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

第 32卷 第 3期 地 震 地 质 Vo1．32，No．3 2010年 9月 SEISMOLOGY AND GE0LOGY Sept．，2010 doi：10．3969／j．issn．0253—4967．2010．03．012 具有干预机制的 CSAMT数据一维最优化反演 周俊杰 强建科 汤井田 王海青 1)中南大学信息物理工程学院，长沙 410083 2)安徽省勘查技术院，蚌埠 233000 摘 要 在 CSAMT测深 中，测深曲线往往受到 电偶极源 的影响而 出现畸变 ，这给CSAMT数据 的反 演解释带来很大的困惑 。文 中依据 CSAMT法在波区与MT理论相似 的原理，提出用Bostick半定量结 果作为 CSAMT一维反演 的初始模型，采用传统最小二乘算法 ，结合多种常识性物性变化特点，优化层 参数并控制反演过程 的迭代方 向，使纯粹 的数值迭代转化为按地球物理特征 自动迭代 ，其结果可更好 地满足物探解释需要。算例和实测数据反演表明，优化后的CSAMT一维反演方法，反演精度较高、计 算速度较快，其反演结果的正演响应曲线和实测 曲线拟合得很好 。 关键词 可控源音频大地 电磁法 一维反演 最小二乘法 初始模 型 中图分类号 ：P631 文献标识码 ：A 文章编号 ：0253～4967(2010)03—0453—12 O 引言 可控源音频大地 电磁法 (ControlledSourceAudio—frequencyMagnetotelluric，简称 CSAMT)是 一 种人工源 电磁测深法 ，是为克服大地 电磁测深法 (MT)信号强度弱 的缺点而发展起来 的。 CSAMT法 自20世纪70年代问世 以来 ，在油气、地热资源勘探方面取得 了一系列的应用性成果 (何继善，1990)。由于 CSAMT法沿用了MT中卡尼亚视电阻率的定义，因此严格来讲只能在波 区进行测量。在过渡区和近区，卡尼亚视电阻率由于非平面波效应而出现畸变 ，使资料处理难 度加大 ，要求有成熟的正反演技术来支持。在正演模拟方面 ，Ward等 (1988)研究 了一维水平电 偶极子场的计算 。汤井田等 (1994，2005)通过修改含贝赛尔积分式 中的核 函数 ，改善 了数值计 算的收敛性。殷长春 (1994)提出直接积分法 ，可对双极源场进行正演计算。Anderson(1979)以 及国内龚强等 (2008)对快速汉克尔变换 的滤波系数的合理选取做 了研究。对于 CSAMT法 的 一 维反演而言，传统做法是对观测数据做近场校正 ，或只采用波区数据，再利用 MT理论方法进 行反演 (Routheta1．，1999；王若等，2003)。也可以不做近场校正，直接对观测数据进行反演 ，通 常采取的方法有最小二乘反演 、阻尼最小二乘反演、奥克姆反演 以及各类完全非线性反演方法 (尚通晓等 ，2007；王若等，2007；李帝铨等，2008；何梅兴等 ，2008)，但至今未在实践中普遍开 [收稿 13期】 2010—01—07收稿 ，2010—07—02改回。 [基金项 目] 国家 自然科学基金 、国家高科技发展计划项 目(863—2006AA06Z105)和深部探测技 术与实验研究专项(Sinoprobe一03)共 同资助 。 通讯作者 ：强建科 ，男，陕西岐山人 ，副教授 ，从事地球物理 电磁法正、反演研究，电话 ： E—mail：qiangjianke@163．con。 454 地 震 地 质 展。最小二乘反演具有简单有效的特点，但其对初值依赖较大，以及反演参数太多导致多解性 影响加大等问题仍有待合理解决。 本文拟采用最小二乘法对 CSAMT观测数据直接进行全 区反演。对于正演计算 ，采用双极 源的电磁场模拟可使之更加精确；对于初始模型的选取 ，依据 CSAMT波区与MT观测数据的相 似性 ，利用 Bostick反演获取一个多层模型，并将其层数优化后作为初