首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 171 毫秒
1.
改进的模拟退火-单纯形综合反演方法   总被引:19,自引:6,他引:13       下载免费PDF全文
实际中的大量地球物理反演是一个多参数、非线性优化问题,所采用的目标函数,即度量由参数化的理论模型得出的预测值与观测值的吻合程度,往往具有多个局部极值.针对这类问题,本文综合全局反演方法具有的全域搜索能力强、局部方法收敛速度快和“均匀设计”布点效率高的特点,提出了模拟退火-单纯性综合反演方法,并通过一维声波非线性反演验证了这种综合方法的搜索能力和效率.  相似文献   

2.
磁法反演属于非线性最优化问题,具有多变量、目标函数多极值、反演多解性等特点,因此,需要稳定的和高效的优化反演算法.粒子群优化已开始被用于地球物理反演计算,但是对于高维数、多峰值函数,粒子群的收敛精度不高,容易陷入局部极值.如果将混沌局部搜索和粒子群优化的优势相结合,通过将种群搜索过程对应为混沌轨道的遍历过程,可使标准粒子群优化的搜索过程具有避免陷入局部极小的能力.本文利用混沌-粒子群优化用于磁法反演计算.数据试验结果表明,该方法可以用于磁法数据的地球物理非线性反演,并且在一定程度上优于标准粒子群优化方法.  相似文献   

3.
从CSAMT信号中提取激电信息有利于提高频率域电磁法反演与解释的精度.目前的研究多以线性反演方法为主,存在依赖初始模型、易陷入局部极值的问题.针对CSAMT信号IP提取问题的非线性和非凸特征,本文提出了一种基于柯西分布和惯性权重的二阶段最小构造混合蛙跳反演方法来提取IP信息.该方法首先利用柯西算子取代随机算子来提高算法的全局搜索能力,并通过引入混沌震荡惯性权重来均衡进化过程中的个体经验和群体经验,保证算法后期的稳定收敛;然后通过引入第二阶段反演过程来强化极化率对观测数据的影响,同时将正则化参数引入混合蛙跳算法的适应度函数来改善反演的多解性问题;最后利用CPU并行计算加速了算法的模因组搜索过程.反演结果表明,上述方法能够较好地重构地电结构和提取激电信息,在加噪环境下具有较强的鲁棒性.相比其他非线性算法(标准混合蛙跳算法SFLA,差分进化算法DE和粒子群优化算法PSO)的反演结果,本文算法具有更强的全局搜索能力和更高的计算效率,适合对微弱的激电信息进行提取.  相似文献   

4.
基于不等式约束的三维电阻率探测混合反演方法   总被引:2,自引:2,他引:0       下载免费PDF全文
三维电阻率探测的线性反演和非线性反演中均存在着多解性的固有难题.电阻率线性反演方法的效率较高,但反演结果对初始模型的依赖性较强,易陷入局部极小;而非线性反演方法不依赖初始模型,但搜索效率极低,尚未见到关于三维电阻率非线性反演的文献.针对上述问题,融合线性与非线性反演方法的互补优势,提出了最小二乘法(线性方法)与改进遗传算法(非线性方法)相结合的混合反演方法的概念和思想.首先,提出了将介质电阻率变化范围作为不等式约束引入反演方程的思路,以实现压制多解性、提高可靠性的目标.提出了宽松不等式约束和基于钻孔推断的局部严格不等式约束的获取及定义方法.在此基础上,分别提出了基于不等式约束的最小二乘线性反演方法和遗传算法非线性反演方法.其次,对于遗传算法在变异搜索方向控制、初始群体产生等方面进行了改进,优化了其搜索方向和初始群体多样性.然后,提出了混合反演方法及其实现方案,利用改进遗传算法进行第一阶段反演,发挥其对初始模型的依赖程度低的优势,搜索到最优解附近的空间,输出当前最优个体;利用最小二乘法进行第二阶段反演,将遗传算法得到的当前最优个体作为初始模型,在最优解附近空间执行高效率的局部线性搜索,最终实现地电结构的三维成像.最后,开展了合成数据与实际工程算例验证,与传统最小二乘方法进行了对比,发现混合反演方法在压制多解性、摆脱初始模型依赖和提高反演效果方面有较好效果.  相似文献   

5.
京津唐张地区速度结构和震源位置联合反演的遗传算法   总被引:12,自引:0,他引:12  
震源参数和速度结构的联合反演是一个典型的非线性名参数最优化问题,常规的局部线性化反演方法往往易于陷入局部极值,且严重依赖于初始模型的选取。模拟生物界进化的遗传算法则是一种简单而高效的全局性搜索方法。  相似文献   

6.
自适应全局混合反演   总被引:27,自引:0,他引:27  
地球物理学反演多是一个多参数、非线性、多极值、非连续函数的优化问题 .在模拟退火与单纯形、均匀设计中引入自适应的退火规则 ,形成自适应全局混合反演方法 .数值检验与模型计算表明 ,该方法具有搜索速度快、效率高的特点  相似文献   

7.
反演瑞雷波频散曲线能有效获取地层横波速度和厚度.但由于其高度的非线性、多参数、多极值等特点,传统的全局搜索方法易出现收敛速度慢、早熟收敛及搜索精度低的问题.鉴于此,本文提出并测试了基于萤火虫优化算法(FA)和带惯性权重的蝙蝠优化算法(WBA)的新的瑞雷波频散曲线反演策略.在瑞雷波频散曲线反演中,FA全局搜索能力强,但后期搜索精度低,而WBA局部搜索能力强,搜索精度高,但易出现早熟收敛.故本文将二者结合,提出了一种新的优化策略,称其为WFBA,即在反演前期使用FA,后期使用WBA,很好地解决了FA后期搜索精度低及WBA早熟收敛的问题.本文首先反演了三个典型理论模型的无噪声、含噪声的数据,验证了WFBA对瑞雷波数据反演的有效性与稳定性.然后将WFBA与WBA、FA单独反演以及不含惯性权重的FBA和粒子群优化算法(PSO)反演的结果进行了对比,说明了WFBA相对于WBA、FA、FBA和PSO具有更稳定、收敛速度更快、求解精度更高等优点.最后,反演了来自美国怀俄明地区的实测资料,检验了WFBA对瑞雷波数据反演的实用性.理论模型试算和实测资料分析表明,WFBA很适用于瑞雷波频散曲线的定量解释,具有很高的实用性价值.  相似文献   

8.
一维波动方程波阻抗反演的同伦方法   总被引:10,自引:3,他引:7       下载免费PDF全文
文中从地震勘探一维波动方程反问题出发,研究了一种反演地层参数的同伦方法,该方法把非线性方程组的求解转化成常微分方程初值问题的数值求解,从而给出一种稳定的计算速度快、抗噪能力强的全局收敛的反演方法.理论模型和实例试算的结果表明了同伦方法是一种有效的反演算法,特别适用于非线性、多极值的地球物理反演问题,在地球物理非线性反演中具有广泛的应用前景.  相似文献   

9.
瑞雷面波技术是浅地层地质勘探的重要手段之一,通过瑞雷面波频散曲线反演可有效获得地下横波速度模型.然而,瑞雷面波频散曲线反演具有多参数、多极值和非线性等特点,面对复杂地震-地质条件下的瑞雷面波资料处理,传统方法难以快速、精确地对地层参数进行反演和重建.本文将自适应对数螺旋路径萤火虫算法(Adaptive Logarithmic Spiral-Lévy Firefly Algorithm,简称ALSL-FA)引入到瑞雷面波频散曲线反演中,有效解决了经典萤火虫算法(Firefly Algorithm,简称FA)精于探索,疏于开发的缺点.ALSL-FA集成了对数螺旋引导萤火虫路径,并在搜索过程中通过自适应切换因子实现了全局搜索与局部开发自适应切换,在增强局部开发能力的同时,确保了全局搜索的能力.通过Rastrigin函数求解,对比分析了FA、Lévy飞行萤火虫算法(Lévy Flying Firefly Algorithm,简称LF-FA)和ALSL-FA的运算性能;通过速度递增和含高速硬夹层地质模型的瑞雷面波频散曲线反演,比较了人工蜂群算法(Artificial Bee Colony alg...  相似文献   

10.
复杂地形条件下静校正的综合寻优   总被引:19,自引:2,他引:19       下载免费PDF全文
复杂地形,特别是山地条件下地震勘探中的大静校正问题是 一个非线性的、具有多参数多极值的全局优化难题. 本文综合利用最大能量法、模拟退火与 遗传算法的各自优势,提出一种具有局部收敛速度快、全局搜索能力强的综合寻优反演方法 . 其中对产生伪静校正量的零空间现象进行空间平滑滤波剔除. 理论数值计算及实际资料的 处理说明本方法适应能力强、静校正效果较好.  相似文献   

11.
研究了基于矢量有限元方法的大地电磁带地形三维反演算法并开发了三维反演计算程序代码.在大地电磁场正演数值模拟方面,采用并行直接稀疏求解器PARDISO且无需进行散度校正的快速正演方案,对典型地形模型,在中等规模计算条件下,与双共轭梯度法(BICG)计算结果比较,发现PARDISO比BICG快10倍以上;通过理论模型试算,并与前人的有限元法计算结果对比,验证了带地形三维正演计算程序的正确性.在反演方面,本研究基于共轭梯度方法编写了大地电磁带地形三维反演代码,为了避免直接求取雅可比矩阵,将反演中的雅可比矩阵计算问题转为求解两次“拟正演”问题,进而将PARDISO的快速正演方案应用于“拟正演”问题的求解,以提高反演计算效率.利用开发的反演算法对多个带地形地电模型的合成数据进行了三维反演,反演结果能很好地重现理论模型的电性结构,验证了本文开发的三维反演算法的正确性和可靠性.最后,利用该算法反演了某矿区大地电磁实测数据,反演得到的三维电性结构清晰地反映了研究区的地电特征,将反演结果与该区已有地质资料结合进行解释,应用效果明显,进一步验证了本文算法的有效性.  相似文献   

12.
井地有限线源三维电阻率反演研究   总被引:2,自引:1,他引:2       下载免费PDF全文
有限线源的三维电阻率反演,是地学探测研究的重要课题之一.本文主要讨论了以下几个方面的内容:(1)结合反演方程式改进的共轭梯度算法;(2)基于扰动法的线源Jacobi矩阵的近似计算方法及光滑系数矩阵在各个方向上的光滑因子的计算公式;(3)通过分析阻尼系数与修改量校正量对反演结果的影响,提出了利用校正系数对模型修正量进行校正;(4)较为系统的讨论了阻尼系数λ对反演分辨能力的影响,较好的改善了三维电阻率反演中的电性异常体重心的"上漂"现象.数值实验表明,改进后的共轭梯度法反演对初始模型的依赖程度较小,反演能够稳定收敛,对模型的空间位置分辨率较高.  相似文献   

13.
基于阵列感应与自然电位联合反演地层水电阻率   总被引:2,自引:2,他引:0       下载免费PDF全文
原状地层水电阻率是重要的储层参数,也是进行精细储层评价的基础.基于泥浆侵入数值模拟与侵入过程中井周岩石物理特征分析,确定了薄膜电位的产生位置,针对储层高、低侵等不同侵入特征,提出了可适用于包括存在"低阻环"等不同侵入特征时储层电阻率分布的"五参数"地层模型,基于几何因子理论与有限元方法,建立了阵列感应与自然电位测井联合反演算法,实现了地层电阻率参数反演,重构了地层径向电阻率剖面,进而精确求取了地层水电阻率.通过对实际资料处理表明:反演算法稳定可靠,"五参数"模型能很好地表征储层电阻率分布形态,重构储层电阻率剖面,确定薄膜电位产生位置;基于阵列感应与自然电位的联合反演,能精确计算原状地层水电阻率,为储层评价与流体性质识别提供依据.  相似文献   

14.
可控源音频大地电磁三维共轭梯度反演研究   总被引:9,自引:5,他引:4       下载免费PDF全文
可控源音频大地电磁法在资源勘探等领域中发挥着重要的作用.我们把有限差分数值模拟方法用于可控源音频大地电磁三维正演,结合正则化反演方案和共轭梯度反演的思路,将反演中的雅可比矩阵计算问题转为求解两次"拟正演"问题,得到模型参数的更新步长,形成反演迭代,实现了可控源音频大地电磁三维共轭梯度反演算法.该反演算法可用于对有限长度电偶源激发下采集到的可控源音频大地电磁全区(近区、过渡区和远区)视电阻率和相位资料进行三维反演定量解释,获得地下三维模型的电阻率结构.理论模型合成数据的反演算例验证了所实现的可控源音频大地电磁三维共轭梯度反演算法的有效性和稳定性.  相似文献   

15.
大地电磁的人工鱼群最优化约束反演   总被引:3,自引:2,他引:1       下载免费PDF全文
大地电磁的反演问题是非线性,如果采用线性反演方法容易陷入局部极小,使得反演结果非唯一性严重.本文将人工鱼群算法引入到地球物理反演之中,提出了非线性的大地电磁人工鱼群最优化反演.该方法不需要进行偏导数的求取,可以对反演的范围进行约束,以减小反演结果的非唯一性.同时我们对搜索步长进行了改进,给出适用于大地电磁反演的人工鱼群参数.大量的理论数据试算表明,人工鱼群反演算法能够较好地寻找到全局最优解.实测数据的处理结果表明,该方法可以用来处理实际资料,并且能够取得很好的应用效果.  相似文献   

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

17.
叠前地质统计学反演将随机模拟与叠前反演相结合,不仅可以反演各种储层弹性参数,还提高了反演结果的分辨率.基于联合概率分布的直接序贯协模拟方法可以在原始数据域对数据进行模拟,不需要对数据进行高斯变换,拓展了地质统计学反演的应用范围;而联合概率分布的应用确保了反演参数之间相关性,提高了反演的精度.本文将基于联合概率分布的直接序贯协模拟方法与蒙特卡洛抽样算法相结合,参考全局随机反演策略,提出了基于蒙特卡洛优化算法的全局迭代地质统计学反演方法.为了提高反演的稳定性,我们修改了局部相关系数的计算公式,提出了一种新的基于目标函数的优化局部相关系数计算公式并应用到协模拟之中.模型测试及实际数据应用表明,该方法可以很好的应用于叠前反演之中.  相似文献   

18.
For density inversion of gravity anomaly data, once the inversion method is determined, the main factors affecting the inversion result are the inversion parameters and subdivision scheme. A set of reasonable inversion parameters and subdivision scheme can, not only improve the inversion process efficiency, but also ensure inversion result accuracy. The gravity inversion method based on correlation searching and the golden section algorithm is an effective potential field inversion method. It can be used to invert 2D and 3D physical properties with potential data observed on flat or rough surfaces. In this paper, we introduce in detail the density inversion principles based on correlation searching and the golden section algorithm. Considering that the gold section algorithm is not globally optimized, we present a heuristic method to ensure the inversion result is globally optimized. With a series of model tests, we systematically compare and analyze the inversion result efficiency and accuracy with different parameters. Based on the model test results, we conclude the selection principles for each inversion parameter with which the inversion accuracy can be obviously improved.  相似文献   

19.
二维地震资料波动方程非线性反演   总被引:3,自引:3,他引:3       下载免费PDF全文
针对反演的要求和实际问题的需要,提出利用地震资料叠前数据进行二维波动方程反演,采用最小平方拟合修正模型参数的非线性反演方法,构造了问题的加速迭代算法.反演算法充分利用了冗余的叠前数据和多道相关性,可以分离噪声和信号,使噪声不参与或很少参与反演,算法抗噪能力强.数值模拟例子表明算法有效和稳定,得到了令人满意的结果.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号