MULTI-LINE SOURCE GROUND-BOREHOLE TRANSIENT ELECTROMAGNETIC DETECTION METHOD AND MULTI-LINE SOURCE GROUND-BOREHOLE TRANSIENT ELECTROMAGNETIC DETECTION

Abstract

Disclosed are a multi-line source ground-borehole transient electromagnetic detection method and a multi-line source ground-borehole transient electromagnetic detection device. The method includes following steps: S1, constructing a multi-line source ground-borehole transient electromagnetic forward model; S2, obtaining multi-line source ground-borehole transient electromagnetic responses of an underground target layer by the multi-line source ground-borehole transient electromagnetic forward model; and S3, recognizing the underground target layer according to electromagnetic diffusion characteristics of the multi-line source ground-borehole transient electromagnetic responses. A resolution and a detection capability of the underground target by a ground-borehole transient electromagnetic method is greatly improved by adopting technical schemes.

Claims

1. A multi-line source ground-borehole transient electromagnetic detection method, comprising: S1, constructing a multi-line source ground-borehole transient electromagnetic forward model; S2, obtaining multi-line source ground-borehole transient electromagnetic responses of an underground target layer by the multi-line source ground-borehole transient electromagnetic forward model; and S3, recognizing the underground target layer according to electromagnetic diffusion characteristics of the multi-line source ground-borehole transient electromagnetic responses.

2. The multi-line source ground-borehole transient electromagnetic detection method according to claim 1, wherein the S1 is specifically as follows: when an electromagnetic field is excited in a layered earth surface and received in a borehole, an electromagnetic field component excited by a single grounding conductor source in an x-direction is expressed as follows: $\begin{array}{cc}{E}_x=-\frac{{\hat{z}}_{0}I}{4\pi }{\int }_{-L}^L{\int }_0^{\infty }\left(1+r_{\mathrm{TE}}\right)\frac{\lambda }{u_0}{e}^{u_{0}z}{J}_0\left(\lambda r\right)d\lambda {\mathrm{dx}}^{\prime }-{\frac{I}{4\pi }\left[\frac{x}{r}{\int }_0^{\infty }\left[\left(1-r_{\mathrm{TM}}\right)\frac{u_0}{{\hat{y}}_0}-\left(1+r_{\mathrm{TE}}\right)\frac{{\hat{z}}_0}{u_0}\right]e^{u_{0}z}{J}_1\left(\lambda r\right)d\lambda \right]}_{R_1}^{R_2}& \left(1\right)\end{array}$ $E_y=-{\frac{I}{4\pi }\left[\frac{y}{r}{\int }_0^{\infty }\left[\left(1-r_{\mathrm{TM}}\right)\frac{u_0}{{\hat{y}}_0}-\left(1+r_{\mathrm{TE}}\right)\frac{{\hat{z}}_0}{u_0}\right]e^{u_{0}z}{J}_1\left(\lambda r\right)d\lambda \right]}_{R_1}^{R_2}$ $E_z=-{\frac{I}{4\pi }\left[{\int }_0^{\infty }\left(1-r_{\mathrm{TM}}\right)\frac{\lambda }{u_0}{e}^{u_{0}z}{J}_0\left(\lambda r\right)d\lambda \right]}_{R_1}^{R_2}$ $H_x=-\frac{I}{4\pi }\frac{y}{r}{\int }_0^{\infty }\left(r_{\mathrm{TM}}+r_{\mathrm{TE}}\right)e^{u_{0}z}{J}_1\left(\lambda r\right)d\lambda {.Math.}_{R_1}^{R_2}$ $H_y=-\frac{I}{4\pi }\frac{x}{r}{\int }_0^{\infty }\left(r_{\mathrm{TM}}+r_{\mathrm{TE}}\right)e^{u_{0}z}{J}_1\left(\lambda r\right)d\lambda {.Math.}_{R_1}^{R_2}-\frac{I}{4\pi }{\int }_{-L}^L{\int }_0^{\infty }\left(1-r_{\mathrm{TE}}\right)e^{u_{0}z}\lambda J_0\left(\lambda r\right)d\lambda {\mathrm{dx}}^{\prime }$ $H_z=\frac{I}{4\pi }{\int }_{-L}^L\frac{y}{r}{\int }_0^{\infty }\left(1+r_{\mathrm{TE}}\right)e^{u_{0}z}\frac{{\lambda }^2}{u_0}{J}_1\left(\lambda r\right)d\lambda {\mathrm{dx}}^{\prime },$ wherein J.sub.1(λr) and J.sub.0(λr) are first-order and zero-order Bessel functions, respectively, λ represents a horizontal wave number, r.sub.i represents receiving and sending distances of measuring points in coordinate systems to which different line sources belong; h.sub.j represents a thickness of a j-th layer; {circumflex over (z)}.sub.j represents an impedance of the j-th layer, ŷ.sub.j represents an admittance of the j-th layer, and z represents a height of a receiving point; $r_{\mathrm{TE}}=\frac{Y_0-{\hat{Y}}_1}{Y_0+{\hat{Y}}_1}$  indicates a reflection coefficient of an electromagnetic wave when an electric field component is incident on a plane, $r_{\mathrm{TM}}=\frac{Z_0-{\stackrel{⌢}{Z}}_1}{Z_0+{\stackrel{⌢}{Z}}_1}$  indicates the reflection coefficient of the electromagnetic wave when a vertical magnetic field is incident on the plane, $Y_0=\frac{u_0}{{\hat{z}}_0}$  indicates an intrinsic admittance of a free space, Ŷ.sub.1 indicates a surface admittance, $Z_0=\frac{u_0}{{\hat{y}}_0}$  indicates an intrinsic impedance of the free space, and {circumflex over (Z)}.sub.1 indicates a surface impedance, where {circumflex over (z)}.sub.0=iωμ.sub.0 and ŷ.sub.0=iωε.sub.0; ω represents an angular frequency, i represents a complex number, ε.sub.0 represents a dielectric coefficient of an underground uniform half space, and μ.sub.0 represents a magnetic permeability of the underground uniform half space; for a multi-line source ground-borehole transient electromagnetic field, z is positive; z is equal to 0 when observed on the ground, but negative when observed in the air; a solution of the responses is obtained by a superposition calculation of the responses of various line sources; Ex, EY, Ez, Hx, HY, and Hz represent electric field components and magnetic field components in an x direction, a y direction, and a z direction respectively.

3. The multi-line source ground-borehole transient electromagnetic detection method according to claim 1, wherein in the S2, the multi-line source ground-borehole transient electromagnetic responses in the underground target layer are obtained by designing following parameters: eight line sources are set, an included angle between adjacent line sources is 45 degrees, a length of each line source is set to 500 meters, an emission current is 10 A, and depths of receiving points are 100, 200, 300, 400 and 500 meters respectively.

4. The multi-line source ground-borehole transient electromagnetic detection method according to claim 2, wherein in the S1, a coordinate system is established with each line source as an x axis and a midpoint of each line source as a source point; for the same measuring point P(x.sub.i,y.sub.1), coordinates in different coordinate systems are obtained by relationships between translations and rotations of the coordinate systems; (a) when the coordinate system X.sub.1OY.sub.1 is translated to X′O′Y′, the coordinates are:
x′=x.sub.1+a
y′=y.sub.1±b, wherein a=L−L cos θ.sub.i, b=L sin θ.sub.i; (b) when the coordinate system X′O′Y′ is rotated to the coordinate system X.sub.iO′Y.sub.i to which the i-th source belongs:
x.sub.i=x′ cos θ.sub.i−y′ sin θ.sub.i
y.sub.i=y′ cos θ.sub.i+x′ sin θ.sub.i, each component of the magnetic field generated by a multi-line source excitation is obtained by the superposition.

5. The multi-line source ground-borehole transient electromagnetic detection method as claimed in claim 1, wherein in the S3, a resolution and a recognition capability of an underground high resistance layer are judged by analyzing response diffusion properties of the underground target layer.

6. A multi-line source ground-borehole transient electromagnetic detection device, comprising: a construction module for constructing a multi-line source ground-borehole transient electromagnetic forward model; an analysis module for obtaining multi-line source ground-borehole transient electromagnetic responses of an underground target layer through the multi-line source ground-borehole transient electromagnetic forward model; and a recognition module for recognizing the underground target layer according to electromagnetic diffusion characteristics of the multi-line source ground-borehole transient electromagnetic responses.

7. The multi-line source ground-borehole transient electromagnetic detection device according to claim 6, wherein a specific process of the construction module comprises: when an electromagnetic field is excited in a layered earth surface and received in a borehole, an electromagnetic field component excited by a single grounding conductor source in an x-direction is expressed as follows: $\begin{array}{cc}{E}_x=-\frac{{\hat{z}}_{0}I}{4\pi }{\int }_{-L}^L{\int }_0^{\infty }\left(1+r_{\mathrm{TE}}\right)\frac{\lambda }{u_0}{e}^{u_{0}z}{J}_0\left(\lambda r\right)d\lambda {\mathrm{dx}}^{\prime }-{\frac{i}{4\pi }\left[\frac{x}{r}{\int }_0^{\infty }\left[\left(1-r_{\mathrm{TM}}\right)\frac{u_0}{{\hat{y}}_0}-\left(1+r_{\mathrm{TE}}\right)\frac{{\hat{z}}_0}{u_0}\right]e^{u_{0}z}{J}_1\left(\lambda r\right)d\lambda \right]}_{R_1}^{R_2}& \left(1\right)\end{array}$ $E_y=-{\frac{I}{4\pi }\left[\frac{y}{r}{\int }_0^{\infty }\left[\left(1-r_{\mathrm{TM}}\right)\frac{u_0}{{\hat{y}}_0}-\left(1+r_{\mathrm{TE}}\right)\frac{{\hat{z}}_0}{u_0}\right]e^{u_{0}z}{J}_1\left(\lambda r\right)d\lambda \right]}_{R_1}^{R_2}$ $E_z=-{\frac{I}{4\pi }\left[{\int }_0^{\infty }\left(1-r_{\mathrm{TM}}\right)\frac{\lambda }{u_0}{e}^{u_{0}z}{J}_0\left(\lambda r\right)d\lambda \right]}_{R_1}^{R_2}$ $H_x=-\frac{I}{4\pi }\frac{y}{r}{\int }_0^{\infty }\left(r_{\mathrm{TM}}+r_{\mathrm{TE}}\right)e^{u_{0}z}{J}_1\left(\lambda r\right)d\lambda {.Math.}_{R_1}^{R_2}$ $H_y=-\frac{I}{4\pi }\frac{x}{r}{\int }_0^{\infty }\left(r_{\mathrm{TM}}+r_{\mathrm{TE}}\right)e^{u_{0}z}{J}_1\left(\lambda r\right)d\lambda {.Math.}_{R_1}^{R_2}-\frac{I}{4\pi }{\int }_{-L}^L{\int }_0^{\infty }\left(1-r_{\mathrm{TE}}\right)e^{u_{0}z}\lambda J_0\left(\lambda r\right)d\lambda {\mathrm{dx}}^{\prime }$ $H_z=\frac{I}{4\pi }{\int }_{-L}^L\frac{y}{r}{\int }_0^{\infty }\left(1+r_{\mathrm{TE}}\right)e^{u_{0}z}\frac{{\lambda }^2}{u_0}{J}_1\left(\lambda r\right)d\lambda {\mathrm{dx}}^{\prime },$ wherein J.sub.1(λr) and J.sub.0(λr) are first-order and zero-order Bessel functions, respectively, λ represents a horizontal wave number, r.sub.i represents receiving and sending distances of measuring points in coordinate systems to which different line sources belong; h.sub.j represents a thickness of a j-th layer; {circumflex over (z)}.sub.j represents an impedance of the j-th layer, ŷ.sub.j represents an admittance of the j-th layer, and z represents a height of a receiving point; $r_{\mathrm{TE}}=\frac{Y_0-{\hat{Y}}_1}{Y_0+{\hat{Y}}_1}$  indicates a reflection coefficient of an electromagnetic wave when an electric field component is incident on a plane, $r_{\mathrm{TM}}=\frac{Z_0-{\stackrel{⌢}{Z}}_1}{Z_0+{\stackrel{⌢}{Z}}_1}$  indicates the reflection coefficient of the electromagnetic wave when a vertical magnetic field is incident on the plane, $Y_0=\frac{u_0}{{\hat{z}}_0}$  indicates an intrinsic admittance of a free space, Ŷ.sub.1 indicates a surface admittance, $Z_0=\frac{u_0}{{\hat{y}}_0}$  indicates an intrinsic impedance of the free space, and {circumflex over (Z)}.sub.1 indicates a surface impedance, where {circumflex over (z)}.sub.0=iωμ.sub.0 and ŷ.sub.0=iωε.sub.0; ω represents an angular frequency, i represents a complex number, ε.sub.0 represents a dielectric coefficient of an underground uniform half space, and μ.sub.0 represents a magnetic permeability of the underground uniform half space; for a multi-line source ground-borehole transient electromagnetic field, z is positive; z is equal to 0 when observed on the ground, but negative when observed in the air; a solution of the responses is obtained by a superposition calculation of the responses of various line sources; Ex, EY, Ez, Hx, HY, and Hz represent electric field components and magnetic field components in an x direction, a y direction, and a z direction respectively.

