论文部分内容阅读
[摘 要]研究工区受沉积环境、构造因素及成岩作用的综合影响下，储层非均质性强、层内与层间变化大、综合含水差异大、流动层带复杂等问题突出，进行储层流动单元研究显得尤为重要。本文通过主因子分析从能够反映沉积相、储层岩石物理特征、成岩作用以及岩石的微观孔隙结构等15个参数中优选出符合研究区实际的孔隙度、渗透率、动层指数、渗透率极比、渗透率均质系数、砂体厚度、隔夹层分布密度等7个参数作为流动单元划分依据，建立判别函数。将研究区流动单元划为3类，结果显示I类流动单元具有良好的综合特性，与沉积微相展布特征相一致。
[关键词]流动单元；非均质性；因子分析；聚类分析
中图分类号：TH862 文献标识码：A 文章编号：1009-914X（2015）03-0354-01
1 研究区地质概况
研究区虎狼峁位于鄂尔多斯盆地一级构造单元陕北斜坡中东部，虎狼峁区长6储层主要发育中生界三叠统延长组地层。结合传统旋回厚度法及地层旋回的等时对比法，在对比标志层控制下，将长6油组划为3个砂组，6个小层。储层岩性以细砂为主，属于三角洲前缘亚相沉积。经岩性、电性特征统计，长6物性较差，平均渗透率2.7×10-3um2，平均孔隙度11.1%；无论是平面，还是纵向上，非均质性都是非常严重。
2 流动单元划分方法
2.1 定性选择流动单元参数基础
流动单元的划分应充分考虑影响的多种因素，主要包括：沉积相、储层岩石物理特征、成岩作用以及岩石的微观孔隙结构等。针对研究区特质选用了孔隙度（？）、渗透率（K）、动层指数（FZI）、砂层厚度（H）、有效厚度（h）、储层品质系数[（？/K ）1/2]、地层系数（KH）、储油能力（？*So*H）、含油饱和度（So）、渗透率变异系数（Vk）、渗透率突进系数（Tk）、渗透率极差（Jk）、渗透率极比（Bk）、渗透率均值系数（kp）、隔夹层分布密度（Dk）共15个参数做为划分基础数据。
2.2 参数优化
2.2.1 因子分析对流动单元降维
因子分析[1]是多参数分析的基础，利用斜方差逼近技术，将n维数据归纳主成分达到降维之目的， 再通过主成分坐标旋转，将变量之间的相关性表示更为清晰明了。本次统计了研究区461口井中15个类型各6241个流动单元的参数，在小层内各参数的均值构成1597*15个样品。通过主因子—R型分析，形成了4个主因子和两副因子载荷图。为了进一步验证所选参数的可靠性， 又对15项参数用逐步判别分析方法进行了挑选， 不断改变引入和剔出参数的临界值 F 1、F 2 ，挑选出孔隙度、渗透率、动层指数、渗透率极比、渗透率均质系数、砂体厚度、隔夹层分布密度等7個参数[2-4]，以此作为研究区提取的实际地质参数[5]。
2.2.2 流动单元聚类分析
聚类分析是根据研究对象之间亲疏关系。本次研究从储层沉积微相、岩性、物性和孔隙度、渗透率特征等反复进行对比分析，逐层对比解释是否合理， 如不合理则需选择不同的分类统计量进行再次分类， 直至每个小层样品对比分类与沉积、岩性、物性、渗流特性吻合，才能算完成聚类分析。
本次提出的方法如下：假设我们想划分为3种简单类型，即： I类、II类和III类。I类表示好的流动单元；II类表示中等的流动单元； III类表示差的流动单元。
2.3 建立流动单元判别函数
在已知分类的基础上，利用逐步判别分析方法建立基于贝叶斯准则下的多组判别函数进而使之推广到那些未得出分类的区域。通过调整引入7个参数，引入顺序为渗透率、动层指数、渗透率极比、渗透率均质系数、砂体厚度、孔隙度、隔夹层分布密度得到如下公式：
F1=-252.761+134.611FZI-77.802（？/K）1/2+1.565H+29.409Dk+3372.389？-18.562kp-0.105Bk
F2=-237.109+132.513FZI-75.465（？/K）1/2+1.596H+26.495Dk+3315.719？-23.117kp-0.080Bk
F3=-223.109+120.237FZI-77.003（？/K）1/2+1.778H+23.462Dk+3228.674？-24.753kp-0.052Bk
3 应用效果
流动单元体现的是一种整体的效果，必然会得出相互叠置的参数范围，但体现的却是多种因素的交互影响。其实际规律与沉积相态有强相关性，因沉积相也是多种地质因素共同作用的结果，所以其实际规律与沉积相态有强相关性。
由表1可以看出：I类流动单元具有良好的综合特性，孔隙几何特征和渗流、储集能力。渗透率均值3.18×10-3μm2，孔隙度均值为11.95%分布在14%—10%之间，动层指数均值为1.15；主要发育水下分流河道，叠置的河口砂坝沉积。此类孔隙分布均匀且连通性好。它们控制了研究区主要的含油有利区范围，水淹程度强。II类流动单元有较好的综合特性，渗透率均值2.00×10-3μm2，孔隙度均值为11.20%在 13.8% -9.4%之间，动层指数均值为0.972 之间；主要发育布在河口坝边部、水下分流河道侧缘；存在一定的剩余油。III类流动单元综合特性较差，渗透率均值1.39×10-3μm2，孔隙度均值为10.43%，动层指数均值为0.82，主要发育水下分流间湾水下分流河道侧缘和前缘席状砂。从层面上来看I类流动单元主要分布在长611、612之上。综合因子分析结果与基本沉积相研究一致。
4 结语
（1）本文通过因子分析和逐步判别分析对参数进行优选，挑选出孔隙度、渗透率、动层指数、渗透率极比、渗透率均质系数、砂体厚度、隔夹层分布密度7个参数作为研究区提取的实际地质参数。利用这7个参数，建立了3个流动单元划分函数，并将流动单元划分为I类、II类和III类。
（2）通过应用效果分析认为，I类流动单元具有良好的综合特性，孔隙几何特征和渗流、储集能力，是主要的含油有利区范围，水淹程度强。II类流动单元有较好的综合特性存在一定的剩余油。III类流动单元综合特性较差。
[关键词]流动单元；非均质性；因子分析；聚类分析
中图分类号：TH862 文献标识码：A 文章编号：1009-914X（2015）03-0354-01
1 研究区地质概况
研究区虎狼峁位于鄂尔多斯盆地一级构造单元陕北斜坡中东部，虎狼峁区长6储层主要发育中生界三叠统延长组地层。结合传统旋回厚度法及地层旋回的等时对比法，在对比标志层控制下，将长6油组划为3个砂组，6个小层。储层岩性以细砂为主，属于三角洲前缘亚相沉积。经岩性、电性特征统计，长6物性较差，平均渗透率2.7×10-3um2，平均孔隙度11.1%；无论是平面，还是纵向上，非均质性都是非常严重。
2 流动单元划分方法
2.1 定性选择流动单元参数基础
流动单元的划分应充分考虑影响的多种因素，主要包括：沉积相、储层岩石物理特征、成岩作用以及岩石的微观孔隙结构等。针对研究区特质选用了孔隙度（？）、渗透率（K）、动层指数（FZI）、砂层厚度（H）、有效厚度（h）、储层品质系数[（？/K ）1/2]、地层系数（KH）、储油能力（？*So*H）、含油饱和度（So）、渗透率变异系数（Vk）、渗透率突进系数（Tk）、渗透率极差（Jk）、渗透率极比（Bk）、渗透率均值系数（kp）、隔夹层分布密度（Dk）共15个参数做为划分基础数据。
2.2 参数优化
2.2.1 因子分析对流动单元降维
因子分析[1]是多参数分析的基础，利用斜方差逼近技术，将n维数据归纳主成分达到降维之目的， 再通过主成分坐标旋转，将变量之间的相关性表示更为清晰明了。本次统计了研究区461口井中15个类型各6241个流动单元的参数，在小层内各参数的均值构成1597*15个样品。通过主因子—R型分析，形成了4个主因子和两副因子载荷图。为了进一步验证所选参数的可靠性， 又对15项参数用逐步判别分析方法进行了挑选， 不断改变引入和剔出参数的临界值 F 1、F 2 ，挑选出孔隙度、渗透率、动层指数、渗透率极比、渗透率均质系数、砂体厚度、隔夹层分布密度等7個参数[2-4]，以此作为研究区提取的实际地质参数[5]。
2.2.2 流动单元聚类分析
聚类分析是根据研究对象之间亲疏关系。本次研究从储层沉积微相、岩性、物性和孔隙度、渗透率特征等反复进行对比分析，逐层对比解释是否合理， 如不合理则需选择不同的分类统计量进行再次分类， 直至每个小层样品对比分类与沉积、岩性、物性、渗流特性吻合，才能算完成聚类分析。
本次提出的方法如下：假设我们想划分为3种简单类型，即： I类、II类和III类。I类表示好的流动单元；II类表示中等的流动单元； III类表示差的流动单元。
2.3 建立流动单元判别函数
在已知分类的基础上，利用逐步判别分析方法建立基于贝叶斯准则下的多组判别函数进而使之推广到那些未得出分类的区域。通过调整引入7个参数，引入顺序为渗透率、动层指数、渗透率极比、渗透率均质系数、砂体厚度、孔隙度、隔夹层分布密度得到如下公式：
F1=-252.761+134.611FZI-77.802（？/K）1/2+1.565H+29.409Dk+3372.389？-18.562kp-0.105Bk
F2=-237.109+132.513FZI-75.465（？/K）1/2+1.596H+26.495Dk+3315.719？-23.117kp-0.080Bk
F3=-223.109+120.237FZI-77.003（？/K）1/2+1.778H+23.462Dk+3228.674？-24.753kp-0.052Bk
3 应用效果
流动单元体现的是一种整体的效果，必然会得出相互叠置的参数范围，但体现的却是多种因素的交互影响。其实际规律与沉积相态有强相关性，因沉积相也是多种地质因素共同作用的结果，所以其实际规律与沉积相态有强相关性。
由表1可以看出：I类流动单元具有良好的综合特性，孔隙几何特征和渗流、储集能力。渗透率均值3.18×10-3μm2，孔隙度均值为11.95%分布在14%—10%之间，动层指数均值为1.15；主要发育水下分流河道，叠置的河口砂坝沉积。此类孔隙分布均匀且连通性好。它们控制了研究区主要的含油有利区范围，水淹程度强。II类流动单元有较好的综合特性，渗透率均值2.00×10-3μm2，孔隙度均值为11.20%在 13.8% -9.4%之间，动层指数均值为0.972 之间；主要发育布在河口坝边部、水下分流河道侧缘；存在一定的剩余油。III类流动单元综合特性较差，渗透率均值1.39×10-3μm2，孔隙度均值为10.43%，动层指数均值为0.82，主要发育水下分流间湾水下分流河道侧缘和前缘席状砂。从层面上来看I类流动单元主要分布在长611、612之上。综合因子分析结果与基本沉积相研究一致。
4 结语
（1）本文通过因子分析和逐步判别分析对参数进行优选，挑选出孔隙度、渗透率、动层指数、渗透率极比、渗透率均质系数、砂体厚度、隔夹层分布密度7个参数作为研究区提取的实际地质参数。利用这7个参数，建立了3个流动单元划分函数，并将流动单元划分为I类、II类和III类。
（2）通过应用效果分析认为，I类流动单元具有良好的综合特性，孔隙几何特征和渗流、储集能力，是主要的含油有利区范围，水淹程度强。II类流动单元有较好的综合特性存在一定的剩余油。III类流动单元综合特性较差。