渤海海域河流相油藏精细描述技术及应用研究

Aiming at the character of Bohaii Sea area and the heterogeneity of fluvial facies reservoir, litho-geophysics experiments and integrated research of geophysical technologies are carried out. To deal with practical problems in oil fields of Bohai area, such as QHD32-6, Southern BZ25-1 and NP35-2 et al., technology of reservoir description based on seismic data and reservoir geophysical methods is built. In this dissertation, three points are emphasized: ①the integration of multidiscipline; ②the application of new methods and technologies; ③the integration of quiescent and dynamic data. At last, research of geology modeling and reservoir numerical simulation based on geophysical data are integrated. There are several innovative results and conclusion in this dissertation: (1)To deal with problems in shallow sea area where seismic data is the key data, a set of technologies for fine reservoir description based on seismic data in Bohai Sea area are built. All these technologies, including technologies of stratigraphic classification, sedimentary facies identification, structure fine characterization, reservoir description, fluid recognition and integration of geological modeling& reservoir numerical simulation, play an important role in the hydrocarbon exploration and development. In the research of lithology and hydrocarbon-bearing condition, petrophysical experiment is carried out. Outdoors inspection and experiment test data are integrated in seismic forward modeling& inversion research. Through the research, the seismic reflection rules of fluid in porosity are generated. Based on all the above research, seismic data is used to classify rock association, identify sedimentary facies belts and recognition hydrocarbon-bearing condition of reservoir. In this research, the geological meaning of geophysical information is more clear and the ambiguity of geophysical information is efficiently reduced, so the reliability in hydrocarbon forecasting is improved. The methods of multi-scales are developed in microfacies research aiming at the condition of shallow sea area in Bohai Sea: ① make the transformation from seismic information to sedimentary facies reality by discriminant analysis; ②in research of planar sedimentary facies, make microfacies research on seismic scale by technologies integration of seismic multi-attributes analysis& optimization, strata slicing and seismic waveform classification; ③descript the sedimentary facies distribution on scales below seismic resolution with the method of stochastic modeling. In the research of geological modeling and reservoir numerical simulation, the way of bilateral iteration between modeling and numerical simulation is carried out in the geological model correction. This process include several steps: ①make seismic forward modeling based on the reservoir numerical simulation results and geological models; ②get trend residual of forward modeling and real seismic data; ③make dynamic correction of the model according to the above trend residual. The modern integration technology of reservoir fine description research in Bohai Sea area, which is developed in this dissertation, is successfully used in (1)the reserve volume evaluation and development research in BZ25-1 oil field and (2)the tracing while drilling research in QHD32-6 oil field. These application researches show wide application potential in hydrocarbon exploration and development research in other oil fields.

叠前叠后联合反演技术及其在QHD32-6油田的应用研究

What geophysical inversion studied includes the common mathematics physical property of inversion and the constitution and appraisal method of solution in geophysics domain, i.e. using observed physical phenomenon from the earth surface to infer space changing and physical property structure of medium within the earth. Seismic inversion is a branch of geophysical inversion. The basic purpose of seismic inversion is to utilizing seismic wave propagating law in the medium underground to infer stratum structure and space distribution of physical property according to data acquisition, processing and interpretation, and then offer the vital foundation for exploratory development. Poststack inversion is convenient and swift, its acoustic impedance inversion product can reflect reservoir interior changing rule to a certain degree, but poststack data lack abundant amplitude and travel time information included in prestack data because of multiple superimpose and weaken the sensitiveness which reflecting reservoir property. Compared with poststack seismic inversion, prestack seismic inversion has better fidelity and more adequate information. Prestack seismic inversion, including waveform inversion, not only suitable for thin strata physical property inversion, it can also inverse reservoir oil-bearing ability. Prestack seismic inversion and prestack elastic impedance inversion maintain avo information, sufficiently applying seismic gathering data with different incident angle, partial angle stack, gradient and intercept seismic data cube. Prestack inversion and poststack inversion technology were studied in this dissertation. A joint inversion method which synthesize prestack elastic wave waveform inversion, prestack elastic impedance inversion and poststack inversion was proposed by making fully use of prestack inversion multiple information and relatively fast and steady characteristic of poststack inversion. Using the proposed method to extract rock physics attribute cube with clear physical significance and reflecting reservoir characterization, such as P-wave and S-wave impedance, P-wave and S-wave velocity, velocity ratio, density, Poisson ratio and Lame’s constant. Regarding loose sand reservoir in lower member of Minghuazhen formation, 32-6 south districts in Qinhuangdao,as the research object, be aimed at the different between shallow layer loose sand and deep layer tight sand, first of all, acquire physical property parameters suitable for this kind of heavy oil pool according to experimental study, establishing initial pressure and shear wave relational model; Afterwards, performing prestack elastic wave forward and inversion research, summarizing rules under the guidance of theoretical research and numerical simulation, performing elastic impedance inversion, calculating rock physics attributes; Finally, predicting sand body distribution according to rock physics parameters, and predicting favorable oil area combine well-logging materials and made good results.

松辽盆地长岭断陷火山岩储层识别与勘探目标评价

Changling fault depression is the biggest fault subsidence in south of Songliao Basin. In its Lower Cretaceous Yingcheng and Shahezi formations developed thick source rocks of deep lake facies and developed poly-phase volcanic rock reservoirs as well. In recent years, significant breakthroughs have been obtained in hydrocarbon exploration of volcanic rock reservoir in the different fault depressions in Songliao basin. Lately, I have been involved in hydrocarbon exploration in the Changling rift depression, especially volcanic rock reservoirs and exploration targets research, participating in the deployment of well Yaoshen 1 which gained over 40 × 104m3 natural gas flow. As quick changes of lithology and facies in Changling area in the south of Songliao basin, and the volcanic rock interludes distribution in continental clastic rock and shale in 3D space, so the identification of volcanic rock types and distribution become a difficult problem. Thus, based on the integrated research of the wild outcrop observation, gravity, magnetic and seismic data, geophysical logging, drilling and coring, laboratory test, this paper carried out the reservoir identification, description and prediction of volcanic rocks in Changling fault depression. In this area, this paper analyzed the volcanic rocks litho-facies, the eruption period, and characteristics of cycles. At the same time, tried to know how to use logging, seismic data to separate volcanic rocks from sandstone and shale, distinguish between volcanic reservoir and non-reservoir, distinguish between intermediate-basic and acidic volcanic rocks, and how to identify traps of volcanic rocks and its gas-bearing properties, etc. Also it is summarized forming conditions and distribution of traps, and possible gas-bearing traps were optimized queuing management. Conclusions as follows: There are two faulted basements in Changling fault depression, granite basement in the southeast and upper paleozoic epimetamorphic basement in the northwest. The main volcanic reservoirs developed in Yingcheng period, which was the intermediate-basic and acidic volcanic eruptions, from the south to north by the intermediate-basic to acid conversion. The volcanic vents are gradually young from south to north. According to information of the re-processing 3D seismic data and gravity-magnetic data, the large volcanic vent or conduit was mainly beaded-distributed along the main fault. The volcanic rocks thickness in Yingcheng formation was changed by the deep faults and basement boundary line. Compared with the clastic rocks, volcanic rocks in Changling area are with high resistance and velocity (4900-5800), abnormal Gamma. All kinds of volcanic rocks are with abnormal strong amplitude reflection on the seismic stacked section except tuff. By analyzing the seismic facies characteristics of volcanic rocks, optimizing seismic attributes constrained by logging, using seismic amplitude and waveforms and other attributes divided volcanic rocks of Yingcheng formation into four seismic zones in map. Currently, most volcanic gas reservoirs are fault-anticline and fault-nose structure. But the volcanic dome lithologic gas reservoirs with large quantity and size are the main gas reservoir types to be found.

多波多分量地震资料走时近似方法的研究与应用

