CN114528780A - Numerical simulation method considering non-Newtonian characteristics of thickened oil - Google Patents

Numerical simulation method considering non-Newtonian characteristics of thickened oil Download PDF

Info

Publication number
CN114528780A
CN114528780A CN202210152007.XA CN202210152007A CN114528780A CN 114528780 A CN114528780 A CN 114528780A CN 202210152007 A CN202210152007 A CN 202210152007A CN 114528780 A CN114528780 A CN 114528780A
Authority
CN
China
Prior art keywords
oil
numerical simulation
flow
viscosity
thickened
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.)
Pending
Application number
CN202210152007.XA
Other languages
Chinese (zh)
Inventor
姜宇
陈秀华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Qingdao Shida Hengrui Petroleum New Technology Development Co ltd
Original Assignee
Qingdao Shida Hengrui Petroleum New Technology Development Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Qingdao Shida Hengrui Petroleum New Technology Development Co ltd filed Critical Qingdao Shida Hengrui Petroleum New Technology Development Co ltd
Priority to CN202210152007.XA priority Critical patent/CN114528780A/en
Publication of CN114528780A publication Critical patent/CN114528780A/en
Pending legal-status Critical Current

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
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • 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)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a numerical simulation method considering the non-Newtonian characteristic of cold recovery of thickened oil, which specifically comprises the following steps of: inputting grid geological parameters; inputting rock and fluid parameters; establishing a numerical model and realizing numerical simulation software; and outputting and displaying the calculation result. The invention provides two methods for representing the flow rule of thickened oil: characterizing the flow rule of the thickened oil based on the phenomenon that the viscosity of the thickened oil is shear-thinned; and characterizing the thick oil flowing rule based on the nonlinear seepage phenomenon of the thick oil flowing in the displacement experiment. The method can quantitatively represent the physical phenomenon that the viscosity of the heavy oil reservoir changes along with displacement and the nonlinear seepage rule thereof, and has high numerical simulation precision and stability; the accuracy of history fitting in the numerical simulation process is obviously improved, the phenomena of non-Newtonian fluid characteristics and non-linear seepage in the heavy oil reservoir are reasonably described, and theoretical and technical supports are provided for guiding the next step of regulation and control of the heavy oil reservoir.

Description

