申请/专利权人：中国地质大学(武汉)

申请日：2021-03-08

公开（公告）日：2024-03-22

公开（公告）号：CN112949133B

主分类号：G06F30/23

分类号：G06F30/23;G06T17/05;G06F119/14

优先权：

专利状态码：有效-授权

法律状态：2024.03.22#授权;2021.07.02#实质审查的生效;2021.06.11#公开

摘要：本发明公开了一种基于PDE的重磁联合反演方法，包括：给定反演初始模型和反演网格边界，利用网格生成方法生成非结构化四面体网格或结构化六面体网格，作为反演网格空间；根据反演初始模型，以及物理场的背景场参数进行PDE重力场与磁场三维正演计算，得到物理场数据d；根据物理场数据d,以及与S2相对应的反演计算，进行反演模型的PDE三维联合反演迭代计算，迭代过程中，根据迭代结果，对当前反演模型周围的反演网格空间进行加密剖分动态调整；当反演目标函数求解误差达到预先设定误差要求时，输出最终反演模型，否则返回S3继续进行三维联合反演迭代计算。本发明对于计算空间中物理场分布、描述几何形状具有良好的适用性，更高的稳定性和精度。

主权项：1.一种基于PDE的重磁联合反演方法，其特征在于：具体包括以下步骤：S1、给定反演初始模型和反演网格边界，利用已知的网格生成方法生成非结构化四面体网格或结构化六面体网格，作为反演网格空间；S2、根据反演初始模型，以及物理场的背景场参数进行PDE重力场与磁场三维正演计算，得到物理场数据d；所述物理场包括重力场和磁场；S3、根据物理场数据d,以及与S2中三维正演计算相对应的反演计算，进行反演磁化率模型的PDE三维联合反演迭代计算，迭代过程中，根据迭代结果，对当前反演模型周围的反演网格空间进行加密剖分动态调整；S4、当三维联合反演目标函数求解误差达到预先设定误差要求时，输出最终反演模型，否则返回S3继续进行三维联合反演迭代计算；步骤S2中，具体为：根据反演初始模型，采用三维有限体积、三维标量有限元或三维矢量有限元方法，基于偏微分方程PDE的物理场的三维正演公式，计算物理场数据d；其中物理场包括重力场和磁场；对于联合反演中的重力场，则所述重力场数据为重力梯度张量数据；对于联合反演中的磁场，所述磁场数据为磁梯度张量数据；步骤S3中，重磁数据的PDE三维联合反演的目标函数如式1： 其中Φd可以表示为如式2： 其中为包含深度权重的权重矩阵；dobs为实际观测数据，参数β1和β2为常数值；Φcg为交叉梯度函数，可以表示为如式3： 其中，m1和m2表示密度和磁化率两种不同的模型参数，D表示联合反演空间；Φm1,m2表示三维联合反演目标函数，采用非线性优化求解；η为动态规整系数，在每一次迭代过程中不断调整η的值，当目标函数最小化达到阈值时，完成非线性优化求解过程。

全文数据：

权利要求：

百度查询： 中国地质大学(武汉) 一种基于PDE的重磁联合反演方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD