最优烃类指标方法的研究与应用

Nowadays, with the development of reservoir exploration, the method of exploration is growing. Together with well information and laboratory results, seismic information with high quality can predict reservoir successfully. Hydrocarbon Indicator is a method, which picks the most sensitive rock properties of hydrocarbons, scans the aim area with rock physics tools, and then indicates the area of reservoir. Obviously, the more is the difference between brine and oil/gas, the better this method works. Which parameter can be used as the Optimal Hydrocarbon Indicators is still in discussion. The author introduced several kinds of Hydrocarbon Indicators in this thesis. After analyzing the response of different parameter to reservoir, together with seismic information, the reservoir can be predicted. In this paper, the reservoir of Zhunge’er is studied to prove this kind of method is suitable for real exploration in China. Besides, the author chose Haila’er reservoir to testify whether this method could be used in metamorphic reservoir other than sandstone reservoir. The results highlighted the meaning of Optimal Hydrocarbon Indicators on reservoir identification. The author also mentioned some thoughts for the development of hydrocarbon indicators in the future.

贝尔凹陷苏德尔特构造带古潜山构造与储层研究

With the develop ment of oil and gas exploration, the conventional struc ture exploration era has gradually been substituted by the concealed reser voir exploration technology. Hill poll becomes one of the most important areas in the future exploration. This paper is based on the three-dimensional seismic interpretation of Sudeerte structure. In terms of the overall character istics of Sudeerte structure, we use the coherent cube and the time slice to interpret the fault in plane. For the interpretation，we especially used the well to adjust the interpretation results. The results of seismic attribution analysis, spectrum decomposition and post- stack seismic inversion forecast that hill pools reservoir are dist ributed in several bands along the north-northeast to northeast-east. Xing'anling Group shows that the potential reser voirs are mainly distributed along Bei 14 –Bei 40 and De (99-212) - Bei (16-1) and Budate Group distributed along Bei 14 –Bei 40 in northwestern direction and De (99-212) - Bei (16-1) in north eastern direction. At the same time, by analyzing the structure and the reservoir, and combining with other data, three models are built. The characteristics of reservoirs dist ribution are concl uded，and potent ional favorable exploration dire ctions are predi cted.

数据驱动型多次波衰减方法的研究

In the practical seismic profile multiple reflections tend to impede the task of even the experienced interpreter in deducing information from the reflection data. Surface multiples are usually much stronger, more broadband, and more of a problem than internal multiples because the reflection coefficient at the water surface is much larger than the reflection coefficients found in the subsurface. For this reason most attempts to remove multiples from marine data focus on surface multiples, as will I. A surface-related multiple attenuation method can be formulated as an iterative procedure. In this essay a fully data-driven approach which is called MPI —multiple prediction through inversion (Wang, 2003) is applied to a real marine seismic data example. This is a pretty promising scheme for predicting a relative accurate multiple model by updating the multiple model iteratively, as we usually do in a linearized inverse problem. The prominent characteristic of MPI method lie in that it eliminate the need for an explicit surface operator which means it can model the multiple wavefield without any knowledge of surface and subsurface structures even a source signature. Another key feature of this scheme is that it can predict multiples not only in time but also in phase and in amplitude domain. According to the real data experiments it is shown that this scheme for multiple prediction can be made very efficient if a good initial estimate of the multiple-free data set can be provided in the first iteration. In the other core step which is multiple subtraction we use an expanded multi-channel matching filter to fulfil this aim. Compared to a normal multichannel matching filter where an original seismic trace is matched by a group of multiple-model traces, in EMCM filter a seismic trace is matched by not only a group of the ordinary multiple-model traces but also their adjoints generated mathematically. The adjoints of a multiple-model trace include its first derivative, its Hilbert transform and the derivative of the Hilbert transform. The third chapter of the thesis is the application for the real data using the previous methods we put forward from which we can obviously find the effectivity and prospect of the value in use. For this specific case I have done three group experiments to test the effectiveness of MPI method, compare different subtraction results with fixed filter length but different window length, invest the influence of the initial subtraction result for MPI method. In terms of the real data application, we do fine that the initial demultiple estimate take on a great deal of influence for the MPI method. Then two approaches are introduced to refine the intial demultiple estimate which are first arrival and masking filter respectively. In the last part some conclusions are drawn in terms of the previous results I have got.

