Logging Interpretation Method, Device and Storage Medium Based on Array Induction

Abstract

An embodiment of the present application provides a logging interpretation method, device, and storage medium based on array induction. Logging data of a horizontal well to be measured is collected through array induction, and the logging data is subjected to environmental correction to obtain preliminary logging data. There is an association relationship between the preliminary logging data and the measurement frequency. The logging data of the horizontal well to be measured includes initial conductivity. The preliminary logging data is subjected to inversion correction to obtain the logging result of the horizontal well to be measured after inversion correction. The technical solution provided by the present application can obtain the logging result of the horizontal well to be measured by performing inversion correction on the preliminary logging data after environmental correction, thereby improving the accuracy of array induction logging.

Description

Technical Field

[0001] The present application relates to the field of logging technology, and in particular, to a logging interpretation method, device, and storage medium based on array induction. Background Art

[0002] In fields such as geological exploration and well logging, array induction is usually used for logging operations. Among them, during measurement, environmental correction is generally performed on the acquired logging data, and software focusing and resolution matching operations are carried out using software to obtain logging results.

[0003] Currently, array induction is usually used for logging, mainly by using multiple sub-arrays and multi-frequency measurements, which can provide resistivity curves with different detection depths for users. The processing of the data measured by the array induction logging instrument mainly includes two modules, an environmental correction module and a software focusing module. Among them, the environmental correction module can perform environmental correction on the logging data, including skin effect correction and borehole correction. The software focusing module can perform software focusing and longitudinal resolution matching. During the operation process, it is often aimed at the measurement of vertical wells on a horizontal ground, that is, the instrument is usually operated with its axis perpendicular to the formation interface. The acquired logging data is processed by the environmental correction module to obtain preliminarily processed data, and the preliminarily processed data is further processed by the software focusing module to obtain resistivity curves with different depths of longitudinal resolution, so as to interpret the situation of the measured well.

[0004] However, the formation distribution may not be longitudinally distributed on a horizontal plane. For example, the well to be measured may be a horizontal well. If the current array induction method is used to measure a horizontal well, there is a large error between the obtained logging result and the true result of the horizontal well, resulting in a low accuracy of the array induction logging result. Summary of the Invention

[0005] The embodiments of the present application provide a logging interpretation method, device, and storage medium based on array induction, which can obtain the logging result of a horizontal well, thereby improving the accuracy of array induction logging.

[0006] In a first aspect, the embodiments of the present application provide a logging interpretation method based on array induction, and the logging interpretation method based on array induction includes:

[0007] Collect logging data of a horizontal well to be measured through array induction, and perform environmental correction on the logging data to obtain preliminary logging data. The preliminary logging data has an association relationship with the measurement frequency, and the logging data of the horizontal well to be measured includes initial conductivity;

[0008] Perform inversion correction on the preliminary logging data to obtain the logging results of the horizontal well to be measured after inversion correction.

[0009] In a possible implementation, the performing inversion correction on the preliminary logging data to obtain the logging results after inversion correction includes:

[0010] Determine the constraint conditions of the logging data.

[0011] Perform inversion correction on the preliminary logging data through a preset inversion formula and the constraint conditions to obtain the logging results after inversion correction.

[0012] The preset inversion formula is:

[0013]

[0014] where δ is the error minimization objective function, is the preliminary logging data, The sub-array response vector calculated from the initial value of the preliminary logging data The initial value Constrants(1) represents the first constraint condition, Constrants(2) represents the second constraint condition, Constrants(N) represents the Nth constraint condition, and N is a natural number.

[0015] In a possible implementation, the constraint conditions at least include: the mutual constraint condition of the upper surrounding rock conductivity, the mutual constraint condition of the target layer conductivity, the mutual constraint condition of the lower surrounding rock conductivity, the constraint condition of the upper surrounding rock formation interface position, the constraint condition of the lower surrounding rock formation interface position, the prior value constraint condition of the upper surrounding rock formation conductivity, the prior value constraint condition of the target layer conductivity, the prior value constraint condition of the lower surrounding rock formation conductivity, and the prior value constraint condition of the upper and lower surrounding rock formation interface positions and formation thickness.

[0016] In a possible implementation, the constraint conditions include the mutual constraint condition of the upper surrounding rock conductivity; the determining the constraint conditions of the logging data includes:

[0017] Calculate the mutual constraint condition of the upper surrounding rock conductivity through the first constraint formula.

[0018] The first constraint formula is:

[0019]

[0020] where Npoint represents the number of inversion points taken, is the mutual constraint coefficient between the resistivities of the upper surrounding rock formations inverted at the kth and lth points, is the conductivity of the upper surrounding rock formation corresponding to the kth point, is the electrical conductivity of the upper surrounding rock formation corresponding to the l-th point.

[0021] The constraint conditions include the mutual constraint conditions of the conductivity of the target layer; determining the constraint conditions of the logging data includes:

[0022] Calculate the mutual constraint conditions of the conductivity of the target layer through the second constraint formula.

[0023] The second constraint formula is:

[0024]

[0025] where is the mutual constraint coefficient between the resistivities of the target layer inverted at the k-th point and the l-th point, is the resistivity of the target layer corresponding to the k-th point, is the resistivity of the target layer corresponding to the l-th point.

[0026] The constraint conditions include the mutual constraint conditions of the conductivity of the lower surrounding rock; determining the constraint conditions of the logging data includes:

[0027] Calculate the mutual constraint conditions of the conductivity of the lower surrounding rock through the third constraint formula.

[0028] The third constraint formula is:

[0029]

[0030] where is the mutual constraint coefficient between the resistivities of the lower surrounding rock formation inverted at the k-th point and the l-th point, is the resistivity of the lower surrounding rock formation corresponding to the k-th point, is the resistivity of the lower surrounding rock formation corresponding to the l-th point.

[0031] The constraint conditions include the constraint conditions of the position of the upper surrounding rock formation interface; determining the constraint conditions of the logging data includes:

[0032] Calculate the constraint conditions of the position of the upper surrounding rock formation interface through the fourth constraint formula.

[0033] The fourth constraint formula is:

[0034]

[0035] where is the mutual constraint coefficient between the positions of the upper surrounding rock interfaces inverted at the k-th point and the l-th point, is the position of the upper surrounding rock interface corresponding to the k-th point, is the position of the upper surrounding rock interface corresponding to the l-th point.

[0036] The constraint conditions include the constraint condition for the position of the lower surrounding rock formation interface; determining the constraint conditions for the logging data includes:

[0037] Calculate the constraint condition for the position of the lower surrounding rock formation interface through the fifth constraint formula.

[0038] The fifth constraint formula is:

[0039]

[0040] Wherein, is the mutual constraint coefficient between the positions of the lower surrounding rock interfaces inverted at the k-th point and the l-th point, is the position of the lower surrounding rock interface corresponding to the k-th point, is the position of the lower surrounding rock interface corresponding to the l-th point.

[0041] The constraint conditions include the prior value constraint condition for the conductivity of the upper surrounding rock formation; determining the constraint conditions for the logging data includes:

[0042] Calculate the prior value constraint condition for the conductivity of the upper surrounding rock formation through the sixth constraint formula.

[0043] The sixth constraint formula is:

[0044]

[0045] Wherein, is the constraint coefficient between the resistivity of the upper surrounding rock formation inverted at the l-th point and the prior value, is the resistivity of the upper surrounding rock formation corresponding to the l-th point, is the reference value of the resistivity of the upper surrounding rock formation corresponding to the l-th point, and this reference value can be an existing prior value.

[0046] The constraint conditions include the prior value constraint condition for the conductivity of the target layer; determining the constraint conditions for the logging data includes:

[0047] Calculate the prior value constraint condition for the conductivity of the target layer through the seventh constraint formula.

[0048] The seventh constraint formula is:

[0049]

[0050] Wherein, is the constraint coefficient between the resistivity of the target layer inverted at the l-th point and the prior value, is the resistivity of the target layer corresponding to the l-th point, is the reference value of the resistivity of the target layer corresponding to the l-th point, and this reference value can be an existing prior value.

[0051] The constraint conditions include the prior value constraint condition of the conductivity of the lower surrounding rock formation; determining the constraint conditions of the logging data includes:

[0052] Calculate the prior value constraint condition of the conductivity of the lower surrounding rock formation through the eighth constraint formula.

[0053] The eighth constraint formula is:

[0054]

[0055] Wherein, is the constraint coefficient between the resistivity of the lower surrounding rock formation inverted at the l-th point and the prior value, is the resistivity of the lower surrounding rock formation corresponding to the l-th point, is the reference value of the resistivity of the lower surrounding rock formation corresponding to the l-th point, and this reference value can be the existing prior value.

[0056] The constraint conditions include the prior value constraint conditions of the positions of the upper and lower surrounding rock formation interfaces and the formation thickness; determining the constraint conditions of the logging data includes:

[0057] Use the ninth constraint formula to calculate the prior value constraint conditions of the positions of the upper and lower surrounding rock formation interfaces and the formation thickness.

[0058] The ninth constraint formula is:

[0059]

[0060] Wherein, is the constraint coefficient between the reference value of the lower surrounding rock formation interface and the formation thickness inverted at the l-th point, is the position of the upper surrounding rock interface corresponding to the l-th point, is the position of the lower surrounding rock interface corresponding to the l-th point, H is the reference value of the formation thickness corresponding to the l-th point, and the reference value is the prior value that has been collected.

[0061] In a possible implementation manner, before calculating each constraint condition in the constraint conditions through multiple constraint condition formulas respectively, the method further includes:

[0062] Collect the reference value of the conductivity of the upper surrounding rock formation, the reference value of the conductivity of the target layer, the reference value of the conductivity of the lower surrounding rock formation, and the reference value of the formation thickness from the adjacent wellbore of the current logging wellbore.

[0063] In a possible implementation manner, the environmental correction of the logging data to obtain preliminary logging data includes:

[0064] Perform borehole correction on the logging data, and determine the corrected logging data as the preliminary logging data;

[0065] Alternatively,

[0066] perform skin effect correction and borehole correction on the logging data to obtain preliminary corrected data; perform a frequency conversion operation on the preliminary corrected data, and determine the converted preliminary corrected data as the preliminary logging data.

[0067] In a possible implementation, the performing a frequency conversion operation on the preliminary corrected data includes:

[0068] Perform a frequency conversion operation on the preliminary corrected data through the following preset conversion formula:

[0069]

[0070] where: k = iωμσ, σ is the preliminary corrected data, L R is the main coil spacing, L B is the shield coil spacing, σ aR is the conductivity corresponding to the corrected sub-array, μ is the magnetic permeability, and i is the imaginary unit.

[0071] where each sub-array corresponds to a σ aR , and determine the values of all sub-arrays as the preliminary logging data.

[0072] In a second aspect, an embodiment of the present application provides a logging interpretation device based on array induction. The logging interpretation device based on array induction includes:

[0073] An environment correction module 1101, configured to collect logging data of a horizontal well to be measured through array induction, and perform environment correction on the logging data to obtain preliminary logging data. The preliminary logging data has an associated relationship with the measurement frequency. The logging data of the horizontal well to be measured includes an initial conductivity;

[0074] An inversion correction module 1102, configured to perform inversion correction on the preliminary logging data to obtain a logging result of the horizontal well to be measured after inversion correction.

[0075] In a third aspect, an embodiment of the present application further provides an electronic device. The electronic device includes: a processor, and a memory communicatively connected to the processor;

[0076] The memory stores computer-executable instructions;

[0077] The processor executes the computer-executable instructions stored in the memory to implement the logging interpretation method based on array induction in any possible implementation manner of the first aspect above.

[0078] Fourthly, an embodiment of the present application further provides a computer-readable storage medium storing computer-executable instructions, which, when executed by a processor, implement the logging interpretation method based on array induction described in any possible implementation manner of the first aspect above.

[0079] Fifthly, an embodiment of the present application further provides a computer program product including a computer program, which, when executed by a processor, implements the logging interpretation method based on array induction described in any possible implementation manner of the first aspect above.

[0080] It can be seen that an embodiment of the present application provides a logging interpretation method, device, and storage medium based on array induction. Logging data of a horizontal well to be measured is collected through array induction, and the logging data is subjected to environmental correction to obtain preliminary logging data. There is an association between the preliminary logging data and the measurement frequency. The logging data of the horizontal well to be measured includes initial conductivity. The preliminary logging data is subjected to inversion correction to obtain the logging result of the horizontal well to be measured after inversion correction. The technical solution provided by the present application performs environmental correction on the logging data of the horizontal well to be measured and performs inversion correction on the preliminary logging data after environmental correction, so as to obtain the logging result of the horizontal well to be measured, making the obtained logging result closer to the true result of the horizontal well, enabling the array induction logging method to be applicable to the measurement of horizontal wells, and thus improving the accuracy of array induction logging. Description of the Drawings

[0081] Figure 1 is the principle of measuring formation conductivity by a dual-coil array provided by an embodiment of the present application;

[0082] Figure 2 is a schematic diagram of the principle of measuring formation conductivity by a three-coil array with a shielded receiving coil provided by an embodiment of the present application;

[0083] Figure 3 is a schematic diagram of the position of an array induction coil placed horizontally provided by an embodiment of the present application;

[0084] Figure 4 is a schematic diagram of the longitudinal differential geometric factor of two sub-arrays provided by an embodiment of the present application;

[0085] Figure 5 is a schematic diagram of a formation model of a vertical well and a horizontal formation provided by an embodiment of the present application;

[0086] Figure 6 is a schematic diagram of the radial integral geometric factor of two sub-arrays provided by an embodiment of the present application;

[0087] Figure 7Schematic diagram of a formation model of a horizontal well in a horizontal formation provided by an embodiment of the present application;

[0088] Figure 8 Schematic flow chart of a logging interpretation method based on array induction provided by an embodiment of the present application;

[0089] Figure 9 Schematic flow chart of another logging interpretation method based on array induction provided by an embodiment of the present application;

[0090] Figure 10 Schematic flow chart of an inversion correction method provided by an embodiment of the present application;

[0091] Figure 11 Schematic diagram of the structure of a logging interpretation device 110 based on array induction provided by an embodiment of the present application;

[0092] Figure 12 Schematic diagram of the structure of an electronic device provided by the present application.

[0093] Through the above-mentioned drawings, specific embodiments of the present disclosure have been shown, and there will be more detailed descriptions hereinafter. These drawings and the written description are not intended to limit the scope of the concept of the present disclosure in any way, but to illustrate the concept of the present disclosure to those skilled in the art by referring to specific embodiments. Specific Embodiments

[0094] Here, exemplary embodiments will be described in detail, and examples are shown in the drawings. When the following description refers to the drawings, unless otherwise indicated, the same numbers in different drawings represent the same or similar elements. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with the present disclosure. On the contrary, they are merely examples of devices and methods consistent with some aspects of the present disclosure as detailed in the appended claims.

[0095] In the embodiments of the present application, "at least one" means one or more, and "a plurality" means two or more. "And / or" describes the association relationship of associated objects and indicates that three relationships may exist. For example, A and / or B may represent: A exists alone, A and B exist simultaneously, and B exists alone. Here, A and B may be singular or plural. In the written description of the present application, the character " / " generally represents an "or" relationship between the associated objects before and after.

[0096] The technical solutions provided by the embodiments of the present application can be applied to the scenario of logging. During the logging process, the method of array induction logging can be used to interpret the situation of the well to be logged. Hereinafter, the principle of array induction logging will be explained. Figure 1This is the principle of measuring formation conductivity using a dual-coil system provided by an embodiment of the present application. According to Figure 1 As shown, the transmitting and receiving coils are on the instrument rod. When the transmitting coil emits electromagnetic waves, a primary field will be generated in the surrounding medium. This primary field will generate eddy current rings in the formation, and these eddy current rings will generate a secondary field. The receiving coil will receive the generated primary field and secondary field. Since the received secondary field is related to the formation conductivity, the formation conductivity can be calculated using the received signal, that is, the secondary field signal. However, the primary field is much larger than the secondary field. Therefore, the logging instrument must overcome the primary field to obtain the secondary field.

[0097] Furthermore, for eliminating the primary field, refer to Figure 2 As shown, Figure 2 This is a schematic diagram of the principle of measuring formation conductivity using a three-coil system with a shielded receiving coil provided by an embodiment of the present application. In Figure 2 there are two receiving coils, namely the outermost main receiving coil and the innermost auxiliary receiving coil, that is, the shielded receiving coil. The long winding directions of the main receiving coil and the shielded receiving coil are opposite, and the connecting coils are in series. This enables the purpose of eliminating the primary field to be achieved by appropriately adjusting the number of turns of the two receiving coils and the distance from the transmitting coil. Currently, array induction logging usually uses an array induction logging instrument composed of three coils as shown in Figure 2 As shown.

[0098] Since the wellbore direction is different during different measurement processes, as shown in Figure 3 As shown, Figure 3 This is a schematic diagram of the position of an array induction coil placed horizontally provided by an embodiment of the present application. In Figure 3 the transmitting coil is represented by T, there is a sub-array of N. The shielded receiving coil of the first sub-array is represented by RB-1, the main receiving coil of the first sub-array is represented by RM-1, and so on. The shielded receiving coil of the Nth sub-array is represented by RB-N, and the main receiving coil of the Nth sub-array is represented by RM-N.

[0099] Exemplarily, Figure 3 the longitudinal geometric factors of two sub-arrays can be referred to Figure 4 As shown, Figure 4 This is a schematic diagram of the longitudinal differential geometric factors of two sub-arrays provided by an embodiment of the present application. In Figure 4 the solid line represents the geometric factor when the coil spacing from the main receiving coil to the transmitting coil is 1 meter, and the dashed line represents the geometric factor when the coil spacing from the main receiving coil to the transmitting coil is 0.5 meter. According to Figure 4It can be known that for the sub-array with a short coil spacing from the main receiving line to the transmitting coil, the main contribution area in the longitudinal direction is relatively narrow, but the contribution value is relatively large, that is, the longitudinal resolution is relatively high. For the sub-array with a long coil spacing from the main receiving line to the transmitting coil, the main contribution area in the longitudinal direction is relatively large, but the contribution value is relatively small, that is, the longitudinal resolution is relatively low. It can be known that the basic principle of induction logging is: the larger the coil spacing, the lower the longitudinal resolution. It must be pointed out here that the longitudinal resolution refers to the resolution along the axis of the array induction logging tool.

[0100] Furthermore, Figure 5 FIG. is a schematic diagram of a formation model of a vertical well and a horizontal formation provided by an embodiment of the present application. Figure 5 As shown, it is a three-layer formation model. The wellbore of the vertical well is perpendicular to the formation interface, that is, the wellbore trajectory of the vertical well is perpendicular to the first formation interface, the second formation interface and the third formation interface. According to the above principle of array induction, it can be known that the formation with a two-dimensional longitudinal distribution can be measured by the method of array induction. Specifically, the logging data of the array induction logging tool can be processed through an environmental correction module and a software focusing module to obtain the logging result.

[0101] Figure 6 FIG. is a schematic diagram of the radial integral geometric factor of two sub-arrays provided by an embodiment of the present application. In Figure 6 , the dashed line represents the geometric factor when the coil spacing from the main receiving line to the transmitting coil is 1 m, and the solid line represents the geometric factor when the coil spacing from the main receiving line to the transmitting coil is 0.5 m. The 50% contribution line represents that the intersection position of the integral geometric factor curve and this line is a point where the response of a sub-array accounts for 50% inside and outside this point. According to Figure 6 It can be known that if the position of the 50% contribution line is used to determine the radial detection range of a sub-array, the radial detection range of the sub-array with a main coil spacing of 0.5 m is about 0.4 m, while the radial detection range of the sub-array with a main coil spacing of 1 m is about 1.6 m. It can be known that the basic principle of induction logging is: the larger the coil spacing, the larger the radial detection range. Among them, the radial detection depth refers to the radius of the cylinder centered on the axis of the array induction logging tool.

[0102] However, during the logging process, the well to be logged may be a horizontal well. Specifically, see Figure 7 as shown in Figure 7 FIG. is a schematic diagram of a formation model of a horizontal well in a horizontal formation provided by an embodiment of the present application. Figure 7In the second formation interface, the wellbore trajectory of the horizontal well is basically parallel to the second formation interface. At this time, if the existing logging method is used to measure the horizontal well, the axis of the array induction logging instrument also needs to be basically parallel to the second formation interface. Therefore, when the existing array induction logging method is used to measure the horizontal well, there is no longitudinal resolution involved in the vertical well measurement, resulting in a large error between the logging result and the true result of the horizontal well, thus leading to a low accuracy of the array induction logging.

[0103] Considering that the reason for the low accuracy of the array induction logging is that the existing array induction method can only be accurately applied to the case of vertical wells, and the existing array induction logging processes the logging data of vertical wells through an environmental correction module and a software focusing module to obtain the measurement result. Among them, the environmental correction includes skin effect correction and borehole correction, and the output result is independent of the frequency. Therefore, each sub-array during logging is subjected to environmental correction to obtain a result related to the frequency, and then through inversion correction, the data after borehole correction is processed to obtain the logging result of the horizontal well, which can improve the accuracy of the array induction logging.

[0104] Next, the logging interpretation method based on array induction provided by the present application will be described in detail through specific embodiments. It can be understood that the following specific embodiments can be combined with each other, and the same or similar concepts or processes may not be repeated in some embodiments.

[0105] Figure 8 It is a schematic flowchart of a logging interpretation method based on array induction provided by an embodiment of the present application. This logging interpretation method based on array induction can be executed by software and / or a hardware device. For example, the hardware device can be a logging interpretation device based on array induction, and this logging interpretation device based on array induction can be a terminal or a processing chip in the terminal. Exemplarily, please refer to Figure 2 As shown, this logging interpretation method based on array induction can include:

[0106] S801. Collect logging data of the horizontal well to be measured through array induction, and perform environmental correction on the logging data to obtain preliminary logging data, and the preliminary logging data has an associated relationship with the measurement frequency.

[0107] The logging data of the horizontal well to be measured includes the initial conductivity. Exemplarily, the logging data of the horizontal well to be measured can also include other data, and the present application embodiment does not make any limitation on the specific logging data.

[0108] When performing environmental correction on the logging data to obtain preliminary logging data, there are the following two possible implementation manners:

[0109] In a possible implementation, borehole correction can be performed on the logging data, and the corrected logging data is determined as the preliminary logging data. Since the logging data after borehole correction has each frequency of each subarray corresponding to the corrected conductivity, that is, the conductivity obtained after borehole correction is frequency-dependent.

[0110] In another possible implementation, skin effect correction and borehole correction can be performed on the logging data to obtain preliminary correction data; a frequency conversion operation is performed on the preliminary correction data, and the converted preliminary correction data is determined as the preliminary logging data. For example, the preliminary correction data is the preliminary correction data.

[0111] In the embodiments of the present application, by only performing borehole correction, or performing a frequency conversion on the correction data after skin effect correction and borehole correction, the preliminary logging data is obtained, so that the obtained preliminary logging data is all frequency-dependent, which is convenient for inversion correction and further improves the accuracy of array induction logging.

[0112] For example, when performing a frequency conversion operation on the preliminary correction data to obtain the preliminary logging data, the following formula (1), that is, a preset conversion formula, can be used to perform a frequency conversion operation on the data after environmental correction:

[0113]

[0114] where: k = iωμσ, σ is the preliminary correction data, L R is the main coil spacing, L B is the shielding coil spacing, σ aR is the conductivity corresponding to the corrected subarray, μ is the magnetic permeability, μ = 2πF0, i is the imaginary unit. Among them, each subarray corresponds to a σ aR , and the values of all subarrays are determined as the preliminary logging data.

[0115] According to the above formula (1), it can be seen that the conductivity corresponding to the corrected subarray is frequency-dependent.

[0116] In the embodiments of the present application, by performing a frequency conversion operation on the preliminary correction data to obtain the preliminary logging data, that is, obtaining that the conductivity corresponding to the corrected subarray is frequency-dependent, so as to ensure that the conductivity obtained after skin effect and borehole correction is frequency-dependent, which is convenient for inversion correction and greatly improves the accuracy of array induction logging.

[0117] S802. Perform inversion correction on the preliminary logging data to obtain the logging result of the horizontal well to be measured after inversion correction.

[0118] For example, when performing inversion correction on preliminary logging data to obtain the logging results after inversion correction, the constraint conditions of the logging data can be determined; the preliminary logging data is inversely corrected through a preset inversion formula and the constraint conditions to obtain the logging results after inversion correction.

[0119] Among them, the preset inversion formula can be the following formula (2):

[0120]

[0121] Among them, δ is the error minimization objective function, is the preliminary logging data, The sub-array response vector calculated from the initial value of the preliminary logging data, the initial value Constrants(1) represents the first constraint condition, Constrants(2) represents the second constraint condition, Constrants(N) represents the Nth constraint condition, and N is a natural number. Among them, can be the array response vector after environmental correction.

[0122] For example, the number of the above-mentioned constraint conditions can be determined according to the convergence of the difference between the array response vector that can be environmentally corrected and the sub-array response vector calculated from the initial value of each sub-array preliminary logging data during the inversion process. The embodiments of the present application do not make specific limitations on this.

[0123] In the embodiments of the present application, inversely correcting the preliminary logging data according to the preset inversion formula can obtain accurate logging results of the horizontal well to be measured, thereby improving the accuracy of array induction logging.

[0124] For example, the constraint conditions at least include: the mutual constraint condition of the upper surrounding rock conductivity, the mutual constraint condition of the target layer conductivity, the mutual constraint condition of the lower surrounding rock conductivity, the constraint condition of the upper surrounding rock formation interface position, the constraint condition of the lower surrounding rock formation interface position, the prior value constraint condition of the upper surrounding rock formation conductivity, the prior value constraint condition of the target layer conductivity, the prior value constraint condition of the lower surrounding rock formation conductivity, and the prior value constraint conditions of the upper and lower surrounding rock formation interface positions and formation thickness.

[0125] In the embodiments of the present application, inversely correcting the logging data through each constraint condition makes the obtained logging results more in line with the true value of the horizontal well, thereby improving the accuracy of array induction logging.

[0126] For example, when performing one-dimensional inversion correction, the inversion parameters may be the conductivity of the formation where the horizontal well is located, the conductivities of the upper and lower surrounding rocks, and the formation interfaces of the upper and lower surrounding rocks. When performing inversion, the de-electrification inversion method can be used, that is, point-by-point inversion. When performing inversion at each point, the inversion parameters between adjacent two points need to be mutually constrained.

[0127] For example, the constraint conditions may include the mutual constraint conditions of the conductivity of the upper surrounding rock. In this embodiment, when determining the constraint conditions of the logging data, the mutual constraint conditions of the conductivity of the upper surrounding rock can be calculated by the first constraint formula, and the first constraint formula is the following formula (3):

[0128]

[0129] In formula (3), Npoint represents the number of inversion points taken, is the mutual constraint coefficient between the resistivity of the upper surrounding rock formation inverted at the k-th point and the l-th point, is the conductivity of the upper surrounding rock formation corresponding to the k-th point, is the conductivity of the upper surrounding rock formation corresponding to the l-th point.

[0130] The constraint conditions may also include the mutual constraint conditions of the conductivity of the target layer, and the mutual constraint conditions of the conductivity of the target layer can be calculated by the second constraint formula; the second constraint formula is the following formula (4):

[0131]

[0132] In formula (4), is the mutual constraint coefficient between the resistivity of the target layer inverted at the k-th point and the l-th point, is the resistivity of the target layer corresponding to the k-th point, is the resistivity of the target layer corresponding to the l-th point.

[0133] The constraint conditions may also include the mutual constraint conditions of the conductivity of the lower surrounding rock, and the mutual constraint conditions of the conductivity of the lower surrounding rock can be calculated by the third constraint formula, and the third constraint formula is the following formula (5):

[0134]

[0135] In formula (5), is the mutual constraint coefficient between the resistivity of the lower surrounding rock formation inverted at the k-th point and the l-th point, is the resistivity of the lower surrounding rock formation corresponding to the k-th point, is the resistivity of the lower surrounding rock formation corresponding to the l-th point.

[0136] The constraint conditions may also include the constraint condition for the position of the upper surrounding rock formation interface, and the constraint condition for the position of the upper surrounding rock formation interface can be calculated by the fourth constraint formula. The fourth constraint formula is the following formula (6):

[0137]

[0138] In formula (6), is the mutual constraint coefficient between the inverted positions of the upper surrounding rock interfaces at the k-th point and the l-th point, is the position of the upper surrounding rock interface corresponding to the k-th point, is the position of the upper surrounding rock interface corresponding to the l-th point.

[0139] The constraint conditions may also include the constraint condition for the position of the lower surrounding rock formation interface, and the constraint condition for the position of the lower surrounding rock formation interface can be calculated by the fifth constraint formula. The fifth constraint formula is the following formula (7):

[0140]

[0141] In formula (7), is the mutual constraint coefficient between the inverted positions of the lower surrounding rock interfaces at the k-th point and the l-th point, is the position of the lower surrounding rock interface corresponding to the k-th point, is the position of the lower surrounding rock interface corresponding to the l-th point.

[0142] The constraint conditions may also include the prior value constraint condition for the conductivity of the upper surrounding rock formation, and the prior value constraint condition for the conductivity of the upper surrounding rock formation can be calculated by the sixth constraint formula. The sixth constraint formula is the following formula (8):

[0143]

[0144] In formula (8), is the constraint coefficient between the resistivity of the upper surrounding rock formation inverted at the l-th point and the prior value, is the resistivity of the upper surrounding rock formation corresponding to the l-th point, is the reference value of the resistivity of the upper surrounding rock formation corresponding to the l-th point, and this reference value can be the existing prior value.

[0145] The constraint conditions may also include the prior value constraint condition for the conductivity of the target layer, and the prior value constraint condition for the conductivity of the target layer can be calculated by the seventh constraint formula. The seventh constraint formula is the following formula (9):

[0146]

[0147] In formula (9), is the constraint coefficient between the resistivity of the target layer inverted at the l-th point and the prior value, is the resistivity of the target layer corresponding to the l-th point, is the reference value of the resistivity of the target layer corresponding to the l-th point, and this reference value can be an existing prior value.

[0148] The constraint condition can also include the prior value constraint condition of the conductivity of the lower surrounding rock formation. The prior value constraint condition of the conductivity of the lower surrounding rock formation is calculated by the eighth constraint formula, and the eighth constraint formula is the following formula (10):

[0149]

[0150] In formula (10), is the constraint coefficient between the resistivity of the lower surrounding rock formation inverted at the l-th point and the prior value, is the resistivity of the lower surrounding rock formation corresponding to the l-th point, is the reference value of the resistivity of the lower surrounding rock formation corresponding to the l-th point, and this reference value can be an existing prior value.

[0151] The constraint condition can also include the prior value constraint conditions of the positions of the upper and lower surrounding rock formation interfaces and the formation thickness. The ninth constraint formula can be used to calculate the prior value constraint conditions of the positions of the upper and lower surrounding rock formation interfaces and the formation thickness, and the ninth constraint formula is the following formula (11):

[0152]

[0153] In formula (11), is the constraint coefficient between the inverted lower surrounding rock formation interface at the l-th point and the reference value of the formation thickness, is the position of the upper surrounding rock interface corresponding to the l-th point, is the position of the lower surrounding rock interface corresponding to the l-th point, H is the reference value of the formation thickness corresponding to the l-th point, and the reference value is the prior value that has been collected.

[0154] In the embodiments of the present application, inversion correction is performed through multiple constraint conditions, so that the obtained logging results are more in line with the actual situation of the horizontal well, and the problem of low accuracy of the measurement results obtained by using the method of array induction logging to measure the horizontal well is solved, thereby improving the accuracy of the array induction logging.

[0155] Before calculating each constraint condition in the constraint condition through multiple constraint condition formulas respectively, the reference value of the conductivity of the upper surrounding rock formation, the reference value of the conductivity of the target layer, the reference value of the conductivity of the lower surrounding rock formation, and the reference value of the formation thickness can also be collected from the side wellbore of the current logging wellbore to make the result of the inversion correction more accurate.

[0156] It can be seen that the logging interpretation method based on array induction provided by the embodiments of the present application collects logging data of the horizontal well to be measured through array induction, performs environmental correction on the logging data to obtain preliminary logging data, and there is a correlation between the preliminary logging data and the measurement frequency. The logging data of the horizontal well to be measured includes the initial conductivity; performs inversion correction on the preliminary logging data to obtain the logging result of the horizontal well to be measured after inversion correction. The technical solution provided by the embodiments of the present application obtains preliminary logging data that has a correlation with the measurement frequency, and obtains logging data through inversion correction, avoiding the problem of low accuracy of the logging results obtained when using the existing array induction logging processing method to measure horizontal wells, and being able to obtain logging results that are more in line with the true values of horizontal wells, thereby improving the accuracy of the array induction logging processing method.

[0157] The logging interpretation method based on array induction provided by the present application will be described in detail below. Figure 9 It is a schematic flowchart of another logging interpretation method based on array induction provided by the embodiments of the present application. According to Figure 9 As shown, when obtaining the logging result of a horizontal well, environmental correction can be first performed on the logging data measured by the array induction logging instrument, and there are two methods for environmental correction. Method 1: Perform skin effect and borehole correction on the logging data to convert the logging data into conductivity related to frequency. Method 2: Perform borehole correction on the logging data, which can also convert the logging data into conductivity related to frequency. Further, perform inversion correction on the conductivity after environmental correction, specifically invert formation parameters to obtain the logging result. Among them, the methods of environmental correction and inversion correction can refer to those described in the above embodiments, and the embodiments of the present application do not make any limitations on this.

[0158] Furthermore, the inversion correction in the embodiments of the present application can be referred to Figure 10 As shown, Figure 10 It is a schematic flowchart of an inversion correction method provided by the embodiments of the present application. According to Figure 10 As shown, the inversion correction steps are as follows:

[0159] Step 1: Select the initial values of the formation conductivity corresponding to each sub-array and the formation dip data. Step 2: Calculate the response values of each three-coil system one-dimensionally. Step 3: Obtain the preliminary logging data after environmental correction of the sub-array, and calculate the residuals between the preliminary logging data and the response values of each three-coil system in Step 2. For example, the residuals can be expressed as Step 4: Determine whether the residuals satisfy convergence. If they satisfy convergence, output the result; if they do not satisfy convergence, execute the following Step 5. Step 5: Calculate the Jacobian matrix. Step 6: Solve the Jacobian matrix to obtain new initial values. Step 7: Update the initial values, and execute the above Steps 1 to 4 according to the updated initial values. Finally, output the result.

[0160] In the embodiments of the present application, through inversion correction, the obtained logging results are made more accurate.

[0161] Figure 11 The following is a schematic structural diagram of a logging interpretation device 110 based on array induction provided by the embodiments of the present application. For example, please refer to Figure 11 As shown, the logging interpretation device 110 based on array induction may include:

[0162] An environmental correction module 1101, configured to collect logging data of a horizontal well to be measured through array induction, perform environmental correction on the logging data to obtain preliminary logging data, and there is an associated relationship between the preliminary logging data and the measurement frequency. The logging data of the horizontal well to be measured includes an initial conductivity;

[0163] An inversion correction module 1102, configured to perform inversion correction on the preliminary logging data to obtain the logging result of the horizontal well to be measured after inversion correction.

[0164] Optionally, the inversion correction module 1102 is specifically configured to determine the constraint conditions of the logging data; perform inversion correction on the preliminary logging data through a preset inversion formula and the constraint conditions to obtain the logging result after inversion correction; the preset inversion formula is:

[0165]

[0166] where δ is the error minimization objective function, is the preliminary logging data, The sub-array response vector calculated from the initial value of the preliminary logging data, the initial value Constrants(1) represents the first constraint condition, Constrants(2) represents the second constraint condition, Constrants(N) represents the Nth constraint condition, and N is a natural number.

[0167] Optionally, the constraint conditions at least include: the mutual constraint condition of the upper formation conductivity, the mutual constraint condition of the target layer conductivity, the mutual constraint condition of the lower formation conductivity, the constraint condition of the upper formation interface position, the constraint condition of the lower formation interface position, the prior value constraint condition of the upper formation conductivity, the prior value constraint condition of the target layer conductivity, the prior value constraint condition of the lower formation conductivity, and the prior value constraint conditions of the upper and lower formation interface positions and formation thickness.

[0168] Optionally, the constraint conditions include the mutual constraint condition of the upper formation conductivity; the inversion correction module 1102 is specifically configured to calculate the mutual constraint condition of the upper formation conductivity through a first constraint formula.

[0169] The first constraint formula is as follows:

[0170]

[0171] Among them, Npoint represents the number of inversion points taken, is the mutual constraint coefficient between the resistivity of the upper surrounding rock formation for the inversion of the k-th point and the l-th point, is the conductivity of the upper surrounding rock formation corresponding to the k-th point, is the conductivity of the upper surrounding rock formation corresponding to the l-th point.

[0172] The constraint conditions include the mutual constraint conditions of the conductivity of the target layer; the inversion correction module 1102 is specifically used to calculate the mutual constraint conditions of the conductivity of the target layer through the second constraint formula.

[0173] The second constraint formula is as follows:

[0174]

[0175] Among them, is the mutual constraint coefficient between the resistivity of the target layer for the inversion of the k-th point and the l-th point, is the resistivity of the target layer corresponding to the k-th point, is the resistivity of the target layer corresponding to the l-th point.

[0176] The constraint conditions include the mutual constraint conditions of the conductivity of the lower surrounding rock; the inversion correction module 1102 is specifically used to calculate the mutual constraint conditions of the conductivity of the lower surrounding rock through the third constraint formula.

[0177] The third constraint formula is as follows:

[0178]

[0179] Among them, is the mutual constraint coefficient between the resistivity of the lower surrounding rock formation for the inversion of the k-th point and the l-th point, is the resistivity of the lower surrounding rock formation corresponding to the k-th point, is the resistivity of the lower surrounding rock formation corresponding to the l-th point.

[0180] The constraint conditions include the constraint conditions of the position of the upper surrounding rock formation interface; the inversion correction module 1102 is specifically used to calculate the constraint conditions of the position of the upper surrounding rock formation interface through the fourth constraint formula.

[0181] The fourth constraint formula is as follows:

[0182]

[0183] Among them, is the mutual constraint coefficient between the positions of the upper surrounding rock interfaces for the inversion of the k-th point and the l-th point, is the position of the upper surrounding rock interface corresponding to the k-th point, is the position of the upper surrounding rock interface corresponding to the l-th point.

[0184] The constraint conditions include the constraint conditions for the position of the lower surrounding rock formation interface; the inversion correction module 1102 is specifically used to calculate the constraint conditions for the position of the lower surrounding rock formation interface through the fifth constraint formula.

[0185] The fifth constraint formula is:

[0186]

[0187] where, is the mutual constraint coefficient between the inversion results of the lower surrounding rock interface positions at the k-th and l-th points, is the position of the lower surrounding rock interface corresponding to the k-th point, is the position of the lower surrounding rock interface corresponding to the l-th point.

[0188] The constraint conditions include the prior value constraint conditions for the conductivity of the upper surrounding rock formation; the inversion correction module 1102 is specifically used to calculate the prior value constraint conditions for the conductivity of the upper surrounding rock formation through the sixth constraint formula.

[0189] The sixth constraint formula is:

[0190]

[0191] where, is the constraint coefficient between the resistivity of the upper surrounding rock formation obtained by inversion at the l-th point and the prior value, is the resistivity of the upper surrounding rock formation corresponding to the l-th point, is the reference value of the resistivity of the upper surrounding rock formation corresponding to the l-th point, and this reference value can be an existing prior value.

[0192] The constraint conditions include the prior value constraint conditions for the conductivity of the target layer; the inversion correction module 1102 is specifically used to calculate the prior value constraint conditions for the conductivity of the target layer through the seventh constraint formula.

[0193] The seventh constraint formula is:

[0194]

[0195] where, is the constraint coefficient between the resistivity of the target layer obtained by inversion at the l-th point and the prior value, is the resistivity of the target layer corresponding to the l-th point, is the reference value of the resistivity of the target layer corresponding to the l-th point, and this reference value can be an existing prior value.

[0196] The constraint conditions include the prior value constraint condition of the lower surrounding rock formation conductivity; the inversion correction module 1102 is specifically configured to calculate the prior value constraint condition of the lower surrounding rock formation conductivity through the eighth constraint formula.

[0197] The eighth constraint formula is:

[0198]

[0199] Wherein, is the constraint coefficient between the resistivity of the lower surrounding rock formation and the prior value for the inversion at the l-th point, is the resistivity of the lower surrounding rock formation corresponding to the l-th point, is the reference value of the resistivity of the lower surrounding rock formation corresponding to the l-th point, and this reference value can be an existing prior value.

[0200] The constraint conditions include the prior value constraint conditions of the upper and lower surrounding rock formation interface positions and formation thickness; the inversion correction module 1102 is specifically configured to calculate the prior value constraint conditions of the upper and lower surrounding rock formation interface positions and formation thickness through the ninth constraint formula.

[0201] The ninth constraint formula is:

[0202]

[0203] Wherein, is the constraint coefficient between the reference value of the lower surrounding rock formation interface and formation thickness for the inversion at the l-th point, is the position of the upper surrounding rock interface corresponding to the l-th point, is the position of the lower surrounding rock interface corresponding to the l-th point, H is the reference value of the formation thickness corresponding to the l-th point, and the reference value is a prior value that has been collected.

[0204] Optionally, the device further includes a collection module 1103, and the collection module is configured to collect the reference value of the upper surrounding rock formation conductivity, the reference value of the target layer conductivity, the reference value of the lower surrounding rock formation conductivity, and the reference value of the formation thickness from the side wellbore of the current logging wellbore.

[0205] Optionally, the environmental correction module 1101 is specifically configured to perform borehole correction on the logging data and determine the corrected logging data as preliminary logging data; or, perform skin effect correction and borehole correction on the logging data to obtain preliminary corrected data; perform a frequency conversion operation on the preliminary corrected data and determine the converted preliminary corrected data as preliminary logging data.

[0206] Optionally, the environmental correction module 1101 is specifically configured to perform a frequency conversion operation on the preliminary corrected data through the following preset conversion formula:

[0207]

[0208] where: k = iωμσ, σ is the preliminary calibration data, L R is the main coil spacing, L B is the shielding coil spacing, σ aR is the conductivity corresponding to the calibrated sub-array, μ is the magnetic permeability, i is the imaginary unit; among them, each sub-array corresponds to a σ aR , and the values of all sub-arrays are determined as the preliminary logging data.

[0209] The logging interpretation device based on array induction provided by the embodiments of the present application can execute the technical solutions of the logging interpretation method based on array induction in any of the above embodiments. Its implementation principle and beneficial effects are similar to those of the logging interpretation method based on array induction. For details, refer to the implementation principle and beneficial effects of the logging interpretation method based on array induction, which will not be elaborated here.

[0210] Figure 12 This is a schematic structural diagram of an electronic device provided by the present application. As Figure 12 shown, the electronic device 1200 may include: at least one processor 1201 and a memory 1202.

[0211] The memory 1202 is used to store programs. Specifically, the program may include program code, and the program code includes computer operation instructions.

[0212] The memory 1202 may include a high-speed RAM memory, and may also include non-volatile memory, such as at least one disk memory.

[0213] The processor 1201 is configured to execute the computer execution instructions stored in the memory 1202 to implement the interface display method described in the foregoing method embodiments. Among them, the processor 1201 may be a central processing unit (CPU), or an application specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of the present application. Specifically, when implementing the logging interpretation based on array induction described in the foregoing method embodiments, the electronic device may be an electronic device with processing functions such as a terminal or a server. When implementing the logging interpretation method based on array induction described in the foregoing method embodiments, the electronic device may be an electronic control unit on a vehicle.

[0214] Optionally, the electronic device 1200 may further include a communication interface 1203. In a specific implementation, if the communication interface 1203, the memory 1202, and the processor 1201 are implemented independently, the communication interface 1203, the memory 1202, and the processor 1201 may be interconnected through a bus and communicate with each other. The bus may be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, an Extended Industry Standard Architecture (EISA) bus, or the like. The bus may be divided into an address bus, a data bus, a control bus, etc., but it does not mean that there is only one bus or one type of bus.

[0215] Optionally, in a specific implementation, if the communication interface 1203, the memory 1202, and the processor 1201 are integrated on a single chip, the communication interface 1203, the memory 1202, and the processor 1201 may communicate through an internal interface.

[0216] The present application also provides a computer-readable storage medium, which may include: various media capable of storing program codes such as a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disc. Specifically, the computer-readable storage medium stores program instructions for the methods in the foregoing embodiments.

[0217] The present application also provides a program product, which includes execution instructions stored in a readable storage medium. At least one processor of the electronic device may read the execution instructions from the readable storage medium, and the at least one processor executes the execution instructions to enable the electronic device to implement the interface display methods provided by the various embodiments described above.

[0218] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present application and are not intended to limit them. Although the present application has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions described in the foregoing embodiments, or perform equivalent replacements for some or all of the technical features. However, these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the various embodiments of the present application.

Claims

1. A logging interpretation method based on array induction, characterized in that, Including: Collecting logging data of a horizontal well to be measured by array induction, and performing environmental correction on the logging data to obtain preliminary logging data. There is an associated relationship between the preliminary logging data and the measurement frequency. The logging data of the horizontal well to be measured includes initial conductivity; Performing inversion correction on the preliminary logging data to obtain the logging result of the horizontal well to be measured after inversion correction; The performing inversion correction on the preliminary logging data to obtain the logging result after inversion correction includes: Determining the constraint conditions of the logging data; Performing inversion correction on the preliminary logging data through a preset inversion formula and constraint conditions to obtain the logging result after inversion correction; The preset inversion formula is: Among them, is the error minimization objective function, is the preliminary logging data, is the sub-array response vector calculated from the initial value of the preliminary logging data , where the initial value , Constrants(1) represents the first constraint condition, Constrants(2) represents the second constraint condition, Constrants(N) represents the Nth constraint condition, and N is a natural number; The constraint conditions include the mutual constraint conditions of the conductivity of the upper surrounding rock. The determining the constraint conditions of the logging data includes: Calculating the mutual constraint conditions of the conductivity of the upper surrounding rock through a first constraint formula; The first constraint formula is: ， Among them, represents the number of inversion points taken, is the mutual constraint coefficient between the resistivity of the upper surrounding rock formation of the m-th point and the -th point inversion, is the conductivity of the upper surrounding rock formation corresponding to the m-th point, is the conductivity of the upper surrounding rock formation corresponding to the -th point.

2. The method according to claim 1, characterized in that, The constraint conditions further include: the mutual constraint conditions of the conductivity of the target layer, the mutual constraint conditions of the conductivity of the lower surrounding rock, the constraint conditions of the position of the upper surrounding rock formation interface, the constraint conditions of the position of the lower surrounding rock formation interface, the prior value constraint conditions of the conductivity of the upper surrounding rock formation, the prior value constraint conditions of the conductivity of the target layer, the prior value constraint conditions of the conductivity of the lower surrounding rock formation, and the prior value constraint conditions of the position and formation thickness of the upper and lower surrounding rock formations.

