自然-社會水循環模型估算平原-丘陵-濕地區水稻種植潛力

宮興龍，付 強，孫愛華，關英紅，王 斌，李 茉

?

自然-社會水循環模型估算平原-丘陵-濕地區水稻種植潛力

宮興龍，付 強※，孫愛華，關英紅，王 斌，李 茉

（東北農業大學水利與土木工程學院，哈爾濱 150030）

自20世紀中葉以來在平原-丘陵-濕地區隨著井灌水稻熱的興起，水稻種植面積急劇擴大甚至超過此類地區水稻種植潛力，引起了濕地面積萎縮、地下水水位下降和水質惡化等問題。為解決這些問題，該文構建自然-社會二元水循環模型揭示人類影響下的平原-丘陵-濕地區水循環機理，進而推算適宜的水稻種植面積。首先，利用數值法構建了基于柵格的分布式水文模型，然后在此模型中嵌入人類活動影響模型從而構建了二元水循環模型。根據水稻種植潛力閾值抽取河道水量極值、地下水最大埋深、地下水最大降深和最大地下水開采量等因素與旱田作物種類組合了18種情景模式。在滿足控制閾值條件下，依據多年平均日降雨、蒸發和情景模式利用二元模型計算了各種情景的水田種植潛力，由潛力分析得：18種模式地下水最大埋深在5.17~7.49 m之間，地下水最大降深在1.67~3.73 m之間；水田處最大坡度范圍0.028~0.053；河道引水量占河道水量的50%~70%；地下水開采比例在79%~112%間；水稻種植潛力為28.36萬hm2~54.12萬hm2，來自旱地面積為21.05萬hm2~40.32萬hm2，來自未利用地面積為5.68萬hm2~11.09萬hm2。以情景模式17為例在水田生育期內對河道生態需水量、地下水水位和旱地作物蒸發等進行了檢驗，驗證得到：整個水稻生育期內分區流域地下水埋深均小于7.12 m，開發水田的分區流域基流比最小值為33.45%，分區流域旱田平均蒸發與1990年土地類型情景下的分區流域旱田平均蒸發的比值大部分位于0.98~1.05，說明水稻生育期內情景17對河道水量和地下水水位的影響在控制范圍內，對旱田蒸發影響比較小，因此情景模式17的水稻種植潛力是可行的。研究可為描述平原-丘陵-濕地區的水文循環過程和推求水田開發潛力提供依據。

水循環；蒸發；社會；平原-丘陵-濕地區；水稻種植潛力；情景模式；格子波爾茲曼法

0 引 言

20世紀中葉以來三江平原隨著水稻熱的興起，特別是井灌區的發展，抽水量越來越大，引起平原-丘陵-濕地區濕地面積逐漸萎縮，地下水水位下降水質惡化，加重了一些區域水分虧缺態勢和發展不可持續等問題[1-4]。這些問題產生的主要原因是抽取河道和地下水使水循環逐漸偏離抽水前的水文循環，從而使其機理發生了變化[5-6]。為了揭示這種變化和防止地下水水位下降水質惡化等問題的發生，必須給出平原-丘陵-濕地區的水田種植潛力，控制井灌水田面積。針對水田種植潛力或地下水開發潛力學者進行了定量的研究。如王韶華等[7]根據水量平衡原理應用農業需水量和地下水可開采模數2種方法計算了研究區水稻適宜種植面積。趙清[8]根據研究區農業需水量和降水情況計算了水稻適宜種植面積。付強等[9]在灌溉定額已知條件下用水分運動數學模型計算了三江地區水稻發展潛力。雖已進行了富有成效的研究，但存在以下4個問題：1）在計算種植潛力時不論研究區域多大，總是將其作為一個整體研究，對降雨、蒸發和地下水等影響種植潛力因素的空間分布考慮不足；2）在計算種植潛力時大部分研究是利用作物整個生長期水量供需平衡來計算潛力，對于生長期內的各個時期影響因素考慮不足或者是沒有考慮；3）計算種植潛力時采用的模式多數是集總式或基于集總式的分散水文模型，難于描述土地結構時空分布對水文過程量時空分布的影響；4）對于種植的作物種類的不同和耕作制度等人類活動因素的影響考慮不足。

目前反映人類影響水文循環的主要工具之一是自然-社會二元水循環模型，即在自然水循環模型里加入人類活動影響的水文循環模型。針對這一模型學者進行了大量的研究工作，取得了豐碩的研究成果。如王浩等[10-11]依據自然-社會二元水循環模式將下墊面變化和人工取用水作為模型變量，構建水資源評價新方法。謝新民等[12]從二元水循環角度研究了人類影響下的華北平原區地表水與地下水統一評價模型。吳普特等[13]進一步發展了這種水模式提出了實體-虛擬水“二維三元”耦合流動理論。張俊娥等[14]利用分布式水文模型MODCYCLE描述天津市二元水循環，并對天津市水量平衡狀態進了分析。從這些實例可以看出，為了揭示人類影響下的多水源區水循環機理，進而推算出適宜的水稻種植面積，必須建立人類活動影響下的降水、河流水、濕地水、水庫水、灌溉水和地下水等水流運動的數值模型即自然-社會二元水循環模型。

鑒于此，為了考慮種植潛力影響因素的時空分布對種植潛力的影響，本文構建了平原-丘陵-濕地區基于柵格的分布式水文模型。為了反映人類的影響，在旱田和水田處按照管理制度建立人工控制模型，將此控制模型嵌入分布式水文模型中，構成自然-社會二元水循環模型。在滿足研究區域的地下水開發控制水深和流域生態需水量等條件下，利用構建的流域二元水循環模型計算出研究區的水稻種植潛力。

1 自然-社會二元水循環模型構建

二元化流域水循環的重要特征有4個[10-11]：二元化的水循環服務功能、水循環結構和二元化的參數、二元化的水循環路徑、二元化的水循環驅動力，其示意結構見圖1。二元結構模型的本質是將自然與人類影響二元結合起來描述水循環，因此其本質是二元化。二元模型運行的核心是循環解耦股和參數的二元化。描述二元作用下的水流循環路徑或通道是其表征。二元模型運行的驅動力是二元化的基礎（圖1）。從圖1可以看出，二元水資源模型的核心是自然和人類共同影響的循環結構，從建模角度即自然水循環過程和人工水循環過程的耦合結構。

圖1 二元水循環模型結構

1.1 自然水循環模型

研究三江平原地區自然水循環時，必須考慮到地下水運動。本文構建的自然狀態下水文循環模型架構如圖2所示。從自然水循環角度和構建模型需要，參照TOPMODEL描述產流的思路，將垂向水流過程分為降雨、根系區蒸散發、非飽和帶滲流和飽和帶基流，描述水平向匯流時將匯流過程分為地表水坡面匯流、地下水匯流和河道洪水波運動[15]。TOPMODEL模型在描述匯流時采用的方法是時間-面積曲線，時間-面積曲線只能按距離出口等流時的條帶來描述匯流過程，難于準確描述匯流的時空過程。基于此，在構建匯流過程時，本文采用基于格子波爾茲曼法（lattice Boltzmann，LBM）改進了TOPMODEL的匯流和非飽和帶土壤水運移模塊，在柵格上對各模塊進行松散耦合構建基于LBM的分布式TOPMODEL[16-19]，通過調整柵格大小來改變產匯流過程的時空尺度。自然水循環模型的模塊為：面雨量空間分布是根據雨量站實測的降雨采用距離倒數推求；根系區的蒸散發模塊采用TOPMODEL設計時采用的模塊見式（1）；非飽和帶土壤水分運動采用基于LBM求解的Richards方程而構建的水分運動模塊；坡面匯流模塊是利用LBM求解坡面和河道運動方程而構建的模塊；河流/渠道模塊是利用LBM五速模型求解線性擴散波方程而構建的模塊；飽和帶水平向土壤水運動和潛水運動采用達西定律建立的數值模塊。

圖2 自然水循環模型示意圖

垂向水流運動模塊主要是解決根系區、非飽和帶、飽和帶和潛水之間的水流運動。將網格降雨量和鄰近網格地表水流入量作為根系區柵格輸入項。輸出項為根據根系區含水率和水面蒸發量計算的根系區蒸發量。根系區含水率大于根系區容水量時，根系區含水率與根系區容水量的差值量作為下滲水量。根系區下滲水量作為非飽和帶上邊界輸入量，由非飽和帶運移數學模型計算出非飽和帶下滲入飽和帶的水量。飽和帶含水率達到田間持水量時向潛水補給水量。

水平向水流流動是按地表、飽和地下水基流和潛水水流3層來模擬的。地表水匯流時將流域上網格分為平原-丘陵區和濕地區。平原-丘陵區邊界為流域分水線和與濕地交接處，這一區域采用基于LBM的匯流模型。在垂向產流模塊模擬的網格水深基礎上，利用匯流模型模擬得網格的水深。

1）蒸散發模塊

流域內非農田柵格的蒸散發

式中為蒸散發，mm；0為水面蒸發量，mm；max為根系區最大容水量，mm；S為根系區缺水量，mm。

2）非飽和帶土壤水運移數學模型

非飽和帶土壤水運動采用Richards方程描述

式中D為擴散系數，m2/s；為導水系數，m/s；為含水率，cm3/cm3；為距離地面深度，m；為時間，s。

式（2）采用三速模型可得平衡態分布函數

3）飽和帶基流數值模塊

飽和帶土壤水基流

4）坡面匯流數值模型

坡面流采用運動波方程描述

用九速度模型求得式（5）的平衡態分布函數為

坡面流水文要素可由式（7）描述：

式中e、e分別為和方向的粒子運動速度，m/s。

5）河道匯流數學模型

河道水流用線性擴散波方程描述

式中c為河道水流波速，m/s；D為河道水流擴散系數，m2/s；為流量，m3/s。

選用五速模型可推得式（8）的平衡態分布函數

式中A、B、C、D、E為中間變量。

將離散速度和平衡分布函數代入式（7），可得河道水流運動要素。

1.2 人類活動影響模型構建

研究區設定在三江平原下游的撓力河流域。該區域影響水文循環的人類活動比較多，如縣、農場場部生活用水、綠化用水、養殖業、井灌區灌溉用水等。由于本地區生活用水、綠化用水和養殖業用水等用水與農業用水相比非常小，因此著重研究農業用水對水循環的影響。本區農業用水主要包括旱地農作物需水和稻田用水。

1）農作物生態需水量模型

需水量ET采用如下公式

ET=K·ET0（10）

式中K為作物系數；ET0為潛在蒸散量，mm。

潛在蒸散量采用聯合國糧食與農業組織(FAO)推薦的Penman-Monteith（1998）模型[20]

