patents.google.com

CN101067641A - Single-ended fault location method of line against the influence of distributed capacitive current and transition resistance - Google Patents

  • ️Wed Nov 07 2007
Single-ended fault location method of line against the influence of distributed capacitive current and transition resistance Download PDF

Info

Publication number
CN101067641A
CN101067641A CN 200710100215 CN200710100215A CN101067641A CN 101067641 A CN101067641 A CN 101067641A CN 200710100215 CN200710100215 CN 200710100215 CN 200710100215 A CN200710100215 A CN 200710100215A CN 101067641 A CN101067641 A CN 101067641A Authority
CN
China
Prior art keywords
fault
voltage
zero
line
omega
Prior art date
2007-06-06
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
CN 200710100215
Other languages
Chinese (zh)
Other versions
CN100526895C (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.)
Tsinghua University
Original Assignee
Tsinghua University
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.)
2007-06-06
Filing date
2007-06-06
Publication date
2007-11-07
2007-06-06 Application filed by Tsinghua University filed Critical Tsinghua University
2007-06-06 Priority to CNB2007101002150A priority Critical patent/CN100526895C/en
2007-11-07 Publication of CN101067641A publication Critical patent/CN101067641A/en
2009-08-12 Application granted granted Critical
2009-08-12 Publication of CN100526895C publication Critical patent/CN100526895C/en
Status Active legal-status Critical Current
2027-06-06 Anticipated expiration legal-status Critical

Links

  • 238000000034 method Methods 0.000 title claims abstract description 66
  • 230000007704 transition Effects 0.000 title claims abstract description 20
  • 230000005540 biological transmission Effects 0.000 claims description 20
  • 238000009434 installation Methods 0.000 claims description 16
  • 230000004807 localization Effects 0.000 description 26
  • 230000005611 electricity Effects 0.000 description 10
  • 238000004422 calculation algorithm Methods 0.000 description 4
  • 238000005259 measurement Methods 0.000 description 4
  • 238000004088 simulation Methods 0.000 description 2
  • 230000001131 transforming effect Effects 0.000 description 2
  • 238000004364 calculation method Methods 0.000 description 1
  • 238000004891 communication Methods 0.000 description 1
  • 238000010586 diagram Methods 0.000 description 1
  • 230000000694 effects Effects 0.000 description 1
  • 238000005516 engineering process Methods 0.000 description 1
  • 238000012821 model calculation Methods 0.000 description 1

Images

Landscapes

  • Emergency Protection Circuit Devices (AREA)

Abstract

The invention belongs to the electrical power system domain, specially relates to the line single end fault distance measuring method of the anti- distributed capacity electric current and the transition resistance influence. Including: surveys the transformer substation protection installment place fault phase voltage, the phase current, the zero sequence voltage, the zero sequence electric current and the negative sequence electric current as the input value; using the fault voltage cross zero time to measure the resistance characteristic which is not influence by transition resistance, from the protected line begin end start, computing the measure voltage real part computed value of the fault point voltage cross zero point time, compared with the measured value of the measure voltage real part, computing the error; and taking the delta-S as the step size, in turn computing the setting range which is protected by the send tripping signal, if cannot obtain the protect tripping signal, searching the total track length which is protected, taking the minimum point as the fault point. The invention method is not influenced by the distributed capacity electric current and transition resistance, does not have the false root problem of solving equation method and non convergence question of the iterative method, and has the very high practical value.

Description

The line one-end fault ranging method of anti-capacitance current and transition resistance influence

Technical field

The invention belongs to field of power, the line one-end fault ranging method of particularly a kind of anti-capacitance current and transition resistance influence.

Background technology

Transmission line of electricity is being undertaken the vital task that transmits electric energy, and its fault directly threatens the safe operation of electric system.Especially, the electric pressure of China develops to extra-high voltage (is representative with 1000kV) from UHV (ultra-high voltage) (is representative with 500kV), the transmission line of electricity capacity further increases, fault localization is got rid of line fault and is restored electricity as early as possible for quickening accurately, thereby the stability, the stable operation of assurance security of system that improve electric system have great importance.

Fault distance-finding method is divided into two kinds: single-ended method and both-end method from using the angle of electric parameters.The both-end method need be used the electric information at circuit two ends, can realize accurate localization of fault on the principle, but in actual applications, still is subjected to the influence of two ends sample frequency difference or phase deviation.Single-ended method does not need extra equipment, only adopts local electric parameters to carry out fault localization, is not subjected to the restriction of mechanics of communication, is the focus that industry is paid close attention to always.

The one-end fault ranging algorithm is divided into two classes from the angle of principle: impedance method and traveling wave method.Traveling wave method is realized fault localization by the time of measuring row wave-wave head commute between observation station and trouble spot, but the single-ended traveling wave method can not effectively be discerned the reflected traveling wave that comes from the trouble spot and come from the reflected traveling wave of other wave impedance discontinuous point, so the single-ended traveling wave method does not obtain the essence application.Impedance method: calculate observation station to the line impedance value between the trouble spot according to voltage, galvanometer after the fault, lead owing to ignored distributed capacitance, the electricity of transmission line of electricity, the proportional example relation of line impedance value and fault distance, thus realize localization of fault; But the single-ended impedance method can not be eliminated the influence of fault transition resistance and offside system impedance fully, precision is relatively poor than the both-end method, a lot of scholars are at this problem, proposed abundant the improving one's methods of kind, mainly contained: power frequency impedance method (1), separated Differential Equation Algorithm (2), zero-sequence current phase place revised law (3), fault current phase place revised law (4), separate quadratic equation method (5) etc.Method (1), (2) are assumed fault point electric current and observation station electric current same-phase all, method (3), (4), (5) all need known offside system impedance value, accurate localization of fault, yet system operation mode is changeable, be difficult to accurately obtain the peer-to-peer system resistance value, so the fault localization precision is not very high.But single-ended impedance fault localization algorithm does not need the information of offside, and is low to hardware requirement, is easy to realize; Though distance accuracy is not high, algorithm stability is fine, so has obtained at the scene using widely.The microcomputer protecting device of middle operation all adopts the single-ended impedance method as its fault localization element basically at the scene at present; utilize the stability of impedance method range finding simultaneously; the comprehensive combined method range finding of using single-ended impedance method and traveling wave method has also greatly promoted the performance that single-ended traveling wave is found range.

