CN107505536A - The grounding net of transformer substation surface potential distribution calculation method of meter and more metallic conduits - Google Patents

The grounding net of transformer substation surface potential distribution calculation method of meter and more metallic conduits Download PDF

Info

Publication number
CN107505536A
CN107505536A CN201710665777.3A CN201710665777A CN107505536A CN 107505536 A CN107505536 A CN 107505536A CN 201710665777 A CN201710665777 A CN 201710665777A CN 107505536 A CN107505536 A CN 107505536A
Authority
CN
China
Prior art keywords
mtd
mrow
msub
mtr
mtable
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
CN201710665777.3A
Other languages
Chinese (zh)
Other versions
CN107505536B (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.)
Shanghai Jiaotong University
State Grid Shanghai Electric Power Co Ltd
Original Assignee
Shanghai Jiaotong University
State Grid Shanghai 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 Shanghai Jiaotong University, State Grid Shanghai Electric Power Co Ltd filed Critical Shanghai Jiaotong University
Priority to CN201710665777.3A priority Critical patent/CN107505536B/en
Publication of CN107505536A publication Critical patent/CN107505536A/en
Application granted granted Critical
Publication of CN107505536B publication Critical patent/CN107505536B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08—Locating faults in cables, transmission lines, or networks
    • G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
    • G01R31/086—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Prevention Of Electric Corrosion (AREA)

Abstract

The present invention relates to a kind of meter and the grounding net of transformer substation surface potential distribution calculation method of more metallic conduits, this method is on the basis of considering to interact between more metallic conduit itself longitudinal electrical resistances and pipeline stray electrical current and grounded screen diffusing, surface potential distribution to grounded screen carries out calculating analysis, to obtain grounding net of transformer substation surface potential distribution comprehensively and accurately result.Compared with prior art, the present invention has the advantages that to improve accuracy, the reduction amount of calculation that the distribution of grounding net of transformer substation surface potential calculates.

Description

