基于WRF和XFlow耦合的多精细度融合污染物扩散分析方法与流程

本发明属于绿色建筑风环境评估技术领域，具体涉及一种基于wrf和xflow耦合的多精细度融合污染物扩散分析方法。

背景技术：

近年来，随着城市化和工业化不断发展，大气污染引起的公共环境越来越受到关注。影响公共环境安全的污染物主要包括工业废气以及医疗废气，这些污染物的排放及处理不当不仅会造成环境污染问题，甚至会直接威胁到人民的身心健康。

基于cfd(computationalfluiddynamics，即计算流体动力学)的环境分析能够对污染物扩散进行定量分析，为建筑物规划和通风系统的合理布置提供参考，避免污染物长期滞留，最大限度地消除污染物排放对周围环境和周边人群造成的不利影响。然而传统的cfd环境分析存在着大气边界条件和污染物源边界假定不明确、建模难度大、难以考虑真实天气等局限性。

基于中尺度气象预报模式wrf的大气数值模拟方法，能够对城市区域进行模拟分析，给出城市区域的三维气象风场信息。然而中尺度气象预报模式wrf的空间分辨率一般在1km量级以上，未能精细模拟城市近地面和建筑物区域的风场特性。

相比现有的污染物扩散分析方法，本文提出的方法能够同时实现城市级和建筑级的多精细度融合污染物扩散模拟，为化工厂、防疫医院等存在污染物排放的企业的科学选址提供精细化的数据支撑。

技术实现要素：

本发明的目的在于，提供一种基于wrf和xflow耦合的多精细度融合污染物扩散分析方法，既能考虑气象因素的城市级污染物扩散过程，还能够考虑建筑物室内微尺度污染物扩散过程。

本发明提供的技术方案如下：

一种基于wrf和xflow耦合的多精细度融合污染物扩散分析方法，包括以下步骤：

s1、建立建筑物污染物扩散情景库；

s2、利用xflow软件建立建筑物室内cfd计算模型，进行特定情景下的微尺度污染物扩散分析，得到进入建筑物通风系统的污染物浓度；

s3、根据通风系统的净化能力，确定建筑物最终排出的污染物浓度，并作为建筑物微尺度污染物浓度代表值；

s4、建立城市气象情景库；

s5、利用中尺度气象数值模式wrf建立区域多精细度融合网格嵌套模型，进行中尺度三维气象驱动场模拟计算分析；

s6、基于wrf计算结果，判断主导风向；

s7、导入城市几何模型，并根据主导风向确定计算域的大小和方向，建立城市级三维污染物扩散分析模型；

s8、建立多精细度耦合：在城市级三维污染物扩散分析模型中，将建筑物微尺度污染物浓度代表值升尺度为城市级模型的污染源入口边界条件，将wrf中尺度三维气象驱动场降尺度为城市级模型的初始条件和大气边界条件；

s9、进行城市污染物扩散多精细度融合计算，并根据计算结果进行污染物定量评估。

优选地，步骤s1具体为：根据建筑物的类别和使用功能，收集污染物扩散情景，建立建筑物污染物扩散情景库。

优选地，步骤s1包括以下子步骤：

s11、按建筑类别搜集以往发生的污染物扩散案例，按照情景要素对案例进行特征提取，进行情景划分、分析和归类，建立建筑物污染物扩散情景库；

s12、根据污染物扩散情景库选取特定的污染物扩散分析工况，确定污染源的空间分布和相对浓度大小。

优选地，步骤s2包括以下子步骤：

s21、打开xflow软件并创建新的建筑物室内污染物扩散分析工程项目，将xflow计算内核设置为3d，设置流动模型为singlephase，分析类型为internal，开启scalartransport模式，指定diffusivity系数为1e-5m²*s^-1；

s22、导入建筑物几何模型和通风系统几何模型，对计算域进行格子划分；

s23、设定建筑物壁面模型、送风口的入口边界条件以及排风口的出口边界条件，根据污染源空间分布和相对浓度大小，确定质量入口边界；

s24、提交xflow计算任务，对计算结果进行后处理，给出进入通风系统的污染物浓度。

优选地，步骤s3包括以下子步骤：

s31、通过对通风系统进行调研、实测或者数值仿真的方式，给出通风系统能够稀释污染物的倍数；

s32、结合进入通风系统的污染物浓度以及通风系统的稀释倍数，得到最终从建筑物排出的污染物浓度代表值。

优选地，步骤s4具体为：收集城市的气象数据，建立城市气象情景库。

优选地，步骤s4具体包括以下子步骤：

s41、通过实测或者访问气象数据库的方式获取目标城市的气象数据；

s42、通过气象数据特征提取，进行分析和归类，归结为目标城市的气象情景库。

优选地，步骤s5包括以下子步骤：

s51、选取目标日期，下载气象再分析和gis地形地貌数据，设置四层嵌套网格计算域，启动单向网格嵌套模拟，水平网格精度依次为27、9、3和1km，其中目标建筑群位于最内层网格中心区域；

