地震波谱元法数值模拟和偏移研究

With the development of both seismic theory and computer technology, numerical modeling technology of seismic wave has achieved great advancement during the past half century. The current methods under development include finite differentiation method (FDM), finite element method (FEM), pseudospectral method (PSM), integral equation method (IEM) and spectral element method (SEM). They exert their very important roles in every corner of seismology and seismic prospecting. Large quantity of researches towards spectral element method in the end of last century bring this method to a new era, which results in perfect solution of many difficult problems. However, parts of posterior works such as seismic migration and inversion which base on spectral element method have never been studied widely at least up to the present whereas are of importance to seismic imaging and seismic wave propagation. Based on previous work, this paper uses spectral element method to investigate the characteristics and laws of the seismic wave propagation in isotropic and anisotropic media. By thoroughly studying this high-accuracy method, we implement a kind of reverse-time pre- and post-stack migration based on SEM. In order to verify the validity of the SEM method, we have simulated the propagation of seismic wave in several different models. The simulation results show that: (1) spectral element method can be used to model any complex models and the computational results are comparable with the expected results and the analytic results; (2) the optimum accuracy can be achieved when the rank is between 4 and 9. When it is below 4, the dispersion may occur; and when it is above 9, the time step-length will be changed accordingly with the reducing space step-length in order to keep the computation stability. This will exponentially increase the computation time and at the same time the memory even if simulating the same media. This paper also applies explosive reflection surface imaging technology, time constancy principle of wave-filed extrapolation and least travetime raytracing technology of surface source to SEM pre- and post-stack migration of isotropic and anisotropic media. All imaging results derived by the above methods agree well with the real geological models and the position of interface and inflexions can also return to their right location well. This indicates that the method proposed in this paper is a kind of technology with high accuracy and robust stability. It can serve as an alternative method in real seismic data processing. All these work can boost the development of high-accuracy seismic imaging, and therefore have significant inference value.

鄂尔多斯盆地环县地区长6-长8沉积相及储层研究

The exploration in recent years shows that the Yanchang Formation in the southwest of Ordos Basin is of great resource potential and good exploration and exploitation prospect. In the thesis ，sedimentary source analysis，sedimentary system，sedimentary microfacies，sandstones distribution and reservoir characteristic are studied and favorable oil area are forecasted in Chang6-Chang8 of Yanchang formation in HuanXian region， by mainly study on the data of field section observation ，core observation, well logging explaination and routine microscope slice identification，scanning Electron Microscope and reservoir analysis of lithology and physical property , Under the guidance of such advanced theories and methods as sedimentology，reservoir sedimentology，lithological oil pool and so on. The stratum of Chang6-Chang8 of Yanchang formation could be divided into pieces of member following the principles that firstly contrasting the big segments, then contrasting the small segments, being controlled by cycle and consulting the thickness etc.And the characteristic of stratum are detailed discussed , respectively. Based on the source direction of the central basin, heavy and light minerals are used to analyse source direction of Chang6 and Chang8 member, in HuanXian area. Research result shows that the source of Chang6 and Chang8 member is mixed provenance,including west-south,west and east-north. By the study of rock types、 sedimentary conformation、lithology and electromotive curve combination and palaeo-biology,lake、delta and braided delta mianly developed in study area are recognized, Subaqueous distributary channels in delta front and in braided delta front, and sand body in deep-lake turbidite, are the main reservoir.forthermore,the characteristic of depositional system and sandy body in space are discussed. Applied with routine microscope slice identification, Scanning Electron Microscope, reservoir lithology and physical property analysis and other analytic machinery, Feldspar-lithic fine-sandstone and feldspar fine-sandstone are mainly sandstone of Y Chang6-Chang8 in Huanxian area, small pore and tiny pore are the main pore types, tiny throat type and micro-fine throat type are widely developed , secondary dissolution porosity, intercrystal porosity, tiny pore and micro-crack are main pore types.Intergranular porosity and dissolution porosity secondary is the main pore secondary. The dominant diagenesis types in the area are compaction, cementation, replacement and dissolution. Chlorite films cementation facies, carbonate cementation facies ,mud cementation compaction facie, compaction 、pressure solution facies are the main diagenetic facies,in which Chlorite films cementation facies is the best diagenetic facies in study area. Reservoir influence factor analysis ,rock types are the main factor forming this low-pore and low-permeability of Chang6-Chang8 member in study area,and relatively higher permeability area are cortrolled by sedimentary facies distribution, diagenesis improved reservoir physical property. According to the distributing of sedimentary micro-facies and sandy body , and the test oil, favorable region in Chang6-Chang8 are forecasted.

各向同性与各向异性介质中地震波旅行时正演模拟与反演成像

By seismic tomography, interesting results have been achieved not only in the research of the geosphere with a large scale but also in the exploration of resources and projects with a small scale since 80'. Compared with traditional inversion methods, seismic tomography can offer more and detailed information about subsurface and has been being paid attention by more and more geophysicists. Since inversion based on forward modeling, we have studied and improved the methods to calculate seismic traveltimes and raypaths in isotropic and anisotropic media, and applied the improved forward methods to traveltime tomography. There are three main kinds of methods to calculate seismic traveltime field and its ray path distribution, which are ray-tracing theory, eikonal equation by the finite-difference and minimum traveltime tree algorithm. In ray tracing, five methods are introduced in the paper, including analytic ray tracing, ray shooting, ray bending, grid ray tracing and rectangle grid ray perturbation with three points. Finite-difference solution of eikonal equation is very efficient in calculation of seismic first-break, but is awkward in calculation of reflection traveltimes. We have put forward a idea to calculate traveltimes of reflected waves using a combining way of eikonal equation method and other one in order to improve its capability of dealing with reflection waves. The minimum traveltime tree algorithm has been studied with emphases. Three improved algorithms are put forward on the basis of basic algorithm of the minimum traveltime tree. The first improved algorithm is called raypath tracing backward minimum traveltime algorithm, in which not only wavelets from the current source but also wavelets from upper source points are all calculated. The algorithm can obviously improve the speed of calculating traveltimes and raypaths in layered or blocked homogeneous media and keep good accuracy. The second improved algorithm is raypath key point minimum traveltime algorithm in which traveltimes and raypaths are calculated with a view of key points of raypaths (key points of raypths mean the pivotal points which determine raypaths). The raypath key point method is developed on the basis of the first improved algorithm, and has better applicability. For example, it is very efficient even for inhomogeneous media. Another improved algorithm, double grid minimum traveltime tree algorithm, bases upon raypath key point scheme, in which a model is divided with two kinds of grids so that the unnecessary calculation can be left out. Violent undulation of curved interface often results in the phenomenon that there are no reflection points on some parts of interfaces where there should be. One efficacious scheme that curved interfaces are divided into segments, and these segments are treated respectively is presented to solve the problem. In addition, the approximation to interfaces with discrete grids leads to large errors in calculation of traveltimes and raypaths. Noting the point, we have thought a new method to remove the negative effect of mesh and to improve calculation accuracy by correcting the traveltimes with a little of additional calculation, and obtained better results.

均质弹性条件下多洞平面问题解析算法及其应用研究

The scholars in the world have been trying to find an effective analytic algorithm of multiple hole problems usually meet in engineering designs. Though some studies on circular or elliptic holes had been achieved under specific conditions, no efforts were made to any multiple hole problems that is most significant for engineering designs. The author has made further studies on any multiple hole problems, using complex variable function method and Schwarz alternating method. After solving a series of technological difficulties, the author obtains an effective analytic algorithm, and acquires stress field and displacement field with high accuracy, which can be conducted for arbitrary many iterations according to practical accuracy requirements. In addition, th solution of stress and displacement fields, even for multiple holes of complex shapes and smaller distances. Further, the author made preliminary studies on viscoelastic displacement solution for any double holes. In terms of the obtained displacement solution of any multiple holes, this paper studies displacement back-analysis for the excavations of two tunnels, and find that the back-analysis method is accurate. Additionally, the author presents the mathematical prove of inversion uniqueness for ground stresses, elastic modulus and Poisson ratio. The author believes that the accurate analytic algorithm provided in this paper will presents an effective way to stress and displacement analysis for any multiple hole problems, optimal arrangement of multiple holes, hole shape optimization of multiple holes, etc..