复杂介质中地震波场数值模拟研究——基于广义正交多项式迭积微分算子法

In this paper, we propose a new numerical modeling method – Convolutional Forsyte Polynomial Differentiator (CFPD), aimed at simulating seismic wave propagation in complex media with high efficiency and accuracy individually owned by short-scheme finite differentiator and general convolutional polynomial method. By adjusting the operator length and optimizing the operator coefficient, both global and local informations can be easily incorporated into the wavefield which is important to invert the undersurface geological structure. The key issue in this paper is to introduce the convolutional differentiator based on Forsyte generalized orthogonal polynomial in mathematics into the spatial differentiation of the first velocity-stress equation. To match the high accuracy of the spatial differentiator, this method in the time coordinate adopts staggered grid finite difference instead of conventional finite difference to model seismic wave propagation in heterogeneous media. To attenuate the reflection artifacts caused by artificial boundary, Perfectly Matched Layer (PML) absorbing boundary is also being considered in the method to deal with boundary problem due to its advantage of automatically handling large-angle emission. The PML formula for acoustic equation and first-order velocity-stress equation are also derived in this paper. There is little difference to implement the PML boundary condition in all kind of wave equations, but in Biot media, special attenuation factors should be taken. Numerical results demonstrate that the PML boundary condition is better than Cerjan absorbing boundary condition which makes it more suitable to hand the artificial boundary reflection. Based on the theories of anisotropy, Biot two-phase media and viscous-elasticity, this paper constructs the constitutive relationship for viscous-elastic and two-phase media, and further derives the first-order velocity-stress equation for 3D viscous-elastic and two-phase media. Numerical modeling using CFPD method is carried out in the above-mentioned media. The results modeled in the viscous-elastic media and the anisotropic pore elastic media can better explain wave phenomena of the true earth media, and can also prove that CFPD is a useful numerical tool to study the wave propagation in complex media.

尾波干涉法精确测量地震波速变化的方法及实验研究

The coda of seismic waves consists of that part of the signal after the directly arrivials. In a finite medium, or in one that is strongly heterogeneous, the coda is dominated by waves which have repeatedly sampled the medium. Small changes in a medium which may have no detectable influence on the first arrivals are amplified by this repeated sampling and may thus be detectable in the coda. Because of this, coda wave is widely used in detecting micro variations in medium。 In this paper, we give a general view of the theory and application of coda wave, especially coda wave interferometry. We focus on discussing the application of coda wave interferometry on data source of active situ experiment。 First, we apply coda wave interferometry in a short time period situ experiment which last for three days. We also apply the method of coda wave interferometry in a situ experiment which last for one month. Daily circle variations of seismic velocity around the experiment site were obtained, and we also observed that the velocity variations in the experiment site have a significant correlation with the environment factors, including air temperature, barometric pressure, solid earth tide and the level of rainfall. We find that the velocity variation during this period is up to 10-3. The relationship between velocity variation and changes in air temperature, barometric pressure and solid earth tide was analyzed with least square linear fitting .The velocity has no dependence on the air temperature. But velocity has a change of 10-6--10-7 when the barometer or earth tide change per Pa. Generally, we conclude the work and results of previous researchers, and we also display our works and results. We hopes to contribute to the future research of coda wave interferometry.

尺度匹配方法在井震联合重建和多次波衰减中的应用

Scale matching method means adjusting information with different scale to the same level. This thesis focuses on scale unification of information with different frequency bandwidth. Well-seismic cooperate inversion is an important component of reservoir geophysics; multiple prediction & subtraction is a development of multiple attenuation in recent years. The common ground of these two methods is that they both related to different frequency bandwidth unification. Well log、cross-hole seismic、VSP、3D seismic and geological information have different spatial resolution, we can decrease multi-solution of reservoir inversion and enhance the vertical and lateral resolution of the geological object by integrate those information together; Compare the predicted multiple generated by SRME with the real multiple, we find the predicted multiple convolutes at least one wavelet more, which brings frequency bandwidth difference between them. So the subtraction method also relates to multi-scale information unification. This thesis gives a method of well constrained seismic high resolution processing basing on auto gain control modulation. It uses base function method which utilizes original well-seismic match result as initial condition and processed seismic trace as initial model to extrapolate the high frequency information of the well logs to the seismic profiles. In this way we can broaden the bandwidth of the seismic and make the high frequency gain geological meaning. In this thesis we introduce the revised base function method to adaptive subtraction and verify the validity of the method using models. Key words: high frequency reconstruction, scale matching, base function, multiple, SRME prediction & subtraction

离散裂缝诱导各向异性粘弹性地震波模拟与特征分析

The real media always attenuate and distort seismic waves as they propagate in the earth. This behavior can be modeled with a viscoelastic and anisotropic wave equation. The real media can be described as fractured media. In this thesis, we present a high-order staggered grid finite-difference scheme for 2-D viscoelastic wave propagation in a medium containing a large number of small finite length fractures. We use the effective medium approach to compute the anisotropic parameters in each grid cell. By comparing our synthetic seismogram by staggered-grid finite-difference with that by complex-ray parameter ray tracing method, we conclude that the high-order staggered-grid finite-difference technique can effectively used to simulate seismic propagation in viscoelastic-anisotropic media. Synthetic seismograms demonstrate that strong attenuation and significant frequency dispersion due to viscosity are important factors of reducing amplitude and delaying arrival time varying with incidence angle or offset. On the other hand, the amount of scattered energy not only provides an indicator of orientation of fracture sets, but can also provide information about the fracture spacing. Analysis of synthetic seismograms from dry- and fluid-filled fractures indicates that dry-filled fractures show more significant scattering on seismic wavefields than fluid-filled ones, and offset-variations in P-wave amplitude are observable. We also analyze seismic response of an anticlinal trap model that includes a gas-filled fractured reservoir with high attenuation, which attenuates and distorts the so-called bright spot.

中国华北和美国加州地区地震层析成像研究

In recent years seismic tomography has become a powerful tool for studying the three-dimensional crust and mantle structure. In this study, we collected a large number of regional and teleseismic travel-time data and used seismic tomography method to study the relationship between earthquake occurrence and crustal heterogeneity for the 1992 Landers earthquake, heterogeneity and evolution of lithosphere under North China Craton and Southern California, and deep structure and origin of the Changbai intraplate volcano in Northeast China. Our results show a correlation between the seismic rupture zone and crustal heterogeneity. The distribution of the Landers aftershocks is cluster-like and separated or terminated in areas where low-velocity anomalies exist．Most of the large earthquakes with magnitudes >4.0 occurred in or around areas with high P-wave velocity．The possibility is that high-velocity areas are brittle and strong parts which can sustain seismogenic stress，and so can generate earthquakes. Our tomographic images show a very heterogeneous structure in the crust and upper mantle beneath Southern California. Three major anomalies in the upper mantle are revealed clearly beneath the southern Sierra Nevada, Transverse Ranges and Salton Trough. We consider that the high-velocity anomaly beneath the Transverse Ranges was formed through asymmetrical two-side convergence of subcrustal lithosphere and sinking to asthenosphere. Formation of the dense crust root and “drip structure” caused the high-velocity anomaly under the southern Sierra Nevada. The Salton Trough low is the response to the lithospheric extension when the Pacific plate was rifted away from the North American Plate. The tomograpic images beneath the North China Craton show that there exist different lithospheric structures under the different blocks. Complex, prominent low-velocity and high-velocity anomalies are imaged beneath the North China Basin, Trans-North China Orogen (TNCO), and Ordos Block which correspond to rifted, orogenic and cratonic lithospheres, respectively. The thickness of the three-type lithospheres is about 70, 90 and >250 km, respectively. Our results suggest that lithospheric thinning under the eastern part of North China Craton is due to long-term replacement and chemical and thermal erosion of the ancient lithosphere by the hot asthenosphere. The remains of ancient lithosphere exist either in the present upper mantle or have sunk into the mantle transition zone. Our tomographic result of the Changbai volcanic area suggests that the origin of the Changbai volcano is related to the deep dehydration of the subducted Pacific slab and corner flow in the big mantle wedge (BMW) above the stagnant Pacific slab.

波前重建法射线追踪及走时表压缩

As a fast and effective method for approximate calculation of seismic numerical simulation, ray tracing method, which has important theory and practical application value, in terms of seismic theory and seismic simulation, inversion, migration, imaging, simplified from seismic theory according to geometric seismic, means that the main energy of seismic wave field propagates along ray paths in condition of high-frequency asymptotic approximation. Calculation of ray paths and traveltimes is one of key steps in seismic simulation, inversion, migration, and imaging. Integrated triangular grids layout on wavefront with wavefront reconstruction ray tracing method, the thesis puts forward wavefront reconstruction ray tracing method based on triangular grids layout on wavefront, achieves accurate and fast calculation of ray paths and traveltimes. This method has stable and reasonable ray distribution, and overcomes problems caused by shadows in conventional ray tracing methods. The application of triangular grids layout on wavefront, keeps all the triangular grids stable, and makes the division of grids and interpolation of a new ray convenient. This technology reduces grids and memory, and then improves calculation efficiency. It enhances calculation accuracy by accurate and effective description and division on wavefront. Ray tracing traveltime table, which shares the character of 2-D or 3-D scatter data, has great amount of data points in process of seismic simulation, inversion, migration, and imaging. Therefore the traveltime table file will be frequently read, and the calculation efficiency is very low. Due to these reasons, reasonable traveltime table compression will be very necessary. This thesis proposes surface fitting and scattered data compression with B-spline function method, applies to 2-D and 3-D traveltime table compression. In order to compress 2-D (3-D) traveltime table, first we need construct a smallest rectangular (cuboidal) region with regular grids to cover all the traveltime data points, through the coordinate range of them in 2-D surface (3-D space). Then the value of finite regular grids, which are stored in memory, can be calculated using least square method. The traveltime table can be decompressed when necessary, according to liner interpolation method of 2-D (3-D) B-spline function. In the above calculation, the coefficient matrix is stored using sparse method and the liner system equations are solved using LU decomposition based on the multi-frontal method according to the sparse character of the least square method matrix. This method is practiced successfully in several models, and the cubic B-spline function can be the best basal function for surface fitting. It make the construction surface smooth, has stable and effective compression with high approximate accuracy using regular grids. In this way, through constructing reasonable regular grids to insure the calculation efficiency and accuracy of compression and surface fitting, we achieved the aim of traveltime table compression. This greatly improves calculation efficiency in process of seismic simulation, inversion, migration, and imaging.

地震采集观测系统的共聚焦分析

Focal beam analysis is a method for assessment of acquisition geometries that is directly linked to pre-stack migration. About dealing with the complex subsurface structures, the conventional survey design methods which do not take into account the subsurface are no longer valid. Based on the Fourier finite-difference (FFD) large-step wave field extrapolation and Born-Kirchhoff (BK) small-step wavefield interpolation, the thesis presents a rapid resolution analysis of 3D seismic survey design by focal beams in complicated media. Subsequently, The SEG/EAEG salt model is used to illustrate the method. Based on the focal beam resolution definition, each kind of influence factor is discussed. The focal beam analysis usually is carried out in a single frequency, but the actual seismic waves always contain a frequency bandwidth. In this thesis, theoretical relationship between focal beam analysis and frequency is derived. Since the effects of focal beam analysis are linear with frequency simply, the multi-frequency focal beam analysis using interpolation is developed. At the same time, the resolution of different frequency bandwidth is interconvertible in accordance with Signal uncertainty principle. The resolution of all frequency bands can be calculated by using only a few focal beam analysis for a seismic survey. In the last section of this thesis, I propose a new approach to predicting acquisition footprint, based on the assumption of Common-Middle-Point stack without constructing a special velocity model. The approach is a simplistic analytical method in which the acquisition footprint pattern is a weighted, linear summation of limited-offset fold-of-stack plots. Because the value of acquisition can be got by quantificational and rapidly calculating, we can exactly do a comparative analysis among different plans of seismic survey by this method.

