CN107989605B - 一种径向复合型油气藏储层测试探测半径的计算方法 - Google Patents

一种径向复合型油气藏储层测试探测半径的计算方法 Download PDF

Info

Publication number
CN107989605B
CN107989605B CN201711438538.0A CN201711438538A CN107989605B CN 107989605 B CN107989605 B CN 107989605B CN 201711438538 A CN201711438538 A CN 201711438538A CN 107989605 B CN107989605 B CN 107989605B
Authority
CN
China
Prior art keywords
radius
composite
detection
time
radial
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
Application number
CN201711438538.0A
Other languages
English (en)
Other versions
CN107989605A (zh
Inventor
郑立军
徐克彬
甘惠娟
刘朝武
张灵
宫航
白田增
夏克文
王中泽
穆海林
陈艳辉
杜鑫芳
曾子轩
提云
李静宇
蒋婉雪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Petroleum Corp
CNPC Bohai Drilling Engineering Co Ltd
Original Assignee
China National Petroleum Corp
CNPC Bohai Drilling Engineering Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China National Petroleum Corp, CNPC Bohai Drilling Engineering Co Ltd filed Critical China National Petroleum Corp
Priority to CN201711438538.0A priority Critical patent/CN107989605B/zh
Publication of CN107989605A publication Critical patent/CN107989605A/zh
Application granted granted Critical
Publication of CN107989605B publication Critical patent/CN107989605B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • 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
    • 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
    • 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
    • E21B49/08—Obtaining fluid samples or testing fluids, in boreholes or wells
    • E21B49/087—Well testing, e.g. testing for reservoir productivity or formation parameters

Abstract

本发明公开了一种径向复合型油气藏储层测试探测半径的计算方法，包括以下步骤：根据现场提供的基础数据，利用径向流动阶段无因次压力导数判断复合半径探测时间；利用测试半径计算公式计算复合半径距离；将复合半径以外的区域从复合半径开始，以同心圆的方式进行逐级分割，微分时间单元及相对应的渗透率，得到在该微分时间段，该渗透率情况下的微分探测距离；将得到的距离求和得到径向复合型储层测试探测半径。本发明采用以划分时间、以不同渗透率下探测距离叠加的方法，对径向复合型油气藏储层测试探测半径进行计算，解决针对径向复合型油气藏储层测试探测半径无法计算解释的问题。用于径向复合型油气藏储层测试解释。

Description

一种径向复合型油气藏储层测试探测半径的计算方法
技术领域
本发明涉及油气田开发及油气勘探技术领域，具体涉及油气井测试资料解释技术，特别涉及径向复合型油气藏储层测试探测半径的计算方法。
背景技术
在油气井测试过程中，探测半径是一个重要的解释参数，尤其是针对研究测试影响到范围及单井控制储量的计算都是至关重要的。目前针对径向复合型油气藏，分析解释方式是能够得出复合半径，但还没有探测半径的有效计算方法。
发明内容
本发明所要解决的技术问题是，提供一种能够精确计算径向复合型油气藏储层测试探测半径的计算方法，用来合理评估径向复合型油气藏测试的探测距离。
为了解决上述技术问题，本发明采用的技术方案是：一种径向复合型油气藏储层测试探测半径的计算方法，包括以下步骤：
步骤1，根据现场提供的基础数据，利用径向流动阶段无因次压力导数判断复合半径探测时间；
步骤2，利用测试半径计算公式计算复合半径距离；
步骤3：将复合半径以外的区域从复合半径开始，以同心圆的方式进行逐级分割，微分时间单元及相对应的渗透率，得到在该微分时间段，该渗透率情况下的微分探测距离；
步骤4，将步骤1、步骤2得到的距离求和得到径向复合型储层测试探测半径。
具体地说，包括以下步骤：
步骤1，根据现场提供的基础数据,利用径向流动阶段无因次压力导数公式
Figure GDA0002702167350000011
进行判断，取P′D＝0.5的末点时间即为复合半径探测时间T1，复合半径以外的区域探测时间为T2＝T-T1，式中：P′D为径向流动阶段无因次压力导数；PD无因次压力；tD为无因次时间；CD为无因次井筒储集系数；T为测试总时间，h；T1为复合半径探测时间，h；T2为复合半径以外的区域探测时间，h；
步骤2，利用公式
Figure GDA0002702167350000021
计算复合半径距离，式中，R1为复合半径，m；Kn为复合半径以内渗透率，10-3μm2；T1为复合半径探测时间，h；φ为地层孔隙度，无因次；μ为储层流体粘度，mPa·s；Ct为地层综合压缩系数，MPa-1；
在步骤3，将T2近似相等微分成i份得到微分时间Δti，按公式
Figure GDA0002702167350000022
计算复合半径以外的区域微分单元探测距离，式中，Δri为在渗透率Ki、探测时间Δti内的探测距离，m；Ki为复合半径以外的区域微分时间对应的渗透率，10-3μm2；Δti为复合半径以外的区域微分时间，h；φ为地层孔隙度，无因次；μ为储层流体粘度，mPa·s；Ct为地层综合压缩系数，MPa-1；
步骤4，按公式
Figure GDA0002702167350000023
计算径向复合型油气藏储层测试探测半径，式中，R为径向复合型储层测试探测半径，R1为复合半径，m；Δri为在渗透率为Ki、探测时间Δti内的探测距离，m。
本发明的有益效果是：方法简单，外区微分越细密，准确度越高，能够精确计算径向复合型储层测试探测半径，同时，在实际生产工作中，采用简单的Excel表格或者简单的编写程序即可实现，易于操作，使得测试解释工作更为完备、便捷。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细说明：
本发明的径向复合型油气藏储层测试探测半径的计算方法，包括以下步骤：
步骤1，根据现场提供的基础数据，利用径向流动阶段无因次压力导数判断复合半径探测时间；
步骤2，利用测试半径计算公式计算复合半径距离；
步骤3：将复合半径以外的区域从复合半径开始，以同心圆的方式进行逐级分割，微分时间单元及相对应的渗透率，得到在该微分时间段，该渗透率情况下的微分探测距离；
步骤4，将步骤1、步骤2得到的距离求和得到径向复合型储层测试探测半径。
如图1所示，具体地说，包括以下步骤：
步骤101，根据现场提供的基础数据,利用径向流动阶段无因次压力导数公式
Figure GDA0002702167350000031
进行判断，取P′D＝0.5的末点时间即为复合半径探测时间T1，复合半径以外的区域探测时间为T2＝T-T1，式中：P′D为径向流动阶段无因次压力导数；PD无因次压力；tD为无因次时间；CD为无因次井筒储集系数；T为测试总时间，h；T1为复合半径探测时间，h；T2为复合半径以外的区域探测时间，h；
步骤102，利用公式
Figure GDA0002702167350000032
计算复合半径距离，式中，R1为复合半径，m；Kn为复合半径以内渗透率，10-3μm2；T1为复合半径探测时间，h；φ为地层孔隙度，无因次；μ为储层流体粘度，mPa·s；Ct为地层综合压缩系数，MPa-1；
在步骤103，将T2近似相等微分成i份得到微分时间Δti，按公式
Figure GDA0002702167350000033
计算复合半径以外的区域微分单元探测距离，式中，Δri为在渗透率Ki、探测时间Δti内的探测距离，m；Ki为复合半径以外的区域微分时间对应的渗透率，10-3μm2；Δti为复合半径以外的区域微分时间，h；φ为地层孔隙度，无因次；μ为储层流体粘度，mPa·s；Ct为地层综合压缩系数，MPa-1；
步骤104，按公式
Figure GDA0002702167350000034
计算径向复合型油气藏储层测试探测半径，式中，R为径向复合型储层测试探测半径，R1为复合半径，m；Δri为在渗透率为Ki、探测时间Δti内的探测距离，m。
以华北油田某地区某井径向复合型储层测试探测半径的计算为例来说明本发明的具体技术方案：
步骤1，根据现场提供的基础数据，由径向流动阶段无因此压力导数公式：
Figure GDA0002702167350000041
进行判断，取P′D＝0.5的末点时间即为复合半径探测时间T1，复合半径以外的区域探测时间为T2＝T-T1，判断结果如表1所示。
其中：P′D为径向流动阶段无因次压力导数；PD无因次压力；tD为无因次时间；CD为无因次井筒储集系数；T为测试总时间，h；T1为复合半径探测时间，h；T2为复合半径以外的区域探测时间，h。
表1探测半径时间判断结果表
Figure GDA0002702167350000042
步骤2，根据现场提供的基础数据(如表2所示)，测试半径计算公式为：
Figure GDA0002702167350000043
其中：Kn为复合半径以内渗透率，10-3μm2；T1为复合半径探测时间，h；φ为地层孔隙度，无因次；μ为储层流体粘度，mPa·s；Ct为地层综合压缩系数，MPa-1。
表2现场提供的基础数据
Figure GDA0002702167350000044
计算复合半径为0.600m。
步骤3，将复合半径以外的区域从复合半径开始，以同心圆的方式进行逐级分割，微分时间单元，及相对应的渗透率，复合半径以外的区域微分单元探测距离按公式(3)进行计，计算过程如表3所示：
Figure GDA0002702167350000051
其中：Ki为复合半径以外的区域微分时间对应的渗透率，10-3μm2；Δti为复合半径以外的区域微分时间，h；φ为地层孔隙度，无因次；μ为储层流体粘度，mPa·s；Ct为地层综合压缩系数，MPa-1。
表3微分单元数据
Figure GDA0002702167350000052
步骤4，将步骤1、步骤2的得到的距离求和得到径向复合型储层测试探测半径，计算径向复合型储层测试探测半径公式为：
Figure GDA0002702167350000053
其中：R1为复合半径，m；Δri为在渗透率为Ki、探测时间Δti内的探测距离，m。
径向复合型储层测试探测半径为6.056m。
综上所述，本发明的内容并不局限在上述的实施例中，相同领域内的有识之士可以在本发明的技术指导思想之内可以轻易提出其他的实施例，但这种实施例都包括在本发明的范围之内。

Claims (1)

1.一种径向复合型油气藏储层测试探测半径的计算方法，其特征在于，包括以下步骤：
步骤1，根据现场提供的基础数据，利用径向流动阶段无因次压力导数判断复合半径探测时间：
根据现场提供的基础数据,利用径向流动阶段无因次压力导数公式
Figure FDA0002723521270000011
进行判断，取P′D＝0.5的末点时间即为复合半径探测时间T1，复合半径以外的区域探测时间为T2＝T-T1，式中：P′D为径向流动阶段无因次压力导数；PD无因次压力；tD为无因次时间；CD为无因次井筒储集系数；T为测试总时间，h；T1为复合半径探测时间，h；T2为复合半径以外的区域探测时间，h；
步骤2，利用测试半径计算公式计算复合半径距离：
利用公式
Figure FDA0002723521270000012
计算复合半径距离，式中，R1为复合半径，m；Kn为复合半径以内渗透率，10-3μm2；T1为复合半径探测时间，h；φ为地层孔隙度，无因次；μ为储层流体粘度，mPa·s；Ct为地层综合压缩系数，MPa-1；
在步骤3，将复合半径以外的区域从复合半径开始，以同心圆的方式进行逐级分割，微分时间单元及相对应的渗透率，得到在该微分时间段，该渗透率情况下的微分探测距离：
将T2近似相等微分成i份得到微分时间Δti，按公式
Figure FDA0002723521270000013
计算复合半径以外的区域微分单元探测距离，式中，Δri为在渗透率Ki、探测时间Δti内的探测距离，m；Ki为复合半径以外的区域微分时间对应的渗透率，10-3μm2；Δti为复合半径以外的区域微分时间，h；φ为地层孔隙度，无因次；μ为储层流体粘度，mPa·s；Ct为地层综合压缩系数，MPa-1；
步骤4，按将步骤2、步骤3得到的距离求和得到径向复合型储层测试探测半径：
公式
Figure FDA0002723521270000014
计算径向复合型油气藏储层测试探测半径，式中，R为径向复合型储层测试探测半径，R1为复合半径，m；Δri为在渗透率为Ki、探测时间Δti内的探测距离，m。
CN201711438538.0A 2017-12-26 2017-12-26 一种径向复合型油气藏储层测试探测半径的计算方法 Active CN107989605B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711438538.0A CN107989605B (zh) 2017-12-26 2017-12-26 一种径向复合型油气藏储层测试探测半径的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711438538.0A CN107989605B (zh) 2017-12-26 2017-12-26 一种径向复合型油气藏储层测试探测半径的计算方法

Publications (2)

Publication Number Publication Date
CN107989605A CN107989605A (zh) 2018-05-04
CN107989605B true CN107989605B (zh) 2021-01-01

Family

ID=62041899

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711438538.0A Active CN107989605B (zh) 2017-12-26 2017-12-26 一种径向复合型油气藏储层测试探测半径的计算方法

Country Status (1)

Country Link
CN (1) CN107989605B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109555516B (zh) * 2018-10-09 2022-02-01 西南石油大学 一种径向复合多区油藏探测半径的计算方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1438404A (zh) * 2003-02-17 2003-08-27 大庆油田有限责任公司 不依赖径向流直线段的试井分析方法
WO2010132658A1 (en) * 2009-05-14 2010-11-18 Fiberod, Inc. Continuous composite rod and methods
CN103161435A (zh) * 2013-03-13 2013-06-19 中国石油大学(北京) 一种稠油热采直井试井解释方法
CN104695928A (zh) * 2015-01-27 2015-06-10 中国地质调查局油气资源调查中心 一种裂缝性致密油藏直井体积改造产能评价方法
CN105793522A (zh) * 2013-06-04 2016-07-20 沙特阿拉伯石油公司 一种用于计算径向、复合储层系统内的勘查半径的系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8851170B2 (en) * 2009-04-10 2014-10-07 Shell Oil Company Heater assisted fluid treatment of a subsurface formation

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1438404A (zh) * 2003-02-17 2003-08-27 大庆油田有限责任公司 不依赖径向流直线段的试井分析方法
WO2010132658A1 (en) * 2009-05-14 2010-11-18 Fiberod, Inc. Continuous composite rod and methods
CN103161435A (zh) * 2013-03-13 2013-06-19 中国石油大学(北京) 一种稠油热采直井试井解释方法
CN105793522A (zh) * 2013-06-04 2016-07-20 沙特阿拉伯石油公司 一种用于计算径向、复合储层系统内的勘查半径的系统
CN104695928A (zh) * 2015-01-27 2015-06-10 中国地质调查局油气资源调查中心 一种裂缝性致密油藏直井体积改造产能评价方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Interference Testing In a Two-region Composite Reservoir With a Frontal Skin;Ambastha et al;《PETSOC JOURNAL PAPER》;ONEPETRO;19940307;第33卷(第3期);第1-6页 *
径向复合油藏探测半径计算方法;崔迪声等;《油气井测试》;20050228;第14卷(第1期);第15-18页 *

Also Published As

Publication number Publication date
CN107989605A (zh) 2018-05-04

Similar Documents

Publication Publication Date Title
CN109281663B (zh) 一种确定煤层气储层渗透率的方法
CN105893679B (zh) 低产水平井续流修正试井解释方法
CN107526891B (zh) 一种聚合物驱大孔道油藏试井分析方法
CN107725034B (zh) 一种用于多级压裂水平井判别来水方向的压力监测方法
CN110984970B (zh) 一种利用地层测试确定启动压力梯度的方法
CN109441415B (zh) 基于邻井干扰的聚合物驱油藏测试井的试井解释方法
CN107944620B (zh) 一种单井稳态产能的非线性预测方法
CN108843303B (zh) 一种基于泥岩蠕变模型的油水井套管损坏预测方法
CN109242364B (zh) 一种高温高压气井模拟井筒体积置换产能评价法
CN113919240A (zh) 一种基于油气井生产数据的高压气藏储层参数的计算方法
CN104504604A (zh) 一种定性气井井筒积液的方法
CN105781262A (zh) 一种产能试井方法
CN104405374A (zh) 一种致密气藏储层应力敏感性的测量方法
CN107989605B (zh) 一种径向复合型油气藏储层测试探测半径的计算方法
CN112796725B (zh) 一种分段压裂页岩气井压裂段产气贡献率确定方法及系统
CN107605474B (zh) 一种随钻预测气层产量的方法及装置
CN113863920A (zh) 一种气窜通道体积检测方法
CN116455946B (zh) 一种基于云端高频井口压力生产数据分析方法
CN109577891B (zh) 一种深水油气井溢流监测方法
CN111764890A (zh) 一种特低渗储层的试井分析方法
WO2014051838A1 (en) Systems and methods for the determination of gas permeability
CN107725035B (zh) 一种用于非均匀产液水平井判别来水方向的压力监测方法
CN104047591B (zh) 一种利用压力定位钻井微芯片示踪器深度的方法
CN111287741B (zh) 一种致密油藏体积压裂改造区渗透率的快速计算方法
CN110965992B (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
TA01 Transfer of patent application right

Effective date of registration: 20201204

Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No.

Applicant after: CHINA NATIONAL PETROLEUM Corp.

Applicant after: CNPC Bohai Drilling Engineering Co.,Ltd.

Address before: 300457 Tianjin Economic Development Zone, No. 83, No. second main street, Tianjin, China Petroleum Building, Bohai Drilling Engineering Co., Ltd.

Applicant before: CNPC Bohai Drilling Engineering Co.,Ltd.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant