CN103970029B - 一种基于脉冲响应序列的输气管道动态仿真方法 - Google Patents
一种基于脉冲响应序列的输气管道动态仿真方法 Download PDFInfo
- Publication number
- CN103970029B CN103970029B CN201410202027.9A CN201410202027A CN103970029B CN 103970029 B CN103970029 B CN 103970029B CN 201410202027 A CN201410202027 A CN 201410202027A CN 103970029 B CN103970029 B CN 103970029B
- Authority
- CN
- China
- Prior art keywords
- sigma
- pressure
- pipeline
- rate
- discharge
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于脉冲响应模型的输气管道动态仿真方法,首先建立管道脉冲响应序列模型;然后,获取管道压力脉冲响应序列,在管道模型中分别阶跃改变管道的入口流量和出口流量,观察管道的压力变化,将每个阶跃响应序列延迟一个时刻,用原序列减延迟序列获得四个脉冲响应序列;然后,进行管道压力脉冲响应仿真测试,包括:压力阶跃响应测试、单输入压力响应测试、双输入压力响应测试;最后,选择初始参考稳态点,根据管道历史采样数据,用脉冲响应实现在线动态仿真,包括单输入单输出离散***的在线动态仿真和双输入双输出离散***的在线动态仿真。
Description
技术领域
本发明属于动态仿真技术领域,尤其涉及一种基于脉冲响应序列的输气管道动态仿真方法。
背景技术
现有的管道仿真都是基于管道的连续性方程、运动方程、能量方程及气体状态方程等机理模型进行动态仿真计算,由于机理模型是一个多变量偏微分方程组,无法直接求解,需要转化为常微分方程组,并且对长输管道进行剖分,在时间和空间上进行离散化,计算过程需要迭代求解,计算量大,并且当管道初始条件设置不好时可能导致计算发散。
发明内容
本发明的目的在于提供一种基于脉冲响应序列的输气管道动态仿真方法,旨在解决利用管道机理模型计算量大、易发散的问题,提高管道动态仿真计算速度,实现管道在线动态仿真。
本发明是这样实现的,一种基于脉冲响应序列的输气管道动态仿真方法包括:
步骤一、建立管道脉冲响应序列模型,管道的入口压力Pi和出口压力Po按下式计算:
步骤二、获取管道压力脉冲响应序列,在管道模型中分别阶跃改变管道的入口流量xi和出口流量xo,观察管道的压力变化,将每个阶跃响应序列延迟一个时刻,用原序列减延迟序列获得四个脉冲响应序列;入口压力单位脉冲响应hii和hio,以及出口压力单位脉冲响应hoi和hoo;
步骤三、管道压力脉冲响应仿真测试,包括:压力阶跃响应测试、单输入压力响应测试、双输入压力响应测试;
步骤四、用脉冲响应实现在线动态仿真,包括单输入单输出离散***的在线动态仿真和双输入双输出离散***的在线动态仿真。
进一步的,步骤三所述的管道压力脉冲响应仿真测试的具体方法为:
第一步、压力阶跃响应测试,根据管道离散***的压力输出计算公式分别计算压力对入口流量、出口流量的阶跃响应;
第二步、单输入压力响应测试,在管道仿真软件中设计仿真实验,获取实验数据及管道参数,出口流量不变,获得管道入口压力和出口压力变化;对输气管道入口流量进行采样,根据获取的脉冲响应序列计算压力输出响应;然后,在仿真软件中对同一管道模型进行管道出口流量对压力影响的仿真实验,入口流量不变,获取管道入口压力和出口压力变化;对输气管道出口流量进行采样,根据脉冲响应序列计算压力输出响应;
第三步、双输入压力响应测试,在管道仿真软件中对同一管道模型同时改变入口流量和出口流量,获取管道入口压力和出口压力变化;对输气管道入口、出口流量进行采样,根据脉冲响应分别计算入口压力、出口压力的输出响应,同时进行线性叠加作为***的输出。
进一步的,步骤四所述的在线动态仿真的具体方法为:
第一步、单输入单输出离散***的在线动态仿真,假设***采样值xc(i),***相对稳态的输入输出为x(i)、y(i),指定某一时刻s的输入输出为参考稳态,记相对参考稳态输入输出为x'(i)、y'(i),则:
x'(s)=0
x'(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分钟后,入口压力、出口压力几乎线性递减,类似的改变出口流量到220000M3/H,15分钟后,入口压力、出口压力几乎线性递减,设定采样间隔为18秒,采样100个点即采样30分钟内的阶跃响应,这样可以获得较高的精度,绘制阶跃响应曲线如图1和图2。
S102:获取管道压力脉冲响应序列,在管道模型中分别阶跃改变管道的入口流量和出口流量,观察管道的压力变化,将每个阶跃响应序列延迟一个时刻,用原序列减延迟序列获得四个脉冲响应序列,如图3所示;
S103:管道压力脉冲响应仿真测试,包括:
第一步、压力阶跃响应测试:
根据管道离散***的压力输出计算公式分别计算压力对入口流量、出口流量的阶跃响应,如图4和图5所示,其中Pi、Po是TGNET仿真结果,Pical、Pocal是理论计算结果,二者几乎重合;
第二步、单输入压力响应测试:
在管道仿真软件中设计仿真实验,获取实验数据及管道参数,管道参数同步骤S101,入口流量在200000M3/H附近变化,出口流量不变,管道入口压力和出口压力变化如图6所示;
对输气管道入口流量进行采样,根据获取的脉冲响应序列计算压力输出响应,与仿真输出数据的对比曲线如图7所示,入口压力的输出响应与仿真数据的最大偏差0.0042MPa,出口压力的输出响应与仿真数据的最大偏差0.0001MPa;
然后,在仿真软件中对同一管道模型进行管道出口流量对压力影响的仿真实验,出口流量在200000M3/H附近变化,入口流量不变,获取管道入口压力和出口压力变化,如图8所示;
对输气管道出口流量进行采样,根据脉冲响应序列计算压力输出响应,与仿真输出数据的对比曲线如图9所示,对变化的出口流量,入口压力的输出响应与仿真数据的最大压力偏差为0.0026MPa,出口压力的输出响应与仿真数据的最大压力偏差为0.0166MPa;
第三步、双输入压力响应测试:
在管道仿真软件中对同一管道模型同时改变入口流量和出口流量,获取管道入口压力和出口压力变化,如图10所示;
对输气管道入口、出口流量进行采样,根据脉冲响应分别计算入口压力、出口压力的输出响应,同时进行线性叠加作为***的输出,与仿真输出数据的对比曲线如图11所示,入口压力的输出响应与仿真数据的最大压力偏差为0.005MPa,出口压力的输出响应与仿真数据的最大压力偏差为0.0153MPa,出现在第31点即0.155小时,此时流量变化剧烈采样间隔相对较大,在流量变化较为缓慢的时间里最大偏差一般小于0.003MPa,足以满足工业现场的精度要求。
S104:用脉冲响应实现在线动态仿真,包括单输入单输出离散***的在线动态仿真和双输入双输出离散***的在线动态仿真;
第一步、单输入单输出离散***的在线动态仿真,假设***相对稳态的输入输出为x(i)、y(i),***参数采样值为xc(i),如果指定某一时刻s的输入输出为参考稳态xs,记相对参考稳态输入输出为x'(i)、y'(i),则:
x'(s)=0
x'(i)=xc(i)-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时刻(未来的)实际工程单位下的仿真输出。
在L+k时刻开始仿真计算,可以将所有的仿真输出Pi'(L+k+m)减去Pi'(L+k),同时以L+k时刻的采样输出Pic(L+k)为基准进行叠加,计算结果将是未来实际工程单位下的仿真输出Pif(L+k+m)。
以半小时数据动态仿真为例,管道仿真数据和脉冲响应序列同前,在k=120、130、140点开始进行仿真测试,分别利用上述方法计算实际工程单位下的仿真输出绘制动态仿真曲线和管道动态数据曲线,入口压力动态仿真曲线如图12所示,出口压力动态仿真曲线如图13所示,入口压力最大误差分别为0.0037MPa、0.0012MPa、0.0003MPa,出口压力最大误差分别为0.0019MPa、0.0013MPa、0.0022Mpa。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性的劳动即可做出的各种修改或变形仍在本发明的保护范围之内。
Claims (3)
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)=0
x'(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时刻实际工程单位下的仿真输出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410202027.9A CN103970029B (zh) | 2014-05-14 | 2014-05-14 | 一种基于脉冲响应序列的输气管道动态仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410202027.9A CN103970029B (zh) | 2014-05-14 | 2014-05-14 | 一种基于脉冲响应序列的输气管道动态仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103970029A CN103970029A (zh) | 2014-08-06 |
CN103970029B true CN103970029B (zh) | 2017-03-29 |
Family
ID=51239671
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410202027.9A Expired - Fee Related CN103970029B (zh) | 2014-05-14 | 2014-05-14 | 一种基于脉冲响应序列的输气管道动态仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103970029B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109697308B (zh) * | 2018-11-30 | 2023-04-07 | 天津大学 | 一种考虑管网延时特性的天然气输气***稳态建模方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4559610A (en) * | 1983-05-04 | 1985-12-17 | Southwest Research Corporation | Gas pumping system analog |
CN102779198A (zh) * | 2011-05-13 | 2012-11-14 | 中国石油天然气股份有限公司 | 天然气管网输送***动态仿真***建立方法 |
CN102779199A (zh) * | 2011-05-13 | 2012-11-14 | 中国石油天然气股份有限公司 | 天然气管网输送***动态仿真过程的实现方法 |
-
2014
- 2014-05-14 CN CN201410202027.9A patent/CN103970029B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4559610A (en) * | 1983-05-04 | 1985-12-17 | Southwest Research Corporation | Gas pumping system analog |
CN102779198A (zh) * | 2011-05-13 | 2012-11-14 | 中国石油天然气股份有限公司 | 天然气管网输送***动态仿真***建立方法 |
CN102779199A (zh) * | 2011-05-13 | 2012-11-14 | 中国石油天然气股份有限公司 | 天然气管网输送***动态仿真过程的实现方法 |
Non-Patent Citations (2)
Title |
---|
煤层气输送管道仿真技术研究;李长俊 等;《2010年中国天然气产业及***管网建设发展论坛》;20101231;132-139 * |
输气管道仿真教学***的研发;王明阳 等;《实验室研究与探索》;20060630;第25卷(第6期);650-652,656 * |
Also Published As
Publication number | Publication date |
---|---|
CN103970029A (zh) | 2014-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108304684B (zh) | 一种火箭发动机尾喷射流仿真方法及*** | |
US20150261893A1 (en) | Method and apparatus for determining pipeline flow status parameter of natural gas pipeline network | |
Legoll et al. | A micro-macro parareal algorithm: application to singularly perturbed ordinary differential equations | |
CN104408744A (zh) | 一种用于目标跟踪的强跟踪容积卡尔曼滤波方法 | |
CN102435844A (zh) | 一种频率无关的正弦信号相量计算方法 | |
CN104019817A (zh) | 一种用于卫星姿态估计的范数约束强跟踪容积卡尔曼滤波方法 | |
CN103728882A (zh) | 微陀螺仪的自适应反演非奇异终端滑模控制方法 | |
CN107066712A (zh) | 基于导叶开度—功率特性的水轮机模型参数辨识方法及装置 | |
WO2014181183A9 (en) | System and method for multi-phase fluid measurement | |
Zhang et al. | Modeling and simulation of time-varying delays | |
CN104199307B (zh) | 硬件在环仿真方法和*** | |
CN105425589B (zh) | 提高航天器惯性参数辨识精度的输入信号设计方法 | |
CN108459087A (zh) | 一种应用于板结构损伤检测的多模态Lamb波模态分离方法 | |
CN103218482A (zh) | 一种动力学***中不确定参数的估计方法 | |
Fetisov et al. | Aggregative simulation method for implementing mathematical models for gas transmission systems | |
CN105989205B (zh) | 飞行器表面脉动压力的确定方法 | |
CN103970029B (zh) | 一种基于脉冲响应序列的输气管道动态仿真方法 | |
CN101479613B (zh) | 瞬时确定交流电网上的信号畸变率的方法和相关装置 | |
CN102654773B (zh) | 一种基于zvdd和pwm混合输入成型器的挠性航天器控制方法 | |
CN106292289B (zh) | 流程工业控制回路对象的混合精英随机搜索优化方法 | |
DE502007006897D1 (de) | Verfahren und vorrichtung zur identifikation einer verzögerungsbehafteten regelstrecke, regeleinrichtung und computerprogrammprodukt | |
EP2738630B1 (en) | Method for high fidelity modeling of an aircraft electrical power system | |
EP2887227A2 (en) | Information processing apparatus and method | |
US11610034B2 (en) | Method and device for synchronizing a simulation with a real-time system | |
CN109858170A (zh) | 一种强非线性Modelica***模型的频率特性估算***及方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170329 Termination date: 20180514 |
|
CF01 | Termination of patent right due to non-payment of annual fee |