，，

(南京大學(xué)環(huán)境規(guī)劃設(shè)計(jì)研究院股份公司，江蘇 南京 210093)

隨著我國(guó)經(jīng)濟(jì)的快速發(fā)展，沿海地區(qū)加大了濱海地區(qū)的開(kāi)發(fā)強(qiáng)度，化工、醫(yī)藥、農(nóng)藥等行業(yè)在濱海地區(qū)迅猛發(fā)展，人工填海地區(qū)成為該行業(yè)重要的選擇[1-3]。濱海地區(qū)的過(guò)度開(kāi)發(fā)導(dǎo)致污染問(wèn)題日益突出，并逐漸成為了建設(shè)生態(tài)文明和諧社會(huì)的一個(gè)主要矛盾。尤其是對(duì)于地下水環(huán)境而言，地下水污染具有隱蔽性、滯后性以及一旦污染便難以清除的特點(diǎn)，因此備受關(guān)注[4]。而針對(duì)建設(shè)項(xiàng)目的地下水環(huán)境影響預(yù)測(cè)評(píng)價(jià)，能預(yù)先評(píng)估和預(yù)防項(xiàng)目可能造成的地下水污染，起到源頭防治和管理的作用，對(duì)保護(hù)地下水環(huán)境意義重大[5-7]。此外，濱海地區(qū)地下水、地表水以及海水之間的水力聯(lián)系密切且復(fù)雜，濱海地區(qū)地下水環(huán)境評(píng)價(jià)成為項(xiàng)目環(huán)評(píng)中的難點(diǎn)。這類(lèi)項(xiàng)目地下水環(huán)境評(píng)價(jià)研究較少，目前可參考的文獻(xiàn)非常有限。本文以某化工企業(yè)為例，基于數(shù)值模擬模型研究濱海地區(qū)地下水環(huán)境影響評(píng)價(jià)。

1 研究思路與背景

1.1 研究思路

通過(guò)資料收集和分析，結(jié)合現(xiàn)場(chǎng)水文地質(zhì)環(huán)境調(diào)查進(jìn)行綜合研究，確定合理的模擬范圍；分析研究區(qū)地下水流動(dòng)特征、含隔水層系統(tǒng)結(jié)構(gòu)和邊界條件等，建立研究區(qū)水文地質(zhì)概念模型；通過(guò)對(duì)研究區(qū)進(jìn)行空間離散、水文地質(zhì)參數(shù)賦值、地下水流場(chǎng)識(shí)別與校正，構(gòu)建研究區(qū)地下水水流模型；研究項(xiàng)目污染源特征，確定合理的污染因子，設(shè)計(jì)地下水泄露污染情景，在地下水水流模型的基礎(chǔ)上耦合污染物運(yùn)移模型，對(duì)研究區(qū)地下水進(jìn)行數(shù)值模擬預(yù)測(cè)，選擇合理的水質(zhì)標(biāo)準(zhǔn)，針對(duì)模擬結(jié)果開(kāi)展評(píng)價(jià)[8-11]。

1.2 研究背景

1.2.1 研究區(qū)概況

研究區(qū)為南通某沿海經(jīng)濟(jì)開(kāi)發(fā)區(qū)，累計(jì)落戶項(xiàng)目100多個(gè)，形成集醫(yī)藥、農(nóng)藥、高分子材料為主導(dǎo)的三大產(chǎn)業(yè)板塊群。開(kāi)發(fā)區(qū)東枕黃海，占地面積約9.7km2，區(qū)位條件得天獨(dú)厚。研究區(qū)地貌特征為典型的濱海平原，場(chǎng)地較平緩開(kāi)闊，地形坡度3°以內(nèi)，自然地面標(biāo)高3.79～5.69 m(黃海高程系)，總體呈西南高東北低狀。

1.2.2 水文地質(zhì)條件

本區(qū)大地構(gòu)造處于揚(yáng)子準(zhǔn)地臺(tái)東部，地層屬于揚(yáng)子地層區(qū)。地表第四系廣泛覆蓋，無(wú)基巖出露。區(qū)內(nèi)地下水類(lèi)型主要為松散巖類(lèi)孔隙水，具有分布廣、層次多、水量豐富，水質(zhì)復(fù)雜等特征。區(qū)內(nèi)潛水含水層地層巖性主要為粉土、粉砂，粉土層厚1.35～1.60 m，平均1.51 m，粉砂層厚19.65～25.60 m，平均22.61 m；潛水含水層下隔水層層厚7.00～17.00 m，平均12.06 m，巖性主要為粉質(zhì)粘土，滲透性較弱。區(qū)內(nèi)淺層地下水埋深較淺，一般在1.55～2.83 m左右，地下水流向主要由西南向東北，研究區(qū)地下水流場(chǎng)圖見(jiàn)圖1。

圖1 研究區(qū)地下水流場(chǎng)圖

2 數(shù)值模型

2.1 水文地質(zhì)概念模型

研究區(qū)水文地質(zhì)條件較單一，主要含水層為第四系松散巖類(lèi)含水層，含水層分布較均勻，在潛水含水層和承壓含水層之間有一層良好的粉質(zhì)粘土隔水層，使得潛水含水層與承壓含水層之間水力聯(lián)系較弱。

結(jié)合研究區(qū)水資源分區(qū)、水系分布，考慮區(qū)域地質(zhì)、水文地質(zhì)、環(huán)境水文地質(zhì)條件，研究區(qū)西部以振洋河為河流界，北部以黃海為水頭界，東部以匡河為河流界，南部以一內(nèi)陸河為河流界，整個(gè)研究區(qū)范圍面積約16 km2。含水層上邊界為地面，通過(guò)該邊界含水層系統(tǒng)與大氣降水、地表水等產(chǎn)生垂向上的水量交換；下邊界為透水性差的粉質(zhì)粘土弱透水層，該層阻斷了潛水含水層與下伏承壓含水層的水力聯(lián)系，故定義為隔水邊界。根據(jù)模擬區(qū)地層條件，污染進(jìn)入地下主要污染潛水含水層。

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

根據(jù)水文地質(zhì)概念模型，評(píng)價(jià)范圍內(nèi)地下水流運(yùn)動(dòng)的數(shù)學(xué)模型可以表示為潛水含水層非均質(zhì)、各向異性三維非穩(wěn)定流數(shù)學(xué)模型，其控制方程及定解條件如下：

H(x，y，z，0)=H0(x，y，z)(x，y，z)∈Ω

H(x，y，z，t)|Γ1=H1(x，y，z，t)

式中：Kxx,Kyy,Kzz為主坐標(biāo)軸方向多孔介質(zhì)的滲透系數(shù)，[LT-1]；h為水頭，[L]；w為單位面積垂向流量，[LT-1]，用以表示源匯項(xiàng)；μ為多孔介質(zhì)的給水度(或飽和差)；z為潛水含水層的底板標(biāo)高，[L]；t為時(shí)間，[T]；Ω表示滲流區(qū)域；Γ1表示第一類(lèi)給定水頭邊界。

污染物在地下水中的運(yùn)移包括對(duì)流、彌散以及溶質(zhì)本身的物理、化學(xué)變化等過(guò)程，可表示為：

本文采用GMS軟件求解，用MODFLOW計(jì)算模塊求解地下水水流運(yùn)動(dòng)數(shù)學(xué)模型，用MT3DMS模塊求解地下水污染物運(yùn)移數(shù)學(xué)模型。

2.3 模型離散

本文數(shù)值模擬求解采用有限差分法，需將評(píng)價(jià)區(qū)剖分為四邊形網(wǎng)格，在廠區(qū)污水處理站處進(jìn)行加密處理，最小網(wǎng)格空間長(zhǎng)度達(dá)到2 m，網(wǎng)格垂向上剖分依據(jù)場(chǎng)區(qū)建設(shè)特點(diǎn)以及評(píng)價(jià)區(qū)內(nèi)潛水含水層特征劃分為兩層：第一層為粉土，厚度1.5 m左右；第二層為粉砂，含水層厚度22.6 m左右；在含水層下部設(shè)置一層相對(duì)隔水層，厚度為12.1 m左右；整個(gè)模型在垂向上一共三層。一共剖分148 512個(gè)網(wǎng)格，其中有效網(wǎng)格為89 733個(gè)。詳細(xì)見(jiàn)圖2。

圖2 研究區(qū)三維網(wǎng)格剖分示意圖

2.4 水文地質(zhì)參數(shù)

水文地質(zhì)參數(shù)的選擇主要依據(jù)該項(xiàng)目水文地質(zhì)調(diào)查所進(jìn)行的各類(lèi)野外及室內(nèi)水文地質(zhì)試驗(yàn)的結(jié)果，并結(jié)合以往在附近區(qū)域各類(lèi)水文地質(zhì)試驗(yàn)數(shù)據(jù)資料確定，其水文地質(zhì)參數(shù)取值如表1所示。

表1 研究區(qū)水文地質(zhì)參數(shù)取值

彌散度是研究污染物在地下水中遷移轉(zhuǎn)化規(guī)律的重要參數(shù)之一，由于彌散度的不確定性，本文將通過(guò)對(duì)比不同彌散度取值條件下污染物遷移規(guī)律，幫助企業(yè)科學(xué)管理地下水。

2.5 模型校正與識(shí)別

將數(shù)值模型進(jìn)行求解，對(duì)比模型計(jì)算結(jié)果與實(shí)際觀測(cè)數(shù)據(jù)差異程度，從而校正檢驗(yàn)?zāi)Ｐ偷暮侠硇?、?zhǔn)確性和穩(wěn)定性。

2.5.1 地下水水位擬合

模擬計(jì)算含水層地下水水位與實(shí)測(cè)地下水水位關(guān)系見(jiàn)圖3。圖中可以看出實(shí)際觀測(cè)水位與計(jì)算水位誤差均在0.2 m以內(nèi)，模擬誤差小，說(shuō)明模型是合理的。

圖3 研究區(qū)地下水水位擬合圖

2.5.2 水均衡

模擬計(jì)算得到的模擬范圍內(nèi)水均衡結(jié)果如表2所示。

表2 研究區(qū)水均衡結(jié)果表 m3/a

根據(jù)水均衡結(jié)果，評(píng)價(jià)區(qū)每年地下水排泄進(jìn)入地表水1 757 080.8 m3，地表水補(bǔ)給地下水的量為785 582.2 m3，表明地下水和地表水存在密切的水力聯(lián)系。綜上，根據(jù)對(duì)地下水水位及水均衡計(jì)算結(jié)果的分析，模型能較好反映該地區(qū)地下水流運(yùn)動(dòng)特征，可以用于地下水環(huán)境影響的預(yù)測(cè)評(píng)價(jià)。

3 地下水模擬預(yù)測(cè)與評(píng)價(jià)

3.1 地下水污染源

建設(shè)項(xiàng)目工程可能存在的地下水污染源主要是各生產(chǎn)裝置產(chǎn)生的生產(chǎn)廢水、各生產(chǎn)廢水的輸送管線、污水處理池等，根據(jù)對(duì)本次改擴(kuò)建項(xiàng)目的生產(chǎn)工藝、設(shè)備設(shè)施分析后，項(xiàng)目裝置區(qū)內(nèi)的污水處理池防治措施較薄弱，也是最重要的地下水污染源。非正常狀況下污水處理池池壁因系統(tǒng)老化、防滲措施失效等原因發(fā)生滲漏可能對(duì)地下水產(chǎn)生影響。本工程排放的廢水主要產(chǎn)生于噻蟲(chóng)胺、乙螨唑、蟲(chóng)胺等原藥生產(chǎn)過(guò)程中的廢水。污水處理池廢水源強(qiáng)進(jìn)水情況見(jiàn)表3。

表3 污水處理池進(jìn)水源強(qiáng)情況表

非正常狀況下，污水處理池發(fā)生滲漏，廢水經(jīng)包氣帶進(jìn)入潛水含水層。挑選污水處理池調(diào)節(jié)池底部及側(cè)壁面積約為220 m2，滲漏面積按面積的5‰計(jì)算，根據(jù)《給水排水構(gòu)筑物工程施工及驗(yàn)收規(guī)范》(GB50141-2008)，鋼筋混凝土結(jié)構(gòu)水池滲水量不得超過(guò)2 L/(m2·d)，非正常狀況按照正常狀況的100倍考慮，則污水處理池調(diào)解池滲水量為0.22 m3/d。

3.2 預(yù)測(cè)因子及時(shí)段

項(xiàng)目地下水的主要污染源來(lái)自污水處理池中的廢水，根據(jù)對(duì)項(xiàng)目廢水的主要成分分析，確定廢水中主要種類(lèi)和標(biāo)準(zhǔn)值屬較大的COD和甲苯作為預(yù)測(cè)的對(duì)象，此類(lèi)物質(zhì)有較大危害性，如進(jìn)入地下水易對(duì)環(huán)境造成危害，因此選取該指標(biāo)進(jìn)行地下水溶質(zhì)遷移預(yù)測(cè)。

考慮項(xiàng)目建設(shè)、運(yùn)營(yíng)和退役期，將地下水環(huán)境影響預(yù)測(cè)時(shí)段擬定為10 000 d。結(jié)合工程特征與環(huán)境特征，預(yù)測(cè)污染發(fā)生100 d、1 000 d及10 000 d后污染物遷移情況。

3.3 評(píng)價(jià)標(biāo)準(zhǔn)

COD標(biāo)準(zhǔn)限值參照《地下水質(zhì)量標(biāo)準(zhǔn)》(GB/T14848-93)Ⅲ類(lèi)標(biāo)準(zhǔn)值[12]，甲苯超標(biāo)范圍參照《地下水質(zhì)標(biāo)準(zhǔn)》(DZ/T0290-2015)Ⅲ類(lèi)標(biāo)準(zhǔn)值[13]。因COD與高錳酸鉀指數(shù)存在換算系數(shù)，根據(jù)南通地區(qū)地表水監(jiān)測(cè)數(shù)據(jù)顯示換算系數(shù)在2.5～3之間，為保守起見(jiàn)，本文COD濃度與高錳酸鉀指數(shù)按4倍進(jìn)行換算。各指標(biāo)具體情況見(jiàn)表4。

表4 地下水污染物評(píng)價(jià)標(biāo)準(zhǔn)

圖4 不同時(shí)刻下COD運(yùn)移分布圖(D=10m)

3.4 預(yù)測(cè)結(jié)果分析

3.4.1 彌散度為10 m時(shí)，污染物遷移規(guī)律

將污染物泄露量和水文地質(zhì)參數(shù)帶入校正模型進(jìn)行模擬預(yù)測(cè)，彌散度為10 m時(shí)，100 d、1 000 d、10 000 d時(shí)刻下COD和甲苯運(yùn)移規(guī)律見(jiàn)圖4和圖5。

從圖中可以看出，隨著時(shí)間不斷增加，污染物在中心處的濃度逐漸升高，往東北方向不斷遷移，距離逐漸增加，污染范圍不斷擴(kuò)大；污染物詳細(xì)遷移規(guī)律見(jiàn)表5。

表5 D=10 m時(shí)不同時(shí)刻污染物運(yùn)移特征表

3.4.2 不同彌散度條件下污染物遷移規(guī)律

溶質(zhì)運(yùn)移模型的彌散度參數(shù)具有“尺度效應(yīng)”(隨研究尺度的增大而增大)，因此彌散度具有很大的不確定性。為了更好的了解污染物可能的運(yùn)移規(guī)律，可以通過(guò)分析不同彌散度條件下污染物遷移規(guī)律，得到污染物遷移距離、污染范圍、污染物濃度等參數(shù)，從而方便企業(yè)后期科學(xué)管理判斷可能存在的污染現(xiàn)象。本文選取彌散度為2 m、10 m、50 m進(jìn)行模擬預(yù)測(cè)，COD和甲苯模擬結(jié)果分別見(jiàn)圖6和圖7。

圖5 不同時(shí)刻下甲苯運(yùn)移分布圖(D=10m)

圖6 不同彌散度下COD運(yùn)移分布圖(T=10 000 d)

圖7 不同彌散度下甲苯運(yùn)移分布圖(T=10 000 d)

圖8 不同彌散度下污染物中心濃度變化趨勢(shì)圖

圖9 不同彌散度下遷移距離與超標(biāo)范圍變化趨勢(shì)圖

模擬結(jié)果顯示，模擬時(shí)間為10 000 d時(shí)，隨著彌散度的不斷增加，污水處理池中心濃度逐漸下降，遷移距離和超標(biāo)范圍不斷擴(kuò)大，見(jiàn)圖8和圖9。彌散度為2 m、10 m、50 m時(shí)，污染物COD和甲苯中心處濃度與初始濃度的比值從分別為0.97、0.94、0.90和0.96、0.91、0.87，遷移距離分別為32.39 m、52.12 m、88.24 m和36.89 m、58.38 m、102.49 m，超標(biāo)范圍分別為2 263.3 m2、4 007.4 m2、10 080.3 m2和2 585.5 m2、4 763.1 m2、12 851.1 m2。COD比甲苯的濃度下降的更慢，而甲苯比COD的遷移距離和超標(biāo)范圍更大。

發(fā)生這一現(xiàn)象的主要原因是，彌散度的增加有利于污染物的擴(kuò)散，同時(shí)隨著對(duì)流作用，污染物越容易向外擴(kuò)散，從而導(dǎo)致遷移距離不斷增加，污染物擴(kuò)散范圍及污染范圍變大；與此同時(shí)，中心處的污染物濃度隨著擴(kuò)散作用的增強(qiáng)會(huì)逐漸下降。

由于彌散度的尺度效應(yīng)，在目前無(wú)法有效解決尺度效應(yīng)情況下，通過(guò)模擬不同彌散度下污染物遷移規(guī)律，確定彌散度對(duì)溶質(zhì)運(yùn)移的影響。與確定的彌散度模型相比，本文可以更好的弄清污染物可能的遷移距離和超標(biāo)范圍，有利于企業(yè)加強(qiáng)地下水管理和判斷可能發(fā)生的污染現(xiàn)象。

4 結(jié)語(yǔ)

(1)本文在對(duì)濱海地區(qū)水文地質(zhì)條件詳細(xì)了解和調(diào)查的情況下，構(gòu)建了地下水?dāng)?shù)值模型，選擇合理的水文地質(zhì)參數(shù)校正模型，所建立的模型能夠真實(shí)地反應(yīng)研究區(qū)實(shí)際的水文地質(zhì)條件。

(2)針對(duì)本項(xiàng)目特點(diǎn)，設(shè)計(jì)了合理的污水處理池泄露情景，分析了非正常工況下污染物對(duì)地下水的影響，模擬預(yù)測(cè)了污染物的時(shí)空遷移與擴(kuò)散規(guī)律，結(jié)果表明隨時(shí)間的增加，中心處污染物濃度逐漸升高，遷移距離逐漸增加，污染范圍不斷擴(kuò)大。

(3)由于彌散度的尺度效應(yīng)，研究了不同彌散度條件下污染物遷移規(guī)律，結(jié)果表明彌散度的不斷增加，污水處理池中心濃度逐漸下降，遷移距離和超標(biāo)范圍不斷擴(kuò)大。通過(guò)模擬不同彌散度下污染物遷移規(guī)律，有利于企業(yè)加強(qiáng)地下水管理和判斷可能發(fā)生的污染現(xiàn)象。

[1]薛禹群，謝春紅，吳吉春，等．山東龍口-萊州地區(qū)的海水入侵[J]．地質(zhì)學(xué)報(bào).1992.66(3):280-291.

[2]薛禹群，謝春紅，吳吉春．海水入侵研究[J]．水文地質(zhì)工程地質(zhì).1992.19(6):29-33.

[3]左銳，王金生，楊潔，等．濱海石化項(xiàng)目地下水環(huán)境影響評(píng)價(jià)的關(guān)鍵問(wèn)題[J]．水文地質(zhì)工程地質(zhì).2010.37(3)：97-101.

[4]余期沖，祝曉彬，吳吉春，等．死端孔隙對(duì)溶質(zhì)運(yùn)移影響的實(shí)驗(yàn)研究[J]．水文地質(zhì)工程地質(zhì).2017.44(4)：160-164.

[5]盧曉華．基于數(shù)值模擬的企業(yè)地下水重金屬污染的環(huán)境影響預(yù)測(cè)評(píng)價(jià)[J]．安全與環(huán)境工程.2014.21(1)：93-97.

[6]鄭玉虎，吳明洲，徐愛(ài)蘭，等．考慮土壤吸附作用的地下水污染物運(yùn)移特征研究[J]．地下水.2017.39(3)：4-7.

[7]岳強(qiáng)，范亞民，耿磊，等．地下水環(huán)境影響評(píng)價(jià)導(dǎo)則執(zhí)行過(guò)程中遇到的問(wèn)題及建議[J]．環(huán)境科學(xué)與管理.2012.37(10):174-177.

[8]薛禹群，吳吉春．地下水動(dòng)力學(xué)[M]．北京:水利水電出版社．2009．

[9]王焰新．地下水污染與防治[M]．北京：高等教育出版社.2007.

[10]朱學(xué)愚，錢(qián)孝星．地下水環(huán)境影響評(píng)價(jià)的工作要點(diǎn)[J]．水資源保護(hù).1998(4):48-53.

[11]余期沖，祝曉彬，吳吉春，等．場(chǎng)地地下水源熱泵系統(tǒng)抽灌實(shí)驗(yàn)及數(shù)值模擬研究[J]．地下水.2015.37(3)：24-27.

[12]GB/T14848-93 地下水質(zhì)量標(biāo)準(zhǔn) [S].

[13]DZ/T0290-2015 地下水質(zhì)標(biāo)準(zhǔn) [S].