一种多孔介质渗透率的确定方法
【技术领域】
[0001] 本发明属于多孔介质技术领域，具体设及一种多孔介质渗透率的确定方法。
【背景技术】
[0002] 多孔介质广泛存在于生产生活中，渗透率是描述多孔介质内流体流动能力的关键 指标，对渗透率的准确测量和预测是多孔介质研究中的重要内容，不同类型的多孔介质通 常有不同的内部微观结构，进而导致了不同多孔介质的渗透率确定方法之间有很大差异。 通常，对于多孔介质渗透率的确定方法，除了直接测量外还可W测量多孔介质部分结构参 数，再通过计算得到。
【发明内容】
[0003] 有鉴于此，本发明提供了一种多孔介质渗透率的确定方法，该方法计算简单、精度 较高。
[0004] 为了达到上述目的，本发明的技术方案为:一种多孔介质渗透率的确定方法，多孔 介质的孔隙率为Φ，组成所述多孔介质的玻璃微珠的平均半径为^，玻璃微珠的粒径分布 偏度为丫；则多孔介质样品的渗透率为：
[0005]
[0006] 进一步地，多孔介质由玻璃微珠烧结制备而成，玻璃微珠的粒径范围为20WI1~200 μη;制备方法为：制作圆柱形的金属模具，将玻璃微珠填满模具并放入恒溫炉进行加热，将 微珠烧结在一起，获得圆柱形的多孔介质，多孔介质的圆形端面直径为为6mm~15mm，厚度 为1mm~5mm。
[0007] 进一步地，多孔介质的孔隙率为Φ的测量方法为：测量获得多孔介质的总质量Μ、 总体积VW及玻璃微珠的密度Ρ，则孔隙率。
[000引进一步地，测量玻璃微珠的粒径分布参数:包括玻璃微珠的平均半径东，中位径 dpMed，粒径累计分布参数dp(0.25)和dp(0.75);其中粒径累计分布参数表示小于该粒径的微 珠占微珠总体积的体积分数;dp(0.25)为小于该粒径的所有微珠的体积总和占微珠体积总 和的25%，dp(0.75)为小于该粒径的所有微珠的体积总和占微珠体积总和的75%。
[0009]所述多孔介质中玻璃微珠的分选系数为：
[0010]当玻璃微珠的粒径分布满足对数正态分布时，其均值为疗=-10貧乂·?和标准差分 别关
[0011] 则所述玻璃微珠粒径分布偏度丫计算方法为
[0012] 其中Ei (i = 1，2，3)为玻璃微珠粒径的i阶原点矩。
[0016] 有益效果：
[0017] 本发明利用由玻璃微珠制成的多孔介质的孔隙率及其玻璃微珠粒径分布参数实 现该多孔介质渗透率的确定，因此运种渗透率的确定方式计算简单，将该方法的渗透率计 算结果与实验测量的渗透率结果进行比对，发现该方法的渗透率计算结果精度较高。对于 微珠粒径分布参数已知的样品，只需称量其质量和体积即可确定渗透率，无需使用专用的 渗透率测量设备进行实验测量。
【附图说明】
[0018] 图1有玻璃微珠制成多孔介质的微观照片
[0019] 图2星形孔隙；
[0020] 图3迂曲度与孔隙率的关系图；
[0021 ]图4迂曲度与偏度的关系图。
【具体实施方式】
[0022] 下面结合附图并举实施例，对本发明进行详细描述。
[0023] 实施例1，本方法步骤如下：
[0024] 步骤一:多孔介质样品制备:该种多孔介质由大量玻璃微珠烧结而成，玻璃微珠的 粒径范围为20μπι~200μπι。制作圆柱形的金属模具，将微珠填满模具并放入恒溫炉进行加 热，将微珠烧结在一起，获得圆柱形的多孔介质样品。样品圆形端面直径通常为6mm~15mm， 厚度为1mm~5mm。
[0025] 步骤二:多孔介质样品孔隙率Φ的测量:测量样品的总质量Μ、总体积VW及玻璃微 珠的密度Ρ，则样品孔隙率为：
[0026]
(1)
[0027] 步骤玻璃微珠粒径分布参数测量:利用激光粒度分析仪对玻璃微珠的粒径分 布进行测量。测量玻璃微珠的平均半径克，中位径dpMed，粒径累计分布参数dp(0.25)和dp (0.75)。其中粒径累计分布参数表示小于该粒径的微珠占微珠总体积的体积分数。即dp (0.25)为小于该粒径的所有微珠的体积总和占微珠体积总和的25%，dp(0.75)为小于该粒 径的所有微珠的体积总和占微珠体积总和的75%。
[0028] 步骤四：玻璃微珠粒径分布参数计算
[0029] 该多孔介质中玻璃微珠的分选系数为：
[0030]
(2)
[0031 ]认为玻璃微珠的粒径分布满足对数正态分布，其均值和标准差分别为：
[0034]其中dpMed单位为mm。再通过下列计算公式得到玻璃微珠粒径分布偏度丫，其中Ei(i =1，2，3)为玻璃微珠粒径的i阶原点矩。
[0041 ]式(9)的原理：
[0042] 对于多孔介质渗透率的计算，有经典算法Κο Z en厂化man化C)方程：
[0043]
(10)
[0044] 其中k = k〇T2称为KC常数，ko为形状因子，与多孔介质中孔隙的形状有关，τ为多孔 介质中的迂曲度，与多孔介质中流道的弯曲程度有关。KC方程是一个半经验方程，通常将多 孔介质中的孔隙直接看成圆形，此时ko取为2，并将KC常数直接取为5。
[0045] 但是对于该种玻璃微珠制成的多孔介质，该KC常数的选择并不准确，本专利测试 了 5种不同粒径分布共18片多孔介质的渗透率，根据实测数据，利用最小二乘法计算式（10) 中的KC常数，应为3.5686,与5相差较大。同时利用该模型计算该种由微珠构成的多孔介质 时误差较大，KC常数取3.5686时，相对误差（（理论值-实测值)/实测值*100%)最小，约为 36.8%。
[0046] 因此分别讨论形状因子ko和迂曲度τ的取值，W提高计算精度。由玻璃微珠构成的 多孔介质，微珠的球形特征保存的较好，如图1所示。因此，微珠之间的孔隙形状应该看成星 形，如图2所示，而不是圆形。利用Ansys-Fluent仿真环境计算星形孔隙的形状因子ko,得 到：
[0047] ko = 0.75 (11)
[004引对于迂曲度τ的取值，本专利认为其应该与多孔介质的孔隙率Φ和微珠粒径分布 偏度丫运两个因素有关。利用式（12)计算迂曲度，并分别分析迂曲度和孔隙率ΦΚ及微珠 粒径分布偏度丫之间的关系，实测实验结果如图3，图4所示。
[0049]
(1巧
[0050] 由实验数据通过最小二乘法拟合分别得到迂曲度和孔隙率W及微珠粒径分布偏 度的定量关系，如式(13)，（14)。
[0051] τι( Φ ) = 9.523Φ-1.029 (13)
[0052] Τ2( 丫）= 8.539 丫 3+33.9 丫 2-44.62 丫+21.15 (14)
[0053] 认为孔隙率与微珠粒径分布偏度对迂曲度的影响程度相同，如式(15)
[0054]
(15)
[0055] 将式(11)、（13)、（14)和(15)代入式(10)即可得式(9)。
[0056] 通过式(9)计算该种多孔介质的渗透率，理论值与实测值的相对误差为16.6%，明 显优于KC方程。与一些其他的渗透率计算方法相比，对于该种由玻璃微珠制成的多孔介质 的渗透率的计算，该种方法得到的理论值最准确。
[0057] 实施例2、步骤一:利用筛选好的玻璃微珠烧制5个圆柱形多孔介质样品，样品直径 为10mm，厚度为1~3mm不等。玻璃微珠密度为P = 2.46g/cm3。
[005引步骤二:测量5个多孔介质样品的质量和体积，利用式（1)计算得其孔隙率分别为： Φ 1 = 33.48%，(62 = 38.56%，(63 = 39.17%，(64=40.73%和 (65 = 43.49%。
[0059] 步骤测量玻璃微珠的基本参数：.1 =巧.12,歷，dpMed = 34.39皿= 0.03439mm，dp (0.25 )=28.25皿和dp (0.75 )=42.08皿。
[0060] 步骤四：利用式（2)~（8)计算玻璃微珠粒径分布参数：So = 1.061，寅=4 86i7， σ =0.42()1 ,Ei = 27.7皿瓜= 8.：MX102皿2，Ε3 = 2.75Χ1〇4μπι3, 丫 =0.9333。
[0061] 步骤五:利用式(9)分别计算出多孔介质样品渗透率：
[0062] Κι = 7.598 X l〇i3m2，Κ2 = 1.097 X l〇i2m2，Κ3 = 1.145 X l〇i2m2，Κ4= 1.275 X l〇i2m2，Κ已 = 1.540Xl〇i2m2。至此可W确定各个多孔介质样品的渗透率。
[0063] 综上，W上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在 本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护 范围之内。
【主权项】
1. 一种多孔介质渗透率的确定方法，其特征在于，所述多孔介质的孔隙率为Φ，组成所 述多孔介质的玻璃微珠的平均半径为A，玻璃微珠的粒径分布偏度为γ ;则多孔介质样品 的渗透率为：2. 如权利要求1所述的一种多孔介质渗透率的确定方法，其特征在于，所述多孔介质由 玻璃微珠烧结制备而成，玻璃微珠的粒径范围为20μπι~200μπι;制备方法为:制作圆柱形的 金属模具，将玻璃微珠填满模具并放入恒温炉进行加热，将微珠烧结在一起，获得圆柱形的 多孔介质，多孔介质的圆形端面直径为为6mm~15mm，厚度为1mm~5_。3. 如权利要求1所述的一种多孔介质渗透率的确定方法，其特征在于，所述多孔介质的 孔隙率为φ的测量方法为:测量获得多孔介质的总质量M、总体积V以及玻璃微珠的密度P， Μ 则孔隙率# = pV4. 如权利要求1所述的一种多孔介质渗透率的确定方法，其特征在于，测量玻璃微珠的 粒径分布参数:包括玻璃微珠的平均半径I，中位径d pMed，粒径累计分布参数dP(0.25WPdP (0.75);其中粒径累计分布参数表示小于该粒径的微珠占微珠总体积的体积分数；dP (0.25)为小于该粒径的所有微珠的体积总和占微珠体积总和的25%，dP(0.75)为小于该粒 径的所有微珠的体积总和占微珠体积总和的75% ; 所述多孔介质中玻璃微珠的分选系数为：当玻璃微珠的粒径分布满足对数正态分布时，其均值为歹=-和标准差分别为则所述玻璃微珠粒径分布偏度γ计算方法爻其中Ei (i = 1，2，3)为玻璃微珠粒径的i阶原点矩；
【专利摘要】本发明公开了一种多孔介质渗透率的确定方法，属于多孔介质技术领域，多孔介质的孔隙率为φ，组成所述多孔介质的玻璃微珠的平均半径为玻璃微珠的粒径分布偏度为γ；则多孔介质样品的渗透率为：<maths num="0001"></maths>该方法计算简单、精度较高。
【IPC分类】G01N15/08
【公开号】CN105628586
【申请号】CN201610094655
【发明人】付梦印, 程思源, 王美玲, 明丽, 李响, 肖梅峰
【申请人】北京理工大学
【公开日】2016年6月1日
【申请日】2016年2月19日