安雙財

摘要：在我國西南基巖山區(qū)存在許多巖溶含水層-裂隙弱透水層交替分布區(qū)，這些地區(qū)巖溶含水層巖溶較為發(fā)育，而裂隙弱透水層隔水性相對較好，因此，開展這類基巖山區(qū)地下水污染風險分析對巖溶含水層水源保護及項目選址具有重要指導意義。本文以云南富源縣某擬建化工項目為例，通過污染項目對其下伏及周邊含水層的影響進行預測分析，為基巖山區(qū)此類項目選址提供參考。

Abstract： There are many alternating distribution of karst aquifer and fissure aquitard in the Bedrock Mountains in southwest of China. The development of karst in the karst aquifer is very well and the permeable of fissure aquitard is relatively well. Therefore， the risk analysis for groundwater contamination in this kind of karst aquifer is very important for the water resource protection and project site selection. Taking a proposed chemical project in Fuyuan county， Yunnan province as an example， this study predicts and analyzes the influence of the pollution project on the underlying and surrounding aquifers， so as to provide reference for the site selection of such projects in bedrock mountain area.

關(guān)鍵詞：巖溶含水層;弱透水層;地下水污染;風險分析

Key words： karst aquifer;aquitard;groundwater contamination;risk analysis

中圖分類號：X824? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻標識碼：A? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文章編號：1006-4311（2020）14-0274-04

0? 引言

地下水污染風險評估是地下水污染管理的關(guān)鍵環(huán)節(jié)。準確評估一個水源地或污染源周邊的地下水污染風險對地下水污染管理至關(guān)重要。因此，近年來，地下水污染風險評估受到了許多學者的廣泛關(guān)注[1-7]。針對不同場地不同尺度的地下水污染風險評估采用的方法各有差異。對于平原區(qū)而言，李紹飛等[2]從污染源荷載、滲透補給條件、含水層巖性、社會經(jīng)濟等方面構(gòu)建評價地下水污染風險指標體系，根據(jù)突變理論，建立平原區(qū)淺層地下水污染風險突變評價模型。而對于巖溶區(qū)而言，崔亞豐等[8]構(gòu)建了適合巖溶區(qū)地下水污染風險評價體系并以地蘇地下河系流域為例開展巖溶區(qū)域地下水污染風險評價，研究表明地下水污染風險空間分布特征與巖性、坡度、巖溶網(wǎng)絡(luò)發(fā)育等自然條件密切相關(guān)。而對于構(gòu)造對地下水污染風險的影響方面，梅向陽等[9，10]針對向斜構(gòu)造對地下水污染風險進行了分析，而王仁敏等[11]則探討了巨厚包氣帶巖溶區(qū)垃圾填埋場選址合理性及其對下游敏感保護目標的影響。

綜上所述，場地地質(zhì)及水文地質(zhì)條件對地下水污染風險評估結(jié)果影響明顯，而本文重點考慮巖溶含水層和裂隙弱透水層分布區(qū)某擬建化工項目對地下水污染風險的影響進行分析評價，為我國基巖山區(qū)重點污染類項目選址及地下水污染風險評估提供借鑒。

1? 研究區(qū)水文地質(zhì)條件特征

1.1 研究區(qū)水文地質(zhì)概況

研究區(qū)位于區(qū)域雨汪旋卷構(gòu)造東南邊界，整體上，研究區(qū)地下水主要受雨汪旋卷構(gòu)造及北東南西向褶皺構(gòu)造的控制，不同次級構(gòu)造區(qū)域地下水富集性差異較大。根據(jù)水文地質(zhì)初步調(diào)查發(fā)現(xiàn)研究區(qū)主要出露三疊系中統(tǒng)個舊組（T2ga）碳酸鹽巖巖溶含水層和三疊系永寧鎮(zhèn)組（T1yb）細砂巖及頁巖弱透水層，為一個相對獨立完整的水文地質(zhì)單元，水文地質(zhì)單元范圍內(nèi)巖溶含水層主要為三疊系中統(tǒng)個舊組灰、深灰色灰?guī)r巖溶含水層，而弱透水層主要為三疊系永寧鎮(zhèn)組黃綠、紫色細砂巖夾粉砂巖頁巖、泥灰?guī)r夾泥巖弱透水層，經(jīng)數(shù)值模擬分析顯示，受永寧鎮(zhèn)組弱透水層的影響，數(shù)值模擬評價范圍內(nèi)地下水主要沿著北東及南西兩個方向徑流，南西方向地下水主要向南側(cè)喜舊溪河排泄，沿著北東方向徑流地下水最終排泄于北東側(cè)黃泥河?？傮w上，研究區(qū)為一個相對獨立的完整水文地質(zhì)單元，因此，本次圈定如圖1黑色虛線范圍約28.49km2作為本次地下水風險分析的調(diào)查范圍，而白色虛線范圍11.72km2則為本次數(shù)值模擬預測分析范圍，根據(jù)水文地質(zhì)調(diào)查發(fā)現(xiàn)模擬計算區(qū)北西側(cè)以F1導水斷層斷層為定流量邊界，南東側(cè)以三疊系永寧鎮(zhèn)組黃綠、紫色細砂巖夾粉砂巖頁巖為隔水邊界，南西側(cè)以個舊組分水嶺為隔水邊界，與區(qū)域地下水流向相符，由此，整個數(shù)值模擬研究區(qū)構(gòu)成一個相對獨立完整的水文地質(zhì)單元。區(qū)域地下水主要受大氣降水入滲補給，沿著三疊系中統(tǒng)（T2ga）巖溶及裂隙介質(zhì)徑流，根據(jù)本次水文地質(zhì)調(diào)查，分析范圍內(nèi)，北側(cè)區(qū)域的地下水總體上沿著北東方向徑流，而南側(cè)箐門前一帶的地下水往南徑流，最終排泄于南側(cè)740號區(qū)域泉群，整體上，研究區(qū)地下水補徑排條件較為清楚。

1.2 研究區(qū)構(gòu)造水文地質(zhì)特征

F1逆斷層：位于場地北側(cè)（圖1），地表部分被覆土掩蓋，走向北東，傾向南東，傾角42～60°，斷裂兩側(cè)均為灰?guī)r，富水性強，透水性強，項目區(qū)地下水沿北東向排泄。

F2逆斷層：位于項目區(qū)北側(cè)，走向北東，傾向南東，傾角25～63°，項目區(qū)地段，斷層北側(cè)為灰?guī)r，南側(cè)為頁巖，頁巖具有一定的阻水作用，故導致北側(cè)的巖溶裂隙水在斷裂F1地段被南側(cè)頁巖阻斷，沿斷裂F1向北東向徑流，根據(jù)圖2落水洞位置推測沿著場地北西側(cè)F2斷層可能存在巖溶發(fā)育帶，尤其是該斷層與F3夾持的個舊組巖溶含水層三角區(qū)域3口深井（100m）均為揭露地下水位，表明該區(qū)域存在巖溶導水帶，F(xiàn)2在場地北西側(cè)導水性較強，場地西側(cè)地表水及地下水沿著溝谷向北側(cè)徑流，至F2處順著落水洞進入地下，沿著地下水暗河向北東側(cè)徑流排泄。

研究區(qū)受區(qū)域雨汪旋卷構(gòu)造影響，各斷層導水性和富水性差異較大，在雨汪旋卷構(gòu)造作用下，地下水主要沿著北東及南西向徑流，同時受旋卷構(gòu)造兩側(cè)褶皺構(gòu)造影響，研究區(qū)范圍內(nèi)形成局部富水盆地，整體上，雨汪旋卷構(gòu)造控制著研究區(qū)地下水補徑排條件特征。

1.3 研究區(qū)地下水補徑排條件

根據(jù)本次調(diào)查和區(qū)域地質(zhì)資料，區(qū)域地下水類型以碳酸鹽巖巖溶水、碳酸鹽巖與非碳酸鹽巖成夾層型的巖溶水，受雨汪旋卷構(gòu)造及北西南東兩側(cè)褶皺構(gòu)造的影響，根據(jù)水文地質(zhì)調(diào)查及勘察顯示項目場地所處三疊系中統(tǒng)（T2ga）碳酸鹽巖巖溶水含水層較為豐富，而場地中部分布的三疊系下統(tǒng)永寧鎮(zhèn)組頁巖隔水層透水性較差，受含水層及隔水層分布特征的影響，研究區(qū)范圍內(nèi)主要受近北東南西向斷層構(gòu)造的控制，地下水主要受大氣降水入滲補給，順著溶蝕裂隙及層間裂隙徑流，受東側(cè)個舊組第二段（T2gb）頁巖隔水底板的影響，研究區(qū)南西側(cè)箐門前水庫一帶形成相對較為富集的富水盆地，根據(jù)雨汪電廠門口散戶水井揭露地下水位埋深約為35m，而穩(wěn)定水位埋深約為15m，富水層位巖溶含水層具有一定承壓性，穩(wěn)定水位標高與項目場地西南側(cè)出露的Q1和GW3龍?zhí)兑约澳衔鱾?cè)箐門前水庫西側(cè)Q2龍?zhí)冻雎兜臉烁呦喈?，由此推斷，雨汪電廠門口及南側(cè)區(qū)域以箐門前水庫為中心形成較為豐富的富水盆地。除了雨汪電廠門口散戶抽取地下水作為水源排泄外，研究區(qū)地下水主要沿著兩個方向排泄，箐門前水庫西南側(cè)主要往西南方向徑流排泄，而項目區(qū)西側(cè)出露的GW3龍?zhí)都癚1排泄進入西側(cè)溝谷河流，往北徑流經(jīng)落水洞進入F2導水斷層地下暗河向北東方向徑流，最終排泄于北東側(cè)黃泥河，受不同斷層構(gòu)造導水性和阻水性差異影響，研究區(qū)北側(cè)地下水埋深及富集程度差異較大，整體上，本次研究區(qū)地下水主要沿著北東側(cè)黃泥河溝谷排泄，而南西側(cè)主要往南側(cè)發(fā)乃村南側(cè)泉群（圖1）及喜舊溪河排泄，地下水排泄主要以泉及抽水井抽水形式排泄，經(jīng)調(diào)查分析，研究區(qū)地下水補徑排條件較為清楚。

1.4 區(qū)域及場地巖溶發(fā)育及分布及特征

根據(jù)本次調(diào)查和區(qū)域地質(zhì)資料，區(qū)域地下水類型以三疊系中統(tǒng)（T2ga）碳酸鹽巖巖溶含水層為主（圖1），受區(qū)域雨汪旋卷構(gòu)造影響，研究區(qū)范圍內(nèi)各斷層導水性及富水性差異較大，根據(jù)調(diào)查項目北西側(cè)落水洞及3個深孔未見地下水可以推斷研究區(qū)范圍內(nèi)主要發(fā)育的巖溶為F2斷層一帶可能存在地下暗河沿著F2斷層向北東側(cè)發(fā)育，最終以泉的形式排泄于北東側(cè)黃泥河，受東南側(cè)（T2gb）頁巖隔水底板影響，項目南東側(cè)以箐門前水庫為中心，地下水形成較為富集的富水盆地。而南側(cè)地下水主要向發(fā)乃一帶的740號區(qū)域泉群排泄，研究區(qū)范圍內(nèi)巖溶發(fā)育分帶特征較為明顯，巖溶強發(fā)育區(qū)主要受雨汪旋卷構(gòu)造的控制。

2? 研究區(qū)地下水污染風險分析

根據(jù)前述水文地質(zhì)條件調(diào)查及分析概化，依據(jù)相關(guān)規(guī)范要求進行模型概化并進行數(shù)值模擬預測分析，下面將給出邊界條件概化及氨氮特征污染物在7200天后污染羽狀物擴散范圍圖如圖3和圖4所示，圖中坐標網(wǎng)格為每個間隔為500m。

從上述預測圖4可以看出，重點污染工程區(qū)氨氮污染羽狀物主要沿著北側(cè)和北西方向遷移擴散，雨水收集池及事故水池區(qū)域的氨氮污染羽狀物在7200天后向北西方向最遠遷移距離約為200m，未能到達GW3龍?zhí)都兴?，而液氨站超標濃度氨氮污染羽狀物?200天后向北側(cè)最遠遷移距離約為300m。

3? 結(jié)果與討論

通過上述數(shù)值模擬分析顯示，對于巖溶區(qū)化工項目的選址，如果巖溶含水層中間夾有裂隙弱透水層，那么盡可能將重度污染工程分部在裂隙弱透水層中，可大大減少項目風險條件下對地下水污染造成的風險。以本次擬建項目為例，位于中部裂隙弱透水層中的重點污染工程羽狀物遷移擴散速度較慢，對周邊的敏感目標影響較小，而如果將項目落在高滲透性的巖溶含水層中或巖溶管道分布區(qū)，那么項目對地下水造成的污染風險將不可控。但對于本項目而言，把重點污染工程區(qū)分布于裂隙弱透水層中，對于選址困難地區(qū)，即使是巖溶區(qū)，也不應是化工類項目的禁區(qū)。此外，從上述污染物遷移擴散方向及距離可以看出，受不同含水層/弱透水層分布影響，羽狀物擴散方向及距離差異較大，但最終未對項目周邊敏感目標造成污染影響，即超標濃度羽狀物在7200天內(nèi)未能到達周邊泉點及溶洞入口。綜合以上分析，對于我國西南巖溶區(qū)項目選址，應因地制宜，而不是一概而論，重點對項目區(qū)地質(zhì)及水文地質(zhì)條件進行深入調(diào)查及分析至關(guān)重要。只有在分析清楚水文地質(zhì)條件的基礎(chǔ)上，開展預測及評價工作，才有可能獲取更合理的評價結(jié)論，從而為項目區(qū)選址及項目環(huán)評提供更為可靠的依據(jù)。

另一方面，對于此類復雜水文地質(zhì)條件地區(qū)化工項目建設(shè)后續(xù)還應該做好地下水污染跟蹤監(jiān)測及應急替代措施，尤其是西側(cè)水源泉點，風險條件下，一旦監(jiān)測到該點受項目建設(shè)影響導致污染，應及時為其做好替代方案。此外，充分利用項目區(qū)地下水西側(cè)溶洞入口做好風險預案，一旦項目區(qū)污染地下水向該處匯集，可在該點采取收集處理措施。

4? 結(jié)論

本文以云南富源縣某巖溶含水層/裂隙弱透水層分布區(qū)某擬建化工項目地下水污染風險分析為例，通過水文地質(zhì)調(diào)查及勘察，在查明擬建污染項目區(qū)水文地質(zhì)條件基礎(chǔ)上進行項目區(qū)地下水污染風險分析及預測，通過研究主要獲得以下認識：

①通過解析計算及數(shù)值模擬分析表明，通過將高污染項目放在裂隙弱透水層中可降低項目建設(shè)對地下水污染風險。

②對于巖溶區(qū)化工項目的選址及地下水污染風險評估過程中，應根據(jù)項目區(qū)地質(zhì)及水文地質(zhì)條件進行充分對比分析，因地制宜，根據(jù)項目性質(zhì)及含水層/隔水層分布特征進行有效布局以及采取有效防控措施至關(guān)重要。

③巖溶區(qū)落水洞作為地表水/地下水匯集區(qū)，在極端情況下可直接在這些地區(qū)采取污染防治工程，減少污染地下水進一步擴散。

參考文獻：

[1]覃榮高，邱仁敏，黎明，曹廣祝，王金生，仵彥卿.包氣帶—含水層地下水污染風險評估研究進展[J].地球科學進展，2020，35（02）：111-123.

[2]李紹飛，陳伏龍，余萍，孫書洪.平原區(qū)淺層地下水污染風險評價[J].武漢大學學報（工學版），2018，51（12）：1035-1040，1049.

[3]張博，李國秀，程品，孫法圣，王洪義.基于隨機理論的地下水環(huán)境風險評價[J].水科學進展，2016（01）：100-106.

[4]孫才志，陳相濤，陳雪姣，鄭德鳳.地下水污染風險評價研究進展[J].水利水電科技進展，2015，35（05）：152-161.

[5]Locatelli L， Binning PJ， Sanchez-Vila X， S？ndergaard GL， Rosenberg L， Bjerg PL. A simple contaminant fate and transport modelling tool for management and risk assessment of groundwater pollution from contaminated sites. J Contam Hydrol 2019， 221： 35-49.

[6]Henri CV， Fernandez-Garcia D， de Barros FPJ. Probabilistic human health risk assessment of degradation-related chemical mixtures in heterogeneous aquifers： Risk statistics， hot spots， and preferential channels. Water Resour Res 2015， 51（6）： 4086-4108.

[7]de Barros FPJ， Fiori A. First-order based cumulative distribution function for solute concentration in heterogeneous aquifers： Theoretical analysis and implications for human health risk assessment. Water Resour Res 2014， 50（5）： 4018-4037.

[8]崔亞豐，何江濤，王曼麗，趙閱坤，王菲.巖溶地區(qū)地下水污染風險評價方法探究——以地蘇地下河系流域為例[J].中國巖溶 2016（04）：372-383.

[9]梅向陽，王菁，覃榮高，曹廣祝.云南曲靖市沾益區(qū)某向斜構(gòu)造巖溶區(qū)地下水污染數(shù)值模擬研究[J].價值工程，2019，38（05）：138-142.

[10]梅向陽，王菁，覃榮高，曹廣祝.云南富寧縣某冶煉項目所處背斜構(gòu)造區(qū)地下水污染風險分析[J].價值工程，2019，38（06）： 125-129.

[11]王仁敏，梅向陽，覃榮高，曹廣祝.含巨厚非飽和帶巖溶含水層地下水污染風險分析[J].價值工程，2019，38（12）：135-139.