一种基于积分方程的退磁效应下磁测数据的正反演方法技术

本发明专利技术提供了一种基于积分方程的退磁效应下磁测数据的正反演方法，包括以下步骤：根据地下目标区域，确定网格剖分的空间范围，使用正交网格对所述空间范围进行网格剖分；设定地磁场的磁倾角、磁偏角和磁感应强度；设定观测位置，利用基于积分方程的在退磁效应下的磁测数据正演计算方法，能够获得考虑退磁效应下在观测位置的磁测数据；输入磁场观测数据，通过基于积分方程在退磁效应下的磁测数据反演方法，结合正则化项进行模型最优化计算，对地下目标进行反演计算，得到反演结果。本发明专利技术能计算考虑退磁效应下的磁测数据正演，并能够根据磁测数据进行考虑退磁效应的反演。

【技术实现步骤摘要】
一种基于积分方程的退磁效应下磁测数据的正反演方法
本专利技术涉及地球物理探测领域，尤其涉及一种基于积分方程的退磁效应下磁测数据的正反演方法。
技术介绍
在地磁探测中，人们通常为了减小计算的时间，将异常体的磁化强度进行近似，这样探测出来的异常体的磁异常数据往往与实际的数据解释存在一定的误差。在低磁背景下，这种误差可以忽略不计，但在高磁背景，该误差对磁场异常探测有很大的影响，不考虑退磁效应的处理结果精度较低，无法满足实际应用的需求。因此考虑退磁效应下的磁梯度张量数据不仅能很好的描述异常体的磁化方向、空间几何结构、地质体的走向等信息，而且在高磁背景下精度准确，能很好地解释地下地质信息。现有文献1“SharmaPV.Rapidcomputationofmagneticanomaliesanddemagnetizationeffectscausedbybodiesofarbitraryshape[J].Pure&AppliedGeophysics,1966,64(1):89-109.”提出了积分方程考虑退磁效应下的磁场正演方法，但不包括反演方法。文献2“Lelièvre,P.G.andOldenburg,D.W.Magneticforwardmodellingandinversionforhighsusceptibility.2006,166:76-90.”提出了偏微分方程下的磁总场退磁方法。文献3“Lelievre,Peter.(2021).Forwardmodelingandinversionofgeophysicalmagneticdata.”使用有限体积、有限元和磁偶极子等方法比较了退磁效应下磁总场和磁三分量之间的磁异常特征。文献4“KrahenbuhlRA,LiY.Investigationofmagneticinversionmethodsinhighlymagneticenvironmentsunderstrongself-demagnetizationeffect[J].2017,82(6):1-60.”提出了积分方程的退磁效应下的磁总场磁异常正反演方法。文献5“ZhdanovMS,CaiH,WilsonGA.3DInversionofSQUIDMagneticTensorData[J].2012,1(1).”提出了磁梯度张量数据的三维正反演方法，显著提高了反演的分辨率，但未考虑退磁效应。文献6“ZhdanovMS,WilsonGA,Technoimaging.3Dinversionoffulltensormagneticgradiometry(FTMG)data[C].2011.”利用全张量磁梯度测量来反演澳大利亚塔拉旺地区的一个磁铁矿，反演结果与当地的地质情况非常吻合。目前在全球范围内，地球物理行业内已有偏微分方程和积分方程的退磁效应下的磁总场和磁三分量三维正演研究，但是积分方程下考虑退磁效应的磁梯度张量数据三维正反演和磁总场、磁三分量数据三维反演方法未见提出。
技术实现思路
有鉴于此，针对以上技术问题，本专利技术提供的一种基于积分方程的退磁效应下磁测数据的正反演方法，该方法能够进行磁测数据为磁梯度张量数据的正反演计算以及磁测数据为磁总场、磁三分量数据的反演计算，具体包括以下步骤：S1、根据地下目标区域，确定网格剖分的空间范围，使用正交网格对所述空间范围进行网格剖分，得到剖分后的网格；S2、设定地磁场的磁倾角、磁偏角和磁感应强度；S3、在剖分后的网格中设定观测位置，利用基于积分方程的退磁效应下的磁测数据正演计算方法，获得在退磁效应下观测位置的磁测数据；S4、根据磁测数据，通过基于积分方程的退磁效应下磁测数据的反演方法，结合正则化项进行模型最优化计算，对地下目标进行反演计算，得到反演结果。进一步地，步骤S3中，基于积分方程的退磁效应下的磁测数据正演计算方法，其计算公式为：M＝K(H0+TM)AM＝bA＝I-KTb＝KH0Bij＝μ0TijM,i∈{x,y,z},j∈{x,y,z}其中，M表示考虑退磁效应后的磁化强度矩阵；K表示磁化率矩阵，H0表示磁场强度矢量,Bij异常场的磁梯度张量数据,μ0是真空中的磁导率；I为单位阵。p的观测点位置位于每个正交网格单元的中心。进一步地，步骤S4中，基于积分方程的退磁效应下磁测数据的反演方法的目标函数为：minφ＝||F(κ)-dobs||2+βΔκF(κ)＝GκG＝Tij(B0+μ0TM)其中，φ表示优化目标，F(κ)表示正演过程计算出的磁测数据，磁测数据类型包括磁总场、磁三分量和磁梯度张量数据，dobs表示实际观测磁测数据,β代表规整化因子，Δκ表示模型规整化项，G为灵敏度矩阵，κ表示模型磁化率矩阵，B0为地磁场强度，μ0表示真空磁导率；本专利技术提出的技术方案能使用积分方程方法计算考虑退磁效应下的磁异常梯度张量数据，对高磁化率环境非常适用，同时本专利技术计算出来的磁正反演结果比不考虑退磁效应下计算出来的磁正反演结果误差相对较小，更加符合真实环境。附图说明图1是本专利技术一种基于积分方程的退磁效应下磁测数据的正反演方法的流程图；图2是本专利技术的网格剖分原理示意图；图3是本专利技术实施例提供的一个地下模型网格剖分示意图；图4是本专利技术实施例提供的考虑不退磁效应和考虑退磁效应下的磁总场和磁三分量的磁异常数据示意图；图5是本专利技术实施例提供的不考虑退磁效应和考虑退磁效应下磁梯度张量数据示意图；图6是本专利技术实施例提供的不考虑退磁效应和考虑退磁效应下的反演结果示意图。具体实施方式为使本专利技术的目的、技术方案和优点更加清楚，下面将结合附图对本专利技术实施方式作进一步地描述。请参考图1，一种基于积分方程的退磁效应下磁测数据的正反演方法，包括以下：S1、根据地下目标区域，确定网格剖分的空间范围，使用正交网格对所述空间范围进行网格剖分，得到剖分后的网格。网格剖分原理如图2所示，按照正交网格的形式进行剖分；实施例提供的模型剖分示意图如图3所示，选定X轴、Y轴、Z轴方向为200m×200m×200m范围的区域作为网格剖分的空间范围，用来剖分的网格为10m×10m×10m的六面体形式，示意图中的不规则物体为磁化率为1SI的地下目标；S2、设定地磁场的磁倾角、磁偏角和磁感应强度。常用的参数为磁倾角90度，磁偏角0度，磁感应强度50000nT；S3、在剖分后的网格中设定观测位置，利用基于积分方程的退磁效应下的磁测数据正演计算方法，获得在退磁效应下观测位置的磁测数据；S4、根据磁测数据，通过基于积分方程的退磁效应下磁测数据的反演方法，结合正则化项进行模型最优化计算，对地下目标进行反演计算，得到反演结果。步骤S3中，基于积分方程的退磁效应下的磁测数据正演计算方法，其计算公式为：M＝K(H0+TM)AM＝bA＝I-κTb本文档来自技高网...

