砂砾岩油藏精细油藏描述–以准噶尔盆地西北缘八区下乌尔禾组油藏为例

Glutenite reservoir is one of the most important reservoir types in china. Because of its particularity of rock structure and pore structure, it is usually difficult in development, especially for its serious heterogeneity. On the basis of seismic, well logs, core data and production performance, the lower Wuerhe group can be divided into one second-order sequences, two third-order sequences and twenty two subsequences, corresponding to the five stages and twenty two minlayers. In addition, the fault systems are interpreted and the control action of fault systems to reservoir development is also described. The lower Wuerhe formation of 8th district belongs to fluvial-dominated fan delta sedimentation, according to the analysis of well logs, logging data and core data. It can be subdivided into two kinds of subfacies and nine kinds of microfacies. The fan delta plain subfacies mainly consist of braided channel, unconcentrated flow, mud flow and sieve deposit microfacies. The fan delta front subfacies include subaqueous distributary channel, subaqueous interdistributary channel, debris flow, subaqueous barrier and grain flow microfacies. Combined with the regional geological characteristics, the porosity model of lower Wuerhe formation is performed using core data. A permeability model based on the flow zone index is also formed according to the pore throat characteristics and flow property. Finally, the heterogeneity is analyzed. The result shows that the lower Wuerhe formation has a feature of middle-high heterogeneity, and it is controlled by material sources and sedimentary facies belt.

鄂尔多斯盆地子洲气田二叠系山西组储层特征与天然气富集规律研究

In this thesis, detailed studies on the sedimentology and petrophysical properties of reservoir rocks in the Shan#2 Member of Shanxi Formation, Zizhou gas field of Ordos Basin, are carried out, based on outcrop description, core description, wireline log interpretation and analysis of petrophysical properties. In the context of stratigraphic division scheme of the Upper Paleozoic in Ordos Basin, the Shan#2 Member is further divided into three subintervals: the Shan#23, Shan#22 and Shan#21, based on the marker beds，depositional cycles, wireline log patterns. Subaqueous deltaic-front distributary channels, distributed from the south to north, is identified，which is the main reservoir sand bodies for gases of Shan#2 Member at Zizhou gas field. Quartzose and lithic-quartzose sandstones, commonly with a high volume of cement, but a low volume of matrix, are the major reservoir rocks in the studied area. All sandstones have been evolved into the late diagenetic stage (referred to as diagenetic stage B) during the burial, experiencing compaction, cementation, replacement and dissolution, in which the compaction and cementation could have reduced the porosity, while dissolution could have improved the petrophysical properties. The pore types in the reservoirs are dominated by intergranular-solutional, intergranular-intercrystal and intercrystal-solutional porosity. According to the parameters and capillary pressure curves of test samples, five types of pore texture (I-V) are differentiated, in which types II and III pore textures displayed by low threshold pressure-wide pore throat and moderate threshold pressure-moderately wide pore throat, exist widespread. Sandstone reservoirs in the studied area are characterized by exceptionally low porosity and permeability, in which the petrophysical properties of those in Shan#23 horizon are relatively better. The petrophysical property of reservoirs was influenced both by the sedimentation and diagenesis. In general, the coarse quartzose sandstones deposited in subaqueous distributary channels show the best petrophysical property, which tends to be worse as the grain size decreases and lithic amount increases. Three types of gas reservoirs in Shan#23 horizon are classified according to petrophysical properties (porosity and permeability), which could have been influenced by the initial depositional facies, diagenesis and tectonics. On the basis of the study on the geological conditions of reservoirs in the area, it is concluded that sedimetary facies, diagenesis and tectonic actions can provide an important foundation for gas pool formation, which can also control the accumulation and distribution of gas reservoirs.

鄂尔多斯盆地陇东地区长8段砂岩输导体定量表征研究

Migration carriers act as the “Bridges” connecting source rock and traps and play important roles in petroleum migration and accumulation system. Among various types of carriers, sandstone carrier constitutes the basis of carrier system consisting of connected sandstone bodies, of sand-bodies connected with other carriers, such as faults and/or unconformities. How do we understand sandstone carrier beyond the traditional reservoirs concept? How could we characterize quantitatively this kind of carriers for petroleum migration? Such subjects are important and difficult contents in dynamic studies on hydrocarbon migration and accumulation. Sandstone carrier of Chang 8 member in Longdong area of Ordos Basin is selected as the research target in this thesis. Through conducting integrated reservoir analysis on many single wells, the correlation between single sandstone thickness and oil thickness seems good. Sketch sandstone is defined in this thesis as the principal part of carrier based on systematical analysis on lithology and sandstone thickness. Geometry connectivity of sandstone bodies was identified by the spatial superposition among them and was proved by the oil property features in oilfields. The connectivity between sandstone carriers is also hydrodynamically studied by observing and analyzed various diagenetic phenomena, especially the authigenic minerals and their forming sequence. The results were used to characterize transporting capability of sandstone carriers during the key petroleum migration periods. It was found that compaction and cementation are main causes to reduce pore space, and resolution may but not so importantly increases pore space after the occurrence of first migration. The cements of ferrocalcite and kiesel seem like the efficient index to demonstrate the hydraulic connection among sandy bodies. Diagenetic sequence and its relationship with petroleum migration phases are analyzed. Sandstone carrier of Chang 8 member was then characterized by studying their pore space and permeable properties. The results show an average porosity and permeability of Chang 8 carriers are respectively 8% and 0.50md, belongs to low porosity - low permeability reservoirs. Further, the physical properties of Chang 81 member are commonly better than those of Chang 82 member. Methods to reconstruct property of sandstone carrier during petroleum migration phase (late Jurassic) are built based on diagenetic sequence. Planal porosity, porosity and permeability of sandstone carrier in this period are statistically analyzed. One combining index - product of thickness and ancient porosity - is selected as the idea parameter to characterize sandstone carrier of late Jurassic after contrast with other parameters. Reservoirs of Chang 8 member in Longdong area are lithological reservoir controlled by sand body in which oil layers in middle part are clamped with dry layers in upper and lower parts, in a sandwich way. Based a newly proposed “migration-diagensis-remigration” model in low permeability sandstone of Chang 8 member in Longdong area, oil migration and accumulation processes during different periods are simulated with the reconstructed sandstone carriers system. Results match well with current reservoir distributions. Finally, suggestions for next favorable exploration areas are given based on all research achievements.

海底金属矿床充水条件及矿坑突水的断层效应研究

The Xinli mine area of Sanshandao mine is adjacent to the Bohai Sea and its main exploitable ore deposit occurs in the undersea rock mass. The mine is the biggest undersea gold mine of China after production. The mine area faces a latent danger of water bursting, even sudden seawater inrush. There is no mature experience in undersea mining in China so far. The vein ore deposit is located in the lower wall of a fault; its possible groundwater sources mainly include bittern, Quaternary pore water and modern seawater. To ensure the safety of undersea mining, to survey the flooding conditions of the ore deposit using proper measures and study the potential seawater inrush pattern are the key technical problems. With the Xinli mine area as a case study, the engineering geological conditions of the Xinli mine area are surveyed in situ, the regional structural pattern and rock mass framework characteristics are found out, the distribution of the structural planes are modeled by a Monte Carlo method and the connectivity coefficients of rock mass structural planes are calculated. The regional hydro-geological conditions are analyzed and the in-situ hydro-geological investigation and sampling are performed in detail, the hydrochemistry and isotopes testing and groundwater dynamic monitoring are conducted, the recharge, runoff, discharge conditions are specified and the sources of flooding are distinguished. Some indices are selected from the testing results to calculate the proportion of each source in some water discharge points and in the whole water discharge of the Xinli mine area. The temporal and spatial variations of each water source of the whole ore deposit flooding are analyzed. According to the special project conditions in the Xinli mine area, the permeability coefficient tensors of the rock mass in Xinli mine area are calculated based on a fracture geometry measurement method, in terms of the connectivity and a few hydraulic testing results, a modified synthetic permeability coefficient are calculated. The hydro-geological conceptual and mathematical model are established，the water yield of mine is predicted using Visual Modflow code. The spreading law of surrounding rock mass deformation and secondary stress are studied by numerical analysis; the intrinsic mechanism of the faults slip caused by the excavation of ore deposit is analyzed. The results show that the development of surrounding rock mass deformation and secondary stress of vein ore deposit in the lower wall of a fault, is different from that in a thick-big ore deposit. The secondary stress caused by the excavation of vein ore deposit in the lower wall of a fault, is mainly distributed in the upper wall of the fault, one surface subsidence center will occur. The influences of fault on the rock mass movement, secondary stress and hydro-geological structures are analyzed; the secondary stress is blocked by the fault and the tensile stress concentration occurs in the rock mass near the fault, the original water blocking structure is destructed and the permeable structure is reconstructed, the primary structural planes begin to expand and newborn fissures occur, so the permeability of the original permeable structure is greatly enhanced, so the water bursting will probably occur. Based on this knowledge, the possible water inrush pattern and position of the Xinli mine area are predicted. Some computer programs are developed using object-oriented design method under the development platform Visual Studio.Net. These programs include a Monte Carlo simulation procedure, a joint diagrammatizing procedure, a structural planes connectivity coefficient calculating procedure, a permeability tensor calculating procedure, a water chemical formula edit and water source fixture conditions calculating procedure. A new computer mapping algorithm of joint iso-density diagram is raised. Based on the powerful spatial data management and icon functions of Geographic Information System, the pit water discharge dynamic monitoring data management information systems are established with ArcView.

土石混合体细观结构力学及其边坡稳定性研究

Soil-rock mixture (S-RM) refers to one extremely uneven loose rock and soil materials system with certain stone content. Its formation has started since Quaternary and it is composed of block stone, fine grained soil and pore with certain project scale and high strength. S-RM has extensive distribution in nature, especially in southwest China where the geotectonic background is complicated, the fracture activity is developed and the geomorphological characteristics of high mountain and steep gorge area are protuberant. This kind of complicated geologic body has developed wider in these areas. S-RM has obvious difference with the general soil or rock (rock mass) in physical and mechanical properties because its two components-“soil” and “rock-block” has extreme differences in physical and mechanical properties. The proposition of S-RM and its deep research are needed in the modern engineering construction. It is also the necessity in the modern development of rock and soil mechanics. The dissertation starts from the meso-structural characteristics of soil-rock and takes a systematic research on its meso-structural mechanics, deformation and failure mechanism and the stability of S-RM slope. In summary, it achieves the following innovative results and conclusions. There are various views on the conception of S-RM and its classification system. Based on the large number of field tests, the dissertation makes the conception and classification of S-RM more systematic. It systematically proposed the conception of meso-structural mechanics of S-RM. Thus the dissertation has laid a foundation for its deep study. With the fast development of the computer technology and digital image processing theory, digital image processing technology has been successfully applied in many fields and provided reliable technology support for the quantitative description of the structural characteristics of S-RM. Based on the digital image processing technology, the dissertation systematically proposes and developed the quantitative analysis method and quantitative index for the meso-structure of S-RM. The results indicate that the meso-structure such as its internal soil-rock granularity composition, the soil-rock shape and the orientability has obvious self-organization in the macro statistical level. The dissertation makes a systematic research on the physical mechanical properties, deformation and failure mechanism of S-RM based on large field test. It proposes the field test for the underwater S-RM and deduces the 3D data analysis method of in-situ horizontal push-shear test. The result indicates that S-RM has significant phenomenon of shear dilatancy in the shearing process, and its dilatancy will be more obvious with the increased proportion of rock or the decreased confining pressure. The proportion of rock has great effect on the strength of S-RM and rock-block, especially the spatial position of particles with comparatively big size has great effect on the shape and spatial position of the sample shear zone. The dissertation makes some improvements in the single ring infiltration test equipment and its application on the permeability of S-RM. The results indicate that the increasing of rock-block would make it more difficult for the soil to fill in the vacuity between the rock-block and the proportion would increase which would result in the increased permeability coefficient. The dissertation builds the real meso-structural model of S-RM based on the digital image processing technology. By using geometric reconstruction technology, it transfers the structural mode represented by Binary image into CAD format, which makes it possible to introduce the present finite element analysis software to take research on numerical experimental investigation. It systematically realizes leaping research from the image，geometric mode, to meso-structural mechanics numerical experiment. By using this method, the dissertation takes large scale numerical direct-shear test on the section of S-RM. From the mesoscopic perspective, it reveals three extended modes about the shear failure plane of S-RM. Based on the real meso-structural model and by using the numerical simulation test, the character and mechanics of seepage failure of S-RM are studied. At the same time, it builds the real structural mode of the slope based on the analysis about the slope crosssection of S-RM. By using the strength reduction method, it takes the research on the stability of S-RM and gets great achievements. The three dimensional geometric reconstruction technology of rock block is proposed, which provides technical support for the reconstruction of the 3D meso-structural model of S-RM. For the first time, the dissertation builds the stochastic structure model of two-dimensional and three-dimensional polygons or polyhedron based on the stochastic simulation technique of monte carlo method. It breaks the traditional research which restricted to the random generation method of regular polygon and develops the relevant software system (R-SRM2D/3D) which has great effect on meso-structural mechanics of S-RM. Based on the R-SRM software system which randomly generates the meso-structural mode of S-RM according to the different meso-structural characteristics, the dissertation takes a series of research on numerical test of dual axis and real three-axis, systematically analyses the meso destroy system, the effects of meso-structural characteristics such as on the stone content, size composition and block directionality on the macro mechanical behavior and macro-permeability. Then it proposes the expression of the upper and lower limit for the macro-permeability coefficient of the inhomogeneous geomaterials, such as S-RM. By using the strength reduction FEM, the dissertation takes the research on the stability of the slope structural mode of the randomly formed S-RM. The results indicate that generally, the stability coefficient of S-RM slope increases with the increasing of stone content; on the condition of the same stone content, the stability coefficient of slope will be different with different size composition and the space position of large block at the internal slop has great effect on the stability. It suggests that meso-structural characteristics, especially the space position of large block should be considered when analyzing the stability of this kind of slope and strengthening design. Taking Xiazanri S-RM slope as an example, the dissertation proposes the fine modeling of complicated geologic body based on reverse engineering and the generation method of FLAC3D mode. It resolves the bottleneck problem about building the fine structural mode of three-dimensional geological body. By using FLAC3D, the dissertation takes research on the seepage field and the displacement field of Xiazanri S-RM slope in the process of reservoir water level rising and decreasing. By using strength reduction method, it analyses the three-dimension stability in the process of reservoir water level rising and decreasing. The results indicate that the slope stability firstly show downward trend in the process of reservoir water level rising and then rebound to increase; the sudden drawdown of reservoir water level has great effect on the slope stability and this effect will increase with the sudden drawdown amplitude rising. Based on the result of the rock block size analysis of S-RM, and using R-SRM2D the stochastic structure model of Xiazanri S-RM slope is built. By using strength reduction method, the stability of the stochastic structure model is analysis, the results shows that the stability factor increases significantly after considering the block.

鄂尔多斯盆地志丹地区延长组储层综合评价

Sedimentary provenance direction，sedimentary facies，reservoir geological characteristic，pore structure； physical property characteristic，reservoir classification and evaluation ，forthermore，favorable area area are forecasted of Yanchang formation in ZhiDan region， by mainly study on the data of field section observation ，core observation, well logging explaination and routine microscope slice identification，scanning Electron Microscope and reservoir analysis of lithology and physical property , Under the guidance of such advanced theories and methods as sedimentology，reservoir sedimentology，lithological oil pool and so on，in the thesis. The following fruits are mainly achieved in this paper： Yanchang formation stratum is divided and correlated in this entire region， and the characteristic of oil layer unit is detailed discussed , respectively. According to main marker bed and supplementary ones.and research result shows that the source of provenance direction during Yanchang Formation mianly is north-east. Delta and lake are mainly developed in study area ，sub-facies and micro-facies are divided，distribution of sedimentary micro-facies in plane and palaeogeographic evolution are described，and gentle slope type- shallow water delta depositional model is established. Fine-grain arkose sandstone is the main reservoir，and which have experienced such different degree diagenesis as compaction, cementation, replacement and dissolution, and in which compaction and cementation are mainly factors to reduce sandstone physical property and dissolution effectively improved physical property during burial diagenesis procedure. All reservoirs of Yanchang Formation have entered A period of late diagenetic stage according to scheme of diagenesis period division . Intergranular porosity，dissolution porosity，fissure porosity are main pore types. And porosity structure are analyse by mercury penetration capillary pressure parameter，fine-shortness type and fine- length throat type are mainly developed. as a whole，the reservoir， with the characteristic of porosity and permeability altering apparently，strong inhomogeneity ， is a medium- porosity and medium permeability one. In plane，higher- porosity and higher-permeability are corresponded well with distributary channel area， physical property and inhomogeneity are affected by both deposition and diagenesis，and distributary channel and underwater distributary channel are favorable facies . According to such characteristic as lithology，physical property，pore structure ，diagenesis and sandstone distribution, the sandy reservoir can be classified 4 types, and the main sandy in every oil layer unit are evaluated according to the standard. The analysis result of petroleum concentration rule shows that Yanchang Formation are with not only favourable oil source rock，reservoir，covering combination ，but also good entrapment condition in study area. Lithology and structure-lithology oil pool are mainly developed ，based on condition of favorable reservoir developments，accounting for deliverability and sandstone superface elevation，zone of profitabilitis are forecasted.

鄂尔多斯盆地三叠系延长组特低渗透储层分类评价

On the basis of the character of sedimentation and reservoir researching as well as diagenesis, using conventional and update testing measures, classificati-on and evaluation of the tesla low permeability reservoir in Ordos Basin is pr-esented. From Chang 8 to Chang 4+5 oil formations, four facies developed, includi-ng alluvial fan facies, delta facies, lake facies as well as density current. They were controlled by the northeastern, the southwest, the southern and the northwestern provenances. Distributary channel underwater and mouth bar of delta fr-ont are the main reservoirs. Detrital component has the different character in s-outh and in north. Sedimentary system in the northeastern part has more felds-par and less quartz. Sedimentary system in the southern part has more quartz and less feldspar. Because of sedimentation and diagenesis, the oil formations in region of interest formed the different features of pore array of the tesla l-ow permeability reservoirs. After researching, it is found that the active porosity and the main throat radius of Chang 4+5 are the highest, and they are positive correlation with per-meability. The exponent of flowing interval falls in the sortorder: Chang 8, Chang 4+5, Chang 6, Chang 7. Using clustering procedure and quaternion, the reservoirs of Yanchang for-mation in Ordos Basin are divided into five types. Ⅰ-good reservoirs and Ⅱ-appreciably good reservoirs occur in distributary channel and mouth bar. Ⅲ-poor reservoirs and Ⅳ-poorer reservoirs exist in natural levee, crevasse splay under-water and turbidity fan. It is forecasted that the oil area in Ⅰ-good reservoirs is about 4336.68 square kilometers, and the oil area in Ⅱ-appreciably good reservoirs is 28013.28 square kilometers or so, and the oil area in Ⅲ-poor rese-rvoirs is 28538.05 square kilometers more or less.

松辽盆地南部大情字井地区CO2埋存体地质格架及储层建模

Increasing attentions have been paid to the subsurface geological storage for CO2 in view of the huge storage capacity of subsurface reservoirs. The basic requirement for subsurface CO2 storage is that the CO2 should be sequestrated as supercritical fluids (physical trapping), which may also interact with ambient reservoir rocks and formation waters, forming new minerals (chemical trapping). In order to the effective, durable and safe storage for CO2, enough storage space and stable sealing caprock with strong sealing capacity are necessitated, in an appropriate geological framework. Up till now, hydrocarbon reservoirs are to the most valid and appropriate CO2 storage container, which is well proven as the favorable compartment with huge storage capacity and sealing condition. The thesis focuses on two principal issues related to the storage and sealing capacity of storage compartment for the Qingshankou and Yaojia formations in the Daqingzijing block, Southern Songliao Basin, which was selected as the pilot well site for CO2-EOR storage. In the operation area, three facies, including deltaic plain, deltaic front and subdeep-deep lake facies associations, are recognized, in which 11 subfacies such as subaqueous distributary channel, river- mouth bar, interdistributary bay, sheet sandbody, crevasse splay and overflooding plain are further identified. These subfacies are the basic genetic units in the reservoir and sealing rocks. These facies further comprise the retrogradational and progradational depositional cycles, which were formed base- level rise and fall, respectively. During the regressive or lake lowstand stage, various sands including some turbidites and fans occurred mostly at the bottom of the hinged slope. During the progradation stage, these sands became smaller in size and episodically stepped backwards upon the slope, with greatly expanded and deeped lake. However, most of Cretaceous strata in the study area, localized in the basin centre under this stage, are mainly composed of grey or grizzly siltstones and grey or dark grey mudstones intercalated with minor fine sandstones and purple mudstones. On the base of borehole and core data, these siltstones are widespread, thin from 10 to 50 m thick, good grain sorting, and have relative mature sedimentary structures with graded bedding and cross- lamination or crossbeds such as ripples, which reflect strong hydrodynamic causes. Due to late diagenesis, pores are not widespread in the reservoirs, especially the first member of Qingshankou formation. There are two types of pores: primary pore and secondary cores. The primary pores include intergranular pores and micropores, and the secondary pores include emposieus and fracture pores. Throat channels related to pores is also small and the radius of throat in the first, second and third member of Qingshankou formation is only 0.757 μm, 0.802 μm and 0.631 μm respectively. In addition, based on analyzing the probability plot according to frequency of occurrence of porosity and permeability, they appear single- peaked distribution, which reflects strong hetero- geneity. All these facts indicate that the conditions of physical property of reservoirs are not better. One reason may be provided to interpret this question is that physical property of reservoirs in the study area is strong controlled by the depositional microfacies. From the statistics, the average porosity and permeability of microfacies such as subaqueous distributary channel, channel mouth bar, turbidites, is more than 9 percent and 1md respectively. On the contrary, the average porosity and permeability of microfacies including sand sheet, flagstone and crevasse splay are less than 9 percent and 0.2md respectively. Basically, different hydrodynamic environment under different microfacies can decide different physical property. According to the reservoir models of the first member of Qingshankou formation in the No. well Hei47 block, the character of sedimentary according to the facies models is accord to regional disposition evolution. Meantime, the parameter models of physical property of reservoir indicate that low porosity and low permeability reservoirs widespread widely in the study area, but the sand reservoirs located in the channels are better than other places and they are the main sand reservoirs. The distribution and sealing ability of fault- fractures and caprock are the key aspects to evaluate the stable conditions of compartments to store CO2 in the study area. Based on the core observation, the fractures widespread in the study area, especially around the wells, and most of them are located in the first and second member of Qingshankou formation, almost very few in the third member of Qingshankou formation and Yaojia formation instead. In addition, analyzing the sealing ability of eleven faults in the three-dimensional area in the study area demonstrates that most of faults have strong sealing ability, especially in the No. well Hei56 and Qing90-27. To some extent, the sealing ability of faults in the No. well Hei49, Qing4-6 and Qing84-29 are worse than others. Besides, the deposition environment of most of formations in the study area belongs to moderately deep and deep lake facies, which undoubtedly take advantage to caprocks composed of mudstones widespread and large scale under this deposition environment. In the study area, these mudstones distribute widely in the third member of Qingshankou formation, Yaojia and Nenjiang formation. The effective thickness of mudstone is nearly ~550m on an average with few or simple faults and fractures. In addition, there are many reservoir beds with widely- developed insulated interbeds consist of mudstones or silty mudstone, which can be the valid barrier to CO2 upper movement or leakage through diffusion, dispersion and convection. Above all, the closed thick mud caprock with underdeveloped fractures and reservoir beds can be taken regard as the favorable caprocks to provide stable conditions to avoid CO2 leakage.

蒙古国宗巴音凹陷沉积相研究及优质储层预测

Abstract:Little fundamental work on petroleum exploration and production of Zuunbayan Subbasin, Mongolia has been done before because of the backward economy and petroleum industry techniques in this country, which also results in our little knowledge of reservoir characteristics of this area. This paper focused on the sedimentary system, sedimentary facies, reservoir characteristics and their genesis distribution of Zuunbayan subbasin with various drilling, well logging, seismic, coring and outcrop data, aiming at providing significant guidances for the petroleum exploration and production of Zuunbayan area. Therefore, several conclusions have been achieved as follows: ①In Zuunbayan Subbasin, there are two chief source areas with Tarkhyata and Totoshan Uplifts in the southeast and Saykhandulaan Uplift in the west, respectively, while two subsidiary ones in the northeast and southwest of this subbasin. The sedimentary system of alluvial fan-fan delta is formed in the southeast highland, meanwhile braided river-braided river delta develops in the western ramp region and fan delta in the southern palaeohigh. ②There are middle to high permeability reservoirs in the upper Zuunbayan Formation and the upper member of lower Zuunbayan Formation meanwhile low-porosity and permeability to ultra-low permeability ones in Tsagaan Tsav Formation and the middle and lower members of lower Zuunbayan Formation. Combing with sedimentary facies belt, oil sources conditions and tectonic settings, favorable reservoir belts have been proved to be existing in the fan delta front reservoirs of lower Zuunbayan – Tsagaan Tsav Formation in the central uplift faulted zone as well as the braided river front ones of lower Zuunbayan-Tsagaan Tsav Formation in Zuunbayan nose anticlinal structural belts. ③The reservoir lithologic composition is complex and also related to volcanic activities. Generally, the types of lithologic composition in Zuunbayan Subbasin are chiefly feldspathic litharenites with low compositional maturity and high-middle textural maturity. The rock constituents from upper Cretaceous to lower Zuunbayan Formation are mainly metamorphic rocks including cleaving stone, phyllite, quartzite and schist while volcanic tuffs and acidic extrusive rocks are the secondary; and in the Tsagaan Tsav Formation are mainly volcanic tuffs with subsidiary cleaving stone, phyllite, quartzite and schist. ④In this paper, high-quality reservoirs in the upper member of lower Zuunbayan Formation have been discovered in the drilled high production wells of favorable reservoir facies through sedimentary system and sedimentary facies research, which benefits the prospect and also will bring a new life for petroleum exploration and production of Zuunbayan Subbasin. Key words: sedimentary system, sedimentary facies, superior quality reservoir, Zuunbayan Subbasin, lower Zuunbayan Formation

二氧化碳在孔隙咸水层中多相流动数值模拟研究

There has been a growing concern about the use of fossil fuels and its adverse effects on the atmospheric greenhouse and ecological environment. A reduction in the release rate of CO2 into the atmosphere poses a major challenge to the land ecology of China. The most promising way of achieving CO2 reduction is to dispose of CO2 in deep saline aquifers. Deep aquifers have a large potential for CO2 sequestration in geological medium in terms of volume and duration. Through the numerical simulation of multiphase flow in a porous media, the transformation and motion of CO2 in saline aquifers has been implemented under various temperature and hydrostatic pressure conditions, which plays an important role to the assessment of the reliability and safety of CO2 geological storage. As expected, the calculated results can provide meaningful and scientific information for management purposes. The key problem to the numerical simulation of multiphase flow in a porous media is to accurately capture the mass interface and to deal with the geological heterogeneity. In this study, the updated CE/SE (Space and time conservation element and solution element) method has been proposed, and the Hybrid Particle Level Set method (HPLS) has extended for multiphase flows in porous medium, which can accurately trace the transformation of the mass interface. The benchmark problems have been applied to evaluate and validate the proposed method. In this study, the reliability of CO2 storage in saline aquifers in Daqingzi oil field in Sunlong basin has been discussed. The simulation code developed in this study takes into account the state for CO2 covering the triple point temperature and pressure to the supercritical region. The geological heterogeneity has been implemented, using the well known geostatistical model (GSLIB) on the base of the hard data. The 2D and 3D model have been set up to simulate the CO2 multiphase flow in the porous saline aquifer, applying the CE/SE method and the HPLS method .The main contents and results are summarized as followings. (1) The 2D CE/SE method with first and second –order accuracy has been extended to simulate the multiphase flow in porous medium, which takes into account the contribution of source and sink in the momentum equation. The 3D CE/SE method with the first accuracy has been deduced. The accuracy and efficiency of the proposed CE/SE method have been investigated, using the benchmark problems. (2) The hybrid particle level set method has been made appropriate and extended for capturing the mass interface of multiphase flows in porous media, and the numerical method for level set function calculated has been formulated. (3) The closed equations for multiphase flow in porous medium has been developed, adept to both the Darcy flow and non-Darcy flow, getting over the limitation of Reynolds number to the calculation. It is found that Darcy number has a decisive influence on pressure as well as velocity given the Darcy number. (4) The new Euler scheme for numerical simulations of multiphase flows in porous medium has been proposed, which is efficient and can accurately capture the mass interface. The artificial compressibility method has been used to couple the velocities and pressure. It is found that the Darcy number has determinant effects on the numerical convergence and stability. In terms of the different Darcy numbers, the coefficient of artificial compressibility and the time step have been obtained. (5) The time scale of the critical instability for critical CO2 in the saline aquifer has been found, which is comparable with that of completely CO2 dissolved saline aquifer. (6) The concept model for CO2 multiphase flows in the saline aquifer has been configured, based on the temperature, pressure, porosity as well as permeability of the field site .Numerical simulation of CO2 hydrodynamic trapping in saline aquifers has been performed, applying the proposed CE/SE method. The state for CO2 has been employed to take into account realistic reservoir conditions for CO2 geological sequestration. The geological heterogeneity has been sufficiently treated , using the geostatistical model. (7) It is found that the Rayleigh-Taylor instability phenomenon, which is associated with the penetration of saline fluid into CO2 fluid in the direction of gravity, has been observed in CO2 multiphase flows in the saline aquifer. Development of a mushroom-type spike is a strong indication of the formation of Kelvin-Helmholtz instability due to the developed short wavelength perturbations present along the interface and parallel to the bulk flow. Additional key findings: the geological heterogeneity can distort the flow convection. The ascending of CO2 can induce the persistent flow cycling effects. The results show that boundary conditions of the field site have determinant effects on the transformation and motion of CO2 in saline aquifers. It is confirmed that the proposed method and numerical model has the reliability to simulate the process of the hydrodynamic trapping, which is the controlling mechanism for the initial period of CO2 storage at time scale of 100 years.

软土亚塑性模型与一维非线性固结理论研究

With the great development of Tianjing New Coastal District economy, people need more land to build and live. Land subsidence, which is caused by its special engineering geological conditions, has restricted the further development in the district. Soft soil consolidation is main factor of land subsidence ;thus , on the basis of consolidation theory, the paper make further study on soft soils one-dimension nonlinear consolidation which contains two parts:(1) the nonlinear consolidation of permeability coefficient and compressibility coefficient changing with time and depth, which means real one-dimension nonlinear consolidation;(2) the non-homogeneous consolidation of permeability coefficient and compressibility coefficient only changing with depth. Firstly, nonlinear characteristics of soft soils are elaborated. Hypoplastic theory is introduced to establish a modified soft soils nonlinear constitutive model; the nonlinear governing equation of compressibility coefficient is built, and the nonlinear characteristics of compressibility coefficient are analyzed. Secondly, Considering Load Fluctuation and soil thickness changing ,the consolidation characteristics of single layer is discussed in the paper; meanwhile, on the basis of the Davis and Raymond’s hypothesis and single layer nonlinear consolidation equation, the doubled-layer one-dimension nonlinear consolidation equation is also derived. The solution of the equation is obtained by analytical method, and the consolidation characteristics of doubled-layer soft soil nonlinear theory is also analyzed. Finally, based on assumption that permeability coefficient and compressibility coefficient is varying along depth, single layer soil one-dimension non-homogeneous consolidation differential equation is derived; and the approximate solution is obtained. Furthermore, the single layer non-homogeneous consolidation is extended to double layer non-homogeneous consolidation theory. By using parabolic differential scheme, the matrix equation is established; and the solution of the matrix equation is obtained by chase method. Consolidation characteristics of soil soft single (double) layer non-homogeneous consolidation theory and Terzaghi’s theory are also discussed.

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

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.

基于软粘土微观特性的非线性固结理论研究

The land subsidence of soft clay is including natural and man-made content, which leads to the research on the mechanism of land subsidence constituted by two different aspects, which are studied by geological engineers and geologist. The main major research is focused on the effects of engineering. The land subsidence engineering of soil mechanics is caused by the consolidation and compression of soft clay, the content of which is including the micro-structural characteristics, the stress - strain constitutive relation, porous law, and consolidation theory. In this paper, it is discussed the nonlinear consolidation and compression theory of soft clay. The main studies and conclusions of this thesis are as follows. （1）The micro-structure and its stability are closely related to the engineering characters of soft clay. The stiffness and force connection status of micro-structure plays a controlling influence to its stability. （2）Under saturated state, clay particles remain in a non-full contact or non-contact status, so it is needed to modify the Terzaghi effective stress principle. With the discharge of pore water, the effective stress is increasing, and part of weakly bound-water begins flow, while the porosity and permeability are became lower. （3）It exist non-linear flow in soft clay, which is caused by the shear flow situation of weakly bounded-water. In this case, permeability coefficient is a nonlinear function of hydraulic gradient. （4）In the initial consolidation stage of soft clay in the initial stage, the porous flow is mainly caused by the excretion of free water. With the decrease of free water content, combined bonded-water start to supply free water. At the later stage of consolidation, the flow of fluid is mainly consisted by weakly bounded-water. The exchange between bonded-water and free water is played a role, which slows down the consolidation process.

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

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

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

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