CN109993459A - A kind of complexity multi-aquifer water bursting in mine water source recognition methods - Google Patents

A kind of complexity multi-aquifer water bursting in mine water source recognition methods Download PDF

Info

Publication number
CN109993459A
CN109993459A CN201910298332.5A CN201910298332A CN109993459A CN 109993459 A CN109993459 A CN 109993459A CN 201910298332 A CN201910298332 A CN 201910298332A CN 109993459 A CN109993459 A CN 109993459A
Authority
CN
China
Prior art keywords
water
identification
sample
water sample
distinguishing indexes
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
CN201910298332.5A
Other languages
Chinese (zh)
Other versions
CN109993459B (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.)
China University of Mining and Technology CUMT
Anhui University
China Coal Xinji Energy Co Ltd
Original Assignee
Anhui University
China Coal Xinji Energy 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 Anhui University, China Coal Xinji Energy Co Ltd filed Critical Anhui University
Priority to CN201910298332.5A priority Critical patent/CN109993459B/en
Publication of CN109993459A publication Critical patent/CN109993459A/en
Application granted granted Critical
Publication of CN109993459B publication Critical patent/CN109993459B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/213Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
    • G06F18/2132Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on discrimination criteria, e.g. discriminant analysis
    • G06F18/21322Rendering the within-class scatter matrix non-singular
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • G06F18/243Classification techniques relating to the number of classes
    • G06F18/2431Multiple classes
    • 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/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • 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
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/213Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
    • G06F18/2132Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on discrimination criteria, e.g. discriminant analysis
    • G06F18/21322Rendering the within-class scatter matrix non-singular
    • G06F18/21324Rendering the within-class scatter matrix non-singular involving projections, e.g. Fisherface techniques
    • 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
    • Y02A20/00Water conservation; Efficient water supply; Efficient water use
    • Y02A20/152Water filtration

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Human Resources & Organizations (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Economics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Strategic Management (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Engineering & Computer Science (AREA)
  • General Business, Economics & Management (AREA)
  • Evolutionary Computation (AREA)
  • Marketing (AREA)
  • Evolutionary Biology (AREA)
  • Tourism & Hospitality (AREA)
  • Development Economics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Educational Administration (AREA)
  • Artificial Intelligence (AREA)
  • Agronomy & Crop Science (AREA)
  • Primary Health Care (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Mining & Mineral Resources (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Animal Husbandry (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Mathematical Physics (AREA)
  • Sampling And Sample Adjustment (AREA)
  • Investigating Or Analyzing Non-Biological Materials By The Use Of Chemical Means (AREA)

Abstract

The invention discloses a kind of complicated multi-aquifer water bursting in mine water source recognition methods, comprising the following steps: S1: according to the water chemistry data of known aqueous layer water sample, establishes the water source database of mine;S2: by data detection and exceptional value treated data as modeling sample;S3: distinguishing indexes and its threshold value are determined and using the validity of the distinguishing indexes in the preliminary identification model of back substitution test and judge;S4: according to effective distinguishing indexes and Fisher method of identification, " comprehensive-gradually method of identification " model is established;S5: measuring water sample distinguishing indexes to be sentenced, and is successively determined by " comprehensive-gradually method of identification " model, identifies water source type.It is an advantage of the invention that comprehensive using characteristic ion method of comparison, ion ratio Y-factor method Y and Fisher method of identification etc., distinct methods are used to the identification of different water cut layer water source, it is complicated after first simple, gradually determine the Identification of Water Source in Mines type that hydrogeologic condition is more complex, filled water bearing strata is more.

Description

A kind of complexity multi-aquifer water bursting in mine water source recognition methods
Technical field
The present invention relates to water source identification field more particularly to a kind of complicated multi-aquifer water bursting in mine water source recognition methods.
Background technique
China is one using coal as the country of main energy sources, and the coal deposit complicated hydrogeological conditions multiplicity of China, The situation is tense for Pitwater disaster, is an important factor for influencing mine safety production.Mine once occur water damage not only result in it is huge Huge economic loss, and it is likely to result in casualties.Identification of Water Source in Mines identification is the important foundation work of mine water disaster prevention and treatment, Foundation can be provided for water control measures and Post disaster relief etc..And selecting suitable water source recognition methods is then that Identification of Water Source in Mines identifies It is crucial.
Currently, Identification of Water Source in Mines identification main method have groundwater chemical method, Water level trend observation method, isotope method, Coolant-temperature gage analytic approach etc..Wherein, water chemical method is since basic data is abundant, versatility obtains more by force relatively broad application. Method using water chemical method identification water bursting in mine water source has very much, and dependency number is combined mainly in water chemistry data basis Manage the identification of model realization water source.Such as using Fisher identification model, Bayes identification model, apart from identification model, neural network mould Type can open up identification model, clustering, SVM model, and first carry out the built-up pattern etc. that principal component analysis is identified again.
Above-mentioned single method mine hydrogeology condition compared with simple, filled water bearing strata type is less under conditions of, can take Obtain preferable recognition effect.But when hydrogeologic condition is more complex, filled water bearing strata is more, identification validity is substantially reduced, It is unable to satisfy demand of the mine safety production to water source identification model.
Accordingly, it is more to be badly in need of a kind of complexity quickly and effectively identified for being able to achieve complicated multi-aquifer water bursting source type at present Water-bearing layer water bursting in mine water source recognition methods.
Summary of the invention
Technical problem to be solved by the present invention lies in provide a kind of to be able to achieve complicated multi-aquifer water bursting source type The complicated multi-aquifer water bursting in mine water source recognition methods quickly and effectively identified, that it is suitable for hydrogeologic conditions is more complex, fills The more water bursting in mine water source identification in water water-bearing layer.
The present invention solves above-mentioned technical problem using following technical scheme: a kind of knowledge of complexity multi-aquifer water bursting in mine water source Other method, method includes the following steps:
S1: according to the water chemistry data of known aqueous layer water sample, the water source database of mine is established;
S2: by data detection and exceptional value treated data as modeling sample;
S3: distinguishing indexes and its threshold value and having using the distinguishing indexes in the preliminary identification model of back substitution test and judge are determined Effect property;
S4: according to effective distinguishing indexes and Fisher method of identification, " comprehensive-gradually method of identification " model is established;
S5: measuring water sample distinguishing indexes to be sentenced, and is successively determined by " comprehensive-gradually method of identification " model, identifies water source Type.
As one of preferred embodiment of the invention, the water chemistry data of the water sample are each aquifer water of the affiliated coal mine of water sample Various mass of ion concentration and each ion milliequivalent and ratio;Conventional ion can be Na++K+、 Ca2+、Mg2+、Cl-、 HCO3-、SO4 2-Deng;This water chemistry data builds library index as establish water source database.
As one of preferred embodiment of the invention, the method for the data detection is anion-cation balance inspection, error control System is within ± 5%;The screening of exceptional sample and processing method are using each index histogram, case figure, Q~Q figure and clustering The one or more of of figure combine, and the water sample data after filtering out exceptional sample and rejecting are as modeling sample
As one of preferred embodiment of the invention, the selection method of the distinguishing indexes is according to each water-bearing layer water sample of modeling Water chemistry data, the index in each water-bearing layer can be distinguished as distinguishing indexes by therefrom selecting, and determine its threshold value, and identification water Identification relationship when source;The standard of distinguishing indexes selection is changes of contents of the distinguishing indexes in the water sample of water-bearing layer than other indexs Changes of contents it is big;Water sample distinguishing indexes to be measured are less than or greater than determining knowledge by the identification relationship of a certain identification when identifying The water sample of other metrics-thresholds is classified as one kind, and distinguishing indexes can at least distinguish the water sample in two water-bearing layers;Wherein, distinguishing indexes are The mass concentration of conventional ion and the milliequivalent ratio of conventional ion.
As one of preferred embodiment of the invention, the mass concentration of the conventional ion and the milliequivalent ratio of conventional ion Value is called characteristic ion comparison and ion ratio coefficient respectively.
As one of preferred embodiment of the invention, using in the preliminary identification model of back substitution test and judge in the step S3 The validity of distinguishing indexes method particularly includes: using the water sample of known aqueous channel type as new samples, successively substitute into the first of foundation It walks in identification model, if recognition result and reality are consistent, illustrates the preliminary identification model of selection and the foundation of distinguishing indexes It is effective;If same water-bearing layer water sample recognition result is most inconsistent with reality, new distinguishing indexes should be reselected, Establish new preliminary identification model.
As one of preferred embodiment of the invention, it is comprehensive in the step S4-gradually method of identification be it is comprehensive using feature from Sub- method of comparison, ion ratio Y-factor method Y and Fisher method of identification use distinct methods to the identification of different water cut layer water source, first simple It is complicated afterwards, gradually determine water source type;I.e. for a certain wait sentence single water-bearing layer water sample, by water sample characteristic ion to be measured when identification Or ion ratio coefficient is classified as certain water source type less than the water sample of determining distinguishing indexes threshold value, or the index is greater than threshold The water sample of value is classified as certain water source type, as can judge otherwise to be classified as new sample to be sentenced as a result, then stop identifying, New distinguishing indexes and its threshold value are selected this, are so recycled;When to that can not determine water bursting source type with distinguishing indexes, then adopt With Fisher method of identification.It can be realized the quick identification of complicated multi-aquifer mine water outlet water source type using the present invention.
As one of preferred embodiment of the invention, gradually specific step is as follows for method of identification for the synthesis-:
The first step judges that the water sample is ash difficult to understand using characteristic ion method of comparison for a certain wait sentence single water-bearing layer water sample Water or two aqueous, as can judging otherwise to carry out second step identification as a result, then stop identifying;
Second step judges that the water sample is nappe water or coal measures water using ion ratio Y-factor method Y;
Third step, if the water sample is determined as nappe water, uses characteristic ion method of comparison on the basis of second step, Determine that the water sample is nappe gneiss water and the cold buck of nappe;If the water sample is determined as coal measures water, Fisher is used Method of identification determines that the water sample is Sandstone Water or too buck.
So far, by above-mentioned three step, should wait sentence single water-bearing layer water sample be judged as two aqueous, nappe gneiss water, Nappe is trembled with fear buck, Sandstone Water, too one of buck, Ordovician karst water.
As one of preferred embodiment of the invention, the synthesis-gradually method of identification concrete operation method are as follows: passing through data In water sample after verification and abnormity removing: firstly, working as TDS > 4300mg/L, Cl->2000mg/L、 Na++K+When > 1350mg/L, The water sample is Ordovician karst water, end of identification;Otherwise, further identification, as TDS < 440mg/L, Cl-<50mg/L、Na++K+< When 80mg/L, the water sample is two aqueous, end of identification;Otherwise, further identification, as γ Cl-/γCa2+When < 5.2, the water Sample is nappe water, is further identified, as TDS < 1360mg/L, Cl-When < 474mg/L, the water sample is nappe water gneiss Water, end of identification;Otherwise the water sample is the cold buck of nappe water, end of identification;As γ Cl-/γCa2+When > 5.2, the water Sample is coal measures water, recycles Fisher method of identification, identifies Sandstone Water or too buck, at this time end of identification.
The present invention compared with prior art the advantages of be: it is an advantage of the invention that it is comprehensive using characteristic ion method of comparison, from Sub- proportionality coefficient method and Fisher method of identification etc. use distinct methods to the identification of different water cut layer water source, complicated after first simple, i.e., It can gradually determine the Identification of Water Source in Mines type that hydrogeologic condition is more complex, filled water bearing strata is more.The present invention directly passes through detection Distinguishing indexes used by identification process in unknown water sample can substitute into " comprehensive-gradually method of identification " model established, into The identification of row water source.It is thus possible to improve the rapidity and accuracy of complicated multi-aquifer water bursting in mine water source identification.
Detailed description of the invention
Fig. 1 is the complicated multi-aquifer water bursting source identification process figure in embodiment 1;
Fig. 2 is the outlier processing histogram of the present embodiment;
Fig. 3 is the outlier processing cluster result figure of the present embodiment;
Fig. 4 is the outlier processing index case figure of the present embodiment;
Fig. 5 is " comprehensive-gradually method of identification " model specific flow chart of the present embodiment.
Specific embodiment
It elaborates below to the embodiment of the present invention, the present embodiment carries out under the premise of the technical scheme of the present invention Implement, the detailed implementation method and specific operation process are given, but protection scope of the present invention is not limited to following implementation Example.
Embodiment 1
As shown in Figure 1, a kind of complicated multi-aquifer water bursting in mine water source recognition methods of the present embodiment, this method include with Lower step:
S1: according to the water chemistry data of known aqueous layer water sample, the water source database of mine is established;The aquation of the water sample Learn the milliequivalent and ratio of various mass of ion concentration and each ion that data are each aquifer water of the affiliated coal mine of water sample;It is conventional Ion can be Na++K+、Ca2+、Mg2+、Cl-、HCO3 -、SO4 2-Deng;This water chemistry data builds library as establish water source database Index.
S2: by data detection and exceptional value treated data as modeling sample;The side of the data detection Method is anion-cation balance inspection, and control errors are within ± 5%.In the anion-cation balance of the data detection is examined: The screening of exceptional sample and processing method use the one or more of each index histogram, case figure, Q~Q figure and clustering figure In conjunction with the water sample data after filtering out exceptional sample and rejecting are as modeling sample.
S3: distinguishing indexes and its threshold value and having using the distinguishing indexes in the preliminary identification model of back substitution test and judge are determined Effect property;The selection method of the distinguishing indexes is according to the water chemistry data for modeling each water-bearing layer water sample, and therefrom selecting can distinguish The index in each water-bearing layer determines its threshold value as distinguishing indexes, and the identification relationship at identification water source;Distinguishing indexes selection Standard is that changes of contents of the distinguishing indexes in the water sample of water-bearing layer is bigger than the changes of contents of other indexs;The identification of a certain identification is closed System, that is, the water sample that water sample distinguishing indexes to be measured are less than or greater than to determining distinguishing indexes threshold value when identifying are classified as one kind, identify Index can at least distinguish the water sample in two water-bearing layers;Wherein, distinguishing indexes are the mass concentration and conventional ion of conventional ion Milliequivalent ratio, the mass concentration of the conventional ion and the milliequivalent ratio of conventional ion are called characteristic ion respectively Comparison and ion ratio coefficient.Using the effective of the distinguishing indexes in the preliminary identification model of back substitution test and judge in the step S3 Property method particularly includes: using the water sample of known aqueous channel type as new samples, successively substitute into the preliminary identification model of foundation, if Recognition result and reality are consistent, then illustrate that the preliminary identification model of selection and the foundation of distinguishing indexes is effective;If same Water-bearing layer water sample recognition result is most inconsistent with reality, then should reselect new distinguishing indexes, establish new preliminary knowledge Other model.
S4: according to effective distinguishing indexes and Fisher method of identification, " comprehensive-gradually method of identification " model is established;It is described Gradually method of identification is that synthesis uses characteristic ion method of comparison, ion ratio Y-factor method Y and Fisher method of identification to synthesis-in step S4, Distinct methods are used to the identification of different water cut layer water source, it is complicated after first simple, gradually determine water source type;I.e. for a certain wait sentence Water sample characteristic ion to be measured or ion ratio coefficient are less than determining distinguishing indexes threshold value when identification by single water-bearing layer water sample Water sample be classified as certain water source type, or the water sample that the index is greater than threshold value is classified as certain water source type, as can judging to tie Fruit, then stop identifying, is otherwise classified as new sample to be sentenced, and selects this new distinguishing indexes and its threshold value, so recycles; When to that can not determine water bursting source type with distinguishing indexes, then Fisher method of identification is used.Fisher method of identification divides totality Cloth does not have particular requirement, is a kind of linear identification method.Its feature is that high dimensional data point is projected to lower dimensional space is (such as one-dimensional Straight line) on, such data point can become than comparatively dense, so as to overcome due to caused by dimension height " the dimension seed of trouble ". The principle of projection is will to separate as far as possible between totality and totality, then the smallest according between class distance maximum, inter- object distance Principle determines discriminance analysis function, and then new sample is carried out classification knowledge.It can be realized complicated multi-aquifer using the present invention The quick identification of mine water outlet water source type.
Gradually specific step is as follows for method of identification for the synthesis-:
The first step judges that the water sample is ash difficult to understand using characteristic ion method of comparison for a certain wait sentence single water-bearing layer water sample Water or two aqueous, as can judging otherwise to carry out second step identification as a result, then stop identifying;
Second step judges that the water sample is nappe water or coal measures water using ion ratio Y-factor method Y;
Third step, if the water sample is determined as nappe water, uses characteristic ion method of comparison on the basis of second step, Determine that the water sample is nappe gneiss water and the cold buck of nappe;If the water sample is determined as coal measures water, Fisher is used Method of identification determines that the water sample is Sandstone Water or too buck.
So far, by above-mentioned three step, should wait sentence single water-bearing layer water sample be judged as two aqueous, nappe gneiss water, Nappe is trembled with fear buck, Sandstone Water, too one of buck, Ordovician karst water.
S5: measuring water sample distinguishing indexes to be sentenced, and is successively determined by " comprehensive-gradually method of identification " model, identifies water source Type.
The synthesis-gradually method of identification concrete operation method are as follows: in the water sample after data check and abnormity removing: Firstly, working as TDS > 4300mg/L, Cl->2000mg/L、Na++K+When > 1350mg/L, the water sample is Ordovician karst water, end of identification; Otherwise, further identification, as TDS < 440mg/L, Cl-<50mg/L、 Na++K+When < 80mg/L, the water sample is two aqueous, is known Do not terminate;Otherwise, further identification, as γ Cl-/ γCa2+When < 5.2, the water sample is nappe water, is further identified, when TDS<1360mg/L、Cl-When < 474mg/L, the water sample is nappe water gneiss water, end of identification;Otherwise the water sample For the cold buck of nappe water, end of identification;As γ Cl-/γCa2+When > 5.2, the water sample is coal measures water, recycles Fisher Method of identification identifies Sandstone Water or too buck, at this time end of identification.
The present embodiment is comprehensive using characteristic ion method of comparison, ion ratio Y-factor method Y and Fisher method of identification etc., to difference Aquifer water identifing source uses distinct methods, complicated after first simple, can gradually determine that hydrogeologic condition is more complex, water-filling contains The more Identification of Water Source in Mines type of water layer.The present embodiment directly passes through identification used by identification process in the unknown water sample of detection and refers to Mark can substitute into " comprehensive-gradually method of identification " model established, and carry out water source identification.It is thus possible to improve complicated contain more The rapidity and accuracy of water layer water bursting in mine water source identification.
Selection below is new to collect aqueous Cenozoic two in two water quality accounts, nappe gneiss water, nappe cold buck, sandstone The water-bearing layers such as water and too buck water sample data are as database, to be specifically described the design scheme and theoretical foundation of the present embodiment.
Data detection and exceptional sample processing
When carrying out single characteristics on Aquifer index and judging, it is desirable that the water sample Hydrochemical Composition analyzed need to be the water-bearing layer Typical Representative.Therefore, it before carrying out characteristic index analysis and being established with water source identification model, need to test to each water-bearing layer water sample With exceptional sample processing.
The method that sample data is examined is anion-cation balance inspection, due to Na++K+、Ca2+、Mg2+、Cl-、 HCO3 -、 SO4 2-Index is measured value, therefore it is required that control errors within ± 5%, i.e., test error is qualified water within this range Sample carries out next step analysis as basic data.
Though it is qualified to examine through zwitterion, there are a small number of water samples durings water sampling, the judgement of layer position, test analysis etc. A certain link there are some problems, cause the water sample that can not really reflect a certain water-bearing layer water quality characteristic, in water analysis Outliers are shown as in the process, need to carry out screening rejecting processing.Below by taking Ordovician karst water sample as an example, the sieve to exceptional sample is illustrated Choosing and processing method, such as Fig. 2-4.
To 17 Ordovician karst water samples TDS, Na of collection++K+、Ca2+、Cl-Ion compares and analyzes, as shown in Fig. 2.Label It is obvious for the above-mentioned several indexs of tri- water samples of O-03, O-04, O-05 and remaining 14 water sample difference.Further to judge Ordovician karst water O-04, O-05 and No. O-06 three water samples and remaining water sample relationship, have carried out hierarchial-cluster analysis to Ordovician karst water sample.Using TDS, Na++K+、Ca2+、Mg2+、Cl-、HCO3 -、 SO4 2-Deng 7 indexs as variable, cluster result is as shown in Figure 3.As seen from Figure 3, it removes Outside marked as O-03, O-04, O-05 sample, remaining sample room otherness is minimum, can most represent Ordovician karst water property, and marked as O-03, O-04, O-05 and remaining sample interval finally just complete classification from maximum.And box traction substation, it is that one kind is used to describe number According to the statistical graph of distribution, it can be used to show the descriptive statistics amounts such as median, 4 quantiles and the extreme value of observation data, from The distribution situation of the angular observation variate-value of vision.Ordovician karst water 17 water samples TDS, Na++K+、Ca2+、Mg2+、Cl-Index case figure is such as Shown in Fig. 4.As seen from the figure, marked as tri- water samples of O-3, O-4, O-5 in this 5 indexs, also display has with remaining water sample Significant difference.
Therefore, O-03, O-04, No. O-05 three samples can be considered exceptional sample in 17 samples of Ordovician karst water, cannot represent Ordovician karst water characteristic feature, will reject in subsequent water sample type identification model foundation, participate in modeling not as standard water sample.
, nappe gneiss water aqueous to Cenozoic two, nappe Cambrian system water, Sandstone Water and too buck etc. remaining 5 Single water-bearing layer water sample is all made of the synthesis screening and rejecting that above-mentioned or other a variety of methods carry out outliers.It is examined by data Test and exceptional sample reject after data, the sample as modeling.
Distinguishing indexes and threshold value determine
According to Such analysis, two aqueous and Ordovician karst waters are in TDS, Cl-And Na++K+There is characteristic feature in content, and it is other Water-bearing layer water sample can be distinguished effectively.Therefore, select above-mentioned 3 indexs as two is aqueous and the distinguishing indexes of Ordovician karst water, Its content range is as shown in table 1.
Table 1 two is aqueous and Ordovician karst water characteristic index content range
(2) nappe water and coal measures water
Nappe gneiss water is close with the cold buck conventional ion of nappe and water quality type, and Sandstone Water is conventional with too buck Ion and water quality type are close.Here nappe gneiss water and the cold buck of nappe are first classified as one kind, are referred to as nappe Water;Sandstone Water and too buck is classified as one kind, is referred to as coal measures water.By modeling sample statistics indicate that, nappe water and coal It is characteristic ion proportionality coefficient (the γ Cl of water-/γCa2+) otherness is obvious (table 2), can be used as the nappe water reclassified and The diacritics of coal measures water.
2 nappe water of table and coal measures water characteristic ion proportionality coefficient scope
(3) nappe gneiss water and the cold buck of nappe
Nappe gneiss water and nappe Cambrian system water are in TDS and Cl-There is characteristic feature in content, can be used as differentiation The mark of these two types of aquifer waters.
3 nappe water characteristic index content range of table
(4) Sandstone Water and too buck
The difference in feature example index and ion ratio coefficient is unobvious with too buck for Sandstone Water, only from above-mentioned two Aspect conventional means are not possible to distinguish.For this purpose, comprehensive selection TDS, Ca2+、Mg2+、Na++K+、 HCO3 -、SO4 2-、Cl-Refer to Deng 7 It is denoted as the recognition factor for Sandstone Water and too classified between buck, judgement water source type is carried out using Fisher method of identification.
Fisher method of identification is a kind of linear identification method to the overall no particular requirement of distribution.Its feature be by High dimensional data point projects on lower dimensional space (such as one-dimensional straight line), and such data point can become than comparatively dense, so as to Overcome due to caused by dimension height " the dimension seed of trouble ".The principle of projection is will to separate as far as possible between totality and totality, then According between class distance, maximum, the smallest principle of inter- object distance determines discriminance analysis function, and then new sample is carried out classification knowledge.
Recognition methods selection and identification step
This two mine complexity multi-aquifer water source type recognition methods of new collection is using " comprehensive-gradually method of identification ".
Specific step is as follows:
The first step judges that the water sample is ash difficult to understand using characteristic ion method of comparison for a certain wait sentence single water-bearing layer water sample Water or two aqueous, as can judging otherwise to carry out second step identification as a result, then stop identifying.
Second step judges that the water sample is nappe water or coal measures water using ion ratio Y-factor method Y.
Third step, if the water sample is determined as nappe water, uses characteristic ion method of comparison on the basis of second step, Determine that the water sample is nappe gneiss water and the cold buck of nappe;If the water sample is determined as coal measures water, Fisher is used Method of identification determines that the water sample is Sandstone Water and too buck.
So far, by above-mentioned three step, should wait sentence single water-bearing layer water sample be judged as two aqueous, nappe gneiss water, Nappe is trembled with fear buck, Sandstone Water, too one of buck, Ordovician karst water.
The foundation of " comprehensive-gradually method of identification " model
For a certain wait sentence single water-bearing layer water sample, data detection and exceptional sample processing, distinguishing indexes and boundary are being carried out After limit value, recognition methods selection and identification step determine, the identification model and process of foundation are as shown in Figure 3.
The inspection of recognition effect
Using the validity of back substitution test and judge distinguishing indexes, using the water sample of known aqueous channel type as new samples, according to It is secondary to substitute into the identification model established, if recognition result and practical consistent, illustrate distinguishing indexes selection and foundation it is " comprehensive Closing-gradually method of identification " model is effective;If same water-bearing layer water sample recognition result is most inconsistent with reality, should weigh New distinguishing indexes are newly selected, new water source identification model is established.
After determining " comprehensive-gradually method of identification " model, to unknown water-bearing layer water sample, pass through the content of test distinguishing indexes Afterwards, established model is substituted into, the aqueous channel type of you can get it unknown water sample.
The foregoing is merely illustrative of the preferred embodiments of the present invention, is not intended to limit the invention, all in essence of the invention Made any modifications, equivalent replacements, and improvements etc., should all be included in the protection scope of the present invention within mind and principle.

Claims (9)

1. a kind of complexity multi-aquifer water bursting in mine water source recognition methods, which is characterized in that method includes the following steps:
S1: according to the water chemistry data of known aqueous layer water sample, the water source database of mine is established;
S2: by data detection and exceptional value treated data as modeling sample;
S3: distinguishing indexes and its threshold value are determined and using the effective of the distinguishing indexes in the preliminary identification model of back substitution test and judge Property;
S4: according to effective distinguishing indexes and Fisher method of identification, " comprehensive-gradually method of identification " model is established;
S5: measuring water sample distinguishing indexes to be sentenced, and is successively determined by " comprehensive-gradually method of identification " model, identifies water source type.
2. complexity multi-aquifer water bursting in mine water source recognition methods according to claim 1, which is characterized in that the water sample Water chemistry data be each aquifer water of the affiliated coal mine of water sample each mass of ion concentration and each ion milliequivalent and ratio; This water chemistry data builds library index as establish water source database.
3. complexity multi-aquifer water bursting in mine water source recognition methods according to claim 1, which is characterized in that the data The method of inspection is anion-cation balance inspection, and control errors are within ± 5%;The screening of exceptional sample and processing method use Each index histogram, case figure, Q~Q figure and the one or more of of clustering figure combine, after filtering out exceptional sample and rejecting Water sample data are as modeling sample.
4. complexity multi-aquifer water bursting in mine water source recognition methods according to claim 1, which is characterized in that the identification The selection method of index is therefrom to select the index that can distinguish each water-bearing layer according to the water chemistry data for modeling each water-bearing layer water sample As distinguishing indexes, and determine its threshold value, and identification relationship when identification water source;The standard of distinguishing indexes selection is that identification refers to The changes of contents being marked in the water sample of water-bearing layer is bigger than the changes of contents of other indexs;The identification relationship of a certain identification, that is, when identifying The water sample that water sample distinguishing indexes to be measured are less than or greater than to determining distinguishing indexes threshold value is classified as one kind;Distinguishing indexes at least can be with Distinguish the water sample in two water-bearing layers;Wherein, distinguishing indexes are the mass concentration of conventional ion and the milliequivalent ratio of conventional ion Value.
5. complexity multi-aquifer water bursting in mine water source recognition methods according to claim 4, which is characterized in that the routine The mass concentration of ion and the milliequivalent ratio of conventional ion are called characteristic ion comparison and ion ratio coefficient respectively.
6. complexity multi-aquifer water bursting in mine water source recognition methods according to claim 1, which is characterized in that the step Using the validity of the distinguishing indexes in the preliminary identification model of back substitution test and judge in S3 method particularly includes: by known aqueous layer class The water sample of type successively substitutes into the preliminary identification model of foundation as new samples, if recognition result and reality are consistent, illustrates to know The preliminary identification model of selection and the foundation of other index is effective;If same water-bearing layer water sample recognition result majority and reality It is inconsistent, then new distinguishing indexes should be reselected, new preliminary identification model is established.
7. complexity multi-aquifer water bursting in mine water source recognition methods according to claim 1, which is characterized in that the step Gradually method of identification is comprehensive using characteristic ion method of comparison, ion ratio Y-factor method Y and Fisher method of identification to synthesis-in S4, to not Distinct methods are used with aquifer water identifing source, it is complicated after first simple, gradually determine water source type;It is i.e. single wait sentence for a certain Water sample characteristic ion to be measured or ion ratio coefficient are less than the water of determining distinguishing indexes threshold value when identification by water-bearing layer water sample Sample is classified as certain water source type, or the water sample that the index is greater than threshold value is classified as certain water source type, as can judging as a result, then Stop identification, be otherwise classified as new sample to be sentenced, select this new distinguishing indexes and its threshold value, so recycles;To nothing When method determines water bursting source type with distinguishing indexes, then Fisher method of identification is used.
8. complexity multi-aquifer water bursting in mine water source recognition methods according to claim 7, which is characterized in that described comprehensive Close-gradually specific step is as follows for method of identification:
The first step, for a certain wait sentence single water-bearing layer water sample, using characteristic ion method of comparison judge the water sample be Ordovician karst water also It is two aqueous, as can judging as a result, then stopping identifying, otherwise progress second step identification;
Second step judges that the water sample is nappe water or coal measures water using ion ratio Y-factor method Y;
Third step if the water sample is determined as nappe water, uses characteristic ion method of comparison, determines on the basis of second step The water sample is nappe gneiss water and the cold buck of nappe;If the water sample is determined as coal measures water, identified using Fisher Method determines that the water sample is Sandstone Water or too buck.
9. complexity multi-aquifer water bursting in mine water source recognition methods according to claim 8, which is characterized in that described comprehensive Close-gradually method of identification concrete operation method are as follows: in the water sample after data check and abnormity removing: firstly, work as TDS > 4300mg/L、Cl->2000mg/L、Na++K+When > 1350mg/L, the water sample is Ordovician karst water, end of identification;Otherwise, further Identification, as TDS < 440mg/L, Cl-<50mg/L、Na++K+When < 80mg/L, the water sample is two aqueous, end of identification;Otherwise, Further identification, as γ Cl-/γCa2+When < 5.2, the water sample be nappe water, further identify, when TDS < 1360mg/L, Cl-When < 474mg/L, the water sample is nappe water gneiss water, end of identification;Otherwise the water sample is the cold ash of nappe water Water, end of identification;As γ Cl-/γCa2+When > 5.2, the water sample is coal measures water, recycles Fisher method of identification, and identification is shaked out Rock water or too buck, at this time end of identification.
CN201910298332.5A 2019-04-15 2019-04-15 Method for identifying water bursting source of complex multi-aquifer mine Active CN109993459B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910298332.5A CN109993459B (en) 2019-04-15 2019-04-15 Method for identifying water bursting source of complex multi-aquifer mine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910298332.5A CN109993459B (en) 2019-04-15 2019-04-15 Method for identifying water bursting source of complex multi-aquifer mine

Publications (2)

Publication Number Publication Date
CN109993459A true CN109993459A (en) 2019-07-09
CN109993459B CN109993459B (en) 2022-09-23

Family

ID=67133385

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910298332.5A Active CN109993459B (en) 2019-04-15 2019-04-15 Method for identifying water bursting source of complex multi-aquifer mine

Country Status (1)

Country Link
CN (1) CN109993459B (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110852364A (en) * 2019-10-31 2020-02-28 中国煤炭地质总局勘查研究总院 Method and device for identifying water source of water burst in mine and electronic equipment
CN111562285A (en) * 2020-06-03 2020-08-21 安徽大学 Mine water inrush source identification method and system based on big data and deep learning
CN112381117A (en) * 2020-10-22 2021-02-19 合肥工业大学 Coal mine water inrush water source mixing proportion calculation and dynamic monitoring method based on conventional water chemistry
CN112508330A (en) * 2020-10-29 2021-03-16 中煤科工集团西安研究院有限公司 Method for distinguishing mine water source under disturbance of western mining area mining
CN112633318A (en) * 2020-11-04 2021-04-09 中国地质大学(北京) Water source identification method based on Java and android platform
CN113689915A (en) * 2021-08-25 2021-11-23 合肥工业大学 Abnormal water guide structure detection method based on conventional water chemistry mine
CN114167021A (en) * 2021-12-09 2022-03-11 山西启诚电子科技有限公司 Mine water source rapid identification instrument and control method thereof
CN114264680A (en) * 2021-11-15 2022-04-01 中煤科工集团西安研究院有限公司 Method for predicting fluoride concentration in mine water based on analogy method
CN118033091A (en) * 2024-04-09 2024-05-14 华北理工大学 Water source discrimination method for silica-karite type iron ore deposit based on ore formation analysis

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101923084A (en) * 2010-07-15 2010-12-22 北京华安奥特科技有限公司 A kind of mining water source recognition methods and identification equipment
CN103617147A (en) * 2013-11-27 2014-03-05 中国地质大学(武汉) Method for identifying mine water-inrush source
CN104122319A (en) * 2014-08-13 2014-10-29 北京华安奥特科技有限公司 Method and system for identifying water source in mining area based on ion composite electrode detecting technology and spectrum analysis technology
CN104297308A (en) * 2014-10-23 2015-01-21 淮南矿业(集团)有限责任公司 Device for rapidly distinguishing water inrush source of mine
CN104597516A (en) * 2015-01-19 2015-05-06 天地科技股份有限公司 Quick distinguishing system for water bursting source of mine
WO2016115816A1 (en) * 2015-01-23 2016-07-28 中国矿业大学 Structural discrimination indexes of ordovician limestone top filling zones and determination method
WO2017004552A1 (en) * 2015-07-02 2017-01-05 Abramowitz Rachel Customized personal care management based on water quality and other environmental factors
CN108805357A (en) * 2018-06-13 2018-11-13 安徽理工大学 A kind of Fisher discrimination model water bursting source prediction techniques based on PCA analyses
CN108876030A (en) * 2018-06-11 2018-11-23 安徽理工大学 A kind of water bursting source prediction technique based on Fisher discrimination model
WO2019019709A1 (en) * 2017-07-24 2019-01-31 厦门快商通科技股份有限公司 Method for detecting water leakage of tap water pipe

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101923084A (en) * 2010-07-15 2010-12-22 北京华安奥特科技有限公司 A kind of mining water source recognition methods and identification equipment
CN103617147A (en) * 2013-11-27 2014-03-05 中国地质大学(武汉) Method for identifying mine water-inrush source
CN104122319A (en) * 2014-08-13 2014-10-29 北京华安奥特科技有限公司 Method and system for identifying water source in mining area based on ion composite electrode detecting technology and spectrum analysis technology
CN104297308A (en) * 2014-10-23 2015-01-21 淮南矿业(集团)有限责任公司 Device for rapidly distinguishing water inrush source of mine
CN104597516A (en) * 2015-01-19 2015-05-06 天地科技股份有限公司 Quick distinguishing system for water bursting source of mine
WO2016115816A1 (en) * 2015-01-23 2016-07-28 中国矿业大学 Structural discrimination indexes of ordovician limestone top filling zones and determination method
WO2017004552A1 (en) * 2015-07-02 2017-01-05 Abramowitz Rachel Customized personal care management based on water quality and other environmental factors
WO2019019709A1 (en) * 2017-07-24 2019-01-31 厦门快商通科技股份有限公司 Method for detecting water leakage of tap water pipe
CN108876030A (en) * 2018-06-11 2018-11-23 安徽理工大学 A kind of water bursting source prediction technique based on Fisher discrimination model
CN108805357A (en) * 2018-06-13 2018-11-13 安徽理工大学 A kind of Fisher discrimination model water bursting source prediction techniques based on PCA analyses

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
HUANG PING-HUA等: "Recognition model of groundwater inrush source of coal mine: a case study on Jiaozuo coal mine in China", 《ARABIAN JOURNAL OF GEOSCIENCES》 *
刘杰刚等: "多元统计分析模型在矿井突水水源判别中的应用", 《中国煤炭》 *
姚洁等: "基于聚类分析方法的矿井水源识别", 《煤矿安全》 *
金洲洋: "新集一矿突水水源综合判别模型研究", 《中国优秀硕士学位论文全文数据库(工程科技I辑)》 *
陈红江等: "用Fisher判别法确定矿井突水水源", 《中南大学学报(自然科学版)》 *
黄平华等: "基于多元统计分析的矿井突水水源Fisher识别及混合模型", 《煤炭学报》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110852364A (en) * 2019-10-31 2020-02-28 中国煤炭地质总局勘查研究总院 Method and device for identifying water source of water burst in mine and electronic equipment
CN111562285A (en) * 2020-06-03 2020-08-21 安徽大学 Mine water inrush source identification method and system based on big data and deep learning
CN112381117A (en) * 2020-10-22 2021-02-19 合肥工业大学 Coal mine water inrush water source mixing proportion calculation and dynamic monitoring method based on conventional water chemistry
CN112381117B (en) * 2020-10-22 2023-10-17 合肥工业大学 Coal mine water inrush source mixing proportion calculation and dynamic monitoring method based on conventional water chemistry
CN112508330A (en) * 2020-10-29 2021-03-16 中煤科工集团西安研究院有限公司 Method for distinguishing mine water source under disturbance of western mining area mining
CN112508330B (en) * 2020-10-29 2024-04-23 中煤科工集团西安研究院有限公司 Method for distinguishing mine water sources under disturbance of mining in western mining area
CN112633318B (en) * 2020-11-04 2023-08-11 中国地质大学(北京) Water source identification method based on Java and android platforms
CN112633318A (en) * 2020-11-04 2021-04-09 中国地质大学(北京) Water source identification method based on Java and android platform
CN113689915A (en) * 2021-08-25 2021-11-23 合肥工业大学 Abnormal water guide structure detection method based on conventional water chemistry mine
CN114264680A (en) * 2021-11-15 2022-04-01 中煤科工集团西安研究院有限公司 Method for predicting fluoride concentration in mine water based on analogy method
CN114264680B (en) * 2021-11-15 2023-06-13 中煤科工集团西安研究院有限公司 Method for predicting concentration of fluoride in mine water based on comparison method
CN114167021A (en) * 2021-12-09 2022-03-11 山西启诚电子科技有限公司 Mine water source rapid identification instrument and control method thereof
CN118033091A (en) * 2024-04-09 2024-05-14 华北理工大学 Water source discrimination method for silica-karite type iron ore deposit based on ore formation analysis

Also Published As

Publication number Publication date
CN109993459B (en) 2022-09-23

Similar Documents

Publication Publication Date Title
CN109993459A (en) A kind of complexity multi-aquifer water bursting in mine water source recognition methods
CN103206208B (en) A kind of method of microcosmic different existence state remaining oil macro-level quantitative
CN103630528B (en) A kind of method utilizing constituent content discriminating tea-leaf producing area in tealeaves
CN110737874B (en) Watershed water quality monitoring abnormal value detection method based on spatial relationship
CN108229811B (en) Method for evaluating water injection effect of fractured-vuggy carbonate reservoir
CN107132266A (en) A kind of Classification of water Qualities method and system based on random forest
CN110261560A (en) The water source recognition methods of complex hydrologic geology water bursting in mine and system
CN110458201A (en) A kind of remote sensing image object-oriented classification method and sorter
CN113516228B (en) Network anomaly detection method based on deep neural network
CN110276264A (en) A kind of crowd density estimation method based on foreground segmentation figure
CN103558311B (en) A kind of bitter taste of green tea method of discrimination based on Tea ingredient
CN109801685A (en) A kind of water bursting source recognition methods based on PCA method Yu Bayes discrimination model
CN109779622B (en) Method and device for characterizing low-efficiency water injection zone of oil reservoir in ultrahigh water cut period
CN114813903A (en) Method for distinguishing ore species based on garnet micro-area chemical components
CN103424420A (en) Fitting-based nuclear magnetic resonance signal processing method
CN110334773A (en) Model based on machine learning enters the screening technique of modular character
CN111832389A (en) Counting and analyzing method of bone marrow cell morphology automatic detection system
CN105089642B (en) The recognition methods of oil-gas reservoir exception producing well and device
CN114332856A (en) Feature extraction method and device for blood type card image and blood type interpretation system
CN106019402A (en) Method for comprehensively discriminating and verifying paleoclimatic cycle inside extremely thick coal seam
CN106054859B (en) The double-deck integrated form industrial process fault detection method based on amendment type independent component analysis
Bakraji Application of multivariate statistical methods to classify archaeological pottery from Tel‐Alramad site, Syria, based on x‐ray fluorescence analysis
CN105824785A (en) Rapid abnormal point detection method based on penalized regression
CN113205161B (en) Traditional Chinese medicine producing area distinguishing system and method based on soil parameters
CN106777724B (en) Basic flow segmentation method for semi-humid and semi-arid regions and humid regions

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
TR01 Transfer of patent right

Effective date of registration: 20240607

Address after: 221008 No.1, Jinshan East Road, Quanshan District, Xuzhou City, Jiangsu Province

Patentee after: CHINA University OF MINING AND TECHNOLOGY

Country or region after: China

Patentee after: ANHUI University

Patentee after: CHINA COAL XINJI ENERGY Co.,Ltd.

Address before: No. 111, Kowloon Road, Hefei, Anhui, Anhui

Patentee before: ANHUI University

Country or region before: China

Patentee before: CHINA COAL XINJI ENERGY Co.,Ltd.