王　琛， 張金柱

(北京城建設(shè)計發(fā)展集團股份有限公司, 北京 100025))

?

基于數(shù)值仿真和巖土特性的反演方法

王琛， 張金柱

(北京城建設(shè)計發(fā)展集團股份有限公司, 北京 100025))

【摘要】文章通過以數(shù)值仿真作為參數(shù)反演工具，以巖土參數(shù)對環(huán)境變形的影響模式和參數(shù)測試數(shù)據(jù)的離散性作為反演分析基礎(chǔ)，建立了基于數(shù)值仿真和巖土特性的反演方法，并初步實現(xiàn)了該方法。實例表明，該方法具有較強的實用性，為工程反演提供了一種新思路。

【關(guān)鍵詞】參數(shù)反演;數(shù)值仿真;巖土特性

地鐵等地下工程要求進行信息化設(shè)計和施工，根據(jù)施工監(jiān)測結(jié)果修正設(shè)計，其重點和難點是參數(shù)反演。目前，工程反演方法主要有逆解法、直接法、圖譜法和智能反演法[1-4]。逆解法中由正分析方程反推得到逆方程，利用量測位移得到待定參數(shù)；直接法中把參數(shù)反演問題轉(zhuǎn)化為一個目標函數(shù)的尋優(yōu)問題，直接利用正分析的過程和格式，通過迭代最小誤差函數(shù)，逐次修正未知參數(shù)的試算值，直至獲得“最佳值”；圖譜法中以預先通過由數(shù)值計算得到的對應于各種不同彈性模量和初始地應力與位移的關(guān)系曲線，建立簡便的圖譜和圖表，并根據(jù)相似原理，由現(xiàn)場量測位移通過圖譜和圖表的圖解反推初始地應力和彈性模量，智能反演法是一種源于自然進化的全局搜索優(yōu)化算法，與直接法相近，只是求解速度有所提高。這些方法是參數(shù)反演研究的優(yōu)秀成果，但也存在不足之處。逆解法要求輸入和輸出之間有明確的函數(shù)關(guān)系，直接法和智能反演法計算工作量大，解的穩(wěn)定性差，圖譜法過于依賴地區(qū)經(jīng)驗。

本文在上述成果的基礎(chǔ)上，以數(shù)值仿真作為反演工具，以巖土參數(shù)對環(huán)境變形的影響模式和參數(shù)測試數(shù)據(jù)的離散性作為反演基礎(chǔ)，建立了基于數(shù)值仿真和巖土特性的反演方法。該方法為巖土工程參數(shù)反演提供了一條新思路，對于參數(shù)反演的理論研究和工程應用都具有一定的參考價值。

1本文反演方法綜述

1.1基本思想

(1)基于數(shù)值仿真建立巖土參數(shù)和環(huán)境變形之間的正向映射關(guān)系。即根據(jù)巖土參數(shù)數(shù)值計算得到環(huán)境變形。假定參數(shù)合理時，數(shù)值計算結(jié)果基本正確。

(2)基于巖土參數(shù)對變形的影響模式建立巖土參數(shù)和環(huán)境位移之間的反向映射關(guān)系。即根據(jù)正向映射關(guān)系歸納出參數(shù)對變形的影響模式，具體到個案工程反演中調(diào)整參數(shù)少量試算來標定影響模式中的系數(shù)，以得到針對個案工程的影響關(guān)系式。把現(xiàn)場監(jiān)測的環(huán)境位移結(jié)果代入上述關(guān)系式，就可以求得參數(shù)的反演值。

(3)基于巖土參數(shù)測試數(shù)據(jù)的離散性確定參數(shù)優(yōu)化的合理區(qū)間。通過對大樣本測試數(shù)據(jù)的統(tǒng)計，能夠得到參數(shù)真實值的分布區(qū)間，即離散區(qū)間。對于反演得到的參數(shù)，如果落于離散區(qū)間之內(nèi)，則取該值作為反演值；如果小于離散下界，則取離散下界值作為反演值。反之亦然。

1.2反演步驟

本文方法分前期研究和個案反演兩個階段。在前期研究階段需要開發(fā)參數(shù)反演中使用的數(shù)值仿真程序，分析巖土參數(shù)對環(huán)境變形的影響模式以及分析參數(shù)的離散性。具體到個案工程的參數(shù)反演中，反演的主要操作步驟如下，其流程如圖1所示。

(1)根據(jù)所反演的參數(shù)對環(huán)境變形的影響模式，確定反演分析的最少試算組。例如，某一參數(shù)以三次多項式規(guī)律影響環(huán)境變形，則反演分析中至少需要試算四組來標定影響關(guān)系式；

(2)根據(jù)按參考值的計算結(jié)果與監(jiān)測結(jié)果的比較，確定參數(shù)反演計算方案中參數(shù)的調(diào)整趨勢。例如，按參考值計算的變形小于監(jiān)測結(jié)果，則參數(shù)應該向減小環(huán)境變形的方向調(diào)整，使反演參數(shù)是通過內(nèi)插求得，而不是外延求得；

(3)根據(jù)反演參數(shù)的離散性分析結(jié)果，確定參數(shù)的合理反演區(qū)間；

(4)根據(jù)上述三項制定計算方案；

(5)按計算方案計算，并回歸建立針對本工程的反演參數(shù)與環(huán)境變形之間的定量關(guān)系式；

(6)根據(jù)監(jiān)測數(shù)據(jù)與計算結(jié)果之間的關(guān)系，建立反演優(yōu)化的目標函數(shù)；

(7)根據(jù)目標函數(shù)全局最優(yōu)原則得到反演參數(shù)的優(yōu)化值；

(8)判斷反演優(yōu)化值是否合理。如果反演優(yōu)化值落于反演區(qū)間之內(nèi)，表明結(jié)果合理，該優(yōu)化值就是反演結(jié)果；如果優(yōu)化值落于反演區(qū)間以外(如大于反演上界)，表明結(jié)果不合理，以反演區(qū)間邊界(如上界)作為反演結(jié)果；反演結(jié)束。

圖1　參數(shù)反演流程

2數(shù)值仿真技術(shù)

數(shù)值仿真方法作為本文參數(shù)反演的工具，主要有有限差分法、有限元法、離散元法等。目前巖土工程數(shù)值仿真分析日趨成熟，但大都是由專業(yè)人員進行專題分析，計算周期較長，普及率較低。本文反演方法中要求在短時間內(nèi)進行多次數(shù)值計算，現(xiàn)有計算軟件無法滿足。筆者以EXCEL數(shù)據(jù)庫和FLAC3D[5]及其fish語言作為開發(fā)平臺，探索開發(fā)了土釘支護邊坡數(shù)值仿真模塊建模程序，圖2是參數(shù)輸入界面，圖3是生成的計算模型。該方法中以EXCEL數(shù)據(jù)表和工程性語言的形式輸入?yún)?shù)，其會自動轉(zhuǎn)換并生成FLAC3D模型。該程序適用于土釘?shù)罃?shù)≤20，土層數(shù)≤10，坡度等其他參數(shù)不限。該方法解決了建模問題，一次建模時間不超過10 min，建模效率較高。但計算核心為FLAC3D程序，計算時間較長，數(shù)值仿真普遍存在的計算耗時的問題還需要更多的學者共同努力解決。筆者認為，隨著計算技術(shù)的不斷發(fā)展，尤其是計算機硬件的不斷提高，數(shù)值計算時間長的問題能夠得到有效解決。

圖2　參數(shù)輸入界面

圖3　土釘墻模型

3巖土參數(shù)對環(huán)境變形的影響模式

巖土參數(shù)對環(huán)境變形的影響模式是本文參數(shù)反演方法的理論基礎(chǔ)之一，這里采用實例的方式來闡述其實現(xiàn)過程。

利用本文第2節(jié)提及的土釘支護基坑模塊建模程序來分析土體變形模量對土釘墻頂水平位移和垂直位移的影響模式。選取文獻[6]中介紹的麗水嘉園基坑工程作為分析實例，文獻中提供的參數(shù)值(勘察數(shù)據(jù))作為基準狀態(tài)參數(shù)值(表1)。在基準參數(shù)值的基礎(chǔ)上按比例調(diào)整參數(shù)，通過數(shù)值計算得到基坑位移的變化結(jié)果。從參數(shù)值的實際可能變化范圍考慮，參數(shù)在參考值的基礎(chǔ)上上下浮動20%和40%，計算得到變形模量變化時的5組環(huán)境變形，調(diào)整后的參數(shù)及其計算結(jié)果如表2所示，環(huán)境位移隨變形模量調(diào)整系數(shù)變化的曲線如圖4所示。結(jié)果表明，變形模量與土釘墻位移之間的關(guān)系基本符合三次多項式關(guān)系模式。

其中，三次多項式1為：

y=92.292x3-356.39x2+485.54x-260.86

三次多項式2為：

y=99.063x3-393.28x2+550.58x-306.85

表1　基準狀態(tài)參數(shù)值

表2　坡頂位移計算結(jié)果

圖4　變形模量-位移關(guān)系曲線

4巖土測試數(shù)據(jù)的離散性分析

4.1分析方法

巖土測試數(shù)據(jù)的離散性是本文反演方法理論基礎(chǔ)之二。工程中最主要的巖土測試數(shù)據(jù)是地質(zhì)勘察數(shù)據(jù)，其為設(shè)計的最主要依據(jù)，但是巖土勘察數(shù)據(jù)是不夠準確的，原因主要有：(1)土體在空間上不均勻，而鉆孔空間離散且數(shù)量有限，因此鉆孔資料不能準確描述全部的地層參數(shù)。(2)土體是由顆粒組成的彈塑性材料，其物理力學性能與應力歷史、應力狀態(tài)和應力路徑有關(guān)，土工試驗會對土體產(chǎn)生擾動，導致應力環(huán)境發(fā)生改變。(3)土工試驗儀器不夠精密，人員操作也往往會帶來誤差。

從上面分析可知，巖土參數(shù)勘查數(shù)據(jù)的不準確性由多個因素造成，對每個影響因素都深入研究清楚固然最理想，但很難實現(xiàn)。退而求其次，這里從宏觀上進行統(tǒng)計分析。在地勘報告中，一項指標通常包括測試樣本數(shù)、最大值、最小值、平均值等。這些數(shù)值都體現(xiàn)了測試結(jié)果的離散分布特性，也體現(xiàn)了巖土參數(shù)真實值可能出現(xiàn)的范圍和概率。對于大樣本，統(tǒng)計結(jié)果能夠確定參數(shù)真實值所在的區(qū)間。

為了對勘察數(shù)據(jù)的離散性進行量化分析，這里假設(shè)某一參數(shù)a，勘察試驗樣本數(shù)為N，最大測試結(jié)果為amax，最小測試結(jié)果為amin，測試結(jié)果平均值為aav。定義：相對最大值ρmax=最大值amax/平均值aav≥1，相對最小值ρmin=最小值amin/平均值aav≤1，相對平均值ρmin=平均值aav/平均值aav=1。例如對于壓縮模量，計算每組數(shù)據(jù)的ρmax，統(tǒng)計ρmax在區(qū)間[1+Δn，1+Δ(n+1) )內(nèi)的出現(xiàn)的頻次(Δ為統(tǒng)計步長，n=0,1,2,3…)。按區(qū)間從小到大順序?qū)Τ霈F(xiàn)的頻次進行累加，得到累積出現(xiàn)頻次。用累積出現(xiàn)頻次除以統(tǒng)計樣本數(shù)，就得到累積概率密度。定義累積概率密度的超越概率為5%時對應的區(qū)間點值為離散上界ρmax95。同理，可以對巖土參數(shù)項進行離散下界分析。定義離散區(qū)間為[ρmin95,ρmax95]。本文中暫取離散界限對應的超越概率為5%，該值應該根據(jù)更深入的研究確定。

4.2巖土勘察數(shù)據(jù)的離散性分析的分析實例

本文收集了北京、深圳、天津、成都、青島、合肥6市19個工程87層土的巖土勘察數(shù)據(jù)，通過第4.1節(jié)中的分析方法對結(jié)果進行統(tǒng)計(表3)。

表3　壓縮模量統(tǒng)計結(jié)果

離散下界通過離散小值的統(tǒng)計得到，統(tǒng)計區(qū)間取[1-0.02n，1-0.02(n+1) )(n=0,1,2,3…，Δ=0.02)，結(jié)果如表3所示。同理，離散上界通過離散大值統(tǒng)計得到，取統(tǒng)計區(qū)間為[1+0.1n，1+0.1(n+1) )(n=0,1,2,3…，Δ=0.1)。鑒于篇幅，具體數(shù)據(jù)不再贅述。

根據(jù)上述統(tǒng)計得到的累積概率密度在統(tǒng)計區(qū)間內(nèi)的分布，繪制累積概率密度分布曲線如圖5所示。從圖中可知，超越概率為5%的離散區(qū)間為[0.564，1.641]，即表明壓縮模量真實值與勘察值的比為0.564～1.641的可能性為95%。所以在壓縮模量反演中，可認為當反演值與勘察值的比值落于[0.564，1.641]之內(nèi)時是合理的，落于之外時是不合理的，取相應的邊界值。

圖5　累積概率密度曲線

由于筆者收集到的巖土勘察數(shù)據(jù)的局限性，這里的分析結(jié)果是不夠理想的。筆者認為，巖土測試數(shù)據(jù)的離散性分析應該分地域、土性、測試方法進行分類統(tǒng)計，這有待于后續(xù)深入研究。

5應用實例

5.1實例概況

