多波多分量地震资料走时近似方法的研究与应用

Compared with the conventional P wave, multi-component seismic data can markedly provide more information, thus improve the quality of reservoir evaluation like formation evaluation etc. With PS wave, better imaging result can be obtained especially in areas involved with gas chimney and high velocity formation. However, the signal-to-noise of multi-component seismic data is normally lower than that of the conventional P wave seismic data, while the frequency range of converted wave is always close to that of the surface wave which adds to the difficulty of removing surface wave. To realize common reflection point data stacking from extracted common conversion point data is a hard nut to crack. The s wave static correction of common receiver point PS wave data is not easy neither. In a word, the processing of multi-component seismic data is more complicated than P wave data. This paper shows some work that has been done, addressing those problems mentioned above. (1) Based on the AVO feature of converted wave, this paper has realized the velocity spectrum of converted waves by using Sarkar’s generalized semblance method taking into account of AVO factor in velocity analysis. (2)We achieve a method of smoothly offset division normal method.Firstly we scan the stacking velocities in different offset divisions for a t0, secondly obtain some hyperbolas using these stacking velocities, then get the travel time for every trace using these hyperbolas; in the end we interpolate the normal move out between two t0 for every trace. (3) Here realize a method of stepwise offset division normal moveout.It is similar to the method of smoothly offset division normal moveout.The main difference is using quadratic curve, sixth order curve or fraction curve to fit these hyperbolas. (4)In this paper, 4 types of travel time versus distance functions in inhomogeneous media whose velocity or slowness varies with depth and vertical travel time have been discussed and used to approximate reflection travel time. The errors of ray path and travel time based on those functions in four layered models were analyzed, and it has shown that effective results of NMO in synthetic or real data can be obtained. (5) The velocity model of converted PS-wave can be considered as that of P -wave based on the ghost source theory, thus the converted wave travel time can be approximated by calculation from 4 equivalent velocity functions: velocity or slowness vary linearly with depth or vertical travel time. Then combining with P wave velocity analysis, the converted wave data can be corrected directly to the P-wave vertical travel time. The improvements were shown in Normal Move out of converted waves with numerical examples and real data. (6) This paper introduces the methods to compute conversion point location in vertical inhomogeneous media based on linear functions of velocity or slowness versus depth or vertical travel time, and introduce three ways to choose appropriate equivalent velocity methods, which are velocity fitting, travel time approximation and semblance coefficient methods.

天然气藏地震识别技术研究

The Research of Seismic Recognition Techniques for Gas Reservoir Shang Yong_sheng(Geophysics) Directed by Yang Chang-chun Abstract Gas reservior is one of the most important nature resources. Someone forecast that the output will exceed crude oil in 2015 and become the largest energy source. Recently,more and more gas reservior are discovered as the oil field and gas filed exploration go deep into development. Although the gas proved reserves rise greatly the explorative degree of natural gas resource in our country is still very low. The potential of gas exploration is very great and our task is so hard. How to recognise and discover new gas reservoir is the first task based on the great gas reservior resources foreground. the gas reservior in different oil and gas field have its special gas generation, reservoiring, physical property conditions. However,it may have the same geophysical characters. So,it is very important to analyse, research, summarizing the geophysical characters of the gas reservior and make use of the characters to identify the gas layer effectively. This paper start with modeling,and it discuss the geophysical characters of the gas reservior response. It analyse the seismical wave characters of the gas reservoir. Furthermore, it summarize the method of using the seismica profile to identify the gas reservior directly. The paper discuss the research of extracting diffraction wave for mass diffraction wave grow at the edge of the gas reservoir at the seismic section. Making use of the technique of extracting diffraction wave to identify the gas reservior is the first experiment of the gas reservoir prediction technique. The avo technology is a new geophysical method. From the pre-stack analysis, this paper discuss the technique of using the rich information to identify the gas reservoir. Based on the case study of the Qidam basin and the Hailaer basin it discuss the method of predicating gas reservoir using pre-stack information. It include pre-stack amplitude preserve process, AVO modeling, fluid replacement technique, AVO analysis and interpretation technique. The paper summarize a gas reservoir prediction procedure focusing on the pre-stack information. The seismic wave will cause great attenuation when it pass through the gas layer and the high frequency component loss. This paper discuss the technique of extracting seismic attributes to represent the attenuation degree of seismica wave. Based on the attenuation attributes,it does the research of the gas reservor identification and prediction. At last,the paper discuss the method of calculating the azimuthal anisotropy to predict the fracture reservoir. Keyword: gas reservoir, diffraction wave, AVO, attenuation attribute,fracture prediction

油气勘探中的弱磁异常提取与区域场分析

Petroleum and natural gas is an important strategic resources. The short of the reserves will block the development of economy and threaten the safety of nation, along with the main oil fields of our country coming to the height of power and splendor of the exploitation and exploration. Therefore, it makes a great sense to inaugurate new explorative field and increase the reserves of petroleum and natural gas. Magnetic exploration is a main method of geophysics exploration. the developing observation apparatus and the perfect processing method provide wide space for magnetic exploration in these years. The method of magnetic bright spot is an application of magnetic exploration. The vertical migration of the hydrocarbon changes physical and chemical environment above the hydrocarbon reservoir, the new environment make tervalent iron translate into bivalent iron, that produce small scale magnetic anomaly, that is magnetic bright spot. The method of magnetic bright spot explores oil and gas field by the relation between the hydrocarbon and magnetic anomaly. This paper systemically research to pick-up and identify magnetic bright spot combining an oil field item, then point out advantaged area. In order to test the result, the author use the seismic information to superpose the magnetic bright spot, that prove the magnetic bright spot is reliable. then, the author complete a software to pick and identify the magnetic bright spot. The magnetic basement is very important to research forming and evolvement of the basin, especially, it is a crucial parameter of exploring residual basin in the research on pre-Cenozoic residual. This paper put forward a new method to inverse the interface of the magnetic layer on the basis of previous work, that is the method of separation of magnetic field step by step. The theory of this method is to translate the result of magnetic layer fluctuation to the result of magnetization density change, and the magnetic layer is flat, the paper choose thickness of magnetic layer as unit thickness, and define magnetic layer as a unit-thickness layer in order to convenient calculation, at the same time, define the variational magnetization density as equivalent magnetic density. Then we translate the relation between magnetic field and layer fluctuation to the relation between magnetic field and equivalent magnetic density, then, we can obtain the layer fluctuation through calculating equivalent magnetic density. Contrast to conventional parker method, model experimentation and example checkout prove this method is effective. The merit of this method is to avoid flat result in a strongly fluctuant area because of using a uniform average depth, the result of this method is closer to the fact, and this method is to inverse equivalent magnetic density, then translate equivalent magnetic density to layer fluctuation, this lays a foundation to inverse variational magnetic density in the landscape orientation and portrait.

改进叠前深度偏移的计算效率＆叠前高分辨率成像

An high-resolution prestack imaging technique of seismic data is developed in this thesis. By using this technique, the reflected coefficients of sheet sands can be gained in order to understand and identify thin oil reservoirs. One-way wave equation based migration methods can more accurately model seismic wave propagation effect such as multi-arrivals and obtain almost correct reflected energy in the presence of complex inhomogeneous media, and therefore, achieve more superiorities in imaging complex structure. So it is a good choice to apply the proposed high-resolution imaging to the presatck depth migration gathers. But one of the main shorting of one-way wave equation based migration methods is the low computational efficiency, thus the improvement on computational efficiency is first carried out. The method to improve the computational efficiency of prestack depth migration is first presented in this thesis, that is frequency-dependent varying-step depth exploration scheme plus a table-driven, one-point wavefield interpolation technology for wave equation based migration methods; The frequency-dependent varying-step depth exploration scheme reduces the computational cost of wavefield depth extrapolation, and the a table-driven, one-point wavefield interpolation technology reconstructs the extrapolated wavefield with an equal, desired vertical step with high computational efficiency. The proposed varying-step depth extrapolation plus one-point interpolation scheme results in 2/3 reduction in computational cost when compared to the equal-step depth extrapolation of wavefield, but gives the almost same imaging. The frequency-dependent varying-step depth exploration scheme is presented in theory by using the optimum split-step Fourier. But the proposed scheme can also be used by other wave equation based migration methods of the frequency domain. The proposed method is demonstrated by using impulse response, 2-D Marmousi dataset, 3-D salt dataset and the 3-D field dataset. A method of high-resolution prestack imaging is presented in the 2nd part of this thesis. The seismic interference method to solve the relative reflected coefficients is presented. The high-resolution imaging is obtained by introducing a sparseness- constrained least-square inversion into the reflected coefficient imaging. Gaussian regularization is first imposed and a smoothed solution is obtained by solving equation derived from the least-square inversion. Then the Cauchy regularization is introducing to the least-square inversion , the sparse solution of relative reflected coefficients can be obtained, that is high-resolution solution. The proposed scheme can be used together with other prestack imaging if the higher resolution is needed in a target zone. The seismic interference method in theory and the solution to sparseness-constrained least-square inversion are presented. The proposed method is demonstrated by synthetic examples and filed data.

非均质储层叠前地震反演技术研究及应用

(1) I research on the relationship between elastic parameters, lithology and liquid. It is a physical base for pre-stack seismic inversion. I research all kinds of approximate expressions of Zoeppritz function. Then the relation of all kinds of approximate expressions can be confirmed. The geological model of water sand and gas sand in different depth was designed. Moreover I research on precision of all kinds of approximate expressions. (2) In process of seismic data which aim at amplitude recovery and apply in pre-stack seismic inversion, I advance to adopt double flow chart for different aim. Pre-stack noise elimination, real amplitude recovery and NMO correction of long offset are the key taches. (3) I made a systemic expatiate for the thinking and applicability about all kinds of expressions of elastic impedance. And mathematical model was applied to compare the precision with all kinds of expressions of elastic impedance. I propose a new pre-stack simultaneous inversion which is based on the Zoppritz function and simulated annealing algorithm. This method can ensure calculation precision of reflection coefficient from different incident angle and get a global optimum solution. Therefore this method improves the precision of pre-stack seismic inversion. (4) The object function of P-S wave pre-stack simultaneous inversion was established. I compared the precision and convergence between simultaneous inversion and P-wave inversion. And the results show that simultaneous inversion is superior to P-wave inversion. Through the study of AVO event of transformed wave, AVO characters of different kinds of gas sand were analyzed. (5) I carried out the study work of pre-stack seismic inversion for carbonate reservoir in middle of Tarim basin and sand shale reservoir in Sulige Area of Erdos Basin. The method and technology in this paper was applied to practical work. And I made a prediction for heterogeneous reservoir. Moreover it acquires a good application effect. Key Word: reflection coefficient， amplitude recovery， pre-stack seismic inversion， Heterogeneous reservoir，prediction.

Fourier有限差分法深度偏移与正演模拟

With the development of seismic exploration, the target becomes more and more complex, which leads to a higher demand for the accuracy and efficiency in 3D exploration. Fourier finite-difference (FFD) method is one of the most valuable methods in complex structure exploration, which keeps the ability of finite-differenc method in dealing with laterally varing media and inherits the predominance of the phase-screen method in stablility and efficiency. In this thesis, the accuracy of the FFD operator is highly improved by using simulated annealing algorithm. This method takes the extrapolation step and band width into account, which is more suitable to various band width and discrete scale than the commonely-used optimized method based on velocity contrast alone. In this thesis, the FFD method is extended to viscoacoustic modeling. Based on one-way wave equation, the presented method is implemented in frequency domain; thus, it is more efficient than two-way methods, and is more convenient than time domain methods in handling attenuation and dispersion effects. The proposed method can handle large velocity contrast and has a high efficiency, which is helpful to further research on earth absorption and seismic resolution. Starting from the frequency dispersion of the acoustic VTI wave equation, this thesis extends the FFD migration method to the acoustic VTI media. Compared with the convetional FFD method, the presented method has a similar computational efficiency, and keeps the abilities of dealing with large velocity contrasts and steep dips. The numerical experiments based on the SEG salt model show that the presented method is a practical migration method for complex acoustical VTI media, because it can handle both large velocity contrasts and large anisotropy variations, and its accuracy is relatively high even in strong anisotropic media. In 3D case, the two-way splitting technique of FFD operator causes artificial azimuthal anisotropy. These artifacts become apparent with increasing dip angles and velocity contrasts, which prevent the application of the FFD method in 3D complex media. The current methods proposed to reduce the azimuthal anisotropy significantly increase the computational cost. In this thesis, the alternating-direction-implicit plus interpolation scheme is incorporated into the 3D FFD method to reduce the azimuthal anisotropy. By subtly utilizing the Fourier based scheme of the FFD method, the improved fast algorithm takes approximately no extra computation time. The resulting operator keeps both the accuracy and the efficiency of the FFD method, which is helpful to the inhancements of both the accuracy and the efficiency for prestack depth migration. The general comparison is presented between the FFD operator and the generalized-screen operator, which is valuable to choose the suitable method in practice. The percentage relative error curves and migration impulse responses show that the generalized-screen operator is much sensiutive to the velocity contrasts than the FFD operator. The FFD operator can handle various velocity contrasts, while the generalized-screen operator can only handle some range of the velocity contrasts. Both in large and weak velocity contrasts, the higher order term of the generalized-screen operator has little effect on improving accuracy. The FFD operator is more suitable to large velocity contrasts, while the generalized-screen operator is more suitable to middle velocity contrasts. Both the one-way implicit finite-difference migration and the two-way explicit finite-differenc modeling have been implemented, and then they are compared with the corresponding FFD methods respectively. This work gives a reference to the choosen of proper method. The FFD migration is illustrated to be more attractive in accuracy, efficiency and frequency dispertion than the widely-used implicit finite-difference migration. The FFD modeling can handle relatively coarse grids than the commonly-used explicit finite-differenc modeling, thus it is much faster in 3D modeling, especially for large-scale complex media.

济阳坳陷深层地质构造与演化

Based on geophysical and geological data in Jiyang depression, the paper has identified main unconformity surfaces (main movement surfaces) and tectonic sequences and established tectonic and strata framework for correlation between different sags. Based on different sorts of structural styles and characteristics of typical structures, the paper summarized characteristics and distribution of deep structures, discussed evolution sequence of structure, analyzed the relation between tectonic evolution and generation of petroleum. The major developments are as following: Six tectonic sequences could be divided from bottom to top in the deep zone of Jiyang depression. These tectonic sequences are Cambrian to Ordovician, Carboniferous to Permian, lower to middle Jurassic, upper Jurassic to lower Cretaceous, upper Cretaceous and Kongdian formation to the fourth member of Shahejie formation. The center of sedimentation and subsidence of tectonic sequences distinguished from each other in seismic profiles is controlled by tectonic movements. Six tectonic evolution stages could be summarized in the deep zone in Jiyang depression. Among these stages, Paleozoic stage is croton sedimentation basin; Indosinian stage, open folds of EW direction are controlled by compression of nearly SN direction in early Indosinian (early to middle Triassic) while fold thrust fault of EW – NWW direction and arch protruding to NNE direction are controlled by strong compression in late Indosinian (latter Triassic); early Yanshanian stage (early to middle Jurassic), in relatively weak movement after Indosinian compressional orogeny, fluviolacustrine is deposited in intermontane basins in the beginning of early Yanshanian and then extensively denudated in the main orogenic phase; middle Yanshanian (late Jurassic to early Cretaceous), strike-slipping basins are wide distribution with extension (negative reversion) of NW – SE direction; latter Yanshanian (late Cretaceous), fold and thrust of NE – NNE direction and positive reversion structure of late Jurassic to early Cretaceous strike-slipping basin are formed by strong compression of NW–SE direction; sedimentation stage of Kongdian formation to the fourth member of Shahejie formation of Cenozoic, half graben basins are formed by extension of SN direction early while uplift is resulted from compression of nearly EW direction latterly. Compression system, extension system and strike-slip system are formed in deep zone of Jiyang depression. According to identifying flower structure of seismic profiles and analysis of leveling layer slice of 3D seismic data and tectonic map of deep tectonic interface, strike-slip structures of deep zone in Jiyang depression are distinguished. In the middle of the Jiyang depression, strike-slip structures extend as SN direction, NNW direction in Huimin sag, but NNE in Zhandong area. Based on map of relict strata thickness, main faults activity and regional tectonic setting, dynamic mechanisms of deep structure are preliminary determination. The main reason is the difference of direction and character of the plate’s movement. Development and rework of multi-stage tectonic effects are benefit for favorable reservoir and structural trap. Based on tectonic development, accumulation conditions of deep sub-sags and exploration achievements in recent years, potential zones of oil-gas reservoir are put forward, such as Dongying sag and Bonan sag.

残留盆地分布综合地球物理研究——以大港地区为例

Macro-distribution of residual basins is a basic question in residual basin research，the main object of macro-distribution study is to build strata framework, compute thickness of residual strata and analyze characteristics of residual basins. With the guidance of the theory of integrated geology and geophysical research, the paper assembled series of methods and established the technical chart based on gravity and magnetic data, with restriction of geology, seismic and drilling data. Based on potential field data processing and analysis, forward and inverse computation, region potential field analysis and potential field separation, etc. it computed depth of gravity/magnetic basement and got strata framework. It had got effective results in the research of macro-distribution of residual basin research in the Dagang area. It did the wavelet transform of gravity/magnetic data with multi-kind of wavelet basis using a trou algorithm. From comparison of processing result and their spectral of wavelet analysis, up continuation and filter method, the wavelet approximation is better to fit the regional potential field, and it is an effective method to separate gravity/magnetic effect caused by deep geology bodies. The experiment of matching pursuit shows that te transform domain methods have great advantage in potential data analysis. From the integrated geophysical study of rock property study, gravity/magnetic basement inversion and fault system analysis of the Dagang area, it gets the strata framework and the thickness of pre-Cenozoic residual strata. Comprehensive study with gravity and magnetotelluric profile inversion and interpretation, three prospect plays of macro-distribution of residual basins are fingered out. It has great residual strata thickness in the northern part of Chengning Uplift and there is thrust fault in the deep zone and good up-Paleozoic hydrocarbon source rocks in this area. With integrated analysis, this area will be the most prospective hydrocarbon location of pre-Cenozoic residual basins.

基于显式分形插值的高保真地震数据重建方法研究

Based on the fractal theories, contractive mapping principles as well as the fixed point theory, by means of affine transform, this dissertation develops a novel Explicit Fractal Interpolation Function（EFIF）which can be used to reconstruct the seismic data with high fidelity and precision. Spatial trace interpolation is one of the important issues in seismic data processing. Under the ideal circumstances, seismic data should be sampled with a uniform spatial coverage. However, practical constraints such as the complex surface conditions indicate that the sampling density may be sparse or for other reasons some traces may be lost. The wide spacing between receivers can result in sparse sampling along traverse lines, thus result in a spatial aliasing of short-wavelength features. Hence, the method of interpolation is of very importance. It not only needs to make the amplitude information obvious but the phase information, especially that of the point that the phase changes acutely. Many people put forward several interpolation methods, yet this dissertation focuses attention on a special class of fractal interpolation function, referred to as explicit fractal interpolation function to improve the accuracy of the interpolation reconstruction and to make the local information obvious. The traditional fractal interpolation method mainly based on the randomly Fractional Brown Motion (FBM) model, furthermore, the vertical scaling factor which plays a critical role in the implementation of fractal interpolation is assigned the same value during the whole interpolating process, so it can not make the local information obvious. In addition, the maximal defect of the traditional fractal interpolation method is that it cannot obtain the function values on each interpolating nodes, thereby it cannot analyze the node error quantitatively and cannot evaluate the feasibility of this method. Detailed discussions about the applications of fractal interpolation in seismology have not been given by the pioneers, let alone the interpolating processing of the single trace seismogram. On the basis of the previous work and fractal theory this dissertation discusses the fractal interpolation thoroughly and the stability of this special kind of interpolating function is discussed, at the same time the explicit presentation of the vertical scaling factor which controls the precision of the interpolation has been proposed. This novel method develops the traditional fractal interpolation method and converts the fractal interpolation with random algorithms into the interpolation with determined algorithms. The data structure of binary tree method has been applied during the process of interpolation, and it avoids the process of iteration that is inevitable in traditional fractal interpolation and improves the computation efficiency. To illustrate the validity of the novel method, this dissertation develops several theoretical models and synthesizes the common shot gathers and seismograms and reconstructs the traces that were erased from the initial section using the explicit fractal interpolation method. In order to compare the differences between the theoretical traces that were erased in the initial section and the resulting traces after reconstruction on waveform and amplitudes quantitatively, each missing traces are reconstructed and the residuals are analyzed. The numerical experiments demonstrate that the novel fractal interpolation method is not only applicable to reconstruct the seismograms with small offset but to the seismograms with large offset. The seismograms reconstructed by explicit fractal interpolation method resemble the original ones well. The waveform of the missing traces could be estimated very well and also the amplitudes of the interpolated traces are a good approximation of the original ones. The high precision and computational efficiency of the explicit fractal interpolation make it a useful tool to reconstruct the seismic data; it can not only make the local information obvious but preserve the overall characteristics of the object investigated. To illustrate the influence of the explicit fractal interpolation method to the accuracy of the imaging of the structure in the earth’s interior, this dissertation applies the method mentioned above to the reverse-time migration. The imaging sections obtained by using the fractal interpolated reflected data resemble the original ones very well. The numerical experiments demonstrate that even with the sparse sampling we can still obtain the high accurate imaging of the earth’s interior’s structure by means of the explicit fractal interpolation method. So we can obtain the imaging results of the earth’s interior with fine quality by using relatively small number of seismic stations. With the fractal interpolation method we will improve the efficiency and the accuracy of the reverse-time migration under economic conditions. To verify the application effect to real data of the method presented in this paper, we tested the method by using the real data provided by the Broadband Seismic Array Laboratory, IGGCAS. The results demonstrate that the accuracy of explicit fractal interpolation is still very high even with the real data with large epicenter and large offset. The amplitudes and the phase of the reconstructed station data resemble the original ones that were erased in the initial section very well. Altogether, the novel fractal interpolation function provides a new and useful tool to reconstruct the seismic data with high precision and efficiency, and presents an alternative to image the deep structure of the earth accurately.

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

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

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

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).

博兴洼陷新生代构造-沉积演化及其对鲁西隆起的响应

The West Shandong Uplift and its adjacent basins, with same evolutional history before Mesozoic, are an important basin-orogenic systems in North China. After late Mesozoic, tectonic differentiation between basin and orogenic belt gradually displayed in the study area. The Boxing sag is a part of Jiyang Depression near to West Shandong Uplift, in which the whole Mesozoic and Cenozoic strata are preserved. Based on the analysis of sedimentary records in the Boxing sag, the Cenozoic structural and sedimentary evolutions in Boxing Sag and its response to Western Shandong uplift are discussed in this dissertation. The main conclusions in this research are presented as follows. Based on Seismic and well logging profile interpretation, fault growth index, thickness difference between bottom wall and top wall and fault activity rate from Eocene to Pliocene are studied. Boxing sag had three main faults, NE, NW and NEE trending faults. Research shows that the activity of the NW trending fault in the Boxing sag became weaken from E1-2S4 to N2m gradually. The evolution of NE and the NEE trending fault can be divided into three episodes, from E1-2k to E2s4, from E2s3 to E3s1, from N2m to E3d. The analysis of Paleogene samples of heavy mineral assemblages shows that metamorphic rocks represented by garnet, intermediate-acid igneous rocks represented by the assemblage of apatite, zircon and tourmaline became less from E1-2k to N2g, and sedimentary rocks represented by the assemblage of pyrite, barite and limonite also became less. Intermediate-basic igneous rocks represented by the assemblage of leucoxene, rutile and ilmenite and metamorphic rocks represented by epidote became more and more. Electronic microprobe analysis shows that glaucophane and barroisite are existed in Kongdian Formation and the 4th member of Shahejie Formation, and they demonstrate that Western Shandong and Eastern Shandong are all the source regions of the Boxing Sag, and they also indicate that oceanic crust existed before the collision between the Yangtze and North China continent. The fact that Eastern Shandong is the source region of Boxing Sag also indicates that Western Shandong was not high enough to prevent sediment from Eastern Shandong at E1-2k and E2s4. The results of the dating of five detrital zircons of Boxing Sag show Kongdian Formation and the 4th member of Shahejie Formation have the age peaks of 2800Ma and 700-800. It means that Eastern Shandong is the source region of Boxing Sag at early Paleogene and Western Shandong is not high enough to prevent the sediment from Eastern Shandong. The ages of 160-180 and 220-260 Ma, which exist in the Guantao Formation and Paleogene, are common in Eastern Shandong and rare in Western Shandong,and it implied that Western Shandong is a low uplift at 24Ma. The Paleogene strata have almost same age groups, while the Guantao Formation has significant variations of age groups, and this indicates that Boxing Sag and Western Shandong uplift had taken place tremendous changes. The results of apatite fission track in Boxing sag show that three times uplifts happened at the source region at 60 Ma, 45Ma and 15Ma respectively, and the Boxing sag experienced two subsidences at 60Ma, 45Ma and one uplift at 20Ma.

