本申请公开了一种反演速度模型建模方法、装置、电子设备及介质，其中，该方法包括：步骤1：将多分量地震记录进行P/S分离，获得单分量的P波数据和S波数据，并进一步获得标量的PP反射波数据和标量的PS反射波数据，同时构建反演初始的纵波速度模型和横波速度模型；步骤2：基于PP反射波数据利用反射波波动方程走时反演，更新当前的纵波速度模型，完成纵波速度模型低波数建模；步骤3：基于PS反射波数据利用反射波波动方程走时反演，更新当前的横波速度模型，完成横波速度模型低波数建模；步骤4：将步骤2和步骤3得到的模型作为背景模型，利用弹性反射波波形反演，完成纵波速度模型中波数建模和横波速度模型中波数建模。横波速度模型中波数建模。横波速度模型中波数建模。
【技术实现步骤摘要】
反演速度模型建模方法、装置、电子设备及介质
[0001]本专利技术涉及石油地震勘探
，更具体地，涉及一种反演速度模型建模方法、装置、电子设备及介质。
技术介绍
[0002]地下介质的弹性效应在流体刻画、岩性识别、气云区成像、储层预测、裂缝识别及应力场刻画等方面有着重要作用。同时弹性参数(纵波速度、横波速度、密度等)也是油气勘探、开发以及动态监测等方面重要的地球物理参数。为更逼真地描述地震波的传播过程，有必要将反演方法从声介质RWI拓展到弹性介质中，但随着介质与波动方程的复杂程度增加，计算量也会不断的提升，所对应的问题也就越困难。
[0003]弹性波全波形反演(E
‑
FWI)是目前处理多分量地震数据较为热门的多参数反演方法。它通过最小化预测数据与观测数据的误差泛函，获取地下介质的弹性参数。尽管较声学近似下的FWI增加了一些计算量，E
‑
FWI也渐渐地在少量实际数据中得到应用(Sears et al.,2008；Brossier et al.,2009；Prieux et al.,2013a；Vigh et al.,2014)。随着对深层勘探需求的不断提高，利用反射波信息进行模型中深部更新越来越成为地震反演方法研究的重点。近几年，反射波反演方法也逐渐由单分量数据向多分量数据发展,弹性反射波走时与波形反演方法的研究逐渐开展，例如，Guo and Alkhalifah(2016)提出同时反演弹性参数的背景模型和扰模型来达到拟合多分量反射波数据的目的。Wang et al.(2018)联合地面P/S分离+DIW获取不同波模式的走时残差，通过空间波场波场解耦梯度预条件来压制参数串扰现象；并在此基础上，设计出稳健的“两步法”弹性反射波走时反演流程。
[0004]除了声波RWI面临的问题，弹性反射波波形反演(E
‑
RWI)主要还面临以下困难：1)数据域目标泛函需要同时匹配P、S两种波模式，若初始模型浅层不够精确，伴随源中容易出现波模式的转换现象；2)P/S波模式转换、波路径重叠等使得不同物理参数在梯度上存在严重的串扰，导致反演结果出现多参数耦合问题。3)复杂介质中横波初始速度模型不容易构建。
[0005]因此，期待一种反演速度模型建模方法，可以解决上述问题。
[0006]公开于本专利技术
技术介绍
部分的信息仅仅旨在加深对本专利技术的一般
技术介绍
的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
技术实现思路
[0007]本专利技术的目的是提出一种反演速度模型建模方法、装置、电子设备及介质，实现降低反演中的非线性程度，提高反演精度与收敛效率。
[0008]第一方面，本公开实施例提供了一种反演速度模型建模方法，包括：
[0009]步骤1：将多分量地震记录进行P/S分离，获得单分量的P波数据和S波数据，并进一步获得标量的PP反射波数据和标量的PS反射波数据，同时构建反演初始的纵波速度模型和
横波速度模型；
[0010]步骤2：基于所述PP反射波数据利用反射波波动方程走时反演，更新当前的纵波速度模型，完成纵波速度模型低波数建模；
[0011]步骤3：基于所述PS反射波数据利用反射波波动方程走时反演，更新当前的横波速度模型，完成横波速度模型低波数建模；
[0012]步骤4：将所述步骤2和所述步骤3得到的模型作为背景模型，利用弹性反射波波形反演，完成纵波速度模型中波数建模和横波速度模型中波数建模。
[0013]可选方案中，所述步骤2包括：
[0014]步骤201：输入所述PP反射波数据与初始的纵波速度模型；
[0015]步骤202：根据当前的纵波速度模型，利用声波方程一阶Born正演获取模拟PP波数据，计算出所述模拟的PP波数据与所述PP反射波数据的走时残差；
[0016]步骤203：根据所述走时残差计算纵波速度模型的梯度方向，以更新所述纵波速度模型；
[0017]步骤204：重复所述步骤202，判断所述走时残差是否收敛，若收敛，实现纵波速度模型低波数建模。
[0018]可选方案中，所述步骤3包括：
[0019]步骤301：输入所述PS反射波数据与初始的横波速度模型；
[0020]步骤302：根据当前的横波速度模型，利用声波方程一阶Born正演获取模拟PS波数据，计算出所述模拟的PS波数据与所述PS反射波数据的走时残差；
[0021]步骤303：根据所述走时残差计算横波速度模型的梯度方向，以更新所述横波速度模型；
[0022]步骤304：重复所述步骤302，判断所述走时残差是否收敛，若收敛，实现横波速度模型低波数建模。
[0023]可选方案中，所述步骤4包括：
[0024]步骤401：输入多分量反射波数据和所述步骤2、所述步骤3走时反演得到的所述背景模型；
[0025]步骤402：根据当前的背景模型，利用弹性波Born正演获取模拟的多分量反射数据，计算出输入的所述多分量反射波数据和模拟的多分量反射数据的波形残差；
[0026]步骤403：根据所述波形残差计算出纵波速度模型和横波速度模型的梯度方向，更新当前背景模型；
[0027]步骤404：重复步骤402，判断所述波形残差是否收敛，若收敛，实现纵波速度模型中波数建模和横波速度模型中波数建模。
[0028]可选方案中，所述步骤2中，所述反射波波动方程走时反演的目标泛函为：
[0029][0030]其中，τ(t)表示走时残差；和分别表示模拟的反射波数据和标量PP反射波数据，E表示为反射波波动方程走时反演的目标函数。
[0031]可选方案中，所述弹性反射波波形反演的目标泛函为：
[0032](m
b
)
*
＝arg min C(m
b
)
[0033]其中代表二范数下的目标函数,代表二范数下的目标函数,代表输入的多分量反射数据与模拟的多分量反射数据的波形残差，为采样算子。
[0034]可选方案中，所述步骤403中的纵波速度模型的梯度表达式为：
[0035][0036]其中，α表示纵波速度；
[0037]横波速度模型的梯度表达式为：
[0038][0039]其中，β表示横波速度。
[0040]可选方案中，所述步骤203中，计算纵波速度模型梯度方向的公式为：
[0041][0042]其中，为时间域反射波场对背景模型的Frechet微分，为反射波走时反演的伴随源，为反射波走时反演的伴随源，表示d的时间导数，E表示为反射波波动方程走时反演的目标函数，和分别表示模拟的反射波数据和标量PP反射波数据。
[0043]第二方面，本公开实施例还提供了一种电子设备，包括：
[0044]至少一个处理器；以及，
[0045]与所述至少一个处理器通信连接的存储器；其中，
[0046]所述存储器存储有可本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.一种反演速度模型建模方法，其特征在于，包括：步骤1：将多分量地震记录进行P/S分离，获得单分量的P波数据和S波数据，并进一步获得标量的PP反射波数据和标量的PS反射波数据，同时构建反演初始的纵波速度模型和横波速度模型；步骤2：基于所述PP反射波数据利用反射波波动方程走时反演，更新当前的纵波速度模型，完成纵波速度模型低波数建模；步骤3：基于所述PS反射波数据利用反射波波动方程走时反演，更新当前的横波速度模型，完成横波速度模型低波数建模；步骤4：将所述步骤2和所述步骤3得到的模型作为背景模型，利用弹性反射波波形反演，完成纵波速度模型中波数建模和横波速度模型中波数建模。2.根据权利要求1所述的反演速度模型建模方法，其特征在于，所述步骤2包括：步骤201：输入所述PP反射波数据与初始的纵波速度模型；步骤202：根据当前的纵波速度模型，利用声波方程一阶Born正演获取模拟PP波数据，计算出所述模拟的PP波数据与所述PP反射波数据的走时残差；步骤203：根据所述走时残差计算纵波速度模型的梯度方向，以更新所述纵波速度模型；步骤204：重复所述步骤202，判断所述走时残差是否收敛，若收敛，实现纵波速度模型低波数建模。3.根据权利要求1所述的反演速度模型建模方法，其特征在于，所述步骤3包括：步骤301：输入所述PS反射波数据与初始的横波速度模型；步骤302：根据当前的横波速度模型，利用声波方程一阶Born正演获取模拟PS波数据，计算出所述模拟的PS波数据与所述PS反射波数据的走时残差；步骤303：根据所述走时残差计算横波速度模型的梯度方向，以更新所述横波速度模型；步骤304：重复所述步骤302，判断所述走时残差是否收敛，若收敛，实现横波速度模型低波数建模。4.根据权利要求1所述的反演速度模型建模方法，其特征在于，所述步骤4包括：步骤401：输入多分量反射波数据和所述步骤2、所述步骤3走时反演得到的所述背景模型；步骤402：根据当前的背景模型，利用弹性波Born正演获取模拟的多分量反射数据，计算出输入的所述多分量反射波数据和模拟的多分量反射数据的波形残差；步骤403：根据所述波形残差计算...
【专利技术属性】
技术研发人员：徐文才，胡光辉，孙思雨，何兵红，
申请(专利权)人：中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。