CN107704711A - A kind of tower type solar mirror field shade and the innovatory algorithm for blocking efficiency - Google Patents
A kind of tower type solar mirror field shade and the innovatory algorithm for blocking efficiency Download PDFInfo
- Publication number
- CN107704711A CN107704711A CN201711036023.8A CN201711036023A CN107704711A CN 107704711 A CN107704711 A CN 107704711A CN 201711036023 A CN201711036023 A CN 201711036023A CN 107704711 A CN107704711 A CN 107704711A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- heliostat
- shade
- point
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
Abstract
A kind of tower type solar mirror field shade and the innovatory algorithm for blocking efficiency, by establishing shade and blocking the pre- judgment models of innovatory algorithm, prediction can produce shade or the region blocked, and combine Monte Carlo ray tracing model and calculate mirror field shade in shade or occlusion area and block efficiency, pre- judgment models and the correctness of innovatory algorithm are demonstrated by actual power station data, the Jing Chang for further preferred arrangement high efficiency, low cost provides good theoretical foundation.
Description
Technical field
The present invention relates to photo-thermal power generation technology, more particularly to a kind of tower type solar mirror field shade and the improvement for blocking efficiency
Algorithm.
Background technology
It is becoming tight environment situation day now, the extreme climate such as haze, coal rain emerges in an endless stream, every profession and trade energy-saving and emission-reduction standard day
Benefit improves, as clean energy resource, contribution significant of the solar energy thermal-power-generating station to environment.Tower type solar energy thermal power generation technology because
Collecting efficiency is high, thermal technology's high conversion efficiency, and system overall efficiency is high, and cost reduction space is big, more energy-conserving and environment-protective, suitable big rule
Mould apply the advantages that and as photo-thermal industry future large-scale application Main way.But tower power station mirror field cost is very high
It is expensive, the 30%~50% of whole power station total investment is accounted for, therefore design high efficiency, the Jing Chang of low cost for promoting tower
Formula solar energy power plant is very necessary.And the factor for influenceing mirror field optics efficiency has a lot, wherein on shade and blocking efficiency
Calculating requires that variable is most, and difficulty in computation is also maximum.It is some to produce for compared with Mi Jingchang in traditional flat projection model
The heliostat of raw shade can be ignored, and cause shade with blocking efficiency calculation or judging that mistiming difference can be very big in advance.
The content of the invention
In order to overcome the above-mentioned deficiencies of the prior art, it is an object of the invention to provide a kind of tower type solar mirror field shade
Innovatory algorithm with blocking efficiency, the shade influenceed in tower type solar mirror field optics efficiency can be accurately acquired and block effect
Rate.
In order to achieve the above object, the technical solution adopted by the present invention is:
A kind of tower type solar mirror field shade and the innovatory algorithm for blocking efficiency, comprise the following steps:
(1) establish shade and block the pre- judgment models of innovatory algorithm;
It is assumed that problem heliostat is heliostat 2, the heliostat of meeting shade or occlusion issue heliostat is heliostat 1, is utilized
Ray tracing can obtain shade or block projected position (x of the mirror heart of heliostat 1 in the surface coordinate system of heliostat 21,h2, y1,h2,
z1,h2), for circular Jing Chang, heliostat 1 is not substantially parallel with heliostat 2, when heliostat 1 is projected to heliostat 2
The figure of gained has no longer been a rectangle in itself, but the parallelogram A'B'C'D', a length of H of a rotationh, it is a width of
Hw;
It is a rectangle A " B " C " D " by the projective amplification of heliostat 1, and four summits projected are the four of the rectangle
On bar side, rectangle A " B " C " D " a length of φy×Hh, a width of φx×Hw, wherein φx, φyLong and wide amplification system is represented respectively
Number, when heliostat 1 may produce shade to heliostat 2 or block, condition need to be met:
Also can be reduced to:Wherein, ψx、ψyFor the correction factor of x-axis and y-axis, andWith φ relation
For:
In the process of optimization of mirror field, ψ and φ value change with mirror field shape, heliostat change in location, therefore
To ensure computational accuracy, it is necessary to estimate that each heliostat φ is worth accordingly.φ calculation formula are:
The relation concluded between solar eyepiece 2 and heliostat 1 can be prejudged with reference to above-mentioned formula, when judge two heliostats occur
Shade or when blocking, need to record projection coordinate (x1, h2, y1, h2, z1, h2) and ψ of the mirror heart of heliostat 1 on the minute surface of heliostat 2
Or φ data, used in case calculating link in next step;
(2) precision of prediction of pre- judgment models is obtained;
To obtain the predictablity rate of forecast model, the shade precision of prediction η pr, sd of heliostat mirror field need to be calculated and blocked
Precision of prediction η pr, bl,
Wherein, nn represents neighbouring heliostat sum, and sd represents shade heliostat;Bl is represented and is blocked heliostat, subscript ac and
Pr represents virtual condition and predicted state respectively, nsdac (i), nblac (i) and nsdpr (i), nblpr (i) represent respectively for
Heliostat quantity under i-th of heliostat virtual condition with that can produce shade under predicted state or block.It is fixed that statistics mirror place has
Solar eyepiece, the precision of prediction for producing the shade of whole heliostat mirror field or blocking;
(3) shade is obtained using Monte Carlo ray tracing model and blocks efficiency;
Shade or the heliostat field blocked can be produced in pre- judgement, uses Monte Carlo ray tracing model computational shadowgraph
With blocking efficiency.Monte Carlo ray tracing model receives tower, heat dump surface around center, target heliostat, closes on the settled date
5 sub-models of mirror and sunray build corresponding coordinate system, improve the mutual transformation model of each coordinate system and build up, and the model is based on system
Meter learns principle, and uniform discrete is carried out to each heliostat surface, passes through the ray tracing to sample point, statistics heliostat surface
The sample point quantity of reachable heat dump receiving plane, calculate the optical delivery efficiency of each heliostat, single heliostat and surrounding
Ray tracing logic flow is between certain adjacent heliostat:
(3-1) dispersive target heliostat surface, obtain all sample points position and to center in the heliostat coordinate system
Tower coordinate system is changed;
(3-2) extracts certain sample point, establishes ray along solar incident ray direction in the point, judges ray and adjacent settled date
Mirror parabolic model whether there is intersection point;
Whether (3-3) is changed if intersection point be present, by the intersection point to adjacent heliostat coordinate system, judge intersection point in heliostat table
In the range of face, if intersection point is not present, illustrate that this closes on shadow-free, into step (3-5);
(3-4) illustrates that the point by shade, directly calculates next typical application, is transferred to step (3- if in the range of minute surface
2), if not in the range of, into step (3-5);
(3-5) is established ray along reflection light in sample point, judged by mirror-reflection definite principle reflected ray vector
Whether with close heliostat parabolic model intersection point be present;
Whether (3-6) is changed if intersection point be present, by the intersection point to adjacent heliostat coordinate system, judge intersection point in heliostat table
In the range of face, if intersection point is not present, illustrate that this closes on shadow-free, into step (3-8);
(3-7) if in the range of minute surface, the point has been blocked, and directly calculates next typical application, is transferred to step (3-2),
If not in the range of, carry out in next step;
(3-8) now the point by shade and blocks checking, for the sample point indirect ray, judges itself and heat dump
There is intersection point in surface model, if intersection point is not present, the point is not intercepted, calculates next sample point, is transferred to step (3-2);
Whether (3-9) is changed if intersection point be present, by the point to heat dump surface coordinate system, judge the point in receiving plane scope
Within, if not in the range of, the point is not intercepted, calculates next sample point, is transferred to step (3-2), if in the range of, the point is
It is intercepted, receives number and add 1, calculate the cosine losses and air dissipation loss, be transferred to next sample point, be transferred to step (3-
2);
(3-10) counts all receiving point numbers and the heliostat shade is obtained compared with total sample point and blocks efficiency.
The described target heliostat light gathering efficiency in heliostat field calculates, adjacent heliostat be it is multiple, only need to be
Shade judges that iteration checking each closes on heliostat respectively for region (2-4) and occlusion area (5-8), you can completes shade with blocking
Efficiency calculation.
The beneficial effects of the invention are as follows:
The present invention is for existing tower type solar mirror field shade and the deficiency for blocking efficiency calculation, it is proposed that one kind is tower too
Positive energy mirror field shade and the innovatory algorithm for blocking efficiency, by establishing shade and blocking the pre- judgment models of innovatory algorithm, obtain accurate
True mirror field shade and shielded area, and combine Monte Carlo ray tracing model and obtain mirror field shade and block efficiency, to enter
One-step optimization arrangement high efficiency, the Jing Chang of low cost provide good theoretical foundation.
Brief description of the drawings
Fig. 1 is a kind of tower type solar mirror field shade of the present invention and the innovatory algorithm flow chart for blocking efficiency.
Fig. 2 is heliostat shade of the present invention and blocks projection theory schematic diagram.
Fig. 3 is projection modification schematic diagram of the present invention.
Fig. 4 is shade of the present invention and the precision of prediction figure for blocking pre- judgment models.
Fig. 5 is ray tracing modular concept figure in Monte Carlo of the present invention.
Embodiment
Make narration in detail to the present invention below in conjunction with the accompanying drawings.
Fig. 1 is a kind of tower type solar mirror field shade of the present invention and the innovatory algorithm flow chart for blocking efficiency.
As shown in figure 1, a kind of innovatory algorithm of tower type solar mirror field shade with blocking efficiency, comprises the following steps:
Step S1, establish shade and block the pre- judgment models of innovatory algorithm;;
Step S2, obtain the precision of prediction of pre- judgment models;
Step S3, shade is obtained using Monte Carlo ray tracing model and blocks efficiency;
In a particular embodiment of the present invention, step S1 is implemented as:
It is assumed that problem heliostat is heliostat 2, the heliostat of meeting shade or occlusion issue heliostat is heliostat 1.Utilize
Ray tracing can obtain shade or block projected position (x of the mirror heart of heliostat 1 in the surface coordinate system of heliostat 21,h2,y1,h2,
z1,h2).For circular Jing Chang, heliostat 1 is not substantially parallel with heliostat 2, when heliostat 1 is projected to heliostat 2
The figure of gained has no longer been a rectangle in itself, but the parallelogram A'B'C'D', a length of H of a rotationh, it is a width of
Hw。
It is a rectangle A " B " C " D " by the projective amplification of heliostat 1, and four summits projected are the four of the rectangle
On bar side, rectangle A " B " C " D " a length of φy×Hh, a width of φx×Hw, wherein φx, φyLong and wide amplification system is represented respectively
Number.When heliostat 1 may produce shade to heliostat 2 or block, condition need to be met:
Also can be reduced to:Wherein, ψx、ψyFor the correction factor of x-axis and y-axis, and ψ and φ relation
For:ψxory=φxory+1/2。
In the process of optimization of mirror field, ψ and φ value change with mirror field shape, heliostat change in location.Therefore
To ensure computational accuracy, it is necessary to estimate that each heliostat φ is worth accordingly.φ calculation formula are:
The relation concluded between solar eyepiece 2 and heliostat 1 can be prejudged with reference to above-mentioned formula, judge be between two heliostats
No generation shade blocks.When judging that two heliostats occur shade or blocked, the mirror heart of heliostat 1 need to be recorded in heliostat 2
Projection coordinate (x on minute surface1,h2, y1,h2, z1,h2) with ψ or φ data.
In a particular embodiment of the present invention, step S2 is implemented as:
To obtain the predictablity rate of forecast model, the shade precision of prediction η of heliostat mirror field need to be calculatedpr,sdIt is pre- with blocking
Survey precision ηpr,bl,
Wherein, nnNeighbouring heliostat sum is represented, sd represents shade heliostat;Bl is represented and is blocked heliostat, subscript ac and
Pr represents virtual condition and predicted state respectively.nsdac(i)、nblacAnd nsd (i)pr(i)、nblpr(i) represent respectively for i-th
Heliostat quantity under individual heliostat virtual condition with that can produce shade under predicted state or block.There is the settled date in statistics mirror place
Mirror, the precision of prediction for producing the shade of whole heliostat mirror field or blocking.
In a particular embodiment of the present invention, step S3 is implemented as:
As shown in Figure 5:Shade or the heliostat field blocked can be produced in pre- judgement, uses Monte Carlo ray tracing mould
Type computational shadowgraph is with blocking efficiency.Monte Carlo ray tracing model receives tower, heat dump surface, target settled date around center
Mirror, close on 5 sub-models of heliostat and sunray and build corresponding coordinate system, improve the mutual transformation model of each coordinate system and build up.
The model is based on Principle of Statistics, and uniform discrete is carried out to each heliostat surface, passes through the ray tracing to sample point, system
The sample point quantity that heliostat surface reaches heat dump receiving plane is counted, calculates the optical delivery efficiency of each heliostat.It is single
Ray tracing logic flow is between heliostat and certain adjacent heliostat around:
(3-1) dispersive target heliostat surface, obtain all sample points position and to center in the heliostat coordinate system
Tower coordinate system is changed;
(3-2) extracts certain sample point, establishes ray along solar incident ray direction in the point, judges ray and adjacent settled date
Mirror parabolic model whether there is intersection point;
Whether (3-3) is changed if intersection point be present, by the intersection point to adjacent heliostat coordinate system, judge intersection point in heliostat table
In the range of face, if intersection point is not present, illustrate that this closes on shadow-free, into step (3-5);
(3-4) illustrates that the point by shade, directly calculates next typical application, is transferred to step (3- if in the range of minute surface
2), if not in the range of, into step (3-5);
(3-5) is established ray along reflection light in sample point, judged by mirror-reflection definite principle reflected ray vector
Whether with close heliostat parabolic model intersection point be present;
Whether (3-6) is changed if intersection point be present, by the intersection point to adjacent heliostat coordinate system, judge intersection point in heliostat table
In the range of face, if intersection point is not present, illustrate that this closes on shadow-free, into step (3-8);
(3-7) if in the range of minute surface, the point has been blocked, and directly calculates next typical application, is transferred to step (3-2),
If not in the range of, carry out in next step;
(3-8) now the point by shade and blocks checking, for the sample point indirect ray, judges itself and heat dump
There is intersection point in surface model, if intersection point is not present, the point is not intercepted, calculates next sample point, is transferred to step (3-2);
Whether (3-9) is changed if intersection point be present, by the point to heat dump surface coordinate system, judge the point in receiving plane scope
Within, if not in the range of, the point is not intercepted, calculates next sample point, is transferred to step (3-2), if in the range of, the point is
It is intercepted, receives number and add 1, calculate the cosine losses and air dissipation loss.Next sample point is transferred to, is transferred to step (3-
2)。
(3-10) counts all receiving point numbers and the heliostat shade is obtained compared with total sample point and blocks efficiency.
Calculated for the target heliostat light gathering efficiency in heliostat field, adjacent heliostat is multiple, need to only be sentenced in shade
Iteration checking each closes on heliostat respectively for disconnected region (step (3-2)-(3-4)) and occlusion area (step (3-5)-(3-8)),
Shade can be completed and block efficiency calculation.
It is application of the method for the invention in specific example below, concrete application scene is solar power tower.
Because mirror field heliostat uses closely spaced array, there will be cast shadow during optically focused between each other and reflection hides
Phenomenon is kept off, heliostat 1 is target heliostat, and heliostat 2 is shade or blocks heliostat.In traditional flat projection model,
When projection lens's heart coordinate | xh2,1| H wide less than heliostatwAnd | yh2,1| H high less than heliostathWhen, heliostat 2 may be to the settled date
Mirror 1 produces shade or blocked, as shown in Figure 2.
Conventional model compared with Mi Jingchang for being calculated or being judged in advance mistiming poor meeting very big, some settled dates that can produce shade
Mirror can be ignored.This is primarily due to when heliostat 2 is projected to heliostat 1, and heliostat 1 is not with heliostat 2
It is substantially parallel, projection gained figure has no longer been a rectangle in itself, and seems the parallelogram of a rotation, is such as schemed
Shown in 3.
As shown in figure 3, the projection of heliostat 1 is amplified into A " B " C " D " from A'B'C'D', amplification coefficient φ, modification are introduced
Traditional projection judgment models, and estimate that each heliostat amplification coefficient φ is worth accordingly, establish heliostat 2 and heliostat 1
Anticipation relation, so as to judge shade whether occurs between two heliostats or block.
The circle most Mi Jingchang of Collado F J propositions is chosen as reference mirror field, arrangement is classical radial orders
Ladder arrangement, mirror field parameters selection Gemasolar power station datas, as shown in table 1, can obtain shade or the screening of whole heliostat mirror field
The precision of prediction of gear, as shown in Figure 4.The shade precision of prediction of pre- judgment models proposed by the present invention approaches as seen from Figure 4
100%, and occlusion prediction precision is close to 85%.The relatively low heliostat entity that is primarily due to of occlusion prediction precision is parabolic lens, instead
It is not directional light to penetrate light yet, significantly optical parallax be present.
In the tower mirror field of reality, according to pre- judgment models, projecting figure of the heliostat 1 on heliostat 2 can be approximate
For the rectangle of an amplification, projection coordinate (x of the mirror heart of heliostat 1 on the minute surface of heliostat 2 is being obtainedh2,1, yh2,1, zh2,1)
Shade or the area blocked can be produced between can determine that heliostat afterwards with ψ or φ data.Chased after using Monte Carlo ray
Track model, as shown in figure 5, for the target heliostat computational shadowgraph in heliostat field and blocking efficiency, adjacent heliostat is more
It is individual, only it need to judge that iteration checking each closes on heliostat respectively for region and occlusion area in shade, you can complete shade with blocking
Efficiency calculation can finally obtain mirror field shade and block efficiency.
It is as shown in table 2 to quote PS10 mirror field design data, the shade and block effect that inventive algorithm model is calculated
Rate ηsbFor 0.9216, from Windelsol algorithms result of calculation as reference value, obtained shade and efficiency eta is blockedsbFor
0.9255, it is sufficiently close to inventive algorithm model result, nuance may be caused by position and calculation error, so far this hair
Bright algorithm model is had verified that to be correct.
Table 1 is the first extensive tower molten salt power station Gemasolar Power Plant Design tables of data in the whole world;
Table 2 is Spain's tower type solar business power station PS10 Power Plant Design tables of data;
Claims (2)
1. a kind of tower type solar mirror field shade and the innovatory algorithm for blocking efficiency, it is characterised in that comprise the following steps:
(1) establish shade and block the pre- judgment models of innovatory algorithm;
It is assumed that problem heliostat is heliostat 2, the heliostat of meeting shade or occlusion issue heliostat is heliostat 1, utilizes ray
Tracking can obtain shade or block projected position (x of the mirror heart of heliostat 1 in the surface coordinate system of heliostat 21,h2, y1,h2,
z1,h2), for circular Jing Chang, heliostat 1 is not substantially parallel with heliostat 2, when heliostat 1 is projected to heliostat 2
The figure of gained has no longer been a rectangle in itself, but the parallelogram A'B'C'D', a length of H of a rotationh, it is a width of
Hw;
It is a rectangle A " B " C " D " by the projective amplification of heliostat 1, and four summits projected are in the four edges of the rectangle
On, rectangle A " B " C " D " a length of φy×Hh, a width of φx×Hw, wherein φx, φyLong and wide amplification coefficient is represented respectively, when
When heliostat 1 may produce shade or block to heliostat 2, condition need to be met:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>&phi;</mi>
<mi>x</mi>
</msub>
<mo>&times;</mo>
<msub>
<mi>H</mi>
<mi>w</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
<mo>+</mo>
<mo>|</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>h</mi>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>|</mo>
<mo><</mo>
<mo>=</mo>
<msub>
<mi>H</mi>
<mi>w</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>&phi;</mi>
<mi>y</mi>
</msub>
<mo>&times;</mo>
<msub>
<mi>H</mi>
<mi>h</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
<mo>+</mo>
<mo>|</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>h</mi>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>|</mo>
<mo><</mo>
<mo>=</mo>
<msub>
<mi>H</mi>
<mi>h</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
Also can be reduced to:Wherein, ψx、ψyFor the correction factor of x-axis and y-axis, andRelation with φ is:
In the process of optimization of mirror field, ψ and φ value change with mirror field shape, heliostat change in location, therefore to protect
Demonstrate,prove computational accuracy, it is necessary to estimate that each heliostat φ is worth accordingly, φ calculation formula are:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&phi;</mi>
<mi>x</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>max</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<msup>
<mi>A</mi>
<mo>&prime;</mo>
</msup>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>h</mi>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<msup>
<mi>B</mi>
<mo>&prime;</mo>
</msup>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>h</mi>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<msup>
<mi>C</mi>
<mo>&prime;</mo>
</msup>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>h</mi>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<msup>
<mi>D</mi>
<mo>&prime;</mo>
</msup>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>h</mi>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>H</mi>
<mi>w</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&phi;</mi>
<mi>y</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>max</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>y</mi>
<mrow>
<msup>
<mi>A</mi>
<mo>&prime;</mo>
</msup>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>h</mi>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mrow>
<msup>
<mi>B</mi>
<mo>&prime;</mo>
</msup>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>h</mi>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mrow>
<msup>
<mi>C</mi>
<mo>&prime;</mo>
</msup>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>h</mi>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mrow>
<msup>
<mi>D</mi>
<mo>&prime;</mo>
</msup>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>h</mi>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>H</mi>
<mi>h</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
The relation concluded between solar eyepiece 2 and heliostat 1 can be prejudged with reference to above-mentioned formula, when judge two heliostats occur shade
Or when blocking, projection coordinate (x1, h2, y1, h2, z1, h2) of the mirror heart of heliostat 1 on the minute surface of heliostat 2 and ψ or φ need to be recorded
Data, in case in next step calculate link use;
(2) precision of prediction of pre- judgment models is obtained;
To obtain the predictablity rate of forecast model, the shade precision of prediction η pr, sd and occlusion prediction of heliostat mirror field need to be calculated
Precision η pr, bl,
<mrow>
<msub>
<mi>&eta;</mi>
<mrow>
<mi>p</mi>
<mi>r</mi>
<mo>,</mo>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>n</mi>
<mi>n</mi>
</msub>
</munderover>
<msub>
<mi>nsd</mi>
<mrow>
<mi>a</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>n</mi>
<mi>n</mi>
</msub>
</munderover>
<msub>
<mi>nsd</mi>
<mrow>
<mi>p</mi>
<mi>r</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>,</mo>
<msub>
<mi>&eta;</mi>
<mrow>
<mi>p</mi>
<mi>r</mi>
<mo>,</mo>
<mi>b</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>n</mi>
<mi>n</mi>
</msub>
</munderover>
<msub>
<mi>nbl</mi>
<mrow>
<mi>a</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>n</mi>
<mi>n</mi>
</msub>
</munderover>
<msub>
<mi>nbl</mi>
<mrow>
<mi>p</mi>
<mi>r</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
Wherein, nn represents neighbouring heliostat sum, and sd represents shade heliostat;Bl, which is represented, blocks heliostat, ac and pr points of subscript
Virtual condition and predicted state are not represented, and nsdac (i), nblac (i) and nsdpr (i), nblpr (i) represent for i-th respectively
There is the settled date with the heliostat quantity that shade can be produced under predicted state or is blocked, statistics mirror place under individual heliostat virtual condition
Mirror, the precision of prediction for producing the shade of whole heliostat mirror field or blocking;
(3) shade is obtained using Monte Carlo ray tracing model and blocks efficiency;
Shade or the heliostat field blocked can be produced in pre- judgement, using Monte Carlo ray tracing model computational shadowgraph with hiding
Keep off efficiency, Monte Carlo ray tracing model around center receive tower, heat dump surface, target heliostat, close on heliostat and
5 sub-models of sunray build corresponding coordinate system, improve the mutual transformation model of each coordinate system and build up, the model is based on statistics
Principle, uniform discrete is carried out to each heliostat surface, by the ray tracing to sample point, statistics heliostat surface can be arrived
Up to the sample point quantity of heat dump receiving plane, the optical delivery efficiency of each heliostat, single heliostat and surrounding phase are calculated
Ray tracing logic flow is between adjacent heliostat:
(3-1) dispersive target heliostat surface, obtain all sample points in the heliostat coordinate system position and to central tower sit
The conversion of mark system;
(3-2) extracts certain sample point, and ray is established along solar incident ray direction in the point, judges that ray is thrown with adjacent heliostat
Object plane model whether there is intersection point;
Whether (3-3) is changed if intersection point be present, by the intersection point to adjacent heliostat coordinate system, judge intersection point in heliostat surface model
In enclosing, if intersection point is not present, illustrate that this closes on shadow-free, into step (3-5);
(3-4) illustrates that the point by shade, directly calculates next typical application, is transferred to step (3-2) if in the range of minute surface, if
Not in the range of, into step (3-5);
(3-5) is established ray along reflection light in sample point, judged whether by mirror-reflection definite principle reflected ray vector
Intersection point be present with close heliostat parabolic model;
Whether (3-6) is changed if intersection point be present, by the intersection point to adjacent heliostat coordinate system, judge intersection point in heliostat surface model
In enclosing, if intersection point is not present, illustrate that this closes on shadow-free, into step (3-8);
(3-7) if in the range of minute surface, the point has been blocked, and directly calculates next typical application, step (3-2) is transferred to, if not
In the range of, carry out in next step;
(3-8) now the point by shade and blocks checking, for the sample point indirect ray, judges itself and heat dump surface
There is intersection point in model, if intersection point is not present, the point is not intercepted, calculates next sample point, is transferred to step (3-2);
(3-9) if intersection point be present, by the point to heat dump surface coordinate system change, judge the point whether receiving plane scope with
Interior, if not in the range of, the point is not intercepted, calculates next sample point, is transferred to step (3-2), if in the range of, the point by
Interception, receive number and add 1, calculate the cosine losses and air dissipation loss, be transferred to next sample point, be transferred to step (3-2);
(3-10) counts all receiving point numbers and the heliostat shade is obtained compared with total sample point and blocks efficiency.
2. a kind of tower type solar mirror field shade according to claim 1 and the innovatory algorithm for blocking efficiency, its feature exist
In the target heliostat light gathering efficiency in heliostat field calculates, and adjacent heliostat is multiple, need to only be sentenced in shade
Iteration checking each closes on heliostat respectively for disconnected region (2-4) and occlusion area (5-8), you can completes shade and blocks efficiency meter
Calculate.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711036023.8A CN107704711A (en) | 2017-10-30 | 2017-10-30 | A kind of tower type solar mirror field shade and the innovatory algorithm for blocking efficiency |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711036023.8A CN107704711A (en) | 2017-10-30 | 2017-10-30 | A kind of tower type solar mirror field shade and the innovatory algorithm for blocking efficiency |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107704711A true CN107704711A (en) | 2018-02-16 |
Family
ID=61176977
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711036023.8A Pending CN107704711A (en) | 2017-10-30 | 2017-10-30 | A kind of tower type solar mirror field shade and the innovatory algorithm for blocking efficiency |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107704711A (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108800618A (en) * | 2018-04-28 | 2018-11-13 | 华西能源工业股份有限公司 | Tower type solar energy thermal power generation mirror field method for arranging and its shadow occlusion efficiency calculation method |
CN109992882A (en) * | 2019-04-01 | 2019-07-09 | 浙江大学 | A kind of solar radiant energy density emulation mode based on Monte Carlo ray trace |
CN110209207A (en) * | 2019-05-07 | 2019-09-06 | 中国神华能源股份有限公司 | Determine the method and apparatus and machine readable storage medium of the lost area of heliostat |
CN110414057A (en) * | 2019-06-28 | 2019-11-05 | 浙江大学 | The radiant energy dfensity analogy method of focus type heliostat in tower type solar thermo-power station |
CN110414059A (en) * | 2019-06-28 | 2019-11-05 | 浙江大学 | The radiant energy dfensity analogy method of planar heliostats in tower type solar thermo-power station |
CN110658857A (en) * | 2019-10-08 | 2020-01-07 | 浙江正泰新能源开发有限公司 | Method and device for verifying tracking accuracy of photovoltaic tracking system |
CN117527991A (en) * | 2024-01-08 | 2024-02-06 | 北京智汇云舟科技有限公司 | Video occlusion region detection method and system based on ray tracing |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102930160A (en) * | 2012-11-02 | 2013-02-13 | 浙江大学 | Calculating method of mirror field optical efficiency of tower type solar thermoelectric system |
CN103530697A (en) * | 2013-09-29 | 2014-01-22 | 浙江大学 | Mirror field optimal design method of radiant tower type solar thermoelectric system |
CN105973505A (en) * | 2016-05-10 | 2016-09-28 | 西安交通大学 | Method for determining heat-flow density in opening of solar cavity type heat absorber |
US20160370032A1 (en) * | 2014-07-22 | 2016-12-22 | Esolar Inc. | Variable Density Heliostat Field Layout |
-
2017
- 2017-10-30 CN CN201711036023.8A patent/CN107704711A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102930160A (en) * | 2012-11-02 | 2013-02-13 | 浙江大学 | Calculating method of mirror field optical efficiency of tower type solar thermoelectric system |
CN103530697A (en) * | 2013-09-29 | 2014-01-22 | 浙江大学 | Mirror field optimal design method of radiant tower type solar thermoelectric system |
US20160370032A1 (en) * | 2014-07-22 | 2016-12-22 | Esolar Inc. | Variable Density Heliostat Field Layout |
CN105973505A (en) * | 2016-05-10 | 2016-09-28 | 西安交通大学 | Method for determining heat-flow density in opening of solar cavity type heat absorber |
Non-Patent Citations (3)
Title |
---|
KUN WANG ET AL.: "A novel integrated simulation approach couples MCRT and Gebhart methods to simulate solar radiation transfer in a solar power tower system with a cavity receiver", 《RENEWABLE ENERGY》 * |
张茂龙等: "塔式太阳能镜场阴影与遮挡效率的改进算法", 《太阳能学报》 * |
胡甜等: "基于光线踪迹法的塔式太阳能镜场布置与优化研究", 《工程热物理学报》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108800618A (en) * | 2018-04-28 | 2018-11-13 | 华西能源工业股份有限公司 | Tower type solar energy thermal power generation mirror field method for arranging and its shadow occlusion efficiency calculation method |
CN108800618B (en) * | 2018-04-28 | 2019-11-26 | 华西能源工业股份有限公司 | Tower type solar energy thermal power generation mirror field method for arranging and its shadow occlusion efficiency calculation method |
CN109992882A (en) * | 2019-04-01 | 2019-07-09 | 浙江大学 | A kind of solar radiant energy density emulation mode based on Monte Carlo ray trace |
CN109992882B (en) * | 2019-04-01 | 2020-08-11 | 浙江大学 | Monte Carlo ray tracing-based solar radiation energy density simulation method |
CN110209207A (en) * | 2019-05-07 | 2019-09-06 | 中国神华能源股份有限公司 | Determine the method and apparatus and machine readable storage medium of the lost area of heliostat |
CN110414057A (en) * | 2019-06-28 | 2019-11-05 | 浙江大学 | The radiant energy dfensity analogy method of focus type heliostat in tower type solar thermo-power station |
CN110414059A (en) * | 2019-06-28 | 2019-11-05 | 浙江大学 | The radiant energy dfensity analogy method of planar heliostats in tower type solar thermo-power station |
CN110414059B (en) * | 2019-06-28 | 2020-12-01 | 浙江大学 | Radiation energy density simulation method of planar heliostat in tower type solar thermal power station |
CN110414057B (en) * | 2019-06-28 | 2020-12-11 | 浙江大学 | Radiation energy density simulation method of focusing heliostat in tower type solar thermal power station |
CN110658857A (en) * | 2019-10-08 | 2020-01-07 | 浙江正泰新能源开发有限公司 | Method and device for verifying tracking accuracy of photovoltaic tracking system |
CN117527991A (en) * | 2024-01-08 | 2024-02-06 | 北京智汇云舟科技有限公司 | Video occlusion region detection method and system based on ray tracing |
CN117527991B (en) * | 2024-01-08 | 2024-03-15 | 北京智汇云舟科技有限公司 | Video occlusion region detection method and system based on ray tracing |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107704711A (en) | A kind of tower type solar mirror field shade and the innovatory algorithm for blocking efficiency | |
CN102519152B (en) | Calculating method for mirror field optical efficiency on basis of graphics processing unit (GPU) tower type solar energy thermoelectric system | |
CN102930160B (en) | Determine the method for tower type solar heat and power system mirror field optics efficiency | |
CN111881576B (en) | Heliostat field optimization scheduling control method for solar tower type photo-thermal power station | |
CN109670248A (en) | The analogy method of the radiant energy dfensity distribution of hot spot in tower type solar energy thermal power generation | |
CN106369838B (en) | A kind of slot light collection solar thermal collection system design method | |
Xu et al. | An optimization design method and optical performance analysis on multi-sectioned compound parabolic concentrator with cylindrical absorber | |
CN108090275A (en) | A kind of turbine air cooling blade configuration method towards parameterized model | |
CN110276168A (en) | The tower non-central point focusing modeling method of photo-thermal power station heliostat field | |
CN104034058A (en) | Imaging method of tower-type solar thermoelectric system mirror field based on GPU | |
CN107388599A (en) | A kind of shade of linear Fresnel formula light condenser field and sheltering analysis optimization distribution method | |
CN110705077B (en) | Method for calculating energy flow density distribution of focusing light spots of tower-type solar heat absorber | |
Rungasamy et al. | Comparative study of the optical and economic performance of etendue-conserving compact linear Fresnel reflector concepts | |
CN106910245B (en) | A kind of image modeling method of parabolic trough concentrators | |
Hu et al. | A comprehensive mathematical approach and optimization principle for solar flux distribution and optical efficiency in a solar tower | |
WO2023115740A1 (en) | Method and system for calculating solar radiation energy received by building | |
CN111459194A (en) | Solar thermal power generation aiming point determination method based on heliostat measured light spot | |
CN114329903B (en) | Method for calculating shadow efficiency of heliostat of tower type solar photo-thermal power station | |
CN109460594A (en) | A kind of dish-style Order Triangular Elements split paraboloid film condenser condenser performance prediction technique | |
CN109615693A (en) | A kind of optimized calculation method of flat single-shaft support spacing and following range | |
CN106918313A (en) | A kind of dish-style Stirling Salar light-gatherings minute surface quality determining method | |
CN110414058B (en) | Fitting method of convolution kernel function in tower type solar thermal power station | |
CN111767646B (en) | Tower type solar thermal power station receiver optimization design method | |
CN110414057B (en) | Radiation energy density simulation method of focusing heliostat in tower type solar thermal power station | |
Serrano-Arrabal et al. | Dual-tower CSP plants: optical assessment and optimization with a novel cone-tracing model |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180216 |