石建省, 李國敏, 梁 杏, 陳宗宇, 邵景力, 宋獻方

1)中國地質(zhì)科學院水文地質(zhì)環(huán)境地質(zhì)研究所, 河北石家莊 050061; 2)中國科學院地質(zhì)與地球物理研究所, 北京 100029; 3)中國地質(zhì)大學(武漢), 湖北武漢 430074; 4)中國地質(zhì)大學(北京), 北京 100083; 5)中國科學院地理科學與資源研究所, 北京 100101

華北平原地下水演變機制與調(diào)控

石建省1), 李國敏2), 梁 杏3), 陳宗宇1), 邵景力4), 宋獻方5)

1)中國地質(zhì)科學院水文地質(zhì)環(huán)境地質(zhì)研究所, 河北石家莊 050061; 2)中國科學院地質(zhì)與地球物理研究所, 北京 100029; 3)中國地質(zhì)大學(武漢), 湖北武漢 430074; 4)中國地質(zhì)大學(北京), 北京 100083; 5)中國科學院地理科學與資源研究所, 北京 100101

自20世紀70年代以來, 華北平原長期的地下水過量開采形成了復雜的地下水降落漏斗, 并引發(fā)了一系列環(huán)境地質(zhì)問題。為建立區(qū)域社會經(jīng)濟發(fā)展與地下水資源與環(huán)境之間的協(xié)調(diào)機制, 國家 973項目“華北平原地下水演變機制與調(diào)控”在大量前人研究基礎(chǔ)上, 進一步分析了近50年來地下水動力場演變特征及其對人類活動的響應, 揭示了包氣帶結(jié)構(gòu)變化影響下包氣帶水力參數(shù)的變化, 研究了深層地下水開采引起的地面沉降、咸水移動和資源量轉(zhuǎn)換機理; 利用改進的地下水數(shù)值模擬技術(shù), 建立了2001—2010年華北平原地下水流數(shù)值模型, 計算了最新的地下水資源量。分析了地下水資源承載力的主要影響因素, 進行了華北平原地下水資源承載力定量評價, 給出了華北平原地下水的優(yōu)化調(diào)控途徑, 目前正在以高精度地下水流數(shù)值模型為核心, 模擬各區(qū)的地下水優(yōu)化調(diào)控方案。

華北平原; 地下水降落漏斗; 地下水資源管理

華北平原(圖1)是我國重要的政治、經(jīng)濟和農(nóng)業(yè)區(qū), 20世紀70年代以來, 華北平原地下水經(jīng)歷了長期的超采過程, 形成了大面積的地下水降落漏斗(張兆吉等, 2009a), 導致地下水資源衰減、地質(zhì)環(huán)境惡化(Liu et al., 2001), 直接威脅到該區(qū)水安全、糧食安全和重大工程設(shè)施安全(李國和等, 2007), 已成為經(jīng)濟社會發(fā)展的主要制約瓶頸(Foster et al., 2004)。未來社會經(jīng)濟的發(fā)展仍然面臨著兩難選擇: 一方面迫切需要提高地下水資源對社會、政治、和經(jīng)濟發(fā)展的支撐能力; 另一方面, 和諧社會建設(shè)又亟需有效地調(diào)控地下水的開發(fā)利用, 解決與地下水相關(guān)的地質(zhì)環(huán)境問題(Liu et al., 2008), 因此, 華北平原很需要探索一條地下水可持續(xù)開發(fā)利用的科學路徑(Xia et al., 2007)。

在世界范圍內(nèi), 像華北平原這樣的大型地下水系統(tǒng)均在當?shù)氐纳鐣?jīng)濟中起著重要作用, 而強烈開采引起的水位下降也往往是這些地下水系統(tǒng)面臨的共同問題。如美國西部的高平原地下水系統(tǒng)(Sophocleous, 2005), 澳大利亞的大自流盆地(Habermehl, 2006)等。但是, 華北平原地下水系統(tǒng)比上述系統(tǒng)更為復雜, 經(jīng)受的開采強度更大, 引發(fā)的環(huán)境問題也更為嚴重。這給華北平原地下水可持續(xù)利用問題的解決提出了很大的挑戰(zhàn)。

為此, 國家重點基礎(chǔ)研究發(fā)展計劃(“973”)項目《華北平原地下水演變機制與調(diào)控》聚焦于二個亟需解決的關(guān)鍵科學問題: ①地下水系統(tǒng)結(jié)構(gòu)變異與地下水資源演變機制; ②地下水資源開發(fā)-環(huán)境變化-經(jīng)濟發(fā)展協(xié)調(diào)機制, 通過開展四項內(nèi)容的研究: (a)人類活動條件下區(qū)域地下水系統(tǒng)響應; (b)含水層系統(tǒng)結(jié)構(gòu)變異與地下水可利用資源變化機理; (c)地下水-環(huán)境-社會經(jīng)濟系統(tǒng)耦合機制與評價體系; (d)環(huán)境、經(jīng)濟約束下的地下水調(diào)控來進一步揭示華北平原地下水復合漏斗形成演變機理, 來闡明人類活動對地下水資源數(shù)量與質(zhì)量的影響機制, 揭示地下水資源開發(fā)-環(huán)境變化-經(jīng)濟發(fā)展協(xié)調(diào)機制, 發(fā)展和完善人類活動影響下地下水系統(tǒng)演變及相關(guān)環(huán)境變化的基礎(chǔ)理論。最終建立地下水系統(tǒng)危機臨界識別指標體系和地下水資源承載力評價體系, 提出“地下水復合漏斗”控制與修復理論方法(中國地質(zhì)科學院, 2014)。項目研究思路見圖2。

圖1 華北平原平面位置圖Fig. 1 Location of the North China Plain

經(jīng)過四年多的研究, 項目組已經(jīng)取得了一系列重要進展, 詳情如下。

1 地下水系統(tǒng)結(jié)構(gòu)變異與地下水資源演變機制

1.1 人類活動條件下區(qū)域地下水系統(tǒng)的演變

通過華北平原1960年、1980年和2000年淺層地下水資源評價及均衡模式的對比(張兆吉等, 2009a),建立了大規(guī)模開采前和大規(guī)模開采后的水均衡模式(圖3)。查明了區(qū)域地下水系統(tǒng)的演變特征如下。

大規(guī)模開采前(1960年左右), 淺層地下水位較淺, 降水入滲補給系數(shù)相對較小, 地表水系發(fā)達,地下水蒸發(fā)強烈, 中東部深層地下水向淺層地下水越流補給; 大規(guī)模開采后(2000年), 深、淺層地下水位大幅度下降(圖4), 形成復合降落漏斗。深層地下水不再向淺層地下水補給, 反而接受淺層地下水的越流補給(Shao et al., 2013)。

地下水的水化學組分, 在大規(guī)模開采前有明顯的規(guī)律性: 山前平原以HCO3型水為主, 中部平原則為 SO4、Cl、HCO3混合型水, 濱海平原則以 Cl及Cl-SO4型水為主(張宗祜等, 1997); 在大規(guī)模開采后,水循環(huán)條件的改變導致淺層地下水化學特征趨于復雜, 尤其是中部平原, 由早期的HCO3-SO4、Cl-SO4、SO4-Cl混合型水演化為 HCO3-SO4、HCO3-Cl、Cl-SO4、SO4-Cl、SO4、Cl-SO4、SO4-Cl等混合型水,潛水Cl及Cl-SO4型水分布面積大大減小, 地下水趨于淡化, 水化學類型轉(zhuǎn)化為HCO3或SO4型(張兆吉等, 2009a)。但總體上, 地下水無機化學組分含量從山前至濱海逐漸升高, 水質(zhì)逐漸變差的趨勢仍然沒有改變(費宇紅等, 2014)。

圖2 項目研究思路框圖Fig. 2 Structure of the study logic of this project

圖3 華北平原地下水循環(huán)模式對比(a-1960年; b-2000年)Fig. 3 Comparison of groundwater circulations in the North China Plain (a-1960; b-2000)

1.2 包氣帶結(jié)構(gòu)變異與地下水補給演變

華北平原長期超采地下水導致水位持續(xù)下降,飽水帶向包氣帶轉(zhuǎn)化, 形成厚度大、非均質(zhì)性更為復雜的包氣帶; 同時, 土壤的釋水改變原有的應力狀態(tài), 土壤物理結(jié)構(gòu)、土壤水力參數(shù)發(fā)生了改變,從而影響降雨或灌溉水的入滲補給(Liang et al., 2010)。本項目利用大規(guī)模示蹤試驗, 獲得了一系列入滲補系數(shù), 經(jīng)計算, 山前平原平均補給量13.7～564.3 mm/a , 平均補給系數(shù)為0.17; 中部平原年均補給量 11.9～191.7 mm/a, 平均補給系數(shù)為0.097(Wang et al., 2011)。灌溉是產(chǎn)生垂向補給的主要驅(qū)動力, 灌溉量越大, 對地下水的補給量也越大(Lu et al., 2011)。隨埋深的增大, 補給系數(shù)與地下水位埋深的關(guān)系不再明顯。華北平原平均補給系數(shù)為0.9, 實驗結(jié)果為0.16。近50年來, 根區(qū)土層滲透系數(shù)的變異性對潛在補給量的影響較大, 根區(qū)以下土層(2 m以下)滲透系數(shù)的變異性對潛在補給量影響較小。當包氣帶厚度小于5 m時, 實際補給過程與降雨事件響應明顯。入滲補給隨地下水位埋深的增大, 補給系數(shù)與地下水位埋深的關(guān)系不再明顯(Huo et al., 2014)。在地下水埋深大于土壤蒸發(fā)極限埋深的情況下, 巖性非均質(zhì)對多年平均補給量無顯著影響(Lu et al., 2011; Zha et al., 2013; 霍思遠等, 2013)。示蹤劑方法是確定深厚包氣帶區(qū)地下水補給的可靠技術(shù), 利用該技術(shù)研究了華北平原淺層地下水的影響因素、補給過程和分布特征。根據(jù)示蹤劑的實驗成果和理論研究, 發(fā)現(xiàn)在地下水的深埋區(qū), 包氣帶厚度的變化對地下水垂直補給量影響不大, 地下水補給量主要決定于淺層水分平衡過程和土壤特征,包氣帶厚度增大將使得地下水補給過程產(chǎn)生滯后并變化平緩(譚秀翠等, 2012; Huo et al., 2014)。

