共查询到20条相似文献，搜索用时 250 毫秒
1.
张卜文 《物探化探计算技术》2018,(3)
研究用于分析层状空间目标电磁散射特性的时域有限差分(FDTD)方法。从二维Maxwell方程出发,导出带有入射角的一维修正Maxwell方程组及辅助方程,来计算二维总场-散射场(TF/SF)垂向边界上电磁场,从而在层状空间中实现直接时域形式的均匀平面波入射;为避免复杂的Sommerfeld积分,用互易原理得到简化的外推算法。用该方法对半空间目标电磁散射场和层状有耗空间中的线电流辐射场进行计算,在半空间近场区域得到了良好的均匀平面波,并且远场与已有文献结果和解析解结果都吻合较好,表明本算法和程序是正确和有效的。用该算法对多层状空间中的隧道、隧道内通车辆等目标以及隧道和车辆在无耗层状空间中的电磁散射特性进行研究。结果表明,金属车体对远区散射场贡献较大,且介质层的损耗特性影响着其对车体的屏蔽能力。 相似文献
2.
有耗媒质中线天线激励下三维异常体的散射 总被引：1，自引：0，他引：1
本文研究均匀有耗媒质中,在半波长偶极天线的场激励下,多种形状和参数的异常体的散射特性。用矩量法分别求解良导电球体散射的磁场积分方程和有限导电的长方形异常体散射的电场积分方程。讨论了与数值法有关的近似、误差和检验等问题。通过数值计算得出媒质电参数改变、散射体大小、形状以及空间位置改变时散射场变化的规律。水中金属球和介质块散射的测量结果与数值计算结果吻合较好,验证了数值模型的正确性。 相似文献
3.
《岩土力学》2020,(1)
为研究横观各向同性(TI)介质的凸起地形对平面qP-qSV波的散射问题,采用间接边界元法进行数值模拟并且分别给出该问题在时域和频域内的解答。在求解过程中将计算模型分解为开口层状半空间域和闭合凸起域,同时将波场分解为自由波场和散射波场,自由波场可通过直接刚度法求得,散射波场由斜线荷载动力格林函数来模拟。频域结果表明,TI介质参数的改变导致了凸起和层状半空间动力特性的改变,进而导致了凸起与层状半空间动力相互作用机制的改变,使得不同TI介质参数对应的地表位移幅值显著不同。而时域结果表明,TI介质参数与波的传播方向均对qP波和qSV波的在凸起地形中的传播有着显著影响。另外,时域结果更为清晰的地展示入射波、反射波和散射波的传播过程。 相似文献
4.
《岩土力学》2021,(1)
为研究横观各向同性(TI)介质的凸起地形对平面qP-qSV波的散射问题,采用间接边界元法进行数值模拟并且分别给出该问题在时域和频域内的解答。在求解过程中将计算模型分解为开口层状半空间域和闭合凸起域,同时将波场分解为自由波场和散射波场,自由波场可通过直接刚度法求得,散射波场由斜线荷载动力格林函数来模拟。频域结果表明,TI介质参数的改变导致了凸起和层状半空间动力特性的改变,进而导致了凸起与层状半空间动力相互作用机制的改变,使得不同TI介质参数对应的地表位移幅值显著不同。而时域结果表明,TI介质参数与波的传播方向均对qP波和qSV波的在凸起地形中的传播有着显著影响。另外,时域结果更为清晰的地展示入射波、反射波和散射波的传播过程。 相似文献
5.
起伏地形频率域可控源电磁二维快速正反演 总被引：1，自引：0，他引：1
起伏地形频率域可控源电磁正反演研究受到众多学者的关注和重视,目前普遍实用的为未考虑地形的二维正反演算法程序。本文基于积分方程法正演和多场源、多频率、对比源反演算法,开展起伏地形条件下频率域可控源电磁二维快速正反演研究。通过引入层状参考模型、将地形与目标体整体作为异常场剖分区域,实现起伏地形可控源电磁正反演计算。针对包含地形在内的大尺度剖分区域引起的大型计算代价问题,采用快速傅里叶算法提高正反演计算效率。通过与现有正反演算法进行模型算例对比,验证了本文采用方法的可行性与有效性。 相似文献
6.
双频激电法感应耦合的数值计算与分析 总被引：1，自引：1，他引：0
在吸收前人工作的基础之上,针对均匀半空间情况下,对双频激电测量中的中梯装置电磁耦合规律与特点进行了研究和计算。在该数值的计算中：①以Maxwell方程为基础,从矢量势及边界条件出发,阐述了各向同性、线性均匀半空间时水平电偶极子源产生的谐变电磁场;②针对中间梯度常规测量装置中的平行线布线方式和共线（供电线呈“n”型）布线方式,推导了电磁感应耦合效应的具体计算公式,计算并分析了EM效应引起的视幅频率随供电极距、工作频率、大地电阻率和偏移距等因素的变化规律。最后,对这二种布线方式的电磁耦合效应进行了比较,平行线布线方式的电磁耦合效应比共线布线方式的电磁耦合效应强烈。 相似文献
7.
利用有限元法实现了任意方向偶极子源在二维介质中频率域电磁响应的数值模拟,研究了波数取值对模拟结果的影响.通过对构造走向的Fourier变换,将全三维电磁问题,转化为一系列二维问题,并在波数域求解,极大地减小了计算工作量,导出了波数域耦合适用于二维电性介质中任何方向电或磁偶极子响应计算的电磁场方程.针对每个给定的波数,上述耦合电磁场方程用等参有限元方法在x-z平面内求解.采用Fourier逆变换,将波数域解积分,得到空间域电磁场.针对电磁模拟计算中,源点的奇异性,采用具有一定面积的伪δ函数表达源电流分布,使数值解精度得以提高.另外,采用等参有限元,使地下复杂地质体得到准确表达.利用不同波数值对均匀介质与层状介质的模拟结果与解析解的对比,验证了算法的正确性与精度.利用层状介质模型的解析解与数值计算结果的对比,分析了波数的优化取值范围及取值点数对数值模拟结果的影响,考察了算法对非均匀介质的适应性. 相似文献
8.
平面波的传播问题通常可以归结为一维波动方程的定解问题。在非均匀介质中,即使简单的一维波动方程也需要借助于数值方法获得近似解。3层5点古典差分格式是计算偏微分方程一种常用算法,作为一种显式迭代格式,需要满足稳定性条件 ,其中 为波速, 为空间采样间隔, 为时间采样间隔。当 时, ,古典差分格式达到临界稳定状态。在这种情况下,平面波在 时间内的传播距离恰好等于空间采样间隔,差分格式真实地反映了平面波的传播原理,因而可以得到一维波动方程的精确解。但是,由于在非均匀介质中存在不连续的波阻抗界面,此方法不适于计算非均匀介质的波场。为了将临界稳定情况下的古典差分格式推广应用至非均匀层状介质,提出了一种能够处理波阻抗界面的有限差分格式,并应用傅里叶分析法得到其稳定性条件。模型算例验证了此算法的正确性。 相似文献
9.
《物探与化探》2015,(5)
采用三角网格剖分的有限元法,研究了2.5维航空瞬变电磁法正演模拟问题。利用时频变换数值方法将时间域电磁场转换到拉氏域,再利用傅里叶变换将三维问题降维变为2.5维问题,然后由有限元法求解得到拉氏域二维电磁场,逆拉氏变换后得到时间域航空瞬变响应。为了回避正演模拟中总感应磁场在场源处的奇异性问题,采用异常场算法,场源响应通过在微分方程中施加背景电磁场实现。由于瞬变电磁信号具有较大的动态范围,而且需要经过两次正、逆拉氏变换和傅里叶变换,每个环节的计算精度和速度要严格控制在较高的水平上,否则积累误差会非常大。模型计算表明均匀大地和层状大地模型解析解与数值解吻合很好。这证明该算法是正确可行的,可作为研究二维复杂地质体的方法手段。 相似文献
10.
11.
鬃变电磁场信号在地下的扩散及地面上的时域响应特性 总被引：1，自引：0，他引：1
瞬变电磁测深是近年来发展较快的电磁勘探方法，是实现精确勘探的一个突破点，本文应用时域有限差分算法分析了瞬变场在二维地电模型情况下低阻体对“烟圈”的吸收现象、瞬变场的滞后现象，对进一步的时频分析和高维反演提供了正演基础。 相似文献
12.
拟解析近似方法是一种求解积分方程的一种近似方法,它可以处理强散射或者大扰动的电磁散射问题,在计算过程中避免了传统微分数值方法解决问题时所遇到的大型矩阵或大型代数方程组的求解。孙建国[5]将其引入直流电场的积分方程中,并给出了求解异常电场积分方程的标量拟解析近似公式。在以前的研究中,已经验证了均匀场中异常球体的拟解析近似解的精度,这里对均匀场中的立方体异常体进行数值模拟,得到了直流电场中异常立方体模型的标量拟解析近似解。由于复杂地电模型可以用立方体的组合进行模拟,因此对立方体异常电场拟解析近似解的研究,为三维直流电场中复杂地电模型的快速正反演模拟打下了基础。 相似文献
13.
本文针对地下两孔电磁法的实际,建立多异常体(包括非均匀异常体)散射的数值模型。给出一种使计算机程序具有一定通用性的处理分块序号与坐标关系的方法。用矩量法解体积分方程,并采用合适的近似公式,计算出不同工作频率下各种异常体散射的场分布。由计算结果可得出一些对解释实测数据有参考价值的结论。 相似文献
14.
在面波层析成像和传统面波勘探中,都假定地球或浅层探测对象是一个均匀分层的多层介质模型,而实际上,介质的横向不均匀性广泛存在于从地球直到超声尺度的各种范围内.这里采用模式耦合的方法,研究了均匀半空间和二层介质中(也可以是多层的),三维非均匀体对Rayleigh面波的多次散射,计算了频率和时间域中的散射波场,并分析了散射波场的特点,为散射面波的应用提供了初步的理论基础. 相似文献
15.
巷?孔瞬变电磁法在存在明显电导率各向异性的勘探区会产生较大的解释误差。基于时域有限差分算法,通过引入电导率各向异性张量构建控制方程、将矩形回线源电流密度加入Maxwell方程安培环路定理实现任意电流源的加载,以差分代替微分对控制方程进行离散,实现巷?孔瞬变电磁三维模型正演计算。在与解析解对比验证算法计算精度的基础上,构建全空间模型、层状模型和三维块状模型进行正演并分析电导率各向异性对巷?孔瞬变电磁三分量响应的影响程度与方式。结果表明:垂直轴电导率对巷?孔瞬变电磁三分量响应基本没有影响,水平轴电导率对巷?孔瞬变电磁三分量响应影响较大,其中,?Bx/?t响应主要受y轴电导率影响,?By/?t响应主要受x轴电导率影响;通过三分量响应的形态和幅值的相互关系可以辨别各向异性介质所在方位及主轴电导率方向;当异常体所处全空间介质为电导率各向异性时,异常体产生的异常响应会被全空间介质电导率各向异性产生的异常响应所淹没。电导率轴向各向异性特征在巷?孔瞬变电磁法解释过程中不可忽略,研究成果为巷?孔瞬变电磁法各向异性解释提供指导,也为各向异性反演提供参考。 相似文献
16.
斜阶跃波层状介质航空瞬变电磁响应数值计算 总被引：1，自引：0，他引：1
目前,航空瞬变电磁方法数据解释主要采用阶跃波均匀半空间模型计算视电导率值,而实际航空电磁系统发射波形的下降沿多为斜阶跃,导致解释时计算的视电导率值存在较大偏差。为此,笔者研究了航空瞬变电磁系统发射电流为斜阶跃时的电磁响应数值计算,将发射波形进行拉氏变换,利用G-S逆拉氏变换与241点汉克尔变换相结合的方法,实现斜阶跃波关断后的层状大地模型航空瞬变电磁响应计算;并对均匀半空间和层状大地模型下,不同关断时间和不同飞行高度对电磁响应的影响进行分析。得出结论:不同关断时间,关断后取样延时2 ms时,均匀半空间电磁响应的平均偏差为27.78%,三层模型的平均偏差为32.16%;当飞行高度从20 m增加到60 m时,均匀半空间和三层模型的感应电动势分别减小了43.6%和83.2%。 相似文献
17.
将无网格局部Petrov-Galerkin算法用于大地电磁二维正演。介绍了该方法的基本原理;从大地电磁二维边值问题出发,利用子域法详细推导了与之对应的局部Petrov-Galerkin弱式方程,并用高斯积分法将其离散化。论述了无网格局部Petrov-Galerkin法较无单元Galerkin法及有限元法的优缺点,最后通过二维模型的计算验证了算法的有效性。 相似文献
18.
为了研究地-井瞬变电磁响应特征,获得不同井位低阻薄板异常体的响应规律:首先采用交错网格有限差分技术离散二次场满足的频率域赫姆霍兹方程;然后结合虚框叠加等效和虚拟界面法将发射源和接收位置扩展到任意层位,解决了全空间背景格林函数的计算问题;之后利用MUMPS求解器求解频率域二次场,再经过余弦变换,获得井中任意位置的瞬变电磁响应;最后采用三维模型对本文算法的正确性进行验证,设计了均匀半空间和嵌入低阻薄板的三维模型,获得垂直磁场的三维分布,分析均匀导电半空间中低阻薄板对地-井瞬变电磁三维响应的影响特点.结果表明:本文三维地-井瞬变电磁计算方案的计算精度与前人基本相同;水平导电薄板的存在主要影响板体附近及穿过薄板的井中瞬变响应,最大的特点是在异常体位置附近的中期时间道响应出现变号现象.本文的研究为定性解释地-井瞬变电磁法异常提供一个技术手段,也为地-井瞬变电磁三维反演奠定基础. 相似文献
19.
在井下瞬变电磁探测的时域有限差分(简称FDTD)计算中,将较薄的煤层用局部网格细化,将巷道中的金属支架、铁轨作为细导线处理,使要模拟的结构很小,需要计算场的物理空间又不能相对地减小情况下的计算仍然能够进行. 相似文献
20.
时域瞬变场电磁场有限差分法 总被引：1，自引：2，他引：1
模拟电磁场在地下介质中的传播规律是理解地下介质响应的重要手段，直接应用有限差分法在时间域对二维、三维瞬变场进行数值分析是一种有效方法。应用反映电磁场基本规律的麦克斯韦方程组，导出时域场量的齐次扩散方程，对所研究的空间区域作差分离散，源作为初始条件加入，采用合适的边界条件，利用差分方程进行计算，从而得到场量的数值结果，展现瞬变场在地下随时间扩散的全过程，对电磁场的瞬态响应有更加直观的了解。 相似文献