Forward modeling calculation method and device of pushing type dielectric scanning logging instrument

Abstract

The invention provides a forward modeling calculation method and device for a sidewall contact type dielectric scanning logging instrument. The method comprises the steps that a measurement equivalent calculation model of the sidewall contact type dielectric scanning logging instrument is obtained; determining an electric field spectral domain field and a magnetic field spectral domain field when the signal receiving end receives signals transmitted by different signal transmitting ends; determining a relationship between an electric field harmonic component and a harmonic number and a relationship between a magnetic field harmonic component and a harmonic number according to the electric field spectral domain field and the magnetic field spectral domain field; according to the relationship between the electric field harmonic component and the harmonic number, the interpolation interval of the preset harmonic number and the first preset harmonic number cut-off value, determining the electric field component of the signal receiving end when the signal receiving end receives the transmitted signals of different signal transmitting ends; determining the magnetic field components of the signal receiving end when receiving signals transmitted by different signal transmitting ends according to the relationship between the magnetic field harmonic component and the harmonic number, the interpolation interval of the preset harmonic number and the first preset harmonic number cutoff value; and determining an amplitude ratio and a phase difference of the electric field component and the magnetic field component.

Description

Technical Field

[0001] The invention relates to the technical field of dielectric well logging, in particular to a forward calculation method and device for a push-type dielectric scanning well logging instrument. Background Art

[0003] The dielectric constant is a key parameter for distinguishing oil from water. At the same time, the dielectric constant of natural gas hydrate-containing sediments is an important parameter in the exploration and development of hydrate resources. Dielectric scanning logging is an important logging method for the exploration and development of unconventional oil and gas energy such as natural gas hydrates, tight oil and gas, heavy oil, and carbonate rocks. The push-type dielectric scanning logging instrument uses a push-type plate and uses multi-frequency and multi-source to measure formation information. It is the current mainstream dielectric logging instrument. In order to obtain the dielectric constant, it is necessary to first obtain the amplitude ratio and phase difference of the dielectric scanning logging through forward calculation. For the push-type dielectric scanning logging instrument, the push-type arm, antenna, wellbore, and formation constitute a complex three-dimensional geological model. Three-dimensional numerical methods such as finite element and finite difference are used for forward calculation. The forward calculation efficiency of the existing three-dimensional numerical method is very low and the computer memory occupied is high. Summary of the invention

[0004] In an embodiment of the present invention, a forward modeling method of a push-type dielectric scanning well logging tool is proposed, which can improve the forward modeling efficiency of the push-type dielectric scanning well logging tool and reduce the occupied computer memory, including:

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

[0006] According to the position of the signal receiving end in the measurement equivalent calculation model, the electric field spectrum domain and the magnetic field spectrum domain of the signal receiving end when receiving signals transmitted by different signal transmitting ends are determined;

[0007] Determine the relationship between the electric field harmonic component and the harmonic number, and the relationship between the magnetic field harmonic component and the harmonic number when the signal receiving end receives signals transmitted by different signal transmitting ends according to the electric field spectrum domain field and the magnetic field spectrum domain field;

[0008] Determine the electric field component of the signal receiving end when receiving signals transmitted by different signal transmitting ends according to the relationship between the electric field harmonic component and the harmonic number, the interpolation interval of the preset harmonic number and the first preset harmonic number cutoff value; determine the magnetic field component of the signal receiving end when receiving signals transmitted by different signal transmitting ends according to the relationship between the magnetic field harmonic component and the harmonic number, the interpolation interval of the preset harmonic number and the first preset harmonic number cutoff value;

[0009] Determine the amplitude ratio and phase difference of the electric field components when the signal receiving end receives signals transmitted by different signal transmitting ends; determine the amplitude ratio and phase difference of the magnetic field components when the signal receiving end receives signals transmitted by different signal transmitting ends.

[0010] In an embodiment of the present invention, a forward modeling computing device for a push-type dielectric scanning well logging tool is proposed, which can improve the forward modeling computing efficiency of the push-type dielectric scanning well logging tool and reduce the occupied computer memory, 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-type dielectric scanning logging tool's push arm to a cylindrical metal core, the antenna to a magnetic dipole, and the formation to a cylindrical formation;

[0012] A first determination module is used to determine the electric field spectrum domain and the magnetic field spectrum domain of the signal receiving end when receiving signals transmitted by different signal transmitting ends according to the position of the signal receiving end in the measurement equivalent calculation model;

[0013] The second determination module is used to determine the relationship between the electric field harmonic component and the harmonic number, and the relationship between the magnetic field harmonic component and the harmonic number when the signal receiving end receives signals transmitted by different signal transmitting ends according to the electric field spectrum domain field and the magnetic field spectrum domain field;

[0014] A third determination module is used to determine the electric field component of the signal receiving end when receiving signals transmitted by different signal transmitting ends according to the relationship between the electric field harmonic component 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 component of the signal receiving end when receiving signals transmitted by different signal transmitting ends according to the relationship between the magnetic field harmonic component and the harmonic number, the interpolation interval of the preset harmonic number and the first preset harmonic number cutoff value;

[0015] The fourth determination module is used to determine the amplitude ratio and phase difference of the electric field components when the signal receiving end receives signals transmitted by different signal transmitting ends; determine the amplitude ratio and phase difference of the magnetic field components when the signal receiving end receives signals transmitted by different signal transmitting ends.

[0016] In an embodiment of the present invention, a computer device is proposed, 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, a forward modeling method of a push-type dielectric scanning logging tool is implemented.

[0017] In an embodiment of the present invention, a computer-readable storage medium is provided. The computer-readable storage medium stores a computer program. When the computer program is executed by a processor, a forward modeling method of a push-type dielectric scanning logging tool is implemented.

[0018] A computer program product is provided in an embodiment of the present invention. The computer program product includes a computer program. When the computer program is executed by a processor, a forward modeling method of a push-type dielectric scanning logging tool is implemented.

[0019] The forward modeling method and device of the push-type dielectric scanning logging instrument proposed in the embodiment of the present invention can solve the problems that the forward modeling efficiency of the existing three-dimensional numerical method is very low and the computer memory occupied is high; the embodiment of the present invention obtains the measurement equivalent calculation model of the push-type dielectric scanning logging instrument; wherein the measurement equivalent calculation model is obtained by equivalently treating the push-type dielectric scanning logging instrument's push arm as a cylindrical metal core, the antenna as a magnetic dipole, and the formation as a columnar formation; according to the position of the signal receiving end in the measurement equivalent calculation model, the electric field spectrum domain field and the magnetic field spectrum domain field of the signal receiving end when receiving signals transmitted by different signal transmitting ends are determined; according to the electric field spectrum domain field and the magnetic field spectrum domain field, the signal receiving end is determined when receiving signals transmitted by different signal transmitting ends. The invention relates to a method for determining a relationship between an electric field harmonic component and a harmonic number, and a relationship between a magnetic field harmonic component and a harmonic number when a signal receiving end transmits a signal; determining an electric field component when a signal receiving end receives signals transmitted by different signal transmitting ends according to the relationship between the electric field harmonic component and the harmonic number, an interpolation interval of a preset harmonic number, and a first preset harmonic number cutoff value; determining a magnetic field component when a signal receiving end receives signals transmitted by different signal transmitting ends according to the relationship between the magnetic field harmonic component and the harmonic number, an interpolation interval of a preset harmonic number, and a first preset harmonic number cutoff value; determining an amplitude ratio and a phase difference of an electric field component when the signal receiving end receives signals transmitted by different signal transmitting ends; and determining an amplitude ratio and a phase difference of a magnetic field component when the signal receiving end receives signals transmitted by different signal transmitting ends. The embodiment of the present invention provides a fast forward algorithm for the amplitude ratio and phase difference response of dielectric scanning logging, uses the measurement equivalent calculation model of the push-type dielectric scanning logging instrument as the basic model for forward calculation, and presets the interpolation interval of the harmonic number during the forward calculation process, which greatly reduces the amount of calculation, improves the forward calculation efficiency of the push-type dielectric scanning logging instrument, and reduces the occupied computer memory. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0021] Figure 1 1 is a flow chart of a forward calculation method of a push-type dielectric scanning logging tool according to an embodiment of the present invention;

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

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

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

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

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

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

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

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

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

[0031] Fig.11 Schematic diagram of a computer device in an embodiment of the present invention. DETAILED DESCRIPTION

[0032] To make the purpose, technical solution and advantages of the embodiments of the present invention more clear, the embodiments of the present invention are further described in detail below in conjunction with the accompanying drawings. Here, the exemplary 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] The term "and / or" herein only describes an association relationship, indicating that three relationships may exist. For example, A and / or B may represent the following three situations: A exists alone, A and B exist at the same time, and B exists alone. In addition, the term "at least one" herein represents any combination of at least two of any one or more of a plurality of. For example, including at least one of A, B, and C may represent 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 "include", "including", "have", "contain", etc. are all open terms, which mean including but not limited to. The descriptions with reference to the terms "one embodiment", "a specific embodiment", "some embodiments", "for example", etc. mean that the specific features, structures or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures or characteristics described may be combined in a suitable manner in any one or more embodiments or examples. The order of steps involved in each embodiment is used to schematically illustrate the implementation of the present application, and the order of steps is not limited and can be appropriately adjusted as needed.

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

[0036] Figure 1 FIG. 1 is a flow chart of a forward calculation method of a push-type dielectric scanning logging tool according to an embodiment of the present invention. Figure 1 As shown, the method includes:

[0037] Step 101, obtaining a measurement equivalent calculation model of a push-type dielectric scanning logging tool; wherein the measurement equivalent calculation model is obtained by equivalently treating the push-type dielectric scanning logging tool's push arm as a cylindrical metal core, the antenna as a magnetic dipole, and the formation as a columnar formation;

[0038] Step 102, determining the electric field spectrum domain and the magnetic field spectrum domain of the signal receiving end when the signal receiving end receives signals transmitted by different signal transmitting ends according to the position of the signal receiving end in the measurement equivalent calculation model;

[0039] Step 103, determining the relationship between the electric field harmonic component and the harmonic number, and the relationship between the magnetic field harmonic component and the harmonic number when the signal receiving end receives signals transmitted by different signal transmitting ends according to the electric field spectrum domain field and the magnetic field spectrum domain field;

[0040] Step 104, determining the electric field component of the signal receiving end when receiving signals transmitted by different signal transmitting ends according to the relationship between the electric field harmonic component and the harmonic number, the interpolation interval of the preset harmonic number, and the first preset harmonic number cutoff value; determining the magnetic field component of the signal receiving end when receiving signals transmitted by different signal transmitting ends according to the relationship between the magnetic field harmonic component and the harmonic number, the interpolation interval of the preset harmonic number, and the first preset harmonic number cutoff value;

[0041] Step 105, determining the amplitude ratio and phase difference of the electric field components when the signal receiving end receives signals transmitted by different signal transmitting ends; determining the amplitude ratio and phase difference of the magnetic field components when the signal receiving end receives signals transmitted by different signal transmitting ends.

