一种稠油油藏水平井蒸汽吞吐初期产能预测方法 |
|||||||
申请号 | CN201710009434.1 | 申请日 | 2017-01-06 | 公开(公告)号 | CN106761733A | 公开(公告)日 | 2017-05-31 |
申请人 | 中国海洋石油总公司; 中海油研究总院; | 发明人 | 李延杰; 朱国金; 田冀; 谭先红; 袁忠超; 李南; 丁祖鹏; 卢川; 贾振; 郑强; 张晓亮; | ||||
摘要 | 本 发明 涉及一种稠油油藏 水 平井 蒸汽 吞吐初期产能预测方法,属于油气田开发领域,其包括以下步骤:1)基于渗流理论,建立稠油油藏热 力 采油数学模型;2)基于步骤1)的稠油油藏热力采油数学模型建立无因次产能曲线;3)识别影响步骤2)中无因次产能曲线的重要因素;4)绘制水平井无因次产能图版;5)基于步骤4)的水平井无因次产能图版来预测待测水平井蒸汽吞吐初期的产油量。本发明有着以下优点,根据影响水平井蒸汽吞吐产能的因素进行多元初步回归,得到水平井最大日产油量峰值计算公式;根据本发明给出的无因次产能图版,能够方便准确的预测水平井蒸汽吞吐初期产能。 | ||||||
权利要求 | 1.一种稠油油藏水平井蒸汽吞吐初期产能预测方法,包括以下步骤: |
||||||
说明书全文 | 一种稠油油藏水平井蒸汽吞吐初期产能预测方法技术领域[0001] 本发明涉及一种稠油油藏水平井蒸汽吞吐初期产能预测方法,属于油气田开发领域。 背景技术[0002] 正确预测水平井蒸汽吞吐生产动态是优化注汽工艺参数、高效开发稠油油藏的基础。目前,预测和分析水平井蒸汽吞吐产能的方法主要有数值模拟法和解析法,数值模拟能够较好地模拟蒸汽吞吐开采的动态过程,但该方法应用起来较为繁琐、费时,实用性较差;而研究预测水平井蒸汽吞吐产能的解析方法较少且由于水平井蒸汽吞吐复杂的作用机理和油藏渗流问题的存在,使得目前现有的解析法不成熟。 发明内容[0003] 针对上述问题,本发明的目的是提供一种更为实用且预测较准的稠油油藏水平井蒸汽吞吐初期产能预测方法。 [0004] 为实现上述目的,本发明采用以下技术方案:一种稠油油藏水平井蒸汽吞吐初期产能预测方法,包括以下步骤: [0005] 1)基于渗流理论,建立稠油油藏热力采油数学模型; [0006] 2)基于步骤1)的稠油油藏热力采油数学模型建立无因次产能曲线; [0007] 3)识别影响步骤2)中无因次产能曲线的重要因素; [0008] 4)绘制水平井无因次产能图版; [0009] 5)基于步骤4)的水平井无因次产能图版来预测待测水平井蒸汽吞吐初期的产油量。 [0010] 所述步骤1)中,基于渗流理论,建立稠油油藏热力采油数学模型的过程如下: [0011] ①确定稠油油藏热力采油数学模型的基本假设条件如下: [0013] b、流动满足达西定律,流动过程中各化学组分间没有发生化学反应; [0015] d、考虑重力和毛管力的影响; [0016] ②建立稠油油藏热力采油基本数学模型如下: [0017] a、质量守恒方程 [0018] [0019] 式中,k为储层渗透率,10-3μm2;krj为j相相对渗透率;μj为第j相流体粘度,mPa·s;ρj为第j相流体密度,kg/m3;xij为组分i在j相中的摩尔分数;Pj为j相流体压力,MPa;g为重力加速度,m/s2;D为标高,m;qj地层条件下,单位时间为单位体积储层中注入或采出第j相流体质量,kg/(m3·d);φ为储层孔隙度;Sj为储层中第j相流体饱和度; [0020] b、能量守恒方程 [0021] [0022] 式中,λR为油层有效导热系数,kJ/(s·m·K);ΔT为T-TR,TR是油藏温度,K;Hj为j相流体的焓,kJ/kg;为单位时间内,单位体积中与顶底层损失有关的能量,kJ/(m3·s);为单位时间内,单位体积中输入或输出的能量,kJ/(m3·s);Uj为第j相流体内能,kJ/kg;ρr为岩石密度,kg/m3;Cr为油层岩石比热,kJ/(kg·K);c、相关辅助方程如下: [0023] 饱和度方程:So+Sw+Sg=1 (3) [0024] 毛管力方程:Pc,wo=Po-Pw (4) [0025] Pc,og=Pg-Po (5) [0026] 摩尔分数归一化方程: [0027] 平衡常数方程: [0028] [0029] 式中,So为含油饱和度;Sw为含水饱和度;Sg为含气饱和度;Pc,wo为油水相间毛管力,Pa;Po为油相毛管力,Pa;Pw为水相毛管力,Pa;Pc,og为油气相间毛管力,Pa;Pg为气相毛管力,Pa;xij为第ij组分摩尔分数;T为温度,K;Koi(Pg,T)为组分i在气油相中分配时的相平衡常数;xig为组分i在气相中的摩尔分数;xio为组分i在油相中的摩尔分数;Kwi(Pg,T)为组分i在气水相中分配时的相平衡常数;xiw为组分i在水相中的摩尔分数; [0030] 式(1)~(8)构成了稠油油藏热力采油基本的数学模型,对其进行求解后即可模拟稠油热采开发效果。 [0031] 所述步骤2)中,基于步骤1)的稠油油藏热力采油数学模型建立无因次产能曲线的过程如下: [0033] 其中,油藏地质参数包括油层厚度、储层孔隙度、储层渗透率、净总厚度比、含油饱和度和原油粘度;注采参数包括注汽强度、注汽速度、井底蒸汽干度、蒸汽温度、焖井时间、采液速度和水平井长度; [0034] ②通过引入无因次参数组合将水平井蒸汽吞吐日产油数据无因次化,绘制水平井蒸汽吞吐无因次产油量qD和无因次时间的关系曲线tD,所运用的无因次参数组合见式(9)~(11); [0035] [0036] [0037] [0038] 式中,qd为水平井蒸汽吞吐无因次产油量;qo1为模拟水平井蒸汽吞吐日产油量,m3/d;qpeak1为水平井蒸汽吞吐模拟最大日产油量,m3/d;tD为水平井蒸汽吞吐无因次时间;t为时间,d;Ni为原油地质储量,104m3;A为油藏面积,m2;h为油层厚度,m;φ为储层孔隙度;Soi为油层平均原始含油饱和度;Boi为原始的原油体积系数,m3/m3。 [0039] 所述步骤3)中,识别影响步骤2)中无因次产能曲线的重要影响因素的过程如下: [0040] 油田参数分为油藏地质参数和注采参数,计算油田各参数不同取值下的无因次产能曲线并将与利用油田各参数基础取值绘制的日产油量典型曲线的偏差作为标准,通过油藏数值模拟方法研究各参数对无因次产能曲线的影响程度,以各参数不同取值下无因次产能曲线与典型曲线的偏差10%作为基准,大于10%认为影响大,小于10%认为影响小,具体判断大小的过程如下: [0041] 以典型曲线作为标准,找出各参数不同取值下无因次产能曲线下的无因次时间对应到典型曲线上的无因次产量与各参数不同取值下无因次产能曲线下无因次产量对比,计算其偏差绝对值的平均值,然后舍去对无因次产能曲线影响较小的因素,保留对其影响较大的因素。 [0042] 所述步骤4)中,绘制水平井无因次产能图版的过程如下: [0043] 参照Narayana的研究方法,绘制水平井蒸汽吞吐无因次产油量qd关于水平井蒸汽吞吐无因次时间td的无因次产能图版。 [0044] 所述步骤5)中,基于步骤4)的水平井无因次产能图版来预测待测水平井蒸汽吞吐初期的产油量的过程如下: [0045] ①以待测油井区域的实际地质和动态资料为依据,基于之前建立的无因次产能图版,查找或插值得到未投产井的无因次产能曲线,从曲线中可知随水平井蒸汽吞吐无因次时间tD变化的水平井蒸汽吞吐无因次产油量qD; [0046] ②确定最大日产油量: [0047] 最大日产油量预测值qpeak通过如下式(12)多元回归方法确定: [0048] [0049] 式中,h为油层厚度;NtG为净总厚度比;k为储层渗透率,10-3μm2;μo为原油粘度,mPa·s;qii为注汽强度,t/m;qL为产液速度,m3/d,L为水平井长度,m;a1~a9和b1~b9均为常数; [0050] ③原油地质储量Ni由式(11)计算求得,根据下述式(13)和式(14)得到水平井蒸汽吞吐日产油量qo随时间t变化的关系式 [0051] [0052] [0053] 本发明由于采取以上技术方案,其具有以下优点:1、本发明根据影响水平井蒸汽吞吐产能的因素进行多元初步回归,得到水平井最大日产油量峰值计算公式。2、本发明给出的无因次产能图版,能够方便准确的预测水平井蒸汽吞吐初期产能。附图说明 [0054] 图1是本发明储层渗透率为2000×10-3μm2时不同原油粘度的无因次产能图版; [0055] 图2是本发明储层渗透率为4000×10-3μm2时不同原油粘度的无因次产能图版; [0056] 图3是本发明储层渗透率为6000×10-3μm2时不同原油粘度的无因次产能图版; [0057] 图4是本发明的LD27-2油田A22水平井蒸汽吞吐产能预测图; [0058] 图5是本发明的NB35-2油田B33水平井多元热流体吞吐产能预测图; [0059] 图6是本发明的NB35-2油田B44水平井多元热流体吞吐产能预测图; [0060] 图7是本发明的LD27-2油田A22水平井蒸汽吞吐产能预测结果对比图; [0061] 图8是本发明的NB35-2油田B33水平井多元热流体吞吐产能预测结果对比图; [0062] 图9是本发明的NB35-2油田B44水平井多元热流体吞吐产能预测结果对比图。 具体实施方式[0063] 下面结合附图和实施例对本发明进行详细的描述。 [0064] 本发明提供的一种稠油油藏水平井蒸汽吞吐初期产能预测方法,包括以下步骤: [0065] 1)基于渗流理论,建立稠油油藏热力采油数学模型,具体过程如下: [0066] ①确定稠油油藏热力采油数学模型的基本假设条件如下: [0067] a、稠油油藏内存在油、气、水三相,共有nc种(nc-1个碳氢化合物组分和一个水组分)化学组分,任一组分在相间分配关系中服从相平衡原理; [0068] b、流动满足达西定律,流动过程中各化学组分间没有发生化学反应; [0069] c、稠油油藏内渗流为非等温渗流,温度影响原油的粘度及油、气、水三相相对渗透率; [0070] d、考虑重力和毛管力的影响。 [0071] ②建立稠油油藏热力采油基本数学模型如下: [0072] a、质量守恒方程 [0073] [0074] 式中,k为储层渗透率,10-3μm2;krj为j相相对渗透率;μj为第j相流体粘度,mPa·s;ρj为第j相流体密度,kg/m3;xij为组分i在j相中的摩尔分数;Pj为j相流体压力,MPa;g为重力加速度,m/s2;D为标高,m;qj地层条件下,单位时间为单位体积储层中注入或采出第j相流体 3 质量,kg/(m·d);φ为储层孔隙度;Sj为储层中第j相流体饱和度。 [0075] b、能量守恒方程 [0076] [0077] 式中,λR为油层有效导热系数,kJ/(s·m·K);ΔT为T-TR,TR是油藏温度,K;Hj为j相流体的焓,kJ/kg; 为单位时间内,单位体积中与顶底层损失有关的能量,kJ/(m3·s);3 为单位时间内,单位体积中输入或输出的能量,kJ/(m·s);Uj为第j相流体内能,kJ/kg; ρr为岩石密度,kg/m3;Cr为油层岩石比热,kJ/(kg·K)。 [0078] c、相关辅助方程 [0079] 饱和度方程:So+Sw+Sg=1 (3) [0080] 毛管力方程:Pc,wo=Po-Pw (4) [0081] Pc,og=Pg-Po (5) [0082] 摩尔分数归一化方程: [0083] 平衡常数方程: [0084] [0085] 式中,So为含油饱和度;Sw为含水饱和度;Sg为含气饱和度;Pc,wo为油水相间毛管力,Pa;Po为油相毛管力,Pa;Pw为水相毛管力,Pa;Pc,og为油气相间毛管力,Pa;Pg为气相毛管力,Pa;xij为第ij组分摩尔分数;T为温度,K;Koi(Pg,T)为组分i在气油相中分配时的相平衡常数;xig为组分i在气相中的摩尔分数;xio为组分i在油相中的摩尔分数;Kwi(Pg,T)为组分i在气水相中分配时的相平衡常数;xiw为组分i在水相中的摩尔分数。 [0086] 式(1)~(8)构成了稠油油藏热力采油基本的数学模型,对其进行求解后即可模拟稠油热采开发效果。 [0087] 2)基于步骤1)的稠油油藏热力采油数学模型建立无因次产能曲线,具体过程如下: [0088] ①将稠油油藏热力采油数学模型利用CMG数值模拟软件预测稠油油藏水平井蒸汽吞吐生产第一年中的模拟日产油量qo1和模拟最大日产油量qpeak1。 [0089] ②通过引入无因次参数组合将水平井蒸汽吞吐日产油数据无因次化,绘制水平井蒸汽吞吐无因次产油量qD和无因次时间的关系曲线tD,所运用的无因次参数组合见式(9)~(11)。 [0090] [0091] [0092] [0093] 式中,qD为水平井蒸汽吞吐无因次产油量;qo1为模拟水平井蒸汽吞吐日产油量,m3/d;qpeak1为水平井蒸汽吞吐模拟最大日产油量,m3/d;tD为水平井蒸汽吞吐无因次时间;t为时间,d;Ni为原油地质储量,104m3;A为油藏面积,m2;h为油层厚度,m;φ为储层孔隙度;Soi为油层平均原始含油饱和度;Boi为原始的原油体积系数,m3/m3。 [0094] 3)识别影响步骤2)中无因次产能曲线的重要因素,具体过程如下: [0095] 现有技术中油田参数分为油藏地质参数和注采参数,将油藏地质参数和注采参数纳入考虑范围,根据稠油油田开发实例统计基本油藏和地质参数取值范围,根据目前工艺水平确定注采参数取值范围。 [0096] 计算各参数不同取值下的无因次产能曲线并与利用油田各参数基础取值绘制的日产油量典型曲线的偏差作为标准,通过油藏数值模拟方法研究各参数对无因次产能曲线的影响程度,以各参数不同取值下无因次产能曲线与典型曲线的偏差10%作为基准,大于10%认为影响大,小于10%认为影响小。具体判断大小的过程如下: [0097] 以典型曲线作为标准,找出各参数不同取值下无因次产能曲线下的无因次时间对应到典型曲线上的无因次产量与各参数不同取值下无因次产能曲线下无因次产量对比,计算其偏差绝对值的平均值。然后舍去对无因次产能曲线影响较小的因素,保留对其影响较大的因素。 [0098] 4)绘制水平井无因次产能图版: [0099] 为保证无因次产能图版的准确性及使用的方便性,参照Narayana的研究方法,绘制水平井蒸汽吞吐无因次产油量qD关于水平井蒸汽吞吐无因次时间tD的无因次产能图版。 [0100] 5)基于步骤4)的水平井无因次产能图版来预测待测水平井蒸汽吞吐初期的产油量,具体过程如下: [0101] ①以待测油井区域的实际地质和动态资料为依据,基于之前建立的无因次产能图版,查找或插值得到未投产井的无因次产能曲线,从曲线中可知随水平井蒸汽吞吐无因次时间tD变化的水平井蒸汽吞吐无因次产油量qD。 [0102] ②确定最大日产油量: [0103] 最大日产油量预测值qpeak通过如下式(12)多元回归方法确定。 [0104] [0105] 式中,为NTG为净总厚度比;k为储层渗透率,10-3μm2;μo为原油粘度,mPa·s;qii为注汽强度,t/m;qL为产液速度,m3/d,L为水平井长度,m;a1~a9和b1~b9均为常数。 [0106] ③原油地质储量Ni由式(11)计算求得,根据下述式(13)和式(14)得到水平井蒸汽吞吐日产油量qo随时间t变化的关系式 [0107] [0108] [0109] 下面通过渤海LD27-2油田A22井,NB35-2油田B33井、B44井的产能预测用以说明本发明的效果,其包括以下步骤: [0110] 1)基于渗流理论,建立稠油油藏热力采油数学模型,具体过程如下: [0111] ①确定稠油油藏热力采油数学模型的基本假设条件如下: [0112] a、稠油油藏内存在油、气、水三相,共有nc种(nc-1个碳氢化合物组分和一个水组分)化学组分,任一组分在相间分配关系中服从相平衡原理; [0113] b、流动满足达西定律,流动过程中各化学组分间没有发生化学反应; [0114] c、稠油油藏内渗流为非等温渗流,温度影响原油的粘度及油、气、水三相相对渗透率; [0115] d、考虑重力和毛管力的影响。 [0116] ②建立稠油油藏热力采油基本数学模型如下: [0117] a、质量守恒方程 [0118] [0119] 式中,k为储层渗透率,10-3μm2;krj为j相相对渗透率;μj为第j相流体粘度,mPa·s;ρj为第j相流体密度,kg/m3;xij为组分i在j相中的摩尔分数;Pj为j相流体压力,MPa;g为重力加速度,m/s2;D为标高,m;qj地层条件下,单位时间为单位体积储层中注入或采出第j相流体 3 质量,kg/(m·d);φ为储层孔隙度;Sj为储层中第j相流体饱和度。 [0120] b、能量守恒方程 [0121] [0122] 式中,λR为油层有效导热系数,kJ/(s·m·K);ΔT为T-TR,TR是油藏温度,K;Hj为j相流体的焓,kJ/kg;为单位时间内,单位体积中与顶底层损失有关的能量,kJ/(m3·s);为单位时间内,单位体积中输入或输出的能量,kJ/(m3·s);Uj为第j相流体内能,kJ/kg; ρr为岩石密度,kg/m3;Cr为油层岩石比热,kJ/(kg·K)。c、相关辅助方程 [0123] 饱和度方程:So+Sw+Sg=1 (3) [0124] 毛管力方程:Pc,wo=Po-Pw (4) [0125] Pc,og=Pg-Po (5) [0126] 摩尔分数归一化方程: [0127] 平衡常数方程: [0128] [0129] 式中,So为含油饱和度;Sw为含水饱和度;Sg为含气饱和度;Pc,wo为油水相间毛管力,Pa;Po为油相毛管力,Pa;Pw为水相毛管力,Pa;Pc,og为油气相间毛管力,Pa;Pg为气相毛管力,Pa;xij为第ij组分摩尔分数;T为温度,K;Koi(Pg,T)为组分i在气油相中分配时的相平衡常数;xig为组分i在气相中的摩尔分数;xio为组分i在油相中的摩尔分数;Kwi(Pg,T)为组分i在气水相中分配时的相平衡常数;xiw为组分i在水相中的摩尔分数。 [0130] 式(1)~(8)构成了稠油油藏热力采油基本数学模型,对其进行求解后即可模拟稠油热采开发效果。 [0131] 2)基于步骤1)的稠油油藏热力采油数学模型建立无因次产能曲线,具体过程如下: [0132] ①将稠油油藏热力采油数学模型利用CMG数值模拟软件预测稠油油藏水平井蒸汽吞吐生产第一年中的模拟日产油量qo1和模拟最大日产油量qpeak1。 [0133] ②通过引入无因次参数组合将水平井蒸汽吞吐日产油数据无因次化,绘制水平井蒸汽吞吐无因次产油量qD和无因次时间的关系曲线tD,所运用的无因次参数组合见式(9)~(11)。 [0134] [0135] [0136] [0137] 式中,qD为水平井蒸汽吞吐无因次产油量;qo1为模拟水平井蒸汽吞吐日产油量,m3/d;qpeak1为水平井蒸汽吞吐模拟最大日产油量,m3/d;tD为水平井蒸汽吞吐无因次时间;t4 3 2 为时间,d;Ni为原油地质储量,10m ;A为油藏面积,m ;h为油层厚度,m;φ为储层孔隙度;Soi为油层平均原始含油饱和度;Boi为原始的原油体积系数,m3/m3。 [0138] 3)识别影响步骤2)中无因次产能曲线的重要因素,具体过程如下: [0139] 确定相关参数及取值范围,以各参数不同取值下无因次产能曲线与典型曲线的偏差作为标准,通过油藏数值模拟方法得出各参数对无因次产能曲线的影响程度如表1: [0140] 表1 稠油热采产能影响因素及取值范围 [0141] [0142] [0143] 4)绘制水平井无因次产能图版: [0144] 基于上述无因次产能曲线参数影响程度分析结果可知,对无因次产能曲线影响较大的因素为储层渗透率、原油粘度、产液速度。产液速度对无因次产能曲线的影响相对于原油粘度及储层渗透率两个因素影响较小。因此,为保证无因次产能图版的准确性及使用的方便性,参照Narayana的研究方法,将不同产液速度下的无因次产能曲线取平均值,绘制不同储层渗透率和不同原油粘度下的水平井蒸汽吞吐无因次产油量qD关于水平井蒸汽吞吐无因次时间tD的无因次产能图版。 [0145] 本发明根据已有数据绘制出储层渗透率为2000×10-3μm2、4000×10-3μm2时和[0146] 6000×10-3μm2时的不同原油粘度的无因次产能图版。数据依次如下表2~4: [0147] 表2 k=2000×10-3μm2时无因次产能图版 [0148] [0149] 表3 k=4000×10-3μm2时无因次产能图版 [0150] [0151] 表4 k=6000×10-3μm2时无因次产能图版 [0152] [0153] 基于上述储层渗透率为2000×10-3μm2、4000×10-3μm2和6000×10-3μm2时不同原油粘度的数据绘制出水平井无因次产能图版依次为图1、图2和图3。 [0154] 5)基于步骤4)的水平井无因次产能图版来预测待测水平井蒸汽吞吐初期的产油量,具体过程如下: [0155] 渤海LD27-2油田A22井,NB35-2油田B33井和B44井相关数据统计见下表5所示。 [0156] 表5 相关热采井数据统计 [0157] [0158] ①从表5可知NB35-2油田B33H、B44H井无因次产能曲线通过渗透率4000×10-3μm2、-3 26000×10 μm 和原油粘度500mPa.s、1000mPa.s无因次产能曲线插值得到,其中LD27-2油田A22井无因次产能曲线通过渗透率2000×10-3μm2、4000×10-3μm2和原油粘度1000mPa.s、 3000mPa.s无因次产能曲线插值得到。其得到的无因次产能曲线依次如图4~6,从曲线中可知随水平井蒸汽吞吐无因次时间tD变化的水平井蒸汽吞吐无因次产油量qD。 [0159] ②确定峰值日产油 [0160] 水平井蒸汽吞吐最大日产油量预测值qpeak,其通过多元回归方法确定,其公式如下式(15) [0161] [0162] 式中,NTG为净总厚度比;k为储层渗透率,10-3μm2;μo为原油粘度,mPa·s;qii为注汽强度,t/m;qL为产液速度,m3/d,L为水平井长度,m。 [0163] 通过式(14)计算得到,NB35-2油田B33H、B44H井最大日产油量分别为45.2m3/d、70.3m3/d,LD27-2油田A22井最大日产油量为80.1m3/d。 [0164] ③原油地质储量Ni由式(11)计算求得,根据下述式(13)和式(14)得到水平井蒸汽吞吐日产油量qo随时间t变化的关系式 [0165] [0166] [0167] ④将运用无因次产能图版预测的产油量与实际日产油量进行对比,得到图7~9。 [0168] 上述各实施例仅用于对本发明的目的、技术方案和有益效果进行了进一步详细说明,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。 |