波形反演遗传算法及其在地震测深数据解释中的应用

为了更好地利用地震测深波形数据,提出了地震体波波形反演的遗传算法. 正演使用能精确快速计算互层结构响应的广义反、透射系数理论地震图算法;反演采用遗传算法,实现了地震体波波形反演的遗传算法. 数值试验表明,该算法具有分辨壳内低速层、高低速薄互层结构和一定的抗噪能力. 青藏高原东北缘泽库、夏河、临洮3炮地震测深P波波形反演,得到了上地壳底部低速层和中、下地壳,以及上地幔顶部薄互层的细结构图象.      

Crustal structure of northeastern Tibetan plateau and Ordos block: Waveform interpretation of the Maqen-Jingbian seismic refraction profile

The Maqen-Jingbian wide-angle seismic reflection and refraction experiment was carried out in 1998, which aims at determining detailed structure in the crust and top of the upper mantle and understanding structural relation between the northeastern Tibetan plateau and the Ordos block. The 1-D crustal models inferred by waveform inversion show strong variations in crustal structure, which can be classified into four different types: ① an Ordos platform with the Proterozoic crust and two high-velocity layers in the northeast section, ② a transitional crust between the northeastern Tibetan plateau and the Ordos block across the Haiyuan earthquake zone, ③ the Qilian orogenic zone in the central part, and 4 the Qinling orogenic zone in the southwestern section. The Moho depth increases from ~42 km to ~62 km from the NE part to the SW part of the profile. The crystalline crust consists of the upper crust and lower crust in northeastern Tibetan plateau. There is an obviously low P-wave velocity layer dipping northeastward, which is 12–13 km thick, at the bottom of the upper crust in Qinling orogenic zone and Haiyuan earthquake zone. The lower crust is characterized by alternating high and low P-wave velocity layers. Beneath Ordos block, i.e., the NE part of the profile, the crust shows quite a smooth increase in P-wave velocity down to the Moho at a depth of about 42 km.     

Crustal velocity structure of central Gansu Province from regional seismic waveform inversion using firework algorithm

The firework algorithm (FWA) is a novel swarm intelligence-based method recently proposed for the optimization of multi-parameter, nonlinear functions. Numerical waveform inversion experiments using a synthetic model show that the FWA performs well in both solution quality and efficiency. We apply the FWA in this study to crustal velocity structure inversion using regional seismic waveform data of central Gansu on the northeastern margin of the Qinghai-Tibet plateau. Seismograms recorded from the moment magnitude (M _W) 5.4 Minxian earthquake enable obtaining an average crustal velocity model for this region. We initially carried out a series of FWA robustness tests in regional waveform inversion at the same earthquake and station positions across the study region, inverting two velocity structure models, with and without a low-velocity crustal layer; the accuracy of our average inversion results and their standard deviations reveal the advantages of the FWA for the inversion of regional seismic waveforms. We applied the FWA across our study area using three component waveform data recorded by nine broadband permanent seismic stations with epicentral distances ranging between 146 and 437 km. These inversion results show that the average thickness of the crust in this region is 46.75 km, while thicknesses of the sedimentary layer, and the upper, middle, and lower crust are 3.15, 15.69, 13.08, and 14.83 km, respectively. Results also show that the P-wave velocities of these layers and the upper mantle are 4.47, 6.07, 6.12, 6.87, and 8.18 km/s, respectively.     

Crustal structure of northeastern Tibetan plateau and Ordos block:Waveform interpretation of the Maqen-Jingbian seismic refraction profile

The Maqen-Jingbian wide-angle seismic reflection and refraction experiment was carried out in 1998,which aims at determining detailed structure in the crust and top of the upper mantle and understanding structural relation between the northeastern Tibetan plateau and the Ordos block. The 1-D crustal models inferred by waveform inversion show strong variations in crustal structure,which can be classified into four different types:1 an Ordos platform with the Proterozoic crust and two high-velocity layers in ...     

川滇地区速度结构的区域地震波形反演研究

利用云南数字地震台网的区域地震波形资料,对川滇地区的地壳上地幔速度结构进行了初步研究. 结果表明,川滇地区上地幔顶部P波速度较小,约78 km/s,P波速度在上地幔表现为较小的正速度梯度,S波在100～160 km深度范围内表现为弱低速层. 对于较短的观测路径,不同路径的平均P波和S波速度存在明显的横向变化. 与川滇菱形块体内部的速度结构不同,在块体边界附近可以观测到比较明显的上地壳低速层,我们认为它可能与块体边界的断裂带有关;川滇菱形块体内部存在的下地壳低速层,有利于块体向南滑动,而中上地壳没有明显低速结构,可能表明川滇菱形块体向南滑动的解耦深度至少在下地壳. 根据不同路径的反演结果,给出了云南中部地区地壳内部的平均速度结构.     

逐步波形反演方法及应用

本文针对天然地震波形反演面临的困难及其复杂性，提出了逐步波形反演方法，第一步，运用波形反演中的试错法，求得地震台站下方成层介质的初步结构；第地一步，以第一步结果为初值，令各层厚度不变，反演速度；第三步，以第二步结果为初值，令速度不变，反演厚度。以上各步还可交替进行，直至得到满意结果。     

Joint inversion of teleseismic P waveforms and surface-wave group velocities from ambient seismic noise in the Bohemian Massif

Joint inversion of teleseismic P-waveforms and local group velocities of surface waves retrieved from ambient seismic noise has been performed to model velocity structure of the crust and uppermost mantle of the Bohemian Massif. We analysed P-waveforms of 381 teleseismic earthquakes recorded at 54 broadband seismic stations located on the territory of the Czech Republic and in its close surroundings. Group velocities of Rayleigh and Love surface waves were obtained by cross-correlating long-term recordings of seismic noise. The basis for waveform inversion is the well-known methodology of P-to-S receiver functions constructed from converted phases. Due to instabilities in direct inversion of receiver functions caused by the necessity of applying deconvolution, we propose an alternative formulation to fit observed and calculated radial components of P waveforms. The joint inversion is transformed into a search for the minimum of the cost function defined as a weighted sum of waveform and group velocity misfits. With the use of the robust stochastic optimizer (Differential Evolution Algorithm), neither derivatives nor a starting model are needed. The task was solved for 1D layered isotropic models of the crust and the uppermost mantle. We have performed a sequence of inversions with models containing one, two, three and four layers above a half-space. By using statistical criteria (F-test) we were able to select the simplest velocity models satisfying data and representing local geological structures. Complex crustal models are typical for stations located close to boundaries of major tectonic units. The relatively low average P to S wave-velocity ratio is in agreement with the generally accepted view that the BM crust is predominantly felsic.     

地震波形反演技术在砂泥岩薄互层结构表征中的应用

基于传播矩阵理论开发砂泥岩薄互层地震合成记录算法,与褶积算法、基于界面模型的Zoeppritz方法以及波动方程等方法相比,该方法更适用于具有复杂结构的薄互层模型,能够在充分考虑地震反射波动力学因素的同时不受网格间距的限制.基于正演算法开发了基于波形对比的砂泥岩薄互层地震反演技术,由地震反射波形特征的变化反演薄地层单元中砂体含量与空间位置等参数,进而确定薄互层段砂泥岩组合结构以及砂体的空间展布.理论模型验证了反演方法的有效性.通过测井分析建立薄互层地震地质模型,并将该技术应用于研究区实际地震数据,反演的砂体空间分布与测井资料进行对比分析,验证了反演方法的实用性.     

2-D crustal structure of Changbaishan-Tianchi volcanic region determined by seismic traveltime inversion

IntroductionTrial-and-error forward modeling of wide-angle seismic reflection/refraction traveltimes for 2-D velocity structure is extremely time-consuming, even for experienced data interpreters. For wide-angle seismic reflection/refraction experiments that consist of numerous shots along a single line, it is quite difficult through repeated trial-and-error forward modeling to construct a 2-D model that fits the data within acceptable limits (Cerveny, et al, 1977; ZHANG, et al, 200 . In ad…     

Joint inversion of travel times and waveforms by means of annealing algorithms. Application to a seismic refraction experiment

An algorithm of annealing is applied to a joint inversion of travel times and waveforms belonging to a synthetically generated seismic refraction experiment. The medium (crust and upper mantle) is modelled by a set of plane stratified layers and a halfspace. The obtained structure (elastic parameters and depth of layers) shows that, in spite of contaminating seismic noise and poor knowledge of the seismic source, annealing methods are a good tool in these kinds of inversion problems. We think that many characteristics of the annealing process described here could be used with real data and more sophisticated media for the crust and upper mantle of the earth than the present example.     

甘肃东南部地壳速度结构的区域地震波形反演

利用2007年完成扩建的甘肃东南部及邻近地区的24个宽频带固定地震台站记录到的2008年8月1日汶川地震余震的三分量地震全波形资料,采用小生境遗传算法和反射率法结合的波形反演方法,给出了甘肃东南部两个分区的地壳速度模型.西区和东区分别对应青藏高原块体和它与鄂尔多斯块体之间的过渡带,反演给出的平均模型显示,两个区上、中地壳的平均P波速度非常接近,由酸性岩和中性岩组成,下地壳P波速度差别较大,东区为6.41 km/s,西区为6.96 km/s,厚度相差也较大,东区为9.3 km,西区为19.8 km,地壳厚度由西向东减小,分别为54.6 km和47.9 km.显示西区下地壳由基性岩组成,而东区下地壳由中性岩组成,支持研究区内青藏高原东北缘地壳增厚主要发生在下地壳的观点.西区的上地幔顶部P波速度为7.73 km/s,对应年轻的构造活动区,而东区为8.05 km/s,对应稳定的古老地块.东区在上地壳上部存在厚约6.5 km的沉积层,P波速度为5.78 km/s,但是该沉积层在西部缺失.和PREM模型给出的全球平均地壳速度值相比,两个分区的地壳速度值整体偏低.     

云南思茅—中甸地震剖面的地壳结构

云南思茅—中甸宽角反射/折射地震剖面切割松潘—甘孜、扬子和华南三个构造单元的部分区域. 我们利用初至波和壳内反射波走时层析成像获得地壳纵波速度结构. 在获得新的地壳速度结构模型基础上，利用地震散射成像思想和低叠加次数的叠前深度偏移方法重建了研究区的地壳、上地幔反射结构. 综合分析研究区地壳P波速度模型和壳内地震反射剖面发现：沿测线从北至南地壳厚度从约50 km减薄至35 km左右，地壳厚度的减薄量主要体现在下地壳，剖面北段下地壳厚度约为30 km，剖面南段下地壳厚度仅为15 km左右；上地幔顶部局部位置P波速度值偏低，一般为76～78 km/s，反映出云南地区是典型的构造活动区的特点.剖面沿线地壳内地震反射发育，其中莫霍强反射出现在景云桥下方；在景云桥弧形断裂带8～10 km深处出现宽约50 km的强反射带.     

陕西省数字地震台下方壳幔速度结构研究

利用广义反、透射系数的传递矩阵计算理论地震图，用快速模拟退火法搜索最佳模型，通过与地震台站实际记录的P波波形进行拟合对比，反演了陕西省13个数字地震台下方的P波速度结构。台站分布区域涉及到秦岭造山带、渭河盆地及鄂尔多斯地台三个构造区域。结果显示，3个构造区的地下速度结构差异明显，呈现了各自的构造特点。其中秦岭造山带上地壳乃至地表速度高，中下地壳出现多个低速层，可能与岩石脱水、构造滑脱带等有关；渭河盆地内上中地壳高低速度层交替出现，反映了断层活动在浅部的影响；鄂尔多斯地台有部分台站保持稳定的匀速构造。反映了古老地台地层稳定发展的特性，部分台站受到后期改造，低速层位增多。因为变化剧烈的台站分布在盆地的南部和西部。而且震源机制解显示该区主应力场方向为近东西向，可以说秦岭北缘的断裂作用以及后来来自青藏高原的挤压应力是导致渭河盆地形成和发展的主要因素。最后我们讨论了速度结构与强震、地震分布以及不同构造演化的关系。     

太行山东缘汤阴地堑地壳结构和活动断裂探测

采用深、浅地震反射和钻孔地质剖面相结合的探测方法,对太行山东缘汤阴地堑的地壳结构和隐伏活动断裂进行了研究.结果表明,该区地壳厚度约36~42 km,莫霍面从华北平原区向太行山下倾伏.汤阴地堑是一个受汤东断裂控制的半地堑构造,其基底面形态与莫霍面展布呈"镜像"关系.汤东断裂是1条继承性的隐伏活动断裂,该断裂向上错断了埋深约20 m的中更新世晚期地层,向下延伸至上地壳底部.综合分析深地震反射和已有深地震宽角反射/折射剖面结果,发现深地震反射剖面上的中-下地壳强反射层和壳幔过渡带反射,与深地震宽角反射/折射剖面上出现的中-下地壳正负速度梯度变化层有着较好的对应关系,这表明本区中-下地壳和壳幔过渡带可能为一系列速度递变层或高低速物质的互变层,埋深约15~16 km的强反射带为上地壳与中-下地壳的转换带,壳幔过渡带的底界为地壳与地幔的分界.研究结果为深入理解该区的深部动力学过程、分析研究深浅构造关系、评价断裂的活动性提供了依据.     

大同阳高震区及其邻区壳幔速度结构与深部构造

利用通过本区６条宽角反射／折射剖面资料对大同阳高震区及邻区地壳上地幔速度结构与构造进行了详细的研究。结果表明,地壳上地幔速度结构与构造在纵向和横向上具有明显的不均一性。浅部基底断裂发育,而在其深部,根据波组特征、壳内界面及速度等值线起伏变化和低速异常体的边界等推测有３处地壳深断裂带。本区最明显的上地壳低速体位于大同—阳原附近,其南界存在地壳深断裂,大同阳高地震群与该低速异常体和深断裂有关。     

根据接收函数反演得到的首都圈地壳上地幔三维S波速度结构

利用2002~2003年中国地震局地质研究所台阵实验室以唐山大震区为中心布设的40个流动宽频带地震台站和首都圈数字台网的33个宽频带台站的远震数据,采用接收函数非线性反演方法得到其中72个宽频带台站下方60 km深度范围内的S波速度结构.根据得到的各台站下方地壳上地幔的S波速度结构,并综合刘启元等（1997）用接收函数非线性反演方法得到的延怀盆地15个宽频带流动台站下方的地壳上地幔S波速度结构模型,给出了39°N～41°N,114°E~119.5°E区域内沿不同走向、不同深度S波速度分布.由于综合了利用首都圈数字地震台网的宽频带台站以及流动地震台阵的观测数据,本文给出了较前人同类研究空间分辨率更好的结果.结果表明: (1)研究区的速度结构,特别是怀来以东的速度结构十分复杂.在10~20 km深度范围内,研究区地壳具有高速和低速异常块体的交错结构.研究区中上地壳速度结构主要被与张渤地震带大体重合的NW向高速条带和穿越唐山大震区的NE向高速条带所控制,而其中下地壳的速度结构主要为延怀—三河—唐山地区上地幔隆起所控制.(2)研究区内存在若干壳内S波低速体,它们主要分布在唐山,三河及延怀盆地等地区.在这些地区,壳内低速体伴随着壳幔界面的隆起和上地幔顶部速度结构的横向变化.(3)地表断层分布与地壳速度结构分区有较好的相关性,表明断层对不同块体有明显的控制作用.其中,宝坻断裂,香河断裂和唐山断裂均为超壳断裂.(4)首都圈内大地震的分布与壳内低速体及上地幔顶部的速度结构有密切关系.对于唐山大地震的成因,仅考虑板块作用引起的水平应力场是不够的,有必要充分重视由于上地幔变形引起的地壳垂直变形和上地幔物质侵入造成的热效应.     

Annealing algorithms in computing stress tensors: A set of numerical simulations for local stress tensor recovery

The determination of a local lithospheric stress pattern from low magnitude seismic activity recorded at a loca array is attempted by means of numerical simulations that include waveform inversion and annealing inversion techniques. A set of flat stratified layers are assumed for the propagating medium. The tests show that the annealing algorithm is a useful tool in this kind of study where collected seismic activity at local arrays could help us better understand local stress distributions in the lithosphere. Although real data recorded at local seismic arrays should be processed, assuming more realistic crust and upper mantle models, many characteristics of the annealing process described here could be used when studying real problems.     

基于Love波相速度反演南北地震带地壳上地幔结构

收集了南北地震带区域地震台网中292个地震台站2008年1月至2011年3月期间的地震波形数据,由频时分析方法提取了Love波相速度频散曲线,经过反演得到了研究区内的Love波相速度分布.根据Love波纯路径频散,采用线性反演方法对0.25°×0.25°的网格点进行了一维S波速度结构反演,利用线性插值获取了南北地震带地区的三维S波速度结构.结果显示了松潘—甘孜地体和川滇菱形块体地区的下地壳具有明显的S波低速层分布,该异常分布特征支持解释青藏高原隆升及其地壳物质运移的下地壳流模型.在100至120km深度上,川滇菱形块体西北部呈现较强的S波高速异常,这可能是印度岩石圈板块沿喜马拉雅东构造结下插至该区域所致,该区域下地壳的低速软弱物质与上地幔的高速强硬物质形成了鲜明对比,暗示了地壳和上地幔可能具有不同的构造运动和变形方式,这为该区域的壳幔动力学解耦提供了条件.     

Study on crust-mantle tectonics and its velocity structure along the Beijing-Huailai-Fengzhen profile

ｔｕｄｙｏｎｃｒｕｓｔｍａｎｔｌｅｔｅｃｔｏｎｉｃｓａｎｄｉｔｓｖｅｌｏｃｉｔｙｓｔｒｕｃｔｕｒｅａｌｏｎｇｔｈｅＢｅｉｊｉｎｇＨｕａｉｌａｉＦｅｎｇｚｈｅｎｐｒｏｆｉｌｅＺＨＩＰＩＮＧＺＨＵ（祝治平）ＸＩＡＮＫＡＮＧＺＨＡＮＧ（张先康）...     

苏鲁地区剪切波速度结构研究及与地震关系探讨

利用S波纯波形拟合法以及T函数法反演了苏鲁地区壳幔剪切波速度结构,并利用长周期P波T函数反演得到了连云港和莱阳台下方800km深度的速度结构。结果显示:(1)苏鲁地区大部分台站地壳表层及上地壳浅部速度偏高,分别对应高压、超高压物质和古老基底出露地区;(2)沿郯庐断裂带分布的台站均显示明显低速层,并具有北浅南深的特点;(3)连云港和莱阳台超深度反演结果显示两台均在150km深度下出现高速层,反映扬子板块的俯冲深度为100km以下,俯冲板片厚度在100km以上;板片拆离下沉深度甚至达到300km或者更深;(4)地震深度分布与低速层关系密切,沿郯庐断裂和烟台—五莲断裂的中小地震震源深度都比较深,有的甚至达到地壳的底部,反映这两条断裂目前切割深度都比较大,而且地幔物质相对比较活跃。