A deep metal mine exploration device and method

Abstract

This article provides a deep metal mine exploration device and method. A test well is set up within the work area. The device includes: an excitation module, a collection module and a data processing module; the excitation module is used to emit excitation signals at several excitation points set up in the test well along the well depth direction; the collection module is used to collect excitation information generated by the excitation signal and transmit the excitation information to the data processing module; the data processing module is used to: obtain the layer differential resistivity coefficient and layer differential polarizability coefficient between adjacent excitation points through the excitation information; obtain the development area of ​​deep metal mines in the work area through the layer differential resistivity coefficient and layer differential polarizability coefficient between adjacent excitation points. Thereby improving the recognition accuracy of the distribution range and depth of deep metal minerals.

Description

Technical Field

[0001] The present invention relates to the field of geophysical exploration, and in particular to an exploration device and method for deep metal mines. Background Art

[0002] Among the existing exploration methods, the ground transient electromagnetic method is often used to explore metal minerals. The ground transient electromagnetic method supplies pulse current to the underground through electric dipoles to generate a primary magnetic field in the surrounding space. This magnetic field generates an induced current in the underground conductor. After the current in the power supply circuit is disconnected, the primary magnetic field decays rapidly. In order to prevent this decay, a secondary electromagnetic field is generated in the underground medium, the so-called transient electromagnetic field. The relevant information of the metal mineral is obtained through the change of the transient electromagnetic field. However, as the buried depth of metal minerals increases, the resolution drops sharply, resulting in the inability of conventional ground transient electromagnetic methods to distinguish metal minerals at a depth of 1000m-3000m.

[0003] On this basis, for deep metal mines with a depth of 1000m-3000m, drilling grids and other methods are commonly used in the prior art to explore deep metal mines. The denser the drilling grid density, the higher the exploration accuracy, and the higher the construction cost and the longer the construction period. However, due to various factors such as the actual construction cost, construction period, and construction conditions, the density of the drilling grid cannot meet actual needs. Most deep metal mines are deeply buried underground and have complex structures, which makes it impossible for the commonly used drilling grid method to accurately explore the specific distribution of deep metal mines.

[0004] Therefore, there is an urgent need for a deep metal mine exploration device to improve the accuracy of deep metal mine exploration, and then improve the accuracy of identifying the distribution range and depth of deep metal minerals. Summary of the invention

[0005] The purpose of the embodiments of this article is to provide a deep metal mine exploration device and method to improve the recognition accuracy of the distribution range and depth of deep metal minerals.

[0006] To achieve the above-mentioned purpose, on the one hand, the embodiment of this article provides a deep metal mine exploration device, a test well is set within the work area, and the device includes: an excitation module, a collection module and a data processing module;

[0007] The excitation module is used to emit excitation signals at a plurality of excitation points arranged in the test well along the well depth direction;

[0008] The acquisition module is used to acquire the excitation information generated by the excitation signal and transmit the excitation information to the data processing module;

[0009] The data processing module is used to:

[0010] Obtaining the layer differential resistivity coefficient and layer differential polarizability coefficient between adjacent excitation points through the excitation information;

[0011] The development area of ​​deep metal ore in the work area is obtained through the layer differential resistivity coefficient and the layer differential polarizability coefficient between adjacent excitation points.

[0012] Preferably, the collecting of the excitation information generated by the excitation signal includes:

[0013] After each of the excitation points emits the excitation signal, horizontal radial electric field components generated by the excitation signal are collected at a plurality of observation points above the working area, wherein the horizontal radial electric field components are the excitation information.

[0014] Preferably, obtaining the layer differential resistivity coefficient between adjacent excitation points through the excitation information includes:

[0015] Obtaining a standard electric field component at the observation point by using the horizontal radial electric field components collected at a number of the observation points;

[0016] Obtaining, by using the standard electric field components at the observation points, an attenuation curve of the standard electric field components corresponding to each of the observation points with respect to time;

[0017] Determine the apparent resistivity of each of the excitation points according to the excitation signal related information of each of the excitation points; wherein the excitation signal related information includes an attenuation curve of the standard electric field component corresponding to the observation point with respect to time;

[0018] According to the apparent resistivity of two adjacent excitation points at the same time and the distances from the two adjacent excitation points to the wellhead, the instantaneous layer differential resistivity coefficients of the two adjacent excitation points are determined;

[0019] The layer differential resistivity coefficient between adjacent excitation points is determined by the instantaneous layer differential resistivity coefficient of two adjacent excitation points.

[0020] Preferably, the excitation module comprises:

[0021] A first excitation electrode disposed outside the test well, and a second excitation electrode disposed inside the test well, wherein the plurality of excitation points are located on the second excitation electrode;

[0022] The excitation signal related information also includes:

[0023] The position of the first excitation electrode, the position of the observation point and the position of the excitation point.

[0024] Preferably, determining the apparent resistivity of each of the excitation points based on the excitation signal-related information of each of the excitation points comprises:

[0025] The apparent resistivity of each of the excitation points is determined by the following formula:

[0026] Res x，y，z =K x，y，z ×|Es x，y，z |×10 -6 ;

[0027]

[0028] Where (x, y) is the observation point number obtained from the position of the observation point, z is the excitation point number obtained from the position of the excitation point, and K x，y，z is the device coefficient, r x，y，z is the distance from the first excitation electrode to the observation point obtained by the position of the first excitation electrode and the position of the observation point, R x，y，z Es is the distance from the excitation point to the observation point obtained by the positions of the excitation point and the observation point. x，y，z is the standard electric field component obtained from the attenuation curve of the observation point, Res x，y，z is the apparent resistivity of the excitation point numbered z.

[0029] Preferably, determining the instantaneous layer differential resistivity coefficient of two adjacent excitation points based on the apparent resistivity of the two adjacent excitation points at the same time and the distances from the two adjacent excitation points to the wellhead respectively comprises:

[0030] The instantaneous layer differential resistivity coefficient of two adjacent excitation points is determined by the following formula:

[0031] R x，y，z12 =(Res x，y，z1 ×L z2 ) / (Res x，y，z2 ×L z1 );

[0032] Among them, R x，y，z12 is the instantaneous layer differential resistivity coefficient of two adjacent excitation points numbered z1 and z2, Res x，y，z1 is the apparent resistivity of the excitation point numbered z1, Res x，y，z2 is the apparent resistivity of the excitation point numbered z2, L z1 is the distance from the excitation point numbered z1 to the wellhead, L z2 is the distance from the excitation point numbered z2 to the wellhead;

[0033] The step of determining the layer differential resistivity coefficient between adjacent excitation points by using the instantaneous layer differential resistivity coefficient of two adjacent excitation points comprises:

[0034] Among the instantaneous layer differential resistivity coefficients of two adjacent excitation points, the maximum value of the instantaneous layer differential resistivity coefficient in the time domain is selected as the layer differential resistivity coefficient between the adjacent excitation points, wherein the time domain is from the initial moment of the attenuation curve of the observation point to the moment corresponding to the maximum slope.

[0035] Preferably, obtaining the layer differential polarizability coefficient between adjacent excitation points through the excitation information includes:

[0036] Obtaining a standard electric field component at the observation point by using the horizontal radial electric field components collected at a number of the observation points;

[0037] Obtaining, by using the standard electric field components at the observation points, an attenuation curve of the standard electric field components corresponding to each of the observation points with respect to time;

[0038] Determine the apparent polarizability of each of the excitation points according to the excitation signal related information of each of the excitation points; wherein the excitation signal related information includes the position of the observation point, the position of the excitation point, and the attenuation curve of the standard electric field component corresponding to the observation point with respect to time;

[0039] The layer differential polarizability coefficient between adjacent excitation points is determined by the apparent polarizability of each of the excitation points and the distances from two adjacent excitation points to the wellhead.

[0040] Preferably, determining the apparent polarizability of each of the excitation points based on the excitation signal-related information of each of the excitation points includes:

[0041] The apparent polarizability of each excitation point is determined by the following formula:

[0042] Eta x，y，z =|Es T0 –Es T1 | / |Es T0 |;

[0043] Where (x, y) is the observation point number obtained from the position of the observation point, z is the excitation point number obtained from the position of the excitation point, T0 is the time corresponding to the maximum value of the standard electric field component obtained from the decay curve of the observation point, T1 is the time corresponding to the maximum value of the slope obtained from the decay curve of the observation point, Es T0 is the standard electric field component in the attenuation curve of the observation point corresponding to time T0, Es T1 is the standard electric field component in the attenuation curve of the observation point corresponding to time T1, Eta x，y，z is the apparent polarizability of the excitation point numbered z.

[0044] Preferably, determining the layer differential polarizability coefficient between adjacent excitation points through the apparent polarizability of each of the excitation points and the distances from two adjacent excitation points to the wellhead comprises:

[0045] The layer differential polarizability coefficient between adjacent excitation points is determined by the following formula:

[0046] Eta x，y，z12 =(Eta x，y，z1 ×L z2 ) / (Eta x，y，z2 ×L z1 );

[0047] Among them, Eta x，y，z12 is the layer differential polarizability coefficient between two adjacent excitation points numbered z1 and z2, Eta x，y，z1 is the apparent polarizability of the excitation point numbered z1, Eta x，y，z2 is the apparent polarizability of the excitation point numbered z2, L z1 is the distance from the excitation point numbered z1 to the wellhead, L z2 is the distance from the excitation point numbered z2 to the wellhead.

[0048] On the other hand, the embodiment of this invention also provides a method for exploring deep metal mines, wherein a test well is set within the working area, and the method comprises:

[0049] Transmitting excitation signals at a plurality of excitation points arranged along the depth direction of the well in the test well;

[0050] collecting excitation information generated by the excitation signal;

[0051] Obtaining the layer differential resistivity coefficient and layer differential polarizability coefficient between adjacent excitation points through the excitation information;

[0052] The development area of ​​deep metal ore in the work area is obtained through the layer differential resistivity coefficient and the layer differential polarizability coefficient between adjacent excitation points.

[0053] It can be seen from the technical solution provided in the above embodiments of this invention that the embodiments of this invention can obtain the layer differential resistivity coefficient and layer differential polarizability coefficient between any two adjacent excitation points by using the excitation information through the data processing module, thereby reflecting the stratigraphic differences between the excitation points at different depths and the excitation points adjacent thereto. Among them, compared with non-metallic minerals, the layer differential resistivity coefficient is small and the layer differential polarizability coefficient is large in the stratum where the metal mineral is located, therefore, after obtaining the layer differential resistivity coefficient and layer differential polarizability coefficient between adjacent excitation points, the area where the layer differential resistivity coefficient is small and the layer differential polarizability is large can be selected, that is, the development area of ​​deep metal minerals in the work area.

[0054] In order to make the above and other purposes, features and advantages of this article more obvious and easy to understand, the following specifically cites preferred embodiments and describes them in detail with reference to the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0055] In order to more clearly illustrate the technical solutions in the embodiments of this article or the prior art, the drawings required for use in the embodiments or the prior art descriptions are briefly introduced below. Obviously, the drawings described below are only some embodiments of this article. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0056] Figure 1 A schematic diagram of the module structure of a deep metal mine exploration device provided in an embodiment of this invention is shown;

[0057] Figure 2 Another schematic diagram of the module structure of a deep metal mine exploration device provided in an embodiment of this invention is shown;

[0058] Figure 3 A schematic diagram of another module structure of a deep metal mine exploration device provided in an embodiment of this invention is shown;

[0059] Figure 4 A schematic diagram of a process for exploring a deep metal mine provided in an embodiment of this invention is shown;

[0060] Figure 5 A schematic diagram of a process for obtaining the layer differential resistivity coefficient between adjacent excitation points in a deep metal mine exploration method provided in an embodiment of this invention is shown;

[0061] Figure 6 A schematic diagram of a process for obtaining layer differential polarizability coefficients between adjacent excitation points in a method for exploring deep metal mines provided in an embodiment of this invention is shown;

[0062] Figure 7 A schematic diagram of the structure of a computer device provided in an embodiment of this invention is shown.

[0063] Description of the accompanying symbols:

[0064] 100. Preprocessing module;

[0065] 200, excitation module;

[0066] 210, transmitter;

[0067] 220, a first exciting electrode;

[0068] 230. a second exciting electrode;

[0069] 300, acquisition module;

[0070] 310, collection station;

[0071] 400. Data processing module;

[0072] 410, pre-processing unit;

[0073] 420, computing unit;

[0074] 702. Computer equipment;

[0075] 704, processor;

[0076] 706. Memory;

[0077] 708, driving mechanism;

[0078] 710, input / output module;

[0079] 712. Input devices;

[0080] 714. Output devices;

[0081] 716. Presentation equipment;

[0082] 718. Graphical user interface;

[0083] 720, network interface;

[0084] 722. Communication link;

[0085] 724. Communication bus. DETAILED DESCRIPTION

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

[0087] Among the existing exploration methods, the ground transient electromagnetic method is often used to explore metal minerals. The ground transient electromagnetic method supplies pulse current to the underground through electric dipoles to generate a primary magnetic field in the surrounding space. This magnetic field generates an induced current in the underground conductor. After the current in the power supply circuit is disconnected, the primary magnetic field decays rapidly. In order to prevent this decay, a secondary electromagnetic field is generated in the underground medium, the so-called transient electromagnetic field. The relevant information of the metal minerals is obtained through the changes in the transient electromagnetic field. However, as the buried depth of metal minerals increases, the resolution drops sharply, resulting in the inability of conventional ground transient electromagnetic methods to distinguish metal minerals at a depth of 1000m-3000m.

[0088] In order to solve the above problems, the embodiments of this article provide a deep metal mine exploration device, which can improve the recognition accuracy of the distribution range and depth of deep metal minerals.

[0089] Reference Figure 1 , a deep metal mine exploration device, a test well is set within the work area, the device includes: an excitation module 200, a collection module 300 and a data processing module 400;

[0090] The excitation module 200 is used to emit excitation signals at a plurality of excitation points arranged along the well depth direction in the test well;

[0091] The acquisition module 300 is used to acquire the excitation information generated by the excitation signal and transmit the excitation information to the data processing module 400;

[0092] The data processing module 400 is used to:

[0093] Obtaining the layer differential resistivity coefficient and layer differential polarizability coefficient between adjacent excitation points through the excitation information;

[0094] The development area of ​​deep metal ore in the work area is obtained through the layer differential resistivity coefficient and the layer differential polarizability coefficient between adjacent excitation points.

[0095] In the above embodiment, the data processing module 400 uses the excitation information to obtain the layer differential resistivity coefficient and layer differential polarizability coefficient between any two adjacent excitation points, thereby reflecting the formation differences between the excitation points at different depths and the excitation points adjacent to them. The apparent resistivity calculated at an excitation point reflects the resistivity characteristics of the formation above the excitation point. The ratio of the electrode distances between adjacent excitation points after normalization can reflect the formation resistivity characteristics between the two excitation points, which we define as the layer differential resistivity coefficient; similarly, the apparent polarizability ratio after the electrode distances between adjacent excitation points after normalization is defined as the layer differential polarizability coefficient, which can reflect the formation polarizability characteristics between the two excitation points. Among them, compared with non-metallic minerals, by comparing all the layer differential resistivity coefficients and layer differential polarizability coefficients of the test wells, the layer differential resistivity coefficient of the stratum where the metal mine is located is small, and the layer differential polarizability coefficient is large. Therefore, after obtaining the layer differential resistivity coefficient and layer differential polarizability coefficient between adjacent excitation points, the area with small layer differential resistivity coefficient and large layer differential polarizability can be selected, which is the development area of ​​deep metal mines in the work area.

[0096] In the embodiments of this article, the work area refers to the underground working area for metal mine exploration, and test wells are set up within the work area.

[0097] Reference Figure 2Preferably, the excitation module 200 further includes a pre-processing module 100, and the pre-processing module 100 is used to:

[0098] First, determine the top depth, bottom depth and thickness of the metal mine. If the metal mine is found in the test well, the top depth, bottom depth and thickness of the metal mine are determined based on the logging resistivity; if the metal mine cannot be found in the test well, the top depth, bottom depth and thickness of the metal mine within the work area are inferred through the seismic data or electrical data near the well.

[0099] Then, a formation resistivity model is established based on the logging resistivity. Specifically, the logging resistivity is input into the resistivity model generation software to generate a formation resistivity model. The horizontal coordinate of the formation resistivity model is depth, and the vertical coordinate is simulated resistivity, which can reflect the distribution of the formation as the well depth changes, and then preliminarily distinguish each formation. If a metal mine is found in the test well, the simulated resistivity corresponding to the location of the metal mine can be determined because the top depth, bottom depth and thickness of the metal mine are known; if the metal mine cannot be found in the test well, the simulated resistivity corresponding to the location of the metal mine can also be determined by the known top depth, bottom depth and thickness of the metal mine.

[0100] Finally, the excitation signal period of the positive and negative zero-crossing mode required by the target metal mine is calculated according to the formation resistivity model. Specifically, the excitation signal period of the positive and negative zero-crossing mode required by the target metal mine is calculated by the following formula, t i =2×10 -6 ×H i ×S i , where i is a positive integer, t i is the excitation signal period required for the i-th layer, H i is the thickness of the i-th layer, S i is the conductivity of the i-th layer, the conductivity of the i-th layer is the reciprocal of the resistivity corresponding to the i-th layer in the formation resistivity model, and the excitation signal period required for detecting metal ore targets is the sum of the signal periods required for the metal ore target layer and all layers above it.

[0101] Reference Figure 3 Preferably, the excitation module 200 comprises:

[0102] A transmitter 210 for transmitting an excitation signal, a first excitation electrode 220 arranged outside the test well, and a second excitation electrode 230 arranged inside the test well, wherein a plurality of excitation points are located on the second excitation electrode 230 .

[0103] Specifically, the transmitter 210 includes an A end and a B end. The A end of the transmitter 210 is connected to the first excitation electrode 220, and the B end is connected to the second excitation electrode 230, forming a closed loop, which can transmit an excitation signal. The end of the first excitation electrode 220 away from the A end of the transmitter 210 is buried in the ground within 30m of the test well. The first excitation electrode 220 can be a well-grounded material such as aluminum foil, steel pipe or copper rod. The second excitation electrode 230 is connected to the B end of the transmitter 210 through an excitation cable, wherein the excitation cable is made of a single-core cable that is resistant to high temperature, voltage and current. The second excitation electrode 230 is made of a conductive alloy that is resistant to high temperature, voltage, current and oxidation, and the second excitation electrode 230 is located in the test well.

[0104] In the embodiments of this invention, the spacing between adjacent excitation points is equal and the spacing is not greater than the vertical resolution accuracy. Specifically, if a metal mine is found in the test well, the top depth of the metal mine is known, the depth of the excitation point located at the top end of the well should not be lower than the top depth of the metal mine plus 5 times the spacing between the excitation points, the excitation point located at the bottom end of the well should not be higher than the bottom depth of the metal mine plus 5 times the spacing between the excitation points, and the excitation point located in the middle should be located in the middle of the top and bottom surfaces of the metal mine. If the metal mine cannot be found in the test well, since the top depth, bottom depth and thickness of the metal mine are known, the depth position of the metal mine corresponding to the test well can be obtained, the depth of the excitation point located at the top end of the well should not be lower than the top depth of the metal mine plus 5 times the spacing between the excitation points, the excitation point located at the bottom end of the well should not be higher than the bottom depth of the metal mine plus 5 times the spacing between the excitation points, and the excitation point located in the middle should be located in the middle of the top and bottom surfaces of the metal mine.

[0105] In the embodiment of this article, according to the position of the excitation point in the well, starting from the excitation point at the bottom of the well, the excitation signal is emitted from bottom to top, and each excitation point is excited according to the excitation signal period of the target metal ore when emitting the excitation signal. Specifically, the excitation signal is a square wave signal in a positive and negative zero-crossing mode. For example, the calculated excitation signal period is 4s, and the excitation current is 40A, then the actual excitation signal is 40A for 0-1s, 0A for 1-2s, -40A for 2-3s, and 0A for 3-4s.

[0106] In the embodiment of this article, the step of collecting the excitation information generated by the excitation signal includes:

[0107] After each of the excitation points emits the excitation signal, horizontal radial electric field components generated by the excitation signal are collected at a plurality of observation points above the working area, wherein the horizontal radial electric field components are the excitation information.

[0108] Specifically, the acquisition module 300 includes a plurality of acquisition stations 310 , which are respectively arranged at observation points. The acquisition stations 310 collect horizontal radial electric field components at corresponding observation points and transmit the horizontal radial electric field components to the data processing module 400 .

[0109] In order to ensure the accuracy of receiving the excitation signal, the observation points are all located within 2 kilometers of the test well. When deploying observation points in a three-dimensional observation manner, several observation points are arranged in an array to cover the ground above the work area, and the distance between two adjacent observation points is no greater than the horizontal resolution accuracy.

[0110] In the embodiment of this article, the data processing module 400 includes a pre-processing unit 410 and a computing unit 420 , and data can be transmitted between the pre-processing unit 410 and the computing unit 420 .

[0111] Preferably, the pre-processing unit 410 is used for:

[0112] Obtaining a standard electric field component at the observation point by using the horizontal radial electric field components collected at a number of the observation points;

[0113] The attenuation curve of the standard electric field component corresponding to each observation point with respect to time is obtained through the standard electric field component at the observation point.

[0114] Specifically, first, the horizontal radial electric field component collected at the observation point is current normalized, that is, the horizontal radial electric field component is divided by the corresponding current value to obtain the standard electric field component at the observation point. Then, the standard electric field component at the observation point is simulated by the curve generation software to generate a time-dependent attenuation curve of the standard electric field component, where the horizontal coordinate of the attenuation curve is time and the vertical coordinate is the standard electric field component, reflecting the law of the change of the standard electric field component over time. Finally, the time-dependent attenuation curve of the standard electric field component corresponding to the observation point is transmitted to the calculation unit 420.

[0115] Preferably, the calculation unit 420 is used for:

[0116] Determine the apparent resistivity and apparent polarizability of each of the excitation points according to the excitation signal related information of each of the excitation points; wherein the excitation signal related information includes the position of the first excitation electrode 220, the position of the observation point, the position of the excitation point, and the attenuation curve of the standard electric field component corresponding to the observation point with respect to time;

[0117] According to the apparent resistivity of two adjacent excitation points at the same time and the distances from the two adjacent excitation points to the wellhead, the instantaneous layer differential resistivity coefficients of the two adjacent excitation points are determined;

[0118] Determine the layer differential resistivity coefficient between adjacent excitation points through the instantaneous layer differential resistivity coefficient of two adjacent excitation points;

[0119] Determine the layer differential polarizability coefficient between adjacent excitation points by using the apparent polarizability of each of the excitation points and the distances from two adjacent excitation points to the wellhead;

[0120] The development area of ​​deep metal ore in the work area is obtained through the layer differential resistivity coefficient and the layer differential polarizability coefficient between adjacent excitation points.

[0121] Furthermore, determining the apparent resistivity of each of the excitation points based on the excitation signal-related information of each of the excitation points includes:

[0122] The apparent resistivity of each of the excitation points is determined by the following formula:

[0123] Res x，y，z =K x，y，z ×|Es x，y，z |×10 -6 ;

[0124]

[0125] Where (x, y) is the observation point number obtained from the position of the observation point, z is the excitation point number obtained from the position of the excitation point, and K x，y，z is the device coefficient, r x，y，z is the distance from the first excitation electrode 220 to the observation point obtained from the position of the first excitation electrode 220 and the position of the observation point, R x，y，z Es is the distance from the excitation point to the observation point obtained by the positions of the excitation point and the observation point. x，y，z is the standard electric field component obtained from the attenuation curve of the observation point, Res x，y，z is the apparent resistivity of the excitation point numbered z.

[0126] Specifically, since the observation points are arranged in the form of an array, any observation point can be represented in the form of coordinates (x, y), where x can be the number of rows and y can be the number of columns. Since the excitation point is located in the test well, the excitation point can be represented as 1, 2, 3... from the bottom of the well to the wellhead, z = 1 or 2 or 3... Since the position of the first excitation electrode 220 and the position of the observation point are both known, the distance between the two can be obtained. Since the position of the excitation point and the position of the observation point are both known, the distance between the two can be obtained. That is, the distance from any observation point to the first excitation electrode 220 and the distance from any excitation point to any observation point can be obtained. Since the standard electric field component in the attenuation curve of the observation point is a quantity that varies with time, the apparent resistivity of the excitation point numbered z is also a quantity that varies with time.

[0127] Furthermore, the instantaneous layer differential resistivity coefficient of two adjacent excitation points is determined based on the apparent resistivity of the two adjacent excitation points at the same time and the distances from the two adjacent excitation points to the wellhead, including:

[0128] The instantaneous layer differential resistivity coefficient of two adjacent excitation points is determined by the following formula:

[0129] R x，y，z12 =(Res x，y，z1 ×L z2 ) / (Res x，y，z2 ×L z1 );

[0130] Among them, R x，y，z12 is the instantaneous layer differential resistivity coefficient of two adjacent excitation points numbered z1 and z2, Res x，y，z1 is the apparent resistivity of the excitation point numbered z1, Res x，y，z2 is the apparent resistivity of the excitation point numbered z2, L z1 is the distance from the excitation point numbered z1 to the wellhead, L z2 is the distance from the excitation point numbered z2 to the wellhead.

[0131] Specifically, the apparent resistivities of two adjacent excitation points z1 and z2 are quantities that vary with time, and the instantaneous layer differential resistivity coefficients of z1 and z2 obtained by the above formula are quantities that vary with time.

[0132] Furthermore, determining the layer differential resistivity coefficient between adjacent excitation points by using the instantaneous layer differential resistivity coefficient of two adjacent excitation points includes:

[0133] Among the instantaneous layer differential resistivity coefficients of two adjacent excitation points, the maximum value of the instantaneous layer differential resistivity coefficient in the time domain is selected as the layer differential resistivity coefficient between the adjacent excitation points, wherein the time domain is from the initial moment of the attenuation curve of the observation point to the moment corresponding to the maximum slope.

[0134] Specifically, the layer differential resistivity coefficient of two adjacent excitation points is a quantity that changes with time. Therefore, the maximum value of the instantaneous layer differential resistivity coefficient in the time domain is selected as the layer differential resistivity coefficient between the adjacent excitation points. As time goes by, the excitation signal decays more slowly and the signal-to-noise ratio becomes lower, resulting in inaccurate calculation results. Therefore, it is necessary to select the time domain from the initial moment of the attenuation curve of the observation point to the moment corresponding to the maximum slope.

[0135] In addition, determining the apparent polarizability of each of the excitation points based on the excitation signal related information of each of the excitation points includes:

[0136] The apparent polarizability of each excitation point is determined by the following formula:

[0137] Eta x , y , z =|Es T0 –Es T1 | / |Es T0 |;

[0138] Among them, T0 is the time corresponding to the maximum value of the standard electric field component obtained from the attenuation curve of the observation point, T1 is the time corresponding to the maximum value of the slope obtained from the attenuation curve of the observation point, and Es T0 is the standard electric field component in the attenuation curve of the observation point corresponding to time T0, Es T1 is the standard electric field component in the attenuation curve of the observation point corresponding to time T1, Eta x，y，z is the apparent polarizability of the excitation point numbered z.

[0139] Furthermore, the layer differential polarizability coefficient between adjacent excitation points is determined by the apparent polarizability of each of the excitation points and the distances from two adjacent excitation points to the wellhead, including:

[0140] The layer differential polarizability coefficient between adjacent excitation points is determined by the following formula:

[0141] Eta x，y，z12 =(Eta x，y，z1 ×L z2 ) / (Eta x，y，z2 ×L z1 );

[0142] Among them, Eta x，y，z12is the layer differential polarizability coefficient between two adjacent excitation points numbered z1 and z2, Eta x，y，z1 is the apparent polarizability of the excitation point numbered z1, Eta x，y，z2 is the apparent polarizability of the excitation point numbered z2, L z1 is the distance from the excitation point numbered z1 to the wellhead, L z2 is the distance from the excitation point numbered z2 to the wellhead.

[0143] Specifically, the layer differential resistivity coefficient and the layer differential polarizability coefficient between adjacent excitation points are obtained by the above formula.

[0144] Furthermore, the area with small layer differential resistivity coefficient and large layer differential polarizability is selected as the development area of ​​deep metal mines in the work area.

[0145] Specifically, compared with the non-metallic area in the stratum, the stratum where the metal mine is located has a small layer differential resistivity coefficient and a large layer differential polarizability coefficient. After obtaining the layer differential resistivity coefficient and layer differential polarizability coefficient between adjacent excitation points, the area where the area with a small layer differential resistivity coefficient overlaps with the area with a large layer differential polarizability is the development area of ​​deep metal mines in the work area.

[0146] Based on the above-mentioned deep metal mine exploration device, the embodiments of this article also provide a deep metal mine exploration method.

[0147] Reference Figure 4 , a test well is set within the work area, and the method may include the following steps:

[0148] S101: transmitting excitation signals at a plurality of excitation points arranged along the depth direction of the test well;

[0149] S102: collecting excitation information generated by the excitation signal;

[0150] S103: obtaining the layer differential resistivity coefficient and the layer differential polarizability coefficient between adjacent excitation points through the excitation information;

[0151] S104: The development area of ​​deep metal ore in the work area is obtained through the layer differential resistivity coefficient and the layer differential polarizability coefficient between adjacent excitation points.

[0152] Furthermore, the collecting of the excitation information generated by the excitation signal includes:

[0153] After each of the excitation points emits the excitation signal, horizontal radial electric field components generated by the excitation signal are collected at a plurality of observation points above the working area, wherein the horizontal radial electric field components are the excitation information.

[0154] Reference Figure 5 Further, the layer differential resistivity coefficient between adjacent excitation points is obtained through the excitation information, including:

[0155] S1031: Obtaining a standard electric field component at the observation point through the horizontal radial electric field components collected at the plurality of observation points;

[0156] S1032: Obtaining, by using the standard electric field components at the observation points, an attenuation curve of the standard electric field components corresponding to each of the observation points with respect to time;

