一种计算气井煤层气吸附量的方法
【专利摘要】一种计算气井煤层气吸附量的方法，包括以下步骤：步骤一：对气井穿过的煤层进行测试分析，获取煤质工业分析结果；收集气井煤层数据；步骤二：根据煤的干燥无灰基挥发分计算镜质组最大反射率；步骤三：根据煤层气动态吸附方程计算各煤层的煤层气吸附量；步骤四：按各煤层地层温度，计算各煤层不同压力下的吸附量，将该值求和便得到气井穿过所有煤层的总吸附量。该方法可大幅度提高气井煤层气吸附量计算的精度，对指导煤层气勘探开发具有重要的实际意义。
【专利说明】
一种计算气井煤层气吸附量的方法
技术领域：
[0001] 本发明涉及煤层气勘探与开发领域，特别涉及一种计算气井煤层气吸附量的方 法。
【背景技术】：
[0002] 气井煤层吸附气含量的定量计算方法是煤层气储层评价的关键和主要技术难题， 利用测井资料计算其含量成为重要途径，其计算精度直接影响到煤层气的勘探开发决策。 煤层临界解吸压力是指解吸与吸附达到平衡时对应的压力，即压力降低使吸附在煤微孔隙 表面上的气体开始解吸时的压力。理论上，当储层压力降低到临界解吸压力以下时，煤孔隙 中吸附的气体开始解吸，向裂隙方向扩散，在压力差的作用下，从裂隙向井筒流动，且煤层 的吸附和解析是可逆的。在煤层的解析过程中，当压力降低时，气体则开始解析。目前，煤层 气开采大多是基于这一原理，通过排水降低压力而达到采气的目的。
[0003] 目前气井煤层气吸附量的计算中把各个煤储层当作是均匀的产气储层，用 Langmuir模型、Freundlich经验模型、BP神经网络模型等计算出气井煤层气的吸附量。然而 煤层气井大都穿过多煤层气储层或厚度较大的煤层气储层，并往往统一进行混合排采，把 各个煤储层当作是均匀的产气储层，这样计算出的吸附量显然不符合实际情况，上述计算 模型仅仅考虑压力对煤层气吸附量的影响，并未考虑温度对吸附气量的影响，实际上随着 深度的改变，煤层压力与温度同时改变，且压力与温度对煤层气吸附量都有影响，因此上述 模型计算煤层气吸附量存在一定误差。
【发明内容】
：
[0004] 有鉴于此，本发明提供一种计算气井煤层气吸附量的方法，该方法考虑气井穿过 多煤层的情况，同时给出一个既考虑温度变化又考虑压力变化的动态吸附方程，大幅度提 高气井煤层气吸附量计算的精度，对指导煤层气勘探开发具有重要的实际意义。
[0005] -种计算气井煤层气吸附量的方法，包括以下步骤：
[0006] 步骤一:对气井穿过的煤层进行测试分析，获取煤质工业分析结果;收集气井煤层 数据，包含:气井穿过的煤层数，各煤层深度、煤厚、温度、压力；
[0007]步骤二:根据煤的干燥无灰基挥发分计算镜质组最大反射率；
[0008] 步骤三:根据煤层气动态吸附方程计算各煤层的煤层气吸附量；
[0009] 步骤四：按各煤层地层温度，计算各煤层不同压力下的吸附量，将该值求和便得到 气井穿过所有煤层的总吸附量。
[0010] 计算气井煤层气吸附量的方法具体步骤为：
[0011]步骤一:对气井穿过的煤层进行测试分析，获取煤质工业分析结果;收集气井煤层 数据，包含:气井穿过的煤层数，各煤层深度、煤厚、温度、压力；
[0012] 进一步，地层温度可按地温梯度2.33°C/100m计算，地层压力梯度可按IMPa/lOOm 计算。
[0013] 步骤二:根据煤的干燥无灰基挥发分与镜质组最大反射率之间的线性关系计算各 煤层镜质组最大反射率，该线性关系用式1表示：
[0014] ^max ~ * Vdaf 式.1
[0015] 式中：式_ :镜质组最大反射率；
[0016] ai、bi:待定系数；
[0017] Vdaf:干燥无灰基挥发分；
[0018]用最小二乘法回归出上式的待定系数&1、1^，根据气井各煤层煤的干燥无灰基挥发 分数据，用上式计算各煤层镜质组最大反射率。
[0019] 步骤三:根据煤层气动杰吸附方程计算各煤层的煤层气吸附量，如式2:
式2
[0021] 式中：V:吸附量，m3/t;
[0022] T:地层温度，K;
[0023] B:吸附流量系数；
[0024] A :吸附分子能量差；
[0025] P:地层压力，MPa;
[0026] 吸附常数；
[0027] 进一步，利用步骤二得出的镜质组最大反射率计算煤层气动态吸附方程的常数B、 A、0，如式3,式4,式5:
A. 3 式4 式5
[0031 ]式中：&2、匕2、。2、(12、&3、匕3、03、(13、&4、匕4、04、(14:回归常数；
[0032] 按最小二乘法回归出式3，式4,式5中的32、匕2、02、(12、33、匕3、03、(13、34、匕4、04、(14回归 常数，并根据^：值求煤层气动态吸附方程中B、A的常数。
[0033]步骤四：按各煤层地层温度，计算各煤层不同压力下的吸附量，将该值求和便得到 气井穿过所有煤层的总吸附量;并绘制气井全部煤层吸附曲线。
[0034]进一步，对于气井不同煤层，其最大的吸附压力P应小于P1>max，若第i煤层地层压力 大于Pi,max，则该层煤层气吸附量为零;Pi,max可由式6计算：
[0035] Pi, max = Po+ (hi-h〇) Pg 式 6
[0036]式中:Pi,max:第i煤层最大吸附压力；
[0037] Po:套管封头压力，MPa;
[0038] hi:第i煤层的埋深，m;
[0039] ho:套管内液位高度，m;
[0040] p:液体密度，kg/m3;
[00411 g:重力加速度。
[0042] -种计算气井煤层气吸附量的方法有益效果在于:该方法考虑气井穿过多煤层的 情况，同时考虑煤层气进行混合排采的情况，并给出一个既考虑温度变化又考虑压力变化 的动态吸附方程，该方程根据不同变质程度的煤得出不同的常数，即考虑煤的变质程度对 吸附量的影响，因此该方法可大幅度提高气井煤层气吸附量计算的精度，对指导煤层气勘 探开发具有重要的实际意义。
【具体实施方式】：
[0043]下面结合具体实施例对本发明作进一步说明：
[0044]步骤一:选择某气井孔号为502的勘探井为计算对象，对该气井穿过的煤层进行测 试分析，获取煤质工业分析结果，如下表所示。
[0046]收集气井煤层数据，包含:气井穿过的煤层数，各煤层见煤深度、煤厚、温度、压力， 如下表所示。
[0048] 步骤二:煤的干燥无灰基挥发分与镜质组最大反射率之间的关系如式1:
[0049] R^^^+b*Vdi,f 式工
[0050] 式中：:镜质组最大反射率；
[0051] ai、bi:待定系数；
[0052] Vdaf:干燥无灰基挥发分；
[0053]用最小二乘法回归出上式中的待定系数ai、h，根据气井各煤层煤的干燥无灰基挥 发分数据，用上式计算各煤层镜质组最大反射率;计算结果如下表。
[0055]步骤三:根据煤层气动态吸附方程计算各煤层的煤层气吸附量，如式2:
式2
[0057] 式中：V:吸附量，m3/t;
[0058] T:地层温度，K;
[0059] B:吸附流量系数；
[0060] A :吸附分子能量差；
[0061 ] P:地层压力，MPa;
[0062] 吸附常数；
[0063]进一步，利用步骤二得出的镜质组最大反射率计算煤层气动态吸附方程的常数B、 A、0，如式3,式4,式5:
[0067]式中：82、匕2、02、(12、&3、匕3、03、(13、&4、匕4、04、(14:回归常数；[0068] 按最小二乘法回归出式3，式4,式5中的32、匕2、02、(12、33、匕3、03、(13、34、匕4、04、(14回归常数，并根据^；值求煤层气动态吸附方程中B、A的常数。 式3 式4 式5
[0070] 步骤四：按各煤层地层温度，计算各煤层不同压力下的吸附量，将该值求和便得到 气井穿过所有煤层的总吸附量;并绘制气井全部煤层吸附曲线。
[0071] 进一步，对于气井不同煤层，其最大的吸附压力P应小于Pi,max，若第i煤层地层压力 大于Pi, max，则该层煤层气吸附量为零;Pi,max可由式6计算：
[0072] Pi, max = Po+ (hi-h〇) Pg 式 6
[0073] 式中:Pi,max:第i煤层最大吸附压力；
[0074] Po:套管封头压力，MPa;
[0075] hi:第i煤层的埋深，m;
[0076] ho:套管内液位高度，m;
[0077] P:液体密度，kg/m3;
[0078] g:重力加速度；
[0079]该气井通过排水降低压力得到煤层气的总吸附量计算结果如下表所示。
【主权项】
1. 一种计算气井煤层气吸附量的方法，其特征在于:一种计算气井煤层气吸附量的方 法包括以下步骤： 步骤一：对气井穿过的煤层进行测试分析，获取煤质工业分析结果；收集气井煤层数 据，包含:气井穿过的煤层数，各煤层深度、煤厚、温度、压力； 步骤二:根据煤的干燥无灰基挥发分计算镜质组最大反射率； 步骤三:根据煤层气动态吸附方程计算各煤层的煤层气吸附量； 步骤四：按各煤层地层温度，计算各煤层不同压力下的吸附量，将该值求和便得到气井 穿过所有煤层的总吸附量。2. 如权利要求1所述一种计算气井煤层气吸附量的方法，其特征在于:步骤二中计算各 煤层镜质组最大反射率方程如式1: ^max： = + h * Klaf 式 1 式中：镜质组最大反射率； ai、bi:待定系数； Vdaf :干燥无灰基挥发分； 用最小二乘法回归出上式的待定系数31士1，根据气井各煤层煤的干燥无灰基挥发分数 据，计算各煤层镜质组最大反射率。3. 如权利要求1所述一种计算气井煤层气吸附量的方法，其特征在于:步骤三中煤层气 动态吸附方程如式2:式中:V:吸附量，m3/t; T:地层温度，K; B:吸附流量系数； A :吸附分子能量差； P:地层压力，MPa; β:吸附常数； 利用步骤二得出的镜质组最大反射率计算煤层气动态吸附方程的常数Β、△、β，如式3， 式4,式5:式中：82、匕2、。2、(12、&3、匕3、03、(13、&4、匕4、04、(14:回归常数； 按最小二乘法回归出式3，式4,式5中的82、匕2、02、(12、33、匕3、03、(13、34、匕4、04、(14回归常数， 并根据U值求煤层气动态吸附方程中B、A j的常数。4.如权利要求1所述一种计算气井煤层气吸附量的方法，其特征在于:步骤四中，对于 气井不同煤层，其最大的吸附压力P应小于P1>max，若第i煤层地层压力大于P1>max，则该层煤 层气吸附量为零;Pi,max可由式6计算： Pi,max = P〇+(hi-ho)pg 式 6 式中:Pi,max:第i煤层最大吸附压力； Po:套管封头压力，MPa; hi:第i煤层的埋深，m; ho:套管内液位高度，m; P:液体密度，kg/m3; g:重力加速度。
【文档编号】G06F19/00GK105930667SQ201610269095
【公开日】2016年9月7日
【申请日】2016年4月27日
【发明人】张学梅, 李 东, 魏亚玲, 苗泽凯, 程实
【申请人】中国矿业大学银川学院