欢迎来到知嘟嘟! 联系电话:13095918853 卖家免费入驻,海量在线求购! 卖家免费入驻,海量在线求购!
知嘟嘟
我要发布
联系电话:13095918853
知嘟嘟经纪人
收藏
专利号: 2020115917574
申请人: 中国地质大学(武汉)
专利类型:发明专利
专利状态:已下证
专利领域: 计算;推算;计数
更新日期:2024-01-05
缴费截止日期: 暂无
价格&联系人
年费信息
委托购买

摘要:

权利要求书:

1.一种非结构化网格的重磁数据三维正反演方法,其特征在于:具体包括以下步骤:S1、给定反演初始模型和反演网格边界,利用现有的四面体网格生成方法生成非结构化四面体网格,作为反演网格空间;

S2、根据反演初始模型,以及物理场的背景场参数进行PDE三维正演计算,得到物理场数据d;所述物理场为重力场或者磁场;

S3、根据物理场数据d,以及与S2中正演计算方法相对应的反演计算,进行反演磁化率模型的PDE三维反演迭代计算,迭代过程中,根据迭代结果,对当前反演模型周围的所述反演网格空间进行加密剖分动态调整;

S4、当反演目标函数求解误差达到预先设定误差要求时,输出最终反演模型,否则返回S3继续进行三维反演迭代计算。

2.如权利要求1所述的一种非结构化网格的重磁数据三维正反演方法,其特征在于:步骤S2中,具体为:根据反演初始模型,采用三维有限体积、三维标量有限元或三维矢量有限元方法,基于偏微分方程PDE的物理场的正演公式,计算物理场数据d。

3.如权利要求2所述的一种非结构化网格的重磁数据三维正反演方法,其特征在于:基于偏微分方程PDE的物理场的正演公式中,对于重力场,正演公式为: 对于磁场,正演公式为: 其中,γ表示万有引力常数,ρ表示非结构化网格模型的密度,μ=μ0(1+χ),μ0表示真空磁导率,χ表示非结构化网格模型的磁化率,φ表示重力势或磁势。

4.如权利要求3所述的一种非结构化网格的重磁数据三维正反演方法,其特征在于:若所述物理场为重力场,则所述物理场数据d包括:重力异常数据、重力三分量异常数据和重力梯度张量数据;所述物理场为磁场,则所述物理场数据d包括:磁总场异常数据、磁三分量异常数据和磁梯度张量数据。

5.如权利要求4所述的一种非结构化网格的重磁数据三维正反演方法,其特征在于:所述磁三分量异常数据,表示为式(1):式(1)中,Bs为磁三分量异常数据的向量形式; 为磁三分量异常数据的三个分量;B0表示背景场,即地磁场; 表示地磁场的三个分量;

所述磁总场异常数据,表示为式(2):所述磁梯度张量数据表示为式(3):

6.如权利要求4所述的一种非结构化网格的重磁数据三维正反演方法,其特征在于:所述重力三分量异常数据,表示为式(4):x y z

其中,g、g和g分别表示重力异常的三个分量的数据;

所述重力异常数据表示为式(5):所述重力梯度张量数据表示为式(6):

7.如权利要求1所述的一种非结构化网格的重磁数据三维正反演方法,其特征在于:步骤S3中,反演磁化率模型的PDE三维反演迭代计算的目标函数如式(7):其中,表示三维反演目标函数,采用非线性优化求解;d0表示已有的观测物理场数据;

d=F(B0,m),m≥0,表示正演计算得到的与d0相应物理场数据;F(·)表示PDE三维正演计算;

m表示反演磁化率模型;φreg表示规整化函数。

8.如权利要求7所述的一种非结构化网格的重磁数据三维正反演方法,其特征在于:所述规整化函数,实用深度规整化、最小体积或者最紧支撑函数的任意一种。