CN105470931B - 一种不受母线数据不同步影响的母线差动保护方法 - Google Patents

一种不受母线数据不同步影响的母线差动保护方法 Download PDF

Info

Publication number
CN105470931B
CN105470931B CN201510999678.XA CN201510999678A CN105470931B CN 105470931 B CN105470931 B CN 105470931B CN 201510999678 A CN201510999678 A CN 201510999678A CN 105470931 B CN105470931 B CN 105470931B
Authority
CN
China
Prior art keywords
bus
current
fault
voltage
differential protection
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
CN201510999678.XA
Other languages
English (en)
Other versions
CN105470931A (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
Xian Jiaotong University
State Grid Shaanxi Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Xian Jiaotong University
State Grid Shaanxi 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, Xian Jiaotong University, State Grid Shaanxi Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201510999678.XA priority Critical patent/CN105470931B/zh
Publication of CN105470931A publication Critical patent/CN105470931A/zh
Application granted granted Critical
Publication of CN105470931B publication Critical patent/CN105470931B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02HEMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
    • H02H7/00Emergency protective circuit arrangements specially adapted for specific types of electric machines or apparatus or for sectionalised protection of cable or line systems, and effecting automatic switching in the event of an undesired change from normal working conditions
    • H02H7/22Emergency protective circuit arrangements specially adapted for specific types of electric machines or apparatus or for sectionalised protection of cable or line systems, and effecting automatic switching in the event of an undesired change from normal working conditions for distribution gear, e.g. bus-bar systems; for switching devices
    • H02H7/226Emergency protective circuit arrangements specially adapted for specific types of electric machines or apparatus or for sectionalised protection of cable or line systems, and effecting automatic switching in the event of an undesired change from normal working conditions for distribution gear, e.g. bus-bar systems; for switching devices for wires or cables, e.g. heating wires
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02HEMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
    • H02H7/00Emergency protective circuit arrangements specially adapted for specific types of electric machines or apparatus or for sectionalised protection of cable or line systems, and effecting automatic switching in the event of an undesired change from normal working conditions
    • H02H7/26Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured
    • H02H7/261Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured involving signal transmission between at least two stations

Landscapes

  • Emergency Protection Circuit Devices (AREA)

Abstract

本发明公开一种不受母线数据不同步影响的母线差动保护方法,包括:步骤一、采集与母线相连的每一回出线保护安装处A、B、C三相电压和电流数据;步骤二、对采集得到的A、B、C三相电压和电流进行低通滤波、采样保持和A/D转换后,得到A、B、C三相电压和电流采样值;步骤三、对采集到的电压电流瞬时值进行故障分量提取;步骤四、对得到的电压电流故障分量瞬时值做全周傅氏变换转换成对应的相量值;步骤五、计算结果判断是否发生母线故障。通过本发明方法可以把从各回出线采集到的原本不同步的电流量转化为同一时刻,从而克服了数据传输不同步对母线电流差动保护的影响。

Description

一种不受母线数据不同步影响的母线差动保护方法
技术领域
本发明涉及电力系统继电保护领域,特别涉及一种母线差动保护方法。
背景技术
目前智能变电站在我国逐步推广使用,在220kV以上变电站规程规定母线保护双重化配置,但110kV母线保护只要求单重配置,直采直跳。对于部分地区重要的110kV变电站,上述规程不能完全满足可靠性要求。为此需配置网络式冗余式的母差保护,网采网跳,实现准双重化。如何消除数据不同步影响是解决冗余式母线保护的关键问题。
此外,分布式母线保护的发展也遇到类似问题。常规微机式母线保护,采用集中式处理方式,不存在数据同步问题。但其存在接线复杂,二次电缆引线长,互感器负担重的问题,近年来分布式母线保护成为研究的重点,但分布式母线保护的一个主要问题是如何解决采样不同步引起的保护误动问题。
发明内容
本发明的目的在于提供一种不受母线数据不同步影响的母线差动保护方法,以解决智能变电站内数据传输不同步及分布式母线保护采样不同步这一问题。
为实现上述目的,本发明采用如下技术方案:
一种不受母线数据不同步影响的母线差动保护方法,具体包括以下步骤:
步骤一、采集与母线相连的每一回出线保护安装处A、B、C三相电压和电流数据;
步骤二、对采集得到的A、B、C三相电压和电流进行低通滤波、采样保持和A/D转换后,得到A、B、C三相电压和电流采样值;将采样结果记为uAi(k),uBi(k),uCi(k),iAi(k),iBi(k),iCi(k);其中k为采样点序号,i为与母线相连线路序号;
步骤三、对采集到的电压电流瞬时值uAi(k),uBi(k),uCi(k),iAi(k),iBi(k),iCi(k)进行故障分量提取;
步骤四、对得到的电压电流故障分量瞬时值做全周傅氏变换转换成对应的相量值;
步骤五、按照式(5)计算结果判断是否发生母线故障:
则故障发生在区内;
则故障发生在线区外;
式中,K为比例制动系数;n为与母线相连出线数量;取与母线相连的任意一回线路互感器采集到的电压;同一回线路上采集到的电压和电流信息是同步的,每一回线路的为从母线看向该回线路的等效导纳。
进一步的,步骤二中利用截止频率为300Hz的低通滤波器对采集得到的A、B、C三相电压和电流进行低通滤波;步骤二中采样频率为4kHz,每工频周波采样点数N=80,采样时间间隔Ts=0.25ms。
进一步的,步骤三种进行故障分量提取具体为:
式中:表示A、B、C三相中的任意一相,分别为电压电流的故障分量,u[0],i[0]为故障前电压电流瞬时值。
进一步的,步骤四的具体变换方式为;
进一步的,步骤五中:0<K<1。
进一步的,K的取值为0.5。
进一步的,步骤一中利用间隔层合并单元采集或分布式保护单元采集的每一回出线保护安装处A、B、C三相电压和电流数据。
相对于现有技术,本发明具有以下有益效果:针对智能变电站中过程层合并单元向站控层传输数据时由于延时不同,造成站控层接收到的不同间隔数据不同步,或分布式母线差动保护采样不同步这两种情况,本发明提出了一种不受母线数据不同步影响的差动保护方法,解决了这种不同步对母线差动保护的影响。
附图说明
图1为母线接线示意图;
图2为母线故障对应故障序网图;
图3为线路故障对应故障序网图;
图4为某地变电站110kV母线接线示意图;
图5为区内故障保护仿真计算结果示意图;
图6为区外故障保护仿真计算结果示意图。
具体实施方式
下面结合附图,对本发明应用于母线差动保护过程做进一步详细描述。
本发明一种不受母线数据不同步影响的母线差动保护方法,包括以下步骤:
步骤一、利用电子型传感器采集与母线相连的每一回出线保护安装处A、B、C三相电压和电流数据。
步骤二、利用截止频率为300Hz的低通滤波器对采集得到的A、B、C三相电压和电流进行低通滤波,通过采样保持和A/D转换后,得到A、B、C三相电压和电流采样值。将A、B、C三相电压和电流采样值采样结果分别记为:uAi(k),uBi(k),uCi(k),iAi(k),iBi(k),iCi(k);其中k为采样点序号,i为与母线相连线路序号。设采样频率为4kHz,此时每工频周波采样点数N=80,采样时间间隔Ts=0.25ms。
步骤三、对采集到的电压电流瞬时值uAi(k),uBi(k),uCi(k),iAi(k),iBi(k),iCi(k)进行故障分量提取,具体计算方法如下:
式中:表示A、B、C三相中的任意一相,分别为电压电流的故障分量,u[0],i[0]为故障前电压电流瞬时值。
步骤四、对得到的电压电流故障分量瞬时值做全周傅氏变换转换成对应的相量值。具体变换方式如下:
式中:N为一个工频周期内采样点个数,根据IEC61850规约规定,智能变电站中合并单元采样频率为4kHz,即N取80。
步骤五、按照以下判据判断进行母线差动保护整定:
式中,K为比例制动系数,取0.5;以图1为例,与母线相连出线共5回,即n取5;取与母线相连的任意一回线路互感器采集到的电压。同一回线路上采集到的电压和电流信息是同步的,所以每一回线路的即为从母线看向该回线路的等效导纳。这个等效导纳大小与系统结构有关,在故障发生后几乎不随时间变化而变化。任意一回线路互感器上的电压理论上都是母线电压,但是在数据传输过程中不同的延时可能会造成互不相同。用任意一回线路电压和这个导纳相乘,其结果是一个等效的电流量,这个电流在时间刻上与相同。用这种方法就算出的电流可以认为是所在时刻各回出线的电流,通过这种方法可以把从各回出线采集到的原本不同步的电流量转化为同一时刻,从而克服了数据传输不同步对母线电流差动保护的影响。
分析故障分别出现在区内和区外时,式(5)计算结果的区别。
当出现区内故障时,图1对应的故障序网图如图2所示。规定由母线指向线路的方向为正方向。母线出现故障时,各条出线均向母线注入故障电流,电流方向均同相且为负,因此有此时如果比率制动系数取0<K<1,即式(5)成立。
当出现区外故障时,图1对应的故障序网图如图3所示。规定由母线指向线路的方向为正方向。线路故障对母线来说是区外故障,母线差动保护不应动作。如图3所示,假设在线路5上出现故障,则线路5电流为故障穿越电流,方向为正,其余各回出现电流为母线的穿越性电流,方向为负,且满足等式而此时远大于零。如果比率制动系数取0<K<1,式(5)不成立。
按照式(5)计算结果判断是否发生母线故障。
则说明故障发生在区内上;
则说明故障发生在区外上。
以上即为本发明应用于母线差动保护的实施方法。
利用PSCAD仿真软件建立某地330kV降压到110kV的变电站,如图4所示,将本发明提出整定方案用于该站中110kV母线,仿真采样频率设为4kHz,全周波傅氏算法数据窗长取为20ms。分别针对区内故障和区外故障进行仿真验证。根据IEC61850规定,站内通信延时上限为3ms,本仿真取极端情况,将变电站内延时设置为5ms,验证式(5)在数据不同步时保护计算的准确性,为了增加仿真结果的说服力,将常规电流差动保护计算结果绘制于同一图中作为对比。区内故障取图4中k处故障,区外故障取图4中k1故处障,仿真结果以A相接地故障为例进行说明,其他故障类型仿真结果与A相接地故障类似,故不再赘述。
当差动电流大于制动电流时差动保护动作,反之差动保护不动作。由图4可知,当发生区内故障时,利用本发明提出的差动保护方法不但可以正确选择故障位置,还能正确选择出故障相;但是常规电流差动保护在数据不同步情况下会发生健全相误动情况,即保护会出现故障时误选相。
由图5可知,当发生区外故障时,本方法方法只有很小的差动电流,差动电流值远小于制动电流,故保护不会在区外故障时误动,且有很高的可靠性;但是常规电流差动保护在数据不同步时有区外故障误动现象,不满足了继电保护选择性要求。
由上述分析可知,本发明提出的差动保护方法(即式(5))在数据不同步时比常规差动保护有更高的可靠性。为进一步验证本发明提出算法的有效性,进一步进行了不同故障过渡电阻的仿真。均已A相故障为例进行说明,由上文所述可知,对于A相故障,非故障相误动为数据不同步引起的主要问题,针对这一问题用仿真数据进行说明。故将非故障相(B相)仿真数据列于表1和表2中(电流单位:kA):
表1 区外故障时B相电流仿真结果
表2 区内故障时B相电流仿真结果
由表1和表2中数据可知,本发明方法在出现故障且数据不同步时时,保护判断结果不会收到过渡电阻不同的影响。无论故障发生在区内还是区外,本发明方法对非故障相计算的差动电流结果都远小于制动电流,有较高的可靠性。

Claims (6)

1.一种不受母线数据不同步影响的母线差动保护方法,其特征在于,包括以下步骤:
步骤一、采集与母线相连的每一回出线保护安装处A、B、C三相电压和电流数据;
步骤二、对采集得到的A、B、C三相电压和电流进行低通滤波、采样保持和A/D转换后,得到A、B、C三相电压和电流采样值;将采样结果记为uAi(k),uBi(k),uCi(k),iAi(k),iBi(k),iCi(k);其中k为采样点序号,i为与母线相连线路序号;
步骤三、对采集到的电压电流瞬时值uAi(k),uBi(k),uCi(k),iAi(k),iBi(k),iCi(k)进行故障分量提取;
步骤四、对得到的电压电流故障分量瞬时值 做全周傅氏变换转换成对应的相量值;
步骤五、按照式(5)计算结果判断是否发生母线故障:
<mrow> <mo>|</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <mi>&amp;Delta;</mi> <msub> <mover> <mi>U</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>m</mi> </msub> <mfrac> <mrow> <mi>&amp;Delta;</mi> <msub> <mover> <mi>I</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>i</mi> </msub> </mrow> <mrow> <mi>&amp;Delta;</mi> <msub> <mover> <mi>U</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>i</mi> </msub> </mrow> </mfrac> <mo>|</mo> <mo>&gt;</mo> <mi>K</mi> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <mo>|</mo> <mi>&amp;Delta;</mi> <msub> <mover> <mi>U</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>m</mi> </msub> <mfrac> <mrow> <mi>&amp;Delta;</mi> <msub> <mover> <mi>I</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>i</mi> </msub> </mrow> <mrow> <mi>&amp;Delta;</mi> <msub> <mover> <mi>U</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>i</mi> </msub> </mrow> </mfrac> <mo>|</mo> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>n</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
则故障发生在区内;
则故障发生在线区外;
式中,K为比例制动系数;n为与母线相连出线数量;取与母线相连的任意一回线路互感器采集到的电压;同一回线路上采集到的电压和电流信息是同步的,每一回线路的为从母线看向该回线路的等效导纳。
2.根据权利要求1所述的一种不受母线数据不同步影响的母线差动保护方法,其特征在于,步骤二中利用截止频率为300Hz的低通滤波器对采集得到的A、B、C三相电压和电流进行低通滤波;步骤二中采样频率为4kHz,每工频周波采样点数N=80,采样时间间隔Ts=0.25ms。
3.根据权利要求1所述的一种不受母线数据不同步影响的母线差动保护方法,其特征在于,步骤三中进行故障分量提取具体为:
式中:表示A、B、C三相中的任意一相,分别为电压电流的故障分量,u[0],i[0]为故障前电压电流瞬时值。
4.根据权利要求3所述的一种不受母线数据不同步影响的母线差动保护方法,其特征在于,步骤四的具体变换方式为;
5.根据权利要求1所述的一种不受母线数据不同步影响的母线差动保护方法,其特征在于,步骤五中:0<K<1。
6.根据权利要求1所述的一种不受母线数据不同步影响的母线差动保护方法,其特征在于,步骤五中,比例制动系数K的取值为0.5。
CN201510999678.XA 2015-12-25 2015-12-25 一种不受母线数据不同步影响的母线差动保护方法 Active CN105470931B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510999678.XA CN105470931B (zh) 2015-12-25 2015-12-25 一种不受母线数据不同步影响的母线差动保护方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510999678.XA CN105470931B (zh) 2015-12-25 2015-12-25 一种不受母线数据不同步影响的母线差动保护方法

Publications (2)

Publication Number Publication Date
CN105470931A CN105470931A (zh) 2016-04-06
CN105470931B true CN105470931B (zh) 2017-12-26

Family

ID=55608376

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510999678.XA Active CN105470931B (zh) 2015-12-25 2015-12-25 一种不受母线数据不同步影响的母线差动保护方法

Country Status (1)

Country Link
CN (1) CN105470931B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110661237B (zh) * 2019-07-03 2021-10-29 天津大学 适用于配电网自同步条件下的电流差动保护方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101800441A (zh) * 2010-03-26 2010-08-11 河南省电力公司许昌供电公司 数字化变电站集中式保护装置边界信息缺失的处理方法
CN203135427U (zh) * 2013-03-30 2013-08-14 国家电网公司 变电站母差保护系统
CN103389422A (zh) * 2013-07-19 2013-11-13 国家电网公司 线路重合闸期间非故障相母差保护误动原因的检测方法
CN103872667A (zh) * 2014-03-26 2014-06-18 国家电网公司 一种防合并单元异常大数的线路电流差动保护方法
CN104393579A (zh) * 2014-11-26 2015-03-04 国家电网公司 一种克服汲出电流对母线差动保护影响的方法
CN104638618A (zh) * 2015-02-11 2015-05-20 国网山东省电力公司潍坊供电公司 一种gis母线设备跳闸快速恢复供电方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5645535B2 (ja) * 2010-08-06 2014-12-24 株式会社日立製作所 多端子送電系統の保護継電システム

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101800441A (zh) * 2010-03-26 2010-08-11 河南省电力公司许昌供电公司 数字化变电站集中式保护装置边界信息缺失的处理方法
CN203135427U (zh) * 2013-03-30 2013-08-14 国家电网公司 变电站母差保护系统
CN103389422A (zh) * 2013-07-19 2013-11-13 国家电网公司 线路重合闸期间非故障相母差保护误动原因的检测方法
CN103872667A (zh) * 2014-03-26 2014-06-18 国家电网公司 一种防合并单元异常大数的线路电流差动保护方法
CN104393579A (zh) * 2014-11-26 2015-03-04 国家电网公司 一种克服汲出电流对母线差动保护影响的方法
CN104638618A (zh) * 2015-02-11 2015-05-20 国网山东省电力公司潍坊供电公司 一种gis母线设备跳闸快速恢复供电方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
交流二次回路中性线阻抗对母线差动保护的影响及防范措施;邓旭阳等;《电力系统保护与控制》;20130801;第41卷(第15期);正文第148-153页 *
基于模型参数识别的母线保护原理;索南加乐等;《中国电机工程学报》;20100805;第30卷(第22期);正文第92-99页 *

Also Published As

Publication number Publication date
CN105470931A (zh) 2016-04-06

Similar Documents

Publication Publication Date Title
CN110488154B (zh) 一种用于调度主站端的小电流接地选线方法
CN107332216B (zh) 一种采用自同步技术的母线差动保护方法及装置
CN102768324A (zh) 小电流接地系统单相接地故障区段定位方法
CN107045093B (zh) 基于快速s变换的小电流单相接地故障选线方法
CN106786427B (zh) 基于波形相关的风电场集电线路电流保护方法
CN110783896B (zh) 一种基于弱故障主被动联合检测的配网单相接地保护方法
CN109975661A (zh) 一种基于斯皮尔曼相关系数的输电线路故障检测方法
CN103267927A (zh) 一种利用工频分量小波系数直线拟合检测的小电流接地系统故障选线方法
CN102608499A (zh) 一种差分滤波抑制不平衡电流的小电流选线装置及控制方法
CN106019173B (zh) 一种应用于电压源换流器的实时故障检测方法
CN102122811A (zh) 数字化变电站电子式互感器双a/d采样的保护装置启动元件
CN103532160A (zh) 一种研究网架结构对改善受端交直流相互影响的分析方法
CN113295966B (zh) 早期故障类型的识别方法和装置
CN105470931B (zh) 一种不受母线数据不同步影响的母线差动保护方法
CN114062846A (zh) 一种交流输电线路保护方法及系统
CN106469905B (zh) 基于配网“三道防线”的配网保护方法
CN115513917B (zh) 一种含分布式电源配电网的综合电流幅值差动保护新方法
CN109001589B (zh) 基于非全相运行故障录波数据的线路参数时域计算方法
McDermott et al. Protection of distribution circuits with high penetration of solar pv: Distance, learning, and estimation-based methods
CN115663766A (zh) 考虑电流变化量相似性的风电场送出线路保护方法
CN113315100A (zh) 一种基于卷积神经网络的微电网保护方法与系统
CN107968387B (zh) 基于阻抗平面分析正序极化电压的继电器控制方法
Venkatanagaraju et al. Synchrophasor-Based Wide Area Protection Scheme for System Stressed Conditions to Avoid Power System Blackout: Wide Area Protection Scheme for System Stressed Conditions
CN112345876A (zh) 一种适用于间隔式dtu的故障定位方法及系统
CN104795805B (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
CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: 100031 No. 86 West Chang'an Avenue, Beijing, Xicheng District

Patentee after: STATE GRID CORPORATION OF CHINA

Country or region after: China

Patentee after: STATE GRID SHAANXI ELECTRIC POWER Co.

Patentee after: XI'AN JIAOTONG University

Address before: 100031 No. 86 West Chang'an Avenue, Beijing, Xicheng District

Patentee before: State Grid Corporation of China

Country or region before: China

Patentee before: STATE GRID SHAANXI ELECTRIC POWER Co.

Patentee before: XI'AN JIAOTONG University

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20240808

Address after: 100031 No. 86 West Chang'an Avenue, Beijing, Xicheng District

Patentee after: STATE GRID CORPORATION OF CHINA

Country or region after: China

Patentee after: State Grid Shaanxi Electric Power Co.,Ltd.

Patentee after: XI'AN JIAOTONG University

Address before: 100031 No. 86 West Chang'an Avenue, Beijing, Xicheng District

Patentee before: STATE GRID CORPORATION OF CHINA

Country or region before: China

Patentee before: STATE GRID SHAANXI ELECTRIC POWER Co.

Patentee before: XI'AN JIAOTONG University