1. 项目概况
普洱市是实行最严格水资源管理制度及水生态文明建设试点城市，已被确定为代表云南省迎接国务院2016年实行最严格水资源管理制度工作考核的州市之一 [1] 。目前，全市最严格水资源管理各项相关工作正在有序进行中。但农田灌溉水有效利用系数作为最严格水资源管理制度考核的重点内容之一，也是普洱市绿色经济考核的重要指标，仍缺乏实测资料支撑。基于以上背景，为提高农业用水效率，将最严格水资源管理制度考核的具体指标落到实处，开展普洱市主要作物需水量试验及不同灌区灌溉水有效利用系数的实测分析工作，提出普洱市主要粮经作物需水量、作物系数、灌溉定额及灌溉制度，摸清全市各区域及各类型灌区灌溉水有效利用系数的现状，为推动全市农业节水管理、制定节水规划、合理评价节水灌溉发展成效及分析节水潜力提供重要基础，同时对全省农田灌溉水有效利用系数的测算起到重要示范和推动作用。
2. 技术流程
开展普洱市农田灌溉水有效利用系数测算和主要作物需水量试验工作技术流程如下：
第一，在对全市灌区综合调研的基础上，分类汇总各类型、不同档次灌区的灌溉面积、工程设施与用水状况等，分析全市样点灌区的代表性，确定各规模与类型样点灌区，结合灌区的作物种植面积等，在样点灌区内选择典型田块，在典型田块开展灌溉试验及灌溉水有效利用系数测算工作。
第二，同步开展主要作物需水量试验灌溉定额，试验中旱作的土壤墒情由智能云墒情监测仪监测记录，稻田水位数据由水位计监测记录。由土壤含水率及水位的变化可灌溉定额计算试验区净灌溉定额。根据灌溉定额样点灌区年度净灌溉定额与实测的样点灌区毛灌溉定额计算得到样点灌区灌溉水有效利用系数，并依次方法测算市级区域不同规模与类型灌区的灌溉水有效利用系数。
第三，通过气象要素监测、土壤墒情监测、水位数据监测、农事活动记录等资料，分析计算试验区相应作物的需水量、灌溉定额及作物系数。
第四，对试验所得灌溉水有效利用系数及主要作物需水量进行综合分析，提出提高普洱市灌溉水有效利用系数的对策建议，提出普洱市主要作物的适宜灌溉制度。具体流程图如图1所示。
3. 方案设计
3.1. 灌区规模及用水情况
3.1.1. 灌区规模
截至2015年，普洱市有设计灌溉面积在1.0~30.0万亩之间的中型灌区10个，其中有设计灌溉面积在1.0~5.0万亩之间的中型灌区3个，5.0~15.0万亩之间的中型灌区7个；全市中型灌区设计灌溉面积总和60.91万亩，有效灌溉面积总和49.51万亩。100亩以上1万亩以下的小型灌区有163个，小型灌区设计灌溉面积总和134.6万亩，有效灌溉面积总和113.9万亩，灌区主要种植作物有水稻、玉米、烤烟、茶叶、咖啡、蔬菜等，灌水方式均为自流。
3.1.2. 灌溉用水情况
据普洱市2015年水资源公报统计，2015年普洱市农业灌溉用水量为9.2492亿m3，农业用水消耗量为6.2326亿m3。目前普洱市灌溉水有效利用系数采用经验值，介于0.436至0.550之间，平均为0.498，最低为西盟县0.436，最高为景谷县0.550。因此，测定普洱市灌溉水有效利用系数非常必要。
3.2. 样点灌区选取
3.2.1. 样点灌区选择原则
根据《全国农田灌溉水有效利用系数测算分析技术指导细则》(以下简称《细则》)的规定，样点灌区的选择应具有代表性、可行性及稳定性等原则。在样点灌区、典型田块及外业监测过程中要综合考虑本地区不同规模与类型灌区的分布、各类型灌区配套水平、作物种植结构等情况，力求使农田灌溉水有效利用系数测算分析成果能够真实反映当地灌区整体特点及农业用水水平。其中，中型灌区按有效灌溉面积大小分为3个档次，即1万亩 ≤ A中型 < 5万亩、5万亩 ≤ A中型 < 15万亩、15万亩 ≤ A中型 < 30万亩，每个档次的样点灌区数量
Figure 1. Flow chart of effective utilization of irrigation water for farmland and major crop water demand test techniques
图1. 农田灌溉水有效利用系数测算和主要作物需水量试验技术工作流程图
不应少于本级区域相应档次灌区总数的5%，各档次样点灌区有效灌溉面积不应少于本级区域相应档次灌区有效灌溉面积的10%；小型灌区样点灌区数量按不少于全市小型灌区取样范围内(100亩~10000亩)数量的0.5%，有效灌溉面积不小于该市级区域全部小型灌区有效灌溉面积的1%为原则 [2] 。
3.2.2. 样点灌区的选取结果
通过对普洱市水资源开发利用情况及各类型灌区的现场查勘及调研，了解灌区的自然地理概况及其分布、灌区水利工程设施、灌溉水源和灌溉用水量、有效灌溉面积、种植结构等实际情况，收集了灌区灌水及其他社会经济统计数据、灌区平面布置图等基础资料(基础资料收集以2015年为现状年)。依据云南省用水定额分区，普洱市为滇西南区(III区)，其中属III-1区的有江城县、孟连县、西盟县3个县；属III-2区的有澜沧县1个县；属III-3区的有思茅区、宁洱县、景谷县、墨江县、景东县、镇沅县6个县。按照样点灌区选择原则并结合普洱市灌区实际情况，选取中型灌区5个，全属于自流灌区，所选样点灌区有效灌溉面积占全市中型灌区有效灌溉面积总和的56.69%；综合考虑取水类型及空间分布，选择5个小型样点灌区，小型样点灌区有效灌溉面积占普洱市小型灌区有效灌溉面积总和的2.04%。样点灌区选取结果见表1。
3.2.3. 典型田块的选取
典型田块的选取要求边界清楚、形状规则、面积适中；综合考虑作物种类、灌溉方式、畦田规格、地形、土地平整程度、土壤类型、灌溉制度与方法、地下水埋深等方面的代表性；有固定的进水口和排水口，配备量水设施 [3] 。对于播种面积超过灌区总播种面积10%以上的作物种类，须分别选择典型田块。根据《细则》的规定，中型灌区样点灌区应至少在上、下游有代表性的农渠控制范围内分别选取，每种需观测的作物种类至少选取3个典型田块。小型灌区样点灌区应按照作物种类、耕作和灌溉制度与方法、田面平整程度等因素选取典型田块，每种需观测的作物种类至少选取2个典型田块。依据以上原则，结合普洱农业分局及种植情况，选取普洱市典型田块范围与数量，五个中型灌区各选取典型田块2块，五个小型灌区各选取典型田块1块，全市共计15块典型田块。
3.3. 监测方案及各参数分析测算
3.3.1. 总体监测方案
灌溉水有效利用系数测算方法多，也越来越准确、细致，但实际应用于灌区的并不多，根据《细则》中的
Table 1. Selection results at sample irrigation area in Pu’er City
表1. 普洱市样点灌区选取结果
推荐，本次普洱市以首尾测算法作为主要测算方法。对毛灌溉用水量的计算采用统计加测量的方法进行。根据不同的灌溉用水取水方式确定不同的监测手段，主要采用水位流量关系曲线法和直接测量法。监测方式以田间固定监测和临时采样监测相结合的方式进行。其中，田间固定监测需要安装气象站及土壤水分监测仪，水稻的灌溉定额监测则通过安装压力式自记水位计，临时采样监测采用速测土壤水分仪进行监测。本项目中开展主要作物灌溉试验的15个试验小区均为典型田块，灌溉试验中灌前和灌后土壤墒情或水位数据可直接作为单次灌溉灌溉定额使用。
3.3.2. 样点灌区灌溉水有效利用系数各参数测算
(1) 典型田块灌溉定额监测
根据《灌溉试验规范》(SL 13-2015)有关规定，在每次灌水前后监测典型田块内不同作物年内相应生育期内计划湿润层的土壤质量含水率或体积含水率(或田间水层变化)，计算该次亩均灌溉定额
，得出该典型田块不同作物种类年亩均灌溉定额
[4] 。
旱作物灌水量，根据典型田块灌溉前后计划湿润层土壤含水率的变化确定某次亩均灌溉定额，计算公式如下：
(1)
式中：
为典型田块某次亩均灌溉定额，m3/亩；
为灌水期内典型田块土壤计划湿润层深度，mm；
为典型田块H土层内土壤干容重，g/cm3；
为水的容重，一般可取1，g/cm3；
为某次灌水前典型田块H土层内土壤质量含水率，%；
为某次灌水后典型田块H土层内土壤质量含水率，% [2] 。
水稻灌水量，根据普洱市实际情况，在本次测算中不考虑水稻的湿润灌溉，只考虑淹水灌溉，但在水稻育秧期，应将育秧期某次灌水的亩均灌溉定额按秧田与试验田的面积比例折算到本田，计入水稻年内生育期亩均灌溉定额，淹水灌溉用水量计算如下：
(2)
式中：
为某次灌水前典型田块田面水深，mm；
为某次灌水后典型田块田面水深，mm [2] 。
(2) 样点灌区灌溉定额总量测算
在计算典型田块灌溉定额的基础上，计算某灌区同区域第i种作物的年灌溉定额，计算公式如下：
(3)
式中：
为样点灌区同片区第i种作物的亩均灌溉定额，m3/亩；
为同片区第i种作物第l个典型田块亩均灌溉定额，m3/亩；
为同片区第i种作物第l个典型田块灌溉面积，亩；N为同片区第i种作物典型田块数量，个。
再根据灌区内不同分区不同作物种类灌溉面积，结合不同作物在不同分区的年亩均灌溉定额，计算得出样点灌区年净灌溉用水总量
，计算公式如下：
(4)
式中：
为样点灌区年净灌溉用水总量，m3；
为样点灌区j个片区内第i种作物亩均灌溉定额，m3/亩；
为样点灌区j个片区内第i种作物灌溉面积，亩；m为样点灌区j个片区内的作物种类，种；n为样点灌区片区数量，个；中型灌区n = 2，小型灌区n = 1 [2] 。
(3) 样点灌区毛灌溉用水量的测算
灌区毛灌溉用水总量
是指灌区全年从水源(一个或多个)取用的用于农田灌溉的总水量，该水量通过实测确定。样点灌区年毛灌溉用水总量的计算公式如下：
(5)
式中：
为样点灌区年毛灌溉用水总量，m3；
为样点灌区第i个水源取水量，m3。n样点灌区水源数量，个 [2] 。
(4) 区域灌溉水有效利用系数推算
以中型灌区2个档次样点灌区灌溉水有效利用系数为基础，采用算术平均法分别计算1~5万亩、5~15万亩灌区的灌溉水有效利用系数；然后将汇总得出的1~5万亩、5~15万亩灌区年毛灌溉用水量加权平均得出普洱市中型灌区的灌溉水有效利用系数；小型灌区灌溉水有效利用系数采用算术平均法得到。因此，在已知各规模与类型灌区灌溉用水有效利用系数和年毛灌溉用水量的情况下，灌溉用水有效利用系数按下式计算：
(6)
式中：
、
分别为全市中、小型灌区的年毛灌溉用水量，万m3；
、
分别为全市中、小型灌区灌溉用水有效利用系数 [2] 。
3.3.3. 主要作物需水量试验设计与方法
(1) 试验目的
开展普洱市主要作物需水量试验，观测主要作物在淹水灌溉模式(农民常规灌溉模式)下全生育期的气象特征、灌溉定额、计划湿润层土壤水分的变化、物候及其他生殖生理指标，分析研究作物全生育期需水总量和净灌溉用水总量、不同生育阶段需水规律和作物系数以及淹水灌溉模式下的灌溉制度，为科学计算主要作物灌溉定额提供依据。
(2) 处理设计
以普洱市种植业区划、农业发展规划为依据，综合考虑全市各县区的地理位置、地形、土壤、气候条件、水资源特征及管理等因素，分别在孟连县设1个水稻、蚕豆轮作灌溉试验点，景谷县设1个玉米、小麦轮作灌溉试验点，思茅区设1个蔬菜灌溉及1个水稻、小麦轮作试验点，景东县设1个水稻、大豆论作试验点，澜沧县设1个水稻、玉米、蚕豆轮作试验点。旱作物分别在江城设置1个咖啡试验点，宁洱设置1个茶叶试验点，墨江设置1个烤烟、小麦轮作试验点，西盟设置1个油菜、陆稻、米荞轮作试验点，几个灌溉试验点水肥管理按当地群众习惯模式进行，不做受旱处理试验。
(3) 试验小区建设
在各个试验点分别建设试验小区开展灌溉试验(见图2)，每个试验点建设3个重复试验小区，每个小区面积60 m2 (5 m × 12 m)。试验小区之间用田埂隔开(水稻试验区田埂需做防渗处理)；每个小区在进水口处安装1个水表，在主管道上安装1个水表，共安装4个水表；在每个小区中部靠近作物根系的位置安装1套土壤水分监测仪，每个水稻小区安装一个自计水位计；在每个试验小区附近开阔地带安装云智能气象站1套。
(4) 观测项目及方法
旱作和水稻均需要进行气象监测同时需要记录作物生育期灌溉用水量、生长发育指标、日常农事活动以及收获时的产量。此外旱作需要记录作物计划湿润层内土壤含水率的变化情况，地下水埋深较浅的地方还需要增加地下水位观测；水稻需要记录稻田内水层深度的变化以及稻田的垂向渗漏。普洱市主要作物灌溉试验所需观测的主要项目及方法如表2所示。
(5) 主要作物需水量测算方法
参考作物需水量计算采用Penman-Monteith公式计算参考作物蒸发蒸腾量(ET0)，其计算公式为 [5] ：
Figure 2. Layout of the main crop water demand test plot
图2. 主要作物需水量试验小区布置图
Table 2. Observation items and test methods
表2. 观测项目与试验方法
(7)
式中：
为参考作物蒸发蒸腾量(mm/d)；
为饱和水汽压与温度曲线的斜率(kPa/℃)；为作物冠层表面的净辐射(MJ/(m2∙d))；为土壤热通量(MJ/(m2∙d))，在逐日估算时取；
为日平均气温(℃)，按最高气温(
)和最低气温(
)的算术平均值计算；为2m高度处的风速(m/s)；为饱和水汽压(kPa)；为实际水汽压(kPa)；为饱和水汽压差，(kPa)；为干湿表常数(kPa/℃)。
旱作物实际作物需水量计算采用水量平衡法计算 [2] ，计算公式为：
(8)
式中：为第i天作物实际需水量(mm)；为第i天作物计划湿润层内土壤储水量的变化值，通过体积含水率变化和计划湿润层深度的乘积求得(mm)；为第i天作物计划湿润内储存的有效降水总量(mm)；为第i天有效灌溉水总量，是灌溉前和灌溉后土壤计划湿润层内储水量变化值(mm)；第i天地下水补给量，在有底测坑或地下水埋深大于2.5 m (砂土、沙壤土)、3.5 m (壤土、粘壤土、粘土)条件下，；为第i天内的排水量，本文主要指土壤计划湿润层在降水或灌溉之后向下渗漏的水量(mm)。
(9)
作物系数计算，作物系数(Kc)是计算作物需水量的重要参数，它反映了作物本身的生物学特性、产量水平、土壤耕作条件对作物需水量的影响，作物系数反映了作物在整个生长季节耗水的变化，在不同的生育阶段有不同的值。作物系数的计算方法在FAO-56中主要分为单作物系数法和双作物系数法，本文采用单作物系数法计算普洱市内各种特色农作物的作物系数 [2] 。Kc的计算公式为：式中：Kc为作物系数(无量纲)；ETc为实际作物需水量(mm/d)；ET0为参考作物需水量(mm/d)。
4. 项目实施与计划
为贯彻落实最严格水资源管理制度，圆满完成普洱市农田灌溉水有效利用系数考核目标任务，根据相关技术规程规范要求，本项目实施工作期为3年(2017~2019年)。即通过对所选样点灌区开展为期三年的灌溉试验及灌溉水有效利用系数监测，形成普洱市主要作物需水量、灌溉定额、作物系数计各类型灌区灌溉用水有效利用系数分析年度成果，在3年测算成果的基础上编制《普洱市农田灌溉水有效利用系数测算及灌溉试验成果报告》。该项目主要由普洱市水务局负责此项工作的组织领导及项目实施过程中的监督管理；样点灌区所属各县(区)水务局及灌区管理局负责为本项目的实施提供相关工作配合。普洱市水文水资源局作为技术支撑单位负责项目具体实施及相关报告的编制。
5. 结语
本文紧扣最严格水资源管理制度考核要求，依据《灌溉水利用率测定技术导则》和《灌溉试验规范》的相关技术规定，经过对全市灌区资料的收集并结合选址查勘，选定样点灌区和典型田块，制定合理可行的农田灌溉水有效利用系数测算和主要作物需水量试验技术方案。主要监测水稻、蚕豆、玉米、小麦、大豆、咖啡、茶叶、烤烟、油菜、陆稻、米荞及各类蔬菜全生育期作物需水量、灌溉用水量及作物系数。建成具有70个左右半固定监测点的普洱市农田灌溉用水效率监测网络，为普洱市农田灌溉水有效利用系数测算及特色粮经作物灌溉用水定额制定提供基础支撑。也填补了普洱市农田灌溉水有效利用系数和主要作物需水量的空白，对推动农业节水管理提供科学依据，提高普洱市农田灌溉用水效率，促进全市高原特色现代化农业发展具有重要意义。