CN116611706A - 基于多能源主体的动态碳排放因子测算方法 - Google Patents
基于多能源主体的动态碳排放因子测算方法 Download PDFInfo
- Publication number
- CN116611706A CN116611706A CN202310577830.XA CN202310577830A CN116611706A CN 116611706 A CN116611706 A CN 116611706A CN 202310577830 A CN202310577830 A CN 202310577830A CN 116611706 A CN116611706 A CN 116611706A
- Authority
- CN
- China
- Prior art keywords
- network
- node
- flow
- carbon emission
- power
- 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
Links
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 title claims abstract description 78
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 78
- 238000000034 method Methods 0.000 title claims abstract description 43
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims abstract description 66
- 239000003345 natural gas Substances 0.000 claims abstract description 33
- 238000004364 calculation method Methods 0.000 claims abstract description 30
- 230000008878 coupling Effects 0.000 claims abstract description 10
- 238000010168 coupling process Methods 0.000 claims abstract description 10
- 238000005859 coupling reaction Methods 0.000 claims abstract description 10
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 33
- 239000011159 matrix material Substances 0.000 claims description 31
- 238000012937 correction Methods 0.000 claims description 6
- 239000007789 gas Substances 0.000 claims description 6
- 238000002347 injection Methods 0.000 claims description 5
- 239000007924 injection Substances 0.000 claims description 5
- 238000005265 energy consumption Methods 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 4
- 230000004907 flux Effects 0.000 claims description 3
- 238000010438 heat treatment Methods 0.000 claims description 3
- 230000005611 electricity Effects 0.000 claims description 2
- 241000283690 Bos taurus Species 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 238000004088 simulation Methods 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- 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
-
- 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/0637—Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
- G06Q10/06375—Prediction of business process outcome or impact based on a proposed change
-
- 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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/80—Management or planning
- Y02P90/84—Greenhouse gas [GHG] management systems
Abstract
本发明公开了一种基于多能源主体的动态碳排放因子测算方法,涉及碳排放技术领域。所述方法通过对综合能源系统的运行构造和系统耦合情况,确定综合能源系统的潮流计算顺序,并将电力系统潮流的计算方法,推广至热力网络和天然气网络,进行综合能源系统潮流计算。运用碳排放流理论对综合能源系统的碳排放责任进行分析,得出综合能源系统的动态电碳因子测算模型,并引入典型算例对综合能源系统碳排放情况进行分析验证,确定了计算方案的可行性和有效性,精准构建了综合能源系统碳排放流模型。
Description
技术领域
本发明涉及碳排放技术领域,尤其涉及一种基于多能源主体的动态碳排放因子测算方法。
背景技术
随着气候变化的加剧和全球碳排放的增加,各国政府和企业纷纷采取行动减少碳排放并转型为可持续能源。特别是对于高载能企业(如钢铁、石化等),能源转型更是至关重要。由于碳排放与经济发展和能源使用方式密切相关,在进行动态碳排放因子测算时往往会参考经济发展预测和能源消费情况,不同的经济、能耗参考可以认为是不同的未来发展情景,故也可以使用多情景模拟法结合其他算法进行预测,但是现有技术中出现的情景模拟方法一般具有测量不准确,误差大等缺点。
发明内容
本发明所要解决的技术问题是提供一种有效且测量准确的基于多能源主体的动态碳排放因子测算方法。
为解决上述技术问题,本发明所采取的技术方案是:一种基于多能源主体的动态碳排放因子测算方法,包括如下步骤:
对综合能源系统的网络结构进行分析;
根据网络结构对综合能源系统进行解耦;
引入平衡节点计算热力网络系统与天然气网络系统潮流;
将网络平衡节点功率作为电力网络负荷节点计算电力网络系统潮流;
对综合能源系统网络进行精度判断,若精度未达到要求则修正雅可比矩阵再次进行迭代,若精度达到要求则进行碳流计算;
导出支路潮流分布矩阵、机组注入分布矩阵、负荷分布矩阵、节点有功通量矩阵,计算节点碳势分布向量,得出支路碳流率分布矩阵;
核算负荷碳流率向量,得出动态碳排放因子;
由动态碳排放因子分布情况调整综合能源系统网络结构;
使用综合能源系统网络结构对高载能负荷能源进行决策。
进一步的技术方案在于:所述能源测量维度考虑电、气、热三系统,将综合能源系统进行耦合分析,解耦计算,测量网络各个支路碳流率。
进一步的技术方案在于:识别能源网络碳排放情况影响因素包括负荷情况、能源消耗、能源强度、能源结构和技术水平。
进一步的技术方案在于:所述不同网络能源负荷结构下,进行网络碳排放情况预测,得到不同的能源系统碳排放因子。
进一步的技术方案在于:根据所述碳排放因子,可进一步计算碳排放总量、单位电量碳排放强度、碳源结构和能源碳排放结构。
采用上述技术方案所产生的有益效果在于:所述方法通过对综合能源系统的运行构造和系统耦合情况,确定综合能源系统的潮流计算顺序,并将电力系统潮流的计算方法,推广至热力网络和天然气网络,进行综合能源系统潮流计算。运用碳排放流理论对综合能源系统的碳排放责任进行分析,得出综合能源系统的动态电碳因子测算模型,并引入典型算例对综合能源系统碳排放情况进行分析验证,确定了计算方案的可行性和有效性,精准构建了综合能源系统碳排放流模型。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1是本发明实施例所述方法的流程图;
图2是本发明实施例所述方法中电力系统的潮流计算流程图。
具体实施方式
下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是本发明还可以采用其他不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广,因此本发明不受下面公开的具体实施例的限制。
本申请所述方法从潮流分析的角度出发建立考虑电-热-气三系统网络计算模型。借助电-热-气三系统网络在计算输出参数时均需对节点求解的共性,且系统节点参数关系均可遵从牛顿第三定律推导得到类同力、速度、功等物理关系,可通过能量守恒定律将三系统连接,将电力系统的潮流求解扩展至三系统网络。
电力系统潮流求解基于电路的节点电压方程进行,将功率与电流作为未知量代入节点电压方程计算各个节点的电压幅值与相角,进而得到节点功率。这一过程常采用牛拉法与PQ分解法完成,采用牛拉法进行求解。
参照电力系统潮流求解过程,在热网网络中基于管道压降构建方程,将流速、流量作为求解量代换压降,选取初值进行泰勒展开。由于天然气网络的稳态分析思路与电力网络稳态潮流分析方法类似,对网络节点进行求解,网络处理方式与热力网络相似,仍基于管道压降构建方程。
通过三系统网络的共同点从电力潮流推广至能源潮流,利用牛拉法结合热网模型、气网模型构建电-热-气三系统耦合综合能源系统的潮流计算模型。
如图1所示,本发明公开了一种基于多能源主体的动态碳排放因子测算方法,包括如下步骤:
对综合能源系统的网络结构进行分析;
根据网络结构对综合能源系统进行解耦;
引入平衡节点计算热力网络系统与天然气网络系统潮流;
将网络平衡节点功率作为电力网络负荷节点计算电力网络系统潮流;
对综合能源系统网络进行精度判断,若精度未达到要求则修正雅可比矩阵再次进行迭代,若精度达到要求则进行碳流计算;
导出支路潮流分布矩阵、机组注入分布矩阵、负荷分布矩阵、节点有功通量矩阵,计算节点碳势分布向量,得出支路碳流率分布矩阵;
核算负荷碳流率向量,得出动态碳排放因子;
由动态碳排放因子分布情况调整综合能源系统网络结构;
使用综合能源系统网络结构对高载能负荷能源进行决策。
电力系统潮流模型:
为计算电力系统中各节点信息,基于线路各节点的运行条件和网络连接结构确定整个系统的运行状态利用牛拉法进行潮流计算。
对电力网络进行分析,构建功率方程,可写成:
式中Pi、Qi、Vi、Yij为节点有功功率、无功功率、节点i的电压和导纳,将电压的表达式代入式(1)中,展开并将实部和虚部分开可得:
通过对电力系统节点进行编号调整,可将系统中PV、PQ和平衡节点进行串编。设编号方式为将电力系统中的1,2,...,m号节点编号为PQ类型节点,m+1,m+2,…,n-1号节点为PV节点,n号节点为平衡节点。对PQ节点可列方程:
对PV节点可列方程:
ΔWn、ΔVn、Jnn为功率修正量、电压修正量、雅可比矩阵元素。修正方程可写成分块矩阵的形式:
通过对已知节点数据代入形成导纳矩阵,选取电压初始值并代入节点电压方程,计算雅可比矩阵元素,并计算PV节点、PQ节点未知参数。对结算结果进行误差分析,当误差未满足时修正节点电压代入雅可比矩阵再次计算误差,直至误差符合要求。上述计算流程如图2所示。
热力系统潮流模型:
热力系统网络同电力网络类似,由供给端(热源)、网络(热网)、负荷端(热负荷)构成。其中热网结构不同于电力网络,由于单一管道无法完成室内热交换,需由拓扑结构相同的进水管道与回水管道构成,但其作用可等同于电力网络。热网中变量包括各节点压强hi、各负荷需水量Li、各管道水流量mij。节点压强等同于电力网络中电压,负荷需水量等同于电力网络中电流,而管道水流量则可看作热力网络的功率。即使在管道外包有隔热层,在热量在管道传输的过程中仍存在热量的流失,且水流在管道中需克服管壁摩擦做功,故热力网络中也存在网络损耗。热力供水网络中热量通过管道内水流进行传递,在负荷侧通过热交换器将热量传递给用户,并回流至热源处再次进行加热。为简化计算,认为供水网络与回水网络中各管道流量近似一致。
对供水网络进行分析,基于管道内流量近似相等的假设,参照电力系统中基尔霍夫电流方程和欧姆定律方程,构建了热力网络压力平衡方程。
根据基尔霍夫电流定律在热力系统的扩展,依据当流量一定时,任一通流截面上的通流面积与流速的乘积为定值,对通流截面上的截面积与流速构建方程,,而通流面积与流速可转化为水流量,其具体公式为:
其中A、m、mq依次为支路关联矩阵、管道水流量和该节点注入水流量。
参照欧姆定理,可通过管道内压强损失关系,构建压力损失方程表征管道内压力差与管道内水流量的关系,其具体方程为:
其中K为供水管道的阻抗系数。
基于热网各节点之间温度关系构建热力方程,通过负荷输入节点的温度Ts、负荷输出节点的温度To和返回温度Tr表征热力网络种各节点的温度关系。输出温度的具体公式为:
热负荷与温度、水流量的关系为:
供水网络节点温度之间关系:
回水网络节点温度之间关系:
供水压强可得:
实际热力系统中,因热力系统种热源常独立连接,少有多热源线路的存在,故常按照辐射网络对热力系统进行计算,即可参照电力网络辐射网潮流计算方法,构建热力系统潮流计算方程。
天然气系统潮流分析:
天然气网络系统潮流求解基于基尔霍夫气体流量定律,通过构建管道压力方程,将流量与流速作为未知量代入管道压力方程,得到各节点压力,进而得到各节点流量。电网、天然气网建立的网络方程均为非线性方程组,对网络方程的求解可采用迭代法,两者具有相似性,为此,电力网络牛顿迭代法可类比应用在气网中[15]。
将天然气网络类比电力网络,构建流量方程如下:
式中:
fij为天然气流量,p为节点i处的压力。
通过将天然气网络的各部分计算转换为消耗的天然气流量,类比于电力网络的节点电压方程,利用天然气网络的流量构建节点流量平衡方程:
写成分块矩阵的形式为:
综合能源系统碳排放因子测算模型
在计算碳排放因子时,应先对系统潮流进行计算。为便于系统参数的传递,优先计算天然气系统潮流。对天然气网络进行分析,构建天然气系统的潮流模型,再结合热力系统的流量方程和能量方程,基于热力辐射网结构,构建热力系统潮流模型。计算电力系统潮流时,将天然气网络、热力网络与电力网络的接入端口视作负荷。进行热力系统潮流分析时,将电力网络接入端口视作机组,天然气网络接入端口视作负荷。进行天然气系统潮流分析时,将天然气网络和热力网络接入端口视为机组,分别研究各系统的潮流,然后研究电力系统、热力系统和天然气系统的耦合碳排放因子模型。
研究电热耦合系统时,以电热系统能量损耗最小为优化目标:
式中:SGH、SCHP、PGi、Dsa依次为供水机组集合、回水机组集合、供水机组i的有功出力和供热功率;a2i~a0i、e0a~e5a为拟合常数。
电力系统的有功功率和无功功率应保持平衡,据此构建电力系统线路和节点的有功功率和无功功率约束,确定综合能源系统的上下限:
构建发电机组的出力约束方程:
为保证电力系统稳定运行不被破坏,需根据电力系统网络节点和线路情况,构建节点电压幅值约束、线路的潮流约束以及发电机的容量约束。
热力系统也需满足热力网络潮流约束:
热力系统稳定运行也需满足管道流量和温度约束:
研究气电系统时,以气电系统能量损耗最小为优化目标:
式中,Sq为天然气网络注入机组,b2i~b0i为拟合常数。研究电气网络约束时,除电力网络自身约束外,还有天然气网络约束条件:
所述方法通过对综合能源系统的运行构造和系统耦合情况,确定综合能源系统的潮流计算顺序,并将电力系统潮流的计算方法,推广至热力网络和天然气网络,进行综合能源系统潮流计算。运用碳排放流理论对综合能源系统的碳排放责任进行分析,得出综合能源系统的动态电碳因子测算模型,并引入典型算例对综合能源系统碳排放情况进行分析验证,确定了计算方案的可行性和有效性,精准构建了综合能源系统碳排放流模型。
Claims (9)
1.一种基于多能源主体的动态碳排放因子测算方法,其特征在于包括如下步骤:
对综合能源系统的网络结构进行分析;
根据网络结构对综合能源系统进行解耦;
引入平衡节点计算热力网络系统与天然气网络系统潮流;
将网络平衡节点功率作为电力网络负荷节点计算电力网络系统潮流;
对综合能源系统网络进行精度判断,若精度未达到要求则修正雅可比矩阵再次进行迭代,若精度达到要求则进行碳流计算;
导出支路潮流分布矩阵、机组注入分布矩阵、负荷分布矩阵、节点有功通量矩阵,计算节点碳势分布向量,得出支路碳流率分布矩阵;
核算负荷碳流率向量,得出动态碳排放因子;
由动态碳排放因子分布情况调整综合能源系统网络结构;
使用综合能源系统网络结构对高载能负荷能源进行决策。
2.如权利要求1所述的基于多能源主体的动态碳排放因子测算方法,其特征在于,利用牛拉法结合热网模型、气网模型构建电-热-气三系统耦合综合能源系统的潮流计算模型:
对电力网络进行分析,构建功率方程,可写成:
式中Pi、Qi、Vi、Yij为节点有功功率、无功功率、节点i的电压和导纳,将电压的表达式代入式(1)中,展开并将实部和虚部分开可得:
通过对电力系统节点进行编号调整,将系统中PV、PQ和平衡节点进行串编;设编号方式为将电力系统中的1,2,...,m号节点编号为PQ类型节点,m+1,m+2,…,n-1号节点为PV节点,n号节点为平衡节点;对PQ节点可列方程:
对PV节点可列方程:
ΔWn、ΔVn、Jnn为功率修正量、电压修正量、雅可比矩阵元素,修正方程可写成分块矩阵的形式:
通过对已知节点数据代入形成导纳矩阵,选取电压初始值并代入节点电压方程,计算雅可比矩阵元素,并计算PV节点、PQ节点未知参数;对结算结果进行误差分析,当误差未满足时修正节点电压代入雅可比矩阵再次计算误差,直至误差符合要求。
3.如权利要求2所述的基于多能源主体的动态碳排放因子测算方法,其特征在于,热力系统潮流模型:
根据基尔霍夫电流定律在热力系统的扩展,依据当流量一定时,任一通流截面上的通流面积与流速的乘积为定值,对通流截面上的截面积与流速构建方程,而通流面积与流速可转化为水流量,其具体公式为:
其中A、m、mq依次为支路关联矩阵、管道水流量和该节点注入水流量;
参照欧姆定理,可通过管道内压强损失关系,构建压力损失方程表征管道内压力差与管道内水流量的关系,其具体方程为:
其中K为供水管道的阻抗系数;
基于热网各节点之间温度关系构建热力方程,通过负荷输入节点的温度Ts、负荷输出节点的温度To和返回温度Tr表征热力网络种各节点的温度关系,输出温度的具体公式为:
热负荷与温度、水流量的关系为:
供水网络节点温度之间关系:
回水网络节点温度之间关系:
供水压强可得:
构建热力系统潮流计算方程:
4.如权利要求3所述的基于多能源主体的动态碳排放因子测算方法,其特征在于,天然气网络系统潮流分析:
将天然气网络类比电力网络,构建流量方程如下:
式中:
fij为天然气流量,p为节点i处的压力;
通过将天然气网络的各部分计算转换为消耗的天然气流量,类比于电力网络的节点电压方程,利用天然气网络的流量构建节点流量平衡方程:
写成分块矩阵的形式为:
5.如权利要求4所述的基于多能源主体的动态碳排放因子测算方法,其特征在于,综合能源系统碳排放因子测算模型:
研究电热耦合系统时,以电热系统能量损耗最小为优化目标:
式中:SGH、SCHP、PGi、Dsa依次为供水机组集合、回水机组集合、供水机组i的有功出力和供热功率;a2i~a0i、e0a~e5a为拟合常数;
电力系统的有功功率和无功功率应保持平衡,据此构建电力系统线路和节点的有功功率和无功功率约束,确定综合能源系统的上下限:
构建发电机组的出力约束方程:
为保证电力系统稳定运行不被破坏,需根据电力系统网络节点和线路情况,构建节点电压幅值约束、线路的潮流约束以及发电机的容量约束;
热力系统也需满足热力网络潮流约束:
热力系统稳定运行也需满足管道流量和温度约束:
研究气电系统时,以气电系统能量损耗最小为优化目标:
式中,Sq为天然气网络注入机组,b2i~b0i为拟合常数;研究电气网络约束时,除电力网络自身约束外,还有天然气网络约束条件:
6.如权利要求1所述的基于多能源主体的动态碳排放因子测算方法,其特征在于:所述能源测量维度考虑电、气、热三系统,将综合能源系统进行耦合分析,解耦计算,测量网络各个支路碳流率。
7.根据权利要求1所述的基于多能源主体的动态碳排放因子测算方法,其特征在于:识别能源网络碳排放情况影响因素包括负荷情况、能源消耗、能源强度、能源结构和技术水平。
8.根据权利要求1所述的基于多能源主体的动态碳排放因子测算方法,其特征在于,所述不同网络能源负荷结构下,进行网络碳排放情况预测,得到不同的能源系统碳排放因子。
9.根据权利要求8所述的基于多能源主体的动态碳排放因子测算方法,其特征在于,根据所述碳排放因子,可进一步计算碳排放总量、单位电量碳排放强度、碳源结构和能源碳排放结构。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310577830.XA CN116611706A (zh) | 2023-05-22 | 2023-05-22 | 基于多能源主体的动态碳排放因子测算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310577830.XA CN116611706A (zh) | 2023-05-22 | 2023-05-22 | 基于多能源主体的动态碳排放因子测算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116611706A true CN116611706A (zh) | 2023-08-18 |
Family
ID=87684865
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310577830.XA Pending CN116611706A (zh) | 2023-05-22 | 2023-05-22 | 基于多能源主体的动态碳排放因子测算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116611706A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116894354A (zh) * | 2023-09-11 | 2023-10-17 | 国网浙江省电力有限公司 | 一种多能源耦合环节的碳流计算方法及系统 |
-
2023
- 2023-05-22 CN CN202310577830.XA patent/CN116611706A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116894354A (zh) * | 2023-09-11 | 2023-10-17 | 国网浙江省电力有限公司 | 一种多能源耦合环节的碳流计算方法及系统 |
CN116894354B (zh) * | 2023-09-11 | 2023-12-26 | 国网浙江省电力有限公司 | 一种多能源耦合环节的碳流计算方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107730129B (zh) | 考虑光热热电联产与电锅炉的电-气-热互联系统风险评估方法 | |
CN105910169B (zh) | 基于机理模型预测控制的城市供热系统热网调节方法及系统 | |
CN110175311B (zh) | 一种基于多能耦合模型的优化潮流计算方法 | |
CN111815042B (zh) | 计及精细化热网模型的电热综合能源系统优化调度方法 | |
CN109492325B (zh) | 一种基于扩展能源集线器的多能耦合系统潮流分析方法 | |
CN109830955B (zh) | 考虑柔性约束与全周期成本的电-气配网柔性规划方法 | |
Riederer et al. | Development of a simulation platform for the evaluation of district energy system performances | |
CN107947182B (zh) | 一种动态潮流分析方法及动态潮流分析系统 | |
CN111428351B (zh) | 基于前推回代法的电-热综合能源系统潮流计算方法 | |
CN111082417A (zh) | 一种基于综合能源系统电气热联合网络的状态估计方法 | |
CN111046594B (zh) | 基于交叉迭代原理的热水供热网络动态仿真方法 | |
CN116611706A (zh) | 基于多能源主体的动态碳排放因子测算方法 | |
CN111463795B (zh) | 基于改进苏霍夫降温公式的电-热系统多能流分析方法 | |
CN111898816A (zh) | 一种综合能源系统动态状态估计方法 | |
Larsen et al. | Equivalent models for district heating systems | |
CN113806972B (zh) | 一种考虑供需双侧灵活性的综合能源系统可靠性分析方法 | |
CN112362096B (zh) | 多能流运行状态监测方法、装置、终端设备及存储介质 | |
CN111783309A (zh) | 基于内部守恒的蒸汽供热网络动态仿真方法 | |
CN115374527B (zh) | 一种区域电-热-冷综合能源动态模拟系统构建方法 | |
CN110020506A (zh) | 基于电热型综合能源系统运行优化的差分格式选择方法 | |
CN115062555A (zh) | 一种基于非平衡节点㶲的综合能源系统㶲流直接计算方法 | |
CN115081193A (zh) | 一种电网-热网模型构建方法、融合仿真方法及系统 | |
CN112906220B (zh) | 综合能源微网园区系统状态的估计方法 | |
CN113283077A (zh) | 考虑相关性的电-热-氢综合能源系统区间能流计算方法 | |
Wu et al. | Neural network based feasible region approximation model for optimal operation of integrated electricity and heating system |
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 |