Coal mine parallel electric method data acquisition and processing system

Abstract

The application discloses a coal mine parallel electric method data acquisition and processing system and relates to the technical field of coal mine geological detection.The application comprises n electrodes and m acquisition terminals, further comprises a processing unit and a prediction model, n / 2 electrodes are evenly arranged in the survey range, m acquisition terminals are evenly arranged in the survey range, the processing unit sequentially supplies power to the arranged electrodes, when the electrodes are powered each time, the m acquisition terminals collect the potential information of the 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 each acquisition terminal.The application can refine and enlarge the area with complex underground components by reducing the detection range in batches, can effectively make up the details of the generated model, can feed back and correct the prediction model through the introduction of the deviation coefficient verification link, and can effectively reduce the error probability of the full-space layered model.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of coal mine geological exploration, in particular to a coal mine parallel electric method data acquisition and processing system. BACKGROUND

[0002] Electric method exploration is mainly divided into two categories, one is the conduction type electric method, and the other is the electromagnetic induction method. The electromagnetic induction method is divided into electromagnetic profile method and electromagnetic sounding method, and the conduction type electric method can be divided into resistivity method, charging method, natural electric field method and induced polarization method. The direct current method is a conduction type resistivity method of electric method exploration, and its basic principle is to use artificial field source or natural field source to study the distribution characteristics and laws of direct current electric field related to geological bodies, so as to find ore and solve some geological problems. The electrical property difference used includes the resistivity difference and polarization rate difference of rocks and minerals, and the measured parameters include apparent resistivity and apparent polarization rate. Through analyzing the measured resistivity data, the following information can be obtained: 1. Geological structure: determine the stratification of strata, faults, folds and other geological structures; 2. Ore body distribution: detect the existence, location and size of ore body; 3. Aquifer: determine the distribution of underground water and the characteristics of aquifer. Comprehensive analysis of various factors can comprehensively understand the underground geological conditions.

[0003] The parallel electric method is a kind of direct current resistivity method, and the existing parallel electric method implementation can realize model forming once, but the parallel sampling has the following limitations: 1. The underground exploration is not intuitive and obvious for some complex underground composition models with large resistivity change; 2. The parallel sampling lacks a more reliable verification link to correct the deviation of the sampling model; therefore, how to design a collection and processing system for the above problems is a technical problem to be solved by the person skilled in the art. SUMMARY

[0004] In view of the deficiencies of the prior art, the present application provides a coal mine parallel electric method data acquisition and processing system, which solves the problems raised in the background art.

[0005] To achieve the above purpose, the present application realizes the following technical scheme: a coal mine parallel electric method data acquisition and processing system, comprising n electrodes and m acquisition terminals, further comprising a processing unit and a prediction model.

[0006] The The n electrodes are evenly arranged in the survey range, and the m acquisition terminals are evenly arranged in the survey range;

[0007] The processing unit sequentially supplies power to the arranged electrodes, and the m acquisition terminals collect the potential information of the corresponding positions at the same time when each electrode is powered and transmit it to the processing unit, and the processing unit calculates the resistivity according to the potential difference between each acquisition terminal , the processing unit will transmit the calculated resistivity to the prediction model, the prediction model generates horizontal layered distribution map i of interval q different depth sections according to the n electrode calculated resistivity , and generates smooth curve between each horizontal layered distribution map i according to Bezier curve formula , the prediction model generates horizontal layered distribution map i of interval q different depth sections according to the n electrode calculated resistivity The split line determines the next stage electrode placement coordinate set a, and the specific representation of the set a is:

[0008]

[0009] The is the position of the electrode placement coordinate in the set a;

