技术摘要：
本发明公开了一种砂岩储层内部阻渗条带的预测方法，包括以下步骤：地震资料品质评估；地震资料去噪处理；三维数据体的局部结构熵的计算；预测结果的标定；阻渗条带的展布预测；油田实际生产中井位的部署及优化。本发明采用局部结构熵算法计算得到剖面图和平面图，克服  全部
背景技术：
油田开发过程中出现的层间干扰、单层突进以及注采不均衡等问题都与油气储层 的非均质性有关，因此油气储层的非均质性是制约油田开发效果的重要因素之一。油气储 层内部能够对流体的流动、运移或者聚集产生影响的不渗透或低渗透条带被称为阻渗条 带，阻渗条带是导致储层非均质性的重要因素之一。砂岩储层是最常见的油气储层类型之 一，我国80％以上的油气储层均为砂岩储层。因沉积作用、成岩作用和构造运动等导致油气 储层的岩性和物性变化及小断层形成砂岩储层内部阻渗条带，阻渗条带对流体渗流产生不 同程度遮挡作用，是导致油田开发注采不受效、储量采出程度低和局部剩余油富集的重要 原因。 在地震信号处理中，局部结构熵算法(Local  Strcctural  Entropy，简称LSE)反映 了道与道之间信号的不连线性(参考煤田地质与勘探期刊2007年35卷第1期论文《局部结构 熵值在地震数据不连续性检测中的应用》)。该方法在研究地震资料的相关性问题中有所应 用，但未见其应用于预测砂岩内部阻渗条带的分布。 砂岩储层内部阻渗条带在指导注水井和采油井的设计和施工的方面具有重要意 义，是油气开发的重要数据。但由于砂岩储层内部阻渗条带的尺度通常小于地震资料的分 辨能力的尺寸范围，通过常规的地震反演或储层预测手段很难对其进行准确识别，目前不 能对砂岩储层内部阻渗条带进行准确预测。
技术实现要素：
本发明提供一种砂岩储层内部阻渗条带的预测方法，用以解决砂岩储层内部阻渗 条带尚且不能进行预测的问题。 本发明提供一种砂岩储层内部阻渗条带的预测方法，所述砂岩储层内部阻渗条带 的预测方法包括以下步骤： 地震资料品质评估：分析目标砂岩储层对应地震资料空间范围的信噪比和分辨 率，对地震资料品质进行评估； 地震资料去噪处理：采用去噪或滤波的手段优化地震资料品质，直至提高地震资 料品质中的信噪比和地震分辨率； 三维数据体的局部结构熵的计算； 预测结果的标定； 阻渗条带的展布预测：根据目标油气储层的顶底解释层位，利用局部结构熵体值 计算顶底界面之间的局部结构熵值的总和，并显示为平面图，该平面图即可用于表示储层 内部阻渗条带的平面分布情况； 4 CN 111580158 A 说　明　书 2/5 页 油田实际生产中井位的部署及优化：基于砂岩储层内部阻渗条带的平面展布认 识，合理指导油田实际生产中井位的部署及优化，生产井的井位需要避开阻渗条带，注水井 和采油井之间避开阻渗条带。 优选地，所述三维数据体的局部结构熵的计算包括以下步骤： 采样点的预处理：对原始三维地震数据体的采样点进行预处理，将每个采样点的 地震数据转化为该数据减去该地震道均值得到的地震数据差，地震数据差的计算公式为： 式中，dxyt和 分别为原始地震数据和转化后的地震数据，Nt为每一地震道的地 震采样点总数，Et表示每一地震道各采样点的均值，x、y为地震数据的空间坐标位置，t表示 地震数据的时间变量； 协方差矩阵的计算：针对要计算的样本点，选取以它为中心的2×L1个线数，2×L2 个道数和N个时间样点构成三维分析时窗，将三维分析时窗分为四个L1×L2×N的象限，在每 个象限体内将地震道依次首尾连接，合并成一个向量，得到对应的四个向量{ai|i＝1,2,3, 4}，从而得到计算样本点对应的协方差矩阵S，表示为： 样本点地震数据的局部结构熵的计算：利用构建的协方差矩阵，计算样本点对应 的局部结构熵值ε，计算公式为： 式中，||·||为Hilbert-Schmidt算子，tr指计算矩阵的迹， 整体地震数据的局部结构熵值：将三维分析时窗在三维数据体中滑动，得到三维 数据体每一个样本点处的局部结构熵值，从而得到整体地震数据的局部结构熵值； 优选地，所述预测结果的标定包括以下步骤： 整理目标油气储层的所有已钻水平井上的泥岩段数量和厚度； 剔除由于非地质因素钻遇的泥岩段信息； 根据地震资料能分辨的储层厚度，筛选用于整个地震数据体的局部结构熵值标定 的泥岩段厚度下限； 利用整体地震数据局部结构熵值过水平井剖面，将其中局部结构熵值较大的异常 值与筛选后的水平井泥岩段进行标定，得到单井吻合率；若全水平井段未钻遇泥岩段的水 平井，若局部结构熵值过井剖面且无较大的异常值，则吻合率为100％，否则吻合率为0％。 本发明的有益效果是： 本发明通过计算地震数据的局部结构熵值，并将局部结构熵值较大的异常值与水 5 CN 111580158 A 说　明　书 3/5 页 平井上钻遇的泥岩段进行标定，克服了砂岩储层内部阻渗条带的尺度通常小于地震资料的 分辨能力的尺寸范围，通过常规的地震反演或储层预测手段很难对其进行准确识别的难 度，解决了砂岩储层内部阻渗条带不能进行预测的问题。 附图说明 图1(a)是本发明提供的去噪处理前的某油田的地震剖面图； 图1(b)是本发明提供的去噪处理后的某油田的地震剖面图； 图2(a)是本发明提供的某油田的地震剖面图； 图2(b)是本发明提供的某油田的局部结构熵剖面图； 图3是本发明提供的某油田水平井上的泥岩段数量和厚度的统计表； 图4(a)是本发明提供的水平井A10H的轨迹钻出储层的泥岩段的剖面图； 图4(b)是本发明提供的水平井A10H上的泥岩段数量和厚度的统计表； 图4(c)是本发明提供的水平井A10H上的岩性分布图； 图5是本发明提供的筛选后得到的用于标定的水平井泥岩段统计表； 图6(a)是本发明提供的水平井A14H的局部结构熵值的过井轨迹剖面； 图6(b)是本发明提供的水平井A14H的泥岩段数据统计表； 图6(c)是本发明提供的水平井A15H的局部结构熵值的过井轨迹剖面； 图6(d)是本发明提供的水平井A13H的局部结构熵值的过井轨迹剖面； 图6(e)是本发明提供的水平井A13H的泥岩段数据统计表； 图7是本发明预测的砂岩内部阻渗条带的平面投影图。