强各向异性层状VTI介质高效射线追踪方法、设备和介质与流程

本发明涉及地球物理正演，具体涉及一种强各向异性层状vti介质高效射线追踪方法、设备和介质。

背景技术：

1、各向异性对地震波传播的振幅和传播时间有明显的影响。层状vti介质广泛应用于地震叠前时间偏移和页岩水力压裂中的微震监测。它们的成像质量关键取决于正演计算的准确性。两点射线追踪方法，如射击法和弯曲法，已被证实是一种实用的方法来寻找射线路径和旅行时间。与各向同性介质不同，地震波在vti介质中会分裂为qp波、qsv波和qsh波，三种波均需要精确计算它们的群速度旅行时间和射线路径。在各向异性介质中，程函方程依赖于相速度角，而在射线路径依赖于群速度角。计算过程中很难找到群速度角与其对应的群速度矢量之间的确切关系。一些近似公式，如byun近似式和sena近似式，在弱各向异性介质中表现良好，而在强各向异性地层中误差较大。另一种选择是预先在群速度和群角之间建立字典，然而，计算过程中查找字典中的大量信息大大降低了计算效率。因此，传统方法在计算强各向异性vti介质射线追踪时，存在计算精度低、耗时大的问题，不适用于工业生产应用。研究表明，慢度系数可以在qp波和qsv波中形成四次曲面，而在qsh波中可以形成二次曲面。vti介质中的群速度有其射线参数的解析式，适用于强各向异性介质。因此，与传统方法利用角度作为自变量不同，可以利用慢度作为自变量计算群速度，进而可以推导出射线路径和走时关于慢度的解析公式，最终通过高斯牛顿法快速迭代获得慢度大小，推导获得各层的射线路径和走时大小。相对于传统方法，基于慢度为自变量的方程适用于任意强度各向异性介质，计算效率更高，结果更加精确。

2、因此，现有技术还有待于改进和发展。

技术实现思路

1、针对现有技术中的不足，本发明提供一种强各向异性层状vti介质高效射线追踪方法、设备和介质，其旨在解决基于现有技术计算vti介质射线路径和走时精度差、效率低的问题。

2、为实现上述目的，本发明可以采用以下技术方案进行：

3、第一方面，本发明提供一种强各向异性层状vti介质高效射线追踪方法，其包括如下步骤：

4、基于克里斯托费尔频散理论推导vti介质中高频近似假设条件下程函方程qp波、qsv波和qsh波慢度平面表达公式；

5、通过慢度平面推导水平方向群速度vx和垂直方向群速度vz的解析表达式；

6、通过vti介质水平方向群速度vx和垂直方向群速度vz推导qp波、qsv波和qsh波的射线路径和走时表达式；

7、通过参数变换，将原慢度参数为自变量的射线路径表达式转换为单层最大水平路径为自变量的射线路径表达式，推导在已知射线水平路径条件下目标函数表达式，并通过高斯牛顿法推导反演求解迭代公式；

8、通过寻求离散值拟合射线路径中水平分量路径曲线，推导反演中的初始值大小；

9、通过水平射线路径公式推导出最大水平路径取值，根据最大水平路径取值推导出该射线路径和走时；

10、综合射线路径和走时即可确定该次射线追踪结果。

11、第二方面，本发明提供一种电子设备，所述电子设备包括处理器和存储器，所述存储器中存储有至少一条指令、至少一段程序、代码集或指令集，所述至少一条指令、所述至少一段程序、所述代码集或指令集由所述处理器加载并执行，以实现如上所述的强各向异性层状vti介质高效射线追踪方法。

12、第三方面，本发明提供一种计算机可读存储介质，所述存储介质中存储有至少一条指令、至少一段程序、代码集或指令集，所述至少一条指令、所述至少一段程序、所述代码集或指令集由处理器加载并执行以实现如上所述的强各向异性层状vti介质高效射线追踪方法。

13、本发明与现有技术相比，其有益效果在于：本发明实施例的一种强各向异性层状vti介质高效射线追踪方法，其基于克里斯托费尔频散理论推导vti介质中高频近似假设条件下程函方程qp波、qsv波和qsh波慢度平面表达公式；通过慢度平面推导水平方向群速度vx和垂直方向群速度vz的解析表达式；通过vx和vz推导qp波、qsv波和qsh波的射线路径和走时表达式；通过参数变换转换为单层最大水平路径为自变量的射线路径表达式；通过高斯牛顿法推导反演求解迭代公式；通过离散值拟合射线路径推导反演中的初始值大小。本方法的优势在于可以高效、高精度地获得任意各向异性vti介质中的旅行时间和射线路径。在强各向异性vti介质中，直达波和反射波的射线路径是精确的。该方法的提出，有利于提高地震叠前时间偏移和页岩水力压裂中微震监测正演计算的准确性，为油气资源的低成本高效开采提供了良好的技术基础。

技术特征：

1.一种强各向异性层状vti介质高效射线追踪方法，其特征在于，包括如下步骤：

2.根据权利要求1所述的强各向异性层状vti介质高效射线追踪方法，其特征在于，基于克里斯托费尔频散理论推导vti介质中高频近似假设条件下程函方程qp波、qsv波和qsh波慢度平面表达公式，具体步骤包括：

3.根据权利要求1所述的强各向异性层状vti介质高效射线追踪方法，其特征在于，通过慢度平面推导水平方向群速度vx和垂直方向群速度vz的解析表达式，具体步骤包括：

4.根据权利要求1所述的强各向异性层状vti介质高效射线追踪方法，其特征在于，通过vti介质水平方向群速度vx和垂直方向群速度vz推导qp波、qsv波和qsh波的射线路径和走时表达式，具体步骤包括：

5.根据权利要求1所述的强各向异性层状vti介质高效射线追踪方法，其特征在于，通过参数变换，将原慢度参数为自变量的射线路径表达式转换为单层最大水平路径为自变量的射线路径表达式，具体步骤包括：

6.根据权利要求1所述的强各向异性层状vti介质高效射线追踪方法，其特征在于，推导在已知射线水平路径条件下目标函数表达式，具体步骤包括：

7.根据权利要求1所述的强各向异性层状vti介质高效射线追踪方法，其特征在于，通过寻求离散值拟合射线路径中水平分量路径曲线，推导反演中的初始值大小，具体步骤包括：

8.根据权利要求1所述的强各向异性层状vti介质高效射线追踪方法，其特征在于，

9.一种电子设备，其特征在于，所述电子设备包括处理器和存储器，所述存储器中存储有至少一条指令、至少一段程序、代码集或指令集，所述至少一条指令、所述至少一段程序、所述代码集或指令集由所述处理器加载并执行，以实现如权利要求1至8任一所述的强各向异性层状vti介质高效射线追踪方法。

10.一种计算机可读存储介质，其特征在于，所述存储介质中存储有至少一条指令、至少一段程序、代码集或指令集，所述至少一条指令、所述至少一段程序、所述代码集或指令集由处理器加载并执行以实现权利要求1至8任一所述的强各向异性层状vti介质高效射线追踪方法。

技术总结
本发明公开了强各向异性层状VTI介质高效射线追踪方法、设备和介质，涉及地球物理正演技术领域，该方法包括：基于克里斯托费尔频散理论推导VTI介质中高频近似假设条件下程函方程qP波、qSV波和qSH波慢度平面表达公式；推导水平方向群速度V<subgt;x</subgt;和垂直方向群速度V<subgt;z</subgt;的解析表达式；推导qP波、qSV波和qSH波的射线路径和走时表达式；推导在已知射线水平路径条件下目标函数表达式，并通过高斯牛顿法推导反演求解迭代公式；推导反演中的初始值大小；推导出最大水平路径取值及该射线路径和走时；综合射线路径和走时即可确定该次射线追踪结果。本发明解决基于现有技术计算VTI介质射线路径和走时精度差、效率低的问题。

技术研发人员：蒋星达,杨华勇,李超,杨得厚,李跃金
受保护的技术使用者：南方海洋科学与工程广东省实验室（广州）
技术研发日：
技术公布日：2024/1/13