CN103779846A - Power grid line wide-area current protection method based on admittance matrix fast correction - Google Patents

Power grid line wide-area current protection method based on admittance matrix fast correction Download PDF

Info

Publication number
CN103779846A
CN103779846A CN201410053653.6A CN201410053653A CN103779846A CN 103779846 A CN103779846 A CN 103779846A CN 201410053653 A CN201410053653 A CN 201410053653A CN 103779846 A CN103779846 A CN 103779846A
Authority
CN
China
Prior art keywords
circuit
matrix
centerdot
delta
power grid
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
Application number
CN201410053653.6A
Other languages
Chinese (zh)
Other versions
CN103779846B (en
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 Fujian Electric Power Co Ltd
Putian Power Supply Co of State Grid Fujian Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Fujian Electric Power Co Ltd
Putian Power Supply Co of State Grid Fujian 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 Fujian Electric Power Co Ltd, Putian Power Supply Co of State Grid Fujian Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201410053653.6A priority Critical patent/CN103779846B/en
Publication of CN103779846A publication Critical patent/CN103779846A/en
Application granted granted Critical
Publication of CN103779846B publication Critical patent/CN103779846B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a power grid line wide-area current protection method based on admittance matrix fast correction. Synchronous measurement of voltage and current phasors of each transformer substation of a power grid is achieved through an existing power-grid wide-area measurement system of the power grid, and power-grid synchronous data are acquired; if a single circuit of the power grid is opened, a middle matrix, a variable matrix, the voltage variation matrix of each transformer substation of the power grid and the variation of the current passing through each circuit of the power grid will be calculated by a wide-area back-up protection system, and wide-area current protection criteria are constituted by the variation of the current passing through each circuit of the power grid and the variation, measured by the wide-area measurement system, of the current passing through each circuit of the power grid. The whole calculation process is finished within 200 ms, the circuits can be reliably locked before circuit overload protection false operation, the phenomenon that after a single circuit of the power grid is opened, overload of other circuits is generated and the other circuits are cut off by mistake due to flow transferring is effectively avoided, fault circuits can be quickly opened before back-up circuit protection motion so as to be isolated from the power grid, and operation safety and stability of the power grid are improved.

Description

The power network line wide area current protection method of revising fast based on admittance matrix
Technical field
The present invention relates to Relay Protection Technology in Power System field, specifically relate to a kind of power network line wide area current protection method of revising fast based on admittance matrix.
Background technology
The task of relaying protection is fault and the irregular operating state of reflection protected element, can be divided into main protection and backup protection according to function.Wherein the task of excision protected element fault fast is only being undertaken in main protection, is not substantially subject to the impact of system operation mode.Along with the application of optical-fibre channel, high pressure, extra high voltage system main protection technology based on Two-Terminal Electrical Quantities are ripe day by day.And the task of doing the backup protection far away of adjacent electric component is also being undertaken in backup protection, owing to existing the possibility of transforming plant DC power-fail, this responsibility is to cancel.The information of protection installation place is only reacted in existing electrical network backup protection, is subject to the impact of power network topology annexation and operational mode.For guaranteeing its reliability, have to be configured and adjust according to the harshest situation; For guaranteeing its selectivity, have to sacrifice rapidity and the sensitivity of backup protection.Because electric network composition is increasingly sophisticated, just cause following potential safety hazard simultaneously:
1) backup protection matching relationship complexity, operate time is long.When serious, likely do not meet the desired critical clearing time of stability of power system, and then become large power grid security hidden danger.
2) backup protection configuration is large with the difficulty of adjusting, and can not follow the tracks of the variation of power system operation mode, even may occur protecting mismatch or under-sensitive situation.
3) backup protection can not be distinguished the trend transfer causing after internal fault and failure removal, and this likely causes the backup protection cascading trip in heavy load situation.
In the accident of having a power failure on a large scale, although just occurred that electric power system backup protection, according to design principle correct operation, objectively but accelerated the phenomenon of system crash.Therefore the backup protection system for damping of, considering electrical network global information is studied extremely urgent.
In recent years, the drawing-in system information in backup protection that appears as of wide area synchronized measurement system (WAMS) provides possibility.WAMS mono-is the electric phasor that can synchronously obtain in the whole network, has realized the monitoring of dynamic process of electrical power system; The 2nd, the renewal speed of measurement amount was tapered to a few tens of milliseconds by several seconds, create condition for realizing dynamic process of electrical power system control, make the design of carrying out backup protection from the angle of electrical network total optimization become possibility.
Summary of the invention
The object of the invention is to overcome the deficiency that prior art exists, a kind of power network line wide area current protection method of revising fast based on admittance matrix is provided.The method can be by reliable circuit locking before circuit overload is protected misoperation; after effectively preventing the tripping operation of electrical network uniline, trend transfer causes that All other routes overload is excised by mistake; can be by faulty line fast trip before line backup protection action; in time faulty line and electrical network are kept apart, improve the security and stability of operation of power networks.
The object of the invention is to be achieved through the following technical solutions:
The power network line wide area current protection method of revising fast based on admittance matrix, is characterized in that, comprises following sequential steps:
(1) wide area backup protection system utilizes the existing WAMS of electrical network to realize the synchro measure of the each transformer substation voltage of electrical network, electric current phasor, and Real-time Obtaining synchronized data, if after the circuit breaker tripping at circuit ij two ends, enter step (2); Wherein, circuit ij connects i transformer station and j transformer station; I=1 ..., n; J=1 ..., n; N is electricity grid substation total number;
(2) wide area backup protection system-computed intermediary matrix W=Y-1M; Wherein, Y is grid nodes admittance matrix, for n × n ties up matrix; M=[0...1...-1...0] tfor n × 1 dimension matrix, in Metzler matrix, capable the 1st column element of i is 1, and in Metzler matrix, capable the 1st column element of j is-1, and other elements are 0; I=1 ..., n; J=1 ..., n; N is electricity grid substation total number;
(3) wide area backup protection system-computed matrix of variables wherein, E is unit matrix, is 1 × 1 dimension matrix;
Figure BDA0000466710210000022
be 1 × n dimension matrix, in J matrix, the 1st row i column element is-y ij, in J matrix, the 1st row j column element is y ij, other elements are 0;
Figure BDA0000466710210000023
before circuit breaker tripping for circuit ij two ends, the each transformer substation voltage phase of the electrical network moment matrix that WAMS measures, is n × 1 dimension matrix; y ijfor the line admittance of circuit ij;
(4) after the circuit breaker tripping at wide area backup protection system-computed circuit ij two ends, the each transformer substation voltage variable quantity of electrical network matrix Δ U · = - WZ ; Wherein, Δ U · = Δ U · 1 Δ U · 2 Δ U · 3 · · · Δ U · n - 2 Δ U · n - 1 Δ U · n T For n × 1 dimension matrix, Δ U · 1 , Δ U · 3 , · · · , Δ U · n - 2 , Δ U · n - 1 , Δ U · n Be respectively after the circuit breaker tripping at circuit ij two ends, be numbered the the 1st, 2,3 ..., the voltage variety of n-2, n-1, n transformer station;
(5) current change quantity of circuit lk after the circuit breaker tripping at wide area backup protection system-computed circuit ij two ends
Figure BDA0000466710210000027
and judgement
Figure BDA0000466710210000028
whether set up, if set up, wide area backup protection system judges that circuit lk is faulty line, sends trip signal to the circuit breaker at circuit lk two ends; Wherein, l=1 ..., n; K=1 ..., n; N is electricity grid substation total number; z lkfor the line impedance of circuit lk; I setfor setting current threshold value;
Figure BDA0000466710210000029
before circuit breaker tripping for circuit ij two ends, WAMS measures the electric current flowing through on circuit lk; H=[0...1...-1...0] tfor n × 1 dimension matrix, in H matrix, capable the 1st column element of l is 1, and in H matrix, capable the 1st column element of k is-1, and other elements are all 0; after circuit breaker tripping for circuit ij two ends, WAMS measures the electric current flowing through on circuit lk; Circuit lk connects l transformer station and k transformer station; Circuit lk is different circuits from circuit ij.
The present invention compared with prior art, has following positive achievement:
The inventive method utilizes the existing WAMS of electrical network to realize the synchro measure of the each transformer substation voltage of electrical network, electric current phasor, obtains synchronized data; If there is the tripping operation of electrical network uniline; wide area backup protection system-computed intermediary matrix; calculate matrix of variables; calculate the each transformer substation voltage variable quantity of electrical network matrix; calculate the current change quantity that each line flows of electrical network crossed, utilize current change quantity and the WAMS that each the line flows of electrical network calculating crossed to measure each line current variable quantity formation power network line wide area current protection criterion of electrical network.The inventive method is in the time that the rear topological structure of electric of electrical network uniline tripping operation changes, directly calculate rear each the current change quantity that line flows is crossed of electrical network uniline tripping operation, whole computational process can complete within the 200ms time, protection quick action, faster than the existing line backup protection algorithm based on this locality amount decision-making, can be by reliable circuit locking before circuit overload is protected misoperation, after effectively preventing the tripping operation of electrical network uniline, trend transfer causes that All other routes overload is excised by mistake, can be by faulty line fast trip before line backup protection action, in time faulty line and electrical network are kept apart, improve the security and stability of operation of power networks.
Accompanying drawing explanation
Fig. 1 is application circuit transmission system schematic diagram of the present invention.
Embodiment
According to Figure of description, technical scheme of the present invention is expressed in further detail below.
Wide area backup protection system utilizes the existing WAMS of electrical network to realize the synchro measure of the each transformer substation voltage of electrical network, electric current phasor, Real-time Obtaining synchronized data.
If after the circuit breaker tripping at circuit ij two ends, wide area backup protection system-computed intermediary matrix W=Y -1m; Wherein, Y is grid nodes admittance matrix, for n × n ties up matrix; M=[0...1...-1...0] tfor n × 1 dimension matrix, in Metzler matrix, capable the 1st column element of i is 1, and in Metzler matrix, capable the 1st column element of j is-1, and other elements are 0; I=1 ..., n; J=1 ..., n; N is electricity grid substation total number; Circuit ij connects i transformer station and j transformer station.
Wide area backup protection system-computed matrix of variables
Figure BDA0000466710210000031
wherein, E is unit matrix, is 1 × 1 dimension matrix;
Figure BDA0000466710210000032
be 1 × n dimension matrix, in J matrix, the 1st row i column element is-y ij, in J matrix, the 1st row j column element is yij, other elements are 0;
Figure BDA0000466710210000033
before circuit breaker tripping for circuit ij two ends, the each transformer substation voltage phase of the electrical network moment matrix that WAMS measures, is n × 1 dimension matrix; y ijfor the line admittance of circuit ij.
After the circuit breaker tripping at wide area backup protection system-computed circuit ij two ends, the each transformer substation voltage variable quantity of electrical network matrix Δ U · = - WZ ; Wherein, Δ U · = Δ U · 1 Δ U · 2 Δ U · 3 · · · Δ U · n - 2 Δ U · n - 1 Δ U · n T For n × 1 dimension matrix, Δ U · 1 , Δ U · 3 , · · · , Δ U · n - 2 , Δ U · n - 1 , Δ U · n Be respectively after the circuit breaker tripping at circuit ij two ends, be numbered the the 1st, 2,3 ..., the voltage variety of n-2, n-1, n transformer station.
The current change quantity of circuit lk after the circuit breaker tripping at wide area backup protection system-computed circuit ij two ends
Figure BDA0000466710210000044
and judgement whether set up, if set up, wide area backup protection system judges that circuit lk is faulty line, sends trip signal to the circuit breaker at circuit lk two ends; Wherein, l=1 ..., n; K=1 ..., n; N is electricity grid substation total number; z lkfor the line impedance of circuit lk; I setfor setting current threshold value;
Figure BDA0000466710210000046
before circuit breaker tripping for circuit ij two ends, WAMS measures the electric current flowing through on circuit lk; H=[0...1.. ...-1...0] tfor n × 1 dimension matrix, in H matrix, capable the 1st column element of l is 1, and in H matrix, capable the 1st column element of k is-1, and other elements are all 0;
Figure BDA0000466710210000047
after circuit breaker tripping for circuit ij two ends, WAMS measures the electric current flowing through on circuit lk; Circuit lk connects l transformer station and k transformer station; Circuit lk is different circuits from circuit ij.
The inventive method utilizes the existing WAMS of electrical network to realize the synchro measure of the each transformer substation voltage of electrical network, electric current phasor, obtains synchronized data; If there is the tripping operation of electrical network uniline; wide area backup protection system-computed intermediary matrix; calculate matrix of variables; calculate the each transformer substation voltage variable quantity of electrical network matrix; calculate the current change quantity that each line flows of electrical network crossed, utilize current change quantity and the WAMS that each the line flows of electrical network calculating crossed to measure each line current variable quantity formation power network line wide area current protection criterion of electrical network.
The inventive method is in the time that the rear topological structure of electric of electrical network uniline tripping operation changes, directly calculate rear each the current change quantity that line flows is crossed of electrical network uniline tripping operation, whole computational process can complete within the 200ms time, protection quick action, faster than the existing line backup protection algorithm based on this locality amount decision-making, can be by reliable circuit locking before circuit overload is protected misoperation, after effectively preventing the tripping operation of electrical network uniline, trend transfer causes that All other routes overload is excised by mistake, can be by faulty line fast trip before line backup protection action, in time faulty line and electrical network are kept apart, improve the security and stability of operation of power networks.
The foregoing is only preferred embodiment of the present invention; but protection scope of the present invention is not limited to this; any be familiar with those skilled in the art the present invention disclose technical scope in, the variation that can expect easily or replacement, within all should being encompassed in protection scope of the present invention.

Claims (1)

1. the power network line wide area current protection method of revising fast based on admittance matrix, comprises following sequential steps:
(1) wide area backup protection system utilizes the existing WAMS of electrical network to realize the synchro measure of the each transformer substation voltage of electrical network, electric current phasor, and Real-time Obtaining synchronized data, if after the circuit breaker tripping at circuit ij two ends, enter step (2); Wherein, circuit ij connects i transformer station and j transformer station; I=1 ..., n; J=1 ..., n; N is electricity grid substation total number;
(2) wide area backup protection system-computed intermediary matrix W=Y -1m; Wherein, Y is grid nodes admittance matrix, for n × n ties up matrix; M=[0...1...-1...0] tfor n × 1 dimension matrix, in Metzler matrix, capable the 1st column element of i is 1, and in Metzler matrix, capable the 1st column element of j is-1, and other elements are 0; I=1 ..., n; J=1 ..., n; N is electricity grid substation total number;
(3) wide area backup protection system-computed matrix of variables
Figure FDA0000466710200000011
wherein, E is unit matrix, is 1 × 1 dimension matrix;
Figure FDA0000466710200000012
be 1 × n dimension matrix, in J matrix, the 1st row i column element is-y ij, in J matrix, the 1st row j column element is y ij, other elements are 0;
Figure FDA0000466710200000013
before circuit breaker tripping for circuit ij two ends, the each transformer substation voltage phase of the electrical network moment matrix that WAMS measures, is n × 1 dimension matrix; y ijfor the line admittance of circuit ij;
(4) after the circuit breaker tripping at wide area backup protection system-computed circuit ij two ends, the each transformer substation voltage variable quantity of electrical network matrix Δ U · = - WZ ; Wherein, Δ U · = Δ U · 1 Δ U · 2 Δ U · 3 · · · Δ U · n - 2 Δ U · n - 1 Δ U · n T For n × 1 dimension matrix, Δ U · 1 , Δ U · 3 , · · · , Δ U · n - 2 , Δ U · n - 1 , Δ U · n Be respectively after the circuit breaker tripping at circuit ij two ends, be numbered the the 1st, 2,3 ..., n-2, n-1, n transformer station voltage variety;
(5) current change quantity of circuit lk after the circuit breaker tripping at wide area backup protection system-computed circuit ij two ends
Figure FDA0000466710200000017
and judgement
Figure FDA0000466710200000018
whether set up, if set up, wide area backup protection system judges that circuit lk is faulty line, sends trip signal to the circuit breaker at circuit lk two ends; Wherein, l=1 ..., n; K=1 ..., n; N is electricity grid substation total number; z lkfor the line impedance of circuit lk; I setfor setting current threshold value;
Figure FDA0000466710200000019
before circuit breaker tripping for circuit ij two ends, WAMS measures the electric current flowing through on circuit lk; H=[0...1...-1...0] tfor n × 1 dimension matrix, in H matrix, capable the 1st column element of l is 1, and in H matrix, capable the 1st column element of k is-1, and other elements are all 0; after circuit breaker tripping for circuit ij two ends, WAMS measures the electric current flowing through on circuit lk; Circuit lk connects l transformer station and k transformer station; Circuit lk is different circuits from circuit ij.
CN201410053653.6A 2014-02-18 2014-02-18 The power network line wide area current protection method of revising fast based on admittance matrix Active CN103779846B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410053653.6A CN103779846B (en) 2014-02-18 2014-02-18 The power network line wide area current protection method of revising fast based on admittance matrix

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410053653.6A CN103779846B (en) 2014-02-18 2014-02-18 The power network line wide area current protection method of revising fast based on admittance matrix

Publications (2)

Publication Number Publication Date
CN103779846A true CN103779846A (en) 2014-05-07
CN103779846B CN103779846B (en) 2016-05-11

Family

ID=50571816

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410053653.6A Active CN103779846B (en) 2014-02-18 2014-02-18 The power network line wide area current protection method of revising fast based on admittance matrix

Country Status (1)

Country Link
CN (1) CN103779846B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104065066A (en) * 2014-07-16 2014-09-24 国家电网公司 Power grid load flow calculation method based on admittance matrix quick fix
CN104092194A (en) * 2014-07-15 2014-10-08 国家电网公司 Power grid circuit wide-area current protection method
CN104332972A (en) * 2014-11-12 2015-02-04 国家电网公司 Grid line wide-area backup protection method applicable to single-line tripping
CN104332971A (en) * 2014-11-12 2015-02-04 国家电网公司 Grid line wide-area backup protection method applicable to double-line tripping

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0665625A2 (en) * 1994-01-24 1995-08-02 Licentia Patent-Verwaltungs-GmbH Method for obtaining an impedance value and for processing in a distance protection device
CN101409442A (en) * 2008-11-25 2009-04-15 中国南方电网有限责任公司 Rapid line backup protection method with reliable selection based on wide area measuring system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0665625A2 (en) * 1994-01-24 1995-08-02 Licentia Patent-Verwaltungs-GmbH Method for obtaining an impedance value and for processing in a distance protection device
CN101409442A (en) * 2008-11-25 2009-04-15 中国南方电网有限责任公司 Rapid line backup protection method with reliable selection based on wide area measuring system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
马静等: "快速适应电网结构变化的广域后备保护新方法", 《电网技术》, vol. 35, no. 7, 31 July 2011 (2011-07-31), pages 214 - 220 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104092194A (en) * 2014-07-15 2014-10-08 国家电网公司 Power grid circuit wide-area current protection method
CN104065066A (en) * 2014-07-16 2014-09-24 国家电网公司 Power grid load flow calculation method based on admittance matrix quick fix
CN104332972A (en) * 2014-11-12 2015-02-04 国家电网公司 Grid line wide-area backup protection method applicable to single-line tripping
CN104332971A (en) * 2014-11-12 2015-02-04 国家电网公司 Grid line wide-area backup protection method applicable to double-line tripping

Also Published As

Publication number Publication date
CN103779846B (en) 2016-05-11

Similar Documents

Publication Publication Date Title
CN103762570B (en) A kind of power network line relay protecting method that is applicable to dual line fault
US20200412167A1 (en) Method and system for fast reconfiguration of power supply network in tens of milliseconds after power grid failure
CN101409442B (en) Rapid line backup protection method with reliable selection based on wide area measuring system
CN101741072B (en) Method for fault component-based wide area backup protection
CN101621216B (en) Data sharing type area protection system based on IEC 61850
CN101872964B (en) Wide area measurement system based back-up protection method of multi-terminal high-voltage power transmission area
CN102521667A (en) Probability assessing method for stage type protection operational risk of electric power system
CN101621217B (en) Protection system for area electric network data sharing
CN101350520A (en) Protection control system and method base on electric network synthesis information
CN103779846B (en) The power network line wide area current protection method of revising fast based on admittance matrix
CN103778274A (en) Ultrahigh-voltage circuit relaying protection system reliability assessment method with covert faults taken into consideration
CN104485736A (en) Intelligent protection center-based wide-area backup protection system and method
CN103580009A (en) Self-adaptive overload recognition system based on composite phasor plane and method thereof
CN104253421A (en) Relay protection optimal configuration method for double-circuit lines on same pole based on substation area information sharing
CN103683229B (en) Method for protecting power grid wide-area relay based on unit output power and load change influence
CN103715669A (en) Distributed wide-area self-adaptive inter-phase backup protection system and method
CN101917055B (en) Branch circuit similarity-based wide area backup protection method
CN102570426B (en) Method for realizing back-up protection in wide area protection by using differential protection
CN106385011B (en) A kind of fault coverage identification device of Intelligent region spare power automatic switching system and recognition methods
CN103779853A (en) Grid circuit overload maloperation-preventing relay protection method based on wide-area information
CN104092194A (en) Power grid circuit wide-area current protection method
CN104037744A (en) Power grid wide-area relay protection method based on tide transient component predication
JP5450149B2 (en) Distribution line ground fault protection relay system
CN102611084B (en) Self-adaptive voltage protection method
CN104332972A (en) Grid line wide-area backup protection method applicable to single-line tripping

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

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

Co-patentee after: STATE GRID FUJIAN ELECTRIC POWER Co.,Ltd.

Patentee after: STATE GRID CORPORATION OF CHINA

Co-patentee after: STATE GRID PUTIAN ELECTRIC POWER SUPPLY Co.

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

Co-patentee before: STATE GRID FUJIAN ELECTRIC POWER Co.,Ltd.

Patentee before: State Grid Corporation of China

Co-patentee before: STATE GRID PUTIAN ELECTRIC POWER SUPPLY Co.

CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: 351100 No. 999, South Garden Road, Xia Lin Street, Chengxiang District, Putian, Fujian

Co-patentee after: STATE GRID FUJIAN ELECTRIC POWER Co.,Ltd.

Patentee after: STATE GRID PUTIAN ELECTRIC POWER SUPPLY Co.

Co-patentee after: STATE GRID CORPORATION OF CHINA

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

Co-patentee before: STATE GRID FUJIAN ELECTRIC POWER Co.,Ltd.

Patentee before: State Grid Corporation of China

Co-patentee before: STATE GRID PUTIAN ELECTRIC POWER SUPPLY Co.