CN117662107A

CN117662107A - Borehole correction method for array induction instrument and related equipment

Info

Publication number: CN117662107A
Application number: CN202211043415.8A
Authority: CN
Inventors: 赵彦伟; 梁小兵; 贾占军; 谢昱北; 刘越; 陈文�; 卢炳文; 孙学凯; 朱瑞明; 杜有定; 史广岩; 张国艳; 伍莹
Original assignee: China National Petroleum Corp; China Petroleum Logging Co Ltd
Current assignee: China National Petroleum Corp; China Petroleum Logging Co Ltd
Priority date: 2022-08-29
Filing date: 2022-08-29
Publication date: 2024-03-08

Abstract

The application discloses a borehole correction method for an array induction instrument and related equipment. The method comprises the following steps: performing resolution matching based on the original conductivity curve to obtain multi-frequency multi-subarray stratum apparent conductivity with the same resolution; acquiring an initial value of a borehole correction variable, a borehole geometry factor, an invasion geometry factor and an undisturbed stratum geometry factor coefficient; determining a first target parameter by using a linear optimization algorithm or a method for solving a linear equation set based on an induction measurement response model of the wellbore geometry factor and the invasion geometry factor; constructing an objective function of which the actual logging value is matched with the simulated measurement response value of the induction measurement response model, and inverting and determining a second objective parameter by a nonlinear optimization method according to the first objective parameter; and obtaining a target borehole geometry factor, a target flushing zone geometry factor and a target undisturbed stratum geometry factor according to the second target parameter table lookup interpolation so as to realize borehole correction of the array induction instrument.

Description

Borehole correction method for array induction instrument and related equipment

Technical Field

The present disclosure relates to the field of oil and gas detection, and more particularly, to a wellbore calibration method for an array induction instrument and related equipment.

Background

Geophysical well logging is an important means for exploring and developing mineral resources such as oil gas, metal and the like in a well bore for drilling a stratum by utilizing physical means such as electromagnetic waves, sound waves, radioactivity, nuclear magnetic resonance, optical fibers and the like to distinguish electromagnetic parameters, elastic mechanical parameters, element properties and the like of fluid media in stratum rock and rock pores.

In petroleum logging hydrocarbon evaluation, formation resistivity is an important parameter in estimating hydrocarbon reserves. Array induction logging based on electromagnetic induction principle is one of the most effective logging methods at present. The method measures the formation resistivity information which is abundant in the well, eliminates the environmental influence of the well and the like by using the treatment of the well correction and the like, eliminates the environmental influence of the surrounding rock and the like by using the treatment of the software focusing and the like, and obtains the formation resistivity and the like in the original state by using the treatment of the radial inversion and the like, and has the advantages of high resolution, deep detection depth, obvious intrusion indication and the like.

The traditional borehole correction method is to design an objective function based on four borehole environment variables of borehole diameter, instrument eccentric amount, slurry conductivity and stratum conductivity and measurement vision conductivity matching mode research through skin-seeking correction and resolution matching, obtain borehole environment variables through optimization, calculate borehole signals and eliminate the influence of the borehole to the greatest extent. The wellbore correction method lays a foundation for obtaining the true resistivity of the undisturbed stratum, describing the invasion profile, determining the stratum permeability and the fluid type, and searching and evaluating the potential hydrocarbon reservoir. However, since conventional borehole correction methods treat formations other than the borehole as uniform formations, neither calculation of borehole correction pattern libraries nor adaptive optimization of borehole environment variables take into account the effects of invasion, nor the differences in resistivity in the radial direction.

Disclosure of Invention

In the summary, a series of concepts in a simplified form are introduced, which will be further described in detail in the detailed description. The summary of the invention is not intended to define the key features and essential features of the claimed subject matter, nor is it intended to be used as an aid in determining the scope of the claimed subject matter.

In order to achieve the accuracy of exploratory well logging, in a first aspect, the present invention provides a method for calibrating a borehole of an array induction apparatus, the method comprising:

performing resolution matching based on the original conductivity curve to obtain multi-frequency multi-subarray stratum apparent conductivity with the same resolution;

obtaining an initial value, a borehole geometry factor, an invasion geometry factor and an undisturbed formation geometry factor coefficient of a borehole correction variable, wherein the borehole correction variable comprises a borehole diameter, a slurry conductivity, an undisturbed formation conductivity, an instrument eccentricity, an invasion radius and a flushing charge conductivity;

determining a first target parameter by using a linear optimization algorithm or a method for solving a linear equation set based on an induction measurement response model of a wellbore geometry factor and an invasion geometry factor, wherein the first target parameter is a first sub-optimal solution of mud conductivity, flushing live conductivity and undisturbed formation conductivity, the induction measurement response model is a quasi-linear combination of the measurement apparent conductivity, which is mud conductivity, flushing live conductivity and undisturbed formation conductivity, and a combination coefficient of the quasi-linear combination is the wellbore geometry factor, flushing zone geometry factor and undisturbed formation geometry factor;

Constructing an objective function of which the actual logging value is matched with the simulated measuring response value of the induction measuring response model, and inverting and determining a second objective parameter by a nonlinear optimization method according to the first objective parameter, wherein the second objective parameter is a second sub-optimal solution of the diameter of a well bore, the slurry conductivity, the conductivity of an undisturbed stratum, the eccentricity of an instrument, the invasion radius and the flushing charge conductivity;

and obtaining a target borehole geometry factor, a target flushing zone geometry factor and a target undisturbed stratum geometry factor according to the second target parameter table lookup interpolation so as to realize borehole correction of the array induction instrument.

Optionally, the performing resolution matching based on the original conductivity curve to obtain the apparent conductivity of the multi-frequency multi-subarray stratum with the same resolution includes:

acquiring an original conductivity curve of each subarray under the full frequency;

acquiring resolution difference information of original conductivity curves between adjacent arrays;

and matching a lower resolution curve in the original conductivity curves between adjacent arrays based on the resolution difference information to obtain resolution matching apparent conductivity.

Optionally, the obtaining the initial value of the borehole correction variable, the borehole geometry factor, the invasion geometry factor and the undisturbed formation geometry factor coefficient includes:

Determining an initial value of a borehole correction variable according to the apparent conductivity, instrument geometric feature information, borehole geometric features, borehole diameters and slurry conductivity acquisition conditions and software processing parameter selection logic;

constructing a borehole geometry factor library and an invasion geometry factor library according to borehole correction variables;

and performing multi-slice interpolation calculation based on the borehole geometry factor library and the invasion geometry factor library to obtain the borehole geometry factor, the invasion geometry factor and the undisturbed stratum geometry factor.

Optionally, the initial values of the borehole correction variables include an initial value ecc0 of an instrument eccentricity, the instrument geometry information includes an instrument outer diameter dtool and a centralizer size standoff, and the borehole geometry includes a borehole diameter cal;

determining initial values of borehole correction variables based on the apparent conductivity, instrument geometry information, borehole geometry, and software process parameter selection logic, comprising:

the initial value of the instrument eccentricity, ecc0, is obtained according to the following:

where dtool is the instrument outer diameter, stardoff is the centralizer size, and cal is the wellbore diameter.

Optionally, the initial values of the borehole correction variables include initial values sigma of the flushing charge conductivity _xo0 ；

obtaining initial value sigma of flushing electrification conductivity according to the following formula _xo0 ：

In the method, in the process of the invention,for the weight coefficient corresponding to the apparent resistivity of the corresponding subarray and the corresponding frequency, +.>For apparent conductivity of different subarrays under different frequencies, M is the number of subarrays, and N is the number of measuring frequencies;

and/or the number of the groups of groups,

the initial values of the borehole correction variables include initial values sigma of the conductivity of the undisturbed formation _t0 ；

obtaining an initial value sigma of the conductivity of the undisturbed stratum according to the following formula _t0 ：

In the method, in the process of the invention,for the weight coefficient corresponding to the apparent resistivity of the corresponding subarray and the corresponding frequency, +.>For apparent conductivity of different subarrays at different frequencies, M is the number of subarrays, and N is the number of measuring frequencies.

Optionally, constructing an objective function of matching the actual logging value with the simulated measurement response value of the induction measurement response model, and inverting and determining a second objective parameter by using a nonlinear optimization method according to the first objective parameter, including:

Respectively changing one or more of the borehole diameter, the instrument eccentricity, the slurry conductivity and the undisturbed stratum conductivity, the flushing charge conductivity and the flushing band radius 6 borehole correction parameters, and solving the 6 borehole correction parameters by inversion based on an optimization method;

optimizing a linear equation set constructed based on a geometric factor theoretical model and the optimization formula to obtain the second target parameter;

wherein, the optimization formula is as follows:

wherein O (cal, ecc, sigma _m ,σ _xo ,l _xo ,σ _t ) To optimize the objective function b _jf For the weight coefficient of each frequency of each subarray of the target induction instrument, Q is a punishment function or a regularization term, lambda is a regularization parameter, P is a norm, cal is the borehole diameter, ecc is instrument eccentricity, sigma _m For mud conductivity, sigma _xo To flush the charged conductivity, l _xo To flush the radius sigma _t For the conductivity of the undisturbed stratum, M is the number of subarrays, N is the number of measuring frequencies,2ft resolution matching for the jth sub-array f frequency is based on the conductivity measurement processing result,/for the jth sub-array f frequency>The model response value of the corresponding frequency of the corresponding subarray is calculated according to the geometric factors in a pre-simulation mode;

the above linear system of equations is:

andThe initial values of the borehole geometry factor, the flush zone geometry factor, the undisturbed formation geometry factor, and the intermediate values of each iteration step are represented for the jth sub-array, f, frequency, respectively.

Optionally, the second target parameter includes a target borehole diameter cal _inv Eccentricity ecc of target instrument _inv Target mud conductivity sigma _minv Target flush charge conductivity sigma _xoinv Target flush zone radius l _xoinv Target undisturbed formation conductivity sigma _tinv ；

The obtaining the target borehole geometry factor, the target flushing zone geometry factor and the target undisturbed formation geometry factor according to the second target parameter table lookup interpolation to realize borehole correction of the array induction instrument comprises:

the array induction instrument was calibrated using the following:

in the method, in the process of the invention,borehole corrected rearview conductivity value, cal, representing the jth sub-array f frequency _inv To target borehole diameter, ecc _inv For target instrument eccentricity, σ _minv For target mud conductivity, σ _xoinv To flush the target with electrical conductivity, l _xoinv For target irrigation band radius, sigma _tinv Is the target undisturbed formation conductivity.

In a second aspect, the present invention also provides an array induction instrument wellbore correction device, comprising:

a first acquisition unit for performing resolution matching based on the original conductivity curve to acquire multi-frequency multi-subarray formation apparent conductivity having the same resolution;

a second obtaining unit, configured to obtain an initial value of a borehole correction variable, a borehole geometry factor, an invasion geometry factor, and an undisturbed formation geometry factor coefficient, where the borehole correction variable includes a borehole diameter, a mud conductivity, an undisturbed formation conductivity, an instrument eccentricity, an invasion radius, and a flushing band conductivity;

The first determining unit is used for determining a first target parameter by utilizing a linear optimization algorithm or a method for solving a linear equation set based on an induction measurement response model of a borehole geometry factor and an invasion geometry factor, wherein the first target parameter is a first sub-optimal solution of mud conductivity, flushing charged conductivity and undisturbed formation conductivity, the induction measurement response model is a quasi-linear combination of the measurement apparent conductivity, namely mud conductivity, flushing charged conductivity and undisturbed formation conductivity, and the combination coefficient of the quasi-linear combination is the borehole geometry factor, flushing charged geometry factor and undisturbed formation geometry factor;

the second determining unit is used for constructing an objective function of which the actual logging value is matched with the simulated measurement response value of the induction measurement response model, and inverting and determining a second objective parameter by a nonlinear optimization method according to the first objective parameter, wherein the second objective parameter is a second sub-optimal solution of the borehole diameter, the slurry conductivity, the undisturbed stratum conductivity, the instrument eccentricity, the invasion radius and the flushing charge conductivity;

and the correction unit is used for looking up a table and interpolating to obtain a target borehole geometry factor, a target flushing zone geometry factor and a target undisturbed stratum geometry factor according to the second target parameter so as to realize borehole correction of the array induction instrument.

In a third aspect, an electronic device, comprising: a memory, a processor and a computer program stored in and executable on the processor for performing the steps of the method of borehole correction for an array induction apparatus as described in any one of the first aspects above when the computer program stored in the memory is executed.

In a fourth aspect, the present invention also provides a computer readable storage medium having stored thereon a computer program which when executed by a processor implements the method of borehole correction for an array induction apparatus of any of the first aspects.

In summary, the method provided by the embodiment of the application introduces a new array induction logging simulation calculation technology and a nonlinear optimization technology into the field, so that the establishment of a borehole correction model and a borehole geometry factor library, which are more approximate to the actual situation, and model solving can be realized on a computer by numerical simulation larger than the original model. The traditional non-invasive borehole environment correction model has reasonable correction results on mud rock sections with no invasion or small invasion degree, but when a sandstone target reservoir with obvious invasion exists, particularly when the contrast between the mud conductivity and the stratum conductivity is high and the borehole is large, invasion influences are often regarded as invasion profile errors caused by overcorrection or undercorrection for correcting or partially eliminating borehole influences, and incorrect flushing zone resistivity and stratum true resistivity results obtained by subsequent one-dimensional resistivity inversion. The method provided by the application considers invasion from the establishment of a theoretical model, and fundamentally overcomes the defects of the traditional borehole correction method. Meanwhile, skin-seeking correction is not performed on the multi-frequency measured array induction well data in advance in order to reduce the dimension of a geometric factor library and the time complexity of an algorithm, so that richer well bore, invasion and undisturbed stratum conductivity information is reserved to a certain extent, conditions are created for the inversion process of environment correction with more parameters of an invasion model, and the optimization problem solving of the well bore correction environment variables is more stable and more accurate. Additional advantages, objects, and features of the invention will be set forth in part in the description which follows and in part will become apparent to those having ordinary skill in the art upon examination of the following or may be learned from practice of the invention.

Drawings

Various other advantages and benefits will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiments. The drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the specification. Also, like reference numerals are used to designate like parts throughout the figures. In the drawings:

FIG. 1 is a schematic flow chart of a method for calibrating a borehole of an array induction instrument according to an embodiment of the present application;

FIG. 2 is a schematic diagram of a four-parameter calibration model according to an embodiment of the present disclosure;

FIG. 3 is a schematic diagram of a three-parameter intrusion model according to an embodiment of the present disclosure;

FIG. 4 is a schematic diagram of another three-parameter intrusion model provided in an embodiment of the present application;

FIG. 5 is a schematic diagram of a six-parameter calibration model according to an embodiment of the present disclosure;

FIG. 6 is a schematic flow chart of another method for calibrating a borehole of an array induction apparatus according to an embodiment of the present application;

FIG. 7 is a schematic diagram of a resistivity matching curve according to an embodiment of the present disclosure;

FIG. 8 is a schematic diagram of another resistivity matching curve provided in an embodiment of the present application;

FIG. 9 is a schematic diagram of a high mud conductivity and formation conductivity provided in an embodiment of the present application;

FIG. 10 is a schematic diagram of another high mud conductivity and formation conductivity provided in an embodiment of the present application;

FIG. 11 is a schematic structural diagram of a borehole correction device for an array induction apparatus according to an embodiment of the present disclosure;

fig. 12 is a schematic structural diagram of an electronic device for a borehole calibration method of an array induction apparatus according to an embodiment of the present application.

Detailed Description

The terms "first," "second," "third," "fourth" and the like in the description and in the claims of this application and in the above-described figures, if any, are used for distinguishing between similar objects and not necessarily for describing a particular sequential or chronological order. It is to be understood that the data so used may be interchanged where appropriate such that the embodiments described herein may be implemented in other sequences than those illustrated or otherwise 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. The following description of the embodiments of the present application will be made clearly and fully with reference to the accompanying drawings, in which it is evident that the embodiments described are only some, but not all, of the embodiments of the present application.

As shown in FIG. 2, the four-parameter empirical correction model currently employed assumes a conductivity of σ _t In a uniform formation of (1), the slurry conductivity is sigma under well conditions _m The diameter of the well is cal, the outer diameter of the instrument is d _tool The distance of the instrument from the well wall is x, and the relative eccentricity of the instrument is defined as ecc as shown in formula (1):

obviously, when the instrument is centered, the relative eccentricity is 0; when the instrument is fully stuck to the well wall, the relative eccentricity is 1. As shown in fig. 2, apparent conductivity measured by the ith sub-array after skin-corrected according to induction logging theory in a uniform formation model under logging conditions without taking into account invasionCan be approximated as formula (2):

wherein,the borehole geometry factor is a weight representing the contribution of the borehole medium to the measured signal. The wellbore geometry factor may be approximated as instrument relative eccentricity ecc, borehole diameter cal, and mud conductivity sigma _mud And formation conductivity sigma _t Is a quaternary function of (a). Apparent conductivity->Can be regarded as mud conductivity sigma _mud And formation conductivity sigma _t Is a quasi-linear function of (c). Relative eccentricity ecc, borehole diameter cal and mud conductivity sigma of instrument _mud And formation conductivity sigma _t Four variables are called wellsAn eye correction environment variable, wherein instrument relative eccentricity ecc and borehole diameter cal are referred to as geometric variables; mud conductivity sigma _mud And formation conductivity sigma _t Referred to as medium conductivity variables.

In order to improve logging accuracy, the application provides a 6-parameter-based borehole correction measurement model, which also considers flushing electrification conductivity sigma on the basis of considering the four variables _xo And a flushing belt radius l _xo Two variables. Referring to fig. 3, fig. 4 is a schematic diagram of a formation model for array induction borehole correction, the formation model comprising a borehole (where the medium is a dielectric having a conductivity σ _mud Commonly known as mud), a flushing zone (in which the medium is mud filtrate to fully displace the formation after the original fluid in the rock in the zone, the zone conductivity being sigma _xo ) Undisturbed formation (region not affected by drilling and invasion of mud filtrate in the wellbore, region conductivity is σ _t ) Three regions. Wherein cal and cal/2 denote the diameter of the borehole and the borehole radius, respectively, l _xo Indicating the irrigation band radius. FIG. 5 is a schematic diagram of an array induction measurement response based on geometry factor theory. The horizontal axis of the graph represents the radial (vertical wellbore direction) distance from the wellbore centerline, the upper half is the geometric factor, the ordinate is the dimensionless quantity, the area size below the curve represents the contribution weight of the formation of the region to the induction measurement signal, and the lower half is the formation conductivity, the ordinate is the conductivity. The induction measurement response is a weighted average or linear combination of conductivities of different regions of the formation (including the wellbore), where the weighting or combining coefficients are geometric factors.

Referring to fig. 1, a schematic flow chart of a method for calibrating a borehole of an array induction apparatus according to an embodiment of the present application may specifically include:

s110, performing resolution matching based on an original conductivity curve to obtain multi-frequency multi-subarray stratum apparent conductivity with the same resolution;

illustratively, the original conductivity curve is obtained by the array induction instrument during the single-frequency (or multi-frequency) well logging, and the apparent conductivity is directly subjected to resolution matching, so that M multiplied by N resolution matching apparent conductivities can be obtainedWherein M is the number of subarrays in the array induction instrument, and N is the number of measuring frequencies. The array induction logging instrument is composed of a series of three-coil series subarrays with different source distances, wherein the subarrays are shallow in detection depth and high in longitudinal resolution, and the long subarrays are deep in detection depth and low in longitudinal resolution. The borehole correction process requires that the logs of the different subarrays have the same longitudinal resolution.

S120, acquiring an initial value of a borehole correction variable, a borehole geometry factor, an invasion geometry factor and an undisturbed stratum geometry factor coefficient, wherein the borehole correction variable comprises a borehole diameter, slurry conductivity, undisturbed stratum conductivity, instrument eccentricity, an invasion radius and flushing charge conductivity;

Exemplary borehole correction variables cal, ecc, σ _m 、σ _xo 、l _xo 、σ _t The difficulty level determined according to the initial value can be classified into two types. Wherein cal, ecc, sigma _m Are all variables in the wellbore, in particular the borehole diameter cal and the mud conductivity sigma _m Can be measured relatively accurately by other instruments that are downhole at the same time, while instrument eccentricity ecc, and initial values can be estimated relatively easily based on centralizer size standoff and well deviation conditions. Sigma (sigma) _xo 、l _xo 、σ _t Are both invasion and undisturbed formation related variables whose initial values can be estimated, and are variables that always need inversion solutions. Generally, the initial value cal of the well diameter ₀ The borehole correction control options may be selected as borehole diameter measurements, BS values (bit size used for drilling) or manually entered borehole sizes. Initial value sigma of slurry conductivity _m0 According to the borehole correction control option, the method can select the conductivity measurement value of the array induction instrument such as the self-charged conductivity measurement or three parameters (a logging instrument which simultaneously goes down well) or calculate the mud conductivity of each depth point according to the measurement value of the surface mud test box and the downhole temperature as an initial value.

The wellbore geometry factor, invasion geometry factor and undisturbed formation geometry factor coefficients can be obtained by interpolation of six-parameter wellbore correction measurement model parameters and a wellbore geometry factor library and an invasion geometry factor library. It should be noted that the six-parameter wellbore correction measurement model parameters and the wellbore geometry factor library and invasion geometry factor library settings may be set as follows:

Wellbore diameter cal (25 in): 4.0,4.5,5.0,5.5,6.0,6.5,7.0,7.5,8.0,8.5,9.0,9.5, 10.0, 10.5, 11.0, 12.0, 13.0, 14.0, 16.0, 18.0, 20.0, 22.0, 25.0, 25.5, 26.0

Mud conductivity sigma _m (20, units S/m): 0.001,0.01,0.1,0.2,0.5,0.8,1.0,2.0,3.0,4.0,6.0,8.0, 10.0, 13.0, 16.0, 20.0, 30.0, 40.0, 50.0, 100.0

Conductivity sigma of undisturbed stratum _t (22, units S/m): 0.001,0.01,0.02,0.05,0.08,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0,1.5,2.0,3.0,4.0,5.0,8.0, 10.0

Eccentricity ecc (21, dimensionless): 0.00,0.05,0.10,0.15,0.20,0.25,0.30,0.35,0.40,0.45,0.50,0.55,0.60,0.65,0.70,0.75,0.80,0.85,0.90,0.95,0.99

Radius of intrusion l _xo (10, units in): 4,6,8, 10, 15, 20, 30, 40, 60, 90

Flushing electrification conductivity sigma _xo (16, units S/m): 0.005,0.01,0.02,0.025,0.03,0.05,0.08,0.1,0.2,0.3,0.4,0.5,1.0,1.5,2.0,5.0.

S130, determining a first target parameter by utilizing a linear optimization algorithm or a method for solving a linear equation set based on an induction measurement response model of a borehole geometry factor and an invasion geometry factor, wherein the first target parameter is a first sub-optimal solution of slurry conductivity, flushing charge conductivity and undisturbed formation conductivity, the induction measurement response model is a quasi-linear combination of the measurement apparent conductivity, namely the slurry conductivity, the flushing charge conductivity and the undisturbed formation conductivity, and the combination coefficient of the quasi-linear combination is the borehole geometry factor, the flushing band geometry factor and the undisturbed formation geometry factor;

Exemplary, for an array induction instrument having M subarrays of N measurement frequencies, the measurements of each subarray of each frequency are based on the geometric factor theory in the presence of a borehole and an invasionCan be expressed as:

wherein,a value of the conductivity representing the f-th frequency measurement of the j-th subarray, a>And->The pre-calculated 6-variable borehole geometry factor library and the intrusion geometry factor library, namely the borehole geometry factor and the flushing zone geometry factor, respectively +.>Is the geometric factor coefficient of the undisturbed stratum, sigma _m 、σ _xo 、σ _t Respectively mud conductivity, flushing charge conductivity, undisturbed formation conductivity, cal, ecc, l _xo The wellbore diameter, instrument eccentricity, and flush zone radius, respectively. Equation (3) can be developed in the form of a quasi-linear combination relation (4) written as follows, namely an inductive measurement response model:

s140, constructing an objective function of which the actual logging value is matched with the simulated measurement response value of the induction measurement response model, and inverting and determining a second objective parameter by a nonlinear optimization method according to the first objective parameter, wherein the second objective parameter is a second sub-optimal solution of the diameter of a well bore, the conductivity of slurry, the conductivity of undisturbed stratum, the eccentricity of an instrument, the invasion radius and the conductivity of flushing charge;

The solving process of the nonlinear equation set (4) may be converted into a nonlinear optimization problem through nonlinear optimization until a preset convergence condition or a suspension condition is met to meet the accuracy requirement, so as to determine the second target parameter.

S150, obtaining a target borehole geometry factor, a target flushing zone geometry factor and a target undisturbed stratum geometry factor according to the second target parameter table lookup interpolation so as to realize borehole correction of the array induction instrument.

Illustratively, the array sensing instrument is calibrated using the following formula:

in the method, in the process of the invention,borehole corrected rearview conductivity value, cal, representing the jth sub-array f frequency _inv For a target borehole diameter, i.e. the borehole diameter is optimized or inverted, ecc _inv To optimize or invert the resulting instrument eccentricity, σ _minv To optimize or invert the resulting mud conductivity, σ _xoinv To optimize or invert the resulting flushing charge conductivity, l _xoinv To optimize or invert the resulting flush zone radius, sigma _tinv To optimize or invert the resulting undisturbed formation conductivity.

Since the geometry factor library table is constructed by calculating a few discrete points of the selected 6 borehole correction variables in advance, in order to improve the borehole correction accuracy, the borehole geometry factors corresponding to the six parameters need to be obtained by a multi-component piecewise interpolation method when determining the borehole geometry factors finally used for the borehole correction.

In summary, the method provided by the embodiment of the application introduces the array induction logging simulation calculation technology and the nonlinear optimization technology into the field, so that the establishment of a borehole correction model and a borehole geometric factor library, which are more approximate to the actual situation, and model solving can be realized on a computer by numerical simulation larger than the original standard model. The traditional non-invasive borehole environment correction model has reasonable correction results on mud rock sections with no invasion or small invasion degree, but when a sandstone target reservoir with obvious invasion exists, particularly when the contrast between the mud conductivity and the stratum conductivity is high and the borehole is large, invasion influences are often regarded as invasion profile errors caused by overcorrection or undercorrection for correcting or partially eliminating borehole influences, and incorrect flushing zone resistivity and stratum true resistivity results obtained by subsequent one-dimensional resistivity inversion. The method provided by the application considers invasion from the establishment of a theoretical model, so that the defect of the traditional borehole correction method is fundamentally overcome. Meanwhile, skin-seeking correction is not performed on the multi-frequency measured array induction well data in advance in order to reduce the dimension of a geometric factor library and the time complexity of an algorithm, so that richer well bore, invasion and undisturbed stratum conductivity information is reserved to a certain extent, conditions are created for the inversion process of environment correction with more parameters of an invasion model, and the optimization problem solving of the well bore correction environment variables is more stable and more accurate.

In some examples, step S110 may specifically further include steps S1101 to S1103:

s1101, acquiring an original conductivity curve of each subarray under the full frequency;

s1102, acquiring resolution difference information of an original conductivity curve between adjacent arrays;

and S1103, matching a lower resolution curve in the original conductivity curves between adjacent arrays based on the resolution difference information so as to obtain resolution matching apparent conductivity.

The basic principle of resolution matching is that high-resolution information in short subarray curves with shallow detection depth and high longitudinal resolution in an induction logging curve is extracted and then added to long subarray curves with deep detection depth and low longitudinal resolution, so that the longitudinal resolution of measurement curves of all subarrays is improved and unified on the basis of not changing the detection depth, and the two subarrays can participate in the inversion solving process of the same borehole correction model. For an array sensing instrument with M subarrays, we log the raw conductivity curves for M depths of investigation at different longitudinal resolutions. For the curves of the M detection depths, the resolution difference information of two adjacent curves is filtered, and then the difference information is matched to the corresponding curve according to the required resolution, so that the resolution matching can be completed. The specific method of resolution matching can be described by the mathematical formula as follows:

Wherein C is _l And C _h Low resolution and high resolution logs, respectively;is a logging curve with improved or enhanced resolution after resolution matching; />Is to make curve C of high resolution _h After passing through a filter f _h-l Obtained and C _l Log curves of the same longitudinal resolution, namely:

C _h the curve is obtained after filtering by a filterCurve possession and C _l Longitudinal resolution of the curve is the same, soThen it is the difference information of the resolution. Adding the difference information to C _l On the above, a resolution enhanced curve +.>

From the definition of the geometric factors and the detection depth, it can be known that C _h Curve obtained after filteringAlthough the longitudinal resolution is reduced, the detection depth is unchanged; similarly, when the handle is->Curve is added to C _l Obtained on the curve +.>The radial detection depth of the curve is also unchanged, i.e. +.>Is enhanced with C _h The curve is the same but the radial detection depth is unchanged from curve C _l Is the same. Thus, after resolution matching, the resolution of the logging curves with low resolution is enhanced, but the detection depth is not changed, so that the logging curves with different detection depths and the same longitudinal resolution can be obtained.

The principles and steps of resolution matching described above with abstract mathematical formulas describe how resolution matching is accomplished with actual log data to arrive at the final log. The M original conductivity curves corresponding to the M subarrays, which are different in detection depth and longitudinal resolution, are respectively recorded asWherein the superscript and the subscript denote longitudinal resolution and radial detection depth, respectively; the resolution matching filter is divided into 9 groups according to different background conductivities, and the filtering process of the filter is the same as that of soft focusingFor simplicity of description we do not consider the background conductivity again, but we mark the filters as resolution matched filters at each depth point, respectively, assuming these filters have been selectedThe curve obtained by filtering the high resolution curve with a filter is denoted +.>This shows that the resolution of the high resolution curves after filtering by the filter is reduced but the detection depth is not changed, and these curves are obtained by the following formula (7):

specifically, in practical application, the focusing composite curves can be matched to three kinds of longitudinal resolution of 1 foot, 2 feet and 4 feet. The enhanced resolution curve, i.e., the resulting log with a longitudinal resolution of 1 foot, is referred to herein as The curve with a longitudinal resolution of 2 feet is noted as The curve with longitudinal resolution of 4 feet is noted as The high resolution curves can be synthesized stepwise from the low resolution curves as follows:

first, a curve with a longitudinal resolution of 10 feet is synthesized into a curve with a resolution of 8 feet, namely:

a curve with a resolution of 8 feet was synthesized into a curve with a resolution of 6 feet, namely:

a curve with a resolution of 6 feet was synthesized into a curve with a resolution of 4 feet, namely:

the 4-foot resolution curve was synthesized as a 3-foot curve, namely:

the 3-foot resolution curve was synthesized as a 2-foot curve, namely:

the 2-foot resolution curve was synthesized as a 1-foot curve, namely:

after the steps, the resolution matching processing of the well logging curve is completed, curves with different detection depths of various resolutions are obtained, and in the actual well correction processing, the well logging curve with 2 feet can be taken as the resolution matching curve to participate in the optimization and inversion of the well correction variables at the back, so that the well correction is realized.

In some examples, step S120 may further include step S1201 value step S1203:

s1201, determining an initial value of a borehole correction variable according to the apparent conductivity, instrument geometric feature information, borehole geometric features, borehole diameters and slurry conductivity acquisition conditions and software processing parameter selection logic;

S1202, constructing a borehole geometry factor library and an invasion geometry factor library according to borehole correction variables;

and S1203, performing multi-component slicing interpolation calculation based on the borehole geometry factor library and the invasion geometry factor library to obtain the borehole geometry factor, the invasion geometry factor and the undisturbed stratum geometry factor.

Exemplary, based on the apparent conductivity, tool geometry information and borehole geometry, and software process parameter selection logic designed based primarily on measurements of via borehole diameter and mud conductivity in the tool in the same time borehole, the initial values cal of the 6 borehole correction variables are determined ₀ 、ecc ₀ 、σ _m0 、σ _xo0 、l _xo0 、σ _t0 Calculating the geometric factors by interpolating and looking up a table of a geometric factor library which is calculated in advance And invasion geometry factor->Determination of undisturbed formation factor/>

In some examples, the initial values of the borehole correction variables include an initial value of instrument eccentricity, ecc ₀ The geometric characteristic information of the instrument comprises an instrument outer diameter d _tool And centralizer size standoff, the wellbore geometry includes wellbore diameter cal, step S1201 may specifically include step S1201-a:

s1201-a: obtaining initial value ecc of instrument eccentricity according to ₀ :

Wherein d _tool For instrument outer diameter, stardoff is centralizer size, cal is wellbore diameter.

Exemplary, instrument eccentricity initial value ecc ₀ The instrument position of the borehole correction control option is set to 0 when centered and estimated as shown in equation (14) based on the instrument outer diameter, the borehole diameter and the centralizer size when the instrument position of the borehole correction control option is borehole wall-mounted.

In some examples, the initial values of the borehole correction variables include an initial value σ of the flushing charge conductivity _xo0 Step S1201 may specifically include step S1201-B:

S1201-B: obtaining initial value sigma of flushing electrification conductivity according to the following formula _xo0 ：

In the method, in the process of the invention,weight coefficient corresponding to the first target depth, < ->For different subarraysThe apparent conductivity under different frequencies is that M is the number of subarrays and N is the number of measuring frequencies;

exemplary initial value of the flushing charge conductivity sigma _xo0 Apparent conductivity of N measurement frequencies of M subarrays after resolution matching can be usedThe weighted average is estimated as shown in equation 15. Wherein the weight coefficient->The soft focus filter coefficients for the detection depth of 10in or 20in may be taken.

In some examples, the initial values of the borehole correction variables include an initial value σ of undisturbed formation conductivity _t0 Step S1201 may specifically include steps S1201-C:

S1201-C: obtaining an initial value sigma of the conductivity of the undisturbed stratum according to the following formula _t0 ：

In the method, in the process of the invention,weight coefficient corresponding to the second target depth, +.>For apparent conductivity of different subarrays at different frequencies, M is the number of subarrays, and N is the number of measuring frequencies. />

Exemplary, weight coefficientsThe soft focus filter coefficients for depth of detection depth 90in, or 120in, are typically taken. The initial value of the invasion radius is half of the well diameter _xo0 ＝0.5×cal。

In some examples, step S140 may specifically include: step S1401 to step S1402

S1401, respectively changing one or more of the borehole diameter, the instrument eccentricity, the slurry conductivity, the undisturbed stratum conductivity, the flushing charge conductivity and the flushing band radius 6 borehole correction parameters, and solving the 6 borehole correction parameters in an inversion mode based on an optimization method;

s1402, optimizing a linear equation set constructed based on a geometric factor theoretical model and the optimization formula to obtain the second target parameter;

wherein, the optimization formula is as follows:

the above linear system of equations is:

Illustratively, the penalty function may be constructed by constraints on the continuity and smoothness of 6 variables, for example, in the case of a stable mud system in the wellbore, the mud conductivity should be continuously changing, and during the lift-off logging, the mud conductivity should decrease as the wellbore temperature decreases, so that the continuity and monotonicity of the mud conductivity at the previous depth point can be taken into account to "penalty" the mud optimization result for deviations from the actual situation at the current depth point.

In some examples, the second target parameter includes a target wellbore diameter cal _inv Eccentricity ecc of target instrument _inv Target mud conductivity sigma _minv Target flush charge conductivity sigma _xoinv Target flush zone radius l _xoinv Target undisturbed formation conductivity sigma _tinv ；

Step S150 may include step S1501:

s1501: the array induction instrument was calibrated using the following:

In some examples, referring to fig. 6, the array induction instrument borehole correction step may specifically include steps S210 to S260:

s210, array induction logging apparent conductivity of N measuring frequencies of M subarrays by using the method from step S1101 to step S1103Matching to uniform resolution->

S220: determining initial values cal of 6 borehole correction variables using the methods described in S1201-A, S1201-B and S1201-C ₀ 、ecc ₀ 、σ _m0 、σ _xo0 、l _xo0 、σ _t0 Calculating the geometric factors by interpolating and looking up a table of a geometric factor library which is calculated in advanceAnd->Determine->

