CN109599899A - A kind of setting method of new energy running simulation boundary condition - Google Patents

A kind of setting method of new energy running simulation boundary condition Download PDF

Info

Publication number
CN109599899A
CN109599899A CN201811531549.8A CN201811531549A CN109599899A CN 109599899 A CN109599899 A CN 109599899A CN 201811531549 A CN201811531549 A CN 201811531549A CN 109599899 A CN109599899 A CN 109599899A
Authority
CN
China
Prior art keywords
wind
wind speed
area
parameter
photovoltaic
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
CN201811531549.8A
Other languages
Chinese (zh)
Other versions
CN109599899B (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
Economic and Technological Research Institute of State Grid Shandong Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Economic and Technological Research Institute of State Grid Shandong 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 State Grid Corp of China SGCC, Economic and Technological Research Institute of State Grid Shandong Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201811531549.8A priority Critical patent/CN109599899B/en
Publication of CN109599899A publication Critical patent/CN109599899A/en
Application granted granted Critical
Publication of CN109599899B publication Critical patent/CN109599899B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H02J3/383
    • H02J3/386
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • 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
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects

Abstract

This application discloses a kind of setting methods of new energy running simulation boundary condition, wherein wind-powered electricity generation running simulation boundary condition setting method includes: according to geographical location, all wind power plants in region to be measured are divided into multiple wind areas, wind farm wind velocity statistical distribution parameter is set according to historical wind speed data, the wind speed related coefficient between multiple wind areas is set according to the distance between historical wind speed data and wind area, blower power producing characteristics parameter is set according to blower model.The setting of photovoltaic running boundary condition includes: setting photovoltaic plant parameter, and setting sets clearness index distribution parameter and according to the distance between photovoltaic plant section, sets the related coefficient between photovoltaic plant section.By the application, the reliability of boundary condition setting can be improved, to improve the accuracy for influencing running simulation result.

Description

A kind of setting method of new energy running simulation boundary condition
Technical field
This application involves new energy grid-connected power technical fields, more particularly to a kind of new energy running simulation boundary condition Setting method.
Background technique
With the development of new energy technology and technology of generating electricity by way of merging two or more grid systems, how new energy efficiently, to be reasonably applied to grid-connected Power generation, is the key problem in terms of new energy grid-connected power.Due to new energy power output there are certain randomness, fluctuation with not Certainty, extensive new energy power generation grid-connection can to the power supply of electric system contribute structure, operation form, the energy saving of operation with Economy etc., which all has, to be significantly affected, to influence the planning of entire electric system.
In order to predict influence of the extensive new energy grid-connected power to electric system, is generallyd use in industry and establish mathematical modulo The method of type simulates influence of the new energy to network system, estimates the effect of new energy grid-connected power according to analog result.And it builds When formwork erection type, need to set new energy running simulation boundary condition, and the boundary condition of running simulation determines the knot of running simulation Therefore fruit sets its simulating boundary condition for data model, be a major issue.
Currently, to the method that the new energy simulating boundary condition in mathematical model is set, usually according to current mould Then parameter required for type determines determines simulating boundary condition according to canonical parameter, and canonical parameter is to come based on experience value It determines.
However, in the method set at present to new energy simulating boundary condition, only due to selected canonical parameter It determines based on experience value, it is not intended that correlation between the corresponding practical new energy carrier of model, therefore, using mesh Parameter set by preceding simulating boundary condition enactment method is not accurate enough, thus influence running simulation as a result, causing in turn It is not accurate enough to the assessment result of extensive new energy power generation grid-connection.
Summary of the invention
This application provides a kind of setting methods of new energy running simulation boundary condition, set in the prior art to solve The not accurate enough problem of fixed parameter.
In order to solve the above-mentioned technical problem, the embodiment of the present application discloses following technical solution:
A kind of setting method of new energy running simulation boundary condition, which comprises wind-powered electricity generation running simulation perimeter strip Part setting and the setting of photovoltaic running simulation boundary condition;
Wherein, the method for the wind-powered electricity generation running simulation boundary condition setting, comprising:
According to geographical location, all wind power plants in region to be measured are divided into multiple wind areas, include multiple in each wind area Wind power plant;
According to historical wind speed data, wind farm wind velocity statistical distribution parameter is set, the historical wind speed data is current wind The average value of wind power plant historical wind speed data contained by area;
According to the distance between historical wind speed data and wind area, the wind speed related coefficient between multiple wind areas is set;
According to blower model, blower power producing characteristics parameter is set;
The method of the photovoltaic running simulation boundary condition setting, comprising:
Set photovoltaic plant parameter, the photovoltaic plant parameter be used for calculate it is unobstructed under the conditions of theoretical solar irradiation it is strong Degree and photovoltaic plant power output;
Clearness index distribution parameter is set, the clearness index distribution parameter is for reflecting clear sky relevant to photovoltaic power generation Index;
According to the distance between photovoltaic plant section, the related coefficient between photovoltaic plant section is set.
Optionally, described according to historical wind speed data, set wind farm wind velocity statistical distribution parameter, comprising:
Calculate the fitting wind speed statistical distribution parameter of each wind power plant in any wind area in multiple wind areas;
According to the capacity of each wind power plant, the fitting wind speed statistical distribution parameter is taken into weighted average, is calculated Obtain the wind speed statistical distribution parameter in any wind area.
Optionally, the fitting wind speed statistical distribution parameter for calculating each wind power plant in any wind area in multiple wind areas, Include:
According to the historical wind speed data of wind power plant, using Weibull (Wei Buer) distribution function fitting wind farm wind velocity point Scale parameter, form parameter and the auto-correlation attenuation coefficient of cloth;
According to the historical wind speed data of wind power plant, calculate wind farm wind velocity the moon characteristic mark curve and day characteristic mark it is bent Line.
Optionally, described according to historical wind speed data, set wind farm wind velocity statistical distribution parameter, further includes:
When there is no historical wind speed data in any wind area, using the wind speed statistical distribution parameter in adjacent wind area as institute The wind speed statistical distribution parameter in any wind area is stated, the adjacent sectors are the wind nearest on geographical location with any wind area Area.
Optionally, described according to the distance between historical wind speed data and wind area, set the wind speed between multiple wind areas Related coefficient, comprising:
Judge whether current wind area has historical wind speed data;
If so, being fitted the wind speed related coefficient between any two wind area according to the historical wind speed data in current wind area;
If not, being utilized according to the distance between any two wind areaIt is calculated between any two wind area Wind speed related coefficient, wherein c is wind speed related coefficient, and d is the distance between any two wind area, and M is wind speed related coefficient With the range attenuation factor.
Optionally, the photovoltaic plant parameter includes: geographical location, photovoltaic array type, photovoltaic array tilt angle, light Photovoltaic array deflection, installed capacity, efficiency and prediction absolute error account for installed capacity percentage, the geographical location, photovoltaic Array type, photovoltaic array tilt angle and photovoltaic array deflection, for calculate it is unobstructed under the conditions of theoretical solar irradiation Intensity, the installed capacity, efficiency and prediction absolute error account for installed capacity percentage, for calculating photovoltaic plant power output.
Optionally, the clearness index distribution parameter, comprising: clearness index probability distribution parameters, clearness index probability point Cloth parameter lamda value, clearness index probability distribution parameters theta value, atmospheric scattering coefficient p, atmospheric scattering coefficient q value and bottom Face reflectivity.
It is optionally, described that related coefficient between photovoltaic plant section is set according to the distance between photovoltaic plant section, Include:
Whether judge the distance between any two photovoltaic plant section >=200km;
If so, the related coefficient between any two photovoltaic plant section is set as 0.15;
If not, the related coefficient between any two photovoltaic plant section is set as 0.4.
The technical solution that embodiments herein provides can include the following benefits:
The application provides a kind of setting method of new energy running simulation boundary condition, and this method mainly includes wind-powered electricity generation operation The setting of simulating boundary condition and the setting of photovoltaic running boundary condition.Wherein, in your boundary condition setting method of wind-powered electricity generation running simulation, Wind power plant is divided into multiple wind areas first, in accordance with geographical location, then sets wind farm wind velocity statistical distribution parameter, wind speed phase Relationship number and blower power producing characteristics parameter.In photovoltaic running boundary condition enactment method, main includes setting photovoltaic plant ginseng Related coefficient between number, clearness index distribution parameter and photovoltaic plant section.To wind-powered electricity generation running simulation side in the present embodiment When boundary's condition is set, the geographical location relationship between current wind area and other wind areas is fully considered, according to different geographical positions The relationship of setting determines related coefficient, and acquired setting result is more accurate, and joins for the wind speed statistical distribution in each wind area Number is not simple averaged, it is also contemplated that the capacity of wind power plant, to all wind power plants in wind area using Function Fitting Mode obtains entire wind area wind speed statistical distribution parameter, is conducive to the accuracy for improving boundary condition setting result.The present embodiment When setting to photovoltaic power generation running simulation boundary condition, the photovoltaics such as solar irradiation intensity, photovoltaic plant power output are fully considered Power station parameter and clearness index distribution parameter are covered the various parameters for influencing photovoltaic power generation more fully hereinafter, and are sufficiently examined Consider the geographical location relationship between photovoltaic plant section, can between photovoltaic plant itself and adjacent photovoltaic power station mutually The angle initialization photovoltaic running simulation boundary condition of influence is conducive to the reliability for improving boundary condition setting, to improve shadow Ring the accuracy of running simulation result.
It should be understood that above general description and following detailed description be only it is exemplary and explanatory, not The application can be limited.
Detailed description of the invention
The drawings herein are incorporated into the specification and forms part of this specification, and shows the implementation for meeting the application Example, and together with specification it is used to explain the principle of the application.
In order to illustrate the technical solutions in the embodiments of the present application or in the prior art more clearly, to embodiment or will show below There is attached drawing needed in technical description to be briefly described, it should be apparent that, for those of ordinary skill in the art Speech, without creative efforts, is also possible to obtain other drawings based on these drawings.
Fig. 1 shows for a kind of process of the setting method of new energy running simulation boundary condition provided by the embodiment of the present application It is intended to.
Fig. 2 is the distribution of correlation coefficient schematic diagram between the photovoltaic power output occlusion coefficient of each observation point in the embodiment of the present application.
Specific embodiment
In order to make those skilled in the art better understand the technical solutions in the application, below in conjunction with the application reality The attached drawing in example is applied, the technical scheme in the embodiment of the application is clearly and completely described, it is clear that described implementation Example is merely a part but not all of the embodiments of the present application.Based on the embodiment in the application, this field is common The application protection all should belong in technical staff's every other embodiment obtained without making creative work Range.
The application in order to better understand explains in detail presently filed embodiment with reference to the accompanying drawing.
Referring to Fig. 1, Fig. 1 is a kind of setting method of new energy running simulation boundary condition provided by the embodiment of the present application Flow diagram.As shown in Figure 1, the setting method of new energy running simulation boundary condition specifically includes that wind in the present embodiment The method of electricity operation simulating boundary condition setting and method two parts of photovoltaic running simulation boundary condition setting.
Wherein, the method for wind-powered electricity generation running simulation boundary condition setting specifically includes that
S11: according to geographical location, all wind power plants in region to be measured is divided into multiple wind areas, include in each wind area Multiple wind power plants.
S12: according to historical wind speed data, wind farm wind velocity statistical distribution parameter is set.Wherein, historical wind speed data is to work as The average value of wind power plant historical wind speed data contained by the area Qian Feng.
Wind power plant historical wind speed data can be back-calculated to obtain according to history wind power output data.
Specifically, step S12 is comprised the following processes:
S121: the fitting wind speed statistical distribution parameter of each wind power plant in any wind area in multiple wind areas is calculated.
Step S121 mainly includes the following steps:
S1211: according to the historical wind speed data of wind power plant, using the fitting wind farm wind velocity distribution of Weibull distribution function Scale parameter, form parameter and auto-correlation attenuation coefficient.
S1212: according to the historical wind speed data of wind power plant, the moon characteristic mark curve and day characteristic of wind farm wind velocity is calculated Mark curve.
Weibull distribution function is a kind of unimodal, two parameter distribution function cluster, is divided in the present embodiment using Weibull Cloth function, probability density function areWherein, k is form parameter, and c is scale parameter, the present embodiment By determining scale parameter c and form parameter k, the probability distribution of wind speed in 1 year can be calculated.In the present embodiment preferably Two-parameter Weibull distribution function can quickly and more accurately obtain wind-powered electricity generation using Two-parameter Weibull distribution function Field wind speed boundary condition, to improve the efficiency and accuracy of wind-powered electricity generation running simulation boundary condition setting.
S122: according to the capacity of each wind power plant, fitting wind speed statistical distribution parameter is taken into weighted average, is calculated The wind speed statistical distribution parameter in any wind area.
The present embodiment can using average weighted method to the setting of the wind speed statistical distribution parameter in current wind area It fully considers the influence of the capacity and weight of each wind power plant in current wind area to entire wind area distribution parameter, can obtain more Add accurate boundary parameter, to improve the accuracy and reliability of dry run result.Determine the wind speed statistics in current wind area After distribution parameter, following all wind power plants are all made of the simulation air speed data within the scope of same wind area, it may be assumed that set same wind area The wind speed of interior wind power plant is perfectly correlated.
Further, in this embodiment according to historical wind speed data, the method for setting wind farm wind velocity statistical distribution parameter, Further include step S123: when there is no historical wind speed data in any wind area, being made with the wind speed statistical distribution parameter in adjacent wind area For the wind speed statistical distribution parameter in any wind area.Wherein, adjacent sectors are the wind areas nearest on geographical location with any wind area.
With continued reference to Fig. 1 it is found that the method for the present embodiment apoplexy electricity operation simulating boundary condition setting, further includes step S13: according to the distance between historical wind speed data and wind area, the wind speed related coefficient between multiple wind areas is set.
Specifically, step S13 is comprised the following processes:
S131: judge whether current wind area has historical wind speed data.
Specifically, judge whether there is at least one wind power plant history to go out force data or historical wind speed in current wind area Data.
If having historical wind speed data in current wind area, step S132 is executed: according to the historical wind speed data in current wind area The wind speed related coefficient being fitted between any two wind area.
If there is historical wind speed data in current wind area, step S133 is executed: according to the distance between any two wind area, It utilizesThe wind speed related coefficient between any two wind area is calculated.Wherein, c is wind speed related coefficient, and d is to appoint The distance between area Yi Lianggefeng, M are wind speed related coefficient with the range attenuation factor.It is, all wind-powered electricity generations in current wind area When field all goes out force data or historical wind speed data without history, step S133 is executed.According to being fitted between resulting wind area Related coefficient, it is estimated that wantonly one or two of wind can be calculated then in conjunction with the distance of each wind-powered electricity generation section in the value of M Wind speed related coefficient between area.
Wind-powered electricity generation section part wind speed correlation is mainly related with geographic distance, at a distance of closer wind-powered electricity generation section due to by same The influence of weather conditions, wind speed will show stronger correlation;And wind-powered electricity generation section apart from each other, it is identical due to encountering The probability of weather conditions is smaller, therefore wind speed correlation is weaker, and the present embodiment utilizes between wind speed related coefficient and wind-powered electricity generation section Negative exponent relationship, enable to set boundary condition more accurate.
S14: according to blower model, blower power producing characteristics parameter is set.
The power producing characteristics parameter of blower is mainly determined that the technical documentation that can specifically consult current blower is come by blower model It determines.
Since the treatment characteristic parameter difference of different types of blower is little, representative value as shown in Table 1 can also be taken to make For blower power producing characteristics parameter.
The typical wind electric field blower parameter of table 1
Below by taking the planning of Shandong Power wind power plant as an example, the method for description wind-powered electricity generation running simulation boundary condition setting:
Firstly, wind power plant is divided into 17 areas Ge Feng according to wind power plant affiliated prefecture-level city, it is respectively as follows: the area Ji Nanfeng, Qingdao wind Area, the area Zi Bofeng, the area Zao Zhuanfeng, Dongying Feng Qu, cigarette typhoon zone, the area Wei Fangfeng, the area Ji Ningfeng, the area Tai Anfeng, the area Wei Haifeng, day According to wind area, the area Bin Zhoufeng, the area De Zhoufeng, the area Liao Chengfeng, the area Lin Yifeng, the area He Zefeng and Laiwu wind area.And set each wind-powered electricity generation Each wind farm wind velocity is essentially identical in synchronization inside section.It contributes using random difference equation simulation wind power plant timing When, for each period, while generating and considering that 17 wind speed of correlation correspond to different wind areas.
Then, the wind velocity distributing paremeter in the area Jiang Gefeng is distributed by Weibull, according to the history of Shandong Province's wind power plant power output number According to being fitted the scale parameter c and form parameter k and auto-correlation attenuation coefficient of each wind area wind speed profile, and different wind area wind Related coefficient between speed distribution.Wherein, wind farm wind velocity parameter in Shandong Province's is as shown in table 2, and Shandong Province's wind farm wind velocity is related Coefficient is as shown in table 3, and the selection of fan parameter is as shown in table 1.
2 Shandong Province's wind farm wind velocity parameter of table
3 Shandong Province's wind farm wind velocity related coefficient of table
With continued reference to Fig. 1 it is found that the method that photovoltaic running simulation boundary condition is set in the present embodiment mainly includes as follows Process:
S21: setting photovoltaic plant parameter, wherein photovoltaic plant parameter be used for calculate it is unobstructed under the conditions of the theoretical sun Irradiation intensity and photovoltaic plant power output.
According to the power output model of photovoltaic panel, photovoltaic power output is directly proportional to practical solar irradiation intensity under the conditions of unobstructed, When solar irradiation intensity is standard irradiation intensity, photovoltaic panel is nominal output.
For calculate it is unobstructed under the conditions of the photovoltaic plant parameter of theoretical solar irradiation intensity include: geographical location, light Photovoltaic array type, photovoltaic array tilt angle and photovoltaic array deflection.Geographical location includes: photovoltaic plant position Longitude, latitude and height above sea level, the same photovoltaic plant section uses identical geographical location in the present embodiment.Photovoltaic array Type includes: fixed angle, horizontal uniaxial, inclination angle single shaft and twin shaft.Photovoltaic array tilt angle is set as photovoltaic plant piece Latitude where area.Photovoltaic array deflection is set as 0.
Photovoltaic plant parameter for calculating photovoltaic plant processing includes: installed capacity, efficiency and prediction absolute error Account for installed capacity percentage.Wherein, installed capacity includes the rated power of photovoltaic panel quantity and photovoltaic panel, and efficiency is photovoltaic panel Availability, prediction absolute error account for installed capacity percentage and are set as 15%.
S22: setting clearness index distribution parameter, wherein clearness index distribution parameter is related to photovoltaic power generation for reflecting Clearness index.
Clearness index refers to practical solar radiation slightly and the ratio of theoretical intensity of solar radiation, clearness index usually with The external world such as temperature, cloud cover and Changes in weather uncertain factor is related, and the modeling Simulation of clearness index needs in the present embodiment Set clearness index distribution parameter.Specifically, clearness index distribution parameter includes: situation exponential probability distribution parameter C, situation Exponential probability distribution parameter lamda value, clearness index probability distribution parameters theta value, atmospheric scattering coefficient p, atmospheric scattering system Number q value and bottom reflection rate.
The setting value of clearness index distribution parameter in the present embodiment, distribution are as follows:
Situation exponential probability distribution parameter C is set as 0.2994, and situation exponential probability distribution parameter lamda value is set as 5.062, clearness index probability distribution parameters theta value is set as 0.0343, and atmospheric scattering coefficient p is set as 1.0303, atmosphere Scattering coefficient q value is set as 1.1515, and bottom reflection rate is set as 0.2.
Weather typing parameter can be obtained by being counted to weather history situation in the present embodiment, and the present embodiment is adopted With weather typing parameter as shown in table 4 below:
4 weather typing parameter of table
It is of course also possible to be compared with the photovoltaic of design using hourage by simulation photovoltaic using hourage, to table 4 In set parameter revised, the accuracy of weather typing parameter setting can be further increased.
S23: according to the distance between photovoltaic plant section, the related coefficient between photovoltaic plant section is set.
Specifically, step S23 is comprised the following processes:
S231: whether judge the distance between any two photovoltaic plant section >=200km.
If the distance between any two photovoltaic plant section >=200km, step S232: any two photovoltaic electric is executed The related coefficient stood between section is set as 0.15.
If the distance between any two photovoltaic plant section < 200km, step S233: any two photovoltaic electric is executed The related coefficient stood between section is set as 0.4.
Related coefficient in the present embodiment between the photovoltaic power output occlusion coefficient of each observation point is referring to fig. 2.
Below by taking the planning of Shandong Power photovoltaic power generation as an example, the method for description photovoltaic running simulation boundary condition setting:
Firstly, photovoltaic plant is divided into 17 photovoltaic plant sections according to photovoltaic plant affiliated districts and cities, it may be assumed that light area, respectively Are as follows: the area Ji Nanguang, the area Qing Daoguang, the area Zi Boguang, the area Zao Zhuanguang, Dongying light area, the area Yan Taiguang, the area Wei Fangguang, the area Ji Ningguang, Thailand The area An Guangqu, Wei Haiguang, day illumination area, the area Bin Zhouguang, the area De Zhouguang, the area Liao Chengguang, the area Lin Yiguang, the area He Zeguang, Laiwu light Area.And set in each photovoltaic plant section weather conditions i.e.: photovoltaic occlusion coefficient, geographical position essentially identical in synchronization It sets identical.The power output of photovoltaic plant under theoretical case can be calculated using geographical location, it may be assumed that the photovoltaic plant of fine day state goes out Power when simulating photovoltaic occlusion coefficient using random difference equation, for each period, while generating 17 light for considering correlation Volt occlusion coefficient corresponds to different light areas.
Then, force data, meteorologic parameter representative value are gone out according to the geographical location of Shandong photovoltaic plant and existing history Setting uses clearness index distribution parameter as described above, each smooth area's geographic factor and photovoltaic plant optimum angle of incidence setting such as table 5 Shown, each smooth area's weather pattern distribution probability is as shown in table 6, and the setting of photovoltaic plant technology parameter is as shown in table 7, and each region hides It is as shown in table 8 to keep off factor related coefficient.
Each smooth area's geographic factor of table 5 and the setting of photovoltaic plant optimum angle of incidence
Each smooth area's weather pattern distribution probability of table 6
7 Shandong Province's photovoltaic plant technology parameter setting table of table
Each regional occlusion factor related coefficient of 8 Shandong Province's photovoltaic plant of table
The above is only the specific embodiment of the application, is made skilled artisans appreciate that or realizing this Shen Please.Various modifications to these embodiments will be apparent to one skilled in the art, as defined herein General Principle can be realized in other embodiments without departing from the spirit or scope of the application.Therefore, the application It is not intended to be limited to the embodiments shown herein, and is to fit to and the principles and novel features disclosed herein phase one The widest scope of cause.

Claims (8)

1. a kind of setting method of new energy running simulation boundary condition, which is characterized in that the described method includes: wind-powered electricity generation runs mould Quasi- boundary condition setting and the setting of photovoltaic running simulation boundary condition;
Wherein, the method for the wind-powered electricity generation running simulation boundary condition setting, comprising:
According to geographical location, all wind power plants in region to be measured are divided into multiple wind areas, include multiple wind-powered electricity generations in each wind area ?;
According to historical wind speed data, wind farm wind velocity statistical distribution parameter is set, the historical wind speed data is current institute, wind area The average value of the historical wind speed data containing wind power plant;
According to the distance between historical wind speed data and wind area, the wind speed related coefficient between multiple wind areas is set;
According to blower model, blower power producing characteristics parameter is set;
The method of the photovoltaic running simulation boundary condition setting, comprising:
Set photovoltaic plant parameter, the photovoltaic plant parameter be used for calculate it is unobstructed under the conditions of theoretical solar irradiation intensity and Photovoltaic plant power output;
Clearness index distribution parameter is set, the clearness index distribution parameter is for reflecting that clear sky relevant to photovoltaic power generation refers to Number;
According to the distance between photovoltaic plant section, the related coefficient between photovoltaic plant section is set.
2. a kind of setting method of new energy running simulation boundary condition according to claim 1, which is characterized in that described According to historical wind speed data, wind farm wind velocity statistical distribution parameter is set, comprising:
Calculate the fitting wind speed statistical distribution parameter of each wind power plant in any wind area in multiple wind areas;
According to the capacity of each wind power plant, the fitting wind speed statistical distribution parameter is taken into weighted average, is calculated The wind speed statistical distribution parameter in any wind area.
3. a kind of setting method of new energy running simulation boundary condition according to claim 2, which is characterized in that described Calculate the fitting wind speed statistical distribution parameter of each wind power plant in any wind area in multiple wind areas, comprising:
According to the historical wind speed data of wind power plant, the scale parameter being distributed using Weibull distribution function fitting wind farm wind velocity, Form parameter and auto-correlation attenuation coefficient;
According to the historical wind speed data of wind power plant, the moon characteristic mark curve and day characteristic mark curve of wind farm wind velocity is calculated.
4. a kind of setting method of new energy running simulation boundary condition according to claim 2, which is characterized in that described According to historical wind speed data, wind farm wind velocity statistical distribution parameter is set, further includes:
When there is no historical wind speed data in any wind area, appoint using the wind speed statistical distribution parameter in adjacent wind area as described The wind speed statistical distribution parameter in one wind area, the adjacent sectors are the wind areas nearest on geographical location with any wind area.
5. a kind of setting method of new energy running simulation boundary condition according to claim 1, which is characterized in that described According to the distance between historical wind speed data and wind area, the wind speed related coefficient between multiple wind areas is set, comprising:
Judge whether current wind area has historical wind speed data;
If so, being fitted the wind speed related coefficient between any two wind area according to the historical wind speed data in current wind area;
If not, being utilized according to the distance between any two wind areaThe wind between any two wind area is calculated Fast related coefficient, wherein c is wind speed related coefficient, and d is the distance between any two wind area, M be wind speed related coefficient with away from From decay factor.
6. a kind of setting method of new energy running simulation boundary condition according to claim 1, which is characterized in that described Photovoltaic plant parameter includes: geographical location, photovoltaic array type, photovoltaic array tilt angle, photovoltaic array deflection, installation appearance Amount, efficiency and prediction absolute error account for installed capacity percentage, and the geographical location, photovoltaic array type, photovoltaic array incline Rake angle and photovoltaic array deflection, for calculate it is unobstructed under the conditions of theoretical solar irradiation intensity, the installed capacity, effect Rate and prediction absolute error account for installed capacity percentage, for calculating photovoltaic plant power output.
7. a kind of setting method of new energy running simulation boundary condition according to claim 1, which is characterized in that described Clearness index distribution parameter, comprising: clearness index probability distribution parameters, clearness index probability distribution parameters lamda value, clear sky refer to Number probability distribution parameters theta value, atmospheric scattering coefficient p, atmospheric scattering coefficient q value and bottom reflection rate.
8. a kind of setting method of new energy running simulation boundary condition according to claim 1, which is characterized in that described According to the distance between photovoltaic plant section, the related coefficient between photovoltaic plant section is set, comprising:
Whether judge the distance between any two photovoltaic plant section >=200km;
If so, the related coefficient between any two photovoltaic plant section is set as 0.15;
If not, the related coefficient between any two photovoltaic plant section is set as 0.4.
CN201811531549.8A 2018-12-14 2018-12-14 Setting method for new energy operation simulation boundary conditions Active CN109599899B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811531549.8A CN109599899B (en) 2018-12-14 2018-12-14 Setting method for new energy operation simulation boundary conditions

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811531549.8A CN109599899B (en) 2018-12-14 2018-12-14 Setting method for new energy operation simulation boundary conditions

Publications (2)

Publication Number Publication Date
CN109599899A true CN109599899A (en) 2019-04-09
CN109599899B CN109599899B (en) 2022-04-19

Family

ID=65960748

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811531549.8A Active CN109599899B (en) 2018-12-14 2018-12-14 Setting method for new energy operation simulation boundary conditions

Country Status (1)

Country Link
CN (1) CN109599899B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112395764A (en) * 2020-11-19 2021-02-23 国网江西省电力有限公司建设分公司 Novel simulation calculation method for new energy photovoltaic output
CN112448390A (en) * 2020-11-13 2021-03-05 梅雪峰 Distributed photovoltaic upscale virtual equivalent power station definition method based on power grid structure and meteorological characteristics
CN115392735A (en) * 2022-08-30 2022-11-25 浙江正泰智维能源服务有限公司 Method, system, equipment and medium for monitoring working performance of photovoltaic power station
CN115392735B (en) * 2022-08-30 2024-04-23 浙江正泰智维能源服务有限公司 Photovoltaic power station working performance monitoring method, system, equipment and medium

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104794325A (en) * 2015-03-10 2015-07-22 国家电网公司 Colony wind power plant output timing sequence simulation method based on random difference equation
CN105005937A (en) * 2015-04-28 2015-10-28 国家电网公司 Photovoltaic power station output sequential simulation method based on clearness indexes
US20160203245A1 (en) * 2015-01-14 2016-07-14 State Grid Corporation Of China Method for simulating wind field of extreme arid region based on wrf
CN106372801A (en) * 2016-08-31 2017-02-01 江西省电力设计院 Wind power field wind region division method and method for selecting wind detection tower position of wind power field wind regions
CN108022019A (en) * 2017-12-14 2018-05-11 西南石油大学 Wind power plant dispatching method and system based on Wind turbines classification
CN108429271A (en) * 2017-02-14 2018-08-21 中国电力科学研究院 The micro- energy optimization method of mixed energy storage system and system adjusted for frequency

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160203245A1 (en) * 2015-01-14 2016-07-14 State Grid Corporation Of China Method for simulating wind field of extreme arid region based on wrf
CN104794325A (en) * 2015-03-10 2015-07-22 国家电网公司 Colony wind power plant output timing sequence simulation method based on random difference equation
CN105005937A (en) * 2015-04-28 2015-10-28 国家电网公司 Photovoltaic power station output sequential simulation method based on clearness indexes
CN106372801A (en) * 2016-08-31 2017-02-01 江西省电力设计院 Wind power field wind region division method and method for selecting wind detection tower position of wind power field wind regions
CN108429271A (en) * 2017-02-14 2018-08-21 中国电力科学研究院 The micro- energy optimization method of mixed energy storage system and system adjusted for frequency
CN108022019A (en) * 2017-12-14 2018-05-11 西南石油大学 Wind power plant dispatching method and system based on Wind turbines classification

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112448390A (en) * 2020-11-13 2021-03-05 梅雪峰 Distributed photovoltaic upscale virtual equivalent power station definition method based on power grid structure and meteorological characteristics
CN112395764A (en) * 2020-11-19 2021-02-23 国网江西省电力有限公司建设分公司 Novel simulation calculation method for new energy photovoltaic output
CN115392735A (en) * 2022-08-30 2022-11-25 浙江正泰智维能源服务有限公司 Method, system, equipment and medium for monitoring working performance of photovoltaic power station
CN115392735B (en) * 2022-08-30 2024-04-23 浙江正泰智维能源服务有限公司 Photovoltaic power station working performance monitoring method, system, equipment and medium

Also Published As

Publication number Publication date
CN109599899B (en) 2022-04-19

Similar Documents

Publication Publication Date Title
Li et al. Onshore and offshore wind energy potential assessment near Lake Erie shoreline: A spatial and temporal analysis
Bukala et al. Modern small wind turbine design solutions comparison in terms of estimated cost to energy output ratio
Veigas et al. Wave and offshore wind potential for the island of Tenerife
TWI476430B (en) A wind energy forecasting method with extreme wind speed prediction function
Asaah et al. Optimal placement of wind turbines in wind farm layout using particle swarm optimization
Marinelli et al. Wind and photovoltaic large-scale regional models for hourly production evaluation
Ohunakin Wind resource evaluation in six selected high altitude locations in Nigeria
CN107292514B (en) Site selection method and device for anemometer tower in production operation of wind power plant
Ataei et al. Optimum design of an off-grid hybrid renewable energy system for an office building
Weisser et al. Implications of seasonal and diurnal variations of wind velocity for power output estimation of a turbine: a case study of Grenada
CN103258242A (en) Wind measurement network layout method based on wind power plant layout in large-scale wind power base
Reed Wind power climatology of the United States: supplement
Garcia-Heller et al. Forecast study of the supply curve of solar and wind technologies in Argentina, Brazil, Chile and Mexico
Reja et al. A review of the evaluation of urban wind resources: Challenges and perspectives
Lindberg et al. Review on power-production modeling of hybrid wind and PV power parks
CN109599899A (en) A kind of setting method of new energy running simulation boundary condition
CN109449987A (en) The power system simulation method and system of new energy power generation grid-connection
Winslow Urban wind generation: comparing horizontal and vertical axis wind turbines at Clark University in Worcester, Massachusetts
KR20200040448A (en) Apparatus and method for predicting solar energy generation amount
Cutler et al. Characterizing future large, rapid changes in aggregated wind power using numerical weather prediction spatial fields
CN109636019A (en) Method is determined based on the anemometer tower arrangement of neural network algorithm
Mohammadi et al. Modeling the allocation and economic evaluation of PV panels and wind turbines in urban areas
Karamanis Management of moderate wind energy coastal resources
Muchiri et al. Wind and solar resource complementarity and its viability in wind/PV hybrid energy systems in Machakos, Kenya
Belu Assessment and analysis of offshore wind energy potential

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