Observation and theoretical analysis for the sand-waves migration in the North Gulf of South China Sea

According to the environmental characteristic of the north gulf of South China Sea, a quasi-3D mechanics model has been built for simulating the small scale sand-waves migration in the seas of southwest of Hainan Island. Based on the submarine micro-geomorphic data induced by multi-beam system and hydrographic survey record, the migrations of the sand-waves in the study area are predicted. The results show that calculation is consistent with the observation data in the groove of sand ridge, but not well in the crest of sand ridge. It is indicated that the mechanics model should be used to predict the migration of the small scale sand-waves which are dominated by bed load in the seas. This paper is very meaningful to project the route of submarine pipeline.

Robust Guaranteed Cost Control for Yaw Control of Helicopter

针对具有时变不确定性且不确定性界为椭球的线性系统提出了一种新的具有自适应机制的鲁棒保性能控制器设计方法。首先,引入一个具有可由自适应律在线调整的可调参数的目标模型,通过该参数来保证由目标模型与被控模型所获得的误差系统渐近稳定。结合保证目标模型稳定性的设计,最终形成保证闭环系统稳定且控制器增益仿射依赖于可调参数的鲁棒保性能跟踪控制器。应用于安装在试验平台上的小型直升机航向控制中,仿真试验表明了该方法的有效性。

自适应鲁棒H2控制及其在小型无人直升机航向控制中的应用

提出了一个新的鲁棒控制方法并将其应用到小型无人直升机的航向控制中.给出了直升机的非线性航向动力学模型并进行简化,然后线性化为具有仿射不确定性的线性模型.针对这个线性模型,基于线性矩阵不等式方法提出了具有自适应机制的鲁棒H2反馈控制器.这种设计方法可以降低固定增益控制器所固有的保守性,提高控制效果.数字仿真验证了设计方法的可行性。

模型直升机航向自适应保性能控制

提出了一种新的模型直升机航向控制算法。针对具有模型不确定性的直升机航向线性模型,提出了一种具有自适应机制的最优保性能控制器。该控制策略通过引入自适应机制降低固定增益控制器所固有的保守性,并且控制器的反馈增益应用线性矩阵不等式(LMIs)方法解得。理论分析和数字仿真表明所设计的控制器具有良好的鲁棒稳定性能。

碳酸盐岩缝洞型储层地球物理响应特征及预测方法研究

Proven by the petroleum exploration activities, the karsts-fissure reservoir in carbonate rocks is significant to find out the large scale oil & gas field. They are made up of the four reservoir types: karsts-cave, karsts-crack, crack-cave and fracture-pore-cave. Each reservoir space and each reservoir bed has different features of reservoir heterogeneity and small scale of pore-crack-cave. The fracture-cave reservoir in carbonate rocks is characteristic by multi-types and long oiliness well. The reservoir shape is controlled by the irregular pore-crack-cave. The development level of fracture and karst-cave is the key element of hydrocarbon enriching, high productivity and stable production. However, most of Carbonate formation are buried deeply and the signal-ration-noise of seismic reflection are very low. It is reason why the fracture-cave reservoir are difficult to be predicted effectively. In terms of surveyed and studied lots of the former research outcome, The author applied the methods of synthetical reservoir geophysical prediction from two ways including macrosopic and microcomic technics in terms of the reservoir-cap condition, geophysics and geology feature and difficulty of prediction in carbonate rocks. It is guiden by the new ideas of stratigraphy, sedimentology, sedimentography, reservoir geology and karst geology. The geophysics technology is key technics. In aspects of macroscopic studies, starting off the three efficiencies of controlling the reservoir distribution including sedimental facies, karst and fracture, by means of comprehensive utilization of geology, geophysics, boring well and well log, the study of reservoir features and karst inside story are developed in terms of data of individual well and multiple well. Through establishing the carbonate deposition model, karstic model and fracture model, the macro-distribution laws of carbonatite are carried out by the study of coherence analysis, seismic reflection feature analysis and palaeotectonics analysis. In aspects of microcosmic studies, starting off analysis in reservoir geophysical response feature of fracture and karst-cave model according to guidance of the macroscopic geological model in carbonate reservoir, the methods of the carbonate reservoir prediction are developed by comprehensively utilization of seismic multi-attribution intersection analysis, seismic inversion restricted by log, seismic discontinuity analysis, seimic spectrum attenuation gradient, moniliform reflection feature analysis and multiparameter karst reservoir appraisement.Through application of carbonate reservoir synthetical geophysics prediction, the author r successfully develops the beneficial reservoir distribution province in Ordovician of Katake block 1in middle Tarim basin. The fracture-cave reservoir distributions are delineated. The prospect direction and favorable aims are demonstrated. There are a set of carbonate reservoir prediction methods in middle Tarim basin. It is the favorable basic technique in predicting reservoir of the Ordovician carbonate in middle Tarim. Proven by exploration drilling, the favorable region of moniliform reflection fracture and pore-cave and cave-fracture in lower-middle Ordovician are coincidence with the region of hydrocarbon show. It’s indicated that the reservoir prediction methods described in the study of Ordovician carbonate formation are feasible practicably.

柴达木盆地地震勘探关键技术研究

The rugged surface topography determined the seismic data acquisition construction conditions and the seismic wave explosive and receiver quality in Qaidam Basin. This dissertation systematically researched the seismic acquisition, imaging process and the attribute analysis techniques of complicated oil and gas reservoir. The main research achievements and cognitions are as follows: 1. Through the stimulation effects research and analysis from the aspect of lithologic water-containing differences, it’s specific that stable hydrous sand layer can effectively enhance the stimulation effects combined with the corresponding field tests. The seismic data S/N ratio has been improved due to the combination explosive stimulation. Through the fold number and maximum offset analyses of target horizon, the complicated geometry has been optimized and the S/N ratio of seismic data has been improved, which made an important basis for improvement of 3D seismic data. 2. It has been proved that the first arrival refraction static correction method under the model constraint of fine surface survey is suitable to the Qaidam Basin of western areas by the real seismic data processing. Although the refraction horizon of near surface has some changes in a certain extent, it’s steady basically. The refraction horizon can be continuously traced in sections, so it’s qualified for the refraction static correction method on the whole. 3. The research is based on the curved-ray pre-stack time migration techniques of rough topography, and improved the imaging precision of complex areas. This techniques adopted the constant and variable velocity scanning mode and enhanced the velocity analysis precision. The 3D pre-stack time migration techniques reasonably solved the imaging and velocity multiple solutions problems of steep-dip faults and the intersections of horizontal layers. What’s more, fine velocity analysis and mute are very important to enhance the imaging precision of the seismic data in complicated Wunan areas. 4. The 3D seismic data edge-preserving processing methods have been realized due to the image process techniques. Because this method uses the large range filter, it can attenuate the noise maximally. The faults, break points, lithologic pinchout points and lithologic body of small scale such as river will not be influenced by blur because of the edge-preserving characterization of the method which is really an effective assistant technique of low S/N ratio seismic data attribute analysis. 5. The use of spectral decomposition technique can effectively identify the reservoirs. The special geology body which will not be identified (or without obvious characters) in the seismic profile may be found through the details changes of different frequencies in the amplitude profiles.

地震岩性与物性联合反演技术研究和应用

