首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 140 毫秒
1.
二维波动方程速度的正则化-同伦-测井约束反演   总被引:17,自引:4,他引:13       下载免费PDF全文
傅红笋  韩波 《地球物理学报》2005,48(6):1441-1448
针对二维波动方程反问题,将大范围收敛的同伦方法引入速度参数的反演过程中,并将其与求解不适定问题的Tikhonov正则化有机结合,提出了一种新的、特别适用于非线性的、不适定的、多极值的地震勘探反演问题的反演策略:正则化-同伦方法. 为了充分利用测井资料和地震资料的互补特征,进一步提高反演分辨率并压制噪声,设计了正则化-同伦-测井约束联合反演方法. 大量数值试验结果表明了这两种方法的有效性.  相似文献   

2.
用地质雷达数据资料反演二维地下介质的方法   总被引:7,自引:3,他引:7       下载免费PDF全文
从二维麦克斯韦方程组出发推导出反演介电常数和电导率等二维介质物性参数的反演公式.反演的步骤是: 建立初始猜测模型,利用电磁波时间域有限差分法模拟正演数据,用正演数据与观测数据之间的数据残差建立目标函数,通过引入一个由麦克斯韦方程计算的伴随场,将目标函数对介质参数的导数表示成显式形式,应用最优化理论得出对初始猜测模型的修改,用共轭梯度法迭代,最终得到反演结果.用合成数据反演具有粗糙地表的非导电介质的介电常数,用实验数据同时反演介电常数和电导率,并比较了麦克斯韦方程反演结果与声波方程反演结果、波动方程偏移剖面的差异.  相似文献   

3.
二维波动方程参数反演的微分进化算法   总被引:2,自引:1,他引:1       下载免费PDF全文
利用微分进化算法进行二维波动方程参数反演的研究.微分进化算法属于全局最优化算法,具有鲁棒性强、收敛速度快、计算精度高的优点.把参数反演问题转化为非线性优化问题,通过保持在解空间不同区域中各个点的搜索,以最大的概率找到问题的全局最优解.同时还进行了数值模拟,结果令人满意.  相似文献   

4.
井间电磁场的一维、二维联合反演方法   总被引:17,自引:0,他引:17       下载免费PDF全文
提出了针对轴对称二维井间模型的一维、二维联合反演方法.该方法基于纵向成层背景地层的Green函数,以一维反演结果作为联合反演的迭代初始值,同时反演一维背景地层和二维异常剖面的电导率.采用递推技术计算Green函数的系数,可以很方便地同时得到Green函数对空间坐标的积分和对背景层电导率的微分,大大加快了计算雅可比矩阵〖WTHX〗M〖WTBZ〗的速度并使计算方便、准确.反演实例说明了该方法的有效性.反演实例还显示,背景层电导率较二维成像剖面的电导率收敛速度要快,测量数据的精度对二维剖面电导率成像质量的影响大于对背景层电导率分布的影响.  相似文献   

5.
高分辨率的频率空间域声波全波形速度反演-理论模型   总被引:1,自引:1,他引:0  
使用最速下降法进行二维频率空间域声波波动方程全波形速度反演,讨论了如何快速实现高精度的二维频率空间域声波波动方程全波形速度反演.多尺度的思想耦合在反演框架中.把非线性问题化为逐步线性问题是我们关注的焦点,目的是把整个非线性反演的黑匣子转化成为每一步可控的过程,尽可能得到想要的反演解.仅仅使用3个离散的频率,每个频率迭代...  相似文献   

6.
基于声波方程的井间地震数据快速WTW反演方法   总被引:4,自引:2,他引:4       下载免费PDF全文
WTW(Wave equation traveltime+Waveform inversion)反演是基于波动方程的走时反演(WT反演)和波形反演的联合反演方法.WT反演利用波动方程计算走时和走时关于速度的导数,和传统以射线为基础的走时反演相比,具有不必射线追踪、不必拾取初至、不必高频假设以及初始模型和实际模型差别较大时也能较好收敛等优点,但WT反演与波形反演相比其结果分辨率低.与之互补的是,波形反演的反演结果分辨率高,但是当所给初始模型和实际模型相差太大时,波形反演迭代算法容易陷入局部极小点.可见结合两种方法的WTW反演是一种比较好的联合反演方法.常规WTW迭代算法是首先以WT反演为主反演得到地质模型的整体特征,然后再以波形反演为主反演模型细节,该算法耗时和占用计算机存储空间接近WT反演或波形反演的两倍.为了节省运算耗时和计算机存储空间,往往采取首先单独利用WT反演然后再单独利用波形反演的算法.这样做的缺点是不能紧密结合两种反演方法,使得它们的优缺点在每一次迭代中无法得到互补,从而影响了最终的反演结果.针对以上事实,本文提出一种新的方法实现WTW,使得WTW运算速度和存储空间在任何情况下等同于WT反演或波形反演.模型计算表明新的算法具有更好的收敛性.  相似文献   

7.
黏弹性双相介质参数反演的同伦方法   总被引:4,自引:1,他引:3       下载免费PDF全文
本文以介质的动态时域响应(位移、速度、加速度)作为反演的依据,对黏弹性双相介质的材料参数进行了反分析研究.根据附加条件,并考虑到由于测量误差和测量数据的不充分造成的反演问题的不适定性,将黏弹性双相介质的参数反演问题转化为非线性泛函的极小值问题,然后应用大范围收敛的同伦方法求出使泛函极小的根作为反问题的解.最后,以二维半空间黏弹性双相介质模型为例进行了数值反演分析,数值结果表明本文方法在参数的反演过程中目标函数稳定收敛,且具有一定的抗“噪声”干扰性能;与利用速度和加速度信息反演的结果相比,利用位移信息反演的结果具有较高的精度和抗噪声能力.  相似文献   

8.
基于二维波动方程的物性参数联合反演方法,以常规地震勘探资料为基础、地震地质解释的测井资料为约束,采用线性、非线性反演、神经网络等方法进行地层物性参数预测。其方法保证了反演的成功进行和成果的可靠性,其高分辨率的岩性参数反演成果可有效地用于煤田岩性地震勘探。  相似文献   

9.
二维地电模型电(电磁)测深法参数化反演,宜采用多参数反演方案.文中导出了多参数反演的基本数学模型,并进行了分解处理,使之变成与剖面上反演点数相等的多个相互独立的小数学模型的组合,用广义逆法逐个求解,即可得一次迭代中整个二维模型参数的修改量.大地电磁测深理论模型和实测数据反演结果表明所提出的二维多参数反演方法可行.它具有稳定的收敛性、较快的收敛速度、能适应较复杂构造的反演和计算量小等特点.  相似文献   

10.
为了构建合理可靠的初始模型,使反演稳定快速收敛且达到更高的精度.本文提出用控制转换算法构建初始模型参数,然后使用最优分割法优化初始模型.一维MT最小二乘反演计算结果表明:该方法构建的初始模型使反演稳定快速收敛;使用相同的反演方法,该方法构建的初始模型的反演迭代次数比Bostick法构建的初始模型的反演迭代次数少20%以上,比均匀半空间作为初始模型的反演迭代次数少100%以上,并且反演精度也有了提高.该方法不但对提高MT一维反演效果有意义,而且对构建其他勘探方法的反演初始模型也具有一定的参考价值.  相似文献   

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

12.
地震波形反演的稀疏约束正则化方法   总被引:5,自引:5,他引:0       下载免费PDF全文
本文考虑地震波形反演问题.为了克服传统的Tikhonov正则化方法过度光滑的弊端,引入了非线性稀疏约束正则化方法,并采用对偶方法求解稀疏约束泛函的极小点.基于二维声波方程波形反演问题进行了数值模拟,针对不同模型对稀疏约束正则化方法进行了测试.结果表明,稀疏约束正则化方法对不连续介质模型的介质边缘具有良好的识别能力.  相似文献   

13.
用平面波延拓方程进行地震数据的叠前速度反演   总被引:1,自引:1,他引:1  
本文讨论地震勘探数据的叠前速度反演方法及其在海洋地震勘探数据上的反演试验.反演主要的计算步骤是:1.采用Fourier-Hankel变换把球面波分解为平面谐波;2.用平面谐波的延拓方程将上行波与下行波同时向下延拓,并计算每一层底部的反射系数和下一层的波阻抗;3.用最小二乘法从波阻抗中确定该层的声波速度.重复第2步与第3步,直到某一预定深度时结束.通过反演试验,对地震振幅比例的改变,子波变形,以及第1层速度和密度的误差对反演方法的稳定性及其精度的影响进行了分析.还通过实际海洋地震勘探数据的反演试验,对这一方法在地震勘探中的应用前景作了论述.  相似文献   

14.
寇岚  张进 《地震工程学报》2019,41(5):1259-1265
利用重力异常反演测试三维地震波速度结构,存在解不唯一、可靠性不高的问题。将面波反演充分融合到重力异常反演方程中,降低传统反演方法的非唯一性,并提升可靠性。以川滇地区为例,采用融合后的重力异常反演方法分析三维地震波速度结构。通过速度和密度的关系转换,得到对应的重力异常数据。由于面波频射数据主要对地震波横波速度敏感,因此将重力异常数据和初始横波速度相连,依据地震波速度和岩石密度之间的关系,获取重力异常反演方程,用于分析速度结构。选取21.6°~34.2°N、97.1°~105.9°E范围内的川滇地区活动块体作为实验数据,经过实验分析发现:使用该方法迭代反演川滇地区地壳上地幔顶部横波速度,重力异常数据和面波频射数据的残差值分别是6.24 mGal和0.027 km/s,实际拟合效果较好;分析该地区不同深度切面横波速度发现,在24 km深度处,上地壳中含有相对低速层,在44 km深度处,中下地壳中存在低速层;且该方法分析川滇地区三维地震波速度结构解的分辨率较高。  相似文献   

15.
海洋勘探环境可以抽象为下伏固体与上覆流体相互耦合的介质,本文针对流-固边界耦合介质提出了一种高效、稳定的多参数(速度和密度)全波形反演方法.本文采用弹性波一阶位移-应力方程作为过渡层耦合声波压力方程与弹性波位移方程来模拟耦合环境,相比于传统的交错网格建模方法或者构建连续性条件,本文提出的方法在正演精度和稳定性上凸显出很大优势,极大降低了计算内存.反演策略对多参数全波形反演至关重要,由于不同参数之间的相互耦合使得密度在多参数全波形反演中较难获得,因此本文将非均匀流-固边界耦合介质多参数全波形反演分为两个步骤完成:第一步利用变密度声波方程结合推导出的密度梯度算子进行纵波速度和密度的双参数反演;第二步根据链式法则求取横波速度的梯度,结合第一步的反演结果使用流-固边界耦合方程反演横波速度.最后通过与声波动方程数值模拟结果对比证明正演算法的准确性;上覆流体的Marmousi-2模型的数值试验测试说明反演方法的有效性和适应性.  相似文献   

16.
本文把波动方程反问题与广义Radon 变换的反演相联系。在假定弱散射条件下,把波动方程反问题转化成广义Radon 变换的反演问题,即如何从一系列关于目标函数在某类子流形上的积分值,去求出目标函数.这种转化提供了一种研究波动方程反问题的途径。  相似文献   

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

18.
对三维电阻率反演问题进行了深入研究,提供了一种利用地表观测数据实现三维反演的实用算法.该方法应用有限差分求正演解,并通过对粗糙度矩阵元素进行适当改进,使之适用于各种情况下粗糙度矩阵的求取,进而建立在模型的总粗糙度极小条件下的反演方程.对反演方程采用收敛速度快且稳定的最小二乘正交分解(LSQR)法进行迭代求解,在迭代求解过程中只需利用偏导数矩阵和其转置矩阵乘以一个向量的结果,回避了直接求偏导数矩阵的繁琐计算,节省了内存,加快了反演的计算速度.不同的计算实例表明上述方法是求解大规模三维电阻率反演问题的有效方法.  相似文献   

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

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