Method and device for determining full-range apparent resistivity

Abstract

The application provides a full-range apparent resistivity determination method and device, and belongs to the technical field of exploration geophysics, and comprises the following steps: obtaining a criterion for full-range apparent resistivity; determining calculation basis of the full-range apparent resistivity according to the criterion; the calculation basis comprises vertical induced electromotive force and a vertical magnetic field; and determining target full-range apparent resistivity based on the calculation basis. The full-range apparent resistivity determination method and device provided by the application use the criterion value to select appropriate calculation basis, improve the accuracy of full-range apparent resistivity calculation, and further improve the accuracy of ground transient electromagnetic method detection results.

Description

Technical Field

[0001] This invention relates to the field of exploration geophysics, and in particular to a method and apparatus for determining apparent resistivity over a complete range of conditions. Background Technology

[0002] Apparent resistivity remains one of the most commonly used and fundamental parameters in transient electromagnetic method interpretation. Its calculation methods are divided into three types: early-term, late-term, and full-term. Compared to early-term or late-term apparent resistivity, full-term apparent resistivity, also known as all-time apparent resistivity or all-period apparent resistivity, has advantages such as a wider applicable time range and closer approximation to the true value, and is therefore widely used.

[0003] Currently, the calculation of apparent resistivity over the entire range using the ground transient electromagnetic method is usually based on data from vertical induced electromotive force or vertical magnetic field.

[0004] However, the above methods all have limitations, resulting in inaccurate calculations of apparent resistivity. Summary of the Invention

[0005] The present invention provides a method, apparatus, electronic device, storage medium and product for determining apparent resistivity throughout the entire process, which solves the defects of existing technologies that lead to inaccurate calculated apparent resistivity, thereby improving the accuracy of the calculation of apparent resistivity throughout the entire process and thus improving the accuracy of the ground transient electromagnetic method detection results.

[0006] This invention provides a method for determining apparent resistivity over a complete range, comprising:

[0007] Criteria for obtaining apparent resistivity throughout the entire range;

[0008] Based on the aforementioned criteria, the calculation basis for the apparent resistivity over the entire range is determined; the calculation basis includes: vertical induced electromotive force and vertical magnetic field;

[0009] Based on the aforementioned calculations, the target apparent resistivity throughout the entire range is determined.

[0010] According to a method for determining apparent resistivity over a given period provided by the present invention, the criterion for obtaining the apparent resistivity over a given period includes:

[0011] The criteria for determining the second time window are based on the initial value of the second time window;

[0012] The initial value of the second time window is determined based on the apparent resistivity of the first time window over the entire duration.

[0013] According to a method for determining the apparent resistivity over a given period provided by the present invention, the step of determining the calculation basis for the apparent resistivity over a given period based on the criterion includes:

[0014] If the criterion for the second time window is not less than the criterion threshold, the calculation basis for determining the second time window is the vertical induced electromotive force.

[0015] If the criterion for the second time window is less than the criterion threshold, the calculation basis for determining the second time window is the vertical magnetic field.

[0016] According to a method for determining the apparent resistivity over the entire range provided by the present invention, when the calculation basis for the second time window is the vertical induced electromotive force, the step of determining the calculation basis for the apparent resistivity over the entire range based on the criterion includes:

[0017] Based on the initial value of the second time window, the vertical induced electromotive force of the second time window is determined;

[0018] Based on the first permissible error of the vertical induced electromotive force, the initial value of the second time window is adjusted to determine the apparent resistivity of the second time window throughout its entire duration.

[0019] The target total apparent resistivity is determined based on the total apparent resistivity of the second time window.

[0020] According to a method for determining apparent resistivity over a given time window provided by the present invention, when the calculation basis for the second time window is a vertical magnetic field, determining the target apparent resistivity over the given time window based on the calculation basis includes:

[0021] Based on the initial value of the second time window, the vertical magnetic field of the second time window is determined;

[0022] Based on the second allowable error of the vertical magnetic field, the initial value of the second time window is adjusted to determine the apparent resistivity of the second time window throughout its entire duration.

[0023] The target total apparent resistivity is determined based on the total apparent resistivity of the second time window.

[0024] According to a method for determining the apparent resistivity over a given time window provided by the present invention, determining the target apparent resistivity over a given time window based on the apparent resistivity over a given time window includes:

[0025] If the second time window is the last time window, the apparent resistivity of the entire second time window is determined to be the target apparent resistivity.

[0026] If the second time window is not the last time window, the apparent resistivity of the second time window is used as the initial value of the previous time window to obtain the apparent resistivity of the previous time window, until it is determined that the previous time window is the last time window, and the apparent resistivity of the previous time window is determined to be the target apparent resistivity.

[0027] The present invention also provides a device for determining apparent resistivity over a given period, comprising:

[0028] The acquisition module is used to obtain the criteria for the apparent resistivity throughout the entire process;

[0029] The first determining module is used to determine the calculation basis for the apparent resistivity over the entire range based on the criterion; the calculation basis includes: vertical induced electromotive force and vertical magnetic field;

[0030] The second determining module is used to determine the target apparent resistivity throughout the entire process based on the calculation basis.

[0031] The present invention also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the full-range apparent resistivity determination method as described above.

[0032] The present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the method for determining the apparent resistivity over the entire range as described above.

[0033] The present invention also provides a computer program product, including a computer program that, when executed by a processor, implements the method for determining the apparent resistivity over the entire range as described above.

[0034] The present invention provides a method, apparatus, electronic device, storage medium and product for determining apparent resistivity throughout the entire process. By using criterion values ​​to select appropriate calculation basis, the accuracy of apparent resistivity calculation throughout the entire process is improved, thereby improving the accuracy of ground transient electromagnetic method detection results. Attached Figure Description

[0035] To more clearly illustrate the technical solutions in this invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of this invention. For those skilled in the art, other drawings can be obtained from these drawings without creative effort.

[0036] Figure 1 This is one of the flowcharts illustrating the method for determining apparent resistivity throughout the entire process provided by the present invention;

[0037] Figure 2 This is an integral schematic diagram of the transmitter source circuit provided by the present invention;

[0038] Figure 3 This is the second flowchart of the method for determining apparent resistivity throughout the entire process provided by the present invention;

[0039] Figure 4 This is a schematic diagram of the structure of the apparent resistivity determination device provided by the present invention;

[0040] Figure 5 This is a schematic diagram of the structure of the electronic device provided by the present invention. Detailed Implementation

[0041] To make the objectives, technical solutions, and advantages of this invention clearer, the technical solutions of this invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of this invention. All other embodiments obtained by those skilled in the art based on the embodiments of this invention without creative effort are within the scope of protection of this invention.

[0042] In the description of this application, it should be understood that the terms "first," "second," etc., are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of technical features indicated. Therefore, a feature defined with "first," "second," etc., may explicitly or implicitly include one or more of that feature. In the description of this application, unless otherwise stated, "a plurality of" means two or more.

[0043] The transient electromagnetic method is a method that uses an ungrounded loop or a grounded source to emit a pulsed magnetic field into the ground, and then uses a coil to observe the induced electromotive force during the interval between the pulsed magnetic fields, thereby inferring the distribution of the underground medium.

[0044] The vertical induced electromotive force is not a single-valued function of resistivity, and the apparent resistivity obtained from it is not unique.

[0045] Although the vertical magnetic field is theoretically a single-valued function of resistivity, on the one hand, most instruments measure induced electromotive force (EMF) rather than magnetic field, and the magnetic field is usually derived by integrating the induced EMF over time, which introduces errors. On the other hand, in some cases, the magnetic field of early transient electromagnetic methods tends to be close to a DC magnetic field. In this case, the magnetic field is less sensitive to resistivity, and coupled with the errors generated during numerical integration, the apparent resistivity obtained by it may deviate greatly from the true value.

[0046] The following is combined with Figures 1 to 5 The present invention describes the method, apparatus, electronic device, storage medium, and product for determining apparent resistivity over a given period of time, as provided by embodiments of the present invention.

[0047] Figure 1 This is one of the flowcharts illustrating the method for determining apparent resistivity throughout the entire process provided by this invention, such as... Figure 1 As shown, including but not limited to the following steps:

[0048] First, in step S1, the criterion for the apparent resistivity over the entire range is obtained.

[0049] Based on transient electromagnetic theory, Z is used as the criterion for the apparent resistivity of the entire range, and the basis for calculating the apparent resistivity of the entire range using the ground transient electromagnetic method is selected.

[0050] The criterion is Z, and its calculation method is as follows:

[0051]

[0052] Where i is the sequence number of the current time window; t i T is the delay of the current time window after the current is completely turned off; o ρ is the current-off time in the transmitter; the initial value of ρ is set to the calculated apparent resistivity of the next time window; R is the minimum distance from the measuring point to the transmitter; μ0 is the vacuum permeability.

[0053] Figure 2 This is an integral schematic diagram of the transmitter source circuit provided by the present invention, as shown below. Figure 2 As shown, o is the electric dipole element; o' is the observation point; and r is the distance from the electric dipole element to the observation point. Let dl be the angle between dl and r.

[0054] I, II, III, and IV are the integral trajectories of the entire transmitter circuit. The transmitter circuit is the entire transmitter circuit, which is the combination of I, II, III, and IV. It needs to be divided into countless electric dipole elements. The current multiplied by dl is the electric dipole element, and o is the position of the electric dipole source.

[0055] If i is the last time window, the initial value of the apparent resistivity ρ over the entire time period is calculated using the following formula:

[0056]

[0057] Where s is the equivalent area of ​​the receiving coil; S is the equivalent area of ​​the transmitting coil; V i For t i The measured vertical induced electromotive force at any given time.

[0058] The time windows provided by this invention are calculated in reverse order. For example, if there are a total of 10 time windows, the 10th time window is the last time window, and the initial value is the apparent resistivity of the 10th time window.

[0059] Optionally, the criterion for obtaining the apparent resistivity over the entire range includes:

[0060] The criteria for determining the second time window are based on the initial value of the second time window;

[0061] The initial value of the second time window is determined based on the apparent resistivity of the first time window over the entire duration.

[0062] Since the time windows are calculated in reverse order, the first time window is the next time window after the second time window. For example, if the index of the first time window is i=3, then the index of the second time window is i=2.

[0063] Figure 3 This is the second flowchart illustrating the method for determining apparent resistivity throughout the entire process provided by this invention. Figure 3 As shown, the apparent resistivity of the entire first time window is ρ. i+1 Then the apparent resistivity of the entire second time window is ρ. i The apparent resistivity of the entire first time window is ρ. i+1 , as the initial value of the second window.

[0064] If the first time window is not the last time window, the initial value of the first time window can be obtained from the apparent resistivity of the previous time window.

[0065] Further, in step S2, the calculation basis for the apparent resistivity of the entire range is determined according to the criterion; the calculation basis includes: vertical induced electromotive force and vertical magnetic field.

[0066] The criterion serves as the basis for selecting the calculation of the apparent resistivity over the entire range.

[0067] Optionally, the step of determining the calculation basis for the apparent resistivity over the entire range based on the criterion includes:

[0068] If the criterion for the second time window is not less than the criterion threshold, the calculation basis for determining the second time window is the vertical induced electromotive force.

[0069] If the criterion for the second time window is less than the criterion threshold, the calculation basis for determining the second time window is the vertical magnetic field.

[0070] The criterion threshold can be flexibly selected according to actual needs, such as 2 or 2.4.

[0071] For example, with a criterion threshold of 2.4, if Z ≥ 2.4, then the vertical induced electromotive force is selected as the basis for calculating the apparent resistivity over the entire range; if Z < 2.4, then the vertical magnetic field is selected as the basis for calculating the apparent resistivity over the entire range.

[0072] Furthermore, in step S3, the target apparent resistivity is determined based on the calculation criteria.

[0073] The vertical induced electromotive force or vertical magnetic field is selected as the basis for calculating the apparent resistivity over the entire range, so as to calculate the final target apparent resistivity over the entire range.

[0074] Optionally, when the calculation basis for the second time window is the vertical induced electromotive force, the step of determining the calculation basis for the apparent resistivity over the entire range according to the criterion includes:

[0075] Based on the initial value of the second time window, the vertical induced electromotive force of the second time window is determined;

[0076] Based on the first permissible error of the vertical induced electromotive force, the initial value of the second time window is adjusted to determine the apparent resistivity of the second time window throughout its entire duration.

[0077] The target total apparent resistivity is determined based on the total apparent resistivity of the second time window.

[0078] First, let ρ = ρ i+1 Substitute it into the following formula to calculate the vertical induced electromotive force V:

[0079]

[0080]

[0081]

[0082]

[0083] Where L is the source circuit; dl is the length of the electric dipole element after dividing the source into electric dipole elements; r is the distance from the electric dipole element to the observation point; T o This refers to the current-off duration in the transmitter source; Let dl be the angle between dl and r.

[0084] If V <V i And |1-V / V i If |>Etol, then ρ increases at certain intervals until |1-V / V i |≤Etol;If V>V i And |1-V / V i If |>Etol, then ρ is decreased at certain intervals until |1-V / V i |≤Etol. Where Etol is the first tolerance.

[0085] The adjustment interval of ρ is not fixed and is related to the adjustment accuracy. The smaller Etol is, the higher the adjustment accuracy and the smaller the adjustment interval; the larger Etol is, the lower the adjustment accuracy and the larger the adjustment interval.

[0086] Let ρ i =ρ, where ρ i For t i Apparent resistivity throughout the entire time period.

[0087] Optionally, when the calculation basis for the second time window is the vertical magnetic field, determining the target's apparent resistivity throughout the entire time frame based on the calculation basis includes:

[0088] Based on the initial value of the second time window, the vertical magnetic field of the second time window is determined;

[0089] Based on the second allowable error of the vertical magnetic field, the initial value of the second time window is adjusted to determine the apparent resistivity of the second time window throughout its entire duration.

[0090] The target total apparent resistivity is determined based on the total apparent resistivity of the second time window.

[0091] The first tolerance can be the same as or different from the second tolerance.

[0092] For example, if the first tolerance and the second tolerance are the same, and we take the second tolerance as Etol, first, let ρ = ρ i+1 Substitute it into the following formula to calculate the vertical magnetic field B:

[0093]

[0094]

[0095] If B i And |1-B / B i If |>Etol, then ρ is decreased at certain intervals until |1-B / B i |≤Etol;If B>B i And |1-B / B i If |>Etol, then ρ increases at certain intervals until |1-B / B i |≤Etol. Where B i For t i The actual vertical magnetic field at any given time is usually obtained by integrating the measured vertical induced electromotive force over time.

[0096] Similarly, the adjustment interval of ρ is not fixed and is related to the adjustment accuracy. The smaller Etol is, the higher the adjustment accuracy and the smaller the adjustment interval; the larger Etol is, the lower the adjustment accuracy and the larger the adjustment interval.​

[0097] Let ρ i =ρ, where ρ i For t i Apparent resistivity throughout the entire time period.

[0098] Optionally, determining the target apparent resistivity based on the apparent resistivity over the second time window includes:

[0099] If the second time window is the last time window, the apparent resistivity of the entire second time window is determined to be the target apparent resistivity.

[0100] If the second time window is not the last time window, the apparent resistivity of the second time window is used as the initial value of the previous time window to obtain the apparent resistivity of the previous time window, until it is determined that the previous time window is the last time window, and the apparent resistivity of the previous time window is determined to be the target apparent resistivity.

[0101] like Figure 3 As shown, in step 1, the initial value ρ of any time window is obtained, and the criterion Z of the time window is calculated based on the initial value ρ.

[0102] Step 2: Select the calculation basis according to criterion Z:

[0103] When criterion Z < 2.4, the vertical magnetic field is selected as the basis for calculation, and ρ is calculated at t. i The vertical magnetic field response B at time t, if |B|≤|B i | and |1-B / B i If |>Etol, then ρ is decreased at certain intervals until |1-B / B i |≤Etol;If |B|>|B i | and |1-B / B i If |>Etol, then ρ increases at certain intervals until |1-B / B i |≤Etol, obtain the apparent resistivity over the entire time window.

[0104] Under the criterion Z≥2.4, the vertical induced electromotive force is selected as the calculation basis, and ρ is calculated at t. i The vertical induced electromotive force response V at time t, if |V|≤|V i | and |1-V / V i If |>Etol, then ρ increases at certain intervals until |1-V / V i |≤Etol;If |V|>|V i | and |1-V / V iIf |>Etol, then ρ is decreased at certain intervals until |1-V / V i |≤Etol, obtain the apparent resistivity over the entire time window.

[0105] Step 3: If the sequence number of the time window is 1, end the calculation; if the sequence number of the time window is not 1, use the apparent resistivity of the entire time window as the initial value of the previous time window and execute Step 4.

[0106] Step 4: Let i = i-1, and repeat steps 1 to 3 until i = 1.

[0107] The method for determining apparent resistivity throughout the entire range provided by this invention improves the accuracy of the calculation of apparent resistivity throughout the entire range by using a criterion value to select an appropriate calculation basis, thereby improving the accuracy of the ground transient electromagnetic method detection results.

[0108] The apparent resistivity determination device provided by the present invention is described below. The apparent resistivity determination device described below and the apparent resistivity determination method described above can be referred to in correspondence.

[0109] Figure 4 This is a schematic diagram of the structure of the apparent resistivity determination device provided by the present invention, as shown below. Figure 4 As shown, it includes:

[0110] The acquisition module 401 is used to obtain the criteria for the apparent resistivity throughout the entire process;

[0111] The first determining module 402 is used to determine the calculation basis of the apparent resistivity throughout the entire range according to the criterion; the calculation basis includes: vertical induced electromotive force and vertical magnetic field;

[0112] The second determining module 403 is used to determine the target apparent resistivity throughout the entire process based on the calculation basis.

[0113] During device operation, the acquisition module 401 acquires the criteria for the apparent resistivity throughout the entire range; the first determination module 402 determines the calculation basis for the apparent resistivity throughout the entire range based on the criteria; the calculation basis includes: vertical induced electromotive force and vertical magnetic field; the second determination module 403 determines the target apparent resistivity throughout the entire range based on the calculation basis.

[0114] The apparent resistivity determination device provided by this invention improves the accuracy of apparent resistivity calculation by selecting appropriate calculation basis using criterion values, thereby improving the accuracy of ground transient electromagnetic method detection results.

[0115] Figure 5 This is a schematic diagram of the structure of the electronic device provided by the present invention, such as... Figure 5As shown, the electronic device may include a processor 510, a communication interface 520, a memory 530, and a communication bus 540, wherein the processor 510, the communication interface 520, and the memory 530 communicate with each other via the communication bus 540. The processor 510 can call logical instructions in the memory 530 to execute a method for determining the apparent resistivity over the entire range. This method includes: obtaining a criterion for the apparent resistivity over the entire range; determining the calculation basis for the apparent resistivity over the entire range based on the criterion; the calculation basis includes: vertical induced electromotive force and vertical magnetic field; and determining the target apparent resistivity over the entire range based on the calculation basis.

[0116] Furthermore, the logical instructions in the aforementioned memory 530 can be implemented as software functional units and, when sold or used as independent products, can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, essentially, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0117] On the other hand, the present invention also provides a computer program product, which includes a computer program that can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer can execute the method for determining the apparent resistivity over the entire range provided by the above methods. The method includes: obtaining a criterion for the apparent resistivity over the entire range; determining the basis for calculating the apparent resistivity over the entire range based on the criterion; the basis for calculating the apparent resistivity over the entire range includes: vertical induced electromotive force and vertical magnetic field; and determining the target apparent resistivity over the entire range based on the basis for calculating the apparent resistivity over the entire range.

[0118] In another aspect, the present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon. When executed by a processor, the computer program implements a method for determining the apparent resistivity over the entire range provided by the methods described above. The method includes: obtaining a criterion for the apparent resistivity over the entire range; determining a calculation basis for the apparent resistivity over the entire range based on the criterion; the calculation basis including: vertical induced electromotive force and vertical magnetic field; and determining a target apparent resistivity over the entire range based on the calculation basis.

[0119] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs. Those skilled in the art can understand and implement this without any creative effort.

[0120] Through the above description of the embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus necessary general-purpose hardware platforms, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solutions, in essence or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute the methods described in the various embodiments or some parts of the embodiments.

[0121] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for determining apparent resistivity over a complete range, characterized in that, include: The criteria for determining the second time window are based on the initial value of the second time window; The initial value of the second time window is determined based on the apparent resistivity over the entire duration of the first time window; the criterion for the second time window is... Based on transient electromagnetic theory, the calculation method is as follows: ； Where i is the sequence number of the second time window; This is the delay of the second time window after the current is completely turned off; This refers to the current-off duration in the transmitter source; The initial value is set to the calculated apparent resistivity of the next time window; R is the minimum distance from the measuring point to the emission source; Permeability of free space; If the criterion for the second time window is not less than the criterion threshold, the calculation basis for determining the second time window is the vertical induced electromotive force. If the criterion for the second time window is less than the criterion threshold, the calculation basis for determining the second time window is the vertical magnetic field. Based on the aforementioned calculation, the target total apparent resistivity is determined, which includes: determining the target total apparent resistivity based on the total apparent resistivity of the second time window, which includes: when the second time window is the last time window, determining the total apparent resistivity of the second time window as the target total apparent resistivity; when the second time window is not the last time window, using the total apparent resistivity of the second time window as the initial value of the previous time window to obtain the total apparent resistivity of the previous time window, until the previous time window is the last time window, and then determining the total apparent resistivity of the previous time window as the target total apparent resistivity.

2. The method for determining apparent resistivity throughout the entire range according to claim 1, characterized in that, When the calculation of the second time window is based on the vertical induced electromotive force, the apparent resistivity of the entire second time window is determined based on the following process: Based on the initial value of the second time window, the vertical induced electromotive force of the second time window is determined; Based on the first permissible error of the vertical induced electromotive force, the initial value of the second time window is adjusted to determine the apparent resistivity of the second time window throughout its entire duration.

3. The method for determining apparent resistivity throughout the entire range according to claim 1, characterized in that, When the calculation of the second time window is based on the vertical magnetic field, the apparent resistivity of the second time window is determined based on the following process: Based on the initial value of the second time window, the vertical magnetic field of the second time window is determined; Based on the second permissible error of the vertical magnetic field, the initial value of the second time window is adjusted to determine the apparent resistivity of the second time window throughout its duration.

4. A device for determining apparent resistivity throughout a given process, characterized in that, include: The acquisition module is used to determine the criteria for the second time window based on the initial value of the second time window; The initial value of the second time window is determined based on the apparent resistivity over the entire duration of the first time window; the criterion for the second time window is... Based on transient electromagnetic theory, the calculation method is as follows: ； Where i is the sequence number of the second time window; This is the delay of the second time window after the current is completely turned off; This refers to the current-off duration in the transmitter source; The initial value is set to the calculated apparent resistivity of the next time window; R is the minimum distance from the measuring point to the emission source; Permeability of free space; The first determining module is used to determine that the calculation basis of the second time window is the vertical induced electromotive force, provided that the criterion of the second time window is not less than the criterion threshold. The first determining module is also used to determine that the calculation basis of the second time window is the vertical magnetic field when the criterion of the second time window is less than the criterion threshold. The second determining module is used to determine the target total apparent resistivity based on the calculation basis, which includes: determining the target total apparent resistivity based on the total apparent resistivity of the second time window, which includes: when the second time window is the last time window, determining the total apparent resistivity of the second time window as the target total apparent resistivity; when the second time window is not the last time window, using the total apparent resistivity of the second time window as the initial value of the previous time window to obtain the total apparent resistivity of the previous time window, until the previous time window is the last time window, and determining the total apparent resistivity of the previous time window as the target total apparent resistivity.

5. An electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the program, it implements the method for determining the apparent resistivity over the entire range as described in any one of claims 1-3.

6. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the method for determining the apparent resistivity over the entire range as described in any one of claims 1-3.

7. A computer program product, comprising a computer program, characterized in that, When the computer program is executed by the processor, it implements the method for determining the apparent resistivity over the entire range as described in any one of claims 1-3.

Images