田红军
- 作品数:2 被引量:10H指数:1
- 供职机构:中南大学地球科学与信息物理学院更多>>
- 发文基金:湖南省自然科学基金国家自然科学基金更多>>
- 相关领域:天文地球更多>>
- 带地形的二维各向异性大地电磁测深数值模拟被引量:10
- 2015年
- 以往主要在水平层状地表研究大地电磁二维各向异性系数对观测值的影响。通过有限元数值模拟方法,讨论同时存在起伏地形以及各向异性对二维大地电磁测深法中TM极化模式的响应。首先给出二维各向异性介质中TM极化模式的变分问题;其次用规则矩形单元对研究区域进行剖分,用双线性函插值对变分问题进行求解;最后在地垒、地堑地形下,改变各向异性系数模拟TM极化模式中视电阻率的变化。结果表明:地垒地形对TM极化模式视电阻率曲线表现为低阻异常,地堑地形对TM极化视电阻率曲线呈现高阻异常,两种地形对观测数据的影响要远大于各向异性系数对观测值的影响。
- 田红军佟铁钢
- 关键词:大地电磁测深有限单元法各向异性
- 任意起伏地形下重力异常三维正演及并行计算
- 2024年
- 为了进一步提高空间-波数域三维重力异常正演算法的适用范围和计算效率,本文采用任意傅里叶变换算法实现了空间-波数域三维重力异常正演,且在NVIDIA CUDA平台上进行CPU-GPU并行加速.任意傅里叶变换算法的基本思想是将二维傅里叶变换转化为两个一维傅里叶变换,一维傅里叶变换积分离散为多个单元积分累加和,离散单元中原函数采用二次插值形函数拟合,求出单元积分的解析表达式.相比现有的傅里叶变换算法,新方法具有采样灵活、积分精度高、计算速度快和傅里叶变换的截断效应小等优势.利用空间-波数域算法的高度并行性,采用CPU并行求解常微分方程,GPU并行计算任意傅里叶变换,实现了CPU-GPU并行加速方案,进一步提升了本文算法效率.利用常密度模型,对比数值解和解析解,结果表明本文算法正确;利用变密度模型对比了任意傅里叶变换算法与高斯快速傅里叶变换算法的计算效率与精度,在相近的数值精度下,本文算法波数选取少,效率高;测试CPU-GPU并行效果,结果表明相比CPU串行算法,CPU-GPU并行算法的计算效率大大提升,千万数量级节点数模型正演仅耗时数秒.最后利用实际地形数据进行三维重力异常场数值模拟,证明了新方法的高效性与实用性,对实现大规模复杂条件下重力异常精细化反演成像与综合解释有重要意义.
- 戴世坤朱德祥张莹李昆陈轻蕊凌嘉宣田红军