AUV海底热液喷口探测方法研究

自1979年海底热液喷口被首次发现以来，因其巨大的经济和科研价值引起了科学界的巨大关注。海底热液喷口释放的热液与周围海水混合，形成热液羽流，其范围可以达到数千米。热液羽流的存在使在几千米深的海底定位范围只有几米的热液喷口成为可能。湍流的作用使热液羽流与喷口位置存在不确定性，而在搜索区域中包含多个热液源会增加这种不确定性，这是海底热液喷口探测需要克服的难题之一。 本文主要研究了使用AUV探测海底热液喷口的方法。这个问题从更大范畴来说属于机器人化学羽流源定位问题（也称为移动机器人气源/味源定位），其潜在应用包括污染与环境监测，化学工厂安全，搜索与救援，反恐，麻醉品控制，爆炸物清除，以及热液喷口探测等。 首先，从AUV探测的角度研究了海底热液羽流的特性，分析了海底热液羽流的模型并根据模型对羽流进行了动态仿真。 从化学羽流源定位的角度研究了两种海底热液喷口探测策略──梯度搜索策略和构建占据栅格地图(Occupancy Grid Mapping，OGM)的策略。并利用仿真羽流环境验证了上述两种定位策略的可行性。 梯度搜索策略通过基于行为的方法实现，将梯度搜索任务分解为五个行为，并设计了行为间的转换规则，AUV按此规则在不同的行为间转化，跟踪羽流浓度梯度的方向，最终到达浓度极值点。 通过将每个栅格的二值状态重新定义为是否存在一个活跃的热液源，可以将OGM应用于热液源定位。融合传感器数据得到的后验概率地图可以反映每个栅格中存在热液源的可能性。本文采用基于贝叶斯规则的算法融合传感器数据。由于热液源的数量稀少，使用标准贝叶斯方法往往对栅格的占据概率做出过高估计，无法清晰的定位热液源。为此，又研究了一种精确算法和一种基于后验独立假设（Independence of Posteriors，IP）的近似算法，并分析了三种算法的优缺点。 最后，将占据栅格地图应用于分阶段海底热液喷口探测，利用栅格地图帮助实现探测的自主嵌套。

全信息精细储层预测方法研究与实践-以周清庄油田沙三段油藏为例

With the continuously proceeding of petroleum exploratory development in China, exploratory development becomes more and more difficult. For increasing reserve volume and production, lithologic hydrocarbon reservoir has been the most workable, potential and universality exploration targets. In the past, Dagang Oil Field use the complicated fault reservoir theory as the guide, develop and form a suit of matching construction and instrument in prospecting complicated fault reservoir that reach top of exploration industry in China. But the research of lithologic hydrocarbon reservoir is not much, which affects the exploitation progress of lithologic hydrocarbon reservoir. In this thesis, is object, through the depth study of lithologic deposition in Shasan segment of Zhouqingzhuang Oil Field, a suit of holographic fine reservoir bed forecasting techniques is built up and finally gets following main results: 1. Applying geology, seism, drilling, logging and other information to sensitivity preferences, geological model, inversion and integrated stratum evaluation, realizing the method and flow of refined multi-information stratum forecast. 2. Built up a full three dimensional fine structural interpretation method: in view of r problem of accurately demarcating 90% inclined well, propose a inclined well air space demarcating method, make bed demarcating more exactly; in view of problem of faults demarcating and combination in seismic interpretation, propose a computational method of seismic interference based on wavelet translation, make identify the fault in different level more dependable and reasonable; for exactly identifying structural attitude, propose a velocity modeling method under multi-well restriction, make structural attitude closer to the facts. 3. Built up a high accuracy reservoir bed inversion method: in view of problem in exactly identifying reservoir and nonreservoir with conventional wave impedance inversion method in this place, propose a reservoir log response characteristic analysis and sensible log parameter inversion method. ①analysis log response of reservoir and nonreservoir in region of interest, make definite the most sensible log parameter in identifying reservoir and nonreservoir in this region; ②make sensible log parameter inversion based on wave impedance inversion, to improve inversion accuracy, the thickness of recognizable reservoir bed reach 4-5m. 4. Built up a 4-D reservoir forcasting circuit: in view of difficulty that in lithologic hydrocarbon reservoir making reservoir space characteristic clear by using structural map and reservoir forecasting techniques once only, propose a 4-D reservoir forcasting circuit. In other words, based on development conceptual design, forcast reservoir of different time, namely multiple 3D reservoir forcasting in time queue, each time the accuracy degree of reservoir forcasting is improved since apply the new well material, thereby achieve high quality and highly efficient in exploratory development. During exploratory development lithologic depositin in Shasan segment of Zhouqingzhuang Oil Field, there are thirteen wells get 100% success rate, which sufficiently proves that this suit of method is scientific and effective.

地球介质非线性与各向异性的地震波响应及其成像效应

The practice of geophysical prospecting shows us the complex interior earth. The studies of the complexity play an important role and practical guide for the subsurface structure. At present, the complexity of the earth mainly means lateral and vertical homogeneity, anisotropy and non-linear quality. And the anisotropy and non-linear media studies become the frontier in seismology and exploration seismology. This paper summarizes the development of complexities and presents the forward and inverse in the non-linear and anisotropic media. Firstly, the paper introduces the theory of seismic wave propagation in the non-linear and anisotropic media, the theoretical basis for simulation and inversion research. Secondly, high quality numerical simulation method with little dispersion has been developed to investigate the influence of complexity including anisotropy and non-linear multi-component seismograms. Because most real data in seismology have a single component, we developed two aspects work on anisotropic multi-component imaging. One is prestack reflection migration. The result show that distorted images are obtained if data from anisotropic media are migrated using isotropic extrapolation. Moreover, image quality will be improved greatly after considering anisotropy in subsurface layers. The other one is the we take advantage of multi-component data to inversion of the anisotropic parameters jointly seimic reflection travel time and polarization information. Based on these research works, we get the following results: 1.Combing numerical simulation, systematical studies indicate that anisotropy and non-linear seismograms characters are significant to detect cracked belts in the earth and to understand deformation field and mechanism. 2.Based on anisotropic media models, we developed an efficient prestack migration method for subsurface structure and different observation methods seismic data, which improving the imaging quality with VSP, seismograms and real data. 3.Jointly seismic inversion combining seismic anisotropic reflection traveltimes and polarizations data show that the complete wrong inversion and the following explanation will be resulted by ignoring anisotropy.

南海大陆边缘地壳结构与盆地形成机制--以珠江口－琼东南盆地为例

South China Sea is located in the convergence of Eurasian plate, the Pacific Ocean plate and Indian Ocean-Australia plate. The total area is about 3,500,000 km2, the geologic structure is complicated, and the structure line cut off reciprocal is the marginal sea taking form by that the seafloor spreads during the middle Oligocene. South China Sea continental margin have developed more than 10 large oil-gas bearing basins and a number of medium-small sized basins. These basins contain abundant mineral resources such as oil & gas. The marginal deepwater area in the north part of South China Sea has become our country’s strategic energy prospecting frontier. The deepwater area of Zhujiangkou and Qiongdongnan basins is the research target in this thesis. The thesis studied deep structure and the earth dynamics of the north part of South China Sea margin, and these researches provide scientific basis for oil-gas resources strategic investigation and valuation in deepwater sea area of north part slope of South China Sea. In order to develop the research of rebuilding velocities and density architecture of earth shell in region of interest, in marginal deepwater area in the north part of South China, we adopted 14 long-cable seismic reflection profile data of 3556.41 kilometers in total, the gravity measurement data along profiles (3851.44 kilometers in total), the magnetic observation along profiles (3838.4 kilometers in total) and depth measurement along profile, the logging data of 11 wells in project, the interpreted fault parameter and preexisting geologic and geophysical research achievement. This thesis has carried out concretely studying research as follows: 1. Overlay-velocity data sampling and analysis, interval velocity calculation, time-depth conversion, model building of earth shell velocity and layering character of earth shell are studied on 14 deep sections. Velocity structure in region of interest has revealed: Changchang is the sag with thinnest crust in Qiongdongnan basin; the sedimentary thickness lowers gradually from north to south, and the thickness change from west to east is milder. The sags’ sedimentary velocities in Qiongdongnan basin have obvious demarcation. The velocity of the 8000 meters sedimentary rocks is 4700 m/s in Shunde sag and Baiyun sag, and is the lowest; at that depth, the velocity very different in Liwan sag and Baiyun sag, which is about 800m/s. 2. Extracting gravity data and building of initial crust density model along the section; With Bouguer gravity anomaly data as constraint, revising density distributes of initial model, and building the crust density model. 3. With crust velocity and density as constraint, correcting the effect of thermobaric field and constructing constitution structure of rock in region of interest. By this research, we known that rocks in Zhujiangkou upper crustal layer are chiefly granite-gneiss, quartzite, granodiorite and basalt, however, rocks in Qiongdongnan basin upper earth shell are chiefly composed of granite-gneiss, quartzite, granodiorite, diorite and basalt. 4. Synthetically crust velocity and density structure, gaining expanding factor on crust and entire crust along section. The result is indicated: the expanding factor in every sag rises from northwest to southeast, which have reflected thinning characteristic of crust from continent to ocean. Intra-crustal deformation degree in Changchang and Ledong-Lingshui sag is bigger than that in Songnan-Baodao sag. Entire crust extension factor in Changchang and Songnan-Baodao sag is greater than that in Ledong-Lingshui sag, which can make an explanation of frequently event and longer heating process in middle-east of Qiongdongnan basin. 5. Synthesize multidisciplinary information to discuss the earth dynamics significance of discordogenic seismic profile in deepwater area of Zhujiangkou and Qiongdongnan basins.

地震波谱元法数值模拟和偏移研究

With the development of both seismic theory and computer technology, numerical modeling technology of seismic wave has achieved great advancement during the past half century. The current methods under development include finite differentiation method (FDM), finite element method (FEM), pseudospectral method (PSM), integral equation method (IEM) and spectral element method (SEM). They exert their very important roles in every corner of seismology and seismic prospecting. Large quantity of researches towards spectral element method in the end of last century bring this method to a new era, which results in perfect solution of many difficult problems. However, parts of posterior works such as seismic migration and inversion which base on spectral element method have never been studied widely at least up to the present whereas are of importance to seismic imaging and seismic wave propagation. Based on previous work, this paper uses spectral element method to investigate the characteristics and laws of the seismic wave propagation in isotropic and anisotropic media. By thoroughly studying this high-accuracy method, we implement a kind of reverse-time pre- and post-stack migration based on SEM. In order to verify the validity of the SEM method, we have simulated the propagation of seismic wave in several different models. The simulation results show that: (1) spectral element method can be used to model any complex models and the computational results are comparable with the expected results and the analytic results; (2) the optimum accuracy can be achieved when the rank is between 4 and 9. When it is below 4, the dispersion may occur; and when it is above 9, the time step-length will be changed accordingly with the reducing space step-length in order to keep the computation stability. This will exponentially increase the computation time and at the same time the memory even if simulating the same media. This paper also applies explosive reflection surface imaging technology, time constancy principle of wave-filed extrapolation and least travetime raytracing technology of surface source to SEM pre- and post-stack migration of isotropic and anisotropic media. All imaging results derived by the above methods agree well with the real geological models and the position of interface and inflexions can also return to their right location well. This indicates that the method proposed in this paper is a kind of technology with high accuracy and robust stability. It can serve as an alternative method in real seismic data processing. All these work can boost the development of high-accuracy seismic imaging, and therefore have significant inference value.

均衡重力异常计算与综合解释——以黄海地区为例

Petroleum and Natural Gas is an important strategic resources. The reserves of Petroleum and Natural Gas can’t meet the need of our country, which also blocks the development of economy and threatens the safety of national. Therefore, it makes a great sense to bring “the second round of oil & gas exploration” into effect and study the exploration of oil and gas of Pre-Cenozoic residual basins in China. The integrated geophysical exploration is the main way to research the Pre-Cenozoic residual basins. Gravity exploration is one of the most important exploration methods, which has played an important role in oil and gas prospecting, such as compartmentalizing geotectonic elements, delineating the distribution range of sedimentary basins, searching oil and gas structure, abstracting oil and gas information, and so on, from its naissance. The isostatic gravity anomalies is significant for exploration, which can help us research deep crustal structure, the equilibrium state of earth, the geologic structure of shallow crust, the basement shape of sedimentary basins and the genetic evolution of sedimentary basins. In the paper, we stress the implication and physical meanings systemically, and discuss the calculation theory. On the basis of previous work, we test different isostatic compensation models and parameters to find out their influences to the result of isostatic gravity anomalies. In addition, we improve the method of isostatic gravity anomalies calculation and give a system of isostatic gravity anomalies calculation which is proved has satisfying effect. From the research above, we find that the results of Platt model and Airy model are consistent, which have similar form and almost the same value. However, by contrast, the Airy model is proved has better adaptability than Platt model. The two main parameters——crust thickness and density difference of crust and mantle, both have influence to the isostatic gravity anomalies, but the latter have more. Finally, we adopt the regional field extending edge method to make the result more of actual geologic condition. On the methods above, we calculate the isostatic gravity anomalies field in Yellow Sea area from the Bouguer gravity anomalies and the water depth and altitude data. And then the isostatic gravity anomalies character is analyzed and the integrated geological-geophysical interpretation is made on the basis of summarizing the previous research result systemically and analyzing other geophysical data and geological information. From the research, we find that the Yellow Sea area belongs to continental type crust equilibrium regions, where the isostatic gravity anomalies field is placid and has less fluctuation values, which implies that the area is in equilibrium state to different extends.

北部湾盆地迈陈凹陷地震勘探方法研究

Maichen Depression lie between Leizhou Peninsula and Qiongzhou Strait. Oil and gas have been discovered in Weixinan Depression, Wushi Depression and Fushan Depression, which pertain to a same basin — North Sea Basin along with Maichen Depression.Jiangsu Oil started exploration at 2002. The first well began to drill at November, 2004 after gravity survey, electric method prospecting and 2D seismic exploration had been finished. Generating rock and hydrocarbon shows have been verified by the drilling. Low yield oil stream has been tested. And we started 3D seismic exploration at November, 2005. My thesis topic came from the actual needs of our exploration in the Maichen Depression. In the thesis, I give emphasis to analyse the own seismic geologic conditions of Maichen Depression. By real tests, we choosed the means to overcome or weaken the unfavorably impress owing to the own coditions in Maichen Depression. Finally, we obtained the usable seismic data. 1. Owing to the multiphase eruptive rock during the Quaternary Period, the near surface layers are very inhomogeneous. By simultaneous testing at same point with short refraction, uphole surveys of radial source and of surface source, the most appropriate method had been sorted out. Radial source uphole survey has been regarding the best practicable means in the complex area. Accurate surficial geology was very helpful to choosing of acquirement means and parameters. Basically the appropriate method of seismic acquirement has been built at Maichen area. 2. The seismic primary data has many, very strong and complex noise. By noise characteristic analysis in different domain, many means of denoising had been paralleled individual and joint application researched. As a result, the pre-stack multidomain joint denoise flow was the appropriate method. It can improve the seismic signal-to-noise ratio. 3. The problem of seismic static correction at Maichen Depression is very conspicuous. Many static correction methods had been tested individual and joint researched. The seismic data quality has been improved after choosing the appropriate combination of static correction flows. 4. Although the above-mentioned process are resultful, the seismic profile quality is just passable. Some reflector continuity and fault zone imagery are ambiguity. So it was the useful method to reduce the structural ambiguity during seismic interpretation that built-up geologic model in accord with real geologic character by areal structure study upon backbone seismic profiles. In the same way, traps have been assessed and drill targets have been selected.

基于量子蒙特卡罗的地震波反演方法

The primary approaches for people to understand the inner properties of the earth and the distribution of the mineral resources are mainly coming from surface geology survey and geophysical/geochemical data inversion and interpretation. The purpose of seismic inversion is to extract information of the subsurface stratum geometrical structures and the distribution of material properties from seismic wave which is used for resource prospecting, exploitation and the study for inner structure of the earth and its dynamic process. Although the study of seismic parameter inversion has achieved a lot since 1950s, some problems are still persisting when applying in real data due to their nonlinearity and ill-posedness. Most inversion methods we use to invert geophysical parameters are based on iterative inversion which depends largely on the initial model and constraint conditions. It would be difficult to obtain a believable result when taking into consideration different factors such as environmental and equipment noise that exist in seismic wave excitation, propagation and acquisition. The seismic inversion based on real data is a typical nonlinear problem, which means most of their objective functions are multi-minimum. It makes them formidable to be solved using commonly used methods such as general-linearization and quasi-linearization inversion because of local convergence. Global nonlinear search methods which do not rely heavily on the initial model seem more promising, but the amount of computation required for real data process is unacceptable. In order to solve those problems mentioned above, this paper addresses a kind of global nonlinear inversion method which brings Quantum Monte Carlo (QMC) method into geophysical inverse problems. QMC has been used as an effective numerical method to study quantum many-body system which is often governed by Schrödinger equation. This method can be categorized into zero temperature method and finite temperature method. This paper is subdivided into four parts. In the first one, we briefly review the theory of QMC method and find out the connections with geophysical nonlinear inversion, and then give the flow chart of the algorithm. In the second part, we apply four QMC inverse methods in 1D wave equation impedance inversion and generally compare their results with convergence rate and accuracy. The feasibility, stability, and anti-noise capacity of the algorithms are also discussed within this chapter. Numerical results demonstrate that it is possible to solve geophysical nonlinear inversion and other nonlinear optimization problems by means of QMC method. They are also showing that Green’s function Monte Carlo (GFMC) and diffusion Monte Carlo (DMC) are more applicable than Path Integral Monte Carlo (PIMC) and Variational Monte Carlo (VMC) in real data. The third part provides the parallel version of serial QMC algorithms which are applied in a 2D acoustic velocity inversion and real seismic data processing and further discusses these algorithms’ globality and anti-noise capacity. The inverted results show the robustness of these algorithms which make them feasible to be used in 2D inversion and real data processing. The parallel inversion algorithms in this chapter are also applicable in other optimization. Finally, some useful conclusions are obtained in the last section. The analysis and comparison of the results indicate that it is successful to bring QMC into geophysical inversion. QMC is a kind of nonlinear inversion method which guarantees stability, efficiency and anti-noise. The most appealing property is that it does not rely heavily on the initial model and can be suited to nonlinear and multi-minimum geophysical inverse problems. This method can also be used in other filed regarding nonlinear optimization.

奇异性分析方法及其在储层预测中的应用研究

In research field of oil geophysical prospecting, reservoir prediction is refers to forecasting physical properties of petroleum reservoir by using data of seismic and well logging, it is a research which can guide oil field development. Singularities of seismic and logging data are caused by the heterogeneity of reservoir physical property. It's one of important methods that using singularity characteristics of seismic and logging data to study the reservoir physical property in recently. Among them, realization of reservoir quantitative prediction by analyzing singularity of the data and enhancing transition description of data is difficulty in method research. Based on wavelet transform and the fractal theory, the paper studied the singularity judgment criterion for seismic and logging data, not only analyzed quantitative relation between singularity data and reservoir physical property, but also applied it in practical reservoir prediction. The main achievements are: 1. A new method which provides singular points and their strength information estimation at only one single scale is proposed by Herrmann (1999). Based on that, the dissertation proposed modified algorithm which realized singularity polarity detection. 2. The dissertation introduced onset function to generalize the traditional geologic boundaries variations model which used singularity characteristics to represent the abruptness of the lithologic velocity transition. We show that singularity analysis reveals generic singularity information conducted from velocity or acoustic impedance to seismogram based on the convolution seismic-model theory. Theory and applications indicated that singularity information calculated from seismic data was a natural attribute for delineating stratigraphy boundaries due to its excellent ability in detecting detailed geologic features. We demonstrated that singularity analysis was a powerful tool to delineate stratigraphy boundaries and inverse acoustic impedance and velocity. 3. The geologic significances of logging data singularity information were also presented. According to our analysis, the positions of singularities indicate the sequence stratigraphic boundary, and there is subtle relationship between the singularity strength and sedimentary environment, meanwhile the singularity polarity used to recognize stratigraphic base-level cycle. Based on all those above, a new method which provided sedimentary cycle analysis based on the singularity information of logging data in multiple scales was proposed in this dissertation. This method provided a quantitative tool for judging interface of stratum sequence and achieved good results in the actual application.

束鹿凹陷泥灰岩储层综合评价与加砂压裂技术研究与应用