复杂地质条件下高精度地震勘探技术研究与应用-以泌阳凹陷为例

In order to carry out high-precision three-dimensional "integration" for the characteristics of the secondary seismic exploration for Biyang Depression, in the implementation process, through a combination of scientific research and production, summed up high-precision seismic acquisition, processing and interpretation technologies suitable for the eastern part of the old liberated areas, achieved the following results: 1. high-precision complex three-dimensional seismic exploration technology series suitable for shallow depression Biyang block group. To highlight the shallow seismic signal, apply goal-based observing system design, trail from the small panel to receive and protect the shallow treatment of a range of technologies; to explain the use of three-dimensional visualization and coherent combination of full-body three-dimensional fine interpretation identification of the 50-100 m below the unconformity surface and its formation of about 10 meters of the distribution of small faults and improve the small block and stratigraphic unconformity traps recognition. 2. high-precision series of three-dimensional seismic exploration technology suitable for deep depression Biyang low signal to noise ratio of information. Binding model using forward and lighting technology, wide-angle observation system covering the design, multiple suppression and raise the energy of deep seismic reflection processing and interpretation of detailed, comprehensive reservoir description, such as research and technology, identified a number of different types of traps. 3. high-precision seismic exploration technology series for the southern Biyang Depression high steep three-dimensional structure. The use of new technology of seismic wave scattering theory and high-precision velocity model based on pre-stack time migration and depth migration imaging of seismic data and other high-precision processing technology, in order to identify the southern steep slope of the local structure prediction and analysis of sandstone bedrock surface patterns provide a wealth of information.

地震资料叠前储层预测综合研究--以辽中凹陷东二下段浊积扇气藏为例

Reservoir prediction techniques from prestack seismic are among the most important ones for exploration of lithologic hydrocarbon reservoir. In this paper, we set the turbidite fan sandstone reservoir in Liao-Zhong depress as our researching target, and aims to solve the apllication difficulties on pre-stack inversion in the area, where the drilling data is scarce and the reservoir is lateral varied. Meanwhile, AVO analysis and pre-stack inversion for gas-bearing detection is systematically researched. The seismic reflection characters of gas-bearing sandstone in turbidite fan with different fluid content are defined, after analyzing results from AVO seismic simulation and porous fluid replacement of real log data, and under the guides of the seismic characters from classical gas-bearing sandstone reservoir and numerical simulation for complicate gas-bearing sandstone. It is confirmed that detecting gas-bearing sandstone in turbidite fan via AVO technologies is feasible. In terms of AVO analysis, two AVO characters, fluid detection factor and product of intercept and gradient, can effectively identify top and bottom boundaries and lateral range of tuibidite gas sand by comparing real drilling data. Cross-plotting of near and far angle stack data could avoid the correlation existing in P-G analysis. After comparing the acoustic impedance inversions with routine stacked data and AVO intercept, impedance derived from AVO intercept attribute could reduce the acoustic impedance estimating error which is caused by AVO. On the aspect of elastic impedance inversion, the AVO information in the pre-stack gathers is properly reserved by creating partial angle stack data. By the far angle elastic impedance alone, the gas sand, with abnormally low range of values, can be identified from the background rocks. The boundary of gas sand can also be clearly determined by cross-plotting of near and far angle elastic impedances. The accuracy of far angle elastic impedance is very sensitive to the parameter K, and by taking the statistical average of Vp/Vs on the targeted section in key wells, the accuracy of low frequency trends is gurranteed; the intensive absorsion within the area of the gas sand, which tends to push the spectral of seismic data to the lower end, will cause errors on the inversion result of elastic impedance. The solution is to confine the inversion on the interested area by improving the wavelet. On the aspect of prestack AVA simultaneous inversion, the constraint of local rock-physical trends between velocities of P-wave、S-wave and density successfully removes the instability of inversion, thus improves the precision of the resulting elastic parameters. Plenty of data on rock properties are derived via AVO analysis and prestack seismic data inversion. Based on them, the fluid anomaly is analysized and lithological interpretation are conducted. The distribution of gas sand can be consistently determined via various of ways, such as cross-plotting of P and G attributes, near and far partial angle stack data, near and far angle elastic impedances, λρ and Vp/Vs, etc. The shear modulo and density are also reliable enough to be used for lithological interpretation. We successfully applied the AVO analysis and pre-stack inversion techniques to gas detecting for turbidite fan sand reservoir in Liao-Zhong depression.

反射地震转换波成像方法研究

Multi-waves and multi-component get more and more attentions from oil industry. On the basis of existent research results, My research focuses on some key steps of OBC 4C datum processing. OBC datum must be preprocessed quite well for getting a good image. We show a flow chart of preprocess including attenuation of noise on multi-component datum、elimination ghost by summing P and Z and rotation of horizontal components. This is a good foundation for the coming steps about OBC processing. How to get exact converted point location and to analyze velocity are key points in processing reflection seismic converted wave data. This paper includes computing converted point location, analyzing velocity and nonhyperbolic moveout about converted waves. Anisotropic affects deeply the location of converted wave and the nonhyperbolic moveout. Supposed VTI, we research anisotropic effect on converted wave location and the moveout. Since Vp/Vs is important, we research the compute method of Vp/Vs from post-stack data and pre-stack data. It is a part of the paper that inversing anisotropic parameter by traveltime. Pre-stack time migration of converted wave is an focus, using common-offset Kirchhoff migration, we research the velocity model updating in anisotropic media. I have achieved the following results: 1) using continued Fractions, we proposed a new converted point approximate equation, when the offset is long enough ,the thomsen’s 2 order equation can’t approximate to the exact location of converted point, our equation is a good approximate for the exact location. 2) our new methods about scanning nonhyperbolic velocity and Vp/Vs can get a high quality energy spectrum. And the new moveout can fit the middle and long offset events. Processing the field data get a good result. 3) a new moveout equation, which have the same form as Alkhalifah’s long offset P wave moveout equation, have the same degree preciseness as thomsen’s moveout equation by testing model data. 4) using c as a function of the ratio offset to depth, we can uniform the Li’s and thomsen’s moveout equation in a same equation, the model test tell us choice the reasonable function C can improve the exact degree of Li’s and thomsen’s equation. 5) using traveltime inversion ,we can get anisotropic parameter, which can help to flat the large offset event and propose a model of anisotropic parameter which will useful for converted wave pre-stack time migration in anisotropic media. 6）using our pre-stack time migration method and flow, we can update the velocity model and anisotropic parameter model then get good image. Key words: OBC, Common converted Point (CCP), Nonhyperbolic moveout equation, Normal moveout correction, Velocity analysis, Anisotropic parameters inversion, Kirchhoff anisotropic pre-stack time migration, migration velocity model updating

印度与欧亚板块强烈相互作用区深部结构特征研究

The central-south Tibet is a part of the products of the continental plate collision between Eurasia and India. To study the deep structure of the study area is significant for understanding the dynamics of the continental-continental collision. A 3-D density model matched well with the observations in the central-south Tibet was proposed in this study. In addition, this study has also used numerical simulation method to prove that Quasi-Love (QL) wave is deduced by anisotropy variation but not by lateral heterogeneity. Meanwhile, anisotropy variation in the upper mantle of the Qiangtang terrane and Lhasa terrane is detected by the QL waves observed in recorded seismograms. Based on the gravity modeling, some results are summarized as follows: 1) Under the constrain of geometrical structure detected by seismic data, a 3-D density model and Moho interface are proposed by gravity inversion of the central-south Tibet. 2) The fact that the lower crustal densities are smaller than 3.2 g/cm3, suggests absence of eclogite or partial eclogitization due to delamination under the central-south Tibet. 3) Seismicity will be strong or weak in the most negative Bouguer gravity anomaly. So there is no a certain relationship between seismicity and Bouguer gravity anomaly. 4) Crustal composition are determined after temperature-pressure calibration of seismic P wave velocity. The composition of lower crust might be one or a mixture of: 1. amphibolite and greenschist facies basalt beneath the Qiangtang terrane; 2. gabbro-norite-troctolite and mafic granulite beneath the Lhasa terrane. Because the composition of the middle crust cannot be well constrained by the above data set, the data set published by Rudnick & Fountain (1995) is used for comparison. It indicated the composition of the middle crust is granulite facies and might be pelitic gneisses.Granulite facies used to be interpreted as residues of partial melting, which coincidences with the previous study on partial melting middle crust. Amphibolite facies are thought to be produced after delamination, when underplating works in the rebound of the lower crust and lithospheric mantle. From the seismology study, I have made several followed conclusions: 1) Through the numerical simulation experiment of surface wave propagating in heterogeneity media, we can find that amplitude and polarization of surface wave only change a little when considering heterogeneity. Furthermore, it is proved that QL waves, generated by surface wave scattering, are caused by lateral variation of anisotropy but not by heterogeneity. 2) QL waves are utilized to determine the variation of uppermost mantle anisotropy of the Tibetan plateau. QL waves are identified from the seismograms of the selected paths recorded by the CAD station. The location of azimuth anisotropy gradient is estimated from the group velocities of Rayleigh wave, Love wave and QL wave. It suggests that south-north lateral variation of azimuthal anisotropy locates in Tanggula mountain, and east-west lateral variation in the north of Gandese mountain with 85°E longitude and near the Jinsha river fault with 85°E longitude.

旋转轴对称介质中的qP波一阶射线追踪

Elastic anisotropy is a very common phenomenon in the Earth’s interior, especial for sedimentary rock as important gas and oil reservoirs. But in the processing and interpretation of seismic data, it is assumption that the media in the Earth’s interior is completely elastic and isotropic, and then the methods based on isotropy are used to deal with anisotropic seismic data, so it makes the seismic resolution lower and the error on images is caused. The research on seismic wave simulation technology can improve our understanding on the rules of seismic wave propagation in anisotropic media, and it can help us to resolve problems caused by anisotropy of media in the processing and interpretation of seismic data. So researching on weakly anisotropic media with rotated axis of symmetry, we study systematically the rules of seismic wave propagation in this kind of media, simulate the process with numerical calculation, and get the better research results. The first-order ray tracing (FORT) formulas of qP wave derived can adapt to every anisotropic media with arbitrary symmetry. The equations are considerably simpler than the exact ray tracing equations. The equations allow qP waves to be treated independently from qS waves, just as in isotropic media. They simplify considerably in media with higher symmetry anisotropy. In isotropic media, they reduce to the exact ray tracing equations. In contrast to other perturbation techniques used to trace rays in weakly anisotropic media, our approach does not require calculation of reference rays in a reference isotropic medium. The FORT-method rays are obtained directly. They are computationally more effective than standard ray tracing equations. Moreover the second-order travel time corrections formula derived can be used to reduce effectively the travel time error, and improve the accuracy of travel time calculation. The tensor transformation equations of weak-anisotropy parameters in media with rotated axis of symmetry derived from the Bond transformation equations resolve effectively the problems of coordinate transformation caused by the difference between global system of coordinate and local system of coordinate. The calculated weak-anisotropy parameters are completely suitable to the first-order ray tracing used in this paper, and their forms are simpler than those from the Bond transformation. In the numerical simulation on ray tracing, we use the travel time table calculation method that the locations of the grids in the ray beam are determined, then the travel times of the grids are obtained by the reversed distance interpolation. We get better calculation efficiency and accuracy by this method. Finally we verify the validity and adaptability of this method used in this paper with numerical simulations for the rotated TI model with anisotropy of about 8% and the rotated ORTHO model with anisotropy of about 20%. The results indicate that this method has better accuracy for both media with different types and different anisotropic strength. Keywords: weak-anisotropy, numerical simulation, ray tracing equation, travel time, inhomogeneity