【技术保护点】
1.一种基于积分方程的退磁效应下磁测数据的正反演方法，其特征在于：/nS1、根据地下目标区域，确定网格剖分的空间范围，使用正交网格对所述空间范围进行均匀网格剖分，得到剖分后的网格；/nS2、设定地磁场的磁倾角、磁偏角和磁感应强度；/nS3、在剖分后的网格中设定观测位置，利用基于积分方程的退磁效应下的磁测数据正演计算方法，获得在退磁效应下观测位置的磁测数据；/nS4、根据磁测数据，通过基于积分方程的退磁效应下磁测数据的反演方法，结合正则化项进行模型最优化计算，对地下目标进行反演计算，得到反演结果。/n

【技术特征摘要】
1.一种基于积分方程的退磁效应下磁测数据的正反演方法，其特征在于：
S1、根据地下目标区域，确定网格剖分的空间范围，使用正交网格对所述空间范围进行均匀网格剖分，得到剖分后的网格；
S2、设定地磁场的磁倾角、磁偏角和磁感应强度；
S3、在剖分后的网格中设定观测位置，利用基于积分方程的退磁效应下的磁测数据正演计算方法，获得在退磁效应下观测位置的磁测数据；
S4、根据磁测数据，通过基于积分方程的退磁效应下磁测数据的反演方法，结合正则化项进行模型最优化计算，对地下目标进行反演计算，得到反演结果。


2.如权利要求1所述的一种基于积分方程的退磁效应下磁测数据的正反演方法，其特征在于：步骤S3中，基于积分方程的退磁效应下的磁测数据正演计算方法，其计算公式为：
M＝K(H0+TM)
AM＝b
A＝I-KT
b＝K...

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