基于频点CSAMT一维最小构造反演汤井田张林成肖晓摘要传统.DOC

基于频点CSAMT一维最小构造反演 汤井田 张林成 肖晓 摘要：传统CSAMT线性反演方法往往对初始模型依赖性强和对某些低阻薄层反映不灵敏的弱点，本文在CSAMT最小二乘反演的基础上，在目标函数加入光滑约束条件，构成最小构造反演。另外根据频点个数和地表电阻率来构建初始模型，使反演的初始模型更加灵活，简单，合理，更易应用于实际数据处理中。 关键词：初始模型, CSAMT, 最小构造反演，频点，地表电阻率 One dimension CSAMT minimum structure inversion based on the frequency Tang Jing-Tian, Zhang Lin-Cheng, Xiao Xiao School of Geosciences and Info-Physics, Central South University,Changsha, 410083, china The methods of traditional linear inversion of CSAMT are with more dependence on the initial models, and the thin layer which is with low resistance always can not be reflected by these methods. In this paper, some smoothness constraints conditions were added in the objective function of the least squares inversion , which was called the Minimum structure inversion. Otherwise, we made the number of the frequency equal to the number of the layers, and used the frequency and the surface resistivity to get the thickness of every layer in order to construct the initial model. On these foundations, the initial model of the CSAMT inversion was more flexible, simple and reasonable, and more easily to be applied in the actual data processing. Key words: the initial model, CSAMT, the Minimum structure inversion, frequency, the surface resistivity 1、引言：大多数线性反演中最小二乘反演是最传统的，也是行之有效的方法，一维CSAMT最小二乘反演方法虽然可以精确地模拟电磁场，但它只限于简单的水平层状模型，并且严重依赖于初始模型，对给予的初始模型参数与真实模型相差过大，反演结果会非常差，常常会陷入局部最小，不能达到全局最优，对实际数据资料反演的适用性很小[1] [2]。考虑到实际地下介质的电性分布，在最小二乘反演的目标函数中加入光滑约束，一方面考虑到地质结构的纵向光滑情况，能抑制地电构造的不合理性；另一方面，又可使模型更加灵活，自由，从而构成CSAMT一维最小构造反演，这样就不会引入数据分辨不出的多余构造干扰，具有较好的稳定性，且随着误差的减小，用适当的方法便能得到真实的结构[3] [4]。同时，针对最小构造反演初始模型层参数，提出了基于频点个数分层，勘探深度划分层厚的方法，与传统方法相比，使初始模型更加简单，合理，更趋于真实模型，极大减小了计算量，也更适宜于实际资料的反演。 2、最小构造反演 对于地球物理反演问题，由于实测数据存在误差，数据有限，反演的结果是非唯一的。如果仅仅进行实测数据进行拟合，往往会得到过于复杂的地电模型，或是引入实测数据无法分辨的冗余构造，使得反演过程不稳定，反演结果偏离实际情况[5][6]。因此在反演过程中,为了让模型尽可能的灵活,又抑制地电构造的不合理性,可通过定义模型粗糙度来解决这一问题,也就是在目标函数中加入光滑度约束，即为最小构造反演。 最小构造反演的实质是求取多层地球模型的最光滑解，即在一定的拟合误差标准下模型的粗糙度最小。最小构造反演的优点是不会引入数据分辨不出的多余构造干扰，具有较好的稳定性，随着误差的减小，用适当的方法可得到真实的