首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 140 毫秒
1.
天然气水合物似海底反射层的全波形反演   总被引:14,自引:4,他引:14       下载免费PDF全文
建立了天然气水合物似海底反射层(BSR)研究的全波形反演方法. 这是一种将 水平层状弹性介质的反射共中心点道集转换为截距时间-水平慢度域的反演方法. 反演过程 中采用了全局搜索方法与非线性局部搜索方法. 分两步进行. 第一步是根据走时数据应用非 常快速模拟算法求得速度结构的长波长分量. 第二步,利用波形资料用共轭梯度法求得速度 的短波长扰动分量. 这样,最后反演得到的速度结构模型包含了长波长与短波长分量. 反演 中利用了多网格参数化技术. 日本东南海海槽双BSR的速度结构的反演表明,全波形反演是 天然气水合物BSR研究的重要手段之一.  相似文献   

2.
变密度声波方程多参数全波形反演策略   总被引:4,自引:3,他引:1       下载免费PDF全文
多参数全波形反演中各参数之间的相互耦合增加了反演的非线性程度.通过分析各参数之间的相互影响,提出合理的多参数反演策略是解决该问题的有效途径.本文从变密度声波方程出发,首先研究了密度在速度反演中的重要作用,然后分析了速度对密度反演的影响程度,进而提出了一种有利于速度、密度分步联合反演的策略.第一步,利用给定的初始模型对速度、密度进行同时反演,得到比较可靠的速度反演结果;第二步,利用第一步反演得到的速度和给定的初始密度作为初始模型,继续进行双参数同时反演,这样可以同时得到比较可靠的速度、密度反演结果.为了进一步提高反演精度,将第二步反演得到的速度、密度作为初始模型,再进行下一轮双参数联合反演.二维理论模型实验结果充分说明了本文提出的这种反演策略的有效性.  相似文献   

3.
海洋勘探环境可以抽象为下伏固体与上覆流体相互耦合的介质,本文针对流-固边界耦合介质提出了一种高效、稳定的多参数(速度和密度)全波形反演方法.本文采用弹性波一阶位移-应力方程作为过渡层耦合声波压力方程与弹性波位移方程来模拟耦合环境,相比于传统的交错网格建模方法或者构建连续性条件,本文提出的方法在正演精度和稳定性上凸显出很大优势,极大降低了计算内存.反演策略对多参数全波形反演至关重要,由于不同参数之间的相互耦合使得密度在多参数全波形反演中较难获得,因此本文将非均匀流-固边界耦合介质多参数全波形反演分为两个步骤完成:第一步利用变密度声波方程结合推导出的密度梯度算子进行纵波速度和密度的双参数反演;第二步根据链式法则求取横波速度的梯度,结合第一步的反演结果使用流-固边界耦合方程反演横波速度.最后通过与声波动方程数值模拟结果对比证明正演算法的准确性;上覆流体的Marmousi-2模型的数值试验测试说明反演方法的有效性和适应性.  相似文献   

4.
反射地震数据的逐层波形反演   总被引:8,自引:2,他引:8       下载免费PDF全文
本文针对层状介质并结合梯度法波形反演,提出逐层波形反演的方法. 首先给出介质扰动响应的概念,并在此基础上分析了梯度法波形反演方法. 波形反演实质上是将实测地震记录和预测地震记录的波形残差信息转化为实际地质模型和预测地质模型的模型残差信息. 波形反演的优点是利用大量振幅相位信息得到高分辨率的反演结果, 其缺点是运行耗时大;当初始模型和实际模型相差较大时,迭代算法容易陷入局部极小点,这是因为目标函数和初始模型同实际模型间的差异是非线性的关系. 逐层波形反演方法是使自上而下每一层的目标函数最小,这样总的目标函数也是最小的. 利用二分法速度扫描确定每一层速度不仅提高了运算速度也避免了迭代算法陷入局部极小点的问题. 结合介质扰动响应和目标函数值变化可以更为准确迅速地确定每一层速度和该层界面位置.  相似文献   

5.
接收函数复谱比的最大或然性估计及非线性反演   总被引:32,自引:23,他引:32       下载免费PDF全文
利用Shamway和Der的多道最大或然性反褶积原理,讨论了从单台三分量远震P波波形中分离接收函数径向与垂向分量复谱比的方法.根据Tarantola的波形反演理论发展了接收函数复谱比的非线性反演方法.本文的反演方法除了拟合接收函数的复谱比之外,还需拟合时间域中接收函数垂向与径向分量的初至振幅比.合成波形数据反演的结果表明本文方法的反演结果不依赖初始模型.利用该方法和CDSN台网兰州台记录的宽频带远震P波波形数据,研究了该台站的接收函数及其随方位角的变化.接收函数复谱比非线性反演给出了该台站下方140km深度岩石层的S波速度结构,得到兰州台下方地壳厚度为54ks,岩石层厚度为94km.  相似文献   

6.
全球地震面波相速度变化及其大地构造学意义   总被引:6,自引:4,他引:6       下载免费PDF全文
从1980年到1992年间震级M≥6.0天然地震的全球数字地震台网记录,搜集约30000个长周期地震图波形数据,建立了周期为85-250s的全球Love波和Rayleish波相速度变化模型.全部地震图均通过时间域和频率域的质量控制,震相为短大圆弧和长大圆弧基频振型G1,R1,G2和R2.利用将区域高分辨率反演结果嵌入全球模型框架的混合反演参量化方法,以及地形、地貌和地壳厚度模型,得到校正前、后的两种反演结果.与已有全球相速度模型比较,长波分量的幅度和空间图像符合较好.短周期相速度异常反映了地表板块构造,与大陆地盾、大洋中脊和构造活动区相关;板块构造相关性随面波周期增大而逐渐消失;部分海洋热点与100-200s周期的低速异常相联系.  相似文献   

7.
得益于结合最优化算法和完整的波动方程,全波形反演已成为地震勘探最前沿的研究方向,正逐渐发展成为获取地球内部信息的重要工具.但是,在初始速度模型不够精准的条件下,利用纯反射波进行波形反演来更新光滑的背景速度模型是非常困难的.本文展示了一种利用纯反射波通过交替进行最小二乘逆时偏移和波形反演的方法来成功地同时反演反射系数和背景速度的方法.其波形反演中用于更新背景速度的层析梯度是通过基于Born近似的正演模拟来分离反射波波场得到的,这样可以很好地更新光滑的背景速度模型.之后,常规全波形反演把上一步得到的光滑背景速度模型作为初始模型,利用纯反射波获得最终的高保真、高分辨率的速度模型.进一步分析揭示了高分辨率的速度模型是由于在全波形反演中应用了反褶积成像条件.相对于其他类型的地震波,反射波带有地下更深部的信息,因此,反射波波形反演能够进一步提高我们对地球内部的成像精度.  相似文献   

8.
基于声波方程的井间地震数据快速WTW反演方法   总被引:4,自引:2,他引:4       下载免费PDF全文
WTW(Wave equation traveltime+Waveform inversion)反演是基于波动方程的走时反演(WT反演)和波形反演的联合反演方法.WT反演利用波动方程计算走时和走时关于速度的导数,和传统以射线为基础的走时反演相比,具有不必射线追踪、不必拾取初至、不必高频假设以及初始模型和实际模型差别较大时也能较好收敛等优点,但WT反演与波形反演相比其结果分辨率低.与之互补的是,波形反演的反演结果分辨率高,但是当所给初始模型和实际模型相差太大时,波形反演迭代算法容易陷入局部极小点.可见结合两种方法的WTW反演是一种比较好的联合反演方法.常规WTW迭代算法是首先以WT反演为主反演得到地质模型的整体特征,然后再以波形反演为主反演模型细节,该算法耗时和占用计算机存储空间接近WT反演或波形反演的两倍.为了节省运算耗时和计算机存储空间,往往采取首先单独利用WT反演然后再单独利用波形反演的算法.这样做的缺点是不能紧密结合两种反演方法,使得它们的优缺点在每一次迭代中无法得到互补,从而影响了最终的反演结果.针对以上事实,本文提出一种新的方法实现WTW,使得WTW运算速度和存储空间在任何情况下等同于WT反演或波形反演.模型计算表明新的算法具有更好的收敛性.  相似文献   

9.
电磁测深方法研究   总被引:18,自引:0,他引:18       下载免费PDF全文
为能够实施电磁测深新方法的基本原理,并为研制探测仪器提供一些依据,分析和讨论了源信号为尖脉冲时,入射波与地层反射波通过二阶带通滤波器的输出波形;利用反射波与入射波前沿瞬态叠加的波形,反演地层的电导率和厚度的方法及与探测仪器有关的主要问题等内容. 结果表明,反演仍能给出惟一解,并具有可实践性.  相似文献   

10.
逆时偏移对速度模型有较强的依赖性,速度模型精度高低决定着偏移效果的好坏,速度模型越准确,偏移效果越好.因此,如何获取高精度的速度模型是实现有效逆时偏移的关键步骤.全波形反演方法利用叠前地震波场的运动学和动力学信息重建地下速度结构能够得到高精度的速度模型.为了验证全波形反演能否提高逆时偏移效果,本论文以OVERTHRUST速度模型为例,在频率域进行全波形反演,并实现了速度模型的更新.逐步缩小了初始速度模型与真实速度模型的残差,提高了速度模型的准确度.然后使用初始OVERTHRUST速度模型以及经过反演的速度模型分别进行逆时偏移,应用结果表明,使用经过全波形反演后的速度模型进行逆时偏移,其成像效果较好.因此,逆时偏移结合全波形反演可以有效的提高成像精度.  相似文献   

11.
利用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)首都圈内大地震的分布与壳内低速体及上地幔顶部的速度结构有密切关系.对于唐山大地震的成因,仅考虑板块作用引起的水平应力场是不够的,有必要充分重视由于上地幔变形引起的地壳垂直变形和上地幔物质侵入造成的热效应.  相似文献   

12.
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.  相似文献   

13.
在沙城以东的延庆盆地及其邻近区域内布设了由GDS-1000宽频带数字地震仪组成的流动地震台阵,利用台阵记录的宽频带远震P波波形数据和非线性接收函数反演方法获得了延怀盆地内0-80km深度范围的地壳、上地幔S波速度结构.利用计算机三维彩色剖分显示技术研究了台阵下地壳、上地幔速度结构的横向非均匀变化。结果表明,研究区域内的地壳厚度为40km左右,壳幔界面有4km左右的上下起伏.地表沉积盖层在延庆盆地中心附近厚度约1km,而在向盆地外围延伸的方向上相对变薄.研究区域内上地壳S波速度结构较复杂,而下地壳与上地幔则相对比较均匀.其上地壳最突出的特点是在10km深度附近有明显的S波低速层.在延庆盆地下方,它延伸到6-20km的深度范围.在延庆盆地南侧,该低速层有从西往东逐渐减弱的趋势.研究区域内的地震基本上都发生在延庆盆地下方上地壳低速体外围.据此推断,延庆盆地及其临近区域内的地震活动与该区域地壳内的热状态有密切关系.  相似文献   

14.
甘肃东南部地壳速度结构的区域地震波形反演   总被引:6,自引:1,他引:5       下载免费PDF全文
利用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模型给出的全球平均地壳速度值相比,两个分区的地壳速度值整体偏低.  相似文献   

15.
用转换函数方法研究腾冲—临沧地区地壳结构   总被引:4,自引:1,他引:3       下载免费PDF全文
根据流动数字地震台网提供的三分量地震波形记录资料,应用转换函数及快速模拟退火算法对腾冲-临沧地区30个地震台站下的地壳横波速度结构进行了反演.反演结果说明,研究区壳幔边界清晰、莫霍界面附近速度跳跃明显,由此得出该区地壳厚度在40 km左右、并具有从南向北增厚趋势.一个普遍的现象是,在腾冲-宝山地块下地壳存在明显的低速带,低速带的厚度在10~20 km间.研究结果进一步表明各台站下方上地幔速度结构复杂.这些结果为探讨青藏高原东南缘下地壳的侧向黏性流动、碰撞板块边界处壳幔物质交换等均提供了重要的地球物理证据,为探讨印-藏汇聚过程中青藏高原东构造结岩石圈变形、高原隆升及其深部动力学有一定理论意义.  相似文献   

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

