CN115130395A - 一种水上光伏对下覆封闭水域水温影响的计算方法 - Google Patents
一种水上光伏对下覆封闭水域水温影响的计算方法 Download PDFInfo
- Publication number
- CN115130395A CN115130395A CN202210606357.9A CN202210606357A CN115130395A CN 115130395 A CN115130395 A CN 115130395A CN 202210606357 A CN202210606357 A CN 202210606357A CN 115130395 A CN115130395 A CN 115130395A
- Authority
- CN
- China
- Prior art keywords
- calculating
- photovoltaic
- water
- substep
- radiation
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 176
- 238000000034 method Methods 0.000 title claims abstract description 28
- 230000005855 radiation Effects 0.000 claims abstract description 122
- 238000004364 calculation method Methods 0.000 claims abstract description 39
- 238000001704 evaporation Methods 0.000 claims abstract description 37
- 230000008020 evaporation Effects 0.000 claims abstract description 37
- 230000004907 flux Effects 0.000 claims abstract description 35
- 230000008859 change Effects 0.000 claims abstract description 23
- 230000010354 integration Effects 0.000 claims description 15
- 229920006395 saturated elastomer Polymers 0.000 claims description 13
- 238000003491 array Methods 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 claims description 8
- 238000013461 design Methods 0.000 claims description 6
- 230000009467 reduction Effects 0.000 claims description 4
- 238000009834 vaporization Methods 0.000 claims description 4
- 230000008016 vaporization Effects 0.000 claims description 4
- 239000002352 surface water Substances 0.000 claims description 3
- 238000011161 development Methods 0.000 abstract description 2
- 238000005338 heat storage Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 241000251468 Actinopterygii Species 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000005286 illumination Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000017525 heat dissipation Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000008239 natural water Substances 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000012950 reanalysis Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 230000005476 size effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Fluid Mechanics (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Marketing (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
本发明涉及一种水上光伏对下覆封闭水域水温影响的计算方法,包括:收集信息;计算太阳短波辐射;计算光伏遮挡下水面接受的大气长辐射;计算蒸发散射;计算对流热通量;计算水体向外长波辐射热通量;计算t 0时刻自由表面总热通量;计算受光伏影响的t 0→t 1时刻水体水温变化;计算t 1时刻受光伏影响的水温。本发明通过水体热量平衡计算方法和各通量分项影响因子受光伏阵列影响的计算公式,针对该封闭水域内的桩基水上光伏电站,依据当地的气象条件和即将建设的光伏阵列布设参数,定量计算水上光伏电站对水温变化的影响,根据此方法计算的水温及水温变化可作为评估光伏电站对水环境影响的依据,为将来光伏电站建成后,对该水域的开发利用提供了科学依据。
Description
技术领域
本发明涉及一种水上光伏对下覆封闭水域水温影响的计算方法,是一种水文观测计算方法,是一种用于水面光伏电站的光伏阵列影响的计算方法。
背景技术
水上光伏电站在水面上布设大量密集的光伏板阵列,光伏板阵列通过遮蔽和吸收太阳短波辐射、改变局部风场等方式,影响水面接受的辐射热和水表蒸发、对流等换热过程,进而改变水体温度。作为水生态的重要环境因子,水温是评估水上光伏电站对水环境影响的重要指标。目前主要基于有限的观测进行光伏电站对水温影响规律的定性分析,缺少系统性的定量计算和评价方法。
发明内容
为了克服现有技术的问题,本发明提出了一种水上光伏对下覆封闭水域水温影响的计算方法。所述的方法针封闭水域内的桩基水上光伏电站,提出了依据气象条件和光伏阵列布设参数,定量计算水上光伏电站中水温变化的方法,根据此方法计算的水温及水温变化可作为评估光伏电站对水环境影响的依据。
本发明的目的是这样实现的:一种水上光伏对下覆封闭水域水温影响的计算方法,所述的方法的步骤如下:
步骤1,收集信息:收集当地水文气象资料,以及预建造的光伏阵列的各项设计数据,光伏阵列设计数据包括:光伏板尺寸、光伏板间距、支柱分布间距;当地水文气象资料包括:太阳辐射总辐射量I,通过计算或收集获取散射辐射量Id;大气长波辐射、气温、云量、湿度、气压、水面风速、水温、水深;
步骤2,计算太阳短波辐射:计算光伏遮挡下水面接受的太阳短波辐射,太阳短波辐射的计算包括如下子步骤:
子步骤2.1,根据收集的信息整理获得太阳辐射总辐射量I,通过计算或收集获取散射辐射量Id;
子步骤2.2,计算天空视域比例因子CFsky:
式中θ1、θ2、θ3、θ4分别为计算空隙后板和前板下檐、上檐到积分点的夹角,n为积分区间L的分段数;L为前后阵列之间宽度;mlmax为积分区间左侧最大空隙数;mrmax为积分区间右侧含区分区间本身的最大空隙数;
子步骤2.3,光伏板后檐对应阴影宽度Dh:
Dh=Hpv·cotα·cosγ
式中α为太阳高度角,γ为太阳方位角,Hpv为光伏板上下沿高度差;
子步骤2.4,计算透过比例因子:直射辐射在前后阵列之间宽度为L内的透过比例因子CFdir为:
Dpv为光伏板前檐到后檐在水平面投影面中的距离;
子步骤2.5,计算光伏区下垫面总辐射量Ipv:
Ipv=(I-Id)CFdir+IdCFsky;
步骤3,计算光伏遮挡下水面接受的大气长辐射:
子步骤3.1,根据收集的信息获得实测大气长波辐射或通过实测的气温、云量计算获得无光伏影响的大气长波辐射值;
大气长波辐射计算公式:
Ea=εσTa0 4g(Fc)
式中:ε为发射率;σ为玻尔兹曼常数;Ta0为t0时刻气温;g(Fc)为云量Fc的函数:
g(Fc)=(1+0.17Fc2);
子步骤3.2,根据光伏对天空遮挡的视域因子,收集的云量信息求解光伏区考虑光伏板遮挡后的等效云量系数Fc′:
Fc′=(1-CFsky)+CFsky·Fc;
子步骤3.3,计算有光伏影响下的等效云量函数:
g(Fc′)=(1+0.17Fc′2);
子步骤3.4,计算有光伏影响下的大气长波辐射Eapv:
Eapv=Ea·g(F′C)/g(Fc);
步骤4,计算蒸发散射:
子步骤4.1,根据收集信息资料获得光伏区风速,采用已有光伏实测水面风速u为风速(m/s),或依据类似已建光伏区实测风速折减比rs和拟计算的光伏区附近无光伏影响气象站风速uc计算:
u=rs·uc
其中:rs=u0/uc0,
u0:类似已建光伏区地面或水面风速;uc0:类似已建光伏区外围附近气象站地面或水面风速;
子步骤4.2,计算水面饱和水汽压力、空气中水汽压力:
t0时刻水面饱和水汽压力es:
es=f(Ts0)
其中:ea为t0时刻空气中水汽压力:
ea=f(Ta0)
其中:Ts0为t0时刻水表面温度,Ta0为t0时刻气温;
水汽压力f(T)计算公式:
其中:T为水表面温度或气温;E0为纯水平液面饱和水汽压,计算公式:
其中:Tc为水的三相点温度;T为水表面温度或气温;
子步骤4.3,计算蒸发量:顿蒸发公式形式:
E=f(u)·(es-ea)
其中:es为水面饱和水汽压力;ea为空气中水汽压力;
风速函数f(u)计算:
子步骤4.4,计算蒸发潜热:蒸发潜热Lv(J/kg)计算公式:
Lv=2.5×106-2.3×103·Ts0
其中:Ts0为t0时刻表层水温;
子步骤4.5:蒸发潜热通量QE可通过单位时间蒸发量得到:
QE=LvρE
其中:E为水面单位时间蒸发量;ρ为水体密度(kg/m3);Lv为蒸发潜热;
步骤5,计算对流热通量:
子步骤5.1:计算波文比b:
b=0.62·Pa/1000(Ts0-Ta0)/(es-ea);
子步骤5.2:根据波文比计算对流热通量QH:
QH=b·QE
其中:QE为蒸发热通量;
步骤6,计算水体向外长波辐射热通量EW:水体向外长波辐射热通量计算公式:
Ew=εσTs0 4
步骤7,计算t0时刻自由表面总热通量F:
F=Ipv(1-αs)+Eapv-Ew-QE-QH
其中,αs为太阳短波辐射反射比;
步骤8,计算受光伏影响的t0→t1时刻水体水温变化ΔTs:
ΔTs=F/(ρ·c·H)
其中:c为水的比热,H为水深;
步骤9,计算t1时刻受光伏影响的水温Ts1:
Ts1=Ts0+ΔTs。
进一步的,子步骤2.1中所述的环境散射辐射值的计算方式包括如下子步骤:
子步骤2.11,计算大气层外辐射量I0:
其中,nd为一年中的日序号;
子步骤2.12,计算晴空指数kT:
子步骤2.13,计算散射辐射与总辐射之比Id/I:下垫面散射辐射Id与总辐射I之比的计算公式:
子步骤2.14,计算散射辐射:
Id=I·rd。
本发明的优点和有益效果是:本发明通过水体热量平衡计算方法,以及一系列各通量分项影响因子受光伏阵列影响的计算公式,针对该封闭水域内的桩基水上光伏电站,依据当地的气象条件和即将建设的光伏阵列布设参数,定量计算水上光伏电站对水温变化的影响,根据此方法计算的水温及水温变化可作为评估光伏电站对水环境影响的依据,为将来光伏电站建成后,对该水域的开发利用提供了科学依据。
附图说明
下面结合附图和实施例对本发明作进一步说明。
图1是本发明实施例一所述方法多块光伏板公式的参数标注示意图;
图2是本发明实施例一所述方法的流程图;
图3是本发明实施例二所述方法的应用实例计算赣榆鱼塘光伏区水温与实测值对比。
具体实施方式
实施例一:
本实施例是一种水上光伏对下覆封闭水域水温影响的计算方法。所述方法通过水体热量平衡计算,以及一系列各通量分项影响因子受光伏阵列影响的计算公式,针对该封闭水域内的桩基水上光伏电站,依据当地的气象条件和即将建设的光伏阵列布设参数,定量计算水上光伏电站对水温变化的影响。
(一)水体热量平衡计算:
水体温度变化是水体蓄热量变化的直接反映。自然水体的蓄热量变化过程主要来自太阳辐射以及与大气之间的热交换,具体包括水面太阳短波辐射、水体和大气间的长波辐射、蒸发(潜热)和对流(显热)通量,此外底床与水体之间的热交换对水温也对水体能量变化略有影响。对于与外界无流动交换的封闭水体单元,根据能量守恒定律,水体热量平衡采用下式进行计算:
QG=E*+QH+QE+Qb (1)
式中:E*为水体吸收的净辐射(辐射差额);QH为水-气间的对流(显热)通量,即为水表与低层大气之间的对流热交换;QE为潜热通量,即水面蒸发热交换;QG为水体蓄热变化量,Qb为底床与水体之间的热交换,底部热交换一般认为是小量,可以忽略。
水体吸收的净辐射E*表示为:
E*=Eg↓+El↓-Er↑-El↑ (2)
式中Eg↓为太阳向下短波辐射,Er↑为水面短波反射辐射,是总辐射到达水面后,经水面向上反射的那部分短波辐射(可用垂向向下布设的总辐射表测量)。El↓为向下长波辐射,El↑为向上长波辐射,上述四组分量数据可通过四分量辐射表测量获取。
水体蓄热变化量QG采用实测分层水温计算得到:
QG=[(Σρc·ΔH.T)t2-(Σρc·ΔH·T)t1]/(t2-t1) (3)
其中ρ为水的密度(kg/m3),c为水的比热(J/kg℃),ΔH为水层厚度(m),T为对应层水温(℃),即利用实测水温进行水体能量垂向积分,得到对应时刻的水体内能,t2与t1时刻的水体内能差即为水体蓄热量变化值。
蒸发潜热通量QE可通过单位时间蒸发量得到:
QE=LvρE (4)
其中E为水面单位时间蒸发量(m/s),Lv为蒸发潜热(J/kg),表达式为:
Lv=2.5·106-2.3·103Ts (5)
式中Ts为表皮水温(℃)。
蒸发量E可以通过实际测量获得,但由于受蒸发皿尺寸效应等影响实测资料准确性难以确定,需要通过其它计算方法进行检验修正。本实施例选用经典的道尔顿蒸发公式:
E=f(u)·(es-ea) (6)
式中f(u)为风速函数,u为风速(m/s),es为水面饱和水汽压力(hPa),ea为空气中水汽压力(hPa)。
风速函数可以采用闵骞修正的风速函数,公式为
对流(显热)通量QH可通过波文比方法得到:
QH=0.62·Pa/1000·(Ts-Ta)/(es-ea) (8)
式中Pa为大气压力(hPa),Ta为气温(℃)。
式(1)中所述各组热通量项受光伏遮挡均会发生变化,其中辐射通量E*的遮挡与光伏板阵列的空间布置有关,短波辐射遮挡还会与太阳角度有关。蒸发和对流散热计算跟气温、湿度和风速相关,其中受光伏板影响较为显著的是风速。因此水上光伏对水温变化影响的计算分析,应主要考虑光伏阵列对太阳短波辐射、大气长波辐射的遮挡效应,及对风速的改变效应。
(二)各通量分项影响因子受光伏阵列影响的计算公式:
i.太阳短波辐射:
根据太阳直射辐射和散射辐射资料计算短波辐射动态遮蔽系数。
散射辐射为各向同性漫射,折减系数采用天空视域因子,对多排缝隙判断该缝隙的空域是否入射到计算单元以及计算入射角度的方法如下:
特定缝隙的散射光能否进入积分单元要判断缝隙左侧板的最低入射角是否大于右侧板的最高入射角,如果是,则该二夹角间的空域为积分单元上可视空域,光线可达,否则该区域空隙被前后板互相遮挡,无法到达积分单元,总的天空视域因子CFsky计算公式如下:
式中θ1、θ2、θ3、θ4分别为计算空隙后板和前板下檐、上檐到积分点的夹角,见图1所示,n为积分区间L的分段数,mlmax为积分区间左侧最大空隙数,mrmax为积分区间右侧(含区分区间本身)最大空隙数。图1显示了光伏阵列某个特定单元的垂向立面图,简略的显示了四排光伏板,并将中间的两个光伏板按照由南向北的方向设定为前排光伏板1和后排光伏板2,并将阳光入射方向3设定在入射角α的角度上,将积分点4设在前排光伏板上檐在水面5的投影点(线)上。
直射辐射入射比例可根据太阳高度角和太阳方位角求解光照区比例获得。对于图1显示的特定单元,当单元位于阴影区6时,直射辐射透过比为0,当单元位于光照区7时,直射辐射透过比为1。对于整个光伏区,前后两阵列之间光照区宽度占间距的比值即为直射透过比例。
光伏板后檐对应阴影宽度Dh,如图1所示:
Dh=Hpv·cotα·cosβ (10)
式中:Dh为光伏板上檐和下檐之间的高度差;α为太阳高度角;β为太阳方位角;α、β与时间有关,可采用相关公式求解。
直射辐射在前后阵列宽度L内的透过比例为:
根据实测的无光伏水域(简称“环境区”)短波总辐射I(总辐射包括直射辐射和散射辐射)和散射辐射值Id,按照式(9)和(11)计算方法分别计算光伏区散射辐射和直射辐射的动态遮蔽系数,并将二者折减后的值相加得到光伏区总辐射Ipv;
Ipv=(I-Id)CFdir+IdCFsky (12)
ii.长波辐射:
光伏板阵列遮挡会减少大气向水体的长波辐射,大气长波辐射公式的一般形式为:
Ea=εσTa 4G(Fc) (13)
式中ε为发射率,可通过实测气象资料进行率定,σ为玻尔兹曼常数,取值5.67×10-8(W·m-2·℃-4);g(Fc)为云量Fc的函数,常用的形式为:
g(Fc)=(1+0.17Fc2) (14)
光伏对天空遮挡的比例系数Fpv可以采用太阳短波散射辐射计算的天空视域因子CFsky:
Fpv=1-CFsky (15)
将光伏遮挡区视为云量,光伏区考虑光伏板遮挡后的等效云量系数Fc′为:
Fc′=Fpv+CFsky·Fc (16)
将Fc′带入(14)式,再按(13)式求解得到光伏区受遮挡的大气辐射。
iii.风速:
光伏板对风的影响可用风变系数a表征:
风变系数a与光伏板布置间距、倾角、高度,以及风向等均相关,可通过原型观测统计分析,或者模型试验或数值计算方法获取,已有文献表明,风变系数a取值大约在0.2~1.2之间。之后带入(6)式求解光伏遮挡区的蒸发和对流散热。
在获取上述关键系数后,可确定公式(1)和(2)中各分项的准确计算值,然后带入公式(3)。对(3)式进行向前时间积分,便可获得光伏区水温T的实时变化过程。
本实施例所述方法的具体步骤如下,流程如图2所示:
步骤1,收集信息:收集当地水文气象资料,以及预建造的光伏阵列的各项设计数据,光伏阵列设计数据包括:光伏板尺寸、光伏板间距、支柱分布间距;当地水文气象资料包括:太阳辐射总辐射量I,通过计算或收集获取散射辐射量Id;大气长波辐射、气温、云量、湿度、气压、水面风速、水温、水深。
本实施例是在将要建设光伏电站的封闭区域,进行水文数据收集,包括:太阳辐射量等。同时收集拟建设的光伏电站的各项参数,如光伏板的大小尺寸、间距、支撑柱的高度、间距等。
步骤2,计算太阳短波辐射:计算光伏遮挡下水面接受的太阳短波辐射,太阳短波辐射的计算包括如下子步骤:
子步骤2.1,根据收集的信息整理获得太阳辐射总辐射量I,通过计算或收集获取散射辐射量Id。
太阳辐射总辐射量在一些水文资料比较齐全的地方,可以通过查找获得,但在多数情况下,因此可以通过检测当地的太阳短辐射以及长辐射量计算得到。
子步骤2.2,计算天空视域比例因子CFsky:
式中θ1、θ2、θ3、θ4(见图1)分别为计算空隙后板和前板下檐、上檐到积分点的夹角,n为积分区间L的分段数;L为前后阵列之间宽度,即阵列前后两排光伏板之间从前排前檐到后排前檐;mlmax为积分区间左侧最大空隙数;mrmax为积分区间右侧含区分区间本身的最大空隙数。
子步骤2.3,光伏板后檐对应阴影宽度Dh:
Dh=Hpv·cotα·cosγ
式中α为太阳高度角,γ为太阳方位角,Hpv为光伏板上下沿高度差;
子步骤2.4,计算透过比例因子:直射辐射在前后阵列之间宽度为L内的透过比例因子CFdir为:
其中:Dpv为光伏板前檐到后檐在水平面投影面中的距离;
子步骤2.5,计算光伏区下垫面总辐射量Ipv:
Ipv=(I-Id)CFdir+IdCFsky;
步骤3,计算光伏遮挡下水面接受的大气长辐射:
子步骤3.1,根据收集的信息获得实测大气长波辐射或通过实测的气温、云量计算获得无光伏影响的大气长波辐射值;
大气长波辐射计算公式:
Ea=εσTa0 4g(Fc)
式中:ε为发射率,可通过实测气象资料进行率定,默认值取0.97;σ为玻尔兹曼常数,取值5.67×10-8(W·m-2·℃-4);Ta0为t0时刻气温(K);g(Fc)为云量Fc的函数:
g(Fc)=(1+0.17Fc2);
子步骤3.2,根据光伏对天空遮挡的视域因子,收集的云量信息求解光伏区考虑光伏板遮挡后的等效云量系数Fc′:
Fc′=(1-CFsky)+CFsky·Fc;
子步骤3.3,计算有光伏影响下的等效云量函数:
g(Fc′)=(1+0.17Fc′2);
子步骤3.4,计算有光伏影响下的大气长波辐射Eapv:
Eapv=Ea·g(F′C)/g(Fc);
步骤4,计算蒸发散射:
子步骤4.1,根据收集信息资料获得光伏区风速,采用已有光伏实测水面风速u为风速(m/s),或依据类似已建光伏区实测风速折减比rs和拟计算的光伏区附近无光伏影响气象站风速uc计算:
u=rs·uc
其中:rs=u0/uc0,
u0:类似已建光伏区地面或水面风速;uc0:类似已建光伏区外围附近气象站地面或水面风速;
子步骤4.2,计算水面饱和水汽压力、空气中水汽压力:
t0时刻水面饱和水汽压力es:
es=f(Ts0)
其中:ea为t0时刻空气中水汽压力:
ea=f(Ta0)
其中:Ts0为t0时刻水表面温度(K),Ta0为t0时刻气温(K);t0时刻是指采样开始的时刻。如果研究白昼的影响,可以将采样起始时间设定为早晨太阳升起前,例如早上5点钟,可以设定时间步长为半小时,计算结束时间为下午太阳落山时,例如下午6点。如果研究昼夜或多日的影响过程,采样起始时间可以设定为拟研究的日分界点,无需特殊限定至某一时刻。
水汽压力f(T)计算公式:
其中:T为水表面温度或气温(K);E0为纯水平液面饱和水汽压(hPa),计算公式:
其中:Tc为水的三相点温度,Tc=273.16K;T为水表面温度或气温(K);
子步骤4.3,计算蒸发量:顿蒸发公式形式:
E=f(u)·(es-ea)
其中:es为水面饱和水汽压力(hPa);ea为空气中水汽压力(hPa);
风速函数f(u)计算:
子步骤4.4,计算蒸发潜热:蒸发潜热Lv(J/kg)计算公式:
Lv=2.5×106-2.3×103·Ts0
其中:Ts0为t0时刻表层水温(℃);
子步骤4.5:蒸发潜热通量QE可通过单位时间蒸发量得到:
QE=LvρE
其中:E为水面单位时间蒸发量(m/s);ρ为水体密度(kg/m3);Lv为蒸发潜热(J/kg);
步骤5,计算对流热通量:
子步骤5.1:计算波文比b:
b=0.62·Pa/1000(Ts0-Ta0)/(es-ea);
子步骤5.2:根据波文比计算对流热通量QH:
QH=b·QE
其中:QE为蒸发热通量;
步骤6,计算水体向外长波辐射热通量EW:水体向外长波辐射热通量计算公式:
Ew=εσTs0 4
步骤7,计算t0时刻自由表面总热通量F:
F=Ipv(1-αs)+Eapv-Ew-QE-QH
其中,αs为太阳短波辐射反射比;可通过实测资料率定,默认取值0.06。
步骤8,计算受光伏影响的t0→t1时刻水体水温变化ΔTs:
ΔTs=F/(ρ·c·H)
其中:c为水的比热(J/kg℃),H为水深(m);
步骤9,计算t1时刻受光伏影响的水温Ts1:
Ts1=Ts0+ΔTs。
实施例二:
本实施例是实施例一的改进,是实施例一关于环境辐射值计算方法的细化。子步骤2.1中所述的环境散射辐射值的计算方式包括如下子步骤:
子步骤2.11,计算大气层外辐射量I0:
其中,n为一年中的日序号;
子步骤2.12,计算晴空指数kT:
子步骤2.13,计算散射辐射与总辐射之比Id/I:下垫面散射辐射Id与总辐射I之比的计算公式:
子步骤2.14,计算散射辐射:
Id=I·rd。
应用实例:
应用实例位于江苏省连云港市赣榆区鱼塘上,鱼塘属于封闭静止水域。根据实际光伏板阵列布设,采用对照组鱼塘(未布设光伏阵列的环境区)观测的短波总辐射、短波散射辐射、长波辐射、气压、风速和光伏区鱼塘的风速,通过气象再分析资料CFSv2数据获取的云量系数,利用本实施例所述方法进行光伏区鱼塘水温的计算。计算的水温值与实测值平均在0.5℃以内,如图3所示,计算的水温值与实测值之间的偏差始终保持在特定区间,未随时间出现累积性偏差,说明本实施例的计算公式及参数是合理的,可以用于快速准确的分析封闭或近似封闭水域中光伏电站下垫面水温的影响。
最后应说明的是,以上仅用以说明本发明的技术方案而非限制,尽管参照较佳布置方案对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案(比如光伏阵列的排布方式、各种公式的运用、步骤的先后顺序等)进行修改或者等同替换,而不脱离本发明技术方案的精神和范围。
Claims (2)
1.一种水上光伏对下覆封闭水域水温影响的计算方法,其特征在于,所述的方法的步骤如下:
步骤1,收集信息:收集当地水文气象资料,以及预建造的光伏阵列的各项设计数据,光伏阵列设计数据包括:光伏板尺寸、光伏板间距、支柱分布间距;当地水文气象资料包括:太阳辐射总辐射量I,通过计算或收集获取散射辐射量Id;大气长波辐射、气温、云量、湿度、气压、水面风速、水温、水深;
步骤2,计算太阳短波辐射:计算光伏遮挡下水面接受的太阳短波辐射,太阳短波辐射的计算包括如下子步骤:
子步骤2.1,根据收集的信息整理获得太阳辐射总辐射量I,通过计算或收集获取散射辐射量Id;
子步骤2.2,计算天空视域比例因子CFsky:
式中θ1、θ2、θ3、θ4分别为计算空隙后板和前板下檐、上檐到积分点的夹角,n为积分区间L的分段数;L为前后阵列之间宽度;mlmax为积分区间左侧最大空隙数;mrmax为积分区间右侧含区分区间本身的最大空隙数;
子步骤2.3,光伏板后檐对应阴影宽度Dh:
Dh=Hpv·cotα·cosγ
式中α为太阳高度角,γ为太阳方位角,Hpv为光伏板上下沿高度差;
子步骤2.4,计算透过比例因子:直射辐射在前后阵列之间宽度为L内的透过比例因子CFdir为:
其中:Dpv为光伏板前檐到后檐在水平面投影面中的距离;
子步骤2.5,计算光伏区下垫面总辐射量Ipv:
Ipv=(I-Id)CFdir+IdCFsky;
步骤3,计算光伏遮挡下水面接受的大气长辐射:
子步骤3.1,根据收集的信息获得实测大气长波辐射或通过实测的气温、云量计算获得无光伏影响的大气长波辐射值;
大气长波辐射计算公式:
Ea=εσTa0 4g(Fc)
式中:ε为发射率;σ为玻尔兹曼常数;Ta0为t0时刻气温;g(Fc)为云量Fc的函数:
g(Fc)=(1+0.17Fc2);
子步骤3.2,根据光伏对天空遮挡的视域因子,收集的云量信息求解光伏区考虑光伏板遮挡后的等效云量系数Fc′:
Fc′=(1-CFsky)+CFsky·Fc;
子步骤3.3,计算有光伏影响下的等效云量函数:
g(Fc′)=(1+0.17Fc′2);
子步骤3.4,计算有光伏影响下的大气长波辐射Eapv:
Eapv=Ea·g(F′C)/g(Fc);
步骤4,计算蒸发散射:
子步骤4.1,根据收集信息资料获得光伏区风速,采用已有光伏实测水面风速u为风速(m/s),或依据类似已建光伏区实测风速折减比rs和拟计算的光伏区附近无光伏影响气象站风速uc计算:
u=rs·uc
其中:rs=u0/uc0,
u0:类似已建光伏区地面或水面风速;uc0:类似已建光伏区外围附近气象站地面或水面风速;
子步骤4.2,计算水面饱和水汽压力、空气中水汽压力:
t0时刻水面饱和水汽压力es:
es=f(Ts0)
其中:ea为t0时刻空气中水汽压力:
ea=f(Ta0)
其中:Ts0为t0时刻水表面温度,Ta0为t0时刻气温;
水汽压力f(T)计算公式:
其中:T为水表面温度或气温;E0为纯水平液面饱和水汽压,计算公式:
其中:Tc为水的三相点温度;T为水表面温度或气温;
子步骤4.3,计算蒸发量:顿蒸发公式形式:
E=f(u)·(es-ea)
其中:es为水面饱和水汽压力;ea为空气中水汽压力;
风速函数f(u)计算:
子步骤4.4,计算蒸发潜热:蒸发潜热Lv(J/kg)计算公式:
Lv=2.5×106-2.3×103·Ts0
其中:Ts0为t0时刻表层水温;
子步骤4.5:蒸发潜热通量QE可通过单位时间蒸发量得到:
QE=LvρE
其中:E为水面单位时间蒸发量;ρ为水体密度(kg/m3);Lv为蒸发潜热;
步骤5,计算对流热通量:
子步骤5.1:计算波文比b:
b=0.62·Pa/1000(Ts0-Ta0)/(es-ea);
子步骤5.2:根据波文比计算对流热通量QH:
QH=b·QE
其中:QE为蒸发热通量;
步骤6,计算水体向外长波辐射热通量EW:水体向外长波辐射热通量计算公式:
Ew=εσTs0 4
步骤7,计算t0时刻自由表面总热通量F:
F=Ipv(1-αs)+Eapv-Ew-QE-QH
其中,αs为太阳短波辐射反射比;
步骤8,计算受光伏影响的t0→t1时刻水体水温变化ΔTs:
ΔTs=F/(ρ·c·H)
其中:c为水的比热,H为水深;
步骤9,计算t1时刻受光伏影响的水温Ts1:
Ts1=Ts0+ΔTs。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210606357.9A CN115130395B (zh) | 2022-05-31 | 2022-05-31 | 一种水上光伏对下覆封闭水域水温影响的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210606357.9A CN115130395B (zh) | 2022-05-31 | 2022-05-31 | 一种水上光伏对下覆封闭水域水温影响的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115130395A true CN115130395A (zh) | 2022-09-30 |
CN115130395B CN115130395B (zh) | 2023-04-25 |
Family
ID=83378202
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210606357.9A Active CN115130395B (zh) | 2022-05-31 | 2022-05-31 | 一种水上光伏对下覆封闭水域水温影响的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115130395B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103985060A (zh) * | 2014-04-25 | 2014-08-13 | 云南省电力设计院 | 一种基于测光数据的光伏电站发电量计算方法 |
CN104406686A (zh) * | 2014-12-10 | 2015-03-11 | 西北师范大学 | 复杂地形条件下太阳短波入射辐射估算方法 |
CN112541275A (zh) * | 2020-12-21 | 2021-03-23 | 中国水利水电科学研究院 | 一种冰期河湖与大气热交换线性化模型构建方法 |
CN113987823A (zh) * | 2021-11-08 | 2022-01-28 | 国家气候中心 | 陆上集中式光伏电站气候效应评估的中尺度数值模拟方法 |
-
2022
- 2022-05-31 CN CN202210606357.9A patent/CN115130395B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103985060A (zh) * | 2014-04-25 | 2014-08-13 | 云南省电力设计院 | 一种基于测光数据的光伏电站发电量计算方法 |
CN104406686A (zh) * | 2014-12-10 | 2015-03-11 | 西北师范大学 | 复杂地形条件下太阳短波入射辐射估算方法 |
CN112541275A (zh) * | 2020-12-21 | 2021-03-23 | 中国水利水电科学研究院 | 一种冰期河湖与大气热交换线性化模型构建方法 |
CN113987823A (zh) * | 2021-11-08 | 2022-01-28 | 国家气候中心 | 陆上集中式光伏电站气候效应评估的中尺度数值模拟方法 |
Non-Patent Citations (1)
Title |
---|
石兴鹏等: "干旱区水上太阳能电池板与浮球联合作用对节水效率的影响", 《水电能源科学》 * |
Also Published As
Publication number | Publication date |
---|---|
CN115130395B (zh) | 2023-04-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | A mathematical model of global solar radiation to select the optimal shape and orientation of the greenhouses in southern China | |
Chen et al. | Effects of foliage clumping on the estimation of global terrestrial gross primary productivity | |
Mobtaker et al. | Solar energy conservation in greenhouse: Thermal analysis and experimental validation | |
CN106909985B (zh) | 一种光伏发电预测系统的预测方法 | |
Shi et al. | China land soil moisture EnKF data assimilation based on satellite remote sensing data | |
Jim et al. | Estimating heat flux transmission of vertical greenery ecosystem | |
Huang et al. | An overview of the semi-arid climate and environment research observatory over the Loess Plateau | |
Almorox et al. | Total solar irradiance's effect on the performance of empirical models for estimating global solar radiation: An empirical-based review | |
Zhang et al. | Dynamic energy balance model of a glass greenhouse: An experimental validation and solar energy analysis | |
Jin et al. | New observational evidence for global warming from satellite | |
Heusinger et al. | Introduction, evaluation and application of an energy balance model for photovoltaic modules | |
Sun et al. | The response of surface mass and energy balance of a continental glacier to climate variability, western Qilian Mountains, China | |
Tan et al. | Impact of different roofing mitigation strategies on near-surface temperature and energy consumption over the Chicago metropolitan area during a heatwave event | |
Leng et al. | A method for deriving all‐sky evapotranspiration from the synergistic use of remotely sensed images and meteorological data | |
Jiang et al. | Observed impacts of utility-scale photovoltaic plant on local air temperature and energy partitioning in the barren areas | |
Zou et al. | Comparison of two satellite‐based evapotranspiration models of the Nagqu River Basin of the Tibetan Plateau | |
Gelsor et al. | Solar energy on the Tibetan Plateau: Atmospheric influences | |
Chen et al. | Performance analysis of radiation and electricity yield in a photovoltaic panel integrated greenhouse using the radiation and thermal models | |
Wang et al. | Urban thermal environment and surface energy balance in 3D high-rise compact urban models: Scaled outdoor experiments | |
Bhatt et al. | A climatological‐dynamical analysis associated with precipitation around the southern part of the Himalayas | |
Schroeder et al. | Validation of solar radiation surfaces from MODIS and reanalysis data over topographically complex terrain | |
Wu et al. | A systematic analysis of multiple structural parameters of Chinese solar greenhouse based on the thermal performance | |
Bo et al. | Light environment simulation for a three-span plastic greenhouse based on greenhouse light environment simulation software | |
CN115049520B (zh) | 一种计算光伏阵列对短波辐射动态遮蔽系数的方法 | |
CN115130395B (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 |