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.

冀东油田柳南地区新近系油藏剩余油形成及分布规律研究

This thesis is based on the research project of Study on the Geological Characteristics and Remaining Oil Distribution Law of Neogene Reservoirs in Liunan Area, which is one of the key research projects set by PetroChina Jidong Oilfield Company in 2006. The determination of remaining oil distribution and its saturation changes are the most important research contents for the development and production modification of oilfields in high water-cut phases. Liunan oilfield, located in Tangshan of Hebei Province geographically and in Gaoliu structural belt of Nanpu sag in Bohai Bay Basin structurally, is one of the earliest fields put into production of Jidong oilfield. Focusing on the development problems encountered during the production of the field, this thesis establishes the fine geological reservoir model through the study of reservoir properties such as fine beds correlation, sedimentary facies, micro structures, micro reservoir architecture, flow units and fluid properties. Using routine method of reservoir engineering and technology of reservoir numerical modeling, remaining oil distribution in the target beds of Liunan area is predicted successfully, while the controling factors of remaining oil distribution are illustrated, and the model of remaining oil distribution for fault-block structure reservoirs is established. Using staged-subdivision reservoir correlation and FZI study, the Strata in Liunan Area is subdivided step by step; oil sand body data-list is recompiled; diagram databases are established; plane and section configuration of monolayer sandstone body, and combination pattern of sandstone bodys are summarized. The study of multi-level staged subdivision for sedimentary micro-facies shows that the Lower member of Minghuazhen formation and the whole Guantao formation in Liunan Area belong to meandering river and braided river sedimentary facies respectively, including 8 micro facies such as after point bar, channel bar, channel, natural levee, crevasse splay, abandoned channel, flood plain and flood basin. Fine 3D geological modeling is performed through the application of advanced software and integration of geological, seismic logging and reservoir engineering data. High resolution numerical simulation is performed with a reserve fitting error less than 3%, an average pressure fitting fluctuation range lower than 2Mpa and an accumulate water cut fitting error less than 5%. In this way, the distribution law of the target reservoir in the study area is basically recognized. Eight major remaining oil distribution models are established after analysis of production status and production features in different blocks and different layers. In addition, fuzzy mathematics method is used to the integreted evaluation and prediction of abundant remaining oil accumulation area in major production beds and key sedimentary time units of the shallow strata in Liunan Area and corresponding modification comments are put forward. In summary, the establishment of fine reservoir geological model, reservoir numerical simulation and distribution prediction of remaining oil make a sound foundation for further stimulation of oilfield development performance.

冀东油田南堡滩海地区油气层损害机理与油气层保护方法研究及应用

Study and Application of Damage Mechanism and Protection Method of reservoir in Nanpu Shallow Beach Sea Area is one of the key research projects of Jidong Oilfield Company of PetroChina Company Limited from 2007 to 2008. Located at Nanpu Sag in Huanghua Depression of Bohaiwan Basin, Nanpu Shallow Beach Sea Area with 1000km2 exploration area posseses three sets, shallow Minghuazhen Formation and Guantao Formation of Upper Tertiary, middle-deep Dongying Formation of Lower Tertiary, deep Ordovician, of oil bearing series, according to the achievement of the connecting 3D seismic structure interpretation and the structural geological comprehensive research. Its main reservoir types include Upper Tertiary structural reservoir, Lower Tertiary structural and lithological-structural reservoir, and Ordovician ancient buried hill reservoir. How to protect reservoir, complete well and lift high efficiently is the key to realize high and stable yield of the oil wells during drilling, completing well, testing and repairing well. It is important for reservoir protecting during drilling that directly relate to efficient exploration. Therefore, beginning with basic characteristics and sensitive analysis of reservoir, study of reservoir damage machinism and analysis of reservoir damage potential factor are emphasized when prediction analysis about three-pressure profiles is carried out. The study both of physical and chemical properties and of the strata of the technology of borehole stabilization and reservoir protecting are outstanding. As the conclusions follow: （1）Based on the laboratory experiment about basalt cores, prediction of three- pressure profiles about 30 wells on No.1 and No.2 structure is practiced. The laws of plane pressure distribution are analyzed. （2）According to the analyses about reservoir feature data and about sensitivity evaluation to damage factor in Nanpu oil field, the scheme of reservoir protecting to the sand reservoir of Guantao Formation and the first section of Dongying Formation is put forward. （3）On basis of the analyses on lithological characteristics, mineral compositions, clay minerals, electrical behavior features, physical and chemical properties of basalt of Guantao formation in No.1 and No.2 structure, instability mechanism of basalt sidewall and technical countermeasures are obtained. （4）Aiming at the characteristics of Ordovician dissolution-pore fracture type carbonate reservoir, the scheme of the reservoir protecting to Ordovician is put forward. Creative study of the film forming and sealing and low invasion reservoir protection drilling fluid are successful. In summary, through the study of reservoir heterogeneity and sensitivity, a set of technology and schemes of reservoir protecting is put forward, which is adaptive during drilling the target bed in the research area and establishes the base for efficient exploration. Significant effect has showed in its application in Nanpu oil field.

反假频重建理论及其在储层预测中的应用

To extend the cross-hole seismic 2D data to outside 3D seismic data, reconstructing the low frequency data to high frequency data is necessary. Blind deconvolution method is a key technology. In this paper, an implementation of Blind deconvolution is introduced. And optimized precondition conjugate gradient method is used to improve the stability of the algorithm and reduce the computation. Then high-frequency retrieved Seismic data and the cross-hole seismic data is combined for constraint inversion. Real data processing proved the method is effective. To solve the problem that the seismic data resolution can’t meet the request of reservoir prediction in the river face thin-layers in Chinese eastern oil fields, a high frequency data reconstruction method is proposed. The extrema of the seismic data are used to get the modulation function which operated with the original seismic data to get the high frequency part of the reconstruction data to rebuild the wide band data. This method greatly saves the computation, and easy to adjust the parameters. In the output profile, the original features of the seismic events are kept, the common feint that breaking the events and adding new zeros to produce alias is avoided. And the interbeded details are enhanced compared to the original profiles. The effective band of seismic data is expended and the method is approved by the processing of the field data. Aim to the problem in the exploration and development of Chinese eastern oil field that the high frequency log data and the relative low frequency seismic data can’t be merged, a workflow of log data extrapolation constrained by time-phase model based on local wave decomposition is raised. The seismic instantaneous phase is resolved by local wave decomposition to build time-phase model, the layers beside the well is matched to build the relation of log and seismic data, multiple log info is extrapolated constrained by seismic equiphase map, high precision attributes inverse sections are produced. In the course of resolve the instantaneous phase, a new method of local wave decomposition --Hilbert transform mean mode decomposition(HMMD) is raised to improve the computation speed and noise immunity. The method is applied in the high resolution reservoir prediction in Mao2 survey of Daqing oil field, Multiple attributes profiles of wave impedance, gamma-ray, electrical resistivity, sand membership degree are produced, of which the resolution is high and the horizontal continuous is good. It’s proved to be a effective method for reservoir prediction and estimation.

渤海湾盆地前第三系构造演化与煤成气成藏研究

The disequilibrium between supply and demand the east part of North China accelerated natural gas exploration in Bohai bay basin. Exploration practice showed that coal-derived gas is important resource. In searching of big to middle scaled coal derived gas field, and realize successive gas supply, the paper carried out integrated study on structural evolution of Pre-Tertiary and evaluation of reservoir forming condition of coal-derived gas. Study work of the paper was based on the following condition: available achievement in this field at present, good understanding of multiphase of tectonic movement. Study work was focused on geological evolution, source rock evaluation and dissection key factors controlling reservoir forming. Based on analysis of seismic data, drilling data, tectonic style of Pre-Tertiary was subdivided, with different tectonic style representing different tectonic process. By means of state of the art, such as analysis of balanced cross section, and erosion restoration, the paper reestablished tectonic history and analyzed basin property during different tectonic phase. Dynamic mechanism for tectonic movement and influence of tectonic evolution on tectonic style were discussed. Study made it clear that tectonic movement is intensive since Mesozoic including 2 phase of compressional movement (at the end of Indo-China movement, and Yanshan movement), 2 phase of extensional movement (middle Yanshan movement, and Himalayan movement), 2 phase of strike slip movement, as well as 2 phase of reversal movement (early Yanshan movement, and early Himalayan movement). As a result, three tectonic provinces with different remnant of strata and different tectonic style took shape. Based on afore mentioned study, the paper pointed out that evolution of Bohai bay basin experienced the following steps: basin of rift valley type (Pt2+3)-craton basin at passive continental margin (∈1－2)-craton basin at active continental margin (∈3- O)-convergent craton basin (C-T1+2)-intracontinental basin (J＋K). Superposition of basins in different stage was discussed. Aimed at tectonic feature of multiple phases, the paper put forward concept model of superposition of tectonic unit, and analyzed its significance on reservoir forming. On basis of the difference among 3 tectonic movements in Mesozoic and Cenozoic, superposition of tectonic unit was classified into the following 3 categories and 6 types: continuous subsidence type (I), subsidence in Mesozoic and uplift for erosion in Cenozoic (II1), repeated subsidence and uplift in Mesozoic and subsidence in Cenozoic (II2), repeated subsidence and uplift in Mesozoic and uplift for erosion in Cenozoic (II3), uplift for erosion in Mesozoic and subsidence in Cenozoic (II4), and continuous uplift (III). Take the organic facies analysis as link, the paper established relationship between sedimentary environment and organic facies, as well as organic facies and organic matter abundance. Combined information of sedimentary environment and logging data, the paper estimated distribution of organic matter abundance. Combined with simulation of secondary hydrocarbon generation, dynamic mechanism of hydrocarbon generation, and thermal history, the paper made static and dynamic evaluation of effective source rock, i.e. Taiyuan formation and Shanxi formation. It is also pointed out that superposition of tectonic unit of type II2, type II4, and type I were the most favorable hydrocarbon generation units. Based on dissection of typical primary coal-derived gas reservoir, including reservoir forming condition and reservoir forming process, the paper pointed out key factors controlling reservoir forming for Carboniferous and Permian System: a. remnant thickness and source rock property were precondition; b. secondary hydrocarbon generation during Himalayan period was key factor; c. tectonic evolution history controlling thermal evolution of source rock was main factor that determine reservoir forming; d. inherited positive structural unit was favorable accumulation direction; e. fault activity and regional caprock determined hydrocarbon accumulation horizon. In the end, the paper established reservoir forming model for different superposition of tectonic units, and pointed out promising exploration belts with 11 of the first class, 5 of the second class and 6 of the third class.

滩海过渡带地震资料拼接处理方法研究

The technology of seismic acquirement has been more and more intact in both land area and marine area, however, since the acquirement equipment can not handle the tasks in this both two kinds of prospects at the same time, it is very hard to acquire in the prospect of seashore area, and also the quality of the raw data is poor, the related work in seashore area has been a tough task. With the development of the level of exploration, acquired equipment and technology, a certain series of processing methods has been built up, which enhanced the resolving power, the ratio of signal to noise and the quality of imaging, these methods could meet most of the requirements of exploration and the development. What is more, more and more nowadays prospects are composed by two or more 3D sub-prospect, the differences between frequencies, energy levels and the ratios of signal to noise are definitely obvious, we can not perform finely construct analysis and interpretation of reservoir without solving these problems in matching different raw data. As for the prospect that composed by both seashore and marine ones, through the analysis on the characteristics on frequency, amplitude, polarity and multiples of land prospect and marine prospect, the methods on how to eliminate those differences mentioned above were researched, and this paper gives a series of methods includes the computing of matching factor, recognizing of time difference, recognizing of phase difference, matching frequency, collapsing of abnormal energy, frequency compensating and so on, to fulfill the attenuations of energy difference, frequency difference, time difference and phase difference, and what is more, improve the innovative methods: includes the quantities method to synthetically and automatically recognize time difference and phase difference on the basis of signal match technology and according to statistic rule. And distinct the differences on phase and frequency characteristics of received raw data between velocity seismometer and piezoelectricity seismometer, and deduces that there are 90 degrees gap between data from these two seismometers, thirdly, synthetically apply zone abnormal frequency processing, surface consistent abnormal frequency compensate, residual abnormal frequency compensate, stack time equilibrium and optimize the processing seis-flow, deal with the would-be problems, such as low signal and noise rate, low resolution, poor imaging and poor state of consistent. Gain excellent effect in Da-gang oil field.