CN102809526B - Method for measuring diffusion coefficient of carbon dioxide in saturated oil core - Google Patents

Method for measuring diffusion coefficient of carbon dioxide in saturated oil core Download PDF

Info

Publication number
CN102809526B
CN102809526B CN201210307962.2A CN201210307962A CN102809526B CN 102809526 B CN102809526 B CN 102809526B CN 201210307962 A CN201210307962 A CN 201210307962A CN 102809526 B CN102809526 B CN 102809526B
Authority
CN
China
Prior art keywords
diffusion
pressure
curve
carbon dioxide
eff
Prior art date
Application number
CN201210307962.2A
Other languages
Chinese (zh)
Other versions
CN102809526A (en
Inventor
李兆敏
李宾飞
王舒华
张超
李松岩
Original Assignee
中国石油大学(华东)
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 中国石油大学(华东) filed Critical 中国石油大学(华东)
Priority to CN201210307962.2A priority Critical patent/CN102809526B/en
Publication of CN102809526A publication Critical patent/CN102809526A/en
Application granted granted Critical
Publication of CN102809526B publication Critical patent/CN102809526B/en

Links

Abstract

The invention relates to a method for measuring a diffusion coefficient of carbon dioxide in a saturated oil core. The method comprises the following steps of: performing nondimensionalization by using a convection-diffusion mathematical model, and performing numerical solution, wherein the diffusion coefficient Deff of the carbon dioxide in the saturated oil core is finally obtained by using the method and required to be used in the numerical solution process; and continuously correcting the value of D'eff, and repeatedly executing the step (13), so a curve c of delta Co and t<1/2> is completely coincided with a curve b of delta PTh and the t<1/2>, wherein the D'eff is the Deff and represents the diffusion coefficient of the carbon dioxide in the saturated oil core. The method has the advantages that the influence of the volumetric expansion of crude oil caused by dissolution of the carbon dioxide in the crude oil on the diffusion process and the influence of a porous medium on the diffusion process are taken into consideration, and the obtained effective diffusion coefficient of the carbon dioxide in the saturated oil core can accurately reflect the true diffusion process.

Description

A kind of method of measuring carbon dioxide coefficient of diffusion in saturated oil rock core
Technical field
The present invention relates to a kind of method of measuring carbon dioxide coefficient of diffusion in saturated oil rock core, belong to the technical field of petrochemical complex.
Background technology
The global warming problem causing due to a large amount of discharges of carbon dioxide is increasingly severe.And collecting carbonic anhydride with bury (CO 2capture and Storage) technology thought a kind of potential, alternative carbon dioxide discharge-reduction scheme widely, can reduce the carbon dioxide content in atmosphere, alleviates climate warming problem.At carbon dioxide geological, bury in scheme, improve oil recovery factor (CO with carbon dioxide 2-EOR), not only can reach the effect of burying of carbon dioxide, and can improve the recovery ratio of Reservoir Crude Oil, therefore carbon dioxide raising oil recovery factor technology is widely used in oil field.In the process of employing carbon-dioxide flooding, the quality transfer law of carbon dioxide in saturated crude oil rock core seems particularly important for the prediction of the migration characteristics of injecting carbon dioxide.Therefore the mensuration of the coefficient of diffusion of carbon dioxide in saturated crude oil rock core has great importance for the development of carbon dioxide flooding oil tech.
Measure at present carbon dioxide method of coefficient of diffusion in porous medium and be mainly the PVT method that non-expansion type spreads, the method has only been considered the diffusion of carbon dioxide in porous medium, do not have to consider, because the dissolving of carbon dioxide causes the impact of crude oil volumetric expansion on diffusion, can not well simulate the true diffusion process of carbon dioxide in saturated crude oil rock core.Meanwhile, the PVT of employing methods are measured the coefficient of diffusion of carbon dioxide in crude oil more at present, and the measurement of this coefficient of diffusion does not reflect the impact of porous medium on diffusion process, thereby can not be used for measuring the coefficient of diffusion of carbon dioxide at saturated crude oil rock core.
Summary of the invention
Summary of the invention
For above technical deficiency, the invention provides a kind of method of measuring carbon dioxide coefficient of diffusion in saturated oil rock core.The present invention, under the combined influence effect that considers crude oil volumetric expansion and porous medium, calculates the mathematical model of carbon dioxide coefficient of diffusion in saturated crude oil rock core by having derived in conjunction with one dimension advection diffusion equation and real gas state equation.By measure the Pressure Drop of carbon dioxide in conjunction with its corresponding calculated with mathematical model the coefficient of diffusion of carbon dioxide in saturated crude oil rock core.
The diffusion model of carbon dioxide in the rock core of saturation water and the difference of the diffusion model of carbon dioxide of the present invention in the rock core of saturated oil are: carbon dioxide solubility in water after, water expands hardly, be that the expansion of water is on not impact of diffusion, and the dissolving of carbon dioxide in oil can cause significantly expansion of oil in the present invention, in rock core, the expansion of oil has obvious influence to diffusion, so it is convection current-diffusion mathematical model that the present invention will be spread mathematics model modification, but because the Analytical Solution of convection current diffusion mathematical model is more difficult, so the present invention is directed to convection current-diffusion mathematical model, first carry out nondimensionalization, then carry out numerical solution, and in numerical solution process, need to use the diffusion coefficient D of the carbon dioxide of finally trying to achieve in saturated oil rock core eff, described D effthat the method for the invention solves and obtains.
Terminological interpretation:
Saturated fluid processing: refer to by vacuumizing the whole sucking-offs of air of loading in diffusion barrel in sample rock core interior detail fine porosity, then utilize constant-flux pump that Experimental Flowing Object displacement is entered in diffusion barrel, injection experiments fluid generation build the pressure, it is logical when pressure is suppressed 10MPa, Experimental Flowing Object in diffusion barrel is under high pressure pressed in the hole of filling sample rock core, has realized whole saturated upper Experimental Flowing Objects in the hole that loads sample rock core.
Detailed Description Of The Invention
Technical scheme of the present invention is as follows:
A kind of carbon dioxide diffusion coefficient D in saturated oil rock core of measuring effmethod, comprise that step is as follows:
(1) diffusion barrel is dried;
(2) will sample cylinder rock core or artificial cylinder rock core as filling sample rock core, described filling sample rock core will be dried, after will upper and lower filling sample rock core both ends of the surface sealing with fluid sealant, pack in diffusion barrel;
(3) diffusion barrel is vacuumized;
(4) utilize well heater to heat diffusion barrel, constant temperature is to wanting simulated formation temperature, stand-by;
(5) in diffusion barrel, pump into Experimental Flowing Object, and be forced into 10MPa, to loading core sample, carry out saturated fluid processing; Described Experimental Flowing Object is simulated formation crude oil, and described simulated formation crude oil is the Simulation of Crude Oil identical with formation pore oil property;
(6) continue diffusion barrel inside Experimental Flowing Object to carry out pressurized operation, be forced into and want simulated formation pressure;
(7) in the bottom of described diffusion barrel, be provided with check valve, adjust the set pressure of check valve, described set pressure is compared with the little 0.1MPa of diffusion barrel internal pressure;
(8) to diffusion barrel, pump into carbon dioxide, the pressure of described carbon dioxide with want simulated formation pressure identical, due to the effect of check valve, the Experimental Flowing Object in emptying diffusion barrel;
(9) sealing diffusion barrel, utilizes pressure transducer to gather the original pressure P in diffusion barrel 0, gather n time point t 1~t npressure P 1~P n, n is more than or equal to 2;
(10) Δ P ex, ifor testing the pressure P of n the time point recording 1-P nrespectively with original pressure P 0pressure differential, as, Δ P ex, 1=P 1-P 0, Δ P ex, 2=P 2-P 0..., Δ P ex, n=P n-P 0; By the pressure drop Δ P calculating ex, 1~ Δ P ex, nrespectively with time point t 1~ t nthe corresponding mapping of 1/2 power, obtain the Pressure Drop Δ P that experiment records exwith t 1/2curve a, calculate the slope k of described curve a cathetus section part;
(11) to Δ P exwith t 1/2curve a revise: because carbon dioxide and the saturated crude oil at filling sample core surface place rigidly connect while touching, need to set up a stable liquid-gas interface boundary condition, thereby the pressure drop of experiment incipient stage is very fast, the Pressure Drop of experiment starting stage belongs to the pressure drop Δ P that the unstable stage of diffusion causes measurement exwith t 1/2curve map be upwards offset a segment distance, for the ease of calculating pressure fall-off curve b with theory below, contrast, need to be to Δ P exwith t 1/2curve a revise.Modification method is as follows: first find the flex point of described curve a, i.e. the corresponding time point t of curve a cathetus section part starting point s; Then by formula, 1. obtain average error
&delta; &OverBar; = &Sigma; i = s + 1 n ( &Delta; P Ex , i - k t i ) ( t i - t i - 1 ) t n - t s
Formula 1. in, Δ P ex, i-for testing the pressure P of n the time point recording 1~P nrespectively with original pressure P 0pressure differential, as, Δ P ex, 1=P 1-P 0, Δ P ex, 2=P 2-P 0..., Δ P ex, n=P n-P 0, kPa; K-be the slope of the curve a cathetus section part of trying to achieve in step (10), t i-for testing i the time point recording, min; t s-be the corresponding time point of curve a cathetus section part starting point of trying to achieve in step (10), min;
According to 2. formula to experimental pressure drop Δ P ex, irevise, obtain revising pressure drop values Δ P co, i:
&Delta; P Co , i = &Delta; P Ex , i - &delta; &OverBar; ( i = s + 1 , s + 2 , &Lambda; , n )
Make Δ P cowith t 1/2curve c;
(12) bring the slope k in step (10) into formula 3.,
&Delta; P Ex = 4 M &infin; ZRT D eff &prime; r 0 V &pi; t = k t
Formula 3. in, M -the time levels off to carbon dioxide when infinite and diffuses into the amount in filling sample rock core, mol; Z-compressibility factor is determined by the pressure and temperature of testing; R-universal gas constant, 8.314Pam 3k -1mol -1; T-experimental temperature, K; D ' effcarbon dioxide effective diffusion cofficient in-unexpansive system, m 2/ s; r 0the xsect radius of-filling sample rock core, m; Annular volume between V-filling sample rock core and diffusion barrel, m 3; T-time point, comprises t 1~ t n; K-Δ P exwith t 1/2the slope of the straight-line segment part that forms;
By formula, 3. obtained, obtain the effective diffusion cofficient D ' of carbon dioxide in unexpansive system eff;
(13) by D ' effbring in convection current-diffusion mathematical model of nondimensionalization and carry out numerical solution, obtain nondimensional gas concentration lwevel and distribute with crude oil volumetric expansion velocity distribution
(14) gas concentration lwevel is distributed with crude oil volumetric expansion velocity distribution there is bonding state equation after dimension, obtain corresponding theory and calculate pressure drop Δ P th, utilize described theory to calculate pressure drop Δ P thwith t 1/2curve plotting b, curve b and the Δ P that has just started experiment and record cowith t 1/2curve a have error, the source of error is exactly D ' eff;
(15) constantly revise D ' effvalue, repeating step (13), makes Δ P cowith t 1/2curve c and Δ P thwith t 1/2curve b overlap completely, D ' now effbe D eff, wherein said D efffor the coefficient of diffusion of carbon dioxide in saturated oil rock core.
The invention has the advantages that:
Present invention is directed at the diffusion process of carbon dioxide in saturated crude oil rock core, considered the crude oil volumetric expansion that causes due to the dissolving of carbon dioxide in the crude oil impact on diffusion process, considered again the impact of porous medium on diffusion process, the effective diffusion cofficient of the carbon dioxide of trying to achieve in saturated crude oil rock core can reflect real diffusion process comparatively accurately simultaneously.
Accompanying drawing explanation
Fig. 1 is that embodiment passes through the P in his-and-hers watches 1 1-P ncorresponding time t 1~ t nmapping, i.e. pressure temporal evolution figure;
Fig. 2 is that embodiment passes through the pressure differential deltap P in his-and-hers watches 1 1-Δ P nto time t 1/2mapping, tests pressure differential deltap P exwith t 1/2variation a, and through revised Pressure Drop Δ P cowith t 1/2curve c, by Convection-diffusion model iterative computation Δ P out thwith t 1/2curve b.
Embodiment
According to embodiment and Figure of description, the present invention is described in detail below, but is not limited to this.
Embodiment,
As shown in Figure 1-2.
A kind of carbon dioxide diffusion coefficient D in saturated oil rock core of measuring effmethod, comprise that step is as follows:
(1) diffusion barrel is dried;
(2) will sample cylinder rock core or artificial cylinder rock core as filling sample rock core, described filling sample rock core will be dried, after will upper and lower filling sample rock core both ends of the surface sealing with fluid sealant, pack in diffusion barrel;
(3) diffusion barrel is vacuumized;
(4) utilize well heater to heat diffusion barrel, constant temperature is to wanting simulated formation temperature, stand-by;
(5) in diffusion barrel, pump into Experimental Flowing Object, and be forced into 10MPa, to loading core sample, carry out saturated fluid processing; Described Experimental Flowing Object is simulated formation crude oil, and described simulated formation crude oil is the Simulation of Crude Oil identical with formation pore oil property;
(6) continue diffusion barrel inside Experimental Flowing Object to carry out pressurized operation, be forced into and want simulated formation pressure;
(7) in the bottom of described diffusion barrel, be provided with check valve, adjust the set pressure of check valve, described set pressure is compared with the little 0.1MPa of diffusion barrel internal pressure;
(8) to diffusion barrel, pump into carbon dioxide, the pressure of described carbon dioxide with want simulated formation pressure identical, due to the effect of check valve, the Experimental Flowing Object in emptying diffusion barrel;
(9) sealing diffusion barrel, utilizes pressure transducer to gather the original pressure P in diffusion barrel 0, gather n time point t 1~ t npressure P 1~P n, (n is more than or equal to 2);
(10) Δ P ex, ifor testing the pressure P of n the time point recording 1-P nrespectively with original pressure P 0pressure differential, as, Δ P ex, 1=P 1-P 0, Δ P ex, 2=P 2-P 0..., Δ P ex, n=P n-P 0; By the pressure drop Δ P calculating ex, 1~ Δ P ex, nrespectively with time point t 1~ t nthe corresponding mapping of 1/2 power, obtain the Pressure Drop Δ P that experiment records exwith t 1/2curve a, calculate the slope of described curve a cathetus section part
(11) to Δ P exwith t 1/2curve a revise: because carbon dioxide and the saturated crude oil at filling sample core surface place rigidly connect while touching, need to set up a stable liquid-gas interface boundary condition, thereby the pressure drop of experiment incipient stage is very fast, the Pressure Drop of experiment starting stage belongs to the pressure drop Δ P that the unstable stage of diffusion causes measurement exwith t 1/2curve map be upwards offset a segment distance, for the ease of calculating pressure fall-off curve b with theory below, contrast, need to be to Δ P exwith t 1/2curve a revise.Modification method is as follows: first find the flex point of described curve a, i.e. the corresponding time point t of curve a cathetus section part starting point s=19.27min; Then by formula, 1. obtain average error
&delta; &OverBar; = &Sigma; i = s + 1 n ( &Delta; P Ex , i - k t i ) ( t i - t i - 1 ) t n - t s
, &delta; &OverBar; = 245.31 ( kPa )
Formula 1. in, Δ P ex, i-for testing the pressure P of n the time point recording 1~P nrespectively with original pressure P 0pressure differential, as, Δ P ex, 1=P 1-P 0, Δ P ex, 2=P 2-P 0..., Δ P ex, n=P n-P 0, kPa; K-is the slope of the curve a cathetus section part of trying to achieve in step (10), t i-for testing i the time point recording, min; t s-be the corresponding time point of curve a cathetus section part starting point of trying to achieve in step (10), min;
According to 2. formula to experimental pressure drop Δ P ex, irevise, obtain revising pressure drop values Δ P co, i:
&Delta; P Co , i = &Delta; P Ex , i - &delta; &OverBar; ( i = s + 1 , s + 2 , &Lambda; , n )
Make Δ P cowith t 1/2curve c;
Table 1: utilize pressure transducer to gather the pressure P of each time point 1~P nwith the pressure differential deltap P calculating exwith correction pressure differential deltap P co
By the pressure in his-and-hers watches 1, the time is mapped, obtain pressure temporal evolution figure, as shown in Figure 1.
(12) bring the slope k in step (10) into formula 3.,
&Delta; P Ex = 4 M &infin; ZRT D eff &prime; r 0 V &pi; t = k t
Formula 3. in, M -the time levels off to carbon dioxide when infinite and diffuses into the amount in filling sample rock core, mol; Z-compressibility factor is determined by the pressure and temperature of testing; R-universal gas constant, 8.314Pam 3k -1mol -1; T-experimental temperature, K; D ' effcarbon dioxide effective diffusion cofficient in-unexpansive system, m 2/ s; r 0the xsect radius of-filling sample rock core, m; Annular volume between V-filling sample rock core and diffusion barrel, m 3; T-time point, comprises t 1~ t n; K-Δ P exwith t 1/2the slope of the straight-line segment part that forms;
By formula, 3. obtained, k = 4 M &infin; ZRT D eff &prime; r 0 V &pi; = 24.3033 kPa / min , Obtain the effective diffusion cofficient D ' of carbon dioxide in unexpansive system eff=2.95 × 10 -10m 2/ s;
(13) by D ' effbring in convection current-diffusion mathematical model of nondimensionalization and carry out numerical solution, obtain nondimensional gas concentration lwevel and distribute with crude oil volumetric expansion velocity distribution
&PartialD; c &OverBar; &PartialD; &tau; + c &OverBar; &PartialD; u &OverBar; &PartialD; r &OverBar; + c &OverBar; u &OverBar; r &OverBar; + &lambda; &PartialD; c &OverBar; &PartialD; r &OverBar; = &PartialD; 2 c &OverBar; &PartialD; r &OverBar; 2 c &OverBar; = 1 ( r &OverBar; = 1 , &tau; > 0 ) u &OverBar; = 0 , &PartialD; c &OverBar; &PartialD; r &OverBar; = 0 ( r &OverBar; = 0 , &tau; > 0 ) u &OverBar; = 0 , c &OverBar; = 0 ( &tau; = 0 , r &OverBar; < 1 ) u &OverBar; = 0 , c &OverBar; = 1 ( &tau; = 0 , r &OverBar; = 1 )
Formula 4. in, dimensionless time: dimensionless radius: nondimensional velocity: dimensionless concentration: c &OverBar; = c c 0 ; &lambda; = u &OverBar; - 1 r &OverBar; ;
Consider that dimensionless diffusion-convection equation fully implicit solution finite difference method that crude oil expands carries out solution procedure as follows, in departure process, the single order of speed single order differential and concentration and second-order differential central difference schemes, the single order differential of time adopts forward difference form.Above-mentioned equation is carried out to discretize, discretize form:
&PartialD; c &OverBar; &PartialD; &tau; = 1 &Delta;&tau; ( c &OverBar; i n + 1 - c &OverBar; i n ) + O ( &Delta;&tau; )
&PartialD; c &OverBar; &PartialD; r &OverBar; = 1 2 &Delta; r &OverBar; ( c &OverBar; i + 1 n + 1 - c &OverBar; i - 1 n + 1 ) + O ( &Delta; r &OverBar; 2 )
&PartialD; 2 c &OverBar; &PartialD; r &OverBar; 2 = 1 &Delta; r &OverBar; 2 ( c &OverBar; i + 1 n + 1 - 2 c &OverBar; i n + 1 + c &OverBar; i - 1 n + 1 ) + O ( &Delta; r &OverBar; 2 )
&PartialD; u &OverBar; &PartialD; r &OverBar; = 1 2 &Delta; r &OverBar; ( u &OverBar; i + 1 n + 1 - u &OverBar; i - 1 n + 1 ) + O ( &Delta; r &OverBar; 2 )
Applying finite difference scheme partial differential equation is above write as:
a i c &OverBar; i - 1 n + 1 + b i c &OverBar; i n + 1 + e i c &OverBar; i + 1 n + 1 = f i , i = 0,1,2 &Lambda; , I
Wherein: a i = - &Delta;&tau; &Delta; r &OverBar; 2 - &lambda;&Delta;&tau; 2 &Delta; r &OverBar; , b i = 1 + &Delta;&tau; 2 &Delta; r &OverBar; ( u &OverBar; i + 1 n + 1 - u &OverBar; i - 1 n + 1 ) + &Delta;&tau; r &OverBar; i u &OverBar; i n + 1 + 2 &Delta;&tau; &Delta; r &OverBar; 2 , e i = &lambda;&Delta;&tau; 2 &Delta; r &OverBar; - &Delta;&tau; &Delta; r &OverBar; 2 ,
f i = c i n .
Boundary condition discretize can obtain:
c &OverBar; - 1 n + 1 = c &OverBar; 1 n + 1 , c &OverBar; I n + 1 = 1 .
The system of equations of substitution boundary condition discretize can be write as:
b 0 a 0 + e 0 a 1 b 1 e 1 a 2 b 2 e 2 &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; a I - 1 b I = 1 c &OverBar; 0 n + 1 c &OverBar; 1 n + 1 c &OverBar; 2 n + 1 &CenterDot; &CenterDot; c &OverBar; I - 1 n + 1 = f 0 f 1 f 2 &CenterDot; &CenterDot; &CenterDot; f I - 1 - e I - 1
Your alternative manner of application Gauss-Saden calculates concentration and the velocity distribution of each time step.To each time step, by the concentration of previous time step as the velocity distribution in this time step of calculation of initial value.Adopt velocity distribution in this time step to calculate the CONCENTRATION DISTRIBUTION of next time step.Then by new CONCENTRATION DISTRIBUTION, calculate the velocity distribution in next time step.The condition that above-mentioned iterative process finishes is that the concentration maximum relative error of each radial direction Nodes is less than limits of error e (e=10 in the present invention -4).
(14) gas concentration lwevel is distributed with crude oil volumetric expansion velocity distribution there is bonding state equation after dimension, obtain corresponding theory and calculate pressure drop Δ P th, utilize described theory to calculate pressure drop Δ P thwith t 1/2curve plotting b, curve b and the Δ P that has just started experiment and record cowith t 1/2curve a have error, the source of error is exactly D ' eff:
The dimensionless gas concentration lwevel of trying to achieve in step (13) is distributed and crude oil volumetric expansion velocity distribution Δ P falls in the pressure carbon dioxide that has dimension and obtain intumescent system tht in time 1/2variation relation
&Delta;P Th = qZRT - P 0 &Delta;V V - &Delta;V
Formula 5. in, Z-compressibility factor by the pressure and temperature of testing determine; R-universal gas constant, 8.314Pam 3k -1mol -1; T-experimental temperature, K; P 0carbon dioxide original pressure in-diffusion barrel, MPa; Annular volume between V-rock sample and diffusion barrel, m 3; Δ V-obtained by 6. formula, crude oil volumetric expansion does not cause the reduction of diffusion barrel inner volume, m in the same time 3; Q-do not diffuse into the in the same time CO in saturated oil rock core 2be dissolved in the CO that expands in crude oil 2amount of substance sum, mol.
&Delta;V = 2 &pi; r 0 2 h&phi; &Sigma; &tau; = 0 &tau; u &OverBar; r &OverBar; = 1 , &tau; &Delta;&tau;
Formula 6. in, r 0the xsect radius of-filling sample rock core, m; The length of h-filling sample rock core, m; -nondimensional velocity; τ-dimensionless time.
The Pressure Drop Δ P of the Convection-diffusion model of trying to achieve thwith time relationship as following table 2,
Table 2: Δ P falls in Convection-diffusion model calculating pressure thwith time relationship table
(15) constantly revise D ' effvalue, repeating step (13), makes Δ P cowith t 1/2curve a and Δ P thwith t 1/2curve b overlap completely, D ' now effbe D eff, wherein said D efffor the coefficient of diffusion of carbon dioxide in saturated oil rock core:
By the Δ P doing in step (10) cowith t 1/2curve c and step (14) in do Δ P thwith t 1/2curve b contrast, as Fig. 2.If two curve differs larger, revise effective diffusion cofficient D ' effand get back in step (13).If two curve relative error is less than limits of error e (e=10 in this method -4), stop iterative process, effective diffusion cofficient D now eff==3.19 × 10 -10m 2/ s is the effective diffusion cofficient of carbon dioxide in saturated crude oil rock core.

Claims (1)

1. measure carbon dioxide diffusion coefficient D in saturated oil rock core for one kind effmethod, it is characterized in that, it is as follows that the method comprising the steps of:
(1) diffusion barrel is dried;
(2) will sample cylinder rock core or artificial cylinder rock core as loading core sample, described filling core sample will be dried, after will upper and lower fillings core sample both ends of the surface sealing with fluid sealant, pack in diffusion barrel;
(3) diffusion barrel is vacuumized;
(4) utilize well heater to heat diffusion barrel, constant temperature is to wanting simulated formation temperature, stand-by;
(5) in diffusion barrel, pump into Experimental Flowing Object, and be forced into 10MPa, to loading core sample, carry out saturated fluid processing; Described Experimental Flowing Object is simulated formation crude oil, and described simulated formation crude oil is the Simulation of Crude Oil identical with formation pore oil property;
(6) continue diffusion barrel inside Experimental Flowing Object to carry out pressurized operation, be forced into and want simulated formation pressure;
(7) in the bottom of described diffusion barrel, be provided with check valve, adjust the set pressure of check valve, described set pressure is compared with the little 0.1MPa of diffusion barrel internal pressure;
(8) to diffusion barrel, pump into carbon dioxide, the pressure of described carbon dioxide with want simulated formation pressure identical, due to the effect of check valve, the Experimental Flowing Object in emptying diffusion barrel;
(9) sealing diffusion barrel, utilizes pressure transducer to gather the original pressure P in diffusion barrel 0, gather n time point t 1~t npressure P 1~P n, n is more than or equal to 2;
(10) Δ P ex, ifor testing the pressure P of n the time point recording 1-P nrespectively with original pressure P 0pressure differential, as, Δ P ex, 1=P 1-P 0, Δ P ex, 2=P 2-P 0..., Δ P ex, n=P n-P 0; By the pressure drop Δ P calculating ex, 1~Δ P ex, nrespectively with time point t 1~t nthe corresponding mapping of 1/2 power, obtain the Pressure Drop Δ P that experiment records exwith t 1/2curve a, calculate the slope k of described curve a cathetus section part;
(11) to Δ P exwith t 1/2curve a revise, method is as follows: first find the flex point of described curve a, i.e. the corresponding time point t of curve a cathetus section part starting point s; Then by formula, 1. obtain average error :
&delta; &OverBar; = &Sigma; i = s + 1 n ( &Delta;P Ex , i - k t i ) ( t i - t i - 1 ) t n - t s ?????①
Formula 1. in, Δ P ex, i-for testing the pressure P of n the time point recording 1~P nrespectively with original pressure P 0pressure differential, as, Δ P ex, 1=P 1-P 0, Δ P ex, 2=P 2-P 0..., Δ P ex, n=P n-P 0, kPa; K-be the slope of the curve a cathetus section part of trying to achieve in step (10), t i-for testing i the time point recording, min; t s-be the corresponding time point of curve a cathetus section part starting point of trying to achieve in step (10), min;
According to 2. formula to experimental pressure drop Δ P ex, irevise, obtain revising pressure drop values Δ P co, i:
&Delta; P Co , i = &Delta;P Ex , i - &delta; &OverBar; , ( i = s + 1 , s + 2 , . . . , n ) ???②
Make Δ P cowith t 1/2curve c;
(12) bring the slope k in step (10) into formula 3.,
&Delta;P Ex = 4 M &infin; ZRT D eff &prime; r 0 V &pi; t = k t ????③
Formula 3. in, M -the time levels off to carbon dioxide when infinite and diffuses into the amount of loading in core sample, mol; Z-compressibility factor, is determined by the pressure and temperature of testing; R-universal gas constant, 8.314Pam 3k -1mol -1; T-experimental temperature, K; D ' effcarbon dioxide effective diffusion cofficient in-unexpansive system, m 2/ s; r 0the xsect radius of-filling core sample, m; Annular volume between V-filling core sample and diffusion barrel, m 3; T-time point, comprises t 1~t n; K-Δ P exwith t 1/2the slope of the straight-line segment part that forms;
By formula, 3. obtained, obtain the effective diffusion cofficient D ' of carbon dioxide in unexpansive system eff;
(13) by D ' effbring in convection current-diffusion mathematical model of nondimensionalization and carry out numerical solution, obtain nondimensional gas concentration lwevel and distribute with crude oil volumetric expansion velocity distribution ;
(14) gas concentration lwevel is distributed with crude oil volumetric expansion velocity distribution have after dimension in conjunction with real gas state equation, obtain corresponding theory and calculate pressure drop Δ P th, utilize described theory to calculate pressure drop Δ P thwith t 1/2curve plotting b, curve b and the Δ P that has just started experiment and record cowith t 1/2curve a have error, the source of error is exactly D ' eff;
(15) constantly revise D ' effvalue, repeating step (13), makes Δ P cowith t 1/2curve c and Δ P thwith t 1/2curve b overlap completely, D ' now effbe D eff, wherein said D efffor the coefficient of diffusion of carbon dioxide in saturated oil rock core.
CN201210307962.2A 2012-08-27 2012-08-27 Method for measuring diffusion coefficient of carbon dioxide in saturated oil core CN102809526B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210307962.2A CN102809526B (en) 2012-08-27 2012-08-27 Method for measuring diffusion coefficient of carbon dioxide in saturated oil core

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210307962.2A CN102809526B (en) 2012-08-27 2012-08-27 Method for measuring diffusion coefficient of carbon dioxide in saturated oil core

Publications (2)

Publication Number Publication Date
CN102809526A CN102809526A (en) 2012-12-05
CN102809526B true CN102809526B (en) 2014-04-16

Family

ID=47233287

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210307962.2A CN102809526B (en) 2012-08-27 2012-08-27 Method for measuring diffusion coefficient of carbon dioxide in saturated oil core

Country Status (1)

Country Link
CN (1) CN102809526B (en)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104975827B (en) * 2014-04-03 2018-01-19 中国石油化工股份有限公司 Predict the material balance method of carbon dioxide flooding oil reservoir index
CN104237079B (en) * 2014-10-10 2017-01-18 中国石油大学(华东) Experimental method and device for measuring diffusion coefficient of carbon dioxide in water
CN104502236B (en) * 2014-12-15 2017-03-22 中国石油大学(华东) Method for measuring diffusion coefficient and equilibrium concentration of CO2 in process of diffusion from water phase to oil phase
WO2017181392A1 (en) * 2016-04-21 2017-10-26 深圳市樊溪电子有限公司 Method for measuring diffusion coefficient of carbon dioxide in crude oil
CN107314950B (en) * 2016-04-27 2019-11-05 中国石油化工股份有限公司 A method of measurement carbon dioxide diffusion coefficient in porous media
CN106596353B (en) * 2016-12-19 2019-02-26 中国科学院武汉岩土力学研究所 Gas diffusivity measuring device and method in coal cinder based on differential pressure method
CN106840973B (en) * 2017-02-16 2019-06-04 中国石油大学(华东) CO in a kind of test porous media2The device and its test method of diffusion concentration and diffusion coefficient
CN106872594B (en) * 2017-02-16 2018-11-16 中国石油大学(华东) CO in a kind of test saturated oils porous media2The method of concentration distribution and diffusion coefficient
CN107063935B (en) * 2017-04-21 2018-07-17 中国石油大学(华东) A kind of measurement CO2Diffusion process and the method for diffusion coefficient between water-oil phase
CN108801860B (en) * 2017-05-02 2021-01-01 中国石油化工股份有限公司 Method for measuring diffusion coefficient of carbon dioxide in crude oil
CN108362614B (en) * 2018-01-19 2020-03-27 中国石油大学(华东) Device and method for measuring diffusion coefficient in huff and puff process of shale oil carbon dioxide
CN109269943B (en) * 2018-11-27 2021-04-20 衡阳师范学院 Method for rapidly measuring effective diffusion coefficient of radon in film
CN109283098B (en) * 2018-11-27 2021-04-20 衡阳师范学院 Method for measuring effective diffusion coefficient of radon in film
CN109883891B (en) * 2019-03-25 2020-03-27 中国石油大学(北京) Measuring device and method for injected gas-liquid phase diffusion coefficient of fracture-cavity oil reservoir

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4915903B2 (en) * 2006-01-31 2012-04-11 国立大学法人横浜国立大学 Method and apparatus for measuring oxygen diffusion coefficient of porous material
CN101929939B (en) * 2010-08-06 2013-03-13 中联煤层气国家工程研究中心有限责任公司 Method for detecting diffusion coefficient of coal bed gas in coal matrix
CN101975718B (en) * 2010-08-13 2012-11-14 中国科学院山西煤炭化学研究所 Method for simultaneously measuring high-pressure gas adsorption capacity and adsorption swell capacity of coal petrography and measuring equipment
CN102121891A (en) * 2010-12-16 2011-07-13 西南石油大学 Device and method for testing molecular diffusion coefficients in multi-component oil-gas system under high temperature and high pressure
CN202195998U (en) * 2011-09-16 2012-04-18 东北石油大学 Device for testing diffusion coefficient of natural gas

Also Published As

Publication number Publication date
CN102809526A (en) 2012-12-05

Similar Documents

Publication Publication Date Title
CN104101564B (en) A kind of method of unstable state High Temperature High Pressure test flow in low permeability core free-boundary problem
Gasda et al. Vertical equilibrium with sub-scale analytical methods for geological CO 2 sequestration
Clarkson et al. Production analysis of tight-gas and shale-gas reservoirs using the dynamic-slippage concept
CN104502224B (en) Saturation water Coal Under rock isothermal desorption curve determination device and method
Shook et al. Determining reservoir properties and flood performance from tracer test analysis
CN103323352B (en) Natural gas hydrate deposit dynamic triaxial mechanic-acoustic-electrical synchronous test experimental device and method
CN103246820B (en) A kind of numerical simulation for oil-gas reservoir computing method
CN106021778B (en) A kind of simulation CO2The determination method of displacement dynamic miscible pressure
Li et al. Diffusion coefficients of supercritical CO2 in oil-saturated cores under low permeability reservoir conditions
CN103257151B (en) In a kind of quantitative evaluation oil and gas secondary migration process, pore throat employs the method for rule
CN104018829B (en) A kind of method that utilization coal bed gas well creation data measures air water phase percolation curve
Liu et al. Coal permeability evolution and gas migration under non-equilibrium state
Li et al. Production behavior of methane hydrate in porous media using huff and puff method in a novel three-dimensional simulator
Connell et al. History matching of enhanced coal bed methane laboratory core flood tests
CN103278615B (en) Test method for geological storage of carbon dioxide coal seams
CN202300351U (en) Long rock core model displacement device
CN202330233U (en) Experiment test device for permeability of rock core under condition of formation pressure
CN103279647A (en) Method for calculating tight sandstone reservoir permeability based on hole feature parameters
CN104568694B (en) Method for testing gas-water relative permeability of dense core
CN203630018U (en) Unstable-state determination device for permeability of hyposmosis rock
CN105840187A (en) Method for calculating staged fracturing productivity of compact reservoir horizontal well
CN104895550B (en) A kind of tight gas pressure break horizontal well numerical well testing model establishes method for solving
CN102914494B (en) Device for measuring dynamic leak-off of foam fracturing fluid and working method thereof
CN103527185A (en) Horizontal-well physical simulation experiment device and experimental method thereof
CN102353625B (en) Method for measuring overburden porosity with water in permeation fluid mechanics experiment

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant