水資源約束下區域耕地資源開發利用研究

孟祥玉，雷國平，孫曉兵，曲曉涵

(1.東北農業大學資源與環境學院，哈爾濱 150030；2.東北大學土地管理研究所，沈陽 110004)

0 引 言

耕地資源作為人類生存發展的物質基礎，其對于保障國家糧食安全具有重要的意義，科學合理的開發利用耕地資源對社會經濟與生態環境的健康可持續發展具有舉足輕重的地位[1]。水是農業的命脈，農業的發展、農業生態系統及糧食安全生產與農業水資源的供給能力直接相關，然而，農業水資源的高強度開采已出現嚴重的用水矛盾現象。水資源作為耕地資源有效利用、農作物健康生長的重要限制因素，影響耕地規模尤其是水田規模，如何有效的平衡水資源與耕地資源之間關系以實現耕地資源的可持續利用與保障國家糧食安全已成為耕地資源利用的關鍵。耕地資源的可持續利用保障社會經濟的持續穩定發展，而東北地區作為我國重要的糧食生產基地，對于保障我國糧食安全具有重要意義[2]?！秶壹Z食中長期規劃綱要(2008-2020年)》規劃對東北地區糧食生產提出了更高的要求，而優化土地資源利用結構與提升利用程度顯得意義非凡[3]。水稻是我國主要的口糧作物，在我國糧食生產中占有極其重要的地位，水稻的生產與發展直接關系到我國的糧食安全[4]。然而，在近幾十年的土地資源利用過程中，東北地區耕地資源利用呈現明顯的“旱改水”的水田化趨勢，并呈現由南向北、由西向東逐步擴張的趨勢[5]。水資源對于水田的擴張及其高效利用具有限制性作用，因而水田的規模變化將影響區域耕地規模的穩定性。

目前，針對水資源約束下適宜耕地規模的研究不多。主要是利用水量平衡原理、水熱平衡原理建立水資源與耕地面積的相關關系式計算區域耕地規模。如王韶華等以農業需水量和地下水可開采模數兩種方法，計算了三江平原適宜水稻種植的面積[6]。趙新風等依據水量平衡原理計算塔里木河流域不同水平年的最大灌溉面積[7]。趙清以建三江分局為研究區在分析了水資源開發利用情況的基礎上，通過農業需水量的方法對當地水田規模進行的估算[8]。向長玉等基于GIS技術，考慮自然約束條件，應用定量化模型計算黑龍江省通河縣的最大水田面積[9]。郝麗娜等利用水熱平衡原理建立適宜規模模型確定黑河干流中游地區耕地規模[10]。鄧寶山等在水熱平衡理論基礎上分析了吐魯番綠洲適宜規模及其穩定性[11]。關于水資源約束條件下耕地規模的研究主要都集中在西北干旱區，東北區域研究較少，還沒有意識水土資源的稀缺性、有限性、地域性決定了農業發展和糧食安全只能依靠現有水土資源的合理利用，而針對較小區域的耕地規模的研究就更少。

黑龍江省作為我國重要的商品糧生產基地，耕地資源利用雖呈現多樣性，水田的面積大幅度增加，旱地快速減少已成為其耕地資源利用的主要特征。本文利用水稻的種植面積計算水田的面積[12]，由于農場現有水田面積已占耕地面積的90.27%，水田是農場耕地開發利用的主要方向，因此本文將以計算農場水田規模為耕地規模。大規模的水稻種植必須有充足的水源。因此本文以三江平原寒季水稻種植第一場的位于黑龍江農墾總局建三江分局青龍山農場為研究區，基于遙感影像數據對其耕地利用時空變化特征分析的基礎上，根據水量平衡的原理，探究水資源約束下耕地的適宜規模，以期為水土資源的高效利用，發展適度規模耕地，農業種植結構調整及其耕地資源的可持續發展提供參考。

1 研究區概況

青龍山農場位于黑龍江省三江平原東北部，隸屬于黑龍江省農墾總局建三江分局，地理坐標132°53′～133°55′E， 47°47′～48°13′N，北靠黑龍江，西臨青龍蓮花河。土地面積為601.33 km2，為建三江分局土地面積的4.85%，其地理位置如圖1所示。其氣候屬于寒溫帶大陸性季風氣候，夏季高溫多雨，冬季寒冷干燥，全年溫差大，年均氣溫為2.0 ℃，年氣溫大于10 ℃以上有效積溫為2 350.3 ℃，年均日照時數約為2 600 h，年均無霜期為127 d，年均蒸發量1 159.6 mm。年均降水量約為550 mm，且降水時空分布不均，多集中在6-9月，地表年徑流量為0.48 億m3，但由區域水資源利用現狀可知，利用率極低，幾乎不利用，在區域這種降水量時空分布不均，地表徑流水利用率極低的狀況下，為保證農業的高產穩產，利用地下水灌溉發展水田是主要手段。區域內地形低平，地勢南高北低，東高西低，地面坡降為1/7 000～1/1 000。土壤類型主要有白漿土、草甸土、沼澤土及泛濫地土壤，土壤肥力較高，適宜農業種植。

圖1 青龍山農場地理位置示意圖

作為寒地水稻的發源地、三江平原水稻種植第一場的農場，青龍山農場經濟發展以農業為主，耕地資源的開發利用因農業發展及種植業結構調整而經歷了很大的變化。農場主要種植水稻、大豆、玉米等作物，由于水稻畝產量及其收益高的優勢及國家政策支持，已成為農場耕地的主導作物。農業生產綜合機械化率96.8%。周邊地表水較為豐富，但由于其缺乏水源工程，田間配套工程使得其開發利用率極低。由于農民可根據自己意愿隨意打井開采地下水并采用漫灌方式種植作物，機電井的數量在2000-2014年間就由679眼增加到1 590眼，增加了134%，地下水嚴重超采并伴隨著用水浪費現象。隨著農場水稻種植面積的增加，農業用水也不斷增加，引發了農業用水與生活用水之間的矛盾。

2 數據來源與研究方法

2.1 數據來源

(1)遙感影像數據。本文主要以全國第二次土地調查數據為基礎數據，基于Landsat TM/OLI遙感影像為數據源，獲取青龍山農場2000、2014年土地利用數據集。通過選取2000和2014年云量較少且覆蓋全境的Landsat遙感影像，選取其分辨率為30 m的綠色、紅色、近紅外波段合成30 m分辨率的標準假彩色影像，經過幾何校正、圖像增強等過程，獲取2000、2014年遙感影像數據集?；贏rcgis10.2軟件平臺以二調數據為參考，根據不同的地類的色調、紋理等特征進行人工目視解譯，提取2000、2014年土地利用數據，根據野外實地勘測檢驗其精度在90%以上。由遙感解譯獲得水田與旱地的空間分布數據。

(2)水利數據。水資源作為水田開發利用的重要限制性因素，本文主要利用參數有降水滲入系數、滲透系數、垂直于計算剖面的水利坡度、含水層厚度、河道或河段長度、灌水定額、入滲系數、洪水淹沒區地表水深度等數據。其中，數據主要源自建三江管局水利局《建三江水利志》。

(3)其他數據。灌溉定額數據通過實地調查，結合建三江信息港、建三江管局寒地優質高產水稻生產技術規程、黑龍江省用水定額標準(DB23/T 727-2003)獲得。利用灌溉定額計算區域耕地規模。降雨量等氣象數據則主要來自青龍山農場氣象局。農場信息主要來自《黑龍江史志》、《建三江農墾志》、《黑龍江農墾地名錄》、《青龍山農場史志》、《黑龍江墾區統計年鑒》(2001-2015年)、《黑龍江省水文地質志》等統計數據。

2.2 研究方法

分析區域耕地內部結構變化，利用水田旱地的種植比例計算農業需水量，與降水量比較求出基于降水條件下區域內耕地適宜規模。根據公式計算地下水可開采量，利用灌溉定額計算出基于地下水可采量下耕地的適宜規模，將計算出的耕地面積與實際的耕地面積進行比較，得出區域耕地的超載情況。利用區域內徑流量結合自流提水灌溉定額計算地表水條件下的耕地規模。

2.2.1 農業需水量

青龍山農場的主要作物是水稻，水稻的生長要經歷5-10月的生育期，11-4月的休閑期。根據相關實驗數據及文獻資料可知各主要作物需水量[13,14]。其中大豆的需水量為498 mm，玉米的需水量為478 mm，水稻的需水量為660 mm。在氣候溫度及緯度的影響下，東北區域農作物均為一年一熟，11-4月內蒸發量很小，水稻面積比例高于旱地且需水量大，旱作物基本不灌溉，故可將旱作物及非耕地上植物的蒸散量用旱作物需水量的均值替代，計算作物需水量[15]。農場旱地作物是由農民意愿隨意種植，可將玉米，大豆按輪作的方式進行種植，由此可定旱田的需水量為488 mm。

農業需水主要是指區域土地1 a所需用水量(包括非耕地與非生育期生態用水)[16]。利用農作物的需水量來表征耕地的需水狀況，結合農作物的播種面積計算出農業需水量。即：

W農=W旱(1-Z)+W水Z

(1)

式中：W農為農業需水量，mm；W旱為旱作物平均需水量，mm；W水為水稻需水量，mm；Z為水稻種植比例。

目前針對耕地規模的已有不少學者開始研究，其中根據向長玉等[9]可知利用以下公式即可求得區域內耕地規模，公式如下：

(2)

式中:S水為水田面積，km2；Q總為總耗水量，mm/km2；Q旱為旱地單位面積耗水量，mm/km2；S總為區域耕地面積，km2。

2.2.2 地下水補給量

在4月中旬開始水稻要進行棚內育秧及泡田整地，為抗低溫及以水泡田化凍的作用。到5月進行秧苗移栽，6月初基本完成稻苗移栽進地工作。7月齊穗期需要陸續對水田供應充足的水分，而到9月成熟期、10月末水稻收獲時卻不需要水資源，與之相反的是越干熱越好，為便于機械收割。降水主要集中在6-8月，在水稻的整個生育期內，由于降水的時間不適宜及空間不均性，常有春旱秋澇的現象發生。地表水資源控制性工程的缺乏，青龍山農場水田主要依靠地下水進行農業灌溉。利用地下水的補給量求得地下水的開采量[17]，當地下水的開采量超過地下水補給量，造成地下水水位下降。因此地下水的開采量應控制在地下水的補給量范圍內[18]。地下水補給量計算公式[19-21]如下：

Q綜補=Q降補+Q河補+Q地側+Q灌補+Q洪補

(3)

(1)降水入滲補給量是指降水滲入到土壤中并在重力作用下滲透補給地下水的水量。

Q降補=αPF×10-5

(4)

式中：Q降補為降水滲入補給量，億m3/a；α為降水滲入系數；P為多年平均降水，mm;F為降水滲入補給面積，km2。

(2)河道地表水滲入補給量。河道滲漏補給量，是指當河道水位高于河道岸邊地下水水位時，河水滲漏補給地下水的水量。

Q河渠=KIHLt×10-8

(5)

式中：Q河渠為河渠地表水滲入補給量，億m3/a；K為滲透系數，m/d；I為垂直于計算剖面的水利坡度；H為含水層厚度；L為河道或河段長度；t為河道滲入補給時間，d/a，取365 d。

(3)地下水側向滲入量補給指山前地下水側向徑流補給本區淺層地下水的水量。

Q地側=KIHLt×10-8

(6)

式中：Q地側為地下水側向滲入量補給量，億m3/a。

(4)井灌回歸補給量是指井灌水(系淺層地下水) 進入田間后，入滲補給地下水的水量。計算公式如下：

Q井灌=Q田間β井

(7)

Q田間=qF×10-8

(8)

式中：Q井灌為井灌回歸補給量，億m3/a；Q田間為井灌田間灌水量，億m3/a；β井為井灌回歸系數，取β井=0.1；q為灌水定額，m3/a；F為水田面積，hm2。

(5)洪水淹沒補給量是指水淹沒期，水體滲入補給地下水的水量。

Q洪補=αFΔH×10-2

(9)

式中：Q洪補為洪水淹沒下補給量，億m3/a；α為入滲系數，取α=0.25；ΔH為洪水淹沒區地表水深度，m；F為計算面積，km2。

2.2.3 地下水可開采量

