Experiment of Shale Failure Mechanism Based on Particle Flow Theory
-
摘要: 为探究页岩气藏水力压裂复杂裂缝网络的形成机理,开展了基于颗粒流的页岩破裂数值试验研究。根据页岩的岩石力学性质,采用并行PSO算法反演得到页岩颗粒细观力学参数,并基于PFC离散单元模型分别进行了预制2条裂纹下的巴西劈裂和水力劈裂数值试验,通过改变试验参数分析了裂纹相互作用对其扩展的影响。数值模拟发现:巴西劈裂试验预制裂纹夹角为30°和45°时,垂直裂纹扩展至倾斜裂纹时停止扩展转而从倾斜裂纹两端开裂并扩展,变形由弹性变形变为塑性变形;裂纹夹角为90°时,垂直裂纹贯穿水平裂纹,变形主要是弹性变形;水力劈裂试验中,预制裂纹夹角越小,倾斜裂纹越容易发生剪切破坏,垂直裂纹越趋于改变扩展方向,扩展路径越复杂。研究表明,裂纹倾角对页岩气藏水力压裂缝网的形态具有重要作用,岩体应变达到峰值后轴向应力-应变呈现先减小后增大的余弦波函数关系,且倾角越小压裂缝网越复杂。Abstract: To discover the mechanism that generates complex fracture network in hydraulic fracturing of shale gas reservoirs, a numerical simulation was conducted hydraulic fracturing on shale cracking based on particle flow theory.Shale particle meso-mechanical parameters were retrieved by a parallel PSO algorithm with shale rock mechanical parameters.Brazilian tests and hydraulic fracturing tests with two preset cracks were carried out with the PFC discrete element model, with the next step that the effect of interaction between fractures on propagation was analyzed through varying test parameters.Numerical simulations demonstrated that when the preset crack was 30°or 45°inclined in the Brazilian test, the vertical crack turned into the preset crack and initiated from the sides of the preset crack.At that point, the stress-strain curve turned from elastic deformation to plastic deformation;when the preset crack was inclined 90°, the vertical crack propagated across the preset crack, and stress-strain curve demonstrated mainly elastic deformation.The smaller inclined angle of preset crack in the hydraulic fracturing test, and the easier preset crack shears, the more readily vertical crack changes its propagation direction, and further the more complicated propagation paths.The results demonstrated that the inclined angle of crack exerts significant influence on hydraulic fracture networks in shale gas reservoir;and after rock strain reaches is peak, the axial stress-strain curve decreases and then increases in a cosine wave relation.At that point, fracture networks become more complicated with decreasing inclined angle.
-
Keywords:
- particle flow /
- shale /
- preset crack /
- Brazilian test /
- hydraulic fracturing
-
浅层建井是深水钻完井作业的基础。钻井过程中,水下井口上部支撑防喷器等组件,下部悬挂各类套管串,为后续钻井提供稳定的通道[1]。海洋深水面临土质强度小、温度低、安全密度窗口小和浅层地质灾害等挑战,使深水油气井浅层建井方法的选择和作业面临巨大挑战。随着技术不断发展进步,深水作业中的浅层建井形成了多种方法,但受技术和装备所限,目前我国主要采用喷射法和钻入法。不同建井方式对海况、地质条件、作业工况等参数要求不同[2-3],深水钻井作业具有环境地质条件复杂、风险大和成本高等特点,深水浅层建井方式选择难度较大[4-5]。
国内外很多学者已经针对不同浅层建井方式的适应性开展了评价研究:杨进等人[6]分析了土质强度对不同建井方式承载力的影响,建立了不同建井方式的表层导管下入深度模型;许云锦等人[7]分析了喷射参数与土质强度的规律,指导了现场表层建井施工作业;杨进等人[8]分析了不同浅层建井方式的时效和经济性;周波等人[9]分析了不同油气井建井方法的优缺点和适应性,建立了导管下入深度计算和井口稳定性分析模型,但未形成现场方法的综合定量评价方法;B. Mackenzie等人[10]基于土体性质变化优化了深水喷射作业参数,分析了喷射法的优缺点及适应范围。但总体而言,目前只是对浅层建井方式单一维度或单一因素的适应性进行了分析,对多维度综合评价尤其是定量化评价研究还不多。因此,为科学合理地选择深水油气开发建井方式,避免经验性决策,笔者分析了不同深水浅层建井方式的工艺特点和力学特征,得到了土质强度、建井质量、作业时效、经济性和风险控制等因素对建井方式的影响,建立了单因素适应性分级;采用AHP层次分析、综合因素分析法,计算了各适应性分级权重[11- 12],依据雷达图定量分析各维度适用性,建立了以雷达图面积为评断标准的多维度深水浅层建井方式优选模型,并在多口井进行了现场作业验证。
1. 深水浅层建井方式
深水作业中采用的浅层建井技术主要有喷射法、钻入法和水下打桩法3种[13- 14]。钻入法浅层建井工艺是在井眼内下入表层导管,向环空内注水泥并固井,导管−水泥−地层土系统产生承载力,如图1所示。当海底浅部土质硬度较大时,钻入法相较于其他建井技术适应性更好。不同深度海域对钻入法施工效果影响不大。钻入下表层导管一般采用海水钻进,入泥达到设计深度后将表层导管下入并固井,此时表层导管的外表面与水泥直接接触[14]。
喷射法是采用专门设备将钻头与管柱送入井内,将两者与表层导管连接,使之形成喷射管柱。另外,钻头的尺寸相较于表层导管内径小,凭借钻头射流对井底的冲击力破碎海底土层形成井眼,表层导管凭借喷射管柱的重量下达设计深度,表层导管凭借其外表面与海底土之间的摩擦在海底支撑,如图2(a)所示[15]。喷射法下入表层导管技术有以下优点:1)可以解决钻入法井口难找的问题;2)避免因为固井的缘故压破地层,减少固井候凝时间;3)可以减少起下钻时间,节约作业成本。
水下打桩法与水上打桩法基本相同,利用打桩锤对导管的冲击使导管破开岩土下至一定深度,如图2(b)所示[16]。然而,水下打桩受水深的限制,打桩锤通过液电控制装置获取动力。与钻入法和喷射法相比,水下打桩法无需使用浮式钻井平台,只需在专用深水打桩船上便可完成,适用于小井口间距的水下基盘开发模式,经济效益良好;相较于喷射法,水下打桩法可以减少对海底土壤的扰动,从而增加导管与土层的侧面摩擦力,使表层导管承载力更高,井口更加稳定。
2. 深水浅层建井方式影响因素分析
与浅水建井相比,深水建井作业难度大、风险高[17],对作业时间窗口和海况条件要求更苛刻。深水浅层建井作业时,为降低作业风险,需考虑更多的影响因素。大量工程数据和实践表明,影响深水浅层建井方式的主要因素有海底土质强度、建井质量、作业时效、经济性和风险控制等。
2.1 土质强度
土质强度是影响浅层建井方式选择的重要因素之一[18]。喷射法和水下打桩法适用于松软地层,土质强度较高的地层易发生下入不到位的情况。钻入法适用于较硬地层,对易发生井漏、井塌等复杂情况的浅部疏松地层,钻入法下表层导管难以满足钻井作业要求。土质强度是井口承载力大小的决定因素,根据土质强度可计算出不同建井方式的井口承载力,井口承载力越大,稳定性越高。采用单桩承载力计算公式,计算得到不同建井方式的井口承载力:
Q=Qf+Qp=fAs+quAp (1) 式中:Q为浅层建井后水下井口极限承载力,kN;Qf为导管桩侧壁摩阻力,kN;Qp为导管桩端阻力,kN;As为导管桩侧壁表面积,m2;Ap为导管桩底部截面积,m2;f为导管桩侧壁单位摩擦力,kPa;qu为导管桩底部单位极限阻力,kPa。
由于导管采用空心管,一般不考虑导管桩端阻力。3种建井工艺导管侧壁单位摩擦力计算模型有所不同[13]。
2.2 建井质量
每种浅层建井方法对作业工艺和作业要求都有所不同,因此不同建井方式的建井质量存在较大差距[19],主要考虑井口稳定性、井斜角、等候时间和下入深度等指标评价建井质量 [20- 21]。井口稳定性方面,钻入法受固井质量好坏影响,在固井质量较差的情况下存在井口失稳风险;喷射法对土扰动较大,同样存在承载力不足的情况;水下打桩法承载力较高,井口稳定性风险较小。井斜角方面,钻入法受土质参数影响较小,主要通过设备参数对钻入过程进行控制,方便操作;喷射法遇较硬地层时不易控制井斜角;水下打桩法对垂直度的控制要求极高,一旦发生倾斜,难以进行回调作业,受地质参数影响较大[22]。等候时间方面,钻入法需考虑固井水泥浆的固结时间;喷射法对土壤扰动较大,需等候土壤恢复承载力;水下打桩法导管下入到设计深度的承载力较大,无需等待。下入深度方面,钻入法在深水钻井中地层稳定性较低,固井质量难以保证,因此导管下入深度需更深;喷射法对土扰动较大,同时可能受浅层气和浅水流影响导致承载力降低,需要下至一定深度;水下打桩法承载力较高,导管下入深度相比钻入法和喷射法可更浅[23]。
2.3 作业时效
深水施工作业时效性是进行方案设计与实施的关键指标。根据现场调研资料分析,钻入法受起下钻时间影响与水深有一定关系,水深500 m单井作业时间为24 h,水深1000 m单井的作业时间为27 h,水深1500 m单井作业时间为30 h。喷射法建井无固井带来的多次起下钻,500~1500 m的单井作业时间为5 h,相比钻入法作业时效提升明显。水下打桩法单井作业时间为6~7 h,其中安装专用设备需4 h,浅层建井作业时间为2~3 h,批量作业时可省去重复安装设备的时间,更适合多井批量作业。
2.4 经济性
经济效益同样是建井方式选择必须考虑的因素,与作业时间及作业装备紧密相关。通过现场调研,钻入法单井服务费为60万元,平台费用大约为120万~150万元;喷射法单井服务费为100万元,平台费用大约为50万元;水下打桩法单井服务费为130万元,平台费用大约为120万元。
2.5 风险控制
井漏、井塌、井眼掩埋、固井质量差、导管下入不到位、井口下沉和井口倾斜等对于不同浅层建井作业方式风险等级不同。海底滑坡、古河道等与海洋地质特征息息相关,对浅层建井作业方式影响不大,可通过钻前预测分析降低风险[24]。采用钻入法时,若浅部地层疏松,会存在井漏和井塌等风险;若水深较深、固井温度较低,固井质量无法保证,会导致导管承载力不足,出现井口下沉风险[25];若海况环境恶劣,会出现井眼掩埋、无法找到井眼等风险。采用喷射法时,若地层较硬,易发生表层导管下入不到位的现象,存在更换作业方式或更换井场等风险;若地层极软,下入到位后,井口稳定性无法保证,井口有下沉风险;同时,喷射作业过快,可能存在井口倾斜的风险。采用水下打桩法时,若地层较硬,表层导管下入不到位或作业控制不顺利可能会导致导管倾斜,而水下打桩法难以进行二次作业;但钻遇浅层地质灾害地层时,水下打桩法可较好地隔开地层[26- 27]。
3. 建井方式定量化优选
对深水浅层建井方式选择的主要影响因素进行适应性分级。以承载力大小为依据对海底土质强度进行适应性分级,以建井质量高低对浅层建井方式适应性分级,以作业时效高低为依据进行适应性分级,以作业费用高低为依据进行适应性分级,以风险控制的可操作性高低进行适应性分级,适应性等级均分为10级。采用综合因素分析法对不同影响因素进行权重计算,以雷达图的面积作为评价指标对深水浅层建井方式进行优选。
3.1 土质强度的适应性分级
根据井口承载力大小,建立不同土质强度的建井方式适应性分级,分为0~10级,等级越低(数值越小)表示对此土质强度的适应性越低,0表示完全不适应,10表示完全适应。根据承载力计算模型,由于作业工艺区别,相同土质强度下采用水下打桩法可获得更大的承载力。因此,以相同承载力为依据,深水浅层土质强度分级范围为0~300 kPa,基于海底土质强度的浅层建井方式的适应性分级见表1(Su代表导管下入到位后管斜处浅层土壤抗剪强度,不同土质强度对应不同的适应性等级)。考虑深水浅层多为黏性土,主要根据黏性土抗剪强度进行适应性分级。
表 1 基于土质强度的浅层建井方式适应性分级Table 1. Adaptability classification of well construction methods in shallow formations based on soil strength适应等级 钻入法/kPa 喷射法/kPa 水下打桩法/kPa 0~1 20≥Su>0 Su>300 Su>200 2~3 300≥Su>250 200≥Su>170 4~5 80≥Su>20 250≥Su>200 170≥Su>130 6~7 200≥Su>150 130≥Su>90 8~9 250≥Su>80 150≥Su>100 90≥Su>50 10 300≥Su>250 100≥Su>0 50≥Su>0 3.2 建井质量的适应性分级
以建井质量指标(井斜角、等候时间、下入深度和井口稳定性)为依据评价建井质量,根据现场建井质量作业历史数据,对不同建井方式的建井质量适应性分级,见表2。该表中数字为各类指标在不同浅层建井方式下的建井质量分级,同样分为0~10级,等级越低表示建井质量越低,4项指标的平均值即为建井质量的适应性分级。
表 2 不同建井方式的建井质量适应性分级Table 2. Adaptability classification of well construction quality for different well construction methods浅层建井方式 井斜角 等候时间 下入深度 井口稳定性 钻入法 10 5 4 7 喷射法 6 6 7 7 水下打桩法 4 10 8 8 3.3 作业时间和经济性的适应性分级
根据工艺时间对比可知,钻入法时效较低,单井作业一般在1 d以上。喷射法作业时效较高,一般约为5 h。水下打桩法由于多井作业不需起下工具,多井批钻可节省较多时间。对于作业时效的适应性分级,时效越高适应性级别越高。以浅层建井方式最小作业时间作为适应性等级的最高等级,单井平均作业时间每增加2 h,适应性级别下降1级。对于作业经济性的适应性分级,费用越低,适应性级别越高。以浅层建井方式最低费用作为适应性等级的最高等级,单井平均作业费用每增加30万元,适应性级别降低1级。
3.4 风险控制的适应性分级
根据现场作业能力对每类浅层风险在不同建井方式下的控制能力进行分级,结果见表3。其中,将井漏、固井质量差和浅层地质灾害的风险严重程度分为高、中、低;表中数值代表不同建井方式对各类风险的可控制度,数值越小表示其风险可控制度越低,各项风险可控制度的平均值即为风险控制的适应性分级。
表 3 钻井安全风险可控制度的适应性分级Table 3. Adaptability classification of drilling safety risk control degree浅层建井
方式井漏 固井质量差 浅层地质灾害 导管下入
不到位井眼掩埋 井口下沉 井口倾斜 高 中 低 高 中 低 高 中 低 钻入法 4 6 8 5 7 9 4 6 8 7 6 8 10 喷射法 10 10 10 10 10 10 5 7 9 7 10 7 7 水下打桩法 10 10 10 10 10 10 10 10 10 8 10 9 5 3.5 多维度深水浅层建井方式优选
基于以上5种因素对浅层建井方式的适应性分级评价体系,应用雷达图建立可视化多维度深水浅层建井方式优选模型,如图3所示。
模型以土质强度、建井质量、作业时效、经济性和风险控制能力作为主要考虑因素,以雷达图中5个影响因素适应性等级连线所形成的圈闭面积作为评价指标,形成的圈闭面积越大说明适应性能力越强。在海况条件、井口间距等因素影响较小的情况下,圈闭面积最大的浅层建井方法推荐为最优方法。考虑不同因素对建井影响的程度不同,设置了对应权重。采用综合因素分析法对各适应性等级进行权重计算。
A=[a11...a15⋮...⋮a51...a55] (2) 根据打分表得到不同因素的权重计算公式:
Wi=∑nj=1aij+n2−1n(n−1) (3) 式中:Wi为权重系数;aij是i因素对j因素的权重打分项,取值范围为0.1~0.9,0.5代表同样重要,分数越高相对权重越大;n为维度数,取5。
结合权重系数和适应性分级机制,得到雷达图面积:
S=sinθ2(W1W2,W2W3,W3W4,W4W5,W5W1)⋅(b1b2,b2b3,b3b4,b4b5,b5b1)T (4) 式中:θ取72°;W1~W5分别为不同因素权重系数;b1~b5分别为不同因素适应性分级级数。
4. 实例分析
以中国南海北部典型土壤环境为例进行浅层建井方式优选,目标区块作业水深为1 000~1 500 m,泥线0~100 m,土质抗剪强度为0~100 kPa,采用深水半潜式平台(工程船),完钻井深3 255 m,海底平缓,倾角较小;海底温度低,浅部土体疏松易漏失,浅层地质未有明显浅层气等地质灾害[28]。假设土质强度、建井质量、作业时效、经济性和风险控制权重系数等相同,根据多维度深水浅层建井方式优选模型,得到单口井不同维度的适应性等级(见表4)和建井方式优选雷达图(见图4)。
表 4 不同影响因素适应性分级结果Table 4. Adaptability classification results of different influencing factors建井方式 土质强度 建井质量 时效 经济性 风险可控制度 钻入法 8 6.00 2 8 6.8 喷射法 10 6.25 10 10 8.5 水下打桩法 10 7.00 7 5 8.8 根据雷达图面积,钻入法雷达图面积为18.98,喷射法雷达图面积为50.13,水下打桩法雷达图面积为36.29。对比可知,喷射法建井方式为该井最优建井方式。
作业井数影响建井方式的时效和经济性。当作业井数为10口时,由于喷射法和水下打桩法可进行批钻作业,分析作业时效和经济性,水下打桩法相比喷射法在多口井批量作业时优势更大。根据多维度深水浅层建井方式优选模型,得到10口井建井方式优选雷达图(见图5)。
根据雷达图,该区块10口井批量作业时,钻入法雷达图面积为16.75,喷射法雷达图面积为45.84,水下打桩法雷达图面积为48.37。对比之下,在装备设施满足要求的条件下,水下打桩法批量建井为最优建井方式。
由于实例分析区域为南海北部典型深水区域,可代表大部分南海北部深水区域。由此可知,当单口井或少量井作业时,采用喷射法较为合适;当批量作业10口井或以上时,采用水下打桩法较为合适。
5. 结论与建议
1)钻入法在深水浅层土质强度小于20 kPa时无法适用,建井质量较低,作业时效和作业经济性相对较差,作业安全风险相对较难控制。
2)喷射法对各因素适应性相对较好,能满足目前大部分深水井浅层建井作业,是南海北部深水区域单井或少量井的最优建井方式。
3)水下打桩法受经济性的影响较大,可批量作业,不适用于在单口井作业,但在多口井作业时设备仅需安装一次,单井平均时效和费用相比喷射法和钻入法具有明显优势。
-
[1] 姚军, 孙海, 黄朝琴, 等.页岩气藏开发中的关键力学问题[J].中国科学:物理学 力学 天文学, 2013, 43(12):1527-1547. Yao Jun, Sun Hai, Huang Zhaoqin, et al.Key mechanical problems in the development of shale gas reservoirs[J].Scientia Sinica:Physica, Mechanica Astronomica, 2013, 43(12):1527-1547. [2] 邵尚奇, 田守嶒, 李根生, 等.泥页岩地层水力裂缝延伸方位研究[J].石油钻探技术, 2014, 42(3):27-31. Shao Shangqi, Tian Shouceng, Li Gensheng, et al.Propagating orientation of hydraulic fractures in muddy shale formation[J].Petroleum Drilling Techniques, 2014, 42(3):27-31. [3] Cundall P A, Strack O D L.Particle flow code in 2 dimensions[M].Minnesota:Itasca Consulting Group Inc, 1999:1-15.
[4] Bahaaddini M, Sharrock G, Hebblewhite B K.Numerical investigation of the effect of joint geometrical parameters on the mechanical properties of a non-persistent jointed rock mass under uniaxial compression[J].Computers and Geotechnics, 2013, 49:206-225.
[5] Zhang Xiaoping, Wong Louis Ngai Yuen.Cracking processes in rock-like material containing a single flaw under uniaxial compression:a numerical study based on parallel bonded-particle model approach[J].Rock Mechanics and Rock Engineering, 2012, 45(5):711-737.
[6] 王连庆, 高谦, 王建国, 等.自然崩落采矿法的颗粒流数值模拟[J].北京科技大学学报, 2007, 29(6):557-561. Wang Lianqing, Gao Qian, Wang Jianguo, et al.Numerical simulation of natural caving method based on particle flow code in two dimensions[J].Journal of University of Science and Technology Beijing, 2007, 29(6):557-561. [7] 周健, 王家全, 曾远, 等.土坡稳定分析的颗粒流模拟[J].岩土力学, 2009, 30(1):86-90. Zhou Jian, Wang Jiaquan, Zeng Yuan, et al.Simulation of slope stability analysis by particle flow code[J].Rock and Soil Mechanics, 2009, 30(1):86-90. [8] 周健, 张刚, 孔戈.渗流的颗粒流细观模拟[J].水利学报, 2006, 37(1):28-32. Zhou Jian, Zhang Gang, Kong Ge.Meso-mechanics simulation of seepage with particle flow code[J].Journal of Hydralic Engineering, 2006, 37(1):28-32. [9] 孙锋, 张顶立, 陈铁林, 等.土体劈裂注浆过程的细观模拟研究[J].岩土工程学报, 2010, 32(3):474-480. Sun Feng, Zhang Dingli, Chen Tielin, et al.Meso-mechanical simulation of fracture grouting in soil[J].Chinese Journal of Geotechnical Engineering, 2010, 32(3):474-480. [10] 杨艳, 常晓林, 周伟, 等.裂隙岩体水力劈裂的颗粒离散元数值模拟[J].四川大学学报:工程科学版, 2012, 44(5):78-85. Yang Yan, Chang Xiaolin, Zhou Wei, et al.Simulation of hydraulic fracturing of fractured rock mass by PFC2D[J].Journal of Sichuan University:Engineering Science Edition, 2012, 44(5):78-85. [11] 杨建, 付永强, 陈鸿飞, 等.页岩储层的岩石力学特性[J].天然气工业, 2012, 32(7):12-14. Yang Jian, Fu Yongqiang, Chen Hongfei, et al.Rock mechanical characteristics of shale reservoirs[J].Natural Gas Industry, 2012, 32(7):12-14. [12] Shen B, Stephansson O, Rinne M, et al.Fracod modeling of rock fracturing and permeability change in excavation damaged zones[J].International Journal of Geomechanics, 2010, 11(4):302-313.
-
期刊类型引用(3)
1. 王思维,何淼,许明标,戴白妹. 高温高压气井测试井筒温度应力场耦合分析. 断块油气田. 2024(01): 168-176 . 百度学术
2. 戴一凡,侯冰,廖志豪. 基于相场法的深层干热岩储层水力压裂模拟研究. 石油钻探技术. 2024(02): 229-235 . 本站查看
3. 蒋振新,李军,郭勇,吴德胜,时培忠,杨宏伟,张更. 井下双梯度控压钻井井筒多相流动规律. 断块油气田. 2024(05): 936-944 . 百度学术
其他类型引用(0)
计量
- 文章访问数:
- HTML全文浏览量:
- PDF下载量:
- 被引次数: 3