CN114626587A - 一种考虑天然气管道掺氢的气电联合系统运行优化方法 - Google Patents

一种考虑天然气管道掺氢的气电联合系统运行优化方法 Download PDF

Info

Publication number
CN114626587A
CN114626587A CN202210193832.4A CN202210193832A CN114626587A CN 114626587 A CN114626587 A CN 114626587A CN 202210193832 A CN202210193832 A CN 202210193832A CN 114626587 A CN114626587 A CN 114626587A
Authority
CN
China
Prior art keywords
gas
natural gas
node
hydrogen
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.)
Pending
Application number
CN202210193832.4A
Other languages
English (en)
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.)
Huazhong University of Science and Technology
Economic and Technological Research Institute of State Grid Hubei Electric Power Co Ltd
Original Assignee
Huazhong University of Science and Technology
Economic and Technological Research Institute of State Grid Hubei Electric Power 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 Huazhong University of Science and Technology, Economic and Technological Research Institute of State Grid Hubei Electric Power Co Ltd filed Critical Huazhong University of Science and Technology
Priority to CN202210193832.4A priority Critical patent/CN114626587A/zh
Publication of CN114626587A publication Critical patent/CN114626587A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06313—Resource planning in a project environment
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/04—Constraint-based CAD
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/14—Pipes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Strategic Management (AREA)
  • Marketing (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Health & Medical Sciences (AREA)
  • Development Economics (AREA)
  • Fluid Mechanics (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Computing Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Algebra (AREA)
  • Educational Administration (AREA)
  • Primary Health Care (AREA)
  • General Health & Medical Sciences (AREA)
  • Water Supply & Treatment (AREA)
  • Public Health (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

一种考虑天然气管道掺氢的气电联合系统运行优化方法，该方法包括以下步骤：S1、收集气电联合系统的基本参数以及气电联合系统中电力负荷、天然气负荷和可再生能源出力的预测曲线；S2、建立天然气网络节点处天然气和氢气的质量流量的计算模型；S3、建立考虑管道掺氢的天然气系统动态模型；S4、建立考虑天然气管道掺氢的气电联合优化运行模型，并对模型进行求解得到日前开机方式。本设计不仅降低了系统运行成本，而且减少了碳排放量。

Description

一种考虑天然气管道掺氢的气电联合系统运行优化方法
技术领域
本发明涉及电气工程领域，尤其涉及一种考虑天然气管道掺氢的气电联合系统运行优化方法，主要适用于降低成本与减少碳排放量。
背景技术
为了应对化石能源过度消耗和全球气候变暖等环境问题，可再生能源在世界范围内得到了广泛的应用，未来的能源供给体系将以电力为主，火电等不可再生能源机组逐步退出、可再生能源占比不断增大是未来能源结构的必然趋势。而风电、光伏等可再生能源出力的波动性为电力系统运行的灵活性提出了更高的要求，以有效保障系统运行的安全性。
燃气机组不仅具有良好的动态特性，对于可再生能源出力波动，能够快速响应，而且较之传统的燃煤机组，环境污染更小，因此在电力系统中的应用越来越广泛。与此同时，随着电转气(power to gas，P2G)技术的产生和发展，电力系统中富余的可再生能源电力能够通过转换成天然气(主要指CH4)并泵入天然气管网的方式实现进一步地消纳，同时进一步加强了电力系统与天然气系统之间的耦合，利用燃气机组和P2G机组实现二者之间的双向耦合，可以在促进可再生能源消纳的同时，提高系统运行的灵活性和安全性。
随着氢能的迅速发展，对氢能储运技术也提出了更高的要求，而利用管道的储气能力，可同时实现氢气的输送与储存，但是建设专门的输氢线路存在前期投资大的问题，对此，如果能够发展混氢天然气技术，利用现有的天然气管道，并通过有限改造实现混合气体的规模化输送，便可通过较少的投资，实现氢气的大规模储运，有效推动能源结构转型和“双碳”目标实现。
发明内容
本发明的目的是克服现有技术中存在的成本高、碳排放量高的缺陷与问题，提供一种成本低、碳排放量低的考虑天然气管道掺氢的气电联合系统运行优化方法。
为实现以上目的，本发明的技术解决方案是：一种考虑天然气管道掺氢的气电联合系统运行优化方法，该方法包括以下步骤：
S1、收集气电联合系统的基本参数以及气电联合系统中电力负荷、天然气负荷和可再生能源出力的预测曲线；
S2、建立天然气网络节点处天然气和氢气的质量流量的计算模型；
S3、建立考虑管道掺氢的天然气系统动态模型；
S4、建立考虑天然气管道掺氢的气电联合优化运行模型，并对模型进行求解得到日前开机方式。
步骤S2中，定义管道内天然气掺氢比例为：
Figure BDA0003526097120000021
式中，RV为天然气掺氢体积比，
Figure BDA0003526097120000022
为t时刻管道内的氢气体积，
Figure BDA0003526097120000023
为t时刻管道内的天然气体积，Δt为时间步长，
Figure BDA0003526097120000024
为t时刻所考虑节点处的氢气的密度，
Figure BDA0003526097120000025
为t时刻所考虑节点处的天然气的密度，
Figure BDA0003526097120000026
为t时刻管道内的氢气的质量，
Figure BDA0003526097120000027
为t时刻管道内的天然气的质量，
Figure BDA0003526097120000028
为t时刻所考虑节点处的氢气的质量流量，
Figure BDA0003526097120000029
为t时刻所考虑节点处的天然气的质量流量；
氢气的质量流量和天然气的质量流量之间的关系为：
Figure BDA00035260971200000210
令
Figure BDA00035260971200000211
则：
Figure BDA00035260971200000212
记混合气体t时刻的质量流量为
Figure BDA00035260971200000213
可得：
Figure BDA0003526097120000031
Figure BDA0003526097120000032
Figure BDA0003526097120000033
气体状态方程为：
Figure BDA0003526097120000034
Figure BDA0003526097120000035
式中，
Figure BDA0003526097120000036
为t时刻所考虑节点处的氢气的压强，
Figure BDA0003526097120000037
为t时刻所考虑节点处的天然气的压强，c为音速；
则标准状态下有：
Figure BDA0003526097120000038
Figure BDA0003526097120000039
联列可得：
Figure BDA00035260971200000310
进一步有：
Figure BDA00035260971200000311
在PV确定的情况下，Kt为一常数，记作K。
步骤S3中，通过以下三个方程来描述混合气体运动过程：
混合气体动量方程：
Figure BDA0003526097120000041
式中，Aij为管道ij横截面积，MFij，t+1为管道ij前端在t+1时刻的混合气体质量流量，MEij，t+1为管道ij末端在t+1时刻的混合气体质量流量，MFij，t为管道ij前端在 t时刻的混合气体质量流量，MEij，t为管道ij末端在t时刻的混合气体质量流量，Δt为时间步长，Lij为管道的长度，pj，t+1为管道ij的末端节点j在t+1时刻的混合气体压强，pi，t+1为管道ij的前端节点i在t+1时刻的混合气体压强，pj，t为管道ij的末端节点j在t时刻的混合气体压强，pi，t为管道ij的前端节点i在t时刻的混合气体压强，λ为摩擦系数，
Figure BDA0003526097120000042
为管道ij内混合气体平均流速，dij为管道ij的直径；
管道物质平衡方程：
Figure BDA0003526097120000043
气体状态方程：
p＝c2ρ
式中，c为音速，p为混合气体压强，ρ为混合气体密度。
步骤S3中，天然气系统中各节点的边界条件约束为：
(1)负荷节点混合气体的质量流量为：
Figure BDA0003526097120000044
Figure BDA0003526097120000045
式中，
Figure BDA0003526097120000051
表示不考虑掺氢时负荷节点i在t时刻的气体负荷质量流量，
Figure BDA0003526097120000052
表示掺氢后负荷节点i在t时刻的混合气体的质量流量，
Figure BDA0003526097120000053
表示氢气的热值，
Figure BDA0003526097120000054
表示天然气的热值，KI为负荷节点集合，MEki，t为管道ki末端在t时刻的混合气体质量流量；
(2)气源节点的混合气体压强和混合气体密度约束为：
Figure BDA0003526097120000055
式中，psj为气源节点的混合气体压强，ρsi为气源节点的混合气体密度；
(3)天然气网络中的每一个中间连接节点的物质平衡约束为：
Figure BDA0003526097120000056
式中，(.)k表示末端节点为k的天然气管道集合，k(.)表示前段节点为k的天然气管道集合，Km为连接节点集合；
(4)管道质量流量与节点气压的上下限约束：
Figure BDA0003526097120000057
Figure BDA0003526097120000058
Figure BDA0003526097120000059
步骤S4中，所述气电联合优化运行模型的目标函数为：
Figure BDA00035260971200000510
式中，Ng为传统火电机组台数，ccoal为标准煤当月实时价格，Ngas为天然气系统气源节点数量，
Figure BDA00035260971200000511
为火电机组i在t时刻的煤耗量，
Figure BDA00035260971200000512
为火电机组i单位启停成本；Sudi，t为火电机组i在t时刻的启停状态变化，为0-1变量，当其t时刻的启停状态不同于上一时刻时，Sudi，t＝1，否则，Sudi，t＝0；cM为天然气的价格，cH为氢气的价格，
Figure BDA0003526097120000061
为天然气网络气源节点注入的混合气体中天然气的质量流量，
Figure BDA0003526097120000062
为天然气网络气源节点注入的混合气体中氢气的质量流量，Δt为时间步长。
步骤S4中，所述气电联合优化运行模型的约束条件为：
(1)功率平衡约束
Figure BDA0003526097120000063
式中，
Figure BDA0003526097120000064
为风电场i在t时刻的风电功率，
Figure BDA0003526097120000065
为负荷j在t时刻的负荷功率，
Figure BDA0003526097120000066
为火电机组i在t时刻的出力，
Figure BDA0003526097120000067
为燃气机组在t时刻的出力，
Figure BDA0003526097120000068
为电转气机组i 在t时刻的消耗功率；
(2)支路传输容量约束
-f1im≤SP≤f1im
式中，flim为线路最大传输功率列向量，S为由节点注入功率确定线路传输功率的灵敏度矩阵，P为由各节点注入功率组成的注入功率列向量；
(3)机组出力上下限约束
Figure BDA0003526097120000069
Figure BDA00035260971200000610
Figure BDA00035260971200000611
式中，ui，t为火电机组i在t时刻的开停机状态，为0-1变量，当其t时刻处于升机状态时，ui，t＝1，否则，ui，t＝0；Pi g，max为火电机组最小出力限值，Pi g，min为火电机组最大出力限值，Pi gas，max为燃气机组最大出力限值，Pi p2g，max为火电机组最大消耗功率；
(4)火电机组启停状态约束
Sudi，1＝0
Figure BDA00035260971200000612
Figure BDA0003526097120000071
(5)火电机组最小启停时间约束
ui，t-1≤ui，t≤Ti U
ui，j≤ui，t-1t≤Ti D
Figure BDA0003526097120000072
Figure BDA0003526097120000073
式中，Ti U为机组i需要维持的最小开机时间，Ti D为机组i需要维持的最小停机时间；
(6)机组爬坡约束
Figure BDA0003526097120000074
Figure BDA0003526097120000075
式中，
Figure BDA0003526097120000076
为火电机组i的最大上爬坡功率，
Figure BDA0003526097120000077
为火电机组i的最大下爬坡功率，M为常数；
(7)燃气机组运行约束
Figure BDA0003526097120000078
式中，
Figure BDA0003526097120000079
为燃气机组i对于掺氢天然气的气电转化效率，
Figure BDA00035260971200000710
为燃气机组i的耗气质量流量；
(8)电转气机组运行约束
Figure BDA00035260971200000711
Figure BDA00035260971200000712
Figure BDA00035260971200000713
式中，
Figure BDA00035260971200000714
为电转气节点i处电转甲烷机组在t时刻的耗电功率，
Figure BDA00035260971200000715
为电转气节点i处电制氢机组在t时刻的耗电功率，
Figure BDA0003526097120000081
为电转气节点i处电转气的产气质量流量，
Figure BDA0003526097120000082
为电转气节点i处电转甲烷机组的制气效率，
Figure BDA0003526097120000083
为电转气节点i处电制氢机组的制气效率。
与现有技术相比，本发明的有益效果为：
本发明一种考虑天然气管道掺氢的气电联合系统运行优化方法中，建立考虑管道掺氢的天然气系统动态模型，模型中根据天然气管道中的掺氢比例建立了天然气网络节点处天然气和氢气的质量流量的计算模型，并根据管道动量方程和物质平衡方程，建立了考虑动态能流的天然气系统模型，更为精确地表达天然气系统的动态运行过程；在此基础上，提出了考虑天然气管道掺氢的气电联合系统运行优化方法，计及天然气系统的动态特性，并充分利用燃气机组、电转甲烷机组、电制氢机组的特性，实现电力系统和天然气系统之间的双向能量耦合，在提高系统灵活性的同时，促进可再生能源消纳，减少碳排放，为混氢天然气管道技术的应用提供参考。
附图说明
图1是本发明一种考虑天然气管道掺氢的气电联合系统运行优化方法的流程图。
图2是本发明的实施例中的39节点电力系统拓扑示意图。
图3是本发明的实施例中的27节点天然气系统拓扑示意图。
图4是本发明的实施例中的电力负荷、风电出力以及天然气负荷曲线图。
图5是本发明的实施例中的火电机组开机方式示意图。
图6是本发明的实施例中的气电联合系统中各机组出力示意图。
具体实施方式
以下结合附图说明和具体实施方式对本发明作进一步详细的说明。
参见图1，一种考虑天然气管道掺氢的气电联合系统运行优化方法，该方法包括以下步骤：
S1、收集气电联合系统的基本参数以及气电联合系统中电力负荷、天然气负荷和可再生能源出力的预测曲线；
确定气电联合系统(即研究的电力系统和天然气系统)的拓扑结构以及二者之间的拓扑联系，确定气电联合系统中各元件、装置、设备的相关参数，相关的元件、装置、设备主要包括：电力系统线路、天然气系统管道、火电机组、燃气机组、风电机组、电转气机组；
收集所研究的气电联合系统中的电力负荷、天然气负荷以及风电出力的预测曲线，作为气电联合系统优化运行的已知输入数据；
S2、建立天然气网络节点处天然气和氢气的质量流量的计算模型；
对于掺氢天然气管道，为保证管道安全可靠运行，采用等比例掺氢的方式，并假定管道内天然气(主要为甲烷CH4)和氢气混合均匀，定义管道内天然气掺氢比例为：
Figure BDA0003526097120000091
对上式进行如下变换：
Figure BDA0003526097120000092
式中，RV为天然气掺氢体积比，即管道内掺入氢气的体积与天然气体积之比，量纲为1；
Figure BDA0003526097120000093
为t时刻管道内的氢气体积，
Figure BDA0003526097120000094
为t时刻管道内的天然气体积，Δt 为时间步长，
Figure BDA0003526097120000095
为t时刻所考虑节点处的氢气的密度，
Figure BDA0003526097120000096
为t时刻所考虑节点处的天然气的密度，
Figure BDA0003526097120000097
为t时刻管道内的氢气的质量，
Figure BDA0003526097120000098
为t时刻管道内的天然气的质量，
Figure BDA0003526097120000099
为t时刻所考虑节点处的氢气的质量流量，
Figure BDA00035260971200000910
为t时刻所考虑节点处的天然气的质量流量；
氢气的质量流量和天然气的质量流量之间的关系为：
Figure BDA00035260971200000911
令
Figure BDA00035260971200000912
则：
Figure BDA00035260971200000913
记混合气体t时刻的质量流量为
Figure BDA0003526097120000101
可得：
Figure BDA0003526097120000102
Figure BDA0003526097120000103
Figure BDA0003526097120000104
气体状态方程为：
Figure BDA0003526097120000105
Figure BDA0003526097120000106
式中，
Figure BDA0003526097120000107
为t时刻所考虑节点处的氢气的压强，
Figure BDA0003526097120000108
为t时刻所考虑节点处的天然气的压强，c为音速；
则标准状态下有：
Figure BDA0003526097120000109
Figure BDA00035260971200001010
联列可得：
Figure BDA00035260971200001011
进一步有：
Figure BDA00035260971200001012
在RV确定的情况下，Kt为一常数，记作K；
S3、建立考虑管道掺氢的天然气系统动态模型；
对于管道内的天然气而言，其传输过程由管道两端的压强所驱动，并与自身的温度、密度等因素相关。管道内的天然气流动过程可以视为状态变量为速度、密度、压强的一维流体运动过程。对于掺氢后天然气管道，假设天然气和氢气均匀混合且混氢体积比保持恒定。
考虑天然气网络节点负荷变化引起的管道内混合气体的慢动态过程，可以通过以下三个方程来描述混合气体运动过程：管道动量方程、管道物质平衡方程以及气体状态方程。
管道动量方程，即Navier-Stokes方程，被用来描述管道内混合气体的动量传输：
Figure BDA0003526097120000111
式中，t和x分别表示时间与空间距离；p为气体压强，单位为Pa；ω为气体流速，单位为m/s；ρ和ρα分别表示与水平面平行以及与水平面成α夹角的气体密度，单位为kg/m3；d为管道直径，单位为m；g为重力加速度，单位为m/s2；λ表示摩擦系数；
Figure BDA0003526097120000112
描述了混合气体在管道内的加速度效应；
Figure BDA0003526097120000113
描述了混合气体的对流效应；
Figure BDA0003526097120000114
描述了混合气体的静力学效应；g(ρ-ρa)sinα描述了水平高度对动量方程的影响；
Figure BDA0003526097120000115
表示二阶偏应力张量分量；
天然气管道内的物质平衡方程描述了混合气体在管道内的流动情况：
Figure BDA0003526097120000116
通过音速c，在管道气体压强与密度之间构建气体状态方程：
p＝c2ρ
基于流体力学的管道内气体物质动量传输方程十分复杂，气体动量传输方程、管道物质平衡方程为复杂偏微分方程，本方法通过一定的假设条件，将这两个方程进行简化。
假设管道气体传输为恒温过程，即忽视温度变化的影响，则音速保持不变；由于对流项，在流体速度接近音速时才会存在，对于管道中混合气体流动过程，该项可以忽略；与此同时，认为管道高度不变，则与高度相关的项也为0；在上述假设条件下，气体动量传输方程可以化简为：
Figure BDA0003526097120000117
由于动量方程中偏应力张量是包含气体流速ω二次项和密度ρ乘积的非线性项，为了使模型线性化，采用混合体气体平均流速
Figure BDA0003526097120000121
来近似表示上式中的二次项：
Figure BDA0003526097120000122
同时采用混合气体质量流量M，单位为kg/s，来表示管道中的混合气体流动情况，其与混合气体密度ρ、流速ω及管道横截面积A(m2)存在如下关系：
M＝ρωA
用质量流量M代替动量方程与管道物质平衡方程中的流速ω，同时利用状态方程将式中密度ρ用压强p表示，可进一步化简：
Figure BDA0003526097120000123
Figure BDA0003526097120000124
经过一定简化后的天然气动量方程和管道物质平衡方程仍然是偏微分方程，无法直接在优化领域直接应用，因此，采用Lax-Wendroff差分法将上述偏微分方程进行差分线性化，Lax-Wendroif差分法一般形式如下所示：
Figure BDA0003526097120000125
Figure BDA0003526097120000126
Figure BDA0003526097120000127
式中，Δt和Δx分别为时间与空间距离的步长；
对
Figure BDA0003526097120000128
进行差分线性化，其中对于空间距离x采用变步长的形式，即对于每一个管道ij的动力方程及物质平衡方程，其Δx等于管道长度Lij，同时管道节点下标i+1也对应管道末端节点j；对于任意管道ij都有线性化后的混合气体动量方程与管道物质平衡方程；
混合气体动量方程：
Figure BDA0003526097120000131
式中，Aij为管道ij横截面积，MFij，t+1为管道ij前端在t+1时刻的混合气体质量流量，MEij，t+1为管道ij末端在t+1时刻的混合气体质量流量，MFij，t为管道ij前端在 t时刻的混合气体质量流量，MEij，t为管道ij末端在t时刻的混合气体质量流量，Δt为时间步长，Lij为管道的长度，pj，t+1为管道ij的末端节点j在t+1时刻的混合气体压强，pi，t+1为管道ij的前端节点i在t+1时刻的混合气体压强，pj，t为管道ij的末端节点j在t时刻的混合气体压强，pi，t为管道ij的前端节点i在t时刻的混合气体压强，λ为摩擦系数，
Figure BDA0003526097120000132
为管道ij内混合气体平均流速，dij为管道ij的直径；
管道物质平衡方程：
Figure BDA0003526097120000133
气体状态方程：
p＝c2ρ
式中，c为音速，p为混合气体压强，ρ为混合气体密度；
天然气系统中各节点的边界条件约束为：
(1)假设氢气和天然气混合均匀成具有统一混合热值的气体，负荷节点混合气体的质量流量为：
Figure BDA0003526097120000134
Figure BDA0003526097120000135
式中，
Figure BDA0003526097120000141
表示不考虑掺氢时负荷节点i在t时刻的气体负荷质量流量，
Figure BDA0003526097120000142
表示掺氢后负荷节点i在t时刻的混合气体的质量流量，
Figure BDA0003526097120000143
表示氢气的热值，
Figure BDA0003526097120000144
表示天然气的热值，KI为负荷节点集合，MEki，t为管道ki末端在t时刻的混合气体质量流量；
(2)气源节点的混合气体压强和混合气体密度约束为：
Figure BDA0003526097120000145
式中，psi为气源节点的混合气体压强，ρsi为气源节点的混合气体密度；
(3)天然气网络中的每一个中间连接节点的物质平衡约束为：
Figure BDA0003526097120000146
式中，(.)k表示末端节点为k的天然气管道集合，k(.)表示前段节点为k的天然气管道集合，Km为连接节点集合；
(4)管道质量流量与节点气压的上下限约束：
Figure BDA0003526097120000147
Figure BDA0003526097120000148
Figure BDA0003526097120000149
S4、建立考虑天然气管道掺氢的气电联合优化运行模型，并对模型进行求解得到日前开机方式；
基于步骤S3中提出的考虑管道掺氢的天然气系统动态模型，考虑燃气机组(GFU)和电转气(P2G)机组在电力系统和天然气系统联合运行过程中的能量双向耦合作用，建立了考虑天然气管道掺氢的气电联合优化运行模型；
所提气电联合优化运行模型的目标为最小化联合系统的运行费用，即电力系统和天然气系统的运行费用之和，采用分段线性化方法对火电机组运行费用进行线性化处理，所述气电联合优化运行模型的目标函数为：
Figure BDA0003526097120000151
式中，Ng为传统火电机组台数，ccoal为标准煤当月实时价格，Ngas为天然气系统气源节点数量，
Figure BDA0003526097120000152
为火电机组i在t时刻的煤耗量(t)，
Figure BDA0003526097120000153
为火电机组i单位启停成本(元/次)；Sudi，t为火电机组i在t时刻的启停状态变化，为0-1变量，当其t时刻的启停状态不同于上一时刻时，Sudi，t＝1，否则，Sudi，t＝0；cM为天然气的价格(元/kg)，cH为氢气的价格(元/kg)，
Figure BDA0003526097120000154
为天然气网络气源节点注入的混合气体中天然气的质量流量(kg/m3)，
Figure BDA0003526097120000155
为天然气网络气源节点注入的混合气体中氢气的质量流量(kg/m3)，Δt为时间步长；
对于考虑天然气管道掺氢的气电联合系统，其除了要满足步骤S3中对于天然气系统的约束，还需要满足电力系统的运行和安全约束，以及天然气系统与电力系统之间的能量耦合约束；
所述气电联合优化运行模型的约束条件为：
(1)功率平衡约束
Figure BDA0003526097120000156
式中，
Figure BDA0003526097120000157
为风电场i在t时刻的风电功率(MW)，
Figure BDA0003526097120000158
为负荷j在t时刻的负荷功率(MW)，
Figure BDA0003526097120000159
为火电机组i在t时刻的出力(MW)，
Figure BDA00035260971200001510
为燃气机组在t时刻的出力(MW)，
Figure BDA00035260971200001511
为电转气机组i在t时刻的消耗功率(MW)；
(2)支路传输容量约束
-flim≤SP≤flim
式中，flim为线路最大传输功率列向量，S为由节点注入功率确定线路传输功率的灵敏度矩阵，P为由各节点注入功率组成的注入功率列向量；
(3)机组出力上下限约束
Figure BDA0003526097120000161
Figure BDA0003526097120000162
Figure BDA0003526097120000163
式中，ui，t为火电机组i在t时刻的升停机状态，为0-1变量，当其t时刻处于升机状态时，ui，t＝1，否则，ui，t＝0；Pi g，max为火电机组最小出力限值，Pi g，min为火电机组最大出力限值，Pi gas，max为燃气机组最大出力限值，Pi p2g，max为火电机组最大消耗功率；
(4)火电机组启停状态约束
Sudi，1＝0
Figure BDA0003526097120000164
Figure BDA0003526097120000165
不考虑初始时刻机组的启停状态变化，故令Sudi，1＝0；
(5)火电机组最小启停时间约束
ui，t-1≤ui，tt≤Ti U
ui，t≤ui，t-1t≤Ti D
Figure BDA0003526097120000166
Figure BDA0003526097120000167
式中，Ti U为机组i需要维持的最小开机时间，Ti D为机组i需要维持的最小停机时间；
(6)机组爬坡约束
Figure BDA0003526097120000168
Figure BDA0003526097120000169
式中，
Figure BDA0003526097120000171
为火电机组i的最大上爬坡功率(MW/h)，
Figure BDA0003526097120000172
为火电机组i的最大下爬坡功率(MW/h)，M为常数；
(7)燃气机组运行约束
Figure BDA0003526097120000173
式中，
Figure BDA0003526097120000174
为燃气机组i对于掺氢天然气的气电转化效率(MW/(kg/s))，
Figure BDA0003526097120000175
为燃气机组i的耗气质量流量(kg/s)；
(8)电转气机组运行约束
Figure BDA0003526097120000176
Figure BDA0003526097120000177
Figure BDA0003526097120000178
式中，
Figure BDA0003526097120000179
为电转气节点i处电转甲烷(P2M)机组在t时刻的耗电功率(MW)，
Figure BDA00035260971200001710
为电转气节点i处电制氢(P2H)机组在t时刻的耗电功率，
Figure BDA00035260971200001711
为电转气节点i处电转气的产气质量流量(kg/s)，
Figure BDA00035260971200001712
为电转气节点i处电转甲烷(P2M)机组的制气效率，
Figure BDA00035260971200001713
为电转气节点i处电制氢(P2H)机组的制气效率((kg/s)/MW)。
本实施例针对39节点电力系统与27节点天然气系统组成的气电联合系统展开分析，如图2、图3所示，考虑天然气管道等比例掺氢，掺氢体积比设为5％；该系统共有6台燃煤机组，2台燃气机组(GFU)，2台风电机组，2台P2G机组；电力系统节点 32、33通过燃气机组与天然气系统节点13、22相连，电力系统节点17、22通过P2G 机组与天然气系统节点27、8相连，机组参数如表1所示。气电联合系统中电力负荷、风电预出力以及天然气负荷曲线如图4所示。煤价格为500元/t，天然气价格为2元 /kg，氢气价格为20元/kg。
根据收集到的上述参数，按本方法建立该气电联合系统的模型，进行求解，系统的运行费用如表2所示，求解得到的火电机组开机方式如图5所示，相应的系统运行策略如图6所示，燃气机组和电转气机组的参与使系统运行更加灵活。
天然气折标煤系数取1.674吨标煤/吨天然气，标煤碳排放系数取2.66吨CO2/ 吨标煤，分别计算考虑天然气管道掺氢和天然气管道不掺氢情况下气电联合系统的碳排放量，通过对比可以发现，考虑天然气管道掺氢情况下CO2排放量较之天然气管道不掺氢情况下CO2排放量，减少了18吨，如表3所示；清洁能源氢气的参与使系统碳排放量减少，随着天然气管道掺氢技术的发展与成熟，管道掺氢运行将可以实现，并且掺氢比可以进一步提高，减排效益也会进一步提升，有利于目标的实现。由此可见，较之传统的模型，本设计方法可以在保障系统运行的前提下有效减少碳排放。
表1掺氢气电联合系统机组参数
Figure BDA0003526097120000181
表2气电联合系统运行费用
Figure BDA0003526097120000191
表3二氧化碳排放量对比
Figure BDA0003526097120000192

Claims (6)

1.一种考虑天然气管道掺氢的气电联合系统运行优化方法，其特征在于，该方法包括以下步骤：
S1、收集气电联合系统的基本参数以及气电联合系统中电力负荷、天然气负荷和可再生能源出力的预测曲线；
S2、建立天然气网络节点处天然气和氢气的质量流量的计算模型；
S3、建立考虑管道掺氢的天然气系统动态模型；
S4、建立考虑天然气管道掺氢的气电联合优化运行模型，并对模型进行求解得到日前开机方式。
2.根据权利要求1所述的一种考虑天然气管道掺氢的气电联合系统运行优化方法，其特征在于：
步骤S2中，定义管道内天然气掺氢比例为：
Figure FDA0003526097110000011
式中，RV为天然气掺氢体积比，
Figure FDA0003526097110000012
为t时刻管道内的氢气体积，
Figure FDA0003526097110000013
为t时刻管道内的天然气体积，Δt为时间步长，
Figure FDA0003526097110000014
为t时刻所考虑节点处的氢气的密度，
Figure FDA0003526097110000015
为t时刻所考虑节点处的天然气的密度，
Figure FDA0003526097110000016
为t时刻管道内的氢气的质量，
Figure FDA0003526097110000017
为t时刻管道内的天然气的质量，
Figure FDA0003526097110000018
为t时刻所考虑节点处的氢气的质量流量，
Figure FDA0003526097110000019
为t时刻所考虑节点处的天然气的质量流量；
氢气的质量流量和天然气的质量流量之间的关系为：
Figure FDA0003526097110000021
令
Figure FDA0003526097110000022
则：
Figure FDA0003526097110000023
记混合气体t时刻的质量流量为
Figure FDA0003526097110000024
可得：
Figure FDA0003526097110000025
Figure FDA0003526097110000026
Figure FDA0003526097110000027
气体状态方程为：
Figure FDA0003526097110000028
Figure FDA0003526097110000029
式中，
Figure FDA00035260971100000210
为t时刻所考虑节点处的氢气的压强，
Figure FDA00035260971100000211
为t时刻所考虑节点处的天然气的压强，c为音速；
则标准状态下有：
Figure FDA00035260971100000212
Figure FDA00035260971100000213
联列可得：
Figure FDA00035260971100000214
进一步有：
Figure FDA00035260971100000215
在RV确定的情况下，Kt为一常数，记作K。
3.根据权利要求2所述的一种考虑天然气管道掺氢的气电联合系统运行优化方法，其特征在于：
步骤S3中，通过以下三个方程来描述混合气体运动过程：
混合气体动量方程：
Figure FDA0003526097110000031
式中，Aij为管道ij横截面积，MFij，t+1为管道ij前端在t+1时刻的混合气体质量流量，MEij，t+1为管道ij末端在t+1时刻的混合气体质量流量，MFij，t为管道ij前端在t时刻的混合气体质量流量，MEij，t为管道ij末端在t时刻的混合气体质量流量，Δt为时间步长，Lij为管道的长度，pj，t+1为管道ij的末端节点j在t+1时刻的混合气体压强，pi，t+1为管道ij的前端节点i在t+1时刻的混合气体压强，pj，t为管道ij的末端节点j在t时刻的混合气体压强，pi，t为管道ij的前端节点i在t时刻的混合气体压强，λ为摩擦系数，
Figure FDA0003526097110000032
为管道ij内混合气体平均流速，dij为管道ij的直径；
管道物质平衡方程：
Figure FDA0003526097110000033
气体状态方程：
p＝c2ρ
式中，c为音速，p为混合气体压强，ρ为混合气体密度。
4.根据权利要求3所述的一种考虑天然气管道掺氢的气电联合系统运行优化方法，其特征在于：步骤S3中，天然气系统中各节点的边界条件约束为：
(1)负荷节点混合气体的质量流量为：
Figure FDA0003526097110000041
Figure FDA0003526097110000042
式中，
Figure FDA0003526097110000043
表示不考虑掺氢时负荷节点i在t时刻的气体负荷质量流量，
Figure FDA0003526097110000044
表示掺氢后负荷节点i在t时刻的混合气体的质量流量，
Figure FDA0003526097110000045
表示氢气的热值，
Figure FDA0003526097110000046
表示天然气的热值，Kl为负荷节点集合，MEki，t为管道ki末端在t时刻的混合气体质量流量；
(2)气源节点的混合气体压强和混合气体密度约束为：
Figure FDA0003526097110000047
式中，psi为气源节点的混合气体压强，ρsi为气源节点的混合气体密度；
(3)天然气网络中的每一个中间连接节点的物质平衡约束为：
Figure FDA0003526097110000048
式中，(.)k表示末端节点为k的天然气管道集合，k(.)表示前段节点为k的天然气管道集合，Km为连接节点集合；
(4)管道质量流量与节点气压的上下限约束：
Figure FDA0003526097110000049
Figure FDA00035260971100000410
Figure FDA00035260971100000411
5.根据权利要求4所述的一种考虑天然气管道掺氢的气电联合系统运行优化方法，其特征在于：步骤S4中，所述气电联合优化运行模型的目标函数为：
Figure FDA0003526097110000051
式中，Ng为传统火电机组台数，ccoal为标准煤当月实时价格，Ngas为天然气系统气源节点数量，
Figure FDA0003526097110000052
为火电机组i在t时刻的煤耗量，
Figure FDA0003526097110000053
为火电机组i单位启停成本；Sudi，t为火电机组i在t时刻的启停状态变化，为0-1变量，当其t时刻的启停状态不同于上一时刻时，Sudi，t＝1，否则，Sudi，t＝0；cM为天然气的价格，cH为氢气的价格，
Figure FDA0003526097110000054
为天然气网络气源节点注入的混合气体中天然气的质量流量，
Figure FDA0003526097110000055
为天然气网络气源节点注入的混合气体中氢气的质量流量，Δt为时间步长。
6.根据权利要求5所述的一种考虑天然气管道掺氢的气电联合系统运行优化方法，其特征在于：步骤S4中，所述气电联合优化运行模型的约束条件为：
(1)功率平衡约束
Figure FDA0003526097110000056
式中，
Figure FDA0003526097110000057
为风电场i在t时刻的风电功率，
Figure FDA0003526097110000058
为负荷j在t时刻的负荷功率，
Figure FDA0003526097110000059
为火电机组i在t时刻的出力，
Figure FDA00035260971100000510
为燃气机组在t时刻的出力，
Figure FDA00035260971100000511
为电转气机组i在t时刻的消耗功率；
(2)支路传输容量约束
-flim≤SP≤flim
式中，flim为线路最大传输功率列向量，S为由节点注入功率确定线路传输功率的灵敏度矩阵，P为由各节点注入功率组成的注入功率列向量；
(3)机组出力上下限约束
Figure FDA00035260971100000512
Figure FDA00035260971100000513
Figure FDA0003526097110000061
式中，ui，t为火电机组i在t时刻的开停机状态，为0-1变量，当其t时刻处于开机状态时，ui，t＝1，否则，ui，t＝0；Pi g，max为火电机组最小出力限值，Pi g，min为火电机组最大出力限值，Pi gas，max为燃气机组最大出力限值，Pi p2g，max为火电机组最大消耗功率；
(4)火电机组启停状态约束
Sudi，1＝0
Sudi，t≥ui，t-1-ui，t
Figure FDA0003526097110000062
Sudi，t≥ui，t-ui，t-1
Figure FDA0003526097110000063
(5)火电机组最小启停时间约束
ui，t-1≤ui，t t≤Ti U
ui，t≤ui，t-1 t≤Ti D
Figure FDA0003526097110000064
t≥Ti U+1
Figure FDA0003526097110000065
t≥Ti D+1
式中，Ti U为机组i需要维持的最小开机时间，Ti D为机组i需要维持的最小停机时间；
(6)机组爬坡约束
Figure FDA0003526097110000066
Figure FDA0003526097110000067
式中，
Figure FDA0003526097110000068
为火电机组i的最大上爬坡功率，
Figure FDA0003526097110000069
为火电机组i的最大下爬坡功率，M为常数；
(7)燃气机组运行约束
Figure FDA00035260971100000610
式中，
Figure FDA0003526097110000071
为燃气机组i对于掺氢天然气的气电转化效率，
Figure FDA0003526097110000072
为燃气机组i的耗气质量流量；
(8)电转气机组运行约束
Figure FDA0003526097110000073
Figure FDA0003526097110000074
Figure FDA0003526097110000075
式中，
Figure FDA0003526097110000076
为电转气节点i处电转甲烷机组在t时刻的耗电功率，
Figure FDA0003526097110000077
为电转气节点i处电制氢机组在t时刻的耗电功率，
Figure FDA0003526097110000078
为电转气节点i处电转气的产气质量流量，
Figure FDA0003526097110000079
为电转气节点i处电转甲烷机组的制气效率，
Figure FDA00035260971100000710
为电转气节点i处电制氢机组的制气效率。
CN202210193832.4A 2022-03-01 2022-03-01 一种考虑天然气管道掺氢的气电联合系统运行优化方法 Pending CN114626587A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210193832.4A CN114626587A (zh) 2022-03-01 2022-03-01 一种考虑天然气管道掺氢的气电联合系统运行优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210193832.4A CN114626587A (zh) 2022-03-01 2022-03-01 一种考虑天然气管道掺氢的气电联合系统运行优化方法

Publications (1)

Publication Number Publication Date
CN114626587A true CN114626587A (zh) 2022-06-14

Family

ID=81899967

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210193832.4A Pending CN114626587A (zh) 2022-03-01 2022-03-01 一种考虑天然气管道掺氢的气电联合系统运行优化方法

Country Status (1)

Country Link
CN (1) CN114626587A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024016858A1 (zh) * 2022-07-22 2024-01-25 上海能源建设工程设计研究有限公司 一种基于干旱胁迫场景的农光气共生系统

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024016858A1 (zh) * 2022-07-22 2024-01-25 上海能源建设工程设计研究有限公司 一种基于干旱胁迫场景的农光气共生系统

Similar Documents

Publication Publication Date Title
Zhou et al. Modeling and simulation of the hydrogen blended gas-electricity integrated energy system and influence analysis of hydrogen blending modes
CN109919399B (zh) 一种综合能源系统日前经济调度方法及系统
CN110544025B (zh) 电转气-储气罐组合的气电综合能源系统优化调度方法
CN109086474A (zh) 基于二阶锥优化算法的天然气系统最优潮流求解方法
CN110163415B (zh) 一种变工况特性下的多能流系统多目标模糊协同优化方法
Gao et al. Wind-photovoltaic co-generation prediction and energy scheduling of low-carbon complex regional integrated energy system with hydrogen industry chain based on copula-MILP
CN113792969A (zh) 考虑气网动态特性和电转气综合能源系统优化调度方法
CN110729720A (zh) 一种蓄热罐容量规划方法和电转气参数敏感性分析方法
Zhang et al. Bi-level optimization dispatch of integrated-energy systems with P2G and carbon capture
CN109685291A (zh) 基于多能互补的区域综合能源系统日前优化分析方法
CN114626587A (zh) 一种考虑天然气管道掺氢的气电联合系统运行优化方法
CN111415045A (zh) 一种乡村生物质综合能源系统两阶段运行优化方法
CN113141005B (zh) 一种面向新能源消纳的综合能源系统多时间尺度调度方法
CN112862158B (zh) 一种综合能源系统日前优化调度方法及系统
Su et al. An interval based cost-emissions optimization strategy for gas-electricity integrated energy systems under uncertainty and demand response
CN106100002A (zh) 一种交直流混合微网的优化运行方法
CN110457745B (zh) 一种具有紧密耦合关系的水-能联系网络模型的建立方法
Hong et al. Mechanism of upgrading low-grade solar thermal energy and experimental validation
CN117411049A (zh) 一种用于电力系统进行储能优化运行的配置方法
CN115455709A (zh) 计及碳捕集设备加装的低碳综合能源系统仿真及配置方法
CN114139349A (zh) 计及管道爬坡的气电联合系统日前机组组合优化方法
CN114997662A (zh) 含电动车充电桩的区域分布式多能源系统低碳规划方法
CN113077173A (zh) 一种基于综合能源系统节点能价策略的综合需求响应方法
CN111106627A (zh) 考虑压缩空气储能的风力发电系统鲁棒优化调度方法
Niu et al. Low-Carbon Oriented Optimization of Integrated Energy System Considering Complex Coupling of Carbon and Hydrogen

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