球坐标系下多震相走时三参数同时反演成像 被引量：5

1

作者 黄国娇 白超英 钱卫 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2015年第10期3627-3638,共12页

球坐标系下多震相走时三参数(速度、震源位置和反射界面)同时反演需要解决两个关键问题:(1)球坐标系下3D速度模型中多次透射、反射(折射)及转换波精确、快速的射线追踪;(2)同时反演时三种不同参数间的强耦合问题.为此,我们将直角坐标系... 球坐标系下多震相走时三参数(速度、震源位置和反射界面)同时反演需要解决两个关键问题:(1)球坐标系下3D速度模型中多次透射、反射(折射)及转换波精确、快速的射线追踪;(2)同时反演时三种不同参数间的强耦合问题.为此,我们将直角坐标系下分区多步不规则最短路径算法推广至球坐标系中,进行区域或者全球尺度的多震相射线追踪.然后将其与适合多参数同时反演的子空间算法相结合,形成一种球坐标系下联合多震相走时三参数同时反演的方法技术.与双参数(速度和反射界面或速度和震源位置)同时反演的数值模拟对比分析显示:三参数与双参数的同时反演结果大体接近,并且它们对到时数据中可容许的随机噪声不太敏感.结果说明本文中的同时反演成像为一种提高成像分辨率,同时反演速度、震源位置和反射界面的有效方法. 展开更多

关键词 球坐标系 分区多步不规则最短路径算法 多震相走时 三参数同时反演 子空间法

在线阅读 下载PDF 职称材料

椭球坐标系下多震相地震射线追踪 被引量：2

2

作者 李兴旺 白超英 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2017年第9期3368-3377,共10页

地震射线追踪方法技术在地震学领域有着较为广泛的应用,然而大多数算法建立在直角坐标系或球坐标系下,实际地球并非完美的球体,而是两极略扁的椭球体,因此,球坐标系下计算结果与真实情况存在一定误差.传统的做法一般是在球坐标系下进行... 地震射线追踪方法技术在地震学领域有着较为广泛的应用,然而大多数算法建立在直角坐标系或球坐标系下,实际地球并非完美的球体,而是两极略扁的椭球体,因此,球坐标系下计算结果与真实情况存在一定误差.传统的做法一般是在球坐标系下进行计算,而后进行椭球校正.本文提出了一种直接在椭球体模型中采用分区多步最短路径算法进行多震相地震射线追踪的方法技术,实现了椭球坐标系下多震相地震波射线路径追踪和走时计算.与解析解的对比表明:该算法具有较高的计算精度,适用于任意形状的椭球体,且不需要进行额外的走时校正.数值模拟结果表明,计算所得P波和PcP反射波的走时与AK135走时表的误差小于0.1 s.当震中距较大时,使用球对称模型和椭球体模型计算所得的走时差异显著,说明采用椭球坐标系的必要性. 展开更多

关键词 椭球坐标系 最短路径算法 分区多步 多震相射线追踪

在线阅读 下载PDF 职称材料

一般各向异性介质中qS1和qS2波波前分离

3

作者 李兴旺 王迪 白超英 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2024年第7期2773-2784,共12页

地震波在各向异性介质中传播时存在横波分裂现象,可用于研究岩石中地震波传播特征和了解地下介质各向异性参数.一般各向异性介质中计算三种地震波(qP,qS1,qS2)相速度和群速度的研究日趋成熟,但分离qS1和qS2波的波前仍然是一项具有挑战... 地震波在各向异性介质中传播时存在横波分裂现象,可用于研究岩石中地震波传播特征和了解地下介质各向异性参数.一般各向异性介质中计算三种地震波(qP,qS1,qS2)相速度和群速度的研究日趋成熟,但分离qS1和qS2波的波前仍然是一项具有挑战性的工作.为此,本文基于新近给出的一般各向异性介质中三种波(qP,qS1,qS2)群速度的解析计算公式,依据相速度光滑连续的条件给出相邻方向上qS1和qS2波识别方法;为避免数值计算误差造成识别错误,构建了置信度函数,并借鉴快速行进法的思想逐点进行判别以分离qS1和qS2波的波前,最终得到三种波(qP,qS1,qS2)各自独立且完整的波前面.多个数值模拟试验结果表明,该算法适用于不同类型各向异性介质,且不受波前复杂程度的影响,具有较高的稳定性,可为一般各向异性介质中多次波射线追踪、多震相走时反演以及横波分裂研究等提供算法支持. 展开更多

关键词 一般各向异性介质 群速度 横波分裂 波前分离

在线阅读 下载PDF 职称材料

三维TI和OA介质中优化黏滞纯声波方程及其数值模拟

4

作者 包乾宗 谭崔文 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2024年第6期2415-2428,共14页

各向异性和黏滞性在地球介质中广泛存在,会对地震波的运动学和动力学特征产生影响,因此研究地震波在地下空间中传播时应该对其予以充分考虑.现有的各向异性波场外推方法主要采用基于声学近似的伪声波方程,该类方程在非声学假设条件下容... 各向异性和黏滞性在地球介质中广泛存在,会对地震波的运动学和动力学特征产生影响,因此研究地震波在地下空间中传播时应该对其予以充分考虑.现有的各向异性波场外推方法主要采用基于声学近似的伪声波方程,该类方程在非声学假设条件下容易出现伪横波污染及传播不稳定.针对此问题,本文应用优化方案对原始的耦合各向异性纯声波频散关系式进行近似展开,进而发展了针对三维TI(横向各向同性)介质和OA(正交各向异性)介质的优化纯声波方程,结合有限差分和泊松算法对其进行了有效求解.同时,考虑到介质对地震波的吸收衰减效应,本文结合发展的纯声波方程和标准线性体模型,推导了黏滞各向异性优化纯声波方程.相速度分析表明优化方案能够产生高精度近似结果.数值模型算例证明优化黏滞纯声波方程在产生精确稳定P波响应的前提下,能够有效模拟各向异性介质中纯声波的振幅衰减和相位频散.整体而言,本文推导的方程能够为各向异性衰减介质成像和反演提供有效帮助. 展开更多

关键词 各向异性 黏滞性 波动方程 数值模拟

在线阅读 下载PDF 职称材料

应用遗传算法估计非稳态地震数据的混合相位子波及Q值 被引量：3

5

作者 朱耀旭 包乾宗 《石油地球物理勘探》 EI CSCD 北大核心 2024年第4期763-770,共8页

地震高分辨率处理或反演的目的是获得精确的反射系数或弹性参数模型。然而地层滤波效应模糊了地震记录中的地层反射信息,因此有必要消除这种滤波效应。现有的大部分提高分辨率方法并不能完全摒弃关于地震子波和Q模型的一些假设。为了获... 地震高分辨率处理或反演的目的是获得精确的反射系数或弹性参数模型。然而地层滤波效应模糊了地震记录中的地层反射信息,因此有必要消除这种滤波效应。现有的大部分提高分辨率方法并不能完全摒弃关于地震子波和Q模型的一些假设。为了获得更为切合实际的地震子波,同时自适应获取Q模型,文中将地震子波相位估计与Q模型估计相结合,提出了基于遗传算法的非稳态地震数据混合相位子波及Q值估计方法。首先通过井旁道记录拟合得到初始地震子波的振幅信息,然后依据子波Z变换的根关于单位圆移动与否构建用于遗传算法的编码链条。另一方面,十进制Q模型对应的二进制表示形式同样能够利用编码链条表征,因此利用全局优化算法能够同时估计地震混合相位子波以及Q模型。结合根变换和遗传算法不断改变子波相位的同时自适应生成Q模型,利用子波和Q模型构造的时变子波矩阵与测井反射系数得到合成地震记录,并与井旁道记录进行匹配,最终得到合理的混合相位子波和地层Q模型,进而构造时变子波矩阵进行时变反褶积。通过井旁道记录与测井数据拟合得到的混合相位子波,与实际地震子波相位更为接近,理论数据和实际数据处理结果证实了该方法的有效性。 展开更多

关键词 高分辨率处理 品质因子 非稳态地震数据 混合相位子波 遗传算法

在线阅读 下载PDF 职称材料

矩形/长方体网格单元时间高阶有限差分波场数值模拟 被引量：1

6

作者 徐世刚 黄兴国 +2 位作者 韩丽 包乾宗 任志明 《石油地球物理勘探》 EI CSCD 北大核心 2024年第5期976-988,共13页

有限差分法在求解不同类型地震波动方程中广泛应用,提高有限差分在时间域和空间域的地震波场数值模拟精度十分重要。为了提高现有时间—空间域高阶差分方法的精度和灵活性,提出了基于矩形/长方体网格单元的改进时间—空间高阶有限差分方... 有限差分法在求解不同类型地震波动方程中广泛应用,提高有限差分在时间域和空间域的地震波场数值模拟精度十分重要。为了提高现有时间—空间域高阶差分方法的精度和灵活性,提出了基于矩形/长方体网格单元的改进时间—空间高阶有限差分方案,分别求解二维和三维声波方程。改进差分法主要通过联合旁轴网格节点和传统轴向节点来共同完成偏导数近似,进一步结合平面波理论和泰勒级数展开求取改进模板中的高阶差分系数。数值精度分析表明,改进差分方案相较于传统方法具有更高的数值精度和稳定性。多个模型算例表明,在相同参数条件下,改进差分方法在兼顾空间模拟精度的前提下能够提高时间模拟精度。此外,改进模板沿不同空间方向可以选用不同的网格步长,有效增加了波场模拟的灵活性,能够为后续的地震成像和反演提供有效波场延拓工具。 展开更多

关键词 地震波动方程 正演模拟 高阶有限差分 矩形/长方体网格 差分系数

在线阅读 下载PDF 职称材料

基于改进交错网格差分结构的地震波场正演模拟

7

作者 徐世刚 包乾宗 任志明 《桂林理工大学学报》 CAS 北大核心 2024年第2期225-231,共7页

有限差分离散在地震波场数值模拟中广泛使用,相比于传统时间二阶差分,时间高阶差分具有更高的精度和稳定性,隐式差分比显式差分具有更高的精度。本文结合时间高阶方法和隐式差分的优势,提出了一种改进的交错网格差分结构求解三维一阶声... 有限差分离散在地震波场数值模拟中广泛使用,相比于传统时间二阶差分,时间高阶差分具有更高的精度和稳定性,隐式差分比显式差分具有更高的精度。本文结合时间高阶方法和隐式差分的优势,提出了一种改进的交错网格差分结构求解三维一阶声波方程。在该方案中,时间高阶离散被用于补偿传统时间二阶差分的近似精度,同时采用隐式交错网格差分求解空间导数;给出了关于时间离散格式的泰勒展开差分系数,求取了空间隐式离散格式中的泰勒展开和最小二乘两种差分系数。与现有的几种交错网格差分方案相比,改进方案能够在空间不同方向上采用不同的网格间距和差分算子长度,因此在保证模拟精度的前提下灵活性更高。数值分析和模型算例表明,改进方案能够采用较短的算子长度获得高精度模拟结果。 展开更多

关键词 数值模拟 波动方程 有限差分 矩形交错网格

在线阅读 下载PDF 职称材料

三维复杂层状介质中多震相走时联合反演成像 被引量：21

8

作者 白超英 黄国娇 李忠生 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2011年第1期182-192,共11页

采用新近提出的改进型不规则最短路径多次波射线追踪正演技术,结合共轭梯度法求解带约束的阻尼最小二乘最优化反演问题,讨论了三维复杂层状模型中利用多震相走时资料进行联合反演成像的技术方法.考虑到不同震相种类走时的拾取误差不同,... 采用新近提出的改进型不规则最短路径多次波射线追踪正演技术,结合共轭梯度法求解带约束的阻尼最小二乘最优化反演问题,讨论了三维复杂层状模型中利用多震相走时资料进行联合反演成像的技术方法.考虑到不同震相种类走时的拾取误差不同,反演算法中引入了不同震相种类数据的权系数;另外,考虑到同时反演速度模型和反射界面起伏中不同参数变化对走时影响程度的不同,Jacobi偏导矩阵元素中引入了不同种类参数的归一化因子.几种数值模拟实例表明:多震相走时的联合或同时反演成像是一种提高走时成像空间分辨率,进而降低重建速度模型失真度行之有效的方法技术. 展开更多

关键词 改进型不规则最短路径 多次波射线追踪 多震相走时 联合反演 同时反演 地震层析成像

在线阅读 下载PDF 职称材料

基于连续小波变换的自适应面波压制方法 被引量：35

9

作者 陈文超 高静怀 包乾宗 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2009年第11期2854-2861,共8页

面波干扰的抑制是陆上地震资料处理的主要问题之一.本文根据炮集记录中面波与反射波主要能量在小波域分布区域的不同,以及面波干扰的影响随炮检距变化等特点,提出了一种具有时变、空变特性的自适应面波衰减方法.文中将该方法用于模型及... 面波干扰的抑制是陆上地震资料处理的主要问题之一.本文根据炮集记录中面波与反射波主要能量在小波域分布区域的不同,以及面波干扰的影响随炮检距变化等特点,提出了一种具有时变、空变特性的自适应面波衰减方法.文中将该方法用于模型及实际炮集资料的处理,并与常用的高通滤波方法进行对比,结果表明,该方法在衰减面波干扰的同时,能更好地保持反射波的振幅及相位信息. 展开更多

关键词 连续小波变换 面波压制 自适应滤波

在线阅读 下载PDF 职称材料

二维复杂层状介质中地震多波走时联合反演成像 被引量：35

10

作者 黄国娇 白超英 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2010年第12期2972-2981,共10页

采用新近提出的多次波射线追踪正演算法,结合共轭梯度法求解带约束的阻尼最小二乘最优化反演问题,分析讨论了利用多震相走时资料进行联合反演成像的方法及技术.考虑到不同震相走时的拾取误差不同,反演算法中引入了不同震相种类数据的权... 采用新近提出的多次波射线追踪正演算法,结合共轭梯度法求解带约束的阻尼最小二乘最优化反演问题,分析讨论了利用多震相走时资料进行联合反演成像的方法及技术.考虑到不同震相走时的拾取误差不同,反演算法中引入了不同震相种类数据的权系数;由于同时反演速度模型和反射界面起伏中不同模型参数变化对走时影响程度的不同,Jacobi偏导矩阵元素中引入了不同参数的归一化因子;另外,为了克服射线密度过大(或过小)区域速度模型的过度(或欠)更新问题,反演算法中引入了等权射线密度的概念.几种数值模拟实例表明(含噪声敏感性试验):多波走时的联合或同时反演成像技术是一种提高走时成像空间分辨率,进而降低重建模型失真度行之有效的方法. 展开更多

关键词 多次波射线追踪 多震相走时 联合反演 同时反演 等权射线密度 地震层析成像 噪声敏感性试验

在线阅读 下载PDF 职称材料

最短路径算法下三维层状介质中多次波追踪 被引量：27

11

作者 唐小平 白超英 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2009年第10期2635-2643,共9页

本文使用改进后的最短路径算法(MSPM)结合分区多步计算技术实现了三维复杂层状起伏介质中的多次透射、反射及转换波波前传播的数值模拟,以及相应走时和射线路径的跟踪计算.其原理是将三维复杂层状模型按速度界面分成若干个独立的计算区... 本文使用改进后的最短路径算法(MSPM)结合分区多步计算技术实现了三维复杂层状起伏介质中的多次透射、反射及转换波波前传播的数值模拟,以及相应走时和射线路径的跟踪计算.其原理是将三维复杂层状模型按速度界面分成若干个独立的计算区域,采用分步计算技术进行多次波的跟踪计算.基于多次波是通过速度界而简单的人射、透射、反射及转换波按一定规律及原理的不同组合,因此可实施分区多步计算技术.数值模拟实例及误差分析表明分区多步计算技术具有单步最短路径算法中的诸多优点:算法简单、数值计算稳健、计算精度高、速度快及全球解等,因此是解决多次波跟踪计算行之有效的方法. 展开更多

关键词 MSPM算法 分区多步计算 多波波前传播模拟 多次波追踪

在线阅读 下载PDF 职称材料

一般TI介质中多震相旅行时层析成像——以井间成像为例 被引量：9

12

作者 黄国娇 白超英 钱卫 《石油地球物理勘探》 EI CSCD 北大核心 2016年第1期115-126,21,共12页

传统的各向异性介质中的正、反演研究大多基于弱各向异性假设下的VTI介质,且主要是有关初至波的射线追踪和相应的旅行时反演。然而,随着勘探精度要求的不断提高,一般普适的(含强)TI介质中的多震相射线正、反演方法,则成为急需解决的问... 传统的各向异性介质中的正、反演研究大多基于弱各向异性假设下的VTI介质,且主要是有关初至波的射线追踪和相应的旅行时反演。然而,随着勘探精度要求的不断提高,一般普适的(含强)TI介质中的多震相射线正、反演方法,则成为急需解决的问题。鉴于此,本文将基于各向同性介质多震相分区多步最短路径射线追踪算法推广到一般各向异性TI介质,结合相、群速度导数的一阶旅行时扰动方程,采用共轭梯度求解带约束的阻尼最小二乘问题,进而提出了一种利用多震相旅行时进行一般TI介质弹性参数反演的方法。井间地震数值实验结果表明,多震相旅行时联合反演可有效提高成像分辨率,验证了方法的有效性和正确性。 展开更多

关键词 各向异性TI介质 多震相射线追踪 弹性参数反演 井间旅行时层析

在线阅读 下载PDF 职称材料

复杂层状模型中多次波快速追踪——一种基于非规则网格的最短路径算法 被引量：21

13

作者 赵瑞 白超英 《地震学报》 CSCD 北大核心 2010年第4期433-444,共12页

