Petroleum System in Deepwater Basins of the Northern South China Sea

The northern South China Sea margin has experienced a rifting stage and a post-rifting stage during the Cenozoic. In the rifting stage, the margin received lacustrine and shallow marine facies sediments. In the post-rifting thermal subsidence, the margin accumulated shallow marine facies and hemipelagic deposits, and the deepwater basins formed. Petroleum systems of deepwater setting have been imaged from seismic data and drill wells. Two kinds of source rocks including Paleogene lacustrine black shale and Oligocene-Early Miocene mudstone were developed in the deepwater basin of the South China Sea. The deepwater reservoirs are characterized by the deep sea channel rill, mass flow complexes and drowned reef carbonate platform. Profitable capping rocks on the top are mudstones with huge thickness in the post-rifting stage. Meanwhile, the faults developed during the rifting stage provide a migration path favourable for the formation of reservoirs. The analysis of seismic and drilling data suggests that the joint structural and stratigraphic traps could form giant hydrocarbon fields and hydrocarbon reservoirs including syn-rifting graben subaqueous delta, deepwater submarine fan sandstone and reef carbonate reservoirs.

Seismic stratigraphy of the Qiongdongnan deep sea channel system, northwest South China Sea

Based on more than 4000 km 2D seismic data and seismic stratigraphic analysis, we discussed the extent and formation mechanism of the Qiongdongnan deep sea channel. The Qiongdongnan deep sea channel is a large incised channel which extends from the east boundary of the Yinggehai Basin, through the whole Qiongdongnan and the Xisha trough, and terminates in the western part of the northwest subbasin of South China Sea. It is more than 570 km long and 4-8 km wide. The chaotic (or continuous) middle (or high) amplitude, middle (or high) continuity seismic facies of the channel reflect the different lithological distribution of the channel. The channel formed as a complex result of global sea level drop during early Pliocene, large scale of sediment supply to the Yinggehai Basin, inversion event of the Red River strike-slip fault, and tilted direction of the Qiongdongnan Basin. The large scale of sediment supply from Red River caused the shelf break of the Yinggehai Basin to move torwards the S and SE direction and developed large scale of prograding wedge from the Miocene, and the inversion of the Red River strike-slip fault induced the sediment slump which formed the Qiongdongnan deep sea channel.

Analysis on deep crustal structure along the onshore-offshore seismic profile across the Binhai (littoral) fault zone in northeastern South China Sea

The onshore-offshore deep seismic experiment was carried out for the first time and filled the blankness of the seismic surveys in the transition area between South China and northeastern South China Sea. The seismic data were analyzed and processed. The different seismic phases were identified and their travel time arrivals were modeled by ray-tracing to study the P-wave velocity crustal structure of this area. The crustal structure of this area is the continental crust. The crust thickness is gradually decreasing southward along the on-shore-offshore seismic line. The low-velocity layer (5.5 similar to 5.9 km (.) s(-1)) exists generally in the middle crust (about 10.0 similar to 18.0km)with about 2.5 similar to 4.0 km thickness, which is also thinning seaward. No obvious high-velocity layer appears in the lower crust. The Binhai (littoral) fault zone is a low velocity zone, which is located about 35km southeast to the Nan'ao station and corresponding to the gradient belt of gravity & magnetism anomalies. The depth of the fault zone is close to the Moho discontinuity. The littoral fault zone is a boundary between the normal continental crust of South China and the thinned continental crust of the sea area.

用光学阵列传感器的机器人物体分类系统

本文介绍用光学阵列传感器的机器人物体分类系统。传感器直接安装在机器人的两个手指上。被抓物体的阴影通过光导纤维传到安放在“安全区”的光敏元件上。计算机识别物体的轮廓后命令机器人抓握物体,并把它运送到指定的地点从而达到物体分类的目的。

随钻地震技术的研究及应用

Seismic While Drilling (SWD) is a new wellbore seismic technique. It uses the vibrations produced by a drill-bit while drilling as a downhole seismic energy source. The continuous signals generated by the drill bit are recorded by a pilot sensor attached to the top of the drill-string. Seismic wave receivers positioned in the earth near its surface receive the seismic waves both directly and reflection from the geologic formations. The pilot signal is cross-correlated with the receiver signals to compute travel-times of the arrivals (direct arrival and reflected arrival) and attenuate incoherent noise. No downhole intrusmentation is required to obtain the data and the data recording does not interfere with the drilling process. These characteristics offer a method by which borehole seismic data can be acquired, processed, and interpreted while drilling. As a Measure-While-Drill technique. SWD provides real-time seismic data for use at the well site . This can aid the engineer or driller by indicating the position of the drill-bit and providing a look at reflecting horizons yet to be encountered by the drill-bit. Furthermore, the ease with which surface receivers can be deployed makes multi-offset VSP economically feasible. First, this paper is theoretically studying drill-bit wavefield， interaction mode between drill-bit and formation below drill-bit , the new technique of modern signal process was applied to seismic data, the seismic body wave radiation pattern of a working roller-cone drill-bit can be characterized by theoretical modeling. Then , a systematical analysis about the drill-bit wave was done, time-distance equation of seismic wave traveling was established, the process of seismic while drilling was simulated using the computer software adaptive modeling of SWD was done . In order to spread this technique, I have made trial SWD modeling during drilling. the paper sketches out the procedure for trial SWD modeling during drilling , the involved instruments and their functions, and the trial effect. Subsurface condition ahead of the drill-bit can be predicted drillstring velocity was obtained by polit sensor autocorrelation. Reference decovolution, the drillstring multiples in the polit signal are removed by reference deconvolution, the crosscorrelation process enhance the signal-to-noise power ratio, lithologies. Final, SWD provides real-time seismic data for use at the well site well trajectory control exploratory well find out and preserve reservoirs. intervel velocity was computed by the traveltime The results of the interval velocity determination reflects the pore-pressure present in the subsurface units ahead of the drill-bit. the presences of fractures in subsurface formation was detected by shear wave. et al.

大港油田复杂断块油气藏评价技术研究与应用

In this paper, in view of characteristics and appraisal difficulty of complicate fault block reservoir, based on the theory of reservoir geology, tectonic geology, sequence stratigraphy, logging geology and sedimentology, according to related multidisciplinary data, such as geological, geophysical, logging and testing, taking 3D seismic acquisition processing, delicate structure interpretation, fine reservoir prediction, reservoir engineering and logging evaluation as research methods, on the basis of comprehensive geologic research, innovative appraisal thinking of complicate fault block reservoir is probed into, reservoir evaluation technology series of different reservoirs were created, and appraisal circuit of complicate fault block reservoir was proposed. According to the research on main controlling factors of hydrocarbon and concentration and concealing laws of Dagang exploration area, family assemble patterns of complicate fault block reservoir were proposed firstly, such as twinning, juxtaposed, overlying, concurrent, connected in series and so on, so theory basement was provided to reservoir evaluation. Taking into account of difficulties of low signal to noise ratio and border identification of lithologic trap of complicate fault block reservoir, in the middle of seismic data acquisition, double focusing layout and optimization technology were adopted, geophysical model was set up, and best observation system was ascertained. During object processing, matching processing technology of removing time difference, phase difference, energy difference and frequency spectrum difference generated by different blazing and receiving elements was studied. Imagery, low amplitude structure and identification of subtle lithologic reservoir was interpreted reasonably. On the basis of characteristics and technical description analysis of structural, structural--lithologic and lithologic reservoirs, innovative appraisal thinking of these reservoirs was formatted, appraisal circuit of complicate fault block reservoir was proposed, and this method could be used in other similar oilfields

裂隙/孔隙型弹性介质中地震波场模拟与叠前逆时偏移

With the development of oil/gas seismic exploration, seismic survey for fracture/porosity type reservoir is becoming more and more important. As for China, since it has over 60% store of low porosity and low permeability oil/gas reservoir, it’s more urgent to validly describe fracture/porosity type oil/gas trap and proposing the related, developed seismic technique. To achieve mapping fracture/porosity region and its development status, it demands profound understanding of seismic wave propagation discipline in complex fractured/pored media. Meanwhile, it has profound scientific significance and applied worth to study forward modeling of fracture/porosity type media and pre-stacked reverse time migration. Especially, pre-stacked reverse-time migration is the lead edge technique in the field of seismology and seismic exploration. In this paper, the author has summarized the meaning, history and the present state of numerical simulation of seismic propagation in fractured/pored media and seismic exploration of fractured/pored reservoirs. Extensive Dilatancy Anisotropy (EDA) model is selected as media object in this work. As to forward modeling, due to local limitation of solving spatial partial derivative when using finite-difference and finite-element method, the author turns to pseudo-spectral method (PSM), which is based on the global characteristic of Fourier transform to simulate three-component elastic wave-field. Artifact boundary effect reduction and simulation algorithm stability are also discussed in the work. The author has completed successfully forward modeling coding of elastic wave-field and numerical simulation of two-dimensional and three-dimensional EDA models with different symmetric axis. Seismic dynamic and kinematical properties of EDA media are analyzed from time slices and seismic records of wave propagation. As to pre-stacked reverse-time migration for elastic wave-field in fractured/pored media, based on the successful experience in forward modeling results with PSM, the author has studied pre-stacked reverse-time depth-domain migration technique using PSM of elastic wave-field in two dimensional EDA media induced by preferred fracture/pore distribution. At the same time, different image conditions will bring up what kind of migration result is detailed in this paper. The author has worded out software for pre-stacked reverse-time depth-domain migration of elastic wave-field in EDA media. After migration processing of a series of seismic shot gathers, influences to migration from different isotropic and anisotropy models are described in the paper. In summary, following creative research achievements are obtained:  Realizing two-dimensional and three-dimensional elastic wave-field modeling for fractured/pored media and related software has been completed.  Proposed pre-stacked reverse-time depth-domain migration technique using PSM of elastic wave-field.  Through analysis of the seismic dynamic and kinematical properties of EDA media, the author made a conclusion that collection of multi-component seismic data can provide important data basis for locating and describing the fracture/pore regions and their magnitudes and the preferred directions.  Pre-stacked reverse-time depth-domain migration technique has the ability to reconstruct complex geological object with steep formations and tilt fracture distribution. Neglecting seismic anisotropy induced by the preferred fracture/pore distribution, will lead to the disastrous imaging results.

裂缝性油气藏储层预测方法研究

In this paper, based on the E & P situation in the oilfield and the theory of geophysical exploration, a series researches are conducted on fracture reservoir prediction technology in general，and it especially focus on some difficult points. The technological series which integrated amplitude preserved data processing、interpretation and its comprehensive application research as a whole were developed and this new method can be applied to the other similar oilfield exploration and development. The contents and results in this paper are listed as follows: 1. An overview was given on the status and development of fracture reservoir estimation technique, compare and analyze those geophysical prediction methods. This will be very helpful to the similar reservoir researches. 2. Analyze and conclude the characters of geologies and well logging response of burial hills fracture reservoir, those conclusions are used to steer the geophysical research and get satisfying results. 3. Forward modeling anisotropy seismic response of fracture reservoir. Quantitatively describe the azimuthal amplitude variation. Amplitude ellipse at each incidence angle is used to identify the fracture orientation. 4. Numerical simulation of structure stress based on finite difference method is carried out. Quantitatively describe and analyze the direction and intensity of fracture. 5. Conventional attributes extraction of amplitude preserved seismic data、attributes with different azimuthal angle and different offset are used to determine the relationship between the results and fracture distribution. 6. With spectrum decomposition method based on wavelet transform， the author disclose the reservoir distribution in space. It is a powerful tool to display its anisotropy. 7. Integrated seismic wave impendence、elastic impendence、spectrum decomposition、attribute extraction、fracture analysis result as a whole to identify and evaluate the fracture reservoir. An optimum workflow is constructed. It is used to practical oil&gas production and good results are obtained. This can indicate the wide foreground of this technique series.

储层评价目标驱动的地震叠前叠后综合反演方法研究

The exploration and development of natural gas in the north of Ordos basin have been one important part in China’s energy stratagem. Reservoir in upper Palaeozoic group is of lithological trap and its prediction is a crux in a series of works. Based on foregoing seismic reservoir prediction, seismic data are re-processed with some optical methods and pre-stack information is used in corresponding inversions. Through the application of diverse methods, a series of techniques for reservoir prediction come into being. Several results are achieved as flowing: 1. A set of log processing and interpretation methods is developed. Porosity, permeability and gas saturation models are rebuilt. 2. Based on the petro-physics analysis of reservoirs in upper Palaeozoic group, the equations about lithology, property, hydrocarbon and elastic parameters are established. 3. Forward modeling based on elastic wave theory is first applied in the study area and increases the resolution of modeling results. 4. A series of techniques such as pre-stack time migration and others are combined to improve the data quality. 5. Pre-stack seismic inversion is first employed in the north of Ordos Basin and brings the results of EI, P-impedance, S-impedance and other elastic parameters. 6. In post-stack inversion, logs indicating reservoir parameters are rebuilt and boost the resolution of lithology inversion. 7. Amplitude, coherence, frequency-discomposed amplitude, waveform and other sensitive attributes are extracted to describe sands distribution. Seismic modes standing for sands of P1x3, P1x2 are established. 8. Among 9 proposed wells, 8 wells encountered sands and became production wells. The output of DK13 amounts to 510,000 m3 per day. Keywords：the north of Ordos Basin, reservoir prediction, pre-stack inversion, post-stack inversion, seismic attributes．

叠前叠后联合反演技术及其在QHD32-6油田的应用研究

What geophysical inversion studied includes the common mathematics physical property of inversion and the constitution and appraisal method of solution in geophysics domain, i.e. using observed physical phenomenon from the earth surface to infer space changing and physical property structure of medium within the earth. Seismic inversion is a branch of geophysical inversion. The basic purpose of seismic inversion is to utilizing seismic wave propagating law in the medium underground to infer stratum structure and space distribution of physical property according to data acquisition, processing and interpretation, and then offer the vital foundation for exploratory development. Poststack inversion is convenient and swift, its acoustic impedance inversion product can reflect reservoir interior changing rule to a certain degree, but poststack data lack abundant amplitude and travel time information included in prestack data because of multiple superimpose and weaken the sensitiveness which reflecting reservoir property. Compared with poststack seismic inversion, prestack seismic inversion has better fidelity and more adequate information. Prestack seismic inversion, including waveform inversion, not only suitable for thin strata physical property inversion, it can also inverse reservoir oil-bearing ability. Prestack seismic inversion and prestack elastic impedance inversion maintain avo information, sufficiently applying seismic gathering data with different incident angle, partial angle stack, gradient and intercept seismic data cube. Prestack inversion and poststack inversion technology were studied in this dissertation. A joint inversion method which synthesize prestack elastic wave waveform inversion, prestack elastic impedance inversion and poststack inversion was proposed by making fully use of prestack inversion multiple information and relatively fast and steady characteristic of poststack inversion. Using the proposed method to extract rock physics attribute cube with clear physical significance and reflecting reservoir characterization, such as P-wave and S-wave impedance, P-wave and S-wave velocity, velocity ratio, density, Poisson ratio and Lame’s constant. Regarding loose sand reservoir in lower member of Minghuazhen formation, 32-6 south districts in Qinhuangdao,as the research object, be aimed at the different between shallow layer loose sand and deep layer tight sand, first of all, acquire physical property parameters suitable for this kind of heavy oil pool according to experimental study, establishing initial pressure and shear wave relational model; Afterwards, performing prestack elastic wave forward and inversion research, summarizing rules under the guidance of theoretical research and numerical simulation, performing elastic impedance inversion, calculating rock physics attributes; Finally, predicting sand body distribution according to rock physics parameters, and predicting favorable oil area combine well-logging materials and made good results.

海上时移地震综合解释技术研究与应用

The practical application and development of the time-lapse seismic reservoir monitor technology has indicated which has already become one of most important development technologies in seeking the surplus oil distribution and improving the reservoir recovering. The paper, first obtained the rock physics experiment analysis data according to the Bohai Sea loose sandstone in-situ measure technical, and determined the feasibility research of the S oil-field on the time-lapse seismic reservoir monitoring combining with the time-lapse numeric simulation technology, which was used to analyze the time-lapse seismic respond raw of the reservoir parameters change and pointed out the attentive problems during the real time-lapse seismic processing and interpretation. Next, simply introduced the technical link and the effect of the time-lapse mutual constrained fidelity and match processing aiming at the local complex gathering condition, geological condition, development engineering condition. Third, introduced the time-lapse integrated interpretation and the technical system with the innovative key technology that includes the time-lapse difference explanation technology, the time-lapse seismic multi-attributes integrated interpretation technology, and the time-lapse constrained reservoir parameters inversion technology, and so on. Using the time-lapse difference direct explanation technology, directly obtained the surplus oil macroscopic distribution through the difference seismic data; Using the presenting 8 big principles of the sublayer isochronisms comparison, carried on the time-lapse integrated interpretation analysis on the fine sublayer comparison and the thin oil-layer(group) contrast and the oil layer (group); The paper putted up the research, contrast, applications of the multi-sides sensitive attribute analysis and the RBF nerve network on the nearest study algorithm, and predicted the reservoir parameters and the surplus oil distribution with them; Combining with innovative researches and the time-lapse seismic constrained reservoir parameters inversion technology realized the good combination of the seismic and the reservoir engineering. Fourth, under fully analyzing the geology condition, the reservoir condition, the exploit dynamic data, and the seismic data of the S oil-field, and analyzing the time-lapse difference factors with reservoir dynamic exploit data, calibrated the oil-gas saturation change, the pressure change, the water saturation change, and determined the rationality of the time-lapse seismic difference, and finally obtained the surplus oil distribution, the water flood characteristic understanding, reservoir degasification, and pressure drop raw, and so on, which had been used in the well pattern tightening plan proof of the S oil-field development adjustment plan. Finally, the paper summarized the knowledge and understanding of the marine time-lapse seismic integrated interpretation, also had pointed out the further need researched question.

叠前地震数据规则化、重建及噪音压制

I address of reconstruction of spatial irregular sampling seismic data to regular grids. Spatial irregular sampling data impairs results of prestack migration, multiple attenuations, spectra estimation. Prestack 5-D volumes are often divided into sub-sections for further processing. Shot gathers are easy to obtain from irregular sampling volumes. My strategy for reconstruction is as follows: I resort irregular sampling gathers into a form of easy to bin and perform bin regularization, then utilize F-K inversion to reconstruct seismic data. In consideration of poor ability of F-K regularization to fill in large gaps, I sort regular sampling gathers to CMP and proposed high-resolution parabolic Radon transform to interpolate data and extrapolate offsets. To strong interfering noise--multiples, I use hybrid-domain high-resolution parabolic Radon transform to attenuate it. F-K regularization demand ultimately for lower computing costs. I proposed several methods to further improve efficiency of F-K inversion: first I introduce 1D and 2D NFFT algorithm for a rapid calculation of DFT operators; then develop fast 1D and 2D CG method to solve least-square equations, and utilize preconditioner to accelerate convergence of CG iterations; what’s more, I use Delaunay triangulation for weight calculation and use bandlimit frequency and varying bandwidth technique for competitive computation. Numerical 2D and 3D examples are offered to verify reasonable results and more efficiency. F-K regularization has poor ability to fill in large gaps, so I rearrange data as CMP gathers and develop hybrid-domain high-resolution parabolic Radon transforms which be used ether to interpolate null traces and extrapolate near and far offsets or suppress a strong interfere noise: multiples. I use it to attenuate multiples to verify performances of our algorithm and proposed routines for industrial application. Numerical examples and field data examples show a nice performance of our method.

Pg震相走时的各向异性影响及各向异性结构成像

As the first arrival of seismic phase in deep seismic sounding, Pg is the important data for studying the attributes of the sedimentary layers and the shape of crystalline basement because of its high intensity and reliable detection. Conventionally, the sedimentary cover is expressed as isotropic, linear increasing model in the interpretation of Pg event. Actually, the sedimentary medium should be anisotropic as preferred cracks or fractures and thin layers are common features in the upper crust, so the interpretation of Pg event needs to be taken account of seismic velocity anisotropy. Traveltime calculation is the base of data processing and interpretation. Here, we only study the type of elliptical anisotropy for the poor quality and insufficiency of DSS data. In this thesis, we first investigate the meaning of elliptical anisotropy in the study of crustal structure and attribute, then derive Pg event’s traveltime-offset relationship by assuming a linear increasing velocity model with elliptical anisotropy and present the invert scheme from Pg traveltime-offset dataset to seismic velocity and its anisotropy of shallow crustal structure. We compare the Pg traveltime calculated by our analytic formula with numerical calculating method to test the accuracy. To get the lateral variation of elliptical anisotropy along the profiling, a tomography inversion method with the derived formula is presented, where the profile is divided into rectangles. Anisotropic imaging of crustal structure and attribute is efficient method for crust study. The imaging result can help us interprete the seismic data and discover the attribute of the rock to analyze the interaction between layers. Traveltime calculation is the base of image. Base on the ray tracing equations, the paper present a realization of three dimension of layer model with arbitrary anisotropic type and an example of Pg traveltime calculation in arbitrary anisotropic type is presented. The traveltime calculation method is complex and it only adapts to nonlinear inversion. Perturbation method of travel-time calculation in anisotropy is the linearization approach. It establishes the direct relation between seismic parameters and travetime and it is fit for inversion in anisotropic structural imaging. The thesis presents a P-wave imaging method of layer media for TTI. Southeastern China is an important part of the tectonic framework concerning the continental margin of eastern China and is commonly assumed to comprise the Yangtze block and the Cathaysia block, the two major tectonic units in the region. It’s a typical geological and geophysical zone. In this part, we fit the traveltime of Pg phase by the raytracing numerical method. But the method is not suitable here because the inefficiency of numerical method and the method itself. By the analytic method, we fit the Pg and Sg and get the lateral variation of elliptical anisotropy and then discuss its implication. The northeastern margin of Qinghai-Tibetan plateau is typical because it is the joint area of Eurasian plate and Indian plate and many strong earthquakes have occurred there in recent years．We use the Pg data to get elliptical anisotropic variation and discuss the possible meaning.

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

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.

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

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.