The grounding net of transformer substation surface potential distribution calculation method of meter and more metallic conduits
Technical field
The present invention relates to grounding net of transformer substation technical field, more particularly, to the transformer station of a kind of meter and more metallic conduits Grounded screen surface potential distribution calculation method.
Background technology
Grounded screen is the important component of transformer station, and its reliability has great to the safe and stable operation of power system Meaning.For ensure grounded screen safe and reliable operation, how accurate and effective find grounded screen incipient fault so that targetedly Safeguard procedures are taken to turn into distinct issues the most in power industry operation maintenance work in ground.Based on connecing for surface potential distribution Earth mat method for diagnosing faults is at present using one of more means.But substation field environment is complex, underground metalliferous The presence of pipeline can have an impact to the accuracy of grounded screen surface potential and its fault diagnosis result.Therefore, it is necessary to pass through The influence that analysis and research underground metal pipes are distributed to grounding net of transformer substation surface potential, so as to effectively examining for grounded screen fault It is disconnected that foundation is provided.
The research that domestic and foreign scholars are distributed to grounding net of transformer substation surface potential is broadly divided into two classes：Numerical computations and scene Test.Wherein, numerical computations mainly carry out calculating analysis, conventional method according to ground connection web frame to the distribution of its surface potential Such as finite difference calculus, FInite Element, Analogue charge method, boundary element method.How on-the-spot test is mainly studied to transformer station Actual grounded screen surface potential distribution carries out more accurate test.
But the calculating analysis of grounding net of transformer substation surface potential is counted according to ground grid design drawing mostly at present Calculate.In addition, when analyzing grounding net of transformer substation fieldtesting results, it is believed that test result and grounded screen state and soil Earth structure is relevant, transformer station's underground metal pipes is have ignored to the influence caused by surface potential, so as to cause resulting phase Close the actual state that test result is difficult to reflect the distribution of grounding net of transformer substation surface potential.
The content of the invention
It is an object of the present invention to overcome the above-mentioned drawbacks of the prior art and provide one kind meter and more metals The grounding net of transformer substation surface potential distribution calculation method of pipeline.
The purpose of the present invention can be achieved through the following technical solutions：
The grounding net of transformer substation surface potential distribution calculation method of a kind of meter and more metallic conduits, comprises the following steps：
1) grounded screen Injection Current i is obtainedset；
2) grounded screen and the electric potential relation expression formula of more pipelines are established：
U0=R00I0+RA0I1+RB0I2+…+RK0IK=R00I+RA0TAU1+RB0TBU2+…+RK0TKUK
UA=R0AI0+RAAI1+RBAI2+…+RKAIK=R0AI+RAATAU1+RBATBU2+…+RKATKUK
UB=R0BI0+RABI1+RBBI2+…+RKBIK=R0BI+RABTAU1+RBBTBU2+…+RKBTKUK
UK=R0KI0+RAKI1+RBKI2+…+RKKIK=R0KI+RAKTAU1+RBKTBU2+…+RKKTKUK
In formula：R00For the resistance coefficient matrix of grounded screen itself；U0,U1,…,UKFor grounded screen and Duo Gen metallic conduits Current potential；I0,I1,…,IKIt is distributed for grounded screen and the diffusing of Duo Gen metallic conduits；T1,T2,…,TKFor more metallic conduits etc. Imitate circuit matrix；RIJIt is grounded screen and Duo Gen metallic conduits to itself and the resistance coefficient matrix that affects one another, I, J=0, 1,…,K；K is the total radical of metallic conduit；
Being write the electric potential relation expression formula as matrix form is：
In formula：RT=[RA0TA RB0TB … RK0TK]；RP=[R0A R0B … R0K]； EmFor m rank unit matrix；M is the total subdivision hop count of all metallic conduits；I dissipates for grounded screen Flow distribution；IA=[I1 I2 … IK] be distributed for the stray electrical current of K root metallic conduits；UAT=[U1 UB … UK] it is K root metals Pipeline current potential；U0For grounded screen current potential.
3) electric potential relation matrix equation described in Antisymmetric Iterative Algorithm for Solving is applied, obtains I, IA、UATAnd U0；
4) current potential for obtaining any point above grounding net of transformer substation is calculated according to step 3), grounding net of transformer substation top is any The current potential of point is defined as grounded screen diffusing and metallic conduit stray electrical current in current potential sum caused by the point.
The equivalent circuit matrix of the metallic conduit is calculated using imfinitesimal method and obtained.
The equivalent circuit matrix expression of i-th metallic conduit is：
Wherein, i=1 ..., K, ZL(j-1)For the longitudinal electrical resistance of the jth segmentation of i-th metallic conduit, j=1,2 ..., a-1； A is the segments of i-th metallic conduit.
In the step 3), using described in Antisymmetric Iterative Algorithm for Solving the step of electric potential relation matrix equation it is as follows：
3a) give an arbitrary matrix X1, its line number and columns are equal to matrix [I IA U0 UAT]TColumns, make i =1, the transposition of subscript T representing matrixs；
3b) W is calculated respectively according to following equatioies1、P1And Q1, calculation formula is
If 3c) W1=0, or W1≠ 0 and Q1=0, then stop calculating, otherwise i=i+1, goes to step 3d)；
3d) W is calculated respectively according to following equatioiesi+1、Pi+1And Qi+1, calculation formula is
In formula, | | Wk| | representing matrix WkNorm；The inner product of ＜ ＞ representing matrixs；
If 3e) Wi+1=0, or W1≠ 0 and Q1=0, then stop calculating, be transferred to step 3f), otherwise go to step 3d)；
3f) make Xk+1=[I IA U0 UAT]T, you can obtain I, IA、UATAnd U0。
In the step 2), during opening relationships matrix equation, equipotential volume is regarded as with grounded screen.
Total radical K >=1 of metallic conduit.
During the total radical K=0 of the metallic conduit, the relational matrix equation of foundation replaces with
In formula, u0For grounded screen current potential；
Solution above-mentioned relation matrix equation obtains the grounded screen current potential and diffusing distribution vector when being free of metallic conduit.
Compared with prior art, the present invention considers influence of transformer station's underground metal pipes to grounded screen, can be as far as possible Truly reflect grounding net of transformer substation surface potential distribution situation, can be that the grounded screen fault based on surface potential distribution is examined Disconnected method provides important evidence, has the advantages that：
(1) distribution of grounding net of transformer substation surface potential can more accurately and comprehensively be calculated.
(2) consider to interact between metallic conduit itself longitudinal electrical resistance and pipeline stray electrical current and grounded screen diffusing, improve The accuracy that the distribution of grounding net of transformer substation surface potential calculates.
(3) the equivalent circuit matrix of metallic conduit is obtained using imfinitesimal method, improves grounding net of transformer substation surface potential point The accuracy that cloth calculates.
(4) grounded screen is regarded as equipotential volume, and applies Antisymmetric Iterative Algorithm for Solving electric potential relation matrix equation, reduced Amount of calculation.
(5) it is adapted to grounded screen and nearby more metallic conduits is present, improves grounding net of transformer substation earth's surface Calculation of Potential Distribution it is comprehensive.
Brief description of the drawings
Fig. 1 is the schematic flow sheet of the present invention；
The surface potential distribution for the grounded screen that Fig. 2 is calculated when being in the absence of metallic conduit；
Fig. 3 is the surface potential distribution for the grounded screen being calculated when metallic conduit be present；
Fig. 4 is the grounded screen surface potential distribution that measurement obtains；
Fig. 5 is surface potential calculating and the measurement result for the conductor that metallic conduit passes through.
Embodiment
The present invention is described in detail with specific embodiment below in conjunction with the accompanying drawings.The present embodiment is with technical solution of the present invention Premised on implemented, give detailed embodiment and specific operating process, but protection scope of the present invention is not limited to Following embodiments.
As shown in figure 1, the present embodiment provides the grounding net of transformer substation surface potential distribution of a kind of meter and more metallic conduits Computational methods, comprise the following steps：
1) grounded screen Injection Current i is obtainedset；
2) equipotential volume is regarded as with grounded screen, establishes grounded screen and the electric potential relation expression formula of more pipelines：
U0=R00I0+RA0I1+RB0I2+…+RK0IK=R00I+RA0TAU1+RB0TBU2+…+RK0TKUK
UA=R0AI0+RAAI1+RBAI2+…+RKAIK=R0AI+RAATAU1+RBATBU2+…+RKATKUK
UB=R0BI0+RABI1+RBBI2+…+RKBIK=R0BI+RABTAU1+RBBTBU2+…+RKBTKUK
UK=R0KI0+RAKI1+RBKI2+…+RKKIK=R0KI+RAKTAU1+RBKTBU2+…+RKKTKUK
In formula：R00For the resistance coefficient matrix of grounded screen itself；U0,U1,…,UKFor grounded screen and Duo Gen metallic conduits Current potential；I0,I1,…,IKIt is distributed for grounded screen and the diffusing of Duo Gen metallic conduits；T1,T2,…,TKFor more metallic conduits etc. Imitate circuit matrix；RIJ(I, J=0,1 ..., K) is grounded screen and Duo Gen metallic conduits to itself and the resistance coefficient that affects one another Matrix；K is the total radical of metallic conduit.
Being write the electric potential relation expression formula as matrix form is：
In formula：RT=[RA0TA RB0TB … RK0TK]；RP=[R0A R0B … R0K]； EmFor m rank unit matrix；M is the total subdivision hop count of all metallic conduits；I dissipates for grounded screen Flow distribution；IA=[I1 I2 … IK] be distributed for the stray electrical current of K root metallic conduits；UAT=[U1 UB … UK] it is metallic conduit Potential distribution；U0For grounded screen current potential.
3) electric potential relation matrix equation described in Antisymmetric Iterative Algorithm for Solving is applied, obtains I, IA、UATAnd U0.Basic step It is as follows：
3a) give an arbitrary matrix X1, its line number and columns are equal to matrix [I IA U0 UAT]TColumns, i= 1.Herein, the transposition of subscript T representing matrixs；
3b) W is calculated respectively according to following equatioies1、P1And Q1, calculation formula is
If 3c) W1=0, or W1≠ 0, Q1=0, then stop calculating, otherwise i=i+1, turns 3d)；
3d) W is calculated respectively according to following equatioiesi+1、Pi+1And Qi+1, calculation formula is
In formula, | | Wk| | representing matrix WkNorm；The inner product of ＜ ＞ representing matrixs.
If 3e) Wi+1=0, or W1≠ 0, Q1=0, then stop calculating, be transferred to 3f), otherwise turn 3d)；
3f) make Xk+1=[I IA U0 UAT]T, you can obtain I, IA、UATAnd U0。
4) current potential for obtaining space any point is calculated according to step 3), the current potential of the space any point is defined as grounded screen and dissipated Stream and metallic conduit stray electrical current are in current potential sum caused by the point.
The equivalent circuit matrix of the metallic conduit is calculated using imfinitesimal method and obtained, the equivalent circuit of i-th metallic conduit Matrix expression is：
Wherein, i=1 ..., K, ZL(j-1)The longitudinal electrical resistance being segmented for the jth of i-th metallic conduit, j=1,2 ..., a-1, A is the segments of i-th metallic conduit.
Above-mentioned distribution calculation method is applied to the situation of total radical K >=1 of metallic conduit.
As the total radical K=0 of metallic conduit, the relational matrix equation of foundation replaces with
In formula, u0For grounding grids current potential；
Solution above-mentioned relation matrix equation obtains the grounded screen current potential and diffusing distribution vector when being free of metallic conduit.
Reference picture 1, to verify that the grounding net of transformer substation surface potential proposed by the invention based on field circuit method is distributed meter The accuracy of calculation method, live surface potential distribution measuring is carried out to certain actual 35kV grounding net of transformer substation.
Fig. 2-Fig. 3 show the surface potential distribution for the 35kV grounding net of transformer substation being calculated, and Fig. 2 is actual conditions Under (i.e. there are underground metal pipes in transformer station) surface potential distribution, Fig. 3 is to assume when transformer station is not present metallic conduit Surface potential.Herein, 2 duct lengths are 2m, and pipe radius is 4cm.From Fig. 2-Fig. 3, when transformer station is not present ground During lower metallic conduit, the surface potential distribution curve of grounded screen is more smooth, and the surface potential above conductor is apparently higher than mesh Top, current potential also can be more slightly higher than conductor interlude at the end node of conductor two, meanwhile, it is higher closer to Current injection points current potential；Ground Table current potential maximum is 1.513V, appears in Current injection points；Potential minimum is 0.693V, is appeared in the mesh of lower left The heart.When transformer station has metallic conduit, substantially rule, surface potential maximum and the minimum value of surface potential distribution be not with depositing Situation in pipeline is consistent；But the surface potential curve above pipeline is not very smooth.
Fig. 4 show the grounded screen surface potential distribution that measurement obtains.As seen from the figure, the earth's surface electricity obtained by in-site measurement Bit distribution rule is consistent with result of calculation.Surface potential maximum is equally present in Current injection points, size 1.484V, with meter The error for calculating result is only 1.94%；Surface potential minimum value is equally present in lower left mesh center, size 0.710V, with The error of result of calculation is only 2.50%.Surface potential distribution measuring and the consistency checking of result of calculation are of the invention to be succeeded in one's scheme Calculate the accuracy of result.
In order to more intuitively embody the grounding net of transformer substation surface potential numerical value proposed by the invention based on field circuit method The accuracy of computational methods, the surface potential of OB sections conductor (conductor that i.e. metallic conduit passes through) is calculated Fig. 5 and measurement result Compare.As seen from the figure, the substantially changing rule of 3 surface potential curves is consistent, is simply passed through in metallic conduit Region certain difference be present.Compare 2 result of calculation curves it can be found that there is earth's surface electricity in the region of metallic conduit situation About 6.66% is at most have dropped when position is than without pipeline.In the region, measurement result has higher overlap with corresponding result of calculation Degree, only because the factor such as measurement error and scene interference influences, part measurement result has little bit different with result of calculation, but poor DRS degree is smaller, in tolerance interval.For surface potential measurement result compared with the result of calculation without pipe condition, current potential is most It has dropped about 8.29%.Therefore, the grounding net of transformer substation surface potential numerical computations side proposed by the present invention based on field circuit method Method can more accurately and comprehensively calculate the distribution of grounding net of transformer substation surface potential.Meanwhile consider transformer station's underground metal pipes Influence to grounded screen, can truly reflect grounding net of transformer substation surface potential distribution situation as far as possible, can be based on The Fault Diagnosis for Grounding Grids method of surface potential distribution provides important evidence.
Preferred embodiment of the invention described in detail above.It should be appreciated that one of ordinary skill in the art without Creative work can is needed to make many modifications and variations according to the design of the present invention.Therefore, all technologies in the art Personnel are available by logical analysis, reasoning, or a limited experiment on the basis of existing technology under this invention's idea Technical scheme, all should be in the protection domain being defined in the patent claims.

