[分享]新冠疫情临时医院排风的环境影响快速模拟

 摘要：2019年末爆发的新型冠状病毒肺炎疫情对现有医疗系统造成了巨大压力，各地开始建设新型冠状病毒肺炎的集中治疗临时医院。临时医院排风口排出的气体可能对临时医院新风口或者周围环境造成污染。由于临时医院的设计和建造周期只有6-10天，因此排风环境影响分析一般需要在几小时内完成。基于上述背景，本研究提出了一个临时医院排风环境影响的快速模拟方法。本方法以开源流体力学计算软件FDS为基础，实现了临时医院建筑的快速建模、基于云计算平台的分布式计算、以及有害空气流动的监测和可视化，为临时医院设计阶段的快速分析提供了专门工具。并以武汉雷神山医院为案例，说明了本方法在此次肺炎疫情防疫工作中的应用价值。

关键词：新型冠状病毒肺炎；污染扩散模拟；FDS；快速建模；云计算

Key words:  novel coronavirus pneumonia; pollution dissipation simulation; FDS; fast modeling; cloud computing

        为应对现有医疗救治条件不足的问题，中国各地开始建设新型冠状病毒肺炎患者的集中治疗临时医院。譬如，武汉新建火神山和雷神山两座临时医院专门收治新型冠状病毒确诊患者；北京改造当年用于SARS治疗的小汤山医院专门收治新型冠状病毒患者。这些医院为了防止病毒扩散引起感染，在室内暖通设计时采取了负压病房空气系统，利用专门的排风系统向室外高空排出病房内受污染的空气，以免向医院内其他未污染区域扩散。虽然排风系统针对含病毒空气采取了很好的过滤措施，但是其通过排风口排放仍然存在一定传染风险。有研究表明，当含SARS的病毒空气被稀释到10000倍以上时，就不再具有传染性[1]。很多学者和设计单位建议了室外排风的设计要求[2-4]。在以上成果的基础上，如果可以进一步结合具体医院的实际情况，模拟室外含病毒空气的扩散机制，则可以定量评估临时医院排放出的含病毒空气的污染风险，为排放物的二次污染防控提供更加定量具体的参考和依据。值得注意的是，由于临时医院的设计和建造周期只有6-10天，留给排风的环境影响分析时间一般只有几个小时。因此，对分析方法的准确性和高效性都提出了很高的要求。

        目前，可以采用四种方法研究建筑物周围的污染物扩散现象，包括 (1) 实地测量、(2) 物理试验、(3) 半经验模型、(4) 数值模拟。实地测量只适用于记录已建成区域的污染物浓度，无法为临时医院的设计提供参考。物理试验需要制作所关注建筑或区域的缩尺模型和进行风洞试验，并且需要满足相似性准则，需要大量的时间[5]，无法满足防疫应急需要。因此，物理试验方法难以满足临时医院快速分析的需求。

        常用的半经验模型有高斯模型和ASHRAE模型[6-7]。这类模型假设稳态情形下漂浮污染物的含量在水平和垂直方向上呈高斯分布[8-9]，在实际分析中通常需要依据观测得到一些经验或半经验参数，以此来实现分析问题的简化[10]。虽然半经验模型能够满足快速分析需求，但目前的半经验模型仅仅适用于平坦的空旷地带，而难以考虑建筑甚至 建筑群对污染物扩散机制的影响[11]。

        近年来，计算流体力学 (Computational Fluid Dynamics, CFD) 技术逐步在建筑物周围的污染物扩散研究方面得到应用。依据研究对象的尺度，这类研究大致可以分为独栋建筑、单个街道以及复杂建筑群三个层次。在独栋建筑周围的污染物扩散研究方面，由于计算尺度较小，除了雷诺平均N-S方程 (RANS) 方法[12-15]，大涡模拟 (LES) 方法 [16-17]和直接数值模拟 (NDS) 方法[18]也得到了一定的应用。而在单个街道层次，研究主要采用RANS方法[19-21]和LES方法[22-24]。在复杂建筑群层次，相关研究大多采用RANS方法[25-27]，LES得到少量应用[28]。

        虽然目前已经有很多比较成熟的污染物扩散模拟方法，但将上述半经验模型或者数值模拟方法应用到本研究时，有两点困难需要克服：(1) 由于临时医院建筑群平面布置多样，简化的半经验模型难以适用；(2) 一些商用流体动力学计算平台虽然可以提供较好的分析功能和前后处理功能，但是考虑到中国多地临时医院的紧急建设需求，购买商用流体计算平台无论是时间上还是费用上都无法满足应急阶段的需要。且Fluent等一些常用流体力学软件计算资源消耗较大且计算时间较长。因此，有必要在本研究中引入计算效率高、代码开源、操作简便的CFD计算软件，以满足临时医院的紧急分析需求。

        FDS (Fire Dynamic Simulator) 是美国国家标准技术研究院开发的一款开源流体动力学计算软件。该软件最初是为了火灾模拟而开发，随着其功能的逐步丰富，目前也在空气污染模拟领域得到了应用[29-30]。FDS软件具有以下优点：(1) FDS为开源软件，便于广泛推广；(2) 采用LES模拟污染物扩散过程，能够捕捉流场和污染物的瞬态特征，更好模拟不同污染物的传染机理；(3) 与一般的商用软件相比，FDS软件能快速实现流域的网格划分，操作简便，便于工程师快速对设计方案进行评估。

        基于上述优点，本研究以FDS为软件基础，针对临时医院的快速分析需求，研发有害气体对环境影响的快速模拟方法，为临时医院设计服务。

