石油二次运移机理物理模拟研究

Oil and gas migration is very important for theoretical hydrocarbon geology study and exploration practice, but related research is weak. Physical simulation is a main method to study oil migration. Systematic experiments were done to quantitatively describe the migration patterns, path characters and oil saturation by adjusting the possible dynamic factors respectively. The following conclusions were drawn. 1. Darcy velocity and pore throat diameter were calculated according to seepage cross-sectional area and glass beads arrangement. With such normalized Darcy velocity and pore throat diameter, the date from one and two dimensional experiments can be reasonably drawn in two phase diagrams. It is found that the migration pattern can be identified using only one dimensionless number L which is defined as the ration of capillary number and Bond number. 2. Oil saturated in the pores between glass beads was used as calibration and oil saturation in the path was measured by magnetic resonance imaging. The results show that oil saturation in the center of migration path can reach 100%, is higher than oil saturation in the edge of migration path. 3. Percolation backbone during secondary oil migration was identified experimentally using Hele-Shaw cell. The backbone formed mainly because of the spatial variation of the cluster conductivity caused by oil saturation heterogeneity, main resistant force change, and path shrinkage and snap-off. Percolation backbone improves hydrocarbon migration efficiency and is a favorable factor for reservoir forming. 4. In the three dimensional filling models, the thickness of the secondary migration path is mall. It is only 2.5cm even for the piston pattern. Inclination of the model is the main influencing factor of the secondary path width.

酒泉盆地酒东坳陷下白垩统储层特征及有利勘探方向研究

Jiuquan basin, located in the middle of the Hexizoulang, is one of the major important Mesozoic、cenozoic oil-gas bearing basin in the west of China. Jiuquan basin is composed of Jiuxi depression、Huahai-jinta depression and Jiudong depression. Basement of Jiudong depression is Silurian shallow metamorphic rock. Ying-er sag , focus of this study, is the biggest sag in Jiudong depression and the targeting study object is cretacic strata. Structure evolution and geological background were carefully studied in this research. A series of methods were applied to this research: values of oxygen and carbon isotope and trace elements analysis were used to recover salinity of the palae-lake water of the sag. The evolution and distribution of sedimentary faces were carefully studied. Also, various analysis and tests were made to study the diagenesis of the reservoir sandstones、porosity evolution and porosity distribution. All the studies indicate that sedimentary faces are main macroscopic factor controlling the reservoir quality; Compaction is the main factor destroying reservoir property. Carbonate cements greatly preserved the porosity in eodiagenesis because it had prevented significant early mechanical compaction and its dissolution in the late diagenesis generated secondary pores. Carbonate cements in the late diagenesis occluded primary porosity and played a negative role in the porosity preservation. Source of the carbonate cements were also preliminarily discussed. Feldspar grains and lithic fragments were dissolved by acid fluid and formed a great amount of secondary pores and developed the reservoir quality. Also, sedimentary-diagenesis zones were identified. On basis of these studies, Reservoir forming factors were studied. Keywords: Jiudong depression, sedimentary faces，reservoirs diagenesis reservoir evaluation，secondary pores

乍得西南部中新生代火山作用与盆地演化

Located in the Central and West African, Chad, which is not well geological explored, is characterized by Mesozoic- Cenozoic intra-continental rift basins. The boreholes exposed that, during Mesozoic-Cenozoic times, volcanic activities were intense in these basins, but study on volcanic rocks is very weak, especially on those embedded in rift basins, and so far systematic and detailed work has still no carried out. Based on the project of China National Oil and Gas Exploration and Development Corporation, “The analysis of reservoir condition and the evaluation of exploration targets of seven basins in block H in Chad”, and the cooperative project between Institute of Geology and Geophysics, CAS and CNPC International (Chad) Co. Ltd., “Chronology and geochemistry studies on Mesozoic-Cenozoic volcanic rocks from southwestern Chad Basins”, systematic geochronology, geochemistry and Sr-Nd-Pb isotopic geochemistry studies on volcanic rocks from southwestern Chad basins have been done in the thesis for the first time. Detailed geochronological study using whole-rock K-Ar and Ar-Ar methods shows the mainly eruption ages of these volcanic rocks are Late Cretaceous- Paleogene. Volcanic rocks in the well Nere-1 and Figuier-1 from Doba basin are products of the Late Cretaceous which majority of the K-Ar (Ar-Ar) ages fall in the interval 95-75 Ma, whereas volcanic rocks in the well Ronier-1 from Bongor Basin and the Well Acacia-1 from Lake Chad Basin formed in the Paleogene which the ages concentrated in 66-52Ma. Two main periods of volcanic activity can be recognized in the study area, namely, the Late Cretaceous period and the Paleogene period. Volcanic activities have a general trend of south to north migration, but this may be only a local expression, and farther future studies should be carried on. Petrology study exhibits these volcanic rocks from southwestern Chad basins are mainly tholeiitic basalt. Major- and trace elements as well as Sr-Nd-Pb isotopic geochemistry studies show that the late Cretaceous and the Paleogene basalts have a definitely genetic relationship, and magmas which the basalts in southwestern Chad basins derived from were produced by fractional crystallization of olivine and clinopyroxene and had not do suffered from crustal contamination. These basalts are prominently enriched light rare earth elements (LREE), large-ion lithophile elements (LILE) and high field strength elements (HFSE) and depleted compatible elements. They have positive Ba, Pb, Sr, Nb, Ta, Zr, Hf anomalies and negative Th, U, P,Y anomalies. It is possible that the basalts from southwestern Chad basins mainly formed by mixing of depleted mantle (DM) and enriched mantle (EMⅡ) sources. The late Cretaceous basalts have higher (87Sr/86Sr)i ratios than the Paleogene basalts’, whereas have lower (143Nd/144Nd)i ratios than the latter, showing a significant temporal evolution. The mantle sources of the Late Cretaceous basalts may have more enriched mantle(EMⅡ) compositions, whereas those of the Paleogene basalts are relatively more asthenospheric mantle (DM) components. The mantle components with temporal change observed in basalts from Chad basins were probably correlated with the asthenospheric mantle upwelling and lithospheric thinning in Central and Western Africa since Mesozoic. Mesozoic- Cenozoic Volcanism in Chad basins probably is a product of intra- plate extensional stress regime, corresponded to the tectonic setting of the whole West and Central African during Cretaceous. Volcanism is closely correlated with rifting. As time passed from early period to late, the basaltic magma of Chad basins, characterized with shallower genetic depth, higher density and smaller viscosity, probably indicates the gradual strengthening evolution of the rifting. In the initial rife stage, volcanic activities are absent in the study area. Volcanic activities are basiccally corresponded with the strong extensional period of Chad basins, and the eruption of basalts was slightly lagged behind the extensional period. In the post-rift stage (30-0Ma), these basins shifted to the thermal sag phase, volcanic activities in the study area significantly decreased and then terminated.

井约束外推反演可信度与井震相关性关系研究

With the deepening development of oil-gas exploration and the sharp rise in costs, modern seismic techniques had been progressed rapidly. The Seismic Inversion Technique extracts seismic attribute from the seismic reflection data, inverses the underground distribution of wave impedance or speed, estimates reservoir parameters, makes some reservoir prediction and oil reservoir description as a key technology of Seismic exploration, which provides a reliable basic material for oil-gas exploration. Well-driven SI is essentially an seismic-logging joint inversion. The low, high-frequency information comes from the logging information, while the structural characteristics and medium frequency band depend on the seismic data. Inversion results mainly depend on the quality of raw data, the rationality of the process, the relativity of synthetic and seismic data, etc. This paper mainly research on how the log-to-seismic correlation have affected the well-driven seismic inversion precision. Synthetic, the comparison between middle –frequency borehole impedance and relative seismic impedance and well-attribute crossplots have been taken into account the log-to-seismic correlation. The results verify that the better log-to-seismic correlation, the more reliable the seismic inversion result, through the analysis of three real working area (Qikou Sag, Qiongdongnan basin, Sulige gas field).

裂缝型储层检测方法——P波方位各向异性技术

Nowadays, the exploration of fractured reservoir plays a vital role in the further development of petroleum industry through out the world. Fractured hydrocarbon reservoirs are widely distributed in China. Usually, S-wave technique prevails, but it also has its disadvantage, prohibitive expense in S-wave data acquisition and processing. So directly utilizing P-wave data to detect fractures, comes to our mind. We briefly introduce theoretical model (HTI) for fractured reservoir. Then study Ruger’s reflectivity method to recognize reflection P-wave reflection coefficient of the top and bottom interface of HTI layer respectively, and its azimuth anisotropy character. Base on that study, we gives a review and comparison of two seismic exploration technologies for fractures available in the industry-- P-wave AVO and AVA. They has shown great potential for application to the oil and gas prediction of fractured reservoir and the reservoir fine description．Every technique has its disadvantage, AVO limited to small reflection angle; and AVA just offering relatively results. So that, We can draw a conclusion that a better way to any particular field is using synthesis of multiple data sources including core、outcrop、well-test、image logs、3D VSPs, generally to improve the accuracy.

基于分形和统计理论的地层非均质研究

Basin-scale heterogeneity contains information about the traces of the past sedimentary cycle and tectonic process, and has been a major concern to geophysicists because of its importance in resource exploration and development. In this paper, the sonic data of 30 wells of Sulige field are used to inverse the power-law spectra slope and correlation length which are measures of the heterogeneity of the velocity of the log using fractal and statistic correlation methods. By taking the heterogeneity parameters of different wells interpolated, we get power law spectra slope and correlation length contours reflecting the stratum heterogeneity. Then using correlation and gradient, we inverse the transverse heterogeneity of Sulige field. Reservior-scale heterogeneity influnce the distribution of remaining oil and hydrocarbon accumulation. Using wavelet modulus maximum method to divide the sedimentary cycle using Gr data, therefore we can calculate the heterogeneity parameter in each layer of each log. Then we get the heterogeneity distribution of each layer of Sulige field. Finally, we analyze the relation between the signal sigularity and the strata heterogeneity, and get two different sigularity profiles in different areas.

波动方程偏移速度分析

On the subject of oil and gas exploration, migration is an efficacious technique for imagining structures underground. Wave-equation migration (WEM) dominates over other migration methods in accuracy, despite of higher computational cost. However, the advantages of WEM will emerge as the progress of computer technology. WEM is sensitive to velocity model more than others. Small velocity perturbations result in grate divergence in the image pad. Currently, Kirrchhoff method is still very popular in the exploration industry for the reason of difficult to provide precise velocity model. It is very urgent to figure out a way to migration velocity modeling. This dissertation is mainly devoted to migration velocity analysis method for WEM: 1. In this dissertation, we cataloged wave equation prestack depth migration. The concept of migration is introduced. Then, the analysis is applied to different kinds of extrapolate operator to demonstrate their accuracy and applicability. We derived the DSR and SSR migration method and apply both to 2D model. 2. The output of prestack WEM is in form of common image gathers (CIGs). Angle domain common image gathers (ADCIGs) gained by wave equation are proved to be free of artifacts. They are also the most potential candidates for migration velocity analysis. We discussed how to get ADCIGs by DSR and SSR, and obtained ADCIGs before and after imaging separately. The quality of post stack image is affected by CIGs, only the focused or flattened CIGs generate the correct image. Based on wave equation migration, image could be enhanced by special measures. In this dissertation we use both prestack depth residual migration and time shift imaging condition to improve the image quality. 3. Inaccurate velocities lead to errors of imaging depth and curvature of coherent events in CIGs. The ultimate goal of migration velocity analysis (MVA) is to focus scattered event to correct depth and flatten curving event by updating velocities. The kinematic figures are implicitly presented by focus depth aberration and kinetic figure by amplitude. The initial model of Wave-equation migration velocity analysis (WEMVA) is the output of RMO velocity analysis. For integrity of MVA, we review RMO method in this dissertation. The dissertation discusses the general ideal of RMO velocity analysis for flat and dipping events and the corresponding velocity update formula. Migration velocity analysis is a very time consuming work. Respect to computational convenience, we discus how RMO works for synthetic source record migration. In some extremely situation, RMO method fails. Especially in the areas of poorly illuminated or steep structure, it is very difficult to obtain enough angle information for RMO. WEMVA based on wave extrapolate theory, which successfully overcome the drawback of ray based methods. WEMVA inverses residual velocities with residual images. Based on migration regression, we studied the linearized scattering operator and linearized residual image. The key to WEMVA is the linearized residual image. Residual image obtained by Prestack residual migration, which based on DSR is very inefficient. In this dissertation, we proposed obtaining residual migration by time shift image condition, so that, WEMVA could be implemented by SSR. It evidently reduce the computational cost for this method.

纵波方位属性裂缝检测技术综合研究

