CN106021652A - An earth soil three-dimensional resistance network model establishing method - Google Patents

An earth soil three-dimensional resistance network model establishing method Download PDF

Info

Publication number
CN106021652A
CN106021652A CN201610301434.4A CN201610301434A CN106021652A CN 106021652 A CN106021652 A CN 106021652A CN 201610301434 A CN201610301434 A CN 201610301434A CN 106021652 A CN106021652 A CN 106021652A
Authority
CN
China
Prior art keywords
soil
resistance
earth
grid
network model
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
CN201610301434.4A
Other languages
Chinese (zh)
Other versions
CN106021652B (en
Inventor
赵丹丹
傅晨钊
苏磊
司文荣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Shanghai Electric Power Co Ltd
East China Power Test and Research Institute Co Ltd
Original Assignee
State Grid Shanghai Electric Power Co Ltd
East China Power Test and Research Institute Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Shanghai Electric Power Co Ltd, East China Power Test and Research Institute Co Ltd filed Critical State Grid Shanghai Electric Power Co Ltd
Priority to CN201610301434.4A priority Critical patent/CN106021652B/en
Publication of CN106021652A publication Critical patent/CN106021652A/en
Application granted granted Critical
Publication of CN106021652B publication Critical patent/CN106021652B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • 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
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Economics (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Geometry (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention relates to an earth soil three-dimensional resistance network model establishing method. The method comprises the steps of simplifying earth soil into a plurality of soil grids; acquiring the geological conditions of the earth soil; establishing a finite element calculating model according to the geological conditions of the earth soil; calculating the electric current electric field distribution in the earth soil according to the finite element calculating model; according to the electric current electric field distribution, calculating the normal plane average potential Ur and the normal electric current vector sum (as shown in the specification) in the radial direction of the soil grids; according to the average potential and the electric current vector sum, calculating the radial equivalent resistance Rr, crosswise equivalent resistance R theta and longitudinal depth equivalent resistance Rd corresponding to the soil grids; abstracting the soil grids into a three-dimensional resistance network according to the radial equivalent resistance Rr, crosswise equivalent resistance R theta and longitudinal depth equivalent resistance Rd, and establishing a three-dimensional resistance network model. Compared with conventional methods, the method has the advantages of small calculating quantity and high calculating accuracy.

Description

A kind of method for building up of the earth soil solid resistance network model
Technical field
The present invention relates to D.C. magnetic biasing research and transformator running technology field, especially relate to a kind of the earth soil The method for building up of three-dimensional resistance network model.
Background technology
In the case of straight-flow system debugging, overhauling or break down, D.C. high voltage transmission can use one pole the earth to return The method of operation on road, huge DC current flows into the earth through direct current grounding pole, and causes ground electricity in a big way The significant change of position so that transformer station's earth potential offsets, and causes DC magnetic bias current.At present to D.C. magnetic biasing The computational methods of electric current mainly have two kinds: a kind of and management, on the Poisson's equation characterizing constant-current source field, utilizes business Industry software part CDEGS, the method amount of calculation is huge, and computational efficiency is low, the electrical network small scale that can consider；Additionally A kind of method is to set up underground, region resistor network, by the calculating that the computational short cut of Electromagnetic field is " road ", calculates Measure little, calculate speed fast, it is contemplated that electrical network scale big, but the earth soil equivalence of multilamellar is by the method at present Monolayer resistor network, does not accounts for the longitudinal propagation impact on result of calculation of soil layering and electric current, calculates Accuracy is the poorest.
Summary of the invention
It is an object of the invention to provide the foundation side of a kind of the earth soil solid resistance network model for the problems referred to above Method.
The purpose of the present invention can be achieved through the following technical solutions:
The method for building up of a kind of the earth soil solid resistance network model, the method comprises the following steps:
(1) the earth soil is simplified to multiple soil grid；
(2) geological conditions of the earth soil is obtained；
(3) limited element calculation model is built according to the geological conditions of the earth soil obtained in step (2)；
(4) current electric fields distribution in the earth soil is calculated according to the limited element calculation model obtained in step (3)；
(5) according to the current electric fields distribution obtained in step (4), calculate soil grid Normal plane radially and put down All current potential UrAnd normal direction current phasor and
(6) according to Normal plane average potential U of the soil grid radial direction obtained in step (5)rAnd normal direction Current phasor andCalculate the radial direction substitutional resistance R that soil grid is correspondingr, width is to substitutional resistance RθAnd depth To substitutional resistance Rd；
(7) according to the radial direction substitutional resistance R obtained in step (6)r, width is to substitutional resistance RθAnd depth to Substitutional resistance Rd, by abstract for the soil grid in step (1) be three-dimensional resistor network, set up three-dimensional resistor network Model.
The geological conditions of described the earth soil includes stratification state, landform and soil resistivity.
Described soil grid is fan-shaped.
In described step (5), the radial direction substitutional resistance R of i-th soil gridriResistance be:
R r i = U r i - U r ( i + 1 ) 1 2 ( I · r i + I · r ( i + 1 ) ) ;
Wherein, Uri、Ur(i+1)It is respectively i-th and i+1 soil grid Normal plane average potential radially；Be respectively i-th and i+1 soil grid normal direction current phasor radially with.
In described step (5), the width of i-th soil grid is to substitutional resistance RθiResistance be:
R θ i = 2 πθR r i r 1 i ( r 2 i - r 1 i ) * l n r 2 i r 1 i ;
Wherein, θ is the soil grid angle relative to the center of circle, r1i、r2iIt is respectively the interior outer radius of i-th soil grid.
In described step (5), the depth of i-th soil grid is to substitutional resistance RdiResistance be:
R d i = R r i d 2 ( r 2 i 2 - r 1 i 2 ) * l n r 2 i r 1 i ;
Wherein, d is the degree of depth of i-th soil grid, r1i、r2iIt is respectively the interior outer radius of i-th soil grid.
Compared with prior art, the present invention sets up a kind of energy consideration the earth soil layering and the multidirectional propagation of electric current The convenient computation model of DC current characteristics of the dispersed flow distribution, and then realize on a large scale, transformator under multi-voltage grade DC magnetic bias current simulation calculation, has the advantages that
1) present invention establish the earth soil solid resistance network model, in the calculating of DC magnetic bias current by field The calculating that computational short cut is " road " that road combines, amount of calculation is little, calculates speed fast；
2) in the earth soil solid resistance network model that the present invention sets up, consider respectively radially substitutional resistance, Width to substitutional resistance and depth to substitutional resistance, compared with prior art, it is contemplated that soil layering is vertical with electric current To the impact propagated result of calculation, accuracy in computation is high.
Accompanying drawing explanation
Fig. 1 is the diffusing current characteristics schematic diagram of soil grid；
Fig. 2 is three-dimensional resistance network；
Fig. 3 is the flow chart obtaining three-dimensional resistor network.
Detailed description of the invention
The present invention is described in detail with specific embodiment below in conjunction with the accompanying drawings.The present embodiment is with the technology of the present invention side Implement premised on case, give detailed embodiment and concrete operating process, but the protection model of the present invention Enclose and be not limited to following embodiment.
As it is shown in figure 1, when direct current grounding pole injection current, DC current, is formed at a distance by the earth diffusing DC electric field centered by direct current grounding pole.The geological conditions of the earth is different, directly affects DC current diffusing special Property.In some region of soil, there is radial component, tangential component and downward vertical component in the electric current of diffusing.
In ground, direct current is the permanent steady current field centered by direct current grounding pole in fact.In permanent steady current field, flow into arbitrarily The electric current at the interface comprising current source is equal to the total current of current source, i.e.
J=σ E
V=∫LEdL
It can be seen that the distribution of current density, J in soil and electric field intensity E and electromotive force V is all only and soil Resistivity.Therefore, soil is simplified to three-dimensional resistor network by the present invention, i.e. such as the soil grid in Fig. 1, abstract For such as the three-dimensional resistor network of Fig. 2.Three-dimensional resistor network is made up of the soil grid of multiple abstract mistakes, three-dimensional resistance Each angle of network includes the radial direction resistance R gone out by soil grid computingr, width is to resistance RθAnd depth is to electricity Resistance Rd。
Resistance value in three-dimensional resistor network can obtain through FEM calculation, is considered as actual soil-geological in calculating Condition, the information such as including layering, landform, soil resistivity.Calculation process as it is shown on figure 3, particularly as follows:
(1) the earth soil is simplified to multiple soil grid；
(2) geological conditions of the earth soil is obtained；
(3) limited element calculation model is built according to the geological conditions of the earth soil obtained in step (2)；
(4) current electric fields distribution in the earth soil is calculated according to the limited element calculation model obtained in step (3)；
(5) according to the current electric fields distribution obtained in step (4), i-th soil grid normal direction radially is calculated Plane average potential UriAnd normal direction current phasor and
(6) according to the average potential obtained in step (5) and current phasor and, calculate i-th soil grid pair The radial direction substitutional resistance R answeredri, width is to substitutional resistance RθiAnd depth is to substitutional resistance Rdi；
(7) according to the radial direction substitutional resistance R obtained in (6)ri, width is to substitutional resistance RθiAnd depth is to equivalence Resistance Rdi, by abstract for the soil grid in (1) be three-dimensional resistor network, set up three-dimensional resistance network model.
The radial direction tried to achieve in step (6), width to and depth be respectively to substitutional resistance
R r i = U r i - U r ( i + 1 ) 1 2 ( I · r i + I · r ( i + 1 ) ) R θ i = 2 πθR r i r 1 i ( r 2 i - r 1 i ) * ln r 2 i r 1 i R d i = R r i d 2 ( r 2 i 2 - r 1 i 2 ) * ln r 2 i r 1 i
In formula, Uri、Ur(i+1)It is respectively the average potential on the face vertical with radial direction of i-th soil grid；
Be respectively i-th soil grid the face vertical with radial direction on normal direction current phasor and；
r1i、r2iIt is respectively the interior outer radius of fan-shaped soil grid；
D is the degree of depth of this sector soil grid.
Greatly the three-dimensional resistance network model of soil can not only simulate the characteristics of the dispersed flow of stratified soil, it is possible to simulates into ground The polytropism of DC current, is applied to this resistance network model in the calculating of D.C. magnetic biasing, is greatly improved calculating Accuracy.

Claims (6)

1. the method for building up of a earth soil solid resistance network model, it is characterised in that under the method includes Row step:
(1) the earth soil is simplified to multiple soil grid；
(2) geological conditions of the earth soil is obtained；
(3) limited element calculation model is built according to the geological conditions of the earth soil obtained in step (2)；
(4) current electric fields distribution in the earth soil is calculated according to the limited element calculation model obtained in step (3)；
(5) according to the current electric fields distribution obtained in step (4), calculate soil grid Normal plane radially and put down All current potential UrAnd normal direction current phasor and
(6) according to Normal plane average potential U of the soil grid radial direction obtained in step (5)rAnd normal direction Current phasor andCalculate the radial direction substitutional resistance R that soil grid is correspondingr, width is to substitutional resistance RθAnd depth To substitutional resistance Rd；
(7) according to the radial direction substitutional resistance R obtained in step (6)r, width is to substitutional resistance RθAnd depth to Substitutional resistance Rd, by abstract for the soil grid in step (1) be three-dimensional resistor network, set up three-dimensional resistor network Model.
The method for building up of the earth soil solid resistance network model the most according to claim 1, its feature exists In, the geological conditions of described the earth soil includes stratification state, landform and soil resistivity.
The method for building up of the earth soil solid resistance network model the most according to claim 1, its feature exists In, described soil grid is fan-shaped.
The method for building up of the earth soil solid resistance network model the most according to claim 1, its feature exists In, in described step (5), the radial direction substitutional resistance R of i-th soil gridriResistance be:
R r i = U r i - U r ( i + 1 ) 1 2 ( I · r i + I · r ( i + 1 ) ) ;
Wherein, Uri、Ur(i+1)It is respectively i-th and i+1 soil grid Normal plane average potential radially；Be respectively i-th and i+1 soil grid normal direction current phasor radially with.
The method for building up of the earth soil solid resistance network model the most according to claim 4, its feature exists In, in described step (5), the width of i-th soil grid is to substitutional resistance RθiResistance be:
R θ i = 2 πθR r i r 1 i ( r 2 i - r 1 i ) * ln r 2 i r 1 i ;
Wherein, θ is the soil grid angle relative to the center of circle, r1i、r2iIt is respectively the interior outer radius of i-th soil grid.
The method for building up of the earth soil solid resistance network model the most according to claim 4, its feature exists In, in described step (5), the depth of i-th soil grid is to substitutional resistance RdiResistance be:
R d i = R r i d 2 ( r 2 i 2 - r 1 i 2 ) * ln r 2 i r 1 i ;
Wherein, d is the degree of depth of i-th soil grid, r1i、r2iIt is respectively the interior outer radius of i-th soil grid.
CN201610301434.4A 2016-05-09 2016-05-09 A kind of method for building up of the earth soil solid resistance network model Active CN106021652B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610301434.4A CN106021652B (en) 2016-05-09 2016-05-09 A kind of method for building up of the earth soil solid resistance network model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610301434.4A CN106021652B (en) 2016-05-09 2016-05-09 A kind of method for building up of the earth soil solid resistance network model

Publications (2)

Publication Number Publication Date
CN106021652A true CN106021652A (en) 2016-10-12
CN106021652B CN106021652B (en) 2019-03-08

Family

ID=57098837

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610301434.4A Active CN106021652B (en) 2016-05-09 2016-05-09 A kind of method for building up of the earth soil solid resistance network model

Country Status (1)

Country Link
CN (1) CN106021652B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108388707A (en) * 2018-02-05 2018-08-10 三峡大学 D.C. magnetic biasing computational methods based on field circuit method under a kind of three-dimensional asymmetric structure soil model
CN111756026A (en) * 2019-03-29 2020-10-09 中国能源建设集团江苏省电力设计院有限公司 Transformer direct-current magnetic bias suppression method based on multiple soil models
CN112668212A (en) * 2020-09-02 2021-04-16 国网内蒙古东部电力有限公司检修分公司 Finite element-based method for analyzing overflow characteristics of grounding electrode under different soil models
CN115267337A (en) * 2022-09-22 2022-11-01 国网山西省电力公司电力科学研究院 Method and device for determining wide area soil resistivity and storage medium
CN115823503A (en) * 2022-11-29 2023-03-21 北京市燃气集团有限责任公司 Community gas pipeline corrosion hotspot field test and identification method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20130062435A (en) * 2011-11-18 2013-06-13 (주)유샘인스트루먼트 Soil water management system
CN103593523A (en) * 2013-11-12 2014-02-19 国网上海市电力公司 Finite element theory based direct current magnetic bias suppression method under condition of multiple direct-current falling points
US20140236868A1 (en) * 2013-02-15 2014-08-21 Banyan Water, Inc. System and method for automated, range-based irrigation
CN105005676A (en) * 2015-08-13 2015-10-28 国网上海市电力公司 Three-dimension design method based on cable engineering information model

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20130062435A (en) * 2011-11-18 2013-06-13 (주)유샘인스트루먼트 Soil water management system
US20140236868A1 (en) * 2013-02-15 2014-08-21 Banyan Water, Inc. System and method for automated, range-based irrigation
CN103593523A (en) * 2013-11-12 2014-02-19 国网上海市电力公司 Finite element theory based direct current magnetic bias suppression method under condition of multiple direct-current falling points
CN105005676A (en) * 2015-08-13 2015-10-28 国网上海市电力公司 Three-dimension design method based on cable engineering information model

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
傅晨钊 等: "基于有限元分析的直流偏磁电流计算模型研究", 《华东电力》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108388707A (en) * 2018-02-05 2018-08-10 三峡大学 D.C. magnetic biasing computational methods based on field circuit method under a kind of three-dimensional asymmetric structure soil model
CN111756026A (en) * 2019-03-29 2020-10-09 中国能源建设集团江苏省电力设计院有限公司 Transformer direct-current magnetic bias suppression method based on multiple soil models
CN111756026B (en) * 2019-03-29 2024-04-12 中国能源建设集团江苏省电力设计院有限公司 Transformer direct-current magnetic bias suppression method based on multiple soil models
CN112668212A (en) * 2020-09-02 2021-04-16 国网内蒙古东部电力有限公司检修分公司 Finite element-based method for analyzing overflow characteristics of grounding electrode under different soil models
CN112668212B (en) * 2020-09-02 2023-03-24 国网内蒙古东部电力有限公司检修分公司 Finite element-based method for analyzing overflow characteristics of grounding electrode under different soil models
CN115267337A (en) * 2022-09-22 2022-11-01 国网山西省电力公司电力科学研究院 Method and device for determining wide area soil resistivity and storage medium
CN115267337B (en) * 2022-09-22 2022-12-09 国网山西省电力公司电力科学研究院 Method and device for determining wide area soil resistivity and storage medium
CN115823503A (en) * 2022-11-29 2023-03-21 北京市燃气集团有限责任公司 Community gas pipeline corrosion hotspot field test and identification method

Also Published As

Publication number Publication date
CN106021652B (en) 2019-03-08

Similar Documents

Publication Publication Date Title
CN106021652A (en) An earth soil three-dimensional resistance network model establishing method
CN101900766A (en) Method for determining position of voltage electrode in measurement of grounding resistance in transformer substation by short-distance method
CN104112057B (en) A kind of large scale fractured reservoir method for numerical simulation
CN108763825B (en) Numerical simulation method for simulating wind field of complex terrain
CN102609982B (en) Topology discovery method of space geological data based on unstructured mode
CN103018633B (en) Based on dissymmetrical impedance model short-circuit current fault signature containing the interval decision method of DG distribution network failure
CN101900764A (en) Method for measuring grounded resistance of grounded screen by short range measuring method
CN108021780A (en) A kind of mountain torrents dynamic emulation method based on random unstrctured grid model
El Kasmi et al. Turbulence modeling of atmospheric boundary layer flow over complex terrain: a comparison of models at wind tunnel and full scale
CN110515117A (en) A kind of underground karst cavity detection method to be impulsed based on pile monkey with decision-tree model
CN103777248A (en) TEM one-dimensional forward modeling method applicable to irregular transmitting loop
Ma et al. Analysis of grounding systems in soils with cylindrical soil volumes
CN104361157A (en) Evaluation method for wind environment between buildings
Dong et al. Evaluating the applicability of the finite element method for modelling of geoelectric fields
CN105490041A (en) Electric transmission line tower grounding body and impulse grounding resistance calculation method therefor
CN107451372A (en) A kind of flood of a mountain area numerical simulation method that kinematic wave is combined with dynamic wave
CN105261067A (en) Overground and underground integrated modeling method based on true 3D volumetric display technique and system
CN108446471B (en) Mountain area electromagnetic field prediction technique based on three-dimensional moment method and two-dimentional Fast Multiple Method
CN101900767A (en) Visual method for utilizing short distance measuring method to measure ground resistance of analogue ground screen
CN106202704A (en) A kind of D.C. magnetic biasing impact evaluation method of determining range
Huang et al. Simulation of groundwater flow in fractured rocks using a coupled model based on the method of domain decomposition
CN106326579A (en) Method for generating visualizing of land and sea environment based on electromagnetic scattering technique
CN101900765A (en) Method for building equivalent grounding network computation model
CN104217010A (en) Conducting wire trip-out rate calculating and evaluating method based on segmented topography
CN103678883A (en) Space sampling method oriented to multisource marine environmental monitoring data

Legal Events

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