Claims (7)

  1. A kind of 1. grounding net of transformer substation surface potential distribution calculation method of meter and more metallic conduits, it is characterised in that including Following steps：
    1) grounded screen Injection Current i is obtainedset；
    2) grounded screen and the electric potential relation expression formula of more pipelines are established：
    U0=R00I0+RA0I1+RB0I2+…+RK0IK=R00I+RA0TAU1+RB0TBU2+…+RK0TKUK
    UA=R0AI0+RAAI1+RBAI2+…+RKAIK=R0AI+RAATAU1+RBATBU2+…+RKATKUK
    UB=R0BI0+RABI1+RBBI2+…+RKBIK=R0BI+RABTAU1+RBBTBU2+…+RKBTKUK
    UK=R0KI0+RAKI1+RBKI2+…+RKKIK=R0KI+RAKTAU1+RBKTBU2+…+RKKTKUK
    In formula：R00For the resistance coefficient matrix of grounded screen itself；U0,U1,…,UKFor grounded screen and the current potential of Duo Gen metallic conduits； I0,I1,…,IKIt is distributed for grounded screen and the diffusing of Duo Gen metallic conduits；T1,T2,…,TKFor the equivalent circuit of more metallic conduits Matrix；RIJIt is grounded screen and Duo Gen metallic conduits to itself and the resistance coefficient matrix that affects one another, I, J=0,1 ..., K；K is The total radical of metallic conduit；
    Being write the electric potential relation expression formula as matrix form is：
    <mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>R</mi> <mn>00</mn> </msub> </mtd> <mtd> <msub> <mi>R</mi> <mi>T</mi> </msub> </mtd> <mtd> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>R</mi> <mi>P</mi> </msub> </mtd> <mtd> <mrow> <mi>G</mi> <mi>T</mi> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>E</mi> <mi>m</mi> </msub> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>G</mi> <mi>T</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>I</mi> </mtd> </mtr> <mtr> <mtd> <msub> <mi>I</mi> <mi>A</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>U</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>U</mi> <mrow> <mi>A</mi> <mi>T</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mrow> <mi>s</mi> <mi>e</mi> <mi>t</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow>
    In formula：RT=[RA0TA RB0TB ... RK0TK]；RP=[R0A R0B … R0K]； EmFor m rank unit matrix；M is the total subdivision hop count of all metallic conduits；I dissipates for grounded screen Flow distribution；IA=[I1 I2 … IK] be distributed for the stray electrical current of K root metallic conduits；UAT=[U1 UB … UK] it is K root metals Pipeline current potential；U0For grounded screen current potential.
    3) electric potential relation matrix equation described in Antisymmetric Iterative Algorithm for Solving is applied, obtains I, IA、UATAnd U0；
    4) current potential for obtaining any point above grounding net of transformer substation is calculated according to step 3), grounding net of transformer substation top any point Current potential is defined as grounded screen diffusing and metallic conduit stray electrical current in current potential sum caused by the point.
  2. 2. the grounding net of transformer substation surface potential distribution calculation method of meter according to claim 1 and more metallic conduits, Obtained characterized in that, the equivalent circuit matrix of the metallic conduit is calculated using imfinitesimal method.
  3. 3. the grounding net of transformer substation surface potential distribution calculation method of meter according to claim 2 and more metallic conduits, Characterized in that, the equivalent circuit matrix expression of i-th metallic conduit is：
    Wherein, i=1 ..., K, ZL(j-1)For the longitudinal electrical resistance of the jth segmentation of i-th metallic conduit, j=1,2 ..., a-1；A is The segments of i-th metallic conduit.
  4. 4. the grounding net of transformer substation surface potential computational methods of meter according to claim 1 and more metallic conduits, it is special Sign is, in the step 3), using described in Antisymmetric Iterative Algorithm for Solving the step of electric potential relation matrix equation it is as follows：
    3a) give an arbitrary matrix X1, its line number and columns are equal to matrix [I IA U0 UAT]TColumns, make i=1, The transposition of subscript T representing matrixs；
    3b) W is calculated respectively according to following equatioies1、P1And Q1, calculation formula is
    <mrow> <msub> <mi>W</mi> <mn>1</mn> </msub> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mrow> <mi>s</mi> <mi>e</mi> <mi>t</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>R</mi> <mn>00</mn> </msub> </mtd> <mtd> <msub> <mi>R</mi> <mi>T</mi> </msub> </mtd> <mtd> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>R</mi> <mi>P</mi> </msub> </mtd> <mtd> <mrow> <mi>G</mi> <mi>T</mi> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>E</mi> <mi>m</mi> </msub> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>G</mi> <mi>T</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <msub> <mi>X</mi> <mn>1</mn> </msub> <mi>E</mi> </mrow>
    <mrow> <msub> <mi>P</mi> <mn>1</mn> </msub> <mo>=</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>R</mi> <mn>00</mn> </msub> </mtd> <mtd> <msub> <mi>R</mi> <mi>T</mi> </msub> </mtd> <mtd> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>R</mi> <mi>P</mi> </msub> </mtd> <mtd> <mrow> <mi>G</mi> <mi>T</mi> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>R</mi> <mi>m</mi> </msub> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>G</mi> <mi>T</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <msub> <mi>W</mi> <mn>1</mn> </msub> <mi>E</mi> </mrow>
    <mrow> <msub> <mi>Q</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mn>1</mn> </msub> <mo>-</mo> <msubsup> <mi>P</mi> <mn>1</mn> <mi>T</mi> </msubsup> <mo>)</mo> </mrow> </mrow>
    If 3c) W1=0, or W1≠ 0 and Q1=0, then stop calculating, otherwise i=i+1, goes to step 3d)；
    3d) W is calculated respectively according to following equatioiesi+1、Pi+1And Qi+1, calculation formula is
    <mrow> <msub> <mi>X</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <msub> <mi>X</mi> <mi>k</mi> </msub> <mo>+</mo> <mfrac> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>W</mi> <mi>k</mi> </msub> <mo>|</mo> <msup> <mo>|</mo> <mn>2</mn> </msup> </mrow> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>Q</mi> <mi>k</mi> </msub> <mo>|</mo> <msup> <mo>|</mo> <mn>2</mn> </msup> </mrow> </mfrac> </mrow>
    <mrow> <msub> <mi>W</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mrow> <mi>s</mi> <mi>e</mi> <mi>t</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>R</mi> <mn>00</mn> </msub> </mtd> <mtd> <msub> <mi>R</mi> <mi>T</mi> </msub> </mtd> <mtd> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>R</mi> <mi>P</mi> </msub> </mtd> <mtd> <mrow> <mi>G</mi> <mi>T</mi> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>E</mi> <mi>m</mi> </msub> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>G</mi> <mi>T</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <msub> <mi>X</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mi>E</mi> </mrow>
    <mrow> <msub> <mi>P</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>R</mi> <mn>00</mn> </msub> </mtd> <mtd> <msub> <mi>R</mi> <mi>T</mi> </msub> </mtd> <mtd> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>R</mi> <mi>P</mi> </msub> </mtd> <mtd> <mrow> <mi>G</mi> <mi>T</mi> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>E</mi> <mi>m</mi> </msub> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>G</mi> <mi>T</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <msub> <mi>W</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mi>E</mi> </mrow>
    <mrow> <msub> <mi>Q</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msubsup> <mi>P</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>T</mi> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <mo>&lt;</mo> <msub> <mi>P</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>&amp;CenterDot;</mo> <msub> <mi>Q</mi> <mi>i</mi> </msub> <mo>&gt;</mo> </mrow> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>Q</mi> <mi>i</mi> </msub> <mo>|</mo> <msup> <mo>|</mo> <mn>2</mn> </msup> </mrow> </mfrac> <msub> <mi>Q</mi> <mi>i</mi> </msub> </mrow>
    In formula, | | Wk| | representing matrix WkNorm；The inner product of ＜ ＞ representing matrixs；
    If 3e) Wi+1=0, or W1≠ 0 and Q1=0, then stop calculating, be transferred to step 3f), otherwise go to step 3d)；
    3f) make Xk+1=[I IA U0 UAT]T, you can obtain I, IA、UATAnd U0。
  5. 5. grounding net of transformer substation surface potential distribution calculation method according to claim 1, it is characterised in that the step 2) in, during opening relationships matrix equation, equipotential volume is regarded as with grounded screen.
  6. 6. the grounding net of transformer substation surface potential distribution calculation method of meter according to claim 1 and more metallic conduits, Characterized in that, total radical K >=1 of metallic conduit.
  7. 7. the grounding net of transformer substation surface potential distribution calculation method of meter according to claim 1 and more metallic conduits, Characterized in that, during the total radical K=0 of the metallic conduit, the relational matrix equation of foundation replaces with
    <mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>R</mi> <mn>00</mn> </msub> </mtd> <mtd> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>I</mi> </mtd> </mtr> <mtr> <mtd> <msub> <mi>u</mi> <mn>0</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mrow> <mi>s</mi> <mi>e</mi> <mi>t</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow>
    In formula, u0For grounded screen current potential；
    Solution above-mentioned relation matrix equation obtains the grounded screen current potential and diffusing distribution vector when being free of metallic conduit.
CN201710665777.3A 2017-08-07 2017-08-07 Transformer substation grounding grid earth surface potential distribution calculation method considering multiple metal pipelines Active CN107505536B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710665777.3A CN107505536B (en) 2017-08-07 2017-08-07 Transformer substation grounding grid earth surface potential distribution calculation method considering multiple metal pipelines

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710665777.3A CN107505536B (en) 2017-08-07 2017-08-07 Transformer substation grounding grid earth surface potential distribution calculation method considering multiple metal pipelines

Publications (2)

Publication Number Publication Date
CN107505536A true CN107505536A (en) 2017-12-22
CN107505536B CN107505536B (en) 2020-03-31

Family

ID=60688930

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710665777.3A Active CN107505536B (en) 2017-08-07 2017-08-07 Transformer substation grounding grid earth surface potential distribution calculation method considering multiple metal pipelines

Country Status (1)

Country Link
CN (1) CN107505536B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112632740A (en) * 2020-09-30 2021-04-09 南方电网科学研究院有限责任公司 Earth surface potential measuring method based on stray current model
CN112677833A (en) * 2021-01-18 2021-04-20 中车青岛四方机车车辆股份有限公司 Rail vehicle shielding grounding system, wiring method thereof and rail vehicle

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5596267A (en) * 1994-12-22 1997-01-21 Atlantic Richfield Company Apparatus and method for measuring potentials through pavements for buried pipeline cathodic protection systems
KR20100057972A (en) * 2008-11-24 2010-06-03 한빛이디에스(주) Diagnostic technique of mesh ground by analyzing the ground electric potential
CN102495327A (en) * 2011-12-13 2012-06-13 广东电网公司珠海供电局 Detection method and device for design of grounding grid of substation
CN104898024A (en) * 2015-06-10 2015-09-09 国网上海市电力公司 Transformer station grounding network fault diagnosis method based on synthesis of surface potential and magnetic induction strength
CN104931853A (en) * 2015-06-08 2015-09-23 上海交通大学 Diagnosis method for transformer-station grounding net based on surface potential distribution
CN104965134A (en) * 2015-06-10 2015-10-07 国网上海市电力公司 Echo state network-based transformer station grounding grid earth surface potential distribution calculation method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5596267A (en) * 1994-12-22 1997-01-21 Atlantic Richfield Company Apparatus and method for measuring potentials through pavements for buried pipeline cathodic protection systems
KR20100057972A (en) * 2008-11-24 2010-06-03 한빛이디에스(주) Diagnostic technique of mesh ground by analyzing the ground electric potential
CN102495327A (en) * 2011-12-13 2012-06-13 广东电网公司珠海供电局 Detection method and device for design of grounding grid of substation
CN104931853A (en) * 2015-06-08 2015-09-23 上海交通大学 Diagnosis method for transformer-station grounding net based on surface potential distribution
CN104898024A (en) * 2015-06-10 2015-09-09 国网上海市电力公司 Transformer station grounding network fault diagnosis method based on synthesis of surface potential and magnetic induction strength
CN104965134A (en) * 2015-06-10 2015-10-07 国网上海市电力公司 Echo state network-based transformer station grounding grid earth surface potential distribution calculation method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王丰华等: "采用地表电位和磁感应强度进行变电站接地网故障诊断的效果对比", 《高电压技术》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112632740A (en) * 2020-09-30 2021-04-09 南方电网科学研究院有限责任公司 Earth surface potential measuring method based on stray current model
CN112632740B (en) * 2020-09-30 2024-02-13 南方电网科学研究院有限责任公司 Surface potential measurement method based on stray current model
CN112677833A (en) * 2021-01-18 2021-04-20 中车青岛四方机车车辆股份有限公司 Rail vehicle shielding grounding system, wiring method thereof and rail vehicle

Also Published As

Publication number Publication date
CN107505536B (en) 2020-03-31

Similar Documents

Publication Publication Date Title
CN109256771A (en) Metro stray current and its caused transformer neutral point DC current calculation method
Zhang et al. Diagnosis of breaks in substation's grounding grid by using the electromagnetic method
CN108562782B (en) Stray current acquisition method and stray current detection system
CN105954594A (en) Novel grounding network grounding resistance inverse short-distance measurement method and apparatus
Qi et al. Calculation of overvoltage on nearby underground metal pipeline due to the lightning strike on UHV AC transmission line tower
CN104898024B (en) The Failure Diagnosis of Substation Ground Network method of comprehensive surface potential and magnetic induction intensity
CN103135137B (en) A kind of grounded screen topology detection method based on the differential method
CN108647438A (en) A kind of new soil equivalent resistance model modeling method
CN107142482B (en) The determination and means of defence that shape high voltage DC earthing pole interferes buried metal works
CN103954894A (en) Partial discharge locating method for three-phase crossed and interconnected cables
CN107505536A (en) The grounding net of transformer substation surface potential distribution calculation method of meter and more metallic conduits
CN105116230A (en) Method using multi-frequency combination to measure impulse grounding resistance of grounding device
CN104599193A (en) Rule base based single-phase ground fault positioning method of distribution network
CN104897996A (en) Transformer station grounding network corrosion diagnosis method based on magnetic field inverse problem solving
CN107451744A (en) A kind of substation grounding unit state appraisal procedure based on fuzzy overall evaluation
CN104965134B (en) Grounding net of transformer substation surface potential distribution calculation method based on echo state network
CN105510733A (en) Parameter measurement method for high-voltage and ultra-high-voltage direct-current transmission lines
Wang et al. An edge computing method using a novel mode component for power transmission line fault diagnosis in distribution network
Martins-Britto et al. Inductive interferences between a 500 kV power line and a pipeline with a complex approximation layout and multilayered soil
Liu et al. A magnetic detecting and evaluation method of substation’s grounding grids with break and corrosion
Chen et al. Travelling wave fault location approach for MMC-HVDC transmission line based on frequency modification algorithm
CN206618805U (en) A kind of super-huge Grounding impedance measurement device
Leal et al. A new approach for ground resistance measurements in onshore wind farms based on clamp-on meters and artificial neural network
CN111796165B (en) Power distribution network fault positioning method based on self-adaptive Fourier transform
CN104833883B (en) A kind of earth mat method of testing based on 10 35kV short circuit groundings

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant