专利名称：一种测量二氧化碳在饱和油岩心中扩散系数的方法
技术领域：
本发明涉及一种测量二氧化碳在饱和油岩心中扩散系数的方法，属于石油化工的技术领域。
背景技术：
由于二氧化碳大量排放而弓I发的全球气候变暖问题日趋严峻。而二氧化碳捕集与埋存(CO2 Capture and Storage)技术已被广泛地认为是一种潜在的、可供选择的二氧化碳减排方案，能减少大气中的二氧化碳含量，缓解气候变暖问题。在二氧化碳地质埋存方案中，用二氧化碳提高原油采收率(CO2-EOR),不仅可以达到二氧化碳的埋存效果，而且可以
提高油藏原油的采收率，因此二氧化碳提高原油采收率技术被广泛的应用于油田。在采用二氧化碳驱油的过程中，二氧化碳在饱和原油岩心中的质量传递规律对于注入二氧化碳的运移特征的预测显得尤为重要。因此二氧化碳在饱和原油岩心中的扩散系数的测定对于二氧化碳驱油技术的发展具有重要的意义。目前测量二氧化碳在多孔介质中扩散系数的方法主要为非膨胀型扩散的PVT方法，该方法仅考虑了二氧化碳在多孔介质中的扩散，没有考虑由于二氧化碳的溶解导致原油体积膨胀对扩散的影响，不能很好的模拟二氧化碳在饱和原油岩心中的真实扩散过程。同时，目前多采用PVT方法测量二氧化碳在原油中的扩散系数，该扩散系数的测量没有反映多孔介质对扩散过程的影响，因而也不能用于测量二氧化碳在饱和原油岩心中的扩散系数。

发明内容
发明概述针对以上的技术不足，本发明提供一种测量二氧化碳在饱和油岩心中扩散系数的方法。本发明在综合考虑原油体积膨胀和多孔介质的综合影响作用下，通过结合一维径向对流扩散方程与真实气体状态方程推导了计算二氧化碳在饱和原油岩心中扩散系数的数学模型。通过测量二氧化碳气体的压力降结合其相应的数学模型计算了二氧化碳在饱和原油岩心中的扩散系数。二氧化碳在饱和水的岩心中的扩散模型与本发明所述的二氧化碳在饱和油的岩心中的扩散模型的区别在于二氧化碳溶解在水中之后，水几乎不膨胀，即水的膨胀对扩散没有影响，而本发明中二氧化碳在油中的溶解会导致油明显的膨胀作用，岩心中油的膨胀作用对扩散有明显的影响作用，所以本发明将扩散数学模型修改为对流-扩散数学模型，但是由于对流扩散数学模型的解析求解比较困难，所以本发明针对对流-扩散数学模型先进行无量纲化，然后进行数值求解，而数值求解过程中需要用到最终求得的二氧化碳在饱和油岩心中的扩散系数Drff，所述的Deff是本发明所述方法求解得到的。术语解释饱和流体处理是指通过抽真空将扩散筒中装填样品岩心内部细小孔隙中的空气全部吸出，然后利用平流泵将实验流体驱替入扩散筒中，注入实验流体并产生憋压，通当压力憋到IOMPa时，使得扩散筒中的实验流体在高压下被压入装填样品岩心的孔隙中，实现了装填样品岩心的孔隙中全部饱和上实验流体。发明详述本发明的技术方案如下一种测量二氧化碳在饱和油岩心中扩散系数Deff的方法，包括步骤如下(I)对扩散筒进行烘干； (2)将取样圆柱岩心或人造圆柱岩心作为装填样品岩心，将所述的装填样品岩心烘干，用密封胶将装填样品岩心上下两端面密封后装入扩散筒中；(3)将扩散筒抽真空；(4)利用加热器对扩散筒进行加热，恒温至要模拟地层温度，待用；(5)向扩散筒中泵入实验流体，并加压至lOMPa，对填装岩心样品进行饱和流体处理；所述的实验流体为模拟地层原油，所述模拟地层原油是与地层孔隙原油性质相同的模拟原油；(6)继续对扩散筒内部实验流体进行加压操作，加压至要模拟地层压力；(7)在所述扩散筒的底部设置有回压阀，调整回压阀的设定压力，所述的设定压力较扩散筒内部压力小O. IMPa ；(8)向扩散筒泵入二氧化碳，所述二氧化碳的压力与要模拟地层压力相同，由于回压阀的作用，排空扩散筒中的实验流体；(9)封闭扩散筒，利用压力传感器采集扩散筒内的初始压力P。、采集η个时间点ti tn的压力P1 Pn, η大于等于2 ；(10) APeju为实验测得的η个时间点的压力P1-Pn分别与初始压力P。的压力差，
如，Δ PEx,!=P1-P0, APex 2=P2-P0........ΔΡΕχ η=Ρη-Ρ0 ;将计算得到的压降 Λ ΡΕχ,广 Λ ΡΕχ,η 分别
与时间点t^tn的1/2次幂对应作图，得到实验测得的压力降APex与t1/2的曲线a，计算所述曲线a中直线段部分的斜率k ；(11)对Λ Pex与？/2的曲线a进行修正由于二氧化碳与装填样品岩心表面处的饱和原油刚接触时需要建立一个稳定的气液界面边界条件，因而实验开始阶段的压降很快，实验初始阶段的压力降属于扩散不稳定阶段导致测量的压降APex与t1/2的曲线图向上偏移了一段距离，为了便于后面与理论计算压降曲线b进行对比，需要对APex与七1/2的曲线a进行修正。修正方法如下首先找到所述曲线a的拐点，即曲线a中直线段部分起点所对应的时间点ts ;然后由公式①求出平均误差孑 ,―厂参 r,,')①
K ~(s在公式①中，ΛΡΕχ -为实验测得的η个时间点的压力P1 Pn分别与初始压力P。的压力差，如，APex1=P1-P0, APex 2=P2-P0,……、ΛPExn=Pn-PQ，kPa ;k—为步骤(10)中求得的曲线a中直线段部分的斜率，kPa / Vmin ； ti一为实验测得的第i个时间点，min ；ts一为步骤(10)中求得的曲线a中直线段部分起始点所对应的时间点，min ；根据②式对实验压降APEx i修正，得到修正压力降值APct^i
A/ = AZVv; - δ (7 = λ. + l”s'+ 2，Λ ,/ )②作出Λ Pco与t1/2的曲线c ;(12)将步骤(10)中的斜率k带入式③，APbx JM ^ zrt^ = k^t 在式③中，M00—时间趋近于无穷时二氧化碳扩散进入装填样品岩心中的量，mo I ；Z—压缩因子由实验的压力和温度确定；R—通用气体常数，8 . 314Pa · m3 · Γ1 · moF1 ;T一实验温度，K ;D' eff—非膨胀体系中二氧化碳有效扩散系数，m2/s 装填样品岩心的横截面半径，m ;V-装填样品岩心与扩散筒之间的环空体积，m3 ;t—时间点，包括ftn
t1/2所形成直线段部分的斜率；
π γ JcY由式③得求出非膨胀体系中二氧化碳的有效扩散系数
IO M , /J i I
D' eff ；(13)将D' rff带入无量纲化的对流-扩散数学模型中进行数值求解，得到无量纲的二氧化碳浓度分布[和原油体积膨胀速度分布 ;( 14)对二氧化碳浓度分布[和原油体积膨胀速度分布 进行有量纲化后结合状态方程，求出相对应的理论计算压降Λ PTh，利用所述的理论计算压降八？ 1与沪绘制曲线比曲线b与刚开始实验测得的八匕与七1/2的曲线a是有误差的，误差的来源就是D' eff ；(15)不断修改D' eff的值，重复步骤(13)，使得ΛPc。与t1/2的曲线c与APtt与t1/2的曲线b完全重合，此时的IV eff即为Drff，其中所述的Deff为二氧化碳在饱和油岩心中的扩散系数。本发明的优点在于本发明针对于二氧化碳在饱和原油岩心中的扩散过程，即考虑了由于二氧化碳在原油中的溶解导致的原油体积膨胀对扩散过程的影响，同时又考虑了多孔介质对扩散过程的影响，所求得的二氧化碳在饱和原油岩心中的有效扩散系数能较为准确的反映真实的扩散过程。


图I为实施例通过对表I中的P1-Pn对应时间t^tn作图，即压力随时间变化图；图2为实施例通过对表I中的压差AP1-APn对时间t1/2作图，即实验压差Λ PEx随t1/2的变化a，以及经过修正后的压力降Λ P。。与t1/2的曲线C，通过对流-扩散模型迭代计算出来的的曲线b。
具体实施例方式下面根据实施例和说明书附图对本发明做详细的说明，但不限于此。实施例、如图1-2所示。一种测量二氧化碳在饱和油岩心中扩散系数Deff的方法，包括步骤如下(I)对扩散筒进行烘干；
(2)将取样圆柱岩心或人造圆柱岩心作为装填样品岩心，将所述的装填样品岩心烘干，用密封胶将装填样品岩心上下两端面密封后装入扩散筒中；(3)将扩散筒抽真空；(4)利用加热器对扩散筒进行加热，恒温至要模拟地层温度，待用；(5)向扩散筒中泵入实验流体，并加压至lOMPa，对填装岩心样品进行饱和流体处理；所述的实验流体为模拟地层原油，所述模拟地层原油是与地层孔隙原油性质相同的模拟原油；(6)继续对扩散筒内部实验流体进行加压操作，加压至要模拟地层压力；(7)在所述扩散筒的底部设置有回压阀，调整回压阀的设定压力，所述的设定压力 较扩散筒内部压力小O. IMPa ； (8)向扩散筒泵入二氧化碳，所述二氧化碳的压力与要模拟地层压力相同，由于回压阀的作用，排空扩散筒中的实验流体；(9)封闭扩散筒，利用压力传感器采集扩散筒内的初始压力P。、采集η个时间点t广tn的压力P1 Pn, (η大于等于2)；(10) APeju为实验测得的η个时间点的压力P1-Pn分别与初始压力P。的压力差，
如，APe^1=P1-P0, APex 2=P2-P0........ΔΡΕχ η=Ρη-Ρ0 ;将计算得到的压降 ΛΡΕχ,广ΛΡΕχ,η 分别
与时间点t^tn的1/2次幂对应作图，得到实验测得的压力降APex与t1/2的曲线a，计算所述曲线a中直线段部分的斜率Λ=24.3033 kPa/Vmin ；(11)对Λ Pex与？/2的曲线a进行修正由于二氧化碳与装填样品岩心表面处的饱和原油刚接触时需要建立一个稳定的气液界面边界条件，因而实验开始阶段的压降很快，实验初始阶段的压力降属于扩散不稳定阶段导致测量的压降APex与t1/2的曲线图向上偏移了一段距离，为了便于后面与理论计算压降曲线b进行对比，需要对APex与七1/2的曲线a进行修正。修正方法如下首先找到所述曲线a的拐点，即曲线a中直线段部分起点所对应的时间点ts=19. 27min ;然后由公式①求出平均误差:
Σ(^pExJ — ) — tI-I) ^j __①
t — t
n s得，5= 245.31 (kPa)在公式①中，ΛΡΕχ -为实验测得的η个时间点的压力P1 Pn分别与初始压力P。的压力差，如，APex1=P1-P0, APex 2=P2-P0,……、ΛPEx n=Pn-PQ，kPa ;k_ 为步骤(10)中求得的曲线a中直线段部分的斜率，kPa/vWn ； tj一为实验测得的第i个时间点，min ；ts一为步骤(10)中求得的曲线a中直线段部分起始点所对应的时间点，min ；根据②式对实验压降APEx i修正，得到修正压力降值APcm Δ/,, .. 二 — (,· = S + I, IV + 2, A . n)②作出Λ Pco与t1/2的曲线c ;表I :利用压力传感器采集各时间点的压力P1 Pn和计算得到的压差APex和修正压差ΔΡ。。时 _ I FR力 I im APex I lJii APco | 温度 | 时_ | &力 |tm APex I tk.S APco I 温度
min MPa KPa__KPa qC min MPaKPa__KPaqC
08.170 O O 70 94 7.690480 23570
18.110 60 O 70 103 7.680490 24570
37,980 190__O__70 HO 7.670500__255__70
47.940 230__O__70 119 7.660510__265__70
57.920 250 O 70 128 7.650520 27570
77.890 280 O 70 13B 7.640530 28570
87.880 290 O 70 147 7.630540 29570 LO 7.860 310 O 70 157 7.620550 305 70
12 7.840 330__O__70 168 7.610560__315__70
15 7.830 340__O__70 179 7.600570__325__70
17 7.820 350__O__70 189 7.590580__335__70
19 7.810 360 O 70 201 7.580590 34570
23 7.800 370 125 70 212 7.570600 35570
30 7.790 380 135 70 225 7.560610 36570
34 7.780 390__145__70 237 7.550620__375__70
42 7.770 400 155 70 250 7,540630 38570
46 7.760 410__165__70 264 7.530640__395__70
52 7.750 420__175__70 275 7.520650__405__70
57 7.740 430__185__70 289 7.510660__415__70
64 7.730 440 195 70 304 7.500670 42570
72 7.720 450 205 70 318 7.490680 43570
78 7.710 460 215 70 333 7.480690__44570
86 7.70047022570通过对表I中的压力对时间作图，得到压力随时间变化图，见图I所示。( 12)将步骤(10)中的斜率k带入式③，A/ =4琴& =咖
r0V4π在式③中，M00-时间趋近于无穷时二氧化碳扩散进入装填样品岩心中的量，mo I ；Z-压缩因子由实验的压力和温度确定；R_通用气体常数，8. 314Pa πιοΓ1 ；T-实验温度，K ;D' rff—非膨胀体系中二氧化碳有效扩散系数，m2/s 一装填样品岩心的横截面半径，m;V-装填样品岩心与扩散筒之间的环空体积，m3;t-时间点，包括t广tn;k-APE!^t1/2所形成直线段部分的斜率；
4V, 7 R i [Jy ~_由式③得，灸=:'j ,上“1'1' =24.3033 kPa，求出非膨胀体系中二氧化碳
的有效扩散系数D' eff=2. 95 X 10—1 V/s ；(13)将D' rff带入无量纲化的对流-扩散数学模型中进行数值求解，得到无量纲的二氧化碳浓度分布]和原油体积膨胀速度分布w ;
权利要求
1.一种测量二氧化碳在饱和油岩心中扩散系数Deff的方法，其特征在于，该方法包括步骤如下 (1)对扩散筒进行烘干； (2)将取样圆柱岩心或人造圆柱岩心作为装填样品岩心，将所述的装填样品岩心烘干，用密封胶将装填样品岩心上下两端面密封后装入扩散筒中； (3)将扩散筒抽真空； (4)利用加热器对扩散筒进行加热，恒温至要模拟地层温度，待用； (5)向扩散筒中泵入实验流体，并加压至lOMPa，对填装岩心样品进行饱和流体处理；所述的实验流体为模拟地层原油，所述模拟地层原油是与地层孔隙原油性质相同的模拟原油； (6)继续对扩散筒内部实验流体进行加压操作，加压至要模拟地层压力； (7)在所述扩散筒的底部设置有回压阀，调整回压阀的设定压力，所述的设定压力较扩散筒内部压力小0. IMPa ； (8)向扩散筒泵入二氧化碳，所述二氧化碳的压力与要模拟地层压力相同，由于回压阀的作用，排空扩散筒中的实验流体； (9)封闭扩散筒，利用压力传感器采集扩散筒内的初始压力P。、采集n个时间点压力P1 Pn，n大于等于2 ； (10)APeju为实验测得的n个时间点的压力P1-Pn分别与初始压力Ptl的压力差，如，A PEx,!=P1-Po^ APex^2=P2-P0........APex n=Pn-P0 ;将计算得到的压降 A PEx,广 A PEx,n 分别与时间点t^tn的1/2次幂对应作图，得到实验测得的压力降APex与t1/2的曲线a，计算所述曲线a中直线段部分的斜率k ； (11)对APex与t1/2的曲线a进行修正，方法如下首先找到所述曲线a的拐点，即曲线a中直线段部分起点所对应的时间点ts ;然后由公式①求出平均误差孑 ^(^Pexj Xe - tI-I)。 t -t n s 在公式①中，APex,广为实验测得的n个时间点的压力P1 Pn分别与初始压力P。的压力差，如，APex1=P1-P0, APex 2=P2-P0,……、APEx n=Pn-PQ，kPa ;k_ 为步骤(10)中求得的曲线a中直线段部分的斜率，kPa/V^ : tr为实验测得的第i个时间点，min ；ts-为步骤(10)中求得的曲线a中直线段部分起始点所对应的时间点，min ； 根据②式对实验压降APEx i修正，得到修正压力降值APco i A"( = — S ( = ^ + I, s -I- 2，八 '")② 作出AP。。与#的曲线c ; (12)将步骤(10)中的斜率k带入式③， 在式③中，M00—时间趋近于无穷时二氧化碳扩散进入装填样品岩心中的量，mo I ；Z—压缩因子由实验的压力和温度确定；R—通用气体常数，8. SHPa-Hi3-T1 - mo F1 ;T—实验温度，K ;D' eff—非膨胀体系中二氧化碳有效扩散系数，m2/s ;r(1 一装填样品岩心的横截面半径，m ;V-装填样品岩心与扩散筒之间的环空体积，m3 ;t—时间点，包括t:tn ；k-APEx与t1/2所形成直线段部分的斜率； 由式③得"=告，求出非膨胀体系中二氧化碳的有效扩散系数D' eff; (13)将D'rff带入无量纲化的对流-扩散数学模型中进行数值求解，得到无量纲的二氧化碳浓度分布[和原油体积膨胀速度分布《 ; (14)对二氧化碳浓度分布卩和原油体积膨胀速度分布^进行有量纲化后结合状态方 程，求出相对应的理论计算压降A PTh，利用所述的理论计算压降八？^与沪绘制曲线比曲线b与刚开始实验测得的AP。。与t1/2的曲线a是有误差的，误差的来源就是D' eff ； (15)不断修改D'rff的值，重复步骤(13)，使得AP。。与t1/2的曲线c与APlh与t1/2的曲线b完全重合，此时的D' eff即为Drff，其中所述的Drff为二氧化碳在饱和油岩心中的扩散系数。
全文摘要
本发明涉及一种测量二氧化碳在饱和油岩心中扩散系数的方法，本发明利用对流-扩散数学模型，通过对流-扩散数学模型先进行无量纲化，然后进行数值求解，而数值求解过程中需要用到最终求得的二氧化碳在饱和油岩心中的扩散系数Deff，所述的Deff是本发明所述方法求解得到的。不断修改D′eff的值，重复步骤(13)，使得ΔCo与t1/2的曲线c与ΔPTh与t1/2的曲线b完全重合，此时的D′eff即为Deff，其中所述的Deff为二氧化碳在饱和油岩心中的扩散系数。本发明即考虑了由于二氧化碳在原油中的溶解导致的原油体积膨胀对扩散过程的影响，同时又考虑了多孔介质对扩散过程的影响，所求得的二氧化碳在饱和原油岩心中的有效扩散系数能较为准确的反映真实的扩散过程。
文档编号G01N13/04GK102809526SQ20121030796
公开日2012年12月5日 申请日期2012年8月27日 优先权日2012年8月27日
发明者李兆敏, 李宾飞, 王舒华, 张超, 李松岩 申请人:中国石油大学(华东)