[0157] S1033: determining the apparent resistivity of each of the excitation points according to the excitation signal related information of each of the excitation points; wherein the excitation signal related information includes an attenuation curve of the standard electric field component corresponding to the observation point with respect to time;

[0158] S1034: Determine the instantaneous layer differential resistivity coefficients of the two adjacent excitation points according to the apparent resistivity of the two adjacent excitation points at the same time and the distances from the two adjacent excitation points to the wellhead;

[0159] S1035: Determine the layer differential resistivity coefficient between adjacent excitation points through the instantaneous layer differential resistivity coefficient of two adjacent excitation points.

[0160] Furthermore, a first excitation electrode 220 is arranged outside the test well, and a second excitation electrode 230 is arranged inside the test well, and a plurality of the excitation points are located on the second excitation electrode 230;

[0161] The excitation signal related information also includes:

[0162] The position of the first excitation electrode 220, the position of the observation point and the position of the excitation point.

[0163] Furthermore, determining the apparent resistivity of each of the excitation points based on the excitation signal-related information of each of the excitation points includes:

[0164] The apparent resistivity of each of the excitation points is determined by the following formula:

[0165] Res x，y，z =K x，y，z ×|Es x，y，z |×10 -6 ;

[0166]

[0167] Where (x, y) is the observation point number obtained from the position of the observation point, z is the excitation point number obtained from the position of the excitation point, and K x，y，zis the device coefficient, r x，y，z is the distance from the first excitation electrode 220 to the observation point obtained from the position of the first excitation electrode 220 and the position of the observation point, R x，y，z Es is the distance from the excitation point to the observation point obtained by the positions of the excitation point and the observation point. x，y，z is the standard electric field component obtained from the attenuation curve of the observation point, Res x，y，z is the apparent resistivity of the excitation point numbered z.

[0168] Furthermore, the instantaneous layer differential resistivity coefficient of two adjacent excitation points is determined based on the apparent resistivity of the two adjacent excitation points at the same time and the distances from the two adjacent excitation points to the wellhead, including:

[0169] The instantaneous layer differential resistivity coefficient of two adjacent excitation points is determined by the following formula:

[0170] R x，y，z12 =(Res x，y，z1 ×L z2 ) / (Res x，y，z2 ×L z1 );

[0171] Among them, R x，y，z12 is the instantaneous layer differential resistivity coefficient of two adjacent excitation points numbered z1 and z2, Res x，y，z1 is the apparent resistivity of the excitation point numbered z1, Res x，y，z2 is the apparent resistivity of the excitation point numbered z2, L z1 is the distance from the excitation point numbered z1 to the wellhead, L z2 is the distance from the excitation point numbered z2 to the wellhead;

[0172] The step of determining the layer differential resistivity coefficient between adjacent excitation points by using the instantaneous layer differential resistivity coefficient of two adjacent excitation points comprises:

[0173] Among the instantaneous layer differential resistivity coefficients of two adjacent excitation points, the maximum value of the instantaneous layer differential resistivity coefficient in the time domain is selected as the layer differential resistivity coefficient between the adjacent excitation points, wherein the time domain is from the initial moment of the attenuation curve of the observation point to the moment corresponding to the maximum slope.

[0174] Reference Figure 6 Further, the layer differential polarizability coefficient between adjacent excitation points is obtained through the excitation information, including:

[0175] S1036: Obtaining a standard electric field component at the observation point through the horizontal radial electric field components collected at the plurality of observation points;

[0176] S1037: Obtaining, by using the standard electric field components at the observation points, an attenuation curve of the standard electric field components corresponding to each of the observation points with respect to time;

[0177] S1038: determining the apparent polarizability of each of the excitation points according to the excitation signal related information of each of the excitation points; wherein the excitation signal related information includes the position of the observation point, the position of the excitation point, and the attenuation curve of the standard electric field component corresponding to the observation point with respect to time;

[0178] S1039: Determine the layer differential polarizability coefficient between adjacent excitation points through the apparent polarizability of each of the excitation points and the distances from two adjacent excitation points to the wellhead.

[0179] Furthermore, determining the apparent polarizability of each of the excitation points based on the excitation signal related information of each of the excitation points includes:

[0180] The apparent polarizability of each excitation point is determined by the following formula:

[0181] Eta x，y，z =|Es T0 –Es T1 | / |Es T0 |;

[0182] Where (x, y) is the observation point number obtained from the position of the observation point, z is the excitation point number obtained from the position of the excitation point, T0 is the time corresponding to the maximum value of the standard electric field component obtained from the decay curve of the observation point, T1 is the time corresponding to the maximum value of the slope obtained from the decay curve of the observation point, Es T0 is the standard electric field component in the attenuation curve of the observation point corresponding to time T0, Es T1 is the standard electric field component in the attenuation curve of the observation point corresponding to time T1, Eta x，y，z is the apparent polarizability of the excitation point numbered z.

[0183] Furthermore, the layer differential polarizability coefficient between adjacent excitation points is determined by the apparent polarizability of each of the excitation points and the distances from two adjacent excitation points to the wellhead, including:

[0184] The layer differential polarizability coefficient between adjacent excitation points is determined by the following formula:

[0185] Eta x，y，z12 =(Eta x，y，z1 ×L z2 ) / (Eta x，y，z2 ×L z1 );

[0186] Among them, Etax，y，z12 is the layer differential polarizability coefficient between two adjacent excitation points numbered z1 and z2, Eta x，y，z1 is the apparent polarizability of the excitation point numbered z1, Eta x，y，z2 is the apparent polarizability of the excitation point numbered z2, L z1 is the distance from the excitation point numbered z1 to the wellhead, L z2 is the distance from the excitation point numbered z2 to the wellhead.

[0187] In one embodiment of this invention, reference is made to Figure 7 As shown, a computer device 702 is also provided. The computer device 702 may include one or more processors 704, such as one or more central processing units (CPUs) or graphics processing units (GPUs), and each processing unit may implement one or more hardware threads. The computer device 702 may also include any memory 706, which is used to store any kind of information such as code, settings, data, etc. In a specific embodiment, a computer program on the memory 706 and executable on the processor 704, when the computer program is executed by the processor 704, can execute instructions according to the above method. Non-limiting, for example, the memory 706 may include any one or more combinations of the following: any type of RAM, any type of ROM, flash memory device, hard disk, optical disk, etc. More generally, any memory may use any technology to store information. Further, any memory may provide volatile or non-volatile retention of information. Further, any memory may represent a fixed or removable component of the computer device 702. In one case, when the processor 704 executes an associated instruction stored in any memory or combination of memories, the computer device 702 may perform any operation of the associated instruction. The computer device 702 also includes one or more drive mechanisms 708 for interacting with any storage, such as a hard disk drive mechanism, an optical disk drive mechanism, etc.

[0188] The computer device 702 may also include an input / output module 710 (I / O) for receiving various inputs (via input devices 712) and for providing various outputs (via output devices 714). A specific output mechanism may include a presentation device 716 and an associated graphical user interface 718 (GUI). In other embodiments, the input / output module 710 (I / O), input device 712, and output device 714 may not be included, and the computer device 702 may be used as a computer device in a network. The computer device 702 may also include one or more network interfaces 720 for exchanging data with other devices via one or more communication links 722. One or more communication buses 724 couple the components described above together.

[0189] The communication link 722 may be implemented in any manner, for example, through a local area network, a wide area network (e.g., the Internet), a point-to-point connection, etc., or any combination thereof. The communication link 722 may include any combination of hardwired links, wireless links, routers, gateway functions, name servers, etc. governed by any protocol or combination of protocols.

[0190] Corresponds to Figure 4-Figure 6 The method in this article also provides a computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, the steps of the above method are executed.

[0191] The embodiment of the present invention also provides a computer readable instruction, wherein when the processor executes the instruction, the program therein causes the processor to execute the following Figures 4 to 6 The method shown.

[0192] It should be understood that in the various embodiments of this document, the size of the serial numbers of the above-mentioned processes does not mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of this document.

[0193] It should also be understood that in the embodiments of this article, the term "and / or" is only a description of the association relationship of the associated objects, indicating that three relationships may exist. For example, A and / or B can represent: A exists alone, A and B exist at the same time, and B exists alone. In addition, the character " / " in this article generally indicates that the associated objects before and after are in an "or" relationship.

[0194] Those of ordinary skill in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the composition and steps of each example have been generally described in terms of function in the above description. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of this article.

[0195] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the systems, devices and units described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0196] In the several embodiments provided herein, it should be understood that the disclosed systems, devices and methods can be implemented in other ways. For example, the device embodiments described above are only schematic, for example, the division of the units is only a logical function division, and there may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the mutual coupling or direct coupling or communication connection shown or discussed can be an indirect coupling or communication connection through some interfaces, devices or units, or can be electrical, mechanical or other forms of connection.

[0197] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the embodiments of this article.

[0198] In addition, each functional unit in each embodiment of this invention may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above integrated unit may be implemented in the form of hardware or in the form of software functional unit.

[0199] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this article is essentially or the part that contributes to the prior art, or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions to enable a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of this article. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk and other media that can store program codes.

[0200] Specific embodiments are used in this article to illustrate the principles and implementation methods of this article. The description of the above embodiments is only used to help understand the methods and core ideas of this article. At the same time, for general technicians in this field, according to the ideas of this article, there will be changes in the specific implementation methods and application scopes. In summary, the content of this specification should not be understood as a limitation on this article.

Claims

1. A deep metal mine exploration device, characterized in that: A test well is set up within the work area, and the device includes: an excitation module, a collection module and a data processing module; The excitation module is used to emit excitation signals at a plurality of excitation points arranged in the test well along the well depth direction; The acquisition module is used to acquire the excitation information generated by the excitation signal and transmit the excitation information to the data processing module; The data processing module is used to: Obtaining the layer differential resistivity coefficient and layer differential polarizability coefficient between adjacent excitation points through the excitation information; The development area of ​​deep metal ore in the work area is obtained through the layer differential resistivity coefficient and the layer differential polarizability coefficient between adjacent excitation points; The step of collecting the excitation information generated by the excitation signal includes: After each of the excitation points emits the excitation signal, horizontal radial electric field components generated by the excitation signal are collected at a plurality of observation points above the work area, wherein the horizontal radial electric field components are the excitation information; The step of obtaining the layer differential resistivity coefficient between adjacent excitation points through the excitation information includes: Obtaining a standard electric field component at the observation point by using the horizontal radial electric field components collected at a number of the observation points; Obtaining, by using the standard electric field components at the observation points, an attenuation curve of the standard electric field components corresponding to each of the observation points with respect to time; Determine the apparent resistivity of each of the excitation points according to the excitation signal related information of each of the excitation points; wherein the excitation signal related information includes an attenuation curve of the standard electric field component corresponding to the observation point with respect to time; According to the apparent resistivity of two adjacent excitation points at the same time and the distances from the two adjacent excitation points to the wellhead, the instantaneous layer differential resistivity coefficients of the two adjacent excitation points are determined; Determine the layer differential resistivity coefficient between adjacent excitation points through the instantaneous layer differential resistivity coefficient of two adjacent excitation points; The step of obtaining the layer differential polarizability coefficient between adjacent excitation points through the excitation information includes: Obtaining a standard electric field component at the observation point by using the horizontal radial electric field components collected at a number of the observation points; Obtaining, by using the standard electric field components at the observation points, an attenuation curve of the standard electric field components corresponding to each of the observation points with respect to time; Determine the apparent polarizability of each of the excitation points according to the excitation signal related information of each of the excitation points; wherein the excitation signal related information includes the position of the observation point, the position of the excitation point, and the attenuation curve of the standard electric field component corresponding to the observation point with respect to time; The layer differential polarizability coefficient between adjacent excitation points is determined by the apparent polarizability of each of the excitation points and the distances from two adjacent excitation points to the wellhead.

2. The deep metal mine exploration device according to claim 1, characterized in that: The excitation module comprises: A first excitation electrode disposed outside the test well, and a second excitation electrode disposed inside the test well, wherein the plurality of excitation points are located on the second excitation electrode; The excitation signal related information also includes: The position of the first excitation electrode, the position of the observation point and the position of the excitation point.