This paper studies how to more effectively invert seismic data and predict reservoir under complicated sedimentary environment, complex rock physical relationships and fewer drills in offshore areas of China. Based on rock physical and seismic amplitude-preserving process, and according to depositional system and laws of hydrocarbon reservoir, in the light of feature of seismic inversion methods present applied, series methods were studied. A joint inversion technology for complex geological condition had been presented, at the same time the process and method system for reservoir prediction had been established. This method consists four key parts. 1)We presented the new conception called generalized wave impedance, established corresponding inversion process, and provided technical means for joint inversion lithology and petrophysical on complex geological condition. 2)At the aspect of high-resolution nonlinear seismic wave impedance joint inversion, this method used a multistage nonlinear seismic convolution model rather than conventional primary structure Robinson seismic convolution model, and used Caianiello neural network implement inversion. Based on the definition of multistage positive and negative wavelet, it adopted both deterministic and statistical physical mechanism, direct inversion and indirect inversion. It integrated geological knowledge, rock physical theory, well data, and seismic data, and improved the resolution and anti-noise ability of wave impedence inversion. 3)At the aspect of high-resolution nonlinear reservoir physical property joint inversion, this method used nonlinear rock physical model which introduced convolution model into the relationship between wave impedance and porosity/clay. Through multistage decomposition, it handles separately the large- and small-scale components of the impedance-porosity/clay relationships to achieve more accurate rock physical relationships. By means of bidirectional edge detection with wavelets, it uses the Caianiello neural network to finish statistical inversion with combined applications of model-based and deconvolution-based methods. The resulted joint inversion scheme can integrate seismic data, well data, rock physical theory, and geological knowledge for estimation of high-resolution petrophysical parameters. 4)At the aspect of risk assessment of lateral reservoir prediction, this method integrated the seismic lithology identification, petrophysical prediction, multi-scale decomposition of petrophysical parameters, P- and H-spectra, and the match relationship of data got from seismics, well logging and geology. It could describe the complexity of medium preferably. Through applications of the joint inversion of seismic data for lithologic and petrophysical parameters in several selected target areas, the resulted high-resolution lithologic and petrophysical sections(impedance, porosity, clay) show that the joint inversion can significantly improve the spatial description of reservoirs in data sets involving complex deposits. It proved the validity and practicality of this method adequately.

油气勘探中的弱磁异常提取与区域场分析

Petroleum and natural gas is an important strategic resources. The short of the reserves will block the development of economy and threaten the safety of nation, along with the main oil fields of our country coming to the height of power and splendor of the exploitation and exploration. Therefore, it makes a great sense to inaugurate new explorative field and increase the reserves of petroleum and natural gas. Magnetic exploration is a main method of geophysics exploration. the developing observation apparatus and the perfect processing method provide wide space for magnetic exploration in these years. The method of magnetic bright spot is an application of magnetic exploration. The vertical migration of the hydrocarbon changes physical and chemical environment above the hydrocarbon reservoir, the new environment make tervalent iron translate into bivalent iron, that produce small scale magnetic anomaly, that is magnetic bright spot. The method of magnetic bright spot explores oil and gas field by the relation between the hydrocarbon and magnetic anomaly. This paper systemically research to pick-up and identify magnetic bright spot combining an oil field item, then point out advantaged area. In order to test the result, the author use the seismic information to superpose the magnetic bright spot, that prove the magnetic bright spot is reliable. then, the author complete a software to pick and identify the magnetic bright spot. The magnetic basement is very important to research forming and evolvement of the basin, especially, it is a crucial parameter of exploring residual basin in the research on pre-Cenozoic residual. This paper put forward a new method to inverse the interface of the magnetic layer on the basis of previous work, that is the method of separation of magnetic field step by step. The theory of this method is to translate the result of magnetic layer fluctuation to the result of magnetization density change, and the magnetic layer is flat, the paper choose thickness of magnetic layer as unit thickness, and define magnetic layer as a unit-thickness layer in order to convenient calculation, at the same time, define the variational magnetization density as equivalent magnetic density. Then we translate the relation between magnetic field and layer fluctuation to the relation between magnetic field and equivalent magnetic density, then, we can obtain the layer fluctuation through calculating equivalent magnetic density. Contrast to conventional parker method, model experimentation and example checkout prove this method is effective. The merit of this method is to avoid flat result in a strongly fluctuant area because of using a uniform average depth, the result of this method is closer to the fact, and this method is to inverse equivalent magnetic density, then translate equivalent magnetic density to layer fluctuation, this lays a foundation to inverse variational magnetic density in the landscape orientation and portrait.

基于随机介质理论的多尺度储层非均质性定量描述

Application of long-term exploration for oil and gas shows that the reservoir technology of prediction is one of the most valuable methods. Quantitative analysis of reservoir complexity is also a key technology of reservoir prediction. The current reservoir technologies of prediction are based on the linear assumption of various physical relationships. Therefore, these technologies cannot handle complex reservoirs with thin sands, high heterogeneities in lithological composition and strong varieties in petrophysical properties. Based on the above-mentioned complex reservoir, this paper conducts a series of researches. Both the comprehending and the quantitative analysis of reservoir heterogeneities have been implemented using statistical and non-linear theories of geophysics. At the beginning, the research of random media theories about reservoir heterogeneities was researched in this thesis. One-dimensional (1-D) and two-dimensional (2-D) random medium models were constructed. The autocorrelation lengths of random medium described the mean scale of heterogeneous anomaly in horizontal and deep directions, respectively. The characteristic of random medium models were analyzed. We also studied the corresponding relationship between the reservoir heterogeneities and autocorrelation lengths. Because heterogeneity of reservoir has fractal nature, we described heterogeneity of reservoir by fractal theory based on analyzing of the one-dimensional (1-D) and two-dimensional (2-D) random medium models. We simulated two-dimensional (2-D) random fluctuation medium in different parameters. From the simulated results, we can know that the main features of the two-dimensional (2-D) random medium mode. With autocorrelation lengths becoming larger, scales of heterogeneous geologic bodies in models became bigger. In addition, with the autocorrelation lengths becoming very larger, the layer characteristic of the models is very obvious. It would be difficult to identify sandstone such as gritstone, clay, dense sandstone and gas sandstone and so on in the reservoir with traditional impedance inversion. According to the obvious difference between different lithologic and petrophysical impedance, we studied multi-scale reservoir heterogeneities and developed new technologies. The distribution features of reservoir lithological and petrophysical heterogeneities along vertical and transverse directions were described quantitatively using multi-scale power spectrum and heterogeneity spectrum methods in this paper. Power spectrum (P spectrum) describes the manner of the vertical distribution of reservoir lithologic and petrophysical parameters and the large-scale and small-scale heterogeneities along vertical direction. Heterogeneity spectrum (H spectrum) describes the structure of the reservoir lithologic and petrophysical parameters mainly, that is to say, proportional composition of each lithological and petrophysical heterogeneities are calculated in this formation. The method is more reasonable to describe the degree of transverse multi-scale heterogeneities in reservoir lithological and petrophysical parameters. Using information of sonic logs in Sulige oil field, two spectral methods have been applied to the oil field, and good analytic results have been obtained. In order to contrast the former researches, the last part is the multi-scale character analysis of reservoir based on the transmission character of wave using the wavelet transform. We discussed the method applied to demarcate sequence stratigraphy and also analyzed the reservoir interlayer heterogeneity.

地震层析成像与天山地区的壳幔深部结构研究

The Tien Shan is the most prominent intracontinental mountain belt on the earth. The active crustal deformation and earthquake activities provide an excellent place to study the continental geodynamics of intracontinental mountain belt. The studies of deep structures in crust and upper mantle are significantly meaningful for understanding the geological evolution and geodynamics of global intracontinental mountain belts. This dissertation focuses on the deep structures and geodynamics in the crust and upper mantle in the Tien Shan mountain belt. With the arrival time data from permanent and temporal seismic stations located in the western and central Tien Shan, using seismic travel time tomographic method, we inversed the P-wave velocity and Vp/Vs structures in the crust and uppermost mantle, the Pn and Sn velocities and Pn anisotropic structures in the uppermost mantle, and the P-wave velocity structures in the crust and mantle deep to 690km depth beneath the Tien Shan. The tomographic results suggest that the deep structures and geodynamics have significant impacts not only on the deformations and earthquake activities in the crust, but also on the mountain building, collision, and dynamics of the whole Tien Shan mountain belt. With the strongly collision and deformations in the crust, the 3-D P-wave velocity and Vp/Vs ratio structures are highly complex. The Pn and Sn velocities in the uppermost mantle beneath the Tien Shan, specially beneath the central Tien Shan, are significantly lower than the seismic wavespeed beneath geological stable regions. We infer that the hot upper mantle from the small-scale convection could elevate the temperature in the lower crust and uppermost mantle, and partially melt the materials in the lower crust. The observations of low P-wave and S-wave velocities, high Vp/Vs ratios near the Moho and the absences of earthquake activities in the lower crust are consistent with this inference. Based on teleseismic tomography images of the upper mantle beneath the Tien Shan, we infer that the lithosphere beneath the Tarim basin has subducted under the Tien Shan to depths as great as 500 km. The lithosphere beneath the Kazakh shield may have subducted to similar depths in the opposite direction, but the limited resolution of this data set makes this inference less certain. These images support the plate boundary model of converge for the Tien Shan, as the lithospheres to the north and south of the range both appear to behave as plates.

