Coal mine parallel electrical method data acquisition and processing system

Abstract

The invention discloses a coal mine parallel electrical method data acquisition and processing system, and relates to the technical field of coal mine geological exploration. The system comprises n electrodes and m acquisition terminals, and further comprises a processing unit and a prediction model, the n / 2 electrodes are horizontally and uniformly arranged in a survey range, the m acquisition terminals are horizontally and uniformly arranged in the survey range, the processing unit supplies power to the arranged electrodes in sequence, and the processing unit supplies power to the m acquisition terminals. When the electrodes supply power every time, the m acquisition terminals acquire potential information of corresponding positions at the same time and transmit the potential information to the processing unit, and the processing unit calculates the resistivity according to the potential difference between the acquisition terminals. The detail representation of the generation model can be effectively complemented, the prediction model is subjected to feedback correction by introducing a deviation coefficient verification link, and the error probability of the total space layered model can be effectively reduced.

Description

Technical Field

[0001] The invention relates to the technical field of coal mine geological detection, in particular to a coal mine parallel electrical method data acquisition and processing system. Background Art

[0002] Electrical prospecting is mainly divided into two categories, one is the conductive electrical method, and the other is the electromagnetic induction method. The electromagnetic induction method is divided into the electromagnetic profiling method and the electromagnetic sounding method. The conductive electrical method can be divided into the resistivity method, the charging method, the natural electric field method and the induced polarization method. The direct current method is a conductive resistivity method of electrical prospecting. Its basic principle is to use artificial or natural field sources to study the distribution characteristics and laws of the direct current electric field related to the geological body, so as to find minerals and solve certain geological problems. The electrical differences used include the resistivity difference and polarizability difference of rocks and ores, and the measured parameters include apparent resistivity and apparent polarizability. By analyzing the resistivity data obtained by the measurement, the following information can be obtained: 1. Geological structure: determine the stratification of the strata, faults, folds and other geological structures; 2. Ore body distribution: detect the existence, location and scale of the ore body; 3. Aquifer: determine the distribution of groundwater and the characteristics of the aquifer. Comprehensive analysis of multiple factors can fully understand the underground geological situation.

[0003] The parallel electrical method is a type of DC resistivity method. The existing implementation of the parallel electrical method can realize the one-time model formation, but parallel sampling has the following limitations: 1. The underground exploration probe is not intuitive enough for some complex underground component models with large resistivity changes; 2. Parallel sampling lacks a more reliable verification link to correct the deviation of the sampling model; therefore, how to design an acquisition and processing system for the above problems is a technical problem that technical personnel in this field need to solve. Summary of the invention

[0004] In view of the deficiencies in the prior art, the present invention provides a coal mine parallel electrical method data acquisition and processing system, which solves the problems raised in the above-mentioned background technology.

[0005] To achieve the above objectives, the present invention is implemented through the following technical solutions: A parallel electrical data acquisition and processing system for coal mines includes n electrodes and m acquisition terminals, and also includes a processing unit and a prediction model;

[0006] Will The electrodes are evenly arranged horizontally within the survey range, and m acquisition terminals are evenly arranged horizontally within the survey range;

[0007] The processing unit sequentially supplies power to the laid electrodes. Each time the electrodes are powered, the m acquisition terminals simultaneously collect the potential information of the corresponding positions and transmit it to the processing unit. The processing unit calculates the resistivity ρ according to the potential difference between each acquisition terminal. The processing unit transmits the calculated resistivity ρ to the prediction model. The prediction model generates a horizontal layered distribution map i with different depth sections with a spacing of q according to the resistivity ρ calculated by the n electrodes, and generates a smooth curve k between each horizontal layered distribution map i according to the Bezier curve formula. i The prediction model determines the electrode placement coordinate set a for the next stage according to the resistivity ρ dividing line of the horizontal layered distribution map i of each layer. The specific expression of the set a is:

[0008]

[0009] The (x1, y1) is the position of the electrode placement coordinates in the collection a;

[0010] The processing unit sequentially gives the coordinate collection a Each of the electrodes is powered, and each time the electrodes are powered, the m acquisition terminals simultaneously acquire the potential information of the corresponding positions and transmit it to the processing unit. The processing unit calculates the resistivity ρ according to the potential difference between each acquisition terminal, and the processing unit transmits the calculated resistivity ρ to the prediction model. The prediction model is based on The resistivity ρ calculated for each electrode generates a spacing of The horizontal layered distribution map u of different depth sections is inserted into the horizontal layered distribution map i, and d points in the horizontal layered distribution map u are selected and the corresponding smooth curve k in the range is i Point yk i The deviation coefficient λ1 is analyzed and calculated. The calculation formula of the deviation coefficient λ1 is: where y d is the value of the ordinate of the point selected in the horizontal layered distribution diagram u, yk i is the value of the ordinate of the point corresponding to the deviation coefficient λ1, and the deviation coefficient λ1 is brought into the horizontal layered distribution diagram i to generate a smooth curve k u The prediction model determines the electrode placement coordinate set b for the next stage according to the horizontal layered distribution map u of each layer and the resistivity ρ dividing line. The specific expression of the set b is:

[0011]

[0012] The (x2, y2) is the position of the electrode placement coordinates in the collection b;

[0013] The processing unit sequentially gives the coordinate set b Each of the electrodes is powered, and each time the electrodes are powered, the m acquisition terminals simultaneously acquire the potential information of the corresponding positions and transmit it to the processing unit. The processing unit calculates the resistivity ρ according to the potential difference between each acquisition terminal, and the processing unit transmits the calculated resistivity ρ to the prediction model. The prediction model is based on The resistivity ρ calculated for each electrode generates a spacing of The horizontal layered distribution map o of different depth sections is inserted into the horizontal layered distribution map u, and g points in the horizontal layered distribution map o are selected and the corresponding smooth curve k in the range is u Point yk i The deviation coefficient λ2 is analyzed and calculated. The calculation formula of the deviation coefficient λ2 is: where y g is the value of the ordinate of the point selected in the horizontal layered distribution diagram o, yk u is the value of the ordinate of the point corresponding to the deviation coefficient λ2, and the deviation coefficient λ2 is brought into the Bezier curve formula between each horizontal layered distribution graph o to generate a smooth curve k o The prediction model combines the horizontal layered distribution map i, the horizontal layered distribution map u, and the horizontal layered distribution map o respectively and uses the smooth curve k o The gaps between the layered distribution maps are fitted to generate a full-space layered model.

[0014] Furthermore, in parallel electrical prospecting, the calculation formula of the resistivity ρ is as follows:

[0015] ΔUC is the potential difference measured between the two acquisition terminals, Z is the device coefficient, and L is the power supply current of the electrodes arranged in the calculated horizontal layered distribution diagram i; the device coefficient Z is related to the position between the arranged electrodes. When the position of the electrodes arranged on site is determined, the specific value of the device coefficient Z can be determined.

[0016] Furthermore, the spacing q is determined by the area of ​​the field survey range and the density of the electrode layout spacing. The resistivity ρ dividing line is the dividing line of resistivity ρ in different ranges, which represents the relationship between rock conductivity and mineral composition within the survey range and is determined by specific analysis of the field survey.

[0017] Furthermore, the (x1, y1) is taken as a unit node of one meter. The electrodes are uniformly arranged horizontally in sequence within the range of the placement coordinate set a, and the electrodes are preferentially placed at coordinates with smaller values ​​of (x1, y1) until After all electrodes are placed, the (x2, y2) is taken as the value of the node with half a meter as the unit. The electrodes are uniformly arranged horizontally in sequence within the range of the placement coordinate set b, and the electrodes are preferentially placed at coordinates with smaller values ​​of (x2, y2) until All electrodes are placed.

