近年来各类社区安全事件频频发生^[1]。以北京市为例，受大风天气影响，树倒墙塌、高空坠物、广告牌散落等情况时有发生，严重影响居民安全。大风天气的预报预警一般针对全城，社区范围内难以识别危险区域，风险防范难度大^[2]。同时，社区内火灾、燃气泄漏等突发事件数量近年来维持在较高水平^[3]，如2010年上海社区火灾事故、2011年北京社区燃气爆炸事件、2018年杭州社区纵火事件，这些突发事件对社区居民安全造成了重大影响。

针对社区安全事件频发的现状，传统的应急响应模式由于缺乏灾害现场的背景流场情况，影响事故发展推演及影响范围评估^[4]。在火灾、燃气泄漏、大风等灾害情境下，决策者和救援队伍由于无法快速获得准确的背景流场信息而无法准确定位社区危险区域^[5]。因此，快速获得满足精度要求的社区实际流场，并根据背景流场指导社区风险防范、应急决策与救援，对于保障社区安全具有重要意义。

在社区风险防范、应急决策与救援领域，传统的数值模拟技术计算时间长，在火灾、燃气泄漏、大风等灾害情境下无法做到快速应对^[6]。近年来，随着人工智能技术的发展，机器学习与数值计算相结合的方法得到了快速发展^[7-10]。对于涉及海量数据的数值模拟研究而言，神经网络模型能从大量模拟结果中学习到无法人工提取的统计规律，可以弥补数值模拟计算效率低、消耗资源大等弊端。深度神经网络是指具有2个或2个以上隐藏层的神经网络模型，相较于传统神经网络，其更多的隐藏层为模型提供了更强的抽象能力。

众多学者将各类神经网络模型与流场结合，开展了大量研究。Bhatnagar等^[11]将卷积神经网络运用到低速翼型的速度场和压力场估计中，使用计算流体力学结果作为训练集，为高精度的基于深度学习的流场预测提供了参考。Sekar等^[12]基于大量仿真数据，利用卷积神经网络建立了深度学习模型对翼型绕流进行估计。Wu等^[13]使用生成对抗网络来估计流场，该模型将翼型参数映射到翼型压力场，但是存在训练稳定性问题。Memarzadeh等^[14]基于长短时记忆神经网络模型进行风速预报，将实测风向和气象参数作为输入变量预测风速值。大量神经网络模型的研发^[15-16]为流场快速重构提供了参考。

本研究以社区风场为例，通过数值模拟与神经网络的结合形成社区风场特征点提取方法。利用提取的特征点风速及风场数据训练神经网络模型，实现特征点风速到风场数据的预测，从而快速获得满足精度要求的社区背景风场，辅助危险区域的识别。具体研究步骤如下：

1) 利用数值模拟技术构建社区风场数据库。建立社区三维模型，使用OpenFOAM完成单次模拟，获得社区风场，通过shell脚本设置不同的入流风速条件，批量运行模拟算例，提取并存储稳定风场数据构建社区风场数据库。

2) 利用聚类算法选取一定数量的社区特征点。选取的特征点可以表征该社区风场信息，并根据神经网络模型的训练与预测结果不断优化特征点选取结果。

3) 利用神经网络模型完成特征点数据到风场数据的训练与预测。根据训练与预测的结果优化该社区特征点选取结果，直至神经网络模型的预测精度及特征点数目达到预定要求。

1 构建流场数据集 1.1 场景建模

以北京市石景山区某社区为研究对象。根据该社区三维实景图和建筑平面图，建立社区场景三维模型。该模型三维信息分别为长600 m、宽300 m、高120 m，三维实景图如图 1所示，场景三维建模图如图 2所示。

1.2 数值模拟

采用OpenFOAM进行数值模拟，具体包括生成计算网格、定义物理模型、设置边界条件、设置求解控制、数值计算与后处理5个步骤^[17]。

生成计算网格包括构造面特征、生成背景网格、生成整体网格3个部分：构造面特征用于提取三维建模的几何特征，为生成模型整体网格作准备；背景网格用于定义计算域尺寸及网格大小；最后通过网格切割与网格贴合生成计算区域整体网格。每个网格的大小为5 m×5 m×5 m，边界条件如表 1所示。计算区域整体网格划分如图 3所示。

计算过程中，建筑物附近的计算区域需要更加精确的计算来保证结果的准确性，通过面加密的方式对几何表面相交的网格进行加密。网格细节展示如图 4所示。

本算例计算的是不可压缩湍流流动，未考虑传热，采用k-ε湍流模型^[18-19]。

本算例采用SimpleFoam求解器，迭代计算500步。SimpleFoam为OpenFOAM中一个基于有限体积法的用于求解稳态不可压缩Newton流体Navier-Stokes方程的求解器。在SimpleFoam中，压力和速度的耦合通过SIMPLE算法进行计算，且没有考虑其他体积力，如重力等。

后处理步骤通过paraview软件完成，可实现流场计算结果的可视化。社区风场计算结果示例如图 5所示。

按照上述数值模拟步骤设置模板算例，编写shell脚本重复调用模板算例完成批量仿真任务。shell脚本依次实现以下功能：模板算例定位、网格生成、参数修改赋值、算例执行、风场提取、删除中间过程文件。全部算例运行完毕后风场数据集构建完成，最后从中提取的测试集与训练集风场共9 943个。

2 流场特征点选取 2.1 标准差计算

以社区风场平面为研究对象，针对场内每个点的速度数据，分别计算其东西方向(x方向)和南北方向(y方向)的标准差，加和作为该点的最终标准差。计算得到的6 681个网格点的标准差范围在0.34~14.82 m/s。标准差统计信息如表 2所示，标准差云图如图 6所示。

不同的建筑布局及不同的入流风速，会产生不同的标准差分布，因此可以选择不同大小的标准差所对应的点作为特征点，用这些特征点表示区域建筑布局信息及入流风速信息，进而通过点的速度预测区域场的速度。

2.2 特征点选取

采用k-means聚类算法对6 681个点的标准差聚类，聚类前将点的坐标与其标准差一一对应，并将标准差划分为3级：1级指标对应小标准差0~9.1 m/s，2级指标对应中标准差9.1~10.2 m/s，3级指标对应大标准差10.2~14.8 m/s。

对1级标准差内的所有点进行坐标聚类，得到3个聚类中心；对2级标准差内的所有点进行坐标聚类，得到3个聚类中心；对3级标准差内的所有点进行坐标聚类，得到4个聚类中心。以聚类中心处的速度为输入，以风场整体速度为输出进行神经网络的训练及预测，并根据训练与预测的结果调整标准差级别划分标准与相应的聚类中心位置，直至网络训练与预测结果的精度满足要求。聚类中心点即为特征点。

最终选取的10个特征点的标准差大小、坐标值(x，y)及索引号信息如表 3所示。

特征点选取位置如图 7所示。

3 神经网络模型 3.1 数据预处理

对输入和输出数据作归一化处理，将所有风场内的风速大小映射至[-1, 1]内，映射过程按式(1)进行，

$ \bar{v}=\frac{v-v_{\min }}{v_{\max }-v_{\min }} . $

(1)

其中: v为归一化后的数据，v为原数据，v_min为样本所有速度数据中的最小值，v_max为样本所有速度数据中的最大值。

风场数据库共有9 943个风场，以7 917个风场为训练集、2 026个风场为测试集进行神经网络训练。

3.2 网络结构设计

数值模拟得到不同入流风速下的风场数据。以特征点处的风速为输入，以平面风场为输出，搭建神经网络模型。

网络层数分别设置为3、4、5、6层；激活函数选择Tanh函数；学习率设置为用分数减缓的方式调整，即随着迭代次数的增加逐步减小学习率；选择小批量梯度下降法，每批风场数量设置为16；损失函数选择基于距离度量的损失函数，如式(2)所示，度量样本真实值和模型预测值之间的距离。

$ L=\frac{1}{N} \sum\limits_{i=1}^N\left|\frac{y_i-\hat{y}_i}{y_i}\right| $

(2)

其中: y_i为真实值，$\hat{y}_i$为模型预测值，N为样本总量。

4 结果

经过训练后的模型在测试集中速度预测误差统计如表 4所示。其中：x方向为东西向，y方向为南北向。由表 4可知，对于1~2 m/s分速度区间，模型预测的平均相对误差为11.9%；对于2~5 m/s分速度区间，模型预测的平均相对误差为6.0%；对于5~10 m/s分速度区间，模型预测的平均相对误差为3.2%；对于10 m/s以上分速度区间，模型预测的平均相对误差为3.5%。

从神经网络模型训练与预测结果来看(表 4)，x方向和y方向1 m/s以上速度的预测值平均相对误差分别为5.8%和6.2%；x方向和y方向5 m/s以上速度的预测值平均相对误差分别为3.4%和3.0%。

随机挑选某模拟风场的数值模拟速度与模型预测速度如图 8和9所示。图 8所示为数值模拟产生的风场；图 9所示为训练后的神经网络模型根据特征点风速产生的模型预测的风场。可见，对于该风场，以数值模拟风场为基准，神经网络模型预测的准确率为95.3%。

基于特征点采用神经网络模型对社区风场进行预测的结果，与采用数值计算方法获得的社区风场结果基本吻合，可见神经网络模型的精度能够满足实际需要。

此外，数值模拟需约十几min产生一个风场，而训练好的神经网络模型只需ms量级时间即可产生对应的高精度流场，用时大大缩短，对于紧急情况下快速获取背景风场具有重要应用价值。

5 结论

本文首先完成社区场景建模，包含社区内各类建筑物及构筑物，再利用数值模拟技术产生社区模拟风场，并编写shell脚本完成批量仿真，得到社区风场数据库。采用k-means聚类算法得到10个特征点，以特征点处的速度为输入，以风场为输出进行神经网络模型的训练及预测，得到神经网络模型。

该神经网络模型可以快速获得满足精度要求的社区背景风场，从而辅助社区风险防范、应急决策与救援。与数值模拟约十几min产生一个场的时间相比，训练好的神经网络模型只需ms量级即可产生对应的流场，时间大大缩短。该模型以社区内实际点位为特征点训练模型，下一步将在特征点处安装风速风向传感器采集真实的风速风向数据，利用真实的风速风向数据修正神经网络模型并验证预测结果。