2. 中国石油大学(北京)石油工程学院;
3. 中国石油集团长城钻探工程有限公司;
4. 东北石油大学石油工程学院
2. College of Petroleum Engineering, China University of Petroleum(Beijing);
3. CNPC Great Wall Drilling Company;
4. School of Petroleum Engineering, Northeast Petroleum University
0 引言
页岩气以非常规天然气的形式进入人们的视野,页岩地层基质渗透率和孔隙度极低,需要用水力压裂方法改造地层,提高地层渗流能力[1-6]。然而四川长宁-威远地区的页岩气井在水力压裂过程中套管损坏现象频发。现场资料显示,截至2017年12月,某一区块共计实施52口井压裂施工,套损井28口,占比53.8%,其中的某一井区套损概率极高,为78.9%。套管损坏后,导致压裂封堵桥塞下入遇阻,不能到达预定位置,严重时甚至被迫放弃部分压裂井段,极大地影响了页岩气井的产量。
目前,页岩气水平井压裂过程中地层滑移导致剪切型套管损坏的事故大范围发生,国内外学者对其形成机理进行了多年研究,并取得了一定成果[7-17],但总体上都是采用数值解和数值模拟的方法,而通过室内试验手段研究地层滑移剪切型套损的研究很少。孙桂生[18]采用高强度膨胀剂模拟地层应力,将模拟套管放入装置中间,膨胀剂装入实验装置一侧,通过膨胀剂膨胀来模拟井筒附近地层直接剪切套管;殷凤磊[19]使用TS-1套管损害试验机对J55套管进行了直接剪切试验。这2个试验都模拟了地层剪切套管现象,但试验中的剪切力源都是直接推力,对实际地层中的剪切力来源未做解释,并且室内试验模拟有很大的简化。除此以外,在石油领域未见其他关于剪切套损室内试验装置的报道。
鉴于此,笔者针对页岩气水平井压裂过程中邻近水平段剪切套损的现象,设计了压裂过程中地层滑移室内模拟试验装置,采用该装置模拟了水力压裂过程中产生的地层不均匀隆起诱发水平层面的滑移过程,并基于非常态相似理论确定了试验模型与地层原型各参数之间的对应关系。最后在试验结果分析的基础上得出了地层变形和层间滑移的规律。研究成果对认识页岩气井剪切型套管变形及如何预防套管损坏具有一定的理论指导意义。
1 地层原型与试验模型对应关系 1.1 室内模拟试验的模拟内容室内模拟试验模拟范围有限,主要模拟范围在纵向上为水力压裂裂缝破坏范围及其上覆地层直至地面,横向上为层理断裂面延伸范围,如图 1所示。页岩气井水力压裂引发的地层膨胀变形涉及从地下目的层直至地面纵向上几千米、横向上几百米的范围,体现出的问题是涉及范围广且横纵比例悬殊。为解决以上问题,使用模拟地层材料代替实际地层并采用非常态相似方法进行模拟试验。为了体现目的层层系结构,将模型分为Ⅰ、Ⅱ两部分,上层相似模型Ⅰ采用软质材料模拟上覆地层,下层相似模型Ⅱ采用多层硬质材料模拟龙马溪组龙1d~龙1a小层。
|
| 图 1 室内模拟试验模拟范围示意图 Fig.1 Schematic diagram of simulation range of indoor experiment |
考虑到页岩气井水力压裂的特点,在整个模拟范围内可认为顶部自由,无限远处为固定边界,同时为了找到一种可调控且方便操作的替代这种能使地层抬升的作用力的方法,分析了底部注水加压、底部加膨胀剂、底部垫高度变化的刚性垫及顶部施加均匀垂向力等多种方法,最后确定模拟地层底部的作用力施加方式为底部注水加压,得出试验装置主体与模拟地层的结构示意图如图 2所示。
|
| 图 2 试验装置主体与模拟地层材料的结构示意图 Fig.2 Schematic diagram of the experimental device body and the simulated formation material |
1.2 试验模型与地层原型几何参数、力学参数对应关系
岩石普遍存在弹性模量较大而抗拉强度较小的问题。模拟试验中,在模拟地层的水平断裂面产生极小滑移量时,模拟地层即发生断裂。为了克服这一问题,选择了抗拉能力较强、弹性模量与岩石相近的高分子材料来代替岩石,作为模拟地层。综合考虑各材料弹性模量、强度以及实际加工能力,选择聚氨酯为相似模型Ⅰ的模拟地层材料,聚氯乙烯作为相似模型Ⅱ的模拟地层材料。试验模型与地层原型参数如表 1和表 2所示。
| 项目 | 弹性模量/GPa | 泊松比 | 地层直径/m | 地层厚度/m |
| 相似模型 | 0.022 6 | 0.45 | 0.4 | 0.02 |
| 实际地层 | 15.000 0 | 0.21 | 100.0 | 2 700.00 |
| 项目 | 弹性模量/GPa | 泊松比 | 地层直径/m | 地层厚度/m |
| 相似模型 | 2.04 | 0.31 | 0.4 | 0.020、0.005、 0.005、0.005 |
| 实际地层 | 25.00 | 0.05 | 100.0 | 22.100、4.900、 3.800、5.200 |
对于本试验,在给定横向、垂向尺寸相似比和模拟地层的力学参数后,即可得出其余各参数的相似比以及相应的预测系数。由于试验模型和地层原型的横纵相似比不同,根据非常态相似理论,采用畸变补偿模型[20],得出上、下两层实际地层原型与试验模型的各参数相似比如表 3和表 4所示。
| 参数 | 尺寸 | 应力 | 应变 | 弹性模量 | 泊松比 | |||||||||
| 横向D | 纵向H | 横向σr | 纵向σz | 横向εr | 纵向εz | 横向Er | 横向Ez | 横向μzr | 横向μrz | |||||
| 相似比 | 500.00 | 135 000.00 | 0.009 1 | 0.009 1 | 54 483 | 54 483 | 663.7 | 663.7 | 0.467 | 0.467 | ||||
| 参数 | 尺寸 | 应力 | 应变 | 弹性模量 | 泊松比 | |||||||||
| 横向D | 纵向H | 横向σr | 纵向σz | 横向εr | 纵向εz | 横向Er | 横向Ez | 横向μzr | 横向μrz | |||||
| 相似比 | 500.00 | 1 028.57 | 2.895 9 | 2.895 9 | 3.918 | 3.918 | 12.255 | 12.255 | 0.806 | 0.806 | ||||
相似模型Ⅰ的各预测系数为:横向弹性模量δ3=0.169 65,纵向弹性模量δ5=0.169 65,纵向泊松比为δ7=8.016 2,纵向应力δ4=0.013 8,横向泊松比δ6=0.117 55。
相似模型Ⅱ的各预测系数为:横向弹性模量δ3=0.236 3,纵向弹性模量δ5=0.236 3,纵向泊松比δ7=3.412 8,横向泊松比δ6=0.190 57,纵向应力δ4=0.05 584。
2 地层滑移室内物理模拟试验 2.1 试验装置物模试验设备主要分为底部水槽、加压装置及顶部固定装置。试验设备可固定模拟地层,使其在底部注水压力的作用下保持边界固定。顶部固定装置可采用法兰转换接头,顶部用于固定模拟地层。加工制作出的试验装置如图 3所示。
|
| 图 3 物模试验装置实物图 Fig.3 Photo of the physical model experiment device |
加压装置采用手摇式液体增压泵与平流泵组合的方式向装置内注入液体,达到增压的效果。手摇式液体增压泵体积200 mL,最大许用压力4 MPa。平流泵最大许用压力42 MPa,拥有电子调节增压面板,可免除手动加压带来的试验误差。试验中首先用手摇泵快速加压,泵内液体全部泵入设备,当压力接近试验压力时,改用平流泵缓慢精确增压,直至达到试验压力。图 4为物模试验增压装置图。
|
| 图 4 物模试验增压装置图 Fig.4 Supercharger for physical model experiment |
手摇式液体增压泵和平流泵与地层变形试验装置通过六通连接,同时在六通上安装压力表以测量内部液体压力。为了比较直观地观察到地层抬升,选用精度为0.001 mm、量程15.000 mm的千分表,并将位移测量器测量探头对准岩石中心,以实时监测地层中心垂向高度变化,达到对地层垂直位移进行精确测量的目的。测量装置如图 5所示。
|
| 图 5 地层垂直位移测量装置 Fig.5 Formation vertical displacement measuring device |
2.2 地层滑移量的测量
地层滑移量的测量原理为:在模拟地层钻出孔眼后,在观察孔眼内灌注仍处于流体状的软胶液。在软胶液未凝固前对模拟地层施加试验力,使其变形。保持该作用力恒定,等待软胶凝固成固体。待软胶完全凝固定型后,将试验力卸载,取出软胶,取出后软胶的形态就是其凝固时的形态,即地层变形后的孔眼形态。通过对软胶形态的测来测量层间滑移量。
本试验选择凝固时间为12 h的环氧树脂软胶,如图 6所示。
|
| 图 6 环氧树脂软胶 Fig.6 Epoxy resin glue |
试验中为防止软胶在孔眼内壁完全胶结、取出时发生破坏,设计了在孔眼中插入聚四氟乙烯管模拟套管的方式;另外,为了防止环氧树脂胶水在注入和试验过程中发生泄漏,将聚四氟乙烯管底部1~2 mm预先封固堵死。由于环氧树脂胶是透明的,为了便于观察,在环氧树脂胶中加入染色剂改变颜色。
为了避免孔眼附近应力集中的互相干扰,观察孔眼布局选择螺旋布局方式。模拟材料所钻孔眼均分布在螺旋线上,这样既可以使孔眼之间间隔距离相同,又可以避免孔眼之间直线距离过近。根据设计方案,绘制的观察孔眼分布图如图 7所示。图 8为直径40 cm的模拟地层材料实际钻孔位置。
|
| 图 7 观察孔眼设计方案 Fig.7 Design of the observing hole position |
|
| 图 8 实际钻孔位置 Fig.8 Actual observing hole position |
由于观察孔眼内部的变形量较小,直接通过测量工具测量会产生较大误差。为了实现精确测量井筒变形量,在得到模拟井内部变形图后,首先拍照放大,对所放大的照片进行轮廓数据提取,将提取后的数据进一步处理,完成井筒形变量的测量。轮廓处理后井筒形变如图 9所示。
|
| 图 9 轮廓处理后井筒形变图 Fig.9 Wellbore deformation diagram after contour processing |
2.3 试验结果及分析
本次页岩压裂套损室内模拟试验共进行了3组,试验参数设置如表 5所示。3组试验变形后的模拟井内部形态如图 10所示。
| 试验参数 | 试验1 | 试验2 | 试验3 |
| 压力/MPa | 0.8 | 0.8 | 0.8 |
| 垂直上抬量/cm | 1.5 | 2.5 | 2.5 |
| 上覆岩层厚度/cm | 2.0 | 2.0 | 2.0 |
| 龙1d厚度/cm | 0.5 | 0.7 | 2.0 |
| 龙1c厚度/cm | 0.5 | 0.5 | 0.5 |
| 龙1b厚度/cm | 0.5 | 0.5 | 0.5 |
| 龙1a厚度/cm | 0.5 | 0.5 | 0.5 |
| 注胶管直径/cm | 0.3 | 0.3 | 0.3 |
|
| 图 10 3组试验模拟井内部变形图 Fig.10 3 group experiments for wellbore deformation simulation |
试验后对模拟井内部泥页岩断裂面相对滑移量进行测量,得到3组不同地层厚度条件下地层滑移量数据。基于3组试验数据总体上可以发现,不同地层层间滑移量的变化趋势相同,都随着距中心点距离的增加呈先增大后减小的趋势,并且最大值集中在距中心点14~17 cm这一范围内。在实际地层中,这里的中心点就是水力压裂点。以试验2为例,在距中心点14.7 cm的第19组滑移量最大,达1.60 mm。试验2中上覆岩层-龙1d、龙1d-龙1c层间断裂面产生的滑移量分布如图 11所示。
|
| 图 11 不同层位滑移量分布图 Fig.11 Distribution of different layer slips |
同时从试验数据还可以看到,3组试验套管剪切滑移基本出现在上覆岩层-龙1d以及龙1d-龙1c这两个层间,在龙1c-龙1b少有分布,在龙1b-龙1a未出现滑移。滑移量大小也是从上至下逐渐减小。因此,不同地层埋深影响着地层滑移量的大小,即套管在纵向穿过地层发生多段套损时,其上部剪切滑移程度要大于下部剪切滑移程度。
为了解不同地层厚度下地层滑移量分布的规律,以3组试验数据为基础,绘制了不同地层厚度条件下地层滑移量与距中心点距离的关系曲线,如图 12所示。3组试验上覆地层厚度相同,试验3龙1d层厚度最大,试验1厚度最小。从试验结果上来看,试验3的平均滑移量为1.33 mm,试验1平均滑移量为0.90 mm,试验2平均滑移量为0.94 mm,试验3的滑移量明显高于另外两组试验,试验2的滑移量与试验1的滑移量接近。因此,不同的地层厚度影响着地层滑移量的大小,即当套管穿过层厚较大的页岩地层时,滑移量往往要比低厚度地层大。
|
| 图 12 相同层位不同地层厚度条件下滑移量对比曲线 Fig.12 Comparison of the slip amount of the same layer with different formation thicknesses |
3 结论
(1) 基于非常态相似理论,解决了完全相似模型试验中模型尺寸受限制的问题,并确定了岩层间滑移模拟试验的模拟地层参数。
(2) 地层滑移量主要与距地层变形中心的距离有关,即与距水力压裂点的距离有关,总体上,从地层膨胀中心,滑移量逐渐开始增大,在距中心点14~17 cm范围内达到最大值,之后开始快速下降。
(3) 不同的地层厚度以及不同的地层埋深影响着地层滑移量的大小,当套管穿过层厚较大的页岩地层时,滑移量要比低厚度地层大,当套管在纵向穿过地层发生多段套损时,其上部剪切滑移程度要大于下面部分。
| [1] |
张东晓, 杨婷云. 页岩气开发综述[J]. 石油学报, 2013, 34(4): 792-801. ZHANG D X, YANG T Y. An overview of shale-gas production[J]. Acta Petrolei Sinica, 2013, 34(4): 792-801. |
| [2] |
琚宜文, 卜红玲, 王国昌. 页岩气储层主要特征及其对储层改造的影响[J]. 地球科学进展, 2014, 29(4): 492-506. JU Y W, BU H L, WANG G C. Main characteristics of shale gas reservoirs and its effects on the reservoir reconstruction[J]. Advances in Earth Science, 2014, 29(4): 492-506. |
| [3] |
杨怀成, 毛国扬, 宋其仓, 等. 彭页HF-1井页岩气藏大型压裂工艺技术[J]. 西南石油大学学报(自然科学版), 2014, 36(5): 117-122. YANG H C, MAO G Y, SONG Q C, et al. Large scale fracturing technology of well Pengye HF-1 shale gas[J]. Journal of Southwest Petroleum University(Natural Science Edition), 2014, 36(5): 117-122. |
| [4] |
石晓闪, 刘大安, 崔振东, 等. 页岩气开采压裂技术分析与思考[J]. 天然气勘探与开发, 2015, 38(3): 62-65, 69. SHI X S, LIU D A, CUI Z D, et al. Analysis and thinking on shale gas mining fracturing technology[J]. Natural Gas Exploration and Development, 2015, 38(3): 62-65, 69. DOI:10.3969/j.issn.1673-3177.2015.03.014 |
| [5] |
赵金洲, 任岚, 沈骋, 等. 页岩气储层缝网压裂理论与技术研究新进展[J]. 天然气工业, 2018, 38(3): 1-14. ZHAO J Z, REN L, SHEN C, et al. Latest research progresses in network fracturing theories and technologies for shale gas reservoirs[J]. Natural Gas Industry, 2018, 38(3): 1-14. |
| [6] |
尹丛彬. 页岩压裂裂缝渗透率的测试与分析[J]. 天然气工业, 2018, 38(3): 60-68. YIN C B. Test and analysis on the permeability of fractured fractures in shale reservoirs[J]. Natural Gas Industry, 2018, 38(3): 60-68. |
| [7] |
SCHWALL G H, DENNEY C A. Subsidence induced casing deformation mechanisms in the Ekofisk field[R]. SPE 28091-MS, 1994.
|
| [8] |
DUSSEAULT M B, BRUNO M S, BARRERA J. Casing shear: causes, cases, cures[R]. SPE 72060-PA, 2001.
|
| [9] |
LI X Q, ZHU T G, FANG L T, et al. Effect of sand production on casing integrity[R]. SPE PETSOC-2006-195, 2006.
|
| [10] |
李军, 李玉梅, 张德龙, 等. 页岩气井分段压裂套损影响因素分析[J]. 断块油气田, 2017, 24(3): 387-390. LI J, LI Y M, ZHANG D L, et al. Analysis of casing damage for staged fracturing in shale gas wells[J]. Fault-Block Oil & Gas Field, 2017, 24(3): 388-390. |
| [11] |
WONG R C K, CHAU K T. Casing impairment induced by shear slip along a weak layer in shale due to fluid (steam) injection[J]. Journal of Canadian Petroleum Technology, 2006, 45(12): 60-66. |
| [12] |
WANG T, JIANG R Y, HAN Y L, et al. The casing failure mechanism in the untra-deep horizontal well and research on the technology of the precaution, diagnose and treatment in H Oilfield[J]. Journal of Southwest Petroleum University, 2009, 31(1): 156-161. |
| [13] |
黄小兰, 刘建军, 杨春和. 注水油田剪切套损机理研究与数值模拟[J]. 武汉工业学院学报, 2009, 28(3): 94-97. HUANG X L, LIU J J, YANG C H. Mechanism analysis and numerical simulation of casing's shear failure in water injection oilfield[J]. Journal of Wuhan Polytechnic University, 2009, 28(3): 94-97. DOI:10.3969/j.issn.1009-4881.2009.03.023 |
| [14] |
HAN X H, KHAN S, ANSARI S A, et al. Prediction of injection induced formation shear[R]. SPE 151840-MS, 2012.
|
| [15] |
艾池, 胡超洋, 崔月明. 延缓大庆油田标准层套损的套管优选[J]. 石油钻探技术, 2015, 43(6): 7-11. AI C, HU C Y, CUI Y M. Casing optimization for delaying casing damage in the datum bed of Daqing Oilfield[J]. Petroleum Drilling Techniques, 2015, 43(6): 7-11. |
| [16] |
田中兰, 石林, 乔磊. 页岩气水平井井筒完整性问题及对策[J]. 天然气工业, 2015, 35(9): 70-76. TAIN Z L, SHI L, QIAO L. Reseach of and countermeasure for wellbore integrity of shale gas horizontal well[J]. Natural Gas Industry, 2015, 35(9): 70-76. DOI:10.3787/j.issn.1000-0976.2015.09.010 |
| [17] |
高利军, 乔磊, 柳占立, 等. 页岩储层剪切套损的数值模拟及固井对策研究[J]. 石油机械, 2016, 44(10): 6-10, 16. GAO L J, QIAO L, LIU Z L, et al. Numerical modeling and cementing countermeasure analysis of casing shear demage in shale reservoir[J]. China Petroleum Machinery, 2016, 44(10): 6-10, 16. |
| [18] |
孙桂生.杏北油田套管损坏预测与处理对策研究[D].大庆: 大庆石油学院, 2008. SUN G S. Study on casing damage prediction and treatment countermeasures in Xingbei Oilfield[D]. Daqing: Daqing Petroleum Institute, 2008. http://cdmd.cnki.com.cn/Article/CDMD-10220-2008207656.htm |
| [19] |
殷凤磊.基于光纤传感技术的套管应变及管外压力的在线监测技术[D].哈尔滨: 哈尔滨工业大学, 2010. YIN F L. On-line monitoring technology for casing strain and external pressure based on optical fiber sensing technology[D]. Harbin: Harbin Institute of Technology, 2010. http://cdmd.cnki.com.cn/Article/CDMD-10213-1011260472.htm |
| [20] |
牛恩宽.畸变模型理论在滑坡模型试验中的应用研究[D].宜昌: 三峡大学, 2007. NIU E K. Study on application for distortion and compensated method in landslide model test[D]. Yichang: Sanxia University, 2007. http://cdmd.cnki.com.cn/Article/CDMD-11934-2007110302.htm |