8. The multi-line source ground-borehole transient electromagnetic detection device according to claim 6, wherein the analysis module obtains the multi-line source ground-borehole transient electromagnetic responses in the underground target layer by designing following parameters: eight line sources are set, an included angle between adjacent line sources is 45 degrees, a length of each line source is set to 500 meters, an emission current is 10 A, and depths of receiving points are 100, 200, 300, 400 and 500 meters respectively.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0031] Here some specific embodiments of the present application are described in detail by way of example and not limitation with reference to the accompanying drawings. Same reference numerals in the drawings indicate the same or similar parts. Those skilled in the art should understand that these drawings are not necessarily drawn to scale.

[0032] FIG. 1 is a structural diagram of an emission source device.

[0033] FIG. 2 is a flowchart of a method according to the application.

[0034] FIG. 3 is a schematic diagram of a coordinate transformation relationship between different sources.

[0035] FIGS. 4a-4f are schematic diagrams of response plane distributions (1 ms) at 200 m depth, in which,

[0036] FIG. 4a is a schematic diagram of a horizontal electric field component distribution in an x direction;

[0037] FIG. 4b is a schematic diagram of the horizontal electric field component distribution in a y direction;

[0038] FIG. 4c is a schematic diagram of a vertical electric field component distribution in a z direction;

[0039] FIG. 4d is a schematic diagram of a horizontal magnetic field component distribution in the x direction;

[0040] FIG. 4e is a schematic diagram of the horizontal magnetic field component distribution in the y direction; and

[0041] FIG. 4f is a schematic diagram of a vertical magnetic field component distribution in the z direction.

[0042] FIGS. 5a-5d are Ex response diffusion diagrams, in which,

[0043] FIG. 5a is a schematic diagram of the response diffusion at 1E-5 s;

[0044] FIG. 5b is a schematic diagram of the response diffusion at 1E-4 s;

[0045] FIG. 5c is a schematic diagram of the response diffusion at 1E-3 s; and

[0046] FIG. 5e is a schematic diagram of the response diffusion at 1E-2 s.

[0047] FIGS. 6a-6d are Ey response diffusion diagrams, in which,

[0048] FIG. 6a is a schematic diagram of the response diffusion at 1E-5 s;

[0049] FIG. 6b is a schematic diagram of the response diffusion at 1E-4 s;

[0050] FIG. 6c is a schematic diagram of the response diffusion at 1E-3 s; and

[0051] FIG. 6d is a schematic diagram of the response diffusion at 1E-2 s.

[0052] FIGS. 7a-7d are Ez response diffusion diagrams, in which,

[0053] FIG. 7a is a schematic diagram of the response diffusion at 1E-5 s;

[0054] FIG. 7b is a schematic diagram of the response diffusion at 1E-4 s;

[0055] FIG. 7c is a schematic diagram of the response diffusion at 1E-3 s; and

[0056] FIG. 7d is a schematic diagram of the response diffusion at 1E-2 s.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0057] In order to make a purpose, technical schemes and advantages of the application clearer, the application is described in detail below with reference to the attached drawings.

[0058] As shown in FIG. 1, the application provides a multi-line source ground-borehole transient electromagnetic detection method, which is based on a multi-line ground-borehole transient electromagnetic method (MGBTEM). A ground emitting source is composed of multiple groups of grounding line sources which emit currents in opposite directions, and may be combined with 8, 10 or 16 different sources. FIG. 2 shows an emission source device in a form of 10 line sources. Each emission source has a same length. Current directions diverge from a center to all directions, and included angles between adjacent emission sources are equal. Specifically, a multi-line source emitting device is arranged on the ground; the multi-line source emitting device consists of four pairs of grounding conductor sources with opposite current directions; through transmission lines, a grounding electrode and a magnetic sensor are put into the borehole, and a vertical electric field and magnetic field multi-components are observed in the borehole. One or more receiving magnetic probes and electrode pairs are arranged downward along the borehole, and the receiving magnetic probes and the electrode pairs are connected to a ground receiver. The ground receiver controls the receiving magnetic probes and the electrode pairs and collects received signals.

[0059] As shown in FIG. 1, the multi-line source ground-borehole transient electromagnetic detection method in an embodiment of the application includes following steps.

[0060] S1, constructing a multi-line source ground-borehole transient electromagnetic forward model.

[0061] when an electromagnetic field is excited in a layered earth surface and received in a borehole, an electromagnetic field component excited by a single grounding conductor source in an x-direction is expressed as follows:

[00006]$\begin{array}{cc}{E}_x=-\frac{{\hat{z}}_{0}I}{4\pi }{\int }_{-L}^L{\int }_0^{\infty }\left(1+r_{\mathrm{TE}}\right)\frac{\lambda }{u_0}{e}^{u_{0}z}{J}_0\left(\lambda r\right)d\lambda {\mathrm{dx}}^{\prime }-{\frac{I}{4\pi }\left[\frac{x}{r}{\int }_0^{\infty }\left[\left(1-r_{\mathrm{TM}}\right)\frac{u_0}{{\hat{y}}_0}-\left(1+r_{\mathrm{TE}}\right)\frac{{\hat{z}}_0}{u_0}\right]e^{u_{0}z}{J}_1\left(\lambda r\right)d\lambda \right]}_{R_1}^{R_2}& \left(1\right)\end{array}$$E_y=-{\frac{I}{4\pi }\left[\frac{y}{r}{\int }_0^{\infty }\left[\left(1-r_{\mathrm{TM}}\right)\frac{u_0}{{\hat{y}}_0}-\left(1+r_{\mathrm{TE}}\right)\frac{{\hat{z}}_0}{u_0}\right]e^{u_{0}z}{J}_1\left(\lambda r\right)d\lambda \right]}_{R_1}^{R_2}{E}_z=-{\frac{I}{4\pi }\left[{\int }_0^{\infty }\left(1-r_{\mathrm{TM}}\right)\frac{\lambda }{u_0}{e}^{u_{0}z}{J}_0\left(\lambda r\right)d\lambda \right]}_{R_1}^{R_2}{H}_x=-\frac{I}{4\pi }\frac{y}{r}{\int }_0^{\infty }\left(r_{\mathrm{TM}}+r_{\mathrm{TE}}\right)e^{u_{0}z}{J}_1\left(\lambda r\right)d\lambda {.Math.}_{R_1}^{R_2}$$H_y=-\frac{I}{4\pi }\frac{x}{r}{\int }_0^{\infty }\left(r_{\mathrm{TM}}+r_{\mathrm{TE}}\right)e^{u_{0}z}{J}_1\left(\lambda r\right)d\lambda {.Math.}_{R_1}^{R_2}-\frac{I}{4\pi }{\int }_{-L}^L{\int }_0^{\infty }\left(1-r_{\mathrm{TE}}\right)e^{u_{0}z}\lambda J_0\left(\lambda r\right)d\lambda {\mathrm{dx}}^{\prime }$$H_z=\frac{I}{4\pi }{\int }_{-L}^L\frac{y}{r}{\int }_0^{\infty }\left(1+r_{\mathrm{TE}}\right)e^{u_{0}z}\frac{{\lambda }^2}{u_0}{J}_1\left(\lambda r\right)d\lambda {\mathrm{dx}}^{\prime };$

