劉 莉，邢立亭

（1.濟(jì)南天下第一泉風(fēng)景區(qū)，山東 濟(jì)南 250000；2.濟(jì)南大學(xué) 水利與環(huán)境學(xué)院，山東 濟(jì)南 250022）

1 引言

地下水作為人類生產(chǎn)和生活重要的水源，其水質(zhì)好壞直接關(guān)系到人民身體健康，一旦被污染很難阻止或逆轉(zhuǎn)［1］。巖溶水作為地下水的一個分支，其形成是地表水及地下水對可溶巖綜合作用的結(jié)果，巖溶的發(fā)育程度決定和控制區(qū)域水文地質(zhì)條件的轉(zhuǎn)化，又因它和周圍巖石性質(zhì)、地質(zhì)構(gòu)造、地下水運(yùn)動、水化學(xué)等因素有密切關(guān)系，且諸因素相互影響，導(dǎo)致區(qū)域水文地質(zhì)條件復(fù)雜化［2］。

項(xiàng)目區(qū)位于淄博南部巖溶區(qū)，大武水源地西側(cè)。大武水源地作為我國北方罕見的特大型巖溶-裂隙地下水水源地，多年來大規(guī)模開采，造成區(qū)域地下水位持續(xù)下降［3］。同時，眾多化工企業(yè)直接建設(shè)于水源地補(bǔ)給區(qū)之上，由于防滲不規(guī)范、排污的不合格、管理不到位等多重因素的影響，區(qū)域地下水中污染超標(biāo)嚴(yán)重［4-6］，受其影響，北方最大巖溶水源地報廢。

本項(xiàng)目為新建化工項(xiàng)目，在生產(chǎn)環(huán)節(jié)和最終外排均存在污廢水。由于地下水具有相通性及很強(qiáng)的隱蔽性，地下水的污染很不容易及時發(fā)現(xiàn)，一旦發(fā)現(xiàn)，其后果也難以消除［8］。因此對項(xiàng)目開展地下水環(huán)境影響評價，論證建設(shè)選址的合理性和可行性，為保護(hù)和恢復(fù)大武水源地、保障用水安全具有十分重要的指導(dǎo)作用［9］。

圖1 項(xiàng)目區(qū)位置及地質(zhì)簡圖（據(jù)徐道余［7］等）

2 水文地質(zhì)條件

項(xiàng)目區(qū)位于淄博市張店區(qū)，氣候?qū)倥瘻貛О霛駶櫦撅L(fēng)區(qū)大陸性氣候，四季分明，夏季炎熱多雨，冬季寒冷干燥，常年風(fēng)向以南南西和南西為主。夏季多東南風(fēng)，冬季多西北風(fēng)，多年平均氣溫12.9℃。多年平均降水量為588.7mm，降水量年內(nèi)分配極不平均，多集中在6～9 月份，占全年降水量的78.7%；多年平均陸地水面蒸發(fā)量為2077.3mm。地貌類型屬于低山丘陵，地形變化較大，整體地形南高北低，南北高差大于10m。

2.1 地下水賦存條件

區(qū)域主要分布三大含水巖組，分別為第四系松散巖類孔隙含水巖組、奧陶系碳酸鹽類巖溶～裂隙含水巖組和石炭二疊系碳酸鹽類裂隙含水巖組。

松散巖類孔隙含水巖組包括現(xiàn)代河床沖積孔隙含水巖組和上更新統(tǒng)沖積～洪積孔隙含水巖組，分布于項(xiàng)目區(qū)的東側(cè)淄河河漫灘及其兩側(cè)、淄河沖洪積扇首部和北部平原區(qū)，巖性為砂卵礫石層、粗砂卵礫石層等，富水性良好。

奧陶系碳酸鹽類巖溶～裂隙含水巖組在區(qū)內(nèi)分布廣泛，除低山丘陵地帶巖性裸露外，均被第四系松散沉積物所覆蓋。含水層巖性為中奧陶系第三段至第六段的含泥質(zhì)、白云質(zhì)泥灰?guī)r、角礫狀泥灰?guī)r及厚層狀青灰色豹皮狀灰?guī)r組成?；?guī)r裂隙巖溶極為發(fā)育，其發(fā)育深度在60 ～300m 之間。以南仇-安里-大武一帶富水性最強(qiáng)，單井涌水量大于6000m3/d；丘陵區(qū)，地勢高，水位埋深大，富水性最差，單井涌水量小于1000m3/d，其它地區(qū)單井涌水量在1000 ～5000m3/d 之間。

石炭二疊系碳酸鹽類裂隙含水巖組主要分布于湖田向斜兩翼和金嶺村南側(cè)，該含水層地下水位埋藏較淺，一般小于 5.0m，湖田向斜兩側(cè)富水性一般，單井涌水量在100 ～1000m3/d 之間，金嶺村南側(cè)富水性差，單井涌水量小于100m3/d。

圖2 區(qū)域水文地質(zhì)圖

2.2 地下水補(bǔ)給、徑流、排泄

松散巖類孔隙水主要接受大氣降水、南部山區(qū)地下水徑流、淄河滲漏及灰?guī)r地下水通過第四系“天窗”等多種方式補(bǔ)給，排泄方式為蒸發(fā)、泉和溢出帶形式的地表徑流為主，地下水總體由南向北流動。近年來由于自然因素的變化及人類活動的影響（如淄河斷流等），地下水水位下降明顯，人工開采成為主要排泄方式。

碳酸鹽巖類巖溶～裂隙水補(bǔ)給來源主要為淄河河谷兩側(cè)至東、西地表分水嶺地區(qū)灰?guī)r地下水匯集于淄河斷裂帶后的徑流補(bǔ)給、大氣降水的入滲補(bǔ)給等，排泄方式主要為人工開采。

3 地下水環(huán)境影響評價方法

地下水環(huán)境影響評價的方法以解析法和數(shù)值模擬兩種方法為主。解析法主要用于水文地質(zhì)條件相對簡單，污染較輕，精度要求不高的情況；數(shù)值模擬法以可視化分析研究地下水流和溶質(zhì)運(yùn)移問題，其方法直觀、可視、有效，越來越受到領(lǐng)域內(nèi)學(xué)者的重視以及廣泛應(yīng)用［10］。

本次采用數(shù)值模擬法對地下水環(huán)境影響進(jìn)行評價分析，運(yùn)用Visual MODFLOW 軟件，該系統(tǒng)無縫集成了MODFLOW-96、WinPEST、MT3D99、MODPATH、MT3D 等軟件［11］，建立了合理的Windows 菜單界面，具有可視化功能強(qiáng)大、求解方法簡單、適用范圍廣泛、數(shù)值模擬能力出色且三維建模簡單等特點(diǎn)［12］。

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

結(jié)合區(qū)域地質(zhì)水文地質(zhì)條件，確定模型范圍為整個湖田水文地質(zhì)單元，具體為：東邊界為金嶺斷裂，西邊界為炒米莊斷層，北邊界為湖田向斜軸，南邊界為地表水與地下水分水嶺，總面積約70km2。

（1）含水層結(jié)構(gòu)概化。

區(qū)域地下水流以水平運(yùn)動為主，地下水流速在x，y，z 三個方向都有分量，概化為三維流；各個水文地質(zhì)參數(shù)隨空間變化，體現(xiàn)了非均質(zhì)性；整個地下水系統(tǒng)的輸入和輸出量隨時間和空間變化，地下水流為非穩(wěn)定流。將研究區(qū)含水層概化為非均質(zhì)、各向同性三維非穩(wěn)定流含水層。

圖3 模型剖分立體示意圖

（2）邊界條件概化。

項(xiàng)目所處的湖田水文地質(zhì)單元是一個相對獨(dú)立的水文地質(zhì)單元，本次研究的目的層為奧陶系碳酸鹽類巖溶～裂隙含水層，呈層狀展布，半開放型。南部丘陵裸露區(qū)為無壓區(qū)，節(jié)理、裂隙、巖溶較發(fā)育，為大氣降水的滲入提供了必要的賦存空間。近山前一帶隱伏區(qū)為承壓區(qū)，裂隙巖溶非常發(fā)育，使得從裸露區(qū)到覆蓋區(qū)含水層相互聯(lián)通，成為具統(tǒng)一地下水流場的連續(xù)介質(zhì)。將含水層概化為單層結(jié)構(gòu)，上邊界概化為承壓～半承壓層，下邊界為巖溶含水層底板，為隔水邊界。

根據(jù)區(qū)域地下水流場特征和地層結(jié)構(gòu)特點(diǎn)，設(shè)定北部湖田向斜、西部炒米莊斷裂均為隔水邊界［2］；東部金嶺斷層為弱透水?dāng)鄬樱?］，受堠皋強(qiáng)排井控制，湖田與大武水文地質(zhì)單元沒有流量交換，將金嶺斷層定義為隔水邊界；東部南側(cè)邊界與地下水等水頭線垂直，定為零流量邊界；南部邊界為地表水與分水嶺，為零流量邊界。

（3）網(wǎng)格剖分。

結(jié)合研究區(qū)的地形地貌特征、含水層特征、水文地質(zhì)條件及模型邊界條件，采用統(tǒng)一的平面網(wǎng)格。網(wǎng)格劃分間距20 ～100m，平面上單層共劃分432448 個單元，其中有效單元345591 個。

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

（1）地下水流模型。

根據(jù)研究區(qū)水文地質(zhì)條件，可建立研究區(qū)地下水流模型如下：

式中： h(x，y，z，t)為含水層的水位分布（m）；kx，ky，kz為空間三個方向的滲透系數(shù)（m/d）； Ss為含水層的單位儲水系數(shù)（1/m）；P 為單位體積源匯項(xiàng)[m3/(m3·h) ]；h0(x，y，z)為地下水初始水位分布；h1(x，y，z，t) 為含水層的一類邊界條件；q 為含水層二類邊界單位面積流量[m3/(m2·h) ]； θx，θy，θz為二類邊界的外法線方向與x，y，z 軸的夾角；Ω 為滲流區(qū)域； г1為一類邊界； г2為二類邊界；г3為自由面。

（2）地下水溶質(zhì)運(yùn)移數(shù)學(xué)模型。

根據(jù)研究區(qū)的具體條件，采用下述的溶質(zhì)運(yùn)移模型：

式中： c 為溶解于水中的污染物濃度（M/L3）；Dij為水動力彌散系數(shù)張量（T/L2）； Xi為空間坐標(biāo)（L）； V1為地下水滲透流速（L/T）； t 為時間（T）；n 為孔隙度，無量綱； F 為固相表面的溶質(zhì)濃度（M/L3）； C0為源匯項(xiàng)的濃度（M/L3）； qW為源（正值）或匯（負(fù)值）的單位流量（L/T）； Ω 為研究區(qū)空間區(qū)域；（x，y，z） 為空間位置； г 為研究區(qū)的邊界。

3.3 模型參數(shù)選擇

模型采用的主要參數(shù)設(shè)定如表1。

表1 模型中各參數(shù)取值表

3.4 模型識別檢驗(yàn)

模型識別與檢驗(yàn)是數(shù)值模擬工作及模型建設(shè)過程中最為關(guān)鍵的一個環(huán)節(jié)。在模型識別與檢驗(yàn)過程中，重新認(rèn)識水文地質(zhì)概念模型、分析水文地質(zhì)條件、進(jìn)一步提升對水文地質(zhì)模型的認(rèn)識。根據(jù)實(shí)測的水位及動態(tài)變化數(shù)據(jù)，對模型進(jìn)行識別和調(diào)試。

（1）模型識別。

利用2007 年1 月1 日-2012 年7 月1 日的長測資料，進(jìn)行數(shù)值模型的識別，從地下水位動態(tài)監(jiān)測井的動態(tài)資料的耦合對比分析，模擬的地下水動態(tài)符合實(shí)際，體現(xiàn)了地下水水位動態(tài)規(guī)律，見圖3。

圖4 模型識別期鉆孔水位擬合曲線圖

（2）模型驗(yàn)證。

為了檢驗(yàn)識別后模型的可靠性，利用2012 年7 月1 日-2017 年7 月1 日期間的數(shù)據(jù)進(jìn)行模型驗(yàn)證，將實(shí)測水位資料與計算水位資料進(jìn)行擬合對比，見圖5，對模型識別期間求取的水文地質(zhì)參數(shù)進(jìn)行校正，反復(fù)調(diào)試水文地質(zhì)參數(shù)、各均衡項(xiàng)，識別水文地質(zhì)條件，最終確定了各水文地質(zhì)參數(shù)及均衡項(xiàng)。經(jīng)過識別、驗(yàn)證后的模型基本能夠刻劃地下水系統(tǒng)的滲流特征，反映出水位動態(tài)變化對源、匯項(xiàng)的響應(yīng)關(guān)系，可以進(jìn)行預(yù)測分析。

圖5 模型驗(yàn)證期鉆孔水位擬合曲線圖

4 地下水環(huán)境影響分析

項(xiàng)目設(shè)計考慮了雨污分流、清污分流、污污分治，將罐區(qū)、裝卸車區(qū)等雨水收集系統(tǒng)單獨(dú)設(shè)計，初期雨水與后期清凈雨水采用溢流式自動切換裝置，將初期雨水收集起來匯入收集池，用污水提升泵限量送至污水處理站進(jìn)行處理，后期雨水進(jìn)入廠區(qū)雨水管網(wǎng)收集至雨水監(jiān)控池，由企業(yè)雨水管網(wǎng)外排。生活污水、生產(chǎn)廢水、地面沖洗水、機(jī)泵冷卻水經(jīng)企業(yè)污水處理站處理后排入水質(zhì)凈化廠進(jìn)行處理。項(xiàng)目廢水最終經(jīng)處理達(dá)到《城鎮(zhèn)污水處理廠污染物排放標(biāo)準(zhǔn)》（GB18918-2002）一級A 標(biāo)準(zhǔn)要求，排入河流。

