印度与欧亚板块强烈相互作用区深部结构特征研究

The central-south Tibet is a part of the products of the continental plate collision between Eurasia and India. To study the deep structure of the study area is significant for understanding the dynamics of the continental-continental collision. A 3-D density model matched well with the observations in the central-south Tibet was proposed in this study. In addition, this study has also used numerical simulation method to prove that Quasi-Love (QL) wave is deduced by anisotropy variation but not by lateral heterogeneity. Meanwhile, anisotropy variation in the upper mantle of the Qiangtang terrane and Lhasa terrane is detected by the QL waves observed in recorded seismograms. Based on the gravity modeling, some results are summarized as follows: 1) Under the constrain of geometrical structure detected by seismic data, a 3-D density model and Moho interface are proposed by gravity inversion of the central-south Tibet. 2) The fact that the lower crustal densities are smaller than 3.2 g/cm3, suggests absence of eclogite or partial eclogitization due to delamination under the central-south Tibet. 3) Seismicity will be strong or weak in the most negative Bouguer gravity anomaly. So there is no a certain relationship between seismicity and Bouguer gravity anomaly. 4) Crustal composition are determined after temperature-pressure calibration of seismic P wave velocity. The composition of lower crust might be one or a mixture of: 1. amphibolite and greenschist facies basalt beneath the Qiangtang terrane; 2. gabbro-norite-troctolite and mafic granulite beneath the Lhasa terrane. Because the composition of the middle crust cannot be well constrained by the above data set, the data set published by Rudnick & Fountain (1995) is used for comparison. It indicated the composition of the middle crust is granulite facies and might be pelitic gneisses.Granulite facies used to be interpreted as residues of partial melting, which coincidences with the previous study on partial melting middle crust. Amphibolite facies are thought to be produced after delamination, when underplating works in the rebound of the lower crust and lithospheric mantle. From the seismology study, I have made several followed conclusions: 1) Through the numerical simulation experiment of surface wave propagating in heterogeneity media, we can find that amplitude and polarization of surface wave only change a little when considering heterogeneity. Furthermore, it is proved that QL waves, generated by surface wave scattering, are caused by lateral variation of anisotropy but not by heterogeneity. 2) QL waves are utilized to determine the variation of uppermost mantle anisotropy of the Tibetan plateau. QL waves are identified from the seismograms of the selected paths recorded by the CAD station. The location of azimuth anisotropy gradient is estimated from the group velocities of Rayleigh wave, Love wave and QL wave. It suggests that south-north lateral variation of azimuthal anisotropy locates in Tanggula mountain, and east-west lateral variation in the north of Gandese mountain with 85°E longitude and near the Jinsha river fault with 85°E longitude.

叠前地震多角度三参数同步反演方法

Pre-stack seismic inversion has become the emphasis and hotspot owing to the exploration & exploitation of oil field and the development of seismic technology. Pre-stack seismic inversion has the strongpoint of making the most of amplitude versus offset compared with the post-stack method. In this dissertation, the three parameters were discussed from multi-angle reflectance of P-wave data based on Zoeppritz’s and Aki & Richard’s equation, include P-wave velocity, S-wave velocity, and density. The three parameters are inversed synchronously from the pre-stack multi-angle P-wave data, based on rockphysics model and aimed at the least remnant difference between model simulation and practical data. In order to improve the stability of inversion and resolution to thin bed, several techniques were employed, such as the wavelet transform with multi-scale function, adding the Bayesian soft constraint and hard constraints (the horizon, structure and so on) to the inversion process. Being the result, the uncertainty of the resolution is reduced, the reliability and precision are improved, the significance of parameters becomes clearer. Meeting to the fundamental requirement of pre-stack inversion, some research in rockphysics are carried out which covered the simulation and inversion of S-wave velocity, the influence of pore fluids to geophysical parameters, and the slecting and analyzing of sensitive parameters. The difference between elastic wave equation modeling and Zoeppritz equation method is also compared. A series of key techniques of pre-stack seismic inversion and description were developed, such as attributes optimization, fluid factors, etc. All the techniques mentioned above are assembled to form a technique sets and process of synchronous pre-stack seismic inversion method of the three parameters based on rock physics and model simulation. The new method and technology were applied in many areas with various reservoirs, obtained both geological and economic significance, which proved to be valid and rational. This study will promote the pre-stack inversion technology and it’s application in hidden reservoirs exploration, face good prospects for development and application.

地震波传播模拟的边界元―体积元积分方程法及其近似求解

To improve the efficiency of boundary-volume integral equation technique, this paper is involved in the approximate solutions of boundary-volume integral equation technique. Firstly, based on different interpretations of the self-interaction and extrapolation operators of the resulting boundary integral equation matrix, two different hybrid BEM+Born series modeling schemes are formulated and validated through comparisons with the full-waveform BE numerical solutions for wave propagation simulation in a semicircular alluvial valley and a complex fault model respectively. Numerical experiments indicate that both the BEM+Born series modeling schemes are suitable for complex geological structures and significantly improve computational efficiency especially for the cases of high frequencies and multisource seismic survey. Then boundary-volume integral equation technique is illuminated in detail and verified by modeling wave propagation in complex media. Furthermore, the first-order and second-order Born approximate solutions for the volume-scattering waves are studied and quantified by numerical simulation in different random medium models. Finally, preconditioning generalized minimal residual method is applied to solve boundary-volume integral equation and compared with Gaussian elimination method. Numerical experiments indicate this method makes the calculations more efficient.

准噶尔盆地南缘台13井区侏罗系油藏表征与建模

Junggar Basin has a large amount of recoverable reserves, However, due to the unfavorable factors, such as bad seismic data quality, complex structure with many faults and less wells, the exploration of oil and gas is still relatively limited, so advanced theoretical guidance and effective technical supports are desirable. Based on the theories of sedimentology, as well as comprehensive studies of outcrops, seismic data, drilling data and setting of this area, the paper establishes the isochronous correlation framework, and analyzes the sedimentary facies types and provenance direction, and obtains the profile and plain maps of the sedimentary facies combined with the logging constrained inversion. Then the paper analyzes the reservoir controlling factors, reservoir lithology attribute, 4-property relationship and sensibility based on the sedimentary facies research, and sets up a 3D geological model using facies controlled modeling. Finally, the paper optimizes some target areas with the conclusions of reservoir, structure and reservoir formation.Firstly, the paper establishs the isochronous correlation framework by the seismic data, drilling data and setting of this area. The sedimentary facies in Tai13 well block are braided river and meandering river according to the analysis of the lithology attribute, logging facies and sedimentary structure attribute of outcrop. The concept of “wetland” is put forward for the first time. The provenance direction of Badaowan and Qigu formation is obtained by the geology setting, sedimentary setting and paleocurrent direction. The paper obtains the profile and plain maps of the sedimentary facies from the sand value of the wells and the sand thickness maps from the logging constrained inversion. Then, this paper takes characteristics and control factors of the Jurassic reservoirs analysis on thin section observation, scanning transmission electron microscope observation and find out the petrology characteristics of reservoir, space types of reservoir and lithofacies division. In this area, primary pores dominate in the reservoir pores, which believed that sedimentation played the most important roles of the reservoir quality and diagenesis is the minor factor influencing secondary porosity. Using stochastic modeling technique，the paper builds quantitative 3-D reservoir Parameter. Finally, combined the study of structure and reservoir formation, the reservoir distribution regularity is concluded: (a) structures control the reservoir formation and accumulation. (b) Locating in the favorable sedimentary facies belt. And the area which meets these conditions mentioned above is a good destination for exploration.

地震波阻抗反演的迭代正则化算法

Post-stack seismic impedance inversion is the key technology of reservoir prediction and identification. Geophysicists have done a lot of research for the problem, but the developed methods still cannot satisfy practical requirements completely. The results of different inversion methods are different and the results of one method used by different people are different too. The reasons are due to the quality of seismic data, inaccurate wavelet extraction, errors between normal incidence assumption and real situation, and so on. In addition, there are two main influence factors: one is the band-limited property of seismic data; the other is the ill-posed property of impedance inversion. Thus far, the most effective way to solve the band-limited problem is the constrained inversion. And the most effective way to solve ill-posed problems is the regularization method assisted with proper optimization techniques. This thesis systematically introduces the iterative regularization methods and numerical optimization methods for impedance inversion. A regularized restarted conjugate gradient method for solving ill-posed problems in impedance inversion is proposed. Theoretic simulations are made and field data applications are performed. It reveals that the proposed algorithm possesses the superiority to conventional conjugate gradient method. Finally, non-smooth optimization is proposed as the further research direction in seismic impedance inversion according to practical situation.

波动方程基准面校正方法研究

Datuming which has not been well solved in complex areas is a long-existing problem in seismic processing and imaging. Theoretically, Wave-equation datuming(WED) works well in the areas with substantial surface topography and areas of complex velocity structure. However, many difficulties still exist in practice. There are three main reasons: (1) It’s difficult to obtain the velocity model. (2) The computational cost is high and the efficiency is low. (3) Reflection waveform distortions are introduced by low S/N ratio in seismic data. The second and third problems are involved in the paper. To improve computational efficiency, DP1 proposed by Fu Li-Yun is applied in WED. Some quantitative and semi-quantitative conclusions of assessing the computational accuracy and efficiency have been obtained by comparing the adaptation of three operators( PS, SSF, DP1) to the surface topography and the lateral velocity variation. Moreover, the impacts of near surface scattering associated with complex surface topography on WED is analyzed theoretically. According to the analysis results, the following conclusions have been obtained. WED is stable and effective when the field data has high S/N ratio and velocity model is accurate. However ,it doesn’t work well when S/N ratio of field data is low. So denoising techniques in process of WED is important for low S/N data. The paper presents the theoretical analysis for the issues facing WED, which is expected to provide a useful reference to the further development of this technology.

井约束外推反演可信度与井震相关性关系研究

With the deepening development of oil-gas exploration and the sharp rise in costs, modern seismic techniques had been progressed rapidly. The Seismic Inversion Technique extracts seismic attribute from the seismic reflection data, inverses the underground distribution of wave impedance or speed, estimates reservoir parameters, makes some reservoir prediction and oil reservoir description as a key technology of Seismic exploration, which provides a reliable basic material for oil-gas exploration. Well-driven SI is essentially an seismic-logging joint inversion. The low, high-frequency information comes from the logging information, while the structural characteristics and medium frequency band depend on the seismic data. Inversion results mainly depend on the quality of raw data, the rationality of the process, the relativity of synthetic and seismic data, etc. This paper mainly research on how the log-to-seismic correlation have affected the well-driven seismic inversion precision. Synthetic, the comparison between middle –frequency borehole impedance and relative seismic impedance and well-attribute crossplots have been taken into account the log-to-seismic correlation. The results verify that the better log-to-seismic correlation, the more reliable the seismic inversion result, through the analysis of three real working area (Qikou Sag, Qiongdongnan basin, Sulige gas field).

人工源宽角地震资料中的广角效应及其对地壳速度成像的影响

The seismic wide-angle reflection/refraction method is the one of the most effective method for probing the crustal and upper mantle structure. It mainly uses the wide-angle reflection information from the boundary in the crust and the top boundary of the upper mantle to rebuild the crust and upper mantle structure. Through analyzing the reflection and transmission coefficients of various incident waves on the interface, we think relative to the pre-critical angle reflection information the post critical angle reflection information that received by wide-angle seismic data exists a time-shift effect with the offset variation, and then it must cause the error for velocity analysis and structure image. The feature of the wide-angle seismic wave field of the fourteen representative crust columns tell us that the wide-angle effects in the different representative tectonic units for the interface depth and the interval velocity in crust. We studied the features of the wide-angle seismic wave field through building the crust model and inverse its travel time by GA method to know the wide-angle influence on crustal velocity image. At last we finished the data processing of the Tunxi-Wenzhou wide-angle seismic profile. The results are as following: (1) Through building crust model, we labeled the travel time for all the phases by ray tracing method and remove wide-angle effects method, it revealed the wide-angle effect exists in the seismic data. (2) The travel time inversion by GA method can tell us that the depth by traditional ray tracing method is shallower than the result by remove wide-angle effects method, the latter can recover the crust structure model in effect. (3) We applied the two method mentioned before to the fourteen representative crust columns in China. It indicates that the removed wide-angle effect method in travel time inversion is reasonable and effective. (4) The real data processing from Tunxi-Wenzhou wide-angle seismic profile give us the basic structure through the two ways. The main influence exhibits in the difference of the interval velocity of the curst, and the wide-angle effects in shallow interface are stronger than the deep interface.

中国大陆地区远震体波接收函数与岩石圈厚度的研究

Movements of separation and convergence between the continental plates, as well as the interaction beween the lithosphere and asthenosphere is the dominant factor in plate evolution. Moreover, those phenomena, the formation, enrichment and storage of energy and mineral strorage, as well as intraplate earthquakes are all relate to plate movement and evolution. Therefore, the study of continental lithosphere, is not only helpful to analysing the dynamic model between lithosphere and asthenoshere as well as different plates, but also important to the nation's economy and the people's livelihood. And the lithospheric thickness or Lithosphere-Asthenosphere Boundary (LAB) is one of the most important parameters in study of continental lithospheric formation and evolution. Chinese continent composed by many small plates, possesses diverse type of lithospheric structure. But our knowledge ahout Chinese continental lithosphere, especially the regional research, is almost based on the low-resolution results of surface wave dispersion analysis and seismic wave tomography. Howere, recently a technique employing S-to-P converted body waves (the S receiver function technique) has been developed that can be used to identify the LAB with a higher resolution. This thesis has collected waveform data of 232 broadband seismic stations that are located in China and neighboring regions. Using teleseismic S-wave and P-wave receiver functions have studied the Chinese continental lithospheric structure. The results of this study indicate that, the thickness of Chinese continental lithosphere become thinner from west to east, and obvious difference exists between different blocks. Four types of lithosphere have been detected: (1) Convergence thicking lithosphere in Tibetan Plateau; (2) Stable lithosphere in Tarim basin and upper Yangtze craton; (4) Active lithosphere in Orogenic belts; (4) Break-up thinning lithosphere in east China craton

南海张裂大陆边缘数值模拟研究

The South China Sea (SCS) is one of the largest marginal seas in the western Pacific, which is located at the junction of Eurasian plate, Pacific plate and Indian-Australian plate. It was formed by continent breakup and sea-floor spreading in Cenozoic. The complicated interaction among the three major plates made tectonic movement complex and geological phenomena very rich in this area. The SCS is an ideal place to study the formation and evolution of rifted continental margin and sea-floor spreading since it is old enough to have experienced the major stages of the basin evolution but still young enough to have preserved its original nature. As the demand for energy grows day by day in our country, the deep water region of the northern continental margin in the SCS has become a focus of oil and gas exploration because of its huge hydrocarbon potential. Therefore, to study the rifted continental margin of the SCS not only can improve our understanding of the formation and evolution processes of rifted continental margin, but also can provide theoretical support for hydrocarbon exploration in rifted continental margin. This dissertation mainly includes five topics as follows: (1) Various classic lithosphere stretching models are reviewed, and the continuous non-uniform stretching model is modified to make it suitable for the case where the extension of lithopheric mantle exceeds that of the crust. Then simple/pure shear flexural cantilever model is applied to model the basement geometries of SO49-18 profile in the northern continental margin of the SCS. By fitting the basements obtained by using 2DMove software with modeling results, it is found that the reasonable effective elastic thickness is less than 5km in this region. According to this result, it is assumed that there is weak lower crust in the northern continental margin in the SCS. (2) We research on the methods for stretching factor estimation based on various lithosphere stretching models, and apply the method based on multiple finite rifting model to estimate the stretching factors of several wells and profiles in the northern continental margin of the SCS. (3) We improve one-dimension strain rate inversion method with conjugate gradient method, and apply it to invert the strain rate of several wells in the northern continental margin of the SCS. Two-dimension strain rate forward modeling is carried out, and the modeling results show that effective elastic thickness is a key parameter to control basin’s geometry. (4) We simulate divergent upwelling mantle flow model using finite difference method, and apply this newly developed model to examine the formation mechanism of the northwest and central sub-basin in the SCS. (5) We inverse plate thickness and basal temperature of oceanic lithosphere using sea-floor ages and bathymetries of the North Pacific and the North Atlantic based on varied-parameters plate model, in which the heat conductivity, heat capacity and coefficient of thermal expansion depend on temperature or depth. A new empirical formula is put forward based the inversed parameters, which depicts the relation among sea-floor age, bathymetry and heat flow. Then various similar empirical formulae, including the newly developed one, are applied to examine the sea-floor spread issue in the SCS based on the heat flow and bathymetry data of the abyssal sub-basin.

三维非均匀各向异性介质中的纵波走时反演研究

Many observations show that seismic anisotropy is very common in the crust and upper mantle of the Earth. Seismic anisotropy can provide some clue about the changing and transporting process inside the earth. in recent years, abundant earthquake travel time data are accumulated, computers become more powerful, and these make the inversion of earthquake travel time data practical. In this thesis we studied the theory of elastic wave in anisotropic media, some formule for travel time inversion were derived. We present an iterative procedure to determine 21 elastic parameters from qP wave travel times. No a priori assumptions about heterogeneity and anisotropy of the model are made. The procedure is suitable for the case when we know nothing about the symmetry of anisotropy of the media, as well as for the case of earthquake travel time inversion which may contain various symmetry of anisotropy. The procedure is tested with a synthetic multiple-source offset VSP experiment. The results proved that the formulae are correct, and the procedure is practical. The results and the related theory indicate that the anisotropic inversion needs more rays than isotropic case. For a 2-D weak anisotropic (WA) medium, we need at least 5 rays in different directions to retrieve the elastic parameters on one grid point, and for a 3-D WA medium we need at least 15 rays in different directions to retrieve the elastic parameters on one grid point. The results also indicate that the starting background velocity has no influence on the final results, at least for the model we specified. Our results also show that insufficient illumination coverage will slow down the convergence rate, and make the results more sensitive to noise. We apply the procedure to a set of field travel time data. The data is from an artificial seismic observation. This observation is for locating micro-seismic events around a tunnel, its purpose is to find out if the digging process and the stress condition around the tunnel can generate micro-cracks. The size of this area is around 100m. The anisotropy derived from qP travel times is the same as the anisotropy showed by apparent velocities, and is also consistent with the anisotropy derived from S-wave splitting phenomena.

各向同性与各向异性介质中地震波旅行时正演模拟与反演成像

By seismic tomography, interesting results have been achieved not only in the research of the geosphere with a large scale but also in the exploration of resources and projects with a small scale since 80'. Compared with traditional inversion methods, seismic tomography can offer more and detailed information about subsurface and has been being paid attention by more and more geophysicists. Since inversion based on forward modeling, we have studied and improved the methods to calculate seismic traveltimes and raypaths in isotropic and anisotropic media, and applied the improved forward methods to traveltime tomography. There are three main kinds of methods to calculate seismic traveltime field and its ray path distribution, which are ray-tracing theory, eikonal equation by the finite-difference and minimum traveltime tree algorithm. In ray tracing, five methods are introduced in the paper, including analytic ray tracing, ray shooting, ray bending, grid ray tracing and rectangle grid ray perturbation with three points. Finite-difference solution of eikonal equation is very efficient in calculation of seismic first-break, but is awkward in calculation of reflection traveltimes. We have put forward a idea to calculate traveltimes of reflected waves using a combining way of eikonal equation method and other one in order to improve its capability of dealing with reflection waves. The minimum traveltime tree algorithm has been studied with emphases. Three improved algorithms are put forward on the basis of basic algorithm of the minimum traveltime tree. The first improved algorithm is called raypath tracing backward minimum traveltime algorithm, in which not only wavelets from the current source but also wavelets from upper source points are all calculated. The algorithm can obviously improve the speed of calculating traveltimes and raypaths in layered or blocked homogeneous media and keep good accuracy. The second improved algorithm is raypath key point minimum traveltime algorithm in which traveltimes and raypaths are calculated with a view of key points of raypaths (key points of raypths mean the pivotal points which determine raypaths). The raypath key point method is developed on the basis of the first improved algorithm, and has better applicability. For example, it is very efficient even for inhomogeneous media. Another improved algorithm, double grid minimum traveltime tree algorithm, bases upon raypath key point scheme, in which a model is divided with two kinds of grids so that the unnecessary calculation can be left out. Violent undulation of curved interface often results in the phenomenon that there are no reflection points on some parts of interfaces where there should be. One efficacious scheme that curved interfaces are divided into segments, and these segments are treated respectively is presented to solve the problem. In addition, the approximation to interfaces with discrete grids leads to large errors in calculation of traveltimes and raypaths. Noting the point, we have thought a new method to remove the negative effect of mesh and to improve calculation accuracy by correcting the traveltimes with a little of additional calculation, and obtained better results.

小波域地震资料的处理方法（WSSP）研究

The modeling formula based on seismic wavelet can well simulate zero - phase wavelet and hybrid-phase wavelet, and approximate maximal - phase and minimal - phase wavelet in a certain sense. The modeling wavelet can be used as wavelet function after suitable modification item added to meet some conditions. On the basis of the modified Morlet wavelet, the derivative wavelet function has been derived. As a basic wavelet, it can be sued for high resolution frequency - division processing and instantaneous feature extraction, in acoordance with the signal expanding characters in time and scale domains by each wavelet structured. Finally, an application example proves the effectiveness and reasonability of the method. Based on the analysis of SVD (Singular Value Decomposition) filter, by taking wavelet as basic wavelet and combining SVD filter and wavelet transform, a new de - noising method, which is Based on multi - dimension and multi-space de - noising method, is proposed. The implementation of this method is discussed the detail. Theoretical analysis and modeling show that the method has strong capacity of de - noising and keeping attributes of effective wave. It is a good tool for de - noising when the S/N ratio is poor. To give prominence to high frequency information of reflection event of important layer and to take account of other frequency information under processing seismic data, it is difficult for deconvolution filter to realize this goal. A filter from Fourier Transform has some problems for realizing the goal. In this paper, a new method is put forward, that is a method of processing seismic data in frequency division from wavelet transform and reconstruction. In ordinary seismic processing methods for resolution improvement, deconvolution operator has poor part characteristics, thus influencing the operator frequency. In wavelet transform, wavelet function has very good part characteristics. Frequency - division data processing in wavelet transform also brings quite good high resolution data, but it needs more time than deconvolution method does. On the basis of frequency - division processing method in wavelet domain, a new technique is put forward, which involves 1) designing filter operators equivalent to deconvolution operator in time and frequency domains in wavelet transform, 2) obtaining derivative wavelet function that is suitable to high - resolution seismic data processing, and 3) processing high resolution seismic data by deconvolution method in time domain. In the method of producing some instantaneous characteristic signals by using Hilbert transform, Hilbert transform is very sensitive to high - frequency random noise. As a result, even though there exist weak high - frequency noises in seismic signals, the obtained instantaneous characteristics of seismic signals may be still submerged by the noises. One method for having instantaneous characteristics of seismic signals in wavelet domain is put forward, which obtains directly the instantaneous characteristics of seismic signals by taking the characteristics of both the real part (real signals, namely seismic signals) and the imaginary part (the Hilbert transfom of real signals) of wavelet transform. The method has the functions of frequency division and noise removal. What is more, the weak wave whose frequency is lower than that of high - frequency random noise is retained in the obtained instantaneous characteristics of seismic signals, and the weak wave may be seen in instantaneous characteristic sections (such as instantaneous frequency, instantaneous phase and instantaneous amplitude). Impedance inversion is one of tools in the description of oil reservoir. one of methods in impedance inversion is Generalized Linear Inversion. This method has higher precision of inversion. But, this method is sensitive to noise of seismic data, so that error results are got. The description of oil reservoir in researching important geological layer, in order to give prominence to geological characteristics of the important layer, not only high frequency impedance to research thin sand layer, but other frequency impedance are needed. It is difficult for some impedance inversion method to realize the goal. Wavelet transform is very good in denoising and processing in frequency division. Therefore, in the paper, a method of impedance inversion is put forward based on wavelet transform, that is impedance inversion in frequency division from wavelet transform and reconstruction. in this paper, based on wavelet transform, methods of time - frequency analysis is given. Fanally, methods above are in application on real oil field - Sansan oil field.

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

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.