李金文，楊清云，韓桂武，楊文艷

1.國(guó)家管網(wǎng)集團(tuán)山東省公司（山東濟(jì)南 250101）

2.中國(guó)石油天然氣管道工程有限公司（河北廊坊 065000）

0 引言

山東地區(qū)某成品油長(zhǎng)輸管道設(shè)計(jì)壓力8.0 MPa，管徑219.1 mm，管道采用L245MB 鋼級(jí)鋼管，外防腐采用三層PE。管道附近有當(dāng)?shù)氐V務(wù)局計(jì)劃開采的3308綜放工作面，其中管道中線距離計(jì)劃開采工作面最近位置約316 m。根據(jù)礦區(qū)的開采計(jì)劃，3308工作面擬于2021年7月底開始回采，煤層開采深度580.7 m，平均開采厚度8.1 m；根據(jù)礦產(chǎn)開采分析報(bào)告，回采初期采動(dòng)影響會(huì)對(duì)本成品油長(zhǎng)輸管道造成一定的運(yùn)營(yíng)安全風(fēng)險(xiǎn)[1]。

采用概率積分法對(duì)管道地表位移進(jìn)行估算[2-3]，并將估算值作為位移荷載作用到在役管道節(jié)點(diǎn)中，采用有限元的方法，對(duì)管道在內(nèi)壓、溫差和地表沉降組合工況下的應(yīng)力狀態(tài)進(jìn)行數(shù)值計(jì)算[4]，計(jì)算得到管道應(yīng)力集中發(fā)生在彎頭部位，最大組合應(yīng)力值為314 MPa，超過了《輸油管道工程設(shè)計(jì)規(guī)范》的管材極限允許極限強(qiáng)度220.5 MPa，如圖1 所示,如果不采取任何措施，管道可能處于失效狀態(tài)[5]。

圖1 在役管道有限元計(jì)算應(yīng)力分布圖

項(xiàng)目前期開展了開挖應(yīng)力釋放和應(yīng)力監(jiān)測(cè)等管道保護(hù)措施，在管道周邊設(shè)置了地表位移GNSS監(jiān)測(cè)點(diǎn)3 個(gè)用于監(jiān)測(cè)煤礦開采期間的地表變形情況[6]，其中監(jiān)測(cè)點(diǎn)3 距離開挖面最近，監(jiān)測(cè)點(diǎn)1 距離開挖面距離最遠(yuǎn)，平面布置如圖2所示。

圖2 管道中線、開采工作面及GNSS監(jiān)測(cè)點(diǎn)布置平面圖

根據(jù)開采進(jìn)度工作面每周開采14～19 m，在2021年12月1日，煤礦開采工作面距離管道監(jiān)測(cè)點(diǎn)距離最近，此時(shí)開挖工作面距離管道彎頭處約316 m，如圖3所示。

圖3 采完工作面與管道位置關(guān)系

1 監(jiān)測(cè)方法選取

管道為線性工程，而煤礦采空區(qū)為面狀塌陷區(qū)。為判定采空區(qū)地表沉降對(duì)在役管道的影響，推薦采用多種監(jiān)測(cè)方法，通過對(duì)場(chǎng)區(qū)關(guān)鍵監(jiān)測(cè)點(diǎn)、管道中線及周圍區(qū)域的總體沉降規(guī)律分析，從而對(duì)地表變形穩(wěn)定做出準(zhǔn)確結(jié)論。工程推薦采用GNSS和InSAR兩種監(jiān)測(cè)方法，綜合判定地表變形規(guī)律，這是由兩種監(jiān)測(cè)方法的特點(diǎn)決定的。

1）從精度上，GNSS 測(cè)量精度一般能達(dá)到±5 mm+0.5×10-6D（D為距離參考點(diǎn)的距離）的精度，InSAR 結(jié)果為毫米級(jí)甚至是厘米級(jí)的精度，采用的GNSS測(cè)量精度要高于InSAR監(jiān)測(cè)結(jié)果[7-8]。

2）從監(jiān)測(cè)范圍上，GNSS 數(shù)據(jù)僅能針對(duì)本位置（精確為一點(diǎn)）的時(shí)程曲線；而遙感數(shù)據(jù)可以從監(jiān)測(cè)平面上總體觀測(cè)研究區(qū)域的特性，對(duì)區(qū)域的整體變形進(jìn)行定性和定量分析[9]。

3）從監(jiān)測(cè)頻率上，GNSS數(shù)據(jù)多采用無線傳輸實(shí)時(shí)傳輸，且監(jiān)測(cè)頻率可網(wǎng)絡(luò)調(diào)整，可以定制為密集的監(jiān)測(cè)頻次，滿足變形劇烈區(qū)域的監(jiān)測(cè)要求；而In-SAR 監(jiān)測(cè)頻率依賴所用數(shù)據(jù)的監(jiān)測(cè)周期[10-11]，例如哨兵數(shù)據(jù)的監(jiān)測(cè)周期為12天。

4）從監(jiān)測(cè)時(shí)間上，GNSS數(shù)據(jù)只能在設(shè)備安裝測(cè)試之后才能正常使用，設(shè)備一旦毀壞數(shù)據(jù)將丟失；而InSAR 數(shù)據(jù)具有很強(qiáng)的可追溯性，可以從監(jiān)測(cè)歷史上分析事件過程。

可見兩種監(jiān)測(cè)方法各有優(yōu)勢(shì)，又各有不足，通過采用兩種監(jiān)測(cè)方法結(jié)果的比對(duì)和判別，從而綜合得到場(chǎng)區(qū)的地表變形規(guī)律。

2 GNSS監(jiān)測(cè)

2.1 沉降監(jiān)測(cè)結(jié)果

項(xiàng)目在管道中線附近設(shè)置地表位移監(jiān)測(cè)點(diǎn)3個(gè)，監(jiān)測(cè)日期自2021 年8 月5 日至2022 年6 月1 日監(jiān)測(cè)的地表沉降值結(jié)果，如圖4所示。

圖4 GNSS監(jiān)測(cè)點(diǎn)累積地表沉降值

2.2 GNSS結(jié)果分析

由于3個(gè)地表沉降監(jiān)測(cè)點(diǎn)距離管道很近，3個(gè)點(diǎn)的監(jiān)測(cè)結(jié)果基本代表了管道彎頭附近的地表沉降變形情況，數(shù)據(jù)解析得出如下規(guī)律：

1）根據(jù)礦上反饋信息，工作面開挖在7 月底開始工作，GNSS 點(diǎn)1～點(diǎn)3 反應(yīng)出來地表沉降在9 月中旬左右出現(xiàn)第一次劇烈變化，說明在開采初期地表沉降顯現(xiàn)滯后于工作面開挖工作，滯后時(shí)間約1.5個(gè)月。

2）在2021 年12 月1 日煤礦工作面距離管道彎頭最近，此時(shí)距離監(jiān)測(cè)點(diǎn)距離也是最近的，此時(shí)GNSS監(jiān)測(cè)點(diǎn)1～點(diǎn)3開始發(fā)生第二次劇烈沉降。

3）隨著工作面繼續(xù)開挖，GNSS 監(jiān)測(cè)點(diǎn)1～點(diǎn)3持續(xù)增大，地表出現(xiàn)劇烈的變形階段，變形周期超過6個(gè)月，最大累積沉降值為25～60 mm不等。

4）從不同監(jiān)測(cè)點(diǎn)的變形趨勢(shì)上看，GNSS點(diǎn)1距離開挖面最遠(yuǎn)，從沉降絕對(duì)值上與GNSS 點(diǎn)2 和點(diǎn)3相比，變形絕對(duì)值較小，而且沉降值的收斂性滯后于其他監(jiān)測(cè)點(diǎn)。

3 InSAR監(jiān)測(cè)

本項(xiàng)目所用SAR數(shù)據(jù)為哨兵數(shù)據(jù)，時(shí)間跨度為2021 年8 月至2022 年5 月，周期為12 天，共23 景。SAR數(shù)據(jù)中包含了幅度信息和相位信息，本項(xiàng)目利用幅度信息獲取了沉降湖的變化，利用相位信息獲取了地表形變情況，具體如下。

3.1 沉降湖擴(kuò)展

基于時(shí)序強(qiáng)度影像發(fā)現(xiàn)出現(xiàn)沉降湖，并且隨著工作面開挖，湖的形態(tài)發(fā)生明顯變化，湖的現(xiàn)場(chǎng)沉降照片如圖5所示，形態(tài)的時(shí)間變化如圖6所示。

圖5 地面形成沉降湖泊

圖6 沉降湖形態(tài)變化圖

根據(jù)圖6 所示，沉降湖的外邊沿處于3308 開挖面中間位置，且隨著工作面的延伸，沉降湖總面積逐漸增大，且邊線向東南方向拓展，這與煤礦的開挖方向基本一致。沉降湖的形態(tài)變化可代表了整體區(qū)域的地表沉降區(qū)域的變化形態(tài)，表面地表塌陷區(qū)域在時(shí)間和空間上與地下煤礦工作面的開挖狀態(tài)是一致的。

3.2 形變結(jié)果

本項(xiàng)目采用SBAS 方法對(duì)時(shí)序數(shù)據(jù)進(jìn)行處理得到的累積形變結(jié)果，并選取沉降區(qū)域及關(guān)鍵點(diǎn)如圖7 所示，并在管道中線位置選取關(guān)鍵點(diǎn)1～關(guān)鍵點(diǎn)5 進(jìn)行分析。分析InSAR 關(guān)鍵點(diǎn)位移時(shí)程曲線如圖8 所示?；陉P(guān)鍵點(diǎn)的時(shí)程曲線結(jié)果分析，得出如下結(jié)論：

圖7 沉降區(qū)域及關(guān)鍵點(diǎn)選取

圖8 InSAR關(guān)鍵點(diǎn)位移時(shí)程曲線

1）距離沉降區(qū)域最近的關(guān)鍵點(diǎn)3表現(xiàn)出最大的位移值，最大位移值52 mm；其次是關(guān)鍵點(diǎn)兩側(cè)的關(guān)鍵點(diǎn)4 和關(guān)鍵點(diǎn)2，最大位移值在25～28 mm；變形量最小的是遠(yuǎn)離工作面的關(guān)鍵點(diǎn)5和關(guān)鍵點(diǎn)1，最大位移值在15 mm之內(nèi)。

2）從時(shí)間上來看，距離開挖面最近的關(guān)鍵點(diǎn)3在2021年10月初大幅度沉降，這與GNSS監(jiān)測(cè)點(diǎn)數(shù)據(jù)基本一致；而附近的關(guān)鍵點(diǎn)4 和關(guān)鍵點(diǎn)2 則表現(xiàn)為較平緩的變形；距離工作面較遠(yuǎn)的關(guān)鍵點(diǎn)1 和關(guān)鍵點(diǎn)5變形隨時(shí)間變化不明顯。

3）關(guān)鍵點(diǎn)1 距離開挖工作面最遠(yuǎn)，隨著時(shí)間發(fā)展甚至表現(xiàn)為小幅度的地表抬升，可能的原因是由于地下有剛度較大的地層作用，導(dǎo)致距離較遠(yuǎn)處的地層翹起所致。

4）從關(guān)鍵點(diǎn)2～關(guān)鍵點(diǎn)5的時(shí)程曲線看，地表在2022 年4 月初表現(xiàn)出地層的穩(wěn)定或位移的恢復(fù)（向0點(diǎn)位移方向發(fā)展），地層漸趨穩(wěn)定。

4 結(jié)論

通過對(duì)GNSS 和InSAR 兩種監(jiān)測(cè)方法的數(shù)據(jù)分析和比對(duì)，得出如下結(jié)論：

1）從上部地表與下部開挖關(guān)系上來看，地表監(jiān)測(cè)點(diǎn)的沉降通常滯后于地下開采工作，通常地表沉降滯后下部采空約1.5個(gè)月。

2）從空間關(guān)系上來看，距離開挖面最近的觀測(cè)點(diǎn)在2021 年10 月發(fā)生最大幅度沉降值（最大值約60 mm），而附近較遠(yuǎn)的觀測(cè)點(diǎn)表現(xiàn)為變形平緩。

3）從沉降發(fā)展方向來看，地表沉降發(fā)生與下部工作面開挖方向是一致的，都為從西北向東南的發(fā)展方向，表現(xiàn)為西北側(cè)逐漸穩(wěn)定，東南側(cè)沉降繼續(xù)增大。

4）從地表穩(wěn)定時(shí)間來看，在2022 年4 月管道上所有觀測(cè)點(diǎn)監(jiān)測(cè)值穩(wěn)定或有很小幅度變化，表明附近地層逐漸穩(wěn)定，具備下一步施工的條件。

地表沉降作為一種典型的地質(zhì)災(zāi)害形式，只有通過多種監(jiān)測(cè)手段掌握其變形規(guī)律，才能為后期管道安全性校核和地貌恢復(fù)提供有利的科學(xué)依據(jù)。