准噶尔盆地滴南凸起中生界层序地层及有利勘探目标优选

Based on the theories of sequence stratigraphy and sedimentology, as well as comprehensive studies of seismic data, drilling data, core interpretation and setting of this area, the thesis presents an analysis for Mesozoic formation in Dinan uplift. By means of recognizing the boundary of the sequence, dividing and correlating the systems tract, Mesozoic of Dinan uplift is divided into ten sequences and twenty-five systems tracts during the establishment of the sequence framework. In the framework, some sequences are featured by mature systems of lowstand, water-transgression and highstand, while some undeveloped systems of lowstand or highstand. The main sedimentary facies in Mesozoic of Dinan uplift are braided river, meandering river, delta and lake. The braided river was divided into sandy river and rudaceous river by the lithology of the river channel and was divided into dry climate and wet climate condition by the color of the flood plain. Additionally, The concept of “wetland” is put forward for the first time and regarded as the consequence of wet climate. The analysis includes the classification of six types of traps: (1) stratigraphic overlap trap, (2) lithologic trap with updip pinchout, (3) stratigraphic unconformity trap, (4) fault-lithology trap, (5) fault trap, (6) anticlinal trap, and combining with the research of the characteristics and distribution rules for the known reservoir, it draws out that “fault control” is the petroleum accumulation pattern in this area, in which fault is the key element of the transporting system. Finally the thesis concludes the distribution characteristics and optimized some targets for the potential exploration zone.

松辽盆地长岭断陷火山岩气藏的形成条件与分布规律

Changling fault depression is a compound fault depression complicated by interior fault, with faults in the west and overlap in the west. North of Changling fault depression show NNE strike while south is NW strike. Changling fault depression has undergone twochasmic stage which control the development and distribution of volcanic rock, one depression stage, later inversion and uplift stage which control the formation of natural gas reservoir, and basin atrophic stage. The main boundary faults and main faults in Changling fault depression control three volcanic cycles and the distribution of volcanic rock. Seismic reflection characteristic and logging response characteristic of volcanic rock in study area are obvious, and the distribution characteristic, volcanic cycle and active stage of volcanic rock can be revealed by seismic attribute, conventional logging data can distinguish clastic rock from volcanic rock or distinguish partial different types of volcanic rock. The reservoir property of rhyolite and volcanic tuff are the best. Favorable volcanic reservoir can be preserved in deep zone. Imaging logging and frequency decompostion technology of seismic data act as effective role in the study of reservoir physical property and gas-bearing properties of volcanic rock.. Hydrocarbon gas in study area is high and over mature coal type gas, the origin of CO2 is complex, it is either inorganic origin or organic origin, or mixing origin. Hydrocarbon gas is mainly originate from Shahezi formation and Yingcheng formation source rocks, CO2 is mainly mantle source gas. Hydrocarbon has the characteristics of continuous accumulation with two charging peak. The first peak represent liquid hydrocarbon accumulation time, The second peak stand for the accumulation time of gaseous hydrocarbon.CO2 accumulate approximately in Neocene. The source rock distribution range, volcanic rock and favorable reservoir facies, distribution characteristic of deep fault (gas source fault) and late inversion structure are the major factors to control gas reservoir formation and distribution. All the results show that these traps that consist of big inherited paleo uplift(paleo slope), stratigraphic overlap and thinning out, volcanic rock, are the most advantageous target zone.

渤海湾盆地前第三系构造演化与煤成气成藏研究

The disequilibrium between supply and demand the east part of North China accelerated natural gas exploration in Bohai bay basin. Exploration practice showed that coal-derived gas is important resource. In searching of big to middle scaled coal derived gas field, and realize successive gas supply, the paper carried out integrated study on structural evolution of Pre-Tertiary and evaluation of reservoir forming condition of coal-derived gas. Study work of the paper was based on the following condition: available achievement in this field at present, good understanding of multiphase of tectonic movement. Study work was focused on geological evolution, source rock evaluation and dissection key factors controlling reservoir forming. Based on analysis of seismic data, drilling data, tectonic style of Pre-Tertiary was subdivided, with different tectonic style representing different tectonic process. By means of state of the art, such as analysis of balanced cross section, and erosion restoration, the paper reestablished tectonic history and analyzed basin property during different tectonic phase. Dynamic mechanism for tectonic movement and influence of tectonic evolution on tectonic style were discussed. Study made it clear that tectonic movement is intensive since Mesozoic including 2 phase of compressional movement (at the end of Indo-China movement, and Yanshan movement), 2 phase of extensional movement (middle Yanshan movement, and Himalayan movement), 2 phase of strike slip movement, as well as 2 phase of reversal movement (early Yanshan movement, and early Himalayan movement). As a result, three tectonic provinces with different remnant of strata and different tectonic style took shape. Based on afore mentioned study, the paper pointed out that evolution of Bohai bay basin experienced the following steps: basin of rift valley type (Pt2+3)-craton basin at passive continental margin (∈1－2)-craton basin at active continental margin (∈3- O)-convergent craton basin (C-T1+2)-intracontinental basin (J＋K). Superposition of basins in different stage was discussed. Aimed at tectonic feature of multiple phases, the paper put forward concept model of superposition of tectonic unit, and analyzed its significance on reservoir forming. On basis of the difference among 3 tectonic movements in Mesozoic and Cenozoic, superposition of tectonic unit was classified into the following 3 categories and 6 types: continuous subsidence type (I), subsidence in Mesozoic and uplift for erosion in Cenozoic (II1), repeated subsidence and uplift in Mesozoic and subsidence in Cenozoic (II2), repeated subsidence and uplift in Mesozoic and uplift for erosion in Cenozoic (II3), uplift for erosion in Mesozoic and subsidence in Cenozoic (II4), and continuous uplift (III). Take the organic facies analysis as link, the paper established relationship between sedimentary environment and organic facies, as well as organic facies and organic matter abundance. Combined information of sedimentary environment and logging data, the paper estimated distribution of organic matter abundance. Combined with simulation of secondary hydrocarbon generation, dynamic mechanism of hydrocarbon generation, and thermal history, the paper made static and dynamic evaluation of effective source rock, i.e. Taiyuan formation and Shanxi formation. It is also pointed out that superposition of tectonic unit of type II2, type II4, and type I were the most favorable hydrocarbon generation units. Based on dissection of typical primary coal-derived gas reservoir, including reservoir forming condition and reservoir forming process, the paper pointed out key factors controlling reservoir forming for Carboniferous and Permian System: a. remnant thickness and source rock property were precondition; b. secondary hydrocarbon generation during Himalayan period was key factor; c. tectonic evolution history controlling thermal evolution of source rock was main factor that determine reservoir forming; d. inherited positive structural unit was favorable accumulation direction; e. fault activity and regional caprock determined hydrocarbon accumulation horizon. In the end, the paper established reservoir forming model for different superposition of tectonic units, and pointed out promising exploration belts with 11 of the first class, 5 of the second class and 6 of the third class.

地幔压力下CaSiO3 Perovskite状态方程和弹性常数的第一性原理研究

In the past decade density functional theory (DFT) has made its way from a peripheral position in quantum chemistry to center. Of course the often excellent accuracy of the DFT based methods has provided the primary driving force of this development. This dissertation is devoted to the study of physical and chemical properties of planetary materials by first-principle calculation. The concerned properties include the geometry, elastic constants and anisotropy. In the first chapter, we give a systematic introduction to theoretical background and review its progress. Development of quantum chemistry promotes the establishment of DFT. Theorem of Hohenberg-Kohn is the fundament of DFT and is developed to Kohn-Sham equation, which can be used to perform real calculations. Now, new corrections and extensions, together with developed exchange-correlation, have made DFT more accurate and suitable for larger systems. In the second chapter, we focus on the calculational methods and technical aspects of DFT. Although it is important to develop methods and program, external package are still often used. At the end of this chapter, we briefly some widely used simulation package and the application of DFT. In the third chapter, we begin to focus on properties of real materials by first principles calculation. We study a kind of minerals named Ca perovskite, investigate its possible structure and anisotropy at Earth’s mental condition. By understanding and predicting geo-physically important materials properties at extreme conditions, we can get the most accurate information to interpret seismic data in the context of likely geophysical processes.

积分偏移的计算几何与并行实现

The Second Round of Oil & Gas Exploration needs more precision imaging method, velocity vs. depth model and geometry description on Complicated Geological Mass. Prestack time migration on inhomogeneous media was the technical basic of velocity analysis, prestack time migration on Rugged surface, angle gather and multi-domain noise suppression. In order to realize this technique, several critical technical problems need to be solved, such as parallel computation, velocity algorithm on ununiform grid and visualization. The key problem is organic combination theories of migration and computational geometry. Based on technical problems of 3-D prestack time migration existing in inhomogeneous media and requirements from nonuniform grid, parallel process and visualization, the thesis was studied systematically on three aspects: Infrastructure of velocity varies laterally Green function traveltime computation on ununiform grid, parallel computational of kirchhoff integral migration and 3D visualization, by combining integral migration theory and Computational Geometry. The results will provide powerful technical support to the implement of prestack time migration and convenient compute infrastructure of wave number domain simulation in inhomogeneous media. The main results were obtained as follows: 1. Symbol of one way wave Lie algebra integral, phase and green function traveltime expressions were analyzed, and simple 2-D expression of Lie algebra integral symbol phase and green function traveltime in time domain were given in inhomogeneous media by using pseudo-differential operators’ exponential map and Lie group algorithm preserving geometry structure. Infrastructure calculation of five parts, including derivative, commutating operator, Lie algebra root tree, exponential map root tree and traveltime coefficients , was brought forward when calculating asymmetry traveltime equation containing lateral differential in 3-D by this method. 2. By studying the infrastructure calculation of asymmetry traveltime in 3-D based on lateral velocity differential and combining computational geometry, a method to build velocity library and interpolate on velocity library using triangulate was obtained, which fit traveltime calculate requirements of parallel time migration and velocity estimate. 3. Combining velocity library triangulate and computational geometry, a structure which was convenient to calculate differential in horizontal, commutating operator and integral in vertical was built. Furthermore, recursive algorithm, for calculating architecture on lie algebra integral and exponential map root tree (Magnus in Math), was build and asymmetry traveltime based on lateral differential algorithm was also realized. 4. Based on graph theory and computational geometry, a minimum cycle method to decompose area into polygon blocks, which can be used as topological representation of migration result was proposed, which provided a practical method to block representation and research to migration interpretation results. 5. Based on MPI library, a process of bringing parallel migration algorithm at arbitrary sequence traces into practical was realized by using asymmetry traveltime based on lateral differential calculation and Kirchhoff integral method. 6. Visualization of geological data and seismic data were studied by the tools of OpenGL and Open Inventor, based on computational geometry theory, and a 3D visualize system on seismic imaging data was designed.

北部湾盆地迈陈凹陷地震勘探方法研究

Maichen Depression lie between Leizhou Peninsula and Qiongzhou Strait. Oil and gas have been discovered in Weixinan Depression, Wushi Depression and Fushan Depression, which pertain to a same basin — North Sea Basin along with Maichen Depression.Jiangsu Oil started exploration at 2002. The first well began to drill at November, 2004 after gravity survey, electric method prospecting and 2D seismic exploration had been finished. Generating rock and hydrocarbon shows have been verified by the drilling. Low yield oil stream has been tested. And we started 3D seismic exploration at November, 2005. My thesis topic came from the actual needs of our exploration in the Maichen Depression. In the thesis, I give emphasis to analyse the own seismic geologic conditions of Maichen Depression. By real tests, we choosed the means to overcome or weaken the unfavorably impress owing to the own coditions in Maichen Depression. Finally, we obtained the usable seismic data. 1. Owing to the multiphase eruptive rock during the Quaternary Period, the near surface layers are very inhomogeneous. By simultaneous testing at same point with short refraction, uphole surveys of radial source and of surface source, the most appropriate method had been sorted out. Radial source uphole survey has been regarding the best practicable means in the complex area. Accurate surficial geology was very helpful to choosing of acquirement means and parameters. Basically the appropriate method of seismic acquirement has been built at Maichen area. 2. The seismic primary data has many, very strong and complex noise. By noise characteristic analysis in different domain, many means of denoising had been paralleled individual and joint application researched. As a result, the pre-stack multidomain joint denoise flow was the appropriate method. It can improve the seismic signal-to-noise ratio. 3. The problem of seismic static correction at Maichen Depression is very conspicuous. Many static correction methods had been tested individual and joint researched. The seismic data quality has been improved after choosing the appropriate combination of static correction flows. 4. Although the above-mentioned process are resultful, the seismic profile quality is just passable. Some reflector continuity and fault zone imagery are ambiguity. So it was the useful method to reduce the structural ambiguity during seismic interpretation that built-up geologic model in accord with real geologic character by areal structure study upon backbone seismic profiles. In the same way, traps have been assessed and drill targets have been selected.

砂泥岩薄互层储层预测方法与应用研究

With the deeply development of exploration and development in petroleum in China, new increasing reserves are found in old oil fields and the verge of the old ones through re-study of geological property. It is more and more important to discovery and develop thin layer or thin inter-bedded layers reservoirs. All of the targets are thin sand-shale inter-bedded reservoirs and the core technology is reservoir predictions between wells in thin sand-shale inter-bedded layers. The continuity of the thin sand-shale inter-bedded layers in space or separating and heterogeneity is the key of reservoir geology research. The seismic reflection, high resolution analysis method and inversion method to thin sand-shale inter-bedded layers are thorough discussed and deeply studied in this paper to try to find the methods and resolutions of reservoir geology research. The below is followed. 1. Based on the pre-research of other people, five models are created: the sand sphenoid body, interlay sandstone and interlay shale of the equal thickness, interlay sandstone of the equal thickness and interlay shale of the unequal thickness, interlay sandstone of the unequal thickness and interlay shale of the unequal thickness, interlay sandstone of the changing thickness in sequence and interlay shale of the changing thickness in sequence. Then the study of the forward modeling are conducted on the thin layer and thin inter-bedded layers geological characters and seismic reflections including amplitude, frequency, phase, wave shape and time-frequency responding in the domains of time and frequency. The affect of petro-physics difference of layers, single thin layer thickness, thickness of inter-bedded, layer number of inter-bedded, incident wavelet domain frequency and types, sample interval to seismic reflection characters, frequency spectrum and time-frequency respond of reflectivity is theoretically discussed. 2. Qualitatively analyzing the sedimentary rhythm of the thin inter-bedded layers in vertical orientation and computing the single layer thickness or the average thickness with the method of generalized S transform. Identifying the reflecting interface or lithology interface using the amplitude value of amplitude spectrum domain frequency. 3. Based on the seismic respond of thin sand-shale inter-bedded layers, bring out the high resolution analysis method of seismic data in thin sand-shale inter-bedded layers using wavelet analysis and the idea of affecting low and high frequency with middle frequency. Then analyzing the effect to the method and testing some wavelets in the method. This method is applied to the theoretical models and the field data. 4. Bring forward one improved very fast simulated annealing method (IVFSA) to resolve the problem nonlinearity and multi-parameters of the inversion in thin inter-bedded layers. And IVFSA is more productive and higher precision than general ways. 5. New target constrained function is used in the inversion based on the property of the inversion in thin inter-bedded layers. 6. Making the full use of geological and logging information, IVFSA and the new function are applied in the non-linear inversion to improve reservoir prediction and evaluation in thin inter-bedded formations combined with the idea of logging and seismic inversion. This method was applied to the field data and got good results.

高分辨率非线性岩性物性联合反演方法和应用

In the prediction of complex reservoir with high heterogeneities in lithologic and petrophysical properties, because of inexact data (e.g., information-overlapping, information-incomplete, and noise-contaminated) and ambiguous physical relationship, inversion results suffer from non-uniqueness, instability and uncertainty. Thus, the reservoir prediction technologies based on the linear assumptions are unsuited for these complex areas. Based on the limitations of conventional technologies, the thesis conducts a series of researches on various kernel problems such as inversions from band-limited seismic data, inversion resolution, inversion stability, and ambiguous physical relationship. The thesis combines deterministic, statistical and nonlinear theories of geophysics, and integrates geological information, rock physics, well data and seismic data to predict lithologic and petrophysical parameters. The joint inversion technology is suited for the areas with complex depositional environment and complex rock-physical relationship. Combining nonlinear multistage Robinson seismic convolution model with unconventional Caianiello neural network, the thesis implements the unification of the deterministic and statistical inversion. Through Robinson seismic convolution model and nonlinear self-affine transform, the deterministic inversion is implemented by establishing a deterministic relationship between seismic impedance and seismic responses. So, this can ensure inversion reliability. Furthermore, through multistage seismic wavelet (MSW)/seismic inverse wavelet (MSIW) and Caianiello neural network, the statistical inversion is implemented by establishing a statistical relationship between seismic impedance and seismic responses. Thus, this can ensure the anti-noise ability. In this thesis, direct and indirect inversion modes are alternately used to estimate and revise the impedance value. Direct inversion result is used as the initial value of indirect inversion and finally high-resolution impedance profile is achieved by indirect inversion. This largely enhances inversion precision. In the thesis, a nonlinear rock physics convolution model is adopted to establish a relationship between impedance and porosity/clay-content. Through multistage decomposition and bidirectional edge wavelet detection, it can depict more complex rock physical relationship. Moreover, it uses the Caianiello neural network to implement the combination of deterministic inversion, statistical inversion and nonlinear theory. Last, by combined applications of direct inversion based on vertical edge detection wavelet and indirect inversion based on lateral edge detection wavelet, it implements the integrative application of geological information, well data and seismic impedance for estimation of high-resolution petrophysical parameters (porosity/clay-content). These inversion results can be used to reservoir prediction and characterization. Multi-well constrains and separate-frequency inversion modes are adopted in the thesis. The analyses of these sections of lithologic and petrophysical properties show that the low-frequency sections reflect the macro structure of the strata, while the middle/high-frequency sections reflect the detailed structure of the strata. Therefore, the high-resolution sections can be used to recognize the boundary of sand body and to predict the hydrocarbon zones.

奇异性分析方法及其在储层预测中的应用研究

In research field of oil geophysical prospecting, reservoir prediction is refers to forecasting physical properties of petroleum reservoir by using data of seismic and well logging, it is a research which can guide oil field development. Singularities of seismic and logging data are caused by the heterogeneity of reservoir physical property. It's one of important methods that using singularity characteristics of seismic and logging data to study the reservoir physical property in recently. Among them, realization of reservoir quantitative prediction by analyzing singularity of the data and enhancing transition description of data is difficulty in method research. Based on wavelet transform and the fractal theory, the paper studied the singularity judgment criterion for seismic and logging data, not only analyzed quantitative relation between singularity data and reservoir physical property, but also applied it in practical reservoir prediction. The main achievements are: 1. A new method which provides singular points and their strength information estimation at only one single scale is proposed by Herrmann (1999). Based on that, the dissertation proposed modified algorithm which realized singularity polarity detection. 2. The dissertation introduced onset function to generalize the traditional geologic boundaries variations model which used singularity characteristics to represent the abruptness of the lithologic velocity transition. We show that singularity analysis reveals generic singularity information conducted from velocity or acoustic impedance to seismogram based on the convolution seismic-model theory. Theory and applications indicated that singularity information calculated from seismic data was a natural attribute for delineating stratigraphy boundaries due to its excellent ability in detecting detailed geologic features. We demonstrated that singularity analysis was a powerful tool to delineate stratigraphy boundaries and inverse acoustic impedance and velocity. 3. The geologic significances of logging data singularity information were also presented. According to our analysis, the positions of singularities indicate the sequence stratigraphic boundary, and there is subtle relationship between the singularity strength and sedimentary environment, meanwhile the singularity polarity used to recognize stratigraphic base-level cycle. Based on all those above, a new method which provided sedimentary cycle analysis based on the singularity information of logging data in multiple scales was proposed in this dissertation. This method provided a quantitative tool for judging interface of stratum sequence and achieved good results in the actual application.

