弱各向异性介质中传播的地震波及其速度重建

We begin our studies to make the best of information of seismic data and carry out the description of cracks parameters by extracting anisotropic information. The researching contents are: (1) velocity and polarization anomaly of seismic wave (qP and qSV wave) in weak anisotropic media; (2) reflection seismic synthetic record in anisotropic media; (3) multiple scattering induced by cracks; (4) anisotropic structure inversion and velocity reconstruction with VSP (Vertical Seismic Profile) data; (5) multi-parameters analysis of anisotropy in time-domain and depth-domain. Then we obtain results as follows: (1) We achieve approximate relation of qP and qSV wave's velocity and polarization property in weak anisotropic media. At the same time, we calculate anisotropic velocity factors and polarization anomaly of several typical sedimentary rocks. The results show there are different anisotropic velocity factors and polarization anomaly in different rocks. It is one of the primary theoretical foundation which is expected to identify lithology; (2) We calculate reflection seismic synthetic record with theoretical model; (3) We simulate scattering induced by cracks with Boundary Element Method. Numerical studies show that in the presence of cracks; spatial and scale-length distributions are important and cannot be ignored in modeling cracked solids; (4) From traveltimes information of VSP data, we study the velocity parameter inversion of seismic wave under isotropic and anisotropic models, and its result indicate that the inversion imaging under anisotropic model will not destroy the original features of isotropic model, but it will bring on some bigger error if we adopt the method of isotropic model for anisotropic model data. Further more, basing on the study we develop the CDP mapping technology of reflecting structure under isotropic and anisotropic models, and we process real data as a trial of the methods; (5) We study the problem of initial model reconstruction of anisotropic parameters structure represented by Anderson parameter in depth domain for surface data.

各向异性介质中的地震面波和S波分裂

Surface wave propagation in the anisotropic media and S-wave splitting in China mainland are focused in this M.S. dissertation. We firstly introduced Anderson parameters in the research of surface wave propagation in the anisotropic media were deduced, respectively. By applying the given initial model to the forward calculation of Love wave, we compared dispersion curves of Love wave in the anisotropic media with the one in the isotropic media. the results show that, although the two kind of results are similar with each other, the effect of anisotropy can not be neglected. Furthermore, the variation of anisotropy factors will result in the variation of dispersion curves, especially for high-mode one. The method of grid dispersion inversion was then described for further tectonic inversion. We also deduced inversion equation on the condition that the layered media is anisotropic, and calculated the phase-velocity partial derivatives with respect to the model parameters, P- and S-wave velocities, density, anisotropic parameters for Rayleigh wave and Love wave. Having analyzed the results of phase-velocity partial derivatives, we concluded that the derivatives within each period decreased with the depth increasing, the phase-velocity of surface wave is sensitive to the S-wave velocities and anisotropic factors and is not sensitive to the densities of layers. Dispersion data of Love wave from the events occurred during the period from 1991 to 1998 around the Qinghai and Tibet Plateau, which magnitudes are more than 5.5, have been used in the grid dispersion inversion. Those data have been preprocessed and analyzed in the F-T domain. Then the results of 1°*1° grid dispersion inversion, the pure path dispersion data, in the area of Qianghai and Tibet Plateau were obtained. As an example, dispersion data have been input for the tectonic inversion in the anisotropic media, and the results of anisotropic factors under the region of Qianghai and Tibet Plateau were initially discussed. As for the other part of this dissertation. We first introduced the phenomena of S-wave splitting and the methods for calculation the splitting parameters. Then, We applied Butterworth band-pass filter to S-wave data recorded at 8 stations in China mainland, and analyzed S-wave splitting at different frequency bands. The results show the delay time and the fast polarization directions of S-wave splitting depend upon the frequency bands. There is an absence of S-wave splitting at the station of Wulumuqi (WMQ) for the band of 0.1-0.2Hz. With the frequency band broaden, the delay time of S-wave splitting decreases at the stations of Beijing (BJI), Enshi (ENH), Kunming (KMI) and Mudanjiang (MDJ); the fast polarization direction at Enshi (ENH) changes from westward to eastward, and eastward to westward at Hailaer (HIA). The variations of delay time with bands at Lanzhou (LZH) and qiongzhong (QIZ) are similar, and there is a coherent trend of fast polarization directions at BJI, KMI and MDJ respectively. Initial interpretations to the results of frequency band-dependence of S-wave splitting were also presented.

降雨触发滑坡的非线性模型与斜坡演化的混沌效应

This paper analyzes landsliding process by nonlinear theories, especially the influence mechanism of external factors (such as rainfall and groundwater) on slope evolution. The author investigates landslide as a consequence of the catastrophic slide of initially stationary or creeping slope triggered by a small perturbation. A fully catastrophe analysis is done for all possible scenarios when a continuous change is imposed to the control parameters. As the slip surface continues and erosion due to rainfall occurs, control parameters of the slip surface may evolve such that a previously stable slope may become unstable (e.g. catastrophe occurs), when a small perturbation is imposed. Thus the present analysis offers a plausible explanation to why slope failure occurs at a particular rainfall, which is not the largest in the history of the slope. It is found, by analysis on the nonlinear dynamical model of the evolution process of slope built, that the relationship between the action of external environment factors and the response of the slope system is complicatedly nonlinear. When the nonlinear action of slope itself is equivalent to the acting ability of external environment, the chaotic phenomenon appears in the evolution process of slope, and its route leading to chaos is realized with bifurcation of period-doublings. On the basis of displacement time series of the slope, a nonlinear dynamic model is set up by improved Backus generalized linear inversion theory in this paper. Due to the equivalence between autonomous gradient system and catastrophe model, a standard cusp catastrophe model can be obtained through variable substitution. The method is applied to displacement data of Huangci landslide and Wolongsi landslide, to show how slopes evolve before landsliding. There is convincing statistical evidence to believe that the nonlinear dynamic model can make satisfied prediction results. Most important of all, we find that there is a sudden fall of D, which indicates the occurrence of catastrophe (when D=0).

Lg尾波衰减成像研究——北京及邻区衰减成像

With the development of the technology of earthquake observation, more and more researchers work at many fields' of seismicity using seismic kinetic property, as the result, the study of attenuation has also made great progress, especially in the mechanism of the attenuation and the physical process. Aki put forward single back scattering theory to explain the forming of the seismic coda wave in 1969. Then, researchers started to develop the study in seismic scattering and attenuation. My thesis is also based on that theory. We assume that the Lg wave is a superposition fo high-mode surface waves, the coda of Lg is caused by scattering. Sato proposed Single Isotropic Scattering model (SIS model) to interpreted the scatter property, and he also formulated the geometrical spreading term. Then Xie (1988) developed the single spectral-ratio (SSR) method to obtain the Lg coda Q and the frequency dependent factor n. Later, he get to lateral images in the area of scatter ellipse. SSR method is explored and used in the study of Lg coda waves of regional earthquakes in my thesis. Choosing the earthquakes records with high ratio of signal-noise ,which were recorded at the stations from 1989 to 1999, we obtain the single trace Lg coda Q and its frequency dependent factor n. The results proved that SIS model is the reasonable model to explain the Lg coda wave, and SSR method also can be used to process Lg coda of regional earthquakes to get to the satisfied Lg coda Q. Based on the Lg coda Q we obtained using the former method, we explore the programs to inverse the regional Lg coda Q independently, and then make use of them to inverse the Lg coda Q of Beijing and adjacent area. The inversion result is satisfied. We conclude that the distribution of Qo (Q in lHz) is marked by the inhomogeneity, which is related to the tectonic structure: The value of Qo in uplift area, for example, Yanshan uplift, Taihang uplift, Luxi Uplift, is higher than the depression area, for example, Jizhong depression, Huanghua depression, and Jiyang depression, and the border between the higher Q area and lower Q area is very clear; Lg coda Q is also related to the velocity structure, higher velocity area is also with higher Q, lower velocity area is with lower Q; and higher heat-flow area is companied with lower Q. All in all, the value of Q reflects the difference of characteristics of lithofaces, porosity, the liquid content between the pores and heat flow. So, the Q value difference between uplift area and depression area reveals the difference of tectonic structure, lithology and physical character of the rock. So, the study of Lg coda Q is help to understand the earthquakes propagation mechanism through the inhomogenous medium, the cause of the coda, attenuation mechanism of the coda. Making use of the lateral images of Q, with velocity images, heat flow results, and other experimental result, we will be promoted to understand the complex structure of the crust, its inhomogenous character, and so on.

三维油藏数值模拟正反演算法研究

The dynamic prediction of complex reservoir development is one of the important research contents of dynamic analysis of oil and gas development. With the increase development of time, the permeabilities and porosities of reservoirs and the permeability of block reservoir at its boundaries are dynamically changing. How to track the dynamic change of permeability and porosity and make certain the permeability of block reservoir at its boundary is an important practical problem. To study developing dynamic prediction of complex reservoir, the key problem of research of dynamic prediction of complex reservoir development is realizing inversion of permeability and porosity. To realize the inversion, first of all, the fast forward and inverse method of 3-dimension reservoir simulation must be studied. Although the inversion has been widely applied to exploration and logging, it has not been applied to3-dimension reservoir simulation. Therefore, the study of fast forward and inverse method of 3-dimension reservoir simulation is a cutting-edge problem, takes on important realistic signification and application value. In this dissertation, 2-dimension and 3-dimension fluid equations in porous media are discretized by finite difference, obtaining finite difference equations to meet the inner boundary conditions by Peaceman's equations, giving successive over relaxation iteration of 3-dimension fluid equations in porous media and the dimensional analysis. Several equation-solving methods are compared in common use, analyzing its convergence and convergence rate. The alternating direction implicit procedure of 2-dimension has been turned into successive over relaxation iteration of alternating direction implicit procedure of 3-dimension fluid equations in porous media, which possesses the virtues of fast computing speed, needing small memory of computer, good adaptability for heterogeneous media and fast convergence rate. The geological model of channel-sandy reservoir has been generated with the help of stochastic simulation technique, whose cross sections of channel-sandy reservoir are parabolic shapes. This method makes the hard data commendably meet, very suit for geological modeling of containing complex boundary surface reservoir. To verify reliability of the method, theoretical solution and numerical solution are compared by simplifying model of 3-dimension fluid equations in porous media, whose results show that the only difference of the two pressure curves is that the numerical solution is lower than theoretical at the wellbore in the same space. It proves that using finite difference to solve fluid equations in porous media is reliable. As numerical examples of 3-dimension heterogeneous reservoir of the single-well and multi-well, the pressure distributions have been computed respectively, which show the pressure distributions there are clearly difference as difference of the permeabilities is greater than one order of magnitude, otherwise there are no clearly difference. As application, the pressure distribution of the channel-sandy reservoir have been computed, which indicates that the space distribution of pressure strongly relies on the direction of permeability, and is sensitive for space distributions of permeability. In this dissertation, the Peaceman's equations have been modified into solving vertical well problem and horizontal well problem simultaneously. In porous media, a 3D layer reservoir in which contain vertical wells and horizontal wells has been calculated with iteration. For channel-sandy reservoir in which there are also vertical wells and horizontal wells, a 3D transient heterogeneous fluid equation has been discretized. As an example, the space distribution of pressure has been calculated with iteration. The results of examples are accord with the fact, which shows the modification of Peaceman's equation is correct. The problem has been solved in the space where there are vertical and horizontal wells. In the dissertation, the nonuniform grid permeability integration equation upscaling method, the nonuniform grid 2D flow rate upscaling method and the nonuniform grid 3D flow rate upscaling method have been studied respectively. In those methods, they enhance computing speed greatly, but the computing speed of 3D flow rate upscaling method is faster than that of 2D flow rate upscaling method, and the precision of 3D flow rate upscaling method is better than that of 2D flow rate upscaling method. The results also show that the solutions of upscaling method are very approximating to that of fine grid blocks. In this paper, 4 methods of fast adaptive nonuniform grid upscaling method of 3D fluid equations in porous media have been put forward, and applied to calculate 3D heterogeneous reservoir and channel-sandy reservoir, whose computing results show that the solutions of nonuniform adaptive upscaling method of 3D heterogeneous fluid equations in porous media are very approximating to that of fine grid blocks in the regions the permeability or porosity being abnormity and very approximating to that of coarsen grid blocks in the other region, however, the computing speed of adaptive upscaling method is 100 times faster than that of fine grid block method. The formula of sensitivity coefficients are derived from initial boundary value problems of fluid equations in porous media by Green's reciprocity principle. The sensitivity coefficients of wellbore pressure to permeability parameters are given by Peaceman's equation and calculated by means of numerical calculation method of 3D transient anisotropic fluid equation in porous media and verified by direct method. The computing results are in excellent agreement with those obtained by the direct method, which shows feasibility of the method. In the dissertation, the calculating examples are also given for 3D reservoir, channel-sandy reservoir and 3D multi-well reservoir, whose numerical results indicate: around the well hole, the value of the sensitivity coefficients of permeability is very large, the value of the sensitivity coefficients of porosity is very large too, but the sensitivity coefficients of porosity is much less than the sensitivity coefficients of permeability, so that the effect of the sensitivity coefficients of permeability for inversion of reservoir parameters is much greater than that of the sensitivity coefficients of porosity. Because computing the sensitivity coefficients needs to call twice the program of reservoir simulation in one iteration, realizing inversion of reservoir parameters must be sustained by the fast forward method. Using the sensitivity coefficients of permeability and porosity, conditioned on observed valley erosion thickness in wells (hard data), the inversion of the permeabilities and porosities in the homogeneous reservoir, homogeneous reservoir only along the certain direction and block reservoir are implemented by Gauss-Newton method or conjugate gradient method respectively. The results of our examples are very approximating to the real data of permeability and porosity, but the convergence rate of conjugate gradient method is much faster than that of Gauss-Newton method.

开发早期油藏动态地质模型与油藏工程优化研究——以埕岛油田馆上段油藏为例

In order to developing reservoir of Upper of Ng at high-speed and high-efficient in Chengdao oilfield which is located in the bally shallow sea, the paper builds up a series of theory and means predicting and descripting reservoir in earlier period of oilfield development. There are some conclusions as follows. 1. It is the first time to form a series of technique of fine geological modeling of the channel-sandy reservoir by means of mainly seismic methods. These technique include the logging restriction seismic inversion, the whole three dimension seismic interpretation, seismic properties analysis and so on which are used to the 3-dimension distributing prediction of sandy body, structure and properties of the channel reservoir by a lot of the seismic information and a small quantity of the drilling and the logging information in the earlier stage of the oil-field development. It is the first time that these methods applied to production and the high-speed development of the shallow sea oilfield. The prediction sandy body was modified by the data of new drilling, the new reservoir prediction thinking of traced inversion is built. The applied effect of the technique was very well, according to approximately 200 wells belonging to 30 well groups in Chengdao oilfield, the drilling succeeded rate of the predicting sandy body reached 100%, the error total thickness only was 8%. 2. The author advanced the thinking and methods of the forecasting residual-oil prediction at the earlier stage of production. Based on well data and seismic data, correlation of sediment units was correlated by cycle-correlation and classification control methods, and the normalization and finely interpretation of the well logging and sedimentation micro-facies were acquired. On the region of poor well, using the logging restriction inversion technique and regarding finished drilling production well as the new restriction condition, the sand body distributing and its property were predicted again and derived 3-dimension pool geologic model including structure, reservoir, fluid, reservoir engineering parameter and producing dynamic etc. According to the reservoir geologic model, the reservoir engineering design was optimized, the tracking simulation of the reservoir numerical simulation was done by means of the dynamic data (pressure, yield and water content) of development well, the production rule and oil-water distributing rule was traced, the distributing of the remaining oil was predicted and controlled. The dynamic reservoir modeling method in metaphase of development was taken out. Based on the new drilling data, the static reservoir geologic model was momentarily modified, the research of the flow units was brought up including identifying flow units, evaluating flow units capability and establishing the fine flow units model; according to the dynamic data of production and well testing data, the dynamic tracing reservoir description was realized through the constant modification of the reservoir geologic model restricted these dynamic data by the theory of well testing and the reservoir numerical simulation. It was built the dynamic tracing reservoir model, which was used to track survey of the remaining oil on earlier period. The reservoir engineering tracking analysis technique on shallow sea oilfield was founded. After renewing the structure history since tertiary in Chengdao area by the balance section technique and estimating the activity character of the Chengbei fault by the sealing fault analysis technique, the meandering stream sediment pattern of the Upper of Ng was founded in which the meandering border was the uppermost reservoir unit. Based on the specialty of the lower rock component maturity and the structure maturity, the author founded 3 kinds of pore structure pattern in the Guanshang member of Chengdao oil-field in which the storing space mainly was primary (genetic) inter-granular pore, little was secondary solution pore and the inter-crystal pore tiny pore, and the type of throat mainly distributed as the slice shape and the contract neck shape. The positive rhythmic was briefly type included the simple positive rhythm, the complex positive rhythm and the compound rhythm. Interbed mainly is mudstone widely, the physical properties and the calcite interbed distribute localized. 5. The author synthetically analyzed the influence action of the micro-heterogeneity, the macro-heterogeneity and the structure heterogeneity to the oilfield water flood development. The efficiency of water flood is well in tiny structure of convex type or even type at top and bottom in which the water breakthrough of oil well is soon at the high part of structure when inject at the low part of structure, and the efficiency of water flood is poor in tiny structure of concave type at top and bottom. The remaining oil was controlled by sedimentary facies; the water flooding efficiency is well in the border or channel bar and is bad in the floodplain or the levee. The separation and inter layer have a little influence to the non-obvious positive rhythm reservoir, in which the remaining oil commonly locate within the 1-3 meter of the lower part of the separation and inter layer with lower water flooding efficiency.

复杂断块油藏形成机制及剩余油预测

The central uplift in the Huimin depression is famous for its large amounts of faults and small-scale fault-block area, and it is the famed typical complicated fault-block group oil & gas field in the whole world. After many years of rolling exploration and exploitation, many complex oil &gas field have been discovered in the central uplift, and won the splendent fruit. With the gradual deepening and development of the rolling exploitation, the exploration faces more and more difficulties. Therefore, it is important to reveal the forming mechanism and distributing rule of the complex fault-block reservoir, and to realize the forecast of the complex fault-block reservoir, sequentially, expedite the exploration step. This article applies the new multi-subject theory, method and technique such as structure geometry, kinematics, dynamics, structural stress field, fluid potential field, well logging record and constrained inversion of seismic records, coherence analysis, the seal mold and seal history of oil-bounded fault etc, and try to reveal the forming mechanism and distributing law of the complex fault-block reservoir, in result, implements the forecast of the fault-block reservoir and the remaining oil distributing. In order to do so, this article synthetically carries out structural estimate, reservoir estimate, fault sealing history estimate, oil-bearing properties estimate and residual. This article also synthetically researches, describes and forecast the complex fault-block in Huimin depression by use of the techniques, e.g. seismetic data post-stack processing technique, multi-component demarcating technique, elaborate description technique for the fault-block structure, technique of layer forecasting, fault sealing analysis technique, comprehensive estimate technique of fault-block, comprehensive analysis and estimate technique of remaining oil etc. The activities of the faults varies dramatically in the Huimin depression, and most of the second-class and the third-class faults are contemporaneous faults, which control the macroscopical distribution of the reservoir in the Huimin depression. The fourth-class faults cause the complication between the oil & gas among the fault-blocks. The multi-period strong activities of the Linyi fracture resulted in the vertical migration of large amount of oil & gas along with the faults. This is the main reason for the long vertical distribution properties near the Linyi fracture in the Huimin depression. The sealing ability of the fault is controlled by the property,size and direction of the main stress, the contact relationship of the both sides of the fault, the shale polluting factor, and the configuration relationship between the fault move period and the migration period of oil & gas. The article suggest four fault-sealing modes in the research zone for the first time, which establishes the foundation for the further forecast of the complex fault-block reservoir. Numerical simulation of the structural stress field reveals the distribution law and the evolvement progress of the three-period stress field from the end of the Dongying period to the Guantao period to nowadays. This article puts forward that the Linyi and Shanghe regions are the low value of the maximum main stress data. This is combined with the fault sealing history estimate, then multi-forming-reservoir in the central uplift is put forward. In the Shanghe oilfield, the article establishes six reservoir geological modes and three remaining oil distributing modes(the plane, the inside layer and the interlayer), then puts forward six increase production measure to enhance the remaining oil recovery ratio. Inducting the exploitation of oilfield, it wins notable economic effects and social effects.

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

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.

陡坡带砂砾岩体油气成藏模式研究——以东营凹陷西北部陡坡带为例

Based on the study of sequence stratigraphy, modern sedimentary, basin analysis, and petroleum system in abrupt slop of depression, this paper builds sedimentary system and model, sandy bodies distribution, and pool-forming mechanism of subtle trap. There are some conclusions and views as follows. By a lot of well logging and seismic analysis, the author founded up the sequence stratigraphic of the abrupt slope, systematically illustrated the abrupt slope constructive framework, and pointed out that there was a special characteristics which was that south-north could be divided to several fault block and east-west could be carved up groove and the bridge in studying area. Based all these, the author divided the studying area to 3 fault block zone in which because of the groove became the basement rock channel down which ancient rivers breathed into the lake, the alluvial fan or fan delta were formed. In the paper, the author illustrated the depositional system and depositional model of abrupt slope zone, and distinguished 16 kinds of lithofacies and 3 kinds of depositional systems which were the alluvial fan and fan-delta system, lake system and the turbidite fan or turbidity current deposition. It is first time to expound completely the genetic pattern and distributing rule of the abrupt slope sandy-conglomeratic fan bodies. The abrupt slope sandy-conglomeratic fan bodies distribute around the heaves showing itself circularity shape. In studying area, the sandy-conglomeratic fan bodies mainly distribute up the southern slope of Binxian heave and Chenjiazhuang heave. There mainly are these sandy-conglomeratic fan body colony which distributes at a wide rage including the alluvial fan, sub-water fluvial and the turbidite fan or the other turbidity current deposition in the I fault block of the Wangzhuang area. In the II fault block there are fan-delta front and sub-water fluvial. And in the Binnan area, there mainly are those the alluvial fan (down the basement rock channel) and the sandy-conglomeratic fan body which formed as narrowband sub-water fluvial (the position of bridge of a nose) in the I fault block, the fan-delta front sandy-conglomeratic fan body in the H fault block and the fan-delta front and the turbidity current deposition sandy-conglomeratic fan body in the m fault block. Based on the reservoir outstanding characteristics of complex classic composition and the low texture maturity, the author comparted the reservoir micro-structure of the Sha-III and Sha-IV member to 4 types including the viscous crude cementation type, the pad cementation type, the calcite pore-funds type and the complex filling type, and hereby synthetically evaluated 4 types sandy- conglomeratic fan body reservoir. In the west-north abrupt slope zone of Dongying Depression, the crude oil source is belonging to the Sha-III and Sha-IV member, the deep oil of Lijin oilfield respectively come from the Sha-III and Sha-IV member, which belongs to the autogeny and original deposition type; and the more crude oil producing by Sha-IV member was migrated to the Wangzhuan area and Zhengjia area. The crude oil of Binnan oil-field and Shanjiasi oil-field belongs to mixed genetic. It is the first time to illustrate systematically the genetic of the viscous crude that largely being in the studying area, which are that the dissipation of the light component after pool-forming, the biological gradation action and the bath-oxidation action, these oil accumulation belonging to the secondary viscous crude accumulation. It is also the first time to compart the studying area to 5 pool-forming dynamical system that have the characteristic including the common pressure and abnormal pressure system, the self-fountain and other-fountain system and the closing and half-closing system etc. The 5 dynamical systems reciprocally interconnected via the disappearance or merger of the Ethology and the fluid pressure compartment zone, the fault and the unconformity surface, hereby formed duplicated pattern oil-gas collecting zone. Three oil-gas pool-forming pattern were founded, which included the self-fountain side-direction migrated collecting pattern, the self-fountain side-direction ladder-shape pool-forming pattern and the other-fountain pressure releasing zone migrated collecting pattern. A series of systemic sandy-conglomeratic fan bodies oil-gas predicting theory and method was founded, based on the groove-fan corresponding relation to confirm the favorable aim area, according as the characteristic of seismic-facies to identify qualitatively the sandy-conglomeratic fan bodies or its scale, used the temporal and frequency analysis technique to score the interior structure of the sandy- conglomeratic fan bodies, applied for coherent-data system analysis technology to describe the boundary of the sandy-conglomeratic fan bodies, and utilized the well logging restriction inversion technique to trace quantificational and forecast the sandy-conglomeratic fan bodies. Applied this technique, totally 15 beneficial sandy-conglomeratic fan bodies were predicted, in studying area the exploration was preferably guided, and the larger economic benefit and social benefit was acquired.

超参数控制下的波阻抗约束反演

The seismic survey is the most effective geophysical method during exploration and development of oil/gas. As a main means in processing and interpreting seismic data, impedance inversion takes up a special position in seismic survey. This is because the impedance parameter is a ligament which connects seismic data with well-logging and geological information, while it is also essential in predicting reservoir properties and sand-body. In fact, the result of traditional impedance inversion is not ideal. This is because the mathematical inverse problem of impedance is poor-pose so that the inverse result has instability and multi-result, so it is necessary to introduce regularization. Most simple regularizations are presented in existent literature, there is a premise that the image(or model) is globally smooth. In fact, as an actual geological model, it not only has made of smooth region but also be separated by the obvious edge, the edge is very important attribute of geological model. It's difficult to preserve these characteristics of the model and to avoid an edge too smooth to clear. Thereby, in this paper, we propose a impedance inverse method controlled by hyperparameters with edge-preserving regularization, the inverse convergence speed and result would be improved. In order to preserve the edge, the potential function of regularization should satisfy nine conditions such as basic assumptions edge preservation and convergence assumptions etc. Eventually, a model with clear background and edge-abnormity can be acquired. The several potential functions and the corresponding weight functions are presented in this paper. The potential functionφLφHL andφGM can meet the need of inverse precision by calculating the models. For the local constant planar and quadric models, we respectively present the neighborhood system of Markov random field corresponding to the regularization term. We linearity nonlinear regularization by using half-quadratic regularization, it not only preserve the edge, and but also simplify the inversion, and can use some linear methods. We introduced two regularization parameters (or hyperparameters) λ2 and δ in the regularization term. λ2 is used to balance the influence between the data term and the transcendental term; δ is a calibrating parameter used to adjust the gradient value at the discontinuous position(or formation interface). Meanwhile, in the inverse procedure, it is important to select the initial value of hyperparameters and to change hyperparameters, these will then have influence on convergence speed and inverse effect. In this paper, we roughly give the initial value of hyperparameters by using a trend- curve of φ-(λ2, δ) and by a method of calculating the upper limit value of hyperparameters. At one time, we change hyperparameters by using a certain coefficient or Maximum Likelihood method, this can be simultaneously fulfilled with the inverse procedure. Actually, we used the Fast Simulated Annealing algorithm in the inverse procedure. This method overcame restrictions from the local extremum without depending on the initial value, and got a global optimal result. Meanwhile, we expound in detail the convergence condition of FSA, the metropolis receiving probability form Metropolis-Hasting, the thermal procession based on the Gibbs sample and other methods integrated with FSA. These content can help us to understand and improve FSA. Through calculating in the theoretic model and applying it to the field data, it is proved that the impedance inverse method in this paper has the advantage of high precision practicability and obvious effect.

接收函数方法与中国及邻区上地幔间断面的研究

The receiver function method applied in researching the discontinuities in upper mantle was systematically studied in this paper. Using the theoretical receiver functions, the characteristics of P410S and P660S phases were analyzed, and the influencing factors for detection of these phases were discussed. The stability of receiver function was studied, and a new computational method of receiver function, RFSSMS (Receiver Function of Stack and Smooth of Multi seismic-records at a Single station), was put forward. We built initial reference velocity model for the media beneath each of 18 seismic stations respectively; then estimated the buried depths of 410-km and 660-km discontinuities(simply marked as '410' and '660') under the stations by using the arrive time differences of P410S and P660S with P. We developed a new receiver function inversion method -PGARFI (Peeling-Genetic Algorithm of Receiver Function Inversion), to obtain the whole crust and upper mantle velocity structure and the depths of discontinuities beneath a station. The major works and results could be summarized as follows: (1) By analysis of the theoretical receiver functions with different velocity models and different ray parameters, we obtain the knowledge: The amplitudes of P410S and P660S phases are decreasing with the increasing of epicentral distance A , and the arrival time differences of these phases with P are shorter as A is longer. The multiple refracted and/or reflected waves yielded on Moho and the discontinuities in the crust interfere the identification of P410S. If existing LVZ under the lithosphere, some multiple waves caused by LVZ will interfere the identification of P410S. The multiple waves produced by discontinuity lied near 120km depth will mix with P410s phase in some range of epicentral distance; and the multiple waves concerned with the discontinuity lied near 210km depth will interfere the identification of P660S. The epicentral distance for P4i0s identification is limited, the upper limit is 80° . The identification of P660S is not restricted by the epicenter distance obviously. The identification of P410S and P6gos in the theoretical receiver functions is interfered weakly from the seismic wave attenuation caused by the media absorption if the Q value in a reasonable range. (2) The stability of receiver function was studied by using synthetic seismograms with different kind of noise. The results show that on the condition of high signal-noise-ratio of seismic records, the high frequency background noise and the low frequency microseism noise do not influence the calculating result of receiver function. But the media "scattering noise" influence the stability of receiver function. When the scattering effect reach some level, the identification of P4iOs and P66os is difficult in single receiver function which is yielded from only one seismic record. We provided a new method to calculate receiver function, that is, with a group of earthquake records, stacking the R and Z components respectively in the frequency domain, and weighted smooth the stacked Z component, then compute the complex spectrum ratio of R to Z. This method can improve the stability of receiver function and protrude the P4i0s and P66os in the receiver function curves. (3) 263 receiver functions were provided from 1364 three component broadband seismograms recorded at 18 stations in China and adjacent areas for the tele-earthquakes. The observed arrival time differences of P410S and P660S with P were obtained in these receiver functions. The initial velocity model for every station was built according to the prior research results. The buried depths of '410' and '660' under a station were acquired by the way of adjusting the depths of these two discontinuities in the initial velocity model until the theoretical arrival time differences of P410S and P660S with P well conformed to the observed. The results show an obvious lateral heterogeneity of buried depths of ' 410' and (660' . The depth of '410' is shallower beneath BJI, XAN, LZH and ENH, but deeper under QIZ and CHTO, and the average is 403km . The average depth of '660' is 663km, deeper under MDJ and MAJO, but shallower under QIZ and HYB. (4) For inversing the whole crust and upper mantle velocity structure, a new inversion method -PGARFI (Peeling-Genetic Algorithm of Receiver Function Inversion) has- been developed here. The media beneath a station is divided into segments, then the velocity structure is inversed from receiver function from surface to deep successively. Using PGARFI, the multi reflection / refraction phases of shallower discontinuities are isolated from the first order refraction transform phase of deep discontinuity. The genetic algorithm with floating-point coding was used hi the inversion of every segment, and arithmetical crossover and non-uniform mutation technologies were employed in the genetic optimization. 10 independent inversions are completed for every segment, and 50 most excellent velocity models are selected according to the priority of fitness from all models produced in the inversion process. The final velocity structure of every segment is obtained from the weighted average of these 50 models. Before inversion, a wide range of velocity variation with depth and depth range of the main discontinuities are given according to priori knowledge. PGARFI was verified with numerical test and applied in the inversion of the velocity structure beneath HIA station down to 700km depth.

地震波散射非线性反演问题不动点理论及方法的研究

How to create a new method to solve the problem or reduce the influence of that the result of the seismic waves scattering nonlinear inversion is not uniqueness is a main purpose of this research work in the paper. On the background of research into the seismic inversion, new progress of the nonlinear inversion is introduced at the first chapter in this paper. Especially, the development, basic theories and assumptions on some major theories of seismic inversion are analyzed, discussed and summarized in mathematics and physics. Also, the problems faced by the mathematical basis of investigations of the seismic inversion are discussed, and inverse questions of strongly seismic scattering due to strong heterogeneous media in the Earth interior are analyzed and viewed. What the kernel of paper is that gathers all our attention making a new nonlinear inversion method of seismic scattering. The paper provides a theory and method of how to introduce the fixed-point theory into the nonlinear seismic scattering inversion and how to obtain the solution, and gives the actually method to create a serials of contractive mappings of velocity parameter's in the mapping space of wave. Therefore, the results testify the existence of fixed point of velocity parameter and give the method the find it. Further, the paper proves the conclusion that the value obtained by taking the fixed point of velocity parameter into wave equation is the fixed point of the wave of the contractive mapping. Thence, the fixed point is the global minima since the stabilities quality of the fixed point. Based on the new theory, in the chapter three, many inverse results are obtained in the numerical value test. By analysis the results one could find a basic facts that all the results, which are inversed by the different initial model, are tended to the true value in theoretical true model. In other words, the new method mostly eliminates the non-uniqueness that which is existed in seismic waves scattering nonlinear inversion in degree. But, since the test results are quite finite now, more test is need here to positive our theory. As a new theoretical method, it must be existed many weaken in it. The chapter four points out all the questions which is bother us. We hope more people to join us to solve the problem together.

天然地震转换波偏移成像方法研究

Butovskaya, a scholar of Former Soviet Union, first determined the depth of basalt layer in Tashkent Zone by using converted waves on seismogram in 1952. From then on, more and more scholars developed the comprehensive research that imaged the earth interior structures by applying converted waves information. With the digitalization of earthquake observation, The inversion imaging of complete or partial waveform record can efficiently improve inversion quality and widen its usage scope, therefore great progress is made in converted wave imaging. This paper makes a certain study in converted wave imaging on that basis. Transmitted PP waves and converted PS waves are generated when a P-wave propagates through an interface separating two media with large impedance contracts. A PS converted wave is a seismic body wave, which result from the conversion of an incident parent P wave at a boundary within the crust to a refracted S wave. The thickness of a single crustal layer can theoretically be determined by observing, with three-componented seismometer at a single station, the difference in time of the arrival of the parent P wave and the arrival of the PS converted wave. For a multilayered media, PS converted wave arrivals corresponding to each of the layers can theoretically be observed, provided the station is sufficiently from the source of the parent P wave to allow initial penetration of the P wave beneath the deepest layer considered. To avoid the difficulty of picking up transmitted P-wave and converted wave phases, this paper proposed a converted wave migration method by estimating the travel time difference between PS converted wave and PP transmitted wave. To verify its validity, we apply the converted wave PS migration algorithm to synthetic data generated by three forward modeling. The migration results indicate that PS converted wave may be migrated to reconstruct the transmitting interface. This technique is helpful to investigate the deep earth structures by using earthquake data.

指数矩阵近似计算在地震波场延拓中的应用

Seismic wave field numerical modeling and seismic migration imaging based on wave equation have become useful and absolutely necessarily tools for imaging of complex geological objects. An important task for numerical modeling is to deal with the matrix exponential approximation in wave field extrapolation. For small value size matrix exponential, we can approximate the square root operator in exponential using different splitting algorithms. Splitting algorithms are usually used on the order or the dimension of one-way wave equation to reduce the complexity of the question. In this paper, we achieve approximate equation of 2-D Helmholtz operator inversion using multi-way splitting operation. Analysis on Gauss integral and coefficient of optimized partial fraction show that dispersion may accumulate by splitting algorithms for steep dipping imaging. High-order symplectic Pade approximation may deal with this problem, However, approximation of square root operator in exponential using splitting algorithm cannot solve dispersion problem during one-way wave field migration imaging. We try to implement exact approximation through eigenfunction expansion in matrix. Fast Fourier Transformation (FFT) method is selected because of its lowest computation. An 8-order Laplace matrix splitting is performed to achieve a assemblage of small matrixes using FFT method. Along with the introduction of Lie group and symplectic method into seismic wave-field extrapolation, accurate approximation of matrix exponential based on Lie group and symplectic method becomes the hot research field. To solve matrix exponential approximation problem, the Second-kind Coordinates (SKC) method and Generalized Polar Decompositions (GPD) method of Lie group are of choice. SKC method utilizes generalized Strang-splitting algorithm. While GPD method utilizes polar-type splitting and symmetric polar-type splitting algorithm. Comparing to Pade approximation, these two methods are less in computation, but they can both assure the Lie group structure. We think SKC and GPD methods are prospective and attractive in research and practice.

凝析油气藏描述及开发对策研究-以莫北2油气藏为例

Gas condensate reservoir research involves not only structure sediment reservoir liquid properties characterization but also the change of the temperature field, the change of the pressure field, the change of liquid phase and the reservoir sensitivity. To develop the gas condensate reservoir effectively .we must depict the static properties of the oil and gas system ,build exact and comprehensive parameter field, predict the rule of dynamic change and do the necessary reservoir characterization development plan dynamic prediction direct production. The MoBei Oil and Gas Field is the first gas condensate reservoirs which is found by the Xinjiang Oil Field Company in ZhunGaEr basin belly.it has deserved some knowledge after prospect evaluation, the MoBei Oil and Gas Field start development ,it is one of the important development blocks of Xinjiang Oil Field Company productivity constuction. During its development , it gradually appears some problems, such as complex oil and gas phase, great change of reservoir stretch .uncertain reservoir type and scale, controling its development strategy and plan difficultly. To deserve the high efficient development and long-term stable production of the gas condensate reservoir, it is necessary to characterize it systematically and form a suit of scientific development strategy. This thesis take the MoBei zone SanGongHe sand group reservoir as research object, applied advanced log techniques ,such as the nulear magnetism log ,MDT testing .etc. After comprehensive research of loging geology information, set up a suit of methods to identify oil gas water layer .these methods can identify the gas-oil level and the oil-water level. On the basis of reasonable development object system, according fine structure interpretation and structure modeling. build any oil water column height of the reservoir accurately. Through carefully analysis of the basic theory and method of reservoir seism prediction. optimize a reservoir inversion method .technique. software fitting the research region aiming strata, set up the GR field, porosity field, Rt field, impedence field .permeability field and initial oil saturation field, generating the base of quantity reservoir characterization. Discussing the characteristic of reservoir fluid and the movement and reallocating of muti-phase fluid in reservoir. And according the material of 100 soviet gas condensate reservoir ,build the recognition method and mode of gas condensate reservoir. Building the 3D geology model ,carry on the static and production evaluation, propose the development strategy and improve plan , provide the base of increasing reserves and advancing production and enriching the prospect development theory of the gas condensate reservoi