CN115773101A - Array induction borehole correction method, device, storage medium and processor - Google Patents

Abstract

The invention discloses an array induction borehole correction method, an array induction borehole correction device, a storage medium and a processor. The method comprises the following steps: obtaining a measurement signal obtained by measuring each subarray of the array induction logging instrument; determining borehole environment parameters corresponding to the borehole environment of the target stratum, and acquiring a target borehole geometric factor of the borehole environment from a borehole geometric factor library according to the borehole environment parameters by adopting a multidimensional spatial interpolation method; establishing a chromatography matrix equation based on the radial stratum model and the measurement signal, and solving the chromatography matrix equation to obtain the radial conductivity distribution of the target stratum; determining a measurement value expression according to the radial stratum model and the radial conductivity distribution, and determining a borehole correction formula by adopting a target borehole geometric factor corresponding to a borehole environment; and determining the corrected measured value of the target stratum based on the measured value expression and the borehole correction formula, and completing borehole correction. The invention solves the technical problem of borehole correction errors caused by neglecting mud invasion in the conventional borehole correction method.

Description

Array induction borehole correction method, device, storage medium and processor

Technical Field

The invention relates to the technical field of well logging data processing, in particular to an array induction borehole correction method, an array induction borehole correction device, a storage medium and a processor.

Background

The measurement signals of each sub-array of the array induction logging instrument are greatly influenced by the borehole environment, if the influence cannot be corrected correctly, borehole correction residual errors are propagated and even amplified in the subsequent data processing processes such as soft focusing, and the like, so that the reliability of the final measurement curve is reduced.

The existing array induction borehole correction often has poor effect, and field engineers or processing personnel often need to repeatedly adjust borehole correction processing control parameters for processing so as to obtain a satisfactory final logging curve, but the obtained measurement result is poor. When invasion exists, invasion response is taken as well bore influence and is corrected, curve difference becomes small after well bore correction, the capability of curve radial invasion analysis is reduced, and the calculation of the true resistivity of the stratum is also influenced.

The present invention proposes an effective solution to the above-mentioned problems.

Disclosure of Invention

The embodiment of the invention provides an array induction borehole correction method, an array induction borehole correction device, a storage medium and a processor, which are used for at least solving the technical problem of borehole correction errors caused by neglecting mud invasion in the existing borehole correction method.

According to an aspect of an embodiment of the present invention, there is provided an array induction borehole correction method, including: obtaining a measurement signal obtained by measuring each subarray of the array induction logging instrument, wherein the measurement signal at least comprises: a conductivity curve group; determining a borehole environment parameter corresponding to a borehole environment of a target stratum, and acquiring a target borehole geometric factor of the borehole environment from a borehole geometric factor library according to the borehole environment parameter by adopting a multi-dimensional spatial interpolation method, wherein the borehole environment parameter comprises: well bore cal, mud conductivity σ _m Formation conductivity σ _t And an eccentricity ecc, said wellbore geometry factor database being a pre-established wellbore geometry factor database, said wellbore geometry factor being a function of said wellbore environment parameter; establishing a chromatography matrix equation based on the radial stratum model and the measurement signal, and solving the chromatography matrix equation to obtain the target areaA radial conductivity distribution of the layer, wherein the radial formation model is used to characterize a conductivity change of the target formation; determining a measurement value expression according to the radial stratum model and the radial conductivity distribution, and determining a borehole correction formula by adopting the target borehole geometric factor corresponding to the borehole environment; and determining the corrected measured value of the target stratum based on the measured value expression and the borehole correction formula, and completing borehole correction.

Optionally, before obtaining the target borehole geometry factor of the borehole environment from the borehole geometry factor library according to the borehole environment parameter by using a multidimensional spatial interpolation method, the method further includes: determining the variation range of the well diameter, and determining a preset number of discrete points of the well diameter based on the variation range of the well diameter to form a first-dimension discrete point set; determining the variation range of the mud conductivity, and determining a preset number of discrete points of the mud conductivity based on the variation range of the mud conductivity to form a second-dimension discrete point set; determining the variation range of the stratum conductivity, and determining a preset number of discrete points of the stratum conductivity based on the variation range of the stratum conductivity to form a third-dimension discrete point set; determining the variation range of the eccentricity, and determining a preset number of discrete points of the eccentricity based on the variation range of the eccentricity to form a fourth-dimension discrete point set; and combining the first-dimension discrete point set, the second-dimension discrete point set, the third-dimension discrete point set and the fourth-dimension discrete point set, calculating to obtain the borehole geometric factors of each subarray, and storing the borehole geometric factors of all the subarrays according to a preset storage rule to construct the borehole geometric factor library.

Optionally, the determining a wellbore environment parameter corresponding to a wellbore environment of the target formation, and obtaining a target wellbore geometry factor of the wellbore environment from a wellbore geometry factor library according to the wellbore environment parameter by using a multidimensional spatial interpolation method, includes: determining a value of a borehole environment parameter corresponding to the borehole environment of the target formation; acquiring the closest target discrete point from the borehole geometric factor library according to the value of the borehole environmental parameter; and determining the target borehole geometric factor by adopting the multi-dimensional space interpolation method based on the first-dimension discrete point set, the second-dimension discrete point set, the third-dimension discrete point set and the fourth-dimension discrete point set.

Optionally, before establishing a chromatography matrix equation based on the radial formation model and the measurement signal and solving the chromatography matrix equation to obtain the radial conductivity distribution of the target formation, the method further includes: dividing the target stratum into a plurality of layers in the horizontal radius direction, wherein the borehole is the innermost layer, and the stratum except the borehole is an invasion flushing zone, an invasion transition zone and an undisturbed stratum; and constructing the radial formation model using the wellbore, the invaded wash zone, the invaded transition zone, and the undisturbed formation.

Optionally, establishing a chromatography matrix equation based on the radial formation model and the measurement signal, and solving the chromatography matrix equation to obtain the radial conductivity distribution of the target formation, including: determining the borehole geometric factors of each subarray and the radial geometric factors of each subarray; establishing the chromatography matrix equation based on the radial formation model and the measurement signals, the borehole geometry factor and the radial geometry factor, wherein the left term of the chromatography matrix equation is the measurement signals of each subarray, the unknown quantity of the chromatography matrix equation is the conductivity difference between radially adjacent layers, the coefficient matrix of the chromatography matrix equation is composed of the borehole geometry factor, the radial geometry factor and a constant of each subarray, and the constraint condition of the chromatography matrix equation is determined according to the borehole environment parameters; and solving the chromatography matrix equation to obtain the radial conductivity distribution of the target stratum.

Optionally, the determining the corrected measured value of the target formation based on the measured value expression and the borehole correction formula to complete borehole correction includes: processing the measurement signal by using the measurement expression and the borehole correction formula to determine the target formation conductivity and the target borehole mud conductivity of the target formation; and performing borehole correction on the target measurement value by using the target formation conductivity and the target borehole mud conductivity.

According to another aspect of an embodiment of the present invention, there is provided an array induction borehole correction device, including: the acquisition module is used for acquiring measurement signals obtained by measuring each subarray of the array induction logging instrument, wherein the measurement signals at least comprise: a conductivity curve group; the first determining module is used for determining borehole environment parameters corresponding to the borehole environment of the target stratum and acquiring the target borehole geometric factors of the borehole environment from a borehole geometric factor library according to the borehole environment parameters by adopting a multidimensional space interpolation method, wherein the borehole environment parameters comprise: well bore cal, mud conductivity σ _m Formation conductivity σ _t And an eccentricity ecc, said wellbore geometry factor database being a pre-established wellbore geometry factor database, said wellbore geometry factor being a function of said wellbore environment parameter; the processing module is used for establishing a chromatography matrix equation based on a radial stratum model and the measuring signal, and solving the chromatography matrix equation to obtain radial conductivity distribution of the target stratum, wherein the radial stratum model is used for representing conductivity change of the target stratum; a second determining module, configured to determine a measurement expression according to the radial formation model and the radial conductivity distribution, and determine a borehole correction formula by using the target borehole geometry factor corresponding to the borehole environment; and the correction module is used for determining the corrected measured value of the target stratum based on the measured value expression and the borehole correction formula so as to finish borehole correction.

According to another aspect of embodiments of the present invention, there is also provided a non-volatile storage medium having stored thereon a plurality of instructions adapted to be loaded by a processor and to perform any one of the array induced borehole correction methods described above.

According to another aspect of embodiments of the present invention, there is also provided a processor for executing a program, wherein the program is configured to execute any one of the array induced borehole correction methods described above when executed.

According to another aspect of the embodiments of the present invention, there is also provided an electronic device, including a memory and a processor, wherein the memory stores a computer program, and the processor is configured to execute the computer program to perform any one of the array induction borehole correction methods described above.

In an embodiment of the present invention, a measurement signal obtained by measuring with an array induction logging tool is obtained, where the measurement signal at least includes: a resistivity curve group; determining borehole environment parameters, establishing a borehole geometric factor library based on the borehole environment parameters, and determining a target borehole geometric factor of a target stratum based on the borehole geometric factor library; processing the measurement signal by using a radial stratum model to determine a measurement value expression, and processing the measurement signal and the target borehole geometric factor by using the radial stratum model to determine a borehole correction formula; and determining the corrected measured value of the target stratum based on the measured value expression and the borehole correction formula, completing borehole correction, and achieving the purpose of constructing a new adaptive borehole environment parameter solving algorithm, so that the technical effect of accurately obtaining the true conductivity of the stratum is realized, and the technical problem of borehole correction errors caused by neglecting mud invasion in the conventional borehole correction method is solved.

Drawings

The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the invention and together with the description serve to explain the invention and do not constitute a limitation of the invention. In the drawings:

FIG. 1 is a flow chart of a method for array-induced borehole correction according to an embodiment of the present invention;

FIG. 2 is a schematic diagram of an alternative array induction logging instrument probe configuration in accordance with embodiments of the present invention;

FIG. 3 is a graphical representation of an alternative exemplary array induction logging tool log in accordance with an embodiment of the present invention;

FIG. 4 is a schematic illustration of an alternative array induction tool logging in accordance with an embodiment of the present invention;

FIG. 5 is a schematic illustration of an alternative array induced borehole geometry factor curve according to an embodiment of the present invention;

FIG. 6 is a schematic view of an alternative radial two-layer model according to an embodiment of the invention;

FIG. 7 is a schematic view of an alternative radial trilayer model in accordance with an embodiment of the present invention;

FIG. 8 shows a graph of σ according to an embodiment of the present invention _m >σ _xo >σ _t Then, the visual conductivity diagram after model correction;

FIG. 9 shows a graph of σ according to an embodiment of the present invention _t >σ _xo ≥σ _m Then, the visual conductivity diagram after model correction;

FIG. 10 shows a graph of σ according to an embodiment of the present invention _t >σ _m >σ _xo Then, the visual conductivity diagram after model correction;

FIG. 11 shows a graph of σ according to an embodiment of the present invention _m ≥σ _t >σ _xo Then, the model is corrected to obtain a visual conductivity diagram;

FIG. 12 is a schematic view of an alternative intrusion transition path according to an embodiment of the invention;

FIG. 13 is a schematic view of a conductivity radial profile representative of an alternative borehole corrected pre/post array induction curve in accordance with embodiments of the present invention;

FIG. 14 is a schematic diagram of a wellbore correction device according to an embodiment of the present invention.

Detailed Description

In order to make the technical solutions of the present invention better understood, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.

It should be noted that the terms "first," "second," and the like in the description and claims of the present invention and in the drawings described above are used for distinguishing between similar elements and not necessarily for describing a particular sequential or chronological order. It is to be understood that the data so used is interchangeable under appropriate circumstances such that the embodiments of the invention described herein are capable of operation in other sequences than those illustrated or described herein. Furthermore, the terms "comprises," "comprising," and "having," and any variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed, but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.

Example 1

While a method of array-induced borehole correction is provided according to an embodiment of the present invention, it should be noted that the steps illustrated in the flowchart of the figure may be performed in a computer system such as a set of computer-executable instructions, and that while a logical order is illustrated in the flowchart, in some cases the steps illustrated or described may be performed in an order different than presented herein.

In the prior art, the array induction borehole correction method has two defects: (1) Establishing a borehole correction model which is formed by adding infinite uniform stratums to a borehole, is inconsistent with an actual stratum and is inconsistent with the original intention of array induction design; (2) In the self-adaptive borehole correction, the internal correlation between the array induction subarray measurement values is poor, and the method is a qualitative concept and does not establish quantitative description.

Aiming at the theoretical defects and pain points of the existing array induction borehole correction method, a chromatography concept can be applied, a plurality of step functions are introduced to approach any radial one-dimensional stratum, then approximate processing is carried out according to an approximation method (Born), a mathematical expression with an inline relation is deduced, and the method is applied to borehole environment parameters for self-adaptive solution. And a borehole correction formula is reconstructed based on the invasion model, so that the phenomenon that invasion influence is taken as borehole influence and is corrected is avoided, the curve difference after borehole correction is kept, and the invasion radial analysis capability of the curve is improved.

FIG. 1 is a flow chart of a method for array-induced borehole correction according to an embodiment of the present invention, as shown in FIG. 1, the method comprising the steps of:

step S102, obtaining a measurement signal obtained by measuring each subarray of the array induction logging instrument, wherein the measurement signal at least comprises: a conductivity curve group;

step S104, determining borehole environment parameters corresponding to the borehole environment of the target stratum, and acquiring the target borehole geometric factors of the borehole environment from a borehole geometric factor library according to the borehole environment parameters by adopting a multi-dimensional spatial interpolation method, wherein the borehole environment parameters comprise: well bore cal, mud conductivity σ _m Stratum conductivity sigma _t And an eccentricity ecc, said wellbore geometry factor database being a pre-established wellbore geometry factor database, said wellbore geometry factor being a function of said wellbore environment parameter;

step S106, establishing a chromatography matrix equation based on a radial stratum model and the measurement signal, and solving the chromatography matrix equation to obtain radial conductivity distribution of the target stratum, wherein the radial stratum model is used for representing conductivity change of the target stratum;

step S108, determining a measurement value expression according to the radial stratum model and the radial conductivity distribution, and determining a borehole correction formula by adopting the target borehole geometric factor corresponding to the borehole environment;

and step S110, determining the corrected measured value of the target stratum based on the measured value expression and the borehole correction formula, and completing borehole correction.

In an embodiment of the present invention, an execution main body of the array induction borehole correction method in steps S102 to S110 is an array induction borehole correction system, and the system is adopted to obtain a measurement signal obtained by measuring each sub-array of the array induction logging instrument, where the measurement signal at least includes: a conductivity curve group; determining a borehole corresponding to a borehole environment of a target formationAnd acquiring a target borehole geometric factor of the borehole environment from a borehole geometric factor library according to the borehole environmental parameter by adopting a multidimensional spatial interpolation method, wherein the borehole environmental parameter comprises: well bore cal, mud conductivity σ _m Stratum conductivity sigma _t And an eccentricity ecc, said wellbore geometry factor database being a pre-established wellbore geometry factor database, said wellbore geometry factor being a function of said wellbore environment parameter; establishing a chromatography matrix equation based on a radial stratum model and the measurement signal, and solving the chromatography matrix equation to obtain radial conductivity distribution of the target stratum, wherein the radial stratum model is used for representing conductivity change of the target stratum; determining a measurement value expression according to the radial stratum model and the radial conductivity distribution, and determining a borehole correction formula by adopting the target borehole geometric factor corresponding to the borehole environment; and determining the corrected measured value of the target stratum based on the measured value expression and the borehole correction formula, and completing borehole correction.

It should be noted that the array induction logging tool includes a plurality of sub-measurement arrays, functionally characterized by a set of curves describing the depth of multi-path exploration of the invasion profile. The resulting set of curves may include, but is not limited to: a set of conductivity curves, a set of resistivity curves, etc. As shown in the schematic diagram of the probe structure of the array induction logging instrument in FIG. 2, the probe is composed of a plurality of (5 or 6 or 7 or more) sub-arrays (for example, A1 to A6 in the figure). Each sub-array is formed by a transmitter coil and a series of receiver coils (at least 2 coils, a main receiver coil and a shield receiver coil) connected together in series. One transmit coil is shared between the sub-arrays.

As an alternative embodiment, a series of signal processing is performed on raw measurement data of a plurality of sub-arrays of the array induction logging instrument to obtain a final measurement signal; the series of signal processing described above includes: scale processing, temperature effect correction, skin effect correction, borehole correction, software focusing, resolution matching, radial inversion processing, and the like.

As an alternative embodiment, after obtaining the measurement signal measured by the array induction logging tool, the measurement signal obtained may be subjected to certain pre-processing, for example: missing data processing, erroneous data processing, etc. As shown in fig. 3, a typical array induction logging tool log graph, raw measurement signals of the array induction logging tool are subjected to data processing steps such as preprocessing, skin effect correction, borehole correction and software focusing, and resistivity curve groups with various longitudinal resolutions (1 foot, 2 feet, 4 feet or 0.5 foot) and various radial detection depths (10 in, 20in, 30in, 60in, 90in and/or 120 in) can be obtained.

In an alternative embodiment, before obtaining the target borehole geometry factor of the borehole environment from the borehole geometry factor library according to the borehole environment parameter by using multidimensional spatial interpolation, the method further comprises: determining the variation range of the well diameter, and determining a preset number of discrete points of the well diameter based on the variation range of the well diameter to form a first-dimension discrete point set; determining the variation range of the mud conductivity, and determining a preset number of discrete points of the mud conductivity based on the variation range of the mud conductivity to form a second-dimension discrete point set; determining the variation range of the stratum conductivity, and determining a preset number of discrete points of the stratum conductivity based on the variation range of the stratum conductivity to form a third-dimension discrete point set; determining the variation range of the eccentricity, and determining a preset number of discrete points of the eccentricity based on the variation range of the eccentricity to form a fourth-dimension discrete point set; and performing combined calculation on the first-dimension discrete point set, the second-dimension discrete point set, the third-dimension discrete point set and the fourth-dimension discrete point set to obtain the borehole geometric factors of the sub-arrays, and storing the borehole geometric factors of all the sub-arrays according to a preset storage rule to construct the borehole geometric factor library.

In the embodiment of the invention, the size of the borehole geometry factor value of a certain sub-array of the array sensing instrument is determined by the size of the four borehole environment parameter values. In order to obtain the value of the geometric factor of each sub-array sensed by the array corresponding to a specific borehole environment, a borehole geometric factor library can be established in advance.

Note that the array induction tool logging schematic, σ, shown in FIG. 4 _t Is the formation conductivity (the inverse of which is the formation resistivity R) _t )，σ _m As the borehole mud conductivity (inverse of mud resistivity R) _m ) Cal is well diameter, d _tool The diameter of the instrument and x is the distance between the instrument and the well wall. For convenience, the relative eccentricity of the instrument, ecc, is introduced:

when the instrument is centered, the relative eccentricity is 0; when the instrument is completely attached to the borehole wall, the relative eccentricity is 1.

In the embodiment of the invention, according to the induction logging principle, the array induction borehole influence is controlled by four environmental parameters, which are respectively: well diameter cal, mud conductivity σ _m (or its inverse, mud resistivity R _m ) Eccentricity of the tool in the borehole, ecc, and formation conductivity, σ _t (or its inverse, formation resistivity R) _t ). And determining the four environmental parameters as the borehole environmental parameters.

As an optional embodiment, the variation range of the four borehole environment parameter values is determined, the number of discrete points is determined based on the variation range, four discrete point sets are formed, the discrete point sets are combined, all combinations are calculated to obtain the borehole geometric factors of each sub-array, and the borehole geometric factors of all sub-arrays corresponding to all combinations are stored according to a certain mode to form the borehole geometric factor library. For example, the well diameter is set up one discrete point every 0.5 inches from 4 inches to 32 inches, and 57 discrete points are obtained in total; setting a discrete point every 0.05 in the eccentricity range from 0 to 0.95 to obtain 20 discrete points; selecting 25 discrete points from the conductivity of the mud between 0.0001S/m and 1000S/m; selecting 20 discrete points from 0.0001S/m to 100S/m in formation conductivity; there are 57x20x25x20 discrete points in the four-dimensional discrete point set, and for a 7 sub-array of array sensing instruments, there are 57x20x25x20x7 borehole geometry factor values in the borehole geometry factor library.

Optionally, for the convenience and rapidness of borehole correction, according to the structure of the array induction instrument probe coil system and other relevant construction parameters, the borehole geometric factor can be calculated in advance: at four wellbore environment parameters (ecc, cal, σ) _m ，σ _t ) And (4) performing dispersion, calculating and storing the geometric factors of each subarray for all the dispersion combination points, and forming a borehole geometric factor library. The schematic diagram of the array-induced borehole geometry factor curve shown in fig. 5 corresponds to Rt =100ohmm, rm =0.1ohmm, and the eccentricity ec is equal to 0.0,0.5,0.75,0.85,0.90, and 0.95, respectively. The figure shows that: (1) The borehole geometric factors of the shorter 3 sub-arrays are sensitive to the borehole environment parameters and are relatively reliable; (2) Relative to the conductivity variable (. Sigma.) _m ，σ _t ) The wellbore geometry factor is more sensitive to geometric variables (ecc, cal).

It should be noted that in wells with high contrast in mud and formation conductivity (e.g., mud/formation conductivity ratio greater than 100), the majority of the shallow probing subarrays (subarray A1, subarray A2, subarray A3) readings are from the borehole. In large boreholes with high formation resistivities, the borehole geometry factor may reach above 0.6, the measurement signal is almost completely contributed by the borehole, and the measurement error of the total reading may be greater than the formation signal, which may be prone to overcorrection or undercorrection. Therefore, the effects of the borehole environment must be eliminated as much as possible prior to the synthetic focus process.

In an optional embodiment, the determining the wellbore environment parameter corresponding to the wellbore environment of the target formation, and obtaining the target wellbore geometry factor of the wellbore environment from the wellbore geometry factor library according to the wellbore environment parameter by using a multidimensional spatial interpolation method includes: determining a value of a borehole environment parameter corresponding to the borehole environment of the target formation; according to the value of the borehole environment parameter, obtaining the closest target discrete point from the borehole geometry factor library; and determining the target borehole geometric factor by adopting the multi-dimensional space interpolation method based on the first-dimension discrete point set, the second-dimension discrete point set, the third-dimension discrete point set and the fourth-dimension discrete point set.

As an alternative embodiment, the values of 4 borehole environment parameters corresponding to the target formation are determined; and acquiring target discrete points closest to the values of 4 borehole environment parameters corresponding to the target stratum from the borehole geometry factor library, and determining the target borehole geometry factor of the target stratum from the borehole geometry factor library by adopting a spatial interpolation method based on the 4 borehole environment parameters.

In an optional embodiment, before establishing a chromatography matrix equation based on a radial formation model and the measurement signal, and solving the chromatography matrix equation to obtain a radial conductivity distribution of the target formation, the method further includes: dividing the target stratum into a plurality of layers in the radius direction, wherein the borehole is the innermost layer, and the stratums except the borehole are an invasion flushing zone, an invasion transition zone and an undisturbed stratum; constructing the radial formation model using the wellbore, the invaded wash zone, the invaded transition zone, and the undisturbed formation.

In an alternative embodiment, establishing a chromatography matrix equation based on a radial formation model and the measurement signal, and solving the chromatography matrix equation to obtain a radial conductivity distribution of the target formation includes: determining the borehole geometry factor of each subarray and the radial geometry factor of each subarray; establishing the chromatography matrix equation based on the radial formation model and the measurement signals, the borehole geometry factor and the radial geometry factor, wherein the left term of the chromatography matrix equation is the measurement signals of each subarray, the unknown quantity of the chromatography matrix equation is the conductivity difference between radially adjacent layers, the coefficient matrix of the chromatography matrix equation is composed of the borehole geometry factor, the radial geometry factor and a constant of each subarray, and the constraint condition of the chromatography matrix equation is determined according to the borehole environment parameters; and solving the chromatography matrix equation to obtain the radial conductivity distribution of the target stratum.

In the embodiment of the invention, the array senses that the original measurement data meet a certain relation, and the response of a subarray can be judged according to the trend formed between fuzzy subarraysThere is no significant anomaly. The tendency to form between subarrays is governed by two factors: radial variation of formation conductivity and radial geometry factor of the sub-array. After skin-seeking correction, each subarray measurement may be considered to satisfy the Born approximation. Under Born approximation, the response of sub-array i at a depth of invasion of D can be expressed by: sigma _a (i，D)＝σ _t + Δ σ (D) × GF (i, D); wherein, delta sigma (D) is step function, and the step amplitude is sigma _xo -σ _t The step position is D; GF (i, D) is the radial integral geometry for sub-array i at radial depth D. When the radial variation in formation conductivity is not a single step, the radial tomography concept can be approximated with a multiple radial step diagram as shown in fig. 7. Step function delta sigma _l (r _l ) Denotes the l-th step with a step amplitude of Δ σ _l ＝σ _l -σ _l-1 Step position r _l . After introducing the step function, the conductivity σ of the formation can be expressed as:

σ＝σ ₀ +Δσ ₁ (r ₁ )+Δσ ₂ (r ₂ )+…+Δσ _K (r _K )

alternatively, the measurement value for the sub-array i using Born approximation can be expressed as follows:

σ _a (i)＝σ ₀ +Δσ ₁ *GF(i，r ₁ )+…+Δσ _i *GF(i，r _i )+…+Δσ _k *GF(i，r _k )

optionally, the matrix form is used to describe the array induction subarray measurement value of the radial multi-order formation, so as to obtain an internal quantitative correlation general formula for describing the array induction subarray measurement value:

wherein M is the number of subarrays; k is the number of analytic layers; chromatographic coefficient G _ij Corresponding to the radial integral geometric factor, and G _i0 ＝1。

As an alternative embodiment, after skin correction, the array senses each subarrayThe column measurements have good linearity, and the subarray measurement is equal to the sum of the contributions of the regions of space around the instrument, which can be expressed by the formula: sigma _a ＝∫∫∫g _j σ _j dv _j . Wherein, g _j Called space volume unit dv _j The differential geometry factor of (2). The differential geometric factor satisfies the normalization condition, i.e. the integral [ integral ] g ] _j dv _j =1. The differential geometry factor is a function of spatial coordinates. Under a cylindrical coordinate system, the independent variables of the geometric factors are z, r,

optionally, the radial integral geometric factor GF (i, D) corresponding to any radial depth D may be determined by selecting N pre-selected radial depths D _j Geometric factor GF (i, D) of point (j =1,. Cndot., N) _j ) Interpolation to approximate the calculation, i.e.:

optionally, the above formula is substituted into a measurement value expression of the i-th sub-array obtained by Born approximation, so as to obtain:

wherein, note a _M+1 ＝σ ₀ ；

j＝1,…,N；G(i,j)＝GF(i,D _j ) (ii) a Writing the above formula into a matrix form to obtain:

note that [ a ] ₁ a ₂ … a _M a _M+1 ] ^T Describing each son as a trend relation vectorA trend of the relationship between the array measurements. In the Born approximation, the equation can approximate subarray measurements of any radial one-dimensional formation. The left term of the chromatography matrix equation is the measurement signal of each subarray, the unknown quantity of the chromatography matrix equation is the conductivity difference value between radially adjacent layers, the coefficient array of the chromatography matrix equation is composed of the borehole geometric factor, the radial geometric factor and the constant of each subarray, and the constraint condition of the chromatography matrix equation is determined according to the borehole environment parameters.

Optionally, to analyze the borehole effects of the array induction measurements, a borehole geometry factor G may be introduced _bh And is used for expressing the weight of the contribution of the borehole medium to the measurement signal. Considering the borehole as a radius

A cylinder with z direction from-infinity to + ∞, then:

introducing a radial integral geometric factor G _r Expressing the height to infinity (z direction from-infinity to + ∞) and the outer radius to the weight of the contribution of the r cylinder to the measurement signal, then:

as can be seen from the above, it is shown that,

borehole geometry factor G, same as borehole effects _bh Is composed of four borehole environment parameters (ecc, cal, sigma) _m ，σ _t ) And (6) determining. After the geometric factors are introduced, the logging response in some special cases is expressed explicitly. For example: in an infinite homogeneous formation model without accounting for invasion, in well conditions, the apparent conductivity measured by the ith sub-array can be expressed as:

σ _a (i)＝σ _m G _bh (i)+σ _t [1-G _bh (i)]。

as an alternative embodiment, the existing array induction borehole correction method is based on a radial two-layer model schematic diagram as shown in fig. 6, "borehole + infinite formation", which may be referred to as a (radial) two-layer model. At this time, according to the geometric factor theory, the basic logging response relationship of each subarray is as follows:

σ _a (i)＝σ _m G _bh (i)+σ _t [1-G _bh (i)]

wherein i represents a sub-array number, from 1 to M, the larger i is, the larger the probing depth of the sub-array is; g _bh (i) A borehole geometry factor for sub-array # i; sigma _a (i) Measuring the apparent conductivity of the obtained i-number sub-array; sigma _t Is the undisturbed formation conductivity.

Alternatively, in a two-layer model, the borehole correction is simply "replacing" the borehole fluid with a conductivity σ _t The corresponding borehole correction equation is:

σ _bhc (i)＝σ _a (i)-(σ _m -σ _t )G _bh (i)。

wherein σ _bhc (i) Corrected measurements for the borehole of sub-array # i.

Optionally, the array induction logging instrument describes an invasion profile, i.e., radial variation of resistivity or conductivity, by measuring multiple resistivity curves of detection depths (10 in-, 20in-, 30in-, 60in-, 90in-, or 120 in-) from shallow to deep in the radial direction, so as to accurately obtain the true resistivity Rt of the formation. The stratum model is an invasion model and is simplified into a (radial) three-layer model, such as a radial three-layer model schematic diagram shown in FIG. 7, wherein D _xo The invaded zone diameter. According to the geometry factor theory, the measured values of the subarray can then be expressed as:

σ _a (i)＝σ _m G _bh (i)+σ _xo [G _xo (i)-G _bh (i)]+σ _t [1-G _xo (i)]。

wherein σ _xo Conductivity of the formation for invaded zone; g _xo (i) The radial integral geometry factor of the invaded zone for the # i sub-array, which is a function of the radius of the invaded zone. "Replacing" wellbore fluids as having conductivity σ _xo The well hole correction is realized by invading the stratum with zone, and the corresponding well hole correction is realizedThe formula is as follows:

σ _bhc (i)＝σ _a (i)-(σ _m -σ _xo )G _bh (i)。

optionally, when no intrusion is present, σ _xo ＝σ _t The three-layer model can be transformed into a two-layer model, and the borehole correction formula is correspondingly transformed.

As an alternative embodiment, σ when there is drag reduction intrusion _m >σ _xo >σ _t . Due to the presence of the borehole and invasion, the raw measurement data is related to: sigma _a (1)>…>σ _a (i)>…>σ _a (M), σ as shown in FIG. 8 _m >σ _xo >σ _t And (4) the visual conductivity diagram after model correction. The results of borehole corrections based on the formula of the tri-layer model (taking invasion into account) in the figure are ideal results. When the borehole correction process is performed using a formula based on a two-layer model (without taking invasion into account), the method is based on the equation _t Lower than sigma _xo ，(σ _m -σ _t )>(σ _m -σ _xo ) Overcorrection occurs. This overcorrection amount gradually decreases as the sub-array number increases. Resulting in a reduced difference between subarray curves, resulting in a smaller variation in the intrusion of the resulting curves.

As an alternative embodiment, σ is when there is an increased resistance intrusion _xo <σ _t If the borehole correction process is performed using a formula based on a two-layer model (without taking invasion into account), the result follows σ _m 、σ _t And σ _xo There are relationship changes between them, and the following three cases occur: the first method comprises the following steps: sigma _t >σ _xo ≥σ _m . At this time, (σ) _m -σ _xo ) Is negative, (sigma) _m -σ _t ) Is also negative, and | σ _m -σ _t |>|σ _m -σ _xo I, and therefore overcorrection also occurs, resulting in a smaller intrusion difference in the final curve, σ as shown in fig. 9 _t >σ _xo ≥σ _m And (4) the visual conductivity diagram after model correction. And the second method comprises the following steps: sigma _t >σ _m >σ _xo . At this time (sigma) _m -σ _xo ) Is positive, but because (σ) _m -σ _t ) Negative, the borehole correction direction is the opposite direction, making the invasion difference of the final curve smaller, shown as σ in FIG. 10 _t >σ _m >σ _xo And (3) in time, the visual conductivity diagram after model correction. And the third is that: sigma _m ≥σ _t >σ _xo At this time (σ) _m -σ _t ) And (σ) _m -σ _xo ) Are all positive or zero, since (σ) _m -σ _xo )>(σ _m -σ _t ) Therefore, the borehole correction is not sufficient, making the invasion difference of the final curve smaller, as shown by σ in FIG. 11 _m ≥σ _t >σ _xo And (3) in time, the visual conductivity diagram after model correction.

When a well calibration formula based on a two-layer model is used, the invasion difference of the calibration result curve is always small. In addition, the commonly used borehole correction formula:

it can be understood that: all formation "produced" conductivity "deduces" to the wellbore except for the wellbore; this "produced" conductivity of all formations is the result of the borehole correction. If the formation conforms to the two-layer model, the borehole corrected result is equal to σ _t (ii) a However, if there is intrusion, i.e., a three-layer model, then the result after correction of this equation is:

in an alternative embodiment, the determining the corrected measure of the target formation based on the measure expression and the borehole correction formula to complete borehole correction includes: processing the measurement signal by using the measurement expression and the borehole correction formula to determine the target formation conductivity and the target borehole mud conductivity of the target formation; and performing borehole correction on the target measurement value by using the target formation conductivity and the target borehole mud conductivity.

In the embodiment of the invention, based on a radial invasion model, an expression formula of the measurement values of the subarray and a borehole correction formula are derived; starting from a multi-order stratum model, applying born approximation to deduce a general formula for describing the internal quantitative relevance of the measured values of the array induction subarrays; solving parameters required by borehole correction according to the internal correlation among the array induction subarray measurement values; finally, a borehole correction formula based on the radial invasion model is applied to carry out borehole correction.

As an optional embodiment, in the solving process, parameters that need to be determined according to experiments in the program only have weights when the fitting residual is calculated, and the top 4 sub-arrays that are shallow to be detected need to be concerned, so the weights corresponding to the sub-arrays are selected to be 1.0, and the weights of the other sub-arrays are gradually reduced. Four alternative borehole correction modes are set up in the program and function as follows:

TABLE 1

Borehole correction mode	Description of functions
		General rule	Knowing sigma m Cal and Ecc, adaptive solving for σ xo And σ t
Self-adaptive slurry obtaining method	Knowing Cal and Ecc, the solution to σ is adaptive m 、σ xo And σ t
		Self-adaptive well diameter calculation	Knowing sigma m And Ecc, and adaptively solving Cal and sigma xo And σ t
Adaptive solution to eccentricity	Known as σ m And Cal, solving Ecc and sigma adaptively xo And σ t

As an alternative embodiment, by applying the adaptive solver, the curve difference can reliably describe the invasion profile, and the true resistivity Rt of the formation can be obtained more accurately.

It should also be noted that the radial three-layer model is a great improvement over the two-layer model due to the inclusion of the invaded zone, but it is still far from the actual radial change in formation conductivity. The invasion transition path diagram of fig. 12, in which the dotted line part shows the radial section closer to the actual formation, represents only four possible paths of the transition zone from the borehole wall flushing zone to the undisturbed formation, and there are infinite possibilities of the actual path.

By adopting the method, the inline relation exists among the measuring curves based on the array induction subarray, the chromatography concept is applied, a plurality of step functions are introduced to approach any radial one-dimensional stratum, the mathematical expression of the inline relation is deduced according to the geometric factor theory, and the method is applied to self-adaptive borehole correction. Based on the invasion model, a new borehole correction formula is determined. The parameters required by borehole correction are correctly solved by applying the internal relation between the array induction subarray measurement curves and a new borehole correction formula, so that invasion influence is prevented from being regarded as borehole influence, curve difference is not reduced after borehole correction, and the capacity of curve radial analysis invasion is maintained.

Through the steps, a new self-adaptive borehole environment parameter solving algorithm can be constructed by applying the inline relation among the array induction subarray measurement curves, and the accuracy of a self-adaptive solving result is greatly improved; obtaining a target borehole correction formula based on the invasion model; the invasion influence is prevented from being regarded as the borehole influence to be corrected, the curve difference is not reduced after the borehole is corrected, and the invasion of the curve radial analysis is keptCapability. A schematic conductivity radial profile of a pre/post borehole correction array induction curve representation as shown in fig. 13, wherein prior to borehole correction, the borehole medium conductivity is generally significantly different from the conductivity of the formation (invaded wash zone, invaded transition zone, undisturbed formation), as shown in (a); the sensing of the array by the array at this time is the response of the model (a). Borehole correction is equivalent to "measuring conductivity as sigma _m By changing the borehole medium to a conductivity of sigma _xo The medium "of (a), when the array senses each subarray measurement is the response of the model of (b). The method solves the technical problems that the existing array induction borehole correction method has the defects that the model is not in accordance with the actual stratum and the internal correlation among the array induction subarray measurement values is low.

Example 2

According to an embodiment of the present invention, there is further provided an embodiment of an apparatus for implementing the array induction borehole correction method, fig. 14 is a schematic structural diagram of a borehole correction apparatus according to an embodiment of the present invention, as shown in fig. 14, the borehole correction apparatus includes: an acquisition module 140, a first determination module 142, a processing module 144, a second determination module 146, and a correction module 148, wherein:

an obtaining module 140, configured to obtain a measurement signal obtained by measuring each subarray of the array induction logging instrument, where the measurement signal at least includes: a conductivity curve group;

the first determining module 142 is configured to determine a wellbore environment parameter corresponding to a wellbore environment of a target formation, and obtain a target wellbore geometry factor of the wellbore environment from a wellbore geometry factor library according to the wellbore environment parameter by using a multidimensional spatial interpolation method, where the wellbore environment parameter includes: well bore cal, mud conductivity σ _m Formation conductivity σ _t And an eccentricity ecc, said wellbore geometry factor database being a pre-established wellbore geometry factor database, said wellbore geometry factor being a function of said wellbore environment parameter;

a processing module 144, configured to establish a chromatography matrix equation based on a radial formation model and the measurement signal, and solve the chromatography matrix equation to obtain a radial conductivity distribution of the target formation, where the radial formation model is used to characterize a conductivity change of the target formation;

a second determining module 146, configured to determine a measurement expression according to the radial formation model and the radial conductivity distribution, and determine a borehole correction formula by using the target borehole geometry factor corresponding to the borehole environment;

and a correction module 148, configured to determine a corrected measured value of the target formation based on the measured value expression and the borehole correction formula, so as to complete borehole correction.

It should be noted here that the acquiring module 140, the first determining module 142, the processing module 144, the second determining module 146 and the correcting module 148 correspond to steps S102 to S110 in embodiment 1, and the modules are the same as the corresponding steps in the implementation example and the application scenario, but are not limited to the disclosure in embodiment 1. It should be noted that the modules described above may be executed in a computer terminal as part of an apparatus.

It should be noted that, reference may be made to the relevant description in embodiment 1 for optional or preferred embodiments of this embodiment, and details are not described here again.

The borehole correction device may further include a processor and a memory, and the acquiring module 140, the first determining module 142, the processing module 144, the second determining module 146, the correcting module 148, and the like are stored in the memory as program units, and the processor executes the program units stored in the memory to implement corresponding functions.

The processor comprises a kernel, and the kernel calls a corresponding program unit from the memory, wherein one or more than one kernel can be arranged. The memory may include volatile memory in a computer readable medium, random Access Memory (RAM) and/or nonvolatile memory such as Read Only Memory (ROM) or flash memory (flash RAM), and the memory includes at least one memory chip.

According to an embodiment of the present application, there is also provided an embodiment of a non-volatile storage medium. Optionally, in this embodiment, the non-volatile storage medium includes a stored program, and the apparatus in which the non-volatile storage medium is located is controlled to execute any one of the array induction borehole correction methods when the program is executed.

Optionally, in this embodiment, the nonvolatile storage medium may be located in any one of a group of computer terminals in a computer network, or in any one of a group of mobile terminals, and the nonvolatile storage medium includes a stored program.

Optionally, the apparatus in which the non-volatile storage medium is controlled to perform the following functions when the program is executed: obtaining a measurement signal obtained by measuring each subarray of the array induction logging instrument, wherein the measurement signal at least comprises: a conductivity curve group; determining a borehole environment parameter corresponding to a borehole environment of a target stratum, and acquiring a target borehole geometric factor of the borehole environment from a borehole geometric factor library according to the borehole environment parameter by adopting a multi-dimensional spatial interpolation method, wherein the borehole environment parameter comprises: well bore cal, mud conductivity σ _m Formation conductivity σ _t And an eccentricity ecc, said wellbore geometry factor database being a pre-established wellbore geometry factor database, said wellbore geometry factor being a function of said wellbore environment parameter; establishing a chromatography matrix equation based on a radial stratum model and the measurement signal, and solving the chromatography matrix equation to obtain radial conductivity distribution of the target stratum, wherein the radial stratum model is used for representing conductivity change of the target stratum; determining a measurement value expression according to the radial stratum model and the radial conductivity distribution, and determining a borehole correction formula by adopting the target borehole geometric factor corresponding to the borehole environment; and determining the corrected measured value of the target stratum based on the measured value expression and the borehole correction formula, and completing borehole correction.

Optionally, the apparatus in which the non-volatile storage medium is controlled to perform the following functions when the program is executed: determining a first variation range of the well diameter, and determining a first preset number of discrete points based on the first variation range to form a first-dimension discrete point set; determining a second variation range of the conductivity of the slurry, and determining a second preset number of discrete points based on the second variation range to form a second dimension discrete point set; determining a third variation range of the formation conductivity, and determining a third preset number of discrete points based on the third variation range to form a third-dimension discrete point set; determining a fourth variation range of the eccentricity, and determining a fourth preset number of discrete points based on the fourth variation range to form a fourth-dimensional discrete point set; and performing combined calculation on the first-dimension discrete point set, the second-dimension discrete point set, the third-dimension discrete point set and the fourth-dimension discrete point set to obtain the borehole geometric factors of the sub-arrays, and storing the borehole geometric factors of all the sub-arrays according to a preset storage rule to construct the borehole geometric factor library.

Optionally, the apparatus in which the nonvolatile storage medium is controlled when the program is running performs the following functions: determining a value of a borehole environment parameter corresponding to the borehole environment of the target formation; acquiring the closest target discrete point from the borehole geometric factor library according to the value of the borehole environmental parameter; and determining the geometric factor of the target borehole by adopting the multi-dimensional space interpolation method based on the first-dimension discrete point set, the second-dimension discrete point set, the third-dimension discrete point set and the fourth-dimension discrete point set.

Optionally, the apparatus in which the non-volatile storage medium is controlled to perform the following functions when the program is executed: dividing the target stratum into a plurality of layers in the horizontal radius direction, wherein the borehole is the innermost layer, and the stratums except the borehole are an invasion flushing zone, an invasion transition zone and an undisturbed stratum; and constructing the radial formation model using the wellbore, the invaded wash zone, the invaded transition zone, and the undisturbed formation.

Optionally, the apparatus in which the non-volatile storage medium is controlled to perform the following functions when the program is executed: determining the borehole geometric factors of each subarray and the radial geometric factors of each subarray; establishing the chromatography matrix equation based on the radial formation model and the measurement signals, the borehole geometry factor and the radial geometry factor, wherein the left term of the chromatography matrix equation is the measurement signals of each subarray, the unknown quantity of the chromatography matrix equation is the conductivity difference between radially adjacent layers, the coefficient matrix of the chromatography matrix equation is composed of the borehole geometry factor, the radial geometry factor and a constant of each subarray, and the constraint condition of the chromatography matrix equation is determined according to the borehole environment parameters; and solving the chromatographic matrix equation to obtain the radial conductivity distribution of the target stratum.

Optionally, the apparatus in which the non-volatile storage medium is controlled to perform the following functions when the program is executed: processing the measurement signal by using the measurement expression and the borehole correction formula to determine the target formation conductivity and the target borehole mud conductivity of the target formation; and performing borehole correction on the target measurement value by using the target formation conductivity and the target borehole mud conductivity.

According to an embodiment of the present application, there is also provided an embodiment of a processor. Optionally, in this embodiment, the processor is configured to execute a program, where the program executes any one of the array induction borehole correction methods.

There is also provided, in accordance with an embodiment of the present application, an embodiment of an electronic device, including a memory and a processor, where the memory stores therein a computer program, and the processor is configured to execute the computer program to perform any one of the array-induced borehole correction methods described above.

There is also provided, in accordance with an embodiment of the present application, an embodiment of a computer program product, which, when executed on a data processing device, is adapted to execute a program for initializing the steps of the array induced borehole correction method of any of the above.

The above-mentioned serial numbers of the embodiments of the present invention are only for description, and do not represent the advantages and disadvantages of the embodiments.

In the above embodiments of the present invention, the description of each embodiment has its own emphasis, and reference may be made to the related description of other embodiments for parts that are not described in detail in a certain embodiment.

In the embodiments provided in the present application, it should be understood that the disclosed technical content can be implemented in other manners. The above-described apparatus embodiments are merely illustrative, and for example, the division of the units may be a logical division, and in actual implementation, there may be another division, for example, multiple units or components may be combined or may be integrated into another system, or some features may be omitted, or may not be executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, units or modules, and may be in an electrical or other form.

The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one position, or may be distributed on a plurality of units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.

In addition, functional units in the embodiments of the present invention may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit. The integrated unit may be implemented in the form of hardware, or may also be implemented in the form of a software functional unit.

The integrated unit, if implemented in the form of a software functional unit and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a removable hard disk, a magnetic disk, or an optical disk, and various media capable of storing program codes.

