一种一维水动力模型泵站特性曲线系数的辨识方法技术

本发明专利技术公开了一种一维水动力模型泵站特性曲线系数的辨识方法，涉及水动力模型参数识别技术领域。方法包括：基于圣维南方程，建立目标明渠的一维明渠水动力模型，并率定待辨识泵站前、后连接渠道的糙率；待辨识泵站边界以连续方程和泵站特性曲线为相容方程，与S1中的圣维南方程耦合，形成包含待辨识泵站边界的一维渠道‑泵站水动力模型；构建以待辨识泵站的三个泵站特性曲线系数为状态向量的状态空间方程包括状态方程：X

An identification method for the coefficient of characteristic curve of pump station characteristic of one dimensional hydrodynamic model

The invention discloses a method for identifying the coefficient of characteristic curve of one dimensional hydrodynamic model pump station, and relates to the technical field of parameter identification of the hydrodynamic model. The method includes: Based on Saint Venant equations, one-dimensional hydrodynamic model to establish target channel channel, and the rate of connection channel roughness set to be identified, to be identified after the pump station; pumping station to boundary continuity equation and pump characteristic curve as the compatibility equation, and Saint Venant equations S1 coupling, which contains the identification of the one-dimensional channel boundary station pump hydrodynamic model; construction of pumping station for identification of three pumping stations characteristic curve coefficient as the state space equation of the state vector including state equation: X

