基于分辨率可控包络生成算子的多尺度全波形反演方法技术

本发明专利技术公开了一种基于分辨率可控包络生成算子的多尺度全波形反演方法：首先，建立分辨率可控包络生成算子，通过调整参数控制反演使用地震数据的频带范围；其次，基于分辨率可控包络生成算子构建全波形反演方法，推导目标函数对模型参数的梯度以及迭代步长的计算公式；最后，使用最速下降法对模型参数进行迭代更新，直到满足算法结束条件为止。本发明专利技术具有更强的克服跳周问题的能力，更适用于岩性油气藏速度模型建模。特别是当地震数据中低频信息缺乏或主频过高时，优势更为明显。本发明专利技术无论地震数据中是否具有充足的低频信息，经过分辨率可控包络生成算子作用后的地震数据具有充足的低频信息；本发明专利技术的参数可调，可以实现多尺度反演。

Multi scale full waveform inversion method based on resolution controllable envelope generating operator

The invention discloses a multi-resolution controlled envelope generator operator full waveform inversion method based on: first, the establishment of resolution controllable envelope generation operator, by adjusting the control parameters inversion of seismic data using frequency range; secondly, the resolution of a building envelope controlled full waveform inversion operator based on the gradient of objective function derivation of the model parameters and the calculation formula of the iterative step; finally, using the steepest descent iterative updates the model parameters, until meeting the conditions so far end. The invention has a stronger ability to overcome the cycle problem, and is more suitable for modeling the velocity model of lithologic oil and gas reservoirs. Especially, when the low frequency information in the seismic data is lack or the main frequency is too high, the advantage is more obvious. The invention has enough low frequency information in the seismic data, and the seismic data after the action of the controllable envelope generation operator has sufficient low frequency information. The parameters of the invention can be adjusted, and multi-scale inversion can be realized.

【技术实现步骤摘要】
基于分辨率可控包络生成算子的多尺度全波形反演方法
本专利技术属于地球物理勘探
，涉及一种地球模型P波速度参数反演方法，特别涉及一种基于分辨率可控包络生成算子的多尺度全波形反演方法。
技术介绍
石油和天然气是事关国家经济发展及国家安全的战略资源，进一步推进我国油气资源的勘探与开发属于国家重大战略。岩性油气藏是我国今后相当长时间内最主要的勘探开发对象。地震勘探方法是油气勘探的主要方法，速度参数是油气检测的重要参数。我国岩性油气藏大多在地质上纵向呈薄互层结构、横向上具有强非均质性，如何对岩性油气藏的速度参数进行精细描述是我国油气勘探面临的关键问题。全波形反演是一类直接从叠前地震数据出发，利用叠前炮集地震数据全部信息来得到高精度地下介质速度模型的方法。因此，全波形反演是刻画岩性油气藏速度参数的有效方法。传统的全波形反演是一种以得到地下介质高分辨率参数模型为目标的数据匹配方法。全波形反演通过最小化一个衡量观测数据与计算数据匹配程度的目标函数来得到地下介质参数。最早将全波形反演方法推向实用的是Lailly和Tarantola的工作。他们提出可以使用由震源激发的正传波场以及由数据残差激发的反传波场进行互相关来构建目标函数对模型参数的梯度。这一误差反传技术克服了显式求取Frechet导数需要大量波场正演这一缺陷。传统全波形反演使用L2范数来衡量观测数据与计算数据的匹配程度，其应用于反射地震资料构建高分辨率速度模型时面临很多问题，例如：全波形反演要求观测数据中包含充足的低频信息或者初始模型足够准确来保证算法可以克服周波跳跃问题进而避免陷入局部极值。因此，降低全波形反演对初始模型的依赖程度就成为全波形反演领域的一个研究热点。Bunks在时间域提出了多尺度全波形反演方法，之后Pratt提出了频率域全波形反演方法。这两种方法在反演过程中首先使用低频数据进行反演。由于低频数据构成的目标函数相比于高频数据更加线性，全局极小值的收敛域更大，所以此时算法对初始模型的要求就会更低。基于低频数据的反演结果可以作为基于高频数据的全波形反演的初始模型。经过在反演中逐渐加入高频数据，最终可以得到高精度的速度模型。虽然多尺度方法有一些成功的应用案例，但是这类方法并不能从根本上克服跳周问题。当这类方法应用于实际资料时会遇到如下问题：通常建立一个足够精确的初始模型是很困难的；实际资料的低频信息往往被噪音影响；实际资料的偏移距往往是有限的。很多学者研究了为全波形反演建立比L2范数更优的度量准则这一问题。他们的目标是建立凸性更好，鲁棒性更强并且能够很好克服跳周问题的目标函数。与传统方法进行波形匹配不同，这一方法强调对波场旅行时信息的匹配。这一思路具有代表性的工作就是Luo和Schuster提出的波动方程旅行时反演。与传统旅行时反演方法依赖于人工拾取波场旅行时信息不同，该方法使用波动方程来得到正演数据，可以被视为一种使用互相关目标函数的全波形反演方法。然而，由于波动方程旅行时反演方法产生的速度模型过于光滑，其不能保证可以为全波形反演建立足够精确的初始模型。Chi等以及Luo等提出了在全波形反演中使用波场的瞬时振幅(包络)来构建反演的目标函数。Wu等证明了地震数据包络的振荡和衰减携带了超低频的信息，这一信息可以被用来估计长波长的速度模型。但是这一目标函数仍然会受到周波跳跃问题的困扰。
技术实现思路
本专利技术的目的在于克服上述现有技术的缺点，提供一种基于分辨率可控包络生成算子的多尺度全波形反演方法。本专利技术的目的是通过以下技术方案来实现的：这种基于分辨率可控包络生成算子的多尺度全波形反演方法，包括以下步骤：1)采集原始地震数据，然后对采集到的地震数据进行预处理，处理后得到叠前炮集地震数据，记为dobs(t,xr；xs)，其中t表示时间变量；xr,xs分别为检波器和震源的位置；2)构建长方形网格地质模型，设定离散网格的空间离散间隔、网格规模以及正演模型时间采样间隔；3)分析地震数据，给出反演使用的地震数据频带范围；4)给定地质模型P波速度v(x)的初始模型，并指定待优化的目标函数J(v(x))；5)求取目标函数J(v(x))对模型参数v(x)的梯度，使用误差反传技术进行高效求解；6)求取最速下降法迭代时的迭代步长，在频率域使用线性搜索方法确定；7)通过最速下降法迭代得到目标函数的全局极小值，进而得到在第3)步给定的频带范围内，P波速度的反演结果；8)调整反演时使用的地震数据的频带范围，并使用步骤7)的结果作为初始模型，重复步骤3)到7)，循环完所有的有效频带范围，输出最优P波速度模型。进一步，上述步骤2)中，离散网格的空间离散间隔、网格规模以及正演模型时间采样间隔是根据实际地震数据的炮检距、有效频带范围和采样时间决定的；选定网格参数的标准是使得基于该网格进行有限差分正演模拟时，不仅满足稳定性条件而且有效压制数值频散。进一步，上述步骤3)中，首先分析接收地震数据的有效频带范围，在反演初期先使用相对低频的地震数据，随着反演的进行，逐渐扩大反演使用的频带范围。进一步，上述步骤4)中，地质模型P波速度v(x)的初始模型是对真实P波速度模型的一个粗略估计，初始模型设定为横向均匀、纵向线性增加的速度模型；所述目标函数是刻画最优解的标准，以计算地震数据与观测地震数据之间的拟合程度或误差大小为标准。进一步的，设定衡量观测炮集地震数据与基于模型的计算炮集地震数据之间的分辨率可控包络数据残差的2范数为目标函数，即目标函数为：式中，T为正演的时间长度；G(.)为分辨率可控包络生成算子；Ss,r为一约束算子来保证只有震源和检波器分别位于xs和xr的波场数据被用来构建目标函数；ucal(t，x；xs)为震源位于xs的计算波场；目标函数值J(v(x))越小说明速度模型v(x)越准确。进一步，上述步骤5)中，目标函数J(v(x))对模型参数v(x)的梯度是使用误差反传技术进行高效求解。进一步，上述步骤6)中，迭代步长是在频率域使用线性搜索方法进行确定。进一步，上述步骤8)中，当使用某一频带范围内的地震数据进行迭代反演后，下一轮反演使用的地震数据频带范围将比上一轮迭代时宽，使更多高频的信息被用于进行反演。本专利技术具有以下有益效果：本专利技术所提出方法的步骤可以概括为以下三步：首先，本专利技术建立了一种分辨率可控包络生成算子，该算子可以通过调整参数控制反演使用地震数据的频带范围；其次，本专利技术基于分辨率可控包络生成算子，构建了全波形反演方法，推导了其目标函数对模型参数的梯度以及迭代步长的计算公式；最后，使用最速下降法对模型参数进行迭代更新，直到满足算法结束条件为止。与常用全波形反演方法相比，该方法具有更强的克服跳周问题的能力，更适用于岩性油气藏速度模型建模。特别是当地震数据中低频信息缺乏或主频过高时，该方法的优势更为明显。本专利技术所提出方法具有上述优势的原因在于：首先，无论地震数据中是否具有充足的低频信息，经过分辨率可控包络生成算子作用后的地震数据具有充足的低频信息；其次，分辨率可控包络生成算子的参数可调，可以实现多尺度反演。附图说明图1是本专利技术流程示意图；图2给定三组不同参数时移线性算子的振幅谱；图3是Marmousi模型算例的结果：(a)真实速度模型；(b)初始速度模型；(c)0-8Hz数据反演本文档来自技高网...

【技术保护点】
一种基于分辨率可控包络生成算子的多尺度全波形反演方法，其特征在于，包括以下步骤：1)采集原始地震数据，然后对采集到的地震数据进行预处理，处理后得到叠前炮集地震数据，记为dobs(t,xr；xs)，其中t表示时间变量；xr,xs分别为检波器和震源的位置；2)构建长方形网格地质模型，设定离散网格的空间离散间隔、网格规模以及正演模型时间采样间隔；3)分析地震数据，给出反演使用的地震数据频带范围；4)给定地质模型P波速度v(x)的初始模型，并指定待优化的目标函数J(v(x))；5)求取目标函数J(v(x))对模型参数v(x)的梯度，使用误差反传技术进行高效求解；6)求取最速下降法迭代时的迭代步长，在频率域使用线性搜索方法确定；7)通过最速下降法迭代得到目标函数的全局极小值，进而得到在第3)步给定的频带范围内，P波速度的反演结果；8)调整反演时使用的地震数据的频带范围，并使用步骤7)的结果作为初始模型，重复步骤3)到7)，循环完所有的有效频带范围，输出最优P波速度模型。

【技术特征摘要】
1.一种基于分辨率可控包络生成算子的多尺度全波形反演方法，其特征在于，包括以下步骤：1)采集原始地震数据，然后对采集到的地震数据进行预处理，处理后得到叠前炮集地震数据，记为dobs(t,xr；xs)，其中t表示时间变量；xr,xs分别为检波器和震源的位置；2)构建长方形网格地质模型，设定离散网格的空间离散间隔、网格规模以及正演模型时间采样间隔；3)分析地震数据，给出反演使用的地震数据频带范围；4)给定地质模型P波速度v(x)的初始模型，并指定待优化的目标函数J(v(x))；5)求取目标函数J(v(x))对模型参数v(x)的梯度，使用误差反传技术进行高效求解；6)求取最速下降法迭代时的迭代步长，在频率域使用线性搜索方法确定；7)通过最速下降法迭代得到目标函数的全局极小值，进而得到在第3)步给定的频带范围内，P波速度的反演结果；8)调整反演时使用的地震数据的频带范围，并使用步骤7)的结果作为初始模型，重复步骤3)到7)，循环完所有的有效频带范围，输出最优P波速度模型。2.根据权利要求2所述的基于分辨率可控包络生成算子的多尺度全波形反演方法，其特征在于，所述步骤2)中，离散网格的空间离散间隔、网格规模以及正演模型时间采样间隔是根据实际地震数据的炮检距、有效频带范围和采样时间决定的；选定网格参数的标准是使得基于该网格进行有限差分正演模拟时，不仅满足稳定性条件而且有效压制数值频散。3.根据权利要求2所述的基于分辨率可控包络生成算子的多尺度全波形反演方法，其特征在于，所述步骤3)中，首先分析接收地震数据的有效频带范围，在反演初期先使用相对低频的地震数据，随...

【专利技术属性】
技术研发人员：高照奇，潘志斌，高静怀，张金淼，
申请(专利权)人：西安交通大学，中海油研究总院，
类型：发明
国别省市：陕西,61