CN103277084A - Horizontal well multi-parameter estimation method based on conducting probe array sensor - Google Patents

Abstract

The invention relates to a horizontal well multi-parameter estimation method based on a conducting probe array sensor. The horizontal well multi-parameter estimation method based on the conducting probe array sensor can be used for dynamic monitoring of a low-yield-liquid horizontal well yield section and is characterized by including the following steps that before the conducting probe array sensor is utilized to measure multiple parameters in the horizontal well, firstly the conducting probe array sensor is calibrated, calibrated data are acquired, and a calibrated database is established; when parameter estimation is conducted, distance matrixes are calculated according to actual response, acquired on the spot, of the conducting probe array sensor; lastly a matching vector in the calibrated database is found out according to the position of the least element in the distance matrixes, and then an azimuth angle of the conducting probe array sensor, aqueous phase electrical conductivity and oil-water interface height in the horizontal well are estimated. The horizontal well multi-parameter estimation method based on the conducting probe array sensor has the advantages that when the conducting probe array sensor is utilized to conduct horizontal well yield section well logging, the oil-water interface height and the aqueous phase electrical conductivity in the horizontal well and the azimuth angle of the conducting probe array sensor can be simultaneously measured in an online mode, a measuring lower limit is low, accuracy is high, and reliability is good.

Description

A kind of horizontal well multi-parameter method of estimation based on the conducting probe sensor array

Technical field

The present invention relates to a kind of horizontal well multi-parameter method of estimation based on the conducting probe sensor array, can be used for a plurality of parameters in the dynamic monitoring of low production liquid horizontal well production profile are accurately estimated.

Background technology

Production logging comprises injection well and the production well log technology after the completion, its objective is that the production of estimating oil well itself and oil reservoir is dynamic.In recent years, the task of production logging then expands to the whole exploitation course from each stage at drilling well initial stage to oil well.Different according to measuring object and range of application, the production logging method is broadly divided into three classes: produce dynamic logging, evaluation of producing well logging, engineering well logging.Horizontal well refers to that hole angle reaches or near 90 °, well bore creeps into the well of certain-length along horizontal direction.Horizontal well can significantly improve the gas in catchment building and the underground water unsaturation band, the contact area of underground oil and gas, improves the extraction efficiency of fluid effectively and improves recovery ratio.In horizontal well, the down-hole mixed flow is known from experience because density variation produces gravitational differentiation, has layered effect when making fluid layering mobile or mobile.Horizontal well is the raising oil well output that develops rapidly of China twentieth century end and the technology of development benefit.Common water horizontal well is 2.5～3.5 times of same well depth straight well cost, but its oil production can increase by 3～5 times, so horizontal well technology has obtained the generally attention of oil circle.Horizontal well drilling technology is gradually improved, but is subjected to the logger difficulty of transportation, the influence of factors such as flow pattern complexity, and horizontal well production logging technical development is slow.Domestic many oil fields all successful Application horizontal well, obtained favorable economic benefit, but long-term exploitation makes the production status instability, behind the part well water breakthrough moisture rising very fast, oil production obviously descends.Address these problems, will make the geologist understand the downhole production situation, so that it determines better constructure scheme, just production test has proposed more urgent requirement to horizontal well for this.

Summary of the invention

A kind of horizontal well multi-parameter method of estimation based on the conducting probe sensor array is characterized in that comprising the following steps:

(1) before use conducting probe sensor array is measured a plurality of parameters in the horizontal well, at first the conducting probe sensor array is demarcated, concrete steps are as follows:

A. before demarcating, on the azimuthal value of conducting probe sensor array zone, uniformly-spaced get m point, be designated as θ respectively ₁, θ ₂..., θ _m, on the value zone of oil-water interfaces height, uniformly-spaced get n point, be designated as h respectively ₁, h ₂..., h _n

B. in calibration process, the azimuth of conducting probe sensor array is adjusted into θ _i, be h with the oil-water interfaces Height Adjustment _j, i=1 wherein, 2 ..., m, j=1,2 ..., n measures the response of conducting probe sensor array this moment, is designated as

In calibration process, keeping the water electrical conductivity is constant, is designated as σ ^c, execution in step B repeatedly, the response of measuring conducting probe sensor array under every kind of condition obtains the nominal data storehouse, is expressed in matrix as

$G^c=\left[\begin{array}{cccc}{\mathrm{<figure-callout\; id="1"\; label="G"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">G</figure-callout>}}_{\mathrm{1,1}}^{\mathrm{<figure-callout\; id="1"\; label="c\; G"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">c</figure-callout>}}& G_{}^{}{}_{\mathrm{1,2}}^c& \&CenterDot;\&CenterDot;\&CenterDot;& {\mathrm{<figure-callout\; id="1"\; label="G"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">G</figure-callout>}}_{1,\mathrm{<figure-callout\; id="2"\; label="n\; c\; G"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">n</figure-callout>}}^c& \\ G_{}^{}{}_{\mathrm{2,1}}^{\mathrm{<figure-callout\; id="2"\; label="c\; G"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">c</figure-callout>}}& G_{}^{}{}_{\mathrm{2,2}}^c& \&CenterDot;\&CenterDot;\&CenterDot;& {\mathrm{<figure-callout\; id="2"\; label="G"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">G</figure-callout>}}_{2,n}^c\\ \&CenterDot;\&CenterDot;\&CenterDot;& \&CenterDot;\&CenterDot;\&CenterDot;& \&CenterDot;\&CenterDot;\&CenterDot;& \&CenterDot;\&CenterDot;\&CenterDot;\\ G_{m,1}^c& G_{m,2}^c& \&CenterDot;\&CenterDot;\&CenterDot;& G_{m,n}^c\end{array}\right];$

(2) oil-water interfaces height, conducting probe sensor array azimuth and water electrical conductivity are carried out parameter Estimation, concrete steps are as follows:

A. obtain the real response of conducting probe sensor array in measure field, be designated as G ^m

B. calculate distance matrix D, be expressed as

$D=\left[\begin{array}{cccc}{\mathrm{<figure-callout\; id="1"\; label="d"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">d</figure-callout>}}_{\mathrm{1,1}}& {\mathrm{<figure-callout\; id="1"\; label="d"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">\; <figure-callout\; id="2"\; label="d"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">d</figure-callout>\; </figure-callout>}}_{\mathrm{1,2}}& \&CenterDot;\&CenterDot;\&CenterDot;& {\mathrm{<figure-callout\; id="1"\; label="d"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">d</figure-callout>}}_{1,\mathrm{<figure-callout\; id="2"\; label="n\; d"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">n</figure-callout>}}& \\ d_{}{}_{\mathrm{2,1}}& {\mathrm{<figure-callout\; id="2"\; label="d"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">\; <figure-callout\; id="2"\; label="d"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">d</figure-callout>\; </figure-callout>}}_{\mathrm{2,2}}& \&CenterDot;\&CenterDot;\&CenterDot;& {\mathrm{<figure-callout\; id="2"\; label="d"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">d</figure-callout>}}_{2,n}\\ \&CenterDot;\&CenterDot;\&CenterDot;& \&CenterDot;\&CenterDot;\&CenterDot;& \&CenterDot;\&CenterDot;\&CenterDot;& \&CenterDot;\&CenterDot;\&CenterDot;\\ d_{m,1}& d_{m,2}& \&CenterDot;\&CenterDot;\&CenterDot;& d_{m,n}\end{array}\right],$

Element d in the matrix D wherein _{I, j}Design formulas be:

$d_{i,j}={\left|\right|\frac{G_{i,j}^c}{{\left|\right|G_{i,j}^c\left|\right|}_{\&infin;}}-\frac{G^m}{{\left|\right|G^m\left|\right|}_{\&infin;}}\left|\right|}_2,$

C. determine the position of least member in the matrix D, be designated as (r, s), the row at least member place number among the r representing matrix D wherein, the row at least member place number among the s representing matrix D are then at nominal data storehouse G ^cIn find out that (r, s) locational element is defined as matching vector, is designated as

D. utilize formula

$\left\{\begin{array}{c}\widetilde{\mathrm{\&theta;}}={\mathrm{\&theta;}}_r\\ \widetilde{\mathrm{\&sigma;}}=\frac{{\left|\right|G_{r,s}^c\left|\right|}_1}{{\left|\right|G^m\left|\right|}_1}\&CenterDot;{\mathrm{\&sigma;}}^c\\ \tilde{h}=h_s\end{array}\right.$

Determine each estimates of parameters, wherein

With

Be respectively water electrical conductivity in conducting probe sensor array azimuth, the horizontal well and the estimated value of oil-water interfaces height.

The invention has the beneficial effects as follows when adopting the conducting probe sensor array to carry out the horizontal well production profile logging, can be to the azimuth on-line measurement simultaneously of the oil-water interfaces height in the well, water electrical conductivity and conducting probe sensor array, measurement lower limit is low, precision is high, good reliability.

Description of drawings

Fig. 1 is flow chart of the present invention;

Fig. 2 is the structure chart of the conducting probe sensor array that adopts in the specific embodiment, among the figure: conducting probe (201), support arm (202), bolster (203);

Fig. 3 is the measurement sectional view of conducting probe sensor array in the specific embodiment, among the figure: horizontal well casing pipe (301), back shaft (302), support arm (303), conducting probe (305).

The specific embodiment

Existing accompanying drawings the specific embodiment of the present invention.Fig. 1 is flow chart of the present invention, Fig. 2 is the structure chart of the conducting probe sensor array that adopts in the specific embodiment, the sensor of this structure is by 24 conducting probes (201) that are distributed on the inside loop and outside loop, be installed in respectively on the support arm (202) of 12 equal angles distributions, back shaft (203) is coaxial with horizontal well, is used for fixing and the support whole sensor.Fig. 3 has shown in working order down, the measurement cross section of conducting probe sensor array and the numbering of probe thereof, and wherein the horizontal well inwall is 125mm, the external diameter of back shaft (302) is 50mm.The azimuth of conducting probe sensor array is its angle that turns over counterclockwise around bolster, and specifies Fig. 3 to show situation when the conducting probe sensor array is 15 ° at azimuth; When measuring, sleeve pipe (301) and bolster (302) current potential equate, are considered to the ground of whole measuring system, and the response of conducting probe sensor array is the electric conductivity value between each conducting probe and the ground, and availability vector is expressed as G={g ₁, g ₂..., g ₂₄.

(1) before use conducting probe sensor array is measured a plurality of parameters in the horizontal well, at first the conducting probe sensor array is demarcated, concrete steps are as follows:

A. before demarcating, the azimuthal value of conducting probe sensor array zone [0,2 π) on uniformly-spaced get m point, m=120, so discrete interval in the present embodiment is 3 °, and is designated as θ according to this ₁, θ ₂..., θ ₁₂₀, on the value zone [0,125] of oil-water interfaces height, uniformly-spaced getting n point, n=51, so discrete interval in the present embodiment is 2.5mm, and is designated as h successively ₁, h ₂..., h _n

B. in calibration process, the azimuth of conducting probe sensor array is adjusted into θ _i, be h with the oil-water interfaces Height Adjustment _j, i=1 wherein, 2 ..., 120, j=1,2 ..., 51, measure the response of conducting probe sensor array this moment, be designated as

In calibration process, keeping the water electrical conductivity is constant, is designated as σ ^c, σ after measured ^c=0.078S/m, execution in step B repeatedly, the response of measuring conducting probe sensor array under every kind of condition obtains the nominal data storehouse, is expressed in matrix as

$G^c=\left[\begin{array}{cccc}{\mathrm{<figure-callout\; id="1"\; label="G"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">G</figure-callout>}}_{\mathrm{1,1}}^{\mathrm{<figure-callout\; id="1"\; label="c\; G"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">c</figure-callout>}}& G_{}^{}{}_{\mathrm{1,2}}^c& \&CenterDot;\&CenterDot;\&CenterDot;& {\mathrm{<figure-callout\; id="1"\; label="G"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">G</figure-callout>}}_{\mathrm{1,51}}^{\mathrm{<figure-callout\; id="2"\; label="c\; G"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">c</figure-callout>}}& \\ G_{}^{}{}_{\mathrm{2,1}}^{\mathrm{<figure-callout\; id="2"\; label="c\; G"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">c</figure-callout>}}& G_{}^{}{}_{\mathrm{2,2}}^c& \&CenterDot;\&CenterDot;\&CenterDot;& G_{\mathrm{2,51}}^c\\ \&CenterDot;\&CenterDot;\&CenterDot;& \&CenterDot;\&CenterDot;\&CenterDot;& \&CenterDot;\&CenterDot;\&CenterDot;& \&CenterDot;\&CenterDot;\&CenterDot;\\ G_{\mathrm{120,1}}^c& G_{\mathrm{120,1}}^c& \&CenterDot;\&CenterDot;\&CenterDot;& G_{\mathrm{120,51}}^c\end{array}\right],$

Because m=120, n=51 in the present embodiment, so nominal data storehouse G ^cBe 120 row, 51 row, and the element in the matrix be vector, vectorial figure place is 24, represents under the combination of different orientations and oil-water interfaces height the response of conducting probe array.

(2) oil-water interfaces height, conducting probe sensor array azimuth and water electrical conductivity are carried out parameter Estimation, concrete steps are as follows:

A. obtain the real response of conducting probe sensor array in measure field, be designated as G ^mFor checking effect of the present invention, in the present embodiment, the real response G of conducting probe sensor array ^mInstitute's respective sensor azimuth angle theta, oil-water interfaces height h and water conductivity's actual value is known, is respectively θ=123.6 °, h=62.8mm, σ=0.897S/m.

B. calculate distance matrix D, be expressed as

$D=\left[\begin{array}{cccc}{\mathrm{<figure-callout\; id="1"\; label="d"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">d</figure-callout>}}_{\mathrm{1,1}}& {\mathrm{<figure-callout\; id="1"\; label="d"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">\; <figure-callout\; id="2"\; label="d"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">d</figure-callout>\; </figure-callout>}}_{\mathrm{1,2}}& \&CenterDot;\&CenterDot;\&CenterDot;& {\mathrm{<figure-callout\; id="1"\; label="d"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">d</figure-callout>}}_{1,51}\\ {\mathrm{<figure-callout\; id="2"\; label="d"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">d</figure-callout>}}_{\mathrm{2,1}}& {\mathrm{<figure-callout\; id="2"\; label="d"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">\; <figure-callout\; id="2"\; label="d"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">d</figure-callout>\; </figure-callout>}}_{\mathrm{2,2}}& \&CenterDot;\&CenterDot;\&CenterDot;& {\mathrm{<figure-callout\; id="2"\; label="d"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">d</figure-callout>}}_{2,51}\\ \&CenterDot;\&CenterDot;\&CenterDot;& \&CenterDot;\&CenterDot;\&CenterDot;& \&CenterDot;\&CenterDot;\&CenterDot;& \&CenterDot;\&CenterDot;\&CenterDot;\\ d_{120,1}& {\mathrm{<figure-callout\; id="2"\; label="d"\; filenames="HDA00003234737700022.png"\; state="\{\{state\}\}">d</figure-callout>}}_{120,2}& \&CenterDot;\&CenterDot;\&CenterDot;& d_{120,51}\end{array}\right],$

Element d in the matrix D wherein _{I, j}Design formulas be:

$d_{i,j}={\left|\right|\frac{G_{i,j}^c}{{\left|\right|G_{i,j}^c\left|\right|}_{\&infin;}}-\frac{G^m}{{\left|\right|G^m\left|\right|}_{\&infin;}}\left|\right|}_2,$

Because m=120, n=51 in the present embodiment, so distance matrix D is 120 row, 51 row, and the element in the matrix is scalar, the real response G of representative sensor ^mWith nominal data storehouse matrix G ^cThe Euclidean distance of middle corresponding element.

C. determine the position of least member in the matrix D, be designated as (r, s), the row at least member place number among the r representing matrix D wherein, the row at least member place number among the s representing matrix D are then at nominal data storehouse G ^cIn find out that (r, s) locational element is defined as matching vector, is designated as

R=42 in the present embodiment, s=51, so matching vector is

D. utilize formula

$\left\{\begin{array}{c}\widetilde{\mathrm{\&theta;}}={\mathrm{\&theta;}}_r\\ \widetilde{\mathrm{\&sigma;}}=\frac{{\left|\right|G_{r,s}^c\left|\right|}_1}{{\left|\right|G^m\left|\right|}_1}\&CenterDot;{\mathrm{\&sigma;}}^c\\ \tilde{h}=h_s\end{array}\right.$

Determine each estimates of parameters, wherein With

Be respectively water electrical conductivity in conducting probe sensor array azimuth, the horizontal well and the estimated value of oil-water interfaces height.M=120, n=51, r=42, s=51, matching vector are in the present embodiment

θ then _r=123 °, h _s=62.5mm, therefore

Through calculating

With actual value θ=123.6 °, h=62.8mm, σ=0.897S/m relatively, illustrate that this method is feasible.

More than to the description of the present invention and embodiment thereof, be not limited thereto, only be one of embodiments of the present invention shown in the accompanying drawing.Under the situation that does not break away from the invention aim, without designing and the similar structure of this technical scheme or embodiment, all belong to protection domain of the present invention with creating.

Claims (1)

1. the horizontal well multi-parameter method of estimation based on the conducting probe sensor array is characterized in that comprising the following steps:

(1) before use conducting probe sensor array is measured a plurality of parameters in the horizontal well, at first the conducting probe sensor array is demarcated, concrete steps are as follows:

A. before demarcating, on the azimuthal value of conducting probe sensor array zone, uniformly-spaced get m point, be designated as θ respectively ₁, θ ₂..., θ _m, on the value zone of oil-water interfaces height, uniformly-spaced get n point, be designated as h respectively ₁, h ₂..., h _n

B. in calibration process, the azimuth of conducting probe sensor array is adjusted into θ _i, be h with the oil-water interfaces Height Adjustment _j, i=1 wherein, 2 ..., m, j=1,2 ..., n measures the response of conducting probe sensor array this moment, is designated as

In calibration process, keeping the water electrical conductivity is constant, is designated as σ ^c, execution in step B repeatedly, the response of measuring conducting probe sensor array under every kind of condition obtains the nominal data storehouse, is expressed in matrix as

$G^c=\left[\begin{array}{cccc}{G}_{\mathrm{1,1}}^c& G_{\mathrm{1,2}}^c& \&CenterDot;\&CenterDot;\&CenterDot;& G_{1,n}^c\\ G_{\mathrm{2,1}}^c& G_{\mathrm{2,2}}^c& \&CenterDot;\&CenterDot;\&CenterDot;& G_{2,n}^c\\ \&CenterDot;\&CenterDot;\&CenterDot;& \&CenterDot;\&CenterDot;\&CenterDot;& \&CenterDot;\&CenterDot;\&CenterDot;& \&CenterDot;\&CenterDot;\&CenterDot;\\ G_{m,1}^c& G_{m,2}^c& \&CenterDot;\&CenterDot;\&CenterDot;& G_{m,n}^c\end{array}\right];$

(2) oil-water interfaces height, conducting probe sensor array azimuth and water electrical conductivity are carried out parameter Estimation, concrete steps are as follows:

A. obtain the real response of conducting probe sensor array in measure field, be designated as G ^m

B. calculate distance matrix D, be expressed as

$D=\left[\begin{array}{cccc}{d}_{\mathrm{1,1}}& d_{\mathrm{1,2}}& \&CenterDot;\&CenterDot;\&CenterDot;& d_{1,n}\\ d_{\mathrm{2,1}}& d_{\mathrm{2,2}}& \&CenterDot;\&CenterDot;\&CenterDot;& d_{2,n}\\ \&CenterDot;\&CenterDot;\&CenterDot;& \&CenterDot;\&CenterDot;\&CenterDot;& \&CenterDot;\&CenterDot;\&CenterDot;& \&CenterDot;\&CenterDot;\&CenterDot;\\ d_{m,1}& d_{m,2}& \&CenterDot;\&CenterDot;\&CenterDot;& d_{m,n}\end{array}\right],$

Element d in the matrix D wherein _{I, j}Design formulas be:

$d_{i,j}={\left|\right|\frac{G_{i,j}^c}{{\left|\right|G_{i,j}^c\left|\right|}_{\&infin;}}-\frac{G^m}{{\left|\right|G^m\left|\right|}_{\&infin;}}\left|\right|}_2,$

C. determine the position of least member in the matrix D, be designated as (r, s), the row at least member place number among the r representing matrix D wherein, the row at least member place number among the s representing matrix D are then at nominal data storehouse G ^cIn find out that (r, s) locational element is defined as matching vector, is designated as

D. utilize formula

$\left\{\begin{array}{c}\widetilde{\mathrm{\&theta;}}={\mathrm{\&theta;}}_r\\ \widetilde{\mathrm{\&sigma;}}=\frac{{\left|\right|G_{r,s}^c\left|\right|}_1}{{\left|\right|G^m\left|\right|}_1}\&CenterDot;{\mathrm{\&sigma;}}^c\\ \tilde{h}=h_s\end{array}\right.$

Determine each estimates of parameters, wherein

With Be respectively water electrical conductivity in conducting probe sensor array azimuth, the horizontal well and the estimated value of oil-water interfaces height.

Images