【技术实现步骤摘要】
一种一维水动力模型泵站特性曲线系数的辨识方法
本专利技术涉及水动力模型参数识别
，尤其涉及一种一维水动力模型泵站特性曲线系数的辨识方法。
技术介绍
泵站是常见的渠系水工建筑物之一，其在灌溉排涝、跨流域调水、生产生活供水等方面发挥着重要作用。泵站特性曲线直接反映了泵站各相关水力要素的联系，是一维水动力模型中常用的泵站边界相容条件。因此，精确辨识泵站特性曲线系数，对提高一维明渠水动力模型的精度具有极其重要的价值。在包含泵站的河渠如平原河网、梯级泵站调水工程等一维水动力模拟中，水泵特性曲线系数常由模型水泵试验得到。但是，由于模型水泵和原型水泵之间存在比尺效应，导致水泵特性曲线系数与其真值存在差异。又由于泵站特性曲线随水泵机组开关、运行工况的变化而变化，而导致水泵特性曲线系数具有时变性。目前，现有技术中并没有能够实时辨识泵站特性曲线系数的方法。因此，需要一种研究泵站特性曲线系数的自动辨识方法，以提高含泵站河网水动力模型的精度和适应性。
技术实现思路
本专利技术的目的在于提供一种一维水动力模型泵站特性曲线系数的辨识方法，从而解决现有技术中存在的前述问题。为了实现上述目的，本专利技术所述一维水动力模型泵站特性曲线系数的辨识方法，所述方法包括：S1，基于圣维南方程，建立目标明渠的一维明渠水动力模型，并以此模型为基础，率定目标明渠中待辨识泵站的泵站前连接渠道的糙率和泵站后连接渠道的糙率；S2，待辨识泵站边界以连续方程和泵站特性曲线为相容方程，与S1中的圣维南方程耦合，形成包含待辨识泵站边界的一维渠道-泵站水动力模型，求解所述一维渠道-泵站水动力模型得到所述目标明渠的水动力过程；所述连续方程为公式(1)：Qj＝Qj+1(1)其中，j表示待辨识泵站进水侧的断面编号，j+1表示待辨识泵站出水侧的断面编号；所述泵站特性曲线为公式(2)：其中，Qj为泵站进水侧流量，Qj+1为泵站出水侧流量，Zj为进水侧水位，Zj+1为出水侧水位，A，B，C均为泵站特性曲线系数；S3，所述目标明渠的水动力过程的基础上，构建以待辨识泵站的三个泵站特性曲线系数为状态向量的状态空间方程，所述状态空间方程包括状态方程和量测方程；所述状态方程为公式(3)：Xn+1＝ΦnXn+wn+1(3)；所述量测方程为公式(4)：Yn+1＝HXn+1+vn+1(4)；公式(3)中，Xn+1为n+1时刻的状态向量，Xn+1＝[An+1Bn+1Cn+1]T；Φn为n时刻到n+1时刻的状态转移矩阵；Xn为n时刻的状态向量，Xn＝[AnBnCn]T，T表示行向量转置列向量的转置符号；wn+1为n+1时刻的扰动矩阵；公式(4)中，Yn+1为n+1时刻的观测向量，即n+1时刻待辨识泵站的水力要素观测数据；H为一维渠道-泵站水动力模型，记为观测算子；vn+1为n+1时刻的观测噪声向量；S4，在泵站机组角度或转速不变时段，获取目标明渠中待辨识泵站现场的水力要素观测数据并在该水力要素观测数据的基础上，采用集合卡尔曼滤波器同化状态空间方程，得到同化后的泵站特性曲线系数。优选地，所述Φn为单位矩阵I。优选地，量测方程的观测矩阵H为非线性观测算子。本专利技术所述方法以一维水动力模型中的泵站特性曲线系数为研究对象，构建以所述系数为状态向量的状态空间方程，采用集合卡尔曼滤波技术，获取同化系数，为提高河网水动力数值模拟的精度提供了科学依据和技术支撑。本专利技术的有益效果是：(1)现有技术一般采用模型试验获取泵站特性曲线系数，缺点为：试验周期长。而本专利技术所述方法利用当时水力要素的观测数据，调用同化模型，实现了在一维明渠水动力模拟过程中实时调整泵站特性调整曲线系数。(2)现有技术采用模型试验方法获得的泵站特性曲线系数，存在比尺效应，直接应用于一维渠道-泵站水动力模型中，模拟精度低。而本专利技术采用观测数据反演得到泵站特性曲线系数，显著提高水动力模拟结果的精度。(3)现有技术一般只能得到特定工况的泵站特性曲线系数，而本专利技术可以得到任意工况下的泵站特性曲线系数。附图说明图1为本专利技术方法流程图；图2为渠道-泵站系统示意图；图3为集合卡尔曼滤波同化泵站特性曲线系数流程图；图4为埝头泵站进水侧水位模拟结果；图5为埝头泵站进水侧水位模拟结果；图6为埝头泵站流量模拟结果。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施方式仅仅用以解释本专利技术，并不用于限定本专利技术。关于本申请所述方法更详细的解释说明：图1为本专利技术方法流程图，本专利技术所述一维水动力模型泵站特性曲线系数的辨识方法，具体技术方案如下：S1，基于圣维南方程(SaintVenant方程)，建立目标明渠的一维明渠水动力模型，并率定目标明渠中待辨识曲线系数泵站的泵站前连接渠道的糙率和泵站后连接渠道的糙率；一维明渠水动力模型用于渠道模拟，根据渠道模拟结果判断率定糙率是否合理。不考虑旁侧入流情况下，明渠一维非恒定流的计算通常采用SaintVenant方程组描述，该方程包括连续方程和动量方程：其中，x和t分别为空间坐标和时间坐标；A为渠道处的过流面积；Q为渠道流量；Z为水位；α为动量校正修正；g为重力加速度；C为谢才系数；R为水力半径。采用Pressimann四点隐式差分格式离散SaintVenant方程，并用追赶法求解。该方法属于已有技术。步骤S1中，率定目标明渠中待辨识曲线系数泵站的泵站前连接渠道的糙率和泵站后连接渠道的糙率时，以渠池上游边界水位的模拟值与实测值之差最小为目标。如图2为渠道-泵站系统示意图，假设渠道顺直，无其他局部建筑物。选取上游侧渠池和下游侧渠池，分别进行渠道糙率率定。率定的方法可以是人工试错法、优化算法等。各相邻渠道糙率率定完成后，可以认为渠道-泵站系统水动力模拟误差基本来自泵站边界的参数。S2，以连续方程和泵站特性曲线为相容条件，与SaintVenant方程耦合求解。在不考虑泵站内部流量损失的情况下，泵站入口流量和泵站出口流量基本相等，满足连续性条件。泵站内部多水泵机组运行工况变化时，水泵特性曲线会随之变化，导致泵站综合特性曲线也会随之变化。泵站特性曲线表征的是泵站的水位差和流量关系，属于泵站的固有属性。因此，以连续方程(式2)和泵站特性曲线(式3)作为相容条件，与SaintVenant方程耦合，求解含泵站河网的水动力过程。Qj＝Qj+1(2)式中，Qj为泵站进水侧流量，Qj+1为泵站出水侧流量，Zj为进水侧水位，Zj+1为出水侧水位，A，B，C均为泵站特性曲线系数。若泵站中仅开一台机组，则系数与水泵机组系数相同；若泵站开多台机组，则曲线系数为多水泵的组合转角系数。运行工况变化时，组合系数随之发生变化。步骤3，选取泵站进水侧和出水侧的水力要素观测数据作为观测量，构建以泵站特性曲线系数A，B，C为状态向量的状态空间方程，即状态方程和量测方程。泵站特性曲线系数对水动力模拟结果的精度有重要的影响，特别是泵站作为内边界时，模拟误差会向上下游传递，影响更为显著。水力要素观测数据为相同观测断面的水位或流量或不同观测断面的水位或流量。同时以泵站特性曲线系数A，B，C为状态向量，以X表示，构建状态空间方程。状态方程见式(4)：Xn+1＝ΦnXn+wn+1(4)式中，Xn+1为n+1时刻的状本文档来自技高网...

【技术保护点】
一种一维水动力模型泵站特性曲线系数的辨识方法，其特征在于，所述方法包括：S1，基于圣维南方程，建立目标明渠的一维明渠水动力模型，并以此模型为基础，率定目标明渠中待辨识泵站的泵站前连接渠道的糙率和泵站后连接渠道的糙率；S2，待辨识泵站边界以连续方程和泵站特性曲线为相容方程，与S1中的圣维南方程耦合，形成包含待辨识泵站边界的一维渠道‑泵站水动力模型，求解所述一维渠道‑泵站水动力模型得到所述目标明渠的水动力过程；所述连续方程为公式(1)：Qj＝Qj+1   (1)其中，j表示待辨识泵站进水侧的断面编号，j+1表示待辨识泵站出水侧的断面编号；所述泵站特性曲线为公式(2)：

【技术特征摘要】
1.一种一维水动力模型泵站特性曲线系数的辨识方法，其特征在于，所述方法包括：S1，基于圣维南方程，建立目标明渠的一维明渠水动力模型，并以此模型为基础，率定目标明渠中待辨识泵站的泵站前连接渠道的糙率和泵站后连接渠道的糙率；S2，待辨识泵站边界以连续方程和泵站特性曲线为相容方程，与S1中的圣维南方程耦合，形成包含待辨识泵站边界的一维渠道-泵站水动力模型，求解所述一维渠道-泵站水动力模型得到所述目标明渠的水动力过程；所述连续方程为公式(1)：Qj＝Qj+1(1)其中，j表示待辨识泵站进水侧的断面编号，j+1表示待辨识泵站出水侧的断面编号；所述泵站特性曲线为公式(2)：其中，Qj为泵站进水侧流量，Qj+1为泵站出水侧流量，Zj为进水侧水位，Zj+1为出水侧水位，A，B，C均为泵站特性曲线系数；S3，所述目标明渠的水动力过程的基础上，构建以待辨识泵站的三个泵站特性曲线系数为状态向量的状态空间方程，所述状态空间方程包括状态方程和量测方程；所...

【专利技术属性】
技术研发人员：雷晓辉，张召，田雨，王浩，孟现勇，王超，张忠波，杨明祥，蔡思宇，
申请(专利权)人：中国水利水电科学研究院，
类型：发明
国别省市：北京,11