【CN110161562A】倾斜横向各向同性介质中的裂缝弱度反演方法及系统【专利】

(19)中华人民共和国国家知识产权局

(12)发明专利申请

(10)申请公布号 (43)申请公布日 (21)申请号 201910354239.1

(22)申请日 2019.04.29

(71)申请人 中国石油大学（华东）

地址 266580 山东省青岛市经济技术开发

区长江西路66号

(72)发明人 潘新朋　张广智　周顺心　李林　

印兴耀　

(74)专利代理机构 济南圣达知识产权代理有限

公司 37221

代理人 李琳

(51)Int.Cl.

G01V 1/30(2006.01)

(54)发明名称倾斜横向各向同性介质中的裂缝弱度反演方法及系统(57)摘要本公开提供了一种倾斜横向各向同性介质中的裂缝弱度反演方法及系统，使用垂直横向各向同性介质和Bond矩阵的刚度矩阵，根据裂缝弱度和倾角推导出倾斜横向各向同性介质的刚度矩阵，结合刚度扰动和散射理论，分离两个弱各向异性倾斜横向各向同性介质的弱对比界面，得到该界面下PP波反射系数的方位角，以倾角作为先验信息，利用部分入射角叠加方位地震数据实现方位角叠前反演，根据反演结果估计贝叶斯框架中的裂缝弱度，合理的估计P波和S波模量和裂缝弱度，为地震波的传播研究和地震灾害预防提

供了技术支持。权利要求书1页 说明书12页 附图1页CN 110161562 A 2019.08.23

C N 110161562

A

权　利　要　求　书1/1页CN 110161562 A

1.一种倾斜横向各向同性介质中的裂缝弱度反演方法，其特征是：包括以下步骤：

使用垂直横向各向同性介质和Bond矩阵的刚度矩阵，根据裂缝弱度和倾角推导出倾斜横向各向同性介质的刚度矩阵，结合刚度扰动和散射理论，分离两个弱各向异性倾斜横向各向同性介质的弱对比界面，得到该界面下PP波反射系数的方位角，以倾角作为先验信息，利用部分入射角叠加方位地震数据实现方位角叠前反演，根据反演结果估计贝叶斯框架中的裂缝弱度。

2.如权利要求1所述的一种倾斜横向各向同性介质中的裂缝弱度反演方法，其特征是：根据法向和切向过剩弱度以及Bond变换对弱各向异性垂直横向各向同性介质的刚度矩阵进行积分，根据法向和切向过剩弱度以及倾角推导倾斜横向各向同性介质的刚度矩阵。

3.如权利要求2所述的一种倾斜横向各向同性介质中的裂缝弱度反演方法，其特征是：考虑由平面水平弱相对界面分离的两个均匀弱各向异性倾斜横向各向同性介质，并且通过使用P波和S波速度作为坐标轴来参数化半空间中的每个倾斜横向各向同性介质，融合密度、法向和切向裂缝弱度参数、以及对称轴和z轴之间的角度差和方位角，确定倾斜横向各向同性介质的刚度矩阵。

4.如权利要求1所述的一种倾斜横向各向同性介质中的裂缝弱度反演方法，其特征是：将散射理论与倾斜横向各向同性介质刚度矩阵的扰动相结合，根据背景P-和S-波模量、密度、裂缝弱度和倾斜角得出线性化PP波反射系数，分离两个弱各向异性倾斜横向各向同性介质的弱对比界面。

5.如权利要求1所述的一种倾斜横向各向同性介质中的裂缝弱度反演方法，其特征是：在贝叶斯框架下，进行方位角的反演，得到贝叶斯框架下的裂缝弱度。

6.如权利要求1所述的一种倾斜横向各向同性介质中的裂缝弱度反演方法，其特征是：利用Cauchy稀疏正则化和低频信息正则化以及迭代加权最小二乘法优化算法的约束来实现方位角的叠前反演，根据反演结果估计裂缝弱度。

7.如权利要求6所述的一种倾斜横向各向同性介质中的裂缝弱度反演方法，其特征是：使用Cauchy函数作为先验分布，使用高斯函数作为似然函数，通过Cauchy先验分布和高斯似然函数联合解决后验分布，得到分布最大化方程，将低频信息正则化与最大化方程相结合，使用迭代重加权最小二乘算法来求解，得到模型参数的最终结果，即反演结果。

8.一种倾斜横向各向同性介质中的裂缝弱度反演系统，其特征是：被配置为运行在服务器或处理器上，执行权利要求1-7中任一项所述的倾斜横向各向同性介质中的裂缝弱度反演方法。

9.一种计算机可读存储介质，其特征是：其中存储有多条指令，所述指令适于由终端设备的处理器加载并执行权利要求1-7中任一项所述的倾斜横向各向同性介质中的裂缝弱度反演方法。

10.一种终端设备，其特征是：包括处理器和计算机可读存储介质，处理器用于实现各指令；计算机可读存储介质用于存储多条指令，所述指令适于由处理器加载并执行权利要求1-7中任一项所述的倾斜横向各向同性介质中的裂缝弱度反演方法。

2