[0042] Depend on Figure 1 It can be seen from the process shown that the embodiment of the present invention obtains the measurement equivalent calculation model of the push-type dielectric scanning logging instrument; wherein the measurement equivalent calculation model is obtained by making the push-type dielectric scanning logging instrument's push arm equivalent to a cylindrical metal core, the antenna equivalent to a magnetic dipole, and the formation equivalent to a columnar formation; according to the position of the signal receiving end in the measurement equivalent calculation model, the electric field spectrum domain and the magnetic field spectrum domain of the signal receiving end when receiving signals transmitted by different signal transmitting ends are determined; according to the electric field spectrum domain and the magnetic field spectrum domain, the relationship between the electric field harmonic component and the harmonic number, and the magnetic field harmonic component of the signal receiving end when receiving signals transmitted by different signal transmitting ends are determined. and the harmonic number; according to the relationship between the electric field harmonic component and the harmonic number, the interpolation interval of the preset harmonic number and the first preset harmonic number cutoff value, determine the electric field component of the signal receiving end when receiving signals transmitted by different signal transmitting ends; according to the relationship between the magnetic field harmonic component and the harmonic number, the interpolation interval of the preset harmonic number and the first preset harmonic number cutoff value, determine the magnetic field component of the signal receiving end when receiving signals transmitted by different signal transmitting ends; determine the amplitude ratio and phase difference of the electric field component of the signal receiving end when receiving signals transmitted by different signal transmitting ends; determine the amplitude ratio and phase difference of the magnetic field component of the signal receiving end when receiving signals transmitted by different signal transmitting ends. The embodiment of the present invention provides a fast forward algorithm for the amplitude ratio and phase difference response of dielectric scanning logging, uses the measurement equivalent calculation model of the push-to-pull dielectric scanning logging instrument as the basic model for forward calculation, and presets the interpolation interval of the harmonic number in the forward calculation process, which greatly reduces the amount of calculation, can improve the forward calculation efficiency of the push-to-pull dielectric scanning logging instrument, and reduce the computer memory occupied.

[0043] In order to explain the forward calculation method of the push-type dielectric scanning logging tool more clearly, each step is described in detail below.

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

[0045] Position reference of dielectric scanner during dielectric logging Figure 2 The dielectric scanner is set in the formation wellbore, and the push-to-leave arm (i.e., the cylindrical 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-to-leave dielectric scanning logging instrument is established. The measurement equivalent calculation model of the push-to-leave dielectric scanning logging instrument is obtained by equating the push-to-leave arm of the push-to-leave dielectric scanning logging instrument to the cylindrical metal core, the antenna to the magnetic dipole, and the formation to the columnar formation.

[0046] In one embodiment of the present invention, according to the position of the signal receiving end in the measurement equivalent calculation model, determining the electric field spectrum domain field and the magnetic field spectrum domain field of the signal receiving end when receiving signals transmitted by different signal transmitting ends includes:

[0047] According to the position of the signal receiving end in the measurement equivalent calculation model, the coupling relationship of the field at the radial layer interface of the measurement equivalent calculation model is described by using the generalized reflection array and transmission array to determine the electric field spectrum domain and magnetic field spectrum domain when the signal receiving end receives signals transmitted by different signal transmitting ends.

[0048] In one embodiment of the present invention, before determining the relationship between the electric field harmonic component and the harmonic number, and the relationship between the magnetic field harmonic component and the harmonic number when the signal receiving end receives signals transmitted by different signal transmitting ends, the method includes:

[0049] Obtaining the electric field component and the magnetic field component of the Fourier series expansion of the signal receiving end in the measurement equivalent calculation model;

[0050] According to the electric field spectrum domain and the magnetic field spectrum domain, the relationship between the electric field harmonic component and the harmonic number, and the relationship between the magnetic field harmonic component and the harmonic number when the signal receiving end receives signals transmitted by different signal transmitting ends are determined, including:

[0051] According to the electric field spectral domain field, the magnetic field spectral domain field, and the electric field components and the magnetic field components of the Fourier series expansion, the relationship between the electric field harmonic component and the harmonic number, and the relationship between the magnetic field harmonic component and the harmonic number when the signal receiving end receives signals transmitted by different signal transmitting ends are determined.

[0052] In the specific implementation, the pseudo-analytical solution of the measurement equivalent calculation model is derived, and a fast calculation method for dielectric scanning logging response of a 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 field component and the magnetic field component are expanded in the form of Fourier series. Taking the signal transmitting end (i.e., the transmitting coil) transmitting the signal in the z direction and the signal receiving end (i.e., the receiving coil) receiving the signal in the z direction as an example, the electric field component and the magnetic field component of the Fourier series expansion are expressed as follows:

[0053]

[0054] Where i represents the imaginary unit; E zz represents the electric field component, H zz represents the magnetic field component. 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. By default, the transmitting coil is located at 0. k z represents the wave number in the z direction; e zzν represents the electric field in the spectral domain; h zzν represents the magnetic field in the spectral domain; represents the angle between the tool and the borehole axis, which can be 0 in the present invention; E xx The electric field components and H xx The magnetic field component can be expressed by E zz Indicates the component and H zz The magnetic field components are found jointly.

[0055] In specific implementation, there are analytical solutions for the electric field spectral domain e and the magnetic field spectral domain h in the cylindrical coordinate system.

[0056] zzνzzν

[0057] The cylindrical multi-layer model uses generalized reflection array and transmission array to describe the coupling relationship of the field at the radial layer interface of the measurement equivalent calculation model, and determines the electric field spectrum domain and magnetic field spectrum domain when the signal receiving end receives signals transmitted by different signal transmitting ends.

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

[0059] In one embodiment of the present invention, according to the relationship between the electric field harmonic component and the harmonic number, the interpolation interval of the preset harmonic number and the first preset harmonic number cutoff value, the electric field component of the signal receiving end when receiving the signal transmitted by different signal transmitting ends is determined, and according to the relationship between the magnetic field harmonic component 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 component of the signal receiving end when receiving the signal transmitted by different signal transmitting ends, the method further includes:

[0060] Step 301, determining a first magnetic field component when a signal receiving end receives a signal transmitted by a signal transmitting end according to a relationship between a magnetic field harmonic component and a harmonic number and a first preset harmonic number cutoff value;

[0061] Step 302, determining a second magnetic field component when the signal receiving end receives the signal transmitted by the signal transmitting end according to the relationship between the magnetic field harmonic component 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;

[0062] Step 303, determining an error between the second magnetic field component and the first magnetic field component;

[0063] Step 304: determining an interpolation interval of a preset harmonic number according to an error between the second magnetic field component and the first magnetic field component.

[0064] In one embodiment of the present invention, determining the interpolation interval of the preset harmonic number according to the error between the second magnetic field component and the first magnetic field component includes:

[0065] Determine a local error function according to the error between the second magnetic field component and the first magnetic field component; 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;

[0066] Determining an interpolation interval of a preset harmonic number according to a local error function;

[0067] Wherein, the local error function is as follows:

[0068]

[0069]

[0070] 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; represents the first magnetic field component; H zz represents the second magnetic field component; τ is a preset constant; P p (v) represents the n-th order interpolation polynomial; where P p (v) = a0 + a1v + a2v 2 +...+a n v n ; a0, a1, a2, ..., a n are the constant coefficients of the n-th order interpolation polynomial; n is the order of the interpolation polynomial; and p is the interpolation interval of the harmonic number.

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

[0072] In the specific implementation, for the problem of summing the electric field component and the magnetic field component of the Fourier series expansion, since the tool and the borehole axis are coplanar, let The electric field component and magnetic field component 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:

[0073]

[0074] Among them, H zz Represents the magnetic field component, the subscript zz indicates that the transmitting coil transmits in the z direction and the receiving coil receives in the z direction; H zzvrepresents the harmonic component of the magnetic field at different harmonic numbers, τ is a preset constant; v represents the harmonic number, when v = 0, τ = 1, otherwise τ = 2; in the above formula, the harmonic number cutoff value N can be preset, and the magnetic field component is expressed as:

[0075]

[0076] When simulating different frequencies and different positions of the transmitting coil and the receiving coil, the harmonic components of the magnetic field |Imag(H zzv )|In the exponential decay law, the Imag function is a function that calculates the imaginary part of a complex number, see Figure 4 Part (1), r bh represents the radius of the cylindrical metal core, r1 and r2 represent the distances between the transmitting coil and the receiving coil and the center point, refer to Figure 4 In part (2), the closer the transmitting coil and the receiving coil are to the cylindrical metal core, the greater the N value is, and the greater the N value is when the frequency is high;

[0077] make:

[0078]

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

[0080] Convert the harmonic components of the magnetic field into logarithmic representation, f v It is the intermediate variable of the formula, representing the form of the harmonic component of the magnetic field after taking the age; P p (v) represents the n-th order interpolation polynomial; where P p (v) = a0 + a1v + a2v 2 +...+a n v n ; a0, a1, a2, ..., a n are the constant coefficients of the n-order interpolation polynomial; n is the order of the interpolation polynomial; p is the interpolation interval of the harmonic number; e p (v) represents the error between the second magnetic field component and the first magnetic field component, wherein the second magnetic field component is calculated under the condition of a second preset harmonic number cutoff value. The larger the second preset harmonic number 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 approximate value, that is, when the error is ignored, the magnetic field component can be approximately represented by the first magnetic field component:

[0081]

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

[0083]

[0084]

[0085] 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; represents the first magnetic field component; H zz represents the second magnetic field component; τ is a preset constant; P p (v) represents the n-th order interpolation polynomial; where P p (v) = a0 + a1v + a2v 2 +...+a n v n ; a0, a1, a2, ..., a n are the constant coefficients of the n-order interpolation polynomial; n is the order of the interpolation polynomial; p is the interpolation interval of the harmonic number; reference Figure 5 , indicating the local error (Loc err ) and the interpolation interval (p) of the harmonic number. When the interpolation method is applied, the error is large when the harmonic number is less than 10. By performing logarithmic interpolation after the harmonic number is greater than 10, such as when the interpolation interval p of the harmonic number is 3 under 1 GHz conditions, the amount of calculation can be reduced by half, thus improving the efficiency of forward modeling of dielectric scanning logging.

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

[0087] In one embodiment of the present invention, measurement equivalent calculation models of a centered cylindrical metal core and an eccentric cylindrical metal core are respectively established for forward calculation comparison. The position of the centered cylindrical metal core is referenced Figure 6 (1) of the part, the position of the eccentric cylindrical metal core is referenced Figure 6 Part (2), r bh represents the radius of the cylindrical metal core, r tool Indicates the distance between the transmitting coil and the receiving coil and 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 )) Distance between the receiving coil and the transmitting coil (L TR), and the relationship between the relative error of each calculation method and the distance between the receiving coil and the transmitting coil. The Imag function is a function of calculating the imaginary part of the magnetic field component; by comparing the calculation method of the present invention based on the centered cylindrical metal core, the software simulation calculation based on the centered cylindrical metal core, and the software simulation calculation based on the eccentric cylindrical metal core, it can be seen that the forward calculation method of the present invention under the measurement equivalent calculation model of the centered cylindrical metal core is consistent with the comparison result of the commercial software, verifying its accuracy; and compared with the eccentric cylindrical metal core, the difference in the response results increases with the increase of the source distance (the distance between the receiving coil and the transmitting coil), and the maximum is about 10%, and the stability is good at different frequencies.

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

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

[0090] In the specific implementation, take two transmitting coils and eight receiving coils as an example, refer to Figure 8 In (1), TA and TB represent transmitting coils, and RA1, RA2, RA3, RA4, RB1, RB2, RB3, and RB4 represent receiving coils. The distance between two transmitting coils and the distance between adjacent receiving coils are both 1 inch (in). The distance between a transmitting coil and an adjacent receiving coil is 1.5 inches. Each group of transmitting coils includes two coils, horizontal and vertical. Figure 8 Part (2) of Indicates that a RA receiving coil receives the signal of a TA transmitting coil, Indicates that a RA receiving coil receives the signal of a TA transmitting coil, Indicates that a RB receiving coil receives the signal of a TB transmitting coil, It indicates that a RB receiving coil receives the signal of a TA transmitting coil. The signal includes electric field component and magnetic field component during calculation. The dielectric logging response is expressed by amplitude ratio Att and phase difference PS scale. The amplitude ratio of the signal (electric field component, magnetic field component) is calculated according to the following formula:

[0091]

[0092]

[0093]

[0094] Among them, Att A Att represents the amplitude ratio of the signals received by the RA receiving coil from the TA and TB transmitting coils; B It represents the amplitude ratio of the signal received by the RB receiving coil from the TA and TB transmitting coils; Att represents Att A and Att B The mean of Indicates that a RA receiving coil receives the signal of a TA transmitting coil, Indicates that a RA receiving coil receives the signal of a TA transmitting coil, Indicates that a RB receiving coil receives the signal of a TB transmitting coil, It means that a RB receiving coil receives the signal of a TA transmitting coil, and abs means taking the absolute value;

[0095] The phase difference is calculated as follows:

[0096]

[0097]

[0098]

[0099] Among them, Δφ A Indicates the phase difference between the signals received by the RA receiving coil and the TA and TB transmitting coils; Δφ B The phase difference between the signals received by the RB receiving coil and the signals received by the TA and TB transmitting coils; PS represents Δφ A and Δφ B The mean of , arg represents the phase.

[0100] In one embodiment of the present invention, two transmitting coils and eight receiving coils are used as an example. Since the dielectric scanner needs to calculate the horizontal and vertical components of the dual transmitting and eight receiving coils system at four working frequencies for each logging point, a total of 128 field value component calculations are required, which is computationally intensive and time-consuming. Therefore, the present invention first reduces the number of field component calculations of the dual transmitting and eight receiving coils system from 16 to 5 based on the invariance of the longitudinal properties of the equivalent formation and the symmetry of the receiving coil and the transmitting coil. At the same time, the electric field spectrum domain e zzv 、Electric field spectrum domain h zzv Independent of z, the spectrum domain field and the integral are calculated separately, and the spectrum domain field only needs to be calculated once for the same measurement point, and then five integrals are performed. As can be seen from Table 1, the present invention can improve the calculation efficiency of the dielectric scanning logging response of a single measurement point by more than 8 times.

[0101] Table 1

[0102]

[0103] Fig. 9 It is a specific example diagram of the forward calculation method of the push-type dielectric scanning logging tool according to the embodiment of the present invention.

[0104] In one embodiment of the present invention, the amplitude ratio (unit: dB) and phase difference (unit: degree) responses under different detection modes are calculated, referring to Fig. 9 , the formation resistivity range is set between 1 and 500Ω·m, the relative dielectric constant (ε) range is between 1 and 100, and a response scale chart is established. The amplitude ratio and phase difference signals corresponding to the four frequencies are obviously different. Fig. 9 In part (1), at 20MHz, the dielectric constant increase curve is almost parallel to the vertical axis, that is, the phase difference basically does not reflect the dielectric constant, and the sensitivity to resistivity is high. However, from the numerical point of view, 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 reflection accuracy of electrical parameters is low, and it is difficult to accurately extract the corresponding formation parameters. As the frequency increases, the dielectric constant increase curve approaches the vertical angle with the vertical axis, that is, under high frequency conditions, the dielectric constant is reflected by the phase difference, which is more obvious at high resistance. Fig. 9 In part (4), at 1 GHz, when the relative dielectric constant 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 frequency is generally more favorable, but due to the small distance between the dielectric transmitting coils and the weak signal difference at low frequency, it is generally difficult to determine the resistivity of high-resistance formations under 20 MHz conditions. Fig. 9 Part (2) and part (3) of the above method need to combine 200MHz and 500MHz frequencies to jointly determine the resistivity of high-resistance formations. 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 with the increase of resistivity and dielectric constant, while the accuracy of dielectric constant extraction in high-resistance formations is high.

[0105] It should be noted that, although the operations of the method of the present invention are described in a specific order in the above embodiments and the accompanying drawings, this does not require or imply that the operations must be performed in the specific order, or that all the operations shown must be performed to achieve the desired results. Additionally or alternatively, some steps may be omitted, multiple steps may be combined into one step, and / or one step may be decomposed into multiple steps.

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

[0107] Based on the same inventive concept, the present invention also proposes a forward calculation device for a push-type dielectric scanning logging tool, such as Fig.10 As shown, the device comprises:

[0108] The model acquisition module 1001 is used to acquire a measurement equivalent calculation model of a push-type dielectric scanning logging tool; wherein the measurement equivalent calculation model is obtained by equivalently treating the push-type dielectric scanning logging tool's push arm as a cylindrical metal core, the antenna as a magnetic dipole, and the formation as a columnar formation;

[0109] The first determination module 1002 is used to determine the electric field spectrum domain and the magnetic field spectrum domain when the signal receiving end receives signals transmitted by different signal transmitting ends according to the position of the signal receiving end in the measurement equivalent calculation model;

[0110] The second determination module 1003 is used to determine the relationship between the electric field harmonic component and the harmonic number, and the relationship between the magnetic field harmonic component and the harmonic number when the signal receiving end receives signals transmitted by different signal transmitting ends according to the electric field spectrum domain field and the magnetic field spectrum domain field;

[0111] The third determination module 1004 is used to determine the electric field component of the signal receiving end when receiving signals transmitted by different signal transmitting ends according to the relationship between the electric field harmonic component and the harmonic number, the interpolation interval of the preset harmonic number and the first preset harmonic number cutoff value; determine the magnetic field component of the signal receiving end when receiving signals transmitted by different signal transmitting ends according to the relationship between the magnetic field harmonic component and the harmonic number, the interpolation interval of the preset harmonic number and the first preset harmonic number cutoff value;

[0112] The fourth determination module 1005 is used to determine the amplitude ratio and phase difference of the electric field components when the signal receiving end receives signals transmitted by different signal transmitting ends; determine the amplitude ratio and phase difference of the magnetic field components when the signal receiving end receives signals transmitted by different signal transmitting ends.

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

[0114] According to the position of the signal receiving end in the measurement equivalent calculation model, the coupling relationship of the field at the radial layer interface of the measurement equivalent calculation model is described by using the generalized reflection array and transmission array to determine the electric field spectrum domain and magnetic field spectrum domain when the signal receiving end receives signals transmitted by different signal transmitting ends.

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

[0116] A component representation acquisition module is used to obtain the electric field component and the magnetic field component of the Fourier series expansion of the signal receiving end in the measurement equivalent calculation model before determining the relationship between the electric field harmonic component and the harmonic number and the relationship between the magnetic field harmonic component and the harmonic number when the signal receiving end receives signals transmitted by different signal transmitting ends;

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

[0118] According to the electric field spectral domain field, the magnetic field spectral domain field, and the electric field components and the magnetic field components of the Fourier series expansion, the relationship between the electric field harmonic component and the harmonic number, and the relationship between the magnetic field harmonic component and the harmonic number when the signal receiving end receives signals transmitted by different signal transmitting ends are determined.

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

[0120] An interpolation interval determination module is used to determine the electric field component of a signal receiving end when receiving signals transmitted by different signal transmitting ends according to the relationship between the electric field harmonic component 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 component of the signal receiving end when receiving signals transmitted by different signal transmitting ends according to the relationship between the magnetic field harmonic component and the harmonic number, the interpolation interval of the preset harmonic number and the first preset harmonic number cutoff value, determine the first magnetic field component of the signal receiving end when receiving a signal transmitted by a signal transmitting end according to the relationship between the magnetic field harmonic component and the harmonic number and the first preset harmonic number cutoff value; determine the second magnetic field component of the signal receiving end when receiving the signal transmitted by the signal transmitting end according to the relationship between the magnetic field harmonic component 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; determine the relationship between the error between the second magnetic field component and the first magnetic field component and the harmonic number; and determine the interpolation interval of the preset harmonic number according to the relationship between the error between the second magnetic field component and the first magnetic field component and the harmonic number.

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

[0122] Determine a local error function according to the error between the second magnetic field component and the first magnetic field component; 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;

[0123] Determining an interpolation interval of a preset harmonic number according to a local error function;

[0124] Wherein, the local error function is as follows:

[0125]

[0126]

[0127] 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; represents the first magnetic field component; H zz represents the second magnetic field component; τ is a preset constant; P p (v) represents the n-th order interpolation polynomial; where P p (v) = a0 + a1v + a2v 2 +...+a n v n ; a0, a1, a2, ..., a n are the constant coefficients of the n-th order interpolation polynomial; n is the order of the interpolation polynomial; and p is the interpolation interval of the harmonic number.

[0128] It should be noted that although several modules of the forward computing device of the push-type dielectric scanning logging tool are mentioned in the above detailed description, this division is only exemplary and not mandatory. In fact, according to an embodiment 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 into multiple modules for embodiment.

[0129] Based on the above invention concept, Fig.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, wherein the processor 1102 implements the forward calculation method of the aforementioned push-type dielectric scanning logging instrument when executing the computer program 1103.

[0130] Based on the aforementioned inventive concept, the present invention proposes a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the forward modeling method of the aforementioned push-type dielectric scanning logging tool is implemented.

[0131] Based on the aforementioned inventive concept, the present invention proposes a computer program product, which includes a computer program. When the computer program is executed by a processor, a forward modeling method of a push-type dielectric scanning logging tool is implemented.

[0132] The forward modeling method and device of the push-type dielectric scanning logging instrument proposed in the embodiment of the present invention can solve the problems that the forward modeling efficiency of the existing three-dimensional numerical method is very low and the computer memory occupied is high; the embodiment of the present invention obtains the measurement equivalent calculation model of the push-type dielectric scanning logging instrument; wherein the measurement equivalent calculation model is obtained by equivalently treating the push-type dielectric scanning logging instrument's push arm as a cylindrical metal core, the antenna as a magnetic dipole, and the formation as a columnar formation; according to the position of the signal receiving end in the measurement equivalent calculation model, the electric field spectrum domain field and the magnetic field spectrum domain field of the signal receiving end when receiving signals transmitted by different signal transmitting ends are determined; according to the electric field spectrum domain field and the magnetic field spectrum domain field, the signal receiving end is determined when receiving signals transmitted by different signal transmitting ends. The invention relates to a method for determining a relationship between an electric field harmonic component and a harmonic number, and a relationship between a magnetic field harmonic component and a harmonic number when a signal receiving end transmits a signal; determining an electric field component when a signal receiving end receives signals transmitted by different signal transmitting ends according to the relationship between the electric field harmonic component and the harmonic number, an interpolation interval of a preset harmonic number, and a first preset harmonic number cutoff value; determining a magnetic field component when a signal receiving end receives signals transmitted by different signal transmitting ends according to the relationship between the magnetic field harmonic component and the harmonic number, an interpolation interval of a preset harmonic number, and a first preset harmonic number cutoff value; determining an amplitude ratio and a phase difference of an electric field component when the signal receiving end receives signals transmitted by different signal transmitting ends; and determining an amplitude ratio and a phase difference of a magnetic field component when the signal receiving end receives signals transmitted by different signal transmitting ends. The embodiment of the present invention provides a fast forward algorithm for the amplitude ratio and phase difference response of dielectric scanning logging, uses the measurement equivalent calculation model of the push-type dielectric scanning logging instrument as the basic model for forward calculation, and presets the interpolation interval of the harmonic number during the forward calculation process, which greatly reduces the amount of calculation, improves the forward calculation efficiency of the push-type dielectric scanning logging instrument, and reduces the occupied computer memory.

[0133] Those skilled in the art will appreciate that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented 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.

[0134] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0135] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture including an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0136] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0137] The specific embodiments described above further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above description is only a specific embodiment of the present invention and is 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 in the scope of protection of the present invention.

Claims

1. A forward calculation method for a push-type dielectric scanning logging tool, characterized in that: include: Obtaining a measurement equivalent calculation model of a push-type dielectric scanning logging tool; wherein the measurement equivalent calculation model is obtained by equating the push-type dielectric scanning logging tool's push arm to a cylindrical metal core, the antenna to a magnetic dipole, and the formation to a cylindrical formation; According to the position of the signal receiving end in the measurement equivalent calculation model, the electric field spectrum domain and the magnetic field spectrum domain of the signal receiving end when receiving signals transmitted by different signal transmitting ends are determined; Determine the relationship between the electric field harmonic component and the harmonic number, and the relationship between the magnetic field harmonic component and the harmonic number when the signal receiving end receives signals transmitted by different signal transmitting ends according to the electric field spectrum domain field and the magnetic field spectrum domain field; Determine the electric field component of the signal receiving end when receiving signals transmitted by different signal transmitting ends according to the relationship between the electric field harmonic component and the harmonic number, the interpolation interval of the preset harmonic number and the first preset harmonic number cutoff value; determine the magnetic field component of the signal receiving end when receiving signals transmitted by different signal transmitting ends according to the relationship between the magnetic field harmonic component and the harmonic number, the interpolation interval of the preset harmonic number and the first preset harmonic number cutoff value; Determine the amplitude ratio and phase difference of the electric field components when the signal receiving end receives signals transmitted by different signal transmitting ends; determine the amplitude ratio and phase difference of the magnetic field components when the signal receiving end receives signals transmitted by different signal transmitting ends.

2. The method according to claim 1, characterized in that According to the position of the signal receiving end in the measurement equivalent calculation model, the electric field spectrum domain and the magnetic field spectrum domain of the signal receiving end when receiving signals transmitted by different signal transmitting ends are determined, including: According to the position of the signal receiving end in the measurement equivalent calculation model, the coupling relationship of the field at the radial layer interface of the measurement equivalent calculation model is described by using the generalized reflection array and transmission array to determine the electric field spectrum domain and magnetic field spectrum domain when the signal receiving end receives signals transmitted by different signal transmitting ends.

3. The method according to claim 1, characterized in that Before determining the relationship between the electric field harmonic component and the harmonic number, and the relationship between the magnetic field harmonic component and the harmonic number when the signal receiving end receives signals transmitted by different signal transmitting ends, it includes: Obtaining the electric field component and the magnetic field component of the Fourier series expansion of the signal receiving end in the measurement equivalent calculation model; According to the electric field spectrum domain and the magnetic field spectrum domain, the relationship between the electric field harmonic component and the harmonic number, and the relationship between the magnetic field harmonic component and the harmonic number when the signal receiving end receives signals transmitted by different signal transmitting ends are determined, including: According to the electric field spectral domain field, the magnetic field spectral domain field, and the electric field components and the magnetic field components of the Fourier series expansion, the relationship between the electric field harmonic component and the harmonic number, and the relationship between the magnetic field harmonic component and the harmonic number when the signal receiving end receives signals transmitted by different signal transmitting ends are determined.

4. The method according to claim 1, characterized in that According to the relationship between the electric field harmonic component and the harmonic number, the interpolation interval of the preset harmonic number and the first preset harmonic number cutoff value, the electric field component of the signal receiving end when receiving the signals transmitted by different signal transmitting ends is determined, and according to the relationship between the magnetic field harmonic component and the harmonic number, the interpolation interval of the preset harmonic number and the first preset harmonic number cutoff value, the magnetic field component of the signal receiving end when receiving the signals transmitted by different signal transmitting ends is determined, before the signal receiving end receives the signals transmitted by different signal transmitting ends, the method further includes: Determine, according to the relationship between the magnetic field harmonic component and the harmonic number, and the first preset harmonic number cutoff value, a first magnetic field component when the signal receiving end receives a signal transmitted by a signal transmitting end; Determine the second magnetic field component of the signal receiving end when receiving the signal transmitted by the signal transmitting end according to the relationship between the magnetic field harmonic component 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; determining an error between the second magnetic field component and the first magnetic field component; An interpolation interval of a preset harmonic number is determined according to an error between the second magnetic field component and the first magnetic field component.

5. The method according to claim 4, characterized in that Determining an interpolation interval of a preset harmonic number according to an error between the second magnetic field component and the first magnetic field component includes: Determine a local error function according to the error between the second magnetic field component and the first magnetic field component; 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; Determining an interpolation interval of a preset harmonic number according to a local error function; Wherein, 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; represents the first magnetic field component; H zz represents the second magnetic field component; τ is a preset constant; P p (v) represents the n-th order interpolation polynomial; where P p (v) = a0 + a1v + a2v 2 +...+a n v n ; a0, a1, a2, ..., a n are the constant coefficients of the n-th order interpolation polynomial; n is the order of the interpolation polynomial; and p is the interpolation interval of the harmonic number.

6. A forward 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-type dielectric scanning logging tool's push arm to a cylindrical metal core, the antenna to a magnetic dipole, and the formation to a cylindrical formation; A first determination module is used to determine the electric field spectrum domain and the magnetic field spectrum domain of the signal receiving end when receiving signals transmitted by different signal transmitting ends according to the position of the signal receiving end in the measurement equivalent calculation model; The second determination module is used to determine the relationship between the electric field harmonic component and the harmonic number, and the relationship between the magnetic field harmonic component and the harmonic number when the signal receiving end receives signals transmitted by different signal transmitting ends according to the electric field spectrum domain field and the magnetic field spectrum domain field; A third determination module is used to determine the electric field component of the signal receiving end when receiving signals transmitted by different signal transmitting ends according to the relationship between the electric field harmonic component 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 component of the signal receiving end when receiving signals transmitted by different signal transmitting ends according to the relationship between the magnetic field harmonic component and the harmonic number, the interpolation interval of the preset harmonic number and the first preset harmonic number cutoff value; The fourth determination module is used to determine the amplitude ratio and phase difference of the electric field components when the signal receiving end receives signals transmitted by different signal transmitting ends; determine the amplitude ratio and phase difference of the magnetic field components when the signal receiving end receives signals transmitted by different signal transmitting ends.

7. The device according to claim 6, characterized in that The first determination module is specifically used for: According to the position of the signal receiving end in the measurement equivalent calculation model, the coupling relationship of the field at the radial layer interface of the measurement equivalent calculation model is described by using the generalized reflection array and transmission array to determine the electric field spectrum domain and magnetic field spectrum domain when the signal receiving end receives signals transmitted by different signal transmitting ends.

8. The device according to claim 6, characterized in that Also includes: A component representation acquisition module is used to obtain the electric field component and the magnetic field component of the Fourier series expansion of the signal receiving end in the measurement equivalent calculation model before determining the relationship between the electric field harmonic component and the harmonic number and the relationship between the magnetic field harmonic component and the harmonic number when the signal receiving end receives signals transmitted by different signal transmitting ends; The second determination module is specifically used for: According to the electric field spectral domain field, the magnetic field spectral domain field, and the electric field components and the magnetic field components of the Fourier series expansion, the relationship between the electric field harmonic component and the harmonic number, and the relationship between the magnetic field harmonic component and the harmonic number when the signal receiving end receives signals transmitted by different signal transmitting ends are determined.

9. The device according to claim 6, characterized in that Also includes: An interpolation interval determination module is used to determine the electric field component of the signal receiving end when receiving signals transmitted by different signal transmitting ends according to the relationship between the electric field harmonic component 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 component of the signal receiving end when receiving signals transmitted by different signal transmitting ends according to the relationship between the magnetic field harmonic component and the harmonic number, the interpolation interval of the preset harmonic number and the first preset harmonic number cutoff value, determine the first magnetic field component of the signal receiving end when receiving a signal transmitted by a signal transmitting end according to the relationship between the magnetic field harmonic component and the harmonic number and the first preset harmonic number cutoff value; determine the second magnetic field component of the signal receiving end when receiving the signal transmitted by the signal transmitting end according to the relationship between the magnetic field harmonic component 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; Determine the relationship between the error between the second magnetic field component and the first magnetic field component and the harmonic number; and determine the interpolation interval of the preset 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.

10. The device according to claim 9, characterized in that The interpolation interval determination module is specifically used for: Determine a local error function according to the error between the second magnetic field component and the first magnetic field component; 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; Determining an interpolation interval of a preset harmonic number according to a local error function; Wherein, 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; represents the first magnetic field component; H zz represents the second magnetic field component; τ is a preset constant; P p (v) represents the n-th order interpolation polynomial; where P p (v) = a0 + a1v + a2v 2 +...+a n v n ; a0, a1, a2, ..., a n are the constant coefficients of the n-th order interpolation polynomial; n is the order of the interpolation polynomial; and 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, the method according to any one of claims 1 to 5 is implemented.

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

13. A computer program product, characterized in that The computer program product comprises a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 5 is implemented.

Images