青藏高原东部电性结构特征及地球动力学意义

As the largest and highest plateau on the Earth, the Tibetan Plateau has been a key location for understanding the processes of mountain building and plateau formation during India-Asia continent-continent collision. As the front-end of the collision, the geological structure of eastern Tibetan Plateau is very complex. It is ideal as a natural laboratory for investigating the formation and evolution of the Tibetan Plateau. Institute of Geophysics, Chinese Academy of Sciences (CAS) carried out MT survey from XiaZayii to Qingshuihe in the east part of the plateau in 1998. After error analysis and distortion analysis, the Non-linear Conjugate Gradient inversion(NLCG), Rapid Relaxation Inversin (RRI) and 2D OCCAM Inversion algorithms were used to invert the data. The three models obtained from 3 algorithms provided similar electrical structure and the NLCG model fit the observed data better than the other two models. According to the analysis of skin depth, the exploration depth of MT in Tibet is much more shallow than in stable continent. For example, the Schmucker depth at period 100s is less than 50km in Tibet, but more than 100km in Canadian Shield. There is a high conductivity layer at the depth of several kilometers beneath middle Qiangtang terrane, and almost 30 kilometers beneath northern Qiangtang terrane. The sensitivity analysis of the data predicates that the depth and resistivity of the crustal high conductivity layer are reliable. The MT results provide a high conductivity layer at 20~40km depth, where the seismic data show a low velocity zone. The experiments show that the rock will dehydrate and partially melt in the relative temperature and pressure. Fluids originated from dehydration and partial melting will seriously change rheological characteristics of rock. Therefore, This layer with low velocity and high conductivity layer in the crust is a weak layer. There is a low velocity path at the depth of 90-110 km beneath southeastern Tibetan Plateau and adjacent areas from seismology results. The analysis on the temperature and rheological property of the lithosphere show that the low velocity path is also weak. GPS measurements and the numerical simulation of the crust-mantle deformation show that the movement rate is different for different terranes. The regional strike derived from decomposition analysis for different frequency band and seismic anisotropy indicate that the crust and upper mantle move separately instead of as a whole. There are material flow in the eastern and southeastern Tibetan Plateau. Therefore, the faults, the crustal and upper mantle weak layers are three different boundaries for relatively movement. Those results support the "two layer wedge plates" geodynamic model on Tibetan formation and evolution.

陕甘宁盆地储层物性预测及油气综合判识方法研究

The thesis show you seven interpretation models of erosion of MAs1+2 in the west of prospect and eroded gully of middle and east, form the interpretation technique how to built up eroded gully of the Ordovician top, according of the practical demand of oil-gas exploration in the ShanGanNing basin, using seismic information, combining well logging and well drilling data, Carefully analyzing geologic deposition background and well logging data, through a great quantity forward and inversion for geologic model and combination geologic model with seismic section. Related to research of reservoir absorption in the ShanGanNing basin, it firstly introduces PRONY transformation multidimensional filter. It can simultaneously express relationship of frequency and absorption decay coefficient, better than FUSAIPU analysis method; PRONY filter have obtain the better effect in the gas field of ZhenChuanBao in the ShanBei area after adopting PRONY filtering method to predict reservoir absorption, by analyzing fixed well and prediction of non-well drilling. In the ShanGanNing basin, general seismic inversion method can produce evident different results or misunderstanding because wave impedance and lithology, physical property, gas property are not sole, especially while have little impedance contrast and even have contract direction; the author carefully analyzes multi-parameter inversion technique, add natural gamma ray and natural potential and other parameter combined making model inversion method according of theory of seismic inversion and applying reservoir velocity and wave impedance information at last, we get the more directly reservoir physical property parameter, judging reservoir physical property is more exact. In accordance with geologic, seismic feature of Shan basin, the thesis conclude Ordovician system top erosion interpretation technology with ChangQing character, and reservoir thickness prediction technique combining inversion technique with wave character analysis, Reservoir physical property that is mainly absorption factor analysis and multi-parameter inversion and oil-gas prediction technology. These technologies obtain the better result in the oil-gas field exploration and have formed comprehensive research method and technology series with ShanGanNing character.

商都盆地二氧化碳气成藏条件研究

The bedding sequences, based on the results from others, have been constructed by geological researches. Furthermore, the reservoir, gas-bearing characteristics and reservoir-blanket association have been increasingly understudied by the geological and seismic studies as well as the log data. The deep dynamics for the formation and development of Shangdu basin resulted from complicated fault system and its continued action have been obtained. The studies on the reservoir condition reveal that the mantle-derived magmatism provided the materials for the CO_2 gas reservoir after Paleogene Period and the huge regional fault not only control the evolution of basin and sedimentary but also pay a role as a passage of the CO_2. The sandstone of river course formed in Paleogene System, with very good reservoir condition, are widely developed in the study area. The blanket with good condition is composed by the basalt in Hannuoba Formation and lake facies shale of Shangdou Formation. Local structures and good encirclement are resulted from the different sedimentary environment and later differential sagging. All statements above demonstrate that there is a very good pool-forming condition for the CO_2. In addition, the high abundance of H_2 recognized during drill exploration are also of significance.More than 30 inorganic CO_2 gas reservoirs have been determined during the exploration for the oil-bearing basins in the eastern China, which are developed along the two sides of Tanlu Fault or within it. In which the CO_2 gas reservoir in Shangdou basin is an inorganic gas reservoir far away from Tanlu Fault. Thus the determination of the CO_2 gas reservoir in Shangdou basin is significant for sciences due to the first exploration for the inorganic CO_2 gas reservoir in our country. The geophysical exploration carried on the CO_2 gas reservoir is benefited for the research of prospecting techniques of CO_2 reservoir.

井间地震技术研究和应用

Cross well seismic technique is a new type of geophysical method, which observes the seismic wave of the geologic body by placing both the source and receiver in the wells. By applying this method, it averted the absorption to high-frequency component of seismic signal caused by low weathering layers, thus, an extremely high-resolution seismic signal can be acquired. And extremely fine image of cross well formations, structure, and reservoir can be achieved as well. An integrated research is conducted to the high-frequency S-wave and P-wave data and some other data to determine the small faults, small structure and resolving the issues concerning the thin bed and reservoir's connectivity, fluid distribution, steam injection and fracture. This method connects the high-resolution surface seismic, logging and reservoir engineering. In this paper, based on the E & P situation in the oilfield and the theory of geophysical exploration, a research is conducted on cross well seismic technology in general and its important issues in cross well seismic technology in particular. A technological series of integrated field acquisition, data processing and interpretation and its integrated application research were developed and this new method can be applied to oilfield development and optimizing oilfield development scheme. The contents and results in this paper are as listed follows: An overview was given on the status quo and development of the cross well seismic method and problems concerning the cross well seismic technology and the difference in cross well seismic technology between China and international levels; And an analysis and comparison are given on foreign-made field data acquisition systems for cross-well seismic and pointed out the pros and cons of the field systems manufactured by these two foreign companies and this is highly valuable to import foreign-made cross well seismic field acquisition system for China. After analyses were conducted to the geometry design and field data for the cross well seismic method, a common wave field time-depth curve equation was derived and three types of pipe waves were discovered for the first time. Then, a research was conducted on the mechanism for its generation. Based on the wave field separation theory for cross well seismic method, we believe that different type of wave fields in different gather domain has different attributes characteristics, multiple methods (for instance, F-K filtering and median filtering) were applied in eliminating and suppressing the cross well disturbances and successfully separated the upgoing and downgoing waves and a satisfactory result has been achieved. In the area of wave field numerical simulation for cross well seismic method, a analysis was conducted on conventional ray tracing method and its shortcomings and proposed a minimum travel time ray tracing method based on Feraiat theory in this paper. This method is not only has high-speed calculation, but also with no rays enter into "dead end" or "blinded spot" after numerous iterations and it is become more adequate for complex velocity model. This is first time that the travel time interpolation has been brought into consideration, a dynamic ray tracing method with shortest possible path has been developed for the first arrivals of any complex mediums, such as transmission, diffraction and refraction, etc and eliminated the limitation for only traveling from one node to another node and increases the calculation accuracy for minimum travel time and ray tracing path and derives solution and corresponding edge conditions to the fourth-order differential sonic wave equation. The final step is to calculate cross well seismic synthetics for given source and receivers from multiple geological bodies. Thus, real cross-well seismic wave field can be recognized through scientific means and provides important foundation to guide the cross well seismic field geometry designing. A velocity tomographic inversion of the least square conjugated gradient method was developed for cross well seismic velocity tomopgraphic inversion and a modification has been made to object function of the old high frequency ray tracing method and put forward a thin bed oriented model for finite frequency velocity tomographic inversion method. As the theory model and results demonstrates that the method is simple and effective and is very important in seismic ray tomographic imaging for the complex geological body. Based on the characteristics of the cross well seismic algorithm, a processing flow for cross well seismic data processing has been built and optimized and applied to the production, a good section of velocity tomopgrphic inversion and cross well reflection imaging has been acquired. The cross well seismic data is acquired from the depth domain and how to interprets the depth domain data and retrieve the attributes is a brand new subject. After research was conducted on synthetics and trace integration from depth domain for the cross well seismic data interpretation, first of all, a research was conducted on logging constraint wave impedance of cross well seismic data and initially set up cross well seismic data interpretation flows. After it applied and interpreted to the cross well seismic data and a good geological results has been achieved in velocity tomographic inversion and reflection depth imaging and a lot of difficult problems for oilfield development has been resolved. This powerful, new method is good for oilfield development scheme optimization and increasing EOR. Based on conventional reservoir geological model building from logging data, a new method is also discussed on constraining the accuracy of reservoir geological model by applying the high resolution cross well seismic data and it has applied to Fan 124 project and a good results has been achieved which it presents a bight future for the cross well seismic technology.

基于信息重建技术的地震信号处理方法研究

The dissertation addressed the problems of signals reconstruction and data restoration in seismic data processing, which takes the representation methods of signal as the main clue, and take the seismic information reconstruction (signals separation and trace interpolation) as the core. On the natural bases signal representation, I present the ICA fundamentals, algorithms and its original applications to nature earth quake signals separation and survey seismic signals separation. On determinative bases signal representation, the paper proposed seismic dada reconstruction least square inversion regularization methods, sparseness constraints, pre-conditioned conjugate gradient methods, and their applications to seismic de-convolution, Radon transformation, et. al. The core contents are about de-alias uneven seismic data reconstruction algorithm and its application to seismic interpolation. Although the dissertation discussed two cases of signal representation, they can be integrated into one frame, because they both deal with the signals or information restoration, the former reconstructing original signals from mixed signals, the later reconstructing whole data from sparse or irregular data. The goal of them is same to provide pre-processing methods and post-processing method for seismic pre-stack depth migration. ICA can separate the original signals from mixed signals by them, or abstract the basic structure from analyzed data. I surveyed the fundamental, algorithms and applications of ICA. Compared with KL transformation, I proposed the independent components transformation concept (ICT). On basis of the ne-entropy measurement of independence, I implemented the FastICA and improved it by covariance matrix. By analyzing the characteristics of the seismic signals, I introduced ICA into seismic signal processing firstly in Geophysical community, and implemented the noise separation from seismic signal. Synthetic and real data examples show the usability of ICA to seismic signal processing and initial effects are achieved. The application of ICA to separation quake conversion wave from multiple in sedimentary area is made, which demonstrates good effects, so more reasonable interpretation of underground un-continuity is got. The results show the perspective of application of ICA to Geophysical signal processing. By virtue of the relationship between ICA and Blind Deconvolution , I surveyed the seismic blind deconvolution, and discussed the perspective of applying ICA to seismic blind deconvolution with two possible solutions. The relationship of PC A, ICA and wavelet transform is claimed. It is proved that reconstruction of wavelet prototype functions is Lie group representation. By the way, over-sampled wavelet transform is proposed to enhance the seismic data resolution, which is validated by numerical examples. The key of pre-stack depth migration is the regularization of pre-stack seismic data. As a main procedure, seismic interpolation and missing data reconstruction are necessary. Firstly, I review the seismic imaging methods in order to argue the critical effect of regularization. By review of the seismic interpolation algorithms, I acclaim that de-alias uneven data reconstruction is still a challenge. The fundamental of seismic reconstruction is discussed firstly. Then sparseness constraint on least square inversion and preconditioned conjugate gradient solver are studied and implemented. Choosing constraint item with Cauchy distribution, I programmed PCG algorithm and implement sparse seismic deconvolution, high resolution Radon Transformation by PCG, which is prepared for seismic data reconstruction. About seismic interpolation, dealias even data interpolation and uneven data reconstruction are very good respectively, however they can not be combined each other. In this paper, a novel Fourier transform based method and a algorithm have been proposed, which could reconstruct both uneven and alias seismic data. I formulated band-limited data reconstruction as minimum norm least squares inversion problem where an adaptive DFT-weighted norm regularization term is used. The inverse problem is solved by pre-conditional conjugate gradient method, which makes the solutions stable and convergent quickly. Based on the assumption that seismic data are consisted of finite linear events, from sampling theorem, alias events can be attenuated via LS weight predicted linearly from low frequency. Three application issues are discussed on even gap trace interpolation, uneven gap filling, high frequency trace reconstruction from low frequency data trace constrained by few high frequency traces. Both synthetic and real data numerical examples show the proposed method is valid, efficient and applicable. The research is valuable to seismic data regularization and cross well seismic. To meet 3D shot profile depth migration request for data, schemes must be taken to make the data even and fitting the velocity dataset. The methods of this paper are used to interpolate and extrapolate the shot gathers instead of simply embedding zero traces. So, the aperture of migration is enlarged and the migration effect is improved. The results show the effectiveness and the practicability.

各向异性与非线性介质中弹性波场传播特征研究

The topic of this study is about the propagation features of elastic waves in the anisotropic and nonlinear media by numerical methods with high accuracy and stability. The main achievements of this paper are as followings: Firstly, basing on the third order elastic energy formula, principle of energy conservation and circumvolved matrix method, we firstly reported the equations of non-linear elastic waves with two dimensions and three components in VTI media. Secondly, several conclusions about some numerical methods have been obtained in this paper. Namely, the minimum suitable sample stepth in space is about 1/8-1/12 of the main wavelength in order to distinctly reduce the numerical dispersion resulted from the numerical mehtod, at the same time, the higher order conventional finite difference (CFD) schemes will give little contribution to avoid the numerical solutions error accumulating with time. To get the similar accuracy with the fourth order center finite difference method, the half truncation length of SFFT should be no less than 7. The FDFCT method can present with the numerical solutions without obvious dispersion when the paprameters of FCT is suitable (we think they should be in the scope from 0.0001 to 0.07). Fortunately, the NADM method not only can reported us with the higher order accuracy solutions (higher than that of the fourth order finite difference method and lower than that of the sixth order finite difference method), but also can distinctly reduce the numerical dispersion. Thirdly, basing on the numerial and theoretical analysis, we reported such nonlinear response accumulating with time as waveform aberration, harmonic generation and resonant peak shift shown by the propagation of one- and two-dimensional non-linear elasticwaves in this paper. And then, we drew the conclusion that these nonlinear responses are controlled by the product between nonlinear strength (SN) and the amplitude of the source. At last, the modified FDFCT numerical method presented by this paper is used to model the two-dimensional non-linear elastic waves propagating in VTI media. Subsequently, the wavelet analysis and polarization are adopted to investigate and understand the numerical results. And then, we found the following principles (attention: the nonlinear strength presented by this paper is weak, the thickness of the -nonlinear media is thin (200m), the initial energy of the source is weak and the anisotropy of the media is weak too): The non-linear response shown by the elastic waves in VTI media is anisotropic too; The instantaneous main frequency sections of seismic records resulted from the media with a non-linear layer have about 1/4 to 1/2 changes of the initial main frequency of source with that resulted from the media without non-linear layer; The responses shown by the elasic waves about the anisotropy and nonlinearity have obvious mutual reformation, namely, the non-linear response will be stronger in some directions because of the anisotropy and the anisotropic strength shown by the elastic waves will be stronger when the media is nonlinear.

