CN101776771B - Lithological data acquisition and treatment method - Google Patents

Abstract

The invention relates to a lithological data acquisition and treatment method for geophysical exploration and development, relates to the field of the data acquisition and test analysis of formation lithological parameters, and is suitable for identifying an underground water layer or an oil layer or the oil layer containing movable water based on the lithological data. The method comprises the following steps: using geological exploration and development equipment to acquire the formation lithological data; analyzing the acquired formation lithological data by using analytical equipment to obtain the basic parameters needed for calculating a sandstone apparent formation resistivity, a formation fluid replacement rate and slurry filtrate distribution coefficients, and judging whether the sandstone formation is the water layer, the oil layer or the oil layer containing movable water by comparing the numerical relationship among the depth lateral apparent resistivity of the polluted sandstone formation, the true resistivity and the actually measured resistivity. The invention avoids the waste caused by mistakenly perforating the water layer by analyzing the pollution level of the sandstone formation around wellbore when the oil well is subjected to development perforation, reduces the erroneous judgment of the abnormal layers, effectively improves the single-well production efficiency of oil reservoirs, and is in favor of economically and highly efficiently developing the oil reservoirs.

Description

A kind of lithological data acquisition and treatment method

Technical field

The present invention relates to the geophysical exploration development method, be specially a kind of lithological data acquisition and treatment method, belong to well logging magnanimity information comprehensive interpretation and evaluation category and relate to the data acquisition test analysis field of formation lithology parameter, be particularly useful for discerning water table or oil reservoir or containing movable water oil reservoir based on lithology data.

Background technology

Estimate the sandstone oil reservoir oil-water-layer rapidly and accurately is the technical barrier that oil field prospecting exploitation scientific and technical personnel will solve always.Oil-water-layer evaluation at present mainly contains five class methods: be respectively the method for coring, method of testing, well logging method, chart method and overall approach.

As everyone knows, utilizing and to core or measuring technology is that to estimate oil-water-layer the most direct, also is effectively evaluating technology, but can not a mouthful mouthful well all cores or each payzone is all tested, and this is both uneconomical also unrealistic.Three kinds of methods in back particularly chart method are estimated the oil-water-layer economical and efficient, use more general.It is to be noted that well-log information is a formation information of using all can not cast aside when any method is estimated oil-water-layer the most general, most worthy.

Chart method normally adopts logging technology in conjunction with coring and testing achievement, formulate the oil-water-layer criterion of identification according to the Archie equation, but because drilling period is long, mud is seriously polluted to stratum (comprising oil reservoir and water layer), destroyed the original electrical property feature in stratum, simultaneously because the mud property disunity between well and the well, thereby can't formulate unified oil-water-layer criterion of identification, cause the oil-water-layer discrimination low, bring very big loss to the oil field exploitation,, carried out depth side direction resistivity logging and studied by the resistivity assessment technique behind the mud contamination because the mud contamination adverse effect that evaluation causes to oil-water-layer for solving.

Because oil saturation changes, the conductive fluid character in the stratum is changed behind the mud intrusion virgin zone.The variation of describing formation fluid property is an objective key of asking for formation resistivity and identification characteristic of fluid behind the solution mud intrusion virgin zone how objective and accurately.

Large quantities of experts and scholars have done number of research projects in this regard, but what all come with some shortcomings, as the patent No. is CN1243958A " properties of fluid in bearing stratum measuring method and equipment ", but this patent is the recombination current of the two or more frequencies sent by recording unit produces respective frequencies in reservoir a measuring-signal, by contrast to tracer signal, think that the little well section reservoir fluid of amplitude difference ratio is a water, the well section reservoir fluid that amplitude difference ratio is big is oil or gas.This patent is differentiated the meeting appearance than large deviation for contaminated oil reservoir or water layer, when the original state fluid in the sandstone formation during by the part or all of displacement of mud, the fluid properties of virgin zone (oil reservoir or water layer) all changes, and uses this patented technology and has oil reservoir or water layer misjudgment phenomenon.

For another example, the patent No. be CN 101363315A's " resistivity on quantification thin layer stratum and the method for hydrocarbon saturation ", this patent is the irreducible water cumulative volume according to each layer of the irreducible water saturation of each layer and total porosity and estimation, calculate stratum horizontal resistivity and vertical resistivity, the horizontal resistivity and the vertical resistivity of estimated value and mensuration are compared, adjust the irreducible water saturation of each layer of being estimated and estimate this value repeatedly, difference between estimated value and the vertical resistivity measured drops to selected valve value, by the irreducible water saturation estimation hydrocarbon volume through each layer of adjusting.This patent adopts adjusted repeatedly irreducible water saturation to realize estimation resistivity and the consistance of surveying resistivity, and with adjusted irreducible water saturation estimation hydrocarbonaceous volume, this patent adopts the hydrocarbon volume confidence level of adjusted irreducible water saturation estimation to reduce the formation resistivity that causes raises or reduction is all summed up in the point that on the irreducible water saturation because mud is invaded thereby directly cause.

Summary of the invention

In order to solve the problem that exists in the top background technology, the present invention proposes a kind of lithological data acquisition and treatment method, belong to well logging magnanimity information comprehensive interpretation and evaluation category and relate to the data acquisition test analysis field of formation lithology parameter, fluid properties in formation lithology or the rock stratum is discerned in its calculating based on resistivity.

According to technical solutions according to the invention, a kind of lithological data acquisition and treatment method is provided, it may further comprise the steps:

1) utilizes geologic prospecting equipment, can gather the formation lithology data;

2) based on the formation lithology data that collect, by sandstone total porosity, irreducible water saturation, local water density and the in-place oil density in the lithology data analytical equipment acquisition stratum;

3), obtain to comprise movable formation water resistivity R by the lithology data analytical equipment based on the formation lithology data that collect _WfWith stratum irreducible water resistivity R _WiFormation water resistivity, mud and mud filtrate resistivity and sandstone formation true resistivity;

4) according in step 2) in sandstone total porosity, irreducible water saturation, local water density and the in-place oil density of acquisition, and, obtain resident fluid replacement rate and mud filtrate partition factor according to the formation water resistivity that in step 3), obtains, mud and mud filtrate resistivity and sandstone formation true resistivity; Wherein the resident fluid replacement rate is meant in the specified scope of wellbore (investigation depth that refers to the side direction resistivity logging), invade intrusion liquid in the perviousness sandstone formation long-pending with the original state sandstone formation in the ratio of movable fluid volume, invade long-pending mud and the mud filtrate volume sum of being meant of liquid; The mud filtrate partition factor is meant the mud filtrate and the long-pending ratio of intrusion liquid of invading in the perviousness sandstone formation.

Resident fluid replacement rate expression formula:

$u_f=\frac{V_F}{V}---\left(1\right)$

V _F＝V _mud+W _mf (2)

In the formula: V _FMud and the mud filtrate cumulative volume in the stratum (side direction resistivity logging investigation depth), m are specified in-intrusion ³(cubic meter); V _MudThe mud volume in the stratum (side direction resistivity logging investigation depth), m are specified in-intrusion ³(cubic meter); V _MfThe mud filtrate volume in the stratum (side direction resistivity logging investigation depth), m are specified in-intrusion ³(cubic meter); V-specifies movable fluid volume in the virgin zone (the side direction resistivity logging is surveyed dark), m ³(cubic meter); u _f-resident fluid replacement rate, 0≤u _f≤ 1, zero dimension;

Mud filtrate partition factor expression formula:

$v_f=\frac{V_{\mathrm{mf}}}{V_{\mathrm{mud}}+V_{\mathrm{mf}}}---\left(3\right)$

Mud filtrate partition factor span 0≤v in the formula _f＜1.

5) base area layer fluid replacement rate and mud filtrate partition factor can determine the attribute of sandstone formation, and promptly whether the stratum is pure water layer or net pay zone or contains movable water oil reservoir.

Wherein, after mud and mud filtrate were invaded the pure water layer, with all or part of displacement of movable water in the hole, the mobile water saturation in the stratum was S at this moment _Wc=1-S _Wi, the apparent resistivity explanation formula of depth side direction resistivity logging in the pure water layer is:

$R_{\mathrm{LLso}}=\frac{1}{\frac{{\mathrm{\φ}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+{\mathrm{\φ}}_t(\frac{u_{\mathrm{fs}}}{R_{\mathrm{mfs}}}+\frac{1-u_{\mathrm{fs}}}{R_{\mathrm{wf}}})(1-S_{\mathrm{wi}})}---\left(4\right)$

$R_{\mathrm{mfs}}=\frac{1}{\frac{v_f}{R_{\mathrm{mf}}}+\frac{1-v_f}{R_m}}---\left(5\right)$

$R_{\mathrm{LLdo}}=\frac{1}{\frac{{\mathrm{\φ}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+{\mathrm{\φ}}_t(\frac{u_{\mathrm{fd}}}{R_{\mathrm{mf}}}+\frac{1-u_{\mathrm{fd}}}{R_{\mathrm{wf}}})(1-S_{\mathrm{wi}})}---\left(6\right)$

In the formula: R _LLsoThe shallow side direction apparent resistivity of-100% water bearing sand, Ω m; R _LLdo-100% water bearing sand deep lateral apparent resistivity, Ω m; R _Wi-stratum irreducible water resistivity, Ω m; R _Wf-movable formation water resistivity, Ω m; R _MfMud filtrate resistivity under the-formation temperature, Ω m; R _Mfs-comprehensive the resistivity of intrusion liquid of (investigation depth 0.35m) in shallow side direction resistivity logging investigative range, Ω m; S _Wi-sandstone irreducible water saturation, 0＜S _Wi≤ 1; u _Fs-resident fluid the replacement rate of (investigation depth 0.35m) in shallow side direction resistivity logging investigative range, 0≤u _Fs≤ 1; u _Fd-resident fluid the replacement rate of (investigation depth 1.15m) in dark side direction resistivity logging investigative range, 0≤u _Fd≤ 1; v _f-mud filtrate partition factor, 0≤v _f＜1; φ _t-sandstone total porosity, 0＜φ _t≤ 47.64%.

When adopting the fresh water mud drilling well, invade liquid resistivity greater than formation water resistivity, it is " increasing resistance " (increase progressively or increase) or " partly increasing resistance " mud invasion character that depth side direction apparent resistivity logging shows water layer, i.e. R _o＜R _LLdo＜R _LLso(accompanying drawing 1 is invaded the first synoptic diagram A of the trapezoidal trend of water layer for fresh water mud), or R _o＜R _LLso＜R _LLdo(accompanying drawing 2 is invaded the second synoptic diagram B of the trapezoidal trend of water layer for fresh water mud), or R _LLso＜R _o＜R _LLdo(accompanying drawing 3 is invaded the 3rd synoptic diagram C of the trapezoidal trend of water layer for fresh water mud); When adopting the salty mud drilling well, invade liquid resistivity less than formation water resistivity, it is " drag reduction " mud invasion character that depth side direction resistivity logging shows water layer, i.e. R _o＞R _LLdo＞R _LLso(accompanying drawing 4 is invaded the 4th synoptic diagram D of the trapezoidal trend of water layer for salty mud) also will satisfy when satisfying one of above-mentioned four kinds of situations

Or

Numerical relation, the sandstone formation of need judging is the pure water layer.

R wherein _oBe sandstone water layer true resistivity; R _TsBe shallow side direction resistivity logging actual measurement resistivity; R _TdBe dark side direction resistivity logging actual measurement resistivity.

Wherein, behind mud and mud filtrate intrusion net pay zone, with the part or all of displacement of the oil in the hole, owing to no original movable water existence, i.e. S in the sandstone pores _Wf=0, obtain the apparent resistivity explanation formula of depth side direction apparent resistivity logging in sandstone formation and be:

$R_{\mathrm{LLs}}=\frac{1}{\frac{{\mathrm{\φ}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+\frac{{\mathrm{\φ}}_t{u}_{\mathrm{fs}}(1-S_{\mathrm{wi}})}{R_{\mathrm{mfs}}}}---\left(7\right)$

$R_{\mathrm{LLd}}=\frac{1}{\frac{{\mathrm{\φ}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+\frac{{\mathrm{\φ}}_t{u}_{\mathrm{fd}}(1-S_{\mathrm{wi}})}{R_{\mathrm{mf}}}}---\left(8\right)$

In the formula: R _LLsThe shallow side direction apparent resistivity of-sandstone oil reservoir, Ω m; R _LLd-sandstone oil reservoir deep lateral apparent resistivity, Ω m; R _Wi-stratum irreducible water resistivity, Ω m; R _Wf-movable formation water resistivity, Ω m; R _MfMud filtrate resistivity under the-formation temperature, Ω m; R _Mfs-comprehensive the resistivity of intrusion liquid of (investigation depth 0.35m) in shallow side direction resistivity logging investigative range, Ω m; S _Wi-sandstone irreducible water saturation, 0＜S _Wi≤ 1; u _Fs-resident fluid the replacement rate of (investigation depth 0.35m) in shallow side direction resistivity logging investigative range, 0≤u _Fs≤ 1; u _Fd-resident fluid the replacement rate of (investigation depth 1.15m) in dark side direction resistivity logging investigative range, 0≤u _Fd≤ 1; v _f-mud filtrate partition factor, 0≤v _f＜1; φ _t-sandstone total porosity, 0＜φ _t≤ 47.64%.

According to R _LLs, R _LLd, R _LLso, R _LLdo, R _t, R _TsAnd R _TdBetween numerical relation, whether can determine sandstone formation is net pay zone.

Further, when adopting the fresh water mud drilling well, invade liquid resistivity greater than formation water resistivity, it is " drag reduction " (successively decreasing) mud invasion character that depth side direction resistivity logging shows the sandstone oil reservoir, promptly has R _t＞R _LLd＞R _LLsNumerical relation the time, accompanying drawing 5 is that fresh water mud invades the 5th synoptic diagram E of the trapezoidal trend of oil reservoir, the sandstone formation that needs to judge is a net pay zone; When adopting the salty mud drilling well, invade liquid resistivity less than formation water resistivity, it is " drag reduction " (speed subtracts) mud invasion character that depth side direction resistivity logging shows oil reservoir, promptly has R _t＞R _LLd＞＞R _LLs(accompanying drawing 6 is invaded the 6th synoptic diagram F of the trapezoidal trend of oil reservoir for salty mud) when satisfying one of above-mentioned two kinds of situations, also will satisfy

Or

Numerical relation, the sandstone formation of need judging is a net pay zone.

R wherein _tBe sandstone oil reservoir true resistivity; R _TsBe shallow side direction resistivity logging actual measurement resistivity; R _TdBe dark side direction resistivity logging actual measurement resistivity.

Additionally, for containing movable water oil reservoir, the apparent resistivity explanation formula of depth side direction resistivity logging in sandstone formation is:

$R_{\mathrm{LLs}}^{\′}=\frac{1}{\frac{{\mathrm{\φ}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+\frac{{\mathrm{\φ}}_t{u}_{\mathrm{fs}}(1-S_{\mathrm{wi}})}{R_{\mathrm{mfs}}}+\frac{{\mathrm{\φ}}_t(1-u_{\mathrm{fs}})S_{\mathrm{wf}}}{R_{\mathrm{wf}}}}---\left(9\right)$

$R_{\mathrm{mfs}}=\frac{1}{\frac{v_f}{R_{\mathrm{mf}}}+\frac{1-v_f}{R_m}}---\left(10\right)$

$R_{\mathrm{LLd}}^{\′}=\frac{1}{\frac{{\mathrm{\φ}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+\frac{{\mathrm{\φ}}_t{u}_{\mathrm{fd}}(1-S_{\mathrm{wi}})}{R_{\mathrm{mf}}}+\frac{{\mathrm{\φ}}_t(1-u_{\mathrm{fd}})S_{\mathrm{wf}}}{R_{\mathrm{wf}}}}---\left(11\right)$

In the formula: R ' _LLs-contain the shallow side direction apparent resistivity of movable water oil reservoir, Ω m; R ' _LLd-contain movable water oil reservoir deep lateral apparent resistivity, Ω m; S _Wf-sandstone formation mobile water saturation, 0≤S _Wf＜1; R _Wi-stratum irreducible water resistivity, Ω m; R _Wf-movable formation water resistivity, Ω m; R _MfMud filtrate resistivity under the-formation temperature, Ω m; R _Mfs-comprehensive the resistivity of intrusion liquid of (investigation depth 0.35m) in shallow side direction resistivity logging investigative range, Ω m; S _Wi-sandstone irreducible water saturation, 0＜S _Wi≤ 1; u _Fs-resident fluid the replacement rate of (investigation depth 0.35m) in shallow side direction resistivity logging investigative range, 0≤u _Fs≤ 1; u _Fd-resident fluid the replacement rate of (investigation depth 1.15m) in dark side direction resistivity logging investigative range, 0≤u _Fd≤ 1; v _f-mud filtrate partition factor, 0≤v _f＜1; φ _t-sandstone total porosity, 0＜φ _t≤ 47.64%.

According to R _LLd, R _LLs, R _LLso, R _LLdo, R ' _LLs, R ' _LLd, R _t, R _TsAnd R _TdBetween numerical relation, whether can determine sandstone formation is to contain movable water oil reservoir.

Further, work as R _LLd＞R ' _LLd＞R _LLdoOr R _LLs＞R ' _LLs＞R _LLso, exist simultaneously

Or Numerical relation the time, the sandstone formation that judge is to contain movable water oil reservoir.

Preferably, obtain the sandstone total porosity according to core analysis factor of porosity and actual measurement interval transit time.

Preferably, adopt the pressure mercury data that not influenced by borehole environment to determine irreducible water saturation.

Preferably, local water density is the function of reservoir water salinity, formation temperature and reservoir pressure.

Preferably, movable formation water resistivity R _WfAfter can being converted to equivalent N aCl salinity by reservoir water salinity (SC), concern that by the total mineralization and the equivalent coefficient of corresponding ion concentration plate looks into and get (synoptic diagram that accompanying drawing 7 concerns for local water equivalent N aCl solution resistance rate and temperature intersection for the equivalent coefficient graph of a relation of total mineralization and corresponding ion concentration, accompanying drawing 8).

Preferably, the relation between pure water layer irreducible water saturation and the mobile water saturation:

1-S _wi＝S _wc

In the formula: S _WcMobile water saturation in-100% water bearing sand, 0≤S _Wc＜1; S _Wi-sandstone irreducible water saturation, 0＜S _Wi≤ 1.

Use and of the present inventionly discern underground pure water layer or net pay zone or contain movable water oil reservoir based on a kind of lithological data acquisition and treatment method, and perforated interval selected, reject the pure water layer, select net pay zone that the well section that will produce is implemented perforation, operation can improve the individual well production efficiency of oil reservoir effectively like this, reduce cost of development, help the exploitation of oil reservoir economical and efficient ground.The present invention has assessed the pollution level of mud to the wellbore sandstone formation, has avoided oil well when dropping into the exploitation perforation, has penetrated pure water layer or the high waste that oil reservoir caused that contains movable water unintentionally by mistake, has reduced unusual erroneous judgement.

Brief Description Of Drawings

Fig. 1 invades the first synoptic diagram A of the trapezoidal trend of water layer for fresh water mud;

Fig. 2 invades the second synoptic diagram B of the trapezoidal trend of water layer for fresh water mud;

Fig. 3 invades the 3rd synoptic diagram C of the trapezoidal trend of water layer for fresh water mud;

Fig. 4 invades the 4th synoptic diagram D of the trapezoidal trend of water layer for salty mud;

Fig. 5 invades the 5th synoptic diagram E of the trapezoidal trend of oil reservoir for fresh water mud;

Fig. 6 invades the 6th synoptic diagram F of the trapezoidal trend of oil reservoir for salty mud;

Fig. 7 is a total mineralization and the equivalent coefficient graph of a relation of corresponding ion concentration;

Fig. 8 is the synoptic diagram of local water equivalent N aCl solution resistance rate and temperature intersection relation;

Fig. 9 is the synoptic diagram of oil saturation and oil reservoir height relationships;

Figure 10 is the synoptic diagram of original solution gas-oil ratio and in-place oil density relationship;

Figure 11 is for using the certain herbaceous plants with big flowers east 103 well interpretation results synoptic diagram of technical solution of the present invention;

Figure 12 is for using the brocade 315 well interpretation results synoptic diagram of technical solution of the present invention.

Figure 13 is for using the certain herbaceous plants with big flowers east 102 well interpretation results synoptic diagram of technical solution of the present invention.

Embodiment

Discern underground pure water layer, net pay zone or contain movable water oil reservoir according to lithological data acquisition and treatment method of the present invention, have the advantage that identification is accurate, False Rate is low.This lithological data acquisition and treatment method mainly obtains based on following fact research.

In side direction resistivity logging investigative range, mud and mud filtrate are invaded the original state sandstone formation, and with local water in the hole or the part or all of displacement of oil, local water character can change because of the intrusion of mud and mud filtrate in the hole.For after accurate description mud invades sandstone formation, the variation of original state fluid in the sandstone formation the present invention proposes the thought of " resident fluid replacement rate " and " mud filtrate partition factor "." resident fluid replacement rate " is meant within the specific limits (investigation depth that refers to the side direction resistivity logging), invade intrusion liquid in the perviousness sandstone formation long-pending with the original state sandstone formation in the ratio of movable fluid volume, invade the long-pending mud and the mud filtrate volume sum of invading in the perviousness sandstone formation that be meant of liquid; " mud filtrate partition factor " is meant the mud filtrate and the long-pending ratio of intrusion liquid of invading in the perviousness sandstone formation.

The present invention adopts following technical scheme to discern underground pure water layer, net pay zone or contains movable water oil reservoir, promptly a kind of lithological data acquisition and treatment method, and it may further comprise the steps:

1) utilizes geologic prospecting equipment, can gather the formation lithology data;

2) based on the formation lithology data that collect, by sandstone total porosity, irreducible water saturation, local water density and the in-place oil density in the lithology data analytical equipment acquisition stratum;

3), obtain to comprise movable formation water resistivity R by the lithology data analytical equipment based on the formation lithology data that collect _WfWith stratum irreducible water resistivity R _WiFormation water resistivity, mud and mud filtrate resistivity and sandstone formation true resistivity;

4) according to the sandstone total porosity, irreducible water saturation, local water density and the in-place oil density that in step 1), obtain, and according in step 2) in formation water resistivity, mud and the mud filtrate resistivity and the sandstone formation true resistivity that obtain, obtain resident fluid replacement rate and mud filtrate partition factor; Wherein the resident fluid replacement rate is meant in the specified scope of wellbore (investigation depth that refers to the side direction resistivity logging), invade intrusion liquid in the perviousness sandstone formation long-pending with the original state sandstone formation in the ratio of movable fluid volume, invade long-pending mud and the mud filtrate volume sum of being meant of liquid; The mud filtrate partition factor is meant the mud filtrate and the long-pending ratio of intrusion liquid of invading in the sandstone formation.

Resident fluid replacement rate expression formula:

$u_f=\frac{V_F}{V}---\left(1\right)$

V _F＝V _mud+F _mf (2)

In the formula: V _FMud and the mud filtrate cumulative volume in the stratum (side direction resistivity logging investigation depth), m are specified in-intrusion ³(cubic meter); V _MudThe mud volume in the stratum (side direction resistivity logging investigation depth), m are specified in-intrusion ³(cubic meter); V _MgThe mud filtrate volume in the stratum (side direction resistivity logging investigation depth), m are specified in-intrusion ³(cubic meter); V-specifies movable fluid volume in the virgin zone (the side direction resistivity logging is surveyed dark), m ³(cubic meter); u _f-resident fluid replacement rate, 0≤u _f≤ 1, zero dimension.

Mud filtrate partition factor expression formula:

$v_f=\frac{V_{\mathrm{mf}}}{V_{\mathrm{mud}}+V_{\mathrm{mf}}}---\left(3\right)$

The limit span of mud filtrate partition factor is 0≤v in the formula _f＜1.

5) base area layer fluid replacement rate and mud filtrate partition factor can determine the attribute of sandstone formation, and promptly whether the stratum is pure water layer or net pay zone or contains movable water oil reservoir.

Wherein, after mud and mud filtrate were invaded water layer, with all or part of displacement of movable water in the hole, the mobile water saturation in the stratum was S at this moment _Wf=1-S _Wi, the apparent resistivity explanation formula of depth side direction resistivity logging in water layer is:

$R_{\mathrm{LLso}}=\frac{1}{\frac{{\mathrm{\φ}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+{\mathrm{\φ}}_t(\frac{u_{\mathrm{fs}}}{R_{\mathrm{mfs}}}+\frac{1-u_{\mathrm{fs}}}{R_{\mathrm{wf}}})(1-S_{\mathrm{wi}})}---\left(4\right)$

$R_{\mathrm{mfs}}=\frac{1}{\frac{v_f}{R_{\mathrm{mf}}}+\frac{1-v_f}{R_m}}---\left(5\right)$

$R_{\mathrm{LLdo}}=\frac{1}{\frac{{\mathrm{\φ}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+{\mathrm{\φ}}_t(\frac{u_{\mathrm{fd}}}{R_{\mathrm{mf}}}+\frac{1-u_{\mathrm{fd}}}{R_{\mathrm{wf}}})(1-S_{\mathrm{wi}})}---\left(6\right)$

In the formula: R _LLsoThe shallow side direction apparent resistivity of-100% water bearing sand, Ω m; R _LLdo-100% water bearing sand deep lateral apparent resistivity, Ω m; R _Wi-stratum irreducible water resistivity, Ω m; R _Wf-movable formation water resistivity, Ω m; R _MfMud filtrate resistivity under the-formation temperature, Ω m; R _Mfs-comprehensive the resistivity of intrusion liquid of (investigation depth 0.35m) in shallow side direction resistivity logging investigative range, Ω m; S _Wi-sandstone irreducible water saturation, 0＜S _Wi≤ 1; u _Fs-resident fluid the replacement rate of (investigation depth 0.35m) in shallow side direction resistivity logging investigative range, 0≤u _Fs≤ 1; u _Fd-resident fluid the replacement rate of (investigation depth 1.15m) in dark side direction resistivity logging investigative range, 0≤u _Fd≤ 1; v _f-mud filtrate partition factor, 0≤v _f＜1; φ _t-sandstone total porosity, 0＜φ _t≤ 47.64%.

When adopting the fresh water mud drilling well, invade liquid resistivity greater than formation water resistivity, it is " increasing resistance " (increase progressively or increase) or " partly increasing resistance " mud invasion character that depth side direction apparent resistivity logging shows water layer, i.e. R _o＜R _LLdo＜R _LLso, or R _o＜R _LLso＜R _LLdo, or R _LLso＜R _o＜R _LLdo, exist simultaneously

Or

Numerical relation; When adopting the salty mud drilling well, invade liquid resistivity less than formation water resistivity, it is " drag reduction " mud invasion character that depth side direction resistivity logging shows the pure water layer, i.e. R _o＞R _LLdo＞R _LLso, exist simultaneously

Or

Numerical relation, the sandstone formation of need judging is the pure water layer.

R wherein _oBe sandstone water layer true resistivity; R _TsBe shallow side direction resistivity logging actual measurement resistivity; R _TdBe dark side direction resistivity logging actual measurement resistivity.

Wherein, behind mud and mud filtrate intrusion net pay zone, with all or part of displacement of the oil in the hole, owing to no original movable water existence, i.e. S in the sandstone pores _Wf=0, obtain the apparent resistivity explanation formula of depth side direction apparent resistivity logging in sandstone formation and be:

$R_{\mathrm{LLs}}=\frac{1}{\frac{{\mathrm{\φ}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+\frac{{\mathrm{\φ}}_t{u}_{\mathrm{fs}}(1-S_{\mathrm{wi}})}{R_{\mathrm{mfs}}}}---\left(7\right)$

$R_{\mathrm{LLd}}=\frac{1}{\frac{{\mathrm{\φ}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+\frac{{\mathrm{\φ}}_t{u}_{\mathrm{fd}}(1-S_{\mathrm{wi}})}{R_{\mathrm{mf}}}}---\left(8\right)$

In the formula: R _LLsThe shallow side direction apparent resistivity of-sandstone oil reservoir, Ω m; R _LLd-sandstone oil reservoir deep lateral apparent resistivity, Ω m; R _Wi-stratum irreducible water resistivity, Ω m; R _Wf-movable formation water resistivity, Ω m; R _MfMud filtrate resistivity under the-formation temperature, Ω m; R _Mfs-comprehensive the resistivity of intrusion liquid of (investigation depth 0.35m) in shallow side direction resistivity logging investigative range, Ω m; S _Wi-sandstone irreducible water saturation, 0＜S _Wi≤ 1; u _Fs-resident fluid the replacement rate of (investigation depth 0.35m) in shallow side direction resistivity logging investigative range, 0≤u _Fs≤ 1; u _Fd-resident fluid the replacement rate of (investigation depth 1.15m) in dark side direction resistivity logging investigative range, 0≤u _Fd≤ 1; v _f-mud filtrate partition factor, 0≤v _f＜1; φ _t-sandstone total porosity, 0＜φ _t≤ 47.64%.

R wherein _tBe sandstone oil reservoir true resistivity, R _TsBe the actual measurement resistivity of shallow side direction resistivity logging, R _TdActual measurement resistivity for dark side direction resistivity logging.

According to R _LLs, R _LLd, R _LLso, R _LLdo, R _t, R _TsAnd R _TdBetween numerical relation, whether can determine sandstone formation is net pay zone.

Further, when adopting the fresh water mud drilling well, invade liquid resistivity greater than formation water resistivity, it is " drag reduction " (successively decreasing) mud invasion character that depth side direction resistivity logging shows the sandstone oil reservoir, promptly has R _t＞R _LLd＞R _LLsNumerical relation, exist simultaneously Or

Numerical relation, the sandstone formation of need judging is a net pay zone.

Additionally, for containing movable water oil reservoir, the apparent resistivity explanation formula of depth side direction resistivity logging in sandstone formation is:

$R_{\mathrm{LLs}}^{\′}=\frac{1}{\frac{{\mathrm{\φ}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+\frac{{\mathrm{\φ}}_t{u}_{\mathrm{fs}}(1-S_{\mathrm{wi}})}{R_{\mathrm{mfs}}}+\frac{{\mathrm{\φ}}_t(1-u_{\mathrm{fs}})S_{\mathrm{wf}}}{R_{\mathrm{wf}}}}---\left(9\right)$

$R_{\mathrm{mfs}}=\frac{1}{\frac{v_f}{R_{\mathrm{mf}}}+\frac{1-v_f}{R_m}}---\left(10\right)$

$R_{\mathrm{LLd}}^{\′}=\frac{1}{\frac{{\mathrm{\φ}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+\frac{{\mathrm{\φ}}_t{u}_{\mathrm{fd}}(1-S_{\mathrm{wi}})}{R_{\mathrm{mf}}}+\frac{{\mathrm{\φ}}_t(1-u_{\mathrm{fd}})S_{\mathrm{wf}}}{R_{\mathrm{wf}}}}---\left(11\right)$

In the formula: R ' _LLs-contain the shallow side direction apparent resistivity of movable water oil reservoir, Ω m; R ' _LLd-contain movable water oil reservoir deep lateral apparent resistivity, Ω m; S _Wf-sandstone formation mobile water saturation, 0≤S _Wf＜1; R _Wi-stratum irreducible water resistivity, Ω m; R _Wf-movable formation water resistivity, Ω m; R _MfMud filtrate resistivity under the-formation temperature, Ω m; R _Mfs-comprehensive the resistivity of intrusion liquid of (investigation depth 0.35m) in shallow side direction resistivity logging investigative range, Ω m; S _Wi-sandstone irreducible water saturation, 0＜S _Wi≤ 1; u _Fs-resident fluid the replacement rate of (investigation depth 0.35m) in shallow side direction resistivity logging investigative range, 0≤u _Fs≤ 1; u _Fd-resident fluid the replacement rate of (investigation depth 1.15m) in dark side direction resistivity logging investigative range, 0≤u _Fd≤ 1; v _f-mud filtrate partition factor, 0≤v _f＜1; φ _t-sandstone total porosity, 0＜φ _t≤ 47.64%.

According to R _LLs, R _LLd, R _LLso, R _LLdo, R ' _LLs, R ' _LLd, R _t, R _TsAnd R _TdBetween numerical relation, whether can determine sandstone formation is to contain movable water oil reservoir.

Further, there is R when depth side direction resistivity _LLd＞R ' _LLd＞R _LLdoOr R _LLs＞R ' _LLs＞R _LLsoNumerical relation, exist simultaneously

Or

Numerical relation, the sandstone formation that judge is to contain movable water oil reservoir.

Preferably, obtain the sandstone total porosity according to core analysis factor of porosity and actual measurement interval transit time.

Preferably, adopt the pressure mercury data that not influenced by borehole environment to determine irreducible water saturation.

Preferably, local water density is the function of reservoir water salinity, formation temperature and reservoir pressure.

Preferably, movable formation water resistivity R _WfAfter can being converted to equivalent N aCl salinity by reservoir water salinity (SC), concern that by the total mineralization and the equivalent coefficient of corresponding ion concentration plate looks into and get (synoptic diagram that accompanying drawing 7 concerns for local water equivalent N aCl solution resistance rate and temperature intersection for the equivalent coefficient graph of a relation of total mineralization and corresponding ion concentration, accompanying drawing 8).

Preferably, the relation between pure water layer irreducible water saturation and the mobile water saturation:

1-S _wi＝S _wc

In the formula: S _WcMobile water saturation in-100% water bearing sand, 0≤S _Wc＜1; S _Wi-sandstone irreducible water saturation, 0＜S _Wi≤ 1.

In order further to describe technical scheme of the present invention in detail, describe below how to obtain adopting among the present invention various lithologic parameters.But explanation and explanation to various lithologic parameters below it will be apparent to those skilled in the art that only are exemplary; Those skilled in the art also can obtain sandstone total porosity, irreducible water saturation, local water density and in-place oil density by other approach, and formation water resistivity, mud and mud filtrate resistivity and sandstone formation true resistivity etc.

1 lithologic parameter

1.1 sandstone total porosity

People such as J.P.Martin in 1986 have proposed acoustic Formation Factor Formula on the basis of Raymer-Hunt work:

$\frac{\mathrm{\Δt}}{{\mathrm{\Δt}}_{\mathrm{ma}}}=\frac{1}{{(1-{\mathrm{\φ}}_t)}^x}---\left(12\right)$

In the formula: φ _t-sandstone total porosity, f; Δ t-sandstone interval transit time, μ s/m; Δ t _Ma-sandstone skeleton the time difference, μ s/m; The x-lithology factor.

Because (12) formula and formation resistivity factor Very similar, so Be called " sound wave formation factor ".

(12) formula can be exchanged into:

${\mathrm{\φ}}_t=1-{\left(\frac{{\mathrm{\Δt}}_{\mathrm{ma}}}{\mathrm{\Δt}}\right)}^{1/x}---\left(13\right)$

According to core analysis total porosity (φ _c) and actual measurement interval transit time (Δ t), try to achieve the skeleton time difference (Δ t by the formula form recurrence after the following formula distortion _Ma) and lithology factor (x), finally obtain factor of porosity and explain equation.

logΔt＝logΔt _ma-xlog(1-φ _c) (14)

In the formula: φ _cThe sandstone total porosity of-lab analysis, f.

1.2 irreducible water saturation

Utilize J (S _w) function classifies to capillary pressure curve and average.

$J\left(S_w\right)=\frac{P_c}{\mathrm{\σ}\cos \mathrm{\θ}}\sqrt{\frac{K}{\mathrm{\φ}}}---\left(15\right)$

In the formula: σ, θ and P _c-be respectively interfacial tension, contact angle and capillary pressure in the laboratory; K-sandstone permeability, mD; φ-sandstone net porosity (lab analysis), f.

If the capillary pressure test is mercury-air system, formula can be converted into:

$J\left(S_w\right)=0.086P_c\sqrt{\frac{K}{\mathrm{\φ}}}---\left(16\right)$

J (S _w) function determines the concrete steps of oil saturation:

The first step is calculated the C value of every sample.

For the capillary pressure curve of a certain sample, σ, cos θ, K and φ are constant.Be convenience of calculation, making the constant term in (16) formula is C, that is:

$C=0.086\sqrt{\frac{K}{\mathrm{\φ}}}---\left(17\right)$

Second step, the J (S of each pressure tap of calculating rock sample _w) functional value.

J(S _w)＝CP _c (18)

The 3rd step, the average mercury saturation degree in statistics different saturation interval, average J (S _w) function and average pressure data.

In the 4th step, represent J (S with ordinate _w) function, horizontal ordinate is represented mercury saturation degree (S _Hg), J (S on the point _w) idea.

If idea is concentrated, illustrate that these samples belong to a kind of pore texture type, can comprehensively be a J (S who represents such reservoir _w) function curve.

The 5th goes on foot, and asks for the average capillary pressure curve of such reservoir.

The mean permeability of known such reservoir and average pore, (17) formula of utilization is calculated average C, then with J (S _w) J (S of arbitrary mercurous saturation degree correspondence on the function curve _w) inverse with average C value on duty, can obtain the corresponding average capillary pressure value (P of this point _c), and then draw the average capillary pressure curve of such reservoir.

${\stackrel{\‾}{P}}_c=\frac{1}{\stackrel{\‾}{C}}J\left(S_w\right)---\left(19\right)$

In the 6th step, the average capillary pressure curve under the laboratory condition is scaled capillary pressure curve under the reservoir condition.

The capillary pressure curve of in the laboratory, measuring, though the sample that uses is the actual rock core of oil reservoir, but non-wetting phase and wetting phase fluid can not directly adopt the fluid under the reservoir condition, and the fluid interface tension force of various combination is different with moisten contact angle, and the capillary pressure that records also has nothing in common with each other.Therefore the capillary pressure curve that no matter adopts oil-water or that system of mercury-air to record when calculating original water saturation, must be proofreaied and correct and is the average capillary pressure curve under the reservoir condition.

Laboratory capillary pressure expression formula and oil reservoir capillary pressure expression formula are respectively:

$P_c=\frac{2\mathrm{\σ}\cos \mathrm{\θ}}{r}---\left(20\right)$

$P_{\mathrm{cR}}=\frac{2{\mathrm{\σ}}_R\cos {\mathrm{\θ}}_R}{r}---\left(21\right)$

In the formula: σ, θ and P _c-be respectively interfacial tension, contact angle and capillary pressure in the laboratory;

σ _R, θ _RAnd P _CR-be respectively interfacial tension, contact angle and capillary pressure under the reservoir condition.

Get by (20), (21) formula simultaneous:

$P_{\mathrm{cR}}=\frac{{\mathrm{\σ}}_R\cos {\mathrm{\θ}}_R}{\mathrm{\σ}\cos \mathrm{\θ}}{P}_c---\left(22\right)$

The interfacial tension of common laboratory fluids and moisten contact angle are known, and owing to reservoir fluid is under stratum higher temperature and the reservoir pressure, and in oil and the water solution gas is arranged, so the interfacial tension of reservoir fluid is difficult to ask for.Moisten contact angle is very difficult under the measurement reservoir condition, in actual use, formula (22) is reduced to:

$P_{\mathrm{cR}}=\frac{{\mathrm{\σ}}_R}{\mathrm{\σ}\cos \mathrm{\θ}}{P}_c---\left(23\right)$

Parameter (table 1 is interfacial tension, contact angle experiments chamber analysis results) substitution (22) formula of lab analysis is got:

$P_{\mathrm{cR}}=\frac{P_c}{17.143}---\left(24\right)$

Table 1 interfacial tension, contact angle experiments chamber analysis results

In the 7th step, the capillary pressure under the reservoir condition is scaled oil column height.

The capillary pressure of oil reservoir is come balance by the gravitational difference of profit, and capillary pressure can be expressed as:

P _cR＝H(ρ _wf-ρ _of)g (25)

In the formula: the above height of the free-water level of H-oil reservoir, m; G-acceleration of gravity, 9.80m/s ²

(25) formula is converted into SI system practical unit, finds the solution the oil-containing height and be:

$H=\frac{{102P}_{\mathrm{cR}}}{{\mathrm{\ρ}}_{\mathrm{wf}}-{\mathrm{\ρ}}_{\mathrm{of}}}---\left(26\right)$

In the formula: P _CR-oil reservoir capillary pressure, MPa; ρ _Wf, ρ _Of-be respectively reservoir condition sub-surface water-mass density and oil density, g/cm ³

In the 8th step, ask for oil saturation.

(24) formula substitution (26) formula is got:

$H=\frac{{5.95P}_c}{{\mathrm{\ρ}}_{\mathrm{wf}}-{\mathrm{\ρ}}_{\mathrm{of}}}---\left(27\right)$

Can obtain oil saturation and oil reservoir height relationships by (27) formula and the average capillary pressure curve of oil reservoir, accompanying drawing 9 is the synoptic diagram of oil saturation and oil reservoir height relationships, and then tries to achieve oil saturation and irreducible water saturation under the different oil reservoir height.

1.3 local water density

Local water density is the function of reservoir water salinity, formation temperature and reservoir pressure, determines that the stratum water-mass density adopts following method.

${\mathrm{\&rho;}}_{\mathrm{wf}}=\frac{{\mathrm{\&rho;}}_{\mathrm{<figure-callout\; id="1"\; label="w"\; filenames="HSA00000012998500041.png,HSA00000012998500051.png"\; state="\{\{state\}\}">w</figure-callout>}1}}{62.4{\mathrm{\&rho;}}_{\mathrm{<figure-callout\; id="2"\; label="w"\; filenames="HSA00000012998500041.png,HSA00000012998500051.png"\; state="\{\{state\}\}">w</figure-callout>}2}{\mathrm{\&rho;}}_{w3}}---\left(28\right)$

${\mathrm{\&rho;}}_{\mathrm{<figure-callout\; id="1"\; label="w"\; filenames="HSA00000012998500041.png,HSA00000012998500051.png"\; state="\{\{state\}\}">w</figure-callout>}1}={10}^{3.05\&times;{10}^{-7}S\&CenterDot;C+1.745}---\left(29\right)$

ρ _w2＝1-1.87×10 ^-5t-1.063×10 ^-6t ² (30)

ρ _w3＝1-1.40×10 ^-5t-2.40×10 ^-6P _i (31)

Reservoir water salinity in the formula (SC) is a measured value, and formation temperature (t) is the statistic equation calculating according to the actual measurement formation temperature of research area and vertical depth (D/m) foundation:

$t=10+\frac{3.2D}{100}---\left(32\right)$

Reservoir pressure (P in the formula _i/ MPa) be to calculate according to the statistic equation that Measured formation pressure and vertical depth (D/m) are set up:

P _i＝0.10+0.01D (33)

1.4 in-place oil density

Under formation condition,, make the variation of in-place oil density and solution gas content closely related owing to be dissolved with rock gas to some extent in the crude oil.Analyze data according to the actual sample that the laboratory provides, original solution gas-oil ratio and in-place oil density ASSOCIATE STATISTICS equation (299 group analysis data) have been set up, accompanying drawing 10 is the synoptic diagram of original solution gas-oil ratio and in-place oil density relationship, and related coefficient is 0.9090.

ρ _of＝1.30736-0.29186logR _si (34)

R _si＝321(B _oi-1) ^1.02 (35)

$B_{\mathrm{oi}}=1+0.01176{\mathrm{\&rho;}}_o^{-16.22}+0.01589D^{3.27}---\left(36\right)$

In the formula: ρ _o-ground oil density (actual analysis value), g/cm ³B _Oi-oil volume factor, zero dimension; R _Si-original solution gas-oil ratio (analyze or ask for) by formula by PVT, m ³/ m ³The D-vertical depth, m.

1.5 formation water resistivity

Comprehensive formation water resistivity (R _w ^*) and movable formation water resistivity (R _Wf) can obtain by lab analysis:

$R_w^{*}=R_o{\mathrm{\&phi;}}_t---\left(37\right)$

In the formula: R _o-100% water bearing sand true resistivity, Ω m; φ _t-sandstone total porosity, f.

Movable formation water resistivity (R _Wf):

$R_{\mathrm{wf}}=\frac{R_w^{*}+(1-S_{\mathrm{wi}})R_w^{*}}{2}---\left(38\right)$

Stratum irreducible water resistivity (R _Wi) acquiring method:

$R_{\mathrm{wi}}=\frac{S_{\mathrm{wi}}{R}_w^{*}{R}_{\mathrm{wf}}}{R_{\mathrm{wf}}-(1-S_{\mathrm{wi}})R_w^{*}}---\left(39\right)$

In addition, movable formation water resistivity (R _Wf) also can be converted to equivalent N aCl salinity by reservoir water salinity (SC) after, equivalent coefficient graph of a relation by total mineralization and corresponding ion concentration obtains, and accompanying drawing 7 be a total mineralization with equivalent coefficient graph of a relation, the accompanying drawing 8 of corresponding ion concentration is the synoptic diagram that local water equivalent N aCl solution resistance rate and temperature intersection concern.

1.6 mud and mud filtrate resistivity

Mud and mud filtrate resistivity be according to the actual measurement mud resistivity (R ' _m), mud density (ρ _m), mud temperature parameters such as (t '), the formula that adopts formula that D.W.Hilchie proposed in 1984 and Schlumberger company to provide calculates to be asked for:

$R_m=\left(\frac{t^{\&prime;}+z}{t+z}\right)R_m^{\&prime;}---\left(40\right)$

$z=\frac{160}{9}+{10}^{0.3861545-0.340396\lg R_m^{\&prime;}}---\left(41\right)$

$R_{\mathrm{mf}}={\mathrm{<figure-callout\; id="10"\; label="R\; m"\; filenames="HSA00000012998500011.png,HSA00000012998500012.png"\; state="\{\{state\}\}">R</figure-callout>}}_m{10}^{0.396-0.0475{\mathrm{\&rho;}}_m}---\left(42\right)$

In the formula: R _mMud resistivity under the-formation temperature, Ω m; T '-well head observed temperature, ℃; The t-formation temperature, ℃; R ' _m-well head actual measurement mud resistivity, Ω m; ρ _m-well head mud density, g/cm ³R _MfMud filtrate resistivity under the-formation temperature, Ω m; The z-conversion coefficient.

1.7 sandstone formation true resistivity

(1) pure water layer true resistivity

Relation between pure water layer irreducible water saturation and the mobile water saturation:

1-S _wi＝S _wc (43)

In the formula: S _WcMobile water saturation in-100% water bearing sand stratum, f; S _Wi-sandstone irreducible water saturation, f.Pure water layer resistivity accounting equation:

$R_o=\frac{\mathrm{<figure-callout\; id="2"\; label="N\; N"\; filenames="HSA00000012998500041.png,HSA00000012998500051.png"\; state="\{\{state\}\}">N</figure-callout>}}{\frac{N^{}}{}}---\left(44\right)$

In the formula: the number of slices of N-unit cube sandstone, individual.

(44) get behind the formula abbreviation:

$R_o=\frac{R_{\mathrm{wi}}{R}_{\mathrm{wf}}}{R_{\mathrm{wi}}{\mathrm{\&phi;}}_t(1-S_{\mathrm{wi}})+R_{\mathrm{wf}}{\mathrm{\&phi;}}_t{S}_{\mathrm{wi}}}---\left(45\right)$

In the formula: R _o-100% water bearing sand true formation resistivity, Ω m.

(2) contain movable water oil reservoir true resistivity

Contain the relation between irreducible water saturation, mobile water saturation and the oil saturation of movable water oil reservoir:

1-S _wi-S _of＝S _wf (46)

According to (46) formula, must contain movable water oil reservoir true resistivity equation:

$R_t^{\&prime;}=\frac{R_{\mathrm{wi}}{R}_{\mathrm{wf}}}{R_{\mathrm{wi}}{\mathrm{\&phi;}}_t(1-S_{\mathrm{wi}}-S_{\mathrm{of}})+R_{\mathrm{wf}}{\mathrm{\&phi;}}_t{S}_{\mathrm{wi}}}---\left(47\right)$

In the formula: R ' _t-moisture oil reservoir true resistivity, Ω m.

(3) net pay zone true resistivity

1-S _wi＝S _oi (48)

According to (48) formula, get net pay zone true resistivity equation:

$R_t=\frac{R_{\mathrm{wi}}}{{\mathrm{\&phi;}}_t{S}_{\mathrm{wi}}}---\left(49\right)$

In the formula: R _t-sandstone oil reservoir true resistivity, Ω m.

1.8 resident fluid replacement rate and mud filtrate partition factor

(1) resident fluid replacement rate

$u_f=\frac{V_F}{V}$

V _F＝V _mud+V _mf

In the formula: V _FThe mud and the mud filtrate cumulative volume on stratum (referring to side direction resistivity logging investigation depth), m are specified in-intrusion ³V _MudThe mud volume in the stratum (referring to side direction resistivity logging investigation depth), m are specified in-intrusion ³V _MfThe mud filtrate volume in the stratum (referring to side direction resistivity logging investigation depth), m are specified in-intrusion ³V-specifies movable fluid volume in the virgin zone (referring to side direction resistivity logging investigation depth), m ³u _f-resident fluid replacement rate, 0≤u _f≤ 1, zero dimension.

Intrusion amount (the V of mud and mud filtrate in the stratum _F) be leak-off mud volume (V in the drilling process _L) with well drilling detritus in movable fluid volume (V _Sc) poor:

V _F＝V _L-V _sc (50)

Movable fluid volume (V in the well drilling detritus _Sc):

$V_{\mathrm{sc}}=\mathrm{\&pi;}{\left(\frac{d}{2}\right)}^{2}h{\stackrel{\&OverBar;}{\mathrm{\&phi;}}}_t(1-{\stackrel{\&OverBar;}{S}}_{\mathrm{wi}})---\left(51\right)$

In the formula: V _ScMovable fluid volume in the-well drilling detritus, m ³The d-bit diameter generally adopts the drill bit of diameter 9.5in, i.e. 9.5 * 0.0254=0.2413m; H-permeable strata thickness, m; φ _t-permeable strata thickness balance factor of porosity, f; S _Wi-permeable strata hole thickness balance irreducible water saturation, f.

Leak-off mud volume (V in the drilling process _L):

Because the volume (V of mud and landwaste, formation pore fluid exchange _L) equate, then mix and (composition) mud density of local water and can be expressed as:

${\mathrm{\&rho;}}_{\mathrm{mc}}=\frac{({\mathrm{\&alpha;V}}_m-V_L){\mathrm{\&rho;}}_m+V_L{\mathrm{\&rho;}}_f}{{\mathrm{\&alpha;V}}_m}---\left(52\right)$

In the formula: ρ _Mc-mixed the mud density (actual analysis) of resident fluid, g/cm ³ρ _m-wellbore mud density (on-the site analysis), g/cm ³ρ _f-under surface condition, in the virgin zone by the pore fluid density (actual analysis) of displacement, g/cm ³V _mThe mud volume of-Theoretical Calculation, m ³V _LThe mud volume of leak-off in the-drilling process, m ³The ratio of actual use amount of mud and Theoretical Calculation mud volume in α-drilling process (employed mud volume is 2.0～3.0 times of calculated value in the actual well drilled process), zero dimension.

Change by (52) formula:

$V_L=\frac{{\mathrm{\&alpha;V}}_m({\mathrm{\&rho;}}_m-{\mathrm{\&rho;}}_{\mathrm{mc}})}{{\mathrm{\&rho;}}_m-{\mathrm{\&rho;}}_f}---\left(53\right)$

Three kinds of density in the following formula (in wellbore mud density, the virgin zone by the pore fluid density of displacement, mixed the mud density of resident fluid) all can in time be obtained by the on-site sampling analysis, and the needed mud volume of drilling process can pass through Theoretical Calculation.

$V_m=\mathrm{\&pi;L}{\left(\frac{d}{2}\right)}^2---\left(54\right)$

In the formula: π-circular constant, value 3.1415; L-wellbore trace length (straight well is drilling depth D), m.

Mud intrusion amount in the depth side direction resistivity logging investigative range: shallow side direction resistivity logging investigation depth is 0.35m, and dark side direction resistivity logging investigation depth is 1.15m.

Movable fluid volume in the shallow side direction resistivity logging investigative range:

$V_s=\mathrm{\&pi;h}{\stackrel{\&OverBar;}{\mathrm{\&phi;}}}_t(1-{\stackrel{\&OverBar;}{S}}_{\mathrm{wi}})[r_s^2-{\left(\frac{d}{2}\right)}^2]---\left(55\right)$

In the formula: r _s-shallow side direction resistivity logging radius of investigation (equal investigation depth and add the drill bit radius), m; V _sMovable fluid volume in the-shallow side direction resistivity logging investigative range, m ³

Movable fluid volume in the dark side direction resistivity logging investigative range (containing the movable fluid volume in the shallow side direction resistivity logging investigative range):

$V_d=\mathrm{\&pi;h}{\stackrel{\&OverBar;}{\mathrm{\&phi;}}}_t(1-{\stackrel{\&OverBar;}{S}}_{\mathrm{wi}})[r_d^2-{\left(\frac{d}{2}\right)}^2]---\left(56\right)$

In the formula: r _d-dark side direction resistivity logging radius of investigation (equal investigation depth and add the drill bit radius), m; V _dMovable fluid volume in the-dark side direction resistivity logging investigative range, m ³

The number percent ω of leak-off mud total amount is distributed in the shallow side direction resistivity investigative range in pit shaft, and remaining 1-ω is distributed in the dark side direction resistivity investigative range, that is:

$u_{\mathrm{fs}}=\frac{{\mathrm{\&omega;V}}_F}{\mathrm{\&pi;h}{\stackrel{\&OverBar;}{\mathrm{\&phi;}}}_t(1-{\stackrel{\&OverBar;}{S}}_{\mathrm{wi}})[r_s^2-{\left(\frac{d}{2}\right)}^2]}---\left(57\right)$

$u_{\mathrm{fd}}=\frac{(1-\mathrm{\&omega;})V_F}{\mathrm{\&pi;h}{\stackrel{\&OverBar;}{\mathrm{\&phi;}}}_t(1-{\stackrel{\&OverBar;}{S}}_{\mathrm{wi}})(r_d^2-r_s^2)}---\left(58\right)$

In the formula: u _Fs-resident fluid the replacement rate of (investigation depth 0.35m) in shallow side direction resistivity logging investigative range, 0≤u _Fs≤ 1; u _Fd-resident fluid the replacement rate of (investigation depth 1.15m) in dark side direction resistivity logging investigative range, 0≤u _Fd≤ 1; Intrusion liquid proportional (analysis is obtained) in ω-intrusion depth side direction investigative range, 0＜ω≤1.

(2) mud filtrate partition factor expression formula

$v_f=\frac{V_{\mathrm{mf}}}{V_{\mathrm{mud}}+V_{\mathrm{mf}}}$

In the formula: v _f-mud filtrate partition factor, 0≤v _f＜1.

The mud filtrate partition factor is to be obtained according to core analysis by the laboratory.

1.9 pure water layer apparent resistivity equation

After mud and mud filtrate were invaded sandstone formation, with the movable water section in the hole or all displacements, the mobile water saturation in the stratum was S at this moment _Wc=1-S _Wi, then can derive the apparent resistivity of depth side direction resistivity logging in sandstone formation and explain equation by (44) formula:

$R_{\mathrm{LLso}}=\frac{1}{\frac{{\mathrm{\&phi;}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+{\mathrm{\&phi;}}_t(\frac{u_{\mathrm{fs}}}{R_{\mathrm{mfs}}}+\frac{1-u_{\mathrm{fs}}}{R_{\mathrm{wf}}})(1-S_{\mathrm{wi}})}$

$R_{\mathrm{mfs}}=\frac{1}{\frac{v_f}{R_{\mathrm{mf}}}+\frac{1-v_f}{R_m}}$

$R_{\mathrm{LLdo}}=\frac{1}{\frac{{\mathrm{\&phi;}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+{\mathrm{\&phi;}}_t(\frac{u_{\mathrm{fd}}}{R_{\mathrm{mf}}}+\frac{1-u_{\mathrm{fd}}}{R_{\mathrm{wf}}})(1-S_{\mathrm{wi}})}$

Invade mud and mud filtrate hybrid resistor rate (R in the formula _Mfs) be to ask for according to the mud filtrate partition factor calculating of lab analysis.

1.10 contain movable water oil reservoir apparent resistivity equation

After mud and mud filtrate are invaded and contained movable water oil reservoir,, can derive the apparent resistivity of depth side direction resistivity logging in sandstone formation by (44) formula and explain equation movable water in the hole and the part or all of displacement of oil:

$R_{\mathrm{LLs}}^{\&prime;}=\frac{1}{\frac{{\mathrm{\&phi;}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+\frac{{\mathrm{\&phi;}}_t{u}_{\mathrm{fs}}(1-S_{\mathrm{wi}})}{R_{\mathrm{mfs}}}+\frac{{\mathrm{\&phi;}}_t(1-u_{\mathrm{fs}})S_{\mathrm{wf}}}{R_{\mathrm{wf}}}}$

$R_{\mathrm{LLd}}^{\&prime;}=\frac{1}{\frac{{\mathrm{\&phi;}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+\frac{{\mathrm{\&phi;}}_t{u}_{\mathrm{fd}}(1-S_{\mathrm{wi}})}{R_{\mathrm{mf}}}+\frac{{\mathrm{\&phi;}}_t(1-u_{\mathrm{fd}})S_{\mathrm{wf}}}{R_{\mathrm{wf}}}}$

1.11 net pay zone apparent resistivity equation

After mud and mud filtrate are invaded net pay zone, with the part or all of displacement of the oil in the hole, owing to no original movable water existence, i.e. S in the sandstone pores _Wf=0, can derive the apparent resistivity of depth side direction resistivity logging in sandstone formation according to (44) formula and explain equation:

$R_{\mathrm{LLs}}=\frac{1}{\frac{{\mathrm{\&phi;}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+\frac{{\mathrm{\&phi;}}_t{u}_{\mathrm{fs}}(1-S_{\mathrm{wi}})}{R_{\mathrm{mfs}}}}$

$R_{\mathrm{LLd}}=\frac{1}{\frac{{\mathrm{\&phi;}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+\frac{{\mathrm{\&phi;}}_t{u}_{\mathrm{fd}}(1-S_{\mathrm{wi}})}{R_{\mathrm{mf}}}}$

2 oil-water-layer evaluations

By above-mentioned series resistance rate is explained equation of analysis, (invade liquid resistivity) when adopting the fresh water mud drilling well greater than formation water resistivity, after mud and mud filtrate are invaded water layer, the resistivity of sandstone water layer is increased; Mud is invaded oil reservoir can make the resistivity of sandstone oil reservoir reduce.And when adopting the salty mud drilling well (invading liquid resistivity less than formation water resistivity), mud can make the resistivity of sandstone water layer reduce after invading water layer; Invading oil reservoir can make the resistivity of sandstone oil reservoir reduce significantly.

2.1 water layer evaluation method

When adopting the fresh water mud drilling well, invade liquid resistivity greater than formation water resistivity, it is " increasing resistance " (increase progressively or increase) or " partly increasing resistance " mud invasion character that depth side direction resistivity logging shows water layer, i.e. R _o＜R _LLdo＜R _LLso, accompanying drawing 1 is invaded the first synoptic diagram A of the trapezoidal trend of water layer for fresh water mud, or R _o＜R _LLso＜R _LLdo, accompanying drawing 2 is invaded the second synoptic diagram B of the trapezoidal trend of water layer for fresh water mud, or R _LLso＜R _o＜R _LLDo, accompanying drawing 3 is invaded the 3rd synoptic diagram C of the trapezoidal trend of water layer for fresh water mud; When adopting the salty mud drilling well, invade liquid resistivity less than formation water resistivity, it is " drag reduction " mud invasion character that depth side direction resistivity logging shows water layer, i.e. R _o＞R _LLdo＞R _LLso, accompanying drawing 4 is invaded the 4th synoptic diagram D of the trapezoidal trend of water layer for salty mud; When satisfying one of above-mentioned four kinds of situations, also to satisfy

Or

Numerical relation, this is the characteristic feature of identification pure water layer, subordinate list 2 is the fluid properties recognition mode.

2.2 reservoir evaluation method

When adopting the fresh water mud drilling well, invade liquid resistivity greater than formation water resistivity, it is " drag reduction " (successively decreasing) mud invasion character that depth side direction resistivity logging shows the sandstone oil reservoir, i.e. R _t＞R _LLd＞R _LLs, accompanying drawing 5 is invaded the 5th synoptic diagram E of the trapezoidal trend of oil reservoir for fresh water mud; When adopting the salty mud drilling well, invade liquid resistivity less than formation water resistivity, it is " drag reduction " (speed subtracts) mud invasion character that depth side direction resistivity logging shows oil reservoir, i.e. R _t＞R _LLd＞＞R _LLs, accompanying drawing 6 is invaded the 6th synoptic diagram F of the trapezoidal trend of oil reservoir for salty mud; When satisfying one of above-mentioned two kinds of situations, also to satisfy Or Numerical relation, this be identification net pay zone characteristic feature, subordinate list 2 is the fluid properties recognition mode.

2.3 contain movable water reservoir evaluation method

Adopt fresh water or or during the salty mud drilling well, invade liquid invade contain movable water oil reservoir after, R appears _LLd＞R ' _LLd＞R _LLdoOr R _LLs＞R ' _LLs＞R _LLsoNumerical relation, exist simultaneously

Or Numerical relation, this is that identification contains the characteristic feature of movable water oil reservoir, subordinate list 2 is the fluid properties recognition mode.

Table 2 fluid properties recognition mode

Sandstone formation resistivity calculating that the present invention proposes and pore fluid recognition technology are after surplus the Liaohe Oil Field ten block used, effect is remarkable, make oil-water-layer well logging interpretation precision bring up to 90%～95% by original 75%～80%, interpretation coincidence rate has improved 15 percentage points.Following illustrative example is illustrated.

Example 1

Certain herbaceous plants with big flowers east 103 well finishing drilling well depth 1813.4m, 26.9 ℃ of wellhead temperatures, well head mud density 1.160g/cm ³, mixing mud density 1.145g/cm ³, belong to salty mud.48.38 ℃ of the corresponding formation temperatures of No. 42 floor buried depths, mud use amount 207m ³Mud resistivity 0.187 Ω m under the formation temperature, mud filtrate resistivity 0.410 Ω m under the formation temperature, bit diameter 9.5in (0.2413m), total porosity 31.1%, irreducible water saturation 40.0%, comprehensive formation water resistivity 1.771 Ω m, irreducible water resistivity 2.833 Ω m, movable formation water resistivity 1.417 Ω m, 100% pure water layer resistivity, 5.70 Ω m, mud filtrate partition factor 0.80, the shallow replacement rate of layer fluid laterally 15.1%, deeply laterally the layer fluid replacement rate 1.2%, invade liquid 65% and be distributed in shallow side direction investigative range, true resistivity is 22.80 Ω m (not to be subjected to mud contamination) when calculating this layer oil-containing thus, and calculating is subjected to the water layer depth side direction apparent resistivity behind the mud contamination: R _LLsoBe 4.14 Ω m, R _LLdoBe 5.58 Ω m, have 4.14 (R this moment _LLso)＜5.58 (R _LLdo)＜5.70 (R _o) numerical relation, exist simultaneously

Or

Numerical relation, be " drag reduction " mud invasion character that salty mud is invaded water layer, No. 42 floor are used the present invention and are interpreted as water layer.Confirm it is water layer through formation testing: daily output water 86.1 sides, accompanying drawing 11 is certain herbaceous plants with big flowers east 103 well interpretation results synoptic diagram.

Example 2

Brocade 315 well finishing drilling well depth 1250m, 20 ℃ of wellhead temperatures, well head mud density 1.100g/cm ³, mixing mud density 1.095g/cm ³, belong to inclined to one side salty mud.27.46 ℃ of the corresponding formation temperatures of No. 17 floor buried depths, mud use amount 114m ³Mud resistivity 0.827 Ω m under the formation temperature, mud filtrate resistivity 1.824 Ω m under the formation temperature, bit diameter 9.5in (0.2413m), total porosity 31.0%, irreducible water saturation 40.0%, comprehensive formation water resistivity 2.002 Ω m, irreducible water resistivity 3.204 Ω m, movable formation water resistivity 1.602 Ω m, 100% pure water layer resistivity, 6.45 Ω m, mud filtrate partition factor 0.75, the shallow replacement rate of layer fluid laterally 5.0%, deeply laterally the layer fluid replacement rate 1.7%, invade liquid 30% and be distributed in shallow side direction investigative range, true formation resistivity is 25.80 Ω m (not to be subjected to mud contamination) when calculating this layer oil-containing thus, calculates and pollutes stratum, back depth side direction apparent resistivity: R _LLsBe 22.06 Ω m, R _LLdBe 24.68 Ω m, have 22.06 (R this moment _LLs)＜24.68 (R _LLd)＜25.80 (R _t) numerical relation, exist simultaneously

Or

Numerical relation, promptly pollute back oil reservoir " drag reduction " mud invasion character, No. 17 former water layer that is interpreted as of floor adopt the present invention to be interpreted as oil reservoir, confirm that through formation testing No. 17 floor of bright and beautiful 315 wells are oil reservoirs: tired produce oil 61.632 sides, accompanying drawing 12 is bright and beautiful 315 well interpretation results synoptic diagram.

Example 3

Certain herbaceous plants with big flowers east 102 well finishing drilling well depth 2306m, 18 ℃ of wellhead temperatures, well head mud density 1.190g/cm ³, mixing mud density 1.181g/cm ³, belong to salty mud.60.93 ℃ of the corresponding formation temperatures of No. 50 floor vertical depths, mud use amount 263m ³Mud resistivity 0.087 Ω m under the formation temperature, mud filtrate resistivity 0.190 Ω m under the formation temperature, bit diameter 9.5in (0.2413m), total porosity 28.7%, irreducible water saturation 45.0%, comprehensive formation water resistivity 2.239 Ω m, irreducible water resistivity 3.471 Ω m, movable formation water resistivity 1.736 Ω m, 100% pure water layer resistivity, 7.80 Ω m, mud filtrate partition factor 0.90, the shallow replacement rate of layer fluid laterally 6.9%, deeply laterally the layer fluid replacement rate 1.5%, invade liquid 40% and be distributed in shallow side direction investigative range, true resistivity is 26.87 Ω m (not to be subjected to mud contamination) when calculating this layer oil-containing thus, and stratum depth side direction apparent resistivity: R is calculated in contaminated back _LLsoBe 5.38 Ω m, R _LLdoBe 7.16 Ω m, have 5.38 (R this moment _LLso)＜7.16 (R _LLdo)＜7.80 (R _o) numerical relation, exist simultaneously

Or

Numerical relation, belong to salty mud and invade " drag reduction " mud invasion character of water layer, No. 50 former oil reservoir that is interpreted as of floor, adopting the present invention to explain is water layer, confirm that through the individual layer formation testing No. 50 floor are water layers: survey liquid level 731.3m, daily output water 97.80 sides, accompanying drawing 13 is certain herbaceous plants with big flowers east 102 well interpretation results synoptic diagram.

As above-mentioned, the clear lithological data acquisition and treatment method that has described the present invention's proposition in detail is discerned water table, oil reservoir or is contained movable water oil reservoir.Although the present invention is described and explained to the preferred embodiments of the present invention in detail, but this area those skilled in the art be appreciated that, under the situation of the spirit and scope of the present invention that do not deviate from the claims definition, can in form and details, make multiple modification.

Claims (11)

1. lithological data acquisition and treatment method, it may further comprise the steps:

1) utilizes geologic prospecting equipment, gather the formation lithology data;

2) based on the formation lithology data that collect, by sandstone total porosity, irreducible water saturation, local water density and the in-place oil density in the lithology data analytical equipment acquisition stratum;

3), obtain to comprise movable formation water resistivity R by the lithology data analytical equipment based on the formation lithology data that collect _WfWith stratum irreducible water resistivity R _WiFormation water resistivity, mud and mud filtrate resistivity and sandstone formation true resistivity;

4) according in step 2) in sandstone total porosity, irreducible water saturation, local water density and the in-place oil density of acquisition, and, obtain resident fluid replacement rate and mud filtrate partition factor according to the formation water resistivity that in step 3), obtains, mud and mud filtrate resistivity and sandstone formation true resistivity; Wherein the resident fluid replacement rate is meant in the specified scope of wellbore, invade intrusion liquid in the perviousness sandstone formation long-pending with the original state sandstone formation in the ratio of movable fluid volume, invade long-pending mud and the mud filtrate volume sum that is meant intrusion of liquid; The mud filtrate partition factor is meant the mud filtrate and the long-pending ratio of intrusion liquid of invading in the perviousness sandstone formation;

Resident fluid replacement rate expression formula:

$u_f=\frac{V_F}{V}---\left(1\right)$

V _F＝V _mud+V _mf (2)

In the formula: V _FMud and the mud filtrate cumulative volume in the stratum, m are specified in-intrusion ³V _MudThe mud volume in the stratum, m are specified in-intrusion ³V _MfThe mud filtrate volume in the stratum, m are specified in-intrusion ³V-specifies movable fluid volume in the virgin zone, m ³u _f-resident fluid replacement rate, 0≤u _f≤ 1, zero dimension;

Mud filtrate partition factor expression formula:

$v_f=\frac{V_{\mathrm{mf}}}{V_{\mathrm{mud}}+V_{\mathrm{mf}}}---\left(3\right)$

5) base area layer fluid replacement rate and mud filtrate partition factor can determine the attribute of sandstone formation, and promptly whether the stratum is pure water layer or net pay zone or contains movable water oil reservoir.

2. according to method described in the claim 1, after wherein mud and mud filtrate were invaded the pure water layer, with all or part of displacement of movable water in the hole, the mobile water saturation in the stratum was S at this moment _Wc=1-S _Wi, the apparent resistivity explanation formula of depth side direction resistivity logging in the pure water layer is:

$R_{\mathrm{LLso}}=\frac{1}{\frac{{\mathrm{\&phi;}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+{\mathrm{\&phi;}}_t(\frac{u_{\mathrm{fs}}}{R_{\mathrm{mfs}}}+\frac{1-u_{\mathrm{fs}}}{R_{\mathrm{wf}}})(1-S_{\mathrm{wi}})}---\left(4\right)$

$R_{\mathrm{mfs}}=\frac{1}{\frac{v_f}{R_{\mathrm{mf}}}+\frac{1-v_f}{R_m}}---\left(5\right)$

$R_{\mathrm{LLdo}}=\frac{1}{\frac{{\mathrm{\&phi;}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+{\mathrm{\&phi;}}_t(\frac{u_{\mathrm{fd}}}{R_{\mathrm{mf}}}+\frac{1-u_{\mathrm{fd}}}{R_{\mathrm{wf}}})(1-S_{\mathrm{wi}})}---\left(6\right)$

In the formula: R _LLsoThe shallow side direction apparent resistivity of-100% water bearing sand, Ω m; R _LLdo-100% water bearing sand deep lateral apparent resistivity, Ω m; R _Wi-stratum irreducible water resistivity, Ω m; R _Wf-movable formation water resistivity, Ω m; R _MfMud filtrate resistivity under the-formation temperature, Ω m; R _Mfs-the comprehensive resistivity of intrusion liquid in shallow side direction resistivity logging investigative range, Ω m; S _Wi-sandstone irreducible water saturation, 0＜S _Wi≤ 1; u _Fs-resident fluid replacement rate in shallow side direction resistivity logging investigative range, 0≤u _Fs≤ 1; u _Fd-resident fluid replacement rate in dark side direction resistivity logging investigative range, 0≤u _Fd≤ 1; v _f-mud filtrate partition factor, 0≤v _f＜1; φ _t-sandstone total porosity, 0＜φ _t≤ 47.64%;

When adopting the fresh water mud drilling well, invade liquid resistivity greater than formation water resistivity, it is " increasing resistance " or " partly increasing resistance " mud invasion character that depth side direction apparent resistivity logging shows water layer, i.e. R _o＜R _LLdo＜R _LLso, or R _o＜R _LLso＜R _LLdo, or R _LLso＜R _o＜R _LLdo, exist simultaneously

Or

Numerical relation; When adopting the salty mud drilling well, invade liquid resistivity less than formation water resistivity, it is " drag reduction " mud invasion character that depth side direction resistivity logging shows water layer, i.e. R _o＞R _LLdo＞R _LLso, exist simultaneously

Or

Numerical relation, the sandstone formation of need judging is the pure water layer;

R wherein _oBe sandstone water layer true resistivity, R _TsBe the actual measurement resistivity of shallow side direction resistivity logging, R _TdActual measurement resistivity for dark side direction resistivity logging.

3. according to method described in the claim 1, wherein behind mud and mud filtrate intrusion net pay zone, with the part or all of displacement of the oil in the hole, owing to no original movable water existence, i.e. S in the sandstone pores _Wf=0, obtain the apparent resistivity explanation formula of depth side direction resistivity logging in sandstone formation and be:

$R_{\mathrm{LLs}}=\frac{1}{\frac{{\mathrm{\&phi;}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+\frac{{\mathrm{\&phi;}}_t{u}_{\mathrm{fs}}(1-S_{\mathrm{wi}})}{R_{\mathrm{mfs}}}}---\left(7\right)$

$R_{\mathrm{LLd}}=\frac{1}{\frac{{\mathrm{\&phi;}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+\frac{{\mathrm{\&phi;}}_t{u}_{\mathrm{fd}}(1-S_{\mathrm{wi}})}{R_{\mathrm{mf}}}}---\left(8\right)$

In the formula: R _LLsThe shallow side direction apparent resistivity of-sandstone oil reservoir, Ω m; R _LLd-sandstone oil reservoir deep lateral apparent resistivity, Ω m; R _Wi-stratum irreducible water resistivity, Ω m; R _Wf-movable formation water resistivity, Ω m; R _MfMud filtrate resistivity under the-formation temperature, Ω m; R _Mfs-the comprehensive resistivity of intrusion liquid in shallow side direction resistivity logging investigative range, Ω m; S _Wi-sandstone irreducible water saturation, 0＜S _Wi≤ 1; u _Fs-resident fluid replacement rate in shallow side direction resistivity logging investigative range, 0≤u _Fs≤ 1; u _Fd-resident fluid replacement rate in dark side direction resistivity logging investigative range, 0≤u _Fd≤ 1; v _f-mud filtrate partition factor, 0≤v _f＜1; φ _t-sandstone total porosity, 0＜φ _t≤ 47.64%;

According to R _LLs, R _LLd, R _LLso, R _LLdo, R _t, R _TsAnd R _TdBetween numerical relation, whether can determine sandstone formation is net pay zone.

4. according to method described in the claim 3, when wherein adopting the fresh water mud drilling well, invade liquid resistivity greater than formation water resistivity, it is " drag reduction " mud invasion character that depth side direction resistivity logging shows the sandstone oil reservoir, promptly has R _t＞R _LLd＞R _LLsOr R _t＞R _LLd＞＞R _LLsNumerical relation, exist simultaneously

Or

Numerical relation, the sandstone formation of need judging is a net pay zone.

5. according to method described in the claim 1, wherein for containing movable water oil reservoir, the apparent resistivity explanation formula of depth side direction resistivity logging in sandstone formation is:

$R_{\mathrm{LLs}}^{\&prime;}=\frac{1}{\frac{{\mathrm{\&phi;}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+\frac{{\mathrm{\&phi;}}_t{u}_{\mathrm{fs}}(1-S_{\mathrm{wi}})}{R_{\mathrm{mfs}}}+\frac{{\mathrm{\&phi;}}_t(1-u_{\mathrm{fs}})S_{\mathrm{wf}}}{R_{\mathrm{wf}}}}---\left(9\right)$

$R_{\mathrm{mfs}}=\frac{1}{\frac{v_f}{R_{\mathrm{mf}}}+\frac{1-v_f}{R_m}}---\left(10\right)$

$R_{\mathrm{LLd}}^{\&prime;}=\frac{1}{\frac{{\mathrm{\&phi;}}_t{S}_{\mathrm{wi}}}{R_{\mathrm{wi}}}+\frac{{\mathrm{\&phi;}}_t{u}_{\mathrm{fd}}(1-S_{\mathrm{wi}})}{R_{\mathrm{mf}}}+\frac{{\mathrm{\&phi;}}_t(1-u_{\mathrm{fd}})S_{\mathrm{wf}}}{R_{\mathrm{wf}}}}---\left(11\right)$

In the formula: R ' _LLs-contain the shallow side direction apparent resistivity of movable water oil reservoir, Ω m; R ' _LLd-contain movable water oil reservoir deep lateral apparent resistivity, Ω m; S _Wf-sandstone formation mobile water saturation, 0≤S _Wf＜1; R _Wi-stratum irreducible water resistivity, Ω m; R _Wf-movable formation water resistivity, Ω m; R _MfMud filtrate resistivity under the-formation temperature, Ω m; R _Mfs-the comprehensive resistivity of intrusion liquid in shallow side direction resistivity logging investigative range, Ω m; S _Wi-sandstone irreducible water saturation, 0＜S _Wi≤ 1; u _Fs-resident fluid replacement rate in shallow side direction resistivity logging investigative range, 0≤u _Fs≤ 1; u _Fd-resident fluid replacement rate in dark side direction resistivity logging investigative range, 0≤u _Fd≤ 1; v _f-mud filtrate partition factor, 0≤v _f＜1; φ _t-sandstone total porosity, 0＜φ _t≤ 47.64%;

According to R _LLd, R _LLs, R _LLso, R _LLdo, R ' _LLs, R ' _LLd, R _TsAnd R _TdBetween numerical relation, whether can determine sandstone formation is to contain movable water oil reservoir.

6. according to method described in the claim 5, when there is R in sandstone formation _LLd＞R ' _LLd＞R _LLdoOr R _LLs＞R ' _LLs＞R _LLsoNumerical relation, exist simultaneously

Or

Numerical relation, the sandstone formation that judge is to contain movable water oil reservoir.

7. according to the arbitrary described method among the claim 1-6, wherein obtain the sandstone formation total porosity according to core analysis factor of porosity and actual measurement interval transit time.

8. according to the arbitrary described method among the claim 1-6, wherein adopt the pressure mercury data that not influenced by borehole environment to determine irreducible water saturation.

9. according to the arbitrary described method among the claim 1-6, wherein local water density is the function of reservoir water salinity, formation temperature and reservoir pressure.

10. according to the arbitrary described method among the claim 1-6, wherein movable formation water resistivity R _WfAfter can being converted to equivalent N aCl salinity by reservoir water salinity, concerning by the total mineralization and the equivalent coefficient of corresponding ion concentration that plate is looked into and get.

11. according to the described method of claim 10, the wherein relation between pure water layer irreducible water saturation and the mobile water saturation: 1-S _Wi=S _Wc

In the formula: S _WcMobile water saturation in-100% water bearing sand, 0≤S _Wc＜1; S _Wi-sandstone irreducible water saturation, 0＜S _Wi≤ 1.

Images