范榮洋，高何鳳，狄 瑜，金曉丹，陳何瀟，李 楊

(1.廣西泰能工程咨詢(xún)有限公司，南寧 530023；2.廣西壯族自治區(qū)環(huán)境保護(hù)科學(xué)研究院，南寧 530022)

前 言

隨著國(guó)民經(jīng)濟(jì)的飛速發(fā)展，城市規(guī)模和經(jīng)濟(jì)總量的不斷擴(kuò)大以及人民消費(fèi)水平的提高，危險(xiǎn)廢物的產(chǎn)生量也在不斷的增加，因此相應(yīng)的危險(xiǎn)廢物焚燒、安全填埋、物化、資源化等處理處置工程也應(yīng)著力推進(jìn)配套建設(shè)。其中危險(xiǎn)廢物安全填埋處置是對(duì)現(xiàn)階段諸多資源化成本高的危險(xiǎn)廢物進(jìn)行減量化和無(wú)害化的一種重要方式，也是現(xiàn)階段危險(xiǎn)廢物處置工作一種不可或缺的手段。但危險(xiǎn)廢物安全填埋場(chǎng)在其運(yùn)營(yíng)期和封場(chǎng)期均會(huì)產(chǎn)生大量的滲濾液，滲濾液中含有多種高濃度污染物質(zhì)[1]，一旦發(fā)生滲漏，各種污染物就會(huì)隨滲濾液通過(guò)破裂的防滲膜和包氣帶進(jìn)入地下水含水層，對(duì)地下水環(huán)境造成不同程度的污染[2]。因此，通過(guò)研究滲濾液中污染物在地下水環(huán)境中的運(yùn)移規(guī)律，從而判定危險(xiǎn)廢物填埋處置項(xiàng)目的地下水環(huán)境可行性，對(duì)于評(píng)價(jià)危險(xiǎn)廢物填埋處置項(xiàng)目建設(shè)對(duì)區(qū)域地下水環(huán)境的影響程度有著重要意義[3]。

1 研究區(qū)概況

1.1 擬建項(xiàng)目概況

擬建某危險(xiǎn)廢物填埋處置項(xiàng)目規(guī)劃用地面積297.2畝，設(shè)計(jì)危險(xiǎn)廢物處理規(guī)模為50 540.2t/a，建設(shè)內(nèi)容主要包括接收暫存設(shè)施、穩(wěn)定化/固化處理系統(tǒng)、滲濾液處理站和安全填埋場(chǎng)等生產(chǎn)設(shè)施及配套的輔助生產(chǎn)和辦公管理設(shè)施，安全填埋場(chǎng)總庫(kù)容約102萬(wàn)m3，有效庫(kù)容83.4萬(wàn)m3，設(shè)計(jì)使用年限為15年。滲濾液處理站采用“氣浮+還原反應(yīng)+中和反應(yīng)+絮凝沉淀+MBR+DTRO”工藝，出水水質(zhì)要求達(dá)到《污水綜合排放標(biāo)準(zhǔn)》(GB8978-1996)一級(jí)標(biāo)準(zhǔn)后外排至園區(qū)污水處理廠。

1.2 水文地質(zhì)概況

1.2.1 研究區(qū)地層巖性

研究區(qū)地層主要由第四系人工填土(Qml)，沖洪積成因(Qal+pl)的黏土、殘坡積成因(Qdl+el)的黏土以及中泥盆統(tǒng)納標(biāo)組(D2n)粉砂質(zhì)泥巖組成。

1.2.2 地下水類(lèi)型及富水性

