胜坨油田沙二段油藏流场剩余油形成机理和分布模式

By applying multi-discipline theory and methods comprehensively and with full use of computer, the paper deeps into studying changing rule and control factor of fluid field of ES2 Shengtuo oil field during waterflood development, physical and chemical function, and stress. Matrix field, network field, fluid field, stress field and physical chemistry field and fluid model for dynamic function were established. Macroscopic and microscopic genesis mechanics, distribution rule and control factor of remaining oil were revealed. Remaining oil and emulate model were established. Macroscopic and microscopic distribution rule of mover remaining oil were predicted, several results were achieved as following: The distribution of remaining oil was controlled by micro-structure. At the same development stage, remaining oil saturation of the wells located in higher position of micro-structure is higher than the average saturation in the same layer. The water content ratio has same law. It is the enrichment district that the high position of micro-structure controlled by seal faults. The remaining oil distribution was affected by sedimentary micro-facies, micro-structure, fault sealing, reservoir heterogeneity and affusion-oil extraction. On the plane, the zone owning higher saturation of remaining oil is the area that at the edge miacro-facies and sand-body distribution discontinuously; on the section, the content of waterflood of the upper or middle-upper oil layer of positive rhythm and positive comprehensive rhythm is lower, middle and weak waterflood is main, remaining oil is in enrichment relatively. The remaining oil is relative enrichment at the zones of well network of affusion and oil extraction not affected. 4D dynamic model of reservoir of Es2in Shengtuo oil field was established. Macroscopic and microscopic forming mechanics, distribution rule and control factor were revealed. The emulate model of dynamic function of Shengtuo oil field was established, the space distribution of remaining oil were predicted. Reservoir flow field, matrix field, network field, seep field, physical and chemical field, stress field and fluid field models were established. Reservoir flow field character and distribution were revealed. An improvement of the development geology theory in continental fault depression continental basin was brought on.

电阻率测井反演方法应用研究

Formation resistivity is one of the most important parameters to be evaluated in the evaluation of reservoir. In order to acquire the true value of virginal formation, various types of resistivity logging tools have been developed. However, with the increment of the proved reserves, the thickness of interest pay zone is becoming thinner and thinner, especially in the terrestrial deposit oilfield, so that electrical logging tools, limited by the contradictory requirements of resolution and investigation depth of this kinds of tools, can not provide the true value of the formation resistivity. Therefore, resitivity inversion techniques have been popular in the determination of true formation resistivity based on the improving logging data from new tools. In geophysical inverse problems, non-unique solution is inevitable due to the noisy data and deficient measurement information. I address this problem in my dissertation from three aspects, data acquisition, data processing/inversion and applications of the results/ uncertainty evaluation of the non-unique solution. Some other problems in the traditional inversion methods such as slowness speed of the convergence and the initial-correlation results. Firstly, I deal with the uncertainties in the data to be processed. The combination of micro-spherically focused log (MSFL) and dual laterolog(DLL) is the standard program to determine formation resistivity. During the inversion, the readings of MSFL are regarded as the resistivity of invasion zone of the formation after being corrected. However, the errors can be as large as 30 percent due to mud cake influence even if the rugose borehole effects on the readings of MSFL can be ignored. Furthermore, there still are argues about whether the two logs can be quantitatively used to determine formation resisitivities due to the different measurement principles. Thus, anew type of laterolog tool is designed theoretically. The new tool can provide three curves with different investigation depths and the nearly same resolution. The resolution is about 0.4meter. Secondly, because the popular iterative inversion method based on the least-square estimation can not solve problems more than two parameters simultaneously and the new laterolog logging tool is not applied to practice, my work is focused on two parameters inversion (radius of the invasion and the resistivty of virgin information ) of traditional dual laterolog logging data. An unequal weighted damp factors- revised method is developed to instead of the parameter-revised techniques used in the traditional inversion method. In this new method, the parameter is revised not only dependency on the damp its self but also dependency on the difference between the measurement data and the fitting data in different layers. At least 2 iterative numbers are reduced than the older method, the computation cost of inversion is reduced. The damp least-squares inversion method is the realization of Tikhonov's tradeoff theory on the smooth solution and stability of inversion process. This method is realized through linearity of non-linear inversion problem which must lead to the dependency of solution on the initial value of parameters. Thus, severe debates on efficiency of this kinds of methods are getting popular with the developments of non-linear processing methods. The artificial neural net method is proposed in this dissertation. The database of tool's response to formation parameters is built through the modeling of the laterolog tool and then is used to training the neural nets. A unit model is put forward to simplify the dada space and an additional physical limitation is applied to optimize the net after the cross-validation method is done. Results show that the neural net inversion method could replace the traditional inversion method in a single formation and can be used a method to determine the initial value of the traditional method. No matter what method is developed, the non-uniqueness and uncertainties of the solution could be inevitable. Thus, it is wise to evaluate the non-uniqueness and uncertainties of the solution in the application of inversion results. Bayes theorem provides a way to solve such problems. This method is illustrately discussed in a single formation and achieve plausible results. In the end, the traditional least squares inversion method is used to process raw logging data, the calculated oil saturation increased 20 percent than that not be proceed compared to core analysis.

基于MRF模型的油气储层属性随机建模方法研究

Stochastic reservoir modeling is a technique used in reservoir describing. Through this technique, multiple data sources with different scales can be integrated into the reservoir model and its uncertainty can be conveyed to researchers and supervisors. Stochastic reservoir modeling, for its digital models, its changeable scales, its honoring known information and data and its conveying uncertainty in models, provides a mathematical framework or platform for researchers to integrate multiple data sources and information with different scales into their prediction models. As a fresher method, stochastic reservoir modeling is on the upswing. Based on related works, this paper, starting with Markov property in reservoir, illustrates how to constitute spatial models for catalogued variables and continuum variables by use of Markov random fields. In order to explore reservoir properties, researchers should study the properties of rocks embedded in reservoirs. Apart from methods used in laboratories, geophysical means and subsequent interpretations may be the main sources for information and data used in petroleum exploration and exploitation. How to build a model for flow simulations based on incomplete information is to predict the spatial distributions of different reservoir variables. Considering data source, digital extent and methods, reservoir modeling can be catalogued into four sorts: reservoir sedimentology based method, reservoir seismic prediction, kriging and stochastic reservoir modeling. The application of Markov chain models in the analogue of sedimentary strata is introduced in the third of the paper. The concept of Markov chain model, N-step transition probability matrix, stationary distribution, the estimation of transition probability matrix, the testing of Markov property, 2 means for organizing sections-method based on equal intervals and based on rock facies, embedded Markov matrix, semi-Markov chain model, hidden Markov chain model, etc, are presented in this part. Based on 1-D Markov chain model, conditional 1-D Markov chain model is discussed in the fourth part. By extending 1-D Markov chain model to 2-D, 3-D situations, conditional 2-D, 3-D Markov chain models are presented. This part also discusses the estimation of vertical transition probability, lateral transition probability and the initialization of the top boundary. Corresponding digital models are used to specify, or testify related discussions. The fifth part, based on the fourth part and the application of MRF in image analysis, discusses MRF based method to simulate the spatial distribution of catalogued reservoir variables. In the part, the probability of a special catalogued variable mass, the definition of energy function for catalogued variable mass as a Markov random field, Strauss model, estimation of components in energy function are presented. Corresponding digital models are used to specify, or testify, related discussions. As for the simulation of the spatial distribution of continuum reservoir variables, the sixth part mainly explores 2 methods. The first is pure GMRF based method. Related contents include GMRF model and its neighborhood, parameters estimation, and MCMC iteration method. A digital example illustrates the corresponding method. The second is two-stage models method. Based on the results of catalogued variables distribution simulation, this method, taking GMRF as the prior distribution for continuum variables, taking the relationship between catalogued variables such as rock facies, continuum variables such as porosity, permeability, fluid saturation, can bring a series of stochastic images for the spatial distribution of continuum variables. Integrating multiple data sources into the reservoir model is one of the merits of stochastic reservoir modeling. After discussing how to model spatial distributions of catalogued reservoir variables, continuum reservoir variables, the paper explores how to combine conceptual depositional models, well logs, cores, seismic attributes production history.

准噶尔盆地低阻油藏实验研究及储层评价方法

Low resistivity reservoir is a special reservoir which is different from normal reservoir in identification and evaluation.Through core experiment and analysis, the achievement of which resistivity is resulted from clay additive electric conductivities and high bound water saturation in Junggar basin is gained. For accurately evaluating low resistivity, a good many of experiment have been completed, such as resistivity index and formation factor in hi^jher temperature and higher pressure, semi-permeability board, cation exchange, bound water, NMR (nucleus magnetism response), non-Nad water in different temperature and salinity, the experiments result show that lower resistivity has complex relation with these electric-parameters and chloric ion content in non-NaCl water.Based on comprehensive interpretation of NMR and normal resistivity data, the volume of moved water, bound water, moved oil and residual oil in the strata can be determined quantitatively and which have significant influence on reservoir recognition and perforation optimized.Experiment data (SEM mold, thin section, X ray diffraction, mercury penetration) can be used to analysis low resistivity forming and the relation between low resistivity and pore texture, to set up relation between porosity, permeability and petrophysical property. The reservoir was sorted, evaluated and described. The oil bedding in southern margin of Junggar basin is low porosity, low resistivity reservoir.Based on invasion theory of electric well-logging, modelling and inversion of resistivity well-logging are accomplished. For enhancing low resistivity resulted from higher bound water saturation and cation exchange, invasion period, invasion radius, the relation between fluid distribution in pore and response of laterolog logging have been studied. Virgin zone resistivity, invasion zone resistivity and invasion radius were inversed and which enhanced evaluation accuracy of reservoir. The method was used to process well-logging data in Luliang oilfield and southern margin in Junggar basin, and reservoir resistivity was enhanced effectively, appropriate oil saturation gained and it has better effect on oil exploration.

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

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.

