一种城市水务系统的雨水排放控制方法技术方案

本发明专利技术公开了一种城市水务系统的雨水排放控制方法。目前城市水务系统无法准确地预测降雨量，给城市排水系统带来了很大困难，易引发城市内涝。本发明专利技术方法包括建立被控对象的状态空间模型、静态输出反馈控制器结构、静态输出反馈控制器求解。本发明专利技术基于具有分段常数转移概率的马尔可夫跳变系统，建立城市排水管网雨水排放的数学模型。利用随机分析方法和矩阵奇异值分解方法，设计出相应的静态输出反馈控制器。本发明专利技术中使用的静态输出反馈控制器结构简单，只需要用到外部可测信号，易于实现，同时也能达到较好的控制要求。本发明专利技术为城市排水系统控制提供了有效的方法，有助于减少城市内涝的发生。

A Rainwater Discharge Control Method for Urban Water System

【技术实现步骤摘要】
一种城市水务系统的雨水排放控制方法
本专利技术属于自动化控制
，涉及一种城市水务系统的雨水排放控制方法，具体是一种分段常数转移概率的马儿可夫跳变系统，利用静态输出反馈控制实现对城市水务系统的雨水排放的控制方法，实现城市水务系统强降雨情况下雨水排放的有效调节。
技术介绍
强降雨对现代化城市的排水系统来说是一个严峻挑战，尤其是经常受台风气象灾害影响的城市。例如位于我国东南沿海的福州和温州，常年受到台风天气带来的强降雨影响，市区易发生内涝。据统计，2017年部分沿海省份的主要城市，如广州、福州、杭州、上海及南京的降雨量都超过了1000mm，个别城市如广州甚至超过2000mm。大量雨水短时间内直接排入下水道以及河道排水系统中，给城市排水系统带来了很大困难，极易引发城市内涝，对人们的生活和安全产生了很大影响。参考我国气象部门的标准，并根据降雨量对排水管网中水流量和水位的影响，可将降雨强度划分为四个不同的模态，即无降雨和降雨速率小于0.25mm/h的零星小雨，降雨速率在0.25mm/h和4.0mm/h之间的小雨和中雨，降雨速率在4.0mm/h和16.0mm/h之间的大雨，降雨速率大于16.0mm/h的暴雨。影响降雨量的因素有很多，如温度、季风、大气对流以及地形条件。在这些复杂因素的共同作用下，使得某一地区的降雨过程存在着显著的随机性和不确定性。某地的降雨过程及对排水管网水流状态的影响，可以认为在上述四个不同模态之间进行随机跳变，且模态之间随机跳变符合马儿可夫特性，由此可建立降雨过程的随机跳变模型。此外，在不同季节和温度，降雨的概率会呈现较大的差异，可以将某地的气温进行分段，在同一温度区间内不同降雨模态之间随机跳变的概率保持不变，但是不同降雨模态之间随机跳变的概率随着温度区间改变而变化。因此，可以将某地的气温进行适当分区，利用不同的分段常数转移概率来表示不同温度区间的降雨概率。将某地一天中主要出现的温度段划分为四个温度区间，例如在夏天，分为低于25℃、25℃～30℃、30℃～35℃和高于35℃。在其他季节，也可进行相应的温度区间划分，如冬天，分为低于-5℃、-5℃～0℃、0℃～5℃和高于5℃。按上述方法，针对某地的降雨过程，建立了一个具有分段转移概率的马儿可夫跳变模型。根据四个降雨强度对应的马儿可夫跳变模态，结合四个不同温度区间的切换模态，设计出合适的静态输出反馈控制器以达到对城市排水系统雨水排放的有效控制。
技术实现思路
本专利技术的目的是针对我国目前城市水务系统无法准确地预测降雨量，从而强降雨天气下导致雨水排放难以有效控制，提供一种城市排水管网雨水排放控制的新方法。本专利技术方法包括：(1).建立被控对象的状态空间模型:(1-1).通过数采仪、流量计、水压表采集城市水务系统中排水管道以及河道的水位值、管道排水量、水流速度数据信息，结合实测水务数据建立如下状态方程：x(k)∈R3，表示k时刻排水系统的管道或者河道中水流的状态向量；y(k)∈R1，表示k时刻排水管道系统的测量输出,符号表示n0维的列向量；x(k)＝[x1(k),x2(k),x3(k)]T，x1(k)、x2(k)、x3(k)分别表示k时刻排水管道或河道的水位值、水流速度、排水量，上标T表示矩阵的转置；当x1(k)的值大于允许值时，表示排水管道或河道水位发生溢出；u(k)∈R1，是k时刻的排水系统控制输入量，表示排水系统中进入排水管道的水流量与流出管道的水流量之差；马尔可夫过程{r(k),k≥0}在有限集合S＝{1,2,3,4}中取值，分别表示四个不同降雨量对应的系统模态之间的跳变；Ar(k)∈R3×3、Br(k)∈R3×1和Cr(k)∈R1×3为已知定常矩阵，符号表示n1×n2维的实矩阵；当r(k)＝i(i∈S)时，Ar(k)、Br(k)、Cr(k)和fr(k)(k,x(k))分别简记为Ai、Bi、Ci和fi(k,x(k))；β(k)是伯努利随机序列，表示系统的非线性扰动是随机发生的，其数学期望和方差分别为和β*，即E{}为数学期望的符号；fi(k,x(k))∈R3是随机发生的无污染工业废水和居民生活废水排入河道对水流状态产生的非线性干扰，其满足如下假设：对于任意向量u,v∈R3，非线性函数fi(k,x(k))满足fi(k,0)＝0以及||fi(k,u)-fi(k,v)||≤||Fi(u-v)||；其中，Fi∈R3×3是已知矩阵，||·||表示矩阵或向量的欧几里得范数；(1-2).切换信号σk在有限集合M＝{1,2,3,4}中取值，分别对应四个不同的温度区间；不同降雨条件下，降雨过程的转移概率矩阵表示为：其中转移概率是一个由切换信号σk决定的分段常数函数，即当σk＝m(m∈{1,2,3,4})时，对所有的i,j∈S，都有成立；表示在第m个温度区间的转移概率，Prob(·)表示随机事件的概率；(1-3).利用实测数据和计算机仿真技术，进行模型校验和修正。(2).静态输出反馈控制器结构：利用实际被控对象的测量输出y(k)，构造排水管网系统的静态输出反馈控制律u(k)＝Ki,my(k)，i∈S,m∈M，其中Ki,m为待求的控制器增益矩阵，记将排水系统状态方程转化为闭环系统：(3).静态输出反馈控制器求解：(3-1).构造Lyapunov函数V(r(k),σk)＝xT(k)Pi,mx(k)，Pi,m＞0(i∈S,m∈M)，为正定对称矩阵；根据构造的Lyapunov函数得到：其中，根据非线性假设条件，对于任意正标量εi，得到：等价为：其中I表示维数匹配的单位矩阵，由此得到：其中式中*表示矩阵中的对称项。(3-2).控制器增益矩阵的求解：在闭环系统稳定的基础上，利用线性矩阵不等式和矩阵奇异值分解方法求解控制器增益矩阵Ki,m：首先，根据Schur补引理及矩阵合同变换方法，将Ψi＜0等价转变为：其中，式中diag{}表示对角矩阵；对于列满秩矩阵Bi∈R3×1，通过矩阵奇异值分解算法，可得两个正交矩阵Ui∈R3×3和Vi∈R1×1，使得Bi可以分解为其中Λi为矩阵Bi的非零奇异值，上标-1表示矩阵的逆。因此，存在矩阵Si∈R3×1，使得其中U1i∈R1×3,U2i∈R2×3，并且存在正定对称矩阵P1i,m∈R1×1,P2i,m∈R2×2满足：那么就存在一个非奇异矩阵使得结合上述条件，得到令得到下述线性矩阵不等式：利用Matlab软件的线性矩阵不等式工具箱，求解该线性矩阵不等式可以得到Xi,m的值，最后得到增益矩阵本专利技术方法基于具有分段常数转移概率的马尔可夫跳变系统，对排水管道水流状态进行描述，建立城市排水管网雨水排放的数学模型。利用随机分析方法建立稳定性条件，结合矩阵奇异值分解以及矩阵变换方法，设计出相应的静态输出反馈控制器，进而通过控制泵站网络和排水管道闸门来优化城市排水系统的调度和控制。相比于状态反馈控制器，本专利技术选择的静态输出反馈控制器结构简单，只需要用到外部可测信号，因而易于实现，同时也能达到较好的控制要求。本专利技术为城市排水系统控制提供了有效的方法，减少了城市内涝发生。具体实施方式一种城市水务系统的雨水排放控制方法，包括：(1).建立被控对象的状态空间模型:(1-1).通过数采仪、流量计、水压表采集城市水务系统中排水管道以及河道的水位值、管道排水量、水流本文档来自技高网...

【技术保护点】
1.一种城市水务系统的雨水排放控制方法，其特征在于该方法包括：(1).建立被控对象的状态空间模型:(1‑1).通过数采仪、流量计、水压表采集城市水务系统中排水管道以及河道的水位值、管道排水量、水流速度数据信息，结合实测水务数据建立如下状态方程：

【技术特征摘要】
1.一种城市水务系统的雨水排放控制方法，其特征在于该方法包括：(1).建立被控对象的状态空间模型:(1-1).通过数采仪、流量计、水压表采集城市水务系统中排水管道以及河道的水位值、管道排水量、水流速度数据信息，结合实测水务数据建立如下状态方程：x(k)∈R3，表示k时刻排水系统的管道或者河道中水流的状态向量；y(k)∈R1，表示k时刻排水管道系统的测量输出,符号表示n0维的列向量；x(k)＝[x1(k),x2(k),x3(k)]T，x1(k)、x2(k)、x3(k)分别表示k时刻排水管道或河道的水位值、水流速度、排水量，上标T表示矩阵的转置；当x1(k)的值大于允许值时，表示排水管道或河道水位发生溢出；u(k)∈R1，表示k时刻的排水系统控制输入量，表示排水系统中进入排水管道的水流量与流出管道的水流量之差；马尔可夫过程{r(k),k≥0}在有限集合S＝{1,2,3,4}中取值，分别表示四个不同降雨量对应的系统模态之间的跳变；Ar(k)∈R3×3、Br(k)∈R3×1和Cr(k)∈R1×3为已知定常矩阵，符号表示n1×n2维的实矩阵；当r(k)＝i，i∈S时，Ar(k)、Br(k)、Cr(k)和fr(k)(k,x(k))分别简记为Ai、Bi、Ci和fi(k,x(k))；β(k)是伯努利随机序列，表示系统的非线性扰动是随机发生的，其数学期望和方差分别为和β*，即E{}为数学期望的符号；fi(k,x(k))∈R3是随机发生的无污染工业废水和居民生活废水排入河道对水流状态产生的非线性干扰，其满足如下假设：对于任意向量u,v∈R3，非线性函数fi(k,x(k))满足fi(k,0)＝0以及||fi(k,u)-fi(k,v)||≤||Fi(u-v)||；其中，Fi∈R3×3是已知矩阵，||·||表示矩阵或向量的欧几里得范数；(1-2).切换信号σk在有限集合M＝{1,2,3,4}中取值，分别对应四...

【专利技术属性】
技术研发人员：陈云，陈臻誉，陈张平，薛安克，
申请(专利权)人：杭州电子科技大学，
类型：发明
国别省市：浙江,33