曹鼎，劉迎云，2，，穆述鑫

（1 南華大學(xué) 土木工程學(xué)院，湖南衡陽 421000；2 南華大學(xué)期刊社，湖南衡陽 421000；3 南華大學(xué) 資源環(huán)境與安全工程學(xué)院，湖南衡陽 421000）

0 引言

隨著我國經(jīng)濟(jì)的飛速發(fā)展，汽車成了人們?nèi)粘Ｉ钪胁豢苫蛉钡慕煌üぞ撸奖懔巳藗兊娜粘Ｉ?，但同時(shí)也帶來了交通安全問題與交通源污染的持續(xù)攀升。交通源顆粒物在不利于擴(kuò)散的風(fēng)環(huán)境下聚集，導(dǎo)致顆粒物濃度升高，能見度下降，對城市環(huán)境造成嚴(yán)重?fù)p害[1]。交通源污染物包括氣態(tài)污染物和大量細(xì)微的顆粒狀物體，顆粒物可吸附周圍的金屬微粒及微生物病原體等。這些細(xì)微顆粒的空氣動(dòng)力學(xué)當(dāng)量直徑大多＜10μm（PM10），特別是＜2.5μm（PM2.5）時(shí)，便可隨著呼吸進(jìn)入人們體內(nèi)，并沉積在肺部，對人體局部組織產(chǎn)生堵塞作用，從而引起哮喘、支氣管炎等疾病，提高了人群中呼吸系統(tǒng)的發(fā)病率和死亡率。

近年來，針對交通源顆粒物在住宅小區(qū)與城市街谷中的擴(kuò)散研究已經(jīng)較為深入[2-3]，但是對城市中的醫(yī)院小區(qū)研究較少。醫(yī)院小區(qū)不同于住宅小區(qū)，其相當(dāng)于大型的建筑綜合體，其中包括辦公樓、門診、住院樓、食堂樓、休閑區(qū)、家屬樓等。所以醫(yī)院小區(qū)在設(shè)計(jì)時(shí)大多注重建筑的功能布置、美觀設(shè)計(jì)及空間利用，小區(qū)內(nèi)建筑布局復(fù)雜，建筑高度不一致，很少考慮內(nèi)部風(fēng)環(huán)境與顆粒物分布對人的影響。2020 年新冠疫情肆虐全球，大量研究表明，新冠肺炎病毒顆粒的傳播方式可能與氣溶膠傳播相似，小粒徑氣溶膠可能長時(shí)間保持懸浮狀態(tài)，并受風(fēng)速影響，導(dǎo)致病毒在較大范圍傳播[4]。由于醫(yī)院小區(qū)人流量大，懸浮顆粒物易導(dǎo)致交叉感染，而小區(qū)內(nèi)的主要顆粒物來源于周邊道路，所以研究醫(yī)院小區(qū)內(nèi)風(fēng)環(huán)境與道路顆粒物在其中的分布具有重要的現(xiàn)實(shí)意義。

本文利用計(jì)算流體力學(xué)（CFD）的方法，對常見風(fēng)向與風(fēng)速條件下，醫(yī)院小區(qū)內(nèi)風(fēng)場與道路顆粒物在醫(yī)院小區(qū)濃度場中的表現(xiàn)進(jìn)行計(jì)算，并對醫(yī)院小區(qū)內(nèi)風(fēng)速及PM2.5、PM10 濃度進(jìn)行分析，結(jié)合分析結(jié)果，提出改善醫(yī)院小區(qū)舒適度和減少室外顆粒物濃度的一些建議。研究對醫(yī)院小區(qū)設(shè)計(jì)規(guī)劃具有一定的參考意義。

1 項(xiàng)目概況

某項(xiàng)目位于湖南省中部衡陽市，建設(shè)項(xiàng)目為公立醫(yī)院，地理位置為26.8°N，112.5°E，總用地面積90632.3m2，總建筑面積120235.62m2。相鄰解放大道為主干道，集賢路與青松路為次干道。衡陽地區(qū)常見風(fēng)向?yàn)楸憋L(fēng)、南風(fēng)和東北風(fēng)，常見風(fēng)速為2 級風(fēng)，平均風(fēng)速為2.5m/s。醫(yī)院小區(qū)平面布局圖如圖1 所示。

圖1 平面布局圖

2 數(shù)值模擬計(jì)算

2.1 數(shù)學(xué)模型

