散射波积分方程的Adomian分解解法     

汪燚林  董良国 《地球物理学报》2021,64(10):3701-3717

在背景模型基础上,求解模型扰动后的地震波散射场,这是目前地震反演中的一个关键步骤.本文将计算数学中求解非线性积分方程的Adomian分解方法,应用到求解标量波散射场的Lippmann-Schwinger积分方程和Ricatti积分方程中,分别得到了散射场的Born序列解和Rytov序列解.通过一维和二维数值算例说明:在满足一定的条件下,散射场的这两种序列解稳定收敛,与传统的Born和Rytov近似解相比,引入散射序列中的高阶项可以更精确地描述地震波散射场.  相似文献   

基于Marchenko理论一步法压制层间多次波     

孙红日  王维红  石颖  徐嘉亮 《地球物理学报》2021,64(9):3295-3303

传统基于Marchenko理论自聚焦法压制层间多次波技术需要初始下行聚焦函数的估计,再进行多次波压制的自适应相减.论文基于传统方法的研究成果,提出了一种无需先验信息和自适应相减步骤的层间多次波压制方法.当聚焦点在实际反射层的反射点上时,上行格林函数为聚焦点上的脉冲源在表面所产生的上行波场,该反射层的反射将是上行格林函数的第一个同相轴.通过将格林函数与下行直达波场进行褶积,可将反射层处的所有焦点重新计算到地表,得到只含一次反射波的波场.利用初始下行聚焦函数和初始直达下行格林函数互为逆函数的关系,引入delta函数,将特定深度上所有聚焦点的聚焦函数投影到采集面上的检波器位置,可只利用地表获取的地震记录一步计算得到一次反射波.水平层状介质模型和复杂的含高速薄透镜地质体的介质模型中,进行了多次波压制的数值试验,通过与传统基于Marchenko理论自聚焦方法的压制效果进行对比分析,显示论文给出的一步法压制层间多次波的效果更好,进而验证了所提出方法的有效性和高精度的特点.  相似文献   

瞬变电磁扩散场到虚拟波场的精细积分变换算法   总被引：1，自引：0，他引：1  

鲁凯亮  李貅  戚志鹏  樊亚楠  周建美  李文翰  李贺  张明晶  王扬州 《地球物理学报》2021,64(9):3379-3390

现有的瞬变电磁解释与反演方法主要对电阻率参数进行解释,较难给出准确的地质构造信息.为了能够较多获得地质构造信息,可将瞬变电磁扩散场信号转换为虚拟波场信号,提高瞬变电磁法的分辨率.本文采用精细积分法,实现了瞬变电磁扩散场到虚拟波场的转换.将一个高度病态的线性方程组转换成求积分的过程,大大降低了解决病态问题的难度.积分步长以2的指数增加,在保证精度的前提下,具有极高的计算效率;设置合理的终止迭代条件,进一步提高本文方法的适用性.对典型地电模型的计算,表明本文方法具有较高的精度;减小两个波峰之间的距离,证明了本文方法具有较好的分辨能力;对含有噪声信号的扩散场进行波场反变换,证明了本文方法具有较好的抗噪性.最后计算三维模型和实测数据的虚拟波场,可以很好地反映出三维地质体的界面信息,证明本文方法稳定、可靠,可以对地下地质目标体的界面进行有效识别.本文方法相比预条件正则化共轭梯度法(PRCG),在计算效率方面提升了 4倍左右;本文方法的相对误差在5％以内,而PRCG的最大相对误差可达40％～50％,精度最高可提升10倍;在分辨率与稳定性方面,也表现出较好的效果.  相似文献   

巽他大陆及其邻区的地壳结构及其构造意义:来自远震接收函数的约束     

冯铭业  陈凌  王旭  韦生吉  王新 《地球物理学报》2021,64(12):4364-4377

巽他大陆位于欧亚板块、印度—澳大利亚板块和太平洋板块俯冲汇聚区域,其地壳结构特征是揭示洋陆过渡带演化及物质能量交换机制的重要依据.本文对巽他大陆及其周缘 19 个宽频带地震台站记录的远震波形进行 P波接收函数分析和H-κ叠加处理,获取了每个台站下方的地壳厚度和平均地壳波速比信息.为了减少参数的主观选择对结果带来的不确定性,研究采用了多种参数组合、综合约束策略.将本文结果与前人 146 个宽频带台站接收函数的研究结果进行整合,我们获得了巽他大陆地区地壳厚度和平均地壳波速比分布,并统计分析了两者的相关性.结果显示:巽他大陆地壳总体较薄,平均地壳厚度约为32 km,远低于全球造山带平均值,而与全球拉张型地壳平均厚度较为接近,可能反映研究区地壳整体处于拉张应力状态;而呵叻高原盆地地区地壳相对较厚,平均约38 km,与周缘地区明显不同.火山弧地区平均地壳波速比普遍大于 1.81,甚至达 1.87以上,并且壳内广泛分布低速层,可能受到了火山弧地区熔融物质的影响;非火山弧地区平均地壳波速比则普遍小于 1.76,反映地壳组分以长英质成分为主;局部地区高于 1.81,甚至高达 1.99,表明地壳以铁镁质成分为主或存在部分熔融,可能与铁镁质岩浆底侵作用或地幔热物质上涌有关.中南半岛中西部、婆罗洲西北部和马来半岛中部莫霍面 Ps转换波和多次波不明显而且具有多峰特征,可能表明该区域经历了复杂的壳幔相互作用.巽他大陆地区地壳厚度和平均地壳波速比总体无明显相关性,说明上地壳和下地壳结构和成分横向变化复杂;但中南半岛内部呵叻高原附近和东南部火山区两者均呈负相关性,与周围地区明显不同.综合区域构造背景和其他多种地球物理观测,推测稳定的呵叻高原盆地阻挡了印支地块的侧向挤出,处于挤压应力环境并发生上地壳增厚;而东南部火山区则处于拉张应力环境并存在基性岩浆底侵,可能与地幔物质上涌有关.  相似文献   

多尺度阵列嵌套组合反演宾川气枪源区横波速度结构   总被引：1，自引：0，他引：1  

孙楠  潘磊  王伟涛  叶泵  王彬  陈晓非 《地球物理学报》2021,64(11):4012-4021

本文基于多尺度阵列嵌套组合的方式,利用频率-贝塞尔变换法(Frequency-Bessel,F-J方法)提取背景噪声面波频散信息,通过多个阵列融合的频散曲线反演得到宾川气枪发射台周边不同深度的横波速度结构.结果显示:浅层一阶面波频散信息的加入,使得基阶反演结果更加收敛,反演深度加深到8 km;深度在8 km以下的结构的研究利用多尺度阵列(密集台阵—宾川气枪台网—云南区域地震台网)嵌套组合的方式,面波基阶低频信息从0.55 Hz拓宽到0.008 Hz,使横波速度结构的反演深度显著增加,同时对反演过程提供约束,使得70 km深度以上的横波速度更收敛.由此本文所得的横波速度结构为该区地下结构的探测提供基础,多尺度阵列嵌套组合频散谱的研究方式也为以后区域结构的研究提供一种新的方法和思路.  相似文献   

基于数据增广训练的深度神经网络方法压制地震多次波     

王坤喜  胡天跃  刘小舟  王尚旭  魏建新 《地球物理学报》2021,64(11):4196-4214

地震数据中存在的多次波影响偏移成像,误导地震资料的解释,因此通常视为相干噪声而被去除.为了对多次波进行智能化衰减,本文提出了一种基于数据增广训练的使用深度神经网络的多次波压制方法.设计的深度神经网络包括卷积编码和卷积解码过程,其中卷积编码过程学习全波场数据中的一次波特征,卷积解码过程利用这些特征来重构一次波并压制多次波和随机噪声.在训练阶段,旋转训练集并在输入数据中加入随机噪声构成增广训练数据集来提升神经网络的抗噪稳定性和泛化性,通过迁移学习让深度神经网络具备跨工区压制多次波的能力.简单模型与Sigsbee2B模型三套模拟数据的实例验证了本文方法在一次波重构和多次波压制中的有效性、稳定性和良好泛化性;一套崎岖海底模型地震物理模拟数据的应用实例表明本文方法具有应用于复杂条件下压制地震多次波的能力.  相似文献   

基于范阿伦双星EMFISIS观测数据的NWC和NAA人工甚低频台站信号的内磁层全球统计分布     

项正  林显浩  陈薇  王勇  陆鹏  龚文颖  马文琛  花漫  刘阳希子 《地球物理学报》2021,64(11):3860-3869

地球表面的人工甚低频台站信号可以穿透电离层泄漏进地球磁层导致内辐射带电子沉降到两极大气.因此研究人工甚低频台站信号的空间全球分布特性对于分析辐射带电子的损失具有重要科学意义.本文使用范阿伦双星从2013年到2018年共计6年的高质量的波动观测数据,统计了 NWC(19.8 kHz)、NAA(24.0 kHz)两个人工VLF台站信号的全球分布,分析了台站信号的电场功率谱密度对地理经纬度、磁壳值L、磁地方时MLT、地磁活动水平的依赖性.结果表明,在内磁层中,人工台站VLF信号主要沿着台站位置对应的磁力线传播,夜侧强度高于日侧,冬季高于夏季.这种日夜和夏冬差异的形成是因为夜侧和冬季的日照强度较弱,电离层电子密度较低,VLF信号较容易穿透电离层进入磁层.此外人工VLF台站信号的全球分布受地磁活动的影响很弱.这些统计观测结果给出了 NWC和NAA两个重要人工VLF台站信号强度的全球分布特征,为进一步分析人工VLF台站信号与地球辐射带电子的波粒相互作用提供了关键信息.  相似文献   

复波场的实部和虚部对同时源波形反演的影响     

杨育臣  方金伟  王宁  李松龄 《地球物理学报》2021,64(10):3718-3729

高精度的成像需要准确的速度场信息,波形反演被认为是目前具有最高分辨率的速度反演方法之一.计算效率是目前全波形反演需要考虑的一个主要问题.为了很好地解决计算效率的问题,本文引入了一种高效的无串扰同时源反演方法,并详细介绍了其原理与计算流程.同时,基于此多震源同时反演方法,本文拓展出实波场反演、虚波场反演及复波场反演的三种反演策略,进而分析复波场的实部和虚部对全波形反演的影响.实验表明,相比于复波场全波形反演,无论是实波场反演还是虚波场反演,反演分辨率有所降低;相比于实波场反演,虚波场反演对初始模型的依赖性较小,目标函数的非线性较弱;最后,通过使用组合反演策略,即初始阶段采用虚波场反演,中后期阶段采用复波场反演,不仅可以降低反演的非线性,而且能够保证高精度建模.  相似文献   

基于密集台阵地震背景噪声成像预测煤矿瓦斯分布     

黄宇奇  查华胜  高级  令狐建设  宣金国  周建斌  董润平  霍晶晶  张海江 《地球物理学报》2021,64(11):3997-4011

瓦斯突出是一种常见的煤矿动力灾害现象,随着煤矿矿井开采深度的增加,煤层瓦斯含量、压力都呈上升趋势,发生煤与瓦斯突出的危险性加大.传统的瓦斯测量方法只能测量局部离散点瓦斯含量,难以从矿井及采区尺度对瓦斯含量进行预测.因此需要寻求一种能够在采区及工作面布设前预测煤层瓦斯富集程度的高效地球物理方法.背景噪声成像方法已经在城市地下空间、矿产资源等近地表成像中得到广泛的运用.本文将该方法首次应用到阳泉寺家庄煤矿井田区域,采用96个台站记录的连续背景噪声数据,通过互相关方法获得了台站对之间的瑞利面波经验格林函数,并进一步提取了 5 Hz～1.4 s的基阶瑞利面波的群速度和相速度频散曲线.本研究首先通过区域的平均频散曲线获得该区域的平均一维横波速度结构作为三维反演的初始模型;其次,利用基于射线追踪的面波频散直接成像方法获得研究区1.0 km以浅的三维横波速度结构;最后,结合获得的三维速度结构,以及岩石物理实验获得的瓦斯含量与地震波速度的经验关系,预测了寺家庄井田15号煤的瓦斯含量,预测的瓦斯含量与实际巷道揭露的瓦斯含量具有较好的一致性.本研究成果表明,对于煤矿瓦斯分布预测来说,背景噪声成像方法是一种潜在有效的全新的技术形式.  相似文献   

喜马拉雅东构造结及周边地区地壳各向异性特征   总被引：1，自引：0，他引：1  

黄臣宇  常利军  丁志峰 《地球物理学报》2021,64(11):3970-3982

利用喜马拉雅东构造结及周边地区的48个宽频带地震台站记录的远震波形数据,通过Pms波分裂测量得到了 295对地壳各向异性横波分裂参数,获得了研究区的地壳各向异性图像.喜马拉雅东构造结的快波偏振方向主要为NE-SW方向,周边地区的快波偏振方向呈现出绕东构造结顺时针旋转的趋势.Pms波分裂的慢波延迟时间范围为0.11～0.30 s,平均值为0.24 s.对比分析研究区内Pms波分裂、近震直达S波分裂和远震SKS波分裂的结果发现,上地壳各向异性对Pms波分裂影响有限,地壳各向异性主要来自于中下地壳矿物和熔体的定向排列;地壳各向异性对SKS波分裂影响较小,SKS波分裂主要反映了上地幔的各向异性特征;Pms波分裂的快波偏振方向与近震直达S波分裂和远震SKS波分裂的快波偏振方向表现出较好的一致性,并且与地表构造和变形特征具有较好的相关性,反映了喜马拉雅东构造结及周边地区的岩石圈变形可能为垂直连贯变形模式.  相似文献