CN106761621B - 一种聚合物驱三层窜流油藏试井解释参数的获取方法 - Google Patents
一种聚合物驱三层窜流油藏试井解释参数的获取方法 Download PDFInfo
- Publication number
- CN106761621B CN106761621B CN201710072135.2A CN201710072135A CN106761621B CN 106761621 B CN106761621 B CN 106761621B CN 201710072135 A CN201710072135 A CN 201710072135A CN 106761621 B CN106761621 B CN 106761621B
- Authority
- CN
- China
- Prior art keywords
- layers
- oil reservoir
- pressure
- channelling
- formula
- 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
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Abstract
本发明涉及一种聚合物驱三层窜流油藏试井解释参数的获取方法,包括以下步骤:1)根据聚合物溶液的流变特性,确定聚合物溶液粘度模型修正公式,并根据建立聚合物驱三层窜流油藏物理模型;2)根据聚合物溶液粘度模型的修正公式和聚合物驱三层窜流油藏物理模型,确定聚合物驱三层窜流油藏数学模型;3)采用有限差分方法求解聚合物驱三层窜流油藏数学模型,得到井底压力的数值解;4)对井筒储集系数、地层系数、储容比、时间和井底压力进行无量纲化,并绘制聚合物驱三层窜流油藏典型曲线理论图版;5)将聚合物驱三层窜流油藏典型曲线理论图版与油田实测数据曲线进行拟合,得到聚合物驱三层窜流油藏的储层渗透率、表皮因子和井筒储集系数。
Description
技术领域
本发明涉及一种聚合物驱三层窜流油藏试井解释参数的获取方法,属于试井技术领域。
背景技术
聚合物驱,特别是部分水解聚丙烯酰胺(HPAM)聚合物驱,其由于具有低成本、高效率等优点,因此是最常用的一种提高油田采收率的方法。聚合物驱提高油田采收率的主要机理是通过改善多孔介质的粘性指进来提高波及体积。但是,在聚合物通过储层岩石时,由于其吸附滞留会引起储层伤害或者渗透率下降,因此,聚驱后需要采用相应的评价储层伤害的技术。
试井方法是获取储层渗透率和表皮因子的主要方法,进而可以评价储层伤害程度。1960年,最早提出了水驱多层油藏的试井分析方法。随着计算机技术的发展,不同模型和算法都相继进行了研究,如单层复合模型、双层窜流油藏模型和数值反演技术等。虽然关于复合油藏和双层窜流油藏的相应试井技术有了一定进展,但是到目前为止,绝大部分试井解释模型仅仅把聚合物溶液流变特性简单地描述为幂律型流体,而忽略了聚合物和地层岩石之间的扩散和对流等作用,并且也没有试井解释模型能够解释聚驱三层窜流油藏,造成大庆、辽河和渤海等许多油田测试数据难以解释,导致油田管理人员不能及时了解储层情况并采取相应措施,进而影响聚合物驱油效果。
发明内容
针对上述问题,本发明的目的是提供一种聚合物驱三层窜流油藏试井解释参数的获取方法,该方法能够准确地对油田实测数据进行解释,并为具有类似储层条件的聚合物油藏试井解释、聚合物驱后储层评价提供依据。
为实现上述目的,本发明采取以下技术方案:一种聚合物驱三层窜流油藏试井解释参数的获取方法,包括以下步骤:1)根据聚合物溶液的流变特性,确定聚合物溶液粘度模型的修正公式,并根据聚合物驱三层窜流油藏特征,建立聚合物驱三层窜流油藏物理模型;2)根据步骤1)中获得的聚合物溶液粘度模型的修正公式和聚合物驱三层窜流油藏物理模型,确定聚合物驱三层窜流油藏数学模型;3)采用有限差分方法求解步骤2)中获得的聚合物驱三层窜流油藏数学模型,得到井底压力的数值解;4)对井筒储集系数、地层系数、储容比、时间和步骤3)中获得的井底压力分别进行无量纲化,并绘制聚合物驱三层窜流油藏典型曲线理论图版;5)将步骤4)中获得的聚合物驱三层窜流油藏典型曲线理论图版与油田实测数据曲线进行拟合,得到聚合物驱三层窜流油藏的储层渗透率、表皮因子和井筒储集系数。
在所述步骤1)中,确定聚合物溶液粘度模型的修正公式,具体过程如下:
①建立考虑扩散和对流作用下的聚合物溶液浓度模型的修正公式:
式中,Cp0为聚合物溶液初始浓度,D为扩散系数,t为时间,v为渗流速度,r为径向距离,erf为高斯误差函数;
②建立考虑剪切作用下的聚合物溶液粘度模型的修正公式:
式中,μp为聚合物溶液的粘度,μw为水相的粘度,为零剪切速率下的聚合物溶液粘度,γ1/2为相应粘度对应的剪切速率,γ为粘度对应的剪切速率,Pa为聚合物溶液的非牛顿幂律指数,1.0<Pa<1.8,Cp为聚合物溶液浓度,Csep为水相中有效阳离子浓度,Sp为实验资料确定的聚合物溶液特性参数,A1、A2、A3为多项式的拟合系数,C'为与迂曲度有关的系数,K为渗透率,φ为孔隙度,Q为流量,k为剪切速率指数,h为地层厚度,r为径向距离;
③根据聚合物溶液流变性实验,得到A1=0.718,A2=0.289,A3=0.975,从而确定聚合物溶液粘度模型的修正公式:
在所述步骤2)中,确定聚合物驱三层窜流油藏数学模型,具体过程如下:
①确定各层单相微可压缩液体的不稳定渗流微分方程,其中,第一油层单相微可压缩液体的不稳定渗流微分方程为:
式中,p1、p2分别为第一、第二油层的压力,h1、h2、h3分别为第一、第二、第三油层的厚度,φ1为第一油层的孔隙度,Ct1为第一油层的综合压缩系数,λ12为第一、第二油层之间的窜流系数,rw为井半径;
第二油层单相微可压缩液体的不稳定渗流微分方程为:
式中,p3为第三油层的压力,φ2为第二油层的孔隙度,Ct2为第二油层的综合压缩系数,λ23为第二、第三油层之间的窜流系数;
第三油层单相微可压缩液体的不稳定渗流微分方程为:
式中,Ct3为第三油层的综合压缩系数,φ3为第三油层的孔隙度;
②确定聚合物驱三层窜流油藏数学模型的初始条件方程、内边界条件方程和外边界条件方程,其中,聚合物驱三层窜流油藏数学模型的初始条件方程:
p1(r,0)=p2(r,0)=p3(r,0)=p0j(j=1,2,3) (7)
p1(rw,t)=p2(rw,t)=p3(rw,t)=pw(t) (8)
式中,pw(t)为井底压力;p0j为第j层的原始地层压力;
聚合物驱三层窜流油藏数学模型的内边界条件方程:
式中,pwf为井底流压,B为体积系数,C为井筒储集系数,S1、S2、S3分别为第一、第二、第三油层的表皮系数,K1、K2、K3分别为第一、第二、第三油层的渗透率;
聚合物驱三层窜流油藏数学模型的外边界条件方程:
p1(∞,t)=p2(∞,t)=p3(∞,t)=p0 (11)
式中:p0为原始地层压力。
在所述步骤3)中,获得井底压力的数值解,具体过程如下:
①对空间和时间进行网格划分;
②对渗流方程、初始条件和边界条件进行差分离散化,其中,离散后的第一油层渗流扩散方程:
式中,i为对空间的离散,△x为空间网格大小,n为对时间的离散,△t为时间步长,R12为中间变量,T1中间变量,分别为n+1时刻时一层中第i+1、i和i-1个网格的压力;为n+1时刻时二层中第i个网格的压力;为n时刻时一层中第i个网格的压力;ri为i个网格的径向距离;分别为n+1时刻时第i和i-1个网格的聚合物溶液粘度;
离散后的第二油层渗流扩散方程:
式中,R22、R23分别为中间变量,T2中间变量;分别为n+1时刻时二层中第i+1和i-1个网格的压力;为n+1时刻时三层中第i个网格的压力;为n时刻时二层中第i个网格的压力;
离散后的第三油层渗流扩散方程:
式中,R33为中间变量,T3为中间变量,分别为n+1时刻时三层中第i+1和i-1个网格的压力;为n时刻时三层中第i个网格的压力;
离散后的第一油层内边界:
离散后的第二油层内边界:
离散后的第三油层内边界:
离散后的各油层外边界条件方程:
式中,为n+1时刻时一层第1个网格压力、为n+1时刻时一层初始压力、为n+1时刻时二层初始压力、为n+1时刻时二层第1个网格压力、为n+1时刻时三层初始压力、为n+1时刻时三层第1个网格压力、为n时刻时一层初始压力、为n时刻时一层第1个网格压力、为n时刻时二层初始压力、为n时刻时二层第1个网格压力、为n时刻时三层初始压力、为n时刻时三层第1个网格压力、pi为原始地层压力、为n时刻外边界压力、为n+1时刻时聚合物溶液初始粘度;
求得聚合物驱三层窜流油藏的压力分布公式:
式中,P(ri,tn)为井压力随井半径以及时间的变化,N为空间网格数,k为时间步,为不同时间步及不同空间网格下的压力;
对时间进行离散的计算公式:
式中,为n时刻的井底压力,单位为MPa;为n时刻初始压力、分别为n时刻一层压力;S为表皮系数;
依据公式(26)可以求出聚合物驱三层窜流油藏井底的压力变化,即:
式中,Pwf(tn)井底压力随时间的变化;
③对上述差分方程组进行数值迭代求解,求得聚合物驱三层窜流油藏井底压力的数值解。
在所述步骤4)中,井底压力无量纲化的计算公式:
式中,pwDj为无量纲井底压力;q为注聚井的注入量;
在所述步骤4)中,井筒储集系数无量纲的计算公式:
式中,CD为无量纲井筒储集系数;
在所述步骤4)中,地层系数无量纲的计算公式:
式中,χj为各油层的无量纲地层系数;
在所述步骤4)中,储容比无量纲的计算公式:
式中,ωj为各油层的无量纲储容比;
在所述步骤4)中,时间无量纲的计算公式:
式中,tD为时间的无量纲。
本发明由于采取以上技术方案,其具有以下优点:1、本发明对聚合物驱三层窜流油藏典型曲线理论图版与实测数据曲线进行拟合,得到聚合物驱三层窜流油藏试井解释的储层渗透率、表皮因子和井筒储集系数,并能够利用储层渗透率、表皮因子和井筒储集系数快速且准确地对实井油田实测数据进行解释,同时可以为具有类似储层条件的聚合物油藏试井解释、聚合物驱后储层评价提供依据,进而为具有类似条件的油藏开发提供了有力保障。2、本发明确立了聚合物驱粘度模型修正公式,并根据实验数据求得聚合物溶液参数,综合考虑了聚合物溶液与地层岩石之间的扩散和对流作用,以及聚合物溶液受到的剪切作用,能够真实的反应聚合物溶液在地层中的渗流机理。3、本发明确定了聚合物驱三层窜流油藏数学模型,能够快速地获得井底压力的数值解,绘制聚合物驱三层窜流油藏典型曲线理论图版,将聚合物驱三层窜流油藏典型曲线理论图版与油田实测数据曲线进行拟合,曲线拟合度好,进一步提高了油田实测数据解释结果的准确性,以便于指导油田开发生产。
附图说明
图1是本发明的流程结构示意图
图2是本发明等效模拟机械降解导致聚合物粘度损失的结构示意图
图3是本发明聚合物通过不同直径油管粘度损失实验结果的结构示意图
图4是本发明油藏简化物理模型的结构示意图
图5是本发明油藏试井解释典型曲线图版的结构示意图
图6是本发明不同窜流系数影响下油藏试井解释典型曲线图版的结构示意图
图7是本发明不同地层系数比影响下油藏试井解释典型曲线图版的结构示意图
图8是本发明不同储容比影响下油藏试井解释典型曲线图版的结构示意图
图9是本发明不同初始聚合物浓度影响下油藏试井解释典型曲线图版的结构示意图
图10是本发明油藏试井解释理论图版与实井测试数据拟合曲线的结构示意图
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示,本发明提出的聚合物驱三层窜流油藏试井解释参数的获取方法,包括以下步骤:
1)根据聚合物溶液的流变特性,确定聚合物溶液粘度模型的修正公式,并根据聚合物驱三层窜流油藏特征,建立聚合物驱三层窜流油藏物理模型。
在步骤1)中,确定聚合物溶液粘度模型的修正公式的过程如下:
①建立考虑扩散和对流作用下的聚合物溶液浓度模型的修正公式:
式中,Cp0为聚合物溶液初始浓度,D为扩散系数,t为时间,v为渗流速度,r为径向距离,erf为高斯误差函数。
②建立考虑剪切作用下的聚合物溶液粘度模型的修正公式:
式中,μp为聚合物溶液的粘度,μw为水相的粘度,为零剪切速率下的聚合物溶液粘度,γ1/2为相应粘度对应的剪切速率,γ为粘度对应的剪切速率,Pa为聚合物溶液的非牛顿幂律指数,1.0<Pa<1.8,Cp为聚合物溶液浓度,Csep为水相中有效阳离子浓度,Sp为实验资料确定的聚合物溶液特性参数,A1、A2、A3为多项式的拟合系数,C'为与迂曲度有关的系数,K为渗透率,φ为孔隙度,Q为流量,k为剪切速率指数,h为地层厚度,r为径向距离。
③根据聚合物溶液流变性实验,可得到A1=0.718,A2=0.289,A3=0.975,从而确定聚合物溶液粘度模型的修正公式:
在步骤1)中,建立聚合物驱三层窜流油藏物理模型需要符合以下假设条件:第一、第二、第三油层中均含有聚合物和水两种组分,且聚合物与水完全混溶。各油层聚合物溶液性质相同,各油层的物性不同。每一油层的流动均遵循达西渗流,并忽略重力的影响。每一油层中的流体微可压缩,每一油层的流动过程均为等温渗流,各层间均为拟稳态窜流。各油层井储,表皮恒定。
2)根据步骤1)中获得的聚合物溶液粘度模型的修正公式和聚合物驱三层窜流油藏物理模型,确定聚合物驱三层窜流油藏数学模型;
在步骤2)中,确定聚合物驱三层窜流油藏数学模型,具体过程如下:
①确定各层单相微可压缩液体的不稳定渗流微分方程,其中,第一油层单相微可压缩液体的不稳定渗流微分方程为:
式中,p1、p2分别为第一、第二油层的压力,h1、h2、h3分别为第一、第二、第三油层的厚度,φ1为第一油层的孔隙度,Ct1为第一油层的综合压缩系数,λ12为第一、第二油层之间的窜流系数,rw为井半径。
第二油层单相微可压缩液体的不稳定渗流微分方程为:
式中,p3为第三油层的压力,φ2为第二油层的孔隙度,Ct2为第二油层的综合压缩系数,λ23为第二、第三油层之间的窜流系数。
第三油层单相微可压缩液体的不稳定渗流微分方程为:
式中,Ct3为第三油层的综合压缩系数,φ3为第三油层的孔隙度。
②确定聚合物驱三层窜流油藏数学模型的初始条件方程、内边界条件方程和外边界条件方程,其中,聚合物驱三层窜流油藏数学模型的初始条件方程:
p1(r,0)=p2(r,0)=p3(r,0)=p0j(j=1,2,3) (7)
p1(rw,t)=p2(rw,t)=p3(rw,t)=pw(t) (8)
式中,pw(t)为井底压力,p0j为第j层的原始地层压力。
聚合物驱三层窜流油藏数学模型的内边界条件方程:
式中,pwf为井底流压,B为体积系数,C为井筒储集系数,S1、S2、S3分别为第一、第二、第三油层的表皮系数,K1、K2、K3分别为第一、第二、第三油层的渗透率。
聚合物驱三层窜流油藏数学模型的外边界条件方程:
p1(∞,t)=p2(∞,t)=p3(∞,t)=p0 (11)
式中:p0为原始地层压力。
3)采用有限差分方法求解步骤2)中获得的聚合物驱三层窜流油藏数学模型,得到井底压力的数值解;
在步骤3)中,获得井底压力的数值解,具体过程如下:
①对空间和时间进行网格划分;
②对渗流方程、初始条件和边界条件进行差分离散化,其中,离散后的第一油层渗流扩散方程:
式中,i为对空间的离散,△x为空间网格大小,n为对时间的离散,△t为时间步长,R12为中间变量,T1中间变量,分别为n+1时刻时一层中第i+1、i和i-1个网格的压力;为n+1时刻时二层中第i个网格的压力;为n时刻时一层中第i个网格的压力;ri为i个网格的径向距离;分别为n+1时刻时第i和i-1个网格的聚合物溶液粘度。
离散后的第二油层渗流扩散方程:
式中,R22、R23分别为中间变量,T2中间变量;分别为n+1时刻时二层中第i+1和i-1个网格的压力;为n+1时刻时三层中第i个网格的压力;为n时刻时二层中第i个网格的压力。
离散后的第三油层渗流扩散方程:
式中,R33为中间变量,T3为中间变量,分别为n+1时刻时三层中第i+1和i-1个网格的压力;为n时刻时三层中第i个网格的压力。
离散后的第一油层内边界:
离散后的第二油层内边界:
离散后的第三油层内边界:
离散后的各油层外边界条件方程:
式中,为n+1时刻时一层第1个网格压力、为n+1时刻时一层初始压力、为n+1时刻时二层初始压力、为n+1时刻时二层第1个网格压力、为n+1时刻时三层初始压力、为n+1时刻时三层第1个网格压力、为n时刻时一层初始压力、为n时刻时一层第1个网格压力、为n时刻时二层初始压力、为n时刻时二层第1个网格压力、为n时刻时三层初始压力、为n时刻时三层第1个网格压力、pi为原始地层压力、为n时刻外边界压力、为n+1时刻时聚合物溶液初始粘度。
联立公式(11)~(24),编制计算机程序,求得聚合物驱三层窜流油藏的压力分布公式:
式中,P(ri,tn)为井压力随井半径以及时间的变化,N为空间网格数,k为时间步,为不同时间步及不同空间网格下的压力。
将公式(9)离散,变形整理可得:
式中,为n时刻的井底压力,单位为MPa;为n时刻初始压力、分别为n时刻一层压力;S为表皮系数。
依据公式(26)可以求出聚合物驱三层窜流油藏井底的压力变化,即:
式中,Pwf(tn)为井底压力随时间的变化。
③对上述差分方程组进行数值迭代求解,求得聚合物驱三层窜流油藏井底压力的数值解。
4)对井筒储集系数、地层系数、储容比、时间和步骤3)中获得的井底压力分别进行无量纲化,并绘制聚合物驱三层窜流油藏典型曲线理论图版;
在步骤4)中,井底压力无量纲化的计算公式:
式中,pwDj为无量纲井底压力;q为注聚井的注入量。
在步骤4)中,井筒储集系数无量纲的计算公式:
式中,CD为无量纲井筒储集系数。
在步骤4)中,地层系数无量纲的计算公式:
式中,χj为各油层的无量纲地层系数。
在步骤4)中,储容比无量纲的计算公式:
式中,ωj为各油层的无量纲储容比。
在步骤4)中,时间无量纲的计算公式:
式中,tD为时间的无量纲。
5)将步骤4)中获得的聚合物驱三层窜流油藏典型曲线理论图版与油田实测数据曲线进行拟合,得到聚合物驱三层窜流油藏的储层渗透率、表皮因子和井筒储集系数等三个试井解释参数,接着利用聚合物驱三层窜流油藏试井解释参数,分析聚合物驱效果和储层状况,并采取相应措施指导油田开发生产;
在步骤5)中,首先输入基础数据,包括每一油层的厚度、孔隙度、渗透率、综合压缩系数、表皮系数、井筒储集系数、注聚井注入量、体积系数、水相粘度、聚合物初始浓度、扩散系数、原始地层压力、井径、储容比,以及相邻两油层间的窜流系数,接着通过调整输入的渗透率,表皮系数,井筒储集系数,窜流系数,地层系数,储容比,计算得到理论压力和理论压力导数曲线;然后利用理论压力曲线、理论压力导数曲线和油田压力的实测数据,拟合理论压力曲线和实井压力曲线,以及理论压力导数曲线和实井压力导数曲线,最后根据拟合结果得到聚合物驱三层窜流油藏的渗透率、表皮系数和井筒储集系数,接着利用储层渗透率、表皮因子和井筒储集系数分析聚合物驱效果以及储层状况,及时采取相应措施指导油田开发生产,提高油田采收率。
下面列举两具体实施例:
实施例一
本实施例使用的部分水解聚丙烯酰胺是由中国海油提供,聚丙烯酰胺水解度为27%、分子量为4630,模拟地层水结合氯化镁、氯化钙和硫酸钠,制成离子浓度为43000毫克/升的盐水。聚合物溶液在井筒和射孔段的流速分别为0.0368米/秒、0.3316米/秒,并将聚合物溶液机械降解后静置一天。为了消除机械降解对聚合物溶液流变性的影响,采用不同半径的管流模拟聚合物溶液从井筒到射孔段的流动(如图2所示),利用转速为60000转/分的HAAKE流变仪(哈克流变仪)测量不同浓度聚合物在0.01s-1的剪切速率下的粘度,结果表明机械降解大约使粘度降低20%(如图3所示)。
1)根据聚合物溶液的流变特性,确定聚合物溶液粘度模型的修正公式,并根据聚合物驱三层窜流油藏特征,建立聚合物驱三层窜流油藏物理模型(如图4所示),求得聚合物溶液参数A1=0.718,A2=0.289,A3=0.975。2)根据步骤1)中获得的聚合物溶液粘度模型的修正公式和聚合物驱三层窜流油藏物理模型,确定聚合物驱三层窜流油藏数学模型。3)采用有限差分方法求解步骤2)中获得的聚合物驱三层窜流油藏数学模型,得到井底压力的数值解。4)对井筒储集系数、地层系数、储容比、时间和步骤3)中获得的井底压力分别进行无量纲化,并绘制聚合物驱三层窜流油藏典型曲线理论图版(如图5所示);可以看出,存在窜流的三层油藏聚合物驱典型曲线可以划分为5个流动阶段:第Ⅰ段是纯井筒储集阶段,压力和压力导数重合,反映的是井储阶段的压力响应特征;第Ⅱ段是过渡段,描述纯井筒储集阶段到内区径向流阶段的压力响应特征;第Ⅲ段是第二油层向第一油层的窜流作用段;第Ⅳ段是第三油层向第二油层的窜流作用段;第Ⅴ段是总系统达到径向流作用阶段,由于受到聚合物溶液牛顿流体性质的影响,第Ⅴ段的曲线表现出小幅度上翘。绘制不同窜流系数影响下的典型曲线图版(如图6所示),不同地层系数比影响下的典型曲线图版(如图7所示),不同储容比影响下的典型曲线图版(如图8所示),不同初始聚合物浓度影响下的典型曲线图版(如图9所示)。5)将步骤4)中获得的聚合物驱三层窜流油藏典型曲线理论图版与油田实测数据曲线进行拟合,得到聚合物驱三层窜流油藏试井解释的储层渗透率、表皮因子和井筒储集系数,并利用储层渗透率、表皮因子和井筒储集系数对油田实测数据进行解释。
实施例二
本实施例对渤海油田聚合物驱井进行了试井解释,本实施例中油田实测数据与理论图版拟合曲线(如图10所示),得到储层渗透率、表皮因子和井筒储集系数的解释数据(如表1所示)。表1中聚驱前的数据是该油田井水驱时的测试数据,聚驱后的数据是对该油田井进行试井解释的相关数据,从解释结果可以看出第一油层、第二油层变化渗透率表皮系数变化较小,第三油层渗透率下降较多,表皮系数增加较多,分析该油层发生了聚合物堵塞,于是对该层采取酸化措施,措施后对应油井产量较措施前提高11.2%,可以得出聚合物驱三层窜流油藏试井分析方法能够较为准确地对油田实测数据进行解释,指导油田开发生产。
表1聚合物驱前后储层参数对比
油藏参数 | 聚驱前 | 聚驱后 |
第一油层渗透率(毫达西) | 977 | 940 |
第二油层渗透率(毫达西) | 166 | 145 |
第三油层渗透率(毫达西) | 148 | 46 |
第一油层表皮系数 | 1.01 | 1.11 |
第二油层表皮系数 | 1.08 | 1.77 |
第三油层表皮系数 | 0.84 | 2.95 |
井筒储集系数(立方米/兆帕) | 4.7 | 4.71 |
上述各实施例仅用于说明本发明,其中各部件的结构、连接方式等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (5)
1.一种聚合物驱三层窜流油藏试井解释参数的获取方法,包括以下步骤:
1)根据聚合物溶液的流变特性,确定聚合物溶液粘度模型的修正公式,并根据聚合物驱三层窜流油藏特征,建立聚合物驱三层窜流油藏物理模型;
2)根据步骤1)中获得的聚合物溶液粘度模型的修正公式和聚合物驱三层窜流油藏物理模型,确定聚合物驱三层窜流油藏数学模型;
3)采用有限差分方法求解步骤2)中获得的聚合物驱三层窜流油藏数学模型,得到井底压力的数值解;
4)对井筒储集系数、地层系数、储容比、时间和步骤3)中获得的井底压力分别进行无量纲化,并绘制聚合物驱三层窜流油藏典型曲线理论图版;
5)将步骤4)中获得的聚合物驱三层窜流油藏典型曲线理论图版与油田实测数据曲线进行拟合,得到聚合物驱三层窜流油藏的储层渗透率、表皮因子和井筒储集系数。
2.如权利要求1所述的一种聚合物驱三层窜流油藏试井解释参数的获取方法,其特征在于:在所述步骤1)中,确定聚合物溶液粘度模型的修正公式,具体过程如下:
①建立考虑扩散和对流作用下的聚合物溶液浓度模型的修正公式:
式中,Cp0为聚合物溶液初始浓度,D为扩散系数,t为时间,v为渗流速度,r为径向距离,erf为高斯误差函数;
②建立考虑剪切作用下的聚合物溶液粘度模型的修正公式:
式中,μp为聚合物溶液的粘度,μw为水相的粘度,为零剪切速率下的聚合物溶液粘度,γ1/2为相应粘度对应的剪切速率,γ为粘度对应的剪切速率,Pa为聚合物溶液的非牛顿幂律指数,1.0<Pa<1.8,Cp为聚合物溶液浓度,Csep为水相中有效阳离子浓度,Sp为实验资料确定的聚合物溶液特性参数,A1、A2、A3为多项式的拟合系数,C'为与迂曲度有关的系数,K为渗透率,φ为孔隙度,Q为流量,k为剪切速率指数,h为地层厚度,r为径向距离;
③根据聚合物溶液流变性实验,得到A1=0.718,A2=0.289,A3=0.975,从而确定聚合物溶液粘度模型的修正公式:
3.如权利要求2所述的一种聚合物驱三层窜流油藏试井解释参数的获取方法,其特征在于:在所述步骤2)中,确定聚合物驱三层窜流油藏数学模型,具体过程如下:
①确定各层单相微可压缩液体的不稳定渗流微分方程,其中,第一油层单相微可压缩液体的不稳定渗流微分方程为:
式中,p1、p2分别为第一、第二油层的压力,h1、h2、h3分别为第一、第二、第三油层的厚度,φ1为第一油层的孔隙度,Ct1为第一油层的综合压缩系数,λ12为第一、第二油层之间的窜流系数,rw为井半径;
第二油层单相微可压缩液体的不稳定渗流微分方程为:
式中,p3为第三油层的压力,φ2为第二油层的孔隙度,Ct2为第二油层的综合压缩系数,λ23为第二、第三油层之间的窜流系数;
第三油层单相微可压缩液体的不稳定渗流微分方程为:
式中,Ct3为第三油层的综合压缩系数,φ3为第三油层的孔隙度;
②确定聚合物驱三层窜流油藏数学模型的初始条件方程、内边界条件方程和外边界条件方程,其中,聚合物驱三层窜流油藏数学模型的初始条件方程:
p1(r,0)=p2(r,0)=p3(r,0)=p0j j为1,2,3 (7)
p1(rw,t)=p2(rw,t)=p3(rw,t)=pw(t) (8)
式中,pw(t)为井底压力;p0j为第j层的原始地层压力;
聚合物驱三层窜流油藏数学模型的内边界条件方程:
式中,pwf为井底流压,B为体积系数,C为井筒储集系数,S1、S2、S3分别为第一、第二、第三油层的表皮系数,K1、K2、K3分别为第一、第二、第三油层的渗透率;
聚合物驱三层窜流油藏数学模型的外边界条件方程:
p1(∞,t)=p2(∞,t)=p3(∞,t)=p0 (11)
式中:p0为原始地层压力。
4.如权利要求3所述的一种聚合物驱三层窜流油藏试井解释参数的获取方法,其特征在于:在所述步骤3)中,获得井底压力的数值解,具体过程如下:
①对空间和时间进行网格划分;
②对渗流方程、初始条件和边界条件进行差分离散化,其中,离散后的第一油层渗流扩散方程:
式中,i为对空间的离散,Δx为空间网格大小,n为对时间的离散,Δt为时间步长,R12为中间变量,T1中间变量,分别为n+1时刻时一层中第i+1、i和i-1个网格的压力;为n+1时刻时二层中第i个网格的压力;为n时刻时一层中第i个网格的压力;ri为i个网格的径向距离;分别为n+1时刻时第i和i-1个网格的聚合物溶液粘度;
离散后的第二油层渗流扩散方程:
式中,R22、R23分别为中间变量,T2中间变量;分别为n+1时刻时二层中第i+1和i-1个网格的压力;为n+1时刻时三层中第i个网格的压力;为n时刻时二层中第i个网格的压力;
离散后的第三油层渗流扩散方程:
式中,R33为中间变量,T3为中间变量,分别为n+1时刻时三层中第i+1和i-1个网格的压力;为n时刻时三层中第i个网格的压力;
离散后的第一油层内边界:
离散后的第二油层内边界:
离散后的第三油层内边界:
离散后的各油层外边界条件方程:
式中,为n+1时刻时一层第1个网格压力、为n+1时刻时一层初始压力、为n+1时刻时二层初始压力、为n+1时刻时二层第1个网格压力、为n+1时刻时三层初始压力、为n+1时刻时三层第1个网格压力、为n时刻时一层初始压力、为n时刻时一层第1个网格压力、为n时刻时二层初始压力、为n时刻时二层第1个网格压力、为n时刻时三层初始压力、为n时刻时三层第1个网格压力、pi为原始地层压力、为n时刻外边界压力、为n+1时刻时聚合物溶液初始粘度;
求得聚合物驱三层窜流油藏的压力分布公式:
式中,P(ri,tn)为井压力随井半径以及时间的变化,N为空间网格数,k为时间步,为不同时间步及不同空间网格下的压力;
对时间进行离散的计算公式:
式中,为n时刻的井底压力,单位为MPa;为n时刻初始压力、分别为n时刻一层压力;S为表皮系数;
依据公式(26)可以求出聚合物驱三层窜流油藏井底的压力变化,即:
式中,Pwf(tn)井底压力随时间的变化;
③对上述差分方程组进行数值迭代求解,求得聚合物驱三层窜流油藏井底压力的数值解。
5.如权利要求4所述的一种聚合物驱三层窜流油藏试井解释参数的获取方法,其特征在于:在所述步骤4)中,井底压力无量纲化的计算公式:
式中,pwDj为无量纲井底压力;q为注聚井的注入量;
在所述步骤4)中,井筒储集系数无量纲的计算公式:
式中,CD为无量纲井筒储集系数;
在所述步骤4)中,地层系数无量纲的计算公式:
式中,χj为各油层的无量纲地层系数;
在所述步骤4)中,储容比无量纲的计算公式:
式中,ωj为各油层的无量纲储容比;
在所述步骤4)中,时间无量纲的计算公式:
式中,tD为时间的无量纲。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710072135.2A CN106761621B (zh) | 2017-02-08 | 2017-02-08 | 一种聚合物驱三层窜流油藏试井解释参数的获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710072135.2A CN106761621B (zh) | 2017-02-08 | 2017-02-08 | 一种聚合物驱三层窜流油藏试井解释参数的获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106761621A CN106761621A (zh) | 2017-05-31 |
CN106761621B true CN106761621B (zh) | 2019-08-02 |
Family
ID=58955748
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710072135.2A Active CN106761621B (zh) | 2017-02-08 | 2017-02-08 | 一种聚合物驱三层窜流油藏试井解释参数的获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106761621B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107526891B (zh) * | 2017-08-24 | 2020-11-10 | 中海石油(中国)有限公司 | 一种聚合物驱大孔道油藏试井分析方法 |
CN110322362B (zh) * | 2018-03-29 | 2021-07-20 | 中国石油化工股份有限公司 | 一种聚合物驱早期窜流识别方法及装置 |
CN110175412B (zh) * | 2019-05-28 | 2023-03-28 | 广东海洋大学 | 油藏数据拟合方法、油藏数据分析方法及装置 |
CN112685971B (zh) * | 2020-12-26 | 2022-12-02 | 合肥工业大学 | 亲水油藏的参数解释方法及装置 |
CN113027434B (zh) * | 2021-03-25 | 2022-07-19 | 东北石油大学 | 一种确定射孔完井近井地带含聚类液体流变性参数的方法 |
CN114372398B (zh) * | 2022-01-12 | 2022-10-11 | 东北石油大学 | 一种考虑储层裂缝闭合效应的聚合物驱试井解释方法 |
CN116562179B (zh) * | 2023-03-16 | 2024-01-09 | 东北石油大学 | 一种考虑流变特性的聚合物驱压力动态特征半解析方法 |
CN116341393B (zh) * | 2023-05-26 | 2023-08-04 | 中国石油大学(华东) | 一种自动非稳态试井解释方法、装置、设备及介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4206809A (en) * | 1978-11-02 | 1980-06-10 | Marathon Oil Company | Slug design in micellar-polymer flooding |
CN104533363A (zh) * | 2015-01-23 | 2015-04-22 | 中国石油大学(华东) | 聚合物驱注入采出井动态关联系数的确定方法 |
CN106150461A (zh) * | 2015-02-11 | 2016-11-23 | 中国石油化工股份有限公司 | 化学驱油藏数值模拟方法 |
-
2017
- 2017-02-08 CN CN201710072135.2A patent/CN106761621B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4206809A (en) * | 1978-11-02 | 1980-06-10 | Marathon Oil Company | Slug design in micellar-polymer flooding |
CN104533363A (zh) * | 2015-01-23 | 2015-04-22 | 中国石油大学(华东) | 聚合物驱注入采出井动态关联系数的确定方法 |
CN106150461A (zh) * | 2015-02-11 | 2016-11-23 | 中国石油化工股份有限公司 | 化学驱油藏数值模拟方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106761621A (zh) | 2017-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106761621B (zh) | 一种聚合物驱三层窜流油藏试井解释参数的获取方法 | |
CN102953726B (zh) | 一种水驱油田优势通道识别方法及装置 | |
CN106522928B (zh) | 一种酸化压裂后停泵测井口压降不稳定试井方法 | |
CN111236908A (zh) | 一种适用在低渗透致密气藏中的多段压裂水平井产能预测模型及产能敏感性分析的方法 | |
CN103226641B (zh) | 深水气液两相流循环温度压力耦合计算方法 | |
CN110593865B (zh) | 一种油藏缝洞特征参数试井解释方法 | |
CN107526891B (zh) | 一种聚合物驱大孔道油藏试井分析方法 | |
CN109001438A (zh) | 一种封缝堵气实验模拟装置及测试方法 | |
Izgec et al. | Real-time performance analysis of water-injection wells | |
CN105317407B (zh) | 一种特高含水期表外储层的开发方法 | |
CN107066679A (zh) | 一种用于聚合物驱双层窜流油藏试井分析系统及方法 | |
CN108894777A (zh) | 一种分层压裂多层合采油气藏储层物性及裂缝特性参数的确定方法 | |
CN108442923B (zh) | 致密油气藏水平井椭圆形泄流区域长、短半轴的确定方法 | |
Liu et al. | Field application of polymer microspheres flooding: a pilot test in offshore heavy oil reservoir | |
CN105443093A (zh) | 用于注聚井的井口组合测试装置及其方法 | |
CN106600439A (zh) | 一种低流度油藏物理约束试井解释方法及装置 | |
El Sgher et al. | The Stress Shadowing Impact on the Production Performance of Marcellus Shale | |
Wang et al. | Case studies: Pressure-transient analysis for water injector with the influence of waterflood-induced fractures in tight reservoir | |
Olsen et al. | Waterflood performance evaluation in a chalk reservoir with an ensemble of tools | |
CN108468540B (zh) | 一种低渗-特低渗油气藏单井泄流面积的计算方法 | |
CN106844975B (zh) | 一种确定注气井注入中早期co2等效波及半径的方法及装置 | |
He et al. | Waterflood direction and front characterization with multiple methods: a case study in Changqing oilfield, China | |
Fan et al. | Evaluation of sweep efficiency in flooding process of reservoir development using substitution index | |
CN112943160B (zh) | 一种油井化学堵水效果的评价方法 | |
Fisch et al. | Hydraulic testing and reservoir characterization of the Taufkirchen site in the Bavarian Molasse Basin, Germany |
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 |