4. 孔渗饱和度参数统计与分布特征分析
孔渗饱和度分布图与储层对比的核心,不只是把单井参数画到平面上,更在于先用统计手段弄清每个参数自身的分布形态、离散程度和异常个体,再据此判断平面图与对比剖面是否可信。本节围绕直方图、累积频率曲线、交会图和分区分层统计四类工具展开,目标是回答三个问题:参数总体服从什么分布、哪些井点偏离主体、不同区块或层段的参数是否存在系统性差异。
一、直方图:先看形态,再看尾部
直方图用于观察单参数的概率密度形态。以孔隙度为例,把岩心分析或测井解释孔隙度按 1% 或 2% 的组距分组,统计各组频数并归一化为频率,即可得到分布形态。实际操作中要注意三点:一是组距选择,组距过粗会掩盖双峰,过细则频数抖动大,一般以参数标准差的 1/4~1/3 作为组距初值;二是样本代表性,岩心分析样点常偏向物性好的层段,需与测井解释样本对照;三是分布形态判读,单峰近对称说明参数较均质,双峰往往指示存在两套不同物性的储层或不同沉积微相,长尾则提示存在低值或高值异常群体。
示例性数据(合理假设):某区块 32 口井的测井解释孔隙度,均值 18.6%,标准差 3.4%,直方图呈单峰略左偏,主峰位于 17%~20%,在 10%~12% 区间出现一个孤立小峰,对应 3 口井。这 3 口井需回到测井曲线与岩心描述核查,判断是真实致密段还是解释偏差。
二、累积频率曲线:读分位数与截断值
累积频率曲线由直方图逐组累加得到,横轴为参数值,纵轴为小于某值的累积频率。它的实用价值在于快速读出分位数:P10、P50、P90 分别对应累积频率 10%、50%、90% 处的参数值,可用于确定储层参数分级界限和储量计算中的参数取值。若曲线在中段近似直线,说明主体分布接近正态;若出现明显台阶,说明某区间样本集中,往往是分级界限或截断值所在。
仍以上述示例数据为例,读得孔隙度 P90≈14.2%、P50≈18.5%、P10≈22.8%。可将 14% 作为有效储层孔隙度下限的参考值,将 22% 以上划为相对优质储层。需要强调,分位数只是统计参考,下限最终要结合岩心物性、试油和产能资料标定,不能仅凭曲线截取。
三、交会图:看参数间相关性与异常点
交会图把两个参数作为横纵坐标投点,用于判断相关性、识别异常点和区分不同群体。储层分析中最常用的是孔隙度—渗透率交会图(半对数坐标,渗透率取对数轴),其次为孔隙度—饱和度、渗透率—饱和度交会图。判读时关注三点:一是趋势线的斜率与相关系数,孔渗在半对数坐标下常呈线性关系,斜率反映孔渗结构效率;二是散点带的宽度,带宽大说明孔渗关系受裂缝、微相或黏土含量影响明显;三是离群点,偏离主体趋势带的点往往是解释错误、裂缝段或致密夹层,应逐点回溯。
示例性做法:将 32 口井的孔渗数据投到半对数交会图上,主体趋势带相关系数约 0.82,但有 4 个点明显位于趋势带右下方,即孔隙度不低而渗透率偏低。这类点通常对应孔隙型但喉道细的储层,若平面图上恰位于某微相边界,则提示该区存在储层质量变差的可能,需在储量计算中单独处理。
四、分区分层统计:把总体拆开看
总体统计容易掩盖系统性差异,因此要按区块、层段、沉积微相分别统计均值、标准差、变异系数和分位数,形成对比表。判断差异是否显著,可比较均值差与标准差的关系,也可用变异系数衡量离散程度。若某层段变异系数明显高于其他层段,说明该层段非均质性强,平面图插值时应收紧搜索半径、增加控制点。
| 统计单元 | 样本数 | 孔隙度均值(%) | 渗透率均值(mD) | 变异系数 | 判读 |
|---|---|---|---|---|---|
| 全区 | 32 | 18.6 | 86 | 0.18 | 总体较均质 |
| 北区 | 14 | 20.1 | 132 | 0.15 | 物性好、离散小 |
| 南区 | 18 | 17.4 | 51 | 0.24 | 物性偏差、非均质强 |
| Ⅰ层段 | 12 | 19.8 | 118 | 0.16 | 主力层,相对均质 |
| Ⅱ层段 | 20 | 17.9 | 67 | 0.23 | 非均质强,需加密控制 |
表中数据为教学示例的合理假设。可以看到,南区与Ⅱ层段的变异系数均超过 0.2,是后续平面图与对比剖面中需要重点关注的区域和层段;北区与Ⅰ层段离散小,插值结果相对可靠。
五、异常值识别与处理原则
异常值不等于错误值,处理前必须先定性。常用识别方法包括:箱线图法(超出 Q1−1.5IQR 或 Q3+1.5IQR 视为可疑)、Z 分数法(|Z|>3 视为可疑)、交会图离群点法。识别出可疑点后,按“先核查、再分类、后处理”的顺序推进:核查测井曲线质量、岩心描述与解释参数;分类为真实低值段、裂缝高值段或解释错误;真实值保留并在统计中单独说明,解释错误则修正或剔除,剔除比例一般不宜超过样本的 5%,且须记录剔除依据。