期刊文献+
共找到183篇文章
< 1 2 10 >
每页显示 20 50 100
含激发极化效应的地空瞬变电磁直接时间域三维正演
1
作者 杨双燕 齐彦福 +3 位作者 李梓源 舒臣臣 孙乃泉 李貅 《地球物理学报》 北大核心 2025年第1期312-326,共15页
当使用地空瞬变电磁法对诸如铜镍硫化物矿或石墨矿等含有激发极化效应的目标进行探测时,电磁感应信号和激发极化效应相互耦合,导致地空瞬变电磁晚期观测数据出现反号现象,使用传统的数据解释技术无法处理这种信号畸变.为此,本文针对含... 当使用地空瞬变电磁法对诸如铜镍硫化物矿或石墨矿等含有激发极化效应的目标进行探测时,电磁感应信号和激发极化效应相互耦合,导致地空瞬变电磁晚期观测数据出现反号现象,使用传统的数据解释技术无法处理这种信号畸变.为此,本文针对含激发极化效应的地空瞬变电磁直接时间域三维正演算法开展研究.借助微分欧姆定律将Cole-Cole模型完整地引入到时间域控制方程中,对于在此过程中产生的分数阶导数,采用Caputo定义格式对其进行离散,从而建立含激发极化效应的时间域控制方程.然后,使用基于非结构化网格的有限元方法和后退欧拉离散格式,对时间域电场控制方程进行空间和时间离散,实现复杂介质模型含激发极化效应的地空瞬变电磁直接时间域三维正演.通过将本文三维正演结果与极化均匀半空间模型一维半解析解进行对比,验证该算法的准确性.在此基础上,本文对典型极化体模型和复杂矿体模型进行三维正演模拟,并分析极化参数对地空瞬变电磁响应的影响规律. 展开更多
关键词 电性源地空瞬变电磁法 激发极化效应 三维正演 有限元算法 CAPUTO分数阶导数
在线阅读 下载PDF
内蒙古白云鄂博主东矿三维地电模型与短偏移距瞬变电磁三维正演分析
2
作者 陈卫营 宋婉婷 +3 位作者 薛国强 雷康信 石金晶 朱昱炼 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2024年第10期3931-3945,共15页
内蒙古白云鄂博大型稀土-铌-铁矿床历经近70年开采,浅部铁矿资源已开发殆尽,亟需寻找深部资源.利用电磁法摸清主矿和东矿深部矿体连接状态及赋存特征对进一步认识白云鄂博成矿模式和寻找深部接替资源意义重大.然而,该区复杂的构造、地... 内蒙古白云鄂博大型稀土-铌-铁矿床历经近70年开采,浅部铁矿资源已开发殆尽,亟需寻找深部资源.利用电磁法摸清主矿和东矿深部矿体连接状态及赋存特征对进一步认识白云鄂博成矿模式和寻找深部接替资源意义重大.然而,该区复杂的构造、地层、岩性、地形等条件给电磁法勘探带来极大挑战.采用数值建模与模拟手段对关键观测要素进行分析有助于勘探设计和数据解释.本文以白云鄂博主东矿及其结合部位为主要研究区域,基于已有地质、钻孔等资料,构建了主东矿三维地电模型,利用时域矢量有限元算法开展了电性源短偏移距瞬变电磁法(SOTEM)三维数值模拟,并重点分析了矿坑地形对观测数据的影响程度以及不同收发装置下的深部矿体的三维电磁学响应特征和探测能力.结果表明矿坑地形会导致SOTEM电磁场响应值升高,在主矿西侧布置发射源,更有助于实现对连接部位下方铁矿体的有效探测.该工作可为今后野外实际观测中工作参数优选和数据解释提供依据. 展开更多
关键词 SOTEM 地电模型 矢量有限元 三维正演 白云鄂博
在线阅读 下载PDF
任意起伏地形下重力异常三维正演及并行计算
3
作者 戴世坤 朱德祥 +4 位作者 张莹 李昆 陈轻蕊 凌嘉宣 田红军 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2024年第2期768-780,共13页
为了进一步提高空间-波数域三维重力异常正演算法的适用范围和计算效率,本文采用任意傅里叶变换算法实现了空间-波数域三维重力异常正演,且在NVIDIA CUDA平台上进行CPU-GPU并行加速.任意傅里叶变换算法的基本思想是将二维傅里叶变换转... 为了进一步提高空间-波数域三维重力异常正演算法的适用范围和计算效率,本文采用任意傅里叶变换算法实现了空间-波数域三维重力异常正演,且在NVIDIA CUDA平台上进行CPU-GPU并行加速.任意傅里叶变换算法的基本思想是将二维傅里叶变换转化为两个一维傅里叶变换,一维傅里叶变换积分离散为多个单元积分累加和,离散单元中原函数采用二次插值形函数拟合,求出单元积分的解析表达式.相比现有的傅里叶变换算法,新方法具有采样灵活、积分精度高、计算速度快和傅里叶变换的截断效应小等优势.利用空间-波数域算法的高度并行性,采用CPU并行求解常微分方程,GPU并行计算任意傅里叶变换,实现了CPU-GPU并行加速方案,进一步提升了本文算法效率.利用常密度模型,对比数值解和解析解,结果表明本文算法正确;利用变密度模型对比了任意傅里叶变换算法与高斯快速傅里叶变换算法的计算效率与精度,在相近的数值精度下,本文算法波数选取少,效率高;测试CPU-GPU并行效果,结果表明相比CPU串行算法,CPU-GPU并行算法的计算效率大大提升,千万数量级节点数模型正演仅耗时数秒.最后利用实际地形数据进行三维重力异常场数值模拟,证明了新方法的高效性与实用性,对实现大规模复杂条件下重力异常精细化反演成像与综合解释有重要意义. 展开更多
关键词 空间-波数域 任意傅里叶变换 三维重力异常 CPU-GPU并行
在线阅读 下载PDF
煤矿典型导含水体井下孔中电法三维正演模拟研究
4
作者 张思蔚 《能源与环保》 2024年第3期80-85,共6页
煤矿开采深度逐步加大,隐伏导含水体成为水害的主要隐患。为了提高孔中电法探测精度,采用有限元法通过电位分析对典型导含水体模型进行全空间三维正演模拟,研究了三维正演理论中的电位计算式、边界问题、变分问题及网格剖分,建立了基础... 煤矿开采深度逐步加大,隐伏导含水体成为水害的主要隐患。为了提高孔中电法探测精度,采用有限元法通过电位分析对典型导含水体模型进行全空间三维正演模拟,研究了三维正演理论中的电位计算式、边界问题、变分问题及网格剖分,建立了基础模型,分析了电位分布规律及其误差;针对隐伏正断层、垂直陷落柱、含水采空区典型模型的三维正演模拟结果,总结了各模型的电位等值线特征及变化规律。结果表明,低电阻率使电位等值线变得稀疏,高电阻率使电位等值线变得稠密,且畸变的范围与模型设定的规格相吻合。煤矿井下孔中电法可以有效探测导含水体,探测效率高,靠近靶点,三维正演模拟效果良好,为反演计算提高探测精度奠定了理论基础。 展开更多
关键词 典型导含水体 孔中电法 三维正演模拟 电位分布
在线阅读 下载PDF
复杂场源形态的海洋可控源电磁三维正演 被引量:33
5
作者 韩波 胡祥云 +3 位作者 Adam SCHULTZ 左博新 李建慧 蔡建超 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2015年第3期1059-1071,共13页
在使用电偶极发射源的可控源电磁法(CSEM)勘探中,发射源的方位、长度、形状等对观测数据有重要的影响,然而现有的大部分三维数值模拟方法没有全面地将这些因素考虑进来,很多都只能应对非常简单的场源形态,例如单一方位的点电偶极子,这... 在使用电偶极发射源的可控源电磁法(CSEM)勘探中,发射源的方位、长度、形状等对观测数据有重要的影响,然而现有的大部分三维数值模拟方法没有全面地将这些因素考虑进来,很多都只能应对非常简单的场源形态,例如单一方位的点电偶极子,这有可能显著降低模拟结果的准确性.本文实现了基于交错网格有限体积(FV)离散的海洋CSEM三维正演算法,能够模拟形态相对复杂的场源,包括任意方位的有限长直导线和弯曲导线发射源.该算法使用一次场/二次场方法,只需对二次场使用FV法求解,避免了场源的奇异性问题;一次场的计算为一维正演问题,使用准解析法求解,并且只要在计算一次场时考虑复杂的场源形态便可以实现同样场源的三维正演.通过与一维理论模型的解析解对比验证了三维程序的准确性,并针对三维理论模型进行了一系列正演测试,初步考察了场源形态对三维正演结果的影响. 展开更多
关键词 海洋可控源电磁法 三维正演 有限长线源 场源形态
在线阅读 下载PDF
基于拟态有限体积法的频率域可控源三维正演计算 被引量:21
6
作者 彭荣华 胡祥云 +1 位作者 韩波 蔡建超 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2016年第10期3927-3939,共13页
大规模地球物理电磁数据的定量解释需要发展高效、稳定的三维正反演算法.本文通过求解离散化的三维电场矢量Helmholtz方程,实现了基于有限体积法的频率域可控源电磁(CSEM)三维正演算法.为模拟具有强电性差异的三维电性介质,该算法采用... 大规模地球物理电磁数据的定量解释需要发展高效、稳定的三维正反演算法.本文通过求解离散化的三维电场矢量Helmholtz方程,实现了基于有限体积法的频率域可控源电磁(CSEM)三维正演算法.为模拟具有强电性差异的三维电性介质,该算法采用拟态有限体积法(MFV)对Maxwell方程组进行离散化;另外,为获得稳定、高精度的正演数值结果,采用直接矩阵分解技术来求解离散所得到的大型稀疏线性方程组.对于具有多个发射源的CSEM测量来说,一次矩阵分解结果能够用于同频率下所有场源的正演计算.为降低场源奇异性及边界条件对数值精度的影响,采用虚拟场源校正技术,避免了散射场公式中在构建场源项时所需的大量时间.对于具有多个频率的CSEM的模拟计算,采用分频并行策略来加快三维正演计算.最后,通过与一维层状模型及三维模型的数值结果的对比验证了本文所开发的正演算法对频率域CSEM模拟计算的准确性及有效性,表明该正演算法能够有效应用于三维介质的数值计算.另外,对于多频率CSEM的并行测试结果表明基于分频并行策略的并行计算能够显著地降低正演计算时间. 展开更多
关键词 可控源电磁法 有限体积法 虚拟场源校技术 三维正演 直接分解法
在线阅读 下载PDF
起伏地表频域/时域航空电磁系统三维正演模拟研究 被引量:18
7
作者 张博 殷长春 +1 位作者 刘云鹤 蔡晶 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2016年第4期1506-1520,共15页
由于航空电磁系统具有工作频率低、时间延迟短等特点,地形对航空电磁响应有很大影响,忽略地形影响会给航空电磁数据解释造成很大误差.本文将基于非结构化网格的矢量有限元法应用于模拟起伏地表条件下频域/时域(FD/TD)三维航空电磁系统响... 由于航空电磁系统具有工作频率低、时间延迟短等特点,地形对航空电磁响应有很大影响,忽略地形影响会给航空电磁数据解释造成很大误差.本文将基于非结构化网格的矢量有限元法应用于模拟起伏地表条件下频域/时域(FD/TD)三维航空电磁系统响应.该方法由于采用非结构网格,与传统的结构化网格电磁正演算法相比,能更好地拟合地形和地下不规则异常体,提高对不规则地形和地下介质航空电磁响应的计算精度.通过将计算结果与半空间模型的半解析解及已发表的结果进行对比,检验了本文算法的精度.通过对典型山峰和山谷地形航空电磁响应分析对比,总结了地形对航空电磁响应的影响特征.研究结果对航空电磁地形效应的识别和校正具有指导意义. 展开更多
关键词 航空电磁 三维正演 地形效应 矢量非结构有限元法
在线阅读 下载PDF
隧道衬砌病害探地雷达三维正演模拟及工程验证 被引量:15
8
作者 舒志乐 刘新荣 +1 位作者 刘保县 王道良 《中国铁道科学》 EI CAS CSCD 北大核心 2013年第4期46-53,共8页
从Maxwell旋度方程出发,通过引入复介电常数和复磁导率,推导导电介质UPML吸收边界在考虑和不考虑磁损耗情况下各区的电场和磁场分量的时域差分公式。通过选取各区不同的匹配参数,建立模型区与UPML吸收边界电场和磁场分量的统一离散差分... 从Maxwell旋度方程出发,通过引入复介电常数和复磁导率,推导导电介质UPML吸收边界在考虑和不考虑磁损耗情况下各区的电场和磁场分量的时域差分公式。通过选取各区不同的匹配参数,建立模型区与UPML吸收边界电场和磁场分量的统一离散差分公式。用MATLAB语言编制基于UPML吸收边界的图像可视化的探地雷达三维正演模拟程序。正演模拟衬砌内方形空洞三维地电模型,揭示方形空洞在不同显示剖面上的图谱特征。单剖面的响应特征类似"平锅盖形",z方向切片图定出空洞的形状;反射界面呈强负相反射异常,在立体剖视图中,空洞角点位置产生了2个交叉的双曲线绕射波。将该方法应用于现场隧道衬砌质量的检测,并开孔验证。结果表明:不密实区反射界面多且不规则,同相轴错断、零乱、不连续的较强反射能量团块状异常;在z方向切片图上能清楚反映隧道衬砌密实程度,并且判译直观、准确。 展开更多
关键词 衬砌病害 探地雷达 导电介质 吸收边界 三维正演模拟 隧道
在线阅读 下载PDF
基于二次耦合势的广域电磁法有限体积三维正演计算 被引量:10
9
作者 彭荣华 胡祥云 +2 位作者 李建慧 胡祖志 杨辉 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2018年第10期4160-4170,共11页
与可控源音频大地电磁(CSAMT)相比,广域电磁法通过采用全区视电阻率定义,突破了卡尼亚视电阻率所需的远区条件限制,极大拓展了可控源电磁观测区域和探测深度.考虑到电偶源激发场的三维特征以及地下复杂三维结构,为提高广域电磁数据解释... 与可控源音频大地电磁(CSAMT)相比,广域电磁法通过采用全区视电阻率定义,突破了卡尼亚视电阻率所需的远区条件限制,极大拓展了可控源电磁观测区域和探测深度.考虑到电偶源激发场的三维特征以及地下复杂三维结构,为提高广域电磁数据解释精度,本文实现了基于二次耦合势的广域电磁法三维正演计算.该算法利用Helmholtz定理将麦克斯韦方程转化为库伦规范下的磁矢势和电标势耦合方程,有效改善了离散所得大型线性方程组的谱性质,并通过强加散度条件来消除电场伪解的影响.此外,采用散射场方法,其中一次场使用准解析法求解,二次场使用有限体积法求解,克服了局部激发场源奇异性问题.通过与一维层状模型下电偶源产生的电磁场准解析解对比,验证了本文算法的正确性.在此基础上,利用本文的正演算法对比分析了广域电磁法与CSAMT对典型三维目标体的探测能力,结果表明在相同的观测条件下,广域电磁法能够更准确地反映地下目标体信息,拥有更优的分辨能力. 展开更多
关键词 广域电磁法 可控源音频电磁法 有限体积法 三维正演 二次耦合势
在线阅读 下载PDF
基于电场总场矢量有限元法的接地长导线源三维正演 被引量:23
10
作者 李建慧 Colin G.Farquharson +1 位作者 胡祥云 曾思红 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2016年第4期1521-1534,共14页
电磁场数值模拟的背景场/异常场算法是三维正演的有效策略之一,优点为采用解析法计算电磁场背景场代替场源项、克服了场源奇异性,缺点为不适用于发射源布置于起伏地表或背景模型复杂的情形.总场算法是直接对电磁场总场开展数值模拟,其... 电磁场数值模拟的背景场/异常场算法是三维正演的有效策略之一,优点为采用解析法计算电磁场背景场代替场源项、克服了场源奇异性,缺点为不适用于发射源布置于起伏地表或背景模型复杂的情形.总场算法是直接对电磁场总场开展数值模拟,其难点是有效加载场源、保证近区与过渡区数值解精度.本文以水平电偶源形式分段加载接地长导线源,并以电场总场Helmholtz方程为矢量有限元法控制方程,实现了基于非结构化四面体网格剖分的接地长导线源频率域电磁法三维正演.通过与均匀全空间中水平电偶源产生的电场解析解对比,验证了本文算法的正确性,并分析了四面体外接圆半径与其最短棱边的最大比值和四面体二面角最小值对数值解精度的影响规律.通过与块状高导体地电模型的积分方程法、有限体积法和基于磁矢量势Helmholtz方程的有限元法数值解对比,进一步验证了本文算法正确性,同时说明了非结构化四面体网格能够更加精细地剖分电性异常体,利于获得精确数值解. 展开更多
关键词 接地长导线源 三维正演 矢量有限元法 电场总场 非结构化网格
在线阅读 下载PDF
基于球坐标系下有限差分的地磁测深三维正演 被引量:7
11
作者 李建平 翁爱华 +5 位作者 李世文 李大俊 李斯睿 杨悦 唐裕 张艳辉 《吉林大学学报(地球科学版)》 EI CAS CSCD 北大核心 2018年第2期411-419,共9页
为了计算全球尺度电磁感应的响应,本文介绍地磁测深频率域三维正演。正演算法采用球坐标系下的交错网格有限差分方法,从Maxwell方程的积分形式出发,采用PARDISO对离散后的方程组求解,避免了迭代求解的散度校正。为了验证本文结果的正确... 为了计算全球尺度电磁感应的响应,本文介绍地磁测深频率域三维正演。正演算法采用球坐标系下的交错网格有限差分方法,从Maxwell方程的积分形式出发,采用PARDISO对离散后的方程组求解,避免了迭代求解的散度校正。为了验证本文结果的正确性和精度,与前人的有限元和有限差分方法进行了对比,一维层状模型的三维交错网格有限差分数值结果和解析解相对误差小于5%,双半球模型的计算结果与前人的计算结果完全吻合。三维"棋盘模型"计算表明磁场分量对异常体的大小和位置具有很好的分辨能力。 展开更多
关键词 地磁测深 球坐标 三维正演 有限差分 地磁响应函数
在线阅读 下载PDF
利用积分方程法的大地电磁三维正演 被引量:21
12
作者 徐凯军 李桐林 +1 位作者 张辉 李建平 《西北地震学报》 CSCD 北大核心 2006年第2期104-107,共4页
利用积分方程法实现了均匀导电半空间三维大地电磁响应的数值模拟。求取张量格林函数积分时,采用二次剖分算法解决计算中奇异值问题,对于含有贝塞尔函数的积分项,利用结合连分式展开的高斯求积代替常规的快速汉克尔变换方法,确保了张量... 利用积分方程法实现了均匀导电半空间三维大地电磁响应的数值模拟。求取张量格林函数积分时,采用二次剖分算法解决计算中奇异值问题,对于含有贝塞尔函数的积分项,利用结合连分式展开的高斯求积代替常规的快速汉克尔变换方法,确保了张量格林函数的正确计算并提高了计算精度。最后通过数值模拟结果的对比及模型试算验证了算法的正确性,所实现的三维大地电磁数值模拟算法为理论研究三维地电构造的大地电磁响应的分布规律提供了有效的工具,也为研究三维反演算法奠定了基础。 展开更多
关键词 积分方程 大地电磁 三维正演 张量格林函数
在线阅读 下载PDF
基于Maxwell 3D瞬变电磁法三维正演研究 被引量:11
13
作者 蒋大青 付志红 +1 位作者 侯兴哲 张淮清 《电测与仪表》 北大核心 2012年第6期29-32,共4页
瞬变电磁法(TEM)是一种时间域电磁法,在地球物理勘探领域得到广泛应用。目前三维反演的快速性与精确度均存在着很多问题,这就需要三维正演分析提供支持,但是应用于瞬变电磁法的专业三维正演软件非常缺乏。基于此,本文创新性地将Ansoft M... 瞬变电磁法(TEM)是一种时间域电磁法,在地球物理勘探领域得到广泛应用。目前三维反演的快速性与精确度均存在着很多问题,这就需要三维正演分析提供支持,但是应用于瞬变电磁法的专业三维正演软件非常缺乏。基于此,本文创新性地将Ansoft Maxwell 3D引入瞬变电磁法地质探测的正演分析中,建立了均匀半空间瞬变电磁正演模型,与传统的电磁模型分析软件EMMA进行对比,证实了Maxwell用于瞬变电磁法正演的可行性与有效性;同时,还构建了低阻三维地质异常体的Maxwell正演模型。仿真结果可以很好地反映地下异常体的存在,证明Ansoft Maxwell 3D能很好的应用于TEM正演模型的构建与分析。 展开更多
关键词 瞬变电磁法 三维正演 ANSOFT MAXWELL 二次场
在线阅读 下载PDF
时域多分辨率法在探地雷达三维正演模拟中的应用 被引量:10
14
作者 冯德山 戴前伟 瓮晶波 《中南大学学报(自然科学版)》 EI CAS CSCD 北大核心 2007年第5期975-980,共6页
应用小波伽略金方法,对Maxwell的2个旋度方程进行离散化,导出DB2-MRTD算法的探地雷达3D差分公式、数值稳定性条件。在此基础上,开发探地雷达MRTD法正演模拟程序。所开发的对三维立方体小空洞体模型进行正演模拟。研究结果表明:采用该程... 应用小波伽略金方法,对Maxwell的2个旋度方程进行离散化,导出DB2-MRTD算法的探地雷达3D差分公式、数值稳定性条件。在此基础上,开发探地雷达MRTD法正演模拟程序。所开发的对三维立方体小空洞体模型进行正演模拟。研究结果表明:采用该程序极大地提高了运算速度,改善了三维探地雷达正演方法;相应的正演合成三维剖视图及切片图指示了立方体小空洞体的雷达波空间的三维响应特征;模拟结果表明,采用时域多分辨率法可提高探地雷达探测的可靠性、准确度。 展开更多
关键词 探地雷达 时域多分辨法 三维正演模拟 DB2小波基
在线阅读 下载PDF
起伏地形条件下长偏移距瞬变电磁三维正演 被引量:5
15
作者 王新宇 严良俊 +3 位作者 毛玉蓉 黄鑫 谢兴兵 周磊 《吉林大学学报(地球科学版)》 CAS CSCD 北大核心 2022年第3期754-765,共12页
地球物理勘探区的地形往往是起伏复杂的,传统长偏移距瞬变电磁法(long-offset transient electromagnetic method, LOTEM)三维正演模拟一般不考虑地形的影响,而地形因素会严重影响LOTEM数据的合理解释。因此,本文开展起伏地形条件下非... 地球物理勘探区的地形往往是起伏复杂的,传统长偏移距瞬变电磁法(long-offset transient electromagnetic method, LOTEM)三维正演模拟一般不考虑地形的影响,而地形因素会严重影响LOTEM数据的合理解释。因此,本文开展起伏地形条件下非结构矢量有限元LOTEM三维正演模拟。首先,应用三维建模软件Blender准确建立复杂地电模型,采用非结构四面体网格精细刻画起伏地形。其次,基于偶极子离散的长导线源近似技术,将长导线源离散为若干段等效电偶极子,并引入二阶后退欧拉格式完成时间离散。进而实现LOTEM快速、高精度三维正演。在保证计算精度的前提下,模拟了起伏地形条件下复杂地质模型的瞬变电磁场响应。模拟结果表明,由于地形的影响,LOTEM电磁场响应“畸变”严重,会削弱甚至完全“埋没”地下异常体的瞬变电磁场响应,严重影响地下目标体的有效探测。 展开更多
关键词 长偏移距瞬变电磁法 起伏地形 非结构矢量有限元 Blender 二阶后退欧拉格式 三维正演
在线阅读 下载PDF
瞬变电磁Crank-Nicolson FDTD三维正演 被引量:14
16
作者 孙怀凤 柳尚斌 杨洋 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2021年第1期343-354,共12页
时域有限差分(FDTD)方法使用Yee网格剖分电磁场的空间采样,通过时间步迭代实现电磁场数值模拟,具有内存消耗低、计算简单等特点,常用于瞬变电磁三维正演.然而,常规FDTD方法的时间迭代步长Δt受Courant-Friedrich-Lewy(CFL)条件严格限制... 时域有限差分(FDTD)方法使用Yee网格剖分电磁场的空间采样,通过时间步迭代实现电磁场数值模拟,具有内存消耗低、计算简单等特点,常用于瞬变电磁三维正演.然而,常规FDTD方法的时间迭代步长Δt受Courant-Friedrich-Lewy(CFL)条件严格限制,过多的迭代次数以及过密的采样往往导致计算速度慢、累积误差不断增大.本文提出一种不受CFL条件约束的无条件稳定隐式差分算法Crank-Nicolson FDTD(CN-FDTD)用于瞬变电磁三维正演.基于Crank-Nicolson差分方法对Maxwell方程组重新离散,空间网格仍然采用Yee元胞,时间步进采用在整时间步电场、磁场同时采样的策略,建立无条件稳定FDTD格式,突破CFL条件限制.与常规FDTD交替采样相比,CN-FDTD电场、磁场同时采样的策略构成的隐式差分格式,需要求解大型稀疏矩阵方程组.通常,瞬变电磁三维正演模型中产生的矩阵阶数往往较大,需要占用大量内存和求解时间.为解决上述问题,采用Crank-Nicolson-cycle-sweep-uniform(CNCSU-FDTD)方法近似求解CN-FDTD方程,在保证求解精度的同时,计算效率大幅提高.在边界条件处理上,采用双线性变换推导了复频率参数完全匹配层(CFS-PML)吸收边界.采用均匀半空间模型、四类三层模型进行精度验证,发现CN-FDTD三维正演结果与解析解、线性数字滤波解吻合较好.之后,与接触带上的低阻复杂模型进行对比,结果显示CN-FDTD正演结果与矢量有限元、有限体积法以及FDTD计算结果吻合较好.在此基础上,研究了时间步放大对CN-FDTD计算精度的影响,发现最大时间步放大到常规FDTD的3200倍时才会在晚期出现较明显的误差.在一台CPU为Intel Core i5-7300HQ的笔记本电脑单线程计算条件下,模拟到关断后30 ms仅需要50 min.在进行并行化后,将有望实现复杂模型分钟级的三维正演,从而为三维反演提供可靠、快速的正演方法. 展开更多
关键词 瞬变电磁 三维正演 CN-FDTD 双线性变换 CFS-PML
在线阅读 下载PDF
频率域海底电磁法三维正演研究 被引量:3
17
作者 傅磊 刘四新 +3 位作者 吴俊军 王飞 张丽丽 李金泉 《吉林大学学报(地球科学版)》 EI CAS CSCD 北大核心 2010年第S1期67-69,86,共4页
为了研究频率域海底电磁法在油气勘探中的可行性,根据海底电磁法的观测特点,建立海底三维地电模型,采用频率域FDTD三维正演模拟程序,计算不同参数时频率域电磁响应,获得了不同模型、不同频率下电磁场的响应。
关键词 海洋油气勘探 频率域电磁法 FDTD 三维正演
在线阅读 下载PDF
基于混合网格有限元的直流电阻率法三维正演研究 被引量:4
18
作者 王新宇 王程 +3 位作者 毛玉蓉 严良俊 周磊 高文龙 《煤田地质与勘探》 CAS CSCD 北大核心 2022年第5期136-143,共8页
直流电阻率法因效率高、成本低等特点广泛应用于有色金属、煤田等矿产资源勘探和地质调查行业。但在井–地、地–井电阻率法勘探过程中,钻孔因素(井液电阻率和钻孔孔径)对不同探测模式下视电阻率响应的影响尚未明确,是否影响视电阻率资... 直流电阻率法因效率高、成本低等特点广泛应用于有色金属、煤田等矿产资源勘探和地质调查行业。但在井–地、地–井电阻率法勘探过程中,钻孔因素(井液电阻率和钻孔孔径)对不同探测模式下视电阻率响应的影响尚未明确,是否影响视电阻率资料的解释是值得探讨的一个问题。因此,提出采用混合网格有限元法实现直流电阻率三维正演。给出异常电位法满足的边值问题及有限元变分问题,应用三棱柱和四面体混合网格实现对计算区域的快速离散,并建立2种网格的线性插值基函数和单元系数矩阵;采用SSOR–PCG迭代算法求解异常电位满足的大型线性方程组,得到各观测点的响应。在保证计算精度的前提下,应用混合网格有效地离散钻孔地电模型,探讨钻孔因素对井–地、地–井观测方式视电阻率数据的影响特征。对于井–地观测方式:钻孔附近视电阻率受钻孔因素影响最大,严重影响视电阻率资料的合理解释,但随着测点收发距的增加,视电阻率响应逐渐趋于围岩电阻率;而且随着发射源深度的增加,钻孔对地表视电阻率响应的影响也逐渐减小。对于地–井观测方式:钻孔因素对浅部数据影响较大,对深部数据影响小,相比于井液电阻率,视电阻率响应更易受钻孔孔径的影响。直流电阻率法混合网格有限元三维正演算法的提出,对实际井–地、地–井电阻率法勘探具有理论指导意义,可结合钻孔信息进行正演模拟分析,选择合适的收发距,有效压制钻孔因素对实测资料的影响。 展开更多
关键词 直流电阻率法 井液电阻率 钻孔孔径 混合网格有限元法 三棱柱 四面体 三维正演
在线阅读 下载PDF
探地雷达时域多分辨法(MRTD)三维正演模拟 被引量:8
19
作者 冯德山 戴前伟 《地球物理学进展》 CSCD 北大核心 2008年第5期1621-1625,共5页
应用小波伽略金方法,对Maxwell方程进行离散化,导出了DB2-MRTD算法的探地雷达3D差分公式、数值稳定性条件.在此基础上,开发了探地雷达MRTD(multi-resolution time domain)法正演模拟程序,该程序极大地提高了运算速度,改善了三维探地雷... 应用小波伽略金方法,对Maxwell方程进行离散化,导出了DB2-MRTD算法的探地雷达3D差分公式、数值稳定性条件.在此基础上,开发了探地雷达MRTD(multi-resolution time domain)法正演模拟程序,该程序极大地提高了运算速度,改善了三维探地雷达正演方法,并利用该自制程序,对三角形金属体模型进行了正演模拟,得到了其相应的正演合成三维剖视图及切片图,通过对这些模拟结果进行分析,可以加深对三维雷达反射特征的认识,提高探地雷达探测的可靠性、准确度,同时也说明时域多分辨率法在探地雷达三维正演模拟中的有效性. 展开更多
关键词 探地雷达 时域多分辨法 三维正演模拟 DB2小波基
在线阅读 下载PDF
基于重启多项式Krylov子空间模型降阶的全波形瞬变电磁三维正演 被引量:2
20
作者 周建美 景旭 +2 位作者 鲁凯亮 刘文韬 李貅 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2023年第6期2646-2657,共12页
瞬变电磁三维正演在内存占用和求解时间上都是大的挑战.瞬变电磁响应受到发射波形的影响,正演计算中还需考虑发射波形.基于重启多项式Krylov子空间模型降阶方法,本文实现了一种新的全波形瞬变电磁三维正演算法.采用规则六面体网格的拟... 瞬变电磁三维正演在内存占用和求解时间上都是大的挑战.瞬变电磁响应受到发射波形的影响,正演计算中还需考虑发射波形.基于重启多项式Krylov子空间模型降阶方法,本文实现了一种新的全波形瞬变电磁三维正演算法.采用规则六面体网格的拟态有限体积方法对全波形控制方程进行空间离散,引入一阶φ函数推导得到全时段响应的统一表达式,利用重启多项式Krylov子空间算法求解φ函数表达式,对于给定的重启子空间维度,基于残差公式得到满足给定精度的任意时刻的正演响应.该算法内存占用小,不需要求解大型线性方程组,内存占用主要由空间离散网格数量和重启子空间维度决定.数值模型正演结果验证了本文算法的有效性. 展开更多
关键词 瞬变电磁 三维正演 发射波形 模型降阶
在线阅读 下载PDF
上一页 1 2 10 下一页 到第
使用帮助 返回顶部