双参数混合正则化方法及在电导率反演成像中的应用

通过引入带有二阶正则算子的正则化项,建立了一种双参数混合正则化方法.为确定最佳正则化参数,这里主要应用L-曲线法、偏差原理和广义交叉校验准则的优化组合来确定.首先对理论模型进行了数值模拟,通过与截断奇异值分解法、共轭梯度法及标准Tikhonov正则化法的结果比较,表明该方法不仅精度高,而且对于数据的随机扰动具有稳定性.然后将此方法应用于对电导率反演成像的计算中,得到的电导率反演成像精确细致可靠,也符合实际情况.这表明该方法求解地球物理反问题时,在增强反演的稳定性及减少多解性的同时,还能进一步提高反演计算速度.     

基于Tikhonov正则化的双频电磁波电导率成像反演

本文将Tikhonov正则化方法与active-set算法相结合,利用双频电磁波电导率成像原理,求解其反演成像方程.不仅对现有算法进行了改进,也促进了算法的实际应用.本文研究了在双频电磁波电导率成像方程建立后,如何根据其严重病态性质,选择合适的算法求解矩阵成像方程.针对电导率非负的特性,引入正则化参数,将问题转化为一个非负最小二乘问题,并用active-set算法求解.采用改进后的迭代算法对理论模型进行了数值模拟计算,验证了该方法的有效性.应用到实际电导率成像反演,与常规的LSQR、SP-LSQR、Tikhonov正则化等算法进行比较,取得了满意的结果.     

基于混合差分进化算法的地球物理线性反演

地球物理反问题线性化处理之后, 各种反演算法归结为对病态线性方程组的求解. 为了快速准确地计算出地球物理参数, 本文提出了一种全新的基于LSQR算法的混合差分进化算法(Hybrid Differential Evolution Algorithm, HDE). 该算法利用LSQR算法给出DE算法的初始种群, 提高DE算法的计算速度和稳定性. 在不同噪声水平下, 对四种正则化方法Tikhonov、TSVD、LSQR和HDE的反演结果进行详细比较. 理论模型和实际数据反演的结果都表明: 改进的HDE算法应用于地球物理反问题的求解是成功的: 反演结果与原设定模型具有较高的相关性, 在稳定性和准确性上较常规的反演算法都具有一定的优势; 而且不需要给定正则化参数, 具有更强的实用性.     

高分辨率阵列侧向测井的数学模型及有限元快速正演

对包含井眼、侵入带、围岩和目的层的轴对称地层模型, 推导了无穷远截断边界上的Robin边界条件, 建立了高分辨率阵列侧向测井的等值面边值问题模型. Robin边界条件较Dirichlet边界条件更加精确, 可大大缩小求解区域而不影响计算精度. 考虑到微分方程和边界条件为线性的, 利用叠加原理简化了原微分方程边值问题的计算, 克服了事先屏蔽电极上电流的不确定性. 采用基于地址矩阵的稀疏存贮模式, 大大减小了内存需求, 且地址矩阵物理意义明确, 方便迭代法调用求解有限元方程. 引入预条件共轭梯度(PCG)法求解有限元计算形成的大型线性方程组, 提高了测井响应的计算速度. 利用本文方法定量考察了地层厚度、井径、侵入带等因素对阵列侧向测井响应的影响, 为后续阵列侧向测井反演的研究奠定了基础, 对实际测井工程具有一定的指导意义.