Disclosure of Invention
1. Problems to be solved
Aiming at the problems that the existing underground tensor resistivity measurement method is incomplete and the accuracy and comprehensiveness of the final measurement result are difficult to ensure and errors possibly occur in the prediction of geological disasters, the invention provides the underground tensor resistivity measurement method which can comprehensively measure the resistivity along the trend, the tendency and the vertical direction of the rock stratum and improve the accuracy of the final measurement result so as to accurately predict the geological disasters.
2. Technical proposal
In order to solve the problems, the invention adopts the following technical scheme.
A method of measuring downhole tensor resistivity, the method comprising the steps of:
1. N substation nodes are arranged at intervals along the roadway construction direction on the underground rock layer, and three transmitting electrodes and three collecting electrodes are respectively arranged on each substation node along the X, Y and Z directions;
2. the master station equipment sends an acquisition command to the substation equipment, the transmitting electrode of the substation node transmits a signal, the acquisition electrode receives the signal and records the potential difference between the transmitting electrode and the acquisition electrode, and then the substation equipment uploads the acquired potential difference to the master station equipment;
3. the master station equipment gathers the collected data and uploads the collected data to the background equipment, and the background equipment calculates the resistivity according to the collected potential difference data among all the nodes.
As a further improvement of the technical scheme, the method also comprises substation node working state inspection, and the specific process is as follows:
S1, sorting a patrol queue: numbering substation nodes along the roadway construction direction and setting an inspection-free queue and an inspected queue;
S2, judging the node state: each substation node is patrolled, the patrolled substation node is added into the patrolled queue, and the non-patrolled substation node is added into the non-patrolled queue; then carrying out next inspection, if the substation node which is not inspected to the queue is inspected, transferring the substation node to the inspected queue; and after the inspection for a plurality of times, the substation nodes which are not inspected in the queue are listed as fault nodes, and the substation nodes which are inspected in the queue are listed as effective nodes.
As a further improvement of the technical scheme, the substation node working state inspection further includes:
S3, calculating the mapping relation between the logical address and the physical address: pairing the logical address of the substation node which is patrolled and checked in the queue with the physical address, wherein the logical address is the sorting number of the substation node in the patrolled and checked in the queue, and the physical address is the number of the substation node in step S1.
As a further improvement of the technical scheme, the method also comprises substation node acquisition state inspection, and the specific process is as follows:
S1, collecting state queue classification: setting three state queues of acquisition completion, non-acquisition and error, and adding all substation nodes into the non-acquisition queue at the beginning;
S2, data acquisition: taking out the substation nodes in the non-collected queue for data collection, adding the substation nodes with the collected data within a set range into a collection completion queue, and adding the substation nodes with the collected data not within the set range into an error queue;
S3, inspection of error nodes: the data acquisition is carried out on the substation nodes in the error queue, if the acquired data is within a set range, the substation nodes are added into an unaacquired queue and acquired again, and if the acquired data is still within the set range, the substation nodes are added into an acquisition completion queue;
S4, repeating the steps S1 to S3 to obtain the acquisition states of all the substation nodes.
As a further improvement of the technical scheme, in the second step, the specific process of collection is as follows:
s1, designating a substation node as a transmitting node;
S2, switching on one transmitting electrode in the transmitting node;
s3, taking the rest substation nodes as acquisition nodes, and circularly connecting acquisition electrodes in each acquisition node to acquire data;
S4, sequentially connecting the rest transmitting electrodes in the transmitting node and repeating the step S3;
s5, sequentially designating all the substation nodes as transmitting nodes and repeating the steps S2 to S4 until all the substation nodes are used as primary transmitting nodes.
As a further improvement of the technical scheme, if the substation nodes exist in the error queue, the substation nodes in the error queue are not used as transmitting nodes and collecting nodes in the collecting process.
As a further improvement of the technical solution, the collecting process further includes:
s6, converting the format of the acquired data into a protocol definition data format meeting the communication transmission requirement.
As a further improvement of the technical scheme, the second step further comprises a processing process of collecting data:
s1, respectively inputting data acquired by an acquisition electrode of each acquisition node in each acquisition process into a waveform generator, and inputting a set direct current level into the waveform generator to generate potential difference waveforms of positive level, zero level and negative level with set pulse width;
s2, setting the signal emission times in the primary acquisition process as N, wherein the generated waveform comprises N periods, and each period comprises potential difference waveforms of positive level, zero level and negative level;
S3, collecting 5-20 data in the positive level waveform in each period and taking an average value as a positive level measurement value HV i;
S4, collecting 5-20 data in the negative level waveform in each period and taking an average value as a negative level measured value LV i;
S5, obtaining an average value of HV i in N periods as a positive level result value HV AVG, and obtaining an average value of LV i in N periods as a negative level result value HV AVG;
s6, calculating the component potential difference of the transmitting electrode and the collecting electrode to be U= (|HV AVG|+|HLAVG |) and 2.
As a further improvement of the technical scheme, in the third step, the tensor resistivity is calculated by adopting the following formula;
E=U/d;
ρ=E/J;
Wherein E is the electric field strength, d is the distance from the transmitting electrode to the collecting electrode, J is the current density, and ρ is the tensor resistivity.
As a further improvement of the technical scheme, if the data volume of the substation node collected by the master station device exceeds a set value, the data is transmitted in a subpacket mode:
wherein TRANSDATA is a transmission quantity, ALLDatas is acquired data, and K is a sub-packet number.
3. Advantageous effects
Compared with the prior art, the invention has the beneficial effects that:
(1) According to the method for measuring the resistivity of the underground tensor, the resistivity along the trend, the tendency and the vertical direction of the rock stratum can be comprehensively measured, the accuracy of a final measurement result is improved, so that geological disasters are accurately predicted and monitored in real time, premonitory data of the inoculation and development of the geological disasters can be timely acquired, and the stress peaks are analyzed and predicted in real time, so that early warning is realized, and the reversible stress of dynamic disasters is timely released;
(2) According to the method for measuring the underground tensor resistivity, disclosed by the invention, the fault node and the node with abnormal acquired data can be effectively screened by setting the substation node working state inspection and the substation node acquisition state inspection, so that the fault node and the node with abnormal acquired data are removed during actual data acquisition, effective data are reserved, the accuracy of a final calculation result is ensured, and in addition, the detected sequence number of the workstation substation node corresponds to the actual physical address number of the substation node, so that a worker can clearly master the calculation result corresponding to each substation node, and confusion is avoided;
(3) According to the method for measuring the underground tensor resistivity, the potential difference waveforms are generated by the acquired data through the waveform generator, the average value of the positive level and the negative level is respectively calculated, and the component potential difference between the transmitting electrode and the acquisition electrode is calculated, so that the accuracy of the finally calculated tensor resistivity is effectively improved, and the ground disaster can be accurately predicted and monitored in real time.
Detailed Description
The invention is further described below in connection with specific embodiments and the accompanying drawings.
Example 1
The method is used for measuring and calculating tensor resistivity of each trend of rock stratum in underground mine, especially coal mining process, and timely predicting geological disasters such as coal mine dynamic disasters and the like, ensuring normal running of underground mining work and personnel safety, and describing specific process in detail below.
The method adopts a tensor resistivity complete device, and mainly comprises main station equipment, substation equipment, an industrial personal computer, and matched communication interfaces and communication devices. The master station equipment comprises a power supply, a control chip, a CAN interface circuit, a W5500 chip Ethernet communication interface circuit and an RJ45 network interface module. The substation equipment comprises a power supply, a multi-path emission power supply output, three paths of emission electrodes, three paths of acquisition electrodes, a control chip ARM and a CAN interface circuit. The industrial personal computer is provided with an industrial Ethernet and Internet interface and is used for receiving data transmitted by the master station equipment and controlling the master station equipment and the substation equipment to work.
As shown in fig. 6, the method constructs an integral monitoring system for address disasters, the industrial personal computer controls the main station equipment and the substation equipment to work, the substation equipment collects component potential differences of all substation nodes, the collected data are uploaded to the main station equipment for integration, the background equipment calculates tensor resistivity according to the collected data, and then staff predicts geological disasters according to calculation results to realize real-time monitoring of geological disasters such as coal mine dynamic disasters. Specifically, the method comprises the following steps:
1. N sub-station nodes are arranged on the underground rock layer surface along the roadway construction direction at intervals, and each sub-station node is sequentially numbered, and 32 sub-station nodes are arranged in the embodiment. In general, the intervals among the substation nodes are kept as consistent as possible, and the specific number of the substation nodes is set according to the roadway length. Each substation node is provided with substation equipment, and three transmitting electrodes and three collecting electrodes are respectively arranged along the directions X, Y and Z, and particularly, the directions X, Y and Z are directions along the trend, the trend and the vertical direction of the rock stratum. As shown in fig. 1, is a working schematic diagram of a substation node. Between two sub-station nodes, three transmitting electrodes of one sub-station respectively transmit signals, and three collecting electrodes of the other sub-station node respectively receive signals, so that 9 component potential differences between the two sub-station nodes in different directions, namely { XX, YX, ZX }, { XY, YY, ZY }, { XZ, YZ, ZZ }, can be measured between the three transmitting electrodes and the three collecting electrodes.
For example, the number of the substation node which is acquired by a certain detection and serves as the transmitting node is 3, the number of the substation node which is acquired as the acquiring node is 12, and the acquired component potential difference value is:
{{7.9112361E-02,8.5602672E-02,8.3562645E-02},{6.4245701E-03,6.1099785E-03,5.8294712E-03},{3.5230172E-03,3.0151256E-03,3.2260576E-03}}.
2. The master station device sends an acquisition command to the substation device, the transmitting electrode of the substation node transmits a signal, the acquisition electrode receives the signal and records the potential difference between the transmitting electrode and the acquisition electrode, and then the substation device uploads the acquired potential difference to the master station device.
Before collection starts, the working states of the substation nodes need to be inspected, namely, whether substation equipment of all the substation nodes has faults or not is judged, as shown in fig. 3, the specific process is as follows:
s1, sorting a patrol queue: the substation nodes are numbered along the roadway construction direction and are provided with an unground queue unPollingQueue and a round queue pollingQueue.
S2, judging the node state: each substation node is patrolled, the patrolled substation node is added into the patrolled queue, and the non-patrolled substation node is added into the non-patrolled queue; then carrying out next inspection, if the substation node which is not inspected to the queue is inspected, transferring the substation node to the inspected queue; and after the inspection for a plurality of times, the substation nodes which are not inspected in the queue are listed as fault nodes, and the substation nodes which are inspected in the queue are listed as effective nodes. The number stored at this time is the physical address number of the substation node in S1.
For example, in this step of a certain inspection:
unPollingQueue={3,12}
pollingQueue={1,2,4,5,6,7,8,9,10,11,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32}
S3, calculating the mapping relation between the logical address and the physical address: the logical addresses and the physical addresses of the substation nodes which are patrolled and checked in the queue are paired to form the corresponding relation ADDRMAPPING < logical address, physical address > of the logical and physical addresses of all the substation nodes. The logical address is the sorting number of the substation node in the sorted queue, and the physical address is the number of the substation node in step S1.
addrMapping={<1,1>,<2,4>,<3,5>,<4,6>,<5,7>,<6,8>,<7,9>,<8,10>,<9,11>,<10,13>,<11,14>,<12,15>,<13,16>,<14,17>,<15,18>,<16,19>,<17,20>,<18,21>,<19,22>,<20,23>,<21,24>,<22,25>,<23,26>,<24,27>,<25,28>,<26,29>,<27,30>,<28,31>,<29,32>,<30,33>}.
I.e. the substation equipment with physical address numbers 3 and 12 fails in this patrol.
Next, as shown in fig. 4, the collecting state of the operable substation node is inspected, and it is determined whether an abnormality occurs in the collecting process, and the inspecting process is as follows:
S1, collecting state queue classification: three state queues of acquisition completion finishedQueue, non-acquisition unFinishedQueue and error errorQueue are set, and all substation nodes are added to the non-acquisition queue at the beginning. In a certain inspection process, a substation node with a physical address number of 1 is used as a transmitting node.
Then unFinishedQueue Start to ={2,4,5,6,7,8,9,10,11,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32}
S2, data acquisition: and taking out the substation nodes in the non-collected queue for data collection, adding the substation nodes with the collected data within a set range into a collection completion queue, and adding the substation nodes with the collected data not within the set range into an error queue.
S3, inspection of error nodes: and (3) carrying out data acquisition on the substation nodes in the error queue, if the acquired data is within a set range, adding the substation nodes into an unaacquired queue and acquiring again, and if the acquired data is still within the set range, adding the substation nodes into an acquisition completion queue.
S4, repeating the steps S1 to S3 to obtain the acquisition states of all the substation nodes. The number of the inspection repetition is 6, and finally the method is as follows:
unFinishedQueue Ending ={20}
finishedQueue={2,4,5,6,7,8,9,10,11,13,14,15,16,17,18,19,21,22,23,24,25,26,27,28,29,30,31,32}
errorQueue={}
Namely, the substation node with the physical address number of 20 still does not meet the requirement after 6 times of polling, so the substation node is set as an unfinished node, data can still be uploaded to an industrial personal computer, and further analysis reasons include unstable work, undervoltage, poor grounding and the like.
By setting the substation node working state inspection and substation node acquisition state inspection, the nodes with faults and abnormal acquired data can be effectively screened, so that the faults are discharged during actual data acquisition, the abnormal acquired data are removed during calculation, effective data are reserved, and the accuracy of a final calculation result is ensured. In addition, the detected sequence numbers of the workstation sub-station nodes are corresponding to the actual physical address numbers of the sub-station nodes, so that a worker can clearly grasp the calculation results corresponding to the sub-station nodes, and confusion is avoided.
As shown in fig. 5, the specific process of acquisition is as follows:
s1, designating a substation node as a transmitting node.
S2, one transmitting electrode in the transmitting node is connected.
S3, taking the rest substation nodes as acquisition nodes, and circularly connecting acquisition electrodes in each acquisition node to acquire data.
S4, sequentially switching on the rest transmitting electrodes in the transmitting node and repeating the step S3.
S5, sequentially designating all the substation nodes as transmitting nodes and repeating the steps S2 to S4 until all the substation nodes are used as primary transmitting nodes.
S6, converting the format of the acquired data into a protocol definition data format meeting the communication transmission requirement.
It should be noted that if there is a substation node in the error queue, the substation node in the error queue is not used as a transmitting node and a collecting node in the collecting process.
Specifically, one of the 32 substation nodes is selected as a transmitting node i, and one of the transmitting electrodes in the transmitting node i is turned on as a transmitting electrode h at this time. And then sequentially enabling the rest substation nodes to serve as acquisition nodes j, and sequentially switching on three acquisition electrodes in the acquisition nodes j to serve as acquisition electrodes k at the time. When all the substation nodes except i are used as acquisition nodes, and potential difference data are acquired by the acquisition electrodes of the substation nodes, the transmission electrode h is circularly switched in i, and the acquisition steps of the rest substation nodes are repeated until all the three transmission electrodes finish one-time signal transmission. And switching the substation nodes and repeating the steps until all the substation nodes serve as one-time transmitting nodes, and completing one-time complete data acquisition process.
It should be noted that, not every acquisition process is necessarily complete, but only the potential difference data of some substation nodes can be acquired separately, so that the acquisition process has various situations as follows:
Case 1: when the acquisition is not specified, i, h, j and k are acquired in a full cycle mode in an automatic mode, and a complete acquisition process is adopted.
Case 2: when the transmitting node i is designated, h, j and k are acquired in a full cycle.
Case 3: when the transmitting node i and the transmitting electrode h are designated, j and k are acquired in a full cycle.
Case 4: when the transmitting node i, the transmitting electrode h and the collecting node j are designated, k is collected in a full cycle.
Case 5: when the transmitting node i and the transmitting electrode h, and the collecting node j and the collecting electrode k are designated, only the component potential difference between the transmitting electrode h and the collecting electrode k is collected.
In order to improve the accuracy of the final calculated tensor resistivity, so as to accurately predict and monitor the ground disaster in real time, the embodiment further includes a process of collecting data, as shown in fig. 2, and the specific process is as follows:
S1, respectively inputting data acquired by the acquisition electrodes of each acquisition node in each acquisition process into a waveform generator, and inputting a set direct current level into the waveform generator to generate potential difference waveforms of positive level, zero level and negative level with set pulse width.
S2, setting the signal emission times in the primary acquisition process as N, wherein the generated waveform comprises N periods, and each period comprises potential difference waveforms of positive level, zero level and negative level.
S3, collecting 5-20 data in the positive level waveform in each period and taking an average value as a positive level measurement value HV i.
S4, collecting 5-20 data in the negative level waveform in each period and taking an average value as a negative level measured value LV i.
It should be noted that, the positive level pulses and the negative level pulses in the period T respectively include three time periods of a start time, a measurement time and an end time, and the data in the measurement time is generally accurate, so that the data in the measurement time are obtained in step S3 and step S4.
S5, the average value of HV i in N periods is obtained as a positive level result value HV AVG, and the average value of LV i in N periods is obtained as a negative level result value HV AVG.
S6, calculating the component potential difference of the transmitting electrode and the collecting electrode to be U= (|HV AVG|+|HLAVG |) and 2.
Finally, calculating tensor resistivity by adopting the following formula;
E=U/d;
ρ=E/J;
Wherein E is the electric field strength, d is the distance from the transmitting electrode to the collecting electrode, J is the current density, and ρ is the tensor resistivity.
It should be noted that if the data volume of the substation node collected by the master station device exceeds a set value, the data needs to be transmitted in a subpacket mode, otherwise, a packet loss phenomenon occurs, and the final calculation result is not accurate enough.
Wherein TRANSDATA is a transmission quantity, ALLDatas is acquired data, K is the number of sub-packets, the calculation result is not an integer, and an integer value of a little bigger is taken.
For example, if the acquisition data ALLDatas is 720 bytes and the transmission amount TRANSDATA is 512 bytes, k+|720/512|=2.
In summary, according to the method for measuring the downhole tensor resistivity in the embodiment, the resistivity along the trend, the tendency and the vertical direction of the rock stratum can be comprehensively measured, and the accuracy of the final measurement result is improved, so that geological disasters are accurately predicted, early state data of the inoculation and development of the geological disasters can be timely acquired, and the stress peaks can be analyzed and predicted in real time, so that early warning is realized, and the reversible stress of dynamic disasters is timely released.
The examples of the present invention are merely for describing the preferred embodiments of the present invention, and are not intended to limit the spirit and scope of the present invention, and those skilled in the art should make various changes and modifications to the technical solution of the present invention without departing from the spirit of the present invention.