葛菲媛， 劉景蘭， 李立偉

（天津市地質(zhì)研究和海洋地質(zhì)中心，天津300170）

本文在某典型圍填造陸區(qū)入駐企業(yè)污染源調(diào)查的基礎(chǔ)上，結(jié)合野外實(shí)際勘查結(jié)果，利用Hydrus-1D對(duì)特征污染物泄漏進(jìn)入包氣帶后的運(yùn)動(dòng)過(guò)程進(jìn)行模擬，識(shí)別特征污染物的運(yùn)移規(guī)律，為圍填造陸區(qū)的包氣帶和地下水保護(hù)提供科學(xué)依據(jù)。

1 研究場(chǎng)地概況

研究區(qū)屬于海河平原，地貌類型為海積低平原亞區(qū)，屬溫暖帶半濕潤(rùn)大陸性季風(fēng)氣候。 根據(jù)本次現(xiàn)場(chǎng)勘查及水文地質(zhì)試驗(yàn)，研究場(chǎng)地內(nèi)包氣帶厚度1.81～2.66 m，平均厚度為2.12 m，包氣帶巖性為人工沖填土，垂向平均滲透系數(shù)為0.15 m／ d，研究場(chǎng)地的包氣帶防污性能屬弱等級(jí)。

通過(guò)野外實(shí)際調(diào)查可知，該場(chǎng)地入駐企業(yè)的污水處理站為鋼筋混凝土結(jié)構(gòu)，各池體均為地上式，池底發(fā)生泄漏不易發(fā)現(xiàn)，泄漏的廢水在重力作用下會(huì)對(duì)土壤和地下水造成污染，因此本次模擬的污染源為污水處理站的廢水。 同時(shí)根據(jù)該企業(yè)廢水進(jìn)水水質(zhì)中重金屬鎳的最大濃度為25 mg ／ L，且鎳在自然界中一旦反應(yīng)為羥基鎳能產(chǎn)生很強(qiáng)的毒性，因此本次預(yù)測(cè)因子選擇鎳。

2 模型建立

本次預(yù)測(cè)選擇污染物以點(diǎn)源形式垂直進(jìn)入包氣帶的情形，利用Hydrus-1D 的水流及溶質(zhì)運(yùn)移兩大模塊進(jìn)行預(yù)測(cè)，預(yù)測(cè)模型為一維非連續(xù)點(diǎn)源非飽和溶質(zhì)垂向運(yùn)移模型。 模型設(shè)定時(shí)間單位為d，質(zhì)量單位為mg，長(zhǎng)度單位為cm（后文數(shù)學(xué)模型中各參數(shù)單位的設(shè)定均與此一致）。

2.1 水流模型的選擇及參數(shù)設(shè)定

水流模型選擇發(fā)展已相對(duì)成熟，本次模擬選擇目前應(yīng)用最為廣泛的VG 模型來(lái)進(jìn)行模擬計(jì)算，不考慮水流運(yùn)動(dòng)的滯后現(xiàn)象。

本次模擬污水處理站池體防滲層出現(xiàn)破損發(fā)生跑冒滴漏，污染物進(jìn)入包氣帶的情形，故水流上邊界條件選擇大氣邊界-可積水。 本次模擬不考慮地下水水位變化對(duì)水流及溶質(zhì)運(yùn)移的影響，故選擇自由排水邊界（Free Drainage）作為下邊界條件。

Hydrus-1D 水流模塊中的Soil Catalog 項(xiàng)包含12種典型土壤介質(zhì)及其土壤水分特征曲線相關(guān)參數(shù)，而該場(chǎng)地包氣帶主要巖性為人工沖填土，不包含其內(nèi)。 因此本次模擬基于土工試驗(yàn)成果使用Neural network prediction 來(lái)計(jì)算土壤水分特征曲線參數(shù)。土工試驗(yàn)成果成果見(jiàn)表1，計(jì)算出的土壤水分特征曲線參數(shù)見(jiàn)表2，其中Ks采用場(chǎng)地滲水試驗(yàn)實(shí)測(cè)值，15 cm／ d。

表1 土工試驗(yàn)成果表

表2 水流模型的參數(shù)

2.2 溶質(zhì)運(yùn)移模型的選擇及參數(shù)設(shè)定

本次模擬使用經(jīng)典對(duì)流-彌散方程，不考慮吸附和各種零級(jí)、一級(jí)及其他反應(yīng)。

根據(jù)污水處理池泄漏的實(shí)際情況，溶質(zhì)運(yùn)移上邊界選擇濃度通量邊界，下邊界選擇零濃度梯度邊界。

本次模擬假設(shè)池體發(fā)生泄漏后，企業(yè)在30 d 可以發(fā)生泄漏并切斷污染源，池體的泄漏量參考《給水排水構(gòu)筑物工程施工及驗(yàn)收規(guī)范》（ GB50141—2008）中關(guān)于滿水試驗(yàn)驗(yàn)收的要求，鋼筋混凝土池體滿水試驗(yàn)驗(yàn)收標(biāo)準(zhǔn)為2.0 L／ m2·d，本次滲漏量按照驗(yàn)收標(biāo)準(zhǔn)的10 倍計(jì)算，即20 L／ m2·d，因此上邊界是變化的濃度通量邊界，前30 d 的通量為2 cm／ d（20 L／ m2·d）；30 d 后的通量為0。 池體中鎳的濃度0.025 mg ／ cm3（25 mg ／ L） 。

ρ根據(jù)土工試驗(yàn)的成果取1 500 mg ／ cm3；DL取包氣帶厚度的十分之一，為21.2 cm。

2.3 包氣帶剖分和模擬時(shí)間

根據(jù)場(chǎng)地水文地質(zhì)調(diào)查結(jié)果，本次模擬土壤類型為一種，包氣帶自上而下均勻布設(shè)6 個(gè)觀測(cè)點(diǎn)，埋深分別為5、40、80、120、160、212 cm，以表明溶質(zhì)在垂向上的分布規(guī)律。

本次模擬時(shí)間為100 d，輸出5 個(gè)時(shí)間節(jié)點(diǎn)（1、5、10、20、100 d）的數(shù)據(jù)，以表明包氣帶剖面上溶質(zhì)隨時(shí)間的運(yùn)動(dòng)變化規(guī)律。

3 模擬結(jié)果及分析

（1）不同深度處鎳濃度隨時(shí)間變化曲線如圖1。

圖1 不同深度處土壤中鎳濃度隨時(shí)間變化曲線

由圖1 可知，0 ～30 d 內(nèi)由于污染物的滲漏，隨著時(shí)間的增加不同深度觀測(cè)點(diǎn)位土壤水中鎳的濃度逐漸升高；30 d 時(shí)鎳停止泄漏，不同深度觀測(cè)點(diǎn)位土壤水中鎳的濃度約在30 d 時(shí)也開(kāi)始趨于穩(wěn)定。其中－160 cm 和－212 cm 處30 d 后仍有緩慢的上升趨勢(shì)，這可能是由于少量毛細(xì)水在重力作用下仍有向下遷移的趨勢(shì)。 同時(shí)可以看出不同深度觀測(cè)點(diǎn)的濃度最終保持穩(wěn)定不變，說(shuō)明即使污染源泄漏后被及時(shí)發(fā)現(xiàn)及時(shí)切斷，但污染物一旦進(jìn)入土壤，在只有對(duì)流-彌散的作用下，會(huì)長(zhǎng)久的存在于包氣帶內(nèi)，造成土壤環(huán)境的污染。

（2）不同時(shí)間包氣帶剖面上土壤水中鎳濃度隨深度變化曲線如圖2。

圖2 剖面上不同時(shí)間土壤中鎳濃度隨深度變化曲線

由圖2 可知，不同時(shí)刻，包氣帶剖面由頂?shù)降?，土壤中鎳的濃度逐漸降低，100 d 時(shí)，頂部濃度為0.0241 mg／ cm3（圖中E 點(diǎn)），底部濃度為0.0149 mg／ cm3（圖中e 點(diǎn)）。 同時(shí)結(jié)合圖1 和表3 可以看出，隨著時(shí)間的遷移，包氣帶剖面上的濃度梯度逐漸減小，污染物遷移的速率逐漸也減小，這是因?yàn)榘鼩鈳Т嬖趦?yōu)勢(shì)通道，前期污染物泄漏后在重力勢(shì)的作用下順著優(yōu)勢(shì)通道遷移，遷移速率較大，后期隨著污染物的持續(xù)泄漏，污染物逐漸填充空隙，遷移速率變小。

表3 濃度梯度和遷移速率計(jì)算表

另外由表3 可知，隨著模擬時(shí)間的增加，鎳逐漸向下遷移，第1 d 鎳遷移的最大距離為53 cm （a點(diǎn)），第5 d 鎳遷移的最大距離為137 cm（b點(diǎn)），第10 d 鎳便遷移到包氣帶底部（c點(diǎn)）。 說(shuō)明對(duì)于滲透性能較好、厚度較小的包氣帶截污能力有限，在只有對(duì)流-彌散的作用下，污染物物一旦發(fā)生泄漏，在較短時(shí)間內(nèi)便會(huì)穿透包氣帶進(jìn)入地下水含水層，從而對(duì)場(chǎng)地內(nèi)潛水產(chǎn)生污染，因此企業(yè)在生產(chǎn)活動(dòng)中應(yīng)對(duì)各個(gè)構(gòu)筑物進(jìn)行嚴(yán)格的防滲措施。

4 結(jié)論

隨著污染源的滲漏和切斷，包氣帶中污染物的濃度逐漸升高并趨于穩(wěn)定。 即使污染源泄漏后被及時(shí)發(fā)現(xiàn)并切斷，但污染物一旦進(jìn)入土壤，只有對(duì)流-彌散的作用下，會(huì)長(zhǎng)久的存在于包氣帶內(nèi)，造成土壤環(huán)境的污染。

在只有對(duì)流-彌散的作用下，隨著污染物的泄漏，污染物逐漸向下遷移并到達(dá)包氣帶底部。 對(duì)于滲透性能較好、厚度較小的包氣帶截污能力有限，污染物一旦發(fā)生泄漏，在較短時(shí)間內(nèi)便會(huì)穿透包氣帶進(jìn)入地下水含水層，從而對(duì)場(chǎng)地內(nèi)潛水產(chǎn)生污染。

企業(yè)在生產(chǎn)活動(dòng)中應(yīng)對(duì)各個(gè)污水處理站各構(gòu)筑物進(jìn)行嚴(yán)格的防滲措施，建立土壤和地下水污染防控體系，及時(shí)發(fā)現(xiàn)并截?cái)辔廴驹?，避免?duì)場(chǎng)地內(nèi)土壤和地下水造成污染。