砂泥岩薄互层储层预测方法与应用研究

With the deeply development of exploration and development in petroleum in China, new increasing reserves are found in old oil fields and the verge of the old ones through re-study of geological property. It is more and more important to discovery and develop thin layer or thin inter-bedded layers reservoirs. All of the targets are thin sand-shale inter-bedded reservoirs and the core technology is reservoir predictions between wells in thin sand-shale inter-bedded layers. The continuity of the thin sand-shale inter-bedded layers in space or separating and heterogeneity is the key of reservoir geology research. The seismic reflection, high resolution analysis method and inversion method to thin sand-shale inter-bedded layers are thorough discussed and deeply studied in this paper to try to find the methods and resolutions of reservoir geology research. The below is followed. 1. Based on the pre-research of other people, five models are created: the sand sphenoid body, interlay sandstone and interlay shale of the equal thickness, interlay sandstone of the equal thickness and interlay shale of the unequal thickness, interlay sandstone of the unequal thickness and interlay shale of the unequal thickness, interlay sandstone of the changing thickness in sequence and interlay shale of the changing thickness in sequence. Then the study of the forward modeling are conducted on the thin layer and thin inter-bedded layers geological characters and seismic reflections including amplitude, frequency, phase, wave shape and time-frequency responding in the domains of time and frequency. The affect of petro-physics difference of layers, single thin layer thickness, thickness of inter-bedded, layer number of inter-bedded, incident wavelet domain frequency and types, sample interval to seismic reflection characters, frequency spectrum and time-frequency respond of reflectivity is theoretically discussed. 2. Qualitatively analyzing the sedimentary rhythm of the thin inter-bedded layers in vertical orientation and computing the single layer thickness or the average thickness with the method of generalized S transform. Identifying the reflecting interface or lithology interface using the amplitude value of amplitude spectrum domain frequency. 3. Based on the seismic respond of thin sand-shale inter-bedded layers, bring out the high resolution analysis method of seismic data in thin sand-shale inter-bedded layers using wavelet analysis and the idea of affecting low and high frequency with middle frequency. Then analyzing the effect to the method and testing some wavelets in the method. This method is applied to the theoretical models and the field data. 4. Bring forward one improved very fast simulated annealing method (IVFSA) to resolve the problem nonlinearity and multi-parameters of the inversion in thin inter-bedded layers. And IVFSA is more productive and higher precision than general ways. 5. New target constrained function is used in the inversion based on the property of the inversion in thin inter-bedded layers. 6. Making the full use of geological and logging information, IVFSA and the new function are applied in the non-linear inversion to improve reservoir prediction and evaluation in thin inter-bedded formations combined with the idea of logging and seismic inversion. This method was applied to the field data and got good results.

叠前时间偏移: 基于波动方程的架构、发展和应用

This dissertation starts from the point that the prestack time migration can been considered as an approximation of the prestack depth migration, giving a wave equation based prestack time migration approach. The new approach includes: analytically getting the travel time and amplitude based on the one way wave equation and the stationary-phase theory, using ‘spread’ imaging method and imaging following the prestack depth migration, updating the velocity model with respect to the flats of the events in CRP gathers. Based on this approach, we present a scheme that can image land seismic data without field static correction. We may determine the correct near surface velocities and stack velocities by picking up the residual correction of the events in the CRP gathers. We may get the rational migration section based on the updated velocities and correct the migration section from a floating datum plane to a universal datum plane. We may adaptively determine the migration aperture according to the dips of the imaging structures. This not only speed up the processing, but may suppress the migration noise produce by the extra aperture. We adopt the deconvolution imaging condition of wave equation migration. It may partially compensate the geometric divergence. In this scheme, we use the table-driven technique which may enhance the computational efficiency. If the subsurface is much more complicated, it may be impossible to distinguish the DTS curve. To solve this problem, we proposed a technique to determine the appropriate range of the DTS curve. We synthesize DTS panel in this range using different velocities and depths, and stack the amplitude around the zero time. Determine the correct velocity and location of the considered grid point by comparing the values.

高分辨率非线性岩性物性联合反演方法和应用