根據(jù)項(xiàng)目水文地質(zhì)勘察報(bào)告，研究區(qū)內(nèi)地下水按其賦存條件、水理性質(zhì)、水動(dòng)力特征等，可分為松散巖類(lèi)含水巖組和碎屑巖類(lèi)含水巖組2種，相應(yīng)的地下水類(lèi)型分為松散巖類(lèi)孔隙水、碎屑巖類(lèi)構(gòu)造裂隙水2種。松散巖類(lèi)孔隙水主要分布于項(xiàng)目場(chǎng)地內(nèi)的丘陵緩坡地帶以及地勢(shì)低洼處的溪溝地帶，巖性以黏性土為主，部分與下伏基巖含水層有水力聯(lián)系，水量小，基本不具統(tǒng)一的地下水位。主要接受大氣降水和地表水滲入補(bǔ)給，除地勢(shì)低洼處的溪溝底外，該層枯季一般不含水，雨季則常具季節(jié)性的含水特性。該層透水性弱，賦水空間有限，水量貧乏，富水性弱。碎屑巖類(lèi)基巖裂隙水分布于整個(gè)項(xiàng)目場(chǎng)地，含水層巖性主要為泥盆系中統(tǒng)納標(biāo)組(D2n)粉砂質(zhì)泥巖，地下水賦存運(yùn)移于巖石的風(fēng)化裂隙、構(gòu)造裂隙中，全～強(qiáng)風(fēng)化層主要以風(fēng)化裂隙發(fā)育為主，其透水性較好，呈面狀發(fā)育，微風(fēng)化層的裂隙較發(fā)育，但多呈閉合狀且多被泥質(zhì)礦物充填，透水性極差，儲(chǔ)水空間有限，地下水水量貧乏，富水性弱。

1.2.3 地下水的補(bǔ)給、徑流、排泄條件

研究區(qū)位于近南北向發(fā)育的丘陵沖溝，東、南、西側(cè)三面環(huán)山，呈近U字型，山體寬厚，隔水層為微風(fēng)化粉砂質(zhì)泥巖，巖體完整，滲透性弱，分水嶺水位高于設(shè)計(jì)堆填標(biāo)高，產(chǎn)生越嶺滲漏的可能性小，北側(cè)開(kāi)口緊鄰板究小溪，形成一個(gè)相對(duì)獨(dú)立的水文地質(zhì)單元。

安全填埋場(chǎng)處于地下水補(bǔ)給、徑流區(qū)，大氣降水是地下水補(bǔ)給主要來(lái)源，降水以地表徑流流失為主，少量補(bǔ)給地下水，庫(kù)區(qū)巖層弱富水性，水量貧乏。地下水流向與地表水流向基本一致，即總體由南向北徑流匯入板究小溪，補(bǔ)給后地下水位的變化不明顯，其徑流具有距離短且就近排泄的特點(diǎn)。

2 數(shù)值模型建立

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

研究區(qū)模擬計(jì)算范圍東、南、西三面以擬建項(xiàng)目紅線(xiàn)外相連封閉的山體為界，北面以板究小溪為界，為一個(gè)獨(dú)立的次級(jí)水文地質(zhì)單元，該模擬范圍僅略大于廠界范圍，本次預(yù)測(cè)以非正常工況下泄漏到地下水的污染物運(yùn)移至板究小溪后，是否會(huì)對(duì)該水體Ⅲ類(lèi)水質(zhì)功能造成影響作為評(píng)價(jià)依據(jù)。

模擬范圍主要防滲層微風(fēng)化粉砂質(zhì)泥巖K=5.85×10-6～7.4×10-5cm/s，為相對(duì)隔水層。項(xiàng)目擬建安全填埋場(chǎng)庫(kù)底水位變幅0.28～0.67m，各含水層在整個(gè)評(píng)價(jià)范圍內(nèi)呈現(xiàn)較大的穩(wěn)定性，且松散巖類(lèi)孔隙水與碎屑巖類(lèi)構(gòu)造裂隙水之間具有十分密切的水力聯(lián)系，因此綜合考慮將研究區(qū)內(nèi)含水層概化為一層，模型底部以微風(fēng)化粉砂質(zhì)泥巖為隔水底板，模型頂部為地表。區(qū)域包氣帶也較薄，假定不考慮包氣帶對(duì)污染物的截留作用，污染物泄漏后直接進(jìn)入含水層系統(tǒng)。另根據(jù)研究區(qū)豐、平、枯水期水位動(dòng)態(tài)監(jiān)測(cè)，年內(nèi)水位變化幅度不大，可以將模型概化為穩(wěn)定的水流系統(tǒng)。研究區(qū)模擬計(jì)算范圍地下水系統(tǒng)三維空間模型圖見(jiàn)圖1。

圖1 模擬計(jì)算區(qū)域地下水系統(tǒng)三維空間模型Fig.1 Three-dimensional spatial model of regional groundwater system simulated and calculated

根據(jù)模擬范圍的確定，分水嶺概化為隔水邊界，平行等水位線(xiàn)邊界概化為流量邊界，上邊界為地表，接受大氣降水及蒸發(fā)，下邊界以微風(fēng)化粉砂質(zhì)泥巖為隔水底板，地表河流概化為定水頭邊界。

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

2.2.1 地下水流數(shù)學(xué)模型

參考《環(huán)境影響評(píng)價(jià)技術(shù)導(dǎo)則 地下水環(huán)境》(HJ610-2016)[4]和肖長(zhǎng)來(lái)編著的《地下水動(dòng)力學(xué)講義》(吉林大學(xué)環(huán)境與資源學(xué)院 2009年10月)，根據(jù)研究區(qū)的地質(zhì)及水文地質(zhì)情況，將區(qū)域地層概化為非均質(zhì)各向同性穩(wěn)定流地下水流系統(tǒng)，其水流數(shù)值模型如下：

①水流模型控制方程

②初始條件

h(x,y,z,t)=h0(x,y,z) (x,y,z)∈Ω,t=0

③邊界條件

第一類(lèi)邊界條件：

h(x,y,z,t)|Γ 1=h(x,y,z,t) (x,y,z)∈Γ1,t≥0

第二類(lèi)邊界條件：

式中：Kxx，Kyy，Kzz：x，y，z方向的滲透系數(shù)，m/d；h：壓力水頭，m；W：源匯項(xiàng)，1/d；H0：給定的初始?jí)毫λ^，m；H1：第一類(lèi)邊界給定的壓力水頭，m；q：在第二類(lèi)邊界條件上給定的通過(guò)邊界的水流量，m3/d；Ω：滲流場(chǎng)；Γ1：第一類(lèi)邊界條件；Γ2：第二類(lèi)邊界條件；n：邊界Γ2的外法線(xiàn)方向。

2.2.2 溶質(zhì)運(yùn)移數(shù)學(xué)模型

參考《環(huán)境影響評(píng)價(jià)技術(shù)導(dǎo)則 地下水環(huán)境》(HJ610-2016)[4]，本次溶質(zhì)運(yùn)移模擬，本著最不利因素出發(fā)，不考慮污染物在地下水體中的反應(yīng)及吸附作用，只考慮對(duì)流彌散作用，其溶質(zhì)運(yùn)移模型如下：

①控制方程

②初始條件

C(x,y,z,t)=C0(x,y,z) (x,y,z)∈Ω,t=0

③定解條件

第一類(lèi)邊界條件：給定濃度邊界

C(x,y,z,t)|Γ 1=c(x,y,z,t)(x,y,z)∈Γ1,t≥0

第二類(lèi)邊界條件：給定彌散通量邊界

式中：R：遲滯系數(shù)，無(wú)量綱；b：介質(zhì)密度，mg/(dm)3；θ：介質(zhì)孔隙度，無(wú)量綱；C：組分濃度，g/L；Dij：水動(dòng)力彌散系數(shù)張量，m2/d；W：水流源匯項(xiàng)，1/d；Cs：組分的濃度，mg/L；C0：已知濃度分布，mg/L；f：邊界Γ2上的已知彌散通量函數(shù)；Ω：滲流場(chǎng)；Γ1：第一類(lèi)邊界條件；Γ2：第二類(lèi)邊界條件。

2.2.3 模擬軟件的選取

本次評(píng)價(jià)工作采用加拿大Waterloo Hydrogeologic Inc.在美國(guó)地質(zhì)調(diào)查局MODFLOW軟件(1984年)的基礎(chǔ)上應(yīng)用可視化技術(shù)開(kāi)發(fā)研制的Visual MODFLOW軟件。該軟件是目前世界上應(yīng)用最廣泛的三維地下水流和溶質(zhì)運(yùn)移模擬的標(biāo)準(zhǔn)可視化專(zhuān)業(yè)軟件系統(tǒng)。該軟件主要由下列軟件包組成：MODFLOW—水流模擬；Zone Budget—水均衡分析；MODPATH—流線(xiàn)示蹤分析；MT3D—溶質(zhì)運(yùn)移模擬；WinPEST—參數(shù)自動(dòng)識(shí)別；3D-Explorer—模擬結(jié)果的三維顯示。

2.2.4 水文地質(zhì)參數(shù)選取

地下室數(shù)值模擬中水文地質(zhì)參數(shù)非常重要的，其合理性與正確性將直接影響地下水模型的準(zhǔn)確性和可信度。本次模擬水文地質(zhì)參數(shù)是根據(jù)水文地質(zhì)勘察成果確定，通過(guò)試坑滲水試驗(yàn)、鉆孔注水試驗(yàn)、抽水試驗(yàn)、壓水試驗(yàn)及室內(nèi)土樣滲透性試驗(yàn)結(jié)果，得出含水層的垂直滲透系數(shù)和水平滲透系數(shù)、給水度、入滲系數(shù)等主要輸入?yún)?shù)。彌散度的確定相對(duì)比較困難，雖然水文地質(zhì)勘查進(jìn)行了野外彌散試驗(yàn)，但野外彌散試驗(yàn)所求出的彌散度遠(yuǎn)遠(yuǎn)大于在實(shí)驗(yàn)室所測(cè)出的值，相差可達(dá)4～5個(gè)數(shù)量級(jí)，即使是同一含水層，溶質(zhì)運(yùn)移距離越大，所計(jì)算出的彌散度也越大。故即使是進(jìn)行野外或室內(nèi)彌散試驗(yàn)也難以獲得準(zhǔn)確的彌散度值。因此，本次模擬參考李國(guó)敏等人的研究成果，縱向彌散度取10m[5]。根據(jù)經(jīng)驗(yàn)，水平橫向彌散度取值應(yīng)比縱向彌散度小一個(gè)數(shù)量級(jí)，垂直橫向彌散度的取值應(yīng)比縱向彌散度小兩個(gè)數(shù)量級(jí)。

3 地下水流數(shù)值模擬

根據(jù)建立的水文地質(zhì)概念模型及數(shù)學(xué)模型，運(yùn)用基于有限差分法的Visual Modflow軟件建立了項(xiàng)目區(qū)的數(shù)值模型，然后經(jīng)過(guò)參數(shù)識(shí)別及校正，對(duì)項(xiàng)目區(qū)的地下水流模型進(jìn)行評(píng)價(jià)分析，為溶質(zhì)運(yùn)移模擬做基礎(chǔ)。

3.1 模型網(wǎng)格剖分

根據(jù)研究區(qū)的基本情況，由于項(xiàng)目所在水文地質(zhì)單元僅略大于廠界范圍，模擬范圍較小，故沿x、y方向分別以10m細(xì)間距等距剖分，整個(gè)模擬范圍均為加密網(wǎng)格布置，共劃分單元格12 100個(gè)，整個(gè)評(píng)價(jià)區(qū)劃分為二層，模型網(wǎng)格剖分示意圖見(jiàn)圖2。

圖2 模型網(wǎng)格剖分示意圖Fig.2 Schematic diagram of model mesh generation

3.2 源匯項(xiàng)處理

根據(jù)水文地質(zhì)報(bào)告，研究區(qū)含水層地下水類(lèi)型為松散巖類(lèi)孔隙水、碎屑巖類(lèi)構(gòu)造裂隙水2種，主要為基巖裂隙水，水量貧乏,降雨為其主要的補(bǔ)給來(lái)源，側(cè)向徑流為其主要的排泄方式。根據(jù)邊界條件的概化，本次地下水模擬各源匯項(xiàng)均通過(guò)Visual Modflow各子程序包賦到模型里。其中：北側(cè)為流量邊界，采用GeneralHead(GHB)子程序包賦值；降雨采用Recharge子程序包賦值；蒸發(fā)采用Evapotranspiration賦值；開(kāi)采井采用PumpingWell子程序包給模型賦值，模型中不賦值條件即默認(rèn)為隔水邊界。

3.3 參數(shù)識(shí)別與模型校驗(yàn)

采用地下水流場(chǎng)作為識(shí)別方法，將相關(guān)水文地質(zhì)參數(shù)的值及模型的邊界條件等輸入模型，運(yùn)行模型得到計(jì)算流場(chǎng)，與實(shí)際測(cè)量的地下水流場(chǎng)進(jìn)行對(duì)比，不斷調(diào)整相關(guān)水文地質(zhì)參數(shù)，使計(jì)算流場(chǎng)與實(shí)測(cè)流場(chǎng)相吻合。識(shí)別方法采用試估校正法，模型中使用Visual Modflow 程序中自帶的PEST子程序包進(jìn)行。最終得到與實(shí)測(cè)流場(chǎng)較為擬合的計(jì)算流場(chǎng)，見(jiàn)圖3。

由上圖，模擬計(jì)算流場(chǎng)不管從水頭值還是從流場(chǎng)的分布情況來(lái)看，模擬計(jì)算流場(chǎng)均與實(shí)測(cè)流場(chǎng)擬合較好，模型可以反映實(shí)際的地下水情況。同時(shí)，模擬計(jì)算流場(chǎng)跟實(shí)際的補(bǔ)給關(guān)系也一致，符合研究區(qū)的水文地質(zhì)條件，基本反映了研究區(qū)地下水系統(tǒng)的水動(dòng)力特征，可用于本次溶質(zhì)運(yùn)移模擬。

圖3 流場(chǎng)擬合圖Fig.3 Flow field fitting diagram

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

4.1 預(yù)測(cè)情景設(shè)置

擬建危險(xiǎn)廢物安全填埋場(chǎng)在其運(yùn)營(yíng)期和封場(chǎng)期均會(huì)產(chǎn)生大量的滲濾液，滲濾液中含有多種高濃度污染物質(zhì)，一旦進(jìn)入地下水，可能會(huì)對(duì)地下水環(huán)境造成污染。為提前預(yù)知滲漏污染造成的影響范圍、污染程度及其運(yùn)移規(guī)律，本次模擬選取安全填埋場(chǎng)防滲膜破裂滲濾液泄漏的事故情景進(jìn)行預(yù)測(cè)。泄漏區(qū)域選取位置為安全填埋場(chǎng)壩體的最下游，在這種最不利情景下，泄漏的污染物距離板究小溪更近，可以較快的到達(dá)受納水體且污染物的濃度削減較少，從而使本次預(yù)測(cè)情景具有代表性?？紤]到防滲膜可能存在的接縫疏忽或鋪設(shè)不到位等情況，將可能發(fā)生滲漏的面積設(shè)定為安全填埋場(chǎng)底部面積的5%。滲漏量根據(jù)公式Q=A×K×T進(jìn)行計(jì)算，其中：A為滲漏面積，m2；K為垂向滲透系數(shù)，m/d；T為滲漏時(shí)間，d，由于地下水長(zhǎng)期監(jiān)測(cè)點(diǎn)的監(jiān)測(cè)頻率為60 d，整體修復(fù)時(shí)間定為30 d，因此將滲漏時(shí)間定為90 d。

