關(guān)鍵詞：南水北調(diào)；地下水埋深；地下水開采；鄭州市

中圖分類號(hào)：Ｐ６４１；ＴＶ２１３．４ 文獻(xiàn)標(biāo)志碼：Ａ ｄｏｉ：１０．３９６９／ ｊ．ｉｓｓｎ．１０００－１３７９．２０２３．０５．０１６

引用格式：田大永，張春艷，程艷紅，等．南水北調(diào)中線工程通水后鄭州市地下水埋深變化研究［Ｊ］．人民黃河，２０２３，４５（５）：７６－８２，１０２．

地下水是寶貴的淡水資源，特別是在干旱半干旱地區(qū)，地下水是主要的甚至是唯一的可用水源［１－２］ 。

２０１４ 年１２ 月１２ 日，南水北調(diào)工程中線通水，重點(diǎn)解決河南、河北、北京、天津４ 個(gè)省（市）的水資源短缺問(wèn)題［３］ 。目前，南水北調(diào)對(duì)受水區(qū)地下水動(dòng)態(tài)的影響研究多集中在水質(zhì)方面，對(duì)地下水埋深的影響研究存在一定的不足［４－８］ 。劉昌明［９］ 討論了南水北調(diào)工程對(duì)華北平原地下水枯竭產(chǎn)生的潛在生態(tài)效應(yīng)，認(rèn)為南水北調(diào)工程可以實(shí)現(xiàn)對(duì)地下水的回補(bǔ)。劉予等［１０］ 分析了南水北調(diào)引水對(duì)北京市環(huán)境產(chǎn)生的影響，發(fā)現(xiàn)南水北調(diào)引水進(jìn)京后對(duì)北京地區(qū)生態(tài)環(huán)境有正面效應(yīng)，主要是增加水資源量和改善現(xiàn)有生態(tài)環(huán)境條件。馮忠倫等［１１］ 研究發(fā)現(xiàn)南水北調(diào)東線輸水后，對(duì)梁濟(jì)運(yùn)河區(qū)域地下水位的影響超過(guò)了降水的影響。Ｌｏｎｇ 等［７］ 研究發(fā)現(xiàn)南水北調(diào)中線工程對(duì)北京市地下水位的穩(wěn)定恢復(fù)有重要作用，對(duì)地下水恢復(fù)的貢獻(xiàn)為４０％。景兆凱等［１２］ 指出南水北調(diào)中線工程對(duì)河南省淺層地下水位恢復(fù)起到積極作用，主要體現(xiàn)在地下水埋深８ ～１２ ｍ區(qū)向埋深４～８ ｍ 區(qū)及４ ｍ 以下區(qū)轉(zhuǎn)變。Ｙａｎｇ 等［１３］ 系統(tǒng)分析了調(diào)水、用水和氣候?qū)θA北平原地下水儲(chǔ)量的影響，認(rèn)為在用水減少和調(diào)水、降水的共同影響下，華北平原未來(lái)地下水位可上升１５０ ｍｍ。綜合前人研究結(jié)果發(fā)現(xiàn)，南水北調(diào)對(duì)受水區(qū)地下水位回升、儲(chǔ)量增加產(chǎn)生了不同程度的積極影響，為地下水位穩(wěn)定恢復(fù)提供了有利條件。

鄭州市人均水資源占有量?jī)H為全國(guó)人均水平的１／１０，屬于資源型和水質(zhì)型缺水城市［１４－１５］ 。為順應(yīng)人口增長(zhǎng)和經(jīng)濟(jì)發(fā)展需要，鄭州市長(zhǎng)期開采地下水，導(dǎo)致地面沉降、塌陷等一系列問(wèn)題。２０１５ 年１ 月１２ 日，鄭州市劃定了承壓水限采區(qū)和禁采區(qū)，南水北調(diào)水源和黃河調(diào)水水源代替因禁限采地下水而減少的供水量。２０１７—２０１９ 年，南水北調(diào)中線工程向鄭州市年均調(diào)水３．６ 億ｍ３。然而，南水北調(diào)入鄭后，對(duì)鄭州市供水結(jié)構(gòu)及地下水位產(chǎn)生的影響，以往文獻(xiàn)中未見(jiàn)相關(guān)報(bào)道。因此，本文以鄭州市淺層和中深層地下水埋深為研究對(duì)象，應(yīng)用地統(tǒng)計(jì)學(xué)方法和主成分分析法，對(duì)南水北調(diào)中線工程通水后鄭州市城區(qū)供用水結(jié)構(gòu)變化、地下水埋深時(shí)空變化及影響因素進(jìn)行分析，以期為研究區(qū)水資源配置、地下水資源評(píng)價(jià)管理以及城市規(guī)劃等提供參考。

１研究方法

１．１地統(tǒng)計(jì)學(xué)方法

地統(tǒng)計(jì)學(xué)方法是分析地下水動(dòng)態(tài)時(shí)空變化的常用方法［１６］ ，其基本工具為半方差函數(shù)和克里金插值。半方差函數(shù)的表達(dá)式為

利用地統(tǒng)計(jì)學(xué)軟件ＧＳ＋，首先對(duì)地下水原始埋深數(shù)據(jù)進(jìn)行正態(tài)性檢驗(yàn)，發(fā)現(xiàn)原始數(shù)據(jù)并不服從正態(tài)分布；然后對(duì)原始數(shù)據(jù)進(jìn)行對(duì)數(shù)變換后再次進(jìn)行正態(tài)檢驗(yàn)，發(fā)現(xiàn)變換后的數(shù)據(jù)服從正態(tài)分布；最后，分別利用球狀模型、指數(shù)模型、高斯模型和線性模型計(jì)算半方差，經(jīng)分析比較，認(rèn)為球狀模型對(duì)淺層地下水埋深半變異函數(shù)進(jìn)行擬合的效果最好。

１．２主成分分析

法主成分分析法的步驟詳見(jiàn)文獻(xiàn)［２１－２２］。選取研究區(qū)２００３—２０１９ 年的地下水年開采量、年降水量、年蒸發(fā)量、黃河年調(diào)水量、南水北調(diào)年調(diào)水量、地下水年供水量、農(nóng)業(yè)年用水量、工業(yè)年用水量、生活年用水量和生態(tài)環(huán)境年用水量１０ 個(gè)因子進(jìn)行計(jì)算。利用ＳＰＳＳ軟件，首先對(duì)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化，發(fā)現(xiàn)各因子間存在強(qiáng)烈共線性，根據(jù)ＫＭＯ 檢驗(yàn)和Ｂａｒｔｌｅｔｔ 球形檢驗(yàn)結(jié)果，統(tǒng)計(jì)量ＫＭＯ_淺層＝ ０．６７１＞０．５、顯著性水平Ｐ_淺層＝ ０．００６ ＜０．０５，ＫＭＯ_中深層＝ ０．６８３＞０．５、著性水平Ｐ_中深層＝ ０．００６＜０．０５，結(jié)果表明可以進(jìn)行主成分降維分析；然后計(jì)算公因子方差，查看信息濃縮情況，根據(jù)總方差解釋結(jié)果和碎石圖，淺層和中深層地下水埋深影響因子各提取３個(gè)主成分，計(jì)算主成分系數(shù)，得到主成分表達(dá)式；最后進(jìn)行主成分評(píng)價(jià)。

２研究區(qū)概況

鄭州市地處河南中部，北臨黃河，西南部有黃土丘陵與嵩山相接，東部地勢(shì)平坦，屬黃河沖積平原，海拔８０～２００ ｍ，地勢(shì)由西南向東北傾斜。以鄭州市惠濟(jì)區(qū)、金水區(qū)、中原區(qū)、管城回族區(qū)、二七區(qū)為研究區(qū)，研究區(qū)內(nèi)有淺層地下水監(jiān)測(cè)井７６ 眼、中深層地下水監(jiān)測(cè)井６８ 眼，位置分布見(jiàn)圖１。

鄭州市屬溫帶半干旱大陸性氣候區(qū)，四季分明，７—９ 月降水量占全年降水量的５０％～６０％。以京廣鐵路為界，淺層含水層組大致分為東西兩部分：東部淺層含水層組位于黃河沖積平原，主要由第四系全新統(tǒng)和上更新統(tǒng)沖積成因的細(xì)中沙組成，局部為粉細(xì)沙；西部淺層含水層組位于塬前平原和臺(tái)塬崗地，主要由第四系全新統(tǒng)和上更新統(tǒng)沖洪積粉細(xì)沙層組成，其次由中更新統(tǒng)沖洪積成因的含黏粒、中細(xì)沙組成，臺(tái)塬處上部以中更新統(tǒng)沖洪積細(xì)中沙為主，下部為下更新統(tǒng)及新近系湖積相地層。淺層地下水為潛水或微承壓水類型。淺層地下水主要受降雨入滲補(bǔ)給，西南部為山前徑流補(bǔ)給，東北部為黃河側(cè)滲補(bǔ)給，其他地表水補(bǔ)給量較少。排泄方式主要受開采和氣象的影響，有越流排泄和少量蒸發(fā)排泄。淺層地下水天然流向整體由西南向東北，局部地區(qū)因開采而形成了橢圓形漏斗，地下水由周邊向中心匯流。淺層地下水埋深受開采、氣象等多重因素影響，２００６—２０１２ 年逐漸增大、２０１２年之后逐漸減?。ㄒ?jiàn)表１）。

根據(jù)埋藏條件和水文地質(zhì)特征，中深層地下水可分為上部亞層組及下部亞層組：上部亞層組主要由中更新統(tǒng)和下更新統(tǒng)組成，含水層巖性為中細(xì)沙、中粗沙及細(xì)沙，為承壓水，是地下水主要開采層段之一；下部亞層組由中粗沙、中細(xì)沙、細(xì)沙組成。中深層地下水受西南山區(qū)降水－巖溶裂隙水徑流和淺層水越流補(bǔ)給，排泄方式主要為開采和側(cè)向徑流排泄。中深層地下水整體流動(dòng)方向由西南向東北，局部因開采而形成了地下水降落漏斗。中深層地下水埋深主要受開采量的影響，局部地區(qū)受淺層地下水越流補(bǔ)給影響，２００５—２０１８年中深層地下水平均埋深與開采量變化趨勢(shì)大致相同（見(jiàn)表１）。

３數(shù)據(jù)來(lái)源

２０１５—２０１９年研究區(qū)地下水埋深數(shù)據(jù)來(lái)源于河南省地質(zhì)局生態(tài)環(huán)境地質(zhì)服務(wù)中心（其余埋深數(shù)據(jù)來(lái)源于前人研究報(bào)告），地下水開采量來(lái)源于鄭州市節(jié)約用水辦公室，氣象數(shù)據(jù)來(lái)源于鄭州市氣象局。

４供用水結(jié)構(gòu)分析

鄭州市年總用水量從２００３ 年的５．２６ 億ｍ^３ 增至２０１９ 年的９．４２ 億ｍ^３。南水北調(diào)中線工程通水前２００３—２０１４ 年年平均用水量５． ７９ 億ｍ^３， 通水后２０１５—２０１９ 年年平均用水量增至８．１３ 億ｍ^３，增加４０．４％。各類用水量及占比如圖２ 所示，農(nóng)業(yè)多年平均用水量占比由２００３—２０１４ 年的２０．７％降至２０１５—２０１９年的４．７％，生活多年平均用水量占比由２００３—２０１４ 年的４１．５％增至２０１５—２０１９ 年的４９．０％，生態(tài)環(huán)境多年平均用水量占比由２００３—２０１４ 年的１９．１％增至２０１５—２０１９ 年的２８．７％，工業(yè)用水量年際間變化幅度不大（占比２０％左右）。用水結(jié)構(gòu)的變化是城市化、社會(huì)經(jīng)濟(jì)結(jié)構(gòu)改變和水源替代政策等因素綜合作用的結(jié)果。

研究區(qū)供水來(lái)源主要有當(dāng)?shù)氐乇硭?、外調(diào)水（引黃水和南水北調(diào)水）、地下水和其他非常規(guī)水源。分析各類水資源年供水量及占比（見(jiàn)圖３）可知，南水北調(diào)中線工程通水前２００３—２０１４ 年地下水多年平均供水量為２．４０ 億～３．５８ 億ｍ^３（占總供水量的３８．０％ ～５９．８％，平均為４９．８％），通水后２０１５—２０１９ 年地下水多年平均供水量為０．５９ 億～３．０８ 億ｍ^３（占總供水量的６．３％～４０％，平均為１９．１％），減少的地下水供水量由調(diào)水量置換。調(diào)水量由２００３ 年的１．８１ 億ｍ^３（只有引黃水，占比３３．４％）增至２０１９ 年的６．３２ 億ｍ^３（引黃水和南水北調(diào)水，占比６７．１％），成為鄭州市主要的供水水源。２００３—２０１４ 年，地下水供水量占地下水可開采量的１３１．１％～３２２．３％（平均為２４２．７％），地下水超采嚴(yán)重；２０１５—２０１９ 年，地下水供水量占地下水可開采量的７５．２％～２２４．３％（平均為１２９．６％），下降１１３．１％。特別是２０１７—２０１９ 年，地下水開采量低于地下水可開采量，平均占比８２．９％，較南水北調(diào)中線工程通水前下降１５９．８％，為地下水穩(wěn)定恢復(fù)提供了有利條件。地下水開采量的大幅度減少是調(diào)水和限采綜合作用的結(jié)果，特別是南水北調(diào)輸水入鄭后（２０１５—２０１９ 年），地下水年均開采量比２００３—２０１４ 年年均開采量下降５８％。

５地下水埋深分析

５．１地下水埋深時(shí)空變化

２０１５ 年和２０１９ 年枯、豐水期，研究區(qū)淺層和中深層地下水埋深空間分布如圖４、圖５ 所示?？梢钥闯?，淺層地下水埋深由西部（大于４０ ｍ） 向東部（小于１０ ｍ）逐漸減小、中深層地下水埋深由西南（大于６０ ｍ）向東北（小于４０ ｍ）逐漸減小。

２０１５—２０１９ 年，淺層和中深層地下水埋深分區(qū)面積變化如圖６ 所示。與２０１５ 年相比，２０１９ 年淺層地下水埋深小于１０ ｍ 的分區(qū)面積在枯水期減?。常玻ァ⒇S水期增大３％；埋深為１０～２０、２０～３０、３０～４０ ｍ 的分區(qū)面積在枯水期分別增大２１％、２２％、１％，豐水期分別增大２％、１３％、３４％；埋深大于４０ ｍ 的分區(qū)面積在枯水期減?。担常ァ⒇S水期減?。福保?；不同年份地下水埋深為１０～２０ ｍ 的分區(qū)面積比例始終最大（占５０％左右）。與２０１５ 年相比，２０１９ 年中深層地下水埋深小于４０ ｍ的分區(qū)面積明顯增大，枯水期增大４６％、豐水期增大３９％；埋深為４０～６０ ｍ 的分區(qū)面積明顯減小，枯水期減?。常埃?、豐水期減?。玻罚ィ宦裆畲笥冢叮?ｍ 的分區(qū)面積略有減小，枯水期減小９％、豐水期減?。保福?。

研究區(qū)地下水埋深參數(shù)統(tǒng)計(jì)見(jiàn)表２。２０１９ 年，淺層地下水埋深為３．３６～４６．９５ ｍ，平均為１６．２４ ｍ；中深層地下水埋深為１４．８１～１１６．９３ ｍ，平均為４５．４５ ｍ。與２０１５ 年相比，２０１９ 年淺層地下水平均埋深在枯水期減?。埃叮?ｍ（減幅３． ７％）、豐水期減?。玻?３７ｍ （減幅１３．０％）；中深層地下水平均埋深在枯水期減?。担保?ｍ（減幅１０．２％）、豐水期減?。玻保埃?ｍ（減幅３１．６％）?？傮w而言，中深層地下水位恢復(fù)程度高于淺層地下水位的。研究區(qū)地下水埋深空間變異均屬中等變異，淺層變異系數(shù)隨時(shí)間減小，中深層變異系數(shù)隨時(shí)間增大。

５．２地下水埋深空間變異

地下水埋深半變異函數(shù)的參數(shù)值見(jiàn)表３。淺層地下水埋深的Ｃ_０ ／ （Ｃ_０ ＋Ｃ） 波動(dòng)較小，２０１５—２０１７ 年在０．００１～０．０１１范圍內(nèi)變化，２０１８—２０１９ 年增至０．０７５ ～０．１２８。Ｃ_０值呈增大趨勢(shì)，表明淺層地下水埋深有強(qiáng)烈的空間相關(guān)性，但隨著時(shí)間推移，空間相關(guān)性減弱，表明人類活動(dòng)等隨機(jī)因素引起的空間變異性增強(qiáng)，地形、水文地質(zhì)條件等結(jié)構(gòu)因素引起的空間變異性雖有所減弱但仍占主導(dǎo)地位。南水北調(diào)入鄭后，地下水開采量減少，淺層地下水除局部小部分區(qū)域埋深增大外，其余大部分區(qū)域埋深呈減小趨勢(shì)（見(jiàn)圖４），水位整體緩慢回升，地下水埋深變化空間響應(yīng)度較低，導(dǎo)致空間相關(guān)性雖有減弱但仍具有強(qiáng)烈相關(guān)性。

主成分Ｆ_ｓ１中，影響最大的４ 個(gè)因子分別為年開采量（－０．３９）、南水北調(diào)年調(diào)水量（０．３８）、生活年用水量（０．３９）和生態(tài)環(huán)境年用水量（０．３８）；主成分Ｆ_ｓ２ 中，影響最大的２ 個(gè)因子為黃河年調(diào)水量（－０．６５）和年降水量（０．３９）；主成分Ｆ_ｓ３中影響最大的２ 個(gè)因子為年降水量（０．６２）和工業(yè)年用水量（－０．６７）。影響中深層地下水埋深的因子歸納為調(diào)水量、開采量、用水量、降水量。計(jì)算主成分綜合得分并進(jìn)行降序排列，結(jié)果見(jiàn)表５。各因子綜合影響力最強(qiáng)的時(shí)段為２０１５—２０１９ 年，且影響力呈逐年增大趨勢(shì)，可見(jiàn)南水北調(diào)中線工程對(duì)中深層地下水埋深變化同樣產(chǎn)生了主導(dǎo)性的影響。

地下水埋深的影響因素分析結(jié)果表明，調(diào)水和減少開采量等對(duì)淺層和中深層地下水埋深產(chǎn)生了顯著的影響，且在南水北調(diào)中線工程通水之后影響力呈逐年增大趨勢(shì)，這與上文得出的人類活動(dòng)等隨機(jī)因素對(duì)地下水埋深空間變異性的影響越來(lái)越大的結(jié)論一致。

６結(jié)論

（１）２００３—２０１９ 年，鄭州市供水結(jié)構(gòu)發(fā)生了較大變化。２０１５ 年之前地下水超采嚴(yán)重，２０１５—２０１９ 年地下水供水量比２００３—２０１４ 年減少１１３．１％。２０１７—２０１９ 年地下水開采量低于地下水可開采量，較南水北調(diào)中線工程通水前減少１５９．８％。

（２）南水北調(diào)中線工程調(diào)水入鄭后，研究區(qū)淺層和中深層地下水多年平均埋深在枯水期和豐水期均有不同程度的減小。淺層地下水埋深小于１０ ｍ 的分區(qū)面積在枯水期明顯減小、豐水期略有增大，埋深為１０～２０、２０～３０、３０～４０ 的分區(qū)面積有不同程度的增大，埋深大于４０ ｍ 的分區(qū)面積明顯減?。恢猩顚拥叵滤裆钚∮冢矗?ｍ 的分區(qū)面積明顯增大，埋深為４０～６０ ｍ 的分區(qū)面積明顯減小，埋深大于６０ ｍ 的分區(qū)面積略有減小。

（ ３）地統(tǒng)計(jì)法分析結(jié)果表明，受南水北調(diào)中線工程的影響，鄭州市淺層和中深層地下水埋深的空間相關(guān)性均呈減弱趨勢(shì)。對(duì)于淺層地下水埋深，人類活動(dòng)等隨機(jī)因素引起的空間變異性越來(lái)越強(qiáng)，地形、水文地質(zhì)條件等結(jié)構(gòu)因素引起的空間變異性雖有所減弱但仍占主導(dǎo)地位；對(duì)于中深層地下水埋深，結(jié)構(gòu)因素和隨機(jī)因素引起的空間變異性逐漸趨于平衡。

（４）主成分分析結(jié)果表明，影響淺層和中深層地下水埋深變化的影響因子為調(diào)水量、開采量、用水量和降水量。南水北調(diào)中線工程通水后，水源替代政策對(duì)淺層和中深層地下水位恢復(fù)起到了主導(dǎo)性的積極作用。

【責(zé)任編輯 張華興】