大气辐射传输模型及其软件

大气辐射传输模型及其软件
大气辐射传输模型及其软件

大气辐射传输模型及其软件?

焦斌亮 高志强 李素静 白云

燕山大学信息科学与工程学院,河北 秦皇岛 066004

摘 要:本文主要阐述了大气辐射传输模型在大气订正中的应用,介绍了大气辐射传输原理,详细地叙述了6S 、LOWTRAN 、MODTRAN 和 FASCODE 等模型,同时提到了在以上模型基础上发展起来的其它辐射传输模型及软件,并对相应的模型及软件的共同特点和主要区别进行了比较,认为大气辐射传输模型在当前的大气订正模型中依然是比较可靠而常用的方法。

关键词:大气订正 辐射传输 6S MODTRAN

1 引 言

大气订正是遥感技术的重要组成部分,主要包括大气参数估计和地表反射率反演两个方面。如果获得了大气特性参数,进行大气订正就变得相对容易,但是获得准确的大气特性参数通常比较困难。通常有两类方法用辐射传输方程来计算大气订正函数:一种是直接的方法,对于大气透过率函数和反射率函数,通过对模型的积分来得到;另一种是间接的方法,它不是直接计算所需要的大气订正函数,而是通过辐射传输模型输出的表观反射率,结合模型输入的参数来求解。大气订正方法有很多,比如:基于图像特征的相对订正法、基于地面线形回归模型法、大气辐射传输模型法和复合模型法等。它是利用电磁波在大气中的辐射传输原理建立起来的模型对遥感图像进行大气订正的方法。 其中,大气辐射传输模型(Atmospheric Radiative Transfer Model)法是较常用的大气订正方法,它用于模拟大气与地表信息之间耦合作用的结果,其过程可以描述为地表光谱信息与大气耦合以后,在遥感器上所获得的信息,其中考虑了光子与大气相互作用机理,物理意义明确,具有很高的反演精度。

2 大气辐射传输原理

电磁辐射在介质中传输时,通常因其与物质的相互作用而减弱。辐射强度的减弱主要是由物质对辐射的吸收和物质散射所造成的,有时也会因相同波长上物质的发射以及多次散射而增强,多次散射使所有其它方向的一部分辐射进入所研究的辐射方向。当电磁辐射为太阳辐射,而且忽略多次散射产生的漫射辐射时,光谱辐射强度的变化规律可以表述为[1]

λλλρI ds

k dI ?= (1) 式中,I λ 是辐射强度,s 是辐射通过物质的厚度,ρ是物质密度,k λ表示对波长λ辐射的质量消光截面。

令在s=0处的入射强度为I λ(0)

,则在经过一定距离s 1后,其出射强度可由式(1)积分得到

?

作者介绍:焦斌亮(1964—),男(汉族),陕西户县人,燕山大学教授,主要从事光学遥感与CCD 应用技术研究;高志强(1981—),男(汉族),河北鹿泉人,燕山大学硕士研究生,研究方向:大气辐射在光学遥感中的应用。

(2) ()()??

?????=∫101exp 0s ds k I s I ρλλλ假定介质是均匀的,则k λ与距离s 无关,因此定义路径长度

(3) ∫=10s ds ρμ则式(2)可表示为

()()()u k I s I λλλ?=exp 01 (4) 上式就是比尔定律,也称朗伯定律。它指出,通过均匀消光介质传输的辐射强度按简单的指数函数减弱,该指数函数的自变量是质量消光截面和路径长度的乘积。它不仅适用于强度量,而且也适用于通量密度和通量。

根据式(4)我们可以定义单色透过率T λ为

()()(μλλλλk I s I T ?==

exp 01) (5) 式中,μ为θ的余弦值。

一般在大气辐射传输实际应用中,假定局域大气为平面平行的,因此只允许辐射强度和大气参数(温度和气体分布廓线)在垂直方向(即高度和气压)上变化,这种假定在物理意义上是适当的。如果用z 表示距离,则定义的普遍辐射传输方程可化为

()()(?θ?θρ)?θθ,;,;,;cos z J z I dz

k z dI +?= (6) 式中,θ为天顶角,φ为方位角,J 是源函数。

当考虑多次散射问题时,引进由大气上界向下测量的垂直光学厚度

(7)

∫∞=z dz k ρτ于是得到描述平面平行大气中多次散射问题的基本方程

()()(?μτ?μττ

)?μτμ,;,;,;J I d dI ?= (8) 3 典型的大气辐射传输模型

从20世纪80年代起,国外一些学者对遥感影像的大气订正研究做了许多工作,在模拟地—气过程的能力上有了很大提高,发展了一系列辐射传输模型,如6S 、LOWTRAN 、MODTRAN 和FASCODE 模型等,下面分别介绍。

3.1 6S 模型

6S (Second Simulation of the Satellite Signal in the Solar Spectrum )模型估计了0.25- 4.0μm 波长电磁波在晴空无云条件下的辐射特性,是在Tanre 等人[2]提出的5S (Simulation of the Satellite Signal in the Solar Spectrum)基础上发展而来的。它在假设均一地表的前提下,描述了非朗伯反射地表情况下的大气影响理论,而后Vermote [3]又将其改进为6S 模型。目前普遍使用的是1997年的4.1版本(以前是3.1版本)。当前若干

大气订正算法利用6S 模型来计算大气订正函数,一般是通过解近似的辐射传输方程来求得所需要的各种直射、散射透过率、大气的程辐射和大气的半球反照率等参数。其基本公式可以写为

1)()()[,(),,(s

s v s a r v s g v s v s t S T T T ρρθθρθθ??θθρ?+=?+ (9) 此处ρt 是传感器测得的表观反射率,ρs 是地表漫反射率,S 为大气半球反照率,ρr+a 为由分子散射和气溶胶散射所构成的大气路径辐射反射率,T g (θs ,θv )为大气吸收所构成的透过率,T(θs ) 与T(θv )分别代表太阳-目标与目标-传感器路径上的直线透过率,ρt 为大气上界的总反射率。

6S 模型主要包括以下5个部分:太阳、地物与传感器之间的几何关系,大气模式,气溶胶模式,传感器的光谱特性和地表反射率,它考虑了太阳的辐射能量通过大气传递到地表,再经地表反射通过大气传递到传感器的整个传播过程。对于吸收系数的计算公式,采用了吸收线的随机指数分布统计模式,这对于宽带传感器是一种很好的近似。为了考虑多次散射及分子散射与气溶胶散射及其相互作用,6S 采用最新近似 (state-of-the-art) 和连续散射SOS (Successive Order of Scattering) 方法来求解辐射传输方程。

3.2 LOWTRAN

LOWTRAN 是由美国地球物理实验室开发的单参数,谱带模式的大气传输模型,是计算大气透过率及辐射的软件包,其原意是“低谱分辨率大气透过率计算程序”[4],适用于从紫外、可见、红外到微波乃至更宽的电磁波谱范围内,包括云、雾、雨等多种大气状况的大气透过率及背景辐射。目前普遍使用的是1989 年2月公布的LOWTRAN7版本。它以20cm -1的光谱分辨率的单参数带模式计算0cm -1 -50 000cm -1的大气透过率、大气背景辐射、单次散射的阳光和月光辐射亮度、太阳直射辐照度。程序考虑了连续吸收,分子、气溶胶、云、雨的散射和吸收,地球曲率及折射对路径及总吸收物质含量计算的影响。大气模式包括13种微量气体的垂直廓线,六种参考大气模式定义了温度、气压、密度、及水汽、臭氧、甲烷、一氧化碳和一氧化二氮的混合比垂直廓线。程序用带模式计算水汽、臭氧、一氧化二氮、甲烷、一氧化碳、氧气、二氧化碳、一氧化氮和二氧化硫的透过率。多次散射参数化计算使用二流近似和累加法,用k-分布与带模式透过率计算衔接。对于水汽、氮气连续吸收以及紫外和可见波段的臭氧吸收,其平均透过率用朗伯定律计算,对其它气体吸收,则采用了双指数经验公式。

3.3 MODTRAN

MODTRAN 是LOWTRAN 的改进模型,其程序的基本结构和框架保持原样。它覆盖了0-22 600 cm -1 (即波长0.44μm- ∞)的光谱范围,具有2 cm -1的光谱分辨率。它利用二流(two steams)近似模型考虑大气多次散射效应。MODTRAN 是一个中分辨率大气辐射传输模型,吸收带模式参数用最新HITRAN [5],[6]数据库计算而得,采用Curtis-Godson 近似将多层的分层路径近似为等价的均匀路径,而且可以计算热红外的辐射亮度、辐照度等。目前普遍使用的是MODTRAN4版本。

3.4 FASCODE

FASCODE 是一个全世界公认的、以完全的逐线 Beer-Lambert 算法计算大气透过率和辐射的软件,它的分辨率很高,提供了“精确”透过率计算,并且考虑了非局地热力平衡状态的影响,原则上它的应用高度不受限制。因此,FASCODE 通常用作评估遥感系统或参数化带模型的标准,也常用于大气精细化结构的研究。

许多大气订正模型就是在以上模型的基础上发展起来的,如SMAC (Simplified method for the atmospheric correction )、ATREM (Atmospheric Removal)、HATCH (The high accuracy atmospheric correction for hyperspectral data)、ATCOR (Atmospheric and Topographic Correction mode1)、ACORN (Atmospheric

CORrection Now)和FLAASH(Fast line-of-sight atmospheric analysis of spectral hypercubes)模型,PcLnWin 和DISORT软件等。

此外,还有SHARC、UVRAD(Ultraviolet and Visible Radiation),TURNER,UCSB的SBDART、SAMM、SERTRAN模型,三维辐射传输模型(MOD3D)等大气辐射传输模型。

4 四种大气辐射传输模型的共同特点和主要差别

4.1 各种大气辐射传输模型的共同特点

6S、LOWTRAN,MODTRAN和 FASCODE是根据不同应用目的而开发的宽带、窄带和逐线计算的大气辐射传输模型及其相应的应用软件,都是用Fortran语言编写的。这几个大气辐射传输模型具有如下共同特点:它们以大气条件和地表条件作为输入参数,以表观反射率为输出结果,都涉及了复杂大气条件下多种辐射传输量的计算。

这些大气模型可以根据理论计算或实测资料,在这些传输模型实用程序中包括了具有代表性的大气和气溶胶的模式(如大气模式有6种:热带大气、中纬度夏季大气、中纬度冬季大气、亚北极区夏季大气、亚北极区冬季大气、美国标准大气),这些复杂的天气环境使它们具有更广泛的应用。而且,还可以由用户自定义模型大气,使模型在用户指定环境下模拟和使用显得特别灵活。此外,还包括了水平、垂直、倾斜向上和向下传输等各种复杂的几何关系,在计算大气倾斜路径及沿着传输路径衰减量时,都考虑了大气折射和地球的曲率。

提供大量的参数文件查找表(Look up table) ,查找表把全球气溶胶划分为若干类型,每种类型的大气参数由观测者获得。这些辐射传输模型都利用了 HITRAN数据库中的基本分子常数。

在计算地表反射率时,为了描述地表的均一性和朗伯状态,都引入了点扩散函数(PSF)与典型的双向反射率分布函数(BRDF)模型。

主要输入参数也相同,有几何条件、大气模式、气溶胶模式、地面能见度、目标高度、BRDF模型和计算波长位置等。

4.2 各种大气辐射传输模型的主要区别

MODTRAN把LOWTRAN 20 cm-1 (FWHM)的光谱分辨率和在 5 cm-1光谱间隔上做分子吸收计算改进为2 cm-1的光谱分辨率和在l cm-1光谱间隔上做分子吸收计算。在计算分子透过率方法上,LOWTRAN 采用单参数带模式,MODTRAN 采用的带模式使用3个与温度相关参数:吸收系数、线密度参数和平均线宽,使之更精确地服从分子跃迁的温度和压力关系(能级粒子数和V oigt线形),而且可以计算热红外的辐射亮度、辐照度等。

它们在不同的波谱范围内使用,LOWTRAN和MODTRAN可以在非常宽的波谱范围内使用,也可以对雨、云处理;而6S仅对太阳辐射中的紫外、可见光与近红外波段有效,对中远红外以外的波段无法处理,也没有考虑雨、云的情况;LOWTRAN和MODTRAN能够处理朗伯体和均匀表面,而6S还可以处理均一非朗伯地表和不均匀地表的情形。

这些辐射传输模型都利用了HITRAN数据库中的基本分子常数,然而,它们采取了各不相同的处置方法把这些常数换算成透过率和辐射度,使得这些辐射传输模型具有不同的精度与速度。从精度上说,FASCODE最好,6S 比MODTRAN 稍微差些,LOWTRAN最差。从速度上说,6S 比MODTRAN与LOWTRAN快的多。

从应用范围上说,6S仅在辐射信号传输中应用,而MODTRAN等在辐射能量传输和辐射信号传输中均可应用,但在国内遥感研究机构多数采用6S模型;FASCODE多用于激光大气传输中。

就软件运行平台而言,6S 和 LOWTRAN 的平台是DOS,MODTRAN的平台是WINDOWS,故此后者比前两者易于操作。

6S、LOWTRAN 和 MODTRAN 没有对非局地热力平衡状态的分子带进行合理的处理,而 FASCODE 考虑了非局地热力平衡。

5 结束语

综上所述,大气辐射传输模型的特点决定了它在大气订正应用中的重要地位,但是它也有自身的缺点。研究表明大气辐射传输模型的方法计算量较大,而且需要输入很多参数。由于大气参数的随机性和非均匀分布,使得很难确定这些参数。如果测得的参数不正确将直接影响计算的精度,从而产生较大的误差,因此进行实时地遥感大气订正是非常困难的。当地表为非均一面时,就必须考虑邻近效应,但目前各种模型对此效应的处理较为简单或不予考虑。因此,这种方法的广泛应用也受到了一定限制。尽管如此,大气辐射传输模型在当前的大气订正模型中依然是比较可靠而常用的方法。

参考文献

[1]K.N.LIOU. 大气辐射导论[M]. 气象出版社, 2004.

[2] D Tanre, C Deroo, P Duhaut et al. Description of a computer code to simulate the satellite signal in the solar

spectrum: the 5S code[J]. International Journal of Remote Sensing, 1990, 11: 659-668.

[3]Vermote E F, Tanre D, Deuze J L et a1. The Second Simulation of the Satellite Signal in the Solar Spectrum.

[4](6S), 6S User Guide Version 2. France: Laboratoire d’Optique Atmospherique, 1997, 7: 5.

[5]吴健等. 大气中的光传输理论[M]. 北京邮电大学出版社, 2005.

[6]Rothman L S, Gamache R R, Goldman A et a1. The HITRAN database[J]. ApplOpt, 1987, 26: 4058.

[7]Massie S T. Goldman A.. Absorption parameters of very dense molecular spectra for the HITRAN

compilation [J]. Quant Spctrosc Radiant Transfer, 1992, 48(6): 713-719.

Atmospheric Radiative Transfer Model and Its Software

JIAO Bin-liang GAO Zhi-qiang LI Su-jing BAI Yun Institute of Information Science and Engineering, Yanshan University, Qinhuangdao, HeBei 066004, China

Abstract: This article is mainly elaborated the atmospheric radiative transfer model applied in atmospheric correction, introduced the atmospheric radiative transfer theory, in detail narrated 6S, LOWTRAN, MODTRAN and FASCODE model and so on, simultaneously mentioned to develop other radiative transfer models and software on the above model foundation, and has carried on the comparison to the common characteristic and the main difference of the relevant model and software in certain aspects, we think the atmospheric radiative transfer model is still a reliable and popular method in atmospheric correction model at present.

Keywords: atmospheric correction; radiative transfer; 6S;MODTRAN

大气辐射传输模型

[转载]大气辐射传输模型 已有 968 次阅读2010-11-6 14:31|个人分类:未分类|系统分类:科普集锦|关键词:辐射传输 转自https://www.360docs.net/doc/8214323824.html,/s/blog_4b700c4c0100jgl7.html 相对辐射校正和绝对辐射校正 基于物理模型的绝对辐射校是利用一系列参数(例如,卫星过境时的地物反射率,大气的能见度,太阳天顶角和卫星传感器的标定参数等)将遥感图像进行校正的方法。仪器引起的误差畸变一般在数据生产过程中由生产单位根据传感器参数进行了校正。对于用户来所,绝对辐射校正的方法主要是辐射传输模型法,该方法校正精度较高,它是利用电磁波在大气中的辐射传输原理建立起来的模型对遥感图像进行大气校正的方法。由于有不同的不同的假设条件和适用的范围,因此产生很多可选择的大气较正模型,例如 6S模型、LOWTRAN模型、MODTRAN模型、ATCOR模型等。 基于统计模型的相对辐射校正,主要包括不变目标法、黑暗像元法与直方图匹配法等等。不变目标法假定图像上存在具有较稳定反射辐射特性的像元,并且可确定这些像元的地理意义,那么就称这些像元为不变目标,这些不变目标在不同时相的遥感图像上的反射率将存在一种线性关系。当确定了不变目标以及它们在不同时相遥感图像中反射率的这种线性关系,就可以对遥感图像进行大气校正。黑暗像元法的基本原理就是在假定待校正的遥感图像上存在黑暗像元区域、地表朗伯面反射、大气性质均一,忽略大气多次散射辐照作用和邻近像元漫反射作用的前提下,反射率很小的黑暗像元由于大气的影响,而使得这些像元的反射率相对增加,可以认为这部分增加的反射率是由于大气程辐射的影响产生的。利用黑暗像元值计算出程辐射,并代入适当的大气校正模型,获得相应的参数后,通过计算就得到了地物真实的反射率。直方图匹配法是指如果确定某个没有受到大气影响的区域和受到大气影响的区域的反射率是相同的,并且可以确定出不受影响的区域,就可以利用它的直方图对受影响地区的直方图进行匹配处理。此外,还有很多基于统计模型的方法,如有人提出利用小波变换的遥感图像相对辐射校正方法。该方法对源图像小波变换域的低频成分实施辐射变换,并保持高频成分不变,重构的图像具有保持高频信息的特性,因而能够较好地保留原图像中由于地物变化引起的辐射差异;也有人利用主成分分析法把遥感图像中有用的信息和大气影响噪音区分开来。 大气辐射传输模型6S 1986年,法国Université des Sciences et Technologies de Lille(里尔科技大学)大气光学实验室Tanré等人为了简化大气辐射传输方程,开发了太阳光谱波段卫星信号模拟程序5S(SIMULATION OF THE SATELLITE SIGNAL IN THE SOLAR SPECTRUM),用来模拟地气系统中太阳辐射的传输过程并计算卫星入瞳处辐射亮度。1997年,Eric Vemote对5S进行了改进,发展到6S(SECOND SIMULATION OF THE SATELLITE SIGNAL IN THE SOLAR SPECTRUM),6S吸收了最新的散射计算方法,使太阳光谱波段的散射计算精度比5S有所提高。 这种模式是在假定无云大气的情况下,考虑了水汽、CO2、O3和O2的吸收、分子和气溶胶的散射以及非均一地面和双向反射率的问题。6S是对5S的改进,光谱积分的步长从5nm 改进到2.5nm,同5S 相比,它可以模拟机载观测、设置目标高程、解释BRDF作用和临近效应,增加了两种吸收气体的计

高斯扩散模型.

大气污染扩散 第一节大气结构与气象 有效地防止大气污染的途径,除了采用除尘及废气净化装置等各种工程技术手段外,还需充分利用大气的湍流混合作用对污染物的扩散稀释能力,即大气的自净能力。污染物从污染源排放到大气中的扩散过程及其危害程度,主要决定于气象因素,此外还与污染物的特征和排放特性,以及排放区的地形地貌状况有关。下面简要介绍大气结构以及气象条件的一些基本概念。 一、大气的结构 气象学中的大气是指地球引力作用下包围地球的空气层,其最外层的界限难以确定。通常把自地面至1200 km左右范围内的空气层称做大气圈或大气层,而空气总质量的98.2%集中在距离地球表面30 km以下。超过1200 km的范围,由于空气极其稀薄,一般视为宇宙空间。 自然状态的大气由多种气体的混合物、水蒸气和悬浮微粒组成。其中,纯净干空气中的氧气、氮气和氩气三种主要成分的总和占空气体积的99.97%,它们之间的比例从地面直到90km高空基本不变,为大气的恒定的组分;二氧化碳由于燃料燃烧和动物的呼吸,陆地的含量比海上多,臭氧主要集中在55~60km高空,水蒸气含量在4%以下,在极地或沙漠区的体积分数接近于零,这些为大气的可变的组分;而来源于人类社会生产和火山爆发、森林火灾、海啸、地震等暂时性的灾害排放的煤烟、粉尘、氯化氢、硫化氢、硫氧化物、氮氧化物、碳氧化物为大气的不定的组分。 大气的结构是指垂直(即竖直)方向上大气的 密度、温度及其组成的分布状况。根据大气温度在 垂直方向上的分布规律,可将大气划分为四层:对 流层、平流层、中间层和暖层,如图5-1所示。 1. 对流层 对流层是大气圈最靠近地面的一层,集中了大 气质量的75%和几乎全部的水蒸气、微尘杂质。受 太阳辐射与大气环流的影响,对流层中空气的湍流 运动和垂直方向混合比较强烈,主要的天气现象云 雨风雪等都发生在这一层,有可能形成污染物易于 扩散的气象条件,也可能生成对环境产生有危害的 逆温气象条件。因此,该层对大气污染物的扩散、输送和转化影响最大。 大气对流层的厚度不恒定,随地球纬度增高而降低,且与季节的变化有关,赤道附近约

大气污染扩散模型

第一节大气污染物的扩散 一、湍流与湍流扩散理论 1. 湍流 低层大气中的风向是不断地变化,上下左右出现摆动;同时,风速也是时强时弱,形成迅速的阵风起伏。风的这种强度与方向随时间不规则的变化形成的空气运动称为大气湍流。湍流运动是由无数结构紧密的流体微团——湍涡组成,其特征量的时间与空间分布都具有随机性,但它们的统计平均值仍然遵循一定的规律。大气湍流的流动特征尺度一般取离地面的高度,比流体在管道内流动时要大得多,湍涡的大小及其发展基本不受空间的限制,因此在较小的平均风速下就能有很高的雷诺数,从而达到湍流状态。所以近地层的大气始终处于湍流状态,尤其在大气边界层内,气流受下垫面影响,湍流运动更为剧烈。大气湍流造成流场各部分强烈混合,能使局部的污染气体或微粒迅速扩散。烟团在大气的湍流混合作用下,由湍涡不断把烟气推向周围空气中,同时又将周围的空气卷入烟团,从而形成烟气的快速扩散稀释过程。 烟气在大气中的扩散特征取决于是否存在 湍流以及湍涡的尺度(直径),如图5-7所示。 图5-7(a)为无湍流时,烟团仅仅依靠分子 扩散使烟团长大,烟团的扩散速率非常缓慢, 其扩散速率比湍流扩散小5~6个数量级;图5 -7(b)为烟团在远小于其尺度的湍涡中扩散, 由于烟团边缘受到小湍涡的扰动,逐渐与周边 空气混合而缓慢膨胀,浓度逐渐降低,烟流几乎呈直线向下风运动;图5-7(c)为烟团在与其尺度接近的湍涡中扩散,在湍涡的切入卷出作用下烟团被迅速撕裂,大幅度变形,横截面快速膨胀,因而扩散较快,烟流呈小摆幅曲线向下风运动;图5-7(d)为烟团在远大于其尺度的湍涡中扩散,烟团受大湍涡的卷吸扰动影响较弱,其本身膨胀有限,烟团在大湍涡的夹带下作较大摆幅的蛇形曲线运动。实际上烟云的扩散过程通常不是仅由上述单一情况所完成,因为大气中同时并存的湍涡具有各种不同的尺度。 根据湍流的形成与发展趋势,大气湍流可分为机械湍流和热力湍流两种形式。机械湍流是因地面的摩擦力使风在垂直方向产生速度梯度,或者由于地面障碍物(如山丘、树木与建筑物等)导致风向与风速的突然改变而造成的。热力湍流主要是由于地表受热不均匀,或因大气温度层结不稳定,在垂直方向产生温度梯度而造成的。一般近地面的大气湍流总是机械湍流和热力湍流的共同作用,其发展、结构特征及强弱决定于风速的大小、地面障碍物形成的粗糙度和低层大气的温度层结状况。 2. 湍流扩散与正态分布的基本理论 气体污染物进入大气后,一面随大气整体飘移,同时由于湍流混合,使污染物从高浓度区向低浓度区扩散稀释,其扩散程度取决于大气湍流的强度。大气污染的形成及其危害程度在于有害物质的浓度及其持续时间,大气扩散理论就是用数理方法来模拟各种大气污染源在

激光大气传输仿真系统的设计

<> 全部作者: 白林董健业 第1作者单位: 北京邮电大学电子工程学院 论文摘要: 本文讨论了激光大气传输的理论模型,然后给出了仿真系统设计方案。理论模型主要基于布格定律以及散射理论,仿真系统主要分为3个模块即用于参数设定的界面,仿真计算,仿真绘图。 关键词: 激光大气传输仿真系统(浏览全文) 发表日期: 2008年03月17日 同行评议: 该文给出了1个考虑散射和温度场梯度等激光大气传输系统仿真设计的框架,对光束在大气中传输仿真有1定的参考价值,但没有具体的仿真实例,另外对大气随机扰动(如热场扰动)如何进行修改?建议补充1个具体实例加已说明验证。 综合评价: 修改稿: 注:同行评议是由特聘的同行专家给出的评审意见,综合评价是综合专家对论文各要素的评议得出的数值,以1至5颗星显示。 <> 1.绪论

1.1研究背景 网络被认为是互联网发展的第三阶段。网络的设计和实施能够带来切身实际的利益,城域网、企业网、局域网、家庭网和个人网络都是网络发展的体现。网络发明的初衷并不仅仅是表现在它的规模上,而是互联互通,资源共享,消除资源访问的壁垒,让生活更加方便、快捷、高效。随着网络技术的发展,网络在应用方面也体现出了很大的潜力,能够共享和调度成千上万的计算设备协同并发工作,能汇聚数百万计的信息资源加以归类、分析和发布,还可以让世界每一个角落的人们实时沟通交流。在现代高速发展的社会里,企业与企业之间的联系日益密切,大量的、复杂的信息交流显得由为重要。随着电子科技的高速发展,那些如何复杂大量的信息,通过网络技术帮助下,就可以轻而易举的从某一地方传送到另一地方,而且简单、快速、准确,给人们带来了很大的方便。而在现代企业中,网络技术在管理中的应用,已显得举足轻重。随着企业信息化进程的进一步深入和发展,计算机在企业中的应用越来越广泛,而企业对计算机的依赖越来越强。随着网络应用的日益丰富以及人们在日常生活中对网络依赖的日渐紧密,那么对于网络吞吐量,网络延时,网络链路的稳定性以及网络服务的多样性就会产生新的要求,同时也希望网络应用的花销能更加低廉,这样针对电信网络运营商所提供的服务将会产生巨大的挑战,本实时通信系统的成功应用将会给运营商们提供更加方便,快捷,稳定,并且低廉的网络运营成本,本实时通信系统帮助企业实现巨大的商业价值的同时也为用户带来的更加高效,快速,稳定并且廉价的网络服务资源。 1.2 选题理论 更多文章 https://www.360docs.net/doc/8214323824.html,/ mxdwk <>

Fluent辐射传热模型理论以及相关设置

Flue nt辐射传热模型理论以及相关设置 目录 1概述..................................... 2基础理论................................... 2.1专业术语解释: ........................... 2.2FLUENT畐射模型介绍:......................... 2.3辐射模型适用范围总结 ........................ 3Flue nt实际案例操作............................ 3.1Casel-测试external emissivity 使用DC模型计算-2D 模型....... 3.2Case2-测试in ter nal emissivity- ........................... 使用DO模型计算-2D 模型 3.3仿真结论 ..............................

1概述 在传热的仿真中,有时候会不可避免的涉及到辐射传热,而我们对Fluent 中辐射模型的了解甚少,很难得到可靠的计算结果。因此,一直以来,Fluent 中的带辐射的传热仿真是我们的一个难点,本专题重点来学习辐射模型的理论,让我们对辐射计算模型有一个深入的了解,以帮助我们攻克这个仿真难点。 2基础理论 2.1 专业术语解释: 在Fluent 中开启辐射模型时,流体介质以及固体壁面会出现一些专业的参数需要用户来设置。 在Fluent help 中介绍辐射模型时会经常提到一些专业术语。 对这些专业参数以及术语,我们来一一解释: 1、Optical thickness (光学深度,无量纲量):介质层不透明性的量度。即介质吸收辐射的能力的量度,等于入射辐射强度与出射辐射强度之比。设入射到吸收物质层的入射辐射强度为I ,透射的辐射强度为e,则T = l/e,其中T为光学深度。按照此定义,那介质完全透明,对辐射不吸收、也不散射,透射的辐射强度e= 入射辐射强度I,即光学深度为T=1,介质不参与辐射。一摘自百度百 科 而FLUENT中T=a L,其中L为介质的特征长度,a为辐射削弱系数(可理解为介质因吸收和散射引起的光强削弱系数)。如果T=0,说明介质不参与辐射,和百度百科中的定义有出入。但是所表达的意思是接近的,一个是前后辐射量的比值;一个是变化量和入射辐射量的比值(根据Fluent help 里的解释,经过介质的辐射损失量=I*T,个人理解,按照此定义,T不可能大于1啊,矛盾。// Theory Guide :: 0 // 5. Heat Transfer // 5.3. Modeling Radiation // 5.3.2. Radiative Transfer Equation )。该问题的解释为:

遥感辐射传输模型

遥感辐射传输模型 姓名:张超 学院:地球科学与环境工程学院 专业:遥感科学与技术 班级:遥感一班 提交时间:2015年5月10日 大气订正是遥感技术的重要组成部分,主要包括大气参数估计和地表反射率反演两个方面。如果获得了大气特性参数,进行大气订正就变得相对容易,但是

获得准确的大气特性参数通常比较困难。通常有两类方法用辐射传输方程计算大气订正函数:一种是直接的方法,对于大气透过率函数和反射率函数,通过对模型的积分来得到;另一种是间接的方法,他不是直接计算所需要的大气订正函数,而是通过辐射传输模型输出的表观反射率,结合模型输入的参数来求解。大气订正方法有很多,比如:基于图像特征的相对订正法、基于地面线形回归模型法、大气辐射传输模型法和复合模型法等。它是利用电磁波在大气中的 辐射传输原理建立起来的模型对遥感图像进行大气订正的方法。 其中,大气辐射传输模型(Atmospheric Radiative Transfer Model)法是较常用的大气订正方法,它用于模拟大气与地表信息之间耦合作用的结果,其过程可以描述为地表光谱信息与大气耦合以后,在遥感器上所获得的信息,其中考虑了光子与大气相互作用机理,物理意义明确,具有很高的反演精度。 大气辐射传输原理 电磁辐射在介质中传输时,通常因其与物质的相互作用而减弱。辐射强度的减弱主要是由物质对辐射的吸收和物质散射所造成的,有时也会因相同波长上物质的发射以及多次散射而增强,多次散射使所有其它方向的一部分辐射进入所研究的辐射方向。当电磁辐射为太阳辐射,而且忽略多次散射产生的漫射辐射时,光谱辐射强度的变化规律可以表述为[1] (1)式中,IΛ是辐射强度, s是辐射通过物质的厚度,ρ是物质密度,KΛ表示对波长λ辐射的质量消光截面。令在s=0 处的入射强度为Iλ(0),则在经过一定距离s1后,其出射强度可由式(1)积分得到 (2)假定介质是均匀的,则kλ与距离s无关,因此定义路径长度 (3)则式(2)可表示为 (4)上式就是比尔定律,也称朗伯定律。它指出,通过均匀消光介质传输的辐射强度按简单的指数函数减 弱,该指数函数的自变量是质量消光截面和路径长度的乘积。它不仅适用于强度

大气辐射传输理论 第一章..

大气辐射传输理论 引言 学科定义: 1、大气辐射学研究辐射能在地球-大气系统内传输和转换的规律及其应用,属大气物理学的一个分支。大气辐射学是天气学、气候学、动力气象学、应用气象学、大气化学和大气遥感等学科的理论基础之一。 2、地球-大气系统的辐射差额是天气变化和气候形成及其演变的基本因素,可以说辐射过程与动力过程的作用共同决定了地球的气候环境。 学习、研究的意义 辐射是地气系统与宇宙空间能量交换的唯一方式 数值天气预报中需要定量化考察大气辐射过程 辐射传输规律是大气遥感的理论基础 气候问题——辐射强迫 近年来人类活动造成的地球大气气候变迁成为大气科学研究热点,其原因也在于人类活动所排放的某些物质会改变地球大气中的辐射过程所致。 大气辐射学主要研究内容: 一、地-气系统辐射传输的基本物理过程和规律,包括 1、太阳的辐射(97%E在0.3~3μm波段内,λ m=0.5μm附近); 2、地-气系统辐射(绝大部分E在4~80μm波段内,λ m=10μm附近); 3、不同地表状态云、气溶胶、水汽、臭氧、二氧化碳等对辐射传输的影响。 二、大气辐射学还要研究辐射传输方程的求解。 辐射传输方程:是描述辐射传播通过介质时与介质发生相互作用(吸收、散射、发射等)而使辐射能按一定规律传输的方程,在地球大气条件下,求解非常复杂,只能在一些假定下求得解析解,因此辐射传输方程的求解,一直是大气辐射学研究的重要内容。 三、另外,对辐射与天气、气候关系的研究也是大气辐射学的重要内容,它是从地-气系统辐射收支的角度,来研究天气和气候的形成以及气候变迁问题的。 相关内容: 许多复杂的物理动力气候学问题中,涉及到海洋、极冰、陆地表面的辐射和热状况,大气中的云、气溶胶、二氧化碳等因子在辐射过程中对气候所造成的影响,以及这些过程和大气辐射过程之间复杂的相互作用和反馈关系。 第一章用于大气辐射的基本知识 第一节辐射的基本概念 太阳辐射和地球大气辐射虽具有不同的特性,其本质是相同的,它们都是电磁辐射。电磁辐射是以波动和粒子形式表现出的一种能量传送形式。 1.1.1电磁波及其特性 一、波:波是振动在空间的传播。有横波和纵波的形式之分。 二、机械波:机械振动在媒质中的传播,如声波、水波和地震波。 三、电磁波(ElectroMagnetic Spectrum):变化电场和变化磁场在空间的传播。 四、电磁辐射: 电磁能量的传递过程(包括辐射、吸收、反射和投射)称为电磁辐射。 五、电磁波的特性: 1、电磁波是横波 2、在真空中以光速传播 3、电磁波具有波粒二相性: 波动性:表现在电磁辐射以波动方式在大气中传播,并发生反射、折射、衍射和偏振等效应。也就是说电

污染物扩散模型-深圳数学建模

赛区评阅编号(由赛区组委会填写): 2015高教社杯全国大学生数学建模竞赛 承诺书 我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛参赛规则》(以下简称为“竞赛章程和参赛规则”,可从全国大学生数学建模竞赛网站下载)。 我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。 我们知道,抄袭别人的成果是违反竞赛章程和参赛规则的,如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。 我们郑重承诺,严格遵守竞赛章程和参赛规则,以保证竞赛的公正、公平性。如有违反竞赛章程和参赛规则的行为,我们将受到严肃处理。 我们授权全国大学生数学建模竞赛组委会,可将我们的论文以任何形式进行公开展示(包括进行网上公示,在书籍、期刊和其他媒体进行正式或非正式发表等)。 我们参赛选择的题号(从A/B/C/D中选择一项填写): C 我们的报名参赛队号(12位数字全国统一编号): 参赛学校(完整的学校全称,不含院系名):温州医科大学 参赛队员 (打印并签名) :1. 章成俊 2. 杨超 3. 谢锦 指导教师或指导教师组负责人 (打印并签名): 日期:年月日

赛区评阅编号(由赛区组委会填写): 2015高教社杯全国大学生数学建模竞赛 编号专用页 送全国评阅统一编号(由赛区组委会填写): 全国评阅随机编号(由全国组委会填写):

对垃圾处理厂污染的动态监控及居民补偿 摘要 城市垃圾处理问题是一个世界性难题。目前垃圾焚烧正逐步成为中国垃圾处理的主要手段之一。本论文构根据题目设置的垃圾处理厂规模,建立了环境动态监控体系,并根据潜在污染风险对周围居民进行了合理经济补偿的设计。 对于问题(1),为了实现对垃圾焚烧厂烟气排放及相关环境影响状况的动态监控,本论文在高斯烟羽模型的基础上进行改进,引入温度、降雨对污染物扩散的影响,建立了新的污染物扩散模型。本论文创新性的提出了风雨影响指数M,用来衡量风向、降雨对颗粒物扩散的影响。本论文将抽象的污染物含量形象化,利用空气污染指数API描述具体的污染程度及其给周围居民带来的影响。并且从不同角度给出了模型检验,验证了所建模型的准确性。 对于问题(1)具体赔偿方案的制定,在综合考虑了不同方位风向频率、受污染时间、受污染程度的基础上,本论文使用了层次分析法,并且进行了一致性检验,使得赔偿方案具有说服力。通过MATLAB编程,计算出当政府和垃圾处理厂共支付风险赔偿金为N时,得出居住地的每位居民应得的赔偿金额计算公式。对于监测点的设置,经计算共需21个,具体布置情况见后文。 对于问题(2),在题目所述的发生事故的情况下,对污染物的具体含量进行了合理的预测与假设。模拟出酸性物质与颗粒物的影响范围,并根据具体的污染程度设置不同的污染区。对每个污染区的不同情况设置更改监测点的设置,并且在问题(1)的基础上对居民的经济补偿进行合理修改。 关键词:高斯烟羽模型,层次分析法,空气污染指数,烟气抬升公式 一、问题重述 “垃圾围城”是世界性难题,在今天的中国显得尤为突出。数据显示,目前全国三分之二以上的城市面临“垃圾围城”问题,垃圾堆放累计侵占土地75万亩。因此,垃圾焚烧正逐步成为中国垃圾处理的主要手段之一。然而,由于政府监管不力、投资者目光短浅等多方面的原因,致使前些年各地建设的垃圾焚烧电厂在运营中出现了环境污染问题,给垃圾焚烧技术在我国的推广造成了很大阻力,许多城市的新建垃圾焚烧厂选址都出现因居民反对而难以落地的局面。在垃圾焚烧厂运行监管方面,目前主要是在垃圾焚烧厂内进行测量监控,缺少从周边环境视角出发的外围动态监控,因而难以形成为民众所信服的全方位垃圾焚烧厂环境监控体系。 深圳市某地点计划建立一个中型的垃圾焚烧厂,计划处理垃圾量1950吨/天(设置三台可处理垃圾650吨/天的焚烧炉,排烟口高度80米,每天24小时运转)。从构建环境动态监控体系、并根据潜在污染风险对周围居民进行合理经济补偿的需求出发,有关部门希望能综合考虑垃圾焚烧厂对周围带来环境污染以及其他危害的多种因素(例如,焚烧炉的污染物排放量、居住点离开垃圾焚烧厂的距离、风力和风向及降雨等气象条件、地形地貌以及建筑物的遮挡程度等等),在进行科学定量分析的基础

大气污染物扩散高斯模型模拟

大气污染物扩散的高斯模型模拟:可视化模拟点源大气污染的扩散Gaussian Atmospheric Dispersion Model 突发性大气污染事故时有发生,对大气污染扩散进行模拟和分析,有利于减小事故的危害,减轻人员伤亡和财产损失。高斯扩散模型是国际原子能机构(IAEA)推荐使用于重气云扩散模拟的数学模型,该模型在非重气云扩散的应用日益广泛。高斯扩散模型是描述大气对有害气体的输移、扩散和稀释作用的物理或数学模型,是进行灾害预测和救援指挥的有力手段之一。 高斯扩散模型 高斯模型又分为高斯烟团模型和高斯烟羽模型。大气污染物泄漏分为瞬时泄漏和连续泄漏,瞬时泄漏是指污染物泄放的时间相对于污染物扩散的时间较短如突发泄漏等的情形,连续泄漏则是指污染物泄放的时间较长的情形。瞬时泄漏采用高斯烟团模型模拟,而连续泄漏采用高斯模型烟羽模型模拟。高斯模型适用于非重气云气体,包括轻气云和中性气云气体。要求气体在扩散过程中,风速均匀稳定。 在高斯烟团模型中,选择风向建立坐标系统,即取泄漏源为坐标原点,x轴指向风向,y轴表示在水平面内与风向垂直的方向,z轴则指向与水平面垂直的方向,具体公式见式: (mg/s); x、y、z轴上的扩散系数,需根据大气稳定度选择参数计算得到(m);x、y、z表示x、y、z上的坐标值(m);u 表示平均风速(m/s);t表示扩散时间(s);H 表示泄漏源的高度(m)。 同理,高斯烟羽模型的表达式如: 技术方法 若用高斯模型算出空间每一个点在一个时刻的污染浓度,这个计算量是很大的。因此所设计的系统一般都是采用先进行图层网格化,由高斯模型计算出有限个网格点的上的污染物浓度,在进行空间内插得到面上每一个点的污染物浓度,并由此得到污染物浓度的等值线。整个过程的示意图如图所示

大气辐射传输模型6S简介

大气辐射传输模型6S简介 1986年,法国Université des Sciences et Technologies de Lille(里尔科技大学)大气光学实验室Tanré等人为了简化大气辐射传输方程,开发了太阳光谱波段卫星信号模拟程序5S(SIMULATION OF THE SATELLITE SIGNAL IN THE SOLAR SPECTRUM),用来模拟地气系统中太阳辐射的传输过程并计算卫星入瞳处辐射亮度。1997年,Eric Vemote对5S进行了改进,发展到6S(SECOND SIMULATION OF THE SATELLITE SIGNAL IN THE SOLAR SPECTRUM),6S吸收了最新的散射计算方法,使太阳光谱波段的散射计算精度比5S有所提高。 这种模式是在假定无云大气的情况下,考虑了水汽、CO2、O3和O2的吸收、分子和气溶胶的散射以及非均一地面和双向反射率的问题。6S是对5S的改进,光谱积分的步长从5nm改进到2.5nm,同5S相比,它可以模拟机载观测、设置目标高程、解释BRDF作用和临近效应,增加了两种吸收气体的计算(CO、N2O)。采用SOS (successive order of scattering) 方法计算散射作用以提高精度。缺点是不能处理球形大气和limb (临边)观测。 它其中主要包括以下几个部分: (1)太阳、地物与传感器之间的几何关系:用太阳天顶角、太阳方位角、观测天顶角、观测方位角四个变量来描述; (2)大气模式:定义了大气的基本成分以及温湿度廓线,包括7种模式,还可以通过自定义的方式来输入由实测的探空数据,生成局地更为精确、实时的大气模式,此外,还可以改变水汽和臭氧含量的模式; (3)气溶胶模式:定义了全球主要的气溶胶参数,如气溶胶相函数、非对称因子和单次散射反照率等,6S中定义了7种缺省的标准气溶胶模式和一些自定义模式; (4)传感器的光谱特性:定义了传感器的通道的光谱响应函数,6S中自带了大部分主要传感器的可见光近红外波段的通道相应光谱响应函数,如TM,MSS,POLDER和MODIS等; (5)地表反射率:定义了地表的反射率模型,包括均一地表与非均一地表两种情况,在均一地表中又考虑了有无方向性反射问题,在考虑方向性时用了9种不同模型)。 这5个部分便构成了辐射传输模型,考虑了大气顶的太阳辐射能量通过大气传递到地表,以及地表的反射辐射通过大气到达传感器的整个辐射传输过程。 6S的输入参数主要有9个部分组成:

污染空气的扩散模型

放射性气体扩散的预估模型 摘要:由于放射性气体泄漏造成惨重损失的报道在国际屡见不鲜,近日日本福岛核电站的放射性气体的泄漏事件更让我们关注放射性气体泄漏时在环境中的浓度问题,为了今后事故发生后提供积极的补救措施, 所以对放射性气体的扩散作深入的研究是很有必要的。本文结合高斯烟羽模型、线性拟合,以及微分方程模型,运用MA TLAB软件,分析了泄漏源强度、风速、大气稳定度参数、地面粗糙度参数和计算精确度等的因素对放射性气体扩散的影响,预测了放射性气体浓度在不同时间,不同地区的浓度变化,并且本文模型中的数据可以根据不同的实际情况而加以改变,因而使本文的应用范围大大增加,可以适用于具有较强的应用性。文章首先在第一问中利用MA TLAB软件对数据进行线性拟合,采用微分方程模型得到核电站周边放射性气体在不同地区,不同时间段的浓度变化,得出随着离泄漏源距离的延伸,最终放射性物质的浓度越来越小,趋近于零,即当L趋向无穷是,C(x,y,z,t)趋向于零;当时间趋于无穷时,C(x,y,z,t)也趋于无穷。问题二,问题三中,建立以核电站周边不同地区得距离以及风速为因变量,设置各个主要因素的参考数据,同时,利用高斯烟羽模型对核电站周边地区的浓度进行预测,然后,利用MATLAB软件,将相关数据代入程序,我们得到核电站周边地区的浓度分布的等高曲线。问题四中,通过实际收集数据,集合核电站周边地区的浓度等高曲线,可以直观的看出日本福岛核电站对我国东海岸以及美国西海岸的影响。 一.问题的提出 1.1背景的介绍 目前,核电的发展给国家带来了巨大的经济效益和社会效益,但核电正常运行以及发生泄露时不可避免的会有气载放射性核素排出,这样就给周围的环境产生了一定的影响,因此,正确的测出大气中放射性物质的浓度在环境检测以及安全评估中具有重要意义。 1.2需要解决的问题 的放射性气体以匀速排出,设有一座核电站遇自然灾害发生泄漏,浓度为p 速度为m kg/s,在无风的情况下,匀速在大气中向四周扩散, 速度为s m/s. (1)请你建立一个描述核电站周边不同距离地区、不同时段放射性物质浓度的预测模型。 (2)当风速为k m/s时,给出核电站周边放射性物质浓度的变化情况。 (3)当风速为k m/s时,分别给出上风和下风L公里处,放射性物质浓度的预测模型。

大气辐射传输校正模型(5S,modtran,acorn)

在遥感的实际应用中,常用很多简化的手段,如假设地面为朗伯面,排除云的存在,采用有关标准大气模式及大气气溶胶模式等,一次产生了许多不同类型的大气辐射传输模型,主要分为两类, 1)采用大气的光学参数 2)直接采用大气物理参数如lowtran、modtran等大气辐射近似计算模型,而且还增加了多次散射计算 1. 5s模型 该模型的代码模拟计算海平面上的均匀朗伯体目标的反射率,并假定大气吸收作用与散射作用可以耦合,就像吸收粒子位于散射层的上面一样,则大气上层测 量的目标反射率可以表示为, 海平面处朗伯体的反射率 大气透过率 分子、气溶胶层的内在反射率 有太阳到地表再到传感器的大气透过率 S为大气的反射率 大气传输辐射校正模型-3 modtran 该模型是由美国空军地球物理实验室研制的大气辐射模拟计算程序,在遥感领域被广泛应用于图像的大气校正。

lowtran7是一个光谱分辨率20cm-1,的大气辐射传输实用软件,它提供了6种参考大气模式的温度、气压、密度的垂直廓线,水汽、臭氧、甲烷、一氧化碳、一氧化二氮的混合比垂直廓线,其他13种微量气体的垂直廓线,城乡大气气溶胶、雾、沙尘、火山喷发物、云、雨的廓线,辐射参量(如消光系数、吸收系数、非对称因子的光谱分布),以及地外太阳光谱。 lowtran7可以根据用户的需要,设置水平、倾斜、及垂直路径,地对空、空对地等各种探测几何形式,适用对象广泛。lowtran7的基本算法包括透过率计算方法,多次散射处理和几何路径计算。 1)多次散射处理 lowtran 采用改进的累加法,自海平面开始向上直至大气的上界,全面考虑整层大气和地表、云层的反射贡献,逐层确定大气分层每一界面上的综合透过率、吸收率、反射率和辐射通量。再用得到的通量计算散射源函数,用二流近似解求辐射传输方程。 2)透过率计算 该模型在单纯计算透过率或仅考虑单次散射时,采用参数化经验方法计算带平均透过率,在计算多次散射时,采用k-分布法 3)光线几何路径计算 考虑了地球曲率和大气折射效应,将大气看作球面分层,逐层考虑大气折射效应 由于lowtran直接使用大气物理参数,因而需要按照下列方法计算出与 lowtran使用的大气物理参数相对应的大气光学参数179页 4.modtran辐射传输模型 modtran可以计算0到50000cm-1的大气透过率和辐射亮度,它在440nm到无限大的波长范围精度是2cm-1,在22680到50000cm-1紫外波(200-440nm)范围的精度是20cm-1,在给定辐射传输驱动、气溶胶和云参数、光源与遥感器的几何立体对和地面光谱信息的基础上,根据辐射传输方程来计算大气的透过率以及辐射亮度。

数学建模高斯扩散模型培训资料

数学建模高斯扩散模 型

§4-2高斯扩散模式 ū —平均风速; Q—源强是指污染物排放速率。与空气中污染物质的浓度成正比,它是研究空气污染问题的基础数据。通常: (ⅰ)瞬时点源的源强以一次释放的总量表示; (ⅱ)连续点源以单位时间的释放量表示; (ⅲ)连续线源以单位时间单位长度的排放量表示; (ⅳ)连续面源以单位时间单位面积的排放量表示。 δy—侧向扩散参数,污染物在y方向分布的标准偏差,是距离y的函数,m; δz—竖向扩散参数,污染物在z方向分布的标准偏差,是距离z的函数,m; 未知量—浓度c、待定函数A(x)、待定系数a、b; 式①、②、③、④组成一方程组,四个方程式有四个未知数,故方程式可解。 二、高斯扩散模式 (一)连续点源的扩散 连续点源一般指排放大量污染物的烟囱、放散管、通风口等。排放口安置在地面的称为地面点源,处于高空位置的称为高架点源。 1. 大空间点源扩散 高斯扩散公式的建立有如下假设:①风的平均流场稳定,风速均匀,风向平直;②污染物的浓度在y、z轴方向符合正态分布;③污染物在输送扩散中质量守恒; ④污染源的源强均匀、连续。 图5-9所示为点源的高斯扩散模式示意图。有效源位于坐标原点o处,平均风向与x轴平行,并与x轴正向同向。假设点源在没有任何障碍物的自由空间扩散,不考虑下垫面的存在。大气中的扩散是具有y与z两个坐标方向的二维正态分布,当两坐

标方向的随机变量独立时,分布密度为每个坐标方向的一维正态分布密度函数的乘积。由正态分布的假设条件②,参照正态分布函数的基本形式式(5-15),取μ=0,则在点源下风向任一点的浓度分布函数为: (5-16)式中 C—空间点(x,y,z)的污染物的浓度,mg/m3; A(x)—待定函数; σy、σz—分别为水平、垂直方向的标准差,即y、x方向的扩散参数,m。 由守恒和连续假设条件③和④,在任一垂直于x轴的烟流截面上有: (5-17) 式中 q—源强,即单位时间内排放的污染物,μg/s; u—平均风速,m/s。 将式(5-16)代入式(5-17), 由风速稳定假设条件①,A与y、z无关,考虑到③和④,积分可得待定函数A(x): (5-18) 将式(5-18)代入式(5-16),得大空间连续点源的高斯扩散模式 (5-19) 式中,扩散系数σy、σz与大气稳定度和水平距离x有关,并随x的增大而增加。当y=0,z=0时,A(x)=C(x,0,0),即A(x)为x轴上的浓度,也是垂直于x轴截面上污染物的最大浓度点C max。当x→∞,σy及σz→∞,则C→0,表明污染物以在大气中得以完全扩散。 2.高架点源扩散

Fluent辐射传热模型理论以及相关设置

Fluent辐射传热模型理论以及相关设置 目录 1概述 (2) 2基础理论 (2) 2.1专业术语解释: (2) 2.2FLUENT辐射模型介绍: (3) 2.3辐射模型适用范围总结 (4) 3Fluent实际案例操作 (5) 3.1Case1-测试external emissivity 使用DO模型计算-2D模型 (5) 3.2Case2-测试internal emissivity-使用DO模型计算-2D模型 (6) 3.3仿真结论 (10)

1概述 在传热的仿真中,有时候会不可避免的涉及到辐射传热,而我们对Fluent中辐射模型的了解甚少,很难得到可靠的计算结果。因此,一直以来,Fluent中的带辐射的传热仿真是我们的一个难点,本专题重点来学习辐射模型的理论,让我们对辐射计算模型有一个深入的了解,以帮助我们攻克这个仿真难点。2基础理论 2.1专业术语解释: 在Fluent中开启辐射模型时,流体介质以及固体壁面会出现一些专业的参数需要用户来设置。 在Fluent help中介绍辐射模型时会经常提到一些专业术语。 对这些专业参数以及术语,我们来一一解释: 1、Optical thickness(光学深度,无量纲量):介质层不透明性的量度。即介质吸收辐射的能力的量度,等于入射辐射强度与出射辐射强度之比。设入射到吸收物质层的入射辐射强度为I ,透射的辐射强度为e,则T = I/e,其中T为光学深度。按照此定义,那介质完全透明,对辐射不吸收、也不散射,透射的辐射强度e=入射辐射强度I,即光学深度为T=1,介质不参与辐射。—摘自百度百科而FLUENT中T=αL,其中L为介质的特征长度,α为辐射削弱系数(可理解为介质因吸收和散射引起的光强削弱系数)。如果T=0,说明介质不参与辐射,和百度百科中的定义有出入。但是所表达的意思是接近的,一个是前后辐射量的比值;一个是变化量和入射辐射量的比值(根据Fluent help里的解释,经过介质的辐射损失量=I*T,个人理解,按照此定义,T不可能大于1啊,矛盾。// Theory Guide :: 0 // 5. Heat Transfer // 5.3. Modeling Radiation // 5.3.2. Radiative Transfer Equation)。该问题的解释为:其实一点也不矛盾,如果Optical thickness =1,就说明辐射在经过一定特征长度L的介质后被完全吸收。如果>1,就说明辐射根本穿透不了特征长度L的介质,而被早早吸收完了。打个比方,Optical thickness=10,说明辐射在经过L/10距离后已经被吸收(或散射)完。 其中α=αA+αS; 2、Absorption Coefficient(αA吸收系数,单位1/m,见图2-1):因为介质吸收而导致的辐射强度在经过每单位长度介质后改变的量。空气作为流体介质时,一般不吸收热辐射,该系数可近视设为0。而当气体中水蒸气和CO2含量较高时,那对辐射的系数就不能忽略了。 3、Scattering Coefficient(αS散射系数,单位1/m):因为介质散射而导致的辐射强度在经过每单位长度介质后改变的量。空气作为流体介质时,一般情况下,该系数可近视设为0。对于含颗粒物的流体,散射作用不容忽视。

大气辐射传输模型及其软件

大气辐射传输模型及其软件? 焦斌亮 高志强 李素静 白云 燕山大学信息科学与工程学院,河北 秦皇岛 066004 摘 要:本文主要阐述了大气辐射传输模型在大气订正中的应用,介绍了大气辐射传输原理,详细地叙述了6S 、LOWTRAN 、MODTRAN 和 FASCODE 等模型,同时提到了在以上模型基础上发展起来的其它辐射传输模型及软件,并对相应的模型及软件的共同特点和主要区别进行了比较,认为大气辐射传输模型在当前的大气订正模型中依然是比较可靠而常用的方法。 关键词:大气订正 辐射传输 6S MODTRAN 1 引 言 大气订正是遥感技术的重要组成部分,主要包括大气参数估计和地表反射率反演两个方面。如果获得了大气特性参数,进行大气订正就变得相对容易,但是获得准确的大气特性参数通常比较困难。通常有两类方法用辐射传输方程来计算大气订正函数:一种是直接的方法,对于大气透过率函数和反射率函数,通过对模型的积分来得到;另一种是间接的方法,它不是直接计算所需要的大气订正函数,而是通过辐射传输模型输出的表观反射率,结合模型输入的参数来求解。大气订正方法有很多,比如:基于图像特征的相对订正法、基于地面线形回归模型法、大气辐射传输模型法和复合模型法等。它是利用电磁波在大气中的辐射传输原理建立起来的模型对遥感图像进行大气订正的方法。 其中,大气辐射传输模型(Atmospheric Radiative Transfer Model)法是较常用的大气订正方法,它用于模拟大气与地表信息之间耦合作用的结果,其过程可以描述为地表光谱信息与大气耦合以后,在遥感器上所获得的信息,其中考虑了光子与大气相互作用机理,物理意义明确,具有很高的反演精度。 2 大气辐射传输原理 电磁辐射在介质中传输时,通常因其与物质的相互作用而减弱。辐射强度的减弱主要是由物质对辐射的吸收和物质散射所造成的,有时也会因相同波长上物质的发射以及多次散射而增强,多次散射使所有其它方向的一部分辐射进入所研究的辐射方向。当电磁辐射为太阳辐射,而且忽略多次散射产生的漫射辐射时,光谱辐射强度的变化规律可以表述为[1] λλλρI ds k dI ?= (1) 式中,I λ 是辐射强度,s 是辐射通过物质的厚度,ρ是物质密度,k λ表示对波长λ辐射的质量消光截面。 令在s=0处的入射强度为I λ(0) ,则在经过一定距离s 1后,其出射强度可由式(1)积分得到 ? 作者介绍:焦斌亮(1964—),男(汉族),陕西户县人,燕山大学教授,主要从事光学遥感与CCD 应用技术研究;高志强(1981—),男(汉族),河北鹿泉人,燕山大学硕士研究生,研究方向:大气辐射在光学遥感中的应用。

点污染源空气污染扩散模型

8点、中午12点、晚上9点都没有排放气体,该怎么算,是不是需要找到一个关于时间t的函数,来计算多长时间之后污染还剩下多少 c=Q./(2*pi*sigy.*sigz*u+eps)*exp(-0.5*(y.A2)./((sigy+eps).A2)).*(exp(-0.5*(z-H).A2./((sigz+eps).A2))+exp(-0.5*(z+H).A2./((sigz+ eps)A2))); 这个函数对吗?该调用什么函数? 问题: 建立单污染源空气污染扩散模型,描述其对周围空气污染的动态影响规律。 现有河北境内某一工厂废气排放烟囱高50m,主要排放物为氮氧化物。早上9点至下午 3点期间的排放浓度为406.92mg/m3,排放速度为1200m3 /h;晚上10点-凌晨4点期间 的排放浓度为1160mg/m3,排放速度为5700m3 /h;通过你的扩散模型求解该工厂方圆51 公里分别在早上浓度8点、中午12点、晚上9点空气污染分布和空气质量等级。 源代码 clear all clc [x,y]=meshgrid(0:20:5100,0:20:5100); Q=135.64; z=1.5; H=50; u=1.94; sigy=0.3914238*x.A0.865014; sigz=0.0757182*x.A1.00770; c=Q./(2*pi*sigy.*sigz*u+eps).*exp(-0.5*(y.A2)./((sigy+eps).A2)).*(exp(-0.5*(z-H).A2./((sigz+eps).A2))+exp(-0.5*(z+H).A2./((sigz+ eps)A2))); mesh(x,y,c); xlabel('X'),ylabel('Y'),zlabel(C), clear all clc [x,y]=meshgrid(-5100:20:5100,-5100:20:5100); Q=1836.7; z=1.5; H=50; u=1.7; sigy=0.3914238*x.A0.865014; sigz=0.0757182*x.A1.00770; c=Q./(2*pi*sigy.*sigz*u+eps).*exp(-0.5*(y.A2)./((sigy+eps).A2)).*(exp(-0.5*(z-H).A2./((sigz+eps).A2))+exp(-0.5*(z+H).A2./((sigz+ eps)A2))); mesh(x,y,c); xlabel('X'),ylabel('Y'),zlabel(C), 分享到: 2015-05-29 16:32 提问者采纳 clear all [x,y]=meshgrid(-51000:100:51000,-51000:100:51000); Q=135.64; z=1.5; H=50; u=1.94; sigy=0.3914238*x.A0.865014;

相关文档
最新文档