大地电磁二维正则化反演研究.pdf

《大地电磁二维正则化反演研究.pdf》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、2014年4月铁道：亡程学报Apr2014第4期(总187)JOURNALOFRAILWAYENGINEERINGSOCIETYNO．4(Ser．187)文章编号：1006—2106(2014)04—0019—05大地电磁二维正则化反演研究尹银轩张连伟童孝忠2(1．铁道第三勘察设计院集团有限公司，天津300142；2．中南大学，长沙410083)摘要：研究目的：在大地电磁反演方法的研究中，模型与观测数据一般是非线性关系，观测数据的有限性及含噪声使得反演结果具有非唯一性，地球物理建模理论的近似和先验信息的不足使得各反演方法难

2、以兼顾解决全局收敛和高效收敛问题。通过在目标函数中采用正则化方法，可以使得不适定反演问题具有稳定的反演结果，改善解的稳定性和非唯一性问题。因此，寻求一种合适的二维正则化反演方法显得尤为重要。研究结论：(1)求解适定的大规模反演线性方程组时，为了改善高斯一牛顿法的局部收敛，可以采用阻尼高斯一牛顿迭代算法，通过对理论模型的试验，基于阻尼高斯一牛顿算法的大地电磁二维正则化反演计算结果真实地反映了模型的地电参数，准确地把握住异常体的深度及大小，说明该反演算法是有效的；(2)将模型反演计算的预测数据与实际的观测数据进行对比分析，无论

3、是预测视电阻率数据，还是阻抗相位数据，都与其实际观测数据基本吻合，说明了反演算法的准确性，适用于大地电磁的二维资料解释；(3)利用GCV法可以有效地确定出较优的正则化参数，将不适定问题转化为适定问题进行求解，保持数据目标函数和模型约束目标函数之间的平衡关系，进而使模型方差和分辨率之间达到最佳折衷；(4)该反演方法适于近似水平地形条件下进行大地电磁二维正演模拟，可用于铁路、公路等勘察领域二维大地电磁数据资料解释中。关键词：大地电磁；二维反演；阻尼高斯一牛顿；不适定；正则化中图分类号：P631．3文献标识码：AResearch

4、ontheTwo—dimensionalMagnetotelluricRegularizati0nInversionYINYin—xuan，ZHANGLian—wei，TONGXiao—zhong(1．TheThirdRailwaySurveyandDesignInstituteGroupCorporation，Tianjin300142，China；2．CentralSouthUniversity，Changsha，Hunan410083，China)Abstract：Researchpurposes：Inthestud

5、yofmagnetotelluricinversionmethods，modelsandobservationaldataaregenerallynon—linear，limitedandnoisyobservationdatamakestheinversionresultsnon—unique，approximategeophysicalmodelingtheoryandlackofprioriinformationmakesitdificulttosolvetheproblemofglobalconvergencean

6、dconvergenceeficiencysimultaneouslybyallinversionmethods．Astablesolutionoftheill—posedinverseproblemcanbeobtainedbyutilizingtheregularizationmethodsintheobjectivefunction．Therefore，itisparticularlyimportanttofindasuitabledimensionalregularizationinversionmethod．Re

7、searchconclusions：(1)Solvinglargescalelinearequationofinverseproblem，thedampedGauss—Newtonalgorithmwasadopted，whichcanimprovelocalconvergenceofGauss—Newtonmethod．Throughthesyntheticmodelsimulation，theinversionresultsbasedondampedGauss——Newtonalgorithmtrulyreflecte

8、dthegeo——electricalparametersofthemodelandaccuratelyshowedthedepthandsizeoftheabnormalbody，Sotheinversionalgorithmhasbeenproventobeeffective．(2)Compared