Mechanism of Deep Shaft Surrounding Rock Failure and Shaft Wall Structure Design
-
摘要: 为指导赵固二矿井筒建设,综合理论分析、数值模拟等方法,对井筒围岩损伤变形及井壁结构进行研究,结果表明:井筒穿越地层水平围压系数为1.75~1.86,导致软弱地层井筒塑性区蝶化分布,塑性区尺寸与双向主应力关系的RPP特性围岩强度影响较大;当埋深500~600 m时,在冻结前后,及时浇筑井筒围岩损伤形态尺寸差异较大。基于此,确定选择冻结法施工,理论计算设计并校核了井壁结构参数,现场施工监测显示设计参数能够满足工程建设的安全稳定,为类似条件井筒建设提供参考。Abstract: To provide effective guidance for shaft construction, theoretical analysis and numerical simulation were used to study the damage and deformation of the surrounding rock and the structure of the shaft wall. The results show that the horizontal confining pressure coefficient for shaft crossing geological formations is 1.75 to 1.86, leading to a butterfly-shaped distribution of the plastic zone in the weak strata. The RPP characteristics of the relationship between plastic zone size and bidirectional principal stress have a significant impact on the strength of the surrounding rock. The size differences of the damage forms of the wellbore surrounding rock before and after freezing and timely pouring are significant at a burial depth of 500~600 meters. The freezing method was selected for construction before and after freezing based on the vector displacement of the surrounding rock. Based on this, the theoretical calculation and design of the shaft wall structure parameters were verified, and on-site construction monitoring shows that the design parameters can meet the safety and stability requirements of the project construction. The research results have practical significance for the construction of similar shafts under similar conditions.
-
Key words:
- deep shaft /
- surrounding rock damage /
- freezing construction /
- shaft structure
-
0. 引言
浅部资源的枯竭促使煤炭资源开发向深部进军,井筒作为深部开采的咽喉,其稳定性事关矿山安全、正常运营[1-3] 。何满潮[4]从深井复杂地质力学环境造成井筒失稳的角度对深部建井力学研究进展进行了全面论述;刘环宇等[5]通过水文与地质调查,结合相关破坏理论,获得了兖州矿区立井井筒的破坏机制;赵海军等[6]分析了竖井围岩稳定性影响因素,通过数值模拟和地质力学分析研究了竖井变形破坏机理;张 好等[7]基于恒源煤矿地质特征及弹性力学理论,给出了表现井壁破坏过程的应力状态的解析解;郤保平等[8]针对热干岩地质构造,基于黏弹性理论,推导了三维应力下井筒围岩损伤破坏判别式;李 康等[9]研究了邻近设施偏压和流固耦合作用下井筒的破坏机理;王胜开等[10]通过数值模拟划分了深部盲井施工安全风险等级。
深部井筒揭露不同地层,其围岩属性及力学环境也各不相同,实质上,井筒的变形破坏主要是由力学不平衡状态所引起的[11]。本文以赵固二矿西回风立井为背景,基于地质力学条件分析了围岩损伤机理,确定了施工工法,设计了井壁结构及参数,对类似矿井具有指导价值。
1. 井筒工程概况
赵固二矿为解决矿井深部开采的通风问题,决定在矿井西翼建设西回风立井,该立井井筒净直径6.0 m,深度902 m,赵固二矿西风井是我国第四个穿过700 m冲积层的井筒。
井筒依次穿越第四系(10.05 m)、新近系(694.55 m)、二叠系下统石盒子组(196.09 m)、二叠系下统山西组(1.27 m),其中第四系和新近系多为黏土和砂质黏土,石盒子组和山西组围岩多裂隙发育,根据土层和岩石的物理力学试验可知:土层黏聚力为16.6~214.1 kPa,含水率5.3%~30.6%,孔隙比0.369~0.834;二叠系地层岩石总体抗压强度为8.7~60.1 MPa,岩石坚固性等级属不坚固—坚固岩石。根据承压完整井计算公式预测新近系和第四系含水层、基岩风化带裂隙承压含水层、二叠系煤层顶板孔隙裂隙承压含水层的预测井筒涌水量分别为1074.45 m3/h、1.67 m3/h、236.90 m3/h。
如此大埋深的井筒,加之其特殊的地质条件,现既有工程的参考性严重不足,对于适用于厚松散层的冻结法施工,其对本区条件下围岩的损伤机理如何仍有待探索。
2. 井筒破坏机理
井筒穿越地层多、地应力差异大,康红普等[12]统计了我国煤矿地应力分布的一般规律:
$$ \sigma_{{\mathrm{v}}} = 0.0245H $$ (1) $$ \sigma_{{\mathrm{H}}} = 0.0215H + 3.267 $$ (2) $$ \sigma_{{\mathrm{h}}} = 0.0113H + 1.954 $$ (3) 式中:H为埋深,m;σv、σH、σh分别为垂直应力、最大水平主应力、最小水平主应力,MPa。
将上式进行分段计算,获得不同埋深条件下的应力分布基本规律如表1所示(井筒围压比η=σH/σh)。由表1可见,井筒断面始终处于非均匀应力状态下,且埋深100 m至埋深800 m时逐渐增加并稳定于1.86。根据非均匀应力场孔洞围岩蝶形塑性区破坏理论[13- 14],围岩自均匀应力场至非均匀应力场变化期间围岩塑性区将呈现圆形—椭圆形—蝶形过渡趋势,在蝶形状态下围岩塑性区将处于RPP(塑性区半径R与最大主应力P1和最小主应力P3的关系)应力敏感状态,极其容易诱发围岩失稳,此外,围岩塑性区的RPP特征曲线具有围岩强度抵抗性,即围岩强度越大,在相同应力状态下围岩塑性损伤越小,如图1所示。
表 1 地应力分布规律H/m σH/MPa σh/MPa σv/MPa η 100 5.42 3.08 2.45 1.75 200 7.57 4.21 4.90 1.79 300 9.72 5.34 7.35 1.81 400 11.87 6.47 9.80 1.83 500 14.02 7.60 12.25 1.84 600 16.17 8.73 14.70 1.85 700 18.32 9.86 17.15 1.86 800 20.47 10.99 19.60 1.86 研究结果[15]表明:井筒冻结除能够达到阻水目的外,还能够显著提升围岩力学强度。基于前文分析,采用FLAC3D模拟软件开展数值分析,以埋深550~600 m为例,采用应变软化(Strain-Softening)破坏准则,分别模拟该埋深条件下冻结前后(见表2)围岩塑性破坏与矢量位移,为说明关键问题,对井壁一次性支护1000 mm的C80混凝土,且文中数值模拟为一次性开挖。冻结前后井筒围岩塑性区分布如图2、图3所示。
表 2 冻结前后围岩力学属性状态 岩石力学参数 体积模量/GPa 剪切模量/MPa 内摩擦角/(°) 黏聚力/MPa 抗拉强度/MPa 冻结前 0.11 0.12 20 0.81 2.57 冻结后 5.08 2.47 30 3.02 0.57 非冻结时沿最小水平应力(σh)方向井身围岩塑性区尺寸为6.9 m,沿最大水平应力方向(σH)围岩塑性区尺寸为3.5 m,井筒断面围岩塑性区类似于蝶形,如图2所示;冻结后的沿最小水平应力(σh)方向井身围岩塑性区尺寸为4.3 m, 沿最大水平应力方向(σH)围岩塑性区尺寸为1.9 m,井筒断面围岩塑性区类似于椭圆形—蝶形的过渡态,如图3所示。
冻结前后井筒整体及断面围岩的矢量位移由图4—图6所示,结合塑性区的分布特征,着重分析沿最小水平应力(σh)方向井身围岩矢量位移,如图4所示,冻结前后最大矢量位移均集中于井筒施工迎头位置,主要体现为未冻结时底鼓量可达146.3 mm,冻结时底鼓量为81.7 mm;冻结前后井身混凝土支护后变形量相对稳定,其中未冻结时明显数值更大,截取井身断面获得的井筒断面矢量位移云图如图5、图6所示。图5为开挖未浇筑时的井筒断面矢量位移,未冻结时断面位移量最大为120.2 mm,而冻结时井筒断面最大位移为59.5 mm,二者差距2倍;及时混凝土浇筑后非冻结时断面位移量最大为61.4 mm,大于50.0 mm的变形区域范围广,而冻结时井筒断面最大位移为44.0 mm,大于40.0 mm的变形区域较小。
考虑到现浇混凝土初凝及其达到设计强度的时间效应、实际穿层围岩强度差异,数值模拟的结果将比实际情况偏小。基于上述分析,选择冻结法施工,除能够阻断涌水,还能够在井壁浇筑混凝土稳定前保持围岩强度,以弱化围岩径向变形。根据冻结井筒深度设计的相关规范要求[16],并考虑本井筒揭露地层及壁座等情况,确定冻结深度为783 m,冻结深度以下井筒采用普通法施工。
3. 井壁结构设计
由于赵固二矿西回风立井主要是冲积层结构,采用冻结法施工,结合我国冻结井筒绝大多数井壁结构形式,采用钢筋混凝土双层井壁结构形式[17]。为保障外层井壁接茬、低温等因素在解冻后的破裂变形,在内外井壁间增加聚氯乙烯塑料夹层,外层井壁与井帮之间敷设聚苯乙烯泡沫塑料板。
外层井壁承受施工阶段冻胀力,内层井壁承受永久水压力,内、外层井壁整体共同承受永久地压。本次井筒按水土压力校核,考虑摩擦力和不均匀载荷作用时:
$$ P_{{\mathrm{JS}}}=0.01kH $$ (4) $$ P_{{\mathrm{DJ}}}=0.01~0.012H $$ (5) $$ P_{{\mathrm{ST}}}=0.013H $$ (6) 式中:PJS为静水压力,MPa;PDJ为冻结压力,MPa;PST为水土压力,MPa;H为计算处深度。
井壁厚度按式(7)计算:
$$ t={r}_{\text{n}}\left(\sqrt{\frac{{f}_{\rm{c}}+{\rho }_{\mathrm{min}}{f}_{\rm{y}}'}{{f}_{\rm{c}}+{\rho }_{\mathrm{min}}{f}_{\rm{y}}'-2{v}_{\text{k}}{P}_{\text{k}}}}-1\right) $$ (7) 式中:t为井壁厚度,m;rn为计算处井壁内半径,m; fc为混凝土轴心抗压强度设计值,N/mm2;ρmin为最小配筋率;fy′ 为吊挂钢筋强度设计值,N/mm2;vk为钢筋混凝土受压构件载荷系数;Pk为井筒冻结压力,MPa。
外层井壁吊挂力按式(8)—式(10)计算:
$$ {\gamma _{\text{0}}}{N_{\text{d}}} \leqslant {f_{\text{y}}}{{\text{A}}_{\text{z}}}{\text{ }} $$ (8) $$ {N_{{\text{d}},{\text{k}}}} = {\text{π}} {\gamma _{\text{h}}}{h_{\text{d}}}(R_{{\text{ww}}}^2 - R_{{\text{wn}}}^2) $$ (9) $$ {N_{\text{d}}} = {\gamma _{{\text{kh}}}}{N_{{\text{d}},{\text{k}}}} $$ (10) 式中:Nd为井壁吊挂力的计算值,MN;Nd,k为井壁吊挂力的标准值,MN;hd为井壁吊挂段高,取15~20 m;γh为混凝土的重力密度,MN/m3;Az为竖向钢筋截面积,m2。
环向稳定性验算按公式(11)—式(13)计算:
$$ \frac{{{L_0}}}{t} \leqslant 30 $$ (11) $$ \frac{{{E_{\text{c}}}{t^3}}}{{4r_0^3(1 - v_{\text{c}}^{\text{2}})}} \geqslant P $$ (12) $$ {L_0} = 1.814{r_0} $$ (13) 式中:L0为计算处井壁圆环计算长度,m;vc为混凝土泊松比,0.2;Ec为混凝土弹性模量,N/mm2。
在三向应力作用下的内缘承载力验算按公式(14)—式(15)计算:
$$ \begin{aligned} & \sqrt {\sigma _{\text{t}}^2 + \sigma _{\text{r}}^2 + \sigma _{\text{z}}^2 - {\sigma _{\text{t}}}{\sigma _{\text{r}}} - {\sigma _{\text{r}}}{\sigma _{\text{z}}} - {\sigma _{\text{z}}}{\sigma _{\text{t}}}} \leqslant \\& \qquad m{f_{\text{c}}} + \rho f_y' \end{aligned}$$ (14) $$ {\sigma _{\text{z}}} = \frac{{{Q_{{\text{z1}},{\text{k}}}} + {Q_{1,{\text{k}}}} + {Q_{2,{\text{k}}}} + {P_{{\mathrm{f}},{\text{k}}}}{F_{\mathrm{w}}}}}{{{A_0}}} $$ (15) 式中:σz为计算截面井壁纵向应力计算值,MN/m2;σr为计算截面井壁径向应力计算值,MN/m2;σt为井壁圆环截面切向应力,MN/m2;Fw为计算截面以上井壁外表面积,m2;A0为计算截面井壁横截面积,m2;Pf,k为计算截面以上井壁单位外表面积竖向附加力标准值,MPa。
赵固二矿西回风立井冻结深度大,穿越的黏土冲击层厚,采用C80混凝土计算所得井筒垂深680~720 m段井壁总厚为2650 mm,考虑大体积混凝土水化热、掘砌工程量等技术经济问题,采用C90—C100高强度混凝土最厚井壁总厚度为1950 mm,相较于C80混凝土减小壁厚250~700 mm,将节省投资653万元。
综合上述分析,并结合配筋设置计算,获得了西回风立井井壁结构设计参数,如表3所示。
表 3 井壁结构设计参数表深度
/ m内层井壁 外层井壁 总厚/mm 砼等级 壁厚/mm 环筋 竖筋 砼等级 壁厚/mm 环筋 竖筋 0~190 C50 450 2ϕ18@250 2ϕ18@250 C50 450 2ϕ18@250 2ϕ18@250 900 190~298 C60 600 2ϕ22@250 2ϕ22@250 C60 600 2ϕ22@250 2ϕ22@250 1200 298~420 C75 800 2ϕ28@250 2ϕ28@250 C75 700 2ϕ25@250 2ϕ25@250 1500 420~532 C80 800 2ϕ28@250 2ϕ28@250 C80 900 2ϕ28@200 2ϕ28@200 1700 532~590 C80 950 2ϕ28@200 2ϕ28@200 C80 1000 2ϕ28@200 2ϕ28@200 1950 590~640 C90 950 2ϕ28@200 2ϕ28@200 C90 1000 2ϕ28@200 2ϕ28@200 1950 640~680 C95 950 2ϕ28@200 2ϕ28@200 C95 1000 2ϕ28@200 2ϕ28@200 1950 680~720 C100 950 2ϕ28@200 2ϕ28@200 C100 1000 2ϕ28@200 2ϕ28@200 1950 720~752 C100 950 2ϕ28@200 2ϕ28@200 C80 800 2ϕ28@200 2ϕ28@200 1750 752~767 整体浇注段 1750 767~797 过渡段 1000 797~902 普通施工法段 800 4. 工程实践与创新
赵固二矿西回风立井井筒采用设计结构参数进行现场施工[18]。在400 m以深的每一段高进行监测,选取每个层位代表性的测点整理监测结果如表4所示。
表 4 部分测点位移监测序号 土层名称 测点垂深
/m壁温
/℃位移速率
/(mm·h−1)位移量
/mm1 黏土 458.2 −5.9 2.46 71.5 2 砂质黏土 534.7 −7.4 5.59 109.9 3 黏土 592.8 −5.6 3.26 85.7 4 砂质黏土 618.9 −8.4 1.90 91.7 5 铝质黏土 661.5 −9.6 3.22 88.5 6 砂质黏土 664.3 −9.8 2.43 75.1 对冻结壁进行变形监测后发现有57%的监测层位超过了50 mm,最大的变形出现在井深534.7 m处的砂质黏土层,影响因素主要是土质、深度和掘砌段高。
赵固二矿西风井工程2018年6月开始施工,2019年3月外壁掘砌完成,2019年5月套内壁结束,同月底进行壁间注浆。总体而言,掘进期间外层井壁平均每月进尺88 m,内壁平均每天进尺11.6 m,工程进展顺利,现场工作照片如图7所示。
本次设计系行业首次应用C100高强度混凝土,填补了国内高标号混凝土应用的技术空白。采用分段设计,设计分段之间的变断面处上下段井壁厚度变化较小,有利于保证井壁结构稳定;设计除考虑了采取提高混凝土早期强度措施以抵抗冻结压力外还采用敷设聚苯乙烯泡沫塑料板的方式以防止膨胀性土层迅速增长的初期冻结压力对井壁的破坏。此外,井壁设计还充分考虑了负摩擦力对井壁产生的影响,过渡段施工、支护渐变过渡,避免应力集中等。
5. 结论
(1)井筒穿越地层期间其断面原岩应力状态始终处于较大的双向偏应力场中,随着埋深的增加最大偏应力系数为1.86,但随着埋深的增加,应力绝对值增加,促使井筒断面发生非均匀破坏。
(2)基于非均匀应力场下的蝶形理论获得了塑性区最大尺寸R与井筒断面P1、P2之间的RPP特性,随着非均匀应力状态的增加,围岩塑性区在某一突变点将急剧增加,围岩强度越大,该突变点越大。
(3)对埋深500~600 m段进行数值模拟,说明了冻结施工除阻水外将更有助于围岩稳定,基于冻结法井筒设计规范,给出了赵固二矿西风井立井井筒结构参数,现场应用后监测结果显示,能够保正井筒正常掘砌的稳定性。
-
表 1 地应力分布规律
H/m σH/MPa σh/MPa σv/MPa η 100 5.42 3.08 2.45 1.75 200 7.57 4.21 4.90 1.79 300 9.72 5.34 7.35 1.81 400 11.87 6.47 9.80 1.83 500 14.02 7.60 12.25 1.84 600 16.17 8.73 14.70 1.85 700 18.32 9.86 17.15 1.86 800 20.47 10.99 19.60 1.86 表 2 冻结前后围岩力学属性
状态 岩石力学参数 体积模量/GPa 剪切模量/MPa 内摩擦角/(°) 黏聚力/MPa 抗拉强度/MPa 冻结前 0.11 0.12 20 0.81 2.57 冻结后 5.08 2.47 30 3.02 0.57 表 3 井壁结构设计参数表
深度
/ m内层井壁 外层井壁 总厚/mm 砼等级 壁厚/mm 环筋 竖筋 砼等级 壁厚/mm 环筋 竖筋 0~190 C50 450 2ϕ18@250 2ϕ18@250 C50 450 2ϕ18@250 2ϕ18@250 900 190~298 C60 600 2ϕ22@250 2ϕ22@250 C60 600 2ϕ22@250 2ϕ22@250 1200 298~420 C75 800 2ϕ28@250 2ϕ28@250 C75 700 2ϕ25@250 2ϕ25@250 1500 420~532 C80 800 2ϕ28@250 2ϕ28@250 C80 900 2ϕ28@200 2ϕ28@200 1700 532~590 C80 950 2ϕ28@200 2ϕ28@200 C80 1000 2ϕ28@200 2ϕ28@200 1950 590~640 C90 950 2ϕ28@200 2ϕ28@200 C90 1000 2ϕ28@200 2ϕ28@200 1950 640~680 C95 950 2ϕ28@200 2ϕ28@200 C95 1000 2ϕ28@200 2ϕ28@200 1950 680~720 C100 950 2ϕ28@200 2ϕ28@200 C100 1000 2ϕ28@200 2ϕ28@200 1950 720~752 C100 950 2ϕ28@200 2ϕ28@200 C80 800 2ϕ28@200 2ϕ28@200 1750 752~767 整体浇注段 1750 767~797 过渡段 1000 797~902 普通施工法段 800 表 4 部分测点位移监测
序号 土层名称 测点垂深
/m壁温
/℃位移速率
/(mm·h−1)位移量
/mm1 黏土 458.2 −5.9 2.46 71.5 2 砂质黏土 534.7 −7.4 5.59 109.9 3 黏土 592.8 −5.6 3.26 85.7 4 砂质黏土 618.9 −8.4 1.90 91.7 5 铝质黏土 661.5 −9.6 3.22 88.5 6 砂质黏土 664.3 −9.8 2.43 75.1 -
[1] 宋许根, 陈从新, 夏开宗, 等. 竖井变形破坏机制与继续使用可行性探究[J]. 岩土力学,2017,38(S1):331-342. [2] 刘孝孔, 绪瑞华, 赵艳鹏, 等. 邻近厚松散层既有立井井筒地面注浆地层加固技术[J]. 煤炭科学技术,2022,50(7):127-134. [3] 刘志强, 宋朝阳, 纪洪广, 等. 深部矿产资源开采矿井建设模式及其关键技术[J]. 煤炭学报,2021,46(3):826-845. [4] 何满潮. 深部建井力学研究进展[J]. 煤炭学报,2021,46(3):726-746. [5] 刘环宇, 陈卫忠, 王争鸣. 兖州矿区立井井筒破坏机制的理论分析[J]. 岩石力学与工程学报,2007(S1):2620-2626. [6] 赵海军, 马凤山, 徐嘉谟, 等. 金属矿充填开采岩体移动对竖井围岩的变形破坏[J]. 岩土工程学报,2012,34(2):340-348. [7] 张 好, 姚多喜, 鲁海峰, 等. 恒源煤矿井筒非采动破坏判据的弹性力学分析[J]. 煤矿安全,2017,48(2):208-211. [8] 郤保平, 李晓雪, 王 磊, 等. 基于黏弹性理论的干热岩井筒变形与损伤破坏研究[J]. 太原理工大学学报,2019,50(6):820-828. [9] 李 康, 高永涛, 周 喻, 等. 立井井筒失稳机理与加固方法[J]. 中南大学学报(自然科学版),2020,51(4):1068-1076. [10] 王胜开, 朱志根, 余一松, 等. 深部高应力竖井井壁稳定性分析及支护优化[J]. 矿业研究与开发,2022,42(4):38-44. [11] 钟应伟. 深厚第四系流砂层竖井井筒稳定性研究[D]. 湘潭: 湖南科技大学, 2013. [12] 康红普, 伊丙鼎, 高富强, 等. 中国煤矿井下地应力数据库及地应力分布规律[J]. 煤炭学报,2019,44(1):23-33. [13] 李 臣, 霍天宏, 吴 峥, 等. 动压巷道顶板非均匀剧烈变形机理及其稳定性控制[J]. 中南大学学报(自然科学版),2020,51(5):1317-1327. [14] 郭晓菲, 郭林峰, 马念杰, 等. 巷道围岩蝶形破坏理论的适用性分析[J]. 中国矿业大学学报,2020,49(4):646-653. [15] 杨更社, 奚家米, 邵学敏, 等. 冻结条件下岩石强度特性的试验[J]. 西安科技大学学报,2010,30(1):14-18. [16] GB 50511—2010 煤矿井巷工程施工规范[S]. [17] 贾成刚. 深井特厚冲积层冻结法凿井井壁结构设计[J]. 建井技术,2019,40(5):1-4. [18] 王 恒, 赵玉明, 李晓伟. 深厚黏土层冻结井筒冻结壁径向位移分析[J]. 建井技术,2019,40(6):20-23. -