首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 564 毫秒
1.
反射地震数据的逐层波形反演   总被引:8,自引:2,他引:8       下载免费PDF全文
本文针对层状介质并结合梯度法波形反演,提出逐层波形反演的方法. 首先给出介质扰动响应的概念,并在此基础上分析了梯度法波形反演方法. 波形反演实质上是将实测地震记录和预测地震记录的波形残差信息转化为实际地质模型和预测地质模型的模型残差信息. 波形反演的优点是利用大量振幅相位信息得到高分辨率的反演结果, 其缺点是运行耗时大;当初始模型和实际模型相差较大时,迭代算法容易陷入局部极小点,这是因为目标函数和初始模型同实际模型间的差异是非线性的关系. 逐层波形反演方法是使自上而下每一层的目标函数最小,这样总的目标函数也是最小的. 利用二分法速度扫描确定每一层速度不仅提高了运算速度也避免了迭代算法陷入局部极小点的问题. 结合介质扰动响应和目标函数值变化可以更为准确迅速地确定每一层速度和该层界面位置.  相似文献   

2.
常规多波联合反演采用Zoeppritz方程的近似式构建正演方程,反演过程中需要假定背景纵横波速度比为常数,其反演精度不高,稳定性不好。本文提出了一种基于精确Zoeppritz方程的多波联合反演方法,结合贝叶斯方法进行广义线性反演。本方法基于精确Zoeppritz方程构建正演方程,避免了近似式反演在大角度时引起的误差;利用贝叶斯方法引入模型参数的先验分布信息,作为反演的正则化项,降低了反演的不适定性;反演目标函数中引入低频软约束,稳定了反演低频结果,提高了反演的鲁棒性;在求解反演目标函数时,利用快速算法,降低了反演的运算量。经过模型试算,证明了该方法的优越性和抗噪性;并在实际资料的应用中证明了该方法的实用性和有效性。  相似文献   

3.
本文为解决直流电阻率法和瑞雷面波法反演存在的多解性问题,构建了新的电震联合反演目标函数,通过增加不同地球物理响应的反演参数,利用权系数调节各参数对反演结果的影响,实现同步联合反演,从而减少多解性.引入人工蜂群算法求解目标函数,有效解决了传统反演方法将非线性问题线性化处理产生的细节信息丢失问题,提高了反演精度.经过对较为复杂的四层速度递增AA型和四层含低速夹层HA型理论模型和实例数据的反演验证,表明基于人工蜂群算法的非线性联合反演能够有效减少单一方法多解性问题并提高反演精度.  相似文献   

4.
多分量地震数据中低频缺失是弹性波全波形反演中的一大难题,低频的缺失导致全波形反演无法有效恢复介质的长波长成分进而使反演陷入局部极值。为此,本文提出了反演介质纵横波速度长波长分量的弹性波包络反演方法。该方法利用包络算子具有的解调多分量数据中隐含的低频信息的能力,构造多分量地震数据的包络目标函数进行反演,用以恢复地下介质纵横波速度的长波长成分。一系列数值试验表明,即使在多分量地震数据中缺失低频信息、并且初始模型缺少先验信息的情况下,这种弹性波包络反演方法能够有效降低波形反演的非线性,可以为后续的常规弹性波全波形反演或者深度偏移提供足够精确的初始模型,且该方法对横波速度长波长分量的重建尤为有效。Mamousi-2模型的高精度纵横波速度的反演结果表明,利用该方法反演的纵横波速度作为常规弹性波全波形反演的初始模型,可以显著提高反演结果的精度。此外,本文对弹性波包络反演方法的适用性也进行了初步的研究与讨论。  相似文献   

5.
包络信号含有丰富的低频分量,即使在地震数据缺失低频条件下,包络目标函数也能有效缓解全波形反演的周期跳跃现象.但是,当初始速度模型较为平滑时,观测数据中的反射地震事件在模拟数据中没有与之相对应的波形,导致包络反演初期无法很好地利用反射波信号进行速度建模.本文提出基于反射地震数据的时频域包络反演方法,通过结合反射波全波形反演理论,构建反射波时频域包络目标函数,来提高包络反演的速度建模精度.本文首先利用Gabor变换获取时频域地震数据,并提取振幅信息,即为时频域包络信号.然后,推导反射波时频域包络反演的伴随震源和梯度算子.Marmousi模型数据测试结果表明,基于反射地震数据的时频域包络反演方法可以为全波形反演提供一个较好的初始速度模型.  相似文献   