The single-ended impedance fault distance-finding method is widely used on high pressure (is representative with 220kV), supertension line, and its prerequisite is to think to measure impedance and the proportional example relation of fault distance that promptly the distributed capacitance of circuit can be ignored.This is because high pressure, supertension line length are generally less than 400km, and it is little to ignore the error that the influence of distributed capacitance brings, and can satisfy the requirement of rig-site utilization.And for the ultra-high/extra-high voltage long transmission line of the above electric pressure of 750kV, the electric current of the capacitance current of circuit and circuit transmission natural power is compared and can be reached more than 76%, and the error of still ignore distributed capacitance this moment, bringing with lumpy line Model Calculation fault distance can not be accepted by the scene.

Correlative study shows: with the distributed parameter transmission line model modeling, and the line impedance Z between observation station and the trouble spot MeaBe the tanh relation with fault distance, that is:

Z mea = Z z 1 th γ · 1 l ,

Z wherein C1Be the positive sequence wave impedance,

Figure A20071010021500052

Be the positive sequence propagation coefficient, they all are nonlinear functions of unit length line resistance, inductance, capacitance.The hyperbolic tangent function characteristic has determined its anti-transition resistance poor ability, and the additional measurement impedance that transition resistance brings will seriously influence the precision of fault localization.Traditional impedance method improvement project: carry out fault localization as adopting reactive component, partly eliminate the influence of transition resistance, no longer suitable for fault distance-finding method with the distributed parameter transmission line model modeling, mainly be because line impedance Z MeaBe the nonlinear function of unit length line resistance, inductance, capacitance, be difficult to decoupling zero; Other scheme is as hypothesis peer-to-peer system impedance known or assumed fault point electric current and observation station electric current same-phase, on the ultra-high/extra-high voltage long transmission line of the above electric pressure of 750kV, the error of bringing can't be accepted, and therefore traditional single-ended impedance fault distance-finding method is difficult to be applied directly on the ultra-high/extra-high voltage long transmission line.

Summary of the invention

The objective of the invention is for overcoming the weak point of prior art, propose a kind of new anti-capacitance current and the line one-end fault ranging method of transition resistance influence, the physical model of this method adopts the distribution parameter modeling, is not subjected to the influence of capacitance current; This method utilizes fault point voltage zero crossing data computation constantly to measure resistance value, is not subjected to the influence of transition resistance; This method is a kind of method of search type, does not have the pseudo-root problem of the method for solving equation and the not convergence problem of process of iteration, has very high practical value.

The line one-end fault ranging method of a kind of anti-capacitance current that the present invention proposes and transition resistance influence may further comprise the steps:

1, measuring circuit is in transforming plant protecting installation place fault phase voltage phasor Fault phase current phasor The residual voltage phasor

Figure A20071010021500055

The zero-sequence current phasor The negative-sequence current phasor

Figure A20071010021500057

As input quantity;

Wherein  is that fault is separate: A phase, B phase or C phase;

2) calculating the fault point voltage zero crossing is constantly:

Figure A20071010021500058

Wherein α is the initial phase angle of route protection installation place negative-sequence current, and ω is the specified angular frequency value of electric system;

3) calculate the fault point voltage zero crossing constantly, the real part measured value of protection installation place measuring voltage is:

Measuring voltage

Figure A20071010021500059

Wherein: U RelayBe the measuring voltage amplitude, θ is the initial phase angle of measuring voltage; In the fault point voltage zero crossing moment, the real part measured value of measuring voltage is: U Mea_meter=U Relaycos(90 °-α+θ)

4) fault distance is taken as initial value l Fault, calculate the protection installation place and measure electric current:

Measure electric current

Figure A200710100215000510

Wherein: I RelayFor measuring current amplitude, β is for measuring the initial phase angle of electric current;

P = Z c 0 Z c 1 ( T · ch γ · 1 l fault + sh γ · 0 l fault - T · ch γ · 0 l fault sh γ · 1 l fault ) - 1

P is the zero sequence current compensation factor based on the circuit distribution parameter, wherein:

Z C1Be the positive sequence wave impedance:

Z c 1 = ( R 1 + jω L 1 ) / ( G 1 + jω C 1 ) ,

R 1, L 1, G 1, C 1The positive sequence resistance, inductance, the electricity that are respectively the unit length circuit are led and capacitance;

Z C0Be the zero sequence wave impedance:

Z c 0 = ( R 0 + jω L 0 ) / ( G 0 + jω C 0 ) ,

R 0, L 0, G 0, C 0The zero sequence resistance, inductance, the electricity that are respectively the unit length circuit are led and capacitance;

Figure A20071010021500064

Be the positive sequence propagation coefficient:

γ · 1 = ( R 1 + jω L 1 ) ( G 1 + jω C 1 ) ;

Be the zero sequence propagation coefficient:

γ · 0 = ( R 0 + jω L 0 ) ( G 0 + jω C 0 ) ;

T is the system equivalent zero sequence impedance based on distributed parameter model:

T = U · 0 Z c 0 I · 0 ;

5) according to Transmission Line Distributed Parameter computational scheme resistance value Z FaultFor:

Z fault = Z c 1 th γ · 1 l fault = Z relay ( cos ( ωt + τ ) + j sin ( ωt + τ ) )

Wherein: Z RelayBe the line impedance amplitude, τ is the line impedance angle;

6) according to the fault point voltage zero crossing constantly, the real part calculated value of protection installation place measurement electric current and line impedance value calculating measuring voltage is:

U mea_cal=I relayZ relaycos(90°-α+τ+β)

7) calculating the real part calculated value of measuring voltage and the error of measured value is:

Error=|U mea_cal-U mea_meter|

8) fault distance initial value l FaultS increases one by one with step delta, return step 4), calculate the error amount of every bit successively, (setting range is according to the route protection of the relay protection national standard principle value of adjusting until the setting range of the protection of sending out trip signal, for example I section ground distance protection is sent out trip signal, and its setting range is 80% of a protected circuit total length; Total track length is 800km; then search procedure is that step-length increases to 640km one by one with Δ S), if can't be protected trip signal, then search for the protected circuit total length; the point of getting the error minimum is the trouble spot, and this distance to the route protection installation place is a fault distance.

Characteristics of the present invention and technique effect:

The inventive method is based on the proposition of Transmission Line Distributed Parameter model, can accurately describe the physical characteristics of transmission line of electricity, has the ability of natural anti-capacitance current influence; The inventive method is utilized fault point voltage zero crossing data computation fault distance constantly, is not subjected to the influence of transition resistance; The inventive method is a kind of method of search type, does not have the pseudo-root problem of the method for solving equation and the not convergence problem of process of iteration, has very high practical value.The inventive method is applicable to the transmission line of electricity of any electric pressure, especially for 750kV and above ultra-high/extra-high voltage transmission of electricity long transmission line, the fault localization precision of application the inventive method significantly improves than other one-end fault ranging method, can satisfy on-the-spot requirement.

Description of drawings

Fig. 1 is for using the system for ultra-high voltage transmission synoptic diagram of the inventive method.

Fig. 2 compares based on the fault localization error characteristics of application of model the inventive method shown in Figure 1 and traditional power frequency impedance method fault localization error characteristics; Wherein:

(a) for using traditional power frequency impedance method fault localization error characteristics;

(b) for using the fault localization error characteristics of the inventive method;

Embodiment

The line one-end fault ranging method embodiment of anti-capacitance current that the present invention proposes and transition resistance influence is described in detail as follows:

Use a kind of 1000kV system for ultra-high voltage transmission type of the present invention as shown in Figure 1, system is a typical both end power supplying system, and two side bus are respectively M and N, and line length is 800km, and the line parameter circuit value value is as shown in table 1.Both sides system impedance parameter is as follows, and N side power supply angle falls behind M side 44 degree, and M side and N side electromotive force are respectively 1.1062 and 1.1069 times of rated voltages.The fault location device of using the inventive method is installed in the M side, and voltage, electric current are respectively from line side voltage transformer (VT) (PT), current transformer (CT), and the positive dirction of electric current is electric current is flowed to circuit by bus a direction.

Table 11000kV UHV transmission line major parameter

Line parameter circuit value Resistance (Ω/km) Reactance (Ω/km) Capacitive reactance (M Ω/km)
The positive sequence zero sequence 0.00805 0.20489 0.25913 0.74606 0.22688 0.35251

Both sides system impedance parameter is:

M side positive sequence system impedance: Z M1=4.2643+j85.14528 Ω

M side zero sequence system impedance: Z M0=98.533+j260.79 Ω

N side positive sequence system impedance: Z N1=7.9956+j159.6474 Ω

N side zero sequence system impedance: Z N0=184.749+j488.981 Ω

The line one-end fault ranging method that the present invention proposes is applicable to the range finding element and the fault location device independently of fault localization module, the fault oscillograph of any route protection.Present embodiment is an example with the fault localization module of using route protection of the present invention; with ground distance protection I section troubles inside the sample space is evaluating objects; protection domain adjust be 80% of total track length, emulation fault be the A of 420km place through 105 ohm of earth faults, the embodiment concrete steps are as follows:

1) measuring circuit is in transforming plant protecting installation place fault phase voltage phasor, phase current phasor, residual voltage phasor, zero-sequence current phasor, negative-sequence current phasor, and as input quantity, the present embodiment fault is the A phase mutually:

The A phase voltage

U · a = - 0.555 + j · 0.110 MV

The A phase current

I · a = - 2.799 + j · 1.244 kA

Residual voltage

U · 0 = 0.341 - j · 0.128 MV

Zero-sequence current

I · 0 = - 0.004 + j · 1.307 kA

Negative-sequence current

I · 2 = - 0 . 326 + j · 0.827 kA

2) the initial phase angle α of computational scheme protection installation place negative-sequence current is:

Figure A20071010021500086

Then the fault point voltage zero crossing is constantly:

Figure A20071010021500087

3) protection installation place measuring voltage is:

Figure A20071010021500088

Then calculate the fault point voltage zero crossing constantly, the real part measured value of protection measuring voltage is:

U mea_meter=U relaycos(90°-α+θ)=0.566cos(-21.53°+168.76°)=-0.476MV

4) fault distance is taken as initial value l Fault=1km, calculate and protect installation place measurement electric current to be:

The positive sequence wave impedance:

Z c 1 = ( R 1 + jω L 1 ) / ( G 1 + jω C 1 ) = 242.5 - j · 3.766 Ω

The zero sequence wave impedance:

Z c 0 = ( R 0 + jω L 0 ) / ( G 0 + jω C 0 ) = 517.55 - j · 69.775 Ω

The positive sequence propagation coefficient:

γ · 1 = ( R 1 + jω L 1 ) ( G 1 + jω C 1 ) = j · 0.0011

The zero sequence propagation coefficient:

γ · 0 = ( R 0 + jω L 0 ) ( G 0 + jω C 0 ) = 0.0002 + j · 0.0015

Utilize zero-sequence current, voltage and zero sequence wave impedance value, calculate T:

T = U · 0 Z c 0 I · 0 = - 0.121 - j · 0.520

Bring the aforementioned calculation result into P value computing formula, ask for the P value and be:

P = Z c 0 Z c 1 ( T · ch γ · 1 l fault + sh γ · 0 l fault - T · ch γ · 0 l fault sh γ · 1 l fault ) - 1 = 1.901 - j · 0.700

Therefore, obtain measuring electric current:

Figure A20071010021500091

5) according to Transmission Line Distributed Parameter computational scheme resistance value Z FaultFor:

Z fault = Z c 1 th γ · 1 l fault = Z relay ( cos ( ωt + τ ) + j sin ( ωt + τ ) )

Figure A20071010021500093

6, constantly, the real part calculated value of protection installation place measurement electric current and line impedance value calculating measuring voltage is according to the fault point voltage zero crossing:

U mea_cal=I relayZ relaycos(90°-α+τ+β)

=4.182×0.259×cos(-21.53°+116.87°+88.209°)=-0.0011MV

7) calculating the real part calculated value of measuring voltage and the error of measured value is:

Error=|U mea_cal-U mea_meter|=|-0.0011+0.476|=0.4749MV

8) fault distance initial value l Fault1km increases one by one with step-length, returns step 4, calculates the error amount of every bit successively, and until the setting range 640km of I section ground distance protection, then error amount is as shown in table 2:

