多层致密油气藏分层压裂井裂缝非均匀导流能力试井分析方法.pdf
《多层致密油气藏分层压裂井裂缝非均匀导流能力试井分析方法.pdf》由会员分享,可在线阅读,更多相关《多层致密油气藏分层压裂井裂缝非均匀导流能力试井分析方法.pdf(13页珍藏版)》请在文库网上搜索。
1、2023 年 9 月第 28 卷第 5 期中 国 石 油 勘 探CHINA PETROLEUM EXPLORATION多层致密油气藏分层压裂井裂缝非均匀导流能力试井分析方法张锦锋1 吴磊磊1 史文洋2 张城玮3 宋佳忆3 汪洋3(1 陕西延长石油(集团)有限责任公司靖边采油厂;2 常州大学;3 中国石油大学(北京)油气资源与探测国家重点实验室)摘要:压裂直井和多级压裂水平井是致密油气藏目前研究的主要井型,鲜有研究多层致密油气藏分层压裂井合采时非均匀导流能力裂缝参数解释问题。针对多层致密油气藏分层压裂工艺,摒弃致密储层传统压裂缝导流能力均匀的假设,考虑分层压裂缝扩展差异性引起的压裂缝导流能力不均
2、匀现象,通过 Laplace 空间变换、Duhamel 叠加原理和 Stehfest 数值反演方法,建立了具有非均匀导流能力压裂缝的多层致密油气藏分层压裂试井分析模型,研制了识别多层致密油气藏分层压裂井流动特征的试井分析图版,分析了井筒储集系数、裂缝表皮系数、压裂缝导流能力分布特征、储层物性等对油气渗流规律的影响。结果表明:多层致密油气藏分层压裂储层的油气渗流可分为 5 个阶段,早期流动阶段受井筒储集系数和裂缝表皮系数影响,中期流动阶段受压裂缝长度及导流能力控制;在相同生产压差下,增大各层压裂缝长度和导流能力有利于提高油气井产量;忽略压裂缝导流能力和压裂缝扩展的非均匀现象,将低估压裂缝区域压力
3、损耗、高估压裂生产井早期产能。通过实例井试井解释,获取了各层裂缝参数。关键词:致密油气藏;分层压裂;多层合试;裂缝导流能力;试井分析中图分类号:TE348 文献标识码:A DOI:10.3969/j.issn.1672-7703.2023.05.013Well testing analysis method for non-uniform fracture conductivity of multi-layer tight oil and gas reservoirs in multi-layer fracturing wellZhang Jinfeng1,Wu Leilei1,Shi Wen
4、yang2,Zhang Chengwei3,Song Jiayi3,Wang Yang3(1 Jingbian Oil Production Plant,Shaanxi Yanchang Petroleum(Group)Co.,Ltd.;2 Changzhou University;3 State Key Laboratory of Petroleum Resources and Prospecting,China University of Petroleum(Beijing)Abstract:Fracturing of vertical well and multi stage fract
5、uring of horizontal well are the two main methods for developing tight oil/gas reservoirs,but few studies have focused on fracture parameter interpretation with non-uniform conductivity of multi-layer tight oil/gas reservoirs in multi-layer fracturing and commingling production wells.Therefore,the n
6、on-uniform conductivity of hydraulic fractures caused by differential fracture propagation during multi-layer fracturing is considered rather than the traditional scenario of uniform fracture conductivity of tight reservoirs.The well testing analysis model for non-uniform fracture conductivity of mu
7、lti-layer tight oil/gas reservoirs in multi-layer fracturing well is established by using Laplace spatial transformation,Duhamel superposition principle and Stehfest numerical inversion method.In addition,the well testing type curve analyzing the flow regimes in multi-layer fractured well is researc
8、hed and the influence of several factors on fluid flow laws is determined,including the wellbore storage capacity,fracture skin factor,fracture conductivity 基金项目:国家自然科学基金项目“致密油藏注水诱发微裂缝机理及反演方法”(11872073);国家科技重大专项“丛式水平井井网优化设计与产能预测新技术”(2017ZX05009-003)。第一作者简介:张锦锋(1985-),男,陕西靖边人,本科,2009年毕业于西安石油大学,工程师,现主
9、要从事石油地质勘探方面的研究工作。地址:陕西省榆林市靖边县张家畔镇长城路南 10 号延长油田股份公司靖边采油厂,邮政编码:718500。E-mail:收稿日期:2022-12-01;修改日期:2023-05-16引用:张锦锋,吴磊磊,史文洋,等多层致密油气藏分层压裂井裂缝非均匀导流能力试井分析方法 J.中国石油勘探,2023,28(5):145-157.Zhang Jinfeng,Wu Leilei,Shi Wenyang,et al.Well testing analysis method for non-uniform fracture conductivity of multi-laye
10、r tight oil and gas reservoirs in multi-layer fracturing wellJ.China Petroleum Exploration,2023,28(5):145-157.146中 国 石 油 勘 探2023 年第 28 卷and reservoir properties.The results show that:(1)The fluid flow in multi-layer tight reservoir in multi-layer fracturing well is divided into five stages,in which
11、the early flow regime is affected by wellbore storage capacity and fracture skin factor,and the mid-term flow regime is controlled by hydraulic fracture length and conductivity;(2)Given the same production pressure difference,the larger fracture length and higher fracture conductivity are beneficial
12、 to increase well production;(3)The pressure drop in hydraulic fracture zone will be underestimated and the early well production capacity will be overestimated if fracture conductivity and non-uniform fracture propagation are ignored.Finally,a case study of well testing interpretation is given and
13、the hydraulic fracture parameters of each layer are obtained.Key words:tight oil and gas reservoir,multi-layer fracturing,commingling production,fracture conductivity,well testing analysis0 引言试井能够确定压裂井高导流能力的裂缝延伸范围,为改善压裂效果提供技术支撑。1938 年,Muskat1首次研究了垂直压裂井的压力分布规律。1973 年,Gringarten 等2给出了格林源函数和纽曼乘积方法来处理垂直
14、压裂井的非稳态流动,运用点源法推导了一系列垂直压裂井试井模型,含无限导流能力的垂直裂缝、均匀流垂直裂缝模型,这些模型奠定了压裂井试井资料解释的基础。1974 年,Gringarten 等3给出了均匀流和无限导流能力压裂井的压力求解方法。随后,Cinco-Ley 等4通过Fredholm 积分和边界元法提出了有限导流能力压裂井的压力瞬时解,基于边界元法的半解析解求解方法提高了理论压力的计算效率。Cinco-Ley 等提出的垂直裂缝模型至今仍被应用于解释有限导流裂缝参数,是目前应用最广泛的压裂井模型。我国对于压裂井的研究起步较晚,但在20世纪90年代由刘慈群、刘曰武提出的压裂井椭圆流是当时风靡一时
15、的压裂井试井分支5。21 世纪,随着压裂改造工艺和技术的发展,能较好表征压裂储层流动特征的有限导流裂缝渗流模型在我国试井分析领域得到了大力发展。2014 年,程时清等6首先建立了水平井筒多段非均匀流量试井分析模型。同年,李双等7在考虑远离井筒段裂缝的变质量线性流动及裂缝长度不等情况下,建立了考虑裂缝变质量线性流动的压裂水平井计算模型及求解方法。随后,众多学者扩展了多条压裂缝的水平井、鱼骨状多分支水平井模型8-9。何佑伟等10给出了多级压裂水平井试井分析模型,以捕捉油气沿压裂缝非均匀流入井筒的压力特性。虽然这些模型考虑了多条压裂缝之间裂缝导流能力或者流量之间的差异性,但未曾考虑裂缝导流能力沿裂缝
16、延伸方向的变化。对于裂缝导流能力不均匀的表征问题,具有四段、双段导流系数的压裂缝渗流模型被相继建立起来11-13,但该类模型假设地层为双线性流动,只在早期测试数据解释方面具有合理的精度。牟珍宝等14推导了圆形封闭油藏的变导流垂直裂缝井渗流数学模型。近年,Huang 和严谨等15-16针对大牛地气田压裂缝试井解释中裂缝长度普遍小于压裂施工设计裂缝长度的现象16,提出了压裂缝部分闭合的分段产液压裂试井模型,用裂缝导流能力不均匀分布来描述压裂缝不均匀闭合特征。Luo 等17建立了压裂铺砂部分支撑裂缝的压裂井试井分析模型,用于研究非均匀铺砂压裂缝的转向问题,该模型进一步扩展了非均匀导流能力裂缝解释方法
17、在多分支压裂井中的应用。叶义平等18考虑裂缝的变导流能力,推导了页岩油压裂水平井试井模型。王慧萍等19将存在部分闭合裂缝的垂直压裂井解释方法扩展应用到邻井生产干扰试井测试中。Li 等20将非均匀导流能力裂缝模型应用在新疆准噶尔盆地玛湖油田压裂井的压裂效果及压裂参数评价中。在数值试井方面,嵌入式离散裂缝21方法可以通过结构化网格对裂缝和储层进行离散化,后来被扩展到非结构化网格以考虑复杂的裂缝几何形状22-23。因此,具有复杂压裂缝网的压裂井数值试井模型被相应地建立起来24。这些数值模型虽然可以较好地模拟压裂缝的非均匀特征,但压裂缝分布及形态的刻画以及非结构化网格生成较为复杂,压裂缝模型建立的工作
18、量较大,限制了该方法在矿场大规模推广和应用。对于致密油气藏,多级压裂直井和多级压裂水平井是以往研究的主要井型。如苏里格致密气田,直井一般钻遇 24 个产气层,最多能钻遇 67 个气层,目前分层压裂工艺一次最多可以压裂 6 个产层。再如大牛地致密气田,多层压裂合采直井测试显示历年产气剖面变化十分明显,气井各层储量动用程度的落实工作比较困难。因此,致密气藏多层合采和分层压裂改造对气井各产层贡献程度的影响逐渐引起了研究人员的关注。目前为止,鲜有考虑多层致密储层分层压裂试井模型。即使近些年层出不穷的致密油气藏多级张锦锋等:多层致密油气藏分层压裂井裂缝非均匀导流能力试井分析方法147第 5 期压裂合采生
19、产井试井模型也多数采用均匀导流能力的平板压裂缝模型,但均没有考虑压裂过程中支撑剂未能完全填充引起的裂缝导流能力不均匀现象。为了填补这项空白,本文考虑压裂缝导流能力不均匀现象,建立了多层致密油气藏分层压裂合采试井分析模型,研究了非均匀导流能力压裂缝对压裂井的流动特征影响。通过双对数坐标下压力和压力导数的试井典型图版,判别了具有非均匀导流能力的多层致密油气藏分层压裂井的流动特征,分析了井底压力对井筒储集系数和裂缝表皮系数、地层系数比和储容系数比、裂缝长度和导流能力分布的敏感性。研究结果及认识可为多层致密油气藏分层压裂效果评价和增产措施提供理论支撑。1 多层致密油气藏分层压裂合采试井模型的建立如图
20、1a 所示,在多层致密油气藏储层中心有一口井,采用分层压裂工艺,分层压裂数为n,气井生产时多层合采,井筒半径为rw。每层水平等厚状,但各层厚度 (h)、储层物性参数(如储层渗透率K、孔隙度、综合压缩系数 ct)不同,每层流体黏度()相同。任一层压裂缝均为两翼对称垂直缝,各层裂缝半长(xf)和裂缝渗透率(Kf)均不相等,且各层压裂缝高度之和与储层厚度ht相同。如图1b所示,任一层压裂缝宽度(wf)和裂缝导流系数(FC)从井眼到压裂缝末端均是变化的。压裂缝末端裂缝较窄,具有较低的导流系数(FCb),其长度为 xpf。井筒附近裂缝较宽,具有较高的导流系数(FCa),其长度为 xf-xpf。KjjCt
21、j2rw(a)(b)层n层1层jxfjhjxpfjxpfKfwfFCaxfFCb 图 1多层致密油气藏分层压裂合采井物理模型Fig.1 Physical model of multi-layer tight oil/gas reservoir in multi-layer fracturing well(a)储层模型;(b)非均匀导流能力压裂缝模型在裂缝区域,裂缝压力分布可由 Cinco-Ley 裂缝模型进一步推导得到。在储层区域,储层压力分布可以通过无限大储层中垂直生产井的源函数解表示。主要通过裂缝压力与储层压力在裂缝储层接触面相等的数学关系,建立储层压裂缝井筒之间的压力计算模型。无量纲参数
22、定义如表 1 所示。表 1 数学模型无量纲参数定义表Table 1 Definition of dimensionless parameters in mathematical model参数定义无量纲压力D i()()1.842QBjjj jhp m pm p=无量纲流量qjqDj=Q无量纲井筒储集系数()22()CCDjhjctxf=jj无量纲距离xxD,xfj=yyDxfj=无量纲时间()a23.6 10-3()()jjjtDthhjctjxfj=裂缝长度比xfjxfj=fj148中 国 石 油 勘 探2023 年第 28 卷当储层流体为气体时,采用规整化拟时间 ta和规整化拟压力m(p
23、)完成无量纲油气渗流方程的统一25-26。下文不再区分油、气,无量纲渗流方程中统一用压力p 和时间 t 表示。1.1 分层储层压力计算根据Gringarten的研究结果,对于任一层而言,其压力分布可由无限大储层中一口定产井压力点源解表示:144jqjpjpijjctjhjt=()()d()exp)0t22144jqjxxyypjpijjctjhjtt+=()()d()exp(1)其中 Kc=t式中 生产时间,h;x 点源的 x 坐标位置,m;y 点源的 y 坐标位置,m;储层压力传导系数,m2/s。将压裂缝视为线源,通过沿压裂缝方向对公式(1)积分获得具有垂直压裂缝的储层压力解:14xfj-x
24、fjjjjqfjpjpijctht=()0t224jx xyyxt+()()expd d()(2)参数定义储层地层系数比()jKj hjKjhj=储层储容系数比jj hjctjj hjctj=()无量纲裂缝导流系数KfjwfjKj xfjFCDj=裂缝扩展比xpfjRpjxfj=注:p压力,MPa;孔隙度;hj第 j 层厚度,m;Q产量,m3/d;黏度,mPas;B体积系数,m3/m3;m(p)规整化拟压力,MPa;pi初始压力,MPa;q流量,m3/d;C井筒储集系数,m3/MPa;ct综合压缩系数,MPa-1;xf裂缝半长,m;xx 方向距离,m;yy 方向距离,m;t时间,h;ta规整化
25、拟时间,h;K渗透率,D;Kf裂缝渗透率,D;wf裂缝宽度,m;xpf远井段裂缝半长,m。下标 j 表示第 j 层,下标 D 表示无量纲形式。续表式中qfj 第 j 层裂缝流量密度,m2/s。在 Laplace 空间,由表 1 定义的无量纲参数,公式(2)中具有压裂缝的储层压力无量纲形式可以表示为 12jjjjpDjqfDjK0 xDxux=fjfjf()d(3)其中 2jjjq xqfDjq=ff式中K0 零阶修正的贝塞尔函数;uLaplace 空间变量;qfDj 在 Laplace 空间下第 j 层压裂缝区域无 量纲裂缝流量密度。1.2 单层裂缝压力计算考虑到裂缝形态的一维性,裂缝区域内部
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 多层 致密 油气藏 层压 裂缝 均匀 导流 能力 分析 方法