CN104036118A - 一种电力系统并行化轨迹灵敏度获取方法 - Google Patents
一种电力系统并行化轨迹灵敏度获取方法 Download PDFInfo
- Publication number
- CN104036118A CN104036118A CN201410203928.XA CN201410203928A CN104036118A CN 104036118 A CN104036118 A CN 104036118A CN 201410203928 A CN201410203928 A CN 201410203928A CN 104036118 A CN104036118 A CN 104036118A
- Authority
- CN
- China
- Prior art keywords
- partiald
- differential algebraic
- sensitivity
- algebraic equations
- equations
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种电力系统并行化轨迹灵敏度获取方法。轨迹灵敏度在电力系统动态性能评估、稳定校核与控制等领域有着重要的应用价值,但在大规模电力系统分析中,存在着计算时间过长等技术困难,无法在以实时控制为基础的现有电力系统中获得有效应用。针对上述问题,本发明提出了应用并行计算技术加速电力系统轨迹灵敏度新型获取方法。通过挖掘和利用前向灵敏度分析中的牛顿迭代的矩阵结构性质,设计了其专有的不诚实牛顿迭代策略,使得牛顿迭代中的线性方程组能够得到解耦并被部署在不同的计算处理单元上完成计算,从而实现了轨迹灵敏度获取的并行化,极大地提高了其计算效率,可用于电力系统稳定分析与控制等领域。
Description
技术领域
本发明属于电力系统的运行、分析与调度技术领域,尤其涉及一种基于并行计算技术的电力系统轨迹灵敏度获取方法。
背景技术
随着我国坚强智能电网建设进程的推进,以三华电网为核心的特高压电网逐渐形成,大量可再生能源并网发电,电力系统动态性能在电网安全、稳定运行中的地位愈加重要。电力系统的动态性能通常通过电力系统时域仿真进行评估和验证,然而由于现代电力系统存在模型参数不准确、运行点经常变化的特性,一般需要大量的时域仿真方能有效的对电网的动态性能进行准确评估。文献《Sensitivity,Approximation,and Uncertainty in Power System DynamicSimulation》系统化的提出了轨迹灵敏度方法,并将其运用于电力系统仿真分析中。通过使用一阶轨迹灵敏度信息,可以有效的对参数空间中基准轨迹的领域内的轨迹簇进行估计,从而替代了计算复杂、费时的多模型参数时域仿真,大大减少了电力系统动态性能评估的计算量。所获得的轨迹灵敏度也为电力系统稳定控制提供了量化的依据,通过分析对失稳现象影响大的轨迹灵敏度,可以以此设计对应的电力系统稳定控制策略,从而提升系统运行的安全性。
由此可见,轨迹灵敏度是电力系统动态性能分析中的有力工具,但在实际工程应用中存在着轨迹灵敏度计算量大,求解时间长等技术困难。并行计算技术能够利用大规模并行计算机高效求解科学计算问题,突破现有串行计算的计算性能与规模的瓶颈。但轨迹灵敏度的计算本质上属于一个顺序串行计算的过程,无法直接应用并行计算技术予以加速求解。因此需要一种新型的轨迹灵敏度并行分解策略,有效的使用并行计算技术提高其计算效率。
发明内容
本发明的目的在于提出一种基于并行计算技术的电力系统轨迹灵敏度获取方法,该方法可提高电力系统动态性能评估的计算性能,可用于电力系统稳定分析与控制等应用领域。
电力系统并行化轨迹灵敏度获取方法包括如下步骤:
第一步:通过电力系统遥测系统获得电网运行静态数据,从稳定数据库中获得系统动态数据,构建参数向量p下的微分代数方程组D0(p)用以表征电力系统的动态行为。
第二步:对参数向量p中的每个分量构造对应的前向灵敏度微分代数方程组,如针对p中的第i个分量pi,对应的前向灵敏度微分代数方程为Di(p)。
第三步:使用隐式梯形积分法求解拓展微分代数方程组{D0(p),D1(p)…DNP(p)},其中NP为参数向量p的维度。差分化后拓展微分代数方程组的解通过一组非线性方程组E(p)描述。
第四步:使用不诚实牛顿法将非线性方程组E(p)转化为一系列的线性方程组L(p)的迭代格式。忽略雅可比矩阵中的非对角块,从而形成块对角结构。进而将线性方程组L(p)解耦为NP+1个独立线性方程组,记其中第i个线性方程组为Li(p),对每个独立线性方程组Li(p)使用独立的计算处理单元求解,实现拓展微分代数方程组求解的并行化。
第五步:微分代数方程D0(p)的解即为电力系统的时域仿真结果,前向灵敏度微分代数方程组Di(p)的解即为第i个参数pi对应的轨迹灵敏度。检测时域仿真结果中是否存在电网失稳的现象,如有电网失稳,则首先寻找轨迹灵敏度中对电网失稳影响较大的参数,然后基于其灵敏度设计电网稳定控制策略,使电力系统的运行规避失稳风险。
所述的第三步中的拓展微分代数方程组差分得到的第n个仿真时步下的非线性方程组E(p)为:
E(p):
其中h为仿真步长。拓展微分代数方程组和对应的拓展状态变量定义如下:拓展微分代数方程组为:
其中,y和分别为原微分代数方程组D0(p)的状态变量及其导数,f为原微分代数方程组D0(p)的定义函数,和分别为拓展微分代数方程组的拓展状态变量及其导数,si和分别为参数向量p的第i个分量对应的轨迹灵敏度及其导数。相应的,拓展微分代数方程组{D0(p),D1(p)…DNP(p)}中的各分量定义为:
D0(p):y(t0)=y0(p)
Di(p):
其中t0为仿真时段的初始时刻,y0为状态变量的初值。
所述的第四步中第n个仿真时步第j次牛顿迭代的L(p)具有下述形式:
其中和分别为第n个仿真时步第j次牛顿迭代的拓展状态变量及其导数。J为精确雅可比矩阵,具有下述矩阵结构:
J0和Ji(其中i取1…NP)为雅可比矩阵中的对角线和边界分块。
将所有矩阵块Ji忽略,从而得到不精确雅可比矩阵
于是L(p)即可解耦为NP+1个独立线性方程组L0(p)和Li(p)(其中i取1…NP),即:
L0(p):
Li(p):
所述的第五步中计算处理单元是指CPU核心、计算服务器节点或其他具有完整的浮点计算与逻辑处理能力的计算设备。
本发明的有益效果是,本发明提出了一种电力系统并行化轨迹灵敏度获取方法。该方法基于灵敏度分析原理、微分代数方程理论和不诚实牛顿法算法,利用轨迹灵敏度计算中特有的矩阵和方程结构,利用并行计算加速其求解过程。与已有的技术相比,本发明提出的方法主要有以下改进:
1、所提出的并行计算方式属于松耦合并行,具有加速比高、并行效率高的技术优势;
2、由于不诚实牛顿法的特性,本方法的计算结果与串行计算相同,即所有的并行化措施均基于数学等价变换;
3、本方法可使用于多种并行计算设备,对处理单元之间的通信条件要求低,具有大规模并行化的技术潜力。
附图说明
图1是电力系统并行化轨迹灵敏度获取方法流程图;
图2是电力系统发电机功角时域仿真曲线与轨迹灵敏度估计曲线;
图3是电力系统节点电压时域仿真曲线与轨迹灵敏度估计曲线;
图4是CASE678算例并行加速比曲线。
具体实施方式
电力系统并行化轨迹灵敏度获取方法包括如下步骤:
第一步:通过电力系统遥测系统获得电网运行静态数据,从稳定数据库中获得系统动态数据,构建参数向量p下的微分代数方程组D0(p)用以表征电力系统的动态行为。
第二步:对参数向量p中的每个分量构造对应的前向灵敏度微分代数方程组,如针对p中的第i个分量pi,对应的前向灵敏度微分代数方程为Di(p)。
第三步:使用隐式梯形积分法求解拓展微分代数方程组{D0(p),D1(p)…DNP(p)},其中NP为参数向量p的维度。差分化后拓展微分代数方程组的解通过一组非线性方程组E(p)描述。
第四步:使用不诚实牛顿法将非线性方程组E(p)转化为一系列的线性方程组L(p)的迭代格式。忽略雅可比矩阵中的非对角块,从而形成块对角结构。进而将线性方程组L(p)解耦为NP+1个独立线性方程组,记其中第i个线性方程组为Li(p),对每个独立线性方程组Li(p)使用独立的计算处理单元求解,实现拓展微分代数方程组求解的并行化。
第五步:微分代数方程D0(p)的解即为电力系统的时域仿真结果,前向灵敏度微分代数方程组Di(p)的解即为第i个参数pi对应的轨迹灵敏度。检测时域仿真结果中是否存在电网失稳的现象,如有电网失稳,则首先寻找轨迹灵敏度中对电网失稳影响较大的参数,然后基于其灵敏度设计电网稳定控制策略,使电力系统的运行规避失稳风险。
所述的第三步中的拓展微分代数方程组差分得到的第n个仿真时步下的非线性方程组E(p)为:
E(p):
其中h为仿真步长。拓展微分代数方程组和对应的拓展状态变量定义如下:拓展微分代数方程组为:
其中,y和分别为原微分代数方程组D0(p)的状态变量及其导数,f为原微分代数方程组D0(p)的定义函数,和分别为拓展微分代数方程组的拓展状态变量及其导数,si和分别为参数向量p的第i个分量对应的轨迹灵敏度及其导数。相应的,拓展微分代数方程组{D0(p),D1(p)…DNP(p)}中的各分量定义为:
D0(p):y(t0)=y0(p)
Di(p):
其中t0为仿真时段的初始时刻,y0为状态变量的初值。
所述的第四步中第n个仿真时步第j次牛顿迭代的L(p)具有下述形式:
其中和分别为第n个仿真时步第j次牛顿迭代的拓展状态变量及其导数。J为精确雅可比矩阵,具有下述矩阵结构:
J0和Ji(其中i取1…NP)为雅可比矩阵中的对角线和边界分块。
将所有矩阵块Ji忽略,从而得到不精确雅可比矩阵
于是L(p)即可解耦为NP+1个独立线性方程组L0(p)和Li(p)(其中i取1…NP),即:
L0(p):
Li(p):
所述的第五步中计算处理单元是指CPU核心、计算服务器节点或其他具有完整的浮点计算与逻辑处理能力的计算设备。
以下结合附图,对本发明的实施例作详细说明,该发明的流程图如图1所示。
实施例:
电力系统动态可以使用下述微分代数方程组进行描述:
其中,p为参数向量,f为微分代数方程组,状态变量x包括:
x=[Ux Uy x1 … xn] \*MERGEFORMAT(2)
f=[fN fC1 … fCn] \*MERGEFORMAT(3)
Ux和Uy为系统节点电压的实部与虚部,fN为系统节点电压方程,fCi为第i个系统动态元件的动态方程。
G(t)和B(t)为系统节点导纳矩阵的实部与虚部。
设定待分析的参数向量p为x0,即微分代数方程的初值。此时轨迹灵敏度分析将用于评估该动力系统初值对时域中系统动态的影响。
为了测试本发明的有效性,以下两个实施例被用于验证本方法的准确性和高效性。所述算法使用C++开发实现,并在一台配有8个IntelXeonE8-8837CPU的并行计算服务器上进行测试和验证。第三方软件包IDAS、ADC04和KLU被分别用于实现时域仿真、自动微分和独立线性方程组求解的功能。表1展示了所使用的测试电力系统算例信息。
表1:示例电力系统的系统参数
实施例1:
一个简化的九节点电力系统CASE9算例被用于验证该算法的有效性和计算结果的准确性。轨迹灵敏度的应用之一就是利用基准轨迹估计参数空间内邻域的轨迹,即:
x(p0+Δp,t)≈x(p0,t)+Δp·s(p0,t) \*MERGEFORMAT(5)
其中s(p0,t)为基准参数向量p0下的轨迹灵敏度,Δp为参数空间下的扰动,x(p0,t)为基准轨迹。受扰后的轨迹x(p0+Δp,t)可以用上述信息进行估计。对应该测试系统,Δp被设定为某一台发电机的出力减少10MW。图2和图3展示了受扰后轨迹(包括发电机功角曲线和系统节点电压曲线)的两种计算方法,即使用完整的时域仿真和使用轨迹灵敏度估计。图中的计算结果显示,所计算的轨迹灵敏度能够精确地对受扰后的时域轨迹进行估计,与实际时域仿真的曲线非常接近,误差可以忽略不计,从而验证了本发明提出的轨迹灵敏度计算的有效性。
实施例2:
一个标准电力系统测试系统CASE678被用于测试本方法的计算效率和并行加速比。为了验证本方法在大规模计算算例中的计算效率,微分代数方程组的所有初值均被选为待分析的参数向量,从而增加了测试使用的计算量,更加真实地模拟了轨迹灵敏度在实际电力系统动态性能评估中的应用场景。
表2给出了在轨迹灵敏度分析中的两个关键计算环节的计算时间:微分代数方程组的函数调用环节(简称函数调用),和微分代数方程组求解的差分化与牛顿迭代环节(简称时域仿真)。同时,为了测试该并行计算方法在大规模并行计算平台上的可拓展性,不同的CPU核心数量被应用于相同的算例中,由此可以观察到计算时间随CPU核心数增加而减少的趋势。
表2:CASE678测试算例的轨迹灵敏度并行计算时间
将表2中的计算性能数据进一步进行可视化处理,可以得到并行加速比,如图4所示。图4展示了对于轨迹灵敏度计算的两个关键步骤(函数调用和时域仿真)的并行加速比。本发明提出的方法均可以对其进行有效的并行加速,并在使用32个CPU核心时获得最大并行加速性能。如果使用更大规模的计算算例,该并行加速性能将会得到更大的提升,加速饱和现象也将出现在更多的CPU核心情景下。此实施例证明了本发明所提出方法的高效性。
Claims (4)
1.一种电力系统并行化轨迹灵敏度获取方法,其特征在于,包括如下步骤:
第一步:通过电力系统遥测系统获得电网运行静态数据,从稳定数据库中获得系统动态数据,构建参数向量p下的微分代数方程组D0(p)用以表征电力系统的动态行为。
第二步:对参数向量p中的每个分量构造对应的前向灵敏度微分代数方程组,如针对p中的第i个分量pi,对应的前向灵敏度微分代数方程为Di(p)。
第三步:使用隐式梯形积分法求解拓展微分代数方程组{D0(p),D1(p)…DNP(p)},其中NP为参数向量p的维度。差分化后拓展微分代数方程组的解通过一组非线性方程组E(p)描述。
第四步:使用不诚实牛顿法将非线性方程组E(p)转化为一系列的线性方程组L(p)的迭代格式。忽略雅可比矩阵中的非对角块,从而形成块对角结构。进而将线性方程组L(p)解耦为NP+1个独立线性方程组,记其中第i个线性方程组为Li(p),对每个独立线性方程组Li(p)使用独立的计算处理单元求解,实现拓展微分代数方程组求解的并行化。
第五步:微分代数方程D0(p)的解即为电力系统的时域仿真结果,前向灵敏度微分代数方程组Di(p)的解即为第i个参数pi对应的轨迹灵敏度。检测时域仿真结果中是否存在电网失稳的现象,如有电网失稳,则首先寻找轨迹灵敏度中对电网失稳影响较大的参数,然后基于其灵敏度设计电网稳定控制策略,使电力系统的运行规避失稳风险。
2.根据权利要求1所述的电力系统并行化轨迹灵敏度获取方法,其特征在于:所述的第三步中的拓展微分代数方程组差分得到的第n个仿真时步下的非线性方程组E(p)为:
E(p):
其中h为仿真步长。拓展微分代数方程组和对应的拓展状态变量定义如下:拓展微分代数方程组为:
其中,y和分别为原微分代数方程组D0(p)的状态变量及其导数,f为原微分代数方程组D0(p)的定义函数,和分别为拓展微分代数方程组的拓展状态变量及其导数,si和分别为参数向量p的第i个分量对应的轨迹灵敏度及其导数。相应的,拓展微分代数方程组{D0(p),D1(p)…DNP(p)}中的各分量定义为:
D0(p):y(t0)=y0(p)
Di(p):
其中t0为仿真时段的初始时刻,y0为状态变量的初值。
3.根据权利要求1所述的电力系统并行化轨迹灵敏度获取方法,其特征在于:所述的第四步中第n个仿真时步第j次牛顿迭代的L(p)具有下述形式:
其中和分别为第n个仿真时步第j次牛顿迭代的拓展状态变量及其导数。J为精确雅可比矩阵,具有下述矩阵结构:
J0和Ji(其中i取1…NP)为雅可比矩阵中的对角线和边界分块。
将所有矩阵块Ji忽略,从而得到不精确雅可比矩阵
于是L(p)即可解耦为NP+1个独立线性方程组L0(p)和Li(p)(其中i取1…NP),即:
L0(p):
Li(p):
4.根据权利要求1所述的电力系统并行化轨迹灵敏度获取方法,其特征在于:所述的第五步中计算处理单元是指CPU核心、计算服务器节点或其他具有完整的浮点计算与逻辑处理能力的计算设备。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410203928.XA CN104036118B (zh) | 2014-05-14 | 2014-05-14 | 一种电力系统并行化轨迹灵敏度获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410203928.XA CN104036118B (zh) | 2014-05-14 | 2014-05-14 | 一种电力系统并行化轨迹灵敏度获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104036118A true CN104036118A (zh) | 2014-09-10 |
CN104036118B CN104036118B (zh) | 2017-06-13 |
Family
ID=51466888
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410203928.XA Active CN104036118B (zh) | 2014-05-14 | 2014-05-14 | 一种电力系统并行化轨迹灵敏度获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104036118B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105098802A (zh) * | 2015-08-12 | 2015-11-25 | 中国电力科学研究院 | 基于轨迹灵敏度的电力系统切机控制决策方法 |
CN107863772A (zh) * | 2017-11-10 | 2018-03-30 | 浙江大学 | 一种基于广义Galerkin的高阶轨迹灵敏度计算方法 |
CN108629136A (zh) * | 2018-05-14 | 2018-10-09 | 北京理工大学 | 一种连续时间系统的并行仿真及误差补偿方法 |
CN108777485A (zh) * | 2018-06-08 | 2018-11-09 | 国网浙江省电力有限公司电力科学研究院 | 考虑离散控制的高阶轨迹灵敏度计算方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101409447A (zh) * | 2008-11-27 | 2009-04-15 | 浙江大学 | 基于部分自动微分技术的电力系统潮流优化方法 |
CN101478159A (zh) * | 2009-01-12 | 2009-07-08 | 浙江大学 | 基于差分等式约束转化降阶的暂态稳定约束潮流优化方法 |
CN101694940A (zh) * | 2009-10-23 | 2010-04-14 | 国网电力科学研究院 | 一种考虑暂态安全约束的最优潮流实现方法 |
CN103166219A (zh) * | 2013-01-30 | 2013-06-19 | 中国电力科学研究院 | 基于超实时仿真的电网在线稳定决策控制系统及控制方法 |
CN103605852A (zh) * | 2013-11-25 | 2014-02-26 | 国家电网公司 | 一种用于大规模电网机电暂态实时仿真的并行拓扑方法 |
-
2014
- 2014-05-14 CN CN201410203928.XA patent/CN104036118B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101409447A (zh) * | 2008-11-27 | 2009-04-15 | 浙江大学 | 基于部分自动微分技术的电力系统潮流优化方法 |
CN101478159A (zh) * | 2009-01-12 | 2009-07-08 | 浙江大学 | 基于差分等式约束转化降阶的暂态稳定约束潮流优化方法 |
CN101694940A (zh) * | 2009-10-23 | 2010-04-14 | 国网电力科学研究院 | 一种考虑暂态安全约束的最优潮流实现方法 |
CN103166219A (zh) * | 2013-01-30 | 2013-06-19 | 中国电力科学研究院 | 基于超实时仿真的电网在线稳定决策控制系统及控制方法 |
CN103605852A (zh) * | 2013-11-25 | 2014-02-26 | 国家电网公司 | 一种用于大规模电网机电暂态实时仿真的并行拓扑方法 |
Non-Patent Citations (1)
Title |
---|
吕晓祥,王建全: "电力系统暂态稳定预防控制算法综述", 《电力系统保护与控制》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105098802A (zh) * | 2015-08-12 | 2015-11-25 | 中国电力科学研究院 | 基于轨迹灵敏度的电力系统切机控制决策方法 |
CN105098802B (zh) * | 2015-08-12 | 2021-07-13 | 中国电力科学研究院 | 基于轨迹灵敏度的电力系统切机控制决策方法 |
CN107863772A (zh) * | 2017-11-10 | 2018-03-30 | 浙江大学 | 一种基于广义Galerkin的高阶轨迹灵敏度计算方法 |
CN107863772B (zh) * | 2017-11-10 | 2020-05-08 | 浙江大学 | 一种基于广义Galerkin的高阶轨迹灵敏度计算方法 |
CN108629136A (zh) * | 2018-05-14 | 2018-10-09 | 北京理工大学 | 一种连续时间系统的并行仿真及误差补偿方法 |
CN108629136B (zh) * | 2018-05-14 | 2021-07-02 | 北京理工大学 | 一种连续时间系统的并行仿真及误差补偿方法 |
CN108777485A (zh) * | 2018-06-08 | 2018-11-09 | 国网浙江省电力有限公司电力科学研究院 | 考虑离散控制的高阶轨迹灵敏度计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104036118B (zh) | 2017-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lu et al. | Moving extremum surrogate modeling strategy for dynamic reliability estimation of turbine blisk with multi-physics fields | |
CN103646152B (zh) | 一种基于矩阵指数的电力系统电磁暂态仿真方法 | |
CN104750932A (zh) | 一种混合不确定性下基于代理模型的结构可靠性分析方法 | |
CN101232180A (zh) | 一种配电系统负荷模糊建模装置及方法 | |
Huang et al. | Damage identification of bridge structures considering temperature variations-based SVM and MFO | |
CN107069696B (zh) | 一种电力系统状态估计的并行计算方法 | |
CN104123463A (zh) | 一种随机动态载荷的时域识别方法 | |
CN104036118A (zh) | 一种电力系统并行化轨迹灵敏度获取方法 | |
CN107123994A (zh) | 区间无功优化模型的线性化求解方法 | |
CN103199525A (zh) | 一种基于等效节点注入电流变换的配电网潮流计算方法 | |
CN107681685A (zh) | 一种考虑光伏非线性相关性的电力系统概率潮流计算方法 | |
CN102364490B (zh) | 基于层次分析模型的自动同调识别方法 | |
CN102968529A (zh) | 一种供水管网模型计算结果不确定性区间的量化方法 | |
Alkayem et al. | A new self-adaptive quasi-oppositional stochastic fractal search for the inverse problem of structural damage assessment | |
CN105243245A (zh) | 一种基于Petri网的电路模块故障机理相关关系的可靠性建模方法 | |
CN104967114B (zh) | 电网负荷实时数字建模方法及系统 | |
CN105846437A (zh) | 一种计及互动相关性的cholesky分解半不变量潮流计算方法 | |
CN103093079B (zh) | 基于仿射算术的辐射型配电网不确定潮流分析方法 | |
CN105353270A (zh) | 考虑分布式电源并网的电能质量扰动源容错定位方法 | |
Zhang et al. | Measurement‐based dynamic load modelling using time‐domain simulation and parallel‐evolutionary search | |
CN103106332A (zh) | 一种测量不确定度的分析方法 | |
CN106485392A (zh) | 互联发电系统可靠性评估方法 | |
CN113177717B (zh) | 一种基于影响增量灵敏度的输电系统韧性快速评估方法 | |
CN104037756B (zh) | 一种含复杂电力设备模型的电力系统稳定评估方法 | |
CN103713997A (zh) | 一种蜕变关系形式化描述与分解方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |