计算某地的从日出到日落各时刻的水平面直射和散射辐射，倾斜面总辐射的Python程序

已于 2022-10-16 20:32:56 修改
阅读量1.9k 收藏 12
点赞数 1
分类专栏： 数值传热 文章标签： python
于 2021-01-16 08:41:51 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_48615832/article/details/112130737
版权

基本概念

在这里插入图片描述

大气质量

在这里插入图片描述
在这里插入图片描述
大气质量m只是

大气透明度τ

直射辐射透明度τb，散射辐射透明度τd。
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
晴空垂直直射辐射Gcnb
在这里插入图片描述
Gon是一年中第n天的地球外太阳辐射（a式简单，b式更精确）(大气层外法向的辐照度)

from sympy import *
month = 8
day = 22
#n()是一年里的第几天
def n():
    d_num = 0
    for m in range(1, month):
        if m in [1, 3, 5, 7, 8, 10, 12]:
            d_num += 31
        elif m in [2]:
            d_num += 28
        else:
            d_num += 30
    d_num += day
    return d_num
#日子数
n = n()
#计算大气层外法向的辐照度
G_sc = 1367  #单位W/m2
G_on1 = (G_sc * (1+0.033*cos(360*n/365*pi/180))).evalf()
B = (n-1)*360/365
G_on2 = (G_sc * (1.000110+0.034221*cos(B*pi/180)+0.001280*sin(B*pi/180)+0.000719*cos(2*B*pi/180)+0.000077*sin(2*B*pi/180))).evalf()
print(G_on1)
print(G_on2)

结果：

1338.48518301793
1335.56388385576

晴空水平直射辐射Gcb
对于一小时时间段内，晴空水平直射辐射Icb
在这里插入图片描述
Tropical —— 热带，0°—30°南北纬
Midlatitude summer —— 中纬度夏季，30°—60°南北纬
Subarctic summer —— 高纬度夏季，60°—90°南北纬
Midlatitude winter —— 中纬度冬季，30°—60°南北纬

晴空，垂直于辐射方向上的太阳辐射照度Gcnb

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

水平面上的直射辐照度Gcb

在这里插入图片描述

一小时内，水平面上直射辐照量Icb

在这里插入图片描述

水平面上的散射辐照度Gcd

在这里插入图片描述

一小时内，水平面上散射辐照量Icd

在这里插入图片描述

一小时内，水平面上的总辐照量Ic

在这里插入图片描述
在这里插入图片描述

大气层外水平面上的太阳辐射照度G0，W/m2

在这里插入图片描述

G_sc = 1367  #单位W/m2
G0 = G_sc * (1+0.033*cos(360*n/365*pi/180)) * cos(Theta_z*pi/180)

大气层外水平面上，1天内太阳的辐照量H0，J/m2

在这里插入图片描述

from sympy import *
G_sc = 1367  #单位W/m2
H0 = 24*3600*G_sc/pi * (1+0.033*cos(360*n/365*pi/180)) * (2*pi*Omega_set/360*sin(Delta*pi/180)*sin(Phi*pi/180) + cos(Delta*pi/180)*cos(Phi*pi/180)*cos(Omega_set*pi/180))

大气层外水平面上，每小时内太阳的辐照量I0，J/m2

在这里插入图片描述
在这里插入图片描述

from sympy import *
G_sc = 1367  #单位W/m2
I0 = 12*3600*G_sc/pi * (1+0.033*cos(360*n/365*pi/180)) * (2*pi*(Omega_2-Omega_1)/360*sin(Delta*pi/180)*sin(Phi*pi/180) + cos(Delta*pi/180)*cos(Phi*pi/180)*(sin(Omega_2*pi/180)-sin(Omega_1*pi/180)))

晴空指数KT

在这里插入图片描述

小时散射辐射量Id与小时总的辐射量I之比，与kT的关系

在这里插入图片描述

from sympy import *
I = 1    #小时的日辐照量
I0 = 2    #同小时大气层外日辐照量
kT = I/I0    #小时晴空指数
Id_I = None    #Id/I
if kT <= 0.22:
    Id_I  = 1 - 0.09*kT
elif kT <= 0.8:
    Id_I = 0.9511 - 0.1604*kT + 4.388*kT**2 - 16.638*kT**3 + 12.336*kT**4
else:
    Id_I = 0.165

或者可以采用下面的式子，这两个式子的结果都差不多，但前者中心部分精确，后者计算简单。
在这里插入图片描述

from sympy import *
I = 1    #小时的日辐照量
I0 = 2    #同小时大气层外日辐照量
kT = I/I0    #日晴空指数
Id_I = None    #Id/I
if kT <= 0.35:
    Id_I  = 1 - 0.249*kT
elif kT <= 0.75:
    Id_I = 1.557 - 1.84*kT 
else:
    Id_I = 0.177

在这里插入图片描述

小时散射辐射量Id与小时总的辐射量I之比，与k_Tc的关系

在这里插入图片描述

I = 1    #小时的日辐照量
I_c = 2    #同小时标准晴空条件下大气层外日辐照量
kTc = I/I_c    #日标准晴空指数
Id_I = None    #Id/I
if kTc < 0.48:
    Id_I  = 1 - 0.1*kTc
elif kTc < 1.10:
    Id_I = 1.11 + 0.0396*kTc - 0.789*kTc**2
else:
    Id_I = 0.20

日散射辐射量Hd与日总的辐射量H之比，与KT的关系

在这里插入图片描述

from sympy import *
H = 1    #某天的日辐照量
H0 = 2    #同一天大气层外日辐照量
KT = H/H0    #日晴空指数
Hd_H = None    #Hd/H
if KT <= 0.17:
    Hd_H  = 0.99
elif KT < 0.75:
    Hd_H = 1.188 - 2.272*KT + 9.473*KT**2 - 21.865*KT**3 + 14.648*KT**4
elif KT < 0.8:
    Hd_H = -0.54*KT + 0.632
else:
    Hd_H = 0.2

小时总辐射量I与全天总辐射量（月平均每日水平辐射量）之比rt

在这里插入图片描述
I = Ib + Id （小时总辐射量 = 小时直射辐射量 + 小时散射辐射量）

from sympy import *
a = 0.409 + 0.5016*sin((Omega_set-60)*pi/180)
b = 0.6609 - 0.4767*sin((Omega_set-60)*pi/180)
rt = pi/24 * (a+b*cos(Omega*pi/180)) * (cos(Omega*pi/180)-cos(Omega_set*pi/180))/(sin(Omega_set*pi/180)-(2*pi*Omega_set/360)*cos(Omega_set*pi/180))

小时散射辐射量与全天散射辐射之比rd

在这里插入图片描述

from sympy import *
rd = pi/24 * (cos(Omega*pi/180)-cos(Omega_set*pi/180))/(sin(Omega_set*pi/180)-(2*pi*Omega_set/360)*cos(Omega_set*pi/180))

倾斜面和水平面上接受到的直射辐照量之比Rb

在这里插入图片描述
在这里插入图片描述

Beta = 20    #表面倾角
Theta = 60    #入射角
Theta_z = 70    #天顶角
Rb = (cos(Theta*pi/180)/cos(Theta_z*pi/180)).evalf()

表面上的总入射辐射IT

在这里插入图片描述
有必要知道或能够估计太阳辐射入射到倾斜表面，如太阳能集热器、窗户或其他被动系统接收器上。入射太阳辐射是一组辐射流的总和，包括直射辐射、天空散射辐射的三个组成部分，以及倾斜表面“看到”的各个表面反射的辐射。这个表面上的总入射辐射可以写成
在这里插入图片描述
在这里插入图片描述
Liu和Jordan（1963）对此模型进行了改进，即各向同性扩散模型。倾斜表面上的辐射包括三个部分：光束、各向同性漫反射和地面漫反射太阳辐射。方程2.14.3中的第三项和第四项为零，因为所有漫射辐射均假定为各向同性。以坡度β从水平面倾斜的表面具有到天空的视角系数Fc−s=（1+cosβ）/2。（如果漫射辐射是各向同性的，这也是Rd，即倾斜表面上的漫射与水平表面上的漫射之比。）倾斜表面对地面的视角系数Fc−g=（1−cosβ）/2，如果周围环境对总太阳辐射的漫反射系数为ρg，则来自周围环境的反射辐射表面应为Iρg（1−cosβ）/2。因此，对方程2.14.3进行了修改，以给出倾斜表面上一小时的总太阳辐射为三项之和：
在这里插入图片描述

from sympy import *
Ib = 1    #小时水平面的直射辐射量 J/m2
Id = 1    #小时水平面的散射辐射量 J/m2
I = Ib + Id    #小时水平面的总辐射
Beta = 20    #表面倾角
Theta = 60    #入射角
Theta_z = 70    #天顶角
Rb = (cos(Theta*pi/180)/cos(Theta_z*pi/180)).evalf()
Rho_g = 0.2    #地面反射率ρg
#倾斜表面上一个小时的总的太阳辐射量I_T，J/m2
I_T = (Ib*Rb + Id*(1+cos((Beta*pi/180)/2))/2 + I*Rho_g*(1-cos((Beta*pi/180)/2))/2).evalf()

倾斜表面上的总辐射IT与水平表面上的总辐射I之比R

在这里插入图片描述
当IT被确定时，倾斜表面上的总辐射与水平表面上的总辐射之比可以被确定。根据定义，
在这里插入图片描述
根据式子（2.15.1），R=IT/I可转化为
在这里插入图片描述

from sympy import *
Ib = 1    #小时水平面的直射辐射量 J/m2
Id = 1    #小时水平面的散射辐射量 J/m2
I = Ib + Id    #小时水平面的总辐射
Beta = 20    #表面倾角
Theta = 60    #入射角
Theta_z = 70    #天顶角
Rb = (cos(Theta*pi/180)/cos(Theta_z*pi/180)).evalf()
Rho_g = 0.2    #地面反射率ρg
#倾斜表面上一个小时的总的太阳辐射量I_T，J/m2
I_T = (Ib*Rb + Id*(1+cos((Beta*pi/180)/2))/2 + I*Rho_g*(1-cos((Beta*pi/180)/2))/2).evalf()
#总的修正因子
R = I_T/I

求水平面直射与散射的例题

英文版

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

中文版

在这里插入图片描述
在这里插入图片描述

代码

#求晴空条件下水平面的小时总辐照量Ic
#若求11:00-12:00的晴天辐照量，用中点11:30来代表。用类似方法对日照时间内每个小时进行详细计算
#麦迪逊，北纬43°，海拔高度270m，8月22日，太阳时11:30AM
from sympy import *
x = symbols('x')
month = 8
day = 22
hour = 11.5    #太阳时为11:30AM
Phi = 43    #纬度φ为43°
A = 0.27    #海拔高度，单位是km
#由于是水平面，所以平面倾角β=0
Beta = 0
#无特别指出，表面方位角为0
Gamma = 0

#n()是一年里的第几天
def n():
    d_num = 0
    for m in range(1, month):
        if m in [1, 3, 5, 7, 8, 10, 12]:
            d_num += 31
        elif m in [2]:
            d_num += 28
        else:
            d_num += 30
    d_num += day
    return d_num
def angle(Theta,Delta,Phi,Beta,Gamma,Omega):
    angle = cos(Theta*pi/180) - sin(Delta*pi/180)*sin(Phi*pi/180)*cos(Beta*pi/180) + sin(Delta*pi/180)*cos(Phi*pi/180)*sin(Beta*pi/180)*cos(Gamma*pi/180)\
            - cos(Delta*pi/180)*cos(Phi*pi/180)*cos(Beta*pi/180)*cos(Omega*pi/180) - cos(Delta*pi/180)*sin(Phi*pi/180)*sin(Beta*pi/180)*cos(Gamma*pi/180)*cos(Omega*pi/180)\
            - cos(Delta*pi/180)*sin(Beta*pi/180)*sin(Gamma*pi/180)*sin(Omega*pi/180)
    return angle

