好文档就是一把金锄头!
欢迎来到金锄头文库![会员中心]
电子文档交易市场
安卓APP | ios版本
电子文档交易市场
安卓APP | ios版本

第五章预测技术及软件-精品文档资料.ppt

238页
  • 卖家[上传人]:枫**
  • 文档编号:591954431
  • 上传时间:2024-09-19
  • 文档格式:PPT
  • 文档大小:32.07MB
  • / 238 举报 版权申诉 马上下载
  • 文本预览
  • 下载提示
  • 常见问题
    • 油气综合预测技术及应用刘豪2008年2月第五章 主要油气预测技术及软件介绍主要油气预测技术及软件介绍 第五章5.1 叠后波阻抗反演叠后波阻抗反演5.2 叠后地震属性叠后地震属性5.3 利用地震参数进行储层预测的主要方利用地震参数进行储层预测的主要方法法 5.1 叠后波阻抗反演5.1.1 地震反演主要类型5.1.2 叠后地震反演存在的问题5.1.3 反演的技术关键 5.1.4 国内外主要反演软件 地震反演技术在20世纪70年代开始出现,以基于褶积模型的叠后一维波阻抗反演为主;80年 代 得 到 了 蓬 勃 发 展 1983年 ,Cooke介绍了地震资料广义线性反演方法,从而揭开了波阻抗反演技术的新篇章90年代至今围绕一维波阻抗反演的各类算法以及应用成果层出不穷 5.1.1 地震反演的主要类型地震反演的主要类型 道积分反演方法道积分反演方法 递推反演方法递推反演方法 基于模型的反演方法基于模型的反演方法 叠后波阻抗反演叠后波阻抗反演 道积分是利用迭后地震资料计算地层相道积分是利用迭后地震资料计算地层相对波阻抗(速度)的直接反演方法因对波阻抗(速度)的直接反演方法因为它是在地层波阻抗随深度连续可微条为它是在地层波阻抗随深度连续可微条件下推导出来的,因而又称连续反演。

      件下推导出来的,因而又称连续反演 道积分反演方法道积分反演方法地震反演的主要类型地震反演的主要类型Rdt 道积分反演方法道积分反演方法地震反演的主要类型地震反演的主要类型优点:优点:(1)递推时累计误差小;(2)计算简单,不需要反射系数的标定;(3)没有井的控制也能作缺点:缺点:(1)由于这种方法受地震固有频宽的限制,分辨率低,无法适应薄层解释的需要;(2)无法求得地层的绝对波阻抗和绝对速度,不能用于定量计算储集层参数;(3)这种方法在处理过程中不能用地质或测井资料对其进行约束控制,因而其结果比较粗略 递推反演是对地震资料的处理过程,其结果的分辨率、信噪比以及可靠程度主要依赖于地震资料本身的品质,因此用于反演的地震资料应具有较宽的频带、较低的噪声、相对振幅保持和准确成像测井资料,尤其是声波测井和密度测井资料,是地震横向预测的对比标准和解释依据,在反演处理之前应仔细校正,使其能够正确反映岩层的物理特征 递推反演方法递推反演方法地震反演的主要类型地震反演的主要类型 递推反演方法递推反演方法地震反演的主要类型地震反演的主要类型 递推反演的技术核心在于由地震递推反演的技术核心在于由地震资料正确估算地层反射系数(或消除资料正确估算地层反射系数(或消除地震子波的影响),比较典型的实现地震子波的影响),比较典型的实现方法有:基于地层反褶积方法、稀疏方法有:基于地层反褶积方法、稀疏脉冲反演和测井控制地震反演等。

      脉冲反演和测井控制地震反演等 比较典型的实现方法有:基于地比较典型的实现方法有:基于地层反褶积方法、稀疏脉冲反演和测井层反褶积方法、稀疏脉冲反演和测井控制地震反演等控制地震反演等 递推反演方法具有较宽的应用领域在勘探初期递推反演方法具有较宽的应用领域在勘探初期递推反演方法具有较宽的应用领域在勘探初期递推反演方法具有较宽的应用领域在勘探初期只有很少钻井的条件下,通过反演资料进行岩相只有很少钻井的条件下,通过反演资料进行岩相只有很少钻井的条件下,通过反演资料进行岩相只有很少钻井的条件下,通过反演资料进行岩相分析确定地层的沉积体系,根据钻井揭示的储集分析确定地层的沉积体系,根据钻井揭示的储集分析确定地层的沉积体系,根据钻井揭示的储集分析确定地层的沉积体系,根据钻井揭示的储集层特征进行横向预测,确定评价井位到开发前层特征进行横向预测,确定评价井位到开发前层特征进行横向预测,确定评价井位到开发前层特征进行横向预测,确定评价井位到开发前期,在储集层较厚的条件下,递推反演资料可为期,在储集层较厚的条件下,递推反演资料可为期,在储集层较厚的条件下,递推反演资料可为期,在储集层较厚的条件下,递推反演资料可为地质建模提供较可靠的构造、厚度和物性信息,地质建模提供较可靠的构造、厚度和物性信息,地质建模提供较可靠的构造、厚度和物性信息,地质建模提供较可靠的构造、厚度和物性信息,优化方案设计。

      在油藏监测阶段,通过时延地震优化方案设计在油藏监测阶段,通过时延地震优化方案设计在油藏监测阶段,通过时延地震优化方案设计在油藏监测阶段,通过时延地震反演速度差异分析,可帮助确定储集层压力、物反演速度差异分析,可帮助确定储集层压力、物反演速度差异分析,可帮助确定储集层压力、物反演速度差异分析,可帮助确定储集层压力、物性的空间变化,进而推断油气前缘性的空间变化,进而推断油气前缘性的空间变化,进而推断油气前缘性的空间变化,进而推断油气前缘由于受地震频带宽度的限制,递推反演资料的分由于受地震频带宽度的限制,递推反演资料的分由于受地震频带宽度的限制,递推反演资料的分由于受地震频带宽度的限制,递推反演资料的分辨率相对较低,不能满足薄储集层研究的需要辨率相对较低,不能满足薄储集层研究的需要辨率相对较低,不能满足薄储集层研究的需要辨率相对较低,不能满足薄储集层研究的需要 递推反演方法递推反演方法地震反演的主要类型地震反演的主要类型 一些难点 反射剖面的极性问题;反射剖面的极性问题;标定问题;标定问题;低频分量的补偿问题低频分量的补偿问题递推反演方法递推反演方法地震反演的主要类型地震反演的主要类型 基于模型的反演方法的基本思路是:先建立一基于模型的反演方法的基本思路是:先建立一基于模型的反演方法的基本思路是:先建立一基于模型的反演方法的基本思路是:先建立一个初始地层波阻抗模型,然后由此模型进行地个初始地层波阻抗模型,然后由此模型进行地个初始地层波阻抗模型,然后由此模型进行地个初始地层波阻抗模型,然后由此模型进行地震正演,求得合成地震记录,将合成地震记录震正演,求得合成地震记录,将合成地震记录震正演,求得合成地震记录,将合成地震记录震正演,求得合成地震记录,将合成地震记录与实际地震记录相比较,根据比较结果,修改与实际地震记录相比较,根据比较结果,修改与实际地震记录相比较,根据比较结果,修改与实际地震记录相比较,根据比较结果,修改地下波阻抗模型的速度、密度、深度值及子波,地下波阻抗模型的速度、密度、深度值及子波,地下波阻抗模型的速度、密度、深度值及子波,地下波阻抗模型的速度、密度、深度值及子波,再正演求取合成地震记录,与实际地震记录比再正演求取合成地震记录,与实际地震记录比再正演求取合成地震记录,与实际地震记录比再正演求取合成地震记录,与实际地震记录比较后,继续修改波阻抗模型,如此多次反复,较后,继续修改波阻抗模型,如此多次反复,较后,继续修改波阻抗模型,如此多次反复,较后,继续修改波阻抗模型,如此多次反复,从而不断地通过迭代修改,直至合成地震记录从而不断地通过迭代修改,直至合成地震记录从而不断地通过迭代修改,直至合成地震记录从而不断地通过迭代修改,直至合成地震记录与实际地震记录最接近,最终得到地下的波阻与实际地震记录最接近,最终得到地下的波阻与实际地震记录最接近,最终得到地下的波阻与实际地震记录最接近,最终得到地下的波阻抗模型。

      抗模型地震反演的主要类型地震反演的主要类型基于模型的反演方法基于模型的反演方法 优点优点避免了一般反褶积方法对子波的最小相位假设;避免了一般反褶积方法对子波的最小相位假设;避免了一般反褶积方法对子波的最小相位假设;避免了一般反褶积方法对子波的最小相位假设;不需要假设反射系数是白噪;不需要假设反射系数是白噪;不需要假设反射系数是白噪;不需要假设反射系数是白噪;这种方法还可以使随机干扰不参与反演因为该方法要求工作人这种方法还可以使随机干扰不参与反演因为该方法要求工作人这种方法还可以使随机干扰不参与反演因为该方法要求工作人这种方法还可以使随机干扰不参与反演因为该方法要求工作人员在一条剖面上只选择少数员在一条剖面上只选择少数员在一条剖面上只选择少数员在一条剖面上只选择少数““““控制道控制道控制道控制道””””,迭代过程中,修改厚度、,迭代过程中,修改厚度、,迭代过程中,修改厚度、,迭代过程中,修改厚度、子波、密度和速度只是在这些控制道上进行有了控制道参数以子波、密度和速度只是在这些控制道上进行有了控制道参数以子波、密度和速度只是在这些控制道上进行有了控制道参数以子波、密度和速度只是在这些控制道上进行。

      有了控制道参数以后,模型就在这些后,模型就在这些后,模型就在这些后,模型就在这些““““控制道控制道控制道控制道””””之间作内插,再用内插结果去作正之间作内插,再用内插结果去作正之间作内插,再用内插结果去作正之间作内插,再用内插结果去作正演这样一方面加快了运算速度;另一方面,由于程序不要求每演这样一方面加快了运算速度;另一方面,由于程序不要求每演这样一方面加快了运算速度;另一方面,由于程序不要求每演这样一方面加快了运算速度;另一方面,由于程序不要求每一道都与正演结果完全吻合,而只要求整个段上的正演结果与实一道都与正演结果完全吻合,而只要求整个段上的正演结果与实一道都与正演结果完全吻合,而只要求整个段上的正演结果与实一道都与正演结果完全吻合,而只要求整个段上的正演结果与实际剖面误差最小,最终的正演结果减去迭偏剖面所得到的只是一际剖面误差最小,最终的正演结果减去迭偏剖面所得到的只是一际剖面误差最小,最终的正演结果减去迭偏剖面所得到的只是一际剖面误差最小,最终的正演结果减去迭偏剖面所得到的只是一些随机干扰些随机干扰些随机干扰些随机干扰以测井资料丰富的高频信息和完整的低频成分补充地震有限带宽以测井资料丰富的高频信息和完整的低频成分补充地震有限带宽以测井资料丰富的高频信息和完整的低频成分补充地震有限带宽以测井资料丰富的高频信息和完整的低频成分补充地震有限带宽的不足,可获得高分辨率的地层波阻抗资料,其结果的低、高频的不足,可获得高分辨率的地层波阻抗资料,其结果的低、高频的不足,可获得高分辨率的地层波阻抗资料,其结果的低、高频的不足,可获得高分辨率的地层波阻抗资料,其结果的低、高频信息来源于测井资料,构造特征及中频段取决于地震数据,从而信息来源于测井资料,构造特征及中频段取决于地震数据,从而信息来源于测井资料,构造特征及中频段取决于地震数据,从而信息来源于测井资料,构造特征及中频段取决于地震数据,从而解决递推反演中难以解决的低频速度补偿问题解决递推反演中难以解决的低频速度补偿问题解决递推反演中难以解决的低频速度补偿问题解决递推反演中难以解决的低频速度补偿问题。

      地震反演的主要类型地震反演的主要类型基于模型的反演方法基于模型的反演方法 主要技术环节主要技术环节储储储储集集集集层层层层地地地地球球球球物物物物理理理理特特特特征征征征分分分分析析析析测测测测井井井井资资资资料料料料,,,,尤尤尤尤其其其其是是是是声声声声波波波波和和和和密密密密度度度度测测测测井井井井,,,,为为为为建建建建立立立立初初初初始始始始模模模模型型型型的的的的基基基基础础础础资资资资料料料料和和和和地地地地质质质质解解解解释释释释的的的的基基基基本本本本依依依依据据据据用用用用于于于于制制制制做做做做初初初初始始始始波波波波阻阻阻阻抗抗抗抗模模模模型的测井资料必须经过型的测井资料必须经过型的测井资料必须经过型的测井资料必须经过环境校正环境校正环境校正环境校正地地地地震震震震子子子子波波波波提提提提取取取取子子子子波波波波是是是是基基基基于于于于模模模模型型型型反反反反演演演演中中中中的的的的关关关关键键键键因因因因素建建建建立立立立初初初初始始始始波波波波阻阻阻阻抗抗抗抗模模模模型型型型建建建建立立立立尽尽尽尽可可可可能能能能接接接接近近近近实实实实际际际际地地地地层层层层情情情情况况况况的的的的波波波波阻阻阻阻抗抗抗抗模模模模型型型型,,,,是是是是减减减减少少少少其其其其最最最最终终终终结结结结果果果果多多多多解解解解性性性性的根本途径。

      的根本途径的根本途径的根本途径地震反演的主要类型地震反演的主要类型基于模型的反演方法基于模型的反演方法 地震反演的主要类型地震反演的主要类型基于模型的反演方法基于模型的反演方法多解性是基于模型地震反演的固有特性,即地震有效频带以多解性是基于模型地震反演的固有特性,即地震有效频带以外的信息不会影响合成地震资料的最终结果,减小基于模型外的信息不会影响合成地震资料的最终结果,减小基于模型方法多解性问题的关键在于正确建立初始模型方法多解性问题的关键在于正确建立初始模型 国内外主要反演方法国内外主要反演方法地震岩性模拟地震岩性模拟地震岩性模拟地震岩性模拟(SLIM)(SLIM)(SLIM)(SLIM) 递推反演递推反演递推反演递推反演 波恩反演波恩反演波恩反演波恩反演 人工神经网络反演人工神经网络反演人工神经网络反演人工神经网络反演 遗传算法反演遗传算法反演遗传算法反演遗传算法反演 混沌反演混沌反演混沌反演混沌反演 波阻抗多尺度反演波阻抗多尺度反演波阻抗多尺度反演波阻抗多尺度反演 多道反演多道反演多道反演多道反演 测井约束反演测井约束反演测井约束反演测井约束反演 地质统计学反演地质统计学反演地质统计学反演地质统计学反演 存在噪声存在噪声 假设条件与实际不符假设条件与实际不符 子波提取不精确子波提取不精确 测井资料的深时转换不准测井资料的深时转换不准 低频分量求取不准低频分量求取不准 约束条件不准确约束条件不准确 振幅、频率不保真振幅、频率不保真 5.1.2 叠后地震反演存在的问题 1 1、三高地震资料处理、三高地震资料处理 提高信噪比处理技术提高信噪比处理技术提高信噪比处理技术提高信噪比处理技术 – –提提高高信信噪噪比比处处理理分分为为叠叠前前和和叠叠后后两两个个环环节节,,叠叠前前去去噪噪方方法法的的选选取取取取决决于于噪噪声声类类型型及及其其在在各各种种道道集集上上的的表表现现形形式式,,处处理理原原则则是是::针针对对噪噪声声的的分分布布规规律律,,选选取取保保真真度度较较高高的的方方法法来来压压制制噪噪声声,,尽尽可可能能地地减减少少对对有有效效信信号号的的损损失失,,为为后后续续的的分分辨辨率率的的提高打好基础。

      提高打好基础振幅保真处理技术振幅保真处理技术振幅保真处理技术振幅保真处理技术 ––振振幅幅保保真真处处理理就就是是通通过过振振幅幅补补偿偿来来消消除除非非非非地地地地质质质质因因因因素素素素对对地地震震波波振振幅幅的的影影响响,,使使其其真真正正反反映映地地下下波波阻阻抗抗界界面面的的反反射射特特征征的的过过程振幅保真处理技术包括时间和空间两个方面的补偿技术程振幅保真处理技术包括时间和空间两个方面的补偿技术提高分辨率处理技术提高分辨率处理技术提高分辨率处理技术提高分辨率处理技术 ––保护低频、拓展高频是提高分辨率的唯一途径,压缩子波来保护低频、拓展高频是提高分辨率的唯一途径,压缩子波来扩展频宽的反褶积是地震数据处理中一个关键环节扩展频宽的反褶积是地震数据处理中一个关键环节,也是提,也是提高分辨率的主要方法高分辨率的主要方法5.1.3 反演的技术关键 三高地震资料处理的一般流程三高地震资料处理的一般流程三高地震资料处理的一般流程三高地震资料处理的一般流程 – – 1 1)解编;)解编;– – 2 2)球面扩散及吸收补偿;)球面扩散及吸收补偿;– – 3 3)高程静校正(初至折射静校正、微测井静校正);)高程静校正(初至折射静校正、微测井静校正);– – 4 4)叠前噪音衰减(小波变换压制面波);)叠前噪音衰减(小波变换压制面波);– – 5 5)地表一致性振幅补偿;)地表一致性振幅补偿;– – 6 6)组合反褶积;)组合反褶积;– – 7 7)高精度速度分析;)高精度速度分析;– – 8 8)多次剩余静校正;)多次剩余静校正;– – 9 9)多次动态时差校正;)多次动态时差校正;– – 1010)叠加;)叠加;– – 1111)叠后噪音衰减;)叠后噪音衰减;– – 1212)反褶积(预测反褶积、反)反褶积(预测反褶积、反滤波、子波零相位化等);滤波、子波零相位化等);– – 1313)偏移。

      偏移 1 1、三高地震资料处理、三高地震资料处理 5.1.3 反演的技术关键 2 2 2 2、测井资料的预处理及、测井资料的预处理及、测井资料的预处理及、测井资料的预处理及与地震资料的尺度匹配与地震资料的尺度匹配与地震资料的尺度匹配与地震资料的尺度匹配 测井数据的预处理测井数据的预处理 ––测井资料环境校正测井资料环境校正 ––测井资料的数据标准化测井资料的数据标准化 测井资料与地震资料的尺度匹配测井资料与地震资料的尺度匹配 ––尺度匹配和分层问题,一般采用中值滤波方法尺度匹配和分层问题,一般采用中值滤波方法 5.1.3 反演的技术关键 3 3、层位的精确标定和解释、层位的精确标定和解释 层位的精确标定层位的精确标定 ––平均速度曲线(深平均速度曲线(深——时转换尺)法;时转换尺)法;––合成地震记录拟合标定法;合成地震记录拟合标定法;––及垂直剖面法(及垂直剖面法(VSPVSP))标定标定 层位的解释层位的解释 5.1.3 反演的技术关键 3 3、层位的精确标定和解释、层位的精确标定和解释 正极性合成地震记录与井旁道对比图 负极性合成地震记录与井旁道对比图 5.1.3 反演的技术关键 4 4、子波处理、子波处理 时变子波的求取时变子波的求取 ––利用标定好的测井曲线、井旁地震道从浅到深分时利用标定好的测井曲线、井旁地震道从浅到深分时窗求取一组优化子波(时变子波)窗求取一组优化子波(时变子波)空变子波的求取空变子波的求取 ––内插法内插法––子波扫描法子波扫描法 ––基于模型的空变子波技术基于模型的空变子波技术 5.1.3 反演的技术关键 5 5、低频信息构建的经验和技巧、低频信息构建的经验和技巧 地震资料与测井资料的最佳匹配地震资料与测井资料的最佳匹配地震资料与测井资料的最佳匹配地震资料与测井资料的最佳匹配 测井资料的标准化测井资料的标准化测井资料的标准化测井资料的标准化 选择构造平缓方向构建低频分量选择构造平缓方向构建低频分量选择构造平缓方向构建低频分量选择构造平缓方向构建低频分量 岩性横向稳定的标志层作为起始层位岩性横向稳定的标志层作为起始层位岩性横向稳定的标志层作为起始层位岩性横向稳定的标志层作为起始层位 精细解释层位,确定低频分量内插、外推的权函数精细解释层位,确定低频分量内插、外推的权函数精细解释层位,确定低频分量内插、外推的权函数精细解释层位,确定低频分量内插、外推的权函数 考虑断层的空间走滑考虑断层的空间走滑考虑断层的空间走滑考虑断层的空间走滑 正确处理断层面处低频信息正确处理断层面处低频信息正确处理断层面处低频信息正确处理断层面处低频信息 考虑测线外井的贡献率考虑测线外井的贡献率考虑测线外井的贡献率考虑测线外井的贡献率 制作三维低频模型制作三维低频模型制作三维低频模型制作三维低频模型 确定岩性横向突变点确定岩性横向突变点确定岩性横向突变点确定岩性横向突变点 取长补短取长补短取长补短取长补短 注意高、低频带的衔接问题注意高、低频带的衔接问题注意高、低频带的衔接问题注意高、低频带的衔接问题 5.1.3 反演的技术关键 6 6、初始模型构建、初始模型构建 井旁初始波阻抗建立井旁初始波阻抗建立 构造构造框架模型框架模型的建立的建立 构建初始的波阻抗模型构建初始的波阻抗模型 5.1.3 反演的技术关键 6 6、初始模型构建、初始模型构建 反演的技术关键宏模型 微模型 构造构造框架模型框架模型的建立的建立 构建初始的波阻抗模型构建初始的波阻抗模型 6 6、初始模型构建、初始模型构建 反演的技术关键层序场层序场层位控制下的低频模型层位控制下的低频模型层序场控制下的低频模型层序场控制下的低频模型地震数据地震数据 通过古构造的恢复建造初始模型的方法示意图通过古构造的恢复建造初始模型的方法示意图 6 6、初始模型构建、初始模型构建 反演的技术关键 6 6、初始模型构建、初始模型构建 反演的技术关键通过古构造的恢复建造初始模型的方法示意图通过古构造的恢复建造初始模型的方法示意图 5.1.4 国内外主要反演软件美国的美国的美国的美国的SeislogSeislogSeislogSeislog反演软件;反演软件;反演软件;反演软件;丹麦丹麦丹麦丹麦ISISISISISISISIS模拟退火地震反演处理软件;模拟退火地震反演处理软件;模拟退火地震反演处理软件;模拟退火地震反演处理软件;荷兰荷兰荷兰荷兰JasonJasonJasonJason综合储层预测及描述软件;综合储层预测及描述软件;综合储层预测及描述软件;综合储层预测及描述软件;加拿大的加拿大的加拿大的加拿大的StrataStrataStrataStrata地震反演系统;地震反演系统;地震反演系统;地震反演系统;法国地球物理公司法国地球物理公司法国地球物理公司法国地球物理公司Strata VistaStrata VistaStrata VistaStrata Vista三维多道岩性反演软件;三维多道岩性反演软件;三维多道岩性反演软件;三维多道岩性反演软件;以色列帕拉代姆地球物理公司的以色列帕拉代姆地球物理公司的以色列帕拉代姆地球物理公司的以色列帕拉代姆地球物理公司的VanguagdVanguagdVanguagdVanguagd综合油藏描述系统;综合油藏描述系统;综合油藏描述系统;综合油藏描述系统;大庆油田勘探开发研究院的大庆油田勘探开发研究院的大庆油田勘探开发研究院的大庆油田勘探开发研究院的DelogIIDelogIIDelogIIDelogII反演软件;反演软件;反演软件;反演软件;石油勘探开发研究院计算机应用研究所的石油勘探开发研究院计算机应用研究所的石油勘探开发研究院计算机应用研究所的石油勘探开发研究院计算机应用研究所的RICHRICHRICHRICH油藏描述系统;油藏描述系统;油藏描述系统;油藏描述系统;石油勘探开发研究院西北分院的石油勘探开发研究院西北分院的石油勘探开发研究院西北分院的石油勘探开发研究院西北分院的CRISCRISCRISCRIS综合储层反演系统。

      综合储层反演系统综合储层反演系统综合储层反演系统 JASONJASONJASONJASON综合储层预测及描述软件综合储层预测及描述软件综合储层预测及描述软件综合储层预测及描述软件美国美国美国美国APEXAPEX技术公司技术公司技术公司技术公司 地震数据覆盖范围大,但分辨率较低;地震数据覆盖范围大,但分辨率较低;测井数据分辨率高,但是一孔之见测井数据分辨率高,但是一孔之见 JasonJason软件通过地层沉积模式的引入,将地质、地软件通过地层沉积模式的引入,将地质、地震、测井有机的结合起来,互相补充,互相震、测井有机的结合起来,互相补充,互相验证从而得到趋进真实的地质参数模型,验证从而得到趋进真实的地质参数模型,即储层预测和描述结果它主要共包括以下即储层预测和描述结果它主要共包括以下几个部分:几个部分: 地质框架模型地质框架模型(( EarthModEarthMod )) 子波分析子波分析(( Wavelets Wavelets )) 地震反演地震反演(( InverTraceInverTrace )) 测井反演测井反演(( InverMod InverMod )) 地质统计反演地质统计反演(( StatModStatMod ))概述概述 B B BB B BB B BA A AA A AA A AA A AC C CC C CC C C地质框架模型(地质框架模型(地质框架模型(地质框架模型(EarthModEarthMod))))地层框架表:定义单一层与其它层的接触关系,包括整合、断层、上超、地层框架表:定义单一层与其它层的接触关系,包括整合、断层、上超、底超、削截、礁体及河道。

      底超、削截、礁体及河道 地质框架模型(地质框架模型(地质框架模型(地质框架模型( EarthModEarthMod))))B B BB B BB B B实现过程:从建立的地质框架表出发,通过恢复剥蚀的层位、拉平地质层实现过程:从建立的地质框架表出发,通过恢复剥蚀的层位、拉平地质层位、井间插值、再恢复现今构造形态等位、井间插值、再恢复现今构造形态等 B B B地质框架模型(地质框架模型(地质框架模型(地质框架模型( EarthModEarthMod)))) 区域加权区域加权区域加权区域加权反向距离加权反向距离加权反向距离加权反向距离加权三角加权三角加权三角加权三角加权插值方式:插值方式:插值方式:插值方式:地质框架模型(地质框架模型(地质框架模型(地质框架模型(EarthModEarthMod)))) 地质框架模型(地质框架模型(地质框架模型(地质框架模型(EarthModEarthMod))))恢复现今状态恢复现今状态重构的结果重构的结果 ( (单层)单层)重扭曲的层位重扭曲的层位 ( (单层单层)))重新剥蚀的层位重新剥蚀的层位 ( (单层单层)) 地质模型的一个层位地质模型的一个层位恢复剥蚀的地层恢复剥蚀的地层恢复扭曲的地层恢复扭曲的地层地质框架模型(地质框架模型(地质框架模型(地质框架模型(EarthModEarthMod)))) 地质框架模型(地质框架模型(地质框架模型(地质框架模型(EarthModEarthMod))))恢复现今状态恢复现今状态重构的结果重构的结果 ( (单层)单层)重扭曲的层位重扭曲的层位 ( (单层)单层)重新剥蚀的层位重新剥蚀的层位 ( (单层单层)) 地震反演(地震反演(地震反演(地震反演(InverTraceInverTrace))))地震反演提供两种方法:地震反演提供两种方法:约束稀疏脉冲反演约束稀疏脉冲反演 建立在一个快速的趋势约束的脉冲反演算建立在一个快速的趋势约束的脉冲反演算法上。

      波阻抗趋势由解释的层位和井控制波阻抗趋势由解释的层位和井控制约束条件是波阻抗趋势加地质控制,产生一个把束条件是波阻抗趋势加地质控制,产生一个把地质模式融合进去的宽带结果地质模式融合进去的宽带结果递归反演递归反演 它是一种标准的生成地震带宽内的波阻抗它是一种标准的生成地震带宽内的波阻抗数据的递归反演算法数据的递归反演算法 地震反演流程地震反演流程地震反演流程地震反演流程1. 1. 1. 1. 从从从从层速度或测井资料出发,利用层速度或测井资料出发,利用层速度或测井资料出发,利用层速度或测井资料出发,利用地质框架模型,地质框架模型,地质框架模型,地质框架模型,得到低频模型得到低频模型得到低频模型得到低频模型2. 2. 2. 2. 井、层控制约束地震反演井、层控制约束地震反演井、层控制约束地震反演井、层控制约束地震反演得到得到得到得到相对相对相对相对波阻抗波阻抗波阻抗波阻抗3. 3. 3. 3. 低频模型低频模型低频模型低频模型+ + + +相对相对相对相对波阻抗得到绝对波阻抗数据体波阻抗得到绝对波阻抗数据体波阻抗得到绝对波阻抗数据体波阻抗得到绝对波阻抗数据体 约束条件约束条件约束条件约束条件 相对波阻抗相对波阻抗绝对波阻抗绝对波阻抗低频模型低频模型 地震反演结果地震反演结果地震反演结果地震反演结果((((InverTraceInverTrace)))) 测井反演(测井反演(测井反演(测井反演(InverMod InverMod InverMod InverMod )))) 以地质框架构建的地质模型作为基础与地震特征信息以地质框架构建的地质模型作为基础与地震特征信息相结合,生成一个精细的储层模型,该模型同时包含地震、相结合,生成一个精细的储层模型,该模型同时包含地震、测井和地质信息。

      测井和地质信息 该方法避开了褶积模型的局限性,不需要求取子波该方法避开了褶积模型的局限性,不需要求取子波而是采用地震特征进行反演而是采用地震特征进行反演具体算法是针对单层通过奇具体算法是针对单层通过奇异值分解算法(异值分解算法(SVDSVD))来优选出具有独立特征的测井信息来优选出具有独立特征的测井信息(主分量),针对每个反演道逐层求取每个测井信息的权(主分量),针对每个反演道逐层求取每个测井信息的权系数,最后,利用这些权系数对测井信息进行加权,得到系数,最后,利用这些权系数对测井信息进行加权,得到反演的波阻抗或其它测井信息反演的波阻抗或其它测井信息 测井反演流程测井反演流程测井反演流程测井反演流程1. 1. 1. 1. 建立水平层状地质模型建立水平层状地质模型建立水平层状地质模型建立水平层状地质模型2. 2. 2. 2. 产生井段产生井段产生井段产生井段3. 3. 3. 3. 主分量分析主分量分析主分量分析主分量分析 用奇异值分解的方法求解组合井段构成的矩阵用奇异值分解的方法求解组合井段构成的矩阵用奇异值分解的方法求解组合井段构成的矩阵用奇异值分解的方法求解组合井段构成的矩阵4. 4. 4. 4. 剔除无意义特征值对应的特征向量和权系数剔除无意义特征值对应的特征向量和权系数剔除无意义特征值对应的特征向量和权系数剔除无意义特征值对应的特征向量和权系数5. 5. 5. 5. 反演:通过解方程组得到权系数和厚度反演:通过解方程组得到权系数和厚度反演:通过解方程组得到权系数和厚度反演:通过解方程组得到权系数和厚度6. 6. 6. 6. 恢复地层现状恢复地层现状恢复地层现状恢复地层现状 建立水平层状地质模型建立水平层状地质模型建立水平层状地质模型建立水平层状地质模型地质模型的一个层位地质模型的一个层位地质模型的一个层位地质模型的一个层位恢复剥蚀的地层恢复剥蚀的地层恢复剥蚀的地层恢复剥蚀的地层恢复扭曲的地层恢复扭曲的地层恢复扭曲的地层恢复扭曲的地层 InverMod InverMod 算法算法算法算法生成组合井段(单层,单井)生成组合井段(单层,单井)生成组合井段(单层,单井)生成组合井段(单层,单井)PorosityPorosity ResistivityResistivityAIAISonicSonicDensityDensity波阻抗波阻抗 密度密度声波声波孔隙度孔隙度 电阻率电阻率 主分量分析主分量分析主分量分析主分量分析8 8 8 8 8 8口井的单层口井的单层口井的单层口井的单层口井的单层口井的单层奇异值分解奇异值分解奇异值分解针对所需要的某种参数对全区井资料单层段进行主分量分析针对所需要的某种参数对全区井资料单层段进行主分量分析针对所需要的某种参数对全区井资料单层段进行主分量分析针对所需要的某种参数对全区井资料单层段进行主分量分析 剔除无意义的特征值剔除无意义的特征值剔除无意义的特征值剔除无意义的特征值单层单层单层单层8 8 8 8个主分量个主分量个主分量个主分量 InverMod InverMod InverMod InverMod 算法算法算法算法初始权系数图初始权系数图初始权系数图初始权系数图, , , ,单层单层单层单层, 7, 7, 7, 7个个个个主分量主分量主分量主分量 InverMod InverMod InverMod InverMod 算法算法算法算法 s s11= f= f11 (w (w1A1A, , ww2A2A, , ww3A3A, , ..., ..., wwnAnA, , 层厚,主分层厚,主分层厚,主分层厚,主分量量量量))ss22= f= f22 (w (w1A1A, w, w2A2A, w, w3A3A,..., ,..., wwnAnA, , 层厚,主分层厚,主分层厚,主分层厚,主分量量量量)) . ... . .. .s smm= = f fm m (w(w1A1A, w, w2A2A, w, w3A,3A,..., ..., wwnAnA, , 层厚,主分层厚,主分层厚,主分层厚,主分量量量量))单道单层方程式单道单层方程式单道单层方程式单道单层方程式单道单层方程式单道单层方程式 ( ( ( ( ( (A)A)A)A)A)A) 用最小平方法求解方程用最小平方法求解方程用最小平方法求解方程用最小平方法求解方程 F(p)F(p)F(p)F(p)F(F(pp) = ) = FFglobalglobal((pp) + ) + FFsmoothsmooth((pp) + ) + FFseismicseismic((pp) + ) + FFwellswells((pp)) InverMod InverMod InverMod InverMod 算法算法算法算法求解出的权系数图求解出的权系数图求解出的权系数图求解出的权系数图求解出的权系数图求解出的权系数图, , , , , , 单层单层单层单层单层单层, 7, 7, 7, 7, 7, 7个个个个个个主分量主分量主分量主分量主分量主分量 InverModInverModInverModInverMod算法算法算法算法 PCPC1 1WW1 1 * *PCPC2 2WW2 2 * *+ + 重构模型重构模型重构模型重构模型重构模型重构模型 ( ( ( ( ( (单层,单道单层,单道单层,单道单层,单道单层,单道单层,单道) ) ) ) ) )= =Porosity, Porosity, etc.etc.+ +...................... .. 重构的结果重构的结果重构的结果重构的结果 ( ( ( (单层)单层)单层)单层)重扭曲的层位重扭曲的层位重扭曲的层位重扭曲的层位 ( ( ( (单层)单层)单层)单层)重新剥蚀的层位重新剥蚀的层位重新剥蚀的层位重新剥蚀的层位 ( ( ( (单层单层单层单层) )) )恢复现今状态恢复现今状态 测井反演结果(测井反演结果(测井反演结果(测井反演结果(InverModInverModInverModInverMod)))) 地地震震反反演演InverTraceInverTrace的频带的频带测测测测井井井井反反反反演演演演InverModInverMod的频带的频带的频带的频带 AmplitudeAmplitudeAmplitudeSeismic Seismic DataDataFrequencyFrequencyFrequency AmplitudeAmplitudeSeismic Seismic DataData 地质统计反演地质统计反演地质统计反演地质统计反演(((( StatModStatMod )))) 用地质统计学对非均质油藏进行用地质统计学对非均质油藏进行随机模拟,可用于滚随机模拟,可用于滚动勘探开发阶段,结果为孔隙度、渗透率、含油饱和度及动勘探开发阶段,结果为孔隙度、渗透率、含油饱和度及岩性分布。

      岩性分布随机模拟及随机反演分析(随机模拟及随机反演分析( StatmodStatmod-Analysis -Analysis ))随机模拟及随机反演建模随机模拟及随机反演建模 ((StatmodStatmod-Modeling -Modeling )) 反演结果精确反演结果精确反演结果精确反演结果精确地质统计模拟地质统计模拟地质统计模拟地质统计模拟地震反演地震反演地震反演地震反演地质统计反演地质统计反演地质统计反演地质统计反演地质统计反演地质统计反演 弹性反演(弹性反演(弹性反演(弹性反演(RockTraceRockTrace)))) 利用远偏移距或梯度偏移距资料进利用远偏移距或梯度偏移距资料进行地震反演,由于包括了转化波的信息,行地震反演,由于包括了转化波的信息,其结果不仅可以更好地区分岩性,而且其结果不仅可以更好地区分岩性,而且能够预测出流体能够预测出流体 声波阻抗声波阻抗弹性波阻抗弹性波阻抗 Gas Oil Water 砂岩声波阻抗砂岩声波阻抗砂砂砂砂岩岩岩岩弹弹弹弹性性性性波波波波阻阻阻阻抗抗抗抗交绘图交绘图 声波阻抗与频率分布的交绘图声波阻抗与频率分布的交绘图声波阻抗与频率分布的交绘图声波阻抗与频率分布的交绘图弹性波阻抗与频率分布的交绘图弹性波阻抗与频率分布的交绘图弹性波阻抗与频率分布的交绘图弹性波阻抗与频率分布的交绘图气气目标区目标区油油水水页岩页岩阻抗频率 气气目标区目标区油油水水页岩页岩频率阻抗 反演技术评价反演技术评价反演技术评价反演技术评价 JasonJason是一个综合性的地学软件,具有多是一个综合性的地学软件,具有多功能反演手段,通过地质框架模式的引入,将功能反演手段,通过地质框架模式的引入,将地质、地震、测井有机的结合起来,从而得到地质、地震、测井有机的结合起来,从而得到趋进真实的地质参数模型,即储层预测结果,趋进真实的地质参数模型,即储层预测结果,可为储层预测提供了多种地质信息。

      可为储层预测提供了多种地质信息 JasonJason软件可适应于勘探和开发阶段,特软件可适应于勘探和开发阶段,特别是开发阶段油气研究,但对于无井或井少地别是开发阶段油气研究,但对于无井或井少地区勘探阶段,区勘探阶段,JasonJason的主要功能难以发挥的主要功能难以发挥 该反演软件需要具有地质概念,地质框架该反演软件需要具有地质概念,地质框架模式的正确与否,直接影响反演的效果模式的正确与否,直接影响反演的效果 STRATA地震反演系统地震反演系统加拿大加拿大加拿大加拿大HAMPSON-RUSSELLHAMPSON-RUSSELL公司公司公司公司9/19/20249/19/2024 反演概述反演概述 Pre-Stack Inversion 迭前反演迭前反演 Travel Time Inversion (Tomography 旅行时层析成像旅行时层析成像反演反演) Amplitude Inversion ( (AVOAVO分析分析) ) Wave Field Inversion 波场反演波场反演 Bandlimited Inversion 带限反演带限反演 Sparse-spike Inversion 稀疏脉冲稀疏脉冲 Model-based Inversion 基于模型基于模型 Post-Stack Inversion 迭后反演迭后反演 Seismic Inversion Methods 地震反演方法地震反演方法 Amplitude Inversion (振幅反演振幅反演 P波波/S波反演波反演) NeuralNeuralNeural NetworkNetworkNetwork 神经网络神经网络神经网络神经网络神经网络神经网络 方法原理方法原理 StrataStrata地震地震反演的约束反演方法是一种基反演的约束反演方法是一种基于模型的波阻抗反演技术。

      这种方法利用测井于模型的波阻抗反演技术这种方法利用测井资料,以地震解释的层位为控制,从井点出发资料,以地震解释的层位为控制,从井点出发进行内插外推,形成初始波阻抗模型;然后利进行内插外推,形成初始波阻抗模型;然后利用共轭梯度法,对初始波阻抗模型不断进行更用共轭梯度法,对初始波阻抗模型不断进行更新,使得模型的合成记录与实际地新,使得模型的合成记录与实际地 震记录相似震记录相似性最好,此时的波阻抗模型便是反演结果性最好,此时的波阻抗模型便是反演结果 测井数据测井数据地震数据地震数据地震处理地震处理子波提取子波提取利用井提取利用井提取子波子波形成合成记录形成合成记录层位拾取层位拾取建模建模反演反演井相关井相关波阻抗波阻抗反反 演演 流流 程程 假设条件假设条件假设条件假设条件n n地震数据经叠后偏移处理,没有多次波存在地震数据经叠后偏移处理,没有多次波存在地震数据经叠后偏移处理,没有多次波存在地震数据经叠后偏移处理,没有多次波存在n n子波是常子波并且不随时间变化子波是常子波并且不随时间变化子波是常子波并且不随时间变化子波是常子波并且不随时间变化n n噪音是随机的,没有相干噪音噪音是随机的,没有相干噪音。

      噪音是随机的,没有相干噪音噪音是随机的,没有相干噪音数学模型数学模型数学模型数学模型 式式中M(i)为第i次迭代的解向量;G为灵敏度矩阵;S为实测数据 向量;S(i+1)为由第次迭代结果M(i)正演合成的的数据向量 地震处理功能地震处理功能 STRATA STRATA 还提供一系列地震处理功能:还提供一系列地震处理功能:还提供一系列地震处理功能:还提供一系列地震处理功能:n n 滤波滤波滤波滤波n n 反褶积反褶积反褶积反褶积n n 道组合道组合道组合道组合 n n 加权加权加权加权n n 相位转换相位转换相位转换相位转换 测井曲线标定测井曲线标定 根据时深关系将测井曲线做深时转换,根据时深关系将测井曲线做深时转换,在时间域与地震记录的解释层位对齐具体在时间域与地震记录的解释层位对齐具体标定时,可以在时间域和深度域内通过拉伸、标定时,可以在时间域和深度域内通过拉伸、压缩或移动曲线完成标定压缩或移动曲线完成标定 沿斜井轨迹提取地震记录和斜井曲线对沿斜井轨迹提取地震记录和斜井曲线对比标定 井相关井相关提供一系列井编辑功提供一系列井编辑功提供一系列井编辑功提供一系列井编辑功能能能能 。

      其中包括通过其中包括通过其中包括通过其中包括通过将合成记录与地震记将合成记录与地震记将合成记录与地震记将合成记录与地震记录拉齐,进而对井曲录拉齐,进而对井曲录拉齐,进而对井曲录拉齐,进而对井曲线进行拉伸和压缩线进行拉伸和压缩线进行拉伸和压缩线进行拉伸和压缩 子波提取子波提取 1. 1. 利用多道记录的自相关统计方法提取子利用多道记录的自相关统计方法提取子利用多道记录的自相关统计方法提取子利用多道记录的自相关统计方法提取子波振幅谱信息,进而求取子波波振幅谱信息,进而求取子波波振幅谱信息,进而求取子波波振幅谱信息,进而求取子波 2.2.利用测井资料和井旁地震数据,估计子利用测井资料和井旁地震数据,估计子利用测井资料和井旁地震数据,估计子利用测井资料和井旁地震数据,估计子波的振幅谱和相位谱,用最小平方法提取子波的振幅谱和相位谱,用最小平方法提取子波的振幅谱和相位谱,用最小平方法提取子波的振幅谱和相位谱,用最小平方法提取子波 3.3.从外部输入子波从外部输入子波从外部输入子波从外部输入子波 地震层位地震层位 地震层位是建模的基础地震层位是建模的基础地震层位是建模的基础地震层位是建模的基础,可以根据测井曲线标定的结果,,可以根据测井曲线标定的结果,,可以根据测井曲线标定的结果,,可以根据测井曲线标定的结果,在地震剖面上自动或人工拾取目的层位,在地震剖面上自动或人工拾取目的层位,在地震剖面上自动或人工拾取目的层位,在地震剖面上自动或人工拾取目的层位,也可以对拾取的层位也可以对拾取的层位也可以对拾取的层位也可以对拾取的层位进行删除或平滑处理。

      进行删除或平滑处理进行删除或平滑处理进行删除或平滑处理 初始波阻抗模型建立初始波阻抗模型建立 建立初始波阻抗模型的过程实际就是把横向上连续变化的地建立初始波阻抗模型的过程实际就是把横向上连续变化的地建立初始波阻抗模型的过程实际就是把横向上连续变化的地建立初始波阻抗模型的过程实际就是把横向上连续变化的地震界面信息与垂向上高分辨率的测井波阻抗信息相结合的过程震界面信息与垂向上高分辨率的测井波阻抗信息相结合的过程震界面信息与垂向上高分辨率的测井波阻抗信息相结合的过程震界面信息与垂向上高分辨率的测井波阻抗信息相结合的过程 井旁初始波阻抗建立:井旁初始波阻抗建立:井旁初始波阻抗建立:井旁初始波阻抗建立: 首先通过井旁地震道与合成记录的相关性对测井曲线进行纵首先通过井旁地震道与合成记录的相关性对测井曲线进行纵首先通过井旁地震道与合成记录的相关性对测井曲线进行纵首先通过井旁地震道与合成记录的相关性对测井曲线进行纵向的拉伸和压缩,当相关系数达到一定的标准时,就可以获得井向的拉伸和压缩,当相关系数达到一定的标准时,就可以获得井向的拉伸和压缩,当相关系数达到一定的标准时,就可以获得井向的拉伸和压缩,当相关系数达到一定的标准时,就可以获得井的初始波阻抗。

      的初始波阻抗的初始波阻抗的初始波阻抗 全局初始波阻抗模型建立:全局初始波阻抗模型建立:全局初始波阻抗模型建立:全局初始波阻抗模型建立: 在地震层位和地质模式的约束下,选取适当的插值方法,对在地震层位和地质模式的约束下,选取适当的插值方法,对在地震层位和地质模式的约束下,选取适当的插值方法,对在地震层位和地质模式的约束下,选取适当的插值方法,对井的初始波阻抗进行内插和外推,建立全局初始波阻抗模型井的初始波阻抗进行内插和外推,建立全局初始波阻抗模型井的初始波阻抗进行内插和外推,建立全局初始波阻抗模型井的初始波阻抗进行内插和外推,建立全局初始波阻抗模型 整合接触整合接触 削蚀削蚀 上超上超地质模式地质模式 初始模型插值方法初始模型插值方法 Well 2Well 2 WellWell 1 1 Well 3Well 3反距离开方反距离开方反距离开方反距离开方反距离平方反距离平方反距离平方反距离平方克里金克里金克里金克里金三角内插三角内插三角内插三角内插用户定义用户定义用户定义用户定义 波阻抗反演波阻抗反演 以初始波阻抗模型为基础,根据地震记录与合成记录的相以初始波阻抗模型为基础,根据地震记录与合成记录的相以初始波阻抗模型为基础,根据地震记录与合成记录的相以初始波阻抗模型为基础,根据地震记录与合成记录的相关性,选取反演算法,对初始波阻抗模型进行优化,就可以获关性,选取反演算法,对初始波阻抗模型进行优化,就可以获关性,选取反演算法,对初始波阻抗模型进行优化,就可以获关性,选取反演算法,对初始波阻抗模型进行优化,就可以获得最终的波阻抗。

      适用于二维或三维工区得最终的波阻抗适用于二维或三维工区得最终的波阻抗适用于二维或三维工区得最终的波阻抗适用于二维或三维工区 四种不同的反演算法:四种不同的反演算法:四种不同的反演算法:四种不同的反演算法:n n 带限反演带限反演带限反演带限反演n n 稀疏脉冲反演稀疏脉冲反演稀疏脉冲反演稀疏脉冲反演n n约束反演(基于模型反演)约束反演(基于模型反演)约束反演(基于模型反演)约束反演(基于模型反演)n n 人工智能神经网络反演人工智能神经网络反演人工智能神经网络反演人工智能神经网络反演 带限反演带限反演 带限反演是通过地带限反演是通过地层反褶积得到带限的反层反褶积得到带限的反射系数序列,并对反射射系数序列,并对反射系数序列递推计算出地系数序列递推计算出地层波阻抗的地震反演方层波阻抗的地震反演方法 最大相似稀疏脉冲反演最大相似稀疏脉冲反演 稀疏脉冲反演是基于反褶积稀疏脉冲反演是基于反褶积算法的一种反演方法对于每一算法的一种反演方法对于每一道,通过每次增加一个反射系数道,通过每次增加一个反射系数来优化反射系数序列,直到所得来优化反射系数序列,直到所得的合成记录与地震记录达到允许的合成记录与地震记录达到允许误差为止,此时就得到最佳反射误差为止,此时就得到最佳反射序列的估算,进而递推计算出波序列的估算,进而递推计算出波阻抗。

      阻抗 约束反演约束反演约束反演是从初始波阻抗模型约束反演是从初始波阻抗模型约束反演是从初始波阻抗模型约束反演是从初始波阻抗模型出发,利用迭代算法优化初始出发,利用迭代算法优化初始出发,利用迭代算法优化初始出发,利用迭代算法优化初始波阻抗模型,当地震道与合成波阻抗模型,当地震道与合成波阻抗模型,当地震道与合成波阻抗模型,当地震道与合成记录的相关系数达到标准时,记录的相关系数达到标准时,记录的相关系数达到标准时,记录的相关系数达到标准时,改进的模型即为反演结果改进的模型即为反演结果改进的模型即为反演结果改进的模型即为反演结果 人工智能神经网络反演人工智能神经网络反演 人工智能神经网络反演是利用神经网络人工智能神经网络反演是利用神经网络方法,分析测井数据和井旁地震数据的关系方法,分析测井数据和井旁地震数据的关系对这种关系不做预先假定,关系可能是复杂对这种关系不做预先假定,关系可能是复杂的关系,也可能是非线性的神经网络反演的关系,也可能是非线性的神经网络反演法在目标区内有足够井的情况下是一种理想法在目标区内有足够井的情况下是一种理想的方法 反演结果诊断反演结果诊断n n 利用提取的子波和获得的反射系数得到合成地震记录利用提取的子波和获得的反射系数得到合成地震记录利用提取的子波和获得的反射系数得到合成地震记录利用提取的子波和获得的反射系数得到合成地震记录n n 通过误差剖面可以显示原始地震记录与合成记录的差别通过误差剖面可以显示原始地震记录与合成记录的差别通过误差剖面可以显示原始地震记录与合成记录的差别通过误差剖面可以显示原始地震记录与合成记录的差别 反演结果诊断反演结果诊断n n 波阻抗,速度,合成记录和误差可重叠显示波阻抗,速度,合成记录和误差可重叠显示波阻抗,速度,合成记录和误差可重叠显示波阻抗,速度,合成记录和误差可重叠显示n n 在两个窗口中可以比较多个反演结果在两个窗口中可以比较多个反演结果在两个窗口中可以比较多个反演结果在两个窗口中可以比较多个反演结果 AVO分析与反演分析与反演 生成不同偏移距的剖面生成不同偏移距的剖面生成不同偏移距的剖面生成不同偏移距的剖面 进行进行进行进行3-3-D AVOD AVO分析处理分析处理分析处理分析处理 利用利用利用利用AVOAVO剖面进行反演剖面进行反演剖面进行反演剖面进行反演 AVO分析分析 不同偏移距的叠加剖面不同偏移距的叠加剖面不同偏移距的叠加剖面不同偏移距的叠加剖面. . P P P P波反演的波阻抗波反演的波阻抗波反演的波阻抗波反演的波阻抗 S S波反演的波阻抗波反演的波阻抗波反演的波阻抗波反演的波阻抗 反演技术评价反演技术评价 STARTA软件是一套较完整的地震反演软件。

      该软件提软件是一套较完整的地震反演软件该软件提供了四种反演方法:带限反演、稀疏脉冲反演、基于模型约供了四种反演方法:带限反演、稀疏脉冲反演、基于模型约束反演、人工智能神经网络反演束反演、人工智能神经网络反演 目前大家常用是基于模型约束反演,该方法把地震和测井目前大家常用是基于模型约束反演,该方法把地震和测井有机的结合起来,突破了传统意义上的地震分辨率的限制,有机的结合起来,突破了传统意义上的地震分辨率的限制,反演结果既有地震信息,亦有测井信息,适应于薄储层预测,反演结果既有地震信息,亦有测井信息,适应于薄储层预测,可应用于勘探和开发阶段的精细油藏描述可应用于勘探和开发阶段的精细油藏描述 多解性是该方法的固有特性,主要取决于初始模型与实多解性是该方法的固有特性,主要取决于初始模型与实际地质情况的符合程度在反演中要注意解决好两方面问题:际地质情况的符合程度在反演中要注意解决好两方面问题:其一是提供正确地质模型来指导测井资料内插外推建立初始其一是提供正确地质模型来指导测井资料内插外推建立初始模型;其二是提高地震资料的分辨率模型;其二是提高地震资料的分辨率 第五章5.1 叠后波阻抗反演叠后波阻抗反演5.2 叠后地震属性叠后地震属性5.3 利用地震参数进行储层预测的主要方利用地震参数进行储层预测的主要方法法 5.2 叠后地震属性叠后地震属性•5.2.1 地震特征参数的类型及意义地震特征参数的类型及意义•5.2.2 地震地震波形属性分析波形属性分析 •5.2.3 相干体属性相干体属性 •5.2.4 频谱分解频谱分解 •5.2.5 亮点和亮点和“AVO”属性属性 •从20世纪70年代开始研究和使用地震属性,最早主要是以振幅为基础的瞬时属性,用来直接指示油气。

      •80年代,地震属性的数量迅速增加,多元属性分析技术首次被使用•90年代地震属性技术的蓬勃发展同样得益于全三维解释技术的发展 沿层提取的地震参数一览表沿层提取的地震参数一览表 振幅统计类振幅统计类:均方根振幅平均绝对振幅最大波峰振幅平均波峰振幅最大波谷振幅平均波谷振幅最大绝对振幅总绝对振幅平均能量总能量平均振幅振幅变差振幅偏斜度振幅峰态复数道统计类复数道统计类:平均反射强度平均瞬时频率平均瞬时相位反射强度梯度瞬时频率梯度频谱统计类频谱统计类:有效带宽有效带宽波形弧长波形弧长平均过零频度平均过零频度优势频率优势频率峰谱频率峰谱频率频谱峰最梯度频谱峰最梯度相关统计类相关统计类:相邻道协方差相邻道协方差相邻道相关时移相邻道相关时移平均信噪比平均信噪比相关长度相关长度相关分量相关分量卡拉信号复杂度卡拉信号复杂度吸收系数类吸收系数类:吸收系数吸收系数 构造变形类构造变形类 倾角倾角方位角方位角平均曲率平均曲率粗度粗度近反射层几何近反射层几何形状形状 层序统计类:层序统计类:过门槛百分比过门槛百分比欠门槛百分比欠门槛百分比能量过半时长能量过半时长能量半时长梯度能量半时长梯度正负样数比正负样数比波峰数波峰数波谷数波谷数一、地震特征参数的类型及意义一、地震特征参数的类型及意义 地震特征参数的类型及意义地震特征参数的类型及意义 地震属性或指示特征可能反映的地质现象或特征相关参数振幅(瞬时+能量)古地貌、岩性差异、岩层连续性、总孔隙度视极性(瞬时+能量)岩性、反射极性差异,含气性频率(瞬时+能量)岩层厚度及流体性质相位(瞬时+能量)岩层连续性、地层结构振幅极小值与极大值数目比及位置古地貌、岩相结构层速度岩性、孔隙度、压力体反射谱分解的各阶分量横向、垂向分辩率,孔隙度、流体及几何形态AVO岩层中流体性质声阻抗孔隙度及泥质含量曲率、边界增强等现象断层及裂缝分布特征倾角、方位角及人工照明等处理成果构造、断层及由地震资料处理得到的地质特征烃指示属性均方根、最小振幅、最大振幅、最大振幅绝对值、波峰平均值岩性、物性指示属性波谷平均值、平均能量、振幅和、振幅绝对值之和频率—烃类指示属性优势频率、平均瞬时频率频率—岩性、物性指示属性半幅能量、门槛值流体指示属性平均瞬时相位岩相水平、垂向变化特征零值个数、弧长、带宽地震属性可能反映的储层性质地震属性可能反映的储层性质 5.2 叠后地震属性叠后地震属性•5.2.1 地震特征参数的类型及意义地震特征参数的类型及意义•5.2.2 地震地震波形属性分析波形属性分析 •5.2.3 相干体属性相干体属性 •5.2.4 频谱分解频谱分解 •5.2.5 亮点和亮点和“AVO”属性属性 Stratimagic分析原理分析原理 通过对实际地震体层位解释,建立等厚或非等厚的层段,运用神经网络识别技术,建立层段内理论模型道序列(不同频率的波形特征),最终在平面或剖面上显示不同波形的展布,如果研究层段内出现有异常体,其层段波形因其波形特征不同将同样出现异常,在平面或剖面中通过颜色的变化也将得以体现。

      地震地震波形属性分析波形属性分析 地震层段的建立地震地震波形属性分析波形属性分析 使用神经网络对道使用神经网络对道形状进行分类形状进行分类模型道的建立地震地震波形属性分析波形属性分析 一旦形成了相图,即可执行分类 该过程从层段中提取每一个实际地震道和模型道进行相关,并为其分配与模型道有最佳相关的相编码(色彩)地震地震波形属性分析波形属性分析 层段的波形特征在剖面上的响应层段的波形特征在剖面上的响应地震地震波形属性分析波形属性分析 地震地震波形属性分析波形属性分析 5.2 叠后地震属性叠后地震属性•5.2.1 地震特征参数的类型及意义地震特征参数的类型及意义•5.2.2 地震地震波形属性分析波形属性分析 •5.2.3 相干体属性相干体属性 •5.2.4 频谱分解频谱分解 •5.2.5 亮点和亮点和“AVO”属性属性 •相干概念最初只限于相邻道间的相关性,是多道数据间相似程度的一种度量它最有效和最重要的应用之一是同相轴的检测,特别是当同相轴的幅值较小并隐藏在噪声中时 •用相干体技术进行断层解释和组合,避免了解释结果的随意性,补充了因解释人员经验所导致的结果不合理性 相干体属性相干体属性 相干体解释的过程与步骤相干体解释的过程与步骤 •把相干数据体加入到三维可视化系统中 •动画浏览。

      •地震相干数据体的时间切片•检查合理性•把层位解释结果输入,沿层进行相干数据体断层平面组合•把空间断层解释结果输入并加入到可视化的常规地震数据体中,为下步工作和出图等做准备相干体属性相干体属性 相干属性的优势相干属性的优势•相干体技术可分辨出在任意方位上的断层,有助于更精确和更有效的解释 •在描述断层和隐藏构造时,相干体时间切片比常规三维地震数据体切片要清晰地多 •对相干数据体和三维可视化的分析能准确标出断层和地层构造的空间位置,对井位设计大有帮助 •工作效率大大提高,解释结果更客观、合理 相干体属性相干体属性 相干体技术应注意相干体技术应注意的问题的问题•对地震资料要进行评估,只有信噪比高的资料解释出的断层才可信,否则很难对相干进行解释•不相干数据异常不一定都是断层,也可能是因岩性变化或其它地质现象所导致,所以在解释时结合垂直剖面进行具体分析,不要一概而论•相干数据体不是对所有的断层都能识别,当垂直落差在一个视周期或它的整数倍时,在相干体中断层没有反映相干体属性相干体属性 剖面与平面的交互对比解释手段:克服了剖面与平面的交互对比解释手段:克服了断层的多解性,极大地提高了复杂地区断断层的多解性,极大地提高了复杂地区断层的解释精度。

      层的解释精度解释过程中:解释过程中:断层解释断层解释——地震属性分析断层解释地震属性分析断层解释相干体属性相干体属性 1336ms 莱莱711360ms 莱莱71相相 干干 体体 发发 现现 了了 东东 营营 组组 曲曲 流流 河河 砂砂 体体Crossline364Crossline364地震剖面地震剖面东部油田曲流河东部油田曲流河三、相干体属性三、相干体属性 ————地震属性分析地震属性分析 采用相干体切采用相干体切片的方法进行裂缝片的方法进行裂缝油气藏的预测油气藏的预测 其中在相干体其中在相干体中切片暗色代表相中切片暗色代表相干性较差,可以认干性较差,可以认为是断层或裂缝发为是断层或裂缝发育区0ms4ms8ms12ms68ms72ms76ms80ms84ms88ms92ms96ms100ms海坨地区:海坨地区:海坨地区:海坨地区:相干体属性相干体属性 5.2 叠后地震属性叠后地震属性•5.2.1 地震特征参数的类型及意义地震特征参数的类型及意义•5.2.2 地震地震波形属性分析波形属性分析 •5.2.3 相干体属性相干体属性 •5.2.4 频谱分解频谱分解 •5.2.5 亮点和亮点和“AVO”属性属性 •在大型3-D地震勘探中,在利用地震数据和其离散傅里叶变换(DFT)对时间地层厚度和地质不连续体进行成像和作图时,频谱分解为之提供了一种新颖的方法。

      经DFT把地震数据变换到频率域,振幅谱描绘时间地层厚度的变化,而相位谱则显示了地质体的横向不连续性在3—D地震勘探中,用该技术已经确定了地层的形成环境(如河道砂和涉及复杂断层系统的构造环境)•频谱分解基于来自薄层的反射在频率域具有指示时间地层厚度的特征性表现 频谱分解频谱分解 薄层的频谱成像(据Partyka et al, 1999) 频谱分解频谱分解 长时窗的频谱分解及其与褶积模型的关系(据Partyka et al, 1999) 频谱分解频谱分解 短时窗的频谱分解及其与褶积模型的关系(据Partyka et al, 1999) 频谱分解频谱分解 块状楔形模型(据Partyka et al, 1999)(a)反射率;(b)滤波后的反射;(c)谱振幅 频谱分解频谱分解楔形模型响应楔形模型响应 关于频率(a)及关于薄层厚度(b)的振幅对频率的薄层调谐(据Partyka et al, 1999) 频谱分解频谱分解楔形模型响应楔形模型响应 调谐立方体调谐立方体 频谱分解频谱分解目标层调谐立方体(据Partyka et al, 1999) 消除子波叠复消除子波叠复 频谱分解频谱分解在均衡频谱之前,调谐立方体是由薄层干涉、地震子波和随机噪声组成的(据Partyka et al, 1999) 消除了子波叠复(均衡频谱)而没有消除反射率调谐特征的示意图(据Partyka et al, 1999) 频谱分解频谱分解消除子波叠复消除子波叠复 目标层以外目标层以外 频谱分解频谱分解离散频率能量立方体(据Partyka et al, 1999) 安哥拉海上某深水油田采用频谱分解技术所成的地下3000ft目的层图像 频谱分解频谱分解 大庆黑帝庙油层H2-2砂组50Hz振幅谱及其地质解释 频谱分解频谱分解 5.2 叠后地震属性叠后地震属性•5.2.1 地震特征参数的类型及意义地震特征参数的类型及意义•5.2.2 地震地震波形属性分析波形属性分析 •5.2.3 相干体属性相干体属性 •5.2.4 频谱分解频谱分解 •5.2.5 亮点和亮点和“AVO”属性属性 ““亮点亮点””技术技术 •当岩石孔隙中含油气时特别是含气时,其纵波波速将比含水时要低,因此其小阻抗将比含水时低,特别在含气时,其纵波波阻抗将明显降低,使含气砂岩与页岩分界面的反射系数可达0.29。

      这远大于一般反射系数,将产生很强的反射振幅,这就是所谓“亮点技术”的物理基础 五、亮点和五、亮点和“AVO”属性属性 AVO技术技术 Koefoed 50年代年代在计算中假设界面两侧的泊松比是不同在计算中假设界面两侧的泊松比是不同的,所得纵波反射系数随入射角的改变而改变的范围大大增加,的,所得纵波反射系数随入射角的改变而改变的范围大大增加,从而得到了一系列对从而得到了一系列对AVO技术发展有重大意义的结论技术发展有重大意义的结论 ,奠定了奠定了AVO技术理论基础技术理论基础 80年代年代Shuey等人简化了等人简化了Zoeppritz方程,把方程,把AVO技术技术投入到实际应用中,但应用的单一性:利用反射振幅随炮检距投入到实际应用中,但应用的单一性:利用反射振幅随炮检距增加而增加特征直接检测气层(识别真假亮点,实际既有亮点,增加而增加特征直接检测气层(识别真假亮点,实际既有亮点,又有暗点存在)又有暗点存在) 80年代中期年代中期--90年代年代—至今:至今:三分量三分量AVO技术、纵横波技术、纵横波联合定量反演等运用联合定量反演等运用AVO技术检测裂缝,研究地下各向异性技术检测裂缝,研究地下各向异性等。

      技术的进步:等技术的进步:AVO处理技术,处理技术,AVO正演综合解释技术等正演综合解释技术等五、亮点和五、亮点和“AVO”属性属性 地震波的速度地震波的速度第一各岩石速度和波阻抗值都有一定分布范围第一各岩石速度和波阻抗值都有一定分布范围第二各岩石速度和波阻抗值又都有一定的重叠第二各岩石速度和波阻抗值又都有一定的重叠五、亮点和五、亮点和“AVO”属性属性AVO技术地质基础技术地质基础 各种岩石的纵波速度各种岩石的纵波速度m/s五、亮点和五、亮点和“AVO”属性属性AVO技术地质基础技术地质基础 各种岩石的波阻抗各种岩石的波阻抗五、亮点和五、亮点和“AVO”属性属性AVO技术地质基础技术地质基础 各种岩石的泊松比各种岩石的泊松比一般沉积岩石的泊松比值大于一般沉积岩石的泊松比值大于0.2,含气砂岩的泊松比值一般小于,含气砂岩的泊松比值一般小于0.2 五、亮点和五、亮点和“AVO”属性属性AVO技术地质基础技术地质基础   沉积岩纵波速度与泊松比关系图沉积岩纵波速度与泊松比关系图五、亮点和五、亮点和“AVO”属性属性AVO技术地质基础技术地质基础 石英含量石英含量 泊松比泊松比 泥质含量泥质含量 泊松比泊松比碳酸盐岩的孔隙度碳酸盐岩的孔隙度 泊松比泊松比含气砂岩的孔隙度含气砂岩的孔隙度 泊松比泊松比含油或含水砂岩,孔隙度含油或含水砂岩,孔隙度 泊松比泊松比饱和气的岩石,孔隙的纵横波比饱和气的岩石,孔隙的纵横波比 泊松比泊松比含水、含气和含中含水、含气和含中—低气油比油的岩石,孔隙低气油比油的岩石,孔隙 压力压力 泊松比泊松比含高气油比油的岩石,当孔隙压力大大超过其含高气油比油的岩石,当孔隙压力大大超过其 泡点压力,随孔隙压力泡点压力,随孔隙压力 泊松比泊松比围压增围压增 泊松比泊松比温度升温度升 泊松比泊松比影响岩石泊松比的主要因素:影响岩石泊松比的主要因素: 岩石的成分,储层中的流体性质等岩石的成分,储层中的流体性质等 Por交会图交会图GR、、DEN、、AC与与POR交会图交会图W1w2W1w2W1w2AVO资料解释资料解释 对对AVO技术而言,远道叠加剖面对气层分辨率提高技术而言,远道叠加剖面对气层分辨率提高1.随炮检距增加能量增强,反射系数增加,.随炮检距增加能量增强,反射系数增加, 能量增强;能量增强;2.由于目的层埋藏较深,近远道叠加剖面频率相当.由于目的层埋藏较深,近远道叠加剖面频率相当 没有降低分辨率;没有降低分辨率;3.远道比近道叠加剖面信噪比高。

      .远道比近道叠加剖面信噪比高 AVO解释中的陷阱解释中的陷阱 在在AVO正演中,特别要研究储层中正演中,特别要研究储层中不同流体性质的变化,不同流体性质的变化, 研究含气层与背景值的偏离,研究含气层与背景值的偏离, 区别含气层与非含气层产生的区别含气层与非含气层产生的AVO现象的异同,现象的异同, 减小减小AVO气层检测中的误判气层检测中的误判 第五章•5.1 叠后波阻抗反演叠后波阻抗反演•5.2 叠后地震属性叠后地震属性•5.3 利用地震参数进行储层预测的主要方利用地震参数进行储层预测的主要方法法 利用地震参数进行储层预测的主要方法利用地震参数进行储层预测的主要方法•一、基于波阻抗与地质参数间物理关系一、基于波阻抗与地质参数间物理关系的储层预测的储层预测 •二、基于地质体几何特征的单一属性空二、基于地质体几何特征的单一属性空间形态解释间形态解释 •三、基于多地震属性模式识别的储层预三、基于多地震属性模式识别的储层预测测 •四、基于统计学随机建模和正演模拟的四、基于统计学随机建模和正演模拟的储层预测储层预测 •波阻抗与储层参数及其它地质因素间的关系波阻抗与储层参数及其它地质因素间的关系 –储层参数主要指岩性、物性和含油气性。

      它们都对岩石的波阻抗有显著的影响,这也正是可望通过波阻抗进行储层参数预测的根本出发点–除了上述储层参数外,还有许多地质因素对岩石的波阻抗有影响,它们干扰着从波阻抗出发进行储层参数预测,因此需要采用适当措施对其排除–波阻抗是速度与密度的乘积,因此要分析其与储层参数以及其它地质因素之间的关系,应当从速度与密度这两个基本要素的影响因素分析入手 基于波阻抗与地质参数间物理关系的储层预测基于波阻抗与地质参数间物理关系的储层预测 •影响岩石速度、密度和波阻抗的控制因影响岩石速度、密度和波阻抗的控制因素的层次分析素的层次分析 –孔隙度,流体的饱和度和粘度(表现为耦合系数)它们属于起直接控制作用的第一层次的因素 –第二层次的因素通过对第一层次因素的控制作用而间接影响着岩石速度 –三层次的因素:它们直接影响着第二层次的因素,与岩石速度间的关系就更为间接如岩性主要取决于沉积条件,包括沉积环境、物源性质、古气候条件等等 基于波阻抗与地质参数间物理关系的储层预测基于波阻抗与地质参数间物理关系的储层预测 储层参数预测方法储层参数预测方法 Vr=f(Vb、Vf、Ф、q、岩性、成岩程度、流体性质、温度压力、沉积条件、成岩环境、地层时代、深度)一)横向相对变化半定量分析法一)横向相对变化半定量分析法 一)横向相对变化半定量分析法一)横向相对变化半定量分析法 三)单参数解释中的三)单参数解释中的““地质消元地质消元””法法 基于波阻抗与地质参数间物理关系的储层预测基于波阻抗与地质参数间物理关系的储层预测 利用波阻抗进行储层预测的常用方法利用波阻抗进行储层预测的常用方法 •对波阻抗剖面调色进行定性岩性解释对波阻抗剖面调色进行定性岩性解释 •根据波阻抗根据波阻抗- -深度深度- -岩性量板进行定量岩岩性量板进行定量岩性解释性解释 •根据物理关系进行储层泥质含量、孔隙根据物理关系进行储层泥质含量、孔隙度等的定量预测度等的定量预测 •地质统计学方法地质统计学方法 基于波阻抗与地质参数间物理关系的储层预测基于波阻抗与地质参数间物理关系的储层预测 对波阻抗剖面调色进行定性岩性解释对波阻抗剖面调色进行定性岩性解释 对波阻抗剖面调色进行定性岩性解释对波阻抗剖面调色进行定性岩性解释 文昌9-1-1井、9-2-1井波阻抗-GR-深度交会图对比 a. Jasona. Jason反演反演 樊131井沙二下段2200-2400m波阻抗曲线与伽玛曲线交会图 交会点的颜色为密度值,由黄交会点的颜色为密度值,由黄→深黄深黄→红红→绿绿→浅兰浅兰→蓝色,其值逐渐减小,即由砂岩向泥岩区过渡。

      蓝色,其值逐渐减小,即由砂岩向泥岩区过渡 砂岩区砂岩区泥岩区泥岩区伽玛值增大伽玛值增大波阻抗值增大波阻抗值增大根据波阻抗根据波阻抗- -深度深度- -岩性量板进行定量岩性解释岩性量板进行定量岩性解释 根据波阻抗根据波阻抗- -深度深度- -岩性量板进行定量岩性解释岩性量板进行定量岩性解释玛湖地区波阻抗-岩性量版图 •经验公式法经验公式法 –经验公式法求取泥质含量常用的是时间平均方程,如下式所示:–求取孔隙度的时间平均方程为: 根据物理关系进行储层泥质含量、孔隙度等的定量预测根据物理关系进行储层泥质含量、孔隙度等的定量预测 •数理统计法数理统计法 –目前较为流行,一般是利用已知泥质含量或孔隙度与层速度(或波阻抗),进行统计分析回归出某一关系式,然后用此关系式将地震求取的层速度(或波阻抗),转化为泥质含量或孔隙度,达到参数预测的目的 –根据研究区的声波测井资料,设计一速度——岩性间的函数关系: 根据物理关系进行储层泥质含量、孔隙度等的定量预测根据物理关系进行储层泥质含量、孔隙度等的定量预测 •该方法在储层参数预测时,确保预测结果在井点处准确无误,在井点以外则根据测井资料、地震资料、测井地震资料的空间变化规律、参数的分布概率等进行储层参数的预测。

      •地质统计学方法以其用地质规律为约束的先进插值算法和等概率多实现的方式,在对空间分布变量的分析预测中显得独具优势 地质统计学方法地质统计学方法 基于地质体几何特征的单一基于地质体几何特征的单一属性空间形态解释属性空间形态解释 •空间切片分析空间切片分析–地震波振幅切片、相干体切片 、属性类切片、统计类切片 •地震切片的应用地震切片的应用–水平切片 、相干属性类切片 •三维可视化技术三维可视化技术 –体内层位自动追踪、目的层面可视化、等时体可视化、层面与监控剖面可视化及多属性体叠合可视化 基于地质体几何特征的单一属性空间形态解释基于地质体几何特征的单一属性空间形态解释 燕南潜山带燕南潜山带燕南潜山带燕南潜山带L850L850线典型解释剖面线典型解释剖面线典型解释剖面线典型解释剖面 馆陶组中层拉平最大波谷振幅图馆陶组下部层拉平零相位平均频率图StratmigicStratmigic地震相分析图地震相分析图地震相分析图地震相分析图22%33%44%55%67%78%89%100%98318-1-122%33%44%55%67%78%89%100%98318-1-122%33%44%55%67%78%89%100%983上30ms 下10ms馆陶组下部层拉平最大波谷振幅图物源方向物源方向物源方向物源方向 燕南潜山带东二上段底燕南潜山带东二上段底最最 大大 波波 谷谷 振振 幅幅 图图储层预测储层预测(2)东营组 空间切片分析空间切片分析振幅的垂直切片(常规地震剖面) 走滑断层燕东断层燕南断层燕南断层燕东断层走滑断层燕南水平时间切片燕南水平时间切片燕南水平时间切片燕南水平时间切片构造解释构造解释振幅的水平切片(时间切片) 空间切片分析空间切片分析 任意线(面)切片任意线(面)切片 空间切片分析空间切片分析文昌A凹陷层序地层综合连井剖面图 地震相位切片地震相位切片5号断层1号断层3号断层 统计类切片统计类切片 空间切片分析空间切片分析波峰波谷平均振幅图 波峰-波谷比值图 9-1-1井南断层连通性2365ms2375ms2405ms2415ms相干体切片 地震切片的应用地震切片的应用 •水平切片水平切片 ①指导、验证构造解释的合理性和精度;②根据同相轴的中断、错开、振幅突变、拐弯、走向不一致等特点,分析小断裂和断层的细微变化,提高圈闭落实程度和可靠性;③应用叠置的水平切片,可以进行圈闭识别和微幅度构造精细解释;④“剖面定倾向、切片定走向”,强化构造解释的精度;⑤依据同相轴振幅、频率、相位等变化特征,进行地震品质评价。

      •相干属性类切片相干属性类切片 ①利用垂直相干切片可完成断层的自动和半自动解释②应用相干切片,研究断裂、断层空间展布规律,指导和验证平面断裂组合,解释在传统水平切片上显示不清的断层③应用沿层相干切片,分析、研究小断裂和小断层空间的细微变化,指导和验证平面断裂组合④应用沿层相干切片或沿层属性切片进行地震相或沉积相的研究和识别 地震切片的应用地震切片的应用 三维可视化技术三维可视化技术 数据体可视化基础数据体可视化基础 VoxelGeo基本工作流程 PG2PG2输入输入SeisX输入输入ULA输入输入数据体浏览,发现感兴趣数据体浏览,发现感兴趣的目标的目标调整透明度等可视化参数调整透明度等可视化参数小小倾倾角角情情况况Trim volume自动追踪或手动解自动追踪或手动解释层位释层位子体雕刻,三维地子体雕刻,三维地质解释质解释三维井位设计三维井位设计三维可视化技术三维可视化技术 主要特点主要特点 •从后向前快速的数据浏览,使得解释员可以快速侦察数据体内部结构并建立起复杂构造和地层关系的整体概念•·基于波形相关对比的快速、准确自动层位追踪解释模块•·快速的三维分析,快速发现确定地质异常,如河道、扇体、亮点等•·直接的振幅解释、子体雕刻•·断层面、层面和井轨迹与地震数据镶嵌显示•·解释管理•·三维导航•·动画功能•·独特的可变透明度的高质量数据浏览显示与子体雕刻、拉平的结合,使得对三维数据体内部的可视化成为可能,并用于识别高孔隙区带、沉积环境、流体成分以及断层。

      •·与StratmagicTM的结合•·开放的结构•允许用户以VoxelGeo为基础开发自己的应用模块从而扩展其功能VoxelGeo可以与单个应用模块或同时运行•·基于OpenInventor和OpenGL•·支持多个数据体浏览(最多可达5个)三维可视化技术三维可视化技术 软件应用中主要参数的意义软件应用中主要参数的意义及参数选择原则及参数选择原则 三维可视化技术三维可视化技术 透明度(opacity) 软件应用中主要参数的意义软件应用中主要参数的意义及参数选择原则及参数选择原则 三维可视化技术三维可视化技术 调色板调色板 软件应用中主要参数的意义软件应用中主要参数的意义及参数选择原则及参数选择原则 三维可视化技术三维可视化技术 颜色直方图颜色直方图((Histogram)) 子体雕刻参数 软件应用中主要参数的意义软件应用中主要参数的意义及参数选择原则及参数选择原则 三维可视化技术三维可视化技术 软件应用中主要参数的意义软件应用中主要参数的意义及参数选择原则及参数选择原则 三维可视化技术三维可视化技术 剖面上种子点图 由种子点自动追踪结果(第一次) 软件应用中主要参数的意义软件应用中主要参数的意义及参数选择原则及参数选择原则 三维可视化技术三维可视化技术 第一次自动追踪结果(底图) 最终自动追踪结果(剖面) 最终自动追踪结果(底图) 复杂构造的快速解释复杂构造的快速解释 三维可视化技术三维可视化技术 复杂构造的快速解释复杂构造的快速解释 三维可视化技术三维可视化技术 自动追踪 复杂构造的快速解释复杂构造的快速解释 三维可视化技术三维可视化技术 原始地震数据可视化得到清晰的断裂系统 最终的构造图 河道、扇体解释河道、扇体解释 三维可视化技术三维可视化技术 复杂的河道系统复杂的河道系统 河道、扇体解释河道、扇体解释 三维可视化技术三维可视化技术 断层断层扇三角洲扇三角洲河道系统河道系统扇三角洲和河道系统扇三角洲和河道系统 •全三维解释离不开三维可视化技术。

      •对层位和断层面的自动追踪技术是全三维解释的基础•主要方法包括自动拾取、层面切片和体元(Voxel)追踪 全三维解释全三维解释 振幅体、相干体可视化浏览了振幅体、相干体可视化浏览了解构造背景,确定解释方法解构造背景,确定解释方法断裂体系解释断裂体系解释在一个或多个数据体在一个或多个数据体中选择一个或多个种中选择一个或多个种子点进行体素单元自子点进行体素单元自动追踪圈定异常体动追踪圈定异常体全三维可视化解释流程全三维可视化解释流程全三维构造解释全三维构造解释异常体顶底异常体顶底面层位提取面层位提取层位解释层位解释可视化技术可视化技术指导振幅体指导振幅体和相干体的和相干体的断层解释断层解释振幅体、相干体、属性体、储层特振幅体、相干体、属性体、储层特征反演数据体可视化浏览了解沉积征反演数据体可视化浏览了解沉积背景和发现异常体背景和发现异常体储层沉积相解释储层沉积相解释抽稀网格层位抽稀网格层位解释解释选择种子点选择种子点体素单元自体素单元自动追踪和层动追踪和层位提取位提取区域自区域自动追踪动追踪面块切面块切片解释片解释锁定时窗或锁定时窗或锁定层位可锁定层位可视化沉积相视化沉积相解释解释可视化技术检验修改解释成果可视化技术检验修改解释成果三维断层模型三维断层模型综合评价及井位设计综合评价及井位设计 为什么要进行全三维解释?为什么要进行全三维解释?传统的三维资料二维解释方法:传统的三维资料二维解释方法:1.解释效率低解释效率低2.解释精度低解释精度低3.二维平面认识三维空间二维平面认识三维空间,认识局限认识局限 振幅异常振幅异常振幅异常振幅异常全三维解释方法全三维解释方法体元自动追踪技术体元自动追踪技术GeoQuestGeoQuest公司北海实例公司北海实例 可视化技术在全三维解释中的应用可视化技术在全三维解释中的应用i数据体浏览数据体浏览i断裂体系解释断裂体系解释i全三维层位解释全三维层位解释i储层解释储层解释i烃类检测烃类检测i综合评价与井位设计综合评价与井位设计 数据体浏览数据体浏览利用可视化技术进行数据体浏览是全三维解释中最重要的利用可视化技术进行数据体浏览是全三维解释中最重要的工作之一,起作用表现在:工作之一,起作用表现在:i在全三维构造解释中,振幅体和相干体的浏览在全三维构造解释中,振幅体和相干体的浏览可以了解构造形态,断裂空间展布特征,检查数可以了解构造形态,断裂空间展布特征,检查数据体质量,从而据体质量,从而有针对性地选择解释方法和控制有针对性地选择解释方法和控制参数参数i在全三维储层沉积相解释中,利用体可视化显在全三维储层沉积相解释中,利用体可视化显示技术对多属性数据体浏览可以确定储层敏感参示技术对多属性数据体浏览可以确定储层敏感参数,了解沉积背景,发现异常体数,了解沉积背景,发现异常体 数据体浏览数据体浏览方差体方差体和振幅和振幅体综合体综合可视化可视化显示和显示和数据体数据体扫描了扫描了解断裂解断裂系统的系统的空间分空间分布特征布特征 InLine地震剖面地震剖面ABCDABCD950ms处振幅切片处振幅切片ABCD不整合面向下不整合面向下4ms-20ms 振幅体三维可视化图振幅体三维可视化图900-1000ms方差方差体体 三维可视化图三维可视化图断层断层ABCD断距较小的断层在地震断距较小的断层在地震剖面和振幅切片中模糊剖面和振幅切片中模糊不清,难以解释,但是不清,难以解释,但是利用方差数据体并进行利用方差数据体并进行三维可视化显示可以准三维可视化显示可以准确解释断层确解释断层断裂体系解释断裂体系解释 FnFn2052ms-2108ms方差体方差体沙三下油层顶面构造图Fn断裂体系解释断裂体系解释FnFn断断层层地地震震反反射射特特征征和和平平面面分分布布 区域自动追踪技术区域自动追踪技术面块切片解释技术面块切片解释技术体素单元追踪自动技术体素单元追踪自动技术全全三三维维层层位位解解释释方方法法区域自动追踪技术区域自动追踪技术全三维层位解释全三维层位解释 全三维层位解释方法全三维层位解释方法区域自动追踪技术区域自动追踪技术 控控制制条条件件振幅变化范围振幅变化范围相关分析相关分析稀疏解释数据稀疏解释数据倾角变化范围倾角变化范围断层边界限制断层边界限制8080年代初出现,方法日益完善年代初出现,方法日益完善稀疏解释数据稀疏解释数据断层边界限制断层边界限制断层边界限制断层边界限制相关分析相关分析断层边界限制断层边界限制 全三维层位解释方法全三维层位解释方法区域自动追踪技术区域自动追踪技术松辽盆地南部,面积松辽盆地南部,面积144km144km2 2, ,地震资料地震资料信躁比高,解释层位构造简单。

      信躁比高,解释层位构造简单 全三维层位解释方法全三维层位解释方法面块切片解释技术面块切片解释技术面块切片(面块切片(surface slice)是采用是采用全新方式显示全新方式显示的一个的一个薄的地震数据体薄的地震数据体,可以理解为由若干水平切片叠合,可以理解为由若干水平切片叠合而成的数据体其显示方式类似于体可视化显示,而成的数据体其显示方式类似于体可视化显示,所不同的是面块切片所不同的是面块切片突出的是目的层的极性突出的是目的层的极性,凡是,凡是不想解释的极性都不显示不想解释的极性都不显示面块切片的极性有六种选择:波峰(面块切片的极性有六种选择:波峰(peak)、)、波波谷(谷(trough)、、负的波峰(负的波峰(local peak)、、正的波谷正的波谷((local trough)、、从正到负零值点和从负到正零值从正到负零值点和从负到正零值点解释过程中根据目的层波形选择显示极性点解释过程中根据目的层波形选择显示极性 全三维层位解释方法全三维层位解释方法面块切片解释技术面块切片解释技术5 5张切片张切片10ms10ms面块面块切片显示,切片显示,黑色为背黑色为背景色,绿景色,绿色为色为100*100*100100测网解释测网解释数据投影数据投影其它颜色其它颜色为波峰显为波峰显示示 全三维层位解释方法全三维层位解释方法体元自动追踪技术体元自动追踪技术体元自动追踪技术是随着三维可视化的出现而产生的,体元自动追踪技术是随着三维可视化的出现而产生的,其自动追踪过程是从解释人员定义种子体元(其自动追踪过程是从解释人员定义种子体元(seed voxel)开始的,和区域自动追踪技术不同的是,体元追开始的,和区域自动追踪技术不同的是,体元追踪是沿着真正的三维路径追踪数据体的,追踪出来的踪是沿着真正的三维路径追踪数据体的,追踪出来的是数据体,而不是解释层位是数据体,而不是解释层位体元自动追踪的控制参数有两个:属性值变化范围和体元自动追踪的控制参数有两个:属性值变化范围和体元追踪方式。

      如选择体元追踪方式如选择6个方向追踪方式限定了一个体个方向追踪方式限定了一个体元只对最邻近的(面接触)体元进行追踪,而元只对最邻近的(面接触)体元进行追踪,而26个方个方向追踪方式则允许体元追踪沿着两个邻近的面到面、向追踪方式则允许体元追踪沿着两个邻近的面到面、边到边和角到角的体元方向追踪边到边和角到角的体元方向追踪 全三维层位解释方法全三维层位解释方法三种全三维解释技术比较三种全三维解释技术比较假定条件:假定条件:三种解释技术都假定数据体是局部连三种解释技术都假定数据体是局部连续的,相位一致续的,相位一致体素单元自动追踪体素单元自动追踪可以进行多属性体约束追踪,追踪出来的可以进行多属性体约束追踪,追踪出来的是满足约束条件的数据体,适用于异常体自动追踪,是满足约束条件的数据体,适用于异常体自动追踪, 在层位在层位解释方面其优势是可以在相干体(或方差体)和振幅体的共解释方面其优势是可以在相干体(或方差体)和振幅体的共同约束下进行追踪同约束下进行追踪面块切片技术面块切片技术在地层倾角小的情况下进行构造和振幅异常体在地层倾角小的情况下进行构造和振幅异常体解释具有明显优势,解释速度快、质量高,编辑工作量小解释具有明显优势,解释速度快、质量高,编辑工作量小区域自动追踪技术区域自动追踪技术适用于在不同资料品质和不同地层倾角条适用于在不同资料品质和不同地层倾角条件下进行构造解释,约束参数的选择和稀疏测网的解释精度件下进行构造解释,约束参数的选择和稀疏测网的解释精度直接影响自动追踪结果直接影响自动追踪结果 实例分析:实例分析:•地震数据:墨西哥湾的一块小三维数据体地震数据:墨西哥湾的一块小三维数据体•测井数据:测井数据:2828口井口井•构造位置:盐丘构造位置:盐丘测试内容:测试内容:•1.1.对原始数据进行检查并建立断层框架模型对原始数据进行检查并建立断层框架模型•2.2.对断层解释质量进行检查对断层解释质量进行检查•3.3.绘制盐丘顶面构造图绘制盐丘顶面构造图•4.4.确定远景区确定远景区全三维构造解释与常规解释比较全三维构造解释与常规解释比较资料来自资料来自GeoQuest公司公司 解释成果:解释成果:•传统方法共解释出传统方法共解释出7条断层、一层盐丘顶面构造图、圈定一个远景区。

      条断层、一层盐丘顶面构造图、圈定一个远景区•基于基于GeoViz的全三维解释方法共解释出的全三维解释方法共解释出14条断层、一层盐丘顶面构造图、圈条断层、一层盐丘顶面构造图、圈定一个远景区定一个远景区对比结果对比结果全三维构造解释与常规解释比较全三维构造解释与常规解释比较 基于基于GeoVizGeoViz的的解释成果解释成果基于传统方法的基于传统方法的解释成果解释成果全三维构造解释与常规解释比较全三维构造解释与常规解释比较 引自引自李玲等李玲等:全三维解释方法探讨与实践全三维解释方法探讨与实践,石油地球物理勘探石油地球物理勘探1996年年8月月资料背景资料背景:面积面积700km2,地层倾角小地层倾角小,断层少断层少,地震资料品质好地震资料品质好全三维构造解释与常规解释比较全三维构造解释与常规解释比较 引自引自王咸彬等王咸彬等:真(全)三维构造解释技术初探真(全)三维构造解释技术初探,石油物勘探石油物勘探2000年年6月月全三维构造解释与常规解释比较全三维构造解释与常规解释比较 三、基于多地震属性模式识别三、基于多地震属性模式识别的储层预测的储层预测一)、统计模式识别方法一)、统计模式识别方法 二)、神经网络分析二)、神经网络分析 三)、分带压缩和主成分分析三)、分带压缩和主成分分析 四)、其它模式识别方法四)、其它模式识别方法 •油气统计模式识别一种根据含油气与不含油气储层的地震波运动学和动力学特征(如波形、振幅、频率、相位等)的差异,从地震资料中提取多种地震属性,采用多元统计的方法,预测含油气储层的位置与范围的一种技术。

      实现过程可分如下两步: –学习过程 –预测过程 统计模式识别方法统计模式识别方法 自相关函数的第二极大值比极大值付氏谱最大功率谱值与总能量之比修正协方差谱频率带宽自回归2阶AR系数伯格谱75%累计功率谱能量对应的频率伯格谱高频区能量与总能量之比 分形振幅谱维数计盒维数突变次数Lyapunov指 数分 形RS维灰色系统发展系数 KL变换的第1主分量 KL变换的第2、3、4、5主分量 •80年代以来产生了许多利用地震信息预测油气的方法,但由于储层的非均质性和各区块之间的差异越来越大,因此需要有一种适应非线性问题的好的方法来预测油气,而神经网络法多参数预测是一种较理想的方法施继承(1996)年提取了32种参数进行预测,获得了成功 神经网络分析神经网络分析 •BP即反向传播,由于BP算法对多层感知器的权值学习和网络学习比较有效,因而在许多文献中,BP与多层感知器是紧密相连的;•BP网络的应用分为学习和识别两个部分; 神经网络分析-神经网络分析-BPBP神经网络神经网络 •SOM网络是自联想最邻近分类器,能够将任意个连续值模式分成P类网络采用竞争学习规则线学习,既可采用有导师学习,也可进行无导师学习,因而SOM网络可适用于已知样本较少甚至没有的情形下的分类、 聚类问题。

      •SOM网络进行分类时的流程是:①确定输出层维数及大小;②网络进行自学习;③根据自学习的结果,并依据已知样本在输出层的响应点位置对网络的输出层进行标定,形成分类模板;④利用分类模板对未知样本逐一进行分类 神经网络分析-神经网络分析-SOMSOM神经网络神经网络 基于神经网络的烃类检测系统基于神经网络的烃类检测系统 •基于神经网络的烃类检测系统主要包括四大部分:–①储层识别和追踪;–②沿储层的多参数提取;–③特征压缩;–④烃类检测 神经网络分析神经网络分析 基于专家经验地震属性有效性分析方法基于专家经验地震属性有效性分析方法 •地震属性优选应遵循的基本准则地震属性优选应遵循的基本准则–①优选后属性集整体与研究对象具有某种相关性,能够对样本进行有效分类; ②)达到属性结构的最优化,以尽可能相互独立的变量组成尽可能低维的变量空间; ③使有用信息损失为最小,剔除起干扰作用的属性 •地震属性有效性分析的工作程序地震属性有效性分析的工作程序 –①地震属性进行初选;②根据实际工区地质特征建立的地震地质模型对初选地震属性进一步分析和筛选;③井形成初选地震属性集;④达到有效分类;⑤如不能达到有效分类,重复以上各步。

      •地震属性优选主要算法地震属性优选主要算法 –①顺序前进法(SFS);顺序后退法(SBS) ;增l减r法(l-r法) 神经网络分析神经网络分析 A平均反射强度平均反射强度平均瞬时频率平均瞬时频率平均瞬时相位平均瞬时相位气层顶构造图气层顶构造图 (D1) (D1) 第二优势频率第二优势频率(F2)(F2)BAABAA AB地地震震属属性性分分析析预测结果预测结果AB 神神经经网网络络模模式式识识别别 苏里格庙地区盒苏里格庙地区盒8 8段储层烃类检测图段储层烃类检测图神神经经网网络络模模式式识识别别 PSS6准层序组岩性展布的BP 神经网络预测 PSS6准层序组岩性展布的SOM神经网络预测 分带压缩和主成分分析分带压缩和主成分分析三维地震相体三维地震相体 #1振幅振幅, 波阻抗波阻抗, 流体因子流体因子分带压缩分带压缩Zonation谱系谱系 地震相地震相 分类分类三维地震相体三维地震相体 #2振幅包络振幅包络, 波阻抗波阻抗, 流体因子流体因子,相干体相干体主成分分析主成分分析PCA谱系谱系 地震相地震相 分类分类AB三维地震数据体s三维地震数据体图 6-4: A - M1层Time图 (波峰)。

      黑色虚线表示inline 332的位置 B –西-东 inline 332 M1 层以红色显示储层段位于两条蓝线之间 注意拉帕尔玛油田位于主南-北背斜构造的高点 分带压缩和主成分分析分带压缩和主成分分析BLPT-3LPT-4LPT-3LPT-4ALPT-3LPT-4LPT-1XLPT-2CRS-1XM15 kmNWE图 6-5: 三维多属性地震分类工作流程对大套砂岩段进行分类的两种方法A –分带压缩分类; B –主成分分析分类注意两种方法使用了不同的输入地震属性数据体 分带压缩和主成分分析分带压缩和主成分分析LPT-4SPLLDSand baseSand topAWELPT-3SPLLDSand baseSand topWEB图 6-6: 穿过地震相体(由分层压缩方法计算)的西东向主测线,A - Inline 329 经过 LPT-4 生产井附近 B - Inline 333 穿过 LPT-3 含水井位置注意从含水井向东(上倾方向)存在黄色地震相段和变厚的红色地震相层段(图B白色箭头)它表明在这一构造高部位具有良好的储层特性 分带压缩和主成分分析分带压缩和主成分分析2 kmNLPT-3LPT-4LPT-2LPT-1X图 6-7: 地震相分类数据体(用分带压缩法计算)的一个层位切片。

      此切片位于M1层以下32ms注意由黑线圈出的北东-西南向黄色地震相带它被解释为含油砂坝注意由白线圈出的北西-南东向红色地震相带,被解释为河道 分带压缩和主成分分析分带压缩和主成分分析AB图 6-8 A – 用于多属性数据体分类的地层模型 B – 由分带压缩法计算出的地震相体的三维显示数据体显示在一个24ms的时窗内(从M1层以下24ms到48ms)黄色地震相显示为全透明,而其它地震相显示为半透明方式注意北东-南西向黄色-红色地震相带(由黑色虚线圈出的辉亮部分),该相带上完钻三口生产井含水井位于该相带西面的蓝色地震相区 分带压缩和主成分分析分带压缩和主成分分析2 kmNLPT-3LPT-4LPT-2LPT-1X图 6-9 PCA法计算的地震相数据体切片该切片位于M1层以下28ms在井标定的基础上,黄色地震相解释为直接烃类指示在流体预测,地层和构造解释的基础上,获得了更多的储层知识这增加了识别有利区的可信度,起初有利区识别仅以构造和振幅分析为基础见图 10地震相数据体和常规振幅数据体层位切片对比图)作为成果,预测出几个有利区并建议实施钻探 分带压缩和主成分分析分带压缩和主成分分析AB5 kmN图 6-10 对比地震相数据体切片 (A)和常规振幅图切片 (B).两幅图都和M1层等值线图叠合。

      注意多数生产井位于振幅异常体的外面同时注意生产井钻遇的北东-南西相带在常规振幅图上看不到 本项研究对两种不同的多属性数据体分类方法进行本项研究对两种不同的多属性数据体分类方法进行了对比分析;其一包括地震道方波化(分带压缩),而了对比分析;其一包括地震道方波化(分带压缩),而另一种使用另一种使用PCA技术两种分类方法达到的分辨率是不技术两种分类方法达到的分辨率是不同的使用分带压缩法计算出的地震相体垂向分辨率较同的使用分带压缩法计算出的地震相体垂向分辨率较低,而且为了预测岩性和流体内容需要对解释层段内整低,而且为了预测岩性和流体内容需要对解释层段内整个复杂的地震相进行详尽的分析然而,使用交汇图使个复杂的地震相进行详尽的分析然而,使用交汇图使解释更为可信另一方面,包括主成分分析的分类技术解释更为可信另一方面,包括主成分分析的分类技术改善了垂向分辨率而且在井标定基础上可以直接进行烃改善了垂向分辨率而且在井标定基础上可以直接进行烃类预测虽然可以做这项工作,但是缺乏将地震相分类类预测虽然可以做这项工作,但是缺乏将地震相分类结果和输入地震属性进行相关的功能,而且不能进行进结果和输入地震属性进行相关的功能,而且不能进行进一步的定量分析。

      综合两种方法可以使结果更可靠一步的定量分析综合两种方法可以使结果更可靠 分带压缩和主成分分析分带压缩和主成分分析 •模糊模式识模糊模式识 •分形油气预测方法分形油气预测方法 •灰色油气预测方法灰色油气预测方法 •RS理论决策分析理论决策分析 其它模式识别方法其它模式识别方法 四、基于统计学随机建模和正四、基于统计学随机建模和正演模拟的储层预测演模拟的储层预测•多参数地震特征反演多参数地震特征反演 •地质统计学随机建模和储层物性反演地质统计学随机建模和储层物性反演 多参数地震特征反演多参数地震特征反演 技术基础技术基础 •((1)地质模型建立;)地质模型建立;•((2)主组分分析;)主组分分析;•((3)模型估算技术)模型估算技术多参数地震特征反演多参数地震特征反演 创建标准正交基产生一组任意的矢量 形成一假想的合成曲线 对初始参数值的软约束作用 平滑约束的作用 多参数地震特征反演多参数地震特征反演 •地质统计学是处理空间分布变量的应用统计学的一个分支近年来人们综合岩芯、测井、地震和地质等各种软信息和硬资料,不断探索运用地质统计学方法进行储层参数预测,以建立精细的储层地质模型,并且随着计算机技术的发展,地质统计学方法不断地模块化,它已经成为地质学家、地球物理学家和油藏工程师进行空间数据分析的重要工具和共同语言,在储层的参数预测方面起着极为重要的作用。

      地质统计学随机建模和储层物性反演地质统计学随机建模和储层物性反演 地质统计学随机建模和储层物性反演地质统计学随机建模和储层物性反演 变差函数示意图变差函数示意图 变差函数变差函数 ((h))与协方差函数与协方差函数C((h))的关系的关系 地质统计学随机建模和储层物性反演地质统计学随机建模和储层物性反演 克里格法克里格法 地质统计学随机建模和储层物性反演地质统计学随机建模和储层物性反演 克克立立格格法法的的实实质质是是根根据据待待估估点点周周围围的的若若干干已已知知信信息息,,应应用用变变差差函函数数的的性性质质,,对对待待估估点点的的值值作作出出最最优优((即即估估计计方方差差最最小小))、、无无偏偏((即即估估计计值值的的均均值值与与观观测测值值的的均均值值相相同同))的的估估计计更更具具体体地地说说,,克克里里格格法法是是在在考考虑虑了了信信息息样样品品的的形形状状、、大大小小及及其其与与待待估估块块段段相相互互之之间间的的空空间间分分布布位位置置等等几几何何特特征征,,以以及及变变量量的的空空间间结结构构信信息息后后,,为为了了达达到到线线性性、、无无偏偏和和最最小小估估计计方方差差,,而而对对每每个个样样品品值值分分别别赋赋予予一一定定的的权权系系数数,,最最后后用用加加权权平平均均法法对对待待估估计计块块段的未知量进行估计的方法。

      段的未知量进行估计的方法 普通克里格预测泥质含量分布图普通克里格预测泥质含量分布图 地质统计学随机建模和储层物性反演地质统计学随机建模和储层物性反演 协克里格法预测泥质含量分布协克里格法预测泥质含量分布 地质统计学随机建模和储层物性反演地质统计学随机建模和储层物性反演 油藏建模技术油藏建模技术构造模型构造模型岩石物理模型岩石物理模型孔隙度模型孔隙度模型渗透率模型渗透率模型油藏剖面模型油藏剖面模型油藏三维模型油藏三维模型 。

      点击阅读更多内容
      关于金锄头网 - 版权申诉 - 免责声明 - 诚邀英才 - 联系我们
      手机版 | 川公网安备 51140202000112号 | 经营许可证(蜀ICP备13022795号)
      ©2008-2016 by Sichuan Goldhoe Inc. All Rights Reserved.