The foregoing is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, various modifications and decorations can be made without departing from the principle of the present invention, and these modifications and decorations should also be regarded as the protection scope of the present invention.

Claims (10)

1. An array induction borehole correction method, comprising:

obtaining a measurement signal obtained by measuring each subarray of the array induction logging instrument, wherein the measurement signal at least comprises: a conductivity curve group;

determining borehole environment parameters corresponding to the borehole environment of a target stratum, and acquiring target borehole geometric factors of the borehole environment from a borehole geometric factor library according to the borehole environment parameters by adopting a multidimensional spatial interpolation method, wherein the borehole environment parameters comprise: well bore cal, mud conductivity σ _m Formation conductivity σ _t And an eccentricity ecc, the wellbore geometry factor database being a pre-established wellbore geometry factor database, the wellbore geometry factor being a function of the wellbore environment parameter;

establishing a chromatography matrix equation based on a radial stratum model and the measurement signal, and solving the chromatography matrix equation to obtain radial conductivity distribution of the target stratum, wherein the radial stratum model is used for representing conductivity change of the target stratum;

determining a measurement value expression according to the radial stratum model and the radial conductivity distribution, and determining a borehole correction formula by adopting the target borehole geometric factor corresponding to the borehole environment;

and determining the corrected measured value of the target stratum based on the measured value expression and the borehole correction formula, and completing borehole correction.

2. The method of claim 1, wherein prior to obtaining the target wellbore geometry factor for the wellbore environment from a wellbore geometry factor library using multi-dimensional spatial interpolation based on the wellbore environment parameters, the method further comprises:

determining the variation range of the well diameter, and determining a preset number of discrete points of the well diameter based on the variation range of the well diameter to form a first-dimension discrete point set;

determining the variation range of the mud conductivity, and determining a preset number of discrete points of the mud conductivity based on the variation range of the mud conductivity to form a second-dimension discrete point set;

determining the variation range of the stratum conductivity, and determining a preset number of discrete points of the stratum conductivity based on the variation range of the stratum conductivity to form a third-dimension discrete point set;

determining the variation range of the eccentricity, and determining a preset number of discrete points of the eccentricity based on the variation range of the eccentricity to form a fourth-dimensional discrete point set;

and combining the first-dimension discrete point set, the second-dimension discrete point set, the third-dimension discrete point set and the fourth-dimension discrete point set, calculating to obtain the borehole geometric factors of each sub-array, and storing the borehole geometric factors of all the sub-arrays according to a preset storage rule to construct the borehole geometric factor library.

3. The method according to claim 2, wherein the determining of the borehole environment parameter corresponding to the borehole environment of the target formation and the obtaining of the target borehole geometry factor of the borehole environment from the borehole geometry factor library according to the borehole environment parameter by using the multidimensional spatial interpolation method comprises:

determining a value of a wellbore environment parameter corresponding to the wellbore environment of the target formation;

acquiring the closest target discrete point from the borehole geometry factor library according to the value of the borehole environment parameter;

and determining the target borehole geometric factor based on the first dimension discrete point set, the second dimension discrete point set, the third dimension discrete point set and the fourth dimension discrete point set by adopting the multi-dimensional space interpolation method.

4. The method of claim 1, wherein prior to establishing a chromatography matrix equation based on a radial formation model and the measurement signal and solving the chromatography matrix equation to obtain a radial conductivity distribution of the target formation, the method further comprises:

dividing the target stratum into a plurality of layers in the radius direction, wherein a borehole is the innermost layer, and the stratums except the borehole are an invasion flushing zone, an invasion transition zone and an undisturbed stratum;

constructing the radial formation model using the wellbore, the invaded wash zone, the invaded transition zone, and the undisturbed formation.

5. The method of claim 1, wherein establishing a tomographic matrix equation based on a radial formation model and the measurement signal, and solving the tomographic matrix equation to obtain a radial conductivity distribution of the target formation comprises:

determining a borehole geometry factor for each of the sub-arrays and a radial geometry factor for each of the sub-arrays;

establishing the chromatography matrix equation based on the radial stratum model and the measurement signal, the borehole geometry factor and the radial geometry factor, wherein the left term of the chromatography matrix equation is the measurement signal of each subarray, the unknown quantity of the chromatography matrix equation is the conductivity difference value between radially adjacent layers, a coefficient array of the chromatography matrix equation is composed of the borehole geometry factor, the radial geometry factor and a constant of each subarray, and the constraint condition of the chromatography matrix equation is determined according to the borehole environment parameter;

and solving the chromatographic matrix equation to obtain the radial conductivity distribution of the target stratum.

6. The method according to any one of claims 1 to 5, wherein the determining the corrected measure of the target formation based on the measure expression and the borehole correction formula to perform borehole correction comprises:

processing the measurement signal by using the measurement expression and the borehole correction formula to determine target formation conductivity and target borehole mud conductivity of the target formation;

and performing borehole correction on the target measurement value by using the target formation conductivity and the target borehole mud conductivity.

7. An array induction borehole correction device, comprising:

the acquisition module is used for acquiring measurement signals obtained by measuring each subarray of the array induction logging instrument, wherein the measurement signals at least comprise: a conductivity curve group;

the first determination module is used for determining borehole environment parameters corresponding to the borehole environment of the target stratum and acquiring target borehole geometric factors of the borehole environment from a borehole geometric factor library according to the borehole environment parameters by adopting a multi-dimensional spatial interpolation method, wherein the borehole environment parameters comprise: well bore cal, mud conductivity σ _m Formation conductivity σ _t And an eccentricity ecc, the wellbore geometry factor database being a pre-established wellbore geometry factor database, the wellbore geometry factor being a function of the wellbore environment parameter;

the processing module is used for establishing a chromatography matrix equation based on a radial stratum model and the measuring signal, and solving the chromatography matrix equation to obtain radial conductivity distribution of the target stratum, wherein the radial stratum model is used for representing conductivity change of the target stratum;

the second determination module is used for determining a measurement value expression according to the radial stratum model and the radial conductivity distribution and determining a borehole correction formula by adopting the target borehole geometric factor corresponding to the borehole environment;

and the correction module is used for determining the corrected measured value of the target stratum based on the measured value expression and the borehole correction formula so as to finish borehole correction.

8. A non-volatile storage medium storing a plurality of instructions adapted to be loaded by a processor and to perform the array induced borehole correction method of any one of claims 1 to 6 and a borehole geometry factor database.

9. A processor for executing a program, wherein the program is arranged to execute the array induced borehole correction method of any of claims 1 to 6 when executed.

10. An electronic device comprising a memory in which is stored a computer program and a borehole geometry factor database, and a processor arranged to run the computer program to perform the array induced borehole correction method of any of claims 1 to 6.

Images