《一种基于启动压力梯度的煤层瓦斯流动状态快速判识方法.pdf》由会员分享，可在线阅读，更多相关《一种基于启动压力梯度的煤层瓦斯流动状态快速判识方法.pdf（13页珍藏版）》请在专利查询网上搜索。
1、10申请公布号CN103670392A43申请公布日20140326CN103670392A21申请号201310746879X22申请日20131230E21B49/0020060171申请人河南理工大学地址454000河南省焦作市高新区世纪大道2001号河南理工大学能源科学与工程学院72发明人韩颖张飞燕王革学苏现波倪小明郭红玉程虹铭宋德尚朱林剑王博曹文涛周玉军余伟凡杨志龙任艳普张平生74专利代理机构南昌新天下专利商标代理有限公司36115代理人施秀瑾54发明名称一种基于启动压力梯度的煤层瓦斯流动状态快速判识方法57摘要本发明涉及一种基于启动压力梯度的煤层瓦斯流动状态快速判识方法，通过在煤矿。
2、井下采集五类煤样，制成标准煤芯，确定出每个煤芯的地质强度指标值，然后再测算出煤芯的启动压力梯度值，并对多种煤样进行一系列的启动压力梯度测试，得出各类煤样的启动压力梯度值，然后将各类煤样的启动压力梯度值与其对应的地质强度指标值进行回归分析，得出二者的关系；最后得出结论，根据在煤矿井下采集煤样或者直接观测煤壁，获取煤层的煤体结构，确定煤层的GSI值后，即可得出煤层瓦斯的启动压力梯度，将其与瓦斯压力梯度进行比较，即可快速判定煤层瓦斯的流动状态。本发明为抽采难易程度评价及抽采工艺的选择提供指导，并为瓦斯抽采理论的进一步完善提供参考。51INTCL权利要求书1页说明书5页附图6页19中华人民共和国国家知。
3、识产权局12发明专利申请权利要求书1页说明书5页附图6页10申请公布号CN103670392ACN103670392A1/1页21一种基于启动压力梯度的煤层瓦斯流动状态快速判识方法，其特征在于，包括如下步骤第一步，在煤矿井下采集具有代表性的I类非破坏煤、II类破坏煤、III类强烈破坏煤、IV类粉碎煤、V类全粉煤煤样；第二步，在实验室内将各煤样制成规格相同的标准煤芯；第三步，确定每个标准煤芯的地质强度指标GSI值；第四步，根据达西定律，得出存在启动压力梯度时的气体渗流方程，采用启动压力测试系统，测出一系列V与P12P22的值，并将多组V与P12P22的值进行回归分析，计算出常数A、B值，代入推算。
4、出的启动压力梯度的计算公式，得出的值，其中，V为气体通过煤芯的流速，M/S；K为煤芯渗透率，M2；P1为入口气体压力，PA；P2为出口气体压力，PA；P0为大气压力，101325PA；为气体粘度，PAS；L为气体流经长度，即煤芯高度，M；第五步，对多种煤样进行一系列的IV类煤样的启动压力梯度测试，得出各类煤样的启动压力梯度值，然后将各类煤样的启动压力梯度值与其对应的地质强度指标GSI值进行回归分析，得出二者的关系；第六步，得出结论，根据在煤矿井下采集煤样或者直接观测煤壁，获取煤层的煤体结构，确定煤层的GSI值后，即可得出煤层瓦斯的启动压力梯度，将与瓦斯压力梯度P/L进行比较，即可快速判定煤层瓦。
5、斯的流动状态P/L，为低速非线性渗流；P/L，为扩散；0，为线性渗流。2根据权利要求1所述的基于启动压力梯度的煤层瓦斯流动状态快速判识方法，其特征在于在第三步中，将每个煤芯与煤体结构量化的GSI图版进行比对，确定其地质强度指标GSI值。3根据权利要求1所述的基于启动压力梯度的煤层瓦斯流动状态快速判识方法，其特征在于所述标准煤芯为圆柱形煤芯，其直径为50MM，高度为50MM或100MM。4根据权利要求1所述的基于启动压力梯度的煤层瓦斯流动状态快速判识方法，其特征在于所述启动压力梯度测试系统包括一个用于放置煤样的煤样罐，所述煤样罐的两端均设有试验机压头，煤样设置于两个试验机压头之间，还包括一个用于。
6、提供高压N2气源的高压气罐，煤样罐的两端设有气道，进气口与高压气罐连通，排气口与大气连通，所述高压气罐的气路上设有减压稳压阀和气体压力传感器，所述压力传感器与静态电阻应变仪和数据处理仪连接，所述静态电阻应变仪与数据处理仪通讯采集入口气体压力，所述煤样罐的排气管路上设有气体质量流量计，所述气体质量流量计与气体流量积算仪和数据采集仪连接。5根据权利要求4所述的基于启动压力梯度的煤层瓦斯流动状态快速判识方法，其特征在于所述煤样与试验机压头之间还设有透气板。权利要求书CN103670392A1/5页3一种基于启动压力梯度的煤层瓦斯流动状态快速判识方法技术领域0001本发明涉及一种煤层瓦斯流动状态判识方。
7、法，尤其是一种基于启动压力梯度的煤层瓦斯流动状态快速判识方法。背景技术0002瓦斯的高效抽采是灾害治理和资源利用的根本性措施之一，其效果主要取决于煤层与外部空间联系的纽带渗透率。我国煤层渗透率一般在1MD以下，具有“低渗”特点。长期的低渗透油气藏开发实践表明，流体在低渗多孔介质中的流动不再遵循达西定律，而是带有启动压力梯度的低速非线性渗流。启动压力梯度概念由前苏联学者BA费劳林于1951年提出，他认为只有当实际压力梯度大于某一临界值时，流动才能发生，此临界值称为启动压力梯度。0003研究表明当瓦斯在煤层中流动时，若瓦斯压力梯度（煤层埋藏深度每增加1M，煤层瓦斯压力的平均增加值，一般为00050。
8、015MPA/M）大于启动压力梯度，瓦斯将发生低速非线性渗流；反之将不流动，此时瓦斯仅通过扩散途径产出；只有当启动压力梯度为零时，瓦斯才会发生线性渗流。但是，现行瓦斯抽采技术大多遵循线性渗流理论达西定律，往往忽略了低渗煤层内存在的低速非线性渗流及扩散现象，导致抽采难易程度评价及抽采工艺选择出现偏差。针对具有不同煤体结构的煤层，如何快速、准确地评价其抽采难易程度以及选择适宜的抽采工艺，以实现瓦斯高效抽采，是煤矿现场亟待解决的问题。而解决该问题的前提及关键，在于对瓦斯在不同煤体结构煤层内流动状态的准确把握。0004本发明中涉及到的若干专业术语做出如下解释压力梯度沿流体流动方向，单位路程长度上的压力。
9、变化，可用增量形式P/L或微分形式DP/DL表示，式中P为压力；L为距离。0005渗透率压力梯度为1时，动力黏滞系数为L的液体在介质中的渗透速度，渗透率是表征土或岩石本身传导液体能力的参数，其大小与孔隙度、液体渗透方向上空隙的几何形状、颗粒大小以及排列方向等因素有关，而与在介质中运动的液体性质无关。渗透率（K）用来表示渗透性的大小。在一定压差下，岩石允许流体通过的性质称为渗透性；在一定压差下，岩石允许流体通过的能力叫渗透率。0006启动压力梯度只有当实际压力梯度大于某一临界值时，流动才能发生，此临界值称为启动压力梯度。0007达西定律反映水在岩土孔隙中渗流规律的实验定律，其表达式为QKFH/L。
10、，式中Q为单位时间渗流量，F为过水断面，H为总水头损失，L为渗流路径长度，IH/L为水力坡度，K为渗透系数。关系式表明，水在单位时间内通过多孔介质的渗流量与渗流路径长度成反比，与过水断面面积和总水头损失成正比。从水力学已知，通过某一断面的流量Q等于流速V与过水断面F的乘积，即QFV。或，据此，达西定律也可以用另一种形式表达VKI，V为渗流速度。上式表明，渗流速度与水力坡度一次方成正比。说明水力坡度与渗流速度呈线性关系，故又称线性渗流定律。达西定律适用的上限有两种看法一种认为达西定律适用于说明书CN103670392A2/5页4地下水的层流运动；另一种认为并非所有地下水层流运动都能用达西定律来表。
11、述，有些地下水层流运动的情况偏离达西定律，达西定律的适应范围比层流范围小。这个定律说明水通过多孔介质的速度同水力梯度的大小及介质的渗透性能成正比。发明内容0008本发明的目的在于克服现有技术的不足，适应现实需要，提供一种基于启动压力梯度的煤层瓦斯流动状态快速判识方法。0009为了实现本发明的目的，本发明采用的技术方案为一种基于启动压力梯度的煤层瓦斯流动状态快速判识方法，其特征在于，包括如下步骤第一步，在煤矿井下采集具有代表性的I类非破坏煤、II类破坏煤、III类强烈破坏煤、IV类粉碎煤、V类全粉煤煤样；第二步，在实验室内将各煤样制成规格相同的标准煤芯；第三步，确定每个标准煤芯的地质强度指标GS。
12、I值；第四步，根据达西定律，得出存在启动压力梯度时的气体渗流方程，采用启动压力测试系统，测出一系列V与P12P22的值，并将多组V与P12P22的值进行回归分析，计算出常数A、B值，代入推算出的启动压力梯度的计算公式，得出的值，其中，V为气体通过煤芯的流速，M/S；K为煤芯渗透率，M2；P1为入口气体压力，PA；P2为出口气体压力，PA；P0为大气压力，101325PA；为气体粘度，PAS；L为气体流经长度，即煤芯高度，M。0010第五步，对多种煤样进行一系列的IV类煤样的启动压力梯度测试，得出各类煤样的启动压力梯度值，然后将各类煤样的启动压力梯度值与其对应的地质强度指标GSI值进行回归分析，。
13、得出二者的关系；第六步，得出结论，根据在煤矿井下采集煤样或者直接观测煤壁，获取煤层的煤体结构，确定煤层的GSI值后，即可得出煤层瓦斯的启动压力梯度，将与瓦斯压力梯度P/L进行比较，即可快速判定煤层瓦斯的流动状态P/L，为低速非线性渗流；P/L，为扩散；0，为线性渗流。0011在第三步中，将每个煤芯与煤体结构量化的GSI图版进行比对，确定其地质强度指标GSI值。0012所述标准煤芯为圆柱形煤芯，其直径为50MM，高度为50MM或100MM。0013所述启动压力梯度测试系统包括一个用于放置煤样的煤样罐，所述煤样罐的两端均设有试验机压头，煤样设置于两个试验机压头之间，还包括一个用于提供高压N2气源的。
14、高压气罐，煤样罐的两端设有气道，进气口与高压气罐连通，排气口与大气连通，所述高压气罐的气路上设有减压稳压阀和气体压力传感器，所述压力传感器与静态电阻应变仪和数据处理仪连接，所述静态电阻应变仪与数据处理仪通讯采集入口气体压力，所述煤样罐的排气管路上设有气体质量流量计，所述气体质量流量计与气体流量积算仪和数据采集仪连说明书CN103670392A3/5页5接。0014所述煤样与试验机压头之间还设有透气板。本发明的有益效果在于1建立基于启动压力梯度的煤层瓦斯流动状态快速判识方法，可实现不同煤体结构煤层瓦斯流动状态的快速判识，查明瓦斯在低渗煤层内的非达西渗流特征，为抽采难易程度评价及抽采工艺的选择提供。
15、指导，并为瓦斯抽采理论的进一步完善提供参考和借鉴。00152可为地面煤层气井产能预测提供理论依据，提高产能预测精度，对于指导煤层气开发具有重要的现实意义和理论价值。附图说明0016图1为煤的破坏类型分类表；图2为煤体结构量化的GSI图版；图3为启动压力梯度测试系统的结构示意图；图4为本案实施例中经实测得出的若干组V与P12P22的值；图5为对图4中V与P12P22进行回归分析得出的分析结果示意图；图6为IV类煤样启动压力梯度测试数据；图7为将图6中的启动压力梯度值与地质强度指标GSI进行回归分析得到的分析结果示意图。具体实施方式0017下面结合附图和实施例对本发明进一步说明实施例参见图17。0。
16、018一种基于启动压力梯度的煤层瓦斯流动状态快速判识方法，其特征在于，包括如下步骤第一步，在煤矿井下采集具有代表性的I类非破坏煤、II类破坏煤、III类强烈破坏煤、IV类粉碎煤、V类全粉煤煤样；关于五类煤样的判断标准，已在行业标准煤与瓦斯突出矿井鉴定规范（AQ10242006）中明确规定，其他文献也多有公开。图1选自煤与瓦斯突出矿井鉴定规范（AQ10242006），为煤的破坏类型分类表。0019第二步，在实验室内将各煤样制成规格相同的标准煤芯，本实施例中，标准煤芯为圆柱形煤芯，其直径为50MM，高度为50MM或100MM；第三步，确定每个标准煤芯的地质强度指标GSI值，本实施例中，采用的方法是。
17、，将每个煤芯与煤体结构量化的GSI图版进行比对，确定其地质强度指标GSI值；GSI图版是现有技术中已经发表过的东西，在此不赘述。0020选取GSI图版是一种优选的实施方式，还可采用其他现有技术来确定每个标准煤芯的地质强度指标GSI值。0021第四步，根据达西定律，得出存在启动压力梯度时的气体渗流方说明书CN103670392A4/5页6程。采用启动压力测试系统，测出一系列V与P12P22的值，并将多组V与P12P22的值进行回归分析，计算出常数A、B值，代入推算出的启动压力梯度的计算公式，得出的值，其中，V为气体通过煤芯的流速，M/S；K为煤芯渗透率，M2；P1为入口气体压力，PA；P2为出口。
18、气体压力，PA；P0为大气压力，101325PA；为气体粘度，PAS；L为气体流经长度，即煤芯高度，M。0022图3示出了启动压力梯度测试系统的结构及其基本工作原理；所述启动压力梯度测试系统包括一个用于放置煤样9的煤样罐7，所述煤样罐7的两端均设有试验机压头6，煤样9设置于两个试验机压头6之间，还包括一个用于提供高压N2气源的高压气罐1，煤样罐7的两端设有气道，进气口与高压气罐1连通，排气口与大气连通，所述高压气罐1的气路上设有减压稳压阀2和气体压力传感器3，所述压力传感器3与静态电阻应变仪4和数据处理仪5连接，所述静态电阻应变仪4与数据处理仪通讯采集入口气体压力，所述煤样罐7的排气管路上设有。
19、气体质量流量计10，所述气体质量流量计10与气体流量积算仪11和数据采集仪12连接。0023所述煤样9与试验机压头6之间还设有透气板8。002410气体质量流量计即时采集通过煤芯的气体质量流量，并与11气体流量积算仪连接，气体瞬时流量及累计流量同时显示在11气体流量积算仪的电子屏幕上。002512数据采集仪对实验过程中的气体流量变化规律自动进行分析，并以曲线形式展现。0026图中，13为声发射探头，14为声发射数据记录仪，15为煤的应力应变计算机处理系统，这几个部件与本发明内容无直接关系，不赘述。0027实验步骤如下将圆柱形标准煤芯放入煤样罐内，利用RMT150B伺服试验机通过试验机压头6对煤。
20、样9加压，施加2KN轴向力，确保煤样罐密封良好，利用真空泵抽真空8H。0028利用RMT150B伺服试验机通过试验机压头6对煤样9继续加压，施加4KN轴向力，2MPA围压。0029打开高压气罐1，将入口气体压力P1调至较高值（为确保测试系统密封良好，测试过程中气体压力不得超过2MPA），关闭减压稳压阀2，气体质量流量计10与气体流量积算仪11采集的数据，经数据采集仪12换算自动生成流经圆柱形标准煤芯的气体流速V。需要说明的一点是，测试系统排气口端直通大气，故出口气体压力P2P0101325PA。0030按上述方法，逐步降低P1值，测试相应的V与P12P22。0031根据达西定律，当不考虑启动压。
21、力梯度时的气体渗流方程为（1）说明书CN103670392A5/5页7式中，V为气体通过煤芯的流速，M/S；K为煤芯渗透率，M2；P1为入口气体压力，PA；P2为出口气体压力，PA；P0为大气压力，101325PA；为气体粘度，PAS；L为气体流经长度，即煤芯高度，M。0032由式（1）可以看出，当气体渗流符合达西定律时，V与P12P22为通过原点的线性关系，除了V与P12P22外，其它参数都为常数、定值，所以公式（1）可以写成VAP12P22。0033当存在启动压力梯度时，气体渗流方程必然不会通过原点，此时，气体渗流方程为（2）式中，A、B为常数，实验测试一系列V与P12P22，通过回归分析。
22、，确定A、B值。0034令V0，则（3）因此，启动压力梯度计算公式为（4）式中，为启动压力梯度，PA/M，L为气体流经长度，即煤芯高度，M；通过图3所示的启动压力测试系统测出一系列V与P12P22的值，本案中选取一种煤样进行实测，得到如图4所示的结果。0035然后对图4中的结果进行回归分析，得到如图5所示的关系示意图，由图5得出，A1276981015，B246874106，将其代入式（4），得出0183MPA/M。0036第五步，对多种煤样进行一系列的IV类煤样的启动压力梯度测试，得出各类煤样的启动压力梯度值，然后将各类煤样的启动压力梯度值与其对应的地质强度指标GSI值进行回归分析，得出二者。
23、的关系。0037对多种煤样进行启动压力测试得到的值详见图6。0038将各煤样的值与地质强度指标GSI进行回归分析，得到的分析结果如图7所示。由图7可以看出，启动压力梯度与地质强度指标关系显著。当GSI4658时，0；当GSI4658时，随GSI的增大而减小；当GSI4658时，随GSI的增大而增大。0039综上所述，只要在煤矿井下采集煤样或者直接观测煤壁，获取煤层的煤体结构，并与煤体结构量化的GSI图版进行比对，确定煤层的GSI值后，即可通过上述关系得出煤层瓦斯的启动压力梯度，将与瓦斯压力梯度P/L进行比较，即可快速判定煤层瓦斯的流动状态P/L，为低速非线性渗流；P/L，为扩散；0，为线性渗流。说明书CN103670392A1/6页8图1说明书附图CN103670392A2/6页9图2说明书附图CN103670392A3/6页10图3图4说明书附图CN103670392A104/6页11图5说明书附图CN103670392A115/6页12图6说明书附图CN103670392A126/6页13图7说明书附图CN103670392A13。