起伏地形条件下CSAMT资料三维反演研究.pdfVIP

起伏地形条件下CSAMT资料的三维反演研究 林昌洪，谭捍东，佟拓 中国地质大学(北京)地球物理与信息技术学院，北京，100083 引言 可控源音频大地电磁法(CSAMT)在勘探石油、天然气、地热、金属矿产、以及水文和环境工 程中发挥了重要的作用fl，2】。然而，在CSAMT的实际工作中，野外工区很难满足水平地表条件，更 多的情况是场源处和测线处高程和地形起伏差异很大，地形起伏是二维甚至三维的。如果反演CSAMT 数据时不考虑地形的影响，会造成很大的误差。正演采用有限差分法，反演采用共轭梯度法，我们开 发了起伏地形条件下CSAMT三维正反演算法。 带地形三维正演 在计算起伏地形条件下三维地电模型的电磁场响应时，我们采用把电磁场的总场分解成一次场 (背景场)和二次场的叠加策略，分别计算出一次场和二次场，再合成总场。一次场通过一维数值模 拟得到。在起伏地形时，选取起伏地形的最高点处的高程为一维数值模拟计算一次场的地表高程。对 于二次场，采用交错采样有限差分方法来求解。在三维网格剖分时，采用“台阶状”网格来逼近地形的 起伏变化。用交错采样剖分网格在研究区域内对剩余电阻率和二次电磁场所满足的麦克斯韦积分方程 进行离散化，可获得关于研究区域内各网格单元采样点处二次磁场满足的正演方程。解正演方程可获 得二次磁场，然后通过二次电场与二次磁场之间的差分关系计算二次电场。进一步，可计算地表视电 阻率和相位响应。 带地形三维反演 我们采用共轭梯度反演方法求解起伏地形条件下CSAMT瓷料的三维反演问题。反滚的数据选用视电阻率和 相位数据。由于背景电阻率是固定值，在反演中选取剩余电阻率为反演参数。具体的反演流程、目标函数雅可比 矩阵的计算等可参考水平地表条件下CSAMT三维共轭梯度反演算法【31。起伏地形时，反演只更新真实地表下方 的模型参数，其它地方保持不变，重新进行正演计算和反演迭代。 227 合成算例 设计的理论模型为：两个大小为200m×200m×120m，电阻率为10Q·m的低阻棱柱体分别埋藏于 起伏地形的山脊和山谷的下方，围岩的电阻率为100Q．m。 分为42×85×45(含10个空气层)。用可控源音频大地电磁法三维共轭梯度反演程序对地表924个测点 处加入随机误差的9个频率(4000-0．IHz)的视电阻率和相位数据进行三维反演。反演结果如图1一 图3所示。 结论 我们实现了起伏地形条件下CSAMT资料三维反演算法。合成数据的反演算例验证了所实现的三 维反演算法的正确性和稳定性。 该研究由国家自然科学基金项目，中央高校基本科研业务费专项资金资助项目和教育 部新世纪优秀人才计划联合资助。 参考文献 ofcontrolledsol／Ice magnetotellurics [1】Routh，P．S．andOldenburg，D．W．，1999，Inversionaudio—frequency datafora 689·1697． horizontallylayeredearth：Geophysics，64(6)，1 relaxationinversionofCSAMT [2】Lu，X．Y．，Unsworth，J．M．，andBooker,J．，1999，Rapid data：Geophysical International，138(2)，381—392． journai 【3】林昌洪，谭捍东，舒晴，佟拓．可控源音频大地电磁三维共轭梯度反演研究．地球物理学报，2011， 55(1I)：3829～3839 雕 。 伸 j毛 珊 一E儿．tg J 287 啪 啪