3. The method according to claim 2, wherein The constraint conditions further include the mutual constraint conditions of the conductivity of the target layer. The determining the constraint conditions of the logging data includes: Calculating the mutual constraint conditions of the conductivity of the target layer through a second constraint formula; The second constraint formula is: ， Among them, is the mutual constraint coefficient between the resistivity of the target layer for the inversion of the m-th point and the point, is the resistivity of the target layer corresponding to the m-th point, is the resistivity of the target layer corresponding to the point; The constraint conditions include the mutual constraint conditions of the conductivity of the lower surrounding rock. The determining the constraint conditions of the logging data includes: Calculating the mutual constraint conditions of the conductivity of the lower surrounding rock through a third constraint formula; The third constraint formula is: ， Among them, is the mutual constraint coefficient between the resistivity of the lower surrounding rock formation at the m-th point and the resistivity of the lower surrounding rock formation after inversion at the is the resistivity of the lower surrounding rock formation corresponding to the m-th point, is the resistivity of the lower surrounding rock formation corresponding to the The constraint conditions include the constraint conditions of the position of the upper surrounding rock formation interface. The determining the constraint conditions of the logging data includes: Calculating the constraint conditions of the position of the upper surrounding rock formation interface through a fourth constraint formula; The fourth constraint formula is: ， Among them, is the mutual constraint coefficient between the position of the upper surrounding rock interface of the inversion of the m-th point and the point, is the position of the upper surrounding rock interface corresponding to the m-th point, is the position of the upper surrounding rock interface corresponding to the point; The constraint conditions include the constraint conditions of the position of the lower surrounding rock formation interface. The determining the constraint conditions of the logging data includes: Calculating the constraint conditions of the position of the lower surrounding rock formation interface through a fifth constraint formula; The fifth constraint formula is: ， Among them, is the mutual constraint coefficient between the position of the lower surrounding rock interface of the m-th point and the position of the lower surrounding rock interface after inversion of the -th point, is the position of the lower surrounding rock interface corresponding to the m-th point, is the position of the lower surrounding rock interface corresponding to the -th point; The constraint conditions include the prior value constraint conditions of the conductivity of the upper surrounding rock formation. The determining the constraint conditions of the logging data includes: Calculating the prior value constraint conditions of the conductivity of the upper surrounding rock formation through a sixth constraint formula; The sixth constraint formula is: ， Among them, is the constraint coefficient between the resistivity of the upper surrounding rock formation at the point and the prior value, is the resistivity of the upper surrounding rock formation corresponding to the point, is the reference value of the resistivity of the upper surrounding rock formation corresponding to the point, and this reference value is the existing prior value; The constraint conditions include the prior value constraint conditions of the conductivity of the target layer. The determining the constraint conditions of the logging data includes: Calculating the prior value constraint conditions of the conductivity of the target layer through a seventh constraint formula; The seventh constraint formula is: ， Among them, is the constraint coefficient between the resistivity of the target layer for point inversion and the prior value, and is the resistivity of the target layer corresponding to the th point, and is the reference value of the resistivity of the target layer corresponding to the th point, and this reference value is the existing prior value. The constraint conditions include the prior value constraint conditions of the conductivity of the lower surrounding rock formation. The determining the constraint conditions of the logging data includes: Calculating the prior value constraint conditions of the conductivity of the lower surrounding rock formation through an eighth constraint formula; The eighth constraint formula is: ， Among them, is the constraint coefficient between the resistivity of the lower surrounding rock formation at the point and the prior value, is the resistivity of the lower surrounding rock formation corresponding to the point, is the reference value of the resistivity of the lower surrounding rock formation corresponding to the point, and this reference value is the existing prior value; The constraint conditions include the upper and lower surrounding rock formation interface positions and the prior value constraint conditions of the formation thickness; determining the constraint conditions of the logging data includes: Using the ninth constraint formula to calculate the upper and lower surrounding rock formation interface positions and the prior value constraint conditions of the formation thickness; The ninth constraint formula is: ， Among them, is the constraint coefficient between the lower surrounding rock formation interface and the reference value of the formation thickness for the point inversion, is the position of the upper surrounding rock interface corresponding to the point, is the position of the lower surrounding rock interface corresponding to the point, H is the reference value of the formation thickness corresponding to the point, and the reference value is the prior value that has been collected.

4. The method according to claim 3, wherein Before calculating each of the constraint conditions in the constraint conditions respectively through multiple constraint condition formulas, the method further includes: Collecting the upper surrounding rock formation conductivity reference value, the target layer conductivity reference value, the lower surrounding rock formation conductivity reference value, and the formation thickness reference value from the offset wellbore of the current logging wellbore.

5. The method according to any one of claims 1 to 4, characterized in that Performing environmental correction on the logging data to obtain preliminary logging data, including: Performing borehole correction on the logging data and determining the corrected logging data as the preliminary logging data; Or, Performing skin effect correction and borehole correction on the logging data to obtain preliminary corrected data; performing a frequency conversion operation on the preliminary corrected data and determining the converted preliminary corrected data as the preliminary logging data.

6. The method according to claim 5, wherein Performing the frequency conversion operation on the preliminary corrected data includes: Performing a frequency conversion operation on the preliminary corrected data through the following preset conversion formula: ， Wherein: , is the preliminary correction data, is the main coil pitch, is the shielding coil pitch, is the conductivity corresponding to the corrected sub-array, is the magnetic permeability, and i is the imaginary unit; Among them, each sub-array corresponds to a , and the values of all sub-arrays are determined as the preliminary logging data.

7. A logging interpretation device based on array induction, characterized in that, Including: An environmental correction module, configured to collect logging data of a horizontal well to be measured through array induction and perform environmental correction on the logging data to obtain preliminary logging data, where the preliminary logging data has an associated relationship with the measurement frequency, and the logging data of the horizontal well to be measured includes an initial conductivity; An inversion correction module, configured to perform inversion correction on the preliminary logging data to obtain the logging result of the horizontal well to be measured after inversion correction; The inversion correction module is specifically configured to determine the constraint conditions of the logging data; performing inversion correction on the preliminary logging data through a preset inversion formula and the constraint conditions to obtain the logging result after inversion correction; the preset inversion formula is: Among them, is the error minimization objective function, is the preliminary logging data, is the sub-array response vector calculated from the initial value of the preliminary logging data , where the initial value , Constrants(1) represents the first constraint condition, Constrants(2) represents the second constraint condition, Constrants(N) represents the Nth constraint condition, and N is a natural number; The constraint conditions include the mutual constraint conditions of the upper surrounding rock conductivity; the inversion correction module is specifically configured to calculate the mutual constraint conditions of the upper surrounding rock conductivity through the first constraint formula; the first constraint formula is: ， Among them, represents the number of inversion points taken, is the mutual constraint coefficient between the resistivity of the upper surrounding rock formation of the m-th point and the -th point inversion, is the conductivity of the upper surrounding rock formation corresponding to the m-th point, is the conductivity of the upper surrounding rock formation corresponding to the -th point.

8. An electronic device, comprising: A processor and a memory communicatively connected to the processor; The memory stores computer-executable instructions; The processor executes the computer-executable instructions stored in the memory to implement the method according to any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that, Computer-executable instructions are stored in the computer-readable storage medium, and when the computer-executable instructions are executed by the processor, they are used to implement the method according to any one of claims 1 to 6.

Images