CN108039737B - 一种源网荷协调运行模拟系统 - Google Patents

一种源网荷协调运行模拟系统 Download PDF

Info

Publication number
CN108039737B
CN108039737B CN201711472764.0A CN201711472764A CN108039737B CN 108039737 B CN108039737 B CN 108039737B CN 201711472764 A CN201711472764 A CN 201711472764A CN 108039737 B CN108039737 B CN 108039737B
Authority
CN
China
Prior art keywords
unit
output
simulation
power
wind
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
CN201711472764.0A
Other languages
English (en)
Other versions
CN108039737A (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.)
State Grid Corp of China SGCC
State Grid Jiangsu Electric Power Co Ltd
State Grid Energy Research Institute Co Ltd
Economic and Technological Research Institute of State Grid Jiangsu Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Jiangsu Electric Power Co Ltd
State Grid Energy Research Institute Co Ltd
Economic and Technological Research Institute of State Grid Jiangsu 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 State Grid Corp of China SGCC, State Grid Jiangsu Electric Power Co Ltd, State Grid Energy Research Institute Co Ltd, Economic and Technological Research Institute of State Grid Jiangsu Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201711472764.0A priority Critical patent/CN108039737B/zh
Publication of CN108039737A publication Critical patent/CN108039737A/zh
Application granted granted Critical
Publication of CN108039737B publication Critical patent/CN108039737B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/12Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
    • H02J3/14Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by switching loads on to, or off from, network, e.g. progressively balanced loading
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/28Arrangements for balancing of the load in a network by storage of energy
    • H02J3/386
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/003Load forecast, e.g. methods or systems for forecasting future load demand
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02BCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
    • Y02B70/00Technologies for an efficient end-user side electric power management and consumption
    • Y02B70/30Systems integrating technologies related to power network operation and communication or information technologies for improving the carbon footprint of the management of residential or tertiary loads, i.e. smart grids as climate change mitigation technology in the buildings sector, including also the last stages of power distribution and the control, monitoring or operating management systems at local level
    • Y02B70/3225Demand response systems, e.g. load shedding, peak shaving
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E70/00Other energy conversion or management systems reducing GHG emissions
    • Y02E70/30Systems combining energy storage with energy generation of non-fossil origin
    • YGENERAL 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S20/00Management or operation of end-user stationary applications or the last stages of power distribution; Controlling, monitoring or operating thereof
    • Y04S20/20End-user application control systems
    • Y04S20/222Demand response systems, e.g. load shedding, peak shaving

Abstract

本发明公开了一种源网荷协调运行模拟系统,其特征在于,包括数据输入模块:用于输入系统模拟的相关数据;运行计算模块:用于“源‑网‑荷”协调运行过程中的运行模拟与计算;结果输出模块:用于输出运行计算模块计算分析获得的数据及信息。本发明提供的一种源网荷协调运行模拟系统考虑了影响风速和风电场出力的各因素,使风电场模拟更加科学,分别对风电场出力进行模拟和对“源‑网‑荷”进行协调运行模拟,使运营模拟结果更加精确,实现电网能够提高大规模风电的接入能力、降低系统消纳新能源的成本,提高可再生能源利用率,为我国制定受端地区电力发展规划、风电优化布局具有重大的现实意义。

Description

