[发明专利]一种基于BRB系统的供水管网建模方法

权利要求书

1.一种基于BRB系统的供水管网建模方法，其特征在于该方法包括以下步骤：

步骤1.收集管网运行的历史数据，具体是：

从供水管网的SCADA系统和数据库系统中获取管网运行的历史数据，所述的历史数据包括水厂出厂压力、水厂供水量、管网中增压泵站压力、管网中增压泵站流量，管道上可控阀门的开启程度、压力监测点的压力、流量监测点的流量、蓄水池水位和用户用水量；

步骤2.确定供水管网模型的输入变量和输出变量，具体是：

（2-1）从实际获得的管网运行的历史数据出发，选取供水管网模型输入输出变量，输入变量包括：水厂出厂压力、水厂供水量、管网中增压泵站压力、管网中增压泵站流量，管道上可控阀门的开启程度；输出变量包括：压力监测点的压力、流量检测点的流量和蓄水池水位；

（2-2）根据选取的输入变量和输出变量，对管网运行的历史数据进行筛选，得到仅含输入变量和输出变量的历史数据，作为样本数据；同时，将样本数据分成两部分，前一部分作为训练数据集，后一部分作为检验数据集；

步骤3.建立基于BRB的供水管网模型，具体是：

基于BRB系统的模型为MISO模型，供水管网为MIMO模型；因此，基于BRB的供水管网MIMO模型为多个基于BRB的供水管网MISO模型的合成；基于BRB的供水管网MISO模型的建立步骤如下：

（3-1）确定规则库中输入向量和输出向量

规则库中输入向量即为步骤2中选择的输入变量构成的向量，规则库中输出向量为步骤2中选择的输出变量的某一个；并设输入向量U＝(U_i,i＝1,...,T)，输出变量为O；

（3-2）建立初始规则库

BRB系统中规则库中的规则为：

Rk:IFU1isA1kΛU2isA2kΛ...ΛUTkisATkk]]>

THEN{(D1,β1k),(D2,β2k),···,(DN,βNk)},(Σi=1Nβik≤1),k=1,...L]]>

规则的权重：θ_k(k＝1,...,L)，输入变量权重：δ₁,δ₂,...,

其中，L为规则的总数，为输入向量，为第k条规则中第i输入变量的参考值，D＝(D₁,D₂,...,D_N)为推理输出向量，β_ik(i＝1,...,T_k)为推理输出D_i的信任度；如果则称第k条规则是完整的，如果则第k条规则被完全忽略；θ_k(k＝1,...,L)为第k规则在规则库中的相对权重，δ₁,δ₂,...,为第k规则中各输入变量的相对权重；

规则库中可变参数为P＝(β_ik,θ_k,δ_j;i＝1,...,N;k＝1,...,L;j＝1,...,T)；为建立初始化规则库，随机初始β_ik，而θ_k，δ_j初始设定为1；

（3-3）训练规则库

初始规则库通过学习训练数据集中的输入数据集和输出数据集的映射关系，自动调整可变参数P＝(β_ik,θ_k,δ_j;i＝1,...,N;k＝1,...,L;j＝1,...,T)，最终获得描述输入向量U＝(U_i,i＝1,...,T)和输出变量O关系的规则库；

规则库中的推理过程如下：

Step1：将输入向量U＝(U_i,i＝1,...,T)转换成信任分布：

S(U_i)＝{(A_i,j,α_i,j),j＝1,...,J_i}

其中J_i为参考值的个数，α_i,j为输入U_i与参考值A_i,j的相似度；α_i,j计算如下：

如果A_i,j≤U_i≤A_i,j+1，则α_i,j+1＝1-α_i,j且对于j’＝1,...,J_i，j’≠j,j+1，有α_i,j’＝0

Step2：计算第k规则被激活的权重w_k：

wk=θk×Πi=1Tk(αi,jk)δi‾Σj=1L[θj×Πl=1Tk(αi,jl)δl‾],δi‾=δimaxi=1,...,Tk{δi}]]>

Step3：计算输出变量O：

先计算输出变量的信任分布O(U)＝{(D_n,β_n(U)),n＝1,...,N}，其中

{Dn}:βn(U)=μ×[Πk=1K(wk(U)βn,k+1-wk(U)Σi=1Nβi,k)-Πk=1K(1-wk(U)Σi=1Nβi,k)]]]>

{D}:βD(U)=μ×[Πk=1K(1-wk(U)Σi=1Nβi,k)-Πk=1K(1-wk(U))]]]>

μ=Σj=1NΠk=1L(wkβjk+1-wkΣj=1Nβjk)-(N-1)Πk=1L(1-wkΣj=1Nβjk)-Πk=1K(1-wk(U))]]>

β_D表示未知推理结果D_n的信任度，且满足

如果信任规则是完整的，即有β_D(U)＝0，则最终的输出为K条规则的综合，表为O(U)＝{(D_n,β_n(U)),n＝1,...,N}，数值输出为O(U)=Σn=1NDnβn(U);]]>

如果信任规则不完整，即有β_D(U)＞0，则最终的数值输出具有上下界O(U)∈[O_min(U),O_max(U)]，其中

Omin(U)=Σn=1N-1u(Dn)βn(U)+u(DN)(βN(U)+βD(U))]]>

Omax(U)=u(D1)(β1(U)+βD(U))+Σn=2Nu(Dn)βn(U)]]>

至此，给定一个输入U＝(U_i,i＝1,...,T)到规则库就有一个输出O；对于信任规则是完整情形，规则库的训练就是求解如下优化问题，其中H为训练样本个数，为训练样本实际输出；

0≤βik,θk,δj≤1Σi=1Nβik=1]]>

对于信任规则不完整情形，规则库的训练就是求解如下优化问题

（3-4）检验规则库

利用检验数据对（3-1）（3-2）（3-3）建立的规则库进行检验，以评价模型的精度；精度的标准用来衡量；若ε(P)不符合建模精度要求，重复（3-1）（3-2）（3-3）重新建模；

至此，基于BRB的供水管网MISO模型已建立；

步骤4.建立基于BRB的供水管网MIMO模型

基于BRB的供水管网MIMO模型为多个基于BRB的供水管网MISO模型的合成；所以，根据输出变量的个数，为每个输出变量重复MISO模型的建模过程。

2.根据权利要求1所述的一种基于BRB系统的供水管网建模方法，其特征在于：所述的供水管网模型输入还包括用户用水量及系统输入噪声。