王日冉，馬清文,2

古墓垂直防滲帷幕安全深度的優(yōu)化設(shè)計(jì)

王日冉1，馬清文1,2

(1.鄭州大學(xué)水利與環(huán)境學(xué)院，河南 鄭州 450001；2.水利與交通基礎(chǔ)設(shè)施安全防護(hù)

河南省協(xié)同創(chuàng)新中心，河南 鄭州 450001)

古墓作為特殊的文物保護(hù)群體，在實(shí)施保護(hù)時(shí)尤其要注意原狀維護(hù)和病害防治。地下水滲漏問題已經(jīng)發(fā)展成為古墓最主要的病害之一，目前的治理方法主要是在古墓周邊布置垂直防滲帷幕，形成不透水地下連續(xù)墻。對(duì)防滲帷幕的安全深度進(jìn)行優(yōu)化設(shè)計(jì)成為解決影響防滲效果和工程造價(jià)問題的關(guān)鍵。結(jié)合地下水滲流理論，分析防滲帷幕深度的主要影響因素，以墓室底部土體不發(fā)生滲透破壞為控制條件計(jì)算垂直防滲帷幕插入深度，并通過工程實(shí)例利用ANSYS軟件模擬地下水滲流場(chǎng)，對(duì)理論計(jì)算結(jié)果進(jìn)行驗(yàn)證。結(jié)果表明：理論計(jì)算得出的古墓垂直防滲帷幕安全深度是可靠的，且滿足劉崇墓古墓的防滲保護(hù)要求。

古墓；垂直防滲帷幕；安全深度；滲流場(chǎng)

大型墓葬遺址的真實(shí)性和不可再生性決定了保護(hù)工作的迫切性與科學(xué)性。經(jīng)過走訪、調(diào)研及查閱歷史文獻(xiàn)資料發(fā)現(xiàn)，在已經(jīng)發(fā)掘的國(guó)家級(jí)及省級(jí)古墓葬中存在大量地表、地下水滲漏和結(jié)構(gòu)破壞問題，如果任其發(fā)展，長(zhǎng)此下去便會(huì)給文物帶來極大的安全隱患。由于墓葬大都是地下洞室，古墓的滲水來源主要以墓地周邊地表及地下水滲流為主，雨后滲水更多。該問題類似于工程上常見的基坑滲水問題，不同之處就是在治理時(shí)要遵循不改變文物原狀和減少人為干預(yù)的原則。目前主要采用疏、堵、排的治理方法，在墓室周邊建立封閉的垂直防滲帷幕，形成不透水地下連續(xù)墻，用來減少滲流量和降低滲流速度，以達(dá)到保護(hù)墓室的目的。

目前對(duì)于地下防滲墻及其覆蓋層的滲透穩(wěn)定問題已有一些研究成果。如毛昶熙等[1]在實(shí)驗(yàn)室水槽中進(jìn)行了懸掛式防滲墻的砂質(zhì)堤基的滲流模型管涌試驗(yàn)，確定了懸掛式防滲墻及其前后水平段的滲流臨界坡降，并進(jìn)行了相關(guān)的理論分析;沈振中等[2]提出了懸掛式混凝土防滲墻深度優(yōu)化設(shè)計(jì)方案;張家發(fā)等[3]通過砂槽模型模擬了滲透變形發(fā)生和擴(kuò)展的過程，指出懸掛式防滲墻對(duì)滲透變形的發(fā)生條件影響很小;謝興華等[4]對(duì)深厚覆蓋層壩基防滲墻最優(yōu)深度進(jìn)行了研究,并指出最優(yōu)深度的取值取決于覆蓋層的深度;楊波等[5]考慮防滲帷幕強(qiáng)度和彈性模量值對(duì)超深防滲帷幕的影響而進(jìn)行優(yōu)化設(shè)計(jì)，同時(shí)以基坑降水深度和允許滲透坡降來確定防滲帷幕厚度;宋永占等[6]采用了Fluent數(shù)值模擬程序?qū)Ψ罎B帷幕在不同貫入深度的工況下進(jìn)行分析,從而得到了防滲帷幕深度合理值；李建華等[7]根據(jù)邊界元的特點(diǎn)和計(jì)算理論，分析了防滲墻位置和深度對(duì)滲流量的影響；楊秀竹等[8]結(jié)合工程實(shí)例說明懸掛式防滲帷幕在建造前后防滲效果顯著，并指出懸掛式防滲帷幕底部小范圍內(nèi)需要提高抗沖刷能力；徐書平等[9]實(shí)現(xiàn)了滲流分析邊界條件處理的新算法；孫慶春等[10]提出的修復(fù)地下水重金屬污染的滲透反應(yīng)墻技術(shù)對(duì)土遺址地下水環(huán)境保護(hù)具有借鑒意義；陳義等[11]應(yīng)用ANSYS對(duì)裂隙全充填的巖體進(jìn)行了數(shù)值模擬計(jì)算，該研究對(duì)于古墓地下防滲墻體的數(shù)值模擬具有參考價(jià)值?？傊?，防滲問題一直以來是土木工程和水利工程領(lǐng)域內(nèi)的重大問題，且防滲技術(shù)與許多工程實(shí)踐密切相關(guān)。由于滲流而產(chǎn)生的工程問題主要有兩大類：滲透變形問題和滲流量與滲流速率的問題，本文結(jié)合這兩方面問題和上述在防滲工程中對(duì)防滲帷幕設(shè)計(jì)方面的研究，主要考慮墓區(qū)土層參數(shù)和地下水位特征等因素，對(duì)古墓防滲帷幕安全深度進(jìn)行了優(yōu)化設(shè)計(jì)。

1 地下水滲流理論及防滲帷幕安全深度影響因素

1.1 二維穩(wěn)定滲流的基本方程

在研究各向同性土體的二維穩(wěn)定滲流基本方程的建立方法時(shí)，假定從穩(wěn)定滲流場(chǎng)中任取一微小的土單元體，在滲流作用下單元體的體積保持不變且水不可壓縮，單位時(shí)間內(nèi)流入單元體的總水量等于流出的總水量，根據(jù)達(dá)西定律和連續(xù)性方程有：

(1)

式中:k為土體滲透系數(shù)(m/s)；h為總水頭或測(cè)壓管水頭(m)。

式(1)化簡(jiǎn)后為

(2)

式(2)即為各向同性土體二維穩(wěn)定滲流的基本方程式，也就是著名的拉普拉斯(Laplace)方程。滲流速度的分量為

1.2 有限元方程

二維穩(wěn)定滲流屬于帶直邊界條件的問題，可以采用線性形函數(shù)近似獨(dú)立變量的方法來解決地下水滲流水頭的空間變化問題?，F(xiàn)對(duì)二維穩(wěn)定滲流問題采用加重余差的迦遼金(Galerkin)有限元方法分析[12]，表征的基本方程為

(∫V[B]T[D][B]dA){H}=q∫A[N]TdL

(3)

式中：[B]為單元坡降矩陣；[D]為單元滲透系數(shù)矩陣；[N]為單元形函數(shù)矩陣；{H}為水頭列陣；q為單元邊界單寬流量；dL和dA分別為長(zhǎng)度積分變量和面積積分變量；V為滲流區(qū)域體積；A為滲流過流面積。

將式(3)進(jìn)一步化簡(jiǎn)，得

[K]{H}={Q}

(4)

式中:[K]表示單元滲透矩陣,[K]=∫V[B]T[D][B]dA;{Q}表示單元邊界流量列陣,{Q}=q∫A[N]TdL。

分析地下水滲流理論，探究影響防滲帷幕安全深度的主要因素，是古墓防滲帷幕優(yōu)化設(shè)計(jì)的關(guān)鍵。

1.3 影響防滲帷幕安全深度的主要因素

1.3.1 土層信息

多數(shù)古墓地下埋深較大，墓冢及墓室邊墻周邊填土土質(zhì)松軟，且雜草、灌木恣意生長(zhǎng)，根系縱橫交錯(cuò)，造成土層滲透性的增加，另外土質(zhì)墓室的墓道通常以夯土為地基，夯土因自然力作用分化，加劇墓室滲水。因此，土層信息對(duì)于古墓地下水防滲至關(guān)重要。古墓防滲帷幕的防滲效果與帷幕的深度、厚度、位置因素有關(guān)，而關(guān)鍵因素是防滲帷幕的安全深度。帷幕的防滲作用對(duì)土層的滲透系數(shù)和帷幕本身的滲透性很敏感，當(dāng)古墓防滲帷幕采用不透水材料時(shí)，帷幕本身的滲透性可以忽略，所以土層信息成為影響防滲帷幕安全深度的關(guān)鍵因素之一，原因是土層信息直接影響滲透系數(shù)的大小，如土的粒度成分、密實(shí)度、飽和度、結(jié)構(gòu)和構(gòu)造等信息。實(shí)際工程中的天然沉積土往往由滲透性不同的土層組成，對(duì)于與土層平行和垂直的滲流情況，當(dāng)各個(gè)土層的滲透性和厚度為已知時(shí)，據(jù)此可以衡量整個(gè)土層的滲透性。

1.3.2 地下水位特征

多數(shù)古墓在發(fā)掘之前長(zhǎng)期處在水量豐富的地下環(huán)境中，墓墻已嚴(yán)重老化，多處存在酥堿、剝蝕破壞，墓底地磚酥軟破碎，已基本喪失防水功能。因此，地下水位特征是影響防滲帷幕安全深度的另一個(gè)重要因素，原因是地下水位直接決定水頭的大小，而水頭差是滲流發(fā)生的必要條件。

1.3.3 不透水層埋深

根據(jù)古墓地基不透水層埋深情況，并根據(jù)實(shí)際需要，垂直防滲帷幕可完全切斷地基的透水層，深入到不透水層中，徹底解決地基土的滲透變形問題，也可不完全切斷透水層，做成懸掛式防滲帷幕，起到延長(zhǎng)滲流途徑、降低下游的逸出坡降的作用；當(dāng)相對(duì)不透水層埋深大且無規(guī)律時(shí)，應(yīng)根據(jù)滲透分析及防滲要求并結(jié)合工程經(jīng)驗(yàn)決定防滲帷幕的安全深度。

2 工程實(shí)例分析

2.1 工程概況

河南省周口市淮陽縣劉崇墓遺址是一座東漢中期的大型磚石墓葬，古墓東西長(zhǎng)28.3 m，南北寬18.2 m，墓室底深5.26 m，總面積為515 m2，2013年5月公布為第七批全國(guó)重點(diǎn)文物保護(hù)單位。由于地下水位上升，古墓墓室內(nèi)產(chǎn)生嚴(yán)重滲水病害，長(zhǎng)期的滲水導(dǎo)致墓頂、墓壁遭到嚴(yán)重的溶蝕，墓道積水導(dǎo)致墓磚嚴(yán)重破壞等問題。為截?cái)酀B流途徑，降低滲流速度，減少滲流量，在古墓東側(cè)1.5 m、西側(cè)6 m、南側(cè)6 m、北側(cè)6 m處設(shè)置一道地下連續(xù)防滲墻，墻厚20 cm，采用超薄型高聚物防滲材料，整體形成“口”字形防滲墻。經(jīng)勘察，近5年最高水位相對(duì)高程為-3.00 m左右，而該古墓的基礎(chǔ)埋深為5.26 m，地下水類型基本為孔隙型潛水。依據(jù)地勘資料，土層自上而下分為：①素填土、②粉土、③黏土、④粉質(zhì)黏土、⑤黏土。土層參數(shù)信息見表1。

表1 土層參數(shù)信息表

2.2 防滲帷幕插入深度分析

垂直防滲帷幕插人基坑開挖高程以下深度是影響防滲效果和工程造價(jià)的一個(gè)重要因素。經(jīng)調(diào)研發(fā)現(xiàn)，劉崇墓墓室內(nèi)地面有局部隆起現(xiàn)象，地面磚塊被頂起，說明流土現(xiàn)象存在,因此在滿足墓室結(jié)構(gòu)穩(wěn)定的前提下，古墓周圍防滲帷幕安全插入深度的確定是以使得墓室底部土體不發(fā)生流土為控制條件來進(jìn)行分析的。此外，防滲帷幕的最終安全深度要經(jīng)過試算確定，鑒于粉質(zhì)黏土層等弱透水層厚度較大，可以設(shè)置懸掛式豎向防滲帷幕來阻隔地下水。如圖1所示，當(dāng)墓室外的水繞過防滲帷幕向墓室底部滲透時(shí)，只有當(dāng)滲透坡降i≤[i]([i]為墓室底部土體的允許滲透坡降)時(shí)，墓室底部土體才不會(huì)發(fā)生滲透破壞。

圖1 防滲帷幕插入深度分析剖面圖Fig.1 Cross-section of the insert depth analysis of anti-seepage curtain

計(jì)算土體滲透坡降時(shí)，應(yīng)選擇緊貼帷幕的最短滲徑計(jì)算，此處的滲透坡降將最大，只要緊貼防滲帷幕內(nèi)側(cè)的土體不發(fā)生流土，便可認(rèn)為墓室底部土體不會(huì)遭到滲透破壞。其計(jì)算公式如下：

(5)

式中:i為水力坡降(m)；d為墓室底部以下防滲帷幕深度(m)；d+5.26為防滲帷幕的總插人深度(m)；b為防滲帷幕的寬度(m),b=0.2 m；h為地下水位至墓室地面的距離(m)，h=5.26-3=2.26 m;[i]為允許滲透坡降(m)，結(jié)合《工程地質(zhì)手冊(cè)》和工程實(shí)踐經(jīng)驗(yàn)[i]取0.3。

因此，防滲帷幕的插入深度為

(6)

將各數(shù)據(jù)代入(6)式計(jì)算得d?2.53 m,考慮古墓保護(hù)工程的重要性和特殊性，計(jì)算防滲帷幕的總插入深度時(shí)取安全系數(shù)為1.5，即防滲帷幕的總插入深度為(2.53+5.08)×1.5=11.42 m，綜合考慮取為12 m，即深入到粉質(zhì)黏土弱透水層，形成懸掛式防滲帷幕。

2.3 ANSYS滲流分析

ANSYS軟件具有強(qiáng)大的分析功能，能夠模擬分析自然界中的應(yīng)力場(chǎng)、電磁場(chǎng)、溫度場(chǎng)和流場(chǎng)，并且支持多場(chǎng)耦合分析。ANSYS軟件提供了大量的單元用于流體力學(xué)問題的模型分析，這些單元包括FLUID15、FLUID66、FLUID79，由于各控制方程的相似性，除了上述幾種單元外，還可以用熱實(shí)體單元中6節(jié)點(diǎn)的三角形單元PLAIN35、4節(jié)點(diǎn)的四邊形單元PLAIN55和8節(jié)點(diǎn)的四邊形單元PLAIN77，本文采用8節(jié)點(diǎn)Solid單元PLAIN77來模擬分析地下滲流問題[13]。滲流分析時(shí)有以下三點(diǎn)主要假設(shè)：將滲流問題簡(jiǎn)化為平面問題；防滲墻采用不透水的防滲材料；地下水位作為滲透的唯一因素。

對(duì)于典型的穩(wěn)定滲流問題，將已知的固定水位邊界作為水頭邊界，另外還有出滲邊界和防滲墻邊界，在求解時(shí)沒有給出邊界條件的區(qū)域邊界默認(rèn)是不透水的，對(duì)應(yīng)于滲流控制方程中的自然邊界條件。液面水頭通常是已知的，分析時(shí)將已知水頭作為確定的邊界條件。假定滲流是發(fā)生在各向同性土體內(nèi)的二維穩(wěn)定流,將地下水位邊界和墓室底部滲流滲出處作為第一類邊界條件；將墓室底部相對(duì)不透水的黏土層和防滲帷幕及四周不透水的周圍介質(zhì)作為第二類邊界條件。在仿真分析時(shí)，選取靠近墓室邊墻一側(cè)的防滲帷幕進(jìn)行模擬，具體邊界尺寸如圖2所示。

