Scheduling method and system considering building temperature control load aggregator and privacy security thereof

Abstract

The invention discloses a centralized scheduling method and system considering building temperature control load aggregators and privacy security thereof, and relates to the field of privacy protection and demand side response. The method comprises the following steps: S1, establishing a building temperature control load aggregator model; s2, establishing a distribution network economic dispatching model considering a building temperature control load aggregator; s3, establishing a scheduling model calculation method considering privacy security, specifically, S31, implementing an encryption method based on random transformation I; s32, converting the inequality constraint into an equality constraint; s33, implementing constraint equivalent expansion; s34, implementing an encryption method based on random transformation II; and S35, economic dispatching considering the building temperature control load aggregator and the privacy security thereof is implemented. According to the method, an electric power system operator can obtain the building temperature control load aggregation model so as to make an economic dispatching plan, and meanwhile privacy information of a building temperature control load aggregator is guaranteed. The method is beneficial for promoting building heat flexibility excavation and improving operation regulation and control flexibility of a power system.

Description

Technical Field

[0001] The present invention relates to the field of power system operation optimization and privacy security, and specifically to a centralized scheduling method considering building thermostatic load aggregators and their privacy security. Background Art

[0002] With the increasing proportion of intermittent renewable energy, the power system is undergoing a major transformation. The power system requires more flexible resources to support its safe and economic operation. The inherent thermal inertia of buildings brings considerable flexibility to heating and cooling, and is considered a potential demand response resource. The concept of using building thermal inertia for power system operation and control services has attracted extensive attention in the academic and engineering fields.

[0003] Due to the large number of buildings, direct information interaction between the energy system and a vast number of building users will generate huge computing and communication burdens. To alleviate this problem, building thermostatic load aggregators are actually used as agents to participate in power system scheduling. However, the way of direct interaction between building thermostatic load aggregators and the power system will directly expose the aggregation model parameters to the power grid, which may bring privacy leakage risks and economic losses to building thermostatic load aggregators. Summary of the Invention

[0004] The technical problem to be solved by the present invention is how to enable the power system operator to obtain the building thermostatic load aggregation model on the basis of the privacy information security of the building thermostatic load aggregator in order to formulate an economic scheduling plan.

[0005] The object of the present invention can be achieved by the following technical solutions:

[0006] A centralized scheduling method considering building thermostatic load aggregators and their privacy security, comprising the following steps:

[0007] S1. Establish a building thermostatic load aggregator model;

[0008] S2. Establish a distribution network economic scheduling model considering building thermostatic loads;

[0009] S3. Establish a calculation method for a scheduling model considering privacy security, specifically including:

[0010] S31. Implement an encryption method based on random transformation I;

[0011] S32. Convert inequality constraints into equality constraints;

[0012] S33. Implement constraint equivalent extension;

[0013] S34. Implement an encryption method based on random transformation II;

[0014] S35. Implement economic dispatch considering building temperature control load aggregators and their privacy and security.

[0015] Furthermore, the step S1 of establishing a building temperature control load aggregator model includes:

[0016] S11. Establish an algebraic form of the building temperature control load aggregation model:

[0017] The building temperature control load aggregator interacts with the distribution network as an agent of building users, and its model can be established as:

[0018]

[0019] Where, represents the aggregated temperature of building temperature control load aggregator k, represents the cooling / heating power of building temperature control load aggregator k, is the model parameter of building temperature control load aggregator k, and the set M = {0, 1,..., M} represents the model order, and τ k are the upper and lower limits of the aggregated temperature of building temperature control load aggregator k respectively.

[0020] S12. Establish a compact form of the building temperature control load aggregation model:

[0021] First, define the constant matrix where m is the model order, and the element (Λ m ) i,j in the i-th row and j-th column of this matrix is:

[0022]

[0023] Secondly, use the state variable to represent the aggregated temperature of building temperature control load aggregator k, and use the control variable to represent the cooling / heating power of building temperature control load aggregator k. Define the state variable vector x k and the control variable vector u k respectively as:

[0024]

[0025] Thirdly, define the constant vector where the t-th element has the following expression:

[0026]

[0027] Define the matrix

[0028]

[0029] Among them, the constant matrix I T represents the T-dimensional identity matrix. Therefore, the building temperature control load aggregation model established by S11 can be transformed into its compact form:

[0030]

[0031] where 1 T represents the T-dimensional vector of 1.

[0032] Furthermore, the step S2 of establishing the distribution network economic dispatch model considering the building temperature control load aggregator further includes: S21, establishing the distribution network constraint conditions; S22, establishing the coupling constraint conditions between the distribution network and the building temperature control load aggregator; S23, establishing the constraint conditions for the control variables of the building temperature control load aggregator; S24, establishing the distribution network economic dispatch model considering the building temperature control load aggregator.

[0033] Furthermore, S21, establishing the distribution network constraint conditions specifically includes:

[0034] Denote the distribution network economic dispatch decision variable as the vector z (z does not include the cooling / heating power of each building temperature control load aggregator ), and denote the feasible region of z as Z. Then the distribution network constraint conditions can be obtained:

[0035] z ∈ Z

[0036] S22, establishing the coupling constraint conditions between the distribution network and the building temperature control load aggregator specifically includes:

[0037] Denote the vector formed by the control variables of each building temperature control load aggregator k. Then the coupling constraint conditions between the distribution network and the building temperature control load aggregator can be expressed as:

[0038] Az + u = 0,

[0039] where the matrix A represents the connection relationship between the distribution network bus and the building temperature control load aggregator.

[0040] S23, establishing the constraint conditions for the control variables of the building temperature control load aggregator specifically includes:

[0041] Define the control variable u k of the building temperature control load aggregator. According to S22, obtain the feasible region U k of u k as:

[0042] U k = {u k |R k xk +S k u k =d k ,

[0043]

[0044] S24. The establishment of a distribution network economic dispatch model considering building temperature control load aggregators specifically includes:

[0045] The distribution network economic dispatch model considering building temperature control load aggregators can be expressed as:

[0046]

[0047] Among them, c represents the cost vector, U k The expression is as shown in S23, that is:

[0048]

[0049] Furthermore, the calculation method for establishing a dispatch model considering privacy and security in step S3 includes:

[0050] S31. Implement an encryption method based on random transformation I

[0051] First, each building temperature control load aggregator k generates a reversible random matrix Then, the model of each building temperature control load aggregator k in S12 can be equivalently transformed into:

[0052]

[0053] S32. Transform inequality constraints into equality constraints

[0054] First, we denote Denote

[0055] , introduce slack variables Thus, the model of building temperature control load aggregator k in S31 can be equivalently transformed into:

[0056]

[0057] Among them, is a random diagonal matrix with all diagonal elements being positive and is E k is owned by building temperature control load aggregator k.

[0058] S33. Implement constraint equivalent expansion

[0059] First, copy the model of building temperature control load aggregator k in S32 as follows:

[0060]

[0061] Record Thus, the model of the building temperature control load aggregator k can be equivalently transformed into:

[0062]

[0063] S34. Implement the encryption method based on random transformation II

[0064] First, the building temperature control load aggregator k generates a reversible random matrix Then, this random matrix is applied to the model of the building temperature control load aggregator k formed in S33, and we get:

[0065]

[0066] S35. Implement the economic dispatch considering the building temperature control load aggregator and its privacy security

[0067] The building temperature control load aggregator k uploads the matrices V K F k , V k G k , V k H k and V k e k to the power system operator. Thus, the power system operator can formulate the economic dispatch problem considering the building temperature control load aggregator and its privacy security as follows:

[0068]

[0069] where, U k The expression is as follows:

[0070]

[0071] The power system operator solves the above economic dispatch problem and returns the solution result to the building temperature control load aggregator k. The building temperature control load aggregator k can restore the true state variable x according to the following formula k :

[0072]

[0073] The present invention also discloses an economic dispatch system considering the building temperature control load aggregator and its privacy security, including the following steps:

[0074] A building temperature control load aggregator establishment module, used for establishing a building temperature control load aggregator model;

[0075] The distribution network economic dispatch model establishment module is used to establish a distribution network economic dispatch model considering the building temperature control load aggregator;

[0076] The privacy protection model is used to establish a calculation method for the dispatch model considering privacy security, specifically including: the first encryption unit, which is used to implement the encryption method based on random transformation I; the constraint conversion unit, which is used to convert the inequality constraint into an equality constraint; the constraint extension unit, which is used to implement the equivalent extension of the constraint; the second encryption unit, which is used to implement the encryption method based on random transformation II; the implementation unit, which is used to implement the economic dispatch considering the building temperature control load aggregator and its privacy security.

[0077] Furthermore, the building temperature control load aggregator model establishment module includes the following units:

[0078] The algebraic form of the building temperature control load aggregation model unit is used to establish the algebraic form of the building temperature control load aggregation model:

[0079] The building temperature control load aggregator interacts with the distribution network as an agent of the building users, and its algebraic form of the building temperature control load b model is established as:

[0080]

[0081] Among them, represents the aggregated temperature of the building temperature control load aggregator k, represents the cooling / heating power of the building temperature control load aggregator k, is the model parameter of the building temperature control load aggregator k, and the set M = {0, 1,..., M} represents the model order, and τ k are respectively the upper and lower limits of the aggregated temperature of the building temperature control load aggregator k;

[0082] The compact form of the building temperature control load aggregation model establishment unit is used to establish the compact form of the building temperature control load aggregation model:

[0083] First, define the constant matrix where m is the model order, and the element (Λ m ) i,j in the i-th row and j-th column of this matrix is:

[0084]

[0085] Secondly, use the state variable to represent the aggregated temperature of the building temperature control load aggregator k, and use the control variable to represent the cooling / heating power Define the state variable vector \(x\) respectively k and the control variable vector \(u\) k as follows:

[0086]

[0087] Secondly, define the constant vector The \(t\)-th element of which has the following expression:

[0088]

[0089] Define the matrix

[0090]

[0091] Among them, the constant matrix \(I\) T represents the \(T\)-dimensional identity matrix. Therefore, the established building temperature control load aggregation model is transformed into its compact form:

[0092]

[0093] where \(1\) T represents the \(T\)-dimensional all-one vector.

[0094] Furthermore, the distribution network economic dispatch model establishment module further includes: a distribution network constraint condition establishment unit for establishing distribution network constraint conditions; a distribution network and building temperature control load aggregator coupling constraint condition establishment unit for establishing distribution network and building temperature control load aggregator coupling constraint conditions; a building temperature control load aggregator control variable constraint condition establishment unit for establishing building temperature control load aggregator control variable constraint conditions; and a distribution network economic dispatch model establishment unit for establishing a distribution network economic dispatch model considering the building temperature control load aggregator.

[0095] Furthermore, the distribution network constraint condition establishment unit specifically executes the following steps:

[0096] Denote the distribution network economic dispatch decision variable as the vector \(z\), and \(z\) does not include the cooling / heating power of each building temperature control load aggregator Denote the feasible region of \(z\) as \(Z\), then the distribution network constraint conditions are obtained:

[0097] \(z\in Z\)

[0098] The distribution network and building temperature control load aggregator coupling constraint condition establishment unit specifically executes the following steps:

[0099] Denote the control variables of each building temperature control load aggregator constituting the vector Then the distribution network and building temperature control load aggregator coupling constraint conditions are expressed as:

[0100] Az + u = 0,

[0101] where the matrix A represents the connection relationship between the distribution network bus and the building temperature control load aggregator;

[0102] The unit for establishing the control variable constraint conditions of the building temperature control load aggregator specifically executes the following steps:

[0103] Define the control variable u of the building temperature control load aggregator k , and obtain u according to step S22 k The feasible region U k is:

[0104]

[0105] The unit for establishing the distribution network economic dispatch model specifically executes the following steps:

[0106] The distribution network economic dispatch model considering the building temperature control load aggregator is expressed as:

[0107]

[0108] where c represents the cost vector, and U k The expression is as shown in S23, that is:

[0109]

[0110] Furthermore, the privacy protection model includes:

[0111] The first encryption unit, which executes the following steps:

[0112] First, each building temperature control load aggregator k generates a reversible random matrix Then each building temperature control load aggregator k's model is equivalently transformed into:

[0113]

[0114] The constraint transformation unit, which executes the following steps:

[0115] First, denote Denote

[0116] , and introduce a slack variable Thus, the model of the building temperature control load aggregator k in step S31 is equivalently transformed into:

[0117]

[0118] where is a random diagonal matrix with all diagonal elements being positive and is E k is owned by the building temperature control load aggregator k;

[0119] The constraint extension unit performs the following steps:

[0120] First, copy the model of the building temperature control load aggregator k in step S32 as follows:

[0121]

[0122] Denote Then the above model of the building temperature control load aggregator k is equivalently transformed into:

[0123]

[0124] The second encryption unit performs the following steps:

[0125] First, the building temperature control load aggregator k generates a reversible random matrix Then apply this random matrix to the model of the building temperature control load aggregator k formed in step S33 to obtain:

[0126]

[0127] The implementation unit performs the following steps:

[0128] The building temperature control load aggregator k uploads the matrices V K F k , V k G k , V k H k and V k e k to the power system operator. Then the power system operator formulates an economic dispatch problem considering the building temperature control load aggregator and its privacy and security as follows:

[0129]

[0130] where U k The expression is as follows:

[0131]

[0132] The power system operator solves the above economic dispatch problem and returns the solution result to the building temperature control load aggregator k. The building temperature control load aggregator k can restore the true state variable x according to the following formula k :

[0133]

[0134] The advantages of the present invention are:

[0135] A centralized scheduling method considering building temperature control load aggregators and their privacy and security is provided. By adopting random transformation technology and introducing random variables, the privacy information security of building temperature control load aggregators is guaranteed. At the same time, the power system operator can obtain the building temperature control load aggregation model to formulate an economic scheduling plan, thereby promoting the excavation of building thermal flexibility and improving the operation and regulation flexibility of the power system. This method has the advantages of good privacy protection performance and small computational complexity, and has wide practical engineering application value. Brief Description of the Drawings

[0136] The present invention will be further described below with reference to the accompanying drawings.

[0137] Figure 1 is a flowchart of a centralized scheduling method considering building temperature control load aggregators and their privacy and security according to the present invention;

[0138] Figure 2 is a structural diagram of a building temperature control load aggregator - building cluster system in Embodiments 1 and 2 of the present invention;

[0139] Figure 3 is a comparison chart of the cooling / heating power of Building Temperature Control Load Aggregator 1 under privacy protection and non-privacy protection calculation methods in Embodiment 2;

[0140] Figure 4 is matrix G in Embodiment 2 1 and V 1 G 1 The comparison result of the heat maps of the first 24 rows. Detailed Embodiments

[0141] The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all the embodiments. All other embodiments obtained by those 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.

[0142] Embodiment 1

[0143] This embodiment is applied to the power system operator - building system.

[0144] This embodiment provides a centralized scheduling method considering building temperature control load aggregators and their privacy and security, including the following steps:

[0145] S1. Establish a building temperature control load aggregator model, specifically including the following steps:

[0146] S11. Establish an algebraic form of the building temperature control load aggregation model:

[0147] The building temperature control load aggregator interacts with the distribution network as an agent of building users, and its model can be established as follows:

[0148]

[0149] where represents the aggregated temperature of building temperature control load aggregator k, represents the cooling / heating power of building temperature control load aggregator k, is the model parameter of building temperature control load aggregator k, m is the model order, and the set M = {0, 1, …, M} represents the model order, and τ k are the upper and lower limits of the aggregated temperature of building temperature control load aggregator k, respectively.

[0150] S12. Establish a compact form of the building temperature control load aggregation model:

[0151] First, define the constant matrix where m is the model order, and the element (Λ m ) i,j in the i-th row and j-th column of this matrix is:

[0152]

[0153] Second, use the state variable to represent the aggregated temperature of building temperature control load aggregator k, and use the control variable to represent the cooling / heating power of building temperature control load aggregator k. Define the state variable vector x k and the control variable vector u k respectively as:

[0154]

[0155] Third, define the constant vector where the t-th element has the following expression:

[0156]

[0157] Define the matrix

[0158]

[0159] where the constant matrix I T represents the T-dimensional identity matrix. Therefore, the building temperature control load aggregation model established in step S11 can be transformed into its compact form:

[0160]

[0161] Among them, 1 T represents a T-dimensional 1 vector.

[0162] S2. Establish a distribution network economic dispatch model considering building temperature control load aggregators, specifically including:

[0163] S21. Establish distribution network constraint conditions

[0164] Denote the distribution network economic dispatch decision variable as vector z (z does not include the cooling / heating power of each building temperature control load aggregator ), and denote the feasible region of z as Z. Then the distribution network constraint conditions can be obtained as:

[0165] z ∈ Z

[0166] S22. Establish coupling constraint conditions between the distribution network and building temperature control load aggregators

[0167] Denote the control variable of each building temperature control load aggregator k The coupling constraint conditions between the distribution network and building temperature control load aggregators can be expressed as:

[0168] Az + u = 0,

[0169] where the matrix A represents the connection relationship between the distribution network bus and building temperature control load aggregators.

[0170] S23. Establish control variable constraint conditions for building temperature control load aggregators

[0171] Define the control variable u of the building temperature control load aggregator k , and according to S22, obtain the feasible region U k of u k as:

[0172]

[0173] S24. Establish a distribution network economic dispatch model considering building temperature control load aggregators

[0174] The distribution network economic dispatch model considering building temperature control load aggregators can be expressed as:

[0175]

[0176] Among them, c represents the cost vector, and U k is expressed as shown in S23, that is:

[0177]

[0178] S3. Establish a scheduling model calculation method considering privacy and security, specifically including

[0179] S31. Implement an encryption method based on random transformation I

[0180] First, each building temperature control load aggregator k generates a reversible random matrix Then, the model of each building temperature control load aggregator k in step S12 can be equivalently transformed into:

[0181]

[0182] S32. Transform the inequality constraints into equality constraints

[0183] First, denote Denote

[0184] , introduce slack variables Then, the model of the building temperature control load aggregator k in S31 can be equivalently transformed into:

[0185]

[0186] Among them, is a random diagonal matrix with all diagonal elements being positive and is E k is owned by the building temperature control load aggregator k.

[0187] S33. Implement constraint equivalent expansion

[0188] First, copy the model of the building temperature control load aggregator k in step S32 as follows:

[0189]

[0190] Denote Then, the above model of the building temperature control load aggregator k can be equivalently transformed into:

[0191]

[0192] S34. Implement an encryption method based on random transformation II

[0193] First, the building temperature control load aggregator k generates a reversible random matrix Then, apply this random matrix to the model of the building temperature control load aggregator k formed in S33 to obtain:

[0194]

[0195] S35. Implement economic dispatch considering building temperature control load aggregators and their privacy and security

[0196] The building temperature control load aggregator k uploads the matrix V K F k , V k G k , V k H k and V k e k to the power system operator, so that the power system operator can formulate the economic dispatch problem considering the building temperature control load aggregator and its privacy and security as follows:

[0197]

[0198] where U k The expression is as follows:

[0199]

[0200] The power system operator solves the above economic dispatch problem and returns the solution result to the building temperature control load aggregator k, and the building temperature control load aggregator k can restore the true state variable x according to the following formula k :

[0201] Example 2

[0202] The building temperature control load aggregator - building cluster system in this example consists of an IEEE33 - node distribution network system and three building temperature control load aggregators accessing the system. The schematic diagram of the system structure is as Figure 2 shown. The random matrix follows a normal distribution with a mean of 0.1 and a standard deviation of 0.1, and the order M of the thermal dynamic aggregation model is set to 1. The simulation is executed on a personal computer with an Intel i7 core and 32GB of RAM. The program is based on MATLAB R2022b and Yalmip software, and Gurobi 10.0.2 is used to solve the economic dispatch optimization problem. The parameter settings of the building temperature control load aggregator model are shown in the following table:

[0203] Table 1: Parameter table of the building temperature control load aggregator model

[0204]

[0205] Solve the economic dispatch problem considering the building temperature control load aggregator and its privacy and security according to the steps of the present invention. Figure 3 It is the comparison chart of the cooling / heating power of the building temperature control load aggregator 1 under the privacy protection and non - privacy protection methods in Example 2. It can be seen that the economic dispatch result of the privacy - protected model calculation method proposed by the present invention has excellent accuracy. Figure 4 (a) and (b) are the matrix G in Example 21 With V 1 G 1 The comparison result of the heat map of the first 24 lines shows that the real information G 1 is covered, proving that the method proposed in this patent has good privacy protection performance.

[0206] In the description of this specification, the description with reference to terms such as "one embodiment", "example", "specific example", etc. means that the specific features, structures, materials or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described can be combined in any one or more embodiments or examples in a suitable manner.

[0207] The above shows and describes the basic principles, main features and advantages of the present invention. Those skilled in the art should understand that the present invention is not limited by the above embodiments, and what is described in the above embodiments and the specification is only to illustrate the principle of the present invention. Without departing from the spirit and scope of the present invention, the present invention will have various changes and improvements, and these changes and improvements all fall within the scope of the present invention claimed.

Claims

1. A centralized scheduling method considering building temperature control load aggregators and their privacy security, characterized in that: include: S1. Establish a building temperature control load aggregator model; S2. Establish a distribution network economic dispatch model considering building temperature control load aggregators; S3. Establish a scheduling model calculation method that takes privacy security into consideration, specifically including: S31. Implement an encryption method based on random transformation I; S32. Convert inequality constraints into equality constraints; S33. Implement constraint equivalent expansion; S34. Implement an encryption method based on random transformation II; S35. Implement economic scheduling that takes into consideration building temperature control load aggregators and their privacy security.

2. A centralized scheduling method considering building temperature control load aggregators and their privacy security according to claim 1, characterized in that: The step S1 of establishing the building temperature control load aggregator model further includes: S11. Establish an algebraic building temperature control load aggregation model: The building temperature control load aggregator interacts with the distribution network as an agent of the building user, and its model is established as follows: in, represents the aggregate temperature of the building temperature control load aggregator k, represents the cooling / heating power of the building temperature control load aggregate k, is the model parameter of the building temperature control load aggregator k, and the set M = {0, 1, ..., M} represents the model order. and τ k are the upper and lower limits of the aggregation temperature of the building temperature control load aggregator k respectively; S12. Establish a compact building temperature control load aggregation model: First, define the constant matrix Where m is the model order, the element in the i-th row and j-th column of the matrix (Λ m ) i,j for: Second, use the state variable To represent the aggregate temperature of the building temperature control load aggregator k Using control variables To represent the cooling / heating power of the building temperature control load aggregate k Define the state variable vector x separately k and the control variable vector u k for: Again, define the constant vector The tth element With the following expression: Defining the Matrix Among them, the constant matrix I T represents a T-dimensional unit matrix. Therefore, the building temperature control load aggregation model established in step S11 is converted into its compact form: R k x k +S k u k =d k , Among them, 1 T Represents a T-dimensional 1-vector.

3. According to claim 1, a centralized scheduling method considering building temperature control load aggregators and their privacy security is characterized in that: The step S2 of establishing a distribution network economic dispatch model considering building temperature control load aggregators further includes: S21, establishing distribution network constraints; S22, establishing distribution network and building temperature control load aggregator coupling constraints; S23, establishing building temperature control load aggregator control variable constraints; S24, establishing a distribution network economic dispatch model considering building temperature control load aggregators.

4. A centralized scheduling method considering building temperature control load aggregators and privacy security according to claim 3, characterized in that: The step S21, establishing the distribution network constraint conditions specifically includes: The economic dispatch decision variable of the distribution network is vector z, which does not include the cooling / heating power of each building temperature control load aggregator Let the feasible region of z be Z, and then we get the distribution network constraints: z∈Z S22. Establishing the coupling constraints between the distribution network and the building temperature control load aggregator specifically includes: Record the control variables of each building temperature control load aggregator The vector composed The coupling constraint condition between the distribution network and the building temperature control load aggregator is expressed as: Az+u=0, The matrix A represents the connection relationship between the distribution network bus and the building temperature control load aggregator; S23. Establishing the control variable constraints of the building temperature control load aggregator specifically includes: Define the building temperature control load aggregator control variable u k , according to step S22, u k The feasible domain U k for: S24. Establishing a distribution network economic dispatch model considering building temperature control load aggregators specifically includes: The distribution network economic dispatch model considering the building temperature control load aggregator is expressed as: stz∈Z,Az+u=0, in k ∈U k ,k∈K Where c represents the cost vector, U k The expression is shown in S23, that is:

5. A centralized scheduling method considering building temperature control load aggregators and privacy security according to claim 1, characterized in that: The calculation method of establishing the scheduling model considering privacy security in step S3 is: S31. Implementing an encryption method based on random transformation I First, each building temperature control load aggregator k generates a reversible random matrix Then each building temperature control load aggregator k model is equivalently transformed into: S32. Convert inequality constraints into equality constraints First, remember remember Introducing slack variables Therefore, the building temperature control load aggregator k model in step S31 is equivalently transformed into: in, is a random diagonal matrix with positive diagonal elements and E k It is owned by building temperature control load aggregator k; S33. Implement constraint equivalence extension First, the model of the building temperature control load aggregator k in step S32 is copied as follows: remember Therefore, the model of the above building temperature control load aggregator k is equivalently transformed into: S34. Implementing encryption method based on random transformation II First, the building temperature control load aggregator k generates a reversible random matrix Then, the random matrix is ​​applied to the model of the building temperature control load aggregator k formed in step S33 to obtain: S35. Implement economic dispatch considering building temperature control load aggregators and their privacy security The building temperature control load aggregator k converts the matrix V K F k , V k G k , V k H k and V k e k Uploaded to the power system operator, the power system operator then poses the following economic dispatch problem considering the building temperature control load aggregator and its privacy security: stz∈Z,Az+u=0, in k ∈U k ,k∈K Among them, U k The expression is as follows: The power system operator solves the above economic dispatch problem and sends the solution results Returned to the building temperature control load aggregator k, the building temperature control load aggregator k can restore the real state variable x according to the following formula k :

6. An economic dispatch system considering building temperature control load aggregators and their privacy security, characterized in that: include: Building control load aggregator establishment module, used to establish building temperature control load aggregator model; The distribution network economic dispatch model establishment module is used to establish a distribution network economic dispatch model that takes into account the building temperature control load aggregator; The privacy protection model is used to establish a scheduling model calculation method that takes privacy security into consideration, specifically including: a first encryption unit, used to implement an encryption method based on random transformation I; a constraint conversion unit, used to convert inequality constraints into equality constraints; a constraint expansion unit, used to implement constraint equivalent expansion; a second encryption unit, used to implement an encryption method based on random transformation II; an implementation unit, used to implement economic scheduling that takes into consideration building temperature control load aggregators and their privacy security.

7. The economic dispatch system considering building temperature control load aggregators and their privacy security according to claim 6 is characterized in that: The building temperature control load aggregator model establishment module includes the following units: The algebraic building temperature control load aggregation model unit is used to establish the algebraic building temperature control load aggregation model: The building temperature control load aggregator interacts with the distribution network as an agent of the building user, and its model is established as follows: in, represents the aggregate temperature of the building temperature control load aggregator k, represents the cooling / heating power of the building temperature control load aggregate k, is the model parameter of the building temperature control load aggregator k, and the set M = {0, 1, ..., M} represents the model order. and τ k are the upper and lower limits of the aggregation temperature of the building temperature control load aggregator k respectively; A compact building temperature control load aggregation model building unit is used to build a compact building temperature control load aggregation model: First, define the constant matrix Where m is the model order, the element in the i-th row and j-th column of the matrix (Λ m ) i,j for: Second, use the state variable To represent the aggregate temperature of the building temperature control load aggregator k Using control variables To represent the cooling / heating power of the building temperature control load aggregate k Define the state variable vector x separately k and the control variable vector u k for: Again, define the constant vector The tth element With the following expression: Defining the Matrix Among them, the constant matrix I T represents the T-dimensional unit matrix, so the established building temperature control load aggregation model is transformed into its compact form: R k x k +S k u k =d k , Among them, 1 T Represents a T-dimensional 1-vector.

8. The economic dispatch system considering building temperature control load aggregators and their privacy security according to claim 6 is characterized in that: The distribution network economic dispatch model establishment module further includes: a distribution network constraint condition establishment unit, which is used to establish distribution network constraints; a distribution network and building temperature control load aggregator coupling constraint condition establishment unit, which is used to establish distribution network and building temperature control load aggregator coupling constraint conditions; a building temperature control load aggregator control variable constraint condition establishment unit, which is used to establish building temperature control load aggregator control variable constraint conditions; a distribution network economic dispatch model establishment unit, which is used to establish a distribution network economic dispatch model considering the building temperature control load aggregator.

9. The economic dispatch system considering building temperature control load aggregators and their privacy security according to claim 8, characterized in that: The distribution network constraint condition establishment unit specifically performs the following steps: The economic dispatch decision variable of the distribution network is vector z, which does not include the cooling / heating power of each building temperature control load aggregator Let the feasible region of z be Z, and then we get the distribution network constraints: z∈Z The unit for establishing coupling constraints between the distribution network and the building temperature control load aggregator specifically performs the following steps: Record the control variables of each building temperature control load aggregator The vector composed The coupling constraint condition between the distribution network and the building temperature control load aggregator is expressed as: Az+u=0, The matrix A represents the connection relationship between the distribution network bus and the building temperature control load aggregator; The building temperature control load aggregator control variable constraint condition establishment unit specifically performs the following steps: Define the building temperature control load aggregator control variable u k , according to step S22, u k The feasible domain U k for: The distribution network economic dispatch model establishment unit specifically performs the following steps: The distribution network economic dispatch model considering the building temperature control load aggregator is expressed as: stz∈Z,Az+u=0, in k ∈U k ,k∈K Where c represents the cost vector, U k The expression is shown in S23, that is:

10. The economic dispatch system considering building temperature control load aggregators and their privacy security according to claim 6, characterized in that: The privacy protection model includes: The first encryption unit performs the following steps: First, each building temperature control load aggregator k generates a reversible random matrix Then each building temperature control load aggregator k model is equivalently transformed into: To constrain the transformation unit, perform the following steps: First, remember remember Introducing slack variables Therefore, the building temperature control load aggregator k model in step S31 is equivalently transformed into: in, is a random diagonal matrix with positive diagonal elements and E k It is owned by building temperature control load aggregator k; To constrain the expansion unit, perform the following steps: First, the model of the building temperature control load aggregator k in step S32 is copied as follows: remember Therefore, the model of the above building temperature control load aggregator k is equivalently transformed into: The second encryption unit performs the following steps: First, the building temperature control load aggregator k generates a reversible random matrix Then, the random matrix is ​​applied to the model of the building temperature control load aggregator k formed in step S33 to obtain: The implementation unit performs the following steps: The building temperature control load aggregator k converts the matrix V K F k , V k G k , V k H k and V k e k Uploaded to the power system operator, the power system operator then poses the following economic dispatch problem considering the building temperature control load aggregator and its privacy security: stz∈Z,Az+u=0, in k ∈U k ,k∈K Among them, U k The expression is as follows: The power system operator solves the above economic dispatch problem and sends the solution results Returned to the building temperature control load aggregator k, the building temperature control load aggregator k can restore the real state variable x according to the following formula k :

Images