
Citation: | PENG Hua, MA Xiu-min, WANG Zhen, et al., 2013. A FIBRE OPTIC EXTRINSIC FABRY-PEROT INTERFEROMETER WITH TEMPERATURE COMPENSATION FOR FAULT MEASUREMENT. Journal of Geomechanics, 19 (3): 315-324. |
盾构隧道在砂土地层中施工时开挖面失稳可能性很大,一旦发生会造成极其严重的损失。防止开挖面失稳的关键是合理设置不同盾构支护平衡模式下的支护压应力。针对这一问题,目前主要的研究成果集中在临界状态下开挖面的整体失稳上,如工程中广泛应用的筒仓楔形体模型的极限平衡法[1]和塑性上、下限理论极限分析法[2~3]均是假设开挖面整体失稳。然而,工程实践和模型试验都证明盾构隧道开挖面失稳主要以局部失稳为主。Chambon[4]和Oblozinsky[5]通过三维离心试验发现,与刚性活塞支护模式下全断面失稳不同,静水压应力支护模式下开挖面上半部失稳。Berthoz[6]利用大比例土压盾构模型试验研究发现,土压平衡支护模式下开挖面上半部失稳。这说明,不同支护模式下,开挖面失稳区不同。对于发生这种现象的原因,目前还未见文献报道。随着国内外盾构隧道建设的日益增多,合理确定不同支护模式下盾构支护压应力尤为重要,要确定盾构支护压应力首先必须明确开挖面前方土体对刀盘的压应力分布形式。
不同支护模式下,开挖面前方土体对刀盘的压应力分布模式和大小直接决定了支护压应力的大小和作用效果。盾构掘进遇到障碍物如切削桩基时[7]和双洞盾构[8]近距离施工时,合理的支护压应力对隧道的稳定性起着决定性作用。但长期以来相关学者一直假定楔形体为刚性且开挖面前方土体对刀盘的压应力均匀分布[9~10],不考虑刀盘与松动土体之间的摩擦力,而且水平土压应力系数依赖经验确定。武军等[11]在考虑松动土体内的三维土拱效应和开挖土体与刀盘的摩擦力的基础上,改进了筒仓楔形体模型,并通过与模型试验结果的对比分析,验证了该计算方法的合理性。在此基础上,文章进一步通过理论推导,分析开挖面所受应力的分布形式,研究了不同支护方式下的开挖面稳定性,最后给出了公式的简化计算方法,以方便实际工程运用。
根据武军等[11]的推导,开挖面前方松动土体对刀盘的水平压应力σh为
σh=(fz2+gz+σH)k2 | (1) |
公式(1)中:z为以计算点与盾构最低点的垂向距离,m;σH为楔形体上覆土压应力,kPa;可由公式(2)求得:
σH=Rγ−ck1tanφ(1−e−k1tanφHR)+σ0e−k1tanφHR | (2) |
公式(2)中:c为土的粘聚力,kPa;γ为土体容重,kN ·m-3;σ0为地表竖向均布荷载,kPa;φ为土体内摩擦角,(°);H为隧道上覆土厚度,m;R为参数,可按公式(3)确定:
R=BDtanω2(B+Dtanω) | (3) |
公式(3)中:B=(πD)/4;ω= 45°-φ/2;D为盾构隧道开挖直径,m。当σ0=0且c=Rγ时,σH为负值,此时取σH=0;f, g为计算参数,其表达式分别如下:
f=σH2D2[(a+bD)2−a]−r2D(a+bD) | (4) |
g=r−σH(b+aD) | (5) |
r、a、b为计算参数:
r=γ−2cB | (6) |
a=k2tanδtanω | (7) |
b=2k1tanφB | (8) |
公式中δ为松动土体与刀盘间的摩擦角,(°);k1, 2[12]为土拱区内的水平土压力系数,其表达式为:
k1,2=3(cos2θ1,2+kasin2θ1,2)3−(1−ka)cos2θ1,2 | (9) |
公式(9)中:ka为郎肯主动土压力系数,ka=tan2(45°-φ/2);θ1, 2为公式(10)、(11)确定的较大值:
θ1=arctan(1−ka±√(1−ka)2−4tan2δka2tanφka) | (10) |
θ2=arctan(1−ka±√(1−ka)2−4tan2δka2tanδka) | (11) |
将公式(1)沿隧道直径D积分,求得松动土体对刀盘的平均压应力σah:
σah=(13fD2+12gD+σC)k2 | (12) |
目前常用的盾构支护模式有气压平衡支护、泥水平衡支护和土压平衡支护三种。理论上有效支护压应力均匀分布的只有开挖面在地下水位以上时的气压平衡支护模式,如图 1a所示;当开挖面在地下水位以下时,一部分气压支护力用以平衡地下水压应力,由于水压应力沿深度呈线性分布,因此作用在开挖面上的有效支护压应力呈上大下小的梯形分布,如图 1b所示;在泥水平衡支护模式下,当开挖面在地下水位以上时,支护力沿竖直方向呈上小下大的梯形分布,梯度为γs(γs为泥水容重),如图 1c所示;当开挖面在地下水位以下时,一部分泥水支护力用以平衡地下水压应力,作用在开挖面上的有效支护压应力也呈上小下大的梯形分布,但梯度减小为γs-γw(γw为地下水容重),根据膨润土和含泥量的不同,泥水容重也不同,但一般情况下在11~12 kN ·m-3之间,因此有效泥水支护应力梯度很小,近似于均匀分布,如图 1d所示;在土压平衡支护模式下,由于渣土会被改良为具有流动性的塑形状态,可忽略渣土的剪切强度,因此渣土压应力分布可假设为静水压应力分布,当开挖面在地下水位以上时,支护力沿竖直方向线性分布,梯度为γm(γm为渣土容重),当开挖面在地下水位以下时,有效支护力沿竖直方向线性分布,梯度为γm-γw,一般情况下,土压平衡模式支护力梯度大于泥水平衡模式支护力梯度,如图 1e、1f所示。
文章假设不考虑渗透力对开挖面稳定性的影响。为不失一般性,取γ=18 kN ·m-3、γ′=10 kN ·m-3、γw=10 kN ·m-3、γs=12 kN ·m-3、γm=18 kN ·m-3。如图 2所示,在三种平衡支护模式下,开挖面失稳均非全断面整体失稳,而是局部失稳。
气压平衡支护模式下,开挖面在地下水位以上,当φ=0°时,不论粘聚力多大,土体均从刀盘下半部分挤入,开挖面下半部失稳; 当φ>0°时,随着内摩擦角的增大,失稳区逐步从开挖面下半部分向上移动,同时当φ值一定时,随着粘聚力的减小,失稳区逐步上移,如图 2a所示。气压平衡支护模式下,开挖面在地下水位以下时,由于有效气压支护力上大下小分布,开挖面均在下半部分失稳,如图 2b所示。需要说明的是:模型试验中采用的气压支护与实际工程中的并不相同,模型试验中的气体并不直接作用于土体而是被填充在橡胶袋中,通过橡胶薄膜支护开挖面,开挖面局部失稳后,失稳土体挤压橡胶膜,由于橡胶膜的张力使得作用在剩余稳定土体的压应力减小,继而引起失稳部位上部未稳定土体坍塌,所以最终表现为Oblozinsky[5]、Thorpe[14]试验中的整体失稳或上部约3/4断面失稳。实际工程中,开挖面局部土体失稳后形成类似于重力挡土墙的堆积体阻止开挖面的进一步失稳,堆积体的增大,最终能阻止开挖面的失稳。
泥水平衡支护模式下,开挖面在地下水位以上,当φ=0°时,土体从刀盘下半部分挤入,开挖面下半部失稳,但是当φ≠0°时,开挖面均在上半部失稳,如图 2c所示,这与Chambon[4]和Oblozinsky[5]的离心模型试验结果规律性是一致的。开挖面在地下水位以下,当φ=0°时,土体从刀盘下半部分挤入,开挖面下半部失稳,而当φ>0°时,随着内摩擦角的增大,失稳区逐步从开挖面下半部分向上移动,同时当φ值一定时,失稳区随着粘聚力的减小,失稳区逐步上移,如图 2d所示。这与气压平衡支护模式下,开挖面在地下水位以上时,开挖面的失稳规律相同,其原因是该情况下有效泥水支护应力梯度很小,近似于均匀分布。
土压平衡支护模式下,开挖面在地下水位以上时,开挖面均在上半部失稳,如图 2e所示。这与Berthoz[6]的大比例模型试验结果规律性是一致的,造成这种现象的原因是该情况下支护力梯度较大导致开挖面下半部支护压应力大于松动土压应力。土压平衡支护模式下,开挖面在地下水位以下时,由于有效支护力梯度的减小,当φ=0°时,土体从刀盘下半部分挤入,开挖面下半部失稳,而当φ>0°时,随着内摩擦角的增大,失稳区逐步从开挖面下半部分向上移动,同时当φ值一定时,失稳区随着粘聚力的减小,失稳区逐步上移,如图 2f所示,这与图 2a、2d表示的规律相同。
根据武军等[11]的分析,粘聚力不影响开挖土体对刀盘的压应力的分布形式,当C/D≥1.5时,埋深对于松动土体对刀盘的压应力大小和分布形式影响不明显;当φ≥30°时,砂土内摩擦角对于松动土体对刀盘的压应力大小和分布形式影响也不明显;而刀盘与开挖土体间的外摩擦角大小直接影响开挖土体对刀盘的压应力大小和分布形式,因此有必要分析不同支护模式下开挖面失稳和刀盘与开挖土体间外摩擦角的关系。
气压平衡支护模式下,开挖面在地下水位以上时,随着外摩擦角发挥程度的增加,开挖面失稳区逐步从开挖面下半部分向上移动,外摩擦角的发挥程度与地层条件、刀盘开口率有直接关系,如图 3a所示。在透气性小的土体中,气体有可能会在刀盘与土体形成厚度很小的气垫层,减小外摩擦角。刀盘开口率较大时,气体和松动土体直接接触的面积增加,而气体和土体间的摩擦角近似为零,因此开口率增大也会降低外摩擦角。气压平衡支护模式下,开挖面在地下水位以下时,开挖面失稳区位于下半部,与外摩擦的发挥程度关系不大,如图 3b所示。泥水平衡支护模式下,外摩擦角的大小主要与泥膜的质量和厚度有关,泥膜质量好、厚度大,外摩擦角就小。当开挖面在地下水位以上时,外摩擦角的大小对开挖面失稳区域的影响不大,失稳区均位于开挖面上半部分;而开挖面在地下水位以下时,失稳区也位于开挖面上半部分,但失稳程度不同,如图 3c、3d所示。这说明开挖面在地下水位以上时,泥膜的的质量和厚度对开挖面最终的失稳破坏影响不大,但在有地下水的条件下,泥膜的质量和厚度对开挖面最终的失稳破坏程度有较大影响。土压平衡支护模式下,外摩擦角的大小主要与渣土改良的好坏和刀盘开口率的大小有关,渣土改良质量好、刀盘开口率的大,外摩擦角就小。无论开挖面在地下水位以上还是以下,开挖面均在其上部失稳,外摩擦角的大小对其影响不大,如图 3e、3f所示。这说明渣土改良和刀盘开口率的大小对开挖面最终的失稳影响不大。
支护压应力设定后,在盾构停机状态下,支护压应力直接通过刀盘上的开口作用于土体,在盾构掘进状态下,支护压应力对土体的作用还受刀盘转动摩擦力和土体进入土舱过程中刀盘对土体产生的向上滑动摩擦力的影响,但是当刀盘旋转切削土体时,土体对刀盘的转动摩擦力合力为一个力偶,该力偶与刀盘扭矩相互平衡,因此土体与刀盘间的转动摩擦力不对开挖面稳定产生影响[11],刀盘对土体产生的向上滑动摩擦力由外摩擦决定。但通过以上分析发现,对于土拱效应明显的砂土地层来说,除气压平衡支护模式下其开挖面在地下水位以上时,外摩擦角对开挖面失稳有较大影响,其余支护模式下,外摩擦角对开挖面失稳影响较小。因此,在砂性地层中,支护压应力的大小设定对防止开挖面失稳最为关键。需要说明的是:由该节分析可知,在各种不同支护模式下,并未发生如筒仓楔形体模型假设的开挖面整体失稳,而是发生了开挖面的局部失稳,而某些局部失稳土体所受的拉应力很大,开挖面的变形随着支护压应力的变小逐步增大,其所受的应力也是逐步变化的,这说明在未到达筒仓楔形体模型所假设的开挖面整体失稳前,开挖面已经发生了局部失稳。这种机理还有待研究,但是对于盾构隧道来说,以筒仓楔形体模型确定的极限稳定支护力是不安全的。
公式(1)和(11)计算过程较为复杂,因此给出C/D=2.0、c=0 kPa时,工程实践中常遇到的土体参数范围内的设计参数k2、f、g、σC无量纲化关系图,如图 4所示,以方便工程设计使用。根据文献[11],当φ≥15°且C/D≥1.5时,粘聚力c和埋深C对松动土体对刀盘压应力的影响不大,因此,图 4也可用于C/D≥1.5、c>0 kPa的情况,且偏于安全。
进一步分析可知,水平土压应力系数k2随着内摩擦角φ的增大而减小,同时当δ/φ≤0.85时,k2随着δ/φ的增加缓慢增大,当δ/φ≥0.85时,k2随着δ/φ的增加快速增大,如图 4a所示。(fD)/γ随着δ/φ的增加而减小,同时当φ≤40°时,(fD)/γ随着φ的增加而快速减小,当φ≥40°时,减小速率变的较为缓慢,如图 4b所示。当φ≤20°时,g/γ随着φ的增大缓慢减小,当φ≥20°时,g/γ随着φ的增大快速增大,如图 4c所示。当φ≤40°时,σC/(γD)随着φ的增大快速减小,当φ≥40°时,减小速率变的较为缓慢,如图 4d所示。
文章通过分析在气压支护模式、泥水支护模式和土压支护三种模式下,盾构隧道开挖面分别在地下水位以上和地下水位以下时所受应力的分布形式以及稳定性,得出如下结论:
(1) 刀盘与松动土体之间的摩擦力对开挖面前方土体对刀盘的压力大小和分布形式有重要的影响,对开挖面稳定性和失稳部位也有重要影响。
(2) 目前工程中广泛使用的建立在开挖面前方土体对刀盘的压力均匀分布和整体失稳基础上的筒仓楔形体模型无法解释模型试验中发现的静水压应力支护模式下开挖面上半部失稳和土压平衡支护模式下开挖面半部失稳的现象,而改进的筒仓楔形体模型能合理解释这一现象,用该改进模型分析开挖面稳定更合理、更安全。
(3) 在气压、泥水和土压平衡支护模式下,开挖面失稳均非全断面整体失稳,而是局部失稳。有效支护应力均匀分布时,除粘土开挖面下部失稳外,其余土体均为开挖面中下部失稳。有效支护压应力呈上小下大的梯形分布时,除软粘土开挖面下部失稳外,其余土体均为开挖面上部失稳。有效支护应力呈上大下小的梯形分布时,所有土体开挖面均为下部失稳。
(4) 对于土拱效应明显的砂土地层来说,除开挖面在地下水位以上其为气压支护模式时,刀盘和土体间的摩擦角发挥程度对开挖面的失稳区域分布有较大影响外,其余支护模式下,刀盘和土体间的摩擦角发挥程度对开挖面的失稳区域分布影响不大。
[1] |
Bhatia V, Murphy K A, Claus R O, et al. Recent developments in optical-fiber-based extrinsic Fabry-Perot interferometric strain sensing technology[J]. Smart Materials and Structures, 1995, 4(4):246~251. doi: 10.1088/0964-1726/4/4/004
|
[2] |
Lee B. Review of the present status of optical fiber sensors[J]. Optical Fiber Technology, 2003, 9(2):57~79. doi: 10.1016/S1068-5200(02)00527-8
|
[3] |
Murphy K A, Gunther M F, May R G, et al. EFPI sensor manufacturing and applications[C]//Crowe C R. Smart Structures and Materials 1996:Industrial and commercial applications of smart structures technologies. San Diego:The International Society for Optical Engineering, 1996:476~482.
|
[4] |
Xiao H, Deng J D, Wang Z Y, et al. Fiber optic pressure sensor with self-compensation capability for harsh environment applications[J]. Optical Engineering, 2005, 44(5):544~547. https://www.netl.doe.gov/File%20Library/Events/2014/crosscutting/Crosscutting_20140521_1600A_Clemson.pdf
|
[5] |
Ferreira L A, Ribeiro A B L, Santos J L, et al. Simultaneous measurement of displacement and temperature using a low finesse cavity and a fiber Bragg grating[J]. IEEE Photonics Technology Letters, 1996, 8(11):1519~1521. doi: 10.1109/68.541569
|
[6] |
Ferreira L A, Ribeiro A B L, Santos J L, et al. Simultaneous displacement and temperature sensing using a white light interrogated low finesse cavity in line with a fiber Bragg grating[J]. Smart Materials and Structures, 1998, 7(2):189~194. doi: 10.1088/0964-1726/7/2/006
|
[7] |
Liu T, Fernando G F, Rao Y J, et al. Simultaneous strain and temperature measurements in composites using a multiplexed fibre Bragg grating sensor and an extrinsic Fabry-Perot sensor[C]//Claus R O. Smart structures and materials 1997:Smart Sensing, Processing, and Instrumentation. San Diego:The International Society for Optical Engineering, 1997:203~212.
|
[8] |
Bremer K, Lewis E, Moss B, et al. Conception and preliminary evaluation of an optical fibre sensor for simultaneous measurement of pressure and temperature[J]. Journal of Physics:Conference Series, 2009, 178(1):12~16. http://www.dbpia.co.kr/Journal/ArticleDetail/NODE00645252
|
[9] |
Bremer K, Lewis E, Leen G, et al. 2011 Temperature compensated miniature all-glass fibre optic pressure sensor[C]//Institute of Electrical and Electronics Engineers. Sensors. Piscataway:IEEE, 2011:105~108.
|
[10] |
Liu T, Fernando G F, Zhang Z Y, et al. Simultaneous strain and temperature measurements in composites using extrinsic Fabry-Perot interferometric and intrinsic rare-earth doped fiber sensors[J]. Sensors and Actuators A:Physical, 2000, 80(3):208~215. doi: 10.1016/S0924-4247(99)00309-X
|
[11] |
Zhang Z Y, Grattan K T V, Palmer A W, et al. Potential for temperature sensor applications of highly neodymium-doped crystals and fiber at up to approximately 1000℃[J]. Review of Scientific Instruments, 1997, 68(7):2759~2763. doi: 10.1063/1.1148191
|
[12] |
Xu J C, Wang X W, Cooper K L, et al. Miniature temperature-insensitive Fabry-Perot fiber-optic pressure sensor[J]. Photonics Technology Letters, 2006, 18(10):1134~1136. doi: 10.1109/LPT.2006.873361
|
[13] |
Jiang Y. 2008 Fourier transform white-light interferometry for the measurement of fiber-optic extrinsic Fabry-Pérot Interferometric sensors[J]. Photonics Technology Letters, 20(2):75~77. doi: 10.1109/LPT.2007.912567
|
[14] |
Jiang Y, Tang C J. High-finesse micro-lens fiber-optic extrinsic Fabry-Perot interferometric sensors[J]. Smart Materials and Structures, 2008, 17(5):55~67. doi: 10.1007/s13320-010-0014-z
|
1. | 彭丽娟,白伟利,李文超,高龙飞. 2022年9月泸定M_S 6.8地震与花莲M_S 6.9地震震害对比分析. 地震地磁观测与研究. 2024(03): 139-144 . ![]() |