Lg尾波衰减成像研究——北京及邻区衰减成像

With the development of the technology of earthquake observation, more and more researchers work at many fields' of seismicity using seismic kinetic property, as the result, the study of attenuation has also made great progress, especially in the mechanism of the attenuation and the physical process. Aki put forward single back scattering theory to explain the forming of the seismic coda wave in 1969. Then, researchers started to develop the study in seismic scattering and attenuation. My thesis is also based on that theory. We assume that the Lg wave is a superposition fo high-mode surface waves, the coda of Lg is caused by scattering. Sato proposed Single Isotropic Scattering model (SIS model) to interpreted the scatter property, and he also formulated the geometrical spreading term. Then Xie (1988) developed the single spectral-ratio (SSR) method to obtain the Lg coda Q and the frequency dependent factor n. Later, he get to lateral images in the area of scatter ellipse. SSR method is explored and used in the study of Lg coda waves of regional earthquakes in my thesis. Choosing the earthquakes records with high ratio of signal-noise ,which were recorded at the stations from 1989 to 1999, we obtain the single trace Lg coda Q and its frequency dependent factor n. The results proved that SIS model is the reasonable model to explain the Lg coda wave, and SSR method also can be used to process Lg coda of regional earthquakes to get to the satisfied Lg coda Q. Based on the Lg coda Q we obtained using the former method, we explore the programs to inverse the regional Lg coda Q independently, and then make use of them to inverse the Lg coda Q of Beijing and adjacent area. The inversion result is satisfied. We conclude that the distribution of Qo (Q in lHz) is marked by the inhomogeneity, which is related to the tectonic structure: The value of Qo in uplift area, for example, Yanshan uplift, Taihang uplift, Luxi Uplift, is higher than the depression area, for example, Jizhong depression, Huanghua depression, and Jiyang depression, and the border between the higher Q area and lower Q area is very clear; Lg coda Q is also related to the velocity structure, higher velocity area is also with higher Q, lower velocity area is with lower Q; and higher heat-flow area is companied with lower Q. All in all, the value of Q reflects the difference of characteristics of lithofaces, porosity, the liquid content between the pores and heat flow. So, the Q value difference between uplift area and depression area reveals the difference of tectonic structure, lithology and physical character of the rock. So, the study of Lg coda Q is help to understand the earthquakes propagation mechanism through the inhomogenous medium, the cause of the coda, attenuation mechanism of the coda. Making use of the lateral images of Q, with velocity images, heat flow results, and other experimental result, we will be promoted to understand the complex structure of the crust, its inhomogenous character, and so on.

丘陵油田三间房组储层三维地质模型和砂体预测

This paper is belonging to Chinese Petrochemical Corporation's science and technology project. Although it is difficult, it has important theoretical and practical value. The study was aimed to reveal inhomogeneity of two kinds of reservoirs of fan-shaped delta and braided river by using new theories, new methods and new technology about 3-D model building and reservoir knowledge repository throughout the world, and to build reservoir knowledge repository and 3-D geological model which would predict the type of sand body forming reason and distribution rule in order to improve exploration result in Qiuling oil fields. Multi-discipline theories such as petroleum structure geology, reservoir geology, petroleum geology, sequence geology, logging geology, geomathematics and so on are used as guide. The information of geology, seism, logging and production test is combined. Outcrop area and overlap area are combined. By making full use of computer, stable structure, reservoir geometric shape, spatial distribution and inhomogeneity of bed of interest are investigated, described and characterized. Petroleum pool 3-D static geological model of reservoir knowledge repository was built. Sand body distribution was predicted. It has guided oil development, lowed the investment and improved development benefits. Several results are achieved as follows: (1) Strata framework of Sanjianfang group in Qiuling oil field has been established. (2) Geometric shape, spatial distribution and evolve rule of two different forming reason's reservoir of fan-shaped delta and braided river of Sanjianfang group in Qiuling oil field are discussed. (3) The two kinds of reservoirs have lower pore and permeability and very strong inhomogeneity. (4) Reservoir knowledge repository of two different forming reasons has been built of Sanjianfang group, which includes 5 geological knowledge sublibrary. (5) 3-D geological model of two kinds of forming reason's reservoirs has been built. (6) That same sequence instruction a simulation and probability field were used to predict sand body of Sanjianfang group was put forward. Coincidence rate is high after production test. It shows this method has great popularity value. (7) A set of theories, methods and technologies of knowledge repository of two kinds of reservoir of braided river and fan-shaped delta and 3-D geological model building were finished. (8) A set of theories, methods and technologies of investigating, describing, characterizing and predicting two kinds of oil pool were developed. It gets noticeable economic benefit after exploration. Theory and method about extrusion basin are developed.

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

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.

鄂尔多斯盆地苏里格气田储层次生孔隙性质与分布预测

Permian reservoir in Sulige area of Ordos Basin, on which this paper focused, belongs to fluvial-delta lithofacies. The majority formations in this area are complicated channel sand deposit with serious inhomogeneity which makes natural gas exploration be very tough in this area. This inhomogeneity can be found everywhere both in large horizontal area and vertical profile of inner and interbedded formations.This paper studied the inhomogeneity characteristic of Permian formation in sulige area of Ordos Basin according to the logging data.Correlating with core data, a criterion to distinguish different type of reservoirs by using logging data is determined after the study of logging response is done considering the diverse conditions of deposit environments, lithology and reservoir space. The characteristic relationships between the various type formations and logging responses fully and systemically are established.It investigated reservoir parameter calculation methods amply. Combining the conventional and special logging data, basing on the feature of low porosity -permeability formation of sulige area, a set of methods to calculate reservoir parameters was formed including primary porosity, secondary porosity, fracture porosity, permeability and water saturation under the conditions of both low porosity-permeability and inhomogenous reservoirs. One thing should be pay close attention is the parameter M for calculating saturation. It is found that the M in low porosity -permeability formation decreases as the porosity decrease, which is opposite to the law that M increases as the porosity decrease in the formation with intermediate to high porosity and permeability. This view has innovated the traditional theory and offered theory basis for the logging interpretation of low porosity - permeability reservoir. Meanwhile it also improved the Arqi formula theoretically and enhanced the logging interpretation accuracy and rescued a number of formations which has been thought to be hopeless according to the old theory.By using advantage logging interpretation procedure, a territorial synthetic geology evaluation to the inhomogeneous reservoir was completed basing on the single well interpretation. All the reservoir evaluation parameters including sand formation thickness, primary porosity, secondary porosity were calculated and evaluated. The rules of changing and development for sand formation thickness, sand physical properties and secondary porous were found at different formations of upper part of the Member 8 of Shihezi, lower part of the Member 8 of Shihezi, the Member 1 of shanxi and the Member 2 of shanxi individually. Evaluation and Correlation of these five formations were also completed and one conclusion was arrived: upper part of the Member 8 of Shihezi formation has the best performance followed by the lower part of the Member 8 of Shihezi, the Member 1 of shanxi and the Member 2 of shanxi formation.After studied the relationship between reservoir deposition characteristic and the natural gas richness, it is regarded that reservoir inhomogeneity is the key issue of the impaction on the natural gas. Natural gas in Sulige gas field was mainly accumulated in sands of channel bar, distributary channel and debouchure bar. Especially, the quartz sand with rich of secondary porous space has obvious better physical properties than other reservoir and usually can forms the concentration of natural gas.