一种基于分布估计算法的灌溉施肥管网优化方法与流程

本发明涉及农业和园林灌溉技术领域，尤其涉及一种基于分布估计算法的灌溉施肥管网优化方法。

背景技术：

在农业或园林管网设计时，主要考虑首部运行和管网设计的经济性。灌溉的设计一般为树状结构的多段管网组成，在考虑施肥的水肥一体化的灌溉管网设计中，费用与管网投资、动力费用及取水口数相关。管径需求越大，管网投资也随之增加，而对水肥一体机的动力要求越小；如果如取水口一定，缩小管径以减小管网费用时，扬程需增加，动力费增加。因此，在设计考虑施肥的水肥一体系统中，需平衡管径和投资费用之间的关系，设计经济实用的施肥一体化管网系统，节约工程投资运行成本。

在灌溉管网的设计中，考虑水肥一体化模型是一个涉及多种约束条件的复杂优化问题。许多学者在管网优化设计中进行了相关的研究，将遗传算法、模拟退火算法、及人工神经网络算法应用求解该优化问题，取得了一些成果，但仍然存在容易陷入局部最优、没有考虑设备施肥的局限，所求得方案经济效益并不是最优解。

技术实现要素：

本发明所要解决的技术问题是针对现有求解灌溉施肥管网优化中问题，提供一种基于分布估计算法的灌溉施肥管网优化方法，采用岛群算法和分布估计算法相结合的方法，处理考虑施肥的水肥一体化问题，改进管网优化水平，提高管网运行效益。

本发明解决上述技术问题的技术方案如下：

一种基于分布估计算法的灌溉施肥管网优化方法，包括以下步骤：

步骤1：建立灌溉施肥管网优化优化数学模型；

步骤2：对所有的管网进行分配，使得所述步骤1中建立的数学模型的目标函数最小，即使得灌溉施肥系统在满足种植需求时运行费用最低；

所述步骤1中建立的数学模型为，

其中：li,j为第i种轮灌组时第j段管网第k个标准管径的长度；ci,j为第i种轮灌组时第j段管网第k个标准管径的长度；a(m)满足压力要求的第m种型号的水肥机价格；fi为第i种轮灌组管网年费用；x为折旧年限；r为年利率；b为年平均维修费率；e为电价；t为水泵年工作小时数；q为水泵流量；h为水泵扬程；η为水泵效率。

约束条件：

1)管道压力约束：

其中：hcg为管段计算最小压力；hmg为计算管段承受压力；h为管网首端压力；hb为底阀及吸水管水头损失；a为局部水头放大损失；f、m、b为管材相关水头损失系数；qj为第j段流量；dij为第i种轮灌组情况下第j段管径内径；zg为g节点地面高程；hgmin，hgmax为g管段最小、最大压力约束；j(g)为g节点父节点数。

2)管径约束：1≤di,j(k)≤m式中m为标准管径数；

3)管长约束：式中lj为第j段管长度；

4)非负约束：si,j(k)≥0

5)施肥机压力约束：h(m)>h式中h(m)为型号为m的施肥一体机的扬程

进一步：所述步骤2采用岛群算法和分布估计算法相结合的方法进行求解，其具体步骤为，

步骤2.1：初始化分布估计算法参数：设置算法参数，在满足变量非负约束、管网长度约束及管径约束的前提下，采用随机方式生成m个岛群中个体数量，每个岛中有n个个体。

步骤2.2：计算个体的适应度值，分析计算管道压力约束是否满足，并采用罚函数法对个体的适应度值进行修正；根据修正后的适应度值排序。

步骤2.3：随机岛群的顺序{s1,s2,s3,s4,…sm},按生成的顺序进行移民操作，将一个岛屿中随机k个个体移动到下一个岛屿中，如将第s1岛屿中的k个个体移动到s2中，顺序的最后一个sm岛屿中的k个个体移动到s1中。

步骤2.4.分析个岛屿中群体空间中个体，计算得到均值参数和方差参数，采用高斯分布进行采样，生成下一代群体需满足管网长度约束。

步骤2.5.采用随机局部搜索方法，并结合自适应搜索范围参数，动态的调整每个岛屿的局部搜索范围，调整的规律与群体的进化代数相关。在局部搜索时，应考虑管网长度约束和非负约束，一次局部搜索调整的变量为成对出现的两个变量，使其在满足管网长度等式约束的前提下进行局部搜索。每代进行局部搜索20次。

步骤2.6.循环到步骤2.2，直到当进化l代后算法停止。

局部搜索方法具体为，采用爬山搜索算法，在个体的每个维度进行爬山收索，即初始个体为pi时，计算pi+ri或pi+ri时个体适应度值是否改进，如果改进则替换原有个体，其中，ri为局部搜索算子，其值为ri＝r0*e^-αgen/gmax其中r0为初始范围，通常为变量定义域范围，gen为当前进化代数，gmax为最大进行代数，α为调节系统。

本发明的有益效果是：通过建立考虑施肥设备的管网优化数学模型，将施肥需求考虑到管网设计中，更符合现代管网设计的要求；采用岛群算法与分布估计算法相结合的方法，提高求解管网优化问题解的质量，提高管网设计经济效益。

附图说明

图1为本发明的方法流程图。

图2为本发明的实施例中的管网示意图。

具体实施方式

以下结合附图对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。

如图1所示，本发明包括以下步骤：

步骤1：建立灌溉施肥管网优化优化数学模型；

其中：li,j为第i种轮灌组时第j段管网第k个标准管径的长度；ci,j为第i种轮灌组时第j段管网第k个标准管径的长度；a(m)满足压力要求的第m种型号的水肥机价格；fi为第i种轮灌组管网年费用；x为折旧年限；r为年利率；b为年平均维修费率；e为电价；t为水泵年工作小时数；q为水泵流量；h为水泵扬程；η为水泵效率。

约束条件：

1)管道压力约束：

其中：hcg为管段计算最小压力；hmg为计算管段承受压力；h为管网首端压力；hb为底阀及吸水管水头损失；a为局部水头放大损失；f、m、b为管材相关水头损失系数；qj为第j段流量；dij为第i种轮灌组情况下第j段管径内径；zg为g节点地面高程；hgmin，hgmax为g管段最小、最大压力约束；j(g)为g节点父节点数。

2)管径约束：1≤di,j(k)≤m式中m为标准管径数；

3)管长约束：式中lj为第j段管长度；

4)非负约束：si,j(k)≥0

5)施肥机压力约束：h(m)>h式中h(m)为型号为m的施肥一体机的扬程

步骤2：对所有的管网进行分配，使得所述步骤1中建立的数学模型的目标函数最小，即使得灌溉施肥系统在满足种植需求时运行费用最低；具体是采用岛群算法和分布估计算法相结合的方法进行求解，其具体步骤为，

步骤2.1：初始化分布估计算法参数：设置算法参数，在满足变量非负约束、管网长度约束及管径约束的前提下，采用随机方式生成m个岛群中个体数量，每个岛中有n个个体。

步骤2.2：计算个体的适应度值，分析计算管道压力约束是否满足，并采用罚函数法对个体的适应度值进行修正；根据修正后的适应度值排序。

步骤2.3：随机岛群的顺序{s1,s2,s3,s4,…sm},按生成的顺序进行移民操作，将一个岛屿中随机k个个体移动到下一个岛屿中，如将第s1岛屿中的k个个体移动到s2中，顺序的最后一个sm岛屿中的k个个体移动到s1中。

步骤2.4.分析个岛屿中群体空间中个体，计算得到均值参数和方差参数，采用高斯分布进行采样，生成下一代群体需满足管网长度约束。

步骤2.5.采用随机局部搜索方法，并结合自适应搜索范围参数，动态的调整每个岛屿的局部搜索范围，调整的规律与群体的进化代数相关。在局部搜索时，应考虑管网长度约束和非负约束，一次局部搜索调整的变量为成对出现的两个变量，使其在满足管网长度等式约束的前提下进行局部搜索。每代进行局部搜索20次。

步骤2.6.循环到步骤2.2，直到当进化l代后算法停止。

步骤2.5.1局部搜索方法具体为，采用爬山搜索算法，在个体的每个维度进行爬山收索，即初始个体为pi时，计算pi+ri或pi+ri时个体适应度值是否改进，如果改进则替换原有个体。

步骤2.5.2ri为局部搜索算子，其值为ri＝r0*e^-αgen/gmax其中r0为初始范围，通常为变量定义域范围，gen为当前进化代数，gmax为最大进行代数，α为调节系统。

本发明的有益效果是：通过建立考虑施肥设备的管网优化数学模型，将施肥需求考虑到管网设计中，更符合现代管网设计的要求；采用岛群算法与分布估计算法相结合的方法，提高求解管网优化问题解的质量，提高管网设计经济效益。

具体实施方式，如图2所示的管网，

机井管网管径及对应价格如表1，节点对应海拔及相应段流量如表2。

表1机井管网不同规格管道价格表

表2机井管网节点高程及流量表

表3水肥机扬程与费用对照表

管道承压水头为40m，节点出水水头为12m，水源地高程为30m，n为0.013，x为20，e为0.2，η为0.7，r为8％，b为5％，hb为0.2，a为1.05，f为1.312×10⁶，m为2，b为5.33。沿程水头损失取10％，k＝1.528×104，c＝130，α＝1.1m。岛群分布估计算法群体规模取100，岛群规模取5，种群进化到1000代时停止进化。计算结果如表4所示，局部搜索系统α为4。

表4机井管网优化结果表

结果表明采用本发明具有较好的结果，以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。