王卫国 蒋维楣
(南京大学大气科学系，南京 210093)
余 兴
(陕西省气象局，西安 210015)
摘要 利用边界层预报模式提供的非定常的风场和湍流场，采用实用的考虑风切变的烟团模式模拟了海岸复杂地形下(深圳地区)高架点源(电厂)排放SO2的地面浓度变化特征.联接模式细致地考虑了地面反射、混合层穿透和反射.结果反映了风场湍流场的非均匀非定常对浓度分布有较大影响.在海陆风转换期内，地面浓度分布型式较乱，只有在海风盛行时，浓度分布与简单的高斯模型结果接近.浓度变化呈现明显的周期性，最大浓度出现在陆风向海风转换期内，日均最大浓度分布在陆上.这些特征在常用的高斯模式中是无法反映的.
关键词 海岸复杂地形；联接模式；烟团模式；空气污染；深圳地区.
NUMERICAL SIMULATION OF ATMOSPHEERIC TRANSPORT AND DIFFUSION OVER COSTAL COMPLEX TERRAIN IN SHENZHEN REGIONS
Wang Weiguo Jiang Weimei
(Department of Atmospheric Sciences,Nanjing University,Nanjing 210093)
Yu Xing
(Shaanxi Wenther Burean,Xi′an 710015)
ABSTRACT The puff trajectory model was used to simulate the transport and dispersion of SO2 released from elevated point source.Time-dependent three-dimensional distributions of mean and turbulence variables needed in the puff diffusion model were provided by nonhydrostatic PBL model with energy closure scheme.The reflection on both of the ground and the top of mixlayer,penetration on the top of mix-layer were considered in detail.The results indicated that time-variation and nonuniformity of airflow and turbulence obviously affected the concentration distributions.Especially during the period of sea-land breeze transition,the distribution pattern of ground level concentration is irregular.It is close to that obtained from Gauss model only when sea-breeze fully developed.In a day,the periodic variation of ground-level concentration was found,its maximum value appeared when land-breeze turn to sea-breeze.The daily mean of concentration was found on land.The analysis revealed that the above-mentioned characteristics of concentration distribution over coastal complex terrain could not be obtained from simple Gauss model.
Keywords coastal complex terrain,linked model,puff model,air pollution,Shenzhen regions.
1 引言
用边界层预报模式和扩散模式联接对沿海地区复杂地形条件下特殊气流系统中的空气质量进行模拟预测，是近一二十年来才进行的工作［1—3］. 这些成功的工作说明了中尺度大气扩散模拟系统在研究空气质量问题中的潜力.对复杂地形地区的空气质量模拟需要考虑两方面，一是要有较好地反映复杂地形特征的非均匀、非定场的风场和湍流场，另一是适合于在复杂气流中使用的扩散模式.本文针对深圳地区特殊的地形，以及污染源(电厂)附近还有大城市(香港、深圳)和大亚湾核电站、建立了一个精细的中尺度大气模拟系统来模拟该地区特殊气流系统(海陆风)下SO2浓度分布特征.三维非静力能量闭合的边界层预报模式［4］提供具有时、空变化的风场和湍流场，扩散部分采用考虑风切变的7层烟团轨迹模式［5，6］，该模式主要考虑了风切变和湍流场的不均匀性.本文简单介绍了模式结构，对深圳地区进行了模拟，并对模拟结果进行了讨论.
2 模式介绍
2.1 三维边界层非静力模式
采用三维非静力的边界层数值模式，闭合方案用工程上实用的能量(E-ε)闭合，模式方程组由运动预报方程(u,v,w)，位温预报方程(θ)，连续方程以及湍流动能预报方程和湍能耗散率预报等组成.初始条件由模拟域内的实测点资料输入内插模式输出而得.此模式能较好地模拟海岸复杂地形下的流场和湍流场［4］.图1是本文大气扩散模拟系统的框图.
2.2 烟团轨迹模式
实际处理中将烟团看成一个椭球体，水平半轴长3σy，垂直半轴长1.5σz，每个烟团在垂直方向上以7个点作为7个子烟团的中心，各点垂直方向间距为0.5σz，则子烟团中心点位置为z0+0.5(i-4)σz,(i为子烟团代号，i=1，2，…，7，z0为大烟团中心点位置，7个子烟团由椭球积分得到质量分别为：Q1=Q7=0.02QT，Q2=Q6=0.136QT，Q3=Q5=0.220QT,Q4=0.248QT，QT为烟团总物质量.各子烟团具有的性状，包括移速、扩散率等由它中心点处具有的量值决定.烟团扩散参数由它经历的行程和大气湍流扩散能力决定，轨迹由中心所在网格上的风场确定.假设某个烟团被分成7个子烟团，以i表示(i=1,2，…，7)，设t时刻第i个子烟团中心点坐标(xi,yi,zi,t)，该点处气流速度分量为u(xi,yi,zi,t),v(xi,yi,zi,t),w(xi,yi,zi,t),则在t+Δt时刻该子烟团位置为：
xi(t+Δt)=xi(t)+u(xi,yi,zi,t)Δt
yi(t+Δt)=yi(t)+v(xi,yi,zi,t)Δt
zi(t+Δt)=z4(t)+w(x4,y4,z4,t)Δt+wp(x4,y4,z4,t)Δt
+0.5(i-4)σz+β*w(xi,yi,zi,t)Δt
式中wp为烟气抬升引起的烟团垂直运动速率［7］，一般在烟团运动10—15倍的源高后为零(即到达终极抬升).β为系数，当i=4时，β=0，否则为1
|
图1 大气扩散模拟系统流程框图 任一时刻t的水平和垂直扩散参数σy(t),σz(t)由Taylor的扩散统计理论求得［8］，取指数形式的自相关系数，拉格朗日时间尺度TLV,TLw参照Hanna［9］给出的，水平和垂直方向上的速度脉动方差σ2V(= |
|
图2 模拟域示意图(等值线单位：m) 图3给出了海风阶段(15：00)、海风向陆风转换阶段(01：00)、陆风阶段(04：00)、陆风向海风转换阶段(11：00)的地面SO2浓度模拟预测结果.图3a是白天海风盛行期形成的浓度分布(15：00)的状况，此阶段气流由海上吹向陆地，SO2散布在陆地上，地面浓度极大值为72μg/m3，出现在源北陆上10km,由于海风风速较大，一部分烟团移出模拟域，日间对流混合剧烈，地面浓度出现大于25μg/m3的范围较大(远)，地面浓度极值不大，且出现距离较远.在离源下风向约30km的陆上，还出现次大值37μg/m3，这可能是下沉气流和烟团从TIBL以上进入内陆混合层造成的.图3b是海风向陆风转换期间造成的地面浓度(01：00)分布，此阶段气流走向复杂，且风速较小.从图上看到，浓度等值线分布比较乱，高值分布已由白天的陆上移到海上，最大值出现在离源6km的海面上，为231μg/m3，极值比海风阶段大，出现距离也近，这是因为此阶段风速小和层结较稳定造成的.陆上还有一些残留的烟团影响，但受气流的影响，也逐渐向海上移动，陆上极值仅为9.7μg/m3.图3c是夜间陆风阶段造成的浓度(04：00)分布.此阶段气流从陆上吹向海面，风速比过渡阶段要大些，且层结稳定，与过渡期间比，浓度极值减少，为162μg/m3，出现的距离增大，为源下风向10km的海面上，这时陆上烟团已基本移到海面上，陆上浓度很小.图3d是白天早上陆风向海风过渡阶段(或海风初始发展阶段) 造成的地面浓度(11：00)分布.气流已经开始转为由海上向陆上吹，但此时速度还不大，陆上层结基本上处于中性或弱不稳定.SO2浓度大值等值线已逐渐移到陆上，海面上还有由于夜间陆风输送的烟团的影响，海面上浓度还有较大的分布，不过这些烟团已开始向陆上移动，随后的时间里，海面上的SO2分布范围减小，而陆上增大.此过渡阶段极大值为635μg/m3，出现在离源6km的陆上.属于高污染期. |
|
图3 不同时间地面SO2浓度等值线(μg/m3) |
|
图4 地面日平均浓度等值线(μg/m3)
|
|
图5 有效源高和逆温层厚度日变化比较 直线烟流高斯模型一般不可能模拟出复杂气流系统中浓度分布的变化过程和非定常、非均匀的影响.不过，在海风盛行期，采用高斯模式估算结果与联接模式结果相差不大，如15：00，前者估算地面浓度(扩散参数为σy=0.32x0.83；σz=0.39x0.80，引自前面脚注)极大值90μg/m3，出现距离为下风向8km,与后者接近，且联接模式浓度等值线与高斯模式分布类似.这是因为海风盛行期，陆上气流基本一致，近似符合高斯模式的条件.但在气流复杂、物理量随时间变化较大时，高斯模式就不适用了.一些清除过程如云、干湿沉降等在本文中没有能够仔细考虑，所以模拟所得的浓度是保守的(偏大)，边界层参数特别是湍流量对浓度结果影响较大，根据本文的结果分析看，与一般的实测结果的量级相比，估计误差在两倍以内. 参考文献 1 Bornstein R,Klotz S et al.Application of linked three-dimensional PBL and dispersion models to New York city.NATO proc 15th meeting on air pollution modeling and its application,1986:543 1996-04-29收到原稿 1997-02-28收到修改稿 |
+
),σ2W(
)，由边界层预报模式输出量提供脉动速度方差.考虑到烟气抬升对扩散参数的贡献，扩散参数表达式加上(wp Δ t/3.5)2.于洪彬等［10］对水陆交界地区的湍流特征进行了观测研究，得到了三维风速脉动方差在不同稳定度下的比值，本文在这些观测和边界层理论推求得到的湍流脉动速度之间的关系的基础上，作为尝试将湍流动能按下列比例分解(H为混合层高度)：
返回列表
