国家自然科学基金(41274116)
- 作品数：8 被引量：105H指数：7
- 相关作者：刘玉柱董良国杨积忠王义谢春更多>>
- 相关机构：同济大学中国石油更多>>
- 发文基金：国家自然科学基金国家科技重大专项中央高校基本科研业务费专项资金更多>>
- 相关领域：天文地球更多>>
- 声波全波形反演目标函数性态被引量：42
- 2013年
- 地震波传播的复杂性所引起的地震反演中强烈的非线性问题是目前全波形反演在技术上遇到的最大难题,了解全波形反演中不同的目标函数随不同物性参数的不同摄动尺度的变化性态,对选择合理的反演方法和反演策略具有重要意义.本文参照Jannane等对波形反演目标函数性态的分析方法,通过变密度声波方程,分析了多种地震数据子集的不同目标函数随物性参数的摄动尺度的变化关系,重点分析了它们的非线性程度,为进行分步骤、分尺度全波形反演方法和反演策略的选择提供了理论指导.
- 董良国迟本鑫陶纪霞刘玉柱
- 关键词：目标函数非线性
- 基于Born波路径的高斯束初至波波形反演被引量：19
- 2014年
- 为了提高表层速度反演精度,本文提出了一种新的波形反演方法.该方法只利用初至波波形信息以减少波形反演对初始模型的依赖性,降低反演多解性与稳定性.由于只利用初至波波形信息,所以该方法利用高斯束计算格林函数和正演波场,以减少正演计算量.为了避免庞大核函数的存储,该方法基于Born波路径,利用矩阵分解算法实现方向与步长的累加计算.将此基于Born波路径的初至波波形反演方法应用于理论模型实验,并与声波方程全波形反演和初至波射线走时层析方法相对比,发现该方法的反演效果略低于全波形反演方法,但明显优于传统初至波射线走时层析方法,而计算效率却与射线走时层析相当.同时,相对于全波形反演,本文方法对初始模型的依赖性也有所降低.
- 刘玉柱谢春杨积忠
- 关键词：核函数矩阵分解
- 有效利用初至信息的偏移距加权地震层析成像方法被引量：11
- 2014年
- 在基于射线理论的初至波地震走时层析成像方法中,很多因素会影响反演效果,其中观测数据能否有效利用是影响反演结果精度的一个重要因素。数值实验表明,初至层析成像反演表层速度结构时,不同偏移距数据对于提高不同深度反演精度的贡献是不同的。根据这一认识,提出了多偏移距(MOR)层析反演策略,以提高表层建模反演的精度;通过理论模型与实际资料测试证实了该策略的有效性。但MOR反演策略实现过程复杂,且只能实现偏移距的离散利用,考虑到多偏移距层析实质上是对观测数据的一种加权实现,进一步提出了利用偏移距构建协方差矩阵的偏移距加权地震层析成像方法。模型数据与实际资料的反演实验结果表明,偏移距加权地震层析成像方法避免了MOR层析过程中的多次嵌套反演,只需一次反演就可以达到连续多偏移距层析策略的反演效果。
- 刘玉柱杨积忠
- 关键词：初至波层析成像协方差矩阵
- 变密度声波方程多参数全波形反演策略被引量：28
- 2014年
- 多参数全波形反演中各参数之间的相互耦合增加了反演的非线性程度.通过分析各参数之间的相互影响,提出合理的多参数反演策略是解决该问题的有效途径.本文从变密度声波方程出发,首先研究了密度在速度反演中的重要作用,然后分析了速度对密度反演的影响程度,进而提出了一种有利于速度、密度分步联合反演的策略.第一步,利用给定的初始模型对速度、密度进行同时反演,得到比较可靠的速度反演结果;第二步,利用第一步反演得到的速度和给定的初始密度作为初始模型,继续进行双参数同时反演,这样可以同时得到比较可靠的速度、密度反演结果.为了进一步提高反演精度,将第二步反演得到的速度、密度作为初始模型,再进行下一轮双参数联合反演.二维理论模型实验结果充分说明了本文提出的这种反演策略的有效性.
- 杨积忠刘玉柱董良国
- 关键词：非线性
- VTI介质多参数联合走时层析成像方法被引量：17
- 2014年
- 本文基于球谐展开群速度表达式计算走时关于各向异性参数的Fréchet核函数,利用共轭梯度法对两种参数化方法进行了VTI介质中多参数联合反演方法研究.经过理论分析和数值试验发现,与经典的Thomsen参数化方法相比,垂直慢度、水平慢度与动校正慢度的参数化方式更有利于VTI介质多参数联合走时层析反演.为了克服走时对ε参数的不敏感性,我们采用了两步法进行双参数反演,理论模型试验反演得到了与垂直速度精度相当的ε参数.可以将两步法扩展到三步法以同时反演各向异性介质中的三个参数,数值试验展示了该策略的应用潜力.
- 刘玉柱王光银董良国杨积忠
- 关键词：VTI介质联合反演初至波层析成像
- 一种改进的全波形反演混合迭代优化方法(英文)
- 2013年
- 全波形反演中目标函数的Hessian信息对加速收敛起着重要作用,但直接计算Hessian矩阵及其逆通常是不可行的。有限内存BFGS(limited memory Broyden-Fletcher-Goldfarb-Shannon,L-BFGS)法或Hessian-free不精确牛顿(Hessian-free inexact Newton,HFN)法可以使用近似的Hessian信息,但使用程度有限。这两种方法能够互相提供Hessian信息,因而可以混合迭代。混合方法的性能依赖于二者间的有效转换。本文设计了基于目标函数下降比率(下降百分比)的迭代方法动态转换的新方案,得到一种改进的混合迭代优化方法。通过比较相同计算代价下两种方法的目标函数下降比率的大小,新方案使混合方法总是执行下降最快的迭代方法。Marmousi和Overthrust模型的数值试验表明,在保证反演质量的同时,改进方法的收敛速度明显快于L-BFGS法,比Enriched法有小幅提升。它也改进了HFN法效率低的不足。
- 王义董良国刘玉柱
- L-BFGS法时间域全波形反演中初始矩阵的选择方法被引量：7
- 2014年
- 有限内存BFGS(Limited-memory Broyden-Fletcher-Goldfarb-Shanno,L-BFGS)方法是地震全波形反演(Full waveform inversion,FWI)中应用最广泛的优化方法之一。该方法需要提供Hessian逆矩阵的一个初始近似矩阵,初始近似矩阵可以在每次迭代中都进行更新,也可以固定不变。L-BFGS法的收敛性能与初始矩阵的选择以及更新方式紧密相关。在分析了全波形反演中Hessian逆矩阵的几种不同近似方式之后,分别将它们作为初始矩阵应用到了时间域L-BFGS法全波形反演中,详细比较了不同近似方式对反演精度和效率的影响,为LBFGS法全波形反演中初始矩阵的选择提供了依据。数值试验结果表明,在不同的初始矩阵选择以及不同更新方式下,L-BFGS法的收敛速度表现出明显差异。当初始矩阵固定不变时,以虚震源为基础的两种方式收敛最快,采用波场能量矩阵的两种方式要慢一些,其中包含有记录残差反传波场能量的方式收敛最慢。初始矩阵每次迭代更新时,各种方式的收敛速度都慢于固定不变时的情形。此外,同种方式分别作为L-BFGS法的初始矩阵和共轭梯度法的预条件矩阵时,L-BFGS法的性能要优于预条件共轭梯度法。
- 王义董良国
- 关键词：预条件
- 基于地震数据子集的波形反演思路、方法与应用被引量：15
- 2015年
- 地震数据与地下介质物性参数之间的复杂关系,决定了地震全波形反演在理论方法上面临着强烈的非线性难题.地下不同物性参数的不同分量在地震数据上具有不同的表现,勘探的不同阶段对地下介质模型的精度也具有不同的要求,这就决定了在地震全波形反演过程中不必时刻追求地震数据全部信息的匹配,部分信息的匹配就有可能解决现阶段的某些问题,还可以一定程度上规避匹配全部地震信息所遇到的强烈非线性难题.基于这样的考虑,我们提出了利用地震数据子集进行波形反演的思路,给出了统一的反演方法,并通过基于包络数据子集以及反射波数据子集的波形反演的理论模型与实际资料反演试验,证明了所提出的波形反演思路和方法的正确性.
- 董良国黄超迟本鑫刘玉柱
- 关键词：波形反演目标函数核函数
