Forward calculation method and device of push-the-dielectric scanning logging instrument

Abstract

This invention provides a forward modeling method and apparatus for a push-to-feed dielectric scanning logging tool. The method includes: acquiring a measurement equivalent calculation model of the push-to-feed dielectric scanning logging tool; determining the electric field spectral domain and magnetic field spectral domain when the signal receiver receives signals transmitted from different signal transmitters; determining the relationship between the electric field harmonic components and the harmonic number, and the relationship between the magnetic field harmonic components and the harmonic number, based on the electric field spectral domain and the magnetic field spectral domain; determining the electric field components when the signal receiver receives signals transmitted from different signal transmitters, based on the relationship between the electric field harmonic components and the harmonic number, a preset harmonic number interpolation interval, and a first preset harmonic number cutoff value; determining the magnetic field components when the signal receiver receives signals transmitted from different signal transmitters, based on the relationship between the magnetic field harmonic components and the harmonic number, a preset harmonic number interpolation interval, and a first preset harmonic number cutoff value; and determining the amplitude ratio and phase difference between the electric field components and the magnetic field components.

Description

Technical Field

[0001] This invention relates to the field of dielectric logging technology, and more particularly to a forward modeling method and apparatus for a push-type dielectric scanning logging tool. Background Technology

[0002] This section is intended to provide background or context for the embodiments of the invention set forth in the claims. The description herein is not an admission that it is prior art simply because it is included in this section.

[0003] The dielectric constant is a key parameter for distinguishing between oil and water. Furthermore, the dielectric constant of sediments containing natural gas hydrates is an important parameter in the exploration and development of hydrate resources. Dielectric scanning logging (DSL) is currently an important logging method for the exploration and development of unconventional oil and gas resources such as natural gas hydrates, tight oil and gas, heavy oil, and carbonate rocks. Push-type DSL logging tools use a push plate and utilize multi-frequency, multi-source measurements of formation information, making them the mainstream DSL logging instrument. To obtain the dielectric constant, the amplitude ratio and phase difference of the DSL logging must first be calculated through forward modeling. For push-type DSL logging tools, the push arm, antenna, wellbore, and formation constitute a complex three-dimensional geological model. Forward modeling using three-dimensional numerical methods such as finite element method and finite difference method is performed. However, existing three-dimensional numerical methods have low forward modeling efficiency and consume a lot of computer memory. Summary of the Invention

[0004] This invention proposes a forward modeling method for a push-type dielectric scanning logging tool, which can improve the forward modeling efficiency of the push-type dielectric scanning logging tool and reduce the computer memory usage, including:

[0005] Obtain the measurement equivalent calculation model of the push-type dielectric scanning logging tool; wherein, the measurement equivalent calculation model is obtained by equating the push arm of the push-type dielectric scanning logging tool to a columnar metal core, the antenna to a magnetic dipole, and the formation to a columnar formation;

[0006] Based on the position of the signal receiver in the equivalent calculation model, determine the electric field spectral domain and magnetic field spectral domain when the signal receiver receives signals transmitted from different signal transmitters.

[0007] Based on the electric field spectral domain and the magnetic field spectral domain, determine the relationship between the electric field harmonic components and the harmonic number when the signal receiver receives signals transmitted from different signal transmitters, and the relationship between the magnetic field harmonic components and the harmonic number.

[0008] Based on the relationship between electric field harmonic components and harmonic numbers, the interpolation interval of preset harmonic numbers, and the first preset harmonic number cutoff value, the electric field components of the signal receiver when receiving signals transmitted from different signal transmitters are determined; based on the relationship between magnetic field harmonic components and harmonic numbers, the interpolation interval of preset harmonic numbers, and the first preset harmonic number cutoff value, the magnetic field components of the signal receiver when receiving signals transmitted from different signal transmitters are determined.

[0009] Determine the amplitude ratio and phase difference of the electric field components when the signal receiver receives signals transmitted from different signal transmitters; determine the amplitude ratio and phase difference of the magnetic field components when the signal receiver receives signals transmitted from different signal transmitters.

[0010] This invention provides a forward modeling calculation device for a push-type dielectric scanning logging tool, which can improve the forward modeling calculation efficiency of the push-type dielectric scanning logging tool and reduce the computer memory usage, including:

[0011] The model acquisition module is used to acquire the measurement equivalent calculation model of the push-type dielectric scanning logging tool; wherein, the measurement equivalent calculation model is obtained by equating the push arm of the push-type dielectric scanning logging tool to a columnar metal core, the antenna to a magnetic dipole, and the formation to a columnar formation.

[0012] The first determining module is used to determine the electric field spectral domain field and magnetic field spectral domain field of the signal receiving end when receiving signals transmitted from different signal transmitting ends, based on the position of the signal receiving end in the measurement equivalent calculation model.

[0013] The second determining module is used to determine the relationship between the electric field harmonic components and the harmonic number, and the relationship between the magnetic field harmonic components and the harmonic number, when the signal receiver receives signals transmitted by different signal transmitters, based on the electric field spectral domain field and the magnetic field spectral domain field.

[0014] The third determining module is used to determine the electric field components of the signal receiver when receiving signals transmitted from different signal transmitters, based on the relationship between the electric field harmonic components and the harmonic number, the interpolation interval of the preset harmonic number, and the first preset harmonic number cutoff value; and to determine the magnetic field components of the signal receiver when receiving signals transmitted from different signal transmitters, based on the relationship between the magnetic field harmonic components and the harmonic number, the interpolation interval of the preset harmonic number, and the first preset harmonic number cutoff value.

[0015] The fourth determining module is used to determine the amplitude ratio and phase difference of the electric field components when the signal receiver receives signals transmitted from different signal transmitters; and to determine the amplitude ratio and phase difference of the magnetic field components when the signal receiver receives signals transmitted from different signal transmitters.

[0016] An embodiment of the present invention proposes a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements a forward modeling method for a push-type dielectric scanning logging tool.

[0017] This invention provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements a forward modeling method for a push-type dielectric scanning logging tool.

[0018] This invention provides a computer program product comprising a computer program that, when executed by a processor, implements a forward modeling method for a push-type dielectric scanning logging tool.

[0019] The forward modeling method and apparatus for the push-type dielectric scanning logging tool proposed in this invention can solve the problems of low efficiency and high memory consumption of existing three-dimensional numerical methods. This invention obtains a measurement equivalent calculation model for the push-type dielectric scanning logging tool. This model is obtained by equating the push arm of the push-type dielectric scanning logging tool to a columnar metal core, the antenna to a magnetic dipole, and the formation to a columnar formation. Based on the position of the signal receiving end in the measurement equivalent calculation model, the electric field and magnetic field spectral domains of the signal receiving end when receiving signals transmitted from different signal transmitting ends are determined. Based on the electric field and magnetic field spectral domains, the calculation of the signal receiving end's response to different signal transmissions is determined. The relationship between the electric field harmonic components and the harmonic number when the signal is transmitted is determined; the relationship between the magnetic field harmonic components and the harmonic number is also determined; based on the relationship between the electric field harmonic components and the harmonic number, the interpolation interval of the preset harmonic number, and the first preset harmonic number cutoff value, the electric field components of the signal receiver when receiving signals transmitted from different signal transmitters are determined; based on the relationship between the magnetic field harmonic components and the harmonic number, the interpolation interval of the preset harmonic number, and the first preset harmonic number cutoff value, the magnetic field components of the signal receiver when receiving signals transmitted from different signal transmitters are determined; the amplitude ratio and phase difference of the electric field components of the signal receiver when receiving signals transmitted from different signal transmitters are determined; the amplitude ratio and phase difference of the magnetic field components of the signal receiver when receiving signals transmitted from different signal transmitters are determined. This invention provides a fast forward modeling algorithm for amplitude ratio and phase difference response of dielectric scanning logging. It uses the measurement equivalent calculation model of the push-type dielectric scanning logging tool as the basic model for forward modeling and presets the interpolation interval of harmonic numbers during the forward modeling process, which greatly reduces the amount of calculation and can improve the forward modeling efficiency of the push-type dielectric scanning logging tool and reduce the computer memory occupied. Attached Figure Description

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

[0021] Figure 1 This is a flowchart illustrating the forward modeling method of the push-type dielectric scanning logging tool according to an embodiment of the present invention.

[0022] Figure 2 This is a specific example diagram of the forward modeling calculation method of the push-type dielectric scanning logging tool according to an embodiment of the present invention;

[0023] Figure 3 This is a specific example diagram of the forward modeling calculation method of the push-type dielectric scanning logging tool according to an embodiment of the present invention;

[0024] Figure 4 This is a specific example diagram of the forward modeling calculation method of the push-type dielectric scanning logging tool according to an embodiment of the present invention;

[0025] Figure 5 This is a specific example diagram of the forward modeling calculation method of the push-type dielectric scanning logging tool according to an embodiment of the present invention;

[0026] Figure 6 This is a specific example diagram of the forward modeling calculation method of the push-type dielectric scanning logging tool according to an embodiment of the present invention;

[0027] Figure 7 This is a specific example diagram of the forward modeling calculation method of the push-type dielectric scanning logging tool according to an embodiment of the present invention;

[0028] Figure 8 This is a specific example diagram of the forward modeling calculation method of the push-type dielectric scanning logging tool according to an embodiment of the present invention;

[0029] Figure 9 This is a specific example diagram of the forward modeling calculation method of the push-type dielectric scanning logging tool according to an embodiment of the present invention;

[0030] Figure 10 This is a schematic diagram of the forward modeling calculation device of the push-type dielectric scanning logging tool according to an embodiment of the present invention;

[0031] Figure 11 This is a schematic diagram of a computer device in an embodiment of the present invention. Detailed Implementation

[0032] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the embodiments of the present invention will be further described in detail below with reference to the accompanying drawings. Here, the illustrative embodiments of the present invention and their descriptions are used to explain the present invention, but are not intended to limit the present invention.

[0033] In this document, the term "and / or" merely describes a relationship, indicating that three relationships can exist. For example, A and / or B can represent three cases: A alone, A and B simultaneously, and B alone. Furthermore, the term "at least one" in this document means any combination of at least two of any one or more elements. For example, including at least one of A, B, and C can mean including any one or more elements selected from the set consisting of A, B, and C.

[0034] In the description of this specification, the terms "comprising," "including," "having," and "containing" are open-ended terms, meaning that they include but are not limited to. The terms "an embodiment," "a specific embodiment," "some embodiments," and "for example," etc., refer to specific features, structures, or characteristics described in connection with that embodiment or example that are included in at least one embodiment or example of this application. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, or characteristics described can be combined in any suitable manner in one or more embodiments or examples. The order of steps involved in the various embodiments is used to illustrate the implementation of this application, and the order of steps is not limited and can be adjusted appropriately as needed.

[0035] The principles and spirit of the present invention will be explained in detail below with reference to several representative embodiments.

[0036] Figure 1 This is a schematic flowchart illustrating the forward modeling method of the push-type dielectric scanning logging tool according to an embodiment of the present invention. Figure 1 As shown, the method includes:

[0037] Step 101: Obtain the measurement equivalent calculation model of the push-type dielectric scanning logging tool; wherein, the measurement equivalent calculation model is obtained by equating the push arm of the push-type dielectric scanning logging tool to a columnar metal core, the antenna to a magnetic dipole, and the formation to a columnar formation.

[0038] Step 102: Based on the position of the signal receiver in the equivalent calculation model, determine the electric field spectral domain and magnetic field spectral domain when the signal receiver receives signals transmitted from different signal transmitters.

[0039] Step 103: Based on the electric field spectral domain and the magnetic field spectral domain, determine the relationship between the electric field harmonic components and the harmonic number, and the relationship between the magnetic field harmonic components and the harmonic number when the signal receiver receives signals transmitted from different signal transmitters.

[0040] Step 104: Based on the relationship between electric field harmonic components and harmonic numbers, the interpolation interval of preset harmonic numbers, and the first preset harmonic number cutoff value, determine the electric field components of the signal receiver when receiving signals transmitted from different signal transmitters; based on the relationship between magnetic field harmonic components and harmonic numbers, the interpolation interval of preset harmonic numbers, and the first preset harmonic number cutoff value, determine the magnetic field components of the signal receiver when receiving signals transmitted from different signal transmitters.

[0041] Step 105: Determine the amplitude ratio and phase difference of the electric field components when the signal receiver receives signals transmitted from different signal transmitters; determine the amplitude ratio and phase difference of the magnetic field components when the signal receiver receives signals transmitted from different signal transmitters.

[0042] Depend on Figure 1 As shown in the flowchart, this embodiment of the invention obtains the measurement equivalent calculation model of the push-type dielectric scanning logging tool. The measurement equivalent calculation model is obtained by equating the push arm of the push-type dielectric scanning logging tool to a columnar metal core, the antenna to a magnetic dipole, and the formation to a columnar formation. Based on the position of the signal receiving end in the measurement equivalent calculation model, the electric field spectral domain and magnetic field spectral domain are determined when the signal receiving end receives signals transmitted from different signal transmitting ends. Based on the electric field spectral domain and magnetic field spectral domain, the relationship between the electric field harmonic components and harmonic numbers, and the magnetic field harmonic components are determined when the signal receiving end receives signals transmitted from different signal transmitting ends. The relationship between electric field harmonic components and harmonic numbers is determined based on the relationship between electric field harmonic components and harmonic numbers, the preset harmonic number interpolation interval, and the first preset harmonic number cutoff value. Similarly, the magnetic field components are determined based on the relationship between magnetic field harmonic components and harmonic numbers, the preset harmonic number interpolation interval, and the first preset harmonic number cutoff value. The amplitude ratio and phase difference of the electric field components when receiving signals from different signal transmitters are also determined. This invention provides a fast forward modeling algorithm for amplitude ratio and phase difference response in dielectric scanning logging. It utilizes the equivalent calculation model of a push-type dielectric scanning logging tool as the basic model for forward modeling and presets the harmonic number interpolation interval during the forward modeling process, significantly reducing the computational load and improving the forward modeling efficiency of the push-type dielectric scanning logging tool while reducing the amount of computer memory required.

[0043] To provide a clearer explanation of the forward modeling method for the aforementioned push-type dielectric scanning logging tool, a detailed explanation of each step is provided below.

[0044] Figure 2 This is a specific example diagram of the forward modeling calculation method of the push-type dielectric scanning logging tool according to an embodiment of the present invention.

[0045] Reference position of dielectric scanner during dielectric logging Figure 2 The dielectric scanner is installed in the formation wellbore, and the push arm (i.e., columnar metal core) of the dielectric scanner is located in the formation invasion zone. Considering the geological conditions such as the instrument, wellbore, invasion zone, and formation, a measurement equivalent calculation model of the push-type dielectric scanning logging tool is established. The measurement equivalent calculation model of the push-type dielectric scanning logging tool is obtained by equating the push arm of the push-type dielectric scanning logging tool to a columnar metal core, the antenna to a magnetic dipole, and the formation to a columnar formation.

[0046] In one embodiment of the present invention, the electric field spectral domain and magnetic field spectral domain of the signal receiver when receiving signals transmitted from different signal transmitters are determined based on the position of the signal receiver in the equivalent calculation model, including:

[0047] Based on the position of the signal receiver in the equivalent measurement model, the coupling relationship of the field at the radial layer interface of the equivalent measurement model is described by the generalized reflection array and transmission array, and the electric field spectral domain and magnetic field spectral domain are determined when the signal receiver receives signals transmitted from different signal transmitters.

[0048] In one embodiment of the present invention, before determining the relationship between the electric field harmonic components and the harmonic number when the signal receiver receives signals transmitted from different signal transmitters, the method includes:

[0049] Obtain the electric and magnetic field components of the Fourier series expansion of the signal receiver in the equivalent calculation model of the measurement.

[0050] Based on the electric and magnetic field spectral domains, the relationships between the electric field harmonic components and harmonic numbers, and between the magnetic field harmonic components and harmonic numbers, are determined when the signal receiver receives signals transmitted from different signal transmitters. These relationships include:

[0051] Based on the electric field spectral domain field, the magnetic field spectral domain field, and the electric and magnetic field components of the Fourier series expansion, the relationship between the electric field harmonic components and the harmonic number, and the relationship between the magnetic field harmonic components and the harmonic number are determined when the signal receiver receives signals transmitted from different signal transmitters.

[0052] In practical implementation, a pseudo-analytical solution for the measurement equivalent calculation model is derived, and a rapid calculation method for the dielectric scanning logging response of the columnar formation model is established. Based on the rotational symmetry of the measurement equivalent calculation model and the invariance of the vertical properties of the formation, the electric and magnetic field components are expanded in the form of Fourier series. Taking the signal transmitter (i.e., the transmitting coil) transmitting a signal in the z-direction and the signal receiver (i.e., the receiving coil) receiving a signal in the z-direction as an example, the electric and magnetic field components expanded in the Fourier series are expressed by the following formulas:

[0053]

[0054] Where i represents the imaginary unit; E zz H represents the electric field component. zz The subscript 'zz' indicates that the transmitting coil transmits in the z-direction and the receiving coil receives in the z-direction; 'v' represents the harmonic number; 'z' represents the ordinate of the receiving coil, with the transmitting coil defaulting to point 0; 'k' represents the magnetic field components. z Indicates the wave number in the z-direction; e zzv Represents the electric field spectral domain; h zzv Represents the magnetic field spectral domain; This represents the angle between the instrument and the wellbore axis, which can be taken as 0 in this invention; E xx Electric field components and H xx The magnetic field components can be derived from E zz Represents components and H zz The magnetic field components are obtained by combining the results.

[0055] In specific implementation, for the electric field spectral domain field e zzv and magnetic field spectral field h zzv There is an analytical solution in cylindrical coordinates. For the cylindrical multilayer model, the coupling relationship of the field at the radial layer interface of the measurement equivalent calculation model is described by the generalized reflection array and transmission array. The electric field spectral domain and magnetic field spectral domain are determined when the signal receiver receives signals transmitted by different signal transmitters.

[0056] Figure 3 This is a specific example diagram of the forward modeling calculation method of the push-type dielectric scanning logging tool according to an embodiment of the present invention.

[0057] In one embodiment of the present invention, before determining the electric field components of the signal receiver when receiving signals transmitted from different signal transmitters based on the relationship between the electric field harmonic components and the harmonic number, the interpolation interval of the preset harmonic number, and the first preset harmonic number cutoff value, and before determining the magnetic field components of the signal receiver when receiving signals transmitted from different signal transmitters based on the relationship between the magnetic field harmonic components and the harmonic number, the interpolation interval of the preset harmonic number, and the first preset harmonic number cutoff value, the method further includes:

[0058] Step 301: Determine the first magnetic field component when the signal receiver receives a signal transmitted by a signal transmitter based on the relationship between the magnetic field harmonic component and the harmonic number and the first preset harmonic number cutoff value.

[0059] Step 302: Based on the relationship between the magnetic field harmonic components and the harmonic number, and the second preset harmonic number cutoff value, determine the second magnetic field component when the signal receiver receives the signal transmitted by the signal transmitter; wherein, the second preset harmonic number cutoff value is greater than the first preset harmonic number cutoff value.

[0060] Step 303: Determine the error between the second magnetic field component and the first magnetic field component;

[0061] Step 304: Determine the interpolation interval of the preset harmonic number based on the error between the second magnetic field component and the first magnetic field component.

[0062] In one embodiment of the present invention, determining the interpolation interval of a preset harmonic number based on the error between the second magnetic field component and the first magnetic field component includes:

[0063] Based on the error between the second magnetic field component and the first magnetic field component, a local error function is determined; wherein, the local error function is used to represent the relationship between the error between the second magnetic field component and the first magnetic field component and the harmonic number;

[0064] Based on the local error function, determine the interpolation interval for the preset harmonic number;

[0065] The local error function is as follows:

[0066]

[0067]

[0068] Among them, Loc err (v) represents the local error function; e p (v) represents the error between the second magnetic field component and the first magnetic field component; v represents the harmonic number; N1 represents the first preset harmonic number cutoff value; N2 represents the second preset harmonic number cutoff value; H represents the first magnetic field component; zz The second magnetic field component is represented by τ, which is a preset constant; P represents the second magnetic field component. p (v) represents an nth-order interpolation polynomial; where P p (v)=a0+a1v+a2v 2 +...+a n v n a0, a1, a2, ..., a n All are constant coefficients of the nth-order interpolation polynomial; n is the order of the interpolation polynomial; p is the interpolation interval of the harmonic number.

[0069] Figure 4 , Figure 5 This is a specific example diagram of the forward modeling calculation method of the push-type dielectric scanning logging tool according to an embodiment of the present invention.

[0070] In practical implementation, regarding the summation of the electric and magnetic field components from the Fourier series expansion, since the instrument and wellbore axes are coplanar, let... The electric and magnetic field components of the Fourier series expansion can be simplified. Taking the magnetic field component as an example, the magnetic field component of the Fourier series expansion can be simplified to the following formula:

[0071]

[0072] Among them, H zz The subscript zz indicates that the transmitting coil transmits in the z-direction and the receiving coil receives in the z-direction; H zzv Let τ represent the harmonic components of the magnetic field at different harmonic numbers, where τ is a preset constant; v represents the harmonic number, and when v = 0, τ = 1, otherwise τ = 2; a preset harmonic number cutoff value N can be set in the above formula, then the magnetic field components are expressed as:

[0073]

[0074] When simulating different frequencies and different positions of the transmitting and receiving coils, the |Imag(H) of the magnetic field harmonic components zzv The Imag function, which calculates the imaginary part of a complex number, exhibits a decay law on an exponential scale. (See reference...) Figure 4 Part (1), r bh The radius of the cylindrical metal core is represented by r1 and r2, which represent the distances of the transmitting and receiving coils from the center point, respectively. (Refer to...) Figure 4 In part (2), the closer the transmitting coil and receiving coil are to the cylindrical metal core, the more significant the N value becomes; the higher the frequency, the more significant the N value becomes.

[0075] make:

[0076]

[0077] f v =P p (v)+e p (v);

[0078] Convert the harmonic components of the magnetic field into a logarithmic representation, f v It is an intermediate variable in the formula, representing the age-dependent form of the magnetic field harmonic components; P p (v) represents an nth-order interpolation polynomial; where P p (v)=a0+a1v+a2v 2 +...+a nv n a0, a1, a2, ..., a n All are constant coefficients of an nth-order interpolation polynomial; n is the order of the interpolation polynomial; p is the interpolation interval for harmonic numbers; e p (v) represents the error between the second magnetic field component and the first magnetic field component. The second magnetic field component is calculated under the condition of a second preset harmonic cutoff value. The larger the second preset harmonic cutoff value, the closer the calculated second magnetic field component is to the actual value. The first magnetic field component can be regarded as an approximation value. That is, when the error is ignored, the magnetic field component can be approximated by the first magnetic field component.

[0079]

[0080] The local error function is defined as follows:

[0081]

[0082]

[0083] Among them, Loc err (v) represents the local error function; e p (v) represents the error between the second magnetic field component and the first magnetic field component; v represents the harmonic number; N1 represents the first preset harmonic number cutoff value; N2 represents the second preset harmonic number cutoff value; H represents the first magnetic field component; zz The second magnetic field component is represented by τ, which is a preset constant; P represents the second magnetic field component. p (v) represents an nth-order interpolation polynomial; where P p (v)=a0+a1v+a2v 2 +...+a n v n a0, a1, a2, ..., a n All are constant coefficients of an nth-order interpolation polynomial; n is the order of the interpolation polynomial; p is the interpolation interval for harmonic numbers; Reference Figure 5 , representing local error (Loc err The relationship between the harmonic number and the interpolation interval (p) is as follows: when the harmonic number is less than 10, the error is large. By performing logarithmic interpolation after the harmonic number is greater than 10, such as under the condition of 1 GHz, the interpolation interval of the harmonic number p = 3, the amount of calculation can be reduced by half, and the efficiency of dielectric scanning logging forward modeling can be improved.

[0084] Figure 6 , Figure 7 This is a specific example diagram of the forward modeling calculation method of the push-type dielectric scanning logging tool according to an embodiment of the present invention.

[0085] In one embodiment of the present invention, equivalent calculation models for measuring a centrally located columnar metal core and an eccentrically located columnar metal core are established and compared using forward modeling. The position of the centrally located columnar metal core is referenced. Figure 6 Part (1), the position reference of the eccentric columnar metal core Figure 6 Part (2), r bh The radius r represents the cylindrical metal core. tool Indicates the distances of the transmitting and receiving coils from the center point; Reference Figure 7 (1) part and Figure 7 Part (2) shows the imaginary part of the magnetic field component (Imag(H)). zz The distance between the receiving coil and the transmitting coil (L) TR The relationship between the magnetic field components and the relative error of each calculation method and the distance between the receiving coil and the transmitting coil is discussed. The Imag function is a function for calculating the imaginary part of the magnetic field components. Comparing the calculation method of this invention based on a centrally located cylindrical metal core, the software simulation calculation based on a centrally located cylindrical metal core, and the software simulation calculation based on an eccentric cylindrical metal core, it can be seen that the forward modeling calculation method of this invention under the measurement equivalent calculation model of a centrally located cylindrical metal core matches the results of commercial software, verifying its accuracy. Compared with the eccentric cylindrical metal core, the difference in response results increases with the source distance (distance between the receiving coil and the transmitting coil), reaching a maximum of about 10%, and the stability is good at different frequencies.

[0086] In one embodiment of the present invention, after obtaining the electric field components and magnetic field components of the signal receiving end when receiving signals transmitted from different signal transmitting ends, the amplitude ratio and phase difference of the electric field components of the signal receiving end when receiving signals transmitted from different signal transmitting ends are determined; the amplitude ratio and phase difference of the magnetic field components of the signal receiving end when receiving signals transmitted from different signal transmitting ends are also determined.

[0087] Figure 8 This is a specific example diagram of the forward modeling calculation method of the push-type dielectric scanning logging tool according to an embodiment of the present invention.

[0088] In practical implementation, taking two transmitting coils and eight receiving coils as an example, refer to... Figure 8 In part (1), TA and TB represent transmitting coils, and RA1, RA2, RA3, RA4, RB1, RB2, RB3, and RB4 all represent receiving coils; the distance between two transmitting coils and the distance between adjacent receiving coils are both 1 inch (in), and the distance between a transmitting coil and an adjacent receiving coil is 1.5 inches. Each set of transmitting coils includes two coils, one horizontal and one vertical. (Refer to...) Figure 8 Part (2), This indicates that the RA receiving coil receives the signal from the TA transmitting coil. This indicates that the RA receiving coil receives the signal from the TA transmitting coil. This indicates that an RB receiving coil receives a signal from a TB transmitting coil. This indicates that an RB receiving coil receives a signal from a TA transmitting coil. The signal, when calculated, includes both electric and magnetic field components. The dielectric logging response is expressed using the amplitude ratio Att and phase difference PS scale. The amplitude ratio of the signal (electric and magnetic field components) is calculated using the following formula:

[0089]

[0090]

[0091]

[0092] Among them, Att A Att represents the amplitude ratio of the signals received by the receiving coil of RA from the transmitting coils TA and TB; B This represents the amplitude ratio of the signals received by the RB receiving coil from the two transmitting coils TA and TB; Att represents Att. A and Att B The mean; This indicates that the RA receiving coil receives the signal from the TA transmitting coil. This indicates that the RA receiving coil receives the signal from the TA transmitting coil. This indicates that an RB receiving coil receives a signal from a TB transmitting coil. This indicates that an RB receiving coil receives the signal from a TA transmitting coil, and abs indicates taking the absolute value;

[0093] The phase difference is calculated using the following formula:

[0094]

[0095]

[0096]

[0097] Where, Δφ A Δφ represents the phase difference between the signals received by the receiving coil of RA and the transmitting coils of TB. B One RB receiving coil receives the phase difference between the signals from the two transmitting coils TA and TB; PS represents Δφ. A and Δφ B The mean of , where arg represents the phase.

[0098] In one embodiment of the present invention, taking two transmitting coils and eight receiving coils as an example, the dielectric scanner requires calculation of the horizontal and vertical components of the dual-transmitter, eight-receiver coil system at four operating frequencies for each logging point, totaling 128 calculations of field components. This results in a large computational load and high time consumption. Therefore, based on the invariance of the equivalent formation vertical properties and the symmetry between the receiving and transmitting coils, the present invention first reduces the number of field component calculations for the dual-transmitter, eight-receiver coil system from 16 to 5. Simultaneously, it utilizes the electric field spectral domain field e... zzv , electric field spectral domain h zzv Regardless of z, the spectral field and integration are calculated separately. For the same measurement point, only one calculation of the spectral field is needed, followed by five integrations. As shown in Table 1, this invention can improve the calculation efficiency of dielectric scanning logging response at a single measurement point by more than 8 times.

[0099] Table 1

[0100]

[0101] Figure 9 This is a specific example diagram of the forward modeling calculation method of the push-type dielectric scanning logging tool according to an embodiment of the present invention.

[0102] In one embodiment of the present invention, the amplitude ratio (unit: dB) and phase difference (unit: degrees) responses under different detection modes are calculated, with reference to... Figure 9 The formation resistivity was set between 1 and 500 Ω·m, and the relative permittivity (ε) between 1 and 100, to establish a response calibration chart. The amplitude ratio and phase difference signals at the four frequencies showed significant differences, as referenced... Figure 9 In part (1), at 20MHz, the dielectric constant increase curve is almost parallel to the vertical axis, meaning the phase difference does not reflect the dielectric constant at all, and it is highly sensitive to resistivity. However, numerically, the phase difference variation range is within 10, and the amplitude ratio variation is within 0.1, corresponding to a resistivity variation of 1 to 500 Ω·m. At this frequency and source distance, the accuracy of reflecting electrical parameters is low, making it difficult to accurately extract the corresponding formation parameters. As the frequency increases, the dielectric constant increase curve approaches the vertical axis at an angle closer to the vertical axis, meaning that at high frequencies, the dielectric constant is reflected by the phase difference, which is more obvious at high resistivity. (Refer to...) Figure 9 In part (4), at 1 GHz, when the relative permittivity is greater than 20, the dielectric increase curve is basically parallel to the vertical axis, which is beneficial for accurately extracting the dielectric constant. For resistivity, low frequencies are generally more advantageous, but due to the small distance between the dielectric transmitting coils, the signal difference is weak at low frequencies, and it is generally difficult to determine the resistivity of high-resistivity strata under 20 MHz conditions. (Refer to...) Figure 9For parts (2) and (3), the resistivity of high-resistivity formations needs to be determined by combining 200MHz and 500MHz frequencies. When the resistivity is greater than 10Ω·m, the phase difference signal can better reflect the magnitude of the resistivity. The calibration accuracy of formation resistivity decreases as resistivity increases and dielectric constant increases, while the dielectric constant extraction accuracy is high for high-resistivity formations.

[0103] It should be noted that although the operation of the method of the present invention has been described in a specific order in the above embodiments and figures, this does not require or imply that the operations must be performed in that specific order, or that all the operations shown must be performed to achieve the desired result. Additionally or alternatively, certain steps may be omitted, multiple steps may be combined into one step, and / or one step may be broken down into multiple steps.

[0104] The implementation of the forward modeling calculation device for the push-type dielectric scanning logging tool can refer to the implementation of the above method, and the repeated parts will not be described again. The term "module" or "unit" used below can be a combination of software and / or hardware that implements a predetermined function. Although the apparatus described in the following embodiments is preferably implemented in software, hardware implementation, or a combination of software and hardware, is also possible and contemplated.

[0105] Based on the same inventive concept, this invention also proposes a forward modeling calculation device for a push-type dielectric scanning logging tool, such as... Figure 10 As shown, the device includes:

[0106] The model acquisition module 1001 is used to acquire the measurement equivalent calculation model of the push-type dielectric scanning logging tool; wherein, the measurement equivalent calculation model is obtained by equating the push arm of the push-type dielectric scanning logging tool to a columnar metal core, the antenna to a magnetic dipole, and the formation to a columnar formation.

[0107] The first determining module 1002 is used to determine the electric field spectral domain field and magnetic field spectral domain field of the signal receiving end when receiving signals transmitted from different signal transmitting ends, based on the position of the signal receiving end in the measurement equivalent calculation model.

[0108] The second determining module 1003 is used to determine the relationship between the electric field harmonic components and the harmonic number, and the relationship between the magnetic field harmonic components and the harmonic number, when the signal receiving end receives signals transmitted by different signal transmitting ends, based on the electric field spectral domain field and the magnetic field spectral domain field.

[0109] The third determining module 1004 is used to determine the electric field components of the signal receiver when receiving signals transmitted from different signal transmitters, based on the relationship between the electric field harmonic components and the harmonic number, the interpolation interval of the preset harmonic number, and the first preset harmonic number cutoff value; and to determine the magnetic field components of the signal receiver when receiving signals transmitted from different signal transmitters, based on the relationship between the magnetic field harmonic components and the harmonic number, the interpolation interval of the preset harmonic number, and the first preset harmonic number cutoff value.

[0110] The fourth determining module 1005 is used to determine the amplitude ratio and phase difference of the electric field components when the signal receiver receives signals transmitted from different signal transmitters; and to determine the amplitude ratio and phase difference of the magnetic field components when the signal receiver receives signals transmitted from different signal transmitters.

[0111] In one embodiment of the present invention, the first determining module 1002 is specifically used for:

[0112] Based on the position of the signal receiver in the equivalent measurement model, the coupling relationship of the field at the radial layer interface of the equivalent measurement model is described by the generalized reflection array and transmission array, and the electric field spectral domain and magnetic field spectral domain are determined when the signal receiver receives signals transmitted from different signal transmitters.

[0113] In one embodiment of the present invention, it further includes:

[0114] The component representation acquisition module is used to acquire the electric field components and magnetic field components of the Fourier series expansion of the signal receiver in the measurement equivalent calculation model before determining the relationship between the electric field harmonic components and the harmonic number when the signal receiver receives signals transmitted from different signal transmitters.

[0115] The second determining module 1003 is specifically used for:

[0116] Based on the electric field spectral domain field, the magnetic field spectral domain field, and the electric and magnetic field components of the Fourier series expansion, the relationship between the electric field harmonic components and the harmonic number, and the relationship between the magnetic field harmonic components and the harmonic number are determined when the signal receiver receives signals transmitted from different signal transmitters.

[0117] In one embodiment of the present invention, it further includes:

[0118] The interpolation interval determination module is used to determine the electric field components of the signal receiver when receiving signals transmitted from different signal transmitters, based on the relationship between the electric field harmonic components and the harmonic number, the interpolation interval of the preset harmonic number, and the first preset harmonic number cutoff value; and to determine the magnetic field components of the signal receiver when receiving signals transmitted from different signal transmitters, based on the relationship between the magnetic field harmonic components and the harmonic number, the interpolation interval of the preset harmonic number, and the first preset harmonic number cutoff value. Before determining the magnetic field components of the signal receiver when receiving signals transmitted from a signal transmitter, the module determines the first magnetic field component of the signal receiver when receiving a signal transmitted from a signal transmitter, based on the relationship between the magnetic field harmonic components and the harmonic number and the first preset harmonic number cutoff value; and to determine the second magnetic field component of the signal receiver when receiving a signal transmitted from the same signal transmitter, based on the relationship between the magnetic field harmonic components and the harmonic number and the second preset harmonic number cutoff value, wherein the second preset harmonic number cutoff value is greater than the first preset harmonic number cutoff value; to determine the relationship between the error between the second and first magnetic field components and the harmonic number; and to determine the interpolation interval of the preset harmonic number based on the relationship between the error between the second and first magnetic field components and the harmonic number.

[0119] In this embodiment of the invention, the interpolation interval determination module is specifically used for:

[0120] Based on the error between the second magnetic field component and the first magnetic field component, a local error function is determined; wherein, the local error function is used to represent the relationship between the error between the second magnetic field component and the first magnetic field component and the harmonic number;

[0121] Based on the local error function, determine the interpolation interval for the preset harmonic number;

[0122] The local error function is as follows:

[0123]

[0124]

[0125] Among them, Loc err (v) represents the local error function; e p (v) represents the error between the second magnetic field component and the first magnetic field component; v represents the harmonic number; N1 represents the first preset harmonic number cutoff value; N2 represents the second preset harmonic number cutoff value; H represents the first magnetic field component; zz The second magnetic field component is represented by τ, which is a preset constant; P represents the second magnetic field component. p (v) represents an nth-order interpolation polynomial; where P p (v)=a0+a1v+a2v 2 +...+a n v n a0, a1, a2, ..., a nAll are constant coefficients of the nth-order interpolation polynomial; n is the order of the interpolation polynomial; p is the interpolation interval of the harmonic number.

[0126] It should be noted that although several modules of the forward modeling calculation device of the push-type dielectric scanning logging tool have been mentioned in the detailed description above, this division is merely exemplary and not mandatory. In fact, according to embodiments of the present invention, the features and functions of two or more modules described above can be embodied in one module. Conversely, the features and functions of one module described above can be further divided and embodied by multiple modules.

[0127] Based on the aforementioned inventive concept, such as Figure 11 As shown, the present invention also proposes a computer device 1100, including a memory 1101, a processor 1102, and a computer program 1103 stored in the memory 1101 and executable on the processor 1102. When the processor 1102 executes the computer program 1103, it implements the aforementioned forward modeling method of the push-type dielectric scanning logging tool.

[0128] Based on the aforementioned inventive concept, the present invention proposes a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the forward modeling method of the aforementioned push-type dielectric scanning logging tool.

[0129] Based on the aforementioned inventive concept, this invention proposes a computer program product, which includes a computer program that, when executed by a processor, implements a forward modeling method for a push-type dielectric scanning logging tool.

[0130] The forward modeling method and apparatus for the push-type dielectric scanning logging tool proposed in this invention can solve the problems of low efficiency and high memory consumption of existing three-dimensional numerical methods. This invention obtains a measurement equivalent calculation model for the push-type dielectric scanning logging tool. This model is obtained by equating the push arm of the push-type dielectric scanning logging tool to a columnar metal core, the antenna to a magnetic dipole, and the formation to a columnar formation. Based on the position of the signal receiving end in the measurement equivalent calculation model, the electric field and magnetic field spectral domains of the signal receiving end when receiving signals transmitted from different signal transmitting ends are determined. Based on the electric field and magnetic field spectral domains, the calculation of the signal receiving end's response to different signal transmissions is determined. The relationship between the electric field harmonic components and the harmonic number when the signal is transmitted is determined; the relationship between the magnetic field harmonic components and the harmonic number is also determined; based on the relationship between the electric field harmonic components and the harmonic number, the interpolation interval of the preset harmonic number, and the first preset harmonic number cutoff value, the electric field components of the signal receiver when receiving signals transmitted from different signal transmitters are determined; based on the relationship between the magnetic field harmonic components and the harmonic number, the interpolation interval of the preset harmonic number, and the first preset harmonic number cutoff value, the magnetic field components of the signal receiver when receiving signals transmitted from different signal transmitters are determined; the amplitude ratio and phase difference of the electric field components of the signal receiver when receiving signals transmitted from different signal transmitters are determined; the amplitude ratio and phase difference of the magnetic field components of the signal receiver when receiving signals transmitted from different signal transmitters are determined. This invention provides a fast forward modeling algorithm for amplitude ratio and phase difference response of dielectric scanning logging. It uses the measurement equivalent calculation model of the push-type dielectric scanning logging tool as the basic model for forward modeling and presets the interpolation interval of harmonic numbers during the forward modeling process, which greatly reduces the amount of calculation and can improve the forward modeling efficiency of the push-type dielectric scanning logging tool and reduce the computer memory occupied.

[0131] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0132] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0133] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0134] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0135] The specific embodiments described above further illustrate the purpose, technical solution, and beneficial effects of the present invention. It should be understood that the above descriptions are merely specific embodiments of the present invention and are not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A forward modeling method for a push-type dielectric scanning logging tool, characterized in that, include: Obtain the measurement equivalent calculation model of the push-type dielectric scanning logging tool; wherein, the measurement equivalent calculation model is obtained by equating the push arm of the push-type dielectric scanning logging tool to a columnar metal core, the antenna to a magnetic dipole, and the formation to a columnar formation; Based on the position of the signal receiver in the equivalent calculation model, determine the electric field spectral domain and magnetic field spectral domain when the signal receiver receives signals transmitted from different signal transmitters. Based on the electric field spectral domain and the magnetic field spectral domain, determine the relationship between the electric field harmonic components and the harmonic number when the signal receiver receives signals transmitted from different signal transmitters, and the relationship between the magnetic field harmonic components and the harmonic number. Based on the relationship between electric field harmonic components and harmonic numbers, the interpolation interval of preset harmonic numbers, and the first preset harmonic number cutoff value, the electric field components of the signal receiver when receiving signals transmitted from different signal transmitters are determined; based on the relationship between magnetic field harmonic components and harmonic numbers, the interpolation interval of preset harmonic numbers, and the first preset harmonic number cutoff value, the magnetic field components of the signal receiver when receiving signals transmitted from different signal transmitters are determined. Determine the amplitude ratio and phase difference of the electric field components when the signal receiver receives signals transmitted from different signal transmitters; determine the amplitude ratio and phase difference of the magnetic field components when the signal receiver receives signals transmitted from different signal transmitters.

2. The method according to claim 1, characterized in that, Based on the position of the signal receiver in the equivalent measurement model, determine the electric and magnetic field spectral domains of the signal receiver when receiving signals transmitted from different signal transmitters, including: Based on the position of the signal receiver in the equivalent measurement model, the coupling relationship of the field at the radial layer interface of the equivalent measurement model is described by the generalized reflection array and transmission array, and the electric field spectral domain and magnetic field spectral domain are determined when the signal receiver receives signals transmitted from different signal transmitters.

3. The method according to claim 1, characterized in that, Before determining the relationship between the electric field harmonic components and the harmonic number when the signal receiver receives signals transmitted from different signal transmitters, and the relationship between the magnetic field harmonic components and the harmonic number, the following is included: Obtain the electric and magnetic field components of the Fourier series expansion of the signal receiver in the equivalent calculation model of the measurement. Based on the electric and magnetic field spectral domains, the relationships between the electric field harmonic components and harmonic numbers, and between the magnetic field harmonic components and harmonic numbers, are determined when the signal receiver receives signals transmitted from different signal transmitters. These relationships include: Based on the electric field spectral domain field, the magnetic field spectral domain field, and the electric and magnetic field components of the Fourier series expansion, the relationship between the electric field harmonic components and the harmonic number, and the relationship between the magnetic field harmonic components and the harmonic number are determined when the signal receiver receives signals transmitted from different signal transmitters.

4. The method according to claim 1, characterized in that, Before determining the electric field components of the signal receiver when receiving signals transmitted from different signal transmitters, based on the relationship between electric field harmonic components and harmonic numbers, the interpolation interval of preset harmonic numbers, and the first preset harmonic number cutoff value, the method further includes: Based on the relationship between the magnetic field harmonic components and the harmonic number, and the first preset harmonic number cutoff value, the first magnetic field component is determined when the signal receiver receives a signal transmitted by a signal transmitter. Based on the relationship between the magnetic field harmonic components and the harmonic number, and the second preset harmonic number cutoff value, the second magnetic field component is determined when the signal receiver receives the signal transmitted by the signal transmitter; wherein, the second preset harmonic number cutoff value is greater than the first preset harmonic number cutoff value. Determine the error between the second magnetic field component and the first magnetic field component; The interpolation interval for the preset harmonic number is determined based on the error between the second magnetic field component and the first magnetic field component.

5. The method according to claim 4, characterized in that, Based on the error between the second magnetic field component and the first magnetic field component, the interpolation interval for the preset harmonic number is determined, including: Based on the error between the second magnetic field component and the first magnetic field component, a local error function is determined; wherein, the local error function is used to represent the relationship between the error between the second magnetic field component and the first magnetic field component and the harmonic number; Based on the local error function, determine the interpolation interval for the preset harmonic number; The local error function is as follows: Among them, Loc err (v) represents the local error function; e p (v) represents the error between the second magnetic field component and the first magnetic field component; v represents the harmonic number; N1 represents the first preset harmonic number cutoff value; N2 represents the second preset harmonic number cutoff value; H represents the first magnetic field component; zz The second magnetic field component is represented by τ, which is a preset constant; P represents the second magnetic field component. p (v) represents an nth-order interpolation polynomial; where P p (v)=a0+a1v+a2v 2 +...+a n v n a0, a1, a2, ..., a n All are constant coefficients of the nth-order interpolation polynomial; n is the order of the interpolation polynomial; p is the interpolation interval of the harmonic number.

6. A forward modeling calculation device for a push-type dielectric scanning logging tool, characterized in that, include: The model acquisition module is used to acquire the measurement equivalent calculation model of the push-type dielectric scanning logging tool; wherein, the measurement equivalent calculation model is obtained by equating the push arm of the push-type dielectric scanning logging tool to a columnar metal core, the antenna to a magnetic dipole, and the formation to a columnar formation. The first determining module is used to determine the electric field spectral domain field and magnetic field spectral domain field of the signal receiving end when receiving signals transmitted from different signal transmitting ends, based on the position of the signal receiving end in the measurement equivalent calculation model. The second determining module is used to determine the relationship between the electric field harmonic components and the harmonic number, and the relationship between the magnetic field harmonic components and the harmonic number, when the signal receiver receives signals transmitted by different signal transmitters, based on the electric field spectral domain field and the magnetic field spectral domain field. The third determining module is used to determine the electric field components of the signal receiver when receiving signals transmitted from different signal transmitters, based on the relationship between the electric field harmonic components and the harmonic number, the interpolation interval of the preset harmonic number, and the first preset harmonic number cutoff value; and to determine the magnetic field components of the signal receiver when receiving signals transmitted from different signal transmitters, based on the relationship between the magnetic field harmonic components and the harmonic number, the interpolation interval of the preset harmonic number, and the first preset harmonic number cutoff value. The fourth determining module is used to determine the amplitude ratio and phase difference of the electric field components when the signal receiver receives signals transmitted from different signal transmitters; and to determine the amplitude ratio and phase difference of the magnetic field components when the signal receiver receives signals transmitted from different signal transmitters.

7. The apparatus according to claim 6, characterized in that, The first determining module is specifically used for: Based on the position of the signal receiver in the equivalent measurement model, the coupling relationship of the field at the radial layer interface of the equivalent measurement model is described by the generalized reflection array and transmission array, and the electric field spectral domain and magnetic field spectral domain are determined when the signal receiver receives signals transmitted from different signal transmitters.

8. The apparatus according to claim 6, characterized in that, Also includes: The component representation acquisition module is used to acquire the electric field components and magnetic field components of the Fourier series expansion of the signal receiver in the measurement equivalent calculation model before determining the relationship between the electric field harmonic components and the harmonic number when the signal receiver receives signals transmitted from different signal transmitters. The second determining module is specifically used for: Based on the electric field spectral domain field, the magnetic field spectral domain field, and the electric and magnetic field components of the Fourier series expansion, the relationship between the electric field harmonic components and the harmonic number, and the relationship between the magnetic field harmonic components and the harmonic number are determined when the signal receiver receives signals transmitted from different signal transmitters.

9. The apparatus according to claim 6, characterized in that, Also includes: The interpolation interval determination module is used to determine the electric field components of the signal receiver when receiving signals transmitted from different signal transmitters, based on the relationship between the electric field harmonic components and the harmonic number, the interpolation interval of the preset harmonic number, and the first preset harmonic number cutoff value. Before determining the magnetic field components of the signal receiver when receiving signals transmitted from different signal transmitters, based on the relationship between the magnetic field harmonic components and the harmonic number, the interpolation interval of the preset harmonic number, and the first preset harmonic number cutoff value, the module determines the first magnetic field component of the signal receiver when receiving a signal transmitted from a signal transmitter, based on the relationship between the magnetic field harmonic components and the harmonic number and the first preset harmonic number cutoff value. Furthermore, based on the relationship between the magnetic field harmonic components and the harmonic number and the second preset harmonic number cutoff value, the module determines the second magnetic field component of the signal receiver when receiving a signal transmitted from the same signal transmitter. The second preset harmonic number cutoff value is greater than the first preset harmonic number cutoff value. Determine the relationship between the error between the second magnetic field component and the first magnetic field component and the harmonic number; based on the relationship between the error between the second magnetic field component and the first magnetic field component and the harmonic number, determine the interpolation interval of the preset harmonic number.

10. The apparatus according to claim 9, characterized in that, The interpolation interval determination module is specifically used for: Based on the error between the second magnetic field component and the first magnetic field component, a local error function is determined; wherein, the local error function is used to represent the relationship between the error between the second magnetic field component and the first magnetic field component and the harmonic number; Based on the local error function, determine the interpolation interval for the preset harmonic number; The local error function is as follows: Among them, Loc err (v) represents the local error function; e p (v) represents the error between the second magnetic field component and the first magnetic field component; v represents the harmonic number; N1 represents the first preset harmonic number cutoff value; N2 represents the second preset harmonic number cutoff value; H represents the first magnetic field component; zz The second magnetic field component is represented by τ, which is a preset constant; P represents the second magnetic field component. p (v) represents an nth-order interpolation polynomial; where P p (v)=a0+a1v+a2v 2 +...+a n v n a0, a1, a2, ..., a n All are constant coefficients of the nth-order interpolation polynomial; n is the order of the interpolation polynomial; p is the interpolation interval of the harmonic number.

11. A computer 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 computer program, it implements the method of any one of claims 1 to 5.

12. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, implements the method of any one of claims 1 to 5.

13. A computer program product, characterized in that, The computer program product includes a computer program that, when executed by a processor, implements the method of any one of claims 1 to 5.

Images