吳 培 寧
(中國勞動關(guān)系學院計算機網(wǎng)絡(luò)應(yīng)用研究所 北京 100048)
隨著城市交通流量的增加與道路建設(shè)的加快,城市快速路、高架復合道路以及高速公路等交通干道的建設(shè)越來越多,交通噪聲給周邊環(huán)境特別是道路兩側(cè)的高層建筑帶來了嚴重的噪聲污染,是亟待解決的城市主要環(huán)境問題之一。GIS系統(tǒng)具有強大的空間數(shù)據(jù)處理和顯示功能[1],可對空間數(shù)據(jù)及相關(guān)屬性信息進行有效的管理,近年來在社會需求的推動下飛速發(fā)展[2-3]?;贕IS的二維環(huán)評可視化系統(tǒng)在大氣環(huán)境與地面水環(huán)境環(huán)評中已取得了良好的應(yīng)用效果[4-5]。
目前,基于GIS的二維評價區(qū)域交通噪聲環(huán)評可視化是研究的熱點[6-9]。文獻[7-8]分別實現(xiàn)了無路基聲影區(qū)的城市地面道路與高速公路地面路段的交通噪聲模擬計算與可視化。文獻[9]通過MapInfo系統(tǒng)二次開發(fā),實現(xiàn)了有路基聲影區(qū)的高架道路交通噪聲環(huán)評可視化,但仍存在以下問題:1) 僅生成了評價區(qū)域地表平面處二維笛卡爾網(wǎng)格;2) 采用的交通噪聲環(huán)評計算模型中未考慮聲屏障的衰減量計算,無法用于有聲屏障路段的環(huán)評計算與可視化分析。
聲屏障已成為城市道路噪聲降噪的主要手段[10],各種聲屏障插入衰減量計算模型以及新型聲屏障材質(zhì)與結(jié)構(gòu)的研究已成為研究的熱點與難點問題[11-15]。但目前對聲屏障降噪的預(yù)測計算僅局限于二維評價區(qū)域內(nèi)離散點的噪聲預(yù)測計算,適用于區(qū)域噪聲環(huán)評可視化的三維大數(shù)據(jù)場預(yù)測計算目前還未實現(xiàn)。
基于GIS的聲屏障路段三維交通噪聲環(huán)評可視化能夠為環(huán)評提供科學依據(jù)、基礎(chǔ)資料與決策支持,對于高架道路聲屏障的規(guī)劃與建設(shè)也具有十分重要的指導意義。針對文獻[9]研究的不足以及聲屏障路段三維交通噪聲環(huán)評可視化在網(wǎng)格剖分、聲影區(qū)與聲照區(qū)判定、繞射聲衰減量計算及噪聲數(shù)據(jù)場可視化等方面的難點問題,本文實現(xiàn)了以下創(chuàng)新:1) 提出了基于自適應(yīng)分層切割單元法的交通噪聲三維評價區(qū)域網(wǎng)格生成方法。2) 提出了各分層內(nèi)聲照區(qū)長度計算公式,聲影區(qū)與聲照區(qū)判定及其笛卡爾網(wǎng)格生成方法。3) 提出了滿足區(qū)域環(huán)評計算要求的聲屏障插入衰減量計算方法。4) 采用了分層等值線可視化顯示方法。
應(yīng)用上述方法開發(fā)了基于GIS的三維交通噪聲屏障路段環(huán)評可視化系統(tǒng)。
文獻[16]規(guī)定地面評價范圍為道路中心線外兩側(cè)200m以內(nèi),考慮到高架道路本身的高度以及地面反射對高處噪聲的增強作用,本文將評價區(qū)水平寬度設(shè)定為隨高度增加而增加,道路三維評價區(qū)域內(nèi)距地面任一高度處的寬度由下式確定:
(1)
式中:Dh為距地面高度為h(h=0,1,…,h評)處,道路中心線至評價區(qū)域一側(cè)邊界的垂直距離;B為道路寬度,H為道路中心線距地面高度;θ(0<θ<90°)為評價區(qū)域外邊界面與地面的夾角,其取值與評價道路交通流量逆相關(guān)。
規(guī)劃道路中心線段集由下式表示:
R={(x0,y0,H0),…,(xi,yi,Hi),…,
(xn,yn,Hn)}
(2)
式中:R為規(guī)劃道路的中心線段集;(xi,yi)為某一道路中心線段的端點;Hi為該端點距地面高度。
交通噪聲三維評價區(qū)域為沿道路兩側(cè)的兩組梯形多面體。如圖1所示,其道路一側(cè)三維評價區(qū)域多面體的高度為評價區(qū)域高度h評,該多面體的內(nèi)邊界為沿道路一側(cè)邊界線段且垂直于地面的一組平面,外邊界為距道路中心線Dh處與地面夾角為θ的一組平面。
圖1 高架道路YOZ平面噪音評價區(qū)域
如圖2所示,沿道路中心線前進方向,其道路左側(cè)及右側(cè)邊界線段集由Rl及Rr表示,則:
圖2 高架道路XOY平面噪音評價區(qū)域
與道路中心線位于同一水平高度H的評價區(qū)域邊界線段集由RL及RR表示,則:
在距地面高度h處(h=0,1,…,h評)且平行于XOY平面的平面上,沿順時針方向排列的左側(cè)評價區(qū)域邊界線端點集由CL,h表示,右側(cè)評價區(qū)域邊界線端點集由CR,h表示。
三維評價區(qū)域環(huán)評計算網(wǎng)格生成分為兩個步驟:1) 采用一組切割平面切割評價區(qū)域多面體;2) 在各切割平面內(nèi)采用笛卡爾切割單元法生成二維笛卡爾網(wǎng)格。
對于復雜道路系統(tǒng),若采用固定層厚切割平面族,則無法兼顧環(huán)評計算的精度與效率。因此本文采用了基于道路幾何特征的自適應(yīng)變層厚技術(shù),對于水平彎曲系數(shù)及縱坡坡度變化較小的評價道路采用較大的切割層厚度,反之則采用較小的層厚。其算法如下:
步驟1 對于由式(2)表示的規(guī)劃道路,計算其水平彎曲系數(shù)k及縱坡坡度η。
(3)
(4)
步驟2 對于切割平面族,先取一個固定層厚Zh,則自適應(yīng)層厚為:
(5)
由評價區(qū)地面至h評,以ΔZh為層厚生成一組切割平面,每個平面與評價區(qū)域多面體求交,得到兩個評價多邊形。對于第j個切割平面,其沿道路中心線左側(cè)評價多邊形端點集為:
沿道路中心線右側(cè)評價多邊形端點集為:
步驟3 在第j個切割平面內(nèi),采用笛卡爾背景網(wǎng)格與其沿道路中心線左右側(cè)評價多邊形端點集CL,j與CR,j中的每一條邊界線段逐一求交,求交過程為切割單元交點追蹤過程。當完成交點追蹤后,將所有與CL,j與CR,j有交點的背景網(wǎng)格單元標記為切割網(wǎng)格單元,存入切割單元鏈表中。此算法的時間復雜度為O(n2)。
切平面j內(nèi)的背景網(wǎng)格完成與端點集CL,j與CR,j的切割后,被分為三類:1) 評價區(qū)域內(nèi)網(wǎng)格;2) 評價區(qū)域外網(wǎng)格;3) 邊界線切割網(wǎng)格。其中:評價區(qū)域外網(wǎng)格作為無用網(wǎng)格需刪除;對于邊界線切割網(wǎng)格,采用碎片切割單元合并算法將其與相鄰切割單元合并。得到兩種網(wǎng)格單元后,結(jié)合當前切片層層厚信息即得到了第j層三維笛卡爾噪聲環(huán)評計算網(wǎng)格。分別存儲于左側(cè)、右側(cè)評價區(qū)域內(nèi)網(wǎng)格單元鏈表,左側(cè)、右側(cè)評價區(qū)域合并單元鏈表。4個鏈表均為雙向鏈表,其數(shù)據(jù)域由下式描述:
式中:(xi1,yi1,jΔZh),(xi2,yi2,jΔZh),(xi3,yi3,jΔZh),(xi4,yi4,jΔZh),(xi1,yi1,(j+1)ΔZh),(xi2,yi2,(j+1)ΔZh),(xi3,yi3,(j+1)ΔZh),(xi4,yi4,(j+1)ΔZh)為三維笛卡爾網(wǎng)格單元8個節(jié)點坐標;(a1,a2,a3,a4,a5,a6,a7,a8)為節(jié)點屬性,其屬性分為邊界線段切割交點與背景網(wǎng)格點兩類。
聲屏障是降低地面運輸噪音的有效措施之一。依據(jù)文獻[17],安裝聲屏障前后在某特定位置上的聲壓級之差為聲屏障插入衰減量Abar。
Abar=ΔLd-ΔLd有限-ΔLt-ΔLr-(ΔLS,ΔLG)max
(6)
式中:ΔLd為繞射聲衰減;ΔLt為透射聲修正量;ΔLr為反射聲修正量;ΔLS為障礙物聲衰減;ΔLG為地面吸收聲衰減;max表示取ΔLS與ΔLG中的最大值。單位均為dB(A)。
本文的研究對象為僅建有單側(cè)聲屏障的道路,因此反射聲修正量ΔLr為零。
當聲源為無限長聲源且聲屏障為無限長時,其繞射聲衰減由下式計算:
(7)
(8)
式中:f為聲波頻率,單位為Hz;δ=a+b-d為繞射路徑聲程差(如圖3所示)單位為m;c為聲速,單位為m/s。
圖3 交通噪聲繞射及透射路徑
當聲源為一無限長聲源而聲屏障為有限長時,其繞射聲衰減仍由式(7)、式(8)計算,但需要減去一個修正量,則:
(9)
步驟1 對于聲屏障路段聲影區(qū)內(nèi)的每一網(wǎng)格節(jié)點,若該點處聲屏障遮蔽角百分率大于70%,則進入步驟2 。
步驟3 若該線段與9條曲線皆無交點,則該點處的修正量為2 dB(A)。若該線段與一條或多條曲線段有交點,則選取交點序號N最高的曲線為其取值曲線。該點處的修正量為:
計算出繞射聲衰減ΔLd后,透射聲修正量ΔLt、地面吸收聲衰減ΔLG可由文獻[17]中的公式計算得到,代入式(6)即得到了聲屏障插入衰減量Abar。
得到聲屏障衰減量Abar后,采用文獻[9]中的算法即可計算得到道路總車流聲效等級,單位為dB(A)。
對于無聲屏障路段、有聲屏障路段的評價區(qū)域內(nèi)網(wǎng)格,均需判斷哪些為聲影區(qū)網(wǎng)格,哪些為聲照區(qū)網(wǎng)格,再依據(jù)相應(yīng)的噪聲模型完成計算。對于圖3所示聲屏障路段,本文提出了基于道路及聲屏障幾何參數(shù)的各分層內(nèi)聲影區(qū)與聲照區(qū)笛卡爾網(wǎng)格生成算法。
步驟1 人機交互輸入道路幾何參數(shù),包括:路基高度H、路面寬度B、路肩高度h肩、聲屏障起始點坐標[(xbs,ybs)、(xbe,ybe)]、聲屏障有效高度hb等參數(shù),其等效車道高度為距路面1.2 m。將輸入?yún)?shù)存入道路幾何參數(shù)表中。
步驟2 以聲屏障位于沿中心線前進方向右側(cè)為例,計算聲屏障造成的聲影區(qū)長度,沿順時針方向排列的XOY平面內(nèi)聲屏障影響區(qū)域邊界線段集由Cb表示:
讀取1.3節(jié)中生成的右側(cè)評價區(qū)域內(nèi)網(wǎng)格單元鏈表、右側(cè)評價區(qū)域合并單元鏈表,判斷上述兩個鏈表中被Cb包圍的網(wǎng)格節(jié)點,將其存儲于新的聲屏障影響區(qū)網(wǎng)格鏈表中。
(10)
對于無聲屏障的高路堤及高架路段,基于道路幾何參數(shù)的各切割平面內(nèi)路基聲影區(qū)長度為:
(11)
采用與聲屏障路段類似的算法,可生成左側(cè)路基聲影區(qū)、左側(cè)路基聲照區(qū)、右側(cè)路基聲影區(qū)及右側(cè)路基聲照區(qū)4個網(wǎng)格單元鏈表。
評價區(qū)內(nèi)三維交通噪聲預(yù)測計算由各分層內(nèi)的二維噪聲環(huán)評計算組成。對于有聲屏障路段,采用了基于分層的評價區(qū)三維交通噪聲環(huán)評算法。
步驟1 人機交互輸入環(huán)評參數(shù),將參數(shù)存入環(huán)評參數(shù)表中。
步驟2 依次遍歷左側(cè)路基聲照區(qū)、左側(cè)路基聲影區(qū)、右側(cè)路基聲影區(qū)、右側(cè)路基聲照區(qū)、聲屏障聲影區(qū)、聲屏障聲照區(qū)等網(wǎng)格單元鏈表。
針對GIS系統(tǒng)的特點,對評價區(qū)域三維交通噪聲環(huán)評數(shù)據(jù)場,本文采用了分層等值線可視化顯示。
2.4節(jié)生成的三個區(qū)域噪聲值鏈表為噪聲等值線的繪制提供了數(shù)據(jù)基礎(chǔ),每個鏈表均為雙向鏈表,其結(jié)點數(shù)據(jù)域由下式描述:
Si,j={{(xi,j,yi,j,jΔZh),i=1,2,…,n,j=0,1,…,n},a,b,Leq}
式中:(xi,j,yi,j,jΔZh)為網(wǎng)格節(jié)點Si,j的坐標;a為網(wǎng)格節(jié)點背景屬性,分為背景網(wǎng)格點與邊界線段切割交點兩類;b為網(wǎng)格節(jié)點噪聲屬性,分為路基聲影區(qū)、路基聲照區(qū)、聲屏障聲照區(qū)、無線長聲屏障聲影區(qū)、有限長聲屏障聲影區(qū)五類;Leq為Si,j節(jié)點處的總車流等效聲級。
本文提出了基于分層網(wǎng)格序列法的等值點追蹤算法。
步驟1 依次遍歷左側(cè)無聲屏障評價區(qū)域噪聲值鏈表、右側(cè)無聲屏障評價區(qū)域噪聲值鏈表及單側(cè)聲屏障影響區(qū)域噪聲值鏈表,對于任一鏈表中分層高度等于jΔZh(j=0,1,…,n)的所有結(jié)點執(zhí)行下述步驟。
步驟2 按照先X軸后Y軸的順序遍歷,對于每一鏈表中的結(jié)點依照其結(jié)點背景屬性及相鄰節(jié)點,構(gòu)造合并網(wǎng)格單元或背景網(wǎng)格單元。對于節(jié)點背景屬性a為邊界線段切割交點的結(jié)點,構(gòu)造如圖4所示合并網(wǎng)格單元。對于節(jié)點背景屬性a為背景網(wǎng)格點的結(jié)點,構(gòu)造如圖5所示正方形背景網(wǎng)格單元。對于每一單元的頂點,比較其Leq值與給定等值線值的大小,Leq值大于等值線值的標記為“+”頂點,反之為“-”頂點,則:1) 只有一個“+”頂點或“-”頂點的單元有兩個等值點;2) 有兩個“+”頂點與兩個“-”頂點的單元,其等值點個數(shù)為2個或4個。對于由“+”頂點與“-”頂點共同組成的單元邊,采用線性插值計算等值線與該條邊的交點,將交點坐標按所屬鏈表存儲于GIS系統(tǒng)中的j層左側(cè)評價區(qū)等值線TAB表、j層右側(cè)評價區(qū)等值線TAB表或j層聲屏障影響區(qū)等值線TAB表。
圖4 合并單元網(wǎng)格等值點追蹤
圖5 背景單元網(wǎng)格等值點追蹤
步驟3 若背景網(wǎng)格單元頂點Leq值分布如圖5(d)所示,單元內(nèi)必有四個等值點,則評價區(qū)域內(nèi)同一個等值線值存在著兩條等值線。因此將該單元剖分為四個三角形單元,在四個三角內(nèi)進行等值點追蹤。新增網(wǎng)格中心點處的Leq值為單元四個頂點處Leq的平均值。圖6 (a)中心點處為“-”頂點,圖6 (b)中心點處為“+”頂點。采用線性插值計算等值線與剖分單元8條邊的交點,可得到兩組共6個等值點。再通過判斷與哪條等值線有共享等值點,將兩組等值點的坐標分別存入擁有相同等值線值的兩條等值線TAB表中。
圖6 剖分單元等值點追蹤
步驟4 重復執(zhí)行步驟1,直至完成j個平面上的等值點追蹤,分別生成左側(cè)評價區(qū)等值線TAB表、右側(cè)評價區(qū)等值線TAB及或聲屏障影響區(qū)等值線TAB表三類等值線表。此算法的時間復雜度為O(mn3)。
噪聲等值點追蹤算法為噪聲等值線繪制及等值線拓撲填充提供了基礎(chǔ),采用三次樣條函數(shù)對存儲于三類TAB的等值點數(shù)據(jù)進行光滑,將求得的插值點坐標存入光滑TAB表中??梢暬@示階段,在GIS系統(tǒng)中創(chuàng)建j個圖層,讀取與每個圖層高度相對應(yīng)的三類光滑TAB表,并用直線段將插值點依次連接起來即得到了光滑后的等值線。
等值線TAB表中某一等值線值的等值點連接后為評價區(qū)域內(nèi)一條非閉合等值線,將其沿道路邊界線段及評價區(qū)域邊界線段進行封閉后形成等值線區(qū)域,判斷各條等值線區(qū)域的包含與并列關(guān)系后可生成等值線拓撲關(guān)系二叉樹。前序遍歷該二叉樹得到一個填充序列,按此序列采用不同顏色依次填充各條等值線內(nèi)的區(qū)域即可得到等值線拓撲填充圖。
為驗證系統(tǒng)三維交通噪聲環(huán)評計算的準確性,將計算結(jié)果與交通噪聲實際測量數(shù)據(jù)進行了對比分析。選取京港澳高速一段單側(cè)有聲屏障路段為實驗道路,測量地點為一棟建設(shè)中的高樓,測量處距路基水平距離63 m。為便于通過視頻辨識與統(tǒng)計交通流量,監(jiān)測時間選擇為晝間。監(jiān)測儀器為HT-8352聲級計,高度測量采用HT-308激光測距儀,長度測量采用JY-361999測距輪。測量時間為30 min,并同時監(jiān)測交通流量,然后再進行統(tǒng)計。表1為聲屏障路段噪聲垂直分布實際監(jiān)測結(jié)果與環(huán)評計算結(jié)果的比較。表中數(shù)據(jù)顯示:低樓層處計算值高于實際監(jiān)測值。原因為該路段已建成二十多年,兩側(cè)生長的已成林綠化帶降噪效果好于初建成的高速公路兩側(cè)綠化帶。高樓層處環(huán)評計算值低于實際監(jiān)測值,原因為監(jiān)測當天曾有降雨,路面較潮濕導致車輛輪胎噪聲增大,路面噪聲修正量不為零。
表1 聲屏障路段噪聲垂直分布監(jiān)測結(jié)果與環(huán)評計算結(jié)果比較
應(yīng)用本文提出的算法,在MapInfo桌面GIS系統(tǒng)上開發(fā)了三維交通噪聲屏障路段環(huán)評可視化系統(tǒng),編程語言采用該系統(tǒng)自帶的二次開發(fā)語言MapBasic。系統(tǒng)流程如圖7所示,可視化顯示技術(shù)路線如圖8所示。
圖7 系統(tǒng)流程圖
圖8 可視化路線圖
如圖9所示,道路寬度B=60 m,高度H=5.8 m,采用式(3)、式(4)計算水平彎曲系數(shù)k及縱坡坡度η后,由式(5)得到切片層厚度ΔZh=3 m,道路前進方向左側(cè)有一座110 kV變電站,周圍為非居住區(qū)。前進方向右側(cè)為居住區(qū),建有高度為3.6 m的聲屏障,此路段為單側(cè)無限長聲屏障路段。高度jΔZh=3 m分層內(nèi),評價區(qū)域邊界距道路中心線距離由式(1)計算得D3=283.93 m,各條等值線值由MapInfo系統(tǒng)以信息窗口形式輸出,等效聲級大于等于70 dB(A)的等值線采用紅色繪制(圖中以灰度表示),小于70 dB(A)的等值線顏色保持系統(tǒng)默認值。因3 m分層高度小于道路與聲屏障高度之和,聲屏障一側(cè)評價區(qū)域全部為聲屏障聲影區(qū),等效聲級等值線的最大值為69 dB(A),最小值為56 dB(A),符合4a類聲環(huán)境功能區(qū)的環(huán)境噪聲限值[18]。無聲屏障一側(cè)的評價區(qū)域,由式(11)可得該側(cè)的路基聲影區(qū)長度為46.7 m,路基聲影區(qū)內(nèi)的等值線最大值為68 dB(A)。垂直距離大于該長度的路基聲照區(qū)內(nèi)等值線的最大值為78 dB(A),已超出環(huán)境噪聲限值。
圖9 無限長聲屏障路段高度3 m分層內(nèi)噪聲等效聲級等值線
圖10 無限長聲屏障路段高度30 m分層內(nèi)噪聲等效聲級等值線
圖11為單側(cè)有限長聲屏障路段,該路段為互通立交樞紐,寬度B=80 m,道路高度H=3.8 m,依式(3)-式(5)得到切片層厚度ΔZh=2 m。圖中的環(huán)境敏感點用五角星標注,聲屏障一側(cè)側(cè)評價區(qū)域從335 m處開始的主干道及彎道均建有高度為3.6 m的聲屏障。高度jΔZh=6 m分層內(nèi),評價區(qū)域邊界距道路中心線距離由式(1)計算得D6=278.65 m,由式(11)可知,評價區(qū)域內(nèi)無聲屏障的路段全部為路基聲照區(qū),其等效聲級等值線最大值為85 dB(A)。因6 m分層高度小于道路與聲屏障高度之和,此層內(nèi)有聲屏障的路段全部為聲屏障聲影區(qū),等值線最大值為70 dB(A)。
圖11 有限長聲屏障路段高度6 m分層內(nèi)噪聲等效聲級等值線
圖12為將圖11中6條等值線沿邊界線段封閉并建立拓撲關(guān)系二叉樹后,等值線拓撲填充圖與背景地圖的疊加顯示。
圖12 有限長聲屏障路段高度6 m分層內(nèi)等值線拓撲填充圖
圖13為6 m分層高度內(nèi)評價區(qū)域笛卡爾網(wǎng)格單元節(jié)點與背景地圖疊加顯示,從左至右的網(wǎng)格節(jié)點噪聲屬性b分別為無線長聲屏障聲影區(qū)、有限長聲屏障路段聲影區(qū)及路基聲照區(qū)共三類網(wǎng)格單元節(jié)點,節(jié)點間距20 m,采用顏色漸變加以區(qū)分。其中有限長聲屏障路段聲影區(qū)網(wǎng)格單元節(jié)點的遮蔽角百分率大于70%、小于90%。
圖13 有限長聲屏障路段高度6 m分層內(nèi)計算網(wǎng)格單元節(jié)點
圖14 有限長聲屏障路段高度12 m分層內(nèi)計算網(wǎng)格單元節(jié)點
本文首次實現(xiàn)了基于GIS的聲屏障路段三維交通噪聲環(huán)評可視化。實驗結(jié)果表明,基于GIS的三維交通噪聲環(huán)評可視化能夠準確直觀地反映道路周邊評價區(qū)域內(nèi)噪聲的空間分布特征,提高了交通噪聲環(huán)境影響評價的準確性與效率。結(jié)合GIS系統(tǒng)獨特的空間與屬性數(shù)據(jù)管理功能,對經(jīng)過醫(yī)院、學校、住宅等噪聲敏感建筑物集中區(qū)域的高速公路與城市高架道路兩側(cè)聲屏障的設(shè)計與安裝具有非常重要的指導意義。通過對比分析,系統(tǒng)采用的聲屏障衰減量計算模型較準確,環(huán)評預(yù)測計算值與實測噪聲值之間的差值在2~4 dB(A)以內(nèi)。下一步將研究道路兩側(cè)已成林綠化帶的吸收聲衰減以及潮濕瀝青路面的噪聲修正量。