US20210270124A1

US20210270124A1 - Method for distinguishing authenticity of high-pressure physical property parameters of oil reservoirs

Info

Publication number: US20210270124A1
Application number: US17/256,005
Authority: US
Inventors: Haijin Zheng
Original assignee: YANGZHOU JIANGSU OILFIELD RUIDA PETROLEUM ENGINEERING TECHNOLOGY DEVELOPMENT Co Ltd
Current assignee: YANGZHOU JIANGSU OILFIELD RUIDA PETROLEUM ENGINEERING TECHNOLOGY DEVELOPMENT Co Ltd
Priority date: 2019-01-30
Filing date: 2019-08-08
Publication date: 2021-09-02
Also published as: CN109869144A; WO2020155601A1

Abstract

A method for distinguishing an authenticity of high-pressure physical property parameters of oil reservoirs is provided. The method includes: collecting the high-pressure physical property parameters, screening a standard well, calculating a pump efficiency by means of formulas, then calculating an absolute value Δ=|η_theory−η_reality| of each group of absolute errors, and taking a smallest absolute value Δ_minof the absolute error. wherein Δ_min≤ε, ε is a set accuracy, generally, 0≤ε≤0.02, and the high-pressure physical property parameters corresponding to a group of the theoretical pump efficiency η_theoryand the real pump efficiency η_realityare real high-pressure physical property parameters of an oil reservoir to be distinguished. The method provides a fast, simple and practical means for distinguishing correct high-pressure physical property parameters. The method is applicable to an oil well in a certain zone of an oil reservoir to be determined.

Description

TECHNICAL FIELD

The present invention relates to the technical field of petroleum engineering, particularly the technical field of exploration of oil wells.

BACKGROUND

High-pressure physical property parameters are important parameters that are indispensable for determining oil reservoir types, formulating development plans and performing oil reservoir engineering calculation, and are the basis for studying oil field driving types, determining oil field exploitation modes, calculating oil field reserves and selecting oil well working systems.
At present, there are three main methods for obtaining the high-pressure physical property parameters:
1. Laboratory determination: for unsaturated oil reservoirs of which the saturation pressure is lower than the original formation pressure, during oil tests and production tests, when the flowing bottom hole pressure is higher than the saturation pressure, representative formation oil samples are obtained through a bottom hole sampler, and then, the high-pressure physical property parameters are determined in a laboratory.
2. Calculation by means of existing parameter-related charts: if there is no laboratory determination data or representative formation oil samples cannot be obtained, high-pressure physical property parameter values of formation crude oil can be found through charts.
3. Prediction by means of empirical formulas: when sampling conditions are not available and the parameters cannot be found from charts, some empirical formulas published in some countries are generally used to predict high-pressure physical properties.
The laboratory determination requires undeveloped wells and needs to meet a series of harsh conditions, such as bottom hole pressure higher than the expected original saturation pressure, no water or water content of no more than 5%, stable oil and gas flow, and no intermittency. Due to the complex composition of crude oil, the charts are often not accurate enough. The empirical formulas also have the problem of the scope of application. Therefore, the experimental method, the plate method and the empirical formula method all need to meet certain conditions and are not accurate enough. In particular, there is no easy method to determine the authenticity of the high-pressure physical property parameter values.

SUMMARY

An objective of the present invention is to provide a fast, simple and accurate method for distinguishing the authenticity of high-pressure physical property parameters of oil reservoirs.
A technical solution of the present invention includes the following steps:
1) collecting the high-pressure physical property parameters: collecting high-pressure physical property parameter data of the sampled wells of an oil reservoir to be distinguished, including gas-oil ratio GOR, solubility coefficient α, saturation pressure P_b, formation crude oil density ρ_o, and formation crude oil viscosity μ_o;
2) screening a standard well: determining the standard well from multiple pumping wells continuously exploited in the same block and the same zone as the sampled well, the so-called standard well is the one whose indicator diagram of the oil well reflects that the pump does not leak and its oil tube string does not leak, and in the multiple oil wells, the standard well has the lowest rate of water content, the highest liquid yield and the largest submergence in the multiple oil wells;
3) calculating pump efficiency: substituting each group of high-pressure physical property parameters collected in the step 1) into the following formula, and calculating the theoretical pump efficiency η_theoryand the real pump efficiency η_realityof the standard well using the related data of the standard well in the step 2):
$η_{theory} = \frac{S - λ}{S} β - \frac{Q_{leak}}{Q_{theory}}$ $λ = \frac{π D^{2} ρ gh}{4 E} (\frac{L_{1}}{f_{1}} + \frac{L_{2}}{f_{2}} + \dots + \frac{L_{n}}{f_{n}} + \frac{L_{p}}{f_{t}})$ $β = \frac{1}{1 + \frac{(GOR - α P_{s}) (1 - f_{w})}{(10 P_{s} + 1)}}$
when P_s≥P_b, making P_s=P_b, where if its oil tubing string is anchored, L_p/f_tis not included in λ;
η_reality =Q _reality /Q _theory*100%
Q _theory =πD ² ρgSN/4
Q _leak =πDρgδ ³ h/(12L _p/μ)
where
Q_realityrepresents the real output of the oil well;
Q_theoryrepresents the theoretical displacement of the oil well;
Q_leakrepresents the pump leakage;
GOR represents the gas-oil ratio and the unit is m³/m³; α represents the solubility coefficient and the unit is m³/(m³·Mpa); P_brepresents the saturation pressure and the unit is Mpa;
λ represents the stroke loss of the standard well and the unit is m; β represents the gas influence coefficient of the standard well, dimensionless; D represents the pump diameter of the standard well and the unit is m; h represents the pump lift of the standard well and the unit is m; L₁, L₂and L_nrespectively represent the lengths of the first, second and n-level rods of the standard well and the unit is m; f₁, f₂and f_nrespectively represent the cross-sectional areas of the first, second and n-level rods of the standard well and the unit is m²; L_prepresents the pump setting depth of the standard well and the unit is m; f_trepresents the cross-sectional area of the metal part of the tubing of the standard well and the unit is m²; E represents the elastic modulus of the steel of the pumping rod in the standard well, 2.1*10⁷N/cm²; P_srepresents the submergence pressure of the standard well and the unit is Mpa; f_wrepresents the rate of water content in the standard well; S represents the stroke length of the standard well and the unit is m; N represents the stroke frequency of the standard well and the unit is 1/min; g represents the gravity acceleration in the standard well and the unit is m/s²; δ represents the annular gap size between the pump plunger and the pump cylinder of the standard well and the unit is m; L_plrepresents the length of the pump plunger of the standard well and the unit is m; μ represents the hydrodynamic viscosity in the standard well and the unit is Pa. S;
some of the above parameters may be calculated from more basic parameters:
mixed liquid density: ρ=(1−f _w)*ρ_o +f _w*ρ_w
effective head: h=h _dynamic+1000*(p _oil −p _casing)*g/ρ
submergence pressure: P _s =p _casing+(h _setting −h _dynamic)*ρ_o *g/1000
hydrodynamic viscosity: μ=f _w+(1−f _w)*μ_o
where f_wrepresents the water content of the pumped liquid of the standard well and the unit is %; ρ_orepresents the formation crude oil density and the unit is t/m³; ρ_wrepresents the water density and the unit is t/m³; h_dynamicrepresents the dynamic liquid level depth of the standard well and the unit is m; hsettm_grepresents the pump setting depth of the standard well and the unit is m; p_oilrepresents the wellhead oil pressure of the standard well and the unit is Mpa; p_casingrepresents the casing pressure of the standard well and the unit is Mpa; μ_orepresents the formation crude oil viscosity and the unit is Pa. S;
4) screening results: after calculation in the step 3), for a group of standard wells screened in the step 2) and each group of high-pressure physical property parameters, obtaining a group of (η_theory, η_reality) data, and for each group of (η_theory, η_reality) data, calculating the absolute value Δ=|η_theory−η_reality|, and taking the smallest absolute value Δ_minof the absolute error:
Δ_min=min|η_theory−η_reality|
where if Δ_minmeets the following condition:
Δmin≤ε (ε is the set accuracy, generally, 0≤ε≤0.02),
the high-pressure physical property parameters corresponding to this group of (η_theory, η_reality) are the real high-pressure physical property parameters of a certain zone of this oil reservoir. The present invention provides a fast, simple and practical means for distinguishing correct high-pressure physical property parameters. The present invention is applicable to those oil wells in a certain zone of an oil reservoir whose high-pressure property parameters is to be determined.

