CN113283092B - Method and system for constructing side-roof water invasion three-dimensional physical model - Google Patents
Method and system for constructing side-roof water invasion three-dimensional physical model Download PDFInfo
- Publication number
- CN113283092B CN113283092B CN202110598395.XA CN202110598395A CN113283092B CN 113283092 B CN113283092 B CN 113283092B CN 202110598395 A CN202110598395 A CN 202110598395A CN 113283092 B CN113283092 B CN 113283092B
- Authority
- CN
- China
- Prior art keywords
- water
- steam
- gas
- steam cavity
- volume
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a method and a system for constructing a three-dimensional physical model of side-roof water invasion, wherein the method comprises the following steps: injecting a first water body into the reaction vessel; the reaction vessel is a vessel for simulating a steam cavity formed by an oil reservoir in the SAGD production process; heating the reaction vessel containing the first body of water to a first temperature to simulate steam chamber conditions; the vapor chamber conditions include: temperature, pressure and gas-liquid ratio of the vapor chamber; the gas-liquid ratio is the ratio of the gas volume in the steam cavity to the water volume; and simulating a top water invasion process by injecting a second water body into the reaction container meeting the steam cavity condition to construct a top water invasion three-dimensional physical model. The method realizes the same steam cavity condition through the reaction vessel, and simulates the top water invasion process by injecting the second water body.
Description
Technical Field
The invention relates to the technical field of petroleum industry, in particular to a method and a system for constructing a side-top water invasion three-dimensional physical model.
Background
The SAGD oil recovery technology (Steam Assisted Gravity Drainage, abbreviated as SAGD) is an oil recovery method in which Steam is injected into an oil reservoir from a vertical well or a horizontal well above a horizontal production well near the bottom of the oil reservoir, and heated crude oil and Steam condensate are produced from the horizontal well at the bottom of the oil reservoir. SAGD is a common method for developing a heavy oil reservoir, steam is continuously injected through an injection well, the steam rises due to buoyancy, the steam is condensed at a steam-oil interface due to heat exchange, heavy oil is heated to reduce viscosity, the heavy oil is extracted from an extraction well, and a reservoir model with a steam cavity, hot oil and cold oil sequentially formed from top to bottom is formed due to the super-covering effect of the steam along with the production.
A top-side water heavy oil reservoir generally refers to a reservoir in which the crude oil viscosity is greater than 50mPa · s or the de-gassed crude oil viscosity is greater than 100mPa · s under oil layer conditions, and there is an external water body at the top and the sides. In the application process of the SAGD oil extraction technology in the top water heavy oil reservoir, the top water plays a very adverse role in the development process of the top water heavy oil reservoir. The limit top water is located outside the heavy oil reservoir, and the cold thickened oil that solidifies is being separated between SAGD's steam chamber, along with the extension of steam chamber, the temperature is constantly transmitted, and the temperature of the cold thickened oil of wall constantly rises, and thickness constantly reduces, under the influence of action of gravity and inside and outside pressure differential, takes place limit top water and invades. SAGD development is carried out on the side-top water heavy oil reservoir, if side-top water invasion occurs, temperature and pressure changes of a steam cavity are directly influenced, the size and the existence of the steam cavity are determined, and the later development of the oil reservoir is indirectly influenced. Because of the large difference in temperature between steam and the side roof bodies of water, the steam cavity can change dramatically and this adverse change can become more severe over time, with the temperature and pressure of the steam cavity affecting SAGD production. Therefore, the physical simulation of the top water heavy oil reservoir at the SAGD opening side is very important, and the design of a rapid and reasonable laboratory physical simulation system is the premise of preventing the aggravation of top water damage.
At present, the research on the water invasion influence of the side top water reservoir is mainly numerical simulation research, and the physical simulation is mainly based on the existing bottom water reservoir three-dimensional simulation rotation model. The three-dimensional physical simulation is a common means for researching water invasion of the bottom water reservoir, and a three-dimensional physical simulation experiment can artificially reproduce the development process of the reservoir and simulate the seepage characteristics of the actual reservoir and an oil well more truly. The three-dimensional physical simulation system of the bottom water reservoir is disclosed as the following patent: CN10516102A, the same as the three-dimensional physical model in the bottom water reservoir water energy three-dimensional physical simulation device and method, generally comprises a model body, a production metering system, a detection system, a control system and a bottom water simulation system. On the basis of the bottom water reservoir three-dimensional physical simulation system, the model rotates, and three-dimensional physical simulation of top water heavy oil reservoir water invasion can be realized.
However, the method has certain limitation in the process of performing top water invasion simulation by bottom water reservoir three-dimensional simulation, and is suitable for the bottom water reservoir with sufficient water energy. The problems that the temperature and pressure conditions which are the same as those of an actual steam cavity are not easy to construct, the water invasion amount is not well controlled and the like exist; in practical application, the relation between water invasion amount and steam cavity temperature and pressure is difficult to quantify, the simulation experiment is long in operation time, various equivalent treatments of oil deposit data are often required, the related physical parameters are more, the operation is very troublesome, the time consumption is long, and the requirement for rapidly realizing the SAGD edge top water invasion physical simulation indoors cannot be met.
Disclosure of Invention
In view of the above-mentioned problem that the conventional three-dimensional physical simulation rotating model of the bottom water reservoir is not easy to construct the same temperature and pressure conditions as the actual steam cavity, the present invention is proposed to provide a method and a system for constructing a three-dimensional physical model of top water invasion, which overcome or at least partially solve the above-mentioned problems.
According to one aspect of the invention, a method for constructing a three-dimensional physical model of edge-top water invasion is provided, and the method comprises the following steps:
injecting a first water body into the reaction vessel; the reaction vessel is a vessel for simulating a steam cavity formed by an oil reservoir in the SAGD production process;
heating the reaction vessel containing the first body of water to a first temperature to simulate steam chamber conditions; the vapor chamber conditions include: temperature, pressure and gas-liquid ratio of the vapor chamber; the gas-liquid ratio is the ratio of the gas volume in the steam cavity to the water volume;
and simulating a top water invasion process by injecting a second water body into the reaction container meeting the steam cavity condition to construct a top water invasion three-dimensional physical model.
Preferably, simulating steam chamber conditions comprises:
acquiring a gas-liquid ratio of the steam cavity;
controlling the vapor chamber conditions inside the reaction vessel by controlling a first pressure, a first temperature, and a gas-to-liquid ratio of the reaction vessel.
Preferably, obtaining the gas-liquid ratio of the vapor chamber comprises:
the formula I and the formula II are combined to obtain a formula III:
wherein q is o In order to produce oil, L is the length of a steam injection section of the horizontal well, k is the effective permeability of an oil phase, phi is porosity, alpha is thermal diffusion coefficient, g is gravity acceleration, S oi And S ors Respectively representing the initial oil saturation of the oil reservoir and the residual saturation after steam injection, h is the thickness of the oil layer, y is the height of the steam cavity, m is the viscosity constant, v o The kinematic viscosity of the crude oil, x is the advancing distance of a steam front edge in the horizontal direction, and t is the injection time of the steam; wherein, Delta S o ＝S oi －S ors ；
Integrating the steam generation time t in equation three to obtain equation four:
acquiring a formula five according to the formula four to obtain the height y of the steam cavity;
correcting the shape curve of the steam cavity and obtaining the corrected volume of the steam cavity;
determining the gas height and the water height of the steam cavity;
and acquiring the gas-liquid ratio of the steam cavity based on the gas height, the water height and the volume of the steam cavity.
Preferably, the step of correcting the shape curve of the steam cavity and obtaining the corrected volume of the steam cavity further comprises:
placing a shape curve of the steam cavity in a coordinate system;
making a tangent line passing through the origin of the coordinate system on the shape curve to obtain a tangent point on the shape curve of the steam cavity;
obtaining the tangent point (x) of the tangent line according to a formula six and a formula five 0 ，y 0 ) Wherein k is x Is the slope of the tangent:
obtaining a horizontal line passing through the tangent point based on the tangent point, and dividing the shape curve into a first part and a second part through the horizontal line;
respectively obtaining the area of the first part and the area of the second part according to the tangent points;
obtaining half of a side area of the steam cavity, the half of the side area being a sum of an area of the first portion and an area of the second portion;
and obtaining the volume of the steam cavity according to the side area and the length of the steam injection section of the horizontal well.
Preferably, the obtaining the gas-liquid ratio of the steam cavity based on the gas height, the water height and the volume of the steam cavity comprises:
obtaining the water volume V of the steam cavity based on formula seven w ：
Wherein h is s Is the gas height of the steam chamber, h-h s The water height of the steam cavity;
acquiring the gas volume of the steam cavity according to the difference between the volume of the steam cavity and the volume of the water body;
and acquiring a gas-liquid ratio of the steam cavity, wherein the gas-liquid ratio of the steam cavity is the ratio of the gas volume to the water volume.
Preferably, the method further comprises: determining the volume of the first body of water injected according to the first temperature and the gas-liquid ratio.
Preferably, the method further comprises:
obtaining the volume V of the injected first water body according to the formula W0 ：
Wherein, V t Is the volume of the reaction vessel, R sw Gas-liquid ratio, p s1 、ρ w1 Respectively water density, gas density, rho at corresponding temperatures w0 Is the density of the injected first body of water.
Preferably, the method further comprises:
injecting a second body of water into the reaction vessel having a first temperature further comprises: and injecting the second water body at a constant speed through the constant pressure pump.
According to another aspect of the present invention, there is provided a side-roof water invasion three-dimensional physical simulation system, comprising:
the first water injection unit is used for injecting a first water into the reaction vessel; the reaction vessel is a vessel for simulating a steam cavity formed by an oil reservoir in the SAGD production process;
a heating unit for heating the reaction vessel containing the first body of water to a first temperature to simulate steam chamber conditions; the vapor chamber conditions include: temperature, pressure and gas-liquid ratio of the vapor chamber; the gas-liquid ratio is the ratio of the gas volume in the steam cavity to the water volume;
and the second water body injection unit is used for simulating the top water invasion process by injecting second water bodies into the reaction container meeting the steam cavity conditions so as to construct a top water invasion three-dimensional physical model.
According to an aspect of the present invention, there is provided a computing apparatus, including a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor executes the computer program to implement the method for building a three-dimensional physical model of edge-top water invasion as described in any one of the above.
The method changes the thinking, adopts the independent steam cavity, simulates the side top water invasion process by injecting a second water body, and establishes the implementation method of the physical model of the process after the side top water intrudes into the steam cavity.
The foregoing description is only an overview of the technical solutions of the present invention, and the embodiments of the present invention are described below in order to make the technical means of the present invention more clearly understood and to make the above and other objects, features, and advantages of the present invention more clearly understandable.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the embodiments or the prior art descriptions will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
FIG. 1 is a flow chart of a method for constructing a three-dimensional physical model of side-roof water invasion in the embodiment of the invention;
FIG. 2 is a schematic view of a modification of the shape curve of the vapor chamber in an embodiment of the present invention;
FIG. 3A, FIG. 3B and FIG. 3C are the relationship curves of the logarithm value of the temperature variation with the time t, which are obtained by adding different second water bodies at 200 ℃ to carry out the water invasion test respectively;
FIG. 4A, FIG. 4B and FIG. 4C are the relationship curves of the logarithm value of the temperature variation with the time t, which are obtained by adding different second water bodies at 250 ℃ to carry out the water invasion test respectively;
FIGS. 5A, 5B, and 5C are graphs of the relationship between gas volume, temperature, and pressure, respectively, and water intrusion;
FIG. 6 is a schematic structural diagram of a side-roof water invasion three-dimensional physical model construction system in an embodiment of the present invention;
FIG. 7 is a diagram illustrating a computing device in accordance with an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The embodiment of the invention provides a method for constructing a three-dimensional physical model of side-roof water invasion, which comprises the following steps of:
103, simulating a top water invasion process by injecting a second water body into the reaction container meeting the steam cavity condition to construct a top water invasion three-dimensional physical model. In a specific embodiment, after the simulated steam cavity is constructed in step 102, a process of side top water invasion can be simulated through step 103, that is, the construction of a side top water invasion three-dimensional physical model is realized. And in the water invasion process, the influence of side top water invasion on the SAGD development of the side top water heavy oil reservoir is quickly predicted by measuring the changes of the gas volume, the temperature and the pressure, so that guidance is provided for the development of the side top water heavy oil reservoir.
The method changes the thinking, adopts the method that the steam cavity is independent, the same steam cavity condition is realized by heating water through the high-temperature high-pressure reaction kettle, the top water invasion process is simulated by injecting the second water body, and the physical model of the process after the top water intrudes into the steam cavity is established.
According to the method for constructing the side-top water invasion three-dimensional physical model, the conditions of the steam cavity are preferably simulated by the following steps:
acquiring a gas-liquid ratio of the steam cavity;
controlling the vapor chamber conditions inside the reaction vessel by controlling a first pressure, a first temperature, and a gas-to-liquid ratio of the reaction vessel.
In a specific embodiment, the steam cavity conditions mainly include pressure, temperature, and a gas-liquid ratio, wherein the gas-liquid ratio is mainly used for determining the volume injected into the first water body, so that when the reaction vessel meets the steam cavity conditions, the simulated steam cavity configuration can be considered to be completed. Wherein, the pressure condition in the steam chamber condition can satisfy through the size of adjusting first pressure, and the temperature condition can satisfy to first temperature through heating reaction vessel, that is to say that after first pressure, first temperature and the three condition of gas-liquid ratio satisfied, steam chamber in the simulation SAGD production process that can be quick need not to reconstruct the model of whole oil reservoir.
In preferred embodiments, indirect control of the first pressure is also achieved by direct control of the gas-to-liquid ratio and the first temperature.
According to the method for constructing the side-roof water invasion three-dimensional physical model provided by the embodiment of the invention, preferably, the step of obtaining the gas-liquid ratio of the steam cavity comprises the following steps:
the formula I and the formula II are combined to obtain a formula III:
wherein q is o For oil production, m 3 S; l is the length of the steam injection section of the horizontal well, m; k is the effective permeability of the oil phase, m 2 (ii) a Phi is porosity, dimensionless; α is the thermal diffusion coefficient, m 2 S; g is the acceleration of gravity, m/s 2 ；S oi And S ors Respectively representing the initial oil saturation and injected steam of the reservoirThe residual saturation after, h is the oil layer thickness, m; y is the steam cavity height, m; m is viscosity constant, v o Is the kinematic viscosity of the crude oil, m 2 S; x is the propulsion distance of the steam front edge in the horizontal direction, m; t is the injection time of steam; and s. Wherein, Delta S o ＝S oi －S ors . Preferably, m is 3.
Specifically, the formula I and the formula II are both SAGD yield calculation formulas, and the formula I and the formula II are combined to obtain the formula III.
Integrating the generation time t of steam in the formula three, and when t is 0, x is 0 to obtain the formula four:
changing the formula IV to obtain a formula V to obtain the height y of the steam cavity;
and correcting the shape curve of the steam cavity and obtaining the corrected volume of the steam cavity. Specifically, from the fifth formula, at the bottom of the steam cavity, in the SAGD production process, with the continuous injection of steam, the bottom of the steam cavity is farther and farther from the steam injection well above, which is not in accordance with the reality. The main reason for this is the assumption of steady-state accumulation of heat used in the Butler model. It is therefore necessary to modify the shape curve of the steam chamber and obtain a modified volume of said steam chamber.
Determining the gas height and the water height of the steam cavity;
and acquiring the gas-liquid ratio of the steam cavity based on the gas height, the water height and the volume of the steam cavity. Specifically, the gas-liquid ratio of the steam cavity is the ratio of the water volume to the gas volume of the steam cavity, and the ratio of the water volume to the gas volume is related to the gas height and the water height, so that the gas-liquid ratio of the steam cavity can be obtained through the gas height, the water height and the volume of the steam cavity.
In the method for constructing a three-dimensional physical model of side-topped water invasion according to the embodiment of the present invention, preferably, the step of correcting the shape curve of the steam cavity and obtaining the corrected volume of the steam cavity further includes:
placing the shape curve of the steam cavity in a coordinate system.
Specifically, the curve in fig. 2 is a half of the shape curve of the steam cavity, i.e. represents a half of the steam cavity, i.e. the volume of the steam cavity can be obtained as long as a half of the volume of the steam cavity is obtained.
Making a tangent line passing through the origin of the coordinate system on the shape curve to obtain a tangent point on the shape curve of the steam cavity; as shown in fig. 2, line OA is a tangent line passing through the origin of the coordinate system and tangent to the shape curve of the vapor chamber, and a is the tangent point.
Obtaining the tangent point A (x) of the tangent line according to a formula six and a formula five 0 ，y 0 ) Wherein k is x Is the slope of the tangent:
wherein y is the height of the steam cavity, and h is the oil layer degree along with the change of the advancing distance x of the steam front edge in the horizontal direction, namely the maximum position of the steam after the steam is stabilized.
As shown in fig. 2, a horizontal line AC passing through the tangent point a is obtained based on the tangent point, and the shape curve is divided into a first part and a second part by the horizontal line AC. Specifically, the first portion is surrounded by the OAC, that is, the second portion is surrounded by the CABD.
And respectively obtaining the area of the first part and the area of the second part according to the tangent points.
Specifically, point of tangency A (x) 0 ，y 0 ) Substituting into equation six, we can solve the tangent point A (x) 0 ，y 0 ) And slope k x Respectively as follows:
further, the area S1 of the first portion and the area S2 of the second portion can be obtained by the following formulas:
obtaining half of a side area of the steam cavity, the half of the side area being a sum of an area of the first portion and an area of the second portion. Specifically, after obtaining the area S1 of the first portion and the area S2 of the second portion, the available side area S of the half steam chamber is:
accordingly, the side area of the entire steam chamber is 2S.
And obtaining the volume of the steam cavity according to the side area and the length of the steam injection section of the horizontal well.
Specifically, the volume of the steam cavity can be obtained by multiplying the lateral area of the steam cavity by the length of the steam injection section of the horizontal well. That is, the volume V of the steam cavity is obtained by the following formula 0 ：
In the method for constructing a three-dimensional physical model of water invasion by edge top water according to the embodiment of the present invention, preferably, the step of obtaining the gas-liquid ratio of the steam cavity based on the gas height, the water height and the volume of the steam cavity comprises:
obtaining the water volume V of the steam cavity based on formula seven w ：
Wherein h is s Is the gas height of the steam chamber, h-h s The water height of the steam cavity.
Acquiring the gas volume of the steam cavity according to the difference between the volume of the steam cavity and the volume of the water body;
that is, the gas volume V is obtained by the following formula s ：
And acquiring a gas-liquid ratio of the steam cavity, wherein the gas-liquid ratio of the steam cavity is the ratio of the gas volume to the water volume.
Specifically, the gas-liquid ratio R of the vapor chamber can be obtained by the following formula sw ：
The method for constructing the three-dimensional physical model of side-roof water invasion, provided by the embodiment of the invention, preferably comprises the following steps: determining the volume of the first body of water injected according to the first temperature and the gas-liquid ratio.
The method for constructing the edge top water invasion three-dimensional physical model in the embodiment of the invention preferably further comprises the following steps:
obtained according to the formula eightVolume V of the injected first body of water W0 ：
Wherein, V t Is the volume of the reaction vessel, R sw Is the gas-liquid ratio, p s1 、ρ w1 Respectively water density, gas density, rho at corresponding temperatures w0 Is the density of the injected first body of water.
In a specific embodiment, the interior of the reaction vessel satisfies the energy conservation and the volume conservation, so that the following can be obtained according to the energy conservation and the volume conservation:
in the above formula, the volume of the reaction vessel is V t Volume of gas V under saturated vapor chamber conditions S1 Volume of liquid V w1 . Solving the above equation to obtain equation eight, i.e. determining the volume V of the first water body W0 。
The method for constructing the edge top water invasion three-dimensional physical model in the embodiment of the invention preferably further comprises the following steps:
injecting a second body of water into the reaction vessel having a first temperature further comprises: and injecting the second water body at a constant speed through the constant pressure pump. Specifically, the constant pressure pump can inject the second water body into the reaction container for a short time at a constant speed within a certain time to simulate the side top water invasion, and the volume of the second water body can be conveniently calculated through the constant speed injection. The temperature and pressure changes in the reaction vessel can then be observed to simulate those that occur when the steam chamber encounters the edge water.
In a preferred embodiment, the method further comprises:
the heat loss of the model is corrected by the following verification method.
Specifically, the room temperature condition is assumed to be the temperature T S0 Pressure P S0 (ii) a The saturated air cavity condition is temperature T S1 Pressure P of S1 (ii) a The temperature and pressure of the steam cavity under the thermo-dynamic effect of temperature difference are: t is s2 ，P s2 (ii) a Inner cavity gas volume becomes V s2 (ii) a Volume of water V w2 (ii) a Volume of water invasion V wi 。
According to the conservation of energy: i.e. the original energy + the implantation energy is the final energy available
H w1 (T s1 )×ρ w1 (T s1 )×V w1 +H s1 (T s1 )×ρ s1 (T s1 )×V s1 +H w0 (T s0 )×ρ w0 (T s0 )×V wi ＝H w2 (T s2 )×ρ w2 (T s2 )×V w2 +H s2 (T s2 )×ρ s2 (T s2 )×V s2
In the formula, H s ，H w The enthalpy value of the water body and the enthalpy value of the gas at the corresponding temperature are KJ/Kg.
According to the conservation of volume, i.e. the volume of gas phase + the volume of liquid phase being the volume of inner cavity
V s2 +V w2 ＝V t
According to conservation of mass, i.e. original mass + implant mass being the final mass available
ρ w1 (T s1 )×V w1 +ρ s1 (T s1 )×V s1 +ρ w0 (T s0 )×V wi ＝ρ w2 (T s2 )×V w2 +ρ s2 (T s2 )×V s2
Solving the relationship with respect to temperature as:
two constants are defined:
K 1 ＝H w1 (T s1 )×ρ w1 (T s1 )×V w1 +H s1 (T s1 )×ρ s1 (T s1 )×V s1 +H w0 (T s0 )×ρ w0 (T s0 )×V wi
K 2 ＝ρ w1 (T s1 )×V w1 +ρ s1 (T s1 )×V s1 +ρ w0 (T s0 )×V wi
defining a function f (T) s2 ) Wherein, T s2 At T S0 And T S1 And (3) solving the root by a dichotomy:
f(T S2 )＝(K 1 -H s2 (T s2 )×ρ s2 (T s2 )×V t )(ρ w2 (T s2 )-ρ s2 (T s2 ))-(K 2 -ρ s2 (T s2 )×V t )(H w2 (T s2 )×ρ w2 (T s2 )-H s2 (T s2 )×ρ s2 (T s2 ))
f (T) is obtained by judgment S2 ) The monotonicity of the function can be solved by a dichotomy:
(1) defining interval as [ a, b]，a＝T s0 ，b＝T s1 Setting precision xi;
(2) solving a middle point c of the interval;
(3) calculating f (c):
if f (c) is 0, c is the zero point of the function;
if f (c) is less than 0, making b ═ c;
if f (c) is 0, let a be c;
(4) judging whether the precision xi is achieved, namely if | a-b | < xi, obtaining a zero point approximate value a (or b), and otherwise, repeating (2) - (4);
from temperature T s2 Calculating the gas volume and the water volume after water invasion:
the saturation pressure P is obtained from the volume of the water body S2 And heat loss is corrected.
The high-temperature high-pressure reaction kettle is used as a reaction container to carry out multi-temperature and multi-proportion water invasion experiments, two groups of results of 200 ℃ and 250 ℃ are compared, and the change curve of the temperature in the reaction container along with the time is measured by the experiments to obtain (T-T) f ) The logarithmic value of (a) with time t is shown in fig. 3A, fig. 3B, fig. 3C, fig. 4A, fig. 4B, and fig. 4C, wherein the initial conditions of fig. 3A, fig. 3B, and fig. 3C are the logarithmic value with time t of the temperature change obtained by heating a first water body of 20 ℃, 100mL to 200 ℃, and then adding 10mL, 20mL, and 30mL of second water bodies of 20 ℃ respectively for water invasion test; fig. 4A, fig. 4B, and fig. 4C are graphs in which initial conditions are curves of logarithmic values of temperature changes with time t obtained by heating a first water body of 20 ℃ and 100mL to 250 ℃, and then adding 10mL, 20mL, and 30mL of second water bodies of 20 ℃ respectively to perform a water invasion test, wherein an intersection point of a reverse extension line of the curves and a longitudinal axis is an initial temperature. R 2 The fitting degree of the curve with the curve obtained by the model of the embodiment is actually verified.
On the basis, the influence of water invasion on physical parameters of the steam cavity is calculated and analyzed by using the established steam cavity temperature difference thermodynamics and water invasion mathematical verification method, so that a relation graph between gas volume, temperature, pressure and water invasion in the steam cavity is obtained, and the relation graph is shown in fig. 5A, 5B and 5C. Wherein, in FIG. 5A, FIG. 5B and FIG. 5C, the gas-liquid ratio in the heated reaction vessel is 54, and the temperature of the added second water body is 45 ℃.
In fig. 5A, the gas volume is shown as a function of water intrusion, where the vertical axis is the gas volume percentage in the steam chamber, initially 1, and the corresponding gas volume percentage decreases with increasing second body of water. Fig. 5B shows the temperature versus water intrusion, with the temperature of the steam chamber decreasing as the water intrusion increases. Fig. 5C shows the pressure of the steam chamber as a function of water intrusion, the pressure of the steam chamber decreasing with increasing second body of water.
An embodiment of the present invention further provides a side-topped water invasion three-dimensional physical simulation system, as shown in fig. 6, the system includes:
a first water injection unit 601 for injecting a first water into the reaction vessel; the reaction container is a container for simulating a steam cavity formed on the basis of an oil reservoir in the SAGD production process;
a heating unit 602 for heating the reaction vessel containing the first body of water to a first temperature to simulate steam chamber conditions; the vapor chamber conditions include: temperature, pressure and gas-liquid ratio of the vapor chamber; the gas-liquid ratio is the ratio of the gas volume in the steam cavity to the water volume;
and a second water injection unit 603, configured to simulate a top water invasion process by injecting a second water into the reaction vessel satisfying the steam cavity condition, so as to construct a top water invasion three-dimensional physical model.
The embodiment of the present invention further provides a computing device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, where the processor executes the computer program to implement the method for building a three-dimensional physical model of top water invasion as described in any of the above embodiments.
As shown in fig. 7, which is a block diagram of a computing device in embodiments herein, as shown in fig. 7, the computing device 702 may include one or more processing devices 704, such as one or more Central Processing Units (CPUs), each of which may implement one or more hardware threads. The computing device 702 may also include any storage resources 706 for storing any kind of information, such as code, settings, data, and the like. For example, and without limitation, the storage resources 706 may include any one or more of the following in combination: any type of RAM, any type of ROM, flash memory devices, hard disks, optical disks, etc. More generally, any storage resource may use any technology to store information. Further, any storage resource may provide volatile or non-volatile reservation of information. Further, any storage resource may represent a fixed or removable component of computing device 702. In one case, when processing device 704 executes associated instructions that are stored in any storage resource or combination of storage resources, computing device 702 can perform any of the operations of the associated instructions. The computing device 702 also includes one or more drive mechanisms 708, such as a hard disk drive mechanism, an optical disk drive mechanism, or the like, for interacting with any storage resource.
It should be understood that, in various embodiments of the present invention, the sequence numbers of the above-mentioned processes do not mean the execution sequence, and the execution sequence of each process should be determined by its function and inherent logic, and should not constitute any limitation on the implementation process of the embodiments of the present invention.
It should also be understood that, in the embodiment of the present invention, the term "and/or" is only one kind of association relation describing an associated object, and means that three kinds of relations may exist. For example, a and/or B, may represent: a exists alone, A and B exist simultaneously, and B exists alone. In addition, the character "/" herein generally indicates that the former and latter associated objects are in an "or" relationship.
Those of ordinary skill in the art will appreciate that the various illustrative components and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, computer software, or combinations of both, and that the components and steps of the various examples have been described above generally in terms of their functionality in order to clearly illustrate this interchangeability of hardware and software. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the technical solution. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
It is clear to those skilled in the art that, for convenience and brevity of description, the specific working processes of the above-described systems, apparatuses and units may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
In the several embodiments provided in the present application, it should be understood that the disclosed system, apparatus and method may be implemented in other ways. For example, the above-described apparatus embodiments are merely illustrative, and for example, the division of the units is only one logical division, and other divisions may be realized in practice, for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may also be an electric, mechanical or other form of connection.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment of the present invention.
In addition, functional units in the embodiments of the present invention may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit. The integrated unit may be implemented in the form of hardware, or may also be implemented in the form of a software functional unit.
The integrated unit, if implemented in the form of a software functional unit and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present invention essentially or partially contributes to the prior art, or all or part of the technical solution can be embodied in the form of a software product stored in a storage medium and including instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
The principle and the implementation mode of the invention are explained by applying specific embodiments in the invention, and the description of the embodiments is only used for helping to understand the method and the core idea of the invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed, and in summary, the content of the present specification should not be construed as a limitation to the present invention.
Claims (8)
1. A method for constructing a three-dimensional physical model of edge top water invasion, which is characterized by comprising the following steps:
injecting a first water body into the reaction vessel; the reaction container is a container for simulating a steam cavity formed on the basis of an oil reservoir in the SAGD production process;
heating the reaction vessel containing the first body of water to a first temperature to simulate steam chamber conditions; the vapor chamber conditions include: temperature, pressure and gas-liquid ratio of the vapor chamber; the gas-liquid ratio is the ratio of the gas volume in the steam cavity to the water volume; simulating a top water invasion process by injecting a second water body into the reaction vessel meeting the steam cavity conditions to construct a top water invasion three-dimensional physical model; wherein simulating vapor chamber conditions comprises:
acquiring a gas-liquid ratio of the steam cavity;
controlling the interior of the reaction vessel to achieve the vapor cavity conditions by controlling a first pressure, a first temperature, and a gas-to-liquid ratio of the reaction vessel;
wherein obtaining a gas-to-liquid ratio of the vapor chamber comprises:
the formula I and the formula II are combined to obtain a formula III:
wherein q is o In order to produce oil, L is the length of a steam injection section of the horizontal well, k is the effective permeability of an oil phase, phi is porosity, alpha is thermal diffusion coefficient, g is gravity acceleration, S oi And S ors Respectively representing the initial oil saturation of the oil reservoir and the residual saturation after steam injection, h is the oil layer thickness, y is the steam cavity height, m is the viscosity constant, v o The kinematic viscosity of the crude oil, x is the advancing distance of a steam front edge in the horizontal direction, and t is the injection time of the steam; wherein, Delta S o ＝S oi －S ors ；
Integrating the generation time t of steam in equation three to obtain equation four:
acquiring a formula five according to the formula four to obtain the height y of the steam cavity;
correcting the shape curve of the steam cavity and obtaining the corrected volume of the steam cavity;
determining the gas height and the water height of the steam cavity;
and acquiring the gas-liquid ratio of the steam cavity based on the gas height, the water height and the volume of the steam cavity.
2. The method for constructing the three-dimensional physical model of water invasion by edge top water according to claim 1, wherein the step of modifying the shape curve of the steam chamber and obtaining the modified volume of the steam chamber further comprises:
placing the shape curve of the steam cavity in a coordinate system;
making a tangent line passing through the origin of the coordinate system on the shape curve to obtain a tangent point on the shape curve of the steam cavity;
obtaining the tangent point (x) of the tangent line according to a formula six and a formula five 0 ，y 0 ) Wherein k is x Is the slope of the tangent:
obtaining a horizontal line passing through the tangent point based on the tangent point, and dividing the shape curve into a first part and a second part through the horizontal line;
respectively obtaining the area of the first part and the area of the second part according to the tangent points;
obtaining half of a side area of the steam cavity, the half of the side area being a sum of an area of the first portion and an area of the second portion;
and obtaining the volume of the steam cavity according to the side area and the length of the steam injection section of the horizontal well.
3. The method for constructing the three-dimensional physical model of water invasion in the top of the edge of claim 2, wherein the step of obtaining the gas-liquid ratio of the steam cavity based on the gas height, the water height and the volume of the steam cavity comprises the steps of:
obtaining the water volume V of the steam cavity based on formula seven w ：
Wherein h is s Is the gas height of the steam chamber, h-h s The water height of the steam cavity;
acquiring the gas volume of the steam cavity according to the difference between the volume of the steam cavity and the volume of the water body;
and acquiring a gas-liquid ratio of the steam cavity, wherein the gas-liquid ratio of the steam cavity is the ratio of the gas volume to the water volume.
4. The method for constructing the three-dimensional physical model of water invasion of the top of edge water according to claim 3, wherein the method further comprises the following steps: determining the volume of the first body of water injected according to the first temperature and the gas-liquid ratio.
5. The method for constructing the three-dimensional physical model of water invasion of the top of edge water according to claim 4, wherein the method further comprises the following steps:
obtaining the volume V of the injected first water body according to the formula W0 ：
Wherein, V t Is the volume of the reaction vessel, R sw Is the gas-liquid ratio, p s1 、ρ w1 Respectively water density, gas density, rho at corresponding temperatures w0 Is the density of the first body of water injected.
6. The method for constructing the three-dimensional physical model of water invasion of the top of edge of claim 1, wherein the method further comprises the following steps:
injecting a second body of water into the reaction vessel having a first temperature further comprises: and injecting the second water body at a constant speed through the constant pressure pump.
7. A three-dimensional physical simulation system for edge-top water intrusion, the system comprising:
the first water injection unit is used for injecting a first water into the reaction vessel; the reaction vessel is a vessel for simulating a steam cavity formed by an oil reservoir in the SAGD production process;
a heating unit for heating the reaction vessel containing the first body of water to a first temperature to simulate steam chamber conditions; the vapor chamber conditions include: temperature, pressure and gas-liquid ratio of the vapor chamber; the gas-liquid ratio is the ratio of the gas volume in the steam cavity to the water volume;
the second water body injection unit is used for simulating a top water invasion process by injecting second water bodies into the reaction container meeting the steam cavity conditions so as to construct a top water invasion three-dimensional physical model;
the simulated vapor chamber conditions include:
acquiring a gas-liquid ratio of the steam cavity;
controlling the interior of the reaction vessel to achieve the vapor cavity conditions by controlling a first pressure, a first temperature, and a gas-to-liquid ratio of the reaction vessel; wherein the content of the first and second substances,
obtaining a gas-liquid ratio of the vapor chamber comprises:
the formula I and the formula II are combined to obtain a formula III:
wherein q is o In order to produce oil, L is the length of a steam injection section of the horizontal well, k is the effective permeability of an oil phase, phi is porosity, alpha is thermal diffusion coefficient, g is gravity acceleration, and S is oi And S ors Respectively representing the initial oil saturation of the oil reservoir and the residual saturation after steam injection, h is the thickness of the oil layer, y is the height of the steam cavity, m is the viscosity constant, v o The kinematic viscosity of the crude oil, x is the advancing distance of a steam front edge in the horizontal direction, and t is the injection time of the steam; wherein, Delta S o ＝S oi －S ors ；
Integrating the generation time t of steam in equation three to obtain equation four:
acquiring a formula five according to the formula four to obtain the height y of the steam cavity;
correcting the shape curve of the steam cavity and obtaining the corrected volume of the steam cavity;
determining the gas height and the water height of the steam cavity;
and acquiring the gas-liquid ratio of the steam cavity based on the gas height, the water height and the volume of the steam cavity.
8. A computing device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of building a three-dimensional physical model of water intrusion into edge top water of any one of claims 1 to 6 when executing the computer program.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110598395.XA CN113283092B (en) | 2021-05-31 | 2021-05-31 | Method and system for constructing side-roof water invasion three-dimensional physical model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110598395.XA CN113283092B (en) | 2021-05-31 | 2021-05-31 | Method and system for constructing side-roof water invasion three-dimensional physical model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113283092A CN113283092A (en) | 2021-08-20 |
| CN113283092B true CN113283092B (en) | 2022-09-13 |
Family
ID=77282529
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110598395.XA Active CN113283092B (en) | 2021-05-31 | 2021-05-31 | Method and system for constructing side-roof water invasion three-dimensional physical model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113283092B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105156102A (en) * | 2015-09-28 | 2015-12-16 | 中国石油大学(北京) | Water body energy three-dimensional physical simulation device and method of bottom water reservoir |
| CN111075441A (en) * | 2019-12-24 | 2020-04-28 | 中国石油化工股份有限公司 | Three-dimensional physical simulation experiment device and method for cold recovery after thermal recovery of side-bottom water heavy oil reservoir |
-
2021
- 2021-05-31 CN CN202110598395.XA patent/CN113283092B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105156102A (en) * | 2015-09-28 | 2015-12-16 | 中国石油大学(北京) | Water body energy three-dimensional physical simulation device and method of bottom water reservoir |
| CN111075441A (en) * | 2019-12-24 | 2020-04-28 | 中国石油化工股份有限公司 | Three-dimensional physical simulation experiment device and method for cold recovery after thermal recovery of side-bottom water heavy oil reservoir |
Non-Patent Citations (1)
| Title |
|---|
| 高压环境双水平井SAGD三维物理模拟实验;田杰等;《石油学报》;20170425(第04期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113283092A (en) | 2021-08-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111236931B (en) | Method and system for generating unsteady well testing chart of gas reservoir highly-deviated well | |
| CN105156081B (en) | A kind of carbonate rock heavy crude reservoir is acidified simulation evaluation method | |
| CN110472372B (en) | Dual-medium-based permeability prediction method and system | |
| CN106442271A (en) | Rock core permeability simulating method and rock core permeability simulating device | |
| CN110984973A (en) | Determination method for single-well controlled reserve of fracture-cave carbonate gas reservoir | |
| Ge et al. | A novel empirical equation for relative permeability in low permeability reservoirs | |
| CN115828704A (en) | Rapid prediction method for underground water pollution | |
| Liu et al. | A simplified and efficient method for water flooding production index calculations in low permeable fractured reservoir | |
| CN113283092B (en) | Method and system for constructing side-roof water invasion three-dimensional physical model | |
| CN115526086A (en) | Carbonate reservoir water drive breakthrough time prediction method based on proxy model | |
| CN116127859A (en) | High-temperature dynamic water grouting simulation method and system for fractured rock | |
| Pan et al. | Data-Driven Deconvolution Using Echo-State Networks Enhances Production Data Analysis in Unconventional Reservoirs | |
| CN105765161A (en) | Optimizing flow control device properties on injector wells in liquid flooding systems | |
| Ivanov et al. | Acid Treatment Optimization Based on Digital Core Analysis | |
| Singh et al. | Pore-scale Ostwald ripening of gas bubbles in the presence of oil and water in porous media | |
| Meng et al. | Preconditioned multiple-relaxation-time lattice Boltzmann equation model for incompressible flow in porous media | |
| Motta et al. | Effects of converging-diverging pore geometry on the acidizing process with non-Newtonian Carreau-type fluids | |
| CN105899755A (en) | Optimizing flow control device properties on a producer well in coupled injector-producer liquid flooding systems | |
| CN116167302A (en) | Description method of artificial complex cracks in natural gas hydrate yield increase simulation | |
| CN116933670A (en) | Method and system for simulating interaction of fluid rock | |
| Li et al. | Phase-field modeling of coupled reactive transport and pore structure evolution due to mineral dissolution in porous media | |
| CN114153007B (en) | Modeling method for water lock effect damage oil-gas layer, damage degree spatial-temporal evolution 4D quantitative and intelligent diagnosis method and system thereof | |
| CN113887110A (en) | Oil well productivity prediction method and device, electronic equipment and computer storage medium | |
| CN107133373B (en) | Coupling simulation method for shale gas reservoir, shaft and ground pipe network | |
| CN108959764B (en) | Low-permeability reservoir water-drive oil-sweeping area coefficient evaluation method and system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |