Seismic characteristics of a reef carbonate reservoir and implications for hydrocarbon exploration in deepwater of the Qiongdongnan Basin, northern South China Sea

Our analysis of approximately 40,000 km of multichannel 2-D seismic data, reef oil-field seismic data, and data from several boreholes led to the identification of two areas of reef carbonate reservoirs in deepwater areas (water depth >= 500 in) of the Qiongdongnan Basin (QDNB), northern South China Sea. High-resolution sequence stratigraphic analysis revealed that the transgressive and highstand system tracts of the mid-Miocene Meishan Formation in the Beijiao and Ledong-Lingshui Depressions developed reef carbonates. The seismic features of the reef carbonates in these two areas include chaotic bedding, intermittent internal reflections, chaotic or blank reflections, mounded reflections, and apparent amplitude anomalies, similar to the seismic characteristics of the LH11-1 reef reservoir in the Dongsha Uplift and Island Reef of the Salawati Basin, Indonesia, which house large oil fields. The impedance values of reefs in the Beijiao and Ledong-Lingshui Depressions are 8000-9000 g/cc x m/s. Impedance sections reveal that the impedance of the LH11-1 reef reservoir in the northern South China Sea is 800010000 g/cc x m/s, whereas that of pure limestone in BD23-1-1 is > 10000 g/cc x m/s. The mid-Miocene paleogeography of the Beijiao Depression was dominated by offshore and neritic environments, with only part of the southern Beijiao uplift emergent at that time. The input of terrigenous sediments was relatively minor in this area, meaning that terrigenous source areas were insignificant in terms of the Beijiao Depression: reef carbonates were probably widely distributed throughout the depression, as with the Ledong-Lingshui Depression. The combined geological and geophysical data indicate that shelf margin atolls were well developed in the Beijiao Depression, as in the Ledong-Lingshui Depression where small-scale patch or pinnacle reefs developed. These reef carbonates are promising reservoirs, representing important targets for deepwater hydrocarbon exploration. (C) 2008 Elsevier Ltd. All rights reserved.

Seismic stratigraphy of the Qiongdongnan deep sea channel system, northwest South China Sea

Based on more than 4000 km 2D seismic data and seismic stratigraphic analysis, we discussed the extent and formation mechanism of the Qiongdongnan deep sea channel. The Qiongdongnan deep sea channel is a large incised channel which extends from the east boundary of the Yinggehai Basin, through the whole Qiongdongnan and the Xisha trough, and terminates in the western part of the northwest subbasin of South China Sea. It is more than 570 km long and 4-8 km wide. The chaotic (or continuous) middle (or high) amplitude, middle (or high) continuity seismic facies of the channel reflect the different lithological distribution of the channel. The channel formed as a complex result of global sea level drop during early Pliocene, large scale of sediment supply to the Yinggehai Basin, inversion event of the Red River strike-slip fault, and tilted direction of the Qiongdongnan Basin. The large scale of sediment supply from Red River caused the shelf break of the Yinggehai Basin to move torwards the S and SE direction and developed large scale of prograding wedge from the Miocene, and the inversion of the Red River strike-slip fault induced the sediment slump which formed the Qiongdongnan deep sea channel.

Analysis on deep crustal structure along the onshore-offshore seismic profile across the Binhai (littoral) fault zone in northeastern South China Sea

The onshore-offshore deep seismic experiment was carried out for the first time and filled the blankness of the seismic surveys in the transition area between South China and northeastern South China Sea. The seismic data were analyzed and processed. The different seismic phases were identified and their travel time arrivals were modeled by ray-tracing to study the P-wave velocity crustal structure of this area. The crustal structure of this area is the continental crust. The crust thickness is gradually decreasing southward along the on-shore-offshore seismic line. The low-velocity layer (5.5 similar to 5.9 km (.) s(-1)) exists generally in the middle crust (about 10.0 similar to 18.0km)with about 2.5 similar to 4.0 km thickness, which is also thinning seaward. No obvious high-velocity layer appears in the lower crust. The Binhai (littoral) fault zone is a low velocity zone, which is located about 35km southeast to the Nan'ao station and corresponding to the gradient belt of gravity & magnetism anomalies. The depth of the fault zone is close to the Moho discontinuity. The littoral fault zone is a boundary between the normal continental crust of South China and the thinned continental crust of the sea area.

Seismic profiles across the middle Tan-Lu fault zone in Laizhou Bay, Bohai Sea, eastern China

The stratigraphic architecture, structure and Cenozoic tectonic evolution of the Tan-Lu fault zone in Laizhou Bay, eastern China, are analyzed based on interpretations of 31 new 2D seismic lines across Laizhou Bay. Cenozoic strata in the study area are divided into two layers separated by a prominent and widespread unconformity. The upper sedimentary layer is made up of Neogene and Quaternary fluvial and marine sediments, while the lower layer consists of Paleogene lacustrine and fluvial facies. In terms of tectonics, the sediments beneath the unconformity can be divided into four main structural units: the west depression, central uplift, east depression and Ludong uplift. The two branches of the middle Tan-Lu fault zone differ in their geometry and offset: the east branch fault is a steeply dipping S-shaped strike-slip fault that cuts acoustic basement at depths greater than 8 km, whereas the west branch fault is a relatively shallow normal fault. The Tan-Lu fault zone is the key fault in the study area, having controlled its Cenozoic evolution. Based on balanced cross-sections constructed along transverse seismic line 99.8 and longitudinal seismic line 699.0, the Cenozoic evolution of the middle Tan-Lu fault zone is divided into three stages: Paleocene-Eocene transtension, Oligocene-Early Miocene transpression and Middle Miocene to present-day stable subsidence. The reasons for the contrasting tectonic features of the two branch faults and the timing of the change from transtension to transpression are discussed. Crown Copyright (C) 2008 Published by Elsevier Ltd. All rights reserved.

随钻地震技术的研究及应用

Seismic While Drilling (SWD) is a new wellbore seismic technique. It uses the vibrations produced by a drill-bit while drilling as a downhole seismic energy source. The continuous signals generated by the drill bit are recorded by a pilot sensor attached to the top of the drill-string. Seismic wave receivers positioned in the earth near its surface receive the seismic waves both directly and reflection from the geologic formations. The pilot signal is cross-correlated with the receiver signals to compute travel-times of the arrivals (direct arrival and reflected arrival) and attenuate incoherent noise. No downhole intrusmentation is required to obtain the data and the data recording does not interfere with the drilling process. These characteristics offer a method by which borehole seismic data can be acquired, processed, and interpreted while drilling. As a Measure-While-Drill technique. SWD provides real-time seismic data for use at the well site . This can aid the engineer or driller by indicating the position of the drill-bit and providing a look at reflecting horizons yet to be encountered by the drill-bit. Furthermore, the ease with which surface receivers can be deployed makes multi-offset VSP economically feasible. First, this paper is theoretically studying drill-bit wavefield， interaction mode between drill-bit and formation below drill-bit , the new technique of modern signal process was applied to seismic data, the seismic body wave radiation pattern of a working roller-cone drill-bit can be characterized by theoretical modeling. Then , a systematical analysis about the drill-bit wave was done, time-distance equation of seismic wave traveling was established, the process of seismic while drilling was simulated using the computer software adaptive modeling of SWD was done . In order to spread this technique, I have made trial SWD modeling during drilling. the paper sketches out the procedure for trial SWD modeling during drilling , the involved instruments and their functions, and the trial effect. Subsurface condition ahead of the drill-bit can be predicted drillstring velocity was obtained by polit sensor autocorrelation. Reference decovolution, the drillstring multiples in the polit signal are removed by reference deconvolution, the crosscorrelation process enhance the signal-to-noise power ratio, lithologies. Final, SWD provides real-time seismic data for use at the well site well trajectory control exploratory well find out and preserve reservoirs. intervel velocity was computed by the traveltime The results of the interval velocity determination reflects the pore-pressure present in the subsurface units ahead of the drill-bit. the presences of fractures in subsurface formation was detected by shear wave. et al.

碳酸盐岩缝洞型储层地球物理响应特征及预测方法研究

