重力梯度张量的预条件共轭梯度法反演 被引量：10

1

作者 陈少华 朱自强 +1 位作者 鲁光银 曹书锦 《中南大学学报（自然科学版）》 EI CAS CSCD 北大核心 2013年第2期619-625,共7页

针对重力梯度张量反演中的问题,提出基于预条件共轭梯度法的重力梯度张量反演。通过在目标函数中加入粗糙度对模型进行约束以避免反演参数远多于采集点数的欠定问题不稳定,并在目标函数中添加深度加权矩阵对核函数进行补偿,以避免核函... 针对重力梯度张量反演中的问题,提出基于预条件共轭梯度法的重力梯度张量反演。通过在目标函数中加入粗糙度对模型进行约束以避免反演参数远多于采集点数的欠定问题不稳定,并在目标函数中添加深度加权矩阵对核函数进行补偿,以避免核函数随着深度的增大而快速衰减的问题。分别反演、比较各重力梯度张量分量和联合5个独立分量,并将重力梯度张量5个独立联合反演应用于Y型岩脉。研究结果表明:联合反演效果明显优于单一分量的反演效果,且能较好地与原始模型相吻合,证明了本文算法的有效性。 展开更多

关键词 重力梯度张量 深度加权函数 预条件共轭梯度法 反演

在线阅读 下载PDF 职称材料

基于GPU的稀疏线性系统的预条件共轭梯度法 被引量：11

2

作者 张健飞 沈德飞 《计算机应用》 CSCD 北大核心 2013年第3期825-829,共5页

研究了基于GPU的稀疏线性方程组的预条件共轭梯度法加速求解问题,并基于统一计算设备架构(CUDA)平台编制了程序,在NVIDIAGT430 GPU平台上进行了程序性能测试和分析。稀疏矩阵采用压缩稀疏行(CSR)格式压缩存储,针对预条件共轭梯度法的算... 研究了基于GPU的稀疏线性方程组的预条件共轭梯度法加速求解问题,并基于统一计算设备架构(CUDA)平台编制了程序,在NVIDIAGT430 GPU平台上进行了程序性能测试和分析。稀疏矩阵采用压缩稀疏行(CSR)格式压缩存储,针对预条件共轭梯度法的算法特性,研究了基于GPU的稀疏矩阵与向量相乘的性能优化、数据从CPU端传到GPU端的加速传输措施。将编制的稀疏矩阵与向量相乘的kernel函数和CUSPARSE函数库中的cusparseDcsrmv函数性能进行了对比,最优得到了2.1倍的加速效果。对于整个预条件共轭梯度法,通过自编kernel函数来实现的算法较之采用CUBLAS库和CUSPARSE库实现的算法稍具优势,与CPU端的预条件共轭梯度法相比,最优可以得到7.4倍的加速效果。 展开更多

关键词 图形处理器 稀疏线性方程组 预条件共轭梯度法 压缩稀疏行 统一计算设备架构

在线阅读 下载PDF 职称材料

预条件共轭梯度法在线天线阵列RCS分析中的应用 被引量：4

3

作者 周后型 童创明 洪伟 《应用科学学报》 CAS CSCD 2001年第2期145-148,共4页

为计算线天线阵列的雷达散射截面 ,用脉冲函数为基函数的 Galerkin法将线天线阵列的电场积分方程组转化为一个线性方程组 ,其系数矩阵是一个对称的二重复 Toeplitz矩阵 .用预条件共轭梯度法与快速傅里叶变换的组合算法求解该线性方程组 ... 为计算线天线阵列的雷达散射截面 ,用脉冲函数为基函数的 Galerkin法将线天线阵列的电场积分方程组转化为一个线性方程组 ,其系数矩阵是一个对称的二重复 Toeplitz矩阵 .用预条件共轭梯度法与快速傅里叶变换的组合算法求解该线性方程组 ,其中预条件器采用 T.Chan的循环预条件器的分块形式 .这样的处理方法降低了对计算机内存的需求 ,加快了迭代速度 。 展开更多

关键词 线天线阵列 雷达散射截面 预条件共轭梯度法 傅立叶变换 RCS分析 线性方程组 迭代速度

在线阅读 下载PDF 职称材料

求解大型稀疏线性方程组的不完全SAOR预条件共轭梯度法 被引量：4

4

作者 温瑞萍 孟国艳 王川龙 《工程数学学报》 CSCD 北大核心 2007年第4期712-718,共7页

预条件共轭梯度法是求解大型稀疏线性方程组的有效方法之一,SSOR预条件方法是基于矩阵分裂的较有效的预条件共轭梯度法。通过矩阵分裂,本文讨论不完全SAOR预条件方法,研究此方法的预条件因子及系数矩阵的预条件数,并证明了此方法的预条... 预条件共轭梯度法是求解大型稀疏线性方程组的有效方法之一,SSOR预条件方法是基于矩阵分裂的较有效的预条件共轭梯度法。通过矩阵分裂,本文讨论不完全SAOR预条件方法,研究此方法的预条件因子及系数矩阵的预条件数,并证明了此方法的预条件数小于SSOR预条件方法的预条件数。最后通过求解离散化波松(Poisson)方程组表明了该方法的有效性。 展开更多

关键词 不完全SAOR 预条件共轭梯度法 条件数

在线阅读 下载PDF 职称材料

地震波阻抗反演的预条件共轭梯度法 被引量：2

5

作者 黄文松 年静波 刘喜武 《新疆石油地质》 CAS CSCD 北大核心 2006年第4期425-428,共4页

研究地震波阻抗线性反演求解方法。基于Cauchy准则,正则化稀疏反演问题,应用预条件共轭梯度法实现反射系数和子波同时迭代反演。在迭代求解正则化方程时,用共轭梯度法求解相应的原问题。给出初始子波估计的策略和波阻抗求解方法。用实... 研究地震波阻抗线性反演求解方法。基于Cauchy准则,正则化稀疏反演问题,应用预条件共轭梯度法实现反射系数和子波同时迭代反演。在迭代求解正则化方程时,用共轭梯度法求解相应的原问题。给出初始子波估计的策略和波阻抗求解方法。用实际数据检验算法,表明预条件共轭梯度法反演的子波、反射系数和波阻抗要比直接稀疏反演精度高、收敛快、数值稳定。 展开更多

关键词 地震波阻抗 Cauchy准则 迭代反演 预条件共轭梯度

在线阅读 下载PDF 职称材料

预条件共轭梯度法在辐射和散射问题中的应用 被引量：4

6

作者 周后型 洪伟 《电波科学学报》 EI CSCD 2000年第4期501-505,共5页

用矩量法求解一些辐射和散射问题 ,如线天线辐射和线状体散射等问题时 ,可以产生一个 Toeplitz线性方程组 ,采用预条件共轭梯度法 (PCG)与快速富里叶变换 (FFT)的结合方法 (PCGFFT)来求解该方程组 ,其中预条件器采用 T.Chan的优化循环... 用矩量法求解一些辐射和散射问题 ,如线天线辐射和线状体散射等问题时 ,可以产生一个 Toeplitz线性方程组 ,采用预条件共轭梯度法 (PCG)与快速富里叶变换 (FFT)的结合方法 (PCGFFT)来求解该方程组 ,其中预条件器采用 T.Chan的优化循环预条件器。使用 PCGFFT算法 ,可有效地节省内存 ,提高了计算速度。为说明其有效性 ,将 PCGFFT算法与 CGFFT算法以及 Levinson递推算法进行了对比。 展开更多

关键词 预条件共轭梯度法 电磁散射 电磁辐射

在线阅读 下载PDF 职称材料

预条件共轭梯度法在拱坝有限元重分析中的应用 被引量：1

7

作者 孙林松 郭兴文 李春和 《河海大学学报（自然科学版）》 CAS CSCD 北大核心 2009年第2期236-239,共4页

以初始设计的劲度矩阵为预条件矩阵,给出了大型结构有限元重分析的预条件共轭梯度算法.该算法不需要形成和存储修改结构的劲度矩阵,占用内存小,并具有较高的精度和收敛速度.拱坝体形修改有限元分析算例表明,即使设计变量有较大改变时,... 以初始设计的劲度矩阵为预条件矩阵,给出了大型结构有限元重分析的预条件共轭梯度算法.该算法不需要形成和存储修改结构的劲度矩阵,占用内存小,并具有较高的精度和收敛速度.拱坝体形修改有限元分析算例表明,即使设计变量有较大改变时,该方法也能较快地收敛到精确解. 展开更多

关键词 预条件共轭梯度法 有限元重分析 拱坝

在线阅读 下载PDF 职称材料

预条件共轭梯度法在大型振子阵列天线RCS分析中的应用

8

作者 周后型 洪伟 童创明 《电子学报》 EI CAS CSCD 北大核心 2001年第12期1601-1604,共4页

本文提出三重Toeplitz线性方程组预条件共轭梯度法 ,并将该法与快速付里叶变换 (FFT)结合 .这种结合算法称为PCGFFT .将PCGFFT应用于振子阵列天线的RCS分析中 .由于预条件器的使用 ,系数矩阵的条件数得到了很大改善 .数值结果表明 ,PCG... 本文提出三重Toeplitz线性方程组预条件共轭梯度法 ,并将该法与快速付里叶变换 (FFT)结合 .这种结合算法称为PCGFFT .将PCGFFT应用于振子阵列天线的RCS分析中 .由于预条件器的使用 ,系数矩阵的条件数得到了很大改善 .数值结果表明 ,PCGFFT不仅降低了对计算机内存的需求 。 展开更多

关键词 振动阵列天线 雷达散射截面 预条件共轭梯度法

在线阅读 下载PDF 职称材料

预条件共轭梯度法求解三维地电场有限元方程的网格分析

9

作者 王威 李景富 刘洁 《中山大学学报（自然科学版）》 CAS CSCD 北大核心 2018年第3期1-6,共6页

为了优化预条件共轭梯度法求解三维地电场有限元方程的效率,通过系统的网格分析,提出在常规六面体网格基础上二次剖分得到四面体网格的方案。模型分析结果表明,对于均匀网格,不完全Cholesky共轭梯度法(ICCG)和超松弛预条件共轭梯度法(SO... 为了优化预条件共轭梯度法求解三维地电场有限元方程的效率,通过系统的网格分析,提出在常规六面体网格基础上二次剖分得到四面体网格的方案。模型分析结果表明,对于均匀网格,不完全Cholesky共轭梯度法(ICCG)和超松弛预条件共轭梯度法(SORCG)均可成功求解;对于非均匀网格,六面体剖分会导致ICCG的预条件因子不符合条件而求解失败,但采用新的四面体剖分ICCG不仅成功求解而且相对六面体网格可以节省约50%的内存需求。 展开更多

关键词 网格分析 预条件共轭梯度 有限元 三维地电场

在线阅读 下载PDF 职称材料

一种改进的二维MT预条件非线性共轭梯度反演方法 被引量：5

10

作者 相鹏 《中国石油大学学报（自然科学版）》 EI CAS CSCD 北大核心 2014年第4期42-49,共8页

在大地电磁反演方法中反演精度与计算效率问题是一对矛盾,高斯牛顿类方法反演精度高但计算效率低,非线性共轭梯度类方法计算效率高,但是反演精度不如高斯牛顿法高。在前人研究的基础上,提出一种改进的预条件非线性共轭梯度法,通过构建... 在大地电磁反演方法中反演精度与计算效率问题是一对矛盾,高斯牛顿类方法反演精度高但计算效率低,非线性共轭梯度类方法计算效率高,但是反演精度不如高斯牛顿法高。在前人研究的基础上,提出一种改进的预条件非线性共轭梯度法,通过构建性状更接近高斯牛顿Hessian矩阵的预条件算子提高反演精度和计算速度。同时采用正则化参数的自适应更新算法保证反演稳定性和反演精度的平衡。模型实验验证了该方法的正确性。与其他方法的对比结果表明,该方法在保证反演精度的同时,提高了计算效率。对中国西部某地的实测MT数据进行处理解释的结果表明,该方法在解决复杂构造问题方面具有较高的实用价值。 展开更多

关键词 大地电磁 预条件非线性共轭梯度 正则化参数 自适应算法 反演

在线阅读 下载PDF 职称材料

三维地震数据离散光滑插值的共轭梯度法 被引量：12

11

作者 李冰 刘洪 李幼铭 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2002年第5期691-699,共9页

针对三维地震数据插值 ,提出采用Laplacian算子进行光滑约束的插值方法 ,并借鉴Mallet研究的离散光滑插值思路 ,采用预条件共轭梯度法 ,直接生成网格节点上的值 ,从而回避寻求满足插值方程的函数 .为了实现其中Laplacian算子的快速求逆 ... 针对三维地震数据插值 ,提出采用Laplacian算子进行光滑约束的插值方法 ,并借鉴Mallet研究的离散光滑插值思路 ,采用预条件共轭梯度法 ,直接生成网格节点上的值 ,从而回避寻求满足插值方程的函数 .为了实现其中Laplacian算子的快速求逆 ,文中引入Claerbout螺旋坐标系谱因式分解理论 .在螺旋坐标系下 ,Laplacian算子的表示矩阵具有Toeplitz结构 ,其快速求逆可由谱法LU分解实现 .基于二维离散光滑插值 ,文中还给出共轭梯度法与NMO相结合的沿时间切片逐层处理的离散光滑插值流程 .最后 。 展开更多

关键词 三维地震数据 LAPLACIAN算子 离散光滑道插值 预条件共轭梯度 螺旋坐标系谱因式分解 NMO

在线阅读 下载PDF 职称材料

基于GPU的SSOR稀疏近似逆预条件研究 被引量：2

12

作者 高家全 王志超 《浙江工业大学学报》 CAS 北大核心 2016年第2期140-145,共6页

由于SSOR预条件共轭梯度算法中预条件方程求解需要前推和回代,导致算法迁移到GPU平台上并行效率不高.为此,基于诺依曼多项式分解技术,提出了一种GPU加速的SSOR稀疏近似逆预条件子(GSSORSAI).它不仅保持了原线性系统系数矩阵的稀疏和对... 由于SSOR预条件共轭梯度算法中预条件方程求解需要前推和回代,导致算法迁移到GPU平台上并行效率不高.为此,基于诺依曼多项式分解技术,提出了一种GPU加速的SSOR稀疏近似逆预条件子(GSSORSAI).它不仅保持了原线性系统系数矩阵的稀疏和对称正定特性,而且预条件方程求解仅需一次稀疏矩阵矢量乘运算,避免了前推和回代过程.实验结果表明:在NVIDIA Tesla C2050GPU上,对比使用Python在单个CPU上SSOR稀疏近似逆预条件子实现方法,GSSORSAI平均快将近100倍;应用到并行的PCG算法中,相比无预条件的CG算法,平均提高了算法的3倍的收敛速度. 展开更多

关键词 SSOR预条件子 预条件共轭梯度算法 稀疏近似逆 GPU

在线阅读 下载PDF 职称材料

矩阵B^TA^(-1)B的特征值估计及预条件处理

13

作者 李铮 邵新慧 李长军 《东北大学学报（自然科学版）》 EI CAS CSCD 北大核心 2005年第6期603-605,共3页

在矩阵A为对称正定和矩阵B为列满秩的假设下,研究矩阵BTA-1B的特征值上下界估计,进而给出了BTA-1B的谱条件数的估计·基于以上论述,论证了当矩阵A的条件较好时矩阵Q=BTB可作为矩阵BTA-1B的预条件矩阵·在数值实验中,采用预条件... 在矩阵A为对称正定和矩阵B为列满秩的假设下,研究矩阵BTA-1B的特征值上下界估计,进而给出了BTA-1B的谱条件数的估计·基于以上论述,论证了当矩阵A的条件较好时矩阵Q=BTB可作为矩阵BTA-1B的预条件矩阵·在数值实验中,采用预条件共轭梯度算法(PCG)对Stokes方程求解,实验结果表明Q=BTB确实是一类有效的预条件矩阵·这一结果也和其他文献的数值结果相吻合· 展开更多

关键词 特征值 谱条件数 预条件矩阵 广义SOR算法 预条件共轭梯度算法 STOKES方程

在线阅读 下载PDF 职称材料

基于预仿真策略的变电站培训仿真系统 被引量：5

14

作者 张炳达 李银芝 《电力自动化设备》 EI CSCD 北大核心 2013年第1期142-147,共6页

为降低数字物理混合型变电站培训仿真系统的建设成本,对一次系统仿真模型及其求解算法进行了研究。利用贝瑞隆法将整个仿真模型分成站内仿真模块、线路仿真模块和负荷仿真模块,通过短路故障转移法与电流逻辑分析法使站内仿真模块的规模... 为降低数字物理混合型变电站培训仿真系统的建设成本,对一次系统仿真模型及其求解算法进行了研究。利用贝瑞隆法将整个仿真模型分成站内仿真模块、线路仿真模块和负荷仿真模块,通过短路故障转移法与电流逻辑分析法使站内仿真模块的规模大幅缩小,采用多项式预条件共轭梯度法提高仿真模型的求解速度。为解决一次系统网络参数突然变化带来的计算风暴,提出了一种计算模块独立于输出模块的预仿真策略。实例测试表明,基于预仿真策略的数字式一次系统能在廉价的多微处理器系统中运行,且能与真实二次设备很好地配合。 展开更多

关键词 变电站 培训 仿真 预条件共轭梯度法 贝瑞隆电路

在线阅读 下载PDF 职称材料

ω循环型边界条件 被引量：4

15

作者 梅金顺 刘洪 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2003年第6期835-841,共7页

通过利用预条件共轭梯度法对对称正定Toeplitz矩阵系统进行分析 ,重点介绍了一种新的嵌入式预条件矩阵构造方法 ,证明了以前的预条件矩阵构造方法大都是这种方法的特例 .提出了ω循环型边界条件 ,并将其与普通循环型及螺旋型边界条件作... 通过利用预条件共轭梯度法对对称正定Toeplitz矩阵系统进行分析 ,重点介绍了一种新的嵌入式预条件矩阵构造方法 ,证明了以前的预条件矩阵构造方法大都是这种方法的特例 .提出了ω循环型边界条件 ,并将其与普通循环型及螺旋型边界条件作了分析、比较后得到了一种新的边界条件即混合型边界条件 . 展开更多

关键词 ω循环 边界条件 预条件共轭梯度法 TOEPLITZ矩阵

在线阅读 下载PDF 职称材料

海量卫星重力梯度观测数据确定地球重力位模型的数值方法 被引量：1

16

作者 朱广彬 常晓涛 +2 位作者 邹贤才 徐新禹 王建强 《大地测量与地球动力学》 CSCD 北大核心 2011年第6期140-144,共5页

基于空域最小二乘法,对卫星重力梯度数据确定地球重力场中的Cholesky分解法、预条件共轭梯度法以及OpenMP并行算法3种数值方法进行比较与分析。研究表明,在计算机硬件资源有限的情况下,传统的Cholesky分解法已经无法满足求解要求;预条... 基于空域最小二乘法,对卫星重力梯度数据确定地球重力场中的Cholesky分解法、预条件共轭梯度法以及OpenMP并行算法3种数值方法进行比较与分析。研究表明,在计算机硬件资源有限的情况下,传统的Cholesky分解法已经无法满足求解要求;预条件共轭梯度法的求解效率较之Cholesky分解法有改进,但其以损失小量精度为代价;OpenMP并行算法在不损失求解精度的条件下,可提高求解的效率。 展开更多

关键词 卫星重力梯度 CHOLESKY分解 预条件共轭梯度 OpenMP并行算法 数据处理

在线阅读 下载PDF 职称材料

基于不等式约束的最小二乘法三维电阻率反演及其算法优化 被引量：75

17

作者 刘斌 李术才 +5 位作者 李树忱 聂利超 钟世航 李利平 宋杰 刘征宇 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2012年第1期260-268,共9页

基于光滑约束的最小二乘法是三维电阻率反演的主要方法,但该方法在某些情况下存在着多解性较强的问题,且普遍耗时较长,严重制约了三维反演方法的推广与发展.为改善上述问题,将表征模型参数变化范围的不等式约束作为先验信息引入最小二... 基于光滑约束的最小二乘法是三维电阻率反演的主要方法,但该方法在某些情况下存在着多解性较强的问题,且普遍耗时较长,严重制约了三维反演方法的推广与发展.为改善上述问题,将表征模型参数变化范围的不等式约束作为先验信息引入最小二乘线性反演方法中,有效地改善了反演结果的精度,降低了反演的多解性问题.为了解决耗时较长的问题,基于预条件共轭梯度(PCG)算法和Cholesky分解法的特点提出了一套优化三维电阻率反演计算效率的计算方案.在该方案中,Cholesky分解法被用来求解敏感度矩阵计算中的多个点源场的正演问题,Cholesky分解法只需对总体系数矩阵进行一次分解,然后对不同的右端向量进行回代即可.将预条件共轭梯度法引入到三维电阻率反演方程的求解中,将雅可比迭代中的对角阵作为预处理矩阵,其具有求逆方便、无需内存空间的特点,有效地加快了收敛速度.对合成数据以及实测数据的反演算例表明,借助不等式约束和反演效率优化方案,最小二乘反演方法可得到较为精确的反演结果,有效地提高了反演计算效率,具有良好的推广前景. 展开更多

关键词 三维电阻率反演 最小二乘反演 不等式约束 预条件共轭梯度算法 Cholesky分解算法 电阻率层析成像

在线阅读 下载PDF 职称材料

基于有限差分和有限体积法相结合的油浸式变压器三维温度场计算 被引量：55

18

作者 王永强 马伦 +3 位作者 律方成 毕建刚 王柳 万涛 《高电压技术》 EI CAS CSCD 北大核心 2014年第10期3179-3185,共7页

为准确计算油浸变压器的温度场,在对比有限体积法(FVM)和有限差分法(FDM)各自优点的基础上,提出了一种计算变压器温度场的混合计算方法。以驻马店金桥变电站的一台31.5 MVA主变为例,建立了3维模型,利用有限体积和有限差分混合计算方法,... 为准确计算油浸变压器的温度场,在对比有限体积法(FVM)和有限差分法(FDM)各自优点的基础上,提出了一种计算变压器温度场的混合计算方法。以驻马店金桥变电站的一台31.5 MVA主变为例,建立了3维模型,利用有限体积和有限差分混合计算方法,并在绕组和铁芯的热传导计算中采用了高精度的预条件共轭梯度解法(PCGM),计算了变压器的3维温度场。将利用混合方法的平均温度计算结果分别与实际监测数据、使用同样模型的FDM及FVM算法计算出的数据、以及使用IEEE导则模型计算出的数据进行了对比。对比结果表明:混合算法的计算精度明显高于利用IEEE导则模型、FDM算法和FVM算法的计算结果,证明了这种混合计算方法能更准确地计算油浸式变压器的内部温度场。 展开更多

关键词 油浸式变压器 3维温度场 绕组温度 有限差分法 有限体积法 预条件共轭梯度法

在线阅读 下载PDF 职称材料

连续电性介质线源井-地电位三维有限元数值模拟 被引量：17

19

作者 刘海飞 陈德鹏 +2 位作者 戴前伟 柳建新 阮百尧 《桂林理工大学学报》 CAS 北大核心 2011年第1期28-38,共11页

为了更好的研究油藏开发过程中剩余油的分布情况以及注水前沿的推进位置,给出了三维空间任意线电流源且电导率连续变化的井-地电位三维有限元数值模拟方法。在模拟过程中,采用异常电位法消除线电流源奇异性的影响,并且将网格节点间的电... 为了更好的研究油藏开发过程中剩余油的分布情况以及注水前沿的推进位置,给出了三维空间任意线电流源且电导率连续变化的井-地电位三维有限元数值模拟方法。在模拟过程中,采用异常电位法消除线电流源奇异性的影响,并且将网格节点间的电导率设计为连续变化,使其能够模拟实际中更为复杂的地质情况;考虑到采油井或注水井可能存在垂直、倾斜及水平等任意情况,推导了任意线电流源正常电位的计算公式;采用对称超松弛预条件共轭梯度法解大型线性方程组,解决了井-地电位三维模拟速度慢的问题,并编制了实用、高效的井-地电位三维模拟软件。通过对3例地电模型进行模拟计算,分析了水泥环胶结质量、回流井到观测井的距离以及目标层的深度对井-地电位观测结果的影响。 展开更多

关键词 井-地电位 三维有限元数值模拟 异常电位 任意线电流源 预条件共轭梯度法

在线阅读 下载PDF 职称材料

基于自适应加权光滑约束与PCG算法的三维电阻率探测反演成像 被引量：11

20

作者 刘斌 李术才 +5 位作者 聂利超 李利平 刘征宇 宋杰 隋斌 周宗青 《岩土工程学报》 EI CAS CSCD 北大核心 2012年第9期1646-1653,共8页

三维电阻率探测在实际工程中日益受到重视,但存在着深部分辨率低、反演效率低等关键问题,严重制约了三维电阻率探测技术的应用和推广。针对深部分辨率低的问题,提出了一种随三维模型网格尺寸自适应调整的加权光滑约束,改善了深部网格的... 三维电阻率探测在实际工程中日益受到重视,但存在着深部分辨率低、反演效率低等关键问题,严重制约了三维电阻率探测技术的应用和推广。针对深部分辨率低的问题,提出了一种随三维模型网格尺寸自适应调整的加权光滑约束,改善了深部网格的电阻率差异容许程度,实现了对不同深度网格约束的差异化加权处理,有效的提高了深部反演的分辨率和成像效果。针对反演计算耗时长、效率低的问题,基于预条件共轭梯度法求解快速稳定的优势,提出了三维电阻率快速稳定反演成像算法。在该算法中,将雅可比迭代中的对角阵作为预条件矩阵,其具有求逆方便、无需内存空间的特点,显著加快了收敛速度。最后,利用合成算例和隧道导水裂隙探测的工程实例验证了上述反演方法的可行性与有效性,表明借助于自适应调整加权光滑约束和预条件共轭梯度算法,有效的提高了深部分辨率和计算效率,显著改善了反演效果。 展开更多

关键词 三维电阻率反演 自适应加权光滑约束 深部分辨率 预条件共轭梯度算法 导水裂隙探测

在线阅读 下载PDF 职称材料