韓宇平　劉存強　趙雨婷　劉中培　王春穎

摘要：基于MODFLOW模型，以華北平原人民勝利渠灌區(qū)為例，結合研究區(qū)水文地質條件，建立地下水運動數值模擬模型，通過參數校準和模型驗證，表明建立的地下水運動模型能夠合理地反映研究區(qū)2012-2013年的地下水運動狀況，模擬結果表明該灌區(qū)地下水處于負均衡狀態(tài)?；谀M結果，進一步預測了氣候情景（選取RCP4.5情景NorESM1-M模式）下灌區(qū)2030年地下水水位情況。結果表明與1997-2013年相比，2030年灌區(qū)地下水位持續(xù)下降，漏斗面積逐漸擴大。以此為基礎，開展地下水開采量情景分析，將開采量分別增加和減少20%，預測2030年在不同開采量情景下地下水水位變化情況。最后根據預測結果初步提出地下水開采量減少20%的調控方案，以保證地下水水位有所上升，漏斗面積減少。

關鍵詞：MODFLOW;灌區(qū);地下水;氣候變化;情景分析

中圖分類號：P641文獻標志碼：A

Abstract：The MODFLOW model was used to simulate groundwater flow in a typical irrigation area，Peoples' Victory Canal Irrigation District，in North China Plain，based on the hydrogeological conditions of the study area.The model parameters were well calibrated and model verification results showed that the simulation results could reasonably reflect the groundwater dynamics in the study area from 2012 to 2013.Under agricultural groundwater exploitation and climate conditions during 2012-2013 period，groundwater balance was in a negative equilibrium state.The spatial distribution of groundwater flow in 2030 under the future climate change scenarios （RCP4.5 scenario NorESM1-M model） was further predicted.The results showed that the groundwater level would continue to decline from 2013 to 2030，and the low groundwater level funnel area was gradually expanding.The spatial distribution of groundwater flow in 2030 was further expllored under different groundwater exploitation scenarios，i.e.，the exploitation quantity increased and decreased by 20%，respectively.According to the forecast results，a 20% reduction in groundwater exploitation was projected to ensure that the groundwater level could be raised and the funnel area could be reduced.However，an increase in exploitation would lead to an obvious increase in the groundwater funnel area.

Key words：MODFLOW;irrigation area;groundwater;climate change;scenario analysis

農業(yè)用水是我國水資源利用的主體，多年平均用水量約占全國總用水量的70%以上。其中地下水是我國農業(yè)灌溉的重要水源，在城鄉(xiāng)生活用水、社會經濟發(fā)展、生態(tài)平衡等方面也起著非常重要的作用[1]。然而我國地下水采補平衡方面面臨嚴重問題，北方地區(qū)特別是黃淮海流域，地下水用水比例不斷上升[2]，地下水開采量的不斷增加引起了一系列生態(tài)環(huán)境問題[3]。研究表明地下水超采會導致水源地地下水位下降過大，形成地下水降落漏斗[4-6]，在沿海地區(qū)甚至還會造成海水入侵[7]。因此研究農業(yè)用水對地下水水位的影響對于地下水資源的可持續(xù)利用和農業(yè)的可持續(xù)發(fā)展非常必要[8]。

