SEISMOLOGY AND GEOLOGY ›› 2010, Vol. 32 ›› Issue (3): 453-464.DOI: 10.3969/j.issn.0253-4967.2010.03.012

• Brief Report • Previous Articles     Next Articles

1-D OPTIMIZATION INVERSION OF CSAMT DATA WITH INTERVENTION MECHANISM

ZHOU Jun-jie1, QIANG Jian-ke1, TANG Jing-tian1, WANG Hai-qing2   

  1. 1. School of Info-physics and Geomatics Engineering, Central South University, Changsha 410083, China;
    2. Geological Exploration Technologies Institute of Anhui Province, Bengbu 233000, China
  • Received:2010-01-06 Revised:2010-07-01 Online:2010-09-30 Published:2011-01-14

具有干预机制的CSAMT数据一维最优化反演

周俊杰1, 强建科1, 汤井田1, 王海青2   

  1. 1. 中南大学信息物理工程学院, 长沙, 410083;
    2. 安徽省勘查技术院, 蚌埠, 233000
  • 通讯作者: 强建科,男,陕西岐山人,副教授,从事地球物理电磁法正、反演研究,电话:13548985729,E-mail:qiangjianke@163.com。
  • 作者简介:周俊杰,男,1985年生,山西忻州人,现为中南大学硕士研究生,从事地球物理电磁法研究,电话:13548985729,E-mail:qiangjianke@163.com.
  • 基金资助:
    国家自然科学基金(40974076);国家高科技发展计划项目(863-2006AA06Z105);深部探测技术与实验研究专项(Sinoprobe-03)共同资助

Abstract: In the application of CSAMT,the sounding curve is often distorted with the impact of electric dipole source,which produces confusions in data inversion and interpretation.Based on the 1-D layered earth conventional least-square inversion for CSAMT data,this paper introduces an optimized technique to get better result by discussing the forward modeling calculation,mechanism of inversion,selection of initial model,and iteration model controlling.Formulas of electromagnetic field components,Ex and Hy,given as integrations which contain Bessel functions of first kind of order 0 and 1,can be obtained using numerical method.During the calculation,the kernel functions of integrations are changed respectively in order to get better convergence properties.Taking the short transmitter-receiver distance into account,the relatively long grounded linear wire source should be assumed as an accumulation of dipoles,thus the theoretical apparent resistivity can be calculated by the field responses of dipole series.The forward calculation of CSAMT needs configuration parameters such as source location,measuring point position,current amplitude and so on,which are different from MT modeling.A primary objective of inversion is to recover a geologically interpretable model that can reproduce the set of observations.Based on the least square inversion theory,the method iteratively updates the vector of model parameters to get the minimum of the vector of relative mean square errors,where the Jacobi matrix can be obtained by difference method.Singular value decomposition method is adopted in equations solving,and we replace small singular value with zero.To get rid of falling into local minimum point during the iteration,it is necessary to control the model parameter excess when the model becomes abnormal.Actually the action will destroy the convergence,but it also increases the chance to find out the global minimum point.Numerical experiments show that the inversion result is better than the conventional one.The selection of initial model has an impact on the inversion results.Based on the theoretical similarity of CSAMT far-field data and MT data,the semi-quantificational Bostick resistivity-depth curve can be used as the base of initial model.Then,by decreasing the layers number according to the resistivity trend with apparent depth,and abandoning the deformed thin layers produced in layer cutting if necessary,the appropriate initial model is created.Some synthetic studies and a case of field data are followed to validate the inversion mechanism,and their resultant response fits the observed data well with acceptable inversion precision and computing time.In conclusion,it is available to adopt semi-quantificational Bostick resistivity-depth curve as the base of initial model for the CSAMT 1-D inversion.The multiplicity of solutions can be suppressed by the number of inversion parameters cut down using the Bostick curve optimization techniques.Appropriate initial models and iteration controlling are introduced to help enhance the capability of algorithm to find the global minimum,which satisfies the geophysical exploration requirements.

Key words: CSAMT, 1-D inversion, least-squares algorithm, initial model

摘要: 在CSAMT测深中,测深曲线往往受到电偶极源的影响而出现畸变,这给CSAMT数据的反演解释带来很大的困惑。文中依据CSAMT法在波区与MT理论相似的原理,提出用Bostick半定量结果作为CSAMT一维反演的初始模型,采用传统最小二乘算法,结合多种常识性物性变化特点,优化层参数并控制反演过程的迭代方向,使纯粹的数值迭代转化为按地球物理特征自动迭代,其结果可更好地满足物探解释需要。算例和实测数据反演表明,优化后的CSAMT一维反演方法,反演精度较高、计算速度较快,其反演结果的正演响应曲线和实测曲线拟合得很好。

关键词: 可控源音频大地电磁法, 一维反演, 最小二乘法, 初始模型

CLC Number: