CN112989721B - 一种致密油藏体积压裂水平井改造体积的快速计算方法 - Google Patents

一种致密油藏体积压裂水平井改造体积的快速计算方法 Download PDF

Info

Publication number
CN112989721B
CN112989721B CN202110242193.1A CN202110242193A CN112989721B CN 112989721 B CN112989721 B CN 112989721B CN 202110242193 A CN202110242193 A CN 202110242193A CN 112989721 B CN112989721 B CN 112989721B
Authority
CN
China
Prior art keywords
volume
yield
pressure
time
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
Application number
CN202110242193.1A
Other languages
English (en)
Other versions
CN112989721A (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.)
Northeast Petroleum University
Original Assignee
Northeast Petroleum University
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 Northeast Petroleum University filed Critical Northeast Petroleum University
Priority to CN202110242193.1A priority Critical patent/CN112989721B/zh
Publication of CN112989721A publication Critical patent/CN112989721A/zh
Application granted granted Critical
Publication of CN112989721B publication Critical patent/CN112989721B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/25Methods for stimulating production
    • E21B43/26Methods for stimulating production by forming crevices or fractures
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Fluid Mechanics (AREA)
  • General Physics & Mathematics (AREA)
  • Mining & Mineral Resources (AREA)
  • Geology (AREA)
  • Evolutionary Computation (AREA)
  • Algebra (AREA)
  • Computer Hardware Design (AREA)
  • Computing Systems (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及的是一种致密油藏体积压裂水平井改造体积的快速计算方法,其包括:步骤一、收集生产数据和计算所需的基本参数,生产数据包括压力史和流量史;基本参数包括地层的原始压力、综合压缩系数、孔隙度、流体体积系数;步骤二、绘制规整化产量与物质平衡时间的关系曲线;步骤三、利用衰竭流阶段压力变化公式,计算改造体积。本发明利用线性流阶段生产数据进行直线回归,结合了线性流流动方程来确定渗透率。本发明直接利用已有的生产动态资料绘制关系曲线,结合衰竭流阶段压力变化公式确定改造体积,无需关井,无需进行曲线拟合,计算方法简单快捷。

Description

一种致密油藏体积压裂水平井改造体积的快速计算方法
技术领域:
本发明涉及的是非常规油气开发技术领域,具体涉及的是一种致密油藏体积压裂水平井改造体积的快速计算方法。
背景技术:
致密油藏地质条件复杂,渗透率极低,非均质性强,由于常规方法开采致密油藏产能较低,因此体积压裂技术成为了提高采收率一项重要手段。通过体积压裂可使主裂缝附近生成复杂缝网,形成储层性质明显有所改善的体积改造区,极大地提高了储层流体整体渗流速度,可以有效改善油藏的渗流状况,提高油井产能。
储层改造体积是评价体积压裂改造效果的重要参数。通过确定改造体积从而实现体积压裂井改造效果的准确量化评价,可为体积压裂施工设计提供重要的理论指导。常规试井解释方法可利用压力恢复数据确定改造体积,但所需关井时间较长,影响油井产量。生产动态分析需拟合产量递减曲线,拟合过程较慢。
发明内容:
本发明的目的是提供一种致密油藏体积压裂水平井改造体积的快速计算方法,这种致密油藏体积压裂水平井改造体积的快速计算方法用于解决现有技术中确定改造体积关井时间长或者拟合过程较慢的问题。
本发明解决其技术问题所采用的技术方案是:这种致密油藏体积压裂水平井改造体积的快速计算方法:
步骤一、收集生产数据和计算所需的基本参数,生产数据包括压力史和流量史;基本参数包括地层的原始压力、综合压缩系数、孔隙度、流体体积系数;
步骤二、绘制规整化产量与物质平衡时间的关系曲线;
步骤三、利用衰竭流阶段压力变化公式,计算改造体积;
步骤3.1、根据衰竭流阶段曲线特征确定回归直线方程:
根据衰竭流阶段试井曲线斜率为-1的特征,得规整化产量与物质平衡时间关系式为:
Figure GDA0003193560930000021
式中:Δp为压差,MPa;q为产量,m3/d;tc为物质平衡时间,h;C1为直线截距;
由式(5)知回归直线斜率为-1,根据斜率确定回归直线方程;
步骤3.2、计算改造体积:
物质平衡时间下衰竭流阶段压力变化公式为:
Figure GDA0003193560930000022
式中:B为体积系数;Ct为综合压缩系数,1/MPa;pwf为井底压力,MPa;q为产量,m3/d;tc为物质平衡时间,h;Vp为孔隙体积,m3
结合式(5)和式(6),孔隙体积计算公式为:
Figure GDA0003193560930000023
因此改造体积为:
Figure GDA0003193560930000024
式中:Vs为改造体积,m3;Φ为孔隙度。
上述方案中步骤二的具体过程为:
步骤2.1、定义变量:
物质平衡时间:
Figure GDA0003193560930000025
规整化产量:
Figure GDA0003193560930000031
式中:tc为物质平衡时间,h;q为产量,m3/d;Np为累计产量,m3;pi为原始地层压力,MPa;
步骤2.2、绘制实测数据关系曲线:
根据压力史数据和产量史数据可以得到任意时刻物质平衡时间:
Figure GDA0003193560930000032
任意时刻规整化产量:
Figure GDA0003193560930000033
式中:tci为任意时刻物质平衡时间,h;Npi为任意时刻累计产量,m3;qi为任意时刻产量,m3/d;Δpi为任意时刻压差,MPa;pi为原始地层压力,MPa;pwfi为任意时刻井底压力,MPa;下标i表示任意时刻;
在双对数坐标系下,绘制规整化产量关于物质平衡时间的关系曲线。
本发明具有以下有益效果:
1、本发明提供了一种致密油藏体积压裂水平井改造体积的快速计算方法,该方法利用衰竭流阶段生产数据进行直线回归,结合了衰竭流压力变化公式来确定改造体积。该方法通过生产数据即可直接计算改造体积,无需进行曲线拟合,提高了计算速度。
2、本发明直接利用已有的生产动态资料绘制关系曲线,结合衰竭流阶段压力变化公式确定改造体积,无需关井,计算方法简单快捷。
3、本发明利用已有的生产动态资料绘制关系曲线,结合拟稳定流阶段流动方程获取改造区体积,适用性强,而致密油藏常规压力恢复试井分析方法需长时间关井,过程复杂,且试井曲线阶段不完整,难以准确拟合。
附图说明:
图1为物质平衡时间示意图;
图2为体积压裂水平井物理模型示意图;
图3为规整化产量与物质平衡时间实测曲线。
具体实施方式:
下面结合附图对本发明做进一步的说明:
这种致密油藏体积压裂水平井改造体积的快速计算方法,包括以下步骤:
步骤一、收集生产资料和计算所需的基本参数。生产资料包括生产阶段中的压力史和流量史;基本参数包括地层的原始压力、综合压缩系数、孔隙度、流体体积系数;步骤二、绘制规整化产量与物质平衡时间的关系曲线,具体包括以下步骤:
步骤2.1、定义变量:
物质平衡时间为目前累计产量与当前日产量的比值:
Figure GDA0003193560930000041
物质平衡时间几何意义如图1所示,通过物质平衡时间可建立变产量生产与定产量生产之间的等效关系。
规整化产量:
Figure GDA0003193560930000042
式中:tc为物质平衡时间,h;q为产量,m3/d;Np为累计产量,m3;pi为原始地层压力,MPa;pwf为井底流压,MPa。
步骤2.2、绘制实测数据关系曲线:
根据压力史数据和产量史数据可以得到任意时刻物质平衡时间:
Figure GDA0003193560930000051
任意时刻规整化产量:
Figure GDA0003193560930000052
式中:tci为任意时刻物质平衡时间,h;Npi为任意时刻累计产量,m3;qi为任意时刻产量,m3/d;Δpi为任意时刻压差,MPa;pi为原始地层压力,MPa;pwfi为任意时刻井底压力,MPa;下标i表示任意时刻;q1为t1时刻下产量,m3/d;tc1为t1时刻对应的物质平衡时间,h;
在双对数坐标系下,绘制规整化产量关于物质平衡时间的关系曲线。
步骤三、利用衰竭流压力变化公式,计算改造体积,具体包括以下步骤:
步骤3.1、根据衰竭流阶段曲线特征确定回归直线方程:
体积压裂水平井物理模型如图2所示,由于改造区与未改造区性质差异较大,在两区交界处形成衰竭流,与拟稳定流阶段特征类似,双对数坐标下规整化产量与物质平衡时间关系式为
Figure GDA0003193560930000053
由式(3)可知回归直线斜率为-1,根据斜率确定回归直线方程。
式中:Δp为压差,MPa;q为产量,m3/d;tc为物质平衡时间,h;C1为直线截距。
由式(5)可知回归直线斜率为-1,根据斜率确定回归直线方程。
步骤3.2、计算改造体积:
物质平衡时间下衰竭流阶段压力变化公式为
Figure GDA0003193560930000054
式中:B为体积系数;Ct为综合压缩系数,1/MPa;pwf为压力,MPa;q为产量,m3/d;tc为物质平衡时间,h;Vp为孔隙体积,m3
结合式(5)和式(6)可知孔隙体积计算公式为
Figure GDA0003193560930000061
因此改造体积为
Figure GDA0003193560930000062
式中:Vs为改造体积,m3;Φ为孔隙度。
实施例1:
(1)基本参数
收集到的基础参数见表1。
表1基础数据表
Figure GDA0003193560930000063
(2)计算结果
由图3可得到衰竭流阶段回归直线截距为4.2492,因此改造体积为5.02×106m3
压力传播达到改造区边界后发生拟稳定衰竭流动,因此可利用生产资料结合拟稳定流阶段流动方法确定改造体积。

Claims (2)

1.一种致密油藏体积压裂水平井改造体积的快速计算方法,其特征在于包括如下步骤:
步骤一、收集生产数据和计算所需的基本参数,生产数据包括压力史和流量史;基本参数包括地层的原始压力、综合压缩系数、孔隙度、流体体积系数;
步骤二、绘制规整化产量与物质平衡时间的关系曲线;
步骤三、利用衰竭流阶段压力变化公式,计算改造体积;
步骤3.1、根据衰竭流阶段曲线特征确定回归直线方程:
根据衰竭流阶段试井曲线斜率为-1的特征,得规整化产量与物质平衡时间关系式为:
Figure FDA0002962599980000011
式中:Δp为压差,MPa;q为产量,m3/d;tc为物质平衡时间,h;C1为直线截距;
由式(5)知回归直线斜率为-1,根据斜率确定回归直线方程;
步骤3.2、计算改造体积:
物质平衡时间下衰竭流阶段压力变化公式为:
Figure FDA0002962599980000012
式中:B为体积系数;Ct为综合压缩系数,1/MPa;pwf为井底压力,MPa;q为产量,m3/d;tc为物质平衡时间,h;Vp为孔隙体积,m3
结合式(5)和式(6),孔隙体积计算公式为:
Figure FDA0002962599980000013
因此改造体积为:
Figure FDA0002962599980000014
式中:Vs为改造体积,m3;φ为孔隙度。
2.根据权利要求1所述的致密油藏体积压裂水平井改造体积的快速计算方法,其特征在于:所述的步骤二的具体过程为:
步骤2.1、定义变量:
物质平衡时间:
Figure FDA0002962599980000021
规整化产量:
Figure FDA0002962599980000022
式中:tc为物质平衡时间,h;q为产量,m3/d;Np为累计产量,m3;pi为原始地层压力,MPa;
步骤2.2、绘制实测数据关系曲线:
根据压力史数据和产量史数据可以得到任意时刻物质平衡时间:
Figure FDA0002962599980000023
任意时刻规整化产量:
Figure FDA0002962599980000024
式中:tci为任意时刻物质平衡时间,h;Npi为任意时刻累计产量,m3;qi为任意时刻产量,m3/d;Δpi为任意时刻压差,MPa;pi为原始地层压力,MPa;pwfi为任意时刻井底压力,MPa;下标i表示任意时刻;
在双对数坐标系下,绘制规整化产量关于物质平衡时间的关系曲线。
CN202110242193.1A 2021-03-04 2021-03-04 一种致密油藏体积压裂水平井改造体积的快速计算方法 Active CN112989721B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110242193.1A CN112989721B (zh) 2021-03-04 2021-03-04 一种致密油藏体积压裂水平井改造体积的快速计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110242193.1A CN112989721B (zh) 2021-03-04 2021-03-04 一种致密油藏体积压裂水平井改造体积的快速计算方法

Publications (2)

Publication Number Publication Date
CN112989721A CN112989721A (zh) 2021-06-18
CN112989721B true CN112989721B (zh) 2021-09-07

Family

ID=76352890

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110242193.1A Active CN112989721B (zh) 2021-03-04 2021-03-04 一种致密油藏体积压裂水平井改造体积的快速计算方法

Country Status (1)

Country Link
CN (1) CN112989721B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113605879B (zh) * 2021-08-02 2023-03-21 中国石油大学(北京) 煤储层原始地层压力的计算方法及装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105758780A (zh) * 2016-05-04 2016-07-13 西南石油大学 一种低渗透气藏非均质复合压力衰竭程度测试方法
CN106869896A (zh) * 2017-04-26 2017-06-20 西南石油大学 一种低渗透油气藏直井体积压裂储层改造体积预测方法
CN107742020A (zh) * 2017-10-09 2018-02-27 中国石油集团川庆钻探工程有限公司地球物理勘探公司 页岩气储层压裂改造体积预测方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104504230B (zh) * 2014-10-16 2017-05-10 西南石油大学 低渗气井采收率与极限动用半径估计方法
CN106484933B (zh) * 2015-08-31 2019-06-18 中国石油化工股份有限公司 一种用于确定页岩气井井控动态储量的方法及系统
CN105447238A (zh) * 2015-11-16 2016-03-30 中石化石油工程技术服务有限公司 一种页岩气水平井压裂改造体积的估算方法
CN109522634B (zh) * 2018-11-09 2022-08-19 中国石油天然气集团有限公司 一种致密气多段体积压裂水平井数值分析方法
CN111287741B (zh) * 2020-04-12 2020-11-06 东北石油大学 一种致密油藏体积压裂改造区渗透率的快速计算方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105758780A (zh) * 2016-05-04 2016-07-13 西南石油大学 一种低渗透气藏非均质复合压力衰竭程度测试方法
CN106869896A (zh) * 2017-04-26 2017-06-20 西南石油大学 一种低渗透油气藏直井体积压裂储层改造体积预测方法
CN107742020A (zh) * 2017-10-09 2018-02-27 中国石油集团川庆钻探工程有限公司地球物理勘探公司 页岩气储层压裂改造体积预测方法

Also Published As

Publication number Publication date
CN112989721A (zh) 2021-06-18

Similar Documents

Publication Publication Date Title
CN106484933B (zh) 一种用于确定页岩气井井控动态储量的方法及系统
CN107506948B (zh) 一种基于动态泄流体积的页岩油气综合产量分析方法
CN111709847B (zh) 一种顶气边水油藏可采储量预测及开发效果评价方法
CN106204304B (zh) 一种砾岩油藏聚驱相对渗透率曲线的确定方法
CN108490156A (zh) 高温高压条件下混合气体驱油埋存量实验测试方法
CN104632158B (zh) 二氧化碳混相驱条件下油井受效阶段划分及判别方法
CN106547930A (zh) 考虑致密气藏渗流机理的泄气半径计算方法
CN111353205A (zh) 用于致密气藏产水气井地层压力和动态产能的计算方法
CN111810101B (zh) 一种水驱油藏动态分析方法及装置
CN108959767A (zh) 一种窄河道型气藏不同井型凝析油伤害数值模拟方法
CN108133086A (zh) 一种应力敏感储层中产水压裂气井裂缝半长反演方法
CN106951649B (zh) 一种测定水平井sagd蒸汽腔扩展速度的方法
CN109162701A (zh) 一种煤层裸眼井破裂压力预测方法
CN105625990A (zh) 一种体积压裂水平井开采页岩油的方法
CN112989721B (zh) 一种致密油藏体积压裂水平井改造体积的快速计算方法
CN111287740A (zh) 基于真实应变下的异常高压气藏动态储量计算方法
CN106761613B (zh) 一种二氧化碳驱替前缘的试井确定方法
CN113863920B (zh) 一种气窜通道体积检测方法
CN108133087A (zh) 一种气-水两相渗流应力敏感储层原始渗透率反演方法
CN109296363B (zh) 特低渗透油藏二氧化碳驱初期产能预测方法
CN112878987B (zh) 利用生产数据计算页岩气井控制储量的方法
CN111734394A (zh) 一种确定致密油藏压裂井不定常流井底压力的方法
CN107247816B (zh) 一种稠油油藏开展冷采化学降粘最远施工半径的判定方法
CN111287741B (zh) 一种致密油藏体积压裂改造区渗透率的快速计算方法
CN109594965A (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