共反射面元叠加

The seismic survey is the most effective prospecting geophysical method during exploration and development of oil/gas. The structure and the lithology of the geological body become increasingly complex now. So it must assure that the seismic section own upper resolution if we need accurately describe the targets. High signal/noise ratio is the precondition of high-resolution. As one important seismic data processing method, Stacking is an effective means to suppress the records noise. Broadening area of surface stacked is more important to enhance genuine reflection signals and suppressing unwanted energy in the form of coherent and random ambient noise. Common reflection surface stack is a macro-model independent seismic imaging method. Based on the similarity of CRP trace gathers in one coherent zone, CRS stack effectively improves S/N ratio by using more CMP trace gathers to stack. It is regarded as one important method of seismic data processing. Performing CRS stack depends on three attributes. However, the equation of CRS is invalid under condition of great offset. In this thesis, one method based on velocity model in depth domain is put forward. Ray tracing is used to determine the traveltime of CRP in one common reflection surface by the least squares method to regress the equation of CRS. Then we stack in the coherent seismic data set according to the traveltime, and get the zero offset section. In the end of flowchart of implementing CRS stack, one method using the dip angle to enhance the ratio of S/N is used. Application of the method on synthetic examples and field seismic records, the results of this method show an excellent performance of the algorithm both in accuracy and efficiency.

珠江口盆地流花11-1油田储层非均质性及流体分布规律研究

In order to discover the distribution law of the remaining oil, the paper focuses on the quantitative characterization of the reservoir heterogeneity and the distribution law of the fluid barrier and interbed, based on fine geological study of the reservoir in Liuhuall-1 oil field. The refined quantitative reservoir geological model has been established by means of the study of core analysis, logging evaluation on vertical well and parallel well, and seismic interpretation and prediction. Utilizing a comprehensive technology combining dynamic data with static data, the distribution characteristics, formation condition and controlling factors of remaining oil in Liuhuall-1 oil field have been illustrated. The study plays an important role in the enrichment regions of the remaining oil and gives scientific direction for the next development of the remaining oil. Several achievements have been obtained as follows: l.On the basis of the study of reservoir division and correlation,eight lithohorizons (layer A, B_1, B_2, B_3, C, D, E, and F) from the top to the bottom of the reservoir are discriminated. The reef facies is subdivided into reef-core facies, fore-reef facies and backreef facies. These three subfacies are further subdivided into five microfacies: coral algal limestone, coralgal micrite, coral algal clastic limestone, bioclastic limestone and foraminiferal limestone. In order to illustrate the distribution law of remaining oil in high watercut period, the stratigraphic structure model and sedimentary model are reconstructed. 2.1n order to research intra-layer, inter-layer and plane reservoir heterogeneity, a new method to characterize reservoir heterogeneity by using IRH (Index of Reservoir Heterogeneity) is introduced. The result indicates that reservoir heterogeneity is medium in layer B_1 and B_3, hard in layer A, B_2, C, E, poor in layer D. 3.Based on the study of the distribution law of fluid barrier and interbed, the effect of fluid battier and interbed on fluid seepage is revealed. Fluid barrier and interbed is abundant in layer A, which control the distribution of crude oil in reservoir. Fluid barrier and interbed is abundant relatively in layer B_2,C and E, which control the spill movement of the bottom water. Layer B_1, B_3 and D tend to be waterflooded due to fluid barrier and interbed is poor. 4.Based on the analysis of reservoir heterogeneity, fluid barrier and interbed and the distribution of bottom water, four contributing regions are discovered. The main lies on the north of well LH11-1A. Two minors lie on the east of well LH11-1-3 and between well LH11-1-3 and well LH11-1-5. The last one lies in layer E in which the interbed is discontinuous. 5.The parameters of reservoir and fluid are obtained recurring to core analysis, logging evaluation on vertical well and parallel well, and seismic interpretation and prediction. Theses parameters provide data for the quantitative characterization of the reservoir heterogeneity and the distribution law of the fluid barrier and interbed. 6.1n the paper, an integrated method about the distribution prediction of remaining oil is put forward on basis of refined reservoir geological model and reservoir numerical simulation. The precision in history match and prediction of remaining oil is improved greatly. The integrated study embodies latest trend in this research field. 7.It is shown that the enrichment of the remaining oil with high watercut in Liuhua 11-1 oil field is influenced by reservoir heterogeneity, fluid barrier and interbed, sealing property of fault, driving manner of bottom water and exploitation manner of parallel well. 8.Using microfacies, IRH, reservoir structure, effective thickness, physical property of reservoir, distribution of fluid barrier and interbed, the analysis of oil and water movement and production data, twelve new sidetracked holes are proposed and demonstrated. The result is favorable to instruct oil field development and have gotten a good effect.

保重心脉冲压缩滤波及其在速度自动拾取中的应用

To pick velocity automatically is not only helpful to improve the efficiency of seismic data process, but also to provide quickly the initial velocity for prestack depth migration. In this thesis, we use the Viterbi algorithm to do automatic picking, but the velocity picked usually is immoderate. By thorough study and analysis, we think that the Viterbi algorithm has the function to do quickly and effectually automatic picking, but the data provided for picking maybe not continuous on derivative of its curved surface, viz., the curved face on velocity spectrum is not slick. Therefore, the velocity picked may include irrational velocity information. To solve the problem above, we develop a new method to filter signal by performing nonlinear transformation of coordinate and filter of function. Here, we call it as Gravity Center Preserved Pulse Compressed Filter (GCPPCF). The main idea to perform the GCPPCF as follows: separating a curve, such as a pulse, to several subsection, calculating the gravity center (coordinate displacement), and then assign the value (density) on the subsection to gravity center. When gravity center departure away from center of its subsection, the value assigned to gravity center is smaller than the actual one, but non other than gravity center anastomoses fully with its subsection center, the assigned value equal to the actual one. By doing so, the curve shape under new coordinate breadthwise narrows down compare to its original one. It is a process of nonlinear transformation of coordinate, due to gravity center changing with the shape of subsection. Furthermore, the gravity function is filter one, because it is a cause of filtering that the value assigned from subsection center to gravity center is obtained by calculating its weight mean of subsetion function. In addition, the filter has the properties of the adaptive time delay changed filter, owing to the weight coefficient used for weight mean also changes with the shape of subsection. In this thesis, the Viterbi algorithm inducted, being applied to auto pick the stack velocity, makes the rule to integral the max velocity spectrum ("energy group") forward and to get the optimal solution in recursion backward. It is a convenient tool to pick automatically velocity. The GCPPCF above not only can be used to preserve the position of peak value and compress the velocity spectrum, but also can be used as adaptive time delay changed filter to smooth object curved line or curved face. We apply it to smooth variable of sequence observed to get a favourable source data ta provide for achieving the final exact resolution. If there is no the adaptive time delay-changed filter to perform optimization, we can't get a finer source data and also can't valid velocity information, moreover, if there is no the Viterbi algorithm to do shortcut searching, we can't pick velocity automatically. Accordingly, combination of both of algorithm is to make an effective method to do automatic picking. We apply the method of automatic picking velocity to do velocity analysis of the wavefield extrapolated. The results calculated show that the imaging effect of deep layer with the wavefield extrapolated was improved dominantly. The GCPPCF above has achieved a good effect in application. It not only can be used to optimize and smooth velocity spectrum, but also can be used to perform a correlated process for other type of signal. The method of automatic picking velocity developed in this thesis has obtained favorable result by applying it to calculate single model, complicated model (Marmousi model) and also the practical data. The results show that it not only has feasibility, but also practicability.

基于MRF模型的油气储层属性随机建模方法研究

Stochastic reservoir modeling is a technique used in reservoir describing. Through this technique, multiple data sources with different scales can be integrated into the reservoir model and its uncertainty can be conveyed to researchers and supervisors. Stochastic reservoir modeling, for its digital models, its changeable scales, its honoring known information and data and its conveying uncertainty in models, provides a mathematical framework or platform for researchers to integrate multiple data sources and information with different scales into their prediction models. As a fresher method, stochastic reservoir modeling is on the upswing. Based on related works, this paper, starting with Markov property in reservoir, illustrates how to constitute spatial models for catalogued variables and continuum variables by use of Markov random fields. In order to explore reservoir properties, researchers should study the properties of rocks embedded in reservoirs. Apart from methods used in laboratories, geophysical means and subsequent interpretations may be the main sources for information and data used in petroleum exploration and exploitation. How to build a model for flow simulations based on incomplete information is to predict the spatial distributions of different reservoir variables. Considering data source, digital extent and methods, reservoir modeling can be catalogued into four sorts: reservoir sedimentology based method, reservoir seismic prediction, kriging and stochastic reservoir modeling. The application of Markov chain models in the analogue of sedimentary strata is introduced in the third of the paper. The concept of Markov chain model, N-step transition probability matrix, stationary distribution, the estimation of transition probability matrix, the testing of Markov property, 2 means for organizing sections-method based on equal intervals and based on rock facies, embedded Markov matrix, semi-Markov chain model, hidden Markov chain model, etc, are presented in this part. Based on 1-D Markov chain model, conditional 1-D Markov chain model is discussed in the fourth part. By extending 1-D Markov chain model to 2-D, 3-D situations, conditional 2-D, 3-D Markov chain models are presented. This part also discusses the estimation of vertical transition probability, lateral transition probability and the initialization of the top boundary. Corresponding digital models are used to specify, or testify related discussions. The fifth part, based on the fourth part and the application of MRF in image analysis, discusses MRF based method to simulate the spatial distribution of catalogued reservoir variables. In the part, the probability of a special catalogued variable mass, the definition of energy function for catalogued variable mass as a Markov random field, Strauss model, estimation of components in energy function are presented. Corresponding digital models are used to specify, or testify, related discussions. As for the simulation of the spatial distribution of continuum reservoir variables, the sixth part mainly explores 2 methods. The first is pure GMRF based method. Related contents include GMRF model and its neighborhood, parameters estimation, and MCMC iteration method. A digital example illustrates the corresponding method. The second is two-stage models method. Based on the results of catalogued variables distribution simulation, this method, taking GMRF as the prior distribution for continuum variables, taking the relationship between catalogued variables such as rock facies, continuum variables such as porosity, permeability, fluid saturation, can bring a series of stochastic images for the spatial distribution of continuum variables. Integrating multiple data sources into the reservoir model is one of the merits of stochastic reservoir modeling. After discussing how to model spatial distributions of catalogued reservoir variables, continuum reservoir variables, the paper explores how to combine conceptual depositional models, well logs, cores, seismic attributes production history.

东营凹陷沙三中亚段湖相浊积岩储层特征及油气成藏规律研究

Turbidity sandstone reservoirs have been an important field of hydrocarbon exploration and development in the basins all over the world, as well as in China. Lithologic pools are composed of turbidity sandstones and other sandstones are frequently found in the Jiyang Depression that is a Mesozoic-Cenozoic non-marine oil-bearing basin. The Dongying Sag lies in the sedimentary center of the basin. The subtle traps with turbidity reservoirs are generally difficult to be predicted and described by using current techniques. The studies on turbidity reservoirs plays thus an important theoretical and theoretical practical role in exploration and development in the Jiyang Depression. The attention is, in this thesis, focused on the petrologic properties and oil accumulating behaviors in lake turbidity sedimentary systems in the middle part of the third section of Shahejie Formation in the Dongying Sag, especially in Dongxin area, which lies on the central uplift of the Sag. The paper has disclosed the origin types of turbidity sandstones, distribution pattern and controlling factors of turbidity sandstones, and set up hydrocarbon accumulation patterns of the middle part of the third section of Shahejie Formation in Dongxin, based on nonmarine high resolution sequence stratigraphy, event sedimentology and new theories of hydrocarbon forming. By studying prediction method and technology of turbidity sandstone reservoirs, using precise geological model developing, new techniques of high resolution seismic inversion constrained by logging, the paper has forecast low permeability turbidity sandstone reservoirs and pointed out advantage exploration aims to progressive exploration and development. The paper has obtained mainly many productions and acknowledges as follows: 1.Turbidity sandstone reservoirs of the third section of Shahejie Formationin Dongying Sag are formed in such specifical geological background as rift and extension of basin. The inherited Dongying delta and transgression make up many turbidity distribution areas by overlaying and joining together. The hydrocarbon migrates from depression area to adjacent turbidity sandstone continuously. Accumulation area which is sufficient in oil is formed. 2.The paper has confirmed distinguishable sign of sequence boundary , established stratigraphic framework of Dongying Sag and realized isotime stratigraphic correlation. Es3 of Dongying delta is divided into eleven stages. Among them, the second period of the lower section in Es3, the sixth period of the middle section in Es3, the third period of the upper section in Es3 correspond to eleven sedimentary isotime surface in seismic profile, namely Es3 is classified into eleven Formations. 3.According to such the features of turbidity sandstone as deep in burial, small in area, strong in subtle property, overlaying and joining together and occurring in groups, management through fault and space variations of restriction quantum are realized and the forecast precision of turbidity sandstone by using precise geological model developing, new techniques of high resolution seismic inversion constrained by logging, based on the analysis of all kinds of interwell seismic inversion techniques. 4.According to the features of low permeable turbidity sandstone reservoirs, new method of log interpretation model is put forward. At the same time, distinguish technology of familiar low resistivity oil layer in the turbidity sandstone reservoirs is studied based on petrophysical laboratory work and "four properties" interrelationship between lithological physical Jogging and bearing hydrocarbon properties. Log interpretation model and reservoir index interpretation model of low resistivity oil layer are set up. So the log interpretation precision is improved. 5.The evolution law and its difference of the turbidity sandstone are embodies as follows: the source of sediments come from the south and east of the study area in the middle period of Es3. East source of sediments is pushed from west to east. However, the south source supply of sediments in the early and middle period of Es3 is in full, especially in Es3. subsequently, the supply is decreased gradually. Turbidity fan moves back toward the south and the size of fan is minished accordingly. The characteristic of turbidity sandstone in Dongying Sag is different in different structural positions. Dongxin in the middle-east of the central lift and Niuzhuang Sag He in Dongying delta front and prodelta deep lake subfacies. Although the turbidity sandstone of the two areas root in the Dongying delta sedimentary system, the sand body has different remarkably characteristic. 6.The sedimentary model of the turbiditys in study area have three types as follows: (1) collapse turbidity fan in respect of delta; (2) fault trench turbidity fan; (3) other types of microturbidity sandstone. Middle fan and outer fan, can be found mainly in sublacustrine fan. Middle fan includes braided channel microfacies, central microfacies and braided interchannel microfacies, which is main prospecting oil-bearing subfacies. The middle section of the third section of Shahejie Formation in study area (for example the central lift) can be divided into middle-lower and upper part. The middle-lower part is characteristic of turbidity fan. The upper part is sedimented mainly by delta-collapse fan. 7.The turbidity reservoirs of the middle part of the third section of Shahejie Formation in study area characterize by low maturity both in component and texture, strong in diagenesis and low in permeability. The reservoir can be classified into four types. Type III is the body of reservoir and comprises two types of H a and HI b. M a belongs to middle porosity - low permeability reservoir and distributes in the central lift. Hlb belongs to low porosity - low permeability and distributes in Haojia region. 8.A11 single sand body of lens turbidity reservoir of the middle part of the third section of Shahejie Formation in study area are surrounded by thick dark source rocks. The oil-water system is complex and behaves that every sandstone is single seal unit. The water body is 1/3-1-5 of the sand body. The edge water is not active. The gas exists in the top of reservoir in the form of mixed gas. For far-range turbidity fan with big scale channel, the area and volume of sand body is large and the gap is big in oil packing degree. There are lots of edge water and bottom water, and the latter increases rapidly during the course of development. 9.By exerting the modern hydrocarbon forming theories, the third section of Shahejie Formation in study area belongs to abnormally pressured fluid compartment. The lithological reservoir of the third section of Shahejie Formation is formed in the compartment. The reservoir-formed dynamic system belongs to lower self-source enclosed type. The result and the practice indicate that the form and accumulation of lithological oil reservoirs are controlled by the temperature and pressure of stratum, microfacies, thickness of sand body, fault and reservoir heterogeneity. 10. Based on studies above, the emphases focus on in south and north part of Dongying structure, west Dongxin region and south part Xinzhen structure in the application of production. The practice proves that the turbidity sandstone reservoirs in Ying 11 block and the fault-lithological reservoirs in Xin 133 block have been obtained significant breakthrough. The next target is still sandstone groups of the third section of Shahejie Formation in the bordering areas of Dongxin region for instance Xin 149 area, He 89 area, Ying 8 area etc.

利用分频与反演提高地震资料信噪比与分辨率方法研究

This thesis mainly talks about the wavelet transfrom and the frequency division method. It describes the frequency division processing on prestack or post-stack seismic data and application of inversion noise attenuation, frequency division residual static correction and high resolution data in reservoir inversion. This thesis not only describes the frequency division and inversion in theory, but also proves it by model calculation. All the methods are integrated together. The actual data processing demonstrates the applying results. This thesis analyzes the differences and limitation between t-x prediction filter and f-x prediction filter noise attenuation from wavelet transform theory. It considers that we can do the frequency division attenuation process of noise and signal by wavelet frequency division theory according to the differences of noise and signal in phase, amplitude and frequency. By comparison with the f-x coherence noise, removal method, it approves the effects and practicability of frequency division in coherence and random noise isolation. In order to solve the side effects in non-noise area, we: take the area constraint method and only apply the frequency division processing in the noise area. So it can solve the problem of low frequency loss in non-noise area. The residual moveout differences in seismic data processing have a great effect on stack image and resolutions. Different frequency components have different residual moveout differences. The frequency division residual static correction realizes the frequency division and the calculation of residual correction magnitude. It also solves the problems of different residual correction magnitude in different frequency and protects the high frequency information in data. By actual data processing, we can get good results in phase residual moveout differences elimination of pre-stack data, stack image quality and improvement of data resolution. This thesis analyses the characters of the random noises and its descriptions in time domain and frequency domain. Furthermore it gives the inversion prediction solution methods and realizes the frequency division inversion attenuation of the random noise. By the analysis of results of the actual data processing, we show that the noise removed by inversion has its own advantages. By analyzing parameter's about resolution and technology of high resolution data processing, this thesis describes the relations between frequency domain and resolution, parameters about resolution and methods to increase resolution. It also gives the processing flows of the high resolution data; the effect and influence of reservoir inversion caused by high resolution data. Finally it proves the accuracy and precision of the reservoir inversion results. The research results of this thesis reveal that frequency division noise attenuation, frequency residual correction and inversion noise attenuation are effective methods to increase the SNR and resolution of seismic data.

地球化学新方法在隐状矿床找矿预测中的应用研究

With the progress of prospecting, the need for the discovery of blind ore deposits become more and more urgent. To study and find out the method and technology for the discovery of blind and buried ores is now a priority task. New geochemical methods are key technology to discover blind ores. Information of mobile components related to blind ores were extracted using this new methods. These methods were tested and applied based on element' s mobile components migrating and enriched in geophysical-geochemical process. Several kinds of partial extraction techniques have tested based on element' s occurrence in hypergenic zone. Middle-large scale geochemical methods for exploration in forest and swamp have been tested. A serious of methods were tested and applied effetely about evaluation of regional geochemical anomaly, 1:25000 bedrock or soil geochemical methods sampling based on the net in dendritic water system instead of the normal net. 1. Element related with ores can be mobiled to migrate upwards and be absorpted by surface soil. These abnomal components can be concentrated by natural or artificial methods. These trace metalic ions partially exist in dissovlvable ion forms of active state, and partially have been absorbed by Fe-Mn oxide, soil and organic matter in the soil so that a series of reaction such as complex reaction have take place. Employing various partial extraction techniques, metallic ions related with the phase of the blind ores can be extracted, such as the technique of organic complex extraction, Fe-Mn oxide extraction and the extraction technique of metallic ions of various absorption phases. 2.1:200000 regional geochemical evaluation anomaly methods: Advantageous ore-forming areas were selected firstly. Center, concentration, morphological feature, belt of anomaly were choosed then. Geological and geochemical anomalies were combined. And geological and geochemical background information were restrained. Xilekuduke area in Fuyun sheet , Zhaheba area in Qiakuerte sheet, the west-north part in Ertai sheet and Hongshanzui anomaly in Daqiao sheet were selected as target areas, in Alertai, in the north of Xinjiang. in Xilekuduke area, 1:25000 soil geochemical methods sampling based on the net in dendritic water system was carried out. Cu anomaly and copper mineralization were determined in the center area. Au , Cu anomalies and high polarization anomaly were determined in the south part. Prospecting by primary halo and organic complex extraction were used to prognosis blind ore in widely rang outcrop of bedrock. 1:25000 bedrock or soil geochemical methods sampling based on the net in dendritic water system were used in transported overburden outside of mining area. Shallow seismic method and primary halo found a new blind orebody in mining area. A mineralization site was fou and outside of Puziwan gold mine, in the north of Shanxi province. Developing middle-large scale geochemical exploration method is a key technique based 1:200000 regional geochemical exploration. Some conditions were tested as Sampling density , distribution sites of sample, grain size of sample and occurrence of element for exploration. 1:50000 exploration method was advanced to sample clast sediment supplement clast sediment in valley. 1:25000 bedrock or soil geochemical methods sampling based on the net in dendritic water system was applied to sample residual material in A or C horizon. 1:2000 primary or soil halo methods used to check anomalies and determine mineralization. Daliang gold mineralization in the northern Moerdaoga was found appling these methods. Thermomagnetic method was tested in miniqi copper-polymetallic ore. Process methods such as grain size of sample, heated temperature, magnetic separating technique were tested. A suite of Thermomagnetic geochemical method was formed. This method was applied in Xiangshan Cu~Ni deposit which is cover by clast or Gobi in the eastern Xinjiang. Element's content and contrast of anomaly with Thermomagnetic geochemical method were higher than soil anomaly. Susceptibility after samples were heated could be as a assessment conference for anomaly. In some sectors thermo-magnetic Cu, Ni, Ti anomalious were found outside deposits area. There were strong anomal ies response up ore tested by several kind of partial extraction methods include Thermomagnetic, enzyme leach and other partial extractions in Kalatongke Cu-Ni deposit in hungriness area in the northern of Xinjiang. Element's anomalies of meobile were mainly in Fe-Mn oxide and salt. A Copper mineralization site in Xilekuduke anomaly area had been determined. A blind ore was foung by shallow seismic and geochemical method and a mineralization site was found outside this mining area in Puziwan gold deposit in shanxi province. A Gold mineralization site was found by 1:50000 geochemical exploration in Daliang, Inner Mongolia.