(1) 临时医院建筑的快速建模模块；

(2) 基于云计算平台的CFD模拟和参数分析模块；

(3) 有害气体流动监测及结果可视化模块。

        其中，模块(1) 需要完成流场域、临时医院建筑、排风口、整体风环境以及有害气体监测点等的建模。模块(2) 需要完成不同风向、不同风速、不同设计方案下的影响分析。模块(3) 需要对计算结果进行处理，并生成直观的可视化效果，满足工程建设决策需要。为满足临时医院建设的急迫时间需求，本文提出了以下高效建模和分析方法。

        FDS中流场模拟和污染扩散分析需要输入以下信息：(1) 流场计算域网格；(2) 建筑物模型；(3) 风环境模型；(4) 排气口模型；(5) 有害气体监测点。具体介绍如下：

        合理的网格是CFD计算精度和效率的重要保障，与常规CFD软件的网格划分方式不同，FDS采用矩形网格对计算域进行网格划分，并能实现对压力的快速求解。因而保证了计算的高效，可以在数十分钟内得到复杂环境的大涡模拟结果。对临时医院建筑进行分析时，计算域的尺寸和边界条件应参考COST 732 [31]和AIJ [32] 等Best Practice Guidelines的要求。FDS允许不同计算域采用不同的网格密度。因此，对于医院建筑附近的计算域，推荐采用0.5 m到1 m的网格，以提升精度。由于临时医院一般高度不超过三层，因此对排风口以上10m范围以外，推荐采用2 m的网格，保证计算效率和精度的平衡。

图1 建设中的(a)火神山(b)雷神山临时医院[33]

Fig1 (a) Huoshenshan and (b) Leishenshan hospital under construction[33]

&OBST XB=X-MIN, X-MAX, Y-MIN, Y-MAX, Z-MIN, Z-MAX/

        其中，“&OBST XB=”为FDS的关键字，表示建立一个实体。X-MIN, X-MAX, Y-MIN, Y-MAX, Z-MIN, Z-MAX分别是建筑物在X、Y、Z三个方向的边界。由于FDS采用了矩形流体网格，因此实体建模也只能建立矩形实体。对于外形复杂的建筑，FDS只能采用大量矩形实体近似其外形。但是对于临时医院而言，由于其外形本来就是规则的矩形（图1），因此反而充分符合了FDS的特点。

        值得一提的是，如果分析对象建筑已经有了Sketchup模型、BIM模型、CAD模型等三维模型，则可以通过Pyrosim软件[34]，将其转换成FDS的模型，建模效率可以进一步提高。

        排气口建模包括三个内容：

(1) 排气口的空间位置

(2) 排气口的排风量

        根据FDS提供的功能，“&SURF”命令可以设定排气口的排气量（采用“VOLUME_FLOW”参数，单位m3/s），这个数据可以根据排风口联通的空调系统的设计通风量确定。本次肺炎疫情临时医院每个病房的设计排风量VR大概是600-750 m3/h。一般N个病房共用一个排风口，则可以计算得到排风口的排风量为VR∙N。

(3) 排气口有害气体的追踪

        最后，为了定量获取指定位置的有害气体浓度，可以采用FDS的“&DEVC”指令，创建有害气体监测点。具体参数包括：(1) 用“QUANTITY”参数，指定监测内容为示踪气体的比例。(2) 用“XYZ=”命令，指定监测点在模型中的坐标。另外，也可以通过剖面、等值面等方式获取有害气体在空间的宏观分布。

        具体操作为：可以通过在华为云平台（也可采用其他云计算平台）购买64核的云计算虚拟机（每台虚拟机的费用大概13-18元/小时），并安装FDS软件。FDS具有很好的并行计算功能，可以充分利用虚拟机上的并行计算资源，支持MPI、OpenMP等多种并行机制，可运行于Windows或Linux等多种操作系统。这样有多少个工况就可以购买多少台云计算虚拟机开展分布式计算，从而以较低的成本满足突发的计算资源需求。

        FDS计算完成后，可以通过以下2种途径考察有害气体的分布情况。

(1) 有害气体监测点结果

     4.5节指定的有害气体监测点可以输出“devc.cvs”文件，提供了每个监测点各时间点的精确有害气体浓度，对“devc.cvs”文件处理，则可以用曲线方式给出各监测点结果。

(2) 示踪粒子、等值面动画等

Fig2  Trace and concentration isosurface of harmful air
(Southwest wind, 1.9 m/s)

        雷神山医院为武汉市为应对肺炎疫情而建设的一座典型临时医院，位于中国湖北省武汉市江夏区军运村，设计床位数为1600张。每个病房每小时换气12次，设计新风量为500～550 m3/h，排风量为650～700 m3/h。采用上述建模方法，建立雷神山医院的三维FDS模型如图3所示。

      采用华为云计算平台，租用Intel Cascade Lake 3.0GHz 64核服务器，完成一个工况计算需要约45分钟（采用64核并行），计算占用内存约40 G。图4所示为模拟得到的不同排风口高度下1.5 m高程处的有害气体相对浓度。可见将排风口高度提高后，有害气体的相对浓度有了显著降低。

        整个分析过程从拿到设计图纸到完成计算报告，一共用时约6个多小时。其中大约3个多小时用于建模，2个多小时用于计算和结果分析。其他时间用于报告的整理和写作。满足了应急医院建设的工程周期需要。

图5 有害气体在监测点的相对浓度对比

Fig5  Comparison of relative concentrations of harmful air at monitoring points

        (1) 本文提出的分析方法，可以定量提供排风点有害气体在整个医院环境中的空间分布情况，为医院排风设计和总体规划提供参考。

        (2) FDS基于矩形网格进行建模和流体分析，非常适合临时医院建筑分析的需求，且可以极大提升计算效率。通过云计算平台以及FDS的开源免费特点，可以满足临时医院设计的迫切时间需求。

        由于本次新冠肺炎疫情非常紧迫，因此在分析中采取了很多简化假定。随着对新冠病毒传播机制研究的深入，可以引入更加合理的流体计算参数，使得模拟结果更符合实际情况。

       感谢北京科技大学吴元、张芙蓉，北京建筑大学陈越等同学提供的帮助。另外，在中信建筑设计研究总院有限公司、北京市建筑设计研究院有限公司的支持下，本方法在2020年1月24~2月2日期间，还在武汉火神山医院及北京地坛医院等工程中得到应用。