Dynamics of the cold-water event off the southeast coast of the United States in the summer of 2003

The cold-water event along the southeast coast of the United States in the summer of 2003 is studied using satellite data combined with in situ observations. The analysis suggests that the cooling is produced by wind-driven coastal upwelling, which breaks the thermocline barrier in the summer of 2003. The strong and persistent southwesterly winds in the summer of 2003 play an important role of lifting the bottom isotherms up to the surface and away from the coast, generating persistent surface cooling in July-August 2003. Once the thermocline barrier is broken, the stratification in the nearshore region is weakened substantially, allowing further coastal cooling of large magnitudes by episodic southerly wind bursts or passage of coastally trapped waves at periods of a few days. These short-period winds or waves would otherwise have no effects on the surface temperature because of the strong thermocline barrier in summer if not for the low-frequency cooling produced by the persistent southwesterly winds.

大气季节内振荡对东亚夏季风活动的影响研究

利用ERA40逐日再分析资料、NCEP/NCAR2逐日再分析资料、中国740个测站日降水资料、上海台风研究所提供的西太平洋热带气旋资料、Kaplan等重建的月平均SSTA资料、NOAA逐日长波辐射（OLR）等资料，应用离散功率谱分析、带通滤波、EOF分析等统计方法，研究了东亚夏季风(EASM）的移动特征、东亚地区季节内振荡(ISO)的基本特征、季节内振荡对东亚夏季风活动的影响、季节内振荡对东亚夏季风异常活动的影响机理。主要结论如下： （1）综合动力和热力因素定义了可动态描述东亚夏季风移动和强度的指数，并利用该指数研究了东亚夏季风的爆发和移动的季节内变化及其年际和年代际变化特征。研究发现，气候平均东亚夏季风前沿分别在28候、33候、36候、38候、40候、44候出现了明显的跳跃。东亚夏季风活动具有显著的年际变率，主要由于季风前沿在某些区域异常停滞和突然跨越北跳或南撤引起，造成中国东部旱涝灾害频繁发生。东亚夏季风的活动具有明显的年代际变化，在1965年、1980年、1994年发生了突变，造成中国东部降水由“南旱北涝”向“南涝北旱”的转变。 （2）东亚季风区季节内变化具有10~25d和30~60d两个波段的季节内振荡周期，以30-60d为主。存在三个主要低频模态，第一模态主要表征了EASM在长江中下游和华北地区活动期间的低频形势；第二模态印度洋-菲律宾由低频气旋式环流控制，主要表现了ISO在EASM爆发期间的低频形势；第三模态主要出现在EASM在华南和淮河活动期间的低频形势。第一模态和第三模态是代表东亚夏季风活动异常的主要低频形势。 （3）热带和副热带地区ISO总是沿垂直切变风的垂直方向传播。因此，在南海-菲律宾东北风垂直切变和副热带西太平洋北风垂直切变下，大气热源激发菲律宾附近交替出现的低频气旋和低频反气旋不断向西北传播，副热带西太平洋ISO以向西传播为主。中高纬度地区，乌拉尔山附近ISO以向东、向南移动或局地振荡为主；北太平洋中部ISO在某些情况下向南、向西传播。 （4）季风爆发期，伴随着热带东印度洋到菲律宾一系列低频气旋和低频反气旋， 冷空气向南输送，10~25天和30~60天季节内振荡低频气旋同时传入南海加快了南海夏季风的爆发。在气候态下，ISO活动表现的欧亚- 太平洋（EAP）以及太平洋-北美（PNA）低频波列分布特征(本文提出的EAP和PNA低频波列与传统意义上的二维定点相关得到的波列不同)。这种低频分布形式使得欧亚和太平洋中高纬度的槽、脊及太平洋副热带高压稳定、加强，东亚地区的低频波列则成为热带和中高纬度ISO相互作用影响东亚夏季风活动的纽带。不同的阶段表现不同的低频模态，30~60d低频模态的转变加快了EASM推进过程中跳跃性；30-60d低频模态的维持使得EASM前沿相对停滞。 （5）30-60d滤波场，菲律宾海域交替出现的低频气旋和低频反气旋不断向西北传播到南海-西太平洋一带。当南海-西太平洋地区低频气旋活跃时，季风槽加强、东伸，季风槽内热带气旋（TC）频数增加；当南海-西太平洋低频反气旋活跃时，季风槽减弱、西退，TC处于间歇期，生成位置不集中。 （6）在El Nino态下，大气季节内振荡偏弱，北传特征不明显，但ISO由中高纬度北太平洋中部向南和副热带西太平洋向西的传播特征显著，东亚地区ISO活动以第三模态为主，EASM集中停滞在华南和淮河流域，常伴随着持续性区域暴雨的出现，易造成华南和江淮流域洪涝灾害，长江和华北持续干旱。在La Nina态下，大气季节内振荡活跃，且具有明显的向北传播特征，PNA低频波列显著，东亚地区ISO活动以第一模态单峰为主；EASM主要停滞在长江中下游和华北地区，这些地区出现异常持续强降水，华南和淮河流域多干旱；在El Nino态向La Nina态转换期，ISO活动以第一模态双峰为主，长江中下游常常出现二度梅。

热带印度洋低频Rossby波及其对上层热结构的影响

Rossby波是地球物理流体动力学中非常重要的一种波动，海洋斜压长Rossby波在海洋动力过程中起着相当重要的作用。它维持并影响强西边界流，是海盆内能量传播的主要机制，它所携带的变异信号从大洋的东边界传播到内部，对海气耦合系统起到很重要的作用。热带印度洋是季风爆发的源地，对季风的年际变化具有重要影响。研究热带印度洋对理解季风变率和提高季风预测水平有重要的科学和应用价值。 本文利用TOPEX/Poseidon等高度计资料、美国国家海洋数据中心（NODC）的世界海洋图集（WOA05）长期气候态水文资料、美国Scripps海洋研究所的上层海温资料、中国Argo资料中心提供的Argo资料、美国国家环境预测中心（NCEP）的海表面温度、FSU（Florida State University）月平均风场和海气界面热通量等观测数据，全面分析了热带印度洋低频Rossby波的基本特征，并深入研究了低频Rossby波的生成机制及其对上层海洋热结构的影响。 采用相关分析等统计方法，结合1.5层约化重力模式，研究了热带南印度洋低频Rossby波的生成机制。结果表明： （1）热带南印度洋低频Rossby波分为东边界扰动产生的Rossby波和南印度洋中部风强迫Rossby波；东边界激发的为自由Rossby波，沿12°S波速大约13 cm/s，向西最远传播到80°E左右，之后被局地变量调整；强迫Rossby波在西传的过程中不断加强，波速较快，沿12°S能超过20 cm/s； （2）东边界扰动由印度尼西亚贯通流（ITF）导致的地转调整过程引起；内区风强迫Rossby波生成和加强的关键区为（70°E–95°E，15°S–5°S）；显著的西传Rossby波同太平洋上的厄尔尼诺/南方涛动（ENSO）事件紧密相连，ENSO通过大气的遥驱动机制激发热带南印度洋低频Rossby波； （3）作为东边界低频变量扰动的一个重要因子，ITF的变化与ENSO事件密切相关，总的来讲，El Niño年ITF偏弱，La Niña年ITF偏强，这与前人的研究结果一致；但它在ENSO的不同位相时期，存在一定差异，并具有夏季锁相特征：El Niño事件发生年的春季到秋季，ITF偏强，夏季最强；从El Niño盛期（冬季）到次年秋季，ITF持续偏弱，夏季最弱。上述夏季锁相特征与夏季风的强弱变化相对应。La Niña期间情况相反。 西南印度洋（SWIO）（50°E–75°E，15°S–5°S）以及苏门答腊-爪哇沿岸地区是海表面高度异常（SSHA）和海表面温度异常（SSTA）相关显著的区域，表明内部海洋动力过程在次表层和表层变量的相互关系中起重要作用。本文以2006–2008年期间三个连续的同El Niño或者La Niña同时发生的正印度洋偶极子（IOD）事件为背景，基于Argo观测资料研究了表层和次表层IOD的演变以及二者的区别和联系，并首次采用垂直模态分解方法探讨了Rossby波对上层海洋热力结构影响的动力学特征，得到如下主要结论： （1）在热带印度洋，海洋动力过程一般主要由第一和第二低阶垂直斜压模态控制，而第一斜压模态处于主导地位——在SWIO海区，第一斜压模态运动的方差解释率为第二模态的2–3倍，在赤道和东南印度洋也达到2倍左右；另外，赤道印度洋地区高阶斜压模态运动对该地区的海洋动力过程也具有一定的贡献； （2）低频斜压Rossby波能影响海洋的垂直层结，尤其是强暖Rossby波使同第二斜压模态运动紧密相连的海洋上层层结减弱，加强第二斜压模态的贡献量，导致上层各等压线向下垂直位移增大，最终通过垂直混合过程调整上层海洋的热力结构；而低频斜压冷Rossby波会加强上层垂直层结，抑制该层内变量变化，因此第二斜压模态的贡献依然很小； （3）表层IOD和次表层IOD分布形态不同：表层东部冷异常主要集中在东南印度洋Sumatra-Java沿岸，次表层冷异常基本关于赤道对称；表层西部暖异常基本关于赤道对称，而在次表层赤道以南海温扰动强度远远大于赤道以北； （4）正IOD事件中，东南印度洋冷SSTA首先出现于Java沿岸，沿岸东南风引起的潜热释放增加以及沿岸上升流是该初始冷异常建立的主要机制，与之相关的SSTA东西梯度加强大气环流变化，并进一步强迫随后的海洋运动；1–2个月后，SST冷异常中心北跳到Sumatra沿岸并向西扩展，同时不断增强，其中Sumatra沿岸上升流、来自赤道印度洋的冷Kelvin波及其反射的西传冷Rossby波是这一演变过程的动力机制，而沿岸上升流起决定作用。

Phylogeography of Pinus tabulaeformis Carr. (Pinaceae), a dominant species of coniferous forest in northern China

How coniferous trees in northern China changed their distribution ranges in response to Quaternary climatic oscillations remains largely unknown. Here we report a study of the phylogeography of Pinus tabulaeformis, an endemic and dominant species of coniferous forest in northern China. We examined sequence variation of maternally inherited, seed-dispersed mitochondrial DNA (mtDNA) (nad5 intron 1 and nad4/3-4) and paternally inherited, pollen- and seed-dispersed chloroplast DNA (cpDNA) (rpl16 and trnS-trnG) within and among 30 natural populations across the entire range of the species. Six mitotypes and five chlorotypes were recovered among 291 trees surveyed. Population divergence was high for mtDNA variation (G(ST) = 0.738, N-ST = 0.771) indicating low levels of seed-based gene flow and significant phylogeographical structure (N-ST > G(ST), P < 0.05). The spatial distribution of mitotypes suggests that five distinct population groups exist in the species: one in the west comprising seven populations, a second with a north-central distribution comprising 15 populations, a third with a southern and easterly distribution comprising five populations, a fourth comprising one central and one western population, and a fifth comprising a single population located in the north-central part of the species' range. Each group apart from the fourth group is characterized by a distinct mitotype, with other mitotypes, if present, occurring at low frequency. It is suggested, therefore, that most members of each group apart from Group 4 are derived from ancestors that occupied different isolated refugia in a previous period of range fragmentation of the species, possibly at the time of the Last Glacial Maximum. Possible locations for these refugia are suggested. A comparison of mitotype diversity between northern and southern subgroups within the north-central group of populations (Group 2) showed much greater uniformity in the northern part of the range both within and between populations. This could indicate a northward migration of the species from a southern refugium in this region during the postglacial period, although alternative explanations cannot be ruled out. Two chlorotypes were distributed across the geographical range of the species, resulting in lower levels of among-population chlorotype variation. The geographical pattern of variation for all five chlorotypes provided some indication of the species surviving past glaciations in more than one refugium, although differentiation was much less marked, presumably due to the greater dispersal of cpDNA via pollen.

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.

基于复合喷气压电驱动器的柔性机械臂振动控制

空间机器人和大型柔性空间结构在航天器调姿、变轨、外部扰动的情况下将引起振动问题,其低频大幅值振动将持续很长时间,这将影响航天器系统的稳定性和控制精度。为了快速抑制低频大幅值振动及残余振动,提出采用复合可控反作用力幅值的喷气式驱动和压电陶瓷驱动方案进行振动控制。进行基于复合控制的柔性臂系统动力学建模并给出控制算法。设计并建立柔性机械臂试验平台,构建气动驱动控制回路及压电驱动控制回路。进行基于压电陶瓷驱动器、喷气式驱动器及复合喷气和压电驱动器的柔性臂大幅值低频模态振动控制的几种方法试验比较研究。试验结果表明,采用的控制方案和方法既可以快速地抑制柔性机械臂统的低频大幅值振动,又明显地同时抑制高频和低频小幅值残余振动。

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

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.

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

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.

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

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.

青藏高原东部长周期大地电磁测深研究

As a complement to conventional MT, Long-period Magnetotellurics (LMT) has been developed at low frequency for soundings of deep electric structures. Eastern Himalayan Syntaxis (EHS) and surrounding area is a key place for the study of dynamics of the uplift of Tibetan plateau. Experiments in the pioneer studies for EHS3D project showed that the study area shares an unusual low resistive crust and upper mantle. Conventional MT could not provide sufficient information about the deep structures of the study area that requested long period MT measurement to be complemented. This thesis presents the LMT studies in eastern Tibet along the EHS3D-3 Profile from Xiachayu to Yushu including data acquisition, processing, inversion and interpretation. The effective period of the measured LMT signals extend from 10s up to 30000s for the duration more than one week measurement. The resulting model shows that the LMT sounding coincides with the MT data in overlapped periods. Especially the induction arrows and tippers derived from LMT data provide more information about the base of the conductors beneath the plateau with higher resolution. Anomalous induction coefficients and 2-D model suggest extensive conductive bodies beneath Lhasa block and Qiangtang terrain which would be a possible evidence for partial melt and fluids at depth.

地震资料叠前储层预测综合研究--以辽中凹陷东二下段浊积扇气藏为例

Reservoir prediction techniques from prestack seismic are among the most important ones for exploration of lithologic hydrocarbon reservoir. In this paper, we set the turbidite fan sandstone reservoir in Liao-Zhong depress as our researching target, and aims to solve the apllication difficulties on pre-stack inversion in the area, where the drilling data is scarce and the reservoir is lateral varied. Meanwhile, AVO analysis and pre-stack inversion for gas-bearing detection is systematically researched. The seismic reflection characters of gas-bearing sandstone in turbidite fan with different fluid content are defined, after analyzing results from AVO seismic simulation and porous fluid replacement of real log data, and under the guides of the seismic characters from classical gas-bearing sandstone reservoir and numerical simulation for complicate gas-bearing sandstone. It is confirmed that detecting gas-bearing sandstone in turbidite fan via AVO technologies is feasible. In terms of AVO analysis, two AVO characters, fluid detection factor and product of intercept and gradient, can effectively identify top and bottom boundaries and lateral range of tuibidite gas sand by comparing real drilling data. Cross-plotting of near and far angle stack data could avoid the correlation existing in P-G analysis. After comparing the acoustic impedance inversions with routine stacked data and AVO intercept, impedance derived from AVO intercept attribute could reduce the acoustic impedance estimating error which is caused by AVO. On the aspect of elastic impedance inversion, the AVO information in the pre-stack gathers is properly reserved by creating partial angle stack data. By the far angle elastic impedance alone, the gas sand, with abnormally low range of values, can be identified from the background rocks. The boundary of gas sand can also be clearly determined by cross-plotting of near and far angle elastic impedances. The accuracy of far angle elastic impedance is very sensitive to the parameter K, and by taking the statistical average of Vp/Vs on the targeted section in key wells, the accuracy of low frequency trends is gurranteed; the intensive absorsion within the area of the gas sand, which tends to push the spectral of seismic data to the lower end, will cause errors on the inversion result of elastic impedance. The solution is to confine the inversion on the interested area by improving the wavelet. On the aspect of prestack AVA simultaneous inversion, the constraint of local rock-physical trends between velocities of P-wave、S-wave and density successfully removes the instability of inversion, thus improves the precision of the resulting elastic parameters. Plenty of data on rock properties are derived via AVO analysis and prestack seismic data inversion. Based on them, the fluid anomaly is analysized and lithological interpretation are conducted. The distribution of gas sand can be consistently determined via various of ways, such as cross-plotting of P and G attributes, near and far partial angle stack data, near and far angle elastic impedances, λρ and Vp/Vs, etc. The shear modulo and density are also reliable enough to be used for lithological interpretation. We successfully applied the AVO analysis and pre-stack inversion techniques to gas detecting for turbidite fan sand reservoir in Liao-Zhong depression.

含激电信息的人工源频率域电磁法正演研究——大功率人工固定源极低频电磁法激电效应预研究

ABSTRACT Base on Dias-model, this paper study the forward modeling to IP response with MT souece and CSAMT source with the author’s program，that is to provide the theoretical basis for extracting the IP information from the far field, near field and transition field in extremely low frequency electromagnetic sounding with artificial strong fixed source signal. The outline of the paper is as follows. Firstly, review the history of the complex-resistivity-model, and analyses the bases for choosing Dias model; meanwhile, effects and responses of each parameter in Dias model are analyses. Afterward, study the forward modeling to 1D layered model with MT source, numerically simulate the IP effects of classic geoelectric objects; Find the clear anomaly and relationship between the peak value of the amplitude anomaly ratio and phase anomaly difference with and without IP parameters to be considered within a frequency range. On the basis of the modeling of MT, Study the CSAMT modeling with a dipole-source, obtain the anomaly responses and the relationship between anomaly and characters of object. Base on the infinite line source, study the 2D IP effect of geoelectric objects with 2D modeling. Calculate the response with different source distance, object depth and different wall rock resistivity systemically. Finally, conclude the viewpoints and give the discussion of the result, and point out the lack of this research also.

粘弹性介质中的地震波传播与波形反演研究

The real earth is far away from an ideal elastic ball. The movement of structures or fluid and scattering of thin-layer would inevitably affect seismic wave propagation, which is demonstrated mainly as energy nongeometrical attenuation. Today, most of theoretical researches and applications take the assumption that all media studied are fully elastic. Ignoring the viscoelastic property would, in some circumstances, lead to amplitude and phase distortion, which will indirectly affect extraction of traveltime and waveform we use in imaging and inversion. In order to investigate the response of seismic wave propagation and improve the imaging and inversion quality in complex media, we need not only consider into attenuation of the real media but also implement it by means of efficient numerical methods and imaging techniques. As for numerical modeling, most widely used methods, such as finite difference, finite element and pseudospectral algorithms, have difficulty in dealing with problem of simultaneously improving accuracy and efficiency in computation. To partially overcome this difficulty, this paper devises a matrix differentiator method and an optimal convolutional differentiator method based on staggered-grid Fourier pseudospectral differentiation, and a staggered-grid optimal Shannon singular kernel convolutional differentiator by function distribution theory, which then are used to study seismic wave propagation in viscoelastic media. Results through comparisons and accuracy analysis demonstrate that optimal convolutional differentiator methods can solve well the incompatibility between accuracy and efficiency, and are almost twice more accurate than the same-length finite difference. They can efficiently reduce dispersion and provide high-precision waveform data. On the basis of frequency-domain wavefield modeling, we discuss how to directly solve linear equations and point out that when compared to the time-domain methods, frequency-domain methods would be more convenient to handle the multi-source problem and be much easier to incorporate medium attenuation. We also prove the equivalence of the time- and frequency-domain methods by using numerical tests when assumptions with non-relaxation modulus and quality factor are made, and analyze the reason that causes waveform difference. In frequency-domain waveform inversion, experiments have been conducted with transmission, crosshole and reflection data. By using the relation between media scales and characteristic frequencies, we analyze the capacity of the frequency-domain sequential inversion method in anti-noising and dealing with non-uniqueness of nonlinear optimization. In crosshole experiments, we find the main sources of inversion error and figure out how incorrect quality factor would affect inverted results. When dealing with surface reflection data, several frequencies have been chosen with optimal frequency selection strategy, with which we use to carry out sequential and simultaneous inversions to verify how important low frequency data are to the inverted results and the functionality of simultaneous inversion in anti-noising. Finally, I come with some conclusions about the whole work I have done in this dissertation and discuss detailly the existing and would-be problems in it. I also point out the possible directions and theories we should go and deepen, which, to some extent, would provide a helpful reference to researchers who are interested in seismic wave propagation and imaging in complex media.

接收函数方法与中国及邻区上地幔间断面的研究

The receiver function method applied in researching the discontinuities in upper mantle was systematically studied in this paper. Using the theoretical receiver functions, the characteristics of P410S and P660S phases were analyzed, and the influencing factors for detection of these phases were discussed. The stability of receiver function was studied, and a new computational method of receiver function, RFSSMS (Receiver Function of Stack and Smooth of Multi seismic-records at a Single station), was put forward. We built initial reference velocity model for the media beneath each of 18 seismic stations respectively; then estimated the buried depths of 410-km and 660-km discontinuities(simply marked as '410' and '660') under the stations by using the arrive time differences of P410S and P660S with P. We developed a new receiver function inversion method -PGARFI (Peeling-Genetic Algorithm of Receiver Function Inversion), to obtain the whole crust and upper mantle velocity structure and the depths of discontinuities beneath a station. The major works and results could be summarized as follows: (1) By analysis of the theoretical receiver functions with different velocity models and different ray parameters, we obtain the knowledge: The amplitudes of P410S and P660S phases are decreasing with the increasing of epicentral distance A , and the arrival time differences of these phases with P are shorter as A is longer. The multiple refracted and/or reflected waves yielded on Moho and the discontinuities in the crust interfere the identification of P410S. If existing LVZ under the lithosphere, some multiple waves caused by LVZ will interfere the identification of P410S. The multiple waves produced by discontinuity lied near 120km depth will mix with P410s phase in some range of epicentral distance; and the multiple waves concerned with the discontinuity lied near 210km depth will interfere the identification of P660S. The epicentral distance for P4i0s identification is limited, the upper limit is 80° . The identification of P660S is not restricted by the epicenter distance obviously. The identification of P410S and P6gos in the theoretical receiver functions is interfered weakly from the seismic wave attenuation caused by the media absorption if the Q value in a reasonable range. (2) The stability of receiver function was studied by using synthetic seismograms with different kind of noise. The results show that on the condition of high signal-noise-ratio of seismic records, the high frequency background noise and the low frequency microseism noise do not influence the calculating result of receiver function. But the media "scattering noise" influence the stability of receiver function. When the scattering effect reach some level, the identification of P4iOs and P66os is difficult in single receiver function which is yielded from only one seismic record. We provided a new method to calculate receiver function, that is, with a group of earthquake records, stacking the R and Z components respectively in the frequency domain, and weighted smooth the stacked Z component, then compute the complex spectrum ratio of R to Z. This method can improve the stability of receiver function and protrude the P4i0s and P66os in the receiver function curves. (3) 263 receiver functions were provided from 1364 three component broadband seismograms recorded at 18 stations in China and adjacent areas for the tele-earthquakes. The observed arrival time differences of P410S and P660S with P were obtained in these receiver functions. The initial velocity model for every station was built according to the prior research results. The buried depths of '410' and '660' under a station were acquired by the way of adjusting the depths of these two discontinuities in the initial velocity model until the theoretical arrival time differences of P410S and P660S with P well conformed to the observed. The results show an obvious lateral heterogeneity of buried depths of ' 410' and (660' . The depth of '410' is shallower beneath BJI, XAN, LZH and ENH, but deeper under QIZ and CHTO, and the average is 403km . The average depth of '660' is 663km, deeper under MDJ and MAJO, but shallower under QIZ and HYB. (4) For inversing the whole crust and upper mantle velocity structure, a new inversion method -PGARFI (Peeling-Genetic Algorithm of Receiver Function Inversion) has- been developed here. The media beneath a station is divided into segments, then the velocity structure is inversed from receiver function from surface to deep successively. Using PGARFI, the multi reflection / refraction phases of shallower discontinuities are isolated from the first order refraction transform phase of deep discontinuity. The genetic algorithm with floating-point coding was used hi the inversion of every segment, and arithmetical crossover and non-uniform mutation technologies were employed in the genetic optimization. 10 independent inversions are completed for every segment, and 50 most excellent velocity models are selected according to the priority of fitness from all models produced in the inversion process. The final velocity structure of every segment is obtained from the weighted average of these 50 models. Before inversion, a wide range of velocity variation with depth and depth range of the main discontinuities are given according to priori knowledge. PGARFI was verified with numerical test and applied in the inversion of the velocity structure beneath HIA station down to 700km depth.

基于信息重建技术的地震信号处理方法研究

The dissertation addressed the problems of signals reconstruction and data restoration in seismic data processing, which takes the representation methods of signal as the main clue, and take the seismic information reconstruction (signals separation and trace interpolation) as the core. On the natural bases signal representation, I present the ICA fundamentals, algorithms and its original applications to nature earth quake signals separation and survey seismic signals separation. On determinative bases signal representation, the paper proposed seismic dada reconstruction least square inversion regularization methods, sparseness constraints, pre-conditioned conjugate gradient methods, and their applications to seismic de-convolution, Radon transformation, et. al. The core contents are about de-alias uneven seismic data reconstruction algorithm and its application to seismic interpolation. Although the dissertation discussed two cases of signal representation, they can be integrated into one frame, because they both deal with the signals or information restoration, the former reconstructing original signals from mixed signals, the later reconstructing whole data from sparse or irregular data. The goal of them is same to provide pre-processing methods and post-processing method for seismic pre-stack depth migration. ICA can separate the original signals from mixed signals by them, or abstract the basic structure from analyzed data. I surveyed the fundamental, algorithms and applications of ICA. Compared with KL transformation, I proposed the independent components transformation concept (ICT). On basis of the ne-entropy measurement of independence, I implemented the FastICA and improved it by covariance matrix. By analyzing the characteristics of the seismic signals, I introduced ICA into seismic signal processing firstly in Geophysical community, and implemented the noise separation from seismic signal. Synthetic and real data examples show the usability of ICA to seismic signal processing and initial effects are achieved. The application of ICA to separation quake conversion wave from multiple in sedimentary area is made, which demonstrates good effects, so more reasonable interpretation of underground un-continuity is got. The results show the perspective of application of ICA to Geophysical signal processing. By virtue of the relationship between ICA and Blind Deconvolution , I surveyed the seismic blind deconvolution, and discussed the perspective of applying ICA to seismic blind deconvolution with two possible solutions. The relationship of PC A, ICA and wavelet transform is claimed. It is proved that reconstruction of wavelet prototype functions is Lie group representation. By the way, over-sampled wavelet transform is proposed to enhance the seismic data resolution, which is validated by numerical examples. The key of pre-stack depth migration is the regularization of pre-stack seismic data. As a main procedure, seismic interpolation and missing data reconstruction are necessary. Firstly, I review the seismic imaging methods in order to argue the critical effect of regularization. By review of the seismic interpolation algorithms, I acclaim that de-alias uneven data reconstruction is still a challenge. The fundamental of seismic reconstruction is discussed firstly. Then sparseness constraint on least square inversion and preconditioned conjugate gradient solver are studied and implemented. Choosing constraint item with Cauchy distribution, I programmed PCG algorithm and implement sparse seismic deconvolution, high resolution Radon Transformation by PCG, which is prepared for seismic data reconstruction. About seismic interpolation, dealias even data interpolation and uneven data reconstruction are very good respectively, however they can not be combined each other. In this paper, a novel Fourier transform based method and a algorithm have been proposed, which could reconstruct both uneven and alias seismic data. I formulated band-limited data reconstruction as minimum norm least squares inversion problem where an adaptive DFT-weighted norm regularization term is used. The inverse problem is solved by pre-conditional conjugate gradient method, which makes the solutions stable and convergent quickly. Based on the assumption that seismic data are consisted of finite linear events, from sampling theorem, alias events can be attenuated via LS weight predicted linearly from low frequency. Three application issues are discussed on even gap trace interpolation, uneven gap filling, high frequency trace reconstruction from low frequency data trace constrained by few high frequency traces. Both synthetic and real data numerical examples show the proposed method is valid, efficient and applicable. The research is valuable to seismic data regularization and cross well seismic. To meet 3D shot profile depth migration request for data, schemes must be taken to make the data even and fitting the velocity dataset. The methods of this paper are used to interpolate and extrapolate the shot gathers instead of simply embedding zero traces. So, the aperture of migration is enlarged and the migration effect is improved. The results show the effectiveness and the practicability.