CN110414059A - The radiant energy dfensity analogy method of planar heliostats in tower type solar thermo-power station - Google Patents
The radiant energy dfensity analogy method of planar heliostats in tower type solar thermo-power station Download PDFInfo
- Publication number
- CN110414059A CN110414059A CN201910577119.8A CN201910577119A CN110414059A CN 110414059 A CN110414059 A CN 110414059A CN 201910577119 A CN201910577119 A CN 201910577119A CN 110414059 A CN110414059 A CN 110414059A
- Authority
- CN
- China
- Prior art keywords
- heliostat
- radiant energy
- imaging plane
- plane
- dfensity
- 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
Landscapes
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of radiant energy dfensity analogy methods of planar heliostats in tower type solar thermo-power station, belong to tower type solar thermo-power station analogue technique field, comprising the following steps: (1) according to heliostat reflection light direction, define imaging plane;(2) by the uniform subdivision of three-dimensional space, there are the heliostats of shade and hiding relation with heliostat for calculating;(3) rasterisation of imaging plane indicates;(4) the radiant energy dfensity distribution on imaging plane is calculated;It (5), will be on the radiant energy dfensity distribution projection to receiving plane on imaging plane by oblique parallel projection.Solution procedure is accelerated using the correlation technique of graphics, and is realized based on GPU, compared to analytic method before, there is high-efficient, accurate advantage.
Description
Technical field
The present invention relates to tower type solar thermo-power station analogue technique fields, specifically, being related to a kind of tower-type solar thermal
The radiant energy dfensity analogy method of planar heliostats in power station.
Background technique
Effective use solar energy is an important channel for solving energy shortage and problem of environmental pollution.In the benefit of solar energy
In, tower type solar system is one important directions of photo-thermal power generation.Mainly have two for tower type solar system emulation at present
Class method: ray trace and analytic modell analytical model (Garcia, P., Ferriere, A., Bezian, J.J., 2008.Codes for
solar flux calculation dedicated to central receiver system applications:a
comparative review.Solar Energy 82,189–197.)
Ray tracing method can be naturally by heliostat shape, receiver shape, solar model and tracking error etc.
Factor is taken into account, and an accurate simulation result is obtained.Although ray trace models more tower type solar system
It is accurate and flexible to add, but since ray trace computational efficiency is low and the unstable disadvantage of peak value, and it is not suitable for the extensive settled date
Jing Chang emulation.
Analytic modell analytical model is indicated the radiant energy dfensity distribution on receiver using mathematical formulae, is suitable for mirror field cloth
The tasks such as office's optimization, focusing strategy optimization and mirror field radiant energy dfensity distributed simulation.Classical analytic modell analytical model includes HFLCAL
With UNIZAR model.
HFLCAL be heliostat field layout and optimization program (P., Pitz-Paal R.,
Schmitz M.Visual HFLCAL—A Software Tool for Layout and Optimisation of
Heliostat Fields [C] //Proceedings of 15th International SolarPACES Symposium,
Berlin, September.2009:15-18.).HFLCAL carves radiant energy dfensity distribution using mathematical model is simplified
It draws, the radiant energy dfensity distribution of heliostat hot spot is described as isotropism Gaussian Profile.HFLCAL model computational efficiency is high, and one
As be used for tower type solar system high-speed simulation.Since HFLCAL assumes the circular symmetry of radiant energy dfensity distribution, HFLCAL
Result and inaccurate.
UNIZAR model indicates the radiation energy on receiver surface using the convolution in gaussian kernel function and effective reflection area domain
Density Distribution (Collado F., Gomez A., Turegano J.An analytic function for the flux
Density due to sunlight reflected from a heliostat [J] .Solar Energy, 1986,37
(3): 215-234.).Alberto improves the versatility of UNIZAR model through projection, defines an imaging first
Plane solves the radiant energy dfensity distribution obtained on imaging plane, finally by (S á on energy distribution projection to receiving plane later
Nchez-Gonz á lez, Alberto, Santana D.Solar flux distribution on central
Receivers:A projection method from analytic function [J] .Renewable Energy,
2015,74:576-587.).
In tower type solar system emulation, the parameter of HFLCAL and UNIZAR model is about heliostat position and time
Function.The model parameter of the heliostat of different location is different, and the same face heliostat can also become in different moments, model parameter
Change.Callado proposes the method for one-site fit, by be fitted hot spot peak value energy, come determine HFLCAL model and
Parameter (the Collado F.J.One-point fitting of the flux density produced of UNIZAR model
by a heliostat[J].Solar Energy,2010,84(4):673-684.).Cruz proposes the method for an interpolation
The model parameter of heliostat is solved, 30 faces in 541 face heliostats has been used to feature the model ginseng of entire Jing Chang
Number (Cruz N.C.,J.D.,Redondo J.L.,et al.A new methodology for building-
up a robust model for heliostat field flux characterization[J].Energies,2017,
10(5):730.).But this method is analyzed just for the place data at a moment, for different moments, is needed again
It is modeled.
The defect of analytic modell analytical model essentially consists at present: precision is low;Parameter is difficult to determine, needs to pass through peak in practical application
The mode of value fitting and interpolation carries out parameter model.
Summary of the invention
It is an object of the present invention to provide a kind of radiant energy dfensity simulation sides of planar heliostats in tower type solar thermo-power station
Method, this method can be used for the radiant energy dfensity emulation of extensive tower type solar system, and utilize the calculation of computer graphics
Method and technology realize acceleration, have high-efficient, accurate and versatile advantage.
To achieve the goals above, in tower type solar thermo-power station provided by the invention planar heliostats radiant energy dfensity
Analogy method the following steps are included:
(1) according to the geographical location of tower type solar thermo-power station and simulation time, sunlight incident direction is calculated, and is adjusted
The direction of heliostat in thermo-power station defines a virtual image plane according to the direction of heliostat reflection light;
(2) it is a series of equally distributed three-dimensional rectangular parallelepiped grids by the uniform subdivision of three-dimensional space, each grid is known as body
Element saves the corresponding relationship of heliostat and voxel;Incidence and reflection direction along light calculate the voxel intersected with light beam,
And the corresponding heliostat of voxel is saved in shade and is blocked in heliostat set;
Preferably, friendship is asked for the light beam and scene of heliostat reflection, four light emitting by heliostat vertex
Pair it replaces light beam and scene carries out that friendship is asked to judge, and then judges that there are the heliostats of shade and hiding relation with heliostat, i.e.,
The light of four vertex of heliostat transmitting calculates the heliostat there are shade and hiding relation by 3D-DDA respectively, and for
As a result union is sought.
(3) imaging plane by heliostat reflection light direction definition heliostat on receiver surface, to imaging plane into
Line rasterization indicates, by the discrete effective reflection area domain for pixel expression heliostat on imaging plane of imaging plane;
Preferably, after obtaining rasterisation expression of effective view field on imaging plane, need to increase an area
Makeover process:
Wherein:
rcorrectIndicate the amendment ratio of discrete results, i.e., the discrete pixel values on imaging plane are multiplied by the amendment ratio;
SgeometryIndicate the preferred view area in effective reflection area domain;
SdisperseThe uncorrected gross area for indicating the discrete representation of effective reflection area domain is asked by the calculating of GPU specification
Solution.
(4) radiant energy dfensity distribution is calculated on imaging plane.The FDD of planar heliostats hot spot is the effective anti-of heliostat
The convolution of the FDD of region and mirror surface infinitesimal hot spot is penetrated, wherein the FDD of single mirror surface infinitesimal hot spot is considered as convolution kernel function,
The effective reflection area domain of mirror surface is convolution function, and wherein the calculating of discrete convolution has used the fast Flourier convolution of GPU version
It calculates:
Wherein:
FimageIndicate the radiant energy dfensity distribution on imaging plane;
B indicates that the rasterisation in mirror surface effective reflection area domain on imaging plane of heliostat indicates;
C is the radiant energy dfensity distribution of heliostat infinitesimal hot spot;
It (5), will be on the radiant energy dfensity distribution projection to receiving plane on imaging plane by oblique parallel projection:
Freceiver(u, v)=Fimage(fu(u,v),fv(u,v))cosθ
Wherein:
FreceiverIndicate the radiant energy dfensity distribution on receiving plane;
(u, v) is on receiving plane a little in the coordinate of receiving plane local coordinate system oUV plane;
(fu(u,v),fv(u, v)) it is some oU ' V ' plane coordinates under imaging plane local coordinate system on imaging plane,
Point (u, v) and point (fu(u,v),fv(u, v)) corresponding relationship be oblique parallel projection;
The angle of θ expression reflection light and receiving plane normal direction.
Preferably, the step on the radiant energy dfensity distribution projection to receiving plane of imaging plane is as follows:
(5-1) obtains coordinate of the discrete pixel in world coordinate system on receiving plane;
(5-2) obtains on imaging plane corresponding points in the coordinate of world coordinate system by oblique parallel projection;
(5-3) acquires position of the corresponding points in imaging plane by coordinate of the corresponding points in world coordinate system;
The radiant energy dfensity of any is the radiant energy dfensity of corresponding points on imaging plane multiplied by receiving plane on (5-4) receiving plane
Cosine, wherein receiving plane cosine is the cosine value of reflection light and receiving plane normal angle.
Preferably, it when calculating the radiant energy dfensity distribution of corresponding points on imaging plane, is calculated using bilinear interpolation pair
The energy value that should be put indicates the radiant energy dfensity of the point by the weighted average of sampled point adjacent pixel radiant energy dfensity.
Compared with prior art, the invention has the benefit that
(1) present invention replaces light beam and scene to carry out that friendship is asked to judge using four light that heliostat vertex emits, and
Uniform grid subdivision is carried out for scene, quickly calculates with light there are shade and blocks using the 3D-DDA algorithm in graphics
The heliostat of relationship;
(2) present invention uses the effective reflection area domain of the rasterization operation Efficient Solution heliostat in rendering pipeline, and
The process of area correction is increased after rasterisation indicates, and using the convolution kernel function of fitting, improves the precision of emulation;
(3) operations such as effective view field's calculating of the invention, convolutional calculation, oblique parallel projection realize have based on GPU
High-efficient advantage.
Detailed description of the invention
It is the heliostat field for verifying the accuracy design of the model proposed that Fig. 1, which is in the embodiment of the present invention,;
Fig. 2 is the definition of imaging plane and receiving plane and imaging plane local coordinate system in the embodiment of the present invention;
Fig. 3 (a) is the floor map that heliostat field ground grid divides in the embodiment of the present invention, and Fig. 3 (b) is that the present invention is real
It applies in example for storing the array of heliostat Yu spatial voxel corresponding relationship;
Fig. 4 is the schematic diagram that effective reflection area domain calculates in the embodiment of the present invention;
Fig. 5 is the flow chart that effective reflection area domain calculates in the embodiment of the present invention;
Fig. 6 is the exemplary diagram of convolutional calculation on imaging plane in the embodiment of the present invention;
Fig. 7 is HFLCAL model (first row) and UNIZAR model (secondary series) He Benfa in the Spring Equinox day moment at noon example
The contour comparison diagram of method in bright embodiment, every a line respectively correspond 0,4,20, No. 24 heliostat in Fig. 1;
Fig. 8 is HFLCAL model (first row) and UNIZAR model (secondary series) He Benfa in the Spring Equinox day moment at noon example
The contour comparison diagram of method in bright embodiment, every a line respectively correspond 7,9,37, No. 39 heliostats in Fig. 1.
Specific embodiment
To make the object, technical solutions and advantages of the present invention clearer, with reference to embodiments and its attached drawing is to this hair
It is bright to be described further.
Embodiment
Referring to Fig. 1 to Fig. 8, by taking an experiment heliostat field as an example, the heliostats distribution of Jing Chang is as shown in Figure 1.Following table 1
It is the design parameter in place.
Table 1 tests mirror field relevant parameter
The radiant energy dfensity analogy method of planar heliostats includes following step in the tower type solar thermo-power station of the present embodiment
It is rapid:
(1) according to heliostat reflection light direction, imaging plane is defined
Imaging plane was receiving plane center, and normal direction is a virtual plane of sunlight exit direction.Definition imaging is flat
The main reason for face, is: heliostat infinitesimal hot spot can be approximated to be on imaging plane circle it is symmetrical, convenient for calculate.
Fig. 2 indicates the corresponding relationship between imaging plane and plane receiving plane, and UVW is the local coordinate system of receiving plane, sits
Mark the center that origin is receiving plane, V is located on the intersection of receiving plane and horizontal plane, U be located at it is on receiving plane and vertical with V, W
In the normal orientation of receiving plane.U'V'W' is the local coordinate system of imaging plane, and V' is located at the intersection of imaging plane and horizontal plane
On, U' is located on imaging plane and vertical with V', and W' is located at the normal direction of imaging plane, is parallel to the side of the reflection light of heliostat
To.
(2) by the uniform subdivision of three-dimensional space, there are the heliostats of shade and hiding relation with heliostat for calculating
The present embodiment divides place using regular grid.As shown in figure 3, Fig. 3 (a) is the plane that Jing Chang is divided
Schematic diagram, Jing Chang are divided into 6 voxels, respectively number 0-5;There are 5 face heliostats in mirror field, number I-V respectively.Due to fixed
There is the operation such as rotation and need to consider the bounding box of heliostat when calculating heliostat and voxel overlapping relation in solar eyepiece.Single body
Element may correspond to 0,1 or multiple heliostats, and a heliostat may also be in multiple voxels.It is fixed as shown in Fig. 3 (b)
The corresponding relationship of solar eyepiece and voxel is stored in array of indexes.Consider GPU efficient realization, index relative be stored in two it is one-dimensional
In array.Wherein, Match array stores the corresponding heliostat index of all voxels in order, and Index array stores corresponding voxel
Initial position of the index data in Match array.Initial position of the voxel i that Index [i] is indicated in Match array, body
Plain i is that the heliostat of Index [i] to Index [i+1] -1 intersect, such as index is in voxel 2 and Match with index in Match
Heliostat intersects in 2 to 3, and Match [2] is that 2, Match [3] is 3, therefore voxel 2 intersects with heliostat 2 and heliostat 3.
The present embodiment replaces light beam and scene to carry out that friendship is asked to judge by four light that heliostat vertex emits, in turn
There are the heliostats of shade and hiding relation with heliostat for judgement, i.e., pass through respectively to the light of four vertex of heliostat transmitting
There are shade and the heliostats that block calculating by 3D-DDA, and seek union for result.
(3) rasterisation of imaging plane indicates
As shown in figure 4, solving H2At the effective reflection area domain on imaging plane, heliostat is calculated first and is being imaged
Projection H ' in plane2.The view field of heliostat is there may be shade and blocks, and obtains H2After projection on imaging plane,
Needing calculating, there are the projections of shade and the heliostat of hiding relation.For there are the heliostat H of shade relationship1, need along
Incident ray telegoniometer calculates it in the projection H ' of settled date mirror plane1, then by H '1From heliostat plane projection to imaging plane
On H "1.For there may be the heliostat H of hiding relation3, need heliostat H3Along emergent ray direction projection at
As the H ' in plane3。
For heliostat H2The calculating in effective reflection area domain, calculation process are as shown in Figure 5.H2Effective reflection area domain solve
It is largely divided into following two step:
(3-1) calculates heliostat H on imaging plane2View field H '2, judge each pixel center coordinate with
H′2Relationship, for the point in region, the value of respective pixel is set as 1, is otherwise provided as 0;
(3-2) is for there may be the heliostat H of shade and hiding relation1And H3, the throwing on imaging plane is found out respectively
Shadow H "1With H '3, 0 is set by the pixel value in corresponding region.
After obtaining rasterisation expression of effective view field on imaging plane, need to increase to an area correction mistake
Journey:
Wherein:
rcorrectIndicate the amendment ratio of discrete results, i.e., the discrete pixel values on imaging plane are multiplied by the amendment ratio;
SgeometryIndicate the preferred view area in effective reflection area domain;
SdisperseThe uncorrected gross area for indicating the discrete representation of effective reflection area domain is asked by the calculating of GPU specification
Solution.
(4) the radiant energy dfensity distribution on imaging plane is calculated
The convolutional calculation of radiant energy dfensity distribution is defined on imaging plane, and the FDD on imaging plane is heliostat infinitesimal
The convolution of hot spot FDD and effective reflection area domain.As shown in fig. 6, dotted line internal representation heliostat reflects the light on imaging plane
Spot;Solid black lines indicate effective reflection area domain of the heliostat on imaging plane, and view field is by discrete for mirror surface infinitesimal, circle
Region indicates that the hot spot of heliostat infinitesimal, final energy hot spot are the cumulative of heliostat infinitesimal hot spot, i.e. effective reflection area domain
With the convolution of heliostat infinitesimal hot spot.
It (5), will be on the radiant energy dfensity distribution projection to receiving plane on imaging plane by oblique parallel projection
Under world coordinate system, p indicates on imaging plane that a bit, respective coordinates are (px,py,pz);On p ' expression receiving plane
A bit, respective coordinates are (p 'x,p′y,p′z);R=(rx,ry,rZ) indicate the direction of oblique parallel projection, i.e. the outgoing side of light
To, defined by imaging plane it is found that r be imaging plane normal direction;The coordinate origin of o expression imaging plane and receiving plane.In
When calculating projection, need to regard sunlight as directional light.According to light along straightline propagation, the relationship between p and p ' can be obtained:
P+tr=p '
Wherein, t indicates the distance between p and p ', and the calculation of t is as follows:
T=r (p '-o)
In summary two formula can obtain p and p ' corresponding relationship:
P=p '+r (o-p ') r
The coordinate of p ' oUV plane under receiving plane local coordinate system is (ur,vr), p is under imaging plane local coordinate system
The coordinate of oU ' V ' plane is (ui,vi).When solving the corresponding relationship of two o'clock, it is necessary first to which local coordinate system is mapped to the world
In coordinate system.
P '=o+urU+vrV
P=o+uiU′+viV′
(ui,vi) and (ur,vr) two o'clock corresponding relationship it is as follows:
FreceiverIndicate the FDD, F on receiving planeimageIndicate the FDD on imaging plane.According to above-mentioned derivation result,
The FDD distribution on receiving plane can be obtained:
Freceiver(u, v)=Fimage(fu(u,v),fv(u,v))cosθ
Wherein, θ indicates the angle of emergent ray and receiver surface normal.
After obtaining the coordinate correspondence relationship of receiving plane and imaging plane by oblique parallel projection, the respective pixel that is calculated
Coordinate is decimal, not specific pixel coordinate, this example carries out bilinear interpolation using the energy value of sampled point adjacent pixel
Obtain the energy value of the point.
Evaluation index:
The evaluation index of radiant energy dfensity distribution modeling mainly has: gross energy error, peak error, mean square error.
In the present embodiment, legitimate reading is the average result that ray trace carries out 100 experiments, wherein heliostat each time
The light number of infinitesimal transmitting is 1024, and the granularity of division of heliostat is 0.01m.
The error of gross energy error gross energy between fitting result and legitimate reading:
Wherein, Ert_totalFor the gross energy of ray trace result, EtotalFor the gross energy of fitting result.
The error of peak error peak value between fitting result and legitimate reading:
Wherein, Crt_peakFor the peak value of ray trace result, CpeakFor the peak value of fitting result.
The difference of mean square error expression fitting result and legitimate reading energy density distribution:
Wherein, M, N indicate the resolution ratio of receiving plane, SsubIndicate the elemental area on receiver surface, i.e., each discrete unit
Area, C (m, n) indicates radiant energy density of the fitting result at (m, n), and Gt (m, n) indicates light result at (m, n)
Radiant energy density.
For the σ in the gaussian kernel function and s in quasi- Cauchy's kernel function, the present embodiment by minimize loss function come
It is solved, loss function is defined as follows:
Loss=α (| Cpeak-Crt_peak|)+(1-α)Errorrms
Wherein, Loss () is the value of loss function, CpeakIndicate the peak value of fitting kernel function, Crt_peakIndicate light with
The peak value of track result, ErrorrmsIndicate the mean square error of fitting kernel function and ray trace result.α is for weighing peak error
With the importance of mean square error, for convolution kernel function, peak error and fitting mean square error are for the accuracy of model
It has a major impact, therefore value 0.5 in the present embodiment α.
Experimental result:
Fig. 7 is the hot spot fitting result of short distance heliostat, and the horizontal distance of heliostat and receiver is 100m to 500m
Between, wherein each column respectively correspond HFLCAL, UNIZAR and planar heliostats model as a result, the unit of reference axis is m, energy
The unit of metric density is W/m2.From figure 7 it can be seen that fitting effect of the heliostat apart from receiving plane HFLCAL when closer is bad,
The hot spot of heliostat #0 is approximately rectangle, and hot spot is fitted ovalisation by HFLCAL, and has ignored the shape of heliostat, hot spot shape
Shape differs greatly, and there is also the problems of form fit inaccuracy for the hot spot of other heliostats.The fitting result and this reality of UNIZAR
It is close to apply a method, the planar heliostats model generally in the present embodiment is more nearly with legitimate reading, such as heliostat #4
Hot spot be ellipse, the fitting result of UNIZAR is more biased towards in a circular distribution.
For the hot spot of remote heliostat, the hot spot fitting result of planar heliostats model be substantially better than HFLCAL and
UNIZAR.Fig. 8 is the fitting result of remote heliostat hot spot, and the horizontal distance of heliostat and receiver is 500m to 1000m
Between, wherein each column respectively correspond the result of HFLCAL, UNIZAR and planar heliostats model.The light of planar heliostats model
Spot is almost the same with legitimate reading, and the result diffusion of HFLCAL and UNIZAR are obvious, such as the fitting knot of heliostat #39
The contour of fruit, HFLCAL and UNZIAR and true value gap are larger, and the result of planar heliostats model is opposite with legitimate reading
Unanimously.
Following table 2 is the mean error that HFLCAL, UNIZAR and planar heliostats model emulate Jing Chang, is Jing Chang
In all heliostats phantom error percentage average value.For peak error and gross energy error, planar heliostats model
Peak total energy error be better than HFLCAL and UNIZAR, the average peak error of UNIZAR model is better than planar heliostats model
And HFLCAL.For fitting mean square error and loss function, planar heliostats model in the present embodiment better than HFLCAL with
UNIZAR model.
2 Jing Chang of table emulates mean error
Claims (8)
1. the radiant energy dfensity analogy method of planar heliostats in a kind of tower type solar thermo-power station, which is characterized in that including with
Lower step:
(1) according to the geographical location of tower type solar thermo-power station and simulation time, sunlight incident direction is calculated, and adjusts thermoelectricity
The direction of heliostat in standing defines a virtual image plane according to the direction of heliostat reflection light;
(2) it is a series of equally distributed three-dimensional rectangular parallelepiped grids by the uniform subdivision of three-dimensional space, each grid is known as voxel, is protected
Deposit the corresponding relationship of heliostat and voxel;Incidence and reflection direction along light, calculate the voxel that intersects with light beam, and by body
The corresponding heliostat of element is saved in shade and blocks in heliostat set;
(3) imaging plane by heliostat reflection light direction definition heliostat on receiver surface carries out light to imaging plane
Gated indicates, by the discrete effective reflection area domain for pixel expression heliostat on imaging plane of imaging plane;
(4) radiant energy dfensity distribution is calculated on imaging plane:
Wherein:
FimageIndicate the radiant energy dfensity distribution on imaging plane;
B indicates that the rasterisation in mirror surface effective reflection area domain on imaging plane of heliostat indicates;
C is the radiant energy dfensity distribution of heliostat infinitesimal hot spot;
It (5), will be on the radiant energy dfensity distribution projection to receiving plane on imaging plane by oblique parallel projection:
Freceiver(u, v)=Fimage(fu(u,v),fv(u,v))cosθ
Wherein:
FreceiverIndicate the radiant energy dfensity distribution on receiving plane;
(u, v) is on receiving plane a little in the coordinate of receiving plane local coordinate system oUV plane;
(fu(u,v),fv(u, v)) be on imaging plane a little under imaging plane local coordinate system oU ' V ' plane coordinates, point (u,
And point (f v)u(u,v),fv(u, v)) corresponding relationship be oblique parallel projection;
The angle of θ expression reflection light and receiving plane normal direction.
2. the radiant energy dfensity analogy method of planar heliostats in tower type solar thermo-power station according to claim 1,
It is characterized in that, the heliostat is planar heliostats.
3. the radiant energy dfensity analogy method of planar heliostats in tower type solar thermo-power station according to claim 1,
It is characterized in that, in step (2), four light emitted by heliostat vertex replace light beam and scene to carry out asking friendship, and then sentence
Disconnected there are the heliostats of shade and hiding relation with heliostat, i.e., pass through 3D- respectively to the light of four vertex of heliostat transmitting
There are shade and the heliostats that block calculating by DDA, and seek union for result.
4. the radiant energy dfensity analogy method of planar heliostats in tower type solar thermo-power station according to claim 1,
It is characterized in that in step (3), indicating whether certain pixel is effective view field of the heliostat on imaging plane using 0 or 1.
5. the radiant energy dfensity analogy method of planar heliostats in tower type solar thermo-power station according to claim 1,
It is characterized in that, in step (3), after obtaining rasterisation expression of effective view field on imaging plane, increases an area and repair
Positive process:
Wherein:
rcorrectIndicate the correction value of discrete results, i.e., the discrete pixel values on imaging plane are multiplied by the correction value;
SgeometryIndicate the preferred view area in effective reflection area domain;
SdisperseThe uncorrected gross area for indicating the discrete representation of effective reflection area domain is solved by the calculating of GPU specification.
6. the radiant energy dfensity analogy method of planar heliostats in tower type solar thermo-power station according to claim 1,
It is characterized in that, in step (4), uses the radiant energy dfensity of the fast Flourier convolutional calculation heliostat infinitesimal hot spot of GPU version
The discrete values convolution of distribution and effective view field.
7. the radiant energy dfensity analogy method of planar heliostats in tower type solar thermo-power station according to claim 1,
It is characterized in that, in step (5), the step on the radiant energy dfensity distribution projection to receiving plane of imaging plane is as follows:
(5-1) obtains coordinate of the discrete pixel in world coordinate system on receiving plane;
(5-2) obtains on imaging plane corresponding points in the coordinate of world coordinate system by oblique parallel projection;
(5-3) acquires position of the corresponding points in imaging plane by coordinate of the corresponding points in world coordinate system;
The radiant energy dfensity of any is the radiant energy dfensity of corresponding points on imaging plane multiplied by more than receiving plane on (5-4) receiving plane
String, wherein receiving plane cosine is the cosine value of reflection light and receiving plane normal angle.
8. the radiant energy dfensity analogy method of planar heliostats in tower type solar thermo-power station according to claim 7,
It is characterized in that, is calculated as in step (5-4) as when the radiant energy dfensity distribution of corresponding points, being come using bilinear interpolation in plane
The energy value of corresponding points is calculated, i.e., indicates the radiation energy of the point by the weighted average of sampled point adjacent pixel radiant energy dfensity
Density.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910577119.8A CN110414059B (en) | 2019-06-28 | 2019-06-28 | Radiation energy density simulation method of planar heliostat in tower type solar thermal power station |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910577119.8A CN110414059B (en) | 2019-06-28 | 2019-06-28 | Radiation energy density simulation method of planar heliostat in tower type solar thermal power station |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110414059A true CN110414059A (en) | 2019-11-05 |
CN110414059B CN110414059B (en) | 2020-12-01 |
Family
ID=68358510
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910577119.8A Active CN110414059B (en) | 2019-06-28 | 2019-06-28 | Radiation energy density simulation method of planar heliostat in tower type solar thermal power station |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110414059B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111625957A (en) * | 2020-05-25 | 2020-09-04 | 浙江大学 | Radiation energy density distribution simulation method for planar light spot of tower type solar mirror field receiver |
Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101329112A (en) * | 2008-08-04 | 2008-12-24 | 北京理工大学 | Central receiving type high temperature solar heat collecting system of light funnel light gathering orienting light transmission |
US20100006087A1 (en) * | 2008-07-10 | 2010-01-14 | Brightsource Industries (Israel) Ltd. | Systems and methods for control of a solar power tower using infrared thermography |
CN101806502A (en) * | 2010-04-14 | 2010-08-18 | 俞琦 | Solar energy collecting system of heliostat |
CN102331793A (en) * | 2011-07-19 | 2012-01-25 | 南京科远自动化集团股份有限公司 | Method for scheduling heliostat field in tower solar thermal power station |
US20120160324A1 (en) * | 2009-09-02 | 2012-06-28 | 3M Innovative Properties Company | Concentrating solar mirror panel assembly with corrugated stiffener |
CN103250010A (en) * | 2010-08-20 | 2013-08-14 | 菲利普·施拉梅克 | Central receiver solar system comprising a heliostat field |
US20130239952A1 (en) * | 2012-03-14 | 2013-09-19 | Brightsource Industries (Israel) Ltd. | Methods and systems for operating solar tower systems |
US20130327400A1 (en) * | 2011-02-11 | 2013-12-12 | Jaime Caselles Fornés | Direct solar-radiation collection and concentration element and panel |
CN103530518A (en) * | 2013-10-14 | 2014-01-22 | 中国科学院电工研究所 | Computing method of lighting surface light condensation energy flux density distribution of solar tower type electric generation thermal absorber |
CN105973505A (en) * | 2016-05-10 | 2016-09-28 | 西安交通大学 | Method for determining heat-flow density in opening of solar cavity type heat absorber |
CN106096085A (en) * | 2016-05-30 | 2016-11-09 | 浙江大学 | A kind of radiant energy dfensity analogy method of tower-type solar thermal power generating system |
CN106225261A (en) * | 2016-08-16 | 2016-12-14 | 华电电力科学研究院 | A kind of solar power tower receiver architecture and method of reseptance |
CN107704711A (en) * | 2017-10-30 | 2018-02-16 | 中国华能集团清洁能源技术研究院有限公司 | A kind of tower type solar mirror field shade and the innovatory algorithm for blocking efficiency |
CN108829906A (en) * | 2017-12-29 | 2018-11-16 | 西安电子科技大学 | A kind of rescattering calculation method of complex target to background light radiation |
CN109670248A (en) * | 2018-12-21 | 2019-04-23 | 浙江大学 | The analogy method of the radiant energy dfensity distribution of hot spot in tower type solar energy thermal power generation |
CN109690208A (en) * | 2016-06-24 | 2019-04-26 | Sfi株式会社 | Heliostat device and method for electric generation using solar energy |
-
2019
- 2019-06-28 CN CN201910577119.8A patent/CN110414059B/en active Active
Patent Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100006087A1 (en) * | 2008-07-10 | 2010-01-14 | Brightsource Industries (Israel) Ltd. | Systems and methods for control of a solar power tower using infrared thermography |
CN101329112A (en) * | 2008-08-04 | 2008-12-24 | 北京理工大学 | Central receiving type high temperature solar heat collecting system of light funnel light gathering orienting light transmission |
US20120160324A1 (en) * | 2009-09-02 | 2012-06-28 | 3M Innovative Properties Company | Concentrating solar mirror panel assembly with corrugated stiffener |
CN101806502A (en) * | 2010-04-14 | 2010-08-18 | 俞琦 | Solar energy collecting system of heliostat |
CN103250010A (en) * | 2010-08-20 | 2013-08-14 | 菲利普·施拉梅克 | Central receiver solar system comprising a heliostat field |
US20130327400A1 (en) * | 2011-02-11 | 2013-12-12 | Jaime Caselles Fornés | Direct solar-radiation collection and concentration element and panel |
CN102331793A (en) * | 2011-07-19 | 2012-01-25 | 南京科远自动化集团股份有限公司 | Method for scheduling heliostat field in tower solar thermal power station |
US20130239952A1 (en) * | 2012-03-14 | 2013-09-19 | Brightsource Industries (Israel) Ltd. | Methods and systems for operating solar tower systems |
CN103530518A (en) * | 2013-10-14 | 2014-01-22 | 中国科学院电工研究所 | Computing method of lighting surface light condensation energy flux density distribution of solar tower type electric generation thermal absorber |
CN105973505A (en) * | 2016-05-10 | 2016-09-28 | 西安交通大学 | Method for determining heat-flow density in opening of solar cavity type heat absorber |
CN106096085A (en) * | 2016-05-30 | 2016-11-09 | 浙江大学 | A kind of radiant energy dfensity analogy method of tower-type solar thermal power generating system |
CN109690208A (en) * | 2016-06-24 | 2019-04-26 | Sfi株式会社 | Heliostat device and method for electric generation using solar energy |
CN106225261A (en) * | 2016-08-16 | 2016-12-14 | 华电电力科学研究院 | A kind of solar power tower receiver architecture and method of reseptance |
CN107704711A (en) * | 2017-10-30 | 2018-02-16 | 中国华能集团清洁能源技术研究院有限公司 | A kind of tower type solar mirror field shade and the innovatory algorithm for blocking efficiency |
CN108829906A (en) * | 2017-12-29 | 2018-11-16 | 西安电子科技大学 | A kind of rescattering calculation method of complex target to background light radiation |
CN109670248A (en) * | 2018-12-21 | 2019-04-23 | 浙江大学 | The analogy method of the radiant energy dfensity distribution of hot spot in tower type solar energy thermal power generation |
Non-Patent Citations (7)
Title |
---|
CAITOU HE .ETAL: "Fast flux density distribution simulation of central receiver system on GPU", 《SOLAR ENERGY》 * |
F.J. COLLAD .ETAL: "An analytic function for the flux density due to sunlight reflected from a heliostat", 《SOLAR ENERGY》 * |
FRANCISCO J. COLLADO: "One-point fitting of the flux density produced by a heliostat", 《SOLAR ENERGY》 * |
PIERRE GARCIA .ETAL: "Codes for solar flux calculation dedicated to central receiver system applications a comparative review", 《SOLAR ENERGY》 * |
WEIDONGHUANG .ETAL: "An analytical solution for the solar flux density produced by a round focusing heliostat", 《RENEWABLE ENERGY》 * |
徐二树 等: "塔式太阳能热电站系统仿真与分析", 《中国电机工程学报》 * |
盛玲霞 等: "塔式太阳能电站接收器的建模及动态仿真", 《化工学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111625957A (en) * | 2020-05-25 | 2020-09-04 | 浙江大学 | Radiation energy density distribution simulation method for planar light spot of tower type solar mirror field receiver |
CN111625957B (en) * | 2020-05-25 | 2022-06-10 | 浙江大学 | Radiation energy density distribution simulation method for tower type solar mirror field receiver plane |
Also Published As
Publication number | Publication date |
---|---|
CN110414059B (en) | 2020-12-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sánchez-González et al. | Solar flux distribution on central receivers: A projection method from analytic function | |
CN109670248B (en) | Simulation method for radiation energy density distribution of light spots in tower type solar thermal power generation | |
Prahl et al. | Airborne shape measurement of parabolic trough collector fields | |
Cardoso et al. | New functionalities for the Tonatiuh ray-tracing software | |
CN104034058B (en) | The formation method of tower type solar heat and power system Jing Chang based on GPU | |
CN110223389A (en) | The scene modeling method of blending image and laser data, system, device | |
García et al. | Spillage estimation in a heliostats field for solar field optimization | |
CN110276168A (en) | The tower non-central point focusing modeling method of photo-thermal power station heliostat field | |
Gebreiter et al. | sbpRAY–A fast and versatile tool for the simulation of large scale CSP plants | |
CN111856747B (en) | Optimal transmission-based reflector design method | |
He et al. | Fast flux density distribution simulation of central receiver system on GPU | |
Bonanos et al. | Heliostat surface shape characterization for accurate flux prediction | |
CN109697315B (en) | Optimization method of analytical model parameters of radiant energy light spots | |
Andraka et al. | Rapid reflective facet characterization using fringe reflection techniques | |
He et al. | An improved flux density distribution model for a flat heliostat (iHFLCAL) compared with HFLCAL | |
CN110414059A (en) | The radiant energy dfensity analogy method of planar heliostats in tower type solar thermo-power station | |
He et al. | An analytical flux density distribution model with a closed-form expression for a flat heliostat | |
CN111625957B (en) | Radiation energy density distribution simulation method for tower type solar mirror field receiver plane | |
Kaňuk et al. | Testing of V3. sun module prototype for solar radiation modelling on 3D objects with complex geometric structure | |
Richter et al. | Computationally fast analytical ray-tracer for central receiver systems | |
CN110414057A (en) | The radiant energy dfensity analogy method of focus type heliostat in tower type solar thermo-power station | |
Sproul et al. | The development of the heliostat focusing and canting enhancement technique: an optical heliostat alignment tool for the national solar thermal test facility | |
Arrif et al. | Shadowing and blocking factors in heliostats: Comparison between parallel and oblique projections | |
Hausner | Multipole expansion of the light vector | |
Owkes | An optical characterization technique for parabolic trough solar collectors using images of the absorber reflection |
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 |