Analysis method and device for water accumulation area of coal mine goaf and storage medium

Abstract

The invention relates to the technical field of address exploration, and discloses a coal goaf water accumulation area analysis method and device and a storage medium. The method comprises the following steps: acquiring the resistivity of a goaf, the resistivity below the height of three zones above the goaf, an induced voltage value of the goaf and a multi-track curve of the goaf by using a transient electromagnetic method; obtaining the ratio of the induced voltage values of the early stage track and the late stage track of the goaf according to the multi-track curve of the goaf; processing the specific value of the induced voltage values of the early stage channel and the late stage channel of the goaf obtained by the resistivity of the goaf at each measuring point, the resistivity below the height of three zones above the goaf, the induced voltage values of the goaf and the multi-measuring-channel curve of the goaf by utilizing a principal component analysis method; synthesizing the four parameters at each measuring point into an evaluation parameter S; and establishing a slice map of the evaluation parameter S of the goaf. The method can improve the explanation accuracy of the goaf ponding area range.

Description

Technical Field

[0001] The invention relates to the technical field of geological exploration, and in particular to an analysis method, a device and a storage medium for water accumulation in a coal mine goaf. Background Art

[0002] After the coal seams in the strata are mined, certain spaces will be formed between the underground rock layers. These spaces are called goafs. The goafs with lower terrain are prone to water accumulation, which will bring safety hazards to the mining process of coal mines. Therefore, it is necessary to analyze and detect the water accumulation areas in the goafs, and drain the water in the water accumulation areas. In the prior art, the goafs and water accumulation areas in the goafs in the coal seams are generally detected by transient electromagnetic method, and the areas with resistivity significantly lower than other positions in the same layer are interpreted as water accumulation areas in the goafs. However, due to the multi-solution nature of the transient electromagnetic method, the error of analyzing and interpreting the water accumulation areas in the goafs using only one attribute, resistivity, is large. For this reason, it is necessary to propose a new method for analyzing and interpreting the water accumulation areas in the goafs of coal mines to improve the accuracy of exploration of water accumulation areas in the goafs. Summary of the invention

[0003] In order to overcome the deficiencies of the prior art, the purpose of the present invention is to provide an analysis method for water accumulation areas in coal mine goafs, which can improve the accuracy of analysis and interpretation of water accumulation areas in goafs and improve the accuracy of exploration of water accumulation areas in goafs.

[0004] To solve the above problems, the technical solution adopted by the present invention is as follows: A method for analyzing water accumulation areas in coal mine goafs, comprising the following steps:

[0005] The transient electromagnetic method is used to obtain the resistivity of the goaf, the resistivity below the three zones above the goaf, the induced voltage value of the goaf and the multi-channel curve of the goaf;

[0006] According to the multi-channel curve of the goaf, the ratio of the induced voltage value of the early channel and the late channel of the goaf is obtained;

[0007] The principal component analysis method is used to process the resistivity of the goaf at each measuring point, the resistivity below the three zones above the goaf, the induced voltage value of the goaf, and the ratio of the induced voltage values ​​of the early track and the late track of the goaf obtained by the multi-track curve of the goaf, and the above four parameters at each measuring point are synthesized into an evaluation parameter S;

[0008] Establish a slice diagram of the evaluation parameter S of the goaf.

[0009] Compared with the prior art, the beneficial effect of the present invention lies in: the location of the waterlogged area in the goaf is explained by using an evaluation parameter S formed by combining four parameters, namely the resistivity of the goaf, the resistivity below the three zones above the goaf, the induced voltage value of the goaf and the ratio of the induced voltage values ​​of the early channel and the late channel of the goaf obtained by the multi-channel curve of the goaf. Since the four parameters have obvious characteristics in the waterlogged area in the goaf, the accuracy of the interpretation of the waterlogged area can be improved by using the evaluation parameter S formed by combining the above four parameters to explain the waterlogged area in the goaf, thereby more accurately confirming the location of the waterlogged area.

