论文部分内容阅读
矩形大定源是瞬变电磁法中应用最广泛的装置之一,其三维正演和反演研究一直是国内外的热点问题。本论文的研究目的是为了以正则化反演技术的时频转换方法为基础,完成时间域信号到频率域响应转换。利用成熟的频率域数据处理技术,实现矩形大定源瞬变电磁数据的三维频率域非线性共轭梯度反演。为实现时间域到频率域变换,假设下阶跃电流波形,从而时间-频率响应的傅里叶变换关系,可以表示为以频率域电磁场响应为未知变量,时间域响应为右手项的线性方程组。由于该矩阵方程的系数矩阵的奇异性,采用正则化反演技术进行求解,其中最优正则化因子采用L-曲线法选取。基于正则化反演技术的时间域到频率域变换的实现为时间域数据的频率域解释奠定基础。在矩形大定源数据频率三维反演过程中的三维正演模拟计算采用交错网格有限差分法。正演模拟以二次场满足的Helmholtz方程为基础。为了计算矩形大定源在背景介质中的格林函数,将矩形发射回线剖分成等面积的虚拟小方回线,任意一点的电磁场响应等于每个虚拟方形回线在该测点处产生的电磁响应的叠加。采用虚界面法时发射源可位于任意层位,计算点可扩展到发射回线边界外的三维空间任意位置。Hankel积分采用直接数值积分技术计算。基于二次场满足的Helmholtz方程,克服了数值模拟过程中源描述的困难以及源奇异性。三维线性方程组利用MUMPS因式分解软件包进行求解,然后利用插值技术获得全空间任意一点的矩形大定源频率域电磁响应。如果结合余弦变换,可实现矩形大定源瞬变电磁全空间三维正演模拟。通过与解析解和前人的研究成果进行对比,验证了本论文提出的结合频率域三维正演和余弦变换技术计算水平地形条件下三维瞬变电磁响应的模拟方案是可行的。设计均匀半空间嵌入低阻异常体模型,模拟地面、地井和地空全空间矩形大定源观测系统。通过估算早期时间道地下感应涡流场的扩散范围,设计矩形大回线外的观测范围。试算结果表明:(1)地面矩形回线源在地下激发的涡流场以“烟圈”的形式不断向下扩散,随着时间衰减,“烟圈”范围扩大,强度减弱;低阻异常体对电磁场有“吸附”性质,低阻体产生的电磁场随着时间延迟衰减较慢;(2)井中瞬变响应受观测位置影响较大,特别是异常体边缘和异常体内部,电压剖面曲线在异常体位置附近变化较大,且中期时间道响应可能会出现变号现象;(3)对于浅部工程地球物理勘查,在瞬变电磁的数据解释和处理过程中,要充分考虑发射源几何形状的影响。频率域三维反演借助非线性共轭梯度反演方法(NLCG)完成。非线性反演过程中不需要显式地计算并且存储灵敏度矩阵,只需计算灵敏度矩阵与向量的乘积,只需要一次正演和一次伴随正演。分别设计简单和复杂模型进行反演试算。对合成数据的反演结果表明:(1)非线性共轭梯度法具有全局收敛性质,可以较准确的反演出地下电性异常体的位置、大小和电阻率等地球物理信息;(2)在矩形大定源多源观测系统中,不同位置的发射源对异常体具有不同的电磁耦合强度,所观测到的数据具有更丰富的地下异常体信息,三维反演效果明显;(3)对比三维反演结果和理论模型,充分证明了矩形大定源非线性共轭梯度三维频率域反演的稳定性。为了将上述的频率域三维反演技术应用于野外实测的矩形大定源瞬变电磁数据的解释,需要开展如下的预处理工作。(1)对数据进行滤波和关断时间矫正,得到理想阶跃发射波形条件下的感应电动势数据;(2)利用感应电动势和时间域磁场的积分关系,将感应电动势数据转换为时频转换所需的时间域磁场;(3)利用正则化反演技术完成时频变换,获得三维反演所需的频率域电磁响应的实部和虚部数据;(4)采用NLCG对复磁场数据进行三维反演。将上述技术应用到中国山西省沁水盆地东部边缘的泊里井田工区采集的矩形大定源瞬变电磁数据的三维反演中,反演获得的电阻率空间分布与测区实际的地质模型基本一致,证明了本论文采用非线性共轭梯度反演技术对矩形大回线源瞬变电磁数据进行三维频率域解释是可靠的。基于正则化反演技术实现了瞬变电磁数据到频率域响应转换和频率域非线性共轭梯度法的融合,为瞬变电磁三维反演提供新的思路,为野外瞬变电磁数据的三维解释开辟了新的途径。