汉克尔自适应算法在电磁计算中的应用

对于地球物理电磁计算中常出现的贝塞尔函数积分,对比几种高精度的快速汉克尔长滤波系数的精度,选择出一种对于一般情况都表现较好的汉克尔滤波系数,在其基础上,进行一个自适应的截断,来计算相关核函数的积分问题。这样既能保证计算的精度,又能节省计算时间。通过计算实例表明,该算法在保证精度的前提下,相比汉克尔滤波直接计算方法,能够节省大量的时间,尤其对于积分核函数相对复杂的情况,节省时间量尤其显著。在地球物理数据计算量越来越繁琐的今天,该算法对节省计算机资源,具有可行性和实用性。     

电导率连续变化的线源FCSEM有限元正演模拟

为了更好地模拟地下介质连续变化及开展连续介质的反演,对二维电导率分块线性变化的线源频率域可控源电磁法进行了有限元正演模拟,在剖分单元内同时对电场及电导率参数线性插值,使电导率参数在剖分单元之间保持连续变化。首先,提出有限元正演模拟的边值问题及变分问题,并详细论述了有限元的剖分、插值、单元分析及总体合成的各个步骤;其次,采用稀疏存储及基于不完全LU分解的BICGSTAB算法求解复系数方程组,节省了内存并提高了计算速度;然后,对一个均匀半空间模型进行模拟,计算结果表明,低频及高频的有限元数值解都与解析解吻合,证明了算法的正确性;最后,对水平层状模型及垂直断层模型进行正演计算,视电阻率及相位的等值线图均较好地反映出了异常体,说明文中算法能够对电导率连续变化的线源可控源电磁法进行有效地模拟。     

基于MLS形函数处理边界的电阻率法有限单元正演模拟

电阻率法有限单元正演模拟中,采用第三类边界条件时为保证精度仍要求较大范围的计算域.无单元法为地球物理领域的新兴正演模拟方法,其计算效率低,但其中采用的移动最小二乘(MLS)形函数相比于有限单元法形函数具有良好的连续性,模拟精度高.本文将MLS形函数应用于电阻率法有限单元2.5维正演的第三类边界条件处理,提出电阻率法有限单元-移动最小二乘(FEM-MLS)耦合正演方法.通过不同正演方法的模型算例模拟结果对比,验证了本文算法的有效性,并讨论了各个参数选择对模拟结果的影响.本文数值模拟结果表明采用第三类边界条件时,在同等计算精度前提下,FEM-MLS耦合法相比于有限单元法可进一步缩小计算域并提高了计算效率,相比于采用较大计算域满足边界条件的有限单元法计算效率提高了约一倍,相比于采用相同小范围计算域的有限单元法平均精度提高了约一倍.     

水平层状介质电磁法数据时频联合反演

通过一维层状介质电磁法数据的时域和频域正演结果,对不同电阻率大小的电阻层进行了反演计算,说明了在这种情况下电磁法数据时频联合反演的合理性.正演计算了一维层状模型下水平电偶极频率域水平电场响应及利用G-S变换计算时间域垂直磁场导数响应;利用最光滑模型约束反演方法建立了时频联合反演迭代公式.通过理论模型反演试算表明,该反演效果较好.     

一维时域磁场与频域电场联合反演方法研究

多方法、多参数的联合是地球物理反演的一个重要的发展方向.本文中,在水平电偶极发射下,频率域利用快速汉克尔变换法计算一维水平层状大地的水平电场响应,时间域利用G-S变换计算垂直磁场对时间的导数值.设计了四个简单模型进行对比,正演计算表明时域磁场对低阻体反映极为灵敏,对高阻体反映不灵敏;而频域电场对低阻体及高阻体反映都比较...     

基于非结构三角网格的海洋CSEM和MT二维联合反演研究

采用非结构三角形网格结合快速Occam算法对海洋可控源电磁(CSEM)与大地电磁(MT)数据进行二维联合反演试验研究。非结构三角形网格能够准确地模拟起伏地形和复杂地质构造,对反演目标区域采用精细网格剖分,其他区域采用粗网格剖分,在满足精度的前提下减少了不必要的计算量。将CSEM与MT数据加入到同一反演数据集中,通过对联合反演数据权重公式进行推导,构建CSEM及MT数据相关权重因子,控制不同数据的拟合权重来实现联合反演。最后对不同模型进行了反演计算,结果表明联合反演比单一反演对海底构造与异常体的还原度更高,验证了联合反演算法的可靠性。