S230: solving 3 unknown variables sigma containing M×N linear equations in equation set (18) by least square method _m ,σ _xo ,σ _t In particular, σ when the mud conductivity measurement is accurate _m Also known as quantity, only 2 unknown variables sigma _xo ,σ _t Solution is required.

S240: updating initial values of 3 borehole correction conductivity variables by using the three variables obtained in the step S230, repeating the steps of the step S220, the step S230 and the step S240 until a preset convergence or suspension condition is met, and then entering the step S250;

s250: using the initial values of the 6 borehole correction variables after S250 ends, maintaining cal ₀ 、ecc ₀ 、σ _m0 、σ _t0 Invariably, inverting the borehole correction by solving the optimization problem of equation (19)Intrusion parameter sigma within a variable _xo0 、l _xo0 . Then turning to S220; repeatedly solving the linear equation set (18) and the optimization problem (19), and optimizing all 6 borehole correction variables in turn or simultaneously after solving the optimization problem for the second time until a preset convergence or suspension condition is met, and then entering into S260;

s260: inversion results cal of 6 borehole correction variables obtained according to S-optimization _inv 、ecc _inv 、σ _minv 、σ _xoinv 、l _xoinv 、σ _tinv Utilizing type (20)

6 parameter array induction borehole corrections with invasion model were completed.

In some examples, a 2ft resolution matching resistivity curve for the array sensing instrument AFIT-A in a certain well is shown in FIG. 7, and a 2ft resolution matching resistivity curve for the array sensing instrument AFIT-A in a certain well is shown in FIG. 8. In the figure, it can be seen that the resistivity curves of the detection depths of the array induction instrument AFIT-A are overlapped in the formation of the mudstone section, the resistivity curves of the detection depths of the sandstone permeable layer are different, and the curve relationship is correct and reasonable.

In a large borehole, as shown in FIG. 9 below, the well diameter of a well reaches 21in at 2633m, the well diameter of a well reaches 13in at 1306m as shown in FIG. 10, and the well section has high mud conductivity and formation conductivity contrast, the measurement value of the well mud at 20 ℃ is 0.7 ohm, the high-resistance formation resistivity is greater than 200 ohm, the borehole effect can still be well corrected, the curve relationship of different detection depths of the sandstone well section is consistent, the invasion indication is correct, and the curves of the impermeable mudstone sections coincide. Reliable borehole correction allows the tool to have a wider mud and formation adaptability and radial invasion resolution.

Referring to fig. 11, one embodiment of an array induction instrument wellbore correction device in an embodiment of the present application may include:

a first obtaining unit 21, configured to perform resolution matching based on the original conductivity curve to obtain multi-frequency multi-subarray formation apparent conductivity with the same resolution;

a second obtaining unit 22, configured to obtain an initial value of a borehole correction variable, a borehole geometry factor, an invasion geometry factor, and an undisturbed formation geometry factor coefficient, where the borehole correction variable includes a borehole diameter, a mud conductivity, an undisturbed formation conductivity, an instrument eccentricity, an invasion radius, and a flushing band conductivity;

A first determining unit 23, configured to determine a first target parameter by using a linear optimization algorithm or a method for solving a linear equation set based on an inductive measurement response model of a wellbore geometry factor and an invasion geometry factor, where the first target parameter is a first sub-optimal solution of a mud conductivity, a flushing charge conductivity, and an undisturbed formation conductivity, the inductive measurement response model is a pseudo-linear combination of a measurement apparent conductivity being the mud conductivity, the flushing charge conductivity, and the undisturbed formation conductivity, and a combination coefficient of the pseudo-linear combination is the wellbore geometry factor, the flushing band geometry factor, and the undisturbed formation geometry factor;

a second determining unit 24, configured to construct an objective function with an actual logging value matched with the simulated measurement response value of the induction measurement response model, and invert and determine a second objective parameter according to the first objective parameter by using a nonlinear optimization method, where the second objective parameter is a second sub-optimal solution of the wellbore diameter, the mud conductivity, the undisturbed formation conductivity, the instrument eccentricity, the invasion radius, and the flushing conductivity;

and the correction unit 25 is configured to obtain the target borehole geometry factor, the target flushing zone geometry factor, and the target undisturbed formation geometry factor by table lookup interpolation according to the second target parameter so as to implement borehole correction of the array induction instrument.

As shown in fig. 12, the embodiment of the present application further provides an electronic device 300, including a memory 310, a processor 320, and a computer program 311 stored in the memory 320 and capable of running on the processor, where the processor 320 executes the steps of any one of the methods of the array induction instrument wellbore correction methods described above.

Since the electronic device described in this embodiment is an apparatus for implementing the wellbore correction device of the array induction apparatus in this embodiment, based on the method described in this embodiment, those skilled in the art can understand the specific implementation of the electronic device of this embodiment and various modifications thereof, so how the electronic device implements the method in this embodiment will not be described in detail herein, and those skilled in the art should only use the apparatus for implementing the method in this embodiment, which is within the scope of protection of this application.

In a specific implementation, the computer program 311 may implement any of the embodiments corresponding to fig. 1 when executed by a processor.

In the foregoing embodiments, the descriptions of the embodiments are focused on, and for those portions of one embodiment that are not described in detail, reference may be made to the related descriptions of other embodiments.

It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present application may 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, and the like) having computer-usable program code embodied therein.

The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded computer, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

Embodiments of the present application also provide a computer program product comprising computer software instructions that, when run on a processing device, cause the processing device to perform the flow of the array induction instrument borehole correction method as in the corresponding embodiment of fig. 1.

The computer program product includes one or more computer instructions. When the computer program instructions are loaded and executed on a computer, the processes or functions in accordance with embodiments of the present application are produced in whole or in part. The computer may be a general purpose computer, a special purpose computer, a computer network, or other programmable apparatus. The computer instructions may be stored in a computer-readable storage medium or transmitted from one computer-readable storage medium to another computer-readable storage medium, for example, the computer instructions may be transmitted from one website, computer, server, or data center to another website, computer, server, or data center by a wired (e.g., coaxial cable, fiber optic, digital subscriber line (digital subscriber line, DSL)) or wireless (e.g., infrared, wireless, microwave, etc.). Computer readable storage media can be any available media that can be stored by a computer or data storage devices such as servers, data centers, etc. that contain an integration of one or more available media. Usable media may be magnetic media (e.g., floppy disks, hard disks, magnetic tapes), optical media (e.g., DVDs), or semiconductor media (e.g., solid State Disks (SSDs)), among others.

It will be clear to those skilled in the art that, for convenience and brevity of description, specific working procedures of the above-described systems, apparatuses and units may refer to corresponding procedures in the foregoing method embodiments, which are not repeated herein.

In the several embodiments provided in this application, it should be understood that the disclosed systems, apparatuses, and methods may be implemented in other ways. For example, the apparatus embodiments described above are merely illustrative, e.g., the division of elements is merely a logical functional division, and there may be additional divisions of actual implementation, e.g., multiple elements or components may be combined or integrated into another system, or some features may be omitted, or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be an indirect coupling or communication connection via some interfaces, devices or units, which may be in electrical, mechanical or other form.

The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed over a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

In addition, each functional unit in each embodiment of the present application may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit. The integrated units may be implemented in hardware or in software functional units.

The integrated units, if implemented in the form of software functional units and sold or used as stand-alone products, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present application may be embodied in essence or a part contributing to the prior art or all or part of the technical solution in the form of a software product stored in a storage medium, including several instructions to cause a computer device (which may be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the methods of the embodiments of the present application. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a random access Memory (Random Access Memory, RAM), a magnetic disk, or an optical disk, or other various media capable of storing program codes.

The above embodiments are only for illustrating the technical solution of the present application, and are not limiting thereof; although the present application has been described in detail with reference to the foregoing embodiments, it should be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the corresponding technical solutions.

Claims

1. A method of wellbore calibration for an array induction instrument, comprising:

determining a first target parameter by utilizing a linear optimization algorithm or a method for solving a linear equation set based on an induction measurement response model of a borehole geometry factor and an invasion geometry factor, wherein the first target parameter is a first sub-optimal solution of mud conductivity, flushing live conductivity and undisturbed formation conductivity, the induction measurement response model is a quasi-linear combination of the measurement apparent conductivity, namely mud conductivity, flushing live conductivity and undisturbed formation conductivity, and a combination coefficient of the quasi-linear combination is the borehole geometry factor, flushing zone geometry factor and undisturbed formation geometry factor;

2. The method of claim 1, wherein the performing a resolution match based on the raw conductivity curve to obtain multi-frequency multi-subarray formation apparent conductivities having the same resolution comprises:

and matching a lower-resolution curve in the original conductivity curves between adjacent arrays based on the resolution difference information to obtain resolution matching apparent conductivity.

3. The method of claim 1, wherein the obtaining initial values of borehole correction variables, borehole geometry factors, invasion geometry factors, and undisturbed formation geometry factor coefficients comprises:

and performing multi-component slicing interpolation calculation based on the borehole geometry factor library and the invasion geometry factor library to obtain the borehole geometry factor, the invasion geometry factor and the undisturbed stratum geometry factor.

4. The method of claim 3, wherein the initial value of the borehole correction variable comprises an initial value of instrument eccentricity, ecc ₀ The instrument geometric characteristic information comprises an instrument outer diameter d _tool And centralizer size standoff, the wellbore geometry comprising wellbore diameter cal;

the determining initial values of borehole correction variables based on the apparent conductivity, instrument geometry information, and borehole geometry and software process parameter selection logic comprises:

obtaining initial value ecc of instrument eccentricity according to ₀ :

5. The method of claim 3, wherein the initial value of the borehole correction variable comprises an initial value σ of a flush charge conductivity _xo0 ；

and/or the number of the groups of groups,

Determining initial values of borehole correction variables based on the apparent conductivity, instrument geometry information, and borehole geometry and software process parameter selection logic, comprising:

6. The method of claim 1, wherein constructing an objective function for matching the actual log value to the simulated measurement response value of the induction measurement response model, inverting a second objective parameter from the first objective parameter using a nonlinear optimization method, comprises:

Respectively changing one or more of the borehole diameter, the instrument eccentricity, the slurry conductivity, the undisturbed stratum conductivity, the flushing charge conductivity and the flushing band radius 6 borehole correction parameters, and solving the 6 borehole correction parameters in an inversion mode based on an optimization method;

optimizing a linear equation set constructed based on a geometric factor theoretical model and the optimization formula to acquire the second target parameter;

wherein, the optimization formula is as follows:

the system of linear equations is:

7. The method of claim 6, wherein the second target parameter comprises a target wellbore diameter cal _inv Eccentricity ecc of target instrument _inv Target mud conductivity sigma _minv Target flush charge conductivity sigma _xoinv Target flush zone radius l _xoinv Target undisturbed formation conductivity sigma _tinv ；

The obtaining the target borehole geometry factor, the target flushing zone geometry factor and the target undisturbed stratum geometry factor according to the second target parameter table lookup interpolation to realize borehole correction of the array induction instrument comprises the following steps:

the array sensing instrument was calibrated using the following:

8. An array induction instrument wellbore correction control device, comprising:

A first determining unit, configured to determine a first target parameter by using a linear optimization algorithm or a method for solving a linear equation set based on an inductive measurement response model of a wellbore geometry factor and an invasion geometry factor, where the first target parameter is a first sub-optimal solution of a mud conductivity, a flushing charge conductivity, and an undisturbed formation conductivity, the inductive measurement response model is a quasi-linear combination of a measurement apparent conductivity being the mud conductivity, the flushing charge conductivity, and the undisturbed formation conductivity, and a combination coefficient of the quasi-linear combination is the wellbore geometry factor, the flushing band geometry factor, and the undisturbed formation geometry factor;

9. An electronic device, comprising: a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor is configured to implement the steps of the array induction instrument wellbore correction method of any of claims 1 to 7 when executing the computer program stored in the memory.

10. A computer-readable storage medium having stored thereon a computer program, characterized by: the computer program when executed by a processor implements the array induction instrument wellbore correction method of any one of claims 1 to 7.