Array induction logging information anisotropy correction method, device, equipment and medium

Abstract

The invention provides an array induction logging information anisotropy correction method, device and equipment and a medium, is suitable for the technical field of oil exploration and development, and can improve the correction efficiency of array induction logging information. The method comprises the following steps: acquiring parameters of an array induction logging instrument; target interval information and adjacent well interval data corresponding to a target interval are obtained; according to the parameters of the array induction logging instrument and an array induction logging subarray forward modeling algorithm, an array induction logging composite signal is determined; constructing a uniform anisotropic medium inversion model, determining an anisotropic coefficient initial value of the current depth point according to target interval information and adjacent well data, and inputting the anisotropic coefficient initial value into the uniform anisotropic medium inversion model; and performing iterative inversion on the uniform anisotropic medium inversion model by using a gradient iterative inversion algorithm, and determining an anisotropic correction result of the actual measurement data of each detection mode at the current depth point.

Description

Technical Field

[0001] The present invention belongs to the technical field of petroleum exploration and development, and in particular relates to an anisotropy correction method, device, equipment and medium for array induction logging data. Background Art

[0002] Array induction logging is a key logging method for complex oil and gas reservoirs. However, due to resistivity anisotropy, it often exhibits anomalies such as curve uplift, disorder, and apparent resistivity "horizons" in highly deviated / horizontal wells, resulting in low logging interpretation compliance. Furthermore, the core "soft focusing" algorithm used in data processing for commercial array induction logging instruments has not yet been made public, and the filter coefficient libraries for the "soft focusing" algorithms used by instruments from different companies vary. Therefore, to address these issues, the present invention proposes a method for rapid anisotropy correction of array induction logging data based on a direct focusing algorithm for array induction logging synthetic signals. This method not only effectively avoids the issue of the "soft focusing" algorithm being undisclosed, but also increases data processing speed by more than 50 times. Because the radial integration geometric factors of synthetic signals are the same for instruments from different companies, this processing method is applicable to instruments from any company. This method effectively improves the resistivity extraction accuracy of reservoirs such as shale oil and tight sandstone, providing accurate electrical information for precise fluid identification. Summary of the Invention

[0003] The present invention provides a method, device, equipment and medium for anisotropy correction of array induction logging data, which can effectively improve the correction efficiency of array induction logging data.

[0004] In view of the above problems, the present invention adopts the following technical solutions:

[0005] In a first aspect, a method for anisotropy correction of array induction logging data is provided, comprising:

[0006] Obtaining parameters of an array induction logging instrument, including: the number of sub-arrays of the array induction logging instrument, the number of detection mode types, the number of turns of the transmitting coil, the number of turns of the receiving coil, the number of turns of the shielding coil, and the coil radius;

[0007] Obtain target layer information and corresponding adjacent well section data. Target layer information includes: lithology, mud resistivity, array induction logging curves, and wellbore trajectory. Adjacent well section data includes: GR curves, dual lateral curves, and SP curves.

[0008] Determining array induction logging synthetic signals according to array induction logging instrument parameters and array induction logging sub-array forward algorithm;

[0009] Construct a homogeneous anisotropic medium inversion model, and determine the initial value of the anisotropy coefficient at the current depth of the target layer based on the target layer information and adjacent well data, and input it into the homogeneous anisotropic medium inversion model;

[0010] The gradient iterative inversion algorithm is used to iteratively invert the inversion model of uniform anisotropic media to determine the anisotropy correction results of the measured data of each detection mode at the current depth point.

[0011] Optionally, determining the array induction logging synthetic signal according to array induction logging instrument parameters and an array induction logging sub-array forward algorithm includes:

[0012] Based on the measured signals of the M detection modes of the array induction logging tool and the parameters of the array induction logging tool, combined with the array induction logging subarray forward algorithm, the N subarray responses of the uniform medium model are determined, where N is the number of subarrays of the array induction logging tool and M is the number of detection modes of the array induction logging tool.

[0013] Determining radial integration geometric factors of N sub-arrays of array induction logging according to array induction logging instrument parameters;

[0014] According to the radial integration geometric factor objective function of the array induction logging synthetic signal, the radial integration geometric factors of the synthetic signals of the M detection modes are determined;

[0015] Determining N subarray coefficients under M detection modes according to radially integrated geometric factors of N subarrays of array induction logging and radially integrated geometric factors of synthetic signals of M detection modes;

[0016] According to the N sub-array responses and the N sub-array coefficients under the M detection modes, array induction logging synthetic signals of the M detection modes are determined.

[0017] Furthermore, the radial integral geometry factor GRS of the j-th subarray is j satisfy:

[0018]

[0019] Where ρ is the radial distance, r j The coil radius of the jth subarray, N T is the number of turns of the transmitting coil of the jth subarray, N R,j is the number of turns of the receiving coil of the jth subarray, N B,j is the number of turns of the shielding coil of the jth subarray, ρ T is the distance from the unit ring to the center of the transmitting coil, is the distance from the unit ring to the center of the j-th sub-array receiving coil, is the distance from the unit ring to the center of the j-th sub-array shielding coil, is the distance from the transmitting coil to the receiving coil of the jth subarray, is the distance from the transmitting coil to the shielding coil of the j-th subarray, 1≤j≤N;

[0020] Radial integral geometric factor objective function GRF of array induction logging synthetic signal of the i-th detection mode i satisfy:

[0021]

[0022] Among them, a x 、b x and c x is the control parameter corresponding to the radial integral geometric factor of the synthetic signal of the i-th detection mode, 1≤i≤M, 1≤x≤5;

[0023] The coefficients of the N subarrays in the i-th detection mode satisfy:

[0024]

[0025] in, are the coefficients of the N subarrays in the i-th detection mode, ρ1,ρ2,…,ρ N is the N radial distances in the i-th detection mode;

[0026] The synthetic signal R of array induction logging in the i-th detection mode i satisfy:

[0027]

[0028] Among them, SUB j is the j-th subarray response of the uniform medium model.

[0029] Optionally, a uniform anisotropic medium inversion model is constructed, and based on target layer information and adjacent well data, an initial value of the anisotropy coefficient at the current depth point of the target layer is determined and input into the uniform anisotropic medium inversion model, including:

[0030] Constructing a homogeneous anisotropic medium inversion model. The parameters of the homogeneous anisotropic medium inversion model include: horizontal resistivity, vertical resistivity, anisotropy coefficient, and well inclination angle;

[0031] Determine the reference resistivity, well inclination angle and measured resistivity of the current depth point of the target layer based on the target layer information and adjacent well section data;

[0032] Determine the initial value of the anisotropy coefficient at the current depth of the target layer according to the reference resistivity, well inclination and measured resistivity at the current depth of the target layer;

[0033] The initial value of the anisotropy coefficient of the current depth point of the target layer is input into the homogeneous anisotropic medium inversion model.

[0034] Furthermore, based on the target layer information and the adjacent well section data, the reference resistivity, well inclination angle and measured resistivity of the current depth point of the target layer are determined, including:

[0035] The average value of the deep lateral curve in the dual lateral curves is determined as the reference resistivity of the current depth point;

[0036] Determine the well inclination angle at the current depth point based on the wellbore trajectory;

[0037] The curve value with a vertical resolution of 2 ft (feet) and a radial detection depth of 10 inches in the array induction logging measured curve of the target layer is determined as the measured resistivity.

[0038] Furthermore, the anisotropy coefficient λ of the inversion model of homogeneous anisotropic media satisfies:

[0039]

[0040] Among them, R h is the horizontal resistivity, R v is the vertical resistivity;

[0041] The initial value of the anisotropy coefficient λ0 at the current depth point of the target layer satisfies:

[0042]

[0043] Among them, R res is the reference resistivity of the target layer at the current depth, R a is the measured resistivity at the current depth of the target layer, and θ is the well inclination angle at the current depth of the target layer.

[0044] In a second aspect, a device for anisotropy correction of array induction logging data is provided, comprising:

[0045] An acquisition module is used to acquire the parameters of the array induction logging instrument, which include: the number of sub-arrays of the array induction logging instrument, the number of detection mode types, the number of turns of the transmitting coil, the number of turns of the receiving coil, the number of turns of the shielding coil, and the coil radius;

[0046] The acquisition module is also used to obtain target layer information and the corresponding adjacent well section data. The target layer information includes: lithology, mud resistivity, array induction logging measured curve, wellbore trajectory, and the adjacent well section data includes: GR curve, dual lateral curve, SP curve;

[0047] The determination module is further used to determine the array induction logging synthetic signal based on the array induction logging instrument parameters and the array induction logging sub-array forward algorithm;

[0048] The determination module is also used to construct a uniform anisotropic medium inversion model and determine the initial value of the anisotropy coefficient at the current depth point of the target layer based on the target layer information and adjacent well data, and input it into the uniform anisotropic medium inversion model;

[0049] The determination module is also used to perform iterative inversion on the uniform anisotropic medium inversion model using a gradient iterative inversion algorithm to determine the anisotropy correction results of the measured data of each detection mode at the current depth point.

[0050] Optionally, the determination module is further configured to:

[0051] Based on the measured signals of the M detection modes of the array induction logging instrument and the parameters of the array induction logging instrument, combined with the array induction logging subarray forward algorithm, the N subarray responses of the uniform medium model are calculated, where N is the number of subarrays of the array induction logging instrument and M is the number of detection modes of the array induction logging instrument.

[0052] Determining radial integration geometric factors of N sub-arrays of array induction logging according to array induction logging instrument parameters;

[0053] According to the radial integration geometric factor objective function of the array induction logging synthetic signal, the radial integration geometric factors of the synthetic signals of the M detection modes are determined;

[0054] Determining N subarray coefficients under M detection modes according to radially integrated geometric factors of N subarrays of array induction logging and radially integrated geometric factors of synthetic signals of M detection modes;

[0055] According to the N sub-array responses and the N sub-array coefficients under the M detection modes, the array induction logging synthetic signals of the M detection modes are calculated.

[0056] Furthermore, the radial integral geometry factor GRS of the j-th subarray is j satisfy:

[0057]

[0058] Where ρ is the radial distance, r j The coil radius of the jth subarray, N T is the number of turns of the transmitting coil of the jth subarray, N R,j is the number of turns of the receiving coil of the jth subarray, N B,j is the number of turns of the shielding coil of the jth subarray, ρ T is the distance from the unit ring to the center of the transmitting coil, is the distance from the unit ring to the center of the j-th sub-array receiving coil, is the distance from the unit ring to the center of the j-th sub-array shielding coil, is the distance from the transmitting coil to the receiving coil of the jth subarray, is the distance from the transmitting coil to the shielding coil of the j-th subarray, 1≤j≤N;

[0059] Radial integral geometric factor objective function GRF of array induction logging synthetic signal of the i-th detection mode i satisfy:

[0060]

[0061] Among them, a x 、b x and c x is the control parameter corresponding to the radial integral geometric factor of the synthetic signal of the i-th detection mode, 1≤i≤M, 1≤x≤5;

[0062] The coefficients of the N subarrays in the i-th detection mode satisfy:

[0063]

[0064] in, are the coefficients of the N subarrays in the i-th detection mode, ρ2,…,ρ N is the N radial distances in the i-th detection mode;

[0065] The synthetic signal R of array induction logging in the i-th detection mode i satisfy:

[0066]

[0067] Among them, SUB j is the j-th subarray response of the uniform medium model.

[0068] Optionally, the determination module is further configured to:

[0069] Constructing a homogeneous anisotropic medium inversion model. The parameters of the homogeneous anisotropic medium inversion model include: horizontal resistivity, vertical resistivity, anisotropy coefficient, and well inclination angle;

[0070] Determine the reference resistivity, well inclination angle and measured resistivity of the current depth point of the target layer based on the target layer information and adjacent well section data;

[0071] Determine the initial value of the anisotropy coefficient at the current depth of the target layer according to the reference resistivity, well inclination and measured resistivity at the current depth of the target layer;

[0072] The initial value of the anisotropy coefficient of the current depth point of the target layer is input into the homogeneous anisotropic medium inversion model.

[0073] Furthermore, the module is further configured to:

[0074] The average value of the deep lateral curve in the dual lateral curves is determined as the reference resistivity of the current depth point;

[0075] Determine the well inclination angle at the current depth point based on the wellbore trajectory;

[0076] The curve value with a vertical resolution of 2 ft and a radial detection depth of 10 inches in the array induction logging measured curve of the target layer is determined as the measured resistivity.

[0077] Furthermore, the anisotropy coefficient λ of the inversion model of homogeneous anisotropic media satisfies:

[0078]

[0079] Among them, R h is the horizontal resistivity, R v is the vertical resistivity;

[0080] The initial value of the anisotropy coefficient λ0 at the current depth point of the target layer satisfies:

[0081]

[0082] Among them, R res is the reference resistivity of the target layer at the current depth, R a is the measured resistivity at the current depth of the target layer, and θ is the well inclination angle at the current depth of the target layer.

[0083] In a third aspect, an electronic device is provided, comprising: a processor coupled to a memory;

[0084] The processor is used to read and execute the program or instruction stored in the memory, so that the electronic device executes the anisotropy correction method for array induction logging data provided in the first aspect.

[0085] In a fourth aspect, a computer-readable storage medium is provided, which stores a program or instruction. When a computer reads and executes the program or instruction, the computer executes the array induction logging data anisotropy correction method provided in the first aspect.

[0086] Based on the anisotropy correction method, device, equipment and medium for array induction logging data provided by the present invention, the equivalent relationship between the radial integral geometric factor of the sub-array and the radial integral geometric factor of the synthetic signal can be utilized to establish a direct focusing method for the radial synthetic signal of array induction logging, thereby effectively avoiding the "soft focusing" algorithm black box problem of existing array induction logging instruments, increasing the data processing speed by more than 50 times, and thus improving the correction efficiency of array induction logging data.

[0087] In addition, since the synthetic signal integration geometric factors of instruments from different companies are the same, this method is applicable to instruments from any company and has stronger applicability.

[0088] Other features and advantages of the present invention will be described in the following description, and in part will become apparent from the description, or will be understood by practicing the present invention. The purpose and other advantages of the present invention can be realized and obtained by the structures pointed out in the description, claims and drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0089] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following is a brief introduction to the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0090] Figure 1 A schematic diagram of a flow chart of an anisotropy correction method for array induction logging data provided by an embodiment of the present invention;

[0091] Figure 2 An example of an anisotropy correction method for array induction logging data provided by an embodiment of the present invention;

[0092] Figure 3 Schematic diagram of the HDIL instrument structure and specific parameters provided by the embodiment of the present invention;

[0093] Figure 4 The array induction logging curves and adjacent well data for the X-well section provided in the embodiment of the present invention;

[0094] Figure 5 The radial integration geometry factor of the HDIL instrument subarray provided by the embodiment of the present invention;

[0095] Figure 6 The radial integration geometry factor of the synthetic signal of the HDIL instrument provided in the embodiment of the present invention;

[0096] Figure 7 A schematic diagram of an anisotropic inversion model provided in an embodiment of the present invention;

[0097] Figure 8 The anisotropy correction result of the array induction logging data of the X-well section provided by the embodiment of the present invention;

[0098] Figure 9 A schematic structural diagram of an anisotropy correction device for array induction logging data provided by an embodiment of the present invention;

[0099] Figure 10 A schematic structural diagram of an electronic device provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0100] To make the objectives, technical solutions, and advantages of the embodiments of the present invention more clear, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts shall fall within the scope of protection of the present invention.

[0101] First combine Figures 1-8 , detailed description of an embodiment of the present invention provides an anisotropy correction method for array induction logging data.

[0102] Figure 1 The present invention provides a flow chart of an anisotropy correction method for array induction logging data. Figure 1 As shown, the method includes:

[0103] S101, obtaining parameters of an array induction logging tool.

[0104] The parameters of the array induction logging instrument include: the number of sub-arrays of the array induction logging instrument, the number of detection mode types, the number of turns of the transmitting coil, the number of turns of the receiving coil, the number of turns of the shielding coil, and the coil radius.

[0105] S102, obtaining target layer information and adjacent well section data corresponding to the target layer.

[0106] The target layer information includes lithology, mud resistivity, array induction logging curves, and wellbore trajectory; the adjacent well section data includes GR curves, dual lateral curves, and SP curves.

[0107] S103 , determining an array induction logging synthetic signal according to array induction logging instrument parameters and an array induction logging sub-array forward algorithm.

[0108] Optionally, S103, determining an array induction logging synthetic signal according to array induction logging instrument parameters and an array induction logging sub-array forward algorithm, includes:

[0109] Based on the measured signals of the M detection modes of the array induction logging tool and the parameters of the array induction logging tool, combined with the array induction logging subarray forward algorithm, the N subarray responses of the uniform medium model are determined, where N is the number of subarrays of the array induction logging tool and M is the number of detection modes of the array induction logging tool.

[0110] Determining radial integration geometric factors of N sub-arrays of array induction logging according to array induction logging instrument parameters;

[0111] According to the radial integration geometric factor objective function of the array induction logging synthetic signal, the radial integration geometric factors of the synthetic signals of the M detection modes are determined;

[0112] Determining N subarray coefficients under M detection modes according to radially integrated geometric factors of N subarrays of array induction logging and radially integrated geometric factors of synthetic signals of M detection modes;

[0113] According to the N sub-array responses and the N sub-array coefficients under the M detection modes, array induction logging synthetic signals of the M detection modes are determined.

[0114] Furthermore, the radial integral geometry factor GRS of the j-th subarray is j satisfy:

[0115]

[0116] Where ρ is the radial distance, r j The coil radius of the jth subarray, N T is the number of turns of the transmitting coil of the jth subarray, N R,j is the number of turns of the receiving coil of the jth subarray, N B,j is the number of turns of the shielding coil of the jth subarray, ρ T is the distance from the unit ring to the center of the transmitting coil, is the distance from the unit ring to the center of the j-th sub-array receiving coil, is the distance from the unit ring to the center of the j-th sub-array shielding coil, is the distance from the transmitting coil to the receiving coil of the jth subarray, is the distance from the transmitting coil to the shielding coil of the j-th subarray, 1≤j≤N;

[0117] The unit ring can be understood as follows: the medium surrounding the coil can be considered to be composed of countless circular rings with a cross-sectional area of ​​drdz and a radius of r. These rings are usually called unit rings. These unit rings are similar to closed coils. Under the influence of the alternating electromagnetic field of the transmitting coil, they generate induced electromotive force and induced current. The magnitude of the induced electromotive force and induced current is related to the conductivity of the unit rings.

[0118] Radial integral geometric factor objective function GRF of array induction logging synthetic signal of the i-th detection mode i satisfy:

[0119]

[0120] Among them, a x 、b x and c x is the control parameter corresponding to the radial integral geometric factor of the synthetic signal of the i-th detection mode, 1≤i≤M, 1≤x≤5, where 5 refers to the control parameter a in the calculation process of the radial integral geometric factor of the synthetic signal of each detection mode x 、b x and c x There are 5 in each;

[0121] The coefficients of the N subarrays in the i-th detection mode satisfy:

[0122]

[0123] in, are the coefficients of the N subarrays in the i-th detection mode, ρ1,ρ2,…,ρ N is the N radial distances in the i-th detection mode;

[0124] The synthetic signal R of array induction logging in the i-th detection mode i satisfy:

[0125]

[0126] Among them, SUB j is the j-th subarray response of the uniform medium model.

[0127] S104: constructing a uniform anisotropic medium inversion model, and determining an initial value of anisotropy coefficient at the current depth of the target layer based on target layer information and adjacent well data, and inputting the initial value into the uniform anisotropic medium inversion model.

[0128] Optionally, S104, constructing a uniform anisotropic medium inversion model, and determining an initial value of anisotropy coefficient at the current depth of the target layer based on target layer information and adjacent well data, and inputting the initial value into the uniform anisotropic medium inversion model, including:

[0129] Constructing a homogeneous anisotropic medium inversion model. The parameters of the homogeneous anisotropic medium inversion model include: horizontal resistivity, vertical resistivity, anisotropy coefficient, and well inclination angle;

[0130] Determine the reference resistivity, well inclination angle and measured resistivity of the current depth point of the target layer based on the target layer information and adjacent well section data;

[0131] Determine the initial value of the anisotropy coefficient at the current depth of the target layer according to the reference resistivity, well inclination and measured resistivity at the current depth of the target layer;

[0132] The initial value of the anisotropy coefficient of the current depth point of the target layer is input into the homogeneous anisotropic medium inversion model.

[0133] Furthermore, based on the target layer information and the adjacent well section data, the reference resistivity, well inclination angle and measured resistivity of the current depth point of the target layer are determined, including:

[0134] The average value of the deep lateral curve in the dual lateral curves is determined as the reference resistivity of the current depth point;

[0135] Determine the well inclination angle at the current depth point based on the wellbore trajectory;

[0136] The curve value with a vertical resolution of 2 ft and a radial detection depth of 10 inches in the array induction logging measured curve of the target layer is determined as the measured resistivity.

[0137] Furthermore, the anisotropy coefficient λ of the inversion model of homogeneous anisotropic media satisfies:

[0138]

[0139] Among them, R h is the horizontal resistivity, R v is the vertical resistivity;

[0140] The initial value of the anisotropy coefficient λ0 at the current depth point of the target layer satisfies:

[0141]

[0142] Among them, R res is the reference resistivity of the target layer at the current depth, R a is the measured resistivity at the current depth of the target layer, and θ is the well inclination angle at the current depth of the target layer.

[0143] S105 , using a gradient iterative inversion algorithm, iteratively inverts the homogeneous anisotropic medium inversion model to determine anisotropy correction results of the measured data of each detection mode at the current depth point.

[0144] The following combination Figure 2-Figure 8 .

[0145] like Figure 2 As shown, a method for rapid anisotropy correction of array induction logging data includes the following steps:

[0146] Step 1: Obtain the parameters of the array induction logging instrument, such as the number of sub-arrays, instrument operating frequency, number of turns of the transmitting / receiving / shielding coils, the distance between the receiving / shielding coils and the transmitting coil of different sub-arrays, and the coil radius.

[0147] Take Baker-Atlas's high-resolution array induction logging tool HDIL as an example. The tool has 7 sub-array coil systems and 8 operating frequencies (10kHz, 30kHz, 50kHz, 70kHz, 90kHz, 110kHz, 130kHz, and 150kHz). The measured signal has 6 different detection modes, which can be expressed as R 10inch ,R 20inch ,R 30inch ,R 60inch ,R 90inch ,R 120inch , the specific structure and parameters are as follows Figure 3 shown.

[0148] Step 2: Obtain target layer information, such as lithology, mud resistivity, array induction logging curves, wellbore trajectory, etc., and obtain adjacent well section data corresponding to the target layer, such as GR curves, dual lateral curves, SP curves, etc.

[0149] Taking the well section X as an example, the target layer section of this well section and the corresponding adjacent well section information are as follows: Figure 4 The X section of Well X is a shale oil reservoir, which typically exhibits strong anisotropy. Analysis of measured data shows that as the well inclination increases, the array induction logging curve is affected by anisotropy and exhibits an overall uplift, which is consistent with the anisotropic response law.

[0150] Step 3: Using the instrument parameters obtained in step 1 and combining it with the array induction logging sub-array forward algorithm, a direct focusing method for array induction logging synthetic signals is constructed.

[0151] The following describes step 3 in further detail:

[0152] Step 3.1: From step 1, we know that the array induction logging tool has 7 sub-arrays, which can measure signals in 6 detection modes. The number of turns of the transmitting coil is N.T =100, the number of turns of the receiving coil is N R,j , the number of turns of the shielding coil is N B,j , j represents the jth subarray, j=1,2,…,7, and the coil radius is r=0.0031m. Using the above parameters, the array induction logging subarray forward algorithm can be used to calculate the 7 subarray responses of the uniform medium model SUB j .

[0153] Step 3.2, taking subarray 1 as an example, the number of turns of the receiving coil of subarray 1 is N R,1 =12, the distance from the receiving coil to the transmitting coil is Shielding coil turns N B,1 =6, the distance from the shielding coil to the transmitting coil is Substituting the parameters into the radial integral geometric factor calculation formula, we can obtain:

[0154]

[0155] Taking the range of ρ as [0,3]m, calculate GRS1, and get the radial integral geometric factor of subarray 1 in the radial range [0,3]m. Similarly, perform the above calculations on the 7 subarrays respectively, and get the radial integral geometric factors of the 7 subarrays, as shown in the figure below: Figure 5 shown.

[0156] Step 3.3: Use the radial integration geometric factor objective function of the array induction logging synthetic signal Determine the radial integration geometric factor of the synthetic signal of the six detection modes, where i represents the i-th detection mode, i = 1, 2, ..., 6, a x 、b x and c x is the control parameter of the radial integral geometric factor fitting formula of the synthetic signal, which can be determined according to the requirements of different radial synthetic depths. Taking detection mode 1 as an example, its radial integral geometric factor control parameter is a x =[0.961,0.5647,0.1897,0.5842,0.2686], b x =[1.188,0.4252,0.6053,0.737,0.2637], c x =[0.4131,0.1638,0.146,0.2605,0.1198], taking the range of ρ as [0,3]m, and calculating GRF1, we can get the radial integral geometric factor of detection mode 1 in the radial range [0,3]m. Perform the above calculations on the six detection modes respectively, and we can get the radial integral geometric factors of the six detection modes, such as Figure 6 shown.

[0157] Step 3.4: Using the formulas in steps 3.2 and 3.3, select 7 sampling points ρ at equal intervals within the radial detection depth range of [0,3]m. j , and record the radial integral geometric factor values ​​of detection mode 1 and 7 sub-arrays corresponding to the 7 sampling points respectively, and construct the equivalent relationship between the radial integral geometric factor of the 7 sub-arrays and detection mode 1, and the following matrix form can be obtained:

[0158] Through polynomial fitting calculation, the coefficients of the 7 sub-arrays corresponding to detection mode 1 can be obtained: w1 = 0.1175, w2 = 0.1644, w3 = 2.9657, w4 = -2.4836, w5 = 0.1274, w6 = 0.1054, w7 = 0. In order to ensure that the apparent conductivity calculated in the homogeneous medium is equal to the true conductivity of the formation, these coefficients need to meet the normalization condition, that is, Repeat this step for the six detection modes to obtain coefficients of the seven sub-arrays corresponding to the six detection depths.

[0159] In step 3.5, the array induction logging synthetic signal of the required detection mode can be directly calculated by using the 7 sub-array responses calculated in step 3.1 and the coefficients of the 7 sub-arrays corresponding to different detection modes obtained in step 3.4, for example, R 10inch =SUB1*0.1175+SUB2*0.1644+…+SUB7*0.

[0160] Step 4: Construct an inversion model for a uniform anisotropic medium. Based on the target layer and adjacent well data obtained in step 2, determine the initial value of the anisotropy coefficient at the current depth point and input it into the inversion model.

[0161] Step 4 is described in further detail below:

[0162] Step 4.1, build Figure 7 The inversion model of the homogeneous anisotropic medium shown in the figure includes the horizontal resistivity R h , vertical resistivity R v , anisotropy coefficient Well inclination angle θ.

[0163] Step 4.2: Use the average value of the deep lateral curve in the dual lateral curves obtained in step 2 as the reference resistivity R at the current depth point. res =7Ω·m, obtain the well inclination angle at the current depth point according to the well trajectory, assuming that the well inclination angle at the current depth point θ0 = 85°, and use the curve value of the array induction logging measured curve of the target layer with a vertical resolution of 2ft and a radial detection depth of 10inch as the measured resistivity value R a =13Ω·m.

[0164] Step 4.3: Combine the target layer information obtained in step 4.2 and use the anisotropy coefficient calculation formula The initial anisotropy coefficient of the target layer segment can be obtained as λ0 = 1.6287.

[0165] In step 4.4, if the inversion result of the previous depth point of the current depth point is known, the inversion result of the previous depth point and the initial anisotropy coefficient λ0 determined in step 4.3 are sequentially input as the initial inversion values ​​of the current depth point into the inversion model constructed in step 4.1 for inversion, and the result with the smaller error between the two and the measured data is selected as the inversion result of the current depth point; otherwise, only the initial anisotropy coefficient λ0 determined in step 4.3 is input as the initial inversion value of the current depth point into the inversion model constructed in step 4.1 for inversion, thereby obtaining the inversion result of the current depth point.

[0166] Step 5: Use the gradient iterative inversion algorithm to iteratively invert the inversion model constructed in step 4 to obtain the horizontal resistivity R after anisotropy correction of the measured data in the single detection mode at the current depth point. h,10inch and vertical resistivity R v,10inch .

[0167] Similarly, all detection modes can be processed in a loop using step 5 to obtain the horizontal resistivity R after anisotropy correction of the measured data of all detection modes at the current depth point. h,10inch -R h,120inch and vertical resistivity R v,10inch -R v,120inch .

[0168] It is easy to understand that the operation in step 5 can be applied to all depth points of the target layer to achieve a rapid correction of anisotropy of the array induction logging curve of the target layer, such as Figure 8 shown.

[0169] Based on the anisotropy correction method for array induction logging data provided by the present invention, the equivalent relationship between the radial integral geometric factor of the sub-array and the radial integral geometric factor of the synthetic signal can be utilized to establish a direct focusing method for the radial synthetic signal of array induction logging, thereby effectively avoiding the "soft focusing" algorithm black box problem of existing array induction logging instruments, increasing the data processing speed by more than 50 times, and thus improving the correction efficiency of array induction logging data.

[0170] In addition, this method can effectively improve the resistivity extraction accuracy of reservoirs such as shale oil and tight sandstone, and provide accurate electrical information for precise identification of fluids.

[0171] In addition, since the synthetic signal integration geometric factors of instruments from different companies are the same, this method is applicable to instruments from any company and has stronger applicability.

[0172] Combined with the above Figures 1-8 The anisotropy correction method for array induction logging data provided by the embodiment of the present invention is described in detail. Figure 9 and Figure 10 An array induction logging data anisotropy correction device 900 and an electronic device 1000 provided in an embodiment of the present invention are respectively described, and are used to execute the array induction logging data anisotropy correction method provided in the above method embodiment.

[0173] like Figure 9 As shown, the array induction logging data anisotropy correction device 900 includes:

[0174] In a second aspect, an array induction logging data anisotropy correction device 900 is provided, comprising:

[0175] An acquisition module 901 is used to acquire parameters of an array induction logging instrument, including the number of subarrays, the number of detection modes, the number of turns of the transmitting coil, the number of turns of the receiving coil, the number of turns of the shielding coil, and the coil radius.

[0176] The acquisition module 901 is further used to obtain target layer information and adjacent well section data corresponding to the target layer. The target layer information includes: lithology, mud resistivity, array induction logging measured curve, wellbore trajectory, and adjacent well section data includes: GR curve, dual lateral curve, SP curve;

[0177] The determination module 902 is further configured to determine the array induction logging synthetic signal based on the array induction logging instrument parameters and the array induction logging sub-array forward algorithm;

[0178] The determination module 902 is further used to construct a uniform anisotropic medium inversion model and determine the initial value of the anisotropy coefficient at the current depth point of the target layer based on the target layer information and the data of the adjacent wells, and input it into the uniform anisotropic medium inversion model;

[0179] The determination module 902 is further configured to perform iterative inversion on the homogeneous anisotropic medium inversion model using a gradient iterative inversion algorithm to determine anisotropy correction results of the measured data of each detection mode at the current depth point.

[0180] Optionally, the determining module 902 is further configured to:

[0181] Based on the measured signals of the M detection modes of the array induction logging instrument and the parameters of the array induction logging instrument, combined with the array induction logging subarray forward algorithm, the N subarray responses of the uniform medium model are calculated, where N is the number of subarrays of the array induction logging instrument and M is the number of detection modes of the array induction logging instrument.

[0182] Determining radial integration geometric factors of N sub-arrays of array induction logging according to array induction logging instrument parameters;

[0183] According to the radial integration geometric factor objective function of the array induction logging synthetic signal, the radial integration geometric factors of the synthetic signals of the M detection modes are determined;

[0184] Determining N subarray coefficients under M detection modes according to radially integrated geometric factors of N subarrays of array induction logging and radially integrated geometric factors of synthetic signals of M detection modes;

[0185] According to the N sub-array responses and the N sub-array coefficients under the M detection modes, the array induction logging synthetic signals of the M detection modes are calculated.

[0186] Furthermore, the radial integral geometry factor GRS of the j-th subarray is j satisfy:

[0187]

[0188] Where ρ is the radial distance, r j The coil radius of the jth subarray, N T is the number of turns of the transmitting coil of the jth subarray, N R,j is the number of turns of the receiving coil of the jth subarray, N B,j is the number of turns of the shielding coil of the jth subarray, ρ T is the distance from the unit ring to the center of the transmitting coil, is the distance from the unit ring to the center of the j-th sub-array receiving coil, is the distance from the unit ring to the center of the j-th sub-array shielding coil, is the distance from the transmitting coil to the receiving coil of the jth subarray, is the distance from the transmitting coil to the shielding coil of the j-th subarray, 1≤j≤N;

[0189] Radial integral geometric factor objective function GRF of array induction logging synthetic signal of the i-th detection mode i satisfy:

[0190]

[0191] Among them, a x 、b x and c x is the control parameter corresponding to the radial integral geometric factor of the synthetic signal of the i-th detection mode, 1≤i≤M, 1≤x≤5;

[0192] The coefficients of the N subarrays in the i-th detection mode satisfy:

[0193]

[0194] in, are the coefficients of the N subarrays in the i-th detection mode, ρ1,ρ2,…,ρ N is the N radial distances in the i-th detection mode;

[0195] The synthetic signal R of array induction logging in the i-th detection mode i satisfy:

[0196]

[0197] Among them, SUB j is the j-th subarray response of the uniform medium model.

[0198] Optionally, the determining module 902 is further configured to:

[0199] Constructing a homogeneous anisotropic medium inversion model. The parameters of the homogeneous anisotropic medium inversion model include: horizontal resistivity, vertical resistivity, anisotropy coefficient, and well inclination angle;

[0200] Determine the reference resistivity, well inclination angle and measured resistivity of the current depth point of the target layer based on the target layer information and adjacent well section data;

[0201] Determine the initial value of the anisotropy coefficient at the current depth of the target layer according to the reference resistivity, well inclination and measured resistivity at the current depth of the target layer;

[0202] The initial value of the anisotropy coefficient of the current depth point of the target layer is input into the homogeneous anisotropic medium inversion model.

[0203] Furthermore, the determination module 902 is further configured to:

[0204] The average value of the deep lateral curve in the dual lateral curves is determined as the reference resistivity of the current depth point;

[0205] Determine the well inclination angle at the current depth point based on the wellbore trajectory;

[0206] The curve value with a vertical resolution of 2 ft and a radial detection depth of 10 inches in the array induction logging measured curve of the target layer is determined as the measured resistivity.

[0207] Furthermore, the anisotropy coefficient λ of the inversion model of homogeneous anisotropic media satisfies:

[0208]

[0209] Among them, R h is the horizontal resistivity, R vis the vertical resistivity;

[0210] The initial value of the anisotropy coefficient λ0 at the current depth point of the target layer satisfies:

[0211]

[0212] Among them, R res is the reference resistivity of the target layer at the current depth, R a is the measured resistivity at the current depth of the target layer, and θ is the well inclination angle at the current depth of the target layer.

[0213] like Figure 10 As shown, the electronic device 1000 includes: a processor 1001, the processor 1001 is coupled to a memory 1002;

[0214] The processor 1001 is configured to read and execute the program or instruction stored in the memory 1002 , so that the electronic device 1000 executes the anisotropy correction method for array induction logging data provided in the above method embodiment.

[0215] Optionally, the electronic device 1000 may further include a transceiver 1003 for the electronic device 1000 to communicate with other devices.

[0216] It should be noted that, for the sake of convenience, Figure 9 and Figure 10 Only the main components of the array induction logging data anisotropy correction device 900 and the electronic device 1000 are shown. In practical applications, the array induction logging data anisotropy correction device 900 and the electronic device 1000 may also include components or assemblies not shown in the figures.

[0217] An embodiment of the present invention further provides a computer-readable storage medium storing a program or instruction. When a computer reads and executes the program or instruction, the computer executes the array induction logging data anisotropy correction method provided in the above method embodiment.

[0218] Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein; and these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for anisotropy correction of array induction logging data, characterized in that: include: Acquiring parameters of an array induction logging instrument, wherein the parameters of the array induction logging instrument include: the number of sub-arrays of the array induction logging instrument, the number of detection mode types, the number of turns of the transmitting coil, the number of turns of the receiving coil, the number of turns of the shielding coil, and the coil radius; Obtain target layer information and adjacent well section data corresponding to the target layer section, wherein the target layer section information includes: lithology, mud resistivity, array induction logging measured curve, wellbore trajectory, and the adjacent well section data includes: GR curve, dual lateral curve, and SP curve; Determining array induction logging synthetic signals according to array induction logging instrument parameters and array induction logging sub-array forward algorithm; Constructing a uniform anisotropic medium inversion model, and determining an initial value of anisotropy coefficient at a current depth point of the target layer based on the target layer information and the adjacent well data, and inputting the initial value into the uniform anisotropic medium inversion model; The gradient iterative inversion algorithm is used to iteratively invert the homogeneous anisotropic medium inversion model to determine the anisotropy correction results of the measured data of each detection mode at the current depth point.

2. The anisotropy correction method for array induction logging data according to claim 1, characterized in that: Determining the array induction logging synthetic signal according to the array induction logging instrument parameters and the array induction logging sub-array forward algorithm includes: Determining N sub-array responses of a uniform medium model based on measured signals of M detection modes of the array induction logging instrument and parameters of the array induction logging instrument in combination with an array induction logging sub-array forward algorithm, where N is the number of sub-arrays of the array induction logging instrument and M is the number of detection modes of the array induction logging instrument; Determining radial integration geometric factors of N sub-arrays of array induction logging according to the array induction logging instrument parameters; According to the radial integration geometric factor objective function of the array induction logging synthetic signal, the radial integration geometric factors of the synthetic signals of the M detection modes are determined; Determining N subarray coefficients under M detection modes according to radially integrated geometric factors of the N subarrays of the array induction logging and radially integrated geometric factors of the synthetic signals of the M detection modes; The array induction logging synthetic signals of the M detection modes are determined according to the N sub-array responses and the N sub-array coefficients under the M detection modes.

3. The anisotropy correction method for array induction logging data according to claim 2, characterized in that: Radial integral geometric factor GRS of the jth subarray j satisfy: Where ρ is the radial distance, r j The coil radius of the jth subarray, N T is the number of turns of the transmitting coil of the jth subarray, N R,j is the number of turns of the receiving coil of the jth subarray, N B,j is the number of turns of the shielding coil of the jth subarray, ρ T is the distance from the unit ring to the center of the transmitting coil, is the distance from the unit ring to the center of the j-th sub-array receiving coil, is the distance from the unit ring to the center of the j-th sub-array shielding coil, is the distance from the transmitting coil to the receiving coil of the jth subarray, is the distance from the transmitting coil to the shielding coil of the j-th subarray, 1≤j≤N; Radial integral geometric factor objective function GRF of array induction logging synthetic signal of the i-th detection mode i satisfy: Among them, a x 、b x and c x is the control parameter corresponding to the radial integral geometric factor of the synthetic signal of the i-th detection mode, 1≤i≤M, 1≤x≤5; The coefficients of the N subarrays in the i-th detection mode satisfy: in, are the coefficients of the N subarrays in the i-th detection mode, ρ2,…,ρ N is the N radial distances in the i-th detection mode; The synthetic signal R of array induction logging in the i-th detection mode i satisfy: Among them, SUB j is the j-th subarray response of the uniform medium model.

4. The anisotropy correction method for array induction logging data according to any one of claims 1 to 3, characterized in that: The method of constructing a uniform anisotropic medium inversion model and determining an initial value of anisotropy coefficient at a current depth point of the target layer according to the target layer information and the adjacent well data, and inputting the initial value into the uniform anisotropic medium inversion model, includes: Constructing a homogeneous anisotropic medium inversion model, wherein the parameters of the homogeneous anisotropic medium inversion model include: horizontal resistivity, vertical resistivity, anisotropy coefficient, and well inclination angle; Determine the reference resistivity, well inclination angle and measured resistivity of the current depth point of the target layer interval according to the target layer interval information and the adjacent well interval data; Determining an initial value of an anisotropy coefficient at the current depth point of the target layer interval according to a reference resistivity, a well inclination angle, and a measured resistivity at the current depth point of the target layer interval; The initial value of the anisotropy coefficient of the current depth point of the target layer is input into the uniform anisotropic medium inversion model.

5. The anisotropy correction method for array induction logging data according to claim 3, characterized in that: Determining the reference resistivity, well inclination angle, and measured resistivity of the current depth point of the target layer interval based on the target layer interval information and the adjacent well interval data includes: Determining an average value of the deep lateral curves in the dual lateral curves as the reference resistivity at the current depth point; Determine the well inclination angle at the current depth point according to the wellbore trajectory; The curve value of the array induction logging measured curve of the target layer with a vertical resolution of 2 ft and a radial detection depth of 10 inches is determined as the measured resistivity.

