1.一种压力修正式气体渗透率计算方法,包括步骤1)、气体渗透率计算:测量时,将试样(7)放置在压力室(6)的缸筒内部,并用氟橡胶套(8)包裹,之后打开气源与缓冲钢瓶(3)之间的阀门A(2),向缓冲钢瓶(3)中注入一定的压力P1,通过缓冲钢瓶(3)与压力室(6)之间的压力表(5)可以读出其值,之后关闭阀门A(2),打开阀门B(4),开始向试样(7)注入气体,从压力表(5)可以读出这一过程中的压力变化,经过Δt时刻,压力表(5)读数变为P1’,那么在Δt时间内,缓冲钢瓶(3)内的平均压力及压力变化为:ΔP=P1‑P1’ (2),其中:P1为初始注入缓冲钢瓶(3)内的压力;P1’为在T1温度下压力表(5)测量的实时压力值,T1为试验时的实时温度;ΔP为打开阀门B(4)向试样(7)注气后,Δt时间内变化的压力值;Δt为变化的时间;Pmean为平均压力;
根据达西定律,可以知道在此时间段内通过的平均流量为其中:k为渗透率,A为流体的截面面积,μ为动粘度, 为对气压的偏导;Qmean为平均流量;
而根据理想气体状态方程,有ΔPV0=PmeanQmeanΔt (4),即
其中,V0为缓冲钢瓶(3)体积+缓冲钢瓶(3)到压力室(6)之间管路体积;
可压缩性气体在试样(7)中的分布服从如下函数:其中,h为试样高度,t为时间,x为气体注气后在试样中流动的距离;
对式(6)进行求导有:
当测量试样入口端的气体渗透率,即气体进入渗透率,则取x=0,有其中:P0为一个大气压,其余符号与在先定义相同;
带入式(3)有
将式(9)与式(5)联立则可以得到入口端气体有效渗透率最终得到试样(7)的有效渗透率k;其特征是,还包括步骤2)、气体压力修正计算:根据理想气体状态方程有:
PV=nRT (11),其中,P为理想气体的压强,V为理想气体的体积,n为气体物质的量,R为理想气体常数,T为热力学温度;
由此得知,当体积恒定的时候,在T1温度下有:P1’V=nRT1 (12),在T2温度下有:
P2V=nRT2 (13),P2为在T2温度下的气压;T2为试验设定的温度;
由此得到:
则有:
之后将换算得到的压力P2作为P1’代入,重新推导计算渗透率k,即消除温度的影响,得到准确的气体渗透率k。