GF4地震属性

GeoFrame 地震属性分析
内容
GeoFrame 地震属性工具包(SATK) GeoFrame 分频技术 GeoFrame 地震属性应用流程
– 地震相分析流程(SeisClass) – 测井物性推广流程(LPM) – 用分频属性估算净产厚度流程
引言
据美国能源部的统计资料,勘探钻空率从1996年的67%下降 到了2002年的55%. 地震属性的提取和分析技术在探井成功率的显著增长中扮 演了重要角色. 当今的有效技术手段中,实时模型技术提供了最完善的地 震属性提取和分析应用的解决方案.
LPM的优点
只使用井资料获得的净产图（Net Pay Map） 使用井资料与地震资料相结合获得的净产图
LPM 的结果能更精确地表达地质 变化，帮你确定储层的走向
LPM 的流程
岩石物理评价
分层
储层物性集总
地震属性提取
相关分析
物性分布图
LPM 的流程
岩石物理评价
分层
储层物性集总
地震属性提取
相关分析
SeisClass 的流程
提取地震属性
地震属性分析
非督导聚类
训练数据选择
督导聚类
地震相平面图的 解释
SeisClass 的流程
提取地震属性
地震属性分析
非督导聚类
训练数据选择
督导聚类ቤተ መጻሕፍቲ ባይዱ
地震相平面图的 解释
SeisClass 的流程
提取地震属性
地震属性分析
非督导聚类
训练数据选择
督导聚类
地震相平面图的 解释
IESX 的体反射谱（VRS ）属性预览工具
GeoViz Virtual Volumes
什么是分频?
著名的Pink Floyd摇滚乐队的徽标 – 月亮的背面(1973)
我们借用该图像简捷地表达分频的处理过程
什么是分频?
10 Hz
20 Hz
30 Hz
40 Hz
50 Hz
GeoFrame 的分频属性
粉砂岩<20% 泥岩<8% Net/Gross Ratio 70-100%
浊流岩相储集层 叠合扇扇前 扇前 超薄扇扇端 无储层超薄扇扇端
为什么要把测井资料与地震属性进行定量相关？
测井 直接度量岩石的物性 覆盖面太小 三唯地震 间接度量岩石的物性 很好的横向覆盖面 (~ 25 meters)
LPM
Log Property Mapping
估算净厚度流程
分层 地震解释
求取净厚度
上半周持 续时间 下半周持 续时间
上半周 面积 下半周 面积
轨迹属性类（Local Attributes）
在地震剖面上显示的属性值曲线有助于理 解属性变化与地震波形之间的关系
轨迹属性类（Local Attributes）
在地震剖面上显示的属性值曲线有助于理 解属性变化与地震波形之间的关系
轨迹属性类（Local Attributes）
余弦相关变换（CCT）算法
Schlumberger R&D专利算法，该算法
用某一特定频率的非常窄的余弦函 数与地震数据进行互相关。
结果是一个某一特定频率的 相关数据体。用于揭示薄层 岩性变化，指示可能的含烃 地层圈闭。
河道的形状更清 晰，河道内的岩 性细节变化
沿河道顶面层位提取的25 Hertz 频率切片 – Offshore Angola
全频段总能量（Spectral Total Energy）
能量
频率
总能量: 指从用户定义的低截频到高截频之间的能谱的总面积
能谱的十分位属性（Spectral Deciles (1-9)）
能量
频率
1分位频率: 该频率之 8分位频率: 该频率之
前的面积占总面积的
10%
前的面积占总面积的
10%
频谱的十分位是指能谱图上位于总面积 10%、20%、30%。。。处的频率值
密西西比河三角洲（照片）
凸显河道与沉积的关系
等频体（IsoFrequency)
时间偏移体 10 Hz 等频时间体
一旦确定了某一关键频率，可以处理 一个该频率的时间或深度数据体
地质与地球物理的结合
SATK给LPM和SeisClass提供了关键性的数据
SeisClass用神经网络工具分析各种地震属性、划分各种相分布，为钻井 提供最有利的目标。 LPM 是把地震属性和测井物性进行定量相关，获得从井点推广开的储 层物性横向分布图。
按15美元/桶计， = $10.61 M B块:
可采储量= 200000 m² * 50 m * 90% Sand * 10% 空隙度 * (1-50%) * 30% 可采储量= 135000 cubic meters = 849150 桶
按15美元/桶计， = $12.74 M
地震相应用实例 – 巴西海上
地震与地质数据的相关
用已知井的地质数据作为 地震相划分的控制依据
气/ 砂 水/ 砂 泥岩
地震与地质数据的相关
用已知井的地质数据作为 地震相划分的控制依据
气/ 砂 水/ 砂 泥岩
SeisClass 的优点
地震相平面图 砂层概率平面图
其结果可提供 理想的井位 + 潜在的风险
应该钻什么部位?
例如，有二个潜在的有利 地区，但今年只剩一口井 的预算。A块是B块面积 的1.5倍；假设地震上看 到的目的层厚度相似，我 们自然会选择钻A块。
物性分布图
LPM 的流程
岩石物理评价
分层
储层物性集总
地震属性提取
相关分析
物性分布图
LPM 的流程
岩石物理评价
分层
储层物性集总
地震属性提取
相关分析
物性分布图
LPM 的流程
岩石物理评价
分层
储层物性集总
地震属性提取
相关分析
物性分布图
LPM 的流程
岩石物理评价
分层
储层物性集总
地震属性提取
相关分析
物性分布图
有什么好处?
对IESX 的好处 性能提高5倍 (5X) 各种体反射谱属性(VRS Attributes) 独立的应用工具包能使用户在解释的同时运 行属性提取各种运算 取消了一次只能运算20 种属性的限制
反(负)二阶微商变换
右侧的反二阶微商变换显著地提升了连 续性,有助于更快、更准确的层位解释
地震道积分
SeisClass 的流程
提取地震属性
地震属性分析
非督导聚类
训练数据选择
督导聚类
地震相平面图的 解释
SeisClass 的流程
提取地震属性
地震属性分析
非督导聚类
训练数据选择
督导聚类
地震相平面图的 解释
SeisClass 的流程
提取地震属性
地震属性分析
非督导聚类
训练数据选择
督导聚类
地震相平面图的 解释
左: SP 曲线 右: 密度曲线
道积分能起到伪波阻抗剖面的作用. 并不是说用它替代反演, 它可以起到 快速指示空隙度变化的作用.
轨迹属性类（Local Attributes）
下半周持续时间 上半周持续时间 下半周面积 上半周面积 上半周偏斜度 峰度 不对称度
轨迹属性类（Local Attributes）
B
投影面积 200000 m² 投影面积
A
300000 m²
使用 SeisClass后，情况会怎么样?
砂岩 砂泥岩 泥岩
90% 砂岩
50% 砂岩
快速储量计算
可采储量 = 面积 * 厚度 * 净厚度 * 空隙度 * (1-含水饱和度) * 采收率
A块:
可采储量= 300000 m² * 50 m * 50% Sand * 10% 空隙度 * (1-50%) * 30% 可采储量= 112500 立方米 = 707625 桶
RMS 1250 - 1350 ms
Repeated windows work as a quick exploration scanner of your seismic volumes
等分重复时窗实例
等分时窗下的RMS振幅提取
Proportional repeated windows adjust seismic attribute extraction to the local stratigraphy
高级刻度算法中有地质统 计学类算法和神经网络法
估算净厚度流程
估算净厚度流程
分 层
用SP曲线划分砂层的顶、底 注意W-12井的净厚度明显较大
估算净厚度流程
分层
求取净厚度
在ResSum里，使用SP曲线， 设定砂岩门坎值，累积砂层段 的净砂层厚度。 右边是各井点处净厚度在底图 上的显示。 注意：W-4井的净厚度最大。 等值线是对应砂层的反射时间
GF 4.0.4.2的新地震属性处理工具
三个新的 GeoFrame 摸块 地震属性工具包(SATK) GF 分频模块（CCT） 属性浏览工具 - SeisClass Explore
地震属性工具包(SATK)
地震属性工具包 (SATK)的特点 独立的应用工具 GF 4.0.4.2 从地震工具总目录下起动
岩石物理与地球物理相关
地震属性 岩石物性
快速指示两个变量 具有高度相关性
岩石物理与地球物理相关
为什么GQ-4显得有点异常?
有特殊的地质原因吗?例如 GQ-4井与其它井之间是否被 断层隔开? 测井数据是否有问题? LPM 能让你更好地理解该地 区地震属性的响应
用井数据标定(刻度)地震属性
LPM 严格遵守井资料 剩余误差校正法是一种快 速、精确的刻度方法
能谱的十分位属性（Spectral Deciles (1-9)）
频/宽比（Spectral Band Ratios (1-10)）
频/宽比 SR = Ai / At Ai = 某一频率之前的面积 能量 At = 总面积
频率
总共可以求10个频/宽比. The first five ratios are equally divided between the “Frequency at Low End of Spectrum” and the “Center Frequency”.
• • • 各种地震属性计算 轨迹(Local)属性, 持续时间属性(Duration), 转译属性(Cypher)和 各种谱属性 体反射强度属性(VRS)
GeoFeature Mapping 许可证 – 3种属性 方差体许可证 – 1种属性
有什么好处?
并行处理 几十种全新属性 自动时窗移动扫描 分频 属性在GeoViz里预览
在地震剖面上显示的属性值曲线有助于理 解属性变化与地震波形之间的关系
上半周面积
横穿上半周面积为低值 区的地震剖面
上半周面积
横穿上半周面积为高值 区的地震剖面
上半周面积属性与振幅属性的比较
上半周面积
地震振幅
持续时间属性类（Duration Attributes）
与轨迹属性类似, 只是使用了时窗 负半周平均持续时间 正半周平均持续时间 平均持续时间 最小周期持续时间 最大周期持续时间 持续时间的标准偏差
时窗间隔内属性类（Interval Attributes）
振幅标准偏差 时差厚度 平均负振幅 平均正振幅 波峰平均值 波谷平均值 最小振幅时间 最大振幅时间
转译属性类（Cypher Attributes）
转译属性类（Cypher Attributes）
层位之上波峰个数
频谱类属性（Spectral Attributes)
频/宽比（Spectral Band Ratios (1-10)）
频/宽比 SR = Ai / At Ai = 某一频率之前的面积 能量 At = 总面积
频率
总共可以求10个频/宽比。The last five are equally divided between the “Center Frequency” and the “Frequency at High End of Spectrum”.
频/宽比（Spectral Band Ratios (1-10)）
重复时窗属性扫描（Repeated Windows）
等时重复时窗实例
RMS 1000 - 1100 ms
RMS 1050 - 1150 ms
RMS 1100 - 1200 ms
RMS 1150 - 1250 ms
RMS 1200 - 1300 ms
Data courtesy of WesternGeco
怎样运行 GF Spectral Decomposition?
底图上显示的沿层切片（25 Hz）
频率体的动画显示
动画显示是非常重要的分析手段
注意不同的频率突出了不同的地质现象
地质与地球物理的结合
墨西哥湾数据体 25Hz频率体与瞬时相位体的确叠合显示
属性浏览工具书 - SeisClass Explore
功能: 给用户提供一个价格低廉的浏览和 交汇分析多属性的工具
地震属性工具包(SATK) – 可供用户选择的选项
SeisAN 许可证 – 12 种属性
• 地震处理和属性数据体
分频属性 许可证 – 2种属性
• 余弦相关变换分频属性和等频属性
CSA 许可证 – 70多种属性