振幅、频率与动态吸收系数在嘉一段储层预测中的应用
Ｄｅ ．２ １ Ｃ ００
振 幅 、 率 与动 态 吸收 系数 频 在 嘉 一 段储 层 预 测 中的应 用
曹 健， 赵宪生 ， 泽进 ， 施 于静芳
（ 成都理工大学 ， 四川 成都 ６ ０５ ） １０ ９
摘要 ： 绍了频率 、 幅和动 态吸收 系数的计算方法原理 。针对 Ｊ 东南赤水地 区嘉陵江组一段 介 振 Ｉ １
１ 地震属性的分析原理和 计算
地震属性 是指 由叠前 或者叠 后地震 数据 ， 经数 学 变换而得 出 的有 关地震 波 的几 何形 态 、 动学特 运
（ 地震 波长度 ）分析各 层的速度分布 ， ， 上部为嘉二段 中上部 的高速 白云岩 ， 下部为嘉一 段的高速 灰岩 ， 因
此在地震剖 面上的波形特征应该是该 套储层 的总体 振 幅响应 。根据各 层 的速度 ， 建立 如 图 １ 示 的地 所
Ｉ ｓ由图可知 ， ｎ ， ／ 随着厚度的减小 在速度不变条件 下 ． 振幅由强减弱 ， 在厚度小于 ５１后 振 幅增加 ， Ｉ Ｔ 但极性
， ，
玎以看 出 ， Ｊ 该地 区 目的层 振 幅值主 要分布 于西北 和 东南部 ， ４中黄 红 色分别 表示 了 中一强 振 幅 带 图
，
根 据井 资料分 析表 明 ， 区段 为孔 隙发育 较好 的颗 该
第 １ ７卷第 ６ 期
２１ ０ ０年 １ ２月
文 章 编 号 ：１０ ６ ３ （０ ０ ０ ０ ３ ０ ０ ６— ５ ５ ２ １ ）６— ０ ５— ３
特 种 油 气 藏
Ｓ ｅ ｉｌＯ ｌ ｎ ｓＲｅ ｅ ｖ ｉ ｐ ｃａ ｉ ａ ｄ Ｇａ ｓ ｒ ｏｒ ｓ
Ｖ０．１ ． １ ７ Ｎｏ ６
态 吸收系数方 法 即提 取 动态 地 震子 波 的 振 幅谱 和 相位谱 ， 由相邻 层 问子 波谱对 数差 的斜率计 算 吸收
系数 。以上属性 均 与岩层岩 性有关 ， 结合储 层特 征
分析 ， 推断储层 岩性 分布 。 可
收 稿 日期 ：０ ０５ ９； 回 日期 ：０ ０ ８９ ２ １０ １ 改 ２ １０ ０
究砂岩 岩性油 气 藏 ，确定 储 层 的分 布 已经取 得 了
显 著 的效 果 ¨ 。该 方 法 认 为 地震 波属 性 是 地 下
裂缝 孔 隙型储层 ， 储层厚 度小 ， 具有低 孔低 渗 、 非均
质性 强 的特 点 。通过 对 嘉 陵 江组 一段 的钻 井 资料
ห้องสมุดไป่ตู้
地 质结构 、岩性 和所 含 流 体性 质 等 多 种 因素 的 综
区西 北 部 。
关 键 词 ： 幅 ； 率 ； 态 吸 收 系数 ； 振 频 动 地震 波 ； 层 预 测 ； 一 段 储 嘉 中 图分 类号 ：Ｅ ２． Ｔ ｌ２ ２ 文 献标 识 码 ： Ａ
引 言
自２ ０世 纪 ９ ０年代 以来 ，利 用 地 震 波 属性 研
２ 储层特征
川 东南 赤 水地 区下三 叠 统嘉 陵 江组 一段 属 于
基金项 目： 国家 自然科学基金委员会与中国石油化工股份有限公司联合基金 资助项 目“ 川东南地区重点层位储层沉积相研究” ４ ７ ９０ ） （ ０ ３９ ３ 作者简介 ： 曹健（ ９ ５ ， ２０ 年毕业于桂林工学 院勘查技术与工程专业， １８ 一） 男，０ ８ 现为成都理工大学在读硕士研究生 ， 从事储层预测研究 。
粒灰岩 ， 弱振 幅表示致 密岩层 发 育
速度 Ｖ ＝ １ ／ ； 间楔形 层 ， 度从 ｈ＝ ０ｍ １ ６ １０ｍ ｓ 中 厚 ３ 到 ｈ＝ 变化 ，速 度从 左 至 右 分 别 为 ｖ ＝６５０ ０ｍ ２ ９ ｍｓ ／ 到 ＝ １ ｍ／ ；下 部 速 度 ＝６ １０ ｒ ｓ ５９ １ ｓ ４ ／ ， ｎ 由剖面可知 ， 随着厚 度 的变薄 ， 速度 的减 小 ， 幅 由 振
３ ６
特 种 油 气 藏
第 １ ７卷
。
强变弱 。图 ３为变厚度楔形模 型合成 记录 剖面
Ｅ
部速度 ＝ １ ／ ； ６１０ｍ ｓ 中间楔形层厚度从 ： ０ｍ 到 ３
图 ４是嘉 陵 江一段顶 部最 大波峰振 幅平 面 图
＾ Ｏｍ变化 ， 速度 ＝ ７ ／ ； ６１ ｓ下部速度 ： ０ ７ｍ ６１ ４
结 构 、岩 性和所 含 流 体进 行 分 析 和研 究 Ｊ 以川 。
东 南赤水地 区嘉 陵江组 一段 为例 ，根据 已有地质 、
钻 井和测 井 资料 ，结 合储 层 和 岩性 特 征 ，利 用振
幅、频率与动态吸收系数分析预测了该 区块的储
层 分布 。
３ 储层地震响应特征 分析
由于嘉 一段 顶 部储 层较 薄 ， 层 厚 度小 于 １４ 单 ／
征 、 力学 特 征 和统 计 学 特 征 的 特殊 量 值 Ｊ 动 。利
用实 际地震道记 录 ， 过希 尔伯特 变换 可以得 到地 通 震记 录 的正 交道 ， 而计 算 出振 幅 、 率属 性 。动 从 频
质模型， 并用过井剖面提取的子波研究剖面的波形
特征 ， 应用 该 子波 合成 楔 形地 质 模 型 正演 记 录 （ 图 ２３ 。图 ２为变速度楔形模型合成记 录剖 面。上部 、）
储层特性 ， 结合钻 井和测井资料研 究， 建立 了研 究 区嘉一段有利 油气富集带的地震振 幅响应模
式 ， 究认 为 中 高振 幅 为 其 地 震振 幅 响 应模 式 。应 用该 模 式 对提 取 的振 幅 、 率 和 动 态吸 收 系 研 频
数的平面图标 定分析研究 ， 测了研 究区嘉一段 有利储层分布 ， 预 储层分布有利区主要位 于研 究
合 反 映 ，综合分 析 地 震属 性 参数 可 以对地 下 地 质
分析 ， 岩性 对储层 的控 制作 用 比较 明显 ， 岩 心及 在 薄 片观察发 现颗 粒 灰 岩溶 蚀 孔 隙发 育 ， 集性 好 。 储
因此 ， 应将川 东南赤 水地 区嘉 陵江组一段 有利储 层 岩性鲕 粒灰岩 作为重 点勘 探 目标 。Ｊ 。
储层预测
5多元统计方法在储层预测中的应用 5.1多元统计方法原理 5.2多元统计方法的应用 6模式识别技术及其在储层预测中的应用 6.1统计模式识别技术的基本原理与应用 6.2人工神经网络基本原理与应用 7边缘检测技术与应用 7.1边缘检测技术的基本原理 7.2基于小波边缘检测技术与应用 7.3基于分形边缘检测技术与应用
8地震波阻抗反演方法理论与应用 8.1波阻抗反演的基本原理 8.2测井约束地震反演方法的应用 9多波多分量地震储层描述与应用 9.1多波多分量地震勘探的理论基础 9.2多波多分量地震资料采集与处理 9.3多波地震资料解释和储层描述 10 四维地震勘探技术与应用 10.1四维地震的可行性与研究前提 10.2四维地震资料处理方法与原则 10.3四维地震资料解释方法
岩石类型
速度 （米/秒）
砾岩碎石干砂
200~800
砂质粘土
300 ~ 900
湿砂
600 ~ 800
粘土
1200 ~ 2500
砂岩
1400 ~ 4500
泥灰岩
2000 ~ 3500
石灰岩，白云岩
2500 ~ 6100
泥质页岩
2700 ~ 4100
盐岩
4200 ~ 5500
几种沉积岩的波速
三、地震波速度与岩石密度的关系
k i1
i 1 N
( fi f )2
i 1
K:反映吸收系数的大小
求出自相关函数
N
i 2k
X (n) x(k)e N
i 1
N
S X 2 (n) n1
总能量
X max MAX [x(1), x(2), x(3),..., x(n)]
储层预测技术详解
LPM 储层预测技术LPM 是斯伦贝谢公司GeoFrame 地震解释系统中最新推出的储层预测软件，利用地震属性体来指导储层参数（如砂岩厚度）在平面的展布，以此来实现储层参数的准确预测。
LPM 预测储层砂体可分两步进行：首先，它是将提取的地震属性特征参数与井孔处的砂岩厚度、有效厚度进行数据分析，将对储层预测起关键作用的地震属性特征参数优选出来，根据线性相关程度的大小，建立线性或非线性方程。
线性方程的建立主要采用多元线性回归方法；非线性方程的建立主要采用神经网络方法；其次，根据建立的方程，利用网格化的地震属性体来指导储层参数（如砂岩厚度）在平面的成图。
4.1.1多元线性回归基本原理设因变量y 与自变量x 1， x 2 ，…，x m 有线性关系，那么建立y 的ｍ元线性回归模型：ξβββ++++=m m x x y 110其中β0，β1，…，βm 为回归系数；ξ是遵从正态分布N(0，σ2)的随机误差。
在实际问题中，对y 与x 1， x 2 ，…，x m 作n 次观测，即x 1t ， x 2t ，…，x mt ，即有：t mt m t t x x y ξβββ++++= 110建立多元回归方程的基本方法是：(1)由观测值确定回归系数β0，β1，…，βm 的估计b 0，b 1， …，b m 得到y t 对x 1t ，x 2t ，…，x mt ；的线性回归方程：t mt m t t e x x y ++++=βββ 110其中t y 表示t y 的估计；t e 是误差估计或称为残差。
(2)对回归效果进行统计检验。
(3)利用回归方程进行预报。
回归系数的最小二乘法估计根据最小二乘法，要选择这样的回归系数b 0，b 1， …，b m 使∑∑∑===----=-==nt n t mt m t t t t n t tx b x b b y y y e Q 11211012)()( 达到极小。
为此，将Q 分别对b 0，b 1， …，b m 求偏导数，并令0=∂∂bQ ，经化简整理可以得到b 0，b 1， …，b m ，必须满足下列正规方程组：⎪⎪⎩⎪⎪⎨⎧=+++=+++=+++my m mm m m y m m y m m S b S b S b S S b S b S b S S b S b S b S22112222212111212111 m m x b x b x b y b ----= 22110其中∑==nt t y n y 11 m i x n x nt it i ,,2,111==∑= ),,2,1())((1))((1111m i x x n x x x x x x S S nt n t jt it jt n t it j jt i n t it ji ij =-=--==∑∑∑∑==== ),,2,1())((1))((1111m i y x n y x y y x x S nt n t t it n t t it t i n t it iy =-=--=∑∑∑∑====解线性方程组，即可求得回归系数i b ，将i b 代入式可求出常数项0b 。
地震属性含义及其应用
地震属性含义及其应用一、瞬时属性19 假定复数道表示为：u(t) = x(t) • iy(t)，则1. 瞬时实振幅IReAmp ( In sta nta neous Amplitude )是在选定的采样点上地震道时域振动振幅。
是振幅属性的基本参数。
广泛用于构造和地层学解释。
用来圈定高或低振幅异常，即亮点、暗点。
反映不同储集层、含气、油、水情况及厚度预测。
2. 瞬时虚振幅IQuadAmp (I nst. Quadrature Amplitude)是复数地震道的虚部，与复数地震道的相位为90o时的时域振动振幅。
即正交道,为虚振幅。
因它只能在特定的相位观测到，多用来识别与薄储层中的AVO异常。
3. 瞬时相位IPhase ( Instantaneous Phase)(t)二Atan(y(t).x(t)),定义为正切，输出相位已转换为角度，数值范围是[-180o,180°]。
为q(t)/f(t)的一个角，是采样点处地震道的相位。
有助于加强储层内部的弱反射同相轴，但同时也加强了噪声，可用于指示横向连续性；显示与波传播有关的相位部分；用于计算相速度；因为没有振幅信息因此能够显示所有同相轴；用于显示不连续；断层、显示层序边界。
由于烃类聚集常引起局部相位变化，也可以做烃类直接指示之一。
4. 瞬时相位余弦CIP ( Cos ine of In st. Phase )是瞬时相位导出的属性。
其计算式为Cos( (t))常用来改进瞬时相位的变异显示。
并用于相位追踪和检查地震剖面对比、解释的质量。
多与瞬时相位联用。
5. 瞬时频率IFreq (I nst. Freque ney)定义为瞬时相位对时间的函数 d (t) dt (以度/毫秒或弧度/毫秒表示)，其量纲为频率的量纲(Hz)，是地震道在频率方面的瞬时属性。
用来计算、估算地震波的衰减。
油气储层常引起高频成分衰减及杂乱反射显示，所以横向上可用于碳氢指示。
高频成份多显示为尖锐的界面或薄层，亦可反映岩相的粗、细变化及地层旋回。
松辽盆地两江地区扶余油层河道储层描述方法与应用
在沉积微相和层序地层等地质研究的基础上， 地震 积 体 系的控制 ，在 大庆长垣 南部 基本 上位 于两 大物源 信息 多参数 综合 评价 方法是 岩性 油气藏 勘探 的有效 方 的交会处 ， 沉积条件极为复杂。 主要发育三角洲平原分 法 。由于松 辽盆地 扶余油 层河道 砂体属 于三 角洲平 流 河道 和三 角洲分 流 间湾两种 微相 类型 ，分流 河道 在 原和前 缘分流 河道沉 积 ， 震的低 分辨率 ， 道砂 体的 纵 向上具有较明显的正旋 回特征 ，河道宽度一般在 地 河 薄互层、 横向变化快、 储层非均质性强等地质特点 ， 是 ３ ０ ０ ｍ， ０ ～６０ 扶余油 层下部河 道较 宽 ， 可达 ｌ０ｍ ， Ｏ０ 砂 当前 应用 地震技 术识 别扶 余一杨 大城 子油 层河道 砂体 岩单层厚度较薄 ， 一般为 ３ ｍ， ～５ 最厚达 ９ ｍ，电测 曲 的 困难 所在 ［ 近年 来地震沉 积学 的发展 为扶余 油层河 线主要表 现为钟 形 ， ３ ］ 。 其次 为箱 形 ， 曲线 大 多光 滑程 度较 道砂体分布研究提供了新的理论基础［ ¨ １ ］ 。本文在进 差 ，齿化较 明显 ； 流 间湾发育在 分流河 道之 间 ，以泥 分 行细 分小 层 层序地 层研究 分析 基础 上 ，应 用地震 分频 岩和粉砂 质泥岩 沉积 为主 ， 当三角洲 向前推进 时 ， 多被 处理 、各 向异性和叠前弹性波阻抗反演预测技术对河 决 口扇 、 洪泛沉 积充填 ， 测井 曲线呈 线性齿状 或微 齿状 道展 布和沉 积演化 进行 了描述 。 特征 ，有时可 见薄 层砂岩 ，厚度小 于 ｌ ｍ，电测 曲线呈
０ ，式 （）可简化为 ： １
ｌ ｌｐ）｛ ４ ［ｐ ２ 卢 ｓ。 ｎ ｎ ａ＋ｌ 一 １ ＋ ｌ ］ｉ ０ （ Ｅ （ ｎ ｎ ｎ ）ｎ ２ ）
储层参数预测
P1 = Vc1(V – Vc2) / V(Vc1 – Vc2) ……..3 式中Vc1、Vc2是组分速度
如已知A、B两井
A井目的层砂岩含量P1A、页岩含量P2A、层速度VA B井目的层砂岩含量P1B、页岩含量P2B、层速度VB 由2式得： Vc1 = VA P1AVc2 / (Vc2 -VA P2A )
2）尖灭型岩性油气藏
储集层沿上倾方向尖灭或渗透性变差而形成的油藏。
2、地震储层预测技术概述
地震勘探的最早梦想之一，就是期望在地面地震资料 中，提取地下直观的岩性信息。
当储层的岩性、储集层的物性（厚度、孔隙度、渗透 率等）、地层的纵向组合关系发生变化时，都会造成地 震反射特征发生相应的变化。当储层的这些变化达到相 应的限度时，将被记录在地震剖面上。
678000
678475 4246950
2100
2120 1739
4246000
1720
1720
4246000
1700
4245000
1680
1660
4244000
1640
1620
4243000
1600
1580
4242000
1560
1540
SH10 SH10
CB802 CB802
SHG3
CBG4 CBG4
注：不同含砂量a、b值是不同的 取纯砂、泥岩值时应取厚度2米以上的，且声波曲线无畸变，
以保证Vi的计算精度 应取足够多的井，至少在10井以上
实际应用见程序计算和量版的制作
b 散点法（精度较差） 在井少或无井区，
可用地震层速度作岩 性指数量板：
灰质背景下浊积岩储层地震响应特征及识别方法——以东营凹陷董集洼陷为例
灰质背景下浊积岩储层地震响应特征及识别方法——以东营凹陷董集洼陷为例于正军【摘要】近年来浊积岩油藏预测的成功率和勘探效益有所降低,灰质成分的存在是影响浊积岩储层预测精度的主要因素.从岩性组合、振幅、速度、波阻抗等特征入手,明确了东营凹陷董集洼陷含灰质成分岩石的地球物理响应特征,确定了灰质背景下浊积岩储层地震识别方法.利用子波分解技术优选适用子波重构新的地震道来消除灰质背景,以提高地震资料的分辨率;通过基于振幅和频率特征的多属性联合分析技术获得的幅频比能够有效压制灰质能量,放大浊积岩储层的信息并确定其边界;利用叠前反演技术获取纵、横波速度比和纵波阻抗来剔除灰质影响,进而确定浊积岩储层分布规模及与邻区储层的展布关系.【期刊名称】《油气地质与采收率》【年(卷),期】2014(021)002【总页数】3页(P95-97)【关键词】灰质背景;浊积岩;子波分解;幅频比;叠前反演;东营凹陷【作者】于正军【作者单位】中国石化胜利油田分公司物探研究院,山东东营257022【正文语种】中文【中图分类】P631.4随着隐蔽油藏勘探的深入，浊积岩油藏已成为东营凹陷重要的岩性油藏勘探类型，具有较大的勘探空间和资源潜力［1-5］。
但是，经过几十年的勘探开发，浊积岩油藏勘探难度越来越大，勘探效益明显降低。
通过对东营凹陷各地区的浊积岩油藏地震地质条件进行分析，发现影响浊积岩储层预测精度的主要因素有地震资料品质、岩性组合、古地貌等，其中岩性组合是浊积岩储层预测的关键因素之一［6］。
多年的勘探实践成果表明，东营凹陷董集洼陷的浊积岩包括砂岩与泥岩、砂岩与油页岩、砂岩与灰质泥岩、灰质砂岩与泥岩及灰质砂岩与灰质泥岩等5种岩性组合方式。
其中，灰质成分的存在是影响浊积岩储层预测精度的主要因素。
因此，笔者在前人沉积模式、成藏规律研究成果的基础上，通过分析浊积岩中各类岩石的地震响应特征，确定了灰质背景下浊积岩不同岩性组合的地震识别方法，以期为提高浊积岩油藏勘探成功率提供技术支持。
基于90°相移的振幅和频率属性融合法油气检测
基于90°相移的振幅和频率属性融合法油气检测ZHENG Jiangfeng;PENG Gang;SUN Jialin;ZHEN Zongyu【摘要】目前主要利用叠后地震资料的振幅和频率信息进行油气检测.甜点属性预测是其中常用的一种方法,但受储层厚度、相邻层干涉等多种非流体因素影响,检测结果存在多解性.为此研究了基于90°相移的振幅属性和频率属性进行油气检测的原理,提出了融合90°相移的振幅和瞬时频率属性进行流体检测的方法.该方法基于原始零相位地震资料提取90°相移的振幅和频率属性,再利用已钻井油气水层的地震属性约束最优化求解融合参数,最后将基于90°相移的振幅和频率属性融合从而获得油气检测结果.模型正演结果表明,相较于常规甜点预测方法,该方法在主力储层厚度范围内(8～25 m)受储层厚度影响小,能有效压制相邻层干涉,更适合流体检测.渤海PL油田的实际应用结果表明该方法能够有效区分油层与水层,且烃检异常空间分布符合油气成藏规律.【期刊名称】《石油物探》【年(卷),期】2019(058)001【总页数】9页(P130-138)【关键词】地震振幅;甜点;油气检测;90°相移;瞬时频率;储层厚度;储层干涉【作者】ZHENG Jiangfeng;PENG Gang;SUN Jialin;ZHEN Zongyu【作者单位】;;;【正文语种】中文【中图分类】P631随着油田勘探开发不断深入、难度不断加大,不仅需要精细刻画储层空间展布,还需要准确预测储层的含流体性质,从而使油田勘探开发能规避风险且进一步挖掘潜力[1]。
目前储层地震流体检测包括叠前检测和叠后检测,叠前流体检测主要是基于AVO理论利用测井资料和地震道集进行叠前流体参数反演,此类方法计算量大且耗时长,易受测井资料、道集质量、模型和反演方法等因素影响;叠后流体检测则主要基于“亮点”或吸收衰减理论,利用地震振幅或频率属性进行流体检测,此类方法计算量小且耗时短,更适合油田高效勘探评价[2]。
储层预测综述
储层预测综述一、序言储层是储集层的简称，在油气勘探生产中特指地下可供油气聚集、赋存的岩层。
通常从储层的岩性、形态、物性和含油气性四大方面对储层进行表征。
储层岩性是用来描述储层构成成分的要素，它直接或间接地反映了岩层的储集性能和储层特征，一般从储层的岩性、所处相带等方面描述，对于碎屑岩储层还常用砂地比(或砂泥岩百分比)来描述其储集性能;储层形态是对储层的几何形态进行描述的重要参数，常用的描述参数主要有储层的分布范围、储层顶界面构造形态、储层厚度等;描述储层物性参数主要是孔隙度和渗透率;储层含油气性描述主要包括储层是否含有流体、储层含流体的类型和含油气饱和度。
储层地震预测技术是以地震信息为主要依据，综合利用其他资料(地质、测井、岩石物理等)作为约束，对油气储层的几何特征、地质特性、油藏物理特性等进行预测的一门专项技术。
储层地震预测主要是通过分析地震波的速度、振幅、相位、频率、波形等参数的变化来预测储集岩层的分布范围、储层特征等。
岩性、储层物性和充填在其中的流体性质的空间变化，造成了地震反射波速度、振幅、相位、频率、波形等的相应变化。
这些变化是目前储层地震预测的主要依据。
在特定的地震地质条件下，只有这些储层特征参数变化达到一定程度时，才能在地震剖面上反映出来。
随着地震资料采集和处理技术的发展、地震资料品质的不断提高，这些特征参数的变化在地震剖面上的清晰度越来越明显，可信度也越来越高。
运用地震波的运动学特征确定地震波传播时间和传播速度，可以确定地层上下起伏变化的几何形态;而研究岩性时就必须运用波的动力学特征，结合运动学特征确定各种物性参数，来判断地层的岩性成分，以便寻找油气。
在储层预测中，储层的空间追踪和描述借助于提取出的储层的各种参数，包括纵波、横波速度、频率、相位、振幅、阻抗、密度、弹性系数、吸收系数及薪滞系数等。
根据这些参数的差异来分辨、识别、预测岩性，甚至油气层。
二、储层预测技术储层地震预测技术是一门方法繁多、综合性强、相互交叉的技术系列，单项技术不下数十种。
地震储层参数预测方法简述
意到其应用条件与局限性，优选参数，必将为油气勘探解决更多的
问题，提高储层预测的准确率。 地球物理学家希望利用地震资料解决以储层岩性、孔、渗、饱
等物性参数及以孔隙流体性质为内容的储层预测问题。 然而，尽管地震储层预测方法研究取得了较大发展，但 应当指出，就目前而言，它对储层研究的最基本需要还 不能完全满足，如对至关重要的孔、渗、饱三个物性参
年研究已发展了多种用地震资料求取孔隙度的方法。其主要方法大
致可分为四类：
1、用Wyllie时间平均方程或修正了的Wyllie公式求孔隙度
它是从地震速度出发，建立速度与孔隙度的函数关系，然后用 此函数关系求取无井处的孔隙度。此方法原理简单，应用方便，是
目前广为应用的方法之一。孔隙度是速度的多值函数，因此，仅根
也存在一定的缺陷：综合反演中利用的
A Fn
变换只持续
到 3 8 ，大于 3 8 的储层则需要借助其它方法；约束反演法应 用于横向上速度变化大的地区，容易产生误差。
（二）地震储层厚度预测方法 3、 CUSI网络法 目前常用的储层厚度预测方法适用条件不同，预测精度有 别，各有优点和局限性，考虑到地震储层厚度预测的复杂性及 特点（不同储层厚度对应不同的地震特征），与BP网络函数逼 近需要利用全体样本的信息、学习效率低等不足，提出了用完 全利用样本信息（缩写为CUSI：Complete Utilization of
厚度预测研究的不断深入，人们提出了多种储层厚度预测方法。
其主要方法大致可分为三大类：
1、单参数与多参数法 ：
（1）振幅图版法。此方法依据薄层和薄互层的地震振幅
与储层厚度呈线性关系，适用范围为薄层和薄互层。
（2）时差法。适用于储层厚度大于λ /4 （ λ 为波长）的厚 层。
基于储层反演定量预测技术对沉积相应用的研究
[收稿日期]2020-07-10[基金项目]中国石油化工集团有限公司科技攻关项目(P20060-1)[作者简介]杨怀宇(1982 ),男,安徽铜陵人,中国石化胜利油田分公司勘探开发研究院副研究员,博士,主要从事含油气盆地沉积学研究㊂doi :10.3969/j.issn.1673-5935.2020.04.001基于储层反演定量预测技术对沉积相应用的研究杨怀宇(中国石化胜利油田分公司勘探开发研究院,山东东营257015)[摘㊀要]㊀应用岩比法或相比法的传统沉积相研究方法具有不确定性强㊁精度不高的劣势,难以满足岩性油气藏边界的准确刻画㊂以草桥北坡王955井区为例,综合运用基于曲线重构技术㊁地震反演技术对井间和无井区砂岩进行定量预测,最终实现沉积相边界的准确刻画㊂研究表明,测井约束反演技术可较好地反映目标砂体的平面展布规律,特别是针对研究区滩坝砂薄储层的预测效果好,为最终实现草桥王955块井间及无井区沉积相的边界刻画提供保障㊂[关键词]㊀草桥北坡;沉积相;曲线重构;储层反演[中图分类号]TE122㊀㊀[文献标识码]A㊀㊀[文章编号]1673-5935(2020)04-0001-04㊀㊀基于钻井砂岩统计的传统单因素沉积相分析法对于勘探程度较高的地区适用性较好[1]㊂东营凹陷南坡乐安油田是勘探开发近50年的老油田[2],在石村断层储量集中区的北部斜坡带油气发现少㊁勘探程度较低,今年该区完钻的王955㊁王956等井在沙四上首次发现了滩坝砂油藏,但该区钻遇此种类型探井极少,传统的沉积相研究方法难以准确刻画砂体展布㊂因此如何精细刻画滩坝砂展布范围成为制约本区勘探的关键地质问题㊂笔者旨在利用地震反演技术实现井间和无井区砂岩进行定量预测㊂1㊀区域沉积特征草桥北坡沙四上早期主要发育来自广饶凸起的扇三角洲沉积,与前期的沙四下红扇沉积物相对应地称为 蓝扇 沉积,岩性组合以蓝灰色㊁浅灰色含砾砂岩㊁粉细砂岩为主,砂体单层厚度较大,一般3~10m 左右,垂向上表现为弱正旋回序列㊂近期,在草桥东坡完钻的王955㊁王956等井在沙四上纯下3砂组钻遇多套油层,岩性组合以灰色㊁浅灰色粉砂岩㊁白云质砂岩㊁泥岩及薄层泥灰岩沉积为主,垂向上以细 粗的反旋回或细 粗 细的复合旋回为主㊂结合邻井王73井区该套砂体的岩心资料来看,岩心上具有浪成砂纹交错层理(图1(a))㊁脉状层理,层面上发育垂直生物潜穴等沉积构造(图1(b))㊂同时,在显微镜下发现碎屑岩中发育薄皮灰质鲕粒㊁灰质内碎屑等特征(图1(c))㊂这些沉积特征具有滨浅湖岸边砂体被沿岸流改造的浅滩环境的特点[3-6],结合岩心薄片的累积粒度曲线特征,王73-斜1曲线表现为2段跳跃组分特征,反应岸边冲刷回流的水动力特征[6],由此认为王955-王73一带为扇三角洲侧翼的滩坝砂沉积相类型㊂自王955井滩坝砂成功之后,该井区向南抬高的高部位草36㊁草38等井储层基本不发育,随之亟需解决的关键问题在于滩坝砂周缘是如何展布,与周围 蓝扇 的扇三角洲砂体之间的关系一直难以认清,因此需对测井约束地震资料进行反演,确定该套砂体的展布范围㊂2㊀储层反演定量预测技术分析研究区地震资料发现,地震主频和有效频段12020年12月中国石油大学胜利学院学报Dec.2020第34卷㊀第4期Journal of Shengli College China University of PetroleumVol.34㊀No.4为19Hz和5~60Hz,若主频换成最大频率,H(厚度)=V(平均速率)/4ˑf(主频)=3000/240ʈ12.5m,而目的层段储层比较薄,储层厚度一般在10m以下,甚至小于5m,因此单纯靠地震资料无法准确刻画储层的分布㊂波阻抗反演是利用地震资料反演地震波阻抗处理技术,是以钻井㊁测井资料为约束对地下空间结构和物理界面识别并成像,主要应用于储层展布的预测,为油田的勘探和开发提供技术依据[7-9]㊂图1㊀草桥北坡王955井区沙四上沉积特征2.1㊀曲线重构技术常规波阻抗反演是采用测井曲线(主要为声波曲线)进行钻井约束,应用地震资料进行波阻抗反演,但在实际工作中声波资料不能较好地反映地下岩性的变化规律,因此应结合其他测井信息进行曲线重构研究[9-10]㊂首先,针对研究区测井资料分析研究,发现测井曲线需要进行异常点消除㊁深度校正和基线偏移等㊂在此基础上,由于测试仪器及围岩环境的不同造成测井响应值会相差较大,即同一套储层邻井之间波阻抗值会存在较大差异㊂因此,需进一步对曲线进行测井标准化处理,使泥岩层在各井中的响应值一致,标准化后的声波曲线及其他应用的曲线泥岩响应值处理为一致的区间㊂其次,通过草桥北坡沙四上亚段的测井资料分析认为,目的层段井声波曲线对岩性反映不敏感,而自然伽马曲线(GR)却能很好的反映岩性变化,能较好地区分储层与非储层(图2),因此本研究将采用拟声波阻抗反演,即声波时差与自然伽马拟合进行储层精细刻画㊂其原理主要考虑结合声波曲线低频信息(地层背景速度)与自然伽马曲线高频信息(对地层岩性变化敏感),重构拟声波曲线,该构建技术更合理,更能准确地预测储层㊂2.2㊀储层反演定量预测在全区井震标定基础上,按照层位标定结果开展严格精细的层位解释,通过不同频率储层预测对比来看,随频率的提高,其反演结果对储层的刻画能力有所提高,但整体仍是分辨率低,储层刻画效果不理想㊂针对此问题,对地震初始模型进行建模,利用分形建模的值在横向上对地震振幅的变异进行调整,结果更为合理㊂此外,针对采样率尝试发现,2ms采样率基本满足反演需求,能较好地识别砂体横向连续性及纵向展布,1ms采样率反演结果的纵横向分辨率相对高些㊂根据实际需要,可以结合2 ms采样率和1ms采样率反演结果㊂以草桥北坡过王955井的近北东向反演剖面可以看出,纯下3砂组目的层王955-王732储层均较发育,向南北至草38㊁王92井储层发育均较差,可见储层反演预测整体较为准确㊂图2㊀拟声波敏感曲线交汇图3㊀反演效果3.1㊀基于反演结果修正沉积相本次储层反演定量预测采取正演与反演相结合㊁无井约束属性与井约束反演相结合㊁常规阻抗反2第34卷㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀中国石油大学胜利学院学报㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀2020年㊀第4期演与特色阻抗反演相结合,最终达到降低储层识别多解性㊁提高储层分辨率的目的(图3)㊂同时,通过色标调整,提取纯下3顶至底时窗内的均方根振幅㊁变差函数等地震属性,并与井点统计的砂体厚度进行叠合修正(图4),达到无井区的储层定量预测的目的㊂从图4可以看出,草桥北坡整体 蓝扇 发育在西坡的草328-通20-官115一带,平面上呈北西走向,东边的王955井区呈近东西展布的滩坝砂,与西部 蓝扇 主体是分开的,由此修正了草桥北坡的沉积相展布图(图5)㊂图3㊀草桥北坡过王955井储层反演剖面图4㊀王955井区纯下3砂组反演砂体与钻井厚度叠合图3杨怀宇:基于储层反演定量预测技术对沉积相应用的研究图5㊀草桥北坡纯下3砂组沉积相平面图3.2㊀应用效果从应用效果上看,基于上述认识在王955井区滩坝砂岩性油藏的高部位部署了王956(油层4.2m/2层)㊁王955-斜1(油层12.7m/3层)㊁王955-斜2(油层13.5m/8层)分别钻遇多套滩坝砂油藏,表明沉积相相边界刻画的准确性,储层厚度与反演预测结果基本一致㊂因此在无井区沉积相的精细刻画应用储层反演预测技术具有较好的勘探效果,可推广应用㊂[参考文献][1]㊀冯增昭.单因素分析多因素综合作图法:定量岩相古地理重建[J].古地理学报,2004,6(1):3-19.[2]㊀才巨宏.乐安油田草4块沙三段:沙二段储层沉积特征与非均质性研究[J].油气地质与采收率,2011,18(3):24-28. [3]㊀阳孝法,林畅松,刘景彦,等.博兴洼陷沙四段滩坝沉积体系及其主控因素[J].油气地质与采收率,2009,17(1):51-56. [4]㊀张宇.东营凹陷西部沙四段上亚段滩坝砂体的沉积特征[J].油气地质与采收率,2008,15(6):35-38.[5]㊀杨勇强,邱隆伟,姜在兴,等.陆相断陷湖盆滩坝沉积模式:以东营凹陷古近系沙四上亚段为例[J].石油学报,2011(3):423.[6]㊀李安夏,王冠民,庞小军,等.间歇性波浪条件下湖相滩坝砂的结构特征:以东营凹陷南斜坡王73井区沙四段为例[J].油气地质与采收率,2010,17(3):12-14.[7]㊀冯凯,查朝阳,钟德盈.反演技术和频谱成像技术在储层预测中的综合应用[J].石油物探,2006,45(3):262-266. [8]㊀罗士利,罗中华,黄景秋.随机反演技术在薄互层储层预测中的应用[J].吉林大学学报(地球科学版),2006,36(专辑): 58-61.[9]㊀鲜本忠,姜在兴,杨林海,等.测井约束反演在泌阳凹陷下切谷储层预测中的应用[J].石油大学学报(自然科学版),2003,27(5):8-15.[10]㊀杨晓兰,刘小平,赵锡桥.曲线重构反演在火成岩储层预测中的应用[J].西安石油大学学报(自然科学版),2007,22(6):32-35.[责任编辑]㊀胡秋媛4第34卷㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀中国石油大学胜利学院学报㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀2020年㊀第4期。