S/P模式计算的冷却塔烟气抬升高度敏感性试验

德国VDI3784的S/P模式为三维流体动力学积分模式,其方程主要描述了无穷小体积元素的质量、动量、静态污染物质量浓度及能量的守恒。利用德国模式进行了冷却塔烟气排放不同参数、不同大气条件下烟气抬升高度的敏感性试验。结果表明：在影响烟气抬升高度的3个气象要素（风速、气温和湿度）中,风速和气温的变化对结果影响较大,而湿度影响较小。在D类稳定度,当环境风速从0.1 m/s增加到15.0 m/s时,抬升高度从711.7 m变为38.5 m。随着环境温度的升高,抬升高度明显单调变小;当稳定度为A类,环境温度从10升到40时,烟气抬升最大高度从688.9 m降低到45.1 m,降低了14倍多。而环境湿度的变化,对抬升高度的影响不是很明显。对于E类稳定度和F类稳定度,当环境湿度从20%增加到70%,最大抬升高度分别从115.3 m和84.6 m降到112.9 m和81.7 m,分别降低了3.43%和2.08%。在影响烟气抬升高度的其他3个因素（凉水塔直径、烟气出口速度和混合气体温度）中,混合气体温度的变化对结果影响较大,而凉水塔直径和烟气出口速度的影响较小。在各类稳定度条件下,当出口温度从20变到90时,烟气抬升高度增加1.2—13.3倍;在各类稳定度条件下,当凉水塔直径从30 m变到90 m,烟气抬升高度仅增加0.63—1.40倍;在各类稳定度条件下,当出口速度从2.5 m/s变到8.0 m/s,烟气抬升高度增加了0.24—0.74倍。     

冷却塔烟气抬升高度及污染物浓度变化特征

介绍了冷却塔烟气排放原理及德国VDI 3784的S/P模式,结合实例研究了冷却塔烟气排放不同参数、不同大气条件下烟气抬升高度变化特征,并使用德国VDI 3945的Austal 2000模式对冷却塔和烟囱的排放浓度进行了对比。结果表明：烟气抬升高度主要取决于出口温度和环境温度之差;随着风速的增大,抬升高度明显变小;随着直径的增加,烟气抬升高度也增加;随着出口速度的增加,抬升高度增加的相对较多;随着出口温度的增加,烟气抬升高度在逐渐增加,稳定度决定增加的倍数;VDI 3945模拟结果表明：小时浓度差别相对较小,而年均浓度差别较大,日均值的差别介于小时和年均之间。     

冷却塔烟气抬升高度及污染物浓度变化

介绍了冷却塔烟气排放原理及德国VDI 3784的S/P模式,结合实例研究了冷却塔烟气排放不同参数、不同大气条件下烟气抬升高度变化特征,并使用德国VDI 3945的Austal 2000模式分析了冷却塔与烟囱的排放浓度对比。结果表明:烟气抬升高度主要取决于出口温度气环境温度之差;随着风速的增大,抬升高度明显变小;随着直径的增加,烟气抬升高度也增加;随着出口速度的增加,抬升高度增加的相对较多;随着出口温度的增加,烟气抬升高度在逐渐增加,稳定度决定增加的倍数。VDI 3945模拟结果表明:小时浓度值差别相对较小,而年平均浓度值差别较大,日平均值的差别介于小时与年平均值之间。     

基于CFD模型的内陆核电厂厂区流场模拟

以内陆某核电厂为例,简述了利用流体力学软件STAR-CCM+模拟内陆核电厂厂区流场及大型自然通风冷却塔雾羽扩散的实现原理,介绍如何将SolidWorks2010建立的核电厂厂区模型导入到STAR-CCM+,给出了STAR-CCM+划分网格的过程,说明了边界层划分的基本假定条件和参数。将STAR-CCM+模拟的数据与风洞实验数据进行了比较,结果显示了较好的一致性。结果表明：在离地面5m的高度处,大型自然通风冷却塔背风面形成较大的空腔区,空腔区风速较小,只有1-1.5 m/s,部分区域达到静风;在冷却塔两侧风速相比入口速度增大了1.66倍;在离开地面100 m的高空,冷却塔背风面的空腔区依然比较明显,冷却塔两侧风速相比入口风速,其变化趋于平稳;在沿主导风向的轴线上,冷却塔两侧风的扰动依次加强;单台冷却塔雾羽最大的抬升高度出现在下风向距离3300 m处,最大抬升高度为690 m,4台冷却塔雾羽在下风向距离3300 m的抬升高度约为850 m,是单台冷却塔的1.23倍。     

二等皮托管测风误差分析及偏差模拟计算

皮托管是组成风速检定装置的主要计量标准设备,其测量精度对检定质量具有重要影响.从皮托管测风原理出发,详细介绍了影响皮托管测风误差的相关因素,模拟计算了各因素对风速测量产生的影响,分析了影响皮托管测风误差的主要因子.结果表明：温度和阻塞系数是影响二等皮托管测风误差的主要因素.当环境温度偏差为±8℃时,可引起二等皮托管风速测量误差为干0.44 m/s（v=30 m/s）.当阻塞修正系数偏差±0.02时,可引起二等皮托管风速测量误差±o.6 m/s（v=30 m/s）;皮托管系数、大气压力和湿度经修正后对二等皮托管测风精度影响相对较小.     

估算模式 AERSCREEN的参数敏感性研究

《环境影响评价技术导则大气环境》（HJ2.2-2018）推荐的估算模式AERSCREEN在气象和地形资料的处理以及建筑物下洗等多个方面做了改进。利用估算模式AERSCREEN,针对30 m左右高度的点源,进行了不同排放参数、不同气象条件下最大落地浓度的敏感性试验。结果表明：随着烟气出口流速的增大,地面浓度最大值逐渐减小;随着波文比的变化,地面浓度最大值没有明显的变化;随着地面粗糙度的增大,地面浓度最大值逐渐减小;随着烟气出口温度的增高,地面浓度最大值逐渐减小;当烟气温度为75℃,粗糙度达到1.3 m时,地面浓度达到最小;随着反照率的增大,地面浓度最大值逐渐减小;随着烟囱高度的增大,地面浓度最大值逐渐减小;在各种烟囱高度条件下,随着最高环境温度的增高,地面浓度最大值逐渐增大;而在各种环境温度条件下,随着烟囱高度的增高,地面浓度最大值在逐渐减小;模式中,随着最低环境温度的增高,地面浓度最大值没有变化;但随着最小风速的增大,模拟得到的地面浓度最大值会逐渐减小。     

巴丹吉林沙漠北缘拐子湖近地层湍流能谱特征分析

利用2013年6月巴丹吉林沙漠拐子湖地区流沙下垫面的陆气通量观测资料,计算并分析了该研究区不同大气稳定层结条件下的湍流速度各分量谱和温度谱及湍流的局地各向同性特征。结果表明:该研究区风速各分量的湍流强度随平均风速的增加而逐渐减小,当平均风速2 m/s时,风速各分量的湍流强度逐渐稳定且基本限定在0.5以内。在不同大气稳定度下,湍流速度和温度能谱曲线在惯性副区内逐渐有合并趋势且遵循Kolmogorov提出的-2/3定律,速度各分量谱在高频段均满足各向同性且符合低频限制理论。随着稳定度的增加,风速分量的能谱曲线逐渐降低且向高频端移动,风速分量和温度能谱对应的谱峰长度则逐渐减小。该研究区水平湍流尺度范围为9.0~600 m,垂直湍流谱峰波长为10.79~75 m。该结果介于草地和森林下垫面之间,与塔克拉玛干沙漠腹地的试验结果较为接近。     

钟形地形动力抬升和重力波传播与地形云和降水形成关系研究

地形云和降水过程在区域水循环、水资源、生态环境及气候变化中具有十分重要的作用。本文利用中尺度数值模式WRF 数值模拟试验,以及通过引入表示大气层流速度、层结稳定度和地形特征的关系参数——湿Froude 数(F_w),研究了北京2009 年5 月1 日湿条件不稳定大气层结下,地形云和降水形成过程与地形动力抬升和地形重力波传播之间的关系及形成机理。研究表明,在地形最大高度2 km、半宽10 km 的条件下,层流速度从2.5 m/s 逐步增加到25 m/s 时,对应的湿F_w 数从0.19 增加到1.81。当F_w≤1 时,地形的阻挡起主要作用,由地形抬升形成的地形云主要产生在迎风坡一侧。地形重力波主要产生在迎风坡,并向上游传播,先形成层状云,最后演变为准稳定浅对流波状云。最大降水主要发生在紧靠山顶的迎风坡一侧,但当F_w 很小时,地形云不产生降水。当F_w>1 时,地形抬升形成的云主要发生在山顶附近,而地形重力波主要形成在背风坡,并向下游方向传播,形成准稳定波状云。最大降水主要产生在紧靠山顶的背风坡一侧。另外,在弱湿条件不稳定大气层流下,地形降水主要由地形动力抬升造成的暖云微物理过程产生,地形重力波形成的波状云几乎不产生降水。     

北京北郊冬季大风过程湍流通量演变特征的分析研究

利用中国科学院大气物理研究所325 m气象观测塔1993年12月～1994年1月大气边界层实验资料, 计算分析了大风过境过程中47 m和120 m高度湍流通量演变特征及其影响因子, 以及与风速、 稳定度等参数的关系。结果表明: 大风过程对近地面层的物质能量输送有着重要影响, 大风之前出现短时间动量上传和热量下传; 大风过程中的湍流通量数值明显高于过境后, 水平方向湍流通量数值和能量增加幅度大于垂直方向; 当风速大于临界值5 m/s时, 湍流通量与风速、 湍流动能的相关迅速增大; 湍流谱特征表现为湍流能量的低频部分增加、 湍流谱曲线变宽; 大风能强烈影响近地面层的能量收支。     

塔克拉玛干沙漠腹地起沙阈值计算解析

临界起沙风速是判别风沙活动能否发生的关键指标,其变化受到地表状况及大气环境的综合影响。为了进一步认识临界起沙风速在野外条件下的变化规律,选取塔克拉玛干沙漠腹地塔中作为研究区,在综合考虑地表土壤粒径、土壤湿度、空气密度等因素的基础上,利用经验公式计算了该地区每月的临界起沙风速。得出：（1）塔中地区2 m高度的临界摩擦速度值介于0.24～0.36 m/s,均值为0.31 m/s;（2）塔中地区2 m高度的临界起沙风速值介于3.9～5.9 m/s,均值为5.1 m/s;（3）塔中地区起沙阈值,最高值出现在夏季,次高值出现在冬季,春季最小。     

A model for the height of the internal boundary layer over an area with an irregular coastline

A model for the time and space variation of the internal boundary-layer height over a land area with an irregular coastline is presented. It is based on the analytical model of the boundary-layer height proposed by Gryning and Batchvarova (1990) and Batchvarova and Gryning (1991), The model accounts for the temperature jump and the mean vertical air motion at the top of the internal boundary-layer. Four cases from experiments in Nanticoke and Vancouver are used for model validation. The agreement between the calculated and measured internal boundary layer height at the observational sites is fairly good. The input information for the model consist of wind speed and direction, friction velocity and kinematic heat flux in time and space for the area, and the potential temperature gradient and the mean vertical air motion above the internal boundary layer. For the experiments used in the validation the effect of subsidence is relatively important in the afternoon under low wind speed high pressure conditions, lowering the height of the internal boundary layer by up to 10%, and it is negligible in the morning hours. The effect of the mixing height over the sea is found to be negligible.     

太原雾天能见度预报

利用中尺度数值预报模式MM5对山西省2009年发生的几场典型雾个例进行数值模拟。结果表明：模拟2 m温度比观测值偏低约2 ℃,相对湿度模拟结果比观测值偏大约15 %,10 m的模拟风速比观测的偏大0-2 m·s^-1。山西省雾的预报指标为20 m液态水含量大于等于0.13 g·kg^-1而小于0.60g·kg^-1、20-1500 m高度大气层存在逆温层、地面风速小于4 m·s^-1。利用太原测站日平均能见度、日平均相对湿度以及空气污染指数进行拟合建立太原能见度预报模型,并利用实测资料订正MM5、CAPPS模式预报误差,给出订正后的能见度预报方程并以两次实例对区域及太原雾天能见度预报表明该能见度预报模型有一定的适用性。     

基于数值模式的月尺度近地层气象要素预报技术研究

利用WRF模式对美国NCEP发布的CFS气候预测业务产品在中国区域内进行动力降尺度预报,可得到预报时效为45天的逐6小时、30 km分辨率基础气象要素预测产品。再利用全国气象站观测资料和3个风电场70 m高度风速、温度观测资料对2015年冬季预测结果进行检验评估和分析,最后通过线性方法对地面要素预测结果和70 m高度风速、温度预测结果进行统计订正。结果表明:（1）2 m温度和相对湿度的全国预报平均绝对误差分别为4.71 ℃和18.81%,在华东、华中和华南地区误差较小;（2）10 m风速预报平均绝对误差为2.42 m/s,在东北、华北和西北地区误差较小;（3）线性订正后,2 m气温、相对湿度和10 m风速的预报绝对误差分别减小1.05 ℃、5.29%和1.47 m/s,并且订正后误差随时间变化更平稳;（4）订正后70 m高度风速和温度的预报绝对误差均减小,风速平均误差减小最大可达1.29 m/s（B塔）,气温平均绝对误差减小最大可达3 ℃（C塔）。研究结果表明,基于CFS产品和WRF模式的、与月尺度风电预报关系密切的气象要素预报性能较好,未来可将该方法尝试于风电场的月尺度功率预测产品研发。      

Nocturnal boundary layer height prediction from surface routine meteorological data

Summary A number of well known diagnostic equations for the determination of the height,h, of the nocturnal boundary layer. with minimum data requirements of at most surface wind speed, air temperature and total cloud cover, have been tested as to their effectiveness. The computed values have been compared with direct estimation ofh, from temperature or wind profiles of rawinsonde ascents available at 00Z (02h LST). The comparison between computed and observed values shows that best agreement is found when the nocturnal boundary layer height is determined through wind profiles. The ratio of the computed to the observed values reveals a strong dependence on stability, resulting in overestimation by the models for very low stability and underestimation for strong stability. The simple expressions involving the wind speed rather than other stability parameters resulted in a better overall fit to the observed values. A simple prognostic model is shown to provide the best estimates of the NBL height compared to both wind and temperature profile definition.With 5 Figures     

Seasonal variations of spectra of wind speed and air temperature in the mesoscale frequency range

Seasonal variations of the spectra of wind speed and air temperature in the mesoscale frequency range from 1.3 × 10^-4 to 1.5 × 10^-3 Hz (10 min to 2 h periods) have been studied through observations over land for one year. Spectrographs [time series of isopleths of spectral densities, f · S(f) vs f] of wind speed and air temperature contain occasional peaks that are attributed to short-lived mesoscale atmospheric activity with narrow frequency bands. Significant spectral peaks of wind speed were found in 19% of the total observations in winter, and in 15–16% in the other seasons; for air temperature, they occured in 12% of observations in autumn, and in 16–19% in the other seasons. The peaks most often occurred in the period range from 30 min to 1 h; most had durations less than 24 h. Mesoscale fluctuations of wind speed and air temperature were highly correlated, and in most cases, phase differences were 90–180 ° with air temperature leading wind speed. Significant spectral peaks of wind speed often occurred during northerly seasonal cold winds in winter, and accompanied tropical and/or mid-latitude cyclones in the other seasons. When the peaks occurred, wind speed was usually relatively high and the atmospheric surface layer was unstable.     

Sensitivity analysis of planetary boundary layer schemes using the WRF model in Northern Colombia during 2016 dry season

This study conducted meteorological simulations in northern Colombia by analyzing different planetary boundary layer (PBL) schemes available in the numerical Weather Research and Forecasting (WRF) model. The study area included three nested domains with horizontal resolutions of 18 km, 6 km, and 2 km, with 38 vertical levels. The evolution and structure of the PBL were analyzed during the driest months (March, April, and May 2016) and in regions with the highest particulate matter concentrations. Sensitivity analysis of the WRF model was performed with two local and two non-local PBL schemes. The results were validated using observations of the surface air temperature, relative humidity, and surface wind speed collected from three meteorological stations in the area. The PBL heights were experimentally determined using radiosonde data provided by a station located in the center of the study area. Variations in PBL heights were estimated using linear regression methods and minimization of statistical errors for the bulk Richardson number, as well as analysis of vertical temperature and wind profiles. The WRF model reliably reproduced the daily values and diurnal cycles of temperature, relative humidity, and wind speed within the PBL and accounted for the influence of topography and sea breezes. Horizontal heat advection dominates the upwelling of air masses when sea breezes are active. The onshore wind direction starts to change from east to northwest, implying a decay in the land breeze regime. All schemes overestimate the mixing height and tend to underestimate surface air temperature values at night. All show wetter conditions and underestimate wind speed. Although the non-local Yonsei University (YSU) scheme shows the best performance, it also shows the largest sources of errors when determining the behavior of the surface layer during stable conditions. Relative humidity and wind speed estimates provided by the local Mellor‐Yamada‐Nakanishi‐Niino (MYNN) scheme were closer to those recorded at the meteorological stations.     

Wind regime in the lower atmosphere over Moscow from the long-term acoustic sounding data

Presented are the results of the sounding of the lower atmospheric 500-meter layer for the period of 2004–2012 carried out at the Meteorological Observatory of the Moscow State University (MSU) with the MODOS Doppler acoustic radar (sodar) produced by METEK (Germany). Discussed is the methodological basis of the sodar wind data analysis. It is demonstrated that in the air layer up to 200 m the maximum values in the annual course of the wind speed are observed more often in autumn and winter, and the minimum values, in summer; this is associated with the fact that during the cold period of the year Moscow is often located in the zones of intense gradient currents. The diurnal course of the wind speed is characterized by the daytime maximum and night-time minimum in the layer up to 40–60 m from the surface; it is poorly pronounced and characterized by the minimum in the morning in the layer of 80–120 m; and the daytime minimum and night-time maximum are observed above 140–160 m. The layer from 80 to 120 m approximately corresponds to the height of the wind rotation. The amplitude of diurnal variations of the wind speed increases from 0.3 m/s at the height of 7 m and 0.6 m/s at the height of 15 m, to 4.5 m/s at the height of 400 m; however, its secondary minimum (0.5 m/s) associated with the rotation height is registered at the altitude of 80 m. The statistical relationship between the wind speed and surface air temperature is direct during the cold season, inverse during the warm season, and is absent in April and October. The average maximum wind speed over Moscow for ten minutes in the layer up to 500 m from the surface reaches 30–35 m/s in some cases if two conditions concur: the capital is located on the periphery of vast pressure formations (usually of deep cyclones) and the local low-level jet stream is present in the wind profile.