李元杰 姜新慧 陳軍 楊志巖 王貴玲 丁慧君

摘要：在查明內(nèi)蒙古巴彥淖爾市臨河區(qū)水文地質(zhì)條件、地下水開發(fā)及引發(fā)的環(huán)境地質(zhì)問題的基礎(chǔ)上，用均衡法及數(shù)值模擬法進(jìn)行地下水資源量計(jì)算。根據(jù)當(dāng)?shù)厮Y源開發(fā)利用方案，利用模型進(jìn)行研究區(qū)2013年-2033年地下水流場(chǎng)演化模擬，結(jié)果表明該方案下研究區(qū)地下水降落漏斗不再繼續(xù)擴(kuò)大，地下水位有所回升，對(duì)該區(qū)水資源管理有指導(dǎo)意義。

關(guān)鍵詞：巴彥淖爾市臨河區(qū)；地下水資源；開發(fā)利用；地下水漏斗；水均衡法；數(shù)值模擬

中圖分類號(hào)：P641 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)：1672-1683（2016）04-0036-06

Abstract：After the identification of the hydrogeological conditions and the environmental geological problems caused by groundwater exploitation，the water-balance and numerical simulation methods were used to evaluate the amount of groundwater in Linhe district of Bayan Nur City，Inner Mongolia.A model was utilized to simulate the change of groundwater flow field in the study area from 2013 to 2033 according to the local water resources exploration.and utility scenario.and the results showed that the groundwater drawdown funnel will cease to expand and the groundwater level will rise.This provides a meaningful guidance for the water resources management in the study area.

Key words：Linhe district of Bayan Nur City；groundwater resources；development and utilization；groundwater funnel；water-balance method；numerical simulation

隨著經(jīng)濟(jì)社會(huì)發(fā)展和全球環(huán)境變化，我國(guó)水短缺、水污染、水生態(tài)、水災(zāi)害、水管理5個(gè)問題復(fù)雜交叉，已成為我國(guó)可持續(xù)發(fā)展的主要瓶頸[1-4]。正確處理有限的水資源與經(jīng)濟(jì)發(fā)展、生態(tài)-環(huán)境保護(hù)之間的矛盾已成為重要和緊迫的國(guó)家戰(zhàn)略要求，其關(guān)鍵在于科學(xué)的水資源管理[5-7]。目前的管理模式多基于“三條紅線”、生態(tài)環(huán)境綜合治理、地表水與地下水聯(lián)合調(diào)度和流域水資源集成管理[8]。運(yùn)用先進(jìn)的地下水模擬軟件，建立地下水流數(shù)值模型，預(yù)測(cè)各種條件下地下水流場(chǎng)變化，已成為區(qū)域地下水管理的有效方法[9-12]。本文運(yùn)用GMS軟件，對(duì)現(xiàn)狀條件下巴彥淖爾市臨河區(qū)地下水資源進(jìn)行數(shù)值模擬，根據(jù)當(dāng)?shù)厮Y源開發(fā)利用方案，進(jìn)行地下水流場(chǎng)預(yù)報(bào)，為該區(qū)水資源管理提供依據(jù)。

1 研究區(qū)概況

1.1 自然地理概況

內(nèi)蒙古自治區(qū)巴彥淖爾市臨河區(qū)面積約1 400 km²，主體地貌為黃河沖湖積平原，地勢(shì)平坦開闊，地面高程在1 034～1 044 m，微向東北傾斜，地形自然坡度0.15‰～0.18‰。本區(qū)屬中溫帶干旱—半干旱大陸性氣候，多年平均降水量149.04 mm，多集中在7月-9月，多年平均蒸發(fā)量為2 183 mm，主要集中在4月-8月。黃河是唯一的長(zhǎng)年地表水流，自西向東沿研究區(qū)南部邊緣流過。

1.2 水文地質(zhì)條件

研究區(qū)位于河套斷陷盆地臨河凹陷區(qū)的南部邊緣，河套斷陷盆地沉積了厚度較大的第四系沉積物，主要有更新統(tǒng)沖、洪積物以及全新統(tǒng)湖積物、風(fēng)積物。研究區(qū)屬于雙層結(jié)構(gòu)含水層，上部為淺層潛水-微承壓含水層，上更新統(tǒng)至全新統(tǒng)，巖性以沖湖積礫砂、中粗砂、中細(xì)砂、細(xì)粉砂為主；下部為承壓含水層，中更新統(tǒng)下段，巖性為細(xì)砂、粉細(xì)砂。下部承壓水埋藏深，一般在300～500 m以下，且其水質(zhì)差（礦化度大于3～5 g/L）、水量小，開采條件差[13-15]，因此本次工作只研究上部潛水-微承壓水含水層（見圖1）。

