本發(fā)明涉及應(yīng)用磨損比耗指數(shù)預(yù)測滾刀磨損量的分層求和法，屬于地下隧道工程的施工技術(shù)領(lǐng)域。

背景技術(shù)：

盾構(gòu)掘進高強度復(fù)合地層，工程面臨的最大風(fēng)險是掘進速率緩慢或掘進困難，且刀具磨損嚴(yán)重。刀具磨損已經(jīng)成為一個影響工程進度及成本的關(guān)鍵問題。刀具的完善性及刀具的磨損狀態(tài)判斷是盾構(gòu)掘進中的關(guān)鍵因素，如何準(zhǔn)確地定量預(yù)測滾刀的磨損量并判斷磨損狀態(tài)，對盾構(gòu)高效掘進至關(guān)重要。準(zhǔn)確預(yù)測盤形滾刀磨損量，對提高盾構(gòu)施工效益有重要意義。目前國內(nèi)外滾刀磨損預(yù)測方法可歸為兩大技術(shù)路徑：力學(xué)分析法及經(jīng)驗預(yù)測法。力學(xué)模型中的關(guān)鍵系數(shù)受試驗條件及巖體性質(zhì)影響顯著，在復(fù)合地層滾刀磨損預(yù)測應(yīng)用中存在局限。經(jīng)驗預(yù)測法，對地層分布在滾刀軌跡、掘進方向上的區(qū)分不足，影響其在巖層分布多變的復(fù)合地層中的適用性。

技術(shù)實現(xiàn)要素：

本發(fā)明針對現(xiàn)有技術(shù)的不足，提出了應(yīng)用磨損比耗指數(shù)預(yù)測滾刀磨損量的分層求和法。

本發(fā)明是應(yīng)用磨損比耗指數(shù)預(yù)測滾刀磨損量的分層求和法，它包括以下步驟：

步驟一：統(tǒng)計由盾構(gòu)機自動記錄的盾構(gòu)掘進每環(huán)的掘進參數(shù)平均值；

步驟二：進行滾刀刀圈的磨損量量測，分別記錄在每次開倉檢查或換刀時的各滾刀磨損增加量δ^ωi-l-j，滾刀破巖地層為復(fù)合地層時記為δ^δi-l-j；

步驟三：計算并記錄下掘進過程中不同滾刀掘進每一環(huán)中各地層的破巖體積；

步驟四：計算磨損比耗指數(shù)swi(specificwearindex)；

步驟五：將位置各異的滾刀的磨損增加量δ^ωi-l-j按掘進地層重新分類統(tǒng)計，作為回歸分析樣本，將同種地層中的swi所對應(yīng)的掘進參數(shù)集合、swi所對應(yīng)的掘進參數(shù)的若干次冪作為元素集，進行元素集間笛卡爾乘積運算，得到屬于同種地層swi的單項變量及二元復(fù)合變量集合，作為多元回歸分析的自變量集；

步驟六：對滾刀掘進各種均質(zhì)地層的swi樣本及相應(yīng)自變量集元素進行多元逐步線性回歸運算，得到分別獨立屬于各種均質(zhì)地層的swi與掘進參數(shù)間的高階多項式函數(shù)關(guān)系，按照共性提升原則，將屬于不同均質(zhì)地層的swi與掘進參數(shù)間的高階多項式函數(shù)關(guān)系中相同的自變量抽出來，將它們線性組合，作為預(yù)測均質(zhì)地層滾刀swi的通式，即swi預(yù)測方程；

步驟七：計算屬于不同地層的swi預(yù)測方程諸系數(shù)的離均差率mdra、mdrb、mdrc、mdrk；

步驟八：根據(jù)勘察報告，了解掘進地層的分布情況；如果掘進地層是復(fù)合地層，也就是說掌子面由若干個地層組合而成，就根據(jù)地層參數(shù)，從而確定了掌子面里面每一個地層的bq值或[bq]值；

步驟九：確定屬于掌子面里面各地層swi方程系數(shù)；即根據(jù)勘察報告計算得出地層bq值或[bq]值，在橫坐標(biāo)上找到該地層bq值或[bq]值，對應(yīng)的點并過該點作一條垂線，垂線分別和mdra與bq值或[bq]值的關(guān)系曲線、mdrb與bq值或[bq]值的關(guān)系曲線、mdrc與bq值或[bq]值的關(guān)系曲線、mdrk與bq值或[bq]值的關(guān)系曲線的4個交點的縱坐標(biāo)，就是屬于該地層的swi預(yù)測方程諸系數(shù)的離均差率mdra、mdrb、mdrc、mdrk，并在此基礎(chǔ)上確定屬于掌子面里面某一個地層swi方程系數(shù)a、b、c、k的具體值；

步驟十：提出一組掘進速率v和刀盤轉(zhuǎn)速n的試算值，代入到a、b、c、k為具體數(shù)值的swi預(yù)測方程中，就得到了屬于這一地層的、這一組v和n對應(yīng)的磨損比耗指數(shù)預(yù)測值swi^ω*；

步驟十一：提出與掘進速率v的試算值和刀盤轉(zhuǎn)速n的試算值相對應(yīng)的開倉環(huán)號預(yù)選值，再從勘察報告和地層剖面圖中計算出每一環(huán)里各地層的體積把它們和與地層對應(yīng)的swi預(yù)測值相乘，就得到了任意一環(huán)中各地層對應(yīng)的磨損增量，如下式：

步驟十二：

把一環(huán)中各地層對應(yīng)的磨損增量求和，就得到了滾刀在任意一環(huán)的磨損增量w為一環(huán)中地層種類數(shù)量，如下式：

步驟十三：計算兩次開倉之間掘進了q環(huán)(從第p+1環(huán)到第p+q環(huán))的磨損增量，如下式：

步驟十四：滾刀壽命內(nèi)，即從換上新刀到更換下來，中間開了t次倉，根據(jù)下式計算滾刀壽命內(nèi)的磨損量，用△表示，和磨損增量δ相區(qū)別：

步驟十五：盾構(gòu)機在一定刀盤扭矩和掘進速率之下掘進時，可以根據(jù)刀盤扭矩和掘進速率實時預(yù)測滾刀達到預(yù)警值所要掘進的長度或環(huán)數(shù)；當(dāng)選定開倉位置時，即在確定的掘進長度條件下，通過計算掘進速率試算值和刀盤扭矩試算值所對應(yīng)的滾刀壽命內(nèi)的刀盤磨損量將滾刀壽命內(nèi)的磨損量與預(yù)警值(如：20mm±1mm)做比較，如果預(yù)測的磨損量低于20mm，則滾刀利用效率太低；如果磨損量預(yù)測值大于20mm，則滾刀與圍巖間作用過于劇烈；這些情況需要重新選取v和n，直到得到滿足要求的v和n，并將滿足要求的v和n設(shè)定為施工參考值。

原理和優(yōu)勢：

針對復(fù)合地層在滾刀軌跡及掘進方向上的地層分布變化，通過逐刀統(tǒng)計軌跡各異的滾刀的“磨損比耗指數(shù)”，作為磨損量回歸分析目標(biāo)值。最大的好處是消除了滾刀安裝半徑對滾刀間磨損增加量差異的影響，將各滾刀的“磨損比耗指數(shù)”按滾刀掘進地層分類統(tǒng)計及回歸分析，得到分別屬于各種均質(zhì)地層的磨損比耗指數(shù)經(jīng)驗公式。提出適用于復(fù)合地層的應(yīng)用磨損比耗指數(shù)預(yù)測滾刀磨損量的分層求和法，在考慮滾刀軌跡差異對滾刀磨損量影響的同時提高了經(jīng)驗預(yù)測法的地層適用性。

首次提出基于磨損比耗指數(shù)回歸分析與滾刀破巖體積分層統(tǒng)計的應(yīng)用磨損比耗指數(shù)預(yù)測滾刀磨損量的分層求和法，實現(xiàn)了對刀盤上任意位置的正面滾刀及中心滾刀的磨損量及壽命的定量預(yù)測，有效地解決了復(fù)雜地層滾刀磨損量及壽命的定量預(yù)測的難題。應(yīng)用磨損比耗指數(shù)預(yù)測滾刀磨損量的分層求和法對滾刀磨損量與壽命具有良好的預(yù)測精度，豐富了滾刀磨損預(yù)測研究成果，有利于科學(xué)編制預(yù)算、合理安排施工進度、提高盾構(gòu)掘進效率，具有重要的理論意義與實用價值。

附圖說明

圖1滾刀破巖區(qū)域劃分截面示意圖；

圖2滾刀破巖空間立體示意圖；

圖3滾刀swi及破巖體積分層統(tǒng)計示意圖；

圖4j31-2地層swi實測值與回歸曲面圖；

圖5j31-3r地層swi實測值與回歸曲面圖；

圖6j31-3p地層swi實測值與回歸曲面圖；

圖7j31-3地層swi實測值與回歸曲面圖；

圖8不同地層swi預(yù)測方程諸系數(shù)的離均差率與地層bq值或[bq]值的關(guān)系示意圖；

圖9確定屬于掌子面里面某一個地層的swi方程系數(shù)a、b、c、k的具體值的方法；

圖10應(yīng)用磨損比耗指數(shù)預(yù)測滾刀磨損量的分層求和法流程圖；

具體實施方法

本發(fā)明應(yīng)用磨損比耗指數(shù)預(yù)測滾刀磨損量的分層求和法，包括以下步驟：

步驟一：

統(tǒng)計由盾構(gòu)機自動記錄的盾構(gòu)掘進每環(huán)的掘進參數(shù)平均值(vi-l-j-m,ni-l-j-m,fi-l-j-m,ti-l-j-m,pi-l-j-m)，其中，v為推進速率、n為刀盤轉(zhuǎn)速、f為有效推力、t為扭矩、p為土倉壓力，i為滾刀編號，l為滾刀更換次序，j為滾刀壽命內(nèi)的開倉量測次序，m為環(huán)號；

步驟二：

進行滾刀刀圈的磨損量量測，分別記錄在每次開倉檢查或換刀時的各滾刀磨損增加量δ^ωi-l-j，滾刀破巖地層為復(fù)合地層時記為δ^δi-l-j，ω為滾刀破巖地層序號；

步驟三：

計算并記錄下掘進過程中不同滾刀掘進每一環(huán)的破巖體積；

以圖1至圖3為例具體說明。

如圖1所示，圖中標(biāo)號1、2、3、4，分別代表四把不同的滾刀，在各自的軌跡上進行掘進。統(tǒng)計時：取相鄰滾刀外緣間距中點為相鄰滾刀破巖斷面邊界。

如圖2中，ω1、ω2、ω3分別為滾刀掘進三種不同地層。

如圖3中，ω1、ω2、ω3、ω4分別為滾刀掘進四種不同地層。

以滾刀為研究對象，滾刀在環(huán)形軌跡上掘進時，可能遇到不同的地層，滾刀在地層中掘進的空心柱體空間的體積是空心柱體中各地層體積之和。如圖3表格中所示，將滾刀在地層中掘進的空心柱體空間中的各地層體積分層統(tǒng)計。其中，如果某次開倉量測滾刀磨損量與上一次開倉量測滾刀磨損量之間，滾刀掘進破巖空心柱體空間中的地層種類多于一種，即滾刀掘進的是復(fù)合地層，由于滾刀磨損增加量δ^△i-l-j與該階段滾刀在復(fù)合地層中的破巖體積之和相對應(yīng)，無法確定滾刀磨損增加量δ^△i-l-j與該階段滾刀在復(fù)合地層中各地層相應(yīng)的破巖體積之間的對應(yīng)關(guān)系，無法計算得到swi實測值；如果某次開倉量測滾刀磨損量與上一次開倉量測滾刀磨損量之間，滾刀掘進的是均質(zhì)地層，這時滾刀磨損增加量δ^ωi-l-j與該階段滾刀掘進地層的類型和體積均明確對應(yīng)，可以計算得到swi實測值。

具體來說，如式(1)，以單把滾刀為研究對象，取軌跡相鄰的滾刀外緣間距中點為相鄰滾刀破巖斷面邊界，滾刀在環(huán)形軌跡上掘進一環(huán)，滾刀破巖產(chǎn)生的空心柱體空間的體積，即滾刀掘進每環(huán)的空心柱體的總體積vi-l-j-m為：

vi-l-j-m＝l×si-l-j-m(1)

按所遇到的不同的地層，計算并統(tǒng)計各滾刀掘進每環(huán)產(chǎn)生的空心柱體中各地層體積，記為v^ωi-l-j-m，且vi-l-j-m＝σv^ωi-l-j-m；其中：l為盾構(gòu)管片環(huán)寬，si-l-j-m為該滾刀破巖環(huán)形區(qū)域的面積；

步驟四：

進行滾刀的swi統(tǒng)計。各滾刀與刀盤旋轉(zhuǎn)中心之間的距離(滾刀安裝半徑)均不相同，導(dǎo)致掘進相同里程時各滾刀破巖體積不同。定義磨損比耗指數(shù)(specificwearindex，swi)，即單位破巖體積對應(yīng)的滾刀磨損增加量，如式(2)。

swi的意義在于：磨損增加量和半徑正相關(guān)，破巖體積和半徑正相關(guān)，兩者相除，得到swi，實現(xiàn)了對半徑的降次處理。所以，swi受半徑的影響大大小于磨損增加量，可以認(rèn)為“消除了滾刀安裝半徑對滾刀間磨損增加量差異的影響”。

滾刀是螺線前進的，但掘進速率變化太大，無法準(zhǔn)確計算出某一段里滾刀運動的螺線長度。所以，我們把破巖體積這個確定的量作為分母，方便分析與應(yīng)用。

步驟五：

最后，將位置各異的滾刀的磨損增加量按掘進地層重新分類統(tǒng)計，作為回歸分析樣本。將同種地層swi相應(yīng)的掘進參數(shù)集合及掘進參數(shù)若干次冪作為元素集，進行元素集間笛卡爾乘積運算，得到屬于同種地層swi的單項變量及二元復(fù)合變量集合，作為多元回歸分析的自變量集。

步驟六：

對滾刀掘進各種均質(zhì)地層的swi樣本及相應(yīng)自變量集元素進行多元逐步線性回歸運算，得到分別獨立屬于各種均質(zhì)地層的swi與掘進參數(shù)間的高階多項式函數(shù)關(guān)系。按照共性提升原則，將屬于各種均質(zhì)地層的swi與掘進參數(shù)間的高階多項式函數(shù)關(guān)系中相同的自變量抽出來，將它們線性組合成式(3)，作為預(yù)測均質(zhì)地層滾刀磨損速率的通式。

swi^ω*＝kv^an^b+c(3)

式中，v為掘進速率；n為刀盤轉(zhuǎn)速；a、b、c及k均為回歸方程系數(shù)，swi^ω*為swi^ω的預(yù)測值。該通式應(yīng)用于南京市寧高城際軌跡交通二期工程盾構(gòu)隧道項目，如圖4～7為不同地層swi預(yù)測方程計算結(jié)果，其中v/mm·min^-1表示掘進速率v的單位為mm·min^-1；n/r·min^-1表示刀盤轉(zhuǎn)速n的單位為r·min^-1。

如表1及表2，通過對比預(yù)測值和實測值，得到回歸方程系數(shù)及方差分析結(jié)果：因為自變量顯著性在各地層中均小于0.05且相關(guān)系數(shù)較大，所以式(3)具有統(tǒng)計意義，對于滾刀磨損比耗指數(shù)具有較好的預(yù)測效果。

表1均質(zhì)地層swi回歸計算系數(shù)

表2預(yù)測均質(zhì)地層滾刀磨損速率的通式的預(yù)測實例

步驟七：

如圖8，利用式(4)、式(5)、式(6)和式(7)，計算屬于不同地層swi預(yù)測方程諸系數(shù)的離均差率(mdra、mdrb、mdrc、mdrk)。

步驟八：

根據(jù)工程勘察報告，可以了解到掘進地層的分布情況。如果掘進地層是復(fù)合地層，也就是說掌子面由若干個地層組合而成，我們就根據(jù)地層參數(shù)如巖體完整性指數(shù)和巖石飽和單軸抗壓強度，從而確定了掌子面里面每一個地層的bq值(或[bq]值)。

步驟九：

確定屬于掌子面里面各地層swi方程系數(shù)。如圖9，根據(jù)勘察報告計算得出地層bq值(或[bq]值)，在橫坐標(biāo)上找到該地層bq值(或[bq]值)對應(yīng)的點并過該點作一條垂線，垂線分別和mdra與bq值(或[bq]值)的關(guān)系曲線、mdrb與bq值(或[bq]值)的關(guān)系曲線、mdrc與bq值(或[bq]值)的關(guān)系曲線、mdrk與bq值(或[bq]值)的關(guān)系曲線的4個交點的縱坐標(biāo)，就是屬于該地層的swi預(yù)測方程諸系數(shù)的離均差率(mdra、mdrb、mdrc、mdrk)。根據(jù)式(4)、式(5)、式(6)和式(7)，確定屬于掌子面里面某一個地層swi方程系數(shù)a、b、c、k的具體值。

對掌子面內(nèi)各地層均進行步驟九內(nèi)的運算，得到屬于掌子面里面各地層swi方程系數(shù)。

步驟十：

如圖10，提出一組掘進速率v和刀盤轉(zhuǎn)速n的試算值，代入到a、b、c、k為具體數(shù)值的swi預(yù)測方程中，就得到了屬于這一地層的、這一組v和n對應(yīng)的磨損比耗指數(shù)預(yù)測值swi^ω*。

步驟十一：

當(dāng)步驟一所述的掘進參數(shù)改變時，滾刀和圍巖間的相互作用發(fā)生改變，滾刀的磨損狀態(tài)隨之變化，即在滾刀換刀磨損量一定的前提下，掘進環(huán)數(shù)隨掘進參數(shù)的變化而變化。然而，更換刀具時掘進停止，在圍巖條件較差的地段可能存在塌方等隱患，換刀時盾構(gòu)機所處的環(huán)號不能隨意設(shè)置；為了盡可能提高滾刀利用效率，應(yīng)該在滾刀換刀前盡可能掘進更長的距離，且掘進速率要滿足工期要求。所以，可以通過預(yù)測并控制掘進參數(shù)，實現(xiàn)滾刀掘進環(huán)數(shù)的定量預(yù)測和控制。

如圖10，提出與掘進速率v的試算值和刀盤轉(zhuǎn)速n的試算值相對應(yīng)的開倉環(huán)號預(yù)選值，再從勘察報告和地層剖面圖中計算出每一環(huán)里各地層的體積把它們和與地層對應(yīng)的swi預(yù)測值相乘，就得到了任意一環(huán)中各地層對應(yīng)的磨損增量，如式(8)

步驟十二：

如式(9)，把一環(huán)中各地層對應(yīng)的磨損增量求和，就得到了滾刀在任意一環(huán)的磨損增量w為一環(huán)中地層種類數(shù)量。

步驟十三：

如圖4，兩次開倉之間掘進了q環(huán)(從第p+1環(huán)到第p+q環(huán))的磨損增量根據(jù)式(10)計算。

具體來說，公式(10)表示：當(dāng)p環(huán)時開倉，測量磨損量一次；當(dāng)從第p+1環(huán)開始掘進到第p+q環(huán)時，再開倉測量磨損量一次；兩次開倉間的磨損增量，是與第p+1環(huán)到第p+q環(huán)之間滾刀的破巖量相對應(yīng)的；

步驟十四：

如圖10，若在滾刀壽命內(nèi)，即從換上新刀到更換下來，中間開了t次倉，那這就是滾刀壽命內(nèi)的磨損量，用大寫的△表示，和磨損增量小寫的δ相區(qū)別。

步驟十五：

盾構(gòu)掘進施工中，掘進參數(shù)的選取和開倉環(huán)號的安排不是隨意性的。掘進速率過大或刀盤扭矩過大，表明滾刀與圍巖之間作用劇烈，有損壞滾刀、刀座、刀盤、刀盤驅(qū)動設(shè)備的隱患。在圍巖完整程度較低、地下水豐富的地層，盾構(gòu)停工開倉會導(dǎo)致地層異常沉降和掌子面坍塌，嚴(yán)重威脅盾構(gòu)機和操作人員的安全。所以，需要根據(jù)勘察報告，選定地層完整程度較高、地下水較少的區(qū)間安排開倉換刀作業(yè)。為了實現(xiàn)盾構(gòu)機掘進至預(yù)定開倉位置之前的磨損量不超過滾刀磨損量預(yù)警值，需要根據(jù)擬定的掘進距離來選擇合適的掘進參數(shù)。當(dāng)盾構(gòu)機在一定刀盤扭矩和掘進速率之下開展掘進時，可以根據(jù)刀盤扭矩和掘進速率實時預(yù)測滾刀達到預(yù)警值所要掘進的長度或環(huán)數(shù)。

具體應(yīng)用方法如下：

如圖10，滾刀壽命內(nèi)的磨損量達到預(yù)警值，應(yīng)該更換滾刀。如果預(yù)測的滾刀壽命內(nèi)的磨損量低于20mm±1mm，則滾刀利用效率太低；如果磨損量預(yù)測值大于20mm±1mm，則滾刀與圍巖間作用過于劇烈。

當(dāng)選定開倉位置時，即在某確定的掘進長度條件下，通過計算掘進速率試算值和刀盤扭矩試算值所對應(yīng)的刀盤磨損量將滾刀壽命內(nèi)的磨損量與預(yù)警值20mm做比較，如果預(yù)測的磨損量低于20mm，則滾刀利用效率太低；如果磨損量預(yù)測值大于20mm，則滾刀與圍巖間作用過于劇烈。這些情況需要重新選取v和n，直到得到滿足要求的v和n，并將滿足要求的v和n設(shè)定為施工參考值。

當(dāng)實際掘進時需要判定一定掘進速率和刀盤扭矩條件下的滾刀磨損狀態(tài)時，應(yīng)根據(jù)掘進速率和刀盤扭矩實測值，結(jié)合地層分布及地層參數(shù)，可以預(yù)測一定掘進速率和刀盤扭矩條件下滾刀掘進任意距離時的磨損量。當(dāng)接近20mm時，應(yīng)提前安排開倉換刀準(zhǔn)備工作，防止?jié)L刀損壞后更換不及時導(dǎo)致的刀座或刀盤損壞，提高掘進效率。

本專利方法已應(yīng)用于南京市寧高城際軌跡交通二期工程盾構(gòu)隧道項目，如表3，相對于實測值，8具中心滾刀及24具正面滾刀的磨損量預(yù)測值平均相對偏差多在20％～30％內(nèi)，而滾刀壽命預(yù)測值平均相對偏差均低于20％且多在7％～16％內(nèi)，表明swi回歸分析與復(fù)合地層滾刀磨損量分層求和法對磨損量及滾刀壽命均具有較高的預(yù)測精度，且其對滾刀壽命的預(yù)測更為準(zhǔn)確。

表3滾刀磨損量及壽命預(yù)測值相對偏差

表3中，dcw為磨損量預(yù)測值與實測值間相對偏差，根據(jù)式(12)計算；dcl為20mm滾刀更換磨損水平的滾刀壽命預(yù)測值與實測值，間相對偏差根據(jù)式(13)計算。