Fractured oil and gas reservoir is an important type of oil and gas reservoir, which is taking a growing part of current oil and gas production in the whole world. Thus these technologies targeted at exploration of fractured oil and gas reservoirs are drawing vast attentions. It is difficult to accurately predict the fracture development orientation and intensity in oil and gas exploration. Focused on this problem, this paper systematically conducted series study of seismic data processing and P-wave attributes fracture detection based on the structure of ZX buried mountain, and obtained good results. This paper firstly stimulated the propagation of P-wave in weak anisotropic media caused by vertical aligned cracks, and analyzed the rule of P-wave attributes’ variation associated with observed azimuth, such as travel-time, amplitude and AVO gradient and so on, and quantitatively described the sensitive degree of these attributes to anisotropy of fracture medium. In order to further study the sensitive degree of these attributes to anisotropy of fractures, meanwhile, this paper stimulated P-wave propagation through different types and different intensity anisotropic medium respectively and summarized the rule of these attributes’ variation associated with observed azimuth in different anisotropic medium. The results of these studies provided reliable references for predicting orientation, extensity and size of actual complicated cracked medium by P-wave azimuth attributes responses. In the paper, amounts of seismic data processing methods are used to keep and recover all kinds of attributes applied for fracture detection, which guarantee the high accurate of these attributes, thus then improve the accurate of fracture detection. During seismic data processing, the paper adopted the three dimensional F-Kx-Ky field cone filter technique to attenuate ground roll waves and multiple waves, then enhances the S/N ratio of pre-stack seismic data; comprehensively applying geometrical spread compensation, surface consistent amplitude compensation, residual amplitude compensation to recover amplitude; common azimuth processing method effectively preserves the azimuthal characteristics of P-wave attributes; the technique of bend ray adaptive aperture pre-stack time migration insures to obtain the best image in each azimuth. Application of these processing methods guaranteed these attributes’ accuracy, and then improved the accuracy of fracture detection. After comparing and analyzing a variety of attributes, relative wave impedance (relative amplitude) attribute is selected to inverse the orientation of fracture medium; attenuation gradient and corresponding frequency of 85% energy are selected to inverse the intensity of fracture medium; then obtained the fracture distribution characteristics of lower Paleozoic and Precambrian in ZX ancient buried mountains. The results are good accord with the characteristics of faults system and well information in this area.

全局优化的Fourier有限差分偏移方法研究

Seismic exploration is the main method of seeking oil and gas. 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 seismic exploration. Fourier finite-difference (FFD) method is one of the most valuable methods in complex structure exploration, which has obtained good effect. However, in complex media with wider angles, the effect of FFD method is not satisfactory. Based on the FFD operator, we extend the two coefficients to be optimized to four coefficients, then optimize them globally using simulated annealing algorithm. Our optimization method select the solution of one-way wave equation as the objective function. Except the velocity contrast, we consider the effects of both frequency and depth interval. The proposed method can improve the angle of FFD method without additional computation time, which can reach 75° in complex media with large lateral velocity contrasts and wider propagation angles. In this thesis, combinating the FFD method and alternative-direction-implicit plus interpolation(ADIPI） method, we obtain 3D FFD with higher accuracy. On the premise of keeping the efficiency of the FFD method, this method not only removes the azimuthal anisotropy but also optimizes the FFD mehod, which is helpful to 3D seismic exploration. We use the multi-parameter global optimization method to optimize the high order term of FFD method. Using lower-order equation to obtain the approximation effect of higher-order equation, not only decreases the computational cost result from higher-order term, but also obviously improves the accuracy of FFD method. We compare the FFD, SAFFD（multi-parameter simulated annealing globally optimized FFD）, PFFD, phase-shift method(PS）, globally optimized FFD （GOFFD）, and higher-order term optimized FFD method. The theoretical analyses and the impulse responses demonstrate that higher-order term optimized FFD method significantly extends the accurate propagation angle of the FFD method, which is useful to complex media with wider propagation angles.

南海张裂大陆边缘数值模拟研究

The South China Sea (SCS) is one of the largest marginal seas in the western Pacific, which is located at the junction of Eurasian plate, Pacific plate and Indian-Australian plate. It was formed by continent breakup and sea-floor spreading in Cenozoic. The complicated interaction among the three major plates made tectonic movement complex and geological phenomena very rich in this area. The SCS is an ideal place to study the formation and evolution of rifted continental margin and sea-floor spreading since it is old enough to have experienced the major stages of the basin evolution but still young enough to have preserved its original nature. As the demand for energy grows day by day in our country, the deep water region of the northern continental margin in the SCS has become a focus of oil and gas exploration because of its huge hydrocarbon potential. Therefore, to study the rifted continental margin of the SCS not only can improve our understanding of the formation and evolution processes of rifted continental margin, but also can provide theoretical support for hydrocarbon exploration in rifted continental margin. This dissertation mainly includes five topics as follows: (1) Various classic lithosphere stretching models are reviewed, and the continuous non-uniform stretching model is modified to make it suitable for the case where the extension of lithopheric mantle exceeds that of the crust. Then simple/pure shear flexural cantilever model is applied to model the basement geometries of SO49-18 profile in the northern continental margin of the SCS. By fitting the basements obtained by using 2DMove software with modeling results, it is found that the reasonable effective elastic thickness is less than 5km in this region. According to this result, it is assumed that there is weak lower crust in the northern continental margin in the SCS. (2) We research on the methods for stretching factor estimation based on various lithosphere stretching models, and apply the method based on multiple finite rifting model to estimate the stretching factors of several wells and profiles in the northern continental margin of the SCS. (3) We improve one-dimension strain rate inversion method with conjugate gradient method, and apply it to invert the strain rate of several wells in the northern continental margin of the SCS. Two-dimension strain rate forward modeling is carried out, and the modeling results show that effective elastic thickness is a key parameter to control basin’s geometry. (4) We simulate divergent upwelling mantle flow model using finite difference method, and apply this newly developed model to examine the formation mechanism of the northwest and central sub-basin in the SCS. (5) We inverse plate thickness and basal temperature of oceanic lithosphere using sea-floor ages and bathymetries of the North Pacific and the North Atlantic based on varied-parameters plate model, in which the heat conductivity, heat capacity and coefficient of thermal expansion depend on temperature or depth. A new empirical formula is put forward based the inversed parameters, which depicts the relation among sea-floor age, bathymetry and heat flow. Then various similar empirical formulae, including the newly developed one, are applied to examine the sea-floor spread issue in the SCS based on the heat flow and bathymetry data of the abyssal sub-basin.

准噶尔盆地热历史与油气

Based on the temperature data from 196 wells and thermal conductivity measurements of 90 rock samples, altogether 35 heat flow data are obtained. The results show that the Junggar basin is a relatively "cold basin" at present. The thermal gradients vary between 11.6 and 26.5 ℃/km, and the thermal conductivity change from 0.17 to 3.6 W/mK. Heat flow ranges from 23.4 to 53.7 mW/m~2 with a mean of 42.3 ± 7.7 mW/m~2. The heat flow pattern shows that heat flow is higher on the uplifts and lower on the depressions. The overall low present-day heat flow in the Junggar Basin reflects its stable cratonic basement and Cenozoic tectonothermal evolution characterized by lithospheric thickening, thrust and fault at shallow crust as well as consequently quick subsidence during the Late Cenozoic. The study of the basin thermal history, which is one of the important content of the basin analysis, reveals not only the process of the basin's tectonothermal evolution, but also the thermal evolution of the source rocks based on the hydrocarbon generation models. The latter is very helpful for petroleum exploration. The thermal history of the Junggar basin has been reconstructed through the heat flow based method using the VR and Fission track data. The thermal evolutions of main source rocks (Permian and Jurassic) and the formations of the Permian and the Jurassic petroleum systems as well as the influences of thermal fields to petroleum system also have been discussed in this paper. Thermal history reconstruction derived from vitrinite reflectance data indicates that the Paleozoic formations experienced their maximum paleotemperature during Permian to Triassic with the higher paleoheat flow of around 70-85 mW/m~2 and the basin cooled down to the present low heat flow. The thermal evolution put a quite important effect on the formation and evolution of the petroleum system. The Jurassic petroleum system in the Junggar basin is quite limited in space and the source rocks of Middle-Lower Jurassic entered oli-window only along the foreland region of the North Tianshan belt, where the Jurassic is buried to the depth of 5-7 km. By contrast, the Middle-Lower Permian source rocks have initiated oil and gas generation in latter Permian to Triassic, and the major petroleum systems, like Mahu-West Pen 1 Well, was formed prior to Triassic when later Paleozoic formation reached the maximum paleotemperature.

复杂碎屑岩储层动、静态参数计算理论与方法研究

Reservoir characterization and reservoir modeling are two key techniques in petroleum exploration and development. They all are based on the reliable static and dynamic properties of the reservoirs, especially the static and dynamic properties of the reservoirs at each borehole. Without the static and dynamic properties of the reservoir, reservoir characterization and reservoir modeling will pass into nothingness. In fact, the static and dynamic properties of the reservoir are needed in every domain and stage of petroleum exploration & development Today, petroleum industry has reached a stage worldwide that most of the simple & large massive reservoirs have been well explored and developed. As a result, oil companies are paying more and more attention to the exploration & development of the complex & middle to small clastic reservoirs (such as low resistivity sandstone reservoirs, low or no resistivity contrast sandstone reservoirs, conglomerate reservoirs, volcanoclastic reservoirs). In the recent years, oil companies inside and outside China are focusing on the exploration and development elastic reservoirs. Most of the theories & methods being applicable for simple clastic reservoirs can not be used in complex clasic reservoirs. Some theories & methods that are not resolved in the case of simple clasic reservoirs become more impossible to be resolved in the case of complex elastic reservoirs. A set of theories & methods being applicable for computing the static and dynamic properties of the complex elastic reservoirs are developed in this paper and they have been put into practice successfully. These theories & methods are developed by integrating multi-subjects such as geology, well logging and reservoir engineering, in which geology is used as direction and modern well logging technology is used as basis and reservoir engineering is used as assistance and computer technology is used as tool. There are three outstanding breakthroughs in this paper: of the low porosity fractured and/or vuggy carbonate/igneous reservoirs too. A set of practical theories and methods of computing the static properties (such as porosity, saturation, lithology and fluid type) & dynamic properties (such as permeability and production rate) of simple clastic reservoirs have been developed with the hard efforts of many petroleum engineers and scientists in the past 70 years. However, only some of the theories & methods being applicable for simple clastic reservoirs can be used in complex clastic reservoirs after little modification because of the complexity of the complex clastic reservoirs. Most of the theories & methods being applicable for simple clastic reservoirs can not be used in complex clasic reservoirs. Some theories & methods that are not resolved in the case of simple clasic reservoirs become more impossible to be resolved in the case of complex clastic reservoirs. A set of theories & methods being applicable for computing the static and dynamic properties of the complex clastic reservoirs are developed in this paper and they have been put into practice successfully. These theories & methods are developed by integrating multi-subjects such as geology, well logging and reservoir engineering, in which geology is used as direction and modern well logging technology is used as basis and reservoir engineering is used as assistance and computer technology is used as tool. There are three outstanding breakthroughs in this paper:

松辽盆地南部岩性油藏成藏模式与分布规律

According to the basic geologic conditions, the paper is directed by the modem oil-gas accumulation theory and petroleum system in which typical oil pools are analyzed and the shape of lithologic trap and geologic factors are pointed out. The process during which oil and gas migrate from source rock to lithologic trap is rebuilt, and the accumulation model of oil pool is set up. With the comprehensive application of seismic geologic and log data and paying attention to the method and technology which is used to distinguish lithologic accumulation. Promising structural-lithofacies zones are got and the distribution rule of various lithologic accumulation is concluded. With making use of the biologic mark compound, different reservoirs are compared. As a result, the oil and gas in HeiDimiao come from Nenjiang Group's source rocks; in SaErTu from QingShenkou Group's and Nenjiang Group's, and in PuTaohua. GaoTaizi and FuYang from QingShankou Group's. According to the development and distribution of effective source rock, oil distribution and the comparison in the south of SongLiao basin, the characteristic of basin structure and reservoir distribution is considered, and then the middle-upper reservoir of SongLiao basin south are divided into two petroleum system and a complex petroleum system. Because of the characteristic of migration and accumulation, two petroleum systems can furtherly be divided into 6-7 sub-petroleum systems,20 sub-petroleum systems in all. As a result of the difference of the migration characteristic, accumulation conditions and the place in the petroleum system, the accumulation degree and accumulation model are different. So three accumulation mechanism and six basic accumulation model of lithologic trap are concluded. The distribution of lithologic pools is highly regular oil and gas around the generation sag distribute on favorable structural-lithofacies zones, the type of lithological pool vary regularly from the core of sandstone block to the upper zone. On the basic of regional structure and sedimentary evolution, main factors which control the form of trap are discovered, and it is the critical factor method which is used to discern the lithologic trap. After lots of exploration, 700km~2 potential trap is distinguished and 18391.86 * 10~4 tons geologic reserves is calculated. Oil-water distribution rule of pinch-out oil pool is put up on plane which is the reservoirs can be divided into four sections. This paper presented the law of distribution of oil and water in updip pinch-out reservoir, that is, hydrocarbon-bearing formation in plane can be divided into four zones: bottom edge water zone, underside oil and water zone, middle pure oil zone and above residual water zone. The site of the first well should be assigned to be middle or above pure oil zone, thus the exploration value of this type of reservoir can be recognized correctly. In accordance with the characteristics of seism and geology of low permeability thin sandstone and mudstone alternation layer, the paper applied a set of reservoir prediction technology, that is: (1)seism multi-parameter model identification; (2) using stratum's absorbing and depleting information to predict reservoir's abnormal hydrocarbon-bearing range. With the analysis of the residual resource potential and the research of two petroleum system and the accumulation model, promising objective zones are predicted scientifically. And main exploration aim is the DaRngZi bore in the west of ChangLin basin, and YingTai-SiFangZi middle-upper assembly in Honggang terrace.

塔里木盆地轮南地区三叠系油气地质模型及油气预测研究

As an important measure to understand oil and gas accumulation during petroleum exploration and development, Petroleum geological model is an integrated system of theories and methods, which includes sedimentology, reservoir geology, structural geology, petroleum geology and other geological theories, and is used to describe or predict the distribution of oil and gas. Progressive exploration and development for oil and gas is commonly used in terrestrial sedimentary basin in China for the oil and gas generation, accumulation and exploitation are very intricate. It is necessary to establish petroleum geological model, adaptive to different periods of progressive exploration and development practice. Meanwhile there is lack of an integrated system of theories and methods of petroleum geological model suitable for different exploration and development stages for oil and gas, because the current different models are intercrossed, which emphasize their different aspects. According to the characteristics of exploration and development for the Triassic oil and gas pool in Lunnan area, Tarim Basin, the Lunnan horst belt was selected as the major study object of this paper. On the basis of the study of petroleum geological model system, the petroleum geological models for different exploration and development stages are established, which could be applied to predict the distribution of oil and gas distribution. The main results are as follows. (1) The generation-accumulation and exploration-development of hydrocarbon are taken as an integrated system during the course of time, so petroleum exploration and development are closely combined. Under the guidance of some philosophical views that the whole world could be understood, the present writer realizes that any one kind of petroleum geological models can be used to predict and guide petroleum exploration and development practice. The writer do not recognize that any one kind of petroleum geological models can be viewed as sole model for guiding the petroleum exploration and development in the world. Based on the differences of extents and details of research work during various stage of exploration and development for oil and gas, the system of classification for petroleum geological models is established, which can be regarded as theoretical basis for progressive petroleum exploration and development. (2) A petroleum geological model was established based on detailed researches on the Triassic stratigraphy, structure, sedimentology and reservoir rocks in the Lunnan area, northern Tarim Basin. Some sub-belt of hydrocarbon accumulation in the Lunnan area are divided and the predominate controlling factors for oil and gas distribution in the Lunnan area are given out. (3) Geological models for Lunnan and Jiefangqudong oil fields were rebuilt by the combinations of seismology and geology, exploration and development, dynamic and static behavior, thus finding out the distribution of potential zones for oil and gas accumulations. Meanwhile Oil and gas accumulations were considered as the important unit in progressive exploration and development, and the classification was made for Lunnan Triassic pools. Petroleum geological model was created through 3D seismic fine interpretation and detailed description of characteristics of reservoir rocks and the distribution of oil and gas, especially for LN3 and LN26 well zones. The possible distribution of Triassic oil traps and their efficiency in the Lunnan area has been forecasted, and quantitative analysis for original oil(water) saturation in oil pools was performed. (4) The concept of oil cell is proposed by the writer for the first time. It represents the relatively oil-rich zones in oil pool, which were formed by the differences of fluid flows during the middle stage of reservoir development. The classification of oil cells is also given out in this paper. After the studies of physical and numerical modeling, the dominant controlling factors for the formation of various oil cells are analyzed. Oil cells are considered as the most important hydrocarbon potential zones after first recovery, which are main object of progressive development adjustment and improvement oil recovery. An example as main target of analysis was made for various oil cells of Triassic reservoir in the LN2 well area. (5) It is important and necessary that the classification of flow unit and the establishment of geological model of flow unit based on analysis of forecast for inter-well reservoir parameters connected with the statistical analysis of reservoir character of horizontal wells. With the help of self-adaptive interpolation and stochastic simulation, the geological model of flow units was built on the basis of division and correlation of flow units, with which the residual oil distribution in TIII reservoir in the LN2 well area after water flooding can be established.

柴达木盆地北缘地区油气成藏过程模拟与含油气系统评价

As an important part of petroleum exploration areas in the west of China, the north part of Qaidam basin is very promising in making great progress for petroleum discovery. But there are still many obstacles to overcome in understanding the process of petroleum formation and evaluation of oil & gas potential because of the complexity of geological evolution in the study area. Based upon the petroleum system theory, the process of petroleum formation is analyzed and the potential of oil & gas is evaluated in different petroleum systems by means of the modeling approach. The geological background for the formation of petroleum systems and the consisting elements of petroleum systems are described in detail. The thickness of strata eroded is estimated by means of vitrinite reflectance modeling, compaction parameter calculating and thickness extrapolating. The buried histories are reconstructed using the transient compaction model, which combines of forward and reverse modeling. The geo-history evolution consists of four stages - sedimentation in different rates with different areas and slow subsidence during Jurassic, uplifting and erosion during Cretaceous, fast subsidence during the early and middle periods of Tertiary, subsidence and uplifting in alternation during the late period of Tertiary and Quaternary. The thermal gradients in the study area are from 2.0 ℃/100m to 2.6 ℃/100m, and the average of heat flow is 50.6 mW/m~2. From the vitrinite reflectance and apatite fission track data, a new approach based up Adaptive Genetic Algorithms for thermal history reconstruction is presented and used to estimate the plaeo-heat flow. The results of modeling show that the heat flow decreased and the basin got cooler from Jurassic to now. Oil generation from kerogens, gas generation from kerogens and gas cracked from oil are modeled by kinetic models. The kinetic parameters are calculated from the data obtained from laboratory experiments. The evolution of source rock maturation is modeled by means of Easy %Ro method. With the reconstruction of geo-histories and thermal histories and hydrocarbon generation, the oil and gas generation intensities for lower and middle Jurassic source rocks in different time are calculated. The results suggest that the source rocks got into maturation during the time of Xiaganchaigou sedimentation. The oil & gas generation centers for lower Jurassic source rocks locate in Yikeyawuru sag, Kunteyi sag and Eboliang area. The centers of generation for middle Jurassic source rocks locate in Saishenteng faulted sag and Yuka faulted sag. With the evidence of bio-markers and isotopes of carbonates, the oil or gas in Lenghusihao, Lenghuwuhao, Nanbaxian and Mahai oilfields is from lower Jurassic source rocks, and the oil or gas in Yuka is from middle Jurassic source rocks. Based up the results of the modeling, the distribution of source rocks and occurrence of oil and gas, there should be two petroleum systems in the study area. The key moments for these two petroleum, J_1-R(!) and J_2-J_3, are at the stages of Xiaganchaigou-Shangyoushashan sedimentation and Xiayoushashan-Shizigou sedimentation. With the kinetic midels for oil generated from kerogen, gas generated from kerogen and oil cracked to gas, the amount of oil and gas generated at different time in the two petroleum systems is calculated. The cumulative amount of oil generated from kerogen, gas generated from kerogen and gas cracked from oil is 409.78 * 10~8t, 360518.40 * 10~8m~3, and 186.50 * 10~8t in J_1-R(!). The amount of oil and gas generated for accumulation is 223.28 * 10~8t and 606692.99 * 10~8m~3 in J_1-R(!). The cumulative amount of oil generated from kerogen, gas generated from kerogen and gas cracked from oil is 29.05 * 10~8t, 23025.29 * 10~8m~3 and 14.42 * 10~8t in J_2-J_3 (!). The amount of oil and gas generated for accumulation is 14.63 * 10~8t and 42055.44 * 10~8m~3 in J_2-J_3 (!). The total oil and gas potential is 9.52 * 10~8t and 1946.25 * 10~8m~3.