Numerical simulation method considering non-Newtonian characteristics of thickened oil
Technical Field
The invention relates to the technical field of numerical reservoir simulation, in particular to a numerical simulation method considering non-Newtonian characteristics of heavy oil.
Background
The thickened oil is a mixture rich in high carbon number, colloid and asphaltene and better accords with non-Newtonian fluid from the flow angle, but the existing thickened oil underground flow almost supposes to meet Darcy seepage flow, and based on Darcy seepage flow theory, the production capacity of hundreds of thousands of thickened oil reservoirs either has no natural production capacity or ultra-low natural production capacity, which is obviously different from the actual production situation of the thickened oil. And water is injected for water-cooling heavy oil recovery, Darcy flow is assumed, and water-containing abnormity breaks through due to huge oil-water fluidity difference, so that the evaluations have great contradiction with the mine field practice.
Therefore, it is urgently needed to improve the simulation precision of the heavy oil reservoir numerical simulation technology on the basis of the knowledge of the reservoir seepage mechanism. Through experimental research, the increase of the flow speed or the shear rate in the process of developing the heavy oil reservoir can cause the reduction of the viscosity of the heavy oil, and the phenomenon that the viscosity is shear-thinned occurs; in a displacement experiment, the thick oil flow is found to be not satisfactory to Darcy seepage. The numerical simulation technology of the conventional heavy oil reservoir mostly does not consider the change process of the shear thinning of the heavy oil and the nonlinear seepage phenomenon. These phenomena have a great influence on the development of heavy oil reservoirs.
How to describe the relationship in the numerical simulation process becomes an important objective of increasing the numerical simulation precision. At present, no corresponding heavy oil reservoir numerical simulation method based on the phenomenon is proposed by anyone at home and abroad
Disclosure of Invention
The invention discloses a numerical simulation method considering the non-Newtonian characteristic of thickened oil, which can quantitatively represent the physical phenomenon that the viscosity of the thickened oil reservoir changes along with displacement and the nonlinear seepage law method thereof and has high numerical simulation precision and stability.
In order to achieve the purpose, the invention adopts the following technical scheme:
a numerical simulation method considering non-Newtonian characteristics of thickened oil specifically comprises the following steps:
s1, inputting grid geological parameters;
s2, inputting rock and fluid parameters;
s3, establishing a numerical model and realizing numerical simulation software;
and S4, outputting and displaying the calculation result.
Further, the basic parameters include grid geological parameters and rock and fluid parameters, wherein the grid geological parameters in step S1 include grid number, permeability, porosity, and the like, and the rock and fluid parameters in step S2 include PVT parameters, permeability curves, and the like.
Further, in step S3, the following two ways are selected to characterize the flow law of the thick oil:
(1) characterizing the flow rule of the thickened oil based on the phenomenon that the viscosity of the thickened oil is shear-thinned;
(2) and characterizing the thick oil flowing rule based on the nonlinear seepage phenomenon of the thick oil flowing in the displacement experiment.
Further, in step S3, when characterizing the viscous oil flow law based on the phenomenon of viscous oil shear thinning, the specific steps are as follows:
(1.1) processing the experimental data of the non-Newtonian fluid of the thickened oil, and establishing a fitting relation or a data table between the viscosity of the thickened oil and the shear rate or the flow speed;
(1.2) calculating the flow speed of each grid, further calculating the shear rate of each grid, and correcting the viscosity of the thickened oil of each grid according to the shear rate or the flow speed of each grid;
and (1.3) after obtaining the newly-calculated viscosity value, performing matrix construction and iterative calculation.
Further, in the step (1.1), a shear rate or flow speed is used for representing the phenomenon of viscosity change of the thick oil, and a table or a model of the relationship between the shear rate or flow speed and the viscosity of the thick oil is established;
the shear rate is as follows:
Figure BDA0003510810530000021
wherein gamma is the shear rate, 1/s; v. ofoIs the crude oil flow velocity, is the Darcy velocity, true velocity or phase velocity m/s; k isMesh permeability, m2；kroRelative permeability of oil phase without dimension;
Figure BDA0003510810530000022
is the mesh porosity, dimensionless.
Further, in step S3, when characterizing the thick oil flow law based on the nonlinear seepage phenomenon that occurs when the thick oil flows in the displacement experiment, the specific steps are as follows:
(2.1) fitting the thickened oil displacement experiment data, and establishing a nonlinear seepage rule between the flow of the thickened oil and the displacement pressure gradient, wherein the rule is represented by the experiment data or a fitting curve;
and (2.2) reconstructing the matrix parameters aiming at the nonlinear seepage rule, and finishing iterative computation.
The invention has the beneficial effects that the invention provides two methods for characterizing the flow rule of the thickened oil:
characterizing the flow rule of the thickened oil based on the phenomenon that the viscosity of the thickened oil is shear-thinned;
secondly, representing the flowing rule of the thickened oil based on the nonlinear seepage phenomenon of the flowing thickened oil in the displacement experiment.
The method obviously improves the accuracy of history fitting in the numerical simulation process, reasonably describes the phenomenon of non-Newtonian fluid characteristics and non-linear seepage in the heavy oil reservoir, guides the next step of regulation and control of the heavy oil reservoir for developing a development scheme, and provides theoretical and technical support for developing the development scheme.
Drawings
FIG. 1 is a data chart of a thick oil nonlinear seepage experiment;
FIG. 2 is a schematic flow chart of the present invention;
FIG. 3 is a schematic diagram of a fitting relationship between the viscosity and the shear rate of a heavy oil reservoir in the present invention;
FIG. 4 is a graph comparing water cut for heavy oil reservoirs with and without non-Newtonian behavior;
fig. 5 is a graph comparing the cumulative oil production of heavy oil reservoirs with and without non-newtonian behavior.
Detailed Description
The technical solutions in the embodiments of the present invention will be described clearly and completely with reference to the accompanying 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.
A numerical simulation method considering non-Newtonian characteristics of thickened oil is combined with a thickened oil nonlinear seepage experiment data chart shown in figure 1, and specifically comprises the following steps:
step one, inputting grid geological parameters.
Firstly, inputting grid initialization parameters of a numerical model, including permeability, porosity, net-to-gross ratio, layer top depth and the like, as shown in table 1:
TABLE 1
Parameter(s) Numerical value Parameter(s) Numerical value
Permeability rate of penetration 800mD Number of cells 31×31×5
Porosity of 0.2 Size of the grid 10m×10m×5m
Viscosity of oil phase 200mPa·s Viscosity of aqueous phase 0.49mPa·s
And step two, inputting rock and fluid parameters and inputting parameters of different phase fluids in the grid, wherein the parameters comprise PVT parameters, phase seepage curves and the like. As shown in table 2 below.
TABLE 2
Figure BDA0003510810530000031
Figure BDA0003510810530000041
And step three, establishing a numerical model and realizing numerical simulation software.
The flow law of the thickened oil is characterized by selecting the following two modes:
one way is to characterize the flow law of the thickened oil based on the phenomenon that the viscosity of the thickened oil is shear-thinned;
the method comprises the following specific steps:
processing the experimental data of the non-Newtonian fluid of the thickened oil, and establishing a fitting relation or a data table of the viscosity of the thickened oil and the shear rate or the flow speed; the phenomenon of viscosity change of the thick oil is represented by utilizing the shear rate or the flow speed, and a table or a model of the relationship between the shear rate or the flow speed and the viscosity of the thick oil is established; the flow velocity refers to the migration velocity of reservoir fluid in the reservoir, and the shear rate refers to the change velocity of the fluid flow velocity relative to the pore, which can be calculated from the parameters of the fluid flow velocity, the porosity and the like.
The shear rate is as follows:
Figure BDA0003510810530000042
wherein gamma is the shear rate, 1/s; v. ofoIs the fluid flow velocity, is the Darcy velocity, true velocity or phase velocity m/s; k is the mesh permeability, m2；kroRelative permeability of oil phase without dimension;
Figure BDA0003510810530000043
is the mesh porosity, dimensionless.
(1.2) calculating the flow speed of each grid, further calculating the shear rate of each grid, and correcting the viscosity of the thickened oil of each grid according to the shear rate or the flow speed of each grid;
and (1.3) after obtaining the newly-calculated viscosity value, performing matrix construction and iterative calculation.
The other mode is to characterize the thick oil flowing rule based on the nonlinear seepage phenomenon of the thick oil flowing in the displacement experiment.
The method comprises the following specific steps:
(2.1) fitting the thickened oil displacement experiment data, and establishing a nonlinear seepage rule between the flow of the thickened oil and the displacement pressure gradient, wherein the rule is represented by the experiment data or a fitting curve;
and (2.2) reconstructing the matrix parameters aiming at the nonlinear seepage rule, and finishing iterative computation.
And recalculating the oil phase flow velocity or shear rate of each grid in each iteration process, and calculating the viscosity of the heavy oil reservoir on the basis of the flow velocity or shear rate, thereby providing an effective numerical simulation method for the heavy oil reservoir.
According to the viscous oil flowing experimental data, a nonlinear seepage model is established to represent the seepage rule of the viscous oil flowing experimental data, the nonlinear seepage model can be obtained by the experimental data or a fitting theoretical model, and a theoretical calculation curve can represent various forms such as power-law fluid, Bingham fluid and the like.
As can be seen from fig. 4 and 5, after the model considers the change of the viscosity of the thick oil along with the displacement, the water cut rise is reduced, and the accumulated oil yield is higher than that of the model which does not consider the change of the viscosity of the thick oil. After the viscosity change of the heavy oil is considered, the viscosity of the heavy oil on the main flow channel is reduced along with the development, the reduction range near a production well is larger, the oil-water fluidity ratio is reduced, the water phase inrush is slowed down, and the water content rise is reduced after the water breakthrough time of the oil reservoir is late. When the actual block history fitting is started, a viscous oil viscosity change model can be utilized to improve the accuracy of the history fitting.
The method can quantitatively represent the physical phenomenon that the viscosity of the heavy oil reservoir changes along with displacement and the nonlinear seepage law method thereof, and has high numerical simulation precision and stability.
It is to be understood that the above description is not intended to limit the present invention, and the present invention is not limited to the above examples, and those skilled in the art may make various changes, modifications, additions and substitutions within the spirit and scope of the present invention.

Claims (6)

1. A numerical simulation method considering non-Newtonian characteristics of thickened oil is characterized by comprising the following steps:
s1, inputting grid geological parameters;
s2, inputting rock and fluid parameters;
s3, establishing a numerical model and realizing numerical simulation software;
and S4, outputting and displaying the calculation result.
2. The numerical simulation method considering non-Newtonian characteristics of heavy oil of claim 1, wherein the basic parameters include grid geological parameters and rock and fluid parameters, wherein the grid geological parameters in step S1 include grid number, permeability and porosity, and the rock and fluid parameters in step S2 include PVT parameters and phase permeability curves.
3. A numerical simulation method considering non-newtonian characteristics of thick oil according to claim 2, wherein in step S3, the flow law of thick oil is characterized by selecting two ways:
(1) characterizing the flow rule of the thickened oil based on the phenomenon that the viscosity of the thickened oil is shear-thinned;
(2) and characterizing the thick oil flowing rule based on the nonlinear seepage phenomenon of the thick oil flowing in the displacement experiment.
4. A numerical simulation method considering non-newtonian characteristics of thick oil as claimed in claim 3, wherein in step S3, when characterizing the flow law of thick oil based on the phenomenon of viscosity shear thinning of thick oil, the specific steps are as follows:
(1.1) processing the experimental data of the non-Newtonian fluid of the thickened oil, and establishing a fitting relation or a data table of the viscosity of the thickened oil and the shear rate or the flow speed;
(1.2) calculating the flow speed of each grid, further calculating the shear rate of each grid, and correcting the viscosity of the thickened oil of each grid according to the shear rate or the flow speed of each grid;
and (1.3) after obtaining the newly-calculated viscosity value, performing matrix construction and iterative calculation.
5. A numerical simulation method considering non-newtonian characteristics of thick oil according to claim 1, wherein in step (1.1), a shear rate or a flow rate is used to characterize the viscosity change of thick oil, and a table or a model of the relationship between the shear rate or the flow rate and the viscosity of thick oil is established;
the shear rate is as follows:
Figure FDA0003510810520000011
wherein gamma is the shear rate, 1/s; v. ofoIs the crude oil flow velocity, is the Darcy velocity, true velocity or phase velocity m/s; k is the mesh permeability, m2；kroRelative permeability of oil phase without dimension;
Figure FDA0003510810520000012
is the mesh porosity, dimensionless.
6. A numerical simulation method considering non-newtonian characteristics of thick oil according to claim 1, wherein in step S3, when characterizing the thick oil flow law based on the nonlinear seepage phenomenon occurred in the thick oil flow in the displacement experiment, the specific steps are as follows:
(2.1) fitting the thickened oil displacement experiment data, and establishing a nonlinear seepage rule between the flow of the thickened oil and the displacement pressure gradient, wherein the rule is represented by the experiment data or a fitting curve;
and (2.2) reconstructing the matrix parameters aiming at the nonlinear seepage rule, and finishing iterative computation.
CN202210152007.XA 2022-02-18 2022-02-18 Numerical simulation method considering non-Newtonian characteristics of thickened oil Pending CN114528780A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210152007.XA CN114528780A (en) 2022-02-18 2022-02-18 Numerical simulation method considering non-Newtonian characteristics of thickened oil

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210152007.XA CN114528780A (en) 2022-02-18 2022-02-18 Numerical simulation method considering non-Newtonian characteristics of thickened oil

Publications (1)

Publication Number Publication Date
CN114528780A true CN114528780A (en) 2022-05-24

Family

ID=81622754

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210152007.XA Pending CN114528780A (en) 2022-02-18 2022-02-18 Numerical simulation method considering non-Newtonian characteristics of thickened oil

Country Status (1)

Country Link
CN (1) CN114528780A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115479867A (en) * 2022-08-02 2022-12-16 广东工业大学 non-Newtonian fluid rapid discrimination method
CN116822391A (en) * 2022-12-15 2023-09-29 长江大学 Heavy oil reservoir bulk phase fluid nonlinear seepage theory model

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115479867A (en) * 2022-08-02 2022-12-16 广东工业大学 non-Newtonian fluid rapid discrimination method
CN115479867B (en) * 2022-08-02 2024-05-28 广东工业大学 Method for rapidly distinguishing non-Newtonian fluid
CN116822391A (en) * 2022-12-15 2023-09-29 长江大学 Heavy oil reservoir bulk phase fluid nonlinear seepage theory model
CN116822391B (en) * 2022-12-15 2024-05-07 长江大学 Heavy oil reservoir bulk phase fluid nonlinear seepage theory method

Similar Documents

Publication Publication Date Title
CN114528780A (en) Numerical simulation method considering non-Newtonian characteristics of thickened oil
CN103939066B (en) An a kind of note adopts the method that well group determines water injection rate determination oilwell produced fluid amount more
CN105386751A (en) Well logging and productivity prediction method of horizontal well based on oil reservoir seepage flow model
CN107327299B (en) Method and device for determining reservoir compressibility
CN111810101B (en) Dynamic analysis method and device for water-drive reservoir
CN109858177B (en) Bottom water reservoir horizontal well water drive numerical simulation method based on pseudo-streamline method
CN111852417B (en) Rapid inversion method and device for gas channeling channel of carbon dioxide flooding reservoir
CN105426620A (en) Quantitative analysis method and device of main controlling factor of oil layer remaining oil
CN104533519B (en) Management method for water burst water disasters in event of vertical shaft wellbore passing through strong water-containing thick rock layer
CN108536982B (en) Evaluation method for multi-crack oil-water mixing cross displacement
CN107355200B (en) Method for improving water drive well selection by nano-micron particle dispersion system
CN116467958A (en) Salt lake brine numerical model construction and water supplementing and mineral dissolving efficiency calculation method
CN112528573B (en) Plugging agent parameter optimization method and system for oil well in ultrahigh water cut period
Rushton Significance of a seepage face on flows to wells in unconfined aquifers
CN112796749A (en) Oil reservoir water drive numerical simulation method based on number of capillary tubes
CN114737947A (en) Method and device for developing and adjusting horizontal well of carbonate reservoir
CN113326465A (en) Oil reservoir development dynamic analysis method and device
CN114004078B (en) Method for calculating multilayer commingled production seepage simulation phase seepage of thin interbed oil reservoir
CN111322054A (en) Three-three combined excavation and potential optimization mining method for sandstone oil reservoir in chemical flooding stage
CN112651123B (en) Ternary combination flooding injection scheme design method
CN112901126B (en) Yield splitting method based on phase permeation and water flooding characteristics
Kootiani et al. Analysis Fraction Flow of Water versus Cumulative Oil Recoveries Using Buckley Leverett Method
Fan et al. Evaluation of sweep efficiency in flooding process of reservoir development using substitution index
CN112196527B (en) Method for determining water body size of fracture-cavity type oil reservoir
Almukhametova The efficiency of oil reserve recovery in a reservoir of uniform permeability when changing the operating mode of development well

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