1.3 地下水開采引發(fā)的環(huán)境地質(zhì)問題

研究區(qū)位于河套平原黃灌區(qū)，農(nóng)業(yè)灌溉以引黃水為主，部分農(nóng)業(yè)用水、生活用水及工業(yè)生產(chǎn)用水主要依靠開采地下水。隨著城鎮(zhèn)化建設(shè)和工業(yè)的發(fā)展，對(duì)地下水的開采強(qiáng)度逐年加大，引發(fā)的環(huán)境地質(zhì)問題主要是地下水降落漏斗的不斷擴(kuò)大。根據(jù)收集的巴彥淖爾市臨河區(qū)1982年、1989年、1998年、2004年4個(gè)時(shí)期的等水位線、降深圖，以及2013年進(jìn)行的該區(qū)地下水位統(tǒng)測(cè)資料，得出不同時(shí)期降落漏斗中心點(diǎn)坐標(biāo)，漏斗中心水位標(biāo)高、漏斗面積等參數(shù)（表1）地下水位降落漏斗發(fā)展變化圖及剖面圖（圖2）。

通過對(duì)臨河區(qū)城市發(fā)展進(jìn)程的了解及對(duì)圖2的分析，城區(qū)降落漏斗的形成與發(fā)展大致經(jīng)歷了以下三個(gè)時(shí)期。

第一時(shí)期：20世紀(jì)70年代末漏斗形成初期。即在城區(qū)第一個(gè)水源地勘查之前，地下水位受城區(qū)廣泛分布的手壓井和企業(yè)自備井開采影響，開采量約219.28萬m³/a。

第二時(shí)期：中心水位穩(wěn)定，漏斗面積擴(kuò)大期。隨著城市規(guī)模的不斷擴(kuò)大，漏斗中心向東北方向（開發(fā)區(qū)）偏移，分布面積76.0 km²，水位降幅變化不明顯，主要受第一水源地及城區(qū)企業(yè)自備井開采影響，開采量為515.17萬m³/a。

第三時(shí)期：中心水位下降，漏斗面積擴(kuò)大期。至2013年，漏斗面積達(dá)153.72 km²，中心水位標(biāo)高1 032.49 m。地下水位受第一、第二水源地，城區(qū)、開發(fā)區(qū)企業(yè)自備井開采影響，開采量為2 023.35萬m³/a。

降落漏斗演化期的分界點(diǎn)是隨著經(jīng)濟(jì)發(fā)展和城市集中水源地的建設(shè)而出現(xiàn)的，1982年-1998年漏斗中心大致處于相同區(qū)域，啤酒廠和一水廠（王明圪旦供水水源地）是該區(qū)地下水的主要開采點(diǎn)，城區(qū)居民院內(nèi)有部分手壓井。2004年隨著巴彥淖爾市“十一五”工業(yè)項(xiàng)目規(guī)劃的實(shí)施，東部開發(fā)區(qū)一批企業(yè)配備自備井大量開采地下水，降落漏斗向東北方向偏移。2008年后隨著城市建設(shè)的不斷拓展，一水廠停采，沿永濟(jì)渠布設(shè)新的水源井，漏斗范圍擴(kuò)大，中心轉(zhuǎn)移至城區(qū)西南部。

2 地下水?dāng)?shù)值模型的建立

2.1 均衡分析

在模型建立之前先用水均衡法進(jìn)行研究區(qū)水資源量計(jì)算，以便檢驗(yàn)?zāi)M結(jié)果的準(zhǔn)確性、客觀性。均衡方程見式（1），即任一均衡區(qū)任一時(shí)段的補(bǔ)給量和排泄量之差等于孔隙水體積的變化量。

均衡區(qū)面積約為1 400 km²，按照地下水系統(tǒng)將其劃分為永濟(jì)灌域和黃河漫灘兩部分。本次計(jì)算選取一個(gè)水文年（2012年5月-2013年5月）作為均衡期。依托在巴彥淖爾市杭錦后旗建立的包氣帶水鹽運(yùn)移試驗(yàn)場(chǎng)[16]，通過長(zhǎng)期試驗(yàn)觀測(cè)，獲取降水入滲系數(shù)、蒸發(fā)系數(shù)及灌溉入滲系數(shù)；通過抽水試驗(yàn)、土工試驗(yàn)等，確定含水層的滲透系數(shù)及給水度。均衡區(qū)內(nèi)地下水資源計(jì)算結(jié)果見表2。

研究區(qū)均衡期降水量為近20年最高值，為豐水年，地下水總補(bǔ)給量為31 041.29萬m³/a，其中引黃田間灌溉入滲量占總補(bǔ)給量的39%，為主要補(bǔ)給項(xiàng)，渠道入滲量、大氣降水入滲量分別占總補(bǔ)給量的29%、28%。排干排泄量、潛水蒸發(fā)和地下水開采是該均衡區(qū)的主要排泄項(xiàng)，排干排泄量潛水蒸發(fā)、地下水開采分別占總排泄量的41%、31%、24%。均衡期地下水處于正均衡狀態(tài)，均衡差為901.05萬 m³/a。

依據(jù)長(zhǎng)觀井水位變幅數(shù)據(jù)和含水層厚度，采用泰森多邊形分區(qū)計(jì)算得出地下水儲(chǔ)存變化量為931.6萬m³/a，與均衡期補(bǔ)排差進(jìn)行比較，相對(duì)誤差為3.35%，本次均衡計(jì)算的結(jié)果較為合理。

2.2 數(shù)值模擬

自20世紀(jì)70年代末，數(shù)值模擬逐漸成為地下水資源評(píng)價(jià)與管理的主要方法之一[17]。根據(jù)研究區(qū)的水文地質(zhì)條件將南部的黃河概化為定水頭邊界，其余為流量邊界；潛水面為上邊界，下部的淤泥質(zhì)粘土為隔水底部邊界，將地下水流系統(tǒng)概化為非均質(zhì)、各向同性、空間三維、非穩(wěn)定流系統(tǒng)，用如下微分方程的定解問題來描述：

本次模擬期為2012年5月1日至2013年4月30日，以2012年5月1日統(tǒng)測(cè)水位為模型的初始流場(chǎng)，每月為一個(gè)應(yīng)力期，時(shí)間步長(zhǎng)為1 d。各流量邊界的水力梯度，依據(jù)研究區(qū)地下水動(dòng)態(tài)監(jiān)測(cè)資料按時(shí)段賦值，大氣降水入滲補(bǔ)給、蒸發(fā)排泄源匯項(xiàng)依據(jù)臨河氣象站提供的資料，分區(qū)分時(shí)段賦值。

調(diào)用GMS地下水?dāng)?shù)值模擬計(jì)算模塊，不斷調(diào)整參數(shù)及源匯項(xiàng)，得到模擬期末的地下水流場(chǎng)，同時(shí)選取研究區(qū)內(nèi)的觀測(cè)孔在模擬期的水位動(dòng)態(tài)資料進(jìn)行擬合，使模型計(jì)算值逼近觀測(cè)值，從而確定模型各分區(qū)的水文地質(zhì)參數(shù)，最終得出研究區(qū)地下水系統(tǒng)水量均衡結(jié)果（見表3），流場(chǎng)模擬結(jié)果及監(jiān)測(cè)孔擬合情況（見圖3）。

由此可見，研究區(qū)在模擬期內(nèi)的水量補(bǔ)排基本平衡，計(jì)算流場(chǎng)基本上反映了地下水流動(dòng)的趨勢(shì)和規(guī)律。各觀測(cè)孔在各個(gè)時(shí)段的計(jì)算水位與觀測(cè)水位擬合程度較好，擬合差小于0.5 m，說明含水層結(jié)構(gòu)、邊界條件概化、水文地質(zhì)參數(shù)及源匯項(xiàng)的選取是合理的，建立的數(shù)值模型較為真實(shí)地刻畫了研究區(qū)地下水系統(tǒng)的特征，仿真性強(qiáng)，可利用該模型進(jìn)行地下水流場(chǎng)預(yù)報(bào)。

3 地下水流場(chǎng)預(yù)報(bào)

根據(jù)巴彥淖爾市臨河區(qū)的水資源開發(fā)利用方案，本次模擬預(yù)報(bào)時(shí)段分為2013年5月-2018年5月、2018年5月-2023年5月、2023年5月-2033年5月。以2013年5月流場(chǎng)作為初始流場(chǎng)，預(yù)報(bào)模型的結(jié)構(gòu)參數(shù)、水文地質(zhì)參數(shù)均采用識(shí)別后的參數(shù)，源匯項(xiàng)及參數(shù)的確定如下。

補(bǔ)給項(xiàng)中降水入滲、黃河側(cè)滲、地下水側(cè)向流入和井灌回歸量按照多年平均值給定。2012年5月自治區(qū)啟動(dòng)四個(gè)千萬畝節(jié)水灌溉工程，逐步開展節(jié)水改造工作，預(yù)測(cè)2018年、2023年、2033年灌溉水有效利用系數(shù)將由41%分別提升至43.5%、46%、50%；渠系水有效利用率由48%提升至52%、56%、62%。預(yù)報(bào)期間引黃田間入滲量、渠道入滲量的確定以多年平均值為基礎(chǔ)值，按照相應(yīng)的系數(shù)進(jìn)行調(diào)整。

排泄項(xiàng)中潛水蒸發(fā)、湖泊、海子蒸發(fā)和側(cè)向排泄按多年平均值給定。預(yù)報(bào)期間，規(guī)劃運(yùn)行黃河水廠引水工程滿足城區(qū)生活及企業(yè)用水需求，2013年5月起封閉企事業(yè)自備井，預(yù)測(cè)2018年自備井全部關(guān)停，2023年部分水源井可壓減現(xiàn)有開采量的一半，2033年原有的水源井全部關(guān)停，地下水僅選擇水質(zhì)較好的地段進(jìn)行開采。預(yù)報(bào)期間地下水開采量為多年平均開采量減去壓采量；排干年均泄水量=多年排干泄水量平均值×（補(bǔ)給項(xiàng)多年平均值-預(yù)報(bào)期內(nèi)年均補(bǔ)給項(xiàng)值）÷補(bǔ)給項(xiàng)多年平均值。

依據(jù)以上源匯項(xiàng)資料以及參數(shù)序列，進(jìn)行地下水流場(chǎng)模擬，模擬結(jié)果見表4、圖4。

由圖4可見2013年-2023年地下水水位呈緩慢回升趨勢(shì)，市中心的水位標(biāo)高為1 033 m的等水位線所封閉的漏斗中心面積呈逐漸縮小的趨勢(shì)；2023年-2033年，降落漏斗面積呈繼續(xù)縮小的趨勢(shì)，地下水流場(chǎng)整體沒有發(fā)生大的變化。

由此可見，研究區(qū)水資源開發(fā)利用方案中的提高灌溉水利用系數(shù)和封停自備井等措施可以保證研究區(qū)地下水降落漏斗不再繼續(xù)擴(kuò)大，地下水位略有回升，地下水流場(chǎng)不發(fā)生較大變化，模擬結(jié)果比較理想。采用該方案研究區(qū)地下水資源的供需矛盾將有所緩解，若要顯著改善當(dāng)?shù)厮Y源問題，保證地下水資源的可持續(xù)利用，還應(yīng)積極調(diào)整區(qū)內(nèi)產(chǎn)業(yè)結(jié)構(gòu)及布局，提高廢、污水處理率及用水效率，引進(jìn)先進(jìn)的節(jié)水灌溉技術(shù)[18-20]，加強(qiáng)地表水與地下水優(yōu)化配置利用研究，積極摸索水資源高效管理的方式。

4 結(jié)論

（1）分別采用數(shù)值法和水均衡法對(duì)研究區(qū)地下水資源量進(jìn)行計(jì)算，兩種方法較好地評(píng)價(jià)了研究區(qū)的地下水資源狀況，模擬期內(nèi)該區(qū)淺層地下水總補(bǔ)給量為3.12億m³/a，總排泄量為3.02億m³/a，補(bǔ)排差為0.09億m³/a，總體處于正均衡狀態(tài)。

（2）通過對(duì)比模擬水位與觀測(cè)水位，二者擬合程度良好，說明含水層結(jié)構(gòu)、邊界條件概化、水文地質(zhì)參數(shù)及源匯項(xiàng)的選取是合理的，可利用該模型進(jìn)行地下水流場(chǎng)預(yù)報(bào)。

（3）根據(jù)當(dāng)?shù)厮Y源開發(fā)利用方案，利用數(shù)值模型進(jìn)行地下水流場(chǎng)預(yù)報(bào)，結(jié)果顯示，在提高灌溉水利用系數(shù)和封停自備井情況下，2018年-2033年該區(qū)地下水降落漏斗中心面積呈逐漸縮小的趨勢(shì)，地下水位略有回升，若要顯著改善當(dāng)?shù)厮Y源問題，還應(yīng)積極探索水資源高效管理的方式。

參考文獻(xiàn)（References）：

[1] 夏軍，翟金良，占車生.我國(guó)水資源研究與發(fā)展的若干思考[J].地球科學(xué)進(jìn)展，2011，26（9）：905-915.（XIA Jun，ZHAI Jin-liang，ZHAN Che-sheng.Some reflections on the research and development of water resources in China[J].Advances in Earth Science，2011，26（9）：905-915.（in Chinese））

[2] 汪黨獻(xiàn)，王浩，馬靜.中國(guó)區(qū)域發(fā)展的水資源支撐能力[J].水利學(xué)報(bào)，2000，（11）：21-27.（WANG Dang-xian，WANG Hao，MA Jing.Water resources supporting capacity for regional development in China[J].Journal of Hydraulic Engineering，2000，（11）：21-27.（in Chinese））

[3] 李九一，李麗娟.中國(guó)水資源對(duì)區(qū)域社會(huì)經(jīng)濟(jì)發(fā)展的支撐能力[J].地理學(xué)報(bào)，2012，67（3）：410-419.（LI Jiu-yi，LI Li-juan.Water resources supporting capacity toregional socio-economic development of China[J].Acta Geographica Sinica，2012，67（3）：410-419.（in Chinese））

[4] 姜蓓蕾，耿雷華，徐澎波，等.我國(guó)水資源管理實(shí)踐發(fā)展及管理模式演變趨勢(shì)淺析[J].中國(guó)農(nóng)村水利水電，2011 （10）：66-69.（JIANG Bei-lei，GENG Lei-hua，XU Peng-bo，et al.A Tentative analysis of water resources management practice and the development tendency of the water resources management mode[J].China Rural Water and Hydropower，2011 （10）：66-69.（in Chinese））

[5] 趙海莉，張志強(qiáng)，趙銳鋒.黑河流域水資源管理制度歷史變遷及其啟示[J].干旱區(qū)地理，2014，37（1）：45-55.（ZHAO Hai-li， ZHANG Zhi-jiang， ZHAO Rui-feng.Historical changes and the enlightenments of water resources management system of Heihe River Basin[J].Arid Land Geography，2014，37（1）：45-55.（in Chinese））

[6] 薛冰，宋新山，嚴(yán)登華.基于系統(tǒng)動(dòng)力學(xué)的天津市水資源模擬及預(yù)測(cè)[J].南水北調(diào)與水利科技，2011，9（6）：43-47.（XUE Bing，SONG Xin-shan，YAN Deng-hua.Simulation and prediction of water resources carrying capacity based on a system dynamic model in Tianjin[J].South to North Water Transfers and Water Science & Technology，2011，9（6）：43-47.（in Chinese））

[7] 張鴻義.中國(guó)干旱區(qū)地下水資源及開發(fā)潛力分析[J].干旱區(qū)研究，2009，26（ 2）：149-161.（Groundwater resources and their exploitation potential in arid area in China[J].Arid Zone Research，2009，26（ 2）：149-161.（in Chinese））

[8] 孫棟元，金彥兆，李元紅，等.干旱內(nèi)陸河水資源管理模式研究[J].中國(guó)農(nóng)村水利水電，2015 （1）：80-84，89.（SUN Dong-yuan， JIN Yan-zhao， LI Yuan-hong，et al.Patterns of water resources management in Arid Inland River Basins[J].China Rural Water and Hydropower，2015 （1）：80-84，89.（in Chinese））

[9] 孫從軍，韓振波，趙振，等.地下水?dāng)?shù)值模擬的研究與應(yīng)用進(jìn)展[J].環(huán)境工程，2013，31（ 5）：9-13，17.（SUN Cong-jun，HAN Zhen-bo，ZHAO Zhen，et al.Advances in research and application of groundwater numerical simulation[J].Environmental Engineering，2013，31（ 5）：9-13，17.（in Chinese））

[10] 紀(jì)媛媛，周金龍，楊廣焱.GMS在我國(guó)地下水資源評(píng)價(jià)與管理中的應(yīng)用[J].地下水，2013，35（2）：76-79.（JI Yuan-yuan，ZHOU Jin-long，YANG Guang-yan.GMS application for groundwater resources evaluation and management in China[J].Ground Water，2013，35（2）：76-79.（in Chinese））

[11] 剛什婷，鄧英爾.基于GMS在地下水資源評(píng)價(jià)與管理中的應(yīng)用綜述[J].地下水，2015，37（2）：33-36.（GANG Shen-ting，DENG Ying-er.Comprehensive summary of GMS application for groundwater resources evaluation and management[J].Ground Water，2015，37（2）：33-36.（in Chinese））

[12] 張奇.數(shù)值模型在地下水管理中的應(yīng)用[J].水文地質(zhì)工程地質(zhì)，2003，（6）：72-79，88.（ZHANG Qi.Application of numerical models in groundwater management[J].Hydrogeology & Engineering Geology，2003，（6）：72-79，88.（in Chinese））

[13] 李元杰，楊志巖，丁慧君，等.內(nèi)蒙古自治區(qū)巴彥淖爾市臨河區(qū)地下水資源綜合評(píng)價(jià)（1：5萬）[R].內(nèi)蒙古自治區(qū)地質(zhì)環(huán)境監(jiān)測(cè)院，2014.（LI Yuan-jie，YANG Zhi-yan，DING Hui-jun，et al.Comprehensive evaluation of groundwater resources at Linhe District in Bayan Nur City，Inner Mongolia （1：50000）[R].Inner Mongolia Institute of Geological Environmental Monitoring，2014.（in Chinese））

[14] 內(nèi)蒙古自治區(qū)水事監(jiān)理服務(wù)中心.《內(nèi)蒙古巴彥淖爾市臨河區(qū)水資源調(diào)查與評(píng)價(jià)報(bào)告》[R].2007.

（Management and Service Center of Water Resources in Inner Mongolia.Survey and evaluation report of groundwater resources at Linhe District in Bayan Nur City，Inner Mongolia[R].2007.（in Chinese））

[15] 中國(guó)地質(zhì)科學(xué)院水文地質(zhì)環(huán)境地質(zhì)研究所，內(nèi)蒙古地質(zhì)環(huán)境監(jiān)測(cè)院.河套盆地地下水資源及其環(huán)境問題調(diào)查評(píng)價(jià)成果報(bào)告（1：25萬）[R].2013.（The Institute of Hydrogeology and Environment Geology，Inner Mongolia Institute of Geological Environmental Monitoring.Investigation and evaluation results report of groundwater resources and environment in Hetao Basin（1：250000）[R].2013.（in Chinese））

[16] 內(nèi)蒙古地質(zhì)局第一區(qū)調(diào)隊(duì).內(nèi)蒙古自治區(qū)杭錦后旗工業(yè)園區(qū)供水水源地水文地質(zhì)詳查報(bào)告（1：5萬）[R].2008.（Inner Mongolia First Regional Survey team.Hydrogeological investigation report at Industrial Park of Hangjinhou Banner in Inner Mongolia （1：50000）[R].2008.（in Chinese））

[17] 王金生，王長(zhǎng)申，滕彥國(guó).地下水可持續(xù)開采量評(píng)價(jià)方法綜述[J].水利學(xué)報(bào)，2006，37（5）：525-533.（WANG Jin-sheng，WANG Chang-shen，TENG Yan-guo.Review on assessment methods of groundwater sustainable yield[J].Journal of Hydraulic Engineering，2006，37（5）：525-533.（in Chinese））

[18] 曾浩，張中旺，孫小舟，等.湖北漢江流域水資源承載力研究[J].南水北調(diào)與水利科技，2013，11（4）：22-25，30.（ZHENG Hao，ZHANG Zhong-wang，SUN Xiao-zhou.et al.Study of water resources carrying capacity in the Hanjiang River Basin of Hubei[J].South to North Water Transfers and Water Science & Technology，2013，11（ 4）：22-25，30.（ in Chinese） ）

[19] 秦大庸，羅翔宇，陳曉軍，等.西北干旱區(qū)水資源開發(fā)利用潛力分析[J].自然資源學(xué)報(bào)，2004，19（2）：143-150.（QIN Da-yong，LUO Xiang-yu，CHEN Xiao-jun.et al.Analysis of Potentials of water resources development and utilization in arid zone of Northwest China[J].Journal of Natural Resources，2004，19（2）：143-150.（in Chinese））

[20] 彭芳，張新林.河套灌區(qū)節(jié)水灌溉潛力與水資源可持續(xù)利用對(duì)策[J].中國(guó)農(nóng)村水利水電，2005（4）：5-7.（PENG Fang，ZHANG Xin-lin.Potentiality of water-saving irrigation and countermeasures of sustainable utilization of water resources in Hetao irrigation region[J].China Rural Water and Hydropower，2005（4）：5-7.（in Chinese））