3. The deep metal mine exploration device according to claim 2, characterized in that: Determining the apparent resistivity of each of the excitation points based on the excitation signal related information of each of the excitation points includes: The apparent resistivity of each of the excitation points is determined by the following formula: Beef x，y，z =K x，y，z ×|It is x，y，z |×10 -6 ; Where (x, y) is the observation point number obtained from the position of the observation point, z is the excitation point number obtained from the position of the excitation point, and K x，y，z is the device coefficient, r x，y，z is the distance from the first excitation electrode to the observation point obtained by the position of the first excitation electrode and the position of the observation point, R x，y，z Es is the distance from the excitation point to the observation point obtained by the position of the excitation point and the position of the observation point. x，y，z is the standard electric field component obtained from the attenuation curve of the observation point, Res x，y，z is the apparent resistivity of the excitation point numbered z.

4. The deep metal mine exploration device according to claim 3, characterized in that: The method of determining the instantaneous layer differential resistivity coefficient of two adjacent excitation points based on the apparent resistivity of the two adjacent excitation points at the same time and the distances from the two adjacent excitation points to the wellhead comprises: The instantaneous layer differential resistivity coefficient of two adjacent excitation points is determined by the following formula: R x，y，z12 =(Res x，y，z1 ×L z2 ) / (Res x，y，z2 ×L z1 ); Among them, R x，y，z12 is the instantaneous layer differential resistivity coefficient of two adjacent excitation points numbered z1 and z2, Res x，y，z1 is the apparent resistivity of the excitation point numbered z1, Res x，y，z2 is the apparent resistivity of the excitation point numbered z2, L z1 is the distance from the excitation point numbered z1 to the wellhead, L z2 is the distance from the excitation point numbered z2 to the wellhead; The step of determining the layer differential resistivity coefficient between adjacent excitation points by using the instantaneous layer differential resistivity coefficient of two adjacent excitation points comprises: Among the instantaneous layer differential resistivity coefficients of two adjacent excitation points, the maximum value of the instantaneous layer differential resistivity coefficient in the time domain is selected as the layer differential resistivity coefficient between the adjacent excitation points, wherein the time domain is from the initial moment of the attenuation curve of the observation point to the moment corresponding to the maximum slope.

5. The deep metal mine exploration device according to claim 1, characterized in that: Determining the apparent polarizability of each of the excitation points based on the excitation signal related information of each of the excitation points includes: The apparent polarizability of each excitation point is determined by the following formula: Eta x，y，z =|It is T0 -Is T1 | / |It is T0 |; Where (x, y) is the observation point number obtained from the position of the observation point, z is the excitation point number obtained from the position of the excitation point, T0 is the time corresponding to the maximum value of the standard electric field component obtained from the decay curve of the observation point, T1 is the time corresponding to the maximum value of the slope obtained from the decay curve of the observation point, Es T0 is the standard electric field component in the attenuation curve of the observation point corresponding to time T0, Es T1 is the standard electric field component in the attenuation curve of the observation point corresponding to time T1, Eta x，y，z is the apparent polarizability of the excitation point numbered z.

6. The deep metal mine exploration device according to claim 5, characterized in that: Determining the layer differential polarizability coefficient between adjacent excitation points by using the apparent polarizability of each excitation point and the distances from two adjacent excitation points to the wellhead, comprises: The layer differential polarizability coefficient between adjacent excitation points is determined by the following formula: And x，y，z12 =(And x，y，z1 ×L z2 ) / (And x，y，z2 ×L z1 ); Among them, Eta x，y，z12 is the layer differential polarizability coefficient between two adjacent excitation points numbered z1 and z2, Eta x，y，z1 is the apparent polarizability of the excitation point numbered z1, Eta x，y，z2 is the apparent polarizability of the excitation point numbered z2, L z1 is the distance from the excitation point numbered z1 to the wellhead, L z2 is the distance from the excitation point numbered z2 to the wellhead.

7. A method for exploring deep metal mines, characterized in that: A test well is set within the work area, and the method includes: Emitting excitation signals at a plurality of excitation points arranged along the depth direction of the well in the test well; collecting excitation information generated by the excitation signal; Obtaining the layer differential resistivity coefficient and layer differential polarizability coefficient between adjacent excitation points through the excitation information; The development area of ​​deep metal ore in the work area is obtained through the layer differential resistivity coefficient and the layer differential polarizability coefficient between adjacent excitation points; The step of collecting the excitation information generated by the excitation signal includes: After each of the excitation points emits the excitation signal, horizontal radial electric field components generated by the excitation signal are collected at a plurality of observation points above the work area, wherein the horizontal radial electric field components are the excitation information; The step of obtaining the layer differential resistivity coefficient between adjacent excitation points through the excitation information includes: Obtaining a standard electric field component at the observation point by using the horizontal radial electric field components collected at a number of the observation points; Obtaining, by using the standard electric field components at the observation points, an attenuation curve of the standard electric field components corresponding to each of the observation points with respect to time; Determine the apparent resistivity of each of the excitation points according to the excitation signal related information of each of the excitation points; wherein the excitation signal related information includes an attenuation curve of the standard electric field component corresponding to the observation point with respect to time; According to the apparent resistivity of two adjacent excitation points at the same time and the distances from the two adjacent excitation points to the wellhead, the instantaneous layer differential resistivity coefficients of the two adjacent excitation points are determined; Determine the layer differential resistivity coefficient between adjacent excitation points through the instantaneous layer differential resistivity coefficient of two adjacent excitation points; The step of obtaining the layer differential polarizability coefficient between adjacent excitation points through the excitation information includes: Obtaining a standard electric field component at the observation point by using the horizontal radial electric field components collected at a number of the observation points; Obtaining, by using the standard electric field components at the observation points, an attenuation curve of the standard electric field components corresponding to each of the observation points with respect to time; Determine the apparent polarizability of each of the excitation points according to the excitation signal related information of each of the excitation points; wherein the excitation signal related information includes the position of the observation point, the position of the excitation point, and the attenuation curve of the standard electric field component corresponding to the observation point with respect to time; The layer differential polarizability coefficient between adjacent excitation points is determined by the apparent polarizability of each of the excitation points and the distances from two adjacent excitation points to the wellhead.

Images