资源描述:
收稿日期 2018 -04 -15 基金项目 国家自然科学基金资助项目 51604062; 东北大学博士后基金资助项目 20170313; 国家科技重大专项 2017ZX05072; 安徽省科技重大专项17030901023. 作者简介 张鹏海1987 - ,男,辽宁沈阳人,东北大学讲师,博士. 第40卷第5期 2019 年 5 月 东北 大 学 学 报 自 然 科 学 版 Journal of Northeastern UniversityNatural Science Vol. 40,No. 5 May2 0 1 9 doi 10.12068/ j. issn.1005 -3026.2019.05.026 低渗储油层水力压裂裂缝延伸过程及成缝机理 张鹏海1, 张子麟2, 李 明2, 张潦源2 1. 东北大学 资源与土木工程学院, 辽宁 沈阳 110819; 2. 中国石化胜利油田分公司 石油工程技术研究院, 山东 东营 257000 摘 要 为了评价花古 6 井低渗储油层的压裂效果,在该井地表建立了微震监测系统. 基于微震监测结 果,利用裂纹延伸路径分析方法及矩张量分析方法对压裂过程中裂缝延伸的时空过程及成缝机理进行了分 析. 分析结果表明不同黏度液体注入时,裂纹的空间展布具有不同的特征,注入高黏液时储层中形成了裂缝 密度较低的缝网,注入低黏液时储层中形成的缝网具有延伸路径短、分叉明显、密度高、范围大等特征,注入中 黏液会在储层形成一些始于分支裂缝端部的裂纹;裂缝的破裂类型以剪切破裂为主占比 89. 5% ,并且破裂 类型占比不受液体黏性的影响. 关 键 词 低渗储油层;水压致裂;微震监测;裂缝延伸;成缝机理 中图分类号 TD 459 文献标志码 A 文章编号 1005 -3026201905 -0745 -05 ExtensionProcessandFractureMechanismofHydraulic Fractures in Low Permeability Reservoir ZHANG Peng-hai1, ZHANG Zi-lin2, LI Ming2, ZHANG Liao-yuan2 1. School of Resources &Civil Engineering, Northeastern University, Shenyang 110819, China; 2. Institute of Petroleum Engineering, Shengli Oil Field Company, SINOPEC, Dongying 257000, China. Corresponding author ZHANG Peng-hai, E-mail zhangpenghai2005@163. com Abstract In order to evaluate the hydraulic fracturing effect in the low permeability reservoir of Huagu 6 well, a microseismic monitoring system was established on the surface of the well. Based on microseismic data, the spatio-temporal extension process and fracture mechanism of hydraulic fractures are analyzed by means of crack extension path and moment tensor inversion. The results show that the spatial distribution characteristics of fractures changed with the liquid viscosity. The injection of high viscosity liquid formed a fracture network with low fracture density. The injection of low viscosity liquid formed a large-ranged dense fracture network with short extension path, and many crack branches. The injection of medium viscosity liquid formed some fractures that initiated from the branching fracture ends. The fracture type was mainly shear fracture89. 5% , and was not affected by the viscosity of the liquid. Key words low permeability reservoir; hydraulic fracturing; microseismic monitoring; fracture extension; fracture mechanism 水力压裂技术是一项提高低渗透油田采油率 的重要技术[1],是利用高压液体水、细沙、胶黏 剂等混合物克服地应力及岩层强度,使岩层破 裂产生裂缝[2]. 当岩层破坏时积蓄在岩层中的应变能会释放 并向外传播弹性波,这种现象叫做微地震[3]. 利用微震监测系统可以对微震波进行采集及分 析,进而获得压裂裂缝的特征,帮助评价压裂 效果[4]. 起初,由于现场监测方法的限制,对水力压裂 裂缝的研究集中在理论分析、数值模拟[5]以及室 内试验这三方面上. 随着微震监测观测装备的进 步和处理技术的发展,微震数据的质量得到提高 并促使对现场压裂裂纹发展过程、形成机制的研 万方数据 究方法不断完善[6]. Aminzadeh 等[7]通过模糊聚 类方法分析了美国 Geysers 地热田水力压裂过程 中裂缝网随时间的迁移路径;Eaton 等[8]分析了 加拿大一处地下 1 950 m 深的水平井水力压裂过 程中诱发的微震信号,根据频谱特征发现裂纹的 形成机制以张拉为主. 然而截至目前,水力压裂裂 缝形成及扩展的力学机理尚未被完全理解[9]. 因 此,如何更好地基于微震波监测数据,开展震源机 制反演、震源破裂过程分析等深层次的研究,是研 究的重点及难点. 本文结合胜利油田花古 6 井压裂过程中的微 震监测结果,利用建立的裂纹延伸路径分析方法 及矩张量分析方法分析裂缝延伸的时空过程及成 缝机理. 1 工程概况 花古 6 井位于山东省高青县花沟镇花魏村东 约 770 m,为济阳坳陷青城凸起西南部地垒高点. 该井于 2016 年 9 月 19 日开钻,同年 12 月 2 日完 井,完井垂直深度 3 167. 71 m. 储油层为常压低渗 透荧光石英砂岩. 压裂过程中采用速溶型低浓度 瓜胶压裂液体系及变黏度多尺度充填技术低黏 液造分支裂缝;中黏液增加分支裂缝宽度,加强与 主裂缝通道连通性;高黏胶液高砂比加砂,保持主 裂缝高导流通道. 为了评价压裂效果,利用 15 个全内置式地面 微地震监测台站在地面建立了“十字形”的微震 监测系统图 1. 截取 2017 年 2 月 13 日台站采 集到的波形进行滤波去噪、拾取到时、定位及计算 震源参数等数据处理,作为本文压裂裂纹时空演 化过程及产生机理分析的基础数据. 图 1 “十字形”微震监测台站布置示意图 Fig. 1 “Cruciform” microseismic monitoring stations layout 2 微震监测数据的分析 本文利用裂纹延伸路径分析方法及矩张量分 析方法对压裂过程中裂纹的产生机理及裂纹扩展 过程进行分析. 2. 1 裂纹延伸路径分析方法 按照微震事件出现的时间先后顺序及空间距 离,再现压裂裂缝的延伸路径,其过程如下①连 接射孔点及事件 1;②计算事件 2 与事件 1 的距 离,及事件 1,2 与射孔点连线的夹角,若夹角 120,则连接射孔点及事件 2,若夹角 120后不需再判断; ③计算事件 3 与事件 1,2 间的距离并连接最短距 离. 由于花古 6 井的射孔相位角为 180,故采用 120来区分初始的微震事件是由哪一侧的射孔诱 发. 以一侧射孔诱发的微震事件为例,按照以上方 法再现的压裂裂纹延伸路径示意图如图 2 所示. 图 2 压裂裂纹延伸路径示意图 Fig. 2 Schematic diagram of the fracture extension path 2. 2 矩张量分析理论及方法 将地震震源力矩模型以二阶张量的形式表 示. 根据 Ohtsu 的理论[10],经过远场近似,只考虑 P 波初动振幅,将震源简化为激冲函数模型,忽略 时间项,信号的初动振幅 Ax可以写成 Ax = CsRet,r R [r1 r2 r3] m11m12m13 m12m22m23 m13m23m33 r1 r2 r3 .1 式中Ax为传感器 P 波信号初动振幅;Cs为传 感器响应系数,其包含了传感器自身固有属性与 波导介质的材料特性;R 为震源至接收点的距离; Ret,r为考虑传感器朝向 t 与射线传播路径方 向震源与接收点连线方向r 的反射修正系数. 式1仅考虑了几何衰减. 工程尺度上,必须 考虑岩体的黏性对波的吸收衰减,需要将式1 647东北大学学报自然科学版 第 40 卷 万方数据 修正为 Ax = CsRet,re - πf vpQR R [r1 r2 r3] m11m12m13 m12m22m23 m13m23m33 r1 r2 r3 .2 式中Q 为岩体 P 波品质因子;vp为 P 波波速;f 为频率成分,当传播距离较大时,可以将 f 用波拐 角频率或主频替代. 当用于矩张量反演的微震事件传感器触发数 大于 6 时,式2可利用最小二乘法很容易计算 出震源矩张量 mpq. 要进一步解译震源力学特性, 需要对矩张量进行分解. 矩张量的分解有多种方 法,本文采用常规分解方法,将其分解为各向同性 成分ISO、双力偶成分DC与补偿线性矢量偶 极CLVD部分 m11m12m13 m21m22m23 m11m11m3 ⇒ M1 M2 M3 = 100 010 001 M1+ M2+ M3 3 →MISOZ + 100 000 00-1 M1- M3 3 →MDCX + -100 020 00-1 2M2- M1- M3 6 →MCLVDY. 3 式中,M1,M2,M3为 mpq的三个特征值,表示最大 力偶,其对应的三个本征矢量 e1,e2,e3表示最大 力偶方向,且 M1 M2 M3. MISO表示震源区域的 体积膨胀或塌缩,MDC表示震源剪切位错. 由于大 部分研究将震源破裂模型等效为双力偶,即纯剪 切模型,因此用 MDC成分所占比例 X 来判断破裂 源的剪切程度. 设 Z 为各项同性张拉成分 MISO所 占比例,Y 为补偿项 MCLVD所占比例,由式3 可得 X = 3M1- M3 - M1+14M1-7M1, Y =1 - X - Z, Z = 2M1+ M2+ M3 - M1+14M1-7M1. 4 当 X 80% 时,判断震源为剪切破坏; 当 80% ≥X≥50% 时,判断震源为介于剪切与张拉 间的混合型破坏;当 X 50% 时为张拉破坏. 同 时,破裂面法向 n 与位错方向 l 与本征向量之间 关系可以表示为 l = e1+ e3 /2, n = e1- e3 /2.} 5 纯剪切源位错情况下,n 与 l 呈垂直关系;纯 张拉时,n 与 l 呈平行关系. 图 3 不同破裂类型示意图 Fig. 3 Schematic diagram of different fracture types a张拉破裂; b剪切破裂; c混合破裂. 3 微震数据分析结果 3. 1 裂纹扩展的时空发展过程 2017 年2 月13 日1330 1520 间裂纹的动 态发展过程如图 4 所示,图中红色路径表示当前 时段内新增的裂纹延伸路径,青色路径表示初始 及当前时段之前的裂纹延伸路径. 图 4 裂纹延伸路径的动态发展过程俯视图 Fig. 4 Development of fracture extension pathplanform a1330 -1341 高黏液; b1341 -1425 低黏液; c1425 -1455 中黏液; d1455 -1520 高黏液. 第一次注入高黏液过程中,在射孔点附近形 成了一些延伸路径较短的裂缝,但在远离射孔点 的区域延伸路径则较长且裂缝网中存在大面积的 无裂纹延伸区域图 4a,由于注入高黏液的主要 目的是保持主裂缝的宽度以保证其导流能力,因 此在该过程中未形成高密度的裂缝网是完全可以 接受的. 747第 5 期 张鹏海等 低渗储油层水力压裂裂缝延伸过程及成缝机理 万方数据 在注入低黏液过程中,裂纹的延伸路径以短 程为主、具有十分明显的分叉特征,裂缝网的密度 及分布范围均显著增加,裂缝网形状趋于圆形且 裂缝网中无裂纹延伸区域面积明显变小图 4b, 这说明注入低黏液过程中形成了大量分支裂缝, 达到了较好的压裂效果. 在注入中黏液过程中,新生的裂缝数量少于 低黏液,除压裂范围北边界处存在 2 条十分明显 的向外延伸裂纹,多数裂纹的延伸路径仍为短程 图 4c,因为注入中黏液的主要目的为增加分支 裂缝宽度,所以该过程中新生的裂纹也大多始于 低黏液注入期间产生的分支裂纹端部. 由于大量分支裂缝已经形成,所以当为保持 主裂缝的宽度再次注入高黏液时,新生裂纹数量 已明显减少图 4d. 由以上分析可以发现,注入不同黏度的液体 时,裂纹的时空演化过程表现出了不同的特征. 这 既符合不同黏液的注入目的,又验证了不同黏液 的注入效果. 3. 2 微震事件震源机制的时空演化过程 2017 年2 月13 日1330 1520 期间不同黏 液注入情况下微震事件的震源机制时空演化过程 如图 5 所示. 图中圆盘的位置为定位微震事件所 在的空间位置,圆盘的大小表示微震事件破裂尺 度的大小,不同的颜色表示不同的破裂机制,蓝色 为剪切破坏,红色为张拉破坏,绿色为混合破坏. 在所有微震事件中,剪切破裂源为 256 个、占 总事件数的 89. 5% ,张拉破坏源及混合破裂源均 为 15 个、分别占总事件数的 5. 25% . 图 5 不同黏液注入情况下微震事件的震源机制演化过程 Fig. 5 Evolution of the focal mechanism of microseismic events under different viscidity liquid injection a1330 -1341 高黏液; b1341 -1425 低黏液; c1425 -1455 中黏液; d1455 -1520 高黏液. 低黏液注入期间产生微震事件 134 个,其中 剪切破裂源 122 个, 占该期间内总事件数的 91% ,略高于平均水平图 5b. 中黏液注入期间 产生微震事件 56 个,其中剪切破裂源 50 个,占该 期间内总事件数的 89. 3% ,略低于平均水平图 5c. 高黏液注入期间定位微震事件96 个,其中剪 切破裂源 84 个,占该期间内总事件数的 87. 5% , 略低于平均水平图 5a、图 5d. 通过以上数据可 以发现不同的黏液注入对微震事件的震源机制, 即水压致裂的裂纹破裂类型影响不大,可认为裂 纹破裂类型与液体的黏性无关. 张拉型裂纹占比低、剪切型裂纹占比高这一 现象与绝大多数文献及报告中记载的监测井微震 分析结果相符,但实验室尺度下水压致裂的物理 模型实验分析结果表明水压致裂产生的裂纹主要 为张拉型,两种尺度下观测到完全相反的物理现 象,其原因 1 从地应力对深部岩体破裂类型的影响来 看,三向压应力作用下岩石的破裂模式多为剪切 破裂,这已被大量的室内实验所证实并成为学者 们的共识. 一般来说,地下深部受到较高的三向压 应力作用,高压注液难以克服地应力岩体本身的 抗拉强度使深部岩体出现拉破坏. 2 从水对岩石强度的影响来看,浆液被压入 岩体的原生节理裂隙中,通过水的润滑作用可以 降低节理裂隙两侧岩石间的摩擦系数,诱发节理 裂隙两侧岩石的滑移,产生剪切破裂. 3 张拉破裂形成的物理过程主要为原先为 一连续体的岩石介质在应力作用下相互分离,基 本无摩擦现象且张拉裂纹一般不贯穿岩石颗粒, 故释放的能量较小,微震信号的幅值较低. 而剪切 破裂形成的物理过程不仅包括原先为一连续体的 岩石介质在应力作用下相互分离现象,还包括已 形成不连续面的两侧岩石介质进行相对运动,即 摩擦现象. 并且在剪切裂纹经过的路径上常见被 贯穿的岩石颗粒,故释放的能量较大,微震信号的 幅值较高. 在同等技术条件下,剪切破裂释放的微 震信号更容易被识别到. 4 从传播介质对微震波的衰减效应来看,张 拉破裂诱发的微震信号频率较高,剪切破裂诱发 的微震信号频率较低,信号频率越高在传播过程 中衰减得越快,所以当张拉破裂诱发的微震波经 过几千米的路程传播到地表传感器所在位置时, 其信号常衰减至与环境噪声相同水平,极低的信 噪比使其难以被分辨出来. 岩性相同的条件下,传 播路程直接决定了波的衰减程度,衰减效应也是 847东北大学学报自然科学版 第 40 卷 万方数据 引起现场大尺度及实验室内小尺度下观测到不同 物理现象的最主要原因. 其中第 1,2项因素为影响震源机制的本质 因素,而第 3,4项为间接因素,是通过影响微震 波的强弱来影响采集及分析结果. 基于以上分析, 可以发现以剪切破裂为主的破裂机制是符合实际 情况的,但由于较低的能量及较强的干扰导致张 拉破裂诱发的微震波难以被监测和识别到,影响 了拉伸剪切破裂占比的准确性这里暂不考虑分 析方法的影响因素. 也正是由于该原因,花古 6 井水压致裂过程中拉伸破裂的真实占比应略高于 5. 25% ,但真实占比无法通过微震监测的方法直 接获得. 4 结 论 1 花古 6 井压裂过程中不同黏度液体注入 时,裂纹的时空演化过程具有不同特征. 以保持主 裂缝宽度为目的注入高黏液会在储层形成裂纹密 度较低的缝网. 以制造分支裂缝为目的注入低黏 液会在储层形成以延伸路径短、分叉明显、密度 高、范围大为特征的缝网. 以增加分支裂缝宽度为 目的注入中黏液会在储层形成一些始于分支裂缝 端部的裂纹. 2 花古 6 井压裂裂缝的破裂类型以剪切破 裂为主占比 89. 5% ,并且破裂类型占比不受液 体黏性的影响. 剪切型裂纹占比高这一现象产生 的原因可归结为高地应力、水的润滑作用、张拉破 裂释放能量小且衰减效应明显这 4 个因素. 参考文献 [ 1 ] Economides M J,Nolte K G. Reservoir stimulation [ M]. ChichesterWiley,20001 -30. [ 2 ] Ma X, Zou Y, Li N, et al. Experimental study on the mechanism ofhydraulicfracturegrowthinaglutenite reservoir [ J]. Journal of Structural Geology, 2017, 97 37 -47. [ 3 ] Zhang P,Yang T,Yu Q,et al. Study of a seepage channel formation using the combination of microseismic monitoring technique and numerical method in Zhangmatun iron mine [J]. Rock Mechanics and Rock Engineering, 2016, 49 3699 -3708. [ 4 ] 陈海潮,唐有彩,钮凤林,等. 利用微地震参数评估水力压 裂改造效果研究进展[J]. 石油科学通报,2016,12 198 -208. Chen Hai-chao,Tang You-cai,Niu Feng-lin,et al. Recent advances in mlcrosseismic monitoring and implications for hydraulicfracturingmapping [ J ].PetroleumScience Bulletin,2016,12198 -208. [ 5 ] Park S,Kim K I,Kwon S,et al. Development of a hydraulic stimulation simulator toolbox for enhanced geothermal system design[J]. Renewable Energy,2018,118879 -895. [ 6 ] Gischig V S,Doetsch J,Maurer H,et al. On the link between stress field and small-scale hydraulic fracture growth in anisotropic rock derived from micro-seismicity [ J]. Solid Earth Discussions,2018,911 -40. [ 7 ] Aminzadeh F, TaftiTA, MaityD.Anintegrated methodology for sub-surface fracture characterization using microseismic data a case study at the NW Geysers [ J]. Computers &Geosciences,2013,545439 -49. [ 8 ] Eaton D W,Baan M V D,Birkelo B,et al. Scaling relations and spectral characteristics of tensile microseismsevidence for opening/ closing cracks during hydraulic fracturing[J]. GeophysicalJournalInternational, 2014, 196 3 1844 -1857. [ 9 ] Chen Z,WangM.Pore‐scalemodelingof hydromechanicalcoupledmechanicsinhydrofracturing process[J]. Journal of Geophysical Research Solid Earth, 2017,1222410 -3429. [10]Grosse C U, Ohtsu M. Acoustic emission testing [ M ]. BerlinSpringer,2008175 -200. 947第 5 期 张鹏海等 低渗储油层水力压裂裂缝延伸过程及成缝机理 万方数据
展开阅读全文