CN110458314A - A kind of load prediction data decomposition method for power grid Tidal forecasting a few days ago - Google Patents

A kind of load prediction data decomposition method for power grid Tidal forecasting a few days ago Download PDF

Info

Publication number
CN110458314A
CN110458314A CN201910234163.9A CN201910234163A CN110458314A CN 110458314 A CN110458314 A CN 110458314A CN 201910234163 A CN201910234163 A CN 201910234163A CN 110458314 A CN110458314 A CN 110458314A
Authority
CN
China
Prior art keywords
load
active
region
prediction data
main transformer
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910234163.9A
Other languages
Chinese (zh)
Other versions
CN110458314B (en
Inventor
李铁
梁晓赫
姜枫
冯占稳
何晓洋
蔡壮
张凯
姜狄
吴志琪
王亮
朱伟峰
何超军
詹克明
李峰
崔岱
李典阳
张宇时
汤磊
王鹏
刘永锋
刘娟
王磊
曾辉
高梓济
常荣明
张建
孙晨光
唐俊刺
孙明一
王淼
孙文涛
胡景锦
韩秋
王明凯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
State Grid Liaoning Electric Power Co Ltd
Electric Power Research Institute of State Grid Liaoning Electric Power Co Ltd
Beijing King Star Hi Tech System Control Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Liaoning Electric Power Co Ltd
Electric Power Research Institute of State Grid Liaoning Electric Power Co Ltd
Beijing King Star Hi Tech System Control Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, State Grid Liaoning Electric Power Co Ltd, Electric Power Research Institute of State Grid Liaoning Electric Power Co Ltd, Beijing King Star Hi Tech System Control Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201910234163.9A priority Critical patent/CN110458314B/en
Publication of CN110458314A publication Critical patent/CN110458314A/en
Application granted granted Critical
Publication of CN110458314B publication Critical patent/CN110458314B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention belongs to technical field of power systems, in particular to a kind of load prediction data decomposition method for power grid Tidal forecasting a few days ago, is a kind of 110kV substation load data decomposition method for Tidal forecasting a few days ago.Before every end of day, read in the total active prediction data of each 220kV transforming plant main transformer of second day, and it is based on current electric grid model and operating status, on the load that the total burden with power prediction data of 220kV main transformer is decomposed to junior 110kV and the 35kV substation of main transformer institute band, calculate the idle of each load.Decomposition computation result is used for the Tidal forecasting a few days ago of power grid by the present invention, and for power grid, Security Checking, idle pressure optimization a few days ago provide basic data a few days ago, improves grid stability and quality of voltage.

Description

A kind of load prediction data decomposition method for power grid Tidal forecasting a few days ago
Technical field
The invention belongs to technical field of power systems, in particular to a kind of load prediction for power grid Tidal forecasting a few days ago Data resolving method, specifically a kind of 110kV substation load data decomposition method for Tidal forecasting a few days ago.
Background technique
In Operation of Electric Systems, for the ease of schedule ahead grid operation mode, scheduling mode personnel need to mention daily The preceding establishment that operation of power networks scheme will be carried out to upcoming tomorrow, that is, work out operational plan a few days ago.Operational plan a few days ago is electricity The important evidence of scheduling is netted, reasonably management and running plan a few days ago is related to the safe and economic operation in system future.With power grid The expansion of scale and the method for operation complicate, and the plan security check a few days ago based on DC power flow had been unable to satisfy practical need in the past Ask, it is desirable that a few days ago plan carry out include static security, dynamic stability and transient stability etc. comprehensive security check, be This, needs to generate the reasonable AC power flow for meeting and planning a few days ago, and the work completed is needed here it is Tidal forecasting a few days ago.It is a few days ago damp Stream forecast just refers to basis power plant's Active Generation plan a few days ago, a few days ago maintenance plan and a few days ago bus load prediction data, generates AC power flow solution a few days ago, since power plants generating electricity plan a few days ago and bus load prediction data are generally the 96 of 15 minutes intervals and count According to, therefore Tidal forecasting will also generate corresponding 96 AC power flow solutions a few days ago.Traditional method be on the basis of DC power flow with Experience provides generator terminal voltage, then carries out AC power flow calculating, poor astringency, and result is often unreasonable, influences The correctness of Security Checking.
Lin Yi, Sun Hongbin " a few days ago in plan security check plan trend Auto " (Automation of Electric Systems, The 20th phase of volume 36 in October, 2012, pp.68-73) in propose a kind of new Tidal forecasting method in planning a few days ago, The secondary daily planning trend of power grid can be calculated, the operation of power networks state and safe feelings of next day are assessed by plan trend Condition.Former PROBLEM DECOMPOSITION is that active power adjustment subproblem and reactive voltage are distributed subproblem substep and solve by the plan trend method.It is logical Cross solve active power adjustment optimization subproblem coordinate it is inconsistent between a variety of planning datas, by solve reactive voltage be distributed son ask Topic avoids to determine reasonable generator terminal voltage using typical set end voltage bring convergence problem.Using being based on Prediction-correction step modern interior point method solves above-mentioned optimization subproblem, has good convergence.This method is at present in state Interior some large regional grid control centres and provincial power grid dispatching center are applied.
The method provided in document mainly solves the active and idle of generator in Tidal forecasting, in Tidal forecasting Active reactive value on middle load bus is data that are given, providing from bus load forecast system.In provincial power network tune In the actual application at degree center, bus load prediction data mainly provides total burden with power prediction data of 220kV main transformer. As dispatching of power netwoks operational management in recent years is intensive, the change of flattening, saves and adjust in electric network model gradually by 220kV bus institute Subordinate load 110kV, 35kV electric network model of band also increases into, in Tidal forecasting, in order to provide complete exchange tide Stream solution, needs to decompose total burden with power forecast data of 220kV main transformer in junior's network load of its band.The present invention In order to solve this problem the method for proposition is exactly.
It is right in the operation of China's dispatching of power netwoks at present the present invention relates to the operation of power networks region of 220kV-110kV-35kV 330kV, 220kV power grid below, mostly use radiation mode to run, i.e., the junior of the every substation 330kV/220kV band 110kV/35kV substation forms independent region, interknits in region, the area, junior of different 330kV, 220kV institute bands There is no electrical link between domain.This typical schematic diagram in region is as shown in annex map 3.
The present invention relates to the new sensitivity of operation of power networks, the physical significance of burden with power sensitivity is on certain bus After increase injection unit is active, the active variation of each main transformer in power grid.Sun Hongbin, Zhang Baiming, Xiang Niande are " quasi-stationary sensitive Spend analysis method " new sensitivity method is proposed in (Proceedings of the CSEE, in April, 1999 V19N4, pp.9-13), Different from conventional static Sensitivity Analysis Method, new sensitivity method considers the quasi-stationary physics of electric system and rings It answers, total variation between the new and old stable state of meter and system control front and back effectively increases the precision of sensitivity analysis.This method is based on electricity The PQ Decoupled Model of Force system, when generator is equipped with automatic voltage regulator (AVR), it is believed that the generator node is PV Node;And when generator adjusts (AQR) or automatic power factor adjusting (APFR) equipped with automatic reactive power, it is believed that the hair Motor node identical as common load bus is PQ node.In addition, considering static load characteristics at the one of node voltage Secondary or conic section.The tide model established so just naturally takes in these quasi-stationary physical responses, thus It is quasi-stationary sensitivity based on the calculated sensitivity of tide model.Above-mentioned new sensitivity is all made of in calculating herein Method.
Summary of the invention
In order to overcome shortcoming existing for above-mentioned prior art, the present invention proposes that one kind is used for power grid Tidal forecasting a few days ago Load prediction data decomposition method.The present invention uses total active prediction data of 220kV transforming plant main transformer in provincial power network, Based on grid operation mode and topological structure, calculating is decomposited in junior 110kV and the 35kV substation that each main transformer is connected and is born The active prediction data of lotus, and calculate corresponding idle prediction data.The result of decomposition computation is used for the trend a few days ago of power grid Forecast, and then power grid Security Checking, a few days ago idle pressure optimization offer basic data a few days ago are provided, improve grid stability and electricity Press the goal of the invention of quality.
In order to achieve the above-mentioned object of the invention, the present invention is achieved through the following technical solutions:
A kind of load prediction data decomposition method for power grid Tidal forecasting a few days ago reads in second before every end of day Day the total active prediction data of each 220kV transforming plant main transformer, and current electric grid model and operating status are based on, by 220kV main transformer On the load for junior 110kV and the 35kV substation that total burden with power prediction data decomposes main transformer institute band, and further calculate Each load it is idle;From the point of view of 110kV substation, comprising the following steps:
T at the time of step 1. presets daily progress decomposition computation, when T is usually daily 22;
Step 2. T at the time of being calculated daily comes temporarily, reads in current electric grid mould from energy management system EMS Type and calculation of tidal current form the region Z of junior's power grid of 220kV substation x institute bandx, x is 220kV substation in power grid Number, initial value 1;
Step 3. is to region Zx, calculate its internal active extraction factor matrix A of loadx
Step 4. is to region Zx, according to the active prediction data of its 220kV main transformer, calculate its internal 110kV's and 35kV The active prediction data of load;
Step 5. is to region Zx, according to the active prediction data of the load of its 110kV and 35kV, calculate 110kV's and 35kV Idle prediction data;
Step 6.x value is incremented by 1, and return step 2 continues to calculate the region of the next substation 220kV band, until all 220kV substation calculates completion.
The ZxThe step of generation, is as follows:
Step 2.1 automatically generates junior's Grid of x institute, 220kV substation band according to topological structure of electric, generation Regional model are as follows:
Wherein,For the 220kV main transformer in the substation, n platform is amounted to;For substation's band Load in the station junior power grid 110kV, amounts to m;To be the substation with negative in the station junior power grid 35kV Lotus amounts to k;
Step 2.2 read-in area ZxIn, each object at power grid peak load moment on the same day is corresponding active and reactive are as follows:
Wherein,For the active and reactive value in 220kV main transformer high-pressure side in region;For junior's power grid The active and reactive value of load in the station 110kV;For the active and reactive value of load in the station junior power grid 35kV;I is Several equipment.
The step 3. is to region Zx, calculate its internal active extraction factor matrix A of loadx, comprise the following steps that
Step 3.1 is to region Zx, its region internal loading is calculated to the active po wer sensitivity matrix of 220kV main transformer high-pressure side winding Sx, it is as follows:
Wherein,For the sensitivity matrix that 110kV load in region is active to main transformer (m*n),It is that i-th of 110kV load is active to the quasi-steady state spirit that j-th of main transformer high-pressure side winding in region is active in region Sensitivity, physical significance are after i-th of 110kV load increase unit is active, and j-th of 220kV main transformer high-pressure side winding is active Variable quantity;I is that load is less than or equal to m, becomes based on j and is less than or equal to n;
Similarly,For the sensitivity matrix (k* that 35kV load in region is active to main transformer N dimension);Matrix SxTotal dimension are as follows: (m+k) * n;
Step 3.2 is to region Zx, its region internal loading is calculated to the active extraction factor square of 220kV main transformer high-pressure side winding Battle array Ax,
Wherein,For 110kV load in region it is active to main transformer draw matrix (m*n dimension), element are as follows:
WhereinFor in formula (2), region ZxMiddle jth platform main transformer high-pressure side is active;For in formula (2), region ZxIn I 110kV load is active;For submatrixElement;
Similarly,For 35kV load in region it is active to main transformer draw matrix (k*n dimension), element are as follows:
WhereinFor in formula (2), region ZxMiddle jth platform main transformer high-pressure side is active;For in formula (2), region ZxIn I 35kV load is active;For submatrixElement.
The step 4. is to region Zx, according to the active prediction data of its 220kV main transformer, calculate its internal 110kV and The active prediction data of the load of 35kV, comprises the following steps that
Step 4.1 reads in the region Z at t-th of moment of second day from energy management system EMSxIn each 220kV master The active prediction data become, is denoted as:
Wherein, t is the markers of prediction data, and value range is 1~T, i.e. the second day prediction number that adds up to T moment According to ordinary circumstance T=96, t initial value is 1;For the prediction data of t moment;
The active prediction data of load of 110kV and 35kV inside step 4.2 zoningIt is as follows:
WhereinFor be calculated second day t-th when Carve region ZxThe active prediction data of middle 110kV and 35kV load;
Step 4.3t value is incremented by 1, and return step 3.1 continues to calculate subsequent time, until second day all moment is active Prediction data calculates completion.
The step 5. is to region Zx, according to the active prediction data of the load of its 110kV and 35kV, calculate 110kV and The idle prediction data of 35kV, comprises the following steps that
Step 5.1 is calculated to second day t moment, using constant power factoring method with reference to the power factor of current time load 110kV and 35kV reactive load prediction data in regionIt is denoted as:
Wherein:
In formula (10)For power grid same day peak load moment for being provided in formula (2) The active and reactive value of load;The i of formula 10 corresponds respectively to m, k of formula 9;
Step 5.2t value is incremented by 1, and return step 4.1 continues to calculate subsequent time, until second day all moment is idle Prediction data calculates completion.
The total active prediction data of each 220kV transforming plant main transformer for reading in second day, and based on current electric grid model and The total burden with power prediction data of 220kV main transformer is decomposed junior 110kV and the 35kV substation of main transformer institute band by operating status On load, and further calculate the idle of each load;With 220kV stand in main transformer calculated, include two in the subregion 220kV main transformer Tr1, Tr2, load Ld1, Ld2, Ld3, Ld4, Ld5, Ld6, Ld7, Ld8, Ld9;The following steps are included:
T at the time of step 1. presets daily progress decomposition computation, when T is usually daily 22;
Step 2. T at the time of being calculated daily comes temporarily, reads in current electric grid mould from energy management system EMS Type and calculation of tidal current form the region Z of junior's power grid of 220kV substation x institute bandx, x is 220kV substation in power grid Number, initial value 1;
Step 3. is to region Zx, calculate its internal active extraction factor matrix A of loadx
Step 4. is to region Zx, according to the active prediction data of its 220kV main transformer, calculate its internal 110kV's and 35kV The active prediction data of load;
Step 5. is to region Zx, according to the active prediction data of the load of its 110kV and 35kV, calculate 110kV's and 35kV Idle prediction data;
Step 6.x value is incremented by 1, and return step 2 continues to calculate the region of the next substation 220kV band, until all 220kV substation calculates completion.
The generation ZxThe step of it is as follows:
Step 2.1 automatically generates junior's Grid of x institute, 220kV substation band according to topological structure of electric, generation Regional model are as follows:
Wherein,For the 220kV main transformer in the substation, n platform is amounted to;For substation's band Load in the station junior power grid 110kV, amounts to m;To be the substation with negative in the station junior power grid 35kV Lotus amounts to k;
According to can be calculated: Z1={ Tr1, Ld1, Ld2, Ld3, Ld4, Ld5, Ld6 }
Z2={ Tr2, Ld7, Ld8, Ld9 }
Step 2.2 read-in area ZxIn, it is each object at power grid peak load moment on the same day corresponding active Px, max, idle Qx, max are as follows:
Wherein,For the active and reactive value in 220kV main transformer high-pressure side in region;For junior's power grid The active and reactive value of load in the station 110kV;For the active and reactive value of load in the station junior power grid 35kV;I is to work as Preload;
It is read according to the model calculating actually calculated:
The step 3. is to region Zx, calculate its internal active extraction factor matrix A of loadx, steps are as follows:
Step 3.1 is to region Zx, its region internal loading is calculated to the active po wer sensitivity matrix of 220kV main transformer high-pressure side winding Sx, it is as follows:
Wherein,For the sensitivity matrix that 110kV load in region is active to main transformer (m*n),It is that i-th of 110kV load is active to the quasi-steady state spirit that j-th of main transformer high-pressure side winding in region is active in region Sensitivity, physical significance are after i-th of 110kV load increase unit is active, and j-th of 220kV main transformer high-pressure side winding is active Variable quantity;I is that load is less than or equal to m, becomes based on j and is less than or equal to n;
Similarly,For the sensitivity matrix (k*n that 35kV load in region is active to main transformer Dimension);Matrix SxTotal dimension are as follows: (m+k) * n;
The sensitivity being calculated are as follows:
Step 3.2 is to region Zx, its region internal loading is calculated to the active extraction factor square of 220kV main transformer high-pressure side winding Battle array Ax,
Wherein,For 110kV load in region it is active to main transformer draw matrix (m*n dimension), element are as follows:
WhereinFor in formula (2), region ZxMiddle jth platform main transformer high-pressure side is active;For in formula (2), region ZxIn I 110kV load is active;For submatrixElement;
Similarly,For 35kV load in region it is active to main transformer draw matrix (k*n dimension), element are as follows:
WhereinFor in formula (2), region ZxMiddle jth platform main transformer high-pressure side is active;For in formula (2), region ZxIn I 35kV load is active;For submatrixElement;
The active extraction factor being calculated are as follows:
The step 4. is to region Zx, according to the active prediction data of its 220kV main transformer, calculate its internal 110kV and The active prediction data of the load of 35kV, comprises the following steps that
Step 4.1 reads in the region Z at t-th of moment of second day from energy management system EMSxIn each 220kV master The active prediction data become, is denoted as:
Wherein, t is the markers of prediction data, and value range is 1~T, i.e. the second day prediction number that adds up to T moment According to ordinary circumstance T=96, t initial value is 1;
The active prediction data of load of 110kV and 35kV inside step 4.2 zoningIt is as follows:
WhereinFor the second day being calculated T moment region ZxThe active prediction data of middle 110kV and 35kV load;
4.3t value is incremented by 1, and return step 3.1 continues to calculate subsequent time, until the active prediction at second day all moment Data calculate completion;
Calculate 2 the data obtained of 1 the data obtained of region and region.
The step 5. is to region Zx, according to the active prediction data of the load of its 110kV and 35kV, calculate 110kV and The idle prediction data of 35kV, comprises the following steps that
Step 5.1 is calculated to second day t moment, using constant power factoring method with reference to the power factor of current time load 110kV and 35kV reactive load prediction data in regionIt is denoted as:
Wherein:
In formula (10)Load for the power grid same day peak load moment provided in formula (2) Active and reactive value;
Idle numerical value can equally be calculated to obtain;
Step 5.2t value is incremented by 1, and return step 4.1 continues to calculate subsequent time, until second day all moment is idle Prediction data calculates completion.
The features of the present invention and beneficial effect are:
In provincial power network Energy Management System, power prediction module generally only provides the total active pre- of 220kV main transformer Measured data cannot be used directly for the whole network Tidal forecasting comprising juniors' electric network models such as 110kV, 35kV and calculate.The present invention is saving Total active prediction data that 220kV transforming plant main transformer is used in grade power grid, is based on grid operation mode and topological structure, calculates and divides The active prediction data of load in junior 110kV and the 35kV substation that each main transformer is connected is solved, and calculates corresponding nothing Function prediction data.The result of decomposition computation be used for power grid Tidal forecasting a few days ago, and then for power grid a few days ago Security Checking, a few days ago without Function pressure optimization provides basic data, improves grid stability and quality of voltage.
Detailed description of the invention
Understand for the ease of those of ordinary skill in the art and implement the present invention, with reference to the accompanying drawing and specific embodiment The present invention is described in further detail, and the following examples are intended to illustrate the invention, it is to be understood that protection model of the invention It encloses and is not limited by the specific implementation.
Fig. 1 is method flow block diagram of the invention.
Fig. 2 is substation of embodiment of the present invention connection relationship diagram.
Fig. 3 is existing region electrical link figure.
Specific embodiment
The present invention proposes a kind of load prediction data decomposition method for power grid Tidal forecasting a few days ago, as shown in Figure 1, figure 1 is method flow block diagram of the invention.Before every end of day, each 220kV transforming plant main transformer for reading in second day always has this method Function prediction data, and it is based on current electric grid model and operating status, the total burden with power prediction data of 220kV main transformer is decomposed into master Become on the load of junior 110kV and the 35kV substation of institute's band, and further calculates the idle of each load.With 110kV substation For, method includes the following steps:
T at the time of step 1. presets daily progress decomposition computation, when T is usually daily 22;
Step 2. T at the time of being calculated daily comes temporarily, reads in current electric grid mould from energy management system EMS Type and calculation of tidal current form the region Z of junior's power grid of 220kV substation x institute bandx, x is 220kV substation in power grid Number, initial value 1.Generate ZxThe step of it is as follows:
2.1 automatically generate junior's Grid of x institute, 220kV substation band, the region of generation according to topological structure of electric Model are as follows:
Wherein,For the 220kV main transformer in the substation, n platform is amounted to;For substation's band Load in the station junior power grid 110kV, amounts to m;To be the substation with negative in the station junior power grid 35kV Lotus amounts to k.
2.2 read-in area ZxIn, each object at power grid peak load moment on the same day is corresponding active and reactive are as follows:
Wherein,For the active and reactive value in 220kV main transformer high-pressure side in region;For junior's power grid The active and reactive value of load in the station 110kV;For the active and reactive value of load in the station junior power grid 35kV;I is Several equipment.
Step 3. is to region Zx, calculate its internal active extraction factor matrix A of loadx, steps are as follows:
3.1 couples of region Zx, its region internal loading is calculated to the active po wer sensitivity matrix S of 220kV main transformer high-pressure side windingx, such as Under:
Wherein,For the sensitivity matrix that 110kV load in region is active to main transformer (m*n),It is that i-th of 110kV load is active to the quasi-steady state spirit that j-th of main transformer high-pressure side winding in region is active in region Sensitivity, physical significance are after i-th of 110kV load increase unit is active, and j-th of 220kV main transformer high-pressure side winding is active Variable quantity;I is that load is less than or equal to m, becomes based on j and is less than or equal to n.
Similarly,For the sensitivity matrix (k* that 35kV load in region is active to main transformer N dimension).Matrix SxTotal dimension are as follows: (m+k) * n.
3.2 couples of region Zx, its region internal loading is calculated to the active extraction factor matrix A of 220kV main transformer high-pressure side windingx,
Wherein,For 110kV load in region it is active to main transformer draw matrix (m*n dimension), element are as follows:
WhereinFor in formula (2), region ZxMiddle jth platform main transformer high-pressure side is active;For in formula (2), region ZxIn I 110kV load is active;For submatrixElement.
Similarly,For 35kV load in region it is active to main transformer draw matrix (k*n dimension), element are as follows:
WhereinFor in formula (2), region ZxMiddle jth platform main transformer high-pressure side is active;For in formula (2), region ZxIn I 35kV load is active;For submatrixElement.
Step 4. is to region Zx, according to the active prediction data of its 220kV main transformer, calculate its internal 110kV's and 35kV The active prediction data of load, steps are as follows:
4.1 read in the region Z at t-th of moment of second day from energy management system EMSxIn each 220kV main transformer Active prediction data, is denoted as:
Wherein, t is the markers of prediction data, and value range is 1~T, i.e. the second day prediction number that adds up to T moment According to ordinary circumstance T=96, t initial value is 1;For the prediction data of t moment.
The active prediction data of load of 110kV and 35kV inside 4.2 zoningsIt is as follows:
WhereinFor be calculated second day t-th when Carve region ZxThe active prediction data of middle 110kV and 35kV load.
4.3t value is incremented by 1, and return step 3.1 continues to calculate subsequent time, until the active prediction at second day all moment Data calculate completion.
Step 5. is to region Zx, according to the active prediction data of the load of its 110kV and 35kV, calculate 110kV's and 35kV Idle prediction data, steps are as follows:
5.1 pairs of second day t moments, using constant power factoring method, with reference to the power factor of current time load, zoning In 110kV and 35kV reactive load prediction dataIt is denoted as:
Wherein:
In formula (10)For power grid same day peak load moment for being provided in formula (2) The active and reactive value of load.The i of formula 10 corresponds respectively to m, k of formula 9.
5.2t value is incremented by 1, and return step 4.1 continues to calculate subsequent time, until the idle prediction at second day all moment Data calculate completion.
Step 6.x value is incremented by 1, and return step 2 continues to calculate the region of the next substation 220kV band, until all 220kV substation calculates completion.
The working principle of the method for the present invention is:
By carrying out new sensitivity calculating to existing province's tune regional model, main transformer in the current station 220kV is obtained With the sensitivity between the load of the station 110kV, 35kV, the burden with power prediction data at the station 220kV is read, is gone out by Calculation of Sensitivity The load prediction data at the station 110kV, 35kV carries out Tidal forecasting meter by decomposing obtained load data in Tidal forecasting It calculates.
Embodiment 2.
The present embodiment is to calculate the main transformer in a station 220kV, the present embodiment station station connection relationship such as Fig. 2 institute Show, 2 220kV main transformer Tr1, Tr2, load Ld1, Ld2, Ld3, Ld4, Ld5, Ld6, Ld7, Ld8, Ld9 are included in the subregion;
The present invention proposes a kind of load prediction data decomposition method for power grid Tidal forecasting a few days ago, including following step It is rapid:
T at the time of step 1. presets daily progress decomposition computation, when T is usually daily 22;
Step 2. T at the time of being calculated daily comes temporarily, reads in current electric grid mould from energy management system EMS Type and calculation of tidal current form the region Z of junior's power grid of 220kV substation x institute bandx, x is 220kV substation in power grid Number, initial value 1.Generate ZxThe step of it is as follows:
2.1 automatically generate junior's Grid of x institute, 220kV substation band, the region of generation according to topological structure of electric Model are as follows:
Wherein,For the 220kV main transformer in the substation, n platform is amounted to;For substation's band Load in the station junior power grid 110kV, amounts to m;To be the substation with negative in the station junior power grid 35kV Lotus amounts to k.
According to can be calculated: Z1={ Tr1, Ld1, Ld2, Ld3, Ld4, Ld5, Ld6 }
Z2={ Tr2, Ld7, Ld8, Ld9 }
2.2 read-in area ZxIn, each object at power grid peak load moment on the same day corresponding active Px, max, idle Qx, Max are as follows:
Wherein,For the active and reactive value in 220kV main transformer high-pressure side in region;For junior's power grid The active and reactive value of load in the station 110kV;For the active and reactive value of load in the station junior power grid 35kV;I is to work as Preload;
It is read according to the model calculating actually calculated:
Step 3. is to region Zx, calculate its internal active extraction factor matrix A of loadx, steps are as follows:
3.1 couples of region Zx, its region internal loading is calculated to the active po wer sensitivity matrix S of 220kV main transformer high-pressure side windingx, such as Under:
Wherein,For the sensitivity matrix that 110kV load in region is active to main transformer (m*n),It is that i-th of 110kV load is active to the quasi-steady state spirit that j-th of main transformer high-pressure side winding in region is active in region Sensitivity, physical significance are after i-th of 110kV load increase unit is active, and j-th of 220kV main transformer high-pressure side winding is active Variable quantity.I is that load is less than or equal to m, becomes based on j and is less than or equal to n.
Similarly,For the sensitivity matrix (k*n that 35kV load in region is active to main transformer Dimension).Matrix SxTotal dimension are as follows: (m+k) * n.
The sensitivity being calculated are as follows:
3.2 couples of region Zx, its region internal loading is calculated to the active extraction factor matrix A of 220kV main transformer high-pressure side windingx,
Wherein,For 110kV load in region it is active to main transformer draw matrix (m*n dimension), element are as follows:
WhereinFor in formula (2), region ZxMiddle jth platform main transformer high-pressure side is active;For in formula (2), region ZxIn I 110kV load is active;For submatrixElement.
Similarly,For 35kV load in region it is active to main transformer draw matrix (k*n dimension), element are as follows:
WhereinFor in formula (2), region ZxMiddle jth platform main transformer high-pressure side is active;For in formula (2), region ZxIn I 35kV load is active;For submatrixElement.
The active extraction factor being calculated are as follows:
Step 4. is to region Zx, according to the active prediction data of its 220kV main transformer, calculate its internal 110kV's and 35kV The active prediction data of load, steps are as follows:
4.1 read in the region Z at t-th of moment of second day from energy management system EMSxIn each 220kV main transformer Active prediction data, is denoted as:
Wherein, t is the markers of prediction data, and value range is 1~T, i.e. the second day prediction number that adds up to T moment According to ordinary circumstance T=96, t initial value is 1;
The active prediction data of load of 110kV and 35kV inside 4.2 zoningsIt is as follows:
WhereinFor the second day being calculated T moment region ZxThe active prediction data of middle 110kV and 35kV load.
4.3t value is incremented by 1, and return step 3.1 continues to calculate subsequent time, until the active prediction at second day all moment Data calculate completion.
Computed information is as follows: region 1
Region 2:
Step 5. is to region Zx, according to the active prediction data of the load of its 110kV and 35kV, calculate 110kV's and 35kV Idle prediction data, steps are as follows:
5.1 pairs of second day t moments, using constant power factoring method, with reference to the power factor of current time load, zoning In 110kV and 35kV reactive load prediction dataIt is denoted as:
Wherein:
In formula (10)Load for the power grid same day peak load moment provided in formula (2) Active and reactive value.
Idle numerical value can equally be calculated to obtain;
5.2 t values are incremented by 1, and return step 4.1 continues to calculate subsequent time, until the idle prediction at second day all moment Data calculate completion.
Step 6.x value is incremented by 1, and return step 2 continues to calculate the region of the next substation 220kV band, until all 220kV substation calculates completion.

Claims (10)

1. a kind of load prediction data decomposition method for power grid Tidal forecasting a few days ago, it is characterized in that: being read before every end of day Enter each 220kV transforming plant main transformer of second day always active prediction data, and is based on current electric grid model and operating status, it will The total burden with power prediction data of 220kV main transformer decomposes on the load of junior 110kV and the 35kV substation of main transformer institute band, goes forward side by side One step calculates the idle of each load;From the point of view of 110kV substation, comprising the following steps:
T at the time of step 1. presets daily progress decomposition computation, when T is usually daily 22;
Step 2. T at the time of being calculated daily come it is interim, from energy management system EMS read in current electric grid model and Calculation of tidal current forms the region Z of junior's power grid of 220kV substation x institute bandx, x is the volume of 220kV substation in power grid Number, initial value 1;
Step 3. is to region Zx, calculate its internal active extraction factor matrix A of loadx
Step 4. is to region Zx, according to the active prediction data of its 220kV main transformer, the load for calculating its internal 110kV and 35kV has Function prediction data;
Step 5. is to region Zx, according to the active prediction data of the load of its 110kV and 35kV, calculate the idle pre- of 110kV and 35kV Measured data;
Step 6.x value is incremented by 1, and return step 2 continues to calculate the region of the next substation 220kV band, until whole 220kV Substation calculates completion.
2. a kind of load prediction data decomposition method for power grid Tidal forecasting a few days ago according to claim 1, special Sign is: the ZxThe step of generation, is as follows:
Step 2.1 automatically generates junior's Grid of x institute, 220kV substation band, the region of generation according to topological structure of electric Model are as follows:
Wherein,For the 220kV main transformer in the substation, n platform is amounted to;For band junior, the substation Load in the station power grid 110kV, amounts to m;To be the substation with the load in the station junior power grid 35kV, always Count k;
Step 2.2 read-in area ZxIn, each object at power grid peak load moment on the same day is corresponding active and reactive are as follows:
Wherein,For the active and reactive value in 220kV main transformer high-pressure side in region;For junior power grid 110kV The active and reactive value of load in standing;For the active and reactive value of load in the station junior power grid 35kV;I is which sets It is standby.
3. a kind of load prediction data decomposition method for power grid Tidal forecasting a few days ago according to claim 1, special Sign is: the step 3. is to region Zx, calculate its internal active extraction factor matrix A of loadx, comprise the following steps that
Step 3.1 is to region Zx, its region internal loading is calculated to the active po wer sensitivity matrix S of 220kV main transformer high-pressure side windingx, such as Under:
Wherein,For the sensitivity matrix (m* that 110kV load in region is active to main transformer N),It is that i-th of 110kV load is active sensitive to j-th of the active quasi-steady state of main transformer high-pressure side winding in region in region Degree, physical significance are the active changes of j-th of 220kV main transformer high-pressure side winding after i-th of 110kV load increase unit is active Change amount;I is that load is less than or equal to m, becomes based on j and is less than or equal to n;
Similarly,For the sensitivity matrix (k*n that 35kV load in region is active to main transformer Dimension);Matrix SxTotal dimension are as follows: (m+k) * n;
Step 3.2 is to region Zx, its region internal loading is calculated to the active extraction factor matrix A of 220kV main transformer high-pressure side windingx,
Wherein,For 110kV load in region it is active to main transformer draw matrix (m*n dimension), element are as follows:
WhereinFor in formula (2), region ZxMiddle jth platform main transformer high-pressure side is active;For in formula (2), region ZxIn i-th 110kV load is active;For submatrixElement;
Similarly,For 35kV load in region it is active to main transformer draw matrix (k*n dimension), element are as follows:
WhereinFor in formula (2), region ZxMiddle jth platform main transformer high-pressure side is active;For in formula (2), region ZxIn i-th 35kV load is active;For submatrixElement.
4. a kind of load prediction data decomposition method for power grid Tidal forecasting a few days ago according to claim 1, special Sign is: the step 4. is to region Zx, according to the active prediction data of its 220kV main transformer, calculate its internal 110kV's and 35kV The active prediction data of load, comprises the following steps that
Step 4.1 reads in the region Z at t-th of moment of second day from energy management system EMSxIn each 220kV main transformer have Function prediction data, is denoted as:
Wherein, t is the markers of prediction data, and value range is 1~T, i.e. the second day prediction data that adds up to T moment, one As situation T=96, t initial value be 1;For the prediction data of t moment;
The active prediction data of load of 110kV and 35kV inside step 4.2 zoningIt is as follows:
WhereinFor t-th of the moment region Z of second day being calculatedxIn The active prediction data of 110kV and 35kV load;
Step 4.3t value is incremented by 1, and return step 3.1 continues to calculate subsequent time, until the active prediction at second day all moment Data calculate completion.
5. a kind of load prediction data decomposition method for power grid Tidal forecasting a few days ago according to claim 1, special Sign is: the step 5. is to region Zx, according to the active prediction data of the load of its 110kV and 35kV, calculate 110kV's and 35kV Idle prediction data, comprises the following steps that
Step 5.1 is to second day t moment, using constant power factoring method, with reference to the power factor of current time load, zoning In 110kV and 35kV reactive load prediction dataIt is denoted as:
Wherein:
In formula (10)Load for the power grid same day peak load moment provided in formula (2) Active and reactive value;The i of formula 10 corresponds respectively to m, k of formula 9;
Step 5.2t value is incremented by 1, and return step 4.1 continues to calculate subsequent time, until the idle prediction at second day all moment Data calculate completion.
6. a kind of load prediction data decomposition method for power grid Tidal forecasting a few days ago according to claim 1, special Sign is: the total active prediction data of each 220kV transforming plant main transformer for reading in second day, and is based on current electric grid model and operation The total burden with power prediction data of 220kV main transformer is decomposed the load of junior 110kV and the 35kV substation of main transformer institute band by state On, and further calculate the idle of each load;With 220kV stand in main transformer calculated, include two 220kV masters in the subregion Become Tr1, Tr2, load Ld1, Ld2, Ld3, Ld4, Ld5, Ld6, Ld7, Ld8, Ld9;The following steps are included:
T at the time of step 1. presets daily progress decomposition computation, when T is usually daily 22;
Step 2. T at the time of being calculated daily come it is interim, from energy management system EMS read in current electric grid model and Calculation of tidal current forms the region Z of junior's power grid of 220kV substation x institute bandx, x is the volume of 220kV substation in power grid Number, initial value 1;
Step 3. is to region Zx, calculate its internal active extraction factor matrix A of loadx
Step 4. is to region Zx, according to the active prediction data of its 220kV main transformer, the load for calculating its internal 110kV and 35kV has Function prediction data;
Step 5. is to region Zx, according to the active prediction data of the load of its 110kV and 35kV, calculate the idle pre- of 110kV and 35kV Measured data;
Step 6.x value is incremented by 1, and return step 2 continues to calculate the region of the next substation 220kV band, until whole 220kV Substation calculates completion.
7. a kind of load prediction data decomposition method for power grid Tidal forecasting a few days ago according to claim 6, special Sign is: the generation ZxThe step of it is as follows:
Step 2.1 automatically generates junior's Grid of x institute, 220kV substation band, the region of generation according to topological structure of electric Model are as follows:
Wherein,For the 220kV main transformer in the substation, n platform is amounted to;For band junior, the substation Load in the station power grid 110kV, amounts to m;To be the substation with the load in the station junior power grid 35kV, always Count k;
According to can be calculated: Z1={ Tr1, Ld1, Ld2, Ld3, Ld4, Ld5, Ld6 }
Z2={ Tr2, Ld7, Ld8, Ld9 }
Step 2.2 read-in area ZxIn, each object at power grid peak load moment on the same day corresponding active Px, max, idle Qx, Max are as follows:
Wherein,For the active and reactive value in 220kV main transformer high-pressure side in region;For junior power grid 110kV The active and reactive value of load in standing;For the active and reactive value of load in the station junior power grid 35kV;I is current negative Lotus;
It is read according to the model calculating actually calculated:
8. a kind of load prediction data decomposition method for power grid Tidal forecasting a few days ago according to claim 6, special Sign is: the step 3. is to region Zx, calculate its internal active extraction factor matrix A of loadx, steps are as follows:
Step 3.1 is to region Zx, its region internal loading is calculated to the active po wer sensitivity matrix S of 220kV main transformer high-pressure side windingx, such as Under:
Wherein,For the sensitivity matrix (m*n) that 110kV load in region is active to main transformer,It is the active new sensitivity active to j-th of main transformer high-pressure side winding in region of i-th of 110kV load in region, Its physical significance is the active variation of j-th of 220kV main transformer high-pressure side winding after i-th of 110kV load increase unit is active Amount;I is that load is less than or equal to m, becomes based on j and is less than or equal to n;
Similarly,For the sensitivity matrix active to main transformer of 35kV load in region (k*n dimension); Matrix SxTotal dimension are as follows: (m+k) * n;
The sensitivity being calculated are as follows:
Step 3.2 is to region Zx, its region internal loading is calculated to the active extraction factor matrix A of 220kV main transformer high-pressure side windingx,
Wherein,For 110kV load in region it is active to main transformer draw matrix (m*n dimension), element are as follows:
WhereinFor in formula (2), region ZxMiddle jth platform main transformer high-pressure side is active;For in formula (2), region ZxIn i-th 110kV load is active;For submatrixElement;
Similarly,For 35kV load in region it is active to main transformer draw matrix (k*n dimension), element are as follows:
WhereinFor in formula (2), region ZxMiddle jth platform main transformer high-pressure side is active;For in formula (2), region ZxIn i-th 35kV load is active;For submatrixElement;
Active extraction factor is calculated are as follows:
9. a kind of load prediction data decomposition method for power grid Tidal forecasting a few days ago according to claim 6, special Sign is: the step 4. is to region Zx, according to the active prediction data of its 220kV main transformer, calculate its internal 110kV's and 35kV The active prediction data of load, comprises the following steps that
Step 4.1 reads in the region Z at t-th of moment of second day from energy management system EMSxIn each 220kV main transformer have Function prediction data, is denoted as:
Wherein, t is the markers of prediction data, and value range is 1~T, i.e. the second day prediction data that adds up to T moment, one As situation T=96, t initial value be 1;
The active prediction data of load of 110kV and 35kV inside step 4.2 zoningIt is as follows:
WhereinFor be calculated second day t-th when Carve region ZxThe active prediction data of middle 110kV and 35kV load;
4.3t value is incremented by 1, and return step 3.1 continues to calculate subsequent time, until the active prediction data at second day all moment It calculates and completes;
Calculate 2 the data obtained of 1 the data obtained of region and region.
10. a kind of load prediction data decomposition method for power grid Tidal forecasting a few days ago according to claim 6, special Sign is: the step 5. is to region Zx, according to the active prediction data of the load of its 110kV and 35kV, calculate 110kV's and 35kV Idle prediction data, comprises the following steps that
Step 5.1 is to second day t moment, using constant power factoring method, with reference to the power factor of current time load, zoning In 110kV and 35kV reactive load prediction dataIt is denoted as:
Wherein:
In formula (10)The load at the power grid same day peak load moment to provide in formula (2) is active, Without work value;
Idle numerical value can equally be calculated to obtain;
Step 5.2t value is incremented by 1, and return step 4.1 continues to calculate subsequent time, until the idle prediction at second day all moment Data calculate completion.
CN201910234163.9A 2019-03-26 2019-03-26 Load prediction data decomposition method for power grid day-ahead tide prediction Active CN110458314B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910234163.9A CN110458314B (en) 2019-03-26 2019-03-26 Load prediction data decomposition method for power grid day-ahead tide prediction

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910234163.9A CN110458314B (en) 2019-03-26 2019-03-26 Load prediction data decomposition method for power grid day-ahead tide prediction

Publications (2)

Publication Number Publication Date
CN110458314A true CN110458314A (en) 2019-11-15
CN110458314B CN110458314B (en) 2023-07-25

Family

ID=68480598

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910234163.9A Active CN110458314B (en) 2019-03-26 2019-03-26 Load prediction data decomposition method for power grid day-ahead tide prediction

Country Status (1)

Country Link
CN (1) CN110458314B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112288598A (en) * 2020-12-24 2021-01-29 中国电力科学研究院有限公司 Method and system for determining composition of load element of transformer substation

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103823837A (en) * 2013-12-13 2014-05-28 国网安徽省电力公司 Key transmission section searching method based on fuzzy clustering and shortest path
WO2015081771A1 (en) * 2013-12-06 2015-06-11 国家电网公司 Adaptive emergency control method for voltage security and stability based on synchronous measurement information
CN105005940A (en) * 2015-07-09 2015-10-28 河海大学 Correlation-considered GEPOPF calculation method
CN105375513A (en) * 2015-11-06 2016-03-02 国家电网公司 Automatic 110kV wind power field voltage control method based on real-time online equivalence
CN107464048A (en) * 2017-07-26 2017-12-12 广东电网有限责任公司电力调度控制中心 A kind of plan security check method a few days ago based on research state
CN108462179A (en) * 2018-01-16 2018-08-28 国网江苏省电力有限公司电力科学研究院 Power Flow Tracing Method based on breadth first search
CN108631308A (en) * 2018-05-23 2018-10-09 国网天津市电力公司电力科学研究院 A kind of prediction technique of 500kV substations burden with power variation tendency

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015081771A1 (en) * 2013-12-06 2015-06-11 国家电网公司 Adaptive emergency control method for voltage security and stability based on synchronous measurement information
CN103823837A (en) * 2013-12-13 2014-05-28 国网安徽省电力公司 Key transmission section searching method based on fuzzy clustering and shortest path
CN105005940A (en) * 2015-07-09 2015-10-28 河海大学 Correlation-considered GEPOPF calculation method
CN105375513A (en) * 2015-11-06 2016-03-02 国家电网公司 Automatic 110kV wind power field voltage control method based on real-time online equivalence
CN107464048A (en) * 2017-07-26 2017-12-12 广东电网有限责任公司电力调度控制中心 A kind of plan security check method a few days ago based on research state
CN108462179A (en) * 2018-01-16 2018-08-28 国网江苏省电力有限公司电力科学研究院 Power Flow Tracing Method based on breadth first search
CN108631308A (en) * 2018-05-23 2018-10-09 国网天津市电力公司电力科学研究院 A kind of prediction technique of 500kV substations burden with power variation tendency

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
ZHANG S ETC.: ""An Improved Network Partition Method of Voltage Control Considering Coordinated Active and Reactive Power Sources"", 《2018 2ND INTERNATIONAL CONFERENCE ON ENVIRONMENTAL AND ENERGY ENGINEERING》 *
吕颖等: "智能电网调度控制系统的安全校核服务及实用化关键技术", 《电力系统自动化》 *
李传栋等: "潮流追踪解析算法", 《电力系统及其自动化学报》 *
谢开贵等: "电力系统功率分配的解析模型和算法", 《中国电机工程学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112288598A (en) * 2020-12-24 2021-01-29 中国电力科学研究院有限公司 Method and system for determining composition of load element of transformer substation

Also Published As

Publication number Publication date
CN110458314B (en) 2023-07-25

Similar Documents

Publication Publication Date Title
CN107069814B (en) The Fuzzy Chance Constrained Programming method and system that distribution distributed generation resource capacity is layouted
CN109948849B (en) Power distribution network frame planning method considering energy storage access
CN108695857B (en) Automatic voltage control method, device and system for wind power plant
CN104751246A (en) Active distribution network planning method based on stochastic chance constraint
CN112653154B (en) Distributed photovoltaic power distribution network reactive power optimization control method based on edge calculation
CN106026113A (en) Micro-grid system monitoring method having reactive automatic compensation function
CN110543696B (en) Method for small unmodeled unit to participate in electric power market clearing and safety check
CN113489003B (en) Source network coordination planning method considering wind-light-water integrated complementary operation
CN105375513A (en) Automatic 110kV wind power field voltage control method based on real-time online equivalence
CN111461919A (en) Wind-powered electricity generation field power control integration monitored control system
CN115622053B (en) Automatic load modeling method and device for considering distributed power supply
CN106779313A (en) Multiple target distributed power source addressing constant volume method based on mixed integer programming
CN109713716A (en) A kind of chance constraint economic load dispatching method of the wind-electricity integration system based on security domain
CN109214713A (en) Active distribution network planing method containing distributed generation resource
CN109524988A (en) A kind of wind-powered electricity generation based on total active power trend prediction collects station voltage control method
CN115811070A (en) Flywheel energy storage self-adaptive capacity configuration method and system for assisting thermal power generating unit in frequency modulation
CN105958530A (en) Microgrid system with reactive power automatic compensation function
CN110458314A (en) A kind of load prediction data decomposition method for power grid Tidal forecasting a few days ago
CN105071397A (en) Coordinated reactive voltage control method of different reactive compensation devices of wind power delivery
CN104809543A (en) Power grid operation mode generating method based on monthly power transmission and distribution equipment maintenance plan
CN105207255B (en) A kind of power system peak regulation computational methods suitable for wind power output
Zhang et al. Research on active distribution network structure planning for multi distributed generation access
CN105977993B (en) A kind of reactive-load compensation method of the intelligent distribution system based on load
Cheng et al. Analysis of multi-scenario power supply and demand balance in Shandong power grid based on the new generation PSDB platform
CN110619436A (en) Active power distribution network planning method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant