METHOD FOR DETERMINING TARGET POINTS OF HELIOSTATS DURING PREHEATING OF TOWER-TYPE SOLAR PHOTO-THERMAL POWER STATION

Abstract

A method for determining target points of heliostats during preheating of a tower-type solar photo-thermal power station comprises: establishing a coordinate system of a heliostat field of the station; obtaining coordinates of each heliostat according to a layout of the heliostat field; obtaining vertex coordinates of each heat absorbing panel on a heat absorber according to a layout of the heat absorbers; carrying out grid division for each panel to obtain vertex coordinates of each grid; obtaining X and Y coordinates of the target point of each heliostat on the panel; taking a Z coordinate of the target point of each heliostat on the panel as an independent variable and a sum of squares of differences between an actual number and an expected number of target points in each grid as an objective function to establish a non-linear optimization model, and solving the model to obtain the Z coordinate.

Claims

1. A method for determining target points of heliostats during preheating of a tower-type solar photo-thermal power station, comprising the following steps: step 1, establishing a coordinate system of a heliostat field of the tower-type solar photo-thermal power station; step 2, obtaining coordinates of each heliostat according to a layout of the heliostat field of the tower-type solar photo-thermal power station; step 3, obtaining vertex coordinates of each heat absorbing panel on a heat absorber according to a layout of the heat absorbers; step 4, carrying out grid division for each heat absorbing panel of the heat absorber to obtain vertex coordinates of each grid; step 5, obtaining an X coordinate and a Y coordinate of the target point of each heliostat on the heat absorbing panel according to the principle of minimum distance; and step 6, according to the divided grids, taking a Z coordinate of the target point of each heliostat on the heat absorbing panel as an independent variable and a sum of squares of differences between an actual number of target points and an expected number of target points in each grid as an objective function to establish a non-linear optimization model, and solving the model to obtain the Z coordinate of the target point.

2. The method for determining target points of heliostats during preheating of a tower-type solar photo-thermal power station according to claim 1, wherein in the step 1, an East-North-Up coordinate system is established with a center point of a ground where a heat absorption tower of the tower-type solar photo-thermal power station is located as an origin of coordinate.

3. The method for determining target points of heliostats during preheating of a tower-type solar photo-thermal power station according to claim 1, wherein in the step 2, the coordinates of the heliostat refer to coordinates of a center point of a mirror of the heliostat, and coordinates of a ith heliostat are obtained according to a geometric layout of the heliostat field and expressed as (hx.sub.i, hy.sub.i, hz.sub.i), i=1 . . . . W, where W is the number of heliostats in the heliostat field.

4. The method for determining target points of heliostats during preheating of a tower-type solar photo-thermal power station according to claim 1, wherein in the step 3, the heat absorber is a regular n-gonal prism structure which comprises N heat absorbing panels, and each heat absorbing panel is rectangular; and coordinates of a jth vertex of a kth heat absorbing panel are obtained according to a height T of the heat absorption tower, a height H of the heat absorbing panel and a length L of a bottom edge of the heat absorbing panel and expressed as (Tx.sub.kj, Ty.sub.kj, Tz.sub.kj), where k=1 . . . . N, and j=1, 2, 3 and 4.

5. The method for determining target points of heliostats during preheating of a tower-type solar photo-thermal power station according to claim 4, wherein in the step 4, for the kth heat absorbing panel, a maximum value Xmax.sub.k and a minimum value Xmin.sub.k of an X axis and a maximum value Zmax.sub.k and a minimum value Zmin.sub.k of a Z axis are respectively obtained:
Xmax.sub.k=max(Tx.sub.k1, . . . ,Tx.sub.k4)
Xmin.sub.k=min(Tx.sub.k1, . . . ,Tx.sub.k4)
Zmax.sub.k=max(Tz.sub.k1, . . . ,Tz.sub.k4)
Zmin.sub.k=min(Tz.sub.k1, . . . ,Tz.sub.k4) the kth heat absorbing panel is divided into MF grids, and vertex coordinates of all grids in the X-axis direction are respectively $X{\min }_k,X{\min }_k+\frac{X{\max }_k-X{\min }_k}{M-1},.Math.,X{\max }_k,$ and the vertex coordinates of all grids in the Z-axis direction are respectively $Z{\min }_k,Z{\min }_k+\frac{Z{\max }_k-Z{\min }_k}{F-1},.Math.,Z{\max }_k.$

6. The method for determining target points of heliostats during preheating of a tower-type solar photo-thermal power station according to claim 4, wherein the step 5 specifically comprises the following steps: according to the principle of minimum distance, a horizontal distance from each heliostat to the target point is considered as a minimum distance from each heliostat to the heat absorption tower to obtain a circumcircle corresponding to a projection of the heat absorption tower on a horizontal plane: $x^2+y^2=\frac{L^2}{4.Math.{\tan }^2\left(\frac{}{N}\right)}$ where x and y are the X coordinate and the Y coordinate of the target point of the heliostat respectively; for the ith heliostat (hx.sub.i, hy.sub.i, hz.sub.i), a following equation is simultaneous: $\{\begin{array}{c}{x}^2+y^2=\frac{L^2}{4.Math.{\tan }^2\left(\frac{}{N}\right)}\\ y=\frac{{\mathrm{hy}}_i}{{\mathrm{hx}}_i}x\end{array}$ and the X coordinate and the Y coordinate of the target point of the ith heliostat are obtained by solving, $x=\frac{{\mathrm{Lhx}}_i}{2.Math.\tan \left(\frac{}{N}\right)\sqrt{{\mathrm{hx}}_i^2+{\mathrm{hy}}_i^2}}y=\frac{{\mathrm{Lhy}}_i}{2.Math.\tan \left(\frac{}{N}\right)\sqrt{{\mathrm{hx}}_i^2+{\mathrm{hy}}_i^2}}.$

7. The method for determining target points of heliostats during preheating of a tower-type solar photo-thermal power station according to claim 5, wherein in the step 6, the established non-linear optimization model is as follows: $\min {.Math.\stackrel{.Math.}{C}\left(z_1,.Math.,z_W\right)-\stackrel{.Math.}{N}.Math.}^2$ $s.tZ\min z_i\mathrm{Zman}$ where (z.sub.1, . . . , z.sub.W)=(c.sub.1, . . . , c.sub.NMF) refers to an actual number of target points in each grid on the heat absorbing panel under a Z coordinate of a current corresponding target point of the heliostat; where z.sub.1 refers to a Z coordinate of a corresponding target point of a first heliostat, z.sub.W refers to a Z coordinate of a corresponding target point of a wth heliostat, c.sub.1 refers to an actual number of target points in a first grid, C.sub.NMF refers to an actual number of target points in a NMFth grid, and a grid to which the Z coordinate of the target point belongs is determined according to the vertex coordinates of the grids obtained previously, and finally, the number of the target points corresponding to each grid is summarized; NMF refers to a total number of grids; $\stackrel{.Math.}{N}=\left(n_1,.Math.,n_{NMF}\right)=\left(\frac{W}{NMF},.Math.,\frac{W}{NMF}\right)$ where n.sub.1 refers to an expected number of target points in the first grid, n.sub.NMF refers to an expected number of target points in a NMFth grid, and is W a total number of heliostats; z.sub.i refers to a Z coordinate of the corresponding target point of the ith heliostat, Zmin is a minimum value of a Z coordinate of the heat absorbing panel, and Zmax is a maximum value of the Z coordinate of the heat absorbing panel; and a MATLAB tool is used to solve an optimal solution of the non-linear optimization model, namely the Z coordinate of the target point of each heliostat.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0030] In order to illustrate the technical solutions of the embodiments of the present invention or the prior art more clearly, the accompanying drawings required for the description of the embodiments or the prior art will be briefly introduced below.

[0031] FIG. 1 is a schematic view of a heliostat field of a tower-type solar photo-thermal power station;

[0032] FIG. 2 is a flow chart of a method for determining target points of heliostats during preheating of a tower-type solar photo-thermal power station;

[0033] FIG. 3 is a structural schematic view of a heat absorber; and

[0034] FIG. 4 is a schematic diagram of a circumcircle corresponding to a projection of a heat absorption tower on a horizontal plane.

[0035] Reference numerals in the figures: 1. heliostat; 2. heat absorber; 3. heat absorption tower; 4. sun; 5. heat absorbing panel.

DESCRIPTION OF THE EMBODIMENTS

[0036] The technical solutions of the embodiments of the present invention will be described clearly and completely below with reference to the accompanying drawings in the embodiments of the present invention.

[0037] The present invention provides a method for determining target points of heliostats during preheating of a tower-type solar photo-thermal power station. As shown in FIG. 2, the method comprises the following steps:

[0038] Step 1, establishing a coordinate system of a heliostat field of the tower-type solar photo-thermal power station; at which an East-North-Up (ENU) coordinate system is established with a center point of a ground where a heat absorption tower 3 of the tower-type solar photo-thermal power station is located as an origin of coordinate.

[0039] Step 2, obtaining coordinates of each heliostat 1 according to a layout of the heliostat field of the tower-type solar photo-thermal power station; at which the coordinates of the heliostat 1 refer to coordinates of a center point of a mirror of the heliostat 1, and coordinates of a ith heliostat 1 are obtained according to a geometric layout of the heliostat field and expressed as (hx.sub.i, hy.sub.i, hz.sub.i), i=1 . . . . W, where W is the number of heliostats 1 in the heliostat field.

[0040] Step 3, obtaining vertex coordinates of each heat absorbing panel 5 on a heat absorber 2 according to a layout of the heat absorbers 2; as shown in FIG. 3, the heat absorber 2 is a regular n-gonal prism structure which comprises N heat absorbing panels 5, and each heat absorbing panel 5 is rectangular; and coordinates of a jth vertex of a kth heat absorbing panel 5 are obtained according to a height T of the heat absorption tower 3, a height H of the heat absorbing panel 5 and a length L of a bottom edge of the heat absorbing panel 5 and expressed as (Tx.sub.kj, Ty.sub.kj, Tz.sub.kj), where k=1 . . . . N, and j=1, 2, 3 and 4.

[0041] Step 4, carrying out grid division for each heat absorbing panel 5 of the heat absorber 2 to obtain vertex coordinates of each grid; [0042] for the kth heat absorbing panel, a maximum value Xmax.sub.k and a minimum value Xmin.sub.k of an X axis and a maximum value Zmax.sub.k and a minimum value Zmin.sub.k of a Z axis are respectively obtained:

Xmax.sub.k=max(Tx.sub.k1, . . . ,Tx.sub.k4)

Xmin.sub.k=min(Tx.sub.k1, . . . ,Tx.sub.k4)

Zmax.sub.k=max(Tz.sub.k1, . . . ,Tz.sub.k4)

Zmin.sub.k=min(Tz.sub.k1, . . . ,Tz.sub.k4) [0043] the kth heat absorbing panel 5 is divided into MF grids, and vertex coordinates of all grids in the X-axis direction are respectively

[00008]$X{\min }_k,X{\min }_k+\frac{X{\max }_k-X{\min }_k}{M-1},.Math.,X{\max }_k,$ and the vertex coordinates of all grids in the Z-axis direction are respectively

[00009]$Z{\min }_k,Z{\min }_k+\frac{Z{\max }_k-Z{\min }_k}{F-1},.Math.,Z{\max }_k.$

[0044] Step 5, obtaining an X coordinate and a Y coordinate of the target point of each heliostat 1 on the heat absorbing panel 5 according to the principle of minimum distance; as shown in FIG. 4, according to the principle of minimum distance, a horizontal distance from each heliostat 1 to the target point is considered as a minimum distance from each heliostat 1 to the heat absorption tower 3 to obtain a circumcircle corresponding to a projection of the heat absorption tower 3 on a horizontal plane:

[00010]$x^2+y^2=\frac{L^2}{4.Math.{\tan }^2\left(\frac{}{N}\right)}$ [0045] where x and y are the X coordinate and the Y coordinate of the target point of the heliostat 1 respectively; [0046] for the ith heliostat 1 (hx.sub.i, hy.sub.i, hz.sub.i), a following equation is simultaneous:

[00011]$\{\begin{array}{c}{x}^2+y^2=\frac{L^2}{4.Math.{\tan }^2\left(\frac{}{N}\right)}\\ y=\frac{{\mathrm{hy}}_i}{{\mathrm{hx}}_i}x\end{array}$ [0047] the X coordinate and the Y coordinate of the target point of the ith heliostat 1 are obtained by solving,

[00012]$x=\frac{{\mathrm{Lhx}}_i}{2.Math.\tan \left(\frac{}{N}\right)\sqrt{{\mathrm{hx}}_i^2+{\mathrm{hy}}_i^2}}y=\frac{{\mathrm{Lhy}}_i}{2.Math.\tan \left(\frac{}{N}\right)\sqrt{{\mathrm{hx}}_i^2+{\mathrm{hy}}_i^2}}.Math.\mathrm{Step}$

[0048] The established non-linear optimization model is as follows:

[00013]$\min {.Math.\stackrel{.Math.}{C}\left(z_1,.Math.,z_W\right)-\stackrel{.Math.}{N}.Math.}^2$$s.tZ\min z_i\mathrm{Zman}$ [0049] where (z.sub.1, . . . , z.sub.W)=(c.sub.1, . . . , c.sub.NMF) refers to the actual number of target points in each grid on the heat absorbing panel 5 under a Z coordinate of a current corresponding target point of the heliostat 1; where z.sub.1 refers to a Z coordinate of a corresponding target point of a first heliostat 1, z.sub.W refers to a Z coordinate of a corresponding target point of a wth heliostat 1, c.sub.1 refers to the actual number of target points in a first grid, C.sub.NMF refers to the actual number of target points in a NMFth grid, and a grid to which the Z coordinate of the target point belongs is determined according to the vertex coordinates of the grids obtained previously, and finally, the number of the target points corresponding to each grid is summarized; NMF refers to the total number of grids;

[00014]$\stackrel{.Math.}{N}=\left(n_1,.Math.,n_{NMF}\right)=\left(\frac{W}{NMF},.Math.,\frac{W}{NMF}\right)$ [0050] where n.sub.1 refers to an expected number of target points in the first grid, n.sub.NMF refers to an expected number of target points in a NMFth grid, and is W a total number of heliostats 1; [0051] z.sub.i refers to a Z coordinate of the corresponding target point of the ith heliostat 1, Zmin is a minimum value of a Z coordinate of the heat absorbing panel 5, and Zmax is a maximum value of the Z coordinate of the heat absorbing panel 5; and [0052] a MATLAB tool is used to solve an optimal solution of the non-linear optimization model, namely the Z coordinate of the target point of each heliostat 1.

[0053] The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the present invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.