CN107591800A - 基于场景分析的含分布式电源配电网运行状态预测方法 - Google Patents
基于场景分析的含分布式电源配电网运行状态预测方法 Download PDFInfo
- Publication number
- CN107591800A CN107591800A CN201710790471.0A CN201710790471A CN107591800A CN 107591800 A CN107591800 A CN 107591800A CN 201710790471 A CN201710790471 A CN 201710790471A CN 107591800 A CN107591800 A CN 107591800A
- Authority
- CN
- China
- Prior art keywords
- scene
- sequence
- mrow
- scenes
- fragment
- 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.)
- Granted
Links
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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/381—Dispersed generators
-
- 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
- H02J13/00—Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network
- H02J13/00002—Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network characterised by monitoring
-
- 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
-
- 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/001—Methods to deal with contingencies, e.g. abnormalities, faults or failures
-
- 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/001—Methods to deal with contingencies, e.g. abnormalities, faults or failures
- H02J3/0012—Contingency detection
-
- 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
-
- 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/466—Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
-
- 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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- 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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS 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
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/12—Monitoring or controlling equipment for energy generation units, e.g. distributed energy generation [DER] or load-side generation
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS 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
- Y04S40/00—Systems for electrical power generation, transmission, distribution or end-user application management characterised by the use of communication or information technologies, or communication or information technology specific aspects supporting them
- Y04S40/20—Information technology specific aspects, e.g. CAD, simulation, modelling, system security
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于场景分析的配电网运行状态预测方法,包括如下步骤:步骤10)获取配电系统网络架构、历史运行信息;步骤20)根据分布式电源历史出力序列,提取分布式电源出力代表性场景序列片段;步骤30)通过历史相似场景匹配,得到未来单时间断面T0多场景预测结果;步骤40)建立未来多时间断面运行场景树;步骤50)深度遍历未来多时间断面运行场景树中的各场景,针对各场景分别进行配电网潮流分析,计算配电网线路电流越限风险、母线电压越限风险,得到含分布式电源配电网未来运行状态变化趋势。按照该方法进行含分布式电源配电网运行状态预测,通过多场景分析,可以实时预测配电网运行状态的发展态势,及时做出风险预判。
Description
技术领域
本发明属于配电网态势感知领域,涉及一种配电网运行状态预测方法,更具体地,涉及一种基于场景分析的含分布式电源配电网运行状态预测方法。
背景技术
含分布式电源配电网态势感知是配电系统可靠、经济和安全运行的重要基础。对含分布式电源配电网进行运行状态预测是主动配电网态势感知工作中的核心环节。与传统配电网相比,含分布式电源配电网的典型特征之一是分布式电源的加入使电力系统的不确定性提高,因而计及不确定性的分布式电源出力预测技术是其中的关键。现有的分布式电源出力预测技术,无论是点值预测还是概率预测,其结果均未对分布式电源输出功率的时空关联特征进行描述,此外,概率法需要知道概率分布信息,当概率分布未知或难以用确定的概率分布描述时,概率预测结果会产生偏差。
场景分析是解决随机问题的一种有效方法,通过对可能出现的场景进行模拟将模型中的不确定性因素转变成多个确定性场景问题,降低建模及求解难度。与传统的分布式电源出力预测相比,相较时间序列预测得到的单一预测结果,构建场景树能够提供了多种预想场景;此外,采用场景分析法既可体现系统运行中的不确定性,同时可以反映系统运行的时序特性,将场景分析应用到含分布式电源配电网运行状态预测具有可实施性与有效性,能够充分利用分布式电源历史运行信息和实时运行信息,给配电网态势预测提供新的思路。
发明内容
技术问题:本发明提供一种基于场景分析的含分布式电源配电网运行状态预测方法,通过对分布式电源出力信息进行多时间断面的多场景预测,给出配电网未来两小时的运行状态变化趋势。
技术方案:本发明的基于场景分析的含分布式电源配电网运行状态预测方法,包括如下步骤:
步骤10)获取配电系统网络架构、历史运行信息,所述历史运行信息包括分布式电源历史出力序列、各负荷点历史需求信息;
步骤20)根据分布式电源历史出力序列,提取分布式电源出力代表性场景序列片段;
步骤30)通过计算分布式电源实时出力序列片段与代表性场景序列片段的动态弯曲时间距离进行历史相似场景匹配,得到未来单时间断面T0的多场景预测结果;
步骤40)根据所述未来单时间断面多场景预测结果,建立未来多时间断面运行场景树;
步骤50)深度遍历未来多时间断面运行场景树中的各场景,针对各场景分别进行配电网潮流分析,计算配电网线路电流越限风险、母线电压越限风险,得到连续时间断面下的线路电流和母线电压越限风险的变化趋势,即为含分布式电源配电网未来运行状态变化趋势。
进一步的,本发明方法中,步骤10)中通过遍历网络进行节点编号,获取各节点类型,分布式电源接入位置,即得到配电系统网络架构。
进一步的,本发明方法中,步骤20)的具体流程如下:
步骤201)根据配电网运行状态预测范围确定需要提取代表性场景序列片段的分布式电源历史出力序列片段,其长度记作L;确定需要的代表性场景序列片段个数M;
步骤202)在分布式电源历史出力序列中截取拟提取代表性场景序列片段的长度为L的时间序列片段,记其数量为N,构成场景集;
步骤203)根据下式计算场景集中每个场景序列片段的出现概率p(ci):
式中,ci表示场景集中的第i个场景序列片段,i为场景序列片段编号;
步骤204)对于每个场景序列片段ci,根据下式计算其与其他场景序列片段之间的Kantorovich距离,找到距离其最近的场景序列片段并在场景集中标记,形成最小场景距离矩阵KD,KD中对应场景序列片段ci的矩阵元素KD(i)根据下式计算:
KD(i)=min{||ci-cj||2,j∈[1,2,3,...N],j≠i},i∈[1,2,3,...N]
其中,cj表示场景集中的第j个场景序列片段,j为场景序列片段编号;
步骤205)对于每个场景序列片段ci,将其对应的最小场景距离与该场景序列片段的概率相乘,求得场景序列片段ci对应的最小场景概率距离,并找出场景集中最小概率距离最小的场景序列片段作为被剔除场景序列片段c*,将其从场景集中剔除,被剔除场景序列片段c*为:
c*=min{KD(i)*p(i)|i∈[1,2,3,...N]}
步骤206)寻找到距离被剔除场景序列片段c*最近的场景序列片段cn,根据下式更新cn的概率p(cn):
p(cn)=p(c*)+p(cn)
步骤207)令场景序列片段总数N=N-1,如果更新后的场景序列片段总数N=M,则结束步骤20),否则返回步骤204)。
进一步的,本发明方法中,步骤30)的具体流程如下:
步骤301)基于步骤20)提取的分布式电源出力序列代表性场景序列片段,计算分布式电源实时出力序列与第k个代表性场景序列片段的动态弯曲时间距离DTWk;
步骤302)取动态弯曲时间距离的倒数并对其进行归一化处理,得到分布式电源实时出力序列与代表性场景序列片段的相似度,将该相似度作为对应预测场景出现的概率,由第k个代表性场景序列及对应的动态弯曲时间距离DTWk计算分布式电源出力序列的未来预测值Fk,M个未来预测值组成未来单时间断面T0的多场景预测结果。
进一步的,本发明方法中,步骤40)的具体流程如下:
步骤401)将步骤30)产生的未来单时间断面T0的多场景预测结果纳入分布式电源实时出力序列,按照与步骤30)相同的方式得到下一时间断面T’=T0+Δt的多场景预测结果,共U=M2个,Δt为预测间隔;
步骤402)针对时间断面T’的多场景预测结果,进行场景削减,设定时间断面T’削减后的场景序列数M’,分别计算U个场景序列之间的Kantorovich距离,形成最小场景距离矩阵KD’,KD’中对应场景序列cs的矩阵元素KD’(s)根据下式计算:
KD′(s)=min{||cs-ct||2,t∈[1,2,3,...M2],t≠s},s∈[1,2,3,...M2]
其中,cs和ct分别表示包含时间断面T预测值F的分布式电源实时出力序列集中的第s和第t个场景序列,s和t为场景序列编号;
步骤403)对于每个场景序列cs,将其对应的最小场景距离与该场景序列的概率相乘,求得场景序列cs对应的最小场景概率距离,并找出场景集中最小概率距离最小的场景序列作为被剔除场景序列c^,将其从场景集中剔除,被剔除场景序列c^为:
c^=min{KD′(s)*p(s)|s∈[1,2,3,...M2]}
寻找到距离被剔除场景序列c^最近的场景序列cm,根据下式更新cm的概率p(cm):
p(cm)=p(c^)+p(cm)
步骤404)令场景总数U=U-1,如果更新后的场景总数U=M’,则进入步骤405),否则返回步骤402);
步骤405)如果T’=T0+n*Δt,则将各时间断面预测结果依时间先后排列,生成未来多时间断面运行场景树,并结束步骤40),否则令T=T’,T’=T+Δt,M=M’,返回步骤401),其中n为需要预测的时间断面个数。
进一步的,本发明方法中,步骤50)中的具体流程如下:
步骤501)深度遍历未来多时间断面运行场景树中的各场景,即在各场景下,将分布式电源预测出力值视为负值负荷,利用前推回代计算配电网潮流,得到各线路电流和母线电压情况;
步骤502)基于潮流计算结果,分别根据以下公式计算各场景下的线路过载值LOL、线路过载严重度SOL(C/E)、电压越限值LOV、母线过压严重度SOV(C/E),得到配电网线路电流越限风险OLR、母线电压越限风险OVR:
线路过载值LOL为:
LOL=L-0.8
其中,L表示流过线路的电流占其额定电流的比例;
线路过载严重度为:
线路电流越限风险OLR:
其中,NL为全网线路数。
电压越限值LOV为:
LOV=|1.05-V|
其中,V为节点电压标幺值。
母线过压严重度为:
母线电压越限风险OVR:
其中,NP为全网节点数;
步骤503)将步骤502)计算结果自时间断面T0至第nn个时间断面依次排列,得到连续时间断面下的线路电流和母线电压越限风险的变化趋势,即为含分布式电源配电网未来运行状态变化趋势。
有益效果:本发明与现有技术相比,具有以下优点:
本发明提出的场景分析方法,充分利用了分布式电源历史出力信息和实时出力信息,给出了未来两小时内的分布式电源出力超短期多场景预测结果,通过构建未来多时间断面运行场景树并对各单一场景进行潮流分析,提供了配电网运行状态的多种发展趋势。与时间序列单场景预测结果相比,本方面提出的方法关注了小概率场景发生的可能性以发生后的配电网运行状态变化趋势,有利于更全面地开展配电网态势感知与风险预警。
附图说明
图1是本发明实施例的方法流程示意图;
图2是接有分布式电源的IEEE-33节点配电系统结构图。
具体实施方式
如图1所示,本发明的一种基于场景分析的含分布式电源配电网运行状态预测方法,图2是接有分布式电源的IEEE-33节点配电系统,给定了网络中平衡节点的电压幅值及相角大小、PQ节点负荷大小、PV节点的电压幅值大小,接入该系统的分布式电源的历史出力信息已知(出力数据五分钟记录一次)。为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例对本发明进行深入地详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定发明。
步骤10)获取配电系统网架结构,遍历网络将节点编号,获取各节点类型,分布式电源接入位置,如图2所示;获取分布式电源历史出力序列、各负荷点历史需求信息。
步骤20)根据分布式电源历史出力序列,提取分布式电源出力代表性场景序列,具体步骤如下:
步骤201)现需预测未来两个小时的配电网运行状态,预测间隔为十五分钟。假设当前时刻为2017年6月1日上午12:00,需要提取代表性场景片段的分布式电源出力序列片段包括历史三年5月15日至6月18日上午10:05-下午14:00的出力信息,每个时间序列片段长度为48;确定需要的代表性场景片段为M=5个;
步骤202)在分布式电源历史出力序列中截取拟提取代表性场景片段的长度为48的时间序列片段,记其数量N=105,构成场景集;
步骤203)根据下式计算场景集中每个场景序列片段的出现概率p(ci):
式中,ci表示场景集中的第i个场景序列,i为场景序列编号。
步骤204)对于每个场景序列片段ci,根据下式计算其与其他场景序列片段之间的Kantorovich距离,找到距离其最近的场景序列片段并在场景集中标记,形成最小场景距离矩阵KD,KD中对应场景序列片段ci的矩阵元素KD(i):
KD(i)=min{||ci-cj||2,j∈[1,2,3,...N],j≠i},i∈[1,2,3,...N]
其中,cj表示场景集中的第j个场景序列片段,j为场景序列片段编号。
步骤205)对于每个场景序列片段ci,将其对应的最小场景距离与该场景序列片段的概率相乘,求得场景序列片段ci对应的最小场景概率距离,并找出场景集中最小概率距离最小的场景序列片段作为被剔除场景序列片段c*,将其从场景集中剔除,被剔除场景序列片段c*为:
c*=min{KD(i)*p(i)|i∈[1,2,3,...N]}
步骤206)寻找到距离被剔除场景序列片段c*最近的场景序列片段cn,根据下式更新cn的概率p(cn):
p(cn)=p(c*)+p(cn)
步骤207)令场景序列片段总数N=N-1,如果更新后的场景序列片段总数N=M,则结束步骤20),否则返回步骤204)。
步骤30),通过计算分布式电源实时出力序列与代表性场景的动态弯曲时间距离进行历史相似场景匹配,得到未来单时间断面多场景预测结果,具体步骤如下:
步骤301)基于步骤20)提取的5个分布式电源出力序列代表性场景序列片段,计算分布式电源实时出力序列R与第k个代表性场景序列片段Q的动态弯曲时间距离DTWk,具体计算方法如下:
第k个代表性场景序列片段Q长度l=24(只计算前10:05-12:00的时间序列片段),分布式电源实时出力序列R长度p=24,即T={t1,t2,…,tl},R={r1,r2,…,rp}。
构造24行24列的距离矩阵A,即
其中,f=2,3,…,24;g=2,3,…,24;D(24,24)为距离矩阵A的最小累加值,即为分布式电源实时出力序列R和第k个代表性场景序列片段Q的最短距离DTWk。
步骤302)取动态弯曲时间距离的倒数并对其进行归一化处理,得到分布式电源实时出力序列与代表性场景序列片段的相似度,将该相似度作为对应预测场景出现的概率,由第k个代表性场景序列及对应的动态弯曲时间距离DTWk计算分布式电源出力序列中12:15时刻的出力预测值Fk,M个未来预测值组成未来单时间断面(2017年6月1日12:15)多场景预测结果。
步骤40)中,根据多场景预测结果,建立未来多时间断面运行场景树,具体步骤如下:
步骤401)将步骤30)产生的未来单时间断面T=T0=2017年6月1日12:15的多场景预测结果(共5个场景)纳入分布式电源出力序列,重复步骤30)进行下一时间断面T’=2017年6月1日12:30的多场景预测工作,预测间隔Δt=15min;
步骤402)针对时间断面2017年6月1日12:30的多场景预测结果,进行场景削减,削减前共有U=M2=25个场景,设定削减后场景数为M’=5。分别计算25个场景序列之间的Kantorovich距离,形成最小场景距离矩阵KD’,KD’中对应场景序列cs的矩阵元素KD’(s)根据下式计算:
KD′(s)=min{||cs-ct||2,t∈[1,2,3,...25],t≠s},s∈[1,2,3,...25]
其中,cs和ct分别表示包含时间断面2017年6月1日12:30的多场景预测结果的分布式电源实时出力序列集中的第s和第t个场景序列,s和t为场景序列编号;。
步骤403)对于每个场景序列cs,将其对应的最小场景距离与该场景序列的概率相乘,求得场景序列cs对应的最小场景概率距离,并找出场景集中最小概率距离最小的场景序列作为被剔除场景序列c^,将其从场景集中剔除,被剔除场景序列c^为:
c^=min{KD′(s)*p(s)|s∈[1,2,3,...M2]}
寻找到距离被剔除场景序列c^最近的场景序列cm,根据下式更新cm的概率p(cm):
p(cm)=p(c^)+p(cm)
步骤404)令场景总数U=U-1,如果更新后的场景总数U=M’,则进入步骤405),否则返回步骤402);
步骤405)如果T’=T0+8*Δt,则将各时间断面预测结果依时间先后排列,生成未来多时间断面运行场景树,并结束步骤40),否则令T=T’,T’=T+Δt,M=M’,返回步骤401)。
步骤50)中,深度遍历未来多时间断面运行场景树中的各场景,针对各场景依次进行配电网潮流分析,计算配电网线路电流越限风险、母线电压越限风险,得到连续时间断面下的线路电流和母线电压越限风险的变化趋势,即为含分布式电源配电网未来运行状态变化趋势,具体步骤如下:
步骤501)深度遍历未来多时间断面运行场景树中的各场景,由未来多时间断面运行场景树最后一次预测生成的单时间断面多场景预测值为起点,依次寻找其父节点,即上一时刻的预测值,直至根节点,以此路径反向生成连续时间断面。
在各场景下,将分布式电源预测出力值视为负值负荷,利用前推回代计算配电网潮流,得到各线路电流和母线电压情况;
初始化:给定平衡节点电压,并为全网其他PQ节点赋电压初始值)PV节点赋无功注入初始功率Qi (0)。
计算各节点运算功率:
从网络末端开始逐步前推,由节点电压求全网各支路功率分布,前推过程:
从始端出发,逐段回推,由支路功率求各节点电压
利用求得的各节点电压修正PV节点电压和无功功率:
根据收敛判据检查是否收敛,若不满足收敛条件,将各节点电压计算值作为新的初始值代入式(2)开始进入下一次迭代。
步骤502)基于潮流计算结果,计算各场景下的线路过载值LOL、线路过载严重度SOL(C/E)、电压越限值LOV、母线过压严重度SOV(C/E),得到配电网线路电流越限风险OLR、母线电压越限风险OVR:
线路过载值LOL为:
LOL=L-0.8
其中,L表示流过线路的电流占其额定电流的比例。
上式反应了单一线路的过载值,在此基础上定义线路过载风险。定义设备过载风险严重度函数SOL(C/E)。设定流经线路的电流决定线路过载风险严重度。当线路电流小于或者等于额定电流的80%时,SOL(C/E)取为0;随着流过线路电流的增加,SOL(C/E)增大,且增加速率变快。
线路过载严重度为:
线路电流越限风险OLR:
其中,NL为全网线路数。
电压越限值LOV为:
LOV=|1.05-V|
其中,V为节点电压标幺值。
上式反应了单一母线的电压越限值,在此基础上定义电压过载风险,评估区域整体的母线过压风险水平。定义母线的电压越限风险严重度函数为SOV(C/E)。设定当母线电压为1.05p.u.时,严重度函数取值为0;随着电压越限值的增加,节点电压越限风险严重度也增加。
母线过压严重度为:
母线电压越限风险OVR:
其中,NP为全网节点数。
步骤503)将步骤502)计算结果自时间断面T0至第nn个时间断面依次排列,得到连续时间断面下的线路电流和母线电压越限风险的变化趋势,即为含分布式电源配电网未来运行状态变化趋势。
上述实施例仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和等同替换,这些对本发明权利要求进行改进和等同替换后的技术方案,均落入本发明的保护范围。
Claims (6)
1.一种基于场景分析的含分布式电源配电网运行状态预测方法,其特征在于,该方法包括如下步骤:
步骤10)获取配电系统网络架构、历史运行信息,所述历史运行信息包括分布式电源历史出力序列、各负荷点历史需求信息;
步骤20)根据分布式电源历史出力序列,提取分布式电源出力代表性场景序列片段;
步骤30)通过计算分布式电源实时出力序列片段与代表性场景序列片段的动态弯曲时间距离进行历史相似场景匹配,得到未来单时间断面T0的多场景预测结果;
步骤40)根据所述未来单时间断面多场景预测结果,建立未来多时间断面运行场景树;
步骤50)深度遍历未来多时间断面运行场景树中的各场景,针对各场景分别进行配电网潮流分析,计算配电网线路电流越限风险、母线电压越限风险,得到连续时间断面下的线路电流和母线电压越限风险的变化趋势,即为含分布式电源配电网未来运行状态变化趋势。
2.根据权利要求1所述的基于场景分析的含分布式电源配电网运行状态预测方法,其特征在于,所述步骤10)中通过遍历网络进行节点编号,获取各节点类型,分布式电源接入位置,即得到配电系统网络架构。
3.根据权利要求1所述的基于场景分析的含分布式电源配电网运行状态预测方法,其特征在于,所述步骤20)的具体流程如下:
步骤201)根据配电网运行状态预测范围确定需要提取代表性场景序列片段的分布式电源历史出力序列片段,其长度记作L;确定需要的代表性场景序列片段个数M;
步骤202)在分布式电源历史出力序列中截取拟提取代表性场景序列片段的长度为L的时间序列片段,记其数量为N,构成场景集;
步骤203)根据下式计算场景集中每个场景序列片段的出现概率p(ci):
<mrow>
<mi>p</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>c</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3</mn>
<mo>,</mo>
<mn>...</mn>
<mi>N</mi>
</mrow>
式中,ci表示场景集中的第i个场景序列片段,i为场景序列片段编号;
步骤204)对于每个场景序列片段ci,根据下式计算其与其他场景序列片段之间的Kantorovich距离,找到距离其最近的场景序列片段并在场景集中标记,形成最小场景距离矩阵KD,KD中对应场景序列片段ci的矩阵元素KD(i)根据下式计算:
KD(i)=min{||ci-cj||2,j∈[1,2,3,...N],j≠i},i∈[1,2,3,...N]
其中,cj表示场景集中的第j个场景序列片段,j为场景序列片段编号;
步骤205)对于每个场景序列片段ci,将其对应的最小场景距离与该场景序列片段的概率相乘,求得场景序列片段ci对应的最小场景概率距离,并找出场景集中最小概率距离最小的场景序列片段作为被剔除场景序列片段c*,将其从场景集中剔除,被剔除场景序列片段c*为:
c*=min{KD(i)*p(i)|i∈[1,2,3,...N]}
步骤206)寻找到距离被剔除场景序列片段c*最近的场景序列片段cn,根据下式更新cn的概率p(cn):
p(cn)=p(c*)+p(cn)
步骤207)令场景序列片段总数N=N-1,如果更新后的场景序列片段总数N=M,则结束步骤20),否则返回步骤204)。
4.根据权利要求1、2或3所述的基于场景分析的含分布式电源配电网运行状态预测方法,其特征在于,所述步骤30)的具体流程如下:
步骤301)基于步骤20)提取的分布式电源出力序列代表性场景序列片段,计算分布式电源实时出力序列与第k个代表性场景序列片段的动态弯曲时间距离DTWk;
步骤302)取动态弯曲时间距离的倒数并对其进行归一化处理,得到分布式电源实时出力序列与代表性场景序列片段的相似度,将该相似度作为对应预测场景出现的概率,由第k个代表性场景序列及对应的动态弯曲时间距离DTWk计算分布式电源出力序列的未来预测值Fk,M个未来预测值组成未来单时间断面T0的多场景预测结果。
5.根据权利要求1、2或3所述的基于场景分析的含分布式电源配电网运行状态预测方法,其特征在于,所述步骤40)的具体流程如下:
步骤401)将步骤30)产生的未来单时间断面T0的多场景预测结果纳入分布式电源实时出力序列,按照与步骤30)相同的方式得到下一时间断面T’=T0+Δt的多场景预测结果,共U=M2个,Δt为预测间隔;
步骤402)针对时间断面T’的多场景预测结果,进行场景削减,设定时间断面T’削减后的场景序列数M’,分别计算U个场景序列之间的Kantorovich距离,形成最小场景距离矩阵KD’,KD’中对应场景序列cs的矩阵元素KD’(s)根据下式计算:
KD′(s)=min{||cs-ct||2,t∈[1,2,3,...M2],t≠s},s∈[1,2,3,...M2]
其中,cs和ct分别表示包含时间断面T预测值F的分布式电源实时出力序列集中的第s和第t个场景序列,s和t为场景序列编号;
步骤403)对于每个场景序列cs,将其对应的最小场景距离与该场景序列的概率相乘,求得场景序列cs对应的最小场景概率距离,并找出场景集中最小概率距离最小的场景序列作为被剔除场景序列c^,将其从场景集中剔除,被剔除场景序列c^为:
c^=min{KD′(s)*p(s)|s∈[1,2,3,...M2]}
寻找到距离被剔除场景序列c^最近的场景序列cm,根据下式更新cm的概率p(cm):
p(cm)=p(c^)+p(cm)
步骤404)令场景总数U=U-1,如果更新后的场景总数U=M’,则进入步骤405),否则返回步骤402);
步骤405)如果T’=T0+n*Δt,则将各时间断面预测结果依时间先后排列,生成未来多时间断面运行场景树,并结束步骤40),否则令T=T’,T’=T+Δt,M=M’,返回步骤401),其中n为需要预测的时间断面个数。
6.根据权利要求1、2或3所述的基于场景分析的含分布式电源配电网运行状态预测方法,其特征在于,所述步骤50)中的具体流程如下:
步骤501)深度遍历未来多时间断面运行场景树中的各场景,即在各场景下,将分布式电源预测出力值视为负值负荷,利用前推回代计算配电网潮流,得到各线路电流和母线电压情况;
步骤502)基于潮流计算结果,分别根据以下公式计算各场景下的线路过载值LOL、线路过载严重度SOL(C/E)、电压越限值LOV、母线过压严重度SOV(C/E),得到配电网线路电流越限风险OLR、母线电压越限风险OVR:
线路过载值LOL为:
LOL=L-0.8
其中,L表示流过线路的电流占其额定电流的比例;
线路过载严重度为:
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>O</mi>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>/</mo>
<mi>E</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mi>e</mi>
<msub>
<mi>L</mi>
<mrow>
<mi>O</mi>
<mi>L</mi>
</mrow>
</msub>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
线路电流越限风险OLR:
<mrow>
<mi>O</mi>
<mi>L</mi>
<mi>R</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>N</mi>
<mi>L</mi>
</mrow>
</munderover>
<msub>
<mi>S</mi>
<mrow>
<mi>O</mi>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>/</mo>
<mi>E</mi>
<mo>)</mo>
</mrow>
</mrow>
其中,NL为全网线路数。
电压越限值LOV为:
LOV=|1.05-V|
其中,V为节点电压标幺值。
母线过压严重度为:
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>O</mi>
<mi>V</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>/</mo>
<mi>E</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mi>e</mi>
<msub>
<mi>L</mi>
<mrow>
<mi>O</mi>
<mi>V</mi>
</mrow>
</msub>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
母线电压越限风险OVR:
<mrow>
<mi>O</mi>
<mi>V</mi>
<mi>R</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>N</mi>
<mi>P</mi>
</mrow>
</munderover>
<msub>
<mi>S</mi>
<mrow>
<mi>O</mi>
<mi>V</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>/</mo>
<mi>E</mi>
<mo>)</mo>
</mrow>
</mrow>
其中,NP为全网节点数;
步骤503)将步骤502)计算结果自时间断面T0至第nn个时间断面依次排列,得到连续时间断面下的线路电流和母线电压越限风险的变化趋势,即为含分布式电源配电网未来运行状态变化趋势。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710790471.0A CN107591800B (zh) | 2017-09-04 | 2017-09-04 | 基于场景分析的含分布式电源配电网运行状态预测方法 |
US16/639,744 US20200212710A1 (en) | 2017-09-04 | 2018-04-27 | Method for predicting operation state of power distribution network with distributed generations based on scene analysis |
PCT/CN2018/084936 WO2019041857A1 (zh) | 2017-09-04 | 2018-04-27 | 一种基于场景分析的含分布式电源配电网运行状态预测方法 |
US17/978,149 US20230052730A1 (en) | 2017-09-04 | 2022-10-31 | Method for predicting operation state of power distribution network with distributed generations based on scene analysis |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710790471.0A CN107591800B (zh) | 2017-09-04 | 2017-09-04 | 基于场景分析的含分布式电源配电网运行状态预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107591800A true CN107591800A (zh) | 2018-01-16 |
CN107591800B CN107591800B (zh) | 2020-01-17 |
Family
ID=61051770
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710790471.0A Expired - Fee Related CN107591800B (zh) | 2017-09-04 | 2017-09-04 | 基于场景分析的含分布式电源配电网运行状态预测方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20200212710A1 (zh) |
CN (1) | CN107591800B (zh) |
WO (1) | WO2019041857A1 (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108376316A (zh) * | 2018-02-12 | 2018-08-07 | 国家电网公司 | 一种风电功率预测方法和系统 |
CN109301877A (zh) * | 2018-09-13 | 2019-02-01 | 国网江苏省电力有限公司 | 一种分布式电源及节点负荷典型运行场景集生成方法 |
WO2019041857A1 (zh) * | 2017-09-04 | 2019-03-07 | 东南大学 | 一种基于场景分析的含分布式电源配电网运行状态预测方法 |
CN110110815A (zh) * | 2019-05-22 | 2019-08-09 | 国网河北省电力有限公司 | 相似类型时间断面的确定方法、装置以及电子设备 |
CN114336792A (zh) * | 2022-02-14 | 2022-04-12 | 华北电力大学(保定) | 一种电网运行状态预测方法与系统 |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3751699B1 (de) * | 2019-06-13 | 2021-09-15 | Siemens Aktiengesellschaft | Verfahren und anordnung zur schätzung eines netzzustands eines energieverteilungsnetzes |
CN110266000B (zh) * | 2019-06-17 | 2022-07-19 | 国网江苏省电力有限公司 | 一种配电网电压越限原因分析方法、系统及存储介质 |
CN110599006B (zh) * | 2019-08-25 | 2022-08-12 | 南京理工大学 | 基于场景分析的配电网运行风险评估方法 |
CN110601204B (zh) * | 2019-10-14 | 2024-02-02 | 国网辽宁省电力有限公司盘锦供电公司 | 基于随机变量状态时序模拟的光伏并网系统概率潮流分析方法 |
CN112307677A (zh) * | 2020-11-05 | 2021-02-02 | 浙江大学 | 基于深度学习的电网振荡模态评估与安全主动预警方法 |
CN112488367A (zh) * | 2020-11-18 | 2021-03-12 | 国网山西省电力公司晋城供电公司 | 一种基于量子遗传的用户相序降损方法及其系统 |
CN112434874B (zh) * | 2020-12-02 | 2024-05-10 | 华北电力大学 | 一种可再生能源消纳的线路容量优化方法及系统 |
CN112736909B (zh) * | 2020-12-28 | 2023-07-18 | 智光研究院(广州)有限公司 | 储能系统实时控制方法、装置、电子设备和存储介质 |
CN113241793A (zh) * | 2021-05-27 | 2021-08-10 | 国网江苏省电力有限公司经济技术研究院 | 一种计及风电场景的含ipfc电力系统预防控制方法 |
CN113326897A (zh) * | 2021-06-25 | 2021-08-31 | 国网冀北电力有限公司承德供电公司 | 一种架空输电线路的测温计划生成方法、装置和电子设备 |
CN113536206B (zh) * | 2021-07-19 | 2023-12-01 | 国网陕西省电力公司 | 一种配电网区域预警方法、系统、终端设备及可读存储介质 |
CN114077921B (zh) * | 2021-10-15 | 2023-03-31 | 国电南瑞科技股份有限公司 | 变压器感知量趋势预测及状态逐级预警方法、装置及系统 |
CN114237183B (zh) * | 2021-12-20 | 2024-04-30 | 东北大学 | 考虑成品油随机需求的多周期生产计划方案的制定方法 |
CN114336608B (zh) * | 2021-12-30 | 2023-07-11 | 国网浙江省电力有限公司电力科学研究院 | 一种考虑动态增容和重构的机组阻塞优化方法及系统 |
CN116455766B (zh) * | 2023-06-13 | 2023-09-08 | 山东大学 | 基于信号序列分解的架空导线载流容量预测方法及系统 |
CN117614132A (zh) * | 2023-11-27 | 2024-02-27 | 广州航海学院 | 用于配电网的配电变压器电压越限画像方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101901428A (zh) * | 2010-07-21 | 2010-12-01 | 中国电力科学研究院 | 一种采用soa技术的电力市场仿真系统 |
CN103972985A (zh) * | 2014-05-26 | 2014-08-06 | 湖南大学 | 一种配电网的在线安全预警与预防控制方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106355511A (zh) * | 2015-07-22 | 2017-01-25 | 国网浙江省电力公司台州供电公司 | 考虑新能源与电动汽车接入的主动配电网重构方法 |
CN105741193B (zh) * | 2016-04-20 | 2019-10-18 | 河海大学 | 计及分布式电源和负荷不确定性的多目标配网重构方法 |
CN106230020B (zh) * | 2016-08-11 | 2018-10-19 | 浙江工业大学 | 一种微电网下考虑分布式电源消纳的电动汽车互动响应控制方法 |
CN107591800B (zh) * | 2017-09-04 | 2020-01-17 | 国网江苏省电力公司南京供电公司 | 基于场景分析的含分布式电源配电网运行状态预测方法 |
-
2017
- 2017-09-04 CN CN201710790471.0A patent/CN107591800B/zh not_active Expired - Fee Related
-
2018
- 2018-04-27 US US16/639,744 patent/US20200212710A1/en not_active Abandoned
- 2018-04-27 WO PCT/CN2018/084936 patent/WO2019041857A1/zh active Application Filing
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101901428A (zh) * | 2010-07-21 | 2010-12-01 | 中国电力科学研究院 | 一种采用soa技术的电力市场仿真系统 |
CN103972985A (zh) * | 2014-05-26 | 2014-08-06 | 湖南大学 | 一种配电网的在线安全预警与预防控制方法 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019041857A1 (zh) * | 2017-09-04 | 2019-03-07 | 东南大学 | 一种基于场景分析的含分布式电源配电网运行状态预测方法 |
CN108376316A (zh) * | 2018-02-12 | 2018-08-07 | 国家电网公司 | 一种风电功率预测方法和系统 |
CN108376316B (zh) * | 2018-02-12 | 2020-10-27 | 国家电网公司 | 一种风电功率预测方法和系统 |
CN109301877A (zh) * | 2018-09-13 | 2019-02-01 | 国网江苏省电力有限公司 | 一种分布式电源及节点负荷典型运行场景集生成方法 |
CN109301877B (zh) * | 2018-09-13 | 2021-08-17 | 国网江苏省电力有限公司 | 一种分布式电源及节点负荷典型运行场景集生成方法 |
CN110110815A (zh) * | 2019-05-22 | 2019-08-09 | 国网河北省电力有限公司 | 相似类型时间断面的确定方法、装置以及电子设备 |
CN114336792A (zh) * | 2022-02-14 | 2022-04-12 | 华北电力大学(保定) | 一种电网运行状态预测方法与系统 |
Also Published As
Publication number | Publication date |
---|---|
US20200212710A1 (en) | 2020-07-02 |
WO2019041857A1 (zh) | 2019-03-07 |
CN107591800B (zh) | 2020-01-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107591800A (zh) | 基于场景分析的含分布式电源配电网运行状态预测方法 | |
CN106972481B (zh) | 规模化充电设施接入有源配电网的安全性量化评估方法 | |
CN104332996B (zh) | 一种评估电力系统可靠性的方法 | |
CN105449713B (zh) | 考虑分布式电源特性的有源配电网智能软开关规划方法 | |
CN105633948B (zh) | 一种分布式能源接入电力系统随机模糊潮流算法 | |
CN105656031B (zh) | 基于高斯混合分布特征的含风电电力系统安全风险评估方法 | |
CN103426120A (zh) | 基于可靠性的中低压配电网的综合评估方法 | |
CN105207197B (zh) | 包含风电场的电力系统可靠性评估方法 | |
CN104852399B (zh) | 光储微网系统的储能容量动态优化方法 | |
CN109102146A (zh) | 基于多参数线性规划的电力系统风险评估加速方法 | |
CN104657796A (zh) | 配电网规划方法 | |
CN103824234A (zh) | 一种基于分块和层级结构的配电系统可靠性评估方法 | |
CN105005708A (zh) | 一种基于ap聚类算法的广义负荷特性聚类方法 | |
CN105305488A (zh) | 一种考虑新能源并网对输电网利用率影响的评价方法 | |
CN113569411A (zh) | 一种面向灾害天气的电网运行风险态势感知方法 | |
CN103347028A (zh) | 云架构下基于贝叶斯的对等网络信任度量模型 | |
CN110348540A (zh) | 基于聚类的电力系统暂态功角稳定故障筛选方法及装置 | |
CN103514365B (zh) | 电力系统暂态相继失稳模式搜索及关联交流断面识别方法 | |
CN108074048A (zh) | 计入风速相关特性的风电并网电力系统安全风险评估方法 | |
CN105787604A (zh) | 一种输配协同的配网态势快速感知方法 | |
CN105205572A (zh) | 基于EMD分解和Elman算法的光伏出力预测方法 | |
CN107908731A (zh) | 基于psd‑bpa广西电网网区负荷数据批量修改方法 | |
CN105207255B (zh) | 一种适用于风电出力的电力系统调峰计算方法 | |
Han et al. | An assessment approach of the power system vulnerability considering the uncertainties of wind power integration | |
AU2021106383A4 (en) | Method for Calculating Reliability of Distribution System Based on Equipment Reliability Cloud Model |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200117 Termination date: 20200904 |
|
CF01 | Termination of patent right due to non-payment of annual fee |