matlab将坐标轴细分然后网格化,基于Matlab建立抚顺市高分辨率网格化大气污染源清单...

最新推荐文章于 2023-02-22 20:05:35 发布
最新推荐文章于 2023-02-22 20:05:35 发布
阅读量505 收藏
点赞数
文章标签： matlab将坐标轴细分然后网格化

秦思达

摘要：本研究基于Matlab对抚顺市大气污染物排放量进行空间分配、时间分配、物种分配，建立可用于空气质量数值模拟的城市高分辨率网格化清单。该网格化清单的建立可以作为基础数据支撑抚顺市开展空气质量数值模拟、重污染天气预报预警、污染物来源解析等研究工作，对抚顺市制定空气质量达标、改善方案，评估大气治理措施实施效果有重要意义。

关键词：Matlab；网格化；污染源清单；抚顺市

伴随我国社会高速发展，大气环境问题越发受到人们关注。应用排放源处理模式，建立区域全口径的大气污染源清单，并实现其网格化，对开展大气污染研究有重要意义。目前，排放源模式的主要作用是将以年为时间单位、以行政区为空间单位的排放清单进行空间和时间上的分配，根据模型选用的化学机制，对排放清单中的初始污染物进行化学物质种类和质量比例分配，以满足空气质量模型对于排放清单在时空分辨率和化学物种方面的高精度要求。配合第三代空气质量模式应用最广的是由美国北卡罗来纳微电子中心开发的SMOKE(Sparse Matrix Operator Kernel Emissions Modeling System)模式[1]。SMOKE模式是基于美国国家排放清单(NEI)设计的，在我国并不适用。首先，SMOKE模式空间分辨率精确到县，但我国尚未建立覆盖全国范围，精确到县市一级完整的一次污染物排放清单，在城市尺度上也只有少数几个工作基础较好的城市具有较为完善的排放清单[2]。SMOKE模式中的排放源以美国SCC code分类，精确到各种设备，而我国排放源分类体系与SCC code存在较大差异，且由于基础数据缺乏，我国区域排放清单估算时源分类很难如此精确。本研究应用Matlab编译建模，搭建污染源清单处理模型，可以实现多维清单、多物种的嵌套、融合，建立可用于第三代空气质量模式，且符合我国实际情况的网格化污染源清单，并将其应用在抚顺市的网格化清单建立中。

一、空间分配

污染源空间分配总体可以分为点源与面源两种，点源可以按照经纬度坐标分配到对应的网格上，而面源的分配较为复杂，采用代用参数权重法进行分配，不同的面源分配权重不同。抚顺市的点源主要为工业源，包含供暖源、电力源、建材源、石化源、钢铁源等源类，可以直接使用修正后的经纬度坐标分配。移动源、居民源、农业源属于面源，其中移动源按照单位网格路网所占面积的大小分配权重值；居民源(散煤燃烧)使用农村人口进行分配；居民源(溶剂使用、餐饮油烟)使用全口径人口数据进行分配；居民源(秸秆焚烧)、农业源使用土地利用类型中的耕地分布进行分配。

二、时间分配

污染物的时间分配是将污染源的年排放量分配到小时，制作成以小时为单位的清单输入模型。分配基于月分配系数、周分配系数、日分配系数处理完成。工业源与电力源排放受国家政策和市场调控影响，下半年的污染物排放略高于上半年。供暖源排放集中在11月至次年3月期间，居民源中的散煤燃烧主要在冬季供暖，部分用于日常生活，因此采暖期污染物排放量大于非采暖期。移动源的污染物排放季节性变化较小，但日变化特征显著，早晚高峰期的排放量较大。

三、物质分配机制

CBM、SAPRC、MCM以及RADM、RACM機理已经被广泛应用在很多的空气质量模式以及气象模式中。CBM机理在空气质量模式、城市气质模式、区域氧化物模式中的应用，SAPRC机理在空气质量模式和城市尺度的欧拉网格模式，MCM机理在光化学轨迹模式，RADM在EPA区域酸沉降模拟中的应用。CMAQ模型可选用cb05、saprc99、saprc07等多种化学反应机制，本研究考虑了海盐气溶胶排放、水/云化学的气相化学反应，选用cb05_ae6_aq作为模式化学反应机制。

cb05反应机制广泛应用于国内的空气质量数值模拟值，在对PM2.5与O3模拟方面具有良好的性能。Cb05是美国EPA于2005年发布的版本，共包括51个物种和156个反应。CB05根据国际纯粹与应用化学协会(International Union of Pure and Applied Chemistry，IUPAC)和NASA提供的最新反应速率常数对机理进行了更新。无机反应方面CB05的主要修改内容包括：1)增加了涉及H2的化学反应以适应较高的对流层环境；2)增加了奇氧原子的反应；3)增加了反映夜间化学过程的NO3反应；4)增加了NOx的反应来表现NOx的多天变化情况。有机反应方面CB05给出了甲烷和乙烷的详细化学机理，将高阶醛类细分为乙醛(ALD2) 和更高的醛 (ALDx) 并加入了几种萜烯物种(TERP)等。

四、抚顺市网格化清单

研究基于Matlab，对抚顺市大气污染物排放量进行了空间分配、时间分配、物种分配机制选择，最终建立了抚顺市空分辨率的大气污染源网格化清单，如下图1所示。

五、结论

本研究网格化清单的建立可以作为基础数据支撑抚顺市开展空气质量数值模拟、重污染天气预报预警、污染物来源解析等研究工作，对抚顺市制定空气质量达标、改善方案，评估大气治理措施实施效果有重要意义。

参考文献：

[1]CLEVE B M，张志涌编译.MATLAB数值计算[M].北京：北京航空航天大学出版社，2013：1-50.

[2]张礼俊.基于Model-3/CMAQ的珠江三角洲区域空气质量模拟与校验研究[D].华南理工大学，2010.

