STATIC VOLTAGE STABILITY MARGIN EVALUATION METHOD AND SYSTEM, AND TERMINAL DEVICE

20250094663 ยท 2025-03-20

    Inventors

    Cpc classification

    International classification

    Abstract

    A static voltage stability margin evaluation method and system, and a terminal device are related to the field of integrated energy system operation. The method includes the following steps: establishing a thermal dynamic model of a heating system; establishing a thermoelectric coupling device model; establishing a static voltage stability margin model of an electric power system that considers thermal dynamics of the heating system; and solving the model to obtain a voltage stability margin. In the present invention, a static voltage stability margin that considers thermal dynamics of a heating system can be obtained, and a Pareto boundary of the static voltage stability margin that considers the thermal dynamics can be obtained through a dual-objective nonlinear optimization method, so that an impact of thermoelectric coupling on voltage stability and an impact of thermal inertia of the heating system on a voltage stability margin can be revealed.

    Claims

    1. A static voltage stability margin evaluation method, comprising: establishing a thermal dynamic model of a heating system based on temperature transmission of a heating network pipe and thermal dynamics of a building; establishing a thermoelectric coupling device model based on operation and a coupling constraint of a thermoelectric coupling device; establishing, based on an electric load growth mode and a thermal load growth mode, a static voltage stability margin model:
    max(.sub.E,.sub.H)
    s.t. g.sub.E(P,Q,V,,.sub.E)0
    g.sub.H(H,T,m,.sub.H)0
    g.sub.EH(P,H)0 wherein .sub.E and .sub.H are respectively an electric load coefficient and a thermal load coefficient; P and Q are respectively active/reactive power of an electric generator or a branch; V and are respectively a voltage and a phase of a bus; H is thermal energy; T is a water temperature or a room temperature; m is a mass flow rate; and g.sub.E()0, g.sub.H()0, and g.sub.EH()0 are respectively constraints of an electric power system, the heating system, and a thermoelectric coupling relationship; and solving the static voltage stability margin model to obtain a voltage stability margin.

    2. The static voltage stability margin evaluation method according to claim 1, wherein the establishing, based on an electric load growth mode and a thermal load growth mode, a static voltage stability margin model further comprises the following steps: establishing an electric load growth mode model: P L , c 1 i , t = ( 1 + E ) P L .Math. c 0 i , t , Q L , c 1 i , t = ( 1 + E ) Q L .Math. c 0 i , t wherein P.sub.L,c.sub.0.sup.i,t and P.sub.L,c.sub.1.sup.i,t are respectively active power of a current operating point and a safety limit point; Q.sub.L,c.sub.0.sup.i,t and Q.sub.L,c.sub.1.sup.i,t are respectively reactive power of the current operating point and the safety limit point; and .sub.E is the electric load coefficient; establishing a thermal load growth mode model: H L , c 1 k , t = ( 1 + H ) H L .Math. c 0 k , t wherein H.sub.L,c.sub.0.sup.k,t and H.sub.L,c.sub.1.sup.k,t are respectively thermal loads at the current operating point and the safety limit point; and .sub.H is the thermal load coefficient; establishing a hydraulic regulation strategy model: m p , c 1 j = { ( 1 + H ) m p , c 0 j , if ( 1 + H ) m p , c 0 j A p j v _ p , j A p j v _ p else . wherein m.sub.p,c.sub.0.sup.j and m.sub.p,c.sub.1.sup.j are respectively mass flow rates of a pipe j at the current operating point and the safety limit point; and v.sub.p is a water velocity upper limit to avoid hydraulic instability and pipe erosion; and establishing the static voltage stability margin model based on the electric load growth mode model, the thermal load growth mode model, and the hydraulic regulation strategy model.

    3. The static voltage stability margin evaluation method according to claim 2, wherein the establishing a thermal dynamic model of a heating system comprises the following steps: establishing a temperature transmission model of the heating network pipe: T out j , t = p j ( k p j T in j , t - p j + ( 1 - k p j ) T in j , t - p j - 1 ) + ( 1 - p j ) T amb t wherein t is a scheduled period set; j is a pipe set; parameters .sub.p.sup.j and k.sub.p.sup.j are coefficients of the pipe j that are related to a transmission delay; .sub.p.sup.j is a heat preservation coefficient; T.sub.amb.sup.t is a pipe ambient temperature at a moment t; T.sub.in.sup.j,t and T.sub.out.sup.j,t are respectively temperatures of a heat medium at a pipe inlet and outlet at the moment t; and calculation formulas of the parameters .sub.p.sup.j, k.sub.p.sup.j, and .sub.p.sup.j are as follows: { p j = .Math. w A p j L p j / ( m p j t ) .Math. - 1 k p j = p j + 1 - w A p j L p j / ( m p j t h ) p j = 1 - exp ( - p j t / ( w c w A p j ) ( p j + 3 / 2 - k p j ) ) wherein t is a time interval; .sub.w is a density of the heat medium; c.sub.w is a specific heat capacity of the heat medium; A.sub.p.sup.j is a cross-sectional area of the pipe j; L.sub.p.sup.j is a length of the pipe j; .sub.p.sup.j is a pipe heat loss coefficient; m.sub.p.sup.j is a mass flow rate of the heat medium of the pipe; and is a round-up function; and establishing a thermal dynamic model of the building: T b k , t = b k T b k , t - 1 + b k H L k , t + b k T out t wherein k is a building set; and .sub.b.sup.k, .sub.b.sup.k and .sub.b.sup.k are parameters that depend on a heat capacity and thermal resistance of the building.

    4. The static voltage stability margin evaluation method according to claim 1, wherein the solving the static voltage stability margin model to obtain a voltage stability margin comprises the following steps: optimizing the static voltage stability margin model as:
    .sub.E()=max .sub.E
    s.t. .sub.H=
    g.sub.E0,g.sub.H0,g.sub.EH0 wherein is a constant for predefining .sub.H; and initializing =0, Flag=1, k=0, and , and performing solving according to the following cycle step, wherein the cycle step comprises: updating .sub.H.sup.k, and calculating a mass flow rate m.sub.p,c.sub.1.sup.j, and parameters .sub.p.sup.j, k.sub.p.sup.j, and .sub.p.sup.j; solving the optimized static voltage stability margin model; and if there is a solution, updating .sub.E.sup.k.sub.E(.sub.H.sup.k); otherwise, updating Flag0; updating kk+1, +; and if Flag=1, repeating the cycle step; or if Flag=, exiting the cycle step.

    5. The static voltage stability margin evaluation method according to claim 1, wherein the thermoelectric coupling device model comprises an operation constraint of a combined heat and power unit, a coupling constraint between the combined heat and power unit and a heating network, and a water temperature constraint of the heating system.

    6. The static voltage stability margin evaluation method according to claim 5, wherein the coupling constraint between the combined heat and power unit and the heating network is as follows: H chp i , t = c w m src i , t ( T src , s i , t - T src , r i , t ) wherein m.sub.src.sup.i,t is a mass flow rate of a heat source; and T.sub.src,s.sup.i,t and T.sub.src,r.sup.i,t, are respectively a supply water temperature and a return water temperature of the heat source.

    7. The static voltage stability margin evaluation method according to claim 5, wherein the water temperature constraint of the heating system is as follows:
    T.sub.src,s.sup.iT.sub.src,s.sup.i,tT.sub.src,s.sup.i,T.sub.src,r.sup.iT.sub.src,r.sup.i,tT.sub.src,r.sup.i where T.sub.src,s.sup.i and T.sub.src,s.sup.i are respectively a lower limit and an upper limit of a supply water temperature of a heat source; and T.sub.src,r.sup.i and T.sub.src,r.sup.i are respectively a lower limit and an upper limit of a return water temperature of the heat source.

    8. A static voltage stability margin evaluation system using the static voltage stability margin evaluation method according to claim 1, comprising: a thermal dynamic model module of a heating system, configured to establish a thermal dynamic model of the heating system based on temperature transmission of a heating network pipe and thermal dynamics of a building; a thermoelectric coupling device model module, configured to establish a thermoelectric coupling device model based on operation and a coupling constraint of a thermoelectric coupling device; a static voltage stability margin model module, configured to establish, based on an electric load growth mode and a thermal load growth mode, a static voltage stability margin model; and a calculation module, configured to solve the static voltage stability margin model based on the thermal dynamic model module of the heating system, the thermoelectric coupling device model module, and the static voltage stability margin model module, to obtain a voltage stability margin.

    9. A terminal device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, the processor, when loading and executing the computer program, using the static voltage stability margin evaluation method according to claim 1.

    10. A computer-readable storage medium, storing a computer program, the computer program, when being loaded and executed by the processor, using the static voltage stability margin evaluation method according to claim 1.

    Description

    BRIEF DESCRIPTION OF THE DRAWINGS

    [0046] The following further describes the present invention in detail with reference to the accompanying drawings.

    [0047] FIG. 1 is a structural diagram of an integrated energy system according to the present invention;

    [0048] FIG. 2 is a flowchart of a static voltage stability margin evaluation method that considers thermal dynamics of a heating system according to the present invention;

    [0049] FIG. 3 is a structural diagram of an integrated energy system according to Embodiment 2 of the present invention;

    [0050] FIG. 4 shows a Pareto frontier of a static voltage stability margin evaluation result that considers thermal dynamics of a heating system according to Embodiment 2 of the present invention;

    [0051] FIG. 5 shows an operating point of a CHP unit when .sub.H=1.5 in a static voltage stability margin evaluation result that considers thermal dynamics of a heating system according to Embodiment 2 of the present invention;

    [0052] FIG. 6 shows active power output results of a generator in four scenarios when .sub.H=1.5 according to Embodiment 2 of the present invention; and

    [0053] FIG. 7 shows thermal load results in four scenarios when .sub.H=1.5 according to Embodiment 2 of the present invention.

    DETAILED DESCRIPTION

    [0054] The technical solutions of embodiments of the present invention are clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are merely a part rather than all of the embodiments of the present invention. All other embodiments obtained by a person of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

    [0055] In the descriptions of this specification, descriptions of reference terms such as an embodiment, an example, and a specific example mean that specific features, structures, materials, or characteristics that are described with reference to the embodiment or the example are included in at least one embodiment or example of the present invention. In this specification, schematic descriptions of the foregoing terms do not necessarily point at a same embodiment or example. In addition, the described specific features, structures, materials, or characteristics may be combined in a proper manner in any one or more of the embodiments or examples.

    Embodiment 1

    [0056] This embodiment is used in an integrated energy system, and a structure of the integrated energy system is shown in FIG. 1.

    [0057] As shown in FIG. 2, a static voltage stability margin evaluation method that considers thermal dynamics of a heating system includes the following steps.

    [0058] S1: Establish a thermal dynamic model of a heating system.

    [0059] S11: Establish a temperature transmission model of a heating network pipe, where the temperature transmission model of the heating network pipe is as follows:

    [00008] T out j , t = p j ( k p j T in j , t - p j + ( 1 - k p j ) T in j , t - p j - 1 ) + ( 1 - p j ) T amb t [0060] where t is a scheduled period set; j is a pipe set; parameters .sub.p.sup.j and k.sub.p.sup.j are coefficients of the pipe j that are related to a transmission delay; .sub.p.sup.j is a heat preservation coefficient; T.sub.amb.sup.t is a pipe ambient temperature at a moment t; T.sub.in.sup.j,t and T.sub.out.sup.j,t are respectively temperatures of a heat medium at a pipe inlet and outlet at the moment t; and [0061] calculation formulas of the parameters .sub.p.sup.j, k.sub.p.sup.j, and .sub.p.sup.j are as follows:

    [00009] { p j = .Math. w A p j L p j / ( m p j t ) .Math. - 1 k p j = p j + 1 - w A p j L p j / ( m p j t h ) p j = 1 - exp ( - p j t / ( w c w A p j ) ( p j + 3 / 2 - k p j ) ) [0062] where t is a time interval; .sub.w is a density of the heat medium; c.sub.w is a specific heat capacity of the heat medium; A.sub.p.sup.j is a cross-sectional area of the pipe j; L.sub.p.sup.j is a length of the pipe j; .sub.p.sup.j is a pipe heat loss coefficient; m.sub.p.sup.j is a mass flow rate of the heat medium of the pipe; and is a round-up function.

    [0063] S12: Establish a thermal dynamic model of a building:

    [00010] T b k , t = b k T b k , t - 1 + b k H L k , t + b k T out t [0064] where k is a building set; and .sub.b.sup.k, .sub.b.sup.k and .sub.b.sup.k are parameters that depend on a heat capacity and thermal resistance of the building.

    [0065] S2: Establish a thermoelectric coupling device model.

    [0066] S21: Establish an operation constraint of a combined heat and power unit.

    [0067] S211: Establish an operation constraint of a condensing combined heat and power unit:

    [00011] a chp i , j P chp i , t + b chp i , j H chp i , t c chp i , j j J [0068] where i is a unit set; j is a set of limit points of a feasible region; a.sub.chp.sup.i,j, b.sub.chp.sup.i,j, and c.sub.chp.sup.i,j are unit parameters; and P.sub.chp.sup.i,t and H.sub.chp.sup.i,t are respectively electric power and thermal power output by the CHP unit.

    [0069] S212: Establish an operation constraint of a back pressure combined heat and power unit:


    P.sub.chp.sup.i,t=k.sub.chp.sup.i,tH.sub.chp.sup.i,t [0070] where k.sub.chp.sup.i,j is a heat-power ratio of the CHP unit.

    [0071] S22: Establish a coupling constraint between the combined heat and power unit and a heating network:

    [00012] H chp i , t = c w m src i , t ( T src , s i , t - T src , r i , t ) [0072] where m.sub.src.sup.i,t is a mass flow rate of a heat source; and T.sub.src,s.sup.i,t and T.sub.src,r.sup.i,t are respectively a supply water temperature and a return water temperature of the heat source.

    [0073] S23: Establish a water temperature constraint of the heating system:


    T.sub.src,s.sup.iT.sub.src,s.sup.i,tT.sub.src,s.sup.i,T.sub.src,r.sup.iT.sub.src,r.sup.i,tT.sub.src,r.sup.i [0074] where T.sub.src,s.sup.i and T.sub.src,s.sup.i (or T.sub.src,r.sup.i and T.sub.src,r.sup.i) are respectively a lower limit and an upper limit of a supply water temperature (or return water temperature) of a heat source.

    [0075] S3: Establish a static voltage stability margin model that considers thermal dynamics of the heating system.

    [0076] S31: Establish a load growth mode model.

    [0077] S311: Establish an electric load growth mode model:

    [00013] P L , c 1 i , t = ( 1 + E ) P L .Math. c 0 i , t , Q L , c 1 i , t = ( 1 + E ) Q L .Math. c 0 i , t [0078] where P.sub.L,c.sub.0.sup.i,t and P.sub.L,c.sub.1.sup.i,t (or Q.sub.L,c.sub.0.sup.i,t and Q.sub.L,c.sub.1.sup.i,t) are respectively active (or reactive) power of a current operating point and a safety limit point; and .sub.E is an electric load coefficient.

    [0079] S312: Establish a thermal load growth mode model based on a current operating point:

    [00014] H L , c 1 k , t = ( 1 + H ) H L .Math. c 0 k , t [0080] where H.sub.L,c.sub.0.sup.k,t and H.sub.L,c.sub.1.sup.k,t are respectively thermal loads at the current operating point and the safety limit point; and .sub.H is a thermal load coefficient.

    [0081] S313: Establish a hydraulic regulation strategy model:

    [00015] m p , c 1 j = { ( 1 + H ) m p , c 0 j , if ( 1 + H ) m p , c 0 j A p j v _ p , j A p j v _ p else . [0082] where m.sub.p,c.sub.0.sup.j and m.sub.p,c.sub.1.sup.j are respectively mass flow rates of a pipe j at the current operating point and the safety limit point; and v.sub.p is a water velocity upper limit to avoid hydraulic instability and pipe erosion.

    [0083] S32: Establish the static voltage stability margin model that considers the thermal dynamics of the heating system, where a mathematical form thereof is shown in the following formula:


    max(.sub.E,.sub.H)


    s.t. g.sub.E(P,Q,V,,.sub.E)0


    g.sub.H(H,T,m,.sub.H)0


    g.sub.EH(P,H)0 [0084] where .sub.E and .sub.H are respectively an electric load coefficient and a thermal load coefficient; P and Q are respectively active/reactive power of an electric generator or a branch; V and are respectively a voltage and a phase of a bus; H is thermal energy; T is a water temperature or a room temperature; m is a mass flow rate; and g.sub.E()0, g.sub.H()0, and g.sub.EH()0 are respectively constraints of an electric power system, the heating system, and a thermoelectric coupling relationship; and

    [0085] S4: Solve the model to obtain a voltage stability margin.

    [0086] S41: Convert a dual-objective optimization model established in S32 into:


    .sub.E()=max .sub.E


    s.t. .sub.H=


    g.sub.E0,g.sub.H0,g.sub.EH0 [0087] where is a constant for predefining .sub.H.

    [0088] S42: Solve the model established in S41 by using the following algorithm:

    [0089] S421: Initialize =0, Flag=1, k=0, and .

    [0090] S422: Update .sub.H.sup.k, and calculate a mass flow rate m.sub.p,c.sub.1.sup.j, and parameters .sub.p.sup.j, k.sub.p.sup.j, and .sub.p.sup.j.

    [0091] S423: Solve the model in step S41; and if there is one solution, update .sub.E.sup.k.sub.E(.sub.H.sup.k); otherwise, update Flag0.

    [0092] S424: Update kk+1, +.

    [0093] S425: If Flag=1, repeat steps S422 to S424; or if Flag=0, exit the cycle.

    Embodiment 2

    [0094] An integrated energy system in this embodiment is formed by a 9-node electric power system and a 12-node heating system. As shown in FIG. 3, the system includes two 300 MW coal-fired units, one 300 MW CHP unit, and one gas boiler with a capacity of 300 MW, a scheduled time interval is 1 h, total scheduled time is 24 h, a supply temperature of a heat source is set to constant (90 C.), and a room temperature is set to an ideal value (20 C.). A step size of .sub.H, namely, , is set to 0.05.

    [0095] According to the steps of the present invention, a static voltage stability margin that considers thermal dynamics of a heating system is solved. A Pareto frontier obtained by solving is shown in FIG. 4. Distribution of operating points of the CHP unit when .sub.H=1.5 in a static voltage stability margin result that considers the thermal dynamics of the heating system is shown in FIG. 5. Active power output results and thermal load results of a generator in four scenarios when .sub.H=1.5 are shown in FIG. 6 and FIG. 7 respectively. It can be seen that, a thermoelectric coupling relationship in an integrated thermoelectric energy system has a negative impact on voltage stability, but thermal inertia of a heating system (especially a building) can be used to improve an operating point of a combined heat and power apparatus, thereby improving a voltage stability margin of an electric power system.

    [0096] Therefore, in the method, an electric power network constraint, a heating network constraint, and a building thermal load constraint of the integrated energy system can be comprehensively considered, thermal dynamic characteristics, thermoelectric coupling, and a voltage stability margin of the integrated thermoelectric energy system are accurately characterized, and interaction of an electric power network and a thermal network of the integrated energy system can be effectively characterized, to implement evaluation on a static voltage stability margin in the integrated energy system.

    [0097] The embodiments of this application disclose a static voltage stability margin evaluation system, including: [0098] a thermal dynamic model module of a heating system, configured to establish a thermal dynamic model of the heating system based on temperature transmission of a heating network pipe and thermal dynamics of a building; [0099] a thermoelectric coupling device model module, configured to establish a thermoelectric coupling device model based on operation and a coupling constraint of a thermoelectric coupling device; [0100] a static voltage stability margin model module, configured to establish, based on an electric load growth mode and a thermal load growth mode, a static voltage stability margin model; and [0101] a calculation module, configured to solve the static voltage stability margin model based on the thermal dynamic model module of the heating system, the thermoelectric coupling device model module, and the static voltage stability margin model module, to obtain a voltage stability margin.

    [0102] The embodiments of this application further disclose a terminal device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, the processor, when executing the computer program, using any static voltage stability margin evaluation method in the foregoing embodiments.

    [0103] The terminal device may be a computer device such as a desktop computer, a laptop computer, or a cloud server, and the terminal device includes but is not limited to a processor and a memory. For example, the terminal device may alternatively include input and output devices, a network access device, a bus, and the like.

    [0104] The processor may be a central processing unit (CPU). Certainly, according to an actual use situation, the processor may be another general purpose processor, a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field programmable gate array (FPGA) or another programmable logical device, a discrete gate or a transistor logical device, a discrete hardware component, or the like. The general purpose processor may be a microprocessor, any conventional processor, or the like. This is not limited in this application.

    [0105] The memory may be an internal storage unit of the terminal device, such as a hard disk or a memory of the terminal device, or an external storage device of the terminal device, such as a plug-in hard disk, a smart memory card (SMC), a secure digital card (SD), or a flash memory card (FC) equipped on the terminal device. The memory may alternatively be a combination of the internal storage unit of the terminal device and the external storage device. The memory is configured to store a computer program and other programs and data required by the terminal device. The memory may alternatively be configured to temporarily store data that has been output or that is to be output. This is not limited in this application.

    [0106] By using the terminal device, any static voltage stability margin evaluation method in the foregoing embodiments is stored in the memory of the terminal device, and is loaded and executed on the processor of the terminal device for convenient use.

    [0107] The embodiments of this application further disclose a computer-readable storage medium, storing a computer program, the computer program, when being executed by a processor, using any static voltage stability margin evaluation method in the foregoing embodiments.

    [0108] The computer program may be stored in the computer-readable medium. The computer program includes computer program code. The computer program code can be in a form of source code, object code, executable file or some middleware form, etc. The computer-readable medium may include: any entity or apparatus that is capable of carrying the computer program code, a recording medium, a USB flash drive, a removable hard disk, a magnetic disk, an optical disc, a computer memory, a read-only memory (ROM), a random access memory (RAM), an electric carrier signal, a telecommunication signal and a software distribution medium, or the like. It should be noted that, the computer-readable medium includes but is not limited to the foregoing components.

    [0109] By using the computer-readable storage medium, any static voltage stability margin evaluation method in the foregoing embodiments is stored in the computer-readable storage medium, and is loaded and executed on the processor for convenient storage and use of the foregoing method.

    [0110] The foregoing shows and describes basic principles, main features of the present invention and advantages of the present invention. A person skilled in the art may understand that the present invention is not limited to the foregoing embodiments. Descriptions in the embodiments and this specification only illustrate the principles of the present invention. Various modifications and improvements are made in the present invention without departing from the spirit and the scope of the present invention, and these modifications and improvements shall fall within the protection scope of the present invention.