准噶尔盆地三个泉凸起油气成藏条件及动力学研究

Sangequan Uplift in Junggar Basin is an inherited positive structure, which has undergone many times of violent tectonic movements, with high tectonic setting, and far away from the oil-source sag, reservoir forming condition is complex. Combining sequence stratigraphy, depositional facies, reservoir formation theory with seismic and well logging analysis, this paper conducted integrated study on the hydrocarbon migration, accumulation, entrapment conditions, the reservoir forming dynamics and the forming model, and acquired the following recognition: (1) The special reservoir formation conditions that enable Sangequan Uplift to form a giant oil-gas field of over 100 million tons of reserves are as follows: (D Deltaic frontal sandbody is developed in Jurassic Xishanyao Formation, Toutunhe Formation and Lower Cretaceous Hutubihe Formation, with good reservoir quality;? Abundant hydrocarbon resources are found in Western Well Pen-1 Sag, which provides sufficient oil sources for reservoir formation of Sagequan Uplift; ?The unconformity-fault-sandbody system has formed a favorable space transporting system and an open conduit for long-distance hydrocarbon migration; ?fault, low amplitude anticline and lithological traps were well developed, providing a favorable space for hydrocarbon accumulation. (2) The most significant source beds in the Western Well Pen-1 Sag are the Mid-Permian Lower Wuerhe Formation and Lower-Permian Fengcheng Formation. The oil in the Well Block Lu-9 and Shinan Oilfield all originated from the hydrocarbon source beds of Fengcheng Formation and Lower Wuerhe Formation in the Western Well Pen-1 Sag and migrated through Jidong and Jinan deep faults linking unconformity of different regions from sources to structural highs of the uplift and shallow horizons. (3) There were 2 reservoir formation periods in District Sangequan: the first was in late Cretaceous during which the upper part of Xishanyao Formation and Toutunhe Formation; the second was in Triassic, the main resources are high-maturity oil and gas from Fengcheng Formation and Wuerhe Formation in Western Pen-1 Well sag and the gas from coal measure strata of Xishanyao Formation, that were accumulated in Hutubihehe Formation. (4) Model of the hydrocarbon migration, accumulation, reservoir formation of the study area are categorized as three types starting from the hydrocarbon source areas, focusing on the faults and unconformity and aiming at reservoirs: ① Model of accumulation and formation of reservoir through faults or unconformities along the "beam" outside source; ②Model of migration, accumulation and reservoir formation through on-slope near source;③Model of migration, accumulation and reservoir formation of marginal mid-shallow burial biogas-intermediate gas. (5) Pinchout, overlap and lithologic traps are developed in transitional zones between Western Well Pen-1 sag and Luliang uplift. Many faulted blocks and faulted nose-like traps are associated with large structures on Sangequan uplift. Above traps will be new prospecting areas for further hydrocarbon exploration in future.

论大型拆离滑覆构造的特征、形成机制及其与油气成藏关系—以渤海湾盆地冀中坳陷饶阳凹陷为例

The studies of this paper is an important part of the "ninth five" science&technology-tackling project of CNPC -The oil and gas distribution regulation and the aims of explortion in jizhong depression. Basing on the former research results, with the materials of regional structural setting, major tectonic movements, bi-and tri-dimension seismic sections, oil well sections and reservoir sections, this paper involves studies of tectonic evolution, sedimentarv evolution, magma movement and reservoir prediction. The existence of huge stripping and gliding nappe is proved in the RaoYang Sag for the frist time. The properties, development, evolution and the relationship with reservoir of the stripping and gliding nappe are discussed in details in this paper. It is also talked about the affects of stripping and nappes to oil and gas exploration theoretically and practically in the paper. The marking attributes of the stripping and gliding nappe includes stripping and gliding plane, two deformation systems, stratigraphic repeat and hiatus close to the stripping and gliding plane, and the deformation attributes in the front and back of stripping and gliding nappe. The RaoYang stripping and gliding nappes can be divided into different belts in north-south direction and different zones in east-west direction. RaoYang Stripping and gliging nappes took place in the late Paleogene period and before the sedimentation of Neogene period. The sliding direction is NWW. The sliding distance is about 6km. The geothermal gradient in the separating slump area is low and stable. The formation of the stripping and gliding nappes is due to the regional structural setting, the sediments of Paleogene system, the soft roof and the uneven rising movement of structure units. The evolution of the stripping and gliding nappes can be divided into the following stages: regional differential elevation and subsidence, unstable gravity and gravitational sliding, the frist wholly stripping faults and sliding stage, and the following second and third stripping faults and sliding stages. The identification of RaoNan stripping and gliding nappes has an important role on the research of regional structure and oil and gas exploration. Basing on the properties of stripping and gliding nappes, we can identtify the gliding fractures, ductile compressional folds, the front and back structures of gliding nappes and gliding plane covered structures. Combination with different reservoir forming conditions, these structures can lead to different categories of reservoirs.

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

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.

辽河盆地东部凹陷热历史及其与油气关系

This dissertation focuses on the basin geothermal history, tectonothermal evolution and the relationship between geothermal field evolution and hydrocarbon generation. Based on the research of present-day geothermal field, geothermal history of Eastern Subdepression of Liaohe Basin was reconstructed with available data from drillings, loggings, seismic cross-sections, BHTs and thermal indicators. 12 heat flow density values were calculated. Ranging from 45.7 mW/m~2 to 70.0 mW/m~2, the mean value of these determinations exhibits 58.0(±5.83mW/m~22). The heat flow density in the uplift and ramp is greater than that in the sag. The main factors affecting the heat flow density are the morphostructure of basement and thickness of sedimentary cover. The Subdepression experienced a much higher heat flow period from 43 Ma to about 25 Ma. The heat flow increased gradually from Sha3 stage to Dongying stage, and reached the peak at the late of Dongying stage, then cooled down. Structural subsidence analysis shows that the subsidence of Eastern Subdepression can be divided into two phases: earlier (25-43Ma) initial (rift) and late (since 25Ma) thermal subsidence. The lower present-day heat flow and the higher palaeo-heat flow corresponding to structural subsidence stage as well as the typical rift subsidence style in Eastern Subdepression provide with some insights to the tectonic-thermal evolution of the basin. The source rocks of Sha3 member began to generate oil in the Shal stage, and entered oil-window at the late of Dongying stage. The source rocks of Shal began to generate oil at the late of Dongying stage, and being at the stage of lower maturation-maturation now. Most of Dongying source rocks are not mature now. The late of Dongying stage is the critical time for the oil system.

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

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.

微功耗海底地震仪的研制

Our motherland has large area of maritime space. Searching and developing ocean becomes more and more important. So Ocean Bottom Seismometer (OBS) as an absolutely necessary equipment can be used in many oceanic fields. OBS not only is an important instrument for discovering structure of lithosphere of ocean bottom, but also plays a main role of oceanic geophysical exploration. The paper introduces my relational work. The MCI micro-power broad frequency seismometer was developed independently. Its power dissipation is less than 300mW. It has some merits including miniature volumeN light mass and cheap price. It is an ideal device not only for the collection high-resolution natural seismic data, but also for the fields of seismic sounding and engineering seismology. Many new high technique were applied to develop this instrument including over-sampling A/D converter, high performance 32bit Micro Process Unit and Flash memory with smart-media interface. Base on the achievement, I have accomplished the showpiece of OBS, which is applied to the deepwater oil and gas geophysical exploration. Because of micro-power dissipation, the seismograph and the sonar releaser can be integrated into a sphere cabin. By this means, the instrument's frequency of resonance and frequency of couple are improved obviously. The data acquisition system of OBS is improved from MCI seismometer. The capacity of flash memory is enlarge from 1G bytes to 8G bytes. The advance MPU in data acquisition system is used to integrate other function modules such as sonar, GPS, compass and digital transmitter.

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

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.

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

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.