期刊文献+
共找到8篇文章
< 1 >
每页显示 20 50 100
基于PARDISO直接求解器的三维自然电位正反演
1
作者 苏朝阳 沈金松 罗辉 《物探与化探》 CAS 2024年第2期451-460,共10页
近年来自然电位法在海底硫化物资源的勘探和评价中发挥了重要作用。本文开展的是基于PARDISO直接求解器的3D自然电位正反演算法研究。首先,利用有限体积法离散自然电位控制方程,采用PARDISO直接求解器提高正演计算的效率,通过数值解与... 近年来自然电位法在海底硫化物资源的勘探和评价中发挥了重要作用。本文开展的是基于PARDISO直接求解器的3D自然电位正反演算法研究。首先,利用有限体积法离散自然电位控制方程,采用PARDISO直接求解器提高正演计算的效率,通过数值解与解析解对比,验证了正演算法的可靠性。其次,在3D反演算法中考虑了地形因素,同时将最小支撑约束与深度加权加入目标泛函中,理论模型数据的反演结果很好地恢复了矿体的结构。最后,利用该算法对室内沙箱实验获得的自然电位数据进行反演,结果显示得出的电流密度异常与金属棒的位置基本一致。因此,本文提出的反演算法在未来大规模自然电位数据反演中具有重要作用。 展开更多
关键词 自然电位法 3D聚焦反演 pardiso直接求解器
下载PDF
轴对称地层中高分辨率阵列侧向测井信赖域反演法 被引量:10
2
作者 潘克家 汤井田 +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
基于矢量有限元的带地形大地电磁三维反演研究 被引量:10
3
作者 顾观文 李桐林 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2020年第6期2449-2465,共17页
研究了基于矢量有限元方法的大地电磁带地形三维反演算法并开发了三维反演计算程序代码.在大地电磁场正演数值模拟方面,采用并行直接稀疏求解器PARDISO且无需进行散度校正的快速正演方案,对典型地形模型,在中等规模计算条件下,与双共轭... 研究了基于矢量有限元方法的大地电磁带地形三维反演算法并开发了三维反演计算程序代码.在大地电磁场正演数值模拟方面,采用并行直接稀疏求解器PARDISO且无需进行散度校正的快速正演方案,对典型地形模型,在中等规模计算条件下,与双共轭梯度法(BICG)计算结果比较,发现PARDISO比BICG快10倍以上;通过理论模型试算,并与前人的有限元法计算结果对比,验证了带地形三维正演计算程序的正确性.在反演方面,本研究基于共轭梯度方法编写了大地电磁带地形三维反演代码,为了避免直接求取雅可比矩阵,将反演中的雅可比矩阵计算问题转为求解两次"拟正演"问题,进而将PARDISO的快速正演方案应用于"拟正演"问题的求解,以提高反演计算效率.利用开发的反演算法对多个带地形地电模型的合成数据进行了三维反演,反演结果能很好地重现理论模型的电性结构,验证了本文开发的三维反演算法的正确性和可靠性.最后,利用该算法反演了某矿区大地电磁实测数据,反演得到的三维电性结构清晰地反映了研究区的地电特征,将反演结果与该区已有地质资料结合进行解释,应用效果明显,进一步验证了本文算法的有效性. 展开更多
关键词 大地电磁 矢量有限元法 pardiso 三维反演 带地形
下载PDF
矿井风网与采空区流场边界耦合模型与求解技术 被引量:3
4
作者 吴奉亮 高亚超 常心坦 《中南大学学报(自然科学版)》 EI CAS CSCD 北大核心 2020年第8期2333-2342,共10页
为揭示煤矿巷道和采空区中整体风流的流动机理,通过将采空区漏风边界的流量添加到矿井通风网络,建立矿井风网与采空区流场的边界耦合模型。基于矿井通风网络的节点压力解算法与采空区渗流场的有限元求解法,研究矿井风网与采空区流场边... 为揭示煤矿巷道和采空区中整体风流的流动机理,通过将采空区漏风边界的流量添加到矿井通风网络,建立矿井风网与采空区流场的边界耦合模型。基于矿井通风网络的节点压力解算法与采空区渗流场的有限元求解法,研究矿井风网与采空区流场边界耦合模型的迭代求解流程。分析基于线性渗流模型求解采空区非线性渗流模型的方法及控制方程组,用PARDISO作为耦合模型中大型稀疏线性方程组的求解器的软件开发接口。采用VC++与ObjectARX编写相应的模拟软件(简称i-MVS)模拟2个算例,验证模型的有效性与性能。研究结果表明:采空区流场的模拟需要矿井通风网络的解算结果为其提供压力边界条件,对新的矿井通风网络的解算又需要采空区流场的模拟结果来提供漏风边界上的流量,两者之间形成了相互依赖的耦合关系;耦合模型中的新风网是一个含有边界流量的不闭合网络,必须采用矿井通风网络的结点压力法进行解算;耦合模型实现了矿井通风网络与多个采空区流场的集成求解,其计算量是常规采空区流场模拟的若干倍,采用PARDISO求解器实现了对耦合模型的高效求解;相比于常用的模拟煤矿采空区流场的方法,耦合模型的准确度提高了12%。 展开更多
关键词 矿井通风网络 采空区流场 有限元方法 节点压力法 耦合 pardiso
下载PDF
基于局部加密等级网格的2.5D直流电法有限元模拟 被引量:3
5
作者 胡宏伶 肖晓 +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
基于矢量有限元的大地电磁快速三维正演研究 被引量:2
6
作者 顾观文 武晔 石砚斌 《物探与化探》 CAS 北大核心 2020年第6期1387-1398,共12页
本文采用直接求解器PARDISO且无需散度校正的正演方案,求解矢量有限元法对应的大型线性方程组,获得不同地形(水平和起伏)条件下地电模型的大地电磁响应,提高了大地电磁三维正演计算的速度。在中等规模计算条件下,通过本文的计算方法与... 本文采用直接求解器PARDISO且无需散度校正的正演方案,求解矢量有限元法对应的大型线性方程组,获得不同地形(水平和起伏)条件下地电模型的大地电磁响应,提高了大地电磁三维正演计算的速度。在中等规模计算条件下,通过本文的计算方法与带散度校正的迭代求解法对比,计算速度可提高十倍以上。 展开更多
关键词 大地电磁 矢量有限元法 三维正演 pardiso
下载PDF
A Parallel FEA Computing Kernel for Building Structures
7
作者 Jin Duan Yungui Li +2 位作者 Xiaoming Chen Hu Qi Jianyun Sun 《Journal of Applied Mathematics and Physics》 2013年第6期26-30,共5页
With the rapid development of high-rise buildings and long-span structures in the recent years, high performance com- putation (HPC) is becoming more and more important, sometimes even crucial, for the design and cons... With the rapid development of high-rise buildings and long-span structures in the recent years, high performance com- putation (HPC) is becoming more and more important, sometimes even crucial, for the design and construction of com- plex building structures. To satisfy the engineering requirements of HPC, a parallel FEA computing kernel, which is designed typically for the analysis of complex building structures, will be presented and illustrated in this paper. This kernel program is based on the Intel Math Kernel Library (MKL) and coded by FORTRAN 2008 syntax, which is a parallel computer language. To improve the capability and efficiency of the computing kernel program, the parallel concepts of modern FORTRAN, such as elemental procedure, do concurrent, etc., have been applied extensively in coding and the famous PARDISO solver in MKL has been called to solve the Large-sparse system of linear equations. The ultimate objective of developing the computing kernel is to make the personal computer have the ability to analysis large building structures up to ten million degree of freedoms (DOFs). Up to now, the linear static analysis and dynamic analysis have been achieved while the nonlinear analysis, including geometric and material nonlinearity, has not been finished yet. Therefore, the numerical examples in this paper will be concentrated on demonstrating the validity and efficiency of the linear analysis and modal analysis for large FE models, while ignoring the verification of the nonlinear analysis capabilities. 展开更多
关键词 High Performance COMPUTING FINITE ELEMENT Analysis BUILDING Structure pardiso
下载PDF
时间域激发极化法三维有限元正演
8
作者 骆红梅 《世界有色金属》 2020年第23期176-177,共2页
本文采用三维有限元数值模拟的方法,通过在三维有限元模拟直流电的基础上加入COLE-COLE模型对时间域激发极化法进行了正演。其中时间域COLE-COLE模型的加入采用数字滤波的方法,有效的避免了计算γ函数带来的麻烦;对三维有限元单元合成... 本文采用三维有限元数值模拟的方法,通过在三维有限元模拟直流电的基础上加入COLE-COLE模型对时间域激发极化法进行了正演。其中时间域COLE-COLE模型的加入采用数字滤波的方法,有效的避免了计算γ函数带来的麻烦;对三维有限元单元合成后的稀疏对称矩阵采用CSR存储格式和PARDISO并行求解方法。试算了模型,结果证明本文所采用的方法计算速度高,精度可靠。 展开更多
关键词 时间域激发极化法 三维有限元正演 COLE-COLE CSR pardiso
下载PDF
上一页 1 下一页 到第
使用帮助 返回顶部