Disclosure of Invention
The invention aims to solve the technical problem of providing a complex multi-aquifer mine water inrush source identification method capable of realizing quick and effective identification of complex multi-aquifer water inrush source types, which is suitable for identifying mine water inrush sources with complex hydrogeological conditions and more water-filled aquifers.
The invention adopts the following technical scheme to solve the technical problems: a method for identifying a water inrush source of a complex multi-aquifer mine comprises the following steps:
s1: establishing a water source database of a mine according to the water chemistry data of the known aquifer water sample;
s2: taking the data subjected to data detection and abnormal value processing as a modeling sample;
s3: determining an identification index and a threshold value thereof, and judging the effectiveness of the identification index in the primary identification model by adopting back-up test;
s4: establishing a comprehensive-step recognition method model according to the effective recognition index and a Fisher recognition method;
s5: and determining the identification indexes of the water sample to be judged, and sequentially judging through a comprehensive-step identification method model to identify the type of the water source.
As one of the preferable modes of the invention, the water chemistry data of the water sample is the mass concentration of various ions and the milliequivalent and ratio of each ion of the water in each aquifer of the coal mine to which the water sample belongs; the conventional ion may be Na + +K + 、 Ca 2+ 、Mg 2+ 、Cl - 、HCO3 - 、SO 4 2- Etc.; the water chemistry data is used as a database building index for building a water source database.
As one of the preferable modes of the invention, the data inspection method is anion-cation balance inspection, and the error is controlled within +/-5%; the abnormal sample screening and processing method adopts one or more combinations of index histogram, box diagram, Q-Q diagram and cluster analysis diagram to screen the abnormal sample and remove the water sample data as the modeling sample
As one of the preferable modes of the invention, the selection method of the identification index is to select the index capable of distinguishing each aquifer as the identification index according to the water chemistry data of the water samples of each aquifer, determine the threshold value of the index and identify the identification relation when identifying the water source; the standard of the identification index selection is that the content change of the identification index in the aquifer water sample is larger than the content change of other indexes; the identification relation of a certain identification means that water samples of which the identification indexes are smaller than or larger than a determined identification index threshold value are classified into a class during identification, and the identification indexes can at least distinguish the water samples of two aquifers; wherein the identification index is the mass concentration of the conventional ions and the milliequivalent ratio of the conventional ions.
As one of the preferred modes of the present invention, the mass concentration of the conventional ions and the milliequivalent ratio of the conventional ions are called the characteristic ion contrast and the ion proportionality coefficient, respectively.
As one preferable embodiment of the present invention, the specific method for judging the validity of the identification index in the preliminary identification model by using the back substitution test in step S3 includes: taking a water sample containing a known water layer type as a new sample, sequentially substituting the new sample into the established preliminary identification model, and if the identification result is consistent with the actual result, indicating that the selection of the identification index and the established preliminary identification model are effective; and if the identification results of the same aquifer water sample are mostly inconsistent with the actual identification results, reselecting a new identification index and establishing a new initial identification model.
As one of the preferable modes of the invention, the comprehensive-gradual identification method in step S4 is to comprehensively adopt a characteristic ion comparison method, an ion proportion coefficient method and a Fisher identification method, and different methods are adopted for identifying different aquifer water sources, wherein the water source types are determined step by step after being simple and complex; for a single aquifer water sample to be judged, classifying characteristic ions or water samples of which the ion proportionality coefficients are smaller than a determined identification index threshold value into a certain water source type during identification, or classifying the water samples of which the indexes are larger than the threshold value into the certain water source type, stopping identification if a result can be judged, or classifying the water samples into a new sample to be judged, selecting a new identification index and a threshold value thereof, and repeating the steps; and when the type of the water inrush source cannot be judged by using the identification index, a Fisher identification method is adopted. The method can realize the rapid identification of the type of the complex multi-aquifer mine effluent source.
As one of the preferable modes of the invention, the comprehensive-step-by-step identification method comprises the following specific steps:
step one, for a certain single aquifer water sample to be judged, judging whether the water sample is AoHu water or Duhui water by adopting a characteristic ion comparison method, if the result can be judged, stopping identification, otherwise, performing identification in the step two;
secondly, judging whether the water sample is cladding body water or coal series water by adopting an ion proportionality coefficient method;
thirdly, on the basis of the second step, if the water sample is judged to be the cladding body water, judging that the water sample is cladding body gneiss water and cladding body cold grey water by adopting a characteristic ion comparison method; and if the water sample is judged to be coal-based water, judging that the water sample is sandstone water or Taigrey water by adopting a Fisher recognition method.
And judging the single aquifer water sample to be judged into one of secondary water, overlying body gneiss water, overlying body cold grey water, sandstone water, taigrey water and Ordovician grey water through the three steps.
As one of the preferable modes of the invention, the specific operation method of the comprehensive-step-by-step identification method is as follows: in the water sample after data verification and abnormal rejection: first, when TDS>4300mg/L、Cl - >2000mg/L、 Na + +K + >When 1350mg/L, the water sample is the Ordovician limestone water, and the identification is finished; otherwise, further identifying when TDS<440mg/L、Cl - <50mg/L、Na + +K + <When the water sample is 80mg/L, the water sample is secondary water, and the identification is finished; otherwise, further identify when Gamma Cl - /γCa 2+ <5.2, the water sample is the cover-pushing body water, and is further identified when TDS<1360mg/L、Cl - <474mg/L, the water sample is the water of the coating body water sheet gneiss waterRespectively finishing; otherwise, the water sample is the push-coating body water cold grey water, and the identification is finished; when gamma Cl is present - /γCa 2+ >And 5.2, identifying the sandstone water or the taigrey water by using a Fisher identification method, wherein the water sample is coal-based water, and the identification is finished at the moment.
Compared with the prior art, the invention has the advantages that: the method has the advantages that a characteristic ion contrast method, an ion proportion coefficient method, a Fisher identification method and the like are comprehensively adopted, different methods are adopted for identifying different aquifer water sources, and the mine water source type with complicated hydrogeological conditions and more water-filled aquifers can be gradually judged after the method is simple and complicated. The method can substitute the identification indexes adopted by the identification process in the unknown water sample into the established comprehensive-gradual identification method model to identify the water source directly. Therefore, the rapidity and the accuracy of identifying the complex multi-aquifer mine water inrush source can be improved.
Detailed Description
The following examples are given for the detailed implementation and specific operation of the present invention, but the scope of the present invention is not limited to the following examples.
Example 1
As shown in fig. 1, the method for identifying a water inrush source of a complex multi-aquifer mine in the embodiment includes the following steps:
s1: establishing a water source database of a mine according to water chemistry data of a known aquifer water sample; the water chemical data of the water sample are all the coal mines to which the water sample belongsMass concentrations of various ions of the water of the aquifer and milliequivalents and ratios of the ions; the conventional ion may be Na + +K + 、Ca 2+ 、Mg 2+ 、Cl - 、HCO 3 - 、SO 4 2- Etc.; the water chemistry data is used as a database building index for building a water source database.
S2: data after data inspection and abnormal value processing are used as modeling samples; the data inspection method is anion-cation balance inspection, and the error is controlled within +/-5%. In the anion-cation balance test of the data test: the abnormal sample screening and processing method adopts one or more combinations of index histograms, box charts, Q-Q charts and cluster analysis charts to screen the abnormal samples and remove the water sample data as a modeling sample.
S3: determining an identification index and a threshold value thereof, and judging the effectiveness of the identification index in the primary identification model by adopting back substitution test; the selection method of the identification indexes comprises the steps of selecting indexes capable of distinguishing each aquifer from water chemical data of water samples of each aquifer according to modeling as identification indexes, determining threshold values of the indexes and identifying the identification relation of water sources; the standard of the identification index selection is that the content change of the identification index in the aquifer water sample is larger than that of other indexes; the identification relation of a certain identification means that water samples of which the identification indexes are smaller than or larger than a determined identification index threshold value are classified into a class during identification, and the identification indexes can at least distinguish the water samples of two aquifers; the identification indexes are the mass concentration of the conventional ions and the milliequivalent ratio of the conventional ions, and the mass concentration of the conventional ions and the milliequivalent ratio of the conventional ions are respectively called characteristic ion contrast and an ion proportionality coefficient. The specific method for judging the validity of the identification index in the preliminary identification model by adopting the back-substitution test in the step S3 comprises the following steps: taking a water sample containing a known water layer type as a new sample, sequentially substituting the new sample into the established preliminary identification model, and if the identification result is consistent with the actual result, indicating that the selection of the identification index and the established preliminary identification model are effective; and if the identification results of the same aquifer water sample are mostly inconsistent with the actual identification results, reselecting a new identification index and establishing a new initial identification model.
S4: establishing a comprehensive-gradual identification method model according to the effective identification index and a Fisher identification method; the comprehensive-step identification method in the step S4 is to comprehensively adopt a characteristic ion contrast method, an ion proportion coefficient method and a Fisher identification method, adopt different methods for identifying different aquifer water sources, and judge the water source types step by step after simple and complex methods; for a single aquifer water sample to be judged, classifying characteristic ions or water samples of which the ion proportionality coefficients are smaller than a determined identification index threshold value into a certain water source type during identification, or classifying the water samples of which the indexes are larger than the threshold value into the certain water source type, stopping identification if a result can be judged, or classifying the water samples into a new sample to be judged, selecting a new identification index and a threshold value thereof, and repeating the steps; and when the type of the water inrush source cannot be judged by using the identification index, a Fisher identification method is adopted. The Fisher recognition method has no specific requirement on the distribution of the population, and is a linear recognition method. It features that high-dimensional data points are projected to low-dimensional space (one-dimensional straight line) so that data points become denser to overcome "dimensional root" caused by high dimension. The principle of projection is to separate the population from the population as much as possible, then determine the recognition analysis function according to the principle of maximum inter-class distance and minimum intra-class distance, and further classify and recognize the new samples. The method can realize the rapid identification of the type of the complex multi-aquifer mine effluent water source.
The comprehensive-step-by-step identification method comprises the following specific steps:
step one, for a certain single aquifer water sample to be judged, judging whether the water sample is AoHu water or Duhui water by adopting a characteristic ion comparison method, if the result can be judged, stopping identification, otherwise, performing identification in the step two;
secondly, judging whether the water sample is cladding body water or coal series water by adopting an ion proportionality coefficient method;
thirdly, on the basis of the second step, if the water sample is judged to be the cladding body water, judging that the water sample is cladding body gneiss water and cladding body cold grey water by adopting a characteristic ion comparison method; and if the water sample is judged to be coal-based water, judging that the water sample is sandstone water or Taigrey water by adopting a Fisher recognition method.
And judging the single aquifer water sample to be judged into one of secondary water, overlying body gneiss water, overlying body cold grey water, sandstone water, taigrey water and Ordovician grey water through the three steps.
S5: and determining the identification indexes of the water sample to be judged, and sequentially judging through a comprehensive-step identification method model to identify the type of the water source.
The specific operation method of the comprehensive-step-by-step identification method comprises the following steps: in the water sample after data verification and abnormal rejection: first, when TDS>4300mg/L、Cl - >2000mg/L、Na + +K + >When 1350mg/L, the water sample is the Aohui water, and the identification is finished; otherwise, further identifying when TDS<440mg/L、Cl - <50mg/L、 Na + +K + <When the water sample is 80mg/L, the water sample is secondary water, and the identification is finished; otherwise, further identification is carried out when gamma Cl - / γCa 2+ <5.2, the water sample is the cover-pushing body water, and is further identified when TDS<1360mg/L、Cl - <474mg/L, the water sample is the water of the coating body water sheet gneiss, and the identification is finished; otherwise, the water sample is the push-coating body water cold grey water, and the identification is finished; when gamma Cl is present - /γCa 2+ >And 5.2, identifying the sandstone water or the Taigray water by using the Fisher identification method, wherein the water sample is coal-based water, and the identification is finished at the moment.
In the embodiment, a characteristic ion contrast method, an ion proportion coefficient method, a Fisher identification method and the like are comprehensively adopted, different methods are adopted for identifying water sources of different aquifers, and the mine water source types with complicated hydrogeological conditions and more water-filled aquifers can be gradually judged after the methods are simple and complex. In the embodiment, the identification indexes adopted by the identification process in the unknown water sample are directly detected, and the identification indexes can be substituted into the established comprehensive-step identification method model to identify the water source. Therefore, the rapidness and the accuracy of identifying the water bursting source of the mine with complex multi-aquifer can be improved.
Next, aquifer water sample data such as new boundary secondary water, tectorial membrane gneiss water, tectorial membrane cold grey water, sandstone water, taigrey water and the like in the new collection secondary water quality standing book are selected as databases to concretely explain the design scheme and the theoretical basis of the embodiment.
Data inspection and exception sample handling
When a single aquifer characteristic index judgment is carried out, the water chemical component of the water sample to be analyzed is required to be a typical representation of the aquifer. Therefore, before characteristic index analysis and water source identification model establishment, each aquifer water sample needs to be inspected and abnormal sample treatment is carried out.
The sample data inspection method is anion-cation balance inspection, due to Na + +K + 、Ca 2+ 、Mg 2+ 、Cl - 、 HCO 3 - 、SO 4 2- The indexes are measured values, so that the error is required to be controlled within +/-5%, namely, qualified water samples with the testing error within the range are used as basic data for further analysis.
Although the water sample is qualified through anion and cation inspection, a few water samples have certain problems in a certain link in the processes of water sample collection, horizon determination, test analysis and the like, so that the water sample cannot truly reflect the water quality characteristics of a certain aquifer, is represented as an outlier sample in the water quality analysis process, and needs to be screened and rejected. The following will illustrate the screening and processing method of abnormal samples by taking Ordovician gray water sample as an example, as shown in FIGS. 2-4.
To 17 collected Oreochromis water samples TDS and Na + +K + 、Ca 2+ 、Cl - The ions were analyzed in comparison, as shown in FIG. 2. The indexes of the three water samples marked as O-03, O-04 and O-05 are obviously different from the other 14 water samples. In order to further judge the relation between the three water samples of the Ordovician gray water No. O-04, O-05 and O-06 and the rest of the water samples, the Ordovician gray water samples are subjected to system cluster analysis. Adopting TDS and Na + +K + 、Ca 2+ 、Mg 2+ 、Cl - 、HCO 3 - 、 SO 4 2- The 7 indexes are used as variables, and the clustering result is shown in FIG. 3. As can be seen from FIG. 3, the differences among the samples except the samples numbered O-03, O-04 and O-05 are the smallest and most representative of the properties of the Aohu waterThe distance between O-03, O-04, O-05 and the rest of the samples is the largest, and the classification is completed at the end. The boxplot is a statistical graph for describing data distribution, and can be used for expressing descriptive statistics such as median, 4 quantile and extreme value of observed data and observing the distribution condition of variable values from a visual angle. 17 water samples TDS and Na of Aohui water + +K + 、Ca 2+ 、Mg 2+ 、Cl - The box diagram is shown in fig. 4. As can be seen from the figure, the 5 water samples marked with the indexes of O-3, O-4 and O-5 also show significant differences from other water samples.
Therefore, three samples O-03, O-04 and O-05 in 17 samples of the Aohu water can be regarded as abnormal samples which cannot represent the typical characteristics of the Aohu water, and are removed in the establishment of a subsequent water sample type identification model and do not participate in modeling as standard water samples.
And (3) comprehensively screening and eliminating outlier samples from the other 5 single aquifer water samples such as the second aquifer water of the new kingdom, the tectorial membrane schist water, the tectorial membrane hanwu system water, the sandstone water, the taigrey water and the like by adopting the method or other methods. And the data subjected to data inspection and abnormal sample removal is used as a modeling sample.
Identification index and threshold determination
Based on the above analysis, the dihydrate and Ordovician water are in TDS and Cl - And Na + +K + The content has typical characteristics, and can be effectively distinguished from other aquifer water samples. Therefore, the above 3 indexes were selected as the identification indexes of the dihydrate and olanzapine water, and the content ranges thereof are shown in table 1.
TABLE 1 content ranges of characteristic indexes of hydrous and olanzapine
(2) Push-coated body water and coal-series water
The types of the conventional ions and water quality of the cladding body gneiss water and the cladding body cold grey water are similar, and the types of the conventional ions and water quality of the sandstone water and the cladding body taigrey water are similar. The gunite water and the cold grey water of the coating body are firstly classified into one type and are collectively called as coating body water; sandstone water and teragrey water are classified into a group, which is collectively called coal-based water. The data of the modeling samples show that the characteristic ion proportionality coefficients (gamma Cl) of the tectorial water and the coal-series water - /γCa 2+ ) The differences were significant (table 2) and could be used as a differential marker for rebinned inferred waters and coal-based waters.
TABLE 2 range of characteristic ion proportionality coefficients for overburden push water and coal-based water
(3) Gunite water and cold grey water of push-coated body
In TDS and Cl, the water of the jacket body gneiss and the water of the jacket body Hanwu system - The content has typical characteristics and can be used as a mark for distinguishing the two types of aquifer water.
TABLE 3 Water characteristic index content Range of the coating
(4) Sandstone water and taigrey water
The sandstone water and the teragrey water have no obvious difference in characteristic example indexes and ion proportionality coefficients, and the conventional means cannot distinguish the sandstone water and the teragrey water only from the two aspects. For this purpose, TDS and Ca are selected comprehensively 2+ 、Mg 2+ 、Na + +K + 、 HCO 3 - 、SO 4 2- 、Cl - And (5) taking the 7 indexes as identification factors for classification between the sandstone water and the taigrey water, and judging the type of the water source by adopting a Fisher identification method.
The Fisher identification method has no specific requirement on the distribution of the population and is a linear identification method. It features that high-dimensional data points are projected to low-dimensional space (one-dimensional straight line) so that data points become denser to overcome "dimensional root" caused by high dimension. The principle of projection is to separate the population from the population as much as possible, then determine the recognition analysis function according to the principle of maximum inter-class distance and minimum intra-class distance, and further classify and recognize the new samples.
Identification method selection and identification step
The new method for identifying the type of the water source of the second-ore complex multi-aquifer adopts a comprehensive-step identification method.
The method comprises the following specific steps:
firstly, judging whether a certain single aquifer water sample to be judged is the Aohu water or the second aquifer water by adopting a characteristic ion comparison method, if the result can be judged, stopping the identification, otherwise, carrying out the second step of identification.
And secondly, judging whether the water sample is the cladding body water or the coal-series water by adopting an ion proportionality coefficient method.
Thirdly, on the basis of the second step, if the water sample is judged to be the cladding body water, judging that the water sample is cladding body gneiss water and cladding body cold grey water by adopting a characteristic ion comparison method; and if the water sample is judged to be coal-series water, judging that the water sample is sandstone water and Taigrey water by adopting a Fisher identification method.
And judging the single aquifer water sample to be judged into one of secondary water, overlying body gneiss water, overlying body cold grey water, sandstone water, taigrey water and Ordovician grey water through the three steps.
Establishment of 'comprehensive-step-by-step recognition' model
For a single aquifer water sample to be judged, after data inspection and abnormal sample processing, identification indexes and limit values, identification method selection and identification step determination are carried out, an established identification model and a process are shown in fig. 3.
Verification of recognition effects
Judging the effectiveness of the identification index by adopting back-substitution test, taking a water sample with a known water layer type as a new sample, sequentially substituting the new sample into the established identification model, and if the identification result is consistent with the reality, indicating that the selection of the identification index and the established 'comprehensive-gradual identification method' model are effective; and if the identification results of the same aquifer water sample are mostly inconsistent with the actual identification results, a new identification index is reselected, and a new water source identification model is established.
After the comprehensive-gradual identification method model is determined, the unknown aquifer water sample is substituted into the established model after the content of the identification index is tested, and the aquifer type of the unknown water sample can be obtained.
The above description is only exemplary of the present invention and should not be taken as limiting the invention, as any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.