weixin_40001442
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
matlab细化坐标系,MATLAB化坐标系（转载的）
weixin_30983563的博客
03-19 1050
plot，plotyy是常用的画图函数，但是默认的格式通常不能准确的画出你想要的图，因此需要很多细节的修改，才能完整准确的表达你的意图。下面列举了一些常用的指令。[plain] view plain copy%单y轴plot(t*1e+9,abs(iGG)/max(abs(iGG)),‘k‘,‘linewidth‘,2);axis([-5,5,0,1])xlabel(‘时间/ns‘);ylabel...
数值分析Matlab三维数据网格图
Zhang Phil
12-13 7293
主要使用了matlab的meshgrid和mesh网格绘制函数。matlab代码：xi=-20:0.5:20; yi=-20:0.5:20; [x,y]=meshgrid(xi,yi); z=(-x.^2-y.^2); mesh(x,y,z);代码运行结果如图：
ffndgrid Fast 'n' Furious ND 数据网格化：FFNDGRID 网格不均匀间隔的 N 维数据。-matlab开发
06-01
FFNDGRID 网格不均匀间隔的数据。 这是适用于 ND 数据的 FFGRID 的概括。 FFNDGRID 也可用于生成 N 维直方图或密度图。 这在处理使用 plot3 难以查看的数据时非常有用（根据我的经验，所有不完全平滑的 3-D 数据都是这种情况）。 FFNDGRID 是一种 Fast ‘n’ Furious 方式，可以完成与 griddata 为您所做的相同工作。 不同之处在于没有插值。 空点留空，而不是尝试使用相邻点填充它们。 此外，FFNDGRID 对落在同一个网格单元中的多个点没有问题。 除非指定了输出参数，否则数据将使用 IMAGSC 显示，在这种情况下，（规则间隔的）数据矩阵与指定分箱的向量一起给出。
Matlab有限元网格化源程序-dcircle.m
08-13
Matlab有限元网格化源程序-dcircle.m MIT的人写的论文，一个简单的有限元网格化方法 确实可以用，但是我理解的不是很好。 不过，希望对需要的人有所帮助。 简单介绍： function [p,t]=distmesh2d This meshing function produces the following outputs: • The node positions p. This N-by-2 array contains the x, y coordinates for each of the N nodes. • The triangle indices t. The row associated with each triangle has 3 integer entries to specify node numbers in that triangle. The input arguments are as follows: • The geometry is given as a distance function fd. This function returns the signed distance from each node location p to the closest boundary. • The desired edge length function h is given as a function fh, which returns h for all input points. • The parameter h0 is the distance between points in the initial distribution p0. For uniform meshes = constant), the element size in the final mesh will usually be a little larger than this input. • The bounding box for the region is an array bbox=[xmin, ymin; xmax, ymax]. • The fixed node positions are given as an array pfix with two columns. • Additional parameters to the functions fd and fh can be given in the last arguments varargin .
带网格的_热点研究 | 2013至2017年空气污染导致过早死亡的高分辨率网格化分析...
weixin_34518190的博客
01-10 217
点击上方蓝字，关注我们?编者按：大量流行病学队列研究证明PM2.5长期暴露可导致过早死亡。使用高空间精度网格化分析，可减少健康负担的估算误差，为政策制定提供建议。研究背景中国作为全球受PM2.5污染严重的国家之一，空气污染对国民健康的影响更加不容忽视。现有研究大多使用低空间分辨率的PM2.5浓度数据，忽略了暴露响应关系参数的空间差异，使得健康负担估算更容易存在偏差。本研究利用混合遥感-地...
wy_scat2grid3.zip_matlab 点云_三维网格化_点云_点云MATLAB
07-15
wy_scat2grid3.m实现三维点云网格化
loopSubdivision.zip_LOOp细分_loop细分 matlab_matlab网格细分_thise8e_网格细分
07-15
loop网格细分算法的实现，使用matlab编写，直接运行test.m即可
Matlab源代码+用于多尺度有限元方法+网格化（源代码+详解手册+模型示例图片）
12-25
资源内容：源代码+详解手册+模型示例图片 运行环境：Matlab 2019b 功能说明： 包含用于多尺度有限元方法的matlab代码 我们实现了吴晓辉[论文]（docs/MsFEM.pdf）中提出的多尺度有限元方法。 查看我们的最终[演示...
xf.rar_ 三角网格_matlab 网格 图像缩放_三角 matlab_二维 MATLAB_细分
09-19
二维三角网格细分，并同时建立新的三角网格的拓扑结构
基于C++MATLAB跨平台编程，实现计算流体力学中结构化网格划分源码+全部资料.zip
最新发布
04-16
基于C++MATLAB跨平台编程，实现计算流体力学中结构化网格划分源码+全部资料.zip基于C++MATLAB跨平台编程，实现计算流体力学中结构化网格划分源码+全部资料.zip 【备注】 1、该资源内项目代码都经过测试运行成功，...
Matlab有限元网格化源程序-distmesh2d.m
08-13
Matlab有限元网格化源程序-distmesh2d.m MIT的人写的论文，一个简单的有限元网格化方法 确实可以用，但是我理解的不是很好。 不过，希望对需要的人有所帮助。 简单介绍： function [p,t]=distmesh2d This meshing function produces the following outputs: • The node positions p. This N-by-2 array contains the x, y coordinates for each of the N nodes. • The triangle indices t. The row associated with each triangle has 3 integer entries to specify node numbers in that triangle. The input arguments are as follows: • The geometry is given as a distance function fd. This function returns the signed distance from each node location p to the closest boundary. • The desired edge length function h is given as a function fh, which returns h for all input points. • The parameter h0 is the distance between points in the initial distribution p0. For uniform meshes = constant), the element size in the final mesh will usually be a little larger than this input. • The bounding box for the region is an array bbox=[xmin, ymin; xmax, ymax]. • The fixed node positions are given as an array pfix with two columns. • Additional parameters to the functions fd and fh can be given in the last arguments varargin .
Matlab有限元网格化源程序-huniform.m
08-13
Matlab有限元网格化源程序-huniform.m MIT的人写的论文，一个简单的有限元网格化方法 确实可以用，但是我理解的不是很好。 不过，希望对需要的人有所帮助。 简单介绍： function [p,t]=distmesh2d This meshing function produces the following outputs: • The node positions p. This N-by-2 array contains the x, y coordinates for each of the N nodes. • The triangle indices t. The row associated with each triangle has 3 integer entries to specify node numbers in that triangle. The input arguments are as follows: • The geometry is given as a distance function fd. This function returns the signed distance from each node location p to the closest boundary. • The desired edge length function h is given as a function fh, which returns h for all input points. • The parameter h0 is the distance between points in the initial distribution p0. For uniform meshes = constant), the element size in the final mesh will usually be a little larger than this input. • The bounding box for the region is an array bbox=[xmin, ymin; xmax, ymax]. • The fixed node positions are given as an array pfix with two columns. • Additional parameters to the functions fd and fh can be given in the last arguments varargin .
点云数据栅格化matlab代码
06-26
文档中.m文件为将txt文件的点云数据进行栅格化，直接可以运行。
输入ply点云进行三角网格化matlab
05-19
输入ply点云进行三角网格化matlab 输出三角网格化的效果
基于MATLAB的大气污染原型系统开发.pdf
10-30
基于MATLAB的大气污染原型系统开发.pdf
matlab画离散图加入网格,Matlab关于画数据网格图
weixin_39937635的博客
03-17 684
1、surf(Z)；Z为一个矩阵，如果Z是向量，那么就需要先将其转换为矩阵reshape(Z,[length(X),length(Y)]); X和Y为你的横坐标取值个数。具体代码调用格式如下Z=reshape(acc,[length(lambda2),length(lambda1)]);surf(Z);xlabel('\alpha','fontsize',20);set(gca, 'xtick'...
Matlab：插入网格数据
希望我的博客，能帮上你解决学习中工作中所遇到的问题
02-22 419
Matlab：插入网格数据
matlab可以实现克里格,插值(克里格法?)MATLAB中的网格数据
weixin_42511151的博客
03-16 723
我有Lat,Lon和NO2的数据是网格化的。我试着从这些数据中选择一些点,然后插值出一些数据/图,比较插值后的数据/图在插值过程中的改进。我的数据如下:Lat尺寸为336x264Lon的尺寸是336x264NO2(:,:,1,12)(我将使用的变量)的大小是336x264,前两个维度分别是lat和lon。例如,lat(1,1)、lon(1,1)和NO2(1,1,1,12)对应。目前,为了插值这些数...
matlab确定污染源位置,确定污染源的位置
weixin_29571873的博客
03-17 1108
确定污染源的位置 5 污染源的判定 在问题(2)中，通过考虑不同的条件，建立了三个模型，首先单独作分析，然后综合三个模型得出结论，具体情况分析如下： 1、差分方程反演模型的建立与求解 1.1 模型I的建立 问题II研究和分析的是长江干流近一年多来受高锰酸钾盐指数及氨氮的污染源情况，建立该模型主要考虑以下几个因素： 一、浓度是随时间、距离变化的量。 二、长江中的水流量及水中污染物在消耗减少，同时也在...
matlab将xy坐标网格化
05-28
在 Matlab 中，可以使用 `meshgrid` 函数将 x 和 y 坐标向量转换为网格。具体来说，`meshgrid` 函数将 x 向量复制为列向量，y 向量复制为行向量，然后生成两个网格矩阵 X 和 Y，其中 X 的每一列都是 x 向量的复制，Y...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值