复杂介质边界元法数值模拟研究

A major impetus to study the rough surface and complex structure in near surface model is because accuracy of seismic observation and geophysical prospecting can be improved. Wave theory study about fluid-satuated porous media has important significance for some scientific problems, such as explore underground resources, study of earth's internal structure, and structure response of multi-phase porous soil under dynamic and seismic effect. Seismic wave numerical modeling is one of the effective methods which understand seismic propagation rules in complex media. As a numerical simulation method, boundary element methods had been widely used in seismic wave field study. This paper mainly studies randomly rough surface scattering which used some approximation solutions based on boundary element method. In addition, I developed a boundary element solution for fluid saturated porous media. In this paper, we used boundary element methods which based on integral expression of wave equation to study the free rough surface scattering effects of Kirchhoff approximation method, Perturbation approximation method, Rytov approximation method and Born series approximation method. Gaussian spectrum model of randomly rough surfaces was chosen as the benchmark model. The approximation methods result were compared with exact results which obtained by boundary element methods, we study that the above approximation methods were applicable how rough surfaces and it is founded that this depends on and ( here is the wavenumber of the incident field, is the RMS height and is the surface correlation length ). In general, Kirchhoff approximation which ignores multiple scatterings between any two surface points has been considered valid for the large-scale roughness components. Perturbation theory based on Taylor series expansion is valid for the small-scale roughness components, as and are .Tests with the Gaussian topographies show that the Rytov approximation methods improves the Kirchhoff approximation in both amplitude and phase but at the cost of an extra treatment of transformation for the wave fields. The realistic methods for the multiscale surfaces come with the Born series approximation and the second-order Born series approximation might be sufficient to guarantee the accuracy of randomly rough surfaces. It could be an appropriate choice that a complex rough surface can be divided into large-, medium-, and small-scale roughness components with their scattering features be studied by the Kirchhoff or Rytov phase approximations, the Born series approximation, and the perturbation theory, respectively. For this purpose, it is important to select appropriate parameters that separate these different scale roughness components to guarantee the divided surfaces satisfy the physical assumptions of the used approximations, respectively. In addition, in this paper, the boundary element methods are used for solving the porous elastic wave propagation and carry out the numerical simulation. Based on the fluid-saturated porous model, this paper analyses and presents the dynamic equation of elastic wave propagation and boundary integral equation formulation of fluid saturated porous media in frequency domain. The fundamental solutions of the elastic wave equations are obtained according to the similarity between thermoelasticity and poroelasticity. At last, the numerical simulation of the elastic wave propagation in the two-phase isotropic media is carried out by using the boundary element method. The results show that a slow quasi P-wave can be seen in both solid and fluid wave-field synthetic seismograms. The boundary element method is effective and feasible.

特大地震激发地球自由振荡的数值模拟研究

Large earthquakes, such as the Chile earthquake in 1960 and the Sumatra-Andaman earthquake on Dec 26, 2004 in Indonesia, have generated the Earth’s free oscillations. The eigenfrequencies of the Earth’s free oscillations are closely related to the Earth’s internal structures. The conventional methods, which mainly focus on calculating the eigenfrequecies by analytical ways, and the analysis on observations can not easily study the whole processes from earthquake occurrence to the Earth’s free oscillation inspired. Therefore, we try to use numerical method incorporated with large-scale parallel computing to study on the Earth’s free oscillations excited by giant earthquakes. We first give a review of researches and developments of the Earth’s free oscillation, and basical theories under spherical coordinate system. We then give a review of the numerical simulation of seismic wave propagation and basical theories of spectral element method to simulate global seismic wave propagation. As a first step to study the Earth’s free oscillations, we use a finite element method to simulate the propagation of elastic waves and the generation of oscillations of the chime bell of Marquis Yi of Zeng, by striking different parts of the bell, which possesses the oval crosssection. The bronze chime bells of Marquis Yi of Zeng are precious cultural relics of China. The bells have a two-tone acoustic characteristic, i.e., striking different parts of the bell generates different tones. By analysis of the vibration in the bell and the spectrum analysis, we further help the understanding of the mechanism of two-tone acoustic characteristics of the chime bell of Marquis Yi of Zeng. The preliminary calculations have clearly shown that two different modes of oscillation can be generated by striking different parts of the bell, and indicate that finite element numerical simulation of the processes of wave propagation and two-tone generation of the chime bell of Marquis Yi of Zeng is feasible. These analyses provide a new quantitative and visual way to explain the mystery of the two-tone acoustic characteristics. The method suggested by this study can be applied to simulate free oscillations excited by great earthquakes with complex Earth structure. Taking into account of such large-scale structure of the Earth, small-scale low-precision numerical simulation can not simply meet the requirement. The increasing capacity in high-performance parallel computing and progress on fully numerical solutions for seismic wave fields in realistic three-dimensional spherical models, Spectral element method and high-performance parallel computing were incorporated to simulate the seismic wave propagation processes in the Earth’s interior, without the effects of the Earth’s gravitational potential. The numerical simulation shows that, the results of the toroidal modes of our calculation agree well with the theoretical values, although the accuracy of our results is much limited, the calculated peaks are little distorted due to three-dimensional effects. There exist much great differences between our calculated values of spheroidal modes and theoretical values, because we don’t consider the effect the Earth’ gravitation in numerical model, which leads our values are smaller than the theoretical values. When , is much smaller, the effect of the Earth’s gravitation make the periods of spheroidal modes become shorter. However, we now can not consider effects of the Earth’s gravitational potential into the numerical model to simulate the spheroidal oscillations, but those results still demonstrate that, the numerical simulation of the Earth’s free oscillation is very feasible. We make the numerical simulation on processes of the Earth’s free oscillations under spherically symmetric Earth model using different special source mechanisms. The results quantitatively show that Earth’s free oscillations excited by different earthquakes are different, and oscillations at different locations are different for free oscillation excited by the same earthquake. We also explore how the Earth’s medium attenuation will take effects on the Earth’s free oscillations, and take comparisons with the observations. The medium attenuation can make influences on the Earth’s free oscillations, though the effects on lower-frequency fundamental oscillations are weak. At last, taking 2008 Wenchuan earthquake for example, we employ spectral element method incorporated with large-scale parallel computing technology to investigate the characteristics of seismic wave propagation excited by Wenchuan earthquake. We calculate synthetic seismograms with one-point source model and three-point source model respectively. Full 3-D visualization of the numerical results displays the profile of the seismic wave propagation with respect to time. The three-point source, which was proposed by the latest investigations through field observation and reverse estimation, can better demonstrate the spatial and temporal characteristics of the source rupture processes than one-point source. Primary results show that those synthetic signals calculated from three-point source agree well with the observations. This can further reveal that the source rupturing process of Wenchuan earthquake is a multi-rupture process, which is composed by at least three or more stages of rupture processes. In conclusion, the numerical simulation can not only solve some problems concluding the Earth’s ellipticity and anisotropy, which can be easily solved by conventional methods, but also finally solve the problems concluding topography model and lateral heterogeneity. We will try to find a way to fully implement self-gravitation in spectral element method in future, and do our best to continue researching the Earth’s free oscillations using the numerical simulations to see how the Earth’ lateral heterogeneous will affect the Earth’s free oscillations. These will make it possible to bring modal spectral data increasingly to bear on furthering our understanding of the Earth’s three-dimensional structure.

南海北部大陆边缘盆地深部地热特征用地球动力学演化

This report is a conclusion of the major research outcome during my post-doctoral residence of research and work. Its content covers the researches of the deep thermal characteristic and dynamics evolution beneath the northern margin basin of South China Sea. In this report, the each other action and effect between lithosphere ad mantle convection were regarded by the combine of deep and shallow study, subdivision from whole to part, and pay equal attention to determine the nature and fixed quantity. The investigative method we used in this report is geothermal and gravity methods. By the help of geological model and geophysics modeling, we calculated lithosphere thermal structure, rheology structure and mantle convection. Firstly, the report introduces concisely the purpose and the previous achievement to this research. Then, it analyzed the characteristic of heat flow on South China Sea. The structure of deep temperature and thermal has been calculated in some models of heat generation and conduction. The rock rheology structure also was computed by the relationship between temperature and viscosity. All these calculations were finished under the guidelines of combine with geology and geophysics. Meanwhile, the fields both deep mantle convection and small scale upper mantle convection are computed. Beside, the density and temperature disorder resulted by mantle convection were also computed with the convection field. After these, the report bring the contribution of local field of mantle convection, thermal construct and effective viscosity beneath the northern margin basin of South China Sea. And, base on the tectonic background and evolution feature, this report discussion the evolution mechanism of south China Sea and its northern margin basin. The end of this report, the main conclusion of this research was summarized and brings out.

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

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.

复杂断块油藏形成机制及剩余油预测

The central uplift in the Huimin depression is famous for its large amounts of faults and small-scale fault-block area, and it is the famed typical complicated fault-block group oil & gas field in the whole world. After many years of rolling exploration and exploitation, many complex oil &gas field have been discovered in the central uplift, and won the splendent fruit. With the gradual deepening and development of the rolling exploitation, the exploration faces more and more difficulties. Therefore, it is important to reveal the forming mechanism and distributing rule of the complex fault-block reservoir, and to realize the forecast of the complex fault-block reservoir, sequentially, expedite the exploration step. This article applies the new multi-subject theory, method and technique such as structure geometry, kinematics, dynamics, structural stress field, fluid potential field, well logging record and constrained inversion of seismic records, coherence analysis, the seal mold and seal history of oil-bounded fault etc, and try to reveal the forming mechanism and distributing law of the complex fault-block reservoir, in result, implements the forecast of the fault-block reservoir and the remaining oil distributing. In order to do so, this article synthetically carries out structural estimate, reservoir estimate, fault sealing history estimate, oil-bearing properties estimate and residual. This article also synthetically researches, describes and forecast the complex fault-block in Huimin depression by use of the techniques, e.g. seismetic data post-stack processing technique, multi-component demarcating technique, elaborate description technique for the fault-block structure, technique of layer forecasting, fault sealing analysis technique, comprehensive estimate technique of fault-block, comprehensive analysis and estimate technique of remaining oil etc. The activities of the faults varies dramatically in the Huimin depression, and most of the second-class and the third-class faults are contemporaneous faults, which control the macroscopical distribution of the reservoir in the Huimin depression. The fourth-class faults cause the complication between the oil & gas among the fault-blocks. The multi-period strong activities of the Linyi fracture resulted in the vertical migration of large amount of oil & gas along with the faults. This is the main reason for the long vertical distribution properties near the Linyi fracture in the Huimin depression. The sealing ability of the fault is controlled by the property,size and direction of the main stress, the contact relationship of the both sides of the fault, the shale polluting factor, and the configuration relationship between the fault move period and the migration period of oil & gas. The article suggest four fault-sealing modes in the research zone for the first time, which establishes the foundation for the further forecast of the complex fault-block reservoir. Numerical simulation of the structural stress field reveals the distribution law and the evolvement progress of the three-period stress field from the end of the Dongying period to the Guantao period to nowadays. This article puts forward that the Linyi and Shanghe regions are the low value of the maximum main stress data. This is combined with the fault sealing history estimate, then multi-forming-reservoir in the central uplift is put forward. In the Shanghe oilfield, the article establishes six reservoir geological modes and three remaining oil distributing modes(the plane, the inside layer and the interlayer), then puts forward six increase production measure to enhance the remaining oil recovery ratio. Inducting the exploitation of oilfield, it wins notable economic effects and social effects.