A coal seam working face floor detection and analysis method

Abstract

The application discloses a coal seam working face floor detection and analysis method and relates to the technical field of detection and analysis; through parallel electric method observation of potential change conditions of different positions and different elevations, and through three-dimensional electric method inversion, the resistivity distribution conditions of different depths of the working face and the floor thereof can be obtained, the resistivity distribution conditions provide technical support for giving objective and accurate geological interpretation, help better understand the hydrogeological conditions of the working face floor stratum, and thus provide an important reference basis for safe mining; through evaluation and correction of thermal effects and electromagnetic effects existing in the environment, the deviation of the measured conductivity value can be accurately corrected, and the accuracy and reliability of the conductivity measurement data are improved; then the corrected conductivity values SR(xi, yi, zi) of the grid points are input into a three-dimensional model and are labeled to construct a conductivity three-dimensional model, so that the electrical property distribution information of the underground medium can be visually displayed.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of detection analysis, in particular to a coal seam working face floor detection analysis method. BACKGROUND

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

[0003] Traditional geological exploration methods such as drilling and geological profile can provide some geological information, but there are limitations in comprehensively understanding the hydrogeological conditions of the floor. These methods can usually only provide local geological information and are difficult to cover the hydrogeological conditions of the entire working face floor rock stratum. In addition, traditional geological exploration methods require a large amount of drilling, measurement and other operations, which not only consumes time and effort, but also has high cost. SUMMARY

[0004] Therefore, it is necessary to provide a coal seam working face floor detection analysis method in view of the above background technology.

[0005] The purpose of the present application can be achieved by the following technical scheme: a coal seam working face floor detection analysis method, comprising the following steps:

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

[0007] T2: environmental interference evaluation correction, the specific process is:

[0008] T21: recording the devices in the environment in the open state as related devices, selecting one of the grids, and quantitatively analyzing 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 at the grid;

[0009] T22: comprehensively analyzing the temperature transfer index of each related device at the grid to obtain the thermal effect interference index Gb corresponding to the grid, so that the thermal effect interference index of each grid is recorded as Gb(xi, yi, zi);

[0010] T23: The electromagnetic transmission index of each related device at the grid is analyzed comprehensively to obtain the electromagnetic effect interference index Cb of the grid, so that the electromagnetic effect interference index of each grid is recorded as Cb(xi, yi, zi);

[0011] T24: The thermal effect interference index Gb(xi, yi, zi) and the electromagnetic effect interference index Cb(xi, yi, zi) are calculated by a set formula to obtain the correction index CG(xi, yi, zi), wherein b4 and b5 are set proportional coefficients, and γ is a set correction conversion coefficient; the measured conductivity value of each grid is multiplied by the corresponding correction index to obtain the corrected conductivity value SR(xi, yi, zi) of the grid point;

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

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

[0014] The devices in the environment in the open state are recorded as related devices, the temperature and electromagnetic parameters of each related device are obtained, wherein the electromagnetic parameters include electromagnetic intensity and electromagnetic direction; the electromagnetic direction applied by each grid for measuring the conductivity value is obtained and recorded as the reference direction;

[0015] Any one grid is taken, the distance between each related device and the grid is obtained and recorded as the related distance; the environmental temperature, the underground medium type between each related device and the grid is obtained; different medium types are set to correspond to different propagation coefficients, and the underground medium type between each related device and the grid is compared with all set underground medium types to match the corresponding propagation coefficient;

[0016] The device temperature ST, the environmental temperature HT, the related distance A1, the electromagnetic intensity P1, and the propagation coefficient P2 are calculated by a set formula group to obtain the temperature transmission index and the electromagnetic transmission index of each related device, wherein 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 transmission index of each related device at the grid is as follows:

[0018] The temperature transmission index of each related device is compared with the set temperature transmission interval to divide the temperature transmission index into a first transmission index, a second transmission index, and a third transmission index, the number of the first transmission index, the second transmission index, and the third transmission index respectively corresponds to the cumulative number of high thermal effect interference, moderate thermal effect interference, and low thermal effect interference, and is respectively recorded as U1, U2, and U3; the first transmission index, the second transmission index, and the third transmission index are respectively summed to obtain a first transmission total value, a second transmission total value, and a third transmission total value, and are respectively recorded as U4, U5, and U6;

[0019] The cumulative number U1 of high thermal effect interference, the cumulative number U2 of moderate thermal effect interference, the cumulative number U3 of low thermal effect interference, the first transmission total value U4, the second transmission total value U5, and the third transmission total value U6 are substituted into the set formula to obtain the thermal effect interference index Ub of the grid, wherein b1, b2, and b3 are respectively set proportion coefficients, and e is a natural constant.

[0020] In some embodiments, the specific process of comprehensively analyzing the electromagnetic transmission index of each related device at the grid is as follows:

[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 transmission index of each related device at the grid is decomposed on the horizontal axis to obtain a decomposed electromagnetic direction; thus, the electromagnetic transmission index of each related device at the grid and the corresponding decomposed electromagnetic direction are obtained, and the electromagnetic transmission index of each related device at the grid is recorded as APj, wherein j=1, 2, 3…J, J is a positive integer, J represents the total number of related devices at the grid, and j represents the serial number of any related device at the grid; the electromagnetic effect interference index Cb of the grid is calculated by the set formula , wherein n is an odd or even number, n is an even number when the decomposed electromagnetic direction is a positive direction, and n is an odd number when the decomposed electromagnetic direction is a negative direction.

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

[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 electrodes, electrode A and electrode B, and two measurement electrodes, electrode M and electrode N; the apparent resistivity is calculated by using the uniform earth resistivity calculation expression to obtain the observed potential difference at the electrode M and the electrode N;

[0024] Each electrode is equipped with an A / D converter and parallel sampling, and is kept in real-time contact with the host computer according to the network protocol, and the electrode sampling part is disconnected when receiving the power supply state command, so that the electrode is in the AB power supply state, otherwise it always works in the voltage sampling state, and the measured data is sent back to the host computer 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 the dual-lane three-dimensional resistivity imaging are:

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

[0027] The smoothing constraint is added to the least square criterion, and a smoothed 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, wherein C is a model smoothing matrix, and the conductivity measurement value of each three-dimensional grid is recorded as R(xi, yi, zi) by solving the Jacobi matrix G and the calculation of the inverse of the large matrix.

[0028] Compared with the prior art, the beneficial effects of the present application are:

[0029] 1. The intelligent parallel electrode acquisition of electrical data ensures the synchronization of potential measurement and avoids the interference problem of different time measurement data. Such data body is particularly suitable for using full-space three-dimensional resistivity inversion technology to obtain the electrical property distribution between working surfaces;

[0030] 2. The potential change at different positions and elevations is observed by the parallel electrical instrument, and the three-dimensional electrical method inversion can obtain the resistivity distribution at different depths in the working surface and its floor, which provides technical support for objective and accurate geological interpretation, and helps better understand the hydrogeological conditions of the working surface floor strata, thereby providing an important reference for safe mining;

[0031] 3. By evaluating and correcting the thermal effect and electromagnetic effect existing 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 correction conductivity value SR(xi, yi, zi) of each grid point into the three-dimensional model and marking to construct the conductivity three-dimensional model, so as to visualize the distribution information of the electrical properties of the underground medium, thereby providing scientific basis and data support for the fields of geological exploration, resource investigation and the like. BRIEF DESCRIPTION OF DRAWINGS

[0033] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the drawings needed to be used in the embodiments or the prior art description will be briefly introduced. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative labor.

[0034] Figure 1 The method flowchart of the present application;

[0035] Figure 2 The symmetric electrode point map of the present application;

[0036] Figure 3 The artificial electric field distribution map of the present application;

[0037] Figure 4 The network parallel electrode potential map of the present application. DETAILED DESCRIPTION

[0038] In order to make the above-mentioned purposes, features and advantages of the present application more apparent and easy to understand, the specific embodiments of the present application will be described in detail below with reference to the drawings. In the following description, many specific details are set forth in order to provide a thorough understanding of the present application. However, the present application can be implemented in many other ways different from those described herein, and those skilled in the art can make similar improvements without departing from the scope of the present application, therefore the present application is not limited to the specific embodiments disclosed below.

[0039] As shown in Figure 1 A coal seam working face floor detection analysis method, comprising the following steps:

[0040] T1: resistivity measurement method:

[0041] T11: as shown in Figure 3 The two ends of the direct current power supply are supplied to the ground through two electrodes A and B buried underground, and a stable electric field is established in the conductive half-space below the ground; as shown in Figure 2As shown in the figure, a symmetric quadrupole device is arranged on the earth surface, wherein the symmetric quadrupole device comprises two power supply electrodes A, B and two measuring electrodes M, N; when current is sent to the ground through the power supply electrodes A, B, a stable electric field is established in a uniform half-space with a resistivity of p. The potential difference is observed at M and N, and the expression for calculating the resistivity of the uniform earth is

[0042]

[0043] wherein K is a device coefficient, and the unit is m; the size of the device coefficient K is only related to the mutual position of the power supply electrodes A, B and the measuring electrodes M, N; when the electrode position is fixed, the value of K can be determined;

[0044] It should be noted that in a uniform isotropic medium, regardless of the electrode arrangement, the calculated resistivity is always equal to the true resistivity of the medium according to the measurement results; in actual work, the ground electrical section commonly encountered is generally non-uniform and relatively complex; when the bottom plate method is still used, the non-uniform 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 is the result of the comprehensive influence of various rock resistivities within the range of the electric field distribution, which is called apparent resistivity and is represented by the symbol p s , so the expression for apparent resistivity is

[0045] Each electrode is equipped with an A / D converter and performs parallel sampling, and is kept in real-time contact with the host computer according to the network protocol; when receiving a power supply state command, the electrode sampling part is disconnected, so that the electrode is in the AB power supply state, otherwise it always works in the voltage sampling state, and the measurement data is sent back to the host computer in real time through the communication line; in this way, intelligent electrode sampling can be realized, and the electrical property data (conductivity measurement value) of each three-dimensional grid is obtained, denoted as R(xi, yi, zi), wherein (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 the grid number of any one grid;

[0046] Through the time sequence relationship between power supply and measurement, the natural field, primary field, secondary field voltage data and current data are automatically sampled, the sampling process maximally reduces the occurrence of idle electrodes, the intelligent electrodes are combined with the network system, parallel bottom plate method is realized, and the cost of collecting electrical property data is greatly reduced;

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

[0048] T12: Three-dimensional resistivity imaging of double lane:

[0049] The general form of the three-dimensional resistivity inverse problem can be expressed as: Δd=GΔm, where G is the Jacobi matrix; Δd is the residual vector of the observation data, d and the theoretical value of the forward problem, d0; Δm is the modification vector of the initial model m; for a three-dimensional problem, the model is divided into a three-dimensional grid, and the inversion requires parameters, which are the conductivity values in each grid cell. The observation data of three-dimensional inversion are the measured single-pole-single-pole potential values or single-pole-dipole potential difference values. Because of their large range of variation, the inversion data and model parameters are generally calibrated by logarithm, which is conducive to improving the stability of inversion. Because there are too many inversion parameters, the traditional damped least squares inversion often leads to an overly complex model, i.e. the so-called redundant structure, which is a structure information that is not required by the data itself or is not resolvable, bringing difficulties to interpretation. In the least squares criterion, a smoothing constraint is added, and a smooth model is obtained by inversion, improving the stability of the solution. The algorithm for solving the model modification quantity is: (G T G+λC T C)Δm=G T Δd, where C is the model smoothing matrix. The calculation of the Jacobi matrix G and the inverse of the large matrix is used to obtain the electrical data (conductivity measurement value) of each three-dimensional grid, denoted 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 takes a positive integer value, I represents the total number of grids, and i represents the number of any one grid;

[0050] By intelligent acquisition of electrical data by parallel electrodes, the synchronization of potential measurement is ensured, and the interference problem of different time measurement data is avoided. Such data volume is particularly suitable for using full-space three-dimensional resistivity inversion technology, which can obtain the electrical property distribution between working surfaces. By observing the potential change at different positions and elevations by the parallel electrical method instrument, the resistivity distribution at different depths in the working surface and its floor can be obtained by three-dimensional electrical method inversion. These resistivity distribution conditions provide technical support for giving objective and accurate geological interpretation, and help better understand the hydrogeological conditions of the working surface floor strata, thereby providing an important reference for safe mining.

[0051] The specific process of environmental interference evaluation correction is as follows:

[0052] The device in the environment in the open state is recorded as a related device, the temperature of each related device is recorded as ST, and the electromagnetic parameter, wherein the electromagnetic parameter includes electromagnetic intensity P1 and electromagnetic direction; It should be noted that the device in the open state in the environment will generate electromagnetic interference or thermal effect interference, etc., so that the measured conductivity value deviates from the true value; The electromagnetic direction applied to the measured conductivity value is obtained, and is recorded as the reference direction;

[0053] Take one of the grids, obtain the distance between each related device and the grid, and record it as the related distance A1; Obtain the ambient temperature and record it as HT; Obtain the type of underground medium between each related device and the grid, set different medium types to correspond to different propagation coefficients, and it should be noted that different types of underground medium have different effects on electromagnetic wave propagation, for example, the propagation effect of rock type underground medium on electromagnetic wave is better than that of soil type underground medium, because electromagnetic wave can reduce energy loss when propagating in rock, so the propagation distance is farther, the penetration is deeper, and the influence on the conductivity measurement value of the grid is greater; The type of underground medium between each related device and the grid is compared with all the set types of underground medium to match the corresponding propagation coefficient P2, and the set formula group The temperature transfer index and the electromagnetic transfer index of each related device are calculated by using the set formula group, wherein a1, a2, a3, a4 and a5 are the set proportion coefficients; According to the formula, when the ambient temperature is higher, the heat transfer from the related device to the environment is less, the thermal effect on the grid is more obvious, and the temperature transfer index is larger; When the related distance is larger, the heat generated by the operation of the related device has weaker thermal effect on the grid, and the temperature transfer index is smaller; When the propagation coefficient is larger, the electromagnetic intensity is larger, and the electromagnetic transfer index is larger; When the related distance is larger, the electromagnetic transfer index is smaller;

[0054] The temperature transfer index of each related device is compared with the set temperature transfer interval. When the temperature transfer index is greater than the maximum value in the set temperature transfer interval, it is indicated that the thermal effect of the related device has a relatively large influence on the conductivity measurement value of the grid. Therefore, the temperature transfer index is recorded as a 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 a second-level transfer index, and a moderate 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 a third-level transfer index, and a low thermal effect interference is accumulated once. The accumulated number of high, moderate, and low thermal effect interferences is counted and recorded as U1, U2, and U3, respectively. The first, second, and third level transfer indexes are summed to obtain the first, second, and third level transfer total values, which are recorded as U4, U5, and U6, respectively. The thermal effect interference index Ub of the grid is calculated using a set formula where b1, b2, and b3 are set proportionality coefficients, and e is the natural constant. Thus, the thermal effect interference index Ub 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 opposite 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 related 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 either positive or negative. Thus, the electromagnetic transfer index of each related device at the grid and its corresponding decomposed electromagnetic direction are obtained, and the electromagnetic transfer index of each related device at the grid is recorded as APj, where j=1, 2, 3,..., J, J is a positive integer representing the total number of related devices at the grid, and j represents the serial number of any related device at the grid. The electromagnetic effect interference index Cb of the grid is calculated using a set formula where n is an odd or even number. When the decomposed electromagnetic direction is positive, n is an even number. When the decomposed electromagnetic direction is negative, n is an odd number. It should be noted that the direction of the electromagnetic field also has a significant impact on the conductivity value. When the electromagnetic field of the related device is consistent with the applied direction of the conductivity measurement value, the conductivity measurement value increases. Conversely, when the electromagnetic field of the related device is opposite to the applied direction of the conductivity measurement value, the conductivity measurement value decreases. Thus, the electromagnetic effect interference index Cb 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 recorded as a set formula The calculation is performed to obtain a correction index CU(xi, yi, zi), wherein b4 and b5 are respectively a set proportion coefficient, and γ is a set correction conversion coefficient; the measured conductivity value of each grid is multiplied 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 effect and electromagnetic effect existing 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 conductivity three-dimensional model

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

[0060] The technical features of the above-described embodiments can be combined arbitrarily, and in order to make the description concise, all possible combinations of the technical features in the above-described embodiments are not described, however, as long as the combination of the technical features does not exist contradictory, it should be considered as the scope of the present application.

[0061] The above-described embodiments only express several embodiments of the present application, the description is more specific and detailed, but it should not be understood as a limitation on the scope of the patent. It should be pointed out that for ordinary skilled in the art, without departing from the concept of the present application, a number of variations and improvements can be made, which are all within the scope of the present application. Therefore, the scope of the patent protection of the present application should be subject to the appended 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 3D resistivity imaging are used to detect the resistivity value of each 3D grid to obtain the conductivity measurement value of each 3D 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 3D 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: Record the devices in the environment that are turned on as related devices. Randomly select one grid point and quantitatively analyze the thermal effect interference and electromagnetic interference generated by the related devices at that 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. The thermal effect interference index of each grid can be recorded 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. The electromagnetic effect interference index of each grid can be recorded 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 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 the three-dimensional conductivity model.

2. A coal seam working face floor detection and analysis method according to claim 1, characterized in that: The specific process of quantitatively analyzing the thermal effect interference and electromagnetic interference generated by 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 to each grid for measuring the conductivity value and record it as the reference direction; Randomly select one of the grids and obtain the distance between each relevant device and the grid, which is recorded as the relevant distance; obtain the ambient temperature and the type of underground medium between each relevant device and the grid; Different propagation coefficients are set for different medium types, and the underground medium types between each relevant device and grid are compared with all the set underground medium types to match the corresponding propagation coefficients; 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.

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 as follows: The temperature transfer index of each relevant device is compared and analyzed 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. The primary transfer index, the secondary transfer index and the tertiary transfer index are summed up 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 interference, the cumulative number of medium thermal effect interference, the cumulative number of low thermal effect interference, 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 coal seam working face floor detection and analysis method according to claim 3, characterized in that: The specific process of comprehensively analyzing 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 opposite 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. The 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. 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, maintaining real-time communication with the host according to the network protocol. When receiving the power supply status command, the electrode sampling part is disconnected, putting 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 coal seam working face floor detection and analysis method 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 and inversely obtain a smooth model 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, 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