一种源网荷协调运行模拟系统
技术领域
本发明涉及电力系统的技术领域,具体涉及一种源网荷协调运行模拟系统。
背景技术
目前我国风电正蓬勃发展,尤其在受端地区大规模风电接入使电力系统面临新的挑战、孕育着重大的变革。一方面,风电的发展将使电力系统更加低碳、经济、高效,另一方面,风电的间歇性使系统运行方式复杂多变,有可能使系统在短时间内经历多种运行方式,引起系统潮流的大幅度变化,增加系统的运行风险。
针对风电出力随机波动性以及大规模风电并网后电力系统运行复杂问题,杜红卫、鲁文等人在2016年提出“城市配电网源网荷互动优化调度技术研究与应用”,该方案主要集中在配电侧,建立了分布式电源与配电网、配电网与上级电源、配电网与多样性负荷之间的优化调度模型;另外,李逐云、雷霞等人在2017年2月提出“考虑“源-网-荷”三方利益的主动配电网协调规划”,该方案重在处理各个间歇性电源出力波动在区域电网中形成的“多点扰动”问题;陈艳伟在2015年提出“基于博弈理论的大电网源网荷适应性评价模型及算法”,该方案在“源-网-荷”三侧中引入博弈理论,其结果可以指导电网规划建设。
因此如何科学准确地模拟风电出力以及接入后电力系统运行,使未来电网能够提高大规模风电的接入能力、降低系统消纳新能源的成本,为我国制定受端地区电力发展规划、风电优化布局具有重大的现实意义。
发明内容
为了解决上述技术问题,本发明所采用的技术方案是提供一种源网荷协调运行模拟系统,包括
数据输入模块:用于输入系统模拟的相关数据;其中,相关数据包括负荷预测信息、联络线送受电计划、电源规划信息、机组运行技术参数、电网规划信息、线路与断面信息、风电场风资源信息与系统运行参数设定;
运行计算模块:用于“源-网-荷”协调运行过程中的运行模拟与计算;
结果输出模块:用于输出运行计算模块计算分析获得的数据及信息;
运行计算模块包括检修计划安排单元:用于根据系统负荷与系统检修规则安排机组检修,其检修结果作为电力系统运行模拟单元的边界条件;
可再生能源运行模拟单元:用于根据可再生能源资源的随机特性以及装机情况,模拟可再生能源的随机运行方式;
多日运行协调单元:用于根据所述检修计划安排单元确定的机组检修计划,所述可再生能源运行模拟单元确定的可再生能源随机运行方式,以及从所述结果统计分析单元获取的历史运行模拟结果的信息,对电力系统运行模拟的参数进行闭环调整;
电力系统运行模拟单元根据检修计划安排单元、可再生能源运行模拟单元以及多日运行协调单元传输至的信息进行运行模拟,并将运行模拟结果传输至结果统计分析单元;
结果统计分析单元:用于根据电力系统运行模拟单元模拟得到的各类信息,统计各类指标,并将结果反馈至多日运行协调单元与结果输出模块。
在上述方案中,所述可再生能源运行模拟单元包括风速模拟子单元、风速修正子单元、风电场模拟出力子单元与风电场修正出力子单元;
风速模拟子单元根据风速统计特征参数模拟风电场模拟风速曲线;
风速修正子单元根据风速曲线对模拟风速进行修正,其中,风速曲线包括但不限于日平均风速曲线、季平均风速曲线;
风电场模拟出力子单元根据修正的风速曲线和风电机组出力特性曲线模拟风电场模拟出力曲线;
风电场修正出力子单元根据风机出力影响因素对风电场出力进行修正,其中,风机出力影响因素包括风机可靠性和风电场风速尾流效应。
在上述方案中,所述风速模拟子单元实现以下计算过程:
计及概率分布、时序自相关系数,采用随机差分方程方法生成风电场时序出力公式如下:
Figure RE-GDA0001609461900000031
式中,v^it*为风电场t时刻模拟风速;Xt为t时刻的随机过程,计算公式如下:
Figure RE-GDA0001609461900000032
式中,θ为风电场风速序列自相关系数,f(x)是符合Weibull分布的X的概率分布,x为随机变量,c为尺度参数,Γ为伽马函数,k为形状参数。
在上述方案中,所述风速修正子单元实现以下计算过程:
考虑风电场风速的季节性与日内的规律性,对随机生成的风速序列v^it* 进行如下修正:
Figure RE-GDA0001609461900000033
式中,vit*为修正的风速曲线,kim为风电场风速季节因子,kih为风电场日内各小时平均风速曲线。
在上述方案中,所述风电场模拟出力子单元实现以下计算过程:
考虑风电机组出力特性曲线,生成风电场模拟出力曲线如下:
Figure RE-GDA0001609461900000041
式中,Ci为风电机组出力特性曲线,其计算公式如下:
Figure RE-GDA0001609461900000042
其中,v、vin、vrated与vout分别为风电机组风速、切入风速、额定风速与切出风速,R为额定出力。
在上述方案中,所述风电场修正出力子单元实现以下计算过程:
考虑风电场尾流效应与风电场风电机组出力可靠性,风电场修正出力曲线如下:
Figure RE-GDA0001609461900000043
其中,nit为风电场尾流效应系数;η为一随机变量,代表风电场内机组可靠性水平。
在上述方案中,所述电力系统运行模拟单元包括确定日可运行机组子单元与修正多区负荷曲线子单元;
确定日可运行机组子单元具体实现以下内容:
根据装机进度表,考虑机组的投运、退役和技改,确定投运机组;
根据检修计划排除检修机组,确定可运行机组及其参数;
修正多区负荷曲线子单元具体实现以下内容:
安排所有可以确定出力的机组,包括外来的协议送电、核电机组以及认为指定出力的机组,根据电源所在区域,修正对应的负荷曲线;
根据可再生能源运行模拟模块随机模拟生成的可再生能源模拟出力,安排新可再生能源出力,修正对应的负荷曲线;
根据原始负荷曲线和正负备用率确定每个时段系统需要的正负备用。
在上述方案中,所述电力系统运行模拟单元包括对剩余机组进行优化模拟运行子单元;
优化模拟运行子单元具体实现以下内容:
对于抽水蓄能和常规水电机组安排其削峰和填谷,抽水可设定为平抽或满抽方式,并满足机组的容量、电量约束,对柔性负荷,根据经济性目标确立是否切除可中断负荷;
对可中断负荷,根据常规机组优化运行结果安排可中断负荷和可再生能源弃风弃光。
在上述方案中,所述常规机组优化运行过程如下:
a1、以系统发电成本最低为目标,包括机组发电成本、切新能源成本、切负荷成本和机组启停成本,确定目标函数,具体如下:
Figure RE-GDA0001609461900000051
式中:T为优化周期内时段总数;Pc,t、Pf,t、Ph,t、Pp,t、Pw,t分别为t时段日内不可启停火电、日内可启停火电、水电、抽蓄与新能源机组的输出功率; Cc(Pc,t)、Cf(Pf,t)、Ch(Ph,t)、Cp(Pp,t)、Cw(Pw,t)为t时段不可启停火电、日内可启停火电、水电、抽蓄与新能源机组的运行成本;Pwd,t为t时段切新能源功率;Dd,t为t时段切负荷功率;Vd为各节点切负荷损失;Cf、Cc为机组启、停费用;θ、η、γ为加权系数,通常情况下为1,或根据需要调整;
a2、以电力平衡约束、火电机组技术出力约束、新能源机组出力约束、水电及抽蓄出力约束、系统正负备用约束、分区备用约束、网络约束、跨区域送电约束、动态约束等作为系统约束。各约束条件如下:
电力平衡约束如下:
Figure RE-GDA0001609461900000061
式中,Z为分区总数,
Figure RE-GDA0001609461900000062
为第z个分区内的电源数量,
Figure RE-GDA0001609461900000063
为第z个分区里第g个电源的出力,
Figure RE-GDA0001609461900000064
为t时刻接受外来电功率,ΔPz loss_t为t时刻分区z 的网损,
Figure RE-GDA0001609461900000065
为t时刻分区z的负荷;
火电机组的技术出力约束主要包括火电机组出力上下限以及火电机组爬坡约束,其约束如下:
Figure RE-GDA0001609461900000066
式中:Pcmin、Pcmax、Pfmin、Pfmax分别为机组最小出力与最大出力;Ic为日内不可启停机组状态变量,If,t为日内可启停机组t时段状态变量;ΔPf,down、ΔPf,up、ΔPc,down、ΔPc,up分别为机组下爬坡与上爬坡速率;
新能源机组要求,在调度运行中,仅根据日前预测值控制出力在预测值或者以下的范围内运行;同时根据电网运行节能、环保的要求,在电网能够接纳的情况下要求新能源全额上网,其约束如下:
Figure RE-GDA0001609461900000067
式中,Pw,t为t时段新能源出力;Pwd,t为t时段切新能源功率;Pwf,t为 t时段新能源预测出力值;
水电机组根据中长期跨流域梯级水电优化调度结果给出的水电机组出力范围以及日发电量在模型中优化分配各时段出力;抽蓄机组考虑每日的抽水量与发电量平衡;
Figure RE-GDA0001609461900000071
式中:Phmin为水电机组最小出力;Phmax为水电机组最大出力,Qhydro为水电机组日发电量,Pp,pump、Pp,gen为抽蓄机组单位时段最大抽水量与发电量, Ip,pump,t、Ip,gen,t分别为抽蓄机组t时段抽水与发电的状态变量,λp为抽蓄机组的效率;
系统正负备用约束分为系统正备用约束与系统负备用约束,具体如下:
系统正备用约束:
Figure RE-GDA0001609461900000072
式中:ru,t为系统在时段t要求的正备用率,在式中,可再生能源对系统备用的贡献应按其预测出力Pwf,t计算,即使被切除,其被切除部分也计入备用容量;
系统正负备用约束:
Figure RE-GDA0001609461900000073
式中:rd,t为系统在时段t要求的负备用率,此时,可再生能源出力在此不计入,即等效的认为可再生能源最小出力为0;
不同区域需要保证各自区域内开机为区域内负荷预留充足的备用容量,分区备用约束可表示为:
Figure RE-GDA0001609461900000081
式中:Z为区域总数;Dz,t为z区域t时段的负荷;ru,z,t、rd,z,t分别为z区域t时段的正备用率与负备用率;Z+和Z分别为z区域送入与送出的联络线集合;fl,t为l线路的潮流;
网络约束基于直流潮流模型建立线路以及断面传输约束,如下式所示:
Figure RE-GDA0001609461900000082
式中:Fl,t、Fs,t分别为t时段线路、断面潮流矩阵;W为发电机转移分布因子;Angc、Angf、Angh、Angp、Angw分别为不同类型机组的节点-机组关联矩阵;Asl为断面-线路关联矩阵;
为模拟实际电网运行中特定的送电计划,引入跨区域送电约束,即限定断面各时段的输送潮流;
Figure RE-GDA0001609461900000083
式中,Θs为跨区送电断面s包含的线路集合;
动态约束要求启停火电在一天中运行状态最多启停一次,大容量火电机组不允许启停。
在上述方案中,所述结果输出模块输出由结果统计分析单元输至的信息,该信息包括发电负荷平衡信息、风电弃风电量、各类机组出力信息、系统运行成本、系统购电费用、环保与排放指标、线路与断面潮流、风电可信容量。
本发明提供的一种“源-网-荷”协调运行模拟系统考虑了影响风速和风电场出力的各因素,使风电场模拟更加科学,分别对风电场出力进行模拟和对“源-网-荷”进行协调运行模拟,使运营模拟结果更加精确,实现电网能够提高大规模风电的接入能力、降低系统消纳新能源的成本,提高可再生能源利用率,为我国制定受端地区电力发展规划、风电优化布局具有重大的现实意义。
附图说明
图1为本发明的提供的“源-网-荷”协调运行模拟系统的结构框图;
图2为本发明系统中可再生能源运行模拟单元的结构框图;
图3为本发明系统中电力系统运行模拟单元的运行示意图。
具体实施方式
本发明提供的一种源网荷协调运行模拟系统考虑了影响风速和风电场出力的各因素,使风电场模拟更加科学,分别对风电场出力进行模拟和对“源-网-荷”进行协调运行模拟,使运营模拟结果更加精确,实现电网能够提高大规模风电的接入能力、降低系统消纳新能源的成本,提高可再生能源利用率,为我国制定受端地区电力发展规划、风电优化布局具有重大的现实意义。下面结合附图对本发明作出详细的说明。
如图1所示,本发明提供了一种源网荷协调运行模拟系统,包括数据输入模块、运行计算模块与结果输出模块。
数据输入模块:用于输入系统模拟的相关数据;其中,相关数据负荷预测信息、联络线送受电计划、电源规划信息、机组运行技术参数、电网规划信息、线路与断面信息、风电场风资源信息与系统运行参数设定等。
运行计算模块:用于“源-网-荷”协调运行过程中的运行模拟与计算。
结果输出模块:用于输出运行计算模块计算分析获得的数据及信息。
运行计算模块包括:检修计划安排单元、可再生能源运行模拟单元、电力系统运行模拟单元、结果统计分析单元、多日运行协调单元。电力系统运行模拟单元根据检修计划安排单元、可再生能源运行模拟单元以及多日运行协调单元传输至的信息进行运行模拟,并将运行模拟结果传输至结果统计分析单元。
检修计划安排单元:用于根据系统负荷与系统检修规则安排机组检修,其检修结果作为电力系统运行模拟单元的边界条件。
可再生能源运行模拟单元:用于根据可再生能源资源的随机特性以及装机情况,模拟可再生能源的随机运行方式。
如图2所示,可再生能源运行模拟单元包括风速模拟子单元、风速修正子单元、风电场模拟出力子单元与风电场修正出力子单元。其中,
风速模拟子单元根据风速统计特征参数模拟风电场模拟风速曲线,其中,风速统计特征参数为实际风速统计得到的参数,包括但不限于概率分布、时序自相关系数、尺度参数、形状参数、风速曲线等。风速模拟子单元实现以下计算过程:
计及概率分布、时序自相关系数,采用随机差分方程方法生成风电场时序出力公式如下:
Figure RE-GDA0001609461900000101
式中,v^it*为风电场t时刻模拟风速;Xt为t时刻的随机过程,计算公式如下:
Figure RE-GDA0001609461900000111
式中,θ为风电场风速序列自相关系数,f(x)是符合Weibull分布的X的概率分布,x为随机变量,c为尺度参数,Γ为伽马函数,k为形状参数。
风速修正子单元根据风速曲线对模拟风速进行修正,其中,风速曲线包括但不限于日平均风速曲线、季平均风速曲线。风速修正子单元实现以下计算过程:
考虑风电场风速的季节性与日内的规律性,对随机生成的风速序列v^it* 进行如下修正:
Figure RE-GDA0001609461900000112
式中,vit*为修正的风速曲线,kim为风电场风速季节因子,kih为风电场日内各小时平均风速曲线。
风电场模拟出力子单元根据修正的风速曲线和风电机组出力特性曲线模拟风电场模拟出力曲线。风电场模拟出力子单元实现以下计算过程:
考虑风电机组出力特性曲线,生成风电场模拟出力曲线如下:
Figure RE-GDA0001609461900000113
式中,Ci为风电机组出力特性曲线,其计算公式如下:
Figure RE-GDA0001609461900000114
其中,v、vin、vrated与vout分别为风电机组风速、切入风速、额定风速与切出风速,R为额定出力。
风电场修正出力子单元根据风机出力影响因素对风电场出力进行修正,其中,风机出力影响因素包括但不限于风机可靠性和风电场风速尾流效应。风电场修正出力子单元实现以下计算过程:
考虑风电场尾流效应与风电场风电机组出力可靠性,风电场修正出力曲线如下:
Figure RE-GDA0001609461900000121
其中,nit为风电场尾流效应系数;η为一随机变量,代表风电场内机组可靠性水平。
电力系统运行模拟单元:用于根据检修计划安排单元确认的边界条件与机组检修计划以及可再生能源运行模拟单元确定的可再生能源随机运行方式,利用考虑网络约束的机组组合和经济调度模型,以日为单位模拟系统的调度运行,获得机组出力、线路潮流等信息。
如图3所示,电力系统运行模拟单元包括确定日可运行机组子单元、修正多区负荷曲线子单元及对剩余机组进行优化模拟运行子单元几部分构成,其具体情况如下:
电力系统运行模拟主要由确定日可运行机组、修正多区负荷曲线和确定备用、对剩余机组和可中断负荷进行优化模拟运行、机组优化运行几部分构成。
确定日可运行机组子单元具体实现以下内容:
1)、根据装机进度表,考虑机组的投运、退役和技改等,确定投运机组;
2)、根据检修计划排除检修机组,确定可运行机组及其参数;
修正多区负荷曲线子单元具体实现以下内容:
1)、安排所有可以确定出力的机组,包括外来的协议送电、核电机组以及认为指定出力的机组,根据电源所在区域,修正对应的负荷曲线;
2)、根据可再生能源运行模拟模块随机模拟生成的可再生能源模拟出力,安排新可再生能源出力,修正对应的负荷曲线;
3)、根据原始负荷曲线和正负备用率确定每个时段系统需要的正负备用。
优化模拟运行子单元具体实现以下内容:
1)、对于抽水蓄能和常规水电机组安排其削峰和填谷,抽水可设定为平抽或满抽方式,并满足机组的容量、电量等约束,对柔性负荷,根据经济性目标确立是否切除可中断负荷。
2)、对可中断负荷,根据常规机组优化运行结果安排可中断负荷和可再生能源弃风弃光。其中,常规机组优化运行过程如下:
a1、以系统发电成本最低为目标,包括机组发电成本、切新能源成本、切负荷成本和机组启停成本,确定目标函数,具体如下:
Figure RE-GDA0001609461900000131
式中:T为优化周期内时段总数;Pc,t、Pf,t、Ph,t、Pp,t、Pw,t分别为t时段日内不可启停火电、日内可启停火电、水电、抽蓄与新能源机组的输出功率; Cc(Pc,t)、Cf(Pf,t)、Ch(Ph,t)、Cp(Pp,t)、Cw(Pw,t)为t时段不可启停火电、日内可启停火电、水电、抽蓄与新能源机组的运行成本;Pwd,t为t时段切新能源功率;Dd,t为t时段切负荷功率;Vd为各节点切负荷损失;Cf、Cc为机组启、停费用;θ、η、γ为加权系数,通常情况下为1,或根据需要调整。
a2、以电力平衡约束、火电机组技术出力约束、新能源机组出力约束、水电及抽蓄出力约束、系统正负备用约束、分区备用约束、网络约束、跨区域送电约束、动态约束等作为系统约束。各约束条件如下:
系统及分区电力平衡约束如下:
Figure RE-GDA0001609461900000132
式中,Z为分区总数,
Figure RE-GDA0001609461900000133
为第z个分区内的电源数量,
Figure RE-GDA0001609461900000134
为第z个分区里第g个电源的出力,
Figure RE-GDA0001609461900000136
为t时刻接受外来电功率,ΔPz loss_t为t时刻分区z 的网损,
Figure RE-GDA0001609461900000135
为t时刻分区z的负荷。
火电机组的技术出力约束主要包括火电机组出力上下限以及火电机组爬坡约束,其约束如下:
Figure RE-GDA0001609461900000141
式中:Pcmin、Pcmax、Pfmin、Pfmax分别为机组最小出力与最大出力;Ic为日内不可启停机组状态变量,If,t为日内可启停机组t时段状态变量;ΔPf,down、ΔPf,up、ΔPc,down、ΔPc,up分别为机组下爬坡与上爬坡速率。
新能源机组要求,在调度运行中,仅根据日前预测值控制出力在预测值或者以下的范围内运行;同时根据电网运行节能、环保的要求,在电网能够接纳的情况下要求新能源全额上网。本模型在日运行模拟模型中引入切除新能源机制,使模型在系统无法提供调峰容量,系统备用容量不足或新能源送出受阻情况下切除部分可再生能源出力。
Figure RE-GDA0001609461900000142
式中,Pw,t为t时段新能源出力;Pwd,t为t时段切新能源功率;Pwf,t为t 时段新能源预测出力值。
水电机组根据中长期跨流域梯级水电优化调度结果给出的水电机组出力范围以及日发电量在模型中优化分配各时段出力;抽蓄机组考虑每日的抽水量与发电量平衡。
Figure RE-GDA0001609461900000151
式中:Phmin为水电机组最小出力;Phmax为水电机组最大出力,Qhydro为水电机组日发电量,Pp,pump、Pp,gen为抽蓄机组单位时段最大抽水量与发电量, Ip,pump,t、Ip,gen,t分别为抽蓄机组t时段抽水与发电的状态变量,λp为抽蓄机组的效率。
系统正负备用约束分为系统正备用约束与系统负备用约束,具体如下:
系统正备用约束如下:
Figure RE-GDA0001609461900000152
式中:ru,t为系统在时段t要求的正备用率,需要注意的是,在式中,可再生能源对系统备用的贡献应按其预测出力Pwf,t计算,即使被切除,其被切除部分也计入备用容量。
系统正负备用约束如下:
Figure RE-GDA0001609461900000153
式中:rd,t为系统在时段t要求的负备用率,此时,可再生能源出力在此不计入,即等效的认为可再生能源最小出力为0,即可以随时被切除。
不同区域需要保证各自区域内开机为区域内负荷预留充足的备用容量,分区备用约束可表示为:
Figure RE-GDA0001609461900000154
式中:Z为区域总数;Dz,t为z区域t时段的负荷;ru,z,t、rd,z,t分别为z区域t时段的正备用率与负备用率;Z+和Z分别为z区域送入与送出的联络线集合;fl,t为l线路的潮流。
网络约束基于直流潮流模型建立线路以及断面传输约束,如下式所示:
Figure RE-GDA0001609461900000161
式中:Fl,t、Fs,t分别为t时段线路、断面潮流矩阵;W为发电机转移分布因子;Angc、Angf、Angh、Angp、Angw分别为不同类型机组的节点-机组关联矩阵;Asl为断面-线路关联矩阵。
为模拟实际电网运行中特定的送电计划,引入跨区域送电约束,即限定断面各时段的输送潮流。
Figure RE-GDA0001609461900000162
式中,Θs为跨区送电断面s包含的线路集合。
动态约束要求启停火电在一天中运行状态最多启停一次,大容量火电机组不允许启停。
结果统计分析单元:用于根据电力系统运行模拟单元模拟得到的各类机组出力、线路与断面潮流等信息,统计机组发电量、系统运行成本、系统购电费用、环保与排放等指标,并将结果反馈至多日运行协调单元与结果输出模块。
多日运行协调单元:用于根据检修计划安排单元确定的机组检修计划,可再生能源运行模拟单元确定的可再生能源随机运行方式,以及从结果统计分析单元获取的历史运行模拟结果等信息,对电力系统运行模拟的参数进行闭环调整。
结果输出模块输出由结果统计分析单元输至的信息,该信息包括发电负荷平衡信息、风电弃风电量、各类机组出力信息、系统运行成本、系统购电费用、环保与排放指标、线路与断面潮流、风电可信容量等信息。
本系统的有益效果入下:
1)对“源-网-荷”进行协调运行模拟,能够提高可再生能源利用率,提高经济型。
2)系统对风电场出力模型进行详细模拟,使运营模拟结果更加精确。
3)综合考虑了影响风速和风电场出力的各因素,使风电场模拟更加科学。
4)综合统筹了各类电源,使系统更具适应性。
5)在模拟过程中考虑了人工指定因素、自然出力因素以及电网约束和经济性,是模拟效果及结果更加符合实际。
最后应说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (10)

1.一种源网荷协调运行模拟系统,其特征在于,包括
数据输入模块:用于输入系统模拟的相关数据;其中,相关数据包括负荷预测信息、联络线送受电计划、电源规划信息、机组运行技术参数、电网规划信息、线路与断面信息、风电场风资源信息与系统运行参数设定;
运行计算模块:用于“源-网-荷”协调运行过程中的运行模拟与计算;包括检修计划安排单元、可再生能源运行模拟单元、多日运行协调单元、电力系统运行模拟单元、结果统计分析单元;
结果输出模块:用于输出运行计算模块计算分析获得的数据及信息;
检修计划安排单元:用于根据系统负荷与系统检修规则安排机组检修,其检修结果作为电力系统运行模拟单元的边界条件;
可再生能源运行模拟单元:用于根据可再生能源资源的随机特性以及装机情况,模拟可再生能源的随机运行方式;
多日运行协调单元:用于根据所述检修计划安排单元确定的机组检修计划,所述可再生能源运行模拟单元确定的可再生能源随机运行方式,以及从所述结果统计分析单元获取的历史运行模拟结果的信息,对电力系统运行模拟的参数进行闭环调整;
电力系统运行模拟单元根据检修计划安排单元、可再生能源运行模拟单元以及多日运行协调单元传输至的信息进行运行模拟,并将运行模拟结果传输至结果统计分析单元;
结果统计分析单元:用于根据电力系统运行模拟单元模拟得到的各类信息,统计各类指标,并将结果反馈至多日运行协调单元与结果输出模块。
2.如权利要求1所述的模拟系统,其特征在于,所述可再生能源运行模拟单元包括风速模拟子单元、风速修正子单元、风电场模拟出力子单元与风电场修正出力子单元;
风速模拟子单元根据风速统计特征参数模拟风电场模拟风速曲线;
风速修正子单元根据风速曲线对模拟风速进行修正,其中,风速曲线包括但不限于日平均风速曲线、季平均风速曲线;
风电场模拟出力子单元根据修正的风速曲线和风电机组出力特性曲线模拟风电场模拟出力曲线;
风电场修正出力子单元根据风机出力影响因素对风电场出力进行修正,其中,风机出力影响因素包括风机可靠性和风电场风速尾流效应。
3.如权利要求2所述的模拟系统,其特征在于,所述风速模拟子单元实现以下计算过程:
计及概率分布、时序自相关系数,采用随机差分方程方法生成风电场时序出力公式如下:
Figure FDA0002582183430000021
式中,v^it*为风电场t时刻模拟风速;Xt为t时刻的随机过程,计算公式如下:
Figure FDA0002582183430000022
式中,θ为风电场风速序列自相关系数,f(x)是符合Weibull分布的X的概率分布,x为随机变量,c为尺度参数,Γ为伽马函数,k为形状参数,Wt为标准布朗运动。
4.如权利要求2所述的模拟系统,其特征在于,所述风速修正子单元实现以下计算过程:
考虑风电场风速的季节性与日内的规律性,对随机生成的风速序列v^it*进行如下修正:
Figure FDA0002582183430000031
式中,vit*为修正的风速曲线,kim为风电场风速季节因子,kih为风电场日内各小时平均风速曲线。
5.如权利要求2所述的模拟系统,其特征在于,所述风电场模拟出力子单元实现以下计算过程:
考虑风电机组出力特性曲线,生成风电场模拟出力曲线如下:
Figure FDA0002582183430000032
式中,Ci为风电机组出力特性曲线,其计算公式如下:
Figure FDA0002582183430000033
其中,v、vin、vrated与vout分别为风电机组风速、切入风速、额定风速与切出风速,R为额定出力。
6.如权利要求2所述的模拟系统,其特征在于,所述风电场修正出力子单元实现以下计算过程:
考虑风电场尾流效应与风电场风电机组出力可靠性,风电场修正出力曲线如下:
Figure FDA0002582183430000034
其中,ηit为风电场尾流效应系数;η为一随机变量,代表风电场内机组可靠性水平。
7.如权利要求1所述的模拟系统,其特征在于,所述电力系统运行模拟单元包括确定日可运行机组子单元与修正多区负荷曲线子单元;
确定日可运行机组子单元具体实现以下内容:
根据装机进度表,考虑机组的投运、退役和技改,确定投运机组;
根据检修计划排除检修机组,确定可运行机组及其参数;
修正多区负荷曲线子单元具体实现以下内容:
安排所有可以确定出力的机组,包括外来的协议送电、核电机组以及认为指定出力的机组,根据电源所在区域,修正对应的负荷曲线;
根据可再生能源运行模拟模块随机模拟生成的可再生能源模拟出力,安排新可再生能源出力,修正对应的负荷曲线;
根据原始负荷曲线和正负备用率确定每个时段系统需要的正负备用。
8.如权利要求1所述的模拟系统,其特征在于,所述电力系统运行模拟单元包括对剩余机组进行优化模拟运行子单元;
优化模拟运行子单元具体实现以下内容:
对于抽水蓄能和常规水电机组安排其削峰和填谷,抽水可设定为平抽或满抽方式,并满足机组的容量、电量约束,对柔性负荷,根据经济性目标确立是否切除可中断负荷;
对可中断负荷,根据常规机组优化运行结果安排可中断负荷和可再生能源弃风弃光。
9.如权利要求8所述的模拟系统,其特征在于,所述常规机组优化运行过程如下:
a1、以系统发电成本最低为目标,包括机组发电成本、切新能源成本、切负荷成本和机组启停成本,确定目标函数,具体如下:
Figure FDA0002582183430000041
式中:T为优化周期内时段总数;Pc,t、Pf,t、Ph,t、Pp,t、Pw,t分别为t时段日内不可启停火电、日内可启停火电、水电、抽蓄与新能源机组的输出功率;Cc(Pc,t)、Cf(Pf,t)、Ch(Ph,t)、Cp(Pp,t)、Cw(Pw,t)为t时段不可启停火电、日内可启停火电、水电、抽蓄与新能源机组的运行成本;Pwd,t为t时段切新能源功率;Dd,t为t时段切负荷功率;Vd为各节点切负荷损失;Cf、Cc为机组启、停费用;θ、η、γ为加权系数,为1;
a2、以电力平衡约束、火电机组技术出力约束、新能源机组出力约束、水电及抽蓄出力约束、系统正负备用约束、分区备用约束、网络约束、跨区域送电约束、动态约束等作为系统约束;各约束条件如下:
电力平衡约束如下:
Figure FDA0002582183430000051
式中,Z为分区总数,
Figure FDA0002582183430000052
为第z个分区内的电源数量,Pgt z为第z个分区里第g个电源的出力,Ot z为t时刻接受外来电功率,ΔPz loss_t为t时刻分区z的网损,Lt z为t时刻分区z的负荷;
火电机组的技术出力约束主要包括火电机组出力上下限以及火电机组爬坡约束,其约束如下:
Figure FDA0002582183430000053
式中:Pcmin、Pcmax、Pfmin、Pfmax分别为机组最小出力与最大出力;Ic为日内不可启停机组状态变量,If,t为日内可启停机组t时段状态变量;ΔPf,down、ΔPf,up、ΔPc,down、ΔPc,up分别为机组下爬坡与上爬坡速率;
新能源机组要求,在调度运行中,仅根据日前预测值控制出力在预测值或者以下的范围内运行;同时根据电网运行节能、环保的要求,在电网能够接纳的情况下要求新能源全额上网,其约束如下:
Figure FDA0002582183430000061
式中,Pw,t为t时段新能源出力;Pwd,t为t时段切新能源功率;Pwf,t为t时段新能源预测出力值;
水电机组根据中长期跨流域梯级水电优化调度结果给出的水电机组出力范围以及日发电量在模型中优化分配各时段出力;抽蓄机组考虑每日的抽水量与发电量平衡;
Figure FDA0002582183430000062
式中:Phmin为水电机组最小出力;Phmax为水电机组最大出力,Qhydro为水电机组日发电量,Pp,pump、Pp,gen为抽蓄机组单位时段最大抽水量与发电量,Ip,pump,t、Ip,gen,t分别为抽蓄机组t时段抽水与发电的状态变量,λp为抽蓄机组的效率;
系统正负备用约束分为系统正备用约束与系统负备用约束,具体如下:
系统正备用约束:
Figure FDA0002582183430000063
式中:ru,t为系统在时段t要求的正备用率,在式中,可再生能源对系统备用的贡献应按其预测出力Pwf,t计算,即使被切除,其被切除部分也计入备用容量;
系统正负备用约束:
Figure FDA0002582183430000071
式中:rd,t为系统在时段t要求的负备用率,此时,可再生能源出力在此不计入,即等效的认为可再生能源最小出力为0;
不同区域需要保证各自区域内开机为区域内负荷预留充足的备用容量,分区备用约束可表示为:
Figure FDA0002582183430000072
式中:Z为区域总数;Dz,t为z区域t时段的负荷;ru,z,t、rd,z,t分别为z区域t时段的正备用率与负备用率;Z+和Z分别为z区域送入与送出的联络线集合;fl,t为l线路的潮流;
网络约束基于直流潮流模型建立线路以及断面传输约束,如下式所示:
Figure FDA0002582183430000073
式中:Fl,t、Fs,t分别为t时段线路、断面潮流矩阵;W为发电机转移分布因子;Angc、Angf、Angh、Angp、Angw分别为不同类型机组的节点-机组关联矩阵;Asl为断面-线路关联矩阵;
为模拟实际电网运行中特定的送电计划,引入跨区域送电约束,即限定断面各时段的输送潮流;
Figure FDA0002582183430000081
式中,Θs为跨区送电断面s包含的线路集合;
动态约束要求启停火电在一天中运行状态最多启停一次,大容量火电机组不允许启停。
10.如权利要求1所述的模拟系统,其特征在于,所述结果输出模块输出由结果统计分析单元输至的信息,该信息包括发电负荷平衡信息、风电弃风电量、各类机组出力信息、系统运行成本、系统购电费用、环保与排放指标、线路与断面潮流、风电可信容量。
CN201711472764.0A 2017-12-29 2017-12-29 一种源网荷协调运行模拟系统 Active CN108039737B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711472764.0A CN108039737B (zh) 2017-12-29 2017-12-29 一种源网荷协调运行模拟系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711472764.0A CN108039737B (zh) 2017-12-29 2017-12-29 一种源网荷协调运行模拟系统

Publications (2)

Publication Number Publication Date
CN108039737A CN108039737A (zh) 2018-05-15
CN108039737B true CN108039737B (zh) 2020-09-11

Family

ID=62098726

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711472764.0A Active CN108039737B (zh) 2017-12-29 2017-12-29 一种源网荷协调运行模拟系统

Country Status (1)

Country Link
CN (1) CN108039737B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109412162B (zh) * 2018-08-27 2020-11-10 国家电网公司 基于外网运行约束等值的电力系统最优潮流计算方法
CN109636041B (zh) * 2018-12-14 2023-04-07 国网山东省电力公司经济技术研究院 一种适于大规模复杂电网的机组检修模拟方法及系统
CN109494730A (zh) * 2018-12-14 2019-03-19 国网山东省电力公司经济技术研究院 新能源并网下电力系统逐日运行模拟仿真方法
CN109449988A (zh) * 2018-12-14 2019-03-08 国网山东省电力公司经济技术研究院 大规模新能源发电并网的电力系统逐日仿真模拟方法
CN109449987A (zh) * 2018-12-14 2019-03-08 国网山东省电力公司经济技术研究院 新能源发电并网的电力系统模拟方法及系统
CN109787296B (zh) * 2019-03-18 2022-06-14 长沙理工大学 一种考虑源荷协整关系的经济调控方法
CN111709112B (zh) * 2020-04-30 2023-05-16 广东电网有限责任公司电网规划研究中心 一种海上风电运行模拟方法、装置及存储介质
CN113283702B (zh) * 2021-04-22 2022-08-02 国电南瑞科技股份有限公司 融合电力系统安全稳定特性的潮流状态评价方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105071389A (zh) * 2015-08-19 2015-11-18 华北电力大学(保定) 计及源网荷互动的交直流混合微电网优化运行方法及装置
CN106549378A (zh) * 2016-12-09 2017-03-29 国网江苏省电力公司金湖县供电公司 一种针对分布式电源出力不确定性的配网协同调度方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10592833B2 (en) * 2016-04-01 2020-03-17 Enel X North America, Inc. Extended control in control systems and methods for economical optimization of an electrical system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105071389A (zh) * 2015-08-19 2015-11-18 华北电力大学(保定) 计及源网荷互动的交直流混合微电网优化运行方法及装置
CN106549378A (zh) * 2016-12-09 2017-03-29 国网江苏省电力公司金湖县供电公司 一种针对分布式电源出力不确定性的配网协同调度方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
源荷储一体化的广义灵活电源双层统筹规划;李海波等;《电力系统自动化》;20171110;第41卷(第21期);第46-54、104页 *
适应城市配电网源网荷三重不对称的智能配电网运行仿真平台研究与应用;马洲俊等;《供用电 》;20170605;第49-55页 *

Also Published As

Publication number Publication date
CN108039737A (zh) 2018-05-15

Similar Documents

Publication Publication Date Title
CN108039737B (zh) 一种源网荷协调运行模拟系统
Liu et al. Optimal operation of independent regional power grid with multiple wind-solar-hydro-battery power
CN103151803B (zh) 一种含风电系统机组及备用配置的优化方法
CN102694391B (zh) 风光储联合发电系统日前优化调度方法
CN106099993B (zh) 一种适应新能源大规模接入的电源规划方法
CN109449993B (zh) 一种面向多能互补电力系统的多时间尺度生产模拟方法
CN108092324B (zh) 一种风电参与调峰调频的agc控制系统和控制方法
CN107276122B (zh) 适应大规模可再生能源并网的调峰资源调用决策方法
CN104915737A (zh) 基于风险控制含大规模风电的电力系统协调经济调度方法
CN108808740B (zh) 多区域互联电网的多时间尺度的调度方法、装置以及存储介质
CN103699941A (zh) 一种电力系统调度运行年方案的制定方法
CN111342500B (zh) 小水电虚拟电厂的多时间尺度优化调度方法
CN110957717A (zh) 一种多电源电力系统多目标日前优化调度方法
CN104143839B (zh) 基于功率预测的风电场集群限出力有功功率分配方法
US20220029424A1 (en) Hybrid power plant
CN108493992A (zh) 一种含分布式潮流控制器的风电场优化调度方法
CN109146320A (zh) 一种考虑配电网安全性的虚拟电厂优化调度方法
US20240006891A1 (en) Two-stage self-organizing optimized aggregation method and system for distributed resources of virtual power plant (vpp)
CN114744687A (zh) 一种虚拟电厂的能源调控方法及系统
CN102856931A (zh) 风电场有功功率的动态分群控制方法
CN112234604A (zh) 一种多能互补电源基地优化配置方法、存储介质及设备
CN110867907B (zh) 一种基于多类型发电资源同质化的电力系统调度方法
CN108667077A (zh) 一种风蓄联合系统优化调度方法
CN105305485A (zh) 一种消纳大规模间歇性能源的安全约束经济调度方法
CN107769266A (zh) 一种多时间尺度发电与备用联合优化方法

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