圖4 1984—2003年華北平原淺層地下水水位變差圖(張兆吉等, 2009b)Fig. 4 Shallow groundwater level variation graph in the North China Plain from 1984 to 2003(ZHANG et al., 2009b)

1.3 弱透水層壓縮與地面沉降機理

地下水的嚴重超采使華北平原的北京、天津、衡水及滄州等地都出現(xiàn)了嚴重的地面沉降現(xiàn)象(何慶成等, 2006)。同時工程建設(shè)也成為影響地面沉降的一個重要人為因素(龔士良, 2008; 楊少華, 2009),認識和區(qū)分這兩種作用對地面沉降發(fā)生的影響, 對評估地下水開采引起的地面沉降具有重要意義。為此, 我們以天津濱海新區(qū)為典型研究區(qū), 建立了三維地下水滲流場與應力場耦合模型, 對高層荷載和地下水開采雙重作用下地面沉降的演化規(guī)律進行了模擬, 識別出濱海新區(qū)地面沉降演變的時空特征, 即: 濱海新區(qū)地面沉降的發(fā)生與發(fā)展與地下水開采關(guān)系非常密切。持續(xù)開采一年后井300 m以內(nèi)為劇烈沉降區(qū), 距井口超過600 m, 沉降趨于平緩。高層建筑荷載引起的地面沉降范圍很小, 局限于高層建筑物 200 m范圍, 并且與建筑物的高度有關(guān)(Ma et al., 2011)。地面沉降過程中, 在固結(jié)初始階段, 非結(jié)合水首先滲出, 隨著固結(jié)時間的延續(xù), 有效應力相應增高, 當有效應力超過微孔隙層次的結(jié)構(gòu)強度時,孔隙瞬間被壓縮, 孔隙水壓力驟然升高, 弱結(jié)合水突然釋出, 致使固結(jié)速率加快, 土體結(jié)構(gòu)經(jīng)過調(diào)整后, 完成黏性土釋水。

1.4 弱透水層在咸水越流中的作用

華北平原第四系含水層具有上咸下淡的結(jié)構(gòu),長期開采深層淡水, 導致淺層咸水越流量增加, 但目前并未發(fā)現(xiàn)大面積的深層地下淡水礦化度明顯增加的現(xiàn)象, 水文地質(zhì)學家們提出了弱透水層“半透膜效應”的阻鹽作用(王家兵, 2002)。本次研究通過采用密閉的原狀土入滲裝置開展了不同密度的原狀粘土越流試驗得到如下認識: (1)越流過程中粘性土弱透水層存在半透膜超濾鹽分的現(xiàn)象; (2)超濾鹽分的容量隨著弱透水層的持續(xù)壓縮而減小; (3)當超濾容量達到某一臨界值時, 發(fā)生咸水穿透。上述認識說明: 隨著華北平原水位持續(xù)下降, 地面沉降持續(xù)發(fā)展, 粘性土弱透水層的半透膜超濾鹽分的能力將逐漸降低, 最終將有可能發(fā)生咸水突破而影響深層承壓水的開發(fā)利用。

1.5 深層地下水漏斗演變與更新性

以滄州為主要研究區(qū), 詳細研究了深層地下水漏斗的演變過程: 1964年以前地下水開采量很小, 地下水流場基本處于天然狀態(tài), 1965年以后, 隨著地下水開采量的增加, 地下水降落漏斗的擴大速度逐漸加快, 到 2000年左右達到頂峰, 隨后地下水開采量趨于穩(wěn)定或減少, 漏斗水位逐漸回升(韓占濤等, 2013)。

以地面沉降嚴重并且監(jiān)測比較完善的天津市平原區(qū)作為典型區(qū), 建立該地區(qū)地下水流-地面沉降數(shù)值模擬模型(李文運等, 2012; Cui et al., 2014), 模擬區(qū)面積1.1×104km2, 垂向上分為6個含水層組, 網(wǎng)格剖分大小為500 m×500 m。通過11年的地下水位、地面沉降監(jiān)測數(shù)據(jù)識別模型, 運用模型分析人類大規(guī)模開采深層水條件下地下水資源量的組成及其變化。結(jié)果表明: 上部潛水通過弱透水層對深層地下水補給, 占開采量的40%; 壓縮釋水量約占開采量的比例平均為33%, 隨著抽水時間的延長, 壓縮釋水量占總開采量的比例略微減小; 通過外圍邊界側(cè)向流入天津市深層含水組的地下水可以占到開采量的24%。

1.6 華北平原地下水資源評價

在前期工作的基礎(chǔ)上(Wang et al., 2008; Shao et al., 2013), 運用地下水模擬軟件MODFLOEW建立了長序列華北平原地下水流數(shù)值模擬模型, 結(jié)合地下水模擬的并行算法(黃林顯等, 2012; 程湯培等, 2013; Cheng et al., 2014)、區(qū)域地下水開采估算(李玲等, 2013)等研究。以該模型為主要手段, 分析了華北平原地下水均衡、評價了地下水資源量。結(jié)果表明, 2001—2010年華北平原地下水淺層年均總補給量為 243×108m3/a, 總排泄量為 262×108m3/a,目前人類活動條件下, 地下水開采量是主要排泄項, 占總排泄量的 66.0%。多年地下水補排差為-19×108m3/a, 華北平原淺層水總體處于超采狀態(tài)。深層水的補給來自潛水含水層的越流和側(cè)向流入, 為25.5×108m3/a, 排泄量為43.1×108m3/a, 其中地下水開采量為42.7×108m3/a, 是深層水的主要排泄項。地下水補排差為-17. 6×108m3/a, 處于嚴重超采狀態(tài)。此外, 本研究還分行政區(qū)、水文地質(zhì)單元、礦化度對地下水補給資源量和可開采資源量進行了評價, 為地下水資源調(diào)控提供了重要基礎(chǔ)數(shù)據(jù)。

2 地下水資源開發(fā)與環(huán)境、經(jīng)濟協(xié)調(diào)機制

2.1 地下水資源開發(fā)及其資源環(huán)境效應

長期、大規(guī)模的地下水開采不但使地下水流場發(fā)生了劇烈改變, 而且引起了一系列環(huán)境效應, 如濕地面積大幅度減小(喬光英等, 2010), 嚴重的地面沉降和地裂縫(何慶成等, 2006; 李世雄等, 2006),嚴重的海水入侵等(左文喆, 2006)。

為此, 我們抓住華北平原的主要耗水因素, 分析了農(nóng)業(yè)生態(tài)用水對地下水位變化的影響, 闡明了主要農(nóng)作物耗水規(guī)律, 并提出了最佳灌溉模式。經(jīng)計算, 最佳灌溉模式可減少22%左右的農(nóng)業(yè)用水量。

為建立地下水開發(fā)與經(jīng)濟發(fā)展相互關(guān)系, 我們分析了人口、土地利用和工農(nóng)業(yè)發(fā)展歷程與地下水開發(fā)利用數(shù)據(jù), 闡明了地下水開發(fā)利用與近50年社會經(jīng)濟發(fā)展的互動關(guān)系。我們還系統(tǒng)分析了環(huán)境問題(海水入侵、土壤鹽漬化和地下水污染)與地下水的關(guān)系, 闡明了再生水利用對地下水環(huán)境影響程度。

2.2 華北平原地下水資源承載力及可調(diào)控途徑

2.2.1 地下水資源承載力概念

為實現(xiàn)地下水資源的可持續(xù)利用, 應當確定合理的地下水資源承載力。前人根據(jù)不同地區(qū)不同的約束條件, 提出了不同的地下水資源承載力概念和評價方法(張鑫等, 2001; 王順久等, 2004; 萬星等, 2006; 董克剛等, 2008)。然而, 針對華北平原特有的社會經(jīng)濟結(jié)構(gòu)和水文地質(zhì)條件, 前人提出的地下水資源承載力概念和方法并不太適用。

為實現(xiàn)華北平原地下水資源的可持續(xù)利用, 我們提出了綜合考慮地下水資源的自然屬性、環(huán)境屬性和社會屬性的地下水資源承載力的概念。地下水資源承載力即在一定發(fā)展階段下, 以可預見的技術(shù)、經(jīng)濟和社會發(fā)展水平為依據(jù), 以可持續(xù)發(fā)展為原則, 以地下水最大可利用量為前提, 地下水資源對區(qū)域社會經(jīng)濟發(fā)展的最大支撐能力。地下水資源承載力具有動態(tài)性、可調(diào)控性和有限性。

2.2.2 華北平原地下水資源承載力評價

(1)地下水資源承載力評價方法

將地下水資源承載力簡化為地下水最大可利用量和用水效率的函數(shù), 建立了地下水資源承載力靜態(tài)評價方法, 提出了以地下水資源開采程度(RQ)和社會經(jīng)濟發(fā)展程度(RG)關(guān)系為核心、以區(qū)域最優(yōu)用水效率為統(tǒng)一標準的多地區(qū)地下水資源利用與社會經(jīng)濟發(fā)展狀況動態(tài)評價方法。

圖5 現(xiàn)狀條件下(2000年)華北平原各地市地下水對社會經(jīng)濟承載能力Fig. 5 Groundwater resources carrying capacities of cities in the North China Plain (2000)

圖6 華北平原各地市RG-RQ關(guān)系圖(2000年)Fig. 6 The relationship between the extent of groundwater exploitation (RQ) and the economic development (RG) in the North China Plain (2000)

(2)華北平原地下水資源承載力靜態(tài)評價

根據(jù) 2000年的社會經(jīng)濟發(fā)展狀況及耗水量統(tǒng)計數(shù)據(jù)(包括產(chǎn)業(yè)結(jié)構(gòu)、用水結(jié)構(gòu)、耗水水平、地下水在不同行業(yè)之間的分配比例等), 華北平原各地市地下水可安全承載的 GDP規(guī)模以北京市最大, 為1145億元, 東營市最小, 不足3億元。其由大到小的排序見圖5。

(3)華北平原地下水資源利用與社會經(jīng)濟發(fā)展狀況評價

圖 6為2000年華北平原各地市地下水開采程度(RQ)與地下水支撐的經(jīng)濟發(fā)展程度(RG)關(guān)系圖。德州、保定、聊城、濮陽、新鄉(xiāng)等位于 I區(qū), 說明該五個地區(qū)地下水開發(fā)利用程度較低, 經(jīng)濟發(fā)展基本不受水資源制約; 焦作、鶴壁、廊坊、衡水、邢臺、安陽、秦皇島、邯鄲、濱州等位于II區(qū), 表明該九個地區(qū)已發(fā)生地下水超采, 經(jīng)濟發(fā)展實際并未超載, 經(jīng)濟發(fā)展受到水資源制約; 石家莊、滄州、唐山、北京位于III區(qū), 表明該四個地區(qū)不僅發(fā)生了地下水超采, 而且也出現(xiàn)了經(jīng)濟超載現(xiàn)象, 經(jīng)濟發(fā)展受到水資源的嚴重制約。

2.2.3 提高地下水資源承載力的可能途徑

由前所述, 地下水資源的承載力是地下水資源可利用量和用水效率的函數(shù), 因此, 提高地下水資源承載力的途徑從兩方面考慮: 提高地下水可利用量和提高用水效率, 而后者的效果更顯著。對于華北平原而言, 提高地下水可利用量的主要途徑包括:山前人工雨洪調(diào)蓄、中東部平原發(fā)展淺層弱滲透含水層淡水開采技術(shù)和微咸水改造利用技術(shù)等; 提高水資源利用效率的途徑包括: 發(fā)展農(nóng)業(yè)節(jié)水、提高工業(yè)用水重復利用率、調(diào)整產(chǎn)業(yè)結(jié)構(gòu)、改變經(jīng)濟增長模式、提高公眾節(jié)水意識等。

3 環(huán)境、經(jīng)濟約束下的地下水調(diào)控

3.1 華北平原地下水臨界水位

為合理開發(fā)利用、有效保護地下水資源, 在某一地區(qū), 對于特定的地下水開采目的層確定合理的地下水開采控制水位, 即為臨界水位, 以防止由于過量開采地下水對地下水流系統(tǒng)自身及周邊環(huán)境引發(fā)破壞作用。

利用中國地質(zhì)環(huán)境監(jiān)測院華北平原地區(qū)國家級和省級近 20年的水位長序列監(jiān)測數(shù)據(jù), 分析確定了研究區(qū)淺層和深層地下水的臨界水位(圖7, 8)。

圖7 華北平原淺層地下水臨界水位Fig. 7 The critical groundwater level of shallow aquifer in the North China Plain

圖8 華北平原深層地下水臨界水位Fig. 8 The critical groundwater level of deep aquifer in the North China Plain

3.2 地下水調(diào)控

基于上述研究成果, 我們正在以新建立的地下水數(shù)值模型為核心, 以確定的各區(qū)地下水位約束值為約束條件, 考慮《海河流域綜合規(guī)劃(2011—2030年》及南水北調(diào)輸入水量的變動, 以及各種節(jié)水、優(yōu)化開采和回補方案, 模擬和預測 2011—2015年, 2016—2020年, 2021—2030年, 2031—2050年四個時段的地下水開采方案、所能提供的水資源量, 以及地下水流場的變化情況, 并將在此基礎(chǔ)上提出各區(qū)合適的地下水開采量及開采方案。

4 結(jié)論

華北平原地下水系統(tǒng)在經(jīng)歷 40余年的大規(guī)模開采后, 已經(jīng)發(fā)生了劇烈變化。深、淺層地下水均由于過量開采而形成大面積降落漏斗。這一改變伴隨著地下水流系統(tǒng)參數(shù)和水質(zhì)的變化, 以及由此引起的一系列環(huán)境問題。淺層地下水的補給由原來受地下水位影響而變成較穩(wěn)定的補給。深層地下水位下降, 引起淺層地下水向下越流補給, 咸水下移,弱透水層壓縮和地面沉降, 以及海水入侵。為保證華北平原地下水可持續(xù)利用, 建立了高精度地下水流數(shù)值模型, 重新評價了地下水資源, 提出了適用于華北平原的地下水資源承載力評價方法, 正模擬地下水資源優(yōu)化調(diào)控方案, 尋求可持續(xù)發(fā)展的地下水資源利用途徑。

程湯培, 莫則堯, 邵景力. 2013. 基于JASMIN的地下水流大規(guī)模并行數(shù)值模擬[J]. 計算物理, 30(3): 317-325.

董克剛, 王威, 于強, 易長榮, 周俊. 2008. 海河流域沉降區(qū)地下水資源承載力評價指標體系[J]. 地下水, 30(4): 12-15.

龔士良. 2008. 上海地面沉降影響因素綜合分析與地面沉降系統(tǒng)調(diào)控對策研究[D]. 上海: 華東師范大學.

費宇紅, 張兆吉, 郭春艷, 王春曉, 雷廷, 劉瑾. 2014. 區(qū)域地下水質(zhì)量評價及影響因素識別方法研究-以華北平原為例[J].地球?qū)W報, 35(2): 131-138.

韓占濤, 王平, 張威, 夏文則, 馬麗莎, 李亞松, 張發(fā)旺. 2013.滄州地區(qū)深層地下水漏斗演化與供水策略分析[J]. 水文地質(zhì)工程地質(zhì), 40(5): 29-33.

何慶成, 劉文波, 李志明. 2006. 華北平原地面沉降調(diào)查與監(jiān)測[J]. 高校地質(zhì)學報, 12(2): 195-209.

黃林顯, 邵景力, 崔亞莉. 2012. 大區(qū)域地下水模擬的預優(yōu)并行GMRES(m)算法研究[J]. 現(xiàn)代地質(zhì), 6(3): 614-619.

霍思遠, 靳孟貴, 梁杏. 2013. 包氣帶弱滲透性黏土透鏡體對降雨入滲補給影響的數(shù)值模擬[J]. 吉林大學學報(地球科學版), 43(5): 1579-1587.

李國和, 許再良, 孫樹禮, 荊志東. 2007. 華北平原地面沉降對高速鐵路的影響及其對策[J]. 鐵道工程學報, (8): 7-12.

李玲, 邵景力, 崔亞莉. 2013. 華北平原區(qū)地下水開采量估算研究[J]. 南水北調(diào)與水利科技, 11(5): 81-85.

李世雄, 李守定, 郜洪強. 2006. 河北平原地裂縫分布特征及成因機制研究[J]. 工程地質(zhì)學報, 14(2): 178-183.

李文運, 崔亞莉, 蘇晨, 張偉, 邵景力. 2012. 天津市地下水流-地面沉降耦合模型[J]. 吉林大學學報(地球科學版), 42(3): 805-813.

喬光建, 孫梅英, 王斌. 2010. 河北省平原濕地減少原因分析[J].水資源保護, 26(3): 33-37.

譚秀翠, 楊金忠. 2012. 石津灌區(qū)地下水潛在補給量時空分布及影響因素分析[J]. 水利學報, 43(2): 143-152.

萬星, 丁晶, 張曉麗. 2006. 區(qū)域地下水資源承載力綜合評價的集對分析方法[J]. 城市環(huán)境與城市生態(tài), 19(2): 8-10.

王家兵. 2002. 華北平原深層淡水在開采條件下接受上覆咸水越流補給——以天津平原為例[J]. 水文地質(zhì)工程地質(zhì), (6): 35-37.

王順久, 楊志峰, 丁晶. 2004. 關(guān)中平原地下水資源承載力綜合評價的投影尋蹤方法[J]. 資源科學, 26(6): 106-110.

楊少華. 2009. 城市建設(shè)引起地面沉降數(shù)值模擬研究[D]. 北京:中國地質(zhì)大學.

張鑫, 王紀科, 蔡煥杰, 王正興. 2001. 區(qū)域地下水資源承載力綜合評價研究[J]. 水土保持通報, 21(3): 24-27.

張兆吉, 費宇紅, 陳宗宇, 趙宗壯, 謝振華, 王亞斌, 苗普祥,楊麗芝, 邵景力, 靳孟貴, 許廣明, 楊齊青. 2009a. 華北平原地下水可持續(xù)利用調(diào)查評價[M]. 北京: 地質(zhì)出版社.

張兆吉, 費宇紅. 2009b. 華北平原地下水可持續(xù)利用圖集[M].北京: 中國地圖出版社.

張宗祜, 施德鴻, 沈照理, 鐘佐燊, 薛禹群. 1997. 人類活動影響下華北平原地下水環(huán)境的演化與發(fā)展[J]. 地球?qū)W報18(4): 337-344.

中國地質(zhì)科學院. 2014. 中國地質(zhì)科學院2013年度十大科技進展揭曉[J]. 地球?qū)W報, 35(1): 1-5.

左文喆. 2006. 秦皇島洋戴河平原海水入侵調(diào)查與研究[D]. 北京: 中國地質(zhì)大學.

References:

CHENG Tang-pei, MO Ze-yao, SHAO Jing-li. 2013. JASMIN-based Massive Parallel Computing of Large Scale Groundwater Flow[J]. Chinese Journal of Computational Physics, 30(3): 317-325(in Chinese with English abstract).

CHENG Tang-pei, MO Ze-yao, SHAO Jing-li. 2014. Accelerating Groundwater Flow Simulation in MODFLOW Using JASMIN-Based Parallel Computing[J]. Ground Water, 52(2): 194-205.

Chinese Academy of Geological Sciences. 2014. Top Ten Scientific and Technological Progress of Chinese Academy of Geological Sciences in the Year 2013 Unveiled[J]. Acta Geoscientica Sinica, 35(1): 1-5(in Chinese with English abstract).

CUI Ya-li, SU Chen, SHAO Jing-li, WANG Ya-bin, CAO Xiao-yuan. 2014. Development and application of a coupled groundwater flow and land subsidence model for the plain region of Tianjin[J]. Journal of Earth Science(in press).

DONG Ke-gang, WANG Wei, YU Qiang, YI Chang-rong, ZHOU Jun. 2008. Assessment Index System ofGroundwater Resources Carrying Capacity in the Land- subsidence Area of-Haihe Basin[J]. Ground water, 30(4): 12-15(in Chinese with English abstract).

FOSTER S, GARDUNO H, EVANS R, OLSON D, TIAN Yuan, ZHANG Wei-zhen, HAN Zai-sheng. 2004. Quaternary aquifer of the North China Plain—assessing and achieving groundwater resource sustainability[J]. Hydrogeology Journal, 12(1): 81-93.

GONG Tu-liang. 2008. The Study on Acting Factors and Systemic Control of Land Subsidence in Shanghai[D]. Shanghai: East China Normal University(in Chinese with English abstract).

HABERMEHL M. 2006. The great artesian basin, Australia[C]//Into the well from which you drink do not throw stones: 82.

FEI Yu-hong, ZHANG Zhao-ji, GUO Chun-yan, WANG Chun-xiao, LEI Ting, LIU Jin. 2014. Research on the Method for Evaluation and Influence Factors Identification of Regional GroundWater Quality: A Case Study of the North China Plain[J]. Acta Geoscientica Sinica, 35(2): 131-138(in Chinese with English abstract).

HAN Zhan-tao, WANG Ping, ZHANG Wei, XIA Wen-ze, MA Li-sha, LI Ya-song, ZHANG Fa-wang. 2013. Analysis of the evolvement of deep confined groundwater depression cone and water supply strategy for Cangzhou area[J]. Hydrogeology ＆ Engineering Geology, 40(5): 29-33(in Chinese with English abstract).

HE Qing-cheng, LIU Wen-bo, LI Zhi-ming. 2006. Land Subsi-dence Survey andMonitoring in the North China Plain[J]. Geological Journal of China Universities, 12(2): 195-209(in Chinese with English abstract).

HUANG Lin-xian, SHAO Jing-li, CUI Ya-li, CHENG Tang-pei, LI Ling. 2012. Study on Parallel GMRES(m) Algorithm Integrating with Preconditioners Applied to Large-scale Groundwater Flow Simulation[J]. Geoscience, 6(3): 614-619(in Chinese with English abstract).

HUO Si-yuan, JIN Meng-gui, LIANG Xing, LIN Dan. 2014. Changes of vertical groundwater recharge with increase in thickness of vadose zone simulated by one-dimensional variably saturated flow model[J]. Journal of China University of Geosciences(in press).

HUO Si-yuan, JIN Meng-gui, LIANG Xing. 2013. Impacts of Low-Permeability Clay Lensin Vadose Zone onto Rainfall Infiltration and Ground water Recharge Using Numerical Simulation of Variably Saturated Flow[J]. Journal of University (Earth Science Edition), 43(5): 1579-1587(in Chinese with English abstract).

LI Guo-he, XU Zai-liang, SUN Shu-li, JING Zhi-dong. 2007. The Influence of Surface Subsidence on Construction of Highspeed Railway in North China Plain and Its Countermeasures[J]. Journal of Railway Engineering Society, (8): 7-12(in Chinese with English abstract).

LI Ling, SHAO Jing-li, CUI Ya-li. 2013. Estimation of Groundwater Withdrawal in the North China Plain[J]. South-to-North Water Transfers and Water Science & Technology, 11(5): 81-85(in Chinese with English abstract).

LI Shi-xiong, LI Duo-ding, HAO Hong-qiang. 2006. The Distribution Characters and Origin Mechanics of Ground Fissureshazard in Hebei Plain[J]. Journal of Engineering Geology, 14(2): 178-183(in Chinese with English abstract).

LI Wen-yun, CUI Ya-li, SU Chen, ZHANG Wei, SHAO Jing-li. 2012. An Integrated Numerical Groundwater and Land Subsidence Model of Tianjin[J]. Journal of University(Earth Science Edition), 42(3): 805-813(in Chinese with English abstract).

LIANG Xing, LIU Yu, JIN Meng-gui, LU Xin-chen, ZHANG Ren-qun. 2010. Direct observation of complex Tóthian groundwater flow systems in the laboratory[J]. Hydrological Processes, 24(24): 3568-3573.

LIU Chang-ming, YU Jing-jie, KENDY E. 2001. Groundwater exploitation and its impact on the environment in the North China Plain[J]. Water International, 26(2): 265-272.

LIU Jie, ZHENG Chun-miao, ZHENG Li, LEI Yu-ping. 2008. Ground Water Sustainability: Methodology and Application to the North China Plain[J]. Ground Water, 46(6): 897-909.

LU Xiao-hui, JIN Meng-gui, VAN GENUCHTEN M T, WANG Bing-guo. 2011. Groundwater Recharge at Five Representative Sites in the Hebei Plain, China[J]. Ground Water, 49(2): 286-294.

MA Feng-shan, WEI Ai-hua, HAN Zhan-tong, ZHAO Hai-jun, GUO Jie. 2011. The Characteristics and Causes of Land Subsidence in Tanggu Based on the GPS Survey System and Numerical Simulation[J]. Acta Geologica Sinica, 85(6): 1495-1507.

QIAO Guang-jian, SUN Mei-ying, WANG Bin. 2010. Analysis of cause of sharp decrease of wetlands in Hebei Plain[J]. Water Resources Protection, 26(3): 33-37(in Chinese with English abstract).

SHAO Jing-li, LI Ling, CUI Ya-li, ZHANG Zhao-ji. 2013. Groundwater Flow Simulation and its Application in Groundwater Resource Evaluation in the North China Plain, China[J]. Acta Geologica Sinica(English Edition), 87(1): 243-253.

SOPHOCLEOUS M. 2005. Groundwater recharge and sustainability in the High Plains aquifer in Kansas, USA[J]. Hydrogeology Journal, 13: 351-365.

TAN Xiu-cui, YANG Jin-zhong. 2012. Temporal and spatial distribution of the potential recharge and influencing factors in Shijin Irrigation District[J]. Journal of Hydraulic Engineering, 43(2): 143-152(in Chinese with English abstract).

WAN Xing, DING Jing, ZHANG Xiao-li. 2006. Set Pair Analysis for Comprehensive Evaluation on Regional Groundwater Resources Carrying Capacity[J]. Urban Environment & Urban Ecology, 19(2): 8-10(in Chinese with English abstract).

WANG Jia-bing. 2012. Leakage recharge from pores saline groundwater to deep fresh groundwater on the condition of pumping in Huabei Plain—A case of Tianjing Plain[J]. Hydrogeology & Engingeering Geology, (6): 35-37(in Chinese with English abstract).

WANG Kang, ZHANG Ren-duo. 2011. Heterogeneous soil water flow and macropores described with combined tracers of dye and iodine[J]. Journal of Hydrology, 397: 105-117.

WANG Shi-qin, SHAO Jing-li, SONG Xian-feng, ZHANG Yong-bo, HUO Zhi-bin, ZHOU Xiao-yuan. 2008. Application of MODFLOW and geographic information system to groundwater flow simulation in North China Plain, China[J]. Environmental geology, 55: 1449-1462.

WANG Shun-jiu, YANG Zhi-feng, DING Jing. 2004. Projection Pursuit Method of Comprehensive Evaluation on Groundwater Resources Carrying Capacity in Guanzhong plain[J]. Resources Science, 26(6): 106-110(in Chinese with English abstract).

XIA Jun, ZHANG Lu, LIU Chang-ming, YU Jing-jie. 2007. Towards better water security in North China[J]. Water Resources Management, 21(1): 233-247.

YANG Shao-hua. 2009. Numerical simulation study on land subsidence caused by urban construction[D]. Beijing: China University of Geosciences(in Chinese with English abstract).

ZHA Yuan-yuan, SHI Liang-sheng, YE Ming, YANG Jin-zhong. 2013. A generalized Ross method for two- and three-dimensional variably saturated flow[J]. Advances in Water Resources, 54: 67-77.

ZHANG Xin, WANG Ji-ke, CAI Huan-jie, WANG Zheng-xing. 2001. Comprehensive Evaluation on Resources Carrying Capacity of Regional Groundwater[J]. Bulletin of Soil and Water Conservation, 21(3): 24-27(in Chinese with English abstract).

ZHANG Zhao-ji, FEI Yu-hong, CHEN Zong-yu, ZHAO Zong-zhuang, XIE Zhen-hua, WANG Ya-bin, MIAO Pu-xiang, YANG Li-zhi, SHAO Jing-li, JIN Meng-gui, XU Guang-ming, YANG Qi-qing. 2009a. The Sustainable Utilization Evaluation Survey of Groundwater in North China Plain[M]. Beijing: Geological Publishing House(in Chinese).

ZHANG Zhao-ji, FEI Yu-hong. 2009b. Atlas of Groundwater Sustainable Utilization in North China Plain[M]. Beijing: Sino Maps Press(in Chinese).

ZHANG Zong-hu, SHI De-hong, SHEN Zhao-li, ZONG Zou-xing, XUE Yu-qun. 1997. Evolution and Development of Groundwater Environment in North China Plain under Human Activities[J]. Acta Geoscientia Sinica, 18(4): 337-344(in Chinese with English abstract).

ZUO Wen-zhe. 2006. Survey and Research on Seawater Intrusion in the Yandaihe Plain of Qinhuangdao City[D]. Beijing: China University of Geosciences(in Chinese with English abstract).

Evolution Mechanism and Control of Groundwater in the North China Plain

SHI Jian-sheng1), LI Guo-min2), LIANG Xing3), CHEN Zong-yu1), SHAO Jing-li4), SONG Xian-fang5)
1) Institute of Hydrogeology and Environmental Geology, Chinese Academy of Geological Sciences, Shijiazhuang, Hebei 050061; 2) Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029; 3) China University of Geosciences(Wuhan), Wuhan, Hubei 430074; 4) China University of Geosciences(Beijing), Beijing 100083; 5) Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101

The long-term overexploitation of groundwater in the North China Plain has produced complicated groundwater depression funnels and caused a series of environmental problems. To realize the harmonization between the regional economic development and the groundwater resources and environment, researchers have implemented since 2010 the project “Evolution Mechanism and Management of the Groundwater in North China”, which is a National Basic Research Program (973) of China. On the basis of numerous previous studies, the evolution of groundwater flowing field and its response to human activities in the past 50 years were analyzed, and the change of water flowing parameters and the structure of vadose zone were investigated. The evolution of ground depression, the movement of saline water bodies, and the composition of deep groundwater were studied. A high resolution numerical groundwater flow model was built up based on renewed data between 2001 and 2010 and improved groundwater modeling technologies. Based on this model, the authors calculated the up-to-date amount of groundwater resources. The groundwater resources carrying capacities (GRCC) were assessed and the main factors affecting the GRCC were analyzed. The optimizing solutions for the management of the groundwater resources in the North China Plain were put forward and were modeling by the newly-built high resolutiongroundwater flow model.

North China Plain; groundwater depression funnel; groundwater resources management

P588.247; P534.42

A

10.3975/cagsb.2014.05.01

本文由國家重點基礎(chǔ)研究計劃(“973”)項目“華北平原地下水演變機制與調(diào)控”(編號: 2010CB428800)資助。獲中國地質(zhì)科學院 2013年度十大科技進展第四名。

2014-04-14; 改回日期: 2014-05-18。責任編輯: 閆立娟。

石建省, 男, 1962年生。研究員, 博士生導師。長期從事地下水資源調(diào)查評價、中國黃土研究、脆弱帶地質(zhì)生態(tài)學等研究。E-mail: tiger7886@263.net。