CN111612258A - Method for judging and identifying gas abnormality by using gas desorption quantity characteristics - Google Patents

Method for judging and identifying gas abnormality by using gas desorption quantity characteristics Download PDF

Info

Publication number
CN111612258A
CN111612258A CN202010457999.8A CN202010457999A CN111612258A CN 111612258 A CN111612258 A CN 111612258A CN 202010457999 A CN202010457999 A CN 202010457999A CN 111612258 A CN111612258 A CN 111612258A
Authority
CN
China
Prior art keywords
gas
ave
shift
gas emission
jth
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
CN202010457999.8A
Other languages
Chinese (zh)
Other versions
CN111612258B (en
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 CN202010457999.8A priority Critical patent/CN111612258B/en
Publication of CN111612258A publication Critical patent/CN111612258A/en
Application granted granted Critical
Publication of CN111612258B publication Critical patent/CN111612258B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Forestry; Mining
    • GPHYSICS
    • G08SIGNALLING
    • G08BSIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B21/00Alarms responsive to a single specified undesired or abnormal condition and not otherwise provided for
    • G08B21/02Alarms for ensuring the safety of persons
    • G08B21/12Alarms for ensuring the safety of persons responsive to undesired emission of substances, e.g. pollution alarms
    • G08B21/14Toxic gas alarms
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16YINFORMATION AND COMMUNICATION TECHNOLOGY SPECIALLY ADAPTED FOR THE INTERNET OF THINGS [IoT]
    • G16Y10/00Economic sectors
    • G16Y10/20Mining
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16YINFORMATION AND COMMUNICATION TECHNOLOGY SPECIALLY ADAPTED FOR THE INTERNET OF THINGS [IoT]
    • G16Y20/00Information sensed or collected by the things
    • G16Y20/10Information sensed or collected by the things relating to the environment, e.g. temperature; relating to location
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16YINFORMATION AND COMMUNICATION TECHNOLOGY SPECIALLY ADAPTED FOR THE INTERNET OF THINGS [IoT]
    • G16Y40/00IoT characterised by the purpose of the information processing
    • G16Y40/10Detection; Monitoring
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16YINFORMATION AND COMMUNICATION TECHNOLOGY SPECIALLY ADAPTED FOR THE INTERNET OF THINGS [IoT]
    • G16Y40/00IoT characterised by the purpose of the information processing
    • G16Y40/20Analytics; Diagnosis
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16YINFORMATION AND COMMUNICATION TECHNOLOGY SPECIALLY ADAPTED FOR THE INTERNET OF THINGS [IoT]
    • G16Y40/00IoT characterised by the purpose of the information processing
    • G16Y40/50Safety; Security of things, users, data or systems
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A50/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE in human health protection, e.g. against extreme weather
    • Y02A50/20Air quality improvement or preservation, e.g. vehicle emission control or emission reduction by using catalytic converters

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • Computing Systems (AREA)
  • General Health & Medical Sciences (AREA)
  • Economics (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Toxicology (AREA)
  • Development Economics (AREA)
  • Theoretical Computer Science (AREA)
  • Marketing (AREA)
  • Tourism & Hospitality (AREA)
  • Environmental & Geological Engineering (AREA)
  • Game Theory and Decision Science (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Biomedical Technology (AREA)
  • Emergency Management (AREA)
  • Accounting & Taxation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Agronomy & Crop Science (AREA)
  • Animal Husbandry (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Mining & Mineral Resources (AREA)
  • Primary Health Care (AREA)
  • Computer Security & Cryptography (AREA)
  • Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
  • Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)

Abstract

The invention discloses a method for judging and identifying gas abnormality by using gas desorption quantity characteristics, which comprises the following steps: s1, determining the ratio W of the maximum gas desorption variation in k minutes to the total gas emission amountj,k(ii) a S2, determining the maximum gas desorption change in m minutesRatio W of volume to total gas emission volumej,m(ii) a S3, determining the ratio S of the maximum gas desorption variation in k minutes to the maximum gas desorption variation in m minutesj,k,m(ii) a S4, determining the ratio T of the maximum value in the k-minute moving average sequence to the maximum value in the m-minute moving average sequencej,k,m(ii) a S5, if the total gas emission amount does not tend to 0, carrying out the next step; s6. if Wj,kNot tend to be 0 or Wj,mNot tending to 0 or Sj,k,mIs not prone to
Figure DDA0002509995290000011
Or Tj,k,mIf the trend is not 1, the next step is carried out; s7, if Wj,k、Wj,m、Sj,k,m、Tj,k,mIf at least one of the two is larger, the gas will flow out abnormally. The method for judging and identifying the gas abnormity by using the gas desorption quantity characteristics can accurately judge and identify whether the gas emission is abnormal or not, and has the advantages of good judgment and identification effect, wide application range and low investment cost.

Description

Method for judging and identifying gas abnormality by using gas desorption quantity characteristics
Technical Field
The invention relates to the field of gas emission, in particular to a method for judging and identifying gas abnormity by using gas desorption quantity characteristics.
Background
The coal sample is positioned in the original coal body, when the gas pressure of the coal sample is equal to the gas pressure of the original coal body, the free gas and the adsorbed gas in the coal sample are in a dynamic balance state, once the balance state is damaged, the gas adsorbed in the coal starts to be desorbed, the gas is released, and the phenomenon of gas emission also occurs; the gas emission quantity has important influence on the mining of the mine, and particularly, the judgment on whether the gas emission is abnormal or not has great significance on the safe production of the mine.
At present, a plurality of methods for judging and identifying the gas emission abnormity of the working surface exist, but some of the methods have low accuracy, limited application range, difficult implementation and high input cost.
Therefore, in order to solve the above problems, a method for identifying gas abnormality by using a gas desorption quantity characteristic is needed, whether gas emission is abnormal or not can be accurately identified, the identification effect is good, the application range is wide, the input cost is low, and a foundation is laid for realizing advanced early warning.
Disclosure of Invention
In view of the above, the present invention provides a method for identifying abnormal gas emission by using a characteristic of gas desorption amount, which can accurately identify whether the gas emission is abnormal, and has the advantages of good identification effect, wide application range, and low investment cost.
The method for judging and identifying gas abnormality by using the gas desorption quantity characteristics comprises the following steps:
s1, determining the maximum gas desorption variable quantity Q in k minutes in the jth shiftj,kRatio W to the total gas emission of shift jj,k(ii) a Wherein, the
Figure BDA0002509995270000011
N·Xj,aveIs the total gas emission amount of the jth shift, N is the total minutes of the shift, Xj,aveThe average value of the gas emission concentration per minute of the jth shift is j, and j is the shift number;
s2, determining the maximum gas desorption variable quantity Q in m minutes in the jth shiftj,mRatio W to the total gas emission of shift jj,m(ii) a Wherein, the
Figure BDA0002509995270000021
S3, determining the maximum gas desorption variable quantity Q in k minutes in the jth shiftj,kMaximum gas desorption variation Q within m minutes of the jth shiftj,mRatio S betweenj,k,m(ii) a Wherein, the
Figure BDA0002509995270000022
S4, determining the maximum value X in the k-minute moving average value sequence in the jth shiftj,k,ave,maxAnd the maximum value X in the m-minute moving average sequence in the jth shiftj,m,ave,maxRatio of T betweenj,k,m(ii) a Wherein, the
Figure BDA0002509995270000023
S5, judging whether the total gas emission amount of the jth shift tends to 0, if so, judging that the gas monitoring sensor has a fault or a working surface does not operate, and if not, entering the step S6;
s6, judging Wj,kWhether or not to tend to 0 or Wj,mWhether or not to go to 0 or Sj,k,mWhether or not to tend to
Figure BDA0002509995270000024
Or Tj,k,mWhether the gas tends to 1 or not, if so, the gas monitoring sensor is in failure or the working surface is not operated, otherwise, the operation goes to step S7;
s7, judging Wj,kWhether or not it is greater than a set threshold lambdakOr Wj,mWhether or not it is greater than a set threshold lambdamOr Sj,k,mWhether it is large or notAt set threshold η or Tj,k,mIf the gas emission is larger than the set threshold value mu, the gas emission is abnormal if the gas emission is larger than the set threshold value mu, and the gas emission is normal if the gas emission is not larger than the set threshold value mu.
Further, in step S1, the average gas emission concentration per minute of the jth shift is determined according to the following steps:
s11, collecting the gas emission concentration of the jth shift to obtain a gas emission concentration sequence: x1,X2,…,Xi,…,XN(ii) a Wherein, XiThe gas emission concentration in the time period from the ith minute to the ith minute is shown, i is a minute number, and the value of i is 1,2, … and N; j is the number of the shift; n is the total minutes of the shift;
s12, calculating the average value of the gas emission concentration per minute of the jth shift
Figure BDA0002509995270000025
Further, in step S1, the maximum gas desorption change Q at k minutes in the jth shift is determined according to the following stepsj,k
S101, determining a k-minute moving average value of the gas emission concentration of the jth shift to obtain a k-minute moving average value sequence of the gas emission concentration: x1,k,ave,X2,k,ave,…,Xi,k,ave,…,XN,k,ave(ii) a Wherein, Xi,k,aveThe gas emission concentration X of the jth shiftiK minutes moving average of (a);
s102, determining the maximum value X in the k-minute moving average sequence of the gas emission concentration of the jth shiftj,k,ave,max
S103, calculating the maximum gas desorption variable quantity Q in k minutes in the jth shiftj,kSaid Q isj,k=k·(Xj,k,ave,max-Xj,ave)。
Further, in step S2, the maximum gas desorption change Q in m minutes in the jth shift is determined according to the following stepsj,m
S201, determining the m-minute moving average value of the gas emission concentration of the jth shift to obtain an m-minute moving average value sequence of the gas emission concentration: x1,m,ave,X2,m,ave,…,Xi,m,ave,…,XN,m,ave(ii) a Wherein, Xi,m,aveThe gas emission concentration X of the jth shiftiM minute moving average of (d);
s202, determining the maximum value X in the m-minute moving average sequence of the gas emission concentration of the jth shiftj,m,ave,max
S203, calculating the maximum gas desorption variable quantity Q in m minutes in the jth shiftj,mSaid Q isj,m=m·(Xj,m,ave,max-Xj,ave)。
Further, in step S101, the method
Figure BDA0002509995270000031
Wherein k is a positive integer and takes a value of 5-180.
Further, in step S201, the method
Figure BDA0002509995270000032
Wherein m is a positive integer, m is greater than k, and the value of m is 30-480.
Further, step S7 includes: determining W of multiple adjacent shiftsj,kWhether or not there is an increasing tendency or Wj,mWhether there is an increasing trend or Sj,k,mWhether there is an increasing trend or Tj,k,mAnd if the gas has an increasing trend, the gas is abnormally discharged, and if the gas does not have an increasing trend, the gas is normally discharged.
The invention has the beneficial effects that: the method for judging and identifying the gas abnormity by using the gas desorption quantity characteristics obtains a plurality of gas desorption characteristic quantities by calculating two groups of maximum gas desorption variable quantities, judges whether the gas emission is abnormal according to the sizes of the plurality of gas desorption characteristic quantities, achieves the aim of accurately judging and identifying the gas emission abnormity, has good judgment and identification effects, wide application range and low investment cost, and lays a foundation for realizing advanced early warning.
Drawings
The invention is further described below with reference to the following figures and examples:
FIG. 1 is a schematic diagram illustrating the principle of identifying abnormal gas emission by using gas desorption characteristics according to the present invention;
Detailed Description
The invention is further described with reference to the accompanying drawings, in which:
the method for judging and identifying gas abnormality by using the gas desorption quantity characteristics comprises the following steps:
s1, determining the maximum gas desorption variable quantity Q in k minutes in the jth shiftj,kRatio W to the total gas emission of shift jj,k(ii) a Wherein, the
Figure BDA0002509995270000041
N·Xj,aveIs the total gas emission amount of the jth shift, N is the total minutes of the shift, Xj,aveThe average value of the gas emission concentration per minute of the jth shift is j, and j is the shift number;
s2, determining the maximum gas desorption variable quantity Q in m minutes in the jth shiftj,mRatio W to the total gas emission of shift jj,m(ii) a Wherein, the
Figure BDA0002509995270000042
S3, determining the maximum gas desorption variable quantity Q in k minutes in the jth shiftj,kMaximum gas desorption variation Q within m minutes of the jth shiftj,mRatio S betweenj,k,m(ii) a Wherein, the
Figure BDA0002509995270000043
S4, determining the maximum value X in the k-minute moving average value sequence in the jth shiftj,k,ave,maxAnd the maximum value X in the m-minute moving average sequence in the jth shiftj,m,ave,maxRatio of T betweenj,k,m(ii) a Wherein, the
Figure BDA0002509995270000044
S5, judging whether the total gas emission amount of the jth shift tends to 0, if so, judging that the gas monitoring sensor has a fault or a working surface does not operate, and if not, entering the step S6;
s6, if Wj,kTend to be0、Wj,mTend to 0, Sj,k,mTend to be
Figure BDA0002509995270000051
Tj,k,mIf at least one of the four items of tendency 1 is true, the gas monitoring sensor is in failure or the working surface is not operated, otherwise, the step S7 is executed;
s7, if Wj,kIs greater than a set threshold lambdak、Wj,mIs greater than a set threshold lambdam、Sj,k,mGreater than set threshold η, Tj,k,mIf at least one of the four items is greater than the set threshold value mu, the gas is abnormal to flow out, otherwise, the gas is normal to flow out. In this embodiment, the threshold λk、λmη and mu differ in their values in different modes of operation in the mine, generally the threshold lambdak、λmGreater than 0, threshold η greater than
Figure BDA0002509995270000052
The threshold μ is greater than 1;
in this embodiment, in step S1, the average gas emission concentration per minute of the jth shift is determined according to the following steps:
s11, collecting the gas emission concentration of the jth shift to obtain a gas emission concentration sequence: x1,X2,…,Xi,…,XN(ii) a Wherein, XiThe gas emission concentration in the time period from the ith minute to the ith minute is shown, i is a minute number, and the value of i is 1,2, … and N; n is the total minutes of the shift; in this embodiment, according to the actual working condition of the downhole operation, one shift represents one working period, and one working period is generally 360 to 480 minutes, where the period of the jth shift is set to 480 minutes, that is, N is 480.
S12, calculating the average value of the gas emission concentration per minute of the jth shift
Figure BDA0002509995270000053
In this embodiment, in step S1, the maximum desorption of gas at k minutes in the jth shift is determined according to the following stepsVariation Qj,k
S101, calculating the gas emission concentration of the jth shift to obtain a gas emission concentration sequence: x1,X2,…,Xi,…,XNObtaining k-minute moving average values of N gas emission concentrations by using the k-minute moving average value of each gas emission concentration; sequencing the k-minute moving average values of the N gas emission concentrations according to the sequence of the original gas emission concentration sequence to obtain a k-minute moving average value sequence of the gas emission concentrations: x1,k,ave,X2,k,ave,…,Xi,k,ave,…,XN,k,ave(ii) a Wherein, Xi,k,aveThe gas emission concentration X of the jth shiftiK minutes moving average of (a);
s102, arranging the k-minute moving average sequence of the gas emission concentration according to the sequence of the sizes, and taking the first value after arrangement as the maximum value X in the k-minute moving average sequence of the gas emission concentration of the jth shiftj,k,ave,max
S103, calculating the maximum gas desorption variable quantity Q in k minutes in the jth shiftj,kSaid Q isj,k=k·(Xj,k,ave,max-Xj,ave)。
In this embodiment, in step S2, the maximum gas desorption change Q in m minutes in the jth shift is determined according to the following stepsj,m
S201, calculating the gas emission concentration of the jth shift to obtain a gas emission concentration sequence: x1,X2,…,Xi,…,XNObtaining the m-minute moving average of N gas emission concentrations by the m-minute moving average of each gas emission concentration; sequencing the m-minute moving average values of the N gas emission concentrations according to the sequence of the original gas emission concentration sequence to obtain an m-minute moving average value sequence of the gas emission concentrations: x1,m,ave,X2,m,ave,…,Xi,m,ave,…,XN,m,ave(ii) a Wherein, Xi,m,aveThe gas emission concentration X of the jth shiftiM minute moving average of (d);
s202, arranging the m-minute moving average value sequence of the gas emission concentration according to the sequence of the magnitude, and taking the first value after arrangement as the m-minute moving average of the gas emission concentration of the jth shiftMaximum value X in the mean sequencej,m,ave,max
S203, calculating the maximum gas desorption variable quantity Q in m minutes in the jth shiftj,mSaid Q isj,m=m·(Xj,m,ave,max-Xj,ave)。
In this embodiment, in step S101, the gas emission concentration X of the jth shiftiK minute moving average of
Figure BDA0002509995270000061
Wherein k is a positive integer and takes a value of 5-180.
In this embodiment, in step S201, the gas emission concentration X of the jth shiftiM minute moving average of
Figure BDA0002509995270000062
Wherein m is a positive integer, m is greater than k, and the value of m is 30-480.
In this embodiment, step S7 further includes: if W of multiple adjacent shiftsj,kW tends to increasej,mHas an increasing tendency, Sj,k,mHas an increasing tendency, Tj,k,mIf at least one of the four items with the increasing trend is true, the gas emission quantity is gradually increased, the gas emission is abnormal, otherwise, the gas emission is normal; in this embodiment, the number of adjacent shifts participating in the judgment is adjusted according to the actually obtained gas monitoring data.
Finally, the above embodiments are only for illustrating the technical solutions of the present invention and not for limiting, although the present invention has been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that modifications or equivalent substitutions may be made to the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention, and all of them should be covered in the claims of the present invention.

Claims (7)

1. A method for judging and identifying gas abnormality by using gas desorption quantity characteristics is characterized by comprising the following steps: the method comprises the following steps:
s1, determining the maximum gas desorption change at k minutes in the jth shiftChemical quantity Qj,kRatio W to the total gas emission of shift jj,k(ii) a Wherein, the
Figure FDA0002509995260000011
N·Xj,aveIs the total gas emission amount of the jth shift, N is the total minutes of the shift, Xj,aveThe average value of the gas emission concentration per minute of the jth shift is j, and j is the shift number;
s2, determining the maximum gas desorption variable quantity Q in m minutes in the jth shiftj,mRatio W to the total gas emission of shift jj,m(ii) a Wherein, the
Figure FDA0002509995260000012
S3, determining the maximum gas desorption variable quantity Q in k minutes in the jth shiftj,kMaximum gas desorption variation Q within m minutes of the jth shiftj,mRatio S betweenj,k,m(ii) a Wherein, the
Figure FDA0002509995260000013
S4, determining the maximum value X in the k-minute moving average value sequence in the jth shiftj,k,ave,maxAnd the maximum value X in the m-minute moving average sequence in the jth shiftj,m,ave,maxRatio of T betweenj,k,m(ii) a Wherein, the
Figure FDA0002509995260000014
S5, judging whether the total gas emission amount of the jth shift tends to 0, if so, judging that the gas monitoring sensor has a fault or a working surface does not operate, and if not, entering the step S6;
s6, judging Wj,kWhether or not to tend to 0 or Wj,mWhether or not to go to 0 or Sj,k,mWhether or not to tend to
Figure FDA0002509995260000015
Or Tj,k,mWhether the gas tends to 1 or not, if so, the gas monitoring sensor fails orIf the working surface is not operated, the process goes to step S7;
s7, judging Wj,kWhether or not it is greater than a set threshold lambdakOr Wj,mWhether or not it is greater than a set threshold lambdamOr Sj,k,mWhether greater than a set threshold η or Tj,k,mIf the gas emission is larger than the set threshold value mu, the gas emission is abnormal if the gas emission is larger than the set threshold value mu, and the gas emission is normal if the gas emission is not larger than the set threshold value mu.
2. The method for identifying gas abnormality by using a gas desorption amount characteristic according to claim 1, characterized in that: in step S1, the average gas emission concentration per minute for the jth shift is determined according to the following steps:
s11, collecting the gas emission concentration of the jth shift to obtain a gas emission concentration sequence: x1,X2,…,Xi,…,XN(ii) a Wherein, XiThe gas emission concentration in the time period from the ith minute to the ith minute is shown, i is a minute number, and the value of i is 1,2, … and N; j is the number of the shift; n is the total minutes of the shift;
s12, calculating the average value of the gas emission concentration per minute of the jth shift
Figure FDA0002509995260000021
3. The method for identifying gas abnormality by using a gas desorption amount characteristic according to claim 1, characterized in that: in step S1, the maximum gas desorption change Q at k minutes in the jth shift is determined according to the following stepsj,k
S101, determining a k-minute moving average value of the gas emission concentration of the jth shift to obtain a k-minute moving average value sequence of the gas emission concentration: x1,k,ave,X2,k,ave,…,Xi,k,ave,…,XN,k,ave(ii) a Wherein, Xi,k,aveThe gas emission concentration X of the jth shiftiK minutes moving average of (a);
s102, determining the maximum value X in the k-minute moving average sequence of the gas emission concentration of the jth shiftj,k,ave,max
S103, calculatingMaximum gas desorption variation Q in k minutes within j shiftsj,kSaid Q isj,k=k·(Xj,k,ave,max-Xj,ave)。
4. The method for identifying gas abnormality by using a gas desorption amount characteristic according to claim 1, characterized in that: in step S2, the maximum gas desorption change Q in m minutes in the jth shift is determined according to the following stepsj,m
S201, determining the m-minute moving average value of the gas emission concentration of the jth shift to obtain an m-minute moving average value sequence of the gas emission concentration: x1,m,ave,X2,m,ave,…,Xi,m,ave,…,XN,m,ave(ii) a Wherein, Xi,m,aveThe gas emission concentration X of the jth shiftiM minute moving average of (d);
s202, determining the maximum value X in the m-minute moving average sequence of the gas emission concentration of the jth shiftj,m,ave,max
S203, calculating the maximum gas desorption variable quantity Q in m minutes in the jth shiftj,mSaid Q isj,m=m·(Xj,m,ave,max-Xj,ave)。
5. The method for identifying gas abnormality by using a gas desorption amount characteristic according to claim 1, characterized in that: in step S101, the
Figure FDA0002509995260000031
Wherein k is a positive integer and takes a value of 5-180.
6. The method for identifying gas abnormality by using a gas desorption amount characteristic according to claim 1, characterized in that: in step S201, the
Figure FDA0002509995260000032
Wherein m is a positive integer, m is greater than k, and the value of m is 30-480.
7. Desorption with gas according to claim 1The method for judging and identifying gas abnormality by quantity characteristics is characterized by comprising the following steps: in step S7, the method further includes: determining W of multiple adjacent shiftsj,kWhether or not there is an increasing tendency or Wj,mWhether there is an increasing trend or Sj,k,mWhether there is an increasing trend or Tj,k,mAnd if the gas has an increasing trend, the gas is abnormally discharged, and if the gas does not have an increasing trend, the gas is normally discharged.
CN202010457999.8A 2020-05-26 2020-05-26 Method for judging and identifying gas abnormality by using gas desorption quantity characteristics Active CN111612258B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010457999.8A CN111612258B (en) 2020-05-26 2020-05-26 Method for judging and identifying gas abnormality by using gas desorption quantity characteristics

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010457999.8A CN111612258B (en) 2020-05-26 2020-05-26 Method for judging and identifying gas abnormality by using gas desorption quantity characteristics

Publications (2)

Publication Number Publication Date
CN111612258A true CN111612258A (en) 2020-09-01
CN111612258B CN111612258B (en) 2023-04-07

Family

ID=72196344

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010457999.8A Active CN111612258B (en) 2020-05-26 2020-05-26 Method for judging and identifying gas abnormality by using gas desorption quantity characteristics

Country Status (1)

Country Link
CN (1) CN111612258B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112282733A (en) * 2020-10-29 2021-01-29 中煤科工集团重庆研究院有限公司 Method for determining coal bed gas abnormity by gas emission quantity characteristic while drilling

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1654591A1 (en) * 1986-12-09 1991-06-07 Институт Геотехнической Механики Ан Усср Method for preventing coal and gas outburst
CN101858228A (en) * 2010-04-27 2010-10-13 煤炭科学研究总院重庆研究院 Continuous prediction method of gas emission dynamic characteristic outburst of tunneling surface
CN102926810A (en) * 2012-11-16 2013-02-13 天地(常州)自动化股份有限公司 Forecasting method of coal and gas outburst
US20140019047A1 (en) * 2011-03-30 2014-01-16 China University Of Mining And Technology Multi-information coupling prediction method of coal and gas outburst danger
CN205246645U (en) * 2015-12-18 2016-05-18 河南理工大学 Experimental test system of unusual desorption of gas in coal sample destructive process
CN110219692A (en) * 2019-06-18 2019-09-10 中煤科工集团重庆研究院有限公司 Using coal bed drilling construct in the prominent main contral parameter of Gas data inversion method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1654591A1 (en) * 1986-12-09 1991-06-07 Институт Геотехнической Механики Ан Усср Method for preventing coal and gas outburst
CN101858228A (en) * 2010-04-27 2010-10-13 煤炭科学研究总院重庆研究院 Continuous prediction method of gas emission dynamic characteristic outburst of tunneling surface
US20140019047A1 (en) * 2011-03-30 2014-01-16 China University Of Mining And Technology Multi-information coupling prediction method of coal and gas outburst danger
CN102926810A (en) * 2012-11-16 2013-02-13 天地(常州)自动化股份有限公司 Forecasting method of coal and gas outburst
CN205246645U (en) * 2015-12-18 2016-05-18 河南理工大学 Experimental test system of unusual desorption of gas in coal sample destructive process
CN110219692A (en) * 2019-06-18 2019-09-10 中煤科工集团重庆研究院有限公司 Using coal bed drilling construct in the prominent main contral parameter of Gas data inversion method

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
刘琼等: "瓦斯涌出智能分析及预警技术研究与应用", 《采矿技术》 *
戴林超等: "采空区瓦斯涌出强度对其流动规律的影响研究", 《安全与环境学报》 *
武福生等: "基于瓦斯涌出量的煤体瓦斯解吸特性研究", 《煤炭工程》 *
范磊等: "瓦斯解吸量在突出预测中的可靠性研究", 《山西焦煤科技》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112282733A (en) * 2020-10-29 2021-01-29 中煤科工集团重庆研究院有限公司 Method for determining coal bed gas abnormity by gas emission quantity characteristic while drilling

Also Published As

Publication number Publication date
CN111612258B (en) 2023-04-07

Similar Documents

Publication Publication Date Title
CN109101749B (en) Common cause failure system reliability assessment method considering environmental factors
JP4927271B2 (en) Sensor failure detection, isolation and adjustment methods
CN114024837B (en) Fault root cause positioning method of micro-service system
CN105067262B (en) A kind of state monitoring of rolling bearing method
CN106404441B (en) A kind of failure modes diagnostic method based on non-linear similarity index
CN1867876A (en) Detecting faults of system components in a continuous process
CN112284440B (en) Sensor data deviation self-adaptive correction method
CN111612258B (en) Method for judging and identifying gas abnormality by using gas desorption quantity characteristics
CN111881603B (en) Mechanical structure fatigue reliability assessment method considering failure correlation
CN106570259A (en) Gross error elimination method for dam displacement data
CN103093078B (en) A kind of data checking method improving 53H algorithm
KR102285117B1 (en) Flow control method for absorption type air dryer device using multi dew point meter
CN111307483B (en) Long-period fatigue test data processing and trend prejudging method for mechanical products
CN110362064B (en) Fault diagnosis method for aero-engine exhaust temperature sensor
CN116771611A (en) Wind turbine generator power monitoring method and device
CN113283089B (en) Product reliability assessment method based on double variable threshold
CN115659653A (en) Correlation matrix test point optimization method considering test point reliability
Steiner et al. Effective monitoring of processes with parts per million defective. A hard problem!
WO2019041732A1 (en) Evaluation method and apparatus for manufacturing process capability
CN110729018B (en) Memory diagnosis data compression method based on dynamic fault mode identification
Noble Variations in conventional control charts
CN112015619A (en) Method for optimizing and screening core key indexes of system through parameters
Love et al. Assessment of DESeq2 performance through simulation
CN116543822B (en) Flash memory test and processing method and system applied to solid state disk
Kapur et al. Distribution based change-point problem with two types of imperfect debugging in software reliability

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