CN110457832A - 一种稳态条件下双玻双面光伏组件各层温度的计算方法 - Google Patents
一种稳态条件下双玻双面光伏组件各层温度的计算方法 Download PDFInfo
- Publication number
- CN110457832A CN110457832A CN201910753453.4A CN201910753453A CN110457832A CN 110457832 A CN110457832 A CN 110457832A CN 201910753453 A CN201910753453 A CN 201910753453A CN 110457832 A CN110457832 A CN 110457832A
- Authority
- CN
- China
- Prior art keywords
- glass
- temperature
- double
- photovoltaic module
- plate
- 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
-
- 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
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
Landscapes
- Photovoltaic Devices (AREA)
Abstract
本发明公开了一种稳态条件下双玻双面光伏组件各层温度的计算方法,根据能量守衡定律建立双玻双面光伏组件每一层的稳态热平衡公式;将所述双玻双面光伏组件每一层的稳态热平衡公式进行整合,得到以玻璃盖板温度,上层EVA温度,电池片温度,下层EVA温度和玻璃背板温度为未知量的方程组;根据所得到的方程组编写程序求解双玻双面光伏组件玻璃盖板温度,上层EVA温度,电池片温度,下层EVA温度和玻璃背板温度。通过环境条件参数准确计算组件各层材料温度,有利于评估组件在不同使用环境下的可靠性能与耐候性能。
Description
技术领域
本发明公开了一种稳态条件下双玻双面光伏组件各层温度的计算方法,属于光伏组件发电技术领域。
背景技术
太阳能电池片的转换效率会随着工作温度的升高而降低,降低电池片产生的热量是改善组件性能的关键因素。双玻双面光伏组件成本与普通组件越来越接近,由于发电增益明显,尤其在高反射地区增益达20%以上,其相关技术引起了产业界与研究者们的广泛关注。不同电池片技术及组件设计都会对组件的散热性能有一定的影响,双玻双面组件主要是通过外界环境将太阳光线反射到组件背面从而加以利用,怎样准确计算双玻双面组件工作温度对评估双面组件设计与应用可行性十分重要。
评估双面组件工作温度最重要的是计算中心电池片的温度,因为电池片的温度直接影响组件的发电功率,而现有的多数方法均忽略组件内部温度差异,将组件看成温度一致的整体。
发明内容
本发明的目的在于提供一种稳态条件下双玻双面光伏组件各层温度的计算方法,通过将组件各层进行拆分,建立各层的能量守恒方程,进而准确计算出双面组件包括电池在内的各层工作温度。
为达到上述目的,本发明采用的技术方案如下:
一种稳态条件下双玻双面光伏组件各层温度的计算方法,包括:
根据稳态条件下双玻双面光伏组件的能量守衡定律,建立双玻双面光伏组件每一层的稳态热平衡公式;
将所述双玻双面光伏组件每一层的稳态热平衡公式进行整合,得到以玻璃盖板温度,上层EVA温度,电池片温度,下层EVA温度和玻璃背板温度为未知量的方程组;
根据所得到的方程组编写程序求解双玻双面光伏组件玻璃盖板温度,上层EVA温度,电池片温度,下层EVA温度和玻璃背板温度。
进一步的,根据稳态条件下双玻双面光伏组件的能量守衡定律,建立双玻双面光伏组件每一层的稳态热平衡公式,包括:
玻璃盖板的稳态热平衡公式为:
I+Rcond=Rfs+Rfg+Rfa+Iρf+Iτf
上层EVA的稳态热平衡公式为:
电池片的稳态热平衡公式为:
下层EVA的稳态热平衡公式为:
玻璃背板的稳态热平衡公式为:
其中,I为太阳辐射值,Rcond为玻璃盖板与上层EVA之间的热传导,ρf为玻璃盖板的反射率,τf为玻璃盖板的透射率,Rfs为玻璃盖板与天空之间的热辐射,Rfg为玻璃盖板与地面之间的热辐射,Rfa为玻璃盖板与空气之间的热对流,Tc为电池片温度,TE1为上层EVA温度,R2为上层EVA与电池片之间的热阻,ρE1为上层EVA的反射率,τE1为上层EVA的透射率,R1为玻璃盖板与上层EVA之间的热阻,Tf为玻璃盖板温度,Gr1为双玻双面光伏组件背面接收到的总辐射,P为双玻双面光伏组件输出的功率,A为双玻双面光伏组件的表面积,TE2为下层EVA温度,ρc为电池片的反射率,τb为玻璃背板的透射率,τE2为下层EVA的透射率,R3为电池片与下层EVA之间的热阻,Tb为玻璃背板温度,R4为下层EVA与玻璃背板之间的热阻,ρE2为下层EVA的反射率,ρb为玻璃背板的反射率,Rbs为玻璃背板与天空的辐射换热,Rbg为玻璃背板与地面的辐射换热,Rba为玻璃背板与地面的对流换热。
进一步的,所述玻璃盖板与上层EVA之间的热传导计算如下:
进一步的,
所述玻璃盖板与天空之间的热辐射计算如下:
其中,σ为波兹曼常数,εf为玻璃盖板的发射率,εs为天空的发射率,Tf为玻璃盖板的温度,Ts为天空的温度,
Ffs为玻璃盖板与天空之间的角系数:
其中,θ为双玻双面光伏组件的安装角度;
所述玻璃盖板与地面之间的热辐射,计算如下:
其中,εg为地面的发射率,Tg为地面的温度,
Ffg为玻璃盖板与地面之间的角系数:
所述玻璃盖板与空气之间的热对流,计算如下:
Rfa=hfa(Tf-Ta)
其中,hfa为玻璃盖板与空气之间的热对流系数,Ta为环境温度。
进一步的,
所述玻璃盖板与上层EVA之间的热阻,计算如下:
其中,δf为玻璃盖板的厚度,δE1为上层EVA的厚度,λf为玻璃盖板的导热系数,λE1为上层EVA的导热系数;
所述上层EVA与电池片之间的热阻,计算如下:
其中,δc为电池片的厚度,λc为电池片的导热系数;
所述电池片与下层EVA之间的热阻,计算如下:
其中,δE2为下层EVA的厚度,λE2为下层EVA的导热系数;
所述下层EVA与玻璃背板之间的热阻,计算如下:
其中,δb为玻璃背板的厚度,λb为玻璃背板的导热系数。
进一步的,通过视角系数法计算双玻双面光伏组件背面接收到的总辐射,计算如下:
每排双玻双面光伏组件背面收到的总辐射包括:
地面反射:
其中,G1为地面对每排双面电池片背面的反射辐射度,ρg为地面反射率,DHI为水平散射辐照度,A1为每排双面电池片地面阴影区域面积,An为每排双面电池片的面积,F1为地面阴影区对每排双面电池片背面的视角系数,GHI为水平总辐照度,A2为每排双面电池片地面非阴影区域面积,F2为地面非阴影区对每排双面电池片背面的视角系数;
天空散射:
其中,G2表示天空散射对每排双面电池片背面的辐射度,θ为双玻双面光伏组件的安装倾角;
以及后排组件反射天空散射:
其中,G3表示后排双面电池片对前排双面电池片背面的辐射度,Gm为后排双面电池片表面接收到的辐照度,ρm为双玻双面光伏组件的反射率,Am为后排双面电池片面积,F3为后排双面电池片对前排双面电池片背面的视角系数;
每排双面电池片背面接收到的总辐照度为:Gr1=G1+G2+G3;
所有排数的双面电池片背面接收到的总辐照度之和即为整个双玻双面光伏组件背面收到的总辐射。
进一步的,
所述玻璃背板与天空的辐射换热,计算如下:
其中,εb为玻璃背板的发射率,
Fbs为玻璃背板与天空的角系数:
玻璃背板与地面的辐射换热,计算如下:
其中,Fbg为背板与地面的角系数:
玻璃背板与地面的对流换热,计算如下:
Rba=hba(Tb-Ta)
其中,hba为背板与空气的对流换热系数。
进一步的,所述天空的温度计算如下:
Ts=0.0552Ta 1.5。
进一步的,所述将所述双玻双面光伏组件每一层的稳态热平衡公式进行整合,得到以玻璃盖板温度,上层EVA温度,电池片温度,下层EVA温度和玻璃背板温度为未知量的方程组,包括:
进一步的,所述根据所得到的方程组编写程序求解双玻双面光伏组件玻璃盖板温度,上层EVA温度,电池片温度,下层EVA温度和玻璃背板温度,包括:
设:A1=σεf,
D2=Iτf(τE1+ρE1-1)
D3=Gr1τb(τE2+ρE2-1)
A5=-σεb,
得到:
A2Tf+B2TE1+C2Tc=D2
A3TE1+B3TE2+C3Tc=D3
A4Tc+B4TE2+C4Tb=D4
运用matlab软件行求解以上的五元四次方程组,得到玻璃盖板温度Tf、上层EVA温度TE1、电池片温度Tc、下层EVA温度TE2和玻璃背板温度Tb;
Matlab程序如下:
本发明的优点在于,本发明运用视角系数法计算双面组件背面辐照强度,综合考虑双面组件正、反面的发电功率、热功率以及外部散热,并将组件各层进行拆分,建立各层的能量守恒方程,可达到根据环境因素、组件自身条件则能够准确计算出双面组件包括玻璃、电池、EVA、背板在内的各层的工作温度;通过环境条件参数准确计算组件各层材料温度,有利于评估组件在不同使用环境下的可靠性能与耐候性能。
通过该方法计算得到太阳电池的准确温度,根据电池温度系数可以计算组件实际工作条件下因为温度引起的功率与发电量损失,准确预测组件实际发电性能,为电站经济性分析,投资可行性评估等提供支撑。
附图说明
图1为本发明双玻双面光伏组件分层工作温度计算的原理图;
图2为视角系数求解示意图;
图3为双玻双面光伏组件玻璃盖板的能量守恒示意图;
图4为双玻双面光伏组件上层EVA的能量守恒示意图;
图5为双玻双面光伏组件中电池片的能量守恒示意图;
图6为双玻双面光伏组件下层EVA的能量守恒示意图;
图7为双玻双面光伏组件玻璃背板的能量守恒示意图;
图8为本发明实施例中双玻双面光伏组件安装示意图。
具体实施方式
下面对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
本发明通过视角系数法,结合典型双面电池组件应用场景中地面反射率与组件间距,计算双面电池背面辐照;根据背面辐照,电池背面效率,吸收率与发射率,计算电池背面热功率;根据太阳辐照强度,电池正面效率,吸收率,反射率计算电池正面热功率;电池正、反面热功率和作为双玻双面光伏组件内热源功率,设定材料导热系数、对流换热系数、发射率,环境温度等,根据内部产生热量等于外部散热量的热平衡方程,建立热稳态模型,列出计算得到双面组件各层工作温度。
参见图1,本发明的基本原理如下:
组件内部热源:电池正(背)面功率=正(背)面辐照-反射的辐照-发电功率。
正(背)面电池对光的吸收率与转化率直接影响发电功率。
组件外部散热:与环境的对流换热、辐射换热。
稳态时,组件内部热源与外部散失的热量相等。
组件整体的能量守恒方程为:
等式左边为组件内部热源,等式右边为外部散热。
其中,I为组件正面接收到的太阳辐射值,Gr1为组件背面接收到的总辐射,ρf为玻璃盖板的反射率,ρb为玻璃背板的反射率,Rfs为玻璃盖板与天空之间的热辐射,Rfg为玻璃盖板与地面之间的热辐射,Rfa为玻璃盖板与空气之间的热对流,Rbs为玻璃背板与天空的辐射换热,Rbg为玻璃背板与地面的辐射换热,Rba为玻璃背板与地面的对流换热,A为组件的表面积,P包含两个部分,分别为组件正面发电功率和背面发电功率。
组件的正面热功率为:
组件的背面热功率为:
P1=ηI
P2=γGr1
其中,P1为组件正面发电功率,P2为组件反面发电功率,η和γ均为光电转化效率,与电池吸收率有关。
本发明实施林的稳态条件下双玻双面光伏组件工作温度的计算方法,具体实施过程如下:
1、通过视角系数法计算双玻双面光伏组件背面辐照
双玻双面光伏组件背面辐照度模型分为视图因子模型和光线跟踪模拟两类。与射线跟踪仿真相比,视图因子模型的计算要求较低,所需参数较少,更易于计算。视图因子,也称为形状和配置因子,量化从一个表面反射到接收表面的辐射能的百分比。参见图2,用F1→2表示平面1发出的辐射能中落到平面2的百分数,计算公式为:
其中,r是平面1与平面2之间的距离,A1表示平面1,A2表示平面2,与为平面1与平面2上的单元面积,θ1为平面1法向量与连接和的线之间的角度,θ2为平面2法向量与连接和的线之间的角度。
到达双玻双面光伏组件背面的辐射来自三个区域,即地面反射,天空散射以及后排组件反射天空散射,由于双玻双面光伏组件多为阵列安装,为了便于计算,本发明中,忽略水平方向上的影响,认为双玻双面光伏组件背面辐照的不均匀度只发生在垂直方向上,则本发明以水平方向上每排电池片为一个单元进行建模计算,把双玻双面光伏组件进行多排划分。然后计算双玻双面光伏组件背面辐射:
地面反射:包括阴影区域和非阴影区域对每排双面电池片背面的反射辐射,表示为:
其中,G1为地面对每排双面电池片背面的反射辐射度,单位(W/m2),ρg为地面反射率,DHI为水平散射辐照度,单位(W/m2),A1为每排双面电池片地面阴影区域面积,单位(m2);An为每排双面电池片的面积,单位(m2),F1为地面阴影区对每排双面电池片背面的视角系数,GHI为水平总辐照度,单位(W/m2),A2为每排双面电池片地面非阴影区域面积,单位(W/m2),F2为地面非阴影区对每排双面电池片背面的视角系数。
天空散射:目前,普遍用来计算倾斜面空中散射辐射分为天空各向同性模型和天空各向异性模型,由于天空各向同性模型容易理解,计算简单,在本发明中采用天空各向同性模型计算:
其中,G2表示天空散射对每排双面电池片背面的辐射度,单位(W/m2),θ为双玻双面光伏组件的安装倾角。
后排组件反射天空散射:
其中,G3表示后排双面电池片对前排双面电池片背面的辐射度,单位(W/m2),Gm为后排双面电池片表面接收到的辐照度,单位(W/m2),ρm为双玻双面光伏组件的反射率,Am为后排双面电池片面积,单位(m2),F3为后排双面电池片对前排双面电池片背面的视角系数。
因此,每排双面电池片背面接收到的总辐照度为:
Gr1=G1+G2+G3 (5)
所有排数的双面电池片背面接收到的总辐照度之和即为整个双玻双面光伏组件背面接收到的总辐照度。
然而,r,θ1和θ2是取决于微分区域的变量,求解方程(1)十分复杂。由于光伏组件多为阵列安装,因此我们可以忽略水平方向上的影响,把光伏组件,地面阴影区和地面非阴影区假设为一个无限长(水平方向)的平面,分别计算视角系数F1,F2,F3。
2、建立双玻双面光伏组件稳态传热模型
根据双玻双面光伏组件的传热机理将组件分为五个部分,由组件正面至背面依次为:玻璃盖板、上层EVA、电池片、下层EVA和玻璃背板。分析各层之间的传热机理,假设正、背面入射的太阳光均不穿透电池片且不考虑时间对温度的影响,即忽略温度的滞后性。基于这五个部分之间的热对流,热辐射和热传导,根据能量守衡定律,建立双玻双面光伏组件每一层的稳态热平衡公式。
2.1玻璃盖板的能量守恒模型
如图3所示,当太阳辐射照射到玻璃盖板面,会有一部分辐射被反射,一部分透射过玻璃进入下层EVA,同时玻璃盖板会和天空、地面存在一定的辐射换热,会和空气产生对流换热,也会和上层EVA存在热传导。根据能量守恒定律,进入玻璃盖板的能量等于出去的能量,可以列出玻璃盖板的能量守恒公式:
I+Rcond=Rfs+Rfg+Rfa+Iρf+Iτf (6)
其中,I为太阳辐射值,单位(W/m2),ρf为玻璃盖板的反射率,τf为玻璃盖板的透射率,Rcond为玻璃盖板与上层EVA之间的热传导,单位(W/m2)。
其中,TE1为上层EVA的温度,单位(K),Tf为玻璃盖板的温度,单位(K),R1为玻璃盖板与上层EVA之间的热阻,单位(K/W):
其中,δf为玻璃盖板的厚度,单位(m),δE1为上层EVA的厚度,单位(m),λf为玻璃盖板的导热系数,单位(W/m·K),λE1为上层EVA的导热系数,单位(W/m·K)。
Rfs为玻璃盖板与天空之间的热辐射,单位(W/m2):
其中,σ为波兹曼常数,εf为玻璃盖板的发射率,εs为天空的发射率,Tf为玻璃盖板的温度,单位(K),Ts为天空的温度,单位(K),
Ffs为玻璃盖板与天空之间的角系数:
其中,θ为双玻双面光伏组件的安装角度。
Rfg为玻璃盖板与地面之间的热辐射,单位(W/m2):
其中,εg为地面的发射率,Tg为地面的温度,单位(K),
Ffg为玻璃盖板与地面之间的角系数:
Rfa为玻璃盖板与空气之间的热对流,单位(W/m2):
Rfa=hfa(Tf-Ta) (13)
其中,hfa为玻璃盖板与空气之间的热对流系数,单位(W/m2·K),Ta为环境温度,单位(℃)。
2.2上层EVA的能量守恒模型
如图4所示,对于上层EVA,进入的能量有到达它的辐射,即透射过玻璃盖板的辐射以及与电池片的热传导,出去的能量有被上层EVA反射的辐射,与玻璃盖板的热传导,透过上层EVA的辐射。根据能量守恒定律,进入上层EVA的能量等于出去的能量,可以列出上层EVA的能量守恒公式:
其中,Tc为电池片的温度,单位(K),ρE1为上层EVA的反射率,τE1为上层EVA的透射率,TE1为上层EVA的温度,
R2为上层EVA与电池片之间的热阻,单位(K/W):
其中,δc为电池片的厚度,单位(m),λc为电池片的导热系数,单位(W/m·K)。
2.3电池片的能量守恒模型
如图5所示,对于电池片,进入的能量有从正、背面入射到达它的辐射,即透射过EVA的辐射,出去的能量有输出的电能、被反射的辐射、与上层EVA的热传导以及与下层EVA的热传导,根据能量守恒定律,可以列出电池片的能量守恒公式:
其中,TE2为下层EVA的温度,单位(K),ρc为电池片的反射率,τb为玻璃背板的透射率,τE2为下层EVA的透射率,Gr1为从双玻双面光伏组件背面收到的总辐射,单位(K/W),P为电池片输出的功率,单位(W),A为双玻双面光伏组件的表面积,单位(m2),
R3为电池片与下层EVA之间的热阻,单位(K/W):
其中,δE2为下层EVA的厚度,单位(m),λE2为下层EVA的导热系数,单位(W/m·K)。
2.4下层EVA的能量守恒模型
如图6所示,对于下层EVA,进入的能量有与电池片的热传导和从背面入射到达下层EVA的辐射,出去的能量有与玻璃背板的热传导、反射的辐射和透射的辐射,根据能量守恒定律,可以列出下层EVA的能量守恒公式:
其中,ρE2为下层EVA的反射率,Tb为背板的温度,单位(K),
R4为下层EVA与背板之间的热阻,单位(K/W):
其中,δb为背板的厚度,单位(m),λb为背板的导热系数,单位(W/m·K)。
2.5玻璃背板的能量守恒模型
如图7所示,对于玻璃背板,进入的能量有与下层EVA的热传导和从背板入射的辐射,出去的能量有与天空的辐射换热、与地面的辐射换热、与空气的对流换热、被反射的太阳辐照和透射的太阳辐射,根据能量守恒定律,可以列出玻璃背板的能量守恒公式:
其中,ρb为玻璃背板的反射率,
Rbs为玻璃背板与天空的辐射换热,单位(W/m2):
其中,εb为背板的发射率,
Fbs为背板与天空的角系数:
Rbg为玻璃背板与地面的辐射换热,单位(W/m2),
其中,Fbg为背板与地面的角系数:
Rba为玻璃背板与地面的对流换热,单位(W/m2):
Rba=hba(Tb-Ta) (25)
其中,hba为背板与空气的对流换热系数,单位(W/m2·K)。
将式(6)-(25)整理合并后得到双玻双面光伏组件稳态传热模型式(26)-(30):
3、根据双玻双面光伏组件稳态传热模型计算双玻双面光伏组件各层工作温度
设:A1=σεf,
D2=Iτf(τE1+ρE1-1)
D3=Gr1τb(τE2+ρE2-1)
A5=-σεb,
以上系数和常数项均可根据实际的组件安装情况、辐照强度、风速和组件阵列数目计算得到,得到:
A2Tf+B2TE1+C2Tc=D2 (32)
A3TE1+B3TE2+C3Tc=D3 (33)
A4Tc+B4TE2+C4Tb=D4 (34)
运用matlab软件行求解以上的五元四次方程组,可得到玻璃盖板温度Tf、上层EVA温度TE1、电池片温度Tc、下层EVA温度TE2和玻璃背板温度Tb。
Matlab程序如下:
实施例
双玻双面光伏组件在户外工作条件下的二维结构示意图如图8所示,组件安装倾角θ为37.5°,组件末端距离地面的距离0.6m。以环境参数:总辐照强度I800W/m2,环境温度Ta20℃,地面温度Tg为30℃,风速1m/s的工作环境为例,计算组件各层材料的工作温度。
假设每排仅放置单一的光伏组件(仅用于简化计算,不影响实施过程),每排之间间隔1m。
组件的基本参数如下:组件长度1.954m、组件宽度0.986m、组件表面积A为1.927m2、玻璃盖板的发射率εf为0.85、玻璃背板的发射率εb为0.85、天空的发射率εs为0.85、地面的发射率εg为0.88,组件标称输出功率P为320W。
采用视角系数法计算背面的辐照强度:设置地面反射率ρg为0.5,水平散射辐照强度为水平总辐射强度的15%(DHI为15%*800=120)阴影区面积A1与非阴影区面积A2均设为1m2,后排双面电池片表面接收到的辐照度为800W/m2。则地面阴影区与非阴影区对玻璃背板的视图因子F1、F2均为0.53。后派组件与前排组件除安装位置外,其余均一致,前后排间距1.5m。组件整体的反射率ρm为0.165,后排双面电池对前排双面电池的视角系数根据公式1计算得0.85。则:
Gr1=126.52+7.9+112.33=246.75
组件各层的导热系数以及厚度设置如表1所示。
表1各层参数设置
根据上表可得:
玻璃盖板与上层EVA之间的热阻R1为:2.38×10-3K/W;
上层EVA与电池片之间的热阻R2为:7.85×10-4K/W;
电池片与下层EVA之间的热阻R3为:7.31×10-4K/W;
下层EVA与玻璃背板之间的热阻R4为:1.73×10-3K/W。
光伏组件的对流换热可以分为两种形式,一种是强制对流换热,一种是自然对流换热。强制对流换热主要是因为组件上表面和下表面有风经过从而引起的空气流动产生的热量的传动,自然对流换热主要是因为组件上下表面温度分布不均匀从而引起的热量的传动。用自然对流换热系数和强制对流换热系数两者中的最大值作为组件对流换热系数。在风速大于1m/s时,强制对流换热系数大于自然对流换热系数,而风速小于1m/s时,自然对流换热系数大于强制对流换热系数。此次将强制对流换热系数作为组件上下表面的对流换热系数。光伏组件的强制对流换热公式如式(36)
hfa=5.7+3.8v (36)
玻璃背板的对流换热系数为玻璃盖板面的一半,风速v=1m/s,则组件上下表面的对流换热系数如表2所示:
表2对流换热系数
光伏组件的辐射换热主要分为两个部分,一个是组件上表面即玻璃面与天空之间的辐射换热,另一个就是组件下表面即背板面与地面之间的辐射换热。天空的温度计算公式为:
Ts=0.0552Ta 1.5 (37)
则当环境温度为20℃时,天空温度为5℃。
在计算组件各层实际吸收的太阳辐射能时,需考虑材料的吸收率、透射率和反射率,其具体数值如表3所示,
表3各层的吸收率、透射率与反射率表
本发明实施例中,玻璃背板的透射率和反射率取一样的,上层EVA和下层EVA的透射率和反射率取一样的。
将以上数值及计算结果带入式(31)-(35),则可以方程的系数如表4所示:
表4方程组系数
运用matlab求解solve函数即可得到组件各层温度。
玻璃盖板温度Tf为316.06K、上层EVA温度TE1为317.00K、电池片温度Tc为317.29K、下层EVA温度TE2为317.16K、玻璃背板温度Tb为316.86K。
该实施例仅做参考,具体实施过程中的参数细节以实际情况为准。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (10)
1.一种稳态条件下双玻双面光伏组件各层温度的计算方法,其特征在于,包括:
根据稳态条件下双玻双面光伏组件的能量守衡定律,建立双玻双面光伏组件每一层的稳态热平衡公式;
将所述双玻双面光伏组件每一层的稳态热平衡公式进行整合,得到以玻璃盖板温度,上层EVA温度,电池片温度,下层EVA温度和玻璃背板温度为未知量的方程组;
根据所得到的方程组编写程序求解双玻双面光伏组件玻璃盖板温度,上层EVA温度,电池片温度,下层EVA温度和玻璃背板温度。
2.根据权利要求1所述的一种稳态条件下双玻双面光伏组件各层温度的计算方法,其特征在于,根据稳态条件下双玻双面光伏组件的能量守衡定律,建立双玻双面光伏组件每一层的稳态热平衡公式,包括:
玻璃盖板的稳态热平衡公式为:
I+Rcond=Rfs+Rfg+Rfa+Iρf+Iτf
上层EVA的稳态热平衡公式为:
电池片的稳态热平衡公式为:
下层EVA的稳态热平衡公式为:
玻璃背板的稳态热平衡公式为:
其中,I为太阳辐射值,Rcond为玻璃盖板与上层EVA之间的热传导,ρf为玻璃盖板的反射率,τf为玻璃盖板的透射率,Rfs为玻璃盖板与天空之间的热辐射,Rfg为玻璃盖板与地面之间的热辐射,Rfa为玻璃盖板与空气之间的热对流,Tc为电池片温度,TE1为上层EVA温度,R2为上层EVA与电池片之间的热阻,ρE1为上层EVA的反射率,τE1为上层EVA的透射率,R1为玻璃盖板与上层EVA之间的热阻,Tf为玻璃盖板温度,Gr1为双玻双面光伏组件背面接收到的总辐射,P为双玻双面光伏组件输出的功率,A为双玻双面光伏组件的表面积,TE2为下层EVA温度,ρc为电池片的反射率,τb为玻璃背板的透射率,τE2为下层EVA的透射率,R3为电池片与下层EVA之间的热阻,Tb为玻璃背板温度,R4为下层EVA与玻璃背板之间的热阻,ρE2为下层EVA的反射率,pb为玻璃背板的反射率,Rbs为玻璃背板与天空的辐射换热,Rbg为玻璃背板与地面的辐射换热,Rba为玻璃背板与地面的对流换热。
3.根据权利要求2所述的一种稳态条件下双玻双面光伏组件各层温度的计算方法,其特征在于,所述玻璃盖板与上层EVA之间的热传导计算如下:
4.根据权利要求2所述的一种稳态条件下双玻双面光伏组件各层温度的计算方法,其特征在于,
所述玻璃盖板与天空之间的热辐射计算如下:
其中,σ为波兹曼常数,εf为玻璃盖板的发射率,εs为天空的发射率,Tf为玻璃盖板的温度,Ts为天空的温度,
Ffs为玻璃盖板与天空之间的角系数:
其中,θ为双玻双面光伏组件的安装角度;
所述玻璃盖板与地面之间的热辐射,计算如下:
其中,εg为地面的发射率,Tg为地面的温度,
Ffg为玻璃盖板与地面之间的角系数:
所述玻璃盖板与空气之间的热对流,计算如下:
Rfa=hfa(Tf-Ta)
其中,hfa为玻璃盖板与空气之间的热对流系数,Ta为环境温度。
5.根据权利要求2所述的一种稳态条件下双玻双面光伏组件各层温度的计算方法,其特征在于,
所述玻璃盖板与上层EVA之间的热阻,计算如下:
其中,δf为玻璃盖板的厚度,δE1为上层EVA的厚度,λf为玻璃盖板的导热系数,λE1为上层EVA的导热系数;
所述上层EVA与电池片之间的热阻,计算如下:
其中,δc为电池片的厚度,λc为电池片的导热系数;
所述电池片与下层EVA之间的热阻,计算如下:
其中,δE2为下层EVA的厚度,λE2为下层EVA的导热系数;
所述下层EVA与玻璃背板之间的热阻,计算如下:
其中,δb为玻璃背板的厚度,λb为玻璃背板的导热系数。
6.根据权利要求2所述的一种稳态条件下双玻双面光伏组件各层温度的计算方法,其特征在于,通过视角系数法计算双玻双面光伏组件背面接收到的总辐射,计算如下:
每排双玻双面光伏组件背面收到的总辐射包括:
地面反射:
其中,G1为地面对每排双面电池片背面的反射辐射度,ρg为地面反射率,DHI为水平散射辐照度,A1为每排双面电池片地面阴影区域面积,An为每排双面电池片的面积,F1为地面阴影区对每排双面电池片背面的视角系数,GHI为水平总辐照度,A2为每排双面电池片地面非阴影区域面积,F2为地面非阴影区对每排双面电池片背面的视角系数;
天空散射:
其中,G2表示天空散射对每排双面电池片背面的辐射度,θ为双玻双面光伏组件的安装倾角;
以及后排组件反射天空散射:
其中,G3表示后排双面电池片对前排双面电池片背面的辐射度,Gm为后排双面电池片表面接收到的辐照度,ρm为双玻双面光伏组件的反射率,Am为后排双面电池片面积,F3为后排双面电池片对前排双面电池片背面的视角系数;
每排双面电池片背面接收到的总辐照度为:Gri=G1+G2+G3;
所有排数的双面电池片背面接收到的总辐照度之和即为整个双玻双面光伏组件背面收到的总辐射。
7.根据权利要求4所述的一种稳态条件下双玻双面光伏组件各层温度的计算方法,其特征在于,
所述玻璃背板与天空的辐射换热,计算如下:
其中,εb为玻璃背板的发射率,
Fbs为玻璃背板与天空的角系数:
玻璃背板与地面的辐射换热,计算如下:
其中,Tbg为背板与地面的角系数:
玻璃背板与地面的对流换热,计算如下:
Rba=hDa(Tb-Ta)
其中,hba为背板与空气的对流换热系数。
8.根据权利要求4所述的一种稳态条件下双玻双面光伏组件各层温度的计算方法,其特征在于,所述天空的温度计算如下:
Ts=0.0552Ta 1.5。
9.根据权利要求4所述的一种稳态条件下双玻双面光伏组件各层温度的计算方法,其特征在于,所述将所述双玻双面光伏组件每一层的稳态热平衡公式进行整合,得到以玻璃盖板温度,上层EVA温度,电池片温度,下层EVA温度和玻璃背板温度为未知量的方程组,包括:
10.根据权利要求9所述的一种稳态条件下双玻双面光伏组件各层温度的计算方法,其特征在于,所述根据所得到的方程组编写程序求解双玻双面光伏组件玻璃盖板温度,上层EVA温度,电池片温度,下层EVA温度和玻璃背板温度,包括:
设:A1=σεf,
D2=Iτf(τE1+ρE1-1)
D3=Gr1τb(τE2+ρE2-1)
A5=-σεb,
得到:
A2Tf+B2TE1+c2Tc=D2
A3TE1+B3TE2+C3Tc=D3
A4Tc+B4TE2+C4Tb=D4
运用matlab软件行求解以上的五元四次方程组,得到玻璃盖板温度Tf、上层EVA温度TE1、电池片温度Tc、下层EVA温度TE2和玻璃背板温度Tb;
Matlab程序如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910753453.4A CN110457832B (zh) | 2019-08-15 | 2019-08-15 | 一种稳态条件下双玻双面光伏组件各层温度的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910753453.4A CN110457832B (zh) | 2019-08-15 | 2019-08-15 | 一种稳态条件下双玻双面光伏组件各层温度的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110457832A true CN110457832A (zh) | 2019-11-15 |
CN110457832B CN110457832B (zh) | 2022-08-16 |
Family
ID=68486771
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910753453.4A Active CN110457832B (zh) | 2019-08-15 | 2019-08-15 | 一种稳态条件下双玻双面光伏组件各层温度的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110457832B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112541152A (zh) * | 2020-10-10 | 2021-03-23 | 天合光能股份有限公司 | 双面光伏组件工作温度的计算方法 |
CN112671336A (zh) * | 2020-12-16 | 2021-04-16 | 广东电网有限责任公司电力调度控制中心 | 检测光伏组件工作温度异常的方法、装置和计算机设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104467663A (zh) * | 2014-11-21 | 2015-03-25 | 河海大学常州校区 | 一种检测光伏热斑损害的太阳电池漏电流的系统及方法 |
CN106546626A (zh) * | 2016-11-04 | 2017-03-29 | 无锡市产品质量监督检验中心 | 一种求解光伏组件温度系数的电池片温度修正方法 |
CN107181463A (zh) * | 2017-06-20 | 2017-09-19 | 河海大学常州校区 | 一种精确预测太阳能电池组件工作温度的方法 |
-
2019
- 2019-08-15 CN CN201910753453.4A patent/CN110457832B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104467663A (zh) * | 2014-11-21 | 2015-03-25 | 河海大学常州校区 | 一种检测光伏热斑损害的太阳电池漏电流的系统及方法 |
CN106546626A (zh) * | 2016-11-04 | 2017-03-29 | 无锡市产品质量监督检验中心 | 一种求解光伏组件温度系数的电池片温度修正方法 |
CN107181463A (zh) * | 2017-06-20 | 2017-09-19 | 河海大学常州校区 | 一种精确预测太阳能电池组件工作温度的方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112541152A (zh) * | 2020-10-10 | 2021-03-23 | 天合光能股份有限公司 | 双面光伏组件工作温度的计算方法 |
CN112541152B (zh) * | 2020-10-10 | 2024-03-22 | 天合光能股份有限公司 | 双面光伏组件工作温度的计算方法 |
CN112671336A (zh) * | 2020-12-16 | 2021-04-16 | 广东电网有限责任公司电力调度控制中心 | 检测光伏组件工作温度异常的方法、装置和计算机设备 |
CN112671336B (zh) * | 2020-12-16 | 2022-01-25 | 广东电网有限责任公司电力调度控制中心 | 检测光伏组件工作温度异常的方法、装置和计算机设备 |
Also Published As
Publication number | Publication date |
---|---|
CN110457832B (zh) | 2022-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Aste et al. | Design, modeling and performance monitoring of a photovoltaic–thermal (PVT) water collector | |
Sellami et al. | Optical efficiency study of PV crossed compound parabolic concentrator | |
Kaplani et al. | Thermal modelling and experimental assessment of the dependence of PV module temperature on wind velocity and direction, module orientation and inclination | |
Corbin et al. | Experimental and numerical investigation on thermal and electrical performance of a building integrated photovoltaic–thermal collector system | |
Rosell et al. | Design and simulation of a low concentrating photovoltaic/thermal system | |
Jones et al. | A thermal model for photovoltaic systems | |
Yousef et al. | Performance investigation of low–Concentration photovoltaic systems under hot and arid conditions: Experimental and numerical results | |
Agrawal et al. | Optimizing the energy and exergy of building integrated photovoltaic thermal (BIPVT) systems under cold climatic conditions | |
Tripanagnostopoulos et al. | Hybrid photovoltaic/thermal solar systems | |
Feng et al. | A novel solar multifunctional PV/T/D system for green building roofs | |
Ji et al. | Thermal characteristics of a building-integrated dual-function solar collector in water heating mode with natural circulation | |
Chaabane et al. | Performance evaluation of concentrating solar photovoltaic and photovoltaic/thermal systems | |
Lu et al. | Energy analysis of silicon solar cell modules based on an optical model for arbitrary layers | |
Zhang et al. | The mathematical and experimental analysis on the steady-state operating temperature of bifacial photovoltaic modules | |
Tina et al. | Assessment of the electrical and thermal performances of building integrated bifacial photovoltaic modules | |
Lim et al. | Design and evaluation of passive concentrator and reflector systems for bifacial solar panel on a highly cloudy region–A case study in Malaysia | |
Singh et al. | Experimental investigations into low concentrating line axis solar concentrators for CPV applications | |
Simonetti et al. | Development and validation of a comprehensive dynamic mathematical model for hybrid PV/T solar collectors | |
CN105260528B (zh) | 光伏组件的输出特性计算方法及系统 | |
CN110457832A (zh) | 一种稳态条件下双玻双面光伏组件各层温度的计算方法 | |
Abdolzadeh et al. | Optical and thermal modeling of a photovoltaic module and experimental evaluation of the modeling performance | |
Baig et al. | Modelling and experimental analysis of a seasonally tracked V-trough PV/T system in India | |
CN106227981B (zh) | 一种双玻光伏组件热电耦合分析方法 | |
Dey et al. | Numerical modeling and performance assessment of elongated compound parabolic concentrator based LCPVT system | |
Wu et al. | Coupled optical-electrical-thermal analysis of a semi-transparent photovoltaic glazing façade under building shadow |
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 |