External characteristic aggregation equivalent interaction method and device for new energy convergence area, and medium

Abstract

The invention discloses an external characteristic aggregation equivalent interaction method and device for a new energy gathering area, and a medium, and belongs to the technical field of new energy external characteristic equivalent interaction. According to the method, boundary plant stations and adjacent plant stations are searched, and marking processing is performed according to whether new energy components are included, so that a mathematical model of a power grid can be constructed more accurately, and existence and influence of the new energy plant stations can be reflected accurately when provincial and regional interface equivalence is calculated; through arrangement of a transition resistor and iterative computation, the change of new energy fault current can be simulated, so that the influence of new energy short-circuit current on an equivalent model and a relay protection constant value is fully considered, and the equivalent model is closer to the actual power grid operation condition; by fully considering the influence of the new energy short-circuit current, an equivalent interaction mode is improved, the accuracy of an equivalent model and the adaptability of a relay protection constant value are improved, and a powerful guarantee is provided for safe and stable operation of a power grid system.

Description

Technical Field

[0001] The present invention relates to the technical field of equivalent interaction of new energy external characteristics, and in particular to an equivalent interaction method for external characteristic aggregation in a new energy gathering area. Background Art

[0002] Multi-level control agencies rely on unified data models, graphical models, setting calculation basic platforms and setting calculation principles to achieve mutual connection and efficient coordination of power grid relay protection setting calculation. Under the current technical architecture, both the provincial control center (provincial control center) and the regional control center (local control center) have adopted an integrated setting calculation system, which adopts a distributed deployment strategy. Specifically, the provincial control center and the local control center are respectively equipped with independent servers, and each is responsible for maintaining the relevant basic data for setting calculation. In terms of the interaction of power grid models, the provincial and local control centers exchange data at the interface in an equivalent form, in which the provincial control center sends the equivalent model of the deployed area to the local control center, and the local control center reports the corresponding equivalent model to the provincial control center. This mechanism ensures that the provincial control center and the local control center can operate based on a unified and complete power grid model when performing setting calculations, thereby effectively ensuring the accuracy of short-circuit current and branch coefficient calculations, as well as the correctness of relay protection settings.

[0003] However, as the penetration rate of renewable energy in local power grids continues to rise, the behavior of the grid in the event of a fault has changed significantly. Renewable energy stations, with their low voltage ride-through characteristics, will inject short-circuit current into the grid when a grid fault occurs. This injection behavior has a significant impact on the sensitivity and selectivity of relay protection settings, which may cause protection failure or misoperation. At present, the equivalent model of the provincial-local interface has not fully considered this impact of renewable energy, which has weakened the accuracy and reliability of provincial relay protection settings to a certain extent.

[0004] In view of the increasingly significant impact of new energy short-circuit current on the equivalent model of the provincial-earth interface and the setting value of relay protection, the constant impedance equivalent interaction mode in the existing technology can no longer meet the setting calculation requirements under the new power system. Therefore, how to fully consider the impact of new energy short-circuit current when calculating the equivalent value of the provincial-earth interface, so as to improve the existing equivalent interaction mode, has become a technical problem that the current technicians in this field need to overcome urgently. Summary of the invention

[0005] The purpose of the present invention is to provide an equivalent interaction method, device and medium for external characteristic aggregation in a new energy convergence area, aiming to overcome the problem in the prior art that the accuracy and reliability of relay protection settings are weakened due to the failure to fully incorporate new energy factors into the equivalent interaction model.

[0006] The present invention solves the above technical problems through the following technical solutions: A method for external characteristic aggregation equivalent interaction of a new energy gathering area includes the following steps: S1. Find all boundary plants and stations based on the voltage level of the local power grid; S2, starting from the boundary plant and station, search for adjacent plants and stations along the voltage level of the local adjustment grid; S3. Based on whether the plant station contains new energy components, the busbar of the plant station is marked. If the plant station contains new energy components, the busbar of the plant station is marked as a new energy busbar; if the plant station does not contain new energy components, the busbar of the plant station is marked as a conventional power busbar, and the plant station includes a boundary plant station and an adjacent plant station; based on the result of the marking process, a mathematical model is generated; S4. When the bus to be calculated is a conventional power bus, according to the generated mathematical model, the equivalent model of each conventional power bus fault node is calculated in turn, and reported to the provincial dispatching; when the bus to be calculated is a new energy bus, according to the generated mathematical model, by setting the transition resistance, the new energy fault current injected into the new energy bus is changed, and the external characteristic equivalent model of each new energy bus fault node is iteratively calculated, and reported to the provincial dispatching; S5. The provincial dispatching agency aggregates the equivalent model of the conventional power bus and the external characteristic equivalent model of the new energy bus to obtain an external characteristic aggregated equivalent model. The provincial dispatching agency interactively controls the new energy gathering area based on the external characteristic aggregated equivalent model.

[0007] A further improvement of the present invention is that the new energy component is a photovoltaic element or a wind turbine element.

[0008] A further improvement of the present invention is that: S3 specifically includes the following steps: traversing all power plants and stations, marking the bus of the power plants and stations, if the power plants and stations contain new energy components, marking the bus of the power plants and stations as a new energy bus; if the power plants and stations do not contain new energy components, marking the bus of the power plants and stations as a conventional power bus, and recording the conventional power bus and the new energy bus in a conventional power bus set A and a new energy bus set B, respectively, wherein the conventional power bus set A contains N conventional power buses, and the new energy bus set B contains M new energy buses; based on the result of the marking process, a mathematical model is generated in combination with the topological connection relationship between the boundary power plants and stations, adjacent power plants and stations and the ground-dispatching power grid, as well as the impedance information and operating status of the ground-dispatching equipment.

[0009] A further improvement of the present invention is that when the bus to be calculated is a conventional power bus, the present invention further includes the following steps: setting r=0, where r represents the position of the current conventional power bus in the conventional power bus set A; when the bus to be calculated is a new energy bus, the present invention further includes the following steps: setting s=0, where s represents the position of the current new energy bus in the new energy bus set B.

[0010] The further improvement of the present invention is that: the equivalent model of each conventional power bus fault node is calculated in sequence, specifically including the following steps: let r=r+1, make the current conventional power bus operate in the maximum operation mode, and calculate the first equivalent information of the current conventional power bus fault node f based on the mathematical model of S3, including the positive sequence impedance Z f1max and zero sequence impedance Z f0max ; Make the current conventional power bus in the minimum operating mode, based on the mathematical model of S3, calculate the second equivalent information of the current conventional power bus fault node f, including the positive sequence impedance Z f1min and zero sequence impedance Z f0min ; Determine whether r is equal to N. If it is not, execute S5 again; if it is yes, report all the calculated first equivalent information and second equivalent information as the equivalent model of the conventional power bus to the provincial dispatching department; Among them, the maximum operating mode refers to an operating mode in which the current conventional power bus has the smallest short-circuit impedance value when operating in this operating mode, and the short-circuit current generated after a short circuit occurs is the largest; the minimum operating mode refers to an operating mode in which the current conventional power bus has the largest short-circuit impedance when operating in this operating mode, and the short-circuit current generated after a short circuit occurs is the smallest.

[0011] A further improvement of the present invention is that the iterative calculation of the external characteristic equivalent model of each new energy bus fault node specifically includes the following steps: S41, let s=s+1, i=1, initialize the transition resistance , set the upper limit of transition resistance to R lim , set the transition resistance increment to , where i represents the number of iterations of the transition resistance; k is set to 0, k represents the number of times the new energy fault current is injected into the current new energy bus; S42, let k = k + 1, obtain the terminal voltage, inject the kth new energy fault current into the current new energy bus, and calculate the first current parameters of all new energy nodes p associated with the plant corresponding to the current new energy bus in the mathematical model based on the terminal voltage, including the positive sequence current and negative sequence current ; S43, obtain the positive sequence impedance, negative sequence impedance and zero sequence impedance of the conventional power supply equivalent of the current new energy bus, where the positive sequence impedance and negative sequence impedance are both Z f1 , zero sequence impedance is Z f0 , according to the first current parameter, calculate the second current parameter of the current new energy bus fault node f, including the positive sequence current and negative sequence current ; S44, based on the first current parameter and the second current parameter, calculate the first voltage parameters of all new energy nodes p associated with the plant corresponding to the current new energy bus in the mathematical model, including the positive sequence voltage , negative sequence voltage and zero sequence voltage ; Determine whether the difference between the current first voltage parameter and the last first voltage parameter meets the convergence accuracy. When k=1, the terminal voltage is the last first voltage parameter. If the difference meets the convergence accuracy, execute S45; if the difference does not meet the convergence accuracy, return to S42; S45, let i=i+1, current transition resistance , is the last transition resistance, judge Is it greater than or equal to the upper limit of transition resistance R lim If the judgment is no, then return to S42; if the judgment is yes, then enter S46; S46, obtaining the second current parameter of the current fault node f, and calculating the second voltage parameter of the current fault node f according to the second current parameter, including the positive sequence voltage , negative sequence voltage , based on the second current parameter and the second voltage parameter, generate a fault amplitude curve and a phase angle curve; determine whether s is equal to M, if not, return to S41; if yes, report all fault amplitude curves and phase angle curves to the provincial dispatching department as the external characteristic equivalent model of the new energy bus.

[0012] A further improvement of the present invention is that the first current parameter is specifically:

[0013] Where j is the symbol of the imaginary part of the complex number; , , are the steady-state positive-sequence short-circuit current, active component and reactive component corresponding to the new energy node p in the kth cycle under the i-th transition resistance; , , They are respectively the steady-state negative-sequence short-circuit current amplitude, active component and reactive component of the plant corresponding to the new energy node p after the current new energy bus fails; the superscript i represents the number of iterations of the i-th transition resistance; the superscript k represents the number of times the new energy fault current is injected into the current new energy bus; the subscript d represents the active component of the current; the subscript q represents the reactive component of the current; the subscript 1 represents the positive-sequence current; the subscript 2 represents the negative-sequence current; To reflect the new energy current and positive sequence voltage controlled by its terminal voltage Function of To reflect the new energy current and negative sequence voltage controlled by its terminal voltage Function of The second current parameter is specifically: If the fault is a single-phase grounding fault, then:

[0014] in, is the positive sequence short-circuit current of the fault node in the kth cycle of the new energy under the i-th transition resistance, is the negative sequence short-circuit current of the fault node in the kth cycle of the new energy under the i-th transition resistance, is the positive sequence impedance of the fault node f, is the zero-sequence impedance of the fault node f, is the grounding resistance of the single-phase short-circuit fault in the i-th cycle, is the positive sequence mutual impedance between the fault node f and the new energy node p; If the fault is a two-phase fault, then:

[0015] If the fault is a two-phase grounding fault, then:

[0016]

[0017] Among them, / / represents AND operation; If the fault is a three-phase phase-to-phase fault, the positive sequence current and negative sequence current of the fault node are:

[0018] .

[0019] A further improvement of the present invention is that the first voltage parameter is specifically:

[0020]

[0021]

[0022] in, is the positive sequence voltage of the new energy node p in the kth cycle of the new energy under the i-th transition resistance, is the negative sequence voltage of the new energy node p in the kth cycle of the new energy under the i-th transition resistance, is the zero-sequence voltage of the new energy node p in the kth cycle of the new energy under the i-th transition resistance, , are the positive-sequence and negative-sequence self-impedances of the new energy node p, respectively. is the negative sequence short-circuit current provided by the new energy node p for the kth time in the i-th cycle, is the zero-sequence short-circuit current provided by the fault node f in the kth cycle of the new energy under the i-th transition resistance, is the negative sequence mutual impedance between the fault node f and the new energy node p, is the zero-sequence mutual impedance between the fault node f and the new energy node p; The second voltage parameter is specifically:

[0023]

[0024] in, , Faulty nodes The positive and negative sequence self-impedance.

[0025] A terminal device comprises a memory, a processor and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the steps of the above method are implemented.

[0026] A computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the above method are implemented.

[0027] Compared with the prior art, the positive and progressive effects of the present invention are: The external characteristic aggregation equivalent interaction method of the new energy gathering area provided by the present invention can more accurately construct the mathematical model of the power grid by searching for boundary plants and adjacent plants and stations and marking them according to whether they contain new energy components, which helps to accurately reflect the existence and influence of new energy sites when calculating the equivalent values ​​of provincial and local interfaces; by setting transition resistance and iterative calculation, it can simulate the change of new energy fault current, so as to fully consider the influence of new energy short-circuit current on the equivalent model and relay protection setting, so that the equivalent model is closer to the actual power grid operation; by fully considering the influence of new energy short-circuit current, the equivalent interaction mode is improved, the accuracy of the equivalent model and the adaptability of the relay protection setting are improved, and a strong guarantee is provided for the safe and stable operation of the power grid system. BRIEF DESCRIPTION OF THE DRAWINGS

[0028] The drawings in the specification are used to provide further understanding of the present invention and constitute a part of the present invention. The schematic embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute improper limitations on the present invention.

[0029] Figure 1 It is a schematic diagram of the process of the present invention; Figure 2 It is a flow chart of conventional power bus; Figure 3 This is a flow chart of the new energy bus. DETAILED DESCRIPTION

[0030] In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Generally, the components of the embodiments of the present invention described and shown in the drawings here can be arranged and designed in various different configurations.

[0031] Therefore, the following detailed description of the embodiments of the present invention provided in the accompanying drawings is not intended to limit the scope of the invention claimed for protection, but merely represents selected embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0032] It should be noted that similar reference numerals and letters denote similar items in the following drawings, and therefore, once an item is defined in one drawing, further definition and explanation thereof is not required in subsequent drawings.

[0033] In the description of the embodiments of the present invention, it should be noted that if the terms "upper", "lower", "horizontal", "inner", etc. indicate an orientation or positional relationship based on the orientation or positional relationship shown in the drawings, or the orientation or positional relationship in which the product of the invention is usually placed when in use, it is only for the convenience of describing the present invention and simplifying the description, and does not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as a limitation on the present invention. In addition, the terms "first", "second", etc. are only used to distinguish the description, and cannot be understood as indicating or implying relative importance.

[0034] In addition, it should be noted that, unless otherwise clearly specified and limited, the terms "set", "install", "connect", and "connect" should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection, or it can be indirectly connected through an intermediate medium, or it can be the internal connection of two components. For ordinary technicians in this field, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0035] The present invention is further described in detail below in conjunction with the accompanying drawings and specific embodiments, which are intended to explain the present invention rather than to limit it.

[0036] See also Figure 1The present invention provides a method for exogenous characteristic aggregation equivalent interaction of a new energy convergence area, comprising the following steps: S1. Find all boundary plants and stations based on the voltage level of the local power grid; S2, starting from the boundary plant and station, search for adjacent plants and stations along the voltage level of the local adjustment grid; S3. Based on whether the plant station contains new energy components, the busbar of the plant station is marked. If the plant station contains new energy components, the busbar of the plant station is marked as a new energy busbar; if the plant station does not contain new energy components, the busbar of the plant station is marked as a conventional power busbar, and the plant station includes a boundary plant station and an adjacent plant station; based on the result of the marking process, a mathematical model is generated; S4. When the bus to be calculated is a conventional power bus, according to the generated mathematical model, the equivalent model of each conventional power bus fault node is calculated in turn, and reported to the provincial dispatching; when the bus to be calculated is a new energy bus, according to the generated mathematical model, by setting the transition resistance, the new energy fault current injected into the new energy bus is changed, and the external characteristic equivalent model of each new energy bus fault node is iteratively calculated, and reported to the provincial dispatching; S5. The provincial dispatching agency aggregates the equivalent model of the conventional power bus and the external characteristic equivalent model of the new energy bus to obtain an external characteristic aggregated equivalent model. The provincial dispatching agency interactively controls the new energy gathering area based on the external characteristic aggregated equivalent model.

[0037] By finding boundary plants and adjacent plants and marking them according to whether they contain new energy components, the mathematical model of the power grid can be constructed more accurately, which helps to accurately reflect the existence and influence of new energy sites when calculating the equivalent values ​​of provincial and local interfaces; by setting transition resistance and iterative calculation, the change of new energy fault current can be simulated, so as to fully consider the influence of new energy short-circuit current on the equivalent model and relay protection settings, making the equivalent model closer to the actual power grid operation; by fully considering the influence of new energy short-circuit current, the equivalent interaction mode is improved, the accuracy of the equivalent model and the adaptability of the relay protection settings are improved, and a strong guarantee is provided for the safe and stable operation of the power grid system.

[0038] In a specific embodiment of the present invention, S2, taking the boundary plant as the starting point, searching for adjacent plants along the voltage level of the ground-controlled power grid, specifically includes the following steps: taking each boundary plant as the starting point, along the voltage level of the ground-controlled power grid, based on the breadth of the electrical distance characteristics, searching for all adjacent plants L1 within the first-level range of the boundary plant; taking the adjacent plant L1 as the starting point, expanding outward along the voltage level of the ground-controlled power grid, searching for the adjacent plant L2 corresponding to the adjacent plant L1, the adjacent plant L2 being the plant within the first-level range to the second-level range of the boundary plant; taking the adjacent plant L2 as the starting point, expanding outward along the voltage level of the ground-controlled power grid, searching for the adjacent plant L3 corresponding to the adjacent plant L2; the level range for searching the boundary plant can be considered as a setting, generally the default is level 3.

[0039] Since the boundary plant stations contain both provincial and local dispatching equipment, during the setting calculation, the provincial dispatching only draws the provincial dispatching model, and the local dispatching model is replaced by an equivalent value; the local dispatching only draws the local dispatching model, and the provincial dispatching model is replaced by an equivalent value, so these plants are called boundary plants; adjacent plants: because the line has two sides, one side is at the R plant station and the other side is at the T plant station, then the T plant station is the adjacent plant station of the R plant station, and the boundary plant station can also be understood as the root node of the tree, the line between the two stations is understood as the branch, and the adjacent plant stations are the tree sub-nodes connected to the tree nodes. These adjacent plant stations are the first-level nodes, and then take the first-level node as the starting point to find the station connected to it, and the station found is the second-level node, and so on.

[0040] Specifically, the new energy components are photovoltaic components or wind turbine components.

[0041] Specifically, S3 specifically includes the following steps: traverse all power plants and stations, mark the bus of the power plants and stations, if the power plants and stations contain new energy components, mark the bus of the power plants and stations as a new energy bus; if the power plants and stations do not contain new energy components, mark the bus of the power plants and stations as a conventional power bus, and record the conventional power bus and the new energy bus in a conventional power bus set A and a new energy bus set B, respectively, wherein the conventional power bus set A contains N conventional power buses, and the new energy bus set B contains M new energy buses; based on the result of the marking process, a mathematical model is generated in combination with the topological connection relationship between the boundary power plants and stations, adjacent power plants and stations and the ground-dispatching power grid, as well as the impedance information and operating status of the ground-dispatching equipment.

[0042] Specifically, when the bus to be calculated is a conventional power bus, the following steps are also included: setting r=0, r represents the position of the current conventional power bus in the conventional power bus set A; when the bus to be calculated is a new energy bus, the following steps are also included: setting s=0, s represents the position of the current new energy bus in the new energy bus set B.

[0043] For details, see Figure 2, the equivalent model of each conventional power bus fault node is calculated in turn, specifically including the following steps: let r = r + 1, make the current conventional power bus operate in the maximum operation mode, based on the mathematical model of S3, calculate the first equivalent information of the current conventional power bus fault node f, including the positive sequence impedance Z f1max and zero sequence impedance Z f0max ; Make the current conventional power bus in the minimum operating mode, based on the mathematical model of S3, calculate the second equivalent information of the current conventional power bus fault node f, including the positive sequence impedance Z f1min and zero sequence impedance Z f0min ; Determine whether r is equal to N. If it is not, execute S5 again; if it is yes, report all the calculated first equivalent information and second equivalent information as the equivalent model of the conventional power bus to the provincial dispatching department; Among them, the maximum operating mode refers to an operating mode in which the current conventional power bus has the smallest short-circuit impedance value when operating in this operating mode, and the short-circuit current generated after a short circuit occurs is the largest; the minimum operating mode refers to an operating mode in which the current conventional power bus has the largest short-circuit impedance when operating in this operating mode, and the short-circuit current generated after a short circuit occurs is the smallest.

[0044] For details, see Figure 3 , iterative calculation of the external characteristic equivalent model of each new energy bus fault node specifically includes the following steps: S41, let s=s+1, i=1, initialize the transition resistance , set the upper limit of transition resistance to R lim , set the transition resistance increment to , where i represents the number of iterations of the transition resistance; k is set to 0, k represents the number of times the new energy fault current is injected into the current new energy bus; S42, let k = k + 1, obtain the terminal voltage, inject the kth new energy fault current into the current new energy bus, and calculate the first current parameters of all new energy nodes p associated with the plant corresponding to the current new energy bus in the mathematical model based on the terminal voltage, including the positive sequence current and negative sequence current ; S43, obtain the positive sequence impedance, negative sequence impedance and zero sequence impedance of the conventional power supply equivalent of the current new energy bus, where the positive sequence impedance and negative sequence impedance are both Z f1 , zero sequence impedance is Z f0 , according to the first current parameter, calculate the second current parameter of the current new energy bus fault node f, including the positive sequence current and negative sequence current ; S44, based on the first current parameter and the second current parameter, calculate the first voltage parameters of all new energy nodes p associated with the plant corresponding to the current new energy bus in the mathematical model, including the positive sequence voltage , negative sequence voltage and zero sequence voltage ; Determine whether the difference between the current first voltage parameter and the last first voltage parameter meets the convergence accuracy. When k=1, the terminal voltage is the last first voltage parameter. If the difference meets the convergence accuracy, execute S45; if the difference does not meet the convergence accuracy, return to S42; S45, let i=i+1, current transition resistance , is the last transition resistance, judge Is it greater than or equal to the upper limit of transition resistance R lim If the judgment is no, then return to S42; if the judgment is yes, then enter S46; S46, obtaining the second current parameter of the current fault node f, and calculating the second voltage parameter of the current fault node f according to the second current parameter, including the positive sequence voltage , negative sequence voltage , based on the second current parameter and the second voltage parameter, generate a fault amplitude curve and a phase angle curve; determine whether s is equal to M, if not, return to S41; if yes, report all fault amplitude curves and phase angle curves to the provincial dispatching department as the external characteristic equivalent model of the new energy bus.

[0045] Specifically, the first current parameter is:

[0046] Where j is the symbol of the imaginary part of the complex number; , , are the steady-state positive-sequence short-circuit current, active component and reactive component corresponding to the new energy node p in the kth cycle under the i-th transition resistance; , , They are respectively the steady-state negative-sequence short-circuit current amplitude, active component and reactive component of the plant corresponding to the new energy node p after the current new energy bus fails; the superscript i represents the number of iterations of the i-th transition resistance; the superscript k represents the number of times the new energy fault current is injected into the current new energy bus; the subscript d represents the active component of the current; the subscript q represents the reactive component of the current; the subscript 1 represents the positive-sequence current; the subscript 2 represents the negative-sequence current; To reflect the new energy current and positive sequence voltage controlled by its terminal voltage Function of To reflect the new energy current and negative sequence voltage controlled by its terminal voltage Function of The second current parameter is specifically: If the fault is a single-phase grounding fault, then:

[0047] in, is the positive sequence short-circuit current of the fault node in the kth cycle of the new energy under the i-th transition resistance, is the negative sequence short-circuit current of the fault node in the kth cycle of the new energy under the i-th transition resistance, is the positive sequence impedance of the fault node f, is the zero-sequence impedance of the fault node f, is the grounding resistance of the single-phase short-circuit fault in the i-th cycle, is the positive sequence mutual impedance between the fault node f and the new energy node p; If the fault is a two-phase fault, then:

[0048] If the fault is a two-phase grounding fault, then:

[0049]

[0050] Among them, / / represents AND operation; If the fault is a three-phase phase-to-phase fault, the positive sequence current and negative sequence current of the fault node are:

[0051] .

[0052] Specifically, the first voltage parameter is:

[0053]

[0054]

[0055] in, is the positive sequence voltage of the new energy node p in the kth cycle of the new energy under the i-th transition resistance, is the negative sequence voltage of the new energy node p in the kth cycle of the new energy under the i-th transition resistance, is the zero-sequence voltage of the new energy node p in the kth cycle of the new energy under the i-th transition resistance, , are the positive-sequence and negative-sequence self-impedances of the new energy node p, respectively. is the negative sequence short-circuit current provided by the new energy node p for the kth time in the i-th cycle, is the zero-sequence short-circuit current provided by the fault node f in the kth cycle of the new energy under the i-th transition resistance, is the negative sequence mutual impedance between the fault node f and the new energy node p, is the zero-sequence mutual impedance between the fault node f and the new energy node p; The second voltage parameter is specifically:

[0056]

[0057] in, , Faulty nodes The positive and negative sequence self-impedance.

[0058] Based on the same inventive concept, an embodiment of the present application provides a terminal device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the steps of the external characteristic aggregation equivalent interaction method of the new energy collection area are implemented, wherein the memory may include a memory, such as a high-speed random access memory, and may also include a non-volatile memory, such as at least one disk memory, etc.; the processor, the network interface, and the memory are interconnected through an internal bus, and the internal bus may be an industrial standard architecture bus, a peripheral component interconnection standard bus, an extended industrial standard structure bus, etc. The bus can be divided into an address bus, a data bus, a control bus, etc. The memory is used to store programs. Specifically, the program may include a program code, and the program code includes computer operation instructions. The memory may include a memory and a non-volatile memory, and provide instructions and data to the processor.

[0059] Based on the same inventive concept, an embodiment of the present application provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, it implements the steps of the external characteristic aggregation equivalent interaction method of the new energy aggregation area. Specifically, the computer-readable storage medium includes, but is not limited to, for example, volatile memory and / or non-volatile memory. The volatile memory may include RAM (Random Access Memory) and / or cache memory, etc. The non-volatile memory may include ROM (Read-Only Memory), hard disk, flash memory, optical disk, magnetic disk, etc.

[0060] It should be understood by those skilled in the art that the embodiments of the present invention may be provided as methods or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM (Compact Disc Read-Only Memory), optical storage, etc.) containing computer-usable program codes.

[0061] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer device or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0062] These computer program instructions may also be stored in a computer readable memory capable of directing a computer device or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer readable memory produce a manufactured product including an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0063] These computer program instructions can also be loaded onto a computer device or other programmable data processing device so that a series of operation steps are executed on the computer device or other programmable device to produce a process implemented by the computer device, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0064] Although the preferred embodiments of the present invention have been described, those skilled in the art may make other changes and modifications to these embodiments once they have learned the basic creative concept. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments and all changes and modifications that fall within the scope of the present invention.

[0065] Obviously, those skilled in the art can make various changes and modifications to the present invention without departing from the spirit and scope of the present invention. Thus, if these modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include these modifications and variations.

Claims

1. A method for external characteristic aggregation and equivalent interaction of a new energy gathering area, characterized in that: The following steps are involved: S1. Find all boundary plants and stations based on the voltage level of the local power grid; S2, starting from the boundary plant and station, search for adjacent plants and stations along the voltage level of the local adjustment grid; S3. Based on whether the plant station contains new energy components, the busbar of the plant station is marked. If the plant station contains new energy components, the busbar of the plant station is marked as a new energy busbar; if the plant station does not contain new energy components, the busbar of the plant station is marked as a conventional power busbar, and the plant station includes a boundary plant station and an adjacent plant station; based on the result of the marking process, a mathematical model is generated; S4. When the bus to be calculated is a conventional power bus, according to the generated mathematical model, the equivalent model of each conventional power bus fault node is calculated in turn, and reported to the provincial dispatching; when the bus to be calculated is a new energy bus, according to the generated mathematical model, by setting the transition resistance, the new energy fault current injected into the new energy bus is changed, and the external characteristic equivalent model of each new energy bus fault node is iteratively calculated, and reported to the provincial dispatching; S5. The provincial dispatching agency aggregates the equivalent model of the conventional power bus and the external characteristic equivalent model of the new energy bus to obtain an external characteristic aggregated equivalent model. The provincial dispatching agency interactively controls the new energy gathering area based on the external characteristic aggregated equivalent model.

2. According to claim 1, a method for external characteristic aggregation equivalent interaction of a new energy gathering area is characterized in that: The new energy component is a photovoltaic element or a wind turbine element.

3. According to claim 1, a new energy gathering area external characteristic aggregation equivalent interaction method is characterized in that: S3 specifically includes the following steps: traverse all power plants and stations, mark the bus of the power plants and stations, if the power plants and stations contain new energy components, mark the bus of the power plants and stations as a new energy bus; if the power plants and stations do not contain new energy components, mark the bus of the power plants and stations as a conventional power bus, and record the conventional power bus and the new energy bus in a conventional power bus set A and a new energy bus set B, respectively, wherein the conventional power bus set A contains N conventional power buses, and the new energy bus set B contains M new energy buses; based on the result of the marking process, a mathematical model is generated in combination with the topological connection relationship between the boundary power plants and stations, adjacent power plants and stations and the ground-dispatching power grid, as well as the impedance information and operating status of the ground-dispatching equipment.

4. According to claim 3, a new energy gathering area external characteristic aggregation equivalent interaction method is characterized in that: When the bus to be calculated is a conventional power bus, the following steps are also included: setting r=0, where r represents the position of the current conventional power bus in the conventional power bus set A; when the bus to be calculated is a new energy bus, the following steps are also included: setting s=0, where s represents the position of the current new energy bus in the new energy bus set B.

5. According to claim 4, a method for external characteristic aggregation equivalent interaction of a new energy gathering area is characterized in that: The equivalent model of each conventional power bus fault node is calculated in turn, which specifically includes the following steps: let r = r + 1, so that the current conventional power bus operates in the maximum operation mode, and based on the mathematical model of S3, calculate the first equivalent information of the current conventional power bus fault node f, including the positive sequence impedance Z f1max and zero sequence impedance Z f0max ; Make the current conventional power bus in the minimum operating mode, based on the mathematical model of S3, calculate the second equivalent information of the current conventional power bus fault node f, including the positive sequence impedance Z f1min and zero sequence impedance Z f0min ; Determine whether r is equal to N. If it is not, execute S5 again; if it is yes, report all the calculated first equivalent information and second equivalent information as the equivalent model of the conventional power bus to the provincial dispatching department; Among them, the maximum operating mode refers to an operating mode in which the current conventional power bus has the smallest short-circuit impedance value when operating in this operating mode, and the short-circuit current generated after a short circuit occurs is the largest; the minimum operating mode refers to an operating mode in which the current conventional power bus has the largest short-circuit impedance when operating in this operating mode, and the short-circuit current generated after a short circuit occurs is the smallest.

6. The method for external characteristic aggregation equivalent interaction of a new energy gathering area according to claim 5 is characterized in that: The iterative calculation of the external characteristic equivalent model of each new energy bus fault node specifically includes the following steps: S41, let s=s+1, i=1, initialize the transition resistance , set the upper limit of transition resistance to R lim , set the transition resistance increment to , where i represents the number of iterations of the transition resistance; k is set to 0, k represents the number of times the new energy fault current is injected into the current new energy bus; S42, let k = k + 1, obtain the terminal voltage, inject the kth new energy fault current into the current new energy bus, and calculate the first current parameters of all new energy nodes p associated with the plant corresponding to the current new energy bus in the mathematical model based on the terminal voltage, including the positive sequence current and negative sequence current ; S43, obtain the positive sequence impedance, negative sequence impedance and zero sequence impedance of the conventional power supply equivalent of the current new energy bus, where the positive sequence impedance and negative sequence impedance are both Z f1 , zero sequence impedance is Z f0 , according to the first current parameter, calculate the second current parameter of the current new energy bus fault node f, including the positive sequence current and negative sequence current ; S44, based on the first current parameter and the second current parameter, calculate the first voltage parameters of all new energy nodes p associated with the plant corresponding to the current new energy bus in the mathematical model, including the positive sequence voltage , negative sequence voltage and zero sequence voltage ; Determine whether the difference between the current first voltage parameter and the last first voltage parameter meets the convergence accuracy. When k=1, the terminal voltage is the last first voltage parameter. If the difference meets the convergence accuracy, execute S45; if the difference does not meet the convergence accuracy, return to S42; S45, let i=i+1, current transition resistance , is the last transition resistance, judge Is it greater than or equal to the upper limit of transition resistance R lim If the judgment is no, then return to S42; if the judgment is yes, then enter S46; S46, obtaining the second current parameter of the current fault node f, and calculating the second voltage parameter of the current fault node f according to the second current parameter, including the positive sequence voltage , negative sequence voltage , based on the second current parameter and the second voltage parameter, generate a fault amplitude curve and a phase angle curve; determine whether s is equal to M, if not, return to S41; if yes, report all fault amplitude curves and phase angle curves to the provincial dispatching department as the external characteristic equivalent model of the new energy bus.

7. The method for external characteristic aggregation equivalent interaction of a new energy gathering area according to claim 6 is characterized in that: The first current parameter is specifically: Where j is the symbol of the imaginary part of the complex number; , , are the steady-state positive-sequence short-circuit current, active component and reactive component corresponding to the new energy node p in the kth cycle under the i-th transition resistance; , , They are respectively the steady-state negative-sequence short-circuit current amplitude, active component and reactive component of the plant corresponding to the new energy node p after the current new energy bus fails; the superscript i represents the number of iterations of the i-th transition resistance; the superscript k represents the number of times the new energy fault current is injected into the current new energy bus; the subscript d represents the active component of the current; the subscript q represents the reactive component of the current; the subscript 1 represents the positive-sequence current; the subscript 2 represents the negative-sequence current; To reflect the new energy current and positive sequence voltage controlled by its terminal voltage Function of To reflect the new energy current and negative sequence voltage controlled by its terminal voltage Function of The second current parameter is specifically: If the fault is a single-phase grounding fault, then: in, is the positive sequence short-circuit current of the fault node in the kth cycle of the new energy under the i-th transition resistance, is the negative sequence short-circuit current of the fault node in the kth cycle of the new energy under the i-th transition resistance, is the positive sequence impedance of the fault node f, is the zero-sequence impedance of the fault node f, is the grounding resistance of the single-phase short-circuit fault in the i-th cycle, is the positive sequence mutual impedance between the fault node f and the new energy node p; If the fault is a two-phase fault, then: If the fault is a two-phase grounding fault, then: Among them, / / represents AND operation; If the fault is a three-phase phase-to-phase fault, the positive sequence current and negative sequence current of the fault node are: 。 8. The method for external characteristic aggregation equivalent interaction of a new energy gathering area according to claim 6 is characterized in that: The first voltage parameter is specifically: in, is the positive sequence voltage of the new energy node p in the kth cycle of the new energy under the i-th transition resistance, is the negative sequence voltage of the new energy node p in the kth cycle of the new energy under the i-th transition resistance, is the zero-sequence voltage of the new energy node p in the kth cycle of the new energy under the i-th transition resistance, , are the positive-sequence and negative-sequence self-impedances of the new energy node p, respectively. is the negative sequence short-circuit current provided by the new energy node p for the kth time in the i-th cycle, is the zero-sequence short-circuit current provided by the fault node f in the kth cycle of the new energy under the i-th transition resistance, is the negative sequence mutual impedance between the fault node f and the new energy node p, is the zero-sequence mutual impedance between the fault node f and the new energy node p; The second voltage parameter is specifically: in, , Faulty nodes The positive and negative sequence self-impedance.

9. A terminal device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the steps of the method according to any one of claims 1 to 8 are implemented.

10. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 8 are implemented.

Images