一种基于动态规划方法的自流灌区单干渠优化配水方法技术

本发明专利技术公开了一种基于动态规划方法的自流灌区单干渠优化配水方法，包括以下步骤：(1)以一次灌水期单干渠所辖各支渠控制的受水区缺水量平方和最小为目标，建立如下目标函数；(2)设置约束条件；(3)模型求解：(a)测定单干渠控制的支渠首节制闸位置、数量，测定各支渠控制的各受水区作物在该灌水期内的需水量；(b)以各支渠首节制闸位置编号为阶段，以前若干个支渠首节制闸配水总量为状态变量，构建状态转移方程，确定各阶段效益费用函数；(c)采用顺序递推法，获得满足各支渠控制面积缺水量平方和最小的各支渠首节制闸最优配水量。本发明专利技术能够获取各支渠首节制闸优化配水方法，获得各支渠首节制闸最优的开度过程。首节制闸最优的开度过程。首节制闸最优的开度过程。

【技术实现步骤摘要】
一种基于动态规划方法的自流灌区单干渠优化配水方法


[0001]本专利技术涉及渠道优化配水的
，具体是涉及一种基于动态规划方法的自流灌区单干渠优化配水方法。

技术介绍

[0002]大中型灌区骨干输水渠道分级多，线路长，沿线控制性建筑物分布广，由于引水水权收到控制，加之管理程度粗放不一，灌溉水利用效率往往较低，影响农业增产和农民增收。如何根据灌区农业种植结构，通过骨干输水渠道沿线控制性节制闸优化调度，实现灌区水资源最大化综合利用，对提高灌溉水利用效率，保障灌区粮食生产，促进农民增收。针对传统优化渠系配水模型主要基于目标函数对作物需水量进行优化配置，而目标函数参数及约束条件较为复杂，难以实现全局优化配水。

技术实现思路

[0003]专利技术目的：针对以上缺点，本专利技术提供一种基于动态规划方法的自流灌区单干渠优化配水方法，以单干渠为研究对象，能够获取以受水区作物缺水量最小为目标的支渠首节制闸优化配水方法，获得单干渠沿线各支渠首节制闸最优的开度过程。
[0004]技术方案：为解决上述问题，本专利技术提供一种基于动态规划方法的自流灌区单干渠优化配水方法，包括以下步骤：
[0005](1)模型建立：以一次灌水期单干渠所辖各支渠控制的受水区缺水量平方和最小为目标，建立如下目标函数；
[0006][0007]式中：F为一次灌水期单干渠所辖各支渠控制的受水区缺水量的最小平方和；f为单干渠所辖各支渠控制的受水区缺水量的平方和；m为干渠控制的支渠数量；j为支渠编号，j＝1,2,
…r/>,m；i为一次灌水期划分的灌水时段编号，i＝1,2,
…
, n，n为时段总数；X
j
(K
ji
)为第j条支渠对应于渠首节制闸开度过程K
ji
的配水量，单位为m3；YS
j
为第j条支渠控制面积对应的作物需水量，单位为m3；
[0008](2)设置约束条件：根据单干渠一次引水总量、灌水时长、干渠流量、干渠水位、支渠首节制闸开度的要求，设置目标函数求解的约束条件；
[0009](3)模型求解：以单干渠沿线各支渠首节制闸配水量为决策变量，以前若干个支渠首节制闸配水总量为状态变量，采用一维动态规划方法，结合约束条件对目标函数求解，具体包括以下步骤：
[0010](31)测定单干渠控制的支渠首节制闸位置、数量，测定各支渠控制的各受水区作物在该灌水期内的需水量；
[0011](32)以单干渠沿线各支渠首节制闸位置编号为阶段，以前若干个支渠首节制闸配水总量为状态变量，构建状态转移方程，确定各阶段效益费用函数；
[0012](33)采用顺序递推法，获得满足各支渠控制面积缺水量平方和最小的各支渠首节制闸最优配水量。
[0013]有益效果：本专利技术相对于现有技术，其显著优点是：结合大中型灌区干渠配水实际情况，设置干渠首引水量约束、干渠配水总时长约束、干渠流量约束、干渠水位约束和支渠首闸门开度约束，提高了模型求解结果的实用性。
[0014]进一步的，本专利技术还包括步骤：
[0015](4)采用一维非恒定流模型模拟各支渠首节制闸闸上水位过程；确定各支渠首节制闸闸上水深过程；再运用闸孔自由出流公式，确定各支渠首节制闸闸门最优开度过程，具体包括以下步骤：
[0016](41)测定单干渠计算参数，包括：渠段数，各渠段渠底高程、渠底宽、边坡系数、糙率；测定各支渠首节制闸数闸底板高程、闸孔净宽、闸孔净高；
[0017](42)根据单干渠一次引水总量及灌水总时长，拟定干渠首引水流量过程；根据获得的各支渠首节制闸最优配水量及灌水总时长，确定各支渠首节制闸配水流量过程；
[0018](43)将获得的各支渠首节制闸配水流量过程作为流量边界，在已知干渠首引水流量过程、干渠初始水位，以及测得的单干渠各渠段断面参数、各支渠首节制闸参数的基础上，采用一维非恒定流模型模拟确定各支渠首节制闸闸上水位过程，并以此确定各支渠首节制闸闸上水深过程；再运用闸孔自由出流公式，计算出各支渠首节制闸闸门最优开度过程。
[0019]有益效果：本专利技术相对于现有技术，其显著优点是：将动态规划法与一维非恒定流数值模拟相结合，能够进一步获得各支渠首节制闸闸门在各时间段内最优开度过程，提高模型求解精度，方便对支渠首闸门开展有效实时调度。
[0020]进一步的，所述步骤(2)中的约束条件包括：
[0021](21)干渠首引水量约束：
[0022][0023](22)干渠灌水总时长约束：
[0024]T
s
≤T0[0025](23)干渠流量约束：
[0026]Q
min
≤Q≤Q
max
[0027](24)干渠水位约束：
[0028]Z
min
≤Z≤Z
max
[0029](25)支渠首闸门开度约束：
[0030]0≤K
ji
≤H
j
[0031]式中，W0为干渠首一次引水总量，单位为m3；T
s
为实际灌水总时长，单位为h；灌水总时长T0，单位为h；Q为干渠流量，单位为m3/s；最大流量Q
max
、最小流量Q
min
，单位为m3/s；Z为干渠输水水位，单位为m；最高水位Z
max
、最低水位Z
min
，单位为m；H
j
为第j个支渠首闸门净高，单位为m。
[0032]进一步的，步骤(32)具体包括以下步骤：
[0033](321)以支渠首节制闸编号j为阶段变量，以前j个支渠首节制闸配水总量为状态
变量λ，构建如下状态转移方程：
[0034]λ
j
＝X
j
(K
ji
)+λ
j
‑1[0035](322)根据目标函数及状态转移方程，确定如下效益费用函数：
[0036]g
j
(λ
j
)＝min{[X
j
(K
ji
)
‑
YS
j
]2+g
j
‑1(λ
j
‑1)}
[0037]进一步的，步骤(33)具体包括以下步骤：
[0038](331)阶段j＝1：
[0039]g1(λ1)＝min[X1(K
1i
)
‑
YS1]2[0040]不考虑节制闸开度K
1i
影响，将该阶段状态变量λ1可在对应可行域内离散：λ1＝0,W1,W2,
…
,W0；对每个离散的λ1，决策变量X1可在对应可行域内离散，应满足：X1≥λ1；将满足要求的X1分别代入上式，根据该阶段已知的YS1，分别得到每个离散λ1值时，最优X
1*
及其对应的g1(λ1)；
[本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于动态规划方法的自流灌区单干渠优化配水方法，其特征在于，包括以下步骤：(1)模型建立：以一次灌水期单干渠所辖各支渠控制的受水区缺水量平方和最小为目标，建立如下目标函数；式中：F为一次灌水期单干渠所辖各支渠控制的受水区缺水量的最小平方和；f为单干渠所辖各支渠控制的受水区缺水量的平方和；m为干渠控制的支渠数量；j为支渠编号，j＝1,2,
…
,m；i为一次灌水期划分的灌水时段编号，i＝1,2,
…
,n，n为时段总数；X
j
(K
ji
)为第j条支渠对应于渠首节制闸开度过程K
ji
的配水量，单位为m3；YS
j
为第j条支渠控制面积对应的作物需水量，单位为m3；(2)设置约束条件：根据单干渠一次引水总量、灌水时长、干渠流量、干渠水位、支渠首节制闸开度的要求，设置求解目标函数的约束条件；(3)模型求解：以单干渠沿线各支渠首节制闸配水量为决策变量，以前若干个支渠首节制闸配水总量为状态变量，采用一维动态规划方法，结合约束条件对目标函数求解，具体包括以下步骤：(3.1)测定单干渠控制的支渠首节制闸位置、数量，测定各支渠控制的各受水区作物在该灌水期内的需水量；(3.2)以单干渠沿线各支渠首节制闸位置编号为阶段，以前若干个支渠首节制闸配水总量为状态变量，构建状态转移方程，确定各阶段效益费用函数；(3.3)采用顺序递推法，获得满足各支渠控制面积缺水量平方和最小的各支渠首节制闸最优配水量。2.根据权利要求1所述的自流灌区单干渠优化配水方法，其特征在于，还包括：(4)采用一维非恒定流模型模拟各支渠首节制闸闸上水位过程；确定各支渠首节制闸闸上水深过程；再运用闸孔自由出流公式，确定各支渠首节制闸闸门最优开度过程，具体包括以下步骤：(4.1)测定单干渠计算参数，包括：渠段数，各渠段渠底高程、渠底宽、边坡系数、糙率；测定各支渠首节制闸数量、闸底板高程、闸孔净宽、闸孔净高；(4.2)根据单干渠一次引水总量及灌水总时长，拟定干渠首引水流量过程；根据获得的各支渠首节制闸最优配水量及灌水总时长，确定各支渠首节制闸配水流量过程；(4.3)将获得的各支渠首节制闸配水流量过程作为流量边界，在已知干渠首引水流量过程、干渠初始水位，以及测得的单干渠各渠段断面参数、各支渠首节制闸参数的基础上，采用一维非恒定流模型模拟确定各支渠首节制闸闸上水位过程，并以此确定各支渠首节制闸闸上水深过程；再运用闸孔自由出流公式，计算出各支渠首节制闸闸门最优开度过程。3.根据权利要求1所述的自流灌区单干渠优化配水方法，其特征在于，所述步骤(2)中的约束条件包括：(2.1)干渠首引水量约束：
(2.2)干渠灌水总时长约束：T
s
≤T0(2.3)干渠流量约束：Q
min
≤Q≤Q
max
(2.4)干渠水位约束：Z
min
≤Z≤Z
max
(2.5)支渠首闸门开度约束：0≤K
ji
≤H
j
式中，W0为干渠首一次引水总量，单位为m3；T
s
为实际灌水总时长，单位为h；灌水总时长T0，单位为h；Q为干渠流量，单位为m3/s；最大流量Q
max
、最小流量Q
min
，单位为m3/s；Z为干渠输水水位，单位为m；最高水位Z
max
、最低水位Z
min
，单位为m；H
j
为第j个支渠首闸门净高，单位为m。4.根据权利要求1所述的自流灌区单干渠优化配水方法，其特征在于，步骤(3.2)具体包括以下步骤：(3.2.1)以支渠首节制闸编号j为阶段变量，以前j个支渠首节制闸配水总量为状态变量λ，构建如下状态转移方程：λ
j
＝X
j
(K
ji
)+λ
j
‑1(3.2.2)根据目标函数及状态转移方程，确定如下效益费用函数：g
j
(λ
j
)＝min{[X
j
(K
ji
)
‑
YS
j
]2+g
j
‑1(λ...

【专利技术属性】
技术研发人员：龚懿，车宇阳，程吉林，朱博文，陈兴，邹文昊，
申请(专利权)人：扬州大学，
类型：发明
国别省市：