6. The anisotropy correction method for array induction logging data according to claim 5, characterized in that: The anisotropy coefficient λ of the homogeneous anisotropic medium inversion model satisfies: Among them, R h is the horizontal resistivity, R v is the vertical resistivity; The initial value λ0 of the initial anisotropy coefficient at the current depth point of the target layer section satisfies: Among them, R res is the reference resistivity of the target layer at the current depth, R a is the measured resistivity at the current depth point of the target layer, and θ is the well inclination angle at the current depth point of the target layer.

7. An array induction logging data anisotropy correction device, characterized in that: include: an acquisition module, configured to acquire parameters of an array induction logging instrument, wherein the parameters of the array induction logging instrument include: the number of subarrays of the array induction logging instrument, the number of detection mode types, the number of turns of the transmitting coil, the number of turns of the receiving coil, the number of turns of the shielding coil, and the coil radius; The acquisition module is further configured to acquire target layer information and adjacent well section data corresponding to the target layer section, wherein the target layer section information includes: lithology, mud resistivity, array induction logging measured curve, and wellbore trajectory, and the adjacent well section data includes: GR curve, dual lateral curve, and SP curve; The determination module is further used to determine the array induction logging synthetic signal based on the array induction logging instrument parameters and the array induction logging sub-array forward algorithm; The determination module is further configured to construct a uniform anisotropic medium inversion model, and determine an initial value of anisotropy coefficient at a current depth point of the target layer based on the target layer information and the adjacent well data, and input the initial value into the uniform anisotropic medium inversion model; The determination module is further configured to perform iterative inversion on the homogeneous anisotropic medium inversion model using a gradient iterative inversion algorithm to determine anisotropy correction results of measured data of each detection mode at the current depth point.

8. The anisotropy correction device for array induction logging data according to claim 7, characterized in that: The determining module is further configured to: Calculating N subarray responses of a uniform medium model based on measured signals of M detection modes of the array induction logging instrument and parameters of the array induction logging instrument in combination with an array induction logging subarray forward algorithm, where N is the number of subarrays of the array induction logging instrument and M is the number of detection modes of the array induction logging instrument; Determining radial integration geometric factors of N sub-arrays of array induction logging according to the array induction logging instrument parameters; According to the radial integration geometric factor objective function of the array induction logging synthetic signal, the radial integration geometric factors of the synthetic signals of the M detection modes are determined; Determining N subarray coefficients under M detection modes according to radially integrated geometric factors of the N subarrays of the array induction logging and radially integrated geometric factors of the synthetic signals of the M detection modes; Calculate array induction logging synthetic signals of the M detection modes based on the N subarray responses and the N subarray coefficients under the M detection modes.

9. The array induction logging data anisotropy correction device according to claim 8, characterized in that: Radial integral geometric factor GRS of the jth subarray j satisfy: Where ρ is the radial distance, r j The coil radius of the jth subarray, N T is the number of turns of the transmitting coil of the jth subarray, N R,j is the number of turns of the receiving coil of the jth subarray, N B,j is the number of turns of the shielding coil of the jth subarray, ρ T is the distance from the unit ring to the center of the transmitting coil, is the distance from the unit ring to the center of the j-th sub-array receiving coil, is the distance from the unit ring to the center of the j-th sub-array shielding coil, is the distance from the transmitting coil to the receiving coil of the jth subarray, is the distance from the transmitting coil to the shielding coil of the j-th subarray, 1≤j≤N; Radial integral geometric factor objective function GRF of array induction logging synthetic signal of the i-th detection mode i satisfy: Among them, a x 、b x and c x is the control parameter corresponding to the radial integral geometric factor of the synthetic signal of the i-th detection mode, 1≤i≤M, 1≤x≤5; The coefficients of the N subarrays in the i-th detection mode satisfy: in, are the coefficients of the N subarrays in the i-th detection mode, ρ2,…,ρ N is the N radial distances in the i-th detection mode; The synthetic signal R of array induction logging in the i-th detection mode i satisfy: Among them, SUB j is the j-th subarray response of the uniform medium model.

10. The anisotropy correction device for array induction logging data according to any one of claims 7 to 9, characterized in that: The determining module is further configured to: Constructing a homogeneous anisotropic medium inversion model, wherein the parameters of the homogeneous anisotropic medium inversion model include: horizontal resistivity, vertical resistivity, anisotropy coefficient, and well inclination angle; Determine the reference resistivity, well inclination angle and measured resistivity of the current depth point of the target layer interval according to the target layer interval information and the adjacent well interval data; Determining an initial value of an anisotropy coefficient at the current depth point of the target layer interval according to a reference resistivity, a well inclination angle, and a measured resistivity at the current depth point of the target layer interval; The initial value of the anisotropy coefficient of the current depth point of the target layer is input into the uniform anisotropic medium inversion model.

11. The array induction logging data anisotropy correction device according to claim 9, characterized in that: The determining module is further configured to: Determining an average value of the deep lateral curves in the dual lateral curves as the reference resistivity at the current depth point; Determine the well inclination angle at the current depth point according to the wellbore trajectory; The curve value of the array induction logging measured curve of the target layer with a vertical resolution of 2 ft and a radial detection depth of 10 inches is determined as the measured resistivity.

12. The array induction logging data anisotropy correction device according to claim 11, characterized in that: The anisotropy coefficient λ of the homogeneous anisotropic medium inversion model satisfies: Among them, R h is the horizontal resistivity, R v is the vertical resistivity; The initial value λ0 of the initial anisotropy coefficient at the current depth point of the target layer section satisfies: Among them, R res is the reference resistivity of the target layer at the current depth, R a is the measured resistivity at the current depth point of the target layer, and θ is the well inclination angle at the current depth point of the target layer.

13. An electronic device, characterized in that: include: a processor coupled to the memory; The processor is configured to read and execute the program or instruction stored in the memory, so that the electronic device executes the anisotropy correction method for array induction logging data according to any one of claims 1 to 6.

14. A computer-readable storage medium, characterized in that A program or instruction is stored, and when a computer reads and executes the program or instruction, the computer executes the anisotropy correction method for array induction logging data according to any one of claims 1 to 6.

Images