劉美英，閔雷雷，沈彥俊

(1.中國(guó)科學(xué)院遺傳與發(fā)育生物學(xué)研究所農(nóng)業(yè)資源研究中心/中國(guó)科學(xué)院農(nóng)業(yè)水資源重點(diǎn)實(shí)驗(yàn)室，河北 石家莊 050021；2.中國(guó)科學(xué)院大學(xué)，北京 100049)

河北平原地處華北平原北部，多年平均降水量為529.1 mm，人均水資源量和畝均水資源量?jī)H占到全國(guó)平均水平的11.18%和9.13%[1]，是海河流域典型的灌溉糧食產(chǎn)區(qū)[2-5]。由于地區(qū)地表水資源較少，灌溉主要靠抽取地下水，農(nóng)用地下水占總地下水開(kāi)采量的74.5%～76.6%，河北平原是我國(guó)當(dāng)前水資源供需矛盾最為嚴(yán)重的地區(qū)之一[2, 6,7, 9]。近年來(lái)，一些學(xué)者在該地區(qū)開(kāi)展了地下水位動(dòng)態(tài)變化分析的研究，張永強(qiáng)等以欒城縣為例，分析了華北平原地下水位動(dòng)態(tài)變化，欒城縣地下水位多年平均下降速率為0.7m/a[8]；許月卿等利用GIS技術(shù)，分析了1975-2000年河北平原(北京、天津以南)淺層、深層地下水位波動(dòng)狀況[9]；劉志國(guó)等在地理信息系統(tǒng)的支持下，運(yùn)用地統(tǒng)計(jì)學(xué)方法研究了1986年到2000年地下水位時(shí)空變化規(guī)律[10]；王仕琴等分析華北平原39處淺層地下水監(jiān)測(cè)結(jié)果，將地下水動(dòng)態(tài)分為六大類型[11]。此外，張光輝等人針對(duì)冀中山前平原的農(nóng)業(yè)區(qū)淺層地下水位不斷下降問(wèn)題，研究了地下水強(qiáng)降弱升機(jī)制[12]。受人口、地下水開(kāi)采量、降水和徑流等諸多因素的影響，河北平原地下水水位不斷變化，并由此衍生了一系列的生態(tài)環(huán)境問(wèn)題，例如地面沉降、塌陷、土壤次生鹽漬化等[13-17]。以往的研究受到數(shù)據(jù)資料的限制，通常是基于月尺度或年尺度水位的數(shù)據(jù)，因而對(duì)水位動(dòng)態(tài)變化的分析不夠深入。本文基于日尺度的自動(dòng)觀測(cè)數(shù)據(jù)分析河北平原區(qū)域內(nèi)的淺層地下水位動(dòng)態(tài)變化，研究結(jié)果有助于深入了解區(qū)域地下水循環(huán)規(guī)律，可為區(qū)域地下水資源的可持續(xù)利用和區(qū)域經(jīng)濟(jì)社會(huì)協(xié)同發(fā)展提供科學(xué)依據(jù)。

1 研究區(qū)概況和監(jiān)測(cè)點(diǎn)分布

1.1 研究區(qū)概況

本文所指的河北平原是石家莊、保定、衡水、滄州以及廊坊的部分區(qū)域(圖1)，從南北方向上，位于河北省中部。地勢(shì)較低平，主要由黃河和海河等不斷沖擊而成，海拔由西部山前平原100 m左右降低到東部濱海平原的3 m左右，屬于溫帶季風(fēng)性氣候，四季變化明顯，冬季干燥寒冷，夏季高溫多雨，春季干旱少雨，蒸發(fā)強(qiáng)烈，年均氣溫和年降水量由南向北隨緯度增加而遞減。春季小麥需水量大，降水量較小，旱情較重，夏季降水量大，常有洪澇。河北平原是中國(guó)嚴(yán)重缺水的地區(qū)之一，人均水資源量?jī)H為335 m3/a，地下水是該地區(qū)的主要用水來(lái)源[18]。

圖1 研究區(qū)地形及地下水監(jiān)測(cè)點(diǎn)分布圖Fig. 1 Location of groundwater monitoring sites in the research area

河北平原是個(gè)大型的沉降盆地，在其地質(zhì)演化過(guò)程中，沉積了較厚層的松散物質(zhì)，形成了河北平原主要的地下水流動(dòng)體系。河北平原自太行山區(qū)向中部平原呈現(xiàn)出水文地質(zhì)條件的水平方向分帶性。自西向東在含水層的巖性、分布、厚度和地下水的化學(xué)成分與運(yùn)動(dòng)特征等方面均有顯著變化，南北方向變化較小。含水層的巖性成分具有以粗顆粒為主變?yōu)橐约?xì)顆粒為主的區(qū)域性變化規(guī)律[19]。區(qū)域內(nèi)，淺層地下水補(bǔ)給的來(lái)源主要是降水入滲補(bǔ)給和灌溉回水補(bǔ)給，此外還有山前側(cè)向徑流補(bǔ)給、河道及渠系滲漏補(bǔ)給等[20]。

1.2 地下水位監(jiān)測(cè)

我國(guó)地下水水位的監(jiān)測(cè)從20世紀(jì)50年代就已經(jīng)開(kāi)始[21]，但是大多數(shù)監(jiān)測(cè)系統(tǒng)由于時(shí)間較長(zhǎng)，數(shù)據(jù)質(zhì)量不高。為深入研究河北平原地下水變化，本文共在河北平原設(shè)置11個(gè)監(jiān)測(cè)點(diǎn)(圖1和表1)，這11個(gè)監(jiān)測(cè)點(diǎn)是從2008年開(kāi)始，依次在河北平原(石家莊、保定、衡水、廊坊)選取了11眼灌溉或生產(chǎn)生活觀測(cè)孔(有些監(jiān)測(cè)點(diǎn)的數(shù)據(jù)存在不同程度缺失)。機(jī)井位置的選擇考慮了數(shù)據(jù)的代表性，監(jiān)測(cè)點(diǎn)遠(yuǎn)離工廠等地下水嚴(yán)重開(kāi)采區(qū)，此外還考慮到儀器管理的安全性和數(shù)據(jù)收集的便利性等，部分監(jiān)測(cè)點(diǎn)在監(jiān)測(cè)過(guò)程中存在數(shù)據(jù)缺失現(xiàn)象。

本研究采用的水位自動(dòng)記錄儀是ONSET公司生產(chǎn)的HOBO Water Level Logger水位計(jì)。地下水位監(jiān)測(cè)頻率為1 h。水位計(jì)可以對(duì)地下水埋深的變化進(jìn)行連續(xù)自動(dòng)監(jiān)測(cè)，進(jìn)而反映地下水位的動(dòng)態(tài)變化。每次采集數(shù)據(jù)時(shí)，均采用加拿大Solinst公司102型地下水位計(jì)對(duì)地下水埋深進(jìn)行校準(zhǔn)驗(yàn)證。

2 數(shù)據(jù)處理

對(duì)于上述11個(gè)HOBO測(cè)得的1 h一次的觀測(cè)數(shù)據(jù)進(jìn)行處理。由于觀測(cè)井多用于生活和農(nóng)業(yè)灌溉，地下水抽取多發(fā)生在白天，觀測(cè)井的水位受人為開(kāi)采地下水的影響存在明顯的日變化，本文選擇凌晨4時(shí)相對(duì)靜態(tài)穩(wěn)定時(shí)間段的水位作為當(dāng)天的水位觀測(cè)值進(jìn)行分析。

3 地下水位動(dòng)態(tài)變化分析

3.1 不同地區(qū)地下水埋深動(dòng)態(tài)變化

河北平原地下水埋深動(dòng)態(tài)變化過(guò)程受到多重因素的影響，例如地形地貌、降水量、土地利用類型、湖泊水庫(kù)等自然因素，人類治水工程、開(kāi)采地下水等人為因素[22]。且這些因素本身也存在時(shí)空變異，變化速率也有差異，導(dǎo)致地下水位呈現(xiàn)出不同的動(dòng)態(tài)變化特征。

3.1.1 山前平原區(qū)地下水埋深動(dòng)態(tài)變化

山前平原區(qū)的監(jiān)測(cè)點(diǎn)主要分布在石家莊和保定。區(qū)域監(jiān)測(cè)點(diǎn)呈現(xiàn)出的整體趨勢(shì)：年內(nèi)存在較大水位波動(dòng)，除7號(hào)監(jiān)測(cè)點(diǎn)(圖2和圖3)外年平均地下水位不斷下降。但由于不同監(jiān)測(cè)點(diǎn)的地貌地形、周圍土地利用現(xiàn)狀以及地下水位埋深等存在空間差異，該區(qū)域各監(jiān)測(cè)點(diǎn)地下水位變化過(guò)程又呈現(xiàn)出不同動(dòng)態(tài)特征。

圖2 山前平原區(qū)1、2號(hào)監(jiān)測(cè)點(diǎn)地下水埋深動(dòng)態(tài)變化曲線Fig.2 The curves of dynamic variation of groundwater depth of No.1, 2 in the region of the piedmont plain.

圖3 山前平原3、4、5號(hào)監(jiān)測(cè)點(diǎn)地下水埋深動(dòng)態(tài)變化曲線Fig.3 The curves of dynamic variation of groundwater depth of No.3, 4, 5 in the region of the piedmont plain

其中石家莊的1、2號(hào)監(jiān)測(cè)點(diǎn)(圖2)，位于河北山前平原最西部。1號(hào)監(jiān)測(cè)點(diǎn)位于山前平原最西部海拔較高的地區(qū)，含水層的滲透性較好，以及地下水位埋深較淺，地下水位受到山區(qū)側(cè)向補(bǔ)給和降雨補(bǔ)給的作用明顯，地下水位下降幅度較小，在降水量較小的2014年(降水量為294.8 mm)地下水位明顯下降。該點(diǎn)的年內(nèi)最低水位出現(xiàn)7月底、8月初，在雨季來(lái)臨之后1個(gè)月，9月份地下水位達(dá)到年內(nèi)最高值，年內(nèi)最高水位與最低水位差在2～3 m之間，9月份達(dá)到最高水位之后又不斷下降，直到翌年7月底8月初。2號(hào)監(jiān)測(cè)點(diǎn)位于石家莊市區(qū)，地下水位埋深較深，但由于靠近山區(qū)，同樣受到山前側(cè)向補(bǔ)給作用明顯，且在雨季來(lái)臨后，接受到山區(qū)側(cè)向補(bǔ)給的時(shí)間較快，年內(nèi)最低水位出現(xiàn)在4-5月份，最高水位出現(xiàn)的時(shí)間與年內(nèi)最大降水量出現(xiàn)的時(shí)間相對(duì)應(yīng)。由于該點(diǎn)周圍生活用水主要靠抽取地下水，開(kāi)采程度較高，地下水位下降的平均速率為1.2 m/a。

3、4、5號(hào)監(jiān)測(cè)點(diǎn)(圖3)分布在山前開(kāi)采程度高的地區(qū)，屬于農(nóng)業(yè)用水集中區(qū)，水位埋深較大，特別是5號(hào)監(jiān)測(cè)點(diǎn)水位埋深大于 45 m。該區(qū)域地下水變化的特點(diǎn)是水位在年際尺度上不斷下降。據(jù)調(diào)查，3號(hào)監(jiān)測(cè)點(diǎn)在2013年之前農(nóng)業(yè)灌溉采用的是污水灌溉，減少了農(nóng)業(yè)灌溉地下水的開(kāi)采，因此在2013年之前年際地下水位無(wú)明顯波動(dòng)，2014年起農(nóng)業(yè)灌溉不再使用污水灌溉，轉(zhuǎn)而依靠抽取地下水，在降水量與前一年基本相同的情況下(2013年降水量596 mm，2014年降水量609 mm)，相比2013年，2014年該監(jiān)測(cè)點(diǎn)水位出現(xiàn)下降趨勢(shì)，下降3 m左右。這3個(gè)監(jiān)測(cè)點(diǎn)所在區(qū)域年內(nèi)地下水位存在較大波動(dòng)，年內(nèi)最高水位和最低水位差值為3～5 m，地下水位變化對(duì)降水量的響應(yīng)比較緩慢，在雨季來(lái)臨之后，水位沒(méi)有立即回升至年內(nèi)最大值，年內(nèi)最高水位出現(xiàn)在2-3月份，最低水位出現(xiàn)在7-8 月份。這主要是由農(nóng)業(yè)灌溉用水引起的，3-7月份是小麥灌溉用水的集中時(shí)間段，小麥生長(zhǎng)季的需水量遠(yuǎn)大于降水量，過(guò)量抽取地下水使得2月份到7月份地下水位持續(xù)下降，7月份之后由于降水量的增加以及灌溉用水的減少，地下水位不斷回升，直到翌年2-3月份。

圖4 山前平原7號(hào)監(jiān)測(cè)點(diǎn)地下水埋深動(dòng)態(tài)變化曲線Fig.4 The curves of dynamic variation of groundwater depth of No.7 in the region of the piedmont plain

7號(hào)監(jiān)測(cè)點(diǎn)位于保定西大洋水庫(kù)附近，且距離山區(qū)較近，海拔相對(duì)較高，地下水位呈現(xiàn)出明顯區(qū)別于其他觀測(cè)點(diǎn)的特征。觀測(cè)年份內(nèi)，7號(hào)監(jiān)測(cè)點(diǎn)(圖4)年際地下水位不斷上升，2011-2013年的年平均地下水埋深分別為8.37、6.28和5.85 m，2011-2012年地下水位大幅度回升。此外由于地下水埋深淺，受到降雨入滲補(bǔ)給和蒸發(fā)排泄作用明顯，同時(shí)受到山區(qū)和水庫(kù)側(cè)向補(bǔ)給作用的影響，年內(nèi)地下水埋深在雨季來(lái)臨之前達(dá)到年內(nèi)最低值，在雨季來(lái)臨之后不斷回升，1～2個(gè)月后達(dá)到年內(nèi)最大值，之后不斷下降。

分析7號(hào)監(jiān)測(cè)點(diǎn)地下水位回升的原因。西大洋水庫(kù)作為保定飲用水水源地，年均供水量為5 000 萬(wàn)m3，但由于水庫(kù)蓄水量嚴(yán)重不足，自2007年底保定市王快水庫(kù)至西大洋水庫(kù)聯(lián)通工程建成通水，正常年份王快水庫(kù)每年可向西大洋水庫(kù)輸水兩億方[23,24]，西大洋水庫(kù)供水不足問(wèn)題得到緩解，這是7號(hào)監(jiān)測(cè)點(diǎn)地下水位在監(jiān)測(cè)年份內(nèi)不斷回升的主要原因。與此同時(shí)，農(nóng)業(yè)用水也是引起河北平原地下水位下降的主要因素，河北平原農(nóng)作物的主要類型是冬小麥和夏玉米，小麥?zhǔn)且鸬叵滤幌陆捣茸畲蟮淖魑锓N類[25]，當(dāng)降水量并不能滿足作物生長(zhǎng)需水時(shí)，必須依靠抽取地下水進(jìn)行灌溉。據(jù)河北農(nóng)村統(tǒng)計(jì)年鑒數(shù)據(jù)，7號(hào)監(jiān)測(cè)點(diǎn)所在縣(曲陽(yáng)縣)的小麥播種面積在不斷減少[圖5(a)]。2000-2006年的年平均降水量為439.17 mm，2007-2013年的年平均降水量為566.05mm，降水量增加了28.9%[圖5(b)]。小麥種植面積減少，農(nóng)業(yè)灌溉用水相對(duì)減少，地下水的開(kāi)采量減少；年平均降水量增加，降水對(duì)地下水的補(bǔ)給量增加，也可能是7號(hào)監(jiān)測(cè)點(diǎn)地下水埋深不斷減小的另一原因[26]。

圖5 曲陽(yáng)縣2000-2013年小麥播種面積和降水量變化圖Fig.5 The variation diagram of wheat planted area and precipitation in Qu yang county in 2000-2013

3.1.2 中部平原區(qū)地下水埋深動(dòng)態(tài)變化

中部平原區(qū)的監(jiān)測(cè)點(diǎn)主要分布在衡水北部、保定西部、滄州和廊坊的東部。地下水的整體特征為地下水埋深較淺，從西到東地下水埋深越來(lái)越小，平均地下水埋深從西邊的9號(hào)監(jiān)測(cè)點(diǎn)24 m到東部11號(hào)監(jiān)測(cè)點(diǎn)8 m。

年際地下水位無(wú)明顯升高降低趨勢(shì)，觀測(cè)年份內(nèi)各水文年水位波動(dòng)小于3 m(表2)，處于動(dòng)態(tài)穩(wěn)定狀態(tài)。年內(nèi)地下水位存在較大波動(dòng)，隨著地下水開(kāi)采量和降水量的變化而不斷變化。由于中部平原地區(qū)地下水埋深較淺，地下水開(kāi)發(fā)利用程度較低，地下水位動(dòng)態(tài)變化主要受到降雨入滲、蒸發(fā)排泄和人為開(kāi)采的影響。最低水位出現(xiàn)在5-6月份，此時(shí)間段降水量較少。進(jìn)入雨季之后，受到降雨入滲的影響，地下水位不斷回升。例如6號(hào)和9號(hào)監(jiān)測(cè)點(diǎn)，地下水位平均埋深分別為15.02、23.14 m，每年最低水位出現(xiàn)在6月底，之后隨著雨季的到來(lái)，降水量不斷增加，地下水位開(kāi)始不斷回升，在10-11月份由于降水量減少，地下水位出現(xiàn)輕微下降，之后不斷回升達(dá)到年內(nèi)最高值。而8號(hào)、10號(hào)和11號(hào)監(jiān)測(cè)點(diǎn)，地下水位平均埋深分別為10.17、6.66和7.69 m，相對(duì)7、10號(hào)監(jiān)測(cè)點(diǎn)，地下水位埋深更淺，地下水位對(duì)降水量的響應(yīng)更快，每年的最低水位出現(xiàn)在5月底、6月初，隨著雨季的到來(lái)，降水量增加，地下水位開(kāi)始快速回升，9月份以后基本達(dá)到年內(nèi)最高值，之后地下水位無(wú)明顯下降，保持著水位高值，直到翌年的3月底、4月初，由于灌溉開(kāi)采量增大和降水量較小，地下水位開(kāi)始持續(xù)下降。在過(guò)去的1980年到2003年，該區(qū)域地下水位也在不斷波動(dòng)，分別經(jīng)歷了穩(wěn)定下降、上升、相對(duì)穩(wěn)定、逐漸下降等階段[20](圖6)。

表2 研究區(qū)各監(jiān)測(cè)點(diǎn)地下水埋深 m

注：監(jiān)測(cè)數(shù)據(jù)在不同年份內(nèi)有數(shù)據(jù)缺失現(xiàn)象。

圖6 中部平原區(qū)6、8、9、10、11號(hào)監(jiān)測(cè)點(diǎn)地下水埋深動(dòng)態(tài)變化曲線Fig.6 The curves of dynamic variation of groundwater depth of No.6, 8, 9, 10, 11 in the region of the central plain

3.2 典型觀測(cè)點(diǎn)年內(nèi)地下水埋深波動(dòng)對(duì)比

不同地區(qū)年內(nèi)地下水位的波動(dòng)變化方式有很大的區(qū)別，選取不同地區(qū)的3個(gè)典型監(jiān)測(cè)點(diǎn)同一時(shí)間段(2010/1/1-2011/7/1)，分析年內(nèi)地下水位變化機(jī)制(圖7)。位于山前平原最西部海拔較高地區(qū)的1號(hào)監(jiān)測(cè)點(diǎn)，水位埋深較淺，年際地下水位下降程度較小(圖2)；位于山前平原的5號(hào)監(jiān)測(cè)點(diǎn)，水位埋深較大，在觀測(cè)年份內(nèi)地下水位不斷下降(圖3)；位于中部平原的6號(hào)監(jiān)測(cè)點(diǎn)，地下水埋深相對(duì)山前平原較淺，在觀測(cè)年份內(nèi)，年際地下水位基本不變(圖6)。地下水位的年際動(dòng)態(tài)變化是由地下水位的年內(nèi)水位升高降低機(jī)制反映出來(lái)的[12]。

在山前平原西部海拔較高的地區(qū)(以1號(hào)監(jiān)測(cè)點(diǎn)為例)，年內(nèi)地下水位變化比較平穩(wěn)，年初地下水位先降低后回升，回升3個(gè)月很快達(dá)到年內(nèi)最高值，之后又不斷下降，水位下降期日變化速率小于水位回升期日變化速率，但是年內(nèi)下降期的時(shí)間遠(yuǎn)大于回升期的時(shí)間，水位在年際上緩慢下降。

在山前低平原區(qū)(以5號(hào)監(jiān)測(cè)點(diǎn)為例)，年內(nèi)地下水位的波動(dòng)相對(duì)于中部平原區(qū)比較平穩(wěn)，年初水位先回升，之后不斷下降，8月左右水位開(kāi)始再次回升。年內(nèi)水位變化趨勢(shì)與山前丘陵區(qū)存在較大的差距，水位下降期日變化速率大于水位回升期日變化速率。在小麥農(nóng)灌時(shí)期，地下水開(kāi)采量較大，補(bǔ)給量較小，日下降速率大于1 cm/d，在雨季來(lái)臨后，降水量增加，并且農(nóng)業(yè)灌溉開(kāi)采用水減少，水位開(kāi)始回升，但是回升速率遠(yuǎn)小于1 cm/d。雖然年內(nèi)回升期的時(shí)間略大于下降期的時(shí)間，但是由于回升期日變化速率遠(yuǎn)小于下降期的日變化速率，下降速度為上升速度的2.1～3.6倍，年際地下水位呈不斷下降的趨勢(shì)。

在中部平原區(qū)(以6號(hào)監(jiān)測(cè)點(diǎn)為例)，年內(nèi)水位波動(dòng)趨勢(shì)與山前平原大致相同，但存在較多的小波動(dòng)，由于水位埋深較淺，水位對(duì)降水量的反應(yīng)較快，春灌期后水位回升的時(shí)間比山前平原早一些，水位下降期日變化速率大于水位回升期日變化速率，下降速度為上升速度的1.7～6.7倍，但是上升時(shí)間約為下降時(shí)間的4倍，年際地下水位基本穩(wěn)定。

4 結(jié) 語(yǔ)

(1)河北平原地下水動(dòng)態(tài)變化呈現(xiàn)出空間分布的整體規(guī)律，平均地下水埋深由山前平原區(qū)到中部平原區(qū)不斷減小。山前平原區(qū)地下水開(kāi)發(fā)利用程度較高，地下水埋深較深，中部平原區(qū)地下水開(kāi)發(fā)利用程度較低，地下水埋深較淺。地下水年際變化分布規(guī)律為：中部平原區(qū)地下水埋深基本不變，觀測(cè)年份內(nèi)各水文年水位波動(dòng)小于3 m，山前平原區(qū)地下水埋深不斷下降，不同監(jiān)測(cè)點(diǎn)下降幅度差異較大。局部地區(qū)，由于靠近水庫(kù)，水庫(kù)蓄水量增加，且小麥種植面積減少和降水量增加，地下水位呈回升趨勢(shì)。地下水年內(nèi)變化規(guī)律為：山前平原最西部海拔較高地區(qū)，受到地貌地形、地下水滲透系數(shù)的影響，年內(nèi)最低水位出現(xiàn)在7-8月份，最高水位出現(xiàn)在9月份；山前平原農(nóng)業(yè)用水集中區(qū)，地下水開(kāi)發(fā)利用程度較高，受到降雨入滲、農(nóng)業(yè)灌溉用水和側(cè)向補(bǔ)給作用的影響，年內(nèi)最高水位出現(xiàn)在2-3月份，最低水位出現(xiàn)在7月份；中部平原，地下水埋深較淺，地下水變化受到蒸發(fā)排泄和降雨入滲的雙重影響，年內(nèi)最低水位出現(xiàn)在4-6月份，最高水位出現(xiàn)在9月份左右。

2)不同地區(qū)年內(nèi)水位波動(dòng)存在較大差異。從山前平原區(qū)到中部平原區(qū)，海拔相對(duì)較高的地區(qū)與低平原區(qū)的年內(nèi)水位下降期和水位上升期分布差異較大，低平原區(qū)年內(nèi)水位下降時(shí)間小于水位回升時(shí)間，而海拔相對(duì)較高的地區(qū)則相反，且兩個(gè)區(qū)域水位下降期和水位上升期日變化速率也不相同。山前低平原區(qū)和中部平原區(qū)，年內(nèi)水位下降期和水位上升期分布大致相同，但是中部平原區(qū)年內(nèi)回升期開(kāi)始的比山前低平原區(qū)早，可能是由于水位埋深較淺，對(duì)降雨的響應(yīng)更及時(shí)，且中部平原區(qū)的年內(nèi)回升期的時(shí)間較長(zhǎng)。

本研究分析了河北平原地下水位動(dòng)態(tài)變化規(guī)律，同時(shí)對(duì)比了區(qū)域內(nèi)的典型觀測(cè)點(diǎn)的同時(shí)間年內(nèi)水位波動(dòng)狀況，是該地區(qū)水位動(dòng)態(tài)變化研究的基礎(chǔ)，下一步將結(jié)合土地利用類型和所處區(qū)位在地下水流場(chǎng)上的不同位置深入分析水位變動(dòng)的規(guī)律和原因。

□

[1] LIU chang-ming, YU jing-jie, Eloise Kendy. Groundwater exploitation and its impact on the environment in the North China Plain [J]. Water International, 2001,26(2):265-272.

[2] 袁再健, 許元?jiǎng)t, 謝櫨樂(lè). 河北平原農(nóng)田耗水與地下水動(dòng)態(tài)及糧食生產(chǎn)相互關(guān)系分析[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 2014,22(8):904-910.

[3] 毛學(xué)森, 劉昌明. 太行山山前平原地下水變化趨勢(shì)與農(nóng)業(yè)持續(xù)發(fā)展[J]. 水土保持研究,2001,(1):147-149.

[4] 張光輝, 費(fèi)宇紅, 王惠軍等. 河北省平原區(qū)農(nóng)田糧食增產(chǎn)與灌溉節(jié)水對(duì)地下水開(kāi)采量的影響[J]. 地質(zhì)通報(bào), 2009,(5):645-650.

[5] 袁再健, 謝櫨樂(lè), 張秉文等. 河北平原農(nóng)田凈灌溉耗水研究[J]. 南水北調(diào)與水利科技,2015,(4):780-784.

[6] 許月卿. 京津以南河北平原地下水位下降驅(qū)動(dòng)因子的定量評(píng)估[J]. 地理科學(xué)進(jìn)展,2003,(5):490-498.

[7] 張衛(wèi)東, 何慶成, 房 浩. 華北平原開(kāi)采地下水產(chǎn)生的正效益估算[J]. 中國(guó)地質(zhì)災(zāi)害與防治學(xué)報(bào),2006,(3):110-113.

[8] 張永強(qiáng), 劉昌明, 沈彥俊. 太行山山前平原淺層地下水位動(dòng)態(tài)分析----以河北省欒城縣為例[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào),2001,(2):42-44.

[9] 許月卿, 蔡運(yùn)龍. 基于 GIS 的河北平原地下水位時(shí)空變化動(dòng)態(tài)分析[J]. 北京大學(xué)學(xué)報(bào): 自然科學(xué)版, 2005,41(2):265-272.

[10] 劉志國(guó), 王恩德, 付建飛, 等. 河北平原地下水水位的時(shí)空變異[J]. 東北大學(xué)學(xué)報(bào) (自然科學(xué)版), 2007,28(5):717-720.

[11] 王仕琴, 宋獻(xiàn)方, 王勤學(xué), 等. 華北平原淺層地下水水位動(dòng)態(tài)變化[J]. 地理學(xué)報(bào), 2008,63(5):462-472.

[12] 張光輝, 田言亮, 王電龍, 等. 冀中山前農(nóng)業(yè)區(qū)地下水位強(qiáng)降弱升特征與機(jī)制[J]. 水科學(xué)進(jìn)展, 2015,26(2):227-232.

[13] 張宗祜, 沈照理. 人類活動(dòng)影響下華北平原地下水環(huán)境的演化與發(fā)展[J]. 地球?qū)W報(bào): 中國(guó)地質(zhì)科學(xué)院院報(bào), 1997,18(4):337-344.

