CN107577831B - 一种缝洞型碳酸盐岩油气藏溶洞规模大小计算方法 - Google Patents
一种缝洞型碳酸盐岩油气藏溶洞规模大小计算方法 Download PDFInfo
- Publication number
- CN107577831B CN107577831B CN201710159508.XA CN201710159508A CN107577831B CN 107577831 B CN107577831 B CN 107577831B CN 201710159508 A CN201710159508 A CN 201710159508A CN 107577831 B CN107577831 B CN 107577831B
- Authority
- CN
- China
- Prior art keywords
- calculating
- karst cave
- straight line
- curve
- line segment
- 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
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Filling Or Discharging Of Gas Storage Vessels (AREA)
Abstract
本发明属于油气田开发技术领域,具体而言,涉及一种缝洞型碳酸盐岩油气藏溶洞规模大小计算方法,包括以下步骤:A、收集地层、流体的基本参数以及关井压力恢复测试数据;B、将压力恢复测试数据做成双对数曲线和半对数曲线;C、在半对数曲线上找到直线段开始和结束的时间点;D、根据半对数曲线计算直线段斜率,再由斜率计算地层系数kh;E、根据步骤C找到的直线段时间区间,在双对数曲线上找到对应时间区间的凹子区间;F、根据数值积分方法计算凹子面积S;G、根据规模溶洞体积公式计算出现场溶洞体积V。该方法具有较高的可行性和准确性,能准确计算出缝洞型碳酸盐岩油气藏溶洞规模的大小,为现场的油气田开发提供了重要的理论依据。
Description
技术领域
本发明涉及一种缝洞型碳酸盐岩油气藏溶洞规模大小计算方法,属于油气田开发技术领域。
背景技术
试井是认识油气藏的一种重要手段和方法。缝洞型碳酸盐岩油气藏储集空间多样,大溶洞作为一种主要的储集空间,裂缝作为一种主要的渗流通道,同时溶洞与裂缝大小不一,相差悬殊,具有明显的多尺度性,且溶洞、裂缝分布不均,连续性差,因此缝洞型碳酸盐岩油气藏流动机理极其复杂。
在缝洞型碳酸盐岩油藏试井方面,前人提出了各种各样的试井数学模型及其样版曲线,然而由于数学模型参数较多,多解性强,矿场实测曲线与理论曲线却往往难以拟合。
目前解决上述问题的理论研究主要有两方面:一是在常规的双重介质渗流理论的基础上,把流动机理和数学模型主要集中在裂缝型介质系统、孔洞-裂缝-基质系统和裂缝-基质介质系统方面,但涉及大型洞穴-裂缝介质系统的研究较少。二是直接采用流体力学中管流和明渠流动相结合方式来描述,但是不能很好的解释相关的地层参数。这两种方法都难以解决缝洞型油气藏试井解释多参数计算和曲线准确拟合问题,导致没有准确可靠的试井解释成果,对缝洞型碳酸盐岩的开发方案制定增加了极大的难度和不确定性。
发明内容
本发明提供了一种缝洞型碳酸盐岩油气藏溶洞规模大小计算方法,其目的在于,解决现有技术中存在的上述问题,能够实现溶洞体积快速、简单、准确的计算。
本发明的技术方案如下:
本发明提供了一种缝洞型碳酸盐岩油气藏溶洞规模大小计算方法,包括下列步骤:
A、收集溶洞储层的地层、流体的基本参数以及关井压力恢复测试数据。
B、将压力恢复测试数据做成半对数horner曲线和双对数曲线。
C、在半对数horner曲线上找到直线段起始时间t1和结束时间t2,在本方法中一个直线段对应一个溶洞,如果半对数horner曲线上有几个直线段,就应该有相应的几个凹子,则代表这段储层存在多个溶洞,这时,t1和t2代表各直线段的起始时间和结束时间。
D、根据半对数horner曲线计算直线段斜率,再由直线斜率计算溶洞所在的相应储层的地层系数kh。
E、根据步骤C找到的直线段时间区间,在双对数曲线上找到对应时间区间的凹子区间。
F、根据数值积分方法计算上述凹子区间内的凹子面积S。
G、根据溶洞体积公式计算出现场溶洞体积V。
在本发明提供的实施例中,上述步骤D中,溶洞储层的地层系数kh的计算公式见式(1)。
式中:kh为溶洞储层的地层系数,mD·m;q为气井产量,104m3/d;μ为气体粘度,mPa·s;B为体积系数,无因次;m为直线段斜率,MPa/h。
在本发明提供的实施例中,上述步骤F中,凹子面积S的计算公式见式(2)。
式中:S为凹子面积,mPa·s;t1为直线段时间起始点,hr;t2为直线段时间结束点,hr;p为相应的储层井底压力,MPa。
在本发明提供的实施例中,上述步骤G中,现场溶洞体积V计算步骤如下:
第一步,先计算实验室模拟条件下溶洞当量直径D,计算公式见式(3)。
式中:D为实验室模拟条件下溶洞的当量直径,m;rw为井筒半径,m。
第二步,计算放大系数f,计算公式见式(4)。
式中:f为放大系数,无因次;Φ为裂缝孔隙度,无因次;Ct为综合压缩系数,MPa-1;C为井筒储集系数,m3/MPa;n为单位压降,MPa/h,其计算公式见式(5)。
式中:P1为t1对应的井底压力,MPa;P2为t2对应的井底压力,MPa。
第三步,计算现场规模溶洞体积,计算公式见式(6)。
V=f×4×π×(D/2)3/3 (6)
式中,V为现场规模溶洞体积,m3。
本发明的有益效果为:本发明提供的一种缝洞型碳酸盐岩油气藏溶洞规模大小计算方法经过大量的实验模拟和现场实践应用,具有较高的可行性和准确性,很好的解决了目前缝洞型碳酸盐岩油气藏在试井方面存在的多解性,理论曲线和实测曲线难以拟合以及无法准确计 算出溶洞体积大小的难题,为现场的油气田开发提供了重要的理论依据。
附图说明
为了更清楚地说明本发明实施方式的技术方案,下面将对实施方式中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1本发明实施例提供的半对数horner曲线图;
图2本发明实施例提供的双对数曲线图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
为使本发明实施方式的目的、技术方案和优点更加清楚,下面将结合本发明实施方式中的附图,对本发明实施方式中的技术方案进行清楚、完整地描述,显然,所描述的实施方式是本发明一部分实施方式,而不是全部的实施方式。因此,以下对在附图中提供的本发明的实施方式的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施方式。基于本发明中的实施方式,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施方式,都属于本发明保护的范围。
本发明提供了一种缝洞型碳酸盐岩油气藏溶洞规模大小计算方法,包括下列步骤:
A、收集溶洞储层的地层、流体的基本参数以及关井压力恢复测试数据。
B、将压力恢复测试数据做成半对数horner曲线和双对数曲线。
C、在半对数horner曲线上找到直线段起始时间t1和结束时间t2。
D、根据半对数horner曲线计算直线段斜率,再由直线斜率计算相应溶洞储层的地层系数kh。
E、根据步骤C找到的直线段时间区间,在双对数曲线上找到对应时间区间的凹子区间。
F、根据数值积分方法计算上述凹子区间内的凹子面积S。
G、根据规模溶洞体积公式计算出现场规模溶洞体积V。
在本发明提供的实施例中,上述步骤A中,储层的基本参数包括:气井产量、气体粘度、气体体积系数、井筒半径、孔隙度、综合压缩系数、井筒储集系数、气井稳定生产时间;所述关井压力恢复测试数据:时间和压力。
在本发明提供的实施例中,上述步骤B中,根据压力恢复测试时间和压力数据,将数据导入编制的软件,以时间为横坐标,压力为纵坐标,做出以横坐标为对数坐标,纵坐标为笛卡尔坐标的半对数horner曲线;再以时间为横坐标,压力和压力导数为纵坐标,做出以横坐 标为对数坐标,纵坐标为对数坐标的双对数曲线。
如图1,在本发明提供的实施例中,上述步骤C中,一个直线段对应一个溶洞,如果半对数horner曲线上有几个直线段,则代表这段储层存在多个溶洞,这时,t1和t2代表各直线段的起始时间和结束时间。则根据半对数horner曲线,找到直线段开始时间与结束时间,并且读出和记录相应的时间,该时间段即为双对数曲线凹子段开始和结束时间段。
在本发明提供的实施例中,上述步骤D中,地层系数kh的计算公式如式(1)所示:
式中:kh为溶洞储层的地层系数,mD·m;q为气井产量,MPa/m3;μ为气体粘度,mPa·s;B为体积系数,无因次;m为直线段斜率,MPa/h。
如图2,在本发明提供的实施例中,上述步骤E中,具体的分为下列几步:
第一步、将双对数曲线上凹子部分开始时间设为T1,凹子部分结束时间设为T2;
第二步、根据步骤3)找到的半对数直线段时间起始点t1与直线段时间结束点t2,由于半对数曲线横坐标为t/(t+tp),气井稳定生产时间tp已知,将t1和t2分别带入式子t1=T1/(T1+tp)和t2=T2/(T2+tp)得到T1与T2时间点;
第三步、在双对数曲线横坐标找到对应的凹子部分开始时间T1和结束时间T2,将两个时间点在曲线上用细线标定,标定后双对数曲线上阴影部分即为凹子。
在本发明提供的实施例中,上述步骤F中,凹子面积S的计算公式如式(2)所示:
式中:S为凹子面积,mPa·s;t1为直线段时间起始点,hr;t2为直线段时间结束点,hr;p为相应的储层井底压力,MPa。
在本发明提供的实施例中,上述步骤G中,现场溶洞体积V计算步骤如下:
第一步,先计算实验室模拟条件下溶洞当量直径D,计算公式如式(3)所示:
式中:D为实验室模拟条件下溶洞的当量直径,m;rw为井筒半径,m。
第二步,计算放大系数f,计算公式如式(4)所示:
式中:f为放大系数,无因次;Φ为裂缝孔隙度,无因次;Ct为综合压缩系数,MPa-1;C为井筒储集系数,m3/MPa;n为单位压降,MPa/h。
式中:P1为t1对应的井底压力,MPa;P2为t2对应的井底压力,MPa。
第三步,计算现场溶洞体积,计算公式如式(6)所示:
V=f×4×π×(D/2)3/3 (6)
式中,V为现场规模溶洞体积,m3。
实施例
如图1和图2所示,以塔中某溶洞性储层井为例,分别用本发明提供的缝洞型碳酸盐岩油气藏溶洞规模大小计算方法和目前井场上常规使用双对数曲线拟合法计算该井溶洞体积的大小,并对两种方法计算的结果与实际结果进行对比。
一、本发明提供的缝洞型碳酸盐岩油气藏溶洞规模大小计算方法具体步骤如下:
A、收集该井溶洞性储层的地层和流体的基本参数如表1所示:
表1塔中某井在溶洞性储层的地层和流体的基本参数
井筒半径(m) | 0.0889 | 气井产量(m^3/d) | 337446 |
地层气粘度(mpa.s) | 0.026 | 地层气体积系数 | 0.004 |
地层气压缩系数(1/Mpa) | 0.018 | 偏差系数 | 1.05 |
井储系数(m^3/Mpa) | 1.19 |
B、将压力恢复测试数据做成半对数horner曲线和双对数曲线。
C、在半对数horner曲线上找到直线段起始时间t1和结束时间t2,因为这段储层的半对数horner曲线上有两个直线段,证明这段储层上存在2个溶洞。
D、根据半对数horner曲线计算直线段斜率,再由直线斜率计算相应溶洞储层的地层系数kh。
E、根据步骤C找到的直线段时间区间,在双对数曲线上找到对应时间区间的凹子区间。
F、根据数值积分方法计算上述凹子区间内的凹子面积S。
G、根据规模溶洞体积公式计算出现场规模溶洞体积V。
其中,通过步骤D至步骤G对每一个溶洞体积进行计算,其计算结果如表2。
表2溶洞1和溶洞2结果及相关计算数据统计表
二、根据常规的双对数曲线拟合法计算的溶洞体积结果如表3。
表3双对数曲线拟合法计算的溶洞体积及相关参数
名称 | 数值 | 单位 |
拟合分析结果 | ||
有效渗透率(K) | 436.83 | mD |
井筒储存系数(C) | 0.3607 | m^3/MPa |
溶洞Ⅰ体积(Vv1) | 366669.75 | m^3 |
溶洞Ⅱ体积(Vv2) | 636.0928 | m^3 |
裂缝Ⅰ长度(XfⅠ) | 84.66 | m |
裂缝Ⅱ长度(XfⅡ) | 148.93 | m |
溶洞总体积(V) | 367305.84 | m^3 |
对于该溶洞性储层的油井,油井的总产量应小于总溶洞体积。目前,该井累计生产原油3.34*104m3,累计产气1.45*108m3,将地面气体体积换算到地下体积,可得目前累计生产油气体积为498360m3,根据双对数曲线拟合方法得到溶洞总体积为367305.84m3,小于当前累计生产油气体积,说明之前双对数曲线拟合算法对于该井计算不准确;根据本发明提供的凹子面积法计算得到溶洞总体积为800080.6m3,结合目前地层压力与气井废弃压力,该方法计算出的结果是可靠的。
本发明提供的缝洞型碳酸盐岩油气藏溶洞规模大小计算方法的有益效果为:本发明提供的一种缝洞型碳酸盐岩油气藏溶洞规模大小计算方法经过大量的实验模拟和现场实践应用,具有较高的可行性和准确性,很好的解决了目前缝洞型碳酸盐岩油藏试井方面存在的油气藏试井解释多参数计算曲线和理论曲线难以拟合,无法准确计算出缝洞型碳酸盐岩油气藏溶洞规模大小的难题,为现场的油气田开发提供了重要的理论依据。
以上所述,并非对本发明作任何形式上的限制,虽然本发明已以实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (4)
1.一种缝洞型碳酸盐岩油气藏溶洞规模大小计算方法,其特征在于,包括下列步骤:
A、收集溶洞性储层的地层、流体的基本参数以及关井压力恢复测试数据;
B、将压力恢复测试数据做成半对数horner曲线和双对数曲线;
C、在半对数horner曲线上找到直线段开始和结束的时间点;
D、根据半对数horner曲线计算该直线段的斜率,再由斜率计算出该直线段所在的地层系数kh;
E、根据步骤C找到的直线段时间区间,在双对数曲线上找到对应该时间区间的凹子区间;
F、根据数值积分方法计算上述凹子区间内的凹子面积S;
G、根据溶洞体积公式计算出现场溶洞体积V。
4.根据权利要求3所述的缝洞型碳酸盐岩油气藏溶洞规模大小计算方法,其特征在于,步骤G中,现场溶洞规模体积大小计算方法如下:
第一步,先计算实验室模拟条件下溶洞当量直径D,计算公式如下:
式中:D实验室模拟条件下溶洞的当量直径,m;rw为井筒半径,m;
第二步,计算放大系数f,计算公式如下:
式中:f为放大系数,无因次;Φ为裂缝孔隙度,无因次;Ct为综合压缩系数,MPa-1;C为井筒储集系数,m3/MPa;n为单位压降,MPa/h;
式中:P1为t1对应的储层井底压力,MPa;P2为t2对应的储层井底压力,MPa;第三步,计算现场溶洞规模体积大小,计算公式如下:
V=f×4×π×(D/2)3/3
式中:V为现场规模溶洞体积,m3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710159508.XA CN107577831B (zh) | 2017-03-17 | 2017-03-17 | 一种缝洞型碳酸盐岩油气藏溶洞规模大小计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710159508.XA CN107577831B (zh) | 2017-03-17 | 2017-03-17 | 一种缝洞型碳酸盐岩油气藏溶洞规模大小计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107577831A CN107577831A (zh) | 2018-01-12 |
CN107577831B true CN107577831B (zh) | 2020-08-11 |
Family
ID=61049370
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710159508.XA Active CN107577831B (zh) | 2017-03-17 | 2017-03-17 | 一种缝洞型碳酸盐岩油气藏溶洞规模大小计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107577831B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110107277B (zh) * | 2018-01-30 | 2021-11-16 | 中国石油化工股份有限公司 | 一种获取碳酸盐岩油藏中井钻遇溶洞的体积的方法 |
CN110096669B (zh) * | 2018-01-30 | 2022-09-23 | 中国石油化工股份有限公司 | 一种获取碳酸盐岩油藏中裂缝体积的方法 |
CN108386187A (zh) * | 2018-03-15 | 2018-08-10 | 中国石油化工股份有限公司 | 基于油井指示曲线分析储层结构的方法 |
CN110297280B (zh) * | 2018-03-23 | 2021-04-09 | 中国石油化工股份有限公司 | 一种碳酸盐岩超压空间分布特征的预测方法 |
CN110593860B (zh) * | 2018-05-25 | 2021-08-13 | 中国科学技术大学 | 一种参数计算方法、装置及设备 |
CN108798634B (zh) * | 2018-07-09 | 2022-04-29 | 合肥辰工科技有限公司 | 一种缝洞型储层单洞-单层均值试井解释方法 |
CN109000599B (zh) * | 2018-09-28 | 2020-03-31 | 山东大学 | 针对封闭充水溶洞的体积测量系统及方法 |
CN111950111B (zh) * | 2019-05-16 | 2024-04-05 | 中国石油化工股份有限公司 | 一种适用于底部开放的碳酸盐岩储层动态分析方法 |
CN112035993B (zh) * | 2019-05-16 | 2024-05-14 | 中国石油化工股份有限公司 | 一种底部定压的碳酸盐岩储层测试评价方法 |
CN110593865B (zh) * | 2019-09-29 | 2022-07-29 | 中国石油集团川庆钻探工程有限公司 | 一种油藏缝洞特征参数试井解释方法 |
CN110984973A (zh) * | 2019-11-29 | 2020-04-10 | 中国石油集团川庆钻探工程有限公司 | 缝洞型碳酸盐岩气藏单井控制储量的确定方法 |
CN111444462B (zh) * | 2020-04-03 | 2022-02-08 | 中国石油大学(北京) | 根据不稳定试井测算串珠体数据的方法及设备 |
CN112232631A (zh) * | 2020-09-11 | 2021-01-15 | 中国石油化工股份有限公司 | 储集体动用状况判断方法及井储关系识别方法 |
CN112267872B (zh) * | 2020-10-27 | 2021-09-14 | 成都理工大学 | 利用生产动态数据确定气井钻遇溶洞系统容积大小的方法 |
CN113027391B (zh) * | 2021-04-22 | 2022-10-28 | 中海油田服务股份有限公司 | 一种增注液波及半径的计算方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2819164C (en) * | 2010-12-16 | 2019-01-15 | Bp Corporation North America Inc. | Method of determining reservoir pressure |
CN102681013B (zh) * | 2012-04-17 | 2013-07-10 | 中国石油天然气股份有限公司 | 碳酸盐岩储层空间模型的建立方法及装置 |
CN106150477B (zh) * | 2015-04-23 | 2019-11-08 | 中国石油化工股份有限公司 | 一种确定缝洞型油藏的单井控制储量的方法 |
-
2017
- 2017-03-17 CN CN201710159508.XA patent/CN107577831B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107577831A (zh) | 2018-01-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107577831B (zh) | 一种缝洞型碳酸盐岩油气藏溶洞规模大小计算方法 | |
CN109594968B (zh) | 一种页岩气多段压裂水平井压后裂缝参数评价方法及系统 | |
RU2723769C1 (ru) | Способ расчета объема обратного потока жидкости для гидроразрыва пласта при гидроразрыве в горизонтальных скважинах в залежах газа трещиноватых песчаников | |
CN107725034B (zh) | 一种用于多级压裂水平井判别来水方向的压力监测方法 | |
CN103912269B (zh) | 页岩气储层地层破裂压力梯度录井确定方法 | |
CN104899411B (zh) | 一种储层产能预测模型建立方法和系统 | |
CN106522928B (zh) | 一种酸化压裂后停泵测井口压降不稳定试井方法 | |
CN111353205B (zh) | 用于致密气藏产水气井地层压力和动态产能的计算方法 | |
CN109242364B (zh) | 一种高温高压气井模拟井筒体积置换产能评价法 | |
CN103912248A (zh) | 水驱油田预测含水率方法 | |
CN105631184A (zh) | 一种油气井产量递减分析方法和系统 | |
CN104504604A (zh) | 一种定性气井井筒积液的方法 | |
CN105781262A (zh) | 一种产能试井方法 | |
CN113239499A (zh) | 一种页岩气井最大油嘴尺寸优选方法 | |
CN104405374A (zh) | 一种致密气藏储层应力敏感性的测量方法 | |
Czarnota et al. | Semianalytical horizontal well length optimization under pseudosteady-state conditions | |
CN112699554B (zh) | 一种基于压裂示踪约束的致密油藏水平井压后分段试井分析方法 | |
CN111950112A (zh) | 一种适用于底部封闭的碳酸盐岩储层动态分析方法 | |
CN113863920A (zh) | 一种气窜通道体积检测方法 | |
CN107725035B (zh) | 一种用于非均匀产液水平井判别来水方向的压力监测方法 | |
CN103498661B (zh) | 一种确定油藏高压物性参数的方法 | |
CN116455946B (zh) | 一种基于云端高频井口压力生产数据分析方法 | |
CN107247816B (zh) | 一种稠油油藏开展冷采化学降粘最远施工半径的判定方法 | |
CN107339096A (zh) | 一种水平井快速找水方法 | |
CN112035993A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |