Inclined shaft anisotropic stratum three-dimensional induction logging data inversion method and device

Abstract

The invention provides an inclined shaft anisotropic stratum three-dimensional induction logging data inversion method and device, and the method comprises the steps: carrying out the azimuth conversion of the pre-obtained three-dimensional induction logging data, so as to generate a three-dimensional induction logging curve under a borehole coordinate system; constructing a window function according to the morphological characteristics of the three-dimensional induction logging curve; segmenting the three-dimensional induction logging curve through the window function; all-parameter inversion is carried out on each section of three-dimensional induction logging curve to obtain the horizontal resistivity, the vertical resistivity, the stratum boundary and the hole drift angle of each horizon; wherein the objective function of all-parameter inversion is that the fitting error of the ZZ component of the three-dimensional induction logging curve is minimum. According to the method, the parameters such as the anisotropic stratum conductivity and the borehole inclination angle of the inclined shaft can be quickly and accurately solved, so that the logging information of the three-dimensional induction instrument can be better quickly processed and applied on site.

Description

Technical Field

[0001] This application belongs to the field of oil and gas field exploration and development technology, particularly the field of well logging data processing technology, specifically relating to a method and device for inverting three-dimensional induction logging data of anisotropic formations in deviated wells. Background Technology

[0002] It is estimated that about 30% of the world's oil and gas reserves are located in thin sandstone-mudstone interbedded layers. These thin interbedded reservoirs can be equivalent to macroscopic uniaxial anisotropic strata (or, laterally isotropic strata, abbreviated as TI strata). Therefore, the detection and identification of such strata are of great significance for the development of oil and gas resources.

[0003] For existing axial induction logging instruments, the vertical resolution is not high enough, so in actual production, such thin interbedded oil reservoirs are often mistaken for high water saturation layers and thus missed.

[0004] Three-dimensional induction logging instruments consist of three perpendicular transmitting coils and three parallel receiving coils. They can detect the horizontal and vertical conductivity of the formation, allowing for the identification of formation characteristics from a three-dimensional perspective. This provides an inherent advantage in detecting thin and complex reservoirs. However, in deviated wells, the logging response of three-dimensional induction instruments is generally related to the formation's horizontal and vertical conductivity as well as the wellbore dip angle, exhibiting significant nonlinearity. Furthermore, the influence of adjacent layers on different components of the logging response varies, all of which increase the difficulty of logging data processing, interpretation, and evaluation.

[0005] In existing technologies, the processing of 3D induction logging data mainly relies on multi-parameter nonlinear iterative inversion methods. However, iterative inversion requires multiple forward calculations to fit the logging data and obtain formation parameters. Inverting an entire segment of logging data—including hundreds or even thousands of meters and covering hundreds of strata—is extremely difficult, mainly for the following reasons:

[0006] The computational load is large and the calculation speed is slow: considering that the three-dimensional induction logging response components are numerous and the relationships are complex; at the same time, it is affected by various environmental factors such as formation, dip angle and adjacent layers, and a single inversion requires multiple forward calculations. Furthermore, the iterative inversion also requires calculating and solving the derivative matrix (Jacobi matrix) of the formation parameters.

[0007] Severe multivaluedness and poor stability: The inversion processing of 3D induction logging data is mathematically a local optimization problem, and multivaluedness and instability are inherent properties of it. The hundreds of formation parameters to be solved can easily have equivalent relationships with each other, causing the inversion program to fail to find the true solution during the solution process, resulting in instability and non-convergence of the inversion program. Summary of the Invention

[0008] One objective of this invention is to provide a method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells. This method enables the rapid and accurate determination of parameters such as the electrical conductivity and wellbore dip angle of anisotropic (TI) formations in deviated wells, allowing for better rapid processing and field application of three-dimensional induction logging data.

[0009] Another object of the present invention is to provide a device for inverting three-dimensional induction logging data of anisotropic formations in deviated wells. A further object of the present invention is to provide an electronic device comprising a memory and a processor, wherein the memory stores a computer program, and the processor, when executing the computer program, implements the steps of the above-described method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells. A further object of the present invention is to provide a readable medium storing a computer program thereon, wherein the computer program, when executed by a processor, implements the steps of the above-described method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells.

[0010] To address the technical problems in the background section of this application, the present invention provides the following technical solutions:

[0011] In a first aspect, the present invention provides a method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells, comprising:

[0012] The azimuth angle of the pre-acquired three-dimensional induction logging data is transformed to generate three-dimensional induction logging curves in the wellbore coordinate system;

[0013] Window functions are constructed based on the morphological characteristics of three-dimensional induction logging curves;

[0014] The three-dimensional induction logging curve is segmented using the window function;

[0015] Full-parameter inversion is performed on each segment of the three-dimensional induction logging curve to obtain the horizontal resistivity, vertical resistivity, formation boundary, and well inclination angle of each layer; wherein, the objective function of the full-parameter inversion is to minimize the fitting error of the ZZ component of the three-dimensional induction logging curve.

[0016] In some embodiments of the present invention, the morphological feature is the morphological feature of the three-dimensional induction logging curve at the formation boundary.

[0017] In some embodiments of the present invention, before performing full-parameter inversion on each segment of the three-dimensional induction logging curve, the method further includes:

[0018] The initial value of the formation boundary is determined based on the XZ and ZX components of the three-dimensional induction logging curve.

[0019] In some embodiments of the present invention, determining the initial value of the formation boundary based on the XZ and ZX components of the three-dimensional induction logging curve includes:

[0020] Determine the cross-component response difference between the XZ component and the ZX component;

[0021] The initial value of the formation boundary is determined based on the local maxima and local minima of the response difference of the cross components.

[0022] In some embodiments of the present invention, a method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells further includes:

[0023] Determine the order of performing full-parameter inversion on multiple segments of three-dimensional induction logging curves.

[0024] In some embodiments of the present invention, determining the order of performing full-parameter inversion on multiple segments of three-dimensional induction logging curves includes:

[0025] Based on the position segmentation of the total window and the position segmentation of the main window of the window function corresponding to the current segment of the three-dimensional induction logging curve, the order of full parameter inversion of the current segment of the three-dimensional induction logging curve is determined.

[0026] In some embodiments of the present invention, a method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells further includes:

[0027] The constraints for the full-parameter inversion are generated based on the relative errors of the three-dimensional induction logging data and the relative errors of the formation parameters, and the vertical resistivity is not less than the horizontal resistivity during the full-parameter inversion process.

[0028] Secondly, the present invention provides a device for inverting three-dimensional induction logging data of anisotropic formations in deviated wells, the device comprising:

