共查询到20条相似文献，搜索用时 109 毫秒
1.
井间电磁波层析成像中的高精度时间域正演计算 总被引：1，自引：0，他引：1
为提高井间电磁波层析成像的分辨率, 正演计算中需保持较高的计算精度.本文正演计算中采用时间域伪谱(PSTD) 法模拟井间电磁波的传播.该算法在时间域有限差分(FDTD) 法的基础上采用快速傅立叶变换(FFT) 计算麦克斯韦方程中的空间导数.该算法除了具备时间域有限差分法的优点外, 在计算条件完全相同的情况下, 计算精度明显高于时间域有限差分法.时间域伪谱法的正演计算为井间电磁波高分辨率层析成像奠定了重要基础. 相似文献
2.
为提高电磁波场数值计算精度,对时间域有限差分法与时间域伪谱法进行了对比研究。时间域有限差分法是一种目前流行的电磁场时域数值计算方法,已被广泛应用于求解与时间有关的偏导数方程。对于大规模数值计算,时间域有限差分法需要较多的内存空间。时间域伪谱法基于时间域有限差分法,该方法使用快速傅里叶度换来计算麦克斯韦方程中的空间导数。由于傅里叶变换的准确性,时间域伪谱法使计算精度提高,数值计算时所需的格子数大大减少,这极大地节省了计算机内存空间,适合于大规模正反演问题的数值计算。 相似文献
3.
为了实现缓变瞬态电磁场的时域有限差分(FDTD)模拟,引入了虚拟位移电流,使FDTD随时间步的递推能够进行。在以往较简单模型的计算中,利用时域瞬变场传播需要时间的性质,可将激励源化为初始条件代入,但是对于复杂3D模型,原有的方法不再适用,应当重新考虑源和边界条件的处理方法。由于实际源带有“斜坡”,故可采取将激励源模拟注入的方式;地下边界根据电磁波在导电媒质中快速衰减的特性设定,将改进的完全匹配层用于空中边界。分析表明,对于由空中和地下区域非均匀网格离散引起的时间稳定性问题,交替方向隐式FDTD是可行的解决办法。 相似文献
4.
时间谱电阻率(TSR)法的二维正演, 是计算可极化二维地电构造上三维电流偶极源的电场瞬变响应, 属所谓2.5维时间域电磁场数值模拟问题, 是目前国际上未妥善解决的计算地球物理疑难问题.针对现有算法的局限性, 建立了新的算法, 其特点是: (1)采用二维有限单元算法, 在矩形网格中增加两对角线形成三角网格剖分, 同时用高斯消元法消除矩形网格中心结点的待求未知量.这样, 既可较准确地模拟任意二维复杂地电断面, 又可节省计算量.(2)采用直接计算二次场的新算法, 只需计算电场和磁场沿地电构造走向两个分量的一次场, 因而, 不但计算精度较高, 而且不显著增加计算量.(3)采用G -S变换法作逆拉氏变换, 并利用拉氏变换延迟定理在倍增的时间间隔中插值, 从而实现对密集采样时间瞬变过程的快速计算.(4)能对可极化和导电大地(即同时包括IP和EM效应). 相似文献
5.
为了研究频率域海底电磁法在油气勘探中的可行性,根据海底电磁法的观测特点,建立海底三维地电模型,采用频率域FDTD三维正演模拟程序,计算不同参数时频率域电磁响应,获得了不同模型、不同频率下电磁场的响应。 相似文献
6.
为了更好地分析探地雷达(GPR)叠前逆时偏移算法对实际地下随机分布介质的成像效果,采用随机过程的谱分解和混合型自相关函数理论构建了不同自相关长度的GPR随机介质模型。基于时域有限差分法(FDTD)构建了GPR叠前逆时偏移成像算法,其中FDTD用于计算正传和反传电磁波场,归一化互相关成像条件用于获取逆时偏移成像剖面。在此基础上,利用该算法对两个随机介质模型的多偏移距正演数据进行计算,并与相应背景介质为均匀介质的逆时偏移结果进行对比。结果表明:电磁波在随机介质中散射强烈,反射波扭曲变形、不连续,形成了明显的随机扰动,致使逆时偏移成像剖面的空间分辨率更低,低频噪声更强;自相关长度是影响随机介质中异常体成像效果的主要因素,自相关长度越小,异常体的成像越清晰,自相关长度越大,异常体的成像效果越差,且不易被识别。 相似文献
7.
为了研究钻孔雷达对金属矿的响应特征,利用时域有限差分法（FDTD）对钻孔雷达探测进行了数值模拟。通过对不同形态的金属矿体,包括等球状体、板状体和一个实际矿体的模拟,得到大量模拟结果。球状目标体在剖面上表现为双曲线形态,通过偏移处理可以精确确定目标体的位置。板状矿体被钻孔穿过时,直达波位置发生错动并变弱,反射同相轴向上下延伸;板状矿体未被钻孔穿过时,反射同相轴呈线状。对于实际矿体,不同的钻孔表现出不同的特征,结合左、中、右3个钻孔的探测结果可以判断出目标体的形态,并且可以判断出此目标矿体具有低阻、高介电常数的特性。 相似文献
8.
9.
三维频散介质中地质雷达信号的FDTD数值模拟 总被引：10，自引：0，他引：10
考虑频散介质的电磁波传播,引入随频率变化的电位移矢量D,并对电场强度E和电位移D进行标准化。开发出一种计算介电常数随频率变化满足Debye关系的频散介质中麦克斯韦方程的时间域有限差分（FDTD）解法,并编写了相应的计算原代码。在边界部分引入假象的介电常数和磁导率,实现了边界的无反射吸收,避免了Berenger完全匹配层中场分裂过程,提高了计算效率。计算实例可以看出频散介质中波的传播情况。 相似文献
10.
地质雷达正演中的频散压制和吸收边界改进方法 总被引：5，自引：0，他引：5
从麦克斯韦方程组出发,建立了地质雷达的时域有限差分法（ＦＤＴＤ）数学模型,导出了理想 散关系和超级吸收边界条件。理想频散关系 考虑了ＦＤＴＤ法的收敛性和稳定性,也考虑了高频电磁波在Ｙｅｅ氏网格中的传播特点;超吸收边界条件则用磁场分量来提高电场分量精度。数值试验表明,理想频散关系能真实地反映雷达波在地下介质中的传播规律,超吸收边界条件能有效减小截断边界的伪反射,提高正演精度。将之应用于实际计算,取得 相似文献
11.
根据不同时域电磁场的频率特性,从频率域到时间域的转换一般有二种方法,即线性滤波算法(余弦变换)和快速傅氏变换(FFT)。利用高密度采样的线性滤波算法计算余弦变换,对层状介质的甚早期瞬变电磁进行模拟,并与均匀半空间的闭合解析式结果进行对比,证明了该方法的有效性。另外,从层状介质的高频电磁场响应出发,利用FFT算法,有效地模拟了不同频率和地电参数的GPR射线,以及层状介质的剖面图。甚早期瞬变电磁和低频GPR的频率,都跨越了传统电磁法中所忽略的中间频率。这里的主要目的就是通过对比选择有效的模拟方法,利用模拟结果来分析该频段时域电磁场的特征。 相似文献
12.
射频大地电磁法(Radio-magnetotelluric, RMT)是浅地表电磁勘探的重要手段之一。由于勘探频段为10～300 kHz,其电磁场传播受地下介质介电常数的影响较大。传统准静态条件下的电磁响应严重制约RMT正演模拟精度,并进一步影响反演成像分辨率。针对这一问题,提出了一种基于Julia并行加速的全电流RMT电磁响应数值模拟方法,利用Julia的分布式计算将各频点的计算发送到不同进程进行求解,从而达到加速计算的目的,同时在计算中考虑位移电流的影响,提升正演的模拟精度。通过计算几种典型高阻/高介电模型的RMT响应,分析并总结了位移电流对射频段电磁场视电阻率及相位响应的影响规律。数值模拟结果表明：当浅部存在高阻覆盖层时,基于准静态假设条件计算的RMT视电阻率和相位响应偏高,且频率越高、覆盖层电阻率越大,响应偏差越大;对于煤炭采空区模型,RMT法能有效反映异常体位置,但忽略位移电流会在采空区及其附近引起较大计算误差;起伏地形的算例表明地形会覆盖地下异常体的RMT数值响应,尤其是地形拐角处;2种不同规模的并行算例对比证明了并行算法的高效性,且随着求解问题规模增大,并行算法效率也随之... 相似文献
13.
井间地震层析成象高精度算法及其在环境工程勘探中的应用 总被引：2，自引：0，他引：2
井间地震层要成象在环境工程勘探中的应用前景广阔，走时层析成象的射线追踪方法和在型稀疏矩阵的方法的方法是其中的关键。文中讨论了ＬＳＱＲ算法系列，数值模拟证明了它的优越性，最后应用于鹏飞大厦场地，取得了良好的效果。 相似文献
14.
局部异常体井间地震走时层析成像研究 总被引：1，自引：0，他引：1
通过数值模拟试验 ,验证了利用井间地震走时层析成像方法技术可以分辨线性尺度为1个像元宽度的目标 ;讨论了表层低速带对基岩中异常体地震走时层析成像分辨率的影响 ;给出了高分辨率探测局部异常体井间地震走时层析成像方法技术 ,即最佳的观测系统、有效的数据采集技术、精细的数据预处理和优化的正反演算法。通过实际工程应用表明 ,该方法技术具有良好的效果 相似文献
15.
钻孔雷达相比传统的地球物理方法具有更大的探测范围,可得到更多的地层信息。在地质雷达正演模拟的方法中,时域有限差分法能提高探测的效果和解释的准确性。将两者相结合,通过单孔反射测量对金属矿进行数值模拟及对结果进行分析,说明了单孔反射测量不仅可以确定矿体的位置,而且对矿体形态也具有很好的识别性。 相似文献
16.
模型的建立与描述是地震资料分析与处理的基础,合理的模型参数化方式对于地震正反演的各个方面都有很好的效果,其中射线追踪在正演模拟、层析成像及偏移等研究领域中都占有很重要的地位。本文针对地震波走时与射线路径计算,设计了一种将离散模型连续化的最小二乘模型参数化方法,分别对速度模型和走时模型进行模型参数化,然后进行梯度速度模型验算及误差分析。计算结果表明,该模型参数化方法使走时计算精度从10-2提高到10-3,使射线路径的计算精度提高了14.33%。最后通过经典Marmousi模型和Sigsbee 2A模型实例验证,充分证明了该模型参数化方法的有效性及普遍适用性。 相似文献
17.
巷?孔瞬变电磁法在存在明显电导率各向异性的勘探区会产生较大的解释误差。基于时域有限差分算法,通过引入电导率各向异性张量构建控制方程、将矩形回线源电流密度加入Maxwell方程安培环路定理实现任意电流源的加载,以差分代替微分对控制方程进行离散,实现巷?孔瞬变电磁三维模型正演计算。在与解析解对比验证算法计算精度的基础上,构建全空间模型、层状模型和三维块状模型进行正演并分析电导率各向异性对巷?孔瞬变电磁三分量响应的影响程度与方式。结果表明:垂直轴电导率对巷?孔瞬变电磁三分量响应基本没有影响,水平轴电导率对巷?孔瞬变电磁三分量响应影响较大,其中,?Bx/?t响应主要受y轴电导率影响,?By/?t响应主要受x轴电导率影响;通过三分量响应的形态和幅值的相互关系可以辨别各向异性介质所在方位及主轴电导率方向;当异常体所处全空间介质为电导率各向异性时,异常体产生的异常响应会被全空间介质电导率各向异性产生的异常响应所淹没。电导率轴向各向异性特征在巷?孔瞬变电磁法解释过程中不可忽略,研究成果为巷?孔瞬变电磁法各向异性解释提供指导,也为各向异性反演提供参考。 相似文献
18.
起伏地形频率域可控源电磁二维快速正反演 总被引：1，自引：0，他引：1
起伏地形频率域可控源电磁正反演研究受到众多学者的关注和重视,目前普遍实用的为未考虑地形的二维正反演算法程序。本文基于积分方程法正演和多场源、多频率、对比源反演算法,开展起伏地形条件下频率域可控源电磁二维快速正反演研究。通过引入层状参考模型、将地形与目标体整体作为异常场剖分区域,实现起伏地形可控源电磁正反演计算。针对包含地形在内的大尺度剖分区域引起的大型计算代价问题,采用快速傅里叶算法提高正反演计算效率。通过与现有正反演算法进行模型算例对比,验证了本文采用方法的可行性与有效性。 相似文献
19.
PPP-BOTDA分布式光纤传感技术及其在试桩中应用 总被引：3，自引：1，他引：2
布里渊时域分析技术（brilliouin optical time domain analyzer, BOTDA）是近年来在国际上研发成功的一种分布式光电传感技术,该技术已开始广泛应用于结构体和工业设备等的温度、应变检测。而新一代分布式应变测量技术--脉冲预泵浦-布里渊时域分析（pulse-prepump-brilliouin optical time domain analyzer, PPP-BOTDA）,实现了10 cm的空间分辨率和±7.5 με的应变测量精度。文中重点介绍了PPP-BOTDA工作原理以及技术验证,并在国内首次将PPP-BOTDA技术应用于管桩桩身应变测量,与同时进行的滑动测微计实测数据进行对比,结果表明PPP-BOTDA可作为一种新型的桩基测试技术加以推广 相似文献
20.
采用 Gaver-Stehfest( G-S)方法模拟水平电偶极源在不同介质模型中地表面的负阶跃波电磁场响应 ,然后引入逆时序偏移的观点 ,把这些电磁场响应作为实际观测值 ,利用 Stratton-Chu积分公式将其向下偏移 ,对偏移后的电磁场值进行归一化处理并对场值求偏导数 ,重建地下拟地震深度剖面图 ,实现时间域电磁场的偏移成像。选用了均匀介质、含良导和高阻薄层水平均匀层状介质模型以及高阻柱体、断层等模型分别进行计算 ,成像的结果表明能够得到较高质量的地下地电构造图 ,从重建的深度剖面图中可以识别出地下反射界面的位置 ,该图比时间剖面直观。 相似文献