Wide-area electric energy quality disturbance problem global cooperative treatment method, system and equipment

Abstract

The invention relates to the technical field of electric energy quality management research, and discloses a global collaborative management method, system and equipment for a wide-area electric energy quality disturbance problem, and the method comprises the steps: constructing a node weight coefficient matrix considering electric energy quality transmission characteristics, responsibility contribution degree, electric energy quality sensitive load and management cost indexes based on all nodes in a power grid; calculating a weighted voltage distortion rate of each node based on the node weight coefficient matrix; constructing a governance scheme optimization model which takes the governance node selection and the adopted governance scheme as optimization objects and takes the minimum sum of weighted voltage distortion rates of all nodes in the power grid as a target function; and solving the governance scheme optimization model by adopting a golden litsea optimization algorithm to obtain an optimal governance node and a governance scheme corresponding to the optimal governance node. According to the invention, the optimization model is constructed by comprehensively considering various treatment schemes, and the proposed technology can integrate the advantages of various treatment schemes to obtain a collaborative treatment method of various treatment schemes.

Description

Technical Field

[0001] The present invention relates to the technical field of power quality management research, and in particular to a global collaborative management method, system and equipment for wide-area power quality disturbance problems. Background Art

[0002] The "double-high" new power system is a typical wide-area large power grid with high voltage levels, wide distribution areas and various types of harmonic emission sources. Although the harmonics of single nonlinear users, new energy sites and converter stations entering the network are suppressed to meet the standards, due to the strong background harmonic fluctuations and the harmonic transmission amplification effect of cables, it is obviously no longer suitable for the harmonic problems of the "double-high" new power system to only suppress harmonics in a single nonlinear user, new energy site or converter station without considering wide-area coordination. Because the harmonics emitted by the harmonic source are coupled with the background harmonics, when it is detected that the harmonics of a certain site exceed the standard, it is considered that the nonlinear users, new energy sites or converter stations connected to the site are the dominant harmonic sources, and it may be a serious mistake to manage them accordingly.

[0003] Since harmonic sources are widely distributed in large power grids and resonance amplification may occur during the transmission process, node harmonic exceeding the standard may be the result of the combined effect of nonlinear users / new energy sites and background harmonics. Governance at only a single node cannot effectively control the wide-area harmonic problem of the power grid. At the same time, when the harmonics in the large power grid exceed the standard in a wide area and the background harmonics are the dominant harmonic source, there is no effective analysis method to select which nodes to govern. In addition, there are various types of node harmonic governance schemes. Only by comprehensively considering various governance schemes and building a wide-area harmonic governance node optimization model can the harmonic problem of the large power grid be fundamentally governed. Summary of the invention

[0004] The present invention provides a global collaborative management method, system and equipment for wide-area power quality disturbance problems to solve the above problems.

[0005] The present invention is achieved through the following technical solutions: A global collaborative governance method for wide-area power quality disturbance problems, comprising: Based on all nodes in the power grid, a node weight coefficient matrix is ​​constructed that takes into account power quality transfer characteristics, responsibility contribution, power quality sensitive loads and governance cost indicators; Calculating the weighted voltage distortion rate of each node based on the node weight coefficient matrix; Construct a governance scheme optimization model with the governance node selection and the adopted governance scheme as the optimization object and the minimum sum of weighted voltage distortion rates of all nodes in the power grid as the objective function; The golden jackal optimization algorithm is used to solve the governance scheme optimization model to obtain the optimal governance node and the governance scheme corresponding to the optimal governance node.

[0006] As an optimization, based on all nodes in the power grid, the specific process of constructing a node weight coefficient matrix that considers power quality transfer characteristics, responsibility contribution, power quality sensitive loads and governance costs is as follows: A1. Calculate the sub-responsibility contribution of all harmonic source nodes in the power grid to each harmonic exceeding node in the power grid: ; In the formula, It represents the sub-responsibility contribution of the harmonic source node m in the power grid to the harmonic exceeding node x, Indicates the harmonic voltage at the node x where the harmonic exceeds the standard, represents the harmonic voltage of the harmonic source node m on the bus where the harmonic exceeding node x is located, express and The angle difference between ; M represents the number of harmonic sources that have been monitored, ,and , represents the measured harmonic injection current of the harmonic source node m, Indicates the bus where the harmonic exceeding node x is located and the harmonic source node m The harmonic transfer impedance between Indicates the background harmonic voltage of the bus where the harmonic exceeding node x is located; A2. Sum up the sub-responsibility contributions of the same harmonic exceeding node m to obtain the responsibility contribution index of the harmonic exceeding node m , ,in, is the total number of nodes exceeding the standard; A3. Set the non-harmonic source nodes Responsibility contribution index Recorded as 0, the node contribution index matrix HI of all nodes in the power grid is finally constructed; A4. Construct a transfer characteristic matrix of all nodes in the power grid regarding the power quality transfer characteristics, thereby obtaining the participation factors of all nodes, wherein the transfer characteristic matrix is ​​composed of the harmonic current of each node divided by the admittance matrix corresponding to each node, and the transfer characteristic matrix is ​​expressed as: ; in, represents the voltage vector of each node, Representation Node j The voltage, , Y is the admittance matrix composed of each node, , represents the eigenvalue diagonal matrix, , They are respectively the admittance matrices Y The left and right eigenvector matrices of is the current matrix composed of the harmonic currents of each node, n is the total number of nodes in the power grid, and the diagonal elements in the admittance matrix are defined as For Node j The participation factor , thus obtaining the participation factor matrix composed of the participation factors of all nodes in the power grid PF , , the participation factor is expressed as a combination of the resonant excitability and observability of mode s on the same bus; A5. Consider power quality sensitive loads to set nodes j The sensitivity factor , let the governance cost index of the node be j The proportion of the governance cost of the whole network to the total governance cost of the whole network is recorded as , thus obtaining the sensitivity factor matrix consisting of the sensitivity factors and governance cost indicators of all nodes in the power grid: and governance cost indicator matrix Cost ; A6. Based on the node contribution index matrix HI , participation factor matrix PF , sensitivity factor matrix and governance cost indicator matrix Cost Construct a node weight coefficient matrix for all nodes in the power grid Q ,in, Q It is expressed as: ; in, represents the i-th node weight coefficient of node j, .

[0007] As an optimization, the specific process of calculating the weighted voltage distortion rate of each node based on the node weight coefficient matrix is: As an optimization, the specific process of calculating the weighted voltage distortion rate of each node based on the node weight coefficient matrix is: B1. Calculate the weight of each node weight coefficient in the node weight coefficient matrix based on the entropy weight method, where i The weight of the node weight coefficient is expressed as: ; in, Indicates i The entropy of the node weight coefficients, ; B2. Calculate the weight of each node based on the weight of each node weight coefficient, and the weight of the node is expressed as: ; in, represents the weight of node j; B3. Calculate the harmonic voltage of each node: ; Where h represents the harmonic order, represents the effective value of the harmonic voltage of the hth harmonic at the node, represents the admittance matrix of the hth harmonic, It represents the effective value of the harmonic injection current of the hth harmonic passing through the node; B4. Calculate the voltage distortion rate of each node: ; in, Representation Node j The voltage distortion rate, Representation Node j The effective value of the fundamental voltage, For Node j The hth harmonic voltage effective value; B5. Based on the weight of each node and the voltage distortion rate, a weighted voltage distortion rate of each node is constructed. The weighted voltage distortion rate is expressed as: ; in, Representation Node j The weighted voltage distortion rate.

[0008] As an optimization, the governance scheme optimization model is expressed as: ; in, represents the preset total governance investment, Indicates whether the node j Adopt governance solutions, Indicates that at the node j Adopt governance solutions, Indicates that at the node j No governance solution is adopted. For Node j The investment required for governance solutions.

[0009] As an optimization, the golden jackal optimization algorithm is used to solve the governance scheme optimization model, and the specific process of obtaining the optimal governance node and the corresponding governance scheme is as follows: C1. Initialize the population position and the maximum number of iterations, wherein the population includes a plurality of prey, each prey consisting of whether a node in the power grid adopts a governance scheme and a corresponding governance scheme; C2. Calculate the fitness value of each individual's location, select the prey with the best fitness value as the male golden jackal, and the prey with the second best fitness value as the female golden jackal; C3. In the exploration phase, the positions of the male and female golden jackals are updated. The specific formula is: ; ; , represent the positions of male and female golden jackals at the tth iteration, , They represent the updated positions of the male and female golden jackals in the tth iteration, E is the escape energy of the prey, is a random number between 0 and 1. represents a random number based on Levy distribution, Represents the prey matrix of the tth iteration ; C4. When the escape energy E of the prey is less than 1, the development phase is entered, and the positions of the male and female golden jackals are updated. The specific formula is: ; ; C5. When the maximum number of iterations is reached When , the position of the golden jackal is updated and the updated position of the golden jackal is used as the output, where the update formula of the position of the golden jackal is: ; in, yes The position of the male golden jackal in the iteration, yes The position of the female golden jackal in iteration .

[0010] As an optimization, the prey matrix It is expressed as: ; in, represents the kth prey, , represents the solution to the e-th dimension optimization problem for the k-th prey, , , d is the dimension of the optimization problem to be solved, and p is the number of prey in the population.

[0011] As an optimization, the prey matrix It is expressed as: ; in, represents the kth prey, , represents the solution to the e-th dimension optimization problem for the k-th prey, , , d is the dimension of the optimization problem to be solved, and p is the number of prey in the population.

[0012] As an optimization, random numbers based on Levy distribution The specific expression is: ; in, is the Levy flight function, , is the intermediate parameter, , and is a random number ranging from 0 to 1. is a constant, represents the Gamma function, .

[0013] The present invention also discloses a global collaborative governance system for wide-area power quality disturbance problems, which is used to execute the aforementioned global collaborative governance method for wide-area power quality disturbance problems, including: The node weight coefficient matrix construction module is used to construct a node weight coefficient matrix based on all nodes in the power grid, taking into account power quality transfer characteristics, responsibility contribution, power quality sensitive loads and governance cost indicators; A weighted voltage distortion rate calculation module, used to calculate the weighted voltage distortion rate of each node based on the node weight coefficient matrix; A governance scheme optimization model building module is used to build a governance scheme optimization model with the governance node selection and the adopted governance scheme as the optimization objects and the minimum sum of weighted voltage distortion rates of all nodes in the power grid as the objective function; The solution module is used to solve the governance scheme optimization model by using the golden jackal optimization algorithm to obtain the best governance node and the governance scheme corresponding to the best governance node.

[0014] The present invention also discloses a data processing device, comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the aforementioned global collaborative governance method for wide-area power quality disturbance problems when executing the computer program.

[0015] Compared with the prior art, the present invention has the following advantages and beneficial effects: The present invention comprehensively considers various governance schemes to construct an optimization model. Compared with the traditional optimization model that only targets a single governance device, the proposed technology can integrate the advantages of various governance schemes. In addition, the existing method does not consider the impact of the access of the governance equipment on the network admittance. The present invention will continuously update the harmonic injection current and admittance matrix after governance during the optimization process. In order to achieve the comprehensive optimal harmonic voltage distortion of the power grid, the present invention comprehensively considers harmonic source governance and harmonic network transmission blocking, and constructs a node weight coefficient that considers harmonic transmission characteristics, responsibility contribution, power quality sensitive load and governance cost. The sum of the weighted harmonic voltage distortion rates of each node in the network is taken as the optimization target, and the total investment in harmonic governance is taken as the constraint condition. The optimal harmonic suppression strategy of the power grid is solved, and finally a collaborative governance method that comprehensively considers various governance schemes is obtained. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] The drawings described herein are used to provide a further understanding of the embodiments of the present invention, constitute a part of this application, and do not constitute a limitation of the embodiments of the present invention. In the drawings: Figure 1 It is a schematic diagram of the process of the present invention; Figure 2 To evaluate the schematic diagram of the power grid topology; Figure 3 The harmonic source current diagram for the setting; Figure 4 Control the harmonic distortion rate before and after each node; Figure 5 Optimize the model optimization process for harmonic mitigation solutions. DETAILED DESCRIPTION

[0017] In order to make the objectives, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with embodiments and drawings. The exemplary embodiments of the present invention and their description are only used to explain the present invention and are not intended to limit the present invention.

[0018] This embodiment 1 provides a global collaborative governance method for wide-area power quality disturbance problems, such as Figure 1 As shown, including: S1. Based on all nodes in the power grid, a node weight coefficient matrix is ​​constructed that takes into account power quality transfer characteristics, responsibility contribution, power quality sensitive loads and governance cost indicators; In some embodiments, the specific process of step S1 is: A1. Calculate the sub-responsibility contribution of all harmonic source nodes in the power grid to each harmonic exceeding node in the power grid: ; In the formula, Indicates the harmonic source node in the power grid mSub-responsibility contribution to the harmonic exceeding node x, i.e. the harmonic source node m Attention Node x The projection percentage of the generated harmonic voltage on the harmonic voltage of the node of interest, Indicates the harmonic voltage at the node x where the harmonic exceeds the standard, represents the harmonic voltage of the harmonic source node m on the bus where the harmonic exceeding node x is located, express and The angle difference between ; M represents the number of harmonic sources that have been monitored, ,and , represents the measured harmonic injection current of the harmonic source node m, Indicates the bus where the harmonic exceeding node x is located and the harmonic source node m The harmonic transfer impedance between (elements of the harmonic transfer impedance matrix), Indicates the background harmonic voltage of the bus where the harmonic exceeding node x is located; It should be noted that voltage can be written as a phasor of amplitude plus angle, and angle difference refers to the subtraction of two angles.

[0019] A2. Sum up the sub-responsibility contributions of the same harmonic exceeding node m to obtain the responsibility contribution index of the harmonic exceeding node m , ,in, is the total number of nodes exceeding the standard; A3. Set the non-harmonic source nodes Responsibility contribution index Recorded as 0, the node contribution index matrix HI of all nodes in the power grid is finally constructed; here, non-harmonic source nodes The harmonic source node m and the harmonic source node m together constitute all the nodes in the power grid.

[0020] A4. Construct a transfer characteristic matrix of all nodes in the power grid regarding the power quality transfer characteristics, thereby obtaining the participation factors of all nodes, wherein the transfer characteristic matrix is ​​composed of the harmonic current of each node divided by the admittance matrix corresponding to each node, and the transfer characteristic matrix is ​​expressed as: ; in, represents the voltage vector of each node, Representation Node j The voltage, , Y is the admittance matrix composed of each node, , represents the eigenvalue diagonal matrix, , They are respectively the admittance matrices Y The left and right eigenvector matrices of is the current matrix composed of the harmonic currents of each node, n is the total number of nodes in the power grid, and the diagonal elements in the admittance matrix are defined as For Node j The participation factor , thus obtaining the participation factor matrix composed of the participation factors of all nodes in the power grid PF , , the participation factor is expressed as a combination of the resonant excitability and observability of mode s on the same bus; A5. Consider power quality sensitive loads to set nodes j The sensitivity factor , let the governance cost index of the node be j The proportion of the governance cost of the whole network to the total governance cost of the whole network is recorded as , thus obtaining the sensitivity factor matrix consisting of the sensitivity factors and governance cost indicators of all nodes in the power grid: and governance cost indicator matrix Cost ; A6. Based on the node contribution index matrix HI , participation factor matrix PF , sensitivity factor matrix and governance cost indicator matrix Cost Construct a node weight coefficient matrix for all nodes in the power grid Q ,in, Q It is expressed as: ; in, represents the i-th node weight coefficient of node j, .

[0021] S2. Calculating the weighted voltage distortion rate of each node based on the node weight coefficient matrix; It should be noted that when solving the optimization model of the governance scheme in the subsequent Golden Jackal algorithm, the weighted voltage distortion rate must be solved. Therefore, the harmonic current injected by the disturbance source after governance is updated according to the governance scheme. With the admittance matrix , including the filter impedance and the harmonic current source equivalent to the active filter. Calculate the harmonic flow after treatment and calculate the weighted voltage distortion rate of each node.

[0022] The specific process of step S2 is: B1. Calculate the weight of each node weight coefficient in the node weight coefficient matrix based on the entropy weight method. A node has a total of 4 node weight coefficients, namely node contribution index, participation factor, sensitivity factor and governance cost index. This step is to calculate the weights of node contribution index, participation factor, sensitivity factor and governance cost index respectively.

[0023] Among them, i The weight of the node weight coefficient is expressed as: ; in, Indicates i The entropy of the node weight coefficients, ; i=1 means that the node weight coefficient is the node contribution index, i=2 means that the node weight coefficient is the participation factor, i=3 means that the node weight coefficient is the sensitivity factor, and i=4 means that the node weight coefficient is the governance cost index.

[0024] B2. Calculate the weight of each node based on the weight of each node weight coefficient, and the weight of the node is expressed as: ; in, represents the weight of node j; B3. Calculate the harmonic voltage of each node: ; Where h represents the harmonic order, represents the effective value of the harmonic voltage of the hth harmonic at the node, represents the admittance matrix of the hth harmonic, Represents the effective value of the harmonic current of the hth harmonic passing through the node; B4. Calculate the voltage distortion rate of each node: ; in, Representation Node j The voltage distortion rate, Representation Node j The effective value of the fundamental voltage, For Node j The hth harmonic voltage effective value; It should be noted that, generally speaking, when performing calculations, h only needs to be considered up to 20, that is, the harmonic frequency order is 19.

[0025] B5. Based on the weight of each node and the voltage distortion rate, a weighted voltage distortion rate of each node is constructed. The weighted voltage distortion rate is expressed as: ; in, Representation Node j The weighted voltage distortion rate.

[0026] S3. Construct a governance scheme optimization model with the governance node selection and the adopted governance scheme as the optimization object and the minimum sum of weighted voltage distortion rates of all nodes in the power grid as the objective function; The governance scheme optimization model is expressed as: ; in, represents the preset total governance investment, Indicates whether the node j Adopt governance solutions, Indicates that at the node j Adopt governance solutions, Indicates that at the node j No governance solution is adopted. For Node j The investment required for governance solutions.

[0027] S4. Use the golden jackal optimization algorithm to solve the governance scheme optimization model to obtain the best governance node and the governance scheme corresponding to the best governance node.

[0028] The specific process of step S4 is: C1. Initialize the population position and the maximum number of iterations, wherein the population (golden jackal population) includes multiple prey, each prey consisting of whether a node in the power grid adopts a governance scheme and the corresponding governance scheme; More specifically, the governance solution is to govern the capacity of the device.

[0029] It should be noted that the specific formula for initializing the population is: ; in, is the initial position of the golden jackal population. , are the upper and lower bounds of the problem to be solved, Indicates taking and A random number between .

[0030] Here, the problem to be solved is the problem that needs to be optimized, that is, the type of governance solution adopted by the node (which can be understood as whether the node adopts the governance solution) and the capacity of the governance device.

[0031] C2. Calculate the fitness value of each individual's location, select the prey with the best fitness value as the male golden jackal, and the prey with the second best fitness value as the female golden jackal; Here, the fitness function may be an objective function. The optimal fitness value refers to the minimum fitness value, and the suboptimal fitness value refers to the second smallest fitness value.

[0032] C3. In the exploration phase, i.e. searching for prey, the male golden jackal is updated and female golden jackals The specific formula is: ; ; , represent the positions of male and female golden jackals at the tth iteration, , They represent the updated positions of the male and female golden jackals in the tth iteration, E is the escape energy of the prey, is a random number between 0 and 1. represents a random number based on Levy distribution, Represents the prey matrix of the tth iteration ; C4. When the escape energy E of the prey is less than 1, the development phase is entered, that is, the phase of surrounding and attacking the prey, and the positions of the male and female golden jackals are updated. The specific formula is: ; ; C5. When the maximum number of iterations is reached When , the position of the golden jackal is updated and the updated position of the golden jackal is used as the output, where the update formula of the position of the golden jackal is: ; in, yes The position of the male golden jackal in the iteration, yes The position of the female golden jackal in iteration .

[0033] In some embodiments, the prey matrix It is expressed as: ; in, represents the kth prey, , represents the solution to the e-th dimension optimization problem for the k-th prey, , , d is the dimension of the optimization problem to be solved, and p is the number of prey in the population, that is, the size of the prey population.

[0034] Here, the dimension is the number of nodes multiplied by the harmonic frequency multiplied by 2. That is, the type of treatment plan and the capacity of the treatment device must be considered for each harmonic and each node.

[0035] In some embodiments, the escape energy E of the prey is expressed as: ; in, represents the decreasing process of prey energy, represents the initial state of prey energy, is a random number in the range of 0-1, is a constant, , t is the current iteration number.

[0036] In some embodiments, a random number based on a Levy distribution The specific expression is: ; Among them, among them, is the Levy flight function, , , and is a random number ranging from 0 to 1. is a constant, represents the Gamma function, .

[0037] Next, the present invention will be described using specific examples.

[0038] The present invention uses IEEE33 nodes for simulation calculation, and the power grid topology is as follows: Figure 2 As shown, 10 harmonic sources are set up, and the harmonic current injected by each harmonic source is randomly 1%~10% of the fundamental current. The node harmonic current is set as follows Figure 3 As shown, the treatment method adopts single tuned filter and active filter, and the cost of treatment device is shown in Table 1-Table 2.

[0039] Table 1 Cost of single tuned filter Single tuned filter Single tuned filter R 396.5 yuan / MW L 55.25 yuan / MVar C 32500 yuan / MVar Table 2 Cost of active filter Active filter 3000 yuan / A The total investment in network governance is set at 50,000 yuan. With the total investment as the constraint condition and the sum of weighted voltage distortion rate as the objective function, a harmonic governance optimization model is established and solved using the golden jackal optimization algorithm. The harmonic distortion rate of each node before and after governance is shown in the figure below: Figure 4 As shown in the figure, it can be seen that the distortion rate of each node is reduced to below 0.2% after treatment. The optimization process of the harmonic control scheme optimization model solved by the golden jackal optimization algorithm is as follows: Figure 5 shown.

[0040] Embodiment 2 further discloses a global collaborative governance system for wide-area power quality disturbance problems, including: The node weight coefficient matrix construction module is used to construct a node weight coefficient matrix based on all nodes in the power grid, taking into account power quality transfer characteristics, responsibility contribution, power quality sensitive loads and governance cost indicators; A weighted voltage distortion rate calculation module, used to calculate the weighted voltage distortion rate of each node based on the node weight coefficient matrix; A governance scheme optimization model building module is used to build a governance scheme optimization model with the governance node selection and the adopted governance scheme as the optimization objects and the minimum sum of weighted voltage distortion rates of all nodes in the power grid as the objective function; The solution module is used to solve the governance scheme optimization model by using the golden jackal optimization algorithm to obtain the best governance node and the governance scheme corresponding to the best governance node.

[0041] Embodiment 3 also discloses a data processing device, characterized in that it comprises: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, it implements the global collaborative governance method for wide-area power quality disturbance problems described in Embodiment 1.

[0042] The specific implementation methods described above further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above description is only a specific implementation method of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A global collaborative governance method for wide-area power quality disturbance problems, characterized in that: include: Based on all nodes in the power grid, a node weight coefficient matrix is ​​constructed that takes into account power quality transfer characteristics, responsibility contribution, power quality sensitive loads and governance cost indicators; Calculating the weighted voltage distortion rate of each node based on the node weight coefficient matrix; Construct a governance scheme optimization model with the governance node selection and the adopted governance scheme as the optimization object and the minimum sum of weighted voltage distortion rates of all nodes in the power grid as the objective function; The golden jackal optimization algorithm is used to solve the governance scheme optimization model to obtain the optimal governance node and the governance scheme corresponding to the optimal governance node.

2. The global collaborative governance method for wide-area power quality disturbance problems according to claim 1 is characterized in that: Based on all nodes in the power grid, the specific process of constructing a node weight coefficient matrix that considers power quality transfer characteristics, responsibility contribution, power quality sensitive loads and governance costs is as follows: A1. Calculate the sub-responsibility contribution of all harmonic source nodes in the power grid to each harmonic exceeding node in the power grid: ； In the formula, It represents the sub-responsibility contribution of the harmonic source node m in the power grid to the harmonic exceeding node x, Indicates the harmonic voltage at the node x where the harmonic exceeds the standard, represents the harmonic voltage of the harmonic source node m on the bus where the harmonic exceeding node x is located, express and The angle difference between ; M represents the number of harmonic sources that have been monitored, ,and , represents the measured harmonic injection current of the harmonic source node m, Indicates the bus where the harmonic exceeding node x is located and the harmonic source node m The harmonic transfer impedance between Indicates the background harmonic voltage of the bus where the harmonic exceeding node x is located; A2. Sum up the sub-responsibility contributions of the same harmonic exceeding node m to obtain the responsibility contribution index of the harmonic exceeding node m , ,in, is the total number of nodes exceeding the standard; A3. Set the non-harmonic source nodes Responsibility contribution index Recorded as 0, the node contribution index matrix HI of all nodes in the power grid is finally constructed; A4. Construct a transfer characteristic matrix of all nodes in the power grid regarding the power quality transfer characteristics, thereby obtaining the participation factors of all nodes, wherein the transfer characteristic matrix is ​​composed of the harmonic current of each node divided by the admittance matrix corresponding to each node, and the transfer characteristic matrix is ​​expressed as: ； in, represents the voltage vector of each node, Representation Node j The voltage, , Y is the admittance matrix composed of each node, , represents the eigenvalue diagonal matrix, , They are respectively the admittance matrices Y The left and right eigenvector matrices of is the current matrix composed of the harmonic currents of each node, n is the total number of nodes in the power grid, and the diagonal elements in the admittance matrix are defined as For Node j The participation factor , thus obtaining the participation factor matrix composed of the participation factors of all nodes in the power grid PF , , the participation factor is expressed as a combination of the resonant excitability and observability of mode s on the same bus; A5. Consider power quality sensitive loads to set nodes j The sensitivity factor , let the governance cost index of the node be j The proportion of the governance cost of the entire network to the total governance cost is recorded as , thus obtaining the sensitivity factor matrix consisting of the sensitivity factors and governance cost indicators of all nodes in the power grid: and governance cost indicator matrix Cost ; A6. Based on the node contribution index matrix HI , participation factor matrix PF , sensitivity factor matrix and governance cost indicator matrix Cost Construct a node weight coefficient matrix for all nodes in the power grid Q ,in, Q It is expressed as: ； in, represents the i-th node weight coefficient of node j, .

3. The global collaborative governance method for wide-area power quality disturbance problems according to claim 2 is characterized in that: The specific process of calculating the weighted voltage distortion rate of each node based on the node weight coefficient matrix is ​​as follows: B1. Calculate the weight of each node weight coefficient in the node weight coefficient matrix based on the entropy weight method, where i The weight of the node weight coefficient is expressed as: ； in, Indicates i The entropy of the node weight coefficients, ; B2. Calculate the weight of each node based on the weight of each node weight coefficient, and the weight of the node is expressed as: ； in, represents the weight of node j; B3. Calculate the harmonic voltage of each node: ； Where h represents the harmonic order, represents the effective value of the harmonic voltage of the hth harmonic at the node, represents the admittance matrix of the hth harmonic, Represents the effective value of the harmonic current of the hth harmonic passing through the node; B4. Calculate the voltage distortion rate of each node: ； in, Representation Node j The voltage distortion rate, Representation Node j The effective value of the fundamental voltage, For Node j The hth harmonic voltage effective value; B5. Based on the weight of each node and the voltage distortion rate, a weighted voltage distortion rate of each node is constructed. The weighted voltage distortion rate is expressed as: ； in, Representation Node j The weighted voltage distortion rate.

4. The global collaborative governance method for wide-area power quality disturbance problems according to claim 3 is characterized in that: The governance scheme optimization model is expressed as: ； in, represents the preset total governance investment, Indicates whether the node j Adopt governance solutions, Indicates that at the node j Adopt governance solutions, Indicates that at the node j No governance solution is adopted. For Node j The investment required for governance solutions.

5. The global collaborative governance method for wide-area power quality disturbance problems according to claim 1 is characterized in that: The specific process of using the Golden Jackal optimization algorithm to solve the governance scheme optimization model and obtain the optimal governance node and the corresponding governance scheme is as follows: C1. Initialize the population position and the maximum number of iterations, wherein the population includes a plurality of prey, each prey consisting of whether a node in the power grid adopts a governance scheme and a corresponding governance scheme; C2. Calculate the fitness value of each individual's location, select the prey with the best fitness value as the male golden jackal, and the prey with the second best fitness value as the female golden jackal; C3. In the exploration phase, the positions of the male and female golden jackals are updated. The specific formula is: ； ； , represent the positions of male and female golden jackals at the tth iteration, , They represent the updated positions of the male and female golden jackals in the tth iteration, E is the escape energy of the prey, is a random number between 0 and 1. represents a random number based on Levy distribution, Represents the prey matrix of the tth iteration ; C4. When the escape energy E of the prey is less than 1, the development phase is entered, and the positions of the male and female golden jackals are updated. The specific formula is: ； ； C5. When the maximum number of iterations is reached When , the position of the golden jackal is updated and the updated position of the golden jackal is used as the output, where the update formula of the position of the golden jackal is: ; in, yes The position of the male golden jackal in the iteration, yes The position of the female golden jackal in iteration .

6. A global collaborative governance method for wide-area power quality disturbance problems according to claim 5, characterized in that: The prey matrix It is expressed as: ； in, represents the kth prey, , represents the solution to the e-th dimension optimization problem for the k-th prey, , , d is the dimension of the optimization problem to be solved, and p is the number of prey in the population.

7. A global collaborative governance method for wide-area power quality disturbance problems according to claim 5, characterized in that: The escape energy E of the prey is expressed as: ； in, represents the decreasing process of prey energy, represents the initial state of prey energy, is a random number in the range of 0-1, is a constant, and t is the current number of iterations.

8. The global collaborative governance method for wide-area power quality disturbance problems according to claim 5 is characterized in that: Random numbers based on Levy distribution The specific expression is: ； in, is the Levy flight function, , is the intermediate parameter, and is a random number ranging from 0 to 1. is a constant, , represents the Gamma function, .

9. A global collaborative governance system for wide-area power quality disturbance problems, used to execute a global collaborative governance method for wide-area power quality disturbance problems as described in any one of claims 1 to 8, characterized in that: include: The node weight coefficient matrix construction module is used to construct a node weight coefficient matrix based on all nodes in the power grid, taking into account power quality transfer characteristics, responsibility contribution, power quality sensitive loads and governance cost indicators; A weighted voltage distortion rate calculation module, used to calculate the weighted voltage distortion rate of each node based on the node weight coefficient matrix; A governance scheme optimization model building module is used to build a governance scheme optimization model with the governance node selection and the adopted governance scheme as the optimization objects and the minimum sum of weighted voltage distortion rates of all nodes in the power grid as the objective function; The solution module is used to solve the governance scheme optimization model by using the golden jackal optimization algorithm to obtain the best governance node and the governance scheme corresponding to the best governance node.

10. A data processing device, characterized in that: include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements a global collaborative management method for wide-area power quality disturbance problems as described in any one of claims 1 to 8 when executing the computer program.

Images