首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
常规AVO三参数反演是通过Zoeppritz方程的近似公式来建立AVO正演模拟的过程,然而在P波入射角过临界角和弹性参数在纵向上变化剧烈的情况下,Zoeppritz方程近似公式精度有限.针对这种情况,可以使用精确的Zoeppritz方程来构建反演目标函数,由于精确Zoeppritz方程中P波反射系数和弹性参数之间是一种复杂的非线性关系,通常解决途径是利用非线性的优化算法来进行数值计算,但是非线性优化算法的缺点是计算量过大;另外一种途径是利用广义线性反演的方法,通过泰勒一阶展开式将P波反射振幅展开后,用线性关系近似表达非线性关系,经过几次迭代后,在理论上可以达到很高的精度,但是广义线性反演算法的核心部分--Jacobian矩阵由于矩阵条件数过大,往往会造成反演算法的不稳定,其应用范围得到了限制.贝叶斯反演方法是通过引入模型参数的先验分布结合噪声的似然函数,生成模型参数的后验分布,通过求取模型参数的最大后验概率分布来得到模型参数的反演解,由于引入模型参数的先验分布信息,可以有效的降低反演的不适定问题.本文将两种反演算法的思想相结合,利用广义线性反演算法的思想,构建AVO正演模拟的过程来提高大角度地震数据反演的精度,同时结合贝叶斯理论,通过引入模型参数的先验分布信息构建反演目标函数的正则化项,可以有效降低由于Jacob矩阵条件数过大带来的反演不适定问题,该算法假设模型参数服从三变量柯西分布.  相似文献   

2.
AVAZ方法是裂缝储层评价的重要手段,然而它受到静态等效地震岩石物理模型的限制,一般用于反演裂缝密度和裂缝走向.近年来,随着裂缝性孔隙介质等效模型的发展,频变AVAZ响应特征被证明可以携带更多的裂缝性质信息.本文针对HTI介质,基于Chapman理论分析流体类型和裂缝长度对频变AVAZ响应特征的影响并建立结合贝叶斯理论和遗传算法反演裂缝密度、裂缝走向、裂缝长度和裂缝充填流体类型的新方法.新反演方法流程为:首先,基于频变AVAZ特征建立反演目标函数并运用遗传算法求解,根据此结果分别建立裂缝参数的似然函数;其次,利用贝叶斯算法将地质、测井等获得的裂缝性质先验信息与地震获得的裂缝性质似然函数结合,计算裂缝性质的后验概率分布函数;最后,通过控制先验裂缝参数的方差,均衡地震数据及地质、测井数据的信息采用量.合成数据的分析和反演结果表明:流体性质和裂缝长度是频变AVAZ响应的敏感参数;本文提出反演裂缝性质方法是可行的.先验裂缝参数影响分析验证了该方法在均衡地质、测井及地震数据的信息采用量方面的有效性.  相似文献   

3.
贝叶斯概率统计方法能够在实际样本空间容量较小的条件下,通过未知参数先验信息进一步对未知参数进行估计,与经典统计频率学存在本质的不同。收集了83组剪跨比在1~2的钢筋混凝土柱在低周反复荷载作用下的受剪承载力试验数据,选定美国ACI318-08规范中钢筋混凝土柱受剪承载力计算公式作为贝叶斯先验模型,综合这两类信息,采用贝叶斯概率统计理论,建立了基于贝叶斯理论的钢筋混凝土柱受剪承载力后验概率模型。利用贝叶斯后验模型参数剔除理论,剔除影响柱抗剪强度的次要因素,对概率模型进行动态更新,得到了柱抗剪强度的简化概率模型。最后,分析验证了贝叶斯概率简化柱受剪承载力计算模型的正确性,同时,对采用不同先验模型条件下所得的后验模型计算值进行了对比分析,表明贝叶斯方法继承了先验信息的完备性和大量试验数据的准确性,能够更准确地预测钢筋混凝土柱的受剪承载力。  相似文献   

4.
含噪声数据反演的概率描述   总被引:5,自引:4,他引:1       下载免费PDF全文
根据贝叶斯理论给出了对含噪声地球物理数据处理的具体流程和方法,主要包括似然函数估计和后验概率计算.我们将数据向量的概念扩展为数据向量的集合,通过引入数据空间内的信赖度,把数据噪声转移到模型空间的概率密度函数上,即获得了反映数据本身的不确定性的似然函数.该方法由于避免了处理阶段数据空间内的人工干预,因而可以保证模型空间中的概率密度单纯反映数据噪声,具有信息保真度高、保留可行解的优点.为了得到加入先验信息的后验分布,本文提出了使用加权矩阵的概率分析法,该方法在模型空间直接引入地质信息,对噪声引起的反演多解性有很强的约束效果.整个处理流程均以大地电磁反演为例进行了展示.  相似文献   

5.
基于弹性阻抗的储层物性参数预测方法   总被引:12,自引:9,他引:3       下载免费PDF全文
储层物性参数是储层描述的重要参数,常规的基于贝叶斯理论的储层物性参数反演方法大多是通过反演获得的弹性参数进一步转换而获得物性参数,本文提出一种基于弹性阻抗数据预测储层物性参数的反演方法.该方法主要通过建立可以表征弹性阻抗与储层物性参数之间关系的统计岩石物理模型,联合蒙特卡罗仿真模拟技术,在贝叶斯理论框架的指导下,应用期望最大化算法估计物性参数的后验概率分布,最终实现储层物性参数反演.经过模型测试和实际资料的处理,其结果表明本文提出的方法具有预测精度高,稳定性强,横向连续性好等优点.  相似文献   

6.
瑞利波具有强能量和频散特性,面波勘探采集瑞利面波数据分析反演得到横波速度结构,在浅地表勘探领域得到了广泛的应用,根据面波频散反演速度结构是面波勘探的重要环节之一.面波频散与地下介质的弹性参数是非线性关系,全局优化方法是解决非线性反问题的有效办法,贝叶斯方法是一种基于统计的全局优化方法.贝叶斯公式反映了先验信息和条件概率乘积与后验概率之间的正比关系,利用马尔科夫链蒙特卡洛采样方法可以获得后验分布的采样.与其他方法相比,该方法不是给出一个模型最优解,而是统计出模型参数的平均值和方差,平均值即是反演的模型解,方差能对该反演结果的不确定性做出评价.本文采用贝叶斯方法对模拟和实测瑞利面波数据频散曲线进行反演,结果显示,该方法能获得比较精确的横波速度和厚度参数,验证了本文方法可行.在高丽营地区采集了多个单炮面波数据,对数据进行频散谱处理拾取频散曲线,通过贝叶斯反演方法获得了该测线二维横波速度剖面图,有助于划分和解释该测区地层及断层地质构造.  相似文献   

7.
随机地震反演关键参数优选和效果分析(英文)   总被引:2,自引:0,他引:2  
随机地震反演技术是将地质统计理论和地震反演相结合的反演方法,它将地震资料、测井资料和地质统计学信息融合为地下模型的后验概率分布,利用马尔科夫链蒙特卡洛(MCMC)方法对该后验概率分布采样,通过综合分析多个采样结果来研究后验概率分布的性质,进而认识地下情况。本文首先介绍了随机地震反演的原理,然后对影响随机地震反演效果的四个关键参数,即地震资料信噪比、变差函数、后验概率分布的样本个数和井网密度进行分析并给出其优化原则。资料分析表明地震资料信噪比控制地震资料和地质统计规律对反演结果的约束程度,变差函数影响反演结果的平滑程度,后验概率分布的样本个数决定样本统计特征的可靠性,而参与反演的井网密度则影响反演的不确定性。最后通过对比试验工区随机地震反演和基于模型的确定性地震反演结果,指出随机地震反演可以给出更符合地下实际情况的模型。  相似文献   

8.
基于MCMC的叠前地震反演方法研究   总被引:6,自引:5,他引:1       下载免费PDF全文
马尔科夫链蒙特卡洛方法(MCMC)是一种启发式的全局寻优算法[1].它在贝叶斯框架下,利用已有资料进行约束,既可使最优解满足参数的统计特性,又通过融入的先验信息,提高解的精度;寻优过程可跳出局部最优,得到全局最优解.利用MCMC方法,可以得到大量来自于后验概率分布的样本,不仅可以得到每个未知参数的估计值,而且可以得到与...  相似文献   

9.
针对常规大地电磁(Magnetotelluric,MT)反演方法对电阻率异常体边界不太敏感的问题,本文尝试基于贝叶斯理论开展二维大地电磁电阻率尖锐边界反演研究.在反演中,模型参数由边界位置及内部电阻率组成,通过贝叶斯理论将模型参数与数据相联系,采用Markov Chain Monte Carlo(MCMC)的Metropolis-Hastings(MH)方法对后验概率密度函数(Posteriori Probability Density,PDD)进行采样.采样过程中无罚值函数约束,完全以数据自身所包含的信息对模型进行约束,同时与有限约束进行比较,并考虑不同起始采样点对结果的影响.以接受率为参考,用模型算例说明MH方法中建议分布函数选择的重要性.当模型参数间相关性较弱时,使用边缘概率分布对采样结果进行分析.该方法能给出模型参数的分布范围,并给出该模型参数范围对应的数据范围.通过与已知模型的对比及数据拟合情况分析检验了该反演方法的有效性.该方法有助于提高大地电磁尖锐边界反演的分辨能力.  相似文献   

10.
通过地震数据获取裂缝储藏中流体的性质并对流体类型进行识别,是地震勘探岩性反演的重要问题之一。由于地震波的速度、储层的密度等弹性参数对某些流体不具有很强的敏感性,使只依赖振幅信息进行流体识别的传统AVO方法面临困境。作为传统叠前振幅反演的一个拓展,频变AVO(FDAVO)技术进一步考虑了振幅对频率的依赖关系,将这种依赖关系与地下裂缝结构、流体填充对应起来,能带来更丰富的流体信息。利用该技术,本文提出了一种基于地震数据参数化Chapman模型的贝叶斯反演新方法(BIDCMP),它包含两步算法,即,FDAVO反演储层的非弹性属性和贝叶斯框架下的流体识别。首先,通过匹配观测数据和模型数据,构造差函数反演裂缝储层非弹性参数。随后,在贝叶斯框架下,使用马尔科夫随机场(MRF)作为先验模型,联合多参数场识别流体。本方法在计算过程中,除综合考虑了弹性参数场、测井资料等常规信息外,还特别地加人了第一步中反演得的非弹性参数的约束,从而充分利用了流体粘性差异,最后在最大后验概率(MAP)准则下输出最佳岩性一流体识别结果。分别对合成地震记录和模拟岩性—流体剖面验证本文方法的有效性,结果证明本文方法获得的流体识别结果准确可信。  相似文献   

11.
秦军 《中国科学D辑》2005,35(8):790-798
遥感反演中使用先验知识, 特别是在观测信息不足情况下, 对于保证反演计算的稳定性和降低反演结果的不确定性具有重要意义. 常用的最优算法, 不能较好提供先验知识的后验分布, 也就不能直接得到反演结果不确定性, 这对于遥感应用是不利的. 引入集合卡曼滤波(EnKF)的方法来实现地表参数的遥感反演, 它能在有效获得反演结果的同时给出先验知识的后验分布. 为了显示其反演优势, 将MODIS标准AMBRALS算法、SCE-UA高效全局最优算法和集合卡曼滤波算法进行了比较. 着重讨论了在观测不足时, 引入先验知识并以反演结果的不确定性为比较指标的情况下, 不同核组合的核驱动模型在不同地表类型上对反照率的反演能力.  相似文献   

12.
遥感反演中使用先验知识, 特别是在观测信息不足情况下, 对于保证反演计算的稳定性和降低反演结果的不确定性具有重要意义. 常用的最优算法, 不能较好提供先验知识的后验分布, 也就不能直接得到反演结果不确定性, 这对于遥感应用是不利的. 引入集合卡曼滤波(EnKF)的方法来实现地表参数的遥感反演, 它能在有效获得反演结果的同时给出先验知识的后验分布. 为了显示其反演优势, 将MODIS标准AMBRALS算法、SCE-UA高效全局最优算法和集合卡曼滤波算法进行了比较. 着重讨论了在观测不足时, 引入先验知识并以反演结果的不确定性为比较指标的情况下, 不同核组合的核驱动模型在不同地表类型上对反照率的反演能力.  相似文献   

13.
基于地质统计先验信息的储层物性参数同步反演   总被引:4,自引:1,他引:3  
本文提出的储层物性参数同步反演是一种高分辨率的非线性反演方法,该方法综合利用岩石物理和地质统计先验信息,在贝叶斯理论框架下,首先通过变差结构分析得到合理的变差函数,进而利用快速傅里叶滑动平均模拟算法(Fast Fourier TransformMoving Average,FFT-MA)和逐渐变形算法(Gradual Deformation Method,GDM)得到基于地质统计学的储层物性参数先验信息,然后根据统计岩石物理模型建立弹性参数与储层物性参数之间的关系,构建似然函数,最终利用Metropolis算法实现后验概率密度的抽样,得到物性参数反演结果。并将此方法处理了中国陆上探区的一块实际资料,本方法的反演结果具有较高的分辨率,与测井数据吻合度较高;由于可以直接反演储层物性参数,避免了误差的累积,大大减少了不确定性的传递,且计算效率较高。  相似文献   

14.
本文提出的储层物性参数同步反演是一种高分辨率的非线性反演方法,该方法综合利用岩石物理和地质统计先验信息,在贝叶斯理论框架下,首先通过变差结构分析得到合理的变差函数,进而利用快速傅里叶滑动平均模拟算法(Fast Fourier TransformMoving Average,FFT-MA)和逐渐变形算法(Gradual Deformation Method,GDM)得到基于地质统计学的储层物性参数先验信息,然后根据统计岩石物理模型建立弹性参数与储层物性参数之间的关系,构建似然函数,最终利用Metropolis算法实现后验概率密度的抽样,得到物性参数反演结果。并将此方法处理了中国陆上探区的一块实际资料,本方法的反演结果具有较高的分辨率,与测井数据吻合度较高;由于可以直接反演储层物性参数,避免了误差的累积,大大减少了不确定性的传递,且计算效率较高。  相似文献   

15.
目前,国内外针对钢筋混凝土梁受剪能力计算主要分为模型化方法和统计计算理论,但仍未形成统一定论。基于贝叶斯统计理论,综合先验模型和先验信息,建立了钢筋混凝土梁贝叶斯后验概率模型,并采用其独特的参数剔除方法,得到简化模型,进而完成了83组试验数据的试验验证。研究表明:在我国规范(GB50010-2010)和美国规范(ACI318-11)基础上建立的梁贝叶斯概率模型,能较好的应用于钢筋混凝土梁受剪承载力的计算,且计算结果与试验值吻合较好。  相似文献   

16.
提出了各向异性页岩储层统计岩石物理反演方法.通过统计岩石物理模型建立储层物性参数与弹性参数的定量关系,使用测井数据及井中岩石物理反演结果作为先验信息,将地震阻抗数据定量解释为储层物性参数、各向异性参数的空间分布.反演过程在贝叶斯框架下求得储层参数的后验概率密度函数,并从中得到参数的最优估计值及其不确定性的定量描述.在此过程中综合考虑了岩石物理模型对复杂地下介质的描述偏差和地震数据中噪声对反演不确定性的影响.在求取最大后验概率过程中使用模拟退火优化粒子群算法以提高收敛速度和计算准确性.将统计岩石物理技术应用于龙马溪组页岩气储层,得到储层泥质含量、压实指数、孔隙度、裂缝密度等物性,以及各向异性参数的空间分布及相应的不确定性估计,为页岩气储层的定量描述提供依据.  相似文献   

17.
针对传统神经网络在电阻率成像反演中存在的过拟合和易陷入局部极值等问题,提出了一种基于剪枝贝叶斯神经网络(PBNN)的非线性反演算法和一种基于K-medoids聚类的样本设计方法。在基于K-medoids聚类的样本设计方法中,利用观测数据的聚类结果提供先验信息构造神经网络的训练样本,从而有针对性地指导神经网络的训练过程;剪枝贝叶斯神经网络是在贝叶斯正则化的基础上,通过评估各隐节点对反演结果的影响来自适应确定神经网络的隐层结构,根据小样本条件下训练样本的先验分布特征,选择了基于广义平均的超参数αk来引导剪枝过程。通过与地球物理领域内其它常用的自适应正则化方法相比较,验证了本文算法的有效性。理论数据和实测数据反演的结果表明:该方法能够较好地抑制神经网络训练过程中噪声的影响,提高网络的泛化能力,其反演结果优于BPNN反演、RBFNN反演和RRBFNN反演以及传统的最小二乘反演。  相似文献   

18.
马尔科夫链蒙特卡洛方法(MCMC)是一种启发式的全局寻优算法,可以用来解决概率反演的问题.基于MCMC方法的反演不依赖于准确的初始模型,可以引入任意复杂的先验信息,通过对先验概率密度函数的采样来获得大量的后验概率分布样本,在寻找最优解的过程中可以跳出局部最优得到全局最优解.MCMC方法由于计算量巨大,应用难度较高,在地...  相似文献   

19.
针对Mogi模型垂直位移与水平位移联合反演中的病态问题,改进火山形变总体最小二乘(Total Least Squares,TLS)联合反演的虚拟观测法,并使用方差分量估计(Variance Components Estimation,VCE)方法确定病态问题的正则化参数.将附有先验信息的参数作为观测方程,与垂直位移和水平位移的观测方程联合解算,推导了三类观测方程联合反演的求解公式及基于总体最小二乘方差分量估计确定正则化参数的表达式,给出了算法的迭代流程.通过算例实验,研究了总体最小二乘联合反演的虚拟观测法在火山Mogi模型形变反演中的应用;算例结果表明,三类数据的联合平差及方差分量估计方法可以确定权比因子并得到修正后的压力源参数,具有一定的实际参考价值.  相似文献   

20.
三维重力反演是地质工作者了解地球深部构造,认知地下结构的重要手段.按照反演单元划分,三维重力反演有离散多面体(Discrete)反演和网格节点(Voxels)反演两种方式.离散多面体反演由于易于吸收先验地质信息得到的理论场能够很好地拟合观测场,因此,在实际重力反演中更受欢迎.目前离散多面体重力反演中初始模型的建立方法繁杂不一,实际应用受到很大的限制.本文本着充分挖掘利用先验信息和重力观测数据得到丰富可靠的反演结果这一原则,以离散多面体反演技术为基础,改进建模过程.在初始模型的建立中,吸收贝叶斯算法优势,采用隐马尔科夫链改善朴素贝叶斯方法的分类效果,通过最大似然函数算法求解,再采取模型降阶技术,固定所建模型中几何体的形态或密度,达到在几何体形态(x,y,z)、密度(σ)和重力值(g)五个参数中降低维数目的,从而减小高维不确定性和正演的计算量,由此反演计算的地质体密度和分布范围相对更准确,更利于重现重力模型结构.通过单位球体和任意形态几何体模拟实验,以及安徽省泥河矿区三维重力反演实践,得到非常接近实际的密度或重力值,大幅提高了三维重力反演的精度和效率,说明该方法是有效、实用的.  相似文献   

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

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