反演速度模型建模方法、装置、电子设备及介质.pdf

(19)国家知识产权局(12)发明专利申请(10)申请公布号CN116029074A(43)申请公布日2023.04.28(21)申请号202111240920.7(22)申请日2021.10.25(71)申请人中国石油化工股份有限公司地址100027北京市朝阳区朝阳门北大街22号申请人中国石油化工股份有限公司石油物探技术研究院(72)发明人徐文才胡光辉孙思雨何兵红(74)专利代理机构北京思创毕升专利事务所11218专利代理师孙向民廉莉莉(51)Int.Cl.G06F30/20(2020.01)G01V1/28(2006.01)权利要求书3页说明书11页附图10页(54)发明名称反演速度模型建模方法、装置、电子设备及介质(57)摘要本申请公开了一种反演速度模型建模方法、装置、电子设备及介质，其中，该方法包括：步骤1：将多分量地震记录进行P/S分离，获得单分量的P波数据和S波数据，并进一步获得标量的PP反射波数据和标量的PS反射波数据，同时构建反演初始的纵波速度模型和横波速度模型；步骤2：基于PP反射波数据利用反射波波动方程走时反演，更新当前的纵波速度模型，完成纵波速度模型低波数建模；步骤3：基于PS反射波数据利用反射波波动方程走时反演，更新当前的横波速度模型，完成横波速度模型低波数建模；步骤4：将步骤2和步骤3得到的模型作为背景模型，利用弹性反射波波形反演，完成纵波速度模型中波数建模和横波速度模型中波数建模。CN116029074ACN116029074A权利要求书1/3页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：根据所述波形残差计算出纵波速度模型和横波速度模型的梯度方向，更新当前背景模型；步骤404：重复步骤402，判断所述波形残差是否收敛，若收敛，实现纵波速度模型中波数建模和横波速度模型中波数建模。5.根据权利要求1所述的反演速度模型建模方法，其特征在于，所述步骤2中，所述反射波波动方程走时反演的目标泛函为：2CN116029074A权利要求书2/3页其中，τ(t)表示走时残差；和分别表示模拟的反射波数据和标量PP反射波数据，E表示为反射波波动方程走时反演的目标函数。6.根据权利要求1所述的反演速度模型建模方法，其特征在于，所述弹性反射波波形反演的目标泛函为：*(mb)＝argminC(mb)其中代表二范数下的目标函数，代表输入的多分量反射数据与模拟的多分量反射数据的波形残差，为采样算子。7.根据权利要求4所述的反演速度模型建模方法，其特征在于，所述步骤403中的纵波速度模型的梯度表达式为：其中，α表示纵波速度；横波速度模型的梯度表达式为：其中，β表示横波速度。8.根据权利要求2所述的反演速度模型建模方法，其特征在于