粘弹性介质中的地震波传播与波形反演研究

The real earth is far away from an ideal elastic ball. The movement of structures or fluid and scattering of thin-layer would inevitably affect seismic wave propagation, which is demonstrated mainly as energy nongeometrical attenuation. Today, most of theoretical researches and applications take the assumption that all media studied are fully elastic. Ignoring the viscoelastic property would, in some circumstances, lead to amplitude and phase distortion, which will indirectly affect extraction of traveltime and waveform we use in imaging and inversion. In order to investigate the response of seismic wave propagation and improve the imaging and inversion quality in complex media, we need not only consider into attenuation of the real media but also implement it by means of efficient numerical methods and imaging techniques. As for numerical modeling, most widely used methods, such as finite difference, finite element and pseudospectral algorithms, have difficulty in dealing with problem of simultaneously improving accuracy and efficiency in computation. To partially overcome this difficulty, this paper devises a matrix differentiator method and an optimal convolutional differentiator method based on staggered-grid Fourier pseudospectral differentiation, and a staggered-grid optimal Shannon singular kernel convolutional differentiator by function distribution theory, which then are used to study seismic wave propagation in viscoelastic media. Results through comparisons and accuracy analysis demonstrate that optimal convolutional differentiator methods can solve well the incompatibility between accuracy and efficiency, and are almost twice more accurate than the same-length finite difference. They can efficiently reduce dispersion and provide high-precision waveform data. On the basis of frequency-domain wavefield modeling, we discuss how to directly solve linear equations and point out that when compared to the time-domain methods, frequency-domain methods would be more convenient to handle the multi-source problem and be much easier to incorporate medium attenuation. We also prove the equivalence of the time- and frequency-domain methods by using numerical tests when assumptions with non-relaxation modulus and quality factor are made, and analyze the reason that causes waveform difference. In frequency-domain waveform inversion, experiments have been conducted with transmission, crosshole and reflection data. By using the relation between media scales and characteristic frequencies, we analyze the capacity of the frequency-domain sequential inversion method in anti-noising and dealing with non-uniqueness of nonlinear optimization. In crosshole experiments, we find the main sources of inversion error and figure out how incorrect quality factor would affect inverted results. When dealing with surface reflection data, several frequencies have been chosen with optimal frequency selection strategy, with which we use to carry out sequential and simultaneous inversions to verify how important low frequency data are to the inverted results and the functionality of simultaneous inversion in anti-noising. Finally, I come with some conclusions about the whole work I have done in this dissertation and discuss detailly the existing and would-be problems in it. I also point out the possible directions and theories we should go and deepen, which, to some extent, would provide a helpful reference to researchers who are interested in seismic wave propagation and imaging in complex media.

湘西地区上震旦统上部－下寒武统底部黑色岩系沉积学及地球化学

Widespread black chert-shales occur in the Ediacaran-Cambrian(E-C) boundary successions along the flank of Yangtze Platform, South China, remarkable changes in sedimentology, geochemistry and biology were recorded. Although extensive studies were carried out upon this boundary succession, the origin of black chert-shales still remain controversial. This paper focuses on the E-C black chert-shales in western Hunan, South China, upon which detailed depositional and geochemical changes are documented, accordingly a depositional model for black chert-shales is proposed. Stratigraphic anatomy across the depositional strike demonstrates that the shallow-water Dengying dolostone along the platform margin sharply pass basinward into the Liuchapo chert successions, which indicate syndepositional extensional faulting at depth could have occurred along the platform margin. The deep-water Niutitang phosphorite-rich black shales are either underlain by the Dengying dolostones on the platform margin toward platform interior or directly by the Liuchaopo chert successions farther basinwards. By detailed investigation, silica chimneys are firsly identified approximately in the chert along platform margin; two types of silica chimneys, including mounded and splayed/funnelized chert(generally brecciated) bodies are further sorted out. The mounded chert are exitbited by domed or hummocky surfaces on the top and irregular spongy to digitiform internal fabrics; within the silica mounds, abundant original vesicles/voids and/or channels were mostly plugged by initial chalcedony, quartze crystals with minor dolomite and bladed barite crystals. Splayed/funnelized brecciated chert “intrusion” cross-cut the uppermost dolostones capping to the horizon underneath, and are directly overlain by the Niutitang phosphorite-rich black shales. Their similarities to the silica chimneys reported from the oceanic spreading centres suggest a similar origin responsible for these unique silica bodies which is also supported by the microthermonmetric data and element geochemistry. High P, Ba, Fe contents and positive correlation between Fe and TOC concentrations in the Niutitang black shales indicate a high palaeo-productivity in the Early Cambrian ocean. The low Th/U and the high V/Cr, V/Sc, V/(V+Ni) ratios in the black shales suggest an anoxic water condition during this interval. Furthermore, Positive Eu anomalies and high Ba contents in the sediments also imply a hydrothermal influence on the formation of Niutitang black shales. To better constrain the placement of deep-water successions straddling the E-C boundary and the timing of hydrothermal silica chimneys, sensitive high-resoluton ion microprobe(SHRIMP) U-Pb dating of zircon grains from tuffs within the chert succession of Liuchapo Formation at Ganziping was conducted and yields a weighted-mean 206Pb/238Pb age at 536.6±5.5Ma, younger than E-C boundary age(542.0±0.3Ma). This age combined with carbon isotopic data is then proposed to correspond to the U-Pb age of zircons(538.2±1.5Ma) from the Zhongyicun member of Meishucun Formation at Meishucun in eastern Yunna, thus, the E-C boundary in Gazngziping was placed between the Dengying formations and Liuchapo formatioms. therefore, the silica chimneys took place at the beginning of the Cambrian period. The temporal coincidence of silica chimneys and negative excursions of δ13C and δ34Spy pairs suggest hydrothermal activities were likely responsible for the isotopic changes. Under such a circumstance, vast amounts of greenhouse gases(CO2, CH4, H2S), with highly 13C-depleted carbon and 34S-depleted sulfur would be released into the ocean and atmosphere. A positive shift in δ34Scas and Δ34S values from the late Ediacaran to the Early Cambrian could be a reflection of enhanced bacterial sulfate reduction(BSR), strengthened by the intensified oceanic anoxia stimulated by hydrothermal activities. Based on the analyses of sedimentology and geochemistry, a model- “oceanic anoxia induced by hydrothermal–volcanic activies” was proposed to responsible for the formation of black chert-shales during this E-C transition. Under this case, hydrothermal-volcanic activies could release large large amount of greenhouse into atmosphere and metal micronutrients into the ocean, which may lead to global warming, stratified ocean, thereby a high palaeoproductivity; on the other hand, the massive releasing of reduced hydrothermal fluids with abundant H2S, could have in turn enhanced the ocean anoxia. All of these were favourable the for preservation of organic matter, and subsequent extensive deposition of black silica-shales.

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

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.

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

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.

小波域地震资料的处理方法（WSSP）研究

The modeling formula based on seismic wavelet can well simulate zero - phase wavelet and hybrid-phase wavelet, and approximate maximal - phase and minimal - phase wavelet in a certain sense. The modeling wavelet can be used as wavelet function after suitable modification item added to meet some conditions. On the basis of the modified Morlet wavelet, the derivative wavelet function has been derived. As a basic wavelet, it can be sued for high resolution frequency - division processing and instantaneous feature extraction, in acoordance with the signal expanding characters in time and scale domains by each wavelet structured. Finally, an application example proves the effectiveness and reasonability of the method. Based on the analysis of SVD (Singular Value Decomposition) filter, by taking wavelet as basic wavelet and combining SVD filter and wavelet transform, a new de - noising method, which is Based on multi - dimension and multi-space de - noising method, is proposed. The implementation of this method is discussed the detail. Theoretical analysis and modeling show that the method has strong capacity of de - noising and keeping attributes of effective wave. It is a good tool for de - noising when the S/N ratio is poor. To give prominence to high frequency information of reflection event of important layer and to take account of other frequency information under processing seismic data, it is difficult for deconvolution filter to realize this goal. A filter from Fourier Transform has some problems for realizing the goal. In this paper, a new method is put forward, that is a method of processing seismic data in frequency division from wavelet transform and reconstruction. In ordinary seismic processing methods for resolution improvement, deconvolution operator has poor part characteristics, thus influencing the operator frequency. In wavelet transform, wavelet function has very good part characteristics. Frequency - division data processing in wavelet transform also brings quite good high resolution data, but it needs more time than deconvolution method does. On the basis of frequency - division processing method in wavelet domain, a new technique is put forward, which involves 1) designing filter operators equivalent to deconvolution operator in time and frequency domains in wavelet transform, 2) obtaining derivative wavelet function that is suitable to high - resolution seismic data processing, and 3) processing high resolution seismic data by deconvolution method in time domain. In the method of producing some instantaneous characteristic signals by using Hilbert transform, Hilbert transform is very sensitive to high - frequency random noise. As a result, even though there exist weak high - frequency noises in seismic signals, the obtained instantaneous characteristics of seismic signals may be still submerged by the noises. One method for having instantaneous characteristics of seismic signals in wavelet domain is put forward, which obtains directly the instantaneous characteristics of seismic signals by taking the characteristics of both the real part (real signals, namely seismic signals) and the imaginary part (the Hilbert transfom of real signals) of wavelet transform. The method has the functions of frequency division and noise removal. What is more, the weak wave whose frequency is lower than that of high - frequency random noise is retained in the obtained instantaneous characteristics of seismic signals, and the weak wave may be seen in instantaneous characteristic sections (such as instantaneous frequency, instantaneous phase and instantaneous amplitude). Impedance inversion is one of tools in the description of oil reservoir. one of methods in impedance inversion is Generalized Linear Inversion. This method has higher precision of inversion. But, this method is sensitive to noise of seismic data, so that error results are got. The description of oil reservoir in researching important geological layer, in order to give prominence to geological characteristics of the important layer, not only high frequency impedance to research thin sand layer, but other frequency impedance are needed. It is difficult for some impedance inversion method to realize the goal. Wavelet transform is very good in denoising and processing in frequency division. Therefore, in the paper, a method of impedance inversion is put forward based on wavelet transform, that is impedance inversion in frequency division from wavelet transform and reconstruction. in this paper, based on wavelet transform, methods of time - frequency analysis is given. Fanally, methods above are in application on real oil field - Sansan oil field.

东海陆架盆地西南部烟囱构造与油气运聚关系研究

In recent years, chimney structure has been proved one of important indicators and a useful guide to major petroleum fields exploration through their exploration history both at home and abroad. Chimney structure, which has been called "gas chimney" or "seismic chimney", is the special fluid-filled fracture swarm, which results from the boiling of active thermal fluid caused by abruptly decreasing of high pressure and high temperature in sedimentary layers of upper lithosphere. Chimney structure is well developed in continental shelf basin of East China Sea, which indicates the great perspectives of petroleum resources there. However, the chimney structure also complicated the petroleum accumulation. So the study of chimney structure on its formation, its effect on occurrence and distribution of petroleum fields is very important not only on theoretical, but also on its applied research. It is for the first time to make a clear definition of chimney structure in this paper, and the existence and practical meaning of chimney structure are illustrated. Firstly, on the viewpoint of exploration, this will amplify exploration area or field, not only in marine, but also on continent. Secondly, this is very important to step-by-step exploration and development of petroleum fields with overpressure. Thirdly, this will provide reference for the study on complex petroleum system with multi-sources, commingled sources and accumulation, multi-stage accumulations, and multi-suits petroleum system in the overlay basin. Fourthly, when the thermal fluid enters the oceanic shallow layer, it can help form gas hydrate under favorable low-temperature and high-pressure conditions. Meanwhile, the thermal fluid with its particular component and thermal content will affect the physical, chemical and ecological environments, which will help solving the problem of global resources and environment. Beginning from the regional tectonic evolution characteristics, this paper discussed the tectonic evolution history of the Taibei depression, then made an dynamical analysis of the tectonic-sedimentary evolution during the Mesozoic and Cenozoic for the East China Sea basin. A numerical model of the tectonic-thermal evolution of the basin via the Basin-Mod technique was carried out and the subsidence-buried history and thermal history of the Taibei depression were inverse calculated: it had undergone a early rapid rift and sag, then three times of uplift and erosion, and finally depressed and been buried. The Taibei depression contains a huge thick clastic sedimentary rock of marine facies, transitional facies and continental facies on the complex basement of ante-Jurassic. It is a part of the back-arc rifting basins occurred during the Mesozoic and Cenozoic. The author analyzed the diagenesis and thermal fluid evolution of this area via the observation of cathodoluminescence, scanning electron microscope and thin section, taking advantage of the evidences of magma activities, paleo-geothermics and structural movement, the author concluded that there were at least three tectonic-thermal events and three epochs of thermal-fluid activities; and the three epochs of thermal-fluid activities were directly relative to the first two tectonic-thermal events and were controlled by the generation and expulsion of hydrocarbon in the source rock simultaneously. Based on these, this paper established the corresponding model between the tectonic-thermal events and the thermal-fluid evolution of the Taibei Depression, which becomes the base for the study on the chimney structures. According to the analyses of the gas-isotope, LAM spectrum component of fluid inclusion, geneses of CO_2 components and geneses of hydrocarbon gases, the author preliminarily verified four sources of the thermal fluid in the Taibei Depression: ① dehydration of mud shale compaction, ② expulsion of hydrocarbon in the source rock; ③ CO_2 gas hydro-thermal decomposition of carbonatite; ④magma-derived thermal fluid including the mantle magma water and volatile components (such as H_2O, CO_2, H_2S, SO_2, N_2 and He etc.). On the basis of the vitrinite reflectance (Ro), homogenization temperature of fluid inclusion, interval transit time of major well-logging, mud density of the wells, measured pressure data and the results of previous studies, this paper analyzed the characteristics of the geothermal fields and geo-pressure fields for the various parts in this area, and discussed the transversal distribution of fluid pressure. The Taibei depression on the whole underwent a temperature-loss process from hot basin to cold basin; and locally high thermal anomalies occurred on the regional background of moderate thermal structure. The seal was primarily formed during the middle and late Paleocene. The overpressured system was formed during the middle and late Eocene. The formation of overpressured system in Lishui Sag underwent such an evolutionary process as "form-weaken-strengthen-weaken". Namely, it was formed during the middle and late Eocene, then was weakened in the Oligocene, even partly broken, then strengthened after the Miocene, and finally weakened. The existence of the thermal fluid rich in volatile gas is a physical foundation for the boiling of the fluid, and sharply pressure depletion was the major cause for the boiling of the fluid, which suggests that there exists the condition for thermal fluid to boil. According to the results of the photoelastic simulation and similarity physical experiments, the geological condition and the formation mechanism of chimnestructures are summarized: well compartment is the prerequisite for chimney formation; the boiling of active thermal fluid is the original physical condition for chimney formation; The local place with low stress by tension fault is easy for chimney formation; The way that thermal fluid migrates is one of the important factors which control the types of chimney structures. Based on where the thermal fluid come from and geometrical characteristics of the chimney structures, this paper classified the genetic types of chimney structures, and concluded that there existed three types and six subtypes chimney structures: organic chimney structures generated by the hydrocarbon-bearing thermal fluid in middle-shallow layers, inorganic and commingling-genetic chimney structures generated by thermal fluid in middle-deep layers. According to the seismic profiles interpretations, well logging response analysis and mineralogical and petrological characteristics in the study area, the author summarized the comprehensive identification marks for chimney structures. Especially the horizon velocity analysis method that is established in this paper and takes advantage of interval velocity anomaly is a semi-quantitative and reliable method of chimney structure s identification. It was pointed out in this paper that the occurrence of the chimney structures in the Taibei depression made the mechanism of accumulation complicated. The author provided proof of episodic accumulation of hydrocarbon in this area: The organic component in the boiling inclusion is the trail of petroleum migration, showing the causality between the boiling of thermal fluid and the chimney structures, meanwhile showing the paroxysmal accumulation is an important petroleum accumulation model. Based on the evolutionary characteristics of various types of chimney structures, this paper discussed their relationships with the migration-accumulation of petroleum respectively. At the same time, the author summarized the accumulating-dynamical models associated with chimney structures. The author analyzed such accumulation mechanisms as the facies state, direction, power of petroleum migration, the conditions of trap, the accumulation, leakage and reservation of petroleum, and the distribution rule of petroleum. The author also provides explanation for such practical problems the existence of a lot of mantle-derived CO_2, and its heterogeneous distribution on plane. By study on and recognition for chimney structure, the existence and distribution of much mantle-derived CO_2 found in this area are explained. Caused by tectonic thermal activities, the deep magma with much CO_2-bearing thermal fluid migrate upward along deep fault and chimney structures, which makes two wells within relatively short distance different gas composition, such as in well LF-1 and well LS36-1-1. Meanwhile, the author predicted the distribution of petroleum accumulation belt in middle-shallow layer for this area, pointed out the three favorable exploration areas in future, and provided the scientific and deciding references for future study on the commingling-genetic accumulation of petroleum in middle-deep layer and the new energy-gas hydrate.

生物－热催化过渡带沉积成岩演化及其成烃作用研究

The biothermocatalytic transitional zone gas is a new type of natural gas genetic theory, and also an clean, effective and high quality energy with shallow burial depth, wide distribution and few investment. Meanwhile, this puts biothermocatalytic transitional zone gas in important position to the energy resource and it is a challenging front study project. This paper introduces the concept, the present situation of study and developmental trend about biothermocatalytic transitional zone gas in detail. Then by using heat simulating of source rocks and catalysis mechanism analysis in the laboratory and studying structural evolution, sedimentation, diagenesis and the conditions of accumulation formation and so on, this paper also discusses catalytic mechanism and evolutionary model of the biothermocatalytic transitional zone gas formation, and establishes the methods of appraisal parameter and resources prediction about the biothermocatalytic transitional zone gas. At last, it shows that geochemical characteristics and differentiated mark of the biothermocatalytic transitional zone gas, and perfect natural gas genetic theory, and points out the conditions of accumulation formation, distribution characteristics and potential distribution region on the biothermocatalytic transitional zone gas m China. The paper mainly focuses on the formation mechanism and the resources potential about the biothermocatalytic transitional zone gas. Based on filed work, it is attached importance to a combination of macroscopic and microcosmic analysis, and the firsthand data are obtained to build up framework and model of the study by applying geologic theory. Based on sedimentary structure, it is expounded that structural actions have an effect on filling space and developmental cource of sediments and evolution of source rocks. Carried out sedimentary environment, sequence stratigraphy, sedimentary system and diagenesis and so on, it is concluded that diagenesis influences developmental evolution of source rocks, and basic geologic conditions of the biothermocatalytic transitional zone gas. Applying experiment simulating and catalytic simulating as well as chemical analysis, catalytic mechanism of clay minerals is discussed. Combined diagenecic dynamics with isotope fractionation dynamics, it is established that basis and method of resource appraisal about the biothermocatalytic transitional zone gas. All these results effectively assess and predict oil&gas resources about the biothermocatalytic transitional zone gas-bearing typical basin in China. I read more than 170 volumes on the biothermocatalytic transitional zone gas and complete the dissertation' summary with some 2.4 ten thousand words, draw up study contents in some detail and set up feasible experimental method and technologic course. 160 pieces of samples are obtained in oilfield such as Liaohe, Shengli, Dagang and Subei and so on, some 86 natural gas samples and more than 30 crude oil samples. Core profiles about 12 wells were observed and some 300 geologic photos were taken. Six papers were published in the center academic journal at home and abroad. Collected samples were analysised more than 1000 times, at last I complete this dissertation with more than 8 ten thousand words, and with 40 figures and 4 plates. According to these studies, it is concluded the following results and understandings. 1. The study indicates structural evolution and action of sedimentary basin influence and control the formation and accumulation the biothermocatalytic transitional zone gas. Then, the structural action can not only control accommodation space of sediments and the origin, migration and accumulation of hydrocarbon matters, but also can supply the origin of energy for hygrocarbon matters foramtion. 2. Sedimentary environments of the biothermocatalytic transitional zone gas are lake, river and swamp delta- alluvial fan sedimentary systems, having a warm, hot and humid climate. Fluctuation of lake level is from low to high., frequency, and piling rate of sedimentary center is high, which reflect a stable depression and rapidly filling sedimentary course, then resulting in source rocks with organic matter. 3. The paper perfects the natural gas genetic theory which is compound and continuous. It expounds the biothermocatalytic transitional zone gas is a special gas formation stage in continuous evolutionary sequence of organic matter, whose exogenic force is temperture and catalysis of clay minerals, at the same time, having decarbxylation, deamination and so on. 4. The methodology is established which is a combination of SEM, TEM and Engery spectrum analysis to identify microstructure of crystal morphology about clay minerals. Using differential thermal-chromatographic analysis, it can understand that hydrocarbon formation potential of different typies kerogens and catalytic method of all kinds of mineral matrix, and improve the surface acidity technology of clay minerals measured by the pyridine analytic method. 5. The experiments confirm catalysis of clay minerals to organic matter hygrocarbon formation. At low temperature (<300 ℃), there is mainly catalysis of montmorillonite, which can improve 2-3 times about produced gas of organic matters and the pyrolyzed temperature decreased 50 ℃; while at the high temperature, there is mainly catalysis of illite which can improve more than 2 times about produced gas of organic matters. 6. It is established the function relationship between organic matter (reactant) concentration and temperature, pressure, time, water and so on, that is C=f (D, t). Using Rali isotope fractionation effect to get methane isotope fractionation formula. According to the relationship between isotope fractionation of diagenesis and depth, and combined with sedimentary rate of the region, it is estimated that relict gas of the biothermocatalytic transitional zone gas in the representative basin. 7. It is revealed that hydrocarbon formation mechanism of the biothermocatalytic transitional zone gas is mainly from montmorillonite to mixed minerals during diagenesis. In interlayer, a lot of Al~(3+) substitute for Si~(4+), resulting in a imbalance between surface charge and interlayer charge of clay minerals and the occurrence of the Lewis and Bronsted acid sites, which promote to form the carbon cation. The cation can form alkene or small carbon cation. 8. It is addressed the comprehensive identification mark of the biothermo - catalytic transitional zone gas. In the temproal-spatial' distribution, its source rocks is mainly Palaeogene, secondly Cretaceous and Jurassic of Mesozoic, Triassic, having mudy rocks and coal-rich, their organic carbon being 0.2% and 0.4% respectively. The vitrinite reflection factor in source rocks Ro is 0.3-0.65%, a few up to 0.2%. The burial depth is 1000-3000m, being characterized by emerge of itself, reservoir of itself, shallow burial depth. In the transitional zone, from shallow to deep, contents of montmorillonites are progressively reduced while contents of illites increasing. Under SEM, it is observed that montmorillonites change into illite.s, firstly being mixed illite/ montmorillonite with burr-like, then itlite with silk-like. Carbon isotope of methane in the biothermocatatytic transitional zone gas , namely δ~(13)C_1-45‰- -60 ‰. 9. From the evolutionary sequence of time, distribution of the biothermocatalytic transitional zone gas is mainly oil&gas bearing basin in the Mesozoic-Neozoic Era. From the distribution region, it is mainly eastern stuctural active region and three large depressions in Bohaiwang basin. But most of them are located in evolutionary stage of the transitional zone, having the better relationship between produced, reservoir and seal layers, which is favorable about forming the biothermocatalytic transitional zone gas reservoir, and finding large gas (oil) field.

东营凹陷下第三系岩性油气藏形成条件与预测

The Dongying depression, located in the northern part of the jiyang Sag in the Buohaiwan Basin, comprises one of the major oil-producing bases of the Shengli oil-field. The prediction and exploration of subtle or litho1ogical oil traps in the oil-field has become the major confronted target. This is also one of the frontier study areas in the highly-explored oil-bearing basins in East China and abroad. Based on the integrated analysis of the geological, seismic and logging data and the theories of sequence stratigraphy, tectono-stratigraphy and petroleum system, the paper has attempted to document the characteristics of the sequence stratigraphic and structural frameworks of the low Tertiary, the syndepositional faults and their control on deposition, and then to investigate the forming conditions and distribution of the tithological oil traps in the depression. The study has set up a set of analysis methods, which can be used to effectively analysis the sequence stratigraphy of inland basins and predict the distribution of sandstone reservoirs in the basins. The major achievements of the study are as follows: 1. The low Tertiary can be divided into 4 second-order sequences and 13 third-order sequences, and the systems tracts in the third-order sequences have been also identified based on the examination and correction of well logging data and seismic profiles. At the same time, the parasequences and their stacking pattern in the deltaic systems of the third member of the Shahejie Formation have been recognized in the key study area. It has been documented that the genetic relation of different order sequences to tectonic, climatic and sediment supply changes. The study suggested that the formation of the second-order sequences was related to multiple rifting, while the activity of the syndepositional faults controlled the stacking pattern of parasequences of the axial deltaic system in the depression. 2. A number of depositional facies have been recognized in the low Tertiary on the basis of seismic facies and well logging analysis. They include alluvial fan, fan delta or braided delta, axial delta, lowstand fan, lacustrine and gravity flow deposits. The lacustrine lowstand fan deposits are firstly recognized in the depression, and their facies architecture and distribution have been investigated. The study has shown that the lowstand fan deposits are the important sandstone reservoirs as lithological oil traps in the depression. 3. The mapping of depositional systems within sequences has revealed the time and special distrbution of depositional systems developed in the basin. It is pointed out that major elastic systems comprise the northern marginal depositional systems consisting of alluvial fan, fan delta and offshore lowstand fan deposits, the southern gentle slope elastic deposits composed of shallow lacustrine, braided delta and lowstand fan deposits and the axial deltaic systems including those from eastern and western ends of the depression. 4. The genetic relationship between the syndepositional faults and the distribution of sandstones has been studied in the paper, upper on the analysis of structural framework and syndepositional fault systems in the depression. The concept of structural slope-break has been firstly introduced into the study and the role of syndepositional faults controlling the development of sequence architecture and distribution of sandstones along the hinged and faulted margins have been widely investigated. It is suggested that structural styles of the structural slope-break controlled the distribution of lowstand fan deposits and formed a favorable zone for the formation of lithological or structure-lithological oil traps in the basin. 5. The paper has made a deep investigation into the forming condition and processes of the lithological traps in the depression, based the analysis of composition of reservoir, seal and resource rocks. It is pointed out that there were two major oil pool-forming periods, namely the end of the Dongying and Guangtao periods, and the later one is the most important. 6. The study has finally predicted a number of favorable targets for exploration of lithologieal traps in the depression. Most of them have been drilled and made great succeed with new discovered thousands tons of raw oil reserves.

弱各向异性介质中传播的地震波及其速度重建

We begin our studies to make the best of information of seismic data and carry out the description of cracks parameters by extracting anisotropic information. The researching contents are: (1) velocity and polarization anomaly of seismic wave (qP and qSV wave) in weak anisotropic media; (2) reflection seismic synthetic record in anisotropic media; (3) multiple scattering induced by cracks; (4) anisotropic structure inversion and velocity reconstruction with VSP (Vertical Seismic Profile) data; (5) multi-parameters analysis of anisotropy in time-domain and depth-domain. Then we obtain results as follows: (1) We achieve approximate relation of qP and qSV wave's velocity and polarization property in weak anisotropic media. At the same time, we calculate anisotropic velocity factors and polarization anomaly of several typical sedimentary rocks. The results show there are different anisotropic velocity factors and polarization anomaly in different rocks. It is one of the primary theoretical foundation which is expected to identify lithology; (2) We calculate reflection seismic synthetic record with theoretical model; (3) We simulate scattering induced by cracks with Boundary Element Method. Numerical studies show that in the presence of cracks; spatial and scale-length distributions are important and cannot be ignored in modeling cracked solids; (4) From traveltimes information of VSP data, we study the velocity parameter inversion of seismic wave under isotropic and anisotropic models, and its result indicate that the inversion imaging under anisotropic model will not destroy the original features of isotropic model, but it will bring on some bigger error if we adopt the method of isotropic model for anisotropic model data. Further more, basing on the study we develop the CDP mapping technology of reflecting structure under isotropic and anisotropic models, and we process real data as a trial of the methods; (5) We study the problem of initial model reconstruction of anisotropic parameters structure represented by Anderson parameter in depth domain for surface data.

山东东营凹陷浊积岩油藏形成机制及预测技术

More and more exploration practises and research results prove that deep-water turbidites and slump turbidites in abrupt slope belt as well as slump turbidites in deposit center are favorable places to find subtle reservoirs. Based on sedimentology, sequence stratigraphy, tectonics and petroleum geology and with the tool of computer, this paper focused a set of questions such as the forming mechanism and m&a of turbidite reservoir in Dongying Sag, and obtained a systematic research result from available information of geology, seismic survey, well logging and testing, etc. Since the structural pattern and tectonic stress are controlling factors on distribution of turbidite reservoirs, deep-water sequence and shallow-water sequence have close relation with various turbidites. The author primarily disclose the regular pattern of the turbidites in the abrupt slope trouth, and predict that the abrupt slope is mainly connected to near-bank slump turbidite fans and deep-water turbidite fans; whereas the depocenter area usually developed with front delta slump fans. According to this, seismic-geology models corresponding to various deposit systems can be set up as a principle for predicting the turbiditea in faulted lacustrine basins. Based on comprehensive study of the turbidite reservoir forming condition and main controlling factors, a differential accumulation model and a multiple migration&accumulation model were also set up for abrupt slope turbidities. Referenced the m

东营凹陷辛镇复杂断块油藏形成机制和剩余油分布规律研究

In this paper, the complex faulted-block oil reservoir of Xinzhen area in Dongying depression is systematically studied from basic conditions forming faulted-block oil and gas reservoir integrating geology, seismic, logging and reservoir engineering information and computer; guided by petroleum geology, geomechanics, structural geology and geophysics and other theories. Based on analysis of background condition such as regional strata, structure and petroleum geology, structural research on geometry, kinemaitcs and dynamics, oil-controlling fault research on the seal features, sealing mechanism and sealing pattern, and research on enrichment rules and controlling factors of complex faulted-block oil reservoir are carried out to give out the formation mechanics of oil reservoir of Xinzhen complex faulted-block oil reservoir. As a result, the reservoir formation pattern is established. At the same time, through dissecting the characteristics and hydrocarbon enrichment law of complex faulted-block oil reservoir, and studying its distribution law of remaining oil after entering extra high water-cut period, a set of technologies are formed to predict complex faulted-block oil reservoir and its remaining oil distribution and to enhance oil recovery (EOR). Based on the time relationship between migration of hydrocarbon and trap formation, accumulating period of Xinzhen oil reservoir is determined. The formation of Xinzhen anticlinal trap was prior to the primary migration. This is favorable to formation of Xinzhen anticlinal hydrocarbon reservoir. Meanwhile, because anticline top caving isn't at the sane time as that of moving or faulted-trap forming inner anticline, oil and gas migrated many times and Xinzhen complex faulted-block oil reservoir formed from ES_3~(upper) to EG. Accumulating law and controlling factors of complex faulted-block reservoir are analyzed from many aspects such as regional structure background controlling hydrocarbon accumulating, plastic arch-open structure controlling oil-bearing series and reservoir types, sealing-opening of fault controlling hydrocarbon distribution and structure pattern controlling enriched trap types. Also, we established the structure pattern in Xinzhen a'ea: the arch-open of underlying strata cause expanding fracture. The main block groups developed here are shovel-like normal fault block group in the north area of Xinzhen and its associated graben block group. Block groups dominate the formation and distribution of reservoirs. We studied qualitatively and quantitatively the sealing characteristics, sealing history and sealing mechanism of faults, too. And, the sealing characteristics are evaluated and the distribution pattern of hydrocarbon controlled by faults is researched. Due to movement intensity of big faults, deep falling of downthrown block, high degree of repture and development of fracture, shallow layers close to the downthrown block of secondary faults are unfavorable to hydrocarbon accumulation. This is confirmed by the exploration practice in Xinzhen anticline. In terms of the downthrown blocks of sencondary contemporaneous faults lied in the south and north area of Xinzhen, hydrocarbon is poor close to fracture belt, while it is relatively abundant in tertiary companion faults. Because of long-term movement of faults that control hydrocarbon, fi'om ES3 to EG, six set of oil-bearing series formed. And their opening causes the inhomogeneity in hydrocarbon abundance among each block--in two flanks of anticline reservoirs are abundant while in the axial area, oil and gas are sporadic. There the sealing characteristics control oil-bearing area of oil/gas accumulation and the height of oil reservoir. Longitudinally, oil and gas are enriched in dip-flat areas in mid-plane of faults. It is established that there are four types of accumulating patterns in complex faulted-block oil reservoirs in Xinzhen. The first is accumulating pattern of lithologic oil reservoirs in E~S_3~(mid-lowwer), that is, self-generating-self-reserving-self-covering lithologic trap pattern. The second is drag-anticline accumulating pattern in Xinzhen. The structure traps are drag anticlines formed by the contemporaneous faults of the second basement in the north of Xinzhen, and the multiple source rocks involve Ek_2, Es_4, Es_3 and Es_1 members. The reservoirs are fluvial-delta sandstones of the upper member of Shahejie formation and Guantao formation, covered by regional thick mudstone of the upper member of Guantao formation and MingHuazhen formation. The third is the accumulating pattern of reverse listric fault, the third-degree fault of Xinzhen anticline limb and the reservoirs form reservoir screened by reverse listric faults. The forth is accumulating pattern of crossing faults which form closing or semi-closing faulted-blocks that accumulate hydrocarbon. The technologies of predicting remaining oil in complex faulted-block reservoir during the mid and late development stage is formed. Remaining oil in simple large faulted-blocks enriches in structural high, structural middle, structural low of thick bottom water reservoirs, points near bent edge-fault oftertiary faults and part the fourth ones with big falling displacement, microstructure high place of oil-sandbodies and areas where local well pattern isn't perfect. While that in small complex faulted-blocks enriches near small nose, small high point, angle of small faults, small oil-bearing faulted-blocks without well and areas with non-perfect well pattern. The technologies of enhancing recovery factor in complex faulted-block reservoir during the mid and late development stage is formed as follows: fine reservoir description, drilling adjust wells, designing directional wells, sub-dividing layer series of development, improving flooding pattern, changing water-injection direction and enhancing swept volume, cyclic waterflooding and gas-injection, etc. Here, directional wells include directional deflecting wells, lateral-drilling wells, lateral-drilling horizontal wells and horizontal wells. The results of this paper have been used in exploration and development of Shengli oilfield, and have achieved great social and economic profit, especially in predicting distribution of complex faulted-block reservoir, remaining oil distribution during middle and late stage of development, and in EOR. Applying the achievement of fault-closure research, new hydrocarbon-bearing blocks are discovered in flanks of Dongying central uplift and in complex blocks with proved reserves 15 million tons. With the study of remaining oil distribution law in complex faulted-block reservoirs, recovery factors are increased greatly in Dongxin, Xianhe and Linpan complex faulted-block reservoirs and accumulated oil production increment is 3 million tons.

青藏高原壳幔形变与层圈耦合数值模拟

Based on multi-principle (such as structures, tectonics and kinematics) exploratory data and related results of continental dynamics in the Tibetan plateau, the author reconstructed the geological-geophysical model of lithospherical structure and tectonic deformation, and the kinetics boundary conditions for the model. Then, the author used the numerical scheme of Fast Lagrangian Analysis of Continua (FLAC), to stimulate the possible process of the stress field and deformational field in the Tibetan plateau and its adjacent area, since the convergence-collision between the Indian continent and Eurasia continent about 50Ma ago. With the above-mentioned results, the author discussed the relationship between crustal movement in shallow layer and the deformational process in interior layers, and its possible dynamic constraints in deep. At the end of the paper, an integrative model has been put forward to explain the outline images of crust-mantle deformation and coupling in the Tibetan Plateau. (1) The characteristics of crust-mantle structure of the Tibetan plateau have been shown to be very complex, and vertical and horizontal difference is significant. The general characteristics of crust-mantle of the Tibetan plateau may be that it's layering in depth direction, and shows blocking from south to north and belting from east to west, mainly according to the results of about 20 seismic sections, such as wide-angle seismic profiles, CMP, seismic tomography and so on. (2) The crust had shortened about 2200km, while the shortening is different for different block from south to north in the Tibetan plateau. It is about 11.5mm/a in Himalayan block, about 9.0mm/a in Lhas-Gangdese block, about 7.0mm/a in Qiangtang block and Songpan-Ganzi-Kekexili block, about 8.0mm/a in Kunlun-Qaidam, and about ll.Omm/a in Qilian block, since the convergence-collision between the Indian continent and Eurasia continent about 50Ma ago. Which - in demonstrates the shortening rate decreases from south to north, but this rate increases near the north edge of the Tibetan plateau. The crust thickening rate is about 0.4mm/a in the whole Tibetan plateau; and this rate is about 0.5mm/a in Himalayan block, about 0.4mm/a in Lhas-Gangdese block, about 0.3mm/a in Qiangtang block, about 0.2mm/a in Songpan-Ganzi-Kekexili block and about O.lmm/a in Kunlun-Qaidam-Qilian block, since the convergence-collision between the Indian continent and Eurasia continent about 50Ma ago. This implies that the thickening rate decreases in the blocks of the Tibetan plateau. From south to north, the displacement of eastern boundary in the Tibetan plateau is about 37mm/a in Himalayan block, about 45mm/a in Lhas-Gangdese block, about 47mm/a in Qiangtang block, about 43mm/a in Songpan-Ganzi-Kekexili block, and about 35mm/a in Kunlun-Qaidam-Qilian block, since the collision-matching between the Indian continent and Eurasia continent had happened about 50Ma ago. This implies that the rate of eastward displacement is biggest in the middle of plateau, and decreases to both sides. The transition of S-N compression stress field in Tibetan Plateau, since about 28Ma+ ago, may be caused by two reasons: On one hand, the movement direction of Eurasia continent changed from northward to southward about 28Ma± ago in the northern plateau. On the other hand, the front belt that is located between India continent's and Eurasia continent's convergence-collision, had moved southward to high Himalayan from Indus-Brahmaputra suture almost at the same time in southern plateau. Affected by the stress field, the earlier tectonics rotated clockwise, NE and NW conjugate strike-slip faults developed, and the SN rift formed. This indicated that the EW movement started. The ratio between upper crust and lower crust of different blocks from south to north in the Tibetan plateau during the process of deformation are as following: about 3.5~5:1 in Himalayan block, about 1~5: 3-4 (which is about 1:3o--4 in south and about 4~5:3 in north) in Lhas-Gangdese block, about 1:3~447mm/a in these blocks: Which is located to the north of Banggong-nujiang suture.

边坡动力响应分析及应用研究

As a marginal subject, dynamic responses of slopes is not only an important problem of engineering geology (Geotechnical problem), but also of other subjects such as seismology, geophysics, seismic engineering and engineering seismic and so on. Owning to the gulf between different subjects, it is arduous to study dynamic responses of slopes and the study is far from ripeness. Studying on the dynamic responses of slopes is very important in theories as well as practices. Supported by hundreds of bibliographies, this paper systemically details the development process of this subject, introduces main means to analyze this subject, and then gives brief remarks to each means respectively. Engineering geology qualitative analysis is the base of slopes dynamic responses study. Because of complexity of geological conditions, engineering geology qualitative analysis is very important in slopes stability study, especially to rock slopes with complex engineering geology conditions. Based on research fruits of forerunners, this paper summarizes factors influencing slopes dynamic stability into five aspects as geology background, stratums, rock mass structure, and topography as well as hydrogeology condition. Based on rock mass structure controlling theory, engineering geology model of the slope is grouped into two typical classes, one is model with obvious controlling discontinuities, which includes horizontal bedded slope, bedding slope, anti-dip slope, slide as well as slope with base rock and weathered crust; the other is model without obvious controlling discontinuities, which includes homogeneous soil slope, joint rock mass slope. Study on slope failure mechanism under dynamic force, the paper concludes that there are two effects will appear in slope during strong earthquake, one is earthquake inertia force, the other is ultra pore pressure buildup. The two effects lead to failure of the slope. To different types of slope failure, the intensity of two effects acting on the slope is different too. To plastic flow failure, pore pressure buildup is dominant; to falling rock failure and toppling failure, earthquake inertia force is dominant in general. This paper briefly introduces the principle of Lagrangian element method. Through a lot of numerical simulations with FLAC3D, the paper comprehensively studies dynamic responses of slopes, and finds that: if the slope is low, displacement, velocity and acceleration are linear enlarging with elevation increasing in vertical direction; if the slope is high enough, displacement, velocity and acceleration are not linear with elevation any more, on the other hand, they fluctuate with certain rhythm. At the same time, the rhythm appears in the horizontal direction in the certain area near surface of the slope. The distribution form of isoline of displacement, velocity and acceleration in the section of the slope is remarkably affected by the slope angle. In the certain area near the slope surface, isoline of displacement, velocity and acceleration is parallel to the surface of the slope, in the mean time, the strike direction of the extreraum area is parallel to the surface of the slope too. Beyond this area, the isoline direction and the strike direction of the extremum area turn to horizontal with invariable distance. But the rhythm appearing or not has nothing to with the slope angle. The paper defines the high slope effect and the low slope effect of slopes dynamic responses, discusses the threshold height H^t of the dynamic high slope effect, and finds that AW is proportional to square root of the dynamic elastic moduli El P , at the same time, it is proportional to period Tof the dynamic input. Thus, the discriminant of H^t is achieved. The discriminant can tell us that to a slope, if its height is larger than one fifth of the wavelength, its response regular will be the dynamic high slope effect; on the other hand, its response regular will be the dynamic low slope effect. Based on these, the discriminant of different slopes taking on same response under the same dynamic input is put forward in this paper. At the same time, the paper studies distribution law of the rhythm extremum point of displacement, velocity and acceleration, and finds that there exists relationship of N = int among the slope height H, the number of the rhythm extremum VHlhro) point N and ffthre- Furthermore, the paper points out that if N^l, the response of the slope will be dynamic high slope effect; \fNdeep fractures are well developed at the left abutment slope, so stability of the slope is a key engineering geological problem. The paper studies the dynamic stability of the slope. In the analysis, considering regional geology and regional seismology of the slope, combining with characteristics of the slope as well as the scale of the project, PGA(peak ground acceleration) of the site is decided (about 197.1cm/s2) by engineering seismology, and then the seismic input used in the slope dynamic analysis is determined. Comprehensive studies are carried out on the slope, especially to its deep fractures, and then the paper concludes that the deep fractures of slope are the result of the combination unloading effect of gravity and tectonic stress. At the same time, the failure model of the slope under dynamic input is attained. Based on these, the stability is comprehensively studied for section IV-IV with numerical simulation method as well as method putting forward in chapter 5. At last, the paper concludes that under dynamic input, the section will slide along fault f9, some deep fractures and fault/5 with certain permanent displacement, and this must be taken into consideration in the engineering.

波动方程辛几何算法三维叠前深度偏移及其并行实现

Oil and scientific groups have been focusing on the 3D wave equation prestack depth migration since it can solve the complex problems of the geologic structure accurately and maintain the wave information, which is propitious to lithology imaging. The symplectic method was brought up by Feng Kang firstly in 1984 and became the hotspot of numerical computation study. It will be widely applied in many scientific field of necessity because of its great virtue in scientific sense. This paper combines the Symplectic method and the 3-D wave equation prestack depth migration to bring up an effectual numerical computation method of wave field extrapolatation technique under the scientific background mentioned above. At the base of deep analysis of computation method and the performance of PC cluster, a seismic prestack depth migration flow considering the virtue of both seismic migration method and Pc cluster has formatted. The software, named 3D Wave Equation Prestack Depth Migration of Symplectic Method, which is based on the flow, has been enrolled in the National Bureau of Copyright (No. 0013767). Dagang and Daqing Oil Field have now put it into use in the field data processing. In this paper, the one way wave equation operator is decompounded into a phase shift operator and a time shift operator and the correct item with high rank Symplectic method when approaching E exponent. After reviewing eliminating alias frequency of operator, computing the maximum angle of migration and the imaging condition, we present the test result of impulse response of the Symplectic method. Taking the imaging results of the SEG/EAGE salt and overthrust models for example and seeing about the imaging ability with complex geologic structure of our software system, the paper has discussed the effect of the selection of imaging parameters and the effectuation on the migration result of the seismic wavelet and compared the 2-D and 3-D prestack depth migration result of the salt mode. We also present the test result of impulse response with the overthrust model. The imaging result of the two international models indicates that the Symplectic method of 3-D prestack depth migration accommodates great transversal velocity variation and complex geologic structure. The huge computing cost is the key obstruction that 3-D prestack depth migration wave equation cannot be adopted by oil industry. After deep analysis of prestack depth migration flow and the character of PC cluster ,the paper put forward :i)parallel algorithms in shot and frequency domain of the common shot gather 3-D wave equation prestack migration; ii)the optimized setting scheme of breakpoint in field data processing; iii)dynamic and static load balance among the nodes of the PC cluster in the 3-D prestack depth migration. It has been proven that computation periods of the 3-D prestack depth migration imaging are greatly shortened given that adopting the computing method mentioned in the paper. In addition,considering the 3-D wave equation prestack depth migration flow in complex medium and examples of the field data processing, the paper put the emphasis on: i)seismic data relative preprocessing, ii) 2.5D prestack depth migration velocity analysis, iii)3D prestack depth migration. The result of field data processing shows satisfied application ability of the flow put forward in the paper.

中国大陆及邻区地震活动背景、环境触发与地震预测研究

This dissertation that includes most of the P. PH.D research work during 2001~2002 covers the large-scale distribution of continental earthquakes in mainland China, the mechanism and statistic features of grouped strong earthquakes related to the tidal triggering, some results in earthquake prediction with correlativity analysis methods, and the flushes from the two strong continental earthquakes in South Asia in 2001. Mainland China is the only continental sub-plate that is compressed by collision boundaries at its two sides, within which earthquakes are dispersive and distributed as seismic belts with different widths. The control capability of the continental block boundaries on the strong earthquakes and seismic hazards is calculated and analyzed in this dissertation. By mapping the distribution of the 31282 ML:3s2,0 earthquakes, I found that the depth of continental earthquakes depend on the tectonic zonings. The events on the boundaries of relatively integrated blocks are deep and those on the new-developed ruptures are shallow. The average depth of earthquakes in the West of China is about 5km deeper than that in the east. The western and southwestern brim of Tarim Basin generated the deepest earthquakes in mainland China. The statistic results from correlation between the grouped M7 earthquakes and the tidal stress show that the strong events were modulated by tidal stress in active periods. Taking Taiwan area as an example, the dependence of moderate events on the moon phase angles (￡>) is analyzed, which shows that the number of the earthquakes in Taiwan when D is 50° ,50° +90° and 50° +180° is more than 2 times of standard deviation over the average frequency at each degree, corresponding to the 4th, 12th and 19th solar day after the new moon. The probability of earthquake attack to the densely populated Taiwan island on the 4th solar day is about 4 times of that on other solar days. On the practice of earthquake prediction, I calculated and analyzed the temporal correlation of the earthquakes in Xinjinag area, Qinghai-Tibet area, west Yunnan area, North China area and those in their adjacent areas, and predicted at the end of 2000 that it would be a special time interval from 2001 to 2003, within which moderate to strong earthquakes would be more active in the west of China. What happened in 2001 partly validated the prediction. Within 10 months, there were 2 great continental earthquakes in south Asia, i.e., the M7.8 event in India on Jan 26 and M8.1 event in China on Nov. 14, 2001, which are the largest earthquake in the past 50 years both for India and China. No records for two great earthquakes in Asia within so short time interval. We should speculate the following aspects from the two incidences: The influence of the fallacious deployment of seismic stations on the fine location and focal mechanism determination of strong earthquakes must be affronted. It is very important to introduce comparative seismology research to seismic hazard analysis and earthquake prediction research. The improvement or changes in real-time prediction of strong earthquakes with precursors is urged. Methods need to be refreshed to protect environment and historical relics in earthquake-prone areas.