In the prediction of complex reservoir with high heterogeneities in lithologic and petrophysical properties, because of inexact data (e.g., information-overlapping, information-incomplete, and noise-contaminated) and ambiguous physical relationship, inversion results suffer from non-uniqueness, instability and uncertainty. Thus, the reservoir prediction technologies based on the linear assumptions are unsuited for these complex areas. Based on the limitations of conventional technologies, the thesis conducts a series of researches on various kernel problems such as inversions from band-limited seismic data, inversion resolution, inversion stability, and ambiguous physical relationship. The thesis combines deterministic, statistical and nonlinear theories of geophysics, and integrates geological information, rock physics, well data and seismic data to predict lithologic and petrophysical parameters. The joint inversion technology is suited for the areas with complex depositional environment and complex rock-physical relationship. Combining nonlinear multistage Robinson seismic convolution model with unconventional Caianiello neural network, the thesis implements the unification of the deterministic and statistical inversion. Through Robinson seismic convolution model and nonlinear self-affine transform, the deterministic inversion is implemented by establishing a deterministic relationship between seismic impedance and seismic responses. So, this can ensure inversion reliability. Furthermore, through multistage seismic wavelet (MSW)/seismic inverse wavelet (MSIW) and Caianiello neural network, the statistical inversion is implemented by establishing a statistical relationship between seismic impedance and seismic responses. Thus, this can ensure the anti-noise ability. In this thesis, direct and indirect inversion modes are alternately used to estimate and revise the impedance value. Direct inversion result is used as the initial value of indirect inversion and finally high-resolution impedance profile is achieved by indirect inversion. This largely enhances inversion precision. In the thesis, a nonlinear rock physics convolution model is adopted to establish a relationship between impedance and porosity/clay-content. Through multistage decomposition and bidirectional edge wavelet detection, it can depict more complex rock physical relationship. Moreover, it uses the Caianiello neural network to implement the combination of deterministic inversion, statistical inversion and nonlinear theory. Last, by combined applications of direct inversion based on vertical edge detection wavelet and indirect inversion based on lateral edge detection wavelet, it implements the integrative application of geological information, well data and seismic impedance for estimation of high-resolution petrophysical parameters (porosity/clay-content). These inversion results can be used to reservoir prediction and characterization. Multi-well constrains and separate-frequency inversion modes are adopted in the thesis. The analyses of these sections of lithologic and petrophysical properties show that the low-frequency sections reflect the macro structure of the strata, while the middle/high-frequency sections reflect the detailed structure of the strata. Therefore, the high-resolution sections can be used to recognize the boundary of sand body and to predict the hydrocarbon zones.

复杂地区地震成像技术研究

Aiming at solving the seismic imaging difficulty in complex area, the static correction methods and the migration imaging techniques taking the anisotropy into account are studied in this dissertation. To solve the static correction problems, a new tomography inversion approach is presented which takes use of the apparent slowness and apparent velocity and inverts both head and diving waves over the complete offset ranges. This approach is also taken practice to the practical seismic data processing of south areas of China and gets ideal effects. There are obvious differences between the actual statics and the statics based on the surface consistency hypothesis. In this dissertation, the exact differences formula is derived. The wave-eqation datuming method based on a single shot gather and the hybrid zero-offset wave-equation datuming algorithm based on f-x domain and f-k domain are presented at the same time. Further more, some forward modelings are made and tested. These methods are also put into practical seismic data processing and good results are made. In this dissertation, the true amplitude Kirchhoff pre-stack time migration fomula in VTI media is presented. The high-dense bispectral scanning technique based on the anelliptical time-shifted hyperbola and the geostatistical filtering are adopted to extract the anellipticity parameter. Simultaneously, combined with the practical seismic data imaging, the anisotropic pre-stack time migration flow is proposed and good processing results are made.

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

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

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

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

利用台阵面波资料研究华北克拉通北缘燕山褶皱带及其邻区的地壳与上地幔结构

