CN113531395B - 一种天然气管网的购销气差率实时监测方法 - Google Patents

一种天然气管网的购销气差率实时监测方法 Download PDF

Info

Publication number
CN113531395B
CN113531395B CN202110727541.4A CN202110727541A CN113531395B CN 113531395 B CN113531395 B CN 113531395B CN 202110727541 A CN202110727541 A CN 202110727541A CN 113531395 B CN113531395 B CN 113531395B
Authority
CN
China
Prior art keywords
natural gas
gas
real
pipe
time
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.)
Active
Application number
CN202110727541.4A
Other languages
English (en)
Other versions
CN113531395A (zh
Inventor
尹祥
郑权明
李建勋
罗东晓
刘华
姚祉百
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Foshan Gas High Pressure Pipe Network Co ltd
Guangdong Foran Technology Co ltd
Fo Ran Energy Group Co ltd
Original Assignee
Foshan Gas High Pressure Pipe Network Co ltd
Guangdong Foran Technology Co ltd
Fo Ran Energy Group Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Foshan Gas High Pressure Pipe Network Co ltd, Guangdong Foran Technology Co ltd, Fo Ran Energy Group Co ltd filed Critical Foshan Gas High Pressure Pipe Network Co ltd
Priority to CN202110727541.4A priority Critical patent/CN113531395B/zh
Publication of CN113531395A publication Critical patent/CN113531395A/zh
Application granted granted Critical
Publication of CN113531395B publication Critical patent/CN113531395B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17D—PIPE-LINE SYSTEMS; PIPE-LINES
    • F17D5/00—Protection or supervision of installations
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17D—PIPE-LINE SYSTEMS; PIPE-LINES
    • F17D1/00—Pipe-line systems
    • F17D1/02—Pipe-line systems for gases or vapours
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17D—PIPE-LINE SYSTEMS; PIPE-LINES
    • F17D3/00—Arrangements for supervising or controlling working operations
    • F17D3/01—Arrangements for supervising or controlling working operations for controlling, signalling, or supervising the conveyance of a product

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Pipeline Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供了一种天然气管网的购销气差率实时监测方法，属于天然气管网技术领域。本发明将整个天然气管网管道划分成若干管段，通过折算系数解决每一个管段及整个天然气管网管道内有效空间体积的获取问题，根据每一个管段内的天然气压力、温度、组分数据，按照相关技术标准及规范得出压缩因子数据，获得充分考虑介质压力、温度、组分及压缩因子影响之后的各管段内天然气的实时储气存量及整个天然气管网管道的实时储气量；按照本发明给出的计算方法，获取预设时间段内整个天然气管网天然气购销气差率。本发明管段划分方法科学，解决了有效空间体积及实时储气量获取的问题，可实时获取准确的天然气管网的购销气差率。

Description

一种天然气管网的购销气差率实时监测方法
技术领域
本发明涉及天然气管网技术领域，尤其涉及一种天然气管网的购销气差率实时监测方法。
背景技术
一个城镇燃气完整的天然气输配管网，一般都涵盖多个天然气门站、调压站，若干长度的天然气输送管道，若干配置有计量仪表的气源输入点，若干配置有计量仪表的气源输出点，以及使用天然气的用户。整个管网的天然气介质流动实时处于有序平衡状态，天然气气源输入与输出时刻处于动态平衡状态。
作为燃气企业，购销气差率是衡量其运营管理水平的一个重要指标，从精细化管理角度考虑，企业迫切需要实时掌握整个天然气管网准确的购销气差率，以便及时发现并解决问题，提高管理水平，提升企业经济效益。
但目前的现实情况是，国内大部分燃气企业一般是采用人工统计核算方式，通过计算某一个时段(通常是间隔一个月甚至更长时段)的气源输入总量与气源输出总量之间的天然气气量差值，在大致核算管网的储气量变化的基础上，得出整个天然气管网该时段的购销气差率。
这种方式存在的明显问题：一是不能实时掌握动态的购销气差率数据；二是因为输配管网系统复杂，管道长，管网各处管道口径、介质压力与温度均不相同，无法准确计算出动态的管道内天然气储气量，采用这种传统的统计方式得出的购销气差率必然准确度不高，特别是时段不够长的情况；三是效率低下，需要花费一定的时间才能大致得出购销气差率数据结果。
可见，这些现有技术至少存在以下不足：
1.不能实时掌握动态的购销气差率数据。
2.无法准确计算出动态的管网管道内天然气储气量，采用这种传统的统计方式得出的购销气差率必然准确度不高。
发明内容
为解决现有技术中存在的技术问题，本发明提供了一种天然气管网的购销气差率实时监测方法，将整个天然气管网管道划分成若干管段，通过折算系数解决每一个管段及整个天然气管网管道内有效空间体积的获取问题，同时解决每一个管段的天然气实时压力、温度、组分的获取问题；根据每一个管段内的天然气压力、温度、组分数据，按照相关技术标准及规范得出压缩因子数据，获得充分考虑管道污垢导致有效空间体积变小、天然气介质压力、温度、组分及压缩因子影响之后的各管段管道内天然气的实时储气存量及整个天然气管网管道的实时储气量；按照本发明给出的计算方法，获取预设时间段内整个天然气管网天然气购销气差率。本发明管段划分方法科学、合理，具有经济性和可操作性，各管段管道受污垢影响条件下的有效空间体积获取方法既可行又方便，圆满解决了各管段管道内天然气实时压力、温度、组分及压缩因子数据获取问题，可实时获取准确的天然气管网的购销气差率。
本发明提供了一种天然气管网的购销气差率实时监测方法，包括如下步骤：
根据管道口径变化、管道长度及传感器安装位置选择难易度，对整个天然气管网进行管段划分；
获得各管段及整个天然气管网的管道内的有效空间体积；
获取每一个管段天然气实时压力和实时温度；
获取每一个管段天然气实时组分；
根据压力、温度和组分获取每一个管段天然气压缩因子；
获取整个天然气管网总储气量；
获取预设时间段内整个天然气管网的天然气储气量的变化量；
获取天然气管网购销气差率。
优选地，根据管道口径变化、管道长度及传感器安装位置选择难易度，进行管段划分包括：
相邻的不同管径的管道，每个管径的管道划分为一个管段；
相邻的管径相同的管道，长度不超过10公里时，划分为一个管段；
相邻的管径相同的管道，长度超过10公里时，划分为2个或多个管段。
划分后的每一个管段内天然气压力、温度相同或连续均衡变化。
优选地，获取各管段的管道内的有效空间体积具体包括：依据各管段的管道的运行年限、材质及运行状态，考虑各管段的管道内壁腐蚀及污垢情况，确定有效空间体积折算系数，根据有效空间体积折算系数和管段的管道内体积确定有效空间体积。
优选地，有效空间体积折算系数具体为：运行时间在5年以内管段的有效空间体积折算系数为1.0；运行时间超过20年且天然气介质含水量大于0.5mg/m3的有效空间体积折算系数为0.9998；其他情况，有效空间体积折算系数在0.9998-1.0之间采用比例分布法推算。
优选地，每隔20公里以上，在管道管径改变处或线路切断阀门位置设置压力监测点，在压力监测点采集压力值，并据此获取各管段的天然气实时压力。
优选地，根据相邻两个压力监测点采集的压力值和管道长度，采用比例分布法获取各管段首尾两端的天然气实时压力，再将各管段首尾两端的天然气实时压力取平均值，获取每一个管段的天然气实时压力，其中比例分布法采用如下公式：
公式如下：
Figure BDA0003139202960000031
其中，
P0为相邻两个压力监测点中距离气源更近的压力监测点处的实时压力；
P1为待求管段端点的实时压力，是一个管段尾端压力，同时也是下一个管段首端压力；
P2为相邻两个压力监测点中距离气源更远的压力监测点处的实时压力；
S1为待求实时压力管段与离气源更近的压力监测点之间的管道长度；
S为相邻两个压力监测点之间的管道长度。
优选地，依据各管段首尾两端的天然气实时压力，采用平均法获取该管段天然气实时压力，公式如下：
PL＝(P0+P1)/2
PL为L管段的天然气实时压力；
P0为L管段距离气源更近的一端(首端)的天然气实时压力；
P1为L管段另外一端(尾端)的天然气实时压力；
优选地，温度监测点设置位置与压力监测点设置位置相同。
优选地，根据相邻两个温度监测点采集的温度值和管道长度，采用比例分布法进行推算。
优选地，获取每一个管段天然气实时组分包括：
对于单一气源管网，在气源输入处设置在线气相色谱分析仪，气相色谱分析仪实时分析数据作为各管段天然气组分；
对于多气源管网，在每个气源输入处设置在线气相色谱分析仪，实时分析天然气组分，根据每一管段管道中各气源输入量的比例，计算出每一个管段的实时天然气组分。
优选地，根据各管段管道内天然气的实时压力、温度、组分、有效空间体积和压缩因子，获得每一管段内天然气实时储气量。
优选地，每一管段内天然气实时储气量采用如下公式计算：
V＝ZmRT/P
V：天然气储气量，单位为m3；
Z：天然气压缩因子，无量纲；
m：气体质量，单位为kg；
R：摩尔气体常数，单位为J/kg﹒K；
T：天然气温度，单位为K；
P：天然气绝对压力：单位为Pa。
优选地，整个管网天然气购销气差率采用如下公式计算：
Figure BDA0003139202960000041
其中：
R：天然气购销气差率；
A：预设时间段内管网的天然气输入气量，单位为m3；
B：预设时间段内管网的天然气储气量变化量，单位为m3；
C：预设时间段内管网的天然气输出气量，单位为m3。
与现有技术相对比，本发明的有益效果如下：
1.本发明全面考虑管道管径变化、同口径管道的长度、管径变化点位置、既有线路切断阀门位置、压力变化曲线、传感器安装难易度、投入成本等因素，通过全面统筹，采用科学合理方法将整个管网划分为若干管段，提出了管段划分的原则和具体技术方案；
2.本发明充分考虑管道内壁腐蚀以及管道内污垢对管道空间体积的影响，通过对历史数据进行统计分析和进行大量的试验测试，结合工程实际经验，提出采用折算系数方式对管道有效空间体积进行修正，并给出了折算系数确定原则和具体技术标准，解决了获取管道内有效空间体积的难题；
3.本发明通过设置有限的压力、温度及组分分析监测点，并未采用每一个管段均设置压力、温度及组分分析监测点，实现了实时获得每一管段天然气压力、温度、组分及压缩因子数据的目标，投入少，且能满足需求；
4.本发明实时采集天然气压力、温度、组分及计量仪表数据，并根据计算得出的每一管段天然气实时压力、温度、组分及压缩因子数据，获得天然气管网实时储气量及预设时间段内管网实时购销气差率，计算过程符合国家相关标准。
附图说明
图1是本发明的一个实施例的管段划分示意图；
图2是本发明的一个实施例的压力监测点设置方法示意图；
图3是本发明的又一个实施例的单一气源管段划分示意图；
图4是本发明的又一个实施例的多气源管段划分示意图；
图5是本发明的一个实施例的天然气管网的购销气差率实时监测方法示意图；
图6是本发明的一个实施例的天然气管网的购销气差率实时监测方法流程图；
图7是本发明的又一个实施例的天然气管网的购销气差率实时监测方法流程图。
具体实施方式
下面结合附图，对本发明的具体实施方式作详细的说明。
本发明提供了一种天然气管网的购销气差率实时监测方法，包括如下步骤：
根据管道口径变化、管道长度及传感器安装位置选择难易度，对整个天然气管网进行管段划分；
获得各管段及整个天然气管网的管道内的有效空间体积；
获取每一个管段天然气实时压力和实时温度；
获取每一个管段天然气实时组分；
根据压力、温度和组分获取每一个管段天然气压缩因子；
获取整个天然气管网总储气量；
获取预设时间段内整个天然气管网的天然气储气量的变化量；
获取天然气管网购销气差率。
根据本发明的一个具体实施方案，根据管道口径变化、管道长度及传感器安装位置选择难易度，进行管段划分包括：
相邻的不同管径的管道，每个管径的管道划分为一个管段；
相邻的管径相同的管道，长度不超过10公里时，划分为一个管段；
相邻的管径相同的管道，长度超过10公里时，划分为2个或多个管段。
划分后的每一个管段内天然气压力、温度相同或连续均衡变化。
根据本发明的一个具体实施方案，获取各管段的管道内的有效空间体积具体包括：依据各管段的管道的运行年限、材质及运行状态，考虑各管段的管道内壁腐蚀及污垢情况，确定有效空间体积折算系数，根据有效空间体积折算系数和管段的管道内体积确定有效空间体积。
根据本发明的一个具体实施方案，有效空间体积折算系数具体为：运行时间在5年以内管段的有效空间体积折算系数为1.0；运行时间超过20年且天然气介质含水量大于0.5mg/m3的有效空间体积折算系数为0.9998；其他情况，有效空间体积折算系数在0.9998-1.0之间采用比例分布法推算。以5年(1.0)和20年(0.9998)为两个临界点，其他情况下有效空间体积折算系数计算采用如下公式进行推算：
C＝1-[(N-5)/(20-5)]*(1-0.9998)＝1-0.0001*(N-5)/2；
其中，
C为折算系数，无量纲；
N为管段运行时间，单位：年；
根据本发明的一个具体实施方案，每隔20公里以上，在管道管径改变处或线路切断阀门位置设置压力监测点，在压力监测点采集压力值，并据此获取各管段的天然气实时压力。
根据本发明的一个具体实施方案，温度监测点设置位置与压力监测点设置位置相同。
根据本发明的一个具体实施方案，根据相邻两个温度监测点采集的温度值和管道长度，采用比例分布法进行推算。不一定离气源近温度就低。由于温度是均衡变化的，通过采集的两个温度监测点温度值，就可以采用比例分布法推算这两个监测点之间某一位置的温度值。
根据本发明的一个具体实施方案，根据相邻两个压力监测点采集的压力值和管道长度，采用比例分布法获取各管段首尾两端的天然气实时压力，再将各管段首尾两端的天然气实时压力取平均值，获取每一个管段的天然气实时压力，其中比例分布法采用如下公式：
Figure BDA0003139202960000061
其中，
P0为相邻两个压力监测点中距离气源更近的压力监测点处的实时压力；
P1为待求管段端点的实时压力，是一个管段尾端压力，同时也是下一个管段首端压力；
P2为相邻两个压力监测点中距离气源更远的压力监测点处的实时压力；
S1为待求实时压力管段与离气源更近的压力监测点之间的管道长度；
S为相邻两个压力监测点之间的管道长度。
根据本发明的一个具体实施方案，依据各管段首尾两端的天然气实时压力，采用平均法获取该管段天然气实时压力，公式如下：
PL＝(P0+P1)/2
PL为L管段的天然气实时压力；
P0为L管段距离气源更近的一端(首端)的天然气实时压力；
P1为L管段另外一端(尾端)的天然气实时压力；
根据本发明的一个具体实施方案，获取每一个管段天然气实时组分包括：
对于单一气源管网，在气源输入处设置在线气相色谱分析仪，气相色谱分析仪实时分析数据作为各管段天然气组分；
对于多气源管网，在每个气源输入处设置在线气相色谱分析仪，实时分析天然气组分，根据每一管段管道中各气源输入量的比例，计算出每一个管段的实时天然气组分。
根据本发明的一个具体实施方案，根据各管段管道内天然气的实时压力、温度、组分、有效空间体积和压缩因子，获得每一管段内天然气实时储气量。
根据本发明的一个具体实施方案，每一管段内天然气实时储气量采用如下公式计算：
V＝ZmRT/P
V：天然气储气量，单位为m3；
Z：天然气压缩因子，无量纲；
m：气体质量，单位为kg；
R：摩尔气体常数，单位为J/kg﹒K；
T：天然气温度，单位为K；
P：天然气绝对压力：单位为Pa。
根据本发明的一个具体实施方案，整个管网天然气购销气差率采用如下公式计算：
Figure BDA0003139202960000071
其中：
R：天然气购销气差率；
A：预设时间段内管网的天然气输入气量，单位为m3；
B：预设时间段内管网的天然气储气量变化量，单位为m3；
C：预设时间段内管网的天然气输出气量，单位为m3。
实施例1
参照附图1，根据本发明的一个具体实施方案，对本发明提供的管段划分的方法进行详细说明。
比如有一段总长度26.9km的天然气管道，中间无调压站或大用户，其中DN800管道有两段，总长度为18km，L1段长度为9km，L2段长度为9km，之后管道管径发生变化，为DN600管道，L3长度为8.9km。
按照如下划分原则，将管道划分为3个管段。划分原则为：
划分后的每一个管段内天然气压力连续均衡变化；
相邻的不同管径的管道，每个管径的管道划分为一个管段；
相邻的管径相同的管道，长度不超过10公里时，划分为一个管段；
相邻的管径相同的管道，长度超过10公里时，划分为2个或多个管段。
即将管径相同但长度较长的18km、DN800管道划分为2个管段L1、L2，长度各为9km；相邻的管径为DN600管道划分为一个管段L3，长度为8.9km。合计3各管段。
实施例2
参照附图2，根据本发明的一个具体实施方案，对本发明提供的压力监测点设置方法进行详细说明。
长度为38公里的DN1000管道划分为4个管段，8.9公里的DN600管道划分为1个管段。
在长度为38公里的DN1000管道与长度为8.9公里的DN600管道连接处(变径处)设置1个压力监测点(P2)；另外在DN1000与DN600管道另外一端各设置1个压力监测点(P1、P3),合计3个压力监测点。只设置3个压力监测点，就可以实时获得5个管段的压力。
具体地：L5管段的实时压力为P2、P3压力监测点采集压力的平均值；L1、L2、L3、L4管段的实时压力则为P1、P2压力监测点采集压力依据管道长度、按照比例分布法分别计算得出每一管段首尾两端的压力，公式如下：
Figure BDA0003139202960000081
其中，
P0为相邻两个压力监测点中距离气源更近的压力监测点处的实时压力；
P1为待求管段端点的实时压力，是一个管段尾端压力，同时也是下一个管段首端压力；
P2为相邻两个压力监测点中距离气源更远的压力监测点处的实时压力；
S1为待求实时压力管段与离气源更近的压力监测点之间的管道长度；
S为相邻两个压力监测点之间的管道长度。
然后，取每一管段首尾两端的实时压力平均值作为该管段实时天然气压力。
在本实施例中，在设置了压力监测点的位置设置温度监测点，各管段实时温度获取方法与压力类似，采用比例分布法获取。
实施例3
参照附图3，根据本发明的一个具体实施方案，为单一气源管网的实施方案，对本发明提供的复杂管网管段划分、压力监测点设置、有效空间体积、管网储气量及管网购销气差率获取进行详细说明。
在图3中，主干管网长度52.5km，管径分别为DN1000、DN800、DN600、DN500、DN400；支线管网9.6km，DN500；一个气源输入点(计量仪表)；两个调压站；一个大客户；3个气源输出点(计量仪表)；管道运行年限8年，天然气介质含水量大于0.5mg/m3。
1、管段划分
顺着气源流动方向，按照本发明的管段划分原则，将整个管网划分为L1～L8共8个管段。划分原则：一般情况下，将不同管径管道分别划分为独立管段，将管径相同但长度较长的管道可划分为多个管段，一个管段长度原则上不超过10km。管段划分越短，计算出来的管网储气量就越准确，但管段越短投资越大，需要统筹考虑。
2、压力监测点的设置与实时压力获取
原则上尽可能减少压力监测点设置数量，以此减少投入。一般情况下，2个压力监测点之间的管道长度必须超过20km；压力相同或连续均衡变化，且能够容易推算出实时压力的管道一般不在中间位置设置压力监测点，依靠相邻压力监测点实时压力推断出来。
按照上述原则，本实施例中，门站及2个调压站、大用户处均本身设置有压力监测点，因此，只需在大用户支线管道与主管道交叉处设置压力监测点即可，即在L3、L4、L5三个管段交汇处设置一个压力监测点，共5个压力监测点(A、B、C、D和E处)。所有管段压力均依据压力检测点采集的实时压力，采用比例分布法进行推算首先得出每一管段首尾两端的压力，公式如下：
Figure BDA0003139202960000091
其中，
P0为相邻两个压力监测点中距离气源更近的压力监测点处的实时压力；
P1为待求管段端点的实时压力，是一个管段尾端压力，同时也是下一个管段首端压力；
P2为相邻两个压力监测点中距离气源更远的压力监测点处的实时压力；
S1为待求实时压力管段与离气源更近的压力监测点之间的管道长度；
S为相邻两个压力监测点之间的管道长度。
然后，取每一管段首尾两端的实时压力平均值作为该管段实时天然气压力。
3、各管段有效空间体积的获取
1)依据各管段管道的长度、管径，按照圆柱体体积计算公式计算出各管段管道内的空间体积；
2)确定各管段折算系数
因为运行年限等条件完全相同，所以8个管段折算系数相同。
按照本发明的折算系数确定原则，折算系数＝1-[(8-5)/(20-5)]×(1-0.9998)＝0.99996；
3)以此折算系数乘以各管段管道的空间体积，即为各管段有效空间体积；
4)将各管段管道有效空间体积合计相加，得出整个管网总的有效空间体积。
4、各管段实时天然气组分
在气源输入处(门站)设置在线气相色谱分析仪，因为是单一气源输入，流经每一个管段内的天然气组分完全相同，因此，采集该分析仪分析结果的组分即为各管段实时天然气组分。
5、各管段压缩因子计算
依据国标《天然气压缩因子的计算》(GB T 17747.2—2011)，用各管段管道内天然气摩尔组成进行计算，获得天然气压缩因子数据；
进而按照如下公式计算获得各管段管道内天然气的实时储存量，并将各管段管道内天然气的实时储存量合计为整个管网的总储气量：
V＝ZmRT/P
V：天然气储气量，单位为m3；
Z：天然气压缩因子，无量纲；
m：气体质量，单位为kg；
R：摩尔气体常数，单位为J/kg﹒K；
T：天然气温度，单位为K；
P：天然气绝对压力：单位为Pa。
6、获取预设时间段内整个天然气管网天然气购销气差率
输入预设时间段(起始时间)，系统采集气源输入与气源输出点计量仪表数据，获得预设时间段内管网的天然气输入气量和天然气输出气量。按照如下计算公式，获取预设时间段内整个天然气管网天然气购销气差率：
Figure BDA0003139202960000111
其中：
R：天然气购销气差率；
A：预设时间段内管网的天然气输入气量，单位为m3；
B：预设时间段内管网的天然气储气量变化量，单位为m3；
C：预设时间段内管网的天然气输出气量，单位为m3。
实施例4
根据本发明的一个具体实施方案，参照附图4，为多气源管网的实施方案，对本发明提供的天然气管网的购销气差率实时监测方法进行详细说明。
图4中，主干管网长度55.1公里，管径分别为DN1000、DN800、DN600、DN500；支线管网21.3km，管径为DN500；两个门站2个气源输入点(计量仪表)；两个调压站；1个大客户；3个气源输出点(计量仪表)；管道运行年限8年，天然气介质含水量大于0.5mg/m3。
1、管段划分
按照本发明的管段划分原则，将整个管网划分为L1～L10共10个管段。
2、压力监测点的设置与实时压力获取
压力监测点设置原则以及实时压力获取的方法与实施例3相同。
本实施例中，2个门站、2个调压站、大用户处均本身设置有压力监测点，因此，只需在大用户支线管道与主管道交叉处设置压力监测点即可，即在L3、L4和L5三个管段交汇处及L7、L8和L9三个管段交汇处设置2个压力监测点，共7个压力监测点(A、B、C、D、E、F和G处)。所有管段压力均依据压力检测点采集的实时压力，采用比例分布法进行推算首先得出每一管段首尾两端的压力，公式如下：
Figure BDA0003139202960000112
其中，
P0为相邻两个压力监测点中距离气源更近的压力监测点处的实时压力；
P1为待求管段端点的实时压力，是一个管段尾端压力，同时也是下一个管段首端压力；
P2为相邻两个压力监测点中距离气源更远的压力监测点处的实时压力；
S1为待求实时压力管段与离气源更近的压力监测点之间的管道长度；
S为相邻两个压力监测点之间的管道长度。
然后，取每一管段首尾两端的实时压力平均值作为该管段实时天然气压力。
3、各管段有效空间体积获取
各管段有效空间体积的获取方法与实施例3相同。
4、各管段实施天然气组分
在气源输入1、气源输入2处设置在线气相色谱分析仪。
L1～L7管段内全部流量均是气源输入1的天然气；L8管段内全部流量均是气源输入2的天然气；L9～L10管段内的气体中，气源输入1的天然气占30％，气源输入2的天然气占70％.
根据本发明的内容：
L1～L7管段内的实时天然气组分就是气源输入1的在线气相色谱分析仪分析结果；
L8管段内的实时天然气组分就是气源输入2的在线气相色谱分析仪分析结果；
L9～L10管段内的实时天然气组分：依据气源输入1和气源输入2的在线气相色谱分析仪分析结果，按照30％与70％比例采用加权平均法得出。
5、各管段压缩因子计算
依据国标《天然气压缩因子的计算》(GB T 17747.2—2011)，用各管段管道内天然气摩尔组成进行计算，获得天然气压缩因子数据；
进而按照如下公式计算获得各管段管道内天然气的实时储存量，并将各管段管道内天然气的实时储存量合计为整个管网的总储气量：
V＝ZmRT/P
V：天然气储气量，单位为m3；
Z：天然气压缩因子，无量纲；
m：气体质量，单位为kg；
R：摩尔气体常数，单位为J/kg﹒K；
T：天然气温度，单位为K；
P：天然气绝对压力：单位为Pa。
6、获取预设时间段内整个天然气管网天然气购销气差率
输入预设时间段(起始时间)，系统采集气源输入与气源输出点计量仪表数据，获得预设时间段内管网的天然气输入气量和天然气输出气量。按照如下公式，获取预设时间段内整个天然气管网天然气购销气差率
Figure BDA0003139202960000131
其中：
R：天然气购销气差率；
A：预设时间段内管网的天然气输入气量，单位为m3；
B：预设时间段内管网的天然气储气量变化量，单位为m3；
C：预设时间段内管网的天然气输出气量，单位为m3。
实施例5
根据本发明的一个具体实施方案，参照附图5-7，对本发明提供的天然气管网的购销气差率实时监测方法进行详细说明。
本发明提供了一种天然气管网的购销气差率实时监测方法，包括如下步骤：
根据管道口径变化、管道长度及传感器安装位置选择难易度，对整个天然气管网进行管段划分，具体包括：
相邻的不同管径的管道，每个管径的管道划分为一个管段；
相邻的管径相同的管道，长度不超过10公里时，划分为一个管段；
相邻的管径相同的管道，长度超过10公里时，划分为2个或多个管段。
划分后的每一个管段内天然气压力、温度相同或连续均衡变化。
获得各管段及整个天然气管网的管道内的有效空间体积，具体包括：依据各管段的管道的运行年限、材质及运行状态，考虑各管段的管道内壁腐蚀及污垢情况，确定有效空间体积折算系数，根据有效空间体积折算系数和管段的管道内体积确定有效空间体积；
有效空间体积折算系数具体可以根据如下规则获取：运行时间在5年以内的有效空间体积折算系数为1.0；运行时间超过20年且天然气介质含水量大于0.5mg/m3的有效空间体积折算系数为0.9998；其他情况，有效空间体积折算系数为0.9998-1.0；
获取每一个管段天然气实时压力和实时温度；压力监测点设置原则具体为：每隔20公里以上，在管道管径改变处或线路切断阀门位置设置压力监测点，在压力监测点采集压力值，并据此获取各管段的天然气实时压力；温度监测点设置点与压力监测点设置点相同；
获取各管段的天然气实时压力，具体为根据相邻两个压力监测点采集的压力值和管道长度，首先采用比例分布法获取各管段两端的天然气实时压力，公式如下：
Figure BDA0003139202960000141
其中，
P0为相邻两个压力监测点中距离气源更近的压力监测点处的实时压力；
P1为待求管段端点的实时压力，是一个管段尾端压力，同时也是下一个管段首端压力；
P2为相邻两个压力监测点中距离气源更远的压力监测点处的实时压力；
S1为待求实时压力管段与离气源更近的压力监测点之间的管道长度；
S为相邻两个压力监测点之间的管道长度。
然后，取每一管段首尾两端的实时压力平均值作为该管段实时天然气压力。
获取每一个管段天然气实时组分，具体包括：
对于单一气源管网，在气源输入处设置在线气相色谱分析仪，气相色谱分析仪实时分析数据作为各管段天然气组分；
对于多气源管网，在每个气源输入处设置在线气相色谱分析仪，实时分析天然气组分，根据每一管段管道中各气源输入量的比例，计算出每一个管段的实时天然气组分。
根据压力、温度和组分获取每一个管段天然气压缩因子；
获取整个天然气管网总储气量，具体为根据各管段管道内天然气的实时压力、温度、组分、有效空间体积和压缩因子，获得每一管段内天然气实时储气量，可以采用如下公式计算：
V＝ZmRT/P
V：天然气储气量，单位为m3；
Z：天然气压缩因子，无量纲；
m：气体质量，单位为kg；
R：摩尔气体常数，单位为J/kg﹒K；
T：天然气温度，单位为K；
P：天然气绝对压力：单位为Pa；
获取预设时间段内整个天然气管网的天然气储气量的变化量；
获取天然气管网购销气差率，具体为整个管网天然气购销气差率采用如下公式计算：
Figure BDA0003139202960000142
其中：
R：天然气购销气差率；
A：预设时间段内管网的天然气输入气量，单位为m3；
B：预设时间段内管网的天然气储气量变化量，单位为m3；
C：预设时间段内管网的天然气输出气量，单位为m3。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均包含在本发明的保护范围之内。

Claims (6)

1.一种天然气管网的购销气差率实时监测方法，其特征在于，包括如下步骤：
根据管道口径变化、管道长度及传感器安装位置选择难易度，对整个天然气管网进行管段划分；
获得各管段及整个天然气管网管道内的有效空间体积；
获取每一个管段天然气实时压力和实时温度；
获取每一个管段天然气实时组分；
根据压力、温度和组分获取每一个管段天然气压缩因子；
根据各管段管道内天然气的实时压力、温度、组分、有效空间体积和压缩因子，获得每一管段内天然气实时储气量；
获取整个天然气管网总储气量；
获取预设时间段内整个天然气管网的天然气储气量的变化量；
获取天然气管网购销气差率；
其中，根据管道口径变化、管道长度及传感器安装位置选择难易度，进行管段划分，包括：
相邻的不同管径的管道，每个管径的管道划分为一个管段；
相邻的管径相同的管道，长度不超过10公里时，划分为一个管段；
相邻的管径相同的管道，长度超过10公里时，划分为2个或多个管段；
其中，获取各管段的管道内的有效空间体积具体包括：依据各管段的管道的运行年限、材质及运行状态，考虑各管段的管道内壁腐蚀及污垢情况，确定有效空间体积折算系数，根据有效空间体积折算系数和管段的管道内体积确定有效空间体积；
其中，有效空间体积折算系数具体为：运行时间在5年以内管段的有效空间体积折算系数为1.0；运行时间超过20年且天然气介质含水量大于0.5mg/m3的有效空间体积折算系数为0.9998；其他情况，有效空间体积折算系数在0.9998-1.0之间采用比例分布法推算。
2.根据权利要求1所述的天然气管网的购销气差率实时监测方法，其特征在于，每隔20公里以上，在管道管径改变处或线路切断阀门位置设置压力监测点，在压力监测点采集压力值，并据此获取各管段的天然气实时压力。
3.根据权利要求2所述的天然气管网的购销气差率实时监测方法，其特征在于，根据相邻两个压力监测点采集的压力值和管道长度，采用比例分布法获取各管段首尾两端的天然气实时压力，再将各管段首尾两端的天然气实时压力取平均值，获取每一个管段的天然气实时压力，其中比例分布法采用如下公式：
Figure 835241DEST_PATH_IMAGE001
其中，
P 0为相邻两个压力监测点中距离气源更近的压力监测点处的实时压力；
P 1为待求管段端点的实时压力，是一个管段尾端压力，同时也是下一个管段首端压力；
P 2为相邻两个压力监测点中距离气源更远的压力监测点处的实时压力；
S 1为待求实时压力管段与离气源更近的压力监测点之间的管道长度；
S为相邻两个压力监测点之间的管道长度。
4.根据权利要求1所述的天然气管网的购销气差率实时监测方法，其特征在于，获取每一个管段天然气实时组分包括：
对于单一气源管网，在气源输入处设置在线气相色谱分析仪，气相色谱分析仪实时分析数据作为各管段天然气实时组分；
对于多气源管网，在每个气源输入处设置在线气相色谱分析仪，实时分析天然气组分，根据每一管段管道中各气源输入量的比例，计算出每一个管段的实时天然气组分。
5.根据权利要求1所述的天然气管网的购销气差率实时监测方法，其特征在于，每一管段内天然气实时储气量采用如下公式计算：
V = ZmRT/P
V：天然气储气量，单位为m3；
Z：天然气压缩因子，无量纲；
m：气体质量，单位为kg；
R：摩尔气体常数，单位为J/kg﹒K；
T：天然气温度，单位为K；
P：天然气绝对压力：单位为Pa。
6.根据权利要求1所述的天然气管网的购销气差率实时监测方法，其特征在于，整个管网天然气购销气差率采用如下公式计算：
Figure 964871DEST_PATH_IMAGE002
其中：
R：天然气购销气差率；
A：预设时间段内管网的天然气输入气量，单位为m3；
B：预设时间段内管网的天然气储气量变化量，单位为m3；
C：预设时间段内管网的天然气输出气量，单位为m3。
CN202110727541.4A 2021-06-29 2021-06-29 一种天然气管网的购销气差率实时监测方法 Active CN113531395B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110727541.4A CN113531395B (zh) 2021-06-29 2021-06-29 一种天然气管网的购销气差率实时监测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110727541.4A CN113531395B (zh) 2021-06-29 2021-06-29 一种天然气管网的购销气差率实时监测方法

Publications (2)

Publication Number Publication Date
CN113531395A CN113531395A (zh) 2021-10-22
CN113531395B true CN113531395B (zh) 2022-11-22

Family

ID=78097134

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110727541.4A Active CN113531395B (zh) 2021-06-29 2021-06-29 一种天然气管网的购销气差率实时监测方法

Country Status (1)

Country Link
CN (1) CN113531395B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6055846A (en) * 1997-03-20 2000-05-02 The Texas A&M University System Method and apparatus for in situ calibration of gas flowmeters
CN108992038A (zh) * 2018-08-10 2018-12-14 清华大学 眼压监测系统及眼压监测方法
CN110765622A (zh) * 2019-10-28 2020-02-07 南方电网科学研究院有限责任公司 一种天然气管道模型的能流获取系统、设备及介质
CN210165192U (zh) * 2019-05-08 2020-03-20 珠海吉泰克燃气设备技术有限公司 燃气智能控制系统
CN111337071A (zh) * 2019-12-14 2020-06-26 中国石油化工股份有限公司天然气分公司计量研究中心 一种天然气计量评价系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112733304A (zh) * 2020-08-25 2021-04-30 濮阳市新宝自动化科技开发有限公司 城市燃气管网多输差协同计算模型及应用方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6055846A (en) * 1997-03-20 2000-05-02 The Texas A&M University System Method and apparatus for in situ calibration of gas flowmeters
CN108992038A (zh) * 2018-08-10 2018-12-14 清华大学 眼压监测系统及眼压监测方法
CN210165192U (zh) * 2019-05-08 2020-03-20 珠海吉泰克燃气设备技术有限公司 燃气智能控制系统
CN110765622A (zh) * 2019-10-28 2020-02-07 南方电网科学研究院有限责任公司 一种天然气管道模型的能流获取系统、设备及介质
CN111337071A (zh) * 2019-12-14 2020-06-26 中国石油化工股份有限公司天然气分公司计量研究中心 一种天然气计量评价系统

Also Published As

Publication number Publication date
CN113531395A (zh) 2021-10-22

Similar Documents

Publication Publication Date Title
CN111486930B (zh) 一种天然气能量计量赋值核查方法及系统
JP6318053B2 (ja) 漏水分布推定装置
NO321700B1 (no) Venturi-stromningsmaler for maling i en fluidstromningskanal
Yıldırım An assessment of hydraulic design of trickle laterals considering effect of minor losses
CN113464850B (zh) 一种天然气管网泄漏监测及应急处置系统
Loureiro et al. A new approach to improve water loss control using smart metering data
Kanshio A review of hydrocarbon allocation methods in the upstream oil and gas industry
Sauer Unit-response method of open-channel flow routing
CN113531395B (zh) 一种天然气管网的购销气差率实时监测方法
US6382032B1 (en) Apparatus and method for measuring flow of gas with entrained liquids
CN116467551A (zh) 一种基于相关系数的供水管网漏损定位方法及系统
Carpenter et al. Applying the IWA approach to water loss performance indicators in Australia
CN108122074B (zh) 一种给水管网模型中漏失水量的分配方法
CN111563332B (zh) 基于夜间调压实测的城镇供水管网漏失定位方法
JP6018970B2 (ja) 配水制御装置および方法
CN109374056B (zh) 3d数字动态模型判评压裂或管道流体性能的方法及系统
Khalifa et al. Analysis and Assessment of Water Losses in Domestic Water Distribution Networks
CN2462377Y (zh) 管道阻力系数测量仪
Schley et al. Gas quality tracking in distribution grids
Langelandsvik et al. Accurate calculations of pipeline transport capacity
Shiu et al. Enhancing the EPANET hydraulic model through genetic algorithm optimization of pipe roughness coefficients
Renaud et al. Introducing pressure and number of connections into water loss indicators for French drinking water supply networks
Halkijević et al. Average pressure in a water supply system
RU2425333C1 (ru) Способ измерения расхода и количества газообразных сред
Trancoso et al. Analysis of interference caused by air blocking valves in residential water meter measurements

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant