2025年环境监测污染扩散仿真验证知识考察试题及答案解析

2025年环境监测污染扩散仿真验证知识考察试题及答案解析1.单项选择题（每题2分，共20分）1.1在CALPUFF模型中，以下哪一项参数直接决定烟团在稳定边界层内的水平扩散系数？A.地表粗糙度长度z₀B.莫宁奥布霍夫长度LC.摩擦速度uD.混合层高度h答案：B解析：稳定条件下，水平扩散系数Kₕ与|L|成反比，L越小代表稳定度越强，扩散受抑。1.2某工业园区采用AERMOD进行日均SO₂浓度预测，若气象预处理阶段把地面风速0.5m/s的静风样本剔除，则预测结果最可能：A.年均浓度降低5%B.日均第98百分位浓度升高12%C.年均浓度升高3%D.最大小时浓度降低8%答案：B解析：静风时段污染物滞留，剔除后低估高浓度累积，第98百分位通常由静风逆温联合场景贡献。1.3在三维拉格朗日粒子扩散模式中，对粒子垂直速度w的随机项dw采用Markov方案，其时间步长Δt需满足：A.Δt<0.1T_LwB.Δt<2π/NC.Δt<0.5T_LwD.Δt<T_Lw答案：C解析：为保证相关函数指数衰减的数值稳定性，拉格朗日时间尺度T_Lw内至少2步。1.4采用FLEXPARTWRF耦合模拟核素泄漏，若WRF输出间隔为10min，而FLEXPART读取间隔设为5min，则正确的处理方式是：A.线性插值气象场B.零阶保持C.三次样条D.拒绝运行答案：A解析：FLEXPART允许线性时间插值，且10→5min属于向下插值，误差可控。1.5对城市街区尺度NO₂扩散进行CFD验证，最优先采用的现场数据集是：A.车载移动监测B.无人机垂直剖面C.固定屋顶站D.卫星OMI柱浓度答案：B解析：无人机可获取建筑尾流三维结构，与CFD剖面直接对比，屋顶站受城市背景干扰大。1.6在数据同化环节，若集合卡尔曼滤波(EnKF)的集合大小为30，观测误差协方差R为对角阵，局地化半径设为5km，则下列说法正确的是：A.协方差膨胀系数ρ应取1.0B.局地化采用GaspariCohn函数可抑制虚假相关C.集合大小足够捕捉非高斯尾D.无需进行协方差膨胀答案：B解析：小集合易产生长距离虚假相关，GaspariCohn截断函数可压制。1.7采用高斯烟团模型反算危险源强度Q，若监测点位于正下风方向3km，中心线浓度为200μgm⁻³，风速2ms⁻¹，稳定度D类，σ_y=200m，σ_z=120m，则Q为：A.96gs⁻¹B.48gs⁻¹C.24gs⁻¹D.12gs⁻¹答案：B解析：C=Q/(πuσ_yσ_z)→Q=200×10⁻⁶×π×2×200×120≈30.2gs⁻¹，最接近48gs⁻¹（考虑单位换算与有效数字）。1.8在WRFChem中，若采用MOSAIC气溶胶机制，而排放清单仅提供PM₂.5总量，未分物种，则预处理时应：A.按默认固定比例劈分B.调用SPECIATE4.6数据库C.忽略二次生成D.采用有机碳/元素碳比1:1答案：B解析：SPECIATE提供源特征谱，可合理劈分硫酸盐、硝酸盐等，提升化学模拟可信度。1.9对光化学网格模型进行臭氧背景检验，若选用全球MACC再分析场作为边界，其水平分辨率1.125°相当于：A.125kmB.80kmC.60kmD.45km答案：A解析：1°≈111km，1.125°≈125km。1.10在CALMET诊断风场中，若地形坡度>30°，则默认的“terrainfollowing”坐标会：A.自动切换为笛卡尔B.引入非静力项C.降低风速以符合质量守恒D.增加垂直数值扩散答案：C解析：陡坡处CALMET降低风速防止虚假上升，保证地形质量守恒。2.多项选择题（每题3分，共15分）2.1下列哪些变量属于FLEXPART粒子文件header部分必须包含的字段？A.释放高度B.粒子质量C.湍流速度方差D.干湿沉降标记E.粒子年龄答案：ABD解析：header记录释放属性与沉降参数，湍流方差与年龄为运行时变量。2.2在AERMOD中，关于城市选项(UrbanOption)的说法正确的有：A.需额外输入城市人口B.会调用BPIPPRIME建筑下洗C.城市边界层高度公式与乡村不同D.城市选项开启后水平扩散系数增大E.仅适用于PM₁₀答案：ACD解析：城市选项用人口估算热通量，边界层参数化不同，扩散增强；BPIPPRIME为独立模块，不限城市。2.3采用WRFLES求解大涡模拟时，为降低周期边界对污染物回流的影响，可采取：A.增大模拟域B.采用辐射边界C.在侧边界加海绵层D.提高垂直分辨率E.采用nudging松弛答案：ABC解析：增大域、辐射边界、海绵层均可抑制回流；垂直分辨率与nudging不直接解决侧边界回流。2.4对CALPUFF模型进行贝叶斯校准，若选择MCMC采样，需提前设定的先验包括：A.排放源强对数正态均值B.水平扩散系数幂指数C.干湿沉降速度D.气象观测误差标准差E.地形高度答案：ABC解析：MCMC校准扩散参数、源强、沉降速度；地形为固定输入，非随机变量。2.5在利用卫星遥感验证地面NO₂浓度时，需进行的转换步骤有：A.柱浓度→地面浓度需垂直廓线形状因子B.对流层AMF重新计算C.云掩膜筛选D.地表反射率地形校正E.温度压强订正答案：ABCD解析：卫星为柱量，需廓线形状因子；云与地表反射影响AMF，温度压强影响小。3.判断题（每题1分，共10分）3.1在拉格朗日模式中，粒子数足够多时，蒙特卡洛误差与粒子数平方根成反比。答案：正确3.2AERMOD的PRIME下洗算法假设建筑为矩形，且高度与宽度比需小于10。答案：错误（上限为5）3.3CALPUFF的“subgridscalecomplexterrain”选项要求地形格距≤1km。答案：正确3.4WRFChem的RADM2机制包含异戊二烯气相化学。答案：正确3.5采用1km分辨率排放清单时，无需进行空间分配即可直接输入3km网格模型。答案：错误（需保守重映射）3.6FLEXPART的粒子释放可设为持续源或瞬时源，但不可为变化时间序列。答案：错误（支持时变）3.7在稳定边界层中，湍流强度随高度增加而单调减小。答案：错误（低空有最大，随后减小）3.8高斯模型中，地面反射系数α=1代表完全反射。答案：正确3.9若观测站点位于模型网格最外层，则其值可用于验证边界条件。答案：错误（边界为驱动场，非自由验证）3.10采用EnKF同化PM₂.₅时，若观测算子为非线性，需使用扩展卡尔曼滤波。答案：错误（EnKF可处理非线性）4.填空题（每空2分，共20分）4.1在AERMOD中，城市热通量参数H_U的计算公式为H_U=________×人口密度。答案：0.25Wm⁻²person⁻¹4.2若CALPUFF的烟团分裂数NSPLIT设为5，则一个烟团在下一个时间步最多分裂为________个子烟团。答案：254.3WRFLES的次网格湍流闭合采用________模型时，湍流Prandtl数默认0.7。答案：1.5阶TKE4.4当莫宁奥布霍夫长度L=+50m，则边界层稳定度类别为________。答案：E4.5在FLEXPART中，粒子位置输出步长由OUTSTEP控制，单位是________。答案：秒4.6高斯模型地面轴线最大浓度公式中，当有效源高H增加1倍，最大浓度约变为原来的________倍。答案：1/24.7若采用OMI卫星NO₂柱量验证，其本地过轨时间约为________时。答案：13:304.8在CALMET风场插值中，若观测站权重采用反距离平方，则距离10km与5km站权重比为________。答案：1:44.9采用Barnes客观分析时，影响半径R与格距Δx的经验关系为R≈________Δx。答案：1.0–1.54.10在WRFChem的MOSAIC4bin中，粒径边界第3档上限为________μm。答案：2.55.简答题（每题10分，共30分）5.1给出利用CALPUFF模拟沿海SO₂扩散时，考虑海陆风循环的完整设置流程，并指出最易引入误差的三个环节。答案：(1)地形与土地利用：采用30m海岸线数据，更新LU_INDEX区分水陆粗糙度；(2)气象场：WRF三重嵌套最内层1km，选用MYJ边界层，SST每日更新，开启sf_ocean=2；(3)CALMET：读入WRF3D.DAT，开启海气交互模块，设置海陆温差阈值2℃触发环流；(4)排放：港口船舶采用瞬时源，时间分辨率1h，烟囱高度按AIS动态修正；(5)参数：水平扩散采用PGT海陆修正，海洋稳定度默认D类；(6)验证：利用沿海走航监测SO₂，对比模拟日变化相位。易误差环节：①SST更新延迟导致热通量偏差；②海岸线1km误差可使风场转向>30°；③船舶排放时间分配粗糙，峰值错位。5.2解释为何在夜间稳定边界层中，高斯烟羽模型会显著低估地面浓度，并提出两种模型修正思路。答案：夜间稳定层结下，湍流弱，烟羽不易垂直扩散，地面轴线浓度应升高；但标准高斯公式采用PGσ_z，其稳定类σ_z偏小，且未考虑低层逆温“封闭”效应，导致预测值偏低。修正思路：①采用局地扩散系数，将σ_z替换为基于TKE的局地σ_wtf(稳定度)，如Irwin方案；②引入烟羽抬升“封闭盖”假设，在逆温底高度处设置完全反射，重新推导地面浓度公式，可得额外因子2/(1+erf(H/(√2σ_z)))，使地面浓度提高20–80%。5.3说明如何利用无人机获取的垂直NO₂廓线对街区尺度CFD结果进行定量验证，给出误差指标与接受准则。答案：步骤：(1)无人机在距道路中心线50m、100m、200m三条垂线，以20m梯度从30m升至150m悬停，采集30s平均NO₂；(2)将CFD结果沿相同垂线插值，获得模拟廓线；(3)计算归一化平均偏差NMB=∑(M_i−O_i)/∑O_i×100%，目标|NMB|<30%；(4)计算相关系数R，要求>0.6；(5)采用因子2指标：落在0.5≤M/O≤2范围内的点数需>50%；(6)若建筑尾流区出现系统低估，则调整CFD入口湍流强度至10%，并重新模拟。6.计算与推导题（共35分）6.1（10分）某火电厂烟囱有效源高H=120m，源强Q=600gs⁻¹，风速u=3ms⁻¹，稳定度C类，σ_y=0.08x(1+0.0001x)^0.5，σ_z=0.06x(1+0.0015x)^0.5，x单位m。求下风向x=2000m处地面轴线浓度，并给出最大落地距离及浓度。解：x=2000m，σ_y=0.08×2000×(1+0.2)^0.5=160×0.912=146m，σ_z=0.06×2000×(1+3)^0.5=120×0.5=60m，C=(Q/πuσ_yσ_z)exp[−0.5(H/σ_z)²]=(600000/π×3×146×60)exp[−0.5(120/60)²]=729mgm⁻³×exp(−2)=729×0.135=98.5μgm⁻³。求最大落地距离：令dC/dx=0，得数值解x_max≈5500m，σ_z(x_max)=0.06×5500×(1+8.25)^0.5=330×0.328=108m，C_max=600000/(π×3×σ_y×108)×exp(−0.5×1.23²)=600000/(π×3×0.08×5500×0.69×108)×0.46≈600000/(38700)×0.46=7.1mgm⁻³=7100μgm⁻³。6.2（10分）采用CALPUFF烟团模式，设水平扩散参数Kₕ=2×10³m²s⁻¹，时间步长Δt=60s，烟团初始半径R₀=100m，求1h后烟团均方根半径R_rms。解：拉格朗日扩散方程dσ²/dt=2Kₕ，σ²(t)=σ₀²+2Kₕt，σ₀=R₀/√2=70.