[发明专利]水深反演方法、装置和计算机可读存储介质

权利要求书

1.一种水深反演方法，其特征在于，包括：

确定预设的当前密度差；

基于重力异常公式和重力正演模型，确定针对研究区的水深模型的目标函数，所述重力异常公式指示重力异常包括长波参考场部分和短波残差场部分，所述短波残差场部分包括密度差参数；

利用所述目标函数，确定当前水深模型，以便获得多个水深模型，其中在所述当前水深模型中，所述当前密度差作为所述密度差参数；

基于所述多个水深模型，确定最优水深模型；

所述基于重力异常公式和重力正演模型，确定针对研究区的水深模型的目标函数，包括：

获取所述研究区中的多个观测点的位置和重力异常，以及多个反演控制点的位置、先验水深和重力异常；

将所述多个反演控制点的先验水深应用于所述重力正演模型，得到所述多个反演控制点的短波重力异常，以便基于所述重力异常公式与所述多个反演控制点的重力异常作差，获得所述多个反演控制点的长波重力异常；

将插值算法应用于所述多个反演控制点的位置和所述多个观测点的位置，基于所述多个反演控制点的长波重力异常，得到所述多个观测点的长波重力异常，以便基于所述重力异常公式与所述多个观测点的重力异常作差，获得所述多个观测点的短波重力异常；

基于所述多个观测点的短波重力异常，确定所述目标函数。

2.根据权利要求1所述的水深反演方法，其特征在于，所述目标函数包括所述水深模型的稳定泛函和失配泛函，所述失配泛函基于所述多个观测点的短波重力异常确定。

3.根据权利要求2所述的水深反演方法，其特征在于，所述目标函数为F(m)＝F_d(m)+αF_s(m)，其中，所述当前水深模型为m，F_s(m)为稳定泛函，F_d(m)为失配泛函，α为正则化因子。

4.根据权利要求3所述的水深反演方法，其特征在于，f_d(m)＝g(m)-g^O，g(m)为所述重力正演值，g^O为重力观测值，N为观测点个数，T表征转置。

5.根据权利要求4所述的水深反演方法，其特征在于，所述利用所述目标函数，确定所述当前水深模型，包括：

通过迭代计算，确定所述当前水深模型，其中，在所述迭代计算中，每次利用α_k进行迭代，确定所述水深模型，其中，α₀＝0,α_k＝α₁q^k-1,k＝1,2,...,itermax；0q1,α₁＝F_d(m₀)/F_s(m₀),itermax为最大迭代次数，m₀为初始水深模型，q控制迭代速率，且0q1。

6.根据权利要求2-5中任一项所述的水深反演方法，其特征在于，还包括：获取多个检核控制点的位置、先验水深和重力异常，其中，所述基于所述多个水深模型，确定最优水深模型，包括：

将所述多个检核控制点的位置和重力异常应用于所述当前水深模型，得到所述多个检核控制点的反演水深；

将所述多个检核控制的反演水深与所述多个检核控制的先验水深进行比较评估，输出作为对应于所述当前水深模型的比较结果，从而得到多个比较结果；

基于所述多个比较结果，选择最优水深模型。

7.根据权利要求1所述的水深反演方法，其特征在于，所述方法还包括：

确定预设的当前步长；

所述获得多个水深模型，包括：

基于所述当前密度差和所述当前步长，确定下一密度差，其中所述当前步长为恒定值。