CN106468160B - A kind of determining CO2Drive the method and CO of foam stream oil ingredient2The analogy method of drive - Google Patents
A kind of determining CO2Drive the method and CO of foam stream oil ingredient2The analogy method of drive Download PDFInfo
- Publication number
- CN106468160B CN106468160B CN201510501873.5A CN201510501873A CN106468160B CN 106468160 B CN106468160 B CN 106468160B CN 201510501873 A CN201510501873 A CN 201510501873A CN 106468160 B CN106468160 B CN 106468160B
- Authority
- CN
- China
- Prior art keywords
- bubble
- phase
- relative permeability
- component
- oil
- 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
Abstract
A kind of determining CO2Drive the method and CO of foam stream oil ingredient2The analogy method of drive, wherein determining CO2The method for driving foam stream oil ingredient content includes: parameter determination, determines thermodynamic parameter, relative permeability and each component of foam stream oil and the mole percent of bubble；Component determines step, the parameter according to obtained in parameter determination, based on default compositional model, determines the content of each component and bubble.This method utilizes note CO2The default compositional model of process " foam oil stream " has been accurately realized in gas injection process " foam oil stream " distribution and the simulation of behavioral characteristics.
Description
Technical field
The present invention relates to oil-gas exploration and development technical fields, specifically, being related to a kind of determining CO2Drive foam stream oil ingredient
Method and CO2The analogy method of drive.
Background technique
China major part oil field belongs to continental deposit, and serious heterogeneity, viscosity of crude are higher.For this kind of oil field, water drive
Recovery ratio can only be up to 33% or so, therefore recovery ratio has the potentiality that are greatly improved.In existing raising recovery efficiency technique, note
Gas drive is that the another important channel of oil recovery is improved after water drive, polymer flooding, steam drive.CO2Drive with it is low in cost,
It achieves noticeable achievement, recycling, asepsis environment-protecting can be recycled and have the advantages that preferable compatibility with crude oil and be listed in raising and adopt
The first choice of yield.
CO2Driving " foam oil stream " is CO2A kind of phenomenon of displacement asphaltenes crude oil.This phenomenon is in Xiang Gaohan
Asphaltic crude injects CO2In the process, CO2It dissolves each other with crude oil and mixes the phenomenon of immiscible region domain formation.In displacement process
In, in displacing front, the gas continuously flowed is formed " foam during similar thickened oil recovery by Asphaltene adsorption function influence
The bubble of oil stream " can delay CO in this way2Coning improves oil recovery.
A large amount of note CO is carried out both at home and abroad at present2Oil recovery research is improved, wherein relating generally to CO2It is mixed phase, immiscible
Phase displacement phase and analysis of seepage flow mechanism have also related to the research of asphaltenes crude oil studies on asphaltene precipitation deposition, but at present
Model research be difficult to simulate CO2The flow characteristics of the formation of bubble and " foam oil stream " in displacement process.
Summary of the invention
To solve the above problems, the present invention provides a kind of application seepage simulations to determine CO2Drive foam stream oil ingredient content
Method, which comprises
Parameter determination determines thermodynamic parameter, relative permeability and each component and bubble of foam stream oil
Mole percent；
Component determines step, is determined each according to parameter obtained in the parameter determination based on default compositional model
The content of a component and bubble.
According to one embodiment of present invention, the default compositional model includes:
The hydrocarbon component equation, water flowing equation, total hydrocarbon equation, fugacity equation, kinetics equation and volume constraint equation.
According to one embodiment of present invention, for the default compositional model, come really by the way of adaptive implicit
The content of fixed each component and bubble.
According to one embodiment of present invention,
The hydrocarbon component equation indicates are as follows:
The water flowing equation indicates are as follows:
The total hydrocarbon equation indicates are as follows:
The fugacity equation indicates are as follows:
lnfig=lnfioI=1,2 ..., nc
The kinetics equation indicates are as follows:
React 1:CO2+ light hydrocarbon component → bubble
Reaction speed x1:x1=F1 × [CO2]+F2 × [light hydrocarbon component]
Reaction 2: bubble → CO2+ light hydrocarbon component
Reaction speed x2:x2=F3 × [bubble]
The volume constraint equation indicates are as follows:
Wherein, Tg、To、TwAnd TbGas phase, oily phase, the conductivity of water phase and bubble are respectively indicated, variable m is equal to n or n+
1, yoi、ygiAnd ybiRespectively indicate molar fraction of i-th component in oily phase, gas phase and bubble, Φg、Φo、ΦwAnd ΦbRespectively
Indicate that the pressure potential of gas phase, oily phase, water phase and bubble, V indicate molal volume, NiIndicate the molal quantity of the i-th component, Ng、No、Nw
And NbThe molal quantity of gas phase, oily phase, water phase and bubble is respectively indicated, F1, F2 and F3 indicate that coefficient of frequency, [A] are rubbing for substance A
That score, ρ indicate molar density, figAnd fioFugacity of i-th component in gas phase and oily phase is respectively indicated, φ indicates porosity,
ncIndicating that the sum of component, w indicate water phase, o indicates oily phase, and g indicates gas phase, and b indicates that bubble, t indicate that total hydrocarbon, q indicate yield,
Time step and new time step before n and n+1 is respectively indicated.
According to one embodiment of present invention, the default compositional model is carried out by the way of adaptive implicit
The process of solution includes:
Equation in the default compositional model is subjected to differencing, obtained difference equation is expanded into iteration surplus
Form, obtain remainder vector；
According to the remainder vector, Jacobian matrix, the content of each component is determined.
According to one embodiment of present invention, the remainder vector includes:
R=[R1,R2,...R2nc+3]-1
Wherein,
According to one embodiment of present invention, the content of each component is determined according to following expression:
J δ x=-R
Wherein, J indicates that Jacobian matrix, p indicate pressure, NiIndicate the molal quantity of the i-th component, NigIt indicates i-th in gas phase
The molal quantity of component, NwIndicate the molal quantity of water phase, i ∈ [1, nc], ncIndicate that the sum of component, R indicate remainder vector.
According to one embodiment of present invention, according to the relative permeability of water phase, relative permeability, the phase of bubble of oily phase
To the relative permeability of permeability and gas phase, conductivity is determined.
According to one embodiment of present invention, the step of determining the relative permeability of the bubble include:
According to critical saturation, the initial flow feature of bubble is determined；
Measure bubble highest relative permeability；
Relative permeability based on gas phase, it is true according to the initial flow feature of the bubble and bubble highest relative permeability
The relative permeability of the fixed bubble.
According to one embodiment of present invention, determine that relative permeability, the relative permeability of oily phase of the water phase are gentle
The step of relative permeability of phase includes:
According to two-phase oil-water system, the relative permeability of the water phase is determined；
The relative permeability that normalization determines the gas phase is demarcated by critical saturation；
Based on three-phase relative permeability mode, according to the relative permeability of the relative permeability of the water phase and gas phase, really
The relative permeability of the fixed oily phase.
The present invention also provides a kind of CO2The analogy method of drive, which comprises
According to reservoir geologic character, geological model is established, and grid dividing is carried out to the geological model；
Each component and the content of bubble in foam stream oil are determined in each grid using method as described above；
According to each component of foam stream oil in each grid and the content of bubble, determine that the oil reservoir is carrying out
CO2During drive, the behavioral characteristics of foam stream oil.
The present invention provides " foam oil stream " distribution characteristics in a kind of simulation gas injection process and the influences to recovery ratio is improved
Method, belong to Eor By Gas Injection analogy method field.Method provided by the present invention based on compositional model, in conjunction with
The advantages of " balance model " and " dynamic model ", regards the flow behavior of " foam oil stream " as the time and introduces flowing item
The function of part.The key of processing is that dynamic model describes " foam oil stream " flow behavior and the variation of time is closed by establishing
System considers that the flowing of bubble in " foam oil stream " is the flowing different from continuous gas phase, also inconsistent with the flowing velocity of crude oil,
The flowing that " bubble " has been handled in foam oil by interpolation and normalized method mutually seeps average value with gas phase by seeking bubble
Inverse oil mutually mutually seeps.Model is solved by adaptive implicit method.The present invention has developed note CO2Asphaltenes crude oil improves harvesting
Rate score analog theory and method have been put forward for the first time note CO2The mathematical model of process " foam oil stream ", science realize gas injection
" foam oil stream " distribution and the simulation of behavioral characteristics in journey.
Other features and advantages of the present invention will be illustrated in the following description, also, partly becomes from specification
It obtains it is clear that understand through the implementation of the invention.The objectives and other advantages of the invention can be by specification, right
Specifically noted structure is achieved and obtained in claim and attached drawing.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is required attached drawing in technical description to do simple introduction:
Fig. 1 is CO according to an embodiment of the invention2The flow chart of the analogy method of drive；
Fig. 2 is the bubble of crude oil to be analyzed according to an embodiment of the invention to permeability and gas phase relative permeability
Compare figure；
Fig. 3 is the curve graph of average gas relative permeability according to an embodiment of the invention；
Fig. 4 is note CO according to an embodiment of the invention2Asphaltenes oil reservoir CO in the process2And the gas of " foam oil stream "
The distribution characteristics schematic diagram of bubble；
Fig. 5 is according to an embodiment of the invention using note CO2Consider foam oil stream model and using note CO2It does not examine
Consider the comparison diagram of the counted crude oil cumulative production of foam oil stream model meter；
Fig. 6 is according to an embodiment of the invention using note CO2Consider foam oil stream model and using note CO2It does not examine
Consider the comparison diagram of the counted gas-oil ratio of foam oil stream model meter.
Specific embodiment
Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings and examples, how to apply to the present invention whereby
Technological means solves technical problem, and the realization process for reaching technical effect can fully understand and implement.It needs to illustrate
As long as not constituting conflict, each feature in each embodiment and each embodiment in the present invention can be combined with each other,
It is within the scope of the present invention to be formed by technical solution.
Meanwhile in the following description, for illustrative purposes and numerous specific details are set forth, to provide to of the invention real
Apply the thorough understanding of example.It will be apparent, however, to one skilled in the art, that the present invention can not have to tool here
Body details or described ad hoc fashion are implemented.
In addition, step shown in the flowchart of the accompanying drawings can be in the department of computer science of such as a group of computer-executable instructions
It is executed in system, although also, logical order is shown in flow charts, and it in some cases, can be to be different from herein
Sequence execute shown or described step.
Fig. 1 shows CO provided by the present embodiment2The flow chart of the analogy method of drive.
As shown in Figure 1, CO provided by the present embodiment2The analogy method of drive is in step s101 according to oil reservoir to be analyzed
Matter feature establishes the geological model of oil reservoir to be analyzed, and carries out grid dividing to geological model.In the present embodiment, to geology
When model carries out grid dividing, roughening treatment has also been carried out.
In step s 102, it is determined in each grid according to the fluid parameter of oil reservoir to be analyzed using default compositional model
The component of foam stream oil.Specifically, in the present embodiment, reservoir fluid well stream object composition to be analyzed is determined by spectrograph, and
Gum asphaltic content in crude oil is determined using asphaltene precipitation experiment.It is determined by out the pitch of asphaltenes crude oil
Matter drop out point and amount of precipitation, and simulate to obtain asphalitine Precipitating pressure according to airwater mist cooling phase equilibrium model and calculate pitch
Matter amount of precipitation.
In the present embodiment, default compositional model includes the hydrocarbon component equation, water flowing equation, total hydrocarbon equation, fugacity equation, moves
Mechanical equation and volume constraint equation.Specifically, the hydrocarbon component equation can indicate are as follows:
Water flowing equation can indicate are as follows:
Total hydrocarbon equation can indicate are as follows:
Fugacity equation can indicate are as follows:
lnfig=ln fioI=1,2 ..., nc (4)
Kinetics equation can indicate are as follows:
React 1:CO2+ light hydrocarbon component → bubble (5)
Reaction speed x1:x1=F1 × [CO2]+F2 × [light hydrocarbon component] (6)
Reaction 2: bubble → CO2+ light hydrocarbon component (7)
Reaction speed x2:x2=F3 × [bubble] (8)
Volume constraint equation can indicate are as follows:
Wherein, Tg、To、TwAnd TbGas phase, oily phase, the conductivity of water phase and bubble are respectively indicated, variable m is equal to n or n+
1, yoi、ygiAnd ybiRespectively indicate molar fraction of i-th component in oily phase, gas phase and bubble, Φg、Φo、ΦwAnd ΦbRespectively
Indicate that the pressure potential of gas phase, oily phase, water phase and bubble, V indicate molal volume, NiIndicate the molal quantity of the i-th component, Ng、No、Nw
And NbThe molal quantity of gas phase, oily phase, water phase and bubble is respectively indicated, F1, F2 and F3 indicate that coefficient of frequency, [A] are rubbing for substance A
That score, ρ indicate molar density, figAnd fioFugacity of i-th component in gas phase and oily phase is respectively indicated, φ indicates porosity,
ncIndicating that the sum of component, w indicate water phase, o indicates oily phase, and g indicates gas phase, and b indicates that bubble, t indicate that total hydrocarbon, q indicate yield,
Time step and new time step before n and n+1 is respectively indicated.
When utilize CO2When carrying out displacement crude oil, CO2, existing uneven conversion process, this meeting between solution gas and free gas
The obvious over-saturation of solution gas in oily phase is caused, to postpone the discharge of solution gas, and forms the low transparent bubble of bubble point more dynamic than really heat
Point.
Above-mentioned non-equilibrium process is influenced by bubble nucleating kinetic factor and gas diffusibility factor, bubble thus at
Core effect depends on the length of nucleation time.Therefore, the macroscopic property for relying solely on state equation to solve bubble is can not
It meets the requirements.
By analysis, it is found that the formation of bubble also receives the adhesion of asphaltene colloid.And this effect is more approximate
In the process of chemical reaction.Therefore, the present invention is analyzed by macroscopic property of the kinetic model to bubble.
Specifically, in the present embodiment, it is assumed that the bubble in foam oil stream was mainly made of three parts, it may be assumed that light hydrocarbon component,
Bubble and CO2, the conversion between them can get up by kinetics equation simultaneous:
In the present embodiment, to CO2Flash calculation is carried out with crude oil, to obtain the composition, molecule of gas phase after flash calculation
The thermodynamic parameters such as amount, the compressed coefficient, critical pressure and critical-temperature and bubble mole percent yb。
In the present embodiment, due to individually handling using bubble as pseudocomponent, and bubble is stream before forming free gas
Dynamic, bubbly flow is different from continuous gas phase, therefore the relative permeability of bubble need to be handled individually.In addition, bubble is critical full
It is lower than the critical saturation of other gas components with spending.Therefore, in the present embodiment, the relative permeability of bubble is seeped into curve
Critical saturation individually defines.
Specifically, in the present embodiment, by defining critical saturation, to determine the initial flow feature of bubble.Due to gas
Bubble mainly flows together with crude oil, therefore the critical saturation of bubble is determined as to the critical saturation of oil.In addition, bubble is to infiltration
Saturating rate is smaller than gas, therefore, in the present embodiment, by measuring bubble highest relative permeability, two-phase flow section in proportion
Normalization is handled to determine the permeability saturation curve of bubble, to obtain the relative permeability of bubble.Specifically, the present embodiment
In, the figure compared with gas phase relative permeability is as shown in Figure 2 to permeability for the bubble of crude oil to be analyzed.
In gas and oil water three-phase system, the relative permeability of water phase can be determined by two-phase oil-water system, water phase
Relative permeability is the function of water phase saturation.But gas phase is made of two parts, therefore in order to solve the relative permeability of oily phase,
Determine that gas phase mean permeability is crucial.Average gas phase relative permeability is determined as bubble pair by method provided by the present embodiment
The weighted average of the mole percent of permeability and gas phase relative permeability.Specifically, in the present embodiment, crude oil to be analyzed
Average gas phase relative permeability krgAs shown in Figure 3.Finally, can determine whole system by three-phase relative permeability model
The relative permeability of oily phase.
It should be noted that in other embodiments of the invention, can also determine foam using other rational methods
The relative permeability of each phase in stream oil, the invention is not limited thereto.
In the present embodiment, above-mentioned default compositional model is solved using adaptive fully implicit method, it is each to obtain
The molal quantity of component.
Specifically, in the present embodiment, difference is carried out to 2n+3 equation included in default compositional model, and will obtain
The expansion of 2n+3 difference, and sequence is carried out in the form of iteration surplus, to obtain remainder vector R.In the present embodiment, remainder
It include 2n+3 element in vector R, i.e. remainder vector R can be indicated are as follows:
Wherein, in the present embodiment, the 1st to the n-th in remainder vector RcA element is the redundancy value obtained according to fugacity equation,
N-thc+ 1 element is the redundancy value obtained according to kinetics equation, n-thc+ 2 to 2ncA element is according to component mobility equation
Obtained redundancy value, 2nc+ 1 element is the redundancy value obtained according to total hydrocarbon equation, 2nc+ 2 elements are according to water flowing
The redundancy value that equation obtains, 2nc+ 3 elements are the redundancy value obtained according to volume constraint equation.
That is, work as z=1,2 ..., ncWhen, exist:
Rz=ln fzg-lnfzo (11)
Work as z=ncWhen+1, exist:
Work as z=nc+2,...,2ncWhen, exist:
Work as z=2ncWhen+1, exist:
Work as z=2ncWhen+2, exist:
Work as z=2ncWhen+3, exist:
In this way, sharing 2n in above-mentioned expression formula (11)~(16)c+ 3 known variables, i.e., the solution of default compositional model.This
In embodiment, the solution vector δ x of default compositional model can be indicated are as follows:
Using Jacobian matrix J, solution vector δ x can be determined according to following expression:
J δ x=-R (18)
In the present embodiment, Jacobian matrix J can be indicated are as follows:
In the present embodiment, the constraint condition and flowing bottom pressure of well have been also contained in remainder vector R and solution vector δ x, together
When, above-mentioned whole system can pass through Newton Algorithm.
Again as shown in Figure 1, analogy method provided by the present embodiment is in step s 103 according to foam in each grid
The component for flowing oil can determine that oil reservoir is carrying out CO2Foam flows behavioral characteristics during drive.
The present embodiment institute providing method solves injection CO2High asphaltenes crude oil forms foam oil stream and is difficult to mould in the process
Quasi- the problem of calculating, foam oil stream and CO can be calculated by this method2Flow characteristics and injection CO2The recovery ratio of crude oil afterwards.
Fig. 4 is that note CO is calculated by above method2Asphaltenes oil reservoir CO in the process2And point of the bubble of " foam oil stream "
Cloth feature schematic diagram.From figure 3, it can be seen that in CO2Injection initial stage, bubble are concentrated mainly on CO2The leading edge of displacement, this is also
Illustrate that the bubble of similar foam oil stream plays the role of stable injection gas front.
Fig. 5 is shown using note CO2Consider foam oil stream model and using note CO2Do not consider that foam oil stream model calculates
The comparison diagram of the crude oil cumulative production obtained.From figure 5 it can be seen that considering compared with the analysis result for not considering foam oil stream
The displacement efficiency of oil stream model significantly improves, and is that there are a large amount of bubbles in displacing front mainly at injection initial stage
When, oil production is apparently higher than that there is no when bubble.
Fig. 6 is shown using note CO2Consider foam oil stream model and using note CO2Do not consider that foam oil stream model calculates
The comparison diagram of the gas-oil ratio obtained.From fig. 6 it can be seen that considering the gas-oil ratio of foam oil stream model obviously than not considering foam oil
It flows low.This illustrates CO2There is foam oil stream and stabilizes displacing front in displacing front, meanwhile, reduce the tempo of penetration of gas, from
And make CO2It is more abundant with crude oil mixed phase, and then improve displacement efficiency.
Method provided in this embodiment is capable of determining that asphaltenes crude oil infuses CO2Foam oil stream and CO2Distribution characteristics,
Different schemes crude output, gas liquid ratio can also be calculated using this method.Shown by numerical simulation study and is utilized now
The calculating crude oil cumulative oil production for having method to analyze is compared, and the crude oil cumulative oil production analyzed using this method is more
Height, this also indicates that displacement efficiency improves very significantly.Meanwhile the gas-oil ratio obtained using this method is obviously more existing than utilizing
The obtained gas-oil ratio of method wants low, this illustrates CO2The appearance of displacing front foam oil stream stabilizes displacing front, meanwhile, subtract
The tempo of penetration for having lacked gas, makes CO2It is more abundant with crude oil mixed phase, to improve displacement efficiency.
The present invention provides " foam oil stream " distribution characteristics in a kind of simulation gas injection process and the influences to recovery ratio is improved
Method, belong to Eor By Gas Injection analogy method field.The present invention based on compositional model, in conjunction with " balance model " and
The advantages of " dynamic model ", regards the flow behavior of " foam oil stream " as the time and introduces the function of flox condition.Place
The key of reason is to describe the variation relation of " foam oil stream " flow behavior and time, consideration " foam by establishing dynamic model
The flowing of bubble is the flowing different from continuous gas phase in oil stream ", also inconsistent with the flowing velocity of crude oil, passes through interpolation and returns
One method changed has handled in foam oil the flowing of " bubble ", and average value inverse oil phase phase is mutually seeped with gas phase by seeking bubble
It seeps.Model is solved by adaptive implicit method.
The present invention has developed note CO2Asphaltenes crude oil improves recovery ratio simulation theory and method, is put forward for the first time
Infuse CO2The mathematical model of process " foam oil stream ", science realize in gas injection process " foam oil stream " distribution and behavioral characteristics
Simulation.
It should be understood that disclosed embodiment of this invention is not limited to specific structure disclosed herein, processing step
Or material, and the equivalent substitute for these features that those of ordinary skill in the related art are understood should be extended to.It should also manage
Solution, term as used herein is used only for the purpose of describing specific embodiments, and is not intended to limit.
" one embodiment " or " embodiment " mentioned in specification means the special characteristic described in conjunction with the embodiments, structure
Or characteristic is included at least one embodiment of the present invention.Therefore, the phrase " reality that specification various places throughout occurs
Apply example " or " embodiment " the same embodiment might not be referred both to.
For convenience, as used herein multiple to may alternatively appear in common list.However, these lists should be interpreted that the list
In each element be respectively identified as independent unique member.Therefore, in the absence of a negative explanation, do not have in the list
One member can be based only upon them and appear in the practical etc. of any other member for being just interpreted same list in common list
Jljl.In addition, can also come together herein together with the substitution for each element referring to various embodiments of the present invention and example.It answers
When understanding, these embodiments, example and substitution are not construed as mutual equivalent, and are considered of the invention independent
Autonomous representative.
Although above-mentioned example is used to illustrate principle of the present invention in one or more application, for the technology of this field
For personnel, without departing from the principles and ideas of the present invention, hence it is evident that can in form, the details of usage and implementation
It is upper that various modifications may be made and does not have to make the creative labor.Therefore, the present invention is defined by the appended claims.
Claims (7)
1. a kind of application seepage simulation determines CO2The method for driving foam stream oil ingredient content, which is characterized in that the described method includes:
Parameter determination, determine foam stream oil thermodynamic parameter, relative permeability and each component and bubble moles hundred
Score；
Component determines step, determines each group based on default compositional model according to parameter obtained in the parameter determination
Point and bubble content；
Wherein, the default compositional model includes: the hydrocarbon component equation, water flowing equation, total hydrocarbon equation, fugacity equation, dynamics
Equation and volume constraint equation；
The hydrocarbon component equation indicates are as follows:
The water flowing equation indicates are as follows:
The total hydrocarbon equation indicates are as follows:
According to the relative permeability of water phase, relative permeability, the relative permeability of the relative permeability of bubble and gas phase of oily phase,
Determine conductivity；
The step of determining the relative permeability of the bubble include:
According to critical saturation, the initial flow feature of bubble is determined；
Measure bubble highest relative permeability；
Relative permeability based on gas phase determines institute according to the initial flow feature of the bubble and bubble highest relative permeability
State the relative permeability of bubble；
The fugacity equation indicates are as follows:
lnfig=lnfioI=1,2 ..., nc
Based on CO2The formation mechenism for driving foam stream bubble in oil establishes the kinetics equation, and by the kinetics equation come
Solve the macroscopic property of bubble；
The CO2Drive the formation mechenism of foam stream bubble in oil are as follows: utilize CO2When carrying out displacement crude oil, CO2, solution gas and free gas
Between there is uneven conversion process, which causes solution gas over-saturation described in oily phase, thus described in postponement
The discharge of solution gas, and form specific heat and move the low transparent bubble point of bubble point；The non-equilibrium process be bullied soak nuclear dynamics factor and
The influence of gas diffusion sexual factor, bubble is so nucleation depends on the length of nucleation time；
The kinetics equation indicates are as follows:
React 1:CO2+ light hydrocarbon component → bubble
Reaction speed x1:x1=F1 × [CO2]+F2 × [light hydrocarbon component]
Reaction 2: bubble → CO2+ light hydrocarbon component
Reaction speed x2:x2=F3 × [bubble]
The volume constraint equation indicates are as follows:
Wherein, Tg、To、TwAnd TbGas phase, oily phase, the conductivity of water phase and bubble are respectively indicated, variable m is equal to n or n+1, yoi、
ygiAnd ybiRespectively indicate molar fraction of i-th component in oily phase, gas phase and bubble, Φg、Φo、ΦwAnd ΦbRespectively indicate gas
The pressure potential of phase, oily phase, water phase and bubble, V indicate molal volume, NiIndicate the molal quantity of the i-th component, Ng、No、NwAnd NbRespectively
Indicate that the molal quantity of gas phase, oily phase, water phase and bubble, F1, F2 and F3 indicate that coefficient of frequency, [A] are the molar fraction of substance A, ρ
Indicate molar density, figAnd fioFugacity of i-th component in gas phase and oily phase is respectively indicated, φ indicates porosity, ncExpression group
The sum divided, w indicate water phase, and o indicates oil phase, and g indicates gas phase, b expression bubble, and t indicates that total hydrocarbon, q indicate yield, and n and n+1 divide
It Biao Shi not preceding time step and new time step.
2. the method as described in claim 1, which is characterized in that for the default compositional model, using adaptive implicit
Mode determines the content of each component and bubble.
3. method according to claim 1 or 2, which is characterized in that come for described default by the way of adaptive implicit
The process that compositional model is solved includes:
Equation in the default compositional model is subjected to differencing, obtained difference equation is expanded into the shape of iteration surplus
Formula obtains remainder vector；
According to the remainder vector, Jacobian matrix, the content of each component is determined.
4. method as claimed in claim 3, which is characterized in that the remainder vector includes:
Wherein,
Wherein, RzIndicate z-th of element in remainder vector R.
5. method as claimed in claim 4, which is characterized in that determine the content of each component according to following expression:
J δ x=-R
Wherein, J indicates that Jacobian matrix, p indicate pressure, NiIndicate the molal quantity of the i-th component, NigIndicate the i-th component in gas phase
Molal quantity, NwIndicate the molal quantity of water phase, i ∈ [1, nc], ncIndicate that the sum of component, R indicate remainder vector.
6. the method as described in claim 1, which is characterized in that determine the relative permeability of the water phase, the opposite infiltration of oily phase
The step of relative permeability of rate and gas phase includes: thoroughly
According to two-phase oil-water system, the relative permeability of the water phase is determined；
The relative permeability that normalization determines the gas phase is demarcated by critical saturation；
Institute is determined according to the relative permeability of the relative permeability of the water phase and gas phase based on three-phase relative permeability mode
State the relative permeability of oily phase.
7. a kind of CO2The analogy method of drive, which is characterized in that the described method includes:
According to reservoir geologic character, geological model is established, and grid dividing is carried out to the geological model；
Using method such as according to any one of claims 1 to 6 determine in each grid in foam stream oil each component and
The content of bubble；
According to each component of foam stream oil in each grid and the content of bubble, determine that the oil reservoir is carrying out CO2It drives
In the process, the behavioral characteristics of foam stream oil.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510501873.5A CN106468160B (en) | 2015-08-14 | 2015-08-14 | A kind of determining CO2Drive the method and CO of foam stream oil ingredient2The analogy method of drive |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510501873.5A CN106468160B (en) | 2015-08-14 | 2015-08-14 | A kind of determining CO2Drive the method and CO of foam stream oil ingredient2The analogy method of drive |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106468160A CN106468160A (en) | 2017-03-01 |
| CN106468160B true CN106468160B (en) | 2019-10-11 |
Family
ID=58213702
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510501873.5A Active CN106468160B (en) | 2015-08-14 | 2015-08-14 | A kind of determining CO2Drive the method and CO of foam stream oil ingredient2The analogy method of drive |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106468160B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108694255B (en) * | 2017-04-06 | 2021-06-15 | 中国石油化工股份有限公司 | Supercritical gas phase displacement medium-condensate oil gas equilibrium phase behavior evaluation method |
| CN107143317B (en) * | 2017-05-31 | 2020-09-22 | 西南石油大学 | Low-permeability reservoir air flooding numerical simulation method and device |
| CN109374824B (en) * | 2017-08-04 | 2021-06-15 | 中国石油化工股份有限公司 | Testing supercritical CO2Method for condensing the phase behavior of oil and gas systems |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103246820A (en) * | 2013-05-21 | 2013-08-14 | 中国石油大学(北京) | Oil and gas reservoir numerical simulation calculation method |
-
2015
- 2015-08-14 CN CN201510501873.5A patent/CN106468160B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103246820A (en) * | 2013-05-21 | 2013-08-14 | 中国石油大学(北京) | Oil and gas reservoir numerical simulation calculation method |
Non-Patent Citations (1)
| Title |
|---|
| 注CO2含沥青质重质原油沉积机理研究;贾英;《中国优秀博硕士学位论文全文数据库（硕士） 工程科技I辑》;20070115;第6.2、6.3节 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106468160A (en) | 2017-03-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110334431B (en) | Single-well control reserve calculation and residual gas analysis method for low-permeability tight gas reservoir | |
| CN105019894B (en) | A kind of stratified reservoir inter well connectivity method for establishing model and system | |
| Alvarez et al. | Unified model for steady-state foam behavior at high and low foam qualities | |
| CN105740563B (en) | Preferential channel identification method for secondary development of mature oil field | |
| CN106547938B (en) | Crack-pore structure dual media coal seam reservoirs air-water two phase flow method for numerical simulation | |
| Li et al. | An analytical method for modeling and analysis gas-water relative permeability in nanoscale pores with interfacial effects | |
| CN106761733A (en) | A kind of horizontal wells in heavy oil reservoir steam soak initial productivity Forecasting Methodology | |
| CN104504230A (en) | Estimation method for recovery ratio and limit drainage radius of low-permeability gas well | |
| CN104915530B (en) | The method for building up of oil reservoir interwell communication relation | |
| CN106468160B (en) | A kind of determining CO2Drive the method and CO of foam stream oil ingredient2The analogy method of drive | |
| Abbaszadeh et al. | Experimentally-based empirical foam modeling | |
| CN103775078A (en) | Method for forecasting effective viscosity of power law fluid in porous medium | |
| CN116306385B (en) | Oil reservoir fracturing imbibition energy increasing numerical simulation method, system, equipment and medium | |
| CN107621350A (en) | One kind simulation supercritical CO2The method of displacement natural gas flow | |
| CN109209307A (en) | A kind of method of quantitative analysis waterflood development of low-permeability reservoirs effect | |
| CN113836695B (en) | Oil reservoir numerical simulation method based on gridless connecting element | |
| Luo et al. | Numerical simulation of the impact of polymer rheology on polymer injectivity using a multilevel local grid refinement method | |
| Hagemann et al. | Multicomponent gas rising through water with dissolution in stratified porous reservoirs–Application to underground storage of H2 and CO2 | |
| Christman et al. | Comparison of laboratory-and field-observed CO2 tertiary injectivity | |
| CN107729587A (en) | Consider the method for numerical simulation of thinner injection | |
| John et al. | A new generation chemical-flooding simulator | |
| Nasriani et al. | Investigation into the effect of capillary number on productivity of a lean gas condensate reservoir | |
| Xu et al. | A fractal physics-based data-driven model for water-flooding reservoir (FlowNet-fractal) | |
| Liu et al. | From molecular dynamics to lattice Boltzmann: a new approach for pore-scale modeling of multi-phase flow | |
| CN110263382A (en) | The prediction technique and system of vertical gas injection crack oil and gas water microscopic seepage rule |
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 |