大气辐射传输模型
[转载]大气辐射传输模型
已有 968 次阅读2010-11-6 14:31|个人分类:未分类|系统分类:科普集锦|关键词:辐射传输
转自https://www.360docs.net/doc/0c3499640.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作用和临近效应，增加了两种吸收气体的计
算（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个部分组成：
（1）几何参数
6S两种输入方法⑴太阳和卫星的天顶角和方位角以及观测时间(月,日)。⑵卫星的接收时间(月,日,年)、像素点数、升交点时间，由程序计算太阳和卫星的天顶角和方位角。特别注意的是这里的时间采用世界时且要精确到1/6秒。
（2）大气模式
6S给出几种可供选择的大气模式，热带、中纬度夏季、中纬度冬季、近极地夏季、近极地冬季、美国62标准大气也可自定义大气模式。
（3）气溶胶模式
三种选择：⑴无气溶胶。⑵自定义气溶胶模式。如，四种基本气溶胶的体积的加权平均；气溶胶的谱分布加光度计测量结果（光学厚度）和复折射指数；直接给出消光系数。⑶提供的三种气溶胶模式大陆型，海洋型和乡村型。
（4）气溶胶浓度
两种选择：⑴在550nm处的光学厚度⑵气象能见度（km）。故它也提供了两者的相互关系。
（5）地面高度
以千米为单位的地面海拔高度（设为负值）。
（6）探测器高度
-1000代表卫星测量，0为地基观测，飞机航测输入以千米为单位的负值。
（7）探测器的光谱条件
给出了常见卫星Meteosat，Goes，NOAA/AVHRR和HRV，Landsat TM 和MSS，Modis Polder的每个通道的光谱响应函数，也可选择自定义。
（8）地表特性
可以选择地表均一或不均一，也可选择地表为郎伯体或双向反射。6S给出了九种比较成熟的BRDF模式供用户选择，也可自定义BRDF函数（输入个角度的反射率及入射强度）
（9）表观反射率
输入反射率或辐射亮度，同时也决定模式是正向还是反向工作。当RAPP<-1时是正向。RAPP>0（辐射亮度）或-1 “湖泊水色遥感大气校正研究”是中科院南京地理与湖泊研究所委托南京大学GIS与遥感实验室的项目，受到中科院领域前沿项目（CXNIGLAS-A02-014）的资助。起止年限：2005.6—2007.6。 2. 研究内容 本研究通过卫星遥感影像获取湖泊水体的离水辐亮度或遥感反射率，即研究如何对湖泊水色遥感中的卫星影象进行大气校正。具体内容包括： 对现有的一般的大气校正方法/模式（如6S、LOWTRAN、MODTRAN等）与海洋水色遥感的大气校正方法/模式进行比较，并分析现有现有软件（如ENVI、PCI等）中的大气校正方法； 提出最适合湖泊水体的大气校正方法/流程/模式，用IDL或MATLAB编程实现。 研究的重点是内陆湖泊水体的遥感大气校正，其目的是找到一种适合内陆湖泊水体遥感大气校正的方法和流程。水体的遥感的目标是通过遥感手段来了解水体的物化参数，就目前而言，能从遥感波谱反射特征所获得的参数主要有叶绿素、悬浮物和黄色物质等。因此，研究这三种物质的光谱特征是水体遥感大气校正的前提。 理想的遥感模型是没有大气存在，地面为朗伯体，遥感接收到的光谱直接反映地面目标的状况。但实际情况中，太阳光在下行到达地面和上行到达传感器的过程中，受到大气的影响而衰减，大气校正就是要研究并消除大气条件对遥感的影响，在传感器位置恢复地面目标的光谱特征。 研究区选在太湖，太湖地处长江三角洲，面积36900平方公里，为流域第一大湖，又是长江中下游五大淡水湖之一。水面面积为2338.1km2。 3. 技术路线 水色遥感是快速获取大面积水体组分参数的重要手段。太阳辐射在太阳-水体目标-传感器的传输过程中，气体的吸收、气溶胶与分子的散射这两种大气过程影响了获得的水体反射率信息。由于水体在可见光和近红外波段的反射率很低，传感器接收到的信号只有5-15%来自水体。因此，利用遥感影像精确反演水体信息，大气校正是至关重要的。 大气校正就是在遥感图像的基础上，消除大气影响，获得地面反射率的过程。许多学者研究大气对辐射传输的影响，提出了不同的大气校正模型，主要包括6S，LOWTRAN/MODT- RAN ，TURNER等。本文主要研究前两种模型。 6S （Second Simulation of the Satellite Signal in the Solar Spectrum）大气校正模型是Eric F. Vermote et al.（1997）在5S模型的基础上发展起来的。6S模型可以很好地模拟太阳光在太阳－地面目标－传感器的传输过程中所受到的大气影响。相对于5S模型，6S模型考虑了地面目标的海拔高度、非朗伯平面的情况和新的吸收气体种类（CH4，N2O，CO）。通过采用the art approximation 近似算法和SOS运算法则，提高了瑞利和气溶胶散射作用的计算精度。光谱步长提高到了2.5nm。6S 模型建立在辐射传输理论基础之上,模型应用范围广，不受研究区特点及目标类型等的影响。 LOWTRAN是一种低分辨率（分辨率≥20cm-1）大气辐射传输模式。在20多年的发展过程中不断扩充和修订基础资料，改进算法，增加可计算的辐射传输结果，从原意义上的“低分辨率大气透过率计算模式”扩展到目前能导出复杂天气条件下多种辐射传输量的“低分辨率大气辐射传输计算模式”，提供了许多新的应用可能性，已被国际上许多应用专家广泛应用于各自的实际问题。MODTRAN（中光谱分辨率大气辐射传输模式）较之LOWTRAN不但提高了光谱分辨率，而且还包括了多次散射辐射传输精确算法——离散纵标法，对有散射大气的辐射传输如太阳短波辐射，比LOWTRAN中的二流近似算法有更高的精度和更大的灵活性。 由于基于辐射传输理论的大气校正需要相应的大气参数，而在一般研究中，这些参数一般很难获得，针对不同研究区的的大气校正方法有时是很实用的。本研究中也尝试使用黑暗象元法、不变目标法、直方图匹配法、大气阻抗植被指数法、主成分分析法等获取水体的遥感反射率或直接获取水体参数信息，并与基本基于辐射传输理论的大气校正结果进行比较。 本文引用地址：https://www.360docs.net/doc/0c3499640.html,/home.php?mod=space&uid=499915&do=blog&id=380974 [转载]大气辐射传输模型 已有 968 次阅读2010-11-6 14:31|个人分类:未分类|系统分类:科普集锦|关键词:辐射传输 转自https://www.360docs.net/doc/0c3499640.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/0c3499640.html,/ mxdwk <> 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简介 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公里处，放射性物质浓度的预测模型。 在遥感的实际应用中，常用很多简化的手段，如假设地面为朗伯面，排除云的存在，采用有关标准大气模式及大气气溶胶模式等，一次产生了许多不同类型的大气辐射传输模型，主要分为两类， 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辐射传热模型理论以及相关设置 目录 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;大气辐射传输模型
高斯扩散模型.
大气污染扩散模型
激光大气传输仿真系统的设计
Fluent辐射传热模型理论以及相关设置
遥感辐射传输模型
大气辐射传输理论 第一章..
污染物扩散模型-深圳数学建模
大气污染物扩散高斯模型模拟
大气辐射传输模型6S简介
污染空气的扩散模型
大气辐射传输校正模型(5S,modtran,acorn)
数学建模高斯扩散模型培训资料
Fluent辐射传热模型理论以及相关设置
大气辐射传输模型及其软件
点污染源空气污染扩散模型