Compared with the conventional P wave, multi-component seismic data can markedly provide more information, thus improve the quality of reservoir evaluation like formation evaluation etc. With PS wave, better imaging result can be obtained especially in areas involved with gas chimney and high velocity formation. However, the signal-to-noise of multi-component seismic data is normally lower than that of the conventional P wave seismic data, while the frequency range of converted wave is always close to that of the surface wave which adds to the difficulty of removing surface wave. To realize common reflection point data stacking from extracted common conversion point data is a hard nut to crack. The s wave static correction of common receiver point PS wave data is not easy neither. In a word, the processing of multi-component seismic data is more complicated than P wave data. This paper shows some work that has been done, addressing those problems mentioned above. (1) Based on the AVO feature of converted wave, this paper has realized the velocity spectrum of converted waves by using Sarkar’s generalized semblance method taking into account of AVO factor in velocity analysis. (2)We achieve a method of smoothly offset division normal method.Firstly we scan the stacking velocities in different offset divisions for a t0, secondly obtain some hyperbolas using these stacking velocities, then get the travel time for every trace using these hyperbolas; in the end we interpolate the normal move out between two t0 for every trace. (3) Here realize a method of stepwise offset division normal moveout.It is similar to the method of smoothly offset division normal moveout.The main difference is using quadratic curve, sixth order curve or fraction curve to fit these hyperbolas. (4)In this paper, 4 types of travel time versus distance functions in inhomogeneous media whose velocity or slowness varies with depth and vertical travel time have been discussed and used to approximate reflection travel time. The errors of ray path and travel time based on those functions in four layered models were analyzed, and it has shown that effective results of NMO in synthetic or real data can be obtained. (5) The velocity model of converted PS-wave can be considered as that of P -wave based on the ghost source theory, thus the converted wave travel time can be approximated by calculation from 4 equivalent velocity functions: velocity or slowness vary linearly with depth or vertical travel time. Then combining with P wave velocity analysis, the converted wave data can be corrected directly to the P-wave vertical travel time. The improvements were shown in Normal Move out of converted waves with numerical examples and real data. (6) This paper introduces the methods to compute conversion point location in vertical inhomogeneous media based on linear functions of velocity or slowness versus depth or vertical travel time, and introduce three ways to choose appropriate equivalent velocity methods, which are velocity fitting, travel time approximation and semblance coefficient methods.

基于窗口精细结构描述的岩体质量评价及力学参数综合优选

This thesis bases on horizontal research project “The research about the fine structure and mechanical parameters of abutment jointed rock mass of high arch dam on Jinping Ⅰ Hydropower Station, Yalong River” and “The research about the fine structure and mechanical parameters of the columnar basalt rock mass on Baihetan Hydropower Station, Jinsha River”. A rounded system about the fine structure description and rock mass classification is established. This research mainly contains six aspects as follow: (1) Methods about fine structure description of the window rock mass; (2) The window rock mass classification about the fine structure; (3) Model test study of intermittent joints; (4) Window rock mass strength theory; (5) Numerical experimentations about window rock mass; (6) The multi-source fusion of mechanical parameters based on Bayes principle. Variation of intact rock strength and joint conditions with the weathering and relaxation degree is studied through the description of window rock mass. And four principal parameters: intact rock point load strength, integration degree of window rock mass, joint conditions, and groundwater condition is selected to assess the window rock mass. Window rock mass is classified into three types using the results of window rock mass fine structure description combined with joints develop model. Scores about intact rock strength, integrality condition, divisional plane condition and groundwater conditions are given based on window rock mass fine structure description. Then quality evaluation about two different types of rock mass: general joint structure and columnar jointing structure are carried out to use this window rock mass classification system. Application results show that the window rock mass classification system is effective and applicable. Aimed at structural features of window structure of “the rock mass damaged by recessive fracture”, model tests and numerical models are designed about intermittent joints. By conducting model tests we get shear strength under different normal stress in integrated samples, through samples and intermittent joints samples. Also, the changing trends of shear strength in various connectivity rates are analyzed. We numerically simulate the entire process of direct shear tests by using PFC2D. In order to tally the stress-strain curve of numerical simulation with experimental tests about both integrated samples and through samples, we adjust mechanical factors between particles. Through adopting the same particle geometric parameter, the numerical sample of intermittent joints in different connective condition is re-built. At the same time, we endow the rock bridges and joints in testing samples with the fixed particle contacting parameters, and conduct a series of direct shear tests. Then the destructive process and mechanical parameters in both micro-prospective and macro-prospective are obtained. By synthesizing the results of numerical and sample tests and analyzing the evolutionary changes of stress and strain on intermittent joints plane, we conclude that the centralization of compressive stress on rock bridges increase the shear strength of it. We discuss the destructive mechanics of intermittent joints rock under direct shear condition, meanwhile, divide the whole shear process into five phases, which are elasticity phase, fracture initiation phase, peak value phase, after-peak phase and residual phase. In development of strength theory, the shear strength mechanisms of joint and rock bridge are analyzed respectively. In order to apply the deducted formulation conveniently in the real projects, a relationship between these formulations and Mohr-Coulomb hypothesis is built up. Some sets of numerical simulation methods, i.e. the distinct element method (UDEC) based on in-situ geology mapping are developed and introduced. The working methods about determining mechanical parameters of intact rock and joints in numerical model are studied. The operation process and analysis results are demonstrated detailed from the research on parameters of rock mass based on numerical test in the Jinping Ⅰ Hydropower Station and Baihetan Hydropower Station. By comparison，the advantages and disadvantages are discussed. Results about numerical simulation study show that we can get the shear strength mechanical parameters by changing the load conditions. The multi-source rock mass mechanical parameters can be fused by the Bayes theory, which are test value, empirical value and theoretical value. Then the value range and its confidence probability of different rock mass grade are induced and these data supports the reliability design.

波前重建法射线追踪及走时表压缩

As a fast and effective method for approximate calculation of seismic numerical simulation, ray tracing method, which has important theory and practical application value, in terms of seismic theory and seismic simulation, inversion, migration, imaging, simplified from seismic theory according to geometric seismic, means that the main energy of seismic wave field propagates along ray paths in condition of high-frequency asymptotic approximation. Calculation of ray paths and traveltimes is one of key steps in seismic simulation, inversion, migration, and imaging. Integrated triangular grids layout on wavefront with wavefront reconstruction ray tracing method, the thesis puts forward wavefront reconstruction ray tracing method based on triangular grids layout on wavefront, achieves accurate and fast calculation of ray paths and traveltimes. This method has stable and reasonable ray distribution, and overcomes problems caused by shadows in conventional ray tracing methods. The application of triangular grids layout on wavefront, keeps all the triangular grids stable, and makes the division of grids and interpolation of a new ray convenient. This technology reduces grids and memory, and then improves calculation efficiency. It enhances calculation accuracy by accurate and effective description and division on wavefront. Ray tracing traveltime table, which shares the character of 2-D or 3-D scatter data, has great amount of data points in process of seismic simulation, inversion, migration, and imaging. Therefore the traveltime table file will be frequently read, and the calculation efficiency is very low. Due to these reasons, reasonable traveltime table compression will be very necessary. This thesis proposes surface fitting and scattered data compression with B-spline function method, applies to 2-D and 3-D traveltime table compression. In order to compress 2-D (3-D) traveltime table, first we need construct a smallest rectangular (cuboidal) region with regular grids to cover all the traveltime data points, through the coordinate range of them in 2-D surface (3-D space). Then the value of finite regular grids, which are stored in memory, can be calculated using least square method. The traveltime table can be decompressed when necessary, according to liner interpolation method of 2-D (3-D) B-spline function. In the above calculation, the coefficient matrix is stored using sparse method and the liner system equations are solved using LU decomposition based on the multi-frontal method according to the sparse character of the least square method matrix. This method is practiced successfully in several models, and the cubic B-spline function can be the best basal function for surface fitting. It make the construction surface smooth, has stable and effective compression with high approximate accuracy using regular grids. In this way, through constructing reasonable regular grids to insure the calculation efficiency and accuracy of compression and surface fitting, we achieved the aim of traveltime table compression. This greatly improves calculation efficiency in process of seismic simulation, inversion, migration, and imaging.

地震波场模拟及叠前深度偏移的李群方法

The processes of seismic wave propagation in phase space and one way wave extrapolation in frequency-space domain, if without dissipation, are essentially transformation under the action of one parameter Lie groups. Consequently, the numerical calculation methods of the propagation ought to be Lie group transformation too, which is known as Lie group method. After a fruitful study on the fast methods in matrix inversion, some of the Lie group methods in seismic numerical modeling and depth migration are presented here. Firstly the Lie group description and method of seismic wave propagation in phase space is proposed, which is, in other words, symplectic group description and method for seismic wave propagation, since symplectic group is a Lie subgroup and symplectic method is a special Lie group method. Under the frame of Hamiltonian, the propagation of seismic wave is a symplectic group transformation with one parameter and consequently, the numerical calculation methods of the propagation ought to be symplectic method. After discrete the wave field in time and phase space, many explicit, implicit and leap-frog symplectic schemes are deduced for numerical modeling. Compared to symplectic schemes, Finite difference (FD) method is an approximate of symplectic method. Consequently, explicit, implicit and leap-frog symplectic schemes and FD method are applied in the same conditions to get a wave field in constant velocity model, a synthetic model and Marmousi model. The result illustrates the potential power of the symplectic methods. As an application, symplectic method is employed to give synthetic seismic record of Qinghai foothills model. Another application is the development of Ray+symplectic reverse-time migration method. To make a reasonable balance between the computational efficiency and accuracy, we combine the multi-valued wave field & Green function algorithm with symplectic reverse time migration and thus develop a new ray+wave equation prestack depth migration method. Marmousi model data and Qinghai foothills model data are processed here. The result shows that our method is a better alternative to ray migration for complex structure imaging. Similarly, the extrapolation of one way wave in frequency-space domain is a Lie group transformation with one parameter Z and consequently, the numerical calculation methods of the extrapolation ought to be Lie group methods. After discrete the wave field in depth and space, the Lie group transformation has the form of matrix exponential and each approximation of it gives a Lie group algorithm. Though Pade symmetrical series approximation of matrix exponential gives a extrapolation method which is traditionally regarded as implicit FD migration, it benefits the theoretic and applying study of seismic imaging for it represent the depth extrapolation and migration method in a entirely different way. While, the technique of coordinates of second kind for the approximation of the matrix exponential begins a new way to develop migration operator. The inversion of matrix plays a vital role in the numerical migration method given by Pade symmetrical series approximation. The matrix has a Toepelitz structure with a helical boundary condition and is easy to inverse with LU decomposition. A efficient LU decomposition method is spectral factorization. That is, after the minimum phase correlative function of each array of matrix had be given by a spectral factorization method, all of the functions are arranged in a position according to its former location to get a lower triangular matrix. The major merit of LU decomposition with spectral factorization (SF Decomposition) is its efficiency in dealing with a large number of matrixes. After the setup of a table of the spectral factorization results of each array of matrix, the SF decomposition can give the lower triangular matrix by reading the table. However, the relationship among arrays is ignored in this method, which brings errors in decomposition method. Especially for numerical calculation in complex model, the errors is fatal. Direct elimination method can give the exact LU decomposition But even it is simplified in our case, the large number of decomposition cost unendurable computer time. A hybrid method is proposed here, which combines spectral factorization with direct elimination. Its decomposition errors is 10 times little than that of spectral factorization, and its decomposition speed is quite faster than that of direct elimination, especially in dealing with a large number of matrix. With the hybrid method, the 3D implicit migration can be expected to apply on real seismic data. Finally, the impulse response of 3D implicit migration operator is presented.

孤北洼陷沙三-沙四段隐蔽油藏形成机制及其预测

Based on the study of sequence stratigraphy, modern sedimentary, basin analysis, and petroleum system in Gubei depression, this paper builds high resolution sequence stratigraphic structure, sedimentary system, sandbody distribution, the effect of tectonic in sequence and sedimentary system evolution and model of tectonic-lithofacies. The pool formation mechanism of subtle trap is developed. There are some conclusions and views as follows. 1.With the synthetic sequence analysis of drilling, seismic, and well log, the highly resolution sequence structure is build in Gubei depression. They are divided two secondary sequences and seven three-order sequences in Shahejie formation. They are include 4 kinds of system traces and 7 kinds of sedimentary systems which are alluvial fan, under water fan, alluvial fan and fan-delta, fan-delta, lacustrine-fan, fluvial-delta-turbidite, lakeshore beach and bar, and deep lake system. Sandbody distribution is show base on third order sequence. 2.Based on a lot of experiment and well log, it is point out that there are many types of pore in reservoir with the styles of corrosion pore, weak cementing, matrix cementing, impure filling, and 7 kinds of diagenetic facies. These reservoirs are evaluated by lateral and profile characteristics of diagenetic facies and reservoir properties. 3.The effect of simultaneous faulting on sediment process is analyzed from abrupt slope, gentle slope, and hollow zone. The 4 kinds of tectonic lithofacies models are developed in several periods in Gubei depression; the regional distribution of subtle trap is predicted by hydro accumulation characteristics of different tectonic lithofacies. 4.There are 4 types of compacting process, which are normal compaction, abnormal high pressure, abnormal low pressure and complex abnormal pressure. The domain type is normal compaction that locates any area of depression, but normal high pressure is located only deep hollow zone (depth more than 3000m), abnormal low pressures are located gentle slope and faulted abrupt slope (depth between 1200~2500m). 5.Two types dynamic systems of pool formation (enclosed and partly enclosed system) are recognized. They are composed by which source rocks are from Es3 and Es4, cap rocks are deep lacustrine shale of Esl and Es3, and sandstone reservoirs are 7 kinds of sedimentary system in Es3 and Es4. According to theory of petroleum system, two petroleum systems are divided in Es3 and Es4 of Gubei depression, which are high or normal pressure self-source system and normal or low pressure external-source system. 6.There are 3 kinds of combination model of pool formation, the first is litholgical pool of inner depression (high or normal pressure self-source type), the second is fault block or fault nose pool in marginal of depression (normal type), the third is fault block-lithological pool of central low lifted block (high or normal pressure type). The lithological pool is located central of depression, other pool are located gentle or abrupt slope that are controlled by lithological, faulting, unconfirmed. 7.This paper raise a new technique and process of exploration subtle trap which include geological modeling, coring description and logging recognition, and well log constrained inversion. These are composed to method and theory of predicting subtle trap. Application these methods and techniques, 6 hydro objects are predicted in three zone of depression.

东营凹陷中央隆起带次生油藏成藏机制和预测

It is the key project of SINOPEC at ninth five years period with a lot of work and very difficult, which the main object are the study of pool-forming mechanism, distribution rule and pool-forming model of complex secondary pool at Dongying formation in high mature exploration area, and building theories and methods of research, description and prediction of secondary fault block pool. This paper apply comprehensively with various theories, method and techniques of geology, seismic, well log, reservoir engineering, meanwhile apply with computer means, then adopt combination of quality and quantitative to develop studies of pool-forming mechanism, model and pool prediction of fault block pool. On the based of stretch, strike-slip, reversal structure theories, integrated the geometry, kinematics, and dynamics of structure, it is show that the structure framework, the structure evolve, formation mechanism of central uplift belt of Dongying depression and control to formation and distribute of secondary complex fault block pool. The opening and sealing properties, sealing mechanism and sealing models of pool-controlling fault are shown by using quality, direction of normal stress, relations between interface and rock of two sides of fault and shale smear factor (SSF), as well as the juxtaposition of fault motion stage and hydrocarbon migration, etc. The sealing history of controlling fault, formation mechanism and distribute the regulation are established by combining together with bury history, structure evolve history, fault growth history stress field evolve history, which can be guide exploration and production oil field. It were bring up for the first time the dynamics mechanism of Dongying central uplift which were the result of compound tress field of stretch, strike-slip and reversal, companion with reversal drag structure, arcogenesis of paste and salt beds. The dual function of migration and sealing of fault were demonstrated in the research area. The ability of migration and sealing oil of pool-controlling fault is controlled by those factors of style of fault combination, activity regulation and intensity of fault at the period of oil migration. The four kinds of sealing model of pool-controlling fault were established in the research area, which the sealing mechanism of fault and distribution regulation of oil in time and space. The sealing ability of fault were controlled by quality, direction of normal stress, relations between interface and rock of two sides of fault and shale smear factor (SSF), as well as the juxtaposition of fault motion stage and hydrocarbon migration, etc. The fuzzy judge of fault sealing is the base of prediction of secondary pool. The pool-forming model of secondary was established in the research area, which the main factors are ability migration and sealing. The transform zone of fault, inner of arc fault and the compound area of multi fault are enrichment region of secondary pool of Dongying formation, which are confirm by exploration with economic performance and social performance.

砂砾岩储层流动单元四维模型研究

This article is an important part of "95" technological subject of SINOPEC. It has a large number of difficulties and workloads, and has significant theoretical meanings and practical value. The study area is composed of sandstone & conglomerate reservoir of alluvial fan & fan delta, which belong to Sha3 lower member and Sha4 upper member of lower tertiary of Yong'an Town Oilfield in Dongying Depression. The target stataum develops in the hanging wall of the synsedimentary fault in the scarp zone of Dongying Depression. The frequently intense movements result in the variation of sandstone and conglomerate reservoir and the evolution of the time and space of Sha3 lower member and Sha4 upper member in Yong'an Town Oilfield. As a result, it is difficult for the individual reservoir correlation at the root of fan, which bring about a tackle problem for the exploitation of oilfield. In this background, the research of fluid units will be more difficult. In this article, the new concepts, the new methods, and the new techniques of sedimentology, petroleum geology, reservoir geology, physics of crystal surface, dynamic & static state reservoir description and well logging geology are synthetically applied, and the computer technology are made full uses of, and the identifying, dividing and appraising of the two-formation-type sandstone & conglomerate reservoir fluid units of Sha3 lower member and Sha4 upper member systemically analyzed in Yong'an Town Oilfield, Dongying Depression. For the first time, the single-well model, the section model, the plane model, the nuclear magnetism log model, the microcosmic network model, the 4-D geology model and the simulation model of the two-formation-type reservoir fluid units of the of sandstone & conglomerate reservoir of Sha3 lower member and Sha4 upper member are established, and the formative mechanism and distributing & enrichment laws of oil-gas of the two type of sandstone and conglomerate reservoir fluid units are revealed. This article established the optimizing, identifying, classifying and appraising standard of the two-formation-type reservoir fluid units of the of sandstone and conglomerate reservoir of Sha3 lower member and Sha4 upper member, which settles the substantial foundations for static state model of the fluid units, reveals the macroscopic & microcosmic various laws of geometrical static state of the fluid units, and instructs the oil exploitation. This article established static state model of the two-formation-type sandstone and conglomerate reservoir fluid units by using the multi-subject theories, information and techniques, and reveals the geometrical configuration, special distribution and the oil-gas enrichment laws of the sandstone and conglomerate reservoir fluid units. For the first time, we established the nuclear magnetism log model of the two-formation-type sandstone and conglomerate reservoir of Sha3 lower member and Sha4 upper member, which reveals not only the character and distributing laws of the porosity and permeability, bat also the formation and distribution of the movable fluid. It established six type of microcosmic net model of the two-formation-type sandstone and conglomerate reservoir of Sha3 lower member and Sha4 upper member in the working area by using the advanced theories, such as rock thin section, SEM, image analysis, intrusive mercury, mold, rock C.T. measure & test image etc., which reveals the microcosmic characteristic of porosity & throat, filterate mode and microcosmic oil-gas enrichment laws of the sandstone and conglomerate reservoir. For the first time, it sets up the 4-D model and mathematic model of the sandstone and conglomerate reservoir, which reveals the distributing and evolving laws of macroscopic & microcosmic parameters of the two-formation-type sandstone and conglomerate reservoir and oil-gas in 4-D space. At the same time, it also forecasts the oil-gas distribution and instructs the oilfield exploitation. It established reservoir simulation model, which reveals the filterate character and distributing laws of oil-gas in different porosity & throat net models. This article established the assistant theories and techniques for researching, describing, indicating and forecasting the sandstone and conglomerate reservoir fluid units, and develops the theories and techniques of the land faces faulted basin exploitation geology. In instructing oilfield exploitation, it had won the notable economic & social benefits.

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

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.

港西油田现代油藏管理研究与实践

Through years of practice, reservoir management has already become the basic mode of foreign oil companies to realize the high-efficient development of the oil field. From the view of reservoir development and technological economy, reservoir management regards the study of the reservoir engineering, designs of reservoir projects and the dynamic analysis of the reservoir's performance as a system. In the fields of reservoir description, the establishment of the geological models and development models, the dynamic simulations of reservoir exploitation and the design of the oil engineering, reservoir management emphasizes the cooperation of the geology and the engineering, the combination of the engineering technology and the economic evaluation. In order to provide the means and basis for the reservoir geology study, reservoir evaluation, reserves calculation, numerical simulation, development plan and risk analysis, it adopts the reservoir management activities(team work) to make and implement the optimized oil field development management strategies so that secientific and democratic decision making can be achieved. Under the planned economic system for a long time, the purpose of Chinese reservoir development has been to fulfill the" mandatory" production task. With the deepening of the reform, the management organization of Chinese petroleum enterprises has been gradually going through the transition and reforms to the operational entity and the establishment of the mode of oil companies under the socialist market economy system. This research aims at introducing the advanced reservoir management technique from foreign countries to further improve the reservoir development results and wholly raise the economic benefits of Chinese mature land facieses sandstone reservoirs in the later stage of the water flooding. We are going to set up a set of modern reservoir management modes according to the reservoir features, current situation and existing problems of GangXi oil field of DaGang oil company. Through the study and implementation of the reservoir description and numerical simulation technology effectively, we plan to work out integrated adjustment projects, to study the related technology of oil recovery; to set up the effective confirmable data procedure and data management system of the reservoir management, to establish the coordinated model and workbench related to geology, engineering and economy in order to realize the real time supervision and evaluation on the process of reservoir development. We hope to stipulate modernization management tools for GangXi oil fields to rationally utilize various kinds of existing technological methods and to realize the economic exploitation and achieve the maximum benefits from the reservoir. The project of the modem reservoir management will be carried out on the GangXi oil field of DaGang oil company for this oil field is typical and has integrated foundamental materials and perfect networks. Besides, it is located in the good geographical position enjoying very convenient traffic. Implementing modern reservoir management will raise the recovery ratio, reduce the production cost and improve the working efficiency. Moreover, the popularization of modern reservoir management will improve the comprehensive benefits of DaGang oil company and even the whole Petro China. Through the reserch of this project, the following technical indicators can be reached: Establishing the concept of modern reservoir management. Establishing a set of integrated data information management system adapt to the features of GangXi reservoir. 3. Forming technical research modes of modern reservoir management suitable for mature reservoirs in the later developing stage. 4. Advancing projects of GangXi reservoir which are maxium optimized in engineering technique and economic benefits of oil exploitation. Besides, this set of technology, research principle and method can guide the mature reservoir of DaGang oil field and even the whole PetroChina to develop the further research of reservoir adjustment and improve the reservoir recovery factor and developing level constantly.

宁夏固原寺口子剖面磁性地层与孢粉植物群初步研究

In Asia continent, several significant environmental events happened during the Cenozoic era, such as uplift of the Tibet Plateau, formation of the Asian Monsoon system, aridification of the temperate inland region in Central Asia. To investigate the history of long-term palaeoclimate evolution during the late Cenozoic, a lacustrine sequence located at Sikouzi, Guyuan county, Ningxia-Hui Autonomous Region, was studied. The Sikouzi section is about 2880 m in thickness and has a general continuous nature according to field observation. Thus this thick lacustrine record is an important archive to further understand those environmental events. In this study, detailed field measurement, layer-after-layer description and sampling, and magnetostratigraphy and palynoflora investigations are conducted at the Sikouzi section and some preliminary results have been achieved as follows. Based on Hipparion fauna, pollen data and long distinctive patterns of the local magnetozones, the confident correlation of the Sikouzi magnetostratigraphic polarity to the GPTS (CK95) is best established, indicating that the top boundary of the Sikouzi formation is dated back to -19.8 Ma B.P. and consequently no Oligocene sediments deposited locally. On the other hand, both the field observation and the correlation to GPTS indicate a nearly continuous nature for the whole sequence. The palynological results show that the grassland has been a dominant vegetation in the Sikouzi area since ~19.8 Ma B.P, although some trees/shrubs were present sparsely during the intervals of relatively warm and wet climatic conditions. This implies that the onset of the aridification in northwestern China is dated back to at least 19.8 Ma B.P. ago. The Neogene Global Climatic Optimum (-16.0 Ma B.P.) occurring between the late Early Miocene and the early Middle Miocene is well documented in the Shanwang formation, Shandong Province. However, such event was not found in the Sikouzi record, and neither in the Guide and the Hualong basins, Qinghai province. This may lead us to the conclusion that the East-Asian Summer Monsoon system remained weak during the period of the Miocene Climatic Optimum although the onset of it was traced back to the Early Oligocene. In the Sikouzi area, it was warm between the latest Miocene and the early Pliocene and then became cool in the late Pliocene. This pattern is consistent with the palaeoclimate record of the Pliocene from other areas in the world.

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

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.

新疆自然区划单元土地利用/覆盖分布图谱研究

Under the direction of Geo-informatic Tupu theory, based on comprehensive natural division in Xinjiang and 1:100000 land use and land cover vector data in 2000 from resource and environment database of CAS, the paper gave out some distribution Tupu of land resource types in different regions of Xinjiang. GIS tools such as ARCTOOLS and ARCV1EW were used to clip the unit of each natural division. Some useful conclusions were established. Then the lower reaches of Tarim River were selected as typical area to analyze the ecological environment evolvement from 2000 to 2004 since the beginning of ecological water delivering. The comprehensive space-time analysis provided a method to monitor the effect dynamically. The main contents of the thesis are listed as follows: (1) Similarities and differences between North Xinjiang and South Xinjiang, 6 second-classified regions and 37 third-classified regions. The conclusion of each region emphasized the percentage of area of main land types, characteristic description and distribution of cultivated land > woodlands grassland-, water area^ urban and rural land and unused land. (2) Thematic Tupu of each region. It concluded artificial oasis Tupu, the proportion between dense-grass% moderate-grass and sparse-grass, pattern of land resources and the land use degree. The artificial oasis mean cultivated land^ woodland^ urban & rural land and other construction land. The proportion of grassland disclosed the quality of all grassland and showed the development of them. Pattern of land resources and the land use degree gave out the detailed development direction and development degree of each region. (3) Ecological environment evolvement of the lower reaches of Tarim River. The effects of each ecological water delivering were compared.