[14] 范鵬飛. 華北平原地下水演化及預(yù)測(cè)[J]. 地球?qū)W報(bào): 中國(guó)地質(zhì)科學(xué)院院報(bào), 1998,19(4):346-352.

[15] 夏 軍, 劉孟雨, 賈紹鳳, 等. 華北地區(qū)水資源及水安全問(wèn)題的思考與研究[J]. 自然資源學(xué)報(bào), 2004,19(5):550-560.

[16] 段永侯,肖國(guó)強(qiáng). 河北平原地下水資源與可持續(xù)利用[J]. 水文地質(zhì)工程地質(zhì),2003,(1):2-8.

[17] 高太忠,羅人明,黃群賢,等. 河北平原地下水開(kāi)發(fā)中的問(wèn)題及可持續(xù)利用[J]. 地質(zhì)災(zāi)害與環(huán)境保護(hù),2002,(3):6-9.

[18] 張光輝, 費(fèi)宇紅, 劉春華, 等. 華北平原灌溉用水強(qiáng)度與地下水承載力適應(yīng)性狀況[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2013,29(1):1-10.

[19] 張宗祜. 河北平原水文地質(zhì)條件及其在水利規(guī)劃中的意義[J]. 中國(guó)地質(zhì)，1961,(10):001.

[20] 王春澤, 張石春, 張建平. 河北平原區(qū)近 30 年地下水動(dòng)態(tài)及可持續(xù)利用對(duì)策[J]. 河北水利, 2005,(2):14-16.

[21] 周仰效, 李文鵬. 區(qū)域地下水位監(jiān)測(cè)網(wǎng)優(yōu)化設(shè)計(jì)方法[J]. 水文地質(zhì)工程地質(zhì), 2007,34(1):1-9.

[22] 費(fèi)宇紅, 張兆吉, 張鳳娥, 等. 華北平原地下水位動(dòng)態(tài)變化影響因素分析[J]. 河海大學(xué)學(xué)報(bào): 自然科學(xué)版, 2005,33(5):538-541.

[23] 袁軍鵬,袁吉棟,史長(zhǎng)瑩. 王快-西大洋水庫(kù)跨流域調(diào)水工程研究[J]. 南水北調(diào)與水利科技,2007,(5):151-153.

[24] 史長(zhǎng)瑩,簡(jiǎn)新平,李麗英. 王快-西大洋水庫(kù)聯(lián)合調(diào)度研究[J]. 黑龍江大學(xué)工程學(xué)報(bào),2012,(1):7-10.

[25] 許月卿. 土地利用對(duì)地下水位下降的影響----以河北平原為例[J]. 地理研究,2005,(2):222-228.

[26] 王 學(xué),李秀彬,辛良杰. 河北平原冬小麥播種面積收縮及由此節(jié)省的水資源量估算[J]. 地理學(xué)報(bào),2013,(5):694-707.