利用变步长网格声波波场模拟的时间域全波形反演方法
【技术领域】
[0001] 本发明涉及石油勘探技术，是一种应用于地震数据深度建模中的利用变步长网格 声波波场模拟的时间域地震波形反演方法。
【背景技术】
[0002] 地震数据处理的重要目的之一是通过深度偏移（D巧thMigration)得到高质量的 深度域图像。深度域图像的质量依赖于地震速度模型的质量，因此建模技术在整个地震数 据处理起着举足轻重的作用。全波形反演（化11Wave化rmInversion)是一种数据域（Data Domain)的地震建模技术，其基本原理是通过比较正演得到的模拟数据和实际观测数据的 差别来确定模型改变量，反复迭代此过程直到模拟数据和观测数据的差别达到最小。和传 统的基于图像域（ImageDomain)的建模方法幼日层析反演)相比，全波形反演方法的优点在 于（1)分辨率高；（2)需要的人工干预少。
[0003]目前业内全波形反演方法的实现方法主要有两个方向；（1)在时间域实现；（2) 在频率域实现。其中时间域算法中，高精度高效的波场模拟算法是其关键因素。在众多模 拟算法中，有限差分方法是应用最为普遍的一个。在有限差分算法中，受计算机内存和计 算量的限制，尽管空间导数会使用高阶逼近，但时间导数一般都使用二阶逼近，由于时间方 向引入的误差偏大，会导致所谓的数字频散现象。该种误差在正演模拟的结果中会表现为 走时误差和相位旋转。为提高时间方向的模拟精度，有学者提出了伪谱法（Tal-Ezeret al，1987，Etgen，1989，aiangandaiang2009)，但该些方法要在毎个时间步中多次运用髙 维付氏变换，比有限差分方法运算量大很多。
【发明内容】
[0004] 本发明的目的是提供一种解决时间域全波形反演实现中由于时间方向模拟精度 不足导致的数值频散现象，并大大降低的运算量，提高效率的利用变步长网格声波波场模 拟的时间域全波形反演方法。
[0005] 本发明实现步骤如下：
[0006] (1)在工区沿炮线方向人工激发地震波，由检波器记录地震波波场；
[0007] (2)根据工区地质信息等信息建立初始速度场；延深度方向用变步长网格进行剖 分，剖分方法如下：根据每一层波场最小传播速度W及模拟数据最大频率确定该层的网格 大小，网格大小不大于保证最短波长的一半。保证稳定性，还要保证与上一层网格大小相 比，增大比例不能太大。
[0008] (3)采用星形格格点分布，在星形權积算子的网格对角线上增加4个计算的格点；
[0009] (4)W速度为常速，在时间-波数域求出声波方程的解析解，得到一个精确的波场 随时间演化关系式：
[0010]戶川，/ +A〇二2cos(vAf |k I.)户(k,〇-/3(k么咕 Q)
[0011] 式中：
[001引 度:为波场，V为当地速度，At为时间步长，k为波数，t为时间；
[0013] (5)在复数域根据空间点的波场传播速度和时间延拓步长用最小二乘寻优的方法 计算星形网格格点的權积系数，使星形网格各点对应的波数域时间演化关系与步骤4)中精 确的波场随时间演化关系式中演化关系间的误差最小化；
[0014] (6)在模拟空间范围内，分析速度变化范围，并建立速度索引，利用步骤5)中最小 二乘寻优的方法，分别求取对应速度的權积算子；
[0015] (7)在时间空间域，利用得到的權积算子和速度索引进行波场正演模拟并记录模 拟的地震波场；
[0016] (8)利用高斯-牛顿优化方法使得野外记录的地震波场和正演模拟的地震波场误 差减小，迭代更新速度模型完成反演。
[0017] 本发明相对于现有技术具有W下优点和技术效果：
[001引（1)解决了现有算法时间导数近似计算的精度不足问题。利用声波方程在时间波 数域的表达形式，得到精确的波场时间演化关系。可W证明波数域中二阶时间差分的有限 差分格式是该个时间演化关系的二阶截断。因此时间演化关系提高了时间导数近似计算的 精度。
[0019] (2)解决了高维化urier变换运算量太大的问题。在空间域求取權积算子，利用 最小二乘法，使之在波数域的有效波数带内拟合时间演化算子，避免多次运用高维化urier 变换，减少运算量。
[0020] (3)解决了波场模拟运算量大计算效率低的问题。尽管在每一个时间步中，时间 演化算子的计算量比相应阶数的有限差分算子略大，但要达到相同的波场模拟精度，有限 差分法必须采用小于数倍的时间步长。因此时间演化算子在精度上和效率上都有很强的优 势。
[0021] (4)采用延深度变步长网格，进一步减少了运算量，解决了波场模拟运算量大计算 效率低的问题。
【附图说明】
[0022] 图1是变步长网格显格式时间演化算子的星形格点分布图；
[0023] 图2是时间演化算法与差分方法分别在近场和远场的单道模拟波形对比图；
[0024] 图3是时间演化算法与差分方法远道模拟效率对比图。
【具体实施方式】
[00巧]W下结合附图和实例详细描述本发明。
[0026] 图1是变步长网格显格式时间演化算子的星形格点分布图；
[0027] 本发明所要解决的技术问题是一种应用于全波场反演的高精度高效率声波波场 模拟技术。具体的实施方案如下：
[0028] (1)在工区沿炮线方向人工激发地震波，由检波器接收并记录地震波场，得到地震 观测数据du(t)=d(X。Xj,t)，其中Xi和Xj,分别代表震源和接收点的空间位置，t代表观测 时间。
[0029] (2)利用采样定理，根据每一层波场传播速度W及模拟数据主频的建立变步长空 间网格，对工区利用延深度的变步长网格进行剖分，根据工区地质信息等建立初始速度场。
[0030] (3)采用星形格格点分布，在星形權积算子的网格对角线上增加4个计算的格点；
[0031] (4)W速度为常速，在时间-波数域求出声波方程的解析解，得到一个精确的波场 随时间演化关系式：
[0032] 例k,/+A〇 二 2cos〇Af|k|)度(!(，。一^(k -A〇 (1)
[0033]式中：
[0034]多:为波场，V为当地速度，At为时间步长，k为波数，t为时间；
[0035] 基于初始速度模型V，在变步长空间网格上应用本发明的波场模拟方法进行数字 正演，模拟地震观测数据，得到地震模拟数据Uu(t)=u(X。Xj,t)。
[0036] (5)在复数域根据空间点的波场传播速度和时间延拓步长用最小二乘寻优的方法 计算星形网格格点的權积系数，使星形网格各点对应的波数域时间演化关系与步骤4)中精 确的波场随时间演化关系式中演化关系间的误差最小化；
[0037] 图2是时间演化算法与差分方法分别在近场和远场的单道模拟波形对比图；在近 场位置，时间演化算法与差分方法模拟波形与参考值接近，误差大致相等；随时间推移，远 场位置时间演化算法的波形与参考值吻合很好，而差分方法发生频散。
[0038] 图3是时间演化算法与差分方法远道模拟效率对比图；如要达到相同的传播精 度，差分算法的时间传播步长需要减少至时间演化算法的1/5。
[0039] (6)在模拟空间范围内，分析速度变化范围，并建立速度索引，利用步骤5)中最小 二乘寻优的方法，分别求取对应速度的權积算子；
[0040] (7)在时间空间域，利用得到的權积算子和速度索引进行波场正演模拟并记录模 拟的地震波场，对观测数据和模拟数据求取残差ry(t)=du(t)-Uu(t)，利用残差ry(t)反演 校正方向g(x。Xj,Xk)，反演过程中应用本发明的波场模拟方法进行炮点正向波场模拟和残 差反向传播的波场模拟；
[0041](8)利用求取的校正方向g(x。Xj,Xk)和合理的步长，求取速度模型修正量，对初始 模型进行修改，利用高斯-牛顿优化方法使得野外记录的地震波场和正演模拟的地震波场 误差ry(t)减小，设立目标函数，多次迭代反演过程，使目标函数最小化，改进速度模型的 精确度。
[0042] 综上所述，本发明提供了一种时间域高效率高精度的波场模拟方法，在时间域全 波形反演技术的实现过程中，可W用此方法代替原来的有限差分波场模拟方法，提高波场 模拟的精度，有利于波形的匹配。
【主权项】
1. 一种利用变步长网格声波波场模拟的时间域全波形反演方法，特点是采用步骤如 下： 1) 在工区沿炮线方向人工激发地震波，由检波器记录地震波波场； 2) 根据工区地质信息建立初始速度场，在延深度方向用变步长网格进行剖分； 3) 采用星形格格点分布，在星形褶积算子的网格对角线上增加4个计算的格点； 4) 以速度为常速，在时间-波数域求出声波方程的解，得到一个精确的波场随时间演 化关系式： 式中：f为波场，V为当地速度，At为时间步长，k为波数，t为时间； 5) 在复数域根据空间点的波场传播速度和时间延拓步长用最小二乘寻优的方法计算 星形网格格点的褶积系数，使星形网格各点对应的波数域时间演化关系与步骤4)中精确的 波场随时间演化关系式中演化关系间的误差最小化； 6) 在模拟空间范围内，分析速度变化范围，并建立速度索引，利用最小二乘寻优的方 法，分别求取对应速度的褶积算子； 7) 在时间空间域，利用得到的褶积算子和速度索引进行波场正演模拟并记录模拟的地 震波场； 8) 利用高斯-牛顿优化方法使得野外记录的地震波场和正演模拟的地震波场误差减 小，迭代更新速度模型完成反演。2. 根据权利要求1的方法，特点是步骤2)所述的剖分是根据每一层波场最小传播速度 以及模拟数据最大频率确定该层的网格大小，网格大小不大于保证最短波长的一半，并且 保证与上一层网格大小相比，增大比例不能过大。
【专利摘要】本发明是一种应用于地震数据深度建模中的利用变步长网格声波波场模拟的时间域地震波形反演方法。延深度方向用变步长网格进行剖分，在星形褶积算子的网格对角线上增加计算格点，在时间-波数域得到精确的波场随时间演化关系式，根据空间点的波场传播速度和时间延拓步长计算网格格点的褶积系数，使各点对应的波数域时间演化关系误差最小化，在模拟空间范围内，并建立速度索引，分别求取对应速度的褶积算子，利用褶积算子和速度索引波场正演模拟，使得野外记录的地震波场和正演模拟的地震波场误差减小，迭代更新速度模型完成反演。本发明解决了现有算法时间导数近似计算的精度不足问题，解决了波场模拟运算量大计算效率低的问题。
【IPC分类】G01V1/28
【公开号】CN104977607
【申请号】CN201410140040
【发明人】戴南浔, 武威
【申请人】中国石油集团东方地球物理勘探有限责任公司
【公开日】2015年10月14日
【申请日】2014年4月9日