克服全波形反演周波跳跃的方法及系统技术方案

本发明专利技术公开了一种克服全波形反演周波跳跃的方法及系统，该方法包括：基于初始速度模型和地震数据，获得深度偏移剖面；计算深度偏移剖面的结构特征张量；基于地震数据的包络信息获得包络反演的速度梯度；基于结构特征张量对速度梯度沿构造平滑预处理，获得预处理后的速度梯度；基于预处理后的速度梯度，更新速度模型。其优点在于：包络反演具有重构低波数背景速度场能力，联合采用地质构造约束的梯度预条件算子，可有效解决传统全波形反演周波跳跃，易于陷入局部极值的问题，进而使实际资料反演更加稳健，更易收敛到全局极小值。

Method and System for Overcoming Cycle Jump in Full Waveform Inversion

The invention discloses a method and system for overcoming cycle jump in full waveform inversion, which includes: obtaining depth migration profile based on initial velocity model and seismic data; calculating structural characteristic tensor of depth migration profile; obtaining velocity gradient of envelope inversion based on envelope information of seismic data; smoothing pre-processing of velocity gradient along structure based on structural characteristic tensor; The velocity gradient after pretreatment is obtained, and the velocity model is updated based on the velocity gradient after pretreatment. Its advantage is that envelope inversion has the ability to reconstruct low-wavenumber background velocity field. Combining with gradient preconditioning operator constrained by geological structure, it can effectively solve the problem of traditional full-waveform inversion cycle jump and easily fall into local extremum, thus making the inversion of actual data more robust and convergent to global minimum.

【技术实现步骤摘要】
克服全波形反演周波跳跃的方法及系统
本专利技术涉及石油物探领域，更具体地，涉及一种克服全波形反演周波跳跃的方法及系统。
技术介绍
随着油气勘探的日益精细化、复杂化，高精度的速度建模在地球物理处理中扮演着越来越重要的角色。由于地震波速度不仅决定偏移成像的质量而且与岩石性质关系密切，能够反映岩石类别和富含流体情况，最终还可以影响地震解释结果的可信性，因而地震波速度是一个非常重要的地层物性参数。地震全波形反演方法在反演过程中会匹配地震波的走时、振幅和相位信息，目前是建模精度最高的一种方法，由于其具有高度非线性性，对初始模型具有高度依赖性，当初始模型较差时，在匹配过程中会出现周波跳跃问题，从而使全波形反演陷入局部极小值的问题。目前业内为了克服全波形反演对初始速度模型的依赖，降低传统全波形反演的周波跳跃问题，采用Laplace域、Laplace-Fourier波形反演方法，利用阻尼波场Fourier变换的零频分量，反演得到低波数度模型；吴如山教授项目组、中国石油大学(华东)、同济大学借鉴天然地震中的包络反演思想，利用信号的包络开展低波数背景速度重建；法国傅里叶大学、吉林大学等知名高校采用旅行时主导匹配反演方法，在波形反演过程中更注重走时的匹配，在此条件下构建更加稳健的低波数背景速度方案。以上方法并不能使全波形反演收敛到全局极小值。因此，有必要开发一种有效解决波形反演的周波跳跃问题，使全波形反演更易收敛到全局极小值的全波形反演周波跳跃的方法及系统。公开于本专利技术
技术介绍
部分的信息仅仅旨在加深对本专利技术的一般
技术介绍
的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
技术实现思路
本专利技术提出了一种克服全波形反演周波跳跃的方法及系统，其能够通过采用地震数据的包络对地下介质模型的低波数背景速度模型成分进行重构，在反演迭代过程中采用地质构造导向约束预条件方法对梯度预处理，解决了波形反演的周波跳跃问题，使全波形反演更易收敛到全局极小值。根据本专利技术的一方面，提出了一种克服全波形反演周波跳跃的方法，所述方法包括：基于初始速度模型和地震数据，获得深度偏移剖面；计算所述深度偏移剖面的结构特征张量；基于所述地震数据的包络信息获得包络反演的速度梯度；基于所述结构特征张量对所述速度梯度沿构造平滑预处理，获得预处理后的速度梯度；基于所述预处理后的速度梯度，更新速度模型。优选地，获得所述深度偏移剖面包括：基于CMP地震数据道集构建所述初始速度模型；基于所述初始速度模型、所述地震数据和叠前深度偏移成像算子，获得所述深度偏移剖面。优选地，所述叠前深度偏移成像算子为Kirchhoff叠前深度偏移算子、单程波深度偏移算子或逆时偏移算子。优选地，通过图像学理论，计算所述深度偏移剖面的结构特征张量。优选地，一阶对称半正定的结构特征张量矩阵为：对上式进行特征值分解可得：T＝λuuuT+λvvvT式中，λu≥λv≥0为特征值；u，v分别为特征值对应的特征向量。优选地，所述速度梯度为：其中，E＝esyn(t)-eobs(t)；式中，y(t)为模拟的地震信号；yH(t)为模拟地震信号的希尔伯特变换；esyn(t)为模拟地震信号的包络eobs(t)为观测地震信号的包络；χ为L2范数约束下的地震信号包络和模拟地震信号包络的残差。优选地，地震信号的包络为：式中，u(t)为地震信号；uH(t)为地震信号的希尔伯特变换。优选地，地震信号包络和模拟地震信号包络的残差L2范数目标泛函数为：其中，优选地，所述更新速度模型为：vn+1＝vn+αGn式中，vn+1为第(n+1)次迭代的低波数背景速度；Gn为基于结构特征张量平滑预处理后的梯度；α为步长。根据本专利技术的另一方面，提出了一种克服全波形反演周波跳跃的系统，所述系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：基于初始速度模型和地震数据，获得深度偏移剖面；计算所述深度偏移剖面的结构特征张量；基于所述地震数据的包络信息获得包络反演的速度梯度；基于所述结构特征张量对所述速度梯度沿构造平滑预处理，获得预处理后的速度梯度；基于所述预处理后的速度梯度，更新速度模型。本法提出了一种克服全波形反演周波跳跃的方法及系统，其优点在于：包络反演具有重构低波数背景速度场能力，联合采用地质构造约束的梯度预条件算子，可有效解决传统全波形反演周波跳跃，易于陷入局部极值的问题，进而使实际资料反演更加稳健，更易收敛到全局极小值。本专利技术的方法和系统具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施例中将是显而易见的，或者将在并入本文中的附图和随后的具体实施例中进行详细陈述，这些附图和具体实施例共同用于解释本专利技术的特定原理。附图说明通过结合附图对本专利技术示例性实施例进行更详细的描述，本专利技术的上述以及其它目的、特征和优势将变得更加明显。图1示出了根据本专利技术的一个示例性实施例的一种克服全波形反演周波跳跃的方法的步骤的流程图。具体实施方式下面将参照附图更详细地描述本专利技术。虽然附图中显示了本专利技术的优选实施例，然而应该理解，可以以各种形式实现本专利技术而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本专利技术更加透彻和完整，并且能够将本专利技术的范围完整地传达给本领域的技术人员。本专利技术提供了一种克服全波形反演周波跳跃的方法，该方法包括：基于初始速度模型和地震数据，获得深度偏移剖面；计算深度偏移剖面的结构特征张量；基于地震数据的包络信息获得包络反演的速度梯度；基于结构特征张量对速度梯度沿构造平滑预处理，获得预处理后的速度梯度；基于预处理后的速度梯度，更新速度模型。作为优选方案，获得深度偏移剖面包括：基于CMP地震数据道集构建初始速度模型；基于初始速度模型、地震数据和叠前深度偏移成像算子，获得深度偏移剖面。其中，获得深度偏移剖面，得到地下介质构造f(x)，其中x表示空间向量。作为优选方案，叠前深度偏移成像算子为Kirchhoff叠前深度偏移算子、是单程波深度偏移算子或逆时偏移算子。其中，通过图像学理论，计算深度偏移剖面的结构特征张量。在图像学领域中可以把深度偏移剖面看做一幅图像，偏移图像的地质构造具有显著的方向性文理结构，求解文理结构得到结构张量以用于表征偏移图像的局部结构信息，即构造特征的主要变化方向及变化速率。作为优选方案，对于二维偏移剖面数据体，它在每点的一阶对称半正定结构张量矩阵为：对上式进行特征值分解可得：T＝λuuuT+λvvvT式中，λu≥λv≥0为特征值；u，v分别为特征值对应的特征向量。其中，特征值及特征向量刻画了局部图像的特征。当两个特征值均为零时，图像区域为常数；当两个特征值相等并都大于零时表示图像为各向同性；当两个特征值不相等并均大于零时表示图像为各向异性，且局部图像区域内存在主方向，较大特征值对应的特征向量表示图像梯度变化最快的方向，其与图像局部的线性特征方向垂直，相对而言较小特征值对应的特征向量代表了图像局部的线性方向。作为优选方案，速度梯度为：其中，E＝esyn(t)-eobs(t)；式中，y(t)为模拟的地震信号；yH(t)为模拟地震信号的希尔伯特变换；esyn(t)为模拟地震信号的包络；eobs(t)为观测地震信号的包络；χ本文档来自技高网...

【技术保护点】
1.一种克服全波形反演周波跳跃的方法，其特征在于，所述方法包括：基于初始速度模型和地震数据，获得深度偏移剖面；计算所述深度偏移剖面的结构特征张量；基于所述地震数据的包络信息获得包络反演的速度梯度；基于所述结构特征张量对所述速度梯度沿构造平滑预处理，获得预处理后的速度梯度；基于所述预处理后的速度梯度，更新速度模型。

【技术特征摘要】
1.一种克服全波形反演周波跳跃的方法，其特征在于，所述方法包括：基于初始速度模型和地震数据，获得深度偏移剖面；计算所述深度偏移剖面的结构特征张量；基于所述地震数据的包络信息获得包络反演的速度梯度；基于所述结构特征张量对所述速度梯度沿构造平滑预处理，获得预处理后的速度梯度；基于所述预处理后的速度梯度，更新速度模型。2.根据权利要求1所述的克服全波形反演周波跳跃的方法，其中，获得所述深度偏移剖面包括：基于CMP地震数据道集构建所述初始速度模型；基于所述初始速度模型、所述地震数据和叠前深度偏移成像算子，获得所述深度偏移剖面。3.根据权利要求2所述的克服全波形反演周波跳跃的方法，其中，所述叠前深度偏移成像算子为Kirchhoff叠前深度偏移算子、单程波深度偏移算子或逆时偏移算子。4.根据权利要求1所述的克服全波形反演周波跳跃的方法，其中，通过图像学理论，计算所述深度偏移剖面的结构特征张量。5.根据权利要求4所述的克服全波形反演周波跳跃的方法，其中，一阶对称半正定的结构特征张量矩阵为：对上式进行特征值分解可得：T＝λuuuT+λvvvT式中，λu≥λv≥0为特征值；u，v分别为特征值对应的特征向量。6.根据权利要求1所述的克服全波形反演周波跳跃的方法，其中，所述速度梯度为：其中，E＝esyn(t)...

【专利技术属性】
技术研发人员：胡光辉，刘定进，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京,11