北京某擬建再生水廠地下水環(huán)境影響研究

李 鵬, 劉 殷, 孫 穎

（北京市水文地質(zhì)工程地質(zhì)大隊(duì), 北京 100195）

北京某擬建再生水廠地下水環(huán)境影響研究

李 鵬, 劉 殷, 孫 穎

（北京市水文地質(zhì)工程地質(zhì)大隊(duì), 北京 100195）

研究區(qū)位于潮白河沖洪積扇中上部，第四系松散孔隙含水層顆粒粗大，易接受補(bǔ)給，地下水固有防污性能處于差和較差，地下水水質(zhì)整體良好。采用數(shù)值模擬方法，建立了三維地下水流數(shù)值模型，對(duì)孔隙度、彌散系數(shù)進(jìn)行了率定，利用GMS中的MT3DMS模塊建立了溶質(zhì)運(yùn)移模型。選取COD作為模擬指標(biāo)，根據(jù)模型對(duì)水廠投產(chǎn)30年后兩種不同情景方案下地下水水質(zhì)狀況進(jìn)行了預(yù)測(cè)分析。結(jié)果顯示，在正常工況和處理流程出現(xiàn)故障兩種條件下，再生水廠對(duì)周圍地下水環(huán)境有一定程度影響，但重點(diǎn)研究區(qū)內(nèi)地下水模擬濃度遠(yuǎn)小于《地下水質(zhì)量標(biāo)準(zhǔn)》中Ⅲ類水標(biāo)準(zhǔn)，污染物濃度雖低，影響距離卻很遠(yuǎn)，因此要建立覆蓋全廠區(qū)的地下水長(zhǎng)期監(jiān)控系統(tǒng)，避免污染事故。

溶質(zhì)運(yùn)移模型；水文地質(zhì)條件；再生水廠；地下水環(huán)境

0 引言

北京市人均水資源占有量為191m3（北京市水務(wù)局，2013），人多水少是北京的基本市情水情。因此，實(shí)現(xiàn)水資源的合理配置、優(yōu)水優(yōu)用、分質(zhì)供水、城市污水回用是解決水資源短缺的一個(gè)重要途徑（呂曉儉等，2015；黃大英，2005；張玲，2007）。擬建的再生水廠經(jīng)污水集中處理后，根據(jù)規(guī)劃將部分退水作為綠化、道路澆灑、工業(yè)冷卻水、河道補(bǔ)水等多種用途。建設(shè)地點(diǎn)位于密云縣經(jīng)濟(jì)開(kāi)發(fā)區(qū)，結(jié)合密云工業(yè)開(kāi)發(fā)區(qū)的發(fā)展規(guī)劃和用地布局，確定項(xiàng)目廠址位于潮白河下游南岸，工業(yè)三期西南角，北側(cè)為左堤路和潮白河，南側(cè)為京承高速公路，東側(cè)為西統(tǒng)路，西側(cè)為綠化帶。

由于該區(qū)第四系顆粒較粗，包氣帶防污性能弱，且位于水源地補(bǔ)給區(qū)上游，下游有八廠水源地。因此，本次確定研究范圍為696km2，西北側(cè)為潮白河，沿潮白河至下游水源八廠，以一個(gè)大的水文地質(zhì)單元為整體；重點(diǎn)評(píng)價(jià)區(qū)范圍為86.8km2，包括潮白河周邊的灌溉區(qū)、水源八廠補(bǔ)給區(qū)等（圖1）。通過(guò)模擬預(yù)測(cè)，評(píng)判再生水廠運(yùn)行過(guò)程中對(duì)地下水環(huán)境可能造成的影響，并針對(duì)可能的影響提出污染防治措施與建議，預(yù)防與控制地下水環(huán)境惡化，保護(hù)區(qū)域地下水安全。

圖1 區(qū)域水文地質(zhì)圖Fig.1 Regional hydrogeological map

1 地下水環(huán)境現(xiàn)狀調(diào)查與評(píng)價(jià)

1.1 研究區(qū)概況

密云縣位于北京市東北部，密云水庫(kù)位于其中，縣境內(nèi)主要河流有潮河、白河、潮白河等。研究區(qū)東部、西部和北部三面環(huán)山，南為平原，總的趨勢(shì)是北部狹窄，南部開(kāi)闊。牛欄山以北及河流兩側(cè)為由全新世沖積、洪積物組成的一級(jí)階地；牛欄山以南的一級(jí)階地外側(cè)分布有更新世沖積、洪積物組成的二級(jí)階地；潮白河等現(xiàn)代河流兩側(cè)有現(xiàn)代沖積層。區(qū)內(nèi)的斷裂構(gòu)造有北東向牛欄山斷裂，南部分布馬昌營(yíng)斷裂。區(qū)內(nèi)第四系為沖洪積成因類型地層，巖性為卵石、礫石、砂類土與粘性土，其中砂類土以中、粗砂為主，細(xì)砂、粉砂次之，粘性土中以砂質(zhì)粘土為主。沉積厚度由北部溪翁莊、東北部河南寨的20～50m，向東南逐漸加大到大胡營(yíng)的200m，在兩河附近最大，超過(guò)300m。下伏基巖地層包括長(zhǎng)城系、薊縣系、青白口系、寒武系、奧陶系、石炭—二疊系和侏羅—白堊系。廠址區(qū)包氣帶巖性平面上主要為砂類土，剖面上從上向下依次為卵石、礫石、砂類土與粘性土，下層粘性土主要是粉質(zhì)粘土和粘質(zhì)粉土，厚度4～8m，且分布連續(xù)、穩(wěn)定。

1.2 水文地質(zhì)條件

研究區(qū)位于潮白河地下水子系統(tǒng)中上部，第四系孔隙水主要賦存于河流沖洪積作用形成的砂及砂卵礫石中，含水層垂直滲透性強(qiáng)，極易接受大氣降水和地表水的補(bǔ)給，地下水和地表水聯(lián)系密切，是補(bǔ)給地下水的極佳地段，但其含水層具有易污染特征。根據(jù)含水層巖性、富水性及埋藏條件將研究區(qū)劃分為三個(gè)主區(qū)，單井出水量在1500～5000 m3/d之間，而山麓坡積裙地帶，第四系沉積物主要為坡積、洪積物，富水性不均一（圖1）。整個(gè)研究區(qū)內(nèi)滲透系數(shù)的變化范圍為3～300m/d，給水度的變化范圍為0.03～0.23，降雨入滲系數(shù)變化范圍為0.15～0.65。根據(jù)場(chǎng)地的滲水試驗(yàn)結(jié)果，其滲透系數(shù)在6.88×10-4～8.55×10-4cm/s，地下水固有防污性能處于差和較差區(qū)（北京市地質(zhì)礦產(chǎn)勘查開(kāi)發(fā)局等，2008）。

研究區(qū)孔隙水的補(bǔ)給來(lái)源主要是大氣降水、農(nóng)田灌溉水、地表水（河道、渠道水）入滲補(bǔ)給和山區(qū)基巖裂隙水側(cè)向補(bǔ)給，近年來(lái)人工回灌也成為地下水補(bǔ)給來(lái)源之一。地下水位變化除受人工開(kāi)采影響外，還受降水量大小影響。孔隙水的排泄以人工開(kāi)采方式為主，其次是側(cè)向徑流排泄。由于受到懷柔應(yīng)急水源地及水源八廠常年抽水影響，加之潮白河已干涸多年，不存在補(bǔ)給潛水的穩(wěn)定地表水，造成了潮白河西岸潛水位標(biāo)高整體低于東岸。地下水徑流方向與地形地貌變化一致，即由山前向平原，由北西向南東流動(dòng)，并隨徑流條件變化，徑流強(qiáng)度發(fā)生變化。現(xiàn)狀條件下，八廠水源地及懷柔應(yīng)急水源地地下水集中開(kāi)采改變了其自然狀態(tài)，在開(kāi)采量較大的地區(qū)，形成了大面積的區(qū)域性降落漏斗（圖2）。

圖2 地下水邊界條件及流場(chǎng)擬合圖Fig.2 Boundary conditions and ftting curve of groundwater

1.3 地下水環(huán)境現(xiàn)狀

根據(jù)重點(diǎn)研究區(qū)內(nèi)近一個(gè)水文年三期潛水含水層的水樣監(jiān)測(cè)結(jié)果，采用單因子指數(shù)法進(jìn)行了水質(zhì)評(píng)價(jià)，標(biāo)準(zhǔn)采用《地下水質(zhì)量標(biāo)準(zhǔn)》（GB/ T14848－93）中的Ⅲ類標(biāo)準(zhǔn)，標(biāo)準(zhǔn)中沒(méi)有的因子參照《生活飲用水衛(wèi)生標(biāo)準(zhǔn)》（GB 5749－2006）。結(jié)果顯示，25項(xiàng)檢測(cè)指標(biāo)中23項(xiàng)濃度均優(yōu)于上述標(biāo)準(zhǔn)值，水質(zhì)良好。總硬度、硝酸鹽（以N計(jì)）等兩個(gè)指標(biāo)水質(zhì)稍差，位于河南寨鎮(zhèn)西南部附近，而其上下游的水質(zhì)監(jiān)測(cè)點(diǎn)兩指標(biāo)的濃度均優(yōu)于標(biāo)準(zhǔn)值。

2 地下水環(huán)境影響預(yù)測(cè)與評(píng)價(jià)

2.1 地下水流數(shù)值模型

數(shù)值模型是某些復(fù)雜的物理實(shí)體及其內(nèi)部流程的簡(jiǎn)化和概化（孫納正，1981）。建立地下水流模型的過(guò)程，就是建立和運(yùn)行能代表實(shí)際含水層行為的模型（James，1981）。本次借助基于有限差分法原理的GMS軟件建立地下水?dāng)?shù)值模型（Brigham Young University-Environmental Modeling Research Laboratory，2002；祝曉彬，2003）。

由地下水位等值線（圖2）可以看出，研究區(qū)地下水流總體由北向南流動(dòng)，在開(kāi)采強(qiáng)度大的地區(qū)，形成區(qū)域地下水位降落漏斗。地下水在第四系松散多孔介質(zhì)中的流動(dòng)符合質(zhì)量守恒定律和達(dá)西定律；考慮到由層間水頭差異引起含水層之間的垂向水量交換，故地下水運(yùn)動(dòng)為三維流；地下水的補(bǔ)排項(xiàng)以及水位是隨時(shí)間變化的，故為非穩(wěn)定流；由于介質(zhì)的非均勻性造成水文地質(zhì)參數(shù)隨空間變化，體現(xiàn)了系統(tǒng)的非均質(zhì)性；由于含水介質(zhì)的成層性，造成垂向和水平水文地質(zhì)參數(shù)的差異，因而可概化為水平各向同性、垂向各向異性介質(zhì)。將重點(diǎn)研究區(qū)及以北區(qū)域模型概化為一層單一的潛水含水層，以南的區(qū)域含水層較厚，未處理為單一含水層。

根據(jù)研究區(qū)水文地質(zhì)條件和地下水開(kāi)發(fā)利用特點(diǎn)，側(cè)向邊界設(shè)置如圖2所示，垂向上第四系松散層的底界處理為隔水邊界。源匯項(xiàng)數(shù)據(jù)均按照鄉(xiāng)鎮(zhèn)分為22個(gè)區(qū)進(jìn)行整理，之后輸入模型中進(jìn)行計(jì)算。將模擬區(qū)各層采用100m×100m的剖分格式。地下水流場(chǎng)模擬期為2010年12月到2011年12月，以2010年12月地下水流場(chǎng)作為模型的初始流場(chǎng)。滲透系數(shù)等水文地質(zhì)參數(shù)處理為shp區(qū)文件直接輸入模型。

利用研究區(qū)2011年12月的實(shí)際流場(chǎng)與模擬的2011年12月的流場(chǎng)，進(jìn)行識(shí)別驗(yàn)證。可以看出，模擬流場(chǎng)與實(shí)際流場(chǎng)總體流動(dòng)方向相同，降落漏斗的影響范圍和最大降深相同，擬合效果良好。

2.2 預(yù)測(cè)指標(biāo)及情景設(shè)定

擬建再生水廠一期（2015年）建設(shè)規(guī)模6.5×104m3/d，二期（2020年）建設(shè)規(guī)模3.5×104m3/d，總規(guī)模10.0×104m3/d。工程處理出水達(dá)到《城市污水再生利用 城市雜用水水質(zhì)》（GB/T18920-2002）中的城市綠化和道路清掃的水質(zhì)標(biāo)準(zhǔn)，《城市污水再生利用 景觀用水水質(zhì)》（GB/T18921-2002）中的觀賞性景觀環(huán)境用水（河道類）的水質(zhì)標(biāo)準(zhǔn)，同時(shí)考慮滿足北京市地方標(biāo)準(zhǔn)《城鎮(zhèn)污水處理廠水污染物排放標(biāo)準(zhǔn)》 （DB11/890-2012）表 1 中的 A 標(biāo)準(zhǔn)。設(shè)計(jì)出水水質(zhì)見(jiàn)表1。

表1 再生水廠進(jìn)出水水質(zhì)指標(biāo)Tab.1 The designed quality of raw wastewater and treated water

擬建再生水廠的進(jìn)水水質(zhì)指標(biāo)中COD的濃度最高，初始濃度500mg/L，故選用COD作為模擬指標(biāo)。設(shè)定出水濃度為20mg/L。研究區(qū)位于潮白河沖洪積扇的中上部，這些區(qū)域地下水循環(huán)快，相對(duì)隨地下水流的排泄來(lái)說(shuō)，在吸附、生物降解等作用下消耗掉的量顯得微弱， 所以本次模擬僅考慮了對(duì)流-彌散作用。

溶質(zhì)運(yùn)移模型的預(yù)測(cè)期為2015年1月至2045年1月，共30年。不同情景污染源強(qiáng)按以下兩種方案確定：

（1）正常工況下即再生水廠正常運(yùn)轉(zhuǎn)

預(yù)測(cè)源強(qiáng)設(shè)定為COD出水濃度20mg/L。

（2）再生水廠處理流程出現(xiàn)故障

再生水廠突發(fā)污染事件，如爆炸等。進(jìn)水未經(jīng)過(guò)處理，直接由再生水廠補(bǔ)給地下水造成污染，模擬其對(duì)地下水的影響。進(jìn)水水質(zhì)的COD濃度為500mg/L，假定以此濃度泄漏持續(xù)時(shí)間15天，通過(guò)計(jì)算污染物滲漏量為4814.8g/d。

2.3 地下水水質(zhì)模型預(yù)測(cè)與評(píng)價(jià)

在地下水流數(shù)值模型的基礎(chǔ)上利用GMS軟件中的MT3DMS建立溶質(zhì)運(yùn)移模型（鄭春苗等，2009），對(duì)污染物的運(yùn)移進(jìn)行預(yù)測(cè)評(píng)價(jià)。溶質(zhì)模型參數(shù)設(shè)定，孔隙度為0.3；縱向彌散系數(shù)，中粗砂設(shè)定為0.2～1.0 m2/d，砂礫石設(shè)定為1～5 m2/d。

（1）正常工況下COD運(yùn)移預(yù)測(cè)

再生水長(zhǎng)期按照出水標(biāo)準(zhǔn)排入潮白河，周邊利用這類水灌溉，做景觀用水。預(yù)測(cè)結(jié)果如圖3所示。污染指標(biāo)COD隨著時(shí)間的增加，濃度越來(lái)越高，污染暈的范圍越來(lái)越大，并且沿著地下水流動(dòng)的方向污染暈的影響距離大于其他方向，最大影響距離為6100m。2020年6月污染物運(yùn)移到水源八廠區(qū)，至2045年1月水源八廠區(qū)污染物濃度最大值為1.66mg/L，小于地下水質(zhì)量標(biāo)準(zhǔn)中Ⅲ類水COD的濃度指標(biāo)3mg/L。故污染物COD對(duì)地下水環(huán)境僅有微弱的影響。

圖3 正常工況下COD污染擴(kuò)散圖Fig.3 Distribution of COD plume in normal operation condition

（2）再生水廠處理流程出現(xiàn)故障

污水COD未經(jīng)過(guò)處理的預(yù)測(cè)結(jié)果如圖4所示。由于COD的初始濃度高，污染物運(yùn)移的速率明顯大于正常工況下污染物的運(yùn)移速率，最大影響距離為8700m。2038年6月污染物運(yùn)移到水源八廠區(qū)，2045年1月水源八廠區(qū)污染物濃度最大值0.22mg/L，遠(yuǎn)小于地下水Ⅲ類水標(biāo)準(zhǔn)。故再生水廠處理流程出現(xiàn)故障時(shí)，污染物COD對(duì)地下水環(huán)境影響甚微。

3 結(jié)語(yǔ)

通過(guò)對(duì)區(qū)域地質(zhì)、水文地質(zhì)條件分析，廠址區(qū)地下水固有防污性能處于差和較差區(qū)，含水層易污染特征為易，采樣分析評(píng)價(jià)地下水環(huán)境現(xiàn)狀，大部分指標(biāo)水質(zhì)良好。利用GMS軟件，建立了地下水流數(shù)值模擬模型及溶質(zhì)運(yùn)移模型。預(yù)測(cè)再生水廠投入使用30年后，不同情景方案下地下水水質(zhì)狀況，結(jié)果顯示再生水廠的使用對(duì)周圍地下水環(huán)境有一定程度影響，但區(qū)內(nèi)主要供水水源地COD濃度均遠(yuǎn)小于《地下水質(zhì)量標(biāo)準(zhǔn)》中Ⅲ類水標(biāo)準(zhǔn)。污染物濃度雖低，影響距離卻很遠(yuǎn)，仍要引起重視。

圖4 處理流程故障條件下COD污染擴(kuò)散圖Fig.4 Distribution of COD plume when process fails

南水北調(diào)來(lái)水后在研究區(qū)進(jìn)行了地下水回灌涵養(yǎng)，地下水位降落漏斗得以恢復(fù)，水動(dòng)力條件發(fā)生變化，污染物運(yùn)移將減緩，從而也可以改善研究區(qū)的地下水環(huán)境。

宜建立覆蓋全廠區(qū)的地下水長(zhǎng)期監(jiān)控系統(tǒng)，對(duì)地下水污染狀況及時(shí)發(fā)現(xiàn)并控制，避免污染事故的發(fā)生。若廠區(qū)發(fā)生地下水污染事故，應(yīng)立即啟動(dòng)應(yīng)急預(yù)案。如控制廠區(qū)事故現(xiàn)場(chǎng)，切斷污染源，進(jìn)行清源處理，同時(shí)上報(bào)相關(guān)部門。對(duì)每個(gè)可能出現(xiàn)事故的污染暈設(shè)置截流帶，布設(shè)截獲井，直至地下水污染物含量低于區(qū)域背景值后停止抽水。

[1]Brigham Young University-Environmental Modeling Research Laboratory. Groundwater Modeling System Tutorials Volume I, 2002.

[2]James W. Mercer，Charles R.Faust. Ground-Water Modeling[M]. National Water Well Association，1981.

[3]北京市水務(wù)局. 北京市水資源公報(bào)(2012)[R]. 北京，2013.

[4]呂曉儉，李宇，李志萍. 北京地下水資源與城市協(xié)調(diào)發(fā)展的思考[J]. 城市地質(zhì)，2015，10(S1)：8～12. [5]黃大英. 淺談北京水資源的合理配置[J]. 水資源保護(hù)，2005，21(4)：67～70.

[6]孫納正. 地下水流的數(shù)學(xué)模型和數(shù)值方法[M]. 北京：地質(zhì)出版社，1981.

[7]北京市地質(zhì)礦產(chǎn)勘查開(kāi)發(fā)局，北京市水文地質(zhì)工程地質(zhì)大隊(duì). 北京地下水[M]. 北京：中國(guó)大地出版社，2008.

[8]張玲. 基于不確定性的水資源優(yōu)化配置模型及案例研究[D]. 北京：中國(guó)地質(zhì)大學(xué)(北京)，2007.

[9]鄭春苗，Bennett G D. 地下水污染物遷移模擬[M].北京：高等教育出版社，2009.

[10]祝曉彬. 地下水模擬系統(tǒng)(GMS)軟件[J]. 水文地質(zhì)工程地質(zhì)，2003，47 (5)：53～55.

Research on the Effect of Groundwater Environment for a Proposed Water Recycling Plant in Beijing

LI Peng, LIU Yin, SUN Ying
(Beijing Institute of Hydrogeology and Engineering Geology, Beijing, 100195)

The study area is located in the upper of the Chaobai River alluvial. Because of the coarse particle size, Quaternary aquifer is easily recharged. Groundwater intrinsic antipollution capability of this area is poor. Groundwater quality is good in general. The index COD is chosen to predict and evaluate the effects of a proposed water recycling plant operation on groundwater in the source area by building up numerical model. The plant may have certain effect on the groundwater environment both in cases of normal operation and during process failures. The forecasting COD concentrations in the groundwater source area are far below the Class III water standard in Standard for Groundwater Quality, but the infuence distance is long, therefore a long-term monitoring system covering the whole plant area should be build to avoid pollution.

Mass transport model; Hydrogeologic conditions; Water recycling plant; Groundwater environment

X523

A

1007-1903（2016）04-0024-05

10.3969/j.issn.1007-1903.2016.04.004

李鵬（1983- ），男，碩士，工程師，主要從事水文地質(zhì)環(huán)境地質(zhì)工作。E-mail：liplip123@163.com.