[0010] The processing unit supplies power to the m electrodes in the coordinate set a in turn, and the m acquisition terminals collect the potential information of the corresponding position at the same time when each electrode is powered, and transmit it to the processing unit, and the processing unit calculates the resistivity according to the potential difference between each acquisition terminal , the processing unit will transmit the calculated resistivity to the prediction model, the prediction model generates horizontal layered distribution map i of interval q different depth sections according to the n electrode calculated resistivity , and generates smooth curve between each horizontal layered distribution map i according to Bezier curve formula , the prediction model generates horizontal layered distribution map u of interval q different depth sections according to the n electrode calculated resistivity , the horizontal layered distribution map u is inserted into the horizontal layered distribution map i, and the points of the d points in the horizontal layered distribution map u and the corresponding smooth curve are analyzed and calculated to obtain the deviation coefficient , the calculation formula of the deviation coefficient is , wherein is the value of the ordinate of the selected point in the horizontal layered distribution map u, is the value of the ordinate of the corresponding point of the deviation coefficient , and the deviation coefficient is brought into the horizontal layered distribution map i to generate smooth curve , the prediction model generates horizontal layered distribution map u of interval q different depth sections according to the n electrode calculated resistivity The split line determines the next stage electrode placement coordinate set b, and the specific representation of the set b is:

[0011]

[0012] The is the position of the electrode placement coordinate in the set b;

[0013] The processing unit sequentially assigns values ​​to the coordinate set b. Each electrode is powered, and each time the electrode is powered, m acquisition terminals simultaneously acquire potential information at the corresponding location and transmit it to the processing unit. The processing unit calculates the resistivity based on the potential difference between each acquisition terminal. The processing unit will calculate the resistivity. Transmitted to the prediction model, the prediction model according to Resistivity calculated for each electrode Generate spacing is A horizontal layered distribution map o at different depths is interspersed within a horizontal layered distribution map u. G points in the horizontal layered distribution map o are selected, and their corresponding smooth curves within that range are compared. point Analysis and calculation of the deviation coefficient Deviation coefficient The calculation formula is: ,in The value of the ordinate of a selected point in the horizontal layered distribution map o. Deviation coefficient The value of the ordinate of the corresponding point, and the deviation coefficient. By substituting the Bézier curve formula between the various horizontal layered distribution maps o, a smooth curve is generated. The prediction model combines horizontal layered distribution map i, horizontal layered distribution map u, and horizontal layered distribution map o, respectively, and uses a smooth curve. A full-space layered model is generated by fitting the gaps between the layered distribution maps.

[0014] Furthermore, in parallel electrical resistivity exploration, the resistivity The specific calculation formula is as follows:

[0015] , The potential difference measured between the two acquisition terminals. For device coefficients, To calculate the supply current of the electrodes arranged in the horizontal layered distribution map i; the device coefficient The device coefficient is related to the position between the electrodes. When the position of the electrodes laid out on site is determined, the device coefficient is... The specific value can then be determined.

[0016] Furthermore, the spacing q is determined by the area of ​​the on-site survey and the density of the electrode spacing, and the resistivity The dividing lines represent different ranges of resistivity. The demarcation line represents the relationship between rock conductivity and mineral composition in the survey area, which is determined by on-site survey analysis.

[0017] Further, the The node values are taken in units of one meter, and The electrodes are evenly arranged in the placement coordinate set a range in sequence, and the electrodes are preferentially placed The numerical value is small until The electrodes are placed in the placement coordinate set a range in sequence, and the electrodes are preferentially placed The node values are taken in units of half a meter, and The electrodes are evenly arranged in the placement coordinate set b range in sequence, and the electrodes are preferentially placed The numerical value is small until The electrodes are placed in the placement coordinate set b range in sequence.

[0018] Further, the values of the corresponding vertical coordinates of the same horizontal coordinates between the three adjacent horizontal layered distribution maps i above and below are taken as three control points, and a smooth curve is obtained using the Bezier curve formula ; the values of the corresponding vertical coordinates of the same horizontal coordinates between the horizontal layered distribution map u and the adjacent horizontal layered distribution maps i above and below are taken as three control points, and a smooth curve is obtained using the Bezier curve formula The difference between the coefficient of the new smooth curve and the smooth curve is taken as the deviation coefficient of the prediction model , and the deviation coefficient is multiplied by the values of the corresponding vertical coordinates of the same horizontal coordinates in the three closest horizontal layered distribution maps i above and below, respectively, to obtain new three control points, and a smooth curve is calculated using the Bezier curve formula ; the values of the corresponding vertical coordinates of the same horizontal coordinates between the horizontal layered distribution map o and the adjacent horizontal layered distribution maps u above and below are taken as three control points, and a smooth curve is obtained using the Bezier curve formula The difference between the coefficient of the new smooth curve and the smooth curve is taken as the deviation coefficient of the prediction model , and the deviation coefficient is multiplied by the values of the corresponding vertical coordinates of the same horizontal coordinates in the three closest horizontal layered distribution maps u above and below, respectively, to obtain new three control points, and a smooth curve is calculated using the Bezier curve formula The Bezier curve formula is a mathematical method for defining a curve through a series of control points, which can accurately control the shape of the curve by adjusting the position of the control points, and is used for data deviation correction.

[0019] Further, the processing unit adopts constant current power supply square wave to sequentially power the arranged electrodes, and calculates to generate horizontal layered distribution map u, wherein the power supply current of the processing unit is , calculates to generate horizontal layered distribution map o, wherein the power supply current of the processing unit is .

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

[0021] 1. The coal mine parallel electrical method data acquisition and processing system can refine and enlarge the area with complex underground composition by reducing the detection range in batches, can effectively make up the details of the generated model, and make the model features more intuitive and obvious.

[0022] 2. The coal mine parallel electrical method data acquisition and processing system can effectively reduce the error probability of the full-space layered model and improve the accuracy and authenticity of the generated model by introducing a deviation coefficient verification link to feed back and correct the prediction model.

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

[0024] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the following will briefly introduce the drawings needed to be used in the embodiment description. 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.

[0025] Fig. 1 The coal mine parallel electrical method data acquisition and processing system block diagram of the present application;

[0026] Fig. 2 The flowchart of the present application;

[0027] Fig. 3 The stacking diagram of the full-space layered model of the present application. DETAILED DESCRIPTION

[0028] The technical solutions in the embodiments of the present application will be described clearly and completely below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only some embodiments of the present application, not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor are within the scope of protection of the present application.

[0029] Please refer to Figs. 1-3 The present application provides a technical solution: a coal mine parallel electrical method data acquisition and processing system, as shown in Fig. 1As shown, it includes n electrodes and m acquisition terminals, as well as a processing unit and a prediction model;

[0030] Will One electrode is horizontally and evenly distributed within the survey area, and m acquisition terminals are horizontally and evenly distributed within the survey area.

[0031] The processing unit sequentially supplies power to the deployed electrodes. Each time an electrode is powered, m acquisition terminals simultaneously collect potential information at the corresponding location and transmit it to the processing unit. The processing unit calculates the resistivity based on the potential difference between each acquisition terminal. The calculated resistivity is then transmitted to the prediction model. The prediction model generates horizontal layered distribution maps i with different depths and spacings of q based on the resistivity calculated from n electrodes. Finally, a smooth curve is generated between each horizontal layered distribution map i using the Bézier curve formula. The prediction model determines the set of coordinates 'a' for electrode placement in the next stage based on the resistivity dividing lines of the horizontal layered distribution map i of each layer. The specific representation of set a is as follows:

[0032]

[0033] The coordinates of the electrodes in set a;

[0034] The processing unit sequentially feeds the coordinate set a into the coordinate set. Each electrode is powered, and each time the electrode is powered, m acquisition terminals simultaneously acquire the potential information at the corresponding location and transmit it to the processing unit. The processing unit calculates the resistivity based on the potential difference between each acquisition terminal, and transmits the calculated resistivity to the prediction model. The prediction model then... The resistivity generation spacing calculated for each electrode is: The horizontal layered distribution map u at different depth sections is interspersed within the horizontal layered distribution map i. The prediction model compares the horizontal layered distribution map u with the corresponding smooth curve within this range. Analysis and calculation of the deviation coefficient and the deviation coefficient Inputting the data into the horizontal layered distribution map i generates a smooth curve. The prediction model determines the coordinate set b for electrode placement in the next stage based on the resistivity dividing lines of the horizontal layered distribution map u of each layer. The specific representation of set b is as follows:

[0035]

[0036] The coordinates of the electrodes in set b;

[0037] The processing unit sequentially feeds the coordinate set b into the coordinate set. electrodes, the m acquisition terminals collect the potential information of the corresponding positions at the same time and transmit the information to the processing unit, the processing unit calculates the potential difference between each acquisition terminal resistivity, the processing unit transmits the calculated resistivity to the prediction model, the prediction model calculates the interval between the two electrodes according to the calculated resistivity The prediction model generates a smooth curve according to the horizontal stratified distribution diagram o and the corresponding smooth curve in the range The prediction model calculates the deviation coefficient and the deviation coefficient is brought into the Bezier curve formula between each horizontal stratified distribution diagram o to generate a smooth curve As shown in Figs. 2-3 The prediction model combines the horizontal stratified distribution diagram i, the horizontal stratified distribution diagram u, and the horizontal stratified distribution diagram o and uses the smooth curve to fit the gap between the stratified distribution diagrams to generate a full-space stratified model; by reducing the detection range in batches, the area with complex underground composition can be refined and enlarged, which can effectively supplement the details of the generated model and make the model features more intuitive and obvious.

[0038] In parallel electrical exploration, the formula for calculating resistivity is as follows:

[0039] , is the resistivity, is the potential difference measured between the two acquisition terminals, is the device coefficient, is the power supply current of the electrodes arranged in the generated horizontal stratified distribution diagram i; the device coefficient is related to the position of the arranged electrodes, when the position of the arranged electrodes is determined, the specific value of the device coefficient can be determined.

[0040] wherein the interval q is determined by the area of the on-site survey range and the electrode arrangement interval density, the resistivity dividing line is the dividing line of different range resistivities, indicating the relationship between rock conductivity and mineral composition in the survey range, which is determined by on-site survey analysis.

[0041] wherein the node is valued in units of one meter, and electrodes are evenly arranged in the placement coordinate set a range, the electrodes are preferentially placed in the numerical value smaller coordinate until all the electrodes are placed, ​The node is taken as a unit of half a meter, and The electrodes are evenly arranged in the placement coordinate set b range in turn, and the electrodes are preferentially placed The numerical value is small, and the coordinate is placed until The electrodes are placed completely.

[0042] Among them, the values of the corresponding vertical coordinates of the same horizontal coordinate between the three adjacent horizontal layer distribution maps i are taken as three control points, and a smooth curve is obtained by using the Bezier curve formula ; the values of the corresponding vertical coordinates of the same horizontal coordinate between the horizontal layer distribution map u and the adjacent horizontal layer distribution map i are taken as three control points, and a smooth curve is obtained by using the Bezier curve formula The coefficient, the new smooth curve The difference between the coefficient and the smooth curve is taken as the deviation coefficient of the prediction model , the deviation coefficient is multiplied by the values of the corresponding vertical coordinates of the same horizontal coordinate in the three closest horizontal layer distribution maps i to obtain new three control points, and a smooth curve is obtained by using the Bezier curve formula ; the values of the corresponding vertical coordinates of the same horizontal coordinate between the horizontal layer distribution map o and the adjacent horizontal layer distribution map u are taken as three control points, and a smooth curve is obtained by using the Bezier curve formula The coefficient, the new smooth curve The difference between the coefficient and the smooth curve is taken as the deviation coefficient of the prediction model , the deviation coefficient is multiplied by the values of the corresponding vertical coordinates of the same horizontal coordinate in the three closest horizontal layer distribution maps u to obtain new three control points, and a smooth curve is obtained by using the Bezier curve formula ; the Bezier curve formula is a mathematical method for defining a curve through a series of control points, and the shape of the curve can be accurately controlled by adjusting the positions of the control points, which is used for correcting the data deviation; the introduction of the deviation coefficient verification link can effectively reduce the error probability of the full-space layer model and improve the accuracy and authenticity of the generated model.

[0043] Among them, the processing unit adopts constant current power supply square wave to supply power to the arranged electrodes in turn, the processing unit power supply current in the calculation and generation of the horizontal layer distribution map u is , the processing unit power supply current in the calculation and generation of the horizontal layer distribution map o is The initial survey range is gradually reduced by comparing the set a and the set b, and the demand for direct current power supply is also correspondingly reduced, which is more conducive to the accuracy of the later model.

[0044] The above merely describes preferred specific embodiments of the present application, but the protection scope of the present application is not limited thereto, and any person skilled in the art, according to the technical solution and inventive concept of the present application, makes equivalent replacement or change within the technical range disclosed by the present application, which should be covered within the protection scope of the present application.

Claims

1. A coal mine parallel electrical method data acquisition and processing system, comprising n electrodes and m acquisition terminals, wherein the n electrodes are evenly arranged in the survey range, and the m acquisition terminals are evenly arranged in the survey range; characterized in that: Also include processing unit and prediction model; ​ The processing unit sequentially supplies power to the electrodes, and each time the electrodes are powered, m acquisition terminals collect potential information at corresponding positions at the same time 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 generates horizontal layered distribution maps i of different depth sections with a spacing of q according to the resistivity calculated by the electrodes , and generates a smooth curve between each horizontal layered distribution map i according to the Bezier curve formula , and determines the next stage electrode placement coordinate set a according to the horizontal layered distribution map i resistivity split line of each layer. The processing unit sequentially supplies power to a in the coordinate set a The acquisition terminal collects the potential information and transmits it to the processing unit, calculates the resistivity according to the potential difference between each acquisition terminal, transmits the calculated resistivity to the prediction model, and the prediction model generates a horizontal layered distribution map u with different depth sections with a distance of The horizontal layered distribution map u is inserted into the horizontal layered distribution map i, selects d points in the horizontal layered distribution map u and the corresponding points of the smooth curve in the range Analyzes and calculates the deviation coefficient The calculation formula of the deviation coefficient is , where is the value of the ordinate of the selected point in the horizontal layered distribution map u, is the value of the ordinate of the corresponding point of the deviation coefficient , and the deviation coefficient is brought into the horizontal layered distribution map i to generate a smooth curve According to the horizontal layered distribution map u resistivity segmentation line of each layer, determine the next stage electrode placement coordinate set b,​ The processing unit sequentially assigns values ​​to the coordinate set b. Each electrode is powered, and the acquisition terminal acquires potential information and transmits it to the processing unit. The resistivity is calculated based on the potential difference between each acquisition terminal, and the calculated resistivity is transmitted to the prediction model. The prediction model then... The resistivity generation spacing calculated for each electrode is: A horizontal layered distribution map o at different depths is interspersed within a horizontal layered distribution map u. G points in the horizontal layered distribution map o are selected, and their corresponding smooth curves within that range are compared. point Analysis and calculation of the deviation coefficient Deviation coefficient The calculation formula is: ,in The value of the ordinate of a selected point in the horizontal layered distribution map o. Deviation coefficient The value of the ordinate of the corresponding point, and the deviation coefficient. By substituting the Bézier curve formula between the various horizontal layered distribution maps o, a smooth curve is generated. Combine the horizontal layered distribution map i, horizontal layered distribution map u, and horizontal layered distribution map o respectively, and use a smooth curve. A full-space layered model is generated by fitting the gaps between the layered distribution maps.

2. The coal mine parallel electrical method data acquisition and processing system according to claim 1, characterized in that, In parallel electrical prospecting, the formula for calculating the resistivity is as follows: , is the resistivity, is the potential difference measured between two acquisition terminals, is the device coefficient, is the supply current of the electrodes laid out in the horizontal stratigraphic profile i calculated. The specific representation of the set a is: ， The To position the electrode placement coordinates in the collection a; The specific representation of the set b is: ， The is the position of the electrode placement coordinates in the collection b.

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

4. The coal mine parallel electrical method data acquisition processing system according to claim 1, characterized in that, The Values ​​are taken at nodes with a unit length of one meter. The electrodes are arranged horizontally and uniformly within the set of coordinates a, with priority given to placing the electrodes in this location. The numerical value is too small, until the coordinates are too small. After all electrodes have been placed, the... Values ​​are taken at nodes with a unit length of half a meter. The electrodes are arranged horizontally and uniformly within the set of coordinates b, with priority given to placing the electrodes in this location. The numerical value is too small until the coordinates are too small. All electrodes have been placed.

5. The coal mine parallel electrical method data acquisition processing system according to claim 1, characterized in that, The processing unit adopts constant current power supply square wave to supply power to the arranged electrodes in turn, and calculates to generate the horizontal layered distribution map u, wherein the power supply current of the processing unit is , and calculates to generate the horizontal layered distribution map o, wherein the power supply current of the processing unit is .

Images