大别（超）包产到户铁岩和阜新火山岩铂族元素斌存规律研究

The platinum-group elements (PGE), including Os, Ir, Ru, Rh, Pt and Pd, axe strongly siderophile and chalcophile. On the basis of melting temperature, the PGE may be divided into two groups: the Ir group (IPGE, >2000°C) consisting of Os, Ir and Ru, and the Pd group (PPGE, <20GO°C) consisting of Rh, Pt and Pd. Because of their unique geochemical properties, PGE provide critical information on global-scale differentiation processes, such as core-mantle segregation, late accretionary history, and core-mantle exchange. In addition, they may be used to identify magma source regions and unravel complex petrogenetic processes including partial melting, melt percolation and metasomatism in the mantle, magma mixing and crustal contamination in magma chambers and melt crystallization.Compared with other rocks, (ultra)mafic rocks have lower REE content but higher PGE content, so PGE have advantages in studying the petrogeneses and evolution of them. In this study, we selected (ultra)mafic rocks collected in Dabie Orogen and volcanic rocks from Fuxin Region. Based on the distribution and behaviour of platinum-group elements, combined with other elements, we speculate the magma evolution and source mantle of these (ultra)mafic rocks and volcanic rocks.Many (ultra)mafic rocks are widely distributed in Dabie Region. According to their deformation and metamorphism, we classed them into three types. One is intrusive (ultra)mafic rocks, which are generally undeformed and show no or little sign of metamorphism, such as (ultra)mafic intrusions in Shacun, zhujiapu, Banzhufan, qingshan, Xiaohekou, Jiaoziyan, Renjiawan and Daoshichong. The other one is ultrahigh pressure metamorphic (ultra)mafic rocks, some of them appeared as eelogites, such as complex in Bixiling and adjacent Maowu. Another one is intense deformed and metamorphic, termed as tectonic slice, alpine-type (ultra)mafic rocks. The most representative is Raobazhai and Dahuapin. However, there are many controversies about the formation of those (ultra)mafic rocks. Here, we select typical rocks of the three types. The PGE were determined by inductively coupled plasma mass spectrometry (ICP-MS) ater NiS fire-assay and tellurium co-precipitation.The PGE tracing shows that three components are needed in the source of the cretaceous (uitra)mafic intrusions. They could be old enriched sub-continental lithospheric mantle, lower crust and depleted asthenospheric mantle. The pattern of PGE also shows the primitive magma of these intrusions underwent S saturation. According to palladium, we can conclude that the mantle enrich in PGE. Distribution of PGE in Bixiiing and Maowu (ultra)mafic rocks display they are products of magmas fractional crystallization. The (ultra)mafic rocks in Bixiiing and Maowu are controlled by various magmatic processes and the source mantle is depleted in PGE. Of interest is that the mantle produced UHP (ultra)mafic rocks are PGE-depleted, whereas the mantle of cretaceous (ultra)mafic intrusions are enrich in PGE. This couldindicate that the mantle change from PGE-enriched to PGE-depleted during120-OOMa, which in accord with the time of tectonic system change in the East China. At the same time, (ultra)mafic intrusions in cretaceous took information of deep mantle, which means the processes in deep mantle arose structural movement in the crust The character of PGE in alpine-type (ultra)mafic rocks declared that the rocks had experienced two types of metasomatic processes - hydrous melt derived from slab and silicate melt. In addition, we analyze the platinum-group elements in volcanic rocks on the northern margin of the North China Craton, Fuxin. The volcanic rocks characterized by negative anomalies of platinum. This indicates that platinum alloys, which may host some Pt resided in the mantle. The PGE patterns also show that Jianguo alkali basalts derived from asthenospheric mantle source, but wulahada high-Mg andesites derived from lithospheric mantle.

CO adsorption and correlation between CO surface coverage and activity/selectivity of preferential CO oxidation over supported Ag catalyst: an in situ FTIR study

In situ IR measurements for CO adsorption and preferential CO oxidation in H-2-rich gases over Ag/SiO2 catalysts are presented in this paper. CO adsorbed on the Ag/SiO2 pretreated with oxygen shows a band centered around 2169 cm(-1), which is assigned to CO linearly bonded to Ag+ sites. The amount of adsorbed CO on the silver particles ( manifested by an IR band at 2169 cm(-1)) depends strongly on the CO partial pressure and the temperature. The steady-state coverage on the Ag surface is shown to be significantly below saturation, and the oxidation of CO with surface oxygen species is probably via a non-competitive Langmuir Hinshelwood mechanism on the silver catalyst which occurs in the high-rate branch on a surface covered with CO below saturation. A low reactant concentration on the Ag surface indicates that the reaction order with respect to Pco is positive, and the selectivity towards CO2 decreases with the decrease of Pco. On the other hand, the decrease of the selectivity with the reaction temperature also reflects the higher apparent activation energy for H-2 oxidation than that for CO oxidation.