付震，張福超，安禹辰

（中水北方勘測(cè)設(shè)計(jì)研究有限責(zé)任公司，天津 300222）

工業(yè)拆遷遺留區(qū)域的土壤污染問(wèn)題日益嚴(yán)重，特別是“水十條”“土十條”頒布以前，部分地區(qū)的一些小型工業(yè)企業(yè)私排偷排嚴(yán)重，很多滲坑滲井不易發(fā)現(xiàn)，成為土壤污染治理的重點(diǎn)難點(diǎn)。南水北調(diào)中線工程線路較長(zhǎng)，天津干線工程穿越華北平原中部。王慶坨水庫(kù)是天津干線的重要調(diào)節(jié)水庫(kù)，占地范圍在歷史上存在多處工業(yè)企業(yè)。水庫(kù)施工過(guò)程中，庫(kù)區(qū)表土剝離后局部區(qū)域發(fā)現(xiàn)了污染土壤的痕跡。掌握工業(yè)遺留區(qū)土壤污染特性并提出處置措施，對(duì)保障施工與水庫(kù)水質(zhì)安全具有重要意義。

傳統(tǒng)的大范圍鉆孔、多點(diǎn)取樣的監(jiān)測(cè)手段費(fèi)時(shí)費(fèi)力，影響工期，同時(shí)極易造成庫(kù)區(qū)底部黏土層穿透風(fēng)險(xiǎn)升高。相比傳統(tǒng)方法，數(shù)學(xué)模型可以模擬土壤中污染物溶質(zhì)運(yùn)移過(guò)程，是土壤污染特性研究的重要方法。HYDRUS-1D軟件是用于模擬飽和-非飽和多孔介質(zhì)中水分運(yùn)移和溶質(zhì)運(yùn)移的數(shù)值模型，在土壤污染物運(yùn)移、地下水污染風(fēng)險(xiǎn)評(píng)價(jià)方面得到了廣泛運(yùn)用。

為確保王慶坨水庫(kù)水質(zhì)安全，在盡量不影響施工工期、避免庫(kù)區(qū)滲漏的前提下，利用地質(zhì)勘察、污染源調(diào)查、HYDRUS模擬等綜合手段，構(gòu)建了一套工業(yè)滲坑土壤重金屬污染模擬方法，并依據(jù)模擬結(jié)果制定處置方案。方案實(shí)施后，取得了顯著環(huán)境效益，為新建水庫(kù)土壤重金屬污染處置提供了借鑒。

1 研究區(qū)域概況

王慶坨水庫(kù)屬新建中型平原水庫(kù)，位于華北平原中北部，庫(kù)區(qū)面積4.21 km2，調(diào)節(jié)庫(kù)容1 500萬(wàn)m3，總庫(kù)容2 000萬(wàn)m3。作為南水北調(diào)天津干線的事故和維修備用水庫(kù)，遇特殊供水事故時(shí)，作為城市備用水源，同時(shí)調(diào)節(jié)南水北調(diào)來(lái)水的不均勻性，保障供水穩(wěn)定。水庫(kù)呈較為規(guī)則的四邊形，占?jí)航ㄔO(shè)用地4.87 hm2，拆遷18家企業(yè)和1家單位。工程于2017年開(kāi)工，2019年完工并投入運(yùn)行。

2 研究方法

2.1 地質(zhì)勘察

施工中發(fā)現(xiàn)污染土壤痕跡后，立即開(kāi)展應(yīng)急修復(fù)，首先進(jìn)行地質(zhì)勘察、現(xiàn)場(chǎng)標(biāo)準(zhǔn)貫入試驗(yàn)及室內(nèi)常規(guī)物理力學(xué)試驗(yàn)。地質(zhì)勘察結(jié)果表明，按地表以下10～18 m范圍內(nèi)地層組合劃分地質(zhì)結(jié)構(gòu)，污染區(qū)為砂、黏雙層結(jié)構(gòu)，上部為全新統(tǒng)上段陸相沖積砂壤土及粉土、粉砂，厚度0.5～3.1 m；下部為全新統(tǒng)上段陸相沖積及全新統(tǒng)中段湖沼相沉積壤土和黏土，厚度1.4～5.6 m。將試驗(yàn)數(shù)據(jù)按地層時(shí)代、成因類型、巖性分類進(jìn)行統(tǒng)計(jì)，土的物理力學(xué)指標(biāo)詳見(jiàn)表1。

表1 土的物理力學(xué)指標(biāo)建議值

2.2 污染源調(diào)查

經(jīng)調(diào)查，水庫(kù)庫(kù)底污染區(qū)原為小型電鍍企業(yè)，約建于20世紀(jì)90年代末，2016年拆除，持續(xù)生產(chǎn)15～16 a。生產(chǎn)產(chǎn)品為自行車金屬加工零件，詳細(xì)生產(chǎn)規(guī)模、生產(chǎn)工藝、污水排放量已經(jīng)無(wú)法調(diào)查。根據(jù)水庫(kù)現(xiàn)場(chǎng)污染分布狀況，可憑經(jīng)驗(yàn)判斷該企業(yè)歷史上存在滲坑或滲井偷排現(xiàn)象。

水庫(kù)庫(kù)區(qū)施工開(kāi)挖后，大氣降雨使土壤污染區(qū)短期內(nèi)殘存有地表積水，相當(dāng)于土壤淋洗溶液，渾濁呈紅色。現(xiàn)場(chǎng)對(duì)污染區(qū)土壤及水質(zhì)進(jìn)行環(huán)境監(jiān)測(cè)，土壤取樣深度為表層0.2 m及下層2 m，水質(zhì)取樣為地表積水，監(jiān)測(cè)結(jié)果顯示：①地表積水強(qiáng)酸性pH值2.78，鐵離子濃度已超過(guò)3×104mg/L；②土壤中鐵含量達(dá)到48.8～76 g/kg。根據(jù)以上結(jié)果，可判斷該區(qū)域已被工業(yè)滲坑污染，且污染持續(xù)時(shí)間長(zhǎng)、積累嚴(yán)重。

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

2.3.1 方法與目的

掌握污染區(qū)的地質(zhì)、表層土壤及水質(zhì)監(jiān)測(cè)數(shù)據(jù)后，尚無(wú)法對(duì)工業(yè)滲坑的污染范圍和程度作出科學(xué)判斷。采用HYDRUS模型，模擬分析工業(yè)污水排放和污染溶質(zhì)在土壤介質(zhì)中的運(yùn)移。一方面，為污染處置提供科學(xué)依據(jù)，縮短因應(yīng)急處置而停擺的施工時(shí)間，保障工程順利實(shí)施；另一方面，避免大范圍、多鉆孔取樣，減少庫(kù)區(qū)底部黏土層的穿透風(fēng)險(xiǎn)，減輕水庫(kù)建成后的污染風(fēng)險(xiǎn)。

2.3.2 水動(dòng)力學(xué)模型

污水在土壤中的下滲采用修改的Richards對(duì)流-彌散方程和迦遼金有限元法進(jìn)行數(shù)值求解?；痉匠虨椋?/p>

式中：t為時(shí)間（d）；h為負(fù)壓水頭（m）；θ為包氣帶體積含水率（%）；K為水力傳導(dǎo)系數(shù)；S為根系吸水項(xiàng)，在本處為0；θe為包氣帶有效含水率（%）；θr為包氣帶殘余含水率（%）；θs為包氣帶飽和含水率（%）；Ks為滲透系數(shù)；n，m，α均為經(jīng)驗(yàn)參數(shù)。初始邊界條件為：

式中：z為空間坐標(biāo)，原點(diǎn)在地面，向下為正；h0為初始包氣帶剖面負(fù)壓值（cm）；ε1為庫(kù)底補(bǔ)給強(qiáng)度（cm/d）；l為地下水埋深（cm）；其余變量含義同上。

2.3.3 溶質(zhì)運(yùn)移模型

溶質(zhì)運(yùn)動(dòng)描述采用CDE對(duì)流-彌散方程，基本方程為：

式中：D為水動(dòng)力彌散系數(shù)；c為土壤溶液中污染質(zhì)濃度（g/cm3）；s為被吸附的固相濃度（%）；q為土體中水的流速（cm/s）；ρ0，ρb為流體密度（g/cm3）；λ1，λ2為經(jīng)驗(yàn)常數(shù)，與包氣帶土壤質(zhì)地和結(jié)構(gòu)有關(guān)；其余變量含義同上。初始條件與邊界條件為：

式中：c0為初始土壤剖面污染質(zhì)濃度（g/cm3）；cm為土壤溶液污染質(zhì)飽和濃度（g/cm3）；cd為地下水中的污染質(zhì)濃度（g/cm3）；其余變量含義同上。

2.4 參數(shù)和邊界條件

模型參數(shù)賦值在地勘實(shí)測(cè)及參考文獻(xiàn)中獲得。以金屬加工工業(yè)酸洗工藝特征污染物Fe2+離子為模擬溶質(zhì)，假設(shè)原工業(yè)企業(yè)在一個(gè)主汛期中連續(xù)排放污水，模擬地表以下20 m土壤從空間及時(shí)間上的污染進(jìn)程。如有必要，調(diào)整模型將模擬時(shí)間尺度最大放寬至15 a。由于項(xiàng)目為平原新建水庫(kù)，地下黏土隔水層對(duì)蓄水后的滲漏、污染阻隔至關(guān)重要，因此在隔水層底部增加觀察點(diǎn)。模型主要參數(shù)詳見(jiàn)表2，邊界及初始條件詳見(jiàn)表3。

表2 土壤水分和溶質(zhì)運(yùn)移參數(shù)

表3 模型邊界及初始條件

3 結(jié)果與分析

假設(shè)歷史上的企業(yè)污水毫無(wú)處理、滲坑連續(xù)偷排，在此最不利條件下，模擬顯示在一個(gè)主汛期內(nèi)，僅需8 d滲坑地表土壤Fe2+離子濃度已超過(guò)1×104mg/m3，主汛期末期地表土壤Fe2+離子濃度已達(dá)3.5×104g/m3，如圖1所示。地表濃度模擬值與應(yīng)急環(huán)境實(shí)測(cè)值相近，說(shuō)明模擬擬合效果較好。水庫(kù)工程最為關(guān)注的原狀土壤黏土層底部，28 d時(shí)被滲坑污染物擊穿，隨著汛期降雨及滲坑偷排的連續(xù)作用，到主汛期末期，黏土層底部Fe2+離子濃度預(yù)計(jì)可達(dá)1×103g/m3，如圖2所示。在地勘鉆孔的20 m土層底部，53 d時(shí)被滲坑污染物擊穿，主汛期末期底部Fe2+離子濃度超過(guò)0.1 g/m3，如圖3所示。整個(gè)土壤剖面污染分布以地表以下10 m為主，濃度由上至下逐步變小，且隨著地表以下深度的增加污染物濃度削減速率逐漸增快，如圖4所示。

圖1 土壤表層濃度時(shí)間關(guān)系曲線

圖2 黏土層底部濃度時(shí)間關(guān)系曲線

圖3 土壤底層濃度時(shí)間關(guān)系曲線

圖4 不同排放周期土壤剖面濃度曲線

4 處置方案

依據(jù)前述HYDRUS模擬結(jié)果，必須對(duì)污染土壤采取處置措施。由于該區(qū)域是飲用水水庫(kù)，水質(zhì)安全十分重要，處置措施選擇污染區(qū)整體換填方案。清理范圍以現(xiàn)有污染區(qū)域?yàn)橹行模蛲鈬由?，? m深度逐層清理。每層清挖后，作業(yè)面必須擱置至少3～4 h，觀察是否發(fā)生顏色變化（注：富含鐵的土壤接觸空氣后，會(huì)由黑色變?yōu)辄S、紅色）。如再無(wú)顏色變化，即可停止清除。清理最大深度不超過(guò)天然黏土層底部，一般清挖控制深度小于7 m。清理后換填4 m厚黏土和3 m厚壤土，機(jī)械夯實(shí)。在污染土方處置的同時(shí)，做好地下水倒排及處理工序。

清理后的污染土壤運(yùn)至指定棄土場(chǎng)，摻加石灰調(diào)節(jié)pH值至堿性。棄土場(chǎng)底部及邊坡鋪設(shè)人工防滲膜，棄土完畢后進(jìn)行覆土，覆土厚度不低于1 m，覆土完畢后撒播草籽。工程完畢后，棄土場(chǎng)嚴(yán)格避免用于農(nóng)業(yè)生產(chǎn)，可在自然凈化一段時(shí)間后，用于制磚、建筑等工程性處置。

5 結(jié)論

掌握工業(yè)拆遷遺留土壤污染特性并提出應(yīng)急處置措施對(duì)確保水庫(kù)水質(zhì)安全具有重要意義。利用地質(zhì)勘察、環(huán)境監(jiān)測(cè)、污染源分析、HYDRUS模型等綜合手段，對(duì)平原區(qū)新建水庫(kù)的污染土壤特性進(jìn)行了研究。結(jié)果表明，在研究區(qū)域特定土壤、降雨條件下，工業(yè)企業(yè)特征污染物在一個(gè)主汛期內(nèi)可擴(kuò)散至黏土層以下，且隨著地表以下深度的增加污染物濃度削減速率逐漸增快，水庫(kù)區(qū)域土壤已被污染。在不破壞庫(kù)底隔水層、確保水庫(kù)無(wú)滲漏風(fēng)險(xiǎn)的前提下，依據(jù)模擬結(jié)果開(kāi)展了以換填土為核心的修復(fù)方案，確保了南水北調(diào)水質(zhì)安全。實(shí)踐表明，運(yùn)用HYDRUS模型數(shù)值模擬預(yù)測(cè)方法能較好地預(yù)測(cè)工業(yè)拆遷遺留土壤污染特性，較傳統(tǒng)的大范圍鉆孔取樣、分層環(huán)境監(jiān)測(cè)的方式時(shí)效性較好，對(duì)減少水庫(kù)滲漏風(fēng)險(xiǎn)、縮短應(yīng)急反應(yīng)時(shí)間、確保施工順利進(jìn)行等均較有利。