随机地震反演关键参数优选和效果分析(英文)

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

基于地质统计先验信息的储层物性参数同步反演

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

基于地质统计先验信息的储层物性参数同步反演(英文)

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

测井、井间地震与地面地震数据联合约束下的地质统计学随机建模方法研究

利用能够整合测井信息与井间地震信息的地质统计学随机模拟方法,结合传统的地质统计学反演思路,得到了一种能够同时整合测井、井间地震与地面地震三种先验信息的地质统计学反演与储层建模方法.由于井间射线信息、测井信息与地面地震数据在随机反演与建模过程当中都得到了尊重,因此与传统地质统计学反演仅利用了测井与地面地震数据相比,本文的地质统计学反演与建模方法更充分地利用了先验信息,有效提高了反演的精度,降低了随机建模中的多解性.基于理论数据的测试证实了上述观点.     

基于统计岩石物理模型的各向异性页岩储层参数反演

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

基于可变临界孔隙度模型的致密砂岩储层参数地震反演方法

储层参数地震反演是通过地震观测数据定量评估储层物性及含油气性的重要途径.中国西部致密砂岩油气分布范围广泛,具有巨大的勘探开发潜力,但其孔隙结构复杂,储层参数的定量预测难度较大.本文提出了一种基于可变临界孔隙度模型的储层参数地震反演方法.为有效描述储层的复杂孔隙结构,本方法将临界孔隙度作为变量,使其具有随深度变化的特征,以表征复杂孔隙结构,进而提高致密砂岩储层岩石物理建模的精度.此外,为直接从观测地震数据反演储层参数,结合岩石物理模型和地震反射系数方程建立正演算子;基于贝叶斯框架构建反演问题,采用混合高斯模型构建储层参数先验约束,用以描述储层的岩性(相)差异特征;基于临界孔隙度、泥质含量与孔隙度的统计关系,结合快速模拟退火算法,在迭代优化中使临界孔隙度与储层参数交替更新,用以提高反演结果稳定性.选用中国西部四川盆地某工区实际资料,对本方法进行应用测试与验证,结果表明本方法对于含气储层具有良好的指示性,相比常规固定临界孔隙度方法,本方法可以有效提升反演结果准确性.     

频变AVAZ响应特征分析及裂缝性质反演方法研究

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

伽马拟声波分频重构反演在储层预测中的应用

在地质和地震资料研究基础上,通过岩心、录井和测井分析,探索不同频率尺度测井曲线的储层意义,优选研究区内钻穿古近系底砂岩段砂泥薄互层储层5口探井的储层敏感曲线.在曲线预处理、标准化和敏感性分析的前提下,利用曲线分频拟合,对储层敏感曲线和声波曲线做小波多尺度分解,提取反映薄互层储层岩性变化的中频中尺度伽马曲线和反映地层速度的低频大尺度声波曲线.利用统计模型信息融合构建具有声波量纲的拟声波曲线.在地质研究和沉积演化模式的指导下,运用井约束反演的方法补充地震资料的高频和低频信息,提高地震资料的质量,预测有利储层分布.研究表明,伽马拟声波分频重构反演能有效的提高地震资料的纵向和横向分辨率,刻画复杂地质条件下的薄层储层分布.     

基于约束最小二乘与信赖域的储层参数反演方法

基于包体岩石物理模型的储层参数地震反演方法面临数学形式复杂、多解性强、适应性差、涉及迭代运算等问题,本文提出一种基于约束最小二乘与信赖域的储层参数地震反演方法.该方法基于储层参数与弹性参数关联岩石物理模型,使用最小二乘方法构建目标函数和信赖域约束全局寻优求解,有效降低了地震反演多解性,极大提高了收敛速度.特别是通过在最小二乘求解中引入垂向约束,有效提高反演结果的抗噪声能力.经过模型测试和实际资料的应用,验证了方法的可行性和适用性.     

地震岩相识别概率表征方法

储层岩相分布信息是油藏表征的重要参数,基于地震资料开展储层岩相识别通常具有较强的不确定性.传统方法仅获取唯一确定的岩相分布信息,无法解析反演结果的不确定性,增加了油藏评价的风险.本文引入基于概率统计的多步骤反演方法开展地震岩相识别,通过在其各个环节建立输入与输出参量的统计关系,然后融合各环节概率统计信息构建地震数据与储层岩相的条件概率关系以反演岩相分布概率信息.与传统方法相比,文中方法通过概率统计关系表征了地震岩相识别各个环节中地球物理响应关系的不确定性,并通过融合各环节概率信息实现了不确定性传递的数值模拟,最终反演的岩相概率信息能够客观准确地反映地震岩相识别结果的不确定性,为油藏评价及储层建模提供了重要参考信息.模型数据和实际资料应用验证了方法的有效性.     

高分辨率非线性储层物性参数反演方法和应用

对于陆相沉积环境下的复杂隐蔽岩性储层,由于观测信息不准确,如信息重叠、信息缺失和噪音污染,以及岩石物理关系模糊等原因,储层横向预测存在不惟一性、不稳定性和不确定性.基于线性假定的常规储层横向预测技术已不适用于复杂隐蔽岩性储层的勘探.本文采用一种非线性储层岩性物性褶积模型,建立波阻抗与孔隙度/泥质含量的函数关系;通过多级结构分解和双向边沿子波检测来刻画复杂岩石物理关系;通过Caianiello褶积神经网络实现确定性反演、统计反演和非线性理论三者有机结合;最后联合应用基于逆算子的反演方法和基于正算子的重建算法实现了综合地质、测井和地震波阻抗信息进行高分辨率储层物性参数反演.非线性储层物性参数反演采用多井约束机制和分频反演方式,在陆相和近海油气勘探资料的实际应用中,取得了明显应用效果.     

岩石物理驱动的正交各向异性方位叠前地震反演方法

基于长波长近似假设,周期性薄互层中发育一组平行排列的垂直裂缝则可视为等效的正交各向异性介质.岩石物理是构建裂缝参数与地震响应之间联系的基础,地震散射理论是各向异性介质参数反演的有效途径.文章提出了一种利用方位叠前地震数据实现正交各向异性裂缝储层Thomsen弱各向异性参数与裂缝弱度参数可靠预测的方法.首先,综合考虑矿物基质、孔隙、裂缝及各向异性岩石中流体替换的影响,通过构建正交各向异性裂缝岩石物理等效模型,实现正交各向异性刚度系数的估测,进而预测储层测井数据的弹性参数、Thomsen弱各向异性参数及裂缝弱度参数,为后续地震反演提供初始模型约束;然后,基于地震散射理论,推导了面向Thomsen弱各向异性参数与裂缝弱度参数反演的正交各向异性介质纵波反射系数方程,为后续地震反演奠定了理论基础;最后,发展了贝叶斯框架下的正交各向异性裂缝储层Thomsen弱各向异性参数与裂缝弱度参数AVAZ反演方法,同时考虑柯西稀疏约束正则化和平滑模型约束正则化约束,使用非线性的迭代重加权最小二乘策略实现正交各向异性特征参数的稳定估算.模型和实际资料处理表明,该方法能够稳定可靠地从方位叠前地震资料中获取正交各向异性特征参数,为正交各向异性介质的特征参数预测提供了一种高可靠性的地震反演方法.     

岩石物理驱动的相约束叠前地震概率化反演方法

地震岩石物理是连接岩石弹性参数与储层物性参数的桥梁,叠前地震反演是实现地下岩石弹性、物性、岩性及含流体性质定量表征的重要方法.文章构建了碎屑岩地震岩石物理高阶近似模型,推导了利用岩石模量高阶近似(Jacobian、Hessian矩阵)表征的叠前地震AVO反射特征方程,并分析了岩石孔隙度、泥质含量及流体饱和度对AVO反射率的贡献度,探讨了此方程在岩石物性参数直接预测方面的可行性.以此为基础,在待反演模型参数服从混合概率先验模型的前提下,文章提出了基于差分进化-马尔可夫链蒙特卡罗随机模型的相约束叠前地震概率化反演方法,兼具差分进化算法的全局寻优特性和马尔可夫链蒙特卡罗模型的不确定性分析能力;通过多条马尔可夫链的交叉并行,可以同步获得待反演模型参数的多个随机解,进而模拟待反演模型的后验概率密度分布,后验均值作为待反演模型的最优解,方差与置信区间用来评价反演结果的不确定性,实现储层弹性、物性、离散岩相及干岩石骨架等参数的同步预测.通过理论试验和实际资料处理验证了该理论方法的有效性.     

带先验约束的碎屑岩储层物性参数反演方法

对弹性参数进行转换而获取物性参数的方法,适用范围有限且预测精度不高.本文基于岩石物理理论,提出一种带先验约束的碎屑岩储层物性参数预测方法,为储层预测、流体识别、储层描述、储量估计、烃类开发方案设计提供重要参数.应用岩石物理理论和物性参数先验信息建立了带先验约束的目标函数;并利用最优化方法对反演目标函数进行优化求解,对储层的物性及孔隙度和泥质含量进行同步迭代反演.经过正演数据及实际地震资料的测试应用,其结果表明本文方法反演得到的孔隙度和泥质含量为储层综合预测提供了愈加丰富的评价指标,同时孔隙度和泥质含量数据的联合应用对减小储层预测的多解性和降低勘探开发的风险有重要意义.     

基于反假频基函数的井震匹配快速算法

河流相储层、碳酸盐岩孔缝洞密集带等地层,地震同相轴不具有层状特征,可形象地称为具有复杂网状特征,利用层状地层模型难以精细描述井外地层.传统的宽带约束反演,其目标函数一般含有一个层状的初始模型、离散光滑约束、稀疏约束和线性拟合函数约束.利用这种层状的初始模型进行井震匹配,其高频信息主要来源于井资料插值得到的层状地层模型,井之间地震信息很难起作用.不仅如此,其计算量较大,且不易调节层状信息与稀疏约束反映的地震资料中的非层状信息的比例,难以适应的复杂网状地层的井震信息融合.本文利用地震资料构造反假频基函数以描述复杂网状地层,避免了对层状模型的过分依赖,使得联合重建的信息在井旁与测井资料拟合较好,在井外能充分地反映地震资料所含非层状信息.本文利用大庆下白垩地层扶余油层的河流相资料进行实验,证实了方法的快速和有效性.     

大斜度井/水平井随钻方位电磁波测井资料实时反演方法

水平井钻进过程中,地层边界的准确预测对地质导向决策至关重要.但水平井测井环境复杂,随钻方位电磁波测井响应不直观,储层性质、地层界面等信息的精确、定量评价困难.本文提出了一种适用于大斜度井/水平井随钻方位电磁波测井资料的电阻率及地层界面信息的实时提取方法.该方法主要基于以下几个技术:(1)引入滑动开窗策略将高维随钻方位电磁波测井资料的反演问题转化为多个连续窗口的一维反演问题;(2)对反演参数进行多初值初始化,结合正则化Levenberg-Marquardt算法保证解的全局最优性;(3)充分利用随钻方位电磁波测井多种探测深度信息,采用分级反演策略提高反演速度.数值和实测资料反演结果表明:降维反演方法可实现对地层界面的实时、准确提取,且适用于多数复杂地层结构和任意井眼轨迹;利用浅探测地质信号可实现随钻方位电磁波测井资料的分级反演,最终反演速度可提升10倍以上.     

基于随机地震反演的Russell流体因子直接估算方法

本文研究了一种基于随机地震反演的Russell流体因子直接估算方法,该方法是一种基于蒙特卡罗的非线性反演,能够有效地融合测井资料中的高频信息,提高反演结果的分辨率.本文应用贝叶斯理论框架,首先通过测井数据计算井位置处的Russell流体因子,利用序贯高斯模拟方法(sequential Gaussian simulation,SGS)得到流体因子的先验信息;然后构建似然函数;最后利用Metropolis抽样算法对后验概率密度进行抽样,得到反演的Russell流体因子.其中对每道数据进行序贯高斯模拟时,采用一种新的逐点模拟方式,具有较高的计算速度.数值试验表明:反演结果与理论模型和实际测井数据吻合较好,具有较高的分辨率,对于判识储层含流体特征具有较好的指示作用.     

基于频变AVO技术对多尺度裂缝内流体属性反演与识别（英文）

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

叠前反演和地震吸收技术在复杂天然气藏地震预测中的应用

针对某复杂断块天然气目标储层,在岩石物理分析的指导下,综合利用地质、地震、测井等资料,提出了一套面向复杂天然气藏的叠前地震预测技术.首先基于地震岩石物理分析得到的初始横波信息,采用叠前贝叶斯非线性三参数反演得到了井旁控制点处精确纵横波速度和密度信息,然后通过叠前/叠后联合反演技术实现了面向目标的弹性阻抗体反演及含气储层敏感参数直接提取,最后结合小波变换时频谱分析的方法从叠前地震资料中估算地层吸收参数值,提高天然气藏识别精度.实际应用表明,综合各种叠前地震预测技术,可以大大提高对复杂天然气藏的识别精度,降低勘探风险.     

基于弹性参数加权统计的地震岩相预测方法

岩相信息能够反映储层岩性及流体特征,在地震储层预测中具有重要作用.常规方法主要利用与岩相信息关系密切的弹性参数定性或定量地转化为岩相信息.在实际应用中,弹性参数的获取主要基于叠前地震反演技术.而不同弹性参数的叠前地震反演精度间存在着差异,势必影响岩相的整体预测精度.本文提出对弹性参数进行加权统计来预测岩相.首先,基于贝叶斯理论,引入权重系数来调节弹性参数信息的采用量,构建出最终的目标反演函数;其次,考虑到勘探初期缺少明确的测井岩相信息,提出利用高斯混合分布函数来自动估算岩相先验概率;最后,根据输入弹性参数的取值,计算每类岩相对应的后验概率密度,将目标反演函数取最大后验概率密度时对应的岩相类别作为最终预测的岩相.新方法旨在减少弹性参数精度间的精度差异对岩相预测结果的影响,以期提高地震岩相的预测精度.模型与实际资料测试均表明该方法可行、有效且预测精度较高.