1. 引言
储层评价工作贯穿于油气勘探开发的全过程中。评价储层的过程，其实就是对储层深化认识的过程。而对储层进行评价，首先要对储层进行合理的分类。目前，国内外研究储层的学者应用不同分析方法对储层评价做出了比较详细的探讨和研究，也给出了一些定量化研究结果，比如刘斌在对朝阳沟油田进行储层评价时应用了层次分析法(Analytic Hierarchy Process) [1] ，高鹏宇等在研究河流相储层时应用了主成分分析法(Principal Component Analysis) [2] ，唐俊等在研究鄂尔多斯盆地的姬塬地区长81储层时应用了Q型聚类分析法(Q Cluster Analysis)和判别分析法(Discriminant analysis) [3] 。但是，针对某一具体研究领域，选用哪些参数、运用哪种评价方法还是一个很难解决的问题。通过对本地区储层的研究和对各种评价方法的对比，本文选用Q型聚类分析方法，对储层进行分类及评价，经分析得出了储层的平面分布特征。
2. 聚类分析
2.1. 聚类分析基本原理
聚类分析(Cluster Analysis)是一种数学方法，是根据研究对象的特征对研究对象进行分类的一种多元统计分析技术 [3] 。聚类分析大致上可分为两类：对样本进行分类处理，称为Q型聚类分析；对变量进行分类处理，则称为R型聚类分析。本文中以应用Q型聚类分析为主。
首先，假设我们所研究的样本之间存在着程度不同的相似性，根据样本可找出一些能够度量样本之间相似程度的统计变量，以此作为划分类别的依据。然后，把一些相似程度较大的样本聚合为一类，把另外一些相似程度较大的聚合为另一类，直到把所有的都聚合完毕，把不同类型一一划出来，形成由小到大的分类系统 [3] 。这个过程，其实就是一个逐步减少类别的过程。最后，再把整个分类系统画成一张谱系图，用这张谱系图就可以把所有样本之间的亲疏关系都表示出来。
2.2. 相似性度量的确定
聚类分析中，样本之间的相似性(亲疏关系)通常用距离描述，而变量之间的相似性(亲疏关系)通常用相似系数来描述。先假设每个样品可由p个变量描述其性质，那么每个样品都可以看成p维空间中的一个点，则n个样品就是p维空间中的n个点 [3] [4] 。将第i个样品xi与第j个样品xj之间的距离记为dij，
把第i个样品xi与第j个样品xj之间的夹角余弦cosθij称为相似系数。则欧氏距离
，相似系数
。
3. 研究区地质概况
3.1. 研究区地理位置
青东5块的地理位置位于山东省东营市莱州湾西部极浅海海域，地势较平坦，平均水深约3.6 m。其构造位置位于济阳坳陷东部，郯庐断裂带西侧青东凹陷西斜坡，西北紧临青坨子凸起，西南为青南洼陷，东临青东凹陷，其面积约20 km2 [5] 。青东5块整体上为一个被断裂系统复杂化的西高东低的斜坡，区内断裂系统以近东西向为主，兼有少量南北向 [5] 。青东5块的构造位置见图1。
3.2. 储层发育特征
青东5块地层自下而上依次为古近系孔店组、沙河街组、新近系馆陶组、明化镇组及第四系平原组，其中，主要的含油气层为沙四段上亚段，地层厚度为140~550 m [5] [6] [7] [8] 。
青东5块储层的主要类型为灰白色粉砂岩以及含砾粗砂岩，储层粒度较粗。岩石的成分成熟度较低，矿物颗粒磨圆度较差，多为次棱角状，分选多为中、中—差，多为颗粒式支撑，点式接触，胶结类型为孔隙式胶结，岩石的结构成熟度也较低 [9] 。储层岩石的矿物成分成熟度及结构成熟度均较低，反映研究区距近物源区较近，属于快速搬运、快速沉积 [10] 。其粒度累积曲线以三段式为主，粒度分布范围宽，累积曲线表现平缓，反映研究区储层分选差 [9] [10] 。其岩性组合以砂泥互层为特征，砂体厚度较薄。垂向上，总体表现为向上粒度变粗的反旋回、进积沉积序列，即表现为扇三角洲沉积特点 [11] [12] 。
Figure 1. Location of Qingdong 5 Block
图1. 青东5块构造位置图
4. 应用
4.1. 储层类型划分标准
利用Q型聚类分析方法进行储层分类评价时，为尽可能合理、准确地划分储层类型，分类样品和储层特征参数的选取是很重要的。所选用的样品应尽可能地包含各种岩性，并且有各种分析资料相对应，这样取得的样品才具有典型性、代表性，有助于储层的完整分类。
不同类型的储层有不同的特征参数值，各特征参数值的大小综合体现了该类储层的内在特征。在测井解释处理过程中总共形成了多种参数，为了选取能够反映本区储层类型的参数，对每个参数进行综合分析，通过相应分析方法与地质经验相结合进行参数选取 [4] [13] ，最终选用孔隙度、渗透率、含油饱和度3个参数作为储层评价的特征参数。
(1) 孔隙度、渗透率是储层研究中的经典参数，反映了储层的总体特征，是表征储层质量优劣的重要参数 [13] 。
(2) 含油饱和度是指油在储层孔隙中的含量分别占总孔隙体积的百分数，反映了油气藏中的含油程度。
综上所述，经过细致分析，筛选出了20口井中的样品，各样品由3个参数定量表征，从而建立了青东5块储层分类的样品集，如表1所示。
Table 1. The reservoir classification standard samples of Qingdong 5 Block
表1. 青东5块井储层分类标准样品集
4.2. 储层分类及结果
储层分类的实质，即是把储层性质相类似的样品划分在同一类中。综合运用Q型聚类分析方法，对所选的20口井的样品进行了分类研究，将这20个样品分成了3种类型的储层，即第1种类型、第2种类型和第3种类型的储层。分类的结果见图2。
结合谱系图，对划分得到的3种类型储层进行分析总结 [3] ：第1种类型储层物性最好，孔隙度均值为34.5%，渗透率均值为2646 × 10−3 μm2，含油饱和度均值为53.8%；第2种类型储层物性较好，孔隙度均值为31.9%，渗透率均值为1650 × 10−3 μm2，含油饱和度均值为43.7%；第3种类型储层物性较差，孔隙度均值为29.6%，渗透率均值为967 × 10−3 μm2，含油饱和度均值为43.1%。参数特征分析表明，由第1种类型到第3种类型的储层，孔隙度、渗透率和含油饱和度基本依次减小，而渗透率的变化幅度较大，基本符合该区储层非均质性较强的地质特征，不同井的储层的质量有所不同。
4.3. 储层平面展布特征
以青东5块Es4上152小层为例，研究该区域储层的平面展布特征。如图3所示。
结合该层沉积微相的平面分布特征进行对比分析 [13] [14] ，可以发现：平面上，青东5块Es4上152小层以第1种类型储层为主，第2种类型储层次之，而第3种类型储层仅分布于较局限的范围内。储集性能较好的第1、2种类型储层主要发育于水下分流河道微相带，主要是沿着三条主河道方向分布；另外，在河口沙坝沉积以及滨浅湖沙坝沉积中也发育有少量的储层，其储层物性较好。第3种类型储层主要发育于水下分流河道的远缘和滨浅湖泥岩带，该区储层的岩性变细，相应地，储层孔渗性降低，储层品质变差。
Figure 2. The dendrogram of cluster analysis to reservoirs in Qingdong 5 Block
图2. 青东5块储层聚类分析谱系图
Figure 3. The classification of Es4 152 reservoirs of Qingdong 5 Block
图3. 青东5块Es4上152储层分类
5. 结论
(1) 青东5块储层具有高孔、中高渗、非均质性强的地质特征，在研究区域上选取孔隙度、渗透率、含油饱和度3个参数作为储层评价的特征参数，利用Q型聚类分析方法对储层进行分类评价，可大致将此储层划分为3种类型。第1种类型储层的物性较好，第2种类型储层次之，第3种类型储层的物性最差。
(2) 该区储层在平面上的分布呈现出规律性变化：第1、2种类型的储层主要分布于水下分流河道微相带及附近，第3种类型储层主要分布于水下分流河道的远缘和滨浅湖泥岩带。
(3) 利用Q型聚类分析方法可以对储层进行定量化的分类以及评价，这在油气藏的勘探开发过程中有非常重要的指示作用。据此得到的评价结果依赖于实际的地质参数，这样可以有效避免由评价者的主观因素所造成的影响，有利于进行更深入的分析研究。