CN110080730A - Fracture hole type series connection reservoir individual well pressure prediction model modeling method and storage medium - Google Patents

Fracture hole type series connection reservoir individual well pressure prediction model modeling method and storage medium Download PDF

Info

Publication number
CN110080730A
CN110080730A CN201910342351.3A CN201910342351A CN110080730A CN 110080730 A CN110080730 A CN 110080730A CN 201910342351 A CN201910342351 A CN 201910342351A CN 110080730 A CN110080730 A CN 110080730A
Authority
CN
China
Prior art keywords
solution cavity
phase
pressure
pressure drop
crack
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
CN201910342351.3A
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.)
Southwest Petroleum University
Original Assignee
Southwest Petroleum University
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 Southwest Petroleum University filed Critical Southwest Petroleum University
Priority to CN201910342351.3A priority Critical patent/CN110080730A/en
Publication of CN110080730A publication Critical patent/CN110080730A/en
Priority to PCT/CN2020/086542 priority patent/WO2020216313A1/en
Pending legal-status Critical Current

Links

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH DRILLING; MINING
    • E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
    • E21B43/166—Injecting a gaseous medium; Injecting a gaseous medium and a liquid medium
    • E21B43/168—Injecting a gaseous medium

Abstract

The invention discloses fracture hole type series connection reservoir individual well pressure prediction model modeling method and storage mediums, solution cavity approximation is taken as cylinder, it is connected between solution cavity by fisstured flow system, flowing of the fluid in fracture-cavity type carbonate reservoir includes vug unit and Crack Element, these two types of units alternatelies occur constituting entire fracture hole coupled system, so that capillary bunch physical model is abstracted as solution cavity, the mathematical model that crack is sequentially connected in series.The pressure drop fundamental equation of vug unit and Crack Element is integrated further according to the boundary condition of solution cavity and crack junction, obtains overall presure drop model by the pressure drop fundamental equation for establishing vug unit, Crack Element respectively.The present invention solves the problems, such as that fractured-cavernous carbonate reservoir modeling is difficult in the prior art, is difficult to predict the pressure drop situation around pit shaft, realizes modeling method, the purpose convenient for being predicted individual well pressure for providing architectonical for fractured-cavernous carbonate reservoir.

Description

Fracture hole type series connection reservoir individual well pressure prediction model modeling method and storage medium
Technical field
The present invention relates to fractured-cavernous carbonate reservoir development fields, and in particular to fracture hole type series connection reservoir individual well pressure is pre- Survey model modelling approach and storage medium.
Background technique
Fracture and vug carbonate reservoir has and buries the spies such as depth, temperature are high, reservoir space is complicated, oil and water zonation is irregular Sign, in China using system in Tahe Oilfield as Typical Representative.The Development Practice of fractured-cavernous carbonate reservoir shows to deposit in fractured-vuggy reservoir In the Coupled Flow of the flowing of solution cavity large scale and crack Porous Media composition, this is the physics of fractured-cavernous carbonate reservoir Model and mathematics model foundation bring great difficulty, cause to be difficult to predict the pressure drop situation around pit shaft, can not Accurate technological parameter foundation is provided for gas injection volume increase by pressure prediction as conventional reservoir.
Summary of the invention
The purpose of the present invention is to provide fracture hole type series connection reservoir individual well pressure prediction model modeling method and storage medium, To solve in the prior art, fractured-cavernous carbonate reservoir modeling is difficult, is difficult to predict the pressure drop situation around pit shaft Problem is embodied as fractured-cavernous carbonate reservoir and provides the modeling method of architectonical, the mesh convenient for being predicted individual well pressure 's.
The present invention is achieved through the following technical solutions:
Fracture hole type series connection reservoir individual well pressure prediction model modeling method, including successively replacing connected solution cavity, crack, packet Include following steps:
(a) solution cavity is established as to cylindrical vug unit, crack is established as to the Crack Element of capillary；Solution cavity list Fluid flowing in member is equivalent to irregular pipe stream, and the fluid flowing in Crack Element is equivalent to the seepage flow of porous media；Its The axis of middle cylindrical model and the axis of capillary model are parallel to each other；
(b) the pressure drop fundamental equation for establishing vug unit establishes the pressure drop fundamental equation of Crack Element；
(c) boundary condition for determining solution cavity and crack junction, makes PVa=Pfb,QV=Qf, in which:
PVa、Pfb--- the pressure of each comfortable tie point of vug unit, Crack Element, MPa；
QV、Qf--- the flow of each comfortable tie point of vug unit, Crack Element, m3/s；
(d) algebraical sum that all solution cavities are generated to pressure drop is added with the algebraical sum that all slits generate pressure drop, obtains stagnation pressure Model drops.
The Development Practice of fractured-cavernous carbonate reservoir shows: there are the flowing of solution cavity large scale and cracks in fractured-vuggy reservoir The Coupled Flow of Porous Media composition.For this feature of reservoir, the invention proposes a kind of pipe stream~seepage flow couplings Solution cavity approximation is taken as cylinder, is connected between solution cavity by fisstured flow system by flow model, and a fracture-cavity units are considered At a kind of netted physical model --- i.e. bundle model.Fluid flowing in solution cavity is approximately considered irregular pipe stream, i.e., Flowing of the fluid in round tube；Fluid flowing in crack is approximately considered the seepage flow of porous media.On this basis, ignore weight The influence of power carrys out founding mathematical models.Flowing of the fluid in fracture-cavity type carbonate reservoir includes two class units: solution cavity list altogether Member and Crack Element, these two types of units alternatelies occur constituting entire fracture hole coupled system, so that capillary bunch physical model be taken out As the mathematical model being sequentially connected in series for solution cavity, crack.Herein on basis, the pressure drop of vug unit, Crack Element is established respectively Fundamental equation, further according to the boundary condition of solution cavity and crack junction, by the pressure drop fundamental equation of vug unit and Crack Element It integrates, is convenient for direct algebraic addition, so that overall presure drop model is obtained, it can be all to pit shaft by obtained overall presure drop model The pressure drop situation enclosed is predicted, and then provides accurate technological parameter foundation for gas injection volume increase.The present invention solves existing skill Fractured-cavernous carbonate reservoir modeling is difficult in art, is difficult to the problem of predicting the pressure drop situation around pit shaft, realizes Modeling method, the purpose convenient for being predicted individual well pressure for providing architectonical for fractured-cavernous carbonate reservoir.
Consider the fluid type in reservoir, it is dynamic that the present invention is divided into oily phase one-way flow, oil-water two-phase flow, oil-gas two-phase flow Three kinds of situations each provide specific mathematical model and establish formula:
When in solution cavity, crack being oily phase one-way flow:
The pressure drop fundamental equation of the vug unit are as follows: Wherein, αa'=αa+ζa, αb'=αb-ζb；Export correction factorEntrance correction factor
In formula: PVa、PVb--- vug unit outlet, the respective pressure of entrance, MPa；LV--- the length of solution cavity, m； The equivalent diameter of d --- solution cavity, m；λ --- loss factor；AVa、AVb--- vug unit entrance, outlet are respective transversal Area, m2；wfαa、αb--- vug unit outlet, the respective kinetic energy correction factor of entrance；ρ --- fluid density；AV--- solution cavity Sectional area, m2；
The pressure drop fundamental equation of the Crack Element are as follows:In formula: Pfa、Pfb--- it splits Stitch unit outlet, the respective pressure of entrance, MPa；Lf--- fracture length, m；Wf--- fracture width, mm；μ --- fluid is viscous Degree, mPas；hf--- fracture height, m.
When in solution cavity, crack being oil-water two-phase flow:
Pressure drop of the water-oil phase in vug unit along journey pressure drop, partial drop of pressure and acceleration pressure by coming down to a lower group into, the solution cavity The pressure drop fundamental equation of unit are as follows: Δ PV=Δ PVw+ΔPVp+ΔPVs；The pressure drop fundamental equation of the Crack Element are as follows:In formula: Δ PV--- vug unit overall presure drop, Pa；ΔPf--- vug unit overall presure drop, Pa； ΔPVw1--- water-oil phase solution cavity is along journey pressure drop, Pa；ΔPVp1--- water-oil phase solution cavity partial drop of pressure, Pa；ΔPVs1--- oil Water two-phase solution cavity accelerates pressure drop, Pa；Qfm--- water-oil phase fluid flow, m in crack3/ks；Lf--- fracture length, m； Wf--- fracture width, mm；H --- fracture height, m；Kro、Krw--- the respective relative permeability of oily phase, water phase, μm2；μo、 μw--- the respective viscosity of oily phase, water phase, mPas.
Wherein, the water-oil phase solution cavity is along journey pressure dropThe water-oil phase solution cavity part Pressure dropThe water-oil phase solution cavity accelerates pressure drop In formula: λm--- oil water mixture loss factor；ρm--- oil water mixture density, kg/m3；VVm--- grease in solution cavity The average speed of mixture, m/s；LV--- the length of solution cavity, m；The equivalent diameter of d --- solution cavity, m；λ --- linear loss system Number；VVam--- solution cavity exports the flow velocity of mixture, m/s；VVbm--- the flow velocity of solution cavity inlet mixture, m/s；αam、αbm—— Solution cavity outlet, the respective kinetic energy correction factor of entrance；Export correction factorEntrance correction factorAVa、AVb--- vug unit entrance, the respective cross-sectional area in outlet, m2。
When being moved in solution cavity, crack for oil-gas two-phase flow:
Pressure drop of the oil gas two-phase in vug unit along journey pressure drop, partial drop of pressure and acceleration pressure by coming down to a lower group into, the solution cavity The pressure drop fundamental equation of unit are as follows: Δ PV=Δ PVw+ΔPVp+ΔPVs；The pressure drop fundamental equation of the Crack Element are as follows:In formula: Δ PV--- vug unit overall presure drop, Pa；ΔPf--- vug unit Overall presure drop, Pa；ΔPVw2--- oil gas two-phase solution cavity is along journey pressure drop, Pa；ΔPVp2--- oil gas two-phase solution cavity partial drop of pressure, Pa； ΔPVs2--- oil gas two-phase solution cavity accelerates pressure drop, Pa；Gfm--- oil-gas two-phase flow weight flow, kg/s in crack；μo、 --- the respective viscosity of oily phase, gas phase, mPas；ρo、 --- the respective relative permeability of oily phase, gas phase, μm2；Lf--- fracture length, m；Wf--- fracture width, mm；hf--- it splits Seam height, m.
Wherein, the oil gas two-phase solution cavity is along journey pressure dropThe oil gas two-phase solution cavity is locally pressed DropThe oil gas two-phase solution cavity accelerates pressure drop Its middle outlet correction factorEntrance correction factorIn formula: λm--- gas mixture Loss factor, the i.e. coefficient of friction resistance；ρm--- gas mixture density, kg/m3；GVm--- the matter of gas mixture in solution cavity Measure flow, kg/s；AV--- solution cavity sectional area, m2；LV--- solution cavity length, m；D --- solution cavity radius, m；AVa、AVb--- it is molten Hole unit entrance, the respective cross-sectional area in outlet, m2；VVam--- solution cavity exports the flow velocity of mixture, m/s；VVbm--- solution cavity The flow velocity of inlet mixture, m/s.
Further, a kind of computer readable storage medium, the computer-readable recording medium storage have computer journey Sequence, the computer program realize the step in this method when being executed by processor.The big meter of specifically used dependence of this method It calculates, it is therefore preferable that above-mentioned calculating process is realized by computer program, so any step comprising being protected in this method Rapid computer program and its storage medium also belongs in the protection scope of the application.
Compared with prior art, the present invention having the following advantages and benefits:
Fracture hole type series connection reservoir individual well pressure prediction model modeling method of the present invention and storage medium, solution cavity approximation is taken as Cylinder is connected by fisstured flow system between solution cavity, and a fracture-cavity units consideration is reticulated physical model --- i.e. hair Restrain model.Fluid flowing in solution cavity is approximately considered irregular pipe stream, i.e. flowing of the fluid in round tube；In crack Fluid flowing be approximately considered the seepage flow of porous media.On this basis, the influence for ignoring gravity carrys out founding mathematical models: fluid Flowing in fracture-cavity type carbonate reservoir includes two class units: vug unit and Crack Element, these two types of units alternatelies altogether Appearance constitutes entire fracture hole coupled system, so that capillary bunch physical model is abstracted as solution cavity, the mathematics that crack is sequentially connected in series Model.Herein on basis, the pressure drop fundamental equation of vug unit, Crack Element is established respectively, is connected further according to solution cavity and crack The boundary condition for meeting place integrates the pressure drop fundamental equation of vug unit and Crack Element, is convenient for direct algebraic addition, from And overall presure drop model is obtained, the pressure drop situation around pit shaft can be predicted by obtained overall presure drop model, Jin Erwei Gas injection volume increase provides accurate technological parameter foundation.It is tired that the present invention solves fractured-cavernous carbonate reservoir modeling in the prior art Hardly possible is difficult to the problem of predicting the pressure drop situation around pit shaft, realizes and provides adult for fractured-cavernous carbonate reservoir The modeling method of system, the purpose convenient for being predicted individual well pressure.
Detailed description of the invention
Attached drawing described herein is used to provide to further understand the embodiment of the present invention, constitutes one of the application Point, do not constitute the restriction to the embodiment of the present invention.In the accompanying drawings:
Fig. 1 is the physical model schematic diagram of the specific embodiment of the invention；
Fig. 2 is the mathematical model schematic diagram of the specific embodiment of the invention.
Specific embodiment
To make the objectives, technical solutions, and advantages of the present invention clearer, below with reference to embodiment and attached drawing, to this Invention is described in further detail, and exemplary embodiment of the invention and its explanation for explaining only the invention, are not made For limitation of the invention.
Embodiment 1:
When individual well is oily phase single-phase flow:
Solution cavity is established as to cylindrical vug unit, crack is established as to the Crack Element of capillary；Vug unit In fluid flowing be equivalent to irregular pipe stream, the flowing of fluid in Crack Element is equivalent to the seepage flow of porous media；Wherein The axis of cylindrical model and the axis of capillary model are parallel to each other；Form physical model as shown in Figure 1.Ignore gravity It influences.Flowing is divided into two class units: vug unit and Crack Element, these two types of units alternatelies occur constituting entire fracture hole coupling Collaboration system.It therefore is mathematical model as shown in Figure 2 by said physical Model Abstraction.Entire fracture-cavity units just by vug unit with Crack Element is interconnected to constitute.When these two types of different units connect, the flowing of whole system keeps continuity, i.e., in adjacent list The pressure and flow of the tie point upper fluid of member are continuous.
(1) vug unit:
Since the present embodiment is oily phase single-phase flow, the energy loss of oil phase flowing is only considered, in the energy of vug unit Loss is mainly along journey energy loss, then energy loss can indicate are as follows:
In formula: LV--- the length of solution cavity, m；The equivalent diameter of d --- solution cavity, m；
VV--- the average speed of fluid, m/s in solution cavity；λ --- loss factor
When fluid is by solution cavity entrance and exit, due to the sudden enlargement of runner size and reduces and generate energy damage It loses, the loss of this portion of energy is called local energy loss, can indicate are as follows:Wherein entrance amendment system NumberExport correction factor
Total energy loss is to be write as along the sum of journey energy loss and local energy loss in solution cavity:
Since the fluid in unit is viscous fluid, the VELOCITY DISTRIBUTION of each cross section upper liquid constitution point of unit is uneven 's.Therefore, the kinetic energy integral term on each section is calculated with kinetic energy correction factor multiplied by the mean flow rate in section, then solution cavity goes out The flow parameter of mouth and entrance are as follows:In formula: Ha、Hb--- solution cavity The head of unit outlet and entrance, m；PVa、PVb--- the pressure of vug unit outlet and entrance, MPa；VVa、VVb--- solution cavity list The section average speed of member outlet and entrance, m/s；α1、α2--- the kinetic energy correction factor of vug unit outlet and entrance.
It can be derived according to bernoulli equation:In formula: αa'=αa+ζa, αb'=αb-ζb。
Had according to continuity equation: QV=VVAV=VVaAVa=VVbAVb.In formula: A1、A2--- vug unit entrance and go out The cross-sectional area of mouth, m2。
To sum up, that is, the pressure drop fundamental equation of vug unit when may make up oily phase single-phase flow:
(2) Crack Element:
Next the pressure drop fundamental equation of analysis Crack Element: the flowing in crack can see linear seepage flow, according to splitting The pungent Lek equation seepage theory of cloth in seam, flow equation may be expressed as:In formula: Pfa、 Pfb--- the pressure of Crack Element outlet and entrance, MPa；μ --- fluid viscosity, mPas；Lf--- fracture length, m； wf--- fracture width, mm；Vf--- percolation flow velocity, m/s.To each Crack Element, seepage flow is considered as one-dimensional incompressible Fluid neuron network, then continuity equation are as follows:That is Vf=const, Qf=const.In formula: L --- the seepage flow of Crack Element Direction length, m.
By above-mentioned flow equation and continuity equation, the pressure drop of Crack Element is substantially square when constituting oily phase single-phase flow Journey:
In formula: Qf--- fluid flow in crack, m3/s；hf--- fracture height, m。
The pressure drop fundamental equation of vug unit when oil phase single-phase flow, Crack Element obtained in the present embodiment is subjected to group It closes, the foundation of complete individual well pressure drop mathematical model can be completed.
It is further:
Consider that fracture hole model has n solution cavity, each solution cavity contributes producing well also different, vacation due to Reservoir Fluid volume difference If i-th of solution cavity is with flow Qi production, then solution cavity internal drop can be expressed as the algebraical sum that n solution cavity generates pressure drop:
Fracture hole model has n crack simultaneously, since crack is not as reservoir space, to producing well without traffic contributions, therefore The flow of i-th crack and i-th of solution cavity are equal with flow, i.e. Qfi=Qi, then crack internal drop can be expressed as n crack Generate the algebraical sum of pressure drop:
Therefore, the present embodiment can also expand the fracture hole model stitched for the hole n n, and stagnation pressure is reduced to solution cavity pressure drop and crack pressure drop The sum of, it can indicate are as follows:
Embodiment 2:
When individual well is oil-water two-phase flow:
Solution cavity is established as to cylindrical vug unit, crack is established as to the Crack Element of capillary；Vug unit In fluid flowing be equivalent to irregular pipe stream, the flowing of fluid in Crack Element is equivalent to the seepage flow of porous media；Wherein The axis of cylindrical model and the axis of capillary model are parallel to each other；Form physical model as shown in Figure 1.Ignore gravity It influences.Flowing is divided into two class units: vug unit and Crack Element, these two types of units alternatelies occur constituting entire fracture hole coupling Collaboration system.It therefore is mathematical model as shown in Figure 2 by said physical Model Abstraction.Entire fracture-cavity units just by vug unit with Crack Element is interconnected to constitute.When these two types of different units connect, the flowing of whole system keeps continuity, i.e., in adjacent list The pressure and flow of the tie point upper fluid of member are continuous.
(1) vug unit:
According to the flow condition of oil-water two-phase flow, it is three parts that the overall presure drop in solution cavity, which is drawn total score, they respectively by Is defined as: along journey pressure drop, partial drop of pressure and accelerate pressure drop.
1. along journey pressure drop in water-oil phase vug unit
It is mainly along journey energy loss in the energy loss of vug unit, then energy loss indicates are as follows:
In formula: λm--- oil water mixture loss factor, decimal；ρm--- grease is mixed Close object density, kg/m3；VVm--- the average speed of oil water mixture, m/s in solution cavity.
2. partial drop of pressure in water-oil phase vug unit
When water-oil phase is by solution cavity entrance and exit, due to the sudden enlargement of runner size and reduces and generate energy Loss, the loss of this portion of energy are local energy loss, are indicated are as follows:Wherein, outlet is repaired Positive coefficientEntrance correction factorIn formula: VVam--- the stream of solution cavity outlet mixture Speed, m/s；VVbm--- the flow velocity of solution cavity inlet mixture, m/s.
3. accelerating pressure drop in water-oil phase vug unit
Accelerate pressure drop Δ PVsBe due to caused by the difference of solution cavity both ends flow velocity, it only and oil water mixture density, flow velocity It is related.According to the speed difference at solution cavity both ends, the present embodiment will accelerate voltage drop meter to be shown as:
Wherein, αam、αbmFor outlet and entrance kinetic energy correction factor, usually Take 1.
4. water-oil phase vug unit overall presure drop
Water-oil phase is equal in the overall presure drop that solution cavity flows along the sum of journey pressure drop, partial drop of pressure and acceleration pressure drop three, molten Hole unit Pressure Drop Model may be expressed as:In formula: αam'= αam+ζa、αbm'=αbm-ζb。
Had according to continuity equation: QVm=VVmAV=VVamAVa=VVbmAVb.In formula: Qfm--- oil-water two-phase flow in solution cavity Body flow, m3/ks；
Aforementioned Pressure Drop Model and continuity equation are combined, considers pressure unit Pa being converted to MPa, grease two can be obtained The pressure drop fundamental equation of vug unit when mutually flowing:
(2) Crack Element:
The present embodiment when the mathematical model of Crack Element, assumes initially that oil, water two-phase exist when establishing oil-water two-phase flow Obedience darcy flow is flowed in crack, in order to provide the modeling systems to match with embodiment 1.It is specific:
Initially set up oily phase flow model:Establish water phase flow model:
Therefore, it may be expressed as: by total liquid measure of any cross-section of river
In formula: Qfm--- water-oil phase fluid flow, m in crack3/ks；Qfm--- oily phase flow rate, m in crack3/ks；
Qfm--- crack inner aqueous phase flow, m3/ks；μo、μw--- oily phase, the viscosity of water phase, mPas；
Kro、Krw--- oily phase, the relative permeability of water phase, μm2；
Total liquid measure formula is integrated to get pressure drop fundamental equation when flowing in Crack Element to water-oil phase::
The pressure drop fundamental equation of vug unit when oil-water two-phase flow, Crack Element obtained in the present embodiment is subjected to group It closes, the foundation of complete individual well pressure drop mathematical model can be completed.
Embodiment 3:
When individual well moves for oil-gas two-phase flow:
Solution cavity is established as to cylindrical vug unit, crack is established as to the Crack Element of capillary；Vug unit In fluid flowing be equivalent to irregular pipe stream, the flowing of fluid in Crack Element is equivalent to the seepage flow of porous media；Wherein The axis of cylindrical model and the axis of capillary model are parallel to each other；Form physical model as shown in Figure 1.Ignore gravity It influences.Flowing is divided into two class units: vug unit and Crack Element, these two types of units alternatelies occur constituting entire fracture hole coupling Collaboration system.It therefore is mathematical model as shown in Figure 2 by said physical Model Abstraction.Entire fracture-cavity units just by vug unit with Crack Element is interconnected to constitute.When these two types of different units connect, the flowing of whole system keeps continuity, i.e., in adjacent list The pressure and flow of the tie point upper fluid of member are continuous.
(1) vug unit:
According to the flow condition of oil-gas two-phase flow, it is three parts that the overall presure drop in solution cavity, which is drawn total score, they respectively by Is defined as: along journey pressure drop, partial drop of pressure and accelerate pressure drop.
1. along journey pressure drop in oil gas two-phase vug unit: it is mainly edge that oil-gas two-phase flow, which moves the energy loss in vug unit, Journey energy loss, then energy loss can indicate are as follows:In formula: λm--- gas mixture is damaged along journey Lose coefficient, the i.e. coefficient of friction resistance；ρm--- gas mixture density, kg/m3；GVm--- the quality stream of gas mixture in solution cavity Amount, kg/s；AV--- solution cavity sectional area, m2；LV--- solution cavity length, m；D --- solution cavity radius, m.
2. partial drop of pressure in oil gas two-phase vug unit: when oil gas two-phase is by solution cavity entrance and exit, due to stream Road size sudden enlargement generates energy loss with reducing, and the loss of this portion of energy is called local energy loss, can indicate are as follows:Its middle outlet correction factorEntrance correction factor
3. accelerating pressure drop in oil gas two-phase vug unit: accelerating pressure drop Δ PVsIt is since the difference of solution cavity both ends flow velocity causes , it is only related with gas mixture density, flow velocity, according to the speed difference at solution cavity both ends, accelerates pressure drop that can indicate are as follows:
4. oil gas two-phase vug unit overall presure drop:
Oil gas two-phase is equal to along the sum of journey pressure drop, partial drop of pressure and acceleration pressure drop three, then in the overall presure drop that solution cavity flows Vug unit Pressure Drop Model may be expressed as:In formula: αam' =αam+ζa, αbm'=αbm-ζb。
Had according to continuity equation: GVm=ρmVVmAV=ρmVVamAVa=ρmVVbmAVb, in formula: Aa、Ab--- vug unit goes out The cross-sectional area of mouth and entrance, m2。
Aforementioned Pressure Drop Model and continuity equation are combined, considers pressure unit Pa being converted to MPa, oil gas two can be obtained The pressure drop fundamental equation of vug unit when mutually flowing:
(2) Crack Element:
It drives a well after oil well gas injection, will occur oil, gas two-phase flow in crack, it is assumed that oil, gas two-phase flow in crack Obey darcy flow.
Oily phase flow model:
Gas phase flow model:In formula: Gfo、Gfg--- it is respectively the quality stream of oil, gas Amount, kg/s.
Then oily, gas two-phase total mass flow rate can be calculated by following formula:
Above formula both sides are integrated:
Define gas and oil two phase pseudo pressure function:
Using gas and oil two phase pseudo pressure function representation, then for
Crack section is divided into n equal part, the ρ in each small Crack Elemento、μo、ρg、μgFor definite value, pass through the pressure of each segment Drop superposition can be in the hope of the pressure drop of entire crack section:
Above formula is carried out to simplify consideration, in pressure P for oily phasefa~PfbInterior ρo、μoTo be thought of as definite value, to gas phase Speech is in pressure Pfa~PfbIt inside takes under average pressureFinal oil gas two-phase then can be obtained to flow in Crack Element When pressure drop fundamental equation:
The pressure drop fundamental equation of vug unit when dynamic of oil-gas two-phase flow obtained in the present embodiment, Crack Element is subjected to group It closes, the foundation of complete individual well pressure drop mathematical model can be completed.
Embodiment 4:
A kind of computer readable storage medium, the computer-readable recording medium storage have computer program, the meter Calculation machine program realizes the step in any of the above-described embodiment when being executed by processor.
Above-described specific embodiment has carried out further the purpose of the present invention, technical scheme and beneficial effects It is described in detail, it should be understood that being not intended to limit the present invention the foregoing is merely a specific embodiment of the invention Protection scope, all within the spirits and principles of the present invention, any modification, equivalent substitution, improvement and etc. done should all include Within protection scope of the present invention.

Claims (7)

  1. The reservoir individual well pressure prediction model modeling method 1. fracture hole type is connected, it is special including successively replacing connected solution cavity, crack Sign is, comprising the following steps:
    (a) solution cavity is established as to cylindrical vug unit, crack is established as to the Crack Element of capillary；In vug unit Fluid flowing be equivalent to irregular pipe stream, the flowing of fluid in Crack Element is equivalent to the seepage flow of porous media；Wherein justify The axis of cylindricality model and the axis of capillary model are parallel to each other；
    (b) the pressure drop fundamental equation for establishing vug unit establishes the pressure drop fundamental equation of Crack Element；
    (c) boundary condition for determining solution cavity and crack junction, makes PVa=Pfb,QV=Qf, in which:
    PVa、Pfb--- the pressure of each comfortable tie point of vug unit, Crack Element, MPa；
    QV、Qf--- the flow of each comfortable tie point of vug unit, Crack Element, m3/s；
    (d) algebraical sum that all solution cavities are generated to pressure drop is added with the algebraical sum that all slits generate pressure drop, obtains overall presure drop mould Type.
  2. The reservoir individual well pressure prediction model modeling method 2. fracture hole type according to claim 1 is connected, which is characterized in that when When being oily phase one-way flow in solution cavity, crack:
    The pressure drop fundamental equation of the vug unit are as follows:Its In, αa'=αa+ζa, αb'=αb-ζb；Export correction factorEntrance correction factor
    In formula:
    PVa、PVb--- vug unit outlet, the respective pressure of entrance, MPa；
    LV--- the length of solution cavity, m；The equivalent diameter of d --- solution cavity, m；λ --- loss factor；
    AVa、AVb--- vug unit entrance, the respective cross-sectional area in outlet, m2；
    wfαa、αb--- vug unit outlet, the respective kinetic energy correction factor of entrance；
    ρ --- fluid density；AV--- solution cavity sectional area, m2；
    The pressure drop fundamental equation of the Crack Element are as follows:In formula:
    Pfa、Pfb--- Crack Element outlet, the respective pressure of entrance, MPa；
    Lf--- fracture length, m；Wf--- fracture width, mm；
    μ --- fluid viscosity, mPas；hf--- fracture height, m.
  3. The reservoir individual well pressure prediction model modeling method 3. fracture hole type according to claim 1 is connected, which is characterized in that when When being oil-water two-phase flow in solution cavity, crack:
    Pressure drop of the water-oil phase in vug unit along journey pressure drop, partial drop of pressure and acceleration pressure by coming down to a lower group into, the vug unit Pressure drop fundamental equation are as follows: Δ PV=Δ PVw+ΔPVp+ΔPVs；
    The pressure drop fundamental equation of the Crack Element are as follows:In formula:
    ΔPV--- vug unit overall presure drop, Pa；ΔPf--- vug unit overall presure drop, Pa；
    ΔPVw1--- water-oil phase solution cavity is along journey pressure drop, Pa；ΔPVp1--- water-oil phase solution cavity partial drop of pressure, Pa；
    ΔPVs1--- water-oil phase solution cavity accelerates pressure drop, Pa；Qfm--- water-oil phase fluid flow, m in crack3/ks；
    Lf--- fracture length, m；Wf--- fracture width, mm；H --- fracture height, m；
    Kro、Krw--- the respective relative permeability of oily phase, water phase, μm2；
    μo、μw--- the respective viscosity of oily phase, water phase, mPas.
  4. The reservoir individual well pressure prediction model modeling method 4. fracture hole type according to claim 3 is connected, which is characterized in that
    The water-oil phase solution cavity is along journey pressure drop
    The water-oil phase solution cavity partial drop of pressure
    The water-oil phase solution cavity accelerates pressure dropIn formula:
    λm--- oil water mixture loss factor；ρm--- oil water mixture density, kg/m3；
    VVm--- the average speed of oil water mixture, m/s in solution cavity；LV--- the length of solution cavity, m；
    The equivalent diameter of d --- solution cavity, m；λ --- loss factor；
    VVam--- solution cavity exports the flow velocity of mixture, m/s；VVbm--- the flow velocity of solution cavity inlet mixture, m/s；
    αam、αbm--- solution cavity outlet, the respective kinetic energy correction factor of entrance；
    Export correction factorEntrance correction factor
    AVa、AVb--- vug unit entrance, the respective cross-sectional area in outlet, m2。
  5. The reservoir individual well pressure prediction model modeling method 5. fracture hole type according to claim 1 is connected, which is characterized in that when When being moved in solution cavity, crack for oil-gas two-phase flow:
    Pressure drop of the oil gas two-phase in vug unit along journey pressure drop, partial drop of pressure and acceleration pressure by coming down to a lower group into, the vug unit Pressure drop fundamental equation are as follows: Δ PV=Δ PVw+ΔPVp+ΔPVs；
    The pressure drop fundamental equation of the Crack Element are as follows:In formula:
    ΔPV--- vug unit overall presure drop, Pa；ΔPf--- vug unit overall presure drop, Pa；
    ΔPVw2--- oil gas two-phase solution cavity is along journey pressure drop, Pa；ΔPVp2--- oil gas two-phase solution cavity partial drop of pressure, Pa；
    ΔPVs2--- oil gas two-phase solution cavity accelerates pressure drop, Pa；Gfm--- oil-gas two-phase flow weight flow, kg/s in crack；
    μo、--- the respective viscosity of oily phase, gas phase, mPas；ρo、--- the respective density of oily phase, gas phase, kg/m3；
    Kro、Krg--- the respective relative permeability of oily phase, gas phase, μm2；
    Lf--- fracture length, m；Wf--- fracture width, mm；hf--- fracture height, m.
  6. The reservoir individual well pressure prediction model modeling method 6. fracture hole type according to claim 5 is connected, which is characterized in that
    The oil gas two-phase solution cavity is along journey pressure drop
    The oil gas two-phase solution cavity partial drop of pressure
    The oil gas two-phase solution cavity accelerates pressure drop
    Its middle outlet correction factorEntrance correction factorIn formula:
    λm--- gas mixture loss factor, the i.e. coefficient of friction resistance；ρm--- gas mixture density, kg/m3；
    GVm--- the mass flow of gas mixture, kg/s in solution cavity；AV--- solution cavity sectional area, m2；
    LV--- solution cavity length, m；D --- solution cavity radius, m；
    AVa、AVb--- vug unit entrance, the respective cross-sectional area in outlet, m2；
    VVam--- solution cavity exports the flow velocity of mixture, m/s；VVbm--- the flow velocity of solution cavity inlet mixture, m/s.
  7. 7. a kind of computer readable storage medium, the computer-readable recording medium storage has computer program, and feature exists In realization is as described in any in claim 1 to 6 the step of the method when the computer program is executed by processor.
CN201910342351.3A 2019-04-26 2019-04-26 Fracture hole type series connection reservoir individual well pressure prediction model modeling method and storage medium Pending CN110080730A (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201910342351.3A CN110080730A (en) 2019-04-26 2019-04-26 Fracture hole type series connection reservoir individual well pressure prediction model modeling method and storage medium
PCT/CN2020/086542 WO2020216313A1 (en) 2019-04-26 2020-04-24 Modeling method and storage medium of fractured-vuggy series reservoir single-well pressure prediction model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910342351.3A CN110080730A (en) 2019-04-26 2019-04-26 Fracture hole type series connection reservoir individual well pressure prediction model modeling method and storage medium

Publications (1)

Publication Number Publication Date
CN110080730A true CN110080730A (en) 2019-08-02

Family

ID=67416788

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910342351.3A Pending CN110080730A (en) 2019-04-26 2019-04-26 Fracture hole type series connection reservoir individual well pressure prediction model modeling method and storage medium

Country Status (2)

Country Link
CN (1) CN110080730A (en)
WO (1) WO2020216313A1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111581786A (en) * 2020-04-19 2020-08-25 东北石油大学 Well testing interpretation model for analyzing fracture-cave series-connection mode double-hole composite reservoir
WO2020216313A1 (en) * 2019-04-26 2020-10-29 西南石油大学 Modeling method and storage medium of fractured-vuggy series reservoir single-well pressure prediction model

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114925632B (en) * 2022-05-26 2023-09-01 西南石油大学 Dynamic simulation method for fracture-cavity type gas reservoir productivity test

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150006131A1 (en) * 2013-06-28 2015-01-01 Schlumberger Technology Corporation Methods utilizing solutal-convection regimes in a porous medium
CN107169680A (en) * 2017-06-29 2017-09-15 中国石油化工股份有限公司 A kind of method of fractured-cavernous carbonate reservoir quantification water filling
CN206972214U (en) * 2017-07-04 2018-02-06 西南石油大学 It is a kind of to simulate fracture-pore reservoir water filling, the experimental provision of gas injection displacement oil

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3616855A (en) * 1970-07-23 1971-11-02 New Mexico Tech Res Found Method of bulking or caving a volume of subsurface material
CN104747180B (en) * 2014-02-26 2016-09-14 中国石油化工股份有限公司 A kind of fracture-pore reservoir for waterflooding extraction analyzes method and application thereof
CN206592111U (en) * 2017-04-05 2017-10-27 长江大学 A kind of new fractured-cavernous carbonate reservoir displacement simulation experimental provision
WO2019103812A1 (en) * 2017-11-22 2019-05-31 Saudi Arabian Oil Company Measuring effective fracture half-length and quantifying flux distribution in and around fractures in petroleum reservoirs
CN110080730A (en) * 2019-04-26 2019-08-02 西南石油大学 Fracture hole type series connection reservoir individual well pressure prediction model modeling method and storage medium

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150006131A1 (en) * 2013-06-28 2015-01-01 Schlumberger Technology Corporation Methods utilizing solutal-convection regimes in a porous medium
CN107169680A (en) * 2017-06-29 2017-09-15 中国石油化工股份有限公司 A kind of method of fractured-cavernous carbonate reservoir quantification water filling
CN206972214U (en) * 2017-07-04 2018-02-06 西南石油大学 It is a kind of to simulate fracture-pore reservoir water filling, the experimental provision of gas injection displacement oil

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
杨瑀: "《缝洞型油藏串联式油井产能及压力预测方法》", 《中国优秀硕士学位论文全文数据库工程科技I辑》 *
陈旭: "《缝洞型油藏两相流动产能方程研究》", 《中国优秀硕士学位论文全文数据库工程科技I辑》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020216313A1 (en) * 2019-04-26 2020-10-29 西南石油大学 Modeling method and storage medium of fractured-vuggy series reservoir single-well pressure prediction model
CN111581786A (en) * 2020-04-19 2020-08-25 东北石油大学 Well testing interpretation model for analyzing fracture-cave series-connection mode double-hole composite reservoir
CN111581786B (en) * 2020-04-19 2021-02-09 东北石油大学 Well test interpretation method of well test interpretation model for analyzing fracture-cave series-mode double-hole composite reservoir

Also Published As

Publication number Publication date
WO2020216313A1 (en) 2020-10-29

Similar Documents

Publication Publication Date Title
CN110080730A (en) Fracture hole type series connection reservoir individual well pressure prediction model modeling method and storage medium
CN105840187B (en) Compactness reservoir horizontal well staged fracturing Productivity
Liu et al. Steady incompressible laminar flow in porous media
Sorbie et al. The extended Washburn equation and its application to the oil/water pore doublet problem
CN104847314B (en) HTHP oil gas straight well single-phase flow perforation completion parameter optimization method
CN104895550A (en) Tight gas fracturing horizontal well numerical value well testing model building and solving method
CN102339326A (en) Method for analyzing and simulating fluid flow of fracture-cavity oil reservoir
CN108266176A (en) A kind of natural gas wellhead method of calculating flux based on wellbore model
CN109948272A (en) Tune based on inter well connectivity blocks up dynamic prediction method and system
CN107145696A (en) A kind of analogy method of coal bed gas above and below ground couple solution
CN104504230A (en) Estimation method for recovery ratio and limit drainage radius of low-permeability gas well
CN104989385B (en) The HTHP oil gas straight well perforating parameter optimization method calculated based on skin factor
CN106150477A (en) A kind of method determining single well controlled reserves
CN104915530B (en) The method for building up of oil reservoir interwell communication relation
CN106528994B (en) A kind of pressure regulation chamber venting hole wind speed simulation method coupled based on gas-liquid interface
CN106869914A (en) The PRODUCTION FORECASTING METHODS that seepage flow is coupled with flowing in pit shaft in a kind of oil reservoir
CN110321575A (en) Gas condensate reservoir Dynamic Productivity Calculation prediction technique
CN106777663B (en) Fracturing fluid filtering stall degree calculation method considering natural fractures
CN107437127A (en) A kind of oil well stop-spraying Formation pressure prediction method
Li et al. Analysis of capillary rise in asymmetric branch-like capillary
CN204371275U (en) A kind of underground steam flow mass dryness fraction distributes regulation device
CN106468160B (en) A kind of determining CO2Drive the method and CO of foam stream oil ingredient2The analogy method of drive
Zhao et al. Characteristics and mechanisms of the tip leakage cavitating flow around a NACA66 (mod) hydrofoil under different cavitation states
CN105569625B (en) Underground steam flow mass dryness fraction distributes regulation device
CN106917614B (en) Injection well downhole flow regulator, injection well downhole flow regulator water nozzle determination method and water injection string

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20190802