[发明专利]一种颗粒样品两级渗透率的获取方法

权利要求书

1.一种颗粒样品两级渗透率的获取方法，其特征在于，所述方法包括：

获取压力随时间变化的曲线P(t)、基准腔体积V_b、放入样品前样品腔体积V_s、样品总体体积V_sb、颗粒平均半径R_a、测试温度T、测试阀开启前基准腔的气体密度ρ₁和样品腔的气体密度ρ₀；

基于氦气状态方程、压力衰减曲线P(t)最终稳定值P₂和温度T，计算对应的气体密度ρ₂、样品孔隙体积V_P、孔隙度φ、颗粒样品外部空间和内部孔隙空间的体积比K_c和测试阀未开启时颗粒外部空间的气体平均密度ρ_c0；

确定气体开始进入样品对应的压力P₃，从所述压力随时间变化的曲线P(t)中确定压力P₃对应的时刻t₃；

选择所有大于所述时刻t₃的测点，基于氦气状态方程计算气体密度ρ(t)和质量分数F_R(t)，对ln(F_R(t))中所有大于所述时刻t₃的测点和时间t进行线性拟合，得到斜率s₁，根据所述斜率s₁计算渗透率k₀；

确定渗透率k′，设定i的初始值为4，以第i个测点为界，将压力衰减曲线分为两段：第i个测点及以前的测点记为A段，压力衰减曲线为P_A(t)；第i个测点及以后的测点记为B段，压力衰减曲线为P_B(t)；

基于氦气状态方程、第i个测点的压力P_i和温度T，计算对应的气体密度ρ_i，计算A段的孔隙体积V_pA、孔隙度φ_A(i)、颗粒样品外部空间和内部孔隙空间的体积比K_cA和测试阀未开启时颗粒外部空间的气体平均密度ρ_c0A，输入(k_A(i)，D_2A(i))＝f(P_A(t)，K_cA，ρ_c0A，ρ₀，R_a，φ_A(i)，k₀，k′)求得对应的渗透率最优解k_A(i)和A段最小偏差平方和D_2A(i)；

计算B段的孔隙体积V_pB、孔隙度φ_B(i)、颗粒样品外部空间和内部孔隙空间的体积比K_cB和测试阀未开启时颗粒外部空间的气体平均密度ρ_c0B，输入(k_B(i)，D_2B(i))＝f(P_B(t)，K_cB，ρ_c0B，ρ₀，R_a，φ_B(i)，k₀，k′)求得对应的渗透率最优解k_B(i)和B段最小偏差平方和D_2B(i)；

确定k_A(i)和k_B(i)为两级渗流通道的渗透率值。

2.根据权利要求1所述的方法，其特征在于，利用下述公式，基于氦气状态方程、压力衰减曲线P(t)最终稳定值P₂和温度T，计算对应的气体密度ρ₂、样品孔隙体积V_P、孔隙度φ、颗粒样品外部空间和内部孔隙空间的体积比K_c和测试阀未开启时颗粒外部空间的气体平均密度ρ_c0：

3.根据权利要求1所述的方法，其特征在于，所述确定气体开始进入样品对应的压力P₃，包括：

利用公式确定

基于氦气状态方程，由ρ₃和温度T计算P₃。