一种页岩气藏气体扩散系数实验测试方法
【技术领域】
[0001] 本发明涉及非常规油气储层气体传质能力实验评价方法技术领域，是一种模拟页 岩气藏气体（甲烷）扩散过程，测试扩散系数的实验方法。
【背景技术】
[0002] 页岩气是典型的非常规天然气，产自极低孔渗，以富有机质页岩为主的储集岩系 中，属经源岩气藏。我国页岩气可米资源量约为25. 08X 1012m3,与常规天然气资源量相当。 在页岩气藏中，页岩气主要以游离态和吸附态赋存。气井开井生产后，游离气首先产出；同 时随着压力波的传递，页岩有机质孔隙中的吸附气会逐渐解吸，并形成一个浓度梯度；在浓 度梯度的作用下，孔隙内部的吸附气会逐渐向已开始解吸的低浓度区域进行扩散传质。因 而，页岩气的产出过程是一个包含了解吸、扩散和渗流的多尺度传质过程。页岩气在有机质 内部的扩散过程不仅控制着吸附气的产出速率，还显著影响吸附气对整个气井产能的贡献 比例。研究页岩基块中气体扩散能力，有助于明确页岩气的多尺度产出机理，建立适合页岩 气井的产能模型。
[0003] 评价气体在岩石中扩散能力强弱的关键参数是气体扩散系数。然而目前关于页岩 中气体扩散系数的获取主要是通过计算机模拟和理论估算的方法，产能模型中对于扩散系 数的取值也缺乏相应的实验支撑。造成这一困境的原因主要是对页岩中气体扩散过程的认 识还不清楚，现有的扩散系数实验测试方法难以有效地反映页岩气在储层中的扩散过程， 并获取可靠的实验结果。
[0004] 现有的测量气体在岩石中扩散系数的方法主要有时滞法、烃浓度法和吸附/解吸 法。时滞法和烃浓度法常被应用于储层岩石中的天然气扩散运移和盖层岩石的封闭性研 究。时滞法测试的是在岩石饱和水条件下，单位岩石中水溶烃浓度梯度下，天然气通过岩石 的扩散系数。但通常认为原地条件下页岩气藏的含水饱和度极低或不含水，因而该方法并 不适用。烃浓度法测试的是压力平衡条件下双组份气体的相对扩散过程中某组分的扩散 系数，该方法通常采用干燥岩样，但在测试过程中需多次取样，极易造成较大的实验操作误 差。同时，时滞法和烃浓度法均未考虑吸附气的影响，且实验模拟过程与页岩气产出中气体 的扩散过程差异较大，难以被运用到页岩中气体扩散系数的实验测试。
[0005] 吸附/解吸法测试扩散系数在煤层气中应用较广，该方法通过分析吸附量与浓度 或压强变化的相应曲线计算扩散系数，能充分体现气体在岩石内部的吸附/解吸-扩散过 程。但吸附/解吸法也存在以下不足：首先，采用粉状样品（60~80目）进行测试，表征尺 度极小，且测得的扩散系数受样品的尺寸影响；其次，粉状样品制样繁琐，且无法反映有效 应力变化对孔隙结构以及扩散过程的影响。
[0006] 此外，在分子筛等多孔材料的研究中，零长柱（ZLC)色谱法、频率响应法、脉冲梯 度场核磁共振和拟弹性种子散射等方法也被应用于扩散系数的测试。但它们或者是实验条 件与工程应用差异极大，或者是仪器设备复杂昂贵，难以普及。
[0007] 综上，页岩气在储层中的扩散系数大小是评价气体在页岩中扩散传质能力强弱， 优化产能模型及预测吸附气产出程度的关键参数，面对现有实验测试方法存在的不足，亟 需探索一种能反映页岩气在储层中实际扩散过程并且易于操作的扩散系数测试方法。
【发明内容】
[0008] 针对现有扩散系数实验测试手段的不足，提供一种新的能模拟有效应力条件，反 映页岩气藏中气体实际扩散过程的扩散系数实验测试方法。该方法首先需要监测恒温定容 条件下气体（甲烷）向有效应力状态下柱塞页岩中的压力衰减过程，再提取其中反映气体 在页岩有机质内部浓度扩散的数据，根据建立的一维扩散模型求取扩散系数。
[0009] 本发明与已有方法相比具有如下优势：
[0010] (1)实验测试过程为页岩气藏气体产出的逆过程，能最大限度地还原吸附气产出 时在页岩有机质内部的解吸-扩散过程；
[0011] (2)相比粉状样品测试，本方法采用柱塞岩样，不仅可施加有效应力，还可以更为 全面地反映微孔结构对扩散过程的影响；
[0012] (3)相比已有的吸附/解吸法，本方法基于岩样参数剔除了气体自由膨胀阶段的 压力衰减数据，消除了其对扩散系数计算的影响。
[0013] 本发明的目的通过以下技术方案实现的：
[0014] 首先钻取、切割页岩岩样，充分烘干（60°C )至恒重后，测量其长度、直径和孔隙 度；实验前先检测装置气密性，再次对管线的死体积进行校准；而后将实验岩样放入夹持 器，对系统施加预设的温度和围压；待温度压力稳定后，对岩样及管线进行抽真空72h以 上；然后通入甲烷气体至微小中间容器，待甲烷气体预热至系统温度，且压力稳定在预设值 后，打开相应控制阀，甲烷开始向岩心中流动，同时开始采集整个腔室的压力衰减数据；待 系统压力24h内压降小于2KPa，认为气体的压力衰减过程达到近似平衡（IOOh以上），数据 采集结束。
[0015] 得到压力衰减曲线后，首先应根据岩样的尺寸和孔隙度求取实验温度压力条件下 岩样孔隙中能储存的甲烷质量，并以此为参照，剔除压力衰减曲线中反映早期气体向岩样 中的自由膨胀阶段的数据。该阶段以后，进入岩样的甲烷主要以吸附态的形式赋存，并且通 过浓度扩散向有机质内部运移。将此过程的压力衰减数据代入一维扩散系数计算模型，即 可求得气体（甲烷）在页岩中的扩散系数。
【附图说明】
[0016] 图1扩散系数实验测试装置示意图。
[0017] 图中：1-6.控制阀，7.高压气瓶，8.微小中间容器，9.真空栗，10.岩心夹持器， 11.围压栗，12.压力传感器，13.恒温温控系统，14.数据采集系统。
[0018] 图2实验测试过程气体（甲烷）的压力衰减曲线。
[0019] 图3气体（甲烷）扩散过程的压力衰减数据线性回归。
【具体实施方式】
[0020] 下面根据实例进一步说明本发明。
[0021] 选取四川盆地龙马溪组页岩柱塞岩样1块，采用本发明方法测试气体（甲烷）在 页岩中的扩散系数。具体步骤如下：
[0022] 步骤1、按照SYT5358-2010岩样制备方法准备柱塞页岩岩样（LMX-I);首先采用氮 气钻取并切割直径为2. 52cm，长度为3. 93cm的页岩岩样，而后放入60°C的烘箱，每隔24h 称重一次，直至相邻两次称重变化小于样品质量的1. 〇%，认为岩样充分烘干；烘干后的岩 样采用CMS300测试岩样孔隙度为4. 06% ;
[0023] 步骤2、将钢岩样（孔隙度为零）装入岩心夹持器中，并施加一定围压，开启恒温 温控系统，待将系统温度加热至设定温度T后，标定钢岩样端面至控制阀1和控制阀3间的 全部空间体积（V d)，标定方法为：关闭控制阀1，打开控制阀2-4,启动真空栗对系统内部抽 真空2h以上，然后关闭控制阀3,并关停真空栗；将已注入一定压力氦气的标准容器（已知 体积为V r)连接至控制阀1，待标准容器内压力稳定后，记录压力P1;打开控制阀1，待系统 内部压力稳定后，记录此时的压力P 2;根据方程（1)可求得岩样端面至控制阀1和控制阀3 间的全部空间体积VdS 58. 8362cm 3;
[0025] 式中Vd-岩样端面至控制阀1和控制阀3间的全部空间体积，cm3;
[0026] Vr一标准容器体积，cm3;
[0027] Z2-分别为压力P2下氦气的偏差因子，无因次；
[0028] 步骤3、将钢岩样装入岩心夹持器，并施加预设的围压IOMPa和温度82°C ;向系统 内通入一定压力为3MPa的氦气（气密性检测压力要高于实验压力），关闭控制阀1和控制 阀3, 24h内压力系统压力下降值小于0· 5KPa，说明气密性良好达到实验要求；
[0029] 步骤4、气密性检测完毕后，取出钢岩样，装入待分析页岩柱塞岩样，再次通过围压 栗11和恒温温控系统13施加预设的围压IOMPa和温度82°C，并保持围压和温度恒定；然 后关闭控制阀1，打开控制阀2-4,开启真空栗9对岩样和管线系统抽真空72h ;然后关闭控 制阀2-4,打开控制阀1，通入一定量的气体（甲烷）至微小中间容器8,待甲烷压力升至预 设压力后，关闭控制阀1 ;当通入的甲烷气体预热至系统温度，且压力保持稳定后，此时气 体压力为1420KPa，打开控制阀2,同时通过压力传感器12和数据采集系统14记录甲烷向 岩样流动过程的气体压力衰减数据；待系统压力24h内压降小于2KPa，认为气体压力达到 平衡；如图2所示，本次实验取压力平衡时间为第IOOh ;
[0030] 步骤5、甲烷向已充分抽真空的岩样中流动，首先气体会通过自由膨胀充满页岩中 的孔隙和微裂缝系统，而后甲烷会进一步吸附在页岩中的有机质纳米孔壁面；吸附气在浓 度梯度的作用下，逐渐向有机质内部进行扩散；鉴于页岩中气体流动的多尺度特征，监测的 压力衰减数据包含了前期甲烷向岩样孔隙中的自由膨胀阶段，而在扩散系数的计算过程中 需要去除该阶段数据；具体处理过程如下，t时刻自由膨胀进入岩样的甲烷的质量可由方 程⑵求得：
[0032] 式中- t时刻自由膨胀进入岩样的甲烷的质量，即此刻的累计储气量，g ;
[0