The real part calculated value of table 2 measuring voltage and the error of measured value

Fault distance Error amount (* 10 5V) Fault distance Error amount (* 10 5V) Fault distance Error amount (* 10 5V) Fault distance Error amount (* 10 5V) Fault distance Error amount (* 10 5V) Fault distance Error amount (* 10 5V) Fault distance Error amount (* 10 5V)
1 4.750 93 3.758 185 2.760 277 1.741 369 0.681 461 0.442 553 0.166
2 4.740 94 3.747 186 2.749 278 1.730 370 0.669 462 0.455 554 0.167
3 4.729 95 3.736 187 2.738 279 1.718 371 0.658 463 0.467 555 0.169
4 4.718 96 3.725 188 2.727 280 1.707 372 0.646 464 0.480 556 0.170
5 4.707 97 3.715 189 2.716 281 1.696 373 0.634 465 0.493 557 0.172
6 4.696 98 3.704 190 2.705 282 1.684 374 0.622 466 0.506 558 0.173
7 4.686 99 3.693 191 2.694 283 1.673 375 0.610 467 0.518 559 0.174
8 4.675 100 3.682 192 2.683 284 1.662 376 0.598 468 0.531 560 0.176
9 4.664 101 3.671 193 2.672 285 1.651 377 0.586 469 0.544 561 0.177
10 4.653 102 3.661 194 2.662 286 1.639 378 0.575 470 0.557 562 0.179
11 4.642 103 3.650 195 2.651 287 1.628 379 0.563 471 0.569 563 0.180
12 4.631 104 3.639 196 2.640 288 1.617 380 0.551 472 0.582 564 0.181
13 4.621 105 3.628 197 2.629 289 1.605 381 0.539 473 0.595 565 0.183
14 4.610 106 3.617 198 2.618 290 1.594 382 0.527 474 0.608 566 0.184
15 4.599 107 3.607 199 2.607 291 1.583 383 0.515 475 0.620 567 0.186
16 4.588 108 3.596 200 2.596 292 1.571 384 0.503 476 0.633 568 0.187
17 4.577 109 3.585 201 2.585 293 1.560 385 0.491 477 0.646 569 0.188
18 4.567 110 3.574 202 2.574 294 1.549 386 0.479 478 0.659 570 0.190
19 4.556 111 3.563 203 2.563 295 1.537 387 0.467 479 0.672 571 0.191
20 4.545 112 3.553 204 2.552 296 1.526 388 0.455 480 0.685 572 0.193
21 4.534 113 3.542 205 2.541 297 1.515 389 0.443 481 0.698 573 0.194
22 4.523 114 3.531 206 2.530 298 1.503 390 0.431 482 0.710 574 0.195
23 4.513 115 3.520 207 2.519 299 1.492 391 0.419 483 0.723 575 0.197
24 4.502 116 3.509 208 2.508 300 1.480 392 0.407 484 0.736 576 0.198
25 4.491 117 3.499 209 2.497 301 1.469 393 0.395 485 0.749 577 0.200
26 4.480 118 3.488 210 2.486 302 1.458 394 0.383 486 0.762 578 0.201
27 4.469 119 3.477 211 2.475 303 1.446 395 0.371 487 0.775 579 0.203

Figure A20071010021500101

The point of error minimum is the 426km place, therefore thinks that the trouble spot is positioned at 426km.Calculating the fault localization relative error is:

426 - 420 800 × 100 % = 0.75 %

Use the inventive method and the fault localization error characteristics of using traditional power frequency impedance method for comparing check, the present invention is based on model shown in Figure 1 and carried out a large amount of Digital Simulations, the trouble spot is selected to be decremented to 10km gradually from 600km, and step-length is 10km; The fault transition resistance is since 5 ohm, with 50 ohm be step-length, increase to 155 ohm gradually.Simulation result as shown in Figure 2.

By Fig. 2 (a) as seen, the fault localization error characteristics of using traditional power frequency impedance method are very poor, even under the less situation of transition resistance (5 ohm, 55 ohm), the fault localization maximum relative error also can reach 30%; Along with the further increase of transition resistance value, the fault localization relative error reaches more than 50%, and the range finding result does not have credibility.

The fault localization error characteristics that adopt the inventive method are shown in Fig. 2 (b), and under the situation of little transition resistance, the fault localization maximum relative error is no more than 0.75%, realize precision ranging; Under the situation of great transition resistance, the precision of fault localization descends to some extent, but maximum relative error can satisfy the requirement of rig-site utilization less than 1.5%.

Claims (1)

1、一种抗分布电容电流和过渡电阻影响的线路单端故障测距方法,包括以下步骤:1. A line single-ended fault location method for resisting the influence of distributed capacitive current and transition resistance, comprising the following steps: 1)测量线路在变电站保护安装处故障相电压相量

Figure A2007101002150002C1

故障相电流相量

Figure A2007101002150002C2

零序电压相量 零序电流相量

Figure A2007101002150002C4

负序电流相量

Figure A2007101002150002C5

作为输入量;其中为故障相别:A相、B相或C相;
1) Measure the fault phase voltage phasor of the line at the protection installation of the substation

Figure A2007101002150002C1

Fault phase current phasor

Figure A2007101002150002C2

zero sequence voltage phasor zero sequence current phasor

Figure A2007101002150002C4

negative sequence current phasor

Figure A2007101002150002C5

As an input; where  is the fault phase: A phase, B phase or C phase;
2)计算故障点电压过零点时刻为:

Figure A2007101002150002C6

其中α为线路保护安装处负序电流的初始相角,ω为电力系统额定角频率值;
2) Calculate the zero-crossing time of the fault point voltage as:

Figure A2007101002150002C6

Where α is the initial phase angle of the negative-sequence current where the line protection is installed, and ω is the rated angular frequency value of the power system;
3)计算故障点电压过零点时刻,保护安装处测量电压的实部实测值为:3) Calculate the time when the voltage at the fault point crosses zero, and the measured value of the real part of the measured voltage at the protection installation is: 测量电压 其中:Urelay为测量电压幅值,θ为测量电压初始相角;故障点电压过零点时刻,测量电压的实部实测值为:Umea_meter=Urelaycos(90°-α+θ)Measuring voltage Among them: U relay is the measured voltage amplitude, θ is the initial phase angle of the measured voltage; when the voltage at the fault point crosses zero, the measured value of the real part of the measured voltage is: U mea_meter = U relay cos(90°-α+θ) 4)故障距离取为初始值lfault,计算保护安装处测量电流:4) The fault distance is taken as the initial value l fault , and the measured current at the protection installation is calculated: 测量电流 其中:Irelay为测量电流幅值,β为测量电流初始相角;measuring current Where: I relay is the measured current amplitude, β is the initial phase angle of the measured current; PP == ZZ cc 00 ZZ cc 11 (( TT ·· chch γγ ·&Center Dot; 11 ll faultfault ++ shsh γγ ·· 00 ll faultfault -- TT ·· chch γγ ·&Center Dot; 00 ll faultfault shsh γγ ·&Center Dot; 11 ll faultfault )) -- 11 P为基于线路分布参数的零序电流补偿系数,其中:P is the zero-sequence current compensation coefficient based on line distribution parameters, where: Zc1为正序波阻抗: Z c 1 = ( R 1 + jω L 1 ) / ( G 1 + jω C 1 ) , R1、L1、G1、C1分别为单位长度线路的正序电阻、电感、电导和电容值;Z c1 is the positive sequence wave impedance: Z c 1 = ( R 1 + jω L 1 ) / ( G 1 + jω C 1 ) , R 1 , L 1 , G 1 , and C 1 are the positive sequence resistance, inductance, conductance, and capacitance of the line per unit length, respectively; Zc0为零序波阻抗: Z c 0 = ( R 0 + jω L 0 ) / ( G 0 + jω C 0 ) , R0、L0、G0、C0分别为单位长度线路的零序电阻、电感、电导和电容值;Z c0 is the zero-sequence wave impedance: Z c 0 = ( R 0 + jω L 0 ) / ( G 0 + jω C 0 ) , R 0 , L 0 , G 0 , and C 0 are the zero-sequence resistance, inductance, conductance, and capacitance of the line per unit length, respectively;

Figure A2007101002150002C12

为正序传播系数: γ · 1 = ( R 1 + j ωL 1 ) ( G 1 + j ωC 1 ) ;

Figure A2007101002150002C12

is the positive sequence propagation coefficient: γ &Center Dot; 1 = ( R 1 + j ωL 1 ) ( G 1 + j ω C 1 ) ;

Figure A2007101002150002C14

为零序传播系数: γ · 0 = ( R 0 + j ωL 0 ) ( G 0 + j ωC 0 ) ;

Figure A2007101002150002C14

is the zero-sequence propagation coefficient: γ · 0 = ( R 0 + j ω L 0 ) ( G 0 + j ω C 0 ) ;
T为基于分布参数模型的系统等值零序阻抗: T = U · 0 Z c 0 I · 0 ; T is the equivalent zero-sequence impedance of the system based on the distributed parameter model: T = u &Center Dot; 0 Z c 0 I · 0 ; 5)根据输电线路分布参数计算线路阻抗值Zfault为:5) Calculate the line impedance value Z fault according to the transmission line distribution parameters as: ZZ faultfault == ZZ cc 11 ththe th γγ ·&Center Dot; 11 ll faultfault == ZZ relayrelay (( coscos (( ωtωt ++ ττ )) ++ jj sinsin (( ωtωt ++ ττ )) )) 其中:Zrelay为线路阻抗幅值,τ为线路阻抗角;Among them: Z relay is the line impedance amplitude, τ is the line impedance angle; 6)根据故障点电压过零点时刻,保护安装处测量电流和线路阻抗值计算测量电压的实部计算值为:6) Calculate the real part of the measured voltage according to the zero-crossing moment of the fault point voltage, the measured current at the protection installation and the line impedance value: Umea_cal=IrelayZrelaycos(90°-α+τ+β)U mea_cal =I relay Z relay cos(90°-α+τ+β) 7)计算测量电压的实部计算值与实测值的误差为:7) Calculate the error between the calculated value of the real part of the measured voltage and the measured value: Error=|Umea_cal-Umea_meter|Error=|U mea_cal -U mea_meter | 8)故障距离初始值lfault以步长ΔS逐次增加,返回步骤4),依次计算每一点的误差值,直至发跳闸信号的保护的整定范围,如果无法得到保护跳闸信号,则搜索被保护线路全长,取误差最小的点为故障点,该点至线路保护安装处的距离为故障距离。8) The initial value of the fault distance l fault is gradually increased with the step size ΔS, return to step 4), and calculate the error value of each point in turn until the setting range of the protection that sends the trip signal. If the protection trip signal cannot be obtained, search for the protected line The overall length, the point with the smallest error is taken as the fault point, and the distance from this point to the line protection installation is the fault distance.
CNB2007101002150A 2007-06-06 2007-06-06 Distributing capacitance current and transition resistance influence resisting line one-end fault ranging method Active CN100526895C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2007101002150A CN100526895C (en) 2007-06-06 2007-06-06 Distributing capacitance current and transition resistance influence resisting line one-end fault ranging method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2007101002150A CN100526895C (en) 2007-06-06 2007-06-06 Distributing capacitance current and transition resistance influence resisting line one-end fault ranging method

Publications (2)

Publication Number Publication Date
CN101067641A true CN101067641A (en) 2007-11-07
CN100526895C CN100526895C (en) 2009-08-12

Family

ID=38880253

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2007101002150A Active CN100526895C (en) 2007-06-06 2007-06-06 Distributing capacitance current and transition resistance influence resisting line one-end fault ranging method

Country Status (1)

Country Link
CN (1) CN100526895C (en)

Cited By (36)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101267107B (en) * 2008-01-08 2010-04-21 江苏省电力试验研究院有限公司 The Method of Suppressing the Single-phase Short-circuit Current of Power Network
CN102023275A (en) * 2010-09-27 2011-04-20 华北电力大学 Single-end ranging method of line based on phase mutation characteristics of positioning function
CN102129014A (en) * 2011-01-28 2011-07-20 福建省电力有限公司福州超高压输变电局 Method for realizing line phase-to-phase fault single-end distance measurement by utilizing distribution parameter model
CN102129011A (en) * 2011-01-20 2011-07-20 福建省电力有限公司福州超高压输变电局 Single-ended phase-to-phase fault location method for distributed capacitance current and fault resistance resistant line
CN102135589A (en) * 2011-02-18 2011-07-27 华北电力大学 Distance measurement method for line phase-to-phase fault distance protection
CN102147443A (en) * 2011-01-13 2011-08-10 国网电力科学研究院 Single-end distance measuring method based on self-adaptive current
CN102175954A (en) * 2011-03-09 2011-09-07 福建省电力有限公司福州超高压输变电局 Circuit inter-phase fault single-end ranging method
CN102967799A (en) * 2012-11-29 2013-03-13 深圳市双合电气股份有限公司 Comprehensive fault distance measuring method for electric power system
CN103116113A (en) * 2013-01-23 2013-05-22 深圳先进技术研究院 Single-ended fault distance detection method of power transmission line based on resonant frequency
CN103163427A (en) * 2013-03-07 2013-06-19 福建省电力有限公司 Method for realizing line single-phase earth fault single-terminal fault locating by using real part of voltage drop along line
CN103163424A (en) * 2013-01-31 2013-06-19 福建省电力有限公司 Line-to-ground fault distance measurement method based on instant sampling value
CN103207353A (en) * 2013-03-13 2013-07-17 福建省电力有限公司 Line inter-phase fault distance measuring method based on discrete sampling value
CN103217628A (en) * 2013-04-15 2013-07-24 国家电网公司 Electric transmission line phase fault single-ended distance measuring method
CN103248026A (en) * 2013-05-10 2013-08-14 国家电网公司 Line single-phase ground fault relay protection method capable of preventing distributed capacitive current and transitional resistance
CN103245876A (en) * 2013-04-15 2013-08-14 国家电网公司 Method for realizing line interphase fault single-ended distance measurement by means of real part characteristics of voltage drop
CN103245878A (en) * 2013-04-15 2013-08-14 国家电网公司 Single-end distance measurement method for electric transmission line single-phase earth fault
CN103245877A (en) * 2013-04-15 2013-08-14 国家电网公司 Method for ranging single-phase earth fault of line by using single-end electric quantity
CN103245890A (en) * 2013-05-10 2013-08-14 国家电网公司 Line single-phase ground fault single-terminal location method capable of preventing influences of both transitional resistance and load current
CN103293439A (en) * 2013-05-19 2013-09-11 国家电网公司 Single-terminal distance measurement method for line single-phase earth faults based on distributed parameter measurement impedance amplitude characteristics
CN103630798A (en) * 2013-09-12 2014-03-12 清华大学 Transmission line single-phase grounding fault single-end ranging method
CN103762571A (en) * 2014-02-18 2014-04-30 国家电网公司 Method for achieving circuit single-phase earth fault relay protection with hyperbolic tangent function amplitude characteristics
CN103779844A (en) * 2014-02-11 2014-05-07 华北电力大学 Self-adaptive distance protection system based on virtual voltage drop and protection method thereof
CN103891077A (en) * 2011-11-01 2014-06-25 Abb技术有限公司 Travelling-wave based fault protection of high-voltage transmission lines
CN104062553A (en) * 2014-07-09 2014-09-24 国家电网公司 Same-tower-double-circuit double-circuit line single-phase earth fault single-end distance measurement method
CN104237742A (en) * 2014-10-10 2014-12-24 华北电力大学 Method for improving fault location tolerance of power distribution network by utilizing fault current relative errors
CN104316842A (en) * 2014-11-14 2015-01-28 国家电网公司 Line phase fault single-ended distance measurement method by means of phase fault position factor phase characteristic
CN105044426A (en) * 2014-04-28 2015-11-11 洛克威尔自动控制技术股份有限公司 Calculating line-to-neutral voltages without a connection to a system neutral or earth ground
CN105606959A (en) * 2016-01-08 2016-05-25 清华大学 Single-end ranging method for arc light high-resistance earth fault of power transmission line
CN105866622A (en) * 2016-03-30 2016-08-17 国网福建省电力有限公司 Method for realizing power transmission line two-end fault distance measurement by using zero-sequence component
CN104466928B (en) * 2015-01-06 2017-11-28 华北电力大学 Distance protecting method and system
CN108535597A (en) * 2018-04-13 2018-09-14 国网福建省电力有限公司 Singlephase earth fault Section Location based on circuit model
CN109100605A (en) * 2018-10-23 2018-12-28 国网江苏省电力有限公司徐州供电分公司 Utilize the single end positioning method of the high-tension cable singlephase earth fault of failure boundary condition
CN110470950A (en) * 2019-08-16 2019-11-19 深圳供电局有限公司 Fault positioning method for transmission line
CN112578225A (en) * 2020-12-07 2021-03-30 广东电网有限责任公司韶关供电局 Method, device and system for determining single-phase earth fault point
CN113064022A (en) * 2021-03-12 2021-07-02 国网河南省电力公司电力科学研究院 A Line Protection Method Based on Transition Resistance Calculation
CN114371361A (en) * 2021-11-01 2022-04-19 天津大学 Single-ended fault location method of transmission line based on multivariate model fusion analysis

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5726574A (en) * 1996-03-11 1998-03-10 Electric Power Research Institute, Inc Method of locating a fault in an electric power cable
CN1508554A (en) * 2001-06-01 2004-06-30 清华大学 Fault line selection method for single-phase grounding in small current grounding system
CN1167176C (en) * 2001-08-17 2004-09-15 清华大学 Distribution Network Relay Protection and Fault Location System
DE10234537B4 (en) * 2002-07-30 2005-05-19 Hagenuk KMT Kabelmeßtechnik GmbH Method and device for locating cable sleeves and cable faults in laid cables
CN1928574A (en) * 2005-09-09 2007-03-14 北京富瑞菲格电力科技有限公司 Traveling wave accidents distance measuring device for hour and minute composite sampling electric transmission line
CN1793995B (en) * 2006-03-09 2010-11-17 保定浪拜迪电气股份有限公司 Measuring method of power transmission line failure distance