地下水作為水循環(huán)的主要組成部分其動態(tài)變化與人類活動和氣候變化密切相關。隨著人類對地下水開采量的增加，地下水的補排平衡被破壞，導致地下水水位逐年下降。與此同時，氣候變化導致的降水、蒸散發(fā)、潛水蒸發(fā)變化也影響著地下水補排關系[9]。為研究其對地下水動態(tài)的影響，國內外學者應用了各種方法和模型。王蕊等[10]建立了地表水地下水耦合模型，研究了南水北調工程實施后對白洋淀水文過程的影響。陳皓銳等[11]應用MODFLOW模擬了未來不同情景下未來的潛水位對氣候和人類活動的響應情況。Scibek等[12]耦合了氣候模式和地下水模型，研究Grand Forks地區(qū)不同氣候情景下地下水的補排情況。Kirshen[13]結合MODFLOW模型研究了全球變暖對馬薩諸塞州東部地下水的潛在影響，研究表明干旱情景下氣候變化對地下水的影響更為嚴重。Oude Essink等[14]建立了地下水運動模型，研究氣候變暖情況下荷蘭沿海地區(qū)地下水動態(tài)變化情況。廖梓龍[15]的研究表明，降雨補給變化對地下水水位變化影響非常大，同時人類開采活動也對地下水位變化造成很大影響。劉路廣等[16]基于MODFLOW模型模擬了柳園口灌區(qū)不同灌溉方式、地下水開采量和種植結構情況下地下水位的響應，結果表明井渠聯(lián)合可以很好地控制地下水位。王電龍等[17]應用了水量平衡原理、統(tǒng)計降尺度和概率統(tǒng)計等方法研究了不同氣候變化情景下華北平原井灌糧區(qū)未來50年糧食生產地下水保障能力，從糧食生產安全用水角度建議了合理的發(fā)展模式。林艷竹[18]應用了GMS及改進后的MODFLOW模型預測了現狀氣候和未來不同模式的長期氣候情景下地下水流場和動態(tài)變化趨勢，分析氣候變化對華北平原地下水流場和均衡的影響。王利書等[19]研究發(fā)現石羊河流域近20年來氣象因素中降水量、蒸發(fā)量的變化對地下水位影響較小，徑流量的變化對地下水位影響較大。大量觀測證據表明，由于大氣中溫室氣體濃度的增加，目前全球正經歷明顯的氣候變化[20-22]。以往研究表明氣候變化與地下水開采可能是造成地下水位變化的主要因素，因此，研究長期氣候變化和農業(yè)用水條件下地下水的演變規(guī)律研究及發(fā)展趨勢具有重要的意義。

2.2 未來氣候情景構建

氣候模式是目前預測未來氣候變化和影響的最主要的工具。氣候模式的輸出結果，通常是在一定溫室氣體和氣溶膠的排放情景（濃度情景或響應輻射強迫情景）下模擬得到[30]。本文選取聯(lián)合國政府間氣候變化專門委員會（IPCC）第五次評估報告中確定使用的“典型濃度路徑”（representative concentration pathways，RCPs）為未來的濃度情景。選取RCP4.5作為氣候模式的模擬情景，因為該情景是當前應用較為廣泛的一種排放情景，該情景下2100年大氣溫室氣體濃度穩(wěn)定在6.5×10-4CO2當量，輻射強迫穩(wěn)定在4.5 W/m2，大氣溫度預計平均升高3.6 ℃。氣候模式模擬數據來源于ISI-MIP提供，中國農科院環(huán)發(fā)所收集整理的中國大陸地區(qū)氣候變化情景數據集，本文選取NorESM1-M模式作為建立模型的氣候模式。該數據集經時間范圍為1951-2050年，空間水平分辨率為0.5°×0.5°。覆蓋研究區(qū)有兩個網格，以這兩個網格的模擬值的算術平均結果作為區(qū)域氣候變化情景。2014-2030年研究區(qū)蒸發(fā)量根據NorESM1-M模式輸出結果計算得到。Noresm1-M模式下2014-2030年灌區(qū)平均氣溫升高1.9 ℃，年平均降水量585 mm，年平均潛在蒸散量為1 335 mm。為闡明氣候變化對地下水的影響，基于2012-2013年地下水開采情況，2014-2030年氣候變化條件下模型模擬結果與1997-2013年歷史氣象條件（年平均降水量608 mm，年平均潛在蒸散量1 078 mm）下灌區(qū)地下水變化情況進行對比。另外，本文開展了氣候變化條件下開采方案情景分析，目的是為了分析未來氣候下人工開采量變化對地下水的影響，為提出合理的地下水開采方案提供依據?；谖磥砣蜃兣榫跋?，大氣溫度升高，作物需水量的增加，同時工業(yè)、生活、生態(tài)的用水量也增加，設置地下水開采量增加20%的情景?；谖磥硐拗频叵滤_采同時采用節(jié)水灌溉設施，設置地下水開采量降低20%的情景[31]。

3 結果與討論

3.1 模型校準與驗證