RNG k—ε 模型在處理小區(qū)中污染物擴(kuò)散問題時(shí)具有很高的精度［5］，因此本文采用RNG k—ε 湍流模型模擬醫(yī)院小區(qū)內(nèi)連續(xù)風(fēng)速場，空氣按不可壓縮流體處理。交通源顆粒物為固體，空氣為氣體，故為氣-固兩項(xiàng)問題。顆粒物在空氣中可視為離散相，本文采用DPM 法追蹤和描述顆粒物在醫(yī)院小區(qū)內(nèi)的擴(kuò)散運(yùn)動(dòng)，空氣流場采用歐拉法（Euler）描述，并用拉格朗日法（Lagrangian）追蹤離散顆粒的運(yùn)動(dòng)軌跡［6］。假設(shè)顆粒相為球體，并不考慮顆粒間的相互碰撞，離散相運(yùn)動(dòng)控制方程為：

式中：mp為顆粒質(zhì)量；up為顆粒速度；wp為顆粒角速度；Ip為質(zhì)量慣性矩，對于球體，Ip=0.1；T 為作用于顆粒上的總旋轉(zhuǎn)矩；FD為顆粒所受曳力；F 為顆粒所受其他作用力。

2.2 幾何建模與網(wǎng)格劃分

本文參考 《湖南省綠色建筑設(shè)計(jì)標(biāo)準(zhǔn)》(DBJ 43/T 006—2017)中的規(guī)定確定計(jì)算區(qū)域的大小。如圖2 所示，H 為醫(yī)院小區(qū)內(nèi)最高建筑高度80m，半徑5H 范圍內(nèi)為水平計(jì)算域，頂部計(jì)算區(qū)域?yàn)?H，建筑小區(qū)距離氣流出口處＞6H。整個(gè)計(jì)算區(qū)域?qū)嶋H尺寸為1490×1258×240m，其中模型面積與迎風(fēng)面計(jì)算區(qū)域截面積的比值小于4%，符合綠色建筑設(shè)計(jì)標(biāo)準(zhǔn)規(guī)范中的規(guī)定。參考當(dāng)?shù)仫L(fēng)向條件，醫(yī)院小區(qū)北部、南部和東部道路為主要污染源。計(jì)算域中，解放大道長305m，寬30m，集賢路長420m，寬16m，青松路長180m，寬16m。

圖2 模擬計(jì)算區(qū)域尺寸圖示

本研究采取非結(jié)構(gòu)化網(wǎng)格技術(shù)實(shí)行網(wǎng)格劃分，為提高數(shù)值計(jì)算精度及計(jì)算速度，對建筑與地面網(wǎng)格進(jìn)行局部加密。數(shù)值模型的網(wǎng)格單元總數(shù)約為807.9 萬，醫(yī)院小區(qū)網(wǎng)格效果圖如圖3 所示。

圖3 醫(yī)院建筑網(wǎng)格效果圖

2.3 邊界條件

2.3.1 入口邊界條件

入口邊界條件主要包括風(fēng)速和風(fēng)向數(shù)據(jù)，其中入口風(fēng)速隨高度發(fā)生變化，采用下列梯度風(fēng)：

式中：Z 表示離地面高度，m；ZR表示參考高度，m；V 表示Z高度處的風(fēng)速，m/s；VR表示參考高度（10m）的平均風(fēng)速，風(fēng)場計(jì)算取用2.5m/s；α 表示不同區(qū)域的地面粗糙度指數(shù)，參考《建筑結(jié)構(gòu)荷載規(guī)范》（GB 50009—2012）規(guī)定，本研究涉及的區(qū)域類型為建筑較密集的大城市區(qū)類，α 取值為0.22。

2.3.2 壁面邊界條件

計(jì)算域的兩側(cè)設(shè)為對稱邊界，地面邊界設(shè)定為無滑移壁面，顆粒相壁面設(shè)置為reflect，出入口為escape，空氣流動(dòng)要受到地面摩擦力的影響。

2.4 顆粒物源強(qiáng)計(jì)算

依據(jù)相關(guān)研究，道路揚(yáng)塵的排放源強(qiáng)遠(yuǎn)大于汽車尾氣的排放源強(qiáng)［7］，故在本次模擬研究中，只考慮道路揚(yáng)塵中PM2.5 與PM10對醫(yī)院小區(qū)的影響。依據(jù)《大氣可吸入顆粒物一次源排放清單編制技術(shù)指南（試行）》（000014672/2014—01379）中的附件《揚(yáng)塵源顆粒物排放清單編制技術(shù)指南》，本次研究道路為鋪裝道路，顆粒物源強(qiáng)的計(jì)算公式為：

式中，WRi為道路揚(yáng)塵源中顆粒物PMi的總排放量，μg/h；ERi為道路揚(yáng)塵源中PMi平均排放因子，g/(km·輛)；LR為道路長度，km；NR為一定時(shí)期內(nèi)車輛在該段道路上的平均車流量，輛/h；nr為不起塵天數(shù)，可使用一年中降水量大于0.25mm/d 的天數(shù)表示。

依據(jù)對衡陽市主城區(qū)道路揚(yáng)塵排放特征的相關(guān)研究［8］，衡陽市主城區(qū)各等級道路揚(yáng)塵排放因子如表1 所示。依據(jù)國家氣象科學(xué)數(shù)據(jù)中心數(shù)據(jù)，2020 年衡陽市區(qū)降雨量＞0.25mm 的天數(shù)為136天。解放大道為主干道，長度為300m，車流量為1430 輛/h，PM2.5源強(qiáng)為2.017×104μg/s，PM10 源強(qiáng)為8.5×104μg/s。集賢路為支路，長度為420m，車流量為750 輛/h，PM2.5 源強(qiáng)為1.81×104μg/s，PM10 源強(qiáng)為7.52×104μg/s。青松路為支路，長度為180m，車流量為750 輛/h，PM2.5 源強(qiáng)為7.76×103μg/s，PM10 源強(qiáng)為3.2×104μg/s。

表1 衡陽市主城區(qū)各等級道路揚(yáng)塵排放系數(shù) ERi/[g·(km·輛)-1]

3 模擬結(jié)果分析

3.1 醫(yī)院小區(qū)室外風(fēng)環(huán)境分析

風(fēng)環(huán)境好壞是評價(jià)空氣質(zhì)量的一個(gè)重要指標(biāo)，影響顆粒物在醫(yī)院小區(qū)中的擴(kuò)散。醫(yī)院小區(qū)應(yīng)具有舒適、衛(wèi)生的生活環(huán)境，所以對于風(fēng)環(huán)境有更高的要求。如圖4 所示，在北風(fēng)影響下，醫(yī)院小區(qū)北部建筑處于相對迎風(fēng)面，由于弧形綜合樓對北風(fēng)的阻擋，綜合樓南部存在風(fēng)影效應(yīng)，并出現(xiàn)大面積無風(fēng)區(qū)與負(fù)壓區(qū)，風(fēng)速小于0.5m/s。兩側(cè)為位移區(qū)，西側(cè)風(fēng)速大于3.5m/s。在住院樓2 與住院樓3 中部由于通風(fēng)區(qū)面積突然減少，此處存在狹管效應(yīng)，因此風(fēng)速明顯增大，局部風(fēng)速大于3.5m/s，不宜病人長期停留，并由于食堂樓的阻擋，氣流在食堂樓北部分流，對休閑區(qū)風(fēng)場產(chǎn)生影響。休閑區(qū)、住院樓2 北部和家屬樓中部出現(xiàn)渦流區(qū)，風(fēng)速范圍為2.3～3.5m/s，易造成渦流區(qū)垃圾堆積，影響環(huán)境。如圖5 所示，在南風(fēng)影響下，醫(yī)院小區(qū)無風(fēng)區(qū)面積較少，小區(qū)內(nèi)通風(fēng)良好。在綜合樓南部與住院樓2 北部存在渦流區(qū)，局部風(fēng)速大于3.5m/s。家屬樓屬于相對迎風(fēng)面，室內(nèi)外壓差較大，不利于病房冬季防冷風(fēng)滲透。醫(yī)院小區(qū)兩側(cè)位移區(qū)風(fēng)速較小，普遍小于1.3m/s。如圖6 所示，東北風(fēng)時(shí)醫(yī)院小區(qū)內(nèi)無風(fēng)區(qū)面積最大，住院樓2 北部存在渦流，風(fēng)速大于3.5m/s，且主要分布在住院樓2 入口處，由于狹管效應(yīng)，住院樓2 與住院樓3 中部風(fēng)速也大于3.5m/s。

圖4 北風(fēng)時(shí)1.5m高度風(fēng)速云圖與矢量圖

圖5 南風(fēng)時(shí)1.5m高度風(fēng)速云圖與矢量圖

圖6 東北風(fēng)時(shí)1.5m高度風(fēng)速云圖與矢量圖

3.2 醫(yī)院小區(qū)室外顆粒物擴(kuò)散分析

風(fēng)向?qū)︶t(yī)院小區(qū)內(nèi)顆粒物濃度分布有顯著影響，顆粒物分布受風(fēng)場影響呈流線型分布。由于PM2.5 與PM10 具有不同的物理特性與排放量，所以在醫(yī)院小區(qū)中的分布也不同。如圖7 所示，北風(fēng)時(shí)由于綜合樓的阻擋作用，PM2.5 大多向醫(yī)院小區(qū)兩側(cè)擴(kuò)散，醫(yī)院小區(qū)內(nèi)相對較少，只在休閑區(qū)有少量聚集，且PM2.5濃度低于5μg/m3；PM10 在醫(yī)院小區(qū)的分布明顯大于PM2.5，主要集中在休閑區(qū)與家屬樓南部，部分區(qū)域PM10 濃度范圍為5～20μg/m3。如圖8 所示，南風(fēng)時(shí)PM2.5 主要分布在家屬樓中部、后勤樓北部與住院樓2 西部，家屬樓與后勤樓北部PM2.5 濃度范圍為6～50μg/m3，住院樓2 西部PM2.5 濃度范圍為6～80μg/m3；PM10 主要分布在家屬樓北部與東部、休閑區(qū)南部，家屬樓北部與東部PM10 濃度范圍為10～20μg/m3，休閑區(qū)南部PM10 濃度范圍為5～15μg/m3。如圖9 所示，東北風(fēng)時(shí)，PM2.5 與PM10 在醫(yī)院小區(qū)中分布更廣，PM2.5 主要分布在休閑區(qū)、家屬樓與后勤樓，濃度范圍為3～15μg/m3；PM10 主要分布在休閑區(qū)與家屬樓，濃度范圍為6～35μg/m3。

圖7 北風(fēng)時(shí)1.5m高度PM2.5與PM10濃度云圖

圖8 南風(fēng)時(shí)1.5m高度PM2.5與PM10濃度云圖

圖9 東北風(fēng)時(shí)1.5m高度PM2.5與PM10濃度云圖

4 結(jié)語

本文通過對醫(yī)院小區(qū)室外風(fēng)環(huán)境與顆粒污染物擴(kuò)散的數(shù)值模擬研究，得出以下結(jié)論與建議：

（1）北風(fēng)與東北風(fēng)時(shí)醫(yī)院小區(qū)無風(fēng)區(qū)面積較大，北風(fēng)時(shí)無風(fēng)區(qū)主要位于醫(yī)院小區(qū)南部，東北風(fēng)時(shí)無風(fēng)區(qū)主要位于醫(yī)院小區(qū)西部休閑區(qū)，且住院樓2 與后勤樓易形成無風(fēng)區(qū)。建議把周圍高大植被換成草坪，增加氣體的流通量。住院樓1 北部易形成渦流區(qū)，風(fēng)速大于3.5m/s，且不利于污染物排放，應(yīng)避免在此處放置垃圾桶或堆放醫(yī)療廢品。食堂樓北部、住院樓1 與住院樓2 中部風(fēng)速較大，風(fēng)速范圍為1～4m/s，說明醫(yī)院小區(qū)內(nèi)建筑的分布與形狀會(huì)影響小區(qū)內(nèi)的風(fēng)速，在小區(qū)內(nèi)風(fēng)速較高區(qū)域可以通過種植樹木來改善小區(qū)風(fēng)環(huán)境；

（2）由于道路揚(yáng)塵中PM10 濃度大于PM2.5，所以醫(yī)院小區(qū)中PM10 分布更廣，濃度相對更高。家屬樓與后勤樓周圍風(fēng)速較小，PM2.5 易聚集且濃度較高，最高濃度達(dá)60μg/m3，病人長期停留可能會(huì)加重病情。PM10 主要集中在休閑區(qū)、家屬樓與后勤樓北部，最高濃度達(dá)35μg/m3，應(yīng)注意顆粒物防治。道路揚(yáng)塵在不同風(fēng)向的影響下會(huì)對醫(yī)院小區(qū)不同區(qū)域產(chǎn)生不同的影響，由于綠植和水對顆粒物有阻擋和吸收作用，所以在顆粒物濃度較高區(qū)域可以通過種植綠植、設(shè)置噴泉或圍欄等方法來降低顆粒物濃度；

（3）在醫(yī)院小區(qū)建設(shè)之初，通過CFD 對風(fēng)環(huán)境與顆粒物濃度進(jìn)行模擬分析，可在醫(yī)院小區(qū)的優(yōu)化設(shè)計(jì)中起到積極作用，為建設(shè)舒適健康的醫(yī)院小區(qū)提供科學(xué)保障。