CN110209207A - Determine the method and apparatus and machine readable storage medium of the lost area of heliostat - Google Patents
Determine the method and apparatus and machine readable storage medium of the lost area of heliostat Download PDFInfo
- Publication number
- CN110209207A CN110209207A CN201910375494.4A CN201910375494A CN110209207A CN 110209207 A CN110209207 A CN 110209207A CN 201910375494 A CN201910375494 A CN 201910375494A CN 110209207 A CN110209207 A CN 110209207A
- Authority
- CN
- China
- Prior art keywords
- heliostat
- area
- vertex
- coordinate system
- determining
- 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
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D3/00—Control of position or direction
- G05D3/12—Control of position or direction using feedback
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Photovoltaic Devices (AREA)
Abstract
The embodiment of the present invention provides the method and apparatus and machine readable storage medium of a kind of lost area of determining heliostat, belongs to thermal control process and fields of measurement.This method comprises: determining the corresponding loss region of the first heliostat, wherein loss region includes shadow region and/or occlusion area, heliostat in shadow region generates shade on the first heliostat, and the heliostat in occlusion area blocks the light reflected via the first heliostat；Determine at least one second heliostat in loss region；Determine at least one second heliostat be incident upon the first heliostat at least one projected area in the plane；Determine union refion of at least one projected area in the surface of the first heliostat；And determine the area of union refion, wherein the area of union refion is lost area.Add the area of the shade of heliostat and/or eclipsing loss under the control mode of the elevation angle in spin angle it is thereby achieved that determining.
Description
Technical field
The present invention relates to thermal control process and fields of measurement, more particularly to a kind of method of the lost area of determining heliostat
With device and machine readable storage medium.
Background technique
For the condenser system (heliostat field) of solar power tower, optical property, the optimization mirror of Jing Chang are improved
Field parameters are conducive to promote power station overall performance and reduce cost, are on the one hand because the cost of investment of condenser system accounts for entirely
The number of the 40%-50% of Jing Chang, another aspect heliostat field reflected sunlight directly determine that the energy of subsequent subsystem is defeated
Enter.
When emulating even power station simulation to condenser system, many times need to know that condenser system inputs to subsequent system
The energy of system, the energy are mainly determined by three factors: the mirror surface gross area, Intensity of the sunlight and mirror field efficiency.The above two compared with
It is easy to get, and mirror field efficiency is related to the optical property of every face heliostat, needs to consider again when seeking between surrounding mirror surface
Mutual alignment relation, thus it is more complicated, and especially when mirror surface quantity is more, calculation amount is often very huge.Specifically seeking
When one face heliostat efficiency, main method is geometric projection and Monte Carlo Ray-tracing Method.The former is using light as parallel
Light simplifies model, and since optical efficiency is made of a variety of component efficiencies, geometric method seeks each component efficiency respectively, then integrates
Obtain heliostat and the overall efficiency of Jing Chang.But due to the simplification of model, when calculating certain component efficiencies, precision is inadequate.For meter
Big problem is measured in calculation, it is thus proposed that mirror field blocking need to only calculate fritter Jing Chang one side or several times settled date by division methods
Mirror efficiency reduces calculation amount in this way, but to sacrifice precision as cost, especially studies different location in mirror field as representing
Heliostat efficiency change when, result error is larger.Ray-tracing Method is more flexible, and can simulate non-ideal optical element,
Its advantage, which also resides in reproduce between photon, simultaneously really interacts, to provide accurate result.But ray trace
Method needs a large amount of calculating time and high computing capability.In this regard, someone have studied the quick light based on computer architecture with
Track method.Somebody proposes the resolving ideas for combining two methods, and two kinds of sides are made full use of when calculating different component efficiencies
The advantage of method.
Jing Chang total light efficiency includes shade and blocks efficiency.It is fixed since Jing Chang arranges difference, and as the sun operates
The posture of solar eyepiece also changes.The gesture stability of heliostat is divided into two kinds at present, respectively the elevation angle+azimuth control mode,
Spin+elevation angle control mode.Due to the difference of control mode, shade and eclipsing loss calculation method be not also identical.The elevation angle+orientation
Shade and eclipsing loss calculating under the control mode of angle have been carried out research.But the yin under spin+elevation angle mode control mode
Shadow and eclipsing loss calculating are not studied.
Summary of the invention
The object of the present invention is to provide the method and apparatus and machine readable storage of a kind of lost area of determining heliostat
Medium can be realized and determine that being in spin angle adds the area of the shade of heliostat and/or eclipsing loss under the control mode of the elevation angle.
To achieve the goals above, one aspect of the present invention provides a kind of for determining the lost area of the first heliostat
Method, the mode for controlling first heliostat is that spin angle adds elevation angle control mode, this method comprises: determining described first
The corresponding loss region of heliostat is in the shadow region wherein the loss area domain includes shadow region and/or occlusion area
Heliostat in domain generates shade on first heliostat, and the heliostat in the occlusion area is blocked via described
The light of first heliostat reflection；Determine at least one second heliostat in the loss region；Described in determining at least
One the second heliostat be incident upon first heliostat at least one projected area in the plane；Determine described at least one
Union refion of a projected area in the surface of first heliostat；And determine the area of the union refion, wherein
The area of the union refion is the lost area.
Optionally it is determined that union refion of at least one the described projected area in the surface of first heliostat and
The area for determining the union refion includes: to draw the surface region of at least one described projected area and first heliostat
It is divided into grid；Lattice point on the grid divided is marked, wherein if a lattice point is at least one described projected area
Preset value is marked in the case where any one of the surface region of first heliostat；By each lattice point it is labeled described in
Preset value is summed to obtain the label summation of the lattice point, wherein in the surface region of first heliostat, label summation is big
In or equal to the region that occupies of twice of lattice point of the preset value be the union refion；And it is determined based on Pick theorem
Label summation in the surface region of first heliostat is greater than or equal to twice of lattice point of the preset value described the
The area in the region occupied in the surface region of one heliostat, wherein the area based on determined by Pick theorem be described in simultaneously
Collect the area in region.
Optionally, at least one second heliostat described in the determination be incident upon first heliostat in the plane
At least one projected area includes: the coordinate value that at least three vertex in first heliostat are determined in the first coordinate system
With the coordinate value on each vertex of every one second heliostat at least one described second heliostat, wherein described first sits
Mark system has the following characteristics that integrating thermal tower location face position as coordinate origin, due east direction is x-axis and direct north as y-axis；Base
Plane where the coordinate value at least three vertex of first heliostat determines first heliostat；Existed based on light
Every one second heliostat in the azimuth of first heliostat and elevation angle and at least one described second heliostat it is every
The coordinate value on one vertex determines straight line where each vertex；Determine every one second settled date at least one described second heliostat
The intersection point of straight line where each vertex of mirror and the plane where first heliostat；And based on it is described at least one second
The intersection point of each vertex correspondence of every one second heliostat in heliostat determines the corresponding projected area of every one second heliostat,
So that it is determined that at least one described projected area.
Optionally, the coordinate value that at least three vertex in first heliostat are determined in the first coordinate system and
The coordinate value on each vertex of every one second heliostat at least one described second heliostat comprises determining that described first is fixed
Coordinate value of the midpoint of solar eyepiece and each heliostat at least one described second heliostat in first coordinate system；?
In corresponding second coordinate system of each heliostat, the coordinate that the vertex of determining projected area is ready to use in each heliostat is determined
Value, wherein second coordinate system has the following characteristics that the center of heliostat for coordinate origin, the method at the center for crossing heliostat
Line direction is y-axis, rotary shaft corresponding with the elevation angle of first heliostat is x-axis；Spin based on first heliostat
Angle, the elevation angle, elevation angle and the transition matrix that second coordinate system is converted to first coordinate system is determined towards angle；And base
Being ready to use in really in identified transition matrix, each heliostat determined in corresponding second coordinate system of each heliostat
Coordinate value of the center of the coordinate value and each heliostat of determining the vertex of projected area in first coordinate system determines each
Coordinate value of the vertex for being ready to use in determining projected area in first coordinate system in heliostat.
Optionally, the transition matrix are as follows:
Wherein, ρ is the spin angle, and θ is the elevation angle, and λ is the elevation angle of first heliostat,For it is described towards
Angle.
Correspondingly, another aspect of the present invention provides a kind of for determining the device of the lost area of the first heliostat, controls
The mode for making first heliostat is that spin angle adds elevation angle control mode, which includes: loss area determination module, is used for
The corresponding loss region of first heliostat is determined, wherein the loss area domain includes shadow region and/or occlusion area, place
Shade is generated on first heliostat in the heliostat in the shadow region, the heliostat in the occlusion area
Block the light reflected via first heliostat；Heliostat determining module is lost, is in the loss region for determining
Second heliostat of at least one interior；Projected area determining module, for determining that at least one described second heliostat is incident upon
First heliostat at least one projected area in the plane；Union refion determining module, it is described at least for determination
Union refion of one projected area in the surface of first heliostat；And lost area determining module, for determining
The area of the union refion, wherein the area of the union refion is the lost area.
Optionally, the union refion determining module determines at least one described projected area in first heliostat
Union refion in surface includes: by the surface segmentation Cheng Fang of at least one described projected area and first heliostat
Lattice；Lattice point on the grid divided is marked, wherein if a lattice point is at least one described projected area and described
Preset value is marked in the case where any one of surface region of first heliostat；And each lattice point is labeled described pre-
If value is summed to obtain the label summation of the lattice point, wherein in the surface region of first heliostat, label summation is greater than
It or equal to the region that occupies of twice of lattice point of the preset value is the union refion；The lost area determining module determines
The area of the union refion includes: to determine the label summation in the surface region of first heliostat based on Pick theorem
More than or equal to the area in the region that twice of lattice point of the preset value occupies in the surface region of first heliostat,
Wherein, the area based on determined by Pick theorem is the area of the union refion.
Optionally, it is fixed to determine that at least one described second heliostat is incident upon described first for the projected area determining module
Solar eyepiece at least one projected area in the plane include: to be determined in first heliostat at least in the first coordinate system
The coordinate value on each vertex of every one second heliostat in the coordinate value and at least one described second heliostat on three vertex,
Wherein, first coordinate system have the following characteristics that integrate thermal tower location face position as coordinate origin, due east direction as x-axis and
Direct north is y-axis；The first heliostat place is determined based on the coordinate value at least three vertex of first heliostat
Plane；It is every in the azimuth of first heliostat and elevation angle and at least one described second heliostat based on light
The coordinate value on each vertex of one second heliostat determines straight line where each vertex；Determine at least one described second heliostat
In every one second heliostat each vertex where straight line and the plane where first heliostat intersection point；And it is based on
The intersection point of each vertex correspondence of every one second heliostat at least one described second heliostat determines every one second settled date
The corresponding projected area of mirror, so that it is determined that at least one described projected area.
Optionally, the coordinate value that at least three vertex in first heliostat are determined in the first coordinate system and
The coordinate value on each vertex of every one second heliostat at least one described second heliostat comprises determining that described first is fixed
Coordinate value of the midpoint of solar eyepiece and each heliostat at least one described second heliostat in first coordinate system；?
In corresponding second coordinate system of each heliostat, the coordinate that the vertex of determining projected area is ready to use in each heliostat is determined
Value, wherein second coordinate system has the following characteristics that the center of heliostat for coordinate origin, the method at the center for crossing heliostat
Line direction is y-axis, rotary shaft corresponding with the elevation angle of first heliostat is x-axis；Spin based on first heliostat
Angle and the elevation angle, elevation angle and the transition matrix that second coordinate system is converted to first coordinate system is determined towards angle；And
Being ready to use in each heliostat determined based on identified transition matrix, in corresponding second coordinate system of each heliostat
Coordinate value of the center of the coordinate value and each heliostat that determine the vertex of projected area in first coordinate system determines every
Coordinate value of the vertex for being ready to use in determining projected area in first coordinate system in one heliostat.
Optionally, the transition matrix are as follows:
Wherein, ρ is the spin angle, and θ is the elevation angle, and λ is the elevation angle of first heliostat,For it is described towards
Angle.
In addition, another aspect of the present invention provides a kind of machine readable storage medium, deposited on the machine readable storage medium
Instruction is contained, which is used for so that machine executes above-mentioned method.
Through the above technical solutions, being under spin angle adds elevation angle control mode, to determine in the mode for controlling the first heliostat
Other heliostats in loss region in the first heliostat the first heliostat at least one projected area in the plane
Domain determines union refion of at least one view field in the surface of the first heliostat, determines the area of union refion, wherein
The area of the union refion is the lost area of the first heliostat, and lost area includes the area of shade and/or eclipsing loss,
Add the area of the shade of heliostat and/or eclipsing loss under the control mode of the elevation angle in spin angle in this way, realizing and determining.
Other features and advantages of the present invention will the following detailed description will be given in the detailed implementation section.
Detailed description of the invention
Attached drawing is to further understand for providing to the embodiment of the present invention, and constitute part of specification, under
The specific embodiment in face is used to explain the present invention embodiment together, but does not constitute the limitation to the embodiment of the present invention.Attached
In figure:
Fig. 1 is the schematic diagram of heliostat gesture stability mode；
Fig. 2 is that spin angle adds elevation angle control mode schematic diagram；
Fig. 3 is that the spin angle of heliostat adds elevational movement pattern diagram；
Fig. 4 is the flow chart of the method for the lost area for determining the first heliostat that one embodiment of the invention provides；
Fig. 5 be another embodiment of the present invention provides the first heliostat receive the side view of shade；
Fig. 6 be another embodiment of the present invention provides the first heliostat receive shade top view；
Fig. 7 be another embodiment of the present invention provides determination shadow region range schematic diagram；
Fig. 8 be another embodiment of the present invention provides exemplary shadow shape schematic diagram；And
Fig. 9 be another embodiment of the present invention provides for determine the first heliostat lost area device structural frames
Figure.
Description of symbols
1 loss area determination module 2 loses heliostat determining module
3 projected area determining module, 4 union refion determining module
5 lost area determining modules
Specific embodiment
It is described in detail below in conjunction with specific embodiment of the attached drawing to the embodiment of the present invention.It should be understood that this
Locate described specific embodiment and be merely to illustrate and explain the present invention embodiment, is not intended to restrict the invention embodiment.
Jing Chang total light efficiency includes atmospheric transmissivity ηat, cosine efficiency etacosin, shade and block efficiency etaS&BAnd spilling
Efficiency etaint.Mirror field efficiency is calculated and analyzed by this system realization.
The optical efficiency of Jing Chang are as follows: ηfield=ηat×ηcosin×ηS&B×ηint
Atmospheric transmissivity ηat:
ηat=exp (- 0.106 × 10-4S0) S0> 1000, wherein S0For the distance of the position Jing Chang to target point.
Cosine efficiency etacosinThe good fortune that heliostat caused by tilting for heliostat is an actually-received is penetrated penetrates less than theoretical maximum good fortune
The phenomenon that: ηcosin=cos θ, wherein θ is acute angle formed by plane where heliostat and the plane vertical with sun incident light.
Shade and block efficiency etaS&BWhat is indicated is that shade does not occur for mirror surface or eclipsing loss " effective " area accounts for the gross area
Ratio situation.
Overflow efficiency etaintIt is that heat collector surface or entrance are not reached by the solar radiant energy of heliostat reflection, is spilled over to
Energy loss resulted in ambient atmosphere is known as overflowing loss.
Since Jing Chang arranges difference, and as the sun operates, the posture of heliostat also changes.Heliostat posture
Control mode can substantially be divided into two kinds: the elevation angle adds azimuth control mode, spin angle to add elevation angle control mode, and control mode is not
Together, shade is different with the calculation method for blocking efficiency.Wherein, spin angle adds elevation angle control mode as shown in Figure 1, wherein in Fig. 1
A figure indicates is that azimuth adds elevation angle control mode, what b in Fig. 1 figure indicated is that spin angle adds elevation angle control mode.Its
In, what spin angle added the control mode at the elevation angle is characterized mainly in that its spin axis is directed toward target heat collector always.
Elevation angle control mode is added to explain spin angle below with reference to Fig. 2.It should be noted that add to spin angle herein
The explanation of elevation angle control mode is only used for understanding that spin angle adds elevation angle control mode, is not intended to restrict the invention.
In Fig. 2, heliostat or its frame are represented with a plane, O is mirror surface central point, and vector OT is directed toward target spot (mesh
Mark heat collector), OS is directed toward the sun, and ON is the normal of mirror surface.The a figure of Fig. 2 is shown the sun and the same orientation of target spot and (for example all exists
Due South) the best time the case where, it is assumed that the posture of mirror and azimuth add elevation angle mode to have no difference at this time, OT, OS and
For ON all in the optical system meridian plane vertical with mirror surface, ON is the angular bisector of the angle of OT and OS, and mirror surface is about at this time
Meridian plane bilateral symmetry.The b figure of Fig. 2 shows the case where sun is in another orientation, and the mirror of new tracking mode is not around vertical
Azimuth axis in ground rotates, but rotates an angle around OT axis, which is indicated with ρ, and angle ρ is spin angle, and spin makes
Normal ON is gone in the plane determined by OT and new OS (OS '), becomes OT ', the plane of symmetry of mirror surface is still meridian at this time
Face.Assuming that the angle of OT and OS ' is 2 θ, mirror surface should also turn while spin around another axis parallel with mirror surface sagitta of arc line
Dynamic, to change the elevation angle of mirror surface, make ON ' and OS ' angle θ.θ is exactly the incidence angle of sunlight at this time, that is, the elevation angle.This
Outside, the direction of motion of spin angle ρ and elevation angle theta is referred to the arrow in Fig. 3 to understand.Wherein, in Fig. 3, around the arrow of T`T
For the direction of motion of spin angle ρ, FF` is the direction of motion of elevation angle theta.
The one aspect of the embodiment of the present invention provides a kind of method for determining the lost area of the first heliostat.Its
In, the mode for controlling the posture of the first heliostat is that spin angle adds elevation angle control mode.
Fig. 4 is the flow chart of the method for the lost area for determining the first heliostat that one embodiment of the invention provides.
As shown in figure 4, this method includes the following contents.
In step s 40, determine the corresponding loss region of the first heliostat, wherein loss region include shadow region and/
Or occlusion area, the heliostat in shadow region generate shade on the first heliostat, the settled date in occlusion area
Mirror blocks the light reflected via the first heliostat.
For example, it is assumed that sunray is parallel rays, the first heliostat is rectangular heliostat, is determined in heliostat field
The height of solar eyepiece is substantially identical.Sunray vertical incidence can determine in the first heliostat generates shade to the first heliostat
The maximum magnitude in region, wherein the farthest of other the first heliostats of heliostat distance of shade can be generated to the first heliostat
Distance may refer to shown in Fig. 5.Fig. 5 indicates that heliostat a (being equal to the first above-mentioned heliostat) receives the side view of shade, d
Indicate can to heliostat a generate shade other heliostats on the front-rear direction along heliostat a apart from the farthest of heliostat a
Distance.Wherein, it according to Fig. 5, can calculateL and W is respectively fixed
The length on two sides of solar eyepiece a, len are the half of the catercorner length of heliostat,In addition,
Other heliostats on the left and right directions of heliostat a can also generate shade to heliostat a.The corresponding shadow region heliostat a
It may refer to Fig. 6, wherein Fig. 6 is that heliostat a receives shade top view.In Fig. 6, A ', B ', C ', the 4 points of areas surrounded D '
Domain A is the corresponding shadow region heliostat a, can determine shadow region by determining A ', B ', C ', 4 points of D ' of seat
Specific range.
Below with reference to Fig. 7 by determine B ' point coordinate for, be illustrated, wherein point A ', C ', D ' coordinate determination
It may refer to determine the seat calibration method of point B '.Assuming that the coordinate of the central point M of heliostat a is (x, y, H), wherein heat collector
The center for integrating thermal tower as coordinate being mounted on, direct north are positive direction of the y-axis, and due east direction is positive direction of the x-axis, because, it is fixed
The height of heliostat in solar eyepiece field is identical, therefore, only determines coordinate of each point on x/y plane.B ' U is perpendicular to UM, ∠
B ' MU=β-270。
Referring to derive point B ' seat calibration method, it can be deduced that point A ', C ', D ' coordinate be respectively as follows:
A ':
C ':
D ':
For the occlusion area of the first heliostat, the method for determination is referred to shadow region.Determine in occlusion area
Solar eyepiece causes to block to the light that the first heliostat reflects, and the light blocked cannot be reflected into heat collector, therefore, can be according to light
The invertibity of line, it is assumed that reflection light is to be issued by heat collector, that is, heat collector is equivalent to the sun in determining shadow region,
And then the occlusion area of the first heliostat is determined referring to the method for above-mentioned determining shadow region.Shade is generated in this way, realizing and determining
The region of the heliostat of loss and/or eclipsing loss.
In step S41, at least one second heliostat in loss region is determined.It is alternatively possible to determine horizontal glass
The coordinate of heliostat in, judges whether the coordinate of heliostat is in loss region.
In step S42, determine at least one second heliostat be incident upon the first heliostat in the plane at least one
A projected area.Specifically, if at least one second heliostat is the heliostat for generating shade to the first heliostat, at least
One heliostat is projected when being projected to plane where the first heliostat along the direction of sunray；If at least one
Second heliostat is to generate the heliostat that blocks to the light of the first heliostat reflection, then, at least one heliostat is to first
Plane where heliostat is projected when being projected along the direction of reflection light；If at least one second heliostat includes difference
Heliostat in shadow region and the heliostat in occlusion area, then respectively along the direction of sunray and reflected light
Plane where from the direction of line to the first heliostat is projected.
In step S43, union refion of at least one projected area in the surface of the first heliostat is determined, that is,
The intersection area of the union refion of at least one projected area and the surface region of the first heliostat, at least one projected area
Union refion is in the part in the surface range of the first heliostat.
In step 44, the area of union refion of at least one projected area in the surface of the first heliostat is determined,
Wherein, the area of identified union refion is lost area.
It is under spin angle adds elevation angle control mode, to determine the damage for being in the first heliostat in the mode for controlling the first heliostat
Lose other heliostats in region the first heliostat at least one view field in the plane, determine at least one projection
Union refion of the region in the surface of the first heliostat, determines the area of union refion, and wherein the area of the union refion is
For the lost area of the first heliostat, lost area includes the area of shade and/or eclipsing loss, in this way, realizing determining place
The shade of heliostat and/or the area of eclipsing loss under spin angle plus elevation angle control mode.
In three dimensions, if the control mode of the posture of heliostat is that spin angle adds elevation angle control mode, heliostat
Motion conditions it is more complicated, then there are many kinds of issuable shadow shapes, as shown in Figure 8.It should be noted that Fig. 8 is only
Illustratively list several possible shadow shapes, it is possible to create shadow shapes be not limited in these types.It can according to Fig. 8
To find out, shadow shapes are not unified, and with the operating of heliostat, and shadow shapes are also variation.Occlusion area is also
So.This just increases difficulty to the area of computational shadowgraph and/or occlusion area.
Optionally, in embodiments of the present invention, determine at least one projected area in the surface of the first heliostat and
The area of collection region and determining union refion, which can be based on lattice Method, to be determined, specifically includes the following contents.
By the surface segmentation of at least one projected area and the first heliostat at grid, the size of medium square can be with
It is determined according to actual needs, grid is smaller, and precision is higher, therefore, can be calculated by controlling the size of grid to control
Precision.
Lattice point on the grid divided is marked, wherein if a lattice point is at least one projected area and
Preset value is marked in the case where any one of surface region of one heliostat.For example, preset value is 1, at least one is second fixed
Solar eyepiece includes 5 heliostats, then at least one projected area includes 5 projected areas, to mark for lattice point Q to be marked
Note lattice point is illustrated.It is lattice point Q label 1 when lattice point Q is in a projected area in 5 projected areas, in this way,
The each projected area whether lattice point Q is located in 5 projected areas is judged one by one, is then lattice when being located in a projected area
Point Q label 1；It when surface region of the lattice point Q also in the first heliostat, then is lattice point Q label 1.
The preset value that each lattice point is labeled is summed to obtain the label summation of the lattice point, wherein in the first heliostat
Surface region in, the region that occupies of twice of lattice point that label summation is greater than or equal to preset value is union refion.Continue with
The label summation for obtaining lattice point to summation for above-mentioned lattice point Q is illustrated.Assuming that lattice point Q is located in 5 projected areas
In two projections, while again in the surface region of the first heliostat, then lattice point Q is labeled 31 in total, then the mark of lattice point Q
Remember that summation is 3.The label summation of each lattice point is determined based on the mode similar with the label summation of lattice point Q is sought.Lattice point Q is in
In the surface region of first heliostat, while the label summation of lattice point Q is greater than 2, then illustrates lattice point Q while being also located at least one
In projected area, therefore, lattice point Q, which is in, can be generated the first heliostat in the region of shade.Similarly, in the first heliostat
In surface region, lattice point of the marked price summation more than or equal to 2 is to be in generate the first heliostat in the region of shade,
Therefore, in the surface region of the first heliostat, the region that lattice point of the label summation more than or equal to 2 occupies is at least one
Union refion of the projected area in the surface of the first heliostat.Preset value be not 1 be 2 in the case where, in the first heliostat
Surface region in, label summation had not only been at least one projected area more than or equal to 4 lattice point, but also was in for the first settled date
In the surface region of mirror.And then can obtain, in the surface region of the first heliostat, label summation is greater than or equal to the two of preset value
The region that lattice point again occupies is union refion of at least one projected area in the surface region of the first heliostat.
Determine that the label summation in the surface region of the first heliostat is greater than or equal to preset value based on Pick theorem
The area in the region that twice of lattice point occupies, wherein the area based on determined by Pick theorem is the area of union refion.After
Continuous to be illustrated with preset value for 1, in the surface region of the first heliostat, lattice point of the label summation equal to 2 is located at least one
The intersection of the surface region of a projected area and the first heliostat, that is, be located on the boundary of union refion, mark summation
Lattice point greater than 2 is located at the inside of union refion.Determine the lattice that the label summation in the surface region of the first heliostat is equal to 2
The number of the lattice point of number and label summation greater than 2 of point, that is, determine the sum for being located at the borderline lattice point of union refion
And the sum of the lattice point in union refion.Then the area of union refion is S=a+b ÷ 2-1, and wherein a indicates to be located at union
The sum of lattice point inside region, b indicate that the sum for being located at the borderline lattice point of union refion, s indicate the face of union refion
Product.
Optionally, in embodiments of the present invention, it in the case where determining union refion, is also based on integration method and determines
The area of union refion out.Specifically, union refion can be divided, is divided into and area is conveniently sought based on integration method
Figure, and then determine according to integration method the area of union refion.
Optionally, in embodiments of the present invention, it determines flat where at least one second heliostat is incident upon the first heliostat
At least one projected area on face includes: the coordinate that at least three vertex in the first heliostat are determined in the first coordinate system
The coordinate value on each vertex of every one second heliostat in value and at least one second heliostat, wherein the first coordinate system tool
Have following characteristics: integrating thermal tower location face position as coordinate origin, due east direction is x-axis and direct north as y-axis；Based on first
The coordinate value at least three vertex of heliostat determines the plane where the first heliostat；Based on light in the side of the first heliostat
The coordinate value on each vertex of every one second heliostat in parallactic angle and elevation angle and at least one second heliostat determines each
Straight line where vertex；Straight line and first where determining each vertex of every one second heliostat at least one second heliostat
The intersection point of plane where heliostat；And each vertex based on every one second heliostat at least one second heliostat
Corresponding intersection point determines the corresponding projected area of every one second heliostat, so that it is determined that at least one projected area.
Below at least one determining second heliostat by taking the heliostat in heliostat field is rectangular heliostat as an example
The first heliostat at least one projected area in the plane be illustrated.Specifically, (as among the above with heliostat b
Second heliostat) it generates on heliostat a (as among the above described the first heliostat) and is illustrated for shade.
Establish coordinate system S ', wherein in coordinate system S ', if the position where the collection thermal tower that heat collector is installed in is seat
Target origin (0,0,0), to integrate the direct north of thermal tower be positive direction of the y-axis, due east direction as positive direction of the x-axis.
Four vertex on three vertex and heliostat b of determining heliostat a coordinate value in coordinate system S '.According to 3 points
The principle for determining a plane determines the plane equation M where heliostat a according to the coordinate on the three of heliostat a vertex.
The intersection point that straight line and plane can be determined according to linear equation and plane equation determines the settled date based on the principle
Four vertex of mirror b be incident upon heliostat a intersection point planar.Determine heliostat b in heliostat a in the plane
When projected area, projected along sunray, determine heliostat b be incident upon heliostat a projection area planar
When domain, the vertex of heliostat b is located on sunray.Determine that four vertex of heliostat b are incident upon plane where heliostat a
Interior intersection point, that is, determine four vertex of heliostat b along plane where straight line where sunray and heliostat a
Intersection point.
Below by determine the vertex P of heliostat b along sunray be incident upon heliostat a carry out for point in the plane
Explanation.Determine other three vertex of heliostat b be incident upon heliostat a point in the plane may refer to determine that vertex P is thrown
Penetrate heliostat a point in the plane method.
Sunray is determined in the elevation angle α and azimuthal angle beta of heliostat a, so that it is determined that the posture of heliostat a out.Too
The elevation angle α of sun is sunray and earth's surface tangent line angle, the azimuthal angle beta of the sun be projection of the sunray on ground level with
Local meridianal angle, can approximatively regard the folder of the shade and Due South established straight line on the ground in the sun as
Angle.Azimuthal angle beta is zero with object direct north, is become larger clockwise, and value range is 0-360 °.Setting top
The coordinate of point P is (xb, yb, zb), then along the linear equation L where the P of sunray vertexPForSimultaneous plane equation M and linear equation LPVertex P can then be found out along the sun
Light be incident upon heliostat a point in the plane coordinate value.Based on same principle, determine heliostat b other three
A vertex is incident upon the intersection point on heliostat a.
Four vertex of heliostat b be incident upon heliostat a intersection point institute area defined in the plane be the settled date
Mirror b be incident upon heliostat a projected area in the plane.Similarly, other for determining can to generate heliostat a shade are fixed
Solar eyepiece be incident upon heliostat a projected area in the plane.
For determine to the first heliostat generate other heliostats for blocking be incident upon the first heliostat in the plane
Projected area can be issued according to the invertibity of light, reflection light by heat collector, determine other heliostats along transmitting light
Be incident upon the first heliostat intersection point in the plane, determine other heliostats be incident upon the first heliostat in the plane
Projected area, specifically, be referred to above-mentioned determining heliostat b be incident upon heliostat a projected area in the plane side
Method.
Optionally, in embodiments of the present invention, at least three vertex in the first heliostat are determined in the first coordinate system
Coordinate value and at least one second heliostat in the coordinate value on each vertex of every one second heliostat include the following contents.
Determine coordinate value of the midpoint of the first heliostat and each heliostat at least one second heliostat in the first coordinate system.
For example, at least second heliostat includes 5 the second heliostats, it is determined that every in the first heliostat and 5 the second heliostats
Coordinate value of the center of one heliostat in the first coordinate system.In corresponding second coordinate system of each heliostat, determine each
The coordinate value on the vertex of determining projected area is ready to use in heliostat, wherein the second coordinate system has the following characteristics that heliostat
Center is coordinate origin, to cross the normal direction at center of heliostat be y-axis, rotation corresponding with the elevation angle of first heliostat
Axis is x-axis.Continue by taking above-mentioned 5 the second heliostats as an example, determines that projected area needs to use at least three in the first heliostat
Each vertex of every one second heliostat in a vertex and 5 the second heliostats.It is based respectively on the first heliostat and 5 second
Heliostat establishes the second coordinate system, in corresponding second coordinate system, determines at least three vertex in the first heliostat
And the coordinate value on each vertex of every one second heliostat in 5 the second heliostats.Spin angle based on the first heliostat is faced upward
Angle, elevation angle and the transition matrix that the second coordinate system is converted to the first coordinate system is determined towards angle；And turn based on determined by
It changes matrix, be ready to use in determining projected area in each heliostat determined in corresponding second coordinate system of each heliostat
Coordinate value of the center of the coordinate value on vertex and each heliostat in the first coordinate system determines being ready to use in each heliostat
Determine coordinate value of the vertex of projected area in the first coordinate system.The elevation angle of first heliostat is relative altitude, for the
For one heliostat, target heat collector and the first heliostat are 0 degree in same level, are equal to 10 expression horizontal lines clockwise 10
Degree.First heliostat is direction relations towards angle, with target heat collector in the positive northern for 0 of the center of the first heliostat of ground
Degree, with target heat collector the first heliostat due west be 90 degree.As shown in Figure 2, target T and O is high in same level
Spending angle is 0 degree, and being equal to 10 indicates 10 degree clockwise of horizontal line；The due north of T point O is 0 degree, and T is 90 degree in the due west of O.
Optionally, in embodiments of the present invention, transition matrix are as follows:
Wherein, ρ is spin angle, and θ is the elevation angle, and λ is the elevation angle of the first heliostat,For towards angle.
Below with the seat on the vertex of determining in the first coordinate system heliostat a (being equal to the first heliostat described above)
It is illustrated for scale value.Where it is assumed that the heliostat in heliostat field is rectangle heliostat.Set heliostat a with
Ground parallel edges is long side, length L；Another side of heliostat a is broadside, width W.
Establish the first coordinate system S ', wherein in the first coordinate system S ', if where the collection thermal tower that heat collector is installed in
Position is origin (0,0,0), to integrate the direct north of thermal tower be positive direction of the y-axis, due east direction as positive direction of the x-axis.Settled date
The height on the centre distance ground of mirror a is H, and in the first coordinate system S ', the coordinate at the center of heliostat a is (x, y, H).
Establish the second coordinate system S, wherein in the second coordinate system S, using the center of heliostat a as coordinate origin (0,0,
0) normal direction for, crossing the center of heliostat a is y-axis, and the corresponding rotary shaft of the elevational movement of heliostat a is x-axis, wherein x-axis
It can be set according to the actual situation with the positive direction of y-axis, when being set at this to the positive direction of x-axis and y-axis, in conjunction with
Fig. 3 is illustrated.As shown in figure 3, with FF ' it is x-axis, FF ' is the corresponding rotary shaft of elevational movement of heliostat a, and is directed toward with F
F ' is positive direction of the x-axis, and with TT ' for y-axis, TT ' was the normal direction at the center of heliostat a, and being directed toward T ' with T is that y-axis is square
To.In addition, using the upward axis for the mirror surface origin for crossing heliostat a as z-axis, and using upwardly direction as positive direction.It needs to illustrate
, only it is illustrated herein with illustrating the positive direction of x-axis and y-axis in corresponding second coordinate system of heliostat a, x
The setting of the positive direction of axis and y-axis is not limited to that, and when setting the positive direction of other x-axis and y-axis, can also be adopted
Coordinate value is determined with the method for determining coordinate value described herein.
Four vertex for setting heliostat a are respectively A, B, C, D, as shown in figure 3, in the second coordinate system S, the seat of A point
It is designated as (L/2,0, W/2), the coordinate of B point is (- L/2,0, W/2), and the coordinate of C point is (- L/2,0 ,-W/2), the coordinate of D point
For (L/2,0 ,-W/2).
Coordinate system is divided into from the second coordinate system S to the coordinate of the first coordinate system S ' conversion one and converts three times.
The x-axis and y-axis for being converted to the second coordinate system S for the first time are obtained centered on y-axis to the negative direction of z-axis rotation angle ρ
To coordinate system s1, wherein ρ is the spin angle of heliostat a, the first transition matrix T are as follows:
Coordinate system s is converted to for the second time1Y-axis and z-axis centered on x-axis, to the negative direction of z-axis rotation angle be θ+λ,
Obtain coordinate system s2, wherein θ is the spin angle of heliostat a, and λ is the elevation angle of heliostat a, which is the collection for collecting thermal tower
The formation elevation angle of hot device and heliostat center.Elevation angle is relative altitude, for example, can be existed with collecting thermal tower and heliostat center
Same level is 0 degree, is equal to 10 degree of 10 degree clockwise of expression horizontal lines.Second transition matrix are as follows:
Third time is converted to coordinate system s2X-axis and y-axis centered on z-axis, to y-axis negative direction (from the negative direction of z-axis
Look counterclockwise to positive direction) rotationObtain coordinate system s3, whereinFor towards angle.It should be collection thermal tower and settled date towards angle
Mirror direction relations.For example, can with collect thermal tower heliostat due north for 0 degree, with collect thermal tower heliostat due west for 90 degree.
Third transition matrix are as follows:
Determining coordinate system s3After, by coordinate system s3The central point (x, y, H) for moving to heliostat a, obtains the first coordinate
It is S ',
According to above-mentioned conversion process it is found that transition matrix from the second coordinate system S to the first coordinate system S ' are as follows:Wherein, ρ is certainly
Swing angle, θ are the elevation angle, and λ is elevation angle,For towards angle.
The coordinate value for the point that heliostat a is in the second coordinate system S is transformed into the transformational relation in the first coordinate system S ' are as follows:
Wherein, s indicates the coordinate value o'clock in the first coordinate system S on heliostat a, and s ' is on heliostat a o'clock the
Coordinate value in two coordinate system S '.
In the second coordinate system S, the coordinate of A point is for (L/2,0, W/2), and the coordinate of B point is (- L/2,0, W/2), C point
Coordinate be (- L/2,0 ,-W/2), the coordinate of D point is (L/2,0 ,-W/2), is brought into above-mentioned transformational relation, and it is fixed to determine
Coordinate value of four vertex of solar eyepiece a in the first coordinate system S ' is respectively as follows:
A point coordinate:
B point coordinate:
C point coordinate:
D point coordinate:
For determining that the point (for example, the point for being ready to use in determining projected area) at least one second heliostat is sat first
Coordinate value in mark system, the method that may refer to the coordinate value o'clock in the first coordinate system on above-mentioned determining heliostat a, for
Every one second heliostat at least one second heliostat establishes corresponding second coordinate system S, and then based on similar to the above
Coordinate conversion, find the transformational relation that the corresponding second coordinate system S of every one second heliostat is transformed into the first coordinate system S ', into
And determine the coordinate value o'clock in the first coordinate system in every one second heliostat.It is calculated in heliostat field in this way, realizing
Heliostat real-time coordinates.
It should be noted that sunray emits on the first heliostat, according to the invertibity of light, reflection light can be with
Regard that the related content of all about determining eclipsing loss can refer to the correlation of determining shadow loss by heat collector sending as
Content understands and realizes.
Correspondingly, the another aspect of the embodiment of the present invention provides a kind of for determining the dress of the lost area of the first heliostat
It sets.Fig. 9 be another embodiment of the present invention provides for determine the first heliostat lost area device structural block diagram.Its
In, the mode for controlling first heliostat is that spin angle adds elevation angle control mode.As shown in figure 9, the device includes loss area
Domain determining module 1, loss heliostat determining module 2, projected area determining module 3, union refion determining module 4, lost area
Determining module 5.Loss area determination module 1 is for determining the corresponding loss region of the first heliostat, wherein loss region includes
Shadow region and/or occlusion area, the heliostat in shadow region generate shade on the first heliostat, are in blocked area
Heliostat in domain blocks the light reflected via the first heliostat；Heliostat determining module 2 is lost for determining in loss
At least one second heliostat in region；Projected area determining module 3 is for determining that at least one second heliostat is incident upon
First heliostat at least one projected area in the plane；Union refion determining module 4 is for determining at least one projection
Union refion of the region in the surface of the first heliostat；Lost area determining module 5 is used to determine the area of union refion,
The area of middle union refion is lost area.
It is under spin angle adds elevation angle control mode, to determine the damage for being in the first heliostat in the mode for controlling the first heliostat
Lose other heliostats in region the first heliostat at least one view field in the plane, determine at least one projection
Union refion of the region in the surface of the first heliostat, determines the area of union refion, and wherein the area of the union refion is
For the lost area of the first heliostat, lost area includes the area of shade and/or eclipsing loss, in this way, realizing determining place
The shade of heliostat and/or the area of eclipsing loss under spin angle plus elevation angle control mode.
Optionally, in embodiments of the present invention, union refion determining module determines that at least one projected area is fixed first
Union refion in the surface of solar eyepiece includes: the surface segmentation Cheng Fang by least one projected area and the first heliostat
Lattice；Lattice point on the grid divided is marked, wherein if a lattice point is at least one projected area and the first settled date
Preset value is marked in the case where any one of surface region of mirror；And the preset value summation for being labeled each lattice point is to obtain
To the label summation of the lattice point, wherein in the surface region of first heliostat, label summation is greater than or equal to preset value
The region that occupies of twice of lattice point be union refion；Lost area determining module determines that the area of the union refion includes:
Determine that the label summation in the surface region of first heliostat is greater than or equal to twice of preset value based on Pick theorem
The area in region that is occupied in the surface region of the first heliostat of lattice point, wherein the area based on determined by Pick theorem
The as area of union refion.
Optionally, in embodiments of the present invention, projected area determining module determines that at least one second heliostat is incident upon
First heliostat at least one projected area in the plane include: to be determined in the first heliostat extremely in the first coordinate system
The coordinate value on each vertex of the coordinate value and every one second heliostat at least one second heliostat on few three vertex,
In, the first coordinate system has the following characteristics that integrating thermal tower location face position is coordinate origin, due east direction as x-axis and the positive north
To for y-axis；The plane where the first heliostat is determined based on the coordinate value at least three vertex of the first heliostat；Based on light
Each vertex of every one second heliostat in the azimuth of the first heliostat and elevation angle and at least one second heliostat
Coordinate value determine straight line where each vertex；Determine each top of every one second heliostat at least one second heliostat
The intersection point of straight line where point and the plane where first heliostat；And based on each at least one second heliostat
The intersection point of each vertex correspondence of second heliostat determines the corresponding projected area of every one second heliostat, so that it is determined that at least one
A projected area.
Optionally, in embodiments of the present invention, at least three vertex in the first heliostat are determined in the first coordinate system
Coordinate value and at least one second heliostat in the coordinate value on each vertex of every one second heliostat comprise determining that first
Coordinate value of the midpoint of heliostat and each heliostat at least one second heliostat in the first coordinate system；Every certain
In corresponding second coordinate system of solar eyepiece, the coordinate value that the vertex of determining projected area is ready to use in each heliostat is determined, wherein
Second coordinate system have the following characteristics that heliostat center be coordinate origin, cross heliostat center normal direction be y-axis,
Rotary shaft corresponding with the elevation angle of the first heliostat is x-axis；Spin angle and the elevation angle, elevation angle based on the first heliostat and towards
Angle determines the transition matrix that the second coordinate system is converted to the first coordinate system；And based on identified transition matrix, every certain
In corresponding second coordinate system of solar eyepiece determine each heliostat in the vertex for being ready to use in determining projected area coordinate value and
Coordinate value of the center of each heliostat in first coordinate system, which determines in each heliostat, is ready to use in determining projection area
Coordinate value of the vertex in domain in the first coordinate system.
Optionally, in embodiments of the present invention, transition matrix are as follows:
Wherein, ρ is spin angle, and θ is the elevation angle, and λ is elevation angle,For towards angle.
It is provided in an embodiment of the present invention for determine the first heliostat lost area device concrete operating principle and
The concrete operating principle and benefit of benefit and the method for the lost area provided in an embodiment of the present invention for being used to determine the first heliostat
Locate similar, will not be described in great detail here.
In addition, the another aspect of the embodiment of the present invention also provides a kind of machine readable storage medium, the machine readable storage
Instruction is stored on medium, which is used for so that machine executes method described in above-described embodiment.
In conclusion being under spin angle adds elevation angle control mode, to determine and be in first in the mode for controlling the first heliostat
Other heliostats in the loss region of heliostat are at least one view field in the plane, the first heliostat institute, and determination is extremely
Few union refion of the view field in the surface of the first heliostat, determines the area of union refion, wherein the union area
The area in domain is the lost area of the first heliostat, and lost area includes the area of shade and/or eclipsing loss, in this way, real
Show and has determined that being in spin angle adds the area of the shade of heliostat and/or eclipsing loss under the control mode of the elevation angle.Furthermore, it is possible to really
The coordinate on the vertex on fixed first heliostat and each vertex of every one second heliostat at least one second heliostat, such as
This, realizes the coordinate value for determining the vertex of heliostat in real time.
The optional embodiment of the embodiment of the present invention is described in detail in conjunction with attached drawing above, still, the embodiment of the present invention is simultaneously
The detail being not limited in above embodiment can be to of the invention real in the range of the technology design of the embodiment of the present invention
The technical solution for applying example carries out a variety of simple variants, these simple variants belong to the protection scope of the embodiment of the present invention.
It is further to note that specific technical features described in the above specific embodiments, in not lance
In the case where shield, it can be combined in any appropriate way.In order to avoid unnecessary repetition, the embodiment of the present invention pair
No further explanation will be given for various combinations of possible ways.
It will be appreciated by those skilled in the art that implementing the method for the above embodiments is that can pass through
Program is completed to instruct relevant hardware, which is stored in a storage medium, including some instructions are used so that single
Piece machine, chip or processor (processor) execute all or part of the steps of each embodiment the method for the application.And it is preceding
The storage medium stated includes: USB flash disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory
The various media that can store program code such as (RAM, Random Access Memory), magnetic or disk.
In addition, any combination can also be carried out between a variety of different embodiments of the embodiment of the present invention, as long as it is not
The thought of the embodiment of the present invention is violated, equally should be considered as disclosure of that of the embodiment of the present invention.
Claims (11)
1. a kind of method for determining the lost area of the first heliostat, which is characterized in that control first heliostat
Mode adds elevation angle control mode for spin angle, this method comprises:
The corresponding loss region of first heliostat is determined, wherein the loss area domain includes shadow region and/or blocked area
Domain, the heliostat in the shadow region generates shade on first heliostat, in the occlusion area
Heliostat blocks the light reflected via first heliostat；
Determine at least one second heliostat in the loss region；
Determine at least one described second heliostat be incident upon first heliostat at least one projection area in the plane
Domain；
Determine union refion of at least one the described projected area in the surface of first heliostat；And
The area of the union refion is determined, wherein the area of the union refion is the lost area.
2. the method according to claim 1, wherein determining that at least one described projected area is fixed described first
Union refion in the surface of solar eyepiece and the area for determining the union refion include:
By the surface segmentation of at least one described projected area and first heliostat at grid；
Lattice point on the grid divided is marked, wherein if a lattice point is at least one described projected area and institute
Preset value is marked in the case where stating any one of surface region of the first heliostat；
The preset value that each lattice point is labeled is summed to obtain the label summation of the lattice point, wherein fixed described first
In the surface region of solar eyepiece, the region that twice of the lattice point that label summation is greater than or equal to the preset value occupies is the union
Region；And
It is described default to determine that the label summation in the surface region of first heliostat is greater than or equal to based on Pick theorem
The area in the region that twice of lattice point of value occupies in the surface region of first heliostat, wherein be based on Pick theorem
Identified area is the area of the union refion.
3. method according to claim 1 or 2, which is characterized in that at least one second heliostat described in the determination is thrown
Penetrate first heliostat at least one projected area in the plane include:
Determined in the first coordinate system at least three vertex in first heliostat coordinate value and it is described at least one the
The coordinate value on each vertex of every one second heliostat in two heliostats, wherein first coordinate system has the following characteristics that
Integrating thermal tower location face position as coordinate origin, due east direction is x-axis and direct north as y-axis；
The plane where first heliostat is determined based on the coordinate value at least three vertex of first heliostat；
It is each in the azimuth of first heliostat and elevation angle and at least one described second heliostat based on light
The coordinate value on each vertex of the second heliostat determines straight line where each vertex；
Straight line and described first where determining each vertex of every one second heliostat at least one described second heliostat
The intersection point of plane where heliostat；And
Based on each vertex correspondence of every one second heliostat at least one described second heliostat intersection point determine it is each
The corresponding projected area of second heliostat, so that it is determined that at least one described projected area.
4. according to the method described in claim 3, it is characterized in that, described determine first heliostat in the first coordinate system
In at least three vertex coordinate value and at least one described second heliostat in every one second heliostat each vertex
Coordinate value include:
Determine the midpoint of first heliostat and each heliostat at least one described second heliostat described first
Coordinate value in coordinate system；
In corresponding second coordinate system of each heliostat, the vertex that determining projected area is ready to use in each heliostat is determined
Coordinate value, wherein second coordinate system has the following characteristics that the center of heliostat is coordinate origin, the center for crossing heliostat
Normal direction be y-axis, rotary shaft corresponding with the elevation angle of first heliostat is x-axis；
Based on the spin angle of first heliostat, the elevation angle, elevation angle and determine second coordinate system to described towards angle
The transition matrix of one coordinate system conversion；And
Based on identified transition matrix, in corresponding second coordinate system of each heliostat determine each heliostat in
For determine the vertex of projected area coordinate value and each heliostat coordinate value of the center in first coordinate system it is true
Coordinate value of the vertex for being ready to use in determining projected area in first coordinate system in fixed each heliostat.
5. according to the method described in claim 4, it is characterized in that, the transition matrix are as follows:
Wherein, ρ is the spin angle, and θ is the elevation angle, and λ is the elevation angle of first heliostat,To be described towards angle.
6. a kind of for determining the device of the lost area of the first heliostat, which is characterized in that control first heliostat
Mode is that spin angle adds elevation angle control mode, which includes:
Area determination module is lost, for determining the corresponding loss region of first heliostat, wherein the loss area domain is wrapped
Shadow region and/or occlusion area are included, the heliostat in the shadow region generates shade on first heliostat,
Heliostat in the occlusion area blocks the light reflected via first heliostat；
Heliostat determining module is lost, for determining at least one second heliostat in the loss region；
Projected area determining module, it is flat where at least one described second heliostat is incident upon first heliostat for determining
At least one projected area on face；
Union refion determining module, for determine at least one described projected area in the surface of first heliostat and
Collect region；And
Lost area determining module, for determining the area of the union refion, wherein the area of the union refion is described
Lost area.
7. device according to claim 6, which is characterized in that
The union refion determining module determine at least one described projected area in the surface of first heliostat and
Collection region includes:
By the surface segmentation of at least one described projected area and first heliostat at grid；
Lattice point on the grid divided is marked, wherein if a lattice point is at least one described projected area and institute
Preset value is marked in the case where stating any one of surface region of the first heliostat；And
The preset value that each lattice point is labeled is summed to obtain the label summation of the lattice point, wherein fixed described first
In the surface region of solar eyepiece, the region that twice of the lattice point that label summation is greater than or equal to the preset value occupies is the union
Region；
The lost area determining module determines that the area of the union refion includes:
It is described default to determine that the label summation in the surface region of first heliostat is greater than or equal to based on Pick theorem
The area in the region that twice of lattice point of value occupies in the surface region of first heliostat, wherein be based on Pick theorem
Identified area is the area of the union refion.
8. device according to claim 6 or 7, which is characterized in that described in the projected area determining module is determining at least
One the second heliostat be incident upon first heliostat at least one projected area in the plane include:
Determined in the first coordinate system at least three vertex in first heliostat coordinate value and it is described at least one the
The coordinate value on each vertex of every one second heliostat in two heliostats, wherein first coordinate system has the following characteristics that
Integrating thermal tower location face position as coordinate origin, due east direction is x-axis and direct north as y-axis；
The plane where first heliostat is determined based on the coordinate value at least three vertex of first heliostat；
It is each in the azimuth of first heliostat and elevation angle and at least one described second heliostat based on light
The coordinate value on each vertex of the second heliostat determines straight line where each vertex；
Straight line and described first where determining each vertex of every one second heliostat at least one described second heliostat
The intersection point of plane where heliostat；And
Based on each vertex correspondence of every one second heliostat at least one described second heliostat intersection point determine it is each
The corresponding projected area of second heliostat, so that it is determined that at least one described projected area.
9. device according to claim 8, which is characterized in that described to determine first heliostat in the first coordinate system
In at least three vertex coordinate value and at least one described second heliostat in every one second heliostat each vertex
Coordinate value include:
Determine the midpoint of first heliostat and each heliostat at least one described second heliostat described first
Coordinate value in coordinate system；
In corresponding second coordinate system of each heliostat, the vertex that determining projected area is ready to use in each heliostat is determined
Coordinate value, wherein second coordinate system has the following characteristics that the center of heliostat is coordinate origin, the center for crossing heliostat
Normal direction be y-axis, rotary shaft corresponding with the elevation angle of first heliostat is x-axis；
Based on the spin angle of first heliostat, the elevation angle, elevation angle and determine second coordinate system to described towards angle
The transition matrix of one coordinate system conversion；And
Based on identified transition matrix, in corresponding second coordinate system of each heliostat determine each heliostat in
For determine the vertex of projected area coordinate value and each heliostat coordinate value of the center in first coordinate system it is true
Coordinate value of the vertex for being ready to use in determining projected area in first coordinate system in fixed each heliostat.
10. device according to claim 9, which is characterized in that the transition matrix are as follows:
Wherein, ρ is the spin angle, and θ is the elevation angle, and λ is the elevation angle of first heliostat,To be described towards angle.
11. a kind of machine readable storage medium, it is stored with instruction on the machine readable storage medium, which is used for so that machine
Perform claim requires method described in any one of 1-5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910375494.4A CN110209207B (en) | 2019-05-07 | 2019-05-07 | Method and apparatus for determining lost area of heliostat and machine-readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910375494.4A CN110209207B (en) | 2019-05-07 | 2019-05-07 | Method and apparatus for determining lost area of heliostat and machine-readable storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110209207A true CN110209207A (en) | 2019-09-06 |
| CN110209207B CN110209207B (en) | 2022-05-24 |
Family
ID=67785589
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910375494.4A Active CN110209207B (en) | 2019-05-07 | 2019-05-07 | Method and apparatus for determining lost area of heliostat and machine-readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110209207B (en) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101859014A (en) * | 2010-06-09 | 2010-10-13 | 中国科学院电工研究所 | Azimuth-elevation tracking method of heliostat with reflecting mirror surface deviated from rotating center |
| CN102519152A (en) * | 2011-11-30 | 2012-06-27 | 浙江大学 | Calculating method for mirror field optical efficiency on basis of graphics processing unit (GPU) tower type solar energy thermoelectric system |
| CN102930160A (en) * | 2012-11-02 | 2013-02-13 | 浙江大学 | Calculating method of mirror field optical efficiency of tower type solar thermoelectric system |
| CN102997453A (en) * | 2012-11-22 | 2013-03-27 | 宁夏光合能源科技有限公司 | Method for controlling heliostats in tower type solar-powered heat collecting device |
| JP2013174363A (en) * | 2010-08-05 | 2013-09-05 | Cosmo Oil Co Ltd | Solar concentrating system and heliostat layout method |
| CN103500277A (en) * | 2013-09-29 | 2014-01-08 | 浙江大学 | Mirror field optimization design method of cornfield and tower type solar thermoelectric system |
| CN106056670A (en) * | 2016-05-30 | 2016-10-26 | 浙江大学 | Occlusion-eliminating radiation energy density simulation method in tower-type solar thermal power generation system |
| CN107168381A (en) * | 2017-05-08 | 2017-09-15 | 东南大学 | A kind of impact analysis method of figure of the earth to heliostat solar tracking error |
| CN107704711A (en) * | 2017-10-30 | 2018-02-16 | 中国华能集团清洁能源技术研究院有限公司 | A kind of tower type solar mirror field shade and the innovatory algorithm for blocking efficiency |
| JP2019047623A (en) * | 2017-09-01 | 2019-03-22 | Ｓｆｉ株式会社 | Photovoltaic power generation method |
| 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 |
-
2019
- 2019-05-07 CN CN201910375494.4A patent/CN110209207B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101859014A (en) * | 2010-06-09 | 2010-10-13 | 中国科学院电工研究所 | Azimuth-elevation tracking method of heliostat with reflecting mirror surface deviated from rotating center |
| JP2013174363A (en) * | 2010-08-05 | 2013-09-05 | Cosmo Oil Co Ltd | Solar concentrating system and heliostat layout method |
| CN102519152A (en) * | 2011-11-30 | 2012-06-27 | 浙江大学 | Calculating method for mirror field optical efficiency on basis of graphics processing unit (GPU) tower type solar energy thermoelectric system |
| CN102930160A (en) * | 2012-11-02 | 2013-02-13 | 浙江大学 | Calculating method of mirror field optical efficiency of tower type solar thermoelectric system |
| CN102997453A (en) * | 2012-11-22 | 2013-03-27 | 宁夏光合能源科技有限公司 | Method for controlling heliostats in tower type solar-powered heat collecting device |
| CN103500277A (en) * | 2013-09-29 | 2014-01-08 | 浙江大学 | Mirror field optimization design method of cornfield and tower type solar thermoelectric system |
| CN106056670A (en) * | 2016-05-30 | 2016-10-26 | 浙江大学 | Occlusion-eliminating radiation energy density simulation method in tower-type solar thermal power generation system |
| CN107168381A (en) * | 2017-05-08 | 2017-09-15 | 东南大学 | A kind of impact analysis method of figure of the earth to heliostat solar tracking error |
| JP2019047623A (en) * | 2017-09-01 | 2019-03-22 | Ｓｆｉ株式会社 | Photovoltaic power generation method |
| CN107704711A (en) * | 2017-10-30 | 2018-02-16 | 中国华能集团清洁能源技术研究院有限公司 | A kind of tower type solar mirror field shade and the innovatory algorithm for blocking efficiency |
| 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 (6)
| Title |
|---|
| K.K. CHONG ⇑, M.H. TAN: "Range of motion study for two different sun-tracking methods in the application of heliostat field", 《SOLAR ENERGY》 * |
| 刘培杰数学工作室: "《Pick定理》", 31 May 2017, 哈尔滨工业大学出版社 * |
| 张立功 等: "塔式太阳能热发电站中定日镜的电源方案研究", 《太阳能》 * |
| 田枫 等: "塔式太阳能热发电相关技术的最新进展", 《南京师范大学学报（工程技术版）》 * |
| 耿其东: "塔式太阳能热发电中的定日镜跟踪系统设计", 《热力发电》 * |
| 谢飞: "塔式太阳能热电系统定日镜场光学仿真与应用研究", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110209207B (en) | 2022-05-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Noone et al. | Heliostat field optimization: A new computationally efficient model and biomimetic layout | |
| Lin et al. | Experimental and theoretical analysis on a linear Fresnel reflector solar collector prototype with V-shaped cavity receiver | |
| CN102519152A (en) | Calculating method for mirror field optical efficiency on basis of graphics processing unit (GPU) tower type solar energy thermoelectric system | |
| CN104034058B (en) | The formation method of tower type solar heat and power system Jing Chang based on GPU | |
| Yousuf et al. | Solar energy potential estimation by calculating sun illumination hours and sky view factor on building rooftops using digital elevation model | |
| CN105320156A (en) | An automatic tracking method for a tower-type solar thermal power generation heliostat | |
| CN106352566A (en) | Tower type solar thermal power generation sun-tracking analysis system | |
| CN109670248A (en) | The analogy method of the radiant energy dfensity distribution of hot spot in tower type solar energy thermal power generation | |
| CN107388599B (en) | A kind of shade of linear Fresnel formula light condenser field and sheltering analysis optimize distribution method | |
| Ulmer et al. | Slope measurements of parabolic dish concentrators using color-coded targets | |
| Lee et al. | Numerical study on optical performances of the first central-receiver solar thermal power plant in Korea | |
| Huang et al. | Development of a new flux density function for a focusing heliostat | |
| Hiller et al. | TRNSHD—a program for shading and insolation calculations | |
| CN110442959A (en) | The prediction of solar energy disc type condensation system condenser performance and keeping method | |
| Bonanos et al. | Heliostat surface shape characterization for accurate flux prediction | |
| He et al. | Fast flux density distribution simulation of central receiver system on GPU | |
| Riaz | A theory of concentrators of solar energy on a central receiver for electric power generation | |
| Zhang et al. | Model construction and performance analysis for asymmetric compound parabolic concentrator with circular absorber | |
| Hou et al. | Optical performance investigation on flat receiver for parabolic trough solar collector based on the MCRT method | |
| CN109684593A (en) | Mountainous region photovoltaic project array pitch calculation method and system | |
| CN106910245B (en) | A kind of image modeling method of parabolic trough concentrators | |
| Skouri et al. | Optical qualification of a solar parabolic concentrator using photogrammetry technique | |
| Hu et al. | Performance analysis and optimization of an integrated azimuth tracking solar tower | |
| CN110209207A (en) | Determine the method and apparatus and machine readable storage medium of the lost area of heliostat | |
| Huang et al. | Solar flux density calculation for a heliostat with an elliptical Gaussian distribution source |
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 |