基于最小支持的2.5维直流电阻率正则化反演研究

地球物理反问题存在多解性和不稳定性,正则化是得到稳定解的重要手段之一,正则化方法中稳定因子的设计与正则化因子的选择是两项重要研究内容.稳定因子的主要功能是对模型的解空间进行限制,以减少多解性,求得稳定解;为更好的定位地下目标体,研究工作采用了最小支持稳定因子,并给出了统一表达形式.正则化因子在模型误差与数据误差之间起到平衡作用,研究工作在“L曲线法”自动选取正则化因子的基础上,提出一种修正的计算正则化因子算法.以2.5维直流电阻率为例,运用非线性共轭梯度法,进行了试算研究,结果表明了修正的“L曲线法”算法快速、稳定,最小支持稳定因子具有较好的聚焦特性,有利于实现陡变边界反演.     

基于二次场算法的时间域激电2.5D有限元正演

二次场算法是提高时域激电正演计算精度的有效方法之一,能很好处理场源奇异问题。研究工作首先通过引入Cole-Cole模型,推导了在此模型下复电阻率二次场满足的边值问题;并利用有限单元法进行数值计算。然后通过Guptasarma滤波算法,实现了频率域激电到时域激电的转换,改变了传统求解时域激电的方法。为了提高计算效率,首先对边界采取近似处理。经处理的刚度矩阵只与供电频率及波数有关,而与供电点位置无关。其次设计了基于图论理论的矩阵重排与填入元分析算法,实现了矩阵的LDLT分解,从而实现了基于二次场算法的时间域2.5维激电正演的快速算法。利用快速算法,计算并分析了视电阻率、视极化率异常特征;其次通过计算频散率异常特征,验证了极化率与频散率对异常体响应具有相似性。     

基于正则化因子与稳定因子的大地电磁反演研究

地球物理反问题存在多解性和不稳定性,正则化是得到稳定解的重要手段。正则化方法中稳定因子的设计与正则化因子的选择是两项重要研究内容。以大地电磁测深法为例,研究了正则化因子的"L"曲线法自动选择算法;针对不同的稳定因子,研究了稳定因子的特点及稳定因子内部参数的选择原则。研究工作为大地电磁反演提供支持。