取文獻[6]中介紹的麗水嘉園土釘支護基坑工程作為反演算例。該工程在文獻中給出的巖土參數(shù)作為參考值，利用參考值計算坡頂水平位移和垂直位移，在該位移值的基礎(chǔ)上增加5%作為目標位移來模擬實測位移值，目標位移值如表4所示。

表4　目標位移值

5.2反演過程

根據(jù)本文第3節(jié)研究結(jié)果，變形模量變化與土釘支護邊坡的位移基本符合三次多項式關(guān)系模式，所以本部分反演計算中取4組數(shù)據(jù)試算。根據(jù)本文第4.2節(jié)中對壓縮模量(壓縮模量與變形模量規(guī)律一致)勘察結(jié)果的統(tǒng)計分析，離散區(qū)間為[0.564，1.641]，同時考慮目標位移值比參考位移值大，即變形模量反演值應該比參考值小，所以調(diào)整參數(shù)取0.6、0.8、1.0和1.2進行試算，計算結(jié)果如表5所示。

表5　位移試算結(jié)果

對表5中數(shù)據(jù)進行三次曲線擬合，得到變形模量調(diào)整系數(shù)和坡頂水平位移的關(guān)系式為：

Wh=18.75t3-179.88t2+367.93t-256.62

變形模量調(diào)整系數(shù)和坡頂垂直位移的關(guān)系式為：

Wh=24.16t3-175.38t2+330.61t-218.26

建立反演優(yōu)化目標函數(shù)為：

F[Eopt]=min{(wh(E)-wh)2+(wv(E)-wv)2}

根據(jù)上述擬合關(guān)系式與上述目標函數(shù)，計算得到當水平位移差和垂直位移差平方和最小時，變形模量調(diào)整系數(shù)為0.965。

5.3反演結(jié)果合理性判斷及結(jié)果驗證

根據(jù)本文第4.2節(jié)統(tǒng)計結(jié)果，壓縮模量超越概率為5%的離散區(qū)間為[0.564，1.641]，本節(jié)參數(shù)反演調(diào)整系數(shù)為0.965，處于離散區(qū)間之內(nèi)，因此認為是合理的。

將參考變形模量參數(shù)值乘以0.965便是反演得到的變形模量值，結(jié)果如表6所示。

表6　變形模量反演值

對上述得到的反演值進行驗證計算，結(jié)果如表7所示。結(jié)果表明，計算值與目標值之間相差較小。

表7　目標值與反演值比較

6結(jié)論

本文提出的基于數(shù)值仿真和巖土特性的參數(shù)反演方法為工程反分析提供了一種新的思路，具有以下優(yōu)點：

(1)以數(shù)值仿真作為工具，以巖土特性作為基礎(chǔ)，該方法具有較強的理論依據(jù)和可操作性。

(2)本文中的參數(shù)對環(huán)境變形的影響，在前期研究階段確定模式，在個案分析中標定關(guān)系式，因此得到的關(guān)系式更具針對性，能夠更好的考慮巖土參數(shù)的普遍性和具體工程的特殊性。

(3)以勘察數(shù)據(jù)的離散性來判斷反演值是否合理，簡單易行，意義明確。

本文反演方法具有很好的發(fā)展前景，但其中的數(shù)值仿真技術(shù)、巖土參數(shù)對環(huán)境變形的影響規(guī)律及巖土測試數(shù)據(jù)的離散性需要進行更加深入的研究。

參考文獻

[1]Sakurai S , Abe S. A design approach to dimensioningunde-rground openings [A] . In : Proc 3rd Int Conf Numerical Methods in Geomechanics[C] .Aachen ,1979： 649-661.

[2]Sakurai ,S ,Takeuchi K.Back analysis of measured displacementof tunnel [J] . Rock Mech and Rock Eng ,1983 ,16 (3) :173-180.

[3]Gioda G, Pandolfi A ,Cividini A. A comparative evaluation ofsome back analysis algorithms and their application to in2situload tests[A]. In :Proc 2nd Int Symp on Field Measurement inGeom[C] . Kobe ,1987： 1131-1144.

[4]楊志法. 有限元法圖譜[M] . 北京:科學出版社,1988.

[5]ITASCA 軟件公司. FLAC程序用戶使用手冊[Z].

[6]段啟偉.土釘支護現(xiàn)場實測與數(shù)值模擬[D]. 北京交通大學,2007.

[作者簡介]王琛(1978～)，男，本科，工程師，主要從事地鐵等地下結(jié)構(gòu)的設(shè)計和研究工作。

【中圖分類號】TU441+.3

【文獻標志碼】A

[定稿日期]2015-10-12