本发明涉及供热管网设计，特别是涉及一种供热管网水力仿真模型参数辨识修正方法。
背景技术：
：：1、随着供热规模日益扩大和环状网络结构日趋复杂，建立严格精确供热管网水力仿真模型变得愈加困难。供热管网阻抗作为管网运行重要参数，对管网运行调节起着至关重要作用。实际工程中若未对管网阻抗辨识，基于理论计算的管网阻抗值多与实际不符，沿程阻力系数管段粗糙度的变化而随管网运行不断变化，而局部阻力系数的估计值则极不准确，对管网水力分析产生较大影响。2、目前国内外对管网阻抗辨识研究多偏向于枝状管网，对环状管网阻抗辨识的研究主要利用稳态水力方程，应用基于最小二乘法或梯度下降方法实现，该方法需要大量稳态工况流量和压力检测值，为管网智能化监控、改造、后期数据处理都有极高要求；且管道内压力值的巨大波动和误差，也导致该方法难以实现工程应用。3、目前对环状管网水力特性研究中，环状管网阻抗辨识问题的难点在于环状管网管段流量、流向难以监测，无法对各管段压差、流量进行辨识。管网成环原因主要是供（回）水网中出现平面的环状结构以及热力站连接了供回水管网形成了“热源-供水网-热力站-回水网-热源”这个拓扑意义上的环状结构。4、对于环状管网各管段流量，有且只有与管网拓扑环数相同的流量值之间是线性无关的，也就是余枝管段的流量。余枝管段中，仅平面环状结构的部分流量是不可检测的，热源及热力站中的监测设备较为全面，能够得到阀门开度与水泵频率的实时参数，上述条件形成了基于动态模型进行状态观测的基础。然而，如何基于该些条件进行动态监测，以满足工程需要，则是须要解决的难题。技术实现思路1、本发明的目的是针对现有技术中存在的技术缺陷，而提供一种供热管网水力仿真模型参数辨识修正方法。2、为实现本发明的目的所采用的技术方案是：3、一种供热管网水力仿真模型参数辨识修正方法，包括：4、s1.基于集总参数假设的集中供热非稳态水力模型，构建线性化模型；5、s2.基于线性化模型，建立集中供热管网动态系统的线性模型；6、s3.判断线性模型是否满足能观性，若是，进入步骤s4，若否进入步骤s8；7、s4.构造状态估计增益矩阵；8、s5.利用状态估计增益矩阵，基于实时传感器数据进行集中供热管网的状态估计，得到余枝管段流量的状态估计值；9、s6.基于余枝管段流量的状态估计值，判断状态估计误差与速度是否满足预设要求，若是则进入步骤s7；若否进入步骤s9；10、s7.基于余枝管段流量的状态估计值，计算管网各管段的阻抗的估计值；11、s8.调整观测的余枝管段，修正基于管网拓扑信息所得到供热管网的基本回路矩阵，然后返回步骤s2；12、s9.调整构造状态估计增益矩阵的特征值及状态估计增益系数，然后返回步骤s4。13、步骤s1中，所述线性化模型表达式如下：14、<mstyledisplaystyle="true"mathcolor="#000000"><msub><moveraccent="true"><mi>q</mi><mo>˙</mo></mover><mi>r</mi></msub><mi>=−</mi><msup><mrow><mo>(</mo><mrow><msub><mi>c</mi><mi>f</mi></msub><msub><mi>k</mi><mi>f</mi></msub><msubsup><mi>c</mi><mi>f</mi><mi>t</mi></msubsup></mrow><mo>)</mo></mrow><mrow><mi>−</mi><mn>1</mn></mrow></msup><msub><mi>c</mi><mi>f</mi></msub><mi>s</mi><msup><mrow><mo>|</mo><mrow><msubsup><mi>c</mi><mstyledisplaystyle="true"mathcolor="#000000"><mi>f</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>t</mi></mstyle></msubsup><msubsup><mi>q</mi><mstyledisplaystyle="true"mathcolor="#000000"><mi>r</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>*</mi></mstyle></msubsup></mrow><mo>|</mo></mrow><mi>d</mi></msup><msubsup><mi>c</mi><mstyledisplaystyle="true"mathcolor="#000000"><mi>f</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>t</mi></mstyle></msubsup><msub><mi>q</mi><mstyledisplaystyle="true"mathcolor="#000000"><mi>r</mi></mstyle></msub><mo>+</mo><msup><mrow><mo>(</mo><mrow><msub><mi>c</mi><mi>f</mi></msub><msub><mi>k</mi><mi>f</mi></msub><msubsup><mi>c</mi><mi>f</mi><mi>t</mi></msubsup></mrow><mo>)</mo></mrow><mrow><mi>−</mi><mn>1</mn></mrow></msup><msub><mi>c</mi><mi>f</mi></msub><msup><mrow><mo>[</mo><mrow><mstyledisplaystyle="true"mathcolor="#000000"><mi>h</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><msub><mrow/><mstyledisplaystyle="true"mathcolor="#000000"><mi>p</mi></mstyle></msub></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>−</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>h</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><msub><mrow/><mstyledisplaystyle="true"mathcolor="#000000"><mi>ϕ</mi></mstyle></msub></mstyle></mrow><mo>]</mo></mrow><mstyledisplaystyle="true"mathcolor="#000000"><mi>d</mi></mstyle></msup></mstyle>，15、其中系数由下式计算：16、，17、其中，表示基于管网拓扑信息所得到供热管网的基本回路矩阵，表示基本回路矩阵的转置，为各余枝管段流量值， 表示对时间的导数，为预定稳态工况下的余枝管段流量值，为管段阻抗矩阵；为水泵扬程；表征调节阀阻力，为管中水流密度；为管段长度；表征管段内径，<mstyledisplaystyle="true"mathcolor="#000000"><msup><mrow><mo>[</mo><mrow/><mo>]</mo></mrow><mstyledisplaystyle="true"mathcolor="#000000"><mi>d</mi></mstyle></msup></mstyle>表示对角阵。18、步骤s2中，所述集中供热管网动态系统的线性模型的表达式如下：19、<mstyledisplaystyle="true"mathcolor="#000000"><msub><moveraccent="true"><mi>q</mi><mo>˙</mo></mover><mi>r</mi></msub><mi>=−</mi><msup><mrow><mo>(</mo><mrow><msub><mi>c</mi><mi>f</mi></msub><msub><mi>k</mi><mi>f</mi></msub><msubsup><mi>c</mi><mi>f</mi><mi>t</mi></msubsup></mrow><mo>)</mo></mrow><mrow><mi>−</mi><mn>1</mn></mrow></msup><msub><mi>c</mi><mi>f</mi></msub><mi>s</mi><msup><mrow><mo>|</mo><mrow><msubsup><mi>c</mi><mi>f</mi><mi>t</mi></msubsup><msubsup><mi>q</mi><mstyledisplaystyle="true"mathcolor="#000000"><mi>r</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>*</mi></mstyle></msubsup></mrow><mo>|</mo></mrow><mi>d</mi></msup><msubsup><mi>c</mi><mi>f</mi><mi>t</mi></msubsup><msub><mi>q</mi><mi>r</mi></msub><mo>+</mo><msup><mrow><mo>(</mo><mrow><msub><mi>c</mi><mi>f</mi></msub><msub><mi>k</mi><mi>f</mi></msub><msubsup><mi>c</mi><mi>f</mi><mi>t</mi></msubsup></mrow><mo>)</mo></mrow><mrow><mi>−</mi><mn>1</mn></mrow></msup><msub><mi>c</mi><mi>f</mi></msub><mi>[</mi><msub><mi>h</mi><mi>p</mi></msub><mi>−</mi><msub><mi>h</mi><mstyledisplaystyle="true"mathcolor="#000000"><mi>ϕ</mi></mstyle></msub><mi>]</mi></mstyle>，20、将集中供热管网动态系统的线性模型简写为如下：21、，22、<mstyledisplaystyle="true"mathcolor="#000000"><mo>{</mo><mtablecolumnalign="left"><mtr><mtd><mi>a</mi><mi>=−</mi><msup><mrow><mo>(</mo><mrow><msub><mi>c</mi><mi>f</mi></msub><msub><mi>k</mi><mi>f</mi></msub><msubsup><mi>c</mi><mi>f</mi><mi>t</mi></msubsup></mrow><mo>)</mo></mrow><mrow><mi>−</mi><mn>1</mn></mrow></msup><msub><mi>c</mi><mi>f</mi></msub><msubsup><mi>sc</mi><mstyledisplaystyle="true"mathcolor="#000000"><mi>f</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>t</mi></mstyle></msubsup></mtd></mtr><mtr><mtd><mstyledisplaystyle="true"mathcolor="#000000"><mi>b</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>=</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mo>(</mo><mrow><msub><mi>c</mi><mi>f</mi></msub><msub><mi>k</mi><mi>f</mi></msub><msubsup><mi>c</mi><mi>f</mi><mi>t</mi></msubsup></mrow><mo>)</mo></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><msup><mrow/><mrow><mi>−</mi><mn>1</mn></mrow></msup></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>c</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><msub><mrow/><mi>f</mi></msub></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>⋅</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mo>[</mo><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mi>−</mi><mn>1</mn></mtd></mtr></mtable><mo>]</mo></mstyle></mtd></mtr><mtr><mtd><mi>u</mi><mi>=</mi><mrow><mo>[</mo><mtablecolumnalign="left"><mtr><mtd><mstyledisplaystyle="true"mathcolor="#000000"><mi>h</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><msub><mrow/><mstyledisplaystyle="true"mathcolor="#000000"><mi>p</mi></mstyle></msub></mstyle></mtd></mtr><mtr><mtd><msub><mi>h</mi><mstyledisplaystyle="true"mathcolor="#000000"><mi>ϕ</mi></mstyle></msub></mtd></mtr></mtable><mo>]</mo></mrow></mtd></mtr></mtable></mstyle>，23、式中，为可监测的余枝管段流量值，表示矩阵，为可监测的余枝管段的数量，矩阵的每列对应余枝管段，各行元素值表示对应余枝管段是否可监测，为的控制向量，矩阵为矩阵，矩阵为矩阵，为环数，为管段数。24、步骤s3中，判断线性模型是否满足能观性，通过构造可观性矩阵判断实现；其中，基于下式构造可观性矩阵，若可观性矩阵列秩为，则认为线性模型满足能观性：25、<mstyledisplaystyle="true"mathcolor="#000000"><mo>[</mo><mtable><mtr><mtd><mi>c</mi></mtd></mtr><mtr><mtd><mi>ca</mi></mtd></mtr><mtr><mtd><msup><mi>ca</mi><mstyledisplaystyle="true"mathcolor="#000000"><mn>2</mn></mstyle></msup></mtd></mtr><mtr><mtd><mi>⋮</mi></mtd></mtr><mtr><mtd><msup><mi>ca</mi><mrow><mstyledisplaystyle="true"mathcolor="#000000"><mi>n</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><msub><mrow/><mstyledisplaystyle="true"mathcolor="#000000"><mi>r</mi></mstyle></msub></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>−</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mn>1</mn></mstyle></mrow></msup></mtd></mtr></mtable><mo>]</mo></mstyle>， 。26、步骤s4中，所述的状态估计增益矩阵的构造过程包括：27、对有负实部的特征值集，构造一个矩阵，满足且，即矩阵和矩阵没有相同特征值；表示维度为的单位矩阵；28、构造矩阵，满足是能控的，使得构造的能控性矩阵的秩为，能控性矩阵由下式构造；29、<mstyledisplaystyle="true"mathcolor="#000000"><mo>[</mo><mtable><mtr><mtd><mstyledisplaystyle="true"mathcolor="#000000"><moveraccent="true"><mi>l</mi><mo>¯</mo></mover></mstyle></mtd><mtd><mstyledisplaystyle="true"mathcolor="#000000"><mi>f</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><moveraccent="true"><mi>l</mi><mo>¯</mo></mover></mstyle></mtd><mtd><mstyledisplaystyle="true"mathcolor="#000000"><mi>f</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><msup><mrow/><mstyledisplaystyle="true"mathcolor="#000000"><mn>2</mn></mstyle></msup></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><moveraccent="true"><mi>l</mi><mo>¯</mo></mover></mstyle></mtd><mtd><mi>⋯</mi></mtd><mtd><mstyledisplaystyle="true"mathcolor="#000000"><mi>f</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><msup><mrow/><mrow><mstyledisplaystyle="true"mathcolor="#000000"><mi>n</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><msub><mrow/><mstyledisplaystyle="true"mathcolor="#000000"><mi>r</mi></mstyle></msub></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>−</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mn>1</mn></mstyle></mrow></msup></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><moveraccent="true"><mi>l</mi><mo>¯</mo></mover></mstyle></mtd></mtr></mtable><mo>]</mo></mstyle>，30、求lyapunov方程唯一解，为矩阵，验证矩阵是非奇异的；31、按下式计算状态估计增益矩阵；32、。33、步骤s5中，通过以下式得到余枝管段流量的状态估计值：34、，35、，36、其中， 表示第时刻的各余枝管段流量的状态估计值对时间的导数，表示各余枝管段流量的状态估计值，是状态估计误差，表示第时刻的状态估计误差对时间的导数，表示时刻的状态估计误差，的初值是；随着时刻的可监测的余枝管段流量值的不断输入，不断得到新的余枝管段流量的状态估计值。37、步骤s6中，判断状态估计误差与速度是否满足预设要求，是通过判断可监测的余枝管段流量值与状态估计值的误差的二范数是否降低到预设值来实现；表达式如下：38、；39、式中，表示时刻的可监测的余枝管段流量值，表示时刻的各余枝管段流量的状态估计值。40、步骤s7，采用以下式计算管网各管段的阻抗的估计值：41、<mstyledisplaystyle="true"mathcolor="#000000"><mo>{</mo><mtablecolumnalign="left"><mtr><mtd><moveraccent="true"><mi>s</mi><mo>^</mo></mover><mi>=</mi><msup><mrow><mo>[</mo><mrow><mi>|</mi><msup><mi>c</mi><mstyledisplaystyle="true"mathcolor="#000000"><mi>t</mi></mstyle></msup><msub><moveraccent="true"><mi>q</mi><mo>^</mo></mover><mstyledisplaystyle="true"mathcolor="#000000"><mi>r</mi></mstyle></msub><mi>|</mi><msup><mi>c</mi><mstyledisplaystyle="true"mathcolor="#000000"><mi>t</mi></mstyle></msup><msub><moveraccent="true"><mi>q</mi><mo>^</mo></mover><mstyledisplaystyle="true"mathcolor="#000000"><mi>r</mi></mstyle></msub></mrow><mo>]</mo></mrow><mrow><mstyledisplaystyle="true"mathcolor="#000000"><mi>−</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mn>1</mn></mstyle></mrow></msup><mrow><mo>[</mo><mrow><mstyledisplaystyle="true"mathcolor="#000000"><moveraccent="true"><mi>h</mi><mo>^</mo></mover></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><msub><mrow/><mi>p</mi></msub></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>−</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><moveraccent="true"><mi>h</mi><mo>^</mo></mover></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><msub><mrow/><mstyledisplaystyle="true"mathcolor="#000000"><mi>ϕ</mi></mstyle></msub></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>−</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>k</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><msub><mrow/><mstyledisplaystyle="true"mathcolor="#000000"><mi>f</mi></mstyle></msub></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><moveraccent="true"><moveraccent="true"><mi>q</mi><mo>^</mo></mover><mo>˙</mo></mover></mstyle></mrow><mo>]</mo></mrow></mtd></mtr><mtr><mtd><mstyledisplaystyle="true"mathcolor="#000000"><moveraccent="true"><mi>h</mi><mo>^</mo></mover></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><msub><mrow/><mi>p</mi></msub></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>=</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mo>[</mo><mrow><mi>a</mi><msup><mrow><mo>|</mo><mrow><msubsup><mi>c</mi><mstyledisplaystyle="true"mathcolor="#000000"><mi>f</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>t</mi></mstyle></msubsup><msub><moveraccent="true"><mi>q</mi><mo>^</mo></mover><mstyledisplaystyle="true"mathcolor="#000000"><mi>r</mi></mstyle></msub></mrow><mo>|</mo></mrow><mi>d</mi></msup><msubsup><mi>c</mi><mstyledisplaystyle="true"mathcolor="#000000"><mi>f</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>t</mi></mstyle></msubsup><mo>+</mo><mfrac><mi>b</mi><msub><mi>r</mi><mstyledisplaystyle="true"mathcolor="#000000"><mi>f</mi></mstyle></msub></mfrac><msubsup><mi>c</mi><mstyledisplaystyle="true"mathcolor="#000000"><mi>f</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>t</mi></mstyle></msubsup></mrow><mo>]</mo></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><msup><mrow/><mstyledisplaystyle="true"mathcolor="#000000"><mi>d</mi></mstyle></msup></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><moveraccent="true"><mi>q</mi><mo>^</mo></mover></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><msub><mrow/><mstyledisplaystyle="true"mathcolor="#000000"><mi>r</mi></mstyle></msub></mstyle></mtd></mtr><mtr><mtd><mstyledisplaystyle="true"mathcolor="#000000"><moveraccent="true"><mi>h</mi><mo>^</mo></mover></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><msub><mrow/><mstyledisplaystyle="true"mathcolor="#000000"><mi>ϕ</mi></mstyle></msub></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>=</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mo>[</mo><mfrac><msup><mi>r</mi><mrow><mstyledisplaystyle="true"mathcolor="#000000"><mn>2</mn></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>(</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mn>1</mn></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>−</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>x</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>)</mi></mstyle></mrow></msup><mrow><mstyledisplaystyle="true"mathcolor="#000000"><mi>c</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><msub><mrow/><mstyledisplaystyle="true"mathcolor="#000000"><mi>v</mi></mstyle></msub></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><msup><mrow/><mstyledisplaystyle="true"mathcolor="#000000"><mn>2</mn></mstyle></msup></mstyle></mrow></mfrac><mo>]</mo></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><msup><mrow/><mstyledisplaystyle="true"mathcolor="#000000"><mi>d</mi></mstyle></msup></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mo>|</mo><mrow><msubsup><mi>c</mi><mstyledisplaystyle="true"mathcolor="#000000"><mi>f</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>t</mi></mstyle></msubsup><msub><moveraccent="true"><mi>q</mi><mo>^</mo></mover><mstyledisplaystyle="true"mathcolor="#000000"><mi>r</mi></mstyle></msub></mrow><mo>|</mo></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><msup><mrow/><mi>d</mi></msup></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>c</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><munderover><mi/><mstyledisplaystyle="true"mathcolor="#000000"><mi>f</mi></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><mi>t</mi></mstyle></munderover></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><moveraccent="true"><mi>q</mi><mo>^</mo></mover></mstyle><mstyledisplaystyle="true"mathcolor="#000000"><msub><mrow/><mstyledisplaystyle="true"mathcolor="#000000"><mi>r</mi></mstyle></msub></mstyle></mtd></mtr></mtable></mstyle>，42、，43、式中，表示管网各管段阻抗的估计值，为水泵流量特性曲线系数，为水泵变频比，为调节阀的阀门开度；为调节阀的阀门可调比，为调节阀的阀门流通能力，<mstyledisplaystyle="true"mathcolor="#000000"><msup><mrow><mo>[</mo><mrow/><mo>]</mo></mrow><mstyledisplaystyle="true"mathcolor="#000000"><mi>d</mi></mstyle></msup></mstyle>表示对角阵，表示调节阀阻力的估计值，表示矩阵c的转置。44、步骤s9中，所述调整构造状态估计增益矩阵的特征值及状态估计增益系数，包括：45、调整矩阵的特征值，修正矩阵以及调整状态估计增益矩阵的状态估计增益系数。46、本发明将集中供热管网视为动态系统，构建以水泵扬程、调节阀阻力为输入，可检测的余枝管段流量值为输出的线性模型，基于该线性模型进行状态估计，对集中供热管网中没有监测设备的管段的不可监测流量以及管网阻抗值进行估计，实现了仅利用连续的流量监测数据，对集中供热管网阻抗的辨识，从而完成了对集中供热管网水力模型的修正与辨识。47、本发明方法仅需利用流量检测值和水泵频率、调节阀的阀门开度等可检测量，不再依赖压力检测值，极大地简化了水力仿真模型的应用要求，规避了压力监测值误差和波动较大的问题。48、本发明方法利用动态线性模型，故仅一个稳态工况附近的，连续的动态的流量检测值，避免了对大量稳态工况数据的统计要求。49、本发明方法基于线性模型的状态估计方法，降低了对现场传感器的布置要求，对大型环状管网，仅需热源、热力站和部分中继泵站的流量值，不再需要在供水网或回水网平面拓扑结构中，每个环都布置传感器，从而节约硬件成本，且方便布置。当前第1页12当前第1页12