古潜山油藏非均质三维地质建模

The discovery of the highly productive Renqiu buried hill reservoir in Bohai Bay Basin in 1975 started the high tide of finding buried hill reservoirs in China and their research. As the advance of E&P technologies, the study of buried hill reservoir in China had a qualitative leap. The reservoir description and some other aspects of development have reached or approached to the international leading level. However, some core techniques for reservoir study such as structure & faulting system study, formation prediction and connection study and heterogeneous model's construction could not completely carry out the quantitative or accurate reservoir description, e. g. the areal distribution of porosity, permeability and oil saturation. Especially, the modeling for reservoir simulation is still wandering in the stage of simplicity. The inaccurate understanding of geology could not derive 3D heterogeneous geological model that can reveal the actual underground situation thus could not design practical and feasible oilfield development plan. Therefore, the problems of low oil recovery rate, low recovery factor and poor development effectiveness have not been solved. The poor connection of the reservoir determined that waterflooding could not get good development effect and the production had to depend on the reservoir elastic energy, and this will bring big difficulty for development modification and improvement of oil recovery. This study formed a series of techniques for heterogeneous model research that can be used to construct heterogeneous model consistent with the reservoir geology. Thus the development effectiveness, success ratio of drilling and percent of producing reserves can be enhanced. This study can make the development of buried hill reservoir be of high recovery rate and high effect. The achievements of this study are as follows: 1. Evaluated the resources, summarized the geological characteristics and carried out the reservoir classification of the buried hill reservoirs in Shengli petroliferous area; 2. Established the markers for stratigraphical correlation and formed the correlation method for complex buried hill reservoirs; 3. Analyzed the structural features of the buried hill reservoirs, finished the structure interpretation and study of faulting system using synthetic seismograms, horizontal slices and coherent analysis, and clarified structural development history of the buried hill reservoirs in Shengli petroliferous area; 4. Determined the 3 classes and 7 types of pore space and the main pore space type, the logging response characteristics and the FMI logging identified difference between artificial and natural fractures by the comprehensive usage of core analysis, other lab analyses, conventional logging, FMI logging and CMR logging; 5. Determined the factors controlled the growth of the fractures, vugs and cavities, proposed the main formation prediction method for buried hill reservoir and analyzed their technical principium and applicability, and formed the seismic method and process for buried hill reservoir description; 6. Established the reserve calculation method for buried hill reservoirs, i. e. the reserves of fractures and matrix are calculated separately; the recoverable reserves are calculated by decline method and are classified by the SPE criteria; 7. Studied restraining barriers and the sealing of the faults thus clarified the oil-bearing formations of the buried hill reservoirs, and verified the multiple reservoir forming theory; 8. Formed reasonable procedure of buried hill reservoir study; 9. Formed the 3 D modeling technology for buried hill reservoirs; 10. Studied a number of buried hill blocks on the aspects of reservoir description, reservoir engineering and development plan optimization based on the above research and the profit and social effect are remarkable.

中-新元古代海洋的硫和碳同位素的演化

Data on seawater carbon isotope in the Mesoproterozoic and Neoproterozoic is abundant. However, the sulfur isotopic age curve of seawater sulfates determined through the analysis of sulfur isotopic composition of marine evaporite is uncertain in the Mesoproterozoic and Neoproterozoic since evaporites are generally rare in Precambrian. The Mesoproterozoic and Neoproterozoic Carbonate Formations preserve not only the carbon isotopic records, but also the sulfur isotopic records of coeval seawater in the Huabei Platform and the Yangtze Platform, China. Sulfur isotopic composition can be determined by the extraction of trace sulfate from carbonate samples. Successive measurements of sulfur and carbon isotopic compositions of carbonate samples from the Mesoproterozoic and Neoproterozoic strata in the Huabei Platform and the Yangtze Platform was accomplished through the extracting of trace sulfate from carbonates. Sulfur and carbon isotopic compositions of coeval seawater were obtained from analytical results of sulfur and carbon isotopes of the same sample without diagenetic alteration. The high-resolution age curve of sulfur isotope given in this paper may reflect the trend of variations in sulfur isotope composition of seawater sulfates during the Mesoproterozoic and Neoproterozoic. It can be correlated with the characteristics of variation in age curve of carbon isotope of coeval seawater carbonates. The δ34S values of seawater varied from +10.3-37.0‰ during the Mesoproterozoic, which took on oscillated variation on the whole. The δ34S values took on high values in the Mesoproterozoic Chuanlinggou stage, Tuanshanzi stage Tieling stage and in Neoproterozoic Jing'eryu stage. The average of those was about +30‰. The sulfates have low δ34S values in the Mesoproterozoic Yangzhuang stage and Hongshuizhuang stage, The average of those was all lower than +20‰. There occured large-amplitude changs in δ34S values of seawater during the Mesoproterozoic. Large-amplitude oscillate of 534S values occured in the intervals of 1600~1400Ma and 1300~1200Ma. The δ13C values of seawater are mostly negative in Changcheng stage of late Paleoproterozoic, -0 ± 1‰ range in Jixian stage of Mesoproterozoic , and the positive 2±2‰ commonly in early Neoproterozoic Jing'eryu stage. From 1000 Ma to 900 Ma, about 108 years interval of oceanic 513C record is shortage. At the end of Paleoproterozoic (1700 - 1600 Ma), the oceanic 813C values change from -3‰ to 0‰, but strongly oscillate near 1600 Ma. Two larger variations of seawater 513C values occur in the Mesoproterozoic: one is a cycle of about 4%o happens at ca. 1400 Ma; another is rise from >2‰ to>5‰ at ca. 1250 Ma and then become stable at the near 1000 Ma. There appears a large positive excursion over +20‰ in 534S value of ancient seawater sulfates in the early Doushantuo stage. Simultaneously, 8 C values of ancient seawater occur a positive excursion reaching 10‰. These allow δ4S values and 513C values to reach high values of+51.7‰ and +6.9‰, respectively. The range of variation in 834S values of seawater is relatively narrow and 513C values are quite high in the middle Doushantuo stage. Then, δ34S values of seawater become oscillating, the same happens in δ13C values. Negative excursions in 834S values and 813C values occur simultaneously at the end of the Doushantuo stage, and the minimum of δ34S values and δ13C values dropped to -11.3‰ and -5.7‰, respectively. The ancient seawater in the Dengying stage has high δS values and δ13C values. Most of the δ34S values of the trace sulfate samples varied between +23.6‰ and +37.9‰ except two boundaries of the Dengying Formation, and the S13C values of the carbonate samples of the Dengying Formation varied between +0.5‰ and +5.0‰. There appeared large negative excursion in 834S values and δ13C values of ancient seawater at the bounder of Precambrian-Cambrian. The isotopic characteristics of sulfur and carbon implicated that the organic productivity and isotopic fractionation caused by biology were low and the palaeoceanic environment was quite unstable during the Mesoproterozoic. The increase and subsequent oscillation of seawater δ13C value occurred from 1700 to 1600 Ma and near 1300 Ma may be responsible to the two global tectonic events happened at coeval time. The characteristics of variation in sulfur and carbon isotopes of ancient seawater imply strong changes in oceanic environment, which became beneficial to inhabitation and propagation of organism. The organic production and the burial rate of organic carbon once reached a quite high level during the Doushantuo stage. However, the state of environment became unstable that means the global climate and the environment possibly were fluctuating and reiterating after the global glaciation. The negative excursions of S34S values and δ13C values occurring at the end of the Doushantuo stage represent a global event, which might be relative to the oxidation of deep seawater. The isotopic characteristics of sulfur and carbon implicated that there were a high organic productivity and a high burial rate of organic carbon in the Dengying stage. It is obvious that the palaeoceanic environment in Dengying stage was stable corresponding and beneficial for biology to inhabit and propagate except for the two boundaries. The tendency of sulfur and carbon isotopic variations maybe resulted from the gradual oxygenation of ocean environment during the Dengying stage. It has been reported that the secular variations of the sulfur isotopic compositions in seawater was negative correlated with that of carbon isotopic compositions. However, our results show that it is not the case. They were negatively correlated in some intervals and positively in some other intervals of the Mesoproterozoic and Neoproterozoic. The difference in correlation may be associated with the changes in conditions of redox in oceanic environment, e.g. sharp change of the oxidation-reduction interface. The strong changes in global environment may induce the abnormality to occur in the biogeo chemical S and C cycles in the ocean and accordingly sharp Variations in isotopic composition of seawater sulfur and carbon during the Mesoproterozoic and Neoproterozoic. Simultaneously, the global tectonism caused large changes of 87Sr/86Sr ratios. The leading factor that causes the variation in isotopic composition is different in the different intervals of the Mesoproterozoic and Neoproterozoic. Thus, there may exist different models of the biogeochemical S and C cycles in the ocean during the Mesoproterozoic and Neoproterozoic.

2200万年来风尘沉积的粒度变化及古气候意义

Three eolian deposit formations, including Quaternary loess (QL, Liu et al.3 1985), Hipparion red earth (HRE, also called red clay, Liu et al., 1985) and Miocene loess (ML, Guo et al., 2002) constitute a set of unique paleoclimatic archives in northern China dated back to at least 22Ma ago. The Miocene loess is a recently discovered loess-soil sequence. Detailed investigation has been made on its origin, chronology and paleoclimatic significance (Guo et al., 2002), but further work is still needed to obtain detailed paleoclimate information, and mechanical links behind paleoclimatic changes. In this study, grain size analysis of QL, HRE and ML has been conducted on two sections: Qinan and Xifeng. The objective is focused on comparison of the grain size distribution characteristics (GSDC) among different eolian deposit formations, and reconstruction of the Asian monsoon circulation in the past 22 Ma. Results show that GDSC of ML resembles that of QL and HRE, and GDSC of ML is especially similar to HRE. Both ML and HRE contain a significant proportion of fine fraction, however, QL has a large amount of coarse sediments. This is mainly due to that the wind system transported aeolian dust was weaker in the late Tertiary than that in the Quaternary. Grain size difference between loess and paleosol in ML is much smaller than that in QL, indicating that the climatic fluctuations during the late Tertiary were much smaller than that happened in the Quaternary The grain size records of the past 22 Ma reveal several evolution phases of the Asia winter monsoon. -2.7 Ma BP is the most important boundary in the process of the winter monsoon evolution: the wind strength have significantly enhanced since 2.7 Ma ago. During a period between -22.0 and -3.6 Ma, three periods with relatively stronger winter monsoon are recorded in the QA-I section, between 21.2 and 19.9, and 16.0 and 13.3, and 8.7 to 6.9 Ma, respectively. From 3.6 to 2.7, the winter monsoon was enhanced gradually. In the Miocene time, the intensified winter monsoon phases (between 21.2 and 19.9, and 16.0 and 13.3, and 8.7 and 6.9 Ma) seemed to have a close relationship with the uplift of the Tibetan Plateau and/or the ongoing global cooling, but the forcing mechanism behind the Asia winter monsoon evolution need to be further investigated. During the Pliocene-Pleistocene time, the Asia winter monsoon strengthened at 3.6 and 2.7Ma ago are in good agreement with the ongoing global cooling and the Arctic ice sheet development. In the mean time, much evidence suggests that an intense uplift of the Tibetan Plateau occurred at ~3.6 Ma, which is synchronous with a major increase in Asia winter monsoon. Therefore, two major factors may be invoked to explain the winter monsoon enhancement: Arctic ice sheet development and Tibetan uplift. We propose that changes in location and intensity of the Siberian-Mongolian high that were caused by the Tibetan uplift and Arctic ice sheet development might be an important factor for Asian winter monsoon evolution in the Pliocene-Pleistocene.

龙滩工程左岸边坡监测资料分析与稳定性评价研究

On the basis of the geological analysis and rock mass toppling deformation and failure mechanism analysis of Longtan engineering left bank slope, the synthetic space-time analysis and influence factors analysis on the surface monitoring data and deep rock mass monitoring data of B-zone of left bank slope are carried on. At the same time, based on the monitoring data analysis in conjunction with the predecessor's mechanics analysis results, the deformation state of B-zone of the left bank slope is discussed and its stability is synthetically evaluated. The detailed research contents and results are as following: According to the monitoring drill histogram analysis of Longtan engineering left bank slope, numerical simulation analysis and model experimentation analysis of bedded counter-inclined steep slope, a new type of toppling deformation and failure mode is proposed, that is "up-slope warping". Then the deformation and failure mode of bedded counter-inclined steep slope is summarized as "down-slope toppling" type, "up-slope warping" type and "complex fold" type. On the basis of synthetic space-time analysis to surface monitoring data and deep rock mass deformation monitoring data of B-zone of Longtan left bank slope;, we can get the conclusion that there exists potential instability rock mass over 520m altitude, especially over 560m altitude of slope B, and the rock mass of around strong-weathering line or creep rock mass breaking band controls the deformation of the whole slope. 1. According to the synthetic space-time analysis and influence factors analysis to the surface monitoring data of B-zone of Longtan left bank slope, a dynamical index, accumulative total acceleration index, which is used to analyze the influence factors of slope surface deformation, is raised. The principle and method of accumulative acceleration index are explained, and the index can be used for the influence factors analysis of the similar slope. 2. Summarize the results of geologic analysis, monitoring analysis and mechanics analysis, the following conclusion can be gotten: the stability of B-zone of the slope is basically good. However, on the condition of drainage and slope toe loading engineering, there is still some creep deformation in the rock mass over 520m altitude, especially over 560m altitude. So, better measures of the monitoring and timely maintenance of the drainage system are suggested in the paper.

库车坳陷第三纪磁性地层和磁组构的研究

The Tertiary Kuche depression, also known as one of the most economically important oil and gas prospecting regions in the Tarim Basin, is a foreland basin formed by flexural subsidence resulting from the southward thrusting of the southern Tianshan. Detailed geological and geophysical studies on the Tertiary sequence of the Kuche region would provide constraints not only on the Mesozoic and Cenozoic deformation patterns, and its controlling on the formation and distribution of oil and gas of the Kuche depression, but also on the India-Asia postcollisional uplift and continental deformation of the Tianshan regions. In this thesis, the Tertiary sequence of the Kuche depression, which is composed of the upper Kumukeliemu Formation, Suweiyi Formation, Jidike Formation, Kangcun Formation, and Kuche Formation, was selected for magnetostratigraphic and anisotropy of magnetic susceptibility (AMS) study. Among a total of 697 collected sites, 688 sites gave magnetic fabric results and 686 sites yielded reliable paleomagnetic results, building up magenetostratigraphy for the Tertiary sequence of the Kuche depression. By correlating with international geomagnetic polarity timescale, the followings were concluded: 1) the time interval of the sampling sections is some 31-8.1 Ma; 2) the boundaries for Kangcun/Jidike Formations, Jidike/Suweiyi Formations, and Suweiyi/Kumukeliemu Formations are at about 13.5 Ma, 26 Ma and 29 Ma respectively; 3) the Jidike Formation might be formed in an interval between the late Upper Oligocene and middle Miocene, whereas the Suweiyi Formation was formed in the Upper Oligocene, in turn, the Paleogene/Neogene boundary is most likely to be at the bottom of the Jidike Formation in the Kuche depression; 4) the dramatic transition from the marine/lacustrine to fluvial/alluvial facies in the Suweiyi Formation and the lower Jidike Formation of the Kuche depression, which may indicate the initiation of Cenozoic thrusting in the Kuche depression, is probably occurred in the late Oligocene, i.e. at about 28-25 Ma. On the other hand, systematic differences of the AMS parameters, such as the corrected anisotropy of magnetic susceptibility, shape parameter of ellipsoids, q factors, angles between the intermediate and minimum magnetic susceptibility directions and bedding, indicate that a dramatic transition of tectonic strain is most likely to occur at about 16 Ma (middle Miocene) in the Kuche depression. Taking into account the initiation o f C enozoic thrusting in the Kuche d epression i s at about 2 8-25 M a, we argue that the intensive Cenozoic thrusting in the Kuche depression is mainly happened during a period between the late Oligocene to middle Miocene.

高含水油藏复合驱剩油分布研究-以羊三木油田馆陶组油藏复合驱为例

Based on the study of the combined flooding test block of Guantao formation in Third faulted block of Yangsanmu oil field, this paper carries out the integration of reservoir precise characterization for very high water cut reservoir, establishes precise 3D geologic model for high water cut development period and states the changing law of the reservoir architecture dtiring development by combined flooding. Then, by subdivided the thick oil reservoir, the study of remaining oil saturation monitoring in fiber glass cased well and tracer monitoring is developed. According the study of multiple constrained combined flooding reservoir numerical simulation, remaining oil distribution are predicted, the methods architecture of predicting remaining oil distribution are established for fluvial facies reservoir at late development stage, develops plan is designed and adjustment associating technologies for enhancing oil recovery. On these base, related measures for tapping the potential are given, it is verified and optimized through the field former test and the good economic effect is achieved . The major achievements of this paper are as follows. The changing law of the reservoir architecture and it's property parameters is revealed, The result indicates that the temperature-pressure of the injecting material and the interaction effect of the injecting material and reservoir petrography are the main factors of the dynamic changes of the reservoir architecture. The quantitative reservoir geologic model, which is tallied with dynamic reservoir parameters of the study area, is established. Subdivided the thick oil reservoir is very important for the study of the remaining oil distribution within the thick oil reservoir. Subdivided the thick oil reservoir technology, which consists of six technologies as follow: micro-cyclic divided, flow unit method, architectural element method, high resolution log technology, high resolution-process technology for normal logging data and using the production data is presented. 3. It is established dynamic monitoring system of remaining oil saturation quantitative research which are inner and interlayer remaining oil saturation from time-lapse logging in fiber glass cased well, inter-well remaining oil saturation from the technology of isotopic tracer monitoring technology, and 4d remaining oil saturation distribution from combined flooding numerical modeling integrated by production datao The forming mechanism of remaining oil for polymer flooding and alkali/polymer combined flooding is clarified, and the plane and vertical distribution law of remaining oil after combined flooding is revealed. Predicting methods and technologies for the combined flooding reservoir of fluvial facies is developed. Combined flooding has been achieved good displacement result in the pilot of Third fault block in Yangsanmu oil field, and accumulated types of important parameters and optimum plans, this technology of combined flooding is expected to increase recovery ratio by 4.77%.

陕甘宁盆地储层物性预测及油气综合判识方法研究

The thesis show you seven interpretation models of erosion of MAs1+2 in the west of prospect and eroded gully of middle and east, form the interpretation technique how to built up eroded gully of the Ordovician top, according of the practical demand of oil-gas exploration in the ShanGanNing basin, using seismic information, combining well logging and well drilling data, Carefully analyzing geologic deposition background and well logging data, through a great quantity forward and inversion for geologic model and combination geologic model with seismic section. Related to research of reservoir absorption in the ShanGanNing basin, it firstly introduces PRONY transformation multidimensional filter. It can simultaneously express relationship of frequency and absorption decay coefficient, better than FUSAIPU analysis method; PRONY filter have obtain the better effect in the gas field of ZhenChuanBao in the ShanBei area after adopting PRONY filtering method to predict reservoir absorption, by analyzing fixed well and prediction of non-well drilling. In the ShanGanNing basin, general seismic inversion method can produce evident different results or misunderstanding because wave impedance and lithology, physical property, gas property are not sole, especially while have little impedance contrast and even have contract direction; the author carefully analyzes multi-parameter inversion technique, add natural gamma ray and natural potential and other parameter combined making model inversion method according of theory of seismic inversion and applying reservoir velocity and wave impedance information at last, we get the more directly reservoir physical property parameter, judging reservoir physical property is more exact. In accordance with geologic, seismic feature of Shan basin, the thesis conclude Ordovician system top erosion interpretation technology with ChangQing character, and reservoir thickness prediction technique combining inversion technique with wave character analysis, Reservoir physical property that is mainly absorption factor analysis and multi-parameter inversion and oil-gas prediction technology. These technologies obtain the better result in the oil-gas field exploration and have formed comprehensive research method and technology series with ShanGanNing character.

胜利油田低渗透砂岩油藏非均质地质建模

Guided by geological theories, the author analyzed factual informations and applied advanced technologies including logging reinterpretation, predicting of fractal-based fracture network system and stochastic modeling to the low permeable sandstone reservoirs in Shengli oilfield. A new technology suitable for precious geological research and 3D heterogeneity modeling was formed through studies of strata precious correlation, relation between tectonic evolution and fractural distribution, the control and modification of reservoirs diagenesis, logging interpretation mathematical model, reservoir heterogeneity, and so on. The main research achievements are as follows: (1) Proposed four categories of low permeable reservoirs, which were preferable, general, unusual and super low permeable reservoir, respectively; (2) Discussed ten geological features of the low permeable reservoirs in Shengli area; (3) Classified turbidite fan of Es_3 member of the Area 3 in Bonan oilfield into nine types of lithological facies, and established the facies sequences and patterns; (4) Recognized that the main diagenesis were compaction, cementation and dissolution, among which the percent compaction was up to 50%~90%; (5) Divided the pore space in ES_3 member reservoir into secondary pores with dissolved carbonate cement and residual intergranular pores strongly compacted and cemented; (6) Established logging interpretation mathematical model guided by facies- control modeling theory; (7) Predicted the fracture distribution in barriers using fractal method; (8) Constructed reservoir structural model by deterministic method and the 3D model of reservoir parameters by stochastic method; (9) Applied permeability magnitudes and directions to describe the fractures' effect on fluid flow, and presented four different fractural configurations and their influence on permeability; (10) Developed 3D modeling technology for the low permeable sandstone reservoirs. The research provided reliable geological foundation for the establishment and modification of development plans in low permeable sandstone reservoirs, improved the development effect and produced more reserves, which provided technical support for the stable and sustained development of Shengli Oilfield.

东营凹陷沙三中亚段湖相浊积岩储层特征及油气成藏规律研究

Turbidity sandstone reservoirs have been an important field of hydrocarbon exploration and development in the basins all over the world, as well as in China. Lithologic pools are composed of turbidity sandstones and other sandstones are frequently found in the Jiyang Depression that is a Mesozoic-Cenozoic non-marine oil-bearing basin. The Dongying Sag lies in the sedimentary center of the basin. The subtle traps with turbidity reservoirs are generally difficult to be predicted and described by using current techniques. The studies on turbidity reservoirs plays thus an important theoretical and theoretical practical role in exploration and development in the Jiyang Depression. The attention is, in this thesis, focused on the petrologic properties and oil accumulating behaviors in lake turbidity sedimentary systems in the middle part of the third section of Shahejie Formation in the Dongying Sag, especially in Dongxin area, which lies on the central uplift of the Sag. The paper has disclosed the origin types of turbidity sandstones, distribution pattern and controlling factors of turbidity sandstones, and set up hydrocarbon accumulation patterns of the middle part of the third section of Shahejie Formation in Dongxin, based on nonmarine high resolution sequence stratigraphy, event sedimentology and new theories of hydrocarbon forming. By studying prediction method and technology of turbidity sandstone reservoirs, using precise geological model developing, new techniques of high resolution seismic inversion constrained by logging, the paper has forecast low permeability turbidity sandstone reservoirs and pointed out advantage exploration aims to progressive exploration and development. The paper has obtained mainly many productions and acknowledges as follows: 1.Turbidity sandstone reservoirs of the third section of Shahejie Formationin Dongying Sag are formed in such specifical geological background as rift and extension of basin. The inherited Dongying delta and transgression make up many turbidity distribution areas by overlaying and joining together. The hydrocarbon migrates from depression area to adjacent turbidity sandstone continuously. Accumulation area which is sufficient in oil is formed. 2.The paper has confirmed distinguishable sign of sequence boundary , established stratigraphic framework of Dongying Sag and realized isotime stratigraphic correlation. Es3 of Dongying delta is divided into eleven stages. Among them, the second period of the lower section in Es3, the sixth period of the middle section in Es3, the third period of the upper section in Es3 correspond to eleven sedimentary isotime surface in seismic profile, namely Es3 is classified into eleven Formations. 3.According to such the features of turbidity sandstone as deep in burial, small in area, strong in subtle property, overlaying and joining together and occurring in groups, management through fault and space variations of restriction quantum are realized and the forecast precision of turbidity sandstone by using precise geological model developing, new techniques of high resolution seismic inversion constrained by logging, based on the analysis of all kinds of interwell seismic inversion techniques. 4.According to the features of low permeable turbidity sandstone reservoirs, new method of log interpretation model is put forward. At the same time, distinguish technology of familiar low resistivity oil layer in the turbidity sandstone reservoirs is studied based on petrophysical laboratory work and "four properties" interrelationship between lithological physical Jogging and bearing hydrocarbon properties. Log interpretation model and reservoir index interpretation model of low resistivity oil layer are set up. So the log interpretation precision is improved. 5.The evolution law and its difference of the turbidity sandstone are embodies as follows: the source of sediments come from the south and east of the study area in the middle period of Es3. East source of sediments is pushed from west to east. However, the south source supply of sediments in the early and middle period of Es3 is in full, especially in Es3. subsequently, the supply is decreased gradually. Turbidity fan moves back toward the south and the size of fan is minished accordingly. The characteristic of turbidity sandstone in Dongying Sag is different in different structural positions. Dongxin in the middle-east of the central lift and Niuzhuang Sag He in Dongying delta front and prodelta deep lake subfacies. Although the turbidity sandstone of the two areas root in the Dongying delta sedimentary system, the sand body has different remarkably characteristic. 6.The sedimentary model of the turbiditys in study area have three types as follows: (1) collapse turbidity fan in respect of delta; (2) fault trench turbidity fan; (3) other types of microturbidity sandstone. Middle fan and outer fan, can be found mainly in sublacustrine fan. Middle fan includes braided channel microfacies, central microfacies and braided interchannel microfacies, which is main prospecting oil-bearing subfacies. The middle section of the third section of Shahejie Formation in study area (for example the central lift) can be divided into middle-lower and upper part. The middle-lower part is characteristic of turbidity fan. The upper part is sedimented mainly by delta-collapse fan. 7.The turbidity reservoirs of the middle part of the third section of Shahejie Formation in study area characterize by low maturity both in component and texture, strong in diagenesis and low in permeability. The reservoir can be classified into four types. Type III is the body of reservoir and comprises two types of H a and HI b. M a belongs to middle porosity - low permeability reservoir and distributes in the central lift. Hlb belongs to low porosity - low permeability and distributes in Haojia region. 8.A11 single sand body of lens turbidity reservoir of the middle part of the third section of Shahejie Formation in study area are surrounded by thick dark source rocks. The oil-water system is complex and behaves that every sandstone is single seal unit. The water body is 1/3-1-5 of the sand body. The edge water is not active. The gas exists in the top of reservoir in the form of mixed gas. For far-range turbidity fan with big scale channel, the area and volume of sand body is large and the gap is big in oil packing degree. There are lots of edge water and bottom water, and the latter increases rapidly during the course of development. 9.By exerting the modern hydrocarbon forming theories, the third section of Shahejie Formation in study area belongs to abnormally pressured fluid compartment. The lithological reservoir of the third section of Shahejie Formation is formed in the compartment. The reservoir-formed dynamic system belongs to lower self-source enclosed type. The result and the practice indicate that the form and accumulation of lithological oil reservoirs are controlled by the temperature and pressure of stratum, microfacies, thickness of sand body, fault and reservoir heterogeneity. 10. Based on studies above, the emphases focus on in south and north part of Dongying structure, west Dongxin region and south part Xinzhen structure in the application of production. The practice proves that the turbidity sandstone reservoirs in Ying 11 block and the fault-lithological reservoirs in Xin 133 block have been obtained significant breakthrough. The next target is still sandstone groups of the third section of Shahejie Formation in the bordering areas of Dongxin region for instance Xin 149 area, He 89 area, Ying 8 area etc.

文33块沙二下段油藏动态模型及剩余油预测

By applying synthetically multi-subject theories, methods and technology, such as petroleum geology, sedimentology, seep mechanics, geochemistry, geophysics and so on; and by making full use of computer; combining quantity and quality, macroscopic and microscopic, intensive static and active description, comprehensive studying and physical modeling, 3 dimension and 4 dimension description; the paper took Wen-33 block of Zhongyuan oil field as an example; and studied reservoir macroscopic and microscopic parameter changing rule and evolve mechanics in different water-blood stage. The reservoir dynamic model and remaining-oil distribution mode was established, and several results were achieved as follows: (1) Three types of parameter gaining, optimizing and whole data body of Wen33th reservoir were established. Strata framework, structure framework, reservoir types and distribution of Wen33th reservoir were discussed. Reservoir genesis types, space distribution law and evolve rule of Wen33th reservoir were explained. 4D dynamic model of macroscopic parameter of reservoir flow dynamic geologic function of Wen33th reservoir was established. The macroscopic remaining-oil distribution and control factor was revealed. The models of the microscopic matrix field, pore-throat network field, fluid field, clay mineral field of Wen-33 block were established. The characters, changing rules and controlled factors in different water stage were revealed. The evolve rule and mechanics of petroleum fluid field in Wen-33 block reservoir were revealed. Macroscopic and microscopic remaining oil distribution mode of Wen-33 block were established. Seven types, namely 12 shapes of dynamic model of microscopic remaining oil were discussed, and the distribution of mover remaining oil was predicted. Emulation model: mathematical model and prediction model of Wen-33 block were established. The changing mechanics of reservoir parameter and distribution of remaining-oil were predicted. Firstly, the paper putting forward that the dynamic geologic function of petroleum development is the factor of controlling remaining-oil, which is the main factor leading to matrix field, network field, clay mineral field, fluid field, physic and chemical field, stress field and fluid field forming and evolving. (10) A set of theories, methods and technologies of investigating, describing, characterizing and predicting complex fault-block petroleum were developed.

油藏数值模拟中的网格粗化和细化技术研究

The grid is a foundation of reservoir description and reservoir simulation. The scale of grid size is vital influence for the precision of reservoir simulation the gridding of reservoir parameters require reasonable interpolation method with computing quickly and accurately. The improved distant weighted interpolation method has many properties, such as logical data points selection, exact interpolation, less calculation and simply programming, and its application can improve the precision of reservoir description and reservoir simulation. The Fractal geologic statistics describes scientifically the distribution law of various geological properties in reservoir. The Fractal interpolation method is applied in grid interpolation of reservoir parameters, and the result more accorded with the geological property and configuration of reservoir, and improved the rationality and quality of interpolation calculation. Incorporating the improved distant weighted interpolation method with Fractal interpolation method during mathematical model of grid-upscaling and grid-downscaling, the softwares of GROUGH(grid-upscaling) and GFINE (grid-downscaling) were developed aiming at the questions of grid-upscaling and grid-downscaling in reservoir description and reservoir simulation. The softwares of GROUGH and GFINE initial applied in the research of fined and large-scale reservoir simulation. It obtained fined distribution of remaining oil applying grid-upscaling and grid-downscaling technique in fined reservoir simulation of Es21-2 Shengtuo oilfield, and provided strongly and scientific basis for integral and comprehensive adjustment. It's a giant tertiary oil recovery pilot area in the alkaline/surfactant/polymer flooding pilot area of west district of Gudao oilfield, and first realized fined reservoir simulation of chemical flooding using grid-upscaling and grid-downscaling technique. It has wide applied foreground and significant research value aiming at the technique of grid-upscaling and grid-downscaling in reservoir description and reservoir simulation.

油气藏自然电位异常与油气赋存的相关性研究

The Stack Spontaneous Potential (SSP) is a direct hydrocarbon location technology and a new hydrocarbon detection method with independent intellectual property. A subsurface hydrocarbon accumulation associated with the upward hydrocarbon micro seepage induces a relatively strong negative potential abnormal zone, of which the anomaly can be measured on the surface with specially designed instruments through careful field measuring procedures. With special software programmed according to a unique geochemical and geophysical model, the original data are analyzed, processed and interpreted on the computer, and then on a series of resulting anomaly distribution maps and/or profiles, the favorable surface locations of the hydrocarbon accumulations can be easily identified. The study of the SSP has been started since 1989, and especially from 1996 to 1997, both profile and area tests were conducted in the Daqing Oilfield. On the testing line of 15kms, there are 6 wells in total, among which some are oil-producing wells, and some are water-producing wells. The final matching ratio of the favorable oil well locations and the possible water well locations predicted by the SSP to those of known wells was up to 83 percent. In the area test, of which the acreage is 800 km2, the matching ratio compared with the existing wells was 87 percent; furthermore, regarding to wells subsequently drilled after the test, the matching ratio was 85 percent. The matching ratio in the development area is more about 10 percent than those of in exploration area. The reason is that, comparing the exploration area, the development area acreage is less and the container rocks are more simplex. In development area there is not so much interference of SSP also. Since 1997 the SSP has been tested and applied all over China to a number of hydrocarbon bearing basins and known oil fields, including the Daqing, Jiangsu, Changqing, Shengli, Nanyang, Jianghan and Zhongyuan Oilfields, only to name a few. The SSP surveys in total areas of over 10,000km2 in more than 30 regions in China so far have been completed in various exploration and development stages, the satisfactory outcomes of which have further evidenced that the dependence between the SP anomaly and abundance of hydrocarbon. Up to date, a substantial amount of successful tests and actual surveys finished in exploration and development practices have evidenced that the SSP is significantly more reliable in comparison with any other similar direct hydrocarbon indication technique generally known to the oil industry, such as the Redox. The SSP can be applied to search for almost all kinds of hydrocarbon accumulations, regardless of the type of traps, such as structural, stratigraphic, buried hill traps, and so on; however, it is interesting to be noted that the SSP seems to be particularly effective in detecting the stratigraphic oil traps according to our practices. On the other hand, there is virtually no surface geographical constrains in terms of field data acquisition, except for those water covered areas, because of the inherent characteristics of the technology itself. Furthermore, utilizing the SSP requires no special considerations to subsurface geological conditions in regard to formation resistivity, since the SSP measurements will not be influenced by either overly high or overly low resistivity of formations lying above the hydrocarbon accumulations. There are two kind of theory, of which, as we know one is called hypbyssal theory such as "Redox"[61 the other is call plutonic theory such as cracking of hydrocarbon [8][9] and natural polarization [3], to describe the mechanism of SP anomaly of oil reservoir and to indicate that the dependence between the SP abnormality and abundance of hydrocarbon has be existed theoretically/The quantitative dependence, which has not been founded due to the complicity of container rocks, be discovered during the exploration and development practices is the crux to the quantitative analysis of SP Anomaly processing. Based on the thorough study of the complex of collector rocks, every kind of thickness of collector rock can be conversed to be a standard effective thickness; the thickness is called apparent effective thickness (AET). The conversation coefficient (ai, 1=1,2,3) could be determined by the variety of every collector rock storability (CRS). The discoveration of quantitative: dependence between AET and the amplitude of SSP, in the practices of exploration and development, is a promotion for the SSP supplied in the oil exploration, and make the data analysis forward to the quantitative stage.

中国东海陆架盆地构造-热演化研究

East China Sea Shelf Basin (ECSSB), as a basin with prospect of oil & gas resource and due to its special geological location on the west margin of the west Pacific, attracts a lot of attention from many geologists in the world.Based on systematic temperature measurements, bottom hole temperature (BHT) and the oil temperature data, the geothermal gradients in the ECS SB are calculated and vary from 25 to 43°C/km, with a mean of 32.7°C/km. The geothermal gradient in Fuzhou Sag has the higher value(40.6°C/km) in Taibei Depression than that in others. The lower value (27.2 °C/km) occurs in in Xihu Depression. The middle values occurs in Jiaojiang and Lishui sags in Taibei Depression with a mean value of 34.8 °C/km. Incorporated with the measured thermal conductivity, heat flow values show that the ECSSB is characterized by present-day heat flow around 70.6mW/m2, varying between 55 and 88 mW/m2. No significant difference in heat flow is observed between the Xihu and the Taibei Depressions. These heat flow data suggest that the ECSSB is geothermally not a modem back-arc basin.Applying the paleogeothermal gradient based method, thermal history is reconstructed using vitrinite reflectance (VR) and apatite fission track (AFT) data. The results suggest that the thermal history was different in the Taibei and the Xihu depressions. Paleo-heat flow values when the pre-Tertiary formations experienced their maximum temperature at the end of the Paleocene reached a mean of 81 mW/m2 in the Taibei Depression, much higher than the present-day value. The lower Tertiary sediments in the Xihu Depression experienced maximum temperatures at the end of Oligocene and reached a mean paleo-heat flow value of 83.4 mW/m2. The time, when the paleo-heat flow reached the maximum value, suggests that the ECSSB rifted eastward.Tectonic subsidence analysis shows that the timing of the major rifting episode was different across the ECSSB. The rifting occurred from the Late Cretaceous to the early Eocene in the Taibei Depression, followed by thermal subsidence from the late Eocene to the end of Miocene. In contrast, in the Xihu Depression the initial subsidence lasted until the early Miocene and thermal subsidence to the end of Miocene. From Pliocene to the present, an accelerated subsidence took place all along the West Pacific margin of the east Asia.The thermal lithosphere thickness is determined by temperature profile in the lithosphere, the mantle adiabat or the dry basalt solidus. It indicates that the thermal lithosphere reached the thinnest thickness at the end of Eocene in the Taibei Depression and the end of Oligocene in the Xihu Depression, respectively, corresponding with a value of 57-66km and 56-64km. In Taibei Depression, the lithosphere thickness decreased 16-22km from the end of Mesozoic to Paleocene. After Paleocene, the thickness increased 13-16km and reached 71-79 km at present-day. In Xihu Depression, From the end of Oligocene to present-day, the thickness increased 10-13km and reached 69-76km at present-day. The evolution of the lithosphere thickness is associated closely with the lithosphere stretching.Combining the reconstructed thermal history and the burial history, the maturation of the Jurassic oil-source rock shows that the main hydrocarbon generation phase was in the mid-Jurassic and a secondary hydrocarbon generation occurred at the end of Paleocene. The secondary generation was controlled mainly by the tectono-thermal background during the Paleocene.

鄂尔多斯盆地苏里格气田储层次生孔隙性质与分布预测

Permian reservoir in Sulige area of Ordos Basin, on which this paper focused, belongs to fluvial-delta lithofacies. The majority formations in this area are complicated channel sand deposit with serious inhomogeneity which makes natural gas exploration be very tough in this area. This inhomogeneity can be found everywhere both in large horizontal area and vertical profile of inner and interbedded formations.This paper studied the inhomogeneity characteristic of Permian formation in sulige area of Ordos Basin according to the logging data.Correlating with core data, a criterion to distinguish different type of reservoirs by using logging data is determined after the study of logging response is done considering the diverse conditions of deposit environments, lithology and reservoir space. The characteristic relationships between the various type formations and logging responses fully and systemically are established.It investigated reservoir parameter calculation methods amply. Combining the conventional and special logging data, basing on the feature of low porosity -permeability formation of sulige area, a set of methods to calculate reservoir parameters was formed including primary porosity, secondary porosity, fracture porosity, permeability and water saturation under the conditions of both low porosity-permeability and inhomogenous reservoirs. One thing should be pay close attention is the parameter M for calculating saturation. It is found that the M in low porosity -permeability formation decreases as the porosity decrease, which is opposite to the law that M increases as the porosity decrease in the formation with intermediate to high porosity and permeability. This view has innovated the traditional theory and offered theory basis for the logging interpretation of low porosity - permeability reservoir. Meanwhile it also improved the Arqi formula theoretically and enhanced the logging interpretation accuracy and rescued a number of formations which has been thought to be hopeless according to the old theory.By using advantage logging interpretation procedure, a territorial synthetic geology evaluation to the inhomogeneous reservoir was completed basing on the single well interpretation. All the reservoir evaluation parameters including sand formation thickness, primary porosity, secondary porosity were calculated and evaluated. The rules of changing and development for sand formation thickness, sand physical properties and secondary porous were found at different formations of upper part of the Member 8 of Shihezi, lower part of the Member 8 of Shihezi, the Member 1 of shanxi and the Member 2 of shanxi individually. Evaluation and Correlation of these five formations were also completed and one conclusion was arrived: upper part of the Member 8 of Shihezi formation has the best performance followed by the lower part of the Member 8 of Shihezi, the Member 1 of shanxi and the Member 2 of shanxi formation.After studied the relationship between reservoir deposition characteristic and the natural gas richness, it is regarded that reservoir inhomogeneity is the key issue of the impaction on the natural gas. Natural gas in Sulige gas field was mainly accumulated in sands of channel bar, distributary channel and debouchure bar. Especially, the quartz sand with rich of secondary porous space has obvious better physical properties than other reservoir and usually can forms the concentration of natural gas.

双聚焦偏移速度估计与可分表示叠前弹性反演

In exploration seismology, the geologic target of oil and gas reservoir in complex medium request the high accuracy image of the structure and lithology of the medium. So the study of the prestack image and the elastic inversion of seismic wave in the complex medium come to the leading edge. The seismic response measured at the surface carries two fundamental pieces of information: the propagation effects of the medium and the reflections from the different layer boundaries in the medium. The propagation represent the low-wavenumber component of the medium, it is so-called the trend or macro layering, whereas the reflections represent the high-wavenumber component of the medium, it is called the detailed or fine layering. The result of migration velocity analysis is the resolution of the low-wavenumber component of the medium, but the prestack elastic inversion provided the resolution of the high-wavvenumber component the medium. In the dissertation, the two aspects about the migration velocity estimation and the elastic inversion have been studied.Firstly, any migration velocity analysis methods must include two basic elements: the criterion that tell us how to know whether the model parameters are correct and the updating that tell us how to update the model parameters when they are incorrect, which are effected on the properties and efficiency of the velocity estimation method. In the dissertation, a migration velocity analysis method based on the CFP technology has been presented in which the strategy of the top-down layer stripping approach are adapted to avoid the difficult of the selecting reduce .The proposed method has a advantage that the travel time errors obtained from the DTS panel are defined directly in time which is the difference with the method based on common image gather in which the residual curvature measured in depth should be converted to travel time errors.In the proposed migration velocity analysis method, the four aspects have been improved as follow:? The new parameterization of velocity model is provided in which the boundaries of layers are interpolated with the cubic spline of the control location and the velocity with a layer may change along with lateral position but the value is calculated as a segmented linear function of the velocity of the lateral control points. The proposed parameterization is suitable to updating procedure.? The analytical formulas to represent the travel time errors and the model parameters updates in the t-p domain are derived under local lateral homogeneous. The velocity estimations are iteratively computed as parametric inversion. The zero differential time shift in the DTS panel for each layer show the convergence of the velocity estimation.? The method of building initial model using the priori information is provided to improve the efficiency of velocity analysis. In the proposed method, Picking interesting events in the stacked section to define the boundaries of the layers and the results of conventional velocity analysis are used to define the velocity value of the layers? An interactive integrate software environment with the migration velocity analysis and prestack migration is built.The proposed method is firstly used to the synthetic data. The results of velocity estimation show both properties and efficiency of the velocity estimation are very good.The proposed method is also used to the field data which is the marine data set. In this example, the prestack and poststack depth migration of the data are completed using the different velocity models built with different method. The comparison between them shows that the model from the proposed method is better and improves obviously the quality of migration.In terms of the theoretical method of expressing a multi-variable function by products of single-variable functions which is suggested by Song Jian (2001), the separable expression of one-way wave operator has been studied. A optimization approximation with separable expression of the one-way wave operator is presented which easily deal with the lateral change of velocity in space and wave number domain respectively and has good approach accuracy. A new prestack depth migration algorithm based on the optimization approximation separable expression is developed and used to testing the results of velocity estimation.Secondly, according to the theory of the seismic wave reflection and transmission, the change of the amplitude via the incident angle is related to the elasticity of medium in the subsurface two-side. In the conventional inversion with poststack datum, only the information of the reflection operator at the zero incident angles can be used. If the more robust resolutions are requested, the amplitudes of all incident angles should be used.A natural separable expression of the reflection/transmission operator is represented, which is the sum of the products of two group functions. One group function vary with phase space whereas other group function is related to elastic parameters of the medium and geological structure.By employing the natural separable expression of the reflection/transmission operator, the method of seismic wave modeling with the one-way wave equation is developed to model the primary reflected waves, it is adapt to a certain extent heterogeneous media and confirms the accuracy of AVA of the reflections when the incident angle is less than 45'. The computational efficiency of the scheme is greatly high.The natural separable expression of the reflection/transmission operator is also used to construct prestack elastic inversion algorithm. Being different from the AVO analysis and inversion in which the angle gathers formed during the prstack migration are used, the proposed algorithm construct a linear equations during the prestack migration by the separable expression of the reflection/transmission operator. The unknowns of the linear equations are related to the elasticity of the medium, so the resolutions of them provided the elastic information of the medium.The proposed method of inversion is the same as AVO inversion in , the difference between them is only the method processing the amplitude via the incident angle and computational domain.