CN111950112A - 一种适用于底部封闭的碳酸盐岩储层动态分析方法 - Google Patents
一种适用于底部封闭的碳酸盐岩储层动态分析方法 Download PDFInfo
- Publication number
- CN111950112A CN111950112A CN201910409037.2A CN201910409037A CN111950112A CN 111950112 A CN111950112 A CN 111950112A CN 201910409037 A CN201910409037 A CN 201910409037A CN 111950112 A CN111950112 A CN 111950112A
- Authority
- CN
- China
- Prior art keywords
- pressure
- dimensionless
- reservoir
- closed
- curve
- 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.)
- Pending
Links
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 title claims abstract description 74
- 238000004458 analytical method Methods 0.000 title claims abstract description 35
- 238000007789 sealing Methods 0.000 title description 2
- 235000019994 cava Nutrition 0.000 claims abstract description 52
- 238000012360 testing method Methods 0.000 claims abstract description 38
- 238000000034 method Methods 0.000 claims abstract description 22
- 238000013178 mathematical model Methods 0.000 claims abstract description 13
- 230000009466 transformation Effects 0.000 claims abstract description 11
- 239000007788 liquid Substances 0.000 claims description 43
- 238000003860 storage Methods 0.000 claims description 40
- 230000007704 transition Effects 0.000 claims description 21
- 230000035699 permeability Effects 0.000 claims description 20
- 239000012530 fluid Substances 0.000 claims description 18
- 230000005484 gravity Effects 0.000 claims description 10
- 238000011084 recovery Methods 0.000 claims description 9
- 230000006835 compression Effects 0.000 claims description 7
- 238000007906 compression Methods 0.000 claims description 7
- 238000004088 simulation Methods 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 2
- 238000011161 development Methods 0.000 abstract description 7
- 239000003208 petroleum Substances 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 7
- 239000003129 oil well Substances 0.000 description 7
- 239000011435 rock Substances 0.000 description 7
- 230000002500 effect on skin Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000011324 bead Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000005206 flow analysis Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明提供了一种适用于底部封闭的碳酸盐岩储层动态分析方法,属于石油勘探开发领域。该方法包括:步骤1,建立底部封闭的碳酸盐岩储层试井数学模型,所述底部封闭的碳酸盐岩储层包括n个溶洞,溶洞的下方为封闭边界;n为大于或等于1的自然数;当n>1时,溶洞之间通过裂缝区连接;步骤2,对所述底部封闭的碳酸盐岩储层试井数学模型进行无量纲化,并求解得到n+1个裂缝区内的压力;步骤3,根据第一个裂缝区内的压力获得无量纲井底压力;步骤4,对所述无量纲井底压力进行反演得到真实空间下的井底压力解;步骤5,根据真实空间下的井底压力解绘制压力曲线和压力导数曲线的双对数曲线;步骤6,进行动态分析。
Description
技术领域
本发明属于石油勘探开发领域,具体涉及一种适用于底部封闭的碳酸盐岩储层动态分析方法。
背景技术
碳酸盐岩油藏在世界已发现油藏中占据重要位置,全世界256个大型油田中,碳酸盐岩油田115个,占45%,已发现的碳酸盐岩油藏储量和产量占世界总量的50%、65%,在油田开发中一直占有重要地位。国外碳酸盐岩油藏主要分布在中东和美洲地区。
国内塔河油田的发现是碳酸盐岩油藏新的重大突破,揭开了我国碳酸盐岩油藏勘探开发新阶段。塔河油田的主体为奥陶系碳酸盐岩缝洞型油藏,主要储层类型为缝洞型储层,是多期岩溶改造作用的结果,储集空间主要有溶洞、孔洞、裂缝等,由这些特征明显不同的储集空间组合形成溶洞型、裂缝-孔洞型、裂缝型以及洞穴型储集体,其在三维空间分布的边界形态极不规则;储集空间分布不连续,孔隙度变化巨大、规律性差,非均质性非常严重。洞穴型储层储集空间为大型洞穴和裂缝,洞穴(包括大洞、巨洞)储集空间巨大,加之裂缝对沟通洞穴和改善渗流性能的作用,形成了储集空间巨大、储渗能力极好的有利储层类型。试井分析可以反演得到缝洞储层的溶洞体积、连通性、储量等关键信息,对确定碳酸盐岩缝洞型油藏储量、指导碳酸盐岩油藏具有重要意义。
目前缝洞型油藏主要有两大类试井分析模型:第一类是连续介质模型,主要包括双重介质模型、三重介质模型或等效三重介质模型、双重复合模型,如“变井筒储存的三重介质油藏试井解释方法研究”、“碳酸盐岩溶洞型储层试井解释新模型”、“井打在大尺度溶洞内的缝洞型油藏试井模型”等文献中公开的模型。第二类是离散介质模型试井模型,如“塔河碳酸盐岩油藏试井解释新技术应用”、“井打在溶洞外的缝洞型油藏数值试井模型”、中国专利公开文献CN106599449A“一种溶洞体积计算的试井解释方法”公开的模型。
但现有模型和方法仅适用于溶洞和裂缝平面径向分布的碳酸盐储层。而很多碳酸盐岩储层为串珠型溶洞,溶洞之间靠垂向裂缝连通,流动主要为垂直方向的流动。因此,现有的试井分析方法对垂向流动的缝洞储层(如串珠型溶洞)适用性较差。目前还没有相关的试井分析方法进行垂向缝洞里流体流动的缝洞储层特征参数描述和压力动态分析。
发明内容
本发明的目的在于解决上述现有技术中存在的难题,提供一种适用于底部封闭的碳酸盐岩储层动态分析方法,形成一种底部封闭储层向n个溶洞垂向供液型缝洞型碳酸盐岩储层试井分析方法,即储层中有1个或多个溶洞,这些溶洞通过其上部垂向的裂缝或高渗带向井筒供液,而溶洞的底部某深处的封闭储层还通过垂向裂缝或高渗带向溶洞供液。从碳酸盐岩缝洞储层中溶洞、裂缝及基岩的地质特征出发,建立耦合溶洞、裂缝垂向流动的试井分析模型,为缝洞型碳酸盐岩油藏确定准确可靠的溶洞体积、高渗或裂缝带流动能力及地质储量提供依据,为该类油藏动态评价提供基础信息,对保障缝洞型碳酸盐岩油藏高效开发、提高经济效益具有重要作用。
本发明是通过以下技术方案实现的:
步骤1,建立底部封闭的碳酸盐岩储层试井数学模型,所述底部封闭的碳酸盐岩储层包括n个溶洞,溶洞的下方为封闭边界;n为大于或等于1的自然数;当n>1时,溶洞之间通过裂缝区连接;
步骤2,对所述底部封闭的碳酸盐岩储层试井数学模型进行无量纲化,并求解得到n+1个裂缝区内的压力;
步骤3,根据第一个裂缝区内的压力获得无量纲井底压力;
步骤4,对所述无量纲井底压力进行反演得到真实空间下的井底压力解;
步骤5,根据真实空间下的井底压力解绘制压力曲线和压力导数曲线的双对数曲线;
步骤6,利用所述压力曲线和压力导数曲线的双对数曲线进行动态分析。
n+1个所述裂缝区均为高渗条带,裂缝区是溶洞与井筒之间、溶洞与溶洞之间、地层与溶洞之间的唯一流动通道。
所述步骤1建立的底部封闭的碳酸盐岩储层试井数学模型如下:
p(z,t=0)=pi(2)
式中:
p为压力,Pa;
pi为第i个裂缝区中的压力,i=1,2,...,n+1,Pa;从上至下依次为第一个裂缝区、第二个裂缝区……、第n+1裂缝区;
z为垂向坐标位置,z轴为以井筒的底端面圆心为原点向下建立的一维坐标轴;
ρ为流体密度,kg/m3;
g为重力加速度,m/s2;
cf为流体压缩系数,1/Pa;
φ为孔隙度,无量纲;
μ为流体粘度,Pa·s;
Ct为综合压缩系数,1/Pa;
k为渗透率,md;
t为时间,s。
pi为原始地层压力,Pa;
H2n-1为第n个溶洞的上边界距离原点的位置,m;
h2n为第n个溶洞的下边界距离原点的位置;
pvn为第n个溶洞中的压力,Pa;
rf为溶洞的半径,m;
Cvn为第n个溶洞的储集系数,m3/Pa;
q为日产量,m3/s;
B为流体体积系数。
所述步骤2中的对所述底部封闭的碳酸盐岩储层试井数学模型进行无量纲化的操作包括:
利用下面的无量纲参数对公式(1)到(6)进行无量纲化得到无量纲化后的模型:
无量纲压力pD:
无量纲时间tD:
无量纲距离zD:
无量纲井筒系数CD:
无量纲溶洞储集系数CvD:
式中:C为井筒存储常数,m3/Pa。
无量纲重力系数GD:
GD=2ρgcf·rf。
所述步骤2中的求解得到n+1个裂缝区内的压力的操作包括:
对所述无量纲化后的模型进行Laplace变换,得到在Laplace空间上的n+1个裂缝区中的无量纲的渗流控制方程;
对所述在Laplace空间上的n+1个裂缝区中的无量纲的渗流控制方程进行求解,得到Laplace空间上无量纲的n+1个裂缝区内的压力。
所述步骤3的操作包括:
在ZD=0时,所述Laplace空间上无量纲的第一个裂缝区内的压力即为拉式空间上的无量纲井底压力。
所述步骤4的操作包括:
采用Stehfest数值反演算法对所述拉式空间上的无量纲井底压力进行反演得到真实空间下的井底压力解,即井底压力与时间的关系。
所述步骤6的操作包括:
判断储层中是否存在溶洞:如果压力曲线和压力导数曲线的双对数曲线中存在过渡流段和溶洞存储段,则判定储层中存在溶洞,否则判定储层中不存在溶洞;
如果n>1,所述步骤6的操作还包括获取储层中的溶洞的数量:在压力曲线和压力导数曲线的双对数曲线中统计出现过渡流段和溶洞存储段的次数,统计得到的次数即为溶洞的数量;
所述过渡流段是指压力导数曲线中骤降的部分;
所述溶洞存储段是指压力导数曲线中与所述过渡流段相连的斜率为1的部分。
所述步骤6的操作包括:
判断溶洞以下的深部储层是否垂向供液:如果压力曲线和压力导数曲线的双对数曲线中存在边界流段,则判定溶洞以下的深部储层垂向供液,否则判定溶洞以下的深部储层不垂向供液;
所述边界流段是指压力曲线和压力导数曲线的斜率都趋向于1的部分。
所述步骤6的操作包括:根据模拟数据和实测数据进行历史拟合得到底部供液的深度、溶洞的大小和高渗带导流能力参数;所述模拟数据是利用所述真实空间下的井底压力解获得的压力和压力导数,所述实测数据是指在压力恢复测试中实际测量得到的压力和压力导数。
与现有技术相比,本发明的有益效果是:本发明形成了一种底部封闭储层向n个溶洞垂向供液型缝洞型碳酸盐岩储层流体流动的试井分析方法,储层中有1个或多个溶洞,这些溶洞通过其上部垂向的裂缝(或高渗透带)向井筒供液,而溶洞的下部深处储层还通过垂向裂缝(或高渗透带)向溶洞供液。从碳酸盐岩缝洞储层中溶洞、裂缝及基岩的地质特征出发,建立耦合溶洞、裂缝垂向流动的试井分析模型。该方法可通过过渡流段和溶洞存储段判断是否存在溶洞及溶洞的个数,对定性判断溶洞的存在提供了技术手段;通过无量纲引入重力系数考虑了垂向流动中的重力影响,可更准确的反映底部封闭储层供液特征;对于1个溶洞的储层,可反映溶洞下部某深处的封闭储层向溶洞供液的碳酸盐岩油井呈现的6个完整流态,对于多个溶洞的储层,可反映底部封闭储层向溶洞供液的碳酸盐岩油井呈现的9个完整流态,有效的反映了井筒存储、溶洞向井筒供液、溶洞里流动、底部地层向溶洞供液等流动的全过程,因此对底部供液型的碳酸盐岩油井适应性好,能用于定量确定该类储层溶洞体积等关键参数。因此该方法反映溶洞中的垂向流体流动和深部地层垂向供液特征,进行裂缝和溶洞垂向分布的缝洞型碳酸盐岩储层流动动态分析,为缝洞型碳酸盐岩油藏确定准确可靠的溶洞体积、高渗带流动能力及地质储量提供依据,为该类油藏动态评价提供基础信息,对保障缝洞型碳酸盐岩油藏高效开发、提高经济效益具有重要作用。
附图说明
图1底部封闭储层向n个溶洞垂向供液型缝洞型碳酸盐岩储层模型示意图;当n=1时,图1中仅仅有高渗区1、溶洞和高渗区2。
图2-1本发明实施例1中得到压力曲线和压力导数曲线的双对数曲线
图2-2本发明实施例2中得到压力曲线和压力导数曲线的双对数曲线;
图3本发明方法的步骤框图。
具体实施方式
下面结合附图对本发明作进一步详细描述:
本发明形成了一种底部封闭储层向n个溶洞垂向供液型缝洞型碳酸盐岩储层试井分析方法,从碳酸盐岩缝洞储层中溶洞、裂缝及基岩的地质特征出发,建立耦合溶洞、裂缝垂向流动的试井分析模型。该模型考虑重力影响、高渗带的流动、溶洞里的流动等流动。定义了所需的无量纲参数,通过Laplace变换和数值反演得到模型的解,形成了缝洞型碳酸盐岩储层试井分析方法,为缝洞型碳酸盐岩油藏确定准确可靠的溶洞体积及地质储量提供依据,为该类油藏动态评价提供基础信息,对保障缝洞型碳酸盐岩油藏高效开发、提高经济效益具有重要作用。封闭边界(或称为封闭储层)表示储层内没有压力供给,储层内的压力只能是衰竭式降低。
(一)底部封闭储层向多个溶洞垂向供液型缝洞型碳酸盐岩储层试井模型物理假设
本发明采用以下假设:
1.地层中有n个溶洞,溶洞之间通过裂缝连接,且溶洞性质保持不变;
2.溶洞为圆柱体,半径为rf,溶洞上边界为h2n-1,下边界为h2n;
3.溶洞处的压力处处相等;
4.高渗条带由裂缝区表示,裂缝区为有限导流能力裂缝,基质渗透率与裂缝渗透率相比可以忽略,裂缝区是溶洞与井筒之间、溶洞与溶洞之间、地层与溶洞之间的唯一流动通道;
5.溶洞间的裂缝系统简化成一个圆柱形区域,区域内的渗透率等效为该区域中的裂缝渗透率;
6.油藏温度不变,原始地层压力处处相等;
7.每个区域内的性质保持不变;
8.地层中的流体为单相流体;
9.流体和岩石均为微可压,且压缩系数均为常数;
10.模型的下边界为垂向上的封闭边界;
11.考虑井储、表皮以及重力效应。
基于上述假设,该模型中的流动可以简化为一维垂向流动。
(二)底部封闭储层向n个溶洞垂向供液型缝洞型碳酸盐岩储层试井数学模型
在一维直角坐标系下,考虑重力的影响,根据质量守恒定律,得出n+1个裂缝区(即图1中的高渗区)中的渗流控制方程分别为:
式中:
p为压力,Pa;
pi为第i个裂缝区中的压力,i=1,2,...,n+1,Pa;从上至下依次为第一个裂缝区、第二个裂缝区……、第n+1个裂缝区;
z为垂向坐标位置,z轴为以井筒的底端面圆心为原点向下建立的一维坐标轴;
ρ为流体密度,kg/m3;
g为重力加速度,m/s2;
cf为流体压缩系数,1/Pa;
φ为孔隙度,无量纲;
μ为流体粘度,Pa·s;
Ct为综合压缩系数,1/Pa;
k为渗透率,md;
t为时间,s。
初始条件为:
p(z,t=0)=pi (2)
式中:pi为原始地层压力,Pa;
外边界条件:垂向上为封闭边界:
假设每个溶洞里的压力处处相等,则有:
式中:
H2n-1为第n个溶洞的上边界距离原点的位置,m;
h2n为第n个溶洞的下边界距离原点的位置;
pvn为第n个溶洞中的压力,Pa;
设溶洞的储集系数为Cv,则根据质量守恒定律,得出溶洞上下边界处的流量关系:
式中:
rf为溶洞的半径,m;
Cvn为第n个溶洞的储集系数,m3/Pa;
根据质量守恒和达西定律,得出在裂缝区域1与井筒相交处的流量关系为:
式中:
q为日产量,m3/s;
B为流体体积系数。
(三)无量纲参数定义
无量纲参数定义如下:
无量纲压力:
无量纲时间:
无量纲距离:
无量纲井筒系数:
无量纲溶洞储集系数:
式中:C为井筒存储常数,m3/Pa。
无量纲重力系数:
GD=2ρgcf·rf
(四)方程无量纲及模型Laplace变换及求解
通过定义的上述无量纲变量,先将试井模型进行无量纲化,再进行Laplace变换,得到在Laplace空间上的n+1个裂缝区中的无量纲的渗流控制方程。
对上述在Laplace空间上的n+1个裂缝区中的无量纲的渗流控制方程进行求解,得到Laplace空间上无量纲的n+1个裂缝区内的压力。
因为在物理建模时是以井底位置(即井筒的底端面、第一裂缝区域的上端面)作为坐标原点的,所以当取ZD=0时,所述Laplace空间上无量纲的第一裂缝区内的压力即为拉式空间上的无量纲井底压力。
(五)模型在真实空间求解
采用Stehfest数值反演算法对所述拉式空间上的无量纲井底压力进行反演得到真实空间下的井底压力解,即井底压力与时间的关系。基于真实空间下的井底压力解,可绘制模型曲线,所述模型曲线包括压力曲线和压力导数曲线。
本发明方法如图3所示,包括:
步骤1,建立底部封闭的碳酸盐岩储层试井数学模型,所述底部封闭的碳酸盐岩储层包括n个溶洞,溶洞的下方为封闭边界;n为大于或等于1的自然数;当n>1时,溶洞之间通过裂缝区连接;
步骤2,对所述底部封闭的碳酸盐岩储层试井数学模型进行无量纲化,并求解得到n+1个裂缝区内的压力;
步骤3,根据第一个裂缝区内的压力获得无量纲井底压力;
步骤4,对所述无量纲井底压力进行反演得到真实空间下的井底压力解;
步骤5,根据真实空间下的井底压力解绘制压力曲线和压力导数曲线的双对数曲线;
步骤6,利用所述压力曲线和压力导数曲线的双对数曲线进行动态分析,包括是否存在溶洞,溶洞的个数,底部储层是否垂向供液,通过压力恢复测试实测数据与模型数据拟合得到溶洞大小、高渗带导流能力等特征参数。
本发明的实施例如下:
实施例1:
以某碳酸盐岩油藏中一口井为例,在距离井底50米处有一个溶洞,溶洞通过高渗条带向井供液,而距离溶洞底部以下500米处的封闭储层还通过高渗条带向溶洞供液,高渗条带的渗透率为1000md。模拟该井压力恢复,得到压力恢复过程中的压力和压力导数的双对数曲线(即根据井底压力解绘制的压力曲线和压力导数曲线),如图2-1所示。从双对数曲线图2-1可以看出,底部储层向溶洞供液的碳酸盐岩油井呈现6个流态段,分别是井筒存储段、表皮效应段、线性流段、过渡流段、溶洞存储段以及边界流。流态1对应纯井筒存储段,其特征为压力曲线和压力导数曲线斜率均为1;流态2对应表皮效应段,其特征为导数曲线上有一个峰,且表皮越大,峰值越大;流态3对应线性流段,代表在高渗条带中的流动,其特征为压力导数曲线斜率为1/2;流态4对应过渡流段,代表从高渗条带向溶洞流动转变,流动通道的突然改变导致导数曲线的骤降;流态5对应溶洞存储段,代表在溶洞1中的流动,压力导数曲线斜率为1;流态6对应边界流段,压力曲线和压力导数曲线斜率都趋向于1,反映了底部封闭储层通过高渗条带向溶洞不断供液的特征。流态4和5是由于溶洞的存在而出现的,因此通过流态4和5的出现可以判断是否有溶洞存在。
通过这6个流态的判断,有效的反映了井筒存储、溶洞向井筒供液、溶洞里流动、底部地层向溶洞供液等流动特征,因此,本发明提出的方法可进行底部供液溶洞型碳酸盐岩的流动动态分析(通过压力曲线和压力导数曲线得到各个流态),可以通过压力恢复测试实测数据与模型数据拟合得到溶洞大小(根据模拟数据和实测数据进行历史拟合,具体就是用实测数据的压力和压力导数,分别去拟合模型数据的压力和压力导数曲线,当得到很好的拟合效果时,就认为模型给定的溶洞大小为实际溶洞大小)、高渗带导流能力等特征参数(与上述求溶洞大小的方法相同,采用历史拟合得到),进而判断储量规模,也可以用于生产过程中碳酸盐岩油气井产量特征、压力特征预测分析。
实施例2:
以某碳酸盐岩油藏中一口井为例,该井在垂向上连通了2个溶洞,分别在距离井底10米和50米的位置,垂向上100米处碰到封闭边界。溶洞通过高渗条带向井供液,而溶洞以下的储层还通过高渗条带向溶洞供液,高渗条带的渗透率为1000md。模拟该井压力恢复,得到压力恢复过程中的压力和压力导数的双对数曲线,即根据真实空间下的井底压力解绘制的压力曲线和压力导数曲线的双对数曲线,如图2-2所示。从双对数曲线图2-2可以看出,底部深部储层向溶洞供液的碳酸盐岩油井呈现9个流态,分别是井筒存储段,表皮效应段,第一次线性流段,第一次过渡流段,第一次溶洞存储段,第二次线性流段,第二次过渡流段,第二次溶洞存储段和边界流段。流态1对应纯井筒存储段,其特征为压力曲线和压力导数曲线斜率均为1;流态2对应表皮效应段,其特征为导数曲线上有一个峰,且表皮越大,峰值越大;流态3对应第一次线性流段,代表在高渗条带1中的流动,其特征为压力导数曲线斜率为1/2;流态4对应第一次过渡流段,代表从高渗条带1向溶洞1流动转变,流动通道的突然改变导致导数曲线的骤降;流态5对应第一次溶洞存储段,代表在溶洞I中的流动,压力导数曲线斜率为1;流态6对应第二次线性流段,代表在高渗条带2中的流动;流态7对应第二次过渡流段,代表流动从高渗条带2向溶洞2中过渡;流态8对应第二次溶洞存储段,即流动到达溶洞2;流态9对应边界流段,压力曲线和压力导数曲线斜率都趋向于1,反映了底部封闭储层通过高渗条带向溶洞供液的特征。
流态4和5是由于溶洞1的存在而出现的,流态7和8是由于溶洞2的存在而出现的,因此通过流态判断是否有溶洞存在,以及有几个溶洞存在。具体的,只有先后出现了过渡流段和溶洞存储段2个流态,才能判断出现了一个溶洞。统计过渡流段和溶洞存储段2个流态出现的次数即可得到溶洞的数量。统计时,第一次出现过度流段和溶洞存储段,则次数加一,第二次出现过度流段和溶洞存储段,则次数再加一,以此类推。
通过这9个流态的判断,有效的反映了井筒存储、溶洞向井筒供液、溶洞里流动、底部地层向溶洞供液等流动特征,因此,本发明提出的方法可进行底部供液溶洞型碳酸盐岩的流动动态分析(通过9个流态的分析来实现),可以通过压力恢复测试实测数据与模型数据拟合得到溶洞大小(根据模拟数据和实测数据进行历史拟合,具体是用实测数据中的压力和压力导数,分别去拟合模型数据的压力和压力导数曲线,当得到很好的拟合效果时,就认为模型给定的溶洞大小为实际溶洞大小)、高渗带导流能力等特征参数(与上述求溶洞大小的方法类似,采用历史拟合获得。),进而判断储量规模,也可以用于生产过程中碳酸盐岩油气井产量特征、压力特征预测分析。
上述实施例实施后,达到了如下效果:(1)可通过过渡流段和溶洞存储段判断是否存在溶洞以及溶洞个数,对定性判断溶洞的存在提供了技术手段;(2)考虑了垂向流动中的重力影响,可更准确的反映底部封闭储层供液特征;(3)该方法可反映底部封闭储层向溶洞供液的碳酸盐岩油井呈现的9个完整流态,有效的反映了井筒存储、溶洞向井筒供液、溶洞里流动、底部地层向溶洞供液等流动的全过程,因此对底部供液型的碳酸盐岩油井适应性好,能用于定量确定该类储层溶洞体积等关键参数。
上述技术方案只是本发明的一种实施方式,对于本领域内的技术人员而言,在本发明公开了应用方法和原理的基础上,很容易做出各种类型的改进或变形,而不仅限于本发明上述具体实施方式所描述的方法,因此前面描述的方式只是优选的,而并不具有限制性的意义。
Claims (10)
1.一种适用于底部封闭的碳酸盐岩储层动态分析方法,其特征在于:所述方法包括:
步骤1,建立底部封闭的碳酸盐岩储层试井数学模型,所述底部封闭的碳酸盐岩储层包括n个溶洞,溶洞的下方为封闭边界;n为大于或等于1的自然数;当n>1时,溶洞之间通过裂缝区连接;
步骤2,对所述底部封闭的碳酸盐岩储层试井数学模型进行无量纲化,并求解得到n+1个裂缝区内的压力;
步骤3,根据第一个裂缝区内的压力获得无量纲井底压力;
步骤4,对所述无量纲井底压力进行反演得到真实空间下的井底压力解;
步骤5,根据真实空间下的井底压力解绘制压力曲线和压力导数曲线的双对数曲线;
步骤6,利用所述压力曲线和压力导数曲线的双对数曲线进行动态分析。
2.根据权利要求1所述的适用于底部封闭的碳酸盐岩储层动态分析方法,其特征在于:n+1个所述裂缝区均为高渗条带,裂缝区是溶洞与井筒之间、溶洞与溶洞之间、地层与溶洞之间的唯一流动通道。
3.根据权利要求1所述的适用于底部封闭的碳酸盐岩储层动态分析方法,其特征在于:所述步骤1建立的底部封闭的碳酸盐岩储层试井数学模型如下:
p(z,t=0)=pi (2)
式中:
p为压力,Pa;
pi为第i个裂缝区中的压力,i=1,2,...,n+1,Pa;从上至下依次为第一个裂缝区、第二个裂缝区……、第n+1裂缝区;
z为垂向坐标位置,z轴为以井筒的底端面圆心为原点向下建立的一维坐标轴;
ρ为流体密度,kg/m3;
g为重力加速度,m/s2;
cf为流体压缩系数,1/Pa;
φ为孔隙度,无量纲;
μ为流体粘度,Pa·s;
Ct为综合压缩系数,1/Pa;
k为渗透率,md;
t为时间,s;
pi为原始地层压力,Pa;
H2n-1为第n个溶洞的上边界距离原点的位置,m;
h2n为第n个溶洞的下边界距离原点的位置;
pvn为第n个溶洞中的压力,Pa;
rf为溶洞的半径,m;
Cvn为第n个溶洞的储集系数,m3/Pa;
q为日产量,m3/s;
B为流体体积系数。
5.根据权利要求4所述的适用于底部封闭的碳酸盐岩储层动态分析方法,其特征在于:所述步骤2中的求解得到n+1个裂缝区内的压力的操作包括:
对所述无量纲化后的模型进行Laplace变换,得到在Laplace空间上的n+1个裂缝区中的无量纲的渗流控制方程;
对所述在Laplace空间上的n+1个裂缝区中的无量纲的渗流控制方程进行求解,得到Laplace空间上无量纲的n+1个裂缝区内的压力。
6.根据权利要求4所述的适用于底部封闭的碳酸盐岩储层动态分析方法,其特征在于:所述步骤3的操作包括:
在ZD=0时,所述Laplace空间上无量纲的第一个裂缝区内的压力即为拉式空间上的无量纲井底压力。
7.根据权利要求6所述的适用于底部封闭的碳酸盐岩储层动态分析方法,其特征在于:所述步骤4的操作包括:
采用Stehfest数值反演算法对所述拉式空间上的无量纲井底压力进行反演得到真实空间下的井底压力解,即井底压力与时间的关系。
8.根据权利要求1所述的适用于底部封闭的碳酸盐岩储层动态分析方法,其特征在于:所述步骤6的操作包括:
判断储层中是否存在溶洞:如果压力曲线和压力导数曲线的双对数曲线中存在过渡流段和溶洞存储段,则判定储层中存在溶洞,否则判定储层中不存在溶洞;
如果n>1,所述步骤6的操作还包括获取储层中的溶洞的数量:在压力曲线和压力导数曲线的双对数曲线中统计出现过渡流段和溶洞存储段的次数,统计得到的次数即为溶洞的数量;
所述过渡流段是指压力导数曲线中骤降的部分;
所述溶洞存储段是指压力导数曲线中与所述过渡流段相连的斜率为1的部分。
9.根据权利要求1所述的适用于底部封闭的碳酸盐岩储层动态分析方法,其特征在于:所述步骤6的操作包括:
判断溶洞以下的封闭储层是否垂向供液:如果压力曲线和压力导数曲线的双对数曲线中存在边界流段,则判定溶洞以下的封闭储层垂向供液,否则判定溶洞以下的封闭储层不垂向供液;
所述边界流段是指压力曲线和压力导数曲线的斜率都趋向于1的部分。
10.根据权利要求1所述的适用于底部封闭的碳酸盐岩储层动态分析方法,其特征在于:所述步骤6的操作包括:根据模拟数据和实测数据进行历史拟合得到底部供液的深度、溶洞的大小和高渗带导流能力参数;所述模拟数据是利用所述真实空间下的井底压力解获得的压力和压力导数,所述实测数据是指在压力恢复测试中实际测量得到的压力和压力导数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910409037.2A CN111950112A (zh) | 2019-05-16 | 2019-05-16 | 一种适用于底部封闭的碳酸盐岩储层动态分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910409037.2A CN111950112A (zh) | 2019-05-16 | 2019-05-16 | 一种适用于底部封闭的碳酸盐岩储层动态分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111950112A true CN111950112A (zh) | 2020-11-17 |
Family
ID=73336687
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910409037.2A Pending CN111950112A (zh) | 2019-05-16 | 2019-05-16 | 一种适用于底部封闭的碳酸盐岩储层动态分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111950112A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112282741A (zh) * | 2019-07-24 | 2021-01-29 | 中国石油化工股份有限公司 | 目标储层试井分析方法、计算机存储介质和计算机设备 |
CN112377183A (zh) * | 2020-11-30 | 2021-02-19 | 中国石油天然气股份有限公司 | 多洞多缝型碳酸盐岩储层产量分析模型构建方法及装置 |
CN112377182A (zh) * | 2020-11-30 | 2021-02-19 | 中国石油天然气股份有限公司 | 大型孔洞裂缝型碳酸盐岩储集体参数确定方法及装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4797821A (en) * | 1987-04-02 | 1989-01-10 | Halliburton Company | Method of analyzing naturally fractured reservoirs |
CN103077558A (zh) * | 2012-05-14 | 2013-05-01 | 中国石油化工股份有限公司 | 碳酸盐岩缝洞型油藏大型溶洞储集体分布模型的建模方法 |
CN105569641A (zh) * | 2015-12-17 | 2016-05-11 | 西南石油大学 | 一种改善近井压降的水平井瞬态压力快速计算模型建立方法 |
US20160363691A1 (en) * | 2015-06-15 | 2016-12-15 | Petrochina Company Limited | Physical simulation method and experiment device of fracture-cavity carbonate reservoir hydrocarbon charge |
US20170191348A1 (en) * | 2016-01-04 | 2017-07-06 | Saudi Arabian Oil Company | Modeling to characterize fractures network in homogeneous petroleum reservoirs |
CN108798634A (zh) * | 2018-07-09 | 2018-11-13 | 合肥辰工科技有限公司 | 一种缝洞型储层单洞-单层均值试井解释方法 |
CN109033519A (zh) * | 2018-06-22 | 2018-12-18 | 中国石油天然气股份有限公司 | 一种异常高压碳酸盐岩挥发性油藏的试井解释方法及装置 |
CN109138974A (zh) * | 2017-06-19 | 2019-01-04 | 中国石油化工股份有限公司 | 一种缝洞型碳酸盐岩油藏离散数值试井分析方法及系统 |
-
2019
- 2019-05-16 CN CN201910409037.2A patent/CN111950112A/zh active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4797821A (en) * | 1987-04-02 | 1989-01-10 | Halliburton Company | Method of analyzing naturally fractured reservoirs |
CN103077558A (zh) * | 2012-05-14 | 2013-05-01 | 中国石油化工股份有限公司 | 碳酸盐岩缝洞型油藏大型溶洞储集体分布模型的建模方法 |
US20160363691A1 (en) * | 2015-06-15 | 2016-12-15 | Petrochina Company Limited | Physical simulation method and experiment device of fracture-cavity carbonate reservoir hydrocarbon charge |
CN105569641A (zh) * | 2015-12-17 | 2016-05-11 | 西南石油大学 | 一种改善近井压降的水平井瞬态压力快速计算模型建立方法 |
US20170191348A1 (en) * | 2016-01-04 | 2017-07-06 | Saudi Arabian Oil Company | Modeling to characterize fractures network in homogeneous petroleum reservoirs |
CN109138974A (zh) * | 2017-06-19 | 2019-01-04 | 中国石油化工股份有限公司 | 一种缝洞型碳酸盐岩油藏离散数值试井分析方法及系统 |
CN109033519A (zh) * | 2018-06-22 | 2018-12-18 | 中国石油天然气股份有限公司 | 一种异常高压碳酸盐岩挥发性油藏的试井解释方法及装置 |
CN108798634A (zh) * | 2018-07-09 | 2018-11-13 | 合肥辰工科技有限公司 | 一种缝洞型储层单洞-单层均值试井解释方法 |
Non-Patent Citations (1)
Title |
---|
邢翠巧 等: ""缝洞型碳酸盐岩油藏试井分析模型研究"", 《河北工业科技》, vol. 35, no. 01, 31 January 2018 (2018-01-31), pages 31 - 36 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112282741A (zh) * | 2019-07-24 | 2021-01-29 | 中国石油化工股份有限公司 | 目标储层试井分析方法、计算机存储介质和计算机设备 |
CN112282741B (zh) * | 2019-07-24 | 2024-05-14 | 中国石油化工股份有限公司 | 目标储层试井分析方法、计算机存储介质和计算机设备 |
CN112377183A (zh) * | 2020-11-30 | 2021-02-19 | 中国石油天然气股份有限公司 | 多洞多缝型碳酸盐岩储层产量分析模型构建方法及装置 |
CN112377182A (zh) * | 2020-11-30 | 2021-02-19 | 中国石油天然气股份有限公司 | 大型孔洞裂缝型碳酸盐岩储集体参数确定方法及装置 |
CN112377182B (zh) * | 2020-11-30 | 2023-08-22 | 中国石油天然气股份有限公司 | 大型孔洞裂缝型碳酸盐岩储集体参数确定方法及装置 |
CN112377183B (zh) * | 2020-11-30 | 2023-09-26 | 中国石油天然气股份有限公司 | 多洞多缝型碳酸盐岩储层产量分析模型构建方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103413030B (zh) | 一种缝洞型碳酸盐岩气藏动态分析方法及系统 | |
CN106150477B (zh) | 一种确定缝洞型油藏的单井控制储量的方法 | |
CN110096718B (zh) | 一种获取碳酸盐岩油藏中溶洞的体积的方法 | |
CN105386751B (zh) | 一种基于油藏渗流模型的水平井测井产能预测方法 | |
CN106599449A (zh) | 一种溶洞体积计算的试井解释方法 | |
CN104533370A (zh) | 压裂水平井油藏、裂缝、井筒全耦合模拟方法 | |
CN111794740B (zh) | 一种适用在缝洞型碳酸盐岩油藏动态储量计算的方法 | |
CN111950112A (zh) | 一种适用于底部封闭的碳酸盐岩储层动态分析方法 | |
CN110107277B (zh) | 一种获取碳酸盐岩油藏中井钻遇溶洞的体积的方法 | |
CN104899411B (zh) | 一种储层产能预测模型建立方法和系统 | |
CN110984973A (zh) | 缝洞型碳酸盐岩气藏单井控制储量的确定方法 | |
CN110362931B (zh) | 一种基于溶洞点源等效原理的油气藏试井解释模型及方法 | |
CN110656915B (zh) | 一种页岩气多段压裂水平井多工作制度产能预测方法 | |
CN108518218B (zh) | 一种非常规油气藏多段压裂水平井单井动态储量确定方法 | |
CN108894777A (zh) | 一种分层压裂多层合采油气藏储层物性及裂缝特性参数的确定方法 | |
CN110439519A (zh) | 一种基于极限限流设计的压裂方法及系统 | |
CN106640021A (zh) | 压后放喷参数的计算方法及装置 | |
CN115587674B (zh) | 油藏改建储气库扩容达产过程气井动态产能预测方法 | |
CN114510854A (zh) | 一种循缝找洞的酸压数值模拟结果准确性评价方法 | |
CN109242364A (zh) | 一种高温高压气井模拟井筒体积置换产能评价法 | |
CN109441415A (zh) | 基于邻井干扰的聚合物驱油藏测试井的试井解释方法 | |
CN114201932A (zh) | 一种复杂情况下的致密油藏压裂井试井模拟方法 | |
CN112576248A (zh) | 一种底水气藏早期产能评价预测方法 | |
CN112035993A (zh) | 一种底部定压的碳酸盐岩储层测试评价方法 | |
CN111950111B (zh) | 一种适用于底部开放的碳酸盐岩储层动态分析方法 |
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 |