[0010] The above-mentioned method for analyzing the water accumulation area in the goaf of a coal mine, wherein the resistivity of the goaf at each measuring point, the resistivity below the three zones above the goaf, the induced voltage value of the goaf, and the multi-channel curve of the goaf are processed by the principal component analysis method to obtain the ratio of the induced voltage values ​​of the early channel and the late channel of the goaf, and the above-mentioned four parameters at each measuring point are synthesized into an evaluation parameter S, which includes the following steps:

[0011] Normalize the resistivity of the goaf, the resistivity below the three zones above the goaf, the induced voltage value of the goaf, and the ratio of the induced voltage values ​​of the early trace to the late trace of the goaf;

[0012] The contribution of the above four normalized parameters at each measuring point is calculated respectively, the principal components are selected according to the contribution of the four parameters, and the principal components are combined into an evaluation parameter S according to the contribution.

[0013] In the above-mentioned method for analyzing water accumulation areas in coal mine goafs, in the step of normalizing the resistivity of the goaf, the resistivity below the three zones above the goaf, the induced voltage value of the goaf, and the ratio of the induced voltage values ​​of the early channel and the late channel of the goaf, normalization is performed by respectively calculating the Z scores of the above four parameters.

[0014] The above-mentioned method for analyzing the water accumulation area in the coal mine goaf, the step of respectively calculating the contribution of the above-mentioned four parameters normalized at each measuring point, calculating the principal components according to the contribution size of the four parameters, and merging the principal components into an evaluation parameter S according to the contribution size includes:

[0015] Establish a covariance matrix for the four normalized parameters of one of the measurement points;

[0016] Solve the characteristic equation |ρ xy -λI|=0, we obtain the eigenvalues ​​λ and eigenunit vectors I of the four parameters, where ρ xy is the item in the covariance matrix;

[0017] The four parameters are sorted from large to small according to the eigenvalue λ, and 2 to 3 parameters whose cumulative variance accounts for more than 70% or parameters whose eigenvalue is greater than 1 are selected as the principal components;

[0018] The principal components are combined into an evaluation parameter S according to their contribution size according to the following formula:

[0019]

[0020] Where P i is the variance proportion of the i-th principal component, S i is the value of the i-th principal component, and P is the total variance share of the n principal components.

[0021] A storage medium is used to store a computer program, which, when called and executed by a processor, implements the above-mentioned method for analyzing water accumulation areas in coal mine goafs.

[0022] A device for analyzing water accumulation areas in coal mine goafs, comprising a processor and the above-mentioned storage medium, wherein the storage medium is electrically connected to the processor, and the processor can implement the above-mentioned method for analyzing water accumulation areas in coal mine goafs by calling a computer program stored in the storage medium.

[0023] An analysis device for water accumulation areas in goafs of coal mines, comprising: an acquisition module, used to acquire the resistivity of the goaf, the resistivity below the height of three zones above the goaf, the induced voltage value of the goaf and the multi-channel curve of the goaf, and obtain the ratio of the induced voltage value of the early channel to the late channel of the goaf according to the multi-channel curve; an evaluation parameter calculation module, used to process the resistivity of the goaf at each measuring point, the resistivity below the height of three zones above the goaf, the induced voltage value of the goaf and the ratio of the induced voltage value of the early channel to the late channel of the goaf obtained by the multi-channel curve of the goaf using a principal component analysis method, and synthesize the above four parameters at each measuring point into an evaluation parameter S; an evaluation parameter graph analysis module, used to establish a slice graph of the evaluation parameter S of the goaf according to the evaluation parameter S of each measuring point calculated by the evaluation parameter calculation module.

[0024] The above-mentioned analysis device for water accumulation area in coal mine goaf, the evaluation parameter calculation module includes: a normalization module, which is used to normalize the resistivity of the goaf, the resistivity below the three zones above the goaf, the induced voltage value of the goaf and the ratio of the induced voltage values ​​of the early track and the late track of the goaf.

[0025] The above-mentioned analysis device for water accumulation area in coal mine goaf, the evaluation parameter calculation module also includes: a covariance calculation module, used to establish the covariance matrix of the four parameters of each measuring point; a characteristic calculation module, used to calculate the eigenvalue and characteristic unit vector of the covariance matrix of each measuring point; a principal component merging module, used to sort the four parameters from large to small according to the eigenvalue, select 2 to 3 parameters with a cumulative variance share of more than 70% or parameters with an eigenvalue greater than 1 as principal components, and merge the principal components of each measuring point into the evaluation parameter S of the measuring point according to the contribution size.

[0026] In the above-mentioned analysis device for water accumulation area in coal mine goaf, the normalization module is used to perform normalization processing by respectively calculating the Z scores of four parameters.

[0027] The present invention is further described in detail below in conjunction with the accompanying drawings and specific embodiments. BRIEF DESCRIPTION OF THE DRAWINGS

[0028] Figure 1 This is a schematic diagram of the goaf and the rock formation above it;

[0029] Figure 2 This is the apparent resistivity cross section of a goaf;

[0030] Figure 3 This is a multi-channel curve diagram of a goaf area;

[0031] Figure 4 This is a transient electromagnetic layer resistivity slice diagram of a coal mine goaf;

[0032] Figure 5 This is a slice diagram of the evaluation parameter S of a coal mine goaf. DETAILED DESCRIPTION

[0033] The following describes an embodiment of the present invention in detail. The embodiment of the present invention provides a method for analyzing water accumulation areas in coal mine goafs, comprising the following steps:

[0034] The transient electromagnetic method is used to obtain the resistivity of the goaf, the resistivity below the three zones above the goaf, the induced voltage value of the goaf and the multi-channel curve of the goaf;

[0035] According to the multi-channel curve of the goaf, the ratio of the induced voltage value of the early channel and the late channel of the goaf is obtained;

[0036] The principal component analysis method is used to process the resistivity of the goaf at each measuring point, the resistivity below the three zones above the goaf, the induced voltage value of the goaf, and the ratio of the induced voltage values ​​of the early track and the late track of the goaf obtained by the multi-track curve of the goaf, and the above four parameters at each measuring point are synthesized into an evaluation parameter S;

[0037] Establish a slice diagram of the evaluation parameter S of the goaf.

[0038] Reference Figure 1 The three zones refer to the collapse of the rock layer above the goaf under the action of gravity, which causes the overlying rock mass of the coal seam to lose its original equilibrium state and undergo a certain degree of rock movement, destroying the integrity and continuity of the rock, causing the rock layer to break and a large number of gaps and cracks to appear, forming a caving zone, a fracture zone and a bending zone, referred to as the three zones. When the gaps in the goaf are filled with water or mud, the resistivity of the area will be significantly lower than the resistivity of the surrounding intact rock, and it will also show a certain low-resistance feature. This feature can be used to divide the scope of the water accumulation area in the goaf. Figure 2 As shown in the figure, there is water accumulation in the goaf of the 9# coal seam in a coal mine. The apparent resistivity at stakes 600 to 1100m near the 9# coal seam is significantly lower than that at other locations in the same layer, and there is an obvious low-resistance area above it. And the mining data of the coal mine shows that the roof collapsed after the goaf, resulting in the existence of three zones at a certain height above the goaf. The increase in rock stratum cracks has led to an increase in water content, and a "shadow zone" with obvious low-resistance characteristics has appeared above, which fully demonstrates that the resistivity above the goaf to the height below the three zones can be used as an obvious feature to explain the scope of the water accumulation area in the goaf.

[0039] According to the above characteristics, when interpreting the goaf, transient electromagnetic apparent resistivity slices can be made in the coal seam and the three zones above the coal seam according to the height of the "three zones" above the coal seam goaf. If the resistivity in the resistivity slices along the layer in the coal seam and above the coal seam is less than the resistivity of the surrounding rock, it is interpreted as a water-containing goaf. If the resistivity in the resistivity slices along the layer in the coal seam is less than the resistivity of the surrounding rock, and the resistivity of the coal seam does not change much compared with the surrounding rock, it is interpreted as a water-free goaf. In this way, by increasing the three-zone resistivity slices of the goaf to interpret the goaf, the accuracy of the interpretation of the coal seam goaf can be effectively improved. Figure 2 Taking coal mines as an example, the impact height of the "three zones" is about 50m. Therefore, in data processing, the layer resistivity of the 9# coal seam and 20m above the 9# coal seam is sliced.

[0040] Since the resistivity of the water-bearing goaf is lower than that of the surrounding rock, the induced voltage will increase under the transient electromagnetic fixed source loop device, which is manifested as a convex curve in the multi-channel curve, and it also affects the late channel data, which is manifested as the greater the change in the late channel data compared with the normal data of the same channel. Therefore, the ratio of the induced voltage of the early channel to the late channel in the multi-channel curve can also be used as a feature to reflect the water accumulation in the goaf. Figure 3As shown, it is a transient electromagnetic multi-channel curve diagram of a known goaf area, among which pile numbers 2960-3040 are known goaf water accumulation areas. The early channel curve does not change much, while the induced voltage of the late channel curve changes significantly, which fully demonstrates that the ratio of the induced voltage values ​​of the early channel and the late channel can be used to reflect the water accumulation in the goaf, and the area with a low ratio is the water accumulation area.

[0041] To sum up, the analysis method of the water accumulation area in the goaf of a coal mine in the embodiment of the present invention uses the principal component analysis method to integrate the resistivity of the goaf, the induced voltage value of the goaf and the above two parameters which also have obvious characteristics in the water accumulation area of ​​the goaf into an evaluation parameter S to explain the water accumulation area in the goaf. Compared with the method of only using the resistivity of the goaf to analyze and explain the scope of the water accumulation area in the goaf, it has higher accuracy. The scope of the water accumulation area in the goaf can be better determined through the slice diagram of the evaluation parameter S, thereby more accurately confirming the location of the water accumulation area in the goaf and improving the safety during coal mining.

[0042] Since the dimensions of the four parameters are different, they usually need to be normalized before being combined using the principal component analysis method. The principal component analysis method includes the following steps:

[0043] Normalize the resistivity of the goaf, the resistivity below the three zones above the goaf, the induced voltage value of the goaf, and the ratio of the induced voltage values ​​of the early trace to the late trace of the goaf;

[0044] The contribution of the above four normalized parameters at each measuring point is calculated respectively, and the principal components are generated according to the contribution sizes of the four parameters, and the principal components are combined into an evaluation parameter S according to the contribution sizes.

[0045] It is understandable that the normalization process may adopt a minimum-maximum normalization method or a Z-score normalization method, etc. In this embodiment, the Z-score is used to normalize the four parameters, and the calculation formula of the Z-score is shown as follows:

[0046]

[0047] Where x′ i is the standardized value of the parameter, x i is the original value of the standardized parameter, σ is the standard deviation of the standardized parameter, is the mean value of the standardized parameter. Z-score normalization can preserve the data distribution characteristics while making each transformed data have the same mean and standard deviation.

[0048] In this embodiment, after the data is normalized, the covariance matrix of the four parameters of each measuring point is established, the characteristic equation is solved to obtain the eigenvalue of each parameter, and the contribution of each parameter is calculated by calculating the ratio of the eigenvalue of each parameter to the sum of the eigenvalues.

[0049] For the four parameters, a 4×4 covariance matrix can be formed, and the eigenvalues ​​λ and eigenunit vectors I of the four parameters can be obtained by solving the following formula:

[0050] |ρ xy -λI|=0 (2)

[0051] Then, the four eigenvalues ​​of the four parameters are sorted from large to small, and 2 to 3 samples whose variance accounts for 70% or samples whose eigenvalues ​​are greater than 1 are selected as principal components. The obtained 2 to 3 principal components are synthesized into a parameter S according to the variance contribution size using the following formula:

[0052]

[0053] Where P i is the variance proportion of the i-th principal component, S i is the value of the i-th principal component, and P is the total variance share of the n principal components.

[0054] Finally, the evaluation parameters S of multiple measuring points are plotted into a slice diagram of the evaluation parameter S, and the area with relatively low S value in the diagram can be interpreted as the water accumulation area in the goaf. Figure 4 and Figure 5 As shown in the figure, they are the transient electromagnetic layer resistivity slice diagram and the evaluation parameter S slice diagram of the same coal mine goaf area. Figure 5 contrast Figure 4 It can be seen that in the transient electromagnetic layer resistivity slice diagram, most of the low-resistance areas are distributed in the goaf, but some areas have little correlation with the goaf; while in the slice diagram of the evaluation parameter S, the areas with relatively low S values ​​are all located in the goaf, and have a better correlation with the goaf, indicating that the method of the embodiment of the present invention can more accurately explain the scope of the water accumulation area in the goaf.

[0055] Based on the same inventive concept, an embodiment of the present invention further provides a computer-readable storage medium, in which a computer program is stored. When the computer program is executed by a processor, the above-mentioned method for analyzing water accumulation areas in coal mine goafs can be implemented.

[0056] In some possible embodiments, various aspects of the method for analyzing water accumulation areas in coal mine goafs provided by the present invention can also be implemented in the form of a program product, which includes program code. When the program product is run on an apparatus, the program code is used to enable the control device to execute the steps of the method for analyzing water accumulation areas in coal mine goafs according to various exemplary embodiments of the present application described above in this specification.

[0057] Based on the same inventive concept, an embodiment of the present invention also provides an analysis device for implementing the above-mentioned analysis method for coal mine waterlogging areas, including a processor and the above-mentioned storage medium, the storage medium is electrically connected to the processor, and the processor is used to execute the computer program stored on the storage medium to implement the above-mentioned analysis method for coal mine waterlogging areas.

[0058] In one possible design, the processor may include one or more processing units, and the processor may integrate an application processor and a modem processor, wherein the application processor mainly processes an operating system, a user interface, and application programs, and the modem processor mainly processes wireless communications. It is understandable that the modem processor may not be integrated into the processor. In some embodiments, the processor and the memory may be implemented on the same chip, and in some embodiments, they may also be implemented separately on separate chips.

[0059] The processor may be a general-purpose processor, such as a central processing unit (CPU), a digital signal processor, an application-specific integrated circuit, a field programmable gate array or other programmable logic device, a discrete gate or transistor logic device, a discrete hardware component, and may implement or execute the methods, steps, and logic block diagrams disclosed in the embodiments of the present application. A general-purpose processor may be a microprocessor or any conventional processor, etc. The steps of the output method of the landing area disclosed in the embodiments of the present application may be directly embodied as being executed by a hardware processor, or may be executed by a combination of hardware and software modules in the processor.

[0060] As a non-volatile computer-readable storage medium, the memory can be used to store non-volatile software programs, non-volatile computer executable programs and modules. The memory can include at least one type of storage medium, for example, it can include flash memory, hard disk, multimedia card, card-type memory, random access memory (Random Access Memory, RAM), static random access memory (Static RandomAccess Memory, SRAM), programmable read-only memory (Programmable Read Only Memory, PROM), read-only memory (Read Only Memory, ROM), electrically erasable programmable read-only memory (Electrically Erasable Programmable Read-Only Memory, EEPROM), magnetic memory, disk, CD, etc. The memory is any other medium that can be used to carry or store the desired program code in the form of instructions or data structures and can be accessed by a computer, but is not limited thereto. The memory in the embodiment of the present application can also be a circuit or any other device that can realize a storage function, for storing program instructions and / or data.

[0061] By designing and programming the processor, the code corresponding to the method for analyzing the water accumulation area in the coal mine goaf area described in the above embodiment can be fixed into the chip, so that the chip can execute the steps of the method for analyzing the water accumulation area in the coal mine goaf area of ​​the embodiment shown in the present invention when running. How to design and program the processor is a technology known to those skilled in the art and will not be described in detail here.

[0062] Based on the same inventive concept, an embodiment of the present invention also provides another analysis device for coal mine goaf water accumulation area, including an acquisition module, an evaluation parameter calculation module and an evaluation parameter map analysis module. Among them, the acquisition module is used to obtain the resistivity of the goaf, the resistivity below the three zones above the goaf, the induced voltage value of the goaf and the multi-channel curve of the goaf, and obtain the ratio of the induced voltage value of the early channel to the late channel of the goaf according to the multi-channel curve. The evaluation parameter calculation module is used to use the principal component analysis method to process the resistivity of the goaf at each measuring point, the resistivity below the three zones above the goaf, the induced voltage value of the goaf and the ratio of the induced voltage value of the early channel to the late channel of the goaf obtained by the multi-channel curve of the goaf, and synthesize the above four parameters at each measuring point into an evaluation parameter S. The evaluation parameter map analysis module is used to establish a slice map of the evaluation parameter S of the goaf according to the evaluation parameter S of each measuring point calculated by the evaluation parameter calculation module.

[0063] In some embodiments, the evaluation parameter calculation module includes a normalization module, a covariance calculation module, a feature calculation module and a principal component merging module. The normalization module is used to normalize the resistivity of the goaf, the resistivity below the three zones above the goaf, the induced voltage value of the goaf, and the ratio of the induced voltage values ​​of the early track to the late track of the goaf. The covariance calculation module is used to establish the covariance matrix of the four parameters of each measuring point. The feature calculation module is used to calculate the eigenvalues ​​and eigenunit vectors of the covariance matrix of each measuring point. The principal component merging module is used to sort the four parameters from large to small according to the eigenvalue, select 2 to 3 parameters with a cumulative variance share of more than 70% or parameters with an eigenvalue greater than 1 as principal components, and merge the principal components of each measuring point into the evaluation parameter S of the measuring point according to the contribution size.

[0064] The present application is described with reference to the flowcharts and / or block diagrams of the methods, apparatus (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of the processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0065] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0066] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0067] In the description of the present invention, "several" means one or more, "more" means two or more, "greater than", "less than", "exceed" etc. are understood to exclude the number itself, and "above", "below", "within" etc. are understood to include the number itself. If there is a description of "first" or "second", it is only used to distinguish the technical features, and cannot be understood as indicating or implying relative importance or implicitly indicating the number of the indicated technical features or implicitly indicating the order of the indicated technical features.

[0068] In the description of the present invention, unless otherwise clearly defined, terms such as setting, installing, connecting, etc. should be understood in a broad sense, and technicians in the relevant technical field can reasonably determine the specific meanings of the above terms in the present invention in combination with the specific content of the technical solution.

[0069] The above-mentioned embodiments are only preferred embodiments of the present invention and cannot be used to limit the scope of protection of the present invention. Any non-substantial changes and substitutions made by technicians in this field on the basis of the present invention shall fall within the scope of protection required by the present invention.

Claims

1. A method for analyzing water accumulation areas in coal mine goafs, characterized in that: The steps include: The transient electromagnetic method is used to obtain the resistivity of the goaf, the resistivity below the three zones above the goaf, the induced voltage value of the goaf and the multi-channel curve of the goaf; According to the multi-channel curve of the goaf, the ratio of the induced voltage value of the early channel and the late channel of the goaf is obtained; The principal component analysis method is used to process the resistivity of the goaf at each measuring point, the resistivity below the three zones above the goaf, the induced voltage value of the goaf, and the ratio of the induced voltage values ​​of the early track and the late track of the goaf obtained by the multi-track curve of the goaf, and the above four parameters at each measuring point are synthesized into an evaluation parameter S; Establish a slice diagram of the evaluation parameter S of the goaf.

2. The method for analyzing water accumulation in coal mine goaf according to claim 1, characterized in that: The step of using the principal component analysis method to process the resistivity of the goaf at each measuring point, the resistivity below the three zones above the goaf, the induced voltage value of the goaf, and the ratio of the induced voltage values ​​of the early track and the late track of the goaf obtained by the multi-track curve of the goaf, and synthesizing the above four parameters at each measuring point into an evaluation parameter S includes: Normalize the resistivity of the goaf, the resistivity below the three zones above the goaf, the induced voltage value of the goaf, and the ratio of the induced voltage values ​​of the early trace to the late trace of the goaf; The contribution of the above four normalized parameters at each measuring point is calculated respectively, the principal components are selected according to the contribution of the four parameters, and the principal components are combined into an evaluation parameter S according to the contribution.

3. The method for analyzing water accumulation in coal mine goaf according to claim 2, characterized in that: In the step of normalizing the resistivity of the goaf, the resistivity below the three zones above the goaf, the induced voltage value of the goaf, and the ratio of the induced voltage values ​​of the early track to the late track of the goaf, normalization is performed by respectively calculating the Z scores of the above four parameters.

4. The method for analyzing water accumulation in coal mine goaf according to claim 2, characterized in that: The step of respectively calculating the contribution of the four normalized parameters at each measuring point, calculating the principal components according to the contribution sizes of the four parameters, and merging the principal components into an evaluation parameter S according to the contribution sizes comprises: Establish a covariance matrix for the four normalized parameters of one of the measurement points; Solve the characteristic equation |ρ xy -λI|=0, we obtain the eigenvalues ​​λ and eigenunit vectors I of the four parameters, where ρ xy is the item in the covariance matrix; The four parameters are sorted from large to small according to the eigenvalue λ, and 2 to 3 parameters whose cumulative variance accounts for more than 70% or parameters whose eigenvalue is greater than 1 are selected as the principal components; The principal components are combined into an evaluation parameter S according to their contribution size according to the following formula: Where P i is the variance proportion of the i-th principal component, S i is the value of the i-th principal component, and P is the total variance share of the n principal components.

5. A storage medium for storing a computer program, characterized in that: When the computer program is called and executed by the processor, the method for analyzing water accumulation areas in coal mine goafs according to any one of claims 1 to 4 is implemented.

6. An analysis device for water accumulation in coal mine goaf, characterized in that: It comprises a processor and a storage medium according to claim 5, wherein the storage medium is electrically connected to the processor, and the processor can implement the analysis method of water accumulation area in coal mine goaf according to any one of claims 1 to 4 by calling a computer program stored in the storage medium.

7. An analysis device for water accumulation in coal mine goaf, characterized in that: include: An acquisition module is used to acquire the resistivity of the goaf, the resistivity below the three zones above the goaf, the induced voltage value of the goaf and the multi-channel curve of the goaf, and obtain the ratio of the induced voltage value of the early channel to the late channel of the goaf according to the multi-channel curve; An evaluation parameter calculation module is used to process the resistivity of the goaf at each measuring point, the resistivity below the three zones above the goaf, the induced voltage value of the goaf, and the ratio of the induced voltage values ​​of the early track and the late track of the goaf obtained by the multi-track curve of the goaf using the principal component analysis method, and synthesize the above four parameters at each measuring point into an evaluation parameter S; The evaluation parameter diagram analysis module is used to establish a slice diagram of the evaluation parameter S of the goaf according to the evaluation parameter S of each measuring point calculated by the evaluation parameter calculation module.

8. The device for analyzing water accumulation in coal mine goaf according to claim 7, characterized in that: The evaluation parameter calculation module includes: The normalization module is used to normalize the resistivity of the goaf, the resistivity below the three zones above the goaf, the induced voltage value of the goaf, and the ratio of the induced voltage values ​​of the early track and the late track of the goaf.

9. The device for analyzing water accumulation in coal mine goaf according to claim 8, characterized in that: The evaluation parameter calculation module also includes: Covariance calculation module, used to establish the covariance matrix of the four parameters of each measuring point; The characteristic calculation module is used to calculate the eigenvalue and characteristic unit vector of the covariance matrix of each measuring point; The principal component merging module is used to sort the four parameters from large to small according to the eigenvalues, select 2 to 3 parameters whose cumulative variance accounts for more than 70% or parameters whose eigenvalues ​​are greater than 1 as principal components, and merge the principal components of each measuring point into the evaluation parameter S of the measuring point according to the contribution size.

10. The device for analyzing water accumulation in coal mine goaf according to claim 8, characterized in that: The normalization module is used to perform normalization processing by respectively calculating the Z scores of the four parameters.

Images