CN111982567A - 深孔反循环取样过程中瓦斯损失量补偿模型的构建方法 - Google Patents

深孔反循环取样过程中瓦斯损失量补偿模型的构建方法 Download PDF

Info

Publication number
CN111982567A
CN111982567A CN202010827687.1A CN202010827687A CN111982567A CN 111982567 A CN111982567 A CN 111982567A CN 202010827687 A CN202010827687 A CN 202010827687A CN 111982567 A CN111982567 A CN 111982567A
Authority
CN
China
Prior art keywords
gas
sampling
desorption
section
particle size
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010827687.1A
Other languages
English (en)
Other versions
CN111982567B (zh
Inventor
陈德敏
隆清明
赵旭生
吕贵春
周锦萱
邱飞
杨娟
张睿
张宪尚
饶家龙
胡杰
常宇
李柏均
刘娟
任文贤
唐勇
周燕
江旭
苏莉
兰祥云
曾银松
张森林
张卫东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CCTEG Chongqing Research Institute Co Ltd
Original Assignee
CCTEG Chongqing Research Institute Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by CCTEG Chongqing Research Institute Co Ltd filed Critical CCTEG Chongqing Research Institute Co Ltd
Priority to CN202010827687.1A priority Critical patent/CN111982567B/zh
Publication of CN111982567A publication Critical patent/CN111982567A/zh
Application granted granted Critical
Publication of CN111982567B publication Critical patent/CN111982567B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N1/00Sampling; Preparing specimens for investigation
    • G01N1/02Devices for withdrawing samples
    • G01N1/04Devices for withdrawing samples in the solid state, e.g. by cutting
    • G01N1/08Devices for withdrawing samples in the solid state, e.g. by cutting involving an extracting tool, e.g. core bit
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D21/00Measuring or testing not otherwise provided for
    • G01D21/02Measuring two or more variables by means not covered by a single other subclass
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
    • G01N15/02Investigating particle size or size distribution
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/15Correlation function computation including computation of convolution operations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N2015/0866Sorption

Abstract

本发明涉及一种深孔反循环取样过程中瓦斯损失量补偿模型的构建方法,属于煤炭开采及煤矿安全领域。该方法包括:S1:根据需要,确定取样位置和取样深度;S2:采用取样时所用钻头,以取样钻进时的钻进速度和钻机转速,在取样煤层施钻,收集煤屑并分析煤屑粒度分布;S3:获得取样管道内各个截面上的温度、压力和粒度分布;S4:建立不同温度、压力及粒度分布条件下的解吸曲线;S5:确定各个截面上瓦斯解吸量极值;S6:计算截面上瓦斯解吸损失量真实值;S7:数据拟合,确定整个取样过程中瓦斯解吸曲线;S8:计算取样过程中瓦斯解吸损失量。本发明使反循环取样过程中瓦斯损失量更加接近真实值,提高煤层瓦斯含量测量的准确性。

Description

深孔反循环取样过程中瓦斯损失量补偿模型的构建方法
技术领域
本发明属于煤炭开采及煤矿安全技术领域,涉及一种深孔反循环取样过程中瓦斯损失量补偿模型的构建方法。
背景技术
煤炭是我国能源的基石,长期以来,煤炭在我国一次性能源消费中占50%以上。但同时我国是世界上瓦斯灾害最严重的国家之一。煤层瓦斯含量是研究煤层瓦斯赋存规律,评估煤层气储量,预测煤与瓦斯突出危险性的核心指标。准确测量煤层瓦斯含量对确保煤矿安全高效生产和煤层气开发利用具有十分重要的意义。然而,“测不准”是目前煤层瓦斯含量测量中公认的难题。
煤层瓦斯含量井下直接测量技术是最常用的煤层瓦斯含量测量方法。煤层瓦斯含量直接测量误差根源是取样过程中的瓦斯损失量。国内外学者对取样过程中瓦斯损失量补偿计算进行了大量的研究,并提出了负指数式等多种计算模型,但这些模型都没有考虑煤样在取样过程中瓦斯解吸规律与煤样在常温常压下解吸规律的差异。现有的研究已经证实,温度、压力和粒度是影响瓦斯解吸逸散的主要因素,而反循环取样过程中,温度、压力并不恒定,煤样在反循环管道内输送过程中,碰撞破碎导致粒径发生显著改变。因此,亟需针对反循环取样的实际,建立更加准确的瓦斯损失量补偿模型,提高煤层瓦斯含量测量的准确性。
发明内容
有鉴于此,本发明的目的在于提供一种深孔反循环取样过程中瓦斯损失量补偿模型的构建方法,通过数学方法计算每个截面上的瓦斯损失量真实值,并通过数值拟合的手段获取取样过程中煤样瓦斯解吸曲线,从而实现了取样过程中瓦斯损失量的更加准确计算。
为达到上述目的,本发明提供如下技术方案:
一种深孔反循环取样过程中瓦斯损失量补偿模型的构建方法,具体包括以下步骤:
S1:根据需要,确定取样位置和取样深度;
S2:采用取样时所用钻头,以取样钻进时的钻进速度和钻机转速等参数,在取样煤层施钻,收集煤屑并分析煤屑粒度分布,从而获取原始粒度分布;
S3:获得取样管道内各个截面上的温度、压力和粒度分布;
S4:建立不同温度、压力及粒度分布条件下的解吸曲线;
S5:确定各个截面上瓦斯解吸量极值;
S6:计算截面上瓦斯解吸损失量真实值;
S7:数据拟合,确定整个取样过程中瓦斯解吸曲线;
S8:计算取样过程中瓦斯解吸损失量。
进一步,步骤S2中,在获取原始粒度分布时,钻头为取样时所用钻头,施钻位置为取样煤层,施钻参数与取样时施钻参数相同。
进一步,步骤S2中,煤屑粒度分布分析方法可采用筛分法或粒度分析仪分析。
进一步,步骤S3中,通过试验方法或数值模拟方法获得取样管道内各个截面上的温度、压力和粒度分布。
进一步,步骤S4中,所述解吸曲线由试验方法获取,依据截面K0、K1、K2、K3、……、Kn上的粒度分布以,以及温度压力值建立解吸曲线D0、D1、D2、D3、……、Dn
所述的解吸曲线D0和D1的坐标原点是t0;解吸曲线D2的坐标原点是D1与t1的交点,解吸曲线D3的坐标原点是D2与t2的交点;依次类推,解吸曲线Dn的坐标原点是Dn-1与tn-1的交点。
进一步,步骤S5中,所述各个截面上瓦斯解吸量极值包含极大值和极小值;煤样在反循环管道内由K0截面运移到K1截面时,瓦斯损失量极小值Q1min为解吸曲线D0与t1的交点,极大值Q1max为解吸曲线D1与t1的交点;在K2截面时,瓦斯损失量极小值Q2min为解吸曲线D1与t2的交点,极大值Q2max为解吸曲线D2与t2的交点;依次类推,在Kn截面时,瓦斯损失量极小值Qnmin为解吸曲线Dn-1与tn的交点,极大值Qnmax为解吸曲线Dn与tn的交点。
进一步,步骤S6中,截面上瓦斯解吸损失量真实值为:Qn=1/2(Qnmin+Qnmax)。
进一步,步骤S7中,取样过程中瓦斯解吸曲线D通过对各个截面上的瓦斯损失量真实值Q1、Q2、Q3、……、Qn通过拟合函数拟合所得。
进一步,步骤S7中,拟合函数采用:f(x)=eat,通过拟合确定参数a。
进一步,步骤S8中,取样过程中瓦斯解吸损失量计算公示为:
Figure BDA0002636820870000021
本发明的有益效果在于:本发明通过获取反循环取样管道各个截面上的煤样瓦斯解吸曲线,以及瓦斯损失量极大值和极小值,来计算各截面上的瓦斯损失量准确值,再通过数值拟合的方法获取反循环取样管道煤样瓦斯解吸曲线,最后通过瓦斯损失量解吸曲线来计算取样过程中的瓦斯损失量。本发明实现了变温变压以及变粒径条件下瓦斯损失量补偿计算,使反循环取样过程中瓦斯损失量更加接近真实值,提高了煤层瓦斯含量测量的准确性。
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述,并且在某种程度上,基于对下文的考察研究对本领域技术人员而言将是显而易见的,或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书来实现和获得。
附图说明
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作优选的详细描述,其中:
图1为本发明的方法原理图。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。需要说明的是,以下实施例中所提供的图示仅以示意方式说明本发明的基本构想,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。
请参阅图1,本发明优选了一种深孔反循环取样过程中瓦斯损失量补偿模型的构建方法,包括以下步骤:
S1:根据需要,确定取样位置和取样深度。
S2:采用取样时所用钻头,以取样钻进时的钻进速度和钻机转速等参数,在取样煤层施钻,收集煤屑后,分析煤屑粒度分布,从而获取原始粒度分布。
其中,钻头为取样时所用钻头,施钻位置为取样煤层,施钻参数与取样时施钻参数相同;在钻取煤样后,采用粒度分析仪分析煤样,获取原始粒度分布。
S3:通过试验方法或数值模拟方法获得取样管道内各个截面上的温度、压力和粒度分布。
S4:建立不同温度、压力及粒度分布条件下的解吸曲线。
将反循环取样管道按照一定的等分距离划分为K0、K1、K2、K3、……、Kn。截面K0上的粒度分布即为原始粒度分布,温度为取样地点的地层温度,压力值由试验或者数值模拟获取,K1、K2、K3、……、Kn截面上的温度由取样地点地层温度递减规律获取,压力值和粒度分布由试验或者数值模拟获取。
在获取各个截面上的温度、压力和粒度分布以后,通过温度、压力和粒度分布条件,由试验获取截面K0、K1、K2、K3……Kn对应的解吸曲线D0、D1、D2、D3……Dn。这里需要说明的是,解吸曲线D0和D1的坐标原点是t0。而解吸曲线D2的坐标原点是D1与t1的交点,解吸曲线D3的坐标原点是D2与t2的交点。依次类推,解吸曲线Dn的坐标原点是Dn-1与tn-1的交点。
S5:确定各个截面上瓦斯解吸量极值。
各个截面上瓦斯解吸量极值包含极大值和极小值。煤样在反循环管道内由K0截面运移到K1截面时,瓦斯损失量极小值Q1min为解吸曲线D0与t1的交点,极大值Q1max为解吸曲线D1与t1的交点。在K2截面时,瓦斯损失量极小值Q2min为解吸曲线D1与t2的交点,极大值Q2max为解吸曲线D2与t2的交点。依次类推,在Kn截面时,瓦斯损失量极小值Qnmin为解吸曲线Dn-1与tn的交点,极大值Qnmax为解吸曲线Dn与tn的交点。
S6:计算截面上瓦斯解吸损失量真实值。取截面Kn上瓦斯解吸损失量真实值Qn=1/2(Qnmin+Qnmax),同样的原理,获取截面K0、K1、K2、K3……Kn瓦斯解吸损失量真实值Q1、Q2、Q3……Qn
S7:数据拟合,确定整个取样过程中瓦斯解吸曲线。采用函数f(x)=eat拟合上述各截面上瓦斯损失量真实值Q1、Q2、Q3、……、Qn得到反循环取样过程中瓦斯解吸曲线和参数a的值。
S8:通过
Figure BDA0002636820870000041
来计算整个反循环取样过程中的瓦斯损失量。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (10)

1.一种深孔反循环取样过程中瓦斯损失量补偿模型的构建方法,其特征在于,该方法具体包括以下步骤:
S1:根据需要,确定取样位置和取样深度;
S2:采用取样时所用钻头,以取样钻进时的钻进速度和钻机转速,在取样煤层施钻,收集煤屑并分析煤屑粒度分布,从而获取原始粒度分布;
S3:获得取样管道内各个截面上的温度、压力和粒度分布;
S4:建立不同温度、压力及粒度分布条件下的解吸曲线;
S5:确定各个截面上瓦斯解吸量极值;
S6:计算截面上瓦斯解吸损失量真实值;
S7:数据拟合,确定整个取样过程中瓦斯解吸曲线;
S8:计算取样过程中瓦斯解吸损失量。
2.根据权利要求1所述的瓦斯损失量补偿模型的构建方法,其特征在于,步骤S2中,在获取原始粒度分布时,钻头为取样时所用钻头,施钻位置为取样煤层,施钻参数与取样时施钻参数相同。
3.根据权利要求1所述的瓦斯损失量补偿模型的构建方法,其特征在于,步骤S2中,煤屑粒度分布分析方法采用筛分法或粒度分析仪分析。
4.根据权利要求1所述的瓦斯损失量补偿模型的构建方法,其特征在于,步骤S3中,通过试验方法或数值模拟方法获得取样管道内各个截面上的温度、压力和粒度分布。
5.根据权利要求1所述的瓦斯损失量补偿模型的构建方法,其特征在于,步骤S4中,所述解吸曲线由试验方法获取,依据截面K0、K1、K2、K3、……、Kn上的粒度分布,以及温度压力值建立解吸曲线D0、D1、D2、D3、……、Dn
所述的解吸曲线D0和D1的坐标原点是t0;解吸曲线D2的坐标原点是D1与t1的交点,解吸曲线D3的坐标原点是D2与t2的交点;依次类推,解吸曲线Dn的坐标原点是Dn-1与tn-1的交点。
6.根据权利要求5所述的瓦斯损失量补偿模型的构建方法,其特征在于,步骤S5中,所述各个截面上瓦斯解吸量极值包含极大值和极小值;煤样在反循环管道的Kn-1截面运移到Kn截面时,瓦斯损失量极小值Qnmin为解吸曲线Dn-1与tn的交点,极大值Qnmax为解吸曲线Dn与tn的交点。
7.根据权利要求6所述的瓦斯损失量补偿模型的构建方法,其特征在于,步骤S6中,截面上瓦斯解吸损失量真实值为:Qn=1/2(Qnmin+Qnmax)。
8.根据权利要求7所述的瓦斯损失量补偿模型的构建方法,其特征在于,步骤S7中,取样过程中瓦斯解吸曲线D通过对各个截面上的瓦斯损失量真实值Q1、Q2、Q3、……、Qn通过拟合函数拟合所得。
9.根据权利要求8所述的瓦斯损失量补偿模型的构建方法,其特征在于,步骤S7中,拟合函数采用:f(x)=eat,通过拟合确定参数a。
10.根据权利要求9所述的瓦斯损失量补偿模型的构建方法,其特征在于,步骤S8中,取样过程中瓦斯解吸损失量计算公示为:
Figure FDA0002636820860000021
CN202010827687.1A 2020-08-17 2020-08-17 深孔反循环取样过程中瓦斯损失量补偿模型的构建方法 Active CN111982567B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010827687.1A CN111982567B (zh) 2020-08-17 2020-08-17 深孔反循环取样过程中瓦斯损失量补偿模型的构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010827687.1A CN111982567B (zh) 2020-08-17 2020-08-17 深孔反循环取样过程中瓦斯损失量补偿模型的构建方法

Publications (2)

Publication Number Publication Date
CN111982567A true CN111982567A (zh) 2020-11-24
CN111982567B CN111982567B (zh) 2024-02-02

Family

ID=73434488

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010827687.1A Active CN111982567B (zh) 2020-08-17 2020-08-17 深孔反循环取样过程中瓦斯损失量补偿模型的构建方法

Country Status (1)

Country Link
CN (1) CN111982567B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113049440A (zh) * 2021-03-24 2021-06-29 中国矿业大学 一种煤层瓦斯含量的井下直接测定方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7429358B1 (en) * 2002-05-20 2008-09-30 Hy-Energy, Llc Method and apparatus for measuring gas sorption and desorption properties of materials
CN102128765A (zh) * 2010-12-09 2011-07-20 中国矿业大学 一种煤矿井下直接快速测定煤层瓦斯含量方法
CN104453721A (zh) * 2014-12-17 2015-03-25 中煤科工集团重庆研究院有限公司 防塌孔反循环取样钻头
CN106093338A (zh) * 2016-07-13 2016-11-09 煤科集团沈阳研究院有限公司 井下反循环取样瓦斯解吸过程模拟测试装置及测试方法
CN106885755A (zh) * 2017-01-10 2017-06-23 中国矿业大学 一种煤矿井下快速测定煤层瓦斯参数的方法与装置
RU2632591C1 (ru) * 2016-04-18 2017-10-06 Лариса Николаевна Прокофьева Способ определения параметров массопереноса метана в угольном пласте
CN107290396A (zh) * 2017-04-06 2017-10-24 河南理工大学 煤岩吸附解吸一体化电阻率测量装置
US20200200004A1 (en) * 2017-10-12 2020-06-25 China University Of Mining And Technology Measurement-while-drilling method and device for assessing outburst risk of coal seam

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7429358B1 (en) * 2002-05-20 2008-09-30 Hy-Energy, Llc Method and apparatus for measuring gas sorption and desorption properties of materials
CN102128765A (zh) * 2010-12-09 2011-07-20 中国矿业大学 一种煤矿井下直接快速测定煤层瓦斯含量方法
CN104453721A (zh) * 2014-12-17 2015-03-25 中煤科工集团重庆研究院有限公司 防塌孔反循环取样钻头
RU2632591C1 (ru) * 2016-04-18 2017-10-06 Лариса Николаевна Прокофьева Способ определения параметров массопереноса метана в угольном пласте
CN106093338A (zh) * 2016-07-13 2016-11-09 煤科集团沈阳研究院有限公司 井下反循环取样瓦斯解吸过程模拟测试装置及测试方法
CN106885755A (zh) * 2017-01-10 2017-06-23 中国矿业大学 一种煤矿井下快速测定煤层瓦斯参数的方法与装置
CN107290396A (zh) * 2017-04-06 2017-10-24 河南理工大学 煤岩吸附解吸一体化电阻率测量装置
US20200200004A1 (en) * 2017-10-12 2020-06-25 China University Of Mining And Technology Measurement-while-drilling method and device for assessing outburst risk of coal seam

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张向阳;郭孟志;宋传杨;高留强;: "解吸法测定煤层瓦斯含量过程中瓦斯损失量3种推算方法对比分析", 煤矿安全, no. 08 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113049440A (zh) * 2021-03-24 2021-06-29 中国矿业大学 一种煤层瓦斯含量的井下直接测定方法
CN113049440B (zh) * 2021-03-24 2022-03-25 中国矿业大学 一种煤层瓦斯含量的井下直接测定方法

Also Published As

Publication number Publication date
CN111982567B (zh) 2024-02-02

Similar Documents

Publication Publication Date Title
CN103334739B (zh) 一种测定煤层瓦斯压力的方法及装置
CN113820078B (zh) 一种示踪法与稀释浓度法联合检测堤坝渗漏的方法
CN102979579A (zh) 实时分析煤与瓦斯突出危险性的方法
CN104695950A (zh) 火山岩油藏产能预测方法
CN109740119B (zh) 一种tbm掘进隧洞围岩单轴抗压强度快速估算方法
CN104636532A (zh) 一种煤矿瓦斯抽采钻孔封孔深度和长度的确定方法
CN110069751B (zh) 一种非欠压实地层孔隙压力预测方法
CN108825226A (zh) 一种采用化学示踪剂评估压后产气量的方法及装置
CN111982567A (zh) 深孔反循环取样过程中瓦斯损失量补偿模型的构建方法
CN110529106B (zh) 一种利用测井资料确定煤层显微组分含量的方法
CN102900432A (zh) 利用微钻时数据随钻计算录井孔隙度进行储层评价的方法
CN102996114A (zh) 一种测井资料的斜井校正方法
CN111950111B (zh) 一种适用于底部开放的碳酸盐岩储层动态分析方法
CN111155980B (zh) 一种水流优势通道识别方法及装置
CN110706344A (zh) 一种碳酸盐岩缝洞储层连通性的3d建模方法
CN205400721U (zh) 一种井口气远程测定装置
CN111206923B (zh) 一种利用钻能确定节理岩体模量比与强度比的测试方法
CN114483043A (zh) 一种用于隧道掘进的探测控制系统及其探测方法
CN110486005B (zh) 一种用于识别页岩气井甜点层的方法
CN114370267B (zh) 一种高温地热田钻进过程中钻遇热储温度的实时计算方法
CN110909311A (zh) 一种计算薄煤层含气量的方法
CN117365642A (zh) 一种基于煤层瓦斯含量的煤层瓦斯抽采半径的测定方法
CN117288587B (zh) 岩体抗拉强度随钻测试方法与系统
CN113627700B (zh) 一种气井生产效果的评价方法及装置
CN110991081B (zh) 一种基于地面钻孔抽采瓦斯确定邻近层抽采瓦斯量的方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant