25ka B.P.以来黑潮流域古环境演化对高频全球变化事件的响应

本文利用黑潮流域主流轴上的两个柱状沉积物岩芯MD05-2908以及PC-1为研究材料，在AMS14C测年的基础上，利用高分辨率的有机地球化学分析记录结合浮游有孔虫氧、碳同位素，恢复并重建了过去25,000 cal a BP以来黑潮流域古海洋环境演化的历史。以及表层海水生产力以及物质输送状况的演化历史。通过利用Uk’37古海水温度以及盐度指标恢复重建了过去25,000 cal a BP以来海洋表层海水温度、盐度；通过机碳、氮含量以及同位素变化、长链正构烷烃以及正构烷醇等指标重建了过去7000a B.P.以来的陆源物质输入历史；通过长链不饱和烯酮含量以及有机碳同位素指标恢复了过去7000a B.P.以来海水表层生产力的历史。此外，通过基于以上各种指标的环境信息与区域以及全球其它气候记录进行对比研究以及不同环境指标的时间系列分析，探讨了该区表层环流系统以及生产力的演化对于全球气候变化的响应，揭示了不同尺度的短周期高频率全球变化事件在黑潮流域的具体作用过程和响应机制。通过这些研究取得了以下的主要认识： 基于有机地球化学指标的古气候环境记录与黑潮流域已有的研究成果有很好的对应关系，从我们高分辨率的有机地球化学记录中可以识别出全新世黑潮强弱变化的几次明显的事件；25000a B.P.以来黑潮流域的环境变化与全球环境变化有着很好的对应性，黑潮强弱演化总体趋势与全球气候背景演化相一致，黑潮对高频气候变化事件的记录与全球记录具有同步性，这种同步性尤其体现在末次冰消期以来的气候快速高频振荡以及全新世以来的气候突变事件上。 全球性的高频气候事件对黑潮主体本身及黑潮流域的相邻区域的大气和海洋环流都具有重要的控制作用，这种控制作用主要通过副高、ITCZ以及季风三种气候要素之间相互关联、彼此影响造成的。具体表现为：太阳辐射量的减少导致热力差异减小，这种相对减小弱化了热带西太平洋的对流活动，造成了西太平洋副热带高压长期偏南、偏东，ITCZ平均北界位置偏南，降雨带长时间集中在南部地区，增强的降雨量提高了风化剥蚀以及沉积物向海洋搬运的能力，陆源物质供应量增加；同时，辐射量以及热力差的减小又与加强的东亚冬季风相联系，增强的冬季风导致了近底层“雾状层”物质的向海传输，物质传输效率增高。这种物源供应以及搬运量的双重增加导致了冲绳海槽流域物质通量的增加。 基于有机地球化学指标的海洋表层生产力的变化与陆源物质供应量以及黑潮流的强弱变化存在着对应关系，通常情况增加的海洋表层生产力对应着高的陆源物质输入以及相对较弱的黑潮。这种变化与东亚夏季风的以及冬季风的强弱都有很好的一致性。陆源物质的输入增加了表层营养物质的含量，导致生产力的勃发；陆源物质的输入增加又对应着减少的太阳辐射量，偏南的ITCZ北界位置以及副热带高压，这些对应于减弱的东亚夏季风（增强的东亚冬季风）。 黑潮流域的高沉积速率事件对应于减弱的黑潮强度和增加的ENSO频度，这些事件与上述的副热带高压、ITCZ位移和强弱的变化相一致。黑潮流域过去25000a B.P.以来南北温度的差异有冰期加大而全新世减小的趋势，但这种冰期与全新世的差异很小，我们认为末次盛冰期的时候黑潮主流轴没有移出冲绳黑潮，只是由于强度的减弱受陆架水体的影响有所加大。 黑潮流域很好的记录到了包括数千年尺度的D/O旋回周期到大气——海洋系统内部振荡所致的PDO、NAO等数十年尺度的高频振荡，说明黑潮流域对过去全球及区域环境变化事件有很好的响应。黑潮流域各古海洋指标所记录周期上的一致性说明这些环境因子控制机理上具有的一致性，即大背景上受太阳活动所引起的辐射量变化控制，局部的高频快速气候波动又受到局域性的气候因素如海气相互作用的放大影响。

Interannual horizontal heat advection in the surface mixed layer over the equatorial Pacific Ocean: assimilation versus TAO analysis

The interannual anomalies of horizontal heat advection in the surface mixed layer over the equatorial Pacific Ocean in an assimilation experiment are studied and compared with existing observational analyses. The assimilation builds upon a hindcast study that has produced a good simulation of the observed equatorial currents and optimizes the simulation of the Reynolds sea surface temperature (SST) data. The comparison suggests that the assimilation has improved the simulation of the interannual horizontal heat advection of the surface mixed layer significantly. During periods of interrupted current measurements, the assimilation is shown to produce more meaningful anomalies of the heat advection than the interpolation of the observational data does. The assimilation also shows that the eddy heat flux due to the correlation between high-frequency current and SST variations, which is largely overlooked by the existing observational analyses, is important for the interannual SST balance over the equatorial Pacific. The interannual horizontal heat advection anomalies are found to be sensitive to SST errors where oceanic currents are strong, which is a challenge for ENSO prediction. The study further suggests that the observational analyses of the tropical SST balance based on the TAO and the Reynolds SST data contain significant errors due to the large gradient errors in the Reynolds SST data, which are amplified into the advection anomalies by the large equatorial currents.

Representation of microscale ocean wavenumber spectrum correct to the second order

In this paper, the analytical representations of four wave source functions in high-frequency spectrum range are given on the basis of ocean wave theory and dimensional analysis, and the perturbation method is used to solve the governing equations of ocean wave high-frequency spectrum on the basis of the temporally stationary and locally homogeneous scale relations of microscale wave. The microscale ocean wavenumber spectrum correct to the second order has an explicit structure, its first order part represents the equilibrium between different source functions, and its second order part represents the contribution of microscale wave propagation.

Ocean wave imaging mechanism by imaging radar

Analytical representations of the high frequency spectra of ocean wave and its variation due to the variation of ocean surface current are derived from the wave-number spectrum balance equation. The ocean surface imaging formulation of real aperture radar (RAR) is given using electromagnetic wave backscattering theory of ocean surface and the modulations of ocean surface winds, currents and their variations to RAR are described. A general representation of the phase modulation induced by the ocean surface motion is derived according to standard synthetic aperture radar (SAR) imaging theory. The detectability of ocean current and sea bottom topography by imaging radar is discussed. The results constitute the theoretical basis for detecting ocean wave fields, ocean surface winds, ocean surface current fields, sea bottom topography, internal wave and so on.

Effects of winds, tides and storm surges on ocean surface waves in the Sea of Japan

Ocean surface waves are strongly forced by high wind conditions associated with winter storms in the Sea of Japan. They are also modulated by tides and storm surges. The effects of the variability in surface wind forcing, tides and storm surges on the waves are investigated using a wave model, a high-resolution atmospheric mesoscale model and a hydrodynamic ocean circulation model. Five month-long wave model simulations are inducted to examine the sensitivity of ocean waves to various wind forcing fields, tides and storm surges during January 1997. Compared with observed mean wave parameters, results indicate that the high frequency variability in the surface wind filed has very great effect on wave simulation. Tides and storm surges have a significant impact on the waves in nearshores of the Tsushima-kaihyo, but not for other regions in the Sea of Japan. High spatial and temporal resolution and good quality surface wind products will be crucial for the prediction of surface waves in the JES and other marginal seas, especially near the coastal regions.

Mitochondrial and chloroplast phylogeography of Picea crassifolia Kom. (Pinaceae) in the Qinghai-Tibetan Plateau and adjacent highlands