6.
包络信号含有丰富的低频分量,即使在地震数据缺失低频条件下,包络目标函数也能有效缓解全波形反演的周期跳跃现象.但是,当初始速度模型较为平滑时,观测数据中的反射地震事件在模拟数据中没有与之相对应的波形,导致包络反演初期无法很好地利用反射波信号进行速度建模.本文提出基于反射地震数据的时频域包络反演方法,通过结合反射波全波形反演理论,构建反射波时频域包络目标函数,来提高包络反演的速度建模精度.本文首先利用Gabor变换获取时频域地震数据,并提取振幅信息,即为时频域包络信号.然后,推导反射波时频域包络反演的伴随震源和梯度算子.Marmousi模型数据测试结果表明,基于反射地震数据的时频域包络反演方法可以为全波形反演提供一个较好的初始速度模型.  相似文献   

7.
速度、密度之间的相互耦合使得密度在多参数全波形反演中较难获得.本文将截断高斯-牛顿法用于声介质速度、密度双参数全波形反演,通过考虑近似Hessian矩阵中反映速度、密度相互作用的非主对角块元素,有效解决了多参数全波形反演中速度、密度之间的耦合问题,在不采用反演策略的情况下,仍能够获得精度较高的速度、密度反演结果.常规的截断牛顿类全波形反演通常利用一阶伴随状态法求取目标函数对模型参数的梯度,利用二阶伴随状态法或有限差分法求解Hessian-向量乘,在每一步内循环迭代过程中需要额外求解两次正演问题,计算量较大.本文基于Born近似,将梯度计算中的核函数-向量乘表示为具有明确物理意义的向量-标量乘的累加运算,同时将Hessian-向量乘转化为两次核函数-向量乘,无需额外求解正演问题,有效降低了计算量.数值实验证明了本文提出的方法的有效性.  相似文献   

8.
单斜介质中方位NMO速度Thomsen参数反演方法研究   总被引:9,自引:0,他引:9       下载免费PDF全文
采用启发式共轭梯度法,即随机爬山法 + 共轭梯度法,利用单斜介质中P波方位NMO速度椭圆轴向偏转角度接近于零这一特性,简化P波方位NMO速度公式,并利用多方位P波NMO速度,反演出某一初始CMP观测线与自然坐标系之间的夹角,作为进一步进行Thomsen各向异性参数反演的基础. 根据各向异性介质中方位NMO速度与Thomsen参数之间的关系,建立了利用三种波的多方位NMO速度及垂直传播速度反演单层单斜各向异性介质Thomsen各向异性参数的目标函数. 对计算的理论值添加具有一定标准差的正态分布的随机噪声,用以模拟实际观测存在的误差,通过对加噪后的数据进行多次反演的误差分析,表明了所建立的目标函数及选用的反演方法是有效可行的,而且相对稳定.  相似文献   

9.
求解非线性最优化问题的遗传算法   总被引:8,自引:1,他引:8  
一、引言地球物理反演是一类典型的非线性最优化问题,其目标是:根据观测得到的地球物理资料推断地下介质模型,即试图寻找一个能最好地解释(说明)观测结果的最佳拟合模型.一般来说,地球物理中的正演模拟计算是模型参数的非线性函数,因而地球物理反演是非线性最优化问题.而常用的目标函数或拟合函数(如归一化互相关函数)具有多个极值(峰或  相似文献   

10.
混合范数下的最优化反演方法   总被引:4,自引:1,他引:4  
在求解地球物理反问题时,通常根据最小二乘准则构造目标函数进行反演,并在实践中得到了广泛的应用.为进一步增强反演的稳健性及减少多解性,不损失反演结果的分辨率,本文提出了混合范数下的最优化反演方法,它根据数据和模型可能服从不同的概率分布,对数据空间和模型空间采用不同的范数来构造目标函数.在给出目标函数的基础上,导出了混合范数下的线性反演方程.由于该线性反演方程的复杂性,我们采用混合范数下迭代再加权共轭梯度法进行求解.最后,通过对模拟的电阻率数据进行反演,验证了本文计算方法是可行的.  相似文献   

11.
Although Genetic Algorithms have found many successful applications in the field of exploration geophysics, the convergence speed remains a big challenge as Genetic Algorithms usually require a huge amount of fitness function evaluations. In this paper, we propose an efficiency-improved Genetic Algorithm, which has both a good global search capability and a good local search capability, and is also capable of robustly handling the premature convergence challenge commonly seen in linear and directed non-linear optimization methods. In our new genetic algorithm, the global search capability is performed via a modified island model, while the local search capability is provided by a novel self-adaptive differential evolution fine tuning scheme. Premature convergence is dealt with via a local exhaustive search method. We first demonstrate the much improved convergence speed of this efficiency-improved Genetic Algorithm over that of our previously proposed advanced Genetic Algorithm on several multimodal functions. We further demonstrate the effectiveness of our efficiency-improved Genetic Algorithm by applying it to a two-dimensional common reflection surface stacking problem, which is a highly nonlinear geophysical optimization problem, to obtain very encouraging results.  相似文献   

12.
改进的地壳速度结构遗传反演方法   总被引:1,自引:0,他引:1  
遗传算法的收敛速度与模型参数初始搜索范围的划分精度有关.本文提出了在低精度的基础上开始迭代,在选代过程中逐步缩小搜索范围的改进措施.这种改进使得遗传算法的收敛速度和反演精度同时得到提高,并且又不增加模型空间的大小.一维地壳速度结构的反演计算结果表明方法是有效的.最后给出了实际应用中应采用的反演策略.  相似文献   

13.
重磁遗传算法三维反演中高速计算及有效存储方法技术   总被引:13,自引:15,他引:13       下载免费PDF全文
将地下场源区域规则划分成很多小长方体单元,并且通过反演确定这些单元的物性变 化,勾画出场源的分布图像,这种方式逐步成为重磁反演,特别是三维反演的重要方向;遗 传算法等非线性技术进行该类反演将逐步成为发展趋势. 本文指出,在应用遗传算法进行该 类反演过程中,隐含着数据量较大时超常规的计算量,它已成为制约该类反演充分发挥作用 的瓶颈问题;同时,本文提出了针对性的分离并存储几何格架的计算策略、以及独特的几何 格架等效压缩存储技术,可以从根本上提高非线性反演计算速度,为该类反演的有效应用奠 定了坚实的基础.  相似文献   

14.
正交各向异性介质P波走时分析及Thomsen参数反演   总被引:10,自引:3,他引:10  
对于包含有垂向裂缝的横向各向同性地层或含有多组正交裂缝的各向同性地层,正交各向异性介质模型是最简单的与实际地层相符的方位各向异性模型.本文对单层水平反射界面正交各向异性模型采用射线追踪法计算了全方位角变化的P波走时,时距曲线表现出强方位各向异性.采用小生境遗传算法,对三条成一定角度的测线的走时信息进行速度和各向异性参数反演.模型算例表明,此方法可以得到高精度的裂缝方位角、P波垂直速度和较高精度的Thomsen各向异性参数.  相似文献   

15.
利用Cole-Cole模型组合和稀释系数理论,可以描述测量得到的视复电阻率频谱与地下多个电性体SIP真参数之间的正演关系。本利用此正演关系进行反演试算,并比较了几种Cole-Cole模型反演方案的效果;提出固定极化目标体和围岩体的极化率同时反演其它各SIP真参数的反演方案。结果表明在固定极化率为真值时,利用此反演方案可以稳定、准确、快速地得出地下电性体SIP真参数。  相似文献   

16.
高分辨率非线性三维整体反演方法是基于非线性理论,在层位控制下,将工区多井(或全部井)的测井数据与井旁地震道数据输入具有多输入多输出的网络,同时进行整体训练,可获得整个工区的自适应权函数,并建立综合非线性映射关系,并根据储层在纵横方向上的地质变化特征更新这种非线性映射关系,这样,就能对反演过程及其反演结果起到约束和控制的作用,从而获得稳定且分辨率高的地震反演剖面(速度反演剖面/波阻抗反演剖面/密度反演剖面),实现整体反演,该方法通过模型试算和实际资料处理,获得较好的地质效果,证明该方法精度高、实用性强,可用于储层的定量分析。  相似文献   

17.
Surface wave dispersion curve inversion is a challenging problem for linear inversion procedures due to its highly non-linear nature and to the large numbers of local minima and maxima of the objective function (multi-modality). In order to improve the reliability of the inversion results, we implemented and tested a two-step inversion scheme based on Genetic Algorithms (GAs). The proposed scheme performs several preliminary “parallel” runs (first step) and a final global run using the previously-determined fittest models as starting population.In this work we focus on the inversion of shear-wave velocity and layer thickness while fixing compressional-wave velocity and density according to user-defined Poisson's ratios and velocity–density relationship respectively. The procedure can nonetheless perform the inversion under different degrees of regularization, depending on the a priori information and the desired degree of freedom of the system.Thanks to the large number of considered models, in addition to the fittest model, a mean model and its accuracy are evaluated by means of a statistical approach based on the estimation of the Marginal Posterior Probability Density (MPPD).We tested the proposed GA-based inversion scheme on three synthetic models reproducing a complex structure with low-to-moderate velocity cover (also including a low-velocity channel) lying over hard bedrock. For all the considered cases the bedrock velocity and depth were properly identified, and velocity inversion was reconstructed with minor uncertainties.The performed tests also investigate the influence of the first higher mode, the reduction of the frequency range of the considered dispersion curve as well as the use of different number of strata. While a limited frequency range of the dispersion curve (maximum frequency reduced from 80 to 40 Hz) does not seem to significantly limit the accuracy of the retrieved model, the adoption of the correct number of strata and the addition of the first higher mode help better focus the final solution.In conclusion, the proposed approach represents an improvement of a purely GA-based optimization scheme and the MPPD-based mean model typically offers a more significant and precise solution than the fittest one.Results of the inversion performed on a field data set were validated by borehole stratigraphy.  相似文献   

18.
使用阻尼最小二乘法进行震源参数和地壳三维速度结构的走时联合反演.所用资料为S波和P波到时差,并用人工地震资料的二维解释结果作为三维速度模型的特定约束条件.为建立初始模型,又利用天然地震构成了准二维剖面.在走时反演基础上,利用遗传算法进行了几个地震事件的波形反演尝试,并对走时反演获得的地壳速度结构模型的局部进行了修正.以34°~42°N,94°~112°E作为研究区域,在该区域中收集了1986年以来大量地震的S波和P波到时差资料,7条人工地震二维速度剖面资料和2个数字化地震台的几个地震的三分向记录资料.对这些资料进行了处理,最后得出了0~25km深度不同截面的速度分布,并对所得结果进行了分析.  相似文献   

19.
粘弹性参数变得越来越重要,其反演算法也逐渐成为众多研究者的研究热点。而遗传算法是一种随机、自适应、启发式的算法, 具有很好的鲁棒性和全局收敛性, 本文基于VSP直达波方程,引入了遗传算法来进行粘弹参数反演, 首先将频率域直达波方程表示为复速度的函数,然后通过遗传算法反演出复速度。而复速度和品质因子又是复速度的函数,从而便可很容易的得出。但若直接反演复速度, 反演参数太多, 不容易实现, 所以又将复速度表示成参数C0和C∞的函数,以减少反演参数数量。最后给出了理论模型实验,以证明该算法的有效性。  相似文献   

20.
A seismic inversion procedure is developed that inverts data available from an unmigrated stacked section to produce an interval velocity model. It attempts to overcome some of the limitations of existing methods by using a generalized linear inversion technique. The inversion process incorporates several features: (i) Lateral interval velocity variations are permitted, (ii) A fast accurate forward model was developed, (iii) Input data is weighted according to the accuracy with which it has been acquired. The procedure is applied to seismic data from the Gippsland Basin, an area offshore South-East Australia.  相似文献   

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

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