陳冬琴

（長(zhǎng)江職業(yè)學(xué)院，湖北武漢 430074）

海綿城市概念提出至今，研究主要集中在“低影響開(kāi)發(fā)”、實(shí)現(xiàn)雨水控制目標(biāo)，如滲、滯、凈、蓄、用、排等多個(gè)方面。很多城市在海綿城市建設(shè)中，缺少對(duì)地質(zhì)環(huán)境研究，對(duì)場(chǎng)地環(huán)境不了解，海綿城市技術(shù)適用性不明確，導(dǎo)致出現(xiàn)很多地質(zhì)環(huán)境問(wèn)題，如地面沉降、開(kāi)裂、特殊土層膨脹等。海綿城市建設(shè)應(yīng)從自然水文條件角度出發(fā)，開(kāi)展城市地質(zhì)調(diào)查工作，分析海綿城市功能與環(huán)境地質(zhì)適宜性的基礎(chǔ)上進(jìn)行規(guī)劃設(shè)計(jì)，實(shí)現(xiàn)自然積存、自然滲透、自然凈化的功能。本文以某城市規(guī)劃區(qū)為研究對(duì)象，采用GMS軟件（表水及地下水耦合）計(jì)算地下水流場(chǎng)結(jié)合Terzaghi固結(jié)理論，進(jìn)行水動(dòng)力學(xué)-力學(xué)耦合數(shù)值模擬，模擬計(jì)算海綿城市滲蓄水時(shí)，對(duì)地下水位的影響，預(yù)測(cè)海綿城市滲蓄水對(duì)地質(zhì)環(huán)境地面沉降量的影響。

1 透水鋪張層及地質(zhì)條件

本次研究透水鋪裝層采用透水混凝土、級(jí)配碎石層、級(jí)配礫石層，共同特點(diǎn)是孔隙率高，具有很好的滲蓄水能力。下墊層地質(zhì)條件為：上部為粉砂土、全新統(tǒng)中期沖積物。中部為細(xì)砂土，分布于西南、西北部，一般頂板埋深1~1.5m；下部粉質(zhì)黏土，頂板埋深1~2m，厚度15~20m；底部為黏土，具水平微層理。

2 數(shù)學(xué)模型

2.1 地下水流動(dòng)數(shù)學(xué)模型

根據(jù)Darcy定律和水均衡原理，忽略密度變化對(duì)地下水影響的三維流動(dòng)數(shù)學(xué)模型為：

式中：Kxx、Kyy、Kzz分別為沿主方向的滲透系數(shù)[L/T]；h(x,y,z,t)為水位[L]；w為單位體積單位時(shí)間內(nèi)注入或抽出的水量[1/t]；為時(shí)間[T]；μs為單位儲(chǔ)水系數(shù)[1/L]；h0(x,y,z)為流場(chǎng)初始水位[L]；B1為一類(lèi)邊界條件；B2為二類(lèi)邊界條件；q(x,y,z,t)單位面積已知流量（l/t）。

2.2 地面沉降與地下水耦合模型

儲(chǔ)水系數(shù)的定義是當(dāng)水位下降一個(gè)單位時(shí)，從空隙介質(zhì)中釋放的水量（體積），空隙介質(zhì)被壓縮，從空隙介質(zhì)中釋放的水量?jī)?chǔ)水系數(shù)表達(dá)式為：

式中：γw為水的重度[M/（L*T）2]；為儲(chǔ)水系數(shù)[1/L]；n為孔隙度[-]；βw為水的體積壓縮系數(shù)[LT2/M]；α為土體體積壓縮系數(shù)[LT2/M]，表示為：

根據(jù)有效應(yīng)力原理，變形只考慮垂向方向，可以把地面沉降和地下水水位耦合起來(lái)：

式中：Em為壓縮模量[M/LT2]；b為厚度[L]。

式中：Δε為應(yīng)變?cè)隽?；Δσ為?yīng)力增量[M/（L*S2）]；Ee和Ev為彈性壓縮模量和塑性壓縮模量[M/（LS2）]。

3 數(shù)值模型構(gòu)建

3.1 研究區(qū)地質(zhì)結(jié)構(gòu)概念模型的構(gòu)建

選擇某城區(qū)海綿城市進(jìn)行數(shù)值模擬如圖1所示，根據(jù)該區(qū)上層透水鋪張、下層水文地質(zhì)條件，對(duì)模擬區(qū)模型進(jìn)行概化，該模型分為5層，總厚度30m，透水鋪張層根據(jù)人工碎石粒徑不同分為3層，第一層透水鋪張層厚度0.3m；第二層透水鋪張層厚度0.25m；第三層透水鋪張層厚0.25m，如表1所示，共0.8m。下部地質(zhì)條件為：上層粉砂土，一般厚度為5~10m；中部細(xì)砂土，厚度5~10m；下部粉質(zhì)黏土，具水平微層理，5~10m；底部為黏土層。

圖1 研究區(qū)概念模型

表1 透水鋪裝層填料情況

3.2 模型參數(shù)取值

3.2.1 參數(shù)基礎(chǔ)值選取

模型網(wǎng)格剖分平面上的網(wǎng)格間距為100m×100m，根據(jù)地層情況確定垂直方向的間距。

透水鋪張層的水文地質(zhì)參數(shù)根據(jù)室內(nèi)實(shí)驗(yàn)確定，模擬區(qū)初始水位根據(jù)2016年3月各層監(jiān)測(cè)孔觀(guān)測(cè)水位確定，模擬邊界根據(jù)模擬區(qū)1 000m附近的監(jiān)測(cè)孔觀(guān)測(cè)水位給定。地面沉降模擬應(yīng)力期2016年10月到2020年11月，時(shí)間間距為1個(gè)月。根據(jù)分析計(jì)算或是工程地質(zhì)分析時(shí)鉆孔中的土樣測(cè)試數(shù)據(jù)確定各層土體的各類(lèi)初始參數(shù)。

3.2.2 模型識(shí)別

模型識(shí)別是通過(guò)反復(fù)調(diào)整參數(shù)，提高數(shù)值模型的精度，使模型的計(jì)算結(jié)果接近實(shí)際值。本次計(jì)算的主要參數(shù)是滲透系數(shù)、儲(chǔ)水系數(shù)、壓縮模量等。模型識(shí)別的初始流場(chǎng)選用2016年10月的統(tǒng)測(cè)水位，擬合檢驗(yàn)以2018年7月的流場(chǎng)、變形監(jiān)測(cè)數(shù)據(jù)為依據(jù)，調(diào)整參數(shù)得識(shí)別后的土體參數(shù)見(jiàn)表2。

表2 識(shí)別后的各層土體參數(shù)表

3.3.3 海綿城市滲蓄水風(fēng)險(xiǎn)數(shù)值模擬

此次用于海綿城市滲蓄水風(fēng)險(xiǎn)數(shù)值模擬的預(yù)測(cè)方法是基于Terzaghi固結(jié)理論和GMS軟件建立。模擬海綿城市不同降雨強(qiáng)度及降雨歷時(shí)、地下水主位補(bǔ)給量、地面沉降度。

根據(jù)模擬區(qū)的降雨特征，豐水期是7月份，基于7月份流場(chǎng)，通過(guò)改變降雨入滲補(bǔ)給量，模擬3種降雨模式（表3），計(jì)算土壓力的變化值，分析降雨對(duì)海綿城市滲蓄水風(fēng)險(xiǎn)的影響如表4所示。

表3 3種降雨預(yù)測(cè)模式

表4 各層土體沉降比重統(tǒng)計(jì)表 單位：%

模擬海綿城市滲水所產(chǎn)生的沉降值。透水鋪張層在飽水情況下，表3中不同的3種降水模式釋放水量到蓄水之前的潛水面，土體沉降的變化值如表5所示。模擬海綿城市滲蓄水交替3次，所產(chǎn)生的沉降值如表6所示。

表5 各層土體沉降比重統(tǒng)計(jì)表 單位：%

表6 各層土體沉降比重統(tǒng)計(jì)表 單位：%

4 模擬結(jié)論分析

在GMS軟件計(jì)算地下水流場(chǎng)的基礎(chǔ)上，采用水動(dòng)力-力學(xué)耦合模型，以海綿城市的滲蓄性對(duì)地下水動(dòng)態(tài)的影響，模擬計(jì)算了土體沉降值，具體結(jié)論如下。

1）在模擬降雨情況下，海綿體在蓄水狀態(tài)時(shí)土體呈現(xiàn)上浮狀態(tài)。降雨量小，降雨歷時(shí)長(zhǎng)時(shí)，下部地層的上部粉砂土變形值最大達(dá)到 30.5%；當(dāng)降雨歷時(shí)增大，強(qiáng)度增強(qiáng)時(shí)，中部細(xì)砂土變形值最大達(dá)到29.9%。

2）海綿城市滲水時(shí)，釋放水量回到蓄水之前的潛水面，土體呈壓縮狀態(tài)。下層土體的變化值都較大，中部細(xì)砂土變形值最大達(dá)到 32.9%。

3）海綿城市在歷經(jīng)滲蓄水交替3次，整個(gè)鋪張層及下層土體都發(fā)生了較大沉降，在滲蓄水頻繁交替時(shí)對(duì)土體的破壞力最大，其中沉降最大的是中部細(xì)砂層達(dá)29%。

以上研究表明，透水鋪張層和下墊層一定要考慮到滲蓄水的風(fēng)險(xiǎn)性，必須經(jīng)過(guò)預(yù)測(cè)與評(píng)價(jià)使沉降值控制在不產(chǎn)生破壞的范圍內(nèi)。下墊層的土體不符合海綿體滲蓄水功能性要求時(shí)，需采取措施加固。