率；
[0054] 步骤105 :计算每个层段中深度点的平均渗透率级差；
[0055] 步骤106 :根据每个层段中深度点的最大渗透率、平均渗透率和平均渗透率级差， 结合渗流单元的分级标准，判断每个层段所属的渗流单元级别。
[0056] 本发明以新疆准噶尔盆地西北缘下克拉玛依组砾岩油藏为例，对本发明方法进行 详细的描述。
[0057] 具体实施时，步骤101为基于测井资料获得待研究地层的多个深度点的渗透率。 具体的，步骤101就是准备数据的过程。所准备的数据还可以包括井名、地层层位和地层孔 隙度。根据表1的形式准备基础数据。以T6014井为例，数据表的内容包括井名、层位（单 砂层）、深度（点）、孔隙度P0R和渗透率PERM。用i表示深度点的个数，i取整数，i大于 等于1，每0. 125m采集一个深度点，即相邻两个深度点的深度值相差为0. 125m。其中泥岩 的孔隙度、渗透率用-999. 25表示，渗透率的单位为毫达西，md。
[0058] 表 1
[0059]
[0060] 具体实施时，步骤102为将待研究地层的多个深度点的渗透率进行分级。具体的， 对渗透率进行分级的标准为：当第i深度点的渗透率PERMi大于等于2000md时，所述第i 深度点的渗透率为1级；当第i深度点的渗透率大于等于500md小于2000md时，所述第i 深度点的渗透率为2级；当第i深度点的渗透率大于等于50md小于500md时，所述第i深 度点的渗透率为3级；当第i深度点的渗透率小于50md时，所述第i深度点的渗透率为4 级；当第i深度点的渗透率等于-999. 25md时，所述第i深度点的渗透率为5级。具体如表 2所示。
[0061] 表 2
[0062]
[0063] 在某一 口井中，在对各个深度点的渗透率进行分级的过程中，可能会出现异常深 度点，所说的异常深度点指的是：若第i深度点Di为泥岩（即渗透率为-999. 25md)，第i+2 深度点Di+2为泥岩，而第i+1深度点Di+1为非泥岩（即渗透率不等于-999. 25md)，那么 Di+1就为异常深度点，如表3中数据点编号为5的行数据就为异常深度点。
[0064]表3
[0065]
[0066] 若出现异常深度点，需将Di+1数据点的岩相修正为泥岩层，即将渗透率的等级修 改正为5。修正后的渗透率分级表如表4所示。
[0067]表 4
[0068]
[0069] 具体实施时，步骤103为将待研究地层的多个深度点按深度值从小到大的顺序进 行排序，形成深度点序列；将深度点序列中深度值相邻且具有相同等级的渗透率的深度点 划分至同一层段中。具体的，一个层段中所包含的深度点的渗透率的等级是相同的，且所包 含的深度点的深度值是相邻的。在分完层段之后，在每一层段的第一个深度点记号1，如表 5所示。
[0070] 表 5
[0071]
[0072] 具体实施时，步骤104为计算每个层段中深度点的平均渗透率。具体的，按照如下 公式计算：
[0073]
[0074] 式（1)中，Kave_j为第j层段中深度点的平均渗透率，j取整数，j大于等于1 ;
[0075] m为第j层段内的相同级别的渗透率的深度点的个数；
[0076] PERMi为第j层段内的第i深度点的渗透率。
[0077] 另外，还可以计算各层段中深度点的厚度H_j，按照如下公式计算：
[0078] H_j= 0. 125Xm (2)
[0079] 式⑵中，H_j为第j层段中深度点的厚度。
[0080] 注：每一层段内泥岩层不参与计算。按照公式（1)和（2)计算所得的数据如表6 所示。
[0081] 表 6
[0082]
[0083] 在具体实施时，若某一层段内所包含的深度点个数小于3的时候，需要将该层段 的渗透率等级修正为与之相邻的层段的渗透率等级。如表6中加粗的部分所示，需要将其 渗透率等级修正为2,在计算该层段中深度点的平均渗透率时，该层段所包含的级别为2的 地层渗透率的深度点按照8个算。
[0084] 具体实施时，步骤105为计算每个层段中深度点的平均渗透率级差。具体的，按照 如下公式计算：
[0085] Kave_dif_j=Kave_j/Kave_min_p(3)
[0086] 式（3)中，Kave_dif_j为第j层段中深度点的平均渗透率级差；
[0087] Kave_min_p为所有层段中深度点的平均渗透率的最小值，其计算公式如下：
[0088] Kave_min_p=Minimum(Kave_j) (4)
[0089] 另外，还可以计算所有层段中深度点的平均渗透率的最大值Kave_maX_p，按照如 下公式计算：
[0090] Kave_max_p=Maximum(Kave_j) (5)
[0091] 检索表6中数据，统计第j层段中深度点的最小渗透率值Kmin_j，第j层段中深度 点的最大渗透率值Kmax_j，各层段中深度点的平均渗透率值Kave_j，按照公式（3)、（4)和 (5)计算之后，获得的数据如表7所示。
[0092] 表 7
[0093]
[0095] 具体实施时，步骤106为根据每个层段中深度点的最大渗透率、平均渗透率和平 均渗透率级差，结合渗流单元的分级标准，判断每个层段所属的渗流单元级别。
[0096] 具体的，首先详细介绍渗流单元的分级标准。以新疆克拉玛依油田砾岩储层为例， 结合地质成因分析，将示踪剂采出曲线按照形态划分为非对称尖峰型、对称尖峰型、箱型和 双峰型等4类。图2为非对称尖峰型曲线示意图，其反映注入水推进快、高渗层薄、等效渗 透率大、渗透率级差大，为I级渗流单元，图2中的纵坐标表示示踪剂含量，"7202"表示井 号，"Er"表示示踪剂含量。图3为对称尖峰型曲线示意图，其反映油井见剂时间较晚，井间 以高渗透条带为主，为II级渗流单元，图3中的纵坐标表示示踪剂含量，"T72311"表示井 号，"Er"表示示踪剂含量。图4为箱型曲线示意图，其表明油井见剂时间略长，井间以中一 高渗透厚层为主，为III级渗流单元，图4中的纵坐标表示示踪剂含量，"T72516"表示井号， "MT10"表示示踪剂含量。图5为双峰型曲线示意图，其表明井间中低渗储层非均质性强， 中渗储层厚度小、多层，为IV级渗流单元，图5中的纵坐标表示示踪剂含量，"T72315"表示 井号，La表示示踪剂含量。对油水井注采单元进行动态模拟，其结果显示，渗流单元的分布 有3个特点：（1)油水井距离小的井区较发育。（2)顺沉积方向的油水井间I级和II级渗 流单元较集中。（3)高产液井附近水流明显增多。根据上述4级渗流单元所在目的层段的 渗流特征参数统计结果表明，I级渗流单元最大渗透率大于2000md，平均渗透率大于50md， 渗透率级差大于17 ;II级渗流单元最大渗透率介于500md至2000md之间，平均渗透率介 于50md至2000md之间，渗透率级差同样大于17 ;111级渗流单元最大渗透率介于500md至 2000md之间，平均渗透率介于50md至2000md之间，渗透率级差介于8至17之间；IV级渗 流单元最大渗透率介于50md至500md之间，平均渗透率介于50md至500md之间，渗透率级 差大于17。另，当渗透率级差小于8,则表明是正常地层。根据上述渗流单元相关参数地质 统计规律，结合储层构型类型、渗流能力、非均质特点、油水井的示踪剂监测、生产动态及油 水井生产模拟结果，建立了该区各级渗流单元划分标准，如表8所示。
[0097] 表 8