VSP克希霍夫法多波联合成像

Geophones being inside the well, VSP can record upgoing and downgoing P waves, upgoing and downgoing S waves simultaneously.Aiming at overcoming the shortages of the known VSP velocity tomography , attenuation tomography , inverse Q filtering and VSP image method , this article mainly do the following jobs:CD; I do the common-source-point raytracing by soving the raytracing equations with Runge-Kutta method, which can provide traveltime , raypath and amplitude for VSP velocity tomography , attenuation tomography and VSP multiwave migration.(D. The velocity distribution can be inversed from the difference between the computed traveltime and the observed traveltime of the VSP downgoing waves. I put forward two methods: A. VSP building-velocity tomography method that doesn't lie on the layered model from which we can derive the slowness of the grids' crunodes . B. deformable layer tomography method from which we can get the location of the interface if the layer's velocity is known..(3). On the basis of the velocity tomography , using the attenuation information shown by the VSP seismic wave , we can derive the attenuation distribution of the subsurface. I also present an algorithm to solve the inverse Q filtering problem directly and accurately from the Q modeling equation . Numerical results presented have shown that our algorithm gives reliable results . ?. According to the theory that the transformed point is the point where the four kinds of wave come into being , and where the stacked energy will be the largest than at other points . This article presents a VSP multiwave Kirchhoff migration method . Application on synthetic examples and field seismic records have shown that the algorithm gives reliable results . (5). When the location of the interface is determined and the velocity of the P wave and S wave is known , we can obtain the transmittivity and reflection coefficient 5 thereby we can gain the elastic parameters . This method is also put into use derive good result.Above all, application on models and field seismic records show that the method mentioned above is efficient and accurate .

多波分量地震资料处理解释方法

To deal with some key problems in multi-component seismic exploration, some methods are introduced in this thesis based on reading amounts of papers about multi-component seismic theories and methods. First, to find a solution for the detection of the fracture density and orientation in igneous, carbonate and shale reservoirs, a large amount of which exist in domestic oil fields with low exploration and development degree, a new fast and slow shear waves separation method called Ratio Method based on S-wave splitting theory is discussed in this thesis, through which the anisotropy coefficient as well as fracture parameters such as density and azimuthal angle can be acquired. Another main point in this thesis involves the application of seismic velocity ratio (Vp/Vs) to predict the Hthological parameters of subsurface medium. To deal with the unfeasibility of velocity ratio calculation method based on time ratio due to the usually low single-noise ratio of S-wave seismic data acquired on land, a new method based on detailed velocity analysis is introduced. Third, pre-stack Kirchhoff integral migration is a new method developed in recent years, through which both S and P component seismic data as well as amplitude ratio of P/S waves can be acquired. In this thesis, the research on untilizing the P and S wave sections as well as amplitude ratio sections to interpret low-amplitude structures and lithological traps is carried out. The fast and slow shear wave separation method is then be applied respectively to detect the density and azimuthal angle of fractures in an igneous rock gas reservoir and the coal formation in a coal field. Two velocity ratio-calculating methods are applied respectively in the lithological prediction at the gas and coal field after summarizing a large amount of experimental results draw domestically and abroad. P and S wave sections as well as amplitude ratio sections are used to identify low-amplitude structures and lithological traps in the slope area of a oil-bearing sedimentary basin. The calculated data concerning fracture density and azimuthal angle through the introduced method matches well with the regional stress and actual drilling data. The predicted lithological data reflects the actual drilling data. Some of the low-amplitude and lithological traps determined by Kirchhoff migration method are verified by the actual drilling data. These results indicate that these methods are very meaningful when dealing with complex oil and gas reservoir, and can be applied in other areas.