CN105888631B - 一种裂缝性稠油油藏热水驱物理模拟相似准数的确定方法 - Google Patents
一种裂缝性稠油油藏热水驱物理模拟相似准数的确定方法 Download PDFInfo
- Publication number
- CN105888631B CN105888631B CN201610305492.4A CN201610305492A CN105888631B CN 105888631 B CN105888631 B CN 105888631B CN 201610305492 A CN201610305492 A CN 201610305492A CN 105888631 B CN105888631 B CN 105888631B
- Authority
- CN
- China
- Prior art keywords
- crack
- hot water
- water flooding
- mathematical model
- heavy crude
- 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
- E21B43/20—Displacing by water
-
- 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
- E21B43/24—Enhanced recovery methods for obtaining hydrocarbons using heat, e.g. steam injection
Abstract
本发明涉及一种裂缝性稠油油藏热水驱物理模拟相似准数的确定方法,包括以下步骤:1)根据裂缝稠油油藏热水驱渗流特点,建立基本假设条件;2)根据裂缝中的物质平衡方程、裂缝中的能量平衡方程,并结合压力、温度初始条件以及内外边界条件,获得裂缝稠油油藏热水驱数学模型;3)将裂缝稠油油藏热水驱数学模型中的相关变量无量纲化,并化简得到裂缝稠油油藏热水驱无量纲数学模型;4)列出裂缝稠油油藏热水驱数学模型、裂缝稠油油藏热水驱无量纲数学模型中的物理量;5)列出步骤4)中全部物理量所包括的基本量纲;6)根据裂缝稠油油藏热水驱无量纲数学模型,结合π定理及相似理论分析,获取相似准数。
Description
技术领域
本发明涉及一种裂缝性稠油油藏热水驱物理模拟相似准数的确定方法。
背景技术
裂缝性稠油油藏注冷水开发效果较差,而注热水则可以降低原油粘度,提高水驱效果。在针对裂缝性稠油油藏热水驱进行物理模拟研究中,首先需要建立一套相似准则,获取相似准数,以指导物理模型与实际油藏满足相似性,但此前尚未发现相关研究报道。
发明内容
针对上述问题,本发明的目的是提供一种裂缝性稠油油藏热水驱物理模拟相似准数的确定方法。
为实现上述目的,本发明采用以下技术方案:一种裂缝性稠油油藏热水驱物理模拟相似准数的确定方法,包括以下步骤:1)根据裂缝稠油油藏热水驱渗流特点,建立基本假设条件;2)根据裂缝中的物质平衡方程、裂缝中的能量平衡方程,并结合压力、温度初始条件以及内外边界条件,获得裂缝稠油油藏热水驱数学模型;3)将裂缝稠油油藏热水驱数学模型中的相关变量无量纲化,并化简得到裂缝稠油油藏热水驱无量纲数学模型;4)列出裂缝稠油油藏热水驱数学模型、裂缝稠油油藏热水驱无量纲数学模型中的物理量;5)列出步骤4)中全部物理量所包括的基本量纲;6)根据裂缝稠油油藏热水驱无量纲数学模型,结合π定理及相似理论分析,获取相似准数。
所述步骤1)所建立的基本假设条件如下:①储层介质包括基质、裂缝,考虑双孔单渗;②流体包括油、水两相;③油藏中任一微元,达到热平衡与相平衡;④忽略裂缝毛管力,考虑基质渗吸,考虑重力;⑤忽略岩石、流体压缩性和热膨胀,忽略烃类热裂解。
所述裂缝中的物质平衡方程如下:
所述裂缝中的能量平衡方程如下:
所述裂缝稠油油藏热水驱数学模型如下:
式中,x、y、z分别表示空间三个方向的坐标;Kx、Ky、Kz分别表示裂缝x、y、z方向渗透率,mD;Sw表示裂缝中水饱和度,无量纲;φ表示裂缝孔隙度,无量纲;μw、μo分别表示地层水、地层油粘度,mPa·s;p+γwz、p+γoz分别表示裂缝中水、油位势,Mpa;γw、γo分别表示水、油重度,106×N/m3;p表示压力,Mpa;R表示单位体积基质岩块可采储量,m3;t*表示基质岩块渗吸半周期,d;λR表示储层热导系数,KJ/m/d/℃;T表示储层温度,℃;ρwCw、ρoCo、ρRCR分别表示水、油、岩石体积比热,KJ/m3/℃;ρw、ρo、ρR分别表示水、油密度Kg/m3;Cw、Co、CR分别表示水、油、岩石质量比热,KJ/Kg/℃;pi表示油藏初始压力,为常数,MPa;Ti、Tinj、T0分别表示油藏初始温度、水井注入水温度和初始注入水温度,其中Ti、T0为常数,℃;lb、ub、n分别表示侧边界、顶底边界和垂直于边界方向;pinj、ppro、△p分别表示水井压力、油井压力和注采压差,△p为常数,MPa;Lx、Ly、Lz分别表示室内模型或油藏模拟区域x、y、z方向的特征长度,m;t'=Lxφ/[(△p/Lx)·(Kx/μw)]表示室内模型或油藏模拟区域的特征时间,d;t表示时间,d;τ是求时间积分时引入符号,表示时间微元;e为常数,约等于2.718281828。
所述步骤3)中将裂缝稠油油藏热水驱数学模型中的相关变量无量纲化的过程如下:
所述步骤4)中的列出的物理量如下:包括4个自变量x、y、z、t,4个因变量p+γwz、p+γoz-pi、Sw、T-Ti,20个参量Kx、Ky、Kz、μw、μo、γw-γo、R、t*、λR、ρwCw、ρoCo、ρRCR、△p、T0-Ti、φ、rw、Lx、Ly、Lz、t'。
所述步骤5)中列出的基本量纲如下:质量M、长度L、温度T、时间t。
所述步骤6)中获取的相似准数为20个,具体如下:
其中,π1~π16为室内模拟常规注水开发所需准数,π17~π20为室内模拟注热水开发时考虑温度及能量变化所附加的准数。
本发明由于采取以上技术方案,其具有以下优点:1、本发明提出了一套技术方法,根据本方法可以获得裂缝性稠油油藏热水驱物理模拟相似准数,可指导物理模型的建立,使建立相似物理模型模拟裂缝性稠油油藏热水驱成为可能。2、本发明给出了定量化、可操作的技术方法和实施步骤。3、本发明不仅适用于油田开发研究领域,还可以供其它与裂缝性介质中两相渗流现象有关的研究领域使用和参考,例如地下污染物运移研究等。
具体实施方式
下面结合实施例对本发明进行详细的描述。
本发明一种裂缝性稠油油藏热水驱物理模拟相似参数的确定方法,包括以下步骤:
1)根据裂缝稠油油藏热水驱渗流特点,建立下述五个基本假设条件:①储层介质包括基质、裂缝,考虑双孔单渗;②流体包括油、水两相;③油藏中任一微元,达到热平衡与相平衡;④忽略裂缝毛管力,考虑基质渗吸,考虑重力;⑤忽略岩石、流体压缩性和热膨胀,忽略烃类热裂解。
2)根据裂缝中的物质平衡方程(1)式、裂缝中的能量平衡方程(2)式,并结合压力、温度初始条件以及内外边界条件,获得裂缝稠油油藏热水驱数学模型(3)式:
上述式中,式中,x、y、z分别表示空间三个方向的坐标;Kx、Ky、Kz分别表示裂缝x、y、z方向渗透率,mD;Sw表示裂缝中水饱和度,无量纲;φ表示裂缝孔隙度,无量纲;μw、μo分别表示地层水、地层油粘度,mPa·s;p+γwz、p+γoz分别表示裂缝中水、油位势,Mpa;γw、γo分别表示水、油重度,106×N/m3;p表示压力,Mpa;R表示单位体积基质岩块可采储量,m3;t*表示基质岩块渗吸半周期,d;λR表示储层热导系数,KJ/m/d/℃;T表示储层温度,℃;ρwCw、ρoCo、ρRCR分别表示水、油、岩石体积比热,KJ/m3/℃;ρw、ρo、ρR分别表示水、油密度Kg/m3;Cw、Co、CR分别表示水、油、岩石质量比热,KJ/Kg/℃;pi表示油藏初始压力,为常数,MPa;Ti、Tinj、T0分别表示油藏初始温度、水井注入水温度和初始注入水温度,其中Ti、T0为常数,℃;lb、ub、n分别表示侧边界、顶底边界和垂直于边界方向;pinj、ppro、△p分别表示水井压力、油井压力和注采压差,△p为常数,MPa;Lx、Ly、Lz分别表示室内模型或油藏模拟区域x、y、z方向的特征长度,m;t'=Lxφ/[(△p/Lx)·(Kx/μw)]表示室内模型或油藏模拟区域的特征时间,d;t表示时间,d;τ是求时间积分时引入符号,表示时间微元;e为常数,约等于2.718281828。
3)将(3)式中的相关变量无量纲化,得到(4)式,再将(4)式代回到(3)式,化简后得到裂缝稠油油藏热水驱无量纲数学模型(5)式:
4)列出热水驱数学模型(3)~(5)式中的物理量,包括4个自变量x、y、z、t,4个因变量p+γwz、p+γoz-pi、Sw、T-Ti,20个参量Kx、Ky、Kz、μw、μo、γw-γo、R、t*、λR、ρwCw、ρoCo、ρRCR、△p、T0-Ti、φ、rw、Lx、Ly、Lz、t'。
5)列出热水驱数学模型中的基本量纲,即4基本量纲:质量M、长度L、温度T、时间t。
6)根据π定理,使用28个物理量描述包含4个基本量纲的裂缝稠油油藏热水驱渗流过程,可等价于用28-4=24个无量纲参数组和描述该渗流过程。为便于物理模拟,这些无量纲参数组合应尽可能少,因此需对无量纲参数组合数进行分析并精简:28个物理量中,Kx、Ky、Kz、μw、μo总是可以推导成Kx/μo、Ky/μo、Kz/μo、μo/μw,φ、t总是可以推导成φ/t,λR、ρwCw、ρoCo、ρRCR总是可以推导成λR/ρRCR、ρwCw/ρRCR、ρoCo/ρRCR,t'=Lxφ/[(△p/Lx)·(Kx/μw)]可视为独立参量。因此,独立物理量共24个,根据π定理,该物理过程等价于使用20个无量纲参数组合来描述,根据裂缝稠油油藏热水驱无量纲数学模型(5)式以及相似理论分析,得到20个无量纲参数组合(即相似准则数组),如表1所示。表1中20个相似准数均具有明确的物理意义,且π1~π16为室内模拟常规注水开发所需准数,π17~π20为室内模拟注热水开发时考虑温度及能量变化所附加的准数。
表1裂缝稠油油藏注热水开发室内比例模拟相似准数组
裂缝性稠油油藏热水驱物理模拟的相似准数确定后,可以以上述相似准数为依据设计物理模拟过程。
上述各实施例仅用于说明本发明,其中方法的实施步骤等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (5)
1.一种裂缝性稠油油藏热水驱物理模拟相似准数的确定方法,包括以下步骤:1)根据裂缝稠油油藏热水驱渗流特点,建立基本假设条件;2)根据裂缝中的物质平衡方程、裂缝中的能量平衡方程,并结合压力、温度初始条件以及内外边界条件,获得裂缝稠油油藏热水驱数学模型;3)将裂缝稠油油藏热水驱数学模型中的相关变量无量纲化,并化简得到裂缝稠油油藏热水驱无量纲数学模型;4)列出裂缝稠油油藏热水驱数学模型、裂缝稠油油藏热水驱无量纲数学模型中的物理量;5)列出步骤4)中全部物理量所包括的基本量纲;6)根据裂缝稠油油藏热水驱无量纲数学模型,结合π定理及相似理论分析,获取相似准数;
所述裂缝中的物质平衡方程如下:
所述裂缝中的能量平衡方程如下:
所述裂缝稠油油藏热水驱数学模型如下:
式中,x、y、z分别表示空间三个方向的坐标;Kx、Ky、Kz分别表示裂缝x、y、z方向渗透率,mD;Sw表示裂缝中水饱和度,无量纲;φ表示裂缝孔隙度,无量纲;μw、μo分别表示地层水、地层油粘度,mPa·s;p+γwz、p+γoz分别表示裂缝中水、油位势,Mpa;γw、γo分别表示水、油重度,106×N/m3;p表示压力,Mpa;R表示单位体积基质岩块可采储量,m3;t*表示基质岩块渗吸半周期,d;λR表示储层热导系数,KJ/m/d/℃;T表示储层温度,℃;ρwCw、ρoCo、ρRCR分别表示水、油、岩石体积比热,KJ/m3/℃;ρw、ρo、ρR分别表示水、油密度Kg/m3;Cw、Co、CR分别表示水、油、岩石质量比热,KJ/Kg/℃;pi表示油藏初始压力,为常数,MPa;Ti、Tinj、T0分别表示油藏初始温度、水井注入水温度和初始注入水温度,其中Ti、T0为常数,℃;lb、ub、n分别表示侧边界、顶底边界和垂直于边界方向;pinj、ppro、△p分别表示水井压力、油井压力和注采压差,△p为常数,MPa;Lx、Ly、Lz分别表示室内模型或油藏模拟区域x、y、z方向的特征长度,m;t'=Lxφ/[(△p/Lx)·(Kx/μw)]表示室内模型或油藏模拟区域的特征时间,d;t表示时间,d;τ是求时间积分时引入符号,表示时间微元;e为常数,约等于2.718281828;
所述步骤6)中获取的相似准数为20个,具体如下:
其中,π1~π16为室内模拟常规注水开发所需准数,π17~π20为室内模拟注热水开发时考虑温度及能量变化所附加的准数。
2.如权利要求1所述的一种裂缝性稠油油藏热水驱物理模拟相似准数的确定方法,其特征在于:所述步骤1)所建立的基本假设条件如下:①储层介质包括基质、裂缝,考虑双孔单渗;②流体包括油、水两相;③油藏中任一微元,达到热平衡与相平衡;④忽略裂缝毛管力,考虑基质渗吸,考虑重力;⑤忽略岩石、流体压缩性和热膨胀,忽略烃类热裂解。
3.如权利要求1所述的一种裂缝性稠油油藏热水驱物理模拟相似准数的确定方法,其特征在于:所述步骤3)中将裂缝稠油油藏热水驱数学模型中的相关变量无量纲化的过程如下:
4.如权利要求3所述的一种裂缝性稠油油藏热水驱物理模拟相似准数的确定方法,其特征在于:所述步骤4)中的列出的物理量如下:包括4个自变量x、y、z、t,4个因变量p+γwz、p+γoz-pi、Sw、T-Ti,20个参量Kx、Ky、Kz、μw、μo、γw-γo、R、t*、λR、ρwCw、ρoCo、ρRCR、△p、T0-Ti、φ、rw、Lx、Ly、Lz、t'。
5.如权利要求4所述的一种裂缝性稠油油藏热水驱物理模拟相似准数的确定方法,其特征在于:所述步骤5)中列出的基本量纲如下:质量M、长度L、温度T、时间t。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610305492.4A CN105888631B (zh) | 2016-05-10 | 2016-05-10 | 一种裂缝性稠油油藏热水驱物理模拟相似准数的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610305492.4A CN105888631B (zh) | 2016-05-10 | 2016-05-10 | 一种裂缝性稠油油藏热水驱物理模拟相似准数的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105888631A CN105888631A (zh) | 2016-08-24 |
CN105888631B true CN105888631B (zh) | 2018-06-26 |
Family
ID=56702792
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610305492.4A Active CN105888631B (zh) | 2016-05-10 | 2016-05-10 | 一种裂缝性稠油油藏热水驱物理模拟相似准数的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105888631B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107741605A (zh) * | 2017-08-17 | 2018-02-27 | 中国海洋石油总公司 | 基于时间推移的微元化导电模型求水淹层相关参数的方法 |
CN108681638B (zh) * | 2018-05-16 | 2020-07-07 | 中国地质大学(武汉) | 一种底水砂岩油藏水平井模型物理模拟实验构建方法 |
CN111812147B (zh) * | 2020-06-24 | 2022-03-22 | 浙江大学 | 一种含热源土壤热湿耦合传递模化实验的设计方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4319635A (en) * | 1980-02-29 | 1982-03-16 | P. H. Jones Hydrogeology, Inc. | Method for enhanced oil recovery by geopressured waterflood |
CN101839127A (zh) * | 2010-04-12 | 2010-09-22 | 北京东方亚洲石油技术服务有限公司 | 一种稠油型油藏开采方法 |
CN101899972A (zh) * | 2010-06-30 | 2010-12-01 | 中国石油大学(北京) | 砂岩油藏注水开发可预测物理模型建立方法 |
CN101942991A (zh) * | 2010-06-30 | 2011-01-12 | 中国石油大学(北京) | 裂缝各向异性油藏注水开发可预测物理模型建立方法 |
CN103835708A (zh) * | 2014-02-25 | 2014-06-04 | 中国海洋石油总公司 | 一种油气田类比开发指标定量预测方法 |
-
2016
- 2016-05-10 CN CN201610305492.4A patent/CN105888631B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4319635A (en) * | 1980-02-29 | 1982-03-16 | P. H. Jones Hydrogeology, Inc. | Method for enhanced oil recovery by geopressured waterflood |
CN101839127A (zh) * | 2010-04-12 | 2010-09-22 | 北京东方亚洲石油技术服务有限公司 | 一种稠油型油藏开采方法 |
CN101899972A (zh) * | 2010-06-30 | 2010-12-01 | 中国石油大学(北京) | 砂岩油藏注水开发可预测物理模型建立方法 |
CN101942991A (zh) * | 2010-06-30 | 2011-01-12 | 中国石油大学(北京) | 裂缝各向异性油藏注水开发可预测物理模型建立方法 |
CN103835708A (zh) * | 2014-02-25 | 2014-06-04 | 中国海洋石油总公司 | 一种油气田类比开发指标定量预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105888631A (zh) | 2016-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108830020B (zh) | 一种模拟海上油田微压裂增注裂缝扩展的方法 | |
CN107762498B (zh) | 一种致密气藏直井体积压裂二区的压力分析方法 | |
CN105888631B (zh) | 一种裂缝性稠油油藏热水驱物理模拟相似准数的确定方法 | |
CN105089595B (zh) | 水平压裂裂缝导流作用下的油藏数值模拟方法及装置 | |
CN101942991A (zh) | 裂缝各向异性油藏注水开发可预测物理模型建立方法 | |
CN104389568B (zh) | 蒸汽辅助重力泄油过程中气体辅助用量的获取方法及装置 | |
CN109025817A (zh) | 干热岩单井双水平人工致裂换热方法 | |
CN105319329B (zh) | 致密砂岩气藏水锁损害及电加热解除水锁模拟装置与方法 | |
CN103590798B (zh) | 一种超稠油注蒸汽采油焖井时间的确定方法及计算装置 | |
Du et al. | Mechanism and feasibility study of nitrogen assisted cyclic steam stimulation for ultra-heavy oil reservoir | |
Han et al. | Analyses of the thermo-hydro-mechanical responses of energy pile subjected to non-isothermal heat exchange condition | |
RU2494242C1 (ru) | Способ разработки залежи высоковязкой нефти с использованием внутрипластового горения | |
CN110348031A (zh) | 水平井压裂近井筒裂缝扭曲形态数值模拟方法 | |
CN109187633A (zh) | 地热藏井网式可变裂缝体积干热岩模拟装置及实验方法 | |
CN106761613B (zh) | 一种二氧化碳驱替前缘的试井确定方法 | |
CN102682195A (zh) | 半潜式平台瞬态钻井井筒温度计算方法 | |
CN110489911B (zh) | 一种考虑工质储存过程的增强型地热系统的数值模拟方法及应用 | |
Liu et al. | Treatment technology of steam breakthrough for offshore heavy oil reservoirs developed by huff and puff with multi-component thermal fluid | |
CN105422086B (zh) | 吸汽剖面获取方法及装置 | |
Gu et al. | Numerical study of dynamic fracture aperture during production of pressure-sensitive reservoirs | |
Rabinovich | An analytical solution for cyclic flow of two immiscible phases | |
He et al. | An integrated model for productivity prediction of cyclic steam stimulation with horizontal well | |
Lavrov et al. | Numerical modeling of tensile thermal stresses in rock around a cased well caused by injection of a cold fluid | |
Yang et al. | Numerical assessments of the effects of injected air temperature and well configuration on the cycle performance in compressed air energy storage in aquifers | |
Wilson | Multiscale simulation of WAG flooding in naturally fractured reservoirs |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information |
Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Applicant after: China Offshore Oil Group Co., Ltd. Applicant after: CNOOC research institute limited liability company Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Applicant before: China National Offshore Oil Corporation Applicant before: CNOOC Research Institute |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |