一种基于PDE的重磁联合反演方法技术

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

【技术实现步骤摘要】
一种基于PDE的重磁联合反演方法
本专利技术涉及地球物理探测
，尤其涉及一种基于PDE的重磁联合反演方法。
技术介绍
根据物理场理论，运用定量描述，由已知的异常体性质求出重力场/磁场分布的过程称为正演问题；反之，由重力场/磁场异常求异常体密度、磁性参数、几何参数的过程被称为反演问题。对于反演问题进行求解时，必须建立在正演问题给出的势场理论数学表达式基础上才能进行。文献“LiYG,OldenburgDW.3-Dinversionofmagneticdata[J].Geophysics,1996,61(2):394-408.”和文献“LiY,OldenburgDW.3-Dinversionofgravitydata[J].Geophysics,1998,63(1):109-119.”提出了一种通过积分方程对结构化六面体(立方体)网格单元分别进行磁场反演与重力场反演的方法。文献“Li,Yaoguo.3-Dinversionofgravitygradiometerdata[J].SegTechnicalProgramExpandedAbstracts,2001:2135.”提出了一种通过重力梯度张量数据进行磁场三维反演的方法。“PeterG.Lelièvre,OldenburgDW.Magneticforwardmodellingandinversionforhighsusceptibility[J].GeophysicalJournalInternational,2006.”提出了一种通过PDE(偏微分方程)进行有限体积的结构化立方体网格的正反演方法。文献“ThomasGünther,CarstenRücker,KlausSpitzer.Three-dimensionalmodellingandinversionofdcresistivitydataincorporatingtopography–II.Inversion[J].GeophysicalJournalInternational,2006,166(2).”提出了基于非结构化网格的电法反演方法，未见基于非结构化网格的三维重磁反演；文献“FarquharsonC.G.,C.R.W.Mosher,2009,Three-dimensionalmodellingofgravitydatausingfinitedifferences,JournalofAppliedGeophysics,68,417-422.”实现了六面体结构化网格的三维差分正演，未实现反演。文献“JahandariH.andC.G.Farquharson,2013,Forwardmodelingofgravitydatausingfinite-volumeandfinite-elementmethodsonunstructuredgrids,Geophysics,78(3),G69–G80”实现了四面体非结构化网格的有限体积、有限元的重力三维正演，未实现反演。文献“GalleyC.G.,P.G.Lelièvre_andC.G.Farquharson,2020,Geophysicalinversionfor3Dcontactsurfacegeometry,Geophysics,85,K27-K45.”采用积分方程实现了非结构化网格的曲面积分反演，与PDE方法在基础理论与方程上显然不同。文献“LuisA.Gallardo,M.A.Pérez-Flores,E.Gómez-Refinementofthree-dimensionalmultilayermodelsofbasinsandcrustalenvironmentsbyinversionofgravityandmagneticdata.2004,397(1):37-54”实现了重磁数据联合反演。文献“GallardoLA.Multiplecross-gradientjointinversionforgeospectralimaging[J].GeophysicalResearchLetters,2007,34(19):5.”实现了电法跟磁法的交叉梯度联合反演，未见采用PDE方法的结构化和非结构的交叉梯度联合反演。然而这些方法存在以下问题：1)结构化六面体(立方体)网格对联合反演模型空间和异常体的剖分精度与效率存在条件限制：一方面，采用立方体单元在对复杂三维结构进行剖分时，容易导致模型离散化误差过大，拟合精度低；另一方面，若采用大量精细立方体单元对复杂三维结构进行剖分，则会导致三维正反演计算量过大，对硬件性能及计算时间要求较高；2)积分方程方法所构建的计算矩阵为非稀疏矩阵，因此在联合反演中求解复杂度高，难以实现大规模并行化计算。
技术实现思路
针对以上技术问题，本专利技术提供了一种基于PDE的重磁联合反演方法，利用现有的网格生成方法构建非结构化四面体网格或结构化六面体网格，并基于非线性PDE(PartialDifferentialEquations，偏微分方程)正反演理论框架，进行三维有限体积、三维标量有限元或三维矢量有限元方法的三维联合正反演计算。本方法所涉及的各个计算矩阵全部为稀疏矩阵，具备并行化计算条件，为大规模重磁数据处理分析提供了基础。本专利技术提出的一种基于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.一种基于PDE的重磁联合反演方法，其特征在于：具体包括以下步骤：/nS1、给定反演初始模型和反演网格边界，利用已知的网格生成方法生成非结构化四面体网格或结构化六面体网格，作为反演网格空间；/nS2、根据反演初始模型，以及物理场的背景场参数进行PDE重力场与磁场三维正演计算，得到物理场数据d；所述物理场包括重力场和磁场；/nS3、根据物理场数据d,以及与S2中三维正演计算相对应的反演计算，进行反演磁化率模型的PDE三维联合反演迭代计算，迭代过程中，根据迭代结果，对当前反演模型周围的反演网格空间进行加密剖分动态调整；/nS4、当三维联合反演目标函数求解误差达到预先设定误差要求时，输出最终反演模型，否则返回S3继续进行三维联合反演迭代计算。/n

【技术特征摘要】
1.一种基于PDE的重磁联合反演方法，其特征在于：具体包括以下步骤：
S1、给定反演初始模型和反演网格边界，利用已知的网格生成方法生成非结构化四面体网格或结构化六面体网格，作为反演网格空间；
S2、根据反演初始模型，以及物理场的背景场参数进行PDE重力场与磁场三维正演计算，得到物理场数据d；所述物理场包括重力场和磁场；
S3、根据物理场数据d,以及与S2中三维正演计算相对应的反演计算，进行反演磁化率模型的PDE三维联合反演迭代计算，迭代过程中，根据迭代结果，对当前反演模型周围的反演网格空间进行加密剖分动态调整；
S4、当三维联合反演目标函数求解误差达到预先设定误差要求时，输出最终反演模型，否则返回S3继续进行三维联合反演迭代计算。


2.如权利要求1所述的一种基于PDE的重磁联合反演方法，其特征在于：步骤S2中，具体为：根据反演初始模型，采用三维有限体积、三维标量有限元或三维矢量有限元方法，基于偏微分方程PDE的物理场的三维正演公式，...

【专利技术属性】
技术研发人员：左博新，
申请(专利权)人：中国地质大学武汉，
类型：发明
国别省市：湖北;42