CN103473478B - 基于能量函数的电网暂态稳定性评估方法 - Google Patents
基于能量函数的电网暂态稳定性评估方法 Download PDFInfo
- Publication number
- CN103473478B CN103473478B CN201310459527.6A CN201310459527A CN103473478B CN 103473478 B CN103473478 B CN 103473478B CN 201310459527 A CN201310459527 A CN 201310459527A CN 103473478 B CN103473478 B CN 103473478B
- Authority
- CN
- China
- Prior art keywords
- energy
- electromotor
- transient
- disturbance
- energy function
- 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.)
- Expired - Fee Related
Links
Abstract
本发明涉及基于广域测量系统与能量函数结合的电网暂态稳定性评估方法。本发明的公开了一种基于能量函数的电网暂态稳定性评估方法。本发明的技术方案是,基于能量函数的电网暂态稳定性评估方法,以广域测量系统采集数据,建立暂态能量函数模型;采集暂态过程中系统数据,构建暂态能量函数模型,得到系统中n台发电机的能量总和;根据二次扰动过程的能量变化轨迹,构建基于二次扰动的暂态稳定性量化指标模型;构建能量裕度指标根据能量裕度指标进行Ts(t)判断。本发明借助广域测量系统的参数在线获取优势和能量函数本身的计算优势,可以提高计算速度,简化计算过程,进一步提高了直接法在电力系统的应用。
Description
技术领域
本发明涉及电力系统暂态稳定分析技术领域,特别涉及基于广域测量系统与能量函数结合的电网暂态稳定性评估方法。
背景技术
近年来,由于电力系统的规模不断扩大,电网结构日益复杂,电力系统安全稳定性问题日趋严重。而在电力系统稳定破坏事故中,暂态稳定破坏事故占据很大部分,可见暂态稳定研究的重要性。能量裕度一直被认为是评价电力系统暂态稳定性的一个重要指标,而能量裕度的在线应用是实现暂态稳定性评估的关键技术,如何更快速、更简便地实现能量裕度的在线应用一直是电力系统暂态稳定分析领域的研究热点。目前在电力系统暂态稳定性评估领域中,研究思路主要集中在两个类别:不确定性的分析方法和确定性的分析法。关于不确定性的分析方法,有研究者基于风险理论来构建暂态稳定性安全风险评估指标,将系统整体划分为几个关键部分来构建暂态安全风险指标,也有基于概率分布来构建暂态稳定性概率模型的。确定性的分析方法主要有基于支持向量机的方法、解轨迹跟踪法和能量函数法。以上分析理论方法包括离线评估、确定不稳定平衡点等,也有基于能量函数,从不同的角度来进行了暂态稳定性的机理及评估指标研究的,包括利用二次扰动法的。但是这些方法计算量大,评估过程复杂,速度低,且容易引起误差。在电力系统发展的新形势下,特别是广域测量系统的大量接入的发展趋势下,基于能量函数的暂态稳定性的在线评估要求向现有的理论方法提出新的挑战,更快速、更简便的评估方法是电力运行人员关注的重点。
发明内容
本发明的所要解决的技术问题,就是克服上述现有技术的不足,提供一种更快速、更简便的基于能量函数的电网暂态稳定性评估方法。
本发明解决所述技术问题,采用的技术方案是,基于能量函数的电网暂态稳定性评估方法,以广域测量系统采集数据,建立暂态能量函数模型,其特征在于,包括以下步骤:
a、采集暂态过程中系统数据,构建暂态能量函数模型,得到系统中n台发电机的能量总和Etot:
其中,为第i台发电机的动能;为第i台发电机的势能;δi、ωi分别为第i台发电机的转子角和转速;为第i台发电机初始稳定时刻的功角;Mi为第i台发电机的惯性时间常数;fi(δi)=Pmi-Pei,Pmi、Pei分别为第i台发电机的机械功率和电磁功率;n为正整数,i≤n;
b、根据二次扰动过程的能量变化轨迹,构建基于二次扰动的暂态稳定性量化指标模型;设时段0-t_start,是一次扰动持续时段,t_start是第一次扰动动能降到最低点时刻,对应系统动能最小点,此时对系统第二次施加一个足够大的扰动,令系统不稳定,系统的动能增加,动能轨迹经过第二次动能最大时刻为t_clear,对应的动能为Eke(t_clear),此时切除故障,待到系统完全恢复稳定,对应时刻为ts2,对应的动能为对应的势能为得到势能极值EPEBS:
c、构建能量裕度指标Ts(t):
d、根据能量裕度指标Ts(t)判断:当Ts(t)∈(0,1],系统是强健的,还具有经受扰动或故障的能力;当Ts(t)=0时,系统已经濒临失稳的临界状态;当Ts(t)<0,系统已失稳。
进一步的,采用相位测量单元采集系统数据,在线获取参数Pmi、Pei、δi、ωi。
进一步的,系统处于稳定平衡状态时,能量函数对时间的导数为零,即:
为第i台发电机的转速对时间的导数。
进一步的,当阻尼存在时
式中,Di为阻尼常数。
本发明的有益效果是,借助广域测量系统的参数在线获取优势和能量函数本身的计算优势,可以提高计算速度,简化计算过程,进一步提高了直接法在电力系统的应用。在计算过程中能够避免直接求取稳定域边界上的不稳定平衡点,还能够适应不同的系统运行条件,为电力系统暂态分析理论和方法的实际应用提供了支撑。
附图说明
图1广域测量系统在三机系统中的应用结构图;
图2三机系统接线图;
图3三机系统的稳定轨迹及初始稳定平衡点;
图4三机系统一次扰动后的运动轨迹及稳定平衡点;
图5三机系统在二次扰动后的运动轨迹及稳定平衡点;
图6两次扰动过程中的系统运动轨迹;
图7三机系统从启动到经历两次扰动的动能变化曲线;
图8三机系统对应不同扰动的势能变化;
图9三机系统从启动到经历两次扰动的能量裕度变化曲线;
具体实施方式
下面结合附图及实施例,详细描述本发明的技术方案。
本发明结合广域测量系统与暂态能量函数模型,分析了基于二次扰动法的势能边界值计算法,并推导了系统暂态稳定性量化指标,其实现方法分为两部分阐述如下:
第一部分:广域测量系统的应用
本发明中,广域测量系统的应用,其具体实现框架包括了9个模块:故障中导纳矩阵模块、故障后导纳矩阵模块、故障中系统参数模块、切除故障时系统参数模块、故障切除后系统稳定时的参数模块、初始稳定系统参数模块、暂态能量计算模块、势能边界值计算模块、暂态稳定性指标计算模块。以上所有模块计算所需参数均由广域测量系统中的相位测量单元(Phase Measurement Unit,PMU)采集传输。暂态能量模块的计算需要故障中的导纳矩阵参数以及故障中的发电机角速度ωf,功角δf和电磁功率Pef。势能边界值的计算需要获取故障切除时刻的发电机角速度ωc,临界切除角δc和切除时刻的电磁功率Pec,故障后系统再次稳定时刻的功角δfs,角速度ωfs,电磁功率Pefs以及初始稳定时刻的功角δs,角速度ωs,电磁功率Pes。暂态能量计算模块和势能边界值计算模块将计算结果提供给暂态稳定性指标计算模块,用于稳定性评价指标的计算。
第二部分:基于能量函数的电网暂态稳定性评估方法
步骤(1):采集暂态过程中系统数据,构建暂态能量函数模型,得到系统中n台发电机的能量总和Etot:
其中,为第i台发电机的动能;为第i台发电机的势能;δi、ωi分别为第i台发电机的转子角和转速;为第i台发电机初始稳定时刻的功角;Mi为第i台发电机的惯性时间常数;fi(δi)=Pmi-Pei,Pmi、Pei分别为第i台发电机的机械功率和电磁功率;n为正整数,i≤n。
该步骤包括:
步骤(1.1):根据单机无穷大系统建立第i台发电机转子运动方程:
式中,δi、ωi分别为发电机的转子角和转速,分别为δi对时间的一阶和二阶导数;Pmi、Pei分别为发电机的机械功率和电磁功率;Mi为惯性时间常数。令
fi(δi)=Pmi-Pei
步骤(1.2):发电机的能量函数表达式:
根据李亚普罗夫直接法,结合步骤(1.1)可得:
其中,
步骤(1.3):多机系统中,发电机的能量函数表达式:
步骤(1.4):故障后稳定平稳时能量函数对时间的导数:
在实际系统中,阻尼是存在的,即
式中,Di为阻尼常数。
步骤(2):根据二次扰动过程的能量变化轨迹,构建基于二次扰动的暂态稳定性量化指标模型;设时段0-t_start,是一次扰动持续时段,t_start是第一次扰动动能降到最低点时刻,对应系统动能最小点,此时对系统第二次施加一个足够大的扰动,令系统不稳定,系统的动能增加,动能轨迹经过第二次动能最大时刻为t_clear,对应的动能为Eke(t_clear),此时切除故障,待到系统完全恢复稳定,对应时刻为ts2,对应的动能为对应的势能为得到势能极值EPEBS:
步骤(3):构建构建能量裕度指标Ts(t):
其中,Etot(t)是系统在在暂态过程中某时刻的总能量,该值越小,说明系统越稳定。
步骤(4)根据能量裕度指标Ts(t)判断:当Ts(t)∈(0,1],系统是强健的,还具有经受扰动或故障的能力;当Ts(t)=0时,系统已经濒临失稳的临界状态;当Ts(t)<0,系统已失稳。
实施例
研究发现,能量函数具有易构建量化指标的优势,能够通过能量函数的优势将多机系统以能量的形式完成量化,给出具体的量化值,但是模型中的部分参数难以快速准确获取,因此结合广域测量系统的参数快速准确获取优势来弥补这一缺陷。基于二次扰动法的势能边界值计算法,能够避免直接求取稳定域边界上的不稳定平衡点。下面以三机系统(n=3)来详细阐述本发明的技术方案。
第一部分:广域测量系统应用框架
本发明中,广域测量系统的应用,其具体实现框架如图1所示。
在图1中,包括了9个模块:故障中导纳矩阵模块、故障后导纳矩阵模块、故障中系统参数模块、切除故障时系统参数模块、故障切除后系统稳定时的参数模块、初始稳定系统参数模块、暂态能量计算模块、势能边界值计算模块、暂态稳定性指标计算模块。以上所有模块计算所需参数均由PMU采集传输。暂态能量模块的计算需要故障中的导纳矩阵参数以及故障中的发电机角速度ωf,功角δf和电磁功率Pef。势能边界值的计算需要获取故障切除时刻的发电机角速度ωc,临界切除角δc和切除时刻的电磁功率Pec,故障后系统再次稳定时刻的功角δfs,角速度ωfs,电磁功率Pefs以及初始稳定时刻的功角δs,角速度ωs,电磁功率Pes。暂态能量计算模块和势能边界值计算模块将计算结果提供给暂态稳定性指标计算模块,用于稳定性评价指标的计算。
第二部分:基于能量函数的电网暂态稳定性评估方法,该部分分为四个步骤展开
第一步:构建三机系统的能量函数模型;
图2为三机系统接线图。其中,G1,G2和G3为发电机,G3为参考发电机,在每台发电机母线处都带有不同大小的负荷,R12、R23和R13为线路电阻,R12=R23=R13,x12、x23和x13为线路电抗,且x12=x23=x13。Load1、Load2、Load3分别为3台发电机的负载。机端电压分别为U1∠δ1,U2∠δ2和U3∠δ3,U1、U2、U3分别为节点1、2、3的电压幅值,δ1、δ2、δ3为节点1、2、3的相角。
三机系统的能量函数模型展开为:
其中,i和j分别为节点编号;Cij=UiUjBij,为支路ij的电导;Dij=UiUjGij,为支路ij的电纳。
三机系统最终的能量函数模型表达式为:
其中,动能部分为:
势能部分为:
第二步:实施二次扰动并利用插值法求取EPEBS。
首先对如图2所示的三机系统实施一次扰动,根据式(3)仿真动能变化,记录系统动能最小Eke(t_start)=0.210以及对应该值的时刻时段t_start=59.133秒,在t_start+=59.133秒时刻,对系统第二次施加一个足够大的扰动,令系统不稳定,系统的动能增加,记录动能轨迹最大动能值Eke(t_clear)=7.976和对应的时刻t_clear=66.156秒,此时切除故障;待到系统第二次受扰切除故障后完全恢复稳定,对应时刻为ts2=78.231秒,对应的动能为根据式(4)可得对应的势能为能量裕度计算如下:
第三步:根据所求取的EPEBS获取能量裕度指标Ts(t);
其表达式为:
第四步,根据能量裕度指标Ts(t)对三机系统进行稳定性分析:
(一)初始稳定平衡点计算。
初始稳定平衡点是指系统受扰前的稳定平衡点,直接通过运行三机系统来获取,其结果如图3所示。
图3展示了从启动到稳定的过程(0s-40s)中系统轨迹变化曲线,从图3可以直接得到三机系统的稳定轨迹及初始稳定平衡点。图3中的三个点Sep1,Sep2和Sep3分别是发电机G1、G2和G3的初始稳定平衡点,分别对应坐标(0.99,54.37°)、(0.99,55.78°)和(0.99,42.13°)。从图中可以看出,三台发电机在启动以后都以很快的速度到达了稳定平衡点。
(二)一次扰动后的运动轨迹及稳定平衡点计算。
一次扰动后的运动轨迹及稳定平衡点如图4所示。图4展示了当一次扰动发生在50s时刻,从故障发生到系统达到稳定平衡点(50s-60s)的轨迹。图4中的Sep1,Sep2和Sep3分别是发电机G1、G2和G3经历一次扰动后的稳定平衡点,分别对应坐标(0.978,54.67°)、(0.978,55.90°)和(0.978,42.15°)。与图3中的稳定平衡点相比,无论是平衡点的角速度还是转子角,都有极小的移动,说明三相短路故障对该系统的整体暂态稳定性有较大影响,导致发电机不能回到初始稳定点。
(三)二次扰动后的运动轨迹及稳定平衡点计算。
图5展示了三机系统在二次扰动后的运动轨迹及各台发电机的稳定平衡点(ω,δ)。在图5的仿真过程中,二次扰动(三相短路故障,发生在支路1-2的中点)发生在70s时刻,SEP1,SEP2和SEP3分别是发电机G1、G2和G3经历二次扰动后的稳定平衡点,分别对应坐标(0.978,54.66°)、(0.978,55.88°)和(0.978,42.14°)。相比较于图4,三台发电机的平衡点位置几乎没有移动,单纯从平衡点的变化来看,第二次三相短路故障没有对系统的稳定性造成太大的影响。
(四)故障过程中运动轨迹及动能和势能的变化曲线仿真。
图6、7、8是系统在故障过程中运动轨迹及动能和势能的变化。图6表示了两次扰动过程中的系统运动轨迹。根据两次扰动后的运动轨迹、动能和势能计算模型,可以得到两次扰动过程中系统的动能和势能变化曲线如图7和图8所示。在系统中,系统任何时刻动能和势能是守恒的,从图7和图8所示的动能和势能变化曲线也验证了这一动能和势能的守恒性。另外,对于系统而言,两次扰动相当于是对系统注入了两次动能,这两次注入的动能在忽略系统其它动能消耗的情况下最终将全部转化为不同形式的势能存在于系统中。因此,在图8中,第一次扰动注入的动能全部转化为势能后,在达到稳定状态后不再变化,在第二次扰动再次注入动能后,势能阶梯型爬升,势能随着时间的变化呈现迭加的态势,一旦动能全部转化为势能以后,系统的势能将维持不变。
(五)三机系统从启动到经历两次扰动的能量裕度变化曲线仿真。
图9是三机系统从启动到经历两次扰动的能量裕度变化曲线。图9中展示了系统在初态以及经受了两次扰动后的能量裕度的变化情况。系统初始状态时的能量裕度指标Ts(t)在1附近,Ts(t)∈(0,1],系统是强健的,还具有经受扰动或故障的能力,而经历了一次扰动以后,能量裕度指标Ts(t)下降到0.5左右,Ts(t)∈(0,1],系统依然是稳定的。在经历了二次扰动后,Ts(t)则下降到0.1左右,此时系统的能量裕度非常小,系统已经开始脆弱化,如果再有一次扰动,系统将濒临失稳。
Claims (4)
1.基于能量函数的电网暂态稳定性评估方法,以广域测量系统采集数据,建立暂态能量函数模型,其特征在于,包括以下步骤:
a、采集暂态过程中系统数据,构建暂态能量函数模型,得到系统中n台发电机的能量总和Etot:
其中,为第i台发电机的动能; 为第i台发电机的势能;δi、ωi分别为第i台发电机的转子角和转速; 为第i台发电机初始稳定时刻的功角;Mi为第i台发电机的惯性时间常数;fi(δi)=Pmi-Pei,Pmi、Pei分别为第i台发电机的机械功率和电磁功率;n为正整数,i≤n;
b、根据二次扰动过程的能量变化轨迹,构建基于二次扰动的暂态稳定性量化指标模型;设时段0-t_start,是一次扰动持续时段,t_start是第一次扰动动能降到最低点时刻,对应系统动能最小点,此时对系统第二次施加一个足够大的扰动,令系统不稳定,系统的动能增加,动能轨迹经过第二次动能最大时刻为t_clear,对应的动能为此时切除故障,待到系统完全恢复稳定,对应时刻为ts2,对应的动能为对应的势能为得到势能极值EPEBS:
c、构建能量裕度指标Ts(t):
d、根据能量裕度指标Ts(t)判断:当Ts(t)∈(0,1],系统是强健的,还具有经受扰动或故障的能力;当Ts(t)=0时,系统已经濒临失稳的临界状态;当Ts(t)<0,系统已失稳。
2.根据权利要求1所述的基于能量函数的电网暂态稳定性评估方法,其特征在于,采用相位测量单元采集系统数据,在线获取参数Pmi、Pei、δi、ωi。
3.根据权利要求1所述的基于能量函数的电网暂态稳定性评估方法,其特征在于,系统处于稳定平衡状态时,能量函数对时间的导数为零,即:
为第i台发电机的转速对时间的导数。
4.根据权利要求1所述的基于能量函数的电网暂态稳定性评估方法,其特征在于,当阻尼存在时
式中,Di为阻尼常数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310459527.6A CN103473478B (zh) | 2013-09-30 | 2013-09-30 | 基于能量函数的电网暂态稳定性评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310459527.6A CN103473478B (zh) | 2013-09-30 | 2013-09-30 | 基于能量函数的电网暂态稳定性评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103473478A CN103473478A (zh) | 2013-12-25 |
CN103473478B true CN103473478B (zh) | 2017-03-01 |
Family
ID=49798325
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310459527.6A Expired - Fee Related CN103473478B (zh) | 2013-09-30 | 2013-09-30 | 基于能量函数的电网暂态稳定性评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103473478B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103887800B (zh) * | 2014-03-19 | 2015-09-30 | 电子科技大学 | 基于能量裕度约束的电压在线协调控制方法 |
CN104680323B (zh) * | 2015-03-02 | 2018-05-25 | 华北电力大学 | 电力系统能量稳定域构建系统及方法 |
CN104732083A (zh) * | 2015-03-23 | 2015-06-24 | 东北电力大学 | 基于支路势能的负荷暂态能量函数构建与评价方法 |
CN108268670B (zh) * | 2016-12-30 | 2023-07-21 | 中国电力科学研究院 | 一种不依赖系统扰动轨迹的电网稳定态势评估方法及系统 |
CN107370180B (zh) * | 2017-07-27 | 2020-11-13 | 华北电力大学 | 一种定量分析风电接入对系统暂态稳定性影响的方法 |
CN110782107A (zh) * | 2018-07-31 | 2020-02-11 | 中国电力科学研究院有限公司 | 一种评估电网暂态稳定性的方法及系统 |
CN109217287B (zh) * | 2018-08-06 | 2020-08-18 | 华南理工大学 | 交直流系统暂态电压稳定安全域的求解方法 |
CN109799456B (zh) * | 2019-02-28 | 2022-12-13 | 中国电力科学研究院有限公司 | 一种确定单机无穷大系统暂态功角稳定性的方法及系统 |
CN110445174B (zh) * | 2019-08-12 | 2022-04-19 | 电子科技大学 | 一种考虑随机风速影响的风电场系统暂态稳定性评估方法 |
CN112819647B (zh) * | 2020-12-08 | 2023-04-07 | 广东电网有限责任公司 | 一种电网故障后暂态稳定性实时预测方法及装置 |
CN112952800B (zh) * | 2021-02-02 | 2022-07-26 | 中国南方电网有限责任公司 | 用于电网暂态稳定水平的分段评估指标获取方法及装置 |
CN114374202B (zh) * | 2021-12-17 | 2023-04-07 | 清华大学 | 平衡点无关的电力系统暂态同步稳定性分析方法及系统 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100403618C (zh) * | 2006-07-17 | 2008-07-16 | 天津大学 | 电力系统暂态稳定分析与控制方法 |
CN102185325A (zh) * | 2011-04-25 | 2011-09-14 | 东北电力大学 | 基于网络量测信息的电力系统暂态稳定性定量评价方法 |
-
2013
- 2013-09-30 CN CN201310459527.6A patent/CN103473478B/zh not_active Expired - Fee Related
Non-Patent Citations (1)
Title |
---|
"网格计算模式下的PSS/E暂态稳定仿真分析";潘爱强;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20080615;第2008年卷(第6期);第C042-93页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103473478A (zh) | 2013-12-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103473478B (zh) | 基于能量函数的电网暂态稳定性评估方法 | |
CN103746370B (zh) | 一种风电场可靠性建模方法 | |
CN102136036B (zh) | 适于电力系统小扰动稳定分析的双馈风电场等值建模方法 | |
CN102799950B (zh) | 基于粒子群算法的船舶电网重构优化方法 | |
CN103973203B (zh) | 一种适用于安全稳定分析的大型光伏电站在线等值建模方法 | |
CN103700036B (zh) | 一种适于电力系统多时间尺度的暂态稳定性投影积分方法 | |
CN103810646B (zh) | 一种基于改进投影积分算法的有源配电系统动态仿真方法 | |
CN102760191B (zh) | 基于转速分群的双馈机组风电场等值建模系统及方法 | |
CN104242300B (zh) | 考虑不同电源形式的配电网中分布式电源的选址定容方法 | |
CN104269867B (zh) | 一种节点扰动功率转移分布均衡度分析方法 | |
CN103488610B (zh) | 一种基于稀疏存储的非零元遍历的电网网络方程求解方法 | |
CN103474992A (zh) | 电力系统节点电压稳定状态的实时在线辨识判据 | |
CN109977541A (zh) | 一种电力信息物理系统联合仿真平台及其同步方法 | |
CN103777525B (zh) | 风电场仿真机与rtds仿真器的自定义接口 | |
CN101794998A (zh) | 基于多机电力系统单一发电机电磁功率简洁表达方式的暂态稳定在线分析方法 | |
CN103995948B (zh) | 一种基于多项式模型的振荡中心电压预测方法 | |
CN102420427A (zh) | 一种考虑外网等值的区域电网电压稳定裕度计算方法 | |
CN103530473A (zh) | 一种含大规模光伏电站的电力系统随机生产模拟方法 | |
CN106169115A (zh) | 一种基于解析‑数值结合的两区域电网最低频率估计算法 | |
CN204314716U (zh) | 一种用于风电场控制系统硬件在环测试平台 | |
CN106229976A (zh) | 基于数据驱动的暂态功角稳定态势预估方法 | |
CN105281371A (zh) | 一种考虑风力发电的可伸缩有功静态安全域 | |
CN105787604A (zh) | 一种输配协同的配网态势快速感知方法 | |
CN105447658A (zh) | 一种含风电随机模糊注入电力系统波动的电压崩溃点求取方法 | |
CN103077291B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170301 Termination date: 20170930 |
|
CF01 | Termination of patent right due to non-payment of annual fee |