首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 265 毫秒
1.
弹性波场的复杂性使得反演问题非线性增强,容易陷入局部极值,需要采用合理的多尺度反演策略降低非线性。在逐频组多尺度反演的基础上引入第二级别的反演策略,即基于阻尼波场的层剥离方法,能够改善反演过程的稳定性。针对频域全波形反演计算效率低、内存占用大的问题,采用双级并行算法:(i)利用多波前大规模并行直接解法(MUMPS)软件包,多节点并行实现波场正演;(ii)基于MPI实现频组内各频率并行计算梯度和步长等,使得多尺度反演算法在提高精度的前提下,保证了计算效率,提高了算法的实用性。Overthrust模型的数值实验表明,本文反演算法能够在有效改善反演稳定性的前提下,高效地获得精度较高的反演结果。  相似文献   

2.
大地电磁三维数据空间反演并行算法研究   总被引:6,自引:4,他引:2  
目前大地电磁三维反演实际应用的主要问题是计算效率低.在对大地电磁三维数据空间反演算法进行深入分析的基础上,本文提出了基于频点和矩阵划分的大粒度并行反演方案和具体实现步骤,并在曙光TC5000A高性能计算平台上实现了基于MPI的大地电磁三维数据空间反演并行算法.该算法实现了包括三维正演、灵敏度矩阵、叉积矩阵以及模型改正量的并行执行,不仅计算效率高,而且每个节点机上灵敏度矩阵的存储空间只需原来微机上的2/N(N是参加并行计算的节点机个数),大大地减少了内存开销.通过两个理论模型合成的数据对实现的三维数据空间反演并行算法进行试算,对比分析了多个节点机下程序的执行效率.测试结果表明,所实现的三维数据空间反演并行算法是可行的、高效的,与单机相比,不仅可以提高运行速度,缩短计算时间,而且还可以扩大计算规模,极大地推动了大地电磁三维反演的实用化.  相似文献   

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

4.
隧道施工期超前探测对于避免突涌水灾害的发生具有重要作用,为满足隧道三维电阻率超前探测快速化解译与成像的要求,本文提出了一种基于GPU并行的蚁群算法与最小二乘方法相结合的混合反演算法.该方法结合线性反演与非线性反演的优点,利用蚁群算法全局搜索能力强的优点为最小二乘反演提供较优的初始模型,以克服最小二乘算法容易陷入局部最优的缺点,提高了隧道三维电阻率反演成像的精度.同时,基于蚁群算法的天然并行性,提出了CUDA环境下的GPU并行策略,实现了三维电阻率反演的快速化成像.其次,开展了基于GPU混合反演的数值算例,与传统最小二乘线性反演进行了对比,基于GPU并行计算的混合反演计算效率得到了显著提高,对含水构造的位置、形态有较好的反映,压制了三维反演的多解性.最后开展了物理模型试验,结果表明基于GPU混合反演探测的低阻异常体与实际含水构造的位置较为相符,发现基于GPU加速的混合反演方法在提高探测精度与加快反演速度方面具有显著优势,为三维电阻率混合反演方法在隧道超前探测实际工程中的应用奠定了基础.  相似文献   

5.
基于FFT-MA谱模拟的快速随机反演方法研究   总被引:3,自引:2,他引:1       下载免费PDF全文
虽然基于地质统计学的随机反演方法能够有效融合测井资料中的高频信息,但计算效率低,占用内存大,限制了它在实际资料中的应用.本文在保留传统随机反演方法优点的基础上,创造性地引入傅里叶滑动平均(Fast Fourier Transform-Moving Average,FFT-MA)谱模拟进行频率域的地质统计模拟,并利用逐步变形算法(Gradual Deformation Method,GDM)确保模拟结果与实际地震数据的匹配,构建了基于FFT-MA谱模拟的新的快速随机反演方法.与常规随机反演相比,新方法不仅分辨率高,而且能够使反演解得到快速收敛,有效提高计算效率,减少内存占用.模型试算获得了与理论模型吻合度较好的高分辨率反演结果.实际资料分析也表明新方法所得到的高分辨率反演结果能够对薄互储层进行良好的展示,为薄储层的识别提供高效可靠的技术支持.  相似文献   

6.
面波频散反演地球内部构造的遗传算法   总被引:41,自引:16,他引:41       下载免费PDF全文
介绍了一种新的算法--遗传算法的基本概念和特点,及其在地震面波反演地球内部构造中的应用,指出了使用遗传算法的注意事项.提出了通过对初步搜索结果参数分布直方图进行分析,从而修改和缩小进一步搜索的范围,逐步搜索以提高搜索效率的方法.并对3层含低速层的理论模型和青藏高原的实际频散资料进行遗传算法反演,获得了满意的结果.讨论了遗传算法在其他地震学问题中进一步应用的可能性.  相似文献   

7.
提出了一种新的非导数搜索方法,用于在多维参数空间中寻找可较好地拟合观测数据的模型,该方法与模拟退火和遗传算法同属一类,一般情况下,模拟退火和遗传算法是用于解决全局最优问题的。这里的目是寻找参数空间中对数据拟合较好域进行优先采样的模型集合,而不是寻找单一的最优模型(相关的文章涉及集合的定量评价)。新的搜索算法利用Voronoi单元的几何建造在参数空间引导搜索,这些单元是根据适当的距离标准定义的最临近区域,算法的概念是简单的,仅需两个调节参数,并且只使用数据的秩拟合准则,而不是数据值,这样避免了与数据失配函数缩放比例相关的一些困难,并能使用任意组合的数据拟合标准,通过不时地使用最相邻的计算代替正演模拟计算,也显示出了Voronoi单元是如何能增强任何现有直接搜索方法的性能的。通过对理论模型的地壳地震结构接收函数反演的应用,对新的直接搜索算法进行了解释,众所周知,这是一个非线性问题,线性化反演技术强烈地依赖于初始解,新的算法生成复杂的自适应搜索行为,据我们所知以前的算法不具有这种功能。  相似文献   

8.
大地电磁测深法是以岩石的电性差异为基础和前提的勘探方法.本文所采用的大地电磁二维反演方法为共轭梯度法,该方法避免了求解雅可比矩阵,效率较高,但是在将模型剖分为较细网格,多频率进行计算时效率有待提高.基于大地电磁根据各频率依次独立处理数据的特点,在本文中采用了MPI并行运算方法,用多个进程同时来计算各频率数据,最后再将数据进行收集,得到最后的计算结果.通过对正演和反演结果的比较,验证了程序的正确性.对并行算法的效率进行了统计,进程数为2~8时,加速比能达到1.63~2.64,验证了并行算法的有效性.  相似文献   

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

10.
基于逆算子估计的AVO反演方法研究   总被引:1,自引:1,他引:0       下载免费PDF全文
传统反演算法以优化算法为主,而基于逆算子估计的AVO反演算法则利用了直接求逆的思路.算法的关键在于寻找存在逆函数的子域,进而可以在子域内直接求逆,这种解决反问题的思路不同于一般的优化类算法所采用的直接搜索解的方式,具有更高的效率.AVO反演利用了振幅随着偏移距的变化特征,反演的精度受到地震资料质量的影响,通过加入L1范数约束以及合理的初始模型有助于提高反演的稳定性以及准确度.模型测算和实际应用表明,基于逆算子估计的AVO反演方法具有较高的精确程度和可靠性.  相似文献   

11.
大地电磁三维快速松弛反演并行算法研究(英文)   总被引:8,自引:1,他引:7  
为了进一步提高大地电磁三维快速松弛反演的计算效率,在深入分析大地电磁三维快速松弛反演算法的基础上,结合MPI自身的优越性,确定了并行计算的思路,实现了三维快速松弛反演的并行计算。通过理论模型合成数据和实测数据对实现的三维快速松弛反演并行程序进行了试算,分析对比了在多种情况下程序的执行效率。测试结果表明,所实现的三维快速松弛反演并行程序运行结果正确,效率提高明显。此成果在普通微机上实现,推动了大地电磁三维反演技术的实用化,可为其它地球物理三维正反演研究所借鉴。  相似文献   

12.
超高密度电法是一种新的地球物理探测技术,它通过多通道数据采集和多装置数据联合反演,极大地提高了电法勘探的成像精度.本文提出一种主成分-正则化极限学习机(PC-RELM)非线性反演方法,该方法针对超高密度电法所获取的高维勘探数据进行反演建模,通过随机设定隐层参数来简化模型的学习过程,通过主成分分析方法来进行高维数据降维,最后引入正则化因子提高反演模型的泛化能力.论文给出了超高密度电法的原理、样本构造方法和非线性反演流程,使用交叉验证方法获得了优化的隐节点数目和正则化参数,构造了优化的反演模型.通过两个经典的超高密度模型的反演结果表明,该方法能够较好地解决超高密度电法反演的高维数据非线性建模问题,能够弥补单一装置数据反演的不足,同时相较其他的非线性反演方法(ELM,BPNN和GRNN)具有更加准确的反演结果.  相似文献   

13.
基于改进粒子群算法的地震标量波方程反演   总被引:4,自引:2,他引:2       下载免费PDF全文
针对标准粒子群优化(PSO)算法存在易出现早熟而陷入局部最优以及进化后期收敛速度慢等缺陷,通过考虑粒子所处位置间相互作用,提出了一种改进的并行粒子群优化算法.由于引入粒子位置间的相互影响,减少了粒子搜索过程盲目性,因此能有效提高算法的收敛速度.数值试验表明,这种改进的粒子群算法适用于二维标量波方程的速度反演,且算法具有...  相似文献   

14.
大地电磁三维交错网格有限差分数值模拟的并行计算研究   总被引:3,自引:3,他引:0  
为了更有效的提高大地电磁三维正演的计算速度,引入了并行处理技术.大地电磁三维交错网格有限差分数值模拟是按照不同频率来计算的,各频率之间求取电磁场值的过程是相互独立的.根据这一特点,可以将多个频率的计算任务平均划分为一个或者几个频率的计算子任务,分配到各个计算节点去并行执行,计算完成后将结果汇总.本文通过采用主从并行模式、分频并行计算的方案,在曙光TC5000A高性能并行平台上实现了基于MPI的大地电磁三维正演的并行计算.通过两个理论模型对实现的大地电磁三维正演并行算法进行试算,对比分析了多个节点机下程序的执行效率.测试结果表明,所实现的三维正演并行算法是正确的、高效的,为进一步的大地电磁三维反演并行算法研究奠定了重要基础.  相似文献   

15.
用斜井钻孔崩落资料求取地壳应力状态是一个非线性反演问题。而遗传算法是解反演问题的有效方法之一。本文由斜井钻孔崩落的理论模型,用人工数据检验了遗传算法应用于此问题的有效性。对冀中坳陷内的饶阳凹陷和深县凹陷两个地区的6口油井的钻孔崩落数据进行了处理,用遗传算法反演了这两个凹陷的上地壳应力状态,与用非线性反演和枚举相结合的方法相比,大大减少了计算量,得到了令人满意的结果。  相似文献   

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

17.
In this work a new algorithm for the fast and efficient 3D inversion of conventional 2D surface electrical resistivity tomography lines is presented. The proposed approach lies on the assumption that for every surface measurement there is a large number of 3D parameters with very small absolute Jacobian matrix values, which can be excluded in advance from the Jacobian matrix calculation, as they do not contribute significant information in the inversion procedure. A sensitivity analysis for both homogeneous and inhomogeneous earth models showed that each measurement has a specific region of influence, which can be limited to parameters in a critical rectangular prism volume. Application of the proposed algorithm accelerated almost three times the Jacobian (sensitivity) matrix calculation for the data sets tested in this work. Moreover, application of the least squares regression iterative inversion technique, resulted in a new 3D resistivity inversion algorithm more than 2.7 times faster and with computer memory requirements less than half compared to the original algorithm. The efficiency and accuracy of the algorithm was verified using synthetic models representing typical archaeological structures, as well as field data collected from two archaeological sites in Greece, employing different electrode configurations. The applicability of the presented approach is demonstrated for archaeological investigations and the basic idea of the proposed algorithm can be easily extended for the inversion of other geophysical data.  相似文献   

18.
在利用地震波数据进行地球物理反演时,地层对地震波的吸收衰减效应会对地层物性参数的准确反演产生较大的影响,因此利用黏弹性声波方程进行反演更符合实际情形.本文在考虑地层衰减效应进行频率空间域正演模拟的基础上,提出基于黏弹性声波方程的频率域逆散射反演算法并对地震波传播速度进行反演重建,在反演过程中分别用地震波传播复速度和实速度来表征是否考虑地层吸收衰减效应.基于反演参数总变差的正则化处理使反演更加稳定,在反演中将低频反演速度模型作为高频反演的背景模型进行逐频反演,由于单频反演过程中背景模型保持不变,故该方法不需要在每次迭代中重新构造正演算子,具有较高的反演效率;此外本文在反演过程中采用了基于MPI的并行计算策略,进一步提高了反演计算的效率.在二维算例中分别对是否考虑地层吸收衰减效应进行了地震波速度反演,反演结果表明考虑衰减效应可以得到与真实模型更加接近的速度分布结果,相反则无法得到正确的地震波速度重建结果.本文算法对复杂地质模型中浅层可以反演得到分辨率较高的速度模型,为其他地震数据处理提供比较准确的速度信息,在地层深部由于地震波能量衰减导致反演分辨率不太理想.  相似文献   

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

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