李 準 張 窮

(中鐵二院工程集團有限責任公司， 成都 610031)

我國西南地區(qū)長期處于地殼活動強烈地帶，地形切割陡峻，地質(zhì)構造復雜，素有中國“地質(zhì)博物館”之稱，在此區(qū)域施工存在巨大的特殊地質(zhì)災害風險[1]。因鐵路隧道定額的普遍性和通用性，未能全面考慮各類特殊復雜地質(zhì)災害引起的安全、資源投入變化情況，西南地區(qū)鐵路隧道施工一旦遇到大變形地質(zhì)災害，現(xiàn)行定額無法充分反映其安全、質(zhì)量、進度成本投入，致使項目投資不足，增加了建設單位、設計單位和施工單位設計變更的工作量和投資控制難度[2]。針對這一情況，本文開展了以新建麗江至香格里拉鐵路軟巖大變形為基礎的定額測定工作，在軟巖三臺階帶仰拱一次開挖工法條件下通過不同大變形段施工斷面(大變形Ⅰ、Ⅱ、ⅡA、ⅢA)參數(shù)分析，對不同變形條件下的資源消耗、成本投入情況進行研究，通過對常規(guī)定額水平以外的特殊復雜地質(zhì)條件下隧道施工資源消耗、成本投入情況等進行研究分析，有針對性的提出概預算編制調(diào)整意見或編制專項補充施工定額，為麗香鐵路投資控制及以后大麗攀、滇藏鐵路等類似重大工程設計概預算編制提供參考。

1 項目背景

新建麗江至香格里拉鐵路位于云南省西北部，南起大麗鐵路麗江車站，向北跨越金沙江，經(jīng)小中甸至香格里拉，全長139.652 km。線路連接大理至麗江鐵路，并通過該鐵路和廣大鐵路與成昆鐵路相聯(lián)，是我國中長期鐵路網(wǎng)規(guī)劃中西部路網(wǎng)的重要組成部分。

新建麗香鐵路地處歐亞板塊和印度洋板塊相互碰撞匯聚形成的青藏高原東南緣之川滇菱形斷塊的西部邊界斷裂帶(金沙江—中甸斷裂帶)內(nèi)，屬我國著名的南北向地震帶南段之滇西地震帶，地質(zhì)構造復雜，新構造運動強烈[3]。受大規(guī)模地質(zhì)構造的強烈擠壓影響，隧道通過之處應力集中度高、多期巖漿侵入的地質(zhì)環(huán)境下存在極高地應力。

2 浮動平均值

浮動平均值求解法是數(shù)理統(tǒng)計學在實際應用中的一種拓展[4]。該方法主要用于改善當數(shù)據(jù)樣本存在一定變動或?qū)崪y數(shù)據(jù)波動較大，存在較大誤差時，采用傳統(tǒng)平均值求解法的適用性和準確性。

傳統(tǒng)的數(shù)理統(tǒng)計方法是以總體分布為正態(tài)分布作為前提推導得出的[5]。從基本概念的推導過程中不難看出，在樣本總體呈正態(tài)分布的情況下，這些統(tǒng)計方法具有良好的穩(wěn)定性，得到結果具有足夠的準確性。但對于大多數(shù)工程應用的實際案例，正態(tài)分布的假設通常只能近似地滿足，有時甚至無法滿足?；谶@一類波動性較強的數(shù)據(jù)樣本，統(tǒng)計方法穩(wěn)健性的研究和應用具有很大的實際意義[6]。

由于在實際觀測中提取的數(shù)據(jù)量通常較大，很難完全甄別出所有的異常數(shù)據(jù)。而如果異常數(shù)據(jù)中有極大或極小值的存在，傳統(tǒng)的算術平均數(shù)的穩(wěn)定性就會降低。為增加平均數(shù)的準確性，通常有兩種途徑：一是采用有效的方法或?qū)?shù)據(jù)盡可能地縮小至更加有效的范圍，以篩選出數(shù)據(jù)中的異常值把它們直接剔除，然后進行傳統(tǒng)平均值的計算；二是設計新的統(tǒng)計方法使得樣本數(shù)據(jù)中的異常值不會對最終結果有過大的影響，一般基于絕對偏差和最小原則，或給中位數(shù)賦予更高的權重，使得樣本的集中趨勢和統(tǒng)計結果的穩(wěn)定性有所改善。

第一種方法作為數(shù)理統(tǒng)計學中的一個重要課題，長久以來一直受到科研工作者們的廣泛關注，針對擁有龐大數(shù)據(jù)量的樣本效果顯著。但是若樣本的數(shù)據(jù)量偏小，則無法利用其自身規(guī)律和趨勢找出異常值或縮小樣本數(shù)據(jù)的有效范圍，第一種方法就難以發(fā)揮作用。而且，即使可以正確地確定某些異常數(shù)據(jù)，在剔除異常數(shù)據(jù)后，如果剩余的樣本數(shù)量太小，同樣可能會導致最終結果的準確性降低。因此，考慮到此課題中實際工程應用背景，本文選用第二種思路對算術平均值進行改進。浮動平均法中考慮樣本中異常數(shù)據(jù)對傳統(tǒng)平均值的影響，引入?yún)?shù)θ，使得所求得的浮動平均值可根據(jù)實際情況(即異常數(shù)據(jù)在整個樣本中的分布和多寡情況)，選擇更加接近中位數(shù)或是算術平均值，能夠比較靈活地在兩者之間浮動切換。

假設數(shù)據(jù)樣本為：x1,x2, …,xn。

樣本數(shù)據(jù)的標準差可通過式(1)求得：

(1)

(2)

浮動平均法不僅依據(jù)樣本中的數(shù)據(jù)對距離算術平均值的遠近相應地賦予了權重，更進一步引入了參數(shù)θ，可根據(jù)數(shù)據(jù)處理的實際背景情況和需要，調(diào)整那些距離算數(shù)平均值較遠數(shù)據(jù)的權重。例如，在某些情況下，一些實測的數(shù)據(jù)雖然距離樣本算術平均值的距離較遠，但仍正確地反應了樣本的某種特性，需賦予與距離較近的數(shù)據(jù)相當?shù)臋嘀匾哉_地反饋出樣本的這種特性。但當我們使用穩(wěn)健平均值進行數(shù)據(jù)分析的時候，這些客觀上距離樣本算術平均值較遠的數(shù)據(jù)會被自然地分配上較低的權重，對最終得到的平均值產(chǎn)生誤導。因此，浮動平均法在工程實測數(shù)據(jù)分析的客觀性和正確性上均優(yōu)于穩(wěn)健平均法。

3 成果運用

3.1 研究樣本

課題以新建麗香鐵路中義、七達里、宗思隧道(設計有效斷面30.97 m2，凈空寬4.54 m，高6.75 m)為樣本進行數(shù)據(jù)采集，對工作條件、步驟相近相似的工序進行合并和細化，重點突出擬編制子目高地應力強烈擠壓變質(zhì)區(qū)大變形的特點[7]。通過確定開挖、出渣、鋼拱架制安、噴射混凝土(機械手)、超前小導管施作、注漿、樹脂錨桿施作、混凝土拌制澆筑運輸、監(jiān)控量測、通風管線路、大變形段拆換鋼拱架及鋼筋網(wǎng)、大變形圬工非爆鑿除及機械拆除等20余個工序的施工組織模型，并以浮動平均法和算數(shù)平均法對測定的數(shù)千個數(shù)據(jù)進行分別處理和對比，以驗證浮動平均法的正確性及優(yōu)越性。受篇幅所限，本文僅以其中中義隧道出口工區(qū)和2號橫洞工區(qū)基礎數(shù)據(jù)為例。

3.2 基于R語言的數(shù)據(jù)處理

由于浮動平均值需通過求解式(2)得到，為得到最優(yōu)Q值，需進行大數(shù)據(jù)測算，為提高工作效率，本文采用R語言編程進行求解。R語言是一款用于統(tǒng)計分析的語言和操作環(huán)境。與之相配套的RGui軟件集統(tǒng)計分析與圖形展示為一體，可自由編寫程序代碼，且完全開源，是統(tǒng)計領域一種非常實用的計算工具[8]。

最終，得到具體計算結果如表1、表2所示，其中參數(shù)θ=0.5。

表1 出口正東基礎延續(xù)時間的兩種平均值對比表

表2 2號橫道正洞基礎延續(xù)時間的兩種平均值對比表

3.3 分析結果

從上表中的各項結果可以看出，求得的浮動平均值與傳統(tǒng)算術平均值大致相等。在對平均值和中位數(shù)給予相等權重的前提條件下(θ=0.5)，浮動平均值與算術平均值之間的差值基本控制在2.5%以內(nèi)。首先，從結果中不難看出，所有數(shù)據(jù)的浮動平均值均與傳統(tǒng)算數(shù)平均值接近，體現(xiàn)出了浮動平均值的有效性；其次，由于浮動平均法將中位數(shù)考慮了進去，遠離中位數(shù)的數(shù)據(jù)對最終結果的影響降低，遠離樣本中心的異常數(shù)據(jù)對最終結果的影響被人為減小，從而使得所得到的結果比傳統(tǒng)算術平均值更加真實。

對于不同的實際工程背景，由于統(tǒng)計方式和工作條件不同，出現(xiàn)異常值的潛在可能情況各有不同。某些場景中，有些數(shù)據(jù)雖距離樣本中位數(shù)較遠，但卻真實地反應了該數(shù)據(jù)項目的某種特性，需給它們賦予與其他距離中位數(shù)較近的數(shù)據(jù)相同的權重，才能使最終求得的平均值能夠更好地體現(xiàn)該組樣本所對應的數(shù)據(jù)特性。因此，在這種情況下，參數(shù)θ的設置可以更加接近1。對于相反的情形，若想對接近中位數(shù)的數(shù)據(jù)賦予更大的權重，則可以將參數(shù)θ設置更加接近0。

4 結束語

傳統(tǒng)的算術平均值可有效地處理大部分較為連續(xù)且分布比較集中的樣本數(shù)據(jù)。但對于許多工程項目中所涉及到的樣本數(shù)據(jù)，由于其數(shù)據(jù)收集方式和收集環(huán)境存在一些難以克服的自然誤差，在處理這些樣本的時候，需尋找一種科學的數(shù)學途徑，有效地減少異常數(shù)據(jù)對最終平均值的影響。本文提出的浮動平均法在處理實際工程數(shù)據(jù)時，可根據(jù)樣本數(shù)據(jù)中異常數(shù)據(jù)分布位置距離樣本中位數(shù)的遠近，通過調(diào)節(jié)參數(shù)θ來有效降低異常數(shù)據(jù)對平均值產(chǎn)生的影響。綜上所述，該方法主要適用于已經(jīng)大致了解異常數(shù)據(jù)分布位置的工程樣本數(shù)據(jù)，根據(jù)所知異常數(shù)據(jù)分布信息調(diào)節(jié)參數(shù)θ求得的浮動平均值的準確性將優(yōu)于傳統(tǒng)算術平均值。