利用固定网格声波波场模拟的时间域全波形反演方法
【技术领域】
[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)对工区利用规则的固定网格进行剖分，根据工区地质信息等建立初始速度场；
[0008] 3)W速度为常速，在时间-波数域求出声波方程的解析解，得到精确的波场随时 间演化关系式：
[0009] 户(k，/ +A/) = 2cos(vA/|l< |)/)(k，〇p{kJ-At). (1)
[0010] 式中：
[0011] 多为波场，V为当地速度，At为时间步长，k为波数，t为时间；
[0012] 4)根据空间点的波场传播速度和时间延拓步长，用最小二乘寻优的方法确定星形 网格各点的權积系数，使星形网格各点对应的波数域时间演化关系与步骤3)中精确的波场 随时间演化关系式中演化关系间的误差最小化；
[0013] 5)在模拟空间范围内，分析速度变化范围，并建立速度索引，利用步骤4)中最小二 乘寻优的方法，分别求取对应速度的權积算子；
[0014] 6)在时间空间域，利用得到的權积算子和速度索引进行波场正演模拟并记录模拟 的地震波场；
[0015] 7)利用高斯-牛顿优化方法使得野外记录的地震波场和正演模拟的地震波场误 差减小，迭代更新速度模型完成反演。
[0016] 本发明相对于现有技术具有W下优点和技术效果：
[0017] 解决了现有方法时间导数近似计算精度不足问题。利用声波方程在时间波数域的 表达形式，得到精确的波场时间演化关系。可W证明波数域中二阶时间差分的有限差分格 式是该个时间演化关系的二阶截断。因此时间演化关系提高了时间导数近似计算的精度。 解决了高维化urier变换运算量太大的问题。在空间域求取權积算子，利用最小二乘法，使 之在波数域的有效波数带内拟合时间演化算子，避免多次运用高维化urier变换，减少运 算量。解决了波场模拟运算量大计算效率低的问题。尽管在每一个时间步中，时间演化算 子的计算量比相应阶数的有限差分算子略大，但要达到相同的波场模拟精度，有限差分法 必须采用小于数倍的时间步长。因此时间演化算子在精度上和效率上都有很强的优势。
【附图说明】
[001引图1是规则网格显格式时间演化算子的星形格点分布图；
[0019] 图2是时间演化算法与差分方法分别在近场和远场的单道模拟波形对比图；
[0020] 图3是时间演化算法与差分方法远道模拟效率对比图。
【具体实施方式】
[0021] W下结合附图和实例详细描述本发明。
[0022] 本发明所要解决的技术问题是一种应用于全波场反演的高精度高效率声波波场 模拟技术。具体的实施方案如下：
[0023] (1)在工区沿炮线方向人工激发地震波，由检波器接收并记录地震波场，得到地震 观测数据错误！未找到引用源。，其中Xi和Xj.，分别代表震源和接收点的空间位置，t代表 观测时间。
[0024] (2)对工区利用规则的固定网格进行剖分，根据工区地质信息等建立初始速度场 V。图1是规则网格显格式时间演化算子的星形格点分布图。
[002引（3)W速度为常速，在时间-波数域求出声波方程的解析解，得到精确的波场随时 间演化关系式：
[0026] 公川，/ +A/) = 2cos〇A/k)公(k，/)_公(k…A/). (1)
[0027] 式中：
[0028] 多为波场，V为当地速度，At为时间步长，k为波数，t为时间；
[0029] 基于初始速度场V，应用本发明的精确的波场随时间演化关系式进行数字正演，模 拟地震观测数据，得到地震模拟数据错误！未找到引用源。；
[0030] (4)根据空间点的波场传播速度和时间延拓步长，用最小二乘寻优的方法确定星 形网格各点的權积系数，使星形网格各点对应的波数域时间演化关系与步骤（3)中精确的 波场随时间演化关系式中演化关系间的误差最小化；
[0031] 图2是时间演化算法与差分方法分别在近场和远场的单道模拟波形对比图；在近 场位置，时间演化算法与差分方法模拟波形与参考值接近，误差大致相等；随时间推移，远 场位置时间演化算法的波形与参考值吻合很好，而差分方法发生频散。
[0032] 图3是时间演化算法与差分方法远道模拟效率对比图；如要达到相同的传播精 度，差分算法的时间传播步长需要减少至时间演化算法的1/5。
[003引（5 )在模拟空间范围内，分析速度变化范围，并建立速度索引，利用步骤（4)中最小 二乘寻优的方法，分别求取对应速度的權积算子；
[0034] (6)在时间空间域，利用得到的權积算子和速度索引进行波场正演模拟并记录模 拟的地震波场，对观测数据和模拟数据求取残差ry(t)=dy(t)-Uu(t)，利用残差错误！未 找到引用源。反演校正方向g(x。Xj,Xk)，反演过程中应用本发明的波场模拟方法进行炮点 正向波场模拟和残差反向传播的波场模拟；
[00巧](7)利用求取的校正方向错误！未找到引用源。和合理的步长，求取速度模型修正 量，对初始模型进行修改，利用高斯-牛顿优化方法使得野外记录的地震波场和正演模拟 的地震波场误差ry(t)减小，设立目标函数，多次迭代反演过程，使目标函数最小化，改进 速度模型的精确度。
[0036] 综上所述，本发明提供了一种时间域高效率高精度的波场模拟方法，在时间域全 波形反演技术的实现过程中，可W用此方法代替原来的有限差分波场模拟方法，提高波场 模拟的精度，有利于波形的匹配。
【主权项】
1. 一利用固定网格声波波场模拟的时间域全波形反演方法，特点是采用以下步骤实 现： 1) 在工区沿炮线方向人工激发地震波，由检波器记录地震波波场； 2) 对工区利用规则的固定网格进行剖分，根据工区地质信息等建立初始速度场； 3) 以速度为常速，在时间-波数域求出声波方程的解析解，得到精确的波场随时间演 化关系式： 式中：?为波场，V为当地速度，At为时间步长，k为波数，t为时间； 4) 根据空间点的波场传播速度和时间延拓步长，用最小二乘寻优的方法确定星形网格 各点的褶积系数，使星形网格各点对应的波数域时间演化关系与步骤3)中精确的波场随时 间演化关系式中演化关系间的误差最小化； 5) 在模拟空间范围内，分析速度变化范围，并建立速度索引，利用步骤4)中最小二乘寻 优的方法，分别求取对应速度的褶积算子； 6) 在时间空间域，利用得到的褶积算子和速度索引进行波场正演模拟并记录模拟的地 震波场； 7) 利用高斯-牛顿优化方法使得野外记录的地震波场和正演模拟的地震波场误差减 小，迭代更新速度模型完成反演。
【专利摘要】本发明是用固定网格声波波场模拟的时间域全波形反演方法，利用规则的固定网格进行剖分建立初始速度场，在时间-波数域得到精确的波场随时间演化关系式，根据空间点的波场传播速度和时间延拓步长，用最小二乘寻优的方法确定星形网格各点的褶积系数，使星形网格各点对应的波数域时间演化关系与精确的波场随时间演化关系式中演化关系间的误差最小化，用最小二乘寻优的方法，分别求取对应速度的褶积算子，利用褶积算子和速度索引进行波场正演模拟并记录模拟的地震波场，使得野外记录的地震波场和正演模拟的地震波场误差减小，迭代更新速度模型完成反演。本发明解决了现有方法时间导数近似计算精度不足问题，运算量小，效率高。
【IPC分类】G01V1/28
【公开号】CN104977608
【申请号】CN201410140367
【发明人】戴南浔, 武威
【申请人】中国石油集团东方地球物理勘探有限责任公司
【公开日】2015年10月14日
【申请日】2014年4月9日