频率域线源CSAET二维反演数值计算.PDF

第 25 卷第 11 期 中国有色金属学报 2015 年 11 月 Volume 25 Number 11 The Chinese Journal of Nonferrous Metals November 2015 文章编号：1004-0609(2015)-11-3196-06 频率域线源CSAET 二维反演数值计算 1, 2 1, 2 1, 2 曹创华 ，童孝忠 ，柳建新 (1. 中南大学 地球科学与信息物理学院，长沙 410083 ； 2. 中南大学 有色资源与地质灾害探查湖南省重点实验室，长沙 410083) 摘 要：由于实测数据包含噪声和其固有的不完备性，决定了频率域电磁测深数据在数学计算上是一个非唯一解 的、不稳定的问题，这就意味着不同的地电模型可以拟合在同一剖面上的视参数数据。面对这种问题可以通过给 定加权的正则化方法使得病态逆问题变为稳定。本研究实现了频率域线源可控源音频大地电流法(CSAET 法)剖面 的二维反演计算。为了使反演过程稳定，采用了模型参数的平滑度约束的方法，并使用自适应正则化参数技术客 观上调整了反演参数。并且在计算过程中利用了伴随方程近似计算方法，减少了计算灵敏度的耗时。数值计算表 明，利用合成的数据可以得到精确的反演结果。 关键词：CSAET 法；线源；二维反演；平滑限制 中图分类号：P631.3 文献标志码：A Inversion numerical calculation of 2D CSAET generated by line current source CAO Chuang-hua1, 2, TONG Xiao-zhong1, 2, LIU Jian-xin1, 2 (1. School of Geosciences and Info-Physics, Central South University, Changsha 410083, China; 2. Key Laboratory of Non-ferrous Resources and Geological Detection, Ministry of Hunan province, Changsha 410083, China) Abstract: The frequency domain electromagnetic sounding calculation is a mathematically non-unique problem with instability because the observed data contain noise and its inherent incompleteness. Thus different models can fit the data from the same cross-section well. This problem can be solved by using regularization method. This study implements the two-dimensional inversion of frequency domain line source CSAET. The method of model parameters smoothness constrain