DETAILED DESCRIPTION OF THE EMBODIMENTS

I. Distinguishing of high-pressure physical property parameters of a certain zone (such as Ef1) of a certain oil reservoir:
1. Collect high-pressure physical property parameters of all oil wells obtained high-pressure physical properties (called sampled wells) in the zone Ef1 of this oil reservoir.


	Gas-oil
	ratio	Solubility	Saturation	Crude oil	Crude oil
	GOR	coefficient α	pressure	density	viscosity
Sampled	(m³/	(m³/	Pb	ρ_o	μ_o
well	m³)	(m³· Mpa))	(Mpa)	(t/m³)	(Pa · S)

Zhuang 2-9	12.8	4.29	2.65	0.8528	15.66
Wei 2-22	18.8	4.24	3.82	0.8633	9.39
Wei 2-15	26	5.61	4.02	0.8662	9.39
Chen 3-7	74.9	2.7	6.96	0.89	20.67

2. Standard wells are determined according to the following steps:
(1) m oil wells with no leakage of the pump and no leakage of the oil tubing string are selected according to the indicator diagram of each of all oil wells of this reservoir zone.
(2) The m oil wells are sorted according to the water content, and the first n oil wells with the lowest water content are selected in the m selected wells.
(3) The n oil wells are sorted according to the liquid yield, and the first p oil wells with the highest liquid yield are selected in the n selected wells.
(4) The p oil wells are sorted according to the submergence, and the first q oil wells with the highest submergence are selected as standard wells in the p selected wells.
The above steps, (2), (3), and (4), are interchangeable.
The results are as follows:


			Dynamic
Standard	Liquid	Water	liquid	Oil	Casing		Stroke	Pump	Pump
well	yield	content	level	pressure	pressure	Stroke	frequency	diameter	depth

Wei 8	17.3	83.1	692.2	0.7	0.6	3	2.96	44	946.64
Ping 3

3. The data collected in step 1 and step 2 is substituted into the following formula, and the theoretical pump efficiency η_theoryand the real pump efficiency η_realityof the standard wells are calculated respectively corresponding to the high-pressure physical property parameters of the sampled wells:
$η_{theory} = \frac{S - λ}{S} β - \frac{Q_{leak}}{Q_{theory}}$ $λ = \frac{π D^{2} ρ gh}{4 E} (\frac{L_{1}}{f_{1}} + \frac{L_{2}}{f_{2}} + \dots + \frac{L_{n}}{f_{n}} + \frac{L_{p}}{f_{t}})$ $β = \frac{1}{1 + \frac{(GOR - α P_{s}) (1 - f_{w})}{(10 P_{s} + 1)}}$
(when P_s≥P_b, making P_s=P_b, where if an oil pipe is anchored, L_p/f_tis not included in λ).
Pump leakage: Q _leak =πDμgδ ³ h/(12L _p/μ).
Theoretical displacement of oil well: Q _theory =πD ² μgSN/4.
Real pump efficiency of oil well: η_reality =Q _reality /Q _theory*100%.
Mixed liquid density: ρ=(1−f _w)*ρ_o +f _w*ρ_w.
Effective head: h=h _dynamic+1000*(p _oil −p _casing)*g/ρ.
Submergence pressure: P _s =p _casing+(h _setting −h _dynamic)*ρ_o *g/1000.
Hydrodynamic viscosity: μ=f _w+(1−f _w)*μ_o
The following results are obtained:


Standard	Sampled	Real pump	Theoretical pump
well	well	efficiency η_reality	efficiency η_theory	Δ

Wei 8 Ping 3	Zhuang 2-9	0.912	0.942	0.030
Wei 8 Ping 3	Wei 2-22	0.912	0.911	0.001
Wei 8 Ping 3	Wei 2-15	0.912	0.894	0.018
Wei 8 Ping 3	Chen 3-7	0.912	0.684	0.228

4. Each group of theoretical pump efficiency η_theoryand real pump efficiency η_realityare substituted into the following formula:
Δ=|η_theory−η_reality|.
The absolute value of the absolute error of each group of theoretical pump efficiency η_theoryand real pump efficiency η_realityis respectively obtained and listed in the 5th column of the above table.
It can be seen from the table that the Δ obtained by the standard well Wei 8 Ping 3 using the high-pressure physical property parameters of the well Wei 2-22 meets the following relationship:
Δ_min=min|η_theory−η_reality|, and Δ_min≤ε (here, assuming that ε=0.001).
It can be confirmed that corresponding to the Wei 2-22, the gas-oil ratio GOR=18.8, the solubility coefficient α=4.24, the saturation pressure P_b=3.82, the formation crude oil density ρ_o=0.8633, and the formation crude oil viscosity μ_o=9.39, which are real high-pressure physical property parameters of the zone Ef1 of this oil reservoir.

Claims

What is claimed is:

1. A method for distinguishing an authenticity of high-pressure physical property parameters of an oil reservoir, comprising the following steps:

1) collecting the high-pressure physical property parameters: collecting the high-pressure physical property parameters of a sampled well of the oil reservoir to be distinguished, wherein the high-pressure physical property parameters comprise a gas-oil ratio GOR, a solubility coefficient α, a saturation pressure P_b, a formation crude oil density ρ_o, and a formation crude oil viscosity μ_o;

2) screening a standard well: determining the standard well from multiple oil wells continuously exploited in a same block and a same zone as the sampled well, wherein an indicator diagram of the standard well reflects a pump of the standard well does not leak and a tubing string of the standard well does not leak, and meanwhile, the standard well has a lowest water content, a highest liquid yield and a largest submergence in the multiple oil wells;

3) calculating a pump efficiency: substituting each group of the high-pressure physical property parameters collected in step 1) and related data of the standard well screened in step 2) into the following formula, and calculating a theoretical pump efficiency η_theoryand a real pump efficiency η_realityof the standard well:

η_{theory} = \frac{S - λ}{S} β - \frac{Q_{leak}}{Q_{theory}}

λ = \frac{π D^{2} ρ gh}{4 E} (\frac{L_{1}}{f_{1}} + \frac{L_{2}}{f_{2}} + \dots + \frac{L_{n}}{f_{n}} + \frac{L_{p}}{f_{t}})

β = \frac{1}{1 + \frac{(GOR - α P_{s}) (1 - f_{w})}{(10 P_{s} + 1)}}

when P_s≥P_b, making P_s=P_b, wherein if an oil pipe is anchored, L_p/f_tis not comprised in λ;

η_reality =Q _reality /Q _theory*100%;

Q _theory =πD ² ρgSN/4;

Q _leak =πDρgδ ³ h/(12L _p/μ);

ρ=(1−f _w)*ρ_o +f _w*ρ_w;

h=h _dynamic+1000*(p _oil −p _casing)*g/ρ;

P _s =p _casing+(h _setting −h _dynamic)*ρ_o *g/1000;

μ=f _w+(1−f _w)*μ_o;

wherein

Q_realityrepresents a real output of an oil well of the multiple oil wells;

Q_theoryrepresents a theoretical displacement of the oil well of the multiple oil wells;

Q_leakrepresents a pump leakage;

the GOR represents the gas-oil ratio and a unit of the GOR is m³/m³; α represents the solubility coefficient and a unit of the solubility coefficient is m³/(m³·Mpa); P_brepresents the saturation pressure and a unit of the saturation pressure is Mpa; ρ_orepresents the formation crude oil density and a unit of the formation crude oil density is t/m³; μ_orepresents the formation crude oil viscosity and a unit of the formation crude oil viscosity is Pa·S;

λ represents a stroke loss of the standard well and a unit of the stroke loss is m; β represents a gas influence coefficient of the standard well, wherein the gas influence coefficient is dimensionless; ρ represents a mixed liquid density of the standard well and a unit of the mixed liquid density is t/m³; h represents a pump lift of the standard well and a unit of the pump lift is m; P_srepresents a submergence pressure of the standard well and a unit of the submergence pressure is Mpa; μ represents a hydrodynamic viscosity of the standard well and a unit of the hydrodynamic viscosity is Pa·S; f_wrepresents a rate of a water content of a pumped liquid in the standard well; h_dynamicrepresents a dynamic liquid level depth of the standard well and a unit of the dynamic liquid level depth is m; D represents a pump diameter of the standard well and a unit of the pump diameter is m; h_settingrepresents a pump setting depth of the standard well and a unit of the pump setting depth is m; p_oilrepresents a wellhead oil pressure of the standard well and a unit of the wellhead oil pressure is Mpa; P_casingrepresents a casing pressure of the standard well and a unit of the casing pressure is Mpa; S represents a stroke length of the standard well and a unit of the stroke length is m; N represents a stroke frequency of the standard well and a unit of the stroke frequency is 1/min; L₁, L₂and L_nrespectively represent a length of a first-level rod, a length of a second-level rod and a length of an n-level rod of the standard well and a unit of the length is m; f₁, f₂and f_nrespectively represent a cross-sectional area of a first-level rod, a cross-sectional area of a second-level rod and a cross-sectional area of an n-level rod of the standard well and a unit of the cross-sectional area is m²; L_prepresents a pump depth of the standard well and a unit of the pump depth is m; f_trepresents a cross-sectional area of a metal part of the oil pipe of the standard well and a unit of the cross-sectional area is m²; E represents an elastic modulus of a steel of a pumping rod in the standard well, and E is 2.1*10⁷N/cm²; g represents a gravity acceleration of the standard well and a unit of the gravity acceleration is m/s²; δ represents an annular gap between a pump plunger and a pump cylinder of the standard well and a unit of the annular gap is m; L_plrepresents a length of the pump plunger of the standard well and a unit of the length of the pump plunger is m; and ρ_wrepresents a water density of the standard well and a unit of the water density is t/m³; and

4) screening results: after the pump efficiency is calculated in the step 3), for the standard well screened in step 2) and the each group of the high-pressure physical property parameters collected in step 1), obtaining a group of the theoretical pump efficiency η_theoryand the real pump efficiency η_reality, and for each group of the theoretical pump efficiency η_theoryand the real pump efficiency η_reality, calculating an absolute value Δ=|η_theory−η_reality| of an absolute error, and taking a smallest absolute value Δ_minof the absolute error:

Δ_min=|η_theory−η_reality|;

wherein if Δ_minmeets a condition: Δ_min≤ε, wherein ε is a set accuracy, and 0≤ε≤0.02,

the high-pressure physical property parameters corresponding to the group of the theoretical pump efficiency η_theoryand the real pump efficiency η_reality, comprising the gas-oil ratio GOR, the solubility coefficient α, the saturation pressure P_b, the formation crude oil density ρ_o, and the formation crude oil viscosity μ_o, are real high-pressure physical property parameters of the oil reservoir to be distinguished.