圖2 地下水滲流場(chǎng)邊界條件Fig.2 Boundary conditions of the groundwater seepage field

本文利用ANSYS軟件，在防滲帷幕深度擬定的情況下，對(duì)地下水滲流場(chǎng)進(jìn)行分析?？紤]到滲流場(chǎng)被簡(jiǎn)化成了規(guī)則圖形，在網(wǎng)格劃分時(shí)采用自由劃分，設(shè)置單元長(zhǎng)度為0.4 m，在應(yīng)力可能集中的部位進(jìn)行了細(xì)化，設(shè)置單元長(zhǎng)度為0.2 m，共劃分1 550個(gè)單元，網(wǎng)格剖分結(jié)果見圖3。土壤中滲流速度矢量分布見圖4，沿路徑A-B上滲流速度的變化見圖5。

圖3 網(wǎng)格剖分圖Fig.3 Grid subdivision graph

圖4 土壤中滲流速度的矢量分布Fig.4 Distribution of the vectors of seepage velocity in the soil

圖5 沿路徑A-B上滲流速度的變化Fig.5 Change of the seepage velocity along the path of A-B

由圖4可見，最大滲流流速分布在防滲帷幕底部，其值為0.002 788 m/s，最小滲流流速為0.999×10-5m/s，其值非常小，即使在最不利情況下，一定時(shí)間內(nèi)的滲透量也是非常有限的。另外，由圖5可見，滲流流速在A點(diǎn)最大，其值為4.24×10-4m/s，并沿路徑A-B下降，在B點(diǎn)滲流流速降為零，最后達(dá)到自由流體狀態(tài)。依據(jù)滲流理論，A點(diǎn)附近的滲流流速要比B點(diǎn)滲流流速大，原因是A點(diǎn)對(duì)水流有較低的阻力，在A點(diǎn)附近有更多的水流經(jīng)過。因此，此次有限元分析結(jié)果驗(yàn)證了理論解的可靠性，從而說明了設(shè)置12 m安全深度的防滲帷幕可滿足古墓防滲的要求。

3 結(jié) 論

古墓保護(hù)工作的特殊性決定了防滲工作的特殊性，因此防滲設(shè)計(jì)必須科學(xué)嚴(yán)謹(jǐn)，工程處理更要做到周全到位，只有這樣才能真正達(dá)到保護(hù)古墓的目的。本文在介紹各向同性土體二維穩(wěn)定地下水滲流基本規(guī)律的基礎(chǔ)上，采用加重余差的迦遼金法的有限元方程進(jìn)一步分析了影響防滲帷幕安全深度的主要因素，并以墓室底部土體的允許滲透坡降為控制因素計(jì)算出防滲帷幕的插入深度，最后通過工程實(shí)例利用ANSYS軟件模擬的滲流場(chǎng)對(duì)理論計(jì)算結(jié)果進(jìn)行了驗(yàn)證。結(jié)果表明：理論計(jì)算出的防滲帷幕安全深度是可靠的，且滿足劉崇墓的防滲保護(hù)要求。因此，本文提出的方法對(duì)形成專門針對(duì)古墓等特殊土遺址群完整而科學(xué)的防滲治理措施具有一定的參考價(jià)值。

[1] 毛昶熙，段祥寶，蔡金傍，等.懸掛式防滲墻控制管涌發(fā)展的試驗(yàn)研究[J].水利學(xué)報(bào),2005,36(1):42-50.

[2] 沈振中,張?chǎng)?陸希,等.西藏老虎嘴水電站左岸滲流控制優(yōu)化[J].水利學(xué)報(bào)，2006，37(10)：1230-1234.

[3] 張家發(fā)，吳昌瑜，朱國(guó)勝.堤基滲透變形擴(kuò)展過程及懸掛式防滲墻控制作用的試驗(yàn)?zāi)M[J].水利學(xué)報(bào)，2002(9)：108-111，116.

[4] 謝興華，王國(guó)慶.深厚覆蓋層壩基防滲墻深度研究[J].巖土力學(xué)，2009,30(9)：2708-2712.

[5] 楊波，劉國(guó)彬，楊光煦.超深防滲帷幕性質(zhì)研究及設(shè)計(jì)優(yōu)化[J].巖土工程技術(shù),2004,18(5):267-270.

[6] 宋永占，趙廷紅，丁睿，等.混凝土重力壩壩基防滲帷幕作用的單因素分析[J].蘭州理工大學(xué)學(xué)報(bào)，2011,37(6)：120-125.

[7] 李建華，侍克斌，王建祥.土石壩防滲墻位置和深度對(duì)滲流量的影響[J].河北工程大學(xué)學(xué)報(bào)，2009，26(2)：36-38.

[8] 楊秀竹，陳福全，雷金山，等.懸掛式帷幕防滲作用的有限元模擬[J].巖土力學(xué)，2005,26(1)：105-107.

[9] 徐書平，喻國(guó)安，楊劍飛.滲流分析邊界條件處理的算法實(shí)現(xiàn)[J].土工基礎(chǔ)，2002,16(2):51-53.

[10]孫慶春，崔康平，許為義,等.重金屬污染地下水修復(fù)的滲透反應(yīng)墻技術(shù)[J].安全與環(huán)境工程，2013,20(5):53-56.

[11]陳義，王傳琪，陳金剛.基于ANSYS的裂隙巖體全充填介質(zhì)膨脹效應(yīng)的數(shù)值模擬[J].安全與環(huán)境工程，2012,19(4)：115-118，128.

[12]潘榮升.高壓噴射注漿防滲帷幕防滲機(jī)理及其應(yīng)用研究[D].長(zhǎng)沙：中南大學(xué)，2012.

[13]Saeed Moaveni.有限元分析—ANSYS理論與應(yīng)用[M].王崧等譯.2版.北京：電子工業(yè)出版社，2005.

Optimal Design about the Safety Depth of Vertical Anti-seepage Curtain in Ancient Tombs

WANG Riran1，MA Qingwen1,2

(1.SchoolofWaterConservancyandEnvironment,ZhengzhouUniversity，Zhengzhou450001，China；2.CollaborativeInnovationCenterofWaterConservancyandTransportationInfrastructureSafetyProtectionofHenanProvince，Zhengzhou450001，China)

Ancient tombs are a special group of cultural relics,of which the original state maintenance and disease control should be paid more attention to.Groundwater seepage problem has become one of the main diseases,and the treatment at present is mainly arrangement of vertical anti-seepage curtain aimed at forming impermeable underground continuous wall.To obtain the optimal design of the safety depth of the anti-seepage curtain becomes the key to solving the problem of seepage control effect and engineering cost.Combined with the groundwater seepage theory,this paper analyzes the main influencing factors of anti-seepage curtain depth,and calculates the safety depth of the vertical anti-seepage curtain by taking the nonoccurrence of the seepage failure as a control condition,and also testifies the theoretical calculation results by adopting the ANSYS software to simulate the groundwater seepage field through a project case study.The results indicate that the theoretical calculation of the safety depth of vertical anti-seepage curtain is reliable,which satisfies the requirements of the anti-seepage prevention of ancient Liuchong tombs.

ancient tomb;vertical anti-seepage curtain;safety depth;seepage field

1671-1556(2016)02-0125-05

2015-07-07

2015-11-20

河南省科技攻關(guān)項(xiàng)目(142102210071)

王日冉(1990—)，女，碩士研究生，主要研究方向?yàn)樗c交通基礎(chǔ)設(shè)施安全防護(hù)等。E-mail：wangriran23409708@163.com

X93;TU761.1+1

A

10.13578/j.cnki.issn.1671-1556.2016.02.023