[0029] The well logging data conversion module is used to convert the azimuth angle of the pre-acquired three-dimensional induction logging data to generate three-dimensional induction logging curves in the wellbore coordinate system.

[0030] The window function construction module is used to construct window functions based on the morphological characteristics of three-dimensional induction logging curves;

[0031] The logging curve segmentation module is used to segment the three-dimensional induction logging curve using the window function;

[0032] The well logging curve inversion module is used to perform full-parameter inversion on each segment of the three-dimensional induction logging curve to obtain the horizontal resistivity, vertical resistivity, formation boundary, and well inclination angle of each layer; wherein, the objective function of the full-parameter inversion is to minimize the fitting error of the ZZ component of the three-dimensional induction logging curve.

[0033] In some embodiments of the present invention, the morphological feature is the morphological feature of the three-dimensional induction logging curve at the formation boundary.

[0034] In some embodiments of the present invention, a three-dimensional induction logging data inversion device for anisotropic formations in deviated wells further includes:

[0035] The formation boundary initial value determination module is used to determine the initial value of the formation boundary based on the XZ component and ZX component of the three-dimensional induction logging curve.

[0036] In some embodiments of the present invention, the formation boundary initial value determination module includes:

[0037] A cross-component response difference determination unit is used to determine the cross-component response difference between the XZ component and the ZX component;

[0038] The formation boundary initial value determination unit is used to determine the initial value of the formation boundary based on the local maxima and local minima of the response difference of the cross components.

[0039] In some embodiments of the present invention, a three-dimensional induction logging data inversion device for anisotropic formations in deviated wells further includes:

[0040] The inversion sequence determination module is used to determine the order in which full-parameter inversions are performed on multiple segments of three-dimensional induction logging curves.

[0041] In some embodiments of the present invention, the inversion order determination module includes:

[0042] The inversion order determination unit is used to determine the order of full parameter inversion of the current segment of the three-dimensional induction logging curve based on the position segmentation of the total window of the window function corresponding to the current segment and the position segmentation of the main window.

[0043] In some embodiments of the present invention, a three-dimensional induction logging data inversion device for anisotropic formations in deviated wells further includes:

[0044] The constraint generation module is used to generate the constraint conditions for the full parameter inversion based on the relative error of the three-dimensional induction logging data and the relative error of the formation parameters, and the vertical resistivity is not less than the horizontal resistivity during the full parameter inversion process.

[0045] Thirdly, the present invention provides a computer program product, including a computer program / instruction, which, when executed by a processor, implements the steps of a method for inverting three-dimensional induction logging data of anisotropic formations in a deviated well.

[0046] Fourthly, the present invention provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the steps of a method for inverting three-dimensional induction logging data of anisotropic formations in a deviated well.

[0047] Fifthly, the present invention provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of a method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells.

[0048] As described above, embodiments of the present invention provide a method and apparatus for inverting three-dimensional induction logging data in anisotropic formations in deviated wells. The corresponding method for inverting three-dimensional induction logging data in anisotropic formations in deviated wells includes: firstly, performing azimuth transformation on pre-acquired three-dimensional induction logging data to generate three-dimensional induction logging curves in the wellbore coordinate system; then, constructing a window function based on the morphological characteristics of the three-dimensional induction logging curves; segmenting the three-dimensional induction logging curves using the window function; and finally, performing full-parameter inversion on each segment of the three-dimensional induction logging curves to obtain the horizontal resistivity, vertical resistivity, formation boundary, and well inclination angle of each layer; wherein, the objective function of the full-parameter inversion is to minimize the fitting error of the ZZ component of the three-dimensional induction logging curve.

[0049] The corresponding anisotropic formation three-dimensional induction logging data inversion device includes: a logging data conversion module, used to convert the azimuth of the pre-acquired three-dimensional induction logging data to generate a three-dimensional induction logging curve in the wellbore coordinate system; a window function construction module, used to construct a window function based on the morphological characteristics of the three-dimensional induction logging curve; a logging curve segmentation module, used to segment the three-dimensional induction logging curve through the window function; and a logging curve inversion module, used to perform full-parameter inversion on each segment of the three-dimensional induction logging curve to obtain the horizontal resistivity, vertical resistivity, formation boundary, and well inclination angle of each layer; wherein, the objective function of the full-parameter inversion is to minimize the fitting error of the ZZ component of the three-dimensional induction logging curve.

[0050] In summary, considering the large volume of measured well logging data in 3D induction logging, this invention performs windowed and segmented inversion on the logging data to be processed, reducing the influence of the surrounding rock and improving the speed and accuracy of 3D induction logging data processing. Furthermore, this invention can rapidly invert 3D induction logging data from anisotropic formations in deviated wells, significantly shortening the logging data processing time and improving the efficiency of the logging instrument. This enables 3D induction logging data to be better applied to field processing and interpretation, and is of great significance for the widespread application of 3D induction logging instruments in the field. Attached Figure Description

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

[0052] Figure 1 This is a flowchart illustrating the method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells, as described in an embodiment of the present invention. Figure 1 ;

[0053] Figure 2 This is a diagram illustrating the coil system structure of a three-dimensional induction logging instrument in an embodiment of the present invention;

[0054] Figure 3 This is a schematic diagram illustrating the relationship between the formation coordinate system, the wellbore coordinate system, and the instrument coordinate system in an embodiment of the present invention. Figure 1 ;

[0055] Figure 4 This is a schematic diagram illustrating the relationship between the formation coordinate system, the wellbore coordinate system, and the instrument coordinate system in an embodiment of the present invention. Figure 2 ;

[0056] Figure 5 This is a flowchart illustrating the method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells, as described in an embodiment of the present invention. Figure 2 ;

[0057] Figure 6 This is a flowchart illustrating step 500 of the method for inverting three-dimensional induction logging data of anisotropic formations in a deviated well, as described in an embodiment of the present invention.

[0058] Figure 7 This is a flowchart illustrating the method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells, as described in an embodiment of the present invention. Figure 3 ;

[0059] Figure 8 This is a flowchart illustrating step 600 of the method for inverting three-dimensional induction logging data of anisotropic formations in a deviated well, as described in an embodiment of the present invention.

[0060] Figure 9 This is a flowchart illustrating the method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells, as described in an embodiment of the present invention. Figure 4 ;

[0061] Figure 10 This is a flowchart illustrating the method for inverting three-dimensional induction logging data of anisotropic formations in a deviated well, according to a specific embodiment of the present invention.

[0062] Figure 11 This is a schematic diagram illustrating the construction of a piecewise inversion window function in a specific embodiment of the present invention;

[0063] Figure 12 This is a schematic diagram of well logging data segmentation in a specific embodiment of the present invention;

[0064] Figure 13This is a graph showing the iterative inversion error in a specific embodiment of the present invention;

[0065] Figure 14 This is a schematic diagram of the well logging data processing results obtained after inversion in a specific embodiment of the present invention.

[0066] Figure 15 A block diagram illustrating a three-dimensional induction logging data inversion device for anisotropic formations in a deviated well, as described in an embodiment of the present invention. Figure 1 ;

[0067] Figure 16 A block diagram illustrating a three-dimensional induction logging data inversion device for anisotropic formations in a deviated well, as described in an embodiment of the present invention. Figure 2 ;

[0068] Figure 17 This is a block diagram of the formation boundary initial value determination module 40 in an embodiment of the present invention;

[0069] Figure 18 A block diagram illustrating a three-dimensional induction logging data inversion device for anisotropic formations in a deviated well, as described in an embodiment of the present invention. Figure 3 ;

[0070] Figure 19 This is a block diagram of the inversion order determination module 60 in an embodiment of the present invention;

[0071] Figure 20 A block diagram illustrating a three-dimensional induction logging data inversion device for anisotropic formations in a deviated well, as described in an embodiment of the present invention. Figure 4 ;

[0072] Figure 21 This is a schematic diagram of the structure of an electronic device in an embodiment of the present invention. Detailed Implementation

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

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

[0075] It should be noted that the terms "comprising" and "having," and any variations thereof, in the specification, claims, and accompanying drawings of this application are intended to cover a non-exclusive inclusion. For example, a process, method, system, product, or apparatus that includes a series of steps or units is not necessarily limited to those explicitly listed, but may include other steps or units not explicitly listed or inherent to these processes, methods, products, or apparatuses. Without conflict, the embodiments and features in the embodiments of this application can be combined with each other. This application will now be described in detail with reference to the accompanying drawings and embodiments.

[0076] The acquisition, storage, use, and processing of data in this application comply with relevant laws and regulations.

[0077] An embodiment of the present invention provides a specific implementation method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells. See [link to implementation details]. Figure 1 The method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells specifically includes the following:

[0078] Step 100: Perform azimuth transformation on the pre-acquired three-dimensional induction logging data to generate three-dimensional induction logging curves in the wellbore coordinate system;

[0079] Step 200: Construct a window function based on the morphological characteristics of the three-dimensional induction logging curve;

[0080] Step 300: Segment the three-dimensional induction logging curve using the window function;

[0081] Step 400: Perform full-parameter inversion on each segment of the three-dimensional induction logging curve to obtain the horizontal resistivity, vertical resistivity, formation boundary, and well inclination angle of each layer; wherein, the objective function of the full-parameter inversion is to minimize the fitting error of the ZZ component of the three-dimensional induction logging curve.

[0082] As described above, this invention provides a method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells, comprising: firstly, performing azimuth transformation on pre-acquired three-dimensional induction logging data to generate three-dimensional induction logging curves in the wellbore coordinate system; then, constructing a window function based on the morphological characteristics of the three-dimensional induction logging curves; segmenting the three-dimensional induction logging curves using the window function; and finally performing full-parameter inversion on each segment of the three-dimensional induction logging curves to obtain the horizontal resistivity, vertical resistivity, formation boundary, and well inclination angle of each layer; wherein, the objective function of the full-parameter inversion is to minimize the fitting error of the ZZ component of the three-dimensional induction logging curve.

[0083] In summary, considering the large volume of measured well logging data in 3D induction logging, this invention performs windowed and segmented inversion on the logging data to be processed, reducing the influence of the surrounding rock and improving the speed and accuracy of 3D induction logging data processing. Furthermore, this invention can rapidly invert 3D induction logging data from anisotropic formations in deviated wells, significantly shortening the logging data processing time and improving the efficiency of the logging instrument. This enables 3D induction logging data to be better applied to field processing and interpretation, and is of great significance for the widespread application of 3D induction logging instruments in the field.

[0084] For step 100, it is understood that the 3D induction logging data is acquired by the 3D induction logging instrument in the wellbore, see [link to relevant documentation]. Figure 2 The coil system of the dimensional induction logging instrument consists of three transmitting coils T, which are centered at a common point and perpendicular to each other. x T y T z Three receiving coils R parallel to it x R y R z (source distance is L1) and three shielding coils B x B y B z Composed of (source distance L2).

[0085] A shielding coil is introduced to counteract the direct coupling component generated by the transmitting coil in the receiving coil. Direct coupling refers to the induced signal directly induced in the receiving coil by the closed-loop characteristic of magnetic flux without passing through a ground layer. The shielding coil and the receiving coil are wound in opposite directions and have unequal numbers of turns, so the direct-coupled electromotive forces generated in the shielding coil and the receiving coil cancel each other out in the air.

[0086] When the transmitting coil system emits sinusoidal alternating current into the surroundings, it can simultaneously measure nine magnetic field components H on the receiving coil system. ij (i=x,y,z; j=x,y,z), where H xyThis represents the magnetic field strength generated by emission in the x direction and reception in the y direction; other components are defined similarly.

[0087] To examine the three-dimensional induction logging response in inclined wells, three coordinate systems need to be introduced, namely, the formation coordinate system OX. f Y f Z f Wellbore coordinate system OX w Y w Z w and instrument coordinate system OX t Y t Z t (like Figure 3 as well as Figure 4 As shown), the magnetic field tensors in these three coordinate systems satisfy the following rotation transformation rule:

[0088]

[0089] in,

[0090]

[0091] In equations (1) and (2), α is the wellbore dip angle, defined as the Z-axis in formation coordinates (Z... f ) and the Z-axis in the wellbore coordinate system (Z w Angle between them; The instrument azimuth angle is defined as the angle between the instrument's projection onto the XY plane and the X-axis. After compensation, the magnetic field strength measured by the receiving coil system... It can be expressed by the following formula (3):

[0092]

[0093] In the formula H ij1 H represents the magnetic field generated by the receiving coil. ij2 This represents the magnetic field generated by the shielding coil. To facilitate comparison between the logging response and formation electrical parameters, the induction logging response is typically normalized to a quantity in the dimension of conductivity. As shown below:

[0094]

[0095] in, K represents the imaginary part of the magnetic field of the coil system. ij For the coil system instrument coefficient,

[0096]

[0097]

[0098] Furthermore, the forward modeling first calculates the electromagnetic field excited by three orthogonal transmitting coils in the formation coordinate system. Then, through a coordinate rotation transformation related to the wellbore inclination angle α, it obtains the three-dimensional induction logging response in the wellbore coordinate system. Finally, through a coordinate rotation transformation related to the instrument azimuth angle φ, it obtains the three-dimensional induction logging response in the instrument coordinate system. Since actual logging data is always obtained in the instrument coordinate system, in data processing, the measured data is inversely transformed by the instrument azimuth angle to obtain the logging response curve in the wellbore coordinate system before further processing.

[0099] For step 200, the morphological characteristics here refer to the morphological characteristics of the three-dimensional induction logging curves at the formation boundaries. That is, the window function is used to reduce the influence of the surrounding rock on the instrument response, and the function shape is constructed based on the instrument response curve shape in the layered formations.

[0100] For step 300, the logging data to be processed is divided into several segments using a window function, and each segment of logging data to be inverted is weighted to reduce the influence of the surrounding rock.

[0101] For step 400, the Gauss-Newton iterative algorithm can be selected for the full parameter inversion. The Gauss-Newton iterative algorithm is a numerical optimization algorithm for nonlinear least squares problems. This nonlinear least squares problem can be formulated as finding a parameter vector that minimizes the sum of squared residuals. The steps of the Gauss-Newton iterative algorithm are as follows:

[0102] Step 1: Initialize the parameter vector, which is generally the initial estimate of the model parameters.

[0103] Step 2: Calculate the residual vector, which is the difference between the actual observed values ​​and the model predictions.

[0104] Step 3: Construct the Jacobian matrix, which contains the partial derivatives of each observation with respect to each parameter.

[0105] Step 4: Calculate the update direction by solving a system of linear equations.

[0106] Step 5: Update the parameter vector.

[0107] Repeat steps 2 to 5 above until the convergence condition is met (e.g., the parameter change is less than a certain threshold).

[0108] In some embodiments of the present invention, the morphological feature is the morphological feature of the three-dimensional induction logging curve at the formation boundary.

[0109] The morphological features in step 200 refer to the morphological features of the three-dimensional induction logging curves corresponding to the bottom boundary.

[0110] In some embodiments of the present invention, see Figure 5Before step 400, a method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells further includes:

[0111] Step 500: Determine the initial value of the formation boundary based on the XZ and ZX components of the three-dimensional induction logging curve.

[0112] In 3D inductive logging technology, a multi-component sensing tool is deployed, capable of measuring electromagnetic field components along different directions. Specifically, the XZ and ZX components are measured separately.

[0113] The XZ component refers to the response of the electromagnetic field generated by the transmitter in the X direction in the Z direction (usually the well axis). This component can provide information about the electrical conductivity of horizontal bedding (the bedding planes perpendicular to the well axis).

[0114] ZX component: This refers to the response of the electromagnetic field generated by the transmitter in the Z direction in the X direction (horizontal direction, possibly the transverse direction of the wellbore). This component helps to assess the conductivity of the layer parallel to the well axis.

[0115] In some embodiments of the present invention, see Figure 6 Step 500 includes:

[0116] Step 501: Determine the cross-component response difference between the XZ component and the ZX component;

[0117] Step 502: Determine the initial value of the formation boundary based on the local maxima and local minima of the cross component response difference.

[0118] Specifically, in steps 501 and 502, the strata are divided according to the local maxima and minima of the cross component response difference (XZ-ZX), and initial values ​​of the strata boundaries are given.

[0119] In some embodiments of the present invention, see Figure 7 A method for inverting three-dimensional induction logging data of anisotropic formations in inclined wells, further comprising:

[0120] Step 600: Determine the order of performing full-parameter inversion on multiple segments of three-dimensional induction logging curves.

[0121] Understandably, for multiple segments of three-dimensional induction logging curves, it is necessary to determine an order and perform full parameter inversion on each segment of the three-dimensional induction logging curve according to that order.

[0122] In some embodiments of the present invention, see Figure 8 Step 600 includes:

[0123] Step 601: Determine the order of full parameter inversion of the current segment of the three-dimensional induction logging curve based on the position segmentation of the total window and the position segmentation of the main window of the window function corresponding to the current segment of the three-dimensional induction logging curve.

[0124] Specifically, given the measurement data to be inverted, it is divided into segments according to the positions of the total window and the main window, and each segment of data is inverted in turn;

[0125] In some embodiments of the present invention, see Figure 9 A method for inverting three-dimensional induction logging data of anisotropic formations in inclined wells, further comprising:

[0126] Step 700: Generate the constraints for the full-parameter inversion based on the relative error of the three-dimensional induction logging data and the relative error of the formation parameters, and ensure that the vertical resistivity is not less than the horizontal resistivity during the full-parameter inversion process.

[0127] During the inversion process, relative errors of well logging response and formation parameters are introduced to monitor the convergence of the inversion program. In addition, certain constraints on formation resistivity and anisotropy coefficients need to be added to the vertical one-dimensional inversion objective function: considering the physical mechanism of equivalent macroscopic anisotropic formations of thin sandstone-mudstone interbedded groups, the vertical resistivity of the formation must be greater than or equal to the horizontal resistivity of the formation during the iterative inversion process.

[0128] As described above, this invention provides a method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells, comprising: firstly, performing azimuth transformation on pre-acquired three-dimensional induction logging data to generate three-dimensional induction logging curves in the wellbore coordinate system; then, constructing a window function based on the morphological characteristics of the three-dimensional induction logging curves; segmenting the three-dimensional induction logging curves using the window function; and finally performing full-parameter inversion on each segment of the three-dimensional induction logging curves to obtain the horizontal resistivity, vertical resistivity, formation boundary, and well inclination angle of each layer; wherein, the objective function of the full-parameter inversion is to minimize the fitting error of the ZZ component of the three-dimensional induction logging curve.

[0129] In summary, considering the large volume of measured well logging data in 3D induction logging, this invention performs windowed and segmented inversion on the logging data to be processed, reducing the influence of the surrounding rock and improving the speed and accuracy of 3D induction logging data processing. Furthermore, this invention can rapidly invert 3D induction logging data from anisotropic formations in deviated wells, significantly shortening the logging data processing time and improving the efficiency of the logging instrument. This enables 3D induction logging data to be better applied to field processing and interpretation, and is of great significance for the widespread application of 3D induction logging instruments in the field.

[0130] In one specific embodiment, the present invention also provides a specific implementation of the method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells, see [link to implementation details]. Figure 10 Specifically, it includes the following steps.

[0131] S1: Acquire the measurement signals from the three-dimensional sensing instrument, wherein the measurement signals include at least: resistivity curve set and well inclination angle.

[0132] Specifically, logging is performed using a three-dimensional induction logging instrument to obtain the three-dimensional induction instrument measurement signal (coil system conductivity curve) and well inclination angle;

[0133] S2: Construct the window function.

[0134] The window function is designed to segment the logging data to be processed into several segments, and each segment of logging data to be inverted is weighted to reduce the influence of the surrounding rock. The main function of the constructed window function is to reduce the influence of the surrounding rock on the instrument response, and the function shape is constructed based on the instrument response curve shape in the layered formation. The length of the main window function is set, preferably one-third of the length of the total window function; specifically:

[0135] Considering the large volume of measured well logging data, the data to be processed is divided into several segments, and windowed inversion is performed on each segment. The presence of surrounding rock will affect the inversion results of the formation parameters of the target layer. Therefore, the following window function is designed:

[0136]

[0137] Weighting is applied to each segment of well logging data to be inverted to reduce the influence of the surrounding rock. The window function expression is shown in equation (6), where w c =1, w f =0.1, L is the length of the main window, z l and z u These are the positions of the top and bottom endpoints of the main window, such as... Figure 11 As shown. Next, the well logging data is divided into three zones according to a window function, as follows: Figure 12 As shown, zones 1 and 3 are the surrounding rock zones, and zone 2 is the main inversion zone. After the logging data inversion for zone 1 is completed, the window function will automatically shift down to the distance of a main window, reaching zone 2. At this point... Figure 12 Zone 2 is designated as the surrounding rock area for this data segment, while Zone 3 is the primary inversion area. After iterative inversion, the inversion results for Zone 3 are output. By moving the window down sequentially, the inversion results for the entire stratigraphic parameters can be obtained.

[0138] S3: Divide the measurement data to be inverted into segments according to the positions of the total window and the main window, and invert the segmented data in sequence.

[0139] Given the measurement data to be inverted, it is divided into segments according to the positions of the total window and the main window, and the inversion is performed on each segment of data in turn;

[0140] Specifically, before the inversion begins, an initial stratigraphic model for the inversion process (step S5) and the measurement data to be inverted are segmented according to the positions of the total window and the main window. Then, each segment of data is inverted. The inversion uses a vertical one-dimensional inversion based on a horizontally layered stratigraphic model to obtain relatively accurate original stratigraphic resistivity, well inclination angle, and the location of the vertical stratigraphic boundaries. An improved Gauss-Newton iterative algorithm is used during the inversion process. The algorithm is as follows:

[0141] Record well logging response data q i (X)=q i (X1,X2,...,X N ), i = 1, 2, ..., M (M is the number of data collection points); X is the parameter to be inverted, q i (X) is a nonlinear function of X, X = X(σ) h ,σ v (d, α) are the parameters to be inverted. N is the number of variables; for example, in a 5-layer model, N represents the 5 horizontal conductivity σ. h 5 vertical conductivity σ v The objective function is defined as follows: The sum of four formation boundary values ​​'d' and one wellbore dip angle 'α', totaling 15 variables.

[0142]

[0143] Where q i (X t ) is the logging response value, which corresponds to X in the absence of measurement noise. t These are the parameter values ​​to be inverted. For ease of writing, a vector matrix can be introduced.

[0144]

[0145]

[0146] The problem of finding the minimum value of the objective function F(X) is usually called a nonlinear least squares problem. Let X... (k) Given the result of the k-th iteration, the modification ΔX of X can be derived using the Gaussian-Newton optimization algorithm. (k) :

[0147] ΔX (k) =-[J T (X (k) )J(X (k) )] -1 J T (X (k) )f(X (k)(10)

[0148] Here, J(X) is the Jacobi matrix of F(X):

[0149]

[0150] Next iteration point X (k+1) for:

[0151] X (k+1) =X (k) +ΔX (k) (12)

[0152] Due to ΔX (k) Is it f i (X) is obtained by performing a Taylor expansion that only retains the first-order terms, only when X... (k) Iterative convergence can only be guaranteed when the solution is sufficiently close to the minimum point of the objective function. However, this is sometimes difficult to guarantee, resulting in poor stability of the Gauss-Newton method solution. The solution is prone to getting trapped in local extrema and may even cause the inversion iteration to fail. To overcome this shortcoming, an N×N dimensional diagonal damping matrix B(β) of the following form is introduced, with the matrix elements of B as follows:

[0153]

[0154] The damping matrix B(β) will gradually approach the identity matrix as the number of iterations increases. That is, the damping effect introduced in this way will automatically decrease and eventually be canceled as the inversion iteration proceeds normally. Then the next iteration point.

[0155] X (k+1) =X (k) +B(β (k) )ΔX (k) (14)

[0156] Among them, it is better for the parameter β to be between 0.5 and 1.

[0157] Introducing relative error of logging response during inversion process Relative error with formation parameters To monitor the convergence of the inversion procedure, the two parameters are defined as follows:

[0158]

[0159]

[0160] In the formula, n is the total number of logging curves, m is the number of measurement points for each logging curve, N is the total number of vertical layers in the formation model, and d i Let R be the boundary location of the i-th layer, α be the wellbore dip angle, and R be the depth of the well. hi Let R be the horizontal resistivity of the i-th stratum.vi Let be the vertical resistivity of the i-th stratum.

[0161] S4: Delineate strata based on the cross-component response difference (XZ-ZX).

[0162] Strata are divided based on the local maxima and minima of the cross component response difference (XZ-ZX), with initial values ​​for the stratigraphic boundaries given.

[0163] S5: The initial stratigraphic model given for the inversion process.

[0164] In addition, the following initial values ​​need to be provided: well inclination angle, formation horizontal resistivity, and formation vertical resistivity;

[0165] S6: Perform full-parameter inversion to obtain the formation horizontal resistivity, vertical resistivity, and well inclination angle.

[0166] The Gauss-Newton iterative algorithm is used for full parameter inversion;

[0167] S7: When the fitting error of the ZZ component is minimized, the iterative inversion ends.

[0168] Specifically, when the fitting error of the ZZ component reaches its minimum, the iterative inversion ends, and the formation horizontal resistivity, vertical resistivity, formation boundary and well inclination angle are obtained.

[0169] In the vertical one-dimensional inversion objective function, certain constraints need to be added to limit the formation resistivity and anisotropy coefficient:

[0170] (1) Considering the physical mechanism of equivalent macroscopic anisotropic strata in thin alternating sandstone and mudstone layers, the vertical resistivity of the strata must be greater than or equal to the horizontal resistivity R during the iterative inversion process. v ≥R h .

[0171] (2) Anisotropy coefficient λ of anisotropic strata 2 =R v / R h The limit is generally set below 20, mainly because three-dimensional induction logging instruments have relatively low sensitivity to vertical resistivity, making it difficult to accurately invert vertical resistivity with anisotropy coefficients greater than 20.

[0172] The following uses a deviated well anisotropic TI formation model as an example to illustrate and verify the method of rapid inversion of three-dimensional induction logging data using windowing functions. Table 1 shows a typical thin interbedded anisotropic (TI) formation model. The three-dimensional induction instrument transmission frequency is f = 25 kHz. A long source-spacing coil system with less influence from the wellbore is preferred: main receiving coil source-spacing L1 = 63 in, shielding coil source-spacing L2 = 49.5 in, where 1 in = 2.54 cm, and the wellbore inclination angle α = 60°. Given the initial inversion values: wellbore inclination angle α0 = 55°, horizontal conductivity Rh0 =1.5R h Vertical conductivity R v0 =1.5R v Considering the actual measurement situation, the simulated response plus 3% random noise is used as the measured signal for inversion.

[0173] Table 1 Anisotropic TI Stratigraphic Model

[0174]

[0175] The variation of inversion iteration error with iteration number n is shown in the figure. Figure 13 In the figure, S is defined by equation (15) as the relative error (%) of logging response fitting; Sp is defined by equation (16) as the relative error (%) of formation parameter fitting; SRh is defined by equation (16) as the relative error (%) of horizontal resistivity fitting; SRv is defined by equation (16) as the relative error (%) of vertical resistivity fitting; Sh is defined by equation (16) as the relative error (%) of longitudinal boundary fitting; and Salpha is defined by equation (16) as the relative error (%) of well inclination angle fitting. Figure 13 It can be seen that the stability and convergence speed of the inversion iteration are very good. Even considering the influence of noise, the relative error of the logging response and the relative error of the formation parameters converge rapidly within 10 steps. The resistivity processing results obtained after iterative inversion are shown in [the figure]. Figure 14 The specific data of the obtained inversion results are shown in Table 2. From Table 2 and... Figure 14 As can be seen, the accuracy of the wellbore dip angle obtained by rapid inversion is very high, and the formation horizontal resistivity R h and vertical resistivity R v The obtained results also match the actual horizontal and vertical conductivity curves of the formation very well, indicating that this rapid inversion method is feasible and yields satisfactory results.

[0176] Table 2. Actual parameters and inversion results of the stratigraphic model.

[0177]

[0178] As described above, the specific embodiments of the present invention provide a method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells, including: firstly, performing azimuth transformation on pre-acquired three-dimensional induction logging data to generate three-dimensional induction logging curves in the wellbore coordinate system; then, constructing a window function based on the morphological characteristics of the three-dimensional induction logging curves; segmenting the three-dimensional induction logging curves through the window function; and finally performing full-parameter inversion on each segment of the three-dimensional induction logging curves to obtain the horizontal resistivity, vertical resistivity, formation boundary, and well inclination angle of each layer; wherein, the objective function of the full-parameter inversion is to minimize the fitting error of the ZZ component of the three-dimensional induction logging curve.

[0179] The above method, considering the large volume of measured well logging data, performs windowed, segmented, one-dimensional inversion on the well logging data to be processed, reducing the influence of the surrounding rock and improving the speed and accuracy of processing three-dimensional induction instrument well logging data. Damping matrices and constraints are introduced during the inversion process to enhance the stability of the inversion procedure.

[0180] In summary, this invention enables rapid inversion of 3D induction logging data from anisotropic formations in deviated wells, significantly shortening data processing time, improving the efficiency of logging instruments, and substantially reducing the cost of on-site testing. This allows 3D induction logging data to be better applied to field processing, interpretation, and evaluation, and is of great significance for the widespread application of 3D induction logging instruments in the field.

[0181] Based on the same inventive concept, this application also provides a device for inverting three-dimensional induction logging data of anisotropic formations in deviated wells, which can be used to implement the method described in the above embodiments, as shown in the following embodiments. Since the principle of solving the problem using the device for inverting three-dimensional induction logging data of anisotropic formations in deviated wells is similar to that of the method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells, the implementation of the device can refer to the implementation of the method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells, and repeated details will not be elaborated further. As used below, the terms "unit" or "module" can refer to a combination of software and / or hardware that performs a predetermined function. Although the system described in the following embodiments is preferably implemented in software, hardware implementation, or a combination of software and hardware, is also possible and contemplated.

[0182] The embodiments of the present invention provide a specific implementation of a three-dimensional induction logging data inversion device for anisotropic formations in deviated wells, capable of realizing a three-dimensional induction logging data inversion method for anisotropic formations in deviated wells, wherein, see... Figure 15 A device for inverting three-dimensional induction logging data of anisotropic formations in inclined wells includes:

[0183] The logging data conversion module 10 is used to convert the azimuth angle of the pre-acquired three-dimensional induction logging data to generate a three-dimensional induction logging curve in the wellbore coordinate system.

[0184] Window function construction module 20 is used to construct window functions based on the morphological characteristics of three-dimensional induction logging curves;

[0185] The logging curve segmentation module 30 is used to segment the three-dimensional induction logging curve using the window function.

[0186] The logging curve inversion module 40 is used to perform full-parameter inversion on each segment of the three-dimensional induction logging curve to obtain the horizontal resistivity, vertical resistivity, formation boundary and well inclination angle of each layer; wherein, the objective function of the full-parameter inversion is to minimize the fitting error of the ZZ component of the three-dimensional induction logging curve.

[0187] In some embodiments of the present invention, the morphological feature is the morphological feature of the three-dimensional induction logging curve at the formation boundary.

[0188] In some embodiments of the present invention, see Figure 16 A device for inverting three-dimensional induction logging data of anisotropic formations in inclined wells, further comprising:

[0189] The formation boundary initial value determination module 50 is used to determine the initial value of the formation boundary based on the XZ component and ZX component of the three-dimensional induction logging curve.

[0190] In some embodiments of the present invention, see Figure 17 The formation boundary initial value determination module 50 includes:

[0191] Cross-component response difference determination unit 50a is used to determine the cross-component response difference between the XZ component and the ZX component;

[0192] The formation boundary initial value determination unit 50b is used to determine the initial value of the formation boundary based on the local maxima and local minima of the response difference of the cross components.

[0193] In some embodiments of the present invention, see Figure 18 A device for inverting three-dimensional induction logging data of anisotropic formations in inclined wells, further comprising:

[0194] The inversion sequence determination module 60 is used to determine the order of full-parameter inversion of multiple three-dimensional induction logging curves.

[0195] In some embodiments of the present invention, see Figure 19 The inversion order determination module 60 includes:

[0196] The inversion order determination unit 60a is used to determine the order of full parameter inversion of the current segment of the three-dimensional induction logging curve based on the position segmentation of the total window of the window function corresponding to the current segment of the three-dimensional induction logging curve and the position segmentation of the main window.

[0197] In some embodiments of the present invention, see Figure 20 A device for inverting three-dimensional induction logging data of anisotropic formations in inclined wells, further comprising:

[0198] The constraint generation module 70 is used to generate the constraint conditions for the full parameter inversion based on the relative error of the three-dimensional induction logging data and the relative error of the formation parameters, and the vertical resistivity is not less than the horizontal resistivity during the full parameter inversion process.

[0199] As described above, this invention provides a device for inverting three-dimensional induction logging data of anisotropic formations in deviated wells, comprising: a logging data conversion module for converting the azimuth of pre-acquired three-dimensional induction logging data to generate a three-dimensional induction logging curve in the wellbore coordinate system; a window function construction module for constructing a window function based on the morphological characteristics of the three-dimensional induction logging curve; a logging curve segmentation module for segmenting the three-dimensional induction logging curve using the window function; and a logging curve inversion module for performing full-parameter inversion on each segment of the three-dimensional induction logging curve to obtain the horizontal resistivity, vertical resistivity, formation boundary, and well inclination angle of each layer; wherein, the objective function of the full-parameter inversion is to minimize the fitting error of the ZZ component of the three-dimensional induction logging curve.

[0200] In summary, considering the large volume of measured well logging data in 3D induction logging, this invention performs windowed and segmented inversion on the logging data to be processed, reducing the influence of the surrounding rock and improving the speed and accuracy of 3D induction logging data processing. Furthermore, this invention can rapidly invert 3D induction logging data from anisotropic formations in deviated wells, significantly shortening the logging data processing time and improving the efficiency of the logging instrument. This enables 3D induction logging data to be better applied to field processing and interpretation, and is of great significance for the widespread application of 3D induction logging instruments in the field.

[0201] This application also provides a specific implementation of an electronic device capable of performing all steps in the above-described method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells. See [link to implementation details]. Figure 21 The electronic devices specifically include the following:

[0202] Processor 1201, memory 1202, communications interface 1203, and bus 1204;

[0203] The processor 1201, memory 1202, and communication interface 1203 communicate with each other via bus 1204; the communication interface 1203 is used to realize information transmission between server-side devices and client-side devices and other related devices.

[0204] The processor 1201 is used to call the computer program in the memory 1202. When the processor executes the computer program, it implements all the steps in the three-dimensional induction logging data inversion method for anisotropic formations in the above embodiment. For example, when the processor executes the computer program, it implements the following steps:

[0205] Step 100: Perform azimuth transformation on the pre-acquired three-dimensional induction logging data to generate three-dimensional induction logging curves in the wellbore coordinate system;

[0206] Step 200: Construct a window function based on the morphological characteristics of the three-dimensional induction logging curve;

[0207] Step 300: Segment the three-dimensional induction logging curve using the window function;

[0208] Step 400: Perform full-parameter inversion on each segment of the three-dimensional induction logging curve to obtain the horizontal resistivity, vertical resistivity, formation boundary, and well inclination angle of each layer; wherein, the objective function of the full-parameter inversion is to minimize the fitting error of the ZZ component of the three-dimensional induction logging curve.

[0209] Embodiments of this application also provide a computer-readable storage medium capable of implementing all steps in the three-dimensional induction logging data inversion method for anisotropic formations in the above embodiments. The computer-readable storage medium stores a computer program, which, when executed by a processor, implements all steps in the three-dimensional induction logging data inversion method for anisotropic formations in the above embodiments. For example, when the processor executes the computer program, it implements the following steps:

[0210] Step 100: Perform azimuth transformation on the pre-acquired three-dimensional induction logging data to generate three-dimensional induction logging curves in the wellbore coordinate system;

[0211] Step 200: Construct a window function based on the morphological characteristics of the three-dimensional induction logging curve;

[0212] Step 300: Segment the three-dimensional induction logging curve using the window function;

[0213] Step 400: Perform full-parameter inversion on each segment of the three-dimensional induction logging curve to obtain the horizontal resistivity, vertical resistivity, formation boundary, and well inclination angle of each layer; wherein, the objective function of the full-parameter inversion is to minimize the fitting error of the ZZ component of the three-dimensional induction logging curve.

[0214] The various embodiments in this specification are described in a progressive manner. Similar or identical parts between embodiments can be referred to interchangeably. Each embodiment focuses on its differences from other embodiments. In particular, hardware + program embodiments are relatively simple in description because they are fundamentally similar to method embodiments; relevant parts can be referred to the descriptions in the method embodiments.

[0215] The foregoing has described specific embodiments of this specification. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recited in the claims may be performed in a different order than that shown in the embodiments and may still achieve the desired result. Furthermore, the processes depicted in the drawings do not necessarily require the specific or sequential order shown to achieve the desired result. In some embodiments, multitasking and parallel processing are possible or may be advantageous.

[0216] While this application provides method operation steps as shown in the embodiments or flowcharts, more or fewer operation steps may be included based on conventional or non-inventive labor. The order of steps listed in the embodiments is merely one possible execution order among many and does not represent the only execution order. In actual device or client product execution, the method can be executed sequentially as shown in the embodiments or drawings, or in parallel (e.g., in a parallel processor or multi-threaded processing environment).

[0217] For ease of description, the above devices are described in terms of function, divided into various modules. Of course, in implementing the embodiments of this specification, the functions of each module can be implemented in one or more software and / or hardware components, or a module that performs the same function can be implemented by a combination of multiple sub-modules or sub-units. The device embodiments described above are merely illustrative. For example, the division of units is only a logical functional division; in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the coupling or direct coupling or communication connection shown or discussed may be through some interfaces, or indirect coupling or communication connection between devices or units, and may be electrical, mechanical, or other forms.

[0218] Those skilled in the art will also know that, besides implementing the controller using purely computer-readable program code, the same functions can be achieved by logically programming the method steps, making the controller function as logic gates, switches, application-specific integrated circuits (ASICs), programmable logic controllers (PLCs), and embedded microcontrollers. Therefore, such a controller can be considered a hardware component, and the devices within it used to implement various functions can also be considered structures within that hardware component. Alternatively, the devices used to implement various functions can be considered as both software modules implementing the method and structures within a hardware component.

[0219] In a typical configuration, a computing device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.

[0220] Memory may include non-persistent storage in computer-readable media, such as random access memory (RAM) and / or non-volatile memory, such as read-only memory (ROM) or flash RAM. Memory is an example of computer-readable media.

[0221] The various embodiments in this specification are described in a progressive manner. Similar or identical parts between embodiments can be referred to mutually. Each embodiment focuses on describing the differences from other embodiments. In particular, system embodiments are basically similar to method embodiments, so the description is relatively simple; relevant parts can be referred to the descriptions in the method embodiments. In the description of this specification, the terms "one embodiment," "some embodiments," "example," "specific example," or "some examples," etc., refer to specific features, structures, materials, or characteristics described in connection with that embodiment or example, which are included in at least one embodiment or example of the embodiments in this specification. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described can be combined in any suitable manner in one or more embodiments or examples. Moreover, without contradiction, those skilled in the art can combine and integrate the different embodiments or examples described in this specification and the features of different embodiments or examples.

[0222] The above description is merely an embodiment of the present specification and is not intended to limit the embodiments of the present specification. For those skilled in the art, various modifications and variations can be made to the embodiments of the present specification. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principle of the embodiments of the present specification should be included within the scope of the claims of the embodiments of the present specification.

Claims

1. A method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells, characterized in that, include: The azimuth angle of the pre-acquired three-dimensional induction logging data is transformed to generate three-dimensional induction logging curves in the wellbore coordinate system; Window functions are constructed based on the morphological characteristics of three-dimensional induction logging curves; The three-dimensional induction logging curve is segmented using the window function; Full-parameter inversion is performed on each segment of the three-dimensional induction logging curve to obtain the horizontal resistivity, vertical resistivity, formation boundary, and well inclination angle of each layer; wherein, the objective function of the full-parameter inversion is to minimize the fitting error of the ZZ component of the three-dimensional induction logging curve.

2. The method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells according to claim 1, characterized in that, The morphological features refer to the morphological features of the three-dimensional induction logging curves at the formation boundary.

3. The method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells according to claim 1, characterized in that, Before performing full-parameter inversion on each segment of the three-dimensional induction logging curve, the following is also included: The initial value of the formation boundary is determined based on the XZ and ZX components of the three-dimensional induction logging curve.

4. The method for inverting three-dimensional induction logging data of anisotropic formations in inclined wells according to claim 3, characterized in that, Based on the XZ and ZX components of the three-dimensional induction logging curve, the initial values ​​of the formation boundary are determined, including: Determine the cross-component response difference between the XZ component and the ZX component; The initial value of the formation boundary is determined based on the local maxima and local minima of the response difference of the cross components.

5. The method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells according to claim 1, characterized in that, Also includes: Determine the order of performing full-parameter inversion on multiple segments of three-dimensional induction logging curves.

6. The method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells according to claim 5, characterized in that, The process of determining the order for full-parameter inversion of multiple three-dimensional induction logging curves includes: Based on the position segmentation of the total window and the position segmentation of the main window of the window function corresponding to the current segment of the three-dimensional induction logging curve, the order of full parameter inversion of the current segment of the three-dimensional induction logging curve is determined.

7. The method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells according to any one of claims 1 to 6, characterized in that, Also includes: The constraints for the full-parameter inversion are generated based on the relative errors of the three-dimensional induction logging data and the relative errors of the formation parameters, and the vertical resistivity is not less than the horizontal resistivity during the full-parameter inversion process.

8. A device for inverting three-dimensional induction logging data of anisotropic formations in inclined wells, characterized in that, include: The well logging data conversion module is used to convert the azimuth angle of the pre-acquired three-dimensional induction logging data to generate three-dimensional induction logging curves in the wellbore coordinate system. The window function construction module is used to construct window functions based on the morphological characteristics of three-dimensional induction logging curves; The logging curve segmentation module is used to segment the three-dimensional induction logging curve using the window function; The well logging curve inversion module is used to perform full-parameter inversion on each segment of the three-dimensional induction logging curve to obtain the horizontal resistivity, vertical resistivity, formation boundary, and well inclination angle of each layer; wherein, the objective function of the full-parameter inversion is to minimize the fitting error of the ZZ component of the three-dimensional induction logging curve.

9. The device for inverting three-dimensional induction logging data of anisotropic formations in inclined wells according to claim 8, characterized in that, The morphological features refer to the morphological features of the three-dimensional induction logging curves at the formation boundary.

10. The device for inverting three-dimensional induction logging data of anisotropic formations in inclined wells according to claim 8, characterized in that, Also includes: The formation boundary initial value determination module is used to determine the initial value of the formation boundary based on the XZ component and ZX component of the three-dimensional induction logging curve.

11. The device for inverting three-dimensional induction logging data of anisotropic formations in inclined wells according to claim 10, characterized in that, The initial value determination module for the formation boundary includes: A cross-component response difference determination unit is used to determine the cross-component response difference between the XZ component and the ZX component; The formation boundary initial value determination unit is used to determine the initial value of the formation boundary based on the local maxima and local minima of the response difference of the cross components.

12. The device for inverting three-dimensional induction logging data of anisotropic formations in inclined wells according to claim 8, characterized in that, Also includes: The inversion sequence determination module is used to determine the order in which full-parameter inversions are performed on multiple segments of three-dimensional induction logging curves.

13. The device for inverting three-dimensional induction logging data of anisotropic formations in inclined wells according to claim 12, characterized in that, The inversion order determination module includes: The inversion order determination unit is used to determine the order of full parameter inversion of the current segment of the three-dimensional induction logging curve based on the position segmentation of the total window of the window function corresponding to the current segment of the three-dimensional induction logging curve and the position segmentation of the main window.

14. The device for inverting three-dimensional induction logging data of anisotropic formations in deviated wells according to any one of claims 8 to 13, characterized in that, Also includes: The constraint generation module is used to generate the constraint conditions for the full parameter inversion based on the relative error of the three-dimensional induction logging data and the relative error of the formation parameters, and the vertical resistivity is not less than the horizontal resistivity during the full parameter inversion process.

15. A computer program product comprising a computer program / instructions, characterized in that, When the computer program / instruction is executed by the processor, it implements the steps of the method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells as described in any one of claims 1 to 7.

16. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the program, it implements the steps of the method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells as described in any one of claims 1 to 7.

17. A computer-readable storage medium having a computer program stored thereon, characterized in that, When executed by a processor, the computer program implements the steps of the method for inverting three-dimensional induction logging data of anisotropic formations in deviated wells as described in any one of claims 1 to 7.

Images