張子涵

(山東農(nóng)業(yè)大學(xué)，山東 泰安 271018)

水是最重要的資源，人類的一切活動(dòng)都離不開(kāi)水的支撐。雖然目前我國(guó)已經(jīng)成為全球第二大經(jīng)濟(jì)體，但這是在“高能耗、高污染、高排放”基礎(chǔ)上實(shí)現(xiàn)的。據(jù)統(tǒng)計(jì)，我國(guó)超過(guò)90%的城市水域污染嚴(yán)重，每年因水污染造成的經(jīng)濟(jì)損失高達(dá)2400億元。因此在保證各行業(yè)正常運(yùn)行的前提下，必須做好水源地劃分和保護(hù)工作，而目前應(yīng)用較為廣泛的是MODFLOW模擬技術(shù)。

1 模擬區(qū)域概況

本文所選擇的模擬區(qū)域位于某市機(jī)場(chǎng)附近，劃定模擬區(qū)域尺寸為：長(zhǎng)×寬=2km×2km，模擬層數(shù)為3層，從上至下依次為：淺砂層(砂礫含水層)、黏土半透水層、深砂層(砂礫含水層)(模型示意圖如圖1所示)。模擬場(chǎng)的污染源為飛機(jī)燃料補(bǔ)充區(qū)，保護(hù)區(qū)為城市供水水源地，2者中間分布有一個(gè)半透水地層。

圖1 模擬區(qū)域水文地質(zhì)模型(單位：m)

2 地下水源保護(hù)區(qū)劃分方法

目前，主要依據(jù)水文地質(zhì)條件來(lái)劃定地下水源保護(hù)區(qū)范圍，需要掌握的資料包括含水層特性、地下水通道、污染物的遷移等。根據(jù)國(guó)內(nèi)外的經(jīng)驗(yàn)和標(biāo)準(zhǔn)，將地下水保護(hù)區(qū)分為3級(jí)：一級(jí)保護(hù)區(qū)、二級(jí)保護(hù)區(qū)和準(zhǔn)保護(hù)區(qū)[1]。

(1)一級(jí)保護(hù)區(qū)。地下水保護(hù)區(qū)的主要?jiǎng)澏ㄒ罁?jù)時(shí)地下水運(yùn)動(dòng)規(guī)律，當(dāng)然這些規(guī)律時(shí)通過(guò)當(dāng)?shù)厮牡刭|(zhì)條件推斷出來(lái)的。相關(guān)規(guī)范規(guī)定：按照地面的細(xì)菌和污染物滲透到地下水源所用時(shí)間在60d以內(nèi)的區(qū)域統(tǒng)一劃定為地下水源一級(jí)保護(hù)區(qū)。該區(qū)域是保證地下水安全的最敏感地帶，如果該區(qū)域被污染，則容易造成大范圍地區(qū)缺水問(wèn)題，該區(qū)域能暫時(shí)保證用水安全。

(2)二級(jí)保護(hù)區(qū)。在一級(jí)保護(hù)區(qū)以外的一定范圍內(nèi)，將其稱為二級(jí)保護(hù)區(qū)，一般地下水水源的補(bǔ)給區(qū)均分布于該區(qū)域，所以其重要性也很大。污染源通過(guò)二級(jí)保護(hù)區(qū)進(jìn)入的一級(jí)保護(hù)區(qū)，會(huì)在中途的被土壤及水分子的吸附作用稀釋很大一部分，因?yàn)閷?duì)保護(hù)水源地安全意義重大。相關(guān)規(guī)范規(guī)定：二級(jí)保護(hù)區(qū)地下水運(yùn)移時(shí)間范圍一般設(shè)定為10a，在這個(gè)范圍內(nèi)的區(qū)域統(tǒng)稱為二級(jí)保護(hù)區(qū)[2]。

(3)準(zhǔn)保護(hù)區(qū)。在二級(jí)保護(hù)區(qū)以外的一定范圍內(nèi)設(shè)置準(zhǔn)保護(hù)區(qū)，其主要作用是防止污染物過(guò)快進(jìn)入二級(jí)保護(hù)區(qū)，主要保護(hù)對(duì)象為集水區(qū)，范圍內(nèi)采空層地下水流入采空井群。相關(guān)規(guī)范規(guī)定：將污染物在地下水中移動(dòng)25a的范圍作為準(zhǔn)保護(hù)區(qū)域范圍[3]。

3 MODFLOW模擬技術(shù)的應(yīng)用分析

MODFLOW是美國(guó)地質(zhì)調(diào)查局在20世紀(jì)80年代開(kāi)發(fā)的一種用于地下水流數(shù)值模擬的三維有限元分析軟件，采用模塊化原理，可模擬的地下水流場(chǎng)和不同時(shí)間水力截流區(qū)來(lái)劃分水源保護(hù)區(qū)。在此規(guī)定：污染物向水源中心移動(dòng)100d和1000d的捕獲區(qū)分別設(shè)為一級(jí)保護(hù)區(qū)和二級(jí)保護(hù)區(qū)，二級(jí)保護(hù)區(qū)外的地下水補(bǔ)給區(qū)設(shè)為準(zhǔn)保護(hù)區(qū)。

3.1 模擬區(qū)域邊界及水文地質(zhì)參數(shù)的確定

3.1.1模擬區(qū)域邊界確定

(1)上部含水層。模擬區(qū)的上部含水層為潛水含水層，均值各向同性。南部有河流，設(shè)計(jì)南北邊界均為給定水頭邊界，其中南部邊界給定水頭6m，北部邊界給定水頭5.2m。上部含水層的東西部均為流量邊界[4]。

(2)下部含水層。模擬區(qū)的下部含水層為承壓含水層；均值各向同性；南北為給定水頭邊界，其中南部邊界給定水頭為18m，北部邊界給定水頭為16.5m；東西部方向?yàn)榻o定流量邊界抽水井P1和P2，取水量分別為400、550m3/d。

3.1.2模擬區(qū)域水文地質(zhì)參數(shù)確定

將模擬區(qū)剖分為900個(gè)有效單元格，每個(gè)單元格尺寸為30m×30m,垂向上為3層。根據(jù)水文地質(zhì)概念模型，設(shè)計(jì)第一層為砂巖，第二層為泥巖，第三層為砂巖。其中水文地質(zhì)初始參數(shù)取值見(jiàn)表1，最終確定的溶質(zhì)運(yùn)移模型參數(shù)見(jiàn)表2[5]。

表1 區(qū)域水文地質(zhì)參數(shù)取值表 單位：m/d

表2 溶質(zhì)運(yùn)移模型參數(shù)表

3.1.3建立溶質(zhì)運(yùn)移模型

(1)調(diào)用MODPATH，在燃料堆放地中投入示蹤劑。進(jìn)行初始化設(shè)置，定義溶質(zhì)運(yùn)移的單位。

(2)對(duì)溶質(zhì)運(yùn)移基本程序進(jìn)行設(shè)置，包括：溶質(zhì)運(yùn)移模擬描述、運(yùn)移程序包選擇、定義組分、溶質(zhì)運(yùn)移問(wèn)題描述、模擬期、輸出控制、初始濃度、孔隙度等的定義或輸出[6]。

(3)在MAP模塊中建立溶質(zhì)運(yùn)移概念模型，包括：縱向彌散度、孔隙度、補(bǔ)給濃度賦值。其中體積密度為1700kg/m3，補(bǔ)給濃度為500.0mg/L。

(4)在3D模塊中對(duì)于對(duì)流彌散程序包進(jìn)行設(shè)置。定義橫向彌散度與縱向彌散度的比值，垂向彌散度與縱向彌散度的比值。

(5)將概念模型導(dǎo)入MT3DMS中，模擬運(yùn)行與結(jié)果查看。

3.2 模型運(yùn)行及結(jié)果分析

3.2.1燃料堆放區(qū)污染物運(yùn)移情況分析

模擬區(qū)域附近機(jī)場(chǎng)的燃料堆放區(qū)是本項(xiàng)目模擬的主要污染源，水動(dòng)力場(chǎng)是帶動(dòng)污染物遷移的主要?jiǎng)恿?。污染物位于模擬區(qū)的正北方，而水動(dòng)力場(chǎng)時(shí)由北向南分布，因此污染物的運(yùn)動(dòng)方向也是逐步向南擴(kuò)散。在此選取100、1000、7300d三個(gè)時(shí)間段對(duì)污染物擴(kuò)散區(qū)域進(jìn)行模擬[7]。

3.2.2模擬結(jié)果分析

圖2—3展示了模型運(yùn)行100、1000、7300d三個(gè)不同時(shí)段下地下水中污染物的遷移擴(kuò)散情況。從圖中可知：①污染物的遷移方向逐步由北向東南部擴(kuò)散，污染源逐步擴(kuò)大。其中1000d的范圍是100d的1.29倍左右，7300d的范圍是100d的3.26倍左右，見(jiàn)表3；②在地下水稀釋作用下，隨著污染物擴(kuò)散范圍越廣，最高濃度逐步降低[8]。

表3 機(jī)場(chǎng)燃料堆放區(qū)污染物運(yùn)移結(jié)果

圖2 第100d(上)和1000d(下)污染物運(yùn)移情況

圖3 第7300d污染物運(yùn)移情況

3.2.3模擬區(qū)水源地保護(hù)區(qū)劃分

在此參考《飲用水水源保護(hù)區(qū)劃分技術(shù)規(guī)劃》，并綜合我國(guó)多地區(qū)相似情況的均值，確定水源地一、二級(jí)保護(hù)區(qū)的降解時(shí)間分別為100、1000d[9]。

本項(xiàng)目利用“示蹤粒子跟蹤法”確定了該模擬區(qū)域地下水運(yùn)動(dòng)遷移情景，并以此確定相關(guān)地下水保護(hù)區(qū)域，如圖4所示[10]，本項(xiàng)目劃定的各保護(hù)區(qū)劃分面積見(jiàn)表4。

表4 水源地保護(hù)區(qū)劃分結(jié)果

圖4 示蹤粒子100、1000d運(yùn)移軌跡圖

4 結(jié)語(yǔ)

水源地保護(hù)區(qū)的劃分對(duì)于保證當(dāng)?shù)氐叵滤|(zhì)安全有著重要意義，MODFLOW模擬技術(shù)可以在宏觀上給出污染源在一定時(shí)間內(nèi)的擴(kuò)散范圍，作為理論參考值。MODFLOW模擬技術(shù)的應(yīng)用步驟文章已經(jīng)詳細(xì)闡述，其中建立模型是首要步驟，應(yīng)當(dāng)緊密結(jié)合當(dāng)?shù)貙?shí)際情況。但在此需要特別說(shuō)明：由于當(dāng)前模擬技術(shù)是根據(jù)事先給定的限定條件進(jìn)行演示，因此缺少水文地質(zhì)條件在實(shí)際情況中復(fù)雜性和動(dòng)態(tài)性，并不能完全將模擬結(jié)果作為最終參考對(duì)象。本項(xiàng)目區(qū)模擬得到了各級(jí)保護(hù)區(qū)面積，在實(shí)際中還需要修定，但遵循的原則是寧可范圍擴(kuò)大，也要避免縮小。