地质导向钻井关键技术及方法

On the basis of analyzing the principle and realization of geo-steering drilling system, the key technologies and methods in it are systematically studied in this paper. In order to recognize lithology, distinguish stratum and track reservoirs, the techniques of MWD and data process about natural gamma, resistivity, inductive density and porosity are researched. The methods for pre-processing and standardizing MWD data and for converting geological data in directional and horizontal drilling are discussed, consequently the methods of data conversion between MD and TVD and those of formation description and adjacent well contrast are proposed. Researching the method of identifying sub-layer yields the techniques of single well explanation, multi-well evaluation and oil reservoir description. Using the extremum and variance clustering analysis realizes logging phase analysis and stratum subdivision and explanation, which provides a theoretical method and lays a technical basis for tracing oil reservoirs and achieving geo-steering drilling. Researching the technique for exploring the reservoir top with a holdup section provides a planning method of wellpath control scheme to trace oil and gas reservoir dynamically, which solves the problem of how to control well trajectory on condition that the layer’s TVD is uncertain. The control scheme and planning method of well path for meeting the demands of target hitting, soft landing and continuous steering respectively provide the technological guarantee to land safely and drill successfully for horizontal, extended-reach and multi-target wells. The integrative design and control technologies are researched based on geology, reservoir and drilling considering reservoir disclosing ratio as a primary index, and the methods for planning and control optimum wellpath under multi-target restriction, thus which lets the target wellpath lie the favorite position in oil reservoir during the process of geo-steering drilling. The BHA (bottomhole assembly) mechanical model is discussed using the finite element method, and the BHA design methods are given on the basis of mechanical analyses according to the shape of well trajectory and the characteristics of BHA’s structure and deformation. The methods for predicting the deflection rate of bent housing motors and designing their assemblies are proposed based on the principle of minimum potential energy, which can clearly show the relation between the BHA’s structure parameters and deflection rate, especially the key factors’ effect to the deflection rate. Moreover, the interaction model between bit and formation is discussed through the process of equivalent formation and equivalent bit considering the formation anisotropy and bit anisotropy on the basis of analyzing the influence factors of well trajectory. Accordingly, the inherence relationship among well trajectory, formation, bit and drilling direction is revealed, which lays the theory basis and technique for predicting and controlling well trajectory.

VTI各向异性叠前深度偏移及速度分析

Theoretical research, laboratory test and field observation show that most of sediment rock has anisotropic features. It will produce some notable errors when applying isotropic methods such as prestack depth migration and velocity analysis to dada acquired under anisotropic condition; it also has a bad effect on geologic interpretation. Generally speaking, the vertical transverse isotropic media is a good approximation to geologic structure, thus it has an important realistic meaning for anisotropic prestack depth migration theory researching and precise complex geologic imaging if considering anisotropic effect of seismic wave propagation. There are two indispensable parts in prestack depth migration of realistic records, one is proper prestack depth migration algorithm, and the other is velocity analysis using prestack seismic data. The paper consists of the two aspects. Based on implicit finite difference research proposed by Dietrich Ristow et al （1997） about VTI media prestack depth migration, the paper proposed split-step Fourier prestack depth migration algorithm （VTISSF） and Fourier finite difference algorithm （VTIFFD） based on wave equation for VTI media, program are designed and the depth migration method are tested using synthetic model. The result shows that VTISSF is a stable algorithm, it generally gets a good result if the reflector dip is not very steep, while undermigration phenomena appeared in steep dips case; the VTIFFD algorithm bring us better result in steep dips with lower efficiency and frequency dispersion. For anisotropic prestack depth migration velocity analysis of VTI media, The paper discussed the basic hypothesis of VTI model in velocity analysis algorithm, basis of anisotropic prestack depth migration velocity analysis and travel time table calculation of VTI media in integral prestack depth migration. Then , analyzed the P-wave common imaging gather in the case of homogeneous velocity and vertically variable velocity . studied the residual correction in common imaging gather produced by media parameter error, analyzed the condition of flat event and correct depth in common imaging gather . In this case, the anisotropic model parameter vector is , is vertical velocity of a point at top surface, is vertical velocity gradient, and are anisotropic parameter. We can get vertical velocity gradient from seismic data; then the P-wave common imaging gather of VTI media whose velocity varies in vertical and horizontal direction, the relationship between media parameter and event residual time shift of common image gather are studied. We got the condition of flattening common imaging gather with correct depth. In this case the anisotropic model parameter vector is , is velocity gradient in horizontal direction. As a result, the vertical velocity grads can be decided uniquely, but horizontal velocity grads and anisotropic parameter can’t be distinguished if no priori information available, our method is to supply parameter by velocity scanning; then, as soon as is supplied we can get another four parameters of VTI media from seismic data. Based on above analysis, the paper discussed the feasibility of migration velocity analysis in vertically and horizontally varied VTI media, synthetic record of three models are used to test the velocity analysis method . Firstly, anisotropic velocity analysis test is done using a simple model with one block, then we used a model with multiple blocks, thirdly, we analyzed the anisotropic velocity using a part of Marmousi model. The model results show that this velocity analysis method is feasible and correct.

构造复杂区变速成图方法研究

In the complex structure areas, velocity field building and structure mapping are important for seismic exploration. With the development of seismic exploration, the methods of structure mapping, reservoir prediction and reservoir description all require high precious velocity field. And more accurate depth-structure maps are required for well site design. Aiming at the problems and defects in velocity analysis and structure mapping in oil seismic exploration, the paper which is based on the studies of real data in several areas combines the theories with practical application, and analyzes the precision and applicability of several methods of velocity model building. After that, the following methods are mainly studied: the coherence inversion methods based on the pre-stack CMP gathers or stacking velocity; the interval velocity inversion methods constrained by multi-well; the Random Simulation method; 3D Image Ray Map Migration method and the structure mapping in floating datum and in fixed datum, and then we conclude the method of building high precious seismic velocity field and structure mapping with variable velocity. Firstly, the paper analyses the distributing rule of the velocity variation in the areas with complex structures in the northwest of China, then points out that velocity is a crucial factor which influences the precision of structure mapping, and the velocity variations have something to do with the shapes of the structures, the variety of lithology and so on. The key point of improving the precision of seismic velocity field is to obtain a structure mapping with high precision. We also describe the range and conditions of these methods. Secondly, by comparing many popular methods of velocity model building, we propose a new method in the use of velocity model building. The new method is more effective in velocity model building under every kind of complex condition and is worthy of spreading. At last, the paper fingers out that it is a system engineering to study variable velocity mapping in every kind of complex structure areas. Every step of the work can affect the final results. So it is important to build high efficient and practical velocity model and the flows of mapping processing. The paper builds the flows and gives some examples. The method has been applied in more than ten exploring surveys. The application proves that this method could bring good effect on researching on low-amplitude trap, reservoir prediction, reservoir description and the integrated research of oil&gas geology. Keywords: structure mapping velocity model building complex structure variable velocity media

海拉尔盆地贝尔凹陷基岩潜山特征及储层预测

The lithology of the buried hill of Triassic Budate group in Beier depression is epimetamorphic clastic rock and volcanic clastic rock stratum. Recently the favorable hydrocarbon show was discovered in buried hill of base rock, and large-duty industrial oil stream was obtained in some wells in Beier depression. Based on the information of seismos and wells, the tectonic framework, tectonic deformation times and faulted system of the Beier depression are comprehensively studied, then configuration, evolutional history, genetic type and distributed regularity of buried hill are defined. According to observing description and analysis of core sample, well logging and interpretive result of FMI, the lithological component, diagenetic type and diagenetic sequence of buried hill reservoir are confirmed, then reservoir space system of buried hill is distinguished, and vegetal feature, genetic mechanism and distributed regularity of buried hill fissure are researched, at the same time the quantitative relationship is build up between core fissures and fissures interpreted by FMI. After that fundamental supervisory action of fault is defined to the vegetal degree of fissure, and the fissure beneficial places are forecasted using fractal theory and approach. At last the beneficial areas of Budate group reservoir are forecasted by reservoir appraisal parameters optimization such as multivariate gradually regression analysis et. al. and reservoir comprehensive appraisal method such as weighing analyze and clustering procedure et. al. which can provide foundation for the next exploratory disposition. Such production and knowledge are obtained in this text as those: 1. Four structural layers and two faulting systems are developed, and four structural layers are carved up by three bed succession boundary surfaces which creates three tectonic distortional times homology. Three types of buried hill are divided, they are ancient physiognomy buried hill, epigenetic buried hill, and contemporaneous buried hill. 2. Reservoir space of Budate buried hill is mainly secondary pore space and fissure, which distributes near the unconformity and/or inside buried hill in sections. The buried hill reservoir experienced multi-type and multi-stage diagenetic reconstruction, which led to the original porosity disappeared, and multi secondary porosity was created by dissolution, superficial clastation and cataclasis et. al. in diagenetic stage, which including middle crystal pore, inter crystal pore, moldic pore, inter particle emposieu, corrosion pore space and fissure et. al. which improved distinctly reservoir capability of buried hill. 3. The inner reservoir of buried hill in Beier depression is not stratigraphic bedded construction, but is fissure developing place formed by inner fault and broken lithogenetic belt. The fissures in inner reservoir of buried hill are developed unequally with many fissure types, which mainly are high angle fissure and dictyonal fissures and its developing degree and distribution is chiefly controlled by faulting. 4. The results of reservoir comprehensive evaluate and reservoir predicting indicates that advantageous areas of reservoir of buried hill chiefly distributes in Sudeerte, Beixi and Huoduomoer, which comprehensive evaluate mainly Ⅱand Ⅲ type reservoir. The clues and results of this text have directive significance for understanding the hydrocarbon reservoir condition of buried hill in Beier depression, for studying hydrocarbon accumulated mechanism and distributed regularity, and for guiding oil and gas exploration. The results of this text also can enrich and improve nonmarine hydrocarbon accumulated theory.

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

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.