CN106777634B - 基于Arcgis的光伏阵列阴影计算方法 - Google Patents
基于Arcgis的光伏阵列阴影计算方法 Download PDFInfo
- Publication number
- CN106777634B CN106777634B CN201611116780.1A CN201611116780A CN106777634B CN 106777634 B CN106777634 B CN 106777634B CN 201611116780 A CN201611116780 A CN 201611116780A CN 106777634 B CN106777634 B CN 106777634B
- Authority
- CN
- China
- Prior art keywords
- points
- photovoltaic
- row
- arrays
- shadow
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Photovoltaic Devices (AREA)
Abstract
本发明涉及光伏电站领域,公开了一种基于Arcgis的光伏阵列阴影计算方法,使得研究人员能够确定任意时间段内前排光伏组件的阴影对后排光伏组件的遮挡情况。本发明通过对太阳的方向向量、组件支架、地面进行建模,从而推导出光伏组件的边缘点到地面投影点的变换矩阵,通过前后组件阴影的重合区域通过变换矩阵求逆准确得到后排组件的被遮挡的区域。本发明适用于对地势非常复杂的山地区域的光伏电站。
Description
技术领域
本发明涉及光伏电站领域,特别涉及基于Arcgis的光伏阵列阴影计算方法。
背景技术
目前按照《光伏发电站设计规范GB50797-2012》中对光伏方阵布置的要求:光伏方阵各排、列的布置间距应保证每天9:00~15:00(当地真太阳时)时间段内前、后、左、右互补遮挡。
由于研究光伏组件的阴影的主要目的为研究前排光伏组件对后排组件的遮挡问题。由此,目前国内的主要研究方向为在保证前排组件对后排组件不遮挡的情况下光伏组件布置的最小间距。此最小间距以当地冬至日真太阳时9:00~15:00时无遮挡为目标,不能具体的分析在其他时间段内前排光伏组件的阴影对后排光伏组件的遮挡情况。
同时,目前的公式仅能对分析在理想的平地、南坡或者北坡下的光伏组件阴影,由于山地光伏地形、地势非常复杂,微地形和微地貌很多。使用上述常规的分析方法无法实现。
发明内容
本发明要解决的技术问题是:提供一种基于Arcgis的光伏阵列阴影计算方法,使得研究人员能够确定任意时间段内前排光伏组件的阴影对后排光伏组件的遮挡情况。
为解决上述问题,本发明采用的技术方案是:基于Arcgis的光伏阵列阴影计算方法,包括步骤:
建立三维坐标系o-xyz;
建立太阳法向量S(Sx,Sy,Sz)在三维坐标系中的数学模型;
通过Arcgis软件获取光伏支架的任意三个角点在三维坐标系中的坐标,根据获取的光伏支架的任意三个角点坐标求得光伏支架的平面方程T(Tx,Ty,Tz),建立光伏支架的平面在三维坐标系中的数学模型;
通过Arcgis软件获取坡面的任意三个点在三维坐标系中的坐标,根据获取的光伏地面的任意三个点的坐标求得坡面的平面方程G(Gx,Gy,Gz),建立坡面在三维坐标系中的数学模型;
根据太阳法向量的数学模型、坡面的的数学模型、光伏支架的平面的数学模型计算光伏组件边缘点T(Tx,Ty,Tz)在坡面上的投影点V(Vx,Vy,Vz)的变化矩阵M;
利用变化矩阵M计算前后两排光伏阵列的边缘点在坡面上的投影点;
根据前后两排光伏阵列的边缘点在坡面上的投影点确定前后两排光伏阵列的阴影交集;
若所述阴影交集不是空集,则获取所述阴影交集的顶点;
根据阴影交集的顶点以及变换矩阵M求逆得到前排阵列在后排阵列上的遮挡区域的顶点,连接所有的前排阵列在后排阵列的遮挡区域的顶点,从而得到前排阵列在后排阵列的遮挡区域。
进一步的,所述变换矩阵M为:
其中,N=Sx·Gx+Sy·Gy+Sz·Gz。
进一步的,获取阴影交集T12′(x,y)的四个顶点的步骤包括:
A.若存在P1、P2、P3、P4四个点互不相等的点,则遍历结束,P1、P2、P3、P4分别为T12′的顶点;
C.若只存在两个不同点,将这两个点做为顶点;同时还需连接这两个点,通过连接的两点求得两点的连线L(x,y)=0,遍历点集T12′(x,y),求出距离连线L(x,y)距离正值和负值最大的两点,并作为T12′的顶点。
本发明的有益效果是:目前对于光伏阵列间距的计算主要是建立在理想的坡面模型的基础上,通过当地纬度、经度、太阳高度角、方位角、前后排支架高差以及支架布置的地形坡度等要素求解光伏布置的最小间距,此方法适用于大型的地面光伏电站的设计。地面光伏电站地形简单,前后排间距大多与固定不变。但山地光伏地形复杂,微地貌较多,通过传统的方法无法计算出前排支架阴影对后排支架的的阴影影响范围。本发明通过对方向向量、组件支架、地面进行建模,推导出光伏组件的边缘点到地面投影点的变换矩阵,可以计算出任意时刻、任意位置的光伏支架的阴影,具有通用性。并且推导出前排支架对地表的投影矩阵,通过Arcgis软件编程可实现光伏阵列的阴影模拟。
2、传统的计算方法,无法计算出前排阵列对后排阵列的影响区域大小。本发明通过算法找到前后两排的阴影区域重合部分,然后利用阴影变换矩阵求逆得到前排阵列在后排阵列的遮挡区域。这有利于准确计算当前后排组件存在遮挡时,发电量的损失的大小。
附图说明
图1为阴影形成的三要素图;
图2为太阳的坐标系示意图;
图3为阴影的公共区域的顶点的判断示意图。
具体实施方式
以下结合附图对本发明做进一步说明。
1、太阳法向量的建模
如图1所示,阴影的行成至少需要三个要素,光线λ、物体P和投影面π,。一个物体只有在光线照射下,遮挡住了物体背面的另一个物体(投影面π),才会在背面的物体上留***影P。根据这三个要素进行分析和表达,可得到光伏支架组件到阴影的变换矩阵。
由于太阳的位置和地球的位置是不断变化的,因此太阳照射下物体的阴影也在不停的改变。在地平坐标系下,太阳相对地面的位置由太阳的高度角α和方位角β决定。建立太阳的方向向量则可以准确的表达出太阳所在位置。
如图2所示,太阳的方向向量S=[Sx,Sy,Sz]对太阳的球坐标转换成直角坐标后,可得到Sx,Sy,Sz与太阳的高度角α和方位角β的关系如式(7)所示:
高度角α和太阳的方位角β的计算公式为:
太阳赤纬角δ
时角ω
ω=(TT-12)×15°
真太阳时TT(h):
其中,T为北京时间(h),λ当地经度。
2、光伏支架的建模
根据光伏支架,将一串光伏组件的四个角点,定义为A、B、C、D点,建立坐标系o-xyz,其中z轴垂直于xoy平面。通过arcgis软件读取一串光伏组件的支架立柱的坐标点,通过支架立柱的坐标计算出A点坐标为(xA,yA,zA),同理得到B、C、D三点的坐标分别为(xB,yB,zB),(xC,yC,zC),(xD,yD,zD)。
任一取三个A、B、C、D中的三个点可求得光伏支架的平面T的方程,例如取A(Ax,Ay,Az)、B(Bx,By,Bz)、C(Cx,Cy,Cz)三个点通过式(8),求得平面方程的表达式为:
Txx+Tyy+Tzz+Tn=0 式(8)
将A(xA,yA,zA)、B(xB,yB,zB)、C(xC,yC,zC)三个点中任意一点带入式(8)即可获得Tn,如将A(xA,yA,zA)带入式(8)可得
Tn=-(TxAx+TyAy+TzAz) 式(10)
3、光伏地面(坡面)的建模
对于光伏电站地面的建模与步骤2中光伏支架平面的建模类似,通过arcgis软件,可以获取光伏地面的任意三个点E(Ex,Ey,Ez)、F(Fx,Fy,Fz)、H(Hx,Hy,Hz)的坐标,同理通过E、F、H三点可得到光伏地面的平面G(Gx,Gy,Gz)的表达为
Gxx+Gyy+Gzz+Gn=0 式(11)
同时,在arcgis中,光伏组件的地表面由等高线生成的三维的平面,记作G(Gx,Gy,Gz)。
4、光伏支架的投影计算
当太阳照射在物体上时,可以看作光线沿着光伏组件的边缘点T(Tx,Ty,Tz)投射在地平面G(Gx,Gy,Gz)上形成阴影点V(Vx,Vy,Vz)的结果。因此,分别计算光伏组件边缘点A(Ax,Ay,Az)、B(Bx,By,Bz)、C(Cx,Cy,Cz)、D(Dx,Dy,Dz)的投影点A′(A′x,A′y,A′z)、B′(B′x,By,B′z)、C′(C′x,C′y,C′z)、D′(D′x,D′y,D′z),再将A′、B′、C′、D′连线得到组件的阴影区域。投影点的V(Vx,Vy,Vz)和T(Tx,Ty,Tz)的关系如下:
V=T+k·S 式(13)
由于投影点位于地平面G上,则有:
G·V=0 式(14)
将式(14)带入式(13)得:
根据式(13)和(16)可得
将式(17)式(18)式(19)写成矩阵得:
N=Sx·Gx+Sy·Gy+Sz·Gz 式(21)
由式(20)可知,光伏组件的边缘点T(Tx,Ty,Tz)在地平面G(Gx,Gy,Gz)的投影点V(Vx,Vy,Vz)主要由太阳的方向向量S=[Sx,Sy,Sz]和地表的向量函数G(Gx,Gy,Gz)确定。其中,地表的方程函数是确定的,可视为常量,而太阳的方向向量随着太阳的移动轨迹变化而变化,在固定时刻,太阳的方向向量也可是为常量。因此,在固定时刻,光伏组件的边缘点T(Tx,Ty,Tz)在地表的投影点的V(Vx,Vy,Vz)的变换矩阵是确定的。其中变换矩阵M为:
5、光伏支架组件前后排遮挡模型
假设在坡面G(Gx,Gy,Gz)上,前后有两排阵列组件T1和T2,通过上述方法建模可计算的T1光伏支架组件的四个角点T1A,T1B,T1C,T1D在G(Gx,Gy,Gz)上的投影点为T1A′,T1B′,T1C′,T1D′,四个角点的在坡面G(Gx,Gy,Gz)投影区域记作集合T1′{T1′∈G(x,y,z)},T2光伏支架组件的四个角点T2A,T2B,T2C,T2D在G(Gx,Gy,Gz)上的投影点为T2A′,T2B′,T2C′,T2D′。投影区域记作集合T2′{T2′∈G(x,y,z)}。
令T1′∩T2′=T12′,点集T12′{T12′∈G(x,y,z)}。点集可以通过Arcgis软件获得。
令T12′(x,y)四个顶点为T12A′,T12B′,T12C′,T12D′,下面为找到四个顶点的步骤如下:
A.若存在P1、P2、P3、P4四个点互不相等的点,则遍历结束,P1、P2、P3、P4分别为T12′的顶点。
B.若存在P1、P2、P3、P4中有一个点同时满足以上四条件中两个,即只存在三个不同点,如P1、P2、P3,或者P1、P2、P4,P1、P3、P4,或者P2、P3、P4,则将这三个作为T12′的顶点;同时,继续根据连接三点中的与
C.若存在P1、P2、P3、P4中有两个点同时满足以上四条件中两个,即只存在两个不同点,将这两个点做为顶点;同时还需连接这两个点,通过连接的两点求得两点的连线L(x,y)=0,求出距离连线L(x,y)距离正值和负值最大的两点,并作为T12′的顶点。如图3所示,满足以上条件的两个不同点位为T12′A(x,y)、T12′C(x,y)(x,y),连接T12′A(x,y)、T12′C(x,y)(x,y),遍历点集T12′(x,y),求出距离连线L(x,y)距离正值和负值最大的两点T12′B(x,y)、T12′D(x,y),最终T12′A(x,y)、T12′C(x,y)、T12′B(x,y)、T12′D(x,y)记为为T12′的顶点。
对公式(20)进行求逆可得处坡面上的投影点V(Vx,Vy,Vz)计算光伏阵列组件产生投影点的点T(Tx,Ty,Tz)的公式,具体公式如下:
利用式(23)可求出T12′所有的顶点在对应后排阵列T2上的阴影点,连接所有的阴影点得到集合T12,则集合T12为后排T2光伏组件被前排光伏组件T1遮挡的区域。
我们根据所述遮挡区域即可计算出后排光伏发阵列电量损失。
以上描述了本发明的基本原理和主要的特征,说明书的描述只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。
Claims (3)
1.基于Arcgis的光伏阵列阴影计算方法,其特征在于,包括步骤:
建立三维坐标系o-xyz;
建立太阳法向量S(Sx,Sy,Sz)在三维坐标系中的数学模型;
通过Arcgis软件获取光伏支架的任意三个角点在三维坐标系中的坐标,根据获取的光伏支架的任意三个角点坐标求得光伏支架的平面方程T(Tx,Ty,Tz),建立光伏支架的平面在三维坐标系中的数学模型;
通过Arcgis软件获取坡面的任意三个点在三维坐标系中的坐标,根据获取的光伏地面的任意三个点的坐标求得坡面的平面方程G(Gx,Gy,Gz),建立坡面在三维坐标系中的数学模型;
根据太阳法向量的数学模型、坡面的的数学模型、光伏支架的平面的数学模型计算光伏组件边缘点T(Tx,Ty,Tz)在坡面上的投影点V(Vx,Vy,Vz)的变化矩阵M;
利用变化矩阵M计算前后两排光伏阵列的边缘点在坡面上的投影点;
根据前后两排光伏阵列的边缘点在坡面上的投影点确定前后两排光伏阵列的阴影交集;
若所述阴影交集不是空集,则获取所述阴影交集的顶点;
根据阴影交集的顶点以及变换矩阵M求逆得到前排阵列在后排阵列上的遮挡区域的顶点,连接所有的前排阵列在后排阵列的遮挡区域的顶点,从而得到前排阵列在后排阵列的遮挡区域。
3.如权利要求2所述的基于Arcgis的光伏阵列阴影计算方法,其特征在于,获取阴影交集T12′(x,y)的四个顶点的步骤包括:
A.若存在P1、P2、P3、P4四个点互不相等的点,则遍历结束,P1、P2、P3、P4分别为T12′的顶点;
C.若只存在两个不同点,将这两个点做为顶点;同时还需连接这两个点,通过连接的两点求得两点的连线L(x,y)=0,遍历点集T12′(x,y),求出距离连线L(x,y)距离正值和负值最大的两点,并作为T12′的顶点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611116780.1A CN106777634B (zh) | 2016-12-07 | 2016-12-07 | 基于Arcgis的光伏阵列阴影计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611116780.1A CN106777634B (zh) | 2016-12-07 | 2016-12-07 | 基于Arcgis的光伏阵列阴影计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106777634A CN106777634A (zh) | 2017-05-31 |
CN106777634B true CN106777634B (zh) | 2020-06-09 |
Family
ID=58882314
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611116780.1A Active CN106777634B (zh) | 2016-12-07 | 2016-12-07 | 基于Arcgis的光伏阵列阴影计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106777634B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108874739B (zh) * | 2018-06-11 | 2022-04-05 | 河海大学常州校区 | 光伏方阵间距遮挡下光伏组件辐照不均匀度计算方法 |
CN109684593B (zh) * | 2018-11-19 | 2023-06-06 | 南京国电南自新能源工程技术有限公司 | 山地光伏项目阵列间距计算方法及*** |
CN109815544B (zh) * | 2018-12-24 | 2024-01-09 | 中国电建集团山东电力建设第一工程有限公司 | 一种基于bim的屋顶光伏布置方法 |
CN110136257B (zh) * | 2019-04-30 | 2023-04-25 | 湖北省电力勘测设计院有限公司 | 一种山地光伏电站的自动化三维建模方法 |
CN110390715B (zh) * | 2019-09-11 | 2023-10-13 | 桂林理工大学 | 一种同时检测建筑物屋顶、建筑物墙体及地面阴影的方法 |
CN112560256B (zh) * | 2020-12-10 | 2023-04-07 | 中国电建集团贵州电力设计研究院有限公司 | 一种光伏组串最优间距计算***和计算方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103823927A (zh) * | 2014-02-12 | 2014-05-28 | 西安建筑科技大学 | 跟踪式光伏设备的阵列式排布方法 |
CN105160166A (zh) * | 2015-08-25 | 2015-12-16 | 许继集团有限公司 | 一种光伏阵列状态判断方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130298959A1 (en) * | 2012-05-09 | 2013-11-14 | Muhammed A. Alam | Shade-tolerant thin film photovoltaic panel |
CN106126887B (zh) * | 2016-06-16 | 2018-07-20 | 河海大学 | 一种基于参数辨识的多峰值光伏阵列最大功率点跟踪方法 |
-
2016
- 2016-12-07 CN CN201611116780.1A patent/CN106777634B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103823927A (zh) * | 2014-02-12 | 2014-05-28 | 西安建筑科技大学 | 跟踪式光伏设备的阵列式排布方法 |
CN105160166A (zh) * | 2015-08-25 | 2015-12-16 | 许继集团有限公司 | 一种光伏阵列状态判断方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106777634A (zh) | 2017-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106777634B (zh) | 基于Arcgis的光伏阵列阴影计算方法 | |
Nguyen et al. | Incorporating shading losses in solar photovoltaic potential assessment at the municipal scale | |
CN110031824A (zh) | 激光雷达联合标定方法及装置 | |
CN105022867B (zh) | 一种基于建筑信息模型的光伏设备自动铺设方法 | |
CN105931234A (zh) | 一种地面三维激光扫描点云与影像融合及配准的方法 | |
CN103823927A (zh) | 跟踪式光伏设备的阵列式排布方法 | |
CN113096250A (zh) | 基于无人机航拍影像序列的三维建筑模型库***构建方法 | |
CN113221222A (zh) | 建筑物屋顶光伏板排布方法 | |
KR20140029882A (ko) | 3차원 영상기반 신재생에너지 현장 평가 방법 및 평가 시스템 | |
Alam et al. | Shadow effect on photovoltaic potentiality analysis using 3D city models | |
CN113221356A (zh) | 复杂建筑物不规则屋顶光伏板排布***及优化方法 | |
CN113223151A (zh) | 建筑物屋顶光伏板排布*** | |
CN106452301A (zh) | 一种针对非正南向坡面的光伏组件安装设计方法 | |
CN105116470A (zh) | 基于地物tin模型的日照时数计算 | |
CN109684593A (zh) | 山地光伏项目阵列间距计算方法及*** | |
Gu et al. | Surveying and mapping of large-scale 3D digital topographic map based on oblique photography technology | |
CN104019800B (zh) | 大侧摆线阵ccd遥感图像对地定位的方法 | |
Deng et al. | Automatic true orthophoto generation based on three-dimensional building model using multiview urban aerial images | |
CN114187399A (zh) | 一种地外星体表面光照仿真影像生成方法和*** | |
CN103217145B (zh) | 一种火星dem制作和航带法空中三角测量方法 | |
CN114912171A (zh) | 一种基于建筑遮挡的日照辐射计算方法 | |
Meller | Analytically calculating shading in regular arrays of sun-pointing collectors | |
CN113420460A (zh) | 基于osg数据天际线的城市建筑限高快速分析方法及*** | |
CN113240725A (zh) | 复杂山体有效光照地面面积确定方法、装置和存储介质 | |
CN102034269B (zh) | 太阳直射光遥感数字图像的辐射变换虚拟仿真方法 |
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 |