劉 丹

(遼寧省河庫管理服務(wù)中心(遼寧省水文局)，遼寧 沈陽 110003)

地下水是區(qū)域水資源的重要組成部分，是經(jīng)濟(jì)社會發(fā)展的重要水源之一，對保障城鎮(zhèn)生活、工業(yè)用水和城鎮(zhèn)經(jīng)濟(jì)社會發(fā)展具有重要作用，也是保障生態(tài)與環(huán)境健康發(fā)展的重要要素[1]。地下水賦存于地質(zhì)介質(zhì)中，具有運動緩慢、補給周期長、循環(huán)更新慢、自我修復(fù)能力差、地下水系統(tǒng)遭到破壞后難以治理和修復(fù)等特點[2]。近年來，部分城鎮(zhèn)地下水源地由于集中、過量的開采，引起水源區(qū)地下水位持續(xù)下降，含水層出水能力降低，許多機電井報廢，致使地下水源地不斷地改建和擴(kuò)建，部分地區(qū)引發(fā)了海(咸)水入侵，出現(xiàn)了地下水水質(zhì)惡化現(xiàn)象[3]。地下水源地開發(fā)利用中存在的上述諸多問題，不僅嚴(yán)重影響了工業(yè)生產(chǎn)和人民生活的正常用水需求，對經(jīng)濟(jì)社會發(fā)展和生態(tài)環(huán)境構(gòu)成極大的危害，而且已嚴(yán)重危及地下水資源的可持續(xù)利用[4]。鞍山西郊區(qū)主要包括西郊水源地和鐵西水源地兩處地下水集中開采區(qū)，地下水開采量占區(qū)域地下水可開采量的90%以上，屬于地下水超采區(qū)[5]。根據(jù)《遼寧省地下水保護(hù)條例》《遼寧省禁止提取地下水規(guī)定》《遼寧省人民政府辦公廳關(guān)于印發(fā)遼寧省實行最嚴(yán)格水資源管理制度“十三五”工作方案和遼寧省“十三五”封閉地下水取水工程總體工作方案的通知》等，亟需對其地下水時空變化的研究。近些年來，信息熵在區(qū)域地下水時空演變分析中得到較為成熟的應(yīng)用[6-15]，但在遼寧地區(qū)還未得到具體應(yīng)用，為此本文基于信息熵原理，以鞍山西郊為研究區(qū)域，對區(qū)域地下水超采區(qū)時空演變特征進(jìn)行分析，為鞍山市西郊地下水保護(hù)區(qū)地下水資源合理開發(fā)、高效利用、有效保護(hù)和科學(xué)管理提供依據(jù)。

1 地下水超采區(qū)現(xiàn)狀位置

西郊區(qū)域地下水超采區(qū)主要由西郊水源地超采區(qū)和鐵西水源地超采區(qū)組成，目前兩個水源地超采區(qū)已形成一個統(tǒng)一的超采區(qū)。西郊地下水超采區(qū)位于運糧河以西，紅旗堡—韭菜臺以南，小鄭家臺—孫青堆子以東，鞍山—騰鰲公路以北。超采區(qū)涉及達(dá)道灣、寧遠(yuǎn)兩個鎮(zhèn)的12個自然村。目前，鞍山市自來水公司在西郊水源地少量集中開采地下水，取水井分布在雙樓臺村、孫周新村、岳家村、孫青村，水源地面積約15km2。

2 超采區(qū)水位動態(tài)評價

2.1 超采區(qū)地下水位時間動態(tài)特征

根據(jù)西郊地下水超采區(qū)形成過程及2005—2018年期間漏斗區(qū)平均最大地下水埋深統(tǒng)計情況可知，超采區(qū)地下水位降落漏斗最大水位埋深出現(xiàn)在2008年，2008—2012年漏斗區(qū)最大水位埋深小幅下降并表現(xiàn)出總體平穩(wěn)的趨勢，2012—2018年漏斗區(qū)最大水位埋深逐年減小，表明地下水位存在明顯上升趨勢，如圖1所示。

圖1 西郊地下水超采區(qū)漏斗中心逐年地下水埋深過程線

據(jù)統(tǒng)計，研究區(qū)2012—2018年漏斗區(qū)最大水位埋深變化率為1.30m/a。根據(jù)超采區(qū)地下水位觀測井的分布特征及數(shù)據(jù)有效性，選取了2012、2015、2016、2018年作為典型年對地下水水位變化特征進(jìn)行時空分布分析，以評價超采區(qū)在2012—2018年期間水位的總體變化特征。