Cited By (60)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101267107B (en) * 2008-01-08 2010-04-21 江苏省电力试验研究院有限公司 The Method of Suppressing the Single-phase Short-circuit Current of Power Network
CN102023275A (en) * 2010-09-27 2011-04-20 华北电力大学 Single-end ranging method of line based on phase mutation characteristics of positioning function
CN102147443A (en) * 2011-01-13 2011-08-10 国网电力科学研究院 Single-end distance measuring method based on self-adaptive current
CN102129011A (en) * 2011-01-20 2011-07-20 福建省电力有限公司福州超高压输变电局 Single-ended phase-to-phase fault location method for distributed capacitance current and fault resistance resistant line
CN102129011B (en) * 2011-01-20 2012-12-05 福建省电力有限公司福州超高压输变电局 Single-ended interphase fault location method for distributed capacitance current and fault resistance resistant line
CN102129014A (en) * 2011-01-28 2011-07-20 福建省电力有限公司福州超高压输变电局 Method for realizing line phase-to-phase fault single-end distance measurement by utilizing distribution parameter model
CN102135589B (en) * 2011-02-18 2013-05-01 华北电力大学 Distance measurement method for line phase-to-phase fault distance protection
CN102135589A (en) * 2011-02-18 2011-07-27 华北电力大学 Distance measurement method for line phase-to-phase fault distance protection
CN102175954B (en) * 2011-03-09 2012-12-05 福建省电力有限公司福州超高压输变电局 Circuit inter-phase fault single-end ranging method
CN102175954A (en) * 2011-03-09 2011-09-07 福建省电力有限公司福州超高压输变电局 Circuit inter-phase fault single-end ranging method
CN103891077A (en) * 2011-11-01 2014-06-25 Abb技术有限公司 Travelling-wave based fault protection of high-voltage transmission lines
CN102967799A (en) * 2012-11-29 2013-03-13 深圳市双合电气股份有限公司 Comprehensive fault distance measuring method for electric power system
CN102967799B (en) * 2012-11-29 2015-02-25 深圳市双合电气股份有限公司 Comprehensive fault distance measuring method for electric power system
CN103116113A (en) * 2013-01-23 2013-05-22 深圳先进技术研究院 Single-ended fault distance detection method of power transmission line based on resonant frequency
CN103163424A (en) * 2013-01-31 2013-06-19 福建省电力有限公司 Line-to-ground fault distance measurement method based on instant sampling value
CN103163427A (en) * 2013-03-07 2013-06-19 福建省电力有限公司 Method for realizing line single-phase earth fault single-terminal fault locating by using real part of voltage drop along line
CN103163427B (en) * 2013-03-07 2015-07-01 福建省电力有限公司 Method for realizing line single-phase earth fault single-terminal fault locating by using real part of voltage drop along line
CN103207353A (en) * 2013-03-13 2013-07-17 福建省电力有限公司 Line inter-phase fault distance measuring method based on discrete sampling value
CN103207353B (en) * 2013-03-13 2015-08-26 福建省电力有限公司 Based on the line inter-phase fault distance measurement method of dis-crete sample values
CN103217628A (en) * 2013-04-15 2013-07-24 国家电网公司 Electric transmission line phase fault single-ended distance measuring method
CN103245878A (en) * 2013-04-15 2013-08-14 国家电网公司 Single-end distance measurement method for electric transmission line single-phase earth fault
CN103245876B (en) * 2013-04-15 2015-08-26 国家电网公司 Voltage drop real part characteristic is utilized to realize circuit inter-phase fault single-end ranging
CN103245876A (en) * 2013-04-15 2013-08-14 国家电网公司 Method for realizing line interphase fault single-ended distance measurement by means of real part characteristics of voltage drop
CN103245877B (en) * 2013-04-15 2015-07-01 国家电网公司 Method for ranging single-phase earth fault of line by using single-end electric quantity
CN103217628B (en) * 2013-04-15 2015-12-02 国家电网公司 A kind of electric transmission line phase fault method of single end distance measurement
CN103245877A (en) * 2013-04-15 2013-08-14 国家电网公司 Method for ranging single-phase earth fault of line by using single-end electric quantity
CN103245878B (en) * 2013-04-15 2015-08-26 国家电网公司 A kind of transmission line one-phase earth fault method of single end distance measurement
CN103248026B (en) * 2013-05-10 2015-11-04 国家电网公司 Relay protection method for line single-phase ground fault against distributed capacitive current and transition resistance
CN103248026A (en) * 2013-05-10 2013-08-14 国家电网公司 Line single-phase ground fault relay protection method capable of preventing distributed capacitive current and transitional resistance
CN103245890B (en) * 2013-05-10 2015-11-04 国家电网公司 Single-end distance measurement method for line single-phase ground faults against the influence of transition resistance and load current
CN103245890A (en) * 2013-05-10 2013-08-14 国家电网公司 Line single-phase ground fault single-terminal location method capable of preventing influences of both transitional resistance and load current
CN103293439B (en) * 2013-05-19 2015-12-02 国家电网公司 Based on distribution parameter measurement impedance magnitude characteristic line single-phase earth fault single-terminal location method
CN103293439A (en) * 2013-05-19 2013-09-11 国家电网公司 Single-terminal distance measurement method for line single-phase earth faults based on distributed parameter measurement impedance amplitude characteristics
CN103630798B (en) * 2013-09-12 2016-02-24 清华大学 Transmission line one-phase earth fault method of single end distance measurement
CN103630798A (en) * 2013-09-12 2014-03-12 清华大学 Transmission line single-phase grounding fault single-end ranging method
CN103779844A (en) * 2014-02-11 2014-05-07 华北电力大学 Self-adaptive distance protection system based on virtual voltage drop and protection method thereof
CN103779844B (en) * 2014-02-11 2016-05-18 华北电力大学 Based on the adaptive distance protection method of virtual voltage landing
CN103762571B (en) * 2014-02-18 2016-08-17 国家电网公司 Hyperbolic tangent function amplitude characteristic is utilized to realize single-phase line earth fault relay protection method
CN103762571A (en) * 2014-02-18 2014-04-30 国家电网公司 Method for achieving circuit single-phase earth fault relay protection with hyperbolic tangent function amplitude characteristics
CN105044426A (en) * 2014-04-28 2015-11-11 洛克威尔自动控制技术股份有限公司 Calculating line-to-neutral voltages without a connection to a system neutral or earth ground
CN105044426B (en) * 2014-04-28 2020-01-31 罗克韦尔自动化技术公司 Calculating line-to-neutral voltage without connecting system neutral or ground
CN104062553A (en) * 2014-07-09 2014-09-24 国家电网公司 Same-tower-double-circuit double-circuit line single-phase earth fault single-end distance measurement method
CN104062553B (en) * 2014-07-09 2016-07-06 国家电网公司 Double-circuit lines on the same pole road singlephase earth fault method of single end distance measurement
CN104237742A (en) * 2014-10-10 2014-12-24 华北电力大学 Method for improving fault location tolerance of power distribution network by utilizing fault current relative errors
CN104316842B (en) * 2014-11-14 2017-04-26 国家电网公司 Line phase fault single-ended distance measurement method by means of phase fault position factor phase characteristic
CN104316842A (en) * 2014-11-14 2015-01-28 国家电网公司 Line phase fault single-ended distance measurement method by means of phase fault position factor phase characteristic
CN104466928B (en) * 2015-01-06 2017-11-28 华北电力大学 Distance protecting method and system
CN105606959B (en) * 2016-01-08 2019-02-22 清华大学 Single-ended ranging method for arc high-resistance grounding faults in transmission lines
CN105606959A (en) * 2016-01-08 2016-05-25 清华大学 Single-end ranging method for arc light high-resistance earth fault of power transmission line
CN105866622A (en) * 2016-03-30 2016-08-17 国网福建省电力有限公司 Method for realizing power transmission line two-end fault distance measurement by using zero-sequence component
CN105866622B (en) * 2016-03-30 2018-10-16 国网福建省电力有限公司 Two ends of electric transmission line fault distance-finding method is realized using zero-sequence component
CN108535597A (en) * 2018-04-13 2018-09-14 国网福建省电力有限公司 Singlephase earth fault Section Location based on circuit model
CN109100605A (en) * 2018-10-23 2018-12-28 国网江苏省电力有限公司徐州供电分公司 Utilize the single end positioning method of the high-tension cable singlephase earth fault of failure boundary condition
CN110470950A (en) * 2019-08-16 2019-11-19 深圳供电局有限公司 Fault positioning method for transmission line
CN110470950B (en) * 2019-08-16 2021-07-16 深圳供电局有限公司 Transmission line fault location method
CN112578225A (en) * 2020-12-07 2021-03-30 广东电网有限责任公司韶关供电局 Method, device and system for determining single-phase earth fault point
CN113064022A (en) * 2021-03-12 2021-07-02 国网河南省电力公司电力科学研究院 A Line Protection Method Based on Transition Resistance Calculation
CN113064022B (en) * 2021-03-12 2022-04-29 国网河南省电力公司电力科学研究院 A Line Protection Method Based on Transition Resistance Calculation
CN114371361A (en) * 2021-11-01 2022-04-19 天津大学 Single-ended fault location method of transmission line based on multivariate model fusion analysis
CN114371361B (en) * 2021-11-01 2025-01-03 天津大学 A single-end fault location method for transmission lines based on multivariate model fusion analysis

