起伏地表频域/时域航空电磁系统三维正演模拟研究

由于航空电磁系统具有工作频率低、时间延迟短等特点,地形对航空电磁响应有很大影响,忽略地形影响会给航空电磁数据解释造成很大误差.本文将基于非结构化网格的矢量有限元法应用于模拟起伏地表条件下频域/时域(FD/TD)三维航空电磁系统响应.该方法由于采用非结构网格,与传统的结构化网格电磁正演算法相比,能更好地拟合地形和地下不规则异常体,提高对不规则地形和地下介质航空电磁响应的计算精度.通过将计算结果与半空间模型的半解析解及已发表的结果进行对比,检验了本文算法的精度.通过对典型山峰和山谷地形航空电磁响应分析对比,总结了地形对航空电磁响应的影响特征.研究结果对航空电磁地形效应的识别和校正具有指导意义.     

2.5维起伏地表条件下频率域地空电磁正演模拟

频率域地空电磁探测方法是指在地面布设人工场源,在空中测量电磁场的一种高效的地球物理勘探技术.该方法具有大范围、高通过性、快速测量的优势,尤其适合崎岖山地、沙漠、沼泽、海陆交互带等复杂地貌区域的资源勘查.但是这些地区的地形起伏通常较大,因此分析地形对地空电磁响应的影响具有重要意义.本文利用有限元法对频率域地空电磁响应进行了正演计算,分析了起伏地表条件下的频率域地空电磁响应特征.首先利用傅里叶变换将2.5维问题转化成二维问题,利用伽辽金加权余量法推导了相应的离散有限元方程组.采用任意四边形单元对区域进行不均匀网格剖分,源和异常体附近网格加密处理,保证计算精度,远离目标区域网格逐渐稀疏,模拟无穷远边界,降低对计算资源的要求.在单元内进行插值,将有限元方程组变换为线性方程组,采用总场算法,利用具有一定面积的伪δ函数表达源电流分布,源项近似为分布在以电偶极源为中心的25个节点上.通过求解线性方程组得到波数域电磁响应,再对波数域电磁场响应进行反傅里叶变换从而获得空间域2.5维频率域电磁场值.通过对比2.5维正演结果与均匀半空间解析解,验证了本文算法的精度,同时本文还对地空电磁场与地面电磁场的响应特...     

三维时间域航空电磁有理Krylov正演研究

常规的三维时间域航空电磁模拟通常采用隐式步长方法进行时间离散,需要几次矩阵分解和上百次右端源项回带,计算效率较低.为了提高正演计算效率,本文提出使用有理Krylov方法求解时间域电场扩散方程.首先使用非结构四面体网格进行空间离散,采用Nédélec矢量基函数近似四面体单元内的电场;然后基于有限元离散给出矩阵指数和矢量乘积表示的电场显式解;最后采用有理Arnoldi算法构造Krylov子空间内的正交基函数并进一步求解矩阵指数与矢量的乘积,直接得到任意时刻的电场解向量,避免步长离散过程.此外,本文还提出一种指数加权偏移参数优化方法,使得有理Arnoldi近似在瞬变衰减晚期具备更高的精度,从而降低Krylov子空间阶数并提高计算效率.通过和层状模型解析解的对比验证了有理Krylov方法的精度.针对三维异常体模型使用全局网格和局部网格剖分并和其他数值方法比较,进一步说明了有理Krylov方法的有效性.

     

起伏地形下线源频率域电磁场二维正演模拟

为了研究地形对线源频率域电磁测深测量结果的影响,本文介绍了起伏地形下的线源频率域电磁场的二维正演模拟方法.在正演模拟过程中采用有限单元法,首先为了模拟起伏地形,用三角形单元对研究区域进行网格剖分;然后再对每个单元进行线性插值、积分、总体合成;最后通过求解线性方程组得出地表的电磁响应.文章最后给出了山谷与山脊模型的数值模...     

时间域航空电磁法一维正演研究

基于电磁勘探理论，导出了层状大地条件下时间域航空电磁法（偶极-偶极装置）的正演计算公式和算法，编制了相应的计算机程序，对若干典型地电断面作了正演计算. 计算结果说明时间域航空电磁法的探测能力和探测条件，进而为时间域直升机航空电磁系统的设计方案提供了依据.     

时间域航空电磁系统瞬变全时响应正演模拟

近年来很多专家学者致力于时间域航空系统正反演的研究.本文针对一维均匀层状介质和三维模型进行正演.不仅计算垂直方向上的电磁响应,还计算了同线方向上的电磁响应,为航空电磁多分量观测提供理论依据.通过比较航空电磁系统的脉冲响应和阶跃响应特征,发现脉冲响应在早期时间存在奇异性,造成计算不稳定.然而,阶跃响应在早期时间没有奇异性,因而利用系统的阶跃响应可得到一种计算时间域航空电磁系统全时响应的稳定算法.该算法具有较高的精度,并很好地保持了磁场强度B和磁感应dB/dt关系的一致性.该算法推广到三维地质体的时间域正演模拟亦取得很好的效果.     

三维起伏地形条件下航空瞬变电磁响应特征研究

航空瞬变电磁法以其速度快、成本低、通行性好等的优势能够有效的应用于地质地形条件复杂的地区.目前对于航空瞬变电磁法的研究主要基于平坦地形的理想情况,对于地形效应的研究相对较少,然而实际应用中地形不可避免,若忽略地形影响将对资料解释造成较大的误差,从而制约航空电磁方法的进一步发展.本文基于交错网格的时域有限差分方法对三维起伏地形条件下航空瞬变电磁进行正演模拟,在保证算法准确性的前提下给出大量模型算例.以经典地形模型为例,利用所给方法计算三维正演响应,结果显示起伏对于航空瞬变电磁数据有着显著的影响且影响主要集中在早期.而后,以实际地质资料为基础,构建起伏地形条件下包含多个异常体的三维复杂模型,计算了复杂模型的航空瞬变电磁响应,并给出三维全域视电阻率曲线,从而对地形效应的影响有了更加直观的认知.最后,通过大量模型讨论了地形的尺寸参数、电性参数、飞行轨迹与飞行高度等因素变化对于航空瞬变电磁数据的影响情况,并得出有价值的结论.

     

基于迭代有限元法的大地电磁二维正演

合理的模型剖分方案是影响大地电磁正演效率的一个重要因素,经典有限元算法为满足控制方程的无穷远边界条件,会在较大的计算空间内进行网格剖分,虽在边界区可以按等比例进行扩展,但依然会形成较高阶的线性方程组,在求解时计算效率较低.针对上述问题,本文开展了基于迭代有限元算法的大地电磁二维正演研究,首先阐述了迭代有限元算法的基本思想及实现过程,建立了基于迭代有限元算法的大地电磁正演模型;其次,结合理论模型的试算,通过与解析解及经典有限元算法的计算结果进行对比分析,验证了迭代有限元算法的准确性及鲁棒性;最后,分析了算法中不同参数对正演精度的影响.结果表明基于迭代有限元算法的大地电磁正演具有计算时间短,占用内存低,能更好的满足远边界条件的优点,可有效提高大地电磁的正演效率,也为后续的反演提供新思路．     

基于瞬时电流脉冲的三维时间域航空电磁全波形正演模拟

传统时间域航空电磁全波形正演模拟主要采用间接法（褶积算法）和直接法（时域有限差分方法等）,然而褶积算法需要获得精确的电流二阶导数,这给发射电流数据采集工作带来极大挑战;时域有限差分方法受到网格和时间步长的严格限制,缺乏灵活性.为解决这些问题,本文采用时域有限元方法,通过直接改变每个时间道上的瞬时电流强度模拟任意发射波形的电磁响应.由于无需计算电流二阶导数,大大提高了正演结果的精度.利用基于非结构四面体网格的矢量有限元方法和后推欧拉技术对时间域电场扩散方程进行空间和时间离散,实现三维航空电磁时间域全波形的直接正演模拟.由此不仅可以模拟复杂的地电结构,而且基于后推欧拉法的无条件稳定性,可以更加灵活地选取时间步长,提高计算效率.通过与1D数值模拟结果进行对比验证了该方法的准确性.本文对三维柱状体模型上HELITEM MULTIPULSE和VTEM系统实际发射波形电磁响应进行模拟,并与褶积算法的结果进行比较,验证了本文算法模拟实际发射波形电磁响应的优越性.对复杂三维地质体模型上不同发射波形电磁响应进行模拟,验证了时间域有限元算法可有效处理复杂地下地质结构.

     

时间域航空电磁各向异性大地三维自适应有限元正演研究

各向异性介质模型电性结构复杂,如何进行合理的网格剖分成为获得高精度正演结果的关键,为此本文开展时间域航空电磁各向异性大地三维自适应有限元正演算法研究.通过结合非结构时间域有限元算法和自适应网格优化技术,实现各向异性介质条件下三维时间域航空电磁自适应正演.考虑到时间域航空电磁响应随时间的衰减特性,为了综合评价不同时刻的后验误差,本文将时间作为加权因子,调整各个时刻后验误差的相对权重,进而实现对浅部和深部网格的同步优化.通过与一维解析结果进行对比验证了本文算法的可靠性.数值实验结果显示电导率各向异性对自适应网格影响严重,其最大主轴电导率的数值及其分布特征直接决定了网格加密效果.此外,各向异性对时间域航空电磁三分量响应的分布形态和异常幅值也会产生严重影响,利用全域视电阻率极性图,可以很好地识别各向异性主轴方向.

     

基于非结构化网格的2.5-D直流电阻率自适应有限元数值模拟

2.5-D直流电阻率有限元数值模拟中,模型的剖分及加密主要通过手动实现.另外,采用的单元类型比较规则如矩形单元等,不易实现复杂模型的模拟.为解决上述问题,文中提出了一种自适应有限元算法.算法中采用稳健的后验误差估计来自动预测下一次网格的单元尺寸,直到设定的迭代条件满足为止.另外,采用非结构化三角形单元实现了任意复杂模型的灵活剖分.基于此,利用垂直接触面模型分析和对比了不同自适应策略的效率.通过对比发现,点源附近的单元得到了加密以消除源的奇异性.另外,对于任意一种策略,有限元结果均能最终收敛到精确解.最后,模拟了两个模型：2-D单个异常体模型和2-D地形模型.     

频率域航空电磁数据变维数贝叶斯反演研究

传统的梯度反演方法已经广泛应用于频率域航空电磁数据处理中,然而此类方法受初始模型影响较大,且容易陷入局部极小.为解决这一问题,本文采用改进的变维数贝叶斯反演方法实现航空电磁数据反演.该方法根据建议分布对反演模型进行随机采样,并依据接受概率筛选合理的候选模型,最终获得反演模型的概率分布和不确定度信息.为解决贝叶斯反演方法对深部低阻层反演效果不佳的问题,本文通过引入合理加权系数,调整对反演模型约束强度,在很大程度上改善了反演效果.通过对模型统计方法进行改进,在遵循原有模型采样方法和接受标准的基础上,将满足数据拟合要求的模型纳入统计范围,削弱不合理模型对统计结果的干扰.本文最后通过对含有高斯噪声的理论数据和实测数据进行反演,并与Occam反演结果进行对比,验证了该方法的有效性.     

时间域航空电磁法激电信息提取研究

时间域航空电磁晚期道数据中常出现符号反转现象,该现象被证明是由激电效应引起.鉴于利用传统实电阻率模型的成像和反演技术处理该类数据无法得到正确的反演结果,本文提出基于Occam理论的反演方法,对层状模型电阻率和充电率参数进行同时反演.本文首先通过正演模拟极化介质中电流密度的分布情况,对时间域航空电磁信号出现多次符号反转现象做出合理解释;然后通过链式法则和波形褶积技术计算Cole-Cole模型参数灵敏度,进而探讨不同Cole-Cole模型参数灵敏度差异及其对反演结果的影响.为了提高多参数反演的稳定性和准确性,在反演过程中对参数采取上下限约束,最后对B_z、dB_z/dt响应单独反演结果和联合反演结果进行对比和分析,发现联合反演可以有效地降低反演结果的多解性.     

基于非结构有限元的时间域海洋电磁三维反演

海洋电磁法是一种有效的钻前储层评价手段,可识别出海底构造是否储油气,从而减少干井率降低勘探成本.近年来,频率域三维海洋电磁反演得到了快速发展,但受到空气波的影响,其在浅海环境中无法取得很好的效果.为解决这一问题,本文研究时间域海洋电磁数据三维正反演.正演模拟和计算中,我们选择基于非结构网格的矢量有限元方法.该算法中空间离散采用非结构四面体网格,可很好地拟合复杂海底地形条件和地下复杂结构;而对于时间离散,我们采用无条件稳定的后推欧拉方法,以确保任意时间步长数值计算的稳定性.反演计算中,灵敏度信息采用伴随正演隐式进行计算;同时,依据时间域反演方法的特点采用L-BFGS方法计算模型修正量.通过利用合成数据反演结果证明了本文提出的三维时间域反演方法可用于复杂海底环境,特别是在浅海环境下的有效性.

     

电导率分块均匀的瞬变电磁25维有限元数值模拟

在吸收并修正部分前人工作的基础上，给出了电导率分块均匀的瞬变电磁25维有限单元法正演算法. 在该算法中，（1）采用在矩形网格中进一步三角细化的剖分方式，同时考虑到矩形单元中心节点的性质，求解线性方程组前用高斯消元法消除矩形网格中心节点上的待求未知量，这样一方面较准确地模拟出任意二维复杂地电断面，另一方面又大大地节省了计算量；（2） 采用计算二次场的算法，事先仅需获得电场和磁场沿地电构造走向方向两个分量的一次场，故而不但计算精度较高，而且不明显增加计算量. 最后，对几例典型模型进行正演，并与解析解对比，验证了本文方法的有效性.