本次模擬時(shí)段設(shè)置為從2020年1月開(kāi)始，連續(xù)模擬預(yù)測(cè)服務(wù)期15年，按5 500 d計(jì)，以10 d或30 d等不同的時(shí)間步長(zhǎng)來(lái)預(yù)測(cè)。預(yù)測(cè)綜合考慮當(dāng)?shù)厮占ｋU(xiǎn)廢物的特征性，同時(shí)結(jié)合等標(biāo)污染負(fù)荷及毒性大小選取有代表性的Cr6+較作為預(yù)測(cè)因子[6]，利用驗(yàn)證后的地下水水流數(shù)值模型預(yù)測(cè)泄漏點(diǎn)污染物隨地下水流擴(kuò)散所造成的影響。

4.2 模擬條件概化

本次模擬將上述情景的污染源設(shè)置為點(diǎn)源濃度邊界，污染源位置按實(shí)際位置概化。本次模擬只考慮污染物在地下水系統(tǒng)中的對(duì)流、彌散作用，不考慮地層的吸附、解吸作用，不考慮化學(xué)反應(yīng)及生物降解等作用。同時(shí)，包氣帶較薄，不考慮包氣帶的阻滯作用，假設(shè)污染物泄漏后直接進(jìn)入含水層中進(jìn)行溶質(zhì)運(yùn)移[7]。

4.3 預(yù)測(cè)結(jié)果及評(píng)價(jià)

將污染源賦給Visual Modflow溶質(zhì)運(yùn)移模型MT3D中，將安全填埋場(chǎng)泄漏設(shè)置為點(diǎn)源濃度邊界條件，在此基礎(chǔ)上運(yùn)行水流及溶質(zhì)模型，得到Cr6+的運(yùn)移擴(kuò)散結(jié)果。

根據(jù)預(yù)測(cè)，污染物Cr6+停止泄漏后5 500天可滿(mǎn)足《地下水質(zhì)量標(biāo)準(zhǔn)》(GB/T14848-2017)中Ⅲ類(lèi)標(biāo)準(zhǔn)鉻(六價(jià))0.05mg/L的限值要求，Cr6+泄漏90天、停止泄漏后1 000 d、2 190 d、3 650 d及5 450 d的預(yù)測(cè)結(jié)果見(jiàn)圖4～圖8，其中圖右側(cè)的剖面圖為縱剖圖，這是由于本項(xiàng)目區(qū)域地下水流向是自南向北，縱剖向?yàn)槲廴疚镞\(yùn)移流動(dòng)方向；各圖中的最外側(cè)邊界線(xiàn)以《地下水質(zhì)量標(biāo)準(zhǔn)》(GB/T14848-2017)中Ⅲ類(lèi)鉻(六價(jià))的濃度0.05mg/L為控制線(xiàn)，預(yù)測(cè)各時(shí)段Cr6+的最大運(yùn)移距離及擴(kuò)散面積分別見(jiàn)下表。

圖4 安全填埋場(chǎng)滲漏90d時(shí)Cr6+運(yùn)移平面及剖面濃度分布Fig.4 Concentration distribution of Cr6+ migration plane and profile at 90 days of leakage in safe landfills

表 各時(shí)段Cr6+運(yùn)移情況Tab. Cr6+ migration in each period

圖5 安全填埋場(chǎng)滲漏1 000d時(shí)Cr6+運(yùn)移平面及剖面濃度分布Fig.5 Concentration distribution of Cr6+ migration plane and profile at 1 000 days of leakage in safe landfills

圖6 安全填埋場(chǎng)滲漏2 190d時(shí)Cr6+運(yùn)移平面及剖面濃度分布Fig.6 Concentration distribution of Cr6+ migration plane and profile at 2 190 days of leakage in safe landfills

圖7 安全填埋場(chǎng)滲漏3 650d時(shí)Cr6+運(yùn)移平面及剖面濃度分布Fig.7 Concentration distribution of Cr6+ migration plane and profile at 3 650 days of leakage in safe landfills

圖8 安全填埋場(chǎng)滲漏5 450d時(shí)Cr6+運(yùn)移平面及剖面濃度分布Fig.8 Concentration distribution of Cr6+ migration plane and profile at 5 450 days of leakage in safe landfills

Cr6+滲漏污染地下水影響評(píng)價(jià)：從上表及圖4～圖8可知，在安全填埋場(chǎng)底部防滲膜破損的情況下，Cr6+會(huì)滲入地下污染地下水，泄漏90 d時(shí)，污染源前鋒距泄漏點(diǎn)最大距離為47.9m，超標(biāo)面積為3 748m2，最大超標(biāo)倍數(shù)為4倍；在停止泄漏后的溶質(zhì)運(yùn)移過(guò)程中，最大超標(biāo)區(qū)域面積為8 101m2。污染物在地下水的擴(kuò)散作用下，濃度逐漸降低，經(jīng)過(guò)5 500 d后滿(mǎn)足《地下水質(zhì)量標(biāo)準(zhǔn)》(GB/T14848-2017)Ⅲ類(lèi)要求，此時(shí)，Cr6+污染暈最遠(yuǎn)運(yùn)行141.7m。

對(duì)下游板究小溪Ⅲ類(lèi)水體水質(zhì)功能的影響：根據(jù)預(yù)測(cè)結(jié)果，污染物Cr6+運(yùn)移5 500 d最遠(yuǎn)運(yùn)行至安全填埋場(chǎng)下游141.7m后即滿(mǎn)足《地下水質(zhì)量標(biāo)準(zhǔn)》(GB/T14848-2017)Ⅲ類(lèi)標(biāo)準(zhǔn)鉻(六價(jià))0.05mg/L的限值要求，而此時(shí)污染暈還在擬建項(xiàng)目紅線(xiàn)范圍內(nèi)，基本不會(huì)對(duì)下游板究小溪水體水質(zhì)功能造成不利影響[8]。

5 結(jié) 語(yǔ)

本文以某危險(xiǎn)廢物填埋處置項(xiàng)目為例，采用Visual MODFLOW軟件進(jìn)行數(shù)值模擬的方法，預(yù)測(cè)項(xiàng)目建成后，在安全填埋場(chǎng)防滲膜破裂導(dǎo)致滲濾液泄漏的事故情景下，污染物隨時(shí)間的運(yùn)移過(guò)程和濃度分布變化[9]。通過(guò)模擬結(jié)果可以看出，發(fā)生事故后，污染物Cr6+由于彌散作用緩慢在地下水流向上進(jìn)行橫向和縱向擴(kuò)散，至填埋場(chǎng)服務(wù)期滿(mǎn)時(shí)，研究區(qū)地下水環(huán)境Cr6+濃度已能滿(mǎn)足《地下水質(zhì)量標(biāo)準(zhǔn)》(GB/T14848-2017)Ⅲ類(lèi)要求，而滲濾液在此前形成的污染暈最遠(yuǎn)運(yùn)移位置仍在項(xiàng)目紅線(xiàn)范圍內(nèi)，因此可以得出擬建危險(xiǎn)廢物填埋處置項(xiàng)目對(duì)區(qū)域地下水環(huán)境影響較小的結(jié)論。

基于Visual MODFLOW軟件進(jìn)行地下水?dāng)?shù)值模擬的方法，可以客觀的分析污染物在地下水中的遷移特征，評(píng)價(jià)其對(duì)地下水環(huán)境及敏感目標(biāo)的影響程度，從而為污染類(lèi)項(xiàng)目是否具有地下水環(huán)境可行性提供決策依據(jù)，并為配套相應(yīng)污染防治措施提供技術(shù)支持。