论文部分内容阅读
在地球物理勘探方法中,常规大地电磁法在深部地壳构造和石油勘测等众多领域有着广泛而深入的应用。而且现阶段三维正反演研究工作在理论和技术都已较为成熟,不过这是在不考虑激电效应而仅仅考虑电磁效应的影响下进行的。然而在地下介质激电效应不容忽视的情况下,这样做从理论上就是有所欠缺,反演结果与实际地质情况也必然存在差异。仅通过更新理论方法进行考虑激电效应的MT三维正反演研究,就可以得到更多地下介质信息,从而使大地电磁方法的应用领域得以扩展。通过借助Cole-Cole模型将复电阻率引入到Maxwell方程中,探究从观测数据中提取Cole-Cole模型中四个参数的可行性,最终编程实现考虑激电效应的大地电磁三维正反演研究。正演选用交错网格有限差分法,通过Cole-Cole模型将Maxwell方程中电导率替换成复电导率,在交错采样网格剖分基础上将Maxwell方程离散化并结合边界条件得到电场值满足的线性方程组。编程实现方程的计算,得到正演电场分量,并由Maxwell方程求出磁场值,进一步计算所需响应值。为了解Cole-Cole模型参数对响应灵敏程度相对大小关系,设计了四组单一参数改变的对比模型。对比各自响应变化情况,也反映了同时反演时各参数反演效果和难易程度。反演采用了数据空间OCCAM法,在拟合实测数据的同时,还包含了对模型水平和垂直方向上的圆滑。由于阻抗对各参数的偏导数矩阵在数值上差异较大,需要对其进行加权处理以改善对响应灵敏程度较低的参数的反演结果。具体加权方式是利用各参数雅克比矩阵二范数的比值来确定各参数拉格朗日因子的大小。进而从目标函数出发推导出灵敏程度相差较大的多参数同时反演时模型参数的迭代公式。此外,针对Cole-Cole模型参数物理特性所决定的取值范围对参数进行相关的限定降低多解性实现其反演。通过理论模型合成数据反演算例证实,开发的正反演算法是正确有效的。