CN110287528B - Electromechanical-electromagnetic transient hybrid simulation power balancing method and device and storage medium - Google Patents

Electromechanical-electromagnetic transient hybrid simulation power balancing method and device and storage medium Download PDF

Info

Publication number
CN110287528B
CN110287528B CN201910431866.0A CN201910431866A CN110287528B CN 110287528 B CN110287528 B CN 110287528B CN 201910431866 A CN201910431866 A CN 201910431866A CN 110287528 B CN110287528 B CN 110287528B
Authority
CN
China
Prior art keywords
electromechanical
transient system
electromagnetic transient
physical quantity
electromagnetic
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
Application number
CN201910431866.0A
Other languages
Chinese (zh)
Other versions
CN110287528A (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.)
CSG Electric Power Research Institute
Guangxi Power Grid Co Ltd
Original Assignee
CSG Electric Power Research Institute
Guangxi 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 CSG Electric Power Research Institute, Guangxi Power Grid Co Ltd filed Critical CSG Electric Power Research Institute
Priority to CN201910431866.0A priority Critical patent/CN110287528B/en
Publication of CN110287528A publication Critical patent/CN110287528A/en
Application granted granted Critical
Publication of CN110287528B publication Critical patent/CN110287528B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Abstract

The invention discloses a power balance method, a device and a storage medium for electromechanical-electromagnetic transient hybrid simulation, which can correct an injection power value (or an injection current value) measured by an electromagnetic transient system in electromechanical-electromagnetic transient hybrid simulation, prevent power impact (or current impact) caused by the inconsistency of models of an electromechanical transient simulation program and an electromagnetic transient simulation program to the simulation, provide a stable initialization state for the electromechanical-electromagnetic transient hybrid simulation, and have simple steps and higher efficiency.

Description

Electromechanical-electromagnetic transient hybrid simulation power balancing method and device and storage medium
Technical Field
The invention relates to the technical field of transient stability research and simulation of a power system, in particular to an electromechanical-electromagnetic transient hybrid simulation power balancing method, device and storage medium.
Background
With the increasing of the interconnection scale of the power system regions, the relevance between the regions is enhanced, and especially with the use of a large number of power electronic devices in recent years, dynamic processes reflecting different physical characteristics are mutually interwoven, so that a high requirement is provided for a simulation program of the power system in order to ensure the safe and stable operation of the power system. Although the traditional electromechanical transient simulation program is rapidly and widely applied to the calculation process, the electromagnetic transient process with rapid local response cannot be accurately simulated, and the traditional electromagnetic transient simulation program is influenced by a computer, is difficult to carry out whole-network electromagnetic simulation on a large-scale power system, and needs equivalent simplification. In order to make up for the above two disadvantages, a conventional method is to adopt electromechanical-electromagnetic transient hybrid simulation, select a plurality of buses as interface buses in an electric power system, and divide the electric power system into an electromechanical transient system and an electromagnetic transient system, wherein a part describing the dynamic response characteristics in detail is divided into the electromagnetic transient system, an electromagnetic transient model is used for performing electromagnetic transient simulation, and the rest is divided into the electromechanical transient system, and a quasi-steady state model is used for performing electromechanical transient simulation.
According to the present computational principle of electromechanical-electromagnetic transient hybrid simulation, the core idea is to simulate an electromagnetic transient system by an electromagnetic transient simulation program, and replace the corresponding part of the electromagnetic transient system by measuring the injection power (or injection current) at the interface bus, in the electromechanical transient simulation program, the electromagnetic transient system is equivalent to an injection current source Ig, the value of the current source is measured from the electromagnetic transient system, and each electromechanical step is measured once.
In the calculation of the actual electromechanical-electromagnetic transient hybrid simulation, the electromechanical transient system is generally executed first, an initial injection power (or an initial injection current) is set at a bus interface position, the electromagnetic transient system is initialized based on the initial state, after the initialization of the electromagnetic transient system is completed, the joint hybrid simulation is started, the equivalent injection power (or the equivalent injection current) at the electromagnetic side is measured in real time, and the initial injection power (or the initial injection current) at the electromechanical side is replaced. However, for the same power system, since the electromechanical side adopts a quasi-steady-state model and the electromagnetic side adopts a detailed model, the results of the two are necessarily slightly different, which means that even if the power system is in a state without fault, the value of the equivalent injection power (or the equivalent injection current) may not completely coincide with the value of the initial injection power (or the initial injection current), and there is a certain power difference between the equivalent injection power and the initial injection power (there is a certain current difference between the equivalent injection current and the initial injection current), so that at a moment when entering the hybrid simulation, a corresponding power difference impact (or current difference impact) may be brought to the power system.
In order to avoid the above problems, the value of the initial injection power (or the initial injection current) is generally continuously modified to be substantially consistent with the equivalent injection power (or the equivalent injection current), however, this method requires a user to modify the basic data, and since the result of the equivalent injection power (or the equivalent injection current) is affected by the initial injection power (or the initial injection current), the value of the equivalent injection power (or the equivalent injection current) needs to be obtained after one electromechanical-electromagnetic transient hybrid simulation is performed in a state where the initial injection power (or the initial injection current) is determined, and once the value of the initial injection power (or the initial injection current) is modified, the corresponding initialization program and the hybrid simulation program need to be re-executed to obtain a new value of the equivalent injection power (or the equivalent injection current), and in order to make the power difference (or the current difference) tend to 0, multiple iterative corrections need to be performed, which is tedious in procedure and inefficient.
Disclosure of Invention
The technical problem to be solved by the embodiments of the present invention is to provide a method, an apparatus, and a storage medium for balancing electromechanical-electromagnetic transient hybrid simulation power, which can correct an injection power value (or an injection current value) measured by an electromagnetic transient system in electromechanical-electromagnetic transient hybrid simulation, prevent power impact (or current impact) caused by model inconsistency between an electromechanical transient simulation program and an electromagnetic transient simulation program, and provide a stable initialization state for electromechanical-electromagnetic transient hybrid simulation, and have simple steps and high efficiency.
In order to solve the above technical problem, an embodiment of the present invention provides a power balancing method for electromechanical-electromagnetic transient hybrid simulation, including:
s11, determining the position of an interface bus of a power system, and dividing the power system into an electromechanical transient system and an electromagnetic transient system according to the position of the interface bus;
s12, obtaining initial injection physical quantity provided by the electromagnetic transient system to the electromechanical transient system according to load flow calculation;
s13, simulating the electromagnetic transient system, and when the electromagnetic transient system reaches a steady state, obtaining equivalent injection physical quantity provided by the electromagnetic transient system to the electromechanical transient system;
s14, obtaining a correction coefficient according to the initial injection physical quantity and the equivalent injection physical quantity;
s15, correcting the equivalent injection physical quantity according to the correction coefficient, and supplying the corrected equivalent injection physical quantity to the electromechanical transient system;
and S16, simulating the electromechanical transient system according to the corrected equivalent injection physical quantity, simulating the electromagnetic transient system after the preset electromechanical simulation time is reached, obtaining the equivalent injection physical quantity which is provided by the electromagnetic transient system to the electromechanical transient system at present after the preset electromagnetic simulation time is reached, and returning to S15 until the simulation process is finished.
Further, the injection physical quantity is injection power or injection current.
Further, the correction coefficient is calculated by the following formula:
Figure BDA0002068534620000031
or
Figure BDA0002068534620000032
Where k represents a correction factor, U represents a voltage at the location of the interface bus, S i Denotes the initial injection power, S g(1) Representing a first acquired equivalent injected power provided by the electromagnetic transient system to the electromechanical transient system; i is i Denotes the initial injection current, I g(1) Representing the equivalent injection current provided by the electromagnetic transient system to the electromechanical transient system acquired for the first time, real () representing the real part of the complex number in brackets.
Further, the corrected equivalent injection physical quantity is obtained by calculating according to the following formula:
Figure BDA0002068534620000033
wherein x is not less than 1,x and represents the number of times of obtaining the equivalent injection physical quantity provided by the electromagnetic transient system to the electromechanical transient system; s. the g(x) Representing the equivalent injection power, S ', provided by the electromagnetic transient system to the electromechanical transient system at the x-th acquisition' g(x) Representing the corrected equivalent injection power corresponding to the equivalent injection power provided by the electromagnetic transient system to the electromechanical transient system for the xth acquisition; i is g(x) Represents an equivalent injection current, l ', provided by the electromagnetic transient system to the electromechanical transient system at the x-th acquisition' g(x) Representing the corrected equivalent injection current corresponding to the equivalent injection current provided by the electromagnetic transient system to the electromechanical transient system acquired at the x-th time, wherein image () represents the imaginary part of a complex number in a bracket, and j is an imaginary factor
Figure BDA0002068534620000041
In order to solve the above technical problem, an embodiment of the present invention further provides an electromechanical-electromagnetic transient hybrid simulation power balancing apparatus, including:
the system dividing module is used for determining the position of an interface bus of a power system and dividing the power system into an electromechanical transient system and an electromagnetic transient system according to the position of the interface bus;
the initial injection physical quantity acquisition module is used for acquiring initial injection physical quantity provided by the electromagnetic transient system to the electromechanical transient system according to load flow calculation;
the electromagnetic transient system simulation module is used for simulating the electromagnetic transient system, and when the electromagnetic transient system reaches a steady state, equivalent injection physical quantity provided by the electromagnetic transient system to the electromechanical transient system is obtained;
a correction coefficient obtaining module, configured to obtain a correction coefficient according to the initial injection physical quantity and the equivalent injection physical quantity;
the equivalent injection physical quantity correction module is used for correcting the equivalent injection physical quantity according to the correction coefficient and providing the corrected equivalent injection physical quantity to the electromechanical transient system; and (c) a second step of,
and the electromechanical transient system simulation module is used for simulating the electromechanical transient system according to the corrected equivalent injection physical quantity, calling the electromagnetic transient simulation module to simulate the electromagnetic transient system after the preset electromechanical simulation time is reached, obtaining the equivalent injection physical quantity which is provided by the electromagnetic transient system to the electromechanical transient system at present after the preset electromagnetic simulation time is reached, and calling the equivalent injection physical quantity correction module to continue processing until the simulation process is finished.
Further, the injection physical quantity is injection power or injection current.
Further, the correction coefficient is calculated by the following formula:
Figure BDA0002068534620000051
or
Figure BDA0002068534620000052
Where k represents a correction factor, U represents a voltage at the location of the interface bus, S i Denotes the initial injection power, S g(1) Representing a first acquired equivalent injection power provided by the electromagnetic transient system to the electromechanical transient system; i is i Denotes the initial injection current, I g(1) Representing the equivalent injection current provided by the electromagnetic transient system to the electromechanical transient system acquired for the first time, real () representing the real part of the complex number in brackets.
Further, the corrected equivalent injection physical quantity is calculated by the following formula:
Figure BDA0002068534620000053
wherein x is not less than 1,x represents the number of times of obtaining the equivalent injection physical quantity provided by the electromagnetic transient system to the electromechanical transient system; s g(x) Representing equivalent injection power, S ', provided by the electromagnetic transient system to the electromechanical transient system for the x-th acquisition' g(x) Representing the corrected equivalent injection power corresponding to the equivalent injection power provided by the electromagnetic transient system to the electromechanical transient system for the xth acquisition; i is g(x) Representing equivalent injection current, l ', provided by the electromagnetic transient system to the electromechanical transient system at the x-th acquisition' g(x) Representing the corrected equivalent injection current corresponding to the equivalent injection current provided by the electromagnetic transient system to the electromechanical transient system acquired at the x-th time, wherein image () represents the imaginary part of a complex number in a bracket, and j is an imaginary factor
Figure BDA0002068534620000054
An embodiment of the present invention further provides a computer-readable storage medium, where the computer-readable storage medium includes a stored computer program; wherein the computer program, when executed, controls an apparatus in which the computer readable storage medium is located to perform any one of the electromechanical-electromagnetic transient hybrid simulation power balancing methods described above.
An embodiment of the present invention further provides an electromechanical-electromagnetic transient hybrid simulation power balancing apparatus, which includes a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, and when the computer program is executed, the processor implements the electromechanical-electromagnetic transient hybrid simulation power balancing method according to any one of the above.
Compared with the prior art, the embodiment of the invention provides an electromechanical-electromagnetic transient hybrid simulation power balancing method, device and storage medium, which can correct the injection power value (or injection current value) measured by an electromagnetic transient system in electromechanical-electromagnetic transient hybrid simulation, prevent power impact (or current impact) brought to simulation due to model inconsistency of an electromechanical transient simulation program and an electromagnetic transient simulation program, provide a stable initialization state for electromechanical-electromagnetic transient hybrid simulation, and have simple steps and higher efficiency.
Drawings
FIG. 1 is a flow chart of a preferred embodiment of a power balancing method for an electromechanical-electromagnetic transient hybrid simulation provided by the present invention;
FIG. 2 is a block diagram of an electromechanical-electromagnetic transient hybrid simulation power balancing apparatus according to a preferred embodiment of the present invention;
fig. 3 is a structural block diagram of another preferred embodiment of the electromechanical-electromagnetic transient hybrid simulation power balancing apparatus provided in 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 any inventive step, are within the scope of the present invention.
An embodiment of the present invention provides a power balancing method for electromechanical-electromagnetic transient hybrid simulation, which is a flowchart of a preferred embodiment of the power balancing method for electromechanical-electromagnetic transient hybrid simulation provided by the present invention, as shown in fig. 1, and the method includes steps S11 to S16:
s11, determining the position of an interface bus of a power system, and dividing the power system into an electromechanical transient system and an electromagnetic transient system according to the position of the interface bus;
s12, obtaining initial injection physical quantity provided by the electromagnetic transient system to the electromechanical transient system according to load flow calculation;
s13, simulating the electromagnetic transient system, and when the electromagnetic transient system reaches a steady state, acquiring equivalent injection physical quantity provided by the electromagnetic transient system to the electromechanical transient system;
s14, obtaining a correction coefficient according to the initial injection physical quantity and the equivalent injection physical quantity;
s15, correcting the equivalent injection physical quantity according to the correction coefficient, and supplying the corrected equivalent injection physical quantity to the electromechanical transient system;
and S16, simulating the electromechanical transient system according to the corrected equivalent injection physical quantity, simulating the electromagnetic transient system after the preset electromechanical simulation time is reached, obtaining the equivalent injection physical quantity which is provided by the electromagnetic transient system to the electromechanical transient system at present after the preset electromagnetic simulation time is reached, and returning to the step S15 until the simulation process is finished.
Specifically, the position of an interface bus is determined in an electric power system, a system on one side of the interface bus is divided into an electromechanical transient system, the simulation is performed by using the existing electromechanical transient simulation program, a system on the other side of the interface bus is divided into an electromagnetic transient system, and the simulation is performed by using the existing electromagnetic transient simulation program; obtaining an initial injection physical quantity provided by the electromagnetic transient system to the electromechanical transient system according to a load flow calculation program of power system calculation software, recording the initial injection physical quantity, starting an electromagnetic transient simulation program, simulating the electromagnetic transient system, measuring the injection physical quantity at the position of an interface bus in real time, judging whether the electromagnetic transient system reaches a steady state or not according to the injection physical quantity obtained by real-time measurement (the steady state means that the injection physical quantity of the electromagnetic transient system at the position of the interface bus does not obviously fluctuate, for example, the electromagnetic transient system reaches the steady state when the proportion of the fluctuation range and the average value of the injection physical quantity is set to be within 1 percent), obtaining an equivalent injection physical quantity provided by the current electromagnetic transient system to the electromechanical transient system when the electromagnetic transient system reaches the steady state, recording the equivalent injection physical quantity obtained for the first time, calculating a correction coefficient according to the recorded initial injection physical quantity and the equivalent injection physical quantity obtained by the first simulation, correcting the equivalent injection physical quantity obtained for the first time according to the correction coefficient obtained by the calculation to obtain the equivalent injection physical quantity after the first correction, supplying the equivalent injection physical quantity after the first correction to the electromechanical transient system, starting an electromechanical transient simulation program, simulating the electromechanical transient system according to the equivalent injection physical quantity after the first correction, starting an electromagnetic transient simulation program after a preset electromechanical simulation time (for example, one electromechanical step length) is reached, simulating the electromagnetic transient system, and obtaining the equivalent injection physical quantity provided by the current electromagnetic transient system to the electromechanical transient system after the preset electromagnetic simulation time (for example, a plurality of electromagnetic step lengths and time corresponding to one electromechanical step length) is reached, and recording the equivalent injection physical quantity obtained for the second time, correcting the equivalent injection physical quantity obtained for the second time according to the correction coefficient to obtain the equivalent injection physical quantity corrected for the second time, supplying the equivalent injection physical quantity corrected for the second time to the electromechanical transient system, starting the electromechanical transient simulation program, simulating the electromechanical transient system according to the equivalent injection physical quantity corrected for the second time, starting the electromagnetic transient simulation program again after the preset electromechanical simulation time is reached, and performing the electromagnetic transient simulation program circularly according to the process until the whole simulation process is finished, for example, when the preset electromechanical-electromagnetic transient hybrid simulation time is reached, the whole simulation process is finished.
It should be noted that, in the initialization stage, the electromechanical transient simulation program is started first, the electromechanical transient system is simulated for a time of one electromechanical step according to the initial injected physical quantity, and then the electromagnetic transient simulation program is started to perform longer-time simulation on the electromagnetic transient system until the electromagnetic transient system reaches a steady state, so as to obtain the corresponding equivalent injected physical quantity (i.e., the equivalent injected physical quantity obtained for the first time).
In the electromechanical-electromagnetic transient hybrid simulation stage, firstly, an electromechanical transient simulation program executes the time simulation of one electromechanical step length according to the corrected equivalent injection physical quantity, then, the electromagnetic transient simulation program executes the time simulation of a plurality of electromagnetic step lengths (corresponding to one electromechanical step length), returns a new equivalent injection physical quantity, and executes the process circularly until the whole simulation process is finished.
As a modification of the above, the injection physical quantity is injection power or injection current.
As a modification of the above, the correction coefficient is calculated by the following formula:
Figure BDA0002068534620000081
or
Figure BDA0002068534620000082
Where k represents a correction factor, U represents a voltage at the location of the interface bus, S i Denotes the initial injection power, S g(1) Representing a first acquired equivalent injected power provided by the electromagnetic transient system to the electromechanical transient system; i is i Denotes the initial injection current, I g(1) Representing the equivalent injection current provided by the electromagnetic transient system to the electromechanical transient system acquired for the first time, real () representing the real part of the complex number in parentheses.
It should be noted that, when the injection power is used for injecting the physical quantity, the correction coefficient is expressed by the formula
Figure BDA0002068534620000083
Calculating to obtain the correction coefficient when the injected physical quantity adopts the injected current
Figure BDA0002068534620000091
And calculating, wherein real () represents the real part of a complex number in brackets, and a symbol denotes the conjugate of a corresponding physical quantity.
As an improvement of the above solution, the modified equivalent injection physical quantity is calculated by the following formula:
Figure BDA0002068534620000092
wherein x is not less than 1,x represents the number of times of obtaining the equivalent injection physical quantity provided by the electromagnetic transient system to the electromechanical transient system; s. the g(x) Representing equivalent injection power, S ', provided by the electromagnetic transient system to the electromechanical transient system for the x-th acquisition' g(x) Representing the corrected equivalent injection power corresponding to the equivalent injection power provided by the electromagnetic transient system to the electromechanical transient system for the xth acquisition; i is g(x) Representing equivalent injection current, l ', provided by the electromagnetic transient system to the electromechanical transient system at the x-th acquisition' g(x) Representing the corrected equivalent injection current corresponding to the equivalent injection current provided by the electromagnetic transient system to the electromechanical transient system acquired at the x-th time, wherein image () represents the imaginary part of a complex number in a bracket, and j is an imaginary factor
Figure BDA0002068534620000093
When the injection power is used as the injection physical quantity, the corrected equivalent injection power is expressed by the formula S' g(x) =real(S g(x) )/k+j·imag(S g(x) ) Calculating to obtain the equivalent injection current after correction when the injection physical quantity adopts the injection current
Figure BDA0002068534620000094
Obtained by calculation, wherein image () represents the imaginary part of the complex number in the brackets, and j is an imaginary factor
Figure BDA0002068534620000095
According to the electromechanical-electromagnetic transient hybrid simulation power balancing method provided by the embodiment of the invention, the value of the equivalent injection power or the equivalent injection current of the electromechanical transient system is corrected by the electromagnetic transient system in the electromechanical-electromagnetic transient hybrid simulation, so that power impact or current impact caused by the inconsistency of models of an electromechanical transient simulation program and an electromagnetic transient simulation program to the hybrid simulation can be prevented, a stable initialization state is provided for the electromechanical-electromagnetic transient hybrid simulation, the initial injection power or the initial injection current does not need to be modified, the initialization program only needs to be executed once, the simulation time is not prolonged, the steps are simple, and the efficiency is higher.
In addition, the method can also ensure that the value of the active part corresponding to the equivalent injection power or the equivalent injection current returned by the electromagnetic transient system for the first time is consistent with the value of the initial injection power or the initial injection current in a steady state, reduces the error between the return value and the initial value of each subsequent time in the case of no fault, effectively avoids the oscillation of the frequency of the power system caused by the mismatching of the active power, and further provides a steady initialization state for the electromechanical-electromagnetic transient hybrid simulation.
The embodiment of the present invention further provides an electromechanical-electromagnetic transient hybrid simulation power balancing apparatus, which is capable of implementing all processes of the electromechanical-electromagnetic transient hybrid simulation power balancing method described in any one of the above embodiments, and the functions and implemented technical effects of each module in the apparatus are respectively the same as those of the electromechanical-electromagnetic transient hybrid simulation power balancing method described in the above embodiment, and are not described herein again.
Referring to fig. 2, it is a block diagram of a preferred embodiment of an electromechanical-electromagnetic transient hybrid simulation power balancing apparatus provided in the present invention, the apparatus includes:
the system dividing module 11 is configured to determine a position of an interface bus of a power system, and divide the power system into an electromechanical transient system and an electromagnetic transient system according to the position of the interface bus;
the initial injection physical quantity acquisition module 12 is configured to obtain an initial injection physical quantity provided by the electromagnetic transient system to the electromechanical transient system according to load flow calculation;
the electromagnetic transient system simulation module 13 is configured to simulate the electromagnetic transient system, and when the electromagnetic transient system reaches a steady state, obtain an equivalent injection physical quantity provided by the electromagnetic transient system to the electromechanical transient system;
a correction coefficient obtaining module 14, configured to obtain a correction coefficient according to the initial injection physical quantity and the equivalent injection physical quantity;
the equivalent injection physical quantity correction module 15 is configured to correct the equivalent injection physical quantity according to the correction coefficient, and provide the corrected equivalent injection physical quantity to the electromechanical transient system; and the number of the first and second groups,
and the electromechanical transient system simulation module 16 is configured to simulate the electromechanical transient system according to the corrected equivalent injected physical quantity, call the electromagnetic transient simulation module 13 to simulate the electromagnetic transient system after a preset electromechanical simulation time is reached, obtain the equivalent injected physical quantity currently provided by the electromagnetic transient system to the electromechanical transient system after the preset electromagnetic simulation time is reached, and call the equivalent injected physical quantity correction module 15 to continue processing until the simulation process is finished.
Preferably, the injected physical quantity is injected power or injected current.
Preferably, the correction coefficient is calculated by the following formula:
Figure BDA0002068534620000111
or
Figure BDA0002068534620000112
Where k represents a correction factor, U represents a voltage at the location of the interface bus, S i Denotes the initial injection power, S g(1) Representing a first acquired equivalent injection power provided by the electromagnetic transient system to the electromechanical transient system; I.C. A i Denotes the initial injection current, I g(1) Representing the equivalent injection current provided by the electromagnetic transient system to the electromechanical transient system acquired for the first time, real () representing the real part of the complex number in parentheses.
Preferably, the corrected equivalent injection physical quantity is obtained by calculating according to the following formula:
Figure BDA0002068534620000113
wherein x is not less than 1,x represents the number of times of obtaining the equivalent injection physical quantity provided by the electromagnetic transient system to the electromechanical transient system; s g(x) Representing equivalent injection power, S ', provided by the electromagnetic transient system to the electromechanical transient system for the x-th acquisition' g(x) Representing the corrected equivalent injection power corresponding to the equivalent injection power provided by the electromagnetic transient system to the electromechanical transient system for the xth acquisition; I.C. A g(x) Represents an equivalent injection current, l ', provided by the electromagnetic transient system to the electromechanical transient system at the x-th acquisition' g(x) Representing the corrected equivalent injection current corresponding to the equivalent injection current provided by the electromagnetic transient system to the electromechanical transient system acquired at the x-th time, wherein image () represents the imaginary part of a complex number in a bracket, and j is an imaginary factor
Figure BDA0002068534620000114
An embodiment of the present invention further provides a computer-readable storage medium, where the computer-readable storage medium includes a stored computer program; wherein the computer program, when running, controls the apparatus on which the computer readable storage medium is located to perform the electromechanical-electromagnetic transient hybrid simulation power balancing method according to any of the above embodiments.
An embodiment of the present invention further provides an electromechanical-electromagnetic transient hybrid simulation power balancing apparatus, which is shown in fig. 3 and is a block diagram of another preferred embodiment of the electromechanical-electromagnetic transient hybrid simulation power balancing apparatus provided by the present invention, the apparatus 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, and the processor 10, when executing the computer program, implements the electromechanical-electromagnetic transient hybrid simulation power balancing method according to any of the above embodiments.
Preferably, the computer program can be divided into one or more modules/units (e.g. computer program 1, computer program 2,) which are stored in the memory 20 and executed by the processor 10 to accomplish the present invention. The one or more modules/units may be a series of computer program instruction segments capable of performing specific functions, which are used to describe the execution of the computer program in the apparatus.
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, a discrete Gate or transistor logic device, a discrete hardware component, etc., the general purpose Processor may be a microprocessor, or the Processor 10 may be any conventional Processor, the Processor 10 is a control center of the apparatus, and various interfaces and lines are used to connect various parts of the apparatus.
The memory 20 mainly includes a program storage area that may store an operating system, an application program required for at least one function, and the like, and a data storage area that may store related data and the like. In addition, the memory 20 may be a high speed random access memory, may also be a non-volatile memory, such as a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash Card (Flash Card), and the like, or the memory 20 may also be other volatile solid state memory devices.
It should be noted that the above-mentioned apparatus may include, but is not limited to, a processor and a memory, and those skilled in the art will understand that the structural block diagram in fig. 3 is only an example of the above-mentioned apparatus, and does not constitute a limitation of the apparatus, and may include more or less components than those shown, or combine some components, or different components.
To sum up, the electromechanical-electromagnetic transient hybrid simulation power balancing method, the electromechanical-electromagnetic transient hybrid simulation power balancing device, and the computer-readable storage medium provided by the embodiments of the present invention modify the value of the equivalent injection power or the equivalent injection current of the electromechanical transient system by the electromagnetic transient system in the electromechanical-electromagnetic transient hybrid simulation, so as to prevent power impact or current impact on the hybrid simulation due to model inconsistency between the electromechanical transient simulation program and the electromagnetic transient simulation program, provide a stable initialization state for the electromechanical-electromagnetic transient hybrid simulation, and do not need to modify the initial injection power or the initial injection current, the initialization program only needs to be executed once, the simulation time is not prolonged, and the steps are simple and the efficiency is high.
The above description is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, several modifications and variations can be made without departing from the technical principle of the present invention, and these modifications and variations should also be regarded as the protection scope of the present invention.

Claims (8)

1. An electromechanical-electromagnetic transient hybrid simulation power balancing method is characterized by comprising the following steps:
s11, determining the position of an interface bus of a power system, and dividing the power system into an electromechanical transient system and an electromagnetic transient system according to the position of the interface bus;
s12, obtaining initial injection physical quantity provided by the electromagnetic transient system to the electromechanical transient system according to load flow calculation;
s13, simulating the electromagnetic transient system, and when the electromagnetic transient system reaches a steady state, obtaining equivalent injection physical quantity provided by the electromagnetic transient system to the electromechanical transient system;
s14, obtaining a correction coefficient according to the initial injection physical quantity and the equivalent injection physical quantity;
s15, correcting the equivalent injection physical quantity according to the correction coefficient, and supplying the corrected equivalent injection physical quantity to the electromechanical transient system;
s16, simulating the electromechanical transient system according to the corrected equivalent injection physical quantity, simulating the electromagnetic transient system after the preset electromechanical simulation time is reached, obtaining the equivalent injection physical quantity which is provided by the electromagnetic transient system to the electromechanical transient system at present after the preset electromagnetic simulation time is reached, and returning to S15 until the simulation process is finished;
the corrected equivalent injection physical quantity is obtained by calculating the following formula:
S′ g(x) =real(S g(x) )/k+j·imag(S g(x) ) Or
Figure FDA0003952097540000011
Wherein x is not less than 1,x represents the number of times of obtaining the equivalent injection physical quantity provided by the electromagnetic transient system to the electromechanical transient system; s g(x) Representing equivalent injection power, S ', provided by the electromagnetic transient system to the electromechanical transient system for the x-th acquisition' g(x) Representing the corrected equivalent injection power corresponding to the equivalent injection power provided by the electromagnetic transient system to the electromechanical transient system for the xth acquisition; i is g(x) Representing equivalent injection current, l ', provided by the electromagnetic transient system to the electromechanical transient system at the x-th acquisition' g(x) Representing the corrected equivalent injection current corresponding to the equivalent injection current provided by the electromagnetic transient system to the electromechanical transient system acquired at the x-th time, wherein image () represents the imaginary part of a complex number in a bracket, and j is an imaginary factor
Figure FDA0003952097540000012
2. The electromechanical-electromagnetic transient hybrid simulation power balancing method of claim 1, wherein the injected physical quantity is an injected power or an injected current.
3. The electromechanical-electromagnetic transient hybrid simulation power balancing method according to claim 1 or 2, wherein the correction factor is calculated by the following formula:
Figure FDA0003952097540000021
or
Figure FDA0003952097540000022
Where k represents a correction factor, U represents a voltage at the location of the interface bus, S i Denotes the initial injection power, S g(1) Representing a first acquired equivalent injected power provided by the electromagnetic transient system to the electromechanical transient system; i is i Denotes the initial injection current, I g(1) Representing the equivalent injection current provided by the electromagnetic transient system to the electromechanical transient system acquired for the first time, real () representing the real part of the complex number in brackets.
4. An electromechanical-electromagnetic transient hybrid simulation power balancing apparatus, comprising:
the system dividing module is used for determining the position of an interface bus of a power system and dividing the power system into an electromechanical transient system and an electromagnetic transient system according to the position of the interface bus;
the initial injection physical quantity acquisition module is used for acquiring initial injection physical quantity provided by the electromagnetic transient system to the electromechanical transient system according to load flow calculation;
the electromagnetic transient system simulation module is used for simulating the electromagnetic transient system, and when the electromagnetic transient system reaches a steady state, equivalent injection physical quantity provided by the electromagnetic transient system to the electromechanical transient system is obtained;
a correction coefficient obtaining module, configured to obtain a correction coefficient according to the initial injection physical quantity and the equivalent injection physical quantity;
the equivalent injection physical quantity correction module is used for correcting the equivalent injection physical quantity according to the correction coefficient and supplying the corrected equivalent injection physical quantity to the electromechanical transient system; and the number of the first and second groups,
the electromechanical transient system simulation module is used for simulating the electromechanical transient system according to the corrected equivalent injected physical quantity, calling the electromagnetic transient simulation module to simulate the electromagnetic transient system after the preset electromechanical simulation time is reached, obtaining the equivalent injected physical quantity which is provided by the electromagnetic transient system to the electromechanical transient system at present after the preset electromagnetic simulation time is reached, and calling the equivalent injected physical quantity correction module to continue processing until the simulation process is finished;
the corrected equivalent injection physical quantity is obtained by calculating the following formula:
S′ g(x) =real(S g(x) )/k+j·imag(S g(x) ) Or
Figure FDA0003952097540000031
Wherein x is not less than 1,x represents the number of times of obtaining the equivalent injection physical quantity provided by the electromagnetic transient system to the electromechanical transient system; s g(x) Representing equivalent injection power, S ', provided by the electromagnetic transient system to the electromechanical transient system for the x-th acquisition' g(x) Representing the corrected equivalent injection power corresponding to the equivalent injection power provided by the electromagnetic transient system to the electromechanical transient system for the xth acquisition; i is g(x) Represents an equivalent injection current, l ', provided by the electromagnetic transient system to the electromechanical transient system at the x-th acquisition' g(x) Said representing the x-th acquisitionThe electromagnetic transient system provides a corrected equivalent injection current corresponding to the equivalent injection current provided by the electromechanical transient system, image () represents the imaginary part of a complex number in a bracket, and j is an imaginary factor
Figure FDA0003952097540000032
5. The electro-mechanical-electro-magnetic transient hybrid simulation power balancing apparatus of claim 4, wherein the injected physical quantity is an injected power or an injected current.
6. The electro-mechanical-electromagnetic transient hybrid simulation power balancing apparatus according to claim 4 or 5, wherein the correction coefficient is calculated by the following formula:
Figure FDA0003952097540000033
or
Figure FDA0003952097540000034
Where k represents a correction factor, U represents a voltage at the location of the interface bus, S i Denotes the initial injection power, S g(1) Representing a first acquired equivalent injected power provided by the electromagnetic transient system to the electromechanical transient system; i is i Denotes the initial injection current, I g(1) Representing the equivalent injection current provided by the electromagnetic transient system to the electromechanical transient system acquired for the first time, real () representing the real part of the complex number in parentheses.
7. A computer-readable storage medium, characterized in that the computer-readable storage medium comprises a stored computer program; wherein the computer program, when executed, controls an apparatus in which the computer-readable storage medium is located to perform the electromechanical-electromagnetic transient hybrid simulation power balancing method of any one of claims 1 to 3.
8. An electro-mechanical-electromagnetic transient hybrid simulation power balancing apparatus, 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 electro-mechanical-electromagnetic transient hybrid simulation power balancing method according to any one of claims 1 to 3 when executing the computer program.
CN201910431866.0A 2019-05-22 2019-05-22 Electromechanical-electromagnetic transient hybrid simulation power balancing method and device and storage medium Active CN110287528B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910431866.0A CN110287528B (en) 2019-05-22 2019-05-22 Electromechanical-electromagnetic transient hybrid simulation power balancing method and device and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910431866.0A CN110287528B (en) 2019-05-22 2019-05-22 Electromechanical-electromagnetic transient hybrid simulation power balancing method and device and storage medium

Publications (2)

Publication Number Publication Date
CN110287528A CN110287528A (en) 2019-09-27
CN110287528B true CN110287528B (en) 2023-02-28

Family

ID=68002230

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910431866.0A Active CN110287528B (en) 2019-05-22 2019-05-22 Electromechanical-electromagnetic transient hybrid simulation power balancing method and device and storage medium

Country Status (1)

Country Link
CN (1) CN110287528B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101957872A (en) * 2010-07-09 2011-01-26 南方电网技术研究中心 Hybrid real-time simulation method for alternating current-direct current power system
CN103279590A (en) * 2013-04-22 2013-09-04 中国南方电网有限责任公司电网技术研究中心 Initial self-correction computation method of interface power in electrical power system hybrid real-time simulation
CN104778317A (en) * 2015-04-07 2015-07-15 中国南方电网有限责任公司电网技术研究中心 Electromechanical-electromagnetic hybrid simulation method based on socket technology for electric power system
CN107622152A (en) * 2017-09-01 2018-01-23 南方电网科学研究院有限责任公司 A kind of electromechanical electromagnetic transient hybrid simulation method of bulk power grid
CN107834541A (en) * 2017-10-19 2018-03-23 中国电力科学研究院有限公司 A kind of half wavelength line electromechanical transient simulation method and system
CN107947155A (en) * 2017-11-14 2018-04-20 南方电网科学研究院有限责任公司 A kind of hybrid simulation method based on electromechanical transient and electro-magnetic transient
CN109245317A (en) * 2017-07-11 2019-01-18 中国电力科学研究院 A kind of electromechanical transient simulation system and method for battery energy storage system
CN109783895A (en) * 2018-12-27 2019-05-21 中国电力科学研究院有限公司 A kind of electric system hybrid simulation method and system containing distributed generation resource

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101957872A (en) * 2010-07-09 2011-01-26 南方电网技术研究中心 Hybrid real-time simulation method for alternating current-direct current power system
CN103279590A (en) * 2013-04-22 2013-09-04 中国南方电网有限责任公司电网技术研究中心 Initial self-correction computation method of interface power in electrical power system hybrid real-time simulation
CN104778317A (en) * 2015-04-07 2015-07-15 中国南方电网有限责任公司电网技术研究中心 Electromechanical-electromagnetic hybrid simulation method based on socket technology for electric power system
CN109245317A (en) * 2017-07-11 2019-01-18 中国电力科学研究院 A kind of electromechanical transient simulation system and method for battery energy storage system
CN107622152A (en) * 2017-09-01 2018-01-23 南方电网科学研究院有限责任公司 A kind of electromechanical electromagnetic transient hybrid simulation method of bulk power grid
CN107834541A (en) * 2017-10-19 2018-03-23 中国电力科学研究院有限公司 A kind of half wavelength line electromechanical transient simulation method and system
CN107947155A (en) * 2017-11-14 2018-04-20 南方电网科学研究院有限责任公司 A kind of hybrid simulation method based on electromechanical transient and electro-magnetic transient
CN109783895A (en) * 2018-12-27 2019-05-21 中国电力科学研究院有限公司 A kind of electric system hybrid simulation method and system containing distributed generation resource

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
宁夏电网多直流送端暂态过电压机电-电磁混合仿真研究;顾雨嘉,曾雪洋,田蓓;《宁夏电力》;20190131;论文第1-4+15页 *
混合仿真电流源等值误差机理分析及改进;肖湘宁,杨洋,陶顺,陈鹏伟;《电力建设》;20161231;论文第38-42页. *

Also Published As

Publication number Publication date
CN110287528A (en) 2019-09-27

Similar Documents

Publication Publication Date Title
CN109342989B (en) Error analysis method and device for electric energy meter
CN108875013B (en) Method and device for processing map data
CN108388509B (en) Software testing method, computer readable storage medium and terminal equipment
CN114626324B (en) FPGA circuit post-simulation verification method and device, electronic equipment and storage medium
CN112966465B (en) Circuit simulation method and device
CN111144082B (en) Multi-table iterative calculation method, device, electronic equipment and storage medium
CN111368259B (en) Distribution network voltage data missing filling method and related device
WO2023202355A1 (en) Soil body state data calculation method and device based on boundary surface plasticity model
CN110516334B (en) Convolution calculation simulation test method and device based on hardware environment and related equipment
CN110287528B (en) Electromechanical-electromagnetic transient hybrid simulation power balancing method and device and storage medium
CN116544999B (en) Method, device, equipment and medium for determining installed capacity of energy storage system
CN112487588A (en) Wind power plant reactive voltage regulation function test method, terminal equipment and storage medium
CN112711739A (en) Data processing method and device, server and storage medium
US20230055464A1 (en) Drc test pattern generation method and apparatus, electronic device, and storage medium
CN115859877A (en) Information synchronization method, device, equipment and storage medium in simulation test
CN110632909B (en) Method and equipment for detecting reasonability of back-checking MAP data
US8527923B1 (en) System, method, and computer program product for hierarchical formal hardware verification of floating-point division and/or square root algorithmic designs using automatic sequential equivalence checking
CN114564336A (en) Data consistency checking method, device, equipment and storage medium
CN104570757B (en) Method and device for operating an integrated control assembly
CN106815136B (en) Unit testing method and device
CN110580517A (en) Feature extraction method and device based on stacked self-encoder and terminal equipment
JP4135210B2 (en) Circuit design verification device
CN113569423B (en) Micro-grid joint fault positioning method, terminal and storage medium
CN116937613A (en) Data processing method and device based on reactive compensation system packaging model
US20070277133A1 (en) Model correspondence method and device

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