使用不规则网格单元划分下的最短路径算法,结合分区多步计算技术实现了二维和三维复杂层状起伏介质中的多次透射、反射及转换波的追踪计算.其原理是将模型按速度界面分成若干个独立的计算区域(在速度界面和起伏地表处采用一种不规则网... 使用不规则网格单元划分下的最短路径算法,结合分区多步计算技术实现了二维和三维复杂层状起伏介质中的多次透射、反射及转换波的追踪计算.其原理是将模型按速度界面分成若干个独立的计算区域(在速度界面和起伏地表处采用一种不规则网格单元划分),采用分步计算技术进行多次波的追踪计算.通过与有限差分下快速行进法的比较,表明无论是计算精度还是CPU时间,不规则最短路径算法均好于快速行进法算法.最后,实例模拟中给出了二维、三维复杂层状模型(包括Marmousi模型及含低速体的模型)中的多次透射、反射及转换波的追踪计算,验证了不规则最短路径算法的功能. 展开更多

关键词 不规则最短路径算法 分区多步计算 多次波追踪 快速行进法

在线阅读 下载PDF 职称材料

多震相走时联合三参数同时反演成像 被引量：9

14

作者 黄国娇 白超英 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2013年第12期4215-4225,共11页

采用新近研制的分区多步不规则最短路径多震相地震射线追踪正演技术,结合流行的子空间反演算法,提出了一种联合多震相走时资料进行地震三参数(速度、反射界面和震源位置)同时反演的方法技术.数值模拟反演实例、以及与双参数(速度和反射... 采用新近研制的分区多步不规则最短路径多震相地震射线追踪正演技术,结合流行的子空间反演算法,提出了一种联合多震相走时资料进行地震三参数(速度、反射界面和震源位置)同时反演的方法技术.数值模拟反演实例、以及与双参数(速度和反射界面或速度和震源位置)同时反演的对比分析表明:三参数同时反演成像结果大体接近双参数同时反演成像的结果.另外,噪声敏感性试验表明:所提算法对到时数据中可容许的随机误差并不敏感,结果说明多震相走时的联合三参数同时反演成像方法技术不失为一种提高走时成像空间分辨率、进而降低重建模型参数失真度、行之有效的方法技术. 展开更多

关键词 分区多步不规则最短路径算法 多震相射线追踪 多震相走时 三参数同时反演 地震层析成像

在线阅读 下载PDF 职称材料

起伏层状TI介质中多次波射线追踪 被引量：12

15

作者 李晓玲 白超英 胡光义 《石油地球物理勘探》 EI CSCD 北大核心 2013年第6期924-931,1016+850,共8页

现有的TI介质中的射线追踪算法很少涉及多次波的追踪计算,加之大多数算法又是建立在规则网格(或单元)模型划分的基础上,不利于追踪不规则界面下的多次波射线路径和相应走时。本文在模型参数化时引入不规则网格(或单元)划分,在主节点上... 现有的TI介质中的射线追踪算法很少涉及多次波的追踪计算,加之大多数算法又是建立在规则网格(或单元)模型划分的基础上,不利于追踪不规则界面下的多次波射线路径和相应走时。本文在模型参数化时引入不规则网格(或单元)划分,在主节点上对TI介质的6个参数(五个弹性参数,外加介质对称轴倾角)进行采样,计算得到三种类型(qP,qSV,qSH)波的群速度,通过单元主节点的群速度值线性插值得到次级节点上的群速度值。然后采用网格(或单元)波前扩展的方法进行分区波前扫描(射线追踪),实现了起伏层状TI介质中多次(透射、反射、及转换)波的追踪计算。均匀各向异性介质中的计算结果表明,数值解与解析解吻合,验证了算法的有效性;模型数值模拟结果表明,该算法可以有效地模拟起伏地表和速度间断面,具有较高的计算精度和较快的CPU运行时间,能够适用于任意倾角的TI介质中的多次波射线追踪。 展开更多

关键词 TI介质 最短路径算法 多次波射线追踪 不规则网格 分区多步计算技术

在线阅读 下载PDF 职称材料

四面体单元剖分下三维各向异性TI介质中多次波射线追踪 被引量：5

16

作者 李兴旺 白超英 李晓玲 《石油地球物理勘探》 EI CSCD 北大核心 2017年第1期48-55,共8页

采用不规则四面体单元做模型剖分,实现了三维复杂层状各向异性TI介质中多次(透射、反射及转换)波的追踪计算,弥补了现行算法主要采用规则网格计算多次波的不足。根据介质的5个弹性参数及其对称轴倾角和走向,计算主节点上qP、qSV和qSH波... 采用不规则四面体单元做模型剖分,实现了三维复杂层状各向异性TI介质中多次(透射、反射及转换)波的追踪计算,弥补了现行算法主要采用规则网格计算多次波的不足。根据介质的5个弹性参数及其对称轴倾角和走向,计算主节点上qP、qSV和qSH波群速度,由此线性插值得到次级节点的群速度值;然后采用不规则最短路径算法,并结合分区多步计算技术实现多次波追踪计算。均匀各向异性介质中计算出的数值解与解析解吻合,验证了算法的有效性。复杂层状模型数值模拟结果表明:采用不规则四面体单元可有效拟合起伏地表及地下不规则速度间断面,同时该算法适用于任意倾角的TI介质,且不受各向异性强度限制。 展开更多

关键词 三维各向异性TI介质 改进型最短路径算法 多次波射线追踪 四面体单元 分区多步计算技术

在线阅读 下载PDF 职称材料

高阶交错网格有限差分弹性波场模拟的精度分析 被引量：2

17

作者 岳晓鹏 白超英 岳崇旺 《煤田地质与勘探》 CAS CSCD 北大核心 2017年第1期125-130,共6页

交错网格波场数值模拟是目前地震正演中广泛使用的方法,为对比分析不同阶数的差分格式下产生的计算效率和精度差异,重新推导了弹性波方程的4种时间4阶、空间2N阶的差分公式及系数,并计算了他们的稳定性条件。利用这4种差分格式进行弹性... 交错网格波场数值模拟是目前地震正演中广泛使用的方法,为对比分析不同阶数的差分格式下产生的计算效率和精度差异,重新推导了弹性波方程的4种时间4阶、空间2N阶的差分公式及系数,并计算了他们的稳定性条件。利用这4种差分格式进行弹性波场数值模拟,对比分析了波场快照、合成地震记录及CPU时间。结果表明:时间4阶、空间6+6阶精度的交错网格有限差分方法在进行地震波场数值模拟时具有较高的计算精度和计算效率。 展开更多

关键词 精度分析 弹性波场 交错网格 高阶有限差分

在线阅读 下载PDF 职称材料

一种基于新差分模板和无穷范数的震源波场重建方法 被引量：2

18

作者 包乾宗 戴雪 梁雪 《石油地球物理勘探》 EI CSCD 北大核心 2022年第6期1384-1394,I0005,共12页

伴随状态法广泛应用于偏移成像和全波形反演中。模型的像或梯度可通过震源波场和伴随波场之间的相互运算得到。但两种波场分别沿时间正向和反向传播,无法同时访问。为避免对震源波场进行存储,可采用边界波场逆向重建震源波场。但现有的... 伴随状态法广泛应用于偏移成像和全波形反演中。模型的像或梯度可通过震源波场和伴随波场之间的相互运算得到。但两种波场分别沿时间正向和反向传播,无法同时访问。为避免对震源波场进行存储,可采用边界波场逆向重建震源波场。但现有的重建方法仍无法平衡精度和存储需求。为此,发展了一种新的交错网格有限差分震源波场重建方法。新方法通过存储边界区域的N层波场和M-N层波场的线性组合来重建内部区域的震源波场(M为差分算子长度参数,0≤N≤M)。推导了基于新差分模板的频散关系,构建了无穷范数型目标函数,采用Remez交换算法优化了重建系数。详细分析了所提方法的精度和稳定性,并将其应用于声波逆时偏移和弹性波全波形反演。数值结果表明:新方法可以获得足够精确的重建波场、偏移剖面和反演结果;其内存需求仅为传统方法的(N+1)/M。 展开更多

关键词 震源波场重建 有限差分 交错网格 逆时偏移 全波形反演

在线阅读 下载PDF 职称材料

三维复杂速度模型中地震的快速精确定位

19

作者 白超英 赵瑞 李忠生 《地震学报》 CSCD 北大核心 2009年第4期385-395,477,共11页

讨论了全局选择震源初始位置下的矩阵反演求取全局解的问题.与流行的全局优化搜寻直接定位方法相比,该反演算法是一种行之有效的定位方法,具有定位精度高、计算时间少以及对噪声数据不敏感等特点.其突出的优点是在不增加计算难度和计算... 讨论了全局选择震源初始位置下的矩阵反演求取全局解的问题.与流行的全局优化搜寻直接定位方法相比,该反演算法是一种行之有效的定位方法,具有定位精度高、计算时间少以及对噪声数据不敏感等特点.其突出的优点是在不增加计算难度和计算时间的前提下,用矩阵反演的方法确保得到全局最小值解,可适应于地震早期预警,海啸早期预警,以及大震速报等实际工作. 展开更多

关键词 射线追踪 全局直接节点搜寻 矩阵反演 全局解 到时残差空间

在线阅读 下载PDF 职称材料

多震相菲涅耳体射线走时同时反演成像

20

作者 李兴旺 白超英 黄国娇 《地震学报》 CSCD 北大核心 2014年第5期895-909,981,共15页

高频假设下的地震射线理论以及相应的地震成像理论表明,在射线稀疏条件下,不可能得到较高分辨率的构造成像;而有限频射线理论更符合实际地震的传播规律,即地震波的走时不仅与中心射线(传统的几何射线)上的速度分布有关,而且与中心射线... 高频假设下的地震射线理论以及相应的地震成像理论表明,在射线稀疏条件下,不可能得到较高分辨率的构造成像;而有限频射线理论更符合实际地震的传播规律,即地震波的走时不仅与中心射线(传统的几何射线)上的速度分布有关,而且与中心射线附近一定范围(称其为第一菲涅耳体)内的速度异常分布有关.鉴于此,本文提出了计算多震相地震波菲涅耳体有限频射线的方法,并定义了走时敏感核函数,同时给出了利用多震相菲涅耳体有限频射线进行速度模型和反射界面同时反演成像的公式.利用多震相走时资料,使用传统射线层析成像方法与有限频射线层析成像方法进行了速度和界面的同时反演成像.结果表明,当射线密度较小时,无论是对速度模型的重建还是对反射界面几何形状的更新,有限频射线层析成像方法均优于传统射线层析成像方法,而变频有限频射线层析成像则是实际地震层析成像的首选反演算法. 展开更多

关键词 第一菲涅耳体 有限频射线 走时敏感核函数 多震相走时 同时反演 菲涅耳体射线成像

在线阅读 下载PDF 职称材料