Proven by the petroleum exploration activities, the karsts-fissure reservoir in carbonate rocks is significant to find out the large scale oil & gas field. They are made up of the four reservoir types: karsts-cave, karsts-crack, crack-cave and fracture-pore-cave. Each reservoir space and each reservoir bed has different features of reservoir heterogeneity and small scale of pore-crack-cave. The fracture-cave reservoir in carbonate rocks is characteristic by multi-types and long oiliness well. The reservoir shape is controlled by the irregular pore-crack-cave. The development level of fracture and karst-cave is the key element of hydrocarbon enriching, high productivity and stable production. However, most of Carbonate formation are buried deeply and the signal-ration-noise of seismic reflection are very low. It is reason why the fracture-cave reservoir are difficult to be predicted effectively. In terms of surveyed and studied lots of the former research outcome, The author applied the methods of synthetical reservoir geophysical prediction from two ways including macrosopic and microcomic technics in terms of the reservoir-cap condition, geophysics and geology feature and difficulty of prediction in carbonate rocks. It is guiden by the new ideas of stratigraphy, sedimentology, sedimentography, reservoir geology and karst geology. The geophysics technology is key technics. In aspects of macroscopic studies, starting off the three efficiencies of controlling the reservoir distribution including sedimental facies, karst and fracture, by means of comprehensive utilization of geology, geophysics, boring well and well log, the study of reservoir features and karst inside story are developed in terms of data of individual well and multiple well. Through establishing the carbonate deposition model, karstic model and fracture model, the macro-distribution laws of carbonatite are carried out by the study of coherence analysis, seismic reflection feature analysis and palaeotectonics analysis. In aspects of microcosmic studies, starting off analysis in reservoir geophysical response feature of fracture and karst-cave model according to guidance of the macroscopic geological model in carbonate reservoir, the methods of the carbonate reservoir prediction are developed by comprehensively utilization of seismic multi-attribution intersection analysis, seismic inversion restricted by log, seismic discontinuity analysis, seimic spectrum attenuation gradient, moniliform reflection feature analysis and multiparameter karst reservoir appraisement.Through application of carbonate reservoir synthetical geophysics prediction, the author r successfully develops the beneficial reservoir distribution province in Ordovician of Katake block 1in middle Tarim basin. The fracture-cave reservoir distributions are delineated. The prospect direction and favorable aims are demonstrated. There are a set of carbonate reservoir prediction methods in middle Tarim basin. It is the favorable basic technique in predicting reservoir of the Ordovician carbonate in middle Tarim. Proven by exploration drilling, the favorable region of moniliform reflection fracture and pore-cave and cave-fracture in lower-middle Ordovician are coincidence with the region of hydrocarbon show. It’s indicated that the reservoir prediction methods described in the study of Ordovician carbonate formation are feasible practicably.

柴达木盆地地震勘探关键技术研究

The rugged surface topography determined the seismic data acquisition construction conditions and the seismic wave explosive and receiver quality in Qaidam Basin. This dissertation systematically researched the seismic acquisition, imaging process and the attribute analysis techniques of complicated oil and gas reservoir. The main research achievements and cognitions are as follows: 1. Through the stimulation effects research and analysis from the aspect of lithologic water-containing differences, it’s specific that stable hydrous sand layer can effectively enhance the stimulation effects combined with the corresponding field tests. The seismic data S/N ratio has been improved due to the combination explosive stimulation. Through the fold number and maximum offset analyses of target horizon, the complicated geometry has been optimized and the S/N ratio of seismic data has been improved, which made an important basis for improvement of 3D seismic data. 2. It has been proved that the first arrival refraction static correction method under the model constraint of fine surface survey is suitable to the Qaidam Basin of western areas by the real seismic data processing. Although the refraction horizon of near surface has some changes in a certain extent, it’s steady basically. The refraction horizon can be continuously traced in sections, so it’s qualified for the refraction static correction method on the whole. 3. The research is based on the curved-ray pre-stack time migration techniques of rough topography, and improved the imaging precision of complex areas. This techniques adopted the constant and variable velocity scanning mode and enhanced the velocity analysis precision. The 3D pre-stack time migration techniques reasonably solved the imaging and velocity multiple solutions problems of steep-dip faults and the intersections of horizontal layers. What’s more, fine velocity analysis and mute are very important to enhance the imaging precision of the seismic data in complicated Wunan areas. 4. The 3D seismic data edge-preserving processing methods have been realized due to the image process techniques. Because this method uses the large range filter, it can attenuate the noise maximally. The faults, break points, lithologic pinchout points and lithologic body of small scale such as river will not be influenced by blur because of the edge-preserving characterization of the method which is really an effective assistant technique of low S/N ratio seismic data attribute analysis. 5. The use of spectral decomposition technique can effectively identify the reservoirs. The special geology body which will not be identified (or without obvious characters) in the seismic profile may be found through the details changes of different frequencies in the amplitude profiles.

大港油田复杂断块油气藏评价技术研究与应用

In this paper, in view of characteristics and appraisal difficulty of complicate fault block reservoir, based on the theory of reservoir geology, tectonic geology, sequence stratigraphy, logging geology and sedimentology, according to related multidisciplinary data, such as geological, geophysical, logging and testing, taking 3D seismic acquisition processing, delicate structure interpretation, fine reservoir prediction, reservoir engineering and logging evaluation as research methods, on the basis of comprehensive geologic research, innovative appraisal thinking of complicate fault block reservoir is probed into, reservoir evaluation technology series of different reservoirs were created, and appraisal circuit of complicate fault block reservoir was proposed. According to the research on main controlling factors of hydrocarbon and concentration and concealing laws of Dagang exploration area, family assemble patterns of complicate fault block reservoir were proposed firstly, such as twinning, juxtaposed, overlying, concurrent, connected in series and so on, so theory basement was provided to reservoir evaluation. Taking into account of difficulties of low signal to noise ratio and border identification of lithologic trap of complicate fault block reservoir, in the middle of seismic data acquisition, double focusing layout and optimization technology were adopted, geophysical model was set up, and best observation system was ascertained. During object processing, matching processing technology of removing time difference, phase difference, energy difference and frequency spectrum difference generated by different blazing and receiving elements was studied. Imagery, low amplitude structure and identification of subtle lithologic reservoir was interpreted reasonably. On the basis of characteristics and technical description analysis of structural, structural--lithologic and lithologic reservoirs, innovative appraisal thinking of these reservoirs was formatted, appraisal circuit of complicate fault block reservoir was proposed, and this method could be used in other similar oilfields

地震岩性与物性联合反演技术研究和应用

This paper studies how to more effectively invert seismic data and predict reservoir under complicated sedimentary environment, complex rock physical relationships and fewer drills in offshore areas of China. Based on rock physical and seismic amplitude-preserving process, and according to depositional system and laws of hydrocarbon reservoir, in the light of feature of seismic inversion methods present applied, series methods were studied. A joint inversion technology for complex geological condition had been presented, at the same time the process and method system for reservoir prediction had been established. This method consists four key parts. 1)We presented the new conception called generalized wave impedance, established corresponding inversion process, and provided technical means for joint inversion lithology and petrophysical on complex geological condition. 2)At the aspect of high-resolution nonlinear seismic wave impedance joint inversion, this method used a multistage nonlinear seismic convolution model rather than conventional primary structure Robinson seismic convolution model, and used Caianiello neural network implement inversion. Based on the definition of multistage positive and negative wavelet, it adopted both deterministic and statistical physical mechanism, direct inversion and indirect inversion. It integrated geological knowledge, rock physical theory, well data, and seismic data, and improved the resolution and anti-noise ability of wave impedence inversion. 3)At the aspect of high-resolution nonlinear reservoir physical property joint inversion, this method used nonlinear rock physical model which introduced convolution model into the relationship between wave impedance and porosity/clay. Through multistage decomposition, it handles separately the large- and small-scale components of the impedance-porosity/clay relationships to achieve more accurate rock physical relationships. By means of bidirectional edge detection with wavelets, it uses the Caianiello neural network to finish statistical inversion with combined applications of model-based and deconvolution-based methods. The resulted joint inversion scheme can integrate seismic data, well data, rock physical theory, and geological knowledge for estimation of high-resolution petrophysical parameters. 4)At the aspect of risk assessment of lateral reservoir prediction, this method integrated the seismic lithology identification, petrophysical prediction, multi-scale decomposition of petrophysical parameters, P- and H-spectra, and the match relationship of data got from seismics, well logging and geology. It could describe the complexity of medium preferably. Through applications of the joint inversion of seismic data for lithologic and petrophysical parameters in several selected target areas, the resulted high-resolution lithologic and petrophysical sections(impedance, porosity, clay) show that the joint inversion can significantly improve the spatial description of reservoirs in data sets involving complex deposits. It proved the validity and practicality of this method adequately.

基于弹性参数分析的碳酸盐岩储层预测方法研究

Prediction of Carbonate Reservoir Based on the Elastic Parameter Analysis Zhang Guangzhi (Solid Geophysics) Directed by Professor Liu Hong Abstract With the exploration and development of Puguang Oilfield, oil-gas exploration of carbonate rock in China has shown good prospects. Research on earthquake prediction methods for carbonate reservoir becomes the key of oil and gas exploration. Starting with analysis of geological characteristics of carbonate rock, prestack AVO inversion method, prestack elastic impedance inversion and parameter calculation method and seismic attribute extraction and optimization method were studied based on the analysis of rock physics in this work. First, variation characteristic and law of carbonate rock reservoir parameters were studied based on experimental data of rock physics, log data, analysis assay data, mud logging data and seismic data, so as to lay a foundation for the further reservoir identification and description. Then, the structure, type and propagation law of seismic wave field were analyzed through seismic forward modeling of the reservoir, and contact between information from log and geology data with elastic parameters, such as compressional wave and shear wave velocity and density were established, so as to provide a standard for reservoir identification and hydrocarbon detection using seismic reflection characteristics of the research area. Starting with the general concept of inverse problem, through analysis of Zoeppritz equation, three kinds of pre-stack inversion methods were derived and analyzed in detail, the AVO 3-parameter inversion based on Bayesian theory, the prestack AVO waveform inversion method and the simultaneous inversion method, based on the statistical hypothesis of inversion parameters and observation data and the Gauss distribution assumption of noise. The three methods were validated by model data and real data. Then, the elastic wave impedance inversion method of carbonate reservoir was investigated and the method of elastic parameter extraction from elastic impedance data was put forward. Based on the analysis of conventional methods of seismic attribute extraction and optimization, the time-frequency attributes and the wavelet attributes with time and amplitude feature were presented, and the prestack seismic attribute calculation method which can characterize the reservoir rock and fluid characteristic was presented. And the optimization of seismic attribute using the nonlinear KPCA method was also put forward. A series of seismic prediction technologies for carbonate reservoir were presented based on analysis of rock physics and seismic forward simulation technology. Practical application of these technologies was implemented in A oil field of Southern China and good effect has been achieved. Key words: carbonate rock; reservoir prediction; rock physics, prestack seismic inversion; seismic attribute

叠前地震数据规则化、重建及噪音压制

I address of reconstruction of spatial irregular sampling seismic data to regular grids. Spatial irregular sampling data impairs results of prestack migration, multiple attenuations, spectra estimation. Prestack 5-D volumes are often divided into sub-sections for further processing. Shot gathers are easy to obtain from irregular sampling volumes. My strategy for reconstruction is as follows: I resort irregular sampling gathers into a form of easy to bin and perform bin regularization, then utilize F-K inversion to reconstruct seismic data. In consideration of poor ability of F-K regularization to fill in large gaps, I sort regular sampling gathers to CMP and proposed high-resolution parabolic Radon transform to interpolate data and extrapolate offsets. To strong interfering noise--multiples, I use hybrid-domain high-resolution parabolic Radon transform to attenuate it. F-K regularization demand ultimately for lower computing costs. I proposed several methods to further improve efficiency of F-K inversion: first I introduce 1D and 2D NFFT algorithm for a rapid calculation of DFT operators; then develop fast 1D and 2D CG method to solve least-square equations, and utilize preconditioner to accelerate convergence of CG iterations; what’s more, I use Delaunay triangulation for weight calculation and use bandlimit frequency and varying bandwidth technique for competitive computation. Numerical 2D and 3D examples are offered to verify reasonable results and more efficiency. F-K regularization has poor ability to fill in large gaps, so I rearrange data as CMP gathers and develop hybrid-domain high-resolution parabolic Radon transforms which be used ether to interpolate null traces and extrapolate near and far offsets or suppress a strong interfere noise: multiples. I use it to attenuate multiples to verify performances of our algorithm and proposed routines for industrial application. Numerical examples and field data examples show a nice performance of our method.

地震数据波动方程反演方法研究

Geophysical inversion is a theory that transforms the observation data into corresponding geophysical models. The goal of seismic inversion is not only wave velocity models, but also the fine structures and dynamic process of interior of the earth, expanding to more parameters such as density, aeolotropism, viscosity and so on. As is known to all, Inversion theory is divided to linear and non-linear inversion theories. In rencent 40 years linear inversion theory has formed into a complete and systematic theory and found extensive applications in practice. While there are still many urgent problems to be solved in non-linear inversion theory and practice. Based on wave equation, this dissertation has been mainly involved in the theoretical research of several non-linear inversion methods: waveform inversion, traveltime inversion and the joint inversion about two methods. The objective of gradient waveform inversion is to find a geologic model, thus synthetic seismograms generated by this geologic model are best fitted to observed seismograms. Contrasting with other inverse methods, waveform inversion uses all characteristics of waveform and has high resolution capacity. But waveform inversion is an interface by interface method. An artificial parameter limit should be provided in each inversion iteration. In addition, waveform information will tend to get stuck in local minima if the starting model is too far from the actual model. Based on velocity scanning in traditional seismic data processing, a layer-by-layer waveform inversion method is developed in this dissertation to deal with weaknesses of waveform inversion. Wave equation is used to calculate the traveltime and derivative (perturbation of traveltime with respect to velocity) in wave-equation traveltime inversion (WT). Unlike traditional ray-based travetime inversion, WT has many advantages. No ray tracing or traveltime picking and no high frequency assumption is necessary and good result can be got while starting model is far from real model. But, comparing with waveform inversion, WT has low resolution. Waveform inversion and WT have complementary advantages and similar algorithm, which proves that the joint inversion is a better inversion method. And another key point which this dissertation emphasizes is how to give fullest play to their complementary advantages on the premise of no increase of storage spaces and amount of calculation. Numerical tests are implemented to prove the feasibility of inversion methods mentioned above in this dissertation. Especially for gradient waveform inversion, field data are inversed. This field data are acquired by our group in Wali park and Shunyi district. Real data processing shows there are many problems for waveform inversion to deal with real data. The matching of synthetic seismograms with observed seismograms and noise cancellation are two primary problems. In conclusion, on the foundation of the former experiences, this dissertation has implemented waveform inversions on the basis of acoustic wave equation and elastic wave equation, traveltime inversion on the basis of acoustic wave equation and traditional combined waveform traveltime inversion. Besides the traditional analysis of inversion theory, there are two innovations: layer by layer inversion of seimic reflection data inversion and rapid method for acoustic wave-equation joint inversion.

多波多分量地震资料走时近似方法的研究与应用

Compared with the conventional P wave, multi-component seismic data can markedly provide more information, thus improve the quality of reservoir evaluation like formation evaluation etc. With PS wave, better imaging result can be obtained especially in areas involved with gas chimney and high velocity formation. However, the signal-to-noise of multi-component seismic data is normally lower than that of the conventional P wave seismic data, while the frequency range of converted wave is always close to that of the surface wave which adds to the difficulty of removing surface wave. To realize common reflection point data stacking from extracted common conversion point data is a hard nut to crack. The s wave static correction of common receiver point PS wave data is not easy neither. In a word, the processing of multi-component seismic data is more complicated than P wave data. This paper shows some work that has been done, addressing those problems mentioned above. (1) Based on the AVO feature of converted wave, this paper has realized the velocity spectrum of converted waves by using Sarkar’s generalized semblance method taking into account of AVO factor in velocity analysis. (2)We achieve a method of smoothly offset division normal method.Firstly we scan the stacking velocities in different offset divisions for a t0, secondly obtain some hyperbolas using these stacking velocities, then get the travel time for every trace using these hyperbolas; in the end we interpolate the normal move out between two t0 for every trace. (3) Here realize a method of stepwise offset division normal moveout.It is similar to the method of smoothly offset division normal moveout.The main difference is using quadratic curve, sixth order curve or fraction curve to fit these hyperbolas. (4)In this paper, 4 types of travel time versus distance functions in inhomogeneous media whose velocity or slowness varies with depth and vertical travel time have been discussed and used to approximate reflection travel time. The errors of ray path and travel time based on those functions in four layered models were analyzed, and it has shown that effective results of NMO in synthetic or real data can be obtained. (5) The velocity model of converted PS-wave can be considered as that of P -wave based on the ghost source theory, thus the converted wave travel time can be approximated by calculation from 4 equivalent velocity functions: velocity or slowness vary linearly with depth or vertical travel time. Then combining with P wave velocity analysis, the converted wave data can be corrected directly to the P-wave vertical travel time. The improvements were shown in Normal Move out of converted waves with numerical examples and real data. (6) This paper introduces the methods to compute conversion point location in vertical inhomogeneous media based on linear functions of velocity or slowness versus depth or vertical travel time, and introduce three ways to choose appropriate equivalent velocity methods, which are velocity fitting, travel time approximation and semblance coefficient methods.

非均质储层叠前地震反演技术研究及应用

(1) I research on the relationship between elastic parameters, lithology and liquid. It is a physical base for pre-stack seismic inversion. I research all kinds of approximate expressions of Zoeppritz function. Then the relation of all kinds of approximate expressions can be confirmed. The geological model of water sand and gas sand in different depth was designed. Moreover I research on precision of all kinds of approximate expressions. (2) In process of seismic data which aim at amplitude recovery and apply in pre-stack seismic inversion, I advance to adopt double flow chart for different aim. Pre-stack noise elimination, real amplitude recovery and NMO correction of long offset are the key taches. (3) I made a systemic expatiate for the thinking and applicability about all kinds of expressions of elastic impedance. And mathematical model was applied to compare the precision with all kinds of expressions of elastic impedance. I propose a new pre-stack simultaneous inversion which is based on the Zoppritz function and simulated annealing algorithm. This method can ensure calculation precision of reflection coefficient from different incident angle and get a global optimum solution. Therefore this method improves the precision of pre-stack seismic inversion. (4) The object function of P-S wave pre-stack simultaneous inversion was established. I compared the precision and convergence between simultaneous inversion and P-wave inversion. And the results show that simultaneous inversion is superior to P-wave inversion. Through the study of AVO event of transformed wave, AVO characters of different kinds of gas sand were analyzed. (5) I carried out the study work of pre-stack seismic inversion for carbonate reservoir in middle of Tarim basin and sand shale reservoir in Sulige Area of Erdos Basin. The method and technology in this paper was applied to practical work. And I made a prediction for heterogeneous reservoir. Moreover it acquires a good application effect. Key Word: reflection coefficient， amplitude recovery， pre-stack seismic inversion， Heterogeneous reservoir，prediction.

反假频重建理论及其在储层预测中的应用

To extend the cross-hole seismic 2D data to outside 3D seismic data, reconstructing the low frequency data to high frequency data is necessary. Blind deconvolution method is a key technology. In this paper, an implementation of Blind deconvolution is introduced. And optimized precondition conjugate gradient method is used to improve the stability of the algorithm and reduce the computation. Then high-frequency retrieved Seismic data and the cross-hole seismic data is combined for constraint inversion. Real data processing proved the method is effective. To solve the problem that the seismic data resolution can’t meet the request of reservoir prediction in the river face thin-layers in Chinese eastern oil fields, a high frequency data reconstruction method is proposed. The extrema of the seismic data are used to get the modulation function which operated with the original seismic data to get the high frequency part of the reconstruction data to rebuild the wide band data. This method greatly saves the computation, and easy to adjust the parameters. In the output profile, the original features of the seismic events are kept, the common feint that breaking the events and adding new zeros to produce alias is avoided. And the interbeded details are enhanced compared to the original profiles. The effective band of seismic data is expended and the method is approved by the processing of the field data. Aim to the problem in the exploration and development of Chinese eastern oil field that the high frequency log data and the relative low frequency seismic data can’t be merged, a workflow of log data extrapolation constrained by time-phase model based on local wave decomposition is raised. The seismic instantaneous phase is resolved by local wave decomposition to build time-phase model, the layers beside the well is matched to build the relation of log and seismic data, multiple log info is extrapolated constrained by seismic equiphase map, high precision attributes inverse sections are produced. In the course of resolve the instantaneous phase, a new method of local wave decomposition --Hilbert transform mean mode decomposition(HMMD) is raised to improve the computation speed and noise immunity. The method is applied in the high resolution reservoir prediction in Mao2 survey of Daqing oil field, Multiple attributes profiles of wave impedance, gamma-ray, electrical resistivity, sand membership degree are produced, of which the resolution is high and the horizontal continuous is good. It’s proved to be a effective method for reservoir prediction and estimation.