塔里木盆地寒武－奥陶系白云岩类型、成因机理及其孔隙性研究

Cambrian-Ordovician dolostones in Tarim Basin are hydrocarbon reservoir rocks of vital importance. Under the guidance of the theories of sedimentology and the sedimentology of carbonate reservoir, based on the first-hand qualitative and quantitative data especially, combined with micro-study, geochemical and reservoir capacity analysis, and precursor research, the origin and reservoir characteristics of the dolostones were discussed. Based on detailed petrographic investigations, four types of dolostone have been recognized, which are, respecitively, mud-silt-sized dolostones, algal laminated dolostones (ALD), prophyritic dolostone, and neomorphic dolostone. Mud-silt-sized dolostones always presents as laminas together with evaporated signatures, its REE patterns and ΣREE are all close to that of the finely crystalline limestone. This kind of dolomite probably experienced relatively low fluid-rock ratio during diagenesis was formed in hypersaline and oxidizing environment and involved fast dolomitization process. It was dolomitized by evaporated seawater in sabkha environment.The main primary fabrics of algal lamination in algal laminated dolomite (ALD) can still be identified and its ΣREE (21.37) is very close to that of algae. This reveals that ALD was dolomitized during early diagenesis and algae possibly played an important role. The ALD was formed under mediation of organic matter and Mg2+ were supplied by magnesium concentrated algal laminites and sea water. Prophyritic dolostones presents mainly as patchy occurrence and yield the lowest δ13C and Z value. Its ΣREE is much less than that of the finely crystalline limestone. These characteristics reveal that the cloudy cores were dolomitized in shallow early diagenetic environments by pore fluids riched in Mg2+. Whereas the clear rims were likely formed in subsequent burial into deeper subsurface environments, and the Mg2+ needed for further dolomitization possibly was supplied by the transformations of clay minerals. Neomorphic dolostones consist of coarse, turbid crystals and exhibits sucrosic and mosaic textures. It has highest Fe2+ contents and average homogeneous temperature (110.2℃). Collectively, these characteristics demonstrate that the neomorphic dolostones was likely formed by recrystallization of pre-existing dolomites during deep burial. The ΣREE of the four types of dolostone distinctly differentiates from each other. However, their REE patterns are all enriched in LREE, depleted in HREE and have Eu negative anomalies. Its ΣREE 13.64 ppm, less than 1/4 of finely crystalline limestone, and ranks the lowest in the 4 types.These characteristics are comparable to those of finely crystalline limestone, and are mainly infuenced by the sea water. These four types of dolostone show similar REE mobility behaviour and no significant fractionation, althouth they have been subjected to evidently different diageneses. Seven main pore types are identified in the dolostones , which are fenestral, moldic, intercrystal, dissolved，breccia, dissolved breccia and stylolite pores. Fenestral pores are primary and the others are secondary. The dissolved pores and intercrystal pores are the most important reservoir spaces and followed by breccias and dissolved breccia pores, and the moldic and fenestral pores are less important. Stylolites can enhance permeability of reservoir rocks in one hand, for the other hand, the capacity of reservoir and permeability are enhanced and then better reservoir rocks can be formed when they are combined with patchy dolostones. The relationship between porosity and the type of dolostones is that the dissolved neomorphic dolostones have the highest porosity of 3.65%, than followed by dissolved Mud-silt-sized dolostones of 3.35%. The mud-silt-sized dolostones without dissolution have the lowest porosity of 0.90%. Moreover, the porosity of prophyritic dolostones and the neomorphic dolostones without dissolution are lower, respectively 1.675% and 1.41%. Although algal laminated dolostones consist of euhedral crystals and riched in intercrystal pores, its porosity just yields 1.20%. The relationship between permeability and the type of dolostones is that that algal laminated dolostones have the highest permeability of 0.462mD and followed by 0.065mD of prophyritic dolostones. Dissolution have no significant influence on the permeability of neomorphic dolostones and this presented by the permeability of dissolved and non-dissolved are very close, respectively 0.043mD and 0.062mD. No matter dissolved or not, mud-silt-sized dolostones are much less permeable. The permeability of non-dissolved and dissolved are 0.051mD and 0.016mD. Collectively, in the study area, neomorphic dolostones can form high quality reservoir.

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

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:

孤东油田馆陶组河流相储层流动单元模型与剩余油分布规律研究

Since C.L. Hearn presented the concept of flow unit in 1984, its generation mechanisms and controlling factors have been studied in many aspects using different methods by researchers. There are some basic methods to do the research, and there are several concepts and classification standards about flow unit. Based on previous achievements and using methodologies from sedimentary geology, geophysics, seismic stratigraphy, and reservoir engineering, the author systemically studies the factors controlling flow unit, puts forward a series of methods for recognition, classification and evaluation of flow unit. The results obtained in this paper have important significance not only for understanding the flow unit, but also for revealing the distribution of remaining oil. As a case, this paper deals with the reservoir rocks in Guantao Group of Gudong Oilfield. Zhanhua Sag, Jiyang Depression in Bohaiwan Basin. Based on the study of stratigraphic, depositional and structural characteristics, the author establishes reservoir geological models, reveals the geological characteristics of oil-bearing reservoir of fluvial facies, points out the factors controlling flow unit and geological parameters for classification of flow unit. and summarizes methods and technologies for flow unit study when geological, well-logging and mathematical methods are used. It is the first attempt in literatures to evaluate reservoir by well-logging data constrained by geological conditions, then a well-logging evaluation model can be built. This kind of model is more precise than ever for calculating physical parameters in flow unit. In a well bore, there are six methods to recognize a flow unit. Among them, the activity function and intra-layer difference methods are the most effective. Along a section, the composition type of flow unit can be located according amplitude and impedance on seismic section. Slice method and other methods are used to distinguish flow unit. In order to reveal the distribution laws of flow unit in space, the author create a new method, named combination and composition of flow unit. Based on microscopic pore structure research, the classification methods of flow unit are developed. There are three types of flow unit in the reservoir of fluvial facies. They have their own lithology, petrophysics and pore structure character. Using judgement method, standard functions are built to determine the class of flow unit of fluvial facies. Combining reservoir engineering methods, the distribution laws of remaining oil in different types, or in different part of a flow unit are studied. It is evident that the remaining oil is controlled by the type of flow unit. The author reveals the relationship between flow unit and remaining oil distribution, builds the flowing models, predicts the variation of reservoir parameters in space, put forward different methods developing remaining oil in different flow unit. Especially, based on the results obtained in this paper, some suggestions for the adjustment of the developing flow units have been applied in Districts No.4 and No.7, and good results have been yielded. So, the results of this paper can guide oil field development. They are useful and significant for developing the remaining oil and enhancing the oil recovery efficiency.

三角洲储层孔隙结构模型与剩余油分布研究——以临南油田沙河街组储层为例

The study of pore structure in reservoir was paid attention to in the early reservoir research and now a systematic research methodology is set up. On the limits of tools and conditions, methodologies and technologies on formation condition and distribution laws of pore structure and the relationship between remaining oil distribution and pore structure are uncertain and some knownage about it is also uncertain. As the development of petroleum industry, the characterization of pore structure and the prediction of remaining oil are the hot spot and difficult point in the research of oil development. The author pays a close attention to this subject and has done much research on it. In a case study in Linnan oilfield Huimin sag Jiyang Depression Bohai Bay basin by using a new method, named varied scale comprehensive modeling of pore structure, the author builds pore structure models for delta reservoir, reveals the remaining oil distribution laws in delta facies, and predicts the distribution of remaining oil in Linnan oilfield. By the application of stratigraphy, sedimentology and structure geology. the author reveals the genetic types of sandbody and its distribution laws, builds the reservoir geological models for delta sandstone reservoir in Shahejie group in Linnan oilfield and points out the geological Factors that control the development of pores and throats. Combining petrology and the reservoir sensitive analysis, the author builds the rock matrix models. It is the first time to state that rocks in different sentimental micro facies have different sensitive .response to fluid pressed into the rocks. Normally. the reservoirs in the delta front have weaker sensitivity to fluid than the reservoirs in delta plain, In same subfacies, the microfacies that have fine grain, such as bank and crevasse splay, have stronger reservoir sensitivity than the microfacies that have coarse grains, such as under-water branched channel and debauch bar. By the application of advanced testing, such as imagine analysis, scan electronic microscope, and morphology method, the author classifies the pore structure and set up the distribution models of pore, throat and pore structure. By the application of advanced theory in well-logging geology, the author finds the relationship between microscope pore structure and macroscopic percolation characteristics, and then builds the well-logging interpretation formulae for calculating pore structure parameters. By using the geostatistics methods, the author reveals the spatial correlative characteristics of pore structure. By application of conditional stochastic simulation methods, the author builds the 3D models of pore structure in delta reservoir. It is the base of predicting remaining oil distribution. By a great deal of experiments and theoretical deduction, The author expounds the laws of percolation flow in different pore structures, and the laws by which the pore structure controls the micro distribution of remaining oil, and then, states the micro mechanism of remaining oil distribution. There are two types of remaining oil. They are by-pass flow caused by micro-fingering and truncation caused by non-piston movement. By new method, the author states the different pore structure has different replacement efficiency, reveals the formation condition and distribution laws of remaining oil. predicts the remaining oil distribution in Linnan oil field, and put forward some idea about how to adjust the oil production. The study yielded good results in the production in Linnan oilfield.

沾化凹陷火成岩油气藏形成条件和分布规律研究

As a kind of special lithologic ones, Igneous rock oil and gas pool is more and more paid attention, and it has different forming condition and distribution from conventional ones, such as various terrane distribution types, serious reservoir anisotropy, complicated hydrocarbon-bearing, so there is not successful experience to follow for exploration and development of this complex subtle oil and gas pool at present. For an example of Igneous oil and gas pool of Luo151 area in Zhanhua seg, Eastern China, this article study the difficult problem, including petrologic nd lithofacies analysis, Origin, invasion age and times of Igneous rock, reservoir anisotropy, Geological Modeling, Igneous reservoir synthesis evaluation. forming condition and distribution are studied synthetically, and an integrated method to predict igneous rock oil and gas pool is formed, which is evaluated by using development data. The Igneous rock is mainly diabase construction in Luo151 area of Zhanhua Sag, and petrologic types include carbonaceous slate, hornfels, and diabases. Based on analyzing synthetically petrologic component, texture and construct, 4 lithofacies zones, such as carbonaceous slate subfacies, hornfels subfacies containing cordierite and grammite, border subfacies and central subfacies, are divided in the diabase and wall rock. By studying on isotopic chronology, terrane configuration and imaging logging data, the diabase intrusion in Zhanhua Sag is formed by tholeiite magma emplacing in Shahejie formation stratum on the rift tension background Lower Tertiary in North China. The diabase intrusion of Luo151 is composed possibly of three periods magma emplacement. There is serious anisotropy in the diabase reservoirs of Luo151 in Zhanhua Sag. Fracture is primary reservoir space, which dominated by tensile fracture in high obliquity, and the fracture zones are mainly developed round joint belt of igneous rock and wall rock and position of terrane thickness changing rapidly. The generation materials of the reservoirs in Luo151 igneous oil pools consist of Intergranular micropore hornfels, condensate blowhole-solution void diabase condensate edge, the edge and center of the condensate seam diabase, of which are divided into horizontal, vertical and reticulated cracks according fracture occurrence. Based on the above research, a conceptual model of igneous rock reservoir is generated, which is vertically divided into 4 belts and horizontally 3 areas. It is built for the first time that classification evaluation pattern of igneous rock reservoir in this area, and 3 key wells are evaluated. The diabase construction is divided into grammite hornfels micropore type and diabase porous-fracture type reservoirs. The heavy mudstone layers in Third Member of Shahejie formation (Es3) provide favorable hydrocarbon source rock and cap formation, diabase and hornfels belts serve as reservoirs, faults and microcracks in the wall rocks as type pathways for oil and gas migration. The time of diabase invasion was about in the later deposition period of Dongying Formation and the middle of that of Guantao Formation, the oil generated from oil source rock of Es3 in the period of the Minghuazhen formation and is earlier more than the period of diabase oil trap and porous space forming. Based on geological and seismic data, the horizon of igneous rocks is demarcated accurately by using VSP and synthetic seismogram, and the shape distribution and continuity of igneous rocks are determined by using cross-hole seismic technology. The reservoir capability is predicted by using logging constraining inversion and neural network technology. An integrated method to predict igneous rock oil and gas pool is formed. The study is appraised by using development data. The result show the reservoir conceptual model can guide the exploration and development of oil pool, and the integrated method yielded marked results in the production.

吐合盆地鲁克沁构造带断块油藏形成机制及油藏预测

Lukeqin arc belt is a compound structure generated by multi-movements and composed of 6 sub-structural zones, which are connected by Huoyanshan Mountain. General characteristics of the arc belt are multi-patterns of structure, multi-phases for petroleum, multi-types of trap and multi-layers for reservoirs. As a part of the eastern Lukeqin arc belt located on the south of Taibei depression, Lukeqin structural zone behaves as a complex faulted-fold zone, in which the formation and distribution of hydrocarbons are controlled by structures. As the dominant source of dynamics for the second migration of hydrocarbon, structure stress field is closely related with the potentials of hydrodynamics. Results derived from the simulations of stress field by finite element method indicate that the northwest tending faults prefer seal to the northeast tending ones. The reason is that the northwest tending faults were squeezed more strongly than the northeast tending ones. Therefor, the northeast tending faults become always the paths for oil to migrate southeastward. Lukeqin structural zone is the main site for oil to concentration because it is surrounded by high stress. Situated on the front of the foreland basin of Turpan-Hami, Lukeqing arc belt is a dam to hold back the southward migrating oil from Shengbei depression. The axis line of Shenquan-Shengnan-Yanmuxi, Lukeqin and Yubei controls the migrating paths and concentrating process of oil and gas. Results derived from stress simulation and structure analyses indicate consistently that both Yubei and Lukeqin structural zones are the favorite areas for oil to migrate. The generally southward paths for oil to migrate out of Taibei depression can be two ways. One of them is from Taibei depression to Yubei structural zone and the other is from Taibei depression to Lukeqin structural zone. By the both ways, oil migrated upward along the faults and southeastward along the structural axis to concentrate in either Permian or Triassic system. The newly ascertained path for oil migration, which is accurately southeastward instead of coarsely southward, indicates the directions for further explorations on the compound Lukeqin block zone. Five kinds of seal models of fault are all found in Lukeqin block zone by studying the seal features of faults occurred in the zone. Having studied the fault seal and their controlling factors by fuzzy set method, the paper deems that the northwest tended faults are better than the northeast tended ones for oil to concentrate. The most important factors to decide the seal extent of faults in this zone are the characteristics of main stress and fluids instead of capillary pressure differences between the two sides of fault and smear mud factors. There exist seal differences not only between the faults of different time but also between the sections within a fault due to the variation of depths, strata and positions. The general distribution rules of reservoirs were dominated by the seal characteristics of a fault during the time reservoirs formed. While the current features of fault seal decide the conservation of reservoirs and heights of oil accumulations. Seal or not of a fault is not absolute because the essential for fault to seal is the distribution of permeability of fault zone. Therefor, the multi cyclical activities of faults create the space-time variation of seal features of the fault. Totally, the seal extent of the faults within the area is not as perfect as to accumulate ordinary crude. Crude oil can only be sealed when it becomes viscous. Process for crude oil to become viscous and viscous happened strongly because of the fault-fold movements. Shallowly burying and even revealing of the objective layers of the reservoirs made the crude oil to be thickened by water washing biologically degradation and oxidation degradation. The northwestward deepening during or after the reservoir formation of the structural zone provided the power for oil to migrate one or more times. The main reason for oil accumulation is the formation of Lukeqin block zone during Xishanyao stage, middle Jurassic Period, Early Yanshanian Movement. While the main reason for reservoir conservation is the placidity of Triassic blocks after the formation of reservoirs. Contrasting to former opinions, it is concluded that the reservoirs in Lukeqin zone, including viscous reservoirs, were formed by one time but not more times. So the author proposes the opinion that the reservoirs of viscous oil were formed by viscous oil migration under the conditions of aptitude sets of fault seals controlled by fluid and other factors. To grope the distribution rules outside Taibei depression and discuss the formation mechanism of Anjurassic reservoirs, it is necessary to study the dominate factors for the formation of reservoirs in Lukeqin structural zone such as structural stress, fault seals and thickening mechanism of crude oil. Also, the necessary studies are the key to break through the Taibei depression and Anjurassic systems. Therefor, they are significant for the future exploration and reserve increasing of hydrocarbon within the Turpan-Hami basin. The paper studied the distribution rules of block reservoirs and forecasted the favorable zones for further exploration in Turpan-Hami basin. Conclusions can be useful for not only the exploration in the area but also the theory consult in the adjacent areas.

岩性油藏储层预测的地震理论与方法研究

At present the main object of the exploration and development (E&D) of oil and gas is not the structural oil-gas pools but the subtle lithological oil-gas reservoir. Since the last 90's, the ratio of this kind of pools in newly-added oil reserves is becoming larger and larger, so is the ratio in the eastern oilfields. The third oil-gas resource evaluation indicates the main exploration object of Jiyang depression is the lithological oil-gas pools in future. However, lack of effective methods that are applied to search for this kind of pool makes E&D difficult and the cost high. In view of the urgent demand of E&D, in this paper we deeply study and analyze the theory and application in which the seismic attributes are used to predict and describe lithological oil-gas reservoirs. The great results are obtained by making full use of abundant physics and reservoir information as well as the remarkable lateral continuity involved in seismic data in combination with well logging, drilling-well and geology. ①Based on a great deal of research and different geological features of Shengli oilfield, the great progresses are made some theories and methods of seismic reservoir prediction and description. Three kinds of extrapolation near well seismic wavelet methods-inverse distance interpolation, phase interpolation and pseudo well reflectivity-are improved; particularly, in sparse well area the method of getting pseudo well reflectivity is given by the application of the wavelet theory. The formulae for seismic attributes and coherent volumes are derived theoretically, and the optimal method of seismic attributes and improved algorithms of picking up coherent data volumes are put forward. The method of making sequence analysis on seismic data is put forward and derived in which the wavelet transform is used to analyze not only qualitatively but also quantitatively seismic characteristics of reservoirs.② According to geologic model and seismic forward simulation, from macro to micro, the method of pre- and post-stack data synthetic analysis and application is put forward using seismic in close combination with geology; particularly, based on making full use of post-stack seismic data, "green food"-pre-stack seismic data is as possible as utilized. ③ In this paper, the formative law and distributing characteristic of lithologic oil-gas pools of the Tertiary in Jiyang depression, the knowledge of geological geophysics and the feasibility of all sorts of seismic methods, and the applied knowledge of seismic data and the geophysical mechanism of oil-gas reservoirs are studied. Therefore a series of perfect seismic technique and software are completed that fit to E&D of different categories of lithologic oil-gas reservoirs. ④ This achievement is different from other new seismic methods that are put forward in the recent years, that is multi-wave multi-component seismic, cross hole seismic, vertical seismic, and time-lapse seismic etc. that need the reacquisition of seismic data to predict and describe the oil-gas reservoir. The method in this paper is based on the conventional 2D/3D seismic data, so the cost falls sharply. ⑤ In recent years this technique that predict and describe lithologic oil-gas reservoirs by seismic information has been applied in E&D of lithologic oil-gas reservoirs on glutenite fans in abrupt slop and turbidite fans in front of abrup slop, slump turbidite fans in front of delta, turbidite fans with channel in low slope and channel sanbody, and a encouraging geologic result has been gained. This achievement indicates that the application of seismic information is one of the most effective ways in solving the present problem of E&D. This technique is significant in the application and popularization, and positive on increasing reserves and raising production as well as stable development in Shengli oilfield. And it will be directive to E&D of some similar reservoirs

吐哈盆地山前带构造岩相带及油气藏预测研究

This paper based on the example of mountain front tectonic belt of Tuha basin - a typical basin in west China, using geological, logging and seismic data, combining the two methods-comprehensive analysis of geological-geophysical data and numerical simulation together, studies the types of tectonic lithofacies and the mechanism of oil/gas accumulation in mountain front tectonic belt, and predicts favorable exploration areas. The tectonic lithofacies in periphery belt of compressional basin is systematically studied for the first time. The distribution and systematically analyzed with geological, logging and seismic data; the bury history of formations in the mountain front are restored with the back-stripping technique; the tectonic depression history of the mountain front is studied according to Airy Balanced Model; the tectonic evolution history is restored with the balanced section technique; according to the space composition of tectonic units and sedimentary systems of different geological time, the tectonic facies belt of the middle and lower Jurassic in the piedmont is divided into 12 types, and the controls of each type's tectonic lithofacies belt on the conditions of oil/gas generation, reservation and sealing are analyzed in depth. 3D numerical simulation and analysis is applied to the tectonic stress field in the hydrocarbon reservoir-forming-period for the first time. Because of the complex evolution of mountain front of Tuha basin, 3D numerical simulation of tectonic stress field in the hydrocarbon reservoir-forming-period helps to study the magnitude of the maximum principal stress, the minimum principal stress and the shear stress, the range and distribution of the principal stress, and controls of the upwards factors on the oil/gas migration and accumulation. Through the study on the oil-controlling fault's evolution, sealing mechanism and sealing history, the coupling of two effects of the fault-passage and sealing screen for oil/gas migration can be defined. Using the basic principle of petroleum system analysis, the paper systematically studies the hydrocarbon reservoir-forming mechanism and the time-space matching of the factors that affect the formation of reservoir, such as the space matching of active oil/gas matching of the active period of fault, the migration period of oil/gas and the formation period of trap. Through comprehensive analysis, the favorable exploration targets are selected in selected in the mountain front.

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

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

火成岩储层裂缝识别、预测理论与方法研究－以济阳坳陷火成岩油藏为例

Since 1970s, igneous reservoirs such as Shang741, Bin674 and Luol51 have been found in Jiyang depression, which are enrichment and heavy-producing. Showing good prospect of exploration and development, igneous reservoirs have been the main part of increasing reserves and production in Shengli oilfield. As fracture igneous reservoir being an extraordinary complex concealed reservoir and showing heavy heterogeneity in spatial distribution, the study of recognition, prediction, formation mechanism and the law of distribution of fracture is essential to develop the reservoir. Guided by multiple discipline theory such as sedimentology, geophysics, mineralogy, petroleum geology, structural geology and reservoir engineering, a set of theories and methods of recognition and prediction of fractured igneous rock reservoir are formed in this paper. Rock data, three-dimensional seismic data, log data, borehole log data, testing data and production data are combined in these methods by the means of computer. Based on the research of igneous rock petrography and reservoir formation mechanism, emphasized on the assessment and forecast of igneous rock reservoir, aimed at establishing a nonhomogeneity quantification model of fractured igneous rock reservoir, the creativity on the fracture recognition, prediction and formation mechanism are achieved. The research result is applied to Jiyang depression, suggestion of exploration and development for fractured igneous rock reservoir is supplied and some great achievement and favourable economic effect are achieved. The main achievements are gained as follows: 1. The main facies models of igneous rock reservoir in JiYang depression are summarized. Based on data and techniques of seism, well log and logging,started from the research of single well rock facies, proceeded by seismic and log facies research, from point to line and line to face, the regional igneous facies models are established. And hypabyssal intrusion allgovite facies model, explosion volcaniclastic rock facies model and overfall basaltic rocks facies model are the main facies models of igneous rock reservoir in JiYang depression. 2. Four nonhomogenous reservoir models of igneous reservoirs are established, which is the base of fracture prediction and recognition. According to characteristics of igneous petrology and spatial types of reservoir, igneous reservoirs of Jiyang depression are divided into four categories: fractured irruptive rock reservoir, fracture-pore thermocontact metamorphic rock and irruptive rock compound reservoir, pore volcanic debris cone reservoir and fracture-pore overfall basaltic rock reservoir. The spatial distribution of each model's reservoir has its features. And reservoirs can be divided into primary ones and secondary ones, whose mechanism of formation and laws of distribution are studied in this paper. 3. Eight geologic factors which dominate igneous reservoirs are presented. The eight geologic factors which dominates igneous reservoirs are igneous facies, epigenetic tectonics deformation, fracture motion, intensity of intrusive effect and adjoining-rock characters, thermo-contact metamorphic rock facies, specific volcano-tectonic position, magmatic cyclicity and epigenetic diagenetic evolution. The interaction of the eight factors forms the four types nonhomogenous reservoir models of igneous reservoirs in Jiyang depression. And igneous facies and fracture motion are the most important and primary factors. 4. Identification patterns of seismic, well log and logging facies of igneous rocks are established. Igneous rocks of Jiyang depression show typical reflecting features on seismic profile. Tabular reflection seismic facies, arc reflection seismic facies and hummocky or mushroom reflection seismic facies are the three main facies. Logging response features of basic basalt and diabase are shown as typical "three low and two high", which means low natural gamma value, low interval transit-time, low neutron porosity, high resistivity and high density. Volcaniclastic rocks show "two high and three low"-high neutron porosity, high interval transit-time, low density, low-resistance and low natural gamma value. Thermo-contact metamorphic rocks surrounding to diabase show "four high and two low" on log data, which is high natural gamma value, high self-potential anomaly, high neutron porosity, high interval transit-time and low density and low-resistance. Based on seismic, well log and logging data, spatial shape of Shang 741 igneous rock is described. 5. The methods of fracture prediction and recognition for fractured igneous reservoir are summarized. Adopting FMI image log and nuclear magnetic resonance log to quantitative analysis of fractured igneous reservoir and according to formation mechanism and shape of fracture, various fractures are recognized, such as high-angle fracture, low-angle fracture, vertical fracture, reticulated fracture, induced fracture, infilling fracture and corrosion vug. Shang 741 intrusive rock reservoir can be divided into pore-vug compound type, pore fracture type, micro-pore and micro-fracture type. Physical properties parameters of the reservoir are computed and single-well fracture model and reservoir parameters model are established. 6. Various comprehensive methods of fracture prediction and recognition for fractured igneous reservoir are put forward. Adopting three-element (igneous facies, fracture motion and rock bending) geologic comprehensive reservoir evaluation technique and deep-shallow unconventional laterolog constrained inversion technique, lateral prediction of fractured reservoir such as Shang 741 is taken and nonhomogeneity quantification models of reservoirs are established.

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

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

塔里木盆地奥陶系不整合面岩溶作用特征与油气储层研究

Based on the principle and methods of carbonate sedimentology and reservoir geology, and guided by the theories of carbonate reservoir geology, the palaeokarst of Ordovician carbonate rocks in Tarim Basin has been comprehensively studied with multiple methods from different branches of geology. It is indicated that the features and distribution of palaeokarstification have developed in Ordovician carbonates. The controlling of karstification to Ordovician carbonate reservoirs has been discussed. Regional distribution of carbonate reservoirs controlled by karstification has been predicted within this basin. The main consents and conclusions of the this dissertation is as follows: Nine key indicators to the recognition of palaeokarst are proposed in terms of careful observation upon the well cores, lithological and geochemical analyses, and drilling and logging responses to the karst caves and fractures. The time and environment of cave filling are documented from careful research of lithofacies, mineralogy, and geochemistry of the physical and chemical fillings within karst caves. The caves in Ordovician carbonates were filled in Early Carboniferous in Lunnan area. The muddy filling in upper caves was deposited under subaerial fresh-water setting, while the muddy filling in lower caves was formed in the mixed water body of fresh-water and dominated sea water. Although most chemical fillings are suggested being precipated in the burial diagenetic environment after karstification but mineralogic and geochemical characteristics of some chemical fillings indicates they formed in meteoric environment during the karstification. It is obvious that the palaeokarst has been zoned in vertical profile. It can be divided into four units from top to bottom: surface karst, vadose karst, phreatic and tranquil flow zones. Between two types of limestone karst and dolostone karst are firstly differentiated in Tarim Basin, based on the comparison of features of each karst zone in limestone and dolostone regions. In Tabei area, the lowest depth of karstification is approximately 300 m below the Upper Ordovician unconformity interface, while the bottom depth of karstification in Tazhong area ranges commonly from 300 to 400 m, in rare cases may be up to 750 m below the upper Ordovician unconformity interface. In Lunnan and Tazhong areas, the palaeokarst morphology and the surface hydrosystem are firstly reconstructed based on the top of carboniferous "Shuangfeng limestone bed (Double-Peaks limestone)" as basal. According to the palaeomorphologic feature, karst topography can be divided into three units: karst upland, karst slope, and karst valley. Vadose zone was well developed in karst upland, and it can be found in a quite depth. Both vadose and phreatic zones were well developed in karst slope and upstream valley. In downstream valley, the karstification is not strong, the vadose and phreatic zones are thin in thickness. In Tazhong and Yingmaili areas, karstification is also developed in relict carbonate palaeo-hills which existed as isolated blocks admits clastic strata.

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

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

三维油藏数值模拟正反演算法研究

The dynamic prediction of complex reservoir development is one of the important research contents of dynamic analysis of oil and gas development. With the increase development of time, the permeabilities and porosities of reservoirs and the permeability of block reservoir at its boundaries are dynamically changing. How to track the dynamic change of permeability and porosity and make certain the permeability of block reservoir at its boundary is an important practical problem. To study developing dynamic prediction of complex reservoir, the key problem of research of dynamic prediction of complex reservoir development is realizing inversion of permeability and porosity. To realize the inversion, first of all, the fast forward and inverse method of 3-dimension reservoir simulation must be studied. Although the inversion has been widely applied to exploration and logging, it has not been applied to3-dimension reservoir simulation. Therefore, the study of fast forward and inverse method of 3-dimension reservoir simulation is a cutting-edge problem, takes on important realistic signification and application value. In this dissertation, 2-dimension and 3-dimension fluid equations in porous media are discretized by finite difference, obtaining finite difference equations to meet the inner boundary conditions by Peaceman's equations, giving successive over relaxation iteration of 3-dimension fluid equations in porous media and the dimensional analysis. Several equation-solving methods are compared in common use, analyzing its convergence and convergence rate. The alternating direction implicit procedure of 2-dimension has been turned into successive over relaxation iteration of alternating direction implicit procedure of 3-dimension fluid equations in porous media, which possesses the virtues of fast computing speed, needing small memory of computer, good adaptability for heterogeneous media and fast convergence rate. The geological model of channel-sandy reservoir has been generated with the help of stochastic simulation technique, whose cross sections of channel-sandy reservoir are parabolic shapes. This method makes the hard data commendably meet, very suit for geological modeling of containing complex boundary surface reservoir. To verify reliability of the method, theoretical solution and numerical solution are compared by simplifying model of 3-dimension fluid equations in porous media, whose results show that the only difference of the two pressure curves is that the numerical solution is lower than theoretical at the wellbore in the same space. It proves that using finite difference to solve fluid equations in porous media is reliable. As numerical examples of 3-dimension heterogeneous reservoir of the single-well and multi-well, the pressure distributions have been computed respectively, which show the pressure distributions there are clearly difference as difference of the permeabilities is greater than one order of magnitude, otherwise there are no clearly difference. As application, the pressure distribution of the channel-sandy reservoir have been computed, which indicates that the space distribution of pressure strongly relies on the direction of permeability, and is sensitive for space distributions of permeability. In this dissertation, the Peaceman's equations have been modified into solving vertical well problem and horizontal well problem simultaneously. In porous media, a 3D layer reservoir in which contain vertical wells and horizontal wells has been calculated with iteration. For channel-sandy reservoir in which there are also vertical wells and horizontal wells, a 3D transient heterogeneous fluid equation has been discretized. As an example, the space distribution of pressure has been calculated with iteration. The results of examples are accord with the fact, which shows the modification of Peaceman's equation is correct. The problem has been solved in the space where there are vertical and horizontal wells. In the dissertation, the nonuniform grid permeability integration equation upscaling method, the nonuniform grid 2D flow rate upscaling method and the nonuniform grid 3D flow rate upscaling method have been studied respectively. In those methods, they enhance computing speed greatly, but the computing speed of 3D flow rate upscaling method is faster than that of 2D flow rate upscaling method, and the precision of 3D flow rate upscaling method is better than that of 2D flow rate upscaling method. The results also show that the solutions of upscaling method are very approximating to that of fine grid blocks. In this paper, 4 methods of fast adaptive nonuniform grid upscaling method of 3D fluid equations in porous media have been put forward, and applied to calculate 3D heterogeneous reservoir and channel-sandy reservoir, whose computing results show that the solutions of nonuniform adaptive upscaling method of 3D heterogeneous fluid equations in porous media are very approximating to that of fine grid blocks in the regions the permeability or porosity being abnormity and very approximating to that of coarsen grid blocks in the other region, however, the computing speed of adaptive upscaling method is 100 times faster than that of fine grid block method. The formula of sensitivity coefficients are derived from initial boundary value problems of fluid equations in porous media by Green's reciprocity principle. The sensitivity coefficients of wellbore pressure to permeability parameters are given by Peaceman's equation and calculated by means of numerical calculation method of 3D transient anisotropic fluid equation in porous media and verified by direct method. The computing results are in excellent agreement with those obtained by the direct method, which shows feasibility of the method. In the dissertation, the calculating examples are also given for 3D reservoir, channel-sandy reservoir and 3D multi-well reservoir, whose numerical results indicate: around the well hole, the value of the sensitivity coefficients of permeability is very large, the value of the sensitivity coefficients of porosity is very large too, but the sensitivity coefficients of porosity is much less than the sensitivity coefficients of permeability, so that the effect of the sensitivity coefficients of permeability for inversion of reservoir parameters is much greater than that of the sensitivity coefficients of porosity. Because computing the sensitivity coefficients needs to call twice the program of reservoir simulation in one iteration, realizing inversion of reservoir parameters must be sustained by the fast forward method. Using the sensitivity coefficients of permeability and porosity, conditioned on observed valley erosion thickness in wells (hard data), the inversion of the permeabilities and porosities in the homogeneous reservoir, homogeneous reservoir only along the certain direction and block reservoir are implemented by Gauss-Newton method or conjugate gradient method respectively. The results of our examples are very approximating to the real data of permeability and porosity, but the convergence rate of conjugate gradient method is much faster than that of Gauss-Newton method.

沾化凹陷泥质岩油气藏成藏特征研究

Mudstone reservoir is a subtle reservoir with extremely inhomogeneous, whose formation is greatly related to the existence of fracture. For this kind of reservoir, mudstone is oil source rock, cover rock and reservoir strata, reservoir type is various, attitude of oil layer changes greatly, and the distribution of oil and gas is different from igneous or clastic rock reservoir as well as from carbonate reservoir of self-producing and self-containing of oil and gas. No mature experience has been obtained in the description, exploration and development of the reservoir by far. Taking Zhanhua depression as an example, we studied in this thesis the tectonic evolution, deposit characteristics, diagenesis, hydrocarbon formation, abnormal formation pressure, forming of fissure in mudstone reservoir, etc. on the basis of core analysis, physical simulation, numerical simulation, integrated study of well logging and geophysical data, and systematically analyzed the developing and distributing of mudstone fissure reservoir and set up a geological model for the formation of mudstone fissure reservoir, and predicted possible fractural zone in studied area. Mudstone reservoir mainly distributed on the thrown side of sedimentary fault along the sloping area of the petroleum generatiion depression in Zhanhua depression. Growing fault controlled subsidence and sedimentation. Both the rate of subsidence and thickness of mudstone are great on the thrown side of growing fault, which result in the formation of surpressure in the area. The unlocking of fault which leads to the pressure discharges and the upward conduct of below stratum, also makes for the surpressure in mudstone. In Zhanhua depression, mudstone reservior mainly developed in sub-compacted stratum in the third segment of Shahejie formation, which is the best oil source rock because of its wide spread in distribution, great in thickness, and rich in organic matter, and rock types of which are oil source mudstone and shale of deep water or semi-deep water sediment in lacustrine facies. It revealed from core analysis that the stratum is rich in limestone, and consists of lamina of dark mudstone and that of light grey limestone alternately, such rock assemblage is in favor of high pressure and fracture in the process of hydrocarbon generation. Fracture of mudstone in the third segment of Shahejie formation was divided into structure fracture, hydrocarbon generation fracture and compound fracture and six secondary types of fracture for the fist time according to the cause of their formation in the thesis. Structural fracture is formed by tectonic movement such as fold or fault, which develops mainly near the faults, especially in the protrude area and the edge of faults, such fracture has obvious directivity, and tend to have more width and extension in length and obvious direction, and was developed periodically, discontinuously in time and successively as the result of multi-tectonic movement in studied area. Hydrocarbon generation fracture was formed in the process of hydrocarbon generation, the fracture is numerous in number and extensively in distribution, but the scale of it is always small and belongs to microfracture. The compound fracture is the result of both tectonic movement and hydrocarbon forming process. The combination of above fractures in time and space forms the three dimension reservoir space network of mudstone, which satellites with abnormal pressure zone in plane distribution and relates to sedimentary faces, rock combination, organic content, structural evolution, and high pressure, etc.. In Zhanhua depression, the mudstone of third segment in shahejie formation corresponds with a set of seismic reflection with better continuous. When mudstone containing oil and gas of abnormal high pressure, the seismic waveform would change as a result of absorb of oil and gas to the high-frequency composition of seismic reflection, and decrease of seismic reflection frequency resulted from the breakage of mudstone structure. The author solved the problem of mudstone reservoir predicting to some degree through the use of coherent data analysis in Zhanhua depression. Numerical modeling of basin has been used to simulate the ancient liquid pressure field in Zhanhua depression, to quantitative analysis the main controlling factor (such as uncompaction, tectonic movement, hydrocarbon generation) to surpressure in mudstone. Combined with factual geologic information and references, we analyzed the characteristic of basin evolution and factors influence the pressure field, and employed numerical modeling of liquid pressure evolution in 1-D and 2-D section, modeled and analyzed the forming and evolution of pressure in plane for main position in different periods, and made a conclusion that the main factors for surpressure in studied area are tectonic movement, uncompaction and hydrocarbon generation process. In Zhanhua depression, the valid fracture zone in mudstone was mainly formed in the last stage of Dongying movement, the mudstone in the third segment of Shahejie formation turn into fastigium for oil generation and migration in Guantao stage, and oil and gas were preserved since the end of the stage. Tectonic movement was weak after oil and gas to be preserved, and such made for the preserve of oil and gas. The forming of fractured mudstone reservoir can be divided into four different stages, i.e. deposition of muddy oil source rock, draining off water by compacting to producing hydrocarbon, forming of valid fracture and collecting of oil, forming of fracture reservoir. Combined with other regional geologic information, we predicted four prior mudstone fracture reservoirs, which measured 18km2 in area and 1200 X 104t in geological reserves.