论文部分内容阅读
可控源音频大地电磁法(CSAMT)是基于大地电磁法(MT)和音频大地电磁法(AMT)的一种电磁法勘探方法,历经几十年的发展,方法从最初的定性分析到较精准的定量解释,而正反演算法的开发对方法的发展有着重大的作用,在实际生产中,不仅电阻率参数对观测数据有较强的影响,局部地区磁导率参数对观测数据的影响也不可忽略,为了保证定量解释的结果更加精确可靠,需要开发相应的正反演算法。本文主要的研究内容是在可控源音频大地电磁法中引入磁导率参数,正演过程中电磁场满足的方程不忽略磁导率参数的影响,选用三维交错网格有限差分法作为正演算子,采用总场分离方法计算二次场,选取两种接收点场值的插值方法,对比两种插值方法的优缺点,选择更适用于实际情况的体积加权计算方法,分别设计一维、二维、三维模型对正演算法精度进行检验,设计模型探究磁导率参数的响应规律及对观测数据的影响,由此得到在部分区域磁导率参数的影响不可忽略的结论。在实现电阻率和磁导率双参数模型正演算法的基础上,进一步实现电阻率和磁导率双参数模型的反演算法。采用有限内存拟牛顿法(LBFGS)作为反演研究手段,反演研究重点为目标函数梯度的计算和不同物性参数灵敏度差异的处理。目标函数梯度采用多元函数求导的法则进行计算,本文先实现了两种物性参数的同步反演,并提出双参数模型正演响应相互抵消时反演存在的问题,进而引入基于交叉梯度约束的两种物性参数同步反演,来实现两种物性参数反演的算法;对于不同参数灵敏度差异的问题,采用基于合作式反演的策略进行处理,即先完成两种物性参数的同步反演,将灵敏度较高的物性参数的同步反演结果作为初始模型,进行引入交叉梯度约束的同步反演,得到最终的双参数反演结果,并通过设计不同模型来检验反演算法的可靠性。两种物性参数理论试算证明,忽略磁导率参数进行单电阻率参数反演,反演结果中会出现额外的虚假异常,由此说明进行两种物性参数反演的必要性。