随钻超深前视电阻率测井是一种新型井中电磁测量技术,能够实现复杂储层探测和开发中井眼轨迹优化与地质导向.本文将基于二次场耦合势三维有限体积法建立一套非均质各向异性地层中随钻超深前视多分量电阻率测井高效数值模拟算法.首先,利用Yee氏交错网格、非均质网格上的等效电导率标准化技术以及三维有限体积法实现二次场耦合势Helmholtz方程离散,并采用传输线法建立背景电磁场并矢Green函数库,同时结合二维Newton插值技术确定背景电磁场,提高异常体内散射电流场计算效率.在此基础上,针对仪器移动测量方式的特点,通过并行直接稀疏求解器PARDISO求解离散方程,结合并行计算同时确定不同位置发射源产生的电磁场.最后,数值结果验证了算法的有效性,显示了二次场方法具有更高的计算精度;建库插值方法显著提高了背景场计算效率,拓展了二次场方法的应用价值.
相似文献为了提高二维大地电磁反演对异常体边界的刻画能力, 我们引入曲波变换建立一种新的稀疏正则化反演方法.与传统的在空间域中对模型电阻率参数求解的方式不同, 我们借助曲波变换将二维电阻率模型转换为曲波系数, 并采用L1范数约束以保证系数的稀疏性.曲波变换是一种多尺度分析方法, 其系数分为粗尺度系数和精细尺度系数, 粗尺度的系数代表电阻率模型的整体概貌, 而精细尺度中较大系数代表目标体的边缘细节.此外, 曲波变换的窗函数满足各向异性尺度关系, 并具有多方向性, 因此曲波变换可以近似最佳地提取目标体的边缘特征信息, 这为我们在反演中恢复边界提供有利条件.通过对大地电磁的理论模型合成数据和实测数据反演, 验证了基于曲波变换稀疏正则化反演对异常体边界的刻画能力优于常规的L2范数和L1范数反演方法.
相似文献为更好地处理与解释复杂海底地形条件下测得的海洋可控源电磁数据,本文提出了一种基于非结构网格剖分的频率域海洋可控源电磁数据三维正则化反演方法.该方法首先对海洋地电模型以非结构四面体单元进行离散,然后基于矢量有限元方法获得海洋可控源电磁响应和灵敏度信息,最后采用共轭梯度法求解高斯-牛顿反演方程计算模型修正量.为提高反演的稳定性,通过在反演过程中采用对数转换方法实现反演模型参数的上下限约束.本文分别测试了单测线水平海底地形反演算例和面积性测量的起伏海底地形反演算例.反演结果表明,本文提出的频率域海洋可控源电磁三维反演能够准确地恢复高阻储油层的位置和电阻率信息,且计算效率较高,可用于实测海洋电磁资料的处理与解释.
相似文献常规的三维时间域航空电磁模拟通常采用隐式步长方法进行时间离散,需要几次矩阵分解和上百次右端源项回带,计算效率较低.为了提高正演计算效率,本文提出使用有理Krylov方法求解时间域电场扩散方程.首先使用非结构四面体网格进行空间离散,采用Nédélec矢量基函数近似四面体单元内的电场;然后基于有限元离散给出矩阵指数和矢量乘积表示的电场显式解;最后采用有理Arnoldi算法构造Krylov子空间内的正交基函数并进一步求解矩阵指数与矢量的乘积,直接得到任意时刻的电场解向量,避免步长离散过程.此外,本文还提出一种指数加权偏移参数优化方法,使得有理Arnoldi近似在瞬变衰减晚期具备更高的精度,从而降低Krylov子空间阶数并提高计算效率.通过和层状模型解析解的对比验证了有理Krylov方法的精度.针对三维异常体模型使用全局网格和局部网格剖分并和其他数值方法比较,进一步说明了有理Krylov方法的有效性.
相似文献传统的均匀分布异常体模型,不能准确描述地下介质的不规则性变化,采用随机介质进行替代,将使电导率这一物性参数更加接近于实际的电导率分布.本文结合地质结构中断裂带特征,采用Von Kármán函数建立三维随机介质模型,通过讨论Hurst指数与自相关长度对三维随机电导率模型建模产生的影响,进行参数优化建立所需的随机介质模型.利用一维傅里叶变换建立三维随机变化的条状断裂带异常体,三维傅里叶变换建立三维随机变化的背景围岩,准确地表征了油气藏断裂带特征.基于时域有限差分方法,实现了磁源激励下的三维随机介质航空时域电磁响应数值模拟.采用均匀半空间模型验证了数值模拟的正确性,分析了随机断裂带与均匀断裂带的电磁响应特征,结果表明随机断裂带可以准确描述地下介质的分布特征,而且与断裂带垂直方向的电磁响应特征清晰地描述了断裂带的倾向、走向与位置,为断裂带结构探测提供了理论依据和技术指导,三维随机断裂带模拟方法同样适用于其它三维随机介质的数值模拟.
相似文献