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

Abstract

The present application provides a logging interpretation method, device, equipment and storage medium based on array induction. Logging data is collected through array induction, and environmental correction is performed on the logging data to obtain preliminary logging data, and inclination correction is performed on the preliminary logging data to obtain the logging result after inclination correction. In this method, through environmental correction and inclination correction, the influence caused by the inclination angle is avoided in the logging result, and the problem that the inclination angle affects the accuracy of the measurement result during logging operations on non-horizontal ground in the prior art is solved, and the logging result is closer to the true value, improving the accuracy of logging.

Description

Technical Field

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

[0002] In fields such as geological survey and wellbore measurement, array induction is usually used for logging operations. Among them, during measurement, environmental correction is generally performed on the collected logging data, and then focusing and resolution matching operations are performed using software to obtain the logging result.

[0003] Currently, in the existing logging process, the array induction method is often used to solve two-dimensional geological problems. During the operation, it is often aimed at the measurement of vertical wells on a horizontal ground, that is, the instrument is often operated with its axis perpendicular to the formation interface.

[0004] However, in practical applications, the formation distribution is not always horizontal, that is, the instrument is not perpendicular to the formation interface. Therefore, when the measurement instrument is not perpendicular to the ground, there will be a large error between the current logging method's measurement result and the true result, resulting in a poor accuracy problem of the logging result when the formation dip angle is large. Summary of the Invention

[0005] The embodiments of the present application provide a logging interpretation method, device, equipment, and storage medium based on array induction to solve the problem of poor accuracy of the logging results of the prior art when the formation dip angle is large.

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

[0007] Collect logging data through array induction, and perform environmental correction on the logging data to obtain preliminary logging data;

[0008] Perform tilt correction on the preliminary logging data to obtain the logging result after tilt correction.

[0009] In a possible design of the first aspect, before performing tilt correction on the preliminary logging data to obtain the logging result after tilt correction, the method further includes:

[0010] Obtain formation dip angle data;

[0011] Correspondingly, performing tilt correction on the preliminary logging data to obtain the logging result after tilt correction includes:

[0012] Perform tilt correction on the preliminary logging data according to the formation dip angle data to obtain the corrected logging result.

[0013] In this possible design, the method further includes:

[0014] Performing a focusing process on the corrected logging result to obtain an optimized logging result, where the focusing process includes software focusing processing and longitudinal resolution matching.

[0015] In this possible design, the step of performing dip correction on the preliminary logging data according to the formation dip data to obtain a corrected logging result includes:

[0016] Performing a conductivity conversion operation on the preliminary logging data to obtain preprocessed data, where the preprocessed data is associated with frequency;

[0017] Performing an inversion process on the preprocessed data to obtain the corrected logging result.

[0018] Optionally, the step of performing a conductivity conversion operation on the preliminary logging data to obtain preprocessed data includes:

[0019] Converting the first conductivity in the preliminary logging data through a preset first formula to obtain a second conductivity having an associated relationship with frequency

[0020] Performing an inversion process on the second conductivity using a constraint condition and a preset second formula to obtain a corrected conductivity, where the preprocessed data includes the corrected conductivity;

[0021] Wherein, the preset first formula is:

[0022]

[0023] Where k = iωμσ, L R is the main coil spacing, L B is the shield coil spacing, σ aR is the apparent conductivity corresponding to the corrected subarray, μ is the magnetic permeability, ω = 2πF0, i is the imaginary unit, and each subarray corresponds to a σ aR , and here the values of all subarrays are put into a vector

[0024] The preset second formula is an inversion penalty function, and the inversion penalty function is:

[0025]

[0026] Where is the initial value of the selected formation conductivity The calculated sub-array response vector, the initial formation conductivity is the conductivity after environmental correction, Constrants(i) is the i-th constraint condition, and i is a natural number.

[0027] Optionally, the constraint conditions at least include the interlayer conductivity constraint condition and the inter-sub-array conductivity constraint condition;

[0028] Correspondingly, the inversion process of the second conductivity using the constraint conditions and the preset second formula to obtain the corrected conductivity includes:

[0029] Inverting the second conductivity through the preset third formula, the interlayer conductivity constraint condition, and the inter-sub-array conductivity constraint condition to obtain the corrected conductivity;

[0030] The preset third formula is:

[0031]

[0032] Wherein, Constrants(1) is the interlayer conductivity constraint condition, and Constrants(2) is the inter-sub-array conductivity constraint condition;

[0033] The interlayer conductivity constraint condition includes:

[0034]

[0035] Wherein, NArry is the number of sub-arrays, Nlayer is the number of formations, σ i，j is the conductivity of the j-th sub-array in the i-th layer, is the constraint value of the conductivity of the j-th sub-array in the i-th layer and the weighted average of the conductivities of the surrounding formations of the i-th layer, and α i，j is the coefficient for adding constraints to the conductivity of the j-th sub-array in the i-th layer;

[0036] The inter-sub-array conductivity constraint condition includes:

[0037]

[0038] Wherein, NArry is the number of sub-arrays, σ k is the conductivity of the k-th sub-array, σ l is the conductivity of the l-th sub-array, and β k，l is the coefficient for adding constraints to the conductivity between the k-th sub-array and the k-th sub-array.

[0039] In a second aspect, an array induction-based logging interpretation device provided by an embodiment of the present application includes: an environmental correction module and a tilt correction module;

[0040] The environmental correction module is used to collect logging data through array induction and perform environmental correction on the logging data to obtain preliminary logging data;

[0041] The tilt correction module is used to perform tilt correction on the preliminary logging data to obtain the logging result after tilt correction.

[0042] In a possible design of the second aspect, before the tilt correction module is used to perform tilt correction on the preliminary logging data to obtain the logging result after tilt correction, it is further used for:

[0043] Obtain formation dip data;

[0044] Correspondingly, the tilt correction module is used to perform tilt correction on the preliminary logging data to obtain the logging result after tilt correction, specifically:

[0045] The tilt correction module is used to perform tilt correction on the preliminary logging data according to the formation dip data to obtain the corrected logging result.

[0046] In this possible design, the tilt correction module is further used to perform a focusing processing operation on the corrected logging result to obtain an optimized logging result, and the focusing processing operation includes software focusing processing and longitudinal resolution matching.

[0047] In this possible design, the tilt correction module is used to perform tilt correction on the preliminary logging data according to the formation dip data to obtain the corrected logging result, specifically:

[0048] The tilt correction module is specifically used for:

[0049] Perform a conductivity conversion operation on the preliminary logging data to obtain preprocessed data, and the preprocessed data is associated with frequency;

[0050] Perform an inversion process on the preprocessed data to obtain the corrected logging result.

[0051] Optionally, the tilt correction module is used to perform a conductivity conversion operation on the preliminary logging data to obtain preprocessed data, specifically:

[0052] The tilt correction module is specifically used for:

[0053] Convert the first conductivity in the preliminary logging data through a preset first formula to obtain a second conductivity that has an associated relationship with frequency

[0054] The second conductivity is inversely processed using the constraint conditions and the preset second formula to obtain the corrected conductivity, and the preprocessed data includes the corrected conductivity;

[0055] Among them, the preset first formula is:

[0056]

[0057] Among them, k = iωμσ, L R is the main coil spacing, L B is the shielding coil spacing, σ aR is the apparent conductivity corresponding to the corrected subarray, μ is the magnetic permeability, ω = 2πF0, i is the imaginary unit, and each subarray corresponds to a σ aR , and here the values of all subarrays are put into a vector

[0058] The preset second formula is an inversion penalty function, and the inversion penalty function is:

[0059]

[0060] Among them, is the initial value of the formation conductivity selected The subarray response vector calculated, the initial value of the formation conductivity is the conductivity after environmental correction, Constrants(i) is the i-th constraint condition, and i is a natural number.

[0061] Optionally, the constraint conditions at least include the interlayer conductivity constraint condition and the inter-subarray conductivity constraint condition;

[0062] Correspondingly, the tilt correction module is specifically used to inversely process the second conductivity using the constraint conditions and the preset second formula to obtain the corrected conductivity, specifically:

[0063] The tilt correction module is specifically used for:

[0064] The second conductivity is inversely processed through the preset third formula, the interlayer conductivity constraint condition and the inter-subarray conductivity constraint condition to obtain the corrected conductivity;

[0065] The preset third formula is:

[0066]

[0067] Among them, Constrants(1) is the interlayer conductivity constraint condition, and Constrants(2) is the inter-subarray conductivity constraint condition;

[0068] The interlayer conductivity constraint condition includes:

[0069]

[0070] Among them, NArry is the number of sub-arrays, Nlayer is the number of formation layers, and σ i，j is the conductivity of the j-th sub-array in the i-th layer, is the constraint value of the conductivity of the j-th sub-array in the i-th layer and the weighted average of the conductivities of the formation layers around the i-th layer, and α i，j is the coefficient for adding constraints to the conductivity of the j-th sub-array in the i-th layer;

[0071] The conductivity constraint conditions between sub-arrays include:

[0072]

[0073] Among them, NArry is the number of sub-arrays, and σ k is the conductivity of the k-th sub-array, and σ l is the conductivity of the l-th sub-array, and β k，l is the coefficient for adding constraints to the conductivity between the k-th sub-array and the k-th sub-array.

[0074] In a third aspect, an embodiment of the present application provides an electronic device, including a processor, a memory, and a computer program stored on the memory and executable on the processor. When the processor executes the program, it implements the method provided in the first aspect and each possible design above.

[0075] In a fourth aspect, an embodiment of the present application may provide a computer-readable storage medium, in which computer-executable instructions are stored. When the computer-executable instructions are executed by a processor, they are used to implement the method provided in the first aspect and each possible design above.

[0076] In a fifth aspect, an embodiment of the present application provides a computer program product, including a computer program. When the computer program is executed by a processor, it is used to implement the method provided in the first aspect and each possible design above.

[0077] For the logging interpretation method, device, equipment, and storage medium provided in the embodiments of the present application, in this method, logging data is collected through array induction, and the logging data is subjected to environmental correction to obtain preliminary logging data, and the preliminary logging data is subjected to tilt correction to obtain the logging result after tilt correction. In this method, through environmental correction and tilt correction, the influence caused by the tilt angle is avoided in the logging result, and the problem that the tilt angle affects the accuracy of the measurement result during logging operations on non-horizontal ground in the prior art is solved, making the logging result closer to the true value and improving the accuracy of logging.

[0078] The above description is only an overview of the technical solution of the present invention. In order to better understand the technical means of the present invention, it can be implemented according to the content of the specification. And in order to make the above and other objects, features and advantages of the present invention more obvious and understandable, the specific embodiments of the present invention are hereinafter specifically exemplified. BRIEF DESCRIPTION OF THE DRAWINGS

[0079] The drawings herein are incorporated into the specification and form a part of the specification, showing embodiments consistent with the present disclosure and used together with the specification to explain the principles of the present disclosure.

[0080] Figure 1A Schematic diagram of the principle of measuring formation conductivity by a dual-coil array in the prior art provided by an embodiment of the present application;

[0081] Figure 1B Schematic diagram of the principle of measuring formation conductivity by a triple-coil array in the prior art provided by an embodiment of the present application;

[0082] Figure 1C Schematic diagram one of the principle of the measuring instrument provided by an embodiment of the present application;

[0083] Figure 1D Schematic diagram of the principle of the measuring instrument provided by an embodiment of the present application Figure 2 ;

[0084] Figure 2 Schematic flowchart of Embodiment 1 of the logging interpretation method based on array induction provided by an embodiment of the present application;

[0085] Figure 3A Schematic diagram of vector decomposition of an equivalent magnetic dipole in a geodetic coordinate system provided by an embodiment of the present application;

[0086] Figure 3B Schematic diagram of the horizontal conductivity and vertical conductivity of a horizontal formation provided by an embodiment of the present application;

[0087] Figure 3C Schematic diagram of the horizontal conductivity and vertical conductivity of an inclined formation provided by an embodiment of the present application;

[0088] Figure 3D Schematic diagram of the logging data response of an inclined formation corresponding to a horizontal formation provided by an embodiment of the present application;

[0089] Figure 4 Schematic flowchart of Embodiment 2 of the logging interpretation method based on array induction provided by an embodiment of the present application;

[0090] Figure 5 Schematic diagram of the structure of the logging interpretation device based on array induction provided by an embodiment of the present application;

[0091] Figure 6 Schematic structural diagram of the electronic device provided by an embodiment of the present application.

[0092] Through the above-mentioned drawings, specific embodiments of the present disclosure have been shown, and more detailed descriptions will be given hereinafter. These drawings and written descriptions 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. Detailed description of specific embodiments

[0093] Here, the exemplary embodiments will be described in detail, and the 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.

[0094] Before introducing the embodiments of the present application, the background technology of the present application will be explained first.

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

[0096] Currently, in the existing well logging process, the array induction-based method is often used to solve two-dimensional geological problems. During the operation, it is often aimed at the vertical well measurement on the horizontal ground, that is, the well logging instrument is often operated with its axis perpendicular to the formation interface.

[0097] Specifically, in the existing well logging process, it is usually carried out in the manner as Figure 1A shown for well logging. Figure 1A Schematic diagram of the principle of measuring formation conductivity by a dual-coil system in the prior art provided by an embodiment of the present application. As Figure 1A shown, the transmitting coil and the receiving coil are on the instrument rod. When the transmitting coil emits electromagnetic waves, a primary field is generated in the surrounding medium. This primary field generates an eddy current ring in the formation, and this eddy current ring will generate a secondary field. The receiving coil will receive these primary fields and secondary fields. Among them, z represents the longitudinal distance of the formation, p represents the transverse distance of the formation, r R represents the distance from the receiving coil to the eddy current ring, r T represents the distance from the transmitting coil to the eddy current ring, and L represents the distance between the receiving coil and the transmitting coil. Since the secondary field is related to the formation conductivity, the received signal can be used to calculate the formation conductivity. However, since the primary field is much larger than the secondary field, the well logging instrument must overcome the influence of the primary field signal.

[0098] To overcome the influence of the primary field signal, Figure 1B The following is a schematic diagram of the principle of measuring formation conductivity using a three-coil array in the prior art provided by the embodiments of the present application. As Figure 1B shown, it includes: a main receiving coil on the upper side, a shielding receiving coil in the middle, and a transmitting coil on the lower side. Among them, P(p, z) represents the coordinates of the observation point, r R represents the distance from the main receiving coil to the observation point, r B represents the distance from the shielding receiving coil to the observation point, r T represents the distance from the transmitting coil to the observation point, L B represents the distance from the shielding receiving coil to the transmitting coil, L T represents the distance from the main receiving coil to the transmitting coil. The winding directions of the main receiving coil and the shielding receiving coil are opposite, and they are connected in series. When the number of turns of the two receiving coils and their distances to the transmitting coil are appropriately adjusted, the primary field signal can be eliminated. And the logging tool used in the current array induction logging process is composed of multiple three-coil systems as Figure 1B shown.

[0099] In an ideal situation, Figure 1C The following is Schematic Diagram 1 of the principle of the measuring instrument provided by the embodiments of the present application. As Figure 1C shown, the eddy current ring is parallel to the formation interface. Among them, r R represents the distance from the receiving coil to the eddy current ring, L T represents the distance from the transmitting coil to the eddy current ring. When the formation is a two-dimensional formation, or the axis of the instrument is perpendicular to the formation interface, the eddy current ring is parallel to the formation interface. At this time, the eddy current does not intersect the formation interface, and the electromagnetic wave emitted by the transmitting coil will not undergo reflection and refraction phenomena.

[0100] However, in practical applications, the formation distribution is not always horizontal, that is, the instrument is not perpendicular to the formation interface. Figure 1D The following is the schematic diagram of the principle of the measuring instrument provided by the embodiments of the present application Figure 2 . As Figure 1D shown, the eddy current ring intersects the formation interface. Among them, r R represents the distance from the receiving coil to the eddy current ring, L T represents the distance from the transmitting coil to the eddy current ring. That is, when the axis of the measuring instrument is not perpendicular to the formation interface and the eddy current ring intersects the formation interface, the eddy current ring will undergo reflection and refraction phenomena on the formation interface. This causes the measurement method not to meet the measurement principle of the array induction instrument, so the measured result cannot obtain the accurate formation resistivity, and the greater the formation dip angle, the greater the difference between the measured formation resistivity and the actual result.

[0101] Therefore, when the current logging method is used and the measuring instrument is not perpendicular to the ground, there will be a large error between the measurement result and the true result.

[0102] In view of the above problems, the inventive concept of the present application is as follows: During the process of using array induction logging, the inventors found that in practical applications, due to the non-horizontal distribution of the formation, there is a certain relationship between the array induction logging data collected and the formation dip angle. The logging data can be corrected by using the formation dip angle-related data, and the obtained logging result can avoid the influence caused by the tilt angle, thus solving the problems existing in the prior art.

[0103] Next, the technical solution of the present application will be described in detail through specific embodiments.

[0104] It should be noted 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 2 FIG. is a schematic flowchart of the first embodiment of the logging interpretation method based on array induction provided by the embodiment of the present application. As Figure 2 shown, the logging interpretation method based on array induction may include the following steps:

[0106] Step 21: Collect logging data through array induction and perform environmental correction on the logging data to obtain preliminary logging data.

[0107] In this step, the logging data collected based on the array method cannot be directly applied in practical applications, and the data needs to be corrected based on the environment to obtain monitoring data that conforms to the actual situation. The environmental correction may include skin effect correction and borehole correction.

[0108] Specifically, during the measurement based on conductivity, there is often interference from the skin effect. The skin effect refers to the phenomenon that when there is alternating current or alternating electromagnetic field in a conductor, the current distribution inside the conductor is uneven, and the current is concentrated in the "skin" part of the conductor, that is, the current is concentrated in a thin layer on the outer surface of the conductor. The closer to the surface of the conductor, the greater the current density, and the actual current inside the conductor is smaller. As a result, the resistance of the conductor increases, and its power loss also increases. This phenomenon is called the skin effect.

[0109] Therefore, the skin effect affects the accuracy of the measurement result and corresponding correction is required. Of course, in the actual application process, a suitable correction method can be selected based on the actual situation, and based on the logging data after environmental correction, preliminary logging data can be obtained.

[0110] Step 22: Perform tilt correction on the preliminary logging data to obtain the logging result after tilt correction.

[0111] In order to explain this scheme more clearly, the transmitting coil and receiving coil in the measuring instrument are first explained. Figure 3A Schematic diagram of vector decomposition of an equivalent magnetic dipole in the geodetic coordinate system provided in the embodiment of the present application. Figure 3A As shown, both the transmitting coil and the receiving coil can be equivalent to magnetic dipoles, and the transmitting coil can be expressed as the magnetic dipole moment M T ,M T is decomposed into components parallel to the bed interface and the component perpendicular to the bed interface The receiving coil can be expressed as the magnetic dipole moment M R ,M R is decomposed into components parallel to the bed interface and the component perpendicular to the bed interface

[0112] At this time, when the transmitting coil M T After the electromagnetic wave is transmitted into the ground, the receiving coil M R The electromagnetic wave signal is received and converted into an induced electromotive force, expressed as V = V ZZ +V ZX +V XX +V xz ,

[0113] Among them, V ZZ yes The component transmits electromagnetic waves in The induced electromotive force on the component, V ZX yes The component transmits electromagnetic waves in The induced electromotive force on the component, V XX yes Components transmit electromagnetic waves in The induced electromotive force on the component, V xz yes The component transmits electromagnetic waves in The induced electromotive force on the component.

[0114] Optionally, in combination with the actual use of the measuring instrument, Figure 3B Schematic diagram of horizontal conductivity and vertical conductivity of horizontal strata provided in the embodiment of the present application. Figure 3B As shown, σ h is the formation conductivity parallel to the formation interface, called horizontal conductivity, σ v The conductivity perpendicular to the formation interface is called vertical conductivity. Figure 3C Schematic diagram of horizontal conductivity and vertical conductivity of inclined strata provided in the embodiment of the present application. his the formation conductivity in the horizontal direction, σ v is the formation conductivity in the vertical direction, r R represents the distance from the receiving coil to the eddy current ring, L T represents the distance from the transmitting coil to the eddy current ring, σ h‘ is the horizontal conductivity relative to the instrument, σ v‘ is the vertical conductivity relative to the instrument. Obviously, when the axis of the instrument is not perpendicular to the formation interface, σ h‘ and σ h are not equal, and σ h‘ is a function of σ h and σ v , and σ v is related to the formation interface.

[0115] Optionally, Figure 3D is a schematic diagram of the logging data response corresponding to the inclined formation and the horizontal formation provided by the embodiment of the present application. As Figure 3D shown, in the case of multiple layers of media, the response of a sub-array when the formation is inclined at 75° (DGR = 75) is compared with the response of a sub-array when the horizontal formation is inclined at 0° (DGR = 0). Among them, the horizontal axis is the depth / foot, and the vertical axis is the imaginary part of the induced electromotive force / V. Among them, the dark gray line shows the logging data of the horizontal formation, which is relatively smooth; the light gray line is the logging data of the formation with an inclination angle of 75°, which is not smooth. It can be seen that since the array induction is developed for the horizontal formation of vertical wells, the response when the formation inclination angle is 75° and the logging results required by the instrument to solve geological problems have a large gap compared with the actual situation. As a result, the current logging method has a problem of low accuracy of the detection results when conducting geological surveys or logging for formations with inclination angles.

[0116] In this step, when the logging process corrects the logging data for the environment, the obtained preliminary logging data is affected by the detected inclination angle on the logging results. Therefore, the preliminary logging data can be corrected for inclination, so as to solve the influence of the inclination angle between the ground and the monitoring equipment on the detection results.

[0117] Specifically, the process of inclination correction can be achieved by obtaining the angle of the inclination angle of the preliminary logging data, and then based on the inclination angle and the influence of the inclination angle on the logging data at different angles for correction, so as to obtain the corrected logging data as the logging data after inclination correction. Specifically, in the process of correction, the correction method can be to perform inversion processing on the logging data based on the constraint conditions by using the inversion formula, so as to achieve the processing result under the constraint conditions, so that the preliminary logging data can be inverted and converged based on the constraint conditions to obtain the logging results that conform to the actual inclination angle.

[0118] The logging interpretation method based on array induction provided by the embodiment of the present application collects logging data through array induction, performs environmental correction on the logging data to obtain preliminary logging data, and performs tilt correction on the preliminary logging data to obtain the logging result after tilt correction. Therefore, when there is a tilt angle during the logging process, tilt correction can be performed to avoid the influence caused by the tilt angle on the logging result, solve the problem that the tilt angle affects the accuracy of the measurement result in the prior art when performing logging operations on non-horizontal ground, make the logging result closer to the true value, and improve the accuracy of logging.

[0119] Figure 4 It is a schematic flowchart of the second embodiment of the logging interpretation method based on array induction provided by the embodiment of the present application. As Figure 4 shown, the logging interpretation method may further include the following steps:

[0120] Step 41: Obtain formation dip data.

[0121] In practical applications, the ground is not always horizontal. Therefore, the logging equipment is not always perpendicular to the ground during the measurement process. Therefore, based on the influence of the dip angle described in the above embodiment, when there is a measurement dip angle during the measurement process, it will interfere with the accuracy of the result. Therefore, it is necessary to correct and process this interference.

[0122] Based on this, after obtaining the preliminary logging data through environmental correction, the preliminary logging data is not sensitive to the formation dip angle. Therefore, when performing tilt correction processing subsequently, it is also necessary to obtain formation dip data from other logging data to provide constraint conditions for tilt correction.

[0123] Optionally, the method and timing for obtaining the formation dip data are not specifically limited.

[0124] Step 42: Perform tilt correction on the preliminary logging data according to the formation dip data to obtain the corrected logging result.

[0125] In this step, perform a preset processing operation on the preliminary logging data to obtain preprocessed data, where the preprocessed data has an associated relationship with the frequency; then perform inversion processing on the preprocessed data to obtain the corrected logging result.

[0126] Optionally, performing a preset processing operation on the preliminary logging data to obtain preprocessed data may be: performing a conductivity conversion operation on the preliminary logging data to obtain preprocessed data. Among them, the preprocessed data is associated with the frequency.

[0127] Specifically, convert the first conductivity in the preliminary logging data through a preset first formula to obtain a second conductivity that has an associated relationship with the frequency The second conductivity is inversely processed using the constraint conditions and a preset second formula to obtain the corrected conductivity. The preprocessed data includes the corrected conductivity.

[0128] Based on the well logging data after environmental correction in the above steps, it is conductivity (or resistivity) well logging data independent of frequency. At this time, it is necessary to convert the conductivity independent of frequency into conductivity (or resistivity) related to frequency. For this purpose, a frequency can be first selected. This frequency can be arbitrarily selected from 0.1 Hz to 1,000,000 Hz and is set as F0. Then, the conductivity after environmental correction can be converted into a response related to frequency using the following first formula.

[0129] Among them, the preset first formula is:

[0130]

[0131] Among them, k = iωμσ, L R is the main coil spacing, L B is the shielding coil spacing, σ aR is the apparent conductivity corresponding to the corrected subarray, μ is the magnetic permeability, ω = 2πF0, i is the imaginary unit, and each subarray corresponds to a σ aR , and here the values of all subarrays are put into a vector

[0132] It should be noted that the selected frequency in the above scheme is only an example. In the specific operation process, multiple frequencies can also be selected, and the above operations are performed separately for the selected multiple frequencies.

[0133] Optionally, the preprocessed data is inversely processed to obtain the corrected well logging result.

[0134] Specifically, the second conductivity is inversely processed using the constraint conditions and a preset second formula to obtain the corrected conductivity.

[0135] The preset second formula is an inversion penalty function, and the inversion penalty function is:

[0136]

[0137] Among them, is the initial value of the selected formation conductivity The subarray response vector calculated, the initial value of the formation conductivity is the conductivity after environmental correction, Constrants(i) is the i-th constraint condition, and i is a natural number.

[0138] Furthermore, based on the actual operation process, there may be multiple constraint conditions, and the most important constraint conditions may include: the interlayer conductivity constraint condition and the inter-subarray conductivity constraint condition; therefore, based on these two constraint conditions, the process of performing inversion in this step can be: using a preset third formula, the interlayer conductivity constraint condition, and the inter-subarray conductivity constraint condition to perform inversion on the second conductivity to obtain the corrected conductivity.

[0139] Correspondingly, using the constraint condition and the preset second formula to perform inversion on the second conductivity to obtain the corrected conductivity includes:

[0140] Using a preset third formula, the interlayer conductivity constraint condition, and the inter-subarray conductivity constraint condition to perform inversion on the second conductivity to obtain the corrected conductivity;

[0141] The preset third formula is:

[0142]

[0143] where Constrants(1) is the interlayer conductivity constraint condition and Constrants(2) is the inter-subarray conductivity constraint condition;

[0144] The interlayer conductivity constraint condition includes:

[0145]

[0146] where NArry is the number of subarrays, Nlayer is the number of formations, σ i，j is the conductivity of the jth subarray in the ith layer, is the constraint value of the conductivity of the jth subarray in the ith layer and the weighted average of the conductivities of the surrounding formations of the ith layer, and α i，j is the coefficient for adding a constraint to the conductivity of the jth subarray in the ith layer;

[0147] The inter-subarray conductivity constraint condition includes:

[0148]

[0149] where NArry is the number of subarrays, σ k is the conductivity of the kth subarray, σ l is the conductivity of the lth subarray, and β k，l is the coefficient for adding a constraint to the conductivity between the kth subarray and the kth subarray.

[0150] In addition, when performing the above inversion, all subarrays share the same tilt angle and the same formation interface, thus ensuring the accuracy of tilt correction.

[0151] Step 43: Perform a focusing process on the corrected logging results to obtain optimized logging results.

[0152] Among them, the focusing process includes software focusing processing and longitudinal resolution matching.

[0153] After the above step 42, the corrected logging results can be output to the user. However, in order to further improve the accuracy and readability of the logging results, in this step, the corrected logging results in the previous step can also be focused through a preset software.

[0154] Specifically, the operation processes of its software focusing processing and longitudinal resolution matching can be processed in an existing manner and will not be elaborated here.

[0155] Figure 5 It is a schematic structural diagram of a logging interpretation device based on array induction provided by an embodiment of the present application. As Figure 5 shown, the logging interpretation device includes: an environmental correction module 51 and an inclination correction module 52.

[0156] The environmental correction module 51 is used to collect logging data through array induction and perform environmental correction on the logging data to obtain preliminary logging data;

[0157] The inclination correction module 52 is used to perform inclination correction on the preliminary logging data to obtain the logging results after inclination correction.

[0158] In a possible design of an embodiment of the present application, before the inclination correction module 52 is used to perform inclination correction on the preliminary logging data to obtain the logging results after inclination correction, it is also used for:

[0159] Obtain formation dip data;

[0160] Correspondingly, the inclination correction module 52 is used to perform inclination correction on the preliminary logging data to obtain the logging results after inclination correction, specifically:

[0161] The inclination correction module 52 is used to perform inclination correction on the preliminary logging data according to the formation dip data to obtain the corrected logging results.

[0162] In this possible design, the inclination correction module 52 is also used to perform a focusing process on the corrected logging results to obtain optimized logging results, and the focusing process includes software focusing processing and longitudinal resolution matching.

[0163] In this possible design, the inclination correction module 52 is used to perform inclination correction on the preliminary logging data according to the formation dip data to obtain the corrected logging results, specifically:

[0164] The tilt correction module 52 is specifically used for:

[0165] Perform conductivity conversion operation on the preliminary logging data to obtain preprocessed data, and the preprocessed data is associated with the frequency;

[0166] Perform inversion processing on the preprocessed data to obtain the corrected logging result.

[0167] Optionally, the tilt correction module 52 is used to perform conductivity conversion operation on the preliminary logging data to obtain preprocessed data. Specifically:

[0168] The tilt correction module 52 is specifically used for:

[0169] Convert the first conductivity in the preliminary logging data through a preset first formula to obtain a second conductivity that has an associated relationship with the frequency

[0170] Perform inversion processing on the second conductivity using the constraint conditions and a preset second formula to obtain the corrected conductivity, and the preprocessed data includes the corrected conductivity;

[0171] Among them, the preset first formula is:

[0172]

[0173] Among them, k = iωμσ, L R is the main coil spacing, L B is the shielding coil spacing, σ aR is the apparent conductivity corresponding to the corrected subarray, μ is the magnetic permeability, ω = 2πF0, i is the imaginary unit, and each subarray corresponds to a σ aR , and here the values of all subarrays are put into a vector

[0174] The preset second formula is an inversion penalty function, and the inversion penalty function is:

[0175]

[0176] Among them, is the initial value of the formation conductivity selected for calculating the subarray response vector. The initial value of the formation conductivity is the conductivity after environmental correction, Constrants(i) is the i-th constraint condition, and i is a natural number.

[0177] Optionally, the constraint conditions at least include the interlayer conductivity constraint condition and the inter-subarray conductivity constraint condition;

[0178] Correspondingly, the tilt correction module 52 is specifically configured to perform inversion processing on the second conductivity by using the constraint conditions and a preset second formula to obtain the corrected conductivity, specifically as follows:

[0179] The tilt correction module 52 is specifically configured to:

[0180] Perform inversion processing on the second conductivity through a preset third formula, the inter-layer conductivity constraint condition, and the inter-subarray conductivity constraint condition to obtain the corrected conductivity;

[0181] The preset third formula is:

[0182]

[0183] where Constrants(1) is the inter-layer conductivity constraint condition, and Constrants(2) is the inter-subarray conductivity constraint condition;

[0184] The inter-layer conductivity constraint condition includes:

[0185]

[0186] where NArry is the number of subarrays, Nlayer is the number of strata, σ i，j is the conductivity of the j-th subarray in the i-th layer, is the constraint value of the conductivity of the j-th subarray in the i-th layer and the weighted average of the conductivities of the surrounding strata of the i-th layer, and α i，j is the coefficient for adding constraints to the conductivity of the j-th subarray in the i-th layer;

[0187] The inter-subarray conductivity constraint condition includes:

[0188]

[0189] where NArry is the number of subarrays, σ k is the conductivity of the k-th subarray, σ l is the conductivity of the l-th subarray, and β k，l is the coefficient for adding constraints to the conductivity between the k-th subarray and the k-th subarray.

[0190] The logging interpretation device provided in this embodiment can be used to execute the solutions in the above method embodiments, and its implementation principle and technical effects are similar, which will not be elaborated here.

[0191] It should be noted that it should be understood that the division of each module of the above device is only a division of logical functions. In actual implementation, it can be fully or partially integrated into a physical entity, or physically separated. And these modules can all be implemented in the form of software called by processing elements; they can also all be implemented in the form of hardware; some modules can also be implemented in the form of software called by processing elements, and some modules can be implemented in the form of hardware. In addition, all or part of these modules can be integrated together or can be independently implemented. In the implementation process, each step of the above logging interpretation method or each of the above modules can be completed by the integrated logic circuit of the hardware in the processor element or the instructions in the form of software.

[0192] Figure 6 FIG. is a schematic structural diagram of an electronic device provided by an embodiment of the present application. As Figure 6 shown, the electronic device may include: a processor 61, a memory 62, and a display 63.

[0193] Among them, the processor 61 executes the computer program instructions stored in the memory 62, so that the processor 61 executes the solutions in the above embodiments. The processor 61 may be a general-purpose processor, including a central processing unit CPU, a network processor (NP), etc.; it may also be a digital signal processor DSP, an application specific integrated circuit ASIC, a field programmable gate array FPGA, or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components.

[0194] The display 63 may be a user interface, which may be used to display logging results, etc. in the above embodiments. The user interface may include graphics, text, icons, videos, and any combination thereof. The display 63 may also be used to provide virtual buttons and / or a virtual keyboard, also referred to as soft buttons and / or a soft keyboard. In some embodiments, the display 63 may be the front panel of the electronic device; in other embodiments, the display 63 may be a flexible display screen, disposed on a curved surface or a folding surface of the electronic device. Even, the display 63 may also be set as a display screen with an irregular non-rectangular shape, that is, a special-shaped screen. The display 63 may be prepared from materials such as a liquid crystal display (LCD), an organic light-emitting diode (OLED), etc.

[0195] Optionally, the above-mentioned various devices of the electronic device may be connected through a system bus.

[0196] The electronic device provided by the embodiment of the present application can be used to execute the solutions in the above embodiments, and its implementation principle and technical effects are similar, and will not be elaborated here.

[0197] An embodiment of this application also provides a chip for running instructions, and the chip is used to execute the solution in the above embodiment.

[0198] An embodiment of this application also provides a computer-readable storage medium, in which computer instructions are stored. When the computer instructions run on a computer, the computer is caused to execute the solution in the above embodiment.

[0199] An embodiment of this application also provides a computer program product, which includes a computer program stored in a computer-readable storage medium. At least one processor can read the computer program from the computer-readable storage medium, and when at least one processor executes the computer program, the solution in the above embodiment can be implemented.

[0200] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of this application, rather than to limit them; although this 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 recorded in the foregoing embodiments, or perform equivalent replacements on some or all of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the scope of the technical solutions of the embodiments of this application.

Claims

1. A logging interpretation method based on array induction, characterized in that, it includes: Collect logging data through array induction and perform environmental correction on the logging data to obtain preliminary logging data; Obtain formation dip data; According to the formation dip data, perform dip correction on the preliminary logging data to obtain the corrected logging result; The performing dip correction on the preliminary logging data according to the formation dip data to obtain the corrected logging result includes: The first conductivity in the preliminary logging data is converted by presetting a first formula to obtain a second conductivity that has an associated relationship with the frequency ; Performing inversion processing on the second conductivity using a constraint condition and a preset second formula to obtain the corrected conductivity; wherein, the preset first formula is: ； Among them, , is the propagation form, is the angular frequency, is the conductivity, is the main coil pitch, is the shielding coil pitch, is the apparent conductivity corresponding to the corrected subarray, is the magnetic permeability, , is the imaginary unit, and each subarray corresponds to a , and here the values of all subarrays are put into a vector , is the real part extraction function; The preset second formula is an inversion penalty function, and the inversion penalty function is: ； Among them, is the initial value of the formation conductivity selected The calculated subarray response vector, and the initial value of the formation conductivity is the conductivity after environmental correction. is the th constraint condition, being a natural number.

2. The method according to claim 1, characterized in that, the method further includes: Performing a focusing processing operation on the corrected logging result to obtain an optimized logging result, and the focusing processing operation includes software focusing processing and longitudinal resolution matching.

3. The method according to claim 1, characterized in that, the constraint condition at least includes an interlayer conductivity constraint condition and an inter-subarray conductivity constraint condition; Correspondingly, the performing inversion processing on the second conductivity using a constraint condition and a preset second formula to obtain the corrected conductivity includes: Performing inversion processing on the second conductivity through a preset third formula, an interlayer conductivity constraint condition and an inter-subarray conductivity constraint condition to obtain the corrected conductivity; The preset third formula is: ； Among them, is the interlayer conductivity constraint condition, is the inter-subarray conductivity constraint condition; The interlayer conductivity constraint condition includes: ； Among them, is the number of sub-arrays, is the number of formations, is the th sub-array's conductivity in the th formation, is the constraint value of the th sub-array's conductivity in the th formation and the weighted average of the conductivities of the formations around the th formation, is the coefficient for adding the constraint to the th sub-array's conductivity in the th formation; The inter-subarray conductivity constraint condition includes: ； Among them, is the number of sub-arrays, is the conductivity of the th sub-array, is the conductivity of the th sub-array, is the coefficient that adds a constraint to the conductivity between the th sub-array and the th sub-array.

4. A logging interpretation device based on array induction, characterized in that, it includes: An environmental correction module, an acquisition module and a dip correction module; The environmental correction module is configured to collect logging data through array induction and perform environmental correction on the logging data to obtain preliminary logging data; The acquisition module is configured to obtain formation dip data; The dip correction module is configured to perform dip correction on the preliminary logging data according to the formation dip data to obtain the corrected logging result; The dip correction module is specifically configured to: The first conductivity in the preliminary logging data is converted by presetting a first formula to obtain a second conductivity associated with the frequency ; Performing inversion processing on the second conductivity using a constraint condition and a preset second formula to obtain the corrected conductivity; wherein, the preset first formula is: ； Among them, , is the propagation form, is the circular frequency, is the conductivity, is the main coil spacing, is the shielding coil spacing, is the apparent conductivity corresponding to the corrected subarray, is the magnetic permeability, , is the imaginary unit, and each subarray corresponds to a , and here the values of all subarrays are put into a vector , is the real part extraction function; The preset second formula is an inversion penalty function, and the inversion penalty function is: ； Among them, is the initial value of the formation conductivity selected The calculated sub-array response vector, and the initial value of the formation conductivity is the conductivity after environmental correction. is the th constraint condition, where is a natural number.

5. An electronic device, including a processor, a memory, and a computer program stored on the memory and executable on the processor, characterized in that, when the processor executes the program, it implements the logging interpretation method based on array induction according to any one of claims 1-3 above.

6. A computer-readable storage medium, characterized in that, the computer-readable storage medium stores computer-executable instructions, and when the computer-executable instructions are executed by a processor, they are used to implement the logging interpretation method based on array induction according to any one of claims 1-3.

7. A computer program product, including a computer program, characterized in that, When the computer program is executed by a processor, it is used to implement the array induction-based logging interpretation method according to any one of claims 1-3.

Images