姚 博

(中國十九冶集團有限公司， 四川成都 611730)

1 概述

東南沿海地區(qū)臺風(fēng)頻發(fā)，無論是高層建筑還是大跨結(jié)構(gòu)建造設(shè)計，其極值風(fēng)速確定較為困難，如果考慮所處地形影響，則極值風(fēng)速分析更為復(fù)雜。傳統(tǒng)極值風(fēng)速分析方法用于建筑極值風(fēng)速分析時，對復(fù)雜地形很難適用，無法直接確定設(shè)計基準(zhǔn)風(fēng)速；如果仍然采用傳統(tǒng)方法在所有的目標(biāo)點都建設(shè)觀測站，則施工費用極大，兩年及其以上的風(fēng)氣候觀測周期也是難以承受的時間成本。因此如何考慮到復(fù)雜地形的影響，精確給出準(zhǔn)確的設(shè)計基準(zhǔn)風(fēng)速對于建筑結(jié)構(gòu)設(shè)計安全也具有十分重要的意義。

在復(fù)雜地形的數(shù)值模擬研究方面。Giebel等[1]對復(fù)雜地形采用最高精度的地形高程數(shù)據(jù)建立網(wǎng)格模型，CFD的計算結(jié)果表明，高精度的地形高程數(shù)據(jù)可以顯著改善分析結(jié)果。Palma等人[2]采用10m精度的地形數(shù)據(jù)，建立高精度的模型網(wǎng)格，研究了復(fù)雜地區(qū)的風(fēng)場特性，結(jié)果表明，非線性的CFD模型計算結(jié)果與實測資料吻合。由上述分析可見，采用CFD來模擬實地風(fēng)場的一個前提條件就是獲得足夠精度的地形數(shù)據(jù)。

對設(shè)計基準(zhǔn)風(fēng)速分析方法的研究。同濟大學(xué)朱樂東[3]在研究壩陵河大橋設(shè)計基準(zhǔn)風(fēng)速時，通過在橋位建造位置安裝觀測塔和風(fēng)速觀測儀器，采集兩年內(nèi)橋位風(fēng)速實測數(shù)據(jù)，建立了與附近氣象站近的風(fēng)速關(guān)系，進而由氣象站長期記錄估算的極值風(fēng)速推算橋位處的極值風(fēng)速作為設(shè)計風(fēng)速，這種方法彌補了橋位位置缺少觀測資料的缺陷，給橋梁建設(shè)提供了較為可靠的理論保障。目前對山區(qū)風(fēng)特性的研究大多數(shù)以Jackson的線性化理論分析為框架，再此基礎(chǔ)上通過實測、風(fēng)洞模擬和數(shù)值模擬等方式，推動了相關(guān)的研究。但存在一個不足就是分析目標(biāo)固定，時間周期較長。

綜上可知，傳統(tǒng)極值方法在分析混合氣候極值風(fēng)速時存在諸多不足。本文通過CFD計算建立氣象站和目標(biāo)點之間的的風(fēng)速比模型，獲取了目標(biāo)位置的極值風(fēng)速序列，基于混合函數(shù)構(gòu)建臺風(fēng)和良態(tài)風(fēng)混合氣候的極值風(fēng)速概率密度函數(shù)，最后采用蒙特卡洛模擬等方法驗證了本文方法的有效性。與傳統(tǒng)方法相比，本文方法能同時考慮復(fù)雜地形影響和混合氣候極值分布特性尾部偏離影響，過程簡單精確，容易被工程實踐采納。

2 地貌選取和風(fēng)速序列生成

2.1 復(fù)雜地形的模型生成和風(fēng)速比計算

復(fù)雜山地模型的生成包括兩個方面的內(nèi)容，即模型的生成和網(wǎng)格的自動化劃分。建立用于分析復(fù)雜目標(biāo)地點CFD模型的地形高程數(shù)據(jù)從GIS地理空間數(shù)據(jù)云下載，數(shù)據(jù)的精度為30×30的DEM數(shù)據(jù)(柵格影像格式)。對下載數(shù)據(jù)首先采用ArcMap軟件進行處理，轉(zhuǎn)化成高程數(shù)據(jù)點，緊接著對高程數(shù)據(jù)點添加實地的(X,Y)坐標(biāo)，這樣就得到了目標(biāo)點真實的地面坐標(biāo)-高度信息，最后將數(shù)據(jù)輸出，并進一步按照建模要求和規(guī)則，選取目標(biāo)地附近適宜區(qū)域進行分析。

建模和網(wǎng)格劃分時，以目標(biāo)地為中心建立一圓形分析區(qū)域模型，模型直徑根據(jù)分析區(qū)域大小以及周邊氣象站的距離來確定，外環(huán)過渡區(qū)寬度取模型半徑的五分之一即可，由于關(guān)心的風(fēng)場特性集中于近地面，所以計算域高度取h=4000m，來流可能任意方向吹來。

風(fēng)速比的計算。將風(fēng)速比rm定義為第m個方向，一定來流位置10m高度和目標(biāo)點10m高度的風(fēng)速比值。計算各個方向的風(fēng)速比，得到風(fēng)速比序列集合{rm}。

2.2 氣象站選取和目標(biāo)點風(fēng)速序列生成

用于CFD計算的模型建立之后，選擇目標(biāo)點四周方向分布較均勻，且具有較長歷史記錄的氣象站，最后結(jié)合各個方向風(fēng)速比和該方向氣象站對應(yīng)的風(fēng)速風(fēng)向歷史記錄，模擬得到目標(biāo)地的極值風(fēng)速序列，具體實現(xiàn)過程為：

(1)方向定義。風(fēng)向定于與氣象站風(fēng)向記錄保持一致，即以正北方向為零度沿順時針劃分成k=16個方向，每個方向βk為22.5 °。

(2)目標(biāo)點風(fēng)速序列的生成。基于周邊氣象站的日風(fēng)速風(fēng)向記錄并結(jié)合風(fēng)速比得到目標(biāo)點的日極值風(fēng)速序列。假設(shè)共有M個氣象站點，風(fēng)速風(fēng)向歷史記錄為N年，在第n=1,2,……N×365天，第m=1,……M個氣象站點的極值風(fēng)速定義為Vn,m,k，其中k為實測記錄的風(fēng)向。極值風(fēng)速可由各氣象站當(dāng)日記錄的風(fēng)速風(fēng)向數(shù)據(jù)乘以對應(yīng)的風(fēng)速比得到，根據(jù)風(fēng)向和氣象站方向的不同，可分為以下情況：

②氣象站記錄風(fēng)向相對于目標(biāo)點的方向在θm=60°之外，則認(rèn)為對目標(biāo)點影響較弱。

2.3 數(shù)據(jù)處理和極值風(fēng)速計算

本文首先基于考慮方向的獨立風(fēng)暴方法進行取樣，再構(gòu)建混合氣候極值風(fēng)速的概率密度函數(shù)，然后基于加權(quán)最小二乘法求解參數(shù)，最終精確計算極值風(fēng)速。該實現(xiàn)過程為：

(1)采用考慮方向的獨立風(fēng)暴法進行取樣。通過選取一定的風(fēng)速閾值，并將超過該閾值的聯(lián)系風(fēng)速樣本記錄為一個獨立風(fēng)暴，取樣的時候取每個獨立風(fēng)暴內(nèi)各個風(fēng)向最大值作為極值樣本。由于每個獨立風(fēng)暴內(nèi)每個方向僅選擇一個樣本，這又保證了每個方向風(fēng)速的獨立性。

(2)確定極值風(fēng)速分布的混合概率密度函數(shù)。將風(fēng)速數(shù)據(jù)分為良態(tài)風(fēng)主導(dǎo)區(qū)和臺風(fēng)主導(dǎo)區(qū)，并分別擬合得到相應(yīng)的概率分布函數(shù)，再此基礎(chǔ)上進一步構(gòu)建如式(1)所示的用于描述混合氣候區(qū)極值風(fēng)速的概率密度函數(shù)，其表達式可以表述為：

(1)

(2)

Z(a,w,Π,σ,ξ)為歸一化函數(shù)，表達式為：

(3)

(4)

(4)一定重現(xiàn)期極值風(fēng)速的確定。由步驟(3)得到混合分布函數(shù)概率密度函數(shù)，進而積分得到其分布函數(shù)，從而得到50a、100a重現(xiàn)期極值風(fēng)速。

3 算例

為了更清楚的展示上述方法過程，這里以某復(fù)雜地形位置極值風(fēng)速的確定為例進行說明。

3.1 風(fēng)場建模

根據(jù)地形特點選取典型地貌區(qū)域，以經(jīng)緯度為中心坐標(biāo)，選取一定半徑L的中心圓區(qū)域和L0的外環(huán)過渡區(qū)進行建模。由于只關(guān)心近地風(fēng)場特性，所以計算域高度取4 000m，根據(jù)選擇分析域的大小選取合適的最大特征長度和加密區(qū)長度分別為100m和40m，粗糙層劃分為10層，加密區(qū)半徑取600m，加密過渡半徑取3 500m，縱向網(wǎng)格生成率為1.2，入口風(fēng)速取10m/s。模型鈍化、網(wǎng)格劃分如圖1所示。

圖1 內(nèi)部鈍化與網(wǎng)格劃分

3.2 風(fēng)速比計算和目標(biāo)點風(fēng)速序列生成

3.3 數(shù)據(jù)處理和極值風(fēng)速計算

結(jié)合各氣象站風(fēng)速風(fēng)向記錄和對應(yīng)氣象站點風(fēng)速比，計算得到目標(biāo)點極值風(fēng)速序列。對于3.2節(jié)得到的極值風(fēng)速，采用獨立風(fēng)暴法進行取樣，再基于2.3節(jié)方法原理對得到數(shù)據(jù)進行分析，得到其混合概率分布函數(shù)，最后計算N年一遇的設(shè)計極值風(fēng)速?；?.3節(jié)原理進行實現(xiàn)的過程如圖2和圖3所示。

圖2 混合模型的分布參數(shù)擬合流程

圖3 混合模型權(quán)函數(shù)的形狀參數(shù)求解流程

3.4 結(jié)果分析

采用上述分析放量對不同位置的兩處風(fēng)速數(shù)據(jù)進行分析，分析結(jié)果如圖4所示。

圖4 臺風(fēng)混合氣候極值風(fēng)速分析結(jié)果

由圖4可知，傳統(tǒng)極值I型擬合方法在對尾部分布考慮欠佳，即在低速良態(tài)風(fēng)主導(dǎo)范圍與實測結(jié)果吻合較好，而在重現(xiàn)期較高的臺風(fēng)主導(dǎo)范圍該方法不再適用。進一步的分析表明，傳統(tǒng)分析方法在處理良態(tài)風(fēng)氣候時具有足夠的計算精度，但是混合分布受尾部偏離的影響，此時仍采用傳統(tǒng)分析方法就有一定的局限性；本文方法與臺風(fēng)MonteCarlo模擬方法[5]和實測數(shù)據(jù)吻合較好，且可以很好的考慮混合極值樣本尾部偏離的現(xiàn)象。

4 結(jié)論

本文通過CFD建立風(fēng)場模型，得到目標(biāo)的極值風(fēng)速時間序列，再基于混合函數(shù)構(gòu)建臺風(fēng)和良態(tài)風(fēng)混合氣候地區(qū)的極值風(fēng)速概率密度函數(shù)。通過考慮復(fù)雜地形影響的臺風(fēng)氣候極值風(fēng)速的研究，得到如下結(jié)論：

(1)內(nèi)陸良態(tài)風(fēng)氣候區(qū)，采用傳統(tǒng)極值I型分布即可獲得精度滿意的計算結(jié)果，但臺風(fēng)混合氣候區(qū)，尾部分布偏離較為明顯，傳統(tǒng)的極值風(fēng)速分析方法會造成不同程度的誤差。

(2)傳統(tǒng)分析方法需要在目標(biāo)點設(shè)立氣象觀測站，建設(shè)成本較高，時間周期較長，本文方法不需要建設(shè)觀測站亦不需要觀測長周期的風(fēng)速樣本，僅需通過CFD數(shù)值模擬建立相關(guān)關(guān)系，因此綜合成本優(yōu)勢較大。

(3)傳統(tǒng)分析方法僅建立目標(biāo)地與附近最近氣象站的關(guān)系，而并不采用周圍其他站點的歷史記錄。本文方法利用了周邊所有的氣象站點數(shù)據(jù)，數(shù)據(jù)信息更豐富，利用率更高。