地下水可開采量是指在技術上可能，經濟上合理和不造成水質惡化，水位持續下降及其他不良后果的條件下，能夠從取水地段中取得的水量[18]。其計算公式如下：

Q開=Q補給P

(10)

式中：Q開為地下水可采量，億m3/a；Q補給為地下水總補給量，億m3/a；P為開采系數，取開采系數P為0.80。

2.2.4 耕地規模計算

根據區域內的水資源量，利用區域灌溉定額計算得到耕地規模。其公式為：

A=W/I

(11)

式中：A為最大可承載面積，km2；W為農業用水量，億m3；I為灌溉定額，m3/km2。

3 結果與分析

3.1 區域耕地變化特征分析

耕地是黑龍江墾區的重要景觀基礎。隨著農業現代化進步與不斷深入，種植業結構調整致使青龍山農場耕地資源利用變化顯著。2000-2014年青龍山農場耕地利用類型發生較大的變化，且呈現出明顯的時間變異性。2000年青龍上農場耕地面積為447.94 km2，土地墾殖率為81.58%，其中，水田面積為103.55 km2，所占比重為23.12%，旱地面積為344.39 km2，所占比重為76.88%。2014年青龍上農場耕地面積為468.84 km2，土地墾殖率為85.38%，耕地面積增長率為4.67%，然而，水田面積為423.21 km2，所占比重為90.27%，增長率為308.70%，年均增長率為22.05%，旱地面積為45.64 km2，所占比重為9.73%，減少比率為86.75%，年均減少比率為6.20%(見圖2)。青龍山耕地利用變化體現出區域經濟發展以農業為主，且水稻已成為區域農業的主要種植作物。然而，隨著耗水作物水稻種植規模的增加，農場的農業水資源量不斷增大，農場早就出現用水矛盾問題。農業用水期間生活用水停水即農業用水占用生活用水及地下水水位下降現象。由此可見，水田規模不能任意擴張，種植面積過大，耗水多，會發生各個行業用水矛盾問題，然而，種植面積過小，會發生限制經濟發展，不能充分利用區域資源問題。

2000-2014年青龍山農場水田面積大幅度增加，且耕地利用狀況呈現出明顯的空間差異性(見表1)。2000年青龍山農場各管理區水田平均面積為14.79 km2，所占比例主要集中于18%～30%之間，其中，第五管理區水田面積較大，所占比例為54.37%，而第四管理區水田面積最小，所占比例僅為12.51%，而2014年各管理區水田平均面積為60.46 km2，除第四管理區以外其余管理區水田所占比例均在92%以上。2000-2014年青龍山農場旱地所占比例快速降低，2014年第四管理區旱地所占比例為37.94%，其余的管理區旱地所占比例均在10%以下。由此可知，2000-2014年青龍山農場耕地利用狀況發生巨大的變化，旱地與水田之間的轉化較為劇烈，水田明顯增加，水資源的合理高效利用對于水田發展具有至關重要的作用。

圖2 青龍山農場耕地利用分布圖

表1 2000-2014年青龍山農場耕地利用變化特征

注：表中一、二、三、四、五、六、七分別為青龍山農場內第一管理區、第二管理區、第三管理區、第四管理區、第五管理區、第六管理區、第七管理區。

3.2 區域水資源

2000-2014年間，青龍山農場耕地規?？傮w變化不大，總體上水田規模變化突出。目前，青龍山農場農業種植結構主要以水田為主，截止到2014年，水田面積達到423.21 km2，占現有耕地面積468.84 km2的90.27%。水稻是農場作物種植結構中耗水最大的作物，因此要有充足的水源保障才能支持水田的發展。農場水資源量由降水、地表水及地下水組成。地表水主要來源于降水及河流形成的地表徑流。由于降水的時空分布不均性及水利工程的缺乏，青龍山農場用于農業灌溉用水量來源于地下水。

