任意起伏地形下重力异常三维正演及并行计算 被引量：1

1

作者 戴世坤 朱德祥 +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 职称材料

探地雷达频率域2.5维正演 被引量：5

2

作者 戴世坤 欧阳振崇 +5 位作者 周印明 张钱江 李昆 赵东东 陈轻蕊 凌嘉宣 《电子与信息学报》 EI CSCD 北大核心 2021年第1期145-153,共9页

该文从频率域电磁法满足的控制方程出发,采用有限单元法实现了频率域2.5维探地雷达(GPR)正演模拟。重点分析了波数域电磁场谱随相对介电常数和收发距变化的规律,探讨了2.5维GPR正演模拟的波数选取问题;基于Open MP并行算法与串行算法的... 该文从频率域电磁法满足的控制方程出发,采用有限单元法实现了频率域2.5维探地雷达(GPR)正演模拟。重点分析了波数域电磁场谱随相对介电常数和收发距变化的规律,探讨了2.5维GPR正演模拟的波数选取问题;基于Open MP并行算法与串行算法的计算效率对比,表明频率域2.5维GPR数值模拟方法具有高效率、高精度和高度并行性的特点,为雷达正演提供重要理论参考依据和技术支撑,是GPR全波形反演的重要基础。 展开更多

关键词 探地雷达 正演 2.5维 频率域 波数选取 Open MP 并行

在线阅读 下载PDF 职称材料

基于矢量位和标量位的空间波数混合域电磁三维正演模拟 被引量：6

3

作者 戴世坤 赵东东 +3 位作者 李昆 张钱江 凌嘉宣 陈轻蕊 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2022年第1期404-416,共13页

高效、高精度电磁三维数值模拟是制约大规模电磁数据精细化三维反演成像、人机交互定量解释的核心问题.针对一问题,本文提出一种基于矢量位和标量位的空间波数混合域电磁场三维数值模拟方法.该方法利用沿水平方向的二维傅里叶变换将空... 高效、高精度电磁三维数值模拟是制约大规模电磁数据精细化三维反演成像、人机交互定量解释的核心问题.针对一问题,本文提出一种基于矢量位和标量位的空间波数混合域电磁场三维数值模拟方法.该方法利用沿水平方向的二维傅里叶变换将空间域矢量位和标量位耦合偏微分方程组转换为波数之间相互独立的常微分耦合方程组,将一个大规模三维问题分解为多个一维小问题,具有高度并行性,由此大大减少了计算量和存储量;保留垂向为空间域,浅层网格剖分适当加密,深层网格剖分适当稀疏,有效兼顾了计算精度与计算效率;采用有限单元法求解不同波数的常微分方程,充分利用追赶法求解定带宽线性方程组的高效性进一步提高数值模拟效率.在模型算例中,设计棱柱体模型验证了本文方法的正确性、计算精度和计算效率.数值试验结果表明本文方法具有数值精度高、并行度高、占用内存小、计算效率高的特性,比传统有限单元法三维数值模拟方法计算效率高1~2个数量级,且网格剖分规模越大,该方法计算效率优势越明显. 展开更多

关键词 空间波数混合域 电磁 三维 数值模拟

在线阅读 下载PDF 职称材料

层状单轴各向异性介质中长导线源电磁响应的快速计算 被引量：7

4

作者 欧阳芳 戴世坤 +2 位作者 张钱江 陈龙伟 陈轻蕊 《石油地球物理勘探》 EI CSCD 北大核心 2016年第5期1012-1020,839,共9页

在场源附近,长导线源不能视作电偶极子源。长导线源电磁响应的常用计算方法是对场源进行离散,再对离散电偶极子电磁场进行叠加,这种方法计算简单,但计算效率较低。针对这一问题,提出了有限单元形函数积分算法。该算法从电偶极子场的解... 在场源附近,长导线源不能视作电偶极子源。长导线源电磁响应的常用计算方法是对场源进行离散,再对离散电偶极子电磁场进行叠加,这种方法计算简单,但计算效率较低。针对这一问题,提出了有限单元形函数积分算法。该算法从电偶极子场的解析表达式出发,将长导线源进行剖分,通过在每个剖分单元上构造形函数来求解积分表达式。通过模型算例对传统电偶极子叠加算法和形函数积分算法的计算精度及计算效率进行了测试。从测试结果可以看出,无论在场源中远区还是近区,形函数积分算法只需要将长导线源剖分成两个单元就能达到很高精度,而电偶极子叠加算法需要将长导线源离散成大量的电偶极子才能达到相同数值精度,相同数值精度条件下形函数积分算法具有更高的计算效率,且离场源越近计算效率优势越明显。 展开更多

关键词 层状模型 单轴各向异性介质 形函数积分 长导线源

在线阅读 下载PDF 职称材料

重力异常场空间-波数混合域三维数值模拟 被引量：6

5

作者 戴世坤 陈轻蕊 +1 位作者 李昆 凌嘉宣 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2020年第5期2107-2119,共13页

重力勘探中复杂条件下的三维正演计算量大存储要求高,使得这种条件下重力勘探高效、精细正反演变得困难.针对这一问题,提出一种空间-波数混合域数值模拟方法,该方法将空间域引力位积分进行水平方向二维傅里叶变换,将三维空间域卷积问题... 重力勘探中复杂条件下的三维正演计算量大存储要求高,使得这种条件下重力勘探高效、精细正反演变得困难.针对这一问题,提出一种空间-波数混合域数值模拟方法,该方法将空间域引力位积分进行水平方向二维傅里叶变换,将三维空间域卷积问题转换为多个不同波数之间相互独立的空间垂向一维积分问题,一维积分垂向可离散为多个单元积分之和,每个单元采用二次形函数表征密度变化,可得出单元积分的解析表达式.该方法计算量和存储需求少,算法高度并行;保留垂向为空间域,优势之一在于可根据实际情况合理调整单元疏密程度,准确模拟任意复杂地形和密度异常体的重力异常,兼顾计算精度与计算效率;优势之二在于用形函数拟合求得积分的解析解,计算精度和效率高;充分利用一维形函数积分的高效和高精度,不同波数之间一维积分高度并行性及快速傅里叶变换的高效性,实现重力异常场三维数值模拟.设计棱柱体模型,通过数值解和解析解对比验证了该方法的正确性、适用性和高效性.针对任意复杂地形条件下的重力场及其张量的模拟问题,提出一种快速算法,对其有效性进行了验证.探究标准FFT法的截断效应对计算精度的影响,对比分析Gauss-FFT法和标准FFT扩边法两种方法的计算精度和效率,总结了二者的选取策略,结果表明选用标准FFT扩边法计算效率更高.实际地形的数值模拟表明本文算法适用于任意复杂地形的高效计算. 展开更多

关键词 重力数值模拟 空间-波数混合域 傅里叶变换 形函数积分

在线阅读 下载PDF 职称材料

各向异性介质空间—波数混合域直流电阻率法三维数值模拟研究 被引量：5

6

作者 戴世坤 凌嘉宣 +2 位作者 陈轻蕊 张莹 李昆 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2022年第7期2729-2740,共12页

将各向异性介质分成各向同性背景介质和各向异性异常介质,并提出一种空间-波数混合域方法实现了各向异性介质下直流电阻率法的三维数值模拟.不同于传统直流电阻率数值模拟方法,本文算法直接对空间域异常电位满足的偏微分方程沿水平方向... 将各向异性介质分成各向同性背景介质和各向异性异常介质,并提出一种空间-波数混合域方法实现了各向异性介质下直流电阻率法的三维数值模拟.不同于传统直流电阻率数值模拟方法,本文算法直接对空间域异常电位满足的偏微分方程沿水平方向进行二维傅里叶变换,使水平方向转换成波数域,保留垂直方向为空间域,可根据地下介质电流密度变化的快慢灵活剖分.这样可把空间域异常电位满足的三维偏微分方程转化成不同波数满足的一维常微分方程,把一个大规模三维数值模拟问题分解为多个一维数值模拟问题,利用一维有限单元法求解方程组,并通过采用压缩算子迭代计算,最终获得较为精确的数值解.与自适应有限单元法对比验证了本文算法的正确性;测试了算法的收敛性,结果表明在满足精度要求的情况下,算法的收敛性只与异常体和围岩之间的电导率差异相关,而与异常体大小和埋深无关;分析了算法计算效率,结果表明算法的计算效率与剖分网格节点成线性关系,算法可在微型计算机中较快计算出剖分节点总数超过千万的各向异性模型的结果;设计简单观测系统并验证其具有反映地下各向异性结构特性的能力;最后模拟异常体沿着不同方向旋转不同角度时的响应特征,对比分析可知异常体为各向异性情况下,沿不同方向旋转相同角度或沿同一个方向旋转不同角度,在同一个位置上测量得到的视电阻率有差异. 展开更多

关键词 空间波数混合域 二维正反傅里叶变换 有限单元法 压缩算子

在线阅读 下载PDF 职称材料

频率和收发距广普适用的均匀层状介质电磁场直接积分算法 被引量：2

7

作者 戴世坤 曾铃 +3 位作者 周印明 李昆 陈轻蕊 凌嘉宣 《石油地球物理勘探》 EI CSCD 北大核心 2020年第6期1364-1372,1394,1167,共11页

在地球物理勘探中,均匀层状介质电磁场解析表达式含有汉克尔积分,其核函数为0阶和1阶Bessel函数,随着其宗量的增大,Bessel函数呈现出快速振荡和慢衰减的特性,由此造成汉克尔积分难以高效、高精度计算,特别是高频和大收发距情况下难度更... 在地球物理勘探中,均匀层状介质电磁场解析表达式含有汉克尔积分,其核函数为0阶和1阶Bessel函数,随着其宗量的增大,Bessel函数呈现出快速振荡和慢衰减的特性,由此造成汉克尔积分难以高效、高精度计算,特别是高频和大收发距情况下难度更大。针对这个问题,提出一种高效、高精度的直接积分方法(Direct integration method),其基本思想是Bessel函数可以在两个区间并分别用不同的多项式展开,每一区间的汉克尔积分可离散为多个单元积分之和,每个单元被积函数可采用三次样条插值函数表示,由此可求得单元积分的解析解,通过叠加,从而求得汉克尔积分的数值解。在此基础上,利用均匀全空间电偶极子电磁场的解析解,正确选取积分范围并合理剖分积分单元。数值解与解析解的对比表明该算法正确可靠。与数字滤波类算法的计算结果对比表明,该算法广泛适用于不同频率和不同收发距条件下电磁场的计算,具有较强的普适性。 展开更多

关键词 汉克尔积分 直接积分法 高频电磁场 层状介质

在线阅读 下载PDF 职称材料

空间—波数域二维强磁场及其梯度张量数值模拟研究 被引量：2

8

作者 戴世坤 冉应强 +3 位作者 张莹 陈轻蕊 凌嘉宣 贾金荣 《石油地球物理勘探》 EI CSCD 北大核心 2023年第1期228-237,共10页

基于磁位泊松方程提出了一种空间—波数混合域二维数值模拟方法。该方法将二维强磁性体磁位满足的偏微分方程通过一维傅里叶变换转化为不同波数之间相互独立的常微分方程,即将二维数值模拟问题转换为一维数值模拟问题,并利用一维有限元... 基于磁位泊松方程提出了一种空间—波数混合域二维数值模拟方法。该方法将二维强磁性体磁位满足的偏微分方程通过一维傅里叶变换转化为不同波数之间相互独立的常微分方程,即将二维数值模拟问题转换为一维数值模拟问题,并利用一维有限元法进行求解,可大大降低存储需求和计算量。对求得的强磁场利用紧算子进行迭代以获得高精度近似解。算法中保留垂向为空间域,可根据需求对网格进行灵活剖分,有利于复杂条件下的模拟。另外,不同波数之间常微分方程相互独立,具有较好的并行性。算法充分利用了傅里叶变换的快速性和迭代算法的稳定性,实现了强磁场及其梯度张量的二维数值模拟。设计三个模型(水平地形下的圆柱体模型、棱柱体模型及起伏地形模型)进行计算,对文中方法数值解与解析解、COMSOL软件所得的数值解进行对比,结果表明该算法计算速度快、精度高,适用于起伏地形条件下的数值计算。 展开更多

关键词 强磁场 梯度张量 空间—波数混合域 二度体正演

在线阅读 下载PDF 职称材料

空间-波数域三维大地电磁场积分方程法数值模拟 被引量：5

9

作者 戴世坤 陈轻蕊 +1 位作者 凌嘉宣 李昆 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2022年第6期2294-2310,共17页

地球物理勘探中,地下复杂地质体电磁场三维数值模拟计算量大、存储要求高.针对这一问题,本文提出一种空间-波数域三维电磁场数值模拟方法,该方法利用电磁场积分方程为卷积的特点,将积分沿水平方向进行二维傅里叶变换,将电磁场三维空间... 地球物理勘探中,地下复杂地质体电磁场三维数值模拟计算量大、存储要求高.针对这一问题,本文提出一种空间-波数域三维电磁场数值模拟方法,该方法利用电磁场积分方程为卷积的特点,将积分沿水平方向进行二维傅里叶变换,将电磁场三维空间域卷积问题转换为不同波数相对独立的多个一维积分问题,由此计算量和存储需求大大减少,易于实现并行计算.采用有限单元法中的形函数进行一维数值积分,一维积分离散为多个单元积分之和,每个单元采用二次形函数表征电流变化,可得出单元积分的解析表达式;保留垂向为空间域,优势之一在于可根据实际情况合理调整单元疏密程度,准确模拟任意复杂地形和导电率异常,兼顾计算精度与计算效率,优势之二是用形函数拟合求得积分的解析解,计算精度和效率高;最后引入压缩算子,采用迭代求得电磁场的数值解.本文方法充分利用不同波数之间一维积分高度并行性、一维形函数积分的高精度及快速傅里叶变换的高效性,实现电磁场三维高效高精度数值模拟.设计模型将本文方法数值模拟结果和软件INTEM3D的数值模拟结果对比,验证了方法的正确性;设计高阻和低阻异常体,研究了异常电导率与背景电导率差异对迭代收敛速度的影响;改变计算规模,随着计算网格增多,算法耗时与存储呈近似线性增长;设计复杂模型与目前主流数值模拟方法对比,本文算法速度快1个数量级以上,且计算规模越大,算法效率的优势越明显.研究结果表明,本文提出的空间-波数域三维电磁场数值模拟方法理论和方法正确,计算效率高,对计算机存储要求低,算法高度并行,适合任意复杂条件大规模三维电磁场高效、高精度数值模拟. 展开更多

关键词 三维电磁场 空间-波数域 数值模拟 形函数积分 压缩算子

在线阅读 下载PDF 职称材料

空间波数混合域三维各向异性强磁性介质磁场数值模拟 被引量：1

10

作者 戴世坤 张莹 +3 位作者 凌嘉宣 陈轻蕊 贾金荣 冉应强 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2023年第2期825-840,共16页

各向异性普遍存在于强磁性体矿物中,为研究各向异性强磁性体磁场响应特征,本文提出一种空间波数混合域三维各向异性磁场数值模拟方法.该方法首先将各向异性强磁性体磁位满足的三维偏微分方程进行水平方向二维傅里叶变换,将其降为不同波... 各向异性普遍存在于强磁性体矿物中,为研究各向异性强磁性体磁场响应特征,本文提出一种空间波数混合域三维各向异性磁场数值模拟方法.该方法首先将各向异性强磁性体磁位满足的三维偏微分方程进行水平方向二维傅里叶变换,将其降为不同波数之间相互独立的一维常微分方程;然后加载准确的上下边界条件,采用二次插值有限单元法计算一维常微分方程,得到五对角方程,采用追赶法进行高效求解;最后采用迭代法求解场分量,引入紧算子保证迭代稳定收敛;综合傅里叶变换的高效性、一维方程求解的快速性和迭代算法的稳定性,实现各向异性强磁性体磁场的三维高效、高精度数值模拟.设计各向异性椭球模型验证算法的正确性,并分析紧算子对不同各向异性磁化率模型的迭代收敛性;与COMSOL Multiphysics软件对比计算效率,表明相同节点下本文算法效率优于常规三维有限元方法,且计算节点总数越多优势越明显.重点研究各向异性参数改变对VTI、HTI、TTI强磁介质异常场响应的影响.最后采用某磁铁矿DEM高程数据模拟起伏地形对各向异性强磁性体磁异常场幅值和形态的影响,体现出本文算法对各向异性强磁性体大规模复杂地形的适应性. 展开更多

关键词 强磁化率 磁化率各向异性 空间波数混合域 紧算子 起伏地形

在线阅读 下载PDF 职称材料

基于Lorenz规范的空间-波数域大地电磁三维正演

11

作者 戴世坤 陈轻蕊 +4 位作者 凌嘉宣 李昆 赵东东 张莹 张钱江 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2022年第12期4963-4980,共18页

大地电磁勘探方法由于其成本低、施工简单、探测深度广等优点,广泛应用于矿产资源普查、油气勘探和深部构造研究等领域.如何提高大规模三维大地电磁数值模拟的精度和效率一直是研究热点.本文基于空间-波数域方法,实现了基于Lorenz规范... 大地电磁勘探方法由于其成本低、施工简单、探测深度广等优点,广泛应用于矿产资源普查、油气勘探和深部构造研究等领域.如何提高大规模三维大地电磁数值模拟的精度和效率一直是研究热点.本文基于空间-波数域方法,实现了基于Lorenz规范的空间-波数域三维大地电磁数值模拟.基于二次场计算原理,引入Lorenz规范,将Maxwell方程组转化为关于二次场矢量位的亥姆霍兹方程;利用水平方向二维傅里叶变换,将空间域三维偏微分方程转换为多个波数下相互独立的常微分方程,方程采用二次插值有限单元法计算,得到定带宽线性方程组,方程计算量小、并行性好,采用追赶法求解,提高了算法效率;引入压缩算子,用迭代法逐次逼近真实解.充分利用了空间-波数域方法数值精度高、内存需求少、效率高的特点.设计棱柱体模型验证了算法的正确性、分析了算法的收敛性,说明算法对不同频率、不同电导率对比度模型均具有很好的适应性.利用Dublin(DTM1)模型进行三维大地电磁数值模拟,结果表明:在满足精度要求的前提下,空间-波数域算法比空间域算法占用内存少、耗时短;相比基于Coulomb规范的空间-波数域算法,基于Lorenz规范的空间-波数域方法耗时更短、占用内存更少,效率提高至少3倍以上,体现了新方法的优势. 展开更多

关键词 空间-波数域 Lorenz规范 三维大地电磁 数值模拟

在线阅读 下载PDF 职称材料

空间波数混合域磁异常场积分解三维数值模拟 被引量：3

12

作者 李昆 戴世坤 +4 位作者 陈轻蕊 张钱江 赵东东 王顺国 凌嘉宣 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2019年第11期4437-4450,共14页

本文提出一种空间波数混合域磁异常场三维数值模拟方法.该方法利用磁位三维空间域积分为卷积的特点,沿水平方向进行二维傅里叶变换,把空间域磁位满足的三维积分问题转化为不同波数之间相互独立的垂向一维积分问题.保留垂向为空间域,优... 本文提出一种空间波数混合域磁异常场三维数值模拟方法.该方法利用磁位三维空间域积分为卷积的特点,沿水平方向进行二维傅里叶变换,把空间域磁位满足的三维积分问题转化为不同波数之间相互独立的垂向一维积分问题.保留垂向为空间域,优势之一在于便于浅层单元剖分可适当加密,随着深度增加,单元剖分适当稀疏,可以准确模拟任意复杂地形和磁性体的磁异常,兼顾了计算精度与计算效率;优势之二在于一维积分垂向可离散为多个单元积分之和,每个单元采用二次形函数表征磁化强度,可得出单元积分的解析表达式,计算精度高、效率高.该方法充分利用一维形函数积分的高效和高精度、快速傅里叶变换的高效性及算法高度并行性,实现了磁异常场高效、高精度的数值模拟.设计棱柱体模型,将模型解析解与空间波数混合域法的数值解对比,结果表明该方法计算精度高、效率高.设计了组合棱柱体复杂模型,对比分析了标准FFT扩边法与Gauss-FFT法的计算精度与计算效率,总结了标准FFT的扩边系数选取策略.针对任意复杂地形条件下的磁异常模拟问题,本文提出一种适用于起伏地形条件下的磁异常场快速计算方法,并对其有效性进行了验证. 展开更多

关键词 磁异常场三维数值模拟 空间波数混合域 傅里叶变换 形函数法

在线阅读 下载PDF 职称材料

基于三维数值模拟的激电法差分装置响应特征分析

13

作者 凌嘉宣 陈志芳 +3 位作者 邓维 戴世坤 周印明 陈轻蕊 《石油地球物理勘探》 EI CSCD 北大核心 2023年第5期1244-1254,共11页

为了深入分析激电法差分装置的响应特征,提出一种高精度、高效率的三维数值模拟方法。该数值模拟方法通过二维傅里叶变换将三维激电数值模拟问题转化为波数域一维数值模拟问题,使用一维有限单元法求解,同时引入算子进行迭代,最终获得高... 为了深入分析激电法差分装置的响应特征,提出一种高精度、高效率的三维数值模拟方法。该数值模拟方法通过二维傅里叶变换将三维激电数值模拟问题转化为波数域一维数值模拟问题,使用一维有限单元法求解,同时引入算子进行迭代,最终获得高精度数值解。该算法可避免进行大型稀疏线性方程组的直接求解,大幅度降低了计算量和存储需求,提高了计算效率。以提出的三维数值模拟方法为工具,分析激电法差分装置响应特征。基于高阻/低阻的高极化异常体模型测试结果表明,单孔井—井测量下差分装置对井旁盲矿的分辨效果优于其他装置,且低阻体较高阻体对电流具有更强的敏感度;基于不同背景介质下的低阻高极化异常体模型的测试结果表明,井—地测量下采用差分装置得到的视电阻率和极化率比采用二极装置得到的视电阻率和极化率能更准确地反映背景介质和异常体在垂向上的空间分布特征。因此,激电法差分装置勘探效果优于其他方法装置,研究结论可为野外勘探测量装置的选择提供参考。 展开更多

关键词 激发极化法 二维傅里叶变换 三维数值模拟 差分装置 电磁响应特征

在线阅读 下载PDF 职称材料