CN112989721A - 一种致密油藏体积压裂水平井改造体积的快速计算方法 - Google Patents
一种致密油藏体积压裂水平井改造体积的快速计算方法 Download PDFInfo
- Publication number
- CN112989721A CN112989721A CN202110242193.1A CN202110242193A CN112989721A CN 112989721 A CN112989721 A CN 112989721A CN 202110242193 A CN202110242193 A CN 202110242193A CN 112989721 A CN112989721 A CN 112989721A
- 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.)
- Granted
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 10
- 239000000463 material Substances 0.000 claims abstract description 35
- 238000004519 manufacturing process Methods 0.000 claims abstract description 19
- 238000000034 method Methods 0.000 claims abstract description 16
- 230000008859 change Effects 0.000 claims abstract description 9
- 239000012530 fluid Substances 0.000 claims abstract description 5
- 230000015572 biosynthetic process Effects 0.000 claims description 6
- 239000011148 porous material Substances 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims description 5
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 208000016253 exhaustion Diseases 0.000 claims description 3
- 238000012986 modification Methods 0.000 abstract description 3
- 230000004048 modification Effects 0.000 abstract description 3
- 230000035699 permeability Effects 0.000 abstract description 2
- 238000012417 linear regression Methods 0.000 abstract 1
- 238000011084 recovery Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011067 equilibration Methods 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- 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
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force 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)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- General Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及的是一种致密油藏体积压裂水平井改造体积的快速计算方法,其包括:步骤一、收集生产数据和计算所需的基本参数,生产数据包括压力史和流量史;基本参数包括地层的原始压力、综合压缩系数、孔隙度、流体体积系数;步骤二、绘制规整化产量与物质平衡时间的关系曲线;步骤三、利用衰竭流阶段压力变化公式,计算改造体积。本发明利用线性流阶段生产数据进行直线回归,结合了线性流流动方程来确定渗透率。本发明直接利用已有的生产动态资料绘制关系曲线,结合衰竭流阶段压力变化公式确定改造体积,无需关井,无需进行曲线拟合,计算方法简单快捷。
Description
技术领域:
本发明涉及的是非常规油气开发技术领域,具体涉及的是一种致密油藏体积压裂水平井改造体积的快速计算方法。
背景技术:
致密油藏地质条件复杂,渗透率极低,非均质性强,由于常规方法开采致密油藏产能较低,因此体积压裂技术成为了提高采收率一项重要手段。通过体积压裂可使主裂缝附近生成复杂缝网,形成储层性质明显有所改善的体积改造区,极大地提高了储层流体整体渗流速度,可以有效改善油藏的渗流状况,提高油井产能。
储层改造体积是评价体积压裂改造效果的重要参数。通过确定改造体积从而实现体积压裂井改造效果的准确量化评价,可为体积压裂施工设计提供重要的理论指导。常规试井解释方法可利用压力恢复数据确定改造体积,但所需关井时间较长,影响油井产量。生产动态分析需拟合产量递减曲线,拟合过程较慢。
发明内容:
本发明的目的是提供一种致密油藏体积压裂水平井改造体积的快速计算方法,这种致密油藏体积压裂水平井改造体积的快速计算方法用于解决现有技术中确定改造体积关井时间长或者拟合过程较慢的问题。
本发明解决其技术问题所采用的技术方案是:这种致密油藏体积压裂水平井改造体积的快速计算方法:
步骤一、收集生产数据和计算所需的基本参数,生产数据包括压力史和流量史;基本参数包括地层的原始压力、综合压缩系数、孔隙度、流体体积系数;
步骤二、绘制规整化产量与物质平衡时间的关系曲线;
步骤三、利用衰竭流阶段压力变化公式,计算改造体积;
步骤3.1、根据衰竭流阶段曲线特征确定回归直线方程:
根据衰竭流阶段试井曲线斜率为-1的特征,得规整化产量与物质平衡时间关系式为:
式中:Δp为压差,MPa;q为产量,m3/d;tc为物质平衡时间,h;C1为直线截距;
由式(5)知回归直线斜率为-1,根据斜率确定回归直线方程;
步骤3.2、计算改造体积:
物质平衡时间下衰竭流阶段压力变化公式为:
式中:B为体积系数;Ct为综合压缩系数,1/MPa;pwf为井底压力,MPa;q为产量,m3/d;tc为物质平衡时间,h;Vp为孔隙体积,m3;
结合式(5)和式(6),孔隙体积计算公式为:
因此改造体积为:
式中:Vs为改造体积,m3;φ为孔隙度。
上述方案中步骤二的具体过程为:
步骤2.1、定义变量:
物质平衡时间:
规整化产量:
式中:tc为物质平衡时间,h;q为产量,m3/d;Np为累计产量,m3;pi为原始地层压力,MPa;
步骤2.2、绘制实测数据关系曲线:
根据压力史数据和产量史数据可以得到任意时刻物质平衡时间:
任意时刻规整化产量:
式中:tci为任意时刻物质平衡时间,h;Npi为任意时刻累计产量,m3;qi为任意时刻产量,m3/d;Δpi为任意时刻压差,MPa;pi为原始地层压力,MPa;pwfi为任意时刻井底压力,MPa;下标i表示任意时刻;
在双对数坐标系下,绘制规整化产量关于物质平衡时间的关系曲线。
本发明具有以下有益效果:
1、本发明提供了一种致密油藏体积压裂水平井改造体积的快速计算方法,该方法利用衰竭流阶段生产数据进行直线回归,结合了衰竭流压力变化公式来确定改造体积。该方法通过生产数据即可直接计算改造体积,无需进行曲线拟合,提高了计算速度。
2、本发明直接利用已有的生产动态资料绘制关系曲线,结合衰竭流阶段压力变化公式确定改造体积,无需关井,计算方法简单快捷。
3、本发明利用已有的生产动态资料绘制关系曲线,结合拟稳定流阶段流动方程获取改造区体积,适用性强,而致密油藏常规压力恢复试井分析方法需长时间关井,过程复杂,且试井曲线阶段不完整,难以准确拟合。
附图说明:
图1为物质平衡时间示意图;
图2为体积压裂水平井物理模型示意图;
图3为规整化产量与物质平衡时间实测曲线。
具体实施方式:
下面结合附图对本发明做进一步的说明:
这种致密油藏体积压裂水平井改造体积的快速计算方法,包括以下步骤:
步骤一、收集生产资料和计算所需的基本参数。生产资料包括生产阶段中的压力史和流量史;基本参数包括地层的原始压力、综合压缩系数、孔隙度、流体体积系数;步骤二、绘制规整化产量与物质平衡时间的关系曲线,具体包括以下步骤:
步骤2.1、定义变量:
物质平衡时间为目前累计产量与当前日产量的比值:
物质平衡时间几何意义如图1所示,通过物质平衡时间可建立变产量生产与定产量生产之间的等效关系。
规整化产量:
式中:tc为物质平衡时间,h;q为产量,m3/d;Np为累计产量,m3;pi为原始地层压力,MPa;pwf为井底流压,MPa。
步骤2.2、绘制实测数据关系曲线:
根据压力史数据和产量史数据可以得到任意时刻物质平衡时间:
任意时刻规整化产量:
式中:tci为任意时刻物质平衡时间,h;Npi为任意时刻累计产量,m3;qi为任意时刻产量,m3/d;Δpi为任意时刻压差,MPa;pi为原始地层压力,MPa;pwfi为任意时刻井底压力,MPa;下标i表示任意时刻。
在双对数坐标系下,绘制规整化产量关于物质平衡时间的关系曲线。
步骤三、利用衰竭流压力变化公式,计算改造体积,具体包括以下步骤:
步骤3.1、根据衰竭流阶段曲线特征确定回归直线方程:
体积压裂水平井物理模型如图2所示,由于改造区与未改造区性质差异较大,在两区交界处形成衰竭流,与拟稳定流阶段特征类似,双对数坐标下规整化产量与物质平衡时间关系式为
由式(3)可知回归直线斜率为-1,根据斜率确定回归直线方程。
式中:Δp为压差,MPa;q为产量,m3/d;tc为物质平衡时间,h;C1为直线截距。
由式(5)可知回归直线斜率为-1,根据斜率确定回归直线方程。
步骤3.2、计算改造体积:
物质平衡时间下衰竭流阶段压力变化公式为
式中:B为体积系数;Ct为综合压缩系数,1/MPa;pwf为压力,MPa;q为产量,m3/d;tc为物质平衡时间,h;Vp为孔隙体积,m3。
结合式(5)和式(6)可知孔隙体积计算公式为
因此改造体积为
式中:Vs为改造体积,m3;φ为孔隙度。
实施例1:
(1)基本参数
收集到的基础参数见表1。
表1基础数据表
(2)计算结果
由图3可得到衰竭流阶段回归直线截距为4.2492,因此改造体积为5.02×106m3。
压力传播达到改造区边界后发生拟稳定衰竭流动,因此可利用生产资料结合拟稳定流阶段流动方法确定改造体积。
Claims (2)
1.一种致密油藏体积压裂水平井改造体积的快速计算方法,其特征在于包括如下步骤:
步骤一、收集生产数据和计算所需的基本参数,生产数据包括压力史和流量史;基本参数包括地层的原始压力、综合压缩系数、孔隙度、流体体积系数;
步骤二、绘制规整化产量与物质平衡时间的关系曲线;
步骤三、利用衰竭流阶段压力变化公式,计算改造体积;
步骤3.1、根据衰竭流阶段曲线特征确定回归直线方程:
根据衰竭流阶段试井曲线斜率为-1的特征,得规整化产量与物质平衡时间关系式为:
式中:Δp为压差,MPa;q为产量,m3/d;tc为物质平衡时间,h;C1为直线截距;
由式(5)知回归直线斜率为-1,根据斜率确定回归直线方程;
步骤3.2、计算改造体积:
物质平衡时间下衰竭流阶段压力变化公式为:
式中:B为体积系数;Ct为综合压缩系数,1/MPa;pwf为井底压力,MPa;q为产量,m3/d;tc为物质平衡时间,h;Vp为孔隙体积,m3;
结合式(5)和式(6),孔隙体积计算公式为:
因此改造体积为:
式中:Vs为改造体积,m3;φ为孔隙度。
2.根据权利要求1所述的致密油藏体积压裂水平井改造体积的快速计算方法,其特征在于:所述的步骤二的具体过程为:
步骤2.1、定义变量:
物质平衡时间:
规整化产量:
式中:tc为物质平衡时间,h;q为产量,m3/d;Np为累计产量,m3;pi为原始地层压力,MPa;
步骤2.2、绘制实测数据关系曲线:
根据压力史数据和产量史数据可以得到任意时刻物质平衡时间:
任意时刻规整化产量:
式中:tci为任意时刻物质平衡时间,h;Npi为任意时刻累计产量,m3;qi为任意时刻产量,m3/d;Δpi为任意时刻压差,MPa;pi为原始地层压力,MPa;pwfi为任意时刻井底压力,MPa;下标i表示任意时刻;
在双对数坐标系下,绘制规整化产量关于物质平衡时间的关系曲线。
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 true CN112989721A (zh) | 2021-06-18 |
CN112989721B 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) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113605879A (zh) * | 2021-08-02 | 2021-11-05 | 中国石油大学(北京) | 煤储层原始地层压力的计算方法及装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104504230A (zh) * | 2014-10-16 | 2015-04-08 | 西南石油大学 | 低渗气井采收率与极限动用半径估计方法 |
CN105447238A (zh) * | 2015-11-16 | 2016-03-30 | 中石化石油工程技术服务有限公司 | 一种页岩气水平井压裂改造体积的估算方法 |
CN105758780A (zh) * | 2016-05-04 | 2016-07-13 | 西南石油大学 | 一种低渗透气藏非均质复合压力衰竭程度测试方法 |
CN106484933A (zh) * | 2015-08-31 | 2017-03-08 | 中国石油化工股份有限公司 | 一种用于确定页岩气井井控动态储量的方法及系统 |
CN106869896A (zh) * | 2017-04-26 | 2017-06-20 | 西南石油大学 | 一种低渗透油气藏直井体积压裂储层改造体积预测方法 |
CN107742020A (zh) * | 2017-10-09 | 2018-02-27 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 页岩气储层压裂改造体积预测方法 |
CN109522634A (zh) * | 2018-11-09 | 2019-03-26 | 中国石油集团川庆钻探工程有限公司 | 一种致密气多段体积压裂水平井数值分析方法 |
CN111287741A (zh) * | 2020-04-12 | 2020-06-16 | 东北石油大学 | 一种致密油藏体积压裂改造区渗透率的快速计算方法 |
-
2021
- 2021-03-04 CN CN202110242193.1A patent/CN112989721B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104504230A (zh) * | 2014-10-16 | 2015-04-08 | 西南石油大学 | 低渗气井采收率与极限动用半径估计方法 |
CN106484933A (zh) * | 2015-08-31 | 2017-03-08 | 中国石油化工股份有限公司 | 一种用于确定页岩气井井控动态储量的方法及系统 |
CN105447238A (zh) * | 2015-11-16 | 2016-03-30 | 中石化石油工程技术服务有限公司 | 一种页岩气水平井压裂改造体积的估算方法 |
CN105758780A (zh) * | 2016-05-04 | 2016-07-13 | 西南石油大学 | 一种低渗透气藏非均质复合压力衰竭程度测试方法 |
CN106869896A (zh) * | 2017-04-26 | 2017-06-20 | 西南石油大学 | 一种低渗透油气藏直井体积压裂储层改造体积预测方法 |
CN107742020A (zh) * | 2017-10-09 | 2018-02-27 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 页岩气储层压裂改造体积预测方法 |
CN109522634A (zh) * | 2018-11-09 | 2019-03-26 | 中国石油集团川庆钻探工程有限公司 | 一种致密气多段体积压裂水平井数值分析方法 |
CN111287741A (zh) * | 2020-04-12 | 2020-06-16 | 东北石油大学 | 一种致密油藏体积压裂改造区渗透率的快速计算方法 |
Non-Patent Citations (5)
Title |
---|
刘雄等: "致密油藏水平井体积压裂产能影响因素分析", 《特种油气藏》 * |
徐鹏: "致密油藏体积压裂水平井试井分析研究", 《中国优秀硕士学位论文全文数据库工程科技I辑》 * |
李宪文等: "基于物质平衡原理解释致密油藏体积压裂有效改造体积的新方法", 《科学技术与工程》 * |
王飞等: "致密气藏压裂水平井反卷积试井模型", 《石油学报》 * |
袁鸿飞等: "致密油藏多层体积压裂水平井压力动态分析", 《IFEDC》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113605879A (zh) * | 2021-08-02 | 2021-11-05 | 中国石油大学(北京) | 煤储层原始地层压力的计算方法及装置 |
CN113605879B (zh) * | 2021-08-02 | 2023-03-21 | 中国石油大学(北京) | 煤储层原始地层压力的计算方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN112989721B (zh) | 2021-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106484933B (zh) | 一种用于确定页岩气井井控动态储量的方法及系统 | |
CN107506948B (zh) | 一种基于动态泄流体积的页岩油气综合产量分析方法 | |
CN102953726B (zh) | 一种水驱油田优势通道识别方法及装置 | |
CN106481332B (zh) | 用于确定页岩气多段压裂水平井内外区动态储量的方法 | |
CN108959767B (zh) | 一种窄河道型气藏不同井型凝析油伤害数值模拟方法 | |
CN106204304B (zh) | 一种砾岩油藏聚驱相对渗透率曲线的确定方法 | |
CN111353205B (zh) | 用于致密气藏产水气井地层压力和动态产能的计算方法 | |
CN104899411B (zh) | 一种储层产能预测模型建立方法和系统 | |
CN104632158B (zh) | 二氧化碳混相驱条件下油井受效阶段划分及判别方法 | |
CN111709847A (zh) | 一种顶气边水油藏可采储量预测及开发效果评价方法 | |
CN111810101B (zh) | 一种水驱油藏动态分析方法及装置 | |
CN108133086A (zh) | 一种应力敏感储层中产水压裂气井裂缝半长反演方法 | |
CN109162701A (zh) | 一种煤层裸眼井破裂压力预测方法 | |
CN111287740A (zh) | 基于真实应变下的异常高压气藏动态储量计算方法 | |
CN111502652A (zh) | 一种三孔介质气藏水平井产量递减及生产动态预测方法 | |
CN112989721B (zh) | 一种致密油藏体积压裂水平井改造体积的快速计算方法 | |
CN107130955B (zh) | 井底流压的确定方法及储集体天然能量的确定方法 | |
CN113863920B (zh) | 一种气窜通道体积检测方法 | |
CN109296363B (zh) | 特低渗透油藏二氧化碳驱初期产能预测方法 | |
CN112878987B (zh) | 利用生产数据计算页岩气井控制储量的方法 | |
CN108664679A (zh) | 一种油气井生产数据分析方法 | |
CN111734394A (zh) | 一种确定致密油藏压裂井不定常流井底压力的方法 | |
CN109594965A (zh) | 一种不同生产压差下二氧化碳驱技术极限井距优化方法 | |
CN111287741B (zh) | 一种致密油藏体积压裂改造区渗透率的快速计算方法 | |
CN111577264A (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 |