The disjunct distribution of forests in the Qinghai-Tibetan Plateau (QTP) and adjacent Helan Shan and Daqing Shan highlands provides an excellent model to examine vegetation shifts, glacial refugia and gene flow of key species in this complex landscape region in response to past climatic oscillations and human disturbance. In this study, we examined maternally inherited mitochondrial DNA (nad1 intron b/c and nad5 intron 1) and paternally inherited chloroplast DNA (trnC-trnD) sequence variation within a dominant forest species, Picea crassifolia Kom. We recovered nine mitotypes and two chlorotypes in a survey of 442 individuals from 32 populations sampled throughout the species' range. Significant mitochondrial DNA population subdivision was detected (G(ST) = 0.512; N-ST = 0.679), suggesting low levels of recurrent gene flow through seeds among populations and significant phylogeographical structure (N-ST > GST, P < 0.05). Plateau haplotypes differed in sequence from those in the adjacent highlands, suggesting a long period of allopatric fragmentation between the species in the two regions and the presence of independent refugia in each region during Quaternary glaciations. On the QTP platform, all but one of the disjunct populations surveyed were fixed for the same mitotype, while most populations at the plateau edge contained more than one haplotype with the mitotype that was fixed in plateau platform populations always present at high frequency. This distribution pattern suggests that present-day disjunct populations on the QTP platform experienced a common recolonization history. The same phylogeographical pattern, however, was not detected for paternally inherited chloroplast DNA haplotypes. Two chlorotypes were distributed throughout the range of the species with little geographical population differentiation (G(ST) = N-ST = 0.093). This provides evidence for highly efficient pollen-mediated gene flow among isolated forest patches, both within and between the QTP and adjacent highland populations. A lack of isolation to pollen-mediated gene flow between forests on the QTP and adjacent highlands is surprising given that the Tengger Desert has been a geographical barrier between these two regions for approximately the last 1.8 million years.

Response of chlorophyll fluorescence to dynamic light in three alpine species differing in plant architecture

A study was carried out to examine the effect of dynamic photosynthetically active photon flux density (PPFD) on photoinhibition and energy use in three herbaceous species, prostrate Saussurea superba, erect-leaved S. katochaete, and half-erect-leaved Gentiana straminea, from the Qinghai-Tibet Plateau. Chlorophyll fluorescence response was measured under each of three sets of high-low PPFD combinations: 1700-0, 1400-300, and 1200-500 mu mol m(-2) s(-1), illuminating in four dynamic frequencies: 1, 5, 15, and 60 cycles per 2 h. The total light exposure time was 2h and the integrated PPFD was the same in all treatments. The highest frequency of PPFD fluctuation resulted in the lowest photochemical activity, the highest level of non-photochemical quenching, and the greatest decrease of F-v/F-m (maximal photochemical efficiency of PSII). The 5 and 15 cycles per 2h treatments resulted in higher photochemical activity than the 1 cycle per 2h treatment. The 1700-0 PPFD combination led to the lowest photochemical activity and more serious photoinhibition in all species. S. superba usually exhibited the highest photochemical activity and CO2 uptake rate, the lowest reduction of F-v/F-m,F- and the smallest fraction of energy in thermal dissipation. With similar fractions of thermal dissipation, S. katochaete had relatively less photoinhibition than G. straminea owing to effective F-o quenching. The results suggest that high frequency of fluctuating PPFD generally results in photoinhibition, which is more serious under periods of irradiation with high light intensity. (c) 2005 Elsevier B.V. All rights reserved.

Phylogeography of the Qinghai-Tibetan Plateau endemic Juniperus przewalskii (Cupressaceae) inferred from chloroplast DNA sequence variation

The vegetation of the northeast Qinghai-Tibetan Plateau is dominated by alpine meadow and desert-steppe with sparse forests scattered within it. To obtain a better understanding of the phylogeography of one constituent species of the forests in this region, we examined chloroplast trnT-trnF and trnS-trnG sequence variation within Juniperus przewalskii, a key endemic tree species. Sequence data were obtained from 392 trees in 20 populations covering the entire distribution range of the species. Six cpDNA haplotypes were identified. Significant population subdivision was detected (G(ST) = 0.772, N-ST = 0.834), suggesting low levels of recurrent gene flow among populations and significant phylogeographic structure (N-ST > G(ST), P < 0.05). Eight of the nine disjunct populations surveyed on the high-elevation northeast plateau were fixed for a single haplotype (A), while the remaining, more westerly population, contained the same haplotype at high frequency together with two low frequency haplotypes (C and F). In contrast, most populations that occurred at lower altitudes at the plateau edge were fixed or nearly fixed for one of two haplotypes, A or E. However, two plateau edge populations had haplotype compositions different from the rest. In one, four haplotypes (A, B, D and E) were present at approximately equivalent frequencies, which might reflect a larger refugium in the area of this population during the last glacial period. Phylogenetic analysis indicated that the most widely distributed haplotype A is not ancestral to other haplotypes. The contrasting phylogeographic structures of the haplotype-rich plateau edge area and the almost haplotype-uniform plateau platform region indicate that the plateau platform was recolonized by J. przewalskii during the most recent postglacial period. This is supported by the findings of a nested clade analysis, which inferred that postglacial range expansion from the plateau edge followed by recent fragmentation is largely responsible for the present-day spatial distribution of cpDNA haplotypes within the species.

UHF RFID 射频前端及其编解码设计

射频识别（Radio Frequency Identification，RFID）技术，是一种利用射频通信实现的非接触式的数据采集和自动识别技术（以下通称RFID技术）。而超高频射频识别技术（Ultra High Frequency RFID，UHF RFID）具有识别距离远、识别准确率高、识别速度快、抗干扰能力强等特点而成为当前研发的热点。UHF RFID读写器的难点就在于射频前端电路和基带编解码的设计，它们设计的好坏直接决定了读写器的性能好坏。 本文首先通过介绍UHF RFID读写器射频前端设计的基本原理，采用射频通用收发模块进行射频前端设计的方法，给出了以ADF7020收发芯片为核心的UHF RFID读写器的射频前端的整体设计和具体的实现电路，设计了包括射频收发电路、射频前端匹配电路、滤波电路、环行器电路、功率放大电路等。 其次根据EPC Gen-2的协议标准进行了UHF RFID读写器的基带编码解码的仿真设计，然后开发了以FPGA为核心的完整的数字基带硬件电路，实际调试表明整个基带编解码软件在硬件基带PCB板上运行状况良好，并能对EPC Gen-2的协议标准的命令进行正确的编码解码。 最后通过研究学习软件无线电的理论和开发方法，把UHF　RFID读写器的射频前端分成射频模拟前端和射频数字前端，给出了一种基于软件无线电思想的UHF RFID射频数字前端设计模型，并借助于SIMULINK中的信号处理工具箱对构建的数字前端的进行仿真验证，仿真结果验证了用软件无线电实现UHF RFID数字前端的可行性。

地震数据波动方程反演方法研究

Geophysical inversion is a theory that transforms the observation data into corresponding geophysical models. The goal of seismic inversion is not only wave velocity models, but also the fine structures and dynamic process of interior of the earth, expanding to more parameters such as density, aeolotropism, viscosity and so on. As is known to all, Inversion theory is divided to linear and non-linear inversion theories. In rencent 40 years linear inversion theory has formed into a complete and systematic theory and found extensive applications in practice. While there are still many urgent problems to be solved in non-linear inversion theory and practice. Based on wave equation, this dissertation has been mainly involved in the theoretical research of several non-linear inversion methods: waveform inversion, traveltime inversion and the joint inversion about two methods. The objective of gradient waveform inversion is to find a geologic model, thus synthetic seismograms generated by this geologic model are best fitted to observed seismograms. Contrasting with other inverse methods, waveform inversion uses all characteristics of waveform and has high resolution capacity. But waveform inversion is an interface by interface method. An artificial parameter limit should be provided in each inversion iteration. In addition, waveform information will tend to get stuck in local minima if the starting model is too far from the actual model. Based on velocity scanning in traditional seismic data processing, a layer-by-layer waveform inversion method is developed in this dissertation to deal with weaknesses of waveform inversion. Wave equation is used to calculate the traveltime and derivative (perturbation of traveltime with respect to velocity) in wave-equation traveltime inversion (WT). Unlike traditional ray-based travetime inversion, WT has many advantages. No ray tracing or traveltime picking and no high frequency assumption is necessary and good result can be got while starting model is far from real model. But, comparing with waveform inversion, WT has low resolution. Waveform inversion and WT have complementary advantages and similar algorithm, which proves that the joint inversion is a better inversion method. And another key point which this dissertation emphasizes is how to give fullest play to their complementary advantages on the premise of no increase of storage spaces and amount of calculation. Numerical tests are implemented to prove the feasibility of inversion methods mentioned above in this dissertation. Especially for gradient waveform inversion, field data are inversed. This field data are acquired by our group in Wali park and Shunyi district. Real data processing shows there are many problems for waveform inversion to deal with real data. The matching of synthetic seismograms with observed seismograms and noise cancellation are two primary problems. In conclusion, on the foundation of the former experiences, this dissertation has implemented waveform inversions on the basis of acoustic wave equation and elastic wave equation, traveltime inversion on the basis of acoustic wave equation and traditional combined waveform traveltime inversion. Besides the traditional analysis of inversion theory, there are two innovations: layer by layer inversion of seimic reflection data inversion and rapid method for acoustic wave-equation joint inversion.

天然气藏地震识别技术研究

The Research of Seismic Recognition Techniques for Gas Reservoir Shang Yong_sheng(Geophysics) Directed by Yang Chang-chun Abstract Gas reservior is one of the most important nature resources. Someone forecast that the output will exceed crude oil in 2015 and become the largest energy source. Recently,more and more gas reservior are discovered as the oil field and gas filed exploration go deep into development. Although the gas proved reserves rise greatly the explorative degree of natural gas resource in our country is still very low. The potential of gas exploration is very great and our task is so hard. How to recognise and discover new gas reservoir is the first task based on the great gas reservior resources foreground. the gas reservior in different oil and gas field have its special gas generation, reservoiring, physical property conditions. However,it may have the same geophysical characters. So,it is very important to analyse, research, summarizing the geophysical characters of the gas reservior and make use of the characters to identify the gas layer effectively. This paper start with modeling,and it discuss the geophysical characters of the gas reservior response. It analyse the seismical wave characters of the gas reservoir. Furthermore, it summarize the method of using the seismica profile to identify the gas reservior directly. The paper discuss the research of extracting diffraction wave for mass diffraction wave grow at the edge of the gas reservoir at the seismic section. Making use of the technique of extracting diffraction wave to identify the gas reservior is the first experiment of the gas reservoir prediction technique. The avo technology is a new geophysical method. From the pre-stack analysis, this paper discuss the technique of using the rich information to identify the gas reservoir. Based on the case study of the Qidam basin and the Hailaer basin it discuss the method of predicating gas reservoir using pre-stack information. It include pre-stack amplitude preserve process, AVO modeling, fluid replacement technique, AVO analysis and interpretation technique. The paper summarize a gas reservoir prediction procedure focusing on the pre-stack information. The seismic wave will cause great attenuation when it pass through the gas layer and the high frequency component loss. This paper discuss the technique of extracting seismic attributes to represent the attenuation degree of seismica wave. Based on the attenuation attributes,it does the research of the gas reservor identification and prediction. At last,the paper discuss the method of calculating the azimuthal anisotropy to predict the fracture reservoir. Keyword: gas reservoir, diffraction wave, AVO, attenuation attribute,fracture prediction

反假频重建理论及其在储层预测中的应用

To extend the cross-hole seismic 2D data to outside 3D seismic data, reconstructing the low frequency data to high frequency data is necessary. Blind deconvolution method is a key technology. In this paper, an implementation of Blind deconvolution is introduced. And optimized precondition conjugate gradient method is used to improve the stability of the algorithm and reduce the computation. Then high-frequency retrieved Seismic data and the cross-hole seismic data is combined for constraint inversion. Real data processing proved the method is effective. To solve the problem that the seismic data resolution can’t meet the request of reservoir prediction in the river face thin-layers in Chinese eastern oil fields, a high frequency data reconstruction method is proposed. The extrema of the seismic data are used to get the modulation function which operated with the original seismic data to get the high frequency part of the reconstruction data to rebuild the wide band data. This method greatly saves the computation, and easy to adjust the parameters. In the output profile, the original features of the seismic events are kept, the common feint that breaking the events and adding new zeros to produce alias is avoided. And the interbeded details are enhanced compared to the original profiles. The effective band of seismic data is expended and the method is approved by the processing of the field data. Aim to the problem in the exploration and development of Chinese eastern oil field that the high frequency log data and the relative low frequency seismic data can’t be merged, a workflow of log data extrapolation constrained by time-phase model based on local wave decomposition is raised. The seismic instantaneous phase is resolved by local wave decomposition to build time-phase model, the layers beside the well is matched to build the relation of log and seismic data, multiple log info is extrapolated constrained by seismic equiphase map, high precision attributes inverse sections are produced. In the course of resolve the instantaneous phase, a new method of local wave decomposition --Hilbert transform mean mode decomposition(HMMD) is raised to improve the computation speed and noise immunity. The method is applied in the high resolution reservoir prediction in Mao2 survey of Daqing oil field, Multiple attributes profiles of wave impedance, gamma-ray, electrical resistivity, sand membership degree are produced, of which the resolution is high and the horizontal continuous is good. It’s proved to be a effective method for reservoir prediction and estimation.

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

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.

KEL-TI介质的地震波场响应与成像方法研究

To deal with the problems in multi-component converted seismic wave exploration in coal fields, the wave propagating features and imaging methods of multi-component converted waves in coal measure strata are researched in this thesis firstly. The relations between viscoelasticity and anisotropy in coal measure strata are analyzed to build KEL-TI model, and which seismic wave propagating and attenuating features are researched. The disadvantages of converted wave imaging methods based on common converted point gather are analyzed and constant velocity no NMO converted wave imaging method based on common scattering point gather is put forward, according to Huygens-Fresnel principle, which applicabilities in the elastic isotropic, elastic TI and KEL-TI situations are discussed. To different model simulation data, the common scattering point gathers’ and stacked profiles’ features are analyzed. The results show that the method can image compressional waves and converted waves with high precision. Secondly, the resolution enhancing theories and methods of converted wave are researched by Rayleigh wave suppressing, converted wave static correction and poststack inverse-Q filtering. 1) The polarization filter is designed by the instantaneous polarization information of seismic waves, and the Rayleigh wave suppressing method is researched. From the spectrum analysis before and after filtering, it can be derived that the amplitudes are kept relatively. 2) To constant velocity no NMO converted wave imaging method, the static correction method based on common equivalent offset point gather is put forward and tested to the actual converted waves. 3) The relation between equivalent quality factor of converted wave, compressional wave quality factor and the ratio of compressional to shear wave velocity is derived. The compressional wave quality inversion method by first arrivals of none-offset VSP is researched, and which is then transformed to the equivalent quality factor to perform inverse-Q filtering of actual converted waves. The result has shown that the method can recover the high frequency energy of converted waves. At last, the theories and methods researched in this thesis are practiced to the 3D3C seismic exploration in Guqiao coal mine in Huainan and achieve good results.