地球介质非线性与各向异性的地震波响应及其成像效应

The practice of geophysical prospecting shows us the complex interior earth. The studies of the complexity play an important role and practical guide for the subsurface structure. At present, the complexity of the earth mainly means lateral and vertical homogeneity, anisotropy and non-linear quality. And the anisotropy and non-linear media studies become the frontier in seismology and exploration seismology. This paper summarizes the development of complexities and presents the forward and inverse in the non-linear and anisotropic media. Firstly, the paper introduces the theory of seismic wave propagation in the non-linear and anisotropic media, the theoretical basis for simulation and inversion research. Secondly, high quality numerical simulation method with little dispersion has been developed to investigate the influence of complexity including anisotropy and non-linear multi-component seismograms. Because most real data in seismology have a single component, we developed two aspects work on anisotropic multi-component imaging. One is prestack reflection migration. The result show that distorted images are obtained if data from anisotropic media are migrated using isotropic extrapolation. Moreover, image quality will be improved greatly after considering anisotropy in subsurface layers. The other one is the we take advantage of multi-component data to inversion of the anisotropic parameters jointly seimic reflection travel time and polarization information. Based on these research works, we get the following results: 1.Combing numerical simulation, systematical studies indicate that anisotropy and non-linear seismograms characters are significant to detect cracked belts in the earth and to understand deformation field and mechanism. 2.Based on anisotropic media models, we developed an efficient prestack migration method for subsurface structure and different observation methods seismic data, which improving the imaging quality with VSP, seismograms and real data. 3.Jointly seismic inversion combining seismic anisotropic reflection traveltimes and polarizations data show that the complete wrong inversion and the following explanation will be resulted by ignoring anisotropy.

东海陆架盆地复杂天然气藏地震勘探判别技术研究

Exploration study proves that East sea shelf basin embeds abundant hydrocarbon resources. However, the exploration knowledge of this area is very low. Many problems in exploration are encountered here. One of them is that the gas reservoir of this area, with rapid lateral variation, is deeply buried. Correlation of Impendence between sandstone, gas sand and shale is very poor. Another problem is that the S/N ratio of the seismic data is very low and multiples are relatively productive which seriously affect reservoir identification. Resolution of the seismic data reflected from 2500-3000 meter is rather low, which seriously affects the application of hydrocarbon direct identification (HDI) technology. This research established a fine geological & geophysical model based on drilling、well logging、geology&seismic data of East sea Lishui area. A Q value extraction method from seismic data is proposed. With this method, Q value inversion from VSP data and seismic data is performed to determine the subsurface absorption of this area. Then wave propagation and absorption rule are in control. Field acquisition design can be directed. And at the same time, with the optimization of source system, the performance of high resolution seismic acquisition layout system is enhanced. So the firm foundation is ensured for east sea gas reservoir exploration. For solving the multiple and amplitude preserving problems during the seismic data processing, wave equation pre-stack amplitude preservation migration and wave equation feedback iteratively multiple attenuation technologies are developed. Amplitude preservation migration technology can preserve the amplitude of imaging condition and wave-field extrapolation. Multiple removing technology is independent of seismic source wavelet and velocity model, which avoiding the weakness of Delft method. Aiming at the complicated formation condition of the gas reservoir in this area, with dissecting typical hydrocarbon reservoir, a series of pertinent advanced gas reservoir seismic identification technologies such as petrophysical properties analyzing and seismic modeling technology、pre-stack/post-stack joint elastic inversion, attribute extraction technology based on seismic non-stationary signal theory and formation absorption characteristic and so on are studied and developed. Integrated analysis of pre-stack/post-stack seismic data, reservoir information, rock physics and attribute information is performed. And finally, a suit of gas reservoir identification technology is built, according to the geological and geophysical characteristics of this area. With developed innovative technologies, practical application and intergrated interpretation appraisal researches are carried out in Lishui 36-1.The validity of these technologies is tested and verified. Also the hydrocarbon charging possibility and position of those three east sea gas exploration targets are clearly pointed out.

油气输导体系研究及应用─以准噶尔盆地东部白家海凸起侏罗系为例

Baijiahai uplift is an important hydrocarbon accumulation belt in eastern Jungger Basin, on which Cainan oilfield and lithologic hydrocarbon reservoir named Cai 43 have been discovered and both of them share the same target formation of Jurassic. However, in the subsequent exploration at this region, several wells that designed for lithologic traps of Jurassic were eventually failed, and that indicates the controlling factors of lithologic reservoir distribution are far more complicated than our previous expectation. This dissertation set the strata of the Jurassic in well Cai 43 region as the target, and based on the integrated analysis of structure evolution、fault sealing ability、simulations of sedimentary microfacies and reservoir beds、distribution analysis of high porosity-high permeability carrier beds、drive forces of hydrocarbons、preferential conduit system and conduit model as well as critical values of the reservoir physical properties for hydrocarbon charging, a special method that different from the conventional way to predict favorable lithologic traps was established. And with this method the controlling factors of the hydrocarbon reservoirs formation are figured out, and further more, the favorable exploration targets are point out. At Baijiahai uplift, fault plays as a crucial factor in the process of the hydrocarbon reservoir formation. In this study, it is found out that the availability of a fault that work as the seal for oil and gas are different. The critical value of the lateral mudstone smear factor (Kssf), which is used to measure the lateral sealing ability of fault, for oil is 3.9 while that for gas is 2.1; and the critical value of vertical sealing factor (F), which similarly a measurement for the vertical sealing ability of fault, for oil is 7.3 while that for gas is 5.1. Dongdaohaizi fault belt that possessed well lateral sealing ability since later Cretaceous have bad vertical sealing ability in later Cretaceous, however, it turns to be well now. Based on the comparison of the physical properties that respectively obtained from electronic log calculating、conventional laboratory rock analysis and the additive-pressure bearing laboratory rock analysis, we established the functions through which the porosity and permeability obtained though conventional method can be converted to the values of the subsurface conditions. With this method, the porosity and permeability of the Jurassic strata at the time of previous Tertiary and that in nowadays are reconstructed respectively, and then the characteristics of the distribution of high porosity-high permeability carrier beds in the evolution processes are determined. With the result of these works, it is found that both well Cai 43 region and Cainan oilfield are located on the preferential conduit direction of hydrocarbon migration. This conclusion is consistent with the result of the fluid potential analysis, in which fluid potential of nowadays and that of later Cretaceous are considered. At the same times, experiment of hydrocarbon injection into the addictive-pressure bearing rock is designed and conducted, from which it is found that, for mid-permeability cores of Jurassic, 0.03MPa is the threshold values for the hydrocarbon charging. And here, the conception of lateral pressure gradient is proposed to describe the lateral driving force for hydrocarbon migration. With this conception, it is found that hydrocarbons largely distributed in the areas where lateral pressure gradient is greater than 0. 03MPa/100m. Analysis of critical physical properties indicated that the value of the critical porosity and critical permeability varied with burial depth, and it is the throat radius of a certain reservoir bed that works as a key factor in controlling hydrocarbon content. Three parameters are proposed to describe the critical physical properties in this dissertation, which composite of effective oil-bearing porosity、effective oil-bearing permeability and preferential flow coefficient. And found that critical physical properties, at least to some extent, control the hydrocarbon distribution of Jurassic in Baijiahai uplift. Synthesize the content discussed above, this dissertation analyzed the key factors i.e., critical physical properties、driving force、conduit system and fluid potential, which controlled the formation of the lithologic reservoir in Baijiahai uplift. In all of which conduit system and fluid potential determined the direction of hydrocarbon migration, and substantially they are critical physical properties of reservoir bed and the lateral pressure gradient that controlled the eventually hydrocarbon distribution. At the same times, sand bodies in the major target formation that are recognized by reservoir bed simulation are appraised, then predict favorite direction of the next step exploration of lithologic reservoir.

徐家围子断陷连片叠前时间偏移处理方法研究

In this paper, the detailed analysis of fundamental seismic data and theoretical method are given, and the tests of some new technologies are performed. For seismic data processing assembly, some key technologies are developed and applied, such as global static correction, amplitude consistency processing, wavelet consistency shaping, fine velocity model establishing and prestack time migration. These technologies can efficiently settle the problems during the course of multiple- block– jointed prestack time migration processing, and it is highly significant for holding the oil output of 40,000,000 tons for Daqing oilfield. Through the research of this dissertation, the following important contributions are shown: (1) The combination of near-surface model method and refraction static correction method is developed, and is applied to solve global static correction for the whole merging area. (2) Prestack amplitude normalization processing method based on fold is developed. The method eliminates the effects of fold on amplitude uniformity, and solves the problem of energy uniformity for tie-area prestack migration processing. (3) Wavelet consistency is investigated. For multiple survey blocks existing in the area, the optimum method of wavelet shaping is developed, which removes the waveform variance between two adjacent blocks. (4) Controlled velocity inversion (CVI) technique is used to establish migration velocity field. It can largely shorten the period of velocity modeling, and improve velocity analysis precision. (5) Float datum level technique is employed, and is able to guarantee prestack migration results of subsurface shallow layers. (6) The static partition of seismic data volume relating to migration aperture is firstly developed. And the precious imaging for huge data volume by prestack time migration is realized. (7) The numerical forward simulation and prestack migration processing is primarily combined to discuss the migration technique for a complex geology structure from practical field information. The combination of numerical simulation and prestack migration is a feasible way to solve the fine imaging of complex volcanic structure. And the combination approach can help to select appropriate migration parameters.

叠前时间偏移: 基于波动方程的架构、发展和应用

This dissertation starts from the point that the prestack time migration can been considered as an approximation of the prestack depth migration, giving a wave equation based prestack time migration approach. The new approach includes: analytically getting the travel time and amplitude based on the one way wave equation and the stationary-phase theory, using ‘spread’ imaging method and imaging following the prestack depth migration, updating the velocity model with respect to the flats of the events in CRP gathers. Based on this approach, we present a scheme that can image land seismic data without field static correction. We may determine the correct near surface velocities and stack velocities by picking up the residual correction of the events in the CRP gathers. We may get the rational migration section based on the updated velocities and correct the migration section from a floating datum plane to a universal datum plane. We may adaptively determine the migration aperture according to the dips of the imaging structures. This not only speed up the processing, but may suppress the migration noise produce by the extra aperture. We adopt the deconvolution imaging condition of wave equation migration. It may partially compensate the geometric divergence. In this scheme, we use the table-driven technique which may enhance the computational efficiency. If the subsurface is much more complicated, it may be impossible to distinguish the DTS curve. To solve this problem, we proposed a technique to determine the appropriate range of the DTS curve. We synthesize DTS panel in this range using different velocities and depths, and stack the amplitude around the zero time. Determine the correct velocity and location of the considered grid point by comparing the values.

基于量子蒙特卡罗的地震波反演方法

The primary approaches for people to understand the inner properties of the earth and the distribution of the mineral resources are mainly coming from surface geology survey and geophysical/geochemical data inversion and interpretation. The purpose of seismic inversion is to extract information of the subsurface stratum geometrical structures and the distribution of material properties from seismic wave which is used for resource prospecting, exploitation and the study for inner structure of the earth and its dynamic process. Although the study of seismic parameter inversion has achieved a lot since 1950s, some problems are still persisting when applying in real data due to their nonlinearity and ill-posedness. Most inversion methods we use to invert geophysical parameters are based on iterative inversion which depends largely on the initial model and constraint conditions. It would be difficult to obtain a believable result when taking into consideration different factors such as environmental and equipment noise that exist in seismic wave excitation, propagation and acquisition. The seismic inversion based on real data is a typical nonlinear problem, which means most of their objective functions are multi-minimum. It makes them formidable to be solved using commonly used methods such as general-linearization and quasi-linearization inversion because of local convergence. Global nonlinear search methods which do not rely heavily on the initial model seem more promising, but the amount of computation required for real data process is unacceptable. In order to solve those problems mentioned above, this paper addresses a kind of global nonlinear inversion method which brings Quantum Monte Carlo (QMC) method into geophysical inverse problems. QMC has been used as an effective numerical method to study quantum many-body system which is often governed by Schrödinger equation. This method can be categorized into zero temperature method and finite temperature method. This paper is subdivided into four parts. In the first one, we briefly review the theory of QMC method and find out the connections with geophysical nonlinear inversion, and then give the flow chart of the algorithm. In the second part, we apply four QMC inverse methods in 1D wave equation impedance inversion and generally compare their results with convergence rate and accuracy. The feasibility, stability, and anti-noise capacity of the algorithms are also discussed within this chapter. Numerical results demonstrate that it is possible to solve geophysical nonlinear inversion and other nonlinear optimization problems by means of QMC method. They are also showing that Green’s function Monte Carlo (GFMC) and diffusion Monte Carlo (DMC) are more applicable than Path Integral Monte Carlo (PIMC) and Variational Monte Carlo (VMC) in real data. The third part provides the parallel version of serial QMC algorithms which are applied in a 2D acoustic velocity inversion and real seismic data processing and further discusses these algorithms’ globality and anti-noise capacity. The inverted results show the robustness of these algorithms which make them feasible to be used in 2D inversion and real data processing. The parallel inversion algorithms in this chapter are also applicable in other optimization. Finally, some useful conclusions are obtained in the last section. The analysis and comparison of the results indicate that it is successful to bring QMC into geophysical inversion. QMC is a kind of nonlinear inversion method which guarantees stability, efficiency and anti-noise. The most appealing property is that it does not rely heavily on the initial model and can be suited to nonlinear and multi-minimum geophysical inverse problems. This method can also be used in other filed regarding nonlinear optimization.

松辽盆地南部腰英台区块低渗透储层整体压裂改造技术研究

The Yaoyingtai Block is located within the northeastern Changling Depression of southern Songliao Basin, where the reservoir sandstones are petrophysically characterized by very low permeability, which results in the low success probability of artificial fracturing, and the low oil yield by water injection in the course of oil production. In order to improve the situations as stated above, this research aims to work out an integral fracturing technology and strategy applicable to the low permeable reservoirs in Yaoyingtai Block. Under the guidance of geological theory, reservoir engineering and technology, the subsurface occurrences of natural and hydraulic fractures in the reservoirs are expected to be delineated, and appropriate fracturing fluids and proppants are to be optimized, based on the data of drilling, well logging, laboratory and field experiments, and geological data. These approaches lay the basis of the integral fracturing technology suitable for the low permeable reservoir in the study area. Based on core sample test, in-situ stress analysis of well logging, and forward and inversion stress field modeling, as well as fluid dynamic analysis, the maximum in-situ stress field is unraveled to be extended nearly along the E-W direction (clustering along N85-135°E) as is demonstrated by the E-W trending tensional fractures. Hydraulic fractures are distributed approximately along the E-W direction as well. Faulting activities could have exerted obvious influences on the distribution of fractures, which were preferentially developed along fault zones. Based on reservoir sensitivity analysis, integrated with studies on rock mechanics, in-situ stress, natural fracture distribution and production in injection-production pilot area, the influences of primary fractures on fracturing operation are analyzed, and a diagnostic technology for primary fractures during depressurization is accordingly developed. An appropriate fracturing fluid (hydroxypropyl guar gum) and a proppant (Yixing ceramsite, with a moderate-density, 0.45-0.9mm in size) applicable to Qingshankou Formation reservoir are worked out through extensive optimization analysis. The fracturing fluid can decrease the damage to the oil reservoir, and the friction in fracturing operation, improving the effect of fracturing operation. Some problems, such as sand-out at early stage and low success rate of fracturing operations, have been effectively solved, through pre-fracturing formation evaluation, “suspension plug” fracturing, real-time monitoring and limited-flow fracturing. Through analysis of fracture-bearing tight reservoir with variable densities and dynamic analysis of influences of well patterns on fracturing by using numerical simulation, a fracturing operation scheme for the Qingshankou Formation reservoir is proposed here as being better to compress the short factures, rather than to compress the long fractures during hydraulic fracturing. It is suggested to adopt the 450m×150m inverted 9-spot well pattern in a diamond shape with wells placed parallel to fractures and a half fracture length of 60-75m.

基于盆地地下水模型与山区地表水模型耦合的济源盆地地下水模拟研究

The aim of this study is to increase the precision of groundwater modeling. The way is use the distributed model calculate the mountain basin groundwater lateral discharge and the river runoff. With appropriate technique help, the groundwater model can couple the distributed model results. This paper’s study object is makeing the distributed hydrological model HEC-HMS coupled to the popular groundwater model Visual MODFLOW. The application example is Jiyuan basin which is a typical basin of North China. HEC-HMS can calculate the surface runoff and subsurface runoff at mountain-pass. The subsurface runoff can turn to recharge well straightly. The water level - runoff course and Trial method is used to back analyze the parameters of surface runoff to Visual MODFLOW. So the distributed hydrological model can coupled to the groundwater model. The research proved that base on couple the distributed surface water model the groundwater model’s results are notability improved. The example is Jiyuan basin where use the distributed model coupled to the groundwater model. On the base of the coupled model applied to Jiyuan basin groundwater modeling. The paper estimates the groundwater change in the study area. Then, by use the water resources integrated planning results, the article calculate the basin groundwater can be development and utilization quantity and potential.

波动方程模拟与偏移的实用化方法研究

Along with the widespread and in-depth applications in petroleum prospecting and development, the seismic modeling and migration technologies are proposed with a higher requirement by oil industrial, and the related practical demand is getting more and more urgent. Based on theories of modeling and migration methods for wave equation, both related with velocity model, I thoroughly research and develop some methods for the goal of highly effective and practical in this dissertation. In the first part, this dissertation probes into the layout designing by wave equations modeling, focusing on the target-oriented layout designing method guided by wave equation modeling in complicated structure areas. It is implemented by using the fourth order staggered grid finite difference (FD) method in velocity-stress 2D acoustic wave equations plus perfectly matched layer (PML) absorbing boundary condition. To design target-oriented layout: (a) match the synthetic record on the surface with events of subsurface structures by analyzing the snapshots of theoretical model; (b) determine the shot-gather distance by tracking the events of target areas and measuring the receiving range when it reaches the surface; (c) restrict the range of valid shot-gather distance by drawing seismic windows in single shot records; (d) choose the best trace distance by comparing the resolution of prospecting targets from the simulated records with different trace distance. Eventually, we obtained the observation system parameters, which achieve the design requirements. In the second part, this dissertation presents the practical method to improve the 3D Fourier Finite Difference (FFD) migration, and carefully analyzes all the factors which influence 3D FFD migration’s efficiency. In which, one of the most important parameters of migration is the extrapolating step. This dissertation presents an efficient 3D FFD migration algorithm, which use FFD propagator to extrapolate wavefields over big layers, and use Born-Kirchhoff interpolator to image wavefields over small layers between the big ones. Finally, I show the effectiveness of this hybrid migration method by comparing migration results from 3D SEG/EAGE model with different methods.

数据驱动型多次波衰减方法的研究

In the practical seismic profile multiple reflections tend to impede the task of even the experienced interpreter in deducing information from the reflection data. Surface multiples are usually much stronger, more broadband, and more of a problem than internal multiples because the reflection coefficient at the water surface is much larger than the reflection coefficients found in the subsurface. For this reason most attempts to remove multiples from marine data focus on surface multiples, as will I. A surface-related multiple attenuation method can be formulated as an iterative procedure. In this essay a fully data-driven approach which is called MPI —multiple prediction through inversion (Wang, 2003) is applied to a real marine seismic data example. This is a pretty promising scheme for predicting a relative accurate multiple model by updating the multiple model iteratively, as we usually do in a linearized inverse problem. The prominent characteristic of MPI method lie in that it eliminate the need for an explicit surface operator which means it can model the multiple wavefield without any knowledge of surface and subsurface structures even a source signature. Another key feature of this scheme is that it can predict multiples not only in time but also in phase and in amplitude domain. According to the real data experiments it is shown that this scheme for multiple prediction can be made very efficient if a good initial estimate of the multiple-free data set can be provided in the first iteration. In the other core step which is multiple subtraction we use an expanded multi-channel matching filter to fulfil this aim. Compared to a normal multichannel matching filter where an original seismic trace is matched by a group of multiple-model traces, in EMCM filter a seismic trace is matched by not only a group of the ordinary multiple-model traces but also their adjoints generated mathematically. The adjoints of a multiple-model trace include its first derivative, its Hilbert transform and the derivative of the Hilbert transform. The third chapter of the thesis is the application for the real data using the previous methods we put forward from which we can obviously find the effectivity and prospect of the value in use. For this specific case I have done three group experiments to test the effectiveness of MPI method, compare different subtraction results with fixed filter length but different window length, invest the influence of the initial subtraction result for MPI method. In terms of the real data application, we do fine that the initial demultiple estimate take on a great deal of influence for the MPI method. Then two approaches are introduced to refine the intial demultiple estimate which are first arrival and masking filter respectively. In the last part some conclusions are drawn in terms of the previous results I have got.

尾波干涉法精确测量地震波速变化的方法及实验研究

The coda of seismic waves consists of that part of the signal after the directly arrivials. In a finite medium, or in one that is strongly heterogeneous, the coda is dominated by waves which have repeatedly sampled the medium. Small changes in a medium which may have no detectable influence on the first arrivals are amplified by this repeated sampling and may thus be detectable in the coda. Because of this, coda wave is widely used in detecting micro variations in medium。 In this paper, we give a general view of the theory and application of coda wave, especially coda wave interferometry. We focus on discussing the application of coda wave interferometry on data source of active situ experiment。 First, we apply coda wave interferometry in a short time period situ experiment which last for three days. We also apply the method of coda wave interferometry in a situ experiment which last for one month. Daily circle variations of seismic velocity around the experiment site were obtained, and we also observed that the velocity variations in the experiment site have a significant correlation with the environment factors, including air temperature, barometric pressure, solid earth tide and the level of rainfall. We find that the velocity variation during this period is up to 10-3. The relationship between velocity variation and changes in air temperature, barometric pressure and solid earth tide was analyzed with least square linear fitting .The velocity has no dependence on the air temperature. But velocity has a change of 10-6--10-7 when the barometer or earth tide change per Pa. Generally, we conclude the work and results of previous researchers, and we also display our works and results. We hopes to contribute to the future research of coda wave interferometry.

地震采集观测系统的共聚焦分析

Focal beam analysis is a method for assessment of acquisition geometries that is directly linked to pre-stack migration. About dealing with the complex subsurface structures, the conventional survey design methods which do not take into account the subsurface are no longer valid. Based on the Fourier finite-difference (FFD) large-step wave field extrapolation and Born-Kirchhoff (BK) small-step wavefield interpolation, the thesis presents a rapid resolution analysis of 3D seismic survey design by focal beams in complicated media. Subsequently, The SEG/EAEG salt model is used to illustrate the method. Based on the focal beam resolution definition, each kind of influence factor is discussed. The focal beam analysis usually is carried out in a single frequency, but the actual seismic waves always contain a frequency bandwidth. In this thesis, theoretical relationship between focal beam analysis and frequency is derived. Since the effects of focal beam analysis are linear with frequency simply, the multi-frequency focal beam analysis using interpolation is developed. At the same time, the resolution of different frequency bandwidth is interconvertible in accordance with Signal uncertainty principle. The resolution of all frequency bands can be calculated by using only a few focal beam analysis for a seismic survey. In the last section of this thesis, I propose a new approach to predicting acquisition footprint, based on the assumption of Common-Middle-Point stack without constructing a special velocity model. The approach is a simplistic analytical method in which the acquisition footprint pattern is a weighted, linear summation of limited-offset fold-of-stack plots. Because the value of acquisition can be got by quantificational and rapidly calculating, we can exactly do a comparative analysis among different plans of seismic survey by this method.

表面相关多次波衰减方法与应用研究

The function of seismic data in prospecting and exploring oil and gas has exceeded ascertaining structural configuration early. In order to determine the advantageous target area more exactly, we need exactly image the subsurface media. So prestack migration imaging especially prestack depth migration has been used increasingly widely. Currently, seismic migration imaging methods are mainly based on primary energy and most of migration methods use one-way wave equation. Multiple will mask primary and sometimes will be regarded as primary and interferes with the imaging of primary, so multiple elimination is still a very important research subject. At present there are three different wavefield prediction and subtraction methods: wavefield extrapolation; feedback loop; and inverse-scattering series. I mainly do research on feedback loop method in this paper. Feedback loop method includs prediction and subtraction.Currently this method has some problems as follows. Firstly, feedback loop method requires the seismic data used to predict multiple is full wavefield data, but usually the original seismic data don’t meet this assumption, so seismic data must be regularized. Secondly, Multiple predicted through feedback loop method usually can’t match the real multiple in seismic data and they are different in amplitude, phase and arrrival time. So we need match the predicted multiple and that in seismic data through estimating filtering factors and subtract multiple from seismic data. It is the key for multiple elimination how to select a correct matching filtering method. There are many matching filtering methods and I put emphasis on Least-square adaptive matching filtering and L1-norm minimizing adaptive matching filtering methods. Least-square adaptive matching filtering method is computationally very fast, but it has two assumptions: the signal has minimum energy and is orthogonal to the noise. When seismic data don’t meet the two assumptions, this method can’t get good matching results and then can’t attenuate multiple correctly. L1-norm adaptive matching filtering methods can avoid these two assumptions and then get good matching results, but this method is computationally a little slow. The results of my research are as follows: 1. Proposed a method that interpolates seismic traces based on F-K migration and demigration. The main advantage of this method is that it can interpolate seismic traces in any offsets. It shows this method is valid through a simple model. 2. Comparing different Least-square adaptive matching filtering methods. The results show that equipose multi-channel adaptive matching filtering methods can get better results of multiple elimination than other matcing methods through three model data and two field data. 3. Proposed equipose multi-channel L1-norm adaptive matching filtering method. Because L1-norm is robust to large amplitude differences, there are no assumption on the signal has minimum energy and orthogonality, this method can get better results of multiple elimination. 4. Research on multiple elimination in inverse data space. The method is a new multiple elimination method and it is different from those methods mentioned above.The advantages of this method is that it is simple in theory and no need for the adaptive subtraction and computationally very fast. The disadvantage of this method is that it is not stabilized in its solution. The results show that equipose multi-channel and equipose pesudo-multi-channel least-square matching filtering and equipose multi-channel and equipose pesudo-multi-channel L1-norm matching filtering methods can get better results of multiple elimination than other matcing methods through three model data and many field data.

叠前时间偏移：各向异性、多次波衰减及应用

At present, in order to image complex structures more accurately, the seismic migration methods has been developed from isotropic media to the anisotropic media. This dissertation develops a prestack time migration algorithm and application aspects for complex structures systematically. In transversely isotropic media with a vertical symmetry axis (VTI media), the dissertation starts from the theory that the prestack time migration is an approximation of the prestack depth migration, based on the one way wave equation and VTI time migration dispersion relation, by combining the stationary-phase theory gives a wave equation based VTI prestack time migration algorithm. Based on this algorithm, we can analytically obtain the travel time and amplitude expression in VTI media, as while conclude how the anisotropic parameter influence the time migration, and by analyzing the normal moveout of the far offset seismic data and lateral inhomogeneity of velocity, we can update the velocity model and estimate the anisotropic parameter model through the time migration. When anisotropic parameter is zero, this algorithm degenerates to the isotropic time migration algorithm naturally, so we can propose an isotopic processing procedure for imaging. This procedure may keep the main character of time migration such as high computational efficiency and velocity estimation through the migration, and, additionally, partially compensate the geometric divergence by adopting the deconvolution imaging condition of wave equation migration. Application of this algorithm to the complicated synthetic dataset and field data demonstrates the effectiveness of the approach. In the dissertation we also present an approach for estimating the velocity model and anisotropic parameter model. After analyzing the velocity and anisotropic parameter impaction on the time migration, and based on the normal moveout of the far offset seismic data and lateral inhomogeneity of velocity, through migration we can update the velocity model and estimate the anisotropic parameter model by combining the advantages of velocity analysis in isotropic media and anisotropic parameter estimation in VTI media. Testing on the synthetic and field data, demonstrates the method is effective and very steady. Massive synthetic dataset、2D sea dataset and 3D field datasets are used for VTI prestack time migration and compared to the stacked section after NMO and prestack isotropic time migration stacked section to demonstrate that VTI prestack time migration method in this paper can obtain better focusing and less positioning errors of complicated dip reflectors. When subsurface is more complex, primaries and multiples could not be separated in the Radon domain because they can no longer be described with simple functions (parabolic). We propose an attenuating multiple method in the image domain to resolve this problem. For a given velocity model，since time migration takes the complex structures wavefield propagation in to account, primaries and multiples have different offset-domain moveout discrepancies, then can be separated using techniques similar to the prior migration with Radon transform. Since every individual offset-domain common-reflection point gather incorporates complex 3D propagation effects, our method has the advantage of working with 3D data and complicated geology. Testing on synthetic and real data, we demonstrate the power of the method in discriminating between primaries and multiples after prestack time migration, and multiples can be attenuated in the image space considerably.

松辽盆地南部大情字井地区CO2埋存体地质格架及储层建模

Increasing attentions have been paid to the subsurface geological storage for CO2 in view of the huge storage capacity of subsurface reservoirs. The basic requirement for subsurface CO2 storage is that the CO2 should be sequestrated as supercritical fluids (physical trapping), which may also interact with ambient reservoir rocks and formation waters, forming new minerals (chemical trapping). In order to the effective, durable and safe storage for CO2, enough storage space and stable sealing caprock with strong sealing capacity are necessitated, in an appropriate geological framework. Up till now, hydrocarbon reservoirs are to the most valid and appropriate CO2 storage container, which is well proven as the favorable compartment with huge storage capacity and sealing condition. The thesis focuses on two principal issues related to the storage and sealing capacity of storage compartment for the Qingshankou and Yaojia formations in the Daqingzijing block, Southern Songliao Basin, which was selected as the pilot well site for CO2-EOR storage. In the operation area, three facies, including deltaic plain, deltaic front and subdeep-deep lake facies associations, are recognized, in which 11 subfacies such as subaqueous distributary channel, river- mouth bar, interdistributary bay, sheet sandbody, crevasse splay and overflooding plain are further identified. These subfacies are the basic genetic units in the reservoir and sealing rocks. These facies further comprise the retrogradational and progradational depositional cycles, which were formed base- level rise and fall, respectively. During the regressive or lake lowstand stage, various sands including some turbidites and fans occurred mostly at the bottom of the hinged slope. During the progradation stage, these sands became smaller in size and episodically stepped backwards upon the slope, with greatly expanded and deeped lake. However, most of Cretaceous strata in the study area, localized in the basin centre under this stage, are mainly composed of grey or grizzly siltstones and grey or dark grey mudstones intercalated with minor fine sandstones and purple mudstones. On the base of borehole and core data, these siltstones are widespread, thin from 10 to 50 m thick, good grain sorting, and have relative mature sedimentary structures with graded bedding and cross- lamination or crossbeds such as ripples, which reflect strong hydrodynamic causes. Due to late diagenesis, pores are not widespread in the reservoirs, especially the first member of Qingshankou formation. There are two types of pores: primary pore and secondary cores. The primary pores include intergranular pores and micropores, and the secondary pores include emposieus and fracture pores. Throat channels related to pores is also small and the radius of throat in the first, second and third member of Qingshankou formation is only 0.757 μm, 0.802 μm and 0.631 μm respectively. In addition, based on analyzing the probability plot according to frequency of occurrence of porosity and permeability, they appear single- peaked distribution, which reflects strong hetero- geneity. All these facts indicate that the conditions of physical property of reservoirs are not better. One reason may be provided to interpret this question is that physical property of reservoirs in the study area is strong controlled by the depositional microfacies. From the statistics, the average porosity and permeability of microfacies such as subaqueous distributary channel, channel mouth bar, turbidites, is more than 9 percent and 1md respectively. On the contrary, the average porosity and permeability of microfacies including sand sheet, flagstone and crevasse splay are less than 9 percent and 0.2md respectively. Basically, different hydrodynamic environment under different microfacies can decide different physical property. According to the reservoir models of the first member of Qingshankou formation in the No. well Hei47 block, the character of sedimentary according to the facies models is accord to regional disposition evolution. Meantime, the parameter models of physical property of reservoir indicate that low porosity and low permeability reservoirs widespread widely in the study area, but the sand reservoirs located in the channels are better than other places and they are the main sand reservoirs. The distribution and sealing ability of fault- fractures and caprock are the key aspects to evaluate the stable conditions of compartments to store CO2 in the study area. Based on the core observation, the fractures widespread in the study area, especially around the wells, and most of them are located in the first and second member of Qingshankou formation, almost very few in the third member of Qingshankou formation and Yaojia formation instead. In addition, analyzing the sealing ability of eleven faults in the three-dimensional area in the study area demonstrates that most of faults have strong sealing ability, especially in the No. well Hei56 and Qing90-27. To some extent, the sealing ability of faults in the No. well Hei49, Qing4-6 and Qing84-29 are worse than others. Besides, the deposition environment of most of formations in the study area belongs to moderately deep and deep lake facies, which undoubtedly take advantage to caprocks composed of mudstones widespread and large scale under this deposition environment. In the study area, these mudstones distribute widely in the third member of Qingshankou formation, Yaojia and Nenjiang formation. The effective thickness of mudstone is nearly ~550m on an average with few or simple faults and fractures. In addition, there are many reservoir beds with widely- developed insulated interbeds consist of mudstones or silty mudstone, which can be the valid barrier to CO2 upper movement or leakage through diffusion, dispersion and convection. Above all, the closed thick mud caprock with underdeveloped fractures and reservoir beds can be taken regard as the favorable caprocks to provide stable conditions to avoid CO2 leakage.