基于SD模型的博斯騰湖水量平衡系統分析與仿真研究

伊麗努爾·阿力甫江，玉素甫江·如素力

（1.新疆師范大學 地理科學與旅游學院，烏魯木齊830054；2.新疆干旱區湖泊環境與資源重點實驗室，烏魯木齊830054；3.中國水利水電科學研究院 流域水循環模擬與調控國家重點實驗室，北京100038）

水是地球生命系統中最基礎的自然資源，也是關系一個國家或地區生態安全、糧食安全及社會經濟健康和可持續發展的重要戰略性資源。干旱地區水資源合理利用是一項世界性的重要課題。其中，水量平衡研究既是水資源評價方法的基礎，也是農田節水調控，農業合理用水的依據［1］。隨著人類活動的加劇和氣候變化的綜合影響，大量湖泊生態環境退化、濕地旱化，湖區沼澤面積逐年縮小，并且湖泊水量失衡［2－4］。已有研究表明近50年來國內湖泊變化顯著［5－6］。由于受氣候暖干化的影響，位于干旱與半干旱地區的博斯騰湖，近幾十年來水位下降，環境質量逐步惡化。湖泊水位由有記錄的2002年最高水位1 049.39m下降到2013年的1 045m左右，已達湖水生態紅線［7］。由于社會、經濟、環境系統及其關系的復雜性和動態性，以及博斯騰湖生態治理和水資源開發利用的不確定性和多目標性，使得各種水資源政策對博斯騰湖水安全格局的影響缺乏定量化的手段和方法［8－10］。

博斯騰湖及其周圍生態環境問題集中表現在博斯騰湖水量平衡和水位變化上。國內學者對博斯騰湖水量平衡的研究始于20世紀50年代。1958年，中國科學院新疆綜合考察隊對博斯騰湖的水量平衡進行過較為詳細的研究［11］；后來，許多學者利用水文、氣象實測資料，運用水量平衡方程定量分析計算博斯騰湖水量平衡要素及其與氣象因素的關系［12－15］。

雖然國內對博斯騰湖水量平衡的研究很多［16－22］，但大多數是在研究水量平衡時用湖泊水量平衡方程式計算，沒有考慮博斯騰湖上游和下游的社會經濟要素對水量平衡的影響以及它們之間的雙向反饋關系。這樣，由湖盆地貌形態決定的湖水位—庫容—水面面積、湖泊水量平衡要素及周圍社會經濟要素等之間的非線性關系不能表達出來，會影響水量平衡模擬結果的正確性。因此本文會在前人研究的基礎上，采用系統動力學方法（SD，System Dynamics），通過構建博斯騰湖水量平衡模型及對湖水未來30年的水位變化進行多種情景分析及趨勢預測，為今后的博斯騰湖及其周圍水資源、生態環境和社會經濟的可持續發展提供科學依據。

1 研究區概況

博斯騰湖流域位于新疆塔里木盆地的西北和西北邊緣。由山區地表水資源形成區（包括開都河上游、黃水溝上游、清水河上游和20條時令河上游）、焉耆盆地綠洲耗水區（包括博斯騰湖）和孔雀河流域耗水區（包括庫爾勒市）三大單元組成，是一個典型的水—生態—經濟復合系統。山區地表水資源形成區人類活動影響較少，水文過程主要受氣候變化、地質地貌和土地覆蓋類型的影響。焉耆盆地綠洲耗水區和孔雀河流域耗水區受人類活動影響顯著。

圖1 博斯騰湖地理位置及水系

博斯騰湖是整個流域的樞紐，是焉耆盆地的最低凹處，氣候干旱，降水稀少，蒸發強烈，曾是我國最大的內陸淡水湖，既是開都河的尾閭，又是孔雀河的源頭，是新疆巴音郭勒蒙古自治州各族人民生活、生產的水源地。博斯騰湖水域遼闊，東西長達55km，南北平均寬約20km。在最高水位1 048.75m時，水面面積1 002.4km2，容積為88億 m3，在水位1 048m 時，平均水深為7.8m，最深為17m［23－25］。1958—2013年期間，水位變化幅度為1987年的1 044.95m到2002年的1 049.39m，非常洪水位和極限死水位之間相差近4m，變化幅度較大。博斯騰湖主要依賴入湖河流補給，入湖河流有開都河、黃水溝、清水河等，常年性入湖河流只有開都河，孔雀河是博斯騰湖唯一的出湖河流，出湖水量自1982年后由人為控制，即揚水泵站投入使用。

2 數據與方法

2.1 數據分析

本文利用的資料主要有大湖水位資料為大河口水文站1983—2013年逐年觀測數據；河水流量資料為1983—2013年的開都河出山口徑流量（大山口水文站）逐年觀測資料；1983—2013年的開都河東、西支入湖水流量（寶浪蘇木分水閘）逐年觀測資料；1983—2013年的揚水泵站逐年水量資料。蒸發與降水數據：大河口水文站1983—2013年逐年降水觀測值；博湖縣氣象站1983—2013年逐年觀測資料；社會經濟年系列資料：主要包括總人口、牲畜頭數、盆地灌區灌溉面積、灌溉定額、工業產值等等；數據來源1983—2013年新疆統計年簽與相關文獻。

2.2 研究方法

2.2.1 模型構建 本文以博斯騰湖水量平衡為核心，全面考慮博斯騰湖上游和下游的社會經濟發展和可能的氣候變化因素，結合系統動力學原理和應用要求，構建以博斯騰湖上游社會經濟發展需水量，博斯騰湖水量平衡和博斯騰湖下游（孔雀河流域）社會經濟發展需水量三部分組成的博斯騰湖水量平衡系統動力學模型，并繪制系統流圖（圖2）。博斯騰湖水量平衡中入湖水量大致分為三個部分：降水、地表水和地下水補給；耗水量包括博斯騰湖水面蒸發損失、蒸騰和博斯騰湖揚水等。系統參數主要由查閱相關資料、參考他人的研究成果和結合研究區實際情況而確定的。在相關的反饋回路里水位—庫容—面積三個因素之間的非線性關系就是通過Lookup Function表達（圖3）。模型模擬時間設置為1983—2013年，共31a，模擬步長為1a，以1982年的統計數據為初始值。

博斯騰湖水量平衡仿真模型中，所選擇的變量共有70多個，各個因素（變量）之間相互作用形成的因果反饋環路構成了模型的基本結構。本文給出了幾個最基本的要素（變量）之間的因果關系回路圖（見圖2）。可以看出，圖中有7個主要反饋環，其中正反饋環有5個：

→→I：博斯騰湖上游：增長率＋農田灌溉面積＋農田灌溉面積增長量。

→→Ⅱ：博斯騰湖上游：人口增長率＋人口數量＋人口增長量。

→→→Ⅲ：博斯騰湖：大湖水面積＋湖面降水量＋大湖地表入流量＋博斯騰湖（大湖）。

→→Ⅳ：博斯騰湖湖下游：農田灌溉面積增長率＋農田灌溉面積＋農田灌溉面積增長量。

→→V：博斯騰湖下游：流域人口增長率＋流域人口數量＋流域人口增長量。

負反饋環有2個：

→→Ⅵ：博斯騰湖：大湖水面積＋博斯騰湖（大湖）－蒸發量。

→→Ⅶ：博斯騰湖：博斯騰湖（大湖）＋模擬水位＋出大湖總水量－博斯騰湖（大湖）

圖2 博斯騰湖水量動態平衡的系統動力學模型

2.2.2 有效性檢驗 根據博斯騰湖1983—2013年歷年水量平衡數據運行模型，進行有效性檢驗，所得模擬數據與歷史監測數據進行比較，驗證其吻合度。由于模型變量較多，筆者主要對博斯騰湖水位、博湖上游的寶浪蘇木分水樞紐（第三水分樞紐）站徑流量和博湖下游的塔什店徑流量數據進行歷史驗證。模型通過相對誤差法（公式1所示）進行驗證。相對誤差通過模擬所造成的絕對誤差與真值之比乘以100%所得的數值，以百分數表示。相對誤差公式為：

式中：δ——實際相對誤差；Δ——絕對誤差；L——真值（其中：絕對誤差＝模擬測得值－真值）。寶浪蘇木分水樞紐（第三水分樞紐）站徑流量、博斯騰湖水位和塔什店徑流量的相對誤差分別3.9%，1.6%，1.2%，實測數據與模型數據基本吻合，最大誤差不超過5%，得出歷史數據與模擬數據基本吻合（圖4），表明模型具有較好的強壯性，能夠反映博斯騰湖水量平衡系統的實際特征。

圖3 博斯騰大湖水位－面積－庫容曲線

圖4 實測數據與模擬數據對比

3 結果與分析

3.1 水量平衡分析

以1983—2013年歷年水文、氣象、社會經濟等資料為基礎，按模型要求將數據存入模型的輸入文件中，應用上述模型進行模擬，得到相應的模擬結果（圖5）。

從1983—2013年的博斯騰湖入（出）湖水量、湖面蒸發耗水量、開都河生態需水量和孔雀河生態需水量的變化曲線（圖5）看出，博斯騰湖水量變化的影響因素是多方面的，其中入湖河水、出湖水量、湖面蒸發量等是影響博斯騰湖湖水位的主要自然因素，開都河和孔雀河生態需水量是影響博斯騰湖水位的主要人為因素。進入博斯騰湖的入湖水量直接受焉耆盆地自然來水的影響，隨自然來水量而變化。1983—2013年，從大湖出的總水量為14.4億m3，多年平均入湖總水量為23.9億m3，湖區總耗水量為9.5億m3，開都河東支入湖河水流量16.1億m3，西支入湖河水流量7.2億m3。西泵站建成后，出湖水量由人為控制，即揚水泵站投入使用，解放一渠直接向孔雀河輸水量也相繼減少，則孔雀河是博斯騰湖唯一的出湖河流，其中大湖區出流占到了50%，多年平均約9億m3，小湖區出流約6.4億m3。

圖5 博斯騰湖入（出）水（蒸發）量及生態需水量變化曲線

出湖水蒸發量受水域面積（水位）影響。1983—2013年蒸發量波動范圍為10.3～20.4億m3，而大湖水位變化范圍為1 045～1 048.62m。SD模型模擬計算表明，在過去的31a間（1983—2013年）博斯騰湖年均蒸發量為12.8億m3，增長速率為0.416億m3/a。其中，1983—1994年博斯騰湖大湖區水面年均蒸發量為12.02億m3，1995—2013年均蒸發量達到13.43億m3，增加了11.73%。而且，隨著湖區蒸發耗水量的增加，博斯騰湖水位還有下降的趨勢。

通過模擬計算得到開都河流域1983年、2000年和2013年生態需水量分別為14.16億m3，12.05億m3，17.41億 m3；孔 雀 河 流 域 1983 年、2000 年 和2013年生態需水量分別為11.87億 m3，26.19億m3，13.26億 m3。1983—2013年開都河、孔雀河流域生態需水量變化發現，近31a來，該流域生態需水量呈增加趨勢，特別是2000年之后增加迅速，特別是孔雀河流域需水量增加更明顯，在此期間，水位處于下降趨勢。

3.2 參數敏感度分析

敏感度分析（Sensitivity Analysis），又稱敏感性分析。參數敏感度分析的目的正是量化各系統參數對系統輸出的影響程度，從而對系統參數進行篩選，所得影響因子亦用于后續參數估計。表達式為［26］：

式中：S——模型參數敏感度；Y——變化后的狀態變量；Y0——基準條件下狀態變量值；ΔY——由參數變化引起的狀態變量的變化幅度；P——變化后的模型參數；P0——基準條件下模型參數值；ΔP——參數的變化幅度。具體含義見表1。

表1 敏感度取值

通過公式（2）對主要參數進行敏感度分析。敏感度較明顯的主要參數及其敏感度值和敏感度級別在表1所示。結果表明（表2），入湖徑流的敏感度值為最高，其次是湖面蒸發量、出湖水量和湖面降水量。其中，湖面蒸發量作為湖泊水量的支出項，與其水位呈負相關關系。表中蒸發量的敏感度值為負，是因為模型中蒸發量與博斯騰湖水位屬于負反饋回路，蒸發量增大水位減少。敏感度值排序為：入湖地表水量＞蒸發量＞出湖水量＞降水量＞灌溉定額＞灌溉面積＞工業產值＞人口。

表2 8個重要參數的敏感度值

3.3 情景分析

情景分析法是考慮影響未來發展諸多因素，系統分析相關問題，構建出多種可能的未來發展態勢，再對系統發展態勢作出自始至終的情景與畫面的描述［27－28］。其基本觀點是雖未來充滿不確定性，但有部分內容可以預測，這是由不確定性決定的。博斯騰湖水位過高或過低都會給生態環境和農業生產帶來影響。因而，在供水管理中最重要的是確定允許的水位范圍，允許水位確定之后，才能對湖泊的可利用水量進行合理分配，如灌溉、城市用水等。如低于此值允許范圍生態系統就會遭到破壞，因此，維護湖泊的合理水位已經成為淡水資源科學配置和永續利用的基本保證。

入湖流量作為湖泊水量的主要輸入項，對湖泊水位的變化有著顯著性的影響。開都河1958—1976年屬平水年，博斯騰湖水位呈波狀緩慢下降趨勢；1976—1985年開都河多年平均年徑流量達到歷史最低點，博斯騰湖水位也達到歷史最低點；1986—1995年多年平均徑流量開始回升，博斯騰湖水位呈急速上升趨勢；1996—2002年開都河多年平均徑流量達到歷史的最高點，博斯騰湖水位也達到歷史最高點。2003—2013年開都河多年平均徑流量開始減少，博斯騰湖水位呈急劇下降趨勢。從以上分析可以得出，入湖流量和水升降的變化規律非常接近。

氣候變化加劇了開都河年徑流變化的不確定性。從短期來看，隨著氣候變暖，冰雪融水增加，近期內河流地表徑流可能呈增加或減少波動狀態，甚至表現為增加趨勢，然而，長時間而言，隨著氣候變化天山山脈雪線上移，山區冰川儲量減少，來水量減少，開都河流域水文情勢可能會發生很大變化，年徑流量在氣候變化的影響下可能會出現重大波動，導致博斯騰湖水位變幅增加。為了避免大規模灌溉引水、跨流域調水以及土地利用變化等人類活動對徑流的影響，本研究選擇大山口水文站徑流數據為研究對象，可以較為單純地反映自然條件下博斯騰湖水位的變化。

農業用水包括灌溉用水和林木漁業用水兩部分，其中，灌溉用水占的比例非常大。農業是高耗水行業，灌溉用水對水資源的需求遠遠超過了工業用水、生活用水和牲畜用水，位居第一位。博斯騰湖流域灌溉用水占農業用水的比例非常大，因此，在研究中主要考慮灌溉用水。2013年，研究區用水總量為17.42億m3，其中灌溉用水量13.93億m3，占總用水量的80%，所以提高農業用水效率，減少灌溉用水量，把農業部門節約出來的水資源轉移到工業用水、生活用水和牲畜用水，從而可以緩解水資源短缺，這在干旱區顯得尤為重要。灌溉用水量取決于灌溉面積、灌溉定額和灌溉面積增長率。因此，根據博斯騰湖水量平衡最敏感影響因子和博斯騰湖流域水資源供需平衡與優化配置［29］，筆者設置兩種情景方案，預測分析博斯騰湖未來30a的水位變化情景。

A方案設置如下：情景1是用2013年的數據，各項參數保持不變進行模擬分析；情景2是除開都河多年平均徑流量34.3億m3（屬平水年）外，各項參數保持不變進行模擬分析；情景3是除開都河多年平均徑流量30.53億m3（屬枯水年）外，其他數據和參數與情景2相同；情景4是除開都河多年平均徑流量45億m3（屬豐水年）外，其他數據和參數與情景2相同。

B方案設置如下：情景1是用2013年的數據，各項參數保持不變進行模擬分析；情景2是開都河年平均徑流量保持不變、灌溉定額、灌溉面積增長率等人為因素數據和參數，進行模擬分析；情景3是除改變人為因素數據和參數外，開都河年平均徑流量與情景2相同；情景4是除改變人為因素數據和參數外，開都河年平均徑流量與情景2相同。B方案的情景指標值見表3。

表3 2014－2044年B方案的情景指標值

盆地綠洲區總耗水量的減少則可增加博斯騰湖總入水量。情景2，3，4都是在提高農業用水效率，減少灌溉用水量的條件下盆地綠洲區實施相同的情景方案，因此它們的農業耗水量也相同，年均約為5.84億m3，比情景1減少1.87億m3。方案A通過開都河平水年、枯水年和豐水年情景，模擬得到博斯騰湖2014—2044年的未來30a的水位變化趨勢。博斯騰湖多年平均水量平衡指標匯總結果見表4。

表4 博斯騰湖水量平衡主要指標匯總

方案A—情景1：開都河多年平均徑流量約為36.3億m3，入湖水量為31.13億m3，其中，東支入大湖水量年均為19.86億m3，西支入小湖水量年均為11.27億m3，出湖水量為12.63億m3，博斯騰湖水位波動范圍為1 043.46～1 047m，多年平均水位是1 045.68m，未來30a博斯騰湖水位持續下降。

情景2：開都河屬平水年，多年平均徑流量約為34.3億m3，入湖水量為28.81億m3，其中，東支入大湖水量年均為18.4億m3，西支入小湖水量年均為10.41億 m3，出湖水量為12.14億 m3，有所減少，但減少幅度不大，博斯騰湖水位波動范圍為1 042.97～1 047m，多年平均水位是1 044.69m，水位比情景1下降1.01m，未來30a博斯騰湖水位還是下降狀態。

情景3：開都河屬枯水年，多年平均徑流量約為30.53億 m3，入湖水量為25.24億 m3，其中，東支入大湖水量年均16.15億m3，西支入小湖水量年均9.09億 m3，出湖水量10.36億 m3，有所減少，減少幅度比情景1，2大，博斯騰湖水位波動范圍為1 041.01～1 047m，多年平均水位是1 043.27m，水位比情景1下降2.41m，平均每年下降0.07m，未來30a博斯騰湖水位還是下降狀態。

情景4：開都河屬豐水年，多年平均徑流量約45億m3，入湖水量40.16億m3，其中，東支入大湖水量年均25.56億m3，西支入小湖水量年均14.6億m3，出湖水量為17.12億m3，有所增加，增加幅度大，博斯騰湖水位波動范圍為1 047～1 048.13m，多年平均水位是1 048.41m，水位比情景1增加2.73m，平均每年增加0.089m，未來30a博斯騰湖水位處于高水位。全球氣候發生重大變化的情況下，開都河也很難一直保持豐水年狀態。因此，情景2適合博斯騰湖的未來實際情況，多年平均水位處于最佳運行水位。

方案B—情景1：開都河多年平均徑流量約36.3億m3，多年平均農業引取水量均值7.71億m3，生態用水量均值11.18億m3，農業引水量較低，博斯騰湖大湖區水位介于1 043.45～1 047m。

情景2：開都河多年平均徑流量約36.3億m3，多年平均農業引取水量均值12.7億m3，生態用水量均值15.82億m3，農業引水量高，博斯騰湖大湖區水位介于1 039.83～1 047m。

情景3：開都河多年平均徑流量約36.3億m3，多年平均農業引取水量均值11.89億m3，生態用水量均值14.60億m3，農業引水量較高，博斯騰湖大湖區水位介于1 041～1 047m。

情景4：開都河多年平均徑流量約36.3億m3，多年平均農業引取水量均值7.15億m3，生態用水量均值9.5億m3，農業引水量低，博斯騰湖大湖區水位介于1 045～1 047m。農業引水量越小，來水量越多，博斯騰湖水位會明顯提升。社會經濟發展需水日益加之的今天，在短時間內大幅度減少農業用水量是不客觀的。因此，情景4適合博斯騰湖的未來實際情況，水位處于最佳運行水位。

從上面的情景分析可以看出，方案A的情景2和方案B的情景4為最適合博斯騰湖及其上下游生態經濟可持續發展的最佳選擇。為了實現如上情景，應盡量保持現在的農業灌溉面積，有效控制農業灌溉面積的再擴大，通過高科技和產業結構的調整減少農業灌溉定額。

4 結 論

博斯騰湖上下游區域人類活動愈發頻繁，加之環境變化下博斯騰湖及其周圍區水—經濟—生態復合系統的復雜性和不確定性，使得水資源變成博斯騰湖流域社會經濟發展的瓶頸。因此，揭示該系統要素之間的復雜反饋關系及其影響機制變成社會經濟和生態環境協調發展的關鍵。本文證明系統動力學方法是研究復雜水—經濟—生態復合系統研究的最佳選擇。結果表明：

（1）環境變化下保持博斯騰湖及其周圍環境良好是一個具有挑戰性的任務。博斯騰湖水位的波動是氣候變化和人類活動耦合作用的結果。其中，近十年來湖水位的下降主要是受博斯騰湖上下游人類社會經濟活動影響的緣故，孔雀河流域社會經濟發展需水量猛增引起的影響更為突出。

（2）從敏感度分析可知，博斯騰湖水量平衡主要是入湖水量、蒸發量、出湖水量與降水量波動等自然因素和灌溉定額、灌溉面積、工業產值與人口等人為因素共同作用的結果，且自然因素的綜合影響程度高于人為因素。

（3）方案A在保持其他參數不變情況下，調節開都河年平均徑流量，方案B在開都河年平均徑流量保持不變的情況下，調節灌溉定額、灌溉面積增長率等人為因素數據和參數，由這兩個方案來產生不同的水位狀況并預測未來30a的湖泊水位變化。結合方案A與方案B分析可知，單一增加入湖徑流量等自然因素或者減少灌溉面積、灌溉定額等人為因素并不能使博斯騰湖流域社會經濟實現可持續發展。例如，加大博斯騰湖上游區域節水力度可以有效保證博斯騰湖入湖徑流量。此外，在干旱區保持較高或較低湖泊水位對上下游區域社會經濟發展造成影響。

綜合考慮博斯騰湖不同水位變化情景，對于氣候變化與人類活動影響均較為突出的博斯騰湖流域，在未來氣候變化這個大環境下，為了維持博斯騰湖水位在生態水位之內，同時實現周圍區社會經濟的可持續發展，必須用嚴格控制農業灌溉面積，加大農業節水技術等方法來有效緩解博斯騰湖水位的劇降和水量失衡問題。調整社會經濟用水與生態用水比例以及各產業用水比例，推進節水型社會建設，提高水資源利用率。

［1］ 張紅亞，呂明輝.水文學概念［M］.北京：北京大學出版社，2007.

［2］ 高華中，姚亦鋒.近50a來人類活動對博斯騰湖水位影響的量化研究［J］.地理科學，2005，25（3）：305－309.

［3］ 馬麗娟，趙景峰，張宏俊，等.氣候變化背景下冰川積雪融水對博斯騰湖水位變化的影響［J］.干旱區地理，2010，33（2）：210－216.

［4］ 李新國，李會志，王影，等.開都河下游灌區土壤鹽漬化研究［J］.水土保持研究，2010，17（6）：111－114.

［5］ 秦伯強.氣候變化對內陸湖泊的影響分析［J］.地理科學，1993，13（3）：212－219.

［6］ 李佳秀，徐長春，王曉，等.開都河流域極端水文事件變化及其對氣候的響應［J］.水土保持研究，2014，21（5）：285－291.

［7］ 陳亞寧.嚴守新疆博斯騰湖生態紅線［EB/OL］.（2013－12－31）［2015－03－05］.中國科學報.

［8］ 李二輝，穆興民，趙廣舉.1919—2010年黃河上中游區徑流量變化分析［J］.水科學進展，2014，25（2）：155－163.

［9］ 汪鷺.博斯騰湖流域主要入湖河流徑流序列分析［D］.成都：四川師范大學，2008.

［10］ 王維霞，王秀君，姜逢清，等.開都河流域上下游過去50a氣溫降水變化特征分析［J］.干旱區地理，2012，35（5）：747－753.

［11］ 巴州國土規劃博湖規劃組.新疆巴音郭楞蒙古自治州國土規劃博斯騰湖保護治理和資源開發規劃報告［R］.1984：4－9.

［12］ Rusuli Y，Li L，Ahmad S，et al.Dynamics model to simulate water and salt balance of Bosten Lake in Xinjiang，China［J］.Environmental Earth Sciences，2013，74（3）：1－12.

［13］ 裴新國.博斯騰湖水鹽動態及其平衡的研究［J］.干旱區地理，1988，1（13）：1－7.

［14］ 裴新國.排入博斯騰湖的農田廢水水量、鹽量的估算［J］.干旱環境監測，1989，3（2）：43－45.

［15］ 成正才，李宇安.博斯騰湖的水鹽平衡與礦化度［J］.干旱區地理，1997，20（3）：43－49.

［16］ 王亞俊，李宇安，王彥國，等.20世紀50年代以來博斯騰湖水鹽變化及趨勢［J］.干旱區研究，2005，22（3）：355－360.

［17］ 孫占東，王潤.氣候背景下冰川在博斯騰湖水量平衡中的作用［J］.湖泊科學，2006，18（5）：484－489.

［18］ 鐘瑞森，董新光，王新菊.博斯騰湖水鹽變化及可調出水量預測［J］.人民黃河，2007，29（11）：45－47.

［19］ 鐘瑞森，董新光.新疆博斯騰湖水鹽平衡及水環境預測［J］.湖泊科學，2008，20（1）：58－64.

［20］ Rusuli Y，Li L，Ahmad S，et al.Dynamics model to simulate water and salt balance of Bosten Lake in Xinjiang，China［J］.Environmental Earth Sciences，2013，74（3）：1－12.

［21］ 劉麗梅，趙景峰，張建平，等.近50a博斯騰湖逐年水量收支估算與水平衡分析［J］.干旱區地理，2013，33（1）：33－40.

［22］ 伊麗努爾·阿力甫江，海米提·依米提，麥麥提·吐爾遜，等.1958—2012年博斯騰湖水位變化驅動力［J］.中國沙漠，2015，35（1）：240－247.

［23］ 李衛紅，袁磊.新疆博斯騰湖水鹽變化及其影響因素探討［J］.湖泊科學，2002，14（3）：223－227.

［24］ 夏軍，左其亭，邵民成.博斯騰湖水資源可持續利用：理論、方法、實踐［M］.北京：科學出版社，2003.

［25］ 孫占東，王潤，黃群.近20a博斯騰湖與岱海水位變化比較分析［J］.干旱區資源與環境，2006，20（5）：56－60.

［26］ 謝媛媛.SWAT模型在黃土丘陵區參數敏感度分析及率正研究［J］.水土保持研究，2012，19（4）：204－206.

［27］ 劉瀏，胡昌偉，徐宗學，等.情景分析技術在未來太湖水位預見中的應用［J］.水利學報，2012，43（4）.

［28］ 吳威，吳松，陳爽.基于SD模型分析的環鄱陽湖地區發展模式探討［J］.湖泊科學，2012，24（2）：252－258.

［29］ 陳亞寧，杜強，陳躍濱，等.博斯騰湖流域水資源可持續利用研究［M］.北京：科學出版社，2012.