一种基于脉冲响应序列的输气管道动态仿真方法技术

本发明专利技术公开了一种基于脉冲响应模型的输气管道动态仿真方法，首先建立管道脉冲响应序列模型；然后，获取管道压力脉冲响应序列，在管道模型中分别阶跃改变管道的入口流量和出口流量，观察管道的压力变化，将每个阶跃响应序列延迟一个时刻，用原序列减延迟序列获得四个脉冲响应序列；然后，进行管道压力脉冲响应仿真测试，包括：压力阶跃响应测试、单输入压力响应测试、双输入压力响应测试；最后，选择初始参考稳态点，根据管道历史采样数据，用脉冲响应实现在线动态仿真，包括单输入单输出离散系统的在线动态仿真和双输入双输出离散系统的在线动态仿真。

【技术实现步骤摘要】
一种基于脉冲响应序列的输气管道动态仿真方法
本专利技术属于动态仿真
，尤其涉及一种基于脉冲响应序列的输气管道动态仿真方法。
技术介绍
现有的管道仿真都是基于管道的连续性方程、运动方程、能量方程及气体状态方程等机理模型进行动态仿真计算，由于机理模型是一个多变量偏微分方程组，无法直接求解，需要转化为常微分方程组，并且对长输管道进行剖分，在时间和空间上进行离散化，计算过程需要迭代求解，计算量大，并且当管道初始条件设置不好时可能导致计算发散。
技术实现思路
本专利技术的目的在于提供一种基于脉冲响应序列的输气管道动态仿真方法，旨在解决利用管道机理模型计算量大、易发散的问题，提高管道动态仿真计算速度，实现管道在线动态仿真。本专利技术是这样实现的，一种基于脉冲响应序列的输气管道动态仿真方法包括：步骤一、建立管道脉冲响应序列模型，管道的入口压力Pi和出口压力Po按下式计算：步骤二、获取管道压力脉冲响应序列，在管道模型中分别阶跃改变管道的入口流量xi和出口流量xo，观察管道的压力变化，将每个阶跃响应序列延迟一个时刻，用原序列减延迟序列获得四个脉冲响应序列；入口压力单位脉冲响应hii和hio，以及出口压力单位脉冲响应hoi和hoo；步骤三、管道压力脉冲响应仿真测试，包括：压力阶跃响应测试、单输入压力响应测试、双输入压力响应测试；步骤四、用脉冲响应实现在线动态仿真，包括单输入单输出离散系统的在线动态仿真和双输入双输出离散系统的在线动态仿真。进一步的，步骤三所述的管道压力脉冲响应仿真测试的具体方法为：第一步、压力阶跃响应测试，根据管道离散系统的压力输出计算公式分别计算压力对入口流量、出口流量的阶跃响应；第二步、单输入压力响应测试，在管道仿真软件中设计仿真实验，获取实验数据及管道参数，出口流量不变，获得管道入口压力和出口压力变化；对输气管道入口流量进行采样，根据获取的脉冲响应序列计算压力输出响应；然后,在仿真软件中对同一管道模型进行管道出口流量对压力影响的仿真实验，入口流量不变，获取管道入口压力和出口压力变化；对输气管道出口流量进行采样，根据脉冲响应序列计算压力输出响应；第三步、双输入压力响应测试，在管道仿真软件中对同一管道模型同时改变入口流量和出口流量，获取管道入口压力和出口压力变化；对输气管道入口、出口流量进行采样，根据脉冲响应分别计算入口压力、出口压力的输出响应，同时进行线性叠加作为系统的输出。进一步的，步骤四所述的在线动态仿真的具体方法为：第一步、单输入单输出离散系统的在线动态仿真，假设系统采样值xc(i)，系统相对稳态的输入输出为x(i)、y(i)，指定某一时刻s的输入输出为参考稳态，记相对参考稳态输入输出为x'(i)、y'(i)，则：x'(s)＝0x'(i)＝xc(i)-xs(xs)相对稳态的输出为：相对参考稳态的输出为：则实际工程单位下的仿真输出：yf(L+k+m)＝yf(L+k)+y'(L+k+m)-y'(L+k)-mh(L)(x0-xs)；第二步、双输入压力响应测试：Pi'(L+k+1)-Pi'(L+k)＝Pic(L+k+1)-Pic(L+k)+[(xio-xis)hii(L)+(xoo-xos)hio(L)]实际工程单位下的仿真输出Pif(L+k+m)为：Pif(L+k+m)＝Pic(L+k)+Pi'(L+k+m)-Pi'(L+k)-m[(xio-xis)hii(L)+(xoo-xos)hio(L)]式中：Pi(L+k)——L+k时刻相对稳态的管道入口压力；xi(n)——n时刻相对稳态的入口流量；xo(n)——n时刻相对稳态的出口流量；xic(n)——n时刻的入口流量采样；xoc(n)——n时刻的出口流量采样；xio——理想稳态的入口流量；xoo——理想稳态的出口流量；hii，hio——分别为入口压力对入口流量和出口流量的脉冲响应序列；Pi'(L+k)——L+k时刻相对参考稳态的管道入口压力；xi'(n)——n时刻相对参考稳态的入口流量；xo'(n)——n时刻相对参考稳态的出口流量；xis——s时刻的参考稳态相对稳态的入口流量；xos——s时刻的参考稳态相对稳态的出口流量；Pic(L+k)——L+k时刻管道入口压力实际工程单位的采样；Pif(L+k+m)——L+k+m时刻(未来的)实际工程单位下的仿真输出。本专利技术只需要通过管道机理模型计算一次管道的阶跃响应序列，通过延时计算可以得到各种采样间隔的脉冲响应序列，在进行动态仿真过程中只需要进行简单的卷积计算，计算量非常小，速度快，不存在计算发散的情况；通过设定参考稳态点，可以计算任意初始状态下的在线动态仿真输出。附图说明图1是本专利技术实施例提供的管道压力对入口流量的阶跃响应；图2是本专利技术实施例提供的管道压力对出口流量的阶跃响应；图3是本专利技术实施例提供的管道压力脉冲响应序列；图4是本专利技术实施例提供的入口流量阶跃变化，出口流量不变时压力对入口流量、出口流量的阶跃响应与实际阶跃响应数据的对比曲线；图5是本专利技术实施例提供的出口流量阶跃变化，入口流量不变时压力对入口流量、出口流量的阶跃响应与实际阶跃响应数据的对比曲线；图中Pi、Po是仿真结果，Pical、Pocal是理论计算结果；图6是本专利技术实施例提供的入口流量在200000M3/H附近变化，出口流量不变时获取的管道入口压力和出口压力变化曲线；图7是本专利技术实施例提供的入口流量在200000M3/H附近变化，出口流量不变时计算压力输出响应，与仿真输出数据的对比曲线；图中，Pi、Po是仿真结果，Pical、Pocal是理论计算结果；图8是本专利技术实施例提供的出口流量在200000M3/H附近变化，入口流量不变时获取的管道入口压力和出口压力变化曲线；图9是本专利技术实施例提供的出口流量在200000M3/H附近变化，入口流量不变时计算压力输出响应，与仿真输出数据的对比曲线；图10是本专利技术实施例提供的对同一管道模型同时改变入口流量和出口流量时获取的管道入口压力和出口压力变化曲线图；图11是本专利技术实施例提供的双输入双输出条件下系统的输出，与仿真输出数据的对比曲线；图12是本专利技术实施例提供的在线动态仿真的入口压力动态曲线；图13是本专利技术实施例提供的在线动态仿真的出口压力动态曲线；图14是本专利技术实施例提供的基于脉冲响应序列的输气管道动态仿真方法流程图。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本专利技术，并不用于限定本专利技术。本专利技术是这样实现的，一种基于脉冲响应序列的输气管道动态仿真方法包括：S101：建立管道脉冲响应序列模型，对于管道入口压力输出需要分别获取入口流量xi和出口流量xo为输入时的单位脉冲响应hii和hio，出口压力输出也需要分别获取入口流量和出口流量为输入时的单位脉冲响应hoi和hoo；在管道仿真软件中建立仿真对象，获取实验数据，管长50km，管道直径900cm，稳态入口压力8MPa，出口压力为6.39MPa，流量200000M3/H，将入口流量阶跃降低到180000M3/H，则入口压力迅速降低，然后指数衰减，经过一定的延时后出口压力缓慢下降，15分钟后，入口压力、出口压力几乎线性递减，类似的改变出口流量到220000本文档来自技高网...

【技术保护点】
一种基于脉冲响应模型的输气管道动态仿真方法，其特征在于，所述的基于脉冲响应模型的输气管道动态仿真方法包括：步骤一、获取管道压力脉冲响应序列，在管道模型中分别阶跃改变管道的入口流量xi和出口流量xo，观察管道的压力变化，将每个阶跃响应序列延迟一个时刻，用原序列减延迟序列获得四个脉冲响应序列；入口压力单位脉冲响应hii和hio，以及出口压力单位脉冲响应hoi和hoo；步骤二、建立管道脉冲响应序列模型，管道的入口压力Pi和出口压力Po可按下式计算：Pi(n)=Σk=1∞xi(k)hii(n-k)+Σk=1∞xo(k)hio(n-k)]]>Po(n)=Σk=1∞xi(k)hoi(n-k)+Σk=1∞xo(k)hoo(n-k)]]>步骤三、管道压力脉冲响应仿真测试，包括：压力阶跃响应测试、单输入压力响应测试、双输入压力响应测试；步骤四、用脉冲响应实现在线动态仿真，包括单输入单输出离散系统的在线动态仿真和双输入双输出离散系统的在线动态仿真。

【技术特征摘要】
1.一种基于脉冲响应模型的输气管道动态仿真方法，其特征在于，所述的基于脉冲响应模型的输气管道动态仿真方法包括：步骤一、建立管道脉冲响应序列模型，管道的入口压力Pi和出口压力Po按下式计算：步骤二、获取管道压力脉冲响应序列，在管道模型中分别阶跃改变管道的入口流量xi和出口流量xo，观察管道的压力变化，将每个阶跃响应序列延迟一个时刻，用原序列减延迟序列获得四个脉冲响应序列；入口压力单位脉冲响应hii和hio，以及出口压力单位脉冲响应hoi和hoo；步骤三、管道压力脉冲响应仿真测试，包括：压力阶跃响应测试、单输入压力响应测试、双输入压力响应测试；步骤四、用脉冲响应实现在线动态仿真，包括单输入单输出离散系统的在线动态仿真和双输入双输出离散系统的在线动态仿真。2.如权利要求1所述的基于脉冲响应模型的输气管道动态仿真方法，其特征在于，步骤三所述的管道压力脉冲响应仿真测试的具体方法为：第一步、压力阶跃响应测试，根据管道离散系统的压力输出计算公式分别计算压力对入口流量、出口流量的阶跃响应；第二步、单输入压力响应测试，在管道仿真软件中设计仿真实验，获取实验数据及管道参数，出口流量不变，获得管道入口压力和出口压力变化；对输气管道入口流量进行采样，根据获取的脉冲响应序列计算压力输出响应；然后,在仿真软件中对同一管道模型进行管道出口流量对压力影响的仿真实验，入口流量不变，获取管道入口压力和出口压力变化；对输气管道出口流量进行采样，根据脉冲响应序列计算压力输出响应；第三步、双输入压力响应测试，在管道仿真软件中对同一管道模型同时改变入口流量和出口流量，获取管道入口压力和出口压力变化；对输气管道入口、出口流量进行采样，根据脉冲响应分别计算入口压力、出口压力的输出响应，同时进行线性叠加作为系统的输出。3.如权利要求1所述的基于脉冲响应模型的输气管道动态仿真方法，其特征在于，步骤四所述的在线动态仿真的具体方法为：第一步、单输入单输出离散系统的在线动态仿真，假设系统采样值xc(i)，系统相对稳态的输入输出为x(i)、y(i)，指定某一时刻s的输入输出为参考稳态，记相对参考稳态输入输出为x'(i)、y'(i)，则：x'(s)＝0x'(i)＝xc(i)-xs(xs)相对稳态的输出为：相对参考稳态的输出为：

【专利技术属性】
技术研发人员：李立刚，戴永寿，
申请(专利权)人：中国石油大学华东，
类型：发明
国别省市：山东;37