磁偶源2.5维瞬变电磁场全空间FDTD数值模拟

为解决模拟瞬变电磁场在3维介质中传播出现的计算量大、对计算机要求高等问题,使用了3维磁偶极源来模拟2维均匀介质和层状介质中方形和薄板低阻体的全空间响应特征。模拟中使用了时域有限差分法,采用非均匀网格和时间步长,分析了瞬变电磁场在均匀介质和层状介质中的传播规律,以及巷道和交界面对场的影响特征,为解释全空间均匀介质和层状介质TEM异常提供参考。     

考虑卷积完全匹配层数值色散的井间电磁三维正演

本文将以卷积完全匹配层为吸收边界条件的时域有限差分法应用到井间电磁的三维正演模拟中.证明了卷积完全匹配层中的数值色散会因为有效延伸因子而产生,列举常规有效延伸因子和网格间距对电磁波相速度各向异性的影响,并通过波场快照验证卷积完全匹配层中数值色散的存在;进而推广三维卷积完全匹配层中最大有效延伸因子、最大网格间距与激励源主频之间的约束,以此完善卷积完全匹配层的最优参数设置方案.在此基础上,通过正演结果中二次场的垂直磁场分量和走时来展示灵敏度分布,以此划定井间电磁勘探的优势区域,提出井间电磁正反演所需要的约束测井系列和最佳取井方案.     

考虑关断时间的回线源激发TEM三维时域有限差分正演

从麦克斯韦旋度方程出发可以直接导出瞬变电磁场扩散方程,然而扩散方程不含电场对时间的一阶导数,不能构成显式的时域有限差分方程,借鉴du Fort-Frankel有限差分离散方法引入虚拟位移电流项构建显式时域有限差分方程.对Wang和Hohmann的经典时域算法进行了两点改进:第一,通过将矩形回线源电流密度加入麦克斯韦方程组的安培环路定理方程,实现回线源瞬变电磁激发源加入;第二,在计算中考虑关断时间.第一点改进使时域有限差分方程考虑了一次场的计算,并且源的计算不再依赖均匀半空间模型响应作为初始条件,使算法能够适应表层电阻率不均匀时的三维复杂模型.由于实际观测中不可能出现阶跃电流的关断形式,第二点改进可以方便设置发射电流下降沿.采用改进的三维时域有限差分正演算法对均匀半空间模型、四类三层模型、均匀半空间中含有低阻块体模型进行了计算并分别与解析解、线性数字滤波解、积分方程解和Wang的三维时域有限差分解进行了对比验证.以H模型为例,采用建立的三维时域有限差分正演算法计算了不同关断时间的斜阶跃脉冲回线源瞬变电磁中心点感应电动势衰减曲线.以实际地质资料为基础,构建包含两层采空区的三维复杂模型,以1 μs的极短关断时间进行了复杂模型定回线源瞬变电磁响应计算,并计算了该复杂模型的视电阻率曲线.     

祁连山老虎沟12号冰川冰内结构特征分析

利用探地雷达(Ground Penetration Radar,GPR)调查冰川冰内结构和冰层厚度是一种监测冰川变化常用手段.应用麦克斯维(Maxwell)方程的二维时域有限差分(FDTD)方法,通过将模拟图像与老虎沟12号冰川的实测图像的对比,分析了雷达剖面中的几种反射特征,如冰裂隙、融洞、暖冰等.对比结果表明,GP...     

Current and Future Applications of 3-D Global Earth-Ionosphere Models Based on the Full-Vector Maxwell’s Equations FDTD Method

Advances in computing technologies in recent decades have provided a means of generating and performing highly sophisticated computational simulations of electromagnetic phenomena. In particular, just after the turn of the twenty-first century, improvements to computing infrastructures provided for the first time the opportunity to conduct advanced, high-resolution three-dimensional full-vector Maxwell’s equations investigations of electromagnetic propagation throughout the global Earth-ionosphere spherical volume. These models, based on the finite-difference time-domain (FDTD) method, are capable of including such details as the Earth’s topography and bathymetry, as well as arbitrary horizontal/vertical geometrical and electrical inhomogeneities and anisotropies of the ionosphere, lithosphere, and oceans. Studies at this level of detail simply are not achievable using analytical methods. The goal of this paper is to provide an historical overview and future prospectus of global FDTD computational research for both natural and man-made electromagnetic phenomena around the world. Current and future applications of global FDTD models relating to lightning sources and radiation, Schumann resonances, hypothesized earthquake precursors, remote sensing, and space weather are discussed.

地质雷达正演中的频散压制和吸收边界改进方法

从麦克斯韦方程组出发,建立了地质雷达的时域有限差分法（ＦＤＴＤ）数学模型,导出了理想 散关系和超级吸收边界条件。理想频散关系 考虑了ＦＤＴＤ法的收敛性和稳定性,也考虑了高频电磁波在Ｙｅｅ氏网格中的传播特点;超吸收边界条件则用磁场分量来提高电场分量精度。数值试验表明,理想频散关系能真实地反映雷达波在地下介质中的传播规律,超吸收边界条件能有效减小截断边界的伪反射,提高正演精度。将之应用于实际计算,取得     

利用时域有限差分方法计算雷电水平电场

雷电水平电场是计算电力架空线雷电感应过电压的重要参数,为了准确获得该参数值,本文利用时域有限差分方法(FDTD)对不同地面电导率情况以及距雷电回击通道不同距离处的雷电水平电场进行了计算和分析,结果表明:(1)地面电导率均匀情况下,地面电导率越大,计算距离越近,雷电水平电场波形愈加接近模式输入的电流波形,这与雷电水平电场中的辐射场分量较小有关,而增大地面电导率或者增加计算距离都会导致计算结果变小甚至波形出现双极性的特征;(2)地面电导率水平分层情况下,计算得到的雷电水平电场受上层土壤的电导率影响更大,但是下层不同电导率的土壤一定程度上改变了土壤整体的等效电导率,使得计算结果与均匀电导率情况下有所差别。     

地面建筑物对地闪回击通道附近电磁辐射环境的影响

采用有限电导率下的时域有限差分（FDTD）方法,分析了建筑物对地闪回击电磁环境的影响.结果发现建筑物对回击电场的影响较为明显,对磁场影响较小,建筑物高度h_b=50 m时对磁场的影响不大于10%.建筑物顶部电场受自身高度的影响较大,屋顶中部的垂直电场随建筑物高度的增加而增大,水平电场随建筑物高度的增加而减小.建筑物高度为20、50和100 m时,屋顶中部的垂直电场分别约为地面处的1.3、1.7和2倍,水平电场约为地面处的0.9、0.8和0.7.电导率分别取0.001和0.1 s/m时对比发现,地表处的水平电场受电导率的影响较大,电导率越大,地表处水平电场的幅值越小,此时建筑物的存在对其影响不明显.地表以上一定高度处的水平电场则受观测点位置和建筑物的影响更加明显.     

频率域海底电磁法三维正演研究

为了研究频率域海底电磁法在油气勘探中的可行性,根据海底电磁法的观测特点,建立海底三维地电模型,采用频率域FDTD三维正演模拟程序,计算不同参数时频率域电磁响应,获得了不同模型、不同频率下电磁场的响应。