4.1 預(yù)測因子與標(biāo)準(zhǔn)

根據(jù)企業(yè)污水中各組分含量，選取COD 作為污染物的代表性因子進(jìn)行預(yù)測。將預(yù)測濃度超過《地下水質(zhì)量標(biāo)準(zhǔn)》（GB/T14848-2017）III 類標(biāo)準(zhǔn)限值定義為超標(biāo)濃度，將指標(biāo)的檢出限作為污染影響判定標(biāo)準(zhǔn)。則COD 超標(biāo)濃度為3.0mg/L，影響濃度為0.05mg/L。

4.2 地下水環(huán)境影響預(yù)測

（1）正常工況下地下水環(huán)境影響預(yù)測。

正常工況下，項(xiàng)目嚴(yán)格按照設(shè)計要求落實(shí)好環(huán)保、防滲措施和管理措施，廢水收集、存放、排放系統(tǒng)密閉，不產(chǎn)生無組織排放。但由于污水收集池防滲層的滲透系數(shù)為10-7cm/s，會有極少量的污水滲入地下。經(jīng)模擬預(yù)測，項(xiàng)目及附近地下水中COD 含量均未出現(xiàn)超標(biāo)，也未出現(xiàn)影響范圍。

（2）事故工況下地下水環(huán)境影響預(yù)測。

“跑、冒、滴、漏”工況：用于儲存和處理污水的收集池池底防滲層出現(xiàn)“跑、冒、滴、漏”等無組織排放，污水泄露量按企業(yè)廢水產(chǎn)生量的1%計。經(jīng)模型預(yù)測，滲入地下的污水在地下水流和堠皋強(qiáng)排井的雙重作用下向東北方向擴(kuò)散。圖6 為事故發(fā)生后5 年和10 年COD 濃度分布圖。項(xiàng)目區(qū)及周邊地下水中COD 含量出現(xiàn)明顯超標(biāo)現(xiàn)象，COD濃度最高達(dá)到200mg/L。事故發(fā)生5 年后，地下水中COD 影響范圍達(dá)到6.35×105m2，最大影響距離為1410m，此時COD 超標(biāo)范圍達(dá)到1.75×105m2，最大超標(biāo)距離為710m。事故發(fā)生10 年后COD的影響范圍達(dá)到9.35×105m2，最大影響距離為2680m，此時超標(biāo)范圍達(dá)到2.95×105m2，最大超標(biāo)距離為1270m。

圖6 “跑、冒、滴、漏”工況不同時段COD 濃度分布圖

圖7 突發(fā)事故工況不同時段COD 濃度分布圖

突發(fā)事故工況：假定事故工況下，廢水收集池池底防滲層出現(xiàn)20m2的破損，事故連續(xù)10 天。預(yù)測結(jié)果見圖7。經(jīng)模擬預(yù)測，事故發(fā)生后1 天地下水中COD 含量就開始超過3mg/L，之后超范圍不斷增大，到577 天超標(biāo)范圍達(dá)到最大，最大超標(biāo)面積約1.5×105m2，此時影響面積為6.55×105m2，最遠(yuǎn)影響距離為1410m。之后污染緩慢消退，超標(biāo)范圍同步減小，但影響范圍和影響距離不斷增大。持續(xù)到第850 天后污染超標(biāo)范圍消失，地下水中COD 含量小于3mg/L，滿足地下水水質(zhì)Ⅲ類標(biāo)準(zhǔn)。事故發(fā)生1000 天后影響范圍達(dá)到最大，為1.13×106m2，最大影響距離為2550m。事故發(fā)生1700 天后，影響范圍消失，此次污染事故對地下水環(huán)境影響消失。

4 結(jié)論

（1）模型預(yù)測結(jié)果表明，在嚴(yán)格落實(shí)地面防滲措施的條件下，項(xiàng)目周邊地下水中COD 含量未出現(xiàn)超標(biāo)，項(xiàng)目對地下水環(huán)境影響較小。“跑、冒、滴、漏”和突發(fā)事故狀態(tài)下，最遠(yuǎn)影響距離超過2km，滲漏污水對項(xiàng)目周邊及下游地下水環(huán)境將構(gòu)成明顯威脅。

（2）由于區(qū)域巖溶地下水埋藏較深，一旦污染，很難修復(fù)。項(xiàng)目日常管理過程中應(yīng)建立污水巡查制度和地下水監(jiān)測系統(tǒng)，如發(fā)現(xiàn)問題及時處理，最大限度保護(hù)區(qū)域地下水資源。

（3）鑒于巖溶水系統(tǒng)的復(fù)雜性，造成數(shù)值模擬的不確定性，從長遠(yuǎn)考慮建議另選廠址。