期刊文献+
共找到4篇文章
< 1 >
每页显示 20 50 100
轴对称地层中高分辨率阵列侧向测井信赖域反演法 被引量:11
1
作者 潘克家 汤井田 +1 位作者 杜华坤 蔡志杰 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2016年第8期3110-3120,共11页
本文研究轴对称地层中高分辨率阵列侧向测井(HRLA)的多参数信赖域反演方法.首先对前期HRLA的有限元正演方法进行改进,提出基于叠加原理和并行直接稀疏求解器PARDISO的快速正演方案,更适合于反演计算.将HRLA反演问题转化为非线性最小二... 本文研究轴对称地层中高分辨率阵列侧向测井(HRLA)的多参数信赖域反演方法.首先对前期HRLA的有限元正演方法进行改进,提出基于叠加原理和并行直接稀疏求解器PARDISO的快速正演方案,更适合于反演计算.将HRLA反演问题转化为非线性最小二乘问题,利用信赖域算法求解.为提高反演速度,推导了目标函数对优化参数偏导数的具体计算公式.对典型地层模型,与已有文献中Jacobi预条件共轭梯度法(JCG)计算结果比较,发现PARDIDO比JCG快10倍以上,验证了本文正演程序的正确性与高效性.利用信赖域算法求解了电阻率四参数反演和传统的三参数反演.研究结果表明:并行直接稀疏求解器PARDISO能有效求解此类HRLA正演问题,对6次不同探测深度的测井模拟,所形成的有限元刚度矩阵完全相同,只须进行一次矩阵分解,大大加快了正反演的速度.信赖域算法收敛速度快,且具有全局收敛性.HRLA的信赖域反演结果几乎不依赖于初值的选取,从较差初值出发仍能得到满意的反演结果.另外信赖域算法抗噪能力比较强,即使对测井数据添加信噪比为10dB(甚至5dB)的高斯白噪声,仍能通过反演得到较为准确的地层参数. 展开更多
关键词 阵列侧向测井 信赖域算法 反演 全局收敛 pardiso
在线阅读 下载PDF
基于矢量有限元的带地形大地电磁三维反演研究 被引量:13
2
作者 顾观文 李桐林 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2020年第6期2449-2465,共17页
研究了基于矢量有限元方法的大地电磁带地形三维反演算法并开发了三维反演计算程序代码.在大地电磁场正演数值模拟方面,采用并行直接稀疏求解器PARDISO且无需进行散度校正的快速正演方案,对典型地形模型,在中等规模计算条件下,与双共轭... 研究了基于矢量有限元方法的大地电磁带地形三维反演算法并开发了三维反演计算程序代码.在大地电磁场正演数值模拟方面,采用并行直接稀疏求解器PARDISO且无需进行散度校正的快速正演方案,对典型地形模型,在中等规模计算条件下,与双共轭梯度法(BICG)计算结果比较,发现PARDISO比BICG快10倍以上;通过理论模型试算,并与前人的有限元法计算结果对比,验证了带地形三维正演计算程序的正确性.在反演方面,本研究基于共轭梯度方法编写了大地电磁带地形三维反演代码,为了避免直接求取雅可比矩阵,将反演中的雅可比矩阵计算问题转为求解两次"拟正演"问题,进而将PARDISO的快速正演方案应用于"拟正演"问题的求解,以提高反演计算效率.利用开发的反演算法对多个带地形地电模型的合成数据进行了三维反演,反演结果能很好地重现理论模型的电性结构,验证了本文开发的三维反演算法的正确性和可靠性.最后,利用该算法反演了某矿区大地电磁实测数据,反演得到的三维电性结构清晰地反映了研究区的地电特征,将反演结果与该区已有地质资料结合进行解释,应用效果明显,进一步验证了本文算法的有效性. 展开更多
关键词 大地电磁 矢量有限元法 pardiso 三维反演 带地形
在线阅读 下载PDF
矿井风网与采空区流场边界耦合模型与求解技术 被引量:6
3
作者 吴奉亮 高亚超 常心坦 《中南大学学报(自然科学版)》 EI CAS CSCD 北大核心 2020年第8期2333-2342,共10页
为揭示煤矿巷道和采空区中整体风流的流动机理,通过将采空区漏风边界的流量添加到矿井通风网络,建立矿井风网与采空区流场的边界耦合模型。基于矿井通风网络的节点压力解算法与采空区渗流场的有限元求解法,研究矿井风网与采空区流场边... 为揭示煤矿巷道和采空区中整体风流的流动机理,通过将采空区漏风边界的流量添加到矿井通风网络,建立矿井风网与采空区流场的边界耦合模型。基于矿井通风网络的节点压力解算法与采空区渗流场的有限元求解法,研究矿井风网与采空区流场边界耦合模型的迭代求解流程。分析基于线性渗流模型求解采空区非线性渗流模型的方法及控制方程组,用PARDISO作为耦合模型中大型稀疏线性方程组的求解器的软件开发接口。采用VC++与ObjectARX编写相应的模拟软件(简称i-MVS)模拟2个算例,验证模型的有效性与性能。研究结果表明:采空区流场的模拟需要矿井通风网络的解算结果为其提供压力边界条件,对新的矿井通风网络的解算又需要采空区流场的模拟结果来提供漏风边界上的流量,两者之间形成了相互依赖的耦合关系;耦合模型中的新风网是一个含有边界流量的不闭合网络,必须采用矿井通风网络的结点压力法进行解算;耦合模型实现了矿井通风网络与多个采空区流场的集成求解,其计算量是常规采空区流场模拟的若干倍,采用PARDISO求解器实现了对耦合模型的高效求解;相比于常用的模拟煤矿采空区流场的方法,耦合模型的准确度提高了12%。 展开更多
关键词 矿井通风网络 采空区流场 有限元方法 节点压力法 耦合 pardiso
在线阅读 下载PDF
基于局部加密等级网格的2.5D直流电法有限元模拟 被引量:3
4
作者 胡宏伶 肖晓 +2 位作者 潘克家 汤井田 谢维 《中南大学学报(自然科学版)》 EI CAS CSCD 北大核心 2014年第7期2259-2267,共9页
合理截取半圆形计算区域,采取局部加密的?-等级网格,结合对称行索引存贮格式(CSR)及并行稀疏直接求解器PARDISO,提出一种高效、高精度的2.5D直流电阻率法有限元正演方案,并编制相应的Fortran程序,对具有解析解的3个典型地电模型进行计... 合理截取半圆形计算区域,采取局部加密的?-等级网格,结合对称行索引存贮格式(CSR)及并行稀疏直接求解器PARDISO,提出一种高效、高精度的2.5D直流电阻率法有限元正演方案,并编制相应的Fortran程序,对具有解析解的3个典型地电模型进行计算与分析。研究结果表明:圆形截断边界不仅便于在径向方向上采取?-等级网格剖分,而且能大大简化有限元模拟中单元刚度矩阵的计算;结构化的等级网格避开了通常非结构化网格有限元计算时繁琐的网格剖分及总体刚度阵的集成过程,且能在不增大问题规模的前提下,显著提高2.5D直流电法正演源点附近的模拟精度;Intel MKL的PARDISO求解器能在普通PC机上5 s内求解电法正演有限元离散得到的100万阶稀疏线性方程组,可广泛用于各种地球物理正演问题。 展开更多
关键词 直流电阻率 等级网格 有限元 局部加密 并行稀疏直接求解器 行索引存贮格式
在线阅读 下载PDF
上一页 1 下一页 到第
使用帮助 返回顶部