一种基于ARX的供水管网在线建模方法技术

本发明专利技术公开了一种基于ARX的供水管网在线建模方法。本发明专利技术首先建立管网系统实时运行数据，建立系统的输入输出过程模型。实时运行数据库包括系统的实时输入数据和实时输出数据；输入数据包括水厂供水压力、水厂供水量、管网中增压泵站压力和流量、管道上可控阀门的开启程度；输出数据包括压力监测点的压力、流量检测点的流量、蓄水池水位。然后以供水管网实时过程运行数据集合为基础，建立离散差分方程形式的ARX模型。最后实时估计模型参数。本发明专利技术可以在线实时估计模型参数，实时跟踪供水管网的状态变化并做出预测；该模型结构简单、鲁棒性强，结合预测控制和自适应控制理论，实现供水管网实时优化调度。

【技术实现步骤摘要】
一种基于ARX的供水管网在线建模方法
本专利技术属于自动控制
，具体是一种基于ARX(AutoRegressiveExogenous)的供水管网在线建模方法。
技术介绍
供水管网系统是一个结构复杂、规模庞大、用水随机性强的动态网络系统，能否快速准确地模拟和预测出管网的运行工况，是进行供水管网优化调度的关键。供水管网系统模型大致可分为两类：微观模型和宏观模型。微观模型为机理模型，需要已知管网的拓扑结构、管长、管径、管材及节点用水量等具体参数，一般为由质量平衡方程和能量平衡方程组成的非线性系统方程，其忽略了管网系统的动力学特性，不能真实地描述实际管网的状态变化过程。宏观模型是经验模型，如回归分析模型，神经网络模型，支持向量机模型等。该类模型不需要已知管网的基本参数，其从系统的角度出发，抓住管网运行的几个宏观变量，如各水厂供水压力、供水量以及监测点的压力等，从大量实际运行数据中提炼出宏观变量之间经验函数关系，建立起供水管网静态输入输出模型。而实际运行中供水管网系统，由于用户用水量时刻变化，管道中水体的惯性及伸缩性作用等，从水厂或泵站的压力和流量调节到管网各节点的响应是一个复杂的动态过程。此外，一般基于传统微观模型及宏观模型的优化调度为离线调度，而离线调度方案实施时，要根据实际管网的状态变化对控制方案进行实时调整，这需要对管网状态的实时跟踪和预测。显然，传统的模型都难以胜任。要实时反映供水管网状态变化，描述其复杂的动态过程，需要建立合适的供水管网的在线预测模型。ARX模型是一种“黑箱”模型，其利用系统的输入、输出信息建立数学模型，结构简单、鲁棒性强，广泛应用于非线性动态系统。且其模型参数较少，容易实现在线建模，非常适用于时变系统。实际运行中的管网系统，管网中节点实时用水量为时变的且难以预测，故其也是一个时变系统。基于ARX供水管网在线预测模型，可以在线实时估计模型参数，实时跟踪供水管网的状态变化并做出预测。
技术实现思路
本专利技术针对现有技术的不足，提供了一种基于ARX的供水管网在线预测模型。具体建模过程如下：(1)建立管网系统实时运行数据建立系统的输入输出过程模型。具体方法是：1)建立供水管网系统实时运行数据库。实时运行数据库包括系统的实时输入数据和实时输出数据。系统的实时输入数据：水厂供水压力，水厂供水量，管网中增压泵站压力和流量，管道上可控阀门的开启程度等。系统实时输出数据：压力监测点的压力，流量检测点的流量，蓄水池水位等。通过SCADA(SupervisoryControlandDataAcquisition)系统采集N组实时过程运行数据，将采集的实时过程运行数据作为数据驱动的样本集合，表示为j＝1,2,…,N，Xj表示第j组系统输入数据，Yj表示第j组系统输出数据。设供水管网系统有p个输入控制变量，r个输出变量。2)以供水管网实时过程运行数据集合为基础，建立离散差分方程形式的ARX模型：Y(k)＝θTX(k)+E(k)其中：X(k)＝[YT(k-1),…,YT(k-n),UT(k-d-1),…,UT(k-d-m)]T，Y(k)＝[Y1(k),Y2(k),…,Yr(k)]T，U(k)＝[U1(k),U2(k),…,Up(k)]T，其中，Y(k)表示当前时刻供水管网系统输出；X表示系统过去时刻的输入和输出数据的集合；U(k)表示当前供水系统的输入控制变量；E(k)表示当前干扰造噪声，在参数估计中又定义为估计误差；k为当前递推步数；θ表示通过辨识得到的模型参数集合；T表示矩阵的转置；n，m，d+1分别为供水管网系统输出变量的阶次、输入变量的阶次、系统的时滞。(2)实时估计模型参数。具体方法为：1)对于具有p个输入，r个输出的ARX模型，对应于模型第i个输出定义均方估计准则函数：2)依据1)中准则函数，估计θi。估计方法为：θi(k)＝θi(k-1)+mi(k)Pi(k-1)X(k-1)εi(k)εi(k)＝Y(k)-θiT(k-1)X其中，Pi(k)，Fi(k)，Hi(k)分别对应于模型第i个输出的增益矩阵、最小增益矩阵及中间矩阵。且Pi(0)＝αI，Fi(k)＝βI，α,β＞0，I为单位矩阵。函数m(k)：其中，δ(k)表示k时刻，系统输入输出给模型估计过程带来的信息量；δ0为信息量阀值。设置δ0可提高估计过程的鲁棒性，且δ0＞0。λi(k)为自适应可变遗忘因子：其中，ρi参数用于控制遗忘因子对误差εi(k)的敏感度。系统通过λi(k)的自适应调整，动态估计模型参数，适应系统的实时变化，且λi(k)∈[λi,min1]，λi,min(k)∈(01]。3)根据θi，构造参数集合本专利技术的有益效果：基于ARX供水管网在线预测模型，可以在线实时估计模型参数，实时跟踪供水管网的状态变化并做出预测；在线建模的一个难点在于网中节点实时用水量难以预测，而本模型将实时用水量变化作为扰动变量，通过实时估计模型参数来适应其变化，从更新模型参数角度来解决该问题；该模型结构简单、鲁棒性强，结合预测控制和自适应控制理论，容易建立供水管网预测控制和自适应控制模型，实现供水管网实时优化调度。附图说明图1实施例供水管网拓扑示意图。具体实施方式为使本专利技术实现的技术手段与创作特征易于明白，下面结合附图和实施例，对本专利技术的实施方式作进一步详述。本实施例考虑一个有两个水厂的供水管网，如图1。建立以2个水厂供水压力，供水量输入控制变量，以6个压力监测点压力为输出的供水管网ARX模型，模型将管网各个节点用水量看作不可测扰动变量。故该供水管网系统有p＝4个输入控制变量，r＝6个输出变量。(1)根据管网系统实时运行数据建立系统的输入输出过程模型。具体方法为：1)建立该供水管网系统实时运行数据库。实时运行数据库包括系统的实时输入数据和实时输出数据。系统的实时输入数据：水厂供水压力，水厂供水量。系统实时输出数据：6个压力监测点的压力。通过SCADA系统采集N组实时过程运行数据，将采集的实时过程运行数据作为数据驱动的样本集合，表示为Xj表示第j组系统输入数据，Yj表示第j组系统输出数据。2)以供水管网实时过程运行数据集合为基础，建立离散差分方程形式的ARX模型：Y(k)＝θTX(k)+E(k)其中：X(k)＝[YT(k-1),…,YT(k-n),UT(k-d-1),…,UT(k-d-m)]T，Y(k)＝[Y1(k),Y2(k),…,Yr(k)]T，U(k)＝[U1(k),U2(k),…,Up(k)]T，其中，Y(k)表示当前时刻供水管网系统输出；X表示系统过去时刻的输入和输出数据的集合；U(k)表示当前供水系统的输入控制变量；E(k)表示当前干扰造噪声，在参数估计中又定义为估计误差；k为当前递推步数；θ表示通过辨识得到的模型参数集合；T表示矩阵的转置；n，m，d+1分别为供水管网系统输出变量的阶次、输入变量的阶次、系统的时滞，其值分别设置为n＝2，m＝2，d＝0。(2)动态估计模型参数。具体方法为：1)对于具有p个输入，r个输出的ARX模型，对应于模型第i个输出定义均方估计准则函数：2)依据1)中准则函数，估计θi。估计方法为：θi(k)＝θi(k-1)+mi(k)Pi(k-1)X(k-1)εi(k)εi(k)＝Y(k)-θiT(k-1)X其中，P本文档来自技高网...

【技术保护点】
一种基于ARX的供水管网在线建模方法，其特征在于该方法包括以下步骤：步骤(1)、建立管网系统实时运行数据，建立系统的输入输出过程模型；具体是：1)建立供水管网系统实时运行数据库：实时运行数据库包括系统的实时输入数据和实时输出数据；系统的实时输入数据包括水厂供水压力、水厂供水量、管网中增压泵站压力和流量、管道上可控阀门的开启程度；系统实时输出数据包括压力监测点的压力、流量检测点的流量、蓄水池水位；通过SCADA系统采集N组实时过程运行数据，将采集的实时过程运行数据作为数据驱动的样本集合，表示为j＝1,2,…,N，Xj表示第j组系统输入数据，Yj表示第j组系统输出数据；设供水管网系统有p个输入控制变量，r个输出变量；2)以供水管网实时过程运行数据集合为基础，建立离散差分方程形式的ARX模型：Y(k)＝θTX(k)+E(k)其中：X(k)＝[YT(k‑1),…,YT(k‑n),UT(k‑d‑1),…,UT(k‑d‑m)]T，Y(k)＝[Y1(k),Y2(k),…,Yr(k)]T，U(k)＝[U1(k),U2(k),…,Up(k)]T，θ=[θ1T,θ2T,...θrT],]]>E(k)＝[ε1(k),ε2(k),…εr(k)]T；其中，Y(k)表示当前时刻供水管网系统输出；X表示系统过去时刻的输入和输出数据的集合；U(k)表示当前供水系统的输入控制变量；E(k)表示当前干扰造噪声，在参数估计中又定义为估计误差；k为当前递推步数；θ表示通过辨识得到的模型参数集合；T表示矩阵的转置；n，m，d+1分别为供水管网系统输出变量的阶次、输入变量的阶次、系统的时滞；步骤(2)实时估计模型参数，具体是：1)对于具有p个输入，r个输出的ARX模型，对应于模型第i个输出定义均方估计准则函数：VN(θi)=1NΣk=1Nλi(k)(Yi(k)-θiTX(k))2,i=1,2,...r]]>2)依据1)中准则函数，估计θi：θi(k)＝θi(k‑1)+mi(k)Pi(k‑1)X(k‑1)εi(k)ei(k)＝Y(k)‑θiT(k‑1)XHi(k)=Pi(k-1)X(k-1)X(k-1)TPi(k-1)λi(k)+X(k-1)TPi(k-1)X(k-1)]]>Pi(k)=1λi(k)[Pi(k-1)-m(k)Hi(k)]+m(k)Fi(k)]]>其中，Pi(k)，Fi(k)，Hi(k)分别对应于模型第i个输出的增益矩阵、最小增益矩阵及中间矩阵；且Pi(0)＝αI，Fi(k)＝βI，α,β＞0，I为单位矩阵；函数m(k)：m(k)=1,δ(k)=X(k-1)TP(k-1)X(k-1)>δ00,δ(k)≤δ0]]>其中，δ(k)表示k时刻，系统输入输出给模型估计过程带来的信息量；δ0为信息量阀值，设置δ0可提高估计过程的鲁棒性，且δ0＞0；λi(k)为自适应可变遗忘因子：λi(k)=λi,min+(1-λi,min)e-ρi|ϵi(k)|,i=1,2,...r]]>其中，ρi参数用于控制遗忘因子对误差εi(k)的敏感度；系统通过λi(k)的自适应调整，动态估计模型参数，适应系统的实时变化，且λi(k)∈[λi,min 1]，λi,min(k)∈(0 1]；3)根据θi，构造参数集合θ=[θ1T,θ2T,...θrT].]]>...

【技术特征摘要】
1.一种基于ARX的供水管网在线建模方法，其特征在于该方法包括以下步骤：步骤(1)、建立管网系统实时运行数据，建立系统的输入输出过程模型；具体是：1)建立供水管网系统实时运行数据库：实时运行数据库包括系统的实时输入数据和实时输出数据；系统的实时输入数据包括水厂供水压力、水厂供水量、管网中增压泵站压力和流量、管道上可控阀门的开启程度；系统实时输出数据包括压力监测点的压力、流量检测点的流量、蓄水池水位；通过SCADA系统采集N组实时过程运行数据，将采集的实时过程运行数据作为数据驱动的样本集合，表示为j＝1,2,…,N，Xj表示第j组系统输入数据，Yj表示第j组系统输出数据；设供水管网系统有p个输入控制变量，r个输出变量；2)以供水管网实时过程运行数据集合为基础，建立离散差分方程形式的ARX模型：Y(k)＝θTX(k)+E(k)其中：X(k)＝[YT(k-1),…,YT(k-n),UT(k-d-1),…,UT(k-d-m)]T，Y(k)＝[Y1(k),Y2(k),…,Yr(k)]T，U(k)＝[U1(k),U2(k),…,Up(k)]T，E(k)＝[ε1(k),ε2(k),…εr(k)]T；其中，Y(k)表示当前时刻供水管网系统输出；X表示系统过去时刻的输入和输出数据的集合；U(k)表示当前供水系统的输入控制变量；E(k)表示当前干扰造噪声，在参数估计中又定义为估计误差；k为当前递推步数；θ表示通过辨识得到的模型参数集合；T表示矩阵的转置；n，m，d+1分别为供水管网系统输出变量的阶次、输入变量的阶次、系统的时滞；步骤(2)实时估计模型参数，具体是：1)对于具有p个输入，r个输出的ARX模型，对应于模型第i个输出定义均方估计准则函数：

【专利技术属性】
技术研发人员：徐哲，杨洁，蔡华强，孔亚广，
申请(专利权)人：杭州电子科技大学，
类型：发明
国别省市：浙江;33