放射性气体扩散浓度的预测模型.doc

《放射性气体扩散浓度的预测模型.doc》由会员分享，可在线阅读，更多相关《放射性气体扩散浓度的预测模型.doc（14页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、【精品文档】如有侵权，请联系网站删除，仅供学习与交流放射性气体扩散浓度的预测模型.精品文档.放 射 性 气 体 扩 散 浓 度 的 预 测 模 型参赛队员：程伟山 15138745435 张 聪 13461028984 赵雅静 15138763238放射性气体扩散浓度的预测模型摘 要关于放射性气体扩散浓度的预测，本文根据日本福岛核电站遭遇自然灾害发生核泄漏的背景进行研究，首先用线性回归模型进行简单预测，其次建立了高斯模型解决以下问题，并根据实际情况对模型进行修正与优化。针对问题一，考虑放射性气体的泄漏是连续不断的。首先利用spss软件根据福岛核电站周围的放射性气体扩散的实际情况建立一个线性回归

2、模型，可以简单预测出福岛核电站西南方的最大影响范围。接着根据“泄漏放射性物质质量守恒定律”和“气体泄漏连续性原理”建立了微积分方程，应用了数学分析中散度、梯度、流量等数学概念，通过Guass公式、四维二阶偏微分方程，因而得到了核电站周边不同距离地区、不同时段放射性物质浓度的预测模型。同时为使模型适用范围更广，本文引入了地面反射系数，考虑由于放射性气体从泄漏口喷出时具有初动量而使其泄漏源有效高度被抬高等因素，所以我们引用了在无风环境中适用范围更广的“高斯修正模型”。针对问题二，主要探究风速对放射性气体浓度分布的影响。本文运用概率论知识，通过图解和数学推导得出“连续点源放射性物质高斯扩散模型”。依

3、次考虑了“重力沉积”、“雨水沉积”、“核衰变”等因素对浓度分布的影响。并通过构建“耗减因子”、“衰变因子”等方法将耗减和衰变的放射性物质“投影”到泄漏源浓度中，得到了经多次合理修正后的“优化高斯模型”，并详细分析了泄漏源周边地区放射性气体的浓度变化情况。针对问题三，本文主要根据问题二的解答模型，结合考虑风速和放射性物质扩散速度在空间中的矢量运算。得出在对上风口、下风口分析时，必须分类讨论风速和自然扩散速度之间的大小关系，当风速大于放射性气体扩散速度时，放射性气体是无法到达上风口的。针对问题四，本文查阅整理了大量气象、地理、新闻等数据资料，以福岛核电站西南方向我国上海作为研究对象，综合考虑对应海

4、域平均风速及风向、地理距离、海水对放射性气体扩散的反射系数等因素，并通过MATLAB编程进行模拟计算，预测出福岛核电站放射性核物质大约经过一周到达我国东海岸，且预测得到131I浓度预测值约为：0.100，即大约与实际情况比较吻合，对人影响已很小。同时对美国西海岸也就行了预估，由于影响因素较多，预估结果相对较差。关键词：放射性气体的扩散 高斯修正模型 浓度预测一、问题重述突发性核泄漏等紧急灾害事件具有爆发性、空间分布不连续性、对周边地形和气象条件的敏感性的特点，所以研究核事故所释放的气体物质的时空扩散需要高度精确的技术，这对于更好地保护环境有着极其重要的意义。现假设有一座核电站遇到自然灾害而发生

5、泄漏，浓度为的放射性气体以匀速排出，速度为，在无风的情况下，匀速在大气中向四周扩散，速度为。建立模型解决以下问题：问题一，建立一个描述核电站周边不同距离地区、不同时段放射性物质浓度的预测模型，这对于研究核污染模式具有重要的意义。 问题二，当风速为时，给出核电站周边放射性物质浓度的变化情况，这对于研究核电站附近放射性气体对实际环境的影响有着重要的作用。 问题三，当风速为时，计算出上风和下风公里处的放射性物质浓度的预估。 问题四，将建立的模型应用于福岛核电站的泄漏，计算出福岛核电站的泄漏对我国东海岸以及美国西海岸的影响，根据这一实际预测，可以明确某一地区实际的污染情况，为我们的核应急决策提供理论数

6、据。二、问题分析针对问题一，根据假设在无风的情况下，放射性气体以的速度，匀速在大气中向四周扩散。在此条件下，建立线性回归和高斯模型来对核电站周边不同距离地区、不同时段放射性物质浓度进行预测。首先我们明确此问题研究的放射性气体扩散是点源连续泄露的扩散问题。虽然题目只是要求考虑在无风情况下放射性气体浓度分布，但为了使模型更贴切实际，更具有科学性，本文还考虑了地面反射、泄漏源有效高度等因素对浓度分布的影响。根据“泄露放射性物质质量守恒定律”和“气体泄漏连续性原理”进行分析，得出核电站周边不同距离地区、不同时段放射性物质浓度的预测模型，最后对于该方程进行分析求解。对于问题二，为了探究风速对发生核泄漏的

7、核电站周边放射性物质浓度分布的影响，运用概率论知识，通过图解和数学推导得出“连续点源放射性物质高斯扩散模型”。应在“连续点源放射性物质高斯扩散模型”的基础上经多次合理修正后得到更好的“优化高斯模型”。对于问题三，该问题要求建立泄漏源上风口和下风口处放射性物质浓度的预测模型，在参考第二问的基础上，主要考虑风速和放射性物质扩散速度在空间中的矢量运算。在对上风口分析时，要分类讨论风速和自然扩散速度之间的大小关系，当风速小于自然扩散速度时，放射性物质是无法到达上风口的。对于问题四，参考了大量气象、地理、新闻等数据资料，选择我国东海岸上海，作为研究对象，综合考虑对应海域平均风速及风向、地理距离、海水对放

8、射性物质扩散的部分反射系数等因素，预测出放射性核物质与实际情况比较。三、模型假设3.1模型假设(1) 放射性气体浓度在y、z风向上分布为正态分布；(2) 整个扩散过程中风速的大小、方向保持不变；(3) 地面对放射性气体的反射是不完全的；(4) 气云中不发生化学反应，地面对气云无吸收；(5) 假设气体的传播服从扩散定律，既单位时间内通过单位法向面积的流量与她的浓度梯度成正比。四、符号说明4.1符号说明:放射性气体的传播速度；:风速，单位；:泄漏点距地面的实际高度；:泄漏点距地面的有效高度；:扩散时间；:空间任意一点在时刻时的放射性物质浓度；:空间任意一点的放射性物质的扩散系数；:空间域；:空间域

9、体积；: 一规则的球面面积；: 单位时间通过单位法向面积的流量；:在内通过的流量；:内放射性物质的增量；:从泄漏源泄漏的放射性物质的总量；:附加高度；:核泄漏出口处的温度；:环境温度；:地面反射系数；:源强；:分别为用浓度标准差表示的轴上的扩散参数；:沉降速度，单位为；:地面干沉积率；:冲洗系数；:放射性核素的半衰期；五、模型的建立与求解5.1问题一模型的建立与求解5.1.1 线性回归模型的建立假设全球大气处于不对流状态，核物质会以恒定的速度向四周扩散，辐射源的核辐射量近似作为核辐射总量。在地理位置上，查的 3月25号水中量在日立-茨城-金町-东京新宿-横滨小雀一线的分布数据如下表1：表1城市

10、放射性离福岛核电站距离日立59.07100.05茨城47.88165.0金町38.04198.78东京新宿32.76215.37横滨小雀23.01270.03利用spss软件得到辐射放射性与距离之间的散点图：图1由散点图得水中I-131量在日立-茨城-金町-东京新宿-横滨小雀一线的分布具有一定的相似性（见图1）。由此可先建立线性模型进行简单的预测。表2Model SummaryModelRR SquareAdjusted R SquareStd. Error of the Estimate1.994a.988.9841.75443a. Predictors: (Constant), 距离Coe

11、fficientsaModelUnstandardized CoefficientsStandardized CoefficientstSig.BStd. ErrorBeta1(Constant)81.7712.76029.625.000距离-.219.014-.994-15.727.001a. Dependent Variable: 辐射可得线性模型为： =-0.219+81.771 (1)其中， 为距离。拟合图为：图2对方程（1）求解可知当距离为373.38时，核辐射量就降为0，而在西南方向，我国上海距离日本辐射源的距离为大约为2273公里，也就是日核泄漏物不会大量的传到我国。5.1.2高

12、斯模型的建立以核泄漏点正下方的地面为坐标原点，指向下风方向为轴，为横风向，垂直方向为轴，建立空间坐标系，则核电站泄漏点距地面的有效高度为，则泄漏点位置坐标为。图3 空间坐标系示意图记时刻，空间任意一点的放射性物质浓度为。根据假设设单位时间通过单位法向面积的流量与浓度梯度成正比，有：是扩散系数，表示梯度，负号表示由浓度高向浓度低的地方扩散。考察空间域，其体积设为，包围的曲面为，在无风的情况下，为一规则的球面，外法线向量为。则在内通过的流量可记为：在内空间域内放射性物质的增量（即存量）为：同时间内从泄漏源泄漏的放射性物质的总量为：根据质量守恒定律和连续性原理，泄漏源在单位时间内向外泄漏的放射性物质

13、等于单位时间内通过所选曲面的向外扩散的放射性物质与曲面内放射性物质增量之和，即有：即：又根据曲面积分的Gauss公式： (其中是散度记号)所以 因为由以上两式得：化简得： 另根据参考文献 解得： (2)上述模型仅是一个最理想化的预测浓度模型，因为它将一切环境视为无边界空间，且没有考虑放射性物质从泄漏口释放出时的初动量、外力、地面反射以及防身性气体自身反应等因素。为使所建立模型更加贴合实际，为实际应用服务，以下将从“有效泄漏源”和“地面反射”这两个方面对所建立模型进行修正。5.1.3考虑热力抬升作用对模型的修正如图4所示，为核泄漏点源的有效源高。它是由两部分构成：一是核泄漏口的实际高度；二是在实

14、际核扩散中核泄漏气团从泄漏口排出时，由于受到热力抬升和本身动力抬升，进而产生的一个附加高度。H为有效高度。图4 热力抬升示意图对于，主要由泄漏的初始动量和浮升力决定，同时还要受到泄漏口温度、大气温度、风速、地形地貌等多种因素的影响。我们直接引用气体污染扩散学中应用较广范的，有关烟气抬升高度的综合分析公式：其中： 泄漏源的实际高度；泄漏源出口处的风速，已知为m/s；泄漏源出口的有效直径；放射性气体的扩散速度，已知为m/s；浮力通量，。由Briggs抬升公式知，浮力通量，其中为核泄漏出口处的温度；为环境温度。由我国“制订地方大气污染物排放标准的技术方法”(GB/T13223-91)中的公式标准：1

15、）在有风()且释放气体温度与环境温度差35K()时，抬升高度：其中为大气压强，2）小风()时，且温度差35K()，抬升高度为：其中为泄漏源的有效高度处上的环境温度梯度()。3）当温度差35K()时，此时的抬升高度：此时，当风速时，取= 结合以上分析，我们将以上结论应用到本题中（带入风速及放射性气体的扩散速度）得：综上所述，泄漏源的有效高度为：5.1.4考虑到地面反射对模型的修正其次考虑到地面会对扩散的放射性气体有一定的反射作用，同时扩散的核素粒子受沉降等作用，核扩散物质又不可能被全部反射回去，因而核辐射物质只能是部分反射回大气，为了便于描述和模拟，设地面反射系数为。这样进入大气的核辐射物质可以

16、看成是两个部分：一是从泄漏源直接扩散到空间点；二是从地面反射进入空间点（见图5）：图5 连续点源扩散地面反射部分示意图（）从几何物理学分析可知，通过地面反射进入大气的核扩散物质反射相当于虚拟泄漏源泄露的核辐射物质在原来空间的一个浓度叠加。如果所设核泄漏点源在距有效地面的高度为的地方，本文求空间任一点的浓度值，则实际泄漏源对点的影响部分可用来表示。在考虑反射系数后，虚泄漏源对点的影响部分可用来表示，于是(2)式所得模型可以修正为：5.2 问题二模型的建立与求解5.2.1高斯连续点源扩散模型的建立问题二，针对在有风速时，选取无界空间高斯连续点源扩散模型进行分析，这也是在所有有风速时，放射性物质扩散

17、模型中最常见也最方便的一种泄漏方式。放射性气体云团在大气中实际迁移和扩散的数值计算基本上可分为二步：第一步：关于所关注区域中风场的计算（理论基础是大气运动方程、连续性方程、状态方程、热力学方程和水汽方程构成的基本方程）。第二步，进行已知风场中放射性气体云团迁移和扩散的计算，可采用类似于处理大气污染的方法，假设放射性气体云团不影响大气流体速度和温度，求解出放射性气体云团的连续性方程。高斯连续点源扩散模型有边界点源，这里设泄漏源有效高度为，取其地面投影为坐标原点，轴指向风向，继续考虑地面反射作用，图6 边界点源示意图可得到高斯连续点源泄漏的浓度分布7：其中，为源强，；为污染物质量浓度，；分别为用浓

18、度标准差表示的轴上的扩散参数；为泄漏高度的平均风速，m/ s；为泄漏有效高度，。考虑到动力抬升的作用，将有效高度（即有效烟云高度）图7 动力抬升示意图此时（18）式可以写成为： (3)5.2.2对高斯浓度计算模型的修正利用标准的高斯模型，计算大气中放射性核素云团的扩散，目前已有较多的研究，而实际的核素扩散过程还存在着粒子的重力沉降、雨洗沉积以及核素衰变等的因素对浓度分布的影响。因此，有必要对高斯浓度扩散模式进行修正，方可较为准确、真实地反映实际的核素扩散规律。由于地面风速对大气扩散的影响起着至关重要的作用，不同的风速其浓度计算方法有较大差异，因此，下面的分析从有风（大于）、小风（小于大于）和静

19、风（小于），三种情况分别论述。1）考虑干沉积时的连续点源扩散(1)粒径小于10的颗粒物可按气体扩散计算(2)大于10的颗粒物有明显的重力沉降，粒子的沉降速度取决于空气阻力和重力平衡，可用斯托克斯公式表示： (4)其中：粒子密度， ；：重力加速度，9. 806 5；：粒子直径，;：空气的动力粘性系数, 可取；：沉降速度，含碘放射性核素的干沉积速度为= 1. 1 , 该值也可由式(4)计算得到。因为在扩散过程中同时有重力沉降的位移迭加到羽流中心线上，中心线就会向下倾斜， 所有粒子相当于在下倾的中心线上扩散。该类扩散和沉降的迭加可认为是羽流运行过程中，实源以的速度向下移动，在处向下移动的高度为，即源

20、高由降到了。实际上，由于大气湍动及其他动力作用， 地面不是全吸收表面，应考虑地面的反射作用。但毕竟存在粒子的沉降作用，又不是全反射，因此，需对反射项乘以一个反射系数，由于反射项的有效源高度也变成了，故相应的浓度计算公式为： (5)式(5)可得到经过干沉积作用后的浓度分布结果。式(5) 中, 反射系数需外部给定， 通常对于放射性核素可取0.5。2）小风和静风（小于）时的计算参考小风和静风时的浓度计算方法，再利用部分反射倾斜烟云模式，可得出该情况下的浓度分布：通常按照每小时排放6个烟团进行积分，即可达到需要的精度。3）地面的干沉积量根据扩散理论和动量传递的普朗特理论，可以导出干沉积的地面沉积为：式

21、中，前面各种情形计算的地面污染物的浓度；干沉积速度，；地面干沉积率，。4）雨洗作用降雨对烟羽中的颗粒物及气溶胶具有清洗作用，可溶性气体与蒸汽亦可溶于雨水中，降雨过程造成的这类湿沉积是导致放射性气溶胶和气体向地面沉积的另一重要机制。通常以冲洗系数，描述降雨对烟羽中污染物清洗作用的大小。与雨强的关系可以表达为： 式中，为雨强() ；为经验系数。式中，按释放物质为含碘、不含碘情况分别取值。对含碘物质，；对于不含碘物质，取。对于湿沉积导致的烟羽耗减，可采用湿沉积耗减因子对源强进行修正，有5）放射性衰变的影响除了干沉积和湿沉积外，放射性物质的衰变也是影响大气中核素浓度分布的主要原因。由于放射性物质服从简

22、单的衰变规律，其浓度随时间的变化可由下式计算：式中， 为初始浓度；衰变常数；经过的时间；可由无衰变的浓度公式计算。参照前面的模式，同样采用衰变耗减因子来对源项进行修正。经过理论推导，可以得出：其中，放射性核素的半衰期。如以上情况都存在，应同时考虑源项修正，即式5.2.3高斯浓度计算模型的求解1）模型参数的确定在事故过程模型中所需参数的选取及确定十分重要，通常情况下气象参数的选取是利用该地区多年气象资料，采取工业安全与环保统计的方法进行有关参数的确定，而其他扩散参数是以实际测定为基础的2）大气稳定度的计算根据国家标准(GB/ T 13201 1991)制定地方大气污染物排放标准的技术方法的规定，

23、划分大气稳定度的级别，共分为6级A F，A为极不稳定；F为极稳定。首先，根据释放源所在地的经度和纬度以及泄漏的日期和时间计算当时的太阳高度角；然后，由太阳高度角和云量查出太阳辐射等级；最后，再根据地面风速确定当时的大气稳定度，计算细节可参考文献。表 3 大气稳定度的级别参考表地面风速白天太阳辐射阴天的白天或夜间有云的夜晚强中弱薄云遮天或低云0.5云量0.46CDDDDD3）扩散参数的计算有风时的扩散参数，的确定采用Briggs给出一套扩散参数，如表4和表5所示。表 4 Briggs 扩散参数（开阔平原田野）大气稳定度ABCDEF表5 Briggs 扩散参数（城市）大气稳定度A-BCDE-F5.

24、2.4模型的求解 图8 模型二计算流程图5.3问题三模型的建立与求解5.3.1高斯连续点源扩散优化模型的建立问题三，针对在有风速时，选取高斯连续点源扩散模型进行分析在给出的上风和下风公里的距离处在结合，在模型二的基础上考虑浓度等效圆的情况。在模型二里面将式(3)中的风速在上风和在下风不同情况下与传播速度之间的比较的分析。1）在上风的情况下，此时满足（）此时若时，即风速小于自然扩散速度，此时的放射性物质是无法到达上风口的，则必须满足风速大于自然扩散速度，仍然根据模型二可以求解，只是将式(3)中的用来代入。2）在上风的情况下，此时无较大的影响和要求此时不论风速小于自然扩散速度，还是风速大于自然扩散

25、速度，此时的放射性物质都可到达下风口的，则仍然根据模型二可以求解，步骤同5.2.2.图9核泄漏对上下风在距离上的影响5.4问题四的求解2011年3月15日下午，日本福岛核电站发生放射性气体泄露，查询日本中北部区域在中低层大气中的风向由西南风转为西北风，高空大气主要以偏西风气流为主，同时近期由于降水发生，有利于核物质沉降，影响范围缩小。未来三天（16日至18日），日本核电站核泄漏产生的放射性污染物主要影响区域为日本中部、北部及其以东的北太平洋区域。在结合问题二、问题三的条件下，在参阅整理大量气象、地理、新闻资料，以我国东海岸上海作为研究对象，综合考虑对应海域平均风速及风向、地理距离、海水对放射性

26、物质扩散的部分反射系数等因素，并通过Matlab编程模拟计算，预测出放射性核物质将经过大约一周到达我国东海岸，且浓度值约为：0.100，与实际情况测得数据比较吻合，验证了模型的准确性。由于到达美国西海岸途径太平洋，距离较远，风速天气等因素具有不确定性，所以预测会产生很大误差，预估将经过20天左右到达，且浓度已不高。与实际情况虽存在误差，但还可以接受。5.5模型的检验本题由于在前四问中均无具体的数据，所以对于模型预估结果，在检验上存在一些麻烦。但是在模型三预测上下风的浓度具有较好的实际意义，同时通过对问题四的解答，对于预测出浓度的变化和对我国的影响也较为合理。故该模型具有必要性、可靠性和参考性，

27、可以进行推广和其它方面的应用。六、模型误差分析 由于本文建立模型前有部分理想化假设，这些假设都会造成模型预测结果与实际情况产生误差，比如放射性气体扩散过程中风向的变化，这是一个很实际的影响因素；同时，在海上扩散时，我们更应该考虑海水与海洋气候的因素；所以我们在建立模型时，尽量加入各方面因素，比如地面反射、雨洗作用、放射性衰变等，对高斯模型进行不断修正，使之与实际情况相差在所能接受的误差范围之内。最后通过对上海的预测，与实际情况的相差不大，所以整体来说，本模型具有一定的准确性，可行性。七、模型的评价与推广通过检验整个模型，有以下优点很突出：该方法可计算核事故中连续点源在不同气象、地形条件下的浓度

28、分布。该结果在核事故的应急救援过程中，对救援人员划定警戒区和确定周围居民的疏散范围具有重要意义，并可为制定救援方案和应急决策提供科学依据。放射性云团在空中迁移和扩散提供浓度的定量描述，对可能发生的核事件的放射性核素浓度监测及监测时间范围提供相关信息。不可避免的，模型也存在一点缺点，模型在建立时忽略的地形的复杂性，实际的扩散情况可能因地形条件的复杂化而有所改变。由于扩散是一个非常复杂、影响因素众多的过程，上述方法还有一定的局限性，有待于进一步检验。另外，次模型课通过浓度的计算结果，换算成人的吸收剂量，这样可以更具体的研究对人的影响。今后将对其计算过程和参数的选取进行多方面验证，使其进一步完善。参

29、考文献1金勇进，统计学M ，北京：中国人民出版社，2010。2梅长林，周家良，实用统计方法M，北京：科学出版社，2002。3罗艾民，魏利军，有毒重气泄露安全距离数值方法J，中国安全科学学报，2005,15(8)：98-100 。4周燕等，Matlab在统计与工程数据分析中的应用M，北京：电子工业出版社，2010。5丁信伟，王淑兰，徐国庆，可燃及毒性气体泄漏扩散研究综述J，化学工业与工程，19992(16)：118122。6胡世明，气体释放源的三维瞬态重气扩散研究J，劳动保护科学技术，20023(20)：2830。7姜启源，谢金星，数学模型(第三版) M，北京：高等教育出版社，2004。8陈宗则

30、, 熊桄清 ,张泽荣 ，核电站安全性分析与发展预测M。 9孙大明，新一代大气扩散模型(ADMS)应用研究M。10冯有前，数学实验M，国防工业出版社。附 录附录一：模型二在上下风L公里的浓度求解源代码function c1,c2=qiunongdu(x,y,z)PI=3.1415926;x=input (L处坐标（米）：);y=input (地点的纵坐标（米）：);z=input (地点的铅直坐标（米）：);t=input (扩散时间（秒）：);a=input (地面反射系数：);s=input (放射性气体扩散速度（m/s）：);k=input (风速（m/s）：);I=input (雨水强度

31、（毫米/小时）：);Ts=input (核物质衰变周期（天）：);p0=input (初始浓度（千克/立方米）：);H=input (有效高度（米）：);daqi=input (大气稳定度等级：);switch(daqi) case A dy=0.22*x/(sqrt(1+0.0001*x); dz=0.2*x; case B dy=0.16*x/(sqrt(1+0.0001*x); dz=0.12*x; case C dy=0.11*x/(sqrt(1+0.0001*x); dz=0.08*x/(sqrt(1+0.0002*x); case D dy=0.08*x/(sqrt(1+0.000

32、1*x); dz=0.06*x/(sqrt(1+0.0015*x); case E dy=0.06*x/(sqrt(1+0.0001*x); dz=0.03*x/(1+0.0003*x); case F dy=0.04*x/(sqrt(1+0.0001*x); dz=0.016*x/(1+0.0003*x);endc1=p0/(2*PI*(s-k)*dy*dz)*exp(-(0.693*x)/(3600*(Ts*24)*(k+s)-(8*1e-5*I0.6*x)/. (k+s)*(exp(-(z-H)+(0.01*x/(s-k)2)/(2*dz*dz)+a*exp(-. (z+H+(0.01*

33、x/(s-k)2)/(2*dz*dz);c2=p0/(2*PI*(s-k)*dy*dz)*exp(-(0.693*(-x)/(3600*(Ts*24)*(k+s)-(8*1e-5*. I0.6*(-x)/(k+s)*(exp(-(z-H)+(0.01*(-x)/(s-k)*(z-H)+(0.01*(-x)/(s-k). /(2*dz*dz)+a*exp(-(z+H+(0.01*(-x)/(s-k)*(z+H+(0.01*(-x)/(s-k)/(2*dz*dz);disp(下风口L处放射性物质浓度);c1disp(上风口L处放射性物质浓度);c2附录二：核泄漏模拟图的绘制Q=1;%输入强源u=2

34、;%输入风速d=1;%步长Z=0.4; %地面粗糙参数x=10:1:100; %下风向距离 y=-100:1:100; %横风向距离x,y=meshgrid(x,y);by0=0.08*x.*(1+0.0001*x).(-1/2);bz0=0.06*x.*(1+0.0015*x).(-1/2);by=by0.*(1+0.38*Z);fz=(2.53-0.13*log(x).*(0.55+0.042*log(x).(-1).*Z.(0.35-0.03*log(x);bz=bz0.*fz;tempy1=-y.*y./by./by./2;tempy2=exp(tempy1);c=Q/pi/u*(by.*bz).(-1).*tempy2Cs=150;%输入求解条数 contour(x,y,c,Cs);shading interp;colorbar;grid;title(核泄漏模拟图)