多尺度分解与相位反褶积提高海上地震资料分辨率和信噪比

Offshore seismic exploration is full of high investment and risk. And there are many problems, such as multiple. The technology of high resolution and high S/N ratio on marine seismic data processing is becoming an important project. In this paper, the technology of multi-scale decomposition on both prestack and poststack seismic data based on wavelet and Hilbert-Huang transform and the theory of phase deconvolution is proposed by analysis of marine seismic exploration, investigation and study of literatures, and integration of current mainstream and emerging technology. Related algorithms are studied. The Pyramid algorithm of decomposition and reconstruction had been given by the Mallat algorithm of discrete wavelet transform In this paper, it is introduced into seismic data processing, the validity is shown by test with field data. The main idea of Hilbert-Huang transform is the empirical mode decomposition with which any complicated data set can be decomposed into a finite and often small number of intrinsic mode functions that admit well-behaved Hilbert transform. After the decomposition, a analytical signal is constructed by Hilbert transform, from which the instantaneous frequency and amplitude can be obtained. And then, Hilbert spectrum. This decomposition method is adaptive and highly efficient. Since the decomposition is based on the local characteristics of the time scale of data, it is applicable to nonlinear and non-stationary processes. The phenomenons of fitting overshoot and undershoot and end swings are analyzed in Hilbert-Huang transform. And these phenomenons are eliminated by effective method which is studied in the paper. The technology of multi-scale decomposition on both prestack and poststack seismic data can realize the amplitude preserved processing, enhance the seismic data resolution greatly, and overcome the problem that different frequency components can not restore amplitude properly uniformly in the conventional method. The method of phase deconvolution, which has overcome the minimum phase limitation in traditional deconvolution, approached the base fact well that the seismic wavelet is phase mixed in practical application. And a more reliable result will be given by this method. In the applied research, the high resolution relative amplitude preserved processing result has been obtained by careful analysis and research with the application of the methods mentioned above in seismic data processing in four different target areas of China Sea. Finally, a set of processing flow and method system was formed in the paper, which has been carried on in the application in the actual production process and has made the good progress and the huge economic benefit.

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

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.

井间地震速度层析成像方法与软件系统

Crosshole Seismic tomography has been broadly studied and applied in the fields of resource exploration and engineering exploration because of its special observing manner and better resolution than normal seismic exploration. This thesis will state the theory and method of Crosshole Seismic tomography. Basing on the previous studies，the thesis studied the initial velocity model，ray-tracing method, and developed the three-dimension tomography software. All the cells that a ray passes through are of the same velocities if the paths from transmitters to receivers are straight. The cells that the each ray passes through are recorded, and rays that pass through each cell are calculated. The ray average velocity which passes through a cell is set as the cell velocity. Analogously we can make a initial node velocity model because the velocity sum is calculated on the all cells which own to a certain node, and the cell number is summed about each nodes，the ratio of the velocity sum to the all cells number is set as the node velocity. The inversion result from the initial node velocity model is better than that of the average velocity model. Ray-bending and Shortest Path for Rays (SPR) have shortcomings and limitations respectively. Using crooked rays obtained from SPR rather than straight lines as the starting point can not only avoid ray bending converging to the local minimum travel time path, but also settle the no smooth ray problem obtained by SPR. The hybrid method costs much computation time, which is roughly equal to the time that SPR expends. The Delphi development tool based on the Object Pascal language standard has an advantage of object-oriented. TDTOM (Three Dimensions Tomography) was developed by using Delphi from the DOS version. Improvement on the part of inversion was made, which bring faster convergence velocity. TDTOM can be used to do velocity tomography from the first arrival travel time of the seismic wave, and it has the good qualities of friendly user interface and convenient operation. TDTOM is used to reconstruct the velocity image for a set of crosshole data from Karamay Oil Field. The geological explanation is then given by comparing the inversion effects of different ray-tracing methods. High velocity zones mean the cover of oil reservoir, and low velocity zones correspond to the reservoir or the steam flooding layer.

晚中新世以来中国北方风成沉积的磁性地层学和沉积学研究

The past two decades have witnessed an unprecedented growth of interest in the palaeoenvironmental significance of the Pleistocene loess deposits in northern China. However, it is only several years ago that the Tertiary red clay sequence underlying Pleistocene loess attracted much attention. One of the major advances in recent studies of eolian deposits on the Loess Plateau is the verification of the eolian origin for the Tertiary red clay sediments. The evidence of the eolian origin for the red clay is mainly from geochemical and sedimentological studies. However, sedimentological studies of the red clay deposits are still few compared with those of the overlying loess sediments. To date, the red clay sections located near Xifeng, Baoji, Lantian, Jiaxian, and Lingtai have been studied, with an emphasis on magnetostratigraphy. These sections have a basal age ranging from ～4.3 Ma to ～7.0 Ma. The thickness of the sections varies significantly, depending perhaps on the development of local geomorphological conditions and the drainage system. Although the stratigraphy of the red clay sections has been recorded in some detail, correlation of the red clay sequences has not yet been undertaken. Geological records (Sun J. et al., 1998) have shown that during glacial periods of the Quaternary the deserts in northem China were greatly expanded compared with modern desert distribution. During interglacial periods, desert areas contracted and retreated mostly to northwestern China because of the increase in inland penetration of monsoonal precipitation. According to pedogenic characteristics of the red clay deposits, the climatic conditions of the Loess Plateau is warmer and wetter generally in the Neogene than in the late Pleistocene. Panicle analyses show that grain size distribution of the red clay sequence is similar to that of the paleosols in the Pleistocene loess record, thus implying a relatively remote provenance of the red clay materials. However, the quantitative or semiquantitative estimates of the distance from the source region to the Loess Plateau during the red clay development remains to be investigated. In this study, magnetostratigraphic and sedimentological studies are conducted at two thick red clay sequences-Jingchuan and Lingtai section. The objectives of these studies are focused on further sedimentological evidence for the eolian origin of the red clay, correlation of red clay sequences, provenance of the red clay, and the palaeoclimate reconstruction in the Neogene. Paleomagnetic studies show that the Jingchuan red clay has a basal age of 8.0 Ma, which is 1 million years older than the previously studied Lingtai section. The Lingtai red clay sequence was divided into five units on the basis of pedogenica characteristics (Ding et al., 1999a). The Jingchuan red clay sequence, however, can be lithologically divided into six units according to field observations. The upper five units of the Jingchuan red clay can generally correlate well with the five units of the Lingtai red clay. Comparison of magnetic susceptibility and color reflectance records of four red clay sections suggests that the Lingtai red clay sequence can be the type-section of the Neogene red clay deposits in northern China. Pleistocene loess and modem dust deposits have a unimodal grain-size distribution. The red clay sediments at Jingchuan and Lingtai also have a unimodal grain-size distribution especially similar to the paleosols in the Pleistocene loess record. Sedimentological studies of a north-south transect of loess deposits above S2 on the Loess Plateau show that loess deposits had distinct temporal and spatial sedimentary differentiation. The characteristics of such sedimentary differentiation can be well presented in a triangular diagram of normalized median grain size, normalized skewness, and normalized kurtosis. The triangular diagrams of the red clay-loess sequence at Lingtai and Jingchuan indicate that loess-paleosol-red clay may be transported and sorted by the same agent wind, thus extending the eolian record in the Loess Plateau from 2.6 Ma back to about 8.0 Ma. It has been recognized that during the last glacial maximum (LGM) the deserts in northern China had a distribution similar to the present, whereas during the Holocene Optimum the deserts retreated to the area west of the Helan Mountains. Advance-retreat cycles of the deserts will lead to changes in the distance of the Loess Plateau to the dust source regions, thereby controlling changes in grain size of the loess deposited in a specific site. To observe spatial changes in sedimentological characteristics of loess during the last glacial-interglacial cycle, the texture of loess was measured along the north-south transect of the Loess Plateau. Since the southern margin of the Mu Us desert during the LGM is already known, several models of grain size parameters versus the minimum distance from the source region to depositional areas were developed. According to these semiquantitative models, the minimum distance from the source region to Lingtai and Jingchuan areas is about 600 km during the Neogene. Therefore the estimated provenance of the Tertiary red clay deposits is the areas now occupied by the Badain Jaran desert and arid regions west of it. The ratio of the free iron to total iron concentration attests to being a good proxy indicator for the summer monsoon evolution. The Lingtai Fe_20_3 ratio record shows high values over three time intervals: 4.8-4.1 Ma, 3.4-2.6 Ma, and during the interglacial periods of the past 0.5 Ma. The increase in summer monsoon intensity over the three intervals also coincides with the well-developed soil characteristics. It is therefore concluded that the East-Asia summer monsoon has experienced a non-linear evolution since the late Miocene. In general, the East Asia summer monsoon was stronger in Neogene than in Quaternary and the strongest East Asia summer monsoon may occur between 4.1 and 4.8 Ma. The relatively small ice volume and high global temperature may be responsible for the strong summer monsoon during the early Pliocene.