CN112052555A - Simulation method and device for lightning electromagnetic transient model of power transmission line - Google Patents

Simulation method and device for lightning electromagnetic transient model of power transmission line Download PDF

Info

Publication number
CN112052555A
CN112052555A CN202010728774.1A CN202010728774A CN112052555A CN 112052555 A CN112052555 A CN 112052555A CN 202010728774 A CN202010728774 A CN 202010728774A CN 112052555 A CN112052555 A CN 112052555A
Authority
CN
China
Prior art keywords
model
equivalent circuit
node
segment
lightning
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
CN202010728774.1A
Other languages
Chinese (zh)
Other versions
CN112052555B (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.)
China South Power Grid International Co ltd
Guangzhou Power Supply Bureau of Guangdong Power Grid Co Ltd
Original Assignee
China South Power Grid International Co ltd
Guangzhou Power Supply Bureau of Guangdong Power Grid 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 China South Power Grid International Co ltd, Guangzhou Power Supply Bureau of Guangdong Power Grid Co Ltd filed Critical China South Power Grid International Co ltd
Priority to CN202010728774.1A priority Critical patent/CN112052555B/en
Publication of CN112052555A publication Critical patent/CN112052555A/en
Application granted granted Critical
Publication of CN112052555B publication Critical patent/CN112052555B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/04—Power grid distribution networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • 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)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a method and a device for simulating a lightning electromagnetic transient model of a power transmission line, wherein the method comprises the following steps: acquiring a preset partial unit equivalent circuit model; the equivalent circuit model of the part of the units is an equivalent circuit model obtained by correcting the equivalent circuit model of the part of the units by utilizing a Green function according to a formula of the equivalent circuit of the part of the units; constructing an overhead line model of the power transmission network according to the equivalent circuit model of the part of the units; and carrying out lightning induced overvoltage simulation on the overhead line model to obtain the lightning resistance level and the spatial electromagnetic field distribution condition of the overhead line model. The method utilizes the Green function to simulate the influence of the lossy earth to correct the existing PEEC method, thereby obtaining an improved partial unit equivalent circuit model and further more accurately simulating the lightning influence of the power transmission line.

Description

Simulation method and device for lightning electromagnetic transient model of power transmission line
Technical Field
The invention relates to the technical field of lightning electromagnetic transient simulation of a power transmission line, in particular to a method and a device for simulating a lightning electromagnetic transient model of the power transmission line.
Background
Lightning is one of the most prominent external factors responsible for distribution line failure. Conductors such as a conducting wire, a ground wire and a pole tower in a distribution line system can be regarded as a multi-conductor system, and lightning transient in the multi-conductor system is a long-standing problem. The multi-conductor system can be run in any direction and interconnected to form a three-dimensional wire structure on lossy ground. When a system is struck by a lightning, a lightning strike current is generated in the conductor. In order to effectively protect against the effects of lightning strikes, it is necessary to build appropriate models to model and evaluate lightning strike transients in the structure.
A partial cell equivalent circuit (PEEC) method is a method for modeling a three-dimensional interconnected fine line structure. It is derived from the mixed potential integral equation and converts the electromagnetic problem into a circuit domain. In recent years, this method has been applied to handling transient currents and voltages in various configurations of lines in any direction, but significant errors may occur because the effect of damaging earth is not taken into account.
Disclosure of Invention
The embodiment of the invention aims to provide a method and a device for simulating a lightning electromagnetic transient model of a power transmission line, which are used for simulating the influence of a lossy earth by utilizing a Green function so as to modify the existing PEEC method, thereby obtaining an improved partial unit equivalent circuit model and further more accurately simulating the lightning influence of the power transmission line.
In order to achieve the above object, an embodiment of the present invention provides a method for simulating a lightning electromagnetic transient model of a power transmission line, including the following steps:
acquiring a preset partial unit equivalent circuit model; the equivalent circuit model of the part of the units is an equivalent circuit model obtained by correcting the equivalent circuit model of the part of the units by utilizing a Green function according to a formula of the equivalent circuit of the part of the units;
constructing an overhead line model of the power transmission network according to the equivalent circuit model of the part of the units;
and carrying out lightning induced overvoltage simulation on the overhead line model to obtain the lightning resistance level and the spatial electromagnetic field distribution condition of the overhead line model.
Preferably, the node voltage of the partial unit equivalent circuit model is Vn＝∑kpnkqkt+∑tjωcntItAnd node capacitance current of
Figure BDA0002598621270000021
Wherein, VnNode voltage segmented for nth node, Ic,nNode capacitance current, p, segmented for the nth nodenkPotential coefficients of conductor segments from the n-th node segment to the k-th node segment, qnkThe amount of charge of the conductor segments for the n-th to k-th node segments, cntIs a preset correction coefficient, omega is angular frequency, the partial unit equivalent circuit model comprises a plurality of node segments, the node segments comprise a field segment and a source segment, ItIs the current of the source section t, pnnAnd (4) segmenting the potential coefficient of the nth node pair, wherein k is not equal to n.
Preferably, the potential coefficient p of the conductor segments of the n-th to k-th node segmentsnkIs calculated by the formula
Figure BDA0002598621270000022
Wherein s isnAnd skSectional areas l corresponding to the nth node segment and the kth node segment, respectivelynAnd lkRespectively corresponding lengths, K, of the nth and the kth node segmentsφAnd (r, r ') is a scalar green function corresponding to the conductor surface node r ═ r ', r is the distance from the nth node segment to a preset observation point, and r ' is the distance from the kth node segment to the observation point.
Preferably, the correction coefficient cntIs calculated by the formula
Figure BDA0002598621270000023
Wherein，snAnd stRespectively corresponding to the nth node segment and the source segment tnAnd ltThe lengths of the n-th node segment and the source segment t are respectively, and P (r, r') is all PnkThe electric potential coefficient matrix is formed by the components,
Figure BDA0002598621270000024
is a unit vector in the z-direction of the vertical field, l'tFor the length l of the source segment ttThe corresponding differentiation.
Another embodiment of the present invention provides a device for simulating a lightning electromagnetic transient model of a power transmission line, the device including:
the model acquisition module is used for acquiring a preset partial unit equivalent circuit model; the equivalent circuit model of the part of the units is an equivalent circuit model obtained by correcting the equivalent circuit model of the part of the units by utilizing a Green function according to a formula of the equivalent circuit of the part of the units;
the line model building module is used for building an overhead line model of the power transmission network according to the part of unit equivalent circuit models;
and the simulation module is used for carrying out lightning induction overvoltage simulation on the overhead line model to obtain the lightning resistance level and the spatial electromagnetic field distribution condition of the overhead line model.
The invention further provides a device using a simulation method of a power transmission line lightning electromagnetic transient model, which comprises a processor, a memory and a computer program stored in the memory and configured to be executed by the processor, wherein the processor implements the simulation method of the power transmission line lightning electromagnetic transient model according to any one of the above items when executing the computer program.
Compared with the prior art, the method and the device for simulating the lightning electromagnetic transient model of the power transmission line provided by the embodiment of the invention simulate the influence of the lossy earth by utilizing the Green function to correct the existing PEEC method, so that an improved partial unit equivalent circuit model is obtained, and the lightning influence of the power transmission line can be simulated more accurately.
Drawings
Fig. 1 is a schematic flow chart of a simulation method of a lightning electromagnetic transient model of a power transmission line according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a structure of a conductor segment in an air-ground medium according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of an equivalent circuit of a segment of conductive wire obtained by the method of the present invention according to an embodiment of the present invention;
FIG. 4 is a schematic diagram of a simulation apparatus for a lightning electromagnetic transient model of a power transmission line according to an embodiment of the present invention;
fig. 5 is a schematic diagram of an apparatus using a simulation method of a lightning electromagnetic transient model of a transmission line according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Referring to fig. 1, which is a schematic flow chart of a method for simulating a lightning electromagnetic transient model of a power transmission line according to an embodiment of the present invention, the method includes steps S1 to S3:
s1, acquiring a preset partial unit equivalent circuit model; the equivalent circuit model of the part of the units is an equivalent circuit model obtained by correcting the equivalent circuit model of the part of the units by utilizing a Green function according to a formula of the equivalent circuit of the part of the units;
s2, constructing an overhead line model of the power transmission network according to the partial unit equivalent circuit model;
and S3, performing lightning induced overvoltage simulation on the overhead line model to obtain the lightning resistance level and the spatial electromagnetic field distribution condition of the overhead line model.
It should be noted that a Partial Element Equivalent Circuit (PEEC) method is a method for modeling a three-dimensional interconnected thin line structure, which is derived from a mixed potential integral equation and converts an electromagnetic problem into a Circuit domain. This method can be applied to handle transient currents and voltages in various configurations of lines in any direction, but does not take into account the effects of damaging the earth and may produce significant errors. The invention utilizes the Green function to correct the mixed potential integral equation in the PEEC formula so as to simulate the influence of the damaged earth, thereby obtaining an improved partial unit equivalent circuit model.
Specifically, a preset partial unit equivalent circuit model is obtained; the equivalent circuit model of the partial unit is obtained by correcting the equivalent circuit model of the partial unit according to a partial equivalent circuit formula by utilizing a Green function;
constructing an overhead line model of the power transmission network according to the partial equivalent circuit model;
and performing lightning induced overvoltage simulation on the overhead line model to obtain the lightning resistance level and the spatial electromagnetic field distribution condition of the overhead line model.
In order to more clearly understand the process of constructing the partial unit equivalent circuit model of the present invention, the following describes the improvement process of the partial unit equivalent circuit model in detail.
The magnetic vector potential a and the scalar potential phi are often used to simplify the mathematical expression of maxwell's equations, satisfying the following at the conductor surface r point:
Figure BDA0002598621270000051
where E is the electric field on the surface of the conductor and ω is the angular frequency.
Assuming the presence of air (0，μ0And σ0) And isotropic lossy floor (1，μ1And σ1) The layered structure of (a) is,0and1dielectric coefficient of air and dielectric coefficient of ground, mu, respectively0And mu1Respectively permeability of air and ground, sigma0And σ1Air conductivity and ground conductivity, respectively, see the figure2, the structure of the conductor segment in the air-ground medium provided by the embodiment of the invention is schematically shown. For conductors located in air, both the magnetic vector potential a and the scalar potential Φ can be expressed in a gray function as shown below:
Figure BDA0002598621270000052
Figure BDA0002598621270000053
where I (r ') and τ (r ') are the conductor current and line charge density, respectively, at r ═ r ' on the conductor. In the formula (I), the compound is shown in the specification,
Figure BDA0002598621270000054
for the dyadic green function of the magnetic vector potential a, there are different components to the current dipole in three orthogonal directions:
Figure BDA0002598621270000055
Figure BDA0002598621270000056
wherein g (r, r ') is a homogeneous medium three-dimensional green's function,
Figure BDA0002598621270000057
wherein R is | R-R '|, R is the distance from the nth node segment to the observation point, R' is the distance from the kth node segment to the observation point, the observation point is any point on the lead and is generally preset manually, d is a calculation coefficient,
Figure BDA0002598621270000058
is expressed as
Figure BDA0002598621270000059
In the formula (I), the compound is shown in the specification,
Figure BDA00025986212700000510
being the unit vector of the current dipole projected on the xy plane or ground,
Figure BDA00025986212700000511
is the unit vector in the x-direction of the horizontal field,
Figure BDA00025986212700000512
being the unit vector in the y-direction of the horizontal field,
Figure BDA0002598621270000061
is a unit vector of the z direction of the vertical field, GxxIs composed of
Figure BDA0002598621270000062
In that
Figure BDA0002598621270000063
Component of dyadic over component, GyyIs composed of
Figure BDA0002598621270000064
In that
Figure BDA0002598621270000065
Component of dyadic over component, GzzIs composed of
Figure BDA0002598621270000066
In that
Figure BDA0002598621270000067
Component of dyadic over component, GzxIs composed of
Figure BDA0002598621270000068
In that
Figure BDA0002598621270000069
Component of dyadic over component, GzyIs composed of
Figure BDA00025986212700000610
In that
Figure BDA00025986212700000611
Component of dyadic over component, GttIs composed of
Figure BDA00025986212700000612
In that
Figure BDA00025986212700000613
Component of dyadic over component, GztIs composed of
Figure BDA00025986212700000614
In that
Figure BDA00025986212700000615
Component of dyadic over the components, and Gxx＝Gyy＝GttL 'is the segment length of the source conductor, and P (r, r') is all PnkA matrix of potential coefficients of composition, KφIs a scalar gray function. For current dipoles oriented in a stratified medium, there is no unique scalar Green's function, K in the present inventionφIs determined from the integration of g (r, r') and Sofaffy.
Dividing the conductor in the air into a plurality of sections, and selecting two sections: field segment i and source segment t, then integrating the formula along the field segment yields the potential equation as follows:
Figure BDA00025986212700000616
Figure BDA00025986212700000617
in the formula, si、liAnd JiRespectively the cross-sectional area, the length and the current density of the field section i,/t、l′t、ItRespectively the length of the source section t, the length differential and the conductor current, phinAnd phimRespectively segment the nth node and the mth node on the field segment iSegmented node voltage, τkAnd l'kThe line charge density and length differential of the kth node segment on field segment i, respectively, σ is the conductivity, n ≠ m, n ≠ k.
From the definition of the average node potential (voltage) V, a circuit equation can be established as follows:
Zs,iI+∑tjωLitIt＝Vn-Vm
Vn＝∑kpnkqkt+∑tjωcntIt
Zs,iis the surface impedance of the conductor field segment i, and is expressed as
Figure BDA00025986212700000618
RaIs the internal impedance of the conductor field section i, for a conductivity σ, a permeability μ0A circular conductor of radius a,
Figure BDA00025986212700000619
i is a first class of higher order Bessel function, I0And I1Respectively, to their 0 th and 1 st order solutions. I istIs the current of the conductor source section t, the parameter LitInductance, p, for the segmentation of the field segments i to the source segments tnkPotential coefficient of segments n to knkAmount of charge for node segment n to node segment k, cntIs the correction factor.
According to the above analysis, there are
Figure BDA0002598621270000071
Figure BDA0002598621270000072
Figure BDA0002598621270000073
Taking the time derivative of the potential to obtain the capacitance current on the node in the frequency domain as:
Figure BDA0002598621270000074
in the formula, qnIs the amount of charge of the node n, Ic,kThe node capacitance current of the kth node.
Node voltage VnAnd branch current IiThe two sets of unknown numbers in the equivalent circuit of the partial unit are both solved by using a conventional circuit analysis tool.
According to the method for simulating the lightning electromagnetic transient model of the power transmission line, which is provided by the embodiment 1 of the invention, the green function is utilized to simulate the influence of the lossy earth to modify the existing PEEC method, so that an improved partial unit equivalent circuit model is obtained, and the lightning influence of the power transmission line can be simulated more accurately.
As an improvement of the scheme, the node voltage of the equivalent circuit model of the partial unit is Vn＝∑kpnkqkt+∑tjωcntItAnd node capacitance current of
Figure BDA0002598621270000075
Wherein, VnNode voltage segmented for nth node, Ic,nNode capacitance current, p, segmented for the nth nodenkPotential coefficients of conductor segments from the n-th node segment to the k-th node segment, qnkThe amount of charge of the conductor segments for the n-th to k-th node segments, cntIs a preset correction coefficient, omega is angular frequency, the partial unit equivalent circuit model comprises a plurality of node segments, the node segments comprise a field segment and a source segment, ItIs the current of the source section t, pnnAnd (4) segmenting the potential coefficient of the nth node pair, wherein k is not equal to n.
Referring to fig. 3, it is a schematic diagram of an equivalent circuit of a segment of conducting wire obtained by the method of the present invention according to the embodiment of the present invention,that is, the equivalent circuit of a part of the unit obtained by the improvement of the present invention is shown in fig. 3. Inductor LitIncluding coupling elements from horizontal (source) segments to vertical (field) segments, which is not present in the conventional definition of inductance. The introduction of this coupling element is purely due to the presence of the lossy earth, a correction term for the mutual inductive coupling. However, due to the asymmetric matrix
Figure BDA0002598621270000076
The vertical current segments do not have any inductive contribution to the horizontal current segments, which is a unidirectional coupling, while the vertical current segments do contribute to the node potential.
Specifically, the node voltage of the partial unit equivalent circuit model is Vn＝∑kpnkqkt+∑tjωcntItAnd node capacitance current of
Figure BDA0002598621270000081
Wherein, VnNode voltage segmented for nth node, Ic,nNode capacitance current, p, segmented for the nth nodenkPotential coefficients of conductor segments from the n-th node segment to the k-th node segment, qnkThe amount of charge of the conductor segments for the n-th to k-th node segments, cntIs a preset correction coefficient, omega is angular frequency, the partial unit equivalent circuit model comprises a plurality of node segments, the node segments comprise a field segment and a source segment, ItIs the current of the source section t, pnnAnd (4) segmenting the potential coefficient of the nth node pair, wherein k is not equal to n.
As an improvement of the above, the potential coefficient p of the conductor segments of the n-th to k-th node segmentsnkIs calculated by the formula
Figure BDA0002598621270000082
Wherein s isnAnd skSectional areas l corresponding to the nth node segment and the kth node segment, respectivelynAnd lkRespectively corresponding lengths, K, of the nth and the kth node segmentsφ(r, r') isAnd the body surface point r is a scalar Green function corresponding to r ', r is the distance from the nth node segment to a preset observation point, and r' is the distance from the kth node segment to the observation point.
In particular, the potential coefficient p of the conductor segments of the n-th to k-th node segmentsnkIs calculated by the formula
Figure BDA0002598621270000083
Wherein s isnAnd skSectional areas l corresponding to the nth node segment and the kth node segment, respectivelynAnd lkRespectively corresponding lengths, K, of the nth and the kth node segmentsφAnd (r, r ') is a scalar Green function corresponding to the conductor surface point r ═ r ', r is the distance from the nth node segment to the preset observation point, and r ' is the distance from the kth node segment to the observation point.
As a modification of the above, the correction coefficient cntIs calculated by the formula
Figure BDA0002598621270000084
Wherein s isnAnd stRespectively corresponding to the nth node segment and the source segment tnAnd ltThe lengths of the n-th node segment and the source segment t are respectively, and P (r, r') is all PnkThe electric potential coefficient matrix is formed by the components,
Figure BDA0002598621270000085
is a unit vector in the z-direction of the vertical field, l'tFor the length l of the source segment ttThe corresponding differentiation.
Specifically, the correction coefficient cntIs calculated by the formula
Figure BDA0002598621270000091
Wherein s isnAnd stRespectively corresponding to the nth node segment and the source segment tnAnd ltThe lengths corresponding to the nth node segment and the source segment t, respectively, and P (r, r') is all PnkThe electric potential coefficient matrix is formed by the components,
Figure BDA0002598621270000092
is a unit vector in the z-direction of the vertical field, l'tIs the length l of the source segment ttThe corresponding differentiation.
Referring to fig. 4, a simulation apparatus for a lightning electromagnetic transient model of a power transmission line according to an embodiment of the present invention is characterized by including:
the model obtaining module 11 is configured to obtain a preset partial unit equivalent circuit model; the equivalent circuit model of the part of the units is an equivalent circuit model obtained by correcting the equivalent circuit model of the part of the units by utilizing a Green function according to a formula of the equivalent circuit of the part of the units;
a line model construction module 12, configured to construct an overhead line model of the power transmission network according to the partial unit equivalent circuit model;
and the simulation module 13 is configured to perform lightning induced overvoltage simulation on the overhead line model to obtain a lightning resistance level and a spatial electromagnetic field distribution condition of the overhead line model.
The simulation device of the power transmission line lightning electromagnetic transient model provided by the embodiment of the invention can realize all processes of the simulation method of the power transmission line lightning electromagnetic transient model described in any one of the embodiments, and the functions and the realized technical effects of each module and unit in the device are respectively the same as those of the simulation method of the power transmission line lightning electromagnetic transient model described in the embodiment and the realized technical effects, and are not repeated here.
Referring to fig. 5, the schematic diagram of an apparatus for a simulation method using a transmission line lightning electromagnetic transient model according to the embodiment of the present invention includes a processor 10, a memory 20, and a computer program stored in the memory 20 and configured to be executed by the processor 10, where the processor 10 implements the simulation method of the transmission line lightning electromagnetic transient model according to any one of the above embodiments when executing the computer program.
Illustratively, the computer program may be divided into one or more modules/units, which are stored in the memory 20 and executed by the processor 10 to implement the present invention. One or more modules/units may be a series of computer program instruction segments capable of performing specific functions, and the instruction segments are used for describing the execution process of the computer program in a simulation method of the lightning electromagnetic transient model of the power transmission line. For example, the computer program may be divided into a model acquisition module, a line model construction module and a simulation module, and the specific functions of each module are as follows:
the model obtaining module 11 is configured to obtain a preset partial unit equivalent circuit model; the equivalent circuit model of the part of the units is an equivalent circuit model obtained by correcting the equivalent circuit model of the part of the units by utilizing a Green function according to a formula of the equivalent circuit of the part of the units;
a line model construction module 12, configured to construct an overhead line model of the power transmission network according to the partial unit equivalent circuit model;
and the simulation module 13 is configured to perform lightning induced overvoltage simulation on the overhead line model to obtain a lightning resistance level and a spatial electromagnetic field distribution condition of the overhead line model.
The device using the simulation method of the lightning electromagnetic transient model of the power transmission line can be computing equipment such as a desktop computer, a notebook computer, a palm computer and a cloud server. The device using the simulation method of the power transmission line lightning electromagnetic transient model can comprise, but is not limited to, a processor and a memory. It will be understood by those skilled in the art that the schematic diagram 5 is merely an example of an apparatus for a simulation method using a transmission line lightning electromagnetic transient model, and does not constitute a limitation of the apparatus for the simulation method using the transmission line lightning electromagnetic transient model, and may include more or less components than those shown, or combine some components, or different components, for example, the apparatus for the simulation method using the transmission line lightning electromagnetic transient model may further include an input and output device, a network access device, a bus, and the like.
The Processor 10 may be a Central Processing Unit (CPU), other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), an off-the-shelf Programmable Gate Array (FPGA) or other Programmable logic device, discrete Gate or transistor logic, discrete hardware components, etc. The general purpose processor may be a microprocessor or the processor 10 may be any conventional processor or the like, and the processor 10 is a control center of the apparatus using the simulation method of the transmission line lightning electromagnetic transient model, and various interfaces and lines are used to connect various parts of the entire apparatus using the simulation method of the transmission line lightning electromagnetic transient model.
The memory 20 may be used to store the computer programs and/or modules, and the processor 10 implements various functions of the apparatus using the simulation method of the transmission line lightning electromagnetic transient model by running or executing the computer programs and/or modules stored in the memory 20 and calling the data stored in the memory 20. The memory 20 may mainly include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required by at least one function (such as a sound playing function, an image playing function, etc.), and the like; the storage data area may store data (such as audio data, a phonebook, etc.) created according to the use of the cellular phone, and the like. In addition, the memory 20 may include high speed random access memory, and may also include non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), at least one magnetic disk storage device, a Flash memory device, or other volatile solid state storage device.
The module integrated by the device using the simulation method of the lightning electromagnetic transient model of the power transmission line can be stored in a computer readable storage medium if the module is realized in the form of a software functional unit and sold or used as an independent product. Based on such understanding, all or part of the flow of the method according to the embodiments of the present invention may also be implemented by a computer program, which may be stored in a computer-readable storage medium and can implement the steps of the embodiments of the method when the computer program is executed by a processor. The computer program includes computer program code, and the computer program code may be in a source code form, an object code form, an executable file or some intermediate form. The computer readable medium may include: any entity or device capable of carrying computer program code, recording medium, U.S. disk, removable hard disk, magnetic disk, optical disk, computer Memory, Read-Only Memory (ROM), Random Access Memory (RAM), electrical carrier wave signals, telecommunications signals, software distribution media, and the like. It should be noted that the computer readable medium may contain other components which may be suitably increased or decreased as required by legislation and patent practice in jurisdictions, for example, in some jurisdictions, in accordance with legislation and patent practice, the computer readable medium does not include electrical carrier signals and telecommunications signals.
The embodiment of the present invention further provides a computer-readable storage medium, where the computer-readable storage medium includes a stored computer program, and when the computer program runs, the device where the computer-readable storage medium is located is controlled to execute the method for simulating the lightning electromagnetic transient model of the power transmission line according to any one of the above embodiments.
To sum up, the method and the device for simulating the lightning electromagnetic transient model of the power transmission line provided by the embodiment of the invention analyze the full-wave PEEC formula of the lightning transient analysis of any directional conductor on the lossy ground, simulate the influence of the lossy foundation by using the dyadic Green function, and finally provide a new PEEC equivalent circuit by considering the correction term caused by the lossy grounding. The improved PEEC equivalent circuit is used for modeling, a distribution network overhead line is established, the electromagnetic coupling effect with the damage to the ground is calculated and processed, the conversion from an electric field to a circuit is realized, the lightning induced overvoltage of the overhead line is modeled and simulated, the lightning resistance level of the overhead line is tested, the space electromagnetic field is analyzed, the induced voltage calculation modeling of a complex network is realized, and the reference for the construction and operation of a distribution network can be realized. The establishment of the invention is beneficial to realizing the full coverage modeling analysis of the power distribution network and provides a new simulation means for lightning protection work. The equivalent circuit considering the influence of the lossy earth can develop a new application scene of the PEEC calculation method, and is applied to the lightning electromagnetic transient problem of overhead conductors, overhead ground wires and towers under the influence of soil electromagnetic coupling.
While the foregoing is directed to the preferred embodiment of the present invention, it will be understood by those skilled in the art that various changes and modifications may be made without departing from the spirit and scope of the invention.

Claims (6)

1. A simulation method of a lightning electromagnetic transient model of a power transmission line is characterized by comprising the following steps:
acquiring a preset partial unit equivalent circuit model; the equivalent circuit model of the part of the units is an equivalent circuit model obtained by correcting the equivalent circuit model of the part of the units by utilizing a Green function according to a formula of the equivalent circuit of the part of the units;
constructing an overhead line model of the power transmission network according to the equivalent circuit model of the part of the units;
and carrying out lightning induced overvoltage simulation on the overhead line model to obtain the lightning resistance level and the spatial electromagnetic field distribution condition of the overhead line model.
2. The method of claim 1, wherein the node voltage of the partial cell equivalent circuit model is Vn＝∑kpnkqkt+∑tjωcntItAnd node capacitance current of
Figure FDA0002598621260000011
Wherein, VnNode voltage segmented for nth node, Ic,nNode capacitance current, p, segmented for the nth nodenkThe potential coefficients of the conductor segments from the n-th node segment to the k-th node segment，qnkThe amount of charge of the conductor segments for the n-th to k-th node segments, cntIs a preset correction coefficient, omega is angular frequency, the partial unit equivalent circuit model comprises a plurality of node segments, the node segments comprise a field segment and a source segment, ItIs the current of the source section t, pnnAnd (4) segmenting the potential coefficient of the nth node pair, wherein k is not equal to n.
3. Method for simulating a lightning electromagnetic transient model of a transmission line according to claim 2, characterized in that the potential coefficient p of the conductor segment from the n-th node segment to the k-th node segmentnkIs calculated by the formula
Figure FDA0002598621260000012
Wherein s isnAnd skSectional areas l corresponding to the nth node segment and the kth node segment, respectivelynAnd lkRespectively corresponding lengths, K, of the nth and the kth node segmentsφAnd (r, r ') is a scalar green function corresponding to the conductor surface node r ═ r ', r is the distance from the nth node segment to a preset observation point, and r ' is the distance from the kth node segment to the observation point.
4. Method for simulating a lightning electromagnetic transient model of a transmission line according to claim 3, characterised in that said correction factor cntIs calculated by the formula
Figure FDA0002598621260000021
Wherein s isnAnd stRespectively corresponding to the nth node segment and the source segment tnAnd ltThe lengths of the n-th node segment and the source segment t are respectively, and P (r, r') is all PnkThe electric potential coefficient matrix is formed by the components,
Figure FDA0002598621260000022
is a unit vector in the z-direction of the vertical field, l'tIs the length of the source segment tDegree ltThe corresponding differentiation.
5. The utility model provides a simulation device of transmission line thunder and lightning electromagnetism transient state model which characterized in that includes:
the model acquisition module is used for acquiring a preset partial unit equivalent circuit model; the equivalent circuit model of the part of the units is an equivalent circuit model obtained by correcting the equivalent circuit model of the part of the units by utilizing a Green function according to a formula of the equivalent circuit of the part of the units;
the line model building module is used for building an overhead line model of the power transmission network according to the part of unit equivalent circuit models;
and the simulation module is used for carrying out lightning induction overvoltage simulation on the overhead line model to obtain the lightning resistance level and the spatial electromagnetic field distribution condition of the overhead line model.
6. An apparatus using a simulation method of a transmission line lightning electromagnetic transient model, characterized by comprising a processor, a memory and a computer program stored in the memory and configured to be executed by the processor, the processor implementing the simulation method of a transmission line lightning electromagnetic transient model according to any one of claims 1 to 4 when executing the computer program.
CN202010728774.1A 2020-07-23 2020-07-23 Simulation method and device for lightning electromagnetic transient model of power transmission line Active CN112052555B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010728774.1A CN112052555B (en) 2020-07-23 2020-07-23 Simulation method and device for lightning electromagnetic transient model of power transmission line

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010728774.1A CN112052555B (en) 2020-07-23 2020-07-23 Simulation method and device for lightning electromagnetic transient model of power transmission line

Publications (2)

Publication Number Publication Date
CN112052555A true CN112052555A (en) 2020-12-08
CN112052555B CN112052555B (en) 2024-02-23

Family

ID=73602779

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010728774.1A Active CN112052555B (en) 2020-07-23 2020-07-23 Simulation method and device for lightning electromagnetic transient model of power transmission line

Country Status (1)

Country Link
CN (1) CN112052555B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112731076A (en) * 2020-12-17 2021-04-30 南方电网科学研究院有限责任公司 Early warning method and device based on insulation flashover rate under overvoltage and storage medium
CN112986659A (en) * 2021-02-01 2021-06-18 中国民航大学 Method for analyzing voltage drop in composite material airplane ground return network
CN113075495A (en) * 2021-02-24 2021-07-06 南方电网科学研究院有限责任公司 Circuit parameter adjusting method, device, equipment and medium of PEEC model
CN113191103A (en) * 2021-03-05 2021-07-30 南方电网科学研究院有限责任公司 Method and device for calculating lightning electromagnetic transient of distribution line
CN113406726A (en) * 2021-06-01 2021-09-17 中国石油大学(北京) Oil and gas station lightning accident early warning method, device, equipment and storage medium
CN117057172A (en) * 2023-10-12 2023-11-14 宁波麦思捷科技有限公司武汉分公司 Method and system for monitoring electric field and magnetic field during lightning

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014173317A1 (en) * 2013-04-27 2014-10-30 国网四川省电力公司电力科学研究院 Power transmission line lightning electromagnetic transient dynamic simulation experimental system
CN106940745A (en) * 2017-03-13 2017-07-11 国家电网公司华中分部 Calculation of lightning over-voltage method based on full-wave electromagnetic transient Model
CN110765666A (en) * 2019-11-26 2020-02-07 国网山东省电力公司电力科学研究院 Simulation method for indirect breakdown fault of power transmission line caused by lightning stroke due to bifurcated lightning
CN110765731A (en) * 2019-11-25 2020-02-07 广州供电局有限公司 Lightning overvoltage simulation calculation method and device of current limiter and computer equipment
CN111079268A (en) * 2019-11-29 2020-04-28 中国电力科学研究院有限公司 Online simulation method and system for LCC direct current transmission electromechanical electromagnetic hybrid system
US11062078B1 (en) * 2018-09-21 2021-07-13 Global Energy Interconnection Research Institute Co., Ltd. Insulation coordination method and system for a series compensation apparatus, storage medium and electronic device

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014173317A1 (en) * 2013-04-27 2014-10-30 国网四川省电力公司电力科学研究院 Power transmission line lightning electromagnetic transient dynamic simulation experimental system
CN106940745A (en) * 2017-03-13 2017-07-11 国家电网公司华中分部 Calculation of lightning over-voltage method based on full-wave electromagnetic transient Model
US11062078B1 (en) * 2018-09-21 2021-07-13 Global Energy Interconnection Research Institute Co., Ltd. Insulation coordination method and system for a series compensation apparatus, storage medium and electronic device
CN110765731A (en) * 2019-11-25 2020-02-07 广州供电局有限公司 Lightning overvoltage simulation calculation method and device of current limiter and computer equipment
CN110765666A (en) * 2019-11-26 2020-02-07 国网山东省电力公司电力科学研究院 Simulation method for indirect breakdown fault of power transmission line caused by lightning stroke due to bifurcated lightning
CN111079268A (en) * 2019-11-29 2020-04-28 中国电力科学研究院有限公司 Online simulation method and system for LCC direct current transmission electromechanical electromagnetic hybrid system

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JOÃO CLAVIO SALARI: "A Methodology for Electromagnetic Transients Calculation—An Application for the Calculation of Lightning Propagation in Transmission Lines", IEEE TRANSACTIONS ON POWER DELIVERY, vol. 22, no. 1 *
何俊佳: "南方电网超／特高压输电线路防雷性能评估技术研究进展", 南方电网技术, vol. 10, no. 9 *
贾磊: "计算输电线路雷电感应过电压的三维电磁场模型", 高电压与绝缘技术 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112731076A (en) * 2020-12-17 2021-04-30 南方电网科学研究院有限责任公司 Early warning method and device based on insulation flashover rate under overvoltage and storage medium
CN112731076B (en) * 2020-12-17 2022-03-15 南方电网科学研究院有限责任公司 Early warning method and device based on insulation flashover rate under overvoltage and storage medium
CN112986659A (en) * 2021-02-01 2021-06-18 中国民航大学 Method for analyzing voltage drop in composite material airplane ground return network
CN113075495A (en) * 2021-02-24 2021-07-06 南方电网科学研究院有限责任公司 Circuit parameter adjusting method, device, equipment and medium of PEEC model
CN113191103A (en) * 2021-03-05 2021-07-30 南方电网科学研究院有限责任公司 Method and device for calculating lightning electromagnetic transient of distribution line
CN113406726A (en) * 2021-06-01 2021-09-17 中国石油大学(北京) Oil and gas station lightning accident early warning method, device, equipment and storage medium
CN117057172A (en) * 2023-10-12 2023-11-14 宁波麦思捷科技有限公司武汉分公司 Method and system for monitoring electric field and magnetic field during lightning
CN117057172B (en) * 2023-10-12 2023-12-29 宁波麦思捷科技有限公司武汉分公司 Method and system for monitoring electric field and magnetic field during lightning

Also Published As

Publication number Publication date
CN112052555B (en) 2024-02-23

Similar Documents

Publication Publication Date Title
CN112052555A (en) Simulation method and device for lightning electromagnetic transient model of power transmission line
KR100709616B1 (en) Line to ground fault location method for under ground cable system
Patel et al. Accurate impedance calculation for underground and submarine power cables using MoM-SO and a multilayer ground model
Tatematsu Development of a surge simulation code VSTL REV based on the 3D FDTD method
CN112711865A (en) Power transmission line inductive electric simulation refined modeling method
CN114298433A (en) Site selection method, device, equipment and medium for 5G base station in transformer substation
CN111598123A (en) Power distribution network line vectorization method and device based on neural network
JP3980811B2 (en) Electromagnetic field intensity calculating device, electromagnetic field intensity calculating method, and computer-readable recording medium storing a program for causing a computer to execute the method
Li et al. Distribution of currents in the lightning protective system of a residential building—Part II: Numerical modeling
Menshov et al. Surface–volume–surface electric field integral equation for magneto-quasi-static analysis of complex 3-D interconnects
CN107895078A (en) Macro model construction method of the strong electromagnetic pulse to long electric power transmission line coupling response
CN108459206B (en) Method for calculating spatial line distribution parameter matrix based on complex mirror image
Sajjadi et al. Lightning Response of Multi-Port Grids Buried in Dispersive Soils: An Approximation versus Full-wave Methods and Experiment
CN113191103A (en) Method and device for calculating lightning electromagnetic transient of distribution line
Sarajcev et al. Interfacing harmonic electromagnetic models of grounding systems with the EMTP-ATP software package
CN113075495A (en) Circuit parameter adjusting method, device, equipment and medium of PEEC model
Ng et al. Simulation of transmission line bend using a non-uniform transmission line model based on scattering theory
CN110909481A (en) Complex cable bundle distribution parameter modeling method based on moment method
Ma et al. Predicting Radiated Emissions from a Complex Transportation System Wiring Harness
Guan et al. A CFIE-based electromagnetic solver for composite objects
CN111859713B (en) Indoor substation power frequency electric field optimizing device
CN111756026A (en) Transformer direct-current magnetic bias suppression method based on multiple soil models
CN115358131B (en) Insulator design method, device, storage medium and system
CN111753391B (en) Power transmission line design method and system based on maximum field intensity of surface of split conductor
CN115828819B (en) High-frequency coupling current calculation method and system for high-precision transmission line

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