共查询到20条相似文献,搜索用时 453 毫秒
1.
2.
3.
在三维电阻率的正演计算中往往涉及到快速、准确求解大型线性方程纽Ax=b的问题。通过采用有限差分法来构造出求解点电源三维地电场的大型稀疏对称线性方程组。并引入Lanczos迭代技术,构造出三对角阵方程组,然后采用正交分解法进行求解,它是Krylov子空间方法中的一种。与传统迭代算法相比,它占用内存少,收敛速度快且稳定。针对大型稀疏矩阵及MATLAB语言的特点,采用简单记录矩阵的非零元素值及其所在行、列值的方法存储大型稀疏矩阵,可大大节省机器内存,提高运算速度。理论分析和计算实例显示,此算法是地电三维正演计算的有效方法,为下一步的反演计算打好基础。 相似文献
4.
在地震勘探中,地震正演模拟是非常重要的技术。与时间域正演相比,频率域正演速度快,计算效率高。如何高效准确地完成频率域正演计算是目前该领域的一个重要问题。数值频散问题和如何提高计算效率降低求解分解阻抗内存占用量一直是频率域正演所需要解决的问题。与传统的直接法求解阻抗矩阵的频率域正演方法不同,本文采用最小二乘共轭梯度法(LSCG法)求解阻抗矩阵进行频率域正演,并提出了一种波数补偿的表达式来压制数值频散现象。经过简单模型和复杂模型的数值测试,采用最小二乘共轭梯度法(LSCG法)求解阻抗矩阵进行频率域正演能够有效降低计算时间,且采用波数补偿的频率域正演方法能够有效压制数值频散现象,提高波场模拟精度。 相似文献
5.
利用ICCG迭代技术加快电阻率三维正演计算 总被引:2,自引:2,他引:0
一般而言,有限差分法求解点源三维地电场正问题所形成的大型稀疏线性方程组Ax=b,直接解法的计算效率极低。本文从系数矩阵A的不完全Cholesky分解及矩阵特征值的特点等角度,说明了不完全Cholesky共轭梯度(ICCG)迭代技术可大大提高电阻率三维正演速度的内在原因。结合矩阵A的稀疏存储模式,使得内存需求也大大减少。 相似文献
6.
一般而言,有限差分法求解点源三维地电场正问题所形成的大型稀疏线必方程组Ax=b,直接解法的计算效率极低。本文从系数矩阵A的不完全Cholesky分解及矩阵特征值的特点等角度,说明了不完全Cholesky共轭梯度(ICCG)迭代技术可大大提高电阻率三维正演速度的内在原因。结合矩阵A的稀疏存储模式,使得内存需求也大大减少。 相似文献
7.
正演是反演技术的基础,正演速度和求解反演问题的系数矩阵存放一直是起伏地形下重、磁三维反演的关键技术问题。这里提出了一种起伏地形下重磁快速正演计算方法,其计算原理是根据反演在垂向的剖分层数,利用水平地形正演计算形成二个不同大小刻度标尺矩阵,然后在模型空间,使用分段线性插值的方式,直接计算出起伏地形观测点的正演值。该方法的主要特点是在保持很高的计算精度下计算速度可提高二倍,且节省计算内存,适合起伏地形下重磁三维反演技术研究。 相似文献
8.
三维电阻率法对反演的精度和速度的要求越来越高,而正演是反演的基础,因此直流电阻率三维正演计算的速度和精度是三维电阻率反演实用化的关键。这里利用对称超松弛预条件共轭梯度法(SSOR-PCG),求解有限差分法离散生成的大型稀疏线性方程组,预条件矩阵的选择大大降低了系数矩阵的条件数,结合矩阵的一维非零元素压缩存储模式,使得正演计算速度得以提高,而内存占用量明显减小。在直流电阻率三维正演中采用异常场法,提高了电源点附近的解的精度。利用编制的有限差分正演程序,对两层模型、垂直接触带模型和低阻异常体模型进行了数值模拟,计算结果表明该算法是可行的,且可以明显提高正演计算的速度和精度。 相似文献
9.
基于有限单元法的二维/三维大地电磁正演模拟策略 总被引:1,自引:0,他引:1
对于二维和三维大地电磁正演问题,有限单元法最后形成了一个线性方程组KX=p。方程组中的K是大型稀疏的带状对称复系数矩阵,其条件数远大于1,为严重病态矩阵,求解其对应方程组会遇到很多困难。不完全LU分解处理的BICGSTAB算法,可用于该线性方程组的求解,并且具有速度快,精度高,稳定性好等优点。为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下,设计了非均匀网格剖分。在程序编制中,因只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间。通过对二维模型和三维模型电磁响应的计算,验证了该算法的正确性。 相似文献
10.
针对三维地电场正演数值计算过程中形成的超大规模稀疏线性方程组,在分析此类线性方程组的一般解法基础上,着重阐述一种适宜求解此类方程组的Lanczos迭代过程与算法原理。同时,当地下介质的电性差异较大时,形成系数矩阵A的条件数就很大,可对算法进行适当改进。讨论采用不完全Cholesky分解方法进行预条件处理,经过条件数改善后,形成新的线性方程组系数矩阵,就会变为一个近似的单位矩阵。经改进后的Lanczos算法,将提高数值计算稳定性,从而加快迭代收敛速度,为提高反演质量提供基础。 相似文献
11.
利用小波变换快速反演大型磁数据 总被引:4,自引:0,他引:4
我们应用小波变换反演大型磁数据,恢复三维磁化率分布问题,快速小波变换通过给出小波系数的门槛值来产生灵敏度矩阵表达式,在波域进行正演模拟计算,通过小波域内稀疏矩阵矢量的乘积运算,和压缩后的矩阵来实施快速正演模拟,所用CPU的时间的减少与矩阵压缩在正比,而且,我们采用对数阻尼法,用共轭梯度解算器来获得反演结果,这种求解手段的组合,直接减少了正演模拟的时间,进而加快了反演速度,产生了更的三维磁反演算法。 相似文献
12.
13.
全波形反演不仅利用相位和振幅信息,还利用波形的细节变化,具有刻画模型精确细节的能力.在对稀疏矩阵直接LU分解求解的基础上,采用梯度预处理方法对声波介质速度模型进行了反射波全波形反演.采用误差反向传播算法计算目标函数梯度以及伪Hessian矩阵的对角线元素来做梯度预处理.数值模型的实验结果表明,利用有效的频率段便能反演出分辨率较高的速度结构,用低频反演出的结果作为高频反演的初始模型,减少了解的非唯一性.二维高斯光滑初始模型提供了有利的低频信息,得到较好的反演结果.伪Hessian矩阵的预处理吸收了高斯牛顿法的二次收敛优势,在不增加计算量的前提下,加快收敛速度. 相似文献
14.
15.
积分方程法求解大型稀疏矩阵制约线源近区三维瞬变电磁正演模拟的速度,采用系统迭代法,把整个异常体区域看成一个系统,每个块状异常体在各自区域内形成子系统,单独处理各个异常体,并建立子系统间的关系,从而得到整个系统的值,这样不仅节省了内存和计算时间,而且降低了离散矩阵的条件数,提高计算速度和精度。通过与解析解的对比证明了该方法的正确性,简单和复杂三维体的二次水平电场和垂直磁场瞬态响应计算结果表明:①二次水平电场值反映了三维异常体的电性特征,同时其极值对应异常体的左右边界;②垂直磁场响应反映了异常体的纵向分辨率,曲线随偏移距增加而变缓。 相似文献
16.
17.
18.
目前有限元技术的开发及在电磁勘查技术中的应用已日趋成熟.然而,有限元正演模拟主要集中于地面和海洋电磁,航空电磁三维有限元模拟尚没有受到足够重视.以前人工作为基础,利用结构化网格实现了航空电磁系统的三维矢量有限元正演模拟.从二次场双旋度矢量非齐次亥姆霍兹方程出发,应用广义变分原理推导出变分方程,并采用六面体单元剖分,将场置于单元棱边上,对每个单元应用线性插值,最后合成含有稀疏矩阵的线性方程组.针对航空电磁多源性问题,利用MUMPS(multifrontal massively parallel sparse direct solver)直接求解器进行求解,在保证精度的前提下大幅度提高计算速度.利用单个异常体模型检验算法的精度和软件的稳定性,进而通过典型地电模型的模拟验证算法的有效性.对不同地下电性结构正演模拟结果进行对比分析,进一步研究了覆盖层和垂直接触带等典型构造对航空电磁响应的影响特征. 相似文献
19.
为了加快大地电磁三维正演的求解速度,本文将一种新型的代数多重网格算法——聚集多重网格(aggregation-based algebraic multigrid, AGMG)算法引入大地电磁三维正演模拟中。首先从准静态条件下的麦克斯韦方程出发,利用交错网格有限体积法进行离散,并采用第一类Dirichlet边界条件形成大型稀疏复线性方程组;然后阐述AGMG算法的粗化策略和套迭代技术,并实施3种不同的AGMG求解算法:1)传统的V循环AGMG算法;2)AGMG预处理共轭梯度(AGMG-CG);3)AGMG预处理广义共轭残差法(AGMG-GCR)。最终实现大地电磁法三维正演模拟。对典型地电模型进行正演模拟,并与已有的大地电磁三维正反演程序(ModEM)进行结果对比,以验证本文算法的准确性。另外,不同剖分网格和极化方式正演模拟结果与准残量最小化(QMR)迭代算法的对比表明,AGMG预处理求解算法(AGMG-CG、AGMG-GCR)不仅能够改善算法的稳定性,而且能够快速有效地求解正演问题;其中AGMG-GCR迭代次数更少,求解速度更快,误差衰减曲线更光滑,在144×152×104网格剖分情况下,相对于现有ModEM程序能够提高十几倍的计算速度,尤其适合大规模大地电磁三维正演问题。 相似文献
20.
系数矩阵存储和线性方程组求解是限制三维电磁积分方程方法发展的主要因素。Zhdanov提出准线性(QL)近似技术,建立了复杂散射场与背景场的线性关系,有效地避免了积分方程中大型线性方程组的求解,但是该算法用于多源问题航空电磁正演模拟时精度不高。因此,本文提出一种基于多重网格准线性(MGQL)近似的算法,并利用系数矩阵的Toeplitz性质存储矩阵和快速傅里叶变换,实现了矩阵与向量的快速乘积、降低了计算复杂度,采用多重网格结合了积分方程方法和准线性近似解法的优点,在保证精度的条件下提高计算速度、减少存储量。针对不同类型网格的模拟实验表明,相比于传统积分方程方法,本文算法在保证计算精度的同时,可以将计算速度极大地提高(>10倍)。 相似文献