以2012年1月實測地下水位為MODFLOW模型的初始水位，模型的校準期選取2012年1月1日至2012年12月31日，校準期間通過反復調整參數使模型模擬值與實測值盡量接近以確定最優(yōu)參數（表2）。模型的驗證期選取2013年1月1日至2013年12月31日。模型模擬期間亢村、夏莊、后河、北翟坡、莊巖、張村6個地下水位觀測井的實測值和模擬值對比見圖4，模型擬合統(tǒng)計參數見表3。結果表明，觀測點在各個時段的模擬水位與觀測水位擬合程度較好（圖4）。模型校準期和驗證期歸一化均方根誤差（NRMS）分別為3.9%和3.1%，均方差（RMS）為0.77 m和0.56 m（表3）。以上結果表明MODFLOW模型能夠較好的模擬各觀測井地下水水位及反映研究區(qū)地下水系統(tǒng)的變化規(guī)律。MODFLOW模型模擬的數據，研究期間水量補給項、排泄項總水量見表4。由表4可看出，降水入滲是主要的補給項，地下水開采灌溉是主要的排泄項。研究期間灌區(qū)地下水量呈負均衡態(tài)勢，即地下水位呈下降趨勢。地下水過度開采是造成下水量呈負均衡態(tài)勢的主要原因。最后，需要指出的是雖然在MODFLOW模型校準與驗證過程中取得較好的模擬效果，但仍存在一些不足。本研究實測地下水位資料的時間序列較短，觀測井數量較少，僅對兩年以月為單位的6個觀測井的地下水變化進行模擬與驗證，今后仍需進一步增加監(jiān)測井數量和延長監(jiān)測時間，補充校準模型所需實測數據，以降低模型參數的不確定性和提高模擬結果準確性。

3.2 氣候變化對地下水位影響及情景分析

為了評價氣候變化對地下水水位影響，選取2012年末地下水位圖（圖5）作為水位的對比圖。氣候變化情景下，基于已校準的MODFLOW模型，以2012年末觀測井的水位作為模型的初始水位，控制開采量、側向流入、河道滲漏、側向流出不變，分別輸入1997-2013年、2014-2030年的氣候變化條件下降水和潛水蒸發(fā)，進而分析氣候變化對地下水的影響。

1997-2013年歷史氣象條件下及2014-2030年氣候變化條件下模型模擬結果見圖6。模擬結果表明2013年末和2030年末灌區(qū)地下水仍然呈現西南多東北少態(tài)勢，地下水流向由西南部流向東北部。

歷史氣象條件下2013年末及氣候變化條件下2030年末地下水水位與初始水位相比，整體有下降趨勢。本文將灌區(qū)地下水水位低于60 m的區(qū)域劃分為地下水漏斗區(qū)域。從圖6的地下水漏斗分布可以看出，灌區(qū)地下水漏斗從夏莊附近逐漸往西南方向延伸，漏斗區(qū)域面積顯著增大，2013年末漏斗面積小于2030年末漏斗面積。2014-2030年和1997-2013年期間地下水補排量分別見表5和表6，從表中可以看出，1997-2013年歷史氣象條件下地下水補給量比2014-2030年氣候變化條件下地下水補給量多。在模擬期間地下水的主要補給項為降水入滲，排泄項中，開采量為主要的排出項。氣候變化條件下補給項中降水入滲量降低，排泄項中潛水蒸發(fā)量增大。2014-2030年期間與1997-2013年期間相比降水入滲補給地下水量降低3.9%，溫度升高，潛水蒸發(fā)量增加23.8%。

3.3 氣候變化條件下地下水開采方案情景分析

為了深入了解2014-2030年未來氣候情景下開采量增加及開采量減少對灌區(qū)地下水位影響，本研究開展地下水開采方案情景分析。地下水開采量變動20%情景對應的各灌溉區(qū)域（圖2）的開采量見表7。地下水開采量分別增加和減少20%，模擬2030年地下水流場變化情況，不同情景下地下水位見圖7。從預測結果來看，開采量增加或減少20%呈現的結果差異較大。地下水開采量增加20%的情景分析結果表明，預測的2030年地下水位中，后河附近地下水水位在64 m左右，莊巖附近水位在68 m左右，亢村附近水位在58 m左右，夏莊為漏斗中心，水位在29 m左右，灌區(qū)地下水下降明顯，地下水漏斗面積增大。地下水開采量減少20%的情景分析結果表明，[HJ2.08mm]預測的2020年地下水位中，后河附近地下水水位在65 m左右，莊巖附近水位在69 m左右，亢村附近水位在66 m左右，夏莊漏斗中心附近水位在50 m左右。整體來講，開采量減少20%與增加20%情景模擬結果對比發(fā)現地下水位有所上升，地下水漏斗面積明顯減小。為了減緩地下水位持續(xù)下降趨勢，降低地下水漏斗的面積，在未來氣候情景下，建議采取地下水開采量減少20%的方案，進而恢復地下水水位和降低漏斗區(qū)面積。為減少農業(yè)灌溉對地下水的開采，可以采取以下措施：（1）優(yōu)化灌區(qū)的渠系工程，增加地表水灌溉量，減少地下水開采量，推廣先進的灌溉技術，比如噴灌，滴灌，以此來提高農業(yè)灌溉用水效率。（2）對灌區(qū)內的自備井加強管理，廢棄一些違反規(guī)定而開掘的用水井[32-33]。

表7 地下水開采量增加及減少20%后各灌溉區(qū)域開采量

Tab.7

The groundwater exploitation of each irrigation area with a ±20% change[JY，1]m3/d

區(qū)域[]1、2[]3[]4[]5[]6[]7

[BHD]增加20%[]40 670[]141 220[]139 830[]121 660[]15 200[]244 210

[BHDW]減少20%[]27 110[]94 150[]93 220[]81 110[]10 130[]162 800

[HT]

圖7 2030年末地下水等水位線（開采量增加20%及減少20%）

Fig.7

Groundwater level contour at the end of 2030 under scenarios of increase and decrease exploitation by 20%

4 結論

氣候變化與地下水開采是造成灌區(qū)地下水位變化的主要因素，基于MODFLOW模型對華北平原典型灌區(qū)人民勝利渠灌區(qū)地下水對農業(yè)用水和氣候變化的響應進行了模擬，研究發(fā)現在2012-2013年氣象條件和農業(yè)用水影響下地下水處于負均衡狀態(tài)。將驗證后的MODFLOW模型進一步與未來氣候變化情景（RCP4.5情景NorESM1-M模式）耦合模擬發(fā)現，灌區(qū)2030年地下水位持續(xù)下降，地下水漏斗逐漸擴大。情景分析表明地下水開采量減少20%可以保證地下水水位有所上升，漏斗面積減少，而開采量增加20%將導致地下水漏斗面積明顯擴大?；诘叵滤\動數值模擬模型MODFLOW開展未來氣候情景下地下水的響應研究可以為防止地下水位的持續(xù)下降提出調控方案，為灌區(qū)地下水資源持續(xù)利用提供指導，進而防止地下水超量開采引起的生態(tài)環(huán)境問題。本文設定的開采量情景仍存在一定的不足，僅考慮了開采量增加和減少20%的情景，建議以后的研究綜合考慮影響開采量的自然和人為因素，更為準確的估算氣候變化條件下的地下水開采情況。

參考文獻（References）：

[1] 吳愛民，荊繼紅，宋博.略論中國水安全問題與地下水的保障作用[J].地質學報，2016，90（10）：2939-2947.（WU A M，JIN J H，SONG B.Water safety issues of China and ensuring roles of groundwater[J].Acta Geologica Sinca，2016，90（10）：2939-2947.（in Chinese）） DOI：10.3969/j.issn.0001-5717.2016.10.029.

[2] 雷鳴，孔祥斌，張雪靚，等.黃淮海平原區(qū)土地利用變化對地下水資源量變化的影響[J].資源科學，2017，39（6）：1099-1116.（LEI M，KONG X B，ZHANG X L，et al.Land use change and impact on groundwater storage in the Huang-Huai-Hai Plain[J].Resources Science，2017，39（6）：1099-1116.（in Chinese）） DOI：10.18402/resci.2017.06.10.

[3] 張旭光，任春春，王浩然.淺談地下水開采引起的生態(tài)環(huán)境問題[J].中國科技博覽，2008（23）：37-37.（ZHANG X G，REN C C，WANG H R.Discussion on ecological environment caused by groundwater exploitation[J].China Science and Technology Review，2008（23）：37-37 （in Chinese））

[4] 呂晨旭，賈紹鳳，季志恒.近30年來白洋淀流域平原區(qū)地下水位動態(tài)變化及原因分析[J].南水北調與水利科技，2010，8（1）：65-68.（LU C X，JIA S F，JI Z H.Dynamics and causes of groundwater table change in Plain Area of Baiyangdian Basin in last 30 years[J].South-to-North Water Transfers and Water Science & Technology，2010，8（1）：65-68.（in Chinese））? DOI：10.3969/j.issn.1672-1683.2010.01.017.

[5] 薛禹群.論地下水超采與地面沉降[J].地下水，2012，34（6）：1-5.（XUE Y Q.Discussion on groundwater overexploitation and ground settlement[J].Ground Water，2012，34（6）：1-5.（in Chinese）） DOI：10.3969/j.issn.1004-1184.2012.06.001.

[6] 周萬畝，齊全，徐敏，等.地下水超采對農業(yè)灌溉的影響及對策研究[J].地下水，2007，29（4）：17-29.（ZHOU W M，QI Q，XU M，et al.The impact of excessive ground water development on agricultural irrigation and the countermeasures[J].Ground Water，2007，29（4）：17-29.（in Chinese）） DOI：10.3969/j.issn.1004-1184.2007.04.008.

[7] 李振函，張春榮，朱偉.日照市沿海地區(qū)海水入侵現狀與分析[J].水文地質工程地質，2009，36（5）：129-132.（LI Z H，ZHANG C R，ZHU W.The present and analysis of the sea water intrusion in coastland of Rizhao[J].Hydrogeology and Engineering Geology，2009，36（5）：129-132.（in Chinese）） DOI：10.3969/j.issn.1000-3665.2009.05.029.

[8] 尚松浩.基于水旱災害的中國農業(yè)水安全情勢評價[J].華北水利水電大學學報 （自然科學版），2018，39（1）：10-14.（SHANG S H.Assessment on agricultural water security regime in China based on flood and drought disasters[J].Journal of North China university of water resources and electric power （Natural science edition），2018，39（1）：10-14.（in Chinese）） DOI：10.3969/j.issn.1002-5634.2018.01.002.

[9] 劉春蓁，劉志雨，謝正輝.地下水對氣候變化的敏感性研究進展[J].水文，2007，27（2）：1-6.（LIU C Z，LIU Z Y，XIE Z H.Recent advances in research on sensitivity of groundwater to climate changes[J].Journal of China Hydrology，2007，27（2）：1-6.（in Chinese）） DOI：10.3969/j.issn.1000-0852.2007.02.001.

[10] [ZK（#]王蕊，王中根，夏軍，等.氣候變化及南水北調中線工程對白洋淀流域水文的影響[J].水力發(fā)電學報，2012，31（2）：14-19.（WANG R，WANG Z G，XIA J，et al.Effects of climate change and the south-north water transfer project on hydrology in Baiyangdian basin[J].Journal of Hydroelectric Engineering，2012，31（2）：14-19.（in Chinese））

[11] 陳皓銳，高占義，王少麗，等.基于 Modflow 的潛水位對氣候變化和人類活動改變的響應[J].水利學報，2012，43（3）：344-353.（CHEN H R，GAO Z Y，WANG S L，et al.Modeling on impacts of climate change and human activities variability on the shallow groundwater level using modflow[J].Journal of Hydraulic Engineering，2012，43（3）：344-353.（in Chinese）） DOI：10.13243/j.cnki.slxb.2012.03.002.

[12] SCIBEK J，ALLEN D M.Modeled impacts of predicted climate change on recharge and groundwater levels[J].Water Resources Research，2006，42（11）.DOI：10.1029/2005WR004742.

[13] KIRSHEN P H.Potential impacts of global warming on groundwater in eastern Massachusetts[J].Journal of Water Resources Planning and Management，2002，128（3）：216-226.DOI：10.1061/（ASCE）0733-9496（2002）128：3（216）.

[14] OUDE ESSINK GHP，VAN BAAREN E S，DE LOUW P G B.Effects of climate change on coastal groundwater systems：A modeling study in the Netherlands[J].Water Resources Research，2010，46（10）.DOI：10.1029/2009WR008719.

[15] 廖梓龍，龍胤慧，劉華琳，等.氣候變化與人類活動對包頭市地下水位的影響[J].干旱區(qū)研究，2014，1：023.（LIAO Z L，LONG Y H，LIU H L，et al.Influence of climate change and human activities on groundwater level in Baotou city[J].Arid Zone Research，2014，1：023.（in Chinese））

[16] 劉路廣，崔遠來，羅玉峰.基于 Modflow 的灌區(qū)地下水管理策略——以柳園口灌區(qū)為例[J].武漢大學學報：工學版，2010，43（1）：25-29.（LIU L G，CUI Y L，LUO Y F.Evaluating groundwater management strategies of irrigation system based on modflow：A case study of Liuyuankou irrigation system[J].Engineering Journal of Wuhan University.2010，43（1）：25-29.（in Chinese））

[17] 王電龍，張光輝，馮慧敏，等.氣候變化下華北平原井灌區(qū)糧食生產地下水保障能力[J].海河水利，2016（02）：34-41.（WANG D L，ZHANG G H，FENG H M，et al.Grain production groundwater ensuring extant of north China plain well irrigation area in climate change scenarios[J].Haihe Water Resources，2016（02）：34-41.（in Chinese））? DOI：10.3969/j.issn.1004-7328.2016.02.014.

[18] 林艷竹.基于數值模擬的長期氣候條件下華北平原地下水響應研究[D].北京：中國地質大學 （北京），2015.（LIN Y Z.Study of the groundwater response in the north China plain under background of long-term climate change based on numerical simulation[D].Beijing：China University of Geosciences（Beijing），2015.（in Chinese））

[19] 王利書，悅琳琳，唐澤軍，等.氣候變化和農業(yè)發(fā)展對石羊河流域地下水位的影響[J].農業(yè)機械學報，2014，45（1）：121-128.（WANG L S，YUE L L，TANG Z J，et al.Influence of climate change and agricultural development on groundwater level in Shiyang river basin[J].Transactions of The Chinese Society of Agricultural Machinery，2014，45（1）：121-128.（in Chinese）） DOI：10.6041/j.issn.1000-1298.2014.01.020

[20] 王焰新，馬騰，郭清海，等.地下水與環(huán)境變化研究[J].地學前緣，2005，12（s）：14-21.（WANG Y X，MA T，GUO Q H，et al.Groundwater and environment change[J].Earth Science Frontiers，2005，12（s）：14-21.（in Chinese））? DOI：10.3321/j.issn：1005-2321.2005.z1.003.

[21] BRADLEY R.1000 years of climate change[J].Science，2000，288（5470）：1353-1355 DOI：10.1126/science.288.5470.1353.

[22] 陳兆波，董文，霍治國，等.中國農業(yè)應對氣候變化關鍵技術研究進展及發(fā)展方向[J].中國農業(yè)科學，2013，46（15）：3097-3104.（CHEN Z B，DONG W，HUO Z G，et al.Research advances and the development direction of the key technique for agriculture addressing climate change in China[J].Scientia Agricultura Sinica，2013，46（15）：3097-3104.（in Chinese）） DOI：10.3864/j.issn.0578-1752.2013.15.003.

[23] 張嘉星，齊學斌，李平，等.人民勝利渠灌區(qū)適宜井渠用水比研究[J].灌溉排水學報，2017，36（2）：58-63.（ZHANG J X，QI X B，LI P，et al.Appropriate well-canal irrigation proportion in irrigation district based on MODFLOW and GIS[J].Journal of irrigation and drainage，2017，36（2）：58-63.（in Chinese））

[24] 姚遙，羅勇，黃建斌.8 個 CMIP5 模式對中國極端氣溫的模擬和預估[J].氣候變化研究進展，2012，8（4）：250-256.（YAO Y，LUO Y，HUANG J B.Evaluation and projection of temperature extremes over China based on 8 modeling data from CMIP5[J].Advances in Climate Change Research，2012，8（4）：250-256.（in Chinese））? DOI：10.3969/j.issn.1673-1717.2012.04.003.

[25] 閆旖君，徐建新，肖恒.2021—2050 年河南省夏玉米凈灌溉需水量對氣候變化的響應[J].氣候變化研究進展，2017，13（2）：138-148.（YAN Y J，XU J X，XIAO H.Summer maize net irrigation water requirement to climate change in Henan province in China during 2021-2050[J].Advances in Climate Change Research，2017，13（2）：138-148.（in Chinese）） DOI：10.12006/j.issn.1673-1719.2016.168.[HJ1.97mm]

[26] 馬黎華.華北區(qū)ENSO年夏冬兩季月降水量的時間序列模擬分析[J].華北水利水電大學學報（自然科學版）.2018.39（2）：69-73.（MA L H，Analysis of the time series of precipitation in the summer and winter months of ENSO in north China[J].Journal of North China University of water resources and electric power （Natural science edition），2018.39（2）：69-73.（in Chinese））

[27] 肖楠森，林鳳勳，華致潔，等.豫北人民勝利渠引黃灌區(qū)的水文地質特征及其與土壤鹽堿化的關系[J].南京大學學報（自然科學），1963（13）：43-60.（XIAO N S，LIN F X，HUA Z J，et al.The hydrogeological characteristics of the Peoples′ Victory canal in the Yellow River in the north of Henan and its relationship with soil salinization[J].Journal of Nanjing University（Natural Sciences），1963（13）：43-60.（in Chinese））

[28] 袁錫泰，龔曉潔，余長合，等.黃河沖積平原灌溉回滲系數研究[J].地下水，2012，34（5）：181-182.（YUAN X T，GONG X J，YU C H，et al.Research on irrigation retardation coefficient of alluvial plain in Yellow River[J].Ground water，2012，34（5）.（in Chinese）） DOI：10.3969/j.issn.1004-1184.2012.05.069

[29] 張巖松，賈海峰.北京典型平原區(qū)地下水環(huán)境模擬及情景分析[J].清華大學學報：自然科學版，2008，48（9）：1436-1440.（ZHANG Y S，JIA H F.Groundwater modeling and scenario analysis in the Beijing typical plain area[J].Journal of Tsinghua University （Science and Technology），2008，48（9）：1436-1440.（in Chinese）） DOI：10.3321/j.issn：1000-0054.2008.09.015.

[30] 張雪芹，彭莉莉，林朝暉.未來不同排放情景下氣候變化預估研究進展[J].地球科學進展，2008，23（2）：174-185.（ZHANG X Q，PENG L L，LIN Z H.Progress on the projections of future climate change with various emission scenarios[J].Advances In Earth Science，2008，23（2）：174-185.（in Chinese））

[31] [ZK（#]王建紅，余啟明，楊俊倉.基于ArcGIS和Visual MODFLOW的黑河流域中游平原區(qū)地下水流數值模擬與預測[J].安全與環(huán)境工程，2016，23（1）：80-87.（WANG J H，YU Q M，YANG J C.Numerical simulation and forecast of groundwater in the middle plains of Heihe River based on ArcGIS and Visual MODFLOW[J].Safety and Environmental Engineering，2016，23（1）：80-87.（in Chinese）） DOI：10.13578/j.cnki.issn.1671-1556.2016.01.016.

[32] 魯東霞，張哲，丁颯，等.河南糧食核心區(qū)建設規(guī)劃地下水環(huán)境影響評價[J].農業(yè)環(huán)境與發(fā)展，2009，26（4）：59-64.（LU D X，ZHANG Z，DING S，et al，Groundwater environment impact assessment of Henan grain core area construction planning[J].Agro-Environment and Development，2009，26（4）：59-64.（in Chinese）） DOI：10.3969/j.issn.1005-4944.2009.04.018.

[33] 太紅禮.寧夏固原北川地下水資源開發(fā)利用探討[J].水利規(guī)劃與設計，2013 （1）：23-24.（TAI H L.Exploitation and utilization of Beichuang groundwater resources in Guyuan，Ningxia[J].Water Resources Planning and Design，2013（1）：23-24.（in Chinese）） DOI：10.3969/j.issn.1672-2469.2013.01.009.

[34] 趙雨婷.黃河中下游典型灌區(qū)地下水位時空演變特征與驅動機制[D].鄭州：華北水利水電大學，2016.（ZHAO Y T.Temporal and spatial characteristics and driving mechanism of ground typical middle and lower reaches of the Yellow River irrigation district[D].Zhengzhou：North China University of Water Resources and Electric Power.（in Chinese））