式中Δ為飽和水汽壓-溫度曲線斜率，kPa/℃；R為凈輻射，MJ/m2；為土壤熱通量，MJ/m2；為干濕計常數；1為空氣平均溫度，℃；2為地面上2 m處風速，m/s；e為飽和水汽壓，kPa；e為實際水汽壓，kPa。

采用聯合國糧食與農業組織推薦的方法[21]和文獻[22-25]提供的方法進行修正作物系數，得撓力河流域作物系數見表1。

表1 撓力河流域作物系數

2）河道內生態用水量

河道內生態用水量采用文獻[26]提出的基流比與河道健康狀況對應關系來描述，見表2。

表2 河道生態流量健康狀況基流比

3）水稻生育期適宜水深及蒸發

統計前人的研究成果[7,27-28]，并根據實地調查、向當地水稻管理部門咨詢和理論分析得到建三江地區的水稻生育期適宜水深和蒸發量（見表3）。

表3 建三江地區水稻生育期及適宜水深

注：0為水面蒸發量。

Note:0is water surface evaporation amount.

4）地下水控制水位

地下水水位觀測井主要集中在八五三、八五二、和五九七等農場，資料年限為1997—2008年。由水位數據分析得撓力河流域4月末—5月初地下水埋深為1.15~ 4.56 m。根據水文地質條件，開采水溫條件、社會生產成本情況、防治鹽堿、提高防洪除澇抗旱庫容和相關學者研究成果[29]，提出“理想埋深（水位）”概念。由地下水水位數據和參考文獻[29]得三江平原下游建三江地區理想埋深應控制在10～17 m范圍內。采用10 m作為控制埋深，由初始埋深最大為4.56 m，得開采地下水時降深應控制在5.44 m以內。

2 研究區概況及模型應用

本文實例為三江平原下游濕地、平原和丘陵交融的區域-撓力河流域，其基本概況如圖3所示。

圖3 研究區基本概況

圖3a為流域數字高程模型（digital elevation model，DEM)。圖3b~圖3d為利用DEM提取的流域水系、流域分區和流域坡度分布。圖3e和圖3f分別為1990年和2013年的土地利用情況。由圖3e和圖3f可得1990—2013年期間撓力河流域耕地（旱地、水田等）面積由121.21萬hm2增加到137.63萬hm2，耕地面積變化不顯著，但是土地類型空間分布發生了顯著變化。變化的主要是水田和未利用地，水田面積由11.73萬hm2增加到41.12萬hm2，未利用地由16.56萬hm2減少到1.61萬hm2，大部分轉為旱地和水田。旱地面積由109.50萬hm2減少到96.59萬hm2，變化不明顯。

撓力河流域內有基巖裂隙孔隙含水層、第四系孔隙承壓含水層、第三系裂隙承壓含水層和第四系孔隙含水層共4個含水層，其中第四系孔隙含水層較厚供水量充沛，從山前向河谷含水層逐漸增厚，地下水埋深逐漸變淺，地下水易于補給和排泄。撓力河流域農業灌溉用水主要來自于第四系含水層。基于此，本次模擬以第四系含水層為主。該區的主要補給方式是河流滲漏、降水入滲和側向地下徑流；人工抽水井開采是地下水排泄的主要方式，在開采量小的時間里還有一部分河流排泄和側向徑流排泄。

3 結果與分析

3.1 模型驗證

降雨量、水面蒸發量和徑流量來源于水文年鑒，資料年限為1964—1989年。地下水資料取自1997—2008年長觀測井的觀測資料。土壤數據來源于寒區旱區科學數據中心(http://westdc.westgis.ac.cn)，其用于估算土壤飽和含水量、田間持水量等水分常數。依據1974—1981年資料率定后的參數，應用1982—1989年7場次洪進行驗證，參加驗證的7場次洪的徑流相對誤差的絕對值區間為1.78%~12.24%，平均值為4.97%；洪峰流量擬合相對誤差的絕對值區間為1.71%~11.55%，平均值為4.52%。洪水過程線的確定性系數區間為0.67~0.85，平均值為0.75。具體驗證過程見文獻[16-18]。

3.2 情景模式設置及水稻種植潛力分析

3.2.1 情景模式

影響水稻種植潛力的因素較多，包括自然供水能力、旱田作物種類和灌溉管理制度等因素，為了能全面反映這些因素和不同因素組合的影響，采用情景模式加以反映。利用二元模型模擬情景模式對環境的影響，優選出適用和可靠的情景模式，其水稻種植面積即為水稻種植潛力。

由1997—2008年觀測井水位記錄信息分析得撓力河區域多年平均4月末—5月初地下水埋深分布見圖4。圖是撓力河流域可耕地區域地下水水位信息和非耕地區域土地類型信息的疊加，非耕地區域土地類型只有4種，因而圖中土地類型只有4種。為了控制地下水超采和防止地下水水位下降過大設定3個指標。指標1采用地下水埋深，要求最大埋深不超過10 m（即最大降深不超過5.44 m）。指標2是水田灌溉末期地下水水位減去水田灌溉初期水位，要求不得超過0.24 m，這是根據灌區多年地下水水位數據統計得到。指標3是地下水開采量。文獻[7,30]在計算三江平原地下水開采能力時，認為本文研究區地下水開采比例不應超過0.8。文獻[31-32]在計算開采能力時，采用地下水開采比例為1時的開采量為極限開采量。這些文獻提出地下水的開采比例是指一年提水量與補給量之比，本文的開采比例是指水田生育期內提水量與補給量之比。水稻生育期開始到次年4月末期間的地下水補給量是生育期內補給量的1.1倍。因此本文地下水開采比例設定0.8、1.0和1.1共3種情景。

為了滿足河道生態用水，本文設置河道引水為河道來水量的50%~70%，此時河道剩余水量為30%~50%，因而河道健康狀態為一般到良好狀態范圍。由文獻[33]數據可得2003－2012年期間，黑龍江省玉米和水稻耕地面積年均增長分別為14.6%和12.9%，這期間大豆耕地面積急劇遞減。小麥和薯類耕地面積與大豆耕地面積比相對較小，根據情景模式將其歸為小麥、玉米和大豆的耕作面積內。撓力河流域旱田作物種植情況與黑龍江省情況相似，主要為小麥、玉米和大豆3種。為了反映農作物需水對水田耕地面積的影響和計算方便，設置流域旱地呈現小麥、玉米和大豆3種極限種植情況。由以上影響因素歸納為18種情景如表4所示。

圖4 灌溉前期（4月末—5月初）地下水埋深多年平均分布

表4 18種情景模式

3.2.2 水稻種植潛力

開墾水田的原則：1）不占用撓力河流域上已有的林地、建筑物用地、草地和水域用地等用地；2）開墾水田處坡度不大于國家規定的最大耕地坡度0.176（國家規定開發耕地最大坡度）；3）水田開墾是按流域坡度從小到大順序進行。當水稻生育期水田處水深達不到表3適宜水深時，水田取水原則：1）取水田處上游坡面來的水；2）取河道用水，各情景河道取水量占河道水量不得高于表4的引河道水量；3）抽取地下水，要求最大降深不得超過5.44 m。具體計算潛力過程如下：1）在圖3c的流域分區內按照開墾水田的原則設定水田分布；2）將多年平均日降雨、蒸發、旱田分布和水田分布等作為輸入數據，利用二元水循環模型模擬流域水循環過程；3）將滿足情景控制條件的最大水田面積作為水田種植潛力。

將18種情景模式作為控制因素，多年平均日降雨和蒸發作為輸入數據，利用二元水循環模型模擬研究區域水文循環，得潛力如表5所示。由表5可知：18種情景地下水最大埋深在5.17~7.49 m之間（<10 m）；地下水最大降深在1.67~3.73 m之間（<5.44 m）；開發水田處最大坡度范圍在0.028~0.053（<0.176）；地下水開采比例在79%~ 112%間，接近110%；時段末水位減去時段初水位（水田生育期末地下水水位減去水田生育期初期地下水水位）范圍<0.24 m；水稻種植潛力為28.36萬hm2~54.12萬hm2，來自旱地轉入水田的面積為21.05萬hm2~40.32萬hm2，來自未利用地轉入水田的面積為5.68萬hm2~11.09萬hm2。

表5 水稻田種植潛力分析

開采地下水比例≥100%的情景有可能產生超采情況，采用情景3、情景6、情景9、情景12和情景15等模式時必須關注這一問題。為了盡可能地利用地下水，不予優先采用地下水開采比例比較小（79%）的情景1、情景4、情景7、情景10、情景13和情景16模式。據文獻[33]可知，近些年來玉米種植面積增加迅速，因此本文推薦情景17模式的水田種植面積為撓力河流域水稻種植潛力如圖5所示。從圖5和圖3可知，流域南部完達山山脈，西南和東南地勢比較高，這些地區為坡度變化較大，被林地覆蓋的低山丘陵帶，不適宜耕地。撓力河流域中部與北部區域為地勢平緩、地下水資源豐富的地區，為土地開發及生產條件優越的地區。

圖5 情景17土地類型分布

3.2.3 河道取水量隨時間變化情況

為了描述水田生育期內河道取水對河道水量影響，以情景17的水田分布作為水稻種植潛力分布，利用1964—1989年間的平均降雨和水面蒸發作為輸入數據模擬得分區流域多年平均河道水量。由情景17模擬的分區流域河道水量除以分區流域多年平均河道水量得基流比如圖6所示。由圖6和表2可知，開發水田的分區流域基流比最小值為33.45%（>30%），即在情景17模式下的河道健康狀況最低程度為一般健康。由圖6可知，分區流域基流比下降比較快和基流比比較小的時間發生在適宜水深相對于前一個生育期增加的時間點后幾天，如水稻返青期和分蘗初開始的最初幾天。

圖6 基流比隨種植水稻天數變化曲線

3.2.4 地下水水位隨時間變化情況

采用情景17情況水田布局，利用構建的二元水循環模型模擬的分區流域地下水埋深隨時間變化過程線如圖7a所示。由圖7a可知：水稻生育期內地下水埋深均小于7.12 m；地下水埋深在水田耕作期的最初1~38 d內呈現下降趨勢，但不是十分明顯；在39~92 d期間地下水埋深下降比較明顯；在93~120 d期間地下水水位緩慢的恢復；121~151 d期間地下水恢復過程明顯加快。由表3可知，水田耕作期的頭38 d是泡田、插秧期和返青期，這一期間水田蒸發量比較小，而這一期間河道引水量和降雨量能滿足水田大部分用水，因此這段時間內地下水埋深下降趨勢不是十分明顯。而在39~92 d期間是水稻生長的主要階段，蒸發量非常大，分蘗期達到9.61 mm，這時期降雨和河道水量滿足不了蒸發量，大量開采地下水導致地下水迅速降低呈現圖7a所示的過程。93~120 d為抽穗期到黃熟期，這一時期水田用水量逐漸減小，地下水進入恢復期。120 d后水稻成熟后不需要水，由于濕地和降雨的補給，地下水水位進入快速恢復期。為了能反映種植水稻期間地下水最大埋深的空間分布情況，按分區流域繪制了地下水最大埋深分布圖如圖7b所示。由圖7b、圖3可知，地下水最大埋深發生在流域上游山丘區與平原區域交接處的分區流域，究其原因為這些區域地下水賦存條件不好和距離流域內較大的河流較遠，導致其利用地下水比例大于其他區域。

圖7 地下水埋深隨時間變化及其最大值空間分布

3.2.5 旱地作物蒸發隨時間變化情況

為了檢驗水田種植面積擴大之后對旱地作物的影響，在相同的降雨和模型參數情況下，利用二元模型模擬情景17土地類型情況下的分區流域旱田平均蒸發與1990年土地類型情景下的分區流域旱田平均蒸發的比值如圖8所示。從圖8可知，從種植水田起到距離種植水田起始日的120 d里，該比值的范圍位于0.95~1.16，多數值位于0.98~1.05。最大值出現于第32天的1.16，其次位于第69 天的1.15，比值>1的部分出現于種植水稻的第30天和第65天周圍。由表3可知，在第32~38天區間時水稻處于返青時期蒸發量比較小，這時的降雨足以滿足蒸發，水田多余的水量被排出，導致水田周圍的旱田蒸發量增大。種植水田的第62~63天是曬田時期，需要將水田水放出來，這導致水田周圍的旱田蒸發增大。

圖8 情景17分區流域旱田平均蒸發與 1990年分區流域旱田平均蒸發比值

綜上，情景17的水田處最大坡度0.051，小于國家規定開發耕地最大坡度0.176；耕地主要來自未利用地和旱地，未開發撓力河流域林業、草業和濕地等資源；情景17抽取地下水灌溉水田時地下水最大埋深為7.12 m，小于研究區控制埋深10 m；情景17水田灌溉末期地下水水位減灌溉開始時地下水水位值為0，說明在灌溉結束時地下水水位已恢復；情景17對河道水量影響顯著，基流比最小值為33.45%（>30%），處于一般健康狀態；情景17的旱地蒸發與1990年土地類型旱地蒸發之比位于0.95~1.16間，大部分值位于0.98~1.05之間。這些數據表明水稻生育期內情景17對河道水量和地下水水位的影響在控制范圍內，而對旱田蒸發影響比較小，因此可得情景17所得水田種植潛力是可行的。

4 討 論

土地利用變化和人類活動對水文循環產生影響，而反過來水文循環變化又影響人類活動和土地利用，他們之間存在復雜關系[18]。本文參照TOPMODEL模型[12]和二元水循環模型[10-12]建模模型的思路和方法，用數值法構建了基于柵格的分布式水文模型，在模型中嵌入人類活動影響模型從而構建了二元水循環模型。模式在考慮人類活動時，沒有考慮經濟社會的影響，僅考慮旱田和水田等的人類活動，對建三江地區居住者生活和綠化等用水沒有考慮。水田適宜水深、蒸散發量和作物系數等受耕作者習慣、氣候和種植區域地形地貌等因素綜合影響[34]，難以準確定量化，本文獲得的這些量可能與實際量存在一定程度的偏差。在使用作物系數時，由于難定量描述農作物播種前的灌溉需水量，使得計算數據與實際可能存在偏差。由于受研究區域地質和地下水資料限制，在分析地下水水位、埋深、降深和開采地下水比例時是按圖3c的流域分區進行的。如資料充足可按照更小的尺度計算種植潛力，這時的潛力數據與本文的計算數據應存在一定程度的偏差。

5 結 論

在撓力河流域建立了二元水循環模型，由洪峰流量擬合相對誤差的絕對值區間為1.71%~11.55%和洪水過程線的確定性系數區間為0.67~0.85可得，模型較好地反映研究區的水文循環過程。采用該模型研究人類影響的平原-丘陵-濕地區的水循環機理，主要結論如下：

1）18種模式地下水最大埋深在5.17~7.49 m之間，地下水最大降深在1.67~3.73 m之間；水田處最大坡度范圍0.028 ~0.053；河道引水量占河道水量的50%~70%；地下水開采比例在79%~112%間；時段末水位減去時段初水位下降范圍<0.24 m；水稻種植潛力為28.36萬hm2~ 54.12萬hm2，來自旱地面積為21.05萬hm2~40.32萬hm2，來自未利用地面積為5.68萬hm2~11.09萬hm2。

2）情景17模式下的分區流域基流比最小值為33.45%，河流健康狀況為一般健康。分區流域水稻生育期內地下水埋深均小于7.12 m，水田種植對旱田蒸發有一定影響，但影響比較小，情景17所得的水田種植潛力是可行的。

在水文氣象、水文地質、灌溉制度和地下水水位等資料不足情況下，構建的二元水循環模型按分區流域盡最大可能利用這些資料，為描述平原-丘陵-濕地區的水文循環過程和推求水田開發潛力提供了參考。

[1] 劉吉平，杜保佳，盛連喜，等. 三江平原沼澤濕地格局變化及影響因素分析[J]. 水科學進展，2017，28(1)：22－31. Liu Jiping, Du Baojia, Sheng Lianxi, et al. Dynamic patterns of change in marshes in the Sanjiang Plain and their influential factors[J]. Advances in Water Science, 2017, 28(1): 22－31. (in Chinese with English abstract)

[2] 周浩，雷國平，楊雪昕，等. RCPs氣候情景下三江平原典型流域耕地動態模擬[J]. 農業機械學報，2017，48(10)：121－133. Zhou Hao, Lei Guoping, Yang Xuexin, et al. Simulation of Cultivated Land under RCPs Scenarios in Typical Basin of Sanjiang Plain[J]. Transactions of the Chinese Society for Agricultural Machinery, 2017, 48(10): 121－133. (in Chinese with English abstract)

[3] 曲藝，羅春雨，張弘強，等. 基于歷史生物多樣性與濕地景觀結構的三江平原濕地恢復優先性研究[J]. 生態學報，2018，38(16)：5709－5716. Qu Yi，Luo Chunyu，Zhang Hongqiang, et al. Prioritization of wetland restoration in Sanjiang Plain, based on historical biodiversity and landscape structure[J]. Acta Ecologica Sinica, 2018, 38(16): 5709－5716. (in Chinese with English abstract)

[4] 危潤初，肖長來，方樟. 黑龍江建三江地區地下水動態趨勢突變點分析. 吉林大學學報：地球科學版，2016，46(1)：202－210.Ｗei Runchu, Xiao Changlai, Fang Zhang. Trends mutation nodes of groundwater dynamic in Jiansanjiang area of Heilongjiang province[J]. Journal of Jilin University: Earth Science Edition, 2016, 46(1): 202－210. (in Chinese with English abstract)

[5] 周浩，雷國平，張博，等. 1990－2013年撓力河流域耕地變化下水土資源平衡效應分析[J]. 農業工程學報，2015，31(1)：272－280. Zhou Hao, Lei Guoping, Zhang Bo, et al. Farmland change induced land and water resource balance in Naoli River Basin from 1990 to 2013[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(1): 272－280. (in Chinese with English abstract)

[6] 任守德，付強，王凱. 基于宏微觀尺度的三江平原區域農業水土資源承載力[J]. 農業工程學報，2011，27(2)：8－14. Ren Shoude, Fu Qiang, Wang Kai. Regional agricultural water and soil resources carrying capacity based on macro- micro scale in Sanjiang Plain[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 27(2): 8－14. (in Chinese with English abstract)

[7] 王韶華，劉文朝，劉群昌. 三江平原農業需水量及適宜水稻種植面積的研究[J]. 農業工程學報，2004，20(4)：50－53. Wang Shaohua, Liu Wenchao, Liu Qunchang. Agricultural water consumption and suitable paddy rice plant areas of the Three-River-Plain[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2004, 20(4): 50－53. (in Chinese with English abstract)

[8] 趙清. 三江平原建三江地區農業需水量及水稻生產潛力研究[J]. 現代農業科技，2009(9)：185－188.

[9] 付強，梁川，楊廣林. 三江平原井灌水稻發展潛力探析[J]. 農業系統科學與綜合研究，2002，18(1)：23－26.

[10] 王浩，王建華，秦大庸，等. 基于二元水循環模式的水資源評價理論方法[J]. 水利學報，2006，37(12)：1496－1502. Wang Hao, Wang Jianhua, Qin Dayong, et al. Theory and methodology of water resources assessment based on dualistic water cycle model[J]. Journal of Hydraulic Engineering, 2006, 37(12): 1496－1502. (in Chinese with English abstract)

[11] 秦大庸，陸垂裕，劉家宏，等. 流域“自然-社會”二元水循環理論框架[J]. 科學通報，2014，59(4/5)：419－427.

[12] 謝新民，郭洪宇，唐克旺，等. 華北平原區地表水與地下水統一評價的二元耦合模型研究[J]. 水利學報，2002，33(12)：95－100. Xie Xinmin, Guo Hongyu, Tang Kewang, et al. Dual coupled model for integrated assessment of surface water and groundwater in North China Plain[J]. Journal of Hydraulic Engineering, 2002, 33(12): 95－100. (in Chinese with English abstract)

[13] 吳普特，高學睿，趙西寧，等. 實體水-虛擬水“二維三元”耦合流動理論基本框架[J]. 農業工程學報，2016，32(12)：1－10. Wu Pute, Gao Xuerui, Zhao Xining, et al. Framework of “two-dimension three-element” coupling flow of real water and virtual water[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(12): 1－10. (in Chinese with English abstract)

[14] 張俊娥，陸垂裕，秦大庸，等. 面向對象模塊化的分布式水文模型MODCYCLEⅡ：模型應用篇[J]. 水利學報，2012，43(11)：1287－1295. Zhang Jun'e, Lu Chuiyu, Qin Dayong, et al. MODCYCLE: An object oriented modularized hydrological model Ⅱ. Application[J]. Journal of Hydraulic Engineering, 2012, 43(11): 1287－1295. (in Chinese with English abstract)

[15] Beven K J, Kirkby M J. A physically based variable contributing area model of basin hydrology[J]. Hydrological Science Bulletin, 1979, 24(1): 43－69.

[16] 宮興龍，芮孝芳，付強，等. 基于格子玻爾茲曼方法的流域匯流數值模型[J]. 水力發電學報，2014，33(2)：19－25. Gong Xinglong, Rui Xiaofang, Fu Qiang, et al. A numerical model of watershed concentration based on lattice boltzmann method[J].Journal of Hydroelectric Engineering,2014, 33(2): 19－25. (in Chinese with English abstract)

[17] 宮興龍，付強，邢貞相，等. 基于格子玻爾茲曼法的TOPMODEL建模與應用[J]. 農業機械學報，2015，46(7)：181－186. Gong Xinglong, Fu Qiang, Xing Zhenxiang, et al. Modelling and application of TOPMODEL based on LBM[J]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(7): 181－186. (in Chinese with English abstract)

[18] 宮興龍，付強，王斌，等. 丘陵-平原-濕地復合區降雨徑流數值模型[J]. 農業機械學報，2016，47(11)：142－149. Gong Xinglong, Fu Qiang, Wang Bin, et al. Numerical model of rainfall runoff in hills-plain-wetland compound area[J]. Transactions of Chinese Society for Agricultural Machinery, 2016, 47(11): 142－149. (in Chinese with English abstract)

[19] 張東輝. 格子玻爾茲曼方法在水文學中的應用研究[D]. 南京：河海大學，2008. Zhang Donghui. Lattice Boltzmann Method and Some Application to Hydrology[D]. Nanjing: Hohai University, 2008. (in Chinese with English abstract)

[20] Allen R G, Smith M, Pereira L S, et al. An update for the calculation of reference evapotranspiration[J]. ICID Bulletin, 1994, 43(2): 35－92.

[21] Allen R G, Pereira L S, Raes D, et al. Crop evapotranspiration- guidelines for computing crop water requirements [M]//FAO Irrigation and Drainage Paper 56. Rome: FAO, 1998.

[22] 劉博，楊曉光，王式功. 東北地區主要糧食作物氣候生產潛力估算與分析[J]. 吉林農業科學，2012，37(3)：57－60. Liu Bo, Yang Xiaoguang, Wang Shigong. Calculation and analysis of climate productive potential of major grain crops in Northeast China[J]. Journal of Jilin Agricultural Sciences, 2012, 37(3): 57－60. (in Chinese with English abstract)

[23] 中國主要農作物需水量等值線圖協作組. 中國主要農作物需水量等值線圖研究[M]. 北京：中國農業科技出版社，1993：19－55.

[24] 劉鈺，汪林，倪廣恒，等. 中國主要作物灌溉需水量空間分布特征[J]. 農業工程學報，2009，25(12)：6－12. Liu Yu, Wang Lin, Ni Guangheng, et al. Spatial distribution characteristics of irrigation water requirement for main crops in China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2009, 25(12): 6－12. (in Chinese with English abstract)

[25] 劉鈺，Pereira L S. 對FAO推薦的作物系數計算方法的驗證[J]. 農業工程學報，2000，16(5)：26－30. Liu Yu, Pereira L S. Validation of FAO methods for estimating crop coefficients[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2000, 16(5): 26－30. (in Chinese with English abstract)

[26] 郭利丹，夏自強，林虹，等. 生態徑流評價中的Tennant法應用[J]. 生態學報，2009，29(4)：1788－1792. Guo Lidan, Xia Ziqiang, Lin Hong, et al. Researches on application of the tenant method in ecological flow evaluention[J]. Acta Ecologica Sinica, 2009, 29(4): 1788－1792. (in Chinese with English abstract)

[27] 彭世彰，艾麗坤，和玉璞，等. 稻田灌排耦合的水稻需水規律研究[J]. 水利學報，2014，45(3)：320－325. Peng Shizhang, Ai Likun, He Yupu, et al. Effect of irrigation and drainage coupling management on rice water requirement [J]. Journal of Hydraulic Engineering, 2014, 45(3): 320－325. (in Chinese with English abstract)

[28] 朱士江，孫愛華，張忠學. 三江平原不同灌溉模式水稻需水規律及水分利用效率試驗研究[J]. 節水灌溉，2009，28(11)：12－14. Zhu Shijiang, Sun Aihua, Zhang Zhongxue. Experimental research on water consumption law and water using efficiency of paddy under different irrigation mode in Sanjiang Plain[J]. Water Saving Irrigation, 2009, 28(11): 12－14. (in Chinese with English abstract) .

[29] 張惠斌，于東，姚章村. 論“打井種稻”與“循環經濟”[J].水利科技與經濟，2006，12(2)：819－821. Zhang Huibin, Yu Dong, Yao Zhangcun. Talk about “planting rice by constructing well for irrigation” and “cycle economy”[J]. Water Conservancy Science and Technology and Economy, 2006, 12(2): 819－821. (in Chinese with English abstract) .

[30] 韓麗偉，付強，劉東，等. 三江平原地下水承載能力綜合評價模型的構建及其應用[J]. 水土保持研究，2010，17(2)：182－185. Han Liwei, Fu Qiang, Liu Dong, et al. Construction and application of the model on groundwater resources carrying capacity in SanJiang plain[J]. Research of Soil and Water Conservation, 2010, 17(2): 182－185. (in Chinese with English abstract)

[31] 楊亞妹，王德文. 三江平原地下水資源開發潛力研究[J]. 水利科技與經濟，2011，17(8)：74－77. Yang Yamei, Wang Dewen. Developed potential study of groundwater resources in Sanjiang Plain[J].Water Conservancy Science and Technology and Economy, 2011, 17(8): 74－77. (in Chinese with English abstract)

[32] 尹喜霖，張烽龍，鄭春曉. 黑龍江省三江平原的地下水開發利用[J]. 黑龍江水專學報，2002，29(4)：26－33. Yin Xilin, Zhang Fenglong, Zheng Chunxiao. Underground water development and utilization in Sanjiang Plain in Heilongjiang province[J]. Journal of Heilongjiang Hydraulic Engineering College, 2002, 29(4): 26－33. (in Chinese with English abstract)

[33] 王斌. 黑龍江省糧食生產與耗水問題探討[J]. 節水灌溉，2015，34(12)：77－80. Wang Bin.Ａ study on grain production and its water consumption in HeiLongJing province[J]. Water Saving Irrigation, 2015, 34(12): 77－80. (in Chinese with English abstract)

[34] 雷志棟，蘇立寧，楊詩秀. 青銅峽灌區水土資源平衡分析的探討[J]. 水利學報，2002，33(6)：9－14. Lei Zhidong, Su Lining, Yang Shixiu. Balance analysis of water resources in Qingtongxia Irrigation Area[J]. Journal of Hydraulic Engineering, 2002, 33(6): 9－14. (in Chinese with English abstract)

Rice planting potential in plain-hill-wetland area estimated by nature-society water cycle model

Gong Xinglong, Fu Qiang※, Sun Aihua, Guan Yinghong, Wang Bin, Li Mo

(150030,)

With the increasing of well irrigation for rice since the middle of 20 Century, planting area of plain-hill-wetland region has expanded rapidly over its planting potential. And this phenomenon has caused a series of problems, such as the shrinking of wetland area, the decline of groundwater level and the deterioration of water quality. In order to prevent these problems, a two element water cycle model should be constructed to reveal the water cycle mechanism of plain-hill-wetland area under the influence of human activities, and to calculate the suitable area of rice planting. A grid-based distributed hydrological model was constructed by numerical method. The absolute value of relative error between fitted and measured flood peak discharge was in the range of 1.71%-11.55%. The determination coefficient of flood process line was in the range of 0.67-0.85. It indicated that the model could well simulate the hydrological cycle process in the study area. And then a two element water cycle model was further constructed by embedding the influence of human activities. According to the threshold of rice planting potential, parameters such as the maximum groundwater depth and the maximum groundwater level drawdown were extracted. And 18 kinds of scenarios were constructed by different combinations of the parameters and crop species in dry land. The planting potential of paddy field under various scenarios was calculated using the two element model under the condition of satisfying the control threshold, according to mean annual rainfall, evaporation and scenarios. It is observed that the maximum groundwater depth was 5.17-7.49 m, the maximum groundwater level drawdown was 1.67-3.73 m, the maximum slope of the paddy field was in the range of 0.028-0.053, the river water division flow was 50%-70%, the ratio of groundwater extraction was 79%-112%, the rice planting potential was 2.836×105hm2-5.412×105hm2, the area from the dry land was 2.105×105hm2- 4.032×105hm2, and the area from the unused land was 5.68×104-1.109×105hm2for the 18 kinds of scenarios. Taking scenario 17 as an example, the river ecological water demand, groundwater depth and evaporation of dry land crops were tested during the growth stage of paddy fields. In the whole rice growing season, the groundwater depth was less than 7.12 m, and the minimum value of base flow ratio in the developed paddy field was 33.45%. The ratio of mean evaporation from dry land of subdivisions in scenario 17 to that under the same land use scenario in 1990 was in the range of 0.98-1.05. It indicated that the effect of scenario 17 on the river water and groundwater level in the rice growth period was within the control range, while the effect on the dry field was relatively small. Hence, the scenario 17 was the optimal one for application. Currently, the parameters such as hydrometeorology, hydrogeology, irrigation system and groundwater depth were not available. However, the water cycle model constructed in the present study could provide valuable information for describing the hydrological cycle process in plain-hill-wetland area and seeking the potential of paddy field development, based on the full use of the limited parameters.

water cycle; evaporation; society; plain-hill-wetland area; rice planting potential; scenario; LBM

2018-07-30

2018-12-10

國家自然科學基金項目（51408107）；黑龍江省科學基金項目（E2017009）；農業水資源高效利用重點實驗室的開放基金（KF201501）

宮興龍，副教授，博士，主要從事流域水文模型研究。 Email：gongxl@neau.edu.cn

付強，教授，博士，主要從事水土資源高效利用研究。 Email：fuqiang@neau.edu.cn

10.11975/j.issn.1002-6819.2019.01.017

S273.4; S165+.27

A

1002-6819(2019)-01-0138-10

宮興龍，付 強，孫愛華，關英紅，王 斌，李 茉. 自然-社會水循環模型估算平原-丘陵-濕地區水稻種植潛力[J]. 農業工程學報，2019，35(1)：138－147. doi：10.11975/j.issn.1002-6819.2019.01.017 http://www.tcsae.org

Gong Xinglong, Fu Qiang, Sun Aihua, Guan Yinghong, Wang Bin, Li Mo. Rice planting potential in plain-hill-wetland area estimated by nature-society water cycle model [J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(1): 138－147. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2019.01.017 http://www.tcsae.org