共查询到20条相似文献，搜索用时 62 毫秒
1.
COMSOL是一个多物理场仿真软件,笔者尝试将其应用于大地电磁法的数值模拟,研究起伏地形对大地电磁响应的影响及解决方案。本文由麦克斯韦方程组出发,施加第一类边界条件,将地形起伏数字化为参数化曲面,导入COMSOL进行计算。在对COMSOL的精度进行了验证的基础上,分析不同频点所需内存以及自由度,然后选择合适的频点进行研究。将三维起伏地形的模拟结果与无地形起伏时的结果进行对比,总结出山谷、山峰地形对三维大地电磁响应的影响,发现在XY模式下,地形对视电阻率的影响大于相位,在视电阻率等值线图中能够清晰地看出起伏地形的轮廓,而在YX模式下,地形对相位的影响大于对视电阻率的影响。最后,利用“比值法”进行地形校正,校正前由于起伏地形影响,在异常体附近会出现部分视电阻率数值跳跃现象,通过校正有效抑制了起伏地形对异常体附近视电阻率的影响,校正后的视电阻率曲线整体上移,并且异常体周围的视电阻率值也趋于背景场值,获得了较好的效果。 相似文献
2.
起伏地形对可控源音频大地电磁(CSAMT)响应具有强烈的影响,因此在CSAMT数据处理解释时需要考虑地形。同时,实际的地下地质情况和地表的地形情况通常比较复杂,地质结构和地形大部分情况下都是三维的。在水平地表三维有限差分CSAMT数值模拟算法的基础上,推导了利用地下交错网格采样点处的总磁场计算起伏地形下空气-地下介质分界面处的总电场和总磁场的表达式,从而实现了起伏地形下三维CSAMT数值模拟算法。在算法实现过程中,采用伪δ函数代替麦克斯韦方程中的场源项和直接计算总场的策略,避免了原有的将总场分离成背景场和二次场的策略在复杂地质条件下难以选择合适背景电阻率的问题。为了直接模拟总场,起伏地形下三维CSAMT数值模拟算法给出了新的三维正演方程的边界条件。将模拟水平地表三维异常体和三维山峰地形两个理论模型得到的响应结果与前人算法的计算结果进行对比,验证了所实现算法的正确性和有效性。 相似文献
3.
采用易于模拟复杂地形起伏和倾斜界面的非结构三角单元剖分网格,并利用对偶加权后验误差估计指导网格自动细化过程,实现了2.5维直流电阻率法自适应有限元数值模拟。在实例模型分析中,分别计算了层状模型和垂直岩脉模型的直流电阻率响应,并与其解析解进行了比较。对比结果表明,该算法所得数值解精度很高,解的相对误差小于0.5%。最后,计算了起伏地形2.5维地电模型视电阻率异常,并利用比较法进行了地形改正。地形改正结果与水平地形时的结果对比表明,比较法可以较好地消除地形影响,突出局部地质体异常。 相似文献
4.
角域是描绘起伏地形的基本单元,通过角域叠加可以获得复杂纯地形产生的视电阻率异常,进而可以进行地形校正。将简单角域视电阻率畸变进行组合叠加得到复杂角域地形的视电阻率畸变,分析了组合叠加误差,并与ANSYS软件数值模拟结果进行了对比。结果表明,在地形起伏不大的情况下,可以通过组合叠加获得复杂角域地形的视电阻率异常。最后给出了利用组合叠加法获得纯地形影响后应用比较法进行地形改正的模型算例;与水平地形时的结果对比表明,比较法可以削弱地形影响,突出局部地质体的异常。 相似文献
5.
2.5维起伏地表条件下坐标变换法直流电场数值模拟 总被引：3，自引：1，他引：2
研究起伏地表对视电阻率分布的影响是进行地形校正的基础。由于很难处理不规则边界,计算简洁且效率高的有限差分方法很少用于解决起伏地表问题。为了解决该问题,引入曲化平思想,实现了一种基于坐标变换法的起伏地表条件下的直流电场数值模拟方法。方法从传统的2.5维基本方程及边界条件出发,通过坐标变换将起伏地表问题转化为水平地表问题,并利用有限差分法求解水平地表问题,最后再通过坐标映射得到起伏地表条件下的电位和视电阻率分布。精度分析及计算实例表明：本方法实现简洁,计算误差主要集中在震源附近,整个计算区域内的平均相对误差为1.39%,计算结果满足起伏地表条件下的电位和视电阻率的分布规律。 相似文献
6.
基于二次场二维起伏地形MT有限元数值模拟 总被引：4，自引：1，他引：3
通过计算二次场来进行二维大地电磁数值模拟;导出了二维大地电磁二次场的微分方程,利用有限单元法来解微分方程;对矩形网格进行对角线的二次剖分,更容易且真实地模拟起伏地形。对几个典型模型进行了试算,与前人总场法的计算结果做了比较,两者视电阻率曲线一致,证明本文算法是正确的;通过2个简单的算例说明复杂地表下2种极化模式的MT观测资料都有明显的异常,视电阻率在TM模式下比TE模式更易受地形影响,TE模式下视电阻率曲线形态与地形呈"正相关",TM模式下反之。 相似文献
7.
高密度电阻率法在岩溶区找水具有较好的效果,但岩溶区地形复杂,高密度电法受地形起伏影响较大。为了提高岩溶区高密度电法找水的勘探效果,掌握地形起伏对高密度电法影响的规律,本文采用高密度电阻率法微测系统,利用铜柱体模拟地下岩溶管道,研究地形起伏对高密度电法探测岩溶管道的影响。研究结果表明:高密度联合剖面法及高密度α2装置受谷地地形影响较大,容易在谷地中间产生虚假的低阻异常,矩形AMN和滚动MNB装置受地形影响较小;受谷地地形影响,高密度电法能勘探的有效异常体深度减小;同一电极距下,平地模型及谷地模型反演异常体顶部埋深均小于真实异常,横向宽度均大于真实异常,平地模型反演异常体的顶部埋深大于谷地模型反演异常体,横向宽度小于谷地模型反演异常体。 相似文献
8.
为了模拟实际的起伏地形,在前人研究的基础上,单元网格设计为任意四边形网格,单元内的场值双线性变化,采用了高斯数值积分法计算单元系数矩阵,并给出二维大地电场的辅助场表达式。通过模型算例表明,模拟结果与解析解的均方误差在1%以内,地形模型与前人的模拟结果相吻合。分析对比了两个不同坡度模型对视电阻率和相位的影响。采用任意四边形网格剖分法,降低了编程难度,可以方便地适应野外地形的起伏变化。 相似文献
9.
为解决大角度起伏地形即陡峭地形(如河流阶地、冲沟、堤坝、山谷、山脊等)对实测视电阻率的影响问题,这里优化了三维有限元数值模拟中的单元积分形式,使之适合大倾角地形条件下对电阻率的数值模拟.为验证方法的正确性,采用保角变换计算了相近地形在均匀电流场情况下,陡峭地形引起的视电阻率变化特征.计算表明,对于缓倾斜地形,有限元数值... 相似文献
10.
射频大地电磁法(Radio-magnetotelluric, RMT)是浅地表电磁勘探的重要手段之一。由于勘探频段为10～300 kHz,其电磁场传播受地下介质介电常数的影响较大。传统准静态条件下的电磁响应严重制约RMT正演模拟精度,并进一步影响反演成像分辨率。针对这一问题,提出了一种基于Julia并行加速的全电流RMT电磁响应数值模拟方法,利用Julia的分布式计算将各频点的计算发送到不同进程进行求解,从而达到加速计算的目的,同时在计算中考虑位移电流的影响,提升正演的模拟精度。通过计算几种典型高阻/高介电模型的RMT响应,分析并总结了位移电流对射频段电磁场视电阻率及相位响应的影响规律。数值模拟结果表明：当浅部存在高阻覆盖层时,基于准静态假设条件计算的RMT视电阻率和相位响应偏高,且频率越高、覆盖层电阻率越大,响应偏差越大;对于煤炭采空区模型,RMT法能有效反映异常体位置,但忽略位移电流会在采空区及其附近引起较大计算误差;起伏地形的算例表明地形会覆盖地下异常体的RMT数值响应,尤其是地形拐角处;2种不同规模的并行算例对比证明了并行算法的高效性,且随着求解问题规模增大,并行算法效率也随之... 相似文献
11.
《物探化探计算技术》2017,(3)
在山区进行音频大地电磁测深(AMT)时,电磁场易受地形起伏影响发生畸变。为研究起伏地形对AMT数据的影响,采用WinGlink软件基于有限差分算法的正演模块,建立二维起伏地形和平地地电模型,并对比分析了两个理论模型的正演响应。结果表明:TE模式的视电阻率曲线在高频段受干扰程度相对较大,在低频段影响减小甚至消失;TM模式的视电阻曲线变化正好相反;二者相位受地形影响则较小。对含异常体的带地形地电模型进行正演计算,在正演响应中加入2.5%的随机误差模拟实测数据,利用WinGlink软件基于非线性共轭梯度(NLCG)算法的反演模块,分别进行带地形和不带地形的TE、TM和TE+TM模式的二维反演。各反演结果的对比表明,带地形的二维反演结果明显优于不带地形的反演结果,能够较好地反映出异常体的位置、形态及电阻率值。依据得到的认识,对青海祁连隧道工程的AMT实测数据进行带地形二维反演,取得了理想的应用效果。 相似文献
12.
高密度电阻率勘测方法分辨率研究与探讨 总被引：2，自引：0，他引：2
随着高密度电阻率法在工程中的应用日益广泛,对分辨率的影响因素分析研究也受到重视。论文以均匀半空间局部直立异常体为探测目标,建立多个基于温纳高密度电法装置的地电模型,采用有限元等方法进行正反演数值模拟分析,对高密度电阻率法探测分辨率的影响规律进行探讨。分析表明,地形起伏会引起异常体的位置与形态发生畸变和位移,分辨率降低; 电阻率差异增大,分辨率提高; 深径比增大,分辨率降低; 当异常体的水平范围一定时,探测分辨率随纵深的增大而提高。 相似文献
13.
地下管线是城市的重要设施,承担着能源输送、信息传递等功能,为城市生活提供了便利和保障。探地雷达(ground penetrating radar, GPR)作为一种高分辨率、高精度、非开挖、非破坏性的探测技术,在管线测量中具有巨大的优势。然而地表地形复杂,多有起伏,对于GPR探测地下管线有很大的影响。因此,本文采用有限单元法对地下管线探测进行了数值模拟,该方法可以与非结构网格结合,更好地拟合地表起伏地形;此外,介绍了如何进行高度校正,将所得剖面数据与地形相吻合,更容易分析异常体特征。最后开展了两个数值实验,分析了起伏地表对于不同埋深、不同间距、不同材质及不同填充物管线探测的影响,为GPR数据解释提供理论基础。实验结果表明,因地表起伏原因,波形和反射波能量将发生畸变,并不能作为判断管线信息的唯一依据。因此,需进行高度校正,利用双曲线的顶点来判断管线的埋深、材质等信息。 相似文献
14.
复电阻率法以岩、矿石的频谱参数特性为勘探依据,多种参数的组合解释能为评价异常源的性质提供更为丰富的信息。但其反演问题一直没有得到较好地解决,已经严重地影响和制约了该方法的应用和发展。本文由频率域电磁场的偏微分方程出发,结合Cole-Cole数学模型进行公式推导,并利用有限单元法实现了复电阻率2.5维正演模拟,能够计算起伏地形条件下包含电磁效应的复电阻率响应。在此基础上,提出了利用多个排列的视电阻率和视相位数据的复电阻率反演方法,并在反演方程中加入了Occam法的光滑模型约束,能够对复电阻率参数直接进行大规模的全区反演,有效减弱了SIP反演中欠定性严重的问题。通过理论模型的反演算例表明,该方法能够快速、准确地反演出二维地质断面上所有单元的4种复电阻率参数,从而验证了该反演算法的正确性与稳定性。 相似文献
15.
16.
水下直流电阻率法数值模拟 总被引：5，自引：0，他引：5
利用有限单元法对水下电阻率测深进行数值模拟。与地面测量不同的是，水下电极所处的深度、水体本身和水底地形都会给电阻率数据造成较大影响，给数据的解释带来困难。若先通过计算，分析不同电阻率装置在水下测量时的视电阻率与水深的关系，再用有限单元法对水下的三维地电断面进行数值模拟，讨论水底地形对视电阻率的影响．即可用比较法消除地形影响。 相似文献
17.
《中国煤炭地质》2019,(7)
倾斜煤层由于受构造挤压,往往地表形成比较陡的单面斜坡,高密度电阻率法在勘探陡倾斜煤层采空区过程中受地形及采空区充填状态的影响,数据反演的多解性给资料解释造成一定的困难。根据倾斜煤层采空区充填状态及断层破碎带设计模型,对模型进行的高密度电阻率法正反演结果表明,电阻率值受背景电阻率值影响较大;高密度电阻率法探测到的倾斜煤系地层采空区实际位置和空间分布范围及形态与实际有一定差异;对断层或采空区的探测,其浅部分辨率好,深部分辨率较低,异常体倾向大于实际目标地质体走向;受体积效应影响,异常响应范围往往大于实际采空区范围。实例证明了高密度电阻率法在倾斜煤层采空区勘探中是可行的。 相似文献
18.
井地电阻率法是在井中供电地表测量的一种方法。在井地电阻率法电法实际工作中电流源通常采用点状电流源和线状电流源,其中点状电流源适用于裸眼井中供电的情况,而线状电流源则适用于套管井供电的情况。结合野外实践工作,主要研究以下几个方面并得到了相应的结论:(1)电流源为点状时,针对电流源分别在异常体的上面、里面和下面3种模型进行正演数值模拟计算,结果表明地表的视电阻率响应变化很大;(2)电流源为垂直线源时,针对电流源分别在异常体的上部、延伸至异常体中间以及贯穿异常体3种模型进行正演数值模拟计算,结果表明地表视电阻率响应的幅值变化逐渐变大;(3)电流源为倾斜线源时,线源倾斜角度不断加大对应地表视电阻率的影响较大,随着线源倾斜角的加大,地表响应的幅值变小。这些结论为野外工作中设置井地电法仪器参数以提高其观测精度提供了理论依据,为初步判断异常体的埋深做出定性的分析,同时为井地电法三维反演奠定了基础。 相似文献
19.
《物探化探计算技术》2021,(2)
航空大地电磁法以其机动灵活、效率高等优点,适用于在地势崎岖的偏远环境中开展快速普查作业。在航空大地电磁探测中,起伏地形会对倾子响应产生影响,另外由于倾子是无量纲传输函数,只包含电性结构的空间变化信息,而与电阻率的数值水平无关,单独利用倾子数据进行反演时,反演多解性更大。为解决上述问题,首先采用四边形网格模拟地形起伏;其次在倾子数据的反演目标函数中加入少量地面大地电磁数据拟合项,克服只拟合倾子的多解性缺陷;最后利用二维地电模型,验证本文有限元算法的正演精度,并分析起伏地形对倾子响应数据的影响。针对倾子数据缺乏背景电阻率信息这一问题,开展了带地形的倾子反演和倾子与地面大地电磁数据反演研究。研究结果表明:向倾子数据中加入少量地面大地电磁数据进行反演,可有效降低单独倾子反演对初始模型的依赖。 相似文献
20.
频率测深二维地形影响的边界元素法正演模拟 总被引：5，自引：0，他引：5
本文用边界元素法模拟二维地形对频率测深的影响，给出了典型地形赤道偶极装置ＥＸ分量的纯地形影响视电阻率曲线，列举了用此值法进行地形改正的实例。 相似文献