[go: up one dir, main page]

CN108233373B - A Probabilistic Harmonic Analysis Method for Distributed Photovoltaic Access to Distribution Network Considering Weather Scenarios - Google Patents

A Probabilistic Harmonic Analysis Method for Distributed Photovoltaic Access to Distribution Network Considering Weather Scenarios Download PDF

Info

Publication number
CN108233373B
CN108233373B CN201711207394.8A CN201711207394A CN108233373B CN 108233373 B CN108233373 B CN 108233373B CN 201711207394 A CN201711207394 A CN 201711207394A CN 108233373 B CN108233373 B CN 108233373B
Authority
CN
China
Prior art keywords
harmonic
photovoltaic
distribution network
harmonic voltage
power
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.)
Active
Application number
CN201711207394.8A
Other languages
Chinese (zh)
Other versions
CN108233373A (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.)
China Electric Power Research Institute Co Ltd CEPRI
Economic and Technological Research Institute of State Grid Ningxia Electric Power Co Ltd
Original Assignee
China Electric Power Research Institute Co Ltd CEPRI
Economic and Technological Research Institute of State Grid Ningxia Electric Power 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 China Electric Power Research Institute Co Ltd CEPRI, Economic and Technological Research Institute of State Grid Ningxia Electric Power Co Ltd filed Critical China Electric Power Research Institute Co Ltd CEPRI
Priority to CN201711207394.8A priority Critical patent/CN108233373B/en
Publication of CN108233373A publication Critical patent/CN108233373A/en
Application granted granted Critical
Publication of CN108233373B publication Critical patent/CN108233373B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for AC mains or AC distribution networks
    • H02J3/01Arrangements for reducing harmonics or ripples
    • H02J3/383
    • 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/50Photovoltaic [PV] energy
    • Y02E10/56Power conversion systems, e.g. maximum power point trackers
    • 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/40Arrangements for reducing harmonics

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明提出的考虑天气场景的分布式光伏接入配电网的概率谐波分析方法,首先由典型天气场景下的分布式光伏系统输出功率曲线通过参数估计方法得到其Beta分布的形状参数,确定光伏系统输出功率的概率分布模型,并进行抽样,得到分布式光伏输出功率值,进行谐波潮流计算,得到配电网的谐波注入电流,进而得到系统各节点的谐波电压,当完成所有抽样次数后,计算得到配电网的谐波电压概率分布和谐波畸变率THD。

Figure 201711207394

The method for probabilistic harmonic analysis of distributed photovoltaics connected to the distribution network considering weather scenarios proposed by the present invention firstly obtains the shape parameters of its Beta distribution from the output power curve of the distributed photovoltaic system under typical weather scenarios through a parameter estimation method, and determines The probability distribution model of the output power of the photovoltaic system, and sampling to obtain the distributed photovoltaic output power value, carry out the harmonic power flow calculation, obtain the harmonic injection current of the distribution network, and then obtain the harmonic voltage of each node of the system. After sampling times, the harmonic voltage probability distribution and harmonic distortion rate THD of the distribution network are calculated.

Figure 201711207394

Description

Probability harmonic analysis method considering weather scene for distributed photovoltaic access power distribution network
Technical Field
The invention relates to the technical field of photovoltaic power generation, in particular to a probability harmonic analysis method for a distributed photovoltaic access power distribution network considering a weather scene.
Background
The distributed power supply has the advantages of environmental protection, high efficiency, safety, flexibility, high cost performance, short construction period and small risk, can be used as supplement of centralized energy, is coordinated with a large power grid for development, can save investment, reduce energy consumption and improve the reliability and flexibility of a power system, and becomes an important direction of a world energy strategy.
Common distributed power sources include wind power generation, solar power generation, small hydropower, micro gas turbines, fuel cells, and the like. Most distributed power generation facilities require power frequency ac power to be delivered to the grid or to the users through utility power electronics. These power electronics devices may cause distortion of the current and voltage waveforms of the power grid, resulting in harmonic pollution of the power grid.
On the other hand, the random change of the wind speed and the illumination intensity is the natural property of renewable energy, so that the electric energy generated by wind power generation and photovoltaic power generation also changes randomly, and further, the harmonic waves generated by the distributed power supply also have randomness. Therefore, it is necessary to perform calculation and analysis by using a probability method. Meteorological conditions have a large influence on the photovoltaic output power, so that the harmonic distortion condition of the distributed photovoltaic access power distribution network is analyzed respectively according to different weather meteorological conditions.
Disclosure of Invention
There is a need to provide a probability harmonic analysis method for a distributed photovoltaic access distribution network considering weather scenarios.
A probability harmonic analysis method for a distributed photovoltaic access power distribution network considering weather scenes is characterized by comprising the following steps:
the first step is as follows: photovoltaic power generation output power P of distributed photovoltaic access power distribution network by utilizing K-MEANS clustering analysis methodpvCounting and clustering to obtain a photovoltaic output power curve and distribution probability of the photovoltaic output power curve in a typical weather scene;
the second step is that: according to the photovoltaic output power curve of the first step of typical weather scene, performing parameter estimation by adopting a least square method, and establishing a photovoltaic power generation power output probability density function of the photovoltaic power station under different weather scenes:
Figure GDA0001672217950000021
in the formula, PPVFor photovoltaic power generation output power, PmaxAlpha and Beta are shape parameters of Beta distribution function, and are the maximum value of the output power of the photovoltaic power generation;
the third step: based on the photovoltaic power generation power output probability density function, the Monte Carlo method is adopted to respectively sample the photovoltaic power generation power under various typical weather scenes, and the photovoltaic output power sampling value under the weather scene is obtained through sampling.
The fourth step: according to the sampling value of the photovoltaic power generation power, carrying out harmonic load flow calculation on the sampling value of the photovoltaic output power of the power distribution network to obtain each subharmonic current and harmonic voltage value injected into the power distribution network;
the fifth step: if the sampling times are less than the sampling set times, returning to the third step to continue sampling calculation, and calculating the average value of the harmonic voltages injected into the power distribution network, otherwise, performing the next step;
and a sixth step: substituting the harmonic voltage mean value of each harmonic into the following harmonic voltage distortion rate formula to calculate the harmonic voltage distortion rate THD of each harmonic under various typical weather scenesu
Harmonic voltage distortion rate formula:
Figure GDA0001672217950000022
in the formula: uh: the harmonic voltage amplitude of the h harmonic; u shape1: a fundamental voltage amplitude; n: the highest order of the harmonic considered;
the seventh step: after the harmonic voltage distortion rate analysis and calculation under various typical weather scenes are completed, the probabilities and the harmonic voltage distortion rates corresponding to the various typical weather scenes are substituted into the following formula for calculation and analysis to obtain the harmonic voltage distortion rate THD of the distributed photovoltaic access power distribution network in a certain region,
Figure GDA0001672217950000031
in the formula: piThe probabilities corresponding to various typical weather scenes; THDUiFor the harmonic voltage distortion rate in various types of typical weather scenes, i is a classification label of various types of typical weather scenes, and S is the number of the typical weather scenes.
According to the method, firstly, the output power curve of the distributed photovoltaic system in a typical weather scene is used for obtaining the Beta distribution shape parameters through a parameter estimation method, a probability distribution model of the output power of the photovoltaic system is determined, sampling is carried out, the distributed photovoltaic output power value is obtained, harmonic load flow calculation is carried out, the harmonic injection current of the power distribution network is obtained, and further the harmonic voltage of each node of the system is obtained. And after all sampling times are finished, calculating to obtain the harmonic voltage probability distribution and the harmonic distortion rate THD of the power distribution network.
Drawings
Fig. 1 is a schematic diagram of photovoltaic power generation.
Fig. 2 is an equivalent circuit of a single harmonic distribution network.
FIG. 3 is a flow chart of the calculation of the present invention.
Fig. 4 shows a medium voltage radiation type distribution network system topology of IEEE33 nodes.
Fig. 5 is a graph of branch circuit parameters and bus load data for an IEEE33 system in an embodiment of the present invention.
FIG. 6 is a graph of the harmonic voltage content of a PCC point for a first type of weather in accordance with an embodiment of the present invention.
FIG. 7 is a graph of the harmonic voltage content of a PCC point on a second typical weather in accordance with an embodiment of the present invention.
FIG. 8 is a graph of the voltage content of each harmonic of a PCC point on a third type of typical weather in accordance with an embodiment of the present invention.
FIG. 9 is a graph of the harmonic voltage content of a PCC point for a fourth type of typical weather in accordance with an embodiment of the present invention.
Fig. 10 is a fifth harmonic voltage probability distribution of PCC during a first typical weather in accordance with an embodiment of the present invention.
Fig. 11 is a seventh harmonic voltage probability distribution of a PCC under a first typical weather in an embodiment of the present invention.
Fig. 12 is a probability distribution of eleven harmonic voltages at the PCC points under a first typical weather condition in accordance with an embodiment of the present invention.
Fig. 13 is a thirteen-harmonic voltage probability distribution of the PCC points under the first typical weather in the embodiment of the present invention.
Fig. 14 is a seventy-seventh harmonic voltage probability distribution of a PCC under a first typical weather in an embodiment of the present invention.
Fig. 15 is a fifth harmonic voltage probability distribution of PCC during a second typical weather in accordance with an embodiment of the present invention.
Fig. 16 is a seventh harmonic voltage probability distribution of a PCC under a second typical weather in an embodiment of the present invention.
Fig. 17 is a probability distribution of eleven harmonic voltages at a PCC in a second typical weather in accordance with an embodiment of the present invention.
Fig. 18 is a thirteen-harmonic voltage probability distribution of a PCC under a second typical weather condition in accordance with an embodiment of the present invention.
Fig. 19 is a seventy-seventh harmonic voltage probability distribution of a PCC under a second typical weather in an embodiment of the present invention.
Fig. 20 is a fifth harmonic voltage probability distribution of PCC during a third typical weather in accordance with an embodiment of the present invention.
Fig. 21 is a seventh harmonic voltage probability distribution of a PCC during a third typical weather in accordance with an embodiment of the present invention.
Fig. 22 is a probability distribution of eleven harmonic voltages at the PCC points under a third typical weather in an embodiment of the present invention.
Fig. 23 is a thirteen-harmonic voltage probability distribution of a PCC under a third typical weather in an embodiment of the present invention.
Fig. 24 is a seventy-seventh harmonic voltage probability distribution of a PCC under a third typical weather in an embodiment of the present invention.
Fig. 25 is a fifth harmonic voltage probability distribution of PCC during a fourth typical weather in accordance with an embodiment of the present invention.
Fig. 26 is a seventh harmonic voltage probability distribution of a PCC under a fourth typical weather in an embodiment of the present invention.
Fig. 27 is a probability distribution of eleven harmonic voltages at the PCC in a fourth exemplary weather in accordance with an embodiment of the present invention.
Fig. 28 is a thirteen-harmonic voltage probability distribution of a PCC under a fourth typical weather in an embodiment of the present invention.
Fig. 29 is a seventy-seventh harmonic voltage probability distribution of a PCC under a fourth typical weather in an embodiment of the present invention.
Detailed Description
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to these drawings without creative efforts.
Referring to fig. 1 to 3, an embodiment of the present invention provides a probability harmonic analysis method for a distributed photovoltaic access distribution network considering a weather scene.
1. Output power probability model for constructing distributed photovoltaic power generation
(1) Photovoltaic power generation power output model
As shown in fig. 1, the principle of solar photovoltaic power generation mainly uses the photoelectric effect to convert the solar light energy into electric energy, and its main components include a photovoltaic array formed by solar panels, and an inverter and a controller used for grid connection. The power output of a photovoltaic cell is mainly related to the intensity of solar radiation, the area of the photovoltaic array and the photoelectric conversion efficiency. The output power of photovoltaic power generation is therefore approximated by
PPV=r·A·ρ (1)
In the formula, r is solar irradiance, A is the area of the photovoltaic array, and rho is photoelectric conversion efficiency.
The output power of the photovoltaic power generation can be obtained by inquiring or counting according to historical monitoring data of photovoltaic power stations in certain areas.
(2) Probability distribution model for solar illumination radiation intensity under different weather scenes
The solar irradiance is generally most commonly distributed by adopting a Beta distribution, and the probability density function of the Beta distribution is as follows:
Figure GDA0001672217950000061
wherein r is solar irradiance, rmaxAlpha and Beta are the shape parameters of the Beta distribution for the solar irradiance maximum.
In the same region, there are many factors affecting the intensity of illumination radiation, wherein the weather scene is one of the main factors. The shape parameters of the probability distribution functions of the corresponding illumination radiation intensities in different weather scenes are also different.
(3) Probability distribution model of photovoltaic power generation power
The probability density function of the photovoltaic power generation power output can be obtained by combining the existing photovoltaic power generation power output model and the probability distribution function of the illumination intensity
Figure GDA0001672217950000062
In the formula, PpvFor photovoltaic power generation output power, PmaxAlpha and Beta are shape parameters of Beta distribution function, and are the maximum value of the output power of the photovoltaic power generation; pmaxCan be calculated by the following formula:
Pmax=rmax·A·ρ (4)
2. clustering analysis and parameter estimation of typical weather scenes
(1) Typical weather scene clustering analysis
In order to make the analyzed weather condition and photovoltaic output more representative, the invention obtains and obtains a photovoltaic output curve under a typical weather scene by adopting a cluster analysis method.
Clustering is the division of a collection of data units into several subsets called clusters, where the data in each cluster has greater similarity. Clustering is mainly used for finding valuable data distribution and data patterns in potential data, and is an important component in data mining. The clustering process can be defined as follows: given N data points of d-dimensional space, the N data points are divided into k clusters according to the degree of similarity between the data points, i.e. it is satisfied that similar samples are in the same cluster and dissimilar samples are in different clusters, so that objects in the same cluster have as much similarity as possible, while objects in different clusters have as much dissimilarity as possible. If a data set X containing N samples X1, …, xN is clustered into c sub-classes X1, …, Xc, X1, …, Xc is required to satisfy:
Figure GDA0001672217950000071
classical cluster analysis methods include the following categories: (1) the marking method typically represents algorithms such as K-MEANS, K-MEDOIDS and CLARANS algorithms. (2) The hierarchy method typically represents algorithms such as CURE, Chameleon, ROCK and BRICH algorithms. (3) Typical algorithms for the density-based clustering method include a density algorithm, a DBSCAN algorithm, an OPTIC algorithm, and the like. (4) A grid-based clustering method. Common grid-based clustering algorithms include a statistical information grid method STING, a clustering high-dimensional space method CIQUE, a wavelet transform-based clustering method Wave Cluster and the like. (5) A model-based clustering method. The method mainly comprises a statistical method and a neural network method. (6) A group intelligence algorithm. Typically representing the main ant colony algorithm and the particle swarm optimization algorithm. The invention adopts the most classical K-MEANS clustering method for analysis.
(2) Parameter estimation
The most widely used method in parameter estimation is the least squares method. The invention adopts a least square method to carry out parameter estimation of two parameters.
For given data (x)i,yi) I is 1, …, n, is
Figure GDA0001672217950000072
Is linearly independent, find
Figure GDA0001672217950000081
Make an error
Figure GDA00016722179500000811
Has the smallest sum of squares, i.e.
Figure GDA0001672217950000082
In a geometric sense, it is sought to associate a given point (x)i,yi) The curve y (p) (x) in which the sum of squared distances of i (1) and … (n) is the minimum, the function p (x) is called a fitting function or a least-squares solution, and the method of obtaining the fitting function p (x) is called a least-squares method of curve fitting.
Is provided with
Figure GDA0001672217950000083
When in use
Figure GDA0001672217950000084
Time, error ri=p(xi)-yiI-1, …, m is the smallest sum of squares, if defined
Figure GDA0001672217950000085
Figure GDA0001672217950000086
Then satisfy
Figure GDA0001672217950000087
The system of linear equations of (1) is:
Figure GDA0001672217950000088
due to the fact that
Figure GDA0001672217950000089
The linear independence, the determinant of the coefficient matrix is not equal to 0, and the equation set has a unique solution.
When the least square method is used for parameter estimation, only a given set of data is required to be substituted into the equation set
Figure GDA00016722179500000810
Obtain the unknown aiThe unknown a can be obtained by left divisioni
3. Principle of harmonic analysis
(1) Harmonic analysis of distributed photovoltaic power distribution network
When the photovoltaic inverter converts direct current into alternating current, certain harmonic waves are generated. The content of the output harmonic is related to the topological structure, the control strategy and the parameter setting of the inverter. When the photovoltaic output is changed, the content of harmonic current generated by photovoltaic power generation is also changed.
The harmonic model of the power distribution network comprises all power elements, the harmonic model comprises distributed photovoltaic, a transformer, a power transmission line and a load, the harmonic model of the power distribution network containing the distributed photovoltaic is built for researching the harmonic current distribution condition of an ideal power distribution network in a steady state, and an equivalent circuit diagram of the single harmonic power distribution network is shown in fig. 2.
In the power distribution network researched by the invention, except for the distributed photovoltaic inverter, other harmonic sources do not exist, namely, no harmonic exists in the power distribution network. After distributed photovoltaic access, h-order harmonic current I is injected into the power distribution network as a harmonic current sourcesh
Figure GDA0001672217950000098
H-harmonic equivalent impedance, I, of a load behind a photovoltaic access Point (PC)LbhFor h-harmonic currents of the load after the point PC, ILahFor h harmonic currents before flowing to the access point,
Figure GDA0001672217950000099
h-harmonic equivalent impedance of the load before the PC point, ILhFor the h harmonic current flowing into the load before the PC point,
Figure GDA00016722179500000910
is h-order harmonic equivalent impedance between a Load access Point (PL) and a Common access Point (PCC), which mainly depends on the harmonic impedance of a transformer and a line, IeqhIs the h harmonic current flowing into the PCC point.
(2) Distributed photovoltaic injection harmonic current distribution analysis
H-order harmonic current is injected into a power distribution network from a PC point of distributed photovoltaic system
Ish-ILbh-ILah=0 (7)
Figure GDA0001672217950000091
Harmonic current I flowing to PCC pointeqhIs composed of
Figure GDA0001672217950000092
As can be seen from the formula (9),harmonic current I flowing into PCC pointeqhBy
Figure GDA0001672217950000093
And
Figure GDA0001672217950000094
are determined jointly by
Figure GDA0001672217950000095
Much less than
Figure GDA0001672217950000096
And
Figure GDA0001672217950000097
the harmonic currents of the distributed photovoltaic injection distribution network mostly flow into the PCC points.
The most common indicator of voltage waveform is Total Harmonic Distortion (THD), which takes into account the effect of each Harmonic component on the signal. The calculation formula is as follows:
Figure GDA0001672217950000101
in the formula: uh: the harmonic voltage amplitude of the h harmonic; u shape1: a fundamental voltage amplitude; n: the highest order of the harmonic considered;
4. probability harmonic analysis method framework of distributed photovoltaic access considering typical weather scene
Sampling photovoltaic power generation power in the photovoltaic power generation power output probability density function under various typical weather scenes by adopting a Monte Carlo method, and performing harmonic load flow calculation on the sampling value of the photovoltaic power generation power to obtain a harmonic current injection value and a harmonic voltage of each subharmonic injected into the power distribution network;
if the sampling times are less than the sampling set times, returning to the previous step to continue sampling calculation, and otherwise, performing the next step; substituting the harmonic voltage of each harmonic into the following harmonic voltage distortion rate formula to calculate to obtain each harmonic voltage under each typical weather sceneHarmonic voltage distortion rate THD of harmonicu
And a sixth step: after the harmonic voltage distortion rate analysis and calculation under various typical weather scenes are completed, the probabilities and the harmonic voltage distortion rates corresponding to the various typical weather scenes are substituted into the following formula for calculation and analysis to obtain the harmonic voltage distortion rate THD of the distributed photovoltaic access power distribution network in a certain region,
Figure GDA0001672217950000102
in the formula: piThe probabilities corresponding to various typical weather scenes; THDUiAnd i is a classification label of each type of typical weather scene.
[ EXAMPLES ]
An IEEE33 node medium voltage radiation type power distribution network system is adopted, and the topological structure is shown in figure 4. The three-phase power reference value Sb of the power distribution system is 10000kVA, the line voltage reference value UB is 12.66kV, the total active load of the system is 3715.0kW, and the total reactive load is 2300.0 kvar. The load data of each node is shown in the table. In the example, node 0 is considered as the balancing node, and the rest of the load nodes are considered as the PQ nodes. And (3) selecting a node 18 to join the distributed photovoltaic system, wherein the installed capacity is 1MW, and the permeability of the distributed photovoltaic is 23%. The branch parameters and bus load data for the system are shown in figure 5.
(1) The photovoltaic output power curve of a certain area all year round is analyzed, and the results of the alpha and beta parameters under four typical weather scenes are obtained by adopting the least square method as shown in the following table 1.
Table 1:
Figure GDA0001672217950000111
(2) monte Carlo simulation calculations show the voltage content of each harmonic of the PCC points for four typical days as shown in FIGS. 6-9.
(3) The results of the probability distribution of the respective harmonic voltages at the PCC points under typical weather are shown in fig. 10-29.
(4) The harmonic distortion rate calculation results are shown in table 2 below.
Table 2:
Figure GDA0001672217950000112
Figure GDA0001672217950000113
the modules or units in the device of the embodiment of the invention can be combined, divided and deleted according to actual needs.
While the invention has been described with reference to a preferred embodiment, it will be understood by those skilled in the art that various changes in form and detail may be made therein without departing from the spirit and scope of the invention.

Claims (1)

1.一种考虑天气场景的分布式光伏接入配电网的概率谐波分析方法,其特征在于包括以下步骤:1. a probabilistic harmonic analysis method for a distributed photovoltaic access distribution network considering weather scenarios, is characterized in that comprising the following steps: 第一步:利用K-MEANS聚类分析方法,对分布式光伏接入配电网的光伏发电输出功率Ppv进行统计和聚类,得到典型天气场景的光伏输出功率曲线及其分布概率;The first step: using the K-MEANS cluster analysis method, statistics and cluster the photovoltaic power generation output power P pv of the distributed photovoltaic connected to the distribution network, and obtain the photovoltaic output power curve and its distribution probability of typical weather scenarios; 第二步:根据第一步典型天气场景的光伏输出功率曲线,采用最小二乘法进行参数估计,建立光伏电站在不同天气场景下的光伏发电功率输出概率密度函数,即Beta分布函数:The second step: According to the photovoltaic output power curve of the typical weather scenario in the first step, the least squares method is used to estimate the parameters, and the output probability density function of the photovoltaic power generation power of the photovoltaic power station in different weather scenarios is established, that is, the Beta distribution function:
Figure FDA0002776920300000011
Figure FDA0002776920300000011
式中,PPV为光伏发电输出功率,Pmax为光伏发电输出功率的最大值,α和β为Beta分布函数的形状参数;In the formula, P PV is the output power of photovoltaic power generation, P max is the maximum value of output power of photovoltaic power generation, and α and β are the shape parameters of the Beta distribution function; 第三步:基于光伏发电功率输出概率密度函数,采用蒙特卡洛法分别对上述各类典型天气场景下光伏发电功率进行采样,采样得到该天气场景下的光伏输出功率采样值;Step 3: Based on the output probability density function of photovoltaic power generation, the Monte Carlo method is used to sample the photovoltaic power generation power in the above various typical weather scenarios, and the sampling value of the photovoltaic output power under the weather scene is obtained by sampling; 第四步:根据光伏发电功率的采样值,对配电网光伏输出功率采样值进行谐波潮流计算,得到注入配电网的各次谐波电流和谐波电压值;Step 4: According to the sampling value of photovoltaic power generation power, perform harmonic power flow calculation on the sampling value of photovoltaic output power of the distribution network, and obtain the harmonic current and harmonic voltage values of each order injected into the distribution network; 第五步:若采样次数小于采样设定次数,返回第三步继续采样计算,计算注入配电网的各次谐波电压均值,否则进行下一步;Step 5: If the sampling times are less than the sampling set times, return to the third step to continue the sampling calculation, and calculate the average value of each harmonic voltage injected into the distribution network, otherwise go to the next step; 第六步:将所述各次谐波的谐波电压均值代入以下谐波电压畸变率公式计算得到各类典型天气场景下各次谐波的谐波电压畸变率THDuStep 6: Substitute the harmonic voltage mean value of each harmonic into the following harmonic voltage distortion rate formula to calculate the harmonic voltage distortion rate THD u of each harmonic under various typical weather scenarios; 谐波电压畸变率公式:Harmonic voltage distortion rate formula:
Figure FDA0002776920300000021
Figure FDA0002776920300000021
式中:Uh:h次谐波的谐波电压幅值;U1:基波电压幅值;n:所考虑的谐波的最高次数;In the formula: Uh: the harmonic voltage amplitude of the h harmonic; U 1 : the fundamental voltage amplitude; n: the highest order of the considered harmonic; 第七步:完成各类典型天气场景下的谐波电压畸变率分析计算后,将所述各类典型天气场景对应的概率及谐波电压畸变率代入下式计算分析,得到某地区的分布式光伏接入配电网的谐波电压畸变率THD,Step 7: After completing the analysis and calculation of the harmonic voltage distortion rate under various typical weather scenarios, substitute the probability and harmonic voltage distortion rate corresponding to the various typical weather scenarios into the following formula for calculation and analysis, and obtain the distributed distribution in a certain area. The harmonic voltage distortion rate THD of photovoltaic connected to the distribution network,
Figure FDA0002776920300000022
Figure FDA0002776920300000022
式中:Pi为各类典型天气场景对应的概率;THDUi为各类典型天气场景下的谐波电压畸变率,i为各类典型天气场景的分类标号,S为典型天气场景数量。where P i is the probability corresponding to various typical weather scenarios; THD Ui is the harmonic voltage distortion rate under various typical weather scenarios, i is the classification label of various typical weather scenarios, and S is the number of typical weather scenarios.
CN201711207394.8A 2017-11-27 2017-11-27 A Probabilistic Harmonic Analysis Method for Distributed Photovoltaic Access to Distribution Network Considering Weather Scenarios Active CN108233373B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711207394.8A CN108233373B (en) 2017-11-27 2017-11-27 A Probabilistic Harmonic Analysis Method for Distributed Photovoltaic Access to Distribution Network Considering Weather Scenarios

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711207394.8A CN108233373B (en) 2017-11-27 2017-11-27 A Probabilistic Harmonic Analysis Method for Distributed Photovoltaic Access to Distribution Network Considering Weather Scenarios

Publications (2)

Publication Number Publication Date
CN108233373A CN108233373A (en) 2018-06-29
CN108233373B true CN108233373B (en) 2021-01-08

Family

ID=62652923

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711207394.8A Active CN108233373B (en) 2017-11-27 2017-11-27 A Probabilistic Harmonic Analysis Method for Distributed Photovoltaic Access to Distribution Network Considering Weather Scenarios

Country Status (1)

Country Link
CN (1) CN108233373B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110957725A (en) * 2019-10-29 2020-04-03 国网江苏省电力有限公司盐城供电分公司 Photovoltaic inverter based on harmonic estimation and control method

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2010295163A1 (en) * 2009-09-18 2012-04-12 Sparq Systems Inc. Distributed power generation interface
CN106548248A (en) * 2016-10-31 2017-03-29 国家电网公司 A kind of harmonic trend probability index Forecasting Methodology of the electrical network of power station containing large-sized photovoltaic
CN107046287A (en) * 2017-04-06 2017-08-15 华北电力大学 A kind of large-sized photovoltaic power station harmonic wave Over-voltage Analysis method

Also Published As

Publication number Publication date
CN108233373A (en) 2018-06-29

Similar Documents

Publication Publication Date Title
Basaran et al. Energy management for on‐grid and off‐grid wind/PV and battery hybrid systems
CN107147152B (en) Collaborative optimization configuration method and system for multi-type active and reactive power sources in new energy distribution network
CN108306303A (en) Voltage stability evaluation method considering load increase and new energy output randomness
CN107301472A (en) Distributed photovoltaic planing method based on scene analysis method and voltage-regulation strategy
CN107910863A (en) Consider the power distribution network dispatching method that photovoltaic is contributed with workload demand forecast interval
CN109274136A (en) A Reactive Power Optimization Method for Photovoltaic System Based on Quantum Particle Swarm Optimization
CN106527568A (en) Maximum power tracker control method based on generalized dynamic fuzzy neural network
CN110909911A (en) Aggregation method for multidimensional time series data considering spatiotemporal correlation
CN105305488B (en) A kind of evaluation method for considering new-energy grid-connected and power transmission network utilization rate being influenced
CN109034587A (en) A kind of active distribution system Optimization Scheduling for coordinating a variety of controllables
Zhu et al. An Energy Storage Performance Improvement Model for Grid‐Connected Wind‐Solar Hybrid Energy Storage System
CN114996908A (en) Active power distribution network extension planning method and system considering intelligent soft switch access
Ali et al. Optimal planning of uncertain renewable energy sources in unbalanced distribution systems by a multi‐objective hybrid PSO–SCO algorithm
Ni et al. A review of line loss analysis of the low-voltage distribution system
Zhao et al. Research on Multiobjective Optimal Operation Strategy for Wind‐Photovoltaic‐Hydro Complementary Power System
CN108233373B (en) A Probabilistic Harmonic Analysis Method for Distributed Photovoltaic Access to Distribution Network Considering Weather Scenarios
CN109149658A (en) Independent micro-capacitance sensor distributed dynamic economic load dispatching method based on congruity theory
Liu et al. Reliability assessment of micro-grid based on Markov model
Yin Fast frequency response technology of photovoltaic power plant based on equivalent modelling
CN103019212A (en) Multi-photovoltaic solar-diesel hybrid system coordination control method based on disturbance observer
CN117748444A (en) Operation simulation method of power distribution system
CN114329857B (en) A distributed power generation planning method based on improved whale algorithm
Xu et al. One-day ahead forecast of PV output based on deep belief network and weather classification
CN114243766B (en) Regional multi-energy system optimal configuration method and system
CN115271445A (en) IIES risk avoidance optimization method and system considering renewable energy uncertainty

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