[0062] where, J.sub.1(λr) and J.sub.0(λr) are first-order and zero-order Bessel functions, respectively, λ represents a horizontal wave number, r.sub.i represents receiving and sending distances of measuring points in coordinate systems to which different line sources belong; h.sub.j represents a thickness of a j-th layer; {circumflex over (z)}.sub.j represents an impedance of the j-th layer, ŷ.sub.j represents an admittance of the j-th layer, and z represents a height of a receiving point;

[00007]$r_{\mathrm{TE}}=\frac{Y_0-{\hat{Y}}_1}{Y_0+{\hat{Y}}_1}$

indicates a reflection coefficient of an electromagnetic wave when an electric field component is incident on a plane,

[00008]$r_{\mathrm{TM}}=\frac{Z_0-{\stackrel{⌢}{Z}}_1}{Z_0+{\stackrel{⌢}{Z}}_1}$

indicates the reflection coefficient of the electromagnetic wave when a vertical magnetic field is incident on the plane,

[00009]$Y_0=\frac{u_0}{{\hat{z}}_0}$

indicates an intrinsic admittance of a free space, Ŷ.sub.1 indicates a surface admittance,

[00010]$Z_0=\frac{u_0}{{\hat{y}}_0}$

indicates an intrinsic impedance of the free space, and {circumflex over (Z)}.sub.1 indicates a surface impedance, where {circumflex over (z)}.sub.0=iωμ.sub.0 and ŷ.sub.0=iωε.sub.0; ω represents an angular frequency, i represents a complex number, ε.sub.0 represents a dielectric coefficient of an underground uniform half space, and μ.sub.0 represents a magnetic permeability of the underground uniform half space.

[0063] For a multi-line source ground-borehole transient electromagnetic field, z is positive; z is equal to 0 when observed on the ground, but negative when observed in the air. A solution of responses is obtained by a superposition calculation of the responses of various line sources, and the coordinate system is established with each line source as an x axis and a midpoint of each line source as a source point, as shown in FIG. 3:

[0064] for the same measuring point P(x.sub.1,y.sub.1), coordinates in different coordinate systems are obtained by relationships between translations and rotations of the coordinate systems;

[0065] (e) when the coordinate system X.sub.1OY.sub.1 is translated to X′O′Y′, the coordinate is:

x′=x.sub.1+a

y′=y.sub.1+b,

[0066] where a=L−L cos θ.sub.i, b=L sin θ.sub.i;

[0067] (f) when the coordinate system X′O′Y′ is rotated to the coordinate system X.sub.iO′Y.sub.i to which the i-th source belongs:

x.sub.i=x′ cos θ.sub.i−y′ sin θ.sub.i

y.sub.i=y′ cos θ.sub.i+x′ sin θ.sub.i,

[0068] each component of the magnetic field generated by a multi-line source excitation is obtained by a superposition.

[0069] S2, obtaining multi-line source ground-borehole transient electromagnetic responses of an underground target by the multi-line source ground-borehole transient electromagnetic forward model.

[0070] Taking eight line sources as an example, the included angle between adjacent line sources is 45 degrees, the length of each line source is 500 meters, and an emission current is 10 A. Multi-line source ground-borehole transient electromagnetic field responses are analyzed by the multi-line source ground-borehole transient electromagnetic forward model (K model), and depths of the receiving points are 100, 200, 300, 400 and 500 meters respectively. Geoelectric parameters of the K model are as follows:

ρ.sub.1=300 Ωm, d.sub.1=200 m; ρ.sub.2=2000 Ωm, d.sub.2=200 m; ρ.sub.3=500 Ωm.

[0071] Taking a source center as a coordinate origin and the line source 1 as an X axis, a planar distribution of multi-components of the electromagnetic field at the depth of 200 m underground is obtained through a calculation, as shown in FIG. 4.

[0072] As shown in FIGS. 4a-4f, horizontal electric fields Ex and Ey show a typical axisymmetric distribution, and their response distribution characteristics are similar. For the Ex, the response of the Ex is 0 at the measuring point where x=0, and similarly, for the Ey, the Ey is 0 at the measuring point where y=0. Their response distribution characteristics should be paid a special attention to in an actual observation. In these special positions, the responses of the actual observations are not zero, but the responses of an underground three-dimensional target. Similar response distribution characteristics also exist in the distribution of horizontal magnetic fields. For the Ez, the response distribution presents a centrosymmetric distribution law, and the responses are centered on the source. With an increase of an offset, the responses gradually decrease. Therefore, if the actual observation is made, try to observe in an area near the source. Compared with a horizontal component, a vertical magnetic field response value is much lower, so the calculated response may be ignored. This is because the vertical magnetic field itself is generated by the line sources, and a combination of the multi-line sources counteracts the vertical magnetic field generated by the line sources. Therefore, in a case of the layered earth, the vertical magnetic field response in a forward calculation is ignored, but a response amplitude is not weak in the actual observation. The response amplitude is generated by the underground three-dimensional target. The observation of the vertical magnetic field realizes a real three-dimensional response observation, while ignoring a background field of the layered earth.

[0073] S3, recognizing an underground target layer according to electromagnetic diffusion characteristics of the multi-line source ground-borehole transient electromagnetic responses.

[0074] A greatest advantage of the multi-line sources is that a fine detection of high resistance targets is realized by observing electric field components, so in a ground-borehole transient electromagnetic observation, the application focuses on observing reflection characteristics of underground electric field components to a high resistance. In order to analyze a detection capability of the observed electric field component in the borehole to a high resistance target layer, response profiles of the high resistance target layer at different times are analyzed. The resolution and the recognition capability of an underground high resistance layer are judged by analyzing response diffusion properties of the underground target layer.

[0075] As shown in FIGS. 5a-5e, FIGS. 6a-6e and FIGS. 7a-7d, for the high resistance layer with the depth of 200-400 m, an Ez component shows an obvious response indication function, and a maximum isoline trap coincides with a top and a bottom of the high resistance layer, so that position information and thickness information of the high resistance layer are well reflected, indicating that the vertical electric field component of the transient electromagnetic field in the ground-borehole has a good resolution for the high resistance target.

[0076] A similar characteristic appears in an Ey component diffusion profile, but a difference is that this indicator characteristic only appears after a certain time, and only in a limited number of the measuring points and when an axial observation is satisfactory relative to the source. When the electric field is observed in the borehole, the vertical electric field should be selected as far as possible; when the electric field is not observed in a roadway, an axial electric field component and the vertical electric field component of the source are considered.

[0077] The application also provides a multi-line source ground-borehole transient electromagnetic detection device, including:

[0078] a construction module, which is used for constructing the multi-line source ground-borehole transient electromagnetic forward model;

[0079] an analysis module, which is used for obtaining the multi-line source ground-borehole transient electromagnetic responses of the underground target layer through the multi-line source ground-borehole transient electromagnetic forward model; and

[0080] a recognition module, which is used for recognizing the underground target layer according to the electromagnetic diffusion characteristics of the multi-line source ground-borehole transient electromagnetic responses.

[0081] The application has the following advantages.

[0082] In order to solve shortcomings of the ground-borehole transient electromagnetic method in detecting high resistance targets, the application provides a ground-borehole transient electromagnetic detection method based on a plurality of ground emission sources. The application analyzes distribution characteristics of various components of the multi-line source ground-borehole transient electromagnetic field and the resolution of an electrical structure, and proposes to realize a more precise recognition of underground high resistance target information by using the vertical electric field component of a pure transverse magnetic field and the horizontal electric field in some areas, so as to further improve the detection capability of the ground-borehole transient electromagnetic method for the high resistance targets.

[0083] Obviously, those skilled in the art may make various changes and modifications to the application without departing from a spirit and a scope of the application. It should be noted that the above description is only a specific embodiment of the application, and does not limit the application. Any modulation and optimization within the spirit and a principle of the application should be covered by the claims of application.