基于局部加密等级网格的2.5D直流电法有限元模拟 |
| |
作者姓名: | 胡宏伶 肖晓 潘克家 汤井田 谢维 |
| |
作者单位: | 湖南师范大学数学与计算机科学学院;中南大学地球科学与信息物理学院;有色金属成矿预测教育部重点实验室;中南大学数学与统计学院;成都理工大学油气藏地质及开发工程国家重点实验室; |
| |
基金项目: | 国家自然科学基金资助项目(41204082,11301176,41104071);中国博士后科学基金特别资助项目(2013T60781);高等学校博士学科点专项科研基金资助项目(20120162120036);成都理工大学油气藏地质及开发工程国家重点实验室基金资助项目(PLC201305) |
| |
摘 要: | 合理截取半圆形计算区域,采取局部加密的?-等级网格,结合对称行索引存贮格式(CSR)及并行稀疏直接求解器PARDISO,提出一种高效、高精度的2.5D直流电阻率法有限元正演方案,并编制相应的Fortran程序,对具有解析解的3个典型地电模型进行计算与分析。研究结果表明:圆形截断边界不仅便于在径向方向上采取?-等级网格剖分,而且能大大简化有限元模拟中单元刚度矩阵的计算;结构化的等级网格避开了通常非结构化网格有限元计算时繁琐的网格剖分及总体刚度阵的集成过程,且能在不增大问题规模的前提下,显著提高2.5D直流电法正演源点附近的模拟精度;Intel MKL的PARDISO求解器能在普通PC机上5 s内求解电法正演有限元离散得到的100万阶稀疏线性方程组,可广泛用于各种地球物理正演问题。
|
关 键 词: | 直流电阻率 等级网格 有限元 局部加密 并行稀疏直接求解器 行索引存贮格式 |
本文献已被 CNKI 等数据库收录! |
|