岩石中地震波速度和衰减的应力变化响应研究

Stress change is one of key factors in seismic nucleating and triggering; therefore for understanding and forecasting earthquakes, it is necessary to research on stress status and its changes in rocks. Propagating in underground structures, wave velocity and attenuation contain information on stress changes of the Earth’s interior. For a better understanding of relationship between seismic data and stress changes, modeling and ultrasonic test supply significant references. In this article, acoustoelastic theory is introduced to explain nonlinear elastic characteristics of rocks. Based on the acoustoelastic theory, a solid-fluid coupled model is given to calculate velocity under different stress for porous and liquid fulfilled rocks. Except for the stress-velocity relationship, effects of pore pressure induced stress changes on ultrasonic coda attenuation are also studied. Intrinsic attenuation quality factors are calculated for a comparison purpose. Finally, the relationship between elastic constants and stress changes is thoroughly investigated, a mixture model from two phases of Hooke media is introduced to explain the differences between dynamic and static moduli, a relation among wave length, wave velocities and elastic moduli considering dimension of microstructure, dimension and state of surface between phases is presented. The most important aspect of this work is exploring and establishing relationships between the seismic properties of rocks and changes of their stress conditions, which will have its application in earthquake forecast and seismic hazard.

煤系地层三维三分量地震数据解释方法研究

Conventional 3D seismic exploration cannot meet the demand of high yield and high efficiency safe production in coal mine any more. Now it is urgent to improve the discovery degree of coal mine geological structures for coal production in China. Based on 3D3C seismic exploration data, multi-component seismic information is fully excavated. First systematic research on 3D3C seismic data interpretation of coal measure strata is carried out. Firstly, by analyzing the coal measure strata, the seismic-geologic model of coal measure strata is built. Shear wave logging is built by using regression analysis. Horizon calibration methods of PP-wave and PS-wave are studied and the multi-wave data are used together to interpret small faults. Using main amplitude analysis technology, small faults which cannot be found from PP-wave sections can be interpreted from the low frequency PS-wave sections. Thus, the purpose to applying PS-wave data to fine structure assistant interpretation is achieved. Secondly, PP- and PS-wave post-stack well constrained inversion methods of coal measure strata are studied. Joint PP- and PS-wave post-stack inversion flow is established. More attribute parameters, which are applied in fine lithology interpretation of coal measure strata, are obtained from combinations of the inversion results. Exploring the relation between rock with negative Poisson’s ratio and anisotropy, fracture development in coal seam are predicted. Petrophysical features of coal measure strata are studied, and the relations between elastic parameters and lithology, fluid and physical properties are established. Inversions of the physical parameters such as porosity, permeability and water saturation, which reflect lithology and fluid property, are obtained. Finally, the approaches of shear wave splitting and Thomsen parameters inversion, which provide new ideas for seismic anisotropy interpretation of coal measure strata, are studied to predict fracture development. The results of practical application indicate that the methods in this paper have good feasibility and applicability. They have positive significance for high yield and high efficiency safe production in coal mine.

南海北部陆缘深水区构造－热演化与烃源岩成熟度史研究

Abstract In order to provide basic data for evaluation of the petroleum potential in the deep water area of the northern margin of the South China Sea (SCS), present-day thermal regime and basin tectonothermal evolution are reconstructed and the maturation history of the Cenozoic major source rocks in the study area is derived. The present-day geothermal regime in the deep water area of the northern margin of SCS is defined according to the geothermal gradient, thermal properties and heat flow data. Tectonic subsidence history is reconstructed based on borehole and seismic data, and accordingly the stretching episodes are determined from the subsidence pattern. Heat flow history in the deep water area of the northern margin of SCS is estimated on a finite time, laterally non-uniform and multi-episode stretching model. Maturation history of the main source rocks in the study area is estimated through EASYRo% kinetic model and thermal history, and the potential of petroleum in the deep water area of the northern margin of SCS is evaluated based on the data above. The results show that the present-day geothermal regime in the deep water area of the northern margin of SCS is characterized by “hot basin” with high geothermal gradient (39.1±7.4℃/km) and high heat flow (77.5±14.8 mW/m2), and that the Qiongdongnan Basin (QDNB) underwent three stretching episodes and consequently suffered three heating episodes (Eocene, Oligocene and Pliocene time) with highest paleo-heat flow of 65～90 mW/m2 at the end of the Pliocene, that the Pearl River Mouth Basin (PRMB) two stretching and two heating episodes (Eocene, Oligocene time) with highest paleo-heat flow of 60～70 mW/m2 at the end of the Oligocene, and that the source rocks matured drastically responding to the heating episodes. There are four hydrocarbon generation kitchens in the deep water area of the northern margin of SCS which are favor of its bright petroleum perspective. Tectonothermal analysis indicates that the present-day geothermal regime which is characterized with “hot basin” in the deep water area of the PRMB resulted mainly from the Cenozoic stretching as well as faulting and magmatic activities during the Neotectonic period, and that the Pliocene heating episode of the QDNB is coupled with the transition from sinistral to dextral gliding of the Red Rive fault, and that the deep water basins in the northern margin of SCS are typical of multiple rifting which caused multi-episode heating process.

微动勘查方法在煤矿采区勘探中的应用研究

The conventional microtremor survey is based on the single point of exploration, which includes collecting field data，estimating the phase velocity, investing the dispersion curve and obtaining the S–wave velocity structure. In the case of large-scale exploration, and when making the two-dimensional velocity section, the inversion is quite time-consuming, laborious and its precision depends on the subjective interpretation, which makes the results differently from person to person. In fact，we do not need the S-wave velocity values but only need the relative variation of velocity. For these reasons, this paper is desired to calculate the apparent S-wave velocity (Vx) to replace the S-wave velocity inversion and to obtain the relative variation of the S-wave velocity. Using this method, we can decrease the analysts’ effect, shorten the data processing time and improve work efficiency. The apparent S-wave velocity is a variable of the surface wave property, which can clearly reflect the downcast columns, mined-out areas and other unusual geological bodies. In this paper, Matlab is used to establish the three-dimensional data volume of the apparent S-wave velocity, from which we can get any apparent S-wave velocity section we need. Through the application case, the designed method is proved to be reliable and effective. The downcast columns, mined-out areas and other unusual geological bodies can be clearly showed in the apparent S-wave velocity section. And from the contour of the apparent S-wave velocity, the interface shape of the major target layers can be controlled basically.

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

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.

一种TEC地图构建方法与全球电子总含量的气候学研究

Global positioning system (GPS) can not only provide precise service for navigation and timing, but also be used to investigate the ionospheric variation. From the GPS observations, we can obtain total electron content (TEC), so-called GPS TEC, which is used to characterize the ionospheric structure. This thesis mainly concerns about GPS TEC data processing and ionospheric climatological analysis as follows. Firstly, develop an algorithm for high-resolution global ionospheric TEC mapping. According to current algorithms in global TEC mapping, we propose a practical way to calibrate the original GPS TEC with the existing GIM results. We also finish global/local TEC mapping by model fitting with the processed GPS TEC data; in practice, we apply it into the local TEC mapping in Southeast of China and obtain some initial results. Next, suggest a new method to calculate equivalent ionospheric global electron content (GEC). We calculate such an equivalent GEC with the TEC data along the geographic longitude 120°E. With the climatological analysis, we can see that GEC climatological variation is mainly composed of three factors: solar cycle, annual and semiannual variations. Solar cycle variation is dominant among them, which indicates the most prominent influence; both annual and semiannual variations play a secondary role and are modulated by solar activity. We construct an empirical GEC model driven by solar activity and seasonal factors on the basis of partial correlation analysis. Generally speaking, our researches not only show that GPS is advantageous in now-casting ionospheric TEC as an important observation, but also show that GEC may become a new index to describe the solar influence on the global ionosphere since the great correlation between GEC and solar activity factor indicates the close relationship between the ionosphere and solar activity.