CN113987823A - Mesoscale numerical simulation method for land centralized photovoltaic power station climate effect evaluation - Google Patents
Mesoscale numerical simulation method for land centralized photovoltaic power station climate effect evaluation Download PDFInfo
- Publication number
- CN113987823A CN113987823A CN202111316011.7A CN202111316011A CN113987823A CN 113987823 A CN113987823 A CN 113987823A CN 202111316011 A CN202111316011 A CN 202111316011A CN 113987823 A CN113987823 A CN 113987823A
- Authority
- CN
- China
- Prior art keywords
- power station
- photovoltaic power
- photovoltaic
- wave radiation
- heat flux
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 50
- 230000000694 effects Effects 0.000 title claims abstract description 47
- 238000004088 simulation Methods 0.000 title claims abstract description 44
- 238000011156 evaluation Methods 0.000 title claims abstract description 7
- 230000005855 radiation Effects 0.000 claims abstract description 62
- 230000004907 flux Effects 0.000 claims abstract description 57
- 239000002689 soil Substances 0.000 claims abstract description 22
- 238000012360 testing method Methods 0.000 claims abstract description 18
- 238000006243 chemical reaction Methods 0.000 claims abstract description 17
- 230000035945 sensitivity Effects 0.000 claims abstract description 13
- 238000011158 quantitative evaluation Methods 0.000 claims abstract description 12
- 238000007781 pre-processing Methods 0.000 claims abstract description 8
- 230000009471 action Effects 0.000 claims abstract description 7
- 238000013016 damping Methods 0.000 claims description 11
- 238000001514 detection method Methods 0.000 claims description 5
- 238000010438 heat treatment Methods 0.000 claims description 4
- 230000010354 integration Effects 0.000 claims description 3
- 238000004379 similarity theory Methods 0.000 claims description 3
- 238000012353 t test Methods 0.000 claims description 3
- 238000000692 Student's t-test Methods 0.000 claims description 2
- 238000007689 inspection Methods 0.000 claims description 2
- 230000007935 neutral effect Effects 0.000 claims description 2
- 238000010248 power generation Methods 0.000 claims description 2
- 238000003908 quality control method Methods 0.000 claims description 2
- 230000001360 synchronised effect Effects 0.000 claims description 2
- 239000000523 sample Substances 0.000 claims 1
- 238000010276 construction Methods 0.000 abstract description 3
- 230000008859 change Effects 0.000 description 5
- YLPGTOIOYRQOHV-UHFFFAOYSA-N Pretilachlor Chemical compound CCCOCCN(C(=O)CCl)C1=C(CC)C=CC=C1CC YLPGTOIOYRQOHV-UHFFFAOYSA-N 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000012827 research and development Methods 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Geometry (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Photovoltaic Devices (AREA)
Abstract
The invention discloses a mesoscale numerical simulation method for onshore centralized photovoltaic power station climate effect evaluation, which comprises the steps of preprocessing multisource meteorological observation data of an onshore centralized photovoltaic power station; carrying out surface short wave radiation parameterization by utilizing the comprehensive albedo and the photoelectric conversion efficiency; directly and explicitly carrying out numerical parameterization on the sensible heat flux by utilizing wind speed, air density and short-wave radiation; latent heat flux, soil heat flux and long wave radiation parameterization are carried out by using a difference contrast method; parameterizing the dragging action of the photovoltaic panel by using the dynamic roughness; and quantitatively evaluating the climate effect of the photovoltaic power station through a sensitivity numerical simulation test. The method solves the problem of objective quantitative evaluation of climate influence of the large photovoltaic power station, and the simulation method has strong universality. By adopting the method, the climate effect of the built photovoltaic power station can be post-evaluated, the climate effect of the proposed power station can be pre-evaluated, and a basis is provided for site selection construction of the climate-friendly power station.
Description
Technical Field
The invention belongs to the technical field of meteorological data analysis and processing, relates to a numerical simulation method for quantitative evaluation of climate effect of a photovoltaic power station, and particularly relates to a land surface process numerical parameterization method based on radiation balance, heat balance and surface dragging action in a photovoltaic power station area.
Background
The quantitative evaluation of the climate effect of the centralized photovoltaic power station not only provides necessary technical support for the site selection design of the climate-friendly photovoltaic power station, but also is the scientific basis for energy structure adjustment and climate change coping in China. At present, the research on the climate effect of the large-scale photovoltaic power station development is only limited to fixed-point observation, theoretical calculation or extensive simulation, and can not provide effective support for the climate adaptation improvement of the photovoltaic power station development and the future development planning of the photovoltaic industry. Therefore, the explicit complete numerical simulation method for the climate effect of the photovoltaic power station is one of core technical tools for scientifically and quantitatively evaluating the climate effect of the photovoltaic power station.
The key link of the mesoscale numerical simulation method for the climate effect of the photovoltaic power station is to perform numerical description on physical processes such as radiation balance, heat balance and surface dragging of the area where the photovoltaic power station is located, establish an explicit and complete mesoscale numerical parameterization tool and further develop a sensitivity numerical simulation test. At present, research and development work aiming at photovoltaic power station climate effect numerical parameterization is started late. Although the american arizona group establishes an energy balance numerical model for different photovoltaic tracking support types from the perspective of micro-scale energy balance of photovoltaic panels, the model is a micro-scale model that needs external meteorological variables to drive and operates independently, and does not realize coupling with a meso-scale meteorological numerical model, and has limited application in quantitative evaluation of climate effects of large-scale photovoltaic power stations. A domestic research and development team preliminarily develops a rudiment tool suitable for carrying out large-scale laying of the climate effect evaluation of the photovoltaic module, but the rudiment tool does not express the daily cycle characteristic of sensible heat flux explicitly, neglects the influence of soil heat flux, latent heat flux and upward long wave on the surface energy balance, and meanwhile, the dynamics rough length in the numerical simulation tool is an empirical value and is difficult to truly reflect the energy balance characteristic of a photovoltaic power station area.
Disclosure of Invention
Aiming at the technical problems in the existing numerical simulation method, the invention provides a mesoscale numerical simulation method for onshore centralized photovoltaic power station climate effect evaluation, which comprises the steps of preprocessing multisource meteorological observation data of an onshore centralized photovoltaic power station; carrying out surface short wave radiation parameterization by utilizing the comprehensive albedo and the photoelectric conversion efficiency; directly and explicitly carrying out numerical parameterization on the sensible heat flux by utilizing wind speed, air density and short-wave radiation; latent heat flux, soil heat flux and long wave radiation parameterization are carried out by using a difference contrast method; parameterizing the dragging action of the photovoltaic panel by using the dynamic roughness; and quantitatively evaluating the climate effect of the photovoltaic power station through a sensitivity numerical test. The numerical simulation method for the climate effect of the photovoltaic power station, which is established based on numerical parameterization and sensitivity numerical simulation tests of the earth surface radiation balance, heat balance and surface drag damping effects of the photovoltaic power station, solves the problem of objective quantitative evaluation of the climate influence of the large photovoltaic power station, and has strong universality. By adopting the numerical simulation method, the climate effect of the built photovoltaic power station can be post-evaluated, the climate effect of the proposed power station can be pre-evaluated, and a basis is provided for the site selection construction of the climate-friendly power station.
The technical scheme adopted by the invention for solving the problem of the numerical simulation method is as follows:
a mesoscale numerical simulation method for onshore centralized evaluation of the climate effects of photovoltaic power plants, said numerical simulation method being developed on the basis of objective quantitative evaluation of the climate effects of large photovoltaic power plants, characterized in that said model comprises at least the following steps:
SS1, preprocessing multi-source meteorological observation data of a terrestrial centralized photovoltaic power station: controlling the quality of original observation data, and selecting synchronous detection data of internal and external measuring points (reference points) of the photovoltaic power station;
and SS2, performing surface short wave radiation parameterization by utilizing the comprehensive albedo and the photoelectric conversion efficiency: calculating the comprehensive albedo of the photovoltaic power station in summer and winter respectively by using the downward and upward short wave radiation observation data after quality control in the photovoltaic power station and adopting a ratio method; meanwhile, the influence of the photovoltaic panel temperature on the power generation efficiency in different seasons is considered, the photoelectric conversion efficiency in summer and the photoelectric conversion efficiency in winter in the photovoltaic power station are respectively calculated by adopting a ratio method according to the generated energy and downward short wave radiation data, and then the short wave radiation balance of the photovoltaic power station area is numerically parameterized;
SS3. direct explicit numerical parameterization of sensible heat flux with wind speed, air density and short wave radiation: the property that the photovoltaic panel heats the atmosphere in the daytime is considered, and the heat-sensitive flux in the photovoltaic power station area is numerically parameterized by combining the heat-sensitive conveying coefficient of the photovoltaic power station, the wind speed after surface damping, the air density, the plate temperature-air temperature difference (described by short-wave radiation), the occupation ratio of the surface coverage type in the photovoltaic power station and the like;
and SS4, carrying out latent heat flux, soil heat flux and long-wave radiation parameterization by using a difference contrast method: calculating the relative changes of latent heat flux, soil heat flux and earth surface temperature observed in the photovoltaic power station at the same period by adopting a difference comparison method and taking an observed value of a reference point as a reference, and carrying out numerical parameterization on the latent heat flux, the soil heat flux and long-wave radiation balance of the photovoltaic power station area by combining the Stefan-Boltzmann thermodynamic law;
ss5. parameterization of the dragging effect of photovoltaic panels with dynamic roughness: calculating the average dynamic roughness Z of the photovoltaic underlying surface by using vortex observation data in the photovoltaic power station according to the Monin-Obukhov similarity theory0Describing the surface drag damping effect of the photovoltaic panel;
SS 6: the climate effect of the photovoltaic power station is quantitatively evaluated through a sensitivity numerical test: based on the steps SS2-SS5, sensitivity numerical simulation tests of the situation of the photovoltaic power station without the photovoltaic power station and the situation of the photovoltaic power station with the photovoltaic power station are respectively carried out, and the time-space difference of key meteorological variables in the two groups of numerical tests is calculated in a comparison mode, so that the quantitative evaluation of the climate effect of the photovoltaic power station is achieved.
Preferably, in step SS1, the data preprocessing includes: removing original observation data such as stiffness values and false values in ground short wave radiation, long wave radiation, sensible heat (latent heat) flux, soil heat flux, average wind speed, surface temperature, turbulence observation data and the like through climate threshold value, space-time consistency inspection and the like; and selecting effective observation data of the photovoltaic power station at the internal and external measuring points in the same period.
Preferably, in the step SS2, the comprehensive albedo of the area where the photovoltaic power station is located is obtained by dividing the upward short-wave radiation observed in winter/summer by the downward short-wave radiation at the detection point in the photovoltaic power station; calculating the photoelectric conversion efficiency of the photovoltaic panel by dividing the generated energy of the photovoltaic power station in winter/summer by the total solar short wave radiation amount reaching the ground surface; and carrying out numerical parameterization on the short-wave radiation balance of the photovoltaic power station area by a mode of superposing the photoelectric conversion efficiency (alpha + epsilon) through comprehensive albedo.
Preferably, in the step SS3, the sensible heat heating process of the photovoltaic panel is explicitly calculated by using the sensible heat transfer coefficient, the average wind speed, the air density, the board temperature-air temperature difference (proportional to the downward short wave radiation reaching the photovoltaic panel) and the air constant pressure specific heat constant in the photovoltaic power station area; on the basis, area weighted average is carried out on the sensible heat flux of the photovoltaic panel and the sensible heat flux of the natural earth surface through the coverage area ratio of the photovoltaic panel and the natural earth surface in the photovoltaic power station area, and a numerical parameterization scheme of the sensible heat flux of the photovoltaic power station area is established.
Preferably, in step SS4, based on the observation value of the reference point outside the photovoltaic power station, the variation ratio of the observation value of the measurement point in the photovoltaic power station with respect to the observation value of the reference point is calculated by using a difference statistical method, and the latent heat flux, the soil heat flux and the long-wave radiation balance in the photovoltaic power station area are parameterized numerically by means of geometric adjustment.
Preferably, in the step SS5, turbulence data observed by the vortex correlation system in the photovoltaic power station is applied to a profile fitting function of the near-formation wind speed, and under the condition of a neutral atmospheric junction, the average dynamic roughness Z of the photovoltaic underlying surface is calculated and obtained0And numerically parameterizing the surface dynamic dragging action of the photovoltaic panel.
Preferably, in the step SS6, in the target simulation area, two sets of numerical tests of photovoltaic power station-free situation and photovoltaic power station-containing situation are respectively designed, and numerical simulation integration is respectively performed based on the mesoscale numerical parameterization tools of photovoltaic power station ground surface radiation balance, heat balance and surface drag damping established by SS2-SS 5. And comparing and calculating the space-time difference of the two groups of simulation results, and performing significance analysis on the difference by using a statistical t test method so as to objectively and quantitatively evaluate the climate effect of the photovoltaic power station.
Compared with the prior art, the invention has the beneficial effects that:
the numerical simulation method for the climate effect of the photovoltaic power station, which is established based on numerical parameterization and sensitivity numerical tests of earth surface radiation balance, heat balance and surface drag damping effects of the photovoltaic power station, solves the problem of objective quantitative evaluation of climate influence of large photovoltaic power stations, adopts an explicit complete technical method to carry out numerical parameterization on core physical processes of short wave radiation balance, sensible heat heating, surface power drag effects and the like of the photovoltaic power station, effectively reduces the high dependence of a model on the position background climate characteristics of the photovoltaic power station, and has strong universality. By adopting the numerical simulation tool, the climate effect of the built photovoltaic power station can be post-evaluated, the climate effect of the proposed power station can be pre-evaluated, and a basis is provided for the site selection construction of the climate-friendly power station.
Drawings
FIG. 1 is a schematic diagram of a photovoltaic power station climate effect mesoscale numerical simulation method established based on photovoltaic power station energy balance numerical parameterization and sensitivity numerical test.
FIG. 2 is a tower gradient observation system of a centralized photovoltaic power station.
Fig. 3 is a graph of the daily change in soil heat flux in winter (a)/summer (b) for a photovoltaic power plant area.
Detailed Description
In order to make the implementation objects, technical solutions and advantages of the present invention clearer, the technical solutions in the embodiments of the present invention will be described in more detail below with reference to the accompanying drawings in the embodiments of the present invention. The described embodiments are only some, but not all embodiments of the invention. The embodiments described below with reference to the drawings are illustrative and intended to be illustrative of the invention and are not to be construed as limiting the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As shown in fig. 1, the present invention relates to a method for numerically simulating a climate effect of a photovoltaic power station, which is established based on numerical parameterization and sensitivity numerical tests of ground surface radiation balance, heat balance and surface drag damping, and the method at least comprises the following steps:
SS1, preprocessing multi-source meteorological observation data of a terrestrial centralized photovoltaic power station;
SS2, performing surface short wave radiation parameterization by utilizing the comprehensive albedo and the photoelectric conversion efficiency;
SS3, directly and explicitly carrying out numerical parameterization on the sensible heat flux by utilizing wind speed, air density and short-wave radiation;
SS4, carrying out latent heat flux, soil heat flux and long-wave radiation parameterization by using a difference contrast method;
SS5. parameterizing the dragging action of the photovoltaic panel by using the dynamic roughness;
SS 6: and quantitatively evaluating the climate effect of the photovoltaic power station through a sensitivity numerical test.
1. In the step SS1, multi-source meteorological observation data preprocessing of the onshore centralized photovoltaic power station
Aiming at meteorological data detected by a tower type gradient observation system (figure 2) of a centralized photovoltaic power station, such as ground short wave radiation, long wave radiation, sensible heat (latent heat) flux, soil heat flux, average wind speed, surface temperature, turbulence observation data and the like, firstly, eliminating false values exceeding reasonable threshold values in original data through climate threshold values of all meteorological variables; for stiffness values with constant detection values for successive times, by ut+1-utScreening and eliminating (0), wherein ut+1Representing the detected value, u, at time t +1tRepresenting the detected value at time t. Finally, by means of a contemporaneous pairAnd selecting effective observation data of the internal and external measuring points of the photovoltaic power station at the same time in a comparison mode.
2. In the step SS2, the comprehensive albedo and the photoelectric conversion efficiency are used for parameterizing the surface short wave radiation
Dividing upward short wave radiation (SWUP, shown in figure 1) observed by a detection point in the photovoltaic power station in winter/summer by downward short wave radiation (SWDOWN, shown in figure 1) to obtain the comprehensive albedo alpha of the area where the photovoltaic power station is located; generating power (S) by winter/summer photovoltaic power stationPVFig. 1) divided by the solar short wave radiation (SWDOWN, fig. 1) reaching the surface, the photoelectric conversion efficiency epsilon of the photovoltaic panel is calculated; and carrying out numerical parameterization on the short-wave radiation balance of the photovoltaic power station area by a mode of superposing the photoelectric conversion efficiency (alpha + epsilon) through comprehensive albedo. The net short-wave radiation (solnet) of a photovoltaic power station area can therefore be parameterized as follows:
solnet＝SWDOWN×(1-α-ε)
taking the photovoltaic power station of five channels in Xinjiang as an example, the comprehensive albedo of the photovoltaic power station in winter/summer is 0.21/0.15, and the photoelectric conversion efficiency in winter/summer is 0.15/0.13.
3. In the above step SS3, direct explicit numerical parameterization of sensible heat flux is performed by using wind speed, air density and short wave radiation
Explicit calculation of sensible heat heating process of photovoltaic panels (HFX) using sensible heat transport coefficient, average wind speed, air density, board temperature-air temperature difference (proportional to the downward short wave radiation reaching the photovoltaic panel) and air constant pressure specific heat constant within the photovoltaic power plant areaPFig. 1); on the basis, the coverage area ratio of the photovoltaic panel and the natural earth surface in the photovoltaic power station area is used for measuring the heat-sensitive flux (HFX) of the photovoltaic panelPFIG. 1) and natural surface sensible Heat Flux (HFX)GFig. 1) area weighted averaging is performed to establish a numerical parameterization scheme of the photovoltaic plant area Heat Flux (HFX):
HFXP＝Ch×Cp×ρ×V×0.35×SWDOWN
HFX＝30％×HFXP×HFXG
under unstable conditions, lightHeat sensitive transport coefficient C of photovoltaic power stationh0.0045; cpIs the constant specific heat constant of air at constant pressure; ρ is the air density and V is the average wind speed, here representing the wind speed after damping of the photovoltaic panel.
4. In the step SS4, latent heat flux, soil heat flux and long wave radiation parameterization is carried out by using a contrast method
Calculating the change proportion of a measuring point observation value in the photovoltaic power station relative to a reference point observation value by using a difference value statistical method based on a reference point observation value outside the photovoltaic power station, and carrying out numerical parameterization on latent heat flux (lambda E, figure 1), soil heat flux (G, figure 1) and upward long-wave radiation (LWUP, figure 1) of a photovoltaic power station area in an equal ratio adjustment mode:
λE＝λE×(1-β)
g ═ G × (1- γ) (winter)
G ═ G × (1+ θ) (summer)
LWUP＝emi×stbolt×(LST-ΔT)4
Due to the physical shielding and wind resistance effects of the photovoltaic panel, the latent heat evaporation in the photovoltaic power station area is reduced, and the latent heat reduction ratio beta is 0.47 on the daily average scale by taking five channels as an example; the soil heat flux of the photovoltaic power station area shows different change characteristics in winter and summer, for example, five ditches are taken as examples, the soil heat flux in winter is reduced, the soil heat flux reduction ratio gamma is 0.22 (fig. 3a) on the average daily scale, the soil heat flux in summer shows an increase characteristic, the increase ratio theta in daytime is 0.34, and the increase ratio theta in night is 0.17 (fig. 3 b); similarly, the surface temperature of the photovoltaic power station area will decrease due to the physical shading of the panels, for example, five channels in Xinjiang, the difference Δ T between the surface temperature in winter/summer and the reference point is 0.4/1.2, and the change in surface temperature causes a decrease in upward long wave according to Stefin-Boltzmann thermodynamics law.
Based on the above steps SS2, SS3 and SS4, the surface energy balance of the centralized photovoltaic power plant area can be fully parameterized by an energy residual equation:
respv＝solnet+LWDOWN-LWUP-HFX-λE-G
it should be noted that the net short wave radiation solnet containsEnergy S of light to electricitypv。
5. In the above step SS5, the dragging action of the photovoltaic panel is parameterized by the dynamic roughness
Calculating the average dynamic roughness Z of the photovoltaic underlying surface by using the vortex observation data of the photovoltaic power station area according to the Monin-Obukhov similarity theory0The method is used for describing the surface drag damping effect of the photovoltaic panel.
6. Quantitative evaluation of climate effects of photovoltaic power stations by sensitivity numerical tests
In a target simulation area, two groups of numerical tests of a photovoltaic power station scene without and a photovoltaic power station scene with are respectively designed, and simulation integration is respectively carried out on a mesoscale numerical parameterization tool of photovoltaic power station earth surface radiation balance, heat balance and surface drag damping based on SS2-SS 5. And comparing and calculating the space-time difference of the two groups of simulation results, and performing significance analysis on the difference by using a statistical t-test method as follows:
wherein the content of the first and second substances,
is the average of the sequence of differences of the two sets of numerical simulation results, SDIs the standard deviation of the sequence of difference values. Respectively calculating the t value at each simulation lattice point by the formula and comparing the t value with the t0.1Corresponding statistical values are compared to identify grid points where spatiotemporal differences pass 90% significance. And finally, realizing the quantitative evaluation of the climate effect of the photovoltaic power station.
The present invention is not limited to the above embodiments, and any changes or substitutions that can be easily made by those skilled in the art within the technical scope of the present invention are also within the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (7)
1. A mesoscale numerical simulation method for onshore centralized evaluation of the climate effects of photovoltaic power plants, said numerical simulation method being developed on the basis of objective quantitative evaluation of the climate effects of large photovoltaic power plants, characterized in that said model comprises at least the following steps:
SS1, preprocessing multi-source meteorological observation data of a terrestrial centralized photovoltaic power station: controlling the quality of original observation data, and selecting synchronous detection data of internal and external measuring points (reference points) of the photovoltaic power station;
and SS2, performing surface short wave radiation parameterization by utilizing the comprehensive albedo and the photoelectric conversion efficiency: calculating the comprehensive albedo of the photovoltaic power station in summer and winter respectively by using the downward and upward short wave radiation observation data after quality control in the photovoltaic power station and adopting a ratio method; meanwhile, the influence of the photovoltaic panel temperature on the power generation efficiency in different seasons is considered, the photoelectric conversion efficiency in summer and the photoelectric conversion efficiency in winter in the photovoltaic power station are respectively calculated by adopting a ratio method according to the generated energy and downward short wave radiation data, and then the short wave radiation balance of the photovoltaic power station area is numerically parameterized;
SS3. direct explicit numerical parameterization of sensible heat flux with wind speed, air density and short wave radiation: the property that the photovoltaic panel heats the atmosphere in the daytime is considered, and the heat-sensitive flux in the photovoltaic power station area is numerically parameterized by combining the heat-sensitive conveying coefficient of the photovoltaic power station, the wind speed after surface damping, the air density, the plate temperature-air temperature difference (described by short-wave radiation), the occupation ratio of the surface coverage type in the photovoltaic power station and the like;
and SS4, carrying out latent heat flux, soil heat flux and long-wave radiation parameterization by using a difference contrast method: calculating the relative changes of latent heat flux, soil heat flux and earth surface temperature observed in the photovoltaic power station at the same period by adopting a difference comparison method and taking an observed value of a reference point as a reference, and carrying out numerical parameterization on the latent heat flux, the soil heat flux and long-wave radiation balance of the photovoltaic power station area by combining the Stefan-Boltzmann thermodynamic law;
ss5. parameterization of the dragging effect of photovoltaic panels with dynamic roughness: by using vortex observation data in the photovoltaic power station according to Monin-Obukhov similarity theory, calculating average dynamic roughness Z of the photovoltaic underlying surface0Describing the surface drag damping effect of the photovoltaic panel;
SS 6: the climate effect of the photovoltaic power station is quantitatively evaluated through a sensitivity numerical test: based on the steps SS2-SS5, sensitivity numerical simulation tests of the situation of the photovoltaic power station without the photovoltaic power station and the situation of the photovoltaic power station with the photovoltaic power station are respectively carried out, and the time-space difference of key meteorological variables in the two groups of numerical tests is calculated in a comparison mode, so that the quantitative evaluation of the climate effect of the photovoltaic power station is achieved.
2. A numerical simulation method according to the preceding claim, characterized in that in said step SS1, said data preprocessing comprises: removing original observation data such as stiffness values and false values in ground short wave radiation, long wave radiation, sensible heat (latent heat) flux, soil heat flux, average wind speed, surface temperature, turbulence observation data and the like through climate threshold value, space-time consistency inspection and the like; and selecting effective observation data of the photovoltaic power station at the internal and external measuring points in the same period.
3. A numerical simulation method according to the preceding claims, characterized in that in step SS2, the integrated albedo of the area in which the photovoltaic plant is located is obtained by dividing the upward short-wave radiation observed in winter/summer by the downward short-wave radiation of the probe points inside the photovoltaic plant; calculating the photoelectric conversion efficiency of the photovoltaic panel by dividing the generated energy of the photovoltaic power station in winter/summer by the total solar short wave radiation amount reaching the ground surface; and carrying out numerical parameterization on the short-wave radiation balance of the photovoltaic power station area by a mode of superposing the photoelectric conversion efficiency (alpha + epsilon) through comprehensive albedo.
4. A numerical simulation method according to the previous claims, characterized in that in step SS3, the sensible heat heating process of the photovoltaic panel is explicitly calculated using the sensible heat transport coefficient, the average wind speed, the air density, the board temperature-air temperature difference (proportional to the short-wave radiation down to the photovoltaic panel) and the air constant pressure specific heat constant in the photovoltaic plant area; on the basis, area weighted average is carried out on the sensible heat flux of the photovoltaic panel and the sensible heat flux of the natural earth surface through the coverage area ratio of the photovoltaic panel and the natural earth surface in the photovoltaic power station area, and a numerical parameterization scheme of the sensible heat flux of the photovoltaic power station area is established.
5. A numerical simulation method according to the preceding claims, characterized in that in step SS4, the variation ratio of the observation value of the measuring point in the photovoltaic power station relative to the observation value of the reference point is calculated by using a difference statistics method based on the observation value of the reference point outside the photovoltaic power station, and the latent heat flux, the soil heat flux and the long wave radiation balance of the photovoltaic power station area are parameterized numerically by means of geometric adjustment.
6. The numerical simulation method of the preceding claims, wherein in step SS5, turbulence data observed by a vortex-related system in a photovoltaic power plant is applied to a profile fitting function of the near-formation wind speed, and under the condition of a neutral atmospheric junction, the average dynamic roughness Z of the photovoltaic underlying surface is calculated0And numerically parameterizing the surface dynamic dragging action of the photovoltaic panel.
7. The numerical simulation method of the previous claims, wherein in the step SS6, two sets of numerical simulation tests of photovoltaic power station-free situation and photovoltaic power station-available situation are respectively designed in a target simulation area, based on the steps SS2-SS5, numerical simulation integration is respectively carried out, the space-time difference of two sets of simulation results is calculated through comparison, and the difference is analyzed significantly by using a statistical t-test method, so that the climate effect of the photovoltaic power station is objectively and quantitatively evaluated.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111316011.7A CN113987823B (en) | 2021-11-08 | 2021-11-08 | Mesoscale numerical simulation method for land centralized photovoltaic power station climate effect evaluation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111316011.7A CN113987823B (en) | 2021-11-08 | 2021-11-08 | Mesoscale numerical simulation method for land centralized photovoltaic power station climate effect evaluation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113987823A true CN113987823A (en) | 2022-01-28 |
| CN113987823B CN113987823B (en) | 2022-04-12 |
Family
ID=79747179
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111316011.7A Active CN113987823B (en) | 2021-11-08 | 2021-11-08 | Mesoscale numerical simulation method for land centralized photovoltaic power station climate effect evaluation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113987823B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114964489A (en) * | 2022-05-31 | 2022-08-30 | 中国水利水电科学研究院 | Observation system and method for influencing water surface to receive short-wave radiation energy through photovoltaic array |
| CN114996943A (en) * | 2022-06-06 | 2022-09-02 | 国家气候中心 | Mesoscale numerical simulation method for reservoir storage climate effect evaluation |
| CN115049519A (en) * | 2022-05-31 | 2022-09-13 | 中国水利水电科学研究院 | Method for calculating influence of overwater photovoltaic on environmental water temperature of underlying open water area |
| CN115130395A (en) * | 2022-05-31 | 2022-09-30 | 中国水利水电科学研究院 | Method for calculating influence of overwater photovoltaic on water temperature of underlying closed water area |
| CN117521282A (en) * | 2023-11-07 | 2024-02-06 | 国家气候中心 | Density-dependent fan parameterization method for wind farm climate characteristic simulation |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109726867A (en) * | 2018-12-27 | 2019-05-07 | 北京恒泰实达科技股份有限公司 | A kind of high-resolution electric power method of meteorological forecast based on multi-mode set |
-
2021
- 2021-11-08 CN CN202111316011.7A patent/CN113987823B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109726867A (en) * | 2018-12-27 | 2019-05-07 | 北京恒泰实达科技股份有限公司 | A kind of high-resolution electric power method of meteorological forecast based on multi-mode set |
Non-Patent Citations (5)
| Title |
|---|
| RUI CHANG ET AL.: "Observed surface radiation and temperature impacts from the largescale deployment of photovoltaics in the barren area of Gonghe, China", 《RENEWABLE ENERGY》 * |
| RUI CHANG.ET AL.: "Simulated local climatic impacts of large-scale photovoltaics over the barren area of Qinghai, China", 《RENEWABLE ENERGY》 * |
| 崔杨等: "光伏电站对局地气候的影响研究进展", 《气候变化研究进展》 * |
| 李芬等: "光伏电站建设运行对气候环境的能量影响", 《气象科技进展》 * |
| 郝晓龙: "太阳能光伏屋顶对城市热环境及能源供需影响的模拟", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114964489A (en) * | 2022-05-31 | 2022-08-30 | 中国水利水电科学研究院 | Observation system and method for influencing water surface to receive short-wave radiation energy through photovoltaic array |
| CN115049519A (en) * | 2022-05-31 | 2022-09-13 | 中国水利水电科学研究院 | Method for calculating influence of overwater photovoltaic on environmental water temperature of underlying open water area |
| CN115130395A (en) * | 2022-05-31 | 2022-09-30 | 中国水利水电科学研究院 | Method for calculating influence of overwater photovoltaic on water temperature of underlying closed water area |
| CN115049519B (en) * | 2022-05-31 | 2023-03-28 | 中国水利水电科学研究院 | Method for calculating influence of overwater photovoltaic on environmental water temperature of underlying open water area |
| CN114996943A (en) * | 2022-06-06 | 2022-09-02 | 国家气候中心 | Mesoscale numerical simulation method for reservoir storage climate effect evaluation |
| CN114996943B (en) * | 2022-06-06 | 2022-11-22 | 国家气候中心 | Mesoscale numerical simulation method for reservoir storage climate effect evaluation |
| CN117521282A (en) * | 2023-11-07 | 2024-02-06 | 国家气候中心 | Density-dependent fan parameterization method for wind farm climate characteristic simulation |
| CN117521282B (en) * | 2023-11-07 | 2024-04-12 | 国家气候中心 | Density-dependent fan parameterization method for wind farm climate characteristic simulation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113987823B (en) | 2022-04-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113987823B (en) | Mesoscale numerical simulation method for land centralized photovoltaic power station climate effect evaluation | |
| Lappalainen et al. | Output power variation of different PV array configurations during irradiance transitions caused by moving clouds | |
| Costanzo et al. | A method of strategic evaluation of energy performance of Building Integrated Photovoltaic in the urban context | |
| CN105375878B (en) | A kind of method of on-line checking and assessment photovoltaic system efficiency | |
| CN109492315B (en) | Spatio-temporal correlation wind-solar sequence simulation method based on Copula function | |
| Lappalainen et al. | Analysis of shading periods caused by moving clouds | |
| Sinapis et al. | A comprehensive study on partial shading response of c-Si modules and yield modeling of string inverter and module level power electronics | |
| CN108399493A (en) | Dust stratification causes photovoltaic power generation quantity loss forecasting method and photovoltaic module to clean judgment method | |
| CN113496311A (en) | Photovoltaic power station generated power prediction method and system | |
| CN103761578A (en) | Solar irradiation predicting method based on multiple linear regression | |
| Simioni et al. | Georeferenced operating-efficiency solar potential maps with local weather conditions–An application to Brazil | |
| Li et al. | Improving the energy production of roof-top solar PV systems through roof design | |
| CN105354655A (en) | Photovoltaic power station group confidence capacity evaluation method considering power correlation | |
| Hsieh et al. | Potential for installing photovoltaic systems on vertical and horizontal building surfaces in urban areas | |
| CN107229817B (en) | Photovoltaic power station dust shielding loss analysis method | |
| Zubair et al. | Assessment of photovoltaic capabilities in urban environments: A case study of Islamabad, Pakistan | |
| Chang et al. | A coupled WRF-PV mesoscale model simulating the near-surface climate of utility-scale photovoltaic plants | |
| Yao et al. | The modified ASHRAE model based on the mechanism of multi-parameter coupling | |
| Mazin et al. | Linear and nonlinear modeling for solar energy prediction on the zone, region and global | |
| Rodrigo et al. | Validation of meso-wake models for array efficiency prediction using operational data from five offshore wind farms | |
| Du et al. | Nowcasting methods for optimising building performance | |
| Wang et al. | Neural network based irradiance mapping model of solar PV power forecasting using sky image | |
| Lestari et al. | Prediction of solar energy radiation using adaptive neuro fuzzy inference system in the tropical region | |
| CN112561181A (en) | Photovoltaic power generation prediction system based on Unet network and foundation cloud picture | |
| Zhou et al. | PV power characteristic modeling based on multi-scale clustering and its application in generation prediction |
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 |