#弦截法迭代求根
def Secant_Method(fx,a,b):
    x0 = a  # 区间下限
    x1 = b  # 区间上限
    x_list = [x1]
    i = 0

    def f(xn):
        f = fx.subs(x,xn)
        return f.evalf()

    while True:
        x2 = x1 - f(x1) * (x1 - x0) / (f(x1) - f(x0))
        x0 = x1
        x1 = x2
        x_list.append(x2.round(1))
        if len(x_list) > 1:
            i += 1
            error = abs((x_list[-1</
最低0.47元/天 解锁文章
Fo*(Bi)
关注
  • 1
    点赞
  • 踩
  • 12
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
专栏目录
PVsyst光伏组件辐照度算法
weixin_64891844的博客
04-18 4200
光伏发电系统—太阳辐照度计算方法 一、太阳辐射能量经光伏组件转化为电能 二、水平面辐射分类 三、水平散射辐射生成算法 四、太阳各辐射分类示意 五、光伏组件直接辐射计算模型 六、 光伏组件倾斜面散射辐射计算-Liu&Jordan模型 七、 光伏组件倾斜面散射辐射计算-Hay模型 八、 光伏组件倾斜面散射辐射计算-Perez模型 九、PVsyst辐射量计算示例 十、 光伏组件地面反射计算模型 ...
基本辐射量总结
zqk的博客
05-24 2194
辐射量 符号 定义 单位 辐射能量 Q 焦耳(J) 辐射通量 Φ ¶ (2) Q/ ¶t(¶ λ) 瓦(W) 辐照度 E ¶ (2) Φ / ¶A (¶ λ) 瓦/米²(W...
【python】计算太阳辐射（GHI、DNI、DHI）
最新发布
a5567899的博客
06-04 590
太阳辐射在许多场景都常有应用，本文将利用pvlib简单计算一下太阳辐射强度。
基于Python的日照时数转太阳辐射计算
03-25
基于Python的日照时数转太阳辐射计算 自动化程度高，计算速度快
Python计算日照总辐射量
毛发浓密的猿工
09-08 783
【代码】Python计算日照总辐射量。
光伏电站+功率+辐照度度+天气
01-18
（1）用于学习研究，间隔15分钟的数据 （2）包含电站率曲线 （3）包含实时天气，温度，湿度，辐照度
地外水平/斜面的太阳理论辐照度计算及其python实现
ysk937111075的博客
05-11 5165
本文讲解了地外太阳辐照度的计算，并进行了python实现
根据经纬度和日期计算日出日落时间
06-28
在IT领域，根据经纬度和日期计算日出日落时间是一项常见的地理信息处理任务，它涉及到天文学、地球物理学和编程技术。这个程序的核心是利用天文学公式来预测太阳在地平线上的位置，从而确定一天中的日出和日落时刻。...
日出日落时间计算
11-09
计算某地某天的日出日落时间、太阳正午时间，以及正午太阳高度角
全国840个气象站点日照时数转日总太阳辐射python代码.zip
09-13
本文档为通过全国气象观测站点观测的日照时数ssd计算日总太阳辐射sr的python脚本，简单易懂。内附输入数据及输出结果文件示例。
全国840个气象站点日照时数转日总太阳辐射python代码
12-10
本文档为通过全国气象观测站点观测的日照时数ssd计算日总太阳辐射sr的python脚本，简单易懂。内附输入数据及输出结果文件示例。
不同倾斜面上太阳辐射量的计算
02-01
采用Klein.S.A和Theilacker.J.C所提出的月平均太阳辐射量的计算方法，计算倾斜面上的太阳辐射量
基于Python实现对房价的预测项目源码.zip
04-21
基于Python实现对房价的预测项目源码.zip项目背景：网上有条段子，某地房价5w每平，月收入刚好过万，掐指一算，命中注定买房是不可能的，这辈子都不可能买房，所以要定个小目标：“我真的还想再活500年······...
pysolar:Pysolar是Python库的集合，用于模拟太阳对地球上任何点的照射。 它包括用于非常精确的星历表计算的代码
05-14
它包含用于极其精确的星历表计算的代码，等等。 注意：目前，Pysolar的最新提交不适用于Python 2.x 版本0.6适用于2.x： ： 但0.7和更高版本有很多更改。 它们已经针对Python 3.4进行了验证，但是3.2版或更早的版本...
基于Python和机器学习开的网络爬虫获取房价信息、数据的预处理和可视化+源代码+文档说明
01-16
基于Python：网络爬虫获取房价信息、数据的预处理和可视化、搭建基于房价预测的机器学习模型、房价预测。 -数据来源：房天下 项目背景：网上有条段子，某地房价5w每平，月收入刚好过万，掐指一算，命中注定买房是不...
python--绘制WRF模式近地面风场以及辐射
【更新于公众号 | 简谱学记】
10-13 4369
读取WRF风速资料U、V，以及辐射（辐射包括三部分）将读取的资料进行绘图，并将绘图后的结果自动 保存到指定路径下，方面后续绘制动图。
2023数学建模美赛E题思路
2201_75432540的博客
02-17 2571
1.光污染相关的数学模型光污染是指人类对夜间光照环境进行不必要、过量、不合理使用而导致的环境问题。为了量化和描述光污染的严重程度，人们提出了多种光污染度量标准数学模型，以下是其中一些常用的模型:光度学模型:该模型考虑到人眼对光线的感受性，通过计算夜间照度值和色温等指标来评估光污染的严重程度。辐射度模型:该模型基于辐射度量标准,通过测量夜间照度和夜间光线的光谱分布来确定光污染水平。
python根据日照天数计算太阳辐射量原理
06-05
太阳辐射量是指太阳对地球各个区域所辐射的能量，而...总之，Python通过采集数据、处理数据和计算输出等步骤，能够很好地实现根据日照天数计算太阳辐射量的功能，而该功能在气象、农业、环保等多个领域有广泛的应用。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值