[0018] Furthermore, the values ​​of the corresponding ordinates under the same abscissa between the three adjacent horizontal layered distribution graphs i are used as three control points, and the smooth curve k is obtained by using the Bezier curve formula. i The values ​​of the corresponding ordinates under the same abscissa between the horizontal layered distribution graph u and the upper and lower adjacent horizontal layered distribution graphs i are used as three control points, and the smooth curve k is obtained using the Bezier curve formula ui coefficient, the newly obtained smooth curve k ui Coefficient and smooth curve k i The difference is used as the deviation coefficient λ1 of the prediction model. The deviation coefficient λ1 is multiplied by the values ​​of the corresponding ordinates under the same abscissa in the three horizontal layered distribution graphs i closest to each other as the three new control points. The smooth curve k is calculated using the Bezier curve formula. u The values ​​of the corresponding ordinates under the same abscissa between the horizontal layered distribution graph o and the upper and lower adjacent horizontal layered distribution graphs u are used as three control points, and the smooth curve k is obtained using the Bezier curve formula ou coefficient, the newly obtained smooth curve k ou Coefficient and smooth curve k u The difference is used as the deviation coefficient λ2 of the prediction model. The deviation coefficient λ2 is multiplied by the values ​​of the corresponding ordinates under the same abscissa in the three horizontal layered distribution graphs u closest to each other as the three new control points. The smooth curve k is calculated using the Bezier curve formula. o ; The Bezier curve formula is a mathematical method to define a curve through a series of control points. By adjusting the positions of the control points, the shape of the curve can be accurately controlled and used to correct data deviations.

[0019] Furthermore, the processing unit uses a constant current power supply square wave to sequentially power the arranged electrodes, and calculates and generates a horizontal layered distribution diagram u. The power supply current of the processing unit is In the horizontal layered distribution diagram o, the power supply current of the processing unit is calculated as follows:

[0020] The present invention has the following beneficial effects:

[0021] 1. The parallel electrical data acquisition and processing system for coal mines can refine and enlarge the area of ​​complex underground components by narrowing the detection range in batches, which can effectively supplement the detailed performance of the generated model and make the model features more intuitive and obvious.

[0022] 2. The parallel electrical data acquisition and processing system of the coal mine can effectively reduce the error probability of the full-space layered model and improve the accuracy and authenticity of the generated model by introducing the deviation coefficient verification link and performing feedback correction on the prediction model.

[0023] Of course, any product implementing the present invention does not necessarily need to achieve all of the advantages described above at the same time. BRIEF DESCRIPTION OF THE DRAWINGS

[0024] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings required for describing the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other accompanying drawings can be obtained based on these accompanying drawings without paying creative work.

[0025] Figure 1 This is a block diagram of the coal mine parallel electrical method data acquisition and processing system of the present invention;

[0026] Figure 2 It is a schematic diagram of the process of the present invention;

[0027] Figure 3 It is a schematic diagram of stacking of the full-space layered model of the present invention. DETAILED DESCRIPTION

[0028] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0029] See also Figure 1-3 The present invention provides a technical solution: a parallel electrical method data acquisition and processing system for coal mines, such as Figure 1 As shown, it includes n electrodes and m acquisition terminals, and also includes a processing unit and a prediction model;

[0030] Will The electrodes are evenly arranged horizontally within the survey range, and the m acquisition terminals are evenly arranged horizontally within the survey range;

[0031] The processing unit supplies power to the laid electrodes in sequence. Each time the electrodes are powered, the m acquisition terminals collect the potential information of the corresponding positions at the same time and transmit it to the processing unit. The processing unit calculates the resistivity according to the potential difference between each acquisition terminal. The processing unit transmits the calculated resistivity to the prediction model. The prediction model generates a horizontal layered distribution map i with a spacing of q different depth sections according to the resistivity calculated by the n electrodes, and generates a smooth curve k between each horizontal layered distribution map i according to the Bezier curve formula. i The prediction model determines the electrode placement coordinate set a in the next stage according to the resistivity segmentation line of the horizontal layered distribution map i of each layer. The specific expression of set a is:

[0032]

[0033] (x1, y1) is the position of the electrode placement coordinates in collection a;

[0034] The processing unit gives the coordinate set a in turn Each time the electrode is powered, the m acquisition terminals collect the potential information of the corresponding position at the same time and transmit it to the processing unit. The processing unit calculates the resistivity according to the potential difference between each acquisition terminal. The processing unit transmits the calculated resistivity to the prediction model. The prediction model calculates the resistivity according to the potential difference between each acquisition terminal. The resistivity generation spacing calculated for each electrode is The horizontal layered distribution map u of different depth sections is interspersed into the horizontal layered distribution map i. The prediction model compares the horizontal layered distribution map u with the corresponding smooth curve k within the range. i Analyze and calculate the deviation coefficient λ1, and bring the deviation coefficient λ1 into the horizontal layered distribution diagram i to generate a smooth curve k u The prediction model determines the electrode placement coordinate set b for the next stage according to the horizontal layered distribution map u of each layer and the resistivity segmentation line. The specific expression of set b is:

[0035]

[0036] (x2, y2) is the position of the electrode placement coordinates in set b;

[0037] The processing unit is given to the coordinate set b in turn Each time the electrode is powered, the m acquisition terminals collect the potential information of the corresponding position at the same time and transmit it to the processing unit. The processing unit calculates the resistivity according to the potential difference ΔUC between each acquisition terminal. The processing unit transmits the calculated resistivity to the prediction model. The prediction model is based on The resistivity generation spacing calculated for each electrode is The horizontal layered distribution map o of different depth sections is interspersed into the horizontal layered distribution map u. The prediction model compares the horizontal layered distribution map o with the corresponding smooth curve k within the range. u The deviation coefficient λ2 is analyzed and calculated, and the deviation coefficient λ2 is brought into the Bezier curve formula between each horizontal layered distribution graph o to generate a smooth curve k o ,like Figure 2-3 As shown, the prediction model combines the horizontal layered distribution map i, the horizontal layered distribution map u, and the horizontal layered distribution map o respectively and uses the smooth curve k o The gaps between the layered distribution maps are fitted to generate a full-space layered model; by gradually narrowing the detection range, the area with complex underground components can be refined and enlarged, which can effectively supplement the detailed performance of the generated model and make the model features more intuitive and obvious.

[0038] Among them, in parallel electrical exploration, the calculation formula of resistivity is as follows:

[0039] ρ is the resistivity, ΔUC is the potential difference measured between the two acquisition terminals, Z is the device coefficient, and L is the power supply current of the electrodes arranged in the calculated horizontal layered distribution diagram i; the device coefficient Z is related to the position between the arranged electrodes. When the position of the electrodes arranged on site is determined, the specific value of the device coefficient Z can be determined.

[0040] Among them, the spacing q is determined by the area of ​​the on-site survey range and the density of the electrode layout spacing. The resistivity dividing line is the dividing line of different resistivity ranges, which represents the relationship between rock conductivity and mineral composition within the survey range and is determined by specific analysis of the on-site survey.

[0041] Among them, (x1, y1) is taken as a unit node of one meter. The electrodes are evenly arranged horizontally in the range of the placement coordinate set a, and the electrodes are placed preferentially at the coordinates with smaller values ​​(x1, y1) until After all electrodes are placed, (x2, y2) is taken as the value of the node with half a meter as the unit. The electrodes are evenly arranged horizontally in the range of the placement coordinate set b, and the electrodes are placed preferentially at the coordinates with smaller values ​​(x2, y2) until All electrodes are placed.

[0042] Among them, the values ​​of the corresponding ordinates under the same abscissa between the three adjacent horizontal layered distribution graphs i are used as three control points, and the smooth curve k is obtained using the Bezier curve formula i ; The values ​​of the corresponding ordinates under the same abscissa between the horizontal layered distribution graph u and the upper and lower adjacent horizontal layered distribution graphs i are used as three control points, and the smooth curve k is obtained using the Bezier curve formula uiCoefficient, the new smooth curve k ui Coefficient and smooth curve k i The difference is used as the deviation coefficient λ1 of the prediction model. The deviation coefficient λ1 is multiplied by the values ​​of the corresponding ordinates under the same abscissa in the three horizontal layered distribution graphs i closest to each other as the new three control points. The smooth curve k is calculated using the Bezier curve formula. u ; The values ​​of the corresponding ordinates under the same abscissa between the horizontal layered distribution graph o and the upper and lower adjacent horizontal layered distribution graphs u are used as three control points, and the smooth curve k is obtained using the Bezier curve formula ou Coefficient, the new smooth curve k ou Coefficient and smooth curve k u The difference is used as the deviation coefficient λ2 of the prediction model. The deviation coefficient λ2 is multiplied by the values ​​of the corresponding ordinates under the same abscissa in the three horizontal layered distribution graphs u closest to each other as the new three control points. The smooth curve k is calculated using the Bezier curve formula. o ; The Bezier curve formula is a mathematical method to define a curve through a series of control points. By adjusting the positions of the control points, the shape of the curve can be precisely controlled and used to correct data deviations. Introducing the deviation coefficient verification link and performing feedback correction on the prediction model can effectively reduce the error probability of the full-space layered model and improve the accuracy and authenticity of the generated model.

[0043] Among them, the processing unit uses a constant current power supply square wave to power the laid electrodes in sequence, and the processing unit power supply current is calculated to generate the horizontal layered distribution diagram u. In the horizontal layered distribution diagram o, the power supply current of the processing unit is calculated as follows: Compared with Collection A and Collection B, the initial survey range is gradually reduced, and the demand for DC power supply will also be reduced accordingly, which is more conducive to the accuracy of the later model.

[0044] The above description is only a preferred specific implementation manner of the present invention, but the protection scope of the present invention is not limited thereto. Any technician familiar with the technical field can make equivalent replacements or changes according to the technical scheme and inventive concept of the present invention within the technical scope disclosed by the present invention, which should be covered by the protection scope of the present invention.

Claims

1. The parallel electrical data acquisition and processing system for coal mines includes n electrodes and m acquisition terminals. The electrodes are evenly arranged horizontally within the survey range, and m acquisition terminals are evenly arranged horizontally within the survey range; the characteristics are: Also included are processing units and prediction models; The processing unit sequentially supplies power to the laid electrodes. Each time the electrodes are powered, the m acquisition terminals simultaneously collect the potential information of the corresponding positions and transmit it to the processing unit. The resistivity is calculated according to the potential difference between each acquisition terminal, and the calculated resistivity is transmitted to the prediction model. The prediction model generates a horizontal layered distribution map i with different depth sections with a spacing of q according to the resistivity calculated by the n electrodes, and generates a smooth curve k between each horizontal layered distribution map i according to the Bezier curve formula. i , according to the resistivity segmentation line of the horizontal layered distribution map i of each layer, the coordinate set a of the electrode placement in the next stage is determined; The processing unit sequentially gives the coordinate set a The electrodes are powered, the acquisition terminals acquire potential information and transmit it to the processing unit, the resistivity is calculated according to the potential difference between each acquisition terminal, the calculated resistivity is transmitted to the prediction model, and the prediction model is based on The resistivity generation spacing calculated for each electrode is The horizontal layered distribution map u of different depth sections is inserted into the horizontal layered distribution map i, and d points in the horizontal layered distribution map u are selected and the corresponding smooth curve k in the range is i Point yk i The deviation coefficient λ1 is analyzed and calculated. The calculation formula of the deviation coefficient λ1 is: where y d is the value of the ordinate of the point selected in the horizontal layered distribution diagram u, yk i is the value of the ordinate of the point corresponding to the deviation coefficient λ1, and the deviation coefficient λ1 is brought into the horizontal layered distribution diagram i to generate a smooth curve k u , according to the horizontal layered distribution map u of each layer, the resistivity segmentation line determines the electrode placement coordinate set b for the next stage, The processing unit sequentially gives the coordinate set b The electrodes are powered, the acquisition terminals acquire potential information and transmit it to the processing unit, the resistivity is calculated according to the potential difference between each acquisition terminal, the calculated resistivity is transmitted to the prediction model, and the prediction model is based on The resistivity generation spacing calculated for each electrode is The horizontal layered distribution map o of different depth sections is inserted into the horizontal layered distribution map u, and g points in the horizontal layered distribution map o are selected and the corresponding smooth curve k in the range is u Point yk i The deviation coefficient λ2 is analyzed and calculated. The calculation formula of the deviation coefficient λ2 is: where y g is the value of the ordinate of the point selected in the horizontal layered distribution diagram o, yk u is the value of the ordinate of the point corresponding to the deviation coefficient λ2, and the deviation coefficient λ2 is brought into the Bezier curve formula between each horizontal layered distribution graph o to generate a smooth curve k o , respectively combine the horizontal layered distribution map i, the horizontal layered distribution map u, and the horizontal layered distribution map o and use the smooth curve k o The gaps between the layered distribution maps are fitted to generate a full-space layered model.

2. The coal mine parallel electrical method data acquisition and processing system according to claim 1 is characterized in that: In parallel electrical prospecting, the resistivity calculation formula is as follows: ρ is the resistivity, ΔUC is the potential difference measured between the two acquisition terminals, Z is the device coefficient, and L is the power supply current of the electrodes arranged in the calculated horizontal layered distribution map i; The specific representation of the collection a is: The (x1, y1) is the position of the electrode placement coordinates in the collection a; The specific representation of the collection b is: The (x2, y2) is the position of the electrode placement coordinates in the collection b.

3. The coal mine parallel electrical method data acquisition and processing system according to claim 1, characterized in that: The spacing q is determined by the area of ​​the field survey range and the density of the electrode layout spacing, and the resistivity dividing line is the dividing line of different resistivity ranges.

4. The coal mine parallel electrical method data acquisition and processing system according to claim 1, characterized in that: The (x1, y1) is taken as a unit node of one meter. The electrodes are uniformly arranged horizontally in sequence within the range of the placement coordinate set a, and the electrodes are preferentially placed at coordinates with smaller values ​​of (x1, y1) until After all electrodes are placed, the (x2, y2) is taken as the value of the node with half a meter as the unit. The electrodes are uniformly arranged horizontally in sequence within the range of the placement coordinate set b, and the electrodes are preferentially placed at coordinates with smaller values ​​of (x2, y2) until All electrodes are placed.

5. The coal mine parallel electrical method data acquisition and processing system according to claim 1, characterized in that: The values ​​of the corresponding ordinates under the same abscissa between the three adjacent horizontal layered distribution graphs i are used as three control points, and the smooth curve k is obtained by using the Bezier curve formula. i The values ​​of the corresponding ordinates under the same abscissa between the horizontal layered distribution graph u and the upper and lower adjacent horizontal layered distribution graphs i are used as three control points, and the smooth curve k is obtained using the Bezier curve formula ui coefficient, the newly obtained smooth curve k ui Coefficient and smooth curve k i The difference is used as the deviation coefficient λ1 of the prediction model. The deviation coefficient λ1 is multiplied by the values ​​of the corresponding ordinates under the same abscissa in the three horizontal layered distribution graphs i closest to each other as the three new control points. The smooth curve k is calculated using the Bezier curve formula. u The values ​​of the corresponding ordinates under the same abscissa between the horizontal layered distribution graph o and the upper and lower adjacent horizontal layered distribution graphs u are used as three control points, and the smooth curve k is obtained using the Bezier curve formula ou coefficient, the newly obtained smooth curve k ou Coefficient and smooth curve k u The difference is used as the deviation coefficient λ2 of the prediction model. The deviation coefficient λ2 is multiplied by the values ​​of the corresponding ordinates under the same abscissa in the three horizontal layered distribution graphs u closest to each other as the three new control points. The smooth curve k is calculated using the Bezier curve formula. o .

6. The coal mine parallel electrical method data acquisition and processing system according to claim 1, characterized in that: The processing unit uses a constant current square wave to supply power to the arranged electrodes in sequence, and calculates and generates a horizontal layered distribution diagram u. The power supply current of the processing unit is In the horizontal layered distribution diagram o, the power supply current of the processing unit is calculated as follows:

Images