吳起帆,吳靜紅,賈立翔,唐柏鑒,龔緒龍
(1.蘇州科技大學(xué) 土木工程學(xué)院,江蘇 蘇州215011;2.蘇州南智傳感科技有限公司,江蘇 蘇州215123;3.江蘇省地質(zhì)調(diào)查研究院 自然資源部地裂縫地質(zhì)災(zāi)害重點(diǎn)實(shí)驗(yàn)室,江蘇 南京210018)
地面沉降是指地面標(biāo)高持續(xù)下降的一種地質(zhì)現(xiàn)象,已成為世界各大城市的一個(gè)主要工程地質(zhì)問(wèn)題[1]。多數(shù)地面沉降是由于大量開采地下水導(dǎo)致的,地面沉降的發(fā)生發(fā)展緩慢而不易察覺,是一種“慢性”災(zāi)害,因此,對(duì)地面沉降的長(zhǎng)期監(jiān)測(cè)就顯得尤為重要。
目前,通過(guò)采用合成孔徑雷達(dá)干涉測(cè)量(InSAR)技術(shù)、全球衛(wèi)星定位系統(tǒng)(GPS)等技術(shù)進(jìn)行地面沉降監(jiān)測(cè)可以獲取地表高程變化。秦曉瓊等[2]利用永久散射體雷達(dá)干涉技術(shù)(PSInSAR技術(shù))監(jiān)測(cè)上海地面沉降,監(jiān)測(cè)結(jié)果與水準(zhǔn)數(shù)據(jù)比較驗(yàn)證,誤差在毫米級(jí),驗(yàn)證了該技術(shù)的有效性和實(shí)用性。潘登等[3]針對(duì)鄭州市地面沉降機(jī)理研究較薄弱的問(wèn)題,利用InSAR和水準(zhǔn)技術(shù)分析出該地區(qū)地面沉降主要分布地區(qū)。羅建忠等、沈玉娜等[4-5]利用GPS監(jiān)測(cè)技術(shù)對(duì)地面沉降區(qū)域地表形變進(jìn)行數(shù)據(jù)處理與分析,驗(yàn)證了該技術(shù)應(yīng)用于地面沉降監(jiān)測(cè)的可行性。然而這些監(jiān)測(cè)技術(shù)主要適用于宏觀的區(qū)域性地面沉降趨勢(shì)分析,不能對(duì)地下不同深度土層的變形量進(jìn)行直接監(jiān)測(cè)。對(duì)于深部土體位移監(jiān)測(cè),現(xiàn)階段主要采用基巖標(biāo)和分層沉降標(biāo),如王美華等[6]以“三維立體化監(jiān)測(cè)”為目標(biāo),提出了土體深層(分層)沉降監(jiān)測(cè)技術(shù)的發(fā)展趨勢(shì),并建立包含空中、地面、土體內(nèi)部的同步監(jiān)測(cè),實(shí)現(xiàn)變形監(jiān)測(cè)智能評(píng)價(jià)和安全預(yù)警。王剛等[7]對(duì)鄭州市地面沉降監(jiān)測(cè)進(jìn)行場(chǎng)地選址、建標(biāo)深度確定等一系列施工,采用三牙輪鉆頭寶塔式鉆具組合不取心鉆進(jìn)施工工藝,對(duì)基巖標(biāo)的建設(shè)有一定的借鑒意義和推廣價(jià)值。但是該方法成本較高且只能對(duì)特定深度土層進(jìn)行變形監(jiān)測(cè),無(wú)法實(shí)現(xiàn)鉆孔全斷面土層變形分析。
分布式光纖感測(cè)技術(shù)作為一種新型的監(jiān)測(cè)技術(shù),能夠獲得地面沉降過(guò)程中多場(chǎng)多參量數(shù)據(jù),實(shí)現(xiàn)精細(xì)化的沉降過(guò)程監(jiān)測(cè),近些年已被試用于地面沉降監(jiān)測(cè)中[8]。顧凱等[9-11]利用它對(duì)一個(gè)300 m深鉆孔進(jìn)行了監(jiān)測(cè),并對(duì)土層結(jié)構(gòu)變形進(jìn)行分析,通過(guò)識(shí)別沉降和變形演化過(guò)程,提供了沿整個(gè)河道的土體變形分布清晰圖像,證明了含水層的變形受孔隙水壓力的顯著控制,認(rèn)為該技術(shù)在地面沉降監(jiān)測(cè)中具有穩(wěn)定的長(zhǎng)期監(jiān)測(cè)性能。
本文結(jié)合全分布式光纖監(jiān)測(cè)技術(shù)布里淵光時(shí)域反射技術(shù)(Brillouin Optical Time Domain Reflectometry,BOTDR),和準(zhǔn)分布式光纖監(jiān)測(cè)技術(shù)布拉格光纖光柵(Fiber Bragg Grating,F(xiàn)BG),對(duì)鹽城市響水縣陳家港中心小學(xué)內(nèi)的一個(gè)350 m深鉆孔土層進(jìn)行了近兩年半的監(jiān)測(cè),分析評(píng)價(jià)了該地區(qū)地下含水層水體變動(dòng)情況及地下土體受水體變動(dòng)影響下的變形特征。
BOTDR技術(shù)是基于入射光波與聲波或傳播的壓力波在光纖中傳播時(shí)產(chǎn)生的非彈性碰撞而出現(xiàn)的光散射[12],可以看成入射光在移動(dòng)光柵上的散射。其應(yīng)變測(cè)量原理如圖1所示。當(dāng)散射光能量達(dá)到最高時(shí),能夠確定位置Z處的布里淵頻率漂移VB,布里淵散射光的頻移VB與光纖應(yīng)變?chǔ)懦示€性關(guān)系,即
其中,VB(ε)是當(dāng)應(yīng)變?yōu)棣艜r(shí)的布里淵頻率的漂移量;VB0為在測(cè)試環(huán)境溫度不變的條件下,光纖自由狀態(tài)時(shí)的布里淵頻率漂移量;為光纖的應(yīng)變系數(shù);ε為光纖的實(shí)際發(fā)生應(yīng)變量。
圖1 BOTDR的應(yīng)變測(cè)量原理圖[12]
FBG技術(shù)是通過(guò)改變光纖芯區(qū)折射率,使其產(chǎn)生小的周期性調(diào)制而形成的,其工作原理如圖2。當(dāng)溫度或應(yīng)變改變后,光纖的折射率發(fā)生了改變,從而引起光柵波長(zhǎng)偏移。利用應(yīng)變和溫度之比與光柵波長(zhǎng)偏移量的線性關(guān)系,得出被測(cè)體應(yīng)變量[13]。波長(zhǎng)漂移量與應(yīng)變和溫度的關(guān)系可表示為?
其中,Δλβ為FBG中心波長(zhǎng)的變化量;λβ為FBG中心波長(zhǎng);Pe為有效光彈系數(shù);ε為光纖軸向應(yīng)變;ΔT為溫度變化量;α為光纖的熱膨脹系數(shù);ζ為光纖的熱光系數(shù)。
準(zhǔn)分布式感測(cè)技術(shù)利用一根信號(hào)傳導(dǎo)光纖,將多個(gè)FBG傳感器串聯(lián)起來(lái),通過(guò)波分復(fù)用和時(shí)分復(fù)用等感測(cè)原理,將多個(gè)傳感器的感測(cè)信號(hào)區(qū)分而獲得各個(gè)傳感器的感測(cè)信息。
圖2 FBG傳感系統(tǒng)原理[13]
鹽城響水地貌為平原,大面積覆蓋第四紀(jì)與新近紀(jì)沉積物,多為淤泥、淤泥質(zhì)亞黏土、亞黏土、粉質(zhì)黏土、粉砂、亞砂土等多層狀結(jié)構(gòu)的松散沉積物。近年來(lái),鹽城響水通過(guò)圍海造陸、工程建設(shè)來(lái)發(fā)展經(jīng)濟(jì),對(duì)區(qū)域內(nèi)大量高壓縮性、高含水量和低承載力的軟土產(chǎn)生的較大擾動(dòng),土體穩(wěn)定性差,易產(chǎn)生地面沉降[14]。截止到2016年全市累計(jì)地面沉降量均大于100 mm,而累計(jì)地面沉降量大于200 mm的區(qū)域面積為10 812 km2,約占鹽城市總面積的63.86%(注:江蘇省地質(zhì)勘察技術(shù)院調(diào)查結(jié)果)。顧晟彥等[15]利用MapGIS中的空間分析功能將分布在鹽城北部的濱海、阜寧、響水射陽(yáng)縣城等定為高危險(xiǎn)區(qū),如圖3所示。據(jù)此,將地面沉降觀測(cè)鉆孔地點(diǎn)選在鹽城響水,鉆孔具體位置處于鹽城市響水縣陳家港中心小學(xué)內(nèi)(見圖3)。
圖3 地面沉降危險(xiǎn)性分區(qū)圖[15]及光纖監(jiān)測(cè)孔位置示意圖
監(jiān)測(cè)鉆孔深度為350 m,傳感器具體布設(shè)安裝方案如圖4所示。為實(shí)現(xiàn)鉆孔全斷面地層巖土體變形監(jiān)測(cè),選用了地層定點(diǎn)雙芯全分布式傳感光纜,其內(nèi)部有兩根光纖,一根為250μm直徑光纖,另一根為900μm直徑光纖,兩根光纖均為局部定點(diǎn)松套光纜,可實(shí)現(xiàn)地層變形分段均一化測(cè)量和較大變形測(cè)量。
根據(jù)相關(guān)地層資料,在0~350 m深度范圍內(nèi),分別在72 m、166 m和197 m含水層布設(shè)1支光纖光柵滲壓計(jì)用于監(jiān)測(cè)含水層內(nèi)滲水壓力變化情況,重點(diǎn)進(jìn)行水位監(jiān)控。此外,對(duì)于抽水含水層內(nèi)及其上部相鄰隔水層分別布設(shè)1支15 m標(biāo)距的光纖光柵位移計(jì),用于監(jiān)測(cè)其變形情況。其中,抽水含水層內(nèi)位移計(jì)頭部深度101 m,尾部深度86 m,相鄰隔水層內(nèi)位移計(jì)頭部深度85 m,尾部深度70 m。
在鉆孔成孔完成后,對(duì)鉆孔進(jìn)行一次掃孔,并用清水進(jìn)行洗孔處理。將光纜和光纖光柵傳感器引線穿入到配重導(dǎo)頭內(nèi),根據(jù)設(shè)計(jì)方案串接連線,將各光纜和引線串接成回路。使用環(huán)氧樹脂膠將光纜與導(dǎo)頭粘結(jié),在其外面使用布基膠帶纏繞固定[16],如圖5所示。將兩條定點(diǎn)光纜上的兩條受力麻繩,交叉綁扎固定在導(dǎo)頭上(見圖5(a))。在鉆孔周圍場(chǎng)地上,將位移計(jì)引標(biāo)安裝固定好(見圖5(b))。將固定好光纜、引線和麻繩等配重導(dǎo)頭放入鉆孔內(nèi)部。將鉆桿放入配重導(dǎo)頭尾部套管內(nèi),通過(guò)下放鉆桿將光纜帶入到鉆孔深部(圖5(c)),待光纜下放到底部后,使用儀器對(duì)光纜進(jìn)行檢測(cè)。最后采用封孔材料進(jìn)行回填封孔。
圖4 光纖傳感器布設(shè)方案
圖5 監(jiān)測(cè)實(shí)施(a:導(dǎo)頭連接;b:位移計(jì)安裝;c:光纜下放)
陳家港監(jiān)測(cè)孔光纖監(jiān)測(cè)點(diǎn)于2014年9月28日布設(shè)完成,封孔材料與周圍土體固結(jié)耦合周期為2~3個(gè)月,初始值采集在封孔材料回填接近3個(gè)月之后的2014年12月14日采集。將采集的初始值數(shù)據(jù)作為數(shù)據(jù)分析的基準(zhǔn)數(shù)值(分布曲線應(yīng)變基準(zhǔn)值),之后監(jiān)測(cè)周期約為3個(gè)月一次,并進(jìn)行長(zhǎng)期監(jiān)測(cè)。
將2014年12月14日的監(jiān)測(cè)數(shù)據(jù)作為初始數(shù)值,之后監(jiān)測(cè)周期所測(cè)得的應(yīng)變數(shù)值減去初始數(shù)值,即可得到各監(jiān)測(cè)周期內(nèi)感測(cè)光纜的應(yīng)變變化值,根據(jù)BOTDR監(jiān)測(cè)所得數(shù)據(jù),截至2017年6月25日,相應(yīng)監(jiān)測(cè)周期內(nèi),雙芯光纜的兩種光纖的應(yīng)變分布曲線見圖6。陳家港監(jiān)測(cè)孔地層定點(diǎn)250μm和900μm光纖監(jiān)測(cè)趨勢(shì)基本相同,在15 m深度范圍以上正應(yīng)變急劇減小為負(fù)應(yīng)變,這是由于淺表土體圍壓較小造成的預(yù)拉狀態(tài)下光纜的回縮,并且該深度地表受到大氣溫度變化而影響。在15~25 m深度范圍壓應(yīng)變隨監(jiān)測(cè)周期逐漸增大,表面該深度范圍內(nèi)土層處于壓縮狀態(tài)。其余層段應(yīng)變基本上處于零值附近,土體未發(fā)生明顯的變形。
將感測(cè)光纜的應(yīng)變值沿深度方向進(jìn)行積分,可以得到土層近似變形量(見圖7)。在深度上,土層變形主要發(fā)生在200 m以上,200 m以下地層變形量較小。其中,0~100 m范圍變形波動(dòng)較大,這主要與地下水的波動(dòng)密切相關(guān)。在時(shí)間上,隨著監(jiān)測(cè)次數(shù)的增多,全斷面土層變形呈逐漸增大的趨勢(shì),說(shuō)明現(xiàn)階段該地區(qū)土層仍有持續(xù)壓縮。圖8給出了陳家港鉆孔布設(shè)的含水層內(nèi)部(頭部深度103 m和尾部深度88 m)及含水層上覆隔水層(頭部深度87 m和尾部深度72 m處)內(nèi)的2支光纖光柵位移計(jì)監(jiān)測(cè)到的土層變形以及相對(duì)應(yīng)位置的全分布式光纖監(jiān)測(cè)到的72~87 m、88~103 m土層變形圖。位移計(jì)監(jiān)測(cè)到的土層變形顯示72~87 m弱透水層呈現(xiàn)出壓縮回彈再壓縮的波動(dòng)變形,整體上有回彈趨勢(shì),壓縮主要發(fā)生在6月至8月的夏季。88~103 m含水層變形同樣存在波動(dòng)變形,但整體上呈壓縮狀態(tài),最大壓縮量在7 mm左右。而全分布式光纖監(jiān)測(cè)到的72~87 m和88~103 m土層變形趨勢(shì)相同,整體呈現(xiàn)出壓縮回彈再壓縮,壓縮主要也發(fā)生在夏季,最大壓縮量在6 mm左右??梢园l(fā)現(xiàn)全分布式監(jiān)測(cè)與準(zhǔn)分布式監(jiān)測(cè)存在一定差異,這是因?yàn)槿植际奖O(jiān)測(cè)得到的應(yīng)變受BOTDR空間分辨率(1 m)的影響,同時(shí)受上下土層應(yīng)變的相互影響。
圖6 陳家港監(jiān)測(cè)孔應(yīng)變分布曲線圖
圖7 陳家港地層變形隨深度的變化關(guān)系圖
圖8 光纖光柵位移計(jì)與全分布式監(jiān)測(cè)土層變形趨勢(shì)圖
在72 m、166 m及197 m處布設(shè)的三支光纖光柵滲壓計(jì)監(jiān)測(cè)了不同含水層的孔隙水壓力變化情況,圖9給出了這三個(gè)含水層孔隙水壓力的變化趨勢(shì)。72 m含水層孔隙水壓力最小,197 m含水層孔隙水壓力最大,即含水層越深孔隙水壓力越大。從圖中可以看出孔隙水壓力總體是先下降后上升的趨勢(shì),三個(gè)含水層孔隙水壓力在2016年6月15日均為最低值,該時(shí)間為夏季,地下水消耗量及蒸發(fā)量大,導(dǎo)致地下水位下降??偟膩?lái)說(shuō),三個(gè)含水層孔隙水壓力呈現(xiàn)波動(dòng)變化,夏季降低,冬季回升。72 m含水層孔隙水壓力在初期是呈下降趨勢(shì),說(shuō)明該區(qū)域含水層監(jiān)測(cè)開始時(shí)在經(jīng)歷排水階段。而166 m和197 m含水層孔隙水壓力初期有一個(gè)上升階段然后再下降,說(shuō)明含水層初期存在補(bǔ)水階段,隨著排水的進(jìn)行,孔隙水壓力慢慢下降直到最低點(diǎn)處。
圖10給出了三個(gè)含水層不同時(shí)期累計(jì)變形量與其孔隙水壓力變化的關(guān)系。從圖中可以看出,72 m深度滲壓計(jì)的壓力整體上先減小后增大再減小,而在72 m含水層則先出現(xiàn)了回彈壓縮再回彈,總體上處于回彈變形,這是由于埋深較小的土層變形除了受水位變化影響外,還受上部建筑荷載等其他因素影響;166 m和197 m深度處滲壓計(jì)的壓力先是緩慢增加后突然減小,后又緩慢增加,呈V字型。166 m和197 m的含水層變形隨時(shí)間的變化也基本一致,呈現(xiàn)出壓縮回彈再壓縮的趨勢(shì),且壓縮量明顯大于回彈量,并且197 m處的含水層壓縮變形量要高于166 m出的含水層壓縮變形量,可以發(fā)現(xiàn)深部土層變形主要受水位變化影響,變形略滯后于水位。
不同的含水層其水位變化不同,出現(xiàn)的土體變形趨勢(shì)也并不完全一致?,F(xiàn)階段開采條件若不變,深部含水層將繼續(xù)發(fā)生壓縮變形。
表1給出了三個(gè)含水層各階段變形及沉降量,從表中可以看出各含水層年沉降量略有差異,其中72 m含水層(Ⅰ)的年沉降量約為1.74 mm,而166 m含水層(Ⅱ)和197 m含水層(Ⅲ)年沉降量分別為5.66 mm和10.4 mm,呈現(xiàn)出隨著含水層的加深而沉降量增大的趨勢(shì)。并且,從每一含水層的月沉降量也可以看出,在每一年的夏季地下水位下降時(shí),沉降量也明顯增大。
圖9 孔隙水壓力隨時(shí)間的變化關(guān)系圖
圖10 孔隙水壓力與各含水層變形關(guān)系圖
表1 三個(gè)含水層各階段變形和沉降量
本文將全分布式監(jiān)測(cè)技術(shù)BOTDR與準(zhǔn)分布式監(jiān)測(cè)技術(shù)FBG相結(jié)合,并應(yīng)用于鹽城響水地面沉降監(jiān)測(cè)中,對(duì)陳家港中心小學(xué)一個(gè)350 m深鉆孔進(jìn)行了兩年半的全斷面精細(xì)化監(jiān)測(cè),得到如下結(jié)論:
(1)陳家港地層應(yīng)變變化較為明顯的是地表至地下25 m范圍內(nèi),0~15 m范圍內(nèi)存在拉應(yīng)變和壓應(yīng)變交替變化的情況,主要是由于大氣溫度的變化較為頻繁及圍壓較小情況下的預(yù)拉回彈;15~25 m土層壓應(yīng)變逐漸增大,為壓縮主要層位,其余深度土層變形相對(duì)穩(wěn)定。
(2)研究區(qū)地下水波動(dòng)具有循環(huán)過(guò)程,呈現(xiàn)出夏季降低,冬季回升的特點(diǎn)。
(3)淺部含水層變形受多種因素影響,與水位變化沒(méi)有直接關(guān)系;深部含水層整體上均呈現(xiàn)出先壓縮再回彈再壓縮趨勢(shì),與水位變化相關(guān),且變形略滯后于水位變化。
(4)三個(gè)監(jiān)測(cè)含水層目前仍有沉降壓縮趨勢(shì),且年沉降量呈現(xiàn)出隨著含水層的加深而增大的趨勢(shì),由于夏季地下水為下降明顯,沉降也明顯增大,其沉降程度較大。