根據公式(3)～(9)，可計算得出區域相關補給量如表2所示。由表2可知，青龍山農場的綜合補給量為1.427 億m3。繼續利用公式(10)可計算出區域內可開采地下水量為1.414 億m3。利用地下水開采模數法并結合區域給水度可算出地下水水位變幅情況[6]，計算區域內地下水水位變幅約為11 m與青龍山農場多年平均地下水位的降深實際情況相吻合，故計算結果較為合理。近些年來，由學者[22,23]對建三江分局地下會可開采量的研究可知農場2000及2006年的地下水補給量及推算可開采量，如表3所示。由表3可知，近14a來，農場地下水總補給量增加，可開采的地下水量也相應地增加。地下水補給量由2000年0.57 億m3到2014年增加到1.427 億m3?？砷_采量也有2000年0.555 億m3增加到2014年的1.141 億m3。從2006-2014年的年均用水量高于2000-2006年的年均用水量。和青龍山農場實際的地下水開采量相比較，地下水均已超采，因此應嚴格控制對地下水的開采。地下水量應先滿足生活、牲畜用水及工業用水后再用于農業灌溉，由于農場工業發展緩慢，用水極低，故本文不考慮工業用水。根據黑龍江省用水定額標準(DB23/T727-2003)可知，生活用灌溉定額人均70L/(人·d)，牲畜用水按豬的用水定額計算為40L/(頭·d)。因此農業地下水量見表6。相比2000年農業地下水量到2014年農業地下水量增長迅猛。

表3青龍山農場地下水可開采量與用水量對比表億m3/a

徑流量主要受降水量的影響。青龍山農場內多年平均徑流量為0.479 億m3。由表4 可知，在不同降水保證率條件下，其徑流量不同。扣除河道內生態用水后即可算出可利用量。根據相關學者對河道內生態需水研究[24,27]可知，河道內生態需水多取多年平均徑流量的20%或40%，20%為最小河道內生態徑流，40%為最適宜河道內生態徑流。根據不同保證率下的徑流總量可估算出地表水可利用量。降水保證率為50%時，其地表徑流量最大0.303 億m3，可利用水量為0.182～0.242 億m3，降水保證率為75%時，地表徑流量最小為0.014 億m3，可利用水量為0.065～0.087 億m3，如表4所示。

表4 地表徑流條件下耕地適宜規模

3.3 區域耕地規模計算

青龍山農場用于耕地灌溉的只有水田，旱地僅用于天然降水即可滿足需求。近些年來，由于水稻單產量、經濟利益高，國家政策的鼓勵等因素的影響，水田規模急速擴增，農業耗用水量也隨著增加。因此，本文地下水的開采量及灌溉定額所求耕地規模為水田規模。

3.3.1 降水條件下耕地規模計算

青龍山農場多年平均降水量約為550mm，降水量能夠滿足旱作物全年的需水量要求，能夠滿足一定規模的水田作物的生長，本文主要通過計算區域農業需水量與多年平均降水量比較，分析降水資源條件下區域內耕地規模尤其是水田規模。

根據公式(1)結合耕地內部不同比例的種植結構計算得到農業需水量如表5所示。由表5可知，隨水稻種植面積的增加農業需水量增大。當水稻種植比例在30%～40%時，農業需水量為539.60～556.80mm，種植比例為90%時，農業需水量為642.80mm。由此可知，到2014年青龍山農場水稻種植比例為92.29%，農業需水量要高于642.80mm，超出多年平均降水量至少98.80mm。考慮地表徑流量即生態環境用水，根據相關文獻資料[28,29]，地表徑流量取降雨量的20%較為合理?？鄢乇韽搅髁考s110mm可利用降水量約為440m，此時已少于農業需水量約202.80mm。由此可知，降水量不能滿足農業需水量。由此可判斷，在當前的降水條件下，青龍山農場耕地規模是超載的，緊靠降水量已無法承載更大規模的耕地開發利用。當水稻種植比例為30%～40%時，農業需水量接近多年平均降水量550mm，此時水稻種植面積介于140.65～187.54km2之間，但此時已充分利用了降水量，對生態環境的可持續發展勢必會有所影響。若水稻種植超過這個區間，便只能利用地表水或開采利用地下水。根據公式(2)可知，區域內適宜耕地規模約為221.15～276.43km2。與現有耕地面積468.84km2相比可知青龍山農場耕地已嚴重超限。超限面積約為192.41～247.69km2。

表5 農業需水量估算

3.3.2 地下水條件下耕地規模計算

地下水開采量過大，就會造成地下水超采，引發供水矛盾。而地下水開采量過小，就難以充分利用地下水資源，制約農業經濟的快速發展。利用機電井開采地下水因其開采工程簡單、利用方便、費用較低的優勢已成為該地區農業主要灌溉水源。在水田種植面積不斷增大的背景下，地下水資源超采現象日益凸顯。根據收集資料顯示農業生產用水占地下水的總用水量比例高達99.5%，而用于生活、工業的總用水量不足地下水總用水量的1%?？芍帻埳睫r場水田主要依靠地下水進行灌溉，因此了解區域內利用地下水灌溉水田規模尤為重要。

地下水總量應在先滿足人類生活及工業用水的基礎上，再用于農業生產，根據表3數據可知，農用地下水量占地下水總量比例很大，且2000-2014年間，農用水用總量逐年增大。根據建三江管局寒地優質高產水稻生產技術規程、黑龍江省用水定額標準(DB23/T727-2003)等相關研究[30]獲得青龍山農場灌溉定額(取灌溉定額為59.997 萬m3/km2)，利用公式(10)計算出耕地規模，見表6。

表6 地下水條件下耕地規模

由計算結果可知，2000-2014年適宜耕地規模明顯增加，耕地年增長速率為6.12km/a，根據區域實際地下水用于農業情況可知此時利用地下水進行計算耕地規模均為水田規模。2000-2006年間耕地規模由91.54km2擴大到116.64km2，耕地年增加速率為4.18km/a。2006-2014年，適宜耕地規模擴張到177.58km2，耕地年增加速率為7.62km/a。2006-2014的耕地增長速率明顯高于2000-2006年間的耕地增長速率，說明在近些年間，水田開發規模增長迅速。與青龍山農場實際水田規模相比較可知，2000年區域耕地適宜規模與實際規模較為均衡。2006年，農場水田規模超載約為36.69km2。到2014年耕地規模已嚴重超載，超載面積約245.63km2。由此，可知區域內超載耕地面積均來源于水田面積的超載。因此，應嚴格控制農場內水田規模的擴增，避免地下水超采利用引起生態環境惡化，資源枯竭的狀況發生。

3.3.3 地表徑流條件下耕地規模估算

青龍山農場地表徑流量與降水量存在著正相關關系。降雨量增加地表徑流量也會相應增加。由于其引水工程的缺乏，地表徑流水利用率很低幾乎不利用，如若能充分利用地表徑流，還能在一定程度上擴大耕地的規模。因此，針對地表徑流數據估算可利用水量，利用公式(10)進而估算出地表徑流條件下耕地規模。取灌溉定額取74.996 3 萬m3/km2，地表徑流條件下的耕地規模(見表4)。當充分利用地表徑流時，在保證率為50%時，耕地規模在最小地表徑流和最適宜地表徑流條件下，利用地表水時耕地規模約為24～32km2，在保證率為75%時，耕地規模在最小徑流和適宜地表徑流條件下，利用地表水時耕地規模約為8.67～11.60km2。此時僅為利用青龍山農場區域內的地表徑流，因數據不全，未計算利用周邊區域的河流引用水量計算區域內耕地規模。

4 結 論

(1)2000-2014年青龍山農場耕地總體面積變化不大，耕地內種植結構發生劇烈變化，水田面積大幅度擴張，旱地急劇減少，且耕地利用變化呈現出明顯的時空差異性。耕地面積由2000年447.94km2發展到2014年468.84km2。水田面積占耕地比重從2000年的23.12%發展到2014年的90.27%。而旱地從占耕地的76.88%降低為9.73%。2000年水田零星分布，旱地集中分布。2014年水田旱地轉變劇烈，水田分布廣泛，且集中連片分布于區域各個管理區，其中，第五管理區水田面積最大，而第四管理區水田面積最小。2000-2014年青龍山農場旱地所占比例快速降低，2014年第四管理區中旱地所占比例最大為37.94%，其余管理區中旱地所占比例均在10%以下。

(2)隨著青龍山農場耕地規模的增加，農業需水量不斷增加，依靠降水、地下水已遠遠不能滿足農業需水的需求。降水量由于其時空分布不均性僅能滿足旱地的需求，受降水資源影響，2000-2014年，青龍山農場地下水總補給量增加，區域地下水用量增加，農業用地下水量也不斷在增加，且農業用地下水量占用地下水總量比例最大。地下水實際開采量均已超采，因此，應嚴格控制地下水的超采。青龍山農場區域內有青龍蓮花河，含有較為豐富的地表徑流，但是由于水利引水工程等的缺乏，利用率極低。白白浪費了地表水。地表徑流量隨著降水保證率增加而減少。在降水保證率為50%時，扣除河道內最小及最適宜生態凈流量后可利用水量為0.182～0.242 億m3，降水保證率為75%時，扣除生態徑流量后可利用水量為0.065～0.087 億m3。

(3)青龍山農場水資源主要由降水，地下水和地表徑流組成，因此，本文主要研究這三種用水條件下的耕地規模。當水稻種植比例在30%～40%時，農業需水量為539.60～556.80mm，與多年平均降水量550mm相近。降水條件下適宜耕地規模約為221.15～276.43km2。與實際規模相比可知實際耕地規模已超限，且超限規模約為192.41～247.69km2。利用地下水計算的耕地規模主要是水田規模，利用地下水的開采量及灌溉定額計算區域耕地規模，并與實際耕地規模比較，得到耕地超載情況。2000-2014年適宜耕地規模明顯增加，2000-2006年間耕地規模由91.54km2擴大到116.64km2，耕地年增加速率為4.18km/a。2006-2014年，適宜耕地規模擴張到177.58km2，耕地年增加速率為7.62km/a。2006-2014的耕地增長速率明顯高于2000-2006年間的耕地增長速率，2000及2006年均算較為平衡，且到2014年耕地規模已嚴重超載。且超采規模由2006年36.69km2增加到2014年的245.63km2。區域內耕地超限規模均來自水田規模的超限，因此應嚴格控制利用開采地下水擴大水田規模。在去除最小河道內生態徑流量、最適宜河道內生態徑流量后得到的河流可利用量利用灌溉定額計算得到耕地規模。在降水保證率為50%，取最小河道內生態徑流量后的地表水可利用水時，耕地規模約為32km2。在降水保證率為75%時，耕地規模約為11.60km2，此時利用地表徑流估算得到的耕地規模均為區域內地表河流，由于無法獲知區域周圍河流水引用量數據，就沒有計算這部分水資源量下的耕地規模。但這確實是區域耕地規模發展的關鍵，應大力發展這部分水資源的利用。

近十多年來，青龍山農場作為三江平原水田擴張的典型區域，其土地利用結構發生了劇烈的變化，水田化已成為三江平原農業現代化發展的趨勢[31]。從宏觀層面，種植業結構調整致使水稻種植面積大幅度增加，對保障國家糧食安全及維持社會穩定具有積極的作用，從微觀層面，種植水稻的經濟驅動促進農戶開發利用水田，有效地提升了耕地利用效率。然而，高耗水作物水稻的種植面積的擴大使得農業用水大量增加，致使水資源利用成為限制其水田可持續發展利用的關鍵。但已有水田開發利用大多是以區域水資源的消耗和地下水位下降為代價，阻礙區域水資源的合理可持續利用。因此，科學合理的開發利用水資源應優先開發利用地表水，控制開采地下水，充分利用區域內地表徑流，大力開發利用地表水及降水資源，提高水資源利用率。區域北靠黑龍江，應大力發展水利引用江水灌溉，可以進一步開發利用耕地?？茖W的在水田面積快速擴張的過程中有效的利用水資源并實現水土資源優化配置，進而推進區域水土資源及社會經濟的可持續發展已成為農業現代化發展的關鍵。

[1] 滿衛東.1990-2013 年東北地區耕地時空變化遙感分析[J].農業工程學報,2016,32(7):1-8.

[2] 周廣云,郭微微.黑龍江省“千億斤糧食產能工程”建設水資源支撐探討[J].安徽農學通報,2014,20(16):1-7.

[3] 于鳳榮,劉 義,劉鳳山.基于RS的友誼農場近10年來水田時空變化分析[J].現代化農業,2012,(11),59-61.

[4] 章秀福,王丹英,方福平.中國糧食安全和水稻生產[J].農業現代化研究,2005,26(2):85-88.

[5] 劉紀遠,張增祥,莊大方.20世紀90年代中國土地利用變化時空特征及其成因分析[J].地理研究,2003,22(1): 1-12.

[6] 王韶華,劉文朝,劉群昌.三江平原農業需水量及適宜水稻種植面積的研究[J].農業工程學報,2004,20(4):50-53.

[7] 趙新風，徐海亮，王 敏，等. 不同水平年塔里木河流域灌溉面積超載分析[J]. 農業工程學報，2015，31(24)77-81.

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

[9] 向長玉,周東興,宋 戈,等.黑龍江省通河縣旱地改水田許可轉換度研究[J].農業工程學報,2012,28(13): 226-230.

[10] 郝麗娜，粟曉玲.黑河干流中游地區適宜綠洲及耕地規模確定[J].農業工程學報，2015,31(10)；262-268.

[11] 鄧寶山，瓦哈甫·哈力克，張玉萍，等. 吐魯番綠洲適宜規模及其穩定性分析[J].干旱區研究，2015,32(4):797-803.

[12] 陶軍德,李艷芳,曾光建.近30年黑龍江省水田規模與分布變化及其對氣候變化的響應[J].中國土地科學, 2011,25(4):26-30.

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

[14] 周 浩,雷國平,張 博.1990-2013 年撓力河流域耕地變化下水土資源平衡效應分析[J].農業工程學報,2015,31(1):272-279.

[15] 王韶華,劉文朝,劉群昌.三江平原農業需水量及適宜水稻種植面積的研究[J].農業工程學報,2004,20(4):50-53.

[16] 田 園. 田園水利文集[M].北京:中國水利水電出版社,1998.

[17] 吳曾生.地下水水文學[M]. 北京：水利水電出版社，1985.

[18] 宋傳英,趙桂芳,孟祥全.利用地下水灌溉農田的水量供需平衡計算[J].黑龍江水專學報,2003,30(1):61-62.

[19] 王國慶,劉冬梅.地下水補給量數學模型解析[J]. 黑龍江水利科技,2012,40(3):166-167.

[20] 梁 冰.水文地質參數確定與地下水資源量計算方法分析[J].黑龍江水利科技,2014,42(12):120.

[21] 孔祥斌,張鳳榮,齊 偉.集約化農區土地利用變化對水資源的影響----以河北省曲周縣為例[J].自然資源學報,2004,19(6):747-753.

[22] 付 強，佟兆生，魏慶利.建三江分局地下水資源計算與評價[J].農機化研究,2002,2(1):47-49.

[23] 李文超.建三江墾區地下水資源評價與預測研究[D]. 哈爾濱：東北農業大學，2007.

[24] 李 捷，夏自強，馬廣慧，等.河流生態徑流計算的逐月頻率計算法[J].生態學報.2007,27(7):2 916-2 921.

[25] 馮夏清,章光新,尹雄瑞.基于生態保護目標的太子河下游河道生態需水量計算[J].環境科學學報，2010,30(7):1 466-1 471.

[26] 李麗娟,鄭紅星.海灤河流域河流系統生態環境需水量計算[J].地理學報,2000,(4):495-501.

[27] 倪晉仁,崔樹彬,李天宏，等. 2002a論河流生態環境需水[J].水利學報,2002,(9):14-26.

[28] 陳小兵,楊勁松,喬曉英.綠洲耕地適宜面積確定與減災研究----以新疆渭干河灌區為例[J].中國地質災害與防治學報,2008,19(1):118-123.

[29] 曹志超,王新平,李衛紅.基于水熱平衡原理的塔里木河下游綠洲適宜規模分析[J].干旱區地理,2012,35(5):807-814.

[30] 魏 然,李 智.青龍山農場灌區灌溉制度分區概述[J].水利科技與經濟，2010,16(2):188-189.

[31] 杜國明,潘 濤,尹哲睿.水田擴張中三江平原北部耕地景觀演化規律及驅動因素分析[J].應用生態學報,2015,26(1):207-214.