摘要: 最大熵谱分析技术是对测井曲线进行合成预测误差滤波分析处理，得到频率趋势线(INPEFA 曲线)，运用 INPEFA 曲线不同幅度的趋势拐点来识别不同级次层序界面。通过对苏北盆地溱潼凹陷南华区块阜三段 GR 曲线进行 INPEFA 处理，将阜三段 III 和 IV 油组划分为 13 个 5 级层序，进而在 INPEFA 曲线所建立的等时框架内对区块内距离相近的 3 口井进行井间的砂体连通性对比，得出在小层中各有 2 个连通性较好的砂层和油层组。INPEFA 技术处理后的测井曲线可以提取出在常规测井曲线中隐藏的旋回趋势特征，从 3 口井变化趋势相近的 INPEFA 曲线来看，该区域的沉积环境变化大致相同，都是从水退到水进的过程。
关键词: 最大熵谱分析;INPEFA;层序划分;溱潼凹陷;苏北盆地
层序地层学自从发展成一门独立学科以来，以其理论上的系统性、综合性，尤其是应用上的可预测性，为油气田的勘探开发作出了重要贡献，并产生了多个学派分支，其中就有以 GALLOWAY 为代表的以最大洪泛面为层序边界的成因地层学。在划分层序界面时，测井曲线是一项重要的识别依据[1]，而在溱潼凹陷南华区块，因缺少地震资料，取心资料有限且不连续，露头资料少，故使得运用测井资料进行层序划分变得十分重要。但由于使用常规测井资料的极限性，且容易穿时，所以在取心资料有限、缺少地震资料的条件下，可以借助 INPEFA 技术，对常规测井资料进行特殊的数学处理，将原始测井曲线中隐藏的地层发育趋势以及沉积旋回等信息提取出来[2]，通过对INPEFA曲线不同级别趋势形态和拐点类型的分析，确定不同级别的层序界面和洪泛面，以及判断沉积环境的变化，进行精确的高频层序对比。
本文对苏北盆地南华区块阜宁组三段测井曲线进行 INPEFA 处理，提取测井曲线所隐藏的周期性特征，完成了对阜三段四、五级高频成因层序的精确划分，判断出该时期的沉积环境变化情况，进而对局部地区相邻的 3 口井进行了精确的高频成因层序及砂体连通性对比。
1 区域地质概况
南华区块所处溱潼凹陷是苏北盆地东台坳陷的一个次级构造单元，南以断阶带与泰州凸起相接，东以小海、梁垛凸起与海安凹陷相隔，北与博镇—吴堡低凸起相接，呈南断北超的箕状特征(图 1)，凹陷面积约 1 200 km2 。
在 2010 年，通过对苏北华东探区资源潜力的重新评价，得出溱潼凹陷油资源量为 10 800×104 t，分布于溱潼凹陷断阶带、东部斜坡带、西部斜坡带上;油气资源在层位上主要分布于阜宁组三段，其中溱潼凹陷阜三段资源量占溱潼凹陷总资源量的 29.86%。所以对阜三段进行精细的层序划分及对比，进而分析该层位的砂体连通性及其中的沉积环境变化，对南华区块的油气资源勘探开发具有重要意义。
2 INPEFA 技术用于层序划分与对比
2.1 基本原理
INPEFA 技术的第一步是进行频谱分析，首先利用最大熵谱分析对测井曲线进行处理，来推算出最大熵频谱分析估计值(MESA)，再由对应的测井曲线真实值(RV)减去最大熵频谱分析估计值得到数据差值(PEFA = RV-MESA)，这个过程就是预测误差滤波分析，对 PEFA 曲线进行积分处理，得到了合成预测误差滤波分析曲线 ( INPEFA)[3]。 INPEFA曲线的正负拐点代表了可能的洪泛面和层序界面。
INPEFA 曲线可以反映水进水退的旋回性变化，而常规测井曲线旋回性变化不明显。INPEFA 曲线的不同趋势代表了不同的地质意义(表 1)[4]，从右到左的负向趋势反映基准面下降、水体变浅、砂质含量增加的反韵律变化，正向拐点(曲线由负趋势变为正趋势的转折点)代表可能的层序界面; 从左到右的正向趋势反映基准面上升、水体变深、泥质含量增加的正韵律变化，表示水进或洪积阶段，负向拐点(曲线由正趋势变为负趋势的转折点)代表可能的洪泛面。
据南华区块阜三段的录井资料，该地区的岩性特征较为单一，均为砂泥岩互层，而在反映泥质含量的变化方面，GR 曲线是测井曲线中效果最显著的，所以本文使用 GR 曲线来进行层序地层研究[5]。
2.2 利用 INPEFA 技术对阜三段进行层序划分
利用INPEFA测井旋回技术进行单井层序地层的识别和划分时，首先确定大级别的层序界面，然后在界定的层序界线内进行更为精细的层序划分。
2.2.1 用 INPEFA 曲线识别阜三段 III、IV 油组界面
INPEFA 曲线不同的变化幅度对应不同级次的旋回，用 INPEFA 技术对南华区块的南华 2-17 井做合成预测误差滤波分析(图 2)。从图 2 中可以得出，频率趋势线旋回形态特征明显，图中的红色箭头代表正向趋势，黑色箭头代表负向趋势。从整体上来看曲线呈现出单一的正趋势，从 III、IV 小层内来看，曲线呈现出由负到正的趋势。
所以从图 2 中可以清楚地识别出 1 768.13 m 和 1 733.47 m 出现相同变化幅度的负拐点，为 Ef3 -III 油组的底部和顶部。界面在 INPEFA 趋势线上非常明显，充分体现了 INPEFA 技术精确度高、准确性强的优点。
2.2.2 阜宁组三段 III、IV 油组精细地层划分
根据阜三段 III、IV 油组的测井、录井资料以及 INPEFA 曲线特征，在 III、IV 油组界线的控制下，结合沉积相序列的演化特征，进行精细的五级层序划分。INPEFA 曲线的拐点预示着沉积环境的变化，由图 2 可知 III、IV 油组的划分界线都是在曲线较大方向上的改变之处。所以运用同样的划分原理，在确定的四级层序界限内，通过识别较小的拐点，分别将阜三段的 III 和 IV 油组(4 级层序)划分为 6 个小层和 7 个小层(5 级层序) (图 3)。图 3 中较大的拐点(粗线)表示 4 级层序的界面，较小的拐点(细线)表示 5 级层序的界面。
2.3 利用 INPEFA 技术对阜三段进行等时对比
在单井高频层序划分的基础上，利用测井旋回曲线的变化特征进行连井高频层序的对比。只有对全区的层序进行划分与对比，并最后闭合，才能达到精细等时地层格架的建立[6]。
将南华区块内 3 口距离相近的井进行精细地层划分与连井对比(图 4)。可以看出 3 口井的趋势线和拐点都较为一致，在 INPEFA 曲线上 4、5 级层序界线都具有良好的可对比性，表现出的旋回特征都明显相似。从测井旋回曲线的特征可以看出，3 口井的目的层在该时期的变化趋势均为负到正，沉积环境大致相近，都是一个从水退到水进的过程。
2.4 利用 INPEFA 技术进行砂体连通性对比
INPEFA曲线能表现出砂泥岩的变化，且具有受区域气候影响的等时性特征。油层的划分是按照砂岩含油对比来进行的，但是在实际的油气田开发中，在同一个沉积环境形成的地层中进行砂泥岩互层对比的时候，只利用测井曲线的形态进行砂对砂、泥对泥的对比，容易造成邻井之间的油层出现穿时现象。因此，利用 INPEFA 测井旋回曲线进行邻井之间曲线形态的特征对比，能在小层等时的框架内划分砂、油层，进而精确地完成砂体连通性的识别[7-9]。
溱潼凹陷南华区块阜三段的岩性较为单一，均为砂岩、泥岩、砂泥岩。图 5 为溱潼凹陷南华区块阜三段 III、IV 油组砂体的连通性对比图，可以看出，在划分的 13 个小层中，只有 III - 4、III - 6、IV- 2、 IV-7 这4 个小层里的砂体是全部连通的，分别为 2 个砂层和油层。在通过 INPEFA 处理后的趋势线建立的等时框架下进行砂体连通性的分析，所得成果与该区块的实际动态开发资料也基本保持一致。
相关期刊推荐：《石油实验地质》(双月刊)曾用刊名：(石油地质实验)本刊是一份报道我国油气普查、勘探及开发成果，反映我国石油地质研究水平的学术性、技术性刊物。创办于1963年，现为双月刊，向国内外公开发行。本刊由中国石化石油勘探开发研究院和中国地质学会石油地质专业委员会主办，中国石化无锡实验地质研究所承办。有投稿需求的，可以咨询期刊天空在线编辑。
3 结论
(1)INPEFA 技术处理得到的趋势线可以很明显地识别出在常规测井曲线上无法识别的旋回趋势特征，且不同的变化幅度对应不同的层序级别，避免了常规人工划分时精度不高、容易穿时、多解性强等问题，使各层序界面更直观、更明显。
(2)由 INPEFA 曲线的旋回趋势确定了溱潼凹陷南华区块阜三段 III、IV 油组的界面，从 4 级层序的曲线变化趋势来看，沉积过程都是从水退到水进，在其框架内进行更精细的层序地层划分，将 III 油组(4 级层序) 划分为 6 个小层(5 级层序)， IV 油组划分为 7 个小层。小层对应的曲线趋势变化十分明显，提高了高频层序划分的精度。
(3)在 INPEFA 旋回趋势建立的等时框架下，对南华区块井间进行精细的小层对比与砂体连通性对比，结合油水关系进行单砂体的连通性预测，刻画出砂体的空间分布，其结果与该区实际的动态开发资料具有很好的一致性。