The continent of eastern China, especially the North China Craton (NCC), has endured intensive tectonic renovation during Mesozoic and Cenozoic, with the presence of widespread magmatism, high heat flow and development of large sedimentary basins and mountain ranges. The cratonic lithosphere of the region has been destroyed remarkably, which is characterized by not only a significant reduction in thickness but also complex modifications in physical and chemical properties of the lithosphere. As for the tectonic regime controlling the evolution of the NCC, various models have been put forward, including the impingement of mantle plumes (“mushroom cloud” model), the collision of south China block and north China block, the subduction of the Pacific plate, etc. Lithosphere delamination and thermal erosion were proposed as the two end-member mechanisms of the lithospheric thinning. However, given the paucity of the data, deep structural evidence is currently still scarce for distinguishing and testifying these models. To better understand the deep structure of the NCC, from 2000 to the present, temporary seismic array observations have been conducted in the NCC by the Seismological Laboratory of the Institute of the Geology and Geophysics, Chinese Academy of Sciences under the North China Interior Structure Project (NCISP). Many arrays extend from the North China Craton and the off-craton regions, and traverse a lot of main tectonic boundaries. A total of more than 300 broadband seismic stations have been deployed along several profiles that traversed the major tectonic units within the craton’s interior, at the boundary areas and in the neighboring off-craton regions. These stations recorded abundant high-quality data, which provides an unprecedented opportunity for us to unravel the deep structural features of the NCC using seismological methods. Among all the seismological methods, the surface wave method appears to be an efficient and widely adopted technique in studying the crustal and upper mantle structures. In particular, it can provide the absolute values of S-wave velocity that are difficult to obtain with other methods. Benefiting from the deployment of dense seismic arrays, progresses have been made in improving the spatial resolution of surface wave imaging, which makes it possible to resolve the fine-scale velocity structures of the crust and upper mantle based on surface wave analysis. Meanwhile, the differences in the S-wave velocities derived from Rayleigh and Love wave data can provide information on the radial anisotropy beneath the seismic arrays. In this thesis, using the NCISP-III broadband data and based on phase velocity dispersion analysis and inversion of fundamental mode Rayleigh and Love waves, I investigated the lateral variations in the S-wave velocity structure of the crust and uppermost mantle beneath the Yanshan Belt and adjacent regions at the northeastern boundary of the NCC. Based on the constructed structural images, I discussed possible deep processes of the craton destruction in the study region.

利用接收函数方位变化研究地壳各向异性

Because of its sensitivity to the velocity discontinuity of the earth, receiver function technique has become a routine procedure used to probe interior structure of the earth. Receiver functions contain anisotropic information of the earth’s interior, however, traditional receiver function techniques such as migration imaging and waveform inversion method, which are based on isotropic media assumption, can not effectively extract the anisotropy information contained in the azimuth variation pattern. Only by using the anisotropic media, e.g. a model with symmetric axis of arbitrary orientation, computing the response, can we obtain the detailed anisotropy information hidden in the radial and transversal receiver function. Focusing on the receiver function variation pattern changing wtih different back azimuths, we introduced different kinds of symmetric systems of seismic anisotropy used often, and summarized some possible causes of anisotropy formation. We show details about how to calculate the response of a stratified anisotropy model with symmetric axis of arbitrary orientation. We also simulated receiver functions among different models and analyzed how the changing of anisotropic parameters influence the azimuth variation pattern of receiver functions. The anisotropy study by receiver function analysis was applied to Taihang Mountain Range (TMR) in North China in this thesis. The maximum entropy spectrum deconvolution technique was used to extract radial and transversal receiver functions from the waveforms of 20 portable seismic stations deployed in TMR. Considering the signal-to-noise ratio and the azimuth coverage, we got the variation pattern of receiver functions for 11 stations. After carefully analyzing the pattern of the receiver functions that we got, we obtained the reliable evidence on the existence of anisotropy in the shallow crust in TMR. Our results show that, although the thickness of the upper crustal layer is only about 1 km, the layer shows a strong anisotropy with magnitude of 8~15%; in the deeper of crust, the magnitudes of anisotropy is about 3%~5%, showing a pattern with fast-symmetric-axis. The crust anisotropy beneath TMR in North China obtained in this study also shows a significant difference in both the lateral and vertical scale, which might imply a regional anisotropy characteristic in the studied region.

地震层析成像与天山地区的壳幔深部结构研究

The Tien Shan is the most prominent intracontinental mountain belt on the earth. The active crustal deformation and earthquake activities provide an excellent place to study the continental geodynamics of intracontinental mountain belt. The studies of deep structures in crust and upper mantle are significantly meaningful for understanding the geological evolution and geodynamics of global intracontinental mountain belts. This dissertation focuses on the deep structures and geodynamics in the crust and upper mantle in the Tien Shan mountain belt. With the arrival time data from permanent and temporal seismic stations located in the western and central Tien Shan, using seismic travel time tomographic method, we inversed the P-wave velocity and Vp/Vs structures in the crust and uppermost mantle, the Pn and Sn velocities and Pn anisotropic structures in the uppermost mantle, and the P-wave velocity structures in the crust and mantle deep to 690km depth beneath the Tien Shan. The tomographic results suggest that the deep structures and geodynamics have significant impacts not only on the deformations and earthquake activities in the crust, but also on the mountain building, collision, and dynamics of the whole Tien Shan mountain belt. With the strongly collision and deformations in the crust, the 3-D P-wave velocity and Vp/Vs ratio structures are highly complex. The Pn and Sn velocities in the uppermost mantle beneath the Tien Shan, specially beneath the central Tien Shan, are significantly lower than the seismic wavespeed beneath geological stable regions. We infer that the hot upper mantle from the small-scale convection could elevate the temperature in the lower crust and uppermost mantle, and partially melt the materials in the lower crust. The observations of low P-wave and S-wave velocities, high Vp/Vs ratios near the Moho and the absences of earthquake activities in the lower crust are consistent with this inference. Based on teleseismic tomography images of the upper mantle beneath the Tien Shan, we infer that the lithosphere beneath the Tarim basin has subducted under the Tien Shan to depths as great as 500 km. The lithosphere beneath the Kazakh shield may have subducted to similar depths in the opposite direction, but the limited resolution of this data set makes this inference less certain. These images support the plate boundary model of converge for the Tien Shan, as the lithospheres to the north and south of the range both appear to behave as plates.

岩石中地震波速度和衰减的应力变化响应研究

Stress change is one of key factors in seismic nucleating and triggering; therefore for understanding and forecasting earthquakes, it is necessary to research on stress status and its changes in rocks. Propagating in underground structures, wave velocity and attenuation contain information on stress changes of the Earth’s interior. For a better understanding of relationship between seismic data and stress changes, modeling and ultrasonic test supply significant references. In this article, acoustoelastic theory is introduced to explain nonlinear elastic characteristics of rocks. Based on the acoustoelastic theory, a solid-fluid coupled model is given to calculate velocity under different stress for porous and liquid fulfilled rocks. Except for the stress-velocity relationship, effects of pore pressure induced stress changes on ultrasonic coda attenuation are also studied. Intrinsic attenuation quality factors are calculated for a comparison purpose. Finally, the relationship between elastic constants and stress changes is thoroughly investigated, a mixture model from two phases of Hooke media is introduced to explain the differences between dynamic and static moduli, a relation among wave length, wave velocities and elastic moduli considering dimension of microstructure, dimension and state of surface between phases is presented. The most important aspect of this work is exploring and establishing relationships between the seismic properties of rocks and changes of their stress conditions, which will have its application in earthquake forecast and seismic hazard.

消除采集脚印的处理方法研究

Seismic technique is in the leading position for discovering oil and gas trap and searching for reserves throughout the course of oil and gas exploration. It needs high quality of seismic processed data, not only required exact spatial position, but also the true information of amplitude and AVO attribute and velocity. Acquisition footprint has an impact on highly precision and best quality of imaging and analysis of AVO attribute and velocity. Acquisition footprint is a new conception of describing seismic noise in 3-D exploration. It is not easy to understand the acquisition footprint. This paper begins with forward modeling seismic data from the simple sound wave model, then processes it and discusses the cause for producing the acquisition footprint. It agreed that the recording geometry is the main cause which leads to the distribution asymmetry of coverage and offset and azimuth in different grid cells. It summarizes the characters and description methods and analysis acquisition footprint’s influence on data geology interpretation and the analysis of seismic attribute and velocity. The data reconstruct based on Fourier transform is the main method at present for non uniform data interpolation and extrapolate, but this method always is an inverse problem with bad condition. Tikhonov regularization strategy which includes a priori information on class of solution in search can reduce the computation difficulty duo to discrete kernel condition disadvantage and scarcity of the number of observations. The method is quiet statistical, which does not require the selection of regularization parameter; and hence it has appropriate inversion coefficient. The result of programming and tentat-ive calculation verifies the acquisition footprint can be removed through prestack data reconstruct. This paper applies migration to the processing method of removing the acquisition footprint. The fundamental principle and algorithms are surveyed, seismic traces are weighted according to the area which occupied by seismic trace in different source-receiver distances. Adopting grid method in stead of accounting the area of Voroni map can reduce difficulty of calculation the weight. The result of processing the model data and actual seismic demonstrate, incorporating a weighting scheme based on the relative area that is associated with each input trace with respect to its neighbors acts to minimize the artifacts caused by irregular acquisition geometry.

KEL-TI介质的地震波场响应与成像方法研究

To deal with the problems in multi-component converted seismic wave exploration in coal fields, the wave propagating features and imaging methods of multi-component converted waves in coal measure strata are researched in this thesis firstly. The relations between viscoelasticity and anisotropy in coal measure strata are analyzed to build KEL-TI model, and which seismic wave propagating and attenuating features are researched. The disadvantages of converted wave imaging methods based on common converted point gather are analyzed and constant velocity no NMO converted wave imaging method based on common scattering point gather is put forward, according to Huygens-Fresnel principle, which applicabilities in the elastic isotropic, elastic TI and KEL-TI situations are discussed. To different model simulation data, the common scattering point gathers’ and stacked profiles’ features are analyzed. The results show that the method can image compressional waves and converted waves with high precision. Secondly, the resolution enhancing theories and methods of converted wave are researched by Rayleigh wave suppressing, converted wave static correction and poststack inverse-Q filtering. 1) The polarization filter is designed by the instantaneous polarization information of seismic waves, and the Rayleigh wave suppressing method is researched. From the spectrum analysis before and after filtering, it can be derived that the amplitudes are kept relatively. 2) To constant velocity no NMO converted wave imaging method, the static correction method based on common equivalent offset point gather is put forward and tested to the actual converted waves. 3) The relation between equivalent quality factor of converted wave, compressional wave quality factor and the ratio of compressional to shear wave velocity is derived. The compressional wave quality inversion method by first arrivals of none-offset VSP is researched, and which is then transformed to the equivalent quality factor to perform inverse-Q filtering of actual converted waves. The result has shown that the method can recover the high frequency energy of converted waves. At last, the theories and methods researched in this thesis are practiced to the 3D3C seismic exploration in Guqiao coal mine in Huainan and achieve good results.

煤系地层三维三分量地震数据解释方法研究

Conventional 3D seismic exploration cannot meet the demand of high yield and high efficiency safe production in coal mine any more. Now it is urgent to improve the discovery degree of coal mine geological structures for coal production in China. Based on 3D3C seismic exploration data, multi-component seismic information is fully excavated. First systematic research on 3D3C seismic data interpretation of coal measure strata is carried out. Firstly, by analyzing the coal measure strata, the seismic-geologic model of coal measure strata is built. Shear wave logging is built by using regression analysis. Horizon calibration methods of PP-wave and PS-wave are studied and the multi-wave data are used together to interpret small faults. Using main amplitude analysis technology, small faults which cannot be found from PP-wave sections can be interpreted from the low frequency PS-wave sections. Thus, the purpose to applying PS-wave data to fine structure assistant interpretation is achieved. Secondly, PP- and PS-wave post-stack well constrained inversion methods of coal measure strata are studied. Joint PP- and PS-wave post-stack inversion flow is established. More attribute parameters, which are applied in fine lithology interpretation of coal measure strata, are obtained from combinations of the inversion results. Exploring the relation between rock with negative Poisson’s ratio and anisotropy, fracture development in coal seam are predicted. Petrophysical features of coal measure strata are studied, and the relations between elastic parameters and lithology, fluid and physical properties are established. Inversions of the physical parameters such as porosity, permeability and water saturation, which reflect lithology and fluid property, are obtained. Finally, the approaches of shear wave splitting and Thomsen parameters inversion, which provide new ideas for seismic anisotropy interpretation of coal measure strata, are studied to predict fracture development. The results of practical application indicate that the methods in this paper have good feasibility and applicability. They have positive significance for high yield and high efficiency safe production in coal mine.

有源电磁扩散场克希霍夫积分偏移

As active electromagnetic method, field data of CSAMT method follow the equation of diffusion. Propagting in solid earth media, diffusion EM signal has strong attenuation and dispersion, otherwise seismic wave shows weak attenuation and dispersion, therefore the resolution power of CSAMT method is not better than seismic reflection method. However, there is consistence and similarity between EM signal and seismic wave in wave equation, we can apply Kirchhoff integral migration technique, a proven one in seismic method in time domain, to carry out seduo-seismic processing for CSAMT signal in frequency domain so that the attenuation and dispersion could be made compensated in some extent, and the resolution power and interpretation precision of active EM wave could be improved. Satisfying passive homogeneous Helmholtz quation, we proceed with Green theorem and combine the active inhomogenous Helmholtz quation, the Kirchhoff integral formula could be derived. Given practical problems, if we only consider the surface integral value, and assume that the intergral value in other interface is zero, combined with Green theorem in uniform half space, the expression could be simplified, and we can obtain frequency-domain Kirchhoff integral formula in surface, which is also called downward continuation of EM field in frequency domain. With image conditions and energy compensation considered, in order to get image conditions in time domain Fourier inverse transformation in frequency domain can be performed, so we can formulate the active Kirchhoff integral migration expression. At first, we construct relative stratified model, with different frequency series taken into account, then we change the distances between transmitter and reciever, the EM response can be obtained. Analyzing the EM properties, we can clarify near and far zone that can instruct us to carry out transmitter layout in practical application. Combined with field data surveyed in far zone, We perform Kirchhoff integral migration and compare the results with model to interpret. Secondly, with far field EM data, we apply TM mode to get EM response of given 2D model, then apply Kirchhoff integral migration on modelling data and interpret the results.

地震波模拟：非规则网格离散的数值解法、吸收边界及其实际应用

This dissertation presents a series of irregular-grid based numerical technique for modeling seismic wave propagation in heterogeneous media. The study involves the generation of the irregular numerical mesh corresponding to the irregular grid scheme, the discretized version of motion equations under the unstructured mesh, and irregular-grid absorbing boundary conditions. The resulting numerical technique has been used in generating the synthetic data sets on the realistic complex geologic models that can examine the migration schemes. The motion equation discretization and modeling are based on Grid Method. The key idea is to use the integral equilibrium principle to replace the operator at each grid in Finite Difference scheme and variational formulation in Finite Element Method. The irregular grids of complex geologic model is generated by the Paving Method, which allow varying grid spacing according to meshing constraints. The grids have great quality at domain boundaries and contain equal quantities of nodes at interfaces, which avoids the interpolation of parameters and variables. The irregular grid absorbing boundary conditions is developed by extending the Perfectly Matched Layer method to the rotated local coordinates. The splitted PML equations of the first-order system is derived by using integral equilibrium principle. The proposed scheme can build PML boundary of arbitrary geometry in the computational domain, avoiding the special treatment at corners in a standard PML method and saving considerable memory and computation cost. The numerical implementation demonstrates the desired qualities of irregular grid based modeling technique. In particular, (1) smaller memory requirements and computational time are needed by changing the grid spacing according to local velocity; (2) Arbitrary surfaces and interface topographies are described accurately, thus removing the artificial reflection resulting from the stair approximation of the curved or dipping interfaces; (3) computational domain is significantly reduced by flexibly building the curved artificial boundaries using the irregular-grid absorbing boundary conditions. The proposed irregular grid approach is apply to reverse time migration as the extrapolation algorithm. It can discretize the smoothed velocity model by irregular grid of variable scale, which contributes to reduce the computation cost. The topography. It can also handle data set of arbitrary topography and no field correction is needed.