一种运用Fischer图解划分砂砾岩沉积期次的方法

摘要

本发明涉及一种利用Fischer图解划分砂砾岩体沉积期次的方法。包括选取测井曲线组合、判断岩性变化面、判断高频基准面旋回、确定高频基准面旋回优势周期、判断高频基准面旋回划分的准确性、制作Fischer图解、划分长期基准面旋回七个步骤。通过分析岩相建立岩相与测井曲线之间的对应关系，选取对岩性反应良好的测井曲线组合进行聚类分析、主成分分析、多项回归分析判断岩性变化；运用测井曲线组合、岩性变化面、离散小波高频成分划分高频旋回基准面，并采用时频分析技术获得测井曲线的高频优势旋回周期，对所划分的高频旋回基准面加以检验。在此基础上上，做出Fischer图解得到反应基准面变化的曲线，划分沉积期次。

说明书

技术领域：

本发明涉及一种砂砾岩体沉积期次划分方法。

背景技术：

砂砾岩油藏，横向上非均质性强，纵向上沉积厚度大，泥岩隔夹层少，古生物标志及标志层缺 乏。、测井响应不明显、地质录井粗糙、地震分辨率低，造成了单井划分与全区对比困难，砂砾岩体沉 积期次的划分对比已成为制约砂砾岩体油气藏勘探与开发的瓶颈之一。

小波变换中，对于非平稳信号而言(测井信号即为非平稳信号)，需要时频窗口具有可调的性质，即 要求在高频部分具有较好的时间分辨率特性，而在低频部分具有较好的频率分辨率特性。对测井曲线进行 小波分析可以识别测井曲线中不同频率的曲线旋回，高频的曲线旋回对应高频的即短期的沉积旋回，低频 的曲线旋回对应低频的即长期的沉积旋回，因此可以用不同频率的测井曲线旋回划分不同周期的沉积旋 回，对应不同级别的层序地层单元。

Fischer图解方法是在对测井曲线高频基准面旋回进行划分的基础上，通过研究沉积旋回在空间上的 叠置规律、沉积基准面变化态势，来为建立沉积基准面变化曲线提供一种客观使用的方法。在砂砾岩地区 制作Fisher图解的基本方法为：(1)运用测井曲线组合、岩性变化面、离散小波高频成分精确划分高频 旋回基准面；(2)计算旋回层的平均厚度(算术平均)；(3)计算各旋回层的厚度偏差和累积厚度偏差；(4)以 横坐标为旋回数，纵坐标为累积偏差，将各旋回数对应的累积偏差连线即得Fischer图解。

发明内容：

本发明涉及一种利用Fischer图解划分砂砾岩体沉积期次的方法。包括选取测井曲线组合、判断 岩性变化面、运用测井曲线组合、岩性变化面、离散小波高频成分判断高频基准面旋回、确定高频基准面 旋回优势周期、判断高频基准面旋回划分的准确性、制作Fischer图解、划分长期基准面旋回七个步骤。 通过分析岩相建立岩相与测井曲线之间的对应关系，选取对岩性反应良好的测井曲线组合进行聚类分 析，主成分分析与多项回归分析判断岩性变化；运用测井曲线组合、岩性变化面、离散小波高频成分划 分高频旋回基准面，并采用时频分析技术获得测井曲线的高频优势旋回周期，对所划分的高频旋回基 准面加以检验，反复修改划分结果，所划分的高频基准面旋回平均厚度与高频优势旋回周期所对应的 厚度一致。在此基础上上，做出Fischer图解得到反应基准面变化的曲线，划分沉积期次。

关键技术要点包括：

(1)识别砂砾岩体岩性变化界面

运用岩心、成像测井资料建立岩性模式，选择对岩性反应良好的测井曲线进行聚类分析、主成 分分析和多项式回归分析判断岩性变化界面。

(2)划分高频基准面旋回

运用测井曲线组合、聚类分析曲线、主成分回归曲线判断高频基准面旋回，当测井曲线组合、聚 类分析曲线、主成分回归曲线叠加趋势发生改变时，认为存在高频基准面旋回。当基准面存在于砂砾 岩内部时，以上组合不具有明显特征，此时对主成分分析所得到的第一主成分进行离散小波分解，当 小波高频成分发生明显的频率变化时，认为存在高频基准面旋回。

附图说明

图1.砂砾岩沉积期次划分流程图

图2.岩性变化面识别示意图

图3.高频旋回基准面划分示意图

图4.高频基准面旋回优势周期能量图

图5.Fischer图解图

具体实施方式：

如图1所示，方法的实施步骤详细介绍如下：

步骤一，运用岩心、成像测井资料建立岩性模式，对不同测井曲线分别做出岩性—测井散点图，通 过观察测井曲线与岩性之间的对应关系选取对岩性反应良好的测井曲线组合；

步骤二，对所选测井曲线组合进行K中心聚类分析，分类数K与岩性数量一致。运用主成份分析对 测井曲线组合降维，选择累计贡献率达到90％以上的m个主成分进行曲线回归分析。当聚类分析曲线呈现 阶梯状变化且主成分回归曲线出现拐点时，判断为岩性变化面(图2)；

步骤三，运用测井曲线组合、岩性变化面判断高频基准面旋回，当测井曲线组合、岩性变化面叠加 趋势发生改变时，认为存在高频基准面旋回。当基准面存在于砂砾岩内部时，以上组合不具有明显特征， 此时对主成分分析所得到的第一主成分进行离散小波分解，当小波高频成分发生明显的频率变化时，认为 存在高频基准面旋回(图3)；

步骤四，对测井曲线运用时频分析技术得到高频基准面旋回优势周期(图4)；

步骤五，判断所划分的高频基准面旋回平均厚度是否符合高频旋回优势周期；若不符合，重复步骤 四直到符合为止，若符合则进行下一步；

步骤六，运用划分好的高频基准面旋回制作Fischer图解(图5)；

步骤七，运用所得到的Fischer图解确定长期基准面旋回，在Fischer图解出现极值时，说明A/S的 (可容纳空间/沉积物供给)变化趋势发生改变，可在极值点处划分长期基准面旋回，确定沉积期次。