基于GMS的某典型石化園區(qū)地下水污染運(yùn)移研究

劉景蘭，李立偉，朱明奕，葛菲媛，秦 磊

(1.天津市地質(zhì)研究和海洋地質(zhì)中心，天津 300170； 2.天津市生態(tài)環(huán)境科學(xué)研究院，天津 300191)

自1984年以來(lái)，我國(guó)石化園區(qū)經(jīng)歷了從有到無(wú)，再發(fā)展壯大的過(guò)程。1990年全國(guó)僅有省級(jí)以上園區(qū)13個(gè)，2014年則發(fā)展至近700家，且其中半數(shù)以上分布于我國(guó)沿海地區(qū)。隨著我國(guó)石油化工行業(yè)的迅速發(fā)展，石化園區(qū)對(duì)周邊水土環(huán)境構(gòu)成的污染風(fēng)險(xiǎn)也成為了管理部門和學(xué)界關(guān)注的焦點(diǎn)。石化園區(qū)中的生產(chǎn)活動(dòng)，往往會(huì)造成較為嚴(yán)重的水土污染問(wèn)題。近些年，隨著我國(guó)石化園區(qū)環(huán)境管理制度的不斷完善，園區(qū)生產(chǎn)所產(chǎn)生的廢水、廢氣和廢渣逐步實(shí)現(xiàn)了統(tǒng)一管理，導(dǎo)致污水處理廠和固廢暫存場(chǎng)地等區(qū)域成為了石化園區(qū)的潛在污染源。園區(qū)污水處理設(shè)施在非正常狀況下可能發(fā)生滲漏事故，在設(shè)施檢修的間隔期，污染物滲入地下水對(duì)周邊環(huán)境所構(gòu)成的風(fēng)險(xiǎn)尤其值得關(guān)注，針對(duì)非正常狀況下污染物在地下水中遷移過(guò)程的預(yù)測(cè)對(duì)園區(qū)環(huán)境管理和園區(qū)周邊敏感目標(biāo)的保護(hù)具有重要的現(xiàn)實(shí)意義。

隨著計(jì)算機(jī)模擬技術(shù)的發(fā)展，數(shù)值模擬方法被廣泛用于地下水污染物遷移的模擬預(yù)測(cè)。目前國(guó)內(nèi)外最為常用的數(shù)值模擬軟件包括GMS(Groundwater Modeling System)、Visual MODFLOW和FEFLOW等，其中GMS和Visual MODFLOW是集成MODFLOW、MODPATH和MT3DMS等模擬程序開發(fā)的商業(yè)軟件，這些模擬程序都是基于有限差分方法的，而FEFLOW軟件則是基于有限元方法開發(fā)的，這些模擬軟件為地下水污染物遷移轉(zhuǎn)化的預(yù)測(cè)研究提供了有效工具。翟勇等[1]利用GMS軟件對(duì)我國(guó)東北某化工場(chǎng)地地下水中揮發(fā)酚類物質(zhì)在未來(lái)10年中的遷移和污染物濃度衰減過(guò)程進(jìn)行了模擬預(yù)測(cè)。何華燕等[2]利用FEFLOW軟件對(duì)氯代烴類污染物在地下水中的稀釋擴(kuò)散過(guò)程進(jìn)行了模擬，發(fā)現(xiàn)鹵代烴類污染物遷移能力較弱，擴(kuò)散范圍較小，對(duì)周邊環(huán)境影響有限。魏傳云等[3]利用Visual MODFLOW軟件的MT3DMS模塊，對(duì)某化工園區(qū)污染物在不同情景下的運(yùn)移規(guī)律進(jìn)行了模擬預(yù)測(cè)，并重點(diǎn)針對(duì)連續(xù)源強(qiáng)和瞬態(tài)源強(qiáng)條件下污染物遷移的差異性進(jìn)行了對(duì)比分析。姜鳳成等[4]采用數(shù)值模擬方法對(duì)某化工場(chǎng)地石油類污染物和重金屬Zn在地下水中的時(shí)空演變規(guī)律進(jìn)行了對(duì)比分析，發(fā)現(xiàn)石油類污染物泄漏后在含水層中可發(fā)生自然衰減，污染物濃度和范圍隨時(shí)間逐漸減小，而重金屬Zn會(huì)隨地下水不斷向下游遷移。王平等[5]針對(duì)某有機(jī)污染場(chǎng)地地下水中污染物對(duì)周邊擬建水庫(kù)的潛在影響進(jìn)行了數(shù)值模擬分析，發(fā)現(xiàn)水庫(kù)建設(shè)會(huì)改變周邊地下水流場(chǎng)，從而使得地下水污染羽范圍有所擴(kuò)大，但并不會(huì)對(duì)水庫(kù)水質(zhì)造成影響。張鳳娟等[6]對(duì)淄博市張店化工園區(qū)建設(shè)對(duì)地下水環(huán)境的影響進(jìn)行了模擬預(yù)測(cè)，發(fā)現(xiàn)在20年中石油類污染物在地下水中遷移了350 m的距離，影響面積達(dá)到3萬(wàn)m2，對(duì)地下水環(huán)境影響嚴(yán)重，有必要采取有效的防滲和監(jiān)測(cè)措施加以防范。鄧一榮等[7]利用GMS軟件的MT3D模塊對(duì)粵港澳大灣區(qū)某典型化工污染場(chǎng)地苯系污染物在地下水中的遷移過(guò)程進(jìn)行了模擬分析，發(fā)現(xiàn)溶解態(tài)苯系污染物沿地下水流向往下游遷移，通過(guò)識(shí)別檢驗(yàn)發(fā)現(xiàn)運(yùn)移過(guò)程主要受對(duì)流彌散作用的控制。

本文以我國(guó)北方沿海某石化園區(qū)為研究區(qū)，以園區(qū)內(nèi)污水處理廠作為潛在污染源，以石油烴為特征污染物，利用GMS軟件的MODFLOW和MT3DMS模塊建立石油烴三維遷移數(shù)值模型，對(duì)非正常狀況下污染物滲入地下水后在未來(lái)30年中的遷移過(guò)程進(jìn)行模擬預(yù)測(cè)，研究可為我國(guó)沿海地區(qū)石化園區(qū)環(huán)境管理提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)范圍如圖1所示，該石化產(chǎn)業(yè)園區(qū)建立于2003年7月，產(chǎn)業(yè)結(jié)構(gòu)以精細(xì)化工和醫(yī)藥制造為主，占地面積約746.84 hm2。研究區(qū)地貌類型屬海積沖積低平原，地勢(shì)平坦，地面坡度小于1/10 000。該石化園區(qū)被河流圍繞，北側(cè)為城排河，南側(cè)依次為荒地排河和獨(dú)流減河，西側(cè)城排明渠自北向南溝通城排河和獨(dú)流減河，東側(cè)是長(zhǎng)青河。

圖1 研究區(qū)范圍和水系分布示意Fig.1 Schematic of range and water system distribution of the study area

1.2 研究區(qū)水文地質(zhì)條件

研究區(qū)場(chǎng)地位于我國(guó)北方海積低平原的咸水及鹽鹵水分布區(qū)。根據(jù)區(qū)域水文地質(zhì)調(diào)查成果[8-10]，研究區(qū)400 m以淺的平原松散地層孔隙水可劃分為4個(gè)含水巖組：第Ⅰ含水組屬于淺層地下水系統(tǒng)；第Ⅱ—Ⅳ含水組屬深層地下水系統(tǒng)。第Ⅰ含水組為咸水及鹽鹵水，底板埋深約100 m，其中頂部潛水含水層隔水底板埋深在18 m左右；第Ⅱ含水組在調(diào)查評(píng)價(jià)區(qū)內(nèi)為咸水，咸水底界埋深在200 m以深；深部的Ⅲ、Ⅳ含水組均為承壓淡水。第Ⅰ含水組屬濱海平原沖海積層咸水及鹽鹵水區(qū)，該部分地下水無(wú)開發(fā)利用情況，不作為居民生活飲用水使用[11-12]。因此，研究將地下水環(huán)境保護(hù)目標(biāo)設(shè)置為第Ⅰ含水巖組的頂部潛水含水層，以防止污染物進(jìn)入該潛水含水層向下游遷移擴(kuò)散對(duì)園區(qū)外部地下水和地表水體造成污染，并對(duì)下游敏感目標(biāo)構(gòu)成環(huán)境風(fēng)險(xiǎn)。

第Ⅰ含水巖組頂部含水層地下水水位埋深在0.24～1.29 m，平均水位埋深為0.70 m。包氣帶主要為人工填土層素填土、粉質(zhì)黏土；含水層為第四系松散巖類，巖性以粉質(zhì)黏土、淤泥質(zhì)黏土為主。潛水含水層厚度為15.01～18.71 m，淤泥質(zhì)黏土單層厚度2.90～9.60 m，粉質(zhì)黏土單層厚度0.90～12.70 m。含水層滲透系數(shù)為0.05～0.20 m/d，單位涌水量為0.015～0.270 m3/(h·m)，富水性弱。地下水化學(xué)類型為Cl-Na型水，pH值為7.45～7.66，礦化度為15～54 g/L。

1.3 概念模型建立

1.3.1 含水層結(jié)構(gòu)

模擬區(qū)內(nèi)主要分布第四系松散巖類孔隙水，含水巖組由上至下分別為人工堆積層、全新統(tǒng)上組陸相沖積層和全新統(tǒng)中組淺海相沉積層，巖性由上至下依次為素填土、粉質(zhì)黏土、粉土和淤泥質(zhì)黏土。該潛水含水巖組富水性較差，其下伏粉質(zhì)黏土厚度較大，分布連續(xù)，且水平滲透系數(shù)和垂直滲透系數(shù)在10-7～10-8cm/s，構(gòu)成了地下水潛水良好的隔水底板。因此，本次將第四系松散巖類孔隙水概化為統(tǒng)一的潛水含水層，同時(shí)考慮垂向上的巖性變化，構(gòu)成非均質(zhì)各向異性的多層結(jié)構(gòu)地下水潛水三維模型。

1.3.2 邊界條件

根據(jù)模擬區(qū)水文地質(zhì)條件及地下水流場(chǎng)特征，模擬區(qū)內(nèi)地下水由區(qū)域中心分別向北西、南西和南東方向的河流進(jìn)行徑流排泄，幾乎無(wú)地下水的側(cè)向流入補(bǔ)給，因此將模型側(cè)向邊界均概化為定水頭邊界。模型的上邊界為潛水面，整個(gè)含水層系統(tǒng)通過(guò)該邊界接收大氣降水補(bǔ)給和蒸發(fā)排泄，從而與外界進(jìn)行垂向上的水量交換。

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

由于區(qū)內(nèi)潛水富水性較差，且水質(zhì)普遍咸化，園區(qū)內(nèi)企業(yè)幾乎不對(duì)區(qū)內(nèi)地下水潛水進(jìn)行開發(fā)利用。因此，模擬過(guò)程中無(wú)需考慮人工開采地下水所導(dǎo)致的地下水潛水位的動(dòng)態(tài)變化?；谏鲜鲈颍舜嗡鶚?gòu)建地下水流數(shù)值模型為穩(wěn)定流模型。

地下水穩(wěn)定流模型所需率定的水文地質(zhì)參數(shù)主要為滲透系數(shù)K，根據(jù)所收集的區(qū)內(nèi)抽水試驗(yàn)計(jì)算結(jié)果，結(jié)合場(chǎng)地含水層結(jié)構(gòu)特征，確定潛水含水層水平滲透系數(shù)Kh為0.05～0.20 m/d，將垂直滲透系數(shù)Kv設(shè)置比水平滲透系數(shù)小一個(gè)數(shù)量級(jí)。根據(jù)經(jīng)驗(yàn)值，有效孔隙度設(shè)置為0.3。由于水動(dòng)力彌散尺度效應(yīng)的存在，難以通過(guò)野外或室內(nèi)彌散試驗(yàn)獲得真實(shí)的彌散度。因此，模擬參照經(jīng)驗(yàn)值，按照保守的評(píng)價(jià)原則，取縱向彌散度值為15 m，橫向彌散度值為1.5 m，垂向彌散度為0.15 m。

1.3.4 源匯項(xiàng)

模擬區(qū)內(nèi)地下水主要接收降水補(bǔ)給，側(cè)向流入補(bǔ)給幾乎可以忽略。根據(jù)包氣帶巖性和地下水埋深特征，結(jié)合經(jīng)驗(yàn)值，確定大氣降水入滲系數(shù)為0.11。由于地下水潛水水位埋深較淺，一般為0.5～2.0 m，蒸騰作用較為強(qiáng)烈，區(qū)內(nèi)地下水主要通過(guò)蒸發(fā)排泄，少部分地下水通過(guò)側(cè)向流出進(jìn)行排泄。此外，根據(jù)地下水潛水監(jiān)測(cè)流場(chǎng)形貌判斷，研究區(qū)內(nèi)存在局部小規(guī)模地下水降落漏斗，推測(cè)該處可能存在少量地下水人工開采。

1.3.5 源強(qiáng)概化

根據(jù)該化工產(chǎn)業(yè)園區(qū)污水處理廠設(shè)計(jì)進(jìn)水水質(zhì)指標(biāo)限值以及各污染因子進(jìn)水濃度和標(biāo)準(zhǔn)限值所計(jì)算的標(biāo)準(zhǔn)指數(shù)，發(fā)現(xiàn)石油類標(biāo)準(zhǔn)指數(shù)最高，對(duì)地下水所構(gòu)成的潛在威脅最大。因此，將研究中污水處理設(shè)施滲漏情境下對(duì)地下水可能造成影響的預(yù)測(cè)因子確定為石油類。

根據(jù)項(xiàng)目工程設(shè)計(jì)，污水處理各構(gòu)筑物的池體均為鋼筋混凝土。池體的泄漏量參考《給水排水構(gòu)筑物工程施工及驗(yàn)收規(guī)范》(GB 50141—2008)中關(guān)于滿水試驗(yàn)驗(yàn)收的要求，鋼筋混凝土池體滿水試驗(yàn)驗(yàn)收標(biāo)準(zhǔn)為2.0 L/(m2·d)，假設(shè)項(xiàng)目在非正常狀況下池底由于地面不均勻沉降或地下水對(duì)池體的腐蝕等多種因素影響下，出現(xiàn)防滲層破裂情況，破裂程度引起的地下水滲漏量按照驗(yàn)收標(biāo)準(zhǔn)的10倍計(jì)算[13-17]，即20 L/(m2·d)。

沉淀池池體長(zhǎng)9.7 m、寬8.5 m、高5.0 m，本次采用池體的表面積作為泄漏面積，即264.45 m2，結(jié)合單位面積滲漏量20 L/(m2·d)，可計(jì)算得到非正常狀況下沉淀池池體滲漏量為5.289 m3/d。由工程分析可知石油類污染物的進(jìn)水濃度為20 mg/L，假定污水處理設(shè)施的地下水防控或檢漏周期為120 d，即發(fā)生非正常狀況后120 d發(fā)現(xiàn)并進(jìn)行修復(fù)切斷滲漏源，因此進(jìn)入含水層中石油類污染物滲漏量為12 213.6 g。

1.4 數(shù)學(xué)模型和軟件選取

1.4.1 地下水流模型

針對(duì)本次所建立的非均質(zhì)各向異性三維穩(wěn)定流地下水流模型，選取下列數(shù)學(xué)模型和定解條件進(jìn)行求解：

(1)

式中，h為地下水潛水位；Ω為地下水滲流區(qū)域；Γ1和Γ2分別為模型第一類邊界和第二類邊界；Kxx、Kyy和Kzz分別為x、y和z主方向的滲透系數(shù)；λ為通過(guò)單位水平面積潛水面注入含水層的水量；ω為源匯項(xiàng)，包括降水入滲補(bǔ)給和井抽水量等；h0(x，y，z)為第一類邊界的地下水水頭函數(shù)；q(x，y，z)為第二類邊界單位面積流量函數(shù)。

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

地下水水質(zhì)模型是定量評(píng)價(jià)和研究地下水水質(zhì)問(wèn)題的數(shù)學(xué)模型。建立地下水水質(zhì)模型時(shí)，將溶質(zhì)運(yùn)移過(guò)程中的對(duì)流、彌散和吸附等反應(yīng)項(xiàng)考慮在內(nèi)，用于描述地下水中溶質(zhì)的運(yùn)移狀態(tài)。模型的偏微分方程式及定解條件如下：

(2)

式中，方程左側(cè)前2項(xiàng)為彌散項(xiàng)，之后2項(xiàng)為對(duì)流項(xiàng)，最后一項(xiàng)為由化學(xué)反應(yīng)或吸附解吸所產(chǎn)生的溶質(zhì)增量；Dxx和Dyy分別為x和y兩個(gè)主方向的彌散系數(shù)；μx和μy為x和y方向上的實(shí)際水流速度；c和c0分別為溶質(zhì)濃度和初始濃度；Ω為溶質(zhì)運(yùn)移的區(qū)域；Γ為三類邊界；φ為邊界溶質(zhì)通量；v為滲流速度；gradc為濃度梯度。

1.4.3 軟件選取

研究所選用的GMS軟件綜合了MODFLOW、MODPATH和MT3DMS等已有的地下水模型基礎(chǔ)上研發(fā)的。GMS軟件相較于Visual MODFLOW等地下水模擬軟件最大的特點(diǎn)是提出了概念模型建模的方法，這為GMS與GIS軟件間提供了良好的數(shù)據(jù)接口，使得建模過(guò)程更為靈活簡(jiǎn)便。除此之外，GMS因其強(qiáng)大的前處理和后處理功能以及優(yōu)良的三維可視化功能，目前已經(jīng)成為了國(guó)際上最受歡迎的地下水模擬軟件之一。

1.5 網(wǎng)格剖分和應(yīng)力期劃分

1.5.1 網(wǎng)格剖分

空間網(wǎng)格剖分是建立地下水水流和溶質(zhì)運(yùn)移數(shù)值模型的基礎(chǔ)。考慮到需要利用所建立的模型對(duì)重點(diǎn)區(qū)域污染物遷移過(guò)程進(jìn)行模擬預(yù)測(cè)，因而在相應(yīng)區(qū)域進(jìn)行了網(wǎng)格的加密剖分。東西向剖分網(wǎng)格64個(gè)，南北向剖分網(wǎng)格74個(gè)，垂向上剖分網(wǎng)格5層，共剖分網(wǎng)格23 680個(gè)，其中有效計(jì)算單元15 960個(gè)。具體空間網(wǎng)格剖分如圖2所示。

圖2 數(shù)值模型空間網(wǎng)格剖分示意Fig.2 Schematic of the spatial grid profiles of the numerical model

1.5.2 應(yīng)力期劃分

此次模擬將地下水非穩(wěn)定流模型和溶質(zhì)運(yùn)移模型模擬期劃分為2個(gè)應(yīng)力期：①第1個(gè)應(yīng)力期設(shè)置為120 d，對(duì)應(yīng)污染物入滲階段；②第2個(gè)應(yīng)力期設(shè)置為11 000 d(30年)，對(duì)應(yīng)污染物擴(kuò)散階段，最長(zhǎng)時(shí)間為污水處理廠最長(zhǎng)服務(wù)年限(30年)。

2 結(jié)果與討論

2.1 地下水流模型識(shí)別檢驗(yàn)

2.1.1 地下水位擬合

通過(guò)模型參數(shù)識(shí)別，得出模型計(jì)算地下水流場(chǎng)和實(shí)際監(jiān)測(cè)地下水流場(chǎng)擬合情況如圖3所示。

圖3 計(jì)算水位與實(shí)際監(jiān)測(cè)水位擬合Fig.3 Fitting of calculated water level and actual monitoring water level

從圖3中可以看到，模型計(jì)算地下水流場(chǎng)和實(shí)際監(jiān)測(cè)地下水流場(chǎng)擬合情況較好，二者之間的水位誤差整體上小于±0.1 m，說(shuō)明所建立的地下水流模型能夠較為準(zhǔn)確地反映評(píng)價(jià)區(qū)地下水潛水的真實(shí)流動(dòng)狀態(tài)，所建立的模型可用于評(píng)價(jià)區(qū)地下水中污染物遷移過(guò)程的模擬預(yù)測(cè)。

2.1.2 地下水均衡分析

通過(guò)模型參數(shù)識(shí)別，得出模型驗(yàn)證期地下水水量均衡結(jié)果見表1。從表1中可以看出，模擬區(qū)地下水潛水主要接受大氣降水入滲補(bǔ)給，并主要通過(guò)潛水蒸騰方式進(jìn)行排泄，并有少量地下水通過(guò)側(cè)向流出進(jìn)行排泄。

表1 模擬區(qū)地下水均衡分析Tab.1 Groundwater balance analysis in simulated area

2.2 污水處理設(shè)施滲漏模擬預(yù)測(cè)

2.2.1 地下水流場(chǎng)動(dòng)態(tài)模擬

滲漏事故情景下，污水滲入地下會(huì)構(gòu)成新的地下水潛水補(bǔ)給來(lái)源，除此之外，模擬區(qū)地下水系統(tǒng)的其他源匯項(xiàng)基本保持不變。因此，非正常狀況下污染物在地下水中的遷移預(yù)測(cè)，有必要在已有地下水流模型的源匯項(xiàng)條件和含水層特征的基礎(chǔ)上，增加入滲補(bǔ)給地下水潛水的污水水量，重新進(jìn)行地下水潛水非穩(wěn)定流模型計(jì)算。污水入滲時(shí)間設(shè)置為120 d，計(jì)算得到的污水處理廠周邊地下水流場(chǎng)變化情況如圖4所示。

圖4 非正常狀況下污水處理廠周邊地下水流場(chǎng)變化情況Fig.4 Changes of groundwater flow field around sewage treatment plant under abnormal conditions

從圖4中可以看到，隨著滲漏時(shí)間延長(zhǎng)，在污水處理設(shè)施周邊開始形成地下水水丘，且規(guī)模不斷增大，在120 d時(shí)達(dá)到最大，且在模擬1 000 d時(shí)也未能完全恢復(fù)。地下水水丘的形成使得泄漏點(diǎn)周邊地下水水力梯度增大，這會(huì)加快污染物向周邊的擴(kuò)散速率。

2.2.2 地下水溶質(zhì)運(yùn)移模擬

污水處理設(shè)施泄漏對(duì)地下水污染預(yù)測(cè)結(jié)果如圖5所示。

圖5 污水處理設(shè)施泄漏對(duì)地下水污染預(yù)測(cè)結(jié)果Fig.5 Prediction results of groundwater pollution caused by leakage of sewage treatment facilities

從圖5中可以看出，在模擬期內(nèi)石油類對(duì)潛水含水層將會(huì)造成污染，并存在超標(biāo)現(xiàn)象。石油類在潛水含水層中不同時(shí)間的超標(biāo)范圍及最大運(yùn)移距離見表2。

表2 非正常狀況下污水處理設(shè)施泄漏對(duì)含水層影響范圍Tab.2 List of influence range of aquifer caused by leakage of sewage treatment facilities under abnormal conditions

從表2中可以看到，1 000 d污染物最遠(yuǎn)遷移41 m左右，超標(biāo)范圍為4 879 m2左右；30年(11 000 d)污染物最遠(yuǎn)遷移56 m左右，超標(biāo)范圍為6 482 m2左右。根據(jù)模擬結(jié)果，在污染物發(fā)生泄漏30年后，地下水中污染物運(yùn)移距離有限，污染羽前緣距離下游地表水體還有120 m左右的距離，這主要與研究區(qū)地下水水力梯度較小、同時(shí)含水層滲透系數(shù)較低有關(guān)。模擬結(jié)果表明，污水處理設(shè)施泄漏對(duì)周邊地表水體影響較小，但由于周邊土壤和地下水遭受污染，對(duì)泄漏點(diǎn)周邊敏感目標(biāo)可能構(gòu)成環(huán)境風(fēng)險(xiǎn)，有必要在污水處理設(shè)施下游進(jìn)行地下水水質(zhì)定期監(jiān)測(cè)。

3 結(jié)論

(1)利用GMS構(gòu)建了研究區(qū)三維非穩(wěn)定地下水流和地下水溶質(zhì)運(yùn)移數(shù)值模型，通過(guò)對(duì)污水處理設(shè)施非正常狀況下廢水的泄漏進(jìn)行模擬預(yù)測(cè)可知：池體泄漏條件下，泄漏點(diǎn)周邊會(huì)形成地下水水丘，在120 d時(shí)規(guī)模達(dá)到最大，加速了污染物向周邊的遷移速率；研究區(qū)屬低滲透性地層，預(yù)測(cè)期內(nèi)污染物運(yùn)移距離較小，尚未對(duì)下游地表水體產(chǎn)生影響。

(2)為防控地下水污染，避免對(duì)泄漏點(diǎn)周邊地下水環(huán)境造成影響，結(jié)合預(yù)測(cè)結(jié)果，建議在污水處理設(shè)施下游30～50 m處布設(shè)跟蹤監(jiān)測(cè)井，定期開展地下水監(jiān)測(cè)工作，以便及時(shí)發(fā)現(xiàn)污染并截?cái)辔廴驹础?/p>

(3)此次模擬遵從保守性原則，未考慮污染物的吸附和降解等生物地球化學(xué)過(guò)程，后期可根據(jù)管理需要，通過(guò)補(bǔ)充相關(guān)反應(yīng)動(dòng)力學(xué)和熱力學(xué)參數(shù)，進(jìn)一步提高模型模擬預(yù)測(cè)精度。