% 功能:烟囱排放污染物烟羽扩散模型 % 邯郸纬度 = 36度 % 太阳高度角 = 90-abs(当地纬度-太阳直射点纬度) % 太阳直射点纬度 = 春分0度——夏至+23.5度[北回归线]——秋分0度——冬至-23.5度[南回归线] % 太阳直射点纬度每天移动的度数 = 23.5*4/当年的天数 % 计算当天的太阳高度角 % 当天的太阳高度角 = 90-abs(36-(23.5-(当前日期-夏至日期)*太阳直射点纬度每天移动的度数)) % 参考《大气污染物浓度估算方法》 % f(太阳高度角、云量)=太阳辐射等级【二维插值或者模糊算法】 % f(太阳辐射等级、风速)=大气稳定度【二维插值或者模糊算法】 % 风速:烟囱出口处的平均风速【区别于地面风速、如果是地面风速{比如在10米测得风速v10、需要校正、校正公式:v10*(烟囱高度/10)*m}、m为某稳定度下的风指数[查表]】 % 参考《第五章 大气扩散参数ppt课件》 % 扩散系数:大气稳定度对应的参数[查表] % 抬升高度 % 参考《大气污染物浓度估算方法》结合《镇江市大气稳定度及抬升高度计算方法的研究》
1
功能:烟囱排放污染物烟羽扩散模型 邯郸纬度 = 36度 太阳高度角 = 90-abs(当地纬度-太阳直射点纬度) 太阳直射点纬度 = 春分0度——夏至+23.5度[北回归线]——秋分0度——冬至-23.5度[南回归线] 太阳直射点纬度每天移动的度数 = 23.5*4/当年的天数 计算当天的太阳高度角 当天的太阳高度角 = 90-abs(36-(23.5-(当前日期-夏至日期)*太阳直射点纬度每天移动的度数)) f(太阳高度角、云量)=太阳辐射等级【二维插值或者模糊算法】 f(太阳辐射等级、风速)=大气稳定度【二维插值或者模糊算法】 风速:烟囱出口处的平均风速【区别于地面风速、如果是地面风速{比如在10米测得风速v10、需要校正、校正公式:v10*(烟囱高度/10)*m}、m为某稳定度下的风指数[查表]】 扩散系数:大气稳定度对应的参数[查表]
1