CN117057284A - Conversion data processing equivalent adjustment method and device for electromagnetic transient model - Google Patents
Conversion data processing equivalent adjustment method and device for electromagnetic transient model Download PDFInfo
- Publication number
- CN117057284A CN117057284A CN202310627934.7A CN202310627934A CN117057284A CN 117057284 A CN117057284 A CN 117057284A CN 202310627934 A CN202310627934 A CN 202310627934A CN 117057284 A CN117057284 A CN 117057284A
- Authority
- CN
- China
- Prior art keywords
- bus
- electromagnetic
- line
- transient model
- electromechanical
- 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.)
- Pending
Links
- 230000001052 transient effect Effects 0.000 title claims abstract description 105
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000012545 processing Methods 0.000 title claims abstract description 34
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 28
- 238000004088 simulation Methods 0.000 claims abstract description 36
- 230000001502 supplementing effect Effects 0.000 claims abstract description 18
- 238000004590 computer program Methods 0.000 claims description 22
- 238000004804 winding Methods 0.000 claims description 10
- 230000007812 deficiency Effects 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 4
- 230000000295 complement effect Effects 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 abstract description 9
- 230000005540 biological transmission Effects 0.000 abstract description 6
- 238000013461 design Methods 0.000 abstract description 6
- 230000007246 mechanism Effects 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 10
- 230000006870 function Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000003993 interaction Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000003238 somatosensory effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J5/00—Circuit arrangements for transfer of electric power between ac networks and dc networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Power Engineering (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The application provides a conversion data processing equivalence adjustment method and device of an electromagnetic transient model, wherein the method comprises the following steps: acquiring an electromagnetic transient model obtained by converting an electromechanical transient model; operating an electromagnetic transient model, and if the line is solved for errors and the zero sequence parameter of the read error line is 0, supplementing default parameters according to a supplementing formula; if the line solving is not wrong, comparing the busbar voltage results of the electromechanical transient model and the electromagnetic transient model, and perfecting the electromagnetic transient model according to the busbar voltage deviation value in the comparison result to obtain the target electromagnetic transient model. The application can quickly perfect electromagnetic transient model parameters without manually checking the parameters of each element. The application is not aimed at specific electromechanical software, can be suitable for data adjustment after different electromechanical software is converted into the same electromagnetic software, can be used for model data adjustment of scenes such as large power grid planning and design, and provides simulation support for limiting factors and mechanism analysis of power transmission and power reception capacity of a planning grid.
Description
Technical Field
The application relates to the technical field of software processing, in particular to a conversion data processing equivalence adjusting method and device of an electromagnetic transient model.
Background
With the sequential operation of extra-high voltage alternating current-direct current engineering, the regional power grids are more and more tightly connected. The power grid presents the phenomenon of strong and weak alternating current, the problem of alternating current-direct current and multi-direct current interaction influence of the transmitting and receiving ends is increasingly outstanding, and the running characteristic of the power grid is increasingly complex. Thus, modeling of multiple regions needs to be performed simultaneously for analysis of large grid characteristics. At present, the common analysis and simulation of the operation characteristics of the large power grid adopts electromechanical transient simulation software, but has the limitation of nonlinear characteristics of networks such as on-off of power electronic devices, has larger working condition errors in the simulation of asymmetrical faults and the like, and cannot accurately simulate the physical characteristics of the power grid containing large-scale direct current; and the electromagnetic transient simulation has high simulation precision, and the alternating current-direct current interaction characteristic and the multi-circuit direct current interaction characteristic can be more accurately researched. The large power grid analysis model is complex, the calculation scale is large, and model data filling is mainly carried out by adopting an electromechanical transient simulation software format. However, with the continuous development of large-scale electromagnetic transient simulation technology, an electromagnetic model of a regional power grid can be realized at present, and the existing electromagnetic transient simulation model has less data, so that an electromechanical model needs to be converted into an electromagnetic model.
The electromechanical simulation software and the electromagnetic simulation software at home and abroad are more in variety, the format difference of different software is larger, and the data quality requirement degree is different. Therefore, the existing model conversion method and model conversion program of the electromechanical-electromagnetic simulation software are often aimed at the appointed simulation software, and most verification examples select a small-sized power grid. The conversion model for large power grids is less comparative. The large power grid model is complex, the parameters are numerous, some element parameters default in the electromechanical simulation software, the parameters are supplemented through a built-in calculation formula, the calculation of the electromechanical simulation is not affected, but waveform errors can occur and even the simulation cannot be performed when the simulation is converted into the electromagnetic simulation software.
Although the model conversion methods of different electromechanical-electromagnetic simulation software are more literature, the method for adjusting the converted data is not blank. Therefore, a great need exists to provide a method for adjusting the equivalent of the conversion data processing of the large power grid, which has strong universality and practicability. To prevent the above-mentioned errors that may occur when converting complex large grid electromagnetic models.
Disclosure of Invention
In view of the above, the present application provides a method and apparatus for adjusting the equivalent of conversion data processing of an electromagnetic transient model, so as to solve at least one of the above-mentioned problems.
In order to achieve the above purpose, the present application adopts the following scheme:
according to a first aspect of the present application, there is provided a conversion data processing equivalence adjustment method of an electromagnetic transient model, the method comprising: acquiring an electromagnetic transient model obtained by converting an electromechanical transient model; and (3) operating the electromagnetic transient model, and if the line solving is in error and the zero sequence parameter of the read error line is 0, supplementing the default parameter according to the following supplementing formula: if the line positive sequence is: g i1 -jB i1 ＝1/(R i1 +jX i1 ) The method comprises the steps of carrying out a first treatment on the surface of the The line zero sequence is: g i0 -jB i0 ＝0.8G i1 -j0.313B i1 ，G ic0 ＝0.7G ic1 ，B ic0 ＝0.68B ic1 The method comprises the steps of carrying out a first treatment on the surface of the In the above, G i1 Is the positive sequence conductance of line i, B i1 Is the positive sequence susceptance of line i, R i1 Is the positive sequence resistance of line i, X i1 Is the positive sequence reactance of line i; g i0 Is the zero sequence conductance of line i, B i0 Is the zero sequence susceptance of line i, G i c 0 Is the zero sequence conductance of the ground to the line i, G i c 1 Is the ground-pair positive sequence conductance of line i, B i c 0 Is the ground-to-ground zero sequence susceptance of line i, B ic1 The positive sequence susceptance of the ground pair of the line i; if the line solving is not wrong, comparing the busbar voltage results of the electromechanical transient model and the electromagnetic transient model, and perfecting the electromagnetic transient model according to the busbar voltage deviation value in the comparison result to obtain a target electromagnetic transient model.
According to a second aspect of the present application, there is provided a conversion data processing equivalence adjustment device of an electromagnetic transient model, the device comprising: the electromagnetic model building unit is used for obtaining an electromagnetic transient model obtained by converting the electromechanical transient model; the electromagnetic model operation unit is used for operating the electromagnetic transient model, and if the line solution is wrong and the zero sequence parameter of the read wrong line is 0, the missing province parameter is complemented according to the following missing complement formula: if the line positive sequence is: g i1 -jB i1 ＝1/(R i1 +jX i1 ) The method comprises the steps of carrying out a first treatment on the surface of the The line zero sequence is: g i0 -jB i0 ＝0.8G i1 -j0.313B i1 ，G ic0 ＝0.7G ic1 ，B ic0 ＝0.68B ic1 The method comprises the steps of carrying out a first treatment on the surface of the In the above, G i1 Is the positive sequence conductance of line i, B i1 Is the positive sequence susceptance of line i, R i1 Is the positive sequence resistance of line i, X i1 Is the positive sequence reactance of line i; g i0 Is the zero sequence conductance of line i, B i0 Is the zero sequence susceptance of line i, G i c 0 Is the zero sequence conductance of the ground to the line i, G i c 1 Is the ground-pair positive sequence conductance of line i, B i c 0 Is the ground-to-ground zero sequence susceptance of line i, B ic1 The positive sequence susceptance of the ground pair of the line i; the bus voltage comparison unit is used for comparing the bus voltage results of the electromechanical transient model and the electromagnetic transient model when the electromagnetic transient model line solution is not wrong by the electromagnetic model operation unit; and the electromagnetic model output unit is used for perfecting the electromagnetic transient model according to the bus voltage deviation value in the comparison result to obtain a target electromagnetic transient model.
According to a third aspect of the present application there is provided an electronic device comprising a memory, a processor and a computer program stored on said memory and executable on said processor, the processor implementing the steps of the above method when executing said computer program.
According to a fourth aspect of the present application there is provided a computer readable storage medium having stored thereon a computer program which when executed by a processor performs the steps of the above method.
According to a fifth aspect of the present application there is provided a computer program product comprising computer programs/instructions which when executed by a processor implement the steps of the above method.
The conversion data processing equivalent adjustment method and device for the electromagnetic transient model can quickly perfect parameters of the electromagnetic transient model, does not need to manually check parameters of each element, and saves a great amount of time cost. The application has strong universality, is not specific to specific electromechanical software, and can be suitable for data adjustment after different electromechanical software are converted into the same electromagnetic software. The simulation support can be used for model data adjustment of scenes such as large power grid planning and design, and the like, and provides simulation support for constraint factors and mechanism analysis of power transmission and reception capacity of a planning grid.
Drawings
In order to more clearly illustrate the embodiments of the application or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the application, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art. In the drawings:
fig. 1 is a schematic flow chart of a conversion data processing equivalence adjustment method of an electromagnetic transient model according to an embodiment of the present application;
FIG. 2 is a schematic flow chart of a method for adjusting the equivalent of conversion data processing of an electromagnetic transient model according to another embodiment of the present application;
fig. 3 is a schematic structural diagram of a conversion data processing equivalent adjusting device of an electromagnetic transient model according to an embodiment of the present application;
fig. 4 is a schematic structural diagram of an electronic device according to an embodiment of the present application.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present application more apparent, the embodiments of the present application will be described in further detail with reference to the accompanying drawings. The exemplary embodiments of the present application and their descriptions herein are for the purpose of explaining the present application, but are not to be construed as limiting the application.
Fig. 1 is a schematic flow chart of a method for adjusting conversion data processing equivalence of an electromagnetic transient model according to an embodiment of the present application, where the method includes:
step S101: and obtaining an electromagnetic transient model obtained by converting the electromechanical transient model.
Step S102: and operating the electromagnetic transient model.
Step S103: if the line solving is in error and the zero sequence parameter of the read error line is 0, the default parameters are complemented according to the following complement formula:
if the line positive sequence is: g i1 -jB i1 ＝1/(R i1 +jX i1 )；
The line zero sequence is: g i0 -jB i0 ＝0.8G i1 -j0.313B i1 ，G ic0 ＝0.7G ic1 ，B ic0 ＝0.68B ic1 ；
In the above, G i1 Is the positive sequence conductance of line i, B i1 Is the positive sequence susceptance of line i, R i1 Is the positive sequence resistance of line i, X i1 Is the positive sequence reactance of line i; g i0 Is the zero sequence conductance of line i, B i0 Is the zero sequence susceptance of line i, G i c 0 Is the zero sequence conductance of the ground to the line i, G i c 1 Is the ground-pair positive sequence conductance of line i, B i c 0 Is the ground-to-ground zero sequence susceptance of line i, B ic1 The positive sequence susceptance of the ground pair of the line i;
step S104: if the line solving is not wrong, comparing the busbar voltage results of the electromechanical transient model and the electromagnetic transient model, and perfecting the electromagnetic transient model according to the busbar voltage deviation value in the comparison result to obtain a target electromagnetic transient model.
In this embodiment, the bus voltage deviation value in this step may include a bus average voltage deviation and a bus extreme voltage deviation, where the bus average voltage deviation satisfies the following index constraint:
ΔU iavg ＝|U iavg machine -U iavg electromagnetic |<ΔU avg threshold
In the formula DeltaU iavg For the average voltage deviation of the ith bus in the preset time, U iavg machine For the average voltage of the ith bus in the preset time in the electromechanical simulation, U iavg electromagnetic For the average voltage delta U of the ith bus in the preset time in electromagnetic simulation avg threshold Is the bus average voltage threshold;
the bus extreme voltage deviation meets the following index constraint:
ΔU iets ＝{|U imax electromechanical -U imax electromagnetic |、|U imin electromechanical -U imin electromagnetism |} max ＜ΔU est threshold
In the formula DeltaU iets U is the extreme voltage deviation of the ith bus in the preset time imax electromechanical 、U imin electromechanical Respectively the maximum voltage and the minimum voltage of the ith bus within the preset time in the electromechanical simulation, U imax electromagnetic 、U imin electromagnetism Respectively the maximum voltage and the minimum voltage of the ith bus within the preset time in electromagnetic simulation, delta U est threshold Is the bus extreme voltage threshold.
Further, this step S104: the electromagnetic transient model is perfected according to the bus voltage deviation value in the comparison result to obtain a target electromagnetic transient model, and the method comprises the following steps of:
step a: and (3) judging whether the average voltage deviation of the bus is less than 1%, if so, entering a step d, and if not, entering a step b.
Step b: searching a bus with the bus voltage of 0, checking parameters of elements connected with the bus, if any one of the element resistance, inductance and capacitance is 0, supplementing the deficiency-saving parameters according to the deficiency supplementing formula, and entering the step c; if the element parameters are not 0, the resistor R<-at 0.001, R is set to 0 and step c is entered when the inductance L<0, modifying the inductance element L into a capacitance element C and entering the step C,if the element parameters are all larger than 0, checking the impedance and entering the step c.
Step c: respectively detecting three-side node voltage U of three-winding transformer 1 、U 2 And U 3 If U is 1k -U 1h ＝U 2k -U 2h ＝U 3k -U 3h =0 (k+.h, k, h=1, …, t), t is the number of transformers in the area, change the three-winding transformer wiring to Y G Y G D 11 And comparing the bus voltage results again, otherwise, directly comparing the bus voltage results again.
Step d: judging whether the bus extreme voltage deviation is smaller than 1%, if so, obtaining a target electromagnetic transient model, and if not, increasing the generator q-axis transient reactance Xq' and comparing the bus voltage results.
The steps a-d are steps for judging the specific position where the problem is generated and modifying the parameters to solve the problem.
According to the technical scheme, the conversion data processing equivalent adjustment method of the electromagnetic transient model can be used for rapidly improving parameters of the electromagnetic transient model, manually checking parameters of each element is not needed, and a large amount of time and cost are saved. The application has strong universality, is not specific to specific electromechanical software, and can be suitable for data adjustment after different electromechanical software are converted into the same electromagnetic software. The simulation support can be used for model data adjustment of scenes such as large power grid planning and design, and the like, and provides simulation support for constraint factors and mechanism analysis of power transmission and reception capacity of a planning grid.
Fig. 2 is a schematic flow chart of a method for adjusting conversion data processing equivalence of an electromagnetic transient model according to another embodiment of the present application, where the method includes the following steps:
step S201: an electromagnetic transient model is generated.
Step S202: and running an electromagnetic transient model.
Step S203: judging whether the line solving is wrong, if so, returning to the step S202 after supplementing the deficiency saving parameters through the deficiency supplementing formula in the embodiment; if not, the process proceeds to step S204.
Step S204: and judging whether the average voltage deviations of the buses are smaller than 1%, if so, proceeding to step S208, and if not, proceeding to step S205.
Step S205: searching a bus with the bus voltage of 0, checking parameters of elements connected with the bus, judging whether any parameter of the element resistance, the inductance and the capacitance is 0, if so, entering step S207 after supplementing the deficiency and province parameters through a deficiency supplementing formula in the embodiment; if not, the process proceeds to step S206.
Step S206: continuing to determine whether the resistance R is less than-0.001 or whether the inductance L is less than 0, if so, modifying the element and element parameters and then proceeding to step S207, specifically, when the resistance R is less than-0.001, turning on the powerThe resistor R is set to 0, and when the inductance L is smaller than 0, the inductance element L is modified to the capacitance element C,if not, the impedance is checked and the process proceeds to step S207.
Step S207: judging whether the wiring modes of the parallel transformers are the same or not, specifically, detecting the voltages U of the three-side nodes of the three-winding transformer respectively 1 、U 2 And U 3 Judging whether different transformers are connected in parallel or not, if so, U 1k -U 1h ＝U 2k -U 2h ＝U 3k -U 3h =0 (k not equal to h, k, h=1, …, t), t is the number of transformers in the area, and indicates that the transformers are in parallel connection, and at this time, whether the parallel connection transformer wiring is Y is compared G Y G D 11 If not, the three-winding transformer wiring is changed into Y G Y G D 11 And returning to the step S202, otherwise, directly returning to the step S202.
Step S208: and judging whether the bus extreme voltage deviation is less than 1%, if so, obtaining a target electromagnetic model, and if not, increasing the generator q-axis transient reactance Xq' and returning to the step S202.
From the above, the conversion data processing equivalent adjustment method of the electromagnetic transient model provided by the application can quickly perfect the parameters of the electromagnetic transient model, does not need to manually check the parameters of each element, and saves a great amount of time and cost. The application has strong universality, is not specific to specific electromechanical software, and can be suitable for data adjustment after different electromechanical software are converted into the same electromagnetic software. The simulation support can be used for model data adjustment of scenes such as large power grid planning and design, and the like, and provides simulation support for constraint factors and mechanism analysis of power transmission and reception capacity of a planning grid.
Fig. 3 is a schematic structural diagram of a conversion data processing equivalence adjusting device for an electromagnetic transient model according to an embodiment of the present application, where the device includes: the electromagnetic model building unit 310, the electromagnetic model running unit 320, the bus voltage comparison unit 330 and the electromagnetic model output unit 340 are sequentially connected with each other. Wherein:
the electromagnetic model building unit 310 is configured to obtain an electromagnetic transient model obtained by converting the electromechanical transient model;
the electromagnetic model operation unit 320 is configured to operate the electromagnetic transient model, and if the line solution is in error and the zero sequence parameter of the read error line is 0, the default parameters are complemented according to the following complement formula:
if the line positive sequence is: g i1 -jB i1 ＝1/(R i1 +jX i1 )；
The line zero sequence is: g i0 -jB i0 ＝0.8G i1 -j0.313B i1 ，G ic0 ＝0.7G ic1 ，B ic0 ＝0.68B ic1 ；
In the above, G i1 Is the positive sequence conductance of line i, B i1 Is the positive sequence susceptance of line i, R i1 Is the positive sequence resistance of line i, X i1 Is the positive sequence reactance of line i; g i0 Is the zero sequence conductance of line i, B i0 Is the zero sequence susceptance of line i, G i c 0 Is the zero sequence conductance of the ground to the line i, G i c 1 Is the ground-pair positive sequence conductance of line i, B i c 0 Is the ground-to-ground zero sequence susceptance of line i, B ic1 The positive sequence susceptance of the ground pair of the line i;
the bus voltage comparison unit 330 is configured to compare bus voltage results of the electromechanical transient model and the electromagnetic transient model when the electromagnetic transient model line solution is not in error by the electromagnetic model operation unit;
the electromagnetic model output unit 340 is configured to perfect the electromagnetic transient model according to the bus voltage deviation value in the comparison result to obtain a target electromagnetic transient model.
Preferably, the bus voltage deviation value includes a bus average voltage deviation and a bus extreme voltage deviation, where the bus average voltage deviation satisfies the following index constraint:
ΔU iavg ＝|U iavg machine -U iavg electromagnetic |<ΔU avg threshold
In the formula DeltaU iavg For the average voltage deviation of the ith bus in the preset time, U iavg machine For the average voltage of the ith bus in the preset time in the electromechanical simulation, U iavg electromagnetic For the average voltage delta U of the ith bus in the preset time in electromagnetic simulation avg threshold Is the bus average voltage threshold;
the bus extreme voltage deviation meets the following index constraint:
ΔU iets ＝{|U imax electromechanical -U imax electromagnetic |、|U imin electromechanical -U imin electromagnetism |} max ＜ΔU est threshold
In the formula DeltaU iets U is the extreme voltage deviation of the ith bus in the preset time imax electromechanical 、U imin electromechanical Respectively the maximum voltage and the minimum voltage of the ith bus within the preset time in the electromechanical simulation, U imax electromagnetic 、U imin electromagnetism Respectively the maximum voltage and the minimum voltage of the ith bus within the preset time in electromagnetic simulation, delta U est threshold Is the bus extreme voltage threshold.
Preferably, the electromagnetic model output unit 340 is specifically configured to implement the following steps: :
a: judging whether the average voltage deviation of the bus is less than 1%, if so, entering a step d, and if not, entering a step b;
b: searching a bus with the bus voltage of 0, checking parameters of elements connected with the bus, if any one of the element resistance, inductance and capacitance is 0, supplementing the deficiency-saving parameters according to the deficiency supplementing formula, and entering the step c; if the element parameters are not 0, the resistor R<-at 0.001, R is set to 0 and step c is entered when the inductance L<0, modifying the inductance element L into a capacitance element C and entering the step C,if the element parameters are all larger than 0, checking the impedance and entering the step c;
c: respectively detecting three-side node voltage U of three-winding transformer 1 、U 2 And U 3 If U is 1k -U 1h ＝U 2k -U 2h ＝U 3k -U 3h =0 (k+.h, k, h=1, …, t), t is the number of transformers in the area, change the three-winding transformer wiring to Y G Y G D 11 And comparing the bus voltage results again;
d: and judging whether the bus extreme voltage deviation is less than 1%, if so, obtaining a target electromagnetic transient model, and if not, increasing the generator q-axis transient reactance Xq' and comparing the bus voltage results again.
From the above, the conversion data processing equivalent adjusting device of the electromagnetic transient model provided by the application can quickly perfect the parameters of the electromagnetic transient model, does not need to manually check the parameters of each element, and saves a great amount of time and cost. The application has strong universality, is not specific to specific electromechanical software, and can be suitable for data adjustment after different electromechanical software are converted into the same electromagnetic software. The simulation support can be used for model data adjustment of scenes such as large power grid planning and design, and the like, and provides simulation support for constraint factors and mechanism analysis of power transmission and reception capacity of a planning grid.
Fig. 4 is a schematic diagram of an electronic device according to an embodiment of the present application. The electronic device shown in fig. 4 is a general-purpose data processing apparatus comprising a general-purpose computer hardware structure including at least a processor 801 and a memory 802. The processor 801 and the memory 802 are connected by a bus 803. The memory 802 is adapted to store one or more instructions or programs executable by the processor 801. The one or more instructions or programs are executed by the processor 801 to implement the steps in the transformation data processing equivalent adjustment method for electromagnetic transient models described above.
The processor 801 may be a separate microprocessor or a collection of one or more microprocessors. Thus, the processor 801 performs the process of processing data and controlling other devices by executing the commands stored in the memory 802, thereby executing the method flow of the embodiment of the present application as described above. The bus 803 connects the above-described components together, while connecting the above-described components to a display controller 804 and a display device and an input/output (I/O) device 805. Input/output (I/O) devices 805 may be a mouse, keyboard, modem, network interface, touch input device, somatosensory input device, printer, and other devices known in the art. Typically, input/output (I/O) devices 805 are connected to the system through input/output (I/O) controllers 806.
The memory 802 may store software components such as an operating system, communication modules, interaction modules, and application programs, among others. Each of the modules and applications described above corresponds to a set of executable program instructions that perform one or more functions and methods described in the embodiments of the application.
The embodiment of the application also provides a computer readable storage medium, on which a computer program is stored, which when being executed by a processor, is used for realizing the steps of the conversion data processing equivalent adjustment method of the electromagnetic transient model.
The embodiment of the application also provides a computer program product, which comprises a computer program/instruction, wherein the computer program/instruction realizes the steps of the conversion data processing equivalent adjustment method of the upper electromagnetic transient model when being executed by a processor.
From the above, the method and the device for adjusting the equivalent of the conversion data processing of the electromagnetic transient model can quickly perfect the parameters of the electromagnetic transient model, do not need to manually check the parameters of each element, and save a great deal of time cost. The method is strong in universality, does not aim at specific electromechanical software, and can be suitable for data adjustment after different electromechanical software is converted into the same electromagnetic software. The simulation support can be used for model data adjustment of scenes such as large power grid planning and design, and the like, and provides simulation support for constraint factors and mechanism analysis of power transmission and reception capacity of a planning grid.
Preferred embodiments of the present application are described above with reference to the accompanying drawings. The many features and advantages of the embodiments are apparent from the detailed specification, and thus, it is intended by the appended claims to cover all such features and advantages of the embodiments which fall within the true spirit and scope thereof. Further, since numerous modifications and changes will readily occur to those skilled in the art, it is not desired to limit the embodiments of the application to the exact construction and operation illustrated and described, and accordingly, all suitable modifications and equivalents may be resorted to, falling within the scope thereof.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The foregoing description of the embodiments has been provided for the purpose of illustrating the general principles of the application, and is not meant to limit the scope of the application, but to limit the application to the particular embodiments, and any modifications, equivalents, improvements, etc. that fall within the spirit and principles of the application are intended to be included within the scope of the application.
Claims (9)
1. A method for converting data processing equivalence adjustment of an electromagnetic transient model, the method comprising:
acquiring an electromagnetic transient model obtained by converting an electromechanical transient model;
and (3) operating the electromagnetic transient model, and if the line solving is in error and the zero sequence parameter of the read error line is 0, supplementing the default parameter according to the following supplementing formula:
if the line positive sequence is: g i1 -jB i1 ＝1/(R i1 +jX i1 )；
The line zero sequence is: g i0 -jB i0 ＝0.8G i1 -j0.313B i1 ，G ic0 ＝0.7G ic1 ，B ic0 ＝0.68B ic1 ；
In the above, G i1 Is the positive sequence conductance of line i, B i1 Is the positive sequence susceptance of line i, R i1 Is the positive sequence resistance of line i, X i1 Is the positive sequence reactance of line i; g i0 Is the zero sequence conductance of line i, B i0 Is the zero sequence susceptance of line i, G i c 0 Is the zero sequence conductance of the ground to the line i, G i c 1 Is the ground-pair positive sequence conductance of line i, B i c 0 Is the ground-to-ground zero sequence susceptance of line i, B ic1 The positive sequence susceptance of the ground pair of the line i;
if the line solving is not wrong, comparing the busbar voltage results of the electromechanical transient model and the electromagnetic transient model, and perfecting the electromagnetic transient model according to the busbar voltage deviation value in the comparison result to obtain a target electromagnetic transient model.
2. The method for converting data processing equivalence adjustment of an electromagnetic transient model according to claim 1, wherein the bus voltage deviation value comprises a bus average voltage deviation and a bus extreme voltage deviation, wherein the bus average voltage deviation satisfies the following index constraint:
ΔU iavg ＝|U iavg machine -U iavg electromagnetic |<ΔU avgthreshold
In the formula DeltaU iavg For the average voltage deviation of the ith bus in the preset time, U iavg machine For the average voltage of the ith bus in the preset time in the electromechanical simulation, U iavg electromagnetic For the average voltage delta U of the ith bus in the preset time in electromagnetic simulation avg threshold Is the bus average voltage threshold;
the bus extreme voltage deviation meets the following index constraint:
ΔU iets ＝{|U imax electromechanical -U imax electromagnetic |、|U imin electromechanical -U imin electromagnetism |} max ＜ΔU est threshold
In the formula DeltaU iets U is the extreme voltage deviation of the ith bus in the preset time imax electromechanical 、U imin electromechanical Respectively the maximum voltage and the minimum voltage of the ith bus within the preset time in the electromechanical simulation, U imax electromagnetic 、U imin electromagnetism Respectively the maximum voltage and the minimum voltage of the ith bus within the preset time in electromagnetic simulation, delta U estthreshold Is the bus extreme voltage threshold.
3. The method for adjusting the equivalent value of the conversion data processing of the electromagnetic transient model according to claim 2, wherein the step of perfecting the electromagnetic transient model according to the bus voltage deviation value in the comparison result to obtain the target electromagnetic transient model comprises the following steps:
a: judging whether the average voltage deviation of the bus is less than 1%, if so, entering a step d, and if not, entering a step b;
b：searching a bus with the bus voltage of 0, checking parameters of elements connected with the bus, if any one of the element resistance, inductance and capacitance is 0, supplementing the deficiency-saving parameters according to the deficiency supplementing formula, and entering the step c; if the element parameters are not 0, the resistor R<-at 0.001, R is set to 0 and step c is entered when the inductance L<0, modifying the inductance element L into a capacitance element C and entering the step C,if the element parameters are all larger than 0, checking the impedance and entering the step c;
c: respectively detecting three-side node voltage U of three-winding transformer 1 、U 2 And U 3 If U is 1k -U 1h ＝U 2k -U 2h ＝U 3k -U 3h =0 (k+.h, k, h=1, …, t), t is the number of transformers in the area, change the three-winding transformer wiring to Y G Y G D 11 And comparing the bus voltage results again;
d: and judging whether the bus extreme voltage deviation is less than 1%, if so, obtaining a target electromagnetic transient model, and if not, increasing the generator q-axis transient reactance Xq' and comparing the bus voltage results again.
4. A conversion data processing equivalence adjusting device of an electromagnetic transient model, characterized in that the device comprises:
the electromagnetic model building unit is used for obtaining an electromagnetic transient model obtained by converting the electromechanical transient model;
the electromagnetic model operation unit is used for operating the electromagnetic transient model, and if the line solution is wrong and the zero sequence parameter of the read wrong line is 0, the missing province parameter is complemented according to the following missing complement formula:
if the line positive sequence is: g i1 -jB i1 ＝1/(R i1 +jX i1 )；
The line zero sequence is: g i0 -jB i0 ＝0.8G i1 -j0.313B i1 ，G ic0 ＝0.7G ic1 ，B ic0 ＝0.68B ic1 ；
In the above, G i1 Is the positive sequence conductance of line i, B i1 Is the positive sequence susceptance of line i, R i1 Is the positive sequence resistance of line i, X i1 Is the positive sequence reactance of line i; g i0 Is the zero sequence conductance of line i, B i0 Is the zero sequence susceptance of line i, G i c 0 Is the zero sequence conductance of the ground to the line i, G i c 1 Is the ground-pair positive sequence conductance of line i, B i c 0 Is the ground-to-ground zero sequence susceptance of line i, B ic1 The positive sequence susceptance of the ground pair of the line i;
the bus voltage comparison unit is used for comparing the bus voltage results of the electromechanical transient model and the electromagnetic transient model when the electromagnetic transient model line solution is not wrong by the electromagnetic model operation unit;
and the electromagnetic model output unit is used for perfecting the electromagnetic transient model according to the bus voltage deviation value in the comparison result to obtain a target electromagnetic transient model.
5. The conversion data processing equivalence adjustment device of an electromagnetic transient model according to claim 4, wherein the bus voltage deviation value comprises a bus average voltage deviation and a bus extreme voltage deviation, wherein the bus average voltage deviation satisfies the following index constraint:
ΔU iavg ＝|U iavg machine -U iavg electromagnetic |<ΔU avg threshold
In the formula DeltaU iavg For the average voltage deviation of the ith bus in the preset time, U iavg machine For the average voltage of the ith bus in the preset time in the electromechanical simulation, U iavg electromagnetic For the average voltage delta U of the ith bus in the preset time in electromagnetic simulation avg threshold Is the bus average voltage threshold;
the bus extreme voltage deviation meets the following index constraint:
ΔU iets ＝{|U imax electromechanical -U imax electromagnetic |、|U imin electromechanical -U imin electromagnetism |} max ＜ΔU est threshold
In the formula DeltaU iets U is the extreme voltage deviation of the ith bus in the preset time imax electromechanical 、U imin electromechanical Respectively the maximum voltage and the minimum voltage of the ith bus within the preset time in the electromechanical simulation, U imax electromagnetic 、U imin electromagnetism Respectively the maximum voltage and the minimum voltage of the ith bus within the preset time in electromagnetic simulation, delta U est threshold Is the bus extreme voltage threshold.
6. The conversion data processing equivalent adjustment device of an electromagnetic transient model according to claim 5, wherein the electromagnetic model output unit is specifically configured to implement the following steps: :
a: judging whether the average voltage deviation of the bus is less than 1%, if so, entering a step d, and if not, entering a step b;
b: searching a bus with the bus voltage of 0, checking parameters of elements connected with the bus, if any one of the element resistance, inductance and capacitance is 0, supplementing the deficiency-saving parameters according to the deficiency supplementing formula, and entering the step c; if the element parameters are not 0, the resistor R<-at 0.001, R is set to 0 and step c is entered when the inductance L<0, modifying the inductance element L into a capacitance element C and entering the step C,if the element parameters are all larger than 0, checking the impedance and entering the step c;
c: respectively detecting three-side node voltage U of three-winding transformer 1 、U 2 And U 3 If U is 1k -U 1h ＝U 2k -U 2h ＝U 3k -U 3h =0 (k+.h, k, h=1, …, t), t is the number of transformers in the area, change the three-winding transformer wiring to Y G Y G D 11 And comparing the bus voltage results again;
d: and judging whether the bus extreme voltage deviation is less than 1%, if so, obtaining a target electromagnetic transient model, and if not, increasing the generator q-axis transient reactance Xq' and comparing the bus voltage results again.
7. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the steps of the method of any of claims 1 to 3 when the computer program is executed by the processor.
8. A computer readable storage medium, on which a computer program is stored, characterized in that the computer program, when being executed by a processor, implements the steps of the method according to any one of claims 1 to 3.
9. A computer program product comprising computer programs/instructions which, when executed by a processor, implement the steps of the method of any one of claims 1 to 3.
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2022108919343 | 2022-07-27 | ||
| CN202210891934.3A CN115293073A (en) | 2022-07-27 | 2022-07-27 | Method and device for converting BPA electromechanical transient model into ADPSS electromagnetic transient model |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117057284A true CN117057284A (en) | 2023-11-14 |
Family
ID=83824033
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210891934.3A Pending CN115293073A (en) | 2022-07-27 | 2022-07-27 | Method and device for converting BPA electromechanical transient model into ADPSS electromagnetic transient model |
| CN202310627934.7A Pending CN117057284A (en) | 2022-07-27 | 2023-05-30 | Conversion data processing equivalent adjustment method and device for electromagnetic transient model |
Family Applications Before (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210891934.3A Pending CN115293073A (en) | 2022-07-27 | 2022-07-27 | Method and device for converting BPA electromechanical transient model into ADPSS electromagnetic transient model |
Country Status (1)
| Country | Link |
|---|---|
| CN (2) | CN115293073A (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117763808B (en) * | 2023-11-30 | 2024-07-02 | 中国电力科学研究院有限公司 | Electromagnetic transient automatic modeling method of direct-current transmission system based on template rule |
-
2022
- 2022-07-27 CN CN202210891934.3A patent/CN115293073A/en active Pending
-
2023
- 2023-05-30 CN CN202310627934.7A patent/CN117057284A/en active Pending
Also Published As
| Publication number | Publication date |
|---|---|
| CN115293073A (en) | 2022-11-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110212521B (en) | Power distribution network line loss calculation method and system | |
| Song et al. | Almost sure stability of switching Markov jump linear systems | |
| CN117057284A (en) | Conversion data processing equivalent adjustment method and device for electromagnetic transient model | |
| CN102946253B (en) | Data sampling method and system and the application process in parameter identification thereof and system | |
| CN109638812B (en) | Self-adaptive distribution line ultra-short term load prediction method and system | |
| Neçaibia et al. | Self-tuning fractional order PI λ D µ controller based on extremum seeking approach | |
| US20080040292A1 (en) | Cost information management system, cost information management method, and cost information management program | |
| CN110429618B (en) | Micro-grid frequency modulation control method, control system and device | |
| CN115562069A (en) | RT-LAB-based energy storage grid-connected test system and test method thereof | |
| CN112383044A (en) | Power grid model comparison method and device based on hierarchical topological structure | |
| CN101154818B (en) | Method and system for parallel connection of UPS derated models with different capacitance grade | |
| US11699013B2 (en) | Computer-implemented method for simulation of an electrical circuit | |
| Su et al. | Distributed load sharing and transmission power loss optimisation for DC microgrids | |
| CN106856334B (en) | A kind of power system state estimation method considering flexible direct current control characteristic | |
| Boulet et al. | The fundamental tradeoff between performance and robustness-a new perspective on loop shaping-classic control revisited part ii | |
| CN105406749B (en) | The robust Controller Design method of gird-connected inverter parallel system | |
| CN116382118A (en) | Simulation method and device of photovoltaic power station and electronic equipment | |
| CN104182566B (en) | For the solution element anode voltage method after part changes to occur in network system | |
| CN114077599B (en) | Method, device and system for constructing transformer sensing quantity database supporting dynamic data repair | |
| Almakhles et al. | Single-bit modulator based controller for capacitive power transfer system | |
| CN109657332B (en) | Method and system for decoupling electromagnetic transient automatic modeling of large-scale power grid | |
| CN106598911B (en) | A kind of column write method of circuit network state equation | |
| CN113111623A (en) | Design method and system of direct current loop impedance | |
| CN111725818A (en) | Three-phase weak current network converter group grid-connected simulation method and simulation terminal | |
| CN111082414A (en) | Transient voltage calculation method and system |
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 |