CN113885085A - Underground direct-current axial dipole dynamic source advanced detection method - Google Patents

Underground direct-current axial dipole dynamic source advanced detection method Download PDF

Info

Publication number
CN113885085A
CN113885085A CN202110896200.XA CN202110896200A CN113885085A CN 113885085 A CN113885085 A CN 113885085A CN 202110896200 A CN202110896200 A CN 202110896200A CN 113885085 A CN113885085 A CN 113885085A
Authority
CN
China
Prior art keywords
distance
electrodes
formula
dipole
receiving
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110896200.XA
Other languages
Chinese (zh)
Other versions
CN113885085B (en
Inventor
石显新
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Research Institute Co Ltd of CCTEG
Original Assignee
Xian Research Institute Co Ltd of CCTEG
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Research Institute Co Ltd of CCTEG filed Critical Xian Research Institute Co Ltd of CCTEG
Priority to CN202110896200.XA priority Critical patent/CN113885085B/en
Publication of CN113885085A publication Critical patent/CN113885085A/en
Application granted granted Critical
Publication of CN113885085B publication Critical patent/CN113885085B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
    • G01V3/26—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Electromagnetism (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses an advanced detection method of an underground direct current axial dipole dynamic source, wherein a receiving electrode M or receiving electrodes M and N are fixedly arranged on a tunneling surface, and transmitting electrodes A and B are arranged from A1And B1The part moves along the central line of the roadway floor to the direction of the driving surface to form a position A1And B1，A2And B2，...，AiAnd Bi...AnAnd BnA dipole moving source of location points; the transmitting electrodes a and B transmit a current once at each location point and the receiving electrode M or receiving electrodes M and N measure the voltage once. The advanced detection maximum detection distance estimation formula provided by the invention establishes the relationship among electrode layout, emission current, formation resistivity and noise level, and provides quantitative basis for selection of underground direct current single-pole moving source advanced detection construction parameters; given a corresponding toiDepth measurement point D ofiAnd apparent resistivity
Figure DDA0003198019140000011
The formula is an algorithm with specificity; and error records in the measured data are used as the judgment standard of abnormal response, so that the reliability guarantee is provided for data interpretation. The method of the invention improves the accuracy of the underground direct current advanced detection.

Description

Underground direct-current axial dipole dynamic source advanced detection method
Technical Field
The invention belongs to the field of electrical and electromagnetic prospecting, relates to an underground direct current advanced detection technology, and particularly relates to an underground direct current axial dipole dynamic source advanced detection method.
Background
The underground direct current advanced detection is one of the main methods for detecting abnormal bodies in front of the driving face of the underground coal mine in China. In the current dc advanced detection method, a monopole transmitting-dipole receiving device called a tripolar method is generally adopted, a transmitting electrode is fixedly arranged on or near a heading face, and a receiving electrode moves away from the heading face. Thereafter, there has been established a multi-monopole-dipole combination method improved on the basis of the three-pole method. Due to the nature of the close-acting DC electric field, the response of the anomaly upon excitation by the emission source is always transmitted outward from the periphery of the anomaly by the field, with the intensity decreasing with increasing distance from the anomaly. By fixing the transmitting electrode near the tunneling surface and moving the receiving electrode to the opposite direction of the tunneling working surface, the obtained abnormity is weak, and the abnormity body behind the tunneling surface can interfere with the advanced detection to cause misjudgment.
Disclosure of Invention
Aiming at the defects and shortcomings in the prior art, the invention provides an advanced detection method for an underground direct current axial dipole dynamic source, which is characterized in that a receiving electrode is fixedly arranged on a tunneling working surface and is close to a front abnormal body as much as possible, and the advanced detection is realized by moving a transmitting electrode to the direction of the tunneling working surface.
In order to achieve the purpose, the invention adopts the following technical scheme:
a method for detecting the active source of underground DC axial dipole includes such steps as fixing the receiving electrodes on the heading surface, locating the receiving electrodes on the vertical bisector of the bottom edge of heading surface, and using the emitting electrodes A and B to emit the active source from A1And B1The part moves along the central line of the roadway floor to the direction of the driving surface to form a position A1And B1Position point, A2And B2Location points, aiAnd BiPoints of position, anAnd BnA dipole moving source of location points; transmitting current once by the transmitting electrodes A and B at each position point, measuring voltage once by the receiving electrode M or the receiving electrodes M and N correspondingly, and performing dipole moving source advanced detection;
when the monopole receiving electrode M is fixedly arranged on the heading face, the estimation formula of the maximum advanced detection distance is as follows:
Figure BDA0003198019120000021
in the above formula, O'1O is the midpoint of emitter electrode distance AB'1Distance to the middle point O of the bottom edge of the driving face, IABmaxIn order to achieve the maximum emission current,
Figure BDA0003198019120000022
the rho is the resistivity of the rock stratum in front of the driving face and is the noise level observed by the receiving electrode M;
when the dipole receiving electrodes M and N are fixedly arranged on the tunneling surface, the estimation formula of the maximum detection distance of the advanced detection is as follows:
Figure BDA0003198019120000023
in the above formula, O'1O is the midpoint of emitter electrode distance AB'1Distance to the midpoint O of the bottom edge of the heading face, O 'O is the distance from the receiving electrode to the midpoint O' of MN, IABmaxIn order to achieve the maximum emission current,
Figure BDA0003198019120000024
ρ is the resistivity of the formation in front of the face for the noise level observed by the receiver electrodes M and N.
The invention also comprises the following technical characteristics:
optionally, A is1And B1Position point, A2And B2Location points, aiAnd BiPoints of position, anAnd BnThe smaller the interval between the position points is, the higher the detection resolution is; emitter electrodes A and B from A1And B1The position of the current field moves to the direction of the heading face at intervals required by resolution along the center line of the roadway bottom plate, and the depth of the current field penetrating into the front of the heading face is from shallow to deep.
Optionally, when the monopole receiving electrode M is fixedly arranged on the heading face, the monopole receiving electrode M and the transmitting electrodes A and B are arranged at AiAnd BiAt the corresponding sounding point DiAnd apparent resistivity
Figure BDA0003198019120000029
The formula of (1) is:
Di≈O′1O-O′iO，(i＝1,2,...,n) (2a)
o 'of the above formula'1O is determined by formula (1a), O'iO is the midpoint of emitter electrode distance AB'iDistance to O;
Figure BDA0003198019120000025
in the above formula
Figure BDA0003198019120000026
Is AiThe distance to M,
Figure BDA0003198019120000027
Is BiA distance to M, wherein AiO is AiDistance to O, BiO is BiDistance to O, MO is the distance from M to O,
Figure BDA0003198019120000028
the emitting electrodes A and B are at AiAnd BiThe emission current of the light source (c),
Figure BDA0003198019120000031
is and
Figure BDA0003198019120000032
a corresponding observed voltage;
when the dipole receiving electrodes M and N are fixedly arranged on the heading face, the dipole receiving electrodes M and N are in A connection with the transmitting electrodes A and BiAnd BiAt the corresponding sounding point DiAnd apparent resistivity
Figure BDA0003198019120000033
The formula of (1) is:
Di≈O′1O-O′iO，(i＝1,2,...,n) (2c)
of which O'1O is determined by formula (1b), O'iO is the midpoint of emitter electrode distance AB'iDistance to O;
Figure BDA0003198019120000034
in the formula
Figure BDA0003198019120000035
AiThe distance to M,
Figure BDA0003198019120000036
Is BiThe distance to M,
Figure BDA0003198019120000037
AiThe distance to N,
Figure BDA0003198019120000038
Is BiA distance to N, wherein AiO is AiDistance to O, BiO is BiDistance to O, MO is the distance from M to O, NO is the distance from N to O,
Figure BDA0003198019120000039
the emitting electrodes A and B are at AiAnd BiThe emission current of the light source (c),
Figure BDA00031980191200000310
is and
Figure BDA00031980191200000311
the corresponding observed voltage.
Optionally, after the detection is finished, taking an error record in the actually measured data as a criterion for the abnormal interpretation;
when the monopole receiving electrode M is fixedly arranged on the heading face, the measured data is interpreted as abnormal if the following formula (3a) is satisfied:
Figure BDA00031980191200000312
when the dipole reception electrodes M and N are fixedly arranged on the heading face, the measured data is interpreted as abnormal if the following equation (3b) is satisfied:
Figure BDA00031980191200000313
in the formula, Mean ± s.d. represents Mean ± standard deviation.
Compared with the prior art, the invention has the beneficial technical effects that:
(1) according to the property of the near-range effect of the direct-current electric field, the receiving electrode is fixedly arranged on the tunneling working face, and response signals of an abnormal body in front of the tunneling face under emission excitation are approached to the maximum extent; the emission electrode moves towards the direction of the heading face, the depth of the current field penetrating into the front of the heading face is observed from shallow to deep, and the signal-to-noise ratio is increased along with the increase of the detection distance.
(2) The maximum detection distance estimation formula provided by the invention establishes the relationship among various elements such as the geometric space of electrode arrangement, the emission current, the rock formation resistivity, the environmental noise level and the like, and provides a quantitative basis for determining the construction parameters of the underground direct current axial dipole dynamic source advanced detection.
(3) The sounding point and apparent resistivity formula corresponding to the movement of the transmitting electrode once is an algorithm designed for the invention, and has specificity; and the error record of the actually measured data is used as a discrimination standard of abnormal response, so that more guarantee is provided for the reliability of data interpretation.
(4) The axial dipole transmission has lower requirement on construction conditions than monopole transmission, can be arranged and constructed without a long roadway, is light and portable in required equipment, and has strong adaptability to the environment and high construction efficiency. The monopole or dipole receiving mode of the device can be selected according to different working conditions and different technical problems.
Drawings
Fig. 1 is a schematic view of the construction layout of the present invention, wherein 1a is monopole reception and 1b is dipole reception.
Fig. 2 is a schematic diagram of the advanced sounding point of the present invention, in which 2a is monopole reception and 2b is dipole reception.
Fig. 3 is a schematic diagram of equipotential between a roadway and an abnormal body under the excitation of an axial dipole source, wherein 3a is a low-resistance body in front of a driving surface, 3b is a high-resistance body in front of the driving surface, and a curve in the diagram is an equipotential line.
Fig. 4 is a graph of experimental calculation results.
In the figure: 1-tunneling surface; 2-a roadway; 3-an anomaly; 4-earth; 5-low resistance body; 6-high resistance.
The invention is described in detail below with reference to the drawings and the detailed description.
Detailed Description
In order to further improve the accuracy of the underground direct current advanced detection, the axial dipole dynamic source advanced detection method of the invention has two forms of axial dipole transmitting-monopole receiving and axial dipole transmitting-dipole receiving: the receiving electrode M or the receiving electrodes M and N are fixed on the excavation face, the transmitting electrodes A and B are arranged along the center line of the roadway floor, and the distance from the excavation face A1B1A series of moving sources A are formed by moving along the center line of the bottom plate and at intervals required by resolution ratio towards the heading faceiBi. The invention provides a maximum detection distance estimation formula of two receiving modes of a dipole dynamic source and an axial source, and AiBiCorresponding depth measuring points and apparent resistivity calculation formulas, an alternative method of electrode arrangement when construction conditions are limited, an acquisition method of noise level in a maximum detection distance estimation formula, a measurement method of signal-to-noise ratio in a detection process, a formula for judging abnormity by using error records in measured data and the like.
Following the above technical solutions, specific embodiments of the present invention are provided below, and it should be noted that the present invention is not limited to the following specific embodiments, and all equivalent changes based on the technical solutions of the present application fall within the protection scope of the present invention. The present invention will be described in further detail with reference to examples.
Example 1:
the embodiment provides an advanced detection method of an underground direct current axial dipole dynamic source, which is respectively shown in a figure 1a and a figure 1b, and the method is characterized in that a receiving electrode M or receiving electrodes M and N are fixedly arranged on a tunneling surface and are positioned on a vertical bisector of the bottom edge of the tunneling surface, so that geological abnormal bodies in front of the tunneling working surface are approached as much as possible, and the influence of the abnormal bodies near a top bottom plate and a side wall of a roadway is avoided; emitter electrodes A and B from A1And B1The position of the probe moves to the direction of the heading face at intervals required by resolution along the center line of the roadway floor, the signal-to-noise ratio increases along with the increase of the detection distance, and the probe is sequentially positioned at A1And B1Position point, A2And B2Location points, aiAnd BiPoints of position, anAnd BnA dipole moving source of location points; transmitting current once by the transmitting electrodes A and B at each position point, measuring voltage once by the receiving electrode M or the receiving electrodes M and N correspondingly, and carrying out dipole dynamic source advanced detection; when the influence of a roadway is neglected, the initial position of the transmitting electrode and the position of the maximum advanced detection distance are symmetrical relative to the tunneling surface, so that the distance between the initial position of the transmitting electrode and the midpoint O of the bottom edge of the tunneling surface is the maximum detection distance;
when the monopole receiving electrode M is fixedly arranged on the heading face, the estimation formula of the maximum advanced detection distance is as follows:
Figure BDA0003198019120000051
in the above formula, O1' O ' is emission electrode distance AB midpoint O '1Distance to the middle point O of the bottom edge of the driving face, IABmaxIn order to achieve the maximum emission current,
Figure BDA0003198019120000052
the rho is the resistivity of the rock stratum in front of the driving face and is the noise level observed by the receiving electrode M;
when the dipole receiving electrodes M and N are fixedly arranged on the tunneling surface, the estimation formula of the maximum detection distance of the advanced detection is as follows:
Figure BDA0003198019120000053
in the above formula, O'1O is the midpoint of emitter electrode distance AB'1The distance from the midpoint O of the bottom edge of the heading face to the point O, O 'O is the distance from the receiving electrode MN to the point O from the midpoint O', IABmaxIn order to achieve the maximum emission current,
Figure BDA0003198019120000054
ρ is the resistivity of the formation in front of the face for the noise level observed by the receiver electrodes M and N.
A1And B1Position point, A2And B2Location points, aiAnd BiPoints of position, anAnd BnThe smaller the interval between the position points is, the higher the detection resolution is; emitter electrodes A and B from A1And B1The part moves to the direction of the heading face at intervals required by resolution along the center line of the roadway bottom plate, and the depth of the current field penetrating into the front of the heading face is from shallow to deep.
In particular, noise levels are determined by the receiver null acquisition when the transmitter is not transmitting current before sounding is initiated
Figure BDA0003198019120000061
Or
Figure BDA0003198019120000062
The maximum detection distance is estimated by substituting into equations (1a) and (1 b).
Each transmitting electrode A in the detecting processiIn position, the noise level is determined by the empty and actual acquisition of the receiver when the transmitter is not transmitting and when it is transmitting current
Figure BDA0003198019120000063
Or
Figure BDA0003198019120000064
And receiving a voltage
Figure BDA0003198019120000065
Or
Figure BDA0003198019120000066
So as to decide whether to make repeated observations or not based on the signal-to-noise ratio. As the transmitting-receiving distance is reduced in the process that the transmitting electrode moves towards the tunneling surface, the signal-to-noise ratio is increased along with the increase of the detection distance. This step may be omitted when the downhole operation time is limited.
Calculating the emitter electrode A from equation (2)iAnd judging whether the abnormal body in front of the driving face is a low resistance body (figure 3a) or a high resistance body (figure 3b) according to the depth measuring point and the apparent resistivity.
When the monopole receiving electrode M is fixedly arranged on the heading face, the monopole receiving electrode M is connected with the AiBiCorresponding depth measuring point DiAnd apparent resistivity
Figure BDA00031980191200000614
The formula of (1) is:
Di≈O′1O-O′iO，(i＝1,2,...,n) (2a)
o 'of the above formula'1O is determined by formula (1a), O'iO is the midpoint of emitter electrode distance AB'iDistance to O;
Figure BDA0003198019120000067
in the above formula
Figure BDA0003198019120000068
Is AiThe distance to M,
Figure BDA0003198019120000069
Is BiA distance to M, wherein AiO is AiDistance to O, BiO is BiDistance to O, MO is the distance from M to O,
Figure BDA00031980191200000610
is AiBiThe emission current of (a) is measured,
Figure BDA00031980191200000611
is and
Figure BDA00031980191200000612
a corresponding observed voltage;
when the dipole reception electrodes M and N are fixedly arranged on the heading face, the dipole reception electrodes A and M are connected with the dipole reception electrodes AiBiCorresponding depth measuring point DiAnd apparent resistivity
Figure BDA00031980191200000613
The formula of (1) is:
Di≈O′1O-O′iO，(i＝1,2,...,n) (2c)
of which O'1O is determined by formula (1b), O'iO is the midpoint of emitter electrode distance AB'iDistance to O;
Figure BDA0003198019120000071
in the formula
Figure BDA0003198019120000072
AiThe distance to M,
Figure BDA0003198019120000073
Is BiThe distance to M,
Figure BDA0003198019120000074
AiThe distance to N,
Figure BDA0003198019120000075
Is BiA distance to N, wherein AiO is AiDistance to O, BiO is BiDistance to O, MO is the distance from M to O, NO is the distance from N to OThe distance of the oxygen is greater than the distance of the oxygen,
Figure BDA0003198019120000076
is AiBiThe emission current of (a) is measured,
Figure BDA0003198019120000077
is and
Figure BDA0003198019120000078
the corresponding observed voltage.
If limited by construction conditions, the transmitting electrode A can also be arranged at other positions of a roadway bottom plate or on a top plate and a side wall of the roadway and moves towards the heading face to form a series of dynamic sources A1B1，A2B2，...，AiBi，...，AnBnThe equations (1a), (1b) can still be used to estimate the maximum detection distance, and the equations (2a), (2b), (2c) and (2d) can still be used to calculate the sounding point and apparent resistivity corresponding to each movement of the transmitting electrode.
If the receiving electrodes are limited by construction conditions, the receiving electrodes M and the receiving electrodes MN can also be fixedly arranged on a roadway top plate, a bottom plate or a side wall close to the tunneling surface; the formulas (1a), (1b) can still be used for estimating the maximum detection distance, and the formulas (2a), (2b), (2c) and (2d) can still be used for calculating the sounding point and the apparent resistivity corresponding to each moving time of the transmitting electrode.
After the detection is finished, taking error records in the actually measured data as a discrimination standard of the abnormal interpretation;
when the monopole receiving electrode M is fixedly arranged on the heading face, the measured data is interpreted as abnormal if the following formula (3a) is satisfied:
Figure BDA0003198019120000079
when the dipole reception electrodes M and N are fixedly arranged on the heading face, the measured data is interpreted as abnormal if the following equation (3b) is satisfied:
Figure BDA00031980191200000710
in the formula, Mean ± s.d. represents Mean ± standard deviation.
And (3) verification experiment:
in order to verify the correctness and the validity of the formula, numerical simulation calculation is carried out by using the formula. Assuming that 1 high-resistance body with the radius of 5m is arranged at the position 20m below the roadway floor, the resistivity of the high-resistance body is 1000 omega.m, the resistivity of the surrounding rock is 100 omega.m, and the calculation is carried out by adopting an axial dipole transmitting and dipole MN receiving device and utilizing a formula (2 d). Fig. 4 shows the corresponding calculation results, with distance in m on the abscissa and apparent resistivity in Ω · m on the ordinate of fig. 4. It can be seen from the figure that the position of the maximum value of the anomaly calculated by the formula corresponds to the actual model position, indicating that the formula is valid and usable.

Claims (4)

1. A method for detecting the moving source of underground DC axial dipole is characterized by that the receiving electrode M or receiving electrodes M and N are fixed on the heading face and positioned on the vertical bisector of the bottom edge of heading face, and the emitting electrodes A and B are arranged from A1And B1The part moves along the central line of the roadway floor to the direction of the driving surface to form a position A1And B1Position point, A2And B2Location points, aiAnd BiPoints of position, anAnd BnA dipole moving source of location points; transmitting current once by the transmitting electrodes A and B at each position point, and measuring voltage once by the receiving electrode M or the receiving electrodes M and N correspondingly, and performing dipole dynamic source advanced detection;
when the monopole receiving electrode M is fixedly arranged on the heading face, the estimation formula of the maximum advanced detection distance is as follows:
Figure FDA0003198019110000011
in the above formula, O'1O is the midpoint of emitter electrode distance AB'1Distance to the middle point O of the bottom edge of the driving face, IABmaxIs at mostThe current is emitted and the current is measured,
Figure FDA0003198019110000012
the rho is the resistivity of the rock stratum in front of the driving face and is the noise level observed by the receiving electrode M;
when the dipole receiving electrodes M and N are fixedly arranged on the tunneling surface, the estimation formula of the maximum detection distance of the advanced detection is as follows:
Figure FDA0003198019110000013
in the above formula, O'1O is the midpoint of emitter electrode distance AB'1The distance from the midpoint O of the bottom edge of the heading face to the point O, O 'O is the distance from the receiving electrode to the point O from the midpoint O' of the MN, IABmaxIn order to achieve the maximum emission current,
Figure FDA0003198019110000014
ρ is the resistivity of the formation in front of the face for the noise level observed by the receiver electrodes M and N.
2. The method of claim 1, wherein A is a direct current axial dipole source advanced detection method1And B1Position point, A2And B2Location points, aiAnd BiPoints of position, anAnd BnThe smaller the interval between the position points is, the higher the detection resolution is; emitter electrodes A and B from A1And B1The part moves to the direction of the heading face at intervals required by resolution along the center line of the roadway bottom plate, and the depth of the current field penetrating into the front of the heading face is from shallow to deep.
3. The method of claim 1 wherein the monopole receiving electrode M is fixedly disposed at the face of the borehole and is aligned with the transmitting electrodes A and B at AiAnd BiAt the corresponding sounding point DiAnd apparent resistivity
Figure FDA0003198019110000021
The formula of (1) is:
Di≈O′1O-O′iO，(i＝1,2,...,n) (2a)
o 'of the above formula'1O is determined by formula (1a), O'iO is the midpoint of emitter electrode distance AB'iDistance to O;
Figure FDA0003198019110000022
in the above formula
Figure FDA0003198019110000023
Is AiThe distance to M,
Figure FDA0003198019110000024
Is BiA distance to M, wherein AiO is AiDistance to O, BiO is BiDistance to O, MO is the distance from M to O,
Figure FDA00031980191100000216
the emitting electrodes A and B are at AiAnd BiThe emission current of the light source (c),
Figure FDA0003198019110000025
is and
Figure FDA00031980191100000213
a corresponding observed voltage;
when the dipole receiving electrodes M and N are fixedly arranged on the heading face, the dipole receiving electrodes M and N are in A connection with the transmitting electrodes A and BiAnd BiAt the corresponding sounding point DiAnd apparent resistivity
Figure FDA0003198019110000026
The formula of (1) is:
Di≈O′1O-O′iO，(i＝1,2,...,n) (2c)
of which O'1O is determined by formula (1b), O'iO is the midpoint of emitter electrode distance AB'iDistance to O;
Figure FDA0003198019110000027
in the formula
Figure FDA0003198019110000028
AiThe distance to M,
Figure FDA0003198019110000029
Is BiThe distance to M,
Figure FDA00031980191100000210
AiThe distance to N,
Figure FDA00031980191100000211
Is BiA distance to N, wherein AiO is AiDistance to O, BiO is BiDistance to O, MO is the distance from M to O, NO is the distance from N to O,
Figure FDA00031980191100000214
the emitting electrodes A and B are at AiAnd BiThe emission current of the light source (c),
Figure FDA00031980191100000212
is and
Figure FDA00031980191100000215
the corresponding observed voltage.
4. The method according to claim 3, wherein after the detection is finished, the error record in the measured data is used as the criterion for the abnormal interpretation;
when the monopole receiving electrode M is fixedly arranged on the heading face, the measured data is interpreted as abnormal if the following formula (3a) is satisfied:
Figure FDA0003198019110000031
when the dipole reception electrodes M and N are fixedly arranged on the heading face, the measured data is interpreted as abnormal if the following equation (3b) is satisfied:
Figure FDA0003198019110000032
in the formula, Mean ± s.d. represents Mean ± standard deviation.
CN202110896200.XA 2021-08-05 2021-08-05 Advanced detection method for underground direct-current axial dipole-motion source Active CN113885085B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110896200.XA CN113885085B (en) 2021-08-05 2021-08-05 Advanced detection method for underground direct-current axial dipole-motion source

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110896200.XA CN113885085B (en) 2021-08-05 2021-08-05 Advanced detection method for underground direct-current axial dipole-motion source

Publications (2)

Publication Number Publication Date
CN113885085A true CN113885085A (en) 2022-01-04
CN113885085B CN113885085B (en) 2023-06-09

Family

ID=79010927

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110896200.XA Active CN113885085B (en) 2021-08-05 2021-08-05 Advanced detection method for underground direct-current axial dipole-motion source

Country Status (1)

Country Link
CN (1) CN113885085B (en)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000017489A2 (en) * 1998-09-19 2000-03-30 Armin Kaus Geoelectric pre-prospecting method
FR2855272A1 (en) * 2003-05-22 2004-11-26 Schlumberger Services Petrol Characterizing underground formation comprises use of diagraph at bottom of borehole which is rotated about azimuth to allow directional measurements to be taken
DE102007021399A1 (en) * 2007-05-18 2008-11-20 Kopp, Thomas Geoelectrical parameter determining method for measuring device for e.g. gas exploration, involves utilizing shield and measuring electrodes for frequency-dependent determination based on analysis of measuring circuit resistances
CN101603423A (en) * 2009-07-09 2009-12-16 煤炭科学研究总院西安研究院 A kind of in coal mine roadway the DC electrical method method of bed-parallel advanced detection of water bearing
CN103278857A (en) * 2013-05-13 2013-09-04 江苏大学 Method for designing underground direct current advanced detection device
CN104007472A (en) * 2014-06-11 2014-08-27 中煤科工集团西安研究院有限公司 In-hole direct current electric method advanced detection method
CN104714254A (en) * 2015-02-02 2015-06-17 中国科学院地质与地球物理研究所 Coal mine water-accumulating gob detection method and device based on transient electromagnetic method
CN106907145A (en) * 2017-02-09 2017-06-30 武汉地大华睿地学技术有限公司 A kind of apparent resistivity measuring system and method with brill advanced prediction
CN107861159A (en) * 2017-10-19 2018-03-30 中煤科工集团西安研究院有限公司 Double Electric Dipole ground well transient electromagnetic detecting methods
US20190339413A1 (en) * 2016-12-14 2019-11-07 China University Of Mining And Technology Mine tem three-component detection method
AU2020100984A4 (en) * 2020-06-11 2020-08-13 Xi’an Northwest Nonferrous Geophysical & Geochemical Exploration Co., Ltd A Method and Apparatus for Ground-tunnel Wide Field Electromagnetic Surveying

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000017489A2 (en) * 1998-09-19 2000-03-30 Armin Kaus Geoelectric pre-prospecting method
FR2855272A1 (en) * 2003-05-22 2004-11-26 Schlumberger Services Petrol Characterizing underground formation comprises use of diagraph at bottom of borehole which is rotated about azimuth to allow directional measurements to be taken
DE102007021399A1 (en) * 2007-05-18 2008-11-20 Kopp, Thomas Geoelectrical parameter determining method for measuring device for e.g. gas exploration, involves utilizing shield and measuring electrodes for frequency-dependent determination based on analysis of measuring circuit resistances
CN101603423A (en) * 2009-07-09 2009-12-16 煤炭科学研究总院西安研究院 A kind of in coal mine roadway the DC electrical method method of bed-parallel advanced detection of water bearing
CN103278857A (en) * 2013-05-13 2013-09-04 江苏大学 Method for designing underground direct current advanced detection device
CN104007472A (en) * 2014-06-11 2014-08-27 中煤科工集团西安研究院有限公司 In-hole direct current electric method advanced detection method
CN104714254A (en) * 2015-02-02 2015-06-17 中国科学院地质与地球物理研究所 Coal mine water-accumulating gob detection method and device based on transient electromagnetic method
US20190339413A1 (en) * 2016-12-14 2019-11-07 China University Of Mining And Technology Mine tem three-component detection method
CN106907145A (en) * 2017-02-09 2017-06-30 武汉地大华睿地学技术有限公司 A kind of apparent resistivity measuring system and method with brill advanced prediction
CN107861159A (en) * 2017-10-19 2018-03-30 中煤科工集团西安研究院有限公司 Double Electric Dipole ground well transient electromagnetic detecting methods
AU2020100984A4 (en) * 2020-06-11 2020-08-13 Xi’an Northwest Nonferrous Geophysical & Geochemical Exploration Co., Ltd A Method and Apparatus for Ground-tunnel Wide Field Electromagnetic Surveying

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张平松;李永盛;胡雄武;: "巷道掘进直流电阻率法超前探测技术应用探讨" *
梁盛军;刘会臣;谭捍东;: "典型地电模型矿井直流电法异常规律研究" *

Also Published As

Publication number Publication date
CN113885085B (en) 2023-06-09

Similar Documents

Publication Publication Date Title
CN106772644B (en) mine transient electromagnetic three-component detection method
WO2020078003A1 (en) Time-domain transient electromagnetic wave well logging far-boundary detection method
CN103675922B (en) Operation period underground pipeline pipe diameter measuring method based on ground penetrating radar
CN102156301B (en) Advanced-prediction observation system while drilling
Ungureanu et al. Use of electric resistivity tomography (ERT) for detecting underground voids on highly anthropized urban construction sites
US20130166215A1 (en) System and method of making environmental measurements
CN109143378B (en) Secondary time difference method for bedding advanced detection of water-containing structure in coal mine tunnel
CN106094043A (en) Transient electromagnetic method ground hole detection method and device
CN104360401B (en) A kind of transient electromagnetic B field descends objective body geological information method definitely
CN101263404A (en) High resolution resistivity earth imager
US20120130641A1 (en) Marine Source To Borehole Electromagnetic Mapping Of Sub-Bottom Electrical Resistivity
NO335899B1 (en) Determination of anisotropy in subsurface formations around a borehole with dipole moment angled with respect to the longitudinal axis of the logging tool
US9810805B2 (en) Method and apparatus to detect a conductive body
CN109343130B (en) Laterally-excited loop source ground well transient electromagnetic detection method and system
CN105116453A (en) Transient electromagnetic exploration method and device for natural gas hydrate in permafrost zone
US10317560B2 (en) Systems and methods of robust determination of boundaries
CN103278855B (en) Method for eliminating influence of roadways and terrains on apparent resistivity in direct-current exploration
CN108828676A (en) A kind of ground-mine laneway transient electromagnetic three-component detection method
CN105842740A (en) Fixed point rotary irradiation large power transient electromagnetic method
CN103821495A (en) Logging method
CN113885083B (en) Underground direct-current axial dipole-motion source abnormity self-display advanced detection method
CN113885085A (en) Underground direct-current axial dipole dynamic source advanced detection method
CN113885086B (en) Underground direct-current equatorial direction dipole dynamic source abnormity self-explicit advanced detection method
CN113900151B (en) Underground direct-current monopole motion source abnormity self-display advanced detection method
CN113885084B (en) Advanced detection method for underground direct-current monopole moving source

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant