CN112052555B - 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 PDFInfo
- Publication number
- CN112052555B CN112052555B CN202010728774.1A CN202010728774A CN112052555B CN 112052555 B CN112052555 B CN 112052555B CN 202010728774 A CN202010728774 A CN 202010728774A CN 112052555 B CN112052555 B CN 112052555B
- Authority
- CN
- China
- Prior art keywords
- segment
- node
- model
- equivalent circuit
- nth
- 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.)
- Active
Links
- 230000005540 biological transmission Effects 0.000 title claims abstract description 52
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000004088 simulation Methods 0.000 title claims abstract description 52
- 230000001052 transient effect Effects 0.000 title claims abstract description 43
- 238000004613 tight binding model Methods 0.000 claims abstract description 16
- 230000006698 induction Effects 0.000 claims abstract description 11
- 230000005672 electromagnetic field Effects 0.000 claims abstract description 10
- 239000004020 conductor Substances 0.000 claims description 40
- 238000004590 computer program Methods 0.000 claims description 20
- 230000006870 function Effects 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 16
- 238000010276 construction Methods 0.000 claims description 7
- 230000004069 differentiation Effects 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 8
- 230000008878 coupling Effects 0.000 description 5
- 238000010168 coupling process Methods 0.000 description 5
- 238000005859 coupling reaction Methods 0.000 description 5
- 230000006872 improvement Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000035699 permeability Effects 0.000 description 3
- 206010042255 Struck by lightning Diseases 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 235000021384 green leafy vegetables Nutrition 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
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
Abstract
The invention discloses a simulation method and a device of 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 partial unit equivalent circuit model is an equivalent circuit model obtained by correcting according to a partial unit equivalent circuit formula and by utilizing a green function; constructing an overhead line model of the power transmission network according to the partial unit equivalent circuit model; and carrying out lightning induction overvoltage simulation on the overhead line model to obtain lightning-resistant horizontal and spatial electromagnetic field distribution conditions of the overhead line model. According to the invention, the influence of the lossy ground is simulated by using the green function so as 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 more accurately simulated.
Description
Technical Field
The invention relates to the technical field of lightning electromagnetic transient simulation of a power transmission line, in particular to a simulation method and device of a lightning electromagnetic transient model of the power transmission line.
Background
Lightning is one of the most dominant external factors responsible for distribution line faults. Conductors such as wires, ground wires, towers, etc. in a distribution line system can be considered as a multi-conductor system where lightning transients are a long standing problem. The multi-conductor system may be run in any direction and interconnected to form a three-dimensional wire structure on the lossy ground. When the system is struck by lightning, lightning current is generated in the conductor. In order to effectively protect against the effects of lightning strikes, it is necessary to build an appropriate model to model and evaluate lightning strike transients in the structure.
The partial cell equivalent circuit (PEEC) method is a method for modeling three-dimensional interconnect fine line structures. It is derived from the mixed potential integration equation and converts the electromagnetic problem into the circuit domain. In recent years, this method has been applied to handle transient currents and voltages in various structures of lines of arbitrary direction, but significant errors may occur since the influence of the lossy ground is not considered.
Disclosure of Invention
The embodiment of the invention aims to provide a simulation method and a device for a lightning electromagnetic transient model of a power transmission line, which utilize green's function to simulate the influence of a lossy ground so as to correct the existing PEEC method, thereby obtaining an improved partial unit equivalent circuit model and further being capable of more accurately simulating the lightning influence of the power transmission line.
In order to achieve the above purpose, 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 partial unit equivalent circuit model is an equivalent circuit model obtained by correcting according to a partial unit equivalent circuit formula and by utilizing a green function;
constructing an overhead line model of the power transmission network according to the partial unit equivalent circuit model;
and carrying out lightning induction overvoltage simulation on the overhead line model to obtain lightning-resistant horizontal and spatial electromagnetic field distribution conditions of the overhead line model.
Preferably, the node voltage of the partial unit equivalent circuit model is V n ＝∑ k p nk q kt +∑ t jωc nt I t And the node capacitance current isWherein V is n Node voltage segmented for nth node, I c,n Node capacitance current segmented for nth node, p nk Potential coefficient of conductor segment for nth node segment to kth node segment, q nk Charge amount of conductor segment for nth node segment to kth node segment, c nt Is 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, I t For the current of the source section t, p nn The potential coefficients for itself are segmented for the nth node, k+.n.
Preferably, the potential coefficient p of the conductor segment from the nth node segment to the kth node segment nk The calculation formula of (2) isWherein s is n Sum s k Respectively the sectional areas corresponding to the nth node segment and the kth node segment, l n And l k Respectively the lengths corresponding to the nth node segment and the kth node segment, K φ (r, r ') is a scalar green' function corresponding to a conductor surface node r=r ', r is a distance from an nth node segment to a preset observation point, and r' is a distance from the kth node segment to the observation point.
Preferably, the correction coefficient c nt The calculation formula of (2) isWherein s is n Sum s t Respectively the sectional areas of the nth node segment and the source segment t, l n And l t Respectively the length corresponding to the nth node segment and the source segment t, and P (r, r') is all P nk A composed potential coefficient matrix->Is the unit vector in the z direction of the vertical field, l' t For the length l of the source segment t t Corresponding differentiation.
The invention further provides a simulation device of a lightning electromagnetic transient model of a power transmission line, which comprises:
the model acquisition module is used for acquiring a preset partial unit equivalent circuit model; the partial unit equivalent circuit model is an equivalent circuit model obtained by correcting according to a partial unit equivalent circuit formula and by utilizing a green function;
the line model construction module is used for constructing an overhead line model of the power transmission network according to the partial unit equivalent circuit model;
and the simulation module is used for carrying out lightning induction overvoltage simulation on the overhead line model to obtain lightning-resistant horizontal and spatial electromagnetic field distribution conditions of the overhead line model.
A further embodiment of the invention correspondingly provides a device for simulating a lightning electromagnetic transient model of a power transmission line, comprising 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 simulating method of the lightning electromagnetic transient model of the power transmission line according to any one of the above when executing the computer program.
Compared with the prior art, the simulation method and the simulation device for the lightning electromagnetic transient model of the power transmission line provided by the embodiment of the invention utilize the green function to simulate the influence of the lossy land so as to correct the existing PEEC method, thereby obtaining an improved partial unit equivalent circuit model and further being capable of more accurately simulating the lightning influence of the power transmission line.
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 conductor segment structure in an air-to-ground medium according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of an equivalent circuit of a length of wire 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 power transmission line according to an embodiment of the invention.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Referring to fig. 1, a 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 is shown, where the method includes steps S1 to S3:
s1, acquiring a preset partial unit equivalent circuit model; the partial unit equivalent circuit model is an equivalent circuit model obtained by correcting according to a partial unit equivalent circuit formula and by utilizing a green function;
s2, constructing an overhead line model of the power transmission network according to the partial unit equivalent circuit model;
s3, carrying out lightning induction overvoltage simulation on the overhead line model to obtain lightning-resistant horizontal and spatial electromagnetic field distribution conditions of the overhead line model.
It should be noted that the partial cell equivalent circuit (Partial Element Equivalent Circuit, PEEC) method is a method for modeling three-dimensional interconnected thin line structures, which is derived from a mixed potential integral equation and converts electromagnetic problems into a circuit domain. This method can be applied to handle transient currents and voltages in various structures of lines in any direction, but without considering the effects of lossy ground, significant errors may occur. The invention corrects the mixed potential integral equation in the PEEC formula by using the green function so as to simulate the influence of the lossy earth, thereby obtaining an improved partial unit equivalent circuit model.
Specifically, a preset partial unit equivalent circuit model is obtained; the partial unit equivalent circuit model is an equivalent circuit model obtained by correcting according to a partial unit equivalent circuit formula and by utilizing a green function;
constructing an overhead line model of the power transmission network according to the partial equivalent circuit model;
and carrying out lightning induction overvoltage simulation on the overhead line model to obtain lightning-resistant horizontal and spatial electromagnetic field distribution conditions of the overhead line model.
In order to more clearly understand the construction process of the partial unit equivalent circuit model of the present invention, the following describes in detail the improvement process of the partial unit equivalent circuit model.
The magnetic vector potential a and scalar potential phi are often used to simplify the mathematical expression of maxwell's equations, satisfying the following equation at the conductor surface r point:where E is the electric field over the conductor surface and ω is the angular frequency.
Assuming air (ε) is present 0 ，μ 0 Sum sigma 0 ) And isotropic lossy ground surface (. Epsilon.) 1 ，μ 1 Sum sigma 1 ) Epsilon of the layered structure of (2) 0 And epsilon 1 Dielectric constant of air and dielectric constant of ground, mu, respectively 0 Sum mu 1 Magnetic permeability of air and magnetic permeability of ground, sigma, respectively 0 Sum sigma 1 The conductivity of air and the conductivity of ground are respectively shown in fig. 2, which is a schematic diagram of a conductor segment structure in an air-ground medium according to the embodiment of the present invention. For conductors located in air, both the magnetic vector potential a and the scalar potential Φ can be expressed as a greens function, as shown in the following equation:
where I (r ') and τ (r ') are the conductor current and line charge density on the conductor at r=r ', respectively. In the method, in the process of the invention,as a phasor green function of the magnetic vector a, there are different components to the current dipole in three orthogonal directions:
wherein g (r, r ') is a three-dimensional spatial green's function of the homogeneous medium,wherein R= |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 wire, the preset is generally performed manually, d is a calculation coefficient, and the calculation coefficient is the calculation coefficient>Is expressed asIn (1) the->For the unit vector of the current dipole projected on the xy plane or ground, +.>Is a unit vector in the x direction of the horizontal field, +.>Is a unit vector in the y-direction of the horizontal field, +.>Is the unit vector in the z direction of the vertical field, G xx Is->At->The vector component on the component, G yy Is->At the position ofThe vector component on the component, G zz Is->At->The vector component on the component, G zx Is->At->The vector component on the component, G zy Is->At->The vector component on the component, G tt Is->At->The vector component on the component, G zt Is->At->The vector component on the component, and G xx ＝G yy ＝G tt L 'is the segmented length of the source conductor and P (r, r') is all P nk Potential coefficient matrix, K φ Is a scalar greens function. For current dipoles oriented in layered media, there is no unique scalar green function, K in the present invention φ Is determined from g (r, r') and Sophine integral.
Dividing the conductor in the air into a plurality of segments, and selecting two segments: the field segment i and the source segment t, and then integrating the formula along the field segment yields the potential equation as follows:
wherein s is i 、l i And J i Respectively the cross-sectional area, the length and the current density of the field section i t 、l′ t 、I t Respectively the length, the length differential and the conductor current of the source section t n And phi m Node voltages, τ, for the nth node segment and the mth node segment, respectively, on field segment i k And l' k The linear charge density and the length differential of the kth node segment on field segment i are respectively, sigma is the conductivity, n noteq m, n noteq k.
From the definition of the average node potential (voltage) V, a circuit equation can be established as follows:
Z s,i I+∑ t jωL it I t ＝V n -V m
V n ＝∑ k p nk q kt +∑ t jωc nt I t
Z s,i the surface impedance of the conductor field section i is expressed asR a For an internal impedance of the conductor field segment i, for an electrical conductivity of σ, the permeability is μ 0 Circular conductor with radius a, +.>I is a higher order Bessel function of the first kind, I 0 And I 1 Solutions to the 0 th and 1 st orders, respectively. I t For the current of the conductor source section t, parameter L it Inductance, p, for the segment from field segment i to source segment t nk Potential coefficients, q, for segments from node segment n to node segment k nk For the amount of charge of node segment n to node segment k, c nt Is a correction coefficient.
According to the above analysis, there are
Taking the time derivative of the potential to obtain the capacitance current on the node in the frequency domain as follows:
wherein q is n For the charge quantity of node n, I c,k The node capacitance current for the kth node.
Node voltage V n And branch current I i Two groups of unknowns in the equivalent circuit of the partial unit can be solved by using a conventional circuit analysis tool.
According to the simulation method of the lightning electromagnetic transient model of the power transmission line, provided by the embodiment 1, the influence of the lossy land is simulated by using the green function, so that the existing PEEC method is corrected, 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 partial unit equivalent circuit model is V n ＝∑ k p nk q kt +∑ t jωc nt I t And the node capacitance current isWherein V is n Node voltage segmented for nth node, I c,n Node capacitance current segmented for nth node, p nk Potential coefficient of conductor segment for nth node segment to kth node segment, q nk Charge amount of conductor segment for nth node segment to kth node segment, c nt Is 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, I t For the current of the source section t, p nn The potential coefficients for itself are segmented for the nth node, k+.n.
Referring to fig. 3, an equivalent circuit diagram of a section of a wire obtained by the method according to the present invention is provided in this embodiment of the present invention, that is, a partial unit equivalent circuit obtained by the improvement of the present invention is shown in fig. 3. Inductance L it Comprising coupling elements from horizontal (source) sections to vertical (field) sections, as is conventionalIs not present in the definition of inductance of (c). The introduction of this coupling element is purely due to the presence of a lossy ground, a correction term for mutual inductance coupling. However, due to the asymmetric matrixThe 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 cell equivalent circuit model is V n ＝∑ k p nk q kt +∑ t jωc nt I t And the node capacitance current isWherein V is n Node voltage segmented for nth node, I c,n Node capacitance current segmented for nth node, p nk Potential coefficient of conductor segment for nth node segment to kth node segment, q nk Charge amount of conductor segment for nth node segment to kth node segment, c nt Is a preset correction coefficient, omega is angular frequency, a part of unit equivalent circuit model comprises a plurality of node segments, each node segment comprises a field segment and a source segment, I t For the current of the source section t, p nn The potential coefficients for itself are segmented for the nth node, k+.n.
As an improvement of the above scheme, the potential coefficient p of the conductor segment from the nth node segment to the kth node segment nk The calculation formula of (2) isWherein s is n Sum s k Respectively the sectional areas corresponding to the nth node segment and the kth node segment, l n And l k Respectively the lengths corresponding to the nth node segment and the kth node segment, K φ (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 a preset observation point, and r' is the distance from the kth node segment to the observation point.
Specifically, the potential coefficient p of the conductor segment from the nth node segment to the kth node segment nk The calculation formula of (2) isWherein s is n Sum s k Respectively the sectional areas corresponding to the nth node segment and the kth node segment, l n And l k Respectively the lengths corresponding to the nth node segment and the kth node segment, K φ (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 a preset observation point, and r' is the distance from the kth node segment to the observation point.
As an improvement of the above scheme, the correction coefficient c nt The calculation formula of (2) isWherein s is n Sum s t Respectively the sectional areas of the nth node segment and the source segment t, l n And l t Respectively the length corresponding to the nth node segment and the source segment t, and P (r, r') is all P nk A composed potential coefficient matrix->Is the unit vector in the z direction of the vertical field, l' t For the length l of the source segment t t Corresponding differentiation.
Specifically, correction coefficient c nt The calculation formula of (2) isWherein s is n Sum s t Respectively the sectional areas of the nth node segment and the source segment t, l n And l t The lengths corresponding to the nth node segment and the source segment t are respectively, and P (r, r') is all P nk A composed potential coefficient matrix->In the z direction of the vertical fieldUnit vector, l' t Length l of source segment t t Corresponding differentiation.
Referring to fig. 4, an apparatus for simulating a lightning electromagnetic transient model of a power transmission line according to an embodiment of the present invention is characterized in that the apparatus includes:
the model acquisition module 11 is used for acquiring a preset partial unit equivalent circuit model; the partial unit equivalent circuit model is an equivalent circuit model obtained by correcting according to a partial unit equivalent circuit formula and by utilizing a green function;
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 used for carrying out lightning induction overvoltage simulation on the overhead line model to obtain lightning-resistant level and space electromagnetic field distribution conditions of the overhead line model.
The simulation device for the lightning electromagnetic transient model of the power transmission line provided by the embodiment of the invention can realize all the processes of the simulation method for the lightning electromagnetic transient model of the power transmission line described in any embodiment, and the functions and the realized technical effects of each module and each unit in the device are respectively the same as those of the simulation method for the lightning electromagnetic transient model of the power transmission line described in the embodiment, and are not repeated here.
Referring to fig. 5, a schematic diagram of an apparatus for a simulation method using a lightning electromagnetic transient model of a power transmission line according to the embodiment of the present invention is provided, where the apparatus for a simulation method using a lightning electromagnetic transient model of a power transmission line 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 using a lightning electromagnetic transient model of a power transmission line according to any one of the embodiments.
By way of example, a computer program may be partitioned into one or more modules/units that are stored in the memory 20 and executed by the processor 10 to perform the present invention. One or more of the modules/units may be a series of computer program instruction segments capable of performing a specific function for describing the execution of the computer program in a method for simulating a lightning electromagnetic transient model of an electrical transmission line. For example, the computer program may be divided into a model acquisition module, a line model construction module and a simulation module, each module having the following specific functions:
the model acquisition module 11 is used for acquiring a preset partial unit equivalent circuit model; the partial unit equivalent circuit model is an equivalent circuit model obtained by correcting according to a partial unit equivalent circuit formula and by utilizing a green function;
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 used for carrying out lightning induction overvoltage simulation on the overhead line model to obtain lightning-resistant level and space electromagnetic field distribution conditions 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 lightning electromagnetic transient model of the power transmission line can comprise, but is not limited to, a processor and a memory. It will be appreciated by those skilled in the art that the schematic diagram 5 is merely an example of an apparatus using a simulation method of a lightning electromagnetic transient model of a power transmission line, and is not limited to the apparatus using the simulation method of a lightning electromagnetic transient model of a power transmission line, and may include more or less components than those illustrated, or may combine some components, or different components, for example, the apparatus using the simulation method of a lightning electromagnetic transient model of a power transmission line may further include an input/output device, a network access device, a bus, and so on.
The processor 10 may be a central processing unit (Central Processing Unit, CPU), but may also be other general purpose processors, digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), off-the-shelf programmable gate arrays (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, or the like. The general processor may be a microprocessor or the processor 10 may also be any conventional processor, etc., and the processor 10 is a control center of the apparatus using the simulation method of the lightning electromagnetic transient model of the transmission line, and connects various parts of the entire apparatus using the simulation method of the lightning electromagnetic transient model of the transmission line by using various interfaces and lines.
The memory 20 may be used to store the computer program and/or module, and the processor 10 implements the various functions of the apparatus using the transmission line lightning electromagnetic transient model simulation method by running or executing the computer program and/or module stored in the memory 20 and invoking data stored in the memory 20. The memory 20 may mainly include a storage program area that may store an operating system, application programs required for at least one function (such as a sound playing function, an image playing function, etc.), and a storage data area; the storage data area may store data (such as audio data, phonebook, etc.) created according to the use of the handset, etc. In addition, the memory 20 may include high-speed random access memory, and may also include nonvolatile memory, such as a hard disk, memory, plug-in hard disk, smart Media Card (SMC), secure Digital (SD) Card, flash Card (Flash Card), at least one disk storage device, flash memory device, or other volatile solid state storage device.
The module integrated by the device of the simulation method using 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 a separate product. Based on such understanding, the present invention may implement all or part of the flow of the method of the above embodiment, or may be implemented by a computer program to instruct related hardware, where the computer program may be stored in a computer readable storage medium, and when the computer program is executed by a processor, the steps of each method embodiment may be implemented. The computer program comprises computer program code, and the computer program code can be in a source code form, an object code form, an executable file or some intermediate form and the like. The computer readable medium may include: any entity or device capable of carrying computer program code, a recording medium, a U disk, a removable hard disk, a magnetic disk, an optical disk, a computer Memory, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), an electrical carrier signal, a telecommunications signal, a software distribution medium, and so forth. It should be noted that the content of the computer readable medium can be appropriately increased or decreased according to the requirements of the jurisdiction's jurisdiction and the patent practice, for example, in some jurisdictions, the computer readable medium does not include electrical carrier signals and telecommunication signals according to the jurisdiction and the patent practice.
The embodiment of the invention also provides a computer readable storage medium, which comprises a stored computer program, wherein when the computer program runs, equipment where the computer readable storage medium is located is controlled to execute the simulation method of the lightning electromagnetic transient model of the power transmission line.
In summary, the simulation method and the simulation device for the lightning electromagnetic transient model of the power transmission line provided by the embodiment of the invention firstly analyze a full-wave PEEC formula of lightning transient analysis of any directional conductor on a lossy ground, simulate the influence of the lossy foundation by using a parallel vector green function, and finally provide a new PEEC equivalent circuit by considering correction items caused by lossy grounding. Modeling is carried out through the improved PEEC equivalent circuit, a distribution network overhead line is established, calculation processing is carried out on the electromagnetic coupling effect of the lossy ground, electric field-circuit conversion is realized, modeling simulation is carried out on the lightning induction overvoltage of the overhead line, the lightning resistance level of the overhead line is tested, the space electromagnetic field is analyzed, the complex network is realized, the induction voltage calculation modeling is realized, and the distribution network construction and operation reference can be provided. The invention is helpful for 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 lossy ground influence 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 the soil electromagnetic coupling.
While the foregoing is directed to the preferred embodiments of the present invention, it will be appreciated by those skilled in the art that changes and modifications may be made without departing from the principles of the invention, such changes and modifications are also intended to be within the scope of the invention.
Claims (3)
1. The simulation method of the lightning electromagnetic transient model of the power transmission line is characterized by comprising the following steps of:
acquiring a preset partial unit equivalent circuit model; the partial unit equivalent circuit model is an equivalent circuit model obtained by correcting according to a partial unit equivalent circuit formula and by utilizing a green function;
constructing an overhead line model of the power transmission network according to the partial unit equivalent circuit model;
carrying out lightning induction overvoltage simulation on the overhead line model to obtain lightning-resistant horizontal and spatial electromagnetic field distribution conditions of the overhead line model;
wherein the node voltage of the partial unit equivalent circuit model is V n ＝∑ k p nk q kt +∑ t jωc nt I t And the node capacitance current isWherein V is n Node voltage segmented for nth node, I c,n Node capacitance current segmented for nth node, p nk Potential coefficient of conductor segment for nth node segment to kth node segment, q nk Charge amount of conductor segment for nth node segment to kth node segment, c nt Is 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, I t For the current of the source section t, p nn Potential coefficients of the n-th node segment pair are k not equal to n;
potential coefficient p of conductor segment from nth node segment to kth node segment nk The calculation formula of (2) isWherein s is n Sum s k Respectively the sectional areas corresponding to the nth node segment and the kth node segment, l n And l k Respectively the lengths corresponding to the nth node segment and the kth node segment, K φ (r, r ') is a scalar green' function corresponding to a 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;
the correction coefficient c nt The calculation formula of (2) isWherein s is n Sum s t Respectively the sectional areas of the nth node segment and the source segment t, l n And l t Respectively the length corresponding to the nth node segment and the source segment t, and P (r, r') is all P nk A composed potential coefficient matrix->Is the unit vector in the z direction of the vertical field, l t ' is the length l of the source segment t t Corresponding differentiation.
2. The simulation device of the lightning electromagnetic transient model of the power transmission line is characterized by comprising:
the model acquisition module is used for acquiring a preset partial unit equivalent circuit model; the partial unit equivalent circuit model is an equivalent circuit model obtained by correcting according to a partial unit equivalent circuit formula and by utilizing a green function;
the line model construction module is used for constructing an overhead line model of the power transmission network according to the partial unit equivalent circuit model;
the simulation module is used for carrying out lightning induction overvoltage simulation on the overhead line model to obtain lightning-resistant level and space electromagnetic field distribution conditions of the overhead line model;
wherein the node voltage of the partial unit equivalent circuit model is V n ＝∑ k p nk q kt +∑ t jωc nt I t And the node capacitance current isWherein V is n Node voltage segmented for nth node, I c,n Node capacitance current segmented for nth node, p nk Potential coefficient of conductor segment for nth node segment to kth node segment, q nk Charge amount of conductor segment for nth node segment to kth node segment, c nt Is 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, I t For the current of the source section t, p nn Potential coefficients of the n-th node segment pair are k not equal to n;
potential coefficient p of conductor segment from nth node segment to kth node segment nk The calculation formula of (2) isWherein s is n Sum s k Respectively the sectional areas corresponding to the nth node segment and the kth node segment, l n And l k Respectively the lengths corresponding to the nth node segment and the kth node segment, K φ (r, r ') is a scalar green' function corresponding to a 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;
the correction coefficient c nt The calculation formula of (2) isWherein s is n Sum s t Respectively the sectional areas of the nth node segment and the source segment t, l n And l t Respectively the length corresponding to the nth node segment and the source segment t, and P (r, r') is all P nk A composed potential coefficient matrix->Is the unit vector in the z direction of the vertical field, l t ' is the length l of the source segment t t Corresponding differentiation.
3. An apparatus for using a simulation method of a lightning electromagnetic transient model of an electrical transmission line, 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 lightning electromagnetic transient model of an electrical transmission line as claimed in claim 1 when the computer program is executed.
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 CN112052555A (en) | 2020-12-08 |
| CN112052555B true 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) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| CN113406726B (en) * | 2021-06-01 | 2022-09-27 | 中国石油大学(北京) | Oil and gas station lightning accident early warning method, device, equipment and storage medium |
| CN117057172B (en) * | 2023-10-12 | 2023-12-29 | 宁波麦思捷科技有限公司武汉分公司 | Method and system for monitoring electric field and magnetic field during lightning |
Citations (6)
| 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 |
| 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 |
| 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 |
-
2020
- 2020-07-23 CN CN202010728774.1A patent/CN112052555B/en active Active
Patent Citations (6)
| 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)
| Title |
|---|
| A Methodology for Electromagnetic Transients Calculation—An Application for the Calculation of Lightning Propagation in Transmission Lines;João Clavio Salari;IEEE TRANSACTIONS ON POWER DELIVERY;第22卷(第1期);全文 * |
| 南方电网超／特高压输电线路防雷性能评估技术研究进展;何俊佳;南方电网技术;第10卷(第9期);全文 * |
| 计算输电线路雷电感应过电压的三维电磁场模型;贾磊;高电压与绝缘技术;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112052555A (en) | 2020-12-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112052555B (en) | Simulation method and device for lightning electromagnetic transient model of power transmission line | |
| Kersting et al. | The application of Carson's equation to the steady-state analysis of distribution feeders | |
| 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 | |
| CN106649946A (en) | Simulation calculation method for power frequency phase parameters of power transmission line | |
| Maslowski et al. | Experimental investigation and numerical modeling of surge currents in lightning protection system of a small residential structure | |
| CN108459206B (en) | Method for calculating spatial line distribution parameter matrix based on complex mirror image | |
| CN109902388B (en) | Optimized wave number selection and application method for infinite long cylinder direct current electric field solution | |
| CN110927484A (en) | Alternating current transmission line radio interference three-dimensional distribution calculation method and system | |
| Valdes et al. | Time-domain single-source integral equations for analyzing scattering from homogeneous penetrable objects | |
| Nuricumbo‐Guillén et al. | Computation of transient voltage and current profiles along illuminated multiconductor lines by means of the numerical Laplace transform | |
| CN113075495A (en) | Circuit parameter adjusting method, device, equipment and medium of PEEC model | |
| CN113191103A (en) | Method and device for calculating lightning electromagnetic transient of distribution line | |
| Wiid et al. | Towards electromagnetic characterization of MeerKAT Telescope | |
| CN108170964B (en) | Method for calculating potential of grounding body of lightning pole tower | |
| De Conti et al. | Lightning transients on branched distribution lines considering frequency-dependent ground parameters | |
| McGahay | Conformal mapping solution for interdigital comb capacitors between ground planes | |
| CN112014815A (en) | Hybrid surface integral analysis method, device, equipment and storage medium | |
| Frei et al. | Calculation of low frequency EMC problems in large systems with a quasi-static approach | |
| Han et al. | A novel field-line-circuit hybrid algorithm for transient responses prediction of transmission lines based on FDTD method | |
| CN115358131B (en) | Insulator design method, device, storage medium and system | |
| CN111756026A (en) | Transformer direct-current magnetic bias suppression method based on multiple soil models | |
| De Silva et al. | On combining classical and numerical techniques for extracting the impedance and admittance of cables and overhead lines | |
| Tikhomirova | Electromagnetic transient modelling of grounding structures | |
| Avula et al. | A novel iterative method for approximating frequency response with equivalent pole/residues |
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 |