街道峽谷對(duì)稱性對(duì)污染物擴(kuò)散的影響

楊 方，鐘 珂，亢燕銘（東華大學(xué)環(huán)境科學(xué)與工程學(xué)院，上海 201620）

街道峽谷對(duì)稱性對(duì)污染物擴(kuò)散的影響

楊 方，鐘 珂，亢燕銘*（東華大學(xué)環(huán)境科學(xué)與工程學(xué)院，上海 201620）

采用數(shù)值模擬的方法，研究了行列式和錯(cuò)列式街谷兩側(cè)建筑物的對(duì)稱性對(duì)街谷內(nèi)和下游建筑尾流區(qū)污染物濃度的影響，并引入3個(gè)用于評(píng)價(jià)室內(nèi)空氣環(huán)境的指標(biāo)來描述街谷中人員區(qū)域內(nèi)的空氣質(zhì)量.結(jié)果表明，平均殘留時(shí)間（ART）和有效通風(fēng)率（PFR）可用于評(píng)價(jià)不對(duì)稱街谷內(nèi)的空氣質(zhì)量.街谷的不對(duì)稱性對(duì)街谷內(nèi)部和下游建筑尾流區(qū)污染物濃度分布均有明顯影響.下游建筑越低，則尾流區(qū)的濃度越高.行列式和錯(cuò)列式街谷中人員活動(dòng)空間的無量綱平均濃度分別在建筑高度比例為7/3和7/2時(shí)最大.為在城市規(guī)劃中盡可能避免出現(xiàn)最不利比例的街谷，根據(jù)平均殘留時(shí)間和有效通風(fēng)率等指標(biāo)，給出了上下游建筑高度比與最不利街谷改善率的關(guān)系.

非對(duì)稱街道峽谷；街谷上下游建筑高度比；最不利比例街谷；有效通風(fēng)率；平均殘留時(shí)間

隨著我國機(jī)動(dòng)車保有量的逐年增加，道路交通排放已成為導(dǎo)致城市空氣環(huán)境質(zhì)量持續(xù)惡化的主因之一.城市街道峽谷（簡稱“街谷”）中的流場(chǎng)性質(zhì)和機(jī)動(dòng)車尾氣擴(kuò)散特征與街谷結(jié)構(gòu)密切相關(guān).街道兩側(cè)的臨街建筑高度相等時(shí)，構(gòu)成對(duì)稱型街谷；道路兩側(cè)建筑高度相差較大時(shí)，則構(gòu)成非對(duì)稱街谷.隨著城市規(guī)劃的多樣性發(fā)展，不對(duì)稱街谷越來越普遍.

大量文獻(xiàn)針對(duì)街谷流場(chǎng)和污染物分布特征進(jìn)行了研究［1-9］，但大多數(shù)都是關(guān)于街道兩側(cè)建筑物高度相等的情況［1-6］，而關(guān)于非對(duì)稱街谷內(nèi)部流場(chǎng)和污染物濃度分布特征的研究較少.Xie等［8］研究了街道兩側(cè)建筑物高度的比值對(duì)街谷內(nèi)渦流結(jié)構(gòu)的影響.何澤能等［9］通過模擬不同高度比的街谷內(nèi)流場(chǎng)，指出峽谷的對(duì)稱性對(duì)其內(nèi)部氣流場(chǎng)有顯著影響.Hang等［10］的研究表明，在非對(duì)稱建筑群中，相對(duì)于湍流作用，水平和垂直方向的平均流對(duì)污染物擴(kuò)散的影響較大.汪立敏等［11］指出，當(dāng)下游建筑高度低于上游建筑時(shí)，街谷內(nèi)的流場(chǎng)及濃度場(chǎng)與對(duì)稱街谷時(shí)相似.當(dāng)下游建筑高于上游建筑時(shí)，街谷內(nèi)的流場(chǎng)較之對(duì)稱街道峽谷有較大的變化.羅昔聯(lián)等［12］的研究表明，上游建筑高于下游建筑的非對(duì)稱街谷中，街谷內(nèi)顆粒物濃度遠(yuǎn)高于同類型其它幾何尺寸街谷.

盡管上述研究考慮了非對(duì)稱街谷中的氣流運(yùn)動(dòng)和污染物空間分布特征，但關(guān)于不對(duì)稱性對(duì)街谷內(nèi)人員活動(dòng)區(qū)空氣品質(zhì)的影響效果尚未涉及，缺乏可供城市規(guī)劃設(shè)計(jì)師參考的合理建議.為深入了解街谷不對(duì)稱性形成的氣流特征對(duì)機(jī)動(dòng)車尾氣擴(kuò)散的影響，以尋求街谷不對(duì)稱程度與街谷內(nèi)空氣品質(zhì)的量化關(guān)系，本文利用計(jì)算流體力學(xué)方法，研究分析三維街谷中上、下游建筑高度比對(duì)街谷內(nèi)以及下游建筑尾流區(qū)氣流和污染物濃度分布的影響.同時(shí)，為了更好地評(píng)價(jià)氣流對(duì)街谷中污染物的清除效果，引入相關(guān)的通風(fēng)效率指標(biāo)對(duì)街谷內(nèi)空氣質(zhì)量進(jìn)行評(píng)價(jià)，為城市的合理規(guī)劃提供理論依據(jù)和參考數(shù)據(jù).

1 物理模型與計(jì)算方法

為獲得街谷不對(duì)稱性對(duì)氣流和污染物擴(kuò)散影響的詳細(xì)信息，對(duì)街道邊界做了簡化，分別對(duì)街道兩側(cè)建筑為規(guī)則排列的行列式布局和錯(cuò)列式布局的街道模型結(jié)構(gòu)進(jìn)行數(shù)值模擬.圖1給出了行列式和錯(cuò)列式建筑布局構(gòu)成的街谷，以及計(jì)算域尺寸.圖中還給出了隨后分析中涉及到的部分幾何參數(shù)（圖中A-A為垂直方向剖面，位于建筑物1/4長度處；L1、L2為街谷呼吸高度處沿街道分布的直線，分別靠近街谷端部和中心）.

圖1 不同建筑布局構(gòu)成的街谷模型和計(jì)算域的尺寸Fig.1 Sketch of street canyon models for different building arrangements and the computational domain

設(shè)所有臨街建筑長60m （L），寬15m （B），街寬為30m （W），街道同側(cè)相鄰建筑間距S=20m，因而臨街建筑物在y方向成周期性排列.考慮街谷內(nèi)污染物擴(kuò)散最不利的情況，設(shè)來流風(fēng)向垂直于街谷軸線.此時(shí)，所形成的街谷流場(chǎng)也呈周期性空間分布.因此，為節(jié)約計(jì)算成本，采用鏡像法，通過設(shè)置對(duì)稱面，僅計(jì)算圖1（a）和（b）中虛線區(qū)域.本文針對(duì)臨街為多層建筑的街谷進(jìn)行研究，建筑物最大高度為21m.數(shù)值模擬中，為盡量消除計(jì)算區(qū)域有限導(dǎo)致的端部效應(yīng)，模擬時(shí)做如下處理［13-14］，迎風(fēng)的首排建筑距求解區(qū)域上游邊界為12H，末排建筑距區(qū)域下游出口17H，建筑頂部距上方自由邊界面5H.如圖1（c）所示.使用非結(jié)構(gòu)網(wǎng)格劃分計(jì)算空間，為保證給出更細(xì)致的近壁面處的計(jì)算結(jié)果.經(jīng)過網(wǎng)格獨(dú)立性測(cè)試后，建筑物附近最細(xì)網(wǎng)格為0.3m，模擬區(qū)域的網(wǎng)格總數(shù)為310～370萬.

設(shè)街谷兩側(cè)上游和下游建筑物高度分別為H1和H2，分別針對(duì)不同的H1和H2組合進(jìn)行模擬計(jì)算，具體組合情況見表1.

設(shè)模擬區(qū)域的入口處流速均勻，并將其邊界設(shè)為速度進(jìn)口；由于出口處可視為充分發(fā)展湍流，故將出口邊界條件設(shè)為outflow.地面和建筑表面為無滑移壁面，計(jì)算域壁面則設(shè)為有滑移壁面.

湍流特性和風(fēng)速的垂直特性受地形的影響，并與地形的粗糙度有著密切的關(guān)系.計(jì)算域入口風(fēng)速的垂直特性由冪指數(shù)律給出，如下式所示：

表1 上、下游建筑的高度Table 1 Heights of upwind and downwind buildings

式中：z為任一高度，m；Vz為任一高度處的平均風(fēng)速，m/s；α為地面粗糙度指數(shù)，地形越粗糙，地表對(duì)氣流的阻滯作用越強(qiáng)，α也越大，根據(jù)相關(guān)文獻(xiàn)［16］，這里α=0.25；h為標(biāo)準(zhǔn)參考高度，m；Vh為在此高度處測(cè)得的標(biāo)準(zhǔn)參考風(fēng)速，m/s，在本文中定義最高建筑屋頂高度處（H=21m）的風(fēng)速為屋頂風(fēng)速Vr.我國長江中下游地區(qū)平均風(fēng)速通常在1.4～5.3m/s［17］，對(duì)應(yīng)的流動(dòng)Reynolds為2×106～7.6×106.此時(shí)流場(chǎng)處于自模區(qū)范圍內(nèi)［18］，即風(fēng)速的變化不會(huì)對(duì)研究結(jié)論有較大的影響.本文對(duì)屋頂來流風(fēng)速Vr=2.0m/s時(shí)進(jìn)行模擬計(jì)算.

為模擬機(jī)動(dòng)車尾氣排放，將交通污染源視作連續(xù)均勻散發(fā)的線源，設(shè)置于街谷底部，它的中心軸線與街谷的中心線重合，寬10m，高0.5m.本文研究粒徑為1μm的顆粒污染物的擴(kuò)散和分布規(guī)律，污染源強(qiáng)度q0=0.02kg/s.為了便于比較，顆粒物濃度用無因次濃度C*的形式表示：

式中：Ci是i點(diǎn)處的顆粒物濃度，kg/m3；Ls為污染源長度，m.

2 模型的合理性驗(yàn)證

為了對(duì)所采用的數(shù)值模型進(jìn)行可靠性驗(yàn)證，分別以Li等［19］和Meroney等［20］關(guān)于街谷流場(chǎng)和濃度分布的實(shí)驗(yàn)參數(shù)（包括街谷幾何參數(shù)、流動(dòng)參數(shù)和污染物參數(shù)）為基本數(shù)據(jù)，利用本文擬采用的離散模式、網(wǎng)格劃分方法和計(jì)算控制等對(duì)二者的實(shí)驗(yàn)分別做了數(shù)值再現(xiàn).其中速度場(chǎng)采用Li等［19］的水槽實(shí)驗(yàn)結(jié)果，并沿水平（x方向）和垂直（z方向）兩個(gè)方向同時(shí)考察風(fēng)速廓線；在速度場(chǎng)得到驗(yàn)證后，繼續(xù)求解濃度場(chǎng)，對(duì)比數(shù)據(jù)采用Meroney等的實(shí)驗(yàn)結(jié)果.圖2給出了數(shù)值計(jì)算結(jié)果與上述實(shí)驗(yàn)研究結(jié)果的比較.由于近年來，有研究建議用大渦模擬分析街谷流動(dòng)問題，圖中也以Liu等［21］的模擬為例，給出了大渦模擬速度場(chǎng)計(jì)算曲線.

圖2 風(fēng)速和濃度數(shù)值計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)的比較Fig.2 Comparisons of wind speed and concentration between numerical results and the experimental data Vzx和Vzz為分z方向的風(fēng)速分量

由圖2（a）和圖2（b）可以看到，本文與Liu等［21］的大渦模擬結(jié)果在街谷內(nèi)大部分區(qū)域均與實(shí)驗(yàn)結(jié)果吻合得較好，這說明本文的湍流模型能夠正確計(jì)算街谷內(nèi)流場(chǎng)情況.圖2（c）中的背風(fēng)區(qū)濃度計(jì)算結(jié)果與實(shí)測(cè)也符合得較好，迎風(fēng)區(qū)則完全吻合.表明本文采用的計(jì)算方法可以很好地預(yù)測(cè)街谷內(nèi)的氣流流動(dòng)和污染物濃度分布.

3 結(jié)果與討論

3.1 街谷對(duì)稱性對(duì)建筑周圍流場(chǎng)和顆粒物平均濃度的影響

街谷兩側(cè)建筑物高度不同時(shí)，來流受到的阻擋作用不同，從而導(dǎo)致街谷內(nèi)渦流結(jié)構(gòu)與對(duì)稱性密切相關(guān).圖3僅以行列式的A-A剖面［圖1（a）］為例，給出了H1/H2不同時(shí)，建筑周圍氣流的流線和顆粒物濃度分布情況.

圖3 行列式街谷中A-A剖面上（y=25m）流場(chǎng)和污染物無量綱濃度分布Fig.3 Flow field and dimensionless concentration distributions at A-A plane （y = 25m） for regular arrangement

從圖3可以看出，當(dāng)H1/H2不同時(shí)，不僅街谷內(nèi)的氣流分布和濃度分布不同，下游建筑尾流區(qū)域的流場(chǎng)和濃度場(chǎng)分布也不同.對(duì)于上階梯型街谷［圖3（a）］，來流在下游建筑的阻擋下，在街谷內(nèi)存在兩個(gè)順時(shí)針旋轉(zhuǎn)的渦流，在橫向方向上幾乎各占街谷的一半.在渦流的作用下，機(jī)動(dòng)車排放的污染物被輸送到背風(fēng)面，再通過背風(fēng)面屋頂處排出街谷，并很快被來流稀釋，已流出污染物幾乎沒有再次流進(jìn)街谷，街谷中機(jī)動(dòng)車排放的污染物流入下游建筑尾流區(qū)的量很小.在完全對(duì)稱的等高街谷中，街谷內(nèi)僅有一個(gè)順時(shí)針渦流［圖3（b）］，在渦流的作用下，污染物在背風(fēng)面和地面產(chǎn)生堆積，濃度隨著高度的遞增而減少.此時(shí)從背風(fēng)面屋頂排出的污染物一部分再次流進(jìn)街谷，一部分流入尾流區(qū).對(duì)于下階梯型街谷［圖3（c）］，街谷中存在一個(gè)范圍相對(duì)較大的渦流，渦流通過下游建筑屋頂延伸至尾流區(qū)，因而街谷中污染物容易被帶出街谷，進(jìn)入尾流區(qū)，使得尾流區(qū)濃度相對(duì)較高.

街谷內(nèi)不同高度處的污染物濃度有很大區(qū)別，因此以不同的高度（或高度范圍）的污染物濃度數(shù)據(jù)進(jìn)行分析討論產(chǎn)生的結(jié)論可能是不同的，降低街谷內(nèi)行人的呼吸濃度是改善街谷空氣品質(zhì)的根本目的，因此本文針對(duì)人員活動(dòng)空間內(nèi)（高度在2m以下）的平均濃度進(jìn)行分析討論.并統(tǒng)一采用人員活動(dòng)空間內(nèi)（高度在2m以下）的平均濃度來評(píng)價(jià)街谷空氣品質(zhì).

本文關(guān)注尾流區(qū)內(nèi)距離下游建筑背風(fēng)面30m以內(nèi)區(qū)域的污染物濃度.圖4給出了街谷中和下游建筑尾流區(qū)的人員活動(dòng)空間內(nèi)污染物濃度的無量綱值與建筑高度比（即不對(duì)稱性）的關(guān)系，同時(shí)也給出了尾流區(qū)和街谷中人員活動(dòng)空間污染物濃度之比k，即：

從圖4（a）可以看出，行列式街谷人員活動(dòng)空間內(nèi)平均濃度受對(duì)稱性影響很大，在H1/H2=7/3時(shí)，街谷人員活動(dòng)空間內(nèi)平均濃度最大.而尾流區(qū)人員活動(dòng)空間中的平均濃度幾乎不隨對(duì)稱性變化而改變，在H1/H2=3/7～7/2，幾乎為常數(shù)，并且相對(duì)街谷人員活動(dòng)空間內(nèi)濃度都很小.這使得尾流區(qū)和街谷中人員活動(dòng)空間污染物濃度的比值k幾乎完全受控于街谷濃度.還可以看到，當(dāng)街谷為上階梯極不對(duì)稱（H1/H2=1/7）時(shí)，下游建筑有效地保護(hù)了尾流區(qū)；當(dāng)街谷為下階梯極不對(duì)稱（H1/H2=7/1）時(shí)，上游建筑對(duì)污染物擴(kuò)散的阻擋作用迫使機(jī)動(dòng)車尾氣在近地面區(qū)域擴(kuò)散開，由于下游建筑過低，造成尾流區(qū)的人員活動(dòng)空間濃度相對(duì)于街谷內(nèi)人員活動(dòng)空間濃度不可忽略.

圖4 下游建筑尾流區(qū)與街谷人員活動(dòng)空間的無量綱濃度及其濃度比Fig.4 Dimensionless pollutant concentrations in the downstream building wake regions and street canyons and the ratios between the two

由圖4（b）可以看到，對(duì)于H1/H2＜1的上階梯式街谷，錯(cuò)列式街谷尾流區(qū)的人員活動(dòng)空間幾乎不受污染；而在下階梯式街谷中，尾流區(qū)濃度隨著下游建筑高度的降低而增大，同樣是因?yàn)橄掠谓ㄖ降停瑢?duì)尾流區(qū)的保護(hù)作用越弱.街谷人員活動(dòng)空間內(nèi)平均濃度最大值出現(xiàn)在H1/H2=7/2.尾流區(qū)和街谷中人員活動(dòng)空間污染物濃度之比k隨著H1/H2單調(diào)增大.

3.2 街谷不對(duì)稱性對(duì)街谷內(nèi)濃度分布規(guī)律的影響離地1.5m高度處是行人平均呼吸高度.為了比較街谷不對(duì)稱性對(duì)街谷內(nèi)顆粒物濃度水平分布的影響，圖5和圖6分別給出了1.5m高度處行列式和錯(cuò)列式街谷在靠近上游建筑端部和中心［具體位置見圖1（a）和圖1（b）中的直線L1和L2］處顆粒物濃度的水平變化曲線.

圖5 行列式街谷內(nèi)行人呼吸高度處（z = 1.5m）顆粒物濃度分布Fig.5 Dimensionless pollutant concentration distributions along the streets at the breathing levels （z = 1.5m） inregular canyons

由圖5（a）可以看出，在靠近行列式街谷缺口處，H1/H2=7/3是臨界點(diǎn)，這時(shí)最大濃度值發(fā)生在街谷中央附近.當(dāng)H1/H2≤7/5時(shí)，濃度峰值發(fā)生在背風(fēng)面附近.這是由于下游建筑高，對(duì)來流產(chǎn)生的較強(qiáng)阻擋作用所誘發(fā)的渦流將機(jī)動(dòng)車排放的污染物輸送至背風(fēng)面，如圖3（a）和（b）所示.此外還可以看到，H1/H2越接近1，即對(duì)稱性越強(qiáng)，街谷內(nèi)濃度最大值越高.當(dāng)H1/H2＞7/3時(shí)，盡管濃度沿街道寬度方向變化較小，但明顯可以看到迎風(fēng)面附近濃度最大.同時(shí)街谷越對(duì)稱，即H1/H2越接近1，濃度越大.

從圖5（b）可以看出，不論對(duì)稱性如何，街谷中心區(qū)的濃度分布特征均為背風(fēng)區(qū)濃度高于迎風(fēng)區(qū)濃度.同樣，H1/H2=7/3是臨界點(diǎn)，此時(shí)街谷中的濃度最大值最高，隨著H1/H2偏離7/3，最大濃度逐漸變小.

圖6 錯(cuò)列式街谷內(nèi)行人呼吸高度處（z = 1.5m）顆粒物濃度分布Fig.6 Dimensionless pollutant concentration distributions along the streets at the breathing levels （z = 1.5m） instaggered canyons

由圖6（a）和6（b）可以看到，錯(cuò)列式街谷中，除了極不對(duì)稱的下階梯街谷中靠近上游建筑端部位置外，街谷不對(duì)稱性不會(huì)改變背風(fēng)區(qū)濃度偏高的特性，但濃度最大值會(huì)受到對(duì)稱性的影響.與圖5比較，可以認(rèn)為，街谷不對(duì)稱性對(duì)街道內(nèi)污染物濃度分布規(guī)律的影響，在行列式街谷中的作用明顯大于錯(cuò)列式街谷.

為了綜合比較建筑布局對(duì)街谷行人呼吸域內(nèi)污染物擴(kuò)散的影響，引入3種原本用于評(píng)價(jià)室內(nèi)空氣質(zhì)量的指標(biāo)來描述街谷呼吸域內(nèi)的空氣質(zhì)量.

有效通風(fēng)率（PFR）：可以理解為將污染物趕出空間所需要的空氣量，它反應(yīng)了空氣清除污染物的能力，其表達(dá)式為［22］：

式中：qp為污染物的產(chǎn)生速率，kg/s；ρ為空氣密度，kg/m3；Cp為空間內(nèi)污染物平均濃度，kg/kg.

訪問頻率（VF） 街谷內(nèi)部分污染物在渦流作用下，可能出現(xiàn)被氣流帶出街谷后有再次進(jìn)入街區(qū)的情況，訪問頻率VF用于表示一個(gè)污染氣體微團(tuán)進(jìn)入空間的次數(shù).空間中所有污染空氣微團(tuán)的平均訪問次數(shù)可用公式（5）計(jì)算，VF值越小，表示較少的污染氣體返回空間，氣流清除污染物的能力較高.

式中：Δqp是顆粒物進(jìn)入空間的流量，kg/s；Si為空間邊界平面i的面積，m2；u是空氣進(jìn)入空間通過邊界平面i時(shí)的平均速度，m/s；c是邊界平面i的污染物濃度，kg/kg；n為空氣進(jìn)入空間通過的邊界平面總數(shù)，u′是速度的波動(dòng)量，m/s；c′是濃度的波動(dòng)量，kg/kg.uc與ρSi相乘代表空間進(jìn)氣量，kg/s與ρSi相乘代表湍流擴(kuò)散部分的進(jìn)氣量，kg/s.

平均殘留時(shí)間（ART） 對(duì)于一個(gè)顆粒物或污染氣體微團(tuán)來說，殘留時(shí)間被定義為污染物從最開始進(jìn)入空間到最后離開空間的時(shí)間.可用下式計(jì)算：

應(yīng)用平均值的概念，將VF和ART相乘表明空間內(nèi)所有顆粒物的平均殘留時(shí)間.V為空間體積，m3.

根據(jù)以上公式分別計(jì)算H1/H2不同時(shí)街谷各項(xiàng)指標(biāo)的變化，圖7分別給出了街谷內(nèi)平均濃度和以上3種指標(biāo)隨著H1/H2的變化曲線.

圖7 H1/H2對(duì)街谷人員活動(dòng)空間內(nèi)各項(xiàng)指標(biāo)的影響Fig.7 Effect of ratios H1/H2on indices in the pedestrian domains

由圖7可以看到，除了訪問頻率（VF），行列式街谷中另外兩個(gè)指標(biāo)和平均濃度的最不利值都出現(xiàn)在7/3，錯(cuò)列式街谷中另外3個(gè)指標(biāo)的最不利值都出現(xiàn)在7/2，表明平均殘留時(shí)間（ART）、有效通風(fēng)率（PFR）和街谷人員空間平均濃度（C）都可以用于評(píng)價(jià)街谷不對(duì)稱性對(duì)濃度分布的影響.本文中污染物產(chǎn)生速率恒定，此時(shí)有效通風(fēng)率（PFR）與濃度（C）成反比關(guān)系，使用濃度（C）評(píng)價(jià)街谷內(nèi)空氣品質(zhì)更直接明了，但若顆粒物產(chǎn)生速率或空氣密度隨時(shí)間發(fā)生變化，則適合采用有效通風(fēng)率（PFR）評(píng)價(jià)街谷內(nèi)的空氣品質(zhì).

圖7表明，行列式街谷和錯(cuò)列式街谷中，人員空間濃度最大值分別在H1/H2=7/3和7/2時(shí)，均為下階梯式街谷.因此在規(guī)劃街道時(shí).應(yīng)結(jié)合當(dāng)?shù)刂鲗?dǎo)風(fēng)向.盡可能避免這種街谷出現(xiàn).

3.3 街谷不對(duì)稱性與最不利街谷空氣品質(zhì)改善率

在城市規(guī)劃中，應(yīng)盡可能偏離最不利街谷比例，為此定義相對(duì)于最不利街谷時(shí)，實(shí)際街谷空氣品質(zhì)的改善率ηi為：

式中：Ai為待規(guī)劃街谷的第i個(gè)評(píng)價(jià)指標(biāo)的預(yù)測(cè)值；Ai0為最不利比例街谷內(nèi)第i個(gè)評(píng)價(jià)指標(biāo)的預(yù)測(cè)值.在本文條件下，有效通風(fēng)率PFR與平均濃度C的改善率完全一致.但為了在污染物強(qiáng)度隨時(shí)間發(fā)生變化時(shí)也能夠全面反應(yīng)空氣品質(zhì)的改善情況，圖8給出街谷人員活動(dòng)空間中指標(biāo)PFR和ART的改善率ηi與建筑高度比的關(guān)系.

圖8 H1/H2對(duì)街谷人員活動(dòng)空間內(nèi)各項(xiàng)指標(biāo)改善率的影響Fig.8 Effect of ratios H1/H2on percentage decrease of the indices in the pedestrian domains

從圖8可以看到，結(jié)合城市主導(dǎo)風(fēng)向形成上階梯式街谷，對(duì)街道內(nèi)空氣品質(zhì)改善作用大于下階梯式街谷.上階梯式街谷兩側(cè)建筑的最佳比例H1/H2為3/7，此時(shí)空氣品質(zhì)改善率高達(dá)70～92%.據(jù)此，城市規(guī)劃設(shè)計(jì)人員可以根據(jù)實(shí)際環(huán)境盡可能地避免出現(xiàn)最不利街谷設(shè)計(jì)，并明確由此對(duì)街谷空氣品質(zhì)的改善程度.

4 結(jié)論

利用三維數(shù)值模擬分別研究了行列式和錯(cuò)列式建筑布局中上、下游建筑高度比例對(duì)街谷內(nèi)和下游建筑尾流區(qū)污染物分布的影響，并引入3種原本用于評(píng)價(jià)室內(nèi)空氣質(zhì)量的指標(biāo)來描述街谷人員活動(dòng)空間內(nèi)的空氣質(zhì)量.

4.1 街谷不對(duì)稱性對(duì)街谷內(nèi)部和下游建筑尾流區(qū)污染物濃度分布均有明顯影響.下游建筑越低矮，下游建筑尾流區(qū)濃度越高.

4.2 盡管平均殘留時(shí)間ART和有效通風(fēng)率PFR是為評(píng)價(jià)對(duì)稱街谷空氣品質(zhì)而提出的，但本文計(jì)算結(jié)果表明，在評(píng)價(jià)不對(duì)稱街谷空氣品質(zhì)時(shí)，平均殘留時(shí)間ART和有效通風(fēng)率PFR是適用的.

4.3 當(dāng)街谷兩側(cè)為多層建筑時(shí)，行列式和錯(cuò)列式街谷中上下游建筑高度比例分別為7/3、7/2時(shí)為最不利于機(jī)動(dòng)車尾氣擴(kuò)散的街谷尺寸，相對(duì)于最不利比例街谷，上階梯式街谷的空氣品質(zhì)改善率高于下階梯式街谷.為便于城市規(guī)劃中盡可能避免出現(xiàn)最不利比例街谷，本文給出了上下游建筑高度比與最不利街谷空氣品質(zhì)間的關(guān)系.

［1］Huang Y， Zhou Z. A numerical study of airflow and pollutant dispersion inside an urban street canyon containing an elevated expressway ［J］. Environmental Modeling and Assessment，2013，18（1）：105-114.

［2］Wang Y， Zhong K， Zhang N， et al. Numerical analysis of solar radiation effects on the flow patterns in street canyons ［J］. Engineering Applications of Computational Fluid Mechanics，2014，8（2）：252－262.

［3］Ng W， Chau C. A modeling investigation of the impact of street and building configuration on personal air pollutant exposure in isolated deep urban canyons ［J］. Science of the Total Environment， 2014，468-469（15）：429-448.

［4］Vardoulakis S， Fisher B E A， Pericleous K， et al. Modelling air quality in street canyons： a review ［J］. Atmospheric Environment，2003，37（2）：155-182.

［5］胡 偉，鐘 秦.壁面加熱作用對(duì)街道峽谷污染物擴(kuò)散的影響［J］. 中國環(huán)境科學(xué)， 2009，29（9）：908-913.

［6］王 樂，張?jiān)苽ィ櫿琢?動(dòng)態(tài)風(fēng)場(chǎng)及交通流量下街道峽谷內(nèi)污染物擴(kuò)散模擬 ［J］. 中國環(huán)境科學(xué)， 2012，32（12）：2161-2167.

［7］邵 毅，張述文，左洪超，等.榆中縣城大氣污染物擴(kuò)散的觀測(cè)實(shí)驗(yàn)與數(shù)值模擬 ［J］. 中國環(huán)境科學(xué)， 2013，33（7）：1223-1230.

［8］Xie X， Huang Z， Wang J. The impact of urban street layout on local atmospheric environment ［J］. Building and Environment，2006，41（10）：1352-1363.

［9］何澤能，高陽華，李永華，等.城市街谷峽谷對(duì)稱性對(duì)內(nèi)部氣流場(chǎng)的影響研究 ［J］. 氣象與環(huán)境學(xué)報(bào)， 2008，24（2）：62-67.

［10］Hang J， Li Y， Sandberg M， et al. The influence of building height variability on pollutant dispersion and pedestrian ventilation in idealized high-rise urban areas ［J］. Building and Environment，2012，56：346-360.

［11］汪立敏，王嘉松，黃 震，等.優(yōu)選數(shù)值計(jì)算方法用于不對(duì)稱街道峽谷的研究 ［J］. 水動(dòng)力學(xué)研究與進(jìn)展， 2005，20（1）：1-7.

［12］羅昔聯(lián)，顧兆林.基于DPM模型的街谷內(nèi)顆粒物擴(kuò)散特性研究［J］. 中國科學(xué)院研究生院學(xué)報(bào)， 2007，24（5）：578-583.

［13］Liu C， Wong C C C. On the pollutant removal， dispersion， and entrainment over two-dimensional idealized street canyons ［J］. Atmospheric Research， 2014，135-136：128-142.

［14］Gromke C， Buccolieri R， Sabatino S D， et al. Dispersion study in a street canyon with tree planting by means of wind tunnel and numerical investigations-Evaluation of CFD data with experimental data ［J］. Atmospheric Environment， 2008，42（37）：8640-8650.

［15］Nazridoust K， Ahmadi G. Airflow and pollutant transport in street canyons ［J］. Journal of Wind Engineering and Industrial Aerodynamics， 2006，94（6）：491-522.

［16］Barratt R. Atmospheric dispersion modelling - an introduction to practical applications ［M］. London： Earthscan Publication， 2001.

［17］GB50736-2012 民用建筑供暖通風(fēng)與空氣調(diào)節(jié)設(shè)計(jì)規(guī)范 ［S］.

［18］Bird R B， Steuart W E. Transport Phenomena ［M］. 2ndEd. New York： John Wiley and Sons， Inc， 2002.

［19］Li X， Leung D Y C， Liu C H， et al. Physical modeling of flow field inside urban street canyons ［J］. Journal of Applied Meteorology and Climatology， 2008，47：2058-2067.

［20］Meroney R N， Pavageau M， Rafadalis S， et al. Study of line source characteristics for 2-D physical modelling of pollutant dispersion in street canyon ［J］. Journal of Wind Engineering and Industrial Aerodynamics， 1996，62（1）：37-56.

［21］Liu C， Barth M C， Leung D Y C. Large-eddy simulation of flow and pollutant transport in street canyons of different buildingheight-to-street-width ratios ［J］. Journal of Applied Meteorology，2004，43（10）：1410-1424.

［22］Kato S， Ito K， Murakami S. Analysis of visitation frequency through particle tracking method based on LES and model experiment ［J］. Indoor Air， 2003，13（2）：182-193.

Effects of street geometric configurations on the pollutant dispersion around the canyons.

YANG Fang， ZHONG Ke，KANG Yan-ming*（School of Environmental Science and Engineering， Donghua University， Shanghai 201620， China）. China Environmental Science， 2015，35（3）：706～713

The effect of street geometric configurations， i.e.， building height layout， on pollutant distributions in both regular and staggered street canyons and the wake region of the downstream buildings are estimated by using the computational fluid dynamic （CFD） simulations. Three ventilation efficiency indices for evaluating indoor environment are applied to assess the air quality in the pedestrian domain. The numerical results show that average residence time（ART） and purging flow rate （PFR） can be used to evaluate the air quality inside the street canyons. The flow field and pollutant distribution inside and outside the canyon are greatly influenced by the building height layout， and the pollutant concentration in the wake region increases with decreasing the downstream building height. The highest values of the average concentration may occur in the pedestrian domain when the height ratios of the upstream building to the downstream building， i.e.， H1/H2， are 7/3 and 7/2 for regular and staggered canyons， respectively. In order to avoid these unfavorable building height layouts during urban planning， the relationships between H1/H2and the percentage reductions in PFR and ART are analyzed for both the regular and staggered arrangements.

asymmetrical street canyon；height ratio of upstream building to downstream building；worst building height layout；purging flow rate；average residence time

X51， P463， TU119， TU984

A

1000-6923（2015）03-0706-08

楊 方（1988-），女，江蘇連云港人，博士研究生，主要從事城市大氣環(huán)境與室內(nèi)空氣品質(zhì)研究.發(fā)表論文4篇.

2014-07-20

國家自然科學(xué)基金（40975093，41275157）；上海市教委科研創(chuàng)新重點(diǎn)項(xiàng)目（14ZZ073）

* 責(zé)任作者， 教授， ymkang@dhu.edu.cn