s52、设置多种描述大气物理变化过程的参数方案，基于wrf最内层网格计算域进行近地面台风风场模拟。

优选地，步骤s6包括以下子步骤：

s61、对wrf气象数据进行脚本批量后处理，计算得到最内层网格嵌套域的三维风速和湍流场信息；

s62、根据wrf气象风场数据，判断主导风向。

优选地，步骤s7包括以下子步骤：

s71、打开xflow软件并创建新的工程项目，将xflow计算内核设置为3d，设置流动模型为singlephase，分析类型为external。开启scalartransport模式，指定diffusivity系数为1e-5m²*s^-1；

s72、导入城市几何模型，根据wrf主导风向确定计算域大小和方向，对计算域进行格子划分，建立城市级三维污染物扩散分析模型。

优选地，步骤s8包括以下子步骤：

s81、在城市级三维污染物扩散分析模型中，选取含污染源的建筑物，根据其排出的污染物浓度设定污染源入口边界条件；

s82、导入wrf气象驱动场，进行网格耦合转换为城市级模型的初始条件，拟合出计算域边界处的大气边界条件，并据此设定计算域侧面的入口边界条件和出口边界条件。

优选地，步骤s9包括以下子步骤：

s91、提交城市级多精细度融合xflow污染物扩散分析计算任务；

s92、对目标区域的污染物浓度值进行评估。

本发明的有益效果为：

本发明的基于wrf和xflow耦合的多精细度融合污染物扩散分析方法，既能考虑气象因素的城市级污染物扩散过程，还能够考虑建筑物室内微尺度污染物扩散过程，实现城市级和建筑级的多精细度融合污染物扩散模拟，为化工厂、防疫医院等存在污染物排放的企业的科学选址提供数据支撑。

附图说明

图1为基于wrf和xflow耦合的多精细度融合污染物扩散分析方法流程图。

图2为本发明实施例中的某防疫医院功能分区示意图。

图3为本发明实施例中的负压病房建筑模型示意图。

图4为本发明实施例中的负压病房含病毒空气浓度三维空间分布云图。

图5为本发明实施例中的多重网格嵌套的wrf风向矢量和风速云图。

图6为本发明实施例中的城市级xflow污染物扩散模型图。

图7为本发明实施例中的隔离病房和icu排风口的污染浓度示意图。

图8为本发明实施例中的多精细度融合分析模型所给出的含病毒空气浓度三维空间分布云图。

具体实施方式

下面将结合附图对本发明作进一步的说明：

为了使本发明的目的、技术方案和优点更加清楚明白，以下结合附图及具体实施案例，对本发明作进一步的说明，但本发明的保护范围不限于下述的实施例。

本发明涉及一种研究污染物扩散的多精细度融合的计算流体力学模拟方法。该方法包括以下步骤：1)建立建筑物污染物扩散情景库。2)应用xflow软件实施特定情景下的微尺度污染物扩散分析，得到进入建筑物排风系统的污染物浓度；3)确定建筑物最终排出的污染物浓度，作为建筑物的污染物浓度代表值；4)建立城市气象情景库；5)利用中尺度气象数值模式wrf分析计算城市建筑群三维气象驱动场；6)处理wrf计算得到的风速和湍流场信息，判断主导风向；7)导入城市几何模型，并根据主导风向确定计算域大小和方向；8)建立多精细度耦合：将建筑物微尺度污染物代表值升尺度为城市级模型的污染源入口边界条件，将wrf中尺度气象驱动场降尺度为城市级模型的初始条件和大气边界条件；9)完成城市污染物扩散多精细度融合计算，根据计算结果对目标区域进行污染物定量评估。本发明既能考虑气象因素的城市级污染物扩散过程，还能够考虑建筑物室内微尺度污染物扩散过程。

本发明实施例的基于wrf和xflow耦合方式的多精细度融合污染物扩散分析方法，如图1所示，包括以下步骤：

针对武汉地区冬季气候的某防疫医院污染物扩散分析为例，这里的污染源为含病毒的空气，具体实施步骤如下：

s1、根据不同建筑物的类别和使用功能，收集可能发生的污染物扩散情景，建立建筑物污染物扩散情景库。具体包括以下子步骤：

s11、该防疫医院是由标准模块装配而成的拼接式建筑，病房区域的负压隔离病房单元采用规格统一的集装箱。如图2所示，医院病房主要分为负压隔离病房和icu病房，根据这两种病房的使用特点，分别进行含病毒空气扩散案例收集，进行情景特征提取。

s12、以负压隔离病房为例，病房内一般包含两位感染病毒的病人，其内部构造如图3所示。负压隔离病房最典型的含病毒空气扩散情景为：含病毒空气从病的口鼻中呼出，在病房内部扩散，然后进入通风系统，经过滤装置净化后排出医院。

s2、应用xflow软件建立建筑物室内cfd计算模型，实施特定情景下的微尺度污染物扩散分析，得到进入建筑物排风系统的污染物浓度。具体包括以下子步骤：

s21、打开xflow软件并创建室内污染物扩散分析工程项目，将xflow计算内核设置为3d，设置流动模型为singlephase，分析类型为internal。开启scalartransport模式，指定diffusivity系数为1e-5m²*s^-1。

s22、导入负压隔离病房、通风系统、病人面部几何模型，综合考虑几何体的复杂程度和计算消耗后，其格子划分尺寸分别设为0.125m、0.03125m和0.0039m，对计算域进行格子划分。

s23、将病房内的墙、地面、床面等几何部件设置为enhancedwall-function壁面模型，根据通风设备的参数分别设置送风口的速度边界条件和出口风口的压力边界条件。由于病人呼出空气含病毒的浓度难以测定，因此这里采用相对浓度来评估通风系统对病人呼出含病毒空气的稀释作用。将病人面部几何体设定为质量入口边界，病人呼出的含病毒的空气相对浓度为1。

s24、提交xflow计算任务，得到如图4所示的含病毒空气扩散过程及其浓度场，对计算结果进行后处理，得到从病房进入通风系统的含病毒空气浓度。

s3、根据排风系统是否具备净化能力，确定建筑物最终排出的污染物浓度，作为建筑物的污染物浓度代表值。具体包括以下子步骤：

s31、根据医院负压隔离病房通风系统采用的过滤装置，通过仿真分析和实测得到通风系统能够稀释含病毒空气的浓度。

s32、结合进入通风系统的含病毒空气的浓度以及通风系统的稀释倍数，得到最终从病房排出的含病毒空气浓度代表值。通过类似的方法，也可以确定其他病房排出的含病毒空气浓度代表值。

s4、收集城市的气象数据，建立城市气象情景库。具体包括以下子步骤：

s41、通过访问气象数据库的方式获取武汉地区的气象数据。

s42、分析武汉地区的季节性气候特征，选出具有代表性的气候类型，归结为武汉的气象情景库。针对本案例的医院，选取冬季气候进行气象分析。

s5、利用中尺度气象数值模式wrf建立区域多精细度融合网格嵌套模型，进行三维气象驱动场模拟计算分析。具体包括以下子步骤：

s51、选取目标日期，下载气象再分析和gis地形地貌数据，设置多层嵌套网格计算域。网格计算域的水平网格精度分别为27km、9km、3km和1km，水平网格数量分别为152×116、142×106、145×109和130×94，垂直方向采用拉伸网格，对近地面高度进行局部加密，一共划分40层，设置最外层网格计算域的时间积分步长为120s，设置时间积分步长嵌套比为3。

s52、为考虑大气环境场的复杂相互作用过程，模拟算例激活了大气微物理方案、边界层方案、积云方案、辐射方案和陆面方案等参数化方案。其中，微物理方案采用wsm6高精度模拟方案，辐射方案采用rrtm长波方案和dudhia短波方案，针对最外两层网格计算域激活kain-fritsch积云方案，采用hong-shin边界层方案，用于求解大气边界层湍流结构。

s6、处理wrf计算得到的风速和湍流场信息，判断主导风向。具体包括以下子步骤：

s61、对wrf气象数据进行后处理，得到最内层网格嵌套域的三维风速和湍流场信息。wrf计算得到的多重网格嵌套的风向和风速矢量图如图5所示。

s62、对wrf气象风向矢量图进行分析，给出冬季主导风向。

s7、导入城市几何模型，并根据主导风向确定计算域的大小和方向，建立城市级三维污染物扩散分析模型。具体包括以下子步骤：

s71、打开xflow软件并创建新的工程项目，将xflow计算内核设置为3d，设置流动模型为singlephase，分析类型为external。开启scalartransport模式，指定diffusivity系数为1e-5m²*s^-1。

s72、导入城市几何模型，根据wrf主导风向和城市几何模型的大小，确定计算域大小为4000m*10000m*300m，远场计算最大格子尺寸为16m，计算域的格子划分形式如图6所示。

s8、建立多精细度融合耦合：在xflow城市级模型中，将建筑物微尺度污染物代表值升尺度为城市级模型的污染源入口边界条件，将wrf中尺度气象驱动场降尺度为城市级模型的初始条件和大气边界条件。具体包括以下子步骤：

s81、根据各病房排出的污染物浓度的代表值，分别选取各病房排风口几何体，设定相应的污染源入口边界条件。各病房排出的污染物浓度代表值如图7所示。

s82、导入wrf气象驱动场，通过网格耦合转换为xflow模型的初始条件，同时拟合出计算域边界处的大气边界条件，并据此设定计算域侧面的入口边界条件和出口边界条件。

s9、完成城市污染物扩散多精细度融合计算，根据计算结果对目标区域进行污染物定量评估。具体包括以下子步骤：

s91、提交城市级多精细度融合xflow污染物扩散分析计算任务。

s92、得到如图8所示的含病毒空气浓度云图，对目标区域的污染物浓度值进行评估，给出相对于病房内病人呼出空气的稀释倍数，从而结合相关研究对医院的防控提出建议。

本领域的技术人员容易理解，以上仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。