郭春生

(上海巖土工程勘察設(shè)計(jì)研究院有限公司，上海 200438)

0 引言

由于盾構(gòu)法通縫隧道的自身結(jié)構(gòu)、周邊水土壓力、行車動(dòng)荷載等因素的影響，隨著時(shí)間的推移，盾構(gòu)隧道大多產(chǎn)生“橫鴨蛋”形變形，為了了解這種變形量的大小，上海地區(qū)多采用全斷面收斂測(cè)量的方法.全斷面掃描收斂測(cè)量作為評(píng)價(jià)圓形盾構(gòu)隧道斷面變形情況的重要指標(biāo)，與沉降、平面位移、滲漏水調(diào)查等方法一起，從不同側(cè)面評(píng)價(jià)盾構(gòu)隧道健康狀況［1］.

目前，全斷面掃描收斂測(cè)量數(shù)據(jù)采集一般采用全站儀、附以機(jī)載程序、沿隧道剖面按20cm左右的步長(zhǎng)自動(dòng)逐點(diǎn)采集，數(shù)據(jù)處理一般采用以橢圓為數(shù)據(jù)模型進(jìn)行全斷面整體擬合，以擬合的橢圓長(zhǎng)軸和短軸與設(shè)計(jì)半徑的較差的大小來(lái)評(píng)價(jià)變形量［2］.本文針對(duì)橢圓模型的不足，提出多弧段擬合法的數(shù)據(jù)處理改進(jìn).

1 橢圓模型擬合算法及不足

地鐵隧道的拱頂下沉、兩側(cè)外擴(kuò)的形變趨勢(shì)，剖面呈近似的橢圓;即使隧道未發(fā)生變形，以橢圓擬合也能得到長(zhǎng)、短相等的特殊橢圓(標(biāo)準(zhǔn)圓).因此，前些年多家測(cè)繪同行單位建立了基于橢圓模型的全斷面擬合算法［3-4］.

圖1基于測(cè)站儀器中心為原點(diǎn)、鉛垂向上為Y軸、剖面內(nèi)水平為X軸的數(shù)據(jù)采集坐標(biāo)系內(nèi)，坐標(biāo)與觀測(cè)量(距離D、垂直角ele)之間的關(guān)系為:

以橢圓圓心為原點(diǎn)、長(zhǎng)軸為X、短軸為Y的坐標(biāo)系下，各極坐標(biāo)點(diǎn)的坐標(biāo)與數(shù)據(jù)采集坐標(biāo)系的轉(zhuǎn)換關(guān)系:

圖1 橢圓模型基于圓心的坐標(biāo)系

亦即(X0+RX)TΛ(X0+RX)-1=0

以所有采集點(diǎn)為觀測(cè)值，以橢圓的參數(shù)a、b及坐標(biāo)轉(zhuǎn)角參數(shù)x0、y0、α為待求量，建立橢圓誤差方程如下:

按最小二乘法擬合，即可求解出圓心坐標(biāo)XO、長(zhǎng)半軸、短半軸及其相對(duì)數(shù)據(jù)采集系的旋轉(zhuǎn)角，進(jìn)而可基于擬合的橢圓計(jì)算剖面上各位置的變形量.

長(zhǎng)期測(cè)量數(shù)據(jù)處理實(shí)踐中發(fā)現(xiàn)，上述處理方式存在以下不足:

(1)觀測(cè)數(shù)據(jù)剔除粗差的閥值設(shè)置困難.擬合前需要剔除管線、手孔等非管片環(huán)上的粗差點(diǎn)，一般以實(shí)測(cè)點(diǎn)偏離擬合橢圓距離為閥值.粗差剔除閥值不能設(shè)得太大，其起不到自動(dòng)數(shù)據(jù)處理的效果，很多點(diǎn)要手工刪除.另一方面，由于盾構(gòu)環(huán)結(jié)構(gòu)特性，變形較大時(shí)，剖面內(nèi)變形通常表現(xiàn)為相鄰管片錯(cuò)臺(tái)、相對(duì)轉(zhuǎn)角，其幾何特性實(shí)際已不是橢圓，較大時(shí)與橢圓偏離2 cm以上，因此粗差剔除閥值不能設(shè)得太小.斷面內(nèi)測(cè)距精度一般能達(dá)到2 mm以內(nèi)，理想的閥值宜設(shè)置為4 mm左右，但基于橢圓模型，對(duì)于變形大的斷面這個(gè)閥值將剔除很多有效觀測(cè)點(diǎn)，顯然不合適，因此觀測(cè)數(shù)據(jù)自動(dòng)剔除粗差的閥值設(shè)置困難.

(2)變形較大時(shí)與水平直徑的直接測(cè)量值相關(guān)較大，不能提供相鄰管片的錯(cuò)臺(tái)、旋轉(zhuǎn)等信息.基于橢圓模型，總體變形較小時(shí)求解的變形量與CAD中展點(diǎn)量測(cè)的變形量基本一致;對(duì)局部變形量較大的區(qū)段，橢圓模型把隧道剖面描述為連續(xù)的曲線，不能反映現(xiàn)狀管片的錯(cuò)臺(tái)、旋轉(zhuǎn)情況，求解的水平直徑常常與實(shí)際值相差超過10 mm，損失了較大的測(cè)量精度.

2 多弧段擬合法

大量收斂測(cè)量斷面測(cè)量數(shù)據(jù)表明，由于封頂塊周圍土壓力的分布特點(diǎn)，相對(duì)拱底塊而言，封頂塊在變形過程中容易發(fā)生豎向位移.在其變形發(fā)展過程中，各管片之間的關(guān)系(見圖2).由于隧道對(duì)稱，取一半隧道剖面考慮，如果封頂塊向下移動(dòng)，隧道管片將繞著圖2所示的A點(diǎn)、B點(diǎn)、C點(diǎn)轉(zhuǎn)動(dòng)(A'、B'、C'為對(duì)稱點(diǎn)，假定對(duì)稱位移).在A點(diǎn)處，隧道接縫外部張開，B點(diǎn)處，隧道接縫內(nèi)部張開，C點(diǎn)處，隧道接縫內(nèi)部張開.由于管片鋼度相對(duì)接口部位強(qiáng)度較大，同一管片的形態(tài)變形較小.

通過對(duì)關(guān)鍵成功領(lǐng)域進(jìn)行解析與細(xì)化，明確了每個(gè)關(guān)鍵成功領(lǐng)域包含的內(nèi)容、最主要因素、達(dá)成每一個(gè)成功領(lǐng)域的標(biāo)準(zhǔn)。如質(zhì)量安全的關(guān)鍵績(jī)效要素為醫(yī)療質(zhì)量和醫(yī)療安全，運(yùn)營(yíng)管理為資產(chǎn)效益、成本控制和管理運(yùn)營(yíng)，創(chuàng)新發(fā)展為學(xué)科人才、新醫(yī)療新技術(shù)、教學(xué)成果、科學(xué)研究。

圖2 通縫隧道變形特點(diǎn)

由于基于橢圓模型的算法存在上述不足，本文提出多弧段斷面擬合的方法.本文將圓形盾構(gòu)隧道看作由若干個(gè)弧段拼接而成，通過弧段各自的擬合，并把擬合后的各弧段連接成一不規(guī)則剖面，與設(shè)計(jì)斷面比較，從而計(jì)算斷面上各位置的變形量、橫向及豎向等典型位置變形量、相鄰弧段間的相對(duì)錯(cuò)臺(tái)量與相對(duì)旋轉(zhuǎn)量，實(shí)現(xiàn)對(duì)圓形盾構(gòu)隧道的精確變形測(cè)量.觀測(cè)和數(shù)據(jù)處理順序如下.

(1)測(cè)量時(shí)，采集觀測(cè)數(shù)據(jù)時(shí)應(yīng)保證每個(gè)管片不少于6點(diǎn)，并加測(cè)相鄰管片拼縫的信息(見圖3).

圖3 數(shù)據(jù)采集時(shí)加測(cè)拼裝縫的信息

(2)把觀測(cè)數(shù)據(jù)以所述管片間的拼接縫作為分界點(diǎn)，將圓形盾構(gòu)隧道斷面劃分為若干弧段;按弧段模型逐一對(duì)各管片采集的數(shù)據(jù)點(diǎn)進(jìn)行最小二乘法擬合［5-7］:

式中:(x0)T是圓心坐標(biāo)，R是圓的半徑對(duì)于測(cè)得點(diǎn)(xi)T，列出誤差方程:

式中:vi相當(dāng)于i點(diǎn)與圓周的距離.

圓心坐標(biāo)的迭代初值可以取為所有測(cè)定點(diǎn)的坐標(biāo)均值，半徑初值可取為任意正數(shù).

(3)對(duì)5個(gè)弧段逐一作為一圓形弧段模型進(jìn)行擬合，計(jì)算每個(gè)采集點(diǎn)與擬合園心的距離Ri，并與擬合半徑比較計(jì)算偏離度δi，以之為統(tǒng)計(jì)量、以置信度95%(大于2倍中誤差)為閥值，剔除偏離度較大的粗差點(diǎn)后進(jìn)行再次擬合.實(shí)踐證明剔除粗差的閥值一般可設(shè)為4 mm，能較好地自動(dòng)完成粗差點(diǎn)剔除;

(4)對(duì)拱底塊則采用兩側(cè)最下部的觀測(cè)點(diǎn)及擬定半徑計(jì)算理論弧段，所有弧段拼接成一個(gè)全斷面;

(5)在由上述5個(gè)實(shí)測(cè)數(shù)據(jù)擬合的弧段、及其拱底塊按理論計(jì)算的弧段組成的剖面內(nèi)，以拱底塊的圓心為原點(diǎn)建立剖面坐標(biāo)系，從而可得到:

①各觀測(cè)點(diǎn)至原點(diǎn)的距離與設(shè)計(jì)半徑的較差可計(jì)算隧道圓周各角度處的徑向變形量;

②根據(jù)6個(gè)拼接縫觀測(cè)點(diǎn)可確定各弧段的分界點(diǎn)，分界點(diǎn)兩側(cè)的兩個(gè)弧段與分界點(diǎn)到原點(diǎn)連線平面內(nèi)相交可求兩個(gè)交點(diǎn)，兩個(gè)交點(diǎn)的距離即為相鄰塊的錯(cuò)臺(tái)量L;

③上步中兩個(gè)交點(diǎn)處兩段弧的切線之間的夾角即為相鄰弧段相對(duì)旋轉(zhuǎn)角a;

④計(jì)算弧段最大橫向距離可求解水平向直徑、計(jì)算弧段最大豎向距離可求解豎向直徑.

圖4 相鄰管片的旋轉(zhuǎn)角、錯(cuò)臺(tái)量計(jì)算

3 算例

某圓形盾構(gòu)隧道的變形量在CAD中展點(diǎn)后，其橫徑變形為0.141，其縱徑變形為 -0.167，針對(duì)此圓形盾構(gòu)隧道分別采用橢圓擬合算法(1)以及本文擬合方法(2)來(lái)計(jì)算其變形程度.

(1)采用傳統(tǒng)橢圓模型解算，按每15°計(jì)算圓周上各點(diǎn)與圓心的距離跟理論半徑比較、以向下為角度零方向、順時(shí)針為角度正方向展開(見圖5).

根據(jù)上圖，采用傳統(tǒng)的橢圓模型時(shí):可解算出橫向變形(90°處變形量+270°處的變形量)、豎向變形(180°處變形量)，且變形量與CAD展點(diǎn)成果相差約6 mm;但不能解算相鄰管片的錯(cuò)臺(tái)量、旋轉(zhuǎn)角，使其無(wú)法全面地測(cè)量出該圓形盾構(gòu)隧道的變形情況.

圖5 基于橢圓模型的斷面變形展開圖

(2)按照本文方法，根據(jù)實(shí)測(cè)的數(shù)據(jù)，可解算出:

①5個(gè)弧段(代表上述的管片2-6)的參數(shù)(見表1).

表1 各擬合弧段參數(shù)表

②根據(jù)表1中各個(gè)弧段的擬合數(shù)據(jù)，計(jì)算相鄰弧段的錯(cuò)臺(tái)量、旋轉(zhuǎn)角(見表2)，其中后管片相對(duì)前管片向外錯(cuò)臺(tái)為正，反之為負(fù);旋轉(zhuǎn)角以外凸為正，反之為負(fù).

表2 相鄰管片錯(cuò)臺(tái)量、旋轉(zhuǎn)角成果表

③按每15°計(jì)算圓周上各點(diǎn)與圓心的距離跟理論半徑比較，同時(shí)納入管縫位置點(diǎn)，同樣以向下為 角度零方向、順時(shí)針為角度正方向展開(見圖5).

圖5 基于多弧段擬合法的斷面變形展開圖

4 結(jié)語(yǔ)

實(shí)際數(shù)據(jù)處理成果表明，采用多弧段擬合法與采用傳統(tǒng)橢圓模型解算成果圖比較，更明顯說(shuō)明每片管片內(nèi)部線性較好，變形、錯(cuò)臺(tái)主要發(fā)生在管片與管片之間的接口部位.采用本文方法的成果與實(shí)際CAD展點(diǎn)圖比較，除突出或凹入隧道表面的粗差點(diǎn)以外，各觀測(cè)點(diǎn)偏離擬合園的距離基本不超過2 mm，并且能夠通過橫徑變形、縱徑變形、相對(duì)錯(cuò)臺(tái)量、相對(duì)旋轉(zhuǎn)角這四個(gè)方面綜合評(píng)定圓形盾構(gòu)隧道多弧段斷面的變形程度.

［1］華錫生，黃 騰.精密工程測(cè)量技術(shù)及應(yīng)用［M］.南京:河海大學(xué)出版社，2005.

［2］王解先，季凱敏.工業(yè)測(cè)量擬合［M］.北京:測(cè)繪出版社，2008.

［3］孫 昊，姚連璧，謝義林.隧道斷面收斂測(cè)量與數(shù)據(jù)處理［J］.工程勘察，2010(11):44-46.

［4］黃茂華，義林謝.隧道收斂測(cè)量方法研究［J］.北京測(cè)繪，2009(4):14-16.

［5］胡龍華，王繼宇，劉雙成.全站儀無(wú)定向自由設(shè)站隧道收斂測(cè)量方法與精度分析［J］.華北科技學(xué)院學(xué)報(bào)，2010，7(3):40-43.

［6］宗玉玉.上海地鐵隧道長(zhǎng)期收斂測(cè)量方法研究［J］.淮海工學(xué)院學(xué)報(bào):自然科學(xué)版，2011(S1):132-134.

［7］趙兵帥，黃 騰，歐 樂.基于橢圓擬合的隧道斷面監(jiān)測(cè)及其應(yīng)用［J］.水利與建筑工程學(xué)報(bào)，2013(1):126-129.