一种多水源灌区灌溉用水量复核方法与流程

本发明属于水利工程灌溉领域，具体涉及到一种多水源灌区农业灌溉用水量复核方法，属于e02b水利工程分类号下。
背景技术：
：按《实行最严格水资源管理制度考核办法》要求，2014年起，水利部组织最严格水资源管理制度考核工作，由流域机构对各省年度用水总量统计成果复核，占总水量一半以上的农业灌溉用水总量复核是工作的重要部分。探讨农业灌溉用水总量复核方法，对促进最严格水资源管理落实具有重要意义。根据《关于印发用水总量统计方案的通知》（办资源【2014】57号）要求，各省区农业灌溉用水量采取取用水大户逐一计量统计、一般用水户抽样调查、综合推算区域灌溉用水的技术方法。到2015年，设计灌溉面积5万亩以上（含5万亩）的灌区、公共供水户和自备水源年取用水量50万吨以上（含50万吨）的工业企业作为重点统计对象，重点统计对象取用水量占用水总量的比例达到55%以上。到2020年设计灌溉面积1万亩以上（含1万亩）的灌区也要纳入重点统计对象，比例进一步提高到65%。由于我国南方灌区以中小灌区居多，灌区位置分散，水源众多，灌溉用水重复性利用现象显著，计量监控设施建设运行成本巨大且难以监测到真实的灌溉用水量，各省区尚未实现《用水总量统计方案》要求的农业用水统计方法。现有灌溉用水量统计复核方法中，南方多水源灌区灌溉用水量统计方法主要有定额法，灌溉水有效利用系数法，趋势分析法，对比法和经验判断法等方法，不同方法不同口径得到的农业灌溉用水量差别较大。目前的复核方法大多受到限制，复核主要以样点校核法和纵向对比法为主。样点校核法依据区域内样点灌区，计算亩均综合灌溉用水量，用该值乘以相应实际灌溉面积，推算区域灌溉用水量，与上报灌溉水量比较。纵向对比法，结合降水量变化，对比当年与上一年亩均毛灌溉用水量:当降水量增加，亩均灌溉水量一般会减少；当降水量减少，亩均灌溉水量一般会增加。对于南方丰水地区，降水量与亩均灌溉用水量的相关关系不明显，该方法对南方多水源灌溉复核适用性低。而样点校核法作为主要复核方法，目前只能做到定性复核。针对现有技术的缺陷，急需适合南方多水源灌区特性，且简便有效的农业灌溉用水量复核方法。目前南方多水源农业灌溉用水量统计主要采用样点灌区灌溉水有效利用系数推算得出，对于设区市以上尺度样点灌区数量多，代表性好，灌溉用水推算量与实际量接近，而对县（市、区）尺度样点灌区数量少且代表性不足，灌溉用水尺度效应明显，目前尚无有效的南方多水源地区县（市、区）灌溉用水量复核方法。技术实现要素：针对现有技术的不足，本发明提供一种南方多水源灌区灌溉用水量复核方法，构建农业灌溉用水量影响指标体系，通过指标体系将设区市灌溉用水量分配到县（市、区），提出了一种基于可变模糊聚类模型的南方多水源地区农业灌溉用水量复核方法，从而进行县（市、区）灌溉用水量的复核工作。本发明的技术方案为。一种多水源灌区农业灌溉用水量复核方法，具体包括以下步骤：第一步，构建南方多水源灌区农业灌溉用水量影响因素指标体系；灌溉用水量的大小与当地的气候、农作物种类、土地利用系数、复种指数、耕作栽培技术、灌溉制度和灌水技术等诸多因素有关，各区域自然条件、农业生产及经济水平不同，其灌溉用水量有较大差异。首先，将直接影响灌溉用水量的年降雨量，播种面积纳入指标体系。其次，由于不同农作物的灌溉用水量相差较大，根据灌溉定额的不同将各设区市统计年鉴中各县（市、区）农作物播种面积指标分为水田播种面积、水浇地播种面积和菜田播种面积三个指标。其中水田指种植早稻、晚稻和单双季稻的田块，菜田指种植水果、蔬菜和花卉苗木的田块，水浇地指种植其他旱作物的田块。最后，由于不同区域灌溉制度和技术的不同，灌溉用水量还与一些人为因素有关，因此将节水灌溉面积和上年度灌溉用水量作为参考指标纳入指标体系。第二步，基于可变模糊理论，构建可变模糊聚类模型，根据各县（市、区）农业灌溉用水量影响因素指标体系，计算各县（市、区）对设区市农业灌溉用水量的综合贡献度，确定各县（市、区）农业灌溉用水量。由于各县（市、区）影响农业灌溉用水量有多个因素，采用相对隶属度法计算各县对农业用水量的综合贡献度，然后进行归一化，根据归一化后的各县（市、区）对农业用水量的综合贡献度乘以全市农业用水量即为该县（市、区）农业用水量复核值。——县（市、区）农业灌溉用水量指标体系相对隶属度，为设区市灌溉用水量，为县（市、区）灌溉用水量。相对隶属度计算过程如下。基础数据无量纲化。假设待分解灌溉用水量的设区市有n个县（市、区），每个县农业用水量受m个指标影响：m个指标对n个县（市、区）的农业用水量预测可用指标特征值矩阵表示为：xij为样本j指标i的特征值，i=1,2,…,m；j=1,2,…,n。作为多指标因素分析法，指标数据的无量纲化处理是必不可少的，无量纲化的方法有极值化方法、标准化方法、均值化方法、标准差化方法，拟采用极值化方法。经过对指标无量纲化处理后，使其对影响农业灌溉用水量的作用效果具有可比性，根据每个指标对农业灌溉用水量影响趋势进行分类。数量值越大，农业灌溉用水量越大的指标无量纲化数值=各县指标值/各县（市、区）最大值；。数量值越大，农业灌溉用水量越小的指标无量纲化数值=1-各县指标值/各县（市、区）最大值。。指标特征值矩阵变换为指标对农业用水量贡献大小的相对隶属度矩阵，即指标特征值规格化矩阵。rij为第j个县（市、区）对于农业用水量指标i的相对隶属度，且0≤rij≤1。计算综合隶属度。样本集的m个指标对优化的影响程度不同，故指标应具有不同的权重，设指标权向量为满足权重应结合当地实际情况，由专家分析和优化方法率定得到。根据各单项的无量纲化数值和指标权重，可以根据层次分析法或者可变模糊优选模型计算各县（市、区）对农业用水量的综合贡献度，拟采用可变模糊优选模型。应用两级模糊优选模型，得到第j个县各指标的对农业综合贡献相对隶属度——相对隶属度，——模型优化准则参数，=1、2分别为最小一、二乘方准则，p——距离参数，通常可取为海明距离p=1，欧氏距离p=2。本发明的有益效果为：该农业灌溉用水量复核方法充分考虑了影响灌溉用水量的主客观信息，采用可变模糊聚类模型计算了不同县（市、区）农业灌溉对设区市农业灌溉用水量的隶属度（综合贡献度），得到了各县（市、区）农业灌溉用水量，为农业灌溉用水量复核提供了定量参考依据，有效提高农业灌溉用水量统计数据合理性和准确性，更好的为最严格水资源管理考核提供支撑。附图说明附图1为2015年杭州市农业灌溉用水量复核模型率定结果。附图2为2016年杭州市农业灌溉用水量复核模型验证结果。具体实施方式本发明考虑了影响灌溉用水量的主客观信息，构建农业灌溉用水量影响指标体系，提出了一种基于可变模糊聚类模型的南方多水源地区农业灌溉用水量复核方法。下面通过浙江省的实施例，对本发明具体实施方法做进一步说明，凡基于本发明的思想做出的任何改进均在本申请保护范围之内。特定的地形地貌决定了浙江省灌区的结构及其特点——多水源型灌区。其中中西部山丘区及金衢盆地多以水库、山塘、堰坝等为水源，采用自流灌溉方式；而沿江、沿湖滩地和沿海平原地区，如杭嘉湖、萧绍甬、温黄、温瑞等平原河网型灌区多提水灌溉方式。全省灌区分布点多面广，水源工程类型主要有水库、塘坝、引水堰闸、小型泵站等，灌溉方式既有蓄水、引水为主，又有提水灌溉，共同发挥作用，确保灌区的用水需求。以浙江省杭州市为实例进行农业灌溉用水量复核，具体步骤如下。第一步，构建灌区农业灌溉用水量影响因素指标体系。考虑县（市、区）级层面数据的可获取性、相关性等，选取对应指标建立农业灌溉用水量分解指标体系，对各设区市统计年鉴中各县（市、区）农作物播种面积进行分类，结合《浙江省农业用水定额标准》（db33/t769-2009），将早稻、晚稻和单季稻灌溉定额最大，其播种面积作为第一个指标，水果、蔬菜和花卉苗木灌溉定额次之，其播种面积作为第二个指标，其他作物均为旱作物，其播种面积作为第三个指标。根据不同设区市特点，综合考虑其用水效率和管理水平，选取以下指标体系。如表1所示。表1农业用水量分解简单型指标体系。【水田播种面积】指当年水稻实际播种面积之和；【水浇地播种面积】指当年其他旱作物实际播种面积之和；【菜田播种面积】指当年蔬菜水果和花卉苗木实际播种面积之和；【年降水量】指由气象或水文部门记录的灌区所在地区的当年实际降水量；【节水灌溉面积率】指当年实施节水改造的灌区面积与灌区总面积比值；【上年度灌溉用水量】指县市区上年度农业灌溉用水量。第二步，率定农业灌溉用水量影响指标体系权重。选取杭州2015年灌溉用水量影响指标体系数据用于模型指标权重率定，指标体系如表2所示，归一化值如表3所示，率定结果如图1。表2杭州市2015年灌溉用水量分解指标体系值。县市区/指标体系水田播种面积（千公顷）水浇地播种面积（千公顷）菜田播种面积（千公顷）降雨量（mm）节水灌溉面积率（%）上年度灌溉用水量杭州市区674379654219330.770.1467萧山8952198434636119330.773.1846余杭1038772123127419330.771.5544富阳市10408149741871420880.391.3699临安市367971791480220590.331.1316桐庐县434999151070423080.500.9678建德市6706130451344621280.860.9423淳安县1619261411067024361.000.5004表3杭州市2015年灌溉用水量分解指标体系归一化值。县市区/指标体系水田播种面积（千公顷）水浇地播种面积（千公顷）菜田播种面积（千公顷）降雨量（mm）节水灌溉面积率（%）上年度灌溉用水量杭州市区0.064760.014500.141111.000000.428570.04606萧山0.860110.759081.000001.000000.428571.00000余杭0.997980.275890.674581.000000.428570.48810富阳市1.000000.572820.403660.925430.846150.43017临安市0.353480.274630.319280.938441.000000.35534桐庐县0.417850.379290.230880.837350.660000.30390建德市0.644310.499020.290030.908180.383720.29590淳安县0.155551.000000.230160.793380.330000.15713权重率定值0.11340.06920.00720.01380.08930.7071第三步，验证农业灌溉用水量复核模型。选取杭州2016年灌溉用水量影响指标体系数据用于模型验证，指标体系如表4所示，归一化值如表5所示，验证结果如图2所示，结果表明本发明能够克服南方多水源灌区农业灌溉用水复杂性带来的用水量复核困难，不但减少复核工作量，也能达到较好的复核效果。表4杭州市2016年灌溉用水量分解指标体系值。县市区/指标体系水田播种面积（千公顷）水浇地播种面积（千公顷）菜田播种面积（千公顷）降雨量（mm）节水灌溉面积率（%）上年度灌溉用水量杭州市区674379656217770.780.1403萧山8952198434635117770.773.0463余杭1048772323125417770.771.4869富阳市10418149741871417810.391.4387临安市367971891480220830.341.0563桐庐县434999151070419380.500.9500建德市6706130651344618640.860.9522淳安县1609260411067022521.000.4350表5杭州市2016年灌溉用水量分解指标体系归一化值。县市区/指标体系水田播种面积（千公顷）水浇地播种面积（千公顷）菜田播种面积（千公顷）降雨量（mm）节水灌溉面积率（%）上年度灌溉用水量杭州市区0.064270.014550.141571.000000.435900.04606萧山0.853630.761991.000001.000000.441561.00000余杭1.000000.277720.674291.000000.441560.48810富阳市0.993420.575020.403750.997870.871790.47228临安市0.350820.276060.319350.853191.000000.34675桐庐县0.414700.380750.230930.917030.680000.31185建德市0.639460.501710.290090.953430.395350.31258淳安县0.153431.000000.230210.789170.340000.14280权重率定值0.11340.06920.00720.01380.08930.7071当前第1页12