Anisotropy preserving DTI processing

Statistical analysis of diffusion tensor imaging (DTI) data requires a computational framework that is both numerically tractable (to account for the high dimensional nature of the data) and geometric (to account for the nonlinear nature of diffusion tensors). Building upon earlier studies exploiting a Riemannian framework to address these challenges, the present paper proposes a novel metric and an accompanying computational framework for DTI data processing. The proposed approach grounds the signal processing operations in interpolating curves. Well-chosen interpolating curves are shown to provide a computational framework that is at the same time tractable and information relevant for DTI processing. In addition, and in contrast to earlier methods, it provides an interpolation method which preserves anisotropy, a central information carried by diffusion tensor data. © 2013 Springer Science+Business Media New York.

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.

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

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

裂隙/孔隙型弹性介质中地震波场模拟与叠前逆时偏移

With the development of oil/gas seismic exploration, seismic survey for fracture/porosity type reservoir is becoming more and more important. As for China, since it has over 60% store of low porosity and low permeability oil/gas reservoir, it’s more urgent to validly describe fracture/porosity type oil/gas trap and proposing the related, developed seismic technique. To achieve mapping fracture/porosity region and its development status, it demands profound understanding of seismic wave propagation discipline in complex fractured/pored media. Meanwhile, it has profound scientific significance and applied worth to study forward modeling of fracture/porosity type media and pre-stacked reverse time migration. Especially, pre-stacked reverse-time migration is the lead edge technique in the field of seismology and seismic exploration. In this paper, the author has summarized the meaning, history and the present state of numerical simulation of seismic propagation in fractured/pored media and seismic exploration of fractured/pored reservoirs. Extensive Dilatancy Anisotropy (EDA) model is selected as media object in this work. As to forward modeling, due to local limitation of solving spatial partial derivative when using finite-difference and finite-element method, the author turns to pseudo-spectral method (PSM), which is based on the global characteristic of Fourier transform to simulate three-component elastic wave-field. Artifact boundary effect reduction and simulation algorithm stability are also discussed in the work. The author has completed successfully forward modeling coding of elastic wave-field and numerical simulation of two-dimensional and three-dimensional EDA models with different symmetric axis. Seismic dynamic and kinematical properties of EDA media are analyzed from time slices and seismic records of wave propagation. As to pre-stacked reverse-time migration for elastic wave-field in fractured/pored media, based on the successful experience in forward modeling results with PSM, the author has studied pre-stacked reverse-time depth-domain migration technique using PSM of elastic wave-field in two dimensional EDA media induced by preferred fracture/pore distribution. At the same time, different image conditions will bring up what kind of migration result is detailed in this paper. The author has worded out software for pre-stacked reverse-time depth-domain migration of elastic wave-field in EDA media. After migration processing of a series of seismic shot gathers, influences to migration from different isotropic and anisotropy models are described in the paper. In summary, following creative research achievements are obtained:  Realizing two-dimensional and three-dimensional elastic wave-field modeling for fractured/pored media and related software has been completed.  Proposed pre-stacked reverse-time depth-domain migration technique using PSM of elastic wave-field.  Through analysis of the seismic dynamic and kinematical properties of EDA media, the author made a conclusion that collection of multi-component seismic data can provide important data basis for locating and describing the fracture/pore regions and their magnitudes and the preferred directions.  Pre-stacked reverse-time depth-domain migration technique has the ability to reconstruct complex geological object with steep formations and tilt fracture distribution. Neglecting seismic anisotropy induced by the preferred fracture/pore distribution, will lead to the disastrous imaging results.

裂缝性油气藏储层预测方法研究

In this paper, based on the E & P situation in the oilfield and the theory of geophysical exploration, a series researches are conducted on fracture reservoir prediction technology in general，and it especially focus on some difficult points. The technological series which integrated amplitude preserved data processing、interpretation and its comprehensive application research as a whole were developed and this new method can be applied to the other similar oilfield exploration and development. The contents and results in this paper are listed as follows: 1. An overview was given on the status and development of fracture reservoir estimation technique, compare and analyze those geophysical prediction methods. This will be very helpful to the similar reservoir researches. 2. Analyze and conclude the characters of geologies and well logging response of burial hills fracture reservoir, those conclusions are used to steer the geophysical research and get satisfying results. 3. Forward modeling anisotropy seismic response of fracture reservoir. Quantitatively describe the azimuthal amplitude variation. Amplitude ellipse at each incidence angle is used to identify the fracture orientation. 4. Numerical simulation of structure stress based on finite difference method is carried out. Quantitatively describe and analyze the direction and intensity of fracture. 5. Conventional attributes extraction of amplitude preserved seismic data、attributes with different azimuthal angle and different offset are used to determine the relationship between the results and fracture distribution. 6. With spectrum decomposition method based on wavelet transform， the author disclose the reservoir distribution in space. It is a powerful tool to display its anisotropy. 7. Integrated seismic wave impendence、elastic impendence、spectrum decomposition、attribute extraction、fracture analysis result as a whole to identify and evaluate the fracture reservoir. An optimum workflow is constructed. It is used to practical oil&gas production and good results are obtained. This can indicate the wide foreground of this technique series.

三维VSP技术的研究及应用

Today, because of high petroleum consumption of our country, society steady development and difficulty increase in new resources exploration, deep exploitation of the existing oilfield is needed. More delicate reservoir imaging and description, such as thin layer identification, interlayer exploitation monitoring, subtle structure imaging, reservoir anisotropy recognition, can provide more detail evidence for new development adjustment scheme and enhanced oil recovery. Now, the people have already realized the 3D VSP technique more effective than the general methods in solving these aspects. But VSP technique especially 3D VSP develop slowly due to some reasons. Carrying out the research of VSP technique, it will be very useful to the EOR service. 3D VSP techniques include acquisition、data processing and interpretation. In this paper, the author carried out some researches around acquisition and processing. The key point of acquisition is the survey design, it is critical to the quality of the data and it will influence the reservoir recognition as follows. The author did detailed researches on the layout pattern of shot point and geophone. Some attributes relate to survey design such as reflectivity, incidence angle, observation area, reflection points distribution, fold, minimum well source distance, azimuth angle and so on are studied seriously. In this geometry design of 3D-VSP exploration in deviated wells, the main problems to be solved are: determining the center position of shots distribution, the effect of shots missing on coverage areas and coverage times，locating the shots and receivers of multi-wells. Through simulating and analyzing, the above problems are discussed and some beneficial conclusions are drawn. These will provide valuable references to actual survey design. In data processing, researches emphasize on those relatively key techniques such as wavefield separation, VSP-CDP imaging, the author carried out deep researches around these two aspects. As a result, variant apparent slowness wavefield separation method developed in this article suit the underground variant velocity field and make wavefield separation well, it can overcome reflection bending shortage aroused by conventional imaging method. The attenuateion range of underground seismic wave is very important for amplitude compensation and oil/gas identification.In this paper, seismic wave attenuateion mechanism is studied by 3D-VSP simulateion and Q-inversion technique. By testing with seismic data, the method of VSP data attenuateion and relationship of attenuateion attribute variant with depth is researched. Also the software of survey design and data processing is developed, it fill the gap of VSP area in our country. The technique developed applied successfully in SZXX-A Oilfield、QKYY-B Oilfield、A area and B area. The good results show that this research is valuable, and it is meaningful to the VSP technique development and application of offshore oil industry and other areas in our country.

储层评价目标驱动的地震叠前叠后综合反演方法研究

The exploration and development of natural gas in the north of Ordos basin have been one important part in China’s energy stratagem. Reservoir in upper Palaeozoic group is of lithological trap and its prediction is a crux in a series of works. Based on foregoing seismic reservoir prediction, seismic data are re-processed with some optical methods and pre-stack information is used in corresponding inversions. Through the application of diverse methods, a series of techniques for reservoir prediction come into being. Several results are achieved as flowing: 1. A set of log processing and interpretation methods is developed. Porosity, permeability and gas saturation models are rebuilt. 2. Based on the petro-physics analysis of reservoirs in upper Palaeozoic group, the equations about lithology, property, hydrocarbon and elastic parameters are established. 3. Forward modeling based on elastic wave theory is first applied in the study area and increases the resolution of modeling results. 4. A series of techniques such as pre-stack time migration and others are combined to improve the data quality. 5. Pre-stack seismic inversion is first employed in the north of Ordos Basin and brings the results of EI, P-impedance, S-impedance and other elastic parameters. 6. In post-stack inversion, logs indicating reservoir parameters are rebuilt and boost the resolution of lithology inversion. 7. Amplitude, coherence, frequency-discomposed amplitude, waveform and other sensitive attributes are extracted to describe sands distribution. Seismic modes standing for sands of P1x3, P1x2 are established. 8. Among 9 proposed wells, 8 wells encountered sands and became production wells. The output of DK13 amounts to 510,000 m3 per day. Keywords：the north of Ordos Basin, reservoir prediction, pre-stack inversion, post-stack inversion, seismic attributes．

叠前叠后联合反演技术及其在QHD32-6油田的应用研究

What geophysical inversion studied includes the common mathematics physical property of inversion and the constitution and appraisal method of solution in geophysics domain, i.e. using observed physical phenomenon from the earth surface to infer space changing and physical property structure of medium within the earth. Seismic inversion is a branch of geophysical inversion. The basic purpose of seismic inversion is to utilizing seismic wave propagating law in the medium underground to infer stratum structure and space distribution of physical property according to data acquisition, processing and interpretation, and then offer the vital foundation for exploratory development. Poststack inversion is convenient and swift, its acoustic impedance inversion product can reflect reservoir interior changing rule to a certain degree, but poststack data lack abundant amplitude and travel time information included in prestack data because of multiple superimpose and weaken the sensitiveness which reflecting reservoir property. Compared with poststack seismic inversion, prestack seismic inversion has better fidelity and more adequate information. Prestack seismic inversion, including waveform inversion, not only suitable for thin strata physical property inversion, it can also inverse reservoir oil-bearing ability. Prestack seismic inversion and prestack elastic impedance inversion maintain avo information, sufficiently applying seismic gathering data with different incident angle, partial angle stack, gradient and intercept seismic data cube. Prestack inversion and poststack inversion technology were studied in this dissertation. A joint inversion method which synthesize prestack elastic wave waveform inversion, prestack elastic impedance inversion and poststack inversion was proposed by making fully use of prestack inversion multiple information and relatively fast and steady characteristic of poststack inversion. Using the proposed method to extract rock physics attribute cube with clear physical significance and reflecting reservoir characterization, such as P-wave and S-wave impedance, P-wave and S-wave velocity, velocity ratio, density, Poisson ratio and Lame’s constant. Regarding loose sand reservoir in lower member of Minghuazhen formation, 32-6 south districts in Qinhuangdao,as the research object, be aimed at the different between shallow layer loose sand and deep layer tight sand, first of all, acquire physical property parameters suitable for this kind of heavy oil pool according to experimental study, establishing initial pressure and shear wave relational model; Afterwards, performing prestack elastic wave forward and inversion research, summarizing rules under the guidance of theoretical research and numerical simulation, performing elastic impedance inversion, calculating rock physics attributes; Finally, predicting sand body distribution according to rock physics parameters, and predicting favorable oil area combine well-logging materials and made good results.

海上时移地震综合解释技术研究与应用

The practical application and development of the time-lapse seismic reservoir monitor technology has indicated which has already become one of most important development technologies in seeking the surplus oil distribution and improving the reservoir recovering. The paper, first obtained the rock physics experiment analysis data according to the Bohai Sea loose sandstone in-situ measure technical, and determined the feasibility research of the S oil-field on the time-lapse seismic reservoir monitoring combining with the time-lapse numeric simulation technology, which was used to analyze the time-lapse seismic respond raw of the reservoir parameters change and pointed out the attentive problems during the real time-lapse seismic processing and interpretation. Next, simply introduced the technical link and the effect of the time-lapse mutual constrained fidelity and match processing aiming at the local complex gathering condition, geological condition, development engineering condition. Third, introduced the time-lapse integrated interpretation and the technical system with the innovative key technology that includes the time-lapse difference explanation technology, the time-lapse seismic multi-attributes integrated interpretation technology, and the time-lapse constrained reservoir parameters inversion technology, and so on. Using the time-lapse difference direct explanation technology, directly obtained the surplus oil macroscopic distribution through the difference seismic data; Using the presenting 8 big principles of the sublayer isochronisms comparison, carried on the time-lapse integrated interpretation analysis on the fine sublayer comparison and the thin oil-layer(group) contrast and the oil layer (group); The paper putted up the research, contrast, applications of the multi-sides sensitive attribute analysis and the RBF nerve network on the nearest study algorithm, and predicted the reservoir parameters and the surplus oil distribution with them; Combining with innovative researches and the time-lapse seismic constrained reservoir parameters inversion technology realized the good combination of the seismic and the reservoir engineering. Fourth, under fully analyzing the geology condition, the reservoir condition, the exploit dynamic data, and the seismic data of the S oil-field, and analyzing the time-lapse difference factors with reservoir dynamic exploit data, calibrated the oil-gas saturation change, the pressure change, the water saturation change, and determined the rationality of the time-lapse seismic difference, and finally obtained the surplus oil distribution, the water flood characteristic understanding, reservoir degasification, and pressure drop raw, and so on, which had been used in the well pattern tightening plan proof of the S oil-field development adjustment plan. Finally, the paper summarized the knowledge and understanding of the marine time-lapse seismic integrated interpretation, also had pointed out the further need researched question.

松辽盆地长岭断陷火山岩储层识别与勘探目标评价

Changling fault depression is the biggest fault subsidence in south of Songliao Basin. In its Lower Cretaceous Yingcheng and Shahezi formations developed thick source rocks of deep lake facies and developed poly-phase volcanic rock reservoirs as well. In recent years, significant breakthroughs have been obtained in hydrocarbon exploration of volcanic rock reservoir in the different fault depressions in Songliao basin. Lately, I have been involved in hydrocarbon exploration in the Changling rift depression, especially volcanic rock reservoirs and exploration targets research, participating in the deployment of well Yaoshen 1 which gained over 40 × 104m3 natural gas flow. As quick changes of lithology and facies in Changling area in the south of Songliao basin, and the volcanic rock interludes distribution in continental clastic rock and shale in 3D space, so the identification of volcanic rock types and distribution become a difficult problem. Thus, based on the integrated research of the wild outcrop observation, gravity, magnetic and seismic data, geophysical logging, drilling and coring, laboratory test, this paper carried out the reservoir identification, description and prediction of volcanic rocks in Changling fault depression. In this area, this paper analyzed the volcanic rocks litho-facies, the eruption period, and characteristics of cycles. At the same time, tried to know how to use logging, seismic data to separate volcanic rocks from sandstone and shale, distinguish between volcanic reservoir and non-reservoir, distinguish between intermediate-basic and acidic volcanic rocks, and how to identify traps of volcanic rocks and its gas-bearing properties, etc. Also it is summarized forming conditions and distribution of traps, and possible gas-bearing traps were optimized queuing management. Conclusions as follows: There are two faulted basements in Changling fault depression, granite basement in the southeast and upper paleozoic epimetamorphic basement in the northwest. The main volcanic reservoirs developed in Yingcheng period, which was the intermediate-basic and acidic volcanic eruptions, from the south to north by the intermediate-basic to acid conversion. The volcanic vents are gradually young from south to north. According to information of the re-processing 3D seismic data and gravity-magnetic data, the large volcanic vent or conduit was mainly beaded-distributed along the main fault. The volcanic rocks thickness in Yingcheng formation was changed by the deep faults and basement boundary line. Compared with the clastic rocks, volcanic rocks in Changling area are with high resistance and velocity (4900-5800), abnormal Gamma. All kinds of volcanic rocks are with abnormal strong amplitude reflection on the seismic stacked section except tuff. By analyzing the seismic facies characteristics of volcanic rocks, optimizing seismic attributes constrained by logging, using seismic amplitude and waveforms and other attributes divided volcanic rocks of Yingcheng formation into four seismic zones in map. Currently, most volcanic gas reservoirs are fault-anticline and fault-nose structure. But the volcanic dome lithologic gas reservoirs with large quantity and size are the main gas reservoir types to be found.

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

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.

Pg震相走时的各向异性影响及各向异性结构成像

As the first arrival of seismic phase in deep seismic sounding, Pg is the important data for studying the attributes of the sedimentary layers and the shape of crystalline basement because of its high intensity and reliable detection. Conventionally, the sedimentary cover is expressed as isotropic, linear increasing model in the interpretation of Pg event. Actually, the sedimentary medium should be anisotropic as preferred cracks or fractures and thin layers are common features in the upper crust, so the interpretation of Pg event needs to be taken account of seismic velocity anisotropy. Traveltime calculation is the base of data processing and interpretation. Here, we only study the type of elliptical anisotropy for the poor quality and insufficiency of DSS data. In this thesis, we first investigate the meaning of elliptical anisotropy in the study of crustal structure and attribute, then derive Pg event’s traveltime-offset relationship by assuming a linear increasing velocity model with elliptical anisotropy and present the invert scheme from Pg traveltime-offset dataset to seismic velocity and its anisotropy of shallow crustal structure. We compare the Pg traveltime calculated by our analytic formula with numerical calculating method to test the accuracy. To get the lateral variation of elliptical anisotropy along the profiling, a tomography inversion method with the derived formula is presented, where the profile is divided into rectangles. Anisotropic imaging of crustal structure and attribute is efficient method for crust study. The imaging result can help us interprete the seismic data and discover the attribute of the rock to analyze the interaction between layers. Traveltime calculation is the base of image. Base on the ray tracing equations, the paper present a realization of three dimension of layer model with arbitrary anisotropic type and an example of Pg traveltime calculation in arbitrary anisotropic type is presented. The traveltime calculation method is complex and it only adapts to nonlinear inversion. Perturbation method of travel-time calculation in anisotropy is the linearization approach. It establishes the direct relation between seismic parameters and travetime and it is fit for inversion in anisotropic structural imaging. The thesis presents a P-wave imaging method of layer media for TTI. Southeastern China is an important part of the tectonic framework concerning the continental margin of eastern China and is commonly assumed to comprise the Yangtze block and the Cathaysia block, the two major tectonic units in the region. It’s a typical geological and geophysical zone. In this part, we fit the traveltime of Pg phase by the raytracing numerical method. But the method is not suitable here because the inefficiency of numerical method and the method itself. By the analytic method, we fit the Pg and Sg and get the lateral variation of elliptical anisotropy and then discuss its implication. The northeastern margin of Qinghai-Tibetan plateau is typical because it is the joint area of Eurasian plate and Indian plate and many strong earthquakes have occurred there in recent years．We use the Pg data to get elliptical anisotropic variation and discuss the possible meaning.

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

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.