[发明专利]基于包络的各向异性参数反演方法、装置及存储介质

权利要求书

1.基于包络的各向异性参数反演方法，其特征在于，所述方法包括：

获取模拟地震数据以及观测地震数据，通过模拟地震数据以及观测地震数据分别获取模拟地震数据包络以及观测地震数据包络，根据模拟地震数据包络以及观测地震数据包络的互相关计算旅行时差；

基于旅行时差定义模型的目标函数，获取由各向异性参数扰动导致的包络变化，通过各向异性参数扰动导致的包络变化获取由各向异性参数扰动导致的旅行时差的变化；

基于旅行时差以及各向异性参数扰动导致的旅行时差的变化得到包络旅行时目标函数的伴随源；

通过伴随源求解伴随波动方程获得逆向传播波场，通过正演模拟获得正向传播波场；

通过逆向传播波场和正向传播波场的互相关获得各向异性参数的梯度，通过线性搜索算法获取迭代步长，通过L-BFGS算法获取每次迭代的近似海森逆矩阵，通过近似海森逆矩阵以及各向异性参数的梯度获得搜索方向，根据搜索方向以及迭代步长对各向异性参数进行迭代反演，直到满足预设的迭代次数为止，得到最终的各向异性参数。

2.根据权利要求1所述的方法，其特征在于，

所述通过模拟地震数据获取模拟地震数据包络的表达式如下所示：

式中，表示模拟地震数据包络，代表希尔伯特变换，表示模拟地震数据；

所述通过观测地震数据获取观测地震数据包络的表达式如下所示：

式中，表示观测地震数据包络，代表希尔伯特变换，表示观测地震数据。

3.根据权利要求2所述的方法，其特征在于，

所述根据模拟地震数据包络以及观测地震数据包络的互相关计算旅行时差的表达式如下所示：

式中，表示旅行时差，为互相关的时间移动量，为观测地震数据采集的最大记录时间。

4.根据权利要求3所述的方法，其特征在于，

所述基于旅行时差定义模型的目标函数的表达式如下所示：

式中，表示各向异性参数。

5.根据权利要求4所述的方法，其特征在于，

所述获取由各向异性参数扰动导致的包络变化的表达式如下所示：

式中，表示由各向异性参数扰动导致的包络变化。

6.根据权利要求5所述的方法，其特征在于，

所述通过各向异性参数扰动导致的包络变化获取由各向异性参数扰动导致的旅行时差的变化的表达式如下所示：

其中，

式中，表示由各向异性参数扰动导致的旅行时差的变化。

7.根据权利要求6所述的方法，其特征在于，

所述基于旅行时差以及各向异性参数扰动导致的旅行时差的变化得到包络旅行时目标函数的伴随源的表达式如下所示：

式中，表示伴随源，将各向异性参数扰动导致的旅行时差的变化的计算公式代入伴随源计算公式中，得到：

其中，

。

8.根据权利要求7所述的方法，其特征在于，

所述通过伴随源求解伴随波动方程获得逆向传播波场的表达式如下所示：

式中，和分别为弹性系数矩阵和应变张量，为密度，为逆向传播波场。