2.2 豐、枯水期超采區(qū)地下水位變化趨勢

采用線性回歸分析法對各典型年1月和8月平均地下水位進(jìn)行計算，得到各年度枯水期和豐水期地下水位總體變化趨勢，如圖2所示。

由圖2可知，從2012年到2018年，枯水期地下水位整體呈上升趨勢，水位上升率為1.20m/a；豐水期地下水位有所波動，但整體仍呈上升趨勢，水位上升率為0.53m/a；枯水期水位上升率高于豐水期水位上升率，至2018年豐水期水位和枯水期水位相差不大，但枯水期水位稍高于豐水期水位，主要受區(qū)域地下水開采情況影響。

2.3 超采區(qū)地下水位時間信息墑變化特征

為更直觀地定量分析西郊超采區(qū)地下水位的時空變化規(guī)律，本次研究引入時間信息熵和時間序列信息熵理論對超采區(qū)地下水位動態(tài)特征進(jìn)行分析，結(jié)果如圖3所示。時間信息熵可以體現(xiàn)出研究時段內(nèi)地下水位的變化強度，熵值越大表示地下水位波動的強度越大，熵值越小表示地下水位波動的強度越?。粫r間序列信息熵可以體現(xiàn)出地下水位的變化趨勢，熵值為正表示地下水位呈上升趨勢，熵值為負(fù)表示地下水位呈下降趨勢，熵值絕對值的大小表示地下水位上升或下降的程度，絕對值越大表示地下水位上升或下降的越明顯，絕對值越小表示地下水位上升或下降的程度越小。

由圖3的時間信息熵分布圖可以看出，枯水期研究區(qū)西部地下水位波動較小，東部地下水位波動較大，變化強度最大的區(qū)域集中在2#、4#、5#、6#、18#、19#、20#、21#井附近，最大熵值為8.69。由時間序列信息熵分布圖可以看出，整個研究區(qū)地下水位多以上升為主，上升區(qū)域占總區(qū)域的99.996%，上升趨勢較大的區(qū)域集中在2#、4#、5#、6#、18#、19#、20#、21#井附近，只有0.004%的區(qū)域地下水位下降，位于3#井附近，可忽略不計。另外，豐水期地下水位波動大的區(qū)域集中在觀測井附近，2#、3#、5#、19#、21#井周圍水位波動較大，最大熵值為9.65，而枯水期最大熵值為8.69，表明豐水期地下水位波動強度大于枯水期；研究區(qū)地下水位上升區(qū)域主要分布在研究區(qū)中部16#井和東北部5#、18#、20#、21#井附近，占整個研究區(qū)域的34.516%，下降區(qū)域主要分布在西部和東南部，占整個研究區(qū)域的65.484%，其中下降最明顯的區(qū)域位于13#、14#井附近。

2.4 超采區(qū)地下水位空間信息墑變化特征

采用各觀測井2018年與2012年地下水位的差值，可以得出地下水位變化情況，如圖4所示。

圖2 各年度枯水期和豐水期地下水位總體變化趨勢

圖3 2012—2018年1月和8月地下水位時間熵及時間序列信息熵分布圖

圖4 2012—2018年1月和8月地下水位變化分布圖

由圖4可以得到，枯水期地下水位上升區(qū)域占整個研究區(qū)域的99.996%，下降區(qū)域可忽略不計，這與從時間序列信息熵得出的結(jié)果一致，枯水期地下水位上升均值為3.250m，最大值為7.190m，出現(xiàn)在5#井附近。豐水期地下水位上升區(qū)域占45.345%，下降區(qū)域占54.655%，與從時間序列信息熵得出的結(jié)果略有不同，這是因為地下水位變化圖只考慮了2012、2018年兩個初末時刻的地下水位，忽略了中間年份地下水位變化導(dǎo)致的水位波動。豐水期地下水位變化均值為0.698m，上升最大值為9.130m，出現(xiàn)在21#井附近，下降最大值為3.293m，出現(xiàn)在3#、13#、14#井附近?？偟膩砜矗S水期地下水位波動大于枯水期。

2.5 超采區(qū)地下水位變化等級信息墑變化特征

利用自然間斷點分級法將枯水期地下水位變化分為5個等級，并求出每個等級對應(yīng)的時間序列信息熵的熵值，得到A1=4.2298、B1=7.1298、C1=8.3597三個閾值，將2012—2018年研究區(qū)域枯水期地下水位變化情況分為輕微上升、上升、明顯上升、劇烈上升4個等級，如圖5所示。

由圖5左圖可以看出，枯水期研究區(qū)地下水位變化程度由西向東逐漸增大，水位上升最大區(qū)域分布在5#、18#、20#、21#井周圍。利用由枯水期確定的地下水位變化范圍對應(yīng)的豐水期時間序列信息熵熵值A(chǔ)2=0.3245、B2=4.588、C2=6.1798，以及由自然間斷點分級法確定D=-2.9105，將2012—2018年研究區(qū)域豐水期地下水位變化情況分為明顯下降、輕微下降、輕微上升、上升、明顯上升、劇烈上升6個等級。由圖5的右圖可以看出，豐水期研究區(qū)13#—15#井周圍地下水位下降最明顯，18#、20#、21#井周圍水位上升最為明顯。在豐水期以地下水位下降為主的范圍內(nèi)，16#井地下水位上升，可能是因為16#井周圍開采井關(guān)停，導(dǎo)致其所在位置處地下水位明顯抬升。

3 地下水超采區(qū)治理措施

對于地下水超采區(qū)，建議采取以下措施進(jìn)行治理：

(1)地下水人工回灌，建立地下水水庫。為了解決鞍山西郊地下水超采問題，在運糧河及楊柳河沖洪積扇建設(shè)地下水回灌場，引太子河流域汛期洪水或大伙房輸水工程的水對地下水進(jìn)行回灌補給，增加地下水補給量，形成地下水庫，從而達(dá)到洪水資源化和實現(xiàn)超采區(qū)地下水資源新的采補平衡，以修復(fù)西郊超采區(qū)，防止生態(tài)環(huán)境地質(zhì)問題發(fā)生，并保證地下水資源的可持續(xù)利用。

(2)大力推進(jìn)工農(nóng)業(yè)節(jié)水建設(shè)，減少現(xiàn)有工農(nóng)業(yè)用水的地下水開采需求。對于不具備替代水源的地下水用水戶，如超采區(qū)周邊工業(yè)企業(yè)用水及干旱期農(nóng)業(yè)灌溉用水等，需大力推廣節(jié)水型企業(yè)建設(shè)，加強農(nóng)業(yè)節(jié)水建設(shè)，減少不必要的地下水開采，以促進(jìn)地下水超采區(qū)恢復(fù)。

(3)建立水權(quán)市場，引導(dǎo)水資源商品化。積極發(fā)展水權(quán)分配制度，實行總量控制與定額管理，將水權(quán)分配到鎮(zhèn)、村、農(nóng)戶，分配到每個用水企業(yè)。外來企業(yè)用水要在原用水戶手中購買用水權(quán)，原用水戶賣掉部分水權(quán)得到的資金，用來建節(jié)水工程。用經(jīng)濟(jì)杠桿調(diào)節(jié)生活、生產(chǎn)和生態(tài)等各方面的水量，促進(jìn)節(jié)水型社會建設(shè)，以減少超采區(qū)地下水分散性開采行為，促進(jìn)超采區(qū)快速恢復(fù)。

4 結(jié)論

(1)為保證鞍山市西郊超采區(qū)恢復(fù)治理、地下水資源合理開發(fā)利用及預(yù)防生態(tài)環(huán)境地質(zhì)問題發(fā)生，建議環(huán)境上限水位埋深不低于2m，工程上限水位埋深不低于5m，下限水位以超采區(qū)地下水位開始明顯回升的節(jié)點時間水位為準(zhǔn)，選擇2012年枯水期水位。

(2)建議在西郊增加地下水補給量，建立地下水庫，從而達(dá)到洪水資源化和實現(xiàn)超采區(qū)地下水資源新的采補平衡，以修復(fù)超采區(qū)，防止生態(tài)環(huán)境地質(zhì)問題發(fā)生，并保證地下水資源的可持續(xù)利用。

圖5 2012—2018年1月和8月地下水位變化等級分布圖

(3)后續(xù)研究應(yīng)重點分析鞍山西郊地下水變化的主因，為制定相對應(yīng)的治理措施提供依據(jù)。