一种多孔介质渗透率的确定方法技术

本发明专利技术公开了一种多孔介质渗透率的确定方法，属于多孔介质技术领域，多孔介质的孔隙率为φ，组成所述多孔介质的玻璃微珠的平均半径为玻璃微珠的粒径分布偏度为γ；则多孔介质样品的渗透率为：该方法计算简单、精度较高。

【技术实现步骤摘要】

本专利技术属于多孔介质
，具体涉及一种多孔介质渗透率的确定方法。
技术介绍
多孔介质广泛存在于生产生活中，渗透率是描述多孔介质内流体流动能力的关键指标，对渗透率的准确测量和预测是多孔介质研究中的重要内容，不同类型的多孔介质通常有不同的内部微观结构，进而导致了不同多孔介质的渗透率确定方法之间有很大差异。通常，对于多孔介质渗透率的确定方法，除了直接测量外还可以测量多孔介质部分结构参数，再通过计算得到。
技术实现思路
有鉴于此，本专利技术提供了一种多孔介质渗透率的确定方法，该方法计算简单、精度较高。为了达到上述目的，本专利技术的技术方案为：一种多孔介质渗透率的确定方法，多孔介质的孔隙率为φ，组成所述多孔介质的玻璃微珠的平均半径为玻璃微珠的粒径分布偏度为γ；则多孔介质样品的渗透率为：K=φ3d‾p227(4.762φ-4.270γ3+16.950γ2-22.310γ+10.061)2(1-φ)2.]]>进一步地，多孔介质由玻璃微珠烧结制备而成，玻璃微珠的粒径范围为20μm～200μm；制备方法为：制作圆柱形的金属模具，将玻璃微珠填满模具并放入恒温炉进行加热，将微珠烧结在一起，获得圆柱形的多孔介质，多孔介质的圆形端面直径为为6mm～15mm，厚度为1mm～5mm。进一步地，多孔介质的孔隙率为φ的测量方法为：测量获得多孔介质的总质量M、总体积V以及玻璃微珠的密度ρ，则孔隙率进一步地，测量玻璃微珠的粒径分布参数：包括玻璃微珠的平均半径中位径dpMed，粒径累计分布参数dp(0.25)和dp(0.75)；其中粒径累计分布参数表示小于该粒径的微珠占微珠总体积的体积分数；dp(0.25)为小于该粒径的所有微珠的体积总和占微珠体积总和的25％，dp(0.75)为小于该粒径的所有微珠的体积总和占微珠体积总和的75％。所述多孔介质中玻璃微珠的分选系数为：当玻璃微珠的粒径分布满足对数正态分布时，其均值为和标准差分别为则所述玻璃微珠粒径分布偏度γ计算方法为其中Ei(i＝1,2,3)为玻璃微珠粒径的i阶原点矩。有益效果：本专利技术利用由玻璃微珠制成的多孔介质的孔隙率及其玻璃微珠粒径分布参数实现该多孔介质渗透率的确定，因此这种渗透率的确定方式计算简单，将该方法的渗透率计算结果与实验测量的渗透率结果进行比对，发现该方法的渗透率计算结果精度较高。对于微珠粒径分布参数已知的样品，只需称量其质量和体积即可确定渗透率，无需使用专用的渗透率测量设备进行实验测量。附图说明图1有玻璃微珠制成多孔介质的微观照片图2星形孔隙；图3迂曲度与孔隙率的关系图；图4迂曲度与偏度的关系图。具体实施方式下面结合附图并举实施例，对本专利技术进行详细描述。实施例1，本方法步骤如下：步骤一：多孔介质样品制备：该种多孔介质由大量玻璃微珠烧结而成，玻璃微珠的粒径范围为20μm～200μm。制作圆柱形的金属模具，将微珠填满模具并放入恒温炉进行加热，将微珠烧结在一起，获得圆柱形的多孔介质样品。样品圆形端面直径通常为6mm～15mm，厚度为1mm～5mm。步骤二：多孔介质样品孔隙率φ的测量：测量样品的总质量M、总体积V以及玻璃微珠的密度ρ，则样品孔隙率为：φ=1-MρV---(1)]]>步骤三：玻璃微珠粒径分布参数测量：利用激光粒度分析仪对玻璃微珠的粒径分布进行测量。测量玻璃微珠的平均半径中位径dpMed，粒径累计分布参数dp(0.25)和dp(0.75)。其中粒径累计分布参数表示小于该粒径的微珠占微珠总体积的体积分数。即dp(0.25)为小于该粒径的所有微珠的体积总和占微珠体积总和的25％，dp(0.75)为小于该粒径的所有微珠的体积总和占微珠体积总和的75％。步骤四：玻璃微珠粒径分布参数计算该多孔介质中玻璃微珠的分选系数为：So=log2(dp(0.25))1og2(dp(0.75))---(2)]]>认为玻璃微珠的粒径分布满足对数正态分布，其均值和标准差分别为：其中dpMed单位为mm。再通过下列计算公式得到玻璃微珠粒径分布偏度γ，其中Ei(i＝1,2,3)为玻璃微珠粒径的i阶原点矩。γ=E3-3E1E2+2E13(E2-E12)32---(8)]]>步骤五：由式(9)即可确定该多孔介质样品的渗透率K=φ3d‾p227(4.762φ-4.270γ3+16.950γ2-22.310γ+10.061)2(1-φ)2---(9)]]>式(9)的原理：对于多孔介质渗透率的计算，有经典算法Kozeny-Carman(KC)方程：K=φ336k(1-φ)2d‾p2---(10)]]>其中k＝k0τ2称为KC常数，k0为形状因子，与多孔介质中孔隙的形状有关，τ为多孔介质中的迂曲度，与多孔介质中流道的弯曲程度有关。KC方程是一个半经验方程，通常将多孔介质中的孔隙直接看成圆形，此时k0取为2，并将KC常数直接取为5。但是对于该种玻璃微珠制成的多孔介质，该KC常数的选择并不准确，本专利测试了5种不同粒径分布共18片多孔介质的渗透率，根据实测数据，利用最小二乘法计算式(10)中的KC常数，应为3.5686，与5相差较大。同时利用该模型计算该种由微珠构成的多孔介质时误差较大，KC常数取3.5686时，相对误差((理论值-实测值)/实测值*100％)最小，约为36.8％。因此分别讨论形状因子k0和迂曲度τ的取值，以提高计算精度。由玻璃微珠构成的多孔介质，微珠的球形特征保存的较好，如图1所示。因此，微珠之间的孔隙形状应该看成星形，如图2所示，而不是圆形。利用Ansys-Fluent仿真环境计算星形孔隙的形状因子k0，得到：k0＝0.75(11)对于迂曲度τ的取值，本专利认为其应该与多孔介质的孔隙率φ和微珠粒径分布偏度γ这两个因素有关。利用式(12)计算迂曲度，并分别分析迂曲度和孔隙率φ以及微珠粒径分布偏度γ之间的关系，实测实验结果如图3，图4所示。τ=φ336k0K(1-φ)2d‾p2---(12)]]>由实验数据通过最小二乘法拟合分别得到迂曲度和孔隙率以及微珠粒径分布偏度的定量关系，如式(13)，(14)。τ1(φ)＝9.523φ-1.029(13)τ2(γ)＝8.539γ3+33.9γ2-44.62γ+21.15(14)...

【技术保护点】
一种多孔介质渗透率的确定方法，其特征在于，所述多孔介质的孔隙率为φ，组成所述多孔介质的玻璃微珠的平均半径为玻璃微珠的粒径分布偏度为γ；则多孔介质样品的渗透率为：K=φ3d‾p227(4.762φ-4.270γ3+16.950γ2-22.310γ+10.061)2(1-φ)2.]]>

【技术特征摘要】
1.一种多孔介质渗透率的确定方法，其特征在于，所述多孔介质的孔隙率
为φ，组成所述多孔介质的玻璃微珠的平均半径为玻璃微珠的粒径分布偏度
为γ；则多孔介质样品的渗透率为：
K=φ3d‾p227(4.762φ-4.270γ3+16.950γ2-22.310γ+10.061)2(1-φ)2.]]>2.如权利要求1所述的一种多孔介质渗透率的确定方法，其特征在于，所
述多孔介质由玻璃微珠烧结制备而成，玻璃微珠的粒径范围为20μm～200μm；
制备方法为：制作圆柱形的金属模具，将玻璃微珠填满模具并放入恒温炉进行
加热，将微珠烧结在一起，获得圆柱形的多孔介质，多孔介质的圆形端面直径
为为6mm～15mm，厚度为1mm～5mm。
3.如权利要求1所述的一种多孔介质渗透率...

【专利技术属性】
技术研发人员：付梦印，程思源，王美玲，明丽，李响，肖梅峰，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：北京;11