一种深度学习优化的变密度声波方程全波形反演方法

本发明涉及一种地下介质的岩石物理参数成像方法，具体是一种深度学习优化的变密度声波方程全波形反演方法，属于深度学习与地震勘探速度建模技术的交叉领域。

背景技术：

1、地下介质的岩石物理属性的准确估计，是发现地下油气资源等靶区的关键。全波形反演是一种可以提供高精度物性参数估计的最高成像精度的方法之一。在过去几十年中，传统的全波形反演取得了巨大的成就，形成了一个成熟的计算过程，无论是在计算精度还是效率上都很高。

2、近年来，深度学习优化全波形反演也逐渐开始研究，目前已经研究了带物理约束(导向)的深度学习全波形反演，但它需要大量的存储空间，因为基于自动微分的梯度计算等同于伴随方法的梯度。目前，在基于深度学习的全波形反演中，无法结合过去几十年累积的自定义目标函数，其原因是：目前深度学习中目标函数的研究主要是针对样本数据的图形处理问题，无法直接适用于多解性的地球物理反问题研究。为了有效地降低全波形反演多解性问题，地球物理学领域在过去几十年积淀了多种行之有效的目标函数，可以克服或者缓解全波形反演中根本性的低频依赖性问题，并且显著地提升了参数建模的精度。另外目前基于深度学习的全波形反演无法实现速度和岩石密度的同时建模，导致反演过程复杂，且速度反演时需要大量的网络参数存储，最终使得整个全波形反演的效率及稳定性还有待提高。

技术实现思路

1、针对上述现有技术存在的问题，本发明提供一种深度学习优化的变密度声波方程全波形反演方法，通过利用传统全波形反演高性能计算方案实现梯度的计算，从而适用于全波形反演的各种自定义目标函数，克服或者缓解全波形反演中根本性的低频依赖性问题；同时其引入一阶变密度速度应力方程，实现了速度和密度的同时建模，从而使得反演过程简化，提高整个全波形反演的效率及稳定性。

2、为了实现上述目的，本发明采用的技术方案是：一种深度学习优化的变密度声波方程全波形反演方法，具体步骤为：

3、a、先确定所需反演的区域，然后采用观测系统获取观测数据；

4、b、构建正传波场方程，所述正传波场方程为二维一阶变密度声波方程；

5、c、构建卷积目标函数，所述卷积目标函数具有良好鲁棒性；

6、d、根据步骤b的正传波场方程和步骤c的卷积目标函数，基于伴随理论确定该卷积目标函数的反传波场方程；

7、e、利用步骤b的正传波场和步骤d的反传波场构建卷积目标函数中模型参数的梯度，所述梯度由速度的梯度和密度的梯度组成；

8、f、将步骤e中模型参数的梯度设置为可训练的变量，然后将可训练的变量输入深度框架内，并进入步骤h进行反演；

9、g、采用深度学习优化方法对步骤f获得的深度框架进行模型参数优化，并将步骤e中计算的梯度输入深度框架中，通过训练深度框架的方式实现全波形反演，最终获得反演后的速度反演结果和密度反演结果。

10、进一步，所述步骤b中二维一阶变密度声波方程具体为：

11、

12、其中，ψ＝[σ,v]t是正传波场变量，σ是正传波场对应的质点震动应力，v＝[vx,vz]表示质点沿着x和z两个方向的振动速度，s是源项，表示对时间t求偏导数，t表示矩阵的转置，链接矩阵ρ表示介质的密度，v表示介质的速度，和表示对x和z求偏导数。

13、进一步，所述步骤c中卷积目标函数具体为：

14、

15、其中，d是观测数据，v是合成波场数据，xr和xref分别表示检波点位置和参考道的位置，空间坐标x＝[x,z]，模型参数m＝[v,ρ]t。

16、进一步，所述步骤d中反传波场方程具体为：

17、

18、其中反传波场其中是反传波场对应的质点震动应力，和表示质点沿着x和z两个方向的振动速度，是反传源项，*表示卷积运算，表示互相关运算。

19、进一步，所述步骤e中模型参数的梯度具体为：

20、

21、其中是速度的梯度，是密度的梯度。

22、进一步，所述步骤f中模型参数的梯度设置为可训练的变量具体为：将模型参数的梯度通过公式设置为可训练的变量，并指定初始速度和密度值及梯度，具体公式为：

23、

24、其中，vini和ρini分别表示初始的速度和密度，接着将可训练的变量输入至深度卷积神经网络模型。

25、进一步，采用如下方法替换步骤f：基于步骤e中模型参数的梯度建立深度卷积神经网络模型，所述深度卷积神经网络模型包括速度模型参数和密度模型参数，并进入步骤g进行反演；

26、其中深度卷积神经网络模型具体为：

27、

28、其中，模型参数γ＝v orρ，用以表示模型参数γ的深度卷积神经网络模型，θ表示网络参数，l表示神经网络的总层数，wl和bl表示卷积层或者全连接层的权重和偏置，其中下标l∈[1,l],l∈z，λ表示网络输入的随机特征向量。

29、进一步，所述步骤g中反演的具体过程为：首先将观测数据作为随机特征变量λ，使用深度卷积神经网络模型表征模型参数γ，传入外部的模型参数梯度通过深度学习优化方法进行优化迭代，修改深度卷积神经网络模型参数θγ，进而表示出更新的模型参数γ，并通过高性能计算获取模型参数的梯度再次优化修正深度卷积网络模型参数θγ，直至迭代次数达到设置次数或目标函数值满足精度要求，输出深度卷积网络模型的表征结果，即为最终反演的模型参数γ。

30、进一步，所述步骤g中深度学习优化方法为rmsprop、adagrad、asgd和adam的其中一种。

31、与现有技术相比，本发明具有如下优点：

32、(1)本发明直接利用传统全波形反演高性能计算方案实现梯度的计算，这些梯度可以对应任何适用于全波形反演的自定义目标函数，通过将外部计算的梯度输入到深度卷积神经网络模型中，以实现多样化的优化算法，并通过优化方法的迭代优化反演代表模型参数的网络参数，从而可以克服或者缓解全波形反演中根本性的低频依赖性问题，并且显著地提升了参数建模的精度。

33、(2)本发明进行全波形反演过程中与现有的全波形反演区别在于，本发明中的深度卷积神经网络模型主要用于描述真实速度与初始速度之间的差异，并不直接表示初始模型参数；即本发明中初始速度的加载采用直接与网络输出叠加的方式，后续进行训练时能有效减少网络预训练次数，同时提高了梯度计算中正演模拟的稳定性。

34、(3)本发明利用一阶变密度速度应力方程作为正传波场，然后构建卷积目标函数，并确定相应的反传波场，最后使用深度学习优化实现了基于卷积目标函数的全波形反演，并在一阶变密度速度应力方程下执行了多尺度反演，实现了速度和密度的同时建模，由于本发明直接采用了一阶变密度速度应力方程，该方程不仅具有高模拟精度，而且还可以自适应地考虑地下介质中密度变化，能为地下介质提供直接的速度和密度的物理信息，从而使得反演过程简化，提高整个全波形反演的效率及稳定性。

技术特征：

1.一种深度学习优化的变密度声波方程全波形反演方法，其特征在于，具体步骤为：

2.根据权利要求1所述深度学习优化的变密度声波方程全波形反演方法，其特征在于，所述步骤b中二维一阶变密度声波方程具体为：

3.根据权利要求2所述深度学习优化的变密度声波方程全波形反演方法，其特征在于，所述步骤c中卷积目标函数具体为：

4.根据权利要求3所述深度学习优化的变密度声波方程全波形反演方法，其特征在于，所述步骤d中反传波场方程具体为：

5.根据权利要求4所述深度学习优化的变密度声波方程全波形反演方法，其特征在于，所述步骤e中模型参数的梯度具体为：

6.根据权利要求5所述深度学习优化的变密度声波方程全波形反演方法，其特征在于，所述步骤f中模型参数的梯度设置为可训练的变量具体为：将模型参数的梯度通过公式设置为可训练的变量，并指定初始速度和密度值及梯度，具体公式为：

7.根据权利要求5所述深度学习优化的变密度声波方程全波形反演方法，其特征在于，采用如下方法替换步骤f：基于步骤e中模型参数的梯度建立深度卷积神经网络模型，所述深度卷积神经网络模型包括速度模型参数和密度模型参数，并进入步骤g进行反演；

8.根据权利要求7所述深度学习优化的变密度声波方程全波形反演方法，其特征在于，所述步骤g中反演的具体过程为：首先将观测数据作为随机特征变量λ，使用深度卷积神经网络模型表征模型参数γ，传入外部的模型参数梯度通过深度学习优化方法进行优化迭代，修改深度卷积神经网络模型参数θγ，进而表示出更新的模型参数γ，并通过高性能计算获取模型参数的梯度再次优化修正深度卷积网络模型参数θγ，直至迭代次数达到设置次数或目标函数值满足精度要求，输出深度卷积网络模型的表征结果，即为最终反演的模型参数γ。

9.根据权利要求1所述深度学习优化的变密度声波方程全波形反演方法，其特征在于，所述步骤g中深度学习优化方法为rmsprop、adagrad、asgd和adam的其中一种。

技术总结
本发明公开了一种深度学习优化的变密度声波方程全波形反演方法，通过直接利用传统全波形反演高性能计算方案实现梯度的计算，这些梯度能对应任何适用于全波形反演的自定义目标函数，并与深度学习网络相结合，从而可以克服或者缓解全波形反演中根本性的低频依赖性问题，并且显著地提升参数建模的精度；另外本发明直接引入一阶变密度速度应力方程作为正传波场，其不仅具有高模拟精度，而且还可以自适应地考虑地下介质中密度变化，能为地下介质提供直接的速度和密度的物理信息，实现了速度和密度的同时建模，并采用优化方法进行迭代优化，从而使得反演过程简化，并在一阶变密度速度应力方程下执行多尺度反演，最终提高了整个全波形反演的效率及稳定性。

技术研发人员：方金伟,李娟娟,章俊,刘盛东
受保护的技术使用者：中国矿业大学
技术研发日：
技术公布日：2024/1/15