Luo Ning ( Mineralogy, Petrology, Deposit Mineralogy) Directed by Fu Liyun With the increase of the level of exploration and development, North China field, as one of the maturing fields in the east, has gradually turned their prospecting targets to frontiers such as deep zones, lithologic hydrocarbon reservoirs, low permeable layers, special lithostromes, etc, which propose new challenges to mating technique of exploration engineering. In it, the special lithostrome of clay carbonate in Shu-Lu cave in Middle Flank exploration area locates in Es_3 generating rock. The area distribution is large, formation thickness is over 100 meters, the oil accumulation condition is excellent, prognostic reserves is over 80,000,000 tons, but how to effectively stimulate the special low permeable and fractured reservoir has become the bottle neck problem of stimulation and stable yields. In this thesis, through comprehensive evaluation and analysis of lithology, lithomechanics, hydrocarbon reservoir characteristics, the characteristics of fluid flow through porous medium and the stimulation measures in the past, we acquire new cognition of clay carbonate reservoirs, in addition, the research and application of first hydraulic fracturing has gained positive effect and formed commensurable comprehensive reservoir evaluation technique and mating engineering technique of hydraulic fracturing. The main cognitions and achievements are as follows: 1.Study of geological information such as lithololy analysis and nuclear magnetic logging, etc, indicates that clay carbonate formation of Shu-Lu cave is anisotropic, low permeable with high shale content, whose accumulation space gives priority to microcracks. 2.The analysis of lithomechanics of clay carbonate indicates that the hardness is moderate, Young’s modulus is between that of sandstone and limestone, clay carbonate presents plastic property and its breakdown pressure is high because of the deep buried depth. 3.The analysis of the drillstem test curves indicates that the flow and build-up pressure curve of clay carbonate of Shu-Lu cave mainly has three types: formation contamination block-up type, low permeable type, formation energy accumulation slowness type; the reservoir characteristics presents double porosity media, radial compounding, uniform flow vertical fracture, isotropy, moniliform reservoir type. The target well Jingu 3 belongs to moniliform reservoir type. 4.Through recognition and re-evaluation of the treatment effect and technologic limitations of acidizing, acid fracturing and gelled acidizing in the past, based on the sufficient survey and study of hydraulic fracturing home and abroad, combined with comprehensive formation study of target well, we launched the study of the optimization of hydraulic fracturing technique, forming the principal clue and commensurable mating technology aimed at clay carbonate formation, whose targets are preventing leak off, preventing sand bridge, preventing embedment, controlling fracture height, forming long fracture. 5. Recognition of stimulation effect evaluation.

波动方程模拟与偏移的实用化方法研究

Along with the widespread and in-depth applications in petroleum prospecting and development, the seismic modeling and migration technologies are proposed with a higher requirement by oil industrial, and the related practical demand is getting more and more urgent. Based on theories of modeling and migration methods for wave equation, both related with velocity model, I thoroughly research and develop some methods for the goal of highly effective and practical in this dissertation. In the first part, this dissertation probes into the layout designing by wave equations modeling, focusing on the target-oriented layout designing method guided by wave equation modeling in complicated structure areas. It is implemented by using the fourth order staggered grid finite difference (FD) method in velocity-stress 2D acoustic wave equations plus perfectly matched layer (PML) absorbing boundary condition. To design target-oriented layout: (a) match the synthetic record on the surface with events of subsurface structures by analyzing the snapshots of theoretical model; (b) determine the shot-gather distance by tracking the events of target areas and measuring the receiving range when it reaches the surface; (c) restrict the range of valid shot-gather distance by drawing seismic windows in single shot records; (d) choose the best trace distance by comparing the resolution of prospecting targets from the simulated records with different trace distance. Eventually, we obtained the observation system parameters, which achieve the design requirements. In the second part, this dissertation presents the practical method to improve the 3D Fourier Finite Difference (FFD) migration, and carefully analyzes all the factors which influence 3D FFD migration’s efficiency. In which, one of the most important parameters of migration is the extrapolating step. This dissertation presents an efficient 3D FFD migration algorithm, which use FFD propagator to extrapolate wavefields over big layers, and use Born-Kirchhoff interpolator to image wavefields over small layers between the big ones. Finally, I show the effectiveness of this hybrid migration method by comparing migration results from 3D SEG/EAGE model with different methods.