准噶尔盆地滴南凸起中生界层序地层及有利勘探目标优选

Based on the theories of sequence stratigraphy and sedimentology, as well as comprehensive studies of seismic data, drilling data, core interpretation and setting of this area, the thesis presents an analysis for Mesozoic formation in Dinan uplift. By means of recognizing the boundary of the sequence, dividing and correlating the systems tract, Mesozoic of Dinan uplift is divided into ten sequences and twenty-five systems tracts during the establishment of the sequence framework. In the framework, some sequences are featured by mature systems of lowstand, water-transgression and highstand, while some undeveloped systems of lowstand or highstand. The main sedimentary facies in Mesozoic of Dinan uplift are braided river, meandering river, delta and lake. The braided river was divided into sandy river and rudaceous river by the lithology of the river channel and was divided into dry climate and wet climate condition by the color of the flood plain. Additionally, The concept of “wetland” is put forward for the first time and regarded as the consequence of wet climate. The analysis includes the classification of six types of traps: (1) stratigraphic overlap trap, (2) lithologic trap with updip pinchout, (3) stratigraphic unconformity trap, (4) fault-lithology trap, (5) fault trap, (6) anticlinal trap, and combining with the research of the characteristics and distribution rules for the known reservoir, it draws out that “fault control” is the petroleum accumulation pattern in this area, in which fault is the key element of the transporting system. Finally the thesis concludes the distribution characteristics and optimized some targets for the potential exploration zone.

松辽盆地长岭断陷火山岩气藏的形成条件与分布规律

Changling fault depression is a compound fault depression complicated by interior fault, with faults in the west and overlap in the west. North of Changling fault depression show NNE strike while south is NW strike. Changling fault depression has undergone twochasmic stage which control the development and distribution of volcanic rock, one depression stage, later inversion and uplift stage which control the formation of natural gas reservoir, and basin atrophic stage. The main boundary faults and main faults in Changling fault depression control three volcanic cycles and the distribution of volcanic rock. Seismic reflection characteristic and logging response characteristic of volcanic rock in study area are obvious, and the distribution characteristic, volcanic cycle and active stage of volcanic rock can be revealed by seismic attribute, conventional logging data can distinguish clastic rock from volcanic rock or distinguish partial different types of volcanic rock. The reservoir property of rhyolite and volcanic tuff are the best. Favorable volcanic reservoir can be preserved in deep zone. Imaging logging and frequency decompostion technology of seismic data act as effective role in the study of reservoir physical property and gas-bearing properties of volcanic rock.. Hydrocarbon gas in study area is high and over mature coal type gas, the origin of CO2 is complex, it is either inorganic origin or organic origin, or mixing origin. Hydrocarbon gas is mainly originate from Shahezi formation and Yingcheng formation source rocks, CO2 is mainly mantle source gas. Hydrocarbon has the characteristics of continuous accumulation with two charging peak. The first peak represent liquid hydrocarbon accumulation time, The second peak stand for the accumulation time of gaseous hydrocarbon.CO2 accumulate approximately in Neocene. The source rock distribution range, volcanic rock and favorable reservoir facies, distribution characteristic of deep fault (gas source fault) and late inversion structure are the major factors to control gas reservoir formation and distribution. All the results show that these traps that consist of big inherited paleo uplift(paleo slope), stratigraphic overlap and thinning out, volcanic rock, are the most advantageous target zone.