17.
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.  相似文献   

18.
喜马拉雅山北部地区的地壳结构模型和速度分布特征   总被引:16,自引:8,他引:16       下载免费PDF全文
本文根据1981年西藏南部喜马拉雅地区的人工地震测深资料进行了震相对比,分辨出t1、t2、t3、t4、t5和t6六组地壳中和莫霍界面的反射波,并用理论走时曲线、绘制速度曲线图、射线跟踪和综合地震图等方法得到了主测线(PP)上各地段的地壳结构模型。初步结果表明,该地区地壳西段较薄(约73公里),东段稍厚(约77公里),平均总厚度约为75公里。地壳的平均P波速度约为6.2-6.3公里/秒。 地壳为高低速相间的多层结构。在中上部有一低速层,其厚度为数公里,速度为5.6-5.7公里/秒,与上层速度差为0.5-0.6公里/秒。低速层在测线东段比较肯定,在西段则不甚明显。结合藏南定日、岗巴一线有强烈水热活动的事实,低速层的存在可能意味着地壳中存在部分熔融的高温物质。下部地壳的速度为6.7-6.8公里/秒,且比较均匀。从莫霍面反射波的特征来看,在紧靠其上方可能有一个速度反转带,其厚度亦为数公里。上部地壳的结构在横向上有较大的差异,这说明在地质历史上,西藏特提斯带曾经历过强烈的地壳变动。  相似文献   

19.
Using the seismic records of 83 temporary and 17 permanent broadband seismic stations deployed in Tangshan earthquake region and its adjacent areas (39°N–41.5°N, 115.5°E–119.5°E), we conducted a nonlinear joint inversion of receiver functions and surface wave dispersion. We obtained some detailed information about the Tangshan earthquake region and its adjacent areas, including sedimentary thickness, Moho depth, and crustal and upper mantle S-wave velocity. Meanwhile, we also obtained the vP/vS structure along two sections across the Tangshan region. The results show that: (1) the Moho depth ranges from 30 km to 38 km, and it becomes shallower from Yanshan uplift area to North China basin; (2) the thickness of sedimentary layer ranges from 0 km to 3 km, and it thickens from Yanshan uplift region to North China basin; (3) the S-wave velocity structure shows that the velocity distribution of the upper crust has obvious correlation with the surface geological structure, while the velocity characteristics of the middle and lower crust are opposite to that of the upper crust. Compared with the upper crust, the heterogeneity of the middle and lower crust is more obvious; (4) the discontinuity of Moho on the two sides of Tangshan fault suggests that Tangshan fault cut the whole crust, and the low vS and high vP/vS beneath the Tangshan earthquake region may reflect the invasion of mantle thermal material through Tangshan fault.  相似文献   

20.
新疆及其毗邻地区地震层析成像   总被引:18,自引:5,他引:18       下载免费PDF全文
利用新疆地震台网及其毗邻地区的国外地震台站提供的地震P波到时数据,重建了新疆及其毗邻地区地壳上地幔的三维速度图像.主要结果是:1.研究区地壳上地幔的地震波速度存在明显的横向不均匀性.这种不均匀性出现在不同构造单元之间,也出现在各构造单元之内.2.一般来说,研究区域内盆地是上地幔的隆起区,褶皱系是上地幔的坳陷区.3.岩石层厚度和软流层的深度与全球的平均结果有较大不同.4.海洋岩石层或大陆岩石层板块在兴都库什和克什米尔与中国交界地区的消减明显存在.5.地壳上地幔物质运动比较剧烈,有部分熔融物质上涌到地壳中部或更浅的深度.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号