基于高斯牛顿法的二维直流电阻率法的快速反演.pdfVIP

第 25卷第 6期 中国有色金属学报 2015年 6月 Volume25Number6 TheChineseJournalofNOHferrousM etals June2O15 文章编号：1004-0609(2015)06—1662—10 基于高斯牛顿法的二维直流电阻率法的 快速反演 赵东东，张钱江，戴世坤，陈龙伟，李 昆 (中南大学 地球科学与信息物理学院，三维电磁反演与成像实验室，长沙 410083) 摘 要：从二维线源问题出发，对二维直流电阻率法高效、高精度正反演方法进行研究。在正演数值模拟中，引 入直接解法求解器求解线性方程组，既保证了起伏地形条件下有限元法正演数值模拟的计算精度和计算效率，又 为反演算法中 “拟正演”快速回代求解提供了条件。结合高效、高精度的正演算法，采用高斯牛顿法对电阻率进 行反演成像。在弱非均匀介质前提下，基于近似海森矩阵主对角线元素严格 占优的特点，采用舍弃海森矩阵非对 角线元素的策略，提高整个反演计算的效率。最后，利用合成数据对反演算法的有效性进行检验。结果表明：给 出的反演算法稳定、快速，结合偶极一偶极装置和三极装置，能有效地反演出异常体的形状、大小和位置。 关键词：直流 电阻率法；高斯牛顿法；起伏地形；反演 中图分类号：P631 文献标志码：A Fastinversionfortwo-dimensionaldirectcurrentresistivitymethod basedonGauss．Newtonmethod ZHAODong—dong，ZHANGQian-jiang，DAIShi—kun，CHENLong—wei，LIKun (TheLaboratoryofThree-dimensionalElectromagneticImaging，SchoolofGeosciencesandInfo-Physics， CentralSouthUniversity,Changsha410083，China) Abstract：Fastandhigh-precisioninversionmethodfortwo—dimensionallinesourceproblemwasstudied．Intheforward numericalsimulation，linearequationssolverwasappliedofrdirectsolution，whichnotonlyimprovedtheprecisionand htespeedofnumericalsimulationoffiniteelementmethodinthecaseofruggedtopography,butalsoprovidedconditions forthe q‘‘uasifowr ard’fastback substitution solution in the inversion algorithm．Combined with high efficient simulationmethod，Gauss-Newtonmethodwasadoptedforinversionofresistiviyt．Inthecaseoflow inhomogeneiyt,the maindiagonalelementsoftheapproximatedHessenmatrixpossessedprioriyt thanohters．Basedonthis，non·diagonal elementsweredeletedwhenGauss—Newtoniterativeequationsweresolved．Thewholeprocessofinversionwasmade moreefficientbythisscheme．F