[發(fā)明專利]三維蒙特卡洛輻射傳輸模式改進(jìn)方法、存儲介質(zhì)及服務(wù)器

有效

• 申請?zhí)枺?/td>: 201810793386.4
• 申請日: 2018-07-18
• 公開（公告）號: CN109039455B
• 公開（公告）日: 2020-12-04
• 發(fā)明（設(shè)計(jì)）人: 王蕾迪;呂達(dá)仁;胡飛
• 申請（專利權(quán)）人: 華南農(nóng)業(yè)大學(xué)
• 主分類號: H04B10/11
• 分類號: H04B10/11
• 代理公司: 佛山市廣盈專利商標(biāo)事務(wù)所(普通合伙) 44339
• 代理人: 李俊

1.一種基于RPV模型的三維蒙特卡洛輻射傳輸模式改進(jìn)方法，其特征在于，所述方法包括：

采用ρ₀,κ,α三個(gè)參數(shù)來對地表進(jìn)行描述，其中ρ₀為地表反射太陽光照射的強(qiáng)度，κ為地表各向異性程度，α為控制太陽光照射相對前向或后向的散射量；

根據(jù)地表描述的ρ₀,κ,α三個(gè)參數(shù)進(jìn)行三個(gè)函數(shù)計(jì)算，獲取計(jì)算結(jié)果；

根據(jù)所述計(jì)算結(jié)果構(gòu)建雙向反射率因子模型，獲取構(gòu)建好的雙向反射率因子模型；

將所述構(gòu)建好的雙向反射率因子模型耦合至三維蒙特卡洛輻射傳輸模式的原始地表模塊的surfaceProperties.f95運(yùn)行文件中，對所述三維蒙特卡洛輻射傳輸模式進(jìn)行改進(jìn)，獲取改進(jìn)的三維蒙特卡洛輻射傳輸模式；

將所述構(gòu)建好的雙向反射率因子模型耦合至三維蒙特卡洛輻射傳輸模式的原始地表模塊的surfaceProperties.f95運(yùn)行文件中，對所述三維蒙特卡洛輻射傳輸模式進(jìn)行改進(jìn)，獲取改進(jìn)的三維蒙特卡洛輻射傳輸模式，包括：

將所述構(gòu)建好的雙向反射率因子模型以FORTRAN代碼的形式進(jìn)行編寫，獲取編寫好的FORTRAN代碼模型；

將所述編寫好的FORTRAN代碼添加至三維蒙特卡洛輻射傳輸模式的原始地表模塊的surfaceProperties.f95運(yùn)行文件中，對所述三維蒙特卡洛輻射傳輸模式進(jìn)行改進(jìn)，獲取改進(jìn)的三維蒙特卡洛輻射傳輸模式；

所述方法還包括：

在所述改進(jìn)的三維蒙特卡洛輻射傳輸模式中的積分器模塊的子程序monteCarloRadiativeTransfer.f95中，以FORTRAN語言形式，加入調(diào)用原始地表模塊中的surfaceProperties.f95運(yùn)行文件的call調(diào)用命令，及修改FORTRAN代碼；

原始地表模塊的surfaceProperties.f95運(yùn)行文件中耦合了相應(yīng)的構(gòu)建好的雙向反射率因子模型，其調(diào)用命令有所不同，需要以FORTRAN語言形式進(jìn)行調(diào)用原始地表模塊中的surfaceProperties.f95運(yùn)行文件中的call調(diào)用命令的相關(guān)編程，在編好之后，加入積分器模塊中monteCarloRadiativeTransfer.f95運(yùn)行程序中，方便后續(xù)積分器模塊中monteCarloRadiativeTransfer.f95運(yùn)行程序?qū)Ω倪M(jìn)的三維蒙特卡洛輻射傳輸模式中的原始地表模塊中的surfaceProperties.f95運(yùn)行文件的調(diào)用。

2.根據(jù)權(quán)利要求1所述的三維蒙特卡洛輻射傳輸模式改進(jìn)方法，其特征在于，所述根據(jù)地表描述的ρ₀,κ,α三個(gè)參數(shù)進(jìn)行三個(gè)函數(shù)計(jì)算的三個(gè)函數(shù)分別如下：

M(κ)＝[μ₀μ(μ₀+μ)]^κ-1；

其中：

其中，θ₀,分別為出射太陽光天頂角和方位角；θ，分別為入射太陽光天頂角和方位角；μ₀，μ分別為出射太陽光天頂角余弦值和入射太陽光天頂角余弦值；ρ₀為地表反射太陽光照射的強(qiáng)度，κ為地表各向異性程度，α為控制太陽光照射相對前向或后向的散射量；M(κ)為與κ相關(guān)的函數(shù)，F(xiàn)(α)為與α相關(guān)的函數(shù)，H(ρ₀)為與ρ₀相關(guān)的函數(shù)。

3.根據(jù)權(quán)利要求2所述的三維蒙特卡洛輻射傳輸模式改進(jìn)方法，其特征在于，所述κ的取值范圍為：0.3≤κ≤1，當(dāng)κ＝1時(shí)，M(κ)函數(shù)值為1；

所述α的取值范圍為：當(dāng)α的取值范圍為：0≤α≤1時(shí)，α為控制太陽光照射相對前向的散射量；當(dāng)α的取值范圍為：-1≤α0時(shí)，α為控制太陽光照射相對后向的散射量。