构造复杂区变速成图方法研究

In the complex structure areas, velocity field building and structure mapping are important for seismic exploration. With the development of seismic exploration, the methods of structure mapping, reservoir prediction and reservoir description all require high precious velocity field. And more accurate depth-structure maps are required for well site design. Aiming at the problems and defects in velocity analysis and structure mapping in oil seismic exploration, the paper which is based on the studies of real data in several areas combines the theories with practical application, and analyzes the precision and applicability of several methods of velocity model building. After that, the following methods are mainly studied: the coherence inversion methods based on the pre-stack CMP gathers or stacking velocity; the interval velocity inversion methods constrained by multi-well; the Random Simulation method; 3D Image Ray Map Migration method and the structure mapping in floating datum and in fixed datum, and then we conclude the method of building high precious seismic velocity field and structure mapping with variable velocity. Firstly, the paper analyses the distributing rule of the velocity variation in the areas with complex structures in the northwest of China, then points out that velocity is a crucial factor which influences the precision of structure mapping, and the velocity variations have something to do with the shapes of the structures, the variety of lithology and so on. The key point of improving the precision of seismic velocity field is to obtain a structure mapping with high precision. We also describe the range and conditions of these methods. Secondly, by comparing many popular methods of velocity model building, we propose a new method in the use of velocity model building. The new method is more effective in velocity model building under every kind of complex condition and is worthy of spreading. At last, the paper fingers out that it is a system engineering to study variable velocity mapping in every kind of complex structure areas. Every step of the work can affect the final results. So it is important to build high efficient and practical velocity model and the flows of mapping processing. The paper builds the flows and gives some examples. The method has been applied in more than ten exploring surveys. The application proves that this method could bring good effect on researching on low-amplitude trap, reservoir prediction, reservoir description and the integrated research of oil&gas geology. Keywords: structure mapping velocity model building complex structure variable velocity media

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.

大陆地壳岩石结构重建：方法与实例

The composition of the continental crust has long been a subject of interest to earth scientists as it can provide key information about the crustal growth and evolution of the continents. In this paper we make a comparative study on the lithological discrimination schemes featuring with the use of different seismic attributes, such as P-wave velocity, P- to S-wave velocity ratio, acoustic or elastic impedances, Lame impedances and high-sensitive identification factors. The results demonstrate that Lame impedances have more powerful constrains than other seismic attributes. In order to fully take the advantage of make the best of the different seismic response of crustal rock, we firstly use seismic attribute that have weak distinguish power to construct loose constrained lithological model, then use seismic attributes that have stronger distinguish power to tighten the constrains of lithological discrimination. We propose a joint scheme (chain constrain technique) by combing all available constrains to reduce the non-uniqueness in mapping rock distribution. We adopt chain constrain technique to construct lithological model beneath Tunxi-Wenzhou transect, Southeastern China, Manzhouli-Suifenhe transect, Northeastern China, and geophysical profile in Bohai Bay Basin, North China. The results can be suumarized as the follows: (1) We compare the sensitivity of different seismic factor constraints on rock types, and conclude that Lame impedances have tighter constrains than seismic velocity, Vp/Vs, density. (2) We propose chain constrains to construct lithological model from integrated geophysical data, and reduce the non-uniqueness in mapping rock distribution. (3) We reconstruct crustal lithological model beneath Tunxi-Wenzhou transect, Southeastern China. The results suggested that Jiangshan-Shaoxing fault is a crust-scale, and it is the boundary between Cathaysia and Yanthze blocks. (4) We construct crustal lithological model beneath Manzhouli-Suifenhe transect, Northeastern China. (5) We map the petrologic distribution along a geophysical profile in Bohai Bay Basin, North China, and construct a three-layered petrology model from the depth 2 km to about 10 km, consisted of basalt (the first layer), pelitic siltstone (the second layer), and silty mudstone and fine sandstone (the third layer).

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

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.

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

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.

云南岩石圈流变学及多震层孕震环境

The Study on rheology of the lithosphere and the environments of the seismogenic layer is currently the basic project of the international earthquake research. Yunnan is the ideal place for studying this project. Through the multi-disciplinary comprehensive study of petrology, geophysics, seismo-geology, rock mechanics, etc., the depth-strength profiles of the lithosphere have been firstly constructed, and the seismogenic layer and its geophysical and tectonic environments in Yunnan have been systematically expounded in this paper. The related results achieved are of the important significances for further understanding the mechanism of strong earthquake generation, dividing the potential foci and exposing recent geodynamical processes in Yunnan. Through the comprehensive contrast of the metamorphic rocks in early and middle Proterozoic outcropping on the surface, DSS data and experimental data of rock seismic velocity under high temperature and high pressure, the petrological structure of the crust and upper mantle has been studied on Yunnan: the upper, middle and lower crust is composed of the metamorphic rocks of greenschist, amphibolite and granulite facies, respectively or granitoids, diorites and gabbros, respectively, and the upper mantle composed of the peridotites. Through the contrast studies of the heat flow and epicenters of the strong earthquakes, the distribution of the geotemperature and the data of focal depth, the relationship of between seismicity and geothermal structure of the lithosphere in Yunnan has been studied: the strong earthquakes with magnitude M ≥ 6.0 mainly take place at the geothermal gradient zone, and the seismic foci densely distribute between 200～500 ℃ isogeotherms. On the basis of studies of the rock properties and constituents of the crust and upper mantle and geothermal structure of the lithosphere, the structure of the rheological stratification of the lithosphere has been studied, and the corresponding depth-strength profiles have been constructed in Yunnan. The lithosphere in majority region of Yunnan has the structure of the rheological stratification, i.e. the brittle regime in the upper crust or upper part of the upper crust, ductile regime in the middle crust or lower part of the upper crust to middle crust, ductile regime in the lower crust and ductile regime in the subcrustal lithosphere. The rheological stratification has the quite marked lateral variations in the various tectonic units. The distributions of the seismogenic layer have been determined by using the high accurate data of focal depth. Through the contrast of the petrological structure, the structure of seismic velocity, electric structure, geotemperature structure, and rheological structure and the study of the focal mechanism in the seismogenic layer, the geophysical environments of the seismogenic layer in Yunnan have been studied. The seismogenic layer in Yunnan is located at the depths of 3 ～ 20 km; the rocks in the seismogenic layer are composed of the metamorphic rocks of greenschist to amphibolite facies (or granites to diorites); the seismogenic layer and its internal focal regions of strong earthquakes have the structure of medium properties with the relatively high seismic velocity, high density and high resistivity; there exists the intracrustal low seismic velocity and high conductivity layer bellow the seismogenic layer, the geotemperature is generally 100～500 ℃ in the depth range in which the seismogenic layer is located. The horizontal stress field predominates in the seismogenic layer, the seismogenic layer corresponds to the brittle regime of the upper crust or brittle regime of the upper crust to semibrittle regime of the middle crust. The formation of the seismogenic layer, preparedness and occurrence of the strong earthquakes is the result of the comprehensive actions of the source fault, rock constituent, structure of the medium properties, distribution of the geotemperature, rheological structure of the seismogenic layer and its external environments. Through the study of the structure, active nature, slip rate, segmentation of the active faults, and seismogenic faults, the tectonic environments of the seismogenic layer in Yunnan have been studied. The source faults of the seismogenic layer in Yunnan are mainly A-type ones and embody mainly the strike slip faults with high dip angle. the source faults are the right-lateral strike slip ones with NW-NNW trend and left-lateral strike slip ones with NE-NEE trend in Southwestern Yunnan, the right-lateral strike slip ones with NNW trend and left-lateral strike slip ones with NNE trend (partially normal ones) in Northwestern Yunnan, the right-lateral strike slip ones with NWW trend in Central Yunnan and left-lateral strike slip ones with NW-NNW trend in Eastern Yunnan. Taking Lijiang earthquake with Ms = 7.0 for example. The generating environments of the strong earthquake and seismogenic mechanical mechanism have been studied: the source region of the strong earthquake has the media structure with the relatively high seismic velocity and high resistivity, there exists the intracrustal low velocity and high conductivity layer bellow it and the strong earthquakes occur near the transitional zone of the crustal brittle to ductile deformation. These characteristics are the generality of the generating environments of strong earthquakes. However, the specific seismogenic tectonic environments and action of the stress field of the seismic source in the various regions, correspondingly constrains the dislocation and rupture mechanical mechanism of source fault of strong earthquake.

内蒙古东部喀喇沁地区早中生代底侵作用

Mafic granulite xenoliths have been extensively concerned over the recent years because they are critical not only to studies of composition and evolution of the deep parts of continental crust but to understanding of the crust-mantle interaction. Detailed petrology, geochemistry and isotope geochronology of the Early Mesozoic mafic-ultramafic cumulate xenoliths and mafic granulite xenoliths and their host diorites from Harqin area, eastern Inner-Mongolia have been studied here. Systematic Rb-Sr isochron, ~(40)Ar-~(39)Ar and K-Ar datings for mafic-ultramafic cumulate xenoliths give ages ranging from 237Ma to 221Ma. Geochemical research and forming temperature and pressure estimates suggest that cumulates are products of the Early Mesozoic mantle-derived magmatic underplating and they formed in the magmatic ponds at the lowermost of the continental crust and are later enclaved by the dioritic magma. Detailed study on the first-discovered mafic granulite xenoliths reveals that their modal composition, mineral chemistry and metamorphic P-T conditions are all different from those of the Precambrian granulite exposed on the earth surface of the North China craton. High-resolution zircon U-Pb dating suggests that the granulite facies metamorphism may take place in 253 ～ 236Ma. Hypersthene single mineral K-Ar dating gives an age of 229Ma, which is believed to represent a cooling age of the granulite. As the host rock of the cumulate and granulite xenoliths, diorites intruded into Archean metamorphic rocks and Permian granite. They are mainly composed of grandodiorite, tonalite and monzogranite and show metaluminous and calc-alkaline features. Whole rock and single mineral K-Ar dating yields age of 221 ～ 223Ma, suggesting a rapid uplift in the forming process of the diorites. Detailed field investigation and geochemical characteristics indicate that these diorites with different rock types are comagmatic rocks, and they have no genetic correlation with cumulate and granulite xenoliths. Geochemical model simulating demonstrates that these diorites in different lithologies are products of highly partial melting of Archean amphibolite. It is considered that the Early Mesozoic underplating induced the intrusion of diorites, and it reflects an extensional geotectonic setting. Compression wave velocity V_P have been measured on 10 representative rock samples from the Early Mesozoic granulite and mafic-ultramafic cumulate xenoliths population as an aid to interpret in-situ seismic velocity data and investigating velocity variation with depth in a mafic lower crust. The experiments have been carried out at constant confining pressures up to 1000MPa and temperatures ranging from 20 ℃ to around 1300 ℃, using the ultrasonic transmission technique. After corrections for estimated in situ crustal pressures and temperatures, elastic wave velocities range from 6.5 ～ 7.4 km s~(-1). On the basis of these experimental data, the Early-Mesozoic continental compression velocity profile has also been reestablished and compared with those of the present and of the different tectonic environments in the world. The result shows that it is similar to the velocity structure of the extensional tectonic area, providing new constraints on the Early Mesozoic continental structure and tectonic evolution of the North-China craton. Combining with some newly advancements about the regional geology, the thesis further proposes some constraints on the Mesozoic geotectonic evolution history, especially the features of deep geology of the North China craton.

泥岩裂缝储层地震识别的理论和方法

During the exploration of fractured reservoirs, worldwide difficult problems will be encountered: how to locate the fractured zones, how to quantitatively determine the azimuth, density, and distribution of the fractures, and how to compute the permeability and porosity of the fractures. In an endeavor to solve these problems, the fractured shale reservoir in SiKou area of ShengLi oil field was chosen as a study area. A study of seismic predictive theory and methods for solving problems encountered in fractured reservoir exploration are examined herein. Building on widely used current fractured reservoir exploration techniques, new seismic theories and methods focusing on wave propagation principles in anisotropic medium are proposed. Additionally, integrated new seismic data acquisition and processing methods are proposed. Based on research and application of RVA and WA methods from earlier research, a new method of acoustic impedance varying with azimuth (IPVA) creatively is put forth. Lastly combining drilling data, well log data, and geologic data, an integrated seismic predictive method for cracked reservoir bed was formed. A summary of the six parts of research work of this paper is outlined below. In part one, conventional geologic and geophysical prediction methods etc. for cracked reservoir exploration are examined, and the weaknesses of these approaches discussed. In part two, seismic wave propagation principles in cracked reservoirs are studied. The wave equation of seismic velocity and attenuation factor in three kinds of fracture mediums is induced, and the azimuth anisotropy of velocity and attenuation in fracture mediums is determined. In part three, building on the research and application of AVA and WA methods by a former researcher, a new method of acoustic impedance creatively varying with azimuth (IPVA) is introduced. A practical software package utilizing this technique is also introduced. In part four, Base on previously discussed theory, first a large full azimuth 3d seismic data (70km~2) was designed and acquired. Next, the volume was processed with conventional processing sequence. Then AVA, WA, and IPVA processing was applied, and finally the azimuth and density of the fractures were quantitatively determined by an integrated method. Predictions were supported by well data that indicate the approach is highly reliable. in part five, geological conditions contributing to cracked reservoir bed formation are analyzed in the LuoJia area resulting in the discovery that the main fractured zones are related to fault distribution in the basin, that also control the accumulation of the oil and gas, the generation mechanisms and types of fractured shale reservoirs are studied. Lastly, by using full 3D seismic attributes, azimuth and density of cracked reservoir zones are successfully quantitative predicted. Using an integrated approach that incorporates seismic, geologic and well log data, the best two fractured oil prospects in LouJia area are proposed. These results herein represent a break through in seismic technology, integrated seismic predictive theory, and production technology for fractured reservoirs. The approach fills a void that can be applied both inside China, and internationally. Importantly, this technique opens a new exploration play in the ShengLi oil field that while difficult has substantial potential. Properly applied, this approach could play an important role toward stabilizing the oil field' production. In addition, this technique could be extended fracture exploration in other oil fields producing substantial economic reward.

叠前深度偏移速度模型的建立方法研究

Seismic exploration is the main tools of exploration for petroleum. as the society needs more petroleum and the level of exploration is going up, the exploration in the area of complex geology construction is the main task in oil industry, so the seismic prestack depth migration appeared, it has good ability for complex construction imaging. Its result depends on the velocity model strongly. So for seismic prestack depth migration has become the main research area. In this thesis the difference in seismic prestack depth migration between our country and the abroad has been analyzed in system. the tomographical method with no layer velocity model, the residual curve velocity analysical method based on velocity model and the deleting method in pre-processing have been developed. In the thesis, the tomographysical method in velocity analysis is been analyzed at first. It characterized with perfection in theory and diffculity in application. This method use the picked first arrivial, compare the difference between the picked first arrival and the calculated arrival in theory velocity model, and then anti-projected the difference along the ray path to get the new velocity model. This method only has the hypothesis of high frequency, no other hypothesis. So it is very effective and has high efficiency. But this method has default still. The picking of first arrival is difficult in the prestack data. The reasons are the ratio of signal to noise is very low and many other event cross each other in prestack data. These phenomenon appear strongly in the complex geology construction area. Based on these a new tomophysical methos in velocity analysis with no layer velocity model is been developed. The aim is to solve the picking problem. It do not need picking the event time contiunely. You can picking in random depending on the reliability. This methos not only need the pick time as the routine tomographysical mehtod, but also the slope of event. In this methos we use the high slope analysis method to improve the precision of picking. In addition we also make research on the residual curve velocity analysis and find that its application is not good and the efficiency is low. The reasons is that the hypothesis is rigid and it is a local optimizing method, it can solve seismic velocity problem in the area with laterical strong velocity variation. A new method is developed to improve the precision of velocity model building . So far the pattern of seismic prestack depth migration is the same as it aborad. Before the work of velocity building the original seismic data must been corrected on a datum plane, and then to make the prestack depth migration work. As we know the successful example is in Mexico bay. It characterized with the simple surface layer construction, the pre-precessing is very simple and its precision is very high. But in our country the main seismic work is in land, the surface layer is very complex, in some area the error of pre-precessing is big, it affect the velocity building. So based on this a new method is developed to delete the per-precessing error and improve the precision of velocity model building. Our main work is, (1) developing a effective tomographical velocity building method with no layer velocity model. (2) a new high resolution slope analysis method is developed. (3) developing a global optimized residual curve velocity buliding method based on velocity model. (4) a effective method of deleting the pre-precessing error is developing. All the method as listed above has been ceritified by the theorical calculation and the actual seismic data.

用宽频带远震体波波形研究中国及邻近地区地壳构造及地质意义初探

The theory and approach of the broadband teleseismic body waveform inversion are expatiated in this paper, and the defining the crust structure's methods are developed. Based on the teleseismic P-wave data, the theoretic image of the P-wave radical component is calculated via the convolution of the teleseismic P-wave vertical component and the transform function, and thereby a P-wavefrom inversion method is built. The applied results show the approach effective, stable and its resolution high. The exact and reliable teleseismic P waveforms recorded by CDSN and IRIS and its geodynamics are utilized to obtain China and its vicinage lithospheric transfer functions, this region ithospheric structure is inverted through the inversion of reliable transfer functions, the new knowledge about the deep structure of China and its vicinage is obtained, and the reliable seismological evidence is provided to reveal the geodynamic evolution processes and set up the continental collisional theory. The major studies are as follows: Two important methods to study crustal and upper mantle structure -- body wave travel-time inversion and waveform modeling are reviewed systematically. Based on ray theory, travel-time inversion is characterized by simplicity, crustal and upper mantle velocity model can be obtained by using 1-D travel-time inversion preliminary, which introduces the reference model for studying focal location, focal mechanism, and fine structure of crustal and upper mantle. The large-scale lateral inhomogeneity of crustal and upper mantle can be obtained by three-dimensional t ravel-time seismic tomography. Based on elastic dynamics, through the fitting between theoretical seismogram and observed seismogram, waveform modeling can interpret the detail waveform and further uncover one-dimensional fine structure and lateral variation of crustal and upper mantle, especially the media characteristics of singular zones of ray. Whatever travel-time inversion and waveform modeling is supposed under certain approximate conditions, with respective advantages and disadvantages, and provide convincing structure information for elucidating physical and chemical features and geodynamic processes of crustal and upper mantle. Because the direct wave, surface wave, and refraction wave have lower resolution in investigating seismic velocity transitional zone, which is inadequate to study seismic discontinuities. On the contrary, both the converse and reflected wave, which sample the discontinuities directly, must be carefully picked up from seismogram to constrain the velocity transitional zones. Not only can the converse wave and reflected wave study the crustal structure, but also investigate the upper mantle discontinuities. There are a number of global and regional seismic discontinuities in the crustal and upper mantle, which plays a significant role in understanding physical and chemical properties and geodynamic processes of crustal and upper mantle. The broadband teleseismic P waveform inversion is studied particularly. The teleseismic P waveforms contain a lot of information related to source time function, near-source structure, propagation effect through the mantle, receiver structure, and instrument response, receiver function is isolated form teleseismic P waveform through the vector rotation of horizontal components into ray direction and the deconvolution of vertical component from the radial and tangential components of ground motion, the resulting time series is dominated by local receiver structure effect, and is hardly irrelevant to source and deep mantle effects. Receiver function is horizontal response, which eliminate multiple P wave reflection and retain direct wave and P-S converted waves, and is sensitive to the vertical variation of S wave velocity. Velocity structure beneath a seismic station has different response to radial and vertical component of an accident teleseismic P wave. To avoid the limits caused by a simplified assumption on the vertical response, the receiver function method is mended. In the frequency domain, the transfer function is showed by the ratio of radical response and vertical response of the media to P wave. In the time domain, the radial synthetic waveform can be obtained by the convolution of the transfer function with the vertical wave. In order to overcome the numerical instability, generalized reflection and transmission coefficient matrix method is applied to calculate the synthetic waveform so that all multi-reflection and phase conversion response can be included. A new inversion method, VFSA-LM method, is used in this study, which successfully combines very fast simulated annealing method (VFSA) with damped least square inversion method (LM). Synthetic waveform inversion test confirms its effectiveness and efficiency. Broadband teleseismic P waveform inversion is applied in lithospheric velocity study of China and its vicinage. According to the data of high quality CDSN and IRIS, we obtained an outline map showing the distribution of Asian continental crustal thickness. Based on these results gained, the features of distribution of the crustal thickness and outline of crustal structure under the Asian continent have been analyzed and studied. Finally, this paper advances the principal characteristics of the Asian continental crust. There exist four vast areas of relatively minor variations in the crustal thickness, namely, northern, eastern southern and central areas of Asian crust. As a byproduct, the earthquake location is discussed, Which is a basic issue in seismology. Because of the strong trade-off between the assumed initial time and focal depth and the nonlinear of the inversion problems, this issue is not settled at all. Aimed at the problem, a new earthquake location method named SAMS method is presented, In which, the objective function is the absolute value of the remnants of travel times together with the arrival times and use the Fast Simulated Annealing method is used to inverse. Applied in the Chi-Chi event relocation of Taiwan occurred on Sep 21, 2000, the results show that the SAMS method not only can reduce the effects of the trade-off between the initial time and focal depth, but can get better stability and resolving power. At the end of the paper, the inverse Q filtering method for compensating attenuation and frequency dispersion used in the seismic section of depth domain is discussed. According to the forward and inverse results of synthesized seismic records, our Q filtrating operator of the depth domain is consistent with the seismic laws in the absorbing media, which not only considers the effect of the media absorbing of the waves, but also fits the deformation laws, namely the frequency dispersion of the body wave. Two post stacked profiles about 60KM, a neritic area of China processed, the result shows that after the forward Q filtering of the depth domain, the wide of the wavelet of the middle and deep layers is compressed, the resolution and signal noise ratio are enhanced, and the primary sharp and energy distribution of the profile are retained.

井间地震技术研究和应用

Cross well seismic technique is a new type of geophysical method, which observes the seismic wave of the geologic body by placing both the source and receiver in the wells. By applying this method, it averted the absorption to high-frequency component of seismic signal caused by low weathering layers, thus, an extremely high-resolution seismic signal can be acquired. And extremely fine image of cross well formations, structure, and reservoir can be achieved as well. An integrated research is conducted to the high-frequency S-wave and P-wave data and some other data to determine the small faults, small structure and resolving the issues concerning the thin bed and reservoir's connectivity, fluid distribution, steam injection and fracture. This method connects the high-resolution surface seismic, logging and reservoir engineering. In this paper, based on the E & P situation in the oilfield and the theory of geophysical exploration, a research is conducted on cross well seismic technology in general and its important issues in cross well seismic technology in particular. A technological series of integrated field acquisition, data processing and interpretation and its integrated application research were developed and this new method can be applied to oilfield development and optimizing oilfield development scheme. The contents and results in this paper are as listed follows: An overview was given on the status quo and development of the cross well seismic method and problems concerning the cross well seismic technology and the difference in cross well seismic technology between China and international levels; And an analysis and comparison are given on foreign-made field data acquisition systems for cross-well seismic and pointed out the pros and cons of the field systems manufactured by these two foreign companies and this is highly valuable to import foreign-made cross well seismic field acquisition system for China. After analyses were conducted to the geometry design and field data for the cross well seismic method, a common wave field time-depth curve equation was derived and three types of pipe waves were discovered for the first time. Then, a research was conducted on the mechanism for its generation. Based on the wave field separation theory for cross well seismic method, we believe that different type of wave fields in different gather domain has different attributes characteristics, multiple methods (for instance, F-K filtering and median filtering) were applied in eliminating and suppressing the cross well disturbances and successfully separated the upgoing and downgoing waves and a satisfactory result has been achieved. In the area of wave field numerical simulation for cross well seismic method, a analysis was conducted on conventional ray tracing method and its shortcomings and proposed a minimum travel time ray tracing method based on Feraiat theory in this paper. This method is not only has high-speed calculation, but also with no rays enter into "dead end" or "blinded spot" after numerous iterations and it is become more adequate for complex velocity model. This is first time that the travel time interpolation has been brought into consideration, a dynamic ray tracing method with shortest possible path has been developed for the first arrivals of any complex mediums, such as transmission, diffraction and refraction, etc and eliminated the limitation for only traveling from one node to another node and increases the calculation accuracy for minimum travel time and ray tracing path and derives solution and corresponding edge conditions to the fourth-order differential sonic wave equation. The final step is to calculate cross well seismic synthetics for given source and receivers from multiple geological bodies. Thus, real cross-well seismic wave field can be recognized through scientific means and provides important foundation to guide the cross well seismic field geometry designing. A velocity tomographic inversion of the least square conjugated gradient method was developed for cross well seismic velocity tomopgraphic inversion and a modification has been made to object function of the old high frequency ray tracing method and put forward a thin bed oriented model for finite frequency velocity tomographic inversion method. As the theory model and results demonstrates that the method is simple and effective and is very important in seismic ray tomographic imaging for the complex geological body. Based on the characteristics of the cross well seismic algorithm, a processing flow for cross well seismic data processing has been built and optimized and applied to the production, a good section of velocity tomopgrphic inversion and cross well reflection imaging has been acquired. The cross well seismic data is acquired from the depth domain and how to interprets the depth domain data and retrieve the attributes is a brand new subject. After research was conducted on synthetics and trace integration from depth domain for the cross well seismic data interpretation, first of all, a research was conducted on logging constraint wave impedance of cross well seismic data and initially set up cross well seismic data interpretation flows. After it applied and interpreted to the cross well seismic data and a good geological results has been achieved in velocity tomographic inversion and reflection depth imaging and a lot of difficult problems for oilfield development has been resolved. This powerful, new method is good for oilfield development scheme optimization and increasing EOR. Based on conventional reservoir geological model building from logging data, a new method is also discussed on constraining the accuracy of reservoir geological model by applying the high resolution cross well seismic data and it has applied to Fan 124 project and a good results has been achieved which it presents a bight future for the cross well seismic technology.

多波分量地震资料处理解释方法

To deal with some key problems in multi-component seismic exploration, some methods are introduced in this thesis based on reading amounts of papers about multi-component seismic theories and methods. First, to find a solution for the detection of the fracture density and orientation in igneous, carbonate and shale reservoirs, a large amount of which exist in domestic oil fields with low exploration and development degree, a new fast and slow shear waves separation method called Ratio Method based on S-wave splitting theory is discussed in this thesis, through which the anisotropy coefficient as well as fracture parameters such as density and azimuthal angle can be acquired. Another main point in this thesis involves the application of seismic velocity ratio (Vp/Vs) to predict the Hthological parameters of subsurface medium. To deal with the unfeasibility of velocity ratio calculation method based on time ratio due to the usually low single-noise ratio of S-wave seismic data acquired on land, a new method based on detailed velocity analysis is introduced. Third, pre-stack Kirchhoff integral migration is a new method developed in recent years, through which both S and P component seismic data as well as amplitude ratio of P/S waves can be acquired. In this thesis, the research on untilizing the P and S wave sections as well as amplitude ratio sections to interpret low-amplitude structures and lithological traps is carried out. The fast and slow shear wave separation method is then be applied respectively to detect the density and azimuthal angle of fractures in an igneous rock gas reservoir and the coal formation in a coal field. Two velocity ratio-calculating methods are applied respectively in the lithological prediction at the gas and coal field after summarizing a large amount of experimental results draw domestically and abroad. P and S wave sections as well as amplitude ratio sections are used to identify low-amplitude structures and lithological traps in the slope area of a oil-bearing sedimentary basin. The calculated data concerning fracture density and azimuthal angle through the introduced method matches well with the regional stress and actual drilling data. The predicted lithological data reflects the actual drilling data. Some of the low-amplitude and lithological traps determined by Kirchhoff migration method are verified by the actual drilling data. These results indicate that these methods are very meaningful when dealing with complex oil and gas reservoir, and can be applied in other areas.

地幔柱与岩浆成矿-峨嵋山大火成岩省地震波速层物相解释及典型PGE矿床成岩成矿机制研究

地幔柱概念在19世纪60至70年代就被提出，但是由于板块构造理论在解释地球上岩浆活动的分布规律时取得了空前的成功，在当时这一理论是被排斥的。板块边界概念可以解释地球上绝大部分的岩浆产出，但在解释板内岩浆的成因时往往显得力不从心，尽管这些岩浆的体积只占地球岩浆总量的2％。地幔柱理论模型发展到现在得到不同学科的支持。地质学、地球化学、地球物理学、古生物学、比较行星学、实验岩石学等等都提供了直接或间接的证据，证明地幔柱几乎存在整个地：质历史时期。当前地幔柱理论中在地球化学领域有两大研究热点：高钦低钦玄武岩的起源以及地幔柱中是否存在循环俯冲洋壳物质。完全解决这些问题才可能深入系统地建立地慢柱成矿作用模型。现在已经建立了一些矿床类型与地慢柱作用的联系：如现在认为赋存在金伯利岩中的金刚石矿床的形成与地慢柱作用密不可分，一些岩浆硫化物矿床和岩浆氧化物矿床很显然是地慢柱岩浆作用形成的，如西伯利亚火成岩省的Noril'sk-Talnakh铜镍铂族元素矿床以及KeweenawaJI大陆裂谷体系的Dultlth杂岩体的Cu-Ni矿床。另外还有赋存在大型基性一超基性层状岩体中的PGE、Ni和cu矿床，如Great Dyke和布什维尔德杂岩体。一些超大型热液矿床也与地慢柱有可能的联系（Pirajno，2000）：如270oMa形成的超大型Kidd Creek火山成因块状硫化物矿床（Bleeker et al.，1 999；Wynan et al.，1999）和南澳大利亚1600Ma形成的超大型olymPicD翻矿床。本文的研究工作包含两方面内容：通过热力学计算峨眉山玄武岩在深部的结晶分异，对峨眉山大火成岩省的岩浆量分布和岩浆氧化物矿床（华Ti磁铁矿矿床）的分布以及下地壳高波速层的物相进行理论解释；对峨眉山大火成岩省金宝山PGE典型矿床进行成岩成矿的地球化学研究，预测整个大火成岩省的岩浆硫化物矿床产出位置。大多数峨眉山玄武岩的 MgO＜7％，Ni为4-232ppm，它们是原始岩浆结晶分异后的产物。峨眉山玄武岩省下地壳和上地幔之间存在厚度为：8-25km1，P彼速为7.1-7.8km／s的附加层（高地震波速层）。滇西地区出露的洲套第三纪富碱斑岩，地球化学和同位素研究表明斑岩的岩浆源是来自“壳一慢混合层”，源区的形成时代为220-25Ma，与峨眉山玄武岩的形成时代一致。所以有理由认为该附加层是由峨眉山玄武岩在此结晶分异形成的。与地慢柱有关的洋岛Hawaii、Marquesas Islands；海底高原Oniong Java、大陆火山岩省ColumbiaRiver Plateaus地震彼研究都表明在上地慢顶部有一高速附加层，Farnetani etal.（1996）的研歼表明高速附加层是由来自地幔柱的岩浆在此结晶分异形成的。玄武岩是一种混合的部分熔融产物，是不同成分的地幔橄榄岩在不同的压力下熔出的。这种降压熔融高温高压实验是做不到的。熔出的熔体成分是温度、压力及橄榄岩成分（源区）的函数，形成的岩浆是一个多压熔融的集合体。热力学计算能够较为精确地计算出生成的岩浆成分和约束岩浆产生的过程。岩浆的结晶分异也是同样的情形，尤其是分离结晶过程，实验岩石学是很精确难模拟其过程的。热力学计算使用的MELTS程序，MELTS适用范围很广，适用于模拟岩石熔融生成岩浆和岩浆的冷却结晶。现今峨眉山大火成岩省的地壳厚度为40恤，这被认为是后期褶皱加厚的缘故。根据峨眉山玄武岩中辉石斑晶成分和玄武岩本身成分计算出分异结晶的压力为6kb，那么当时的地壳厚度约为20km：选择氧逸度为QFM，这一氧逸度范围认为是大多数大陆溢流玄武岩结晶分异时的氧化还原环境。热力学计算结果通过峨眉山玄武岩成分进行约束和验证。Al2O3、NaZO＋K 20、CaO与MgO计算的演化趋势线与实际观察的演化符合较好，橄榄石和斜方辉石的结晶使得CaO随着MgO的降低而增高；当单斜辉石成为液相线矿物时，cao也随着Mgo的降低而降低了。单斜辉石在岩浆演化到MgO＝10．3％时成为液相线矿物。整个计算过程中斜长石未成为液相线矿物，这与大多数玄武岩不具有Eu异常是一致的，并月＿Al2O3随着MgO的减小单调增加也说明了这点。不过大多数峨眉山玄武岩常含有斜长石斑晶，这是低压下结晶分异的结果。由于斜长石密度小，所有很难与高铁玄武岩分离。整个计算的难点也是创新点是波速计算。通过分离的堆晶矿物组合中各种矿物的成分和质量分数计算的附加层波速比观察值高，不过堆积岩体常常会有残留岩浆存在矿物晶粒间，这样会降低岩石的压缩波速。大型基性一超基性岩体常常会残留有或者捕获5-30％的岩浆。假定两个高波速附加层分别捕获7叭，和巧％的残留岩浆，计算的结果就大体等于观察值。热力学和质量平衡计算研究表明：高地震波速层为橄榄辉石岩一辉石岩的巨型侵入岩体；峨眉山中岩区的岩浆量最大也符合含V-Ti磁铁矿矿床只产在中岩区，如太和、白马、攀枝花、红格等岩体；西岩区的岩浆量最小表明几乎没有可能在西岩区形成有规模的V-Ti磁铁矿矿床，实际观察仅仅只见到数量少而小的岩体；东岩区下地壳厚达20灿1的高波速层暗示东岩区上地壳的侵入岩体积也应该具有相当规模，应该是V－Ti磁铁矿矿床成矿区。目前在东岩区很少发现与峨眉山玄武岩有关的岩浆矿床的主要原因是：东岩区的剥蚀深度不够，没有可观的侵入岩体出露，而中岩区侵入岩都侵入在元古代地层中。按照质量平衡的计算方法，最保守的估算整个峨眉地慢柱岩浆事件产生的岩浆量为8.9*106km3，上地壳峨眉山玄武岩和侵入岩体积为3.9*106km3。如果按照初始覆盖面积5x106km2计算（与西伯利亚暗色岩初始覆盖面积相当），喷发高峰期为2Ma，计算的喷发速率为3.9km3／year。这并不亚于西伯利亚暗色岩的喷发速率4km3／year。这对于研究峨眉山大火成岩浆事件与二叠·三叠交界或end-QuadaluPian生物灭绝之间的可能联系具有重要意义。本文另一方面的研究工作是：首先系统地介绍了岩浆硫化物矿床的基本原理，然后通过金宝山PGE矿床实例研究，提出金宝山岩体成岩模式，并且对整个峨眉山大火成岩省的岩浆硫化物矿床产出位置进行理论预测。详细地球化学研究表明金宝山镁铁一超镁铁岩是峨眉山大火成岩省古老火山岩浆房的残留物。岩体主要由底部超镁铁岩和上部镁铁岩组成，两种岩石的质量大致相同。根据超镁铁岩的矿物组合计算的成岩时的氧逸度较高，热力学方法计算的成岩压力为2kb左右。超镁铁岩的包嵌结构和铁铁岩的微晶一细晶结构说明超镁铁岩为镁铁岩结晶的矿物堆积形成的。镁铁一超镁铁岩的蚀变程度不同以及Sc、Sr、Eu等元素在两类岩石中的不同特征指示了整个成岩过程。金宝山岩体的原始岩浆 MgO=8％说明高镁玄武岩并不是形成PGE矿床的必要条件。金宝山的成岩模式是：在火山喷发前，岩浆侵位时橄榄石和少量铬尖晶石先结晶，沉淀在岩浆房底部；随后结晶的是斜方辉石和斜长石，斜方辉石也沉淀在岩浆房底部，斜长石由于密度较小集中中岩浆房上部，岩浆房的中部是：少量的斜长石小斑晶。由于斜方辉石和斜长石的结晶，这样岩浆中的Sc、Sr和Eu就会亏损，也是岩浆房底部堆积岩的原始捕获岩浆。火山喷发后，由于压力的突然降低，岩浆房底部的堆晶会发生再熔融，几乎消耗掉所有的斜方辉石，橄榄石也呈熔蚀状浑圆形态，重新熔融的斜方辉石导致超镁铁岩中残留岩浆比原始捕获岩浆更加富Sc，这种岩浆由于富MgO和在快速冷却的环境下同时结晶，最终形成光性方位一致的单刹辉石。喷发后岩浆房空间的剩余导致围岩-灰岩进入，造成岩浆房中剩余岩浆强烈的碳酸盐化。峨眉山玄武岩Cr-Mg#的相关关系定义一条正常玄武岩演化线。大多数这些玄武岩的Ni也保持了这种演化关系，其中低钦玄武岩和过渡型高钦玄武岩Ni-Mg＃相关关系远离了正常演化线，这些玄武岩的Cu-Mg#相关关系也有类似的情形。峨眉山低钦和过渡类型高钦玄武岩Ni和 Cu的非正常亏损，表明它们在地表下经历了硫饱和事件。金宝山岩浆硫化物矿床成岩模型的建立，为在整个大火成岩省寻找岩浆硫化物矿床提供了一种新认识。低钦和过渡型高钦玄武岩的古老火山口下部是岩浆硫化物矿床的所在地。

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.

Multiphase technique improves mud-pulse velocity calculations

Measurement while drilling (MWD) has become a popular survey technology to monitor directional data, drilling data, formation evaluation data and safety data in the world. And closed loop drilling shows promise in recent years. Obviously, the method of tr