论文部分内容阅读
【摘要】根据Kozeny-Caman 方程，建立了不同流动单元渗透率的表达形式。采用聚类分析方法，选取最大隶属度大于70%的最优样本进行建模，得到流动单元的判别模型； 建立了复杂储层渗透率模型，提高了储层流体渗透性能的预测精度。
【关键词】渗透率 流动单元 Kozeny-Caman方程 聚类分析 判别模型（式6）
由（6）式可知，在RQI与Φz的双对数坐标图上，所有具有相似FZI值的样品点将形成一条斜率为1的直线。FZI的数值可由该斜率直线与Φz=1相交的截距确定。一般认为，属于同一直线的样品点（具有相似FZI值的样品点）具有相似的孔喉性质，属于同一类流动单元。
2 流动单元法与聚类分析结合的渗透率模型的建立
为了验证以上公式的精度，把流动单元计算的渗透率和该工区广泛应用的经验公式计算的渗透率进行对比，可以看出流动单元划分得出的渗透率精度高，符合精细评价的标准。
2.1 应用聚类判别分析对未取心井段渗透率的评价
流动层带指数是表征流动单元的关键参数之一，此外，各种测井方法，如自然伽马、声波时差等都与流动单元的划分有关。基于模糊聚类的思想，将最能反映岩石岩性和物性的自然伽马（ GR ）、声波时差（AC ）、补偿密度（ DEN ）、补偿中子（CNL）等测井数据并结合FZI 参数进行模糊聚类分析，根据最大隶属度原则得到其聚类结果。
由模糊聚类分析结果可以看出，采用模糊聚类算法得出的4种储层类别与图1所划分的4种通过聚类分析得出的类别和我们所划分的流动单元对应性良好，说明这种模型是我们划分流动单元的另一种表现形式。为了对未取芯井进行渗透率计算，选择最大隶属度大于70%的岩芯数据建立标准样本，选取GR、AC、DEN、CNL 等4 种测井曲线进行主成分压缩，建立判别函数，进行流动单元划分，计算不同流动单元的渗透率。
采用FISHER线性判别方法，可以得到4类不同流动单元的判别公式
F1=3.265GR+5.567AC+823.692DEN-4 3 3 . 9 0 1 C N L - 1 7 0 4 . 9 2 5（式1 1）F2=3.768GR+5.672AC+813.0912DEN-4 3 3 . 2 6 9 C N L - 1 7 6 0 . 5 6 8 （式1 2）F2=3.692GR+4.998AC+843.274DEN-446.299CNL-1615.439 （式13）F4=4.584G R+5.829A C+820.432E N-420.80CNL-1859.966 （式14）
通过对样品回带判别，其符合率达到94.8%，说明建立模型的可行性。通过把实测的测井值代入判别公式，其最大值就是流动单元的类别。再利用（7）-（10）可以算出来不同流动单元的渗透率。
3 结论
（1）渗透率和孔隙度之间的关系不是简单的指数关系，它受储层孔隙结构以及流动单元的影响，不同的流动单元可以建立不同的解释模型。
（2）采用实测的测井数据来计算渗透率，克服了未取心井段FZI难以划分的缺点。
（3）根据最大隶属原则和聚类判别分析可以提高渗透率计算的精度。
参考文献
[1] Ebanks WJ . Flow unit concept-integrated approach to reservoir description for engineering projects[J ]. AAPGBulletin，1987，71（5）：551O552
[2] 欧家强，罗明高，王小蓉. 低渗透油藏中储层流动单元划分[J]. 西安石油大学学报： 自然科学版，2008，23（ 5）： 24- 28
[3] 袁新涛，彭仕宓，林承焰，贾瑞高，分流动单元精确求取储层渗透率的方法[J].石油学报，2005，26（6）：21-23