扩散波勘探方法基础及其在石油油气藏开发中的应用研究

As we know, the essence of exploration is objective body determined by getting the information. Such as seismic、electrical and electromagnetic prospecting, they are the common methods of the exploration. Therefore, They have a complete set of theory now. In fact, the effective information can also be got by the diffusion way, it is called diffusion prospecting. The diffusion way prospecting is necessary and important. The way of diffusion prospecting is studied in the paper and main works include below: (1) On the basis of studying basic law of the diffusion, the paper gives the idea of diffusion wave and the formulas of computing diffusion wave function. (2) The paper studies the way of the diffusion prospecting and the methods of data processing. At the same time, it also expounds the characteristics and the applied foreground of the diffusion prospecting. (3) The paper gives the tomography idea and the basic method of diffusion CT. Meanwhile, it also expounds the foreground that the diffusion CT is applied in oil development prospecting. (4) As the inversion of the diffusion equation is a part of the diffusion prospecting way, the methods of diffusion equation inversion are studied and the two formulas are deduced --Laplace transform and polynomial fitting inversion formulas. As the other important result of diffusion equation inversion, the inversion can offer a new analysis method for well Testing in oil development. In order to show a set of methods in the paper feasible, forward、inversion and CT numerical simulation are done in the paper.

沾化凹陷火成岩油气藏形成条件和分布规律研究

As a kind of special lithologic ones, Igneous rock oil and gas pool is more and more paid attention, and it has different forming condition and distribution from conventional ones, such as various terrane distribution types, serious reservoir anisotropy, complicated hydrocarbon-bearing, so there is not successful experience to follow for exploration and development of this complex subtle oil and gas pool at present. For an example of Igneous oil and gas pool of Luo151 area in Zhanhua seg, Eastern China, this article study the difficult problem, including petrologic nd lithofacies analysis, Origin, invasion age and times of Igneous rock, reservoir anisotropy, Geological Modeling, Igneous reservoir synthesis evaluation. forming condition and distribution are studied synthetically, and an integrated method to predict igneous rock oil and gas pool is formed, which is evaluated by using development data. The Igneous rock is mainly diabase construction in Luo151 area of Zhanhua Sag, and petrologic types include carbonaceous slate, hornfels, and diabases. Based on analyzing synthetically petrologic component, texture and construct, 4 lithofacies zones, such as carbonaceous slate subfacies, hornfels subfacies containing cordierite and grammite, border subfacies and central subfacies, are divided in the diabase and wall rock. By studying on isotopic chronology, terrane configuration and imaging logging data, the diabase intrusion in Zhanhua Sag is formed by tholeiite magma emplacing in Shahejie formation stratum on the rift tension background Lower Tertiary in North China. The diabase intrusion of Luo151 is composed possibly of three periods magma emplacement. There is serious anisotropy in the diabase reservoirs of Luo151 in Zhanhua Sag. Fracture is primary reservoir space, which dominated by tensile fracture in high obliquity, and the fracture zones are mainly developed round joint belt of igneous rock and wall rock and position of terrane thickness changing rapidly. The generation materials of the reservoirs in Luo151 igneous oil pools consist of Intergranular micropore hornfels, condensate blowhole-solution void diabase condensate edge, the edge and center of the condensate seam diabase, of which are divided into horizontal, vertical and reticulated cracks according fracture occurrence. Based on the above research, a conceptual model of igneous rock reservoir is generated, which is vertically divided into 4 belts and horizontally 3 areas. It is built for the first time that classification evaluation pattern of igneous rock reservoir in this area, and 3 key wells are evaluated. The diabase construction is divided into grammite hornfels micropore type and diabase porous-fracture type reservoirs. The heavy mudstone layers in Third Member of Shahejie formation (Es3) provide favorable hydrocarbon source rock and cap formation, diabase and hornfels belts serve as reservoirs, faults and microcracks in the wall rocks as type pathways for oil and gas migration. The time of diabase invasion was about in the later deposition period of Dongying Formation and the middle of that of Guantao Formation, the oil generated from oil source rock of Es3 in the period of the Minghuazhen formation and is earlier more than the period of diabase oil trap and porous space forming. Based on geological and seismic data, the horizon of igneous rocks is demarcated accurately by using VSP and synthetic seismogram, and the shape distribution and continuity of igneous rocks are determined by using cross-hole seismic technology. The reservoir capability is predicted by using logging constraining inversion and neural network technology. An integrated method to predict igneous rock oil and gas pool is formed. The study is appraised by using development data. The result show the reservoir conceptual model can guide the exploration and development of oil pool, and the integrated method yielded marked results in the production.

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

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

南海北部大陆边缘盆地热演化史及其与油气关系

This thesis focuses on the present-day thermal field features, evolution and their connections to hydrocarbon generation of the three continental margin basins-the Yinggehai (Yingge Sea), Qiongdongnan(southeast Qiong), and Pear River Mouth basins-in northern South China Sea, based on available data from drillings, loggings, seismic cross-sections, BHTs, thermal indicators (Ro%, inclusion, etc) and geopressure measurements. After studying of present-day distribution of geothermal field and thermal disturbance of fluid in the sedimentary strata, the author discovered that the distribution of gas fields in Yinggehai Basin are closely related to the distribution of anomalously high thermal gradient area, whereas it is not the case for the Pear River Mouse Basin. And detailed processing of the fluid inclusion data indicates that geothermal fluids activated frequently in this area, and they may mainly be derived upward from the overpressure and hydrocarbon-generating beds, 3000-4500 m in depth. Therefore, the abnormal gradients in sedimentary beds were mainly caused by the active geothermal fluids related to hydrocarbon migrating and accumulating in this area. Because of the effect of overpressure retarding on vitrinite reflectance, the thermal indicators for thermal history reconstruction should be assessed before put into use. Although some factors, such as different types of kerogen, heating ratio, activities of thermal fluids and overpressure, may have effects on the vitrinite reflectance, under the circumstance that thermal fluids and overpressure co-exist, overpressure retarding is dominant. And the depth and correction method of overpressure retarding were also determined in this paper. On the basis of reviewing the methods of thermal history studies as well as existing problems, the author believes that the combination of thermal-indicator-inversion and tectono-thermal modeling is an effective method of the thermal history reconstruction for sedimentary basins. Also, a software BaTherMod for modeling thermal history of basins was successfully developed in this work. The Yinggehai Basin has been active since Tertiary, and this was obviously due to its tectonic position-the plate transition zone. Under the background of high thermal flow, long-term quick subsidence and fluid activities were the main reasons that lead to high temperature and overpressure in this basin. The Zhujiangkou Basin, a Tertiary fault-basin within the circum-Pacific tectonic realm, was tectonically controlled by the motion of the Pacific Plate and resembles the other petroliferous basins in eastern China. This basin developed early, and characterized intensive extension in the early stage and weak activity in the later stage of its development. Whereas the Qiongdongnan Basin was in a weak extension early and intensity of extension increased gradually. The relative geographical locations and the extensional histories of three basins ilustrate that the northern continental margin of South China Sea spread from south to north. On the other hand, the Qiongdongnan and Yinggehai Basins may have been controlled by the same tectonic regime since later Tertiary, whereas the Zhujiangkou Basin was not meaningfully influenced. So, the tectono-thermal evolution character of the Qiongdonnan basin should be closely to the other two. It may be concluded that the three basins have been developed within the active continental margin since Tertiary, and the local lithosphere might undergo intensive extension-perhaps two or three times of episodic extension occurred. Extension lead to large tectonoc subsidence and extreme thick Tertiary sediments for hydrocarbon generation in the basins. In response to the periodic extension of the basins, the palaeothermal flow were also periodical. The three basins all have the characteristics of multi-phase thermal evolutions that is good for oil-gas generation. And the overpressure expands the depth range of oil-gas habitat, which is meaningful to petroleum exploration in this region.

均质弹性条件下多洞平面问题解析算法及其应用研究

The scholars in the world have been trying to find an effective analytic algorithm of multiple hole problems usually meet in engineering designs. Though some studies on circular or elliptic holes had been achieved under specific conditions, no efforts were made to any multiple hole problems that is most significant for engineering designs. The author has made further studies on any multiple hole problems, using complex variable function method and Schwarz alternating method. After solving a series of technological difficulties, the author obtains an effective analytic algorithm, and acquires stress field and displacement field with high accuracy, which can be conducted for arbitrary many iterations according to practical accuracy requirements. In addition, th solution of stress and displacement fields, even for multiple holes of complex shapes and smaller distances. Further, the author made preliminary studies on viscoelastic displacement solution for any double holes. In terms of the obtained displacement solution of any multiple holes, this paper studies displacement back-analysis for the excavations of two tunnels, and find that the back-analysis method is accurate. Additionally, the author presents the mathematical prove of inversion uniqueness for ground stresses, elastic modulus and Poisson ratio. The author believes that the accurate analytic algorithm provided in this paper will presents an effective way to stress and displacement analysis for any multiple hole problems, optimal arrangement of multiple holes, hole shape optimization of multiple holes, etc..

地质规律约束下的波阻抗反演

Impedance inversion is very important in seismic technology. It is based on seismic profile. Good inversion result is derived from high quality seismic profile, which is formed using high resolution imaging resolution. High-resolution process demands that signal/noise ratio is high. It is very important for seismic inversion to improve signal/noise ratio. the main idea is that the physical parameter (wave impedance), which describes the stratigraphy directly, is achieved from seismic data expressing structural style indirectly. The solution of impedance inversion technology, which is based on convolution model, is arbitrary. It is a good way to apply the priori information as the restricted condition in inversion. An updated impedance inversion technology is presented which overcome the flaw of traditional model and highlight the influence of structure. Considering impedance inversion restricted by sedimentary model, layer filling style and congruence relation, the impedance model is built. So the impedance inversion restricted by geological rule could be realized. there are some innovations in this dissertation: 1. The best migration aperture is achieved from the included angle of time surface of diffracted wave and reflected wave. Restricted by structural model, the dip of time surface of reflected wave and diffracted wave is given. 2. The conventional method of FXY forcasting noise is updated, and the signal/noise ratio is improved. 3. Considering the characteristic of probability distribution of seismic data and geological events fully, an object function is constructed using the theory of Bayes estimation as the criterion. The mathematics is used here to describe the content of practice theory. 4. Considering the influence of structure, the seismic profile is interpreted to build the model of structure. A series of structure model is built. So as the impedance model. The high frequency of inversion is controlled by the geological rule. 5. Conjugate gradient method is selected to improve resolving process for it fit the demands of geophysics, and the efficiency of algorithm is enhanced. As the geological information is used fully, the result of impedance inversion is reasonable and complex reservoir could be forecasted further perfectly.

火成岩储层裂缝识别、预测理论与方法研究－以济阳坳陷火成岩油藏为例

Since 1970s, igneous reservoirs such as Shang741, Bin674 and Luol51 have been found in Jiyang depression, which are enrichment and heavy-producing. Showing good prospect of exploration and development, igneous reservoirs have been the main part of increasing reserves and production in Shengli oilfield. As fracture igneous reservoir being an extraordinary complex concealed reservoir and showing heavy heterogeneity in spatial distribution, the study of recognition, prediction, formation mechanism and the law of distribution of fracture is essential to develop the reservoir. Guided by multiple discipline theory such as sedimentology, geophysics, mineralogy, petroleum geology, structural geology and reservoir engineering, a set of theories and methods of recognition and prediction of fractured igneous rock reservoir are formed in this paper. Rock data, three-dimensional seismic data, log data, borehole log data, testing data and production data are combined in these methods by the means of computer. Based on the research of igneous rock petrography and reservoir formation mechanism, emphasized on the assessment and forecast of igneous rock reservoir, aimed at establishing a nonhomogeneity quantification model of fractured igneous rock reservoir, the creativity on the fracture recognition, prediction and formation mechanism are achieved. The research result is applied to Jiyang depression, suggestion of exploration and development for fractured igneous rock reservoir is supplied and some great achievement and favourable economic effect are achieved. The main achievements are gained as follows: 1. The main facies models of igneous rock reservoir in JiYang depression are summarized. Based on data and techniques of seism, well log and logging,started from the research of single well rock facies, proceeded by seismic and log facies research, from point to line and line to face, the regional igneous facies models are established. And hypabyssal intrusion allgovite facies model, explosion volcaniclastic rock facies model and overfall basaltic rocks facies model are the main facies models of igneous rock reservoir in JiYang depression. 2. Four nonhomogenous reservoir models of igneous reservoirs are established, which is the base of fracture prediction and recognition. According to characteristics of igneous petrology and spatial types of reservoir, igneous reservoirs of Jiyang depression are divided into four categories: fractured irruptive rock reservoir, fracture-pore thermocontact metamorphic rock and irruptive rock compound reservoir, pore volcanic debris cone reservoir and fracture-pore overfall basaltic rock reservoir. The spatial distribution of each model's reservoir has its features. And reservoirs can be divided into primary ones and secondary ones, whose mechanism of formation and laws of distribution are studied in this paper. 3. Eight geologic factors which dominate igneous reservoirs are presented. The eight geologic factors which dominates igneous reservoirs are igneous facies, epigenetic tectonics deformation, fracture motion, intensity of intrusive effect and adjoining-rock characters, thermo-contact metamorphic rock facies, specific volcano-tectonic position, magmatic cyclicity and epigenetic diagenetic evolution. The interaction of the eight factors forms the four types nonhomogenous reservoir models of igneous reservoirs in Jiyang depression. And igneous facies and fracture motion are the most important and primary factors. 4. Identification patterns of seismic, well log and logging facies of igneous rocks are established. Igneous rocks of Jiyang depression show typical reflecting features on seismic profile. Tabular reflection seismic facies, arc reflection seismic facies and hummocky or mushroom reflection seismic facies are the three main facies. Logging response features of basic basalt and diabase are shown as typical "three low and two high", which means low natural gamma value, low interval transit-time, low neutron porosity, high resistivity and high density. Volcaniclastic rocks show "two high and three low"-high neutron porosity, high interval transit-time, low density, low-resistance and low natural gamma value. Thermo-contact metamorphic rocks surrounding to diabase show "four high and two low" on log data, which is high natural gamma value, high self-potential anomaly, high neutron porosity, high interval transit-time and low density and low-resistance. Based on seismic, well log and logging data, spatial shape of Shang 741 igneous rock is described. 5. The methods of fracture prediction and recognition for fractured igneous reservoir are summarized. Adopting FMI image log and nuclear magnetic resonance log to quantitative analysis of fractured igneous reservoir and according to formation mechanism and shape of fracture, various fractures are recognized, such as high-angle fracture, low-angle fracture, vertical fracture, reticulated fracture, induced fracture, infilling fracture and corrosion vug. Shang 741 intrusive rock reservoir can be divided into pore-vug compound type, pore fracture type, micro-pore and micro-fracture type. Physical properties parameters of the reservoir are computed and single-well fracture model and reservoir parameters model are established. 6. Various comprehensive methods of fracture prediction and recognition for fractured igneous reservoir are put forward. Adopting three-element (igneous facies, fracture motion and rock bending) geologic comprehensive reservoir evaluation technique and deep-shallow unconventional laterolog constrained inversion technique, lateral prediction of fractured reservoir such as Shang 741 is taken and nonhomogeneity quantification models of reservoirs are established.

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

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.