Coal seam working face floor detection and analysis method

Abstract

The invention discloses a coal seam working face floor detection and analysis method, and relates to the technical field of detection and analysis. Potential change conditions of different positions and different elevations are observed through a parallel electrical method instrument, resistivity distribution conditions of different depths in a working face and a bottom plate of the working face can be obtained through three-dimensional electrical method inversion, and the resistivity distribution conditions provide technical support for objective and accurate geological interpretation. The hydrogeological condition of the working face floor rock stratum can be better understood, so that an important reference basis is provided for safe mining; the thermal effect and the electromagnetic effect existing in the environment are evaluated and corrected, so that the deviation of the measured conductivity value can be accurately corrected, and the accuracy and the reliability of conductivity measurement data are improved; and inputting the corrected conductivity value SR (xi, yi, zi) of each grid point into a three-dimensional model, and marking to construct a conductivity three-dimensional model, so as to visually display the electrical property distribution information of the underground medium.

Description

Technical Field

[0001] The present invention relates to the field of detection and analysis technology, and in particular to a method for detecting and analyzing the bottom plate of a coal seam working face. Background Art

[0002] Coal mining is a high-risk industry. The stability of the working face floor is directly related to the safe production of the mine. Through detection and analysis, the unstable factors and potential safety hazards of the floor can be discovered in time, and corresponding measures can be taken to ensure the safety of miners and the integrity of mine facilities. Therefore, for coal mining, working face floor detection is very necessary.

[0003] Although traditional geological exploration methods such as drilling and geological profiling can provide geological information to a certain extent, they have certain limitations in fully understanding the hydrogeological conditions of the bottom plate. These methods can usually only provide local geological information and it is difficult to cover the hydrogeological conditions of the entire bottom plate rock layer of the working face. In addition, traditional geological exploration methods require a large amount of drilling, measurement and other operations, which are not only time-consuming and labor-intensive, but also costly. Summary of the invention

[0004] Based on this, it is necessary to provide a coal seam working face floor detection and analysis method to address the problems mentioned in the above background technology.

[0005] The purpose of the present invention can be achieved by the following technical solution: A method for detecting and analyzing the floor of a coal seam working face comprises the following steps:

[0006] T1: The parallel electrode intelligent acquisition method and double-lane three-dimensional resistivity imaging are used to detect the resistivity value of each three-dimensional grid to obtain the conductivity measurement value of each three-dimensional grid, and the conductivity measurement value is recorded as R(xi, yi, zi), where (xi, yi, zi) represents the specific position of the grid numbered i in the three-dimensional coordinates; i = 1, 2, 3...I, I is a positive integer, I represents the total number of grids, and i represents any grid number among them;

[0007] T2: Environmental interference assessment and correction, the specific process is as follows:

[0008] T21: The devices in the environment that are turned on are recorded as related devices. Randomly select one of the grids and quantitatively analyze the thermal effect interference and electromagnetic interference generated by the related devices at the grid point to obtain the temperature transfer index and electromagnetic transfer index of each related device in the grid;

[0009] T22: Comprehensively analyze the temperature transfer index of each related device at the grid to obtain the thermal effect interference index Gb corresponding to the grid, thereby obtaining the thermal effect interference index of each grid and recording it as Gb(xi, yi, zi);

[0010] T23: Comprehensively analyze the electromagnetic transfer index of each relevant device at the grid to obtain the electromagnetic effect interference index Cb corresponding to the grid, thereby obtaining the electromagnetic effect interference index of each grid and recording it as Cb(xi, yi, zi);

[0011] T24: The thermal effect interference index Gb(xi, yi, zi) and the electromagnetic effect interference index are recorded as Cb(xi, yi, zi) through the set formula Calculate and obtain the correction index CG (xi, yi, zi), where b4 and b5 are respectively the set proportional coefficients, and γ is the set correction conversion coefficient; multiply the conductivity measurement value obtained by measuring each grid by the corresponding correction index to obtain the corrected conductivity value SR (xi, yi, zi) of the grid point;

[0012] T3: Establish an underground three-dimensional model, input the corrected conductivity value SR (xi, yi, zi) of each grid point into the three-dimensional model and mark it to construct a three-dimensional conductivity model.

[0013] In some embodiments, the specific process of quantitatively analyzing the thermal effect interference and electromagnetic interference generated by the relevant equipment at the grid point is:

[0014] Record the devices in the environment that are turned on as relevant devices, obtain the temperature and electromagnetic parameters of each relevant device, where the electromagnetic parameters include electromagnetic intensity and electromagnetic direction; obtain the electromagnetic direction applied by each grid for measuring the conductivity value, and record it as the reference direction;

[0015] Take any one of the grids, obtain the distance between each relevant device and the grid, and record it as the relevant distance; obtain the ambient temperature and the underground medium type between each relevant device and the grid; set different medium types to correspond to different propagation coefficients, and compare the underground medium type between each relevant device and the grid with all the set underground medium types to match the corresponding propagation coefficient;

[0016] The device temperature ST, ambient temperature HT, related distance A1, electromagnetic intensity P1, and propagation coefficient P2 are combined through the set formula The temperature transfer index and electromagnetic transfer index of each relevant device are calculated, where a1, a2, a3, a4, and a5 are set proportional coefficients, and e is a natural constant.

[0017] In some embodiments, the specific process of comprehensively analyzing the temperature transfer index of each related device at the grid is:

[0018] Compare and analyze the temperature transfer index of each relevant equipment with the set temperature transfer range to divide the temperature transfer index into a primary transfer index, a secondary transfer index and a tertiary transfer index. The number of the primary transfer index, the secondary transfer index and the tertiary transfer index corresponds to the cumulative number of high thermal effect interference, medium thermal effect interference and low thermal effect interference, respectively, and they are recorded as U1, U2 and U3, respectively; sum up the primary transfer index, the secondary transfer index and the tertiary transfer index to obtain the primary transfer total value, the secondary transfer total value and the tertiary transfer total value, and record them as U4, U5 and U6, respectively;

[0019] Substitute the cumulative number of high thermal effect interference U1, the cumulative number of medium thermal effect interference U2, the cumulative number of low thermal effect interference U3, the total value of the first-level transfer U4, the total value of the second-level transfer U5 and the total value of the third-level transfer U6 into the set formula The thermal effect interference index Ub of the grid is calculated, where b1, b2, and b3 are respectively set proportional coefficients, and e is a natural constant.

[0020] In some embodiments, the specific process of comprehensively analyzing the electromagnetic transfer index of each relevant device at the grid is:

[0021] The reference direction of the grid is taken as the positive direction of the horizontal axis, and the reverse direction of the reference direction is taken as the negative direction of the horizontal axis; the electromagnetic direction corresponding to the electromagnetic transfer index of each relevant device at the grid is decomposed on the horizontal axis to obtain the decomposed electromagnetic direction; thus, the electromagnetic transfer index of each relevant device at the grid and its corresponding decomposed electromagnetic direction can be obtained, and the electromagnetic transfer index of each relevant device at the grid is recorded as APj, where j = 1, 2, 3...J, J is a positive integer J represents the total number of relevant devices at the grid, and j represents the serial number of any relevant device at the grid; by setting the formula The electromagnetic effect interference index Cb of the grid is calculated, where n takes an odd or even value. When the decomposed electromagnetic direction is a positive direction, n takes an even value; when the decomposed electromagnetic direction is a negative direction, n takes an odd value.

[0022] In some embodiments, the implementation steps of the parallel electrode intelligent acquisition method are:

[0023] A symmetrical quadrupole device is arranged on the surface of the earth to monitor the potential difference of each three-dimensional grid, wherein the symmetrical quadrupole device includes two power supply electrodes, namely electrode A and electrode B, and two measuring electrodes, namely electrode M and electrode N. The observed potential difference at electrode M and electrode N is calculated using the uniform earth resistivity calculation expression to obtain the apparent resistivity;

[0024] Each electrode is equipped with an A / D converter and performs parallel sampling. It maintains real-time contact with the host according to the network protocol. When receiving the power supply status command, the electrode sampling part is disconnected to keep the electrode in the AB power supply state. Otherwise, it always works in the voltage sampling state and sends the measurement data back to the host in real time through the communication line; and the conductivity measurement value of each three-dimensional grid is recorded as R (xi, yi, zi).

[0025] In some embodiments, the implementation steps of double-lane three-dimensional resistivity imaging are:

[0026] The expression of the three-dimensional inverse resistivity problem is: Δd = GΔm, where G is the Jacobi matrix; Δd is the observed data, d and the forward theoretical value d 0 The residual vector of ; Δm is the modification vector of the initial model m;

[0027] Add smooth constraints to the least squares criterion, invert and obtain a smooth model to improve the stability of the solution. The algorithm for solving the model modification is: (G T G+λC T C) Δm=G T Δd, where C is the model smoothing matrix, and the conductivity measurement value of each three-dimensional grid is obtained by solving the Jacobi matrix G and calculating the large matrix inverse, which is recorded as R(xi, yi, zi).

[0028] Compared with the prior art, the present invention has the following beneficial effects:

[0029] 1. Intelligently collect electrical data through parallel electrodes to ensure the synchronization of potential measurement and avoid the interference of data measured at different times. Such data bodies are particularly suitable for using full-space three-dimensional resistivity inversion technology to obtain the electrical distribution between working surfaces;

[0030] 2. Through the parallel electrical instrument to observe the potential changes at different locations and elevations, and through three-dimensional electrical inversion, the resistivity distribution at different depths in the working face and its floor can be obtained. These resistivity distributions provide technical support for objective and accurate geological interpretation, which helps to better understand the hydrogeological conditions of the rock formations at the bottom of the working face, thus providing an important reference for safe mining;

[0031] 3. By evaluating and correcting the thermal and electromagnetic effects in the environment, the deviation of the measured conductivity value can be accurately corrected, the accuracy and reliability of the conductivity measurement data are improved, and a more reliable and accurate basis is provided for subsequent data analysis and application;

[0032] 4. By inputting the corrected conductivity value SR(xi, yi, zi) of each grid point into the three-dimensional model and marking it, a three-dimensional conductivity model is constructed to visualize the distribution information of the electrical properties of the underground medium, thereby providing scientific basis and data support for geological exploration, resource survey and other fields. BRIEF DESCRIPTION OF THE DRAWINGS

[0033] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the conventional technology, the drawings required for use in the embodiments or the conventional technology descriptions are briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0034] Figure 1 is a flow chart of the method of the present invention;

[0035] Figure 2 It is the symmetrical electrode point diagram of the present invention;

[0036] Figure 3 It is the artificial electric field distribution diagram of the present invention;

[0037] Figure 4 The potential map is collected for the network parallel electrodes of the present invention. DETAILED DESCRIPTION

[0038] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the specific embodiments of the present invention are described in detail below in conjunction with the accompanying drawings. In the following description, many specific details are set forth to facilitate a full understanding of the present invention. However, the present invention can be implemented in many other ways different from those described herein, and those skilled in the art can make similar improvements without violating the connotation of the present invention, so the present invention is not limited by the specific embodiments disclosed below.

[0039] like Figure 1 As shown, a method for detecting and analyzing the floor of a coal seam working face comprises the following steps:

[0040] T1: Resistivity measurement method:

[0041] T11: Figure 3 As shown, the two ends of the DC power supply are supplied to the earth through two electrodes A and B buried underground, and a stable electric field is established in the conductive half space below the ground; Figure 2As shown in the figure, a symmetrical quadrupole device is arranged on the surface of the earth, wherein the symmetrical quadrupole device includes two power supply electrodes A and B and two measuring electrodes M and N. When current is sent to the ground through the power supply electrodes A and B, a stable electric field is established in the uniform half space of the underground resistivity ρ. The potential difference is observed at M and N, and the uniform earth resistivity calculation expression is:

[0042]

[0043] in is the device coefficient, and its unit is m. The size of the device coefficient K is only related to the relative positions of the power supply electrodes A, B and the measuring electrodes M, N. When the electrode positions are fixed, the K value can be determined;

[0044] It should be noted that in a uniform isotropic medium, regardless of the pole arrangement, the resistivity calculated based on the measurement results is always equal to the true resistivity of the medium; in actual work, the geoelectric sections usually encountered are generally uneven and complex. When the bottom plate electrical method is still performed using a quadrupole device, the uneven section is replaced by an equivalent uniform section, so formula (1) is still used to calculate the resistivity of the underground medium. The resistivity obtained in this way is not equal to the true resistivity of a certain rock layer, but the result of the combined influence of the resistivity of various rocks within the distribution range of the electric field, which is called the apparent resistivity and is represented by the symbol ρ s Therefore, the expression of apparent resistivity is

[0045] Each electrode is equipped with an A / D converter and performs parallel sampling. It keeps in real-time contact with the host according to the network protocol. When receiving the power supply status command, the electrode sampling part is disconnected to put the electrode in the AB power supply state. Otherwise, it always works in the voltage sampling state and sends the measurement data back to the host in real time through the communication line. In this way, intelligent electrode sampling can be realized, and the electrical data (conductivity measurement value) of each three-dimensional grid is obtained and recorded as R (xi, yi, zi), where (xi, yi, zi) represents the specific position of the grid numbered i in the three-dimensional coordinates; i = 1, 2, 3...I, I is a positive integer, I represents the total number of grids, and i represents any grid number among them;

[0046] The voltage and current data of the natural field, primary field, and secondary field are automatically sampled through the timing relationship between power supply and measurement. The sampling process minimizes the occurrence of idle electrodes. The intelligent electrode is combined with the network system to realize the parallel bottom plate electrical method, thereby greatly reducing the acquisition cost of electrical method data.

[0047] like Figure 4As shown in the figure, according to the different electrode observation devices, the parallel electrical method data acquisition methods are divided into two types: AM method (a) and ABM method (b). The data collected by the parallel electrical method instrument provide technical support for high-density resistivity method, high-resolution georesistivity method, two-dimensional and three-dimensional resistivity imaging interpretation.

[0048] T12: Dual-lane 3D resistivity imaging:

[0049] The general form of the three-dimensional inverse resistivity problem can be expressed as: Δd = GΔm, where G is the Jacobi matrix; Δd is the observed data, d and the forward theoretical value, d 0 The residual vector; Δm is the modification vector of the initial model m; for three-dimensional problems, the model is divided into three-dimensional grids, and the inversion requirement parameter is the conductivity value in each grid unit. The observation data of the three-dimensional inversion is the measured monopole-monopole potential value or monopole-dipole potential difference value. Due to their large range of variation, logarithms are generally used to calibrate the inversion data and model parameters, which is beneficial to improve the stability of the inversion. Due to too many inversion parameters, traditional damped least squares inversion often leads to overly complex models, that is, the so-called redundant structure, which is the structural information that is not required by the data itself or cannot be distinguished, which brings difficulties to the interpretation. Smooth constraints are added to the least squares criterion, and a smooth model is obtained by inversion to improve the stability of the solution. The algorithm for solving the model modification amount is: (G T G+λC T C) Δm=G T Δd, where C is the model smoothing matrix. The electrical data (conductivity measurement value) of each three-dimensional grid is obtained by solving the Jacobi matrix G and calculating the inverse of the large matrix, which is recorded as R(xi, yi, zi), where (xi, yi, zi) represents the specific position of the grid numbered i in the three-dimensional coordinates; i = 1, 2, 3...I, I is a positive integer, I represents the total number of grids, and i represents any grid number;

[0050] Intelligent collection of electrical data by parallel electrodes ensures the synchronization of potential measurement and avoids interference problems of data measured at different times. Such data bodies are particularly suitable for full-space three-dimensional resistivity inversion technology, which can obtain the electrical distribution between working faces. The changes in potential at different locations and elevations are observed by parallel electrical instruments; based on three-dimensional electrical inversion, the resistivity distribution at different depths in the working face and its floor can be obtained. These resistivity distributions provide technical support for objective and accurate geological interpretations, help to better understand the hydrogeological conditions of the rock formations at the bottom of the working face, and thus provide an important reference for safe mining.

[0051] T2: The specific process of environmental interference assessment and correction is as follows:

[0052] The devices in the environment that are turned on are recorded as related devices, and the temperature of each related device is recorded as ST, and the electromagnetic parameters, where the electromagnetic parameters include electromagnetic intensity recorded as P1 and electromagnetic direction; it should be noted that the devices in the environment that are turned on will generate electromagnetic interference or thermal effect interference, etc., which will cause the measured conductivity value to deviate from the true value; the electromagnetic direction applied by each grid for measuring the conductivity value is obtained, and it is recorded as the reference direction;

[0053] Take any one of the grids, obtain the distance between each relevant device and the grid, and record it as the relevant distance A1; obtain the ambient temperature, and record it as HT; obtain the underground medium type between each relevant device and the grid, and set different medium types to correspond to different propagation coefficients. It should be noted that different types of underground media have different propagation effects on electromagnetic waves. For example, the propagation effect of rock-type underground media on electromagnetic waves is better than that of soil-type underground media. This is because electromagnetic waves can reduce energy loss when propagating in rocks, so that the propagation distance is farther and the penetration is deeper, which has a greater impact on the conductivity measurement value of the grid; compare the underground medium type between each relevant device and the grid with all the set underground medium types to match the corresponding propagation coefficient, which is recorded as P2, and use the set formula group The temperature transfer index and electromagnetic transfer index of each relevant device are calculated, where a1, a2, a3, a4, and a5 are respectively the set proportional coefficients; it can be seen from the formula that when the ambient temperature is higher, the relevant device transfers less heat to the environment, and the thermal effect on the grid is more obvious, then the temperature transfer index is higher; when the relevant distance is larger, the heat generated by the operation of the relevant device has a weaker thermal effect on the grid, then the temperature transfer index is smaller; when the transmission coefficient is larger, the electromagnetic intensity is larger, then the electromagnetic transfer index is larger; when the relevant distance is larger, the electromagnetic transfer index is smaller;

[0054] The temperature transfer index of each relevant device is compared and analyzed with the set temperature transfer interval. When the temperature transfer index is greater than the maximum value in the set temperature transfer interval, it means that the thermal effect of the relevant device has a greater impact on the conductivity measurement value of the grid. The temperature transfer index is recorded as the first-level transfer index, and a high thermal effect interference is accumulated once; when the temperature transfer index is within the set temperature transfer interval, the temperature transfer index is recorded as the second-level transfer index, and a medium thermal effect interference is accumulated once; when the temperature transfer index is less than the minimum value in the set temperature transfer interval, the temperature transfer index is recorded as the third-level transfer index, and a low thermal effect interference is accumulated once; the cumulative number of high thermal effect interference, medium thermal effect interference and low thermal effect interference are counted respectively, and they are recorded as U1, U2 and U3 respectively; the first-level transfer index, the second-level transfer index and the third-level transfer index are summed up to obtain the first-level transfer total value, the second-level transfer total value and the third-level transfer total value, and they are recorded as U4, U5 and U6 respectively; using the set formula The thermal effect interference index Ub of the grid is calculated, where b1, b2, and b3 are respectively set proportional coefficients, and e is a natural constant; thus, the thermal effect interference index of each grid is recorded as Ub(xi, yi, zi); it should be noted that the greater the thermal effect, the greater the conductivity measurement value;

[0055] The reference direction of the grid is taken as the positive direction of the horizontal axis, and the reverse direction of the reference direction is taken as the negative direction of the horizontal axis; the electromagnetic direction corresponding to the electromagnetic transfer index of each relevant device at the grid is decomposed on the horizontal axis to obtain the decomposed electromagnetic direction. It should be noted that the decomposed electromagnetic direction is a positive direction or a negative direction; thus, the electromagnetic transfer index of each relevant device at the grid and its corresponding decomposed electromagnetic direction can be obtained, and the electromagnetic transfer index of each relevant device at the grid is recorded as APj, where j = 1, 2, 3...J, J is a positive integer, J represents the total number of relevant devices at the grid, and j represents the serial number of any relevant device at the grid; using the set formula The electromagnetic effect interference index Cb of the grid is calculated, where n is an odd number or an even number. When the decomposed electromagnetic direction is a positive direction, n is an even number; when the decomposed electromagnetic direction is a negative direction, n is an odd number. It should be noted that the direction of the electromagnetic field is also very important for the conductivity value. When the electromagnetic field of the relevant equipment is consistent with the direction of application of the measured conductivity value, the conductivity measurement value increases; conversely, when the electromagnetic field of the relevant equipment is opposite to the direction of application of the measured conductivity value, the conductivity measurement value decreases. Therefore, the electromagnetic effect interference index of each grid is recorded as Cb(xi, yi, zi);

[0056] The thermal effect interference index Ub(xi,yi,zi) and the electromagnetic effect interference index Cb(xi,yi,zi) are calculated by the formula Calculate and obtain the correction index CU (xi, yi, zi), where b4 and b5 are respectively the set proportional coefficients, and γ is the set correction conversion coefficient; multiply the conductivity measurement value obtained by measuring each grid by the corresponding correction index to obtain the corrected conductivity value SR (xi, yi, zi) of the grid point;

[0057] By evaluating and correcting the thermal and electromagnetic effects in the environment, the deviation of the conductivity measurement value can be accurately corrected, the accuracy and reliability of the conductivity measurement data are improved, and a more reliable and accurate basis is provided for subsequent data analysis and application.

[0058] T3: Constructing a three-dimensional conductivity model:

[0059] An underground three-dimensional model is established, and the corrected conductivity value SR(xi, yi, zi) of each grid point is input into the three-dimensional model and annotated to construct a three-dimensional conductivity model, so as to visualize the distribution information of the electrical properties of the underground medium, thereby providing scientific basis and data support for geological exploration, resource survey and other fields.

[0060] The technical features of the above-described embodiments may be arbitrarily combined. To make the description concise, not all possible combinations of the technical features in the above-described embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0061] The above-mentioned embodiments only express several implementation methods of the present invention, and the descriptions thereof are relatively specific and detailed, but they cannot be understood as limiting the scope of the invention patent. It should be pointed out that, for ordinary technicians in this field, several variations and improvements can be made without departing from the concept of the present invention, and these all belong to the protection scope of the present invention. Therefore, the protection scope of the patent of the present invention shall be subject to the attached claims.

Claims

1. A method for detecting and analyzing the floor of a coal seam working face, characterized in that: The following steps are involved: T1: The parallel electrode intelligent acquisition method and double-lane three-dimensional resistivity imaging are used to detect the resistivity value of each three-dimensional grid to obtain the conductivity measurement value of each three-dimensional grid, and the conductivity measurement value is recorded as R(xi, yi, zi), where (xi, yi, zi) represents the specific position of the grid numbered i in the three-dimensional coordinates; i = 1, 2, 3...I, I is a positive integer, I represents the total number of grids, and i represents any grid number among them; T2: Environmental interference assessment and correction, the specific process is as follows: T21: The devices in the environment that are turned on are recorded as related devices. Randomly select one of the grids and quantitatively analyze the thermal effect interference and electromagnetic interference generated by the related devices at the grid point to obtain the temperature transfer index and electromagnetic transfer index of each related device in the grid; T22: Comprehensively analyze the temperature transfer index of each related device at the grid to obtain the thermal effect interference index Gb corresponding to the grid, thereby obtaining the thermal effect interference index of each grid and recording it as Gb(xi, yi, zi); T23: Comprehensively analyze the electromagnetic transfer index of each relevant device at the grid to obtain the electromagnetic effect interference index Cb corresponding to the grid, thereby obtaining the electromagnetic effect interference index of each grid and recording it as Cb(xi, yi, zi); T24: The thermal effect interference index Gb(xi, yi, zi) and the electromagnetic effect interference index are recorded as Cb(xi, yi, zi) through the set formula Calculate and obtain the correction index CG (xi, yi, zi), where b4 and b5 are respectively the set proportional coefficients, and γ is the set correction conversion coefficient; multiply the conductivity measurement value obtained by measuring each grid by the corresponding correction index to obtain the corrected conductivity value SR (xi, yi, zi) of the grid point; T3: Establish an underground three-dimensional model, input the corrected conductivity value SR (xi, yi, zi) of each grid point into the three-dimensional model and mark it to construct a three-dimensional conductivity model.

2. A method for detecting and analyzing the floor of a coal seam working face according to claim 1, characterized in that: The specific process of quantitatively analyzing the thermal effect interference and electromagnetic interference generated by the relevant equipment at this grid point is as follows: Record the devices in the environment that are turned on as relevant devices, obtain the temperature and electromagnetic parameters of each relevant device, where the electromagnetic parameters include electromagnetic intensity and electromagnetic direction; obtain the electromagnetic direction applied by each grid for measuring the conductivity value, and record it as the reference direction; Take any one of the grids, obtain the distance between each relevant device and the grid, and record it as the relevant distance; obtain the ambient temperature and the type of underground medium between each relevant device and the grid; Different medium types are set to correspond to different propagation coefficients, and the underground medium types between the relevant equipment and the grid are compared with all the set underground medium types to match the corresponding propagation coefficients; The equipment temperature ST, ambient temperature HT, related distance A1, electromagnetic intensity P1, and propagation coefficient P2 are combined through the set formula The temperature transfer index and electromagnetic transfer index of each relevant equipment are calculated, where a1, a2, a3, a4, and a5 are set proportional coefficients, and e is a natural constant.

3. A coal seam working face floor detection and analysis method according to claim 2, characterized in that: The specific process of comprehensively analyzing the temperature transfer index of each relevant equipment at the grid is: Compare and analyze the temperature transfer index of each relevant device with the set temperature transfer range to divide the temperature transfer index into a primary transfer index, a secondary transfer index and a tertiary transfer index. The number of the primary transfer index, the secondary transfer index and the tertiary transfer index corresponds to the cumulative number of high thermal effect interference, medium thermal effect interference and low thermal effect interference respectively; sum up the primary transfer index, the secondary transfer index and the tertiary transfer index to obtain the primary transfer total value, the secondary transfer total value and the tertiary transfer total value respectively; The thermal effect interference index of the grid is obtained by formulating and analyzing the cumulative number of high thermal effect interferences, the cumulative number of medium thermal effect interferences, the cumulative number of low thermal effect interferences, the total value of the first-level transfer, the total value of the second-level transfer and the total value of the third-level transfer.

4. A method for detecting and analyzing the floor of a coal seam working face according to claim 3, characterized in that: The specific process of comprehensive analysis of the electromagnetic transfer index of each relevant device at the grid is as follows: The reference direction of the grid is taken as the positive direction of the horizontal axis, and the reverse direction of the reference direction is taken as the negative direction of the horizontal axis; the electromagnetic direction corresponding to the electromagnetic transfer index of each relevant device at the grid is decomposed on the horizontal axis to obtain the decomposed electromagnetic direction; thereby, the electromagnetic transfer index of each relevant device at the grid and its corresponding decomposed electromagnetic direction can be obtained, and the electromagnetic effect interference index of the grid is obtained by cumulative calculation and analysis.

5. A method for detecting and analyzing the floor of a coal seam working face according to claim 1, characterized in that: The implementation steps of the parallel electrode intelligent acquisition method are as follows: A symmetrical quadrupole device is arranged on the surface of the earth to monitor the potential difference of each three-dimensional grid, wherein the symmetrical quadrupole device includes two power supply electrodes, namely electrode A and electrode B, and two measuring electrodes, namely electrode M and electrode N. The observed potential difference at electrode M and electrode N is calculated using the uniform earth resistivity calculation expression to obtain the apparent resistivity; Each electrode is equipped with an A / D converter and performs parallel sampling. It maintains real-time contact with the host according to the network protocol. When receiving the power supply status command, the electrode sampling part is disconnected to keep the electrode in the AB power supply state. Otherwise, it always works in the voltage sampling state and sends the measurement data back to the host in real time through the communication line; and the conductivity measurement value of each three-dimensional grid is recorded as R (xi, yi, zi).

6. A method for detecting and analyzing the floor of a coal seam working face according to claim 1, characterized in that: The implementation steps of double-lane 3D resistivity imaging are as follows: The expression of the three-dimensional inverse resistivity problem is: Δd = GΔm, where G is the Jacobi matrix; Δd is the residual vector between the observed data d and the forward theoretical value d0; Δm is the modification vector of the initial model m; Add smooth constraints to the least squares criterion, invert and obtain a smooth model to improve the stability of the solution. The algorithm for solving the model modification is: (G T G+λC T C) Δm=G T Δd, where C is the model smoothing matrix, and the conductivity measurement value of each three-dimensional grid is obtained by solving the Jacobi matrix G and calculating the large matrix inverse, which is recorded as R(xi, yi, zi).

Images