Electromagnetic full-range apparent resistivity calculation method, device, medium and equipment

Abstract

The present disclosure relates to a method, device, medium and equipment for electromagnetic full-range apparent resistivity calculation, comprising: acquiring the vertical induced electromotive force received by each time window electric dipole element, determining the initial value of the apparent resistivity of each time window; determining the transient vertical magnetic field response and the transient vertical electromotive force response of each time window according to the measurement parameters of the electric dipole element, the initial value corresponding to each time window and the vertical induced electromotive force; determining whether to adjust the initial value according to the transient vertical electromotive force response of each time window, the corresponding vertical induced electromotive force and the preset error threshold; if it is determined to adjust the initial value of any time window, determining the adjustment strategy for the initial value according to the corresponding transient vertical magnetic field response and the vertical induced magnetic field value, recalculating the corresponding transient vertical magnetic field response and the transient vertical electromotive force response according to the adjusted initial value of the apparent resistivity, until the initial value of each time window does not need to be adjusted, and obtaining the electromagnetic full-range apparent resistivity.

Description

TECHNICAL FIELD

[0001] The present disclosure relates to the field of exploration geophysics, in particular, to a method, device, medium and equipment for calculating electromagnetic full-time apparent resistivity. BACKGROUND

[0002] The electric source transient electromagnetic method uses a grounding line source to emit a primary pulse magnetic field to the underground, and uses a coil to observe an induced electromotive force during the intermittent period of the primary pulse magnetic field, so as to infer the underground medium distribution according to the induced electromotive force. The calculation method of the apparent resistivity includes early, late and full-time three modes. Compared with the early or late apparent resistivity, the full-time apparent resistivity has the advantages of large time application range and closer to the true value.

[0003] In related scenarios, the full-time apparent resistivity is determined according to the vertical induced electromotive force or the vertical magnetic field data. First, the vertical induced electromotive force is not a single-valued function of the resistivity, and the apparent resistivity solved according to the vertical induced electromotive force is not unique. Although the vertical magnetic field is theoretically a single-valued function of the resistivity, in practice, the magnetic field is integrated from the induced electromotive force with time, and errors are introduced in the integration process. Moreover, in some scenarios, the early transient magnetic field will tend to be a direct current magnetic field, at which time the sensitivity of the magnetic field to the resistivity is weak, and combined with the errors produced in the numerical integration process, the solved apparent resistivity greatly deviates from the true value. SUMMARY

[0004] The purpose of the present disclosure is to provide a method, device, medium and equipment for calculating electromagnetic full-time apparent resistivity, so as to overcome the technical defects of the deviation of the apparent resistivity from the true value in related scenarios, thereby improving the accuracy of the full-time apparent resistivity calculation and improving the precision of the detection results of the electric source transient electromagnetic method.

[0005] To achieve the above-mentioned purpose, the first aspect of the present disclosure provides a method for calculating electromagnetic full-time apparent resistivity, comprising:

[0006] Obtaining the vertical induced electromotive force received by each time window electric dipole element, and determining the initial value of the apparent resistivity corresponding to each time window;

[0007] According to the measurement parameters of the electric dipole element, the initial value of the apparent resistivity corresponding to each time window and the vertical induced electromotive force corresponding to the time window, determining the transient vertical magnetic field response and the transient vertical electromotive force response of each time window;

[0008] According to the transient vertical electromotive force response of each time window, the corresponding vertical induced electromotive force and the preset error threshold, determining whether to adjust the initial value of the apparent resistivity corresponding to each time window;

[0009] In a case where the initial value of the apparent resistivity corresponding to any of the time windows is determined to be adjusted, an initial value adjustment strategy of the apparent resistivity corresponding to the time window is determined according to the transient vertical magnetic field response and the vertical induced magnetic field value of the time window, and the initial value of the apparent resistivity corresponding to the time window is adjusted, and the corresponding transient vertical magnetic field response and transient vertical electromotive force response are recalculated according to the adjusted initial value of the apparent resistivity, until the initial value of the apparent resistivity corresponding to each of the time windows does not need to be adjusted, and the electromagnetic full-range apparent resistivity is obtained.

[0010] Optionally, the measurement parameters include: a transmitting source length corresponding to the electric dipole element, an electric dipole element length after the grounding line source is divided into electric dipole elements, an observation distance between the electric dipole element and an observation point, an included angle between the electric dipole element and the observation point along the grounding line source, an attenuation delay between adjacent time windows, and a preset vacuum magnetic permeability.

[0011] Optionally, the determination of the transient vertical magnetic field response and the transient vertical electromotive force response of each of the time windows according to the measurement parameters of the electric dipole element, the initial value of the apparent resistivity corresponding to each of the time windows, and the vertical induced electromotive force corresponding to the time window includes:

[0012] determination of a calculation factor of each of the time windows according to the initial value of the apparent resistivity corresponding to each of the time windows, the preset vacuum magnetic permeability in the measurement parameters, the observation distance, and the attenuation delay;

[0013] substitution of the preset vacuum magnetic permeability, the electric dipole element length, the transmitting source length corresponding to the electric dipole element, the calculation factor corresponding to each of the time windows, the observation distance, and the included angle into a first calculation formula to obtain the transient vertical magnetic field response corresponding to each of the time windows;

[0014] substitution of the electric dipole element length, the transmitting source length corresponding to the electric dipole element, the calculation factor corresponding to each of the time windows, the observation distance, the initial value of the apparent resistivity, and the included angle into a second calculation formula to obtain the transient vertical electromotive force response corresponding to each of the time windows.

[0015] Optionally, the first calculation formula is:

[0016]

[0017] wherein Φ (u) and u are the calculation factor, r is the observation distance, μ0 is the preset vacuum magnetic permeability, is the included angle between the electric dipole element and the observation point along the grounding line source, and L x is the transmitting source length corresponding to the electric dipole element, and l is the electric dipole element length after the grounding line source is divided into electric dipole elements.

[0018] Optionally, the second calculation formula is:

[0019]

[0020] wherein, Φ(u) and u are the calculation factor, r is the observation distance, p i is an initial value of apparent resistivity corresponding to the i-th time window, is an angle between the electric dipole element and the observation point along the ground wire source, L x is a length of the transmission source corresponding to the electric dipole element, and l is a length of the electric dipole element after the ground wire source is divided into electric dipole elements.

[0021] Optionally, the determining whether to adjust the initial value of the apparent resistivity corresponding to each of the time windows according to the transient vertical electric motive force response, the corresponding vertical induced electric motive force, and a preset error threshold value of each of the time windows comprises:

[0022] calculating a ratio of the vertical induced electric motive force to the corresponding transient vertical electric motive force response of each of the time windows;

[0023] calculating an absolute value of a difference between the ratio and 1;

[0024] determining whether a size relationship between the absolute value and the preset error threshold value meets a preset condition;

[0025] determining whether to adjust the initial value of the apparent resistivity corresponding to each of the time windows according to whether the preset condition is met.

[0026] Optionally, the preset condition is that the absolute value is less than or equal to the preset error threshold value.

[0027] Optionally, in a case where it is determined to adjust the initial value of the apparent resistivity corresponding to any of the time windows, the initial value adjustment strategy of the apparent resistivity corresponding to the time window is determined according to the transient vertical magnetic field response and the vertical induced magnetic field value of the time window.

[0028] In a case where it is determined to adjust the initial value of the apparent resistivity corresponding to any of the time windows, if the transient vertical magnetic field response of the time window is greater than the vertical induced magnetic field value, the adjustment strategy is to reduce the initial value of the apparent resistivity corresponding to the time window, and if the transient vertical magnetic field response of the time window is less than the vertical induced magnetic field value, the adjustment strategy is to increase the initial value of the apparent resistivity corresponding to the time window.

[0029] Optionally, the determining the initial value of the apparent resistivity corresponding to each of the time windows comprises:

[0030] The equivalent area of the receiving coil of the electric dipole element, the decay delay time between adjacent time windows in the measurement parameter, the preset vacuum magnetic permeability, and the vertical induced electromotive force corresponding to a next time window of a current time window are substituted into a third calculation formula to obtain an initial value of the apparent resistivity corresponding to the current time window.

[0031] Optionally, the third calculation formula is:

[0032]

[0033] wherein S is the equivalent area of the receiving coil of the electric dipole element, t i is the decay delay time between adjacent time windows in the measurement parameter, μ0 is the preset vacuum magnetic permeability, is the vertical induced electromotive force corresponding to a next time window of a current time window, ρ i is the initial value of the apparent resistivity corresponding to the current time window.

[0034] In a second aspect of the present disclosure, an electromagnetic full-range apparent resistivity calculation device is provided, comprising:

[0035] an acquisition module configured to acquire the vertical induced electromotive force received by each time window electric dipole element, and determine the initial value of the apparent resistivity corresponding to each time window;

[0036] a first determination module configured to determine the transient vertical magnetic field response and the transient vertical electromotive force response of each time window according to the measurement parameter of the electric dipole element, the initial value of the apparent resistivity corresponding to each time window, and the vertical induced electromotive force corresponding to the time window;

[0037] a second determination module configured to determine whether to adjust the initial value of the apparent resistivity corresponding to each time window according to the transient vertical electromotive force response of each time window, the corresponding vertical induced electromotive force, and a preset error threshold;

[0038] an adjustment module configured to, in a case where it is determined to adjust the initial value of the apparent resistivity corresponding to any time window, determine an adjustment strategy for the initial value of the apparent resistivity corresponding to the time window according to the transient vertical magnetic field response and the vertical induced magnetic field value of the time window, and recalculate the corresponding transient vertical magnetic field response and the transient vertical electromotive force response according to the adjusted initial value of the apparent resistivity, until the initial value of the apparent resistivity corresponding to each time window does not need to be adjusted, to obtain the electromagnetic full-range apparent resistivity.

[0039] Optionally, the measurement parameters comprise: a length of a transmitting source corresponding to the electric dipole element, a length of an electric dipole element after the grounding line source is divided into the electric dipole element, an observation distance between the electric dipole element and an observation point, an included angle between the electric dipole element and the observation point along the grounding line source, an attenuation delay between adjacent time windows, and a preset vacuum permeability.

[0040] Optionally, the first determining module is configured to:

[0041] determine a calculation factor of each of the time windows according to an initial value of apparent resistivity corresponding to each of the time windows, the preset vacuum permeability in the measurement parameters, the observation distance, and the attenuation delay;

[0042] substitute the preset vacuum permeability, the length of the electric dipole element, the length of the transmitting source corresponding to the electric dipole element, the calculation factor corresponding to each of the time windows, the observation distance, and the included angle into a first calculation formula to obtain the transient vertical magnetic field response corresponding to each of the time windows;

[0043] substitute the length of the electric dipole element, the length of the transmitting source corresponding to the electric dipole element, the calculation factor corresponding to each of the time windows, the observation distance, the initial value of the apparent resistivity, and the included angle into a second calculation formula to obtain the transient vertical electromotive force response corresponding to each of the time windows.

[0044] Optionally, the first calculation formula is:

[0045]

[0046] wherein Φ (u) and u are the calculation factor, r is the observation distance, μ0 is the preset vacuum permeability, is an included angle between the electric dipole element and the observation point along the grounding line source, L x is the length of the transmitting source corresponding to the electric dipole element, and l is the length of the electric dipole element after the grounding line source is divided into the electric dipole element.

[0047] Optionally, the second calculation formula is:

[0048]

[0049] wherein Φ (u) and u are the calculation factor, r is the observation distance, ρ i is an initial value of apparent resistivity corresponding to the i-th time window, is an included angle between the electric dipole element and the observation point along the grounding line source, L x is the length of the transmitting source corresponding to the electric dipole element, and l is the length of the electric dipole element after the grounding line source is divided into the electric dipole element.

[0050] Optionally, the second determining module is configured to:

[0051] calculate a ratio of the vertical induced electromotive force of each of the time windows to the corresponding transient vertical electromotive force response;

[0052] calculate an absolute value of a difference between the ratio and 1;

[0053] determine whether a size relationship between the absolute value and the preset error threshold satisfies a preset condition;

[0054] determine whether to adjust the initial value of the apparent resistivity corresponding to each of the time windows according to whether the preset condition is satisfied.

[0055] Optionally, the preset condition is that the absolute value is less than or equal to the preset error threshold.

[0056] Optionally, the adjusting module is configured to:

[0057] in a case where it is determined to adjust the initial value of the apparent resistivity corresponding to any of the time windows, if the transient vertical magnetic field response of the time window is greater than the vertical induced magnetic field value, the adjustment strategy is to reduce the initial value of the apparent resistivity corresponding to the time window, and if the transient vertical magnetic field response of the time window is less than the vertical induced magnetic field value, the adjustment strategy is to increase the initial value of the apparent resistivity corresponding to the time window.

[0058] Optionally, the obtaining module is configured to:

[0059] substitute the equivalent area of the receiving coil of the electric dipole element, the decay delay time between adjacent time windows in the measurement parameter, the preset vacuum magnetic permeability, and the vertical induced electromotive force corresponding to a time window subsequent to the current time window into a third calculation formula to obtain the initial value of the apparent resistivity corresponding to the current time window.

[0060] Optionally, the third calculation formula is:

[0061]

[0062] wherein S is the equivalent area of the receiving coil of the electric dipole element, t i is the decay delay time between adjacent time windows in the measurement parameter, μ0 is the preset vacuum magnetic permeability, is the vertical induced electromotive force corresponding to a time window subsequent to the current time window, ρ i is the initial value of the apparent resistivity corresponding to the current time window.

[0063] A third aspect of this disclosure provides a non-transitory computer-readable storage medium having a computer program stored thereon that, when executed by a processor, implements the steps of the method described in any of the first aspects.

[0064] A fourth aspect of this disclosure provides an electronic device, comprising:

[0065] A memory on which computer programs are stored;

[0066] A processor for executing the computer program in the memory to implement the steps of the method of any one of the first aspects.

[0067] The above technical solution can achieve at least the following beneficial effects:

[0068] By employing the transient vertical electromotive force response and the corresponding vertical induced electromotive force, it is determined whether the initial value of the apparent resistivity corresponding to the time window needs to be adjusted. If the initial value of the apparent resistivity corresponding to the time window needs to be adjusted, the adjustment strategy for the initial value of the apparent resistivity corresponding to the time window is determined based on the transient vertical magnetic field response and the vertical induced magnetic field value of that time window. In this way, the combined approach of vertical magnetic field and vertical induced electromotive force can improve the accuracy of the apparent resistivity calculation throughout the entire process and improve the accuracy of the detection results of the transient electromagnetic method of the electrical source.

[0069] Other features and advantages of this disclosure will be described in detail in the following detailed description section. Attached Figure Description

[0070] The accompanying drawings are provided to further illustrate the present disclosure and form part of the specification. They are used together with the following detailed description to explain the present disclosure, but do not constitute a limitation thereof. In the drawings:

[0071] Figure 1 This is a flowchart of a method for calculating the apparent resistivity of an electromagnetic circuit according to an embodiment.

[0072] Figure 2 This is a schematic diagram of a grounding source arrangement according to an embodiment.

[0073] Figure 3 This is an implementation provided according to one embodiment. Figure 1 The flowchart of step S2.

[0074] Figure 4 This is an implementation provided according to one embodiment. Figure 1 The flowchart of step S3.

[0075] Figure 5 This is a block diagram of an electromagnetic apparent resistivity calculation device according to an embodiment.

[0076] Figure 6 is a block diagram of an electronic device according to an embodiment. DETAILED DESCRIPTION

[0077] The specific embodiments of the present disclosure will be described in detail below with reference to the accompanying drawings. It should be understood that the specific embodiments described herein are merely intended to illustrate and explain the present disclosure, and are not intended to limit the present disclosure.

[0078] To overcome the technical problem that the full-range apparent resistivity calculated according to the vertical induced electromotive force or the vertical magnetic field data deviates greatly from the true value in the related scenario, an electromagnetic full-range apparent resistivity calculation method is provided in the embodiments of the present disclosure, as shown in FIG. 1, the method comprises the following steps. Figure 1

[0079] In step S1, the vertical induced electromotive force received by each time window electric dipole element is obtained, and the initial value of the apparent resistivity corresponding to each time window is determined.

[0080] It can be explained that the electric dipole element in the embodiments of the present disclosure is only one electric dipole device, instead of the rectangular electric dipole device with four rods in the prior art.

[0081] wherein the vertical induced electromotive force is the normalized induced electromotive force, for example, the normalized vertical induced electromotive force V can be calculated by the following formula z 0

[0082]

[0083] wherein U is the initial vertical induced electromotive force collected, I is the transmitting current, and s is the receiving coil area.

[0084] wherein the initial value of the apparent resistivity corresponding to each time window is determined according to the initial value of the apparent resistivity of the next time window of the current time window based on the reverse order method. For example, the apparent resistivity p of the i+1th time window is the initial value of the apparent resistivity of the ith time window. That is, in the embodiments of the present disclosure, the apparent resistivity of the last time window is calculated first, then the apparent resistivity of the second last time window is calculated, and so on, and finally the apparent resistivity of the first time window is calculated. i+1

[0085] In step S2, the transient vertical magnetic field response and the transient vertical electromotive force response of each time window are determined according to the measurement parameters of the electric dipole element, the initial value of the apparent resistivity corresponding to each time window, and the vertical induced electromotive force of the corresponding time window. ​​​

[0086] In the embodiments of the present disclosure, the measured parameters of the electric dipole element, the initial values of the apparent resistivity corresponding to each time window, and the vertical induced electromotive force corresponding to the time window are directly substituted into the pre-designed calculation formula to obtain the transient vertical magnetic field response and the transient vertical electromotive force response corresponding to each time window.

[0087] In step S3, whether to adjust the initial value of the apparent resistivity corresponding to each time window is determined according to the transient vertical electromotive force response, the corresponding vertical induced electromotive force, and the preset error threshold.

[0088] In a case where it is determined that the initial value of the apparent resistivity corresponding to each time window is not adjusted, the initial value of the apparent resistivity is taken as the final value of the apparent resistivity of the current time window, and the value of the apparent resistivity is directly taken as the initial value of the apparent resistivity of the previous time window.

[0089] In step S4, in a case where it is determined that the initial value of the apparent resistivity corresponding to any time window is adjusted, the initial value adjustment strategy of the apparent resistivity corresponding to the time window is determined according to the transient vertical magnetic field response and the vertical induced magnetic field value of the time window, and the corresponding transient vertical magnetic field response and transient vertical electromotive force response are recalculated according to the adjusted initial value of the apparent resistivity until the initial value of the apparent resistivity corresponding to each time window does not need to be adjusted, to obtain the electromagnetic full-path apparent resistivity.

[0090] The vertical induced magnetic field value of the time window is calculated according to the transient vertical magnetic field response of the time window. The vertical induced magnetic field value of the time window is calculated according to the transient vertical magnetic field response of the time window.

[0091] In the embodiments of the present disclosure, the step of step S2 is performed according to the adjusted initial value of the apparent resistivity, the corresponding transient vertical magnetic field response and transient vertical electromotive force response are recalculated until the initial value of the apparent resistivity corresponding to each time window does not need to be adjusted, and the electromagnetic full-path apparent resistivity is generated according to the apparent resistivity corresponding to each time window.

[0092] The above technical solution adopts the transient vertical electromotive force response and the corresponding vertical induced electromotive force to determine whether the initial value of the apparent resistivity corresponding to the time window needs to be adjusted, and in a case where the initial value of the apparent resistivity corresponding to the time window needs to be adjusted, the initial value adjustment strategy of the apparent resistivity corresponding to the time window is determined according to the transient vertical magnetic field response and the vertical induced magnetic field value of the time window. In this way, the vertical magnetic field and the vertical induced electromotive force are jointly used, which can improve the accuracy of the full-path apparent resistivity calculation and improve the precision of the detection result of the electrical source transient electromagnetic method.

[0093] Optionally, referring to Figure 2The measurement parameters include: a launch source length corresponding to the electric dipole element, an electric dipole element length after the grounding line source is divided into electric dipole elements, an observation distance between the electric dipole element and an observation point, an included angle between the electric dipole element and the observation point along the grounding line source, an attenuation delay between adjacent time windows, and a preset vacuum permeability.

[0094] The observation distance r is a distance from a midpoint of the actual electric dipole element to the observation point, and the included angle between the electric dipole element and the observation point along the grounding line source is an included angle between a line connecting the midpoint of the electric dipole element and the observation point and the X axis.

[0095] Optionally, referring to Figure 3 In step S2, the transient vertical magnetic field response and the transient vertical electromotive force response of each time window are determined according to the measurement parameters of the electric dipole element, the initial value of the apparent resistivity corresponding to each time window, and the vertical induced electromotive force of the corresponding time window, and the determination includes:

[0096] In step S21, the calculation factor of each time window is determined according to the initial value of the apparent resistivity corresponding to each time window, the preset vacuum permeability in the measurement parameters, the observation distance, and the attenuation delay.

[0097] In the embodiments of the present disclosure, the calculation factor includes a first calculation sub-factor and a second calculation sub-factor, wherein the initial value of the apparent resistivity corresponding to the time window, the preset vacuum permeability in the measurement parameters, the observation distance, and the attenuation delay are substituted into a first factor calculation formula to obtain the first calculation sub-factor, and the first factor calculation formula can be:

[0098]

[0099] Wherein, r is the observation distance, μ0 is the preset vacuum permeability, ρ i is the initial value of the apparent resistivity corresponding to the i-th time window, t i is the attenuation delay between adjacent time windows in the measurement parameters.

[0100] Further, the first calculation sub-factor is substituted into a second factor calculation formula to obtain the second calculation sub-factor, and the second factor calculation formula can be:

[0101]

[0102] In step S22, the preset vacuum permeability, the electric dipole element length, the launch source length corresponding to the electric dipole element, the calculation factor corresponding to each time window, the observation distance, and the included angle are substituted into a first calculation formula to obtain the transient vertical magnetic field response corresponding to each time window.

[0103] The first calculation formula is:

[0104]

[0105] Where Φ(u) and u are calculation factors, r is the observation distance, and μ0 is the preset vacuum permeability. L is the angle between the electric dipole element and the observation point along the grounding source. x denoted as the emitter length corresponding to the electric dipole element, and l is the length of the electric dipole element after the grounding wire source is divided into electric dipole elements.

[0106] In step S23, the length of the electric dipole element, the length of the emitter corresponding to the electric dipole element, the calculation factor corresponding to each time window, the observation distance, the initial value of the apparent resistivity, and the included angle are substituted into the second calculation formula to obtain the transient vertical electromotive force response corresponding to each time window.

[0107] The second calculation formula is:

[0108]

[0109] Where Φ(u) and u are calculation factors, r is the observation distance, and ρ is the distance between the observation and ... i Let be the initial value of the apparent resistivity corresponding to the i-th time window. L is the angle between the electric dipole element and the observation point along the grounding source. x denoted as the emitter length corresponding to the electric dipole element, and l is the length of the electric dipole element after the grounding wire source is divided into electric dipole elements.

[0110] Optionally, see Figure 4 As shown, in step S3, based on the transient vertical electromotive force response of each time window, the corresponding vertical induced electromotive force, and the preset error threshold, it is determined whether to adjust the initial value of the apparent resistivity corresponding to each time window, including:

[0111] In step S31, the ratio of the vertical induced electromotive force to the corresponding transient vertical electromotive force response for each time window is calculated.

[0112] In step S32, the absolute value of the difference between the ratio and 1 is calculated.

[0113] In step S33, it is determined whether the relationship between the absolute value and the preset error threshold satisfies the preset condition.

[0114] The following analytical formula is used to determine whether the relationship between the absolute value and the preset error threshold meets the preset conditions: Where Etol is the preset error threshold. In cases where the relationship between the absolute value and the preset error threshold does not meet the preset conditions, In this case, the relationship between the absolute value and the preset error threshold is determined to satisfy the preset conditions.

[0115] In step S34, it is determined whether to adjust the initial value of the apparent resistivity corresponding to each time window according to whether the preset condition is met.

[0116] In the embodiments of the present disclosure, in the case where the size relationship between the absolute value and the preset error threshold does not meet the preset condition, it is determined that the initial value of the apparent resistivity corresponding to each time window needs to be adjusted, and in the case where the size relationship between the absolute value and the preset error threshold meets the preset condition, it is determined that the initial value of the apparent resistivity corresponding to each time window does not need to be adjusted.

[0117] Optionally, the preset condition is that the absolute value is less than or equal to the preset error threshold.

[0118] Optionally, in the case where it is determined to adjust the initial value of the apparent resistivity corresponding to any time window, an initial value adjustment strategy of the apparent resistivity corresponding to the time window is determined according to the transient vertical magnetic field response and the vertical induced magnetic field value of the time window, and the initial value adjustment strategy includes:

[0119] In the case where it is determined to adjust the initial value of the apparent resistivity corresponding to any time window, if the transient vertical magnetic field response of the time window is greater than the vertical induced magnetic field value, the adjustment strategy is to reduce the initial value of the apparent resistivity corresponding to the time window, and if the transient vertical magnetic field response of the time window is less than the vertical induced magnetic field value, the adjustment strategy is to increase the initial value of the apparent resistivity corresponding to the time window.

[0120] For example, if and , the apparent resistivity is reduced i , and V z is recalculated until , the apparent resistivity is output i , if and , the apparent resistivity is increased i , and V z is recalculated until , the apparent resistivity is output i ; and in other cases, the apparent resistivity is directly output i .

[0121] In the embodiments of the present disclosure, in the case where it is determined to adjust the initial value of the apparent resistivity corresponding to any time window, the initial value of the apparent resistivity can be determined in a manner of preset gradient.

[0122] Optionally, in step S1, the initial value of the apparent resistivity corresponding to each time window is determined, including:

[0123] The equivalent area of the receiving coil of the electric dipole element, the decay delay time between adjacent time windows in the measurement parameter, the preset vacuum permeability, and the vertical induced electromotive force corresponding to the next time window of the current time window are substituted into the third calculation formula to obtain the initial value of the apparent resistivity corresponding to the current time window.

[0124] Optionally, the third calculation formula is:

[0125]

[0126] wherein S is the equivalent area of the receiving coil of the electric dipole element, t i is the decay delay time between adjacent time windows in the measurement parameter, μ0 is the preset vacuum permeability, is the vertical induced electromotive force corresponding to the next time window of the current time window, and ρ i is the initial value of the apparent resistivity corresponding to the current time window.

[0127] Referring to Figure 5 The electromagnetic full-path apparent resistivity calculation device 500 comprises:

[0128] An acquisition module 510 is configured to acquire the vertical induced electromotive force received by each time window electric dipole element, and determine the initial value of the apparent resistivity corresponding to each time window;

[0129] A first determination module 520 is configured to determine the transient vertical magnetic field response and the transient vertical electromotive force response of each time window according to the measurement parameter of the electric dipole element, the initial value of the apparent resistivity corresponding to each time window, and the vertical induced electromotive force corresponding to the time window;

[0130] A second determination module 530 is configured to determine whether to adjust the initial value of the apparent resistivity corresponding to each time window according to the transient vertical electromotive force response of each time window, the corresponding vertical induced electromotive force, and a preset error threshold;

[0131] An adjustment module 540 is configured to, in a case where it is determined to adjust the initial value of the apparent resistivity corresponding to any time window, determine an adjustment strategy for the initial value of the apparent resistivity corresponding to the time window according to the transient vertical magnetic field response and the vertical induced magnetic field value of the time window, and recalculate the corresponding transient vertical magnetic field response and transient vertical electromotive force response according to the adjusted initial value of the apparent resistivity, until the initial value of the apparent resistivity corresponding to each time window does not need to be adjusted, to obtain the electromagnetic full-path apparent resistivity.

[0132] Optionally, the measurement parameters comprise: a length of a transmitting source corresponding to the electric dipole element, a length of an electric dipole element after the grounding line source is divided into the electric dipole element, an observation distance between the electric dipole element and an observation point, an included angle between the electric dipole element and the observation point along the grounding line source, an attenuation delay between adjacent time windows, and a preset vacuum permeability.

[0133] Optionally, the first determining module 520 is configured to:

[0134] determine a calculation factor of each time window according to an initial value of apparent resistivity corresponding to each time window, the preset vacuum permeability in the measurement parameters, the observation distance, and the attenuation delay;

[0135] substitute the preset vacuum permeability, the length of the electric dipole element, the length of the transmitting source corresponding to the electric dipole element, the calculation factor corresponding to each time window, the observation distance, and the included angle into a first calculation formula to obtain the transient vertical magnetic field response corresponding to each time window;

[0136] substitute the length of the electric dipole element, the length of the transmitting source corresponding to the electric dipole element, the calculation factor corresponding to each time window, the observation distance, the initial value of the apparent resistivity, and the included angle into a second calculation formula to obtain the transient vertical electromotive force response corresponding to each time window.

[0137] Optionally, the first calculation formula is:

[0138]

[0139] wherein Φ (u) and u are the calculation factor, r is the observation distance, μ0 is the preset vacuum permeability, is an included angle between the electric dipole element and the observation point along the grounding line source, L x is a length of a transmitting source corresponding to the electric dipole element, and l is a length of an electric dipole element after the grounding line source is divided into the electric dipole element.

[0140] Optionally, the second calculation formula is:

[0141]

[0142] wherein Φ (u) and u are the calculation factor, r is the observation distance, ρ i is an initial value of apparent resistivity corresponding to the i-th time window, is an included angle between the electric dipole element and the observation point along the grounding line source, L x is a length of a transmitting source corresponding to the electric dipole element, and l is a length of an electric dipole element after the grounding line source is divided into the electric dipole element.

[0143] Optionally, the second determining module 530 is configured to:

[0144] calculate a ratio of the vertical induced electromotive force of each of the time windows to the corresponding transient vertical electromotive force response;

[0145] calculate an absolute value of a difference between the ratio and 1;

[0146] determine whether a size relationship between the absolute value and the preset error threshold meets a preset condition;

[0147] determine whether to adjust the initial value of the apparent resistivity corresponding to each of the time windows according to whether the preset condition is met.

[0148] Optionally, the preset condition is that the absolute value is less than or equal to the preset error threshold.

[0149] Optionally, the adjusting module 540 is configured to:

[0150] in a case where it is determined to adjust the initial value of the apparent resistivity corresponding to any of the time windows, if the transient vertical magnetic field response of the time window is greater than the vertical induced magnetic field value, the adjustment strategy is to reduce the initial value of the apparent resistivity corresponding to the time window, and if the transient vertical magnetic field response of the time window is less than the vertical induced magnetic field value, the adjustment strategy is to increase the initial value of the apparent resistivity corresponding to the time window.

[0151] Optionally, the obtaining module 510 is configured to:

[0152] substitute the equivalent area of the receiving coil of the electric dipole element, the decay delay time between adjacent time windows in the measurement parameter, the preset vacuum magnetic permeability, and the vertical induced electromotive force corresponding to a next time window of a current time window into a third calculation formula to obtain the initial value of the apparent resistivity corresponding to the current time window.

[0153] Optionally, the third calculation formula is:

[0154]

[0155] wherein S is the equivalent area of the receiving coil of the electric dipole element, t i is the decay delay time between adjacent time windows in the measurement parameter, μ0 is the preset vacuum magnetic permeability, is the vertical induced electromotive force corresponding to a next time window of a current time window, ρ i is the initial value of the apparent resistivity corresponding to the current time window.

[0156] With regard to the apparatus in the above-described embodiments, the specific manner in which the respective modules perform operations has been described in detail in the embodiments related to the method, and will not be described in detail here.

[0157] The embodiments of the present disclosure further provide a non-transitory computer readable storage medium, having stored thereon a computer program, which, when executed by a processor, implements the steps of the electromagnetic full-range apparent resistivity calculation method according to any one of the preceding embodiments.

[0158] The embodiments of the present disclosure further provide an electronic device, comprising:

[0159] a memory having stored thereon a computer program;

[0160] a processor configured to execute the computer program in the memory to implement the steps of the electromagnetic full-range apparent resistivity calculation method according to any one of the preceding embodiments.

[0161] Figure 6 is a block diagram of an electronic device 700 according to an exemplary embodiment. As shown, the electronic device 700 can include a processor 701 and a memory 702. The electronic device 700 can also include one or more of a multimedia component 703, an input / output (I / O) interface 704, and a communication component 705. Figure 6

[0162] ​The processor 701 is configured to control overall operations of the electronic device 700 to complete all or part of the steps of the electromagnetic full-range apparent resistivity calculation method described above. The memory 702 is configured to store various types of data to support operations of the electronic device 700, which can include, for example, instructions for any application or method operating on the electronic device 700, and application-related data. The memory 702 can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as a static random access memory (SRAM), an electrically erasable programmable read-only memory (EEPROM), an erasable programmable read-only memory (EPROM), a programmable read-only memory (PROM), a read-only memory (ROM), a magnetic memory, a flash memory, a magnetic disk, or an optical disk. The multimedia component 703 can include a screen and an audio component. The screen can be, for example, a touch screen, and the audio component is configured to output and / or input audio signals. For example, the audio component can include a microphone configured to receive external audio signals. The received audio signals can be further stored in the memory 702 or transmitted through the communication component 705. The audio component also includes at least one speaker configured to output audio signals. The I / O interface 704 provides an interface between the processor 701 and other interface modules, which can be a keyboard, a mouse, a button, and the like. The buttons can be virtual buttons or physical buttons. The communication component 705 is configured to perform wired or wireless communication between the electronic device 700 and other devices. The wireless communication, such as Wi-Fi, Bluetooth, near field communication (NFC), 2G, 3G, 4G, NB-IOT, eMTC, or other 5G, and the like, or a combination of one or more of them, is not limited herein. Therefore, the corresponding communication component 705 can include a Wi-Fi module, a Bluetooth module, an NFC module, and the like.

[0163] In an exemplary embodiment, the electronic device 700 can be implemented by one or more Application Specific Integrated Circuits (ASICs), Digital Signal Processors (DSPs), Digital Signal Processing Devices (DSPDs), Programmable Logic Devices (PLDs), Field Programmable Gate Arrays (FPGAs), controllers, microcontrollers, microprocessors, or other electronic elements for performing the electromagnetic full-range apparent resistivity calculation method described above.

[0164] In another exemplary embodiment, a computer-readable storage medium including program instructions is also provided, which, when executed by a processor, implement the steps of the electromagnetic full-range apparent resistivity calculation method described above. For example, the computer-readable storage medium can be the memory 702 described above including program instructions, which can be executed by the processor 701 of the electronic device 700 to complete the electromagnetic full-range apparent resistivity calculation method described above.

[0165] The preferred embodiments of the present disclosure are described in detail above with reference to the accompanying drawings, but the present disclosure is not limited to the specific details in the above-described embodiments. Within the technical concept of the present disclosure, various simple modifications can be made to the technical solutions of the present disclosure, and these simple modifications all belong to the protection scope of the present disclosure.

[0166] In addition, it should be noted that each specific technical feature described in the above specific embodiments can be combined in any appropriate manner without contradiction, and in order to avoid unnecessary repetition, the present disclosure will not further describe various possible combinations.

[0167] Furthermore, any combination of the various different embodiments of the present disclosure can also be made, as long as it does not deviate from the idea of the present disclosure, and it should also be considered as disclosed by the present disclosure.

Claims

1. An electromagnetic fullbore resistivity calculation method, characterized in that, The method comprises: acquiring a vertical induced electromotive force received by each time window electric dipole element, and determining an initial value of apparent resistivity corresponding to each time window; determining a transient vertical magnetic field response and a transient vertical electromotive force response of each time window according to a measurement parameter of the electric dipole element, the initial value of apparent resistivity corresponding to each time window, and the vertical induced electromotive force corresponding to the time window; determining whether to adjust the initial value of apparent resistivity corresponding to each time window according to the transient vertical electromotive force response of each time window, the corresponding vertical induced electromotive force, and a preset error threshold value; in a case where it is determined to adjust the initial value of apparent resistivity corresponding to any time window, determining an initial value adjustment strategy of apparent resistivity corresponding to the time window according to the transient vertical magnetic field response and the vertical induced magnetic field value of the time window, and recalculating the corresponding transient vertical magnetic field response and transient vertical electromotive force response according to the adjusted initial value of apparent resistivity until the initial value of apparent resistivity corresponding to each time window does not need to be adjusted, to obtain the electromagnetic full-path apparent resistivity.

2. The method of claim 1, wherein, The measurement parameter comprises: a transmission source length corresponding to the electric dipole element, an electric dipole element length after the grounding line source is divided into electric dipole elements, an observation distance between the electric dipole element and an observation point, an included angle between the electric dipole element and the observation point along the grounding line source, an attenuation delay between adjacent time windows, and a preset vacuum permeability.

3. The method of claim 2, wherein, The determination of the transient vertical magnetic field response and the transient vertical electromotive force response of each time window according to the measurement parameter of the electric dipole element, the initial value of apparent resistivity corresponding to each time window, and the vertical induced electromotive force corresponding to the time window comprises: determining a calculation factor of each time window according to the initial value of apparent resistivity corresponding to each time window, the preset vacuum permeability in the measurement parameter, the observation distance, and the attenuation delay; substituting the preset vacuum permeability, the electric dipole element length, the transmission source length corresponding to the electric dipole element, the calculation factor corresponding to each time window, the observation distance, and the included angle into a first calculation formula to obtain the transient vertical magnetic field response corresponding to each time window; substituting the electric dipole element length, the transmission source length corresponding to the electric dipole element, the calculation factor corresponding to each time window, the observation distance, the initial value of apparent resistivity, and the included angle into a second calculation formula to obtain the transient vertical electromotive force response corresponding to each time window.

4. The method of claim 3, wherein, The first calculation formula is: wherein, and is the calculation factor, r is the observation distance, is the preset vacuum permeability, is the angle between the electric dipole element and the observation point along the ground line source, is the corresponding transmitting source length of the electric dipole element, is the electric dipole element length after the ground line source is divided into electric dipole elements, is the transient vertical magnetic field response corresponding to the i-th time window, t i is the decay delay between adjacent time windows, is the initial value of the apparent resistivity corresponding to the i-th time window.

5. The method of claim 3, wherein, The second calculation formula is: wherein, and is the calculation factor, r is the observation distance, is the initial value of the apparent resistivity corresponding to the i-th time window, is the angle between the electric dipole element and the observation point along the ground line source, is the length of the transmission source corresponding to the electric dipole element, is the length of the electric dipole element after the ground line source is divided into electric dipole elements, is the transient vertical electromotive force response corresponding to the i-th time window, t i is the decay delay time between adjacent time windows.

6. The method according to any one of claims 1-5, characterized in that, The determination of whether to adjust the initial value of apparent resistivity corresponding to each time window according to the transient vertical electromotive force response of each time window, the corresponding vertical induced electromotive force, and the preset error threshold value comprises: calculating a ratio of the vertical induced electromotive force to the corresponding transient vertical electromotive force response of each time window; calculating an absolute value of a difference between the ratio and 1; determining whether a size relationship between the absolute value and the preset error threshold value meets a preset condition. Determine whether to adjust the initial value of the apparent resistivity corresponding to each time window according to whether the preset condition is met.

7. The method of claim 6, wherein, The preset condition is that the absolute value is less than or equal to the preset error threshold.

8. The method according to any one of claims 1-5, characterized in that, In a case where it is determined to adjust the initial value of the apparent resistivity corresponding to any time window, the initial value adjustment strategy for the apparent resistivity corresponding to the time window is determined according to the transient vertical magnetic field response and the vertical induced magnetic field value of the time window, including: In a case where it is determined to adjust the initial value of the apparent resistivity corresponding to any time window, if the transient vertical magnetic field response of the time window is greater than the vertical induced magnetic field value, the adjustment strategy is to reduce the initial value of the apparent resistivity corresponding to the time window, and if the transient vertical magnetic field response of the time window is less than the vertical induced magnetic field value, the adjustment strategy is to increase the initial value of the apparent resistivity corresponding to the time window.

9. The method according to any one of claims 1-5, characterized in that, The initial value of the apparent resistivity corresponding to each time window is determined, including: The equivalent area of the receiving coil of the electric dipole element, the decay delay time between adjacent time windows in the measurement parameter, the preset vacuum permeability, and the vertical induced electromotive force corresponding to the next time window of the current time window are substituted into a third calculation formula to obtain the initial value of the apparent resistivity corresponding to the current time window.

10. The method of claim 9, wherein, The third calculation formula is: wherein, is an equivalent area of a receiving coil of the electric dipole element, is an attenuation delay between adjacent time windows in the measurement parameter, is the preset vacuum permeability, is a vertical induced electromotive force corresponding to a next time window of the current time window, is an initial value of the apparent resistivity corresponding to the current time window.

11. An electromagnetic fullbore resistivity calculation device, characterised in that, including: The acquisition module is configured to acquire the vertical induced electromotive force received by each time window electric dipole element, and determine the initial value of the apparent resistivity corresponding to each time window; The first determination module is configured to determine the transient vertical magnetic field response and the transient vertical electromotive force response of each time window according to the measurement parameter of the electric dipole element, the initial value of the apparent resistivity corresponding to each time window, and the vertical induced electromotive force corresponding to the time window; The second determination module is configured to determine whether to adjust the initial value of the apparent resistivity corresponding to each time window according to the transient vertical electromotive force response of each time window, the corresponding vertical induced electromotive force, and the preset error threshold; The adjustment module is configured to, in a case where it is determined to adjust the initial value of the apparent resistivity corresponding to any time window, determine the initial value adjustment strategy for the apparent resistivity corresponding to the time window according to the transient vertical magnetic field response and the vertical induced magnetic field value of the time window, and recalculate the corresponding transient vertical magnetic field response and transient vertical electromotive force response according to the adjusted initial value of the apparent resistivity, until the initial value of the apparent resistivity corresponding to each time window does not need to be adjusted, to obtain the electromagnetic full-range apparent resistivity.

12. A non-transitory computer-readable storage medium having stored thereon a computer program, characterized in that, The program is executed by the processor to implement the steps of the method of any one of claims 1-10.

13. An electronic device, comprising: including: The memory has a computer program stored thereon; The processor is configured to execute the computer program in the memory to implement the steps of the method of any one of claims 1-10.

Images