磨料空气射流破煤冲蚀模型研究.pdf

第4 5 卷第5 期 2 0 2 0 年5 月 煤炭学报 J O U R N A LO FC H I N AC O A LS O C I E T Y V 0 1 ．4 5N o ．5 M a y 2 0 2 0 移动阅读 刘勇，李志飞，魏建平，等．磨料空气射流破煤冲蚀模型研究[ J ] ．煤炭学报，2 0 2 0 ，4 5 5 1 7 3 3 1 7 4 2 ．d o i i 0 ． 1 3 2 2 5 ／j ．c n k i ．j C C S ．D Y 2 0 ．0 2 8 2 L I UY o n g ，L IZ h i f e i ，W E IJ i a n p i n g ，e ta 1 ．E r o s i o nm o d e lo fa b r a s i v ea i rj e tu s e di nc o a lb r e a k i n g [ J ] ．J o u m a lo fC h i n a C o a lS o c i e t y ，2 0 2 0 ，4 5 5 1 7 3 3 1 7 4 2 ．d o i 1 0 ．1 3 2 2 5 ／j ．c n k i ．j C C ．D Y 2 0 ．0 2 8 2 磨料空气射流破煤冲蚀模型研究 刘勇1 ”，李志飞1 ”，魏建平1 ’2 ，张慧栋1 ’2 1 ．河南理工大学瓦斯地质与瓦斯治理国家重点实验室培育基地，河南焦作4 5 4 0 0 0 ；2 ．煤炭安全生产与清洁高效利用省部共建协同创新中心 河南焦作4 5 4 0 0 0 摘要进行预抽煤矿瓦斯时，采用水力化增透措施容易出现抱钻、塌孔等现象，降低瓦斯抽采率。 为此，提出磨料空气射流破煤增透措施，提高钻孔利用率。为明确磨料空气射流破煤参数，采用 A N S Y SF l u e n t 软件中的离散相模型对磨料粒子在喷嘴中的加速效果进行了计算，分析了空气压力、 磨料密度和磨料粒径对磨料加速的影响，考虑了反弹磨料二次冲蚀作用，建立了磨料空气射流冲蚀 能量方程；基于有限元和光滑粒子模型，采用L S D Y N A 对磨料空气射流冲蚀煤岩体进行数值模 拟，分析了射流扩散角和磨料的形状特性对煤岩体冲蚀的影响规律，建立了磨料空气射流冲蚀能量 转化率方程；采用高压磨料空气射流实验系统，使用 L a v a l 喷嘴，在不同气体压力条件选用标准的 1 8 0I x m 石榴石磨料对煤体、砂岩、灰岩和花岗岩进行磨料空气射流冲蚀实验，获取了单位能量破碎 体积参数，建立了磨料空气射流冲蚀模型；采用控制变量法，进行磨料空气射流冲蚀灰岩的实验，利 用不同影响因素对冲蚀效果的影响实验数据对冲蚀模型进行了修正；进行磨料空气射流冲击煤体 实验，验证了冲蚀模型的准确性。对磨料空气射流冲蚀煤体体积模型进行了分析，结果表明影响冲 蚀效果的因素权重依次为气体压力 扩散角 磨料密度 磨料粒径 磨料形状；磨料空气射流破煤 的最优射流参数为1 8 0 m 石榴石、气体压力1 5M P a 、扩散角1 0 ．5 2 0 。在进行模型建立过程中，未 将受压条件放入模型中，今后将会进行相关实验，把受压条件作为修正系数，进一步完善该冲蚀模 型。 关键词磨料空气射流；冲蚀模型；二次冲蚀；卸压增透；瓦斯抽采 中图分类号T D 7 1 2 文献标志码A文章编号0 2 5 3 - 9 9 9 3 2 0 2 0 0 5 1 7 3 3 - 1 0 E r o s i o nm o d e lo fa b r a s i v ea i rj e tu s e di nc o a lb r e a k i n g L I UY o n 9 1 ．- ，L IZ h i f e i l ”，W E IJ i a n p i n 9 1 ⋯，Z H A N GH u i d o n 9 1 ’2 1 - S t a t eK e yL a b o r a t o r yC u l t i v a t i o nB a s ef o rG a sG e o l o g ya n dG a sC o n t r o l ，H e n a nP o l y t e c h n i cU n i v e r s i t y ，J i a o z u o4 5 4 0 0 0 ，C h i n a ；2 ．S t a t eC o U a b o r a t i v eI n n o v a t i o nC e n t e ro fC o a lW o r kS a f e t ya n dC l e a n e f f i c i e n c yU t i l i z a t i o n ，J i a o z u o4 5 4 0 0 0 ，C h i n a A b s t r a c t I n d u c e db yt h eh y d r a u l i ci m p r o v i n gp e r m e a b i l i t yt e c h n o l o g yi nc o a ls e a m ，t h eb o r e h o l ec o l l a p s e se a s i l y ．T h i s p a p e rd e v e l o p e dan e wt e c h n o l o g yo fi m p r o v i n gp e r m e a b i l i t yw i t ha b r a s i v ea i rj e t ．T oc l a r i f yt h ec o a lb r e a k i n gp a r a m e - t e r so fa b r a s i v ea i rj e t ，u s i n gA N S Y SF l u e n ts o f t w a r ei nad i s c r e t ep h a s em o d e l D P M ，t h ea c c e l e r a t i n ge f f e c to fa b r a s i v ep a r t i c l ei nt h en o z z l ew a sc a l c u l a t e d ，a n dt h ei n f l u e n c eo fa i rp r e s s u r e ，a b r a s i v ed e n s i t ya n dp a r t i c l es i z eo nt h ea c c e l e r a t i o no fa b r a s i v ew a sa n a l y z e d ．B a s e do nt h ef i n i t ee l e m e n t F E M a n ds m o o t hp a r t i c l em e t h o d S P H ，t h eL S D Y N As o f t w a r ew a su s e dt os i m u l a t et h ea b r a s i v eg a sj e te r o s i o nt a r g e tu n d e rc e r t a i ne n e r g yc o n d i t i o n s ．U s i n gt h ee x - 收稿日期2 0 2 0 - 0 2 - 2 6修回日期2 0 2 0 - 0 4 2 9责任编辑郭晓炜 基金项目国家重点研发计划资助项目 2 0 1 7 Y F C 0 8 0 4 2 0 7 ；国家自然科学基金资助项目 5 1 7 0 4 0 9 6 ，5 1 9 7 4 1 0 9 作者简介刘勇 1 9 8 4 一 ，男，山东l 临沂人，副教授，博士。E - m a i l y o o n l i u h p u ．e d u ．e n 通讯作者魏建平 1 9 7 l 一 ，男，河南驻马店人，教授，博士。E - m a i l w e i j i a n p i n g h p u ．e d u ．e n 万方数据 煤炭 学报 2 0 2 0 年第4 5 卷 p e r i m e n ts y s t e mw i t ha i rh i g hp r e s s u r ea b r a s i v ej e t ，a n dt h eL a v a ln o z z l e ，u n d e rt h ec o n d i t i o no fd i f f e r e n tg a sp r e s - s u r e s ，t h ea b r a s i v ej e ta i re r o s i o ne x p e r i m e n t so nc o a l ，s a n d s t o n e ，l i m e s t o n ea n dg r a n i t ew e r ec o n d u c t e du s i n gs t a n d a r d 18 0m i c r o n sg a r n e ta b r a s i v e ，a n do b t a i n e dt h ep a r a m e t e r so fc r u s h i n ge n e r g yp e ru n i tv o l u m e ，t h e nt h ee q u a t i o no fa b - r a s i v ea i rj e te r o s i o ne n e r g yc o n v e r s i o nr a t ew a se s t a b l i s h e d ．U s i n gt h ec o n t r o lv a r i a b l em e t h o d ，t h ei n f l u e n c eo fd i f f e r - e n ti n f l u e n c i n gf a c t o r so nt h ee r o s i o ne f f e c tw a sc o r r e c t e db yt h ee x p e r i m e n t a ld a t ao fa b r a s i v ea i ri e te r o s i o no fl i m e ． s t o n e ．T h ea b r a s i v ea i rj e ti m p a c tc o a le x p e r i m e n tw a sc a r r i e do u tt ov e r i f yt h ea c c u r a c yo ft h ee r o s i o nm o d e l ．T h ev o l u m e t r i cm o d e lo fa b r a s i v ea i rj e te r o s i o nm e d i aw a sa n a l y z e d ．T h er e s u l t ss h o wt h a tt h ew e i g h t so ft h ef a c t o r sa f f e c t i n g t h ee r o s i o ne f f e c ta r e a i rp r e s s u r e d i f f u s i o na n g l e a b r a s i v ed e n s i t y a b r a s i v ep a r t i c l es i z e a b r a s i v es h a p e ．T h eo p t i - m u m j e tp a r a m e t e r so ft h ea b r a s i v ea i rj e tb r e a k i n gc o a la r e18 0 斗mg a r n e t ，a i rp r e s s u r e15M P a ，a n dd i f f u s i o na n g l e 1 0 ．5 2 。．I nt h ep r o c e s so fm o d e lb u i l d i n g ，t h ec o m p r e s s i o nc o n d i t i o n sa r en o tp u ti n t ot h em o d e l ．I nt h ef u t u r e ，s o m e r e l e v a n te x p e r i m e n t sw i l lb ec a r r i e do u tt ot a k et h ec o m p r e s s i o nc o n d i t i o n sa st h ec o r r e c t i o nc o e f f i c i e n tt of u r t h e ri m p r o v et h ee r o s i o nm o d e l ． K e yw o r d s a b r a s i v ea i rj e t ；e r o s i o nm o d e l ；s e c o n d a r ye r o s i o n ；p r e s s u r er e l i e ff o ri m p r o v i n gp e r m e a b i l i t y ；g a sd r a i n a g e 卸压抽采是强化瓦斯抽采和防治瓦斯灾害的重 要手段，其中水力化卸压技术应用最为广泛，如高压 水射流割缝和水力冲孔等，均取得了较好的应用效 果。2J 。但上述技术在松软煤层中应用时，经常出现 塌孔等现象，严重限制了在松软煤层中的推广应用。 除此之外，煤体的低孔、低渗的特点，使水锁效应尤为 突出【3 。4J 。水的侵入使煤体孔隙含水饱和度升高，有 效应力发生改变，煤体渗透率大幅度降低。尤其在俯 孔中，孔底积水容易造成煤体泥化，堵塞瓦斯抽采通 道”J 。随着开采深度的不断增加，高地应力和水的 耦合作用使上述问题更为严重。 采用气体作为卸压抽采的动力，可避免水力化卸 压技术的缺点，是对水力卸压增透的重要补充。如液 态二氧化碳相变致裂技术【6J 、超临界二氧化碳射流 破煤岩技术。列和高压空气爆破冲击等∽J 。基于气力 化措施的思想，笔者提出采用磨料空气射流破煤扩孔 提高瓦斯抽采率。由于磨料空气射流能够以较低压 力使射流达到超声速，使磨料充分加速，其破煤效率 较高，甚至优于水射流破煤效率一。10 | 。磨料空气射流 破煤深度是射流冲击应力波效应和冲蚀磨损的综合 作用结果。应力波效应能够在冲蚀破煤的基础上进 一步扩大破煤范围1 l 。磨料空气射流的冲蚀磨损是 破煤的基础，其本质为固相颗粒冲蚀磨损造成煤岩体 的冲蚀破坏。冲蚀磨损效率决定于磨料冲击动能和 能量转化率。磨料的加速段为磨料提供冲击动能，加 速效果主要受到气体压力，磨料密度和粒径的影 响【1 2 1 ；能量转化率与磨料粒径、形状、密度、扩散角和 靶体力学性质等有关3 | 。目前关于冲蚀机理和冲蚀 模型的研究多针对单颗磨料粒子对塑性材料冲蚀磨 损。而高压磨料空气射流以射流束的形式作用于煤 体，磨料冲击煤体后，发生反弹，反弹磨料对煤体“二 次冲蚀”使冲蚀坑的体积进一步增大4 ‘16 。。且煤体 具有典型的脆性特征，目前的冲蚀模型不能准确计算 冲蚀参数。 为准确计算磨料空气射流扩孔参数，需建立以粒 子束冲击为基础，考虑“二次冲蚀”作用以及气体压 力、磨料粒径、磨料密度、扩散角、磨料形状以及脆性 材料力学性质影响的冲蚀模型。为此，笔者基于L a - v a l 喷嘴，提取数值模拟结果中磨料空气射流冲蚀过 程中各因素的作用效果数据，结合理论推导建立磨料 空气射流冲蚀模型，通过磨料空气射流冲蚀实验计算 模型参数，对模型进行修正。 1 磨料空气射流冲蚀模型的建立 磨料空气射流主要以粒子冲蚀破碎煤体，根据其 作用机理，冲蚀模型的建立主要从磨料空气射流冲蚀 能量方程和能量转化率方程两方面考虑。 1 ．1 磨料空气射流冲蚀能量方程 磨料空气射流对煤体的破碎能量主要来源于磨 料粒子的冲击动能，而冲击动能取决于磨料速度、磨 料密度和磨料粒径。在磨料空气射流中，不同的气体 压力以及磨料的密度和粒径会产生不同的磨料速度。 A N S Y SF l u e n t 数值模拟软件中的离散相模型 D P M 可以较为准确的对两相流中粒子加速进行计算川， 因此本文使用F l u e n t 对4 种磨料在不同粒径 7 5 ， 1 2 5 ，1 8 0 ，4 1 0 和不同压力 5 ，1 0 ，1 5 ，2 0 ，2 5 条件下 磨料粒子的加速效果进行了计算。 数值模拟中使用的喷嘴结构如图1 所示，数值模 拟中喷嘴出口压力0 ．1M P a ，温度均为3 0 0K ，磨料质 量流量为0 ．0 1 6k g ／s 。数值模拟结果中部分磨料粒 子的速度分布云图如图2 所示。 根据动能定理和磨料能量E 与参数之间的耦合 万方数据 第5 期刘勇等磨料空气射流破煤冲蚀模型研究 1 7 3 5 图1喷嘴几何形状及尺寸 F i g ．1 S t r u c t u r eo fL a v a ln o z z l e ， 管 g 世 瑙 蒌 墩 f E 越 删 婪 恤 图2 石榴石磨料速度分布 F i g ．2 A b r a s i v ev e l o c i t yd i s t r i b u t i o nd i a g r a mo fg a r n e t ， 譬 g 赵 删 婪 鼬 关系，即得到以下的数学表达式 E ．厂 d ，P ，P 1 式中，P 为气体压力，M P a ；J D 为磨料密度，k g ／m 3 ；d 为 磨料粒径，I l l 。 分别提取数值模拟结果中喷嘴出口出轴线位置 磨料速度，以磨料速度M m ／s 为因变量，以气体压 力、磨料密度、磨料粒径为自变量分析磨料粒子在不 同压力条件下对冲蚀速度的影响，如图3 a 所示。 通过曲线分析得到，气体压力与磨料速度呈幂函数关 系，相关性系数R 2 0 ．9 9 8 ，气体压力单因素与磨料 速度回归方程为 M 后l P o 3 1 6 2 式中，矗．为受磨料粒径和密度影响的常数。 分别对数值模拟结果中气体压力为2 0M P a 时 不同粒径和密度的影响进行分析，得到在一定磨料粒 磨料密度／ k g m 。3 1 b ， 譬 吕 蜊 删 婪 墩 图3 气体压力、磨料密度和粒径与磨料速度的关系 F i g ．3R e l a t i o n s h i pb e t w e e ng a sp r e s s u r e ．a b r a s i v ep a r t i c a ls i z ea n da b r a s i v ev e l o c i t y 径、密度条件下磨料速度和压力的关系如式 3 所 示。再对图3 b 中粒径分别为7 5 ，1 8 0 ，4 1 0 斗n 条件 下不同密度对速度的影响通过曲线回归的方法分析 得到2 0M P a 下粒径1 8 0 m 磨料密度单因素与磨料 速度的回归方程 式 4 。最后对图3 c 中磨料密 度为26 6 0 ，35 0 0 ，39 5 0k g ／m 3 时不同粒径对速度影 响的分析得到得到2 0M P a 下密度35 0 0k g ／m 3 磨料 粒径单因素与速度的回归方程 式 5 。 坠坐竺f 当。3 1 6r 3 、 配2 0 ，p ，d 七1 2 0 0 3 1 6 1 2 0 ／ 、。7 式中，“％，。为气体压力P 、密度P 、粒径d 下磨料速度 值；“0 D ．，，为气体压力2 0M P a 、密度P 、粒径d 下磨料 速度。 M 一2 ．2 9 1 0 5 p 一35 0 0 2 3 3 0 ．0 2 4 u 一1 ．8 4 1 0 3 d 1 0 6 9 0 2 3 2 8 ．8 5 分析式 4 和式 5 得到 p - 35 0 0 2 ， d 1 0 6 9 0 2 与M 呈一次线性关系，对M 进行二元线性回归 分析，拟定的回归方程为 M b o b lA 6 2 肛 A P 一35 0 0 2 6 肛 d 1 0 6 9 0 2 式中，b 。，b 。，b 为常数。 求解式 6 的二元线性回归方程，得到R 2 0 ．9 9 6 ，说明 p - 35 0 0 2 和 d x l 0 6 9 0 2 对磨料速度 M 具有显著影响，并具有线性相关性，得到了2 0M P a 下关于d 和P 与M 的二元非线性回归方程 M [ 一2 ．2 9 1 0 5 p 一35 0 0 2 1 ．8 4 1 0 3 d 1 0 6 9 0 2 3 4 5 ] 7 联合式 3 ， 7 得到气体压力、磨料粒径、磨料 密度与磨料速度的三元回归方程为 ，D 、03 1 6 E 2 ㈥ [ - 2 2 9 1 0 ‘5 p 一35 0 0 2 1 ．8 4 1 0 3 d 1 0 6 9 0 2 3 4 5 ] 8 根据动能定理，结合磨料速度得到单颗磨料粒子 能量E 。 J ，并推算出磨料空气射流单位时间冲击动 能E J ／s 与磨料密度、磨料粒径、气体压力的关系如 式 1 1 所示。 O 0 O O O 0 O O 3 8 2 7 2 6 ●6 j ，一 钙∞弘引卯勉懈B∞邪0 万方数据 1 7 3 6 煤炭学报 2 0 2 0 年第4 5 卷 E 。， m 。，M 2 9 E E 。．一m 1 0 1 玎Z E 鲁 嘉 。6 3 2 [ 一2 ．2 9 - 。一5 c p 一35 。。，2 1 ．8 4 1 0 - 3 d 1 0 6 9 0 2 3 4 5 ] 11 式中，m 。为单颗磨料粒子质量，k g ；南为磨料的质量 流量，k g ／s 。 射流能量的数学表达式 1 1 中，磨料空气射流 冲蚀效果主要与磨料速度呈2 次方关系，这与F i n n i e 提出的微切削理论中8 ] 磨料入射速度与靶体冲蚀体 积之间的关系为平方关系一致。 1 ．2 磨料空气射流能量转化率方程 磨料空气射流粒子束的冲蚀效果受到射流扩散 角、磨料形状的影响，其中磨料的形状使用球形度表 示，扩散角和磨料球形度影响了磨料空气射流能量分 布以及利用率，即 r ／ 八p ，妒 1 2 式中，叼为磨料空气射流能量的利用率；0 为磨料空 气射流扩散角；妒为磨料球形度。 球形度的计算式为 ∞．掣 1 3 妒2 ‘i 二二二 L， , ／6 、，讯 式中，5 为磨料的表面积；V 为磨料的体积。 能量利用率使用式 1 4 计算，综合能量转化率 如式 1 5 所示。 5 9 ， ∞ 9 4 ； ．3 0 刨 删 ．6 5 7 7 ， 1 i _ 二 1 4 D 1 s f l - ∑叼， 1 5 ，Z i 0 式中，5 ，为磨料的综合能量转化率，其值越高，说明 粒子作用在冲蚀岩石的能量比例越大；凡为所做的差 值数量，0 磨料粒径 磨料形状特性。 通过对冲蚀率方程式 2 7 进行求解计算，得到 当选用磨料密度为35 0 0k g ／m 3 、磨料粒径1 8 0 斗m 、 扩散角1 0 ．9 2 0 、磨料球形度1 ．9 2 时磨料空气射流冲 蚀率达到最大；相较于常用磨料，采用1 8 0 斗m 石榴 石磨料其冲蚀效果比其他磨料冲蚀效果较优，对灰岩 的冲蚀率可达到0 ．5 0 25c m 3 ／s 。 使用修正后的冲蚀模型与试验结果进行对比如 图1 4 所示。 2 ．5 冲蚀煤体试验验证 使用4 , 5 0m m x1 0 0m m 的煤体试样进行磨料气 体射控流冲蚀煤体试验，试验磨料质量流量 0 ．0 1 6k g ／s 、扩散角1 0 ．5 2 0 ，冲蚀时间2 0S ，气体压力 选用1 0M P a ，磨料使用经过测量试验粒径为1 8 0 仙m 的石英砂、棕刚玉、碳化硅和石榴石，试验结果如图 1 5 所示。 c 图1 5 冲蚀煤样试验效果 F i g ．1 5E x p e r i m e n t a le f f e c to fe r o s i o nc o a ls a m p l e 冲蚀试验中，磨料空气射流对煤体试样较高的破 坏程度导致冲蚀效果数据的采集较为困难，笔者通过 冲蚀残留煤样的冲蚀坑进行还原，分别以“V ”字形冲 蚀坑计算其冲蚀最小体积，如图1 5 a 所示，以及冲 蚀残留的冲蚀坑体积加上煤样损失体积为冲蚀最大 体积，如图1 5 b 所示，并结合冲蚀模型的计算数值 进行比较。通过测量煤体冲蚀深度分别为4 ．4 1 ， 5 ．6 0 ，5 ．6 8 ，6 ．1 2c m ，煤样碎裂损失长度为3 ．2 1 ， 2 ．9 5 ，3 ．5 0 ，3 ．3 6c m ，计算得到试验冲蚀体积和冲蚀 模型计算体积见表3 。 表3 试验结果与冲蚀模型计算结果 T a b l e3 E x p e r i m e n t a lr e s u l t sa n dc a l c u l a t i o nr e s u l t so f e r o s i o nm o d e I 通过试验结果对比可以看出，计算的冲蚀体积与 试验冲蚀最小体积相接近，其误差在1 0 ％以内。从 试验煤样破坏形态来看，冲蚀最小体积的倒锥形符合 冲蚀理论特征，为磨料颗粒切削破坏，而冲蚀破坏造 成的最大体积与最小体积的差值部分为断裂破坏，其 破坏形式符合应力波破煤特征l ’驯，试验结果表明 本文构建的磨料空气射流冲蚀煤体模型可以较好的 预测磨料空气射流对煤体的冲蚀体积。 3结论 1 为明确磨料空气射流扩孔破煤参数，建立了 磨料空气射流冲蚀体积计算模型。该模型考虑了粒 子束、磨料形状、磨料密度、射流扩散角和煤体力学性 质对冲蚀效果的影响。模型考虑了能量转化率，较好 的体现二次冲蚀作用以及磨料在冲击过程中的破碎 现象。 2 通过对冲蚀体积模型的计算和冲蚀试验发 现磨料空气射流的影响程度由大到小分别为气体压 力 扩散角 磨料密度 磨料粒径 磨料形状。 3 通过磨料空气射流冲蚀煤体体积计算模型 计算得出磨料空气射流冲蚀煤体的最佳射流参数为 1 8 0 恤m 石榴石、气体压力1 5M P a 、扩散角1 0 ．5 2 。。 参考文献 R e f e r e n c e s [ 1 ] 袁亮，薛俊华，张农，等．煤层气抽采和煤与瓦斯共采关键技术 现状与展望[ J ] ．煤炭科学技术，2 0 1 3 ，4 1 9 6 1 1 ，1 7 ． Y U A NL i a n g ，X U EJ u n h u a ，Z H A N GN o n g ，e ta 1 ．D e v e l o p m e n to r i e n t a t i o na n ds t a t u so fk e yt e c h n o l o g yf o rm i n eu n d e r g r o u n dc o a l b e dm e t h a n ed r a i n a g ea sw a l la sc o a la n dg a ss i m u l t a n e o u sm i n i n g [ J ] ．C o a lS c i e n c ea n dT e c h n o l o g y ，2 0 1 3 ，4 1 9 6 一1 1 ，1 7 ． [ 2 ] 徐冬冬．煤矿井下水力冲孔卸压抽采技术现状及展望[ J ] ．煤炭 技术，2 0 1 6 ，3 5 2 1 7 6 - 1 7 8 ． X UD o n g d o n g ．R e s e a r c hc u r r e n ts i t u a t i o na n dp r o s p e c to fh y d r a u l i c f l u s h i n gt e c h n o l o g y f o r p r e s s u r er e l e a s i n ga n dg a sd r a i n a g e i n 万方数据 1 7 4 2 煤炭 学报 2 0 2 0 年第4 5 卷 c o a lm i n eu n d e r g r o u n d [ J ] ．C o a lT e c h n o l o g y ，2 0 1 6 ，3 5 2 1 7 6 一 1 7 8 ． [ 3 ] H O S S E I NJ ，F A R U KC ．D e t e r m i n a t i o no fm u l t i c o m p o n e n tg a sa n d w a t e re q u i l i b r i u ma n dn o n e q u i l i b r i u ms o r p t i o ni s o t h e r m si nc a r b o n a c e o u ss o l i d sf r o m e a r l y t i m em e a s u r e m e n t s [ J ] ．F u e l ，2 0 0 7 ， 8 6 1 0 ／1 1 1 6 0 1 1 6 1 3 ， [ 4 ] L A W R E N C EO A ．G r o u n d w a t e rf l o wa s s o c i a t e dw i t hc n a l b e dg a sp r o - d u c t i o n ，F e r r o ns a n d s t o n e ，e a s t c e n t r a lU t a h [ J ] ．I n t e r n a t i o n a lJ o u r - n a lo fC o a lG e o l o g y ，2 0 0 3 ，5 6 1 ／2 6 9 9 5 ． [ 5 ] 魏建平，秦恒洁，王登科．基于水分影响的加一卸载围压条件下 含瓦斯煤渗流特性研究[ J ] ．采矿与安全工程学报，2 0 1 4 ， 3 l 6 9 8 7 9 9 4 ． W E IJ i a n p i n g ，Q I NH e n g j i e ，W A N GD e n g k e ．S e e p a g ec h a r a c t e r i s t i c s o fc o a lm a s sc o n t a i n i n gg a sc o n s i d e r i n gm o i s t u r ee f f e c ti nl o a d i n g - ．u n ．． 1 0 a d i n gc o n f i n i n gp r e s s u r et e s t [ J ] ．J o u r n a lo fM i n i n g ＆S a f e t yE n g i n e e r i n g ，2 0 1 4 ，3 l 6 9 8 7 9 9 4 ． [ 6 ] 王兆丰，周大超，李豪君，等．液态C O 相变致裂二次增透技术 [ J ] ．河南理工大学学报 自然科学版 ，2 0 1 6 ，3 5 5 5 9 7 - 6 0 0 ． W A N GZ h a o f o n g ，Z H O UD a c h a o ，L IH a o j u n ，e ta 1 ．R e s e a r c ho ns e c - o n d a r yf r a c t u r i n gt e c h n o l o g yb yp h a s et r a n s f o r m a t i o no fl i q u i dC 0 2 [ J ] ．J o u r n a lo fH e n a nP o l y t e c h n i cU n i v e r s i t y N a t u r a lS c i e n c e ， 2 0 1 6 ，3 5 5 5 9 7 6 0 0 ． [ 7 ] 王磊，梁卫国．超临界C 0 2 压裂下煤岩体裂缝扩展规律试验研 究[ J ] ．煤炭科学技术，2 0 1 9 ，4 7 2 6 5 7 0 ． W A N GL e i ，L I A N GW e i g u o ．E x p e r i m e n t a ls t u d yo nc r a c kp r o p a g a l i o no fc o a l r o c kn l a s u n d e r s u p e r c r i t i c a lC 0 2f r a c t u r i n g [ J ] ． C o a lS c i e n c ea n dT e c h n o l o g y ，2 0 1 9 ，4 7 2 6 5 7 0 ． [ 8 ] 李守国，吕进国，贾宝山，等．高压空气爆破低透气性煤层增透 技术应用研究[ J ] ．中国安全科学学报，2 0 1 6 ，2 6 4 i 1 9 - 1 2 5 ． L IS h o u g u o ，L UJ i n g u o ，J I AB a o s h a n ，e ta 1 ．A n t i r e f l e c t i o nt e c h n o l o g y a p p l i c a t i o ns t u d yo fh i g hp r e s s u r ea i rb l a s t i n gl o wp e r m e a b i l i t yc o a l s e a m [ J ] ．C h i n aS a f e t yS c i e n c eJ o u r n a l ，2 0 1 6 ，2 6 4 1 1 9 1 2 5 ． [ 9 ] 刘