CN111291318B - 一种基于辐射通量密度的三维云层辐射计算方法 - Google Patents
一种基于辐射通量密度的三维云层辐射计算方法 Download PDFInfo
- Publication number
- CN111291318B CN111291318B CN202010129747.2A CN202010129747A CN111291318B CN 111291318 B CN111291318 B CN 111291318B CN 202010129747 A CN202010129747 A CN 202010129747A CN 111291318 B CN111291318 B CN 111291318B
- Authority
- CN
- China
- Prior art keywords
- radiation
- cloud
- cloud layer
- flux density
- horizontal
- 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
- 230000005855 radiation Effects 0.000 title claims abstract description 214
- 230000004907 flux Effects 0.000 title claims abstract description 76
- 238000004364 calculation method Methods 0.000 title claims abstract description 66
- 238000000034 method Methods 0.000 claims abstract description 41
- 230000008859 change Effects 0.000 claims abstract description 34
- 230000000191 radiation effect Effects 0.000 claims abstract description 19
- 230000005540 biological transmission Effects 0.000 claims abstract description 18
- 238000013507 mapping Methods 0.000 claims abstract description 14
- 230000005571 horizontal transmission Effects 0.000 claims abstract description 6
- 230000000694 effects Effects 0.000 claims description 12
- 238000012937 correction Methods 0.000 claims description 8
- 230000003287 optical effect Effects 0.000 claims description 6
- 230000010354 integration Effects 0.000 claims description 3
- 230000003313 weakening effect Effects 0.000 claims description 3
- 210000004027 cell Anatomy 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000008033 biological extinction Effects 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 210000004692 intercellular junction Anatomy 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000032823 cell division Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于辐射通量密度的三维云层辐射计算方法,包括如下步骤:首先利用高斯方向和大气属性构建水平辐射通量密度传输模型,然后利用爱丁顿近似法和严格单次散射计算式计算三维云层基础辐射场;考虑邻云辐射效应对水平传输方向上的三维云层基础辐射场的影响,将其分为邻云影响辐射增强和邻云影响辐射减弱,计算由邻云辐射效应影响所产生的水平辐射通量密度变化;将水平辐射通量密度变化映射回高斯方向,再将三维云层基础辐射场与变化辐射强度结合得到更新后三维云层辐射强度,并由其计算多次散射源函数;结合由初始条件计算的单次散射源函数与热源函数,得到三维云层辐射源函数。
Description
技术领域
本发明涉及大气辐射计算领域,尤其是一种基于辐射通量密度的三维云层辐射计算方法。
背景技术
云是自然界中十分复杂的物理现象之一,在大气辐射传输过程中产生着较大的影响。由于云在三维空间内的分布的不均匀性,使得三维云层场景的辐射传输问题成为大气辐射传输领域中的难题。为了研究有云大气中的辐射传输问题,需要充分考虑云的三维结构和辐射特性。
三维大气辐射传输的精确计算,相对应的是成本的大幅度提高。其中,三维空间内多次散射的精确计算很大程度上影响着最终计算成本。通常可以通过多次迭代方法求解多次散射源函数,复杂的多次迭代过程势必会导致较高的成本和较慢的运算速度。因此,合理地简化多次散射源函数的计算过程来降低所需成本并提高运算速度,成为大气辐射传输领域亟需研究讨论的问题。
发明内容
为了解决上述背景技术中的问题,本发明旨在提供一种基于辐射通量密度的三维云层辐射计算方法,能够在保持较高计算精度的前提下,加快计算三维云层辐射源函数,降低计算成本,并且可以很好地应用于三维多云场景。
为解决上述技术问题,本发明提供一种基于辐射通量密度的三维云层辐射计算方法,包括如下步骤:
(1)利用高斯方向和大气属性构建水平辐射通量密度传输模型;
(2)利用爱丁顿近似法和严格单次散射计算式计算三维云层基础辐射场;
(3)考虑邻云辐射效应对水平传输方向上的三维云层基础辐射场的影响,将其分为邻云影响辐射增强和邻云影响辐射减弱,计算由邻云辐射效应影响所产生的水平辐射通量密度变化;
(4)将水平辐射通量密度变化映射回高斯方向,再将三维云层基础辐射场与变化辐射强度结合得到更新后三维云层辐射强度,并由其计算多次散射源函数;
(5)结合由初始条件计算的单次散射源函数与热源函数,得到三维云层辐射源函数。
优选的,步骤(1)中,利用高斯方向和大气属性构建水平辐射通量密度传输模型,具体步骤如下:
(1)将三维云层场景分为所需大小尺度的网格单元,每个单元都有其相邻的四个单元(如图2),利用大气属性文件,将大气属性根据空间位置存储在每个单元内;
(2)基于离散坐标下的高斯积分,将三维云层辐射场离散到已设置个数的高斯方向上,且同一方位角的高斯方向总和构成一个“三棱柱”空间(如图3);
(3)将高斯方向上的三维云层辐射场映射到水平辐射传输方向上。
优选的,步骤(3)中,考虑邻云辐射效应对水平传输方向上的三维云层基础辐射场的影响,将其分为邻云影响辐射增强和邻云影响辐射减弱,具体步骤如下:
(1)基于高斯方向上的由爱丁顿近似法和严格单次散射计算式计算得到的三维云层基础辐射场,将高斯方向上的辐射强度映射到水平辐射通量密度传输方向上;
(2)将邻云辐射效应分为邻云影响辐射增强和邻云影响辐射减弱即I+(τ;Ω′ts)和I-(τ;Ω′ts),得到水平辐射通量密度的增强值或减弱值;
(3)考虑太阳天顶角变化导致的辐射误差,对水平辐射通量密度的增强值和减弱值进行修正。
优选的,步骤(3)中,计算由邻云辐射效应影响所产生的水平辐射通量密度变化,具体步骤如下:
(1)选择三维云层场景内第一个单元A,再搜索与其水平方向相邻的四个单元B1、B2、B3和B4;
(2)将基础辐射场映射到水平传输方向,计算映射系数ψts;
ψts=cos(Δφts)sin(θts)
其中θts高斯方向的天顶角;Δφts高斯方向在XY平面上的投影与其所处ts空间内水平坐标轴方向的平面夹角,Δφts的取值范围是-45°<Δφts<45°。
(3)计算水平辐射通量密度在单元A与其水平相邻单元B之间的反射比例,即反照率R;
其中,τ为光学厚度,ω为单次散射反照率,g为不对称因子。
(4)计算单元A的水平辐射通量密度受邻云辐射影响的变化I+(τ;Ω′ts)和I-(τ;Ω′ts);
其中I+(τ;Ω′t1)是Ω′t1方向水平辐射通量密度的增加值;I-(τ;Ω′t1)是Ω′t1方向水平辐射通量密度的减弱值;Ω′t2方向是与Ω′t1相同高斯方向相反方位的方向(如图3);R1是水平辐射通量密度在A与B1的交界处发生反射的比例;R2是水平辐射通量密度在A与B2的交界处发生反射的比例。
(5)计算单元A由于邻云辐射效应影响所产生的辐射射通量密度变化ICE(τ;Ω′ts),计算修正系数ACE;
ACE=exp(sinθ0-sinθE)-1
ICE(τ;Ω′ts)=(1+ACE)I+(τ;Ω′ts)-(1-ACE)I-(τ;Ω′ts)
其中θ0为太阳天顶角;θE为二流近似法计算辐射通量所使用的高斯方向天顶角;ICE(τ;Ω′ts)代表某一单元处Ω′ts方向由邻云辐射效应所产生的水平辐射通量密度变化值,I+(τ;Ω′ts)为Ω′ts方向水平辐射通量密度的增加值,I-(τ;Ω′ts)为Ω′ts方向水平辐射通量密度的减弱值,ACE为弥补太阳方向变化对模型带来误差的修正系数。
(6)重复步骤1-5,直到选则到三维云层场景内最后一个单元A,并且单元A的ICE(τ;Ω′ts)都计算完成为止。
优选的,步骤(4)中,将水平辐射通量密度变化映射回高斯方向,再将三维云层基础辐射场与变化辐射强度结合得到更新后三维云层辐射强度,并由其计算多次散射源函数,具体步骤如下:
(1)在三维云层场景内选取第一个单元A;
(2)计算从水平通量密度方向到高斯方向的分配权重P(Ω′ts;Ωs)和高斯积分权重Ats(Ω′ts),并由此计算高斯方向变化辐射强度I'CE(τ;Ω′):
ICE'(τ;Ω′)=P(Ω′ts;Ωs)Ats(Ω′ts)ICE(τ;Ω′)
(3)利用爱丁顿近似法计算和严格单次散射计算式计算得到单元A的三维云层基础辐射场,结合变化辐射强度得到更新后三维云层辐射强度:
其中IEMS(τ;Ω′)为由爱丁顿近似计算得到三维云层近似多次散射强度,ISS(τ;Ω′)为由严格单次散射计算式计算得到的三维云层单次散射强度,I”(τ;Ω′)为以三维云层基础辐射场I′(τ;Ω′)作为输入计算得到的更新后三维云层辐射强度;
(4)利用更新后三维云层辐射亮度计算单元A处的多次散射源函数;
(5)重复步骤1-4,直到选取了三维云层场景内最后一个单元作为A点,且该处的多次散射源函数项计算完成为止。
本发明的有益效果为:本发明可在保持较高计算精度的前提下,加快计算三维云层辐射源函数,降低计算成本,并且可以很好地应用于三维多云场景。
附图说明
图1为本发明的方法流程示意图;
图2为本发明的水平方向单元划分示意图;
图3为本发明的高斯方向构成的三棱柱空间示意图;
图4为本发明的应用于所设云场景1随太阳天顶角变化的改进辐射源函数计算法和爱丁顿近似法与真实辐射值的相对误差结果图;
图5为本发明的应用于所设云场景2随太阳天顶角变化的改进辐射源函数计算法和爱丁顿近似法与真实辐射值的相对误差结果图。
具体实施方式
现将结合附图对本发明的技术方案进行完整的描述。以下描述仅仅是本发明的一部分实施案例而已,并非全部。基于本发明中的实施案例,本领域技术人员在没有作出创造性劳动的前提下所获得的所有其他实施案例,都属于本发明的权利保护范围之内。
如图1所示,一种基于辐射通量密度的三维云层辐射计算方法,包括如下步骤:
(1)利用高斯方向和大气属性构建水平辐射通量密度传输模型;
(2)利用爱丁顿近似法和严格单次散射计算式计算三维云层基础辐射场;
(3)考虑邻云辐射效应对水平传输方向上的三维云层基础辐射场的影响,将其分为邻云影响辐射增强和邻云影响辐射减弱,计算由邻云辐射效应影响所产生的水平辐射通量密度变化;
(4)将水平辐射通量密度变化映射回高斯方向,再将三维云层基础辐射场与变化辐射强度结合得到更新后三维云层辐射强度,并由其计算多次散射源函数;
(5)结合由初始条件计算的单次散射源函数项与热源函数项,得到三维云层辐射源函数。
1.邻云辐射效应
首先将三维云层场景分为所需大小尺度的网格单元,每个单元都有其相邻四个单元(如图2),再由已知的云层含水量分布及米氏散射结果生成大气属性文件,计算对应网格位置处三维云层的反照率、消光系数和相函数,再将以上三类参量的大气属性存储在它们所对应位置的单元中。
通过爱丁顿近似计算得到多次散射项、严格单次散射计算式得到的单次散射项之后求得三维云层基础辐射场,再基于邻云辐射效应,计算得到三维云层变化辐射强度,结合三维云层基础辐射场和变化辐射强度,再通过辐射强度到源函数的转换,可以实现快速计算三维云层场景辐射源函数。
邻云辐射效应是由于光和热等辐射量在云层与邻云之间互相作用和交换而产生的辐射场分布特性改变,该效应可以准确而全面地对非均匀性较强的云层辐射进行描述。故通过计算与辐射相关参量的具体变化数值,可以对邻云辐射效应进行量化表达。三维云层空间内每个单元处都储存着由属性文件计算得到的辐射通量密度,考虑邻云辐射效应对辐射通量密度造成的影响,在网格单元内由影响造成的辐射通量密度变化分为邻云影响辐射增强和邻云影响辐射减弱(如图1):
ICE(τ;Ω′ts)=(1+ACE)I+(τ;Ω′ts)-(1-ACE)I-(τ;Ω′ts)
其中ICE(τ;Ω′ts)代表某一单元处Ω′ts方向由邻云辐射效应所产生的水平辐射通量密度变化值,I+(τ;Ω′ts)为Ω′ts方向水平辐射通量密度的增加值,I-(τ;Ω′ts)为Ω′ts方向水平辐射通量密度的减弱值,ACE为弥补太阳方向变化对模型带来误差的修正系数。
1.1映射到水平方向
对于三维云层辐射问题,邻云辐射效应可以表现为光束在云层之间传播路径的改变以及由此而产生的辐射能量变化,由邻云辐射效应所产生的水平辐射通量密度变化,可以计算多次散射源函数的近似值。辐射通量密度的方向为空间内的已设置高斯方向,分析邻云影响产生的水平辐射通量密度变化,需要将高斯方向辐射映射到水平方向上来,由高斯方向到水平方向的变换系数ψts:
ψts=cos(Δφts)sin(θts)
其中θts为高斯方向的天顶角;Δφts为高斯方向在XY平面上的投影与其所处ts空间内水平坐标轴方向的平面夹角,Δφts的取值范围是-45°<Δφts<45°。
1.2邻云辐射影响增强和减弱
首先计算三维云层空间的基础辐射场,再找到第一个单元A,由其基础辐射强度及其相邻的单元B处的基础辐射强度计算得到受邻云辐射影响增强的辐射通量密度和减弱的辐射通量密度:
其中I+(τ;Ω′t1)是Ω′t1方向水平辐射通量密度的增加值;I-(τ;Ω′t1)是Ω′t1方向水平辐射通量密度的减弱值;Ω′t2方向是与Ω′t1相同高斯方向相反方位的方向(如图3);R1是水平辐射通量密度在A与B1的交界处发生反射的比例;R2是水平辐射通量密度在A与B2的交界处发生反射的比例。
2辐射通量密度变化值计算
2.1修正系数计算
考虑太阳方向变化对模型带来的影响,从而弥补方法模型中存在的输入性误差。修正系数计算式如下:
ACE=exp(sinθ0-sinθE)-1
其中θ0为太阳天顶角;θE为二流近似法计算辐射通量所使用的高斯方向天顶角。
由于高斯积分的标准计算式可以写为:
其中f(γ)被积分函数;γ积分变量;aj高斯积分权重;γj高斯积分点。且在爱丁顿近似法等二流近似方法中,m的取值为1,此时高斯积分权重a1的取值为0.5773503,其为该高斯方向天顶角的余弦,因此可以求得(6)式中θE的数值约为54.74°。
2.2反照率R模型
找到三维空间网格内某一单元A,当在水平方向尺度上对邻云辐射效应进行分析时,考虑与其水平相邻的单元B1、B2、B3和B4,在三维空间含云区域单元交界处会发生反射,反照率R表示单元交界处发生反射的比例。当τ在0-20时,反照率R*的计算参考邻云间反射率计算方法和二流近似中的反照率计算方法得到:
其中,τ为光学厚度;ω为单次散射反照率;g为不对称因子。
当τ不断增大时,R*将逐渐趋近于1,此时无论ω和g取值如何,反照率的计算结果都会发生“饱和”,R*的计算结果将会严重偏离真实值,针对该情况,对R*进行修正得到适用于任意光学厚度的反照率R:
不对称因子可以描述前向散射和后向散射的程度,其计算方法为:
其中,Ω为散射角;P(Ω)为散射相函数。
光学厚度描述辐射传输过程中经过云层或大气吸收和散射产生的总衰弱:
其中,k为消光系数;s为路径距离。
2.3权重系数计算
将水平辐射通量密度分配至相应的高斯方向处,以实现该高斯方向上辐射强度变化的更新。辐射场辐射场从水平通量密度方向到高斯方向的分配权重为P(Ω′ts;Ωs)Ats(Ω′ts)。
其中,P(Ω′ts;Ωs)为水平方向到高斯方向的夹角作为散射角方式的散射相函数;Ats(Ω′ts)为高斯积分权重。
2.4映射回高斯方向
根据得到的邻云辐射影响增强和邻云辐射影响减弱,计算单元A处的水平方向变化辐射通量密度:
ICE(τ;Ω′ts)=(1+ACE)I+(τ;Ω′ts)-(1-ACE)I-(τ;Ω′ts)
其中ICE(τ;Ω′ts)代表某一单元处Ω′ts方向由邻云辐射效应所产生的水平辐射通量密度变化值,I+(τ;Ω′ts)为Ω′ts方向水平辐射通量密度的增加值,I-(τ;Ω′ts)为Ω′ts方向水平辐射通量密度的减弱值。
高斯方向上辐射通量密度更新为:
ICE'(τ;Ω′)=P(Ω′ts;Ωs)Ats(Ω′ts)ICE(τ;Ω′)
于是,可以计算得到单元A处由于邻云辐射效应影响得到的变化辐射通量密度ICE',接着选取新的A点,重复1与2所述计算过程,直到选取了三维云层场景内最后一个单元作为A点,即所有单元的ICE′都计算完成为止。
3.三维云层源函数计算
在所有A点的ICE′都计算完成后,将三维云层基础辐射场与其结合,得到更新后辐射强度,并利用其计算多次散射源函数,进而结合单次散射源函数项与热源函数项,得到三维云层空间总辐射源函数。
3.1辐射强度累加
多次散射源函数近似计算模型以基础辐射场I′(τ;Ω′)作为输入,基础辐射场由通过严格单次散射计算式得到的单次散射强度和利用爱丁顿近似法得到的多次散射强度构成。接下来,计算单元A处的更新后辐射强度:
I′(τ;Ω′)=ISS(τ;Ω′)+IEMS(τ;Ω′)
I″(τ;Ω′)=I′(τ;Ω′)+ICE'(τ;Ω′)
其中ISS(τ;Ω′)通过严格单次散射计算式得到的单次散射强度;IEMS(τ;Ω′)通过爱丁顿近似法计算得到的多次散射强度。
3.2多次散射函数项计算
利用邻云辐射效应影响的变化辐射亮度更新后的辐射强度I”来计算单元A处的多次散射函数JMS(τ;Ω):
其中P(Ω;Ω′)为三维云层空间单元A处的散射相函数。
3.3辐射源函数计算
在多次散射源函数JMS计算完成后,将其与单次散射源函数JSS和热源函数JT结合,得到辐射源函数J:
其中ω为反照率;F0为从辐射源产生并直接照射在单元上的辐射亮度;τ为大气光学厚度;B为普朗克公式;T为温度。
选取新的单元A,重复3.1、3.2与3.3所述计算过程,直到选取了三维云层场景内最后一个单元A,且该处的总辐射源函数计算完成为止。
4.辐射计算方法应用
为了验证本文提出的基于辐射通量密度的三维云层辐射计算方法,选用了真实云场景进行实验。该云场景为I3RC Phase Ⅱ的层积云场,网格尺寸为64*64*16,空间分辨率为55m*55m*25m,所处空间范围为x:0-3.52km,y:0-3.52km,z:10.1-11.6km。为测试方法在云层场景中的适用情况,截取了其中5*5*5的网格点数且云覆盖率较高的两块区域进行实验。考虑实验计算成本,将高斯方向离散数设为4,方位角离散数设为4,即在16个高斯方向上描述整个辐射场。
基于实验环境,将本文改进辐射源函数计算方法与爱丁顿近似法、SHDOM法进行比较分析。爱丁顿近似法是二流近似法的一种,其一维计算模式不考虑水平辐射变化,本方法在其基础上考虑水平辐射交换,与其比较可以明显观察到本方法的改进效果;而SHDOM方法虽然耗费成本较高,但计算结果趋于真实值,可以将其视为自然场景真实辐射值,与其比较能够分析本方法的误差。
在截取的云场景1中,云覆盖率为93.3%,在层积云场所处空间范围为x:1.98-2.2km,y:2.2-2.42km,z:10.9-11.3km。在截取的云场景2中,云覆盖率为69.44%,在层积云场所处空间范围为x:2.585-2.805km,y:2.805-3.025km,z:10.7-11.1km。通过三种方法的辐射计算,得到在不同太阳天顶角情况下爱丁顿近似法和改进辐射源函数计算法同真实辐射值的相对误差结果图。
结果中展示在上行空间本文辐射计算方法的平均相对误差明显低于爱丁顿近似法,云层场景1的平均精度提升为22.5%,云层场景2的平均精度提升为15%。对于下行空间而言,下行空间中存在太阳直射方向的前向方向,前向散射中包含了云层中的大多数单次散射能量,故下行空间内的辐射值大多高于上行空间,因此无论是爱丁顿近似还是本文方法,计算得到的的相对误差都较小,爱丁顿近似法与本方法皆趋近于真实辐射情况,本文方法的平均相对误差基本低于爱丁顿近似法,云层场景1的平均精度提升为2.8%,云层场景2的平均精度提升为1%。
在整个三维云层空间内,云层场景1的平均相对误差为8.75%,云层场景2的平均相对误差为11%,即在两种云场景中,方法误差皆在可接受范围内,可以满足一般应用需求。
本发明可在保持较高计算精度的前提下,加快计算三维云层辐射源函数,降低计算成本,并且可以很好地应用于三维多云场景。
以上实施例仅供说明本发明之用,而非对本发明的限制,有关技术领域的技术人员,在不脱离本发明的精神和范围的情况下,所作出各种变换或变型,均属于本发明的范畴。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前下,本发明还会有各种变化和改进,本发明要求保护范围由所附的权利要求书、说明书及其等效物界定。
Claims (1)
1.一种基于辐射通量密度的三维云层辐射计算方法,其特征在于,包括如下步骤:
1)利用高斯方向和大气属性构建水平辐射通量密度传输模型;步骤1)中,利用高斯方向和大气属性构建水平辐射通量密度传输模型,具体为:
步骤1.1,将三维云层场景分为所需大小尺度的网格单元,每个单元都有其相邻的四个单元,利用大气属性文件,将大气属性根据空间位置存储在每个单元内;
步骤1.2,基于离散坐标下的高斯积分,将三维云层辐射场离散到已设置个数的高斯方向上,且同一方位角的高斯方向总和构成一个三棱柱空间;
步骤1.3,将高斯方向上的三维云层辐射场映射到水平辐射传输方向上;
2)利用爱丁顿近似法和严格单次散射计算式计算三维云层基础辐射场;
3)将所述三维云层基础辐射场分为邻云影响辐射增强和邻云影响辐射减弱,计算由邻云辐射效应影响所产生的水平辐射通量密度变化;步骤3)中,将所述三维云层基础辐射场分为邻云影响辐射增强和邻云影响辐射减弱具体为:
步骤3.1,将高斯方向上的辐射强度映射到水平辐射通量密度传输方向上;
步骤3.2,将邻云辐射效应分为邻云影响辐射增强和邻云影响辐射减弱即I+(τ;Ω′ts)和I-(τ;Ω′ts),得到水平辐射通量密度的增强值或减弱值;
步骤3.3,对水平辐射通量密度的增强值和减弱值进行修正,以消除太阳天顶角变化导致的辐射误差;
步骤3)中,计算由邻云辐射效应影响所产生的水平辐射通量密度变化,具体为:
步骤3.4,选择三维云层场景内任一单元A,再搜索与其水平方向相邻的四个单元B1、B2、B3和B4;
步骤3.5,将基础辐射场映射到水平传输方向,计算映射系数ψts;
步骤3.6,计算水平辐射通量密度在单元A与其水平方向相邻的四个单元B1、B2、B3和B4之间的反射比例,即反照率R;
步骤3.7,计算单元A的水平辐射通量密度受邻云影响辐射增强和邻云影响辐射减弱I+(τ;Ω′ts)和I-(τ;Ω′ts);
步骤3.8,计算单元A由于邻云辐射效应影响所产生的辐射通量密度变化ICE(τ;Ω′ts),计算修正系数ACE;
ACE=exp(sinθ0-sinθE)-1
ICE(τ;Ω′ts)=(1+ACE)I+(τ;Ω′ts)-(1-ACE)I-(τ;Ω′ts)
其中θ0为太阳天顶角;θE为二流近似法计算辐射通量所使用的高斯方向天顶角;ICE(τ;Ω′ts)代表某一单元处Ω′ts方向由邻云辐射效应所产生的水平辐射通量密度变化值,I+(τ;Ω′ts)为Ω′ts方向邻云影响辐射增强,I-(τ;Ω′ts)为Ω′ts方向邻云影响辐射增强和邻云影响辐射减弱,ACE为弥补太阳方向变化对模型带来误差的修正系数,τ为光学厚度,g为不对称因子;
步骤3.9,重复步骤3.4至3.8,直到选择到三维云层场景内最后一个单元A,并且单元的ICE(τ;Ω′ts)都计算完成为止;
4)将水平辐射通量密度变化映射回高斯方向,再将三维云层基础辐射场与变化辐射强度结合得到更新后三维云层辐射强度,并由其计算多次散射源函数;步骤4)中,将水平辐射通量密度变化映射回高斯方向,再将三维云层基础辐射场与变化辐射强度结合得到更新后三维云层辐射强度,并由其计算多次散射源函数,具体为:
步骤4.1,选择三维云层场景内任一个单元A;
步骤4.2,将单元A从水平通量传输方向映射到离散高斯方向,计算分配权重P(Ω′ts;Ωs)和高斯积分权重Ats(Ω′ts),并由此计算高斯方向变化辐射强度I′CE(τ;Ω′):
ICE′(τ;Ω′)=P(Ω′ts;Ωs)Ats(Ω′ts)ICE(τ;Ω′)
步骤4.3,利用爱丁顿近似法和严格单次散射计算式计算得到单元A的三维云层基础辐射场,结合高斯方向变化辐射强度I′CE(τ;Ω′)得到更新后三维云层辐射强度:
其中IEMS(τ;Ω′)为由爱丁顿近似计算得到的三维云层近多次似散射强度,ISS(τ;Ω′)为由严格单次散射计算式计算得到的三维云层单次散射强度,I″(τ;Ω′)为以三维云层基础辐射场I′(τ;Ω′)作为输入计算得到的更新后的三维云层辐射强度;
步骤4.4,利用更新后三维云层辐射强度计算单元A处的多次散射源函数;
步骤4.5,重复步骤4.1至4.4,直到选择到三维云层场景内最后一个单元A,且单元处的多次散射源函数计算完成为止;
5)结合由初始条件计算的单次散射源函数与热源函数,得到三维云层辐射源函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010129747.2A CN111291318B (zh) | 2020-02-28 | 2020-02-28 | 一种基于辐射通量密度的三维云层辐射计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010129747.2A CN111291318B (zh) | 2020-02-28 | 2020-02-28 | 一种基于辐射通量密度的三维云层辐射计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111291318A CN111291318A (zh) | 2020-06-16 |
CN111291318B true CN111291318B (zh) | 2023-12-19 |
Family
ID=71025691
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010129747.2A Active CN111291318B (zh) | 2020-02-28 | 2020-02-28 | 一种基于辐射通量密度的三维云层辐射计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111291318B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112348872B (zh) * | 2020-09-28 | 2023-10-03 | 西安电子科技大学 | 一种基于物理过程的红外积云的建模方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110826015A (zh) * | 2019-09-25 | 2020-02-21 | 南京航空航天大学 | 一种基于球谐离散坐标法的三维含云大气辐射计算方法 |
-
2020
- 2020-02-28 CN CN202010129747.2A patent/CN111291318B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110826015A (zh) * | 2019-09-25 | 2020-02-21 | 南京航空航天大学 | 一种基于球谐离散坐标法的三维含云大气辐射计算方法 |
Non-Patent Citations (1)
Title |
---|
章超 等.三维有云大气辐射变温率数值模拟研究――邻云辐射效应.气象学报.2016,(第04期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN111291318A (zh) | 2020-06-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yoon et al. | Evaluation of hourly solar radiation on inclined surfaces at Seoul by Photographical Method | |
CN111291318B (zh) | 一种基于辐射通量密度的三维云层辐射计算方法 | |
CN107341284B (zh) | 高精度预测低频电波传播特性的双向抛物方程方法 | |
Ferragut et al. | Comparison between 2.5-D and 3-D realistic models for wind field adjustment | |
CN107632964A (zh) | 一种平面地磁异常场向下延拓递归余弦变换法 | |
CN107391880A (zh) | 一种精确高效的天线罩电性能分析方法 | |
Nikolaeva et al. | The influence of neighbouring clouds on the clear sky reflectance studied with the 3-D transport code RADUGA | |
CN104835203B (zh) | 一种基于OptiX的自然场景温度场计算方法及系统 | |
CN106874690B (zh) | 考量水平传播因素的地形重力波拖曳参数化方法 | |
Ghassemi et al. | Multilevel quasidiffusion method with mixed-order time discretization for multigroup thermal radiative transfer problems | |
CN111859704A (zh) | 一种分布式多视角下非刚体目标电磁散射建模方法 | |
CN114898040A (zh) | 基于硬件管线加速的实时红外成像仿真方法及系统 | |
CN106471504A (zh) | 模拟在铸造过程中热耦合表面辐射的影响的方法和算法 | |
Turner et al. | MPAS-Seaice (v1. 0.0): Sea-ice dynamics on unstructured Voronoi meshes | |
CN106548509A (zh) | 一种基于cuda及三维成像载荷的三维影像生成方法 | |
Berkova et al. | Temperature effect of muon component and practical questions of how to take into account in real time | |
Lian et al. | Real‐time temperature estimation method for electromagnetic performance improvement of a large axisymmetric radio telescope under solar radiation | |
CN107748733B (zh) | 适用于云微物理特性连续变化的辐射传输计算方法 | |
CN110826015B (zh) | 一种基于球谐离散坐标法的三维含云大气辐射计算方法 | |
CN114187399A (zh) | 一种地外星体表面光照仿真影像生成方法和系统 | |
Li et al. | A fast piecewise 3D ray tracing algorithm for determining slant total delays | |
CN106796736B (zh) | 用于模拟或计算过程中的射线追踪方法 | |
CN108957454B (zh) | 机载干涉合成孔径雷达的干涉相位仿真方法 | |
KR20180045401A (ko) | 지형효과를 고려한 태양에너지가 반영된 기상 예측 시스템 | |
Meller | Analytically calculating shading in regular arrays of sun-pointing collectors |
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 |