Also Published As

Publication number Publication date
CN100526895C (en) 2009-08-12

Similar Documents

Publication Publication Date Title
CN101067641A (en) 2007-11-07 Single-ended fault location method of line against the influence of distributed capacitive current and transition resistance
CN102175954B (en) 2012-12-05 Circuit inter-phase fault single-end ranging method
CN101074979A (en) 2007-11-21 Method for measuring power-transmission circuit double-ended distance by distributing parameter
CN106093708B (en) 2019-01-04 A kind of modified common-tower double-return direct current transmission line double-end fault distance-finding method of velocity of wave
CN103401238B (en) 2016-01-20 A kind of power load modelling approach based on Measurement-based approach
CN110927510B (en) 2021-10-01 A frequency domain method for double-ended traveling wave fault location of transmission lines
CN101350521A (en) 2009-01-21 Fault distance-finding method for nonuniform zero sequence mutual inductance same-lever aerial multi-back line
CN103792466B (en) 2016-06-08 Overhead transmission line-power cable mixed connection line double-end Fault Locating Method
CN102081132A (en) 2011-06-01 Two-end distance measuring method of transmission line fault under dynamic condition
CN101635457A (en) 2010-01-27 Electric network parameter estimation method based on parameter sensitivity of state estimation residual error
CN102023275A (en) 2011-04-20 Single-end ranging method of line based on phase mutation characteristics of positioning function
CN105606959A (en) 2016-05-25 Single-end ranging method for arc light high-resistance earth fault of power transmission line
CN104914351A (en) 2015-09-16 Area power network fault positioning method based on optimal wave velocity
CN102129014A (en) 2011-07-20 Method for realizing line phase-to-phase fault single-end distance measurement by utilizing distribution parameter model
CN104062550A (en) 2014-09-24 Method for locating non-homonymic phase overline earth fault of double-circuit line based on binary search
CN103176102B (en) 2015-09-02 A kind of range finding yardstick minimum principle that utilizes realizes line single-phase earth fault single-terminal location method
CN103217623B (en) 2015-08-05 The line double-end fault distance-finding method had nothing to do with transition resistance and load current
CN107167703A (en) 2017-09-15 A kind of wind power plant current collection line single phase grounding failure distance-finding method
CN102298105A (en) 2011-12-28 Method for checking lightning shielding failure performance of line
CN103293441B (en) 2016-06-01 Distribution parameter is utilized to realize line single-phase earth fault single-terminal location method
CN105699858B (en) 2018-08-24 A kind of aerial-cable hybrid line fault distance-finding method considering connection resistance
CN103675565A (en) 2014-03-26 Method for measuring power frequency parameters of superhigh/extrahigh-voltage alternating-current (direct-current) power transmission circuit
CN105588984A (en) 2016-05-18 Mixed-pressure bipolar direct-current power transmission line zero-sequence parameter precisely measuring method
CN102147443B (en) 2013-07-17 Single-ended ranging method based on adaptive current
CN107154627B (en) 2019-09-24 A kind of small power station's group of planes Dynamic Equivalence suitable for bulk power grid analysis

Legal Events

Date Code Title Description
2007-11-07 C06 Publication
2007-11-07 PB01 Publication
2008-01-02 C10 Entry into substantive examination
2008-01-02 SE01 Entry into force of request for substantive examination
2009-08-12 C14 Grant of patent or utility model
2009-08-12 GR01 Patent grant