共查询到18条相似文献，搜索用时 140 毫秒
1.
隧道施工期超前探测对于避免突涌水灾害的发生具有重要作用,为满足隧道三维电阻率超前探测快速化解译与成像的要求,本文提出了一种基于GPU并行的蚁群算法与最小二乘方法相结合的混合反演算法.该方法结合线性反演与非线性反演的优点,利用蚁群算法全局搜索能力强的优点为最小二乘反演提供较优的初始模型,以克服最小二乘算法容易陷入局部最优的缺点,提高了隧道三维电阻率反演成像的精度.同时,基于蚁群算法的天然并行性,提出了CUDA环境下的GPU并行策略,实现了三维电阻率反演的快速化成像.其次,开展了基于GPU混合反演的数值算例,与传统最小二乘线性反演进行了对比,基于GPU并行计算的混合反演计算效率得到了显著提高,对含水构造的位置、形态有较好的反映,压制了三维反演的多解性.最后开展了物理模型试验,结果表明基于GPU混合反演探测的低阻异常体与实际含水构造的位置较为相符,发现基于GPU加速的混合反演方法在提高探测精度与加快反演速度方面具有显著优势,为三维电阻率混合反演方法在隧道超前探测实际工程中的应用奠定了基础. 相似文献
2.
基于光滑约束的最小二乘法是三维电阻率反演的主要方法,但该方法在某些情况下存在着多解性较强的问题,且普遍耗时较长,严重制约了三维反演方法的推广与发展.为改善上述问题,将表征模型参数变化范围的不等式约束作为先验信息引入最小二乘线性反演方法中,有效地改善了反演结果的精度,降低了反演的多解性问题.为了解决耗时较长的问题,基于预条件共轭梯度(PCG)算法和Cholesky分解法的特点提出了一套优化三维电阻率反演计算效率的计算方案.在该方案中,Cholesky分解法被用来求解敏感度矩阵计算中的多个点源场的正演问题,Cholesky分解法只需对总体系数矩阵进行一次分解,然后对不同的右端向量进行回代即可.将预条件共轭梯度法引入到三维电阻率反演方程的求解中,将雅可比迭代中的对角阵作为预处理矩阵,其具有求逆方便、无需内存空间的特点,有效地加快了收敛速度.对合成数据以及实测数据的反演算例表明,借助不等式约束和反演效率优化方案,最小二乘反演方法可得到较为精确的反演结果,有效地提高了反演计算效率,具有良好的推广前景. 相似文献
3.
瞬变电磁反演存在高度的非线性特征,常用的最小二乘等线性反演方法往往对初始模型高度依赖,并且极易陷入局部最优解.本文基于观测数据与模拟数据的L1范数建立目标函数,采用模拟退火非线性全局最优化方法实现瞬变电磁一维反演.初始模型完全随机产生,通过指数函数退温机制模拟系统能量最小实现迭代,通过接收概率函数评价当前模型,实现局部最优解的跳出,最终实现全局最优化求解.通过数值算例发现,无论给定的反演层数等于还是大于设计模型,都可以获得较好的反演效果,因而可以在反演初始就设计较多的层数,实现反演模型的自动拟合;同时,利用含噪声数据反演进一步验证算法的稳定性.最后,对实测数据进行了反演测试,结果与钻孔编录基本一致,表明提出的基于L1范数的模拟退火反演可用于实测数据处理. 相似文献
4.
5.
采用非线性最小二乘法结合蒙特卡罗法,实现宽带高频电磁场椭圆极化率数据的精确反演,确定地下层状介质的真实电阻率和介电常数.反演结果表明,对于均匀半空间和二层介质模型,最小二乘法能够很好地实现反演,而对于三层或更多层的介质,首先利用蒙特卡罗法确定拟合初始模型,再进行最小二乘反演,能够避免收敛到局部极小值,提高了反演的稳定性.为了加速正演响应函数的计算和迭代的速度,采用高密度采样的线性滤波算法,大大加快了该精确反演方法的速度.针对如覆盖区地质填图和土壤调查等大面积确定地质体性质的应用,本文还给出了一种近似反演方法（相位矢量图法）,能够快速获取视电阻率和视介电常数,不仅可以为应用提供有用的基础信息,而且可作为精确反演方法的初始模型. 相似文献
6.
7.
8.
利用遗传算法进行不考虑近场校正的全场资料CSAMT反演研究.遗传算法属于全局最优化方法,具有对初始模型依赖小,不易陷入局部极值的优点,然而,当未知数较多时,多解性仍是该方法的瓶颈.为了减小多层反演的多解性,在反演中引入最小构造约束,针对CSAMT的遗传算法反演问题定义了最小构造目标函数,经过模型试验找到了其具体表达式,并找到了适合CSAMT资料反演的拉格朗日乘子的最佳取值μ=0.5,实现了基于遗传算法的CSAMT最小构造反演.利用H、A、K、Q和HKH、KHA模型对方法进行了数值试验,在无噪和加入10％噪声情况下,反演结果与模型一致;加入20％噪声后,反演仍取得良好结果,与理论模型基本吻合.将该方法用于水平层状地层和横向变化地层的实测资料反演,结果与地质资料吻合.不同的计算实例表明了该方法的有效性. 相似文献
9.
混合范数下的最优化反演方法 总被引：4，自引：1，他引：4
在求解地球物理反问题时，通常根据最小二乘准则构造目标函数进行反演，并在实践中得到了广泛的应用.为进一步增强反演的稳健性及减少多解性，不损失反演结果的分辨率，本文提出了混合范数下的最优化反演方法，它根据数据和模型可能服从不同的概率分布，对数据空间和模型空间采用不同的范数来构造目标函数.在给出目标函数的基础上，导出了混合范数下的线性反演方程.由于该线性反演方程的复杂性，我们采用混合范数下迭代再加权共轭梯度法进行求解.最后，通过对模拟的电阻率数据进行反演，验证了本文计算方法是可行的. 相似文献
10.
基于包体岩石物理模型的储层参数地震反演方法面临数学形式复杂、多解性强、适应性差、涉及迭代运算等问题,本文提出一种基于约束最小二乘与信赖域的储层参数地震反演方法.该方法基于储层参数与弹性参数关联岩石物理模型,使用最小二乘方法构建目标函数和信赖域约束全局寻优求解,有效降低了地震反演多解性,极大提高了收敛速度.特别是通过在最小二乘求解中引入垂向约束,有效提高反演结果的抗噪声能力.经过模型测试和实际资料的应用,验证了方法的可行性和适用性. 相似文献
11.
直流电阻率法以其对含水构造响应敏感且适应性强的特点,是隧道施工期超前地质预报的常用方法.针对目前定点源三极观测方式易受到测线附近旁侧异常干扰与聚焦观测方式探测距离过短无法三维成像的问题,提出了一种基于多同性源阵列电阻率的隧道超前探测新方法,利用环形布置的多同性源供电压制隧道掌子面后方的异常体的干扰,同时增强掌子面前方远处异常体的敏感度,阵列电极测量可为含水构造的三维反演成像提供三维数据.首先,开展了含水构造超前探测的数值正演模拟,发现多同性源方法的视电阻率纵向微分曲线极小值以及与横轴的交点位置可对异常体位置进行定性判断,探讨了多同性源方法对后方及侧向异常干扰的响应特征及识别方法.其次,将电阻率光滑约束与轨迹光滑策略引入蚁群反演算法,提出了一种最小二乘方法与改进蚁群算法相结合的混合反演算法,反演数值算例表明混合反演算法能够实现含水构造的三维成像,较好地反演出不同距离的含水构造位置,并使其边界较为清晰.为了验证多同性源探测方法的探测效果开展了物理模型试验,多同性源阵列电阻率法能较好的反映含水构造的存在,与实际含水构造的位置较为相符. 相似文献
12.
13.
针对核磁共振测井双TW观测数据分析和流体识别的需要,研究了基于全局搜索的遗传算法和局部搜索的最小二乘法的联合反演算法,实现了核磁共振双TW观测数据处理.首先,研究了饱和油气水岩石物理模型的核磁共振双TW观测模式的测井响应机理;然后,利用全局搜索性能优良的遗传算法, 对核磁共振回波差数据进行了反演,计算出了流体的核磁弛豫性质及其体积;最后,以遗传算法的反演结果为初值, 利用阻尼最小二乘方法对双TW回波串进行更精细的反演,计算出了双TW的T2分布、孔隙度和流体饱和度.理想模型的合成数据和实际测井资料应用表明,遗传算法与最小二乘方法相结合是一种行之有效的联合反演方法,能较好地实现核磁共振测井双TW观测数据的处理和流体评价. 相似文献
14.
基于交叉梯度结构约束的不同数据的联合反演可以提高地球物理成像的可靠度,但是由于不同观测数据对地下介质的灵敏度不同以及不同反演系统灵敏度矩阵元素的数值可能存在较大的差异,导致构建和求解联合反演系统存在很大的挑战.针对以上问题,本文提出一种新的基于单独反演模型更新量与交叉梯度结构约束相结合的联合反演策略.该策略利用单独反演系统分别确定出各个模型的更新量,然后利用它们约束交叉梯度系统的反演,得出新的模型更新量.通过这样的方式,有效实现了数据拟合与结构约束的平衡,实现了基于交叉梯度联合反演的目标.新的联合反演策略不需要对原来单独反演程序修改即可实现联合反演,减小了联合反演实现的难度,极大地提高了联合反演的易实现性,而且避免了联合反演矩阵存储及结构过于复杂难以求解的问题.基于新的策略,本文首次实现了基于交叉梯度结构约束的三维地震走时与直流电阻率联合反演.合成模型测试表明,与单一成像相比,联合成像减少了地震走时反演中出现的干扰异常并提高了电阻率反演的分辨率. 相似文献
15.
在已有研究成果的基础上,为了适应物性参数剧烈变化的复杂模型并满足联合反演的要求,开发了速度和电阻率随机分布共网格单元模型的建模技术.基于这种统一的物性随机分布的网格介质模型,利用有限元方法和改进的射线追踪法分别正演计算大地电磁场和地震走时,结合改进的模拟退火算法,研究实现了电阻率和速度随机分布条件下的大地电磁与地震资料的同步联合反演.对物性界面不完全一致和物性变化剧烈的带地形复杂模型的试验,表明了该方法在精细反演复杂电阻率和速度结构方面的效果,克服了以往研究局限于简单模型的不足.对地震资料品质差的地区开展的实际资料联合反演,表明了方法的适用性,先验信息约束下的联合反演提高了反演精度. 相似文献
16.
径向基神经网络（RBFNN）具有结构简单、学习速度快、不易陷入局部极小等优点,能够有效地提高电阻率层析成像反演的收敛速度和求解质量.本文针对电阻率层析成像反演的非线性特征,提出了一种基于汉南-奎因信息准则（HQC）的正交最小二乘法（OLS）学习算法（HQOLS）.该算法通过计算HQC的最优值来自动选择RBFNN的网络结构,避免了传统OLS学习算法中阈值参数的设定,保证了网络的泛化性能.通过比较聚类法、梯度法、OLS和HQOLS等学习算法的反演性能,构建了基于RBFNN的电阻率层析成像反演模型.数值仿真和模型反演的结果表明,该方法实现简单,在准确性上优于BP反演,成像质量优于传统最小二乘法反演. 相似文献
17.
Introduction The research on the structure and physical property of ancient hidden hill, igneous rocks and basement is relatively difficult by using seismic data only. If we combine seismic data, magneto-telluric (MT) data and geophysical data together, better results can be obtained for the above problem. A number of geophysicists at home and abroad, such as CHEN and WANG (1990), Siri-punvarapor and Egbert (2000) have tried many methods to solve the problem by the inversion of seismic da… 相似文献
18.
地磁测深研究的周期范围通常为10~5~10~7s,缺少反映浅部电性结构的短周期信息,而C-响应受浅部电阻率影响明显,因此本文提出在反演中增加浅部(约200 km)电阻率约束以提高深部反演的稳定性和可靠性.在磁层环状电流满足P_1~0假设的条件下,球坐标系中一维导电薄球层状地球的C-响应和电导率分布关系由边界条件通过递推的方法计算得到.反演采用有限内存拟牛顿(L-BFGS)法;浅部电阻率约束通过将目标函数对模型参数的梯度设为零来实现;通过置信区间分析评价约束反演结果的可靠性.合成数据的无约束反演虽然最终的拟合效果很好,但浅部电阻率受初始模型影响,差异较大;采用浅部约束后,反演结果对初始模型依赖性明显减小,同时还能显著提高200~600 km范围内反演结果的准确性.对全球近地轨道卫星观测的C-响应数据约束反演后结果与前人一致,表现为地幔电导率整体上随着深度的增加而增加.参数置信区间分析表明,由于约束反演加入了浅部信息,电阻率的变化范围更加紧致,说明反演结果更加可靠.因此,有必要通过其他地球物理方法,如长周期大地电磁测深等获得浅部电阻率分布,作为先验信息参加反演,进行浅部约束的C-响应反演,获得更可靠的一维全深度电性结构,为地磁测深数据解释奠定基础. 相似文献