一种基于积分地震子波的时间域全波形反演方法技术

本发明专利技术属于海洋地震勘探领域，公开了一种基于积分地震子波的时间域全波形反演方法，包括如下步骤：输入野外观测地震数据和初始速度参数模型；根据野外观测地震数据提取地震子波，对地震子波进行时间积分，构建积分地震子波；以积分地震子波作为震源，数值求解时间域波动方程得到模拟地震数据；将野外观测地震数据与模拟地震数据逐点做差，以残差的时间积分作为震源，再次数值求解时间域波动方程，得到残差波场，同时计算梯度算子；计算速度参数模型更新量迭代反演，终止迭代并输出反演结果。本发明专利技术构建积分地震子波生成模拟地震记录，解决了时间域全波形反演对低频数据的压制问题，在梯度中增加了中低波数分量比重，提高了参数反演的精度。反演的精度。反演的精度。

【技术实现步骤摘要】
一种基于积分地震子波的时间域全波形反演方法


[0001]本专利技术属于海洋地震勘探领域，特别涉及一种基于积分地震子波的时间域全波形反演方法。

技术介绍

[0002]海洋地震反演通过观测到的地震数据估算相应的地球物理参数，进而反推海底地下的结构形态及物质成分，可以有效地识别地质构造、预测自然灾害和勘探油气藏。地震波的传播速度参数不仅是处理、解释地震数据的主要依据，而且其本身就是反映底下介质构造和岩性的重要数据。因此，求取地震波在地下传播的精确速度显得尤为重要。
[0003]全波形反演以波动理论为基础，匹配模拟地震数据与野外观测地震数据的波形信息建立目标函数，利用局部最优化的方法迭代寻找目标函数的全局最小值，可以高分辨率地恢复海底底质的速度参数，是目前海洋地震勘探领域的前沿方向之一。全波形反演可以在时间域实现，也可以在频率域实现，频率域全波形反演对计算机的要求更为苛刻，所以适用性更好的当属时间域全波形反演。但是，时间域全波形反演存在对低频数据的压制问题，不利于中低波数参数分量的恢复，容易使反演陷入局部极小值而导致失败。因此，解决时间域全波形反演这种低频压制问题不仅可以增加反演的容错率，同时可以提高反演的精度。

技术实现思路

[0004]为解决上述技术问题，本专利技术提供了一种基于积分地震子波的时间域全波形反演方法，以达到提高反演精度的目的。
[0005]为达到上述目的，本专利技术的技术方案如下：一种基于积分地震子波的时间域全波形反演方法，包括如下步骤：步骤一，输入数据：输入野外观测地震数据d
obs
(x,t)和初始速度参数模型m0；步骤二，构建积分地震子波：根据野外观测地震数据d
obs
(x,t)提取地震子波，对地震子波进行时间积分，构建积分地震子波：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)式中，t表示时间，f(t)表示根据野外观测地震数据提取的地震子波，表示时间积分后的地震子波，即积分地震子波；步骤三，正演模拟：以积分地震子波作为震源，数值求解时间域波动方程得到所有时间点的正演模拟波场u(x,t)，从正演模拟波场u(x,t)中得到接收点处的数据，即模拟地震数据d
cal
(x,t)，时间域波动方程表示形式如下：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)式中，u(x,t)表示正演模拟波场，x表示空间位置矢量，m表示速度参数模型，初始值为m0；表示空间导数之和；步骤四，梯度算子计算：将野外观测地震数据d
obs
(x,t)与模拟地震数据d
cal
(x,t)逐点做差，则残差表示为d
obs
(x,t)
‑
d
cal
(x,t)，以残差的时间积分作为震源替换方程(2)中的积分地震子波，然后再次数值求解时间域波动方程(2)，得到残差波场q(x,t)，同时根据方程(3)计算梯度算子；
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)式中，q(x,t)表示残差波场，s和g分别表示产生野外观测地震数据的震源点和接收点；步骤五，计算速度参数模型更新量迭代反演：先给予一个试探步长，再用Armijo条件的一维线搜索方法求取合适的步长作用于梯度算子计算本次迭代的速度参数模型更新量，保证每次迭代的速度参数模型更新量在20m/s
‑
100m/s之间，利用方程(4)计算更新后的速度参数模型m：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(4)式中，m0表示初始速度参数模型，λ表示步长，表示本次迭代的速度参数模型更新量；返回步骤三，将更新后的速度参数模型m代入方程(2)进行正演模拟，重新计算模拟地震数据，以野外观测地震数据与模拟地震数据残差的平方和[d
obs
(x,t)
‑
d
cal
(x,t)]2作为判定条件，判断反演是否收敛，进而选择继续迭代更新或是终止迭代并输出反演结果。
[0006]上述方案中，步骤四中，当判定条件小于门限值或连续5次增加，则认为反演已收敛。
[0007]优选地，所述门限值为0.001。
[0008]通过上述技术方案，本专利技术提供的一种基于积分地震子波的时间域全波形反演方法具有如下有益效果：常规的时间域全波形反演所用的子波是直接从野外观测地震数据中提取得到的，在计算梯度算子时，地震数据相当于通过了一个高频滤波器，低频信息被压制，进而速度参数的中低波数分量得不到有效的更新。本专利技术构建积分地震子波生成模拟地震记录，解决了时间域全波形反演对低频数据的压制问题，在梯度中增加了中低波数分量比重，提高了参数反演的精度。
附图说明
[0009]为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍。
[0010]图1为本专利技术实施例所公开的一种基于积分地震子波的时间域全波形反演方法流程示意图。
[0011]图2为Marmousi2速度参数模型，其中(a)为真实速度参数模型，(b)为初始速度参数模型。
[0012]图3 为地震子波与积分地震子波频带对比，其中，黑色实线为期待的频带，左边灰色实线为积分地震子波的频带，右边灰色实线为常规地震子波的频带。
[0013]图4 为反演速度参数模型对比，其中(a)为常规地震子波时间域全波形反演结果，(b)为积分地震子波时间域全波形反演结果。
[0014]图5 为单道反演速度曲线对比，其中(a)为常规地震子波时间域全波形反演结果，(b)为积分地震子波时间域全波形反演结果，其中黑色实线为真实速度曲线，灰色实线分别为反演速度曲线和初始速度曲线。
[0015]图6为梯度对比，其中(a)为常规地震子波时间域全波形反演梯度，(b)为积分地震子波时间域全波形反演梯度。
[0016]图7为不同距离处的单道梯度波数谱，其中(a)为距离0.8km处单道梯度波数谱，(b)为距离1.2km处单道梯度波数谱，(c)为距离1.4km处单道梯度波数谱，(d)为距离1.6km处单道梯度波数谱，其中黑色实线为常规子波时间域全波形反演梯度的波数谱，灰色实线为积分地震子波时间域全波形反演梯度的波数谱。
具体实施方式
[0017]下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述。
[0018]本专利技术提供了一种基于积分地震子波的时间域全波形反演方法，如图1所示，包括如下步骤：步骤一，输入数据：输入野外观测地震数据d
obs
(x,t)和初始速度参数模型m0；初始速度参数模型m0是已知的，其速度值从浅层到深层渐变增加，本实施例中选用Marmousi2速度参数模型，如图2中(b)所示。真实速度参数模型是未知的，如图2中(a)所示，需要通过本专利技术的方法反演得到。
[0019]步骤二，构本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于积分地震子波的时间域全波形反演方法，其特征在于，包括如下步骤：步骤一，输入数据：输入野外观测地震数据d
obs
(x,t)和初始速度参数模型m0；步骤二，构建积分地震子波：根据野外观测地震数据d
obs
(x,t)提取地震子波，对地震子波进行时间积分，构建积分地震子波：
ꢀꢀꢀꢀ
(1)式中，t表示时间，f(t)表示根据野外观测地震数据提取的地震子波， 表示时间积分后的地震子波，即积分地震子波；步骤三，正演模拟：以积分地震子波作为震源，数值求解时间域波动方程得到所有时间点的正演模拟波场u(x,t)，从正演模拟波场 u(x,t)中得到接收点处的数据，即模拟地震数据d
cal
(x,t)，时间域波动方程表示形式如下：
ꢀꢀꢀꢀꢀꢀꢀꢀ
(2)式中， u(x,t)表示正演模拟波场，x表示空间位置矢量，m表示速度参数模型，初始值为m0；表示空间导数之和；步骤四，梯度算子计算：将野外观测地震数据d
obs
(x,t)与模拟地震数据d
cal
(x,t)逐点做差，则残差表示为d
obs
(x,t)
‑
d
cal
(x,t)，以残差的时间积分作为震源替换方程(2)中的积...

【专利技术属性】
技术研发人员：梁展源，何传林，李冬，程广森，王振，白一奇，张晓语，于砚廷，龚金龙，尹佳玮，杜金燕，许岩，魏波，邢思宇，
申请(专利权)人：青岛欧谱赛斯海洋科技有限公司，
类型：发明
国别省市：