蔡彥楓，黃穗，黃增浩

（1.中國(guó)能源建設(shè)集團(tuán)廣東省電力設(shè)計(jì)研究院有限公司，廣州 510663；2.廣東科諾勘測(cè)工程有限公司，廣州 510663；3.南方電網(wǎng)科學(xué)研究院，廣州 510663）

0 引言

臺(tái)風(fēng)是造成電力設(shè)施損毀，威脅沿海輸電線路安全穩(wěn)定運(yùn)行的最主要自然災(zāi)害。以南方電網(wǎng)為例，自2003 年以來(lái)，110 kV 等級(jí)及以上的主干電網(wǎng)因登陸臺(tái)風(fēng)造成的倒塔、桿塔受損、斷線等的重大風(fēng)災(zāi)事故共發(fā)生10 起，其中2011—2015 年更是連續(xù)5 年每年出現(xiàn)1 起［1］，國(guó)民經(jīng)濟(jì)蒙受巨大損失。因此，防范和應(yīng)對(duì)臺(tái)風(fēng)引起的風(fēng)災(zāi)事故一直都是電網(wǎng)防災(zāi)減災(zāi)工作的重、難點(diǎn)［2-4］。為了準(zhǔn)確掌握臺(tái)風(fēng)的致災(zāi)機(jī)理，科學(xué)制定防災(zāi)減災(zāi)措施，有必要針對(duì)沿海輸電線路局地風(fēng)場(chǎng)開(kāi)展觀測(cè)、數(shù)值模擬等一系列基礎(chǔ)研究。

沿海輸電線路普遍采用架空線-塔形式，對(duì)風(fēng)荷載的變化較為敏感［5-6］。輸電線路風(fēng)荷載計(jì)算需要輸入的風(fēng)特性參數(shù)有：10 m高度風(fēng)速、水平風(fēng)速廓線、湍流強(qiáng)度、陣風(fēng)系數(shù)等。臺(tái)風(fēng)作為低壓渦旋性天氣系統(tǒng)，其自身的風(fēng)特性參數(shù)顯著區(qū)別于輸電線路結(jié)構(gòu)設(shè)計(jì)所依據(jù)的一般風(fēng)特性參數(shù)［7-8］，因此引起了工程技術(shù)人員的廣泛關(guān)注。目前在結(jié)構(gòu)風(fēng)工程領(lǐng)域中，關(guān)于臺(tái)風(fēng)影響的研究手段主要有現(xiàn)場(chǎng)觀測(cè)、風(fēng)洞實(shí)驗(yàn)、數(shù)值模擬等3 種。其中，現(xiàn)場(chǎng)觀測(cè)利用風(fēng)速傳感器、超聲測(cè)風(fēng)儀等觀測(cè)儀器獲取臺(tái)風(fēng)風(fēng)場(chǎng)內(nèi)實(shí)際風(fēng)特性參數(shù)，并配合加速度傳感器、應(yīng)變計(jì)等傳感器獲取結(jié)構(gòu)的真實(shí)動(dòng)力響應(yīng)進(jìn)行分析［9-10］。風(fēng)洞實(shí)驗(yàn)依托大氣邊界層風(fēng)洞設(shè)施，在實(shí)驗(yàn)室環(huán)境中還原臺(tái)風(fēng)近地面的平均風(fēng)速廓線，針對(duì)不同形式的線-塔氣彈模型，獲取動(dòng)力響應(yīng)特征［11］。數(shù)值模擬則借助計(jì)算流體力學(xué)（computational fluid dynamics，CFD）軟件模擬不同地形下的臺(tái)風(fēng)風(fēng)場(chǎng)空間分布，由于CFD模擬過(guò)程可重復(fù)，模擬結(jié)果具有物理意義，可以與現(xiàn)場(chǎng)觀測(cè)、風(fēng)洞試驗(yàn)等相互印證，互為補(bǔ)充，因此特別適合復(fù)雜地形的研究工作［12-16］。

文獻(xiàn)［17］在1604 號(hào)強(qiáng)臺(tái)風(fēng)“妮妲”登陸期間，在距離臺(tái)風(fēng)中心約15 km 的丘陵地形，利用激光測(cè)風(fēng)雷達(dá)進(jìn)行了輸電線路局地風(fēng)場(chǎng)觀測(cè)。針對(duì)臺(tái)風(fēng)“妮妲”，本文利用中、微尺度兩類(lèi)數(shù)值模擬工具，開(kāi)展臺(tái)風(fēng)期間的沿海輸電線路局地風(fēng)場(chǎng)模擬研究，驗(yàn)證中、微尺度嵌套模擬方法的適用性，分析臺(tái)風(fēng)影響下水平風(fēng)速廓線的差異及成因，并為今后沿海輸電線路臺(tái)風(fēng)防災(zāi)減災(zāi)研究提出建議。

1 模型與資料

1.1 中尺度氣象模式

選用由美國(guó)國(guó)家大氣研究中心（National Center for Atmospheric Research，NCAR）、國(guó)家環(huán)境預(yù)報(bào)中心（National Centers for Environmental Prediction，NCEP）以及多個(gè)大學(xué)、研究院所聯(lián)合研發(fā)的新一代區(qū)域中尺度模式（weather research and forecast，WRF），對(duì)1604 號(hào)強(qiáng)臺(tái)風(fēng)“妮妲”登陸期間的臺(tái)風(fēng)風(fēng)場(chǎng)開(kāi)展數(shù)值模擬。WRF 模式包含三維大氣運(yùn)動(dòng)方程、連續(xù)方程、能量守恒方程、大氣狀態(tài)方程、水汽守恒方程等組成的動(dòng)力框架，以及輻射傳輸、云微物理、行星邊界層、近地層與陸面過(guò)程等物理過(guò)程參數(shù)化方案，重點(diǎn)解決1～10 km 空間分辨率尺度、時(shí)效為60h 以?xún)?nèi)的有限區(qū)域天氣預(yù)報(bào)和模擬問(wèn)題，也常與CFD 軟件聯(lián)合使用［18-20］。本文所使用的WRF版本為v4.1.5。

1.2 微尺度CFD工具

選用挪威WindSim AS 公司開(kāi)發(fā)的CFD 計(jì)算流體力學(xué)軟件WindSim，對(duì)1604 號(hào)強(qiáng)臺(tái)風(fēng)“妮妲”登陸期間的微地形風(fēng)場(chǎng)進(jìn)行數(shù)值模擬。WindSim 軟件內(nèi)置英國(guó)CHAM 公司的Phoenics CFD 求解器，通過(guò)求解地形跟隨坐標(biāo)下的雷諾平均Navier-Stokes 方程（Reynolds-averaged Navier-Stokes equations，RANS），重點(diǎn)解決10～100 m 空間分辨率量級(jí)的各種復(fù)雜地形條件下近地層風(fēng)場(chǎng)模擬問(wèn)題。該軟件計(jì)算收斂速度快，能考慮不同大氣穩(wěn)定度，邊界條件設(shè)置靈活，在風(fēng)電領(lǐng)域應(yīng)用較廣泛［21-22］。本文所使用的WindSim版本為v10.0。

1.3 實(shí)測(cè)資料

1.3.1 地面觀測(cè)

“妮妲”登陸期間，香港橫瀾島氣象站、香港國(guó)際機(jī)場(chǎng)氣象站、澳門(mén)國(guó)際機(jī)場(chǎng)氣象站記錄了較為完整的風(fēng)速、風(fēng)向、氣壓等氣象要素，可用來(lái)對(duì)中尺度氣象模式的模擬結(jié)果進(jìn)行驗(yàn)證。

1.3.2 梯度觀測(cè)

除參考文獻(xiàn)［17］中的激光測(cè)風(fēng)雷達(dá)觀測(cè)資料，深圳氣象局350 m 氣象梯度觀測(cè)鐵塔也獲取了“妮妲”登陸期間的實(shí)測(cè)風(fēng)速、風(fēng)向特征［23］，上述風(fēng)廓線觀測(cè)資料可用來(lái)對(duì)微尺度CFD的模擬結(jié)果進(jìn)行驗(yàn)證。實(shí)測(cè)資料說(shuō)明見(jiàn)表1。

表1 用于驗(yàn)證的實(shí)測(cè)資料Tab.1 Observation data for simulation validation

2 分析方法

2.1 模擬方案

設(shè)置3 組數(shù)值試驗(yàn)，記為T(mén)1、T2 和T3。T1 為基于WRF 的中尺度模擬，T2 為基于WindSim 的微尺度模擬，T3 為WRF +WindSim 的中、微尺度嵌套模擬。

2.1.1 T1試驗(yàn)設(shè)置

在T1 試驗(yàn)中，共設(shè)4 層水平網(wǎng)格，前3 層模擬區(qū)域如圖1（a）所示，最內(nèi)層模擬區(qū)域如圖1（b）所示；空間分辨率依次為27 km、9 km、3 km、1 km，每層網(wǎng)格數(shù)分別為251×201、196×166、196×166、193×163。模式垂直分層為50 層，為了更精確地模擬近地面風(fēng)廓線形態(tài)，200 m 高度以下設(shè)置8 層，模式頂氣壓50 hPa。模擬時(shí)段為8 月1 日00：00 時(shí)～8 月2 日23：00 時(shí)，模式時(shí)間積分步長(zhǎng)135 s（確保積分穩(wěn)定），輸出逐10 min的風(fēng)速、風(fēng)向、氣壓、氣溫等要素（運(yùn)行前12 h 作為模式啟動(dòng)時(shí)間不予分析）。模式運(yùn)行的氣象初始場(chǎng)與邊界場(chǎng)均來(lái)自歐洲中期天氣預(yù)報(bào)中心（European Centre for Medium-Range Weather Forecasts，ECMWF）的ERA5 再分析資料。地形高程與下墊面類(lèi)型等靜態(tài)數(shù)據(jù)來(lái)自模式自帶的地形數(shù)據(jù)集（Global Multiresolution Terrain Elevation Data 2010，GMTED2010）與基于中分辨率成像光譜儀（moderate-resolution imaging spectroradiometer，MODIS）反演數(shù)據(jù)的土地利用類(lèi)型數(shù)據(jù)集，空間分辨率30″（約0.9km）。根據(jù)WRF 模式在以往臺(tái)風(fēng)個(gè)例模擬上的最佳物理過(guò)程參數(shù)化方案選項(xiàng)與組合研究成果［24-26］，并且在本次個(gè)例模擬上經(jīng)過(guò)反復(fù)比選，臺(tái)風(fēng)路徑與強(qiáng)度模擬誤差最小的主要物理過(guò)程參數(shù)化方案如下：長(zhǎng)、短波輻射方案為RRTMG（rapid radiative transfer model for global climate model），云微物理方案為WSM5（WRF-singlemoment-microphysics classes 5），行星邊界層參數(shù)化方案為YSU（Yonsei University），陸面過(guò)程方案為Noah（Noah land model）。

圖1 T1、T2、T3試驗(yàn)中的模擬區(qū)域示意圖Fig.1 Schematic diagrams of simulation area in T1,T2 and T3 case

2.1.2 T2試驗(yàn)設(shè)置

在T2 試驗(yàn)中，模擬區(qū)域如圖1（c）所示，即以激光雷達(dá)為中心的10 km×10 km 范圍；中心區(qū)域3 km×3 km 范圍進(jìn)行水平網(wǎng)格加密，加密區(qū)域內(nèi)空間分辨率相同（600 m、300 m、100 m、30 m 共4 種情形），其余區(qū)域的網(wǎng)格分辨率依次線性增加。垂直分層52 層，為了更精確地模擬近地面風(fēng)廓線形態(tài)，離地10～190 m 高度間隔20 m 設(shè)1 層；模擬區(qū)域頂高程5 000 m（確保網(wǎng)格空區(qū)比大于0.95，避免非物理加速影響）。對(duì)特定風(fēng)向下的定常風(fēng)場(chǎng)進(jìn)行CFD 模擬，地形高程與下墊面類(lèi)型數(shù)據(jù)來(lái)自ASTER GDEM 數(shù)據(jù)集與Globalcover 數(shù)據(jù)集。依照WindSim 軟件在復(fù)雜地形環(huán)境的應(yīng)用經(jīng)驗(yàn)，確保CFD 模擬結(jié)果穩(wěn)定收斂的模型設(shè)置有：求解器GCV（generalized cross validation），湍流模型RNG（renormalization group），迭代殘差收斂閾值0.001，大氣穩(wěn)定度為中性，頂邊界條件為固定壓力，側(cè)邊界條件為軟件默認(rèn)廓線形式。

2.1.3 T3試驗(yàn)設(shè)置

在T2 試驗(yàn)的基礎(chǔ)上，利用WindSim 軟件提供WRF 模式接口程序，將T1 試驗(yàn)中的逐10 min 經(jīng)向風(fēng)、緯向風(fēng)、垂直速度、氣壓和氣溫等輸出結(jié)果插值到圖1（c）的格點(diǎn)上，形成“妮妲”影響期間隨時(shí)間變化的風(fēng)、溫廓線，代替T2 試驗(yàn)中的邊界條件與大氣穩(wěn)定度設(shè)置，進(jìn)行逐10 min 風(fēng)場(chǎng)的CFD模擬。

2.2 模擬誤差統(tǒng)計(jì)

設(shè)置平均誤差E、絕對(duì)誤差M、均方根誤差S和相關(guān)系數(shù)R等4個(gè)統(tǒng)計(jì)指標(biāo)對(duì)T1試驗(yàn)的模擬時(shí)間序列進(jìn)行驗(yàn)證，表達(dá)式如下。

式中：xs，i為模擬序列中的第i個(gè)樣本；xo，i為實(shí)測(cè)序列中的第i個(gè)樣本；n為序列樣本數(shù)量；為模擬序列均值；為實(shí)測(cè)序列均值。

2.3 模擬結(jié)果分析

設(shè)置無(wú)量綱水平風(fēng)速廓線V'（z）與無(wú)量綱垂直速度廓線W'（z），對(duì)各試驗(yàn)中的風(fēng)廓線形態(tài)進(jìn)行驗(yàn)證，表達(dá)式如下。

式中：V（z）為任一高度z的水平風(fēng)速；V100m為離地100 m 高度的水平風(fēng)速；W（z）為任一高度z的垂直速度；W100m為離地100 m高度的垂直速度。

另外，設(shè)置徑向速度Vr，對(duì)T3 試驗(yàn)中的風(fēng)場(chǎng)空間分布形態(tài)進(jìn)行對(duì)比，表達(dá)式如下。

式中：u、v、w分別為WindSim 輸出的x-y-z坐標(biāo)系中的風(fēng)速分量，分別代表東向、北向和垂直方向上的風(fēng)速值；φ為方位角（與正東方向的夾角）；θ為傾角（與天頂之間的夾角）。

3 結(jié)果與討論

3.1 T1試驗(yàn)

如圖2所示，就激光雷達(dá)點(diǎn)位而言，基于WRF的中尺度模擬較好地再現(xiàn)了“妮妲”登陸期間的氣象要素時(shí)間變化過(guò)程，氣壓的“U”型變化、風(fēng)速的“雙峰”型變化、超過(guò)120°的風(fēng)向連續(xù)變化特征均在模擬值中得以體現(xiàn)。此外，模擬值在氣象梯度鐵塔、香港橫瀾島氣象站、香港國(guó)際機(jī)場(chǎng)氣象站、澳門(mén)國(guó)際機(jī)場(chǎng)氣象站等點(diǎn)位也呈現(xiàn)出與實(shí)測(cè)值較為一致的變化過(guò)程。各點(diǎn)位的模擬誤差統(tǒng)計(jì)結(jié)果如表2所示，氣壓的相關(guān)系數(shù)在0.965～0.986 之間，風(fēng)速的相關(guān)系數(shù)在0.682～0.941 之間，風(fēng)向的相關(guān)系數(shù)在0.766～0.988 之間，印證了中尺度模式模擬的可行性。

圖2 T1試驗(yàn)中的激光雷達(dá)點(diǎn)位結(jié)果對(duì)比Fig.2 Comparison of LiDAR locations in T1 case

表2 T1試驗(yàn)的模擬值誤差統(tǒng)計(jì)Tab.2 Simulation error statistics in T1 case

另外，就氣壓而言，各點(diǎn)位的平均誤差、平均絕對(duì)誤差、均方根誤差均很小，數(shù)值不超過(guò)2 hPa，反映了WRF 模式對(duì)氣壓場(chǎng)強(qiáng)度的模擬具有高精度。就風(fēng)速而言，大部分點(diǎn)位的平均誤差為正值，表明WRF 模式對(duì)風(fēng)速有一定程度的高估。上述情況與模式空間分辨率有關(guān)，即使T1 試驗(yàn)中的最高空間分辨率為1 km，但由于WRF模式自身的限制，10～100 m 空間分辨率量級(jí)的微小地形仍然無(wú)法在模式中得到體現(xiàn)，從而造成模式中地形對(duì)風(fēng)場(chǎng)的阻塞效應(yīng)可能被低估，相應(yīng)的絕對(duì)誤差、均方根誤差基本穩(wěn)定在2.2～3.6 m/s之間，模擬精度不及氣壓。就風(fēng)向而言，各點(diǎn)位的平均誤差在14 °～30 °之間，若按16 個(gè)風(fēng)向扇區(qū)而言，誤差控制在左、右1 個(gè)扇區(qū)之內(nèi)，但部分點(diǎn)位的絕對(duì)誤差、均方根誤差較大，其原因也可能與中尺度模式中的地形效應(yīng)偏差，以及本次臺(tái)風(fēng)眼區(qū)影響時(shí)段內(nèi)觀測(cè)風(fēng)向較為離散有關(guān)。

進(jìn)一步利用無(wú)量綱水平風(fēng)速廓線V'（z），對(duì)激光雷達(dá)和氣象鐵塔點(diǎn)位的局地風(fēng)廓線進(jìn)行驗(yàn)證，如圖3所示。依照參考文獻(xiàn)［17］的做法，選取3個(gè)時(shí)段進(jìn)行對(duì)比：時(shí)段A（8 月1 日23：00 時(shí)—8 月2 日02：00 時(shí)，對(duì)應(yīng)“妮妲”登陸前的眼壁強(qiáng)風(fēng)區(qū)影響時(shí)段）、時(shí)段B（8月2日04：00時(shí)—06：00時(shí)，對(duì)應(yīng)“妮妲”眼區(qū)影響時(shí)段）、時(shí)段C（8 月2 日07：00 時(shí)—10：00時(shí)，對(duì)應(yīng)“妮妲”登陸后的眼壁強(qiáng)風(fēng)區(qū)影響時(shí)段）。由圖3可知，時(shí)段B的模擬廓線與實(shí)測(cè)值較為一致，而時(shí)段A、C 的模擬廓線在200 m 高度以上出現(xiàn)明顯偏離。具體表現(xiàn)為實(shí)測(cè)V'（z）更加傾斜，而模擬值的高、低層之間水平風(fēng)速梯度小于實(shí)測(cè)值，原因也與局地風(fēng)場(chǎng)的地形阻塞效應(yīng)可能在模式中被低估有關(guān)，并且這種情況在大風(fēng)影響時(shí)段內(nèi)更為顯著。

圖3 T1試驗(yàn)中的水平風(fēng)速廓線對(duì)比Fig.3 Comparison of horizontal wind speed profiles in T1 case

3.2 T2試驗(yàn)

鑒于中尺度氣象模式在空間分辨率上的限制，通過(guò)基于WindSim 的微尺度模擬得到更高分辨率的風(fēng)場(chǎng)模擬結(jié)果，并對(duì)比不同分辨率情形下的無(wú)量綱水平風(fēng)速廓線V'（z）與無(wú)量綱垂直速度廓線W'（z），揭示地形影響局地風(fēng)場(chǎng)的細(xì)節(jié)。

如圖4 所示，對(duì)于激光雷達(dá)點(diǎn)位，利用335 °～15 °之間的7 個(gè)典型風(fēng)向，模擬時(shí)段A 內(nèi)臺(tái)風(fēng)眼壁強(qiáng)風(fēng)的可能風(fēng)向，模擬結(jié)果均穩(wěn)定收斂（迭代殘差小于或等于閾值），結(jié)果可信。隨著空間分辨率提高，各個(gè)典型風(fēng)向下無(wú)量綱水平風(fēng)速廓線V'（z）與無(wú)量綱垂直速度廓線W'（z）出現(xiàn)顯著差異，反映了不同空間尺度的地形效應(yīng)。另一方面，當(dāng)分辨率提高為100 m 甚至更高的30 m 時(shí)，會(huì)出現(xiàn)某個(gè)風(fēng)向或某兩個(gè)臨近風(fēng)向的廓線模擬值迅速逼近實(shí)測(cè)值的現(xiàn)象，時(shí)段C 則利用150 °～190 °之間的8 個(gè)典型風(fēng)向進(jìn)行模擬，也出現(xiàn)類(lèi)似趨勢(shì)，上述現(xiàn)象與臺(tái)風(fēng)環(huán)流背景下局地風(fēng)廓線的形成機(jī)理密切相關(guān)。

圖4 T2試驗(yàn)中激光雷達(dá)點(diǎn)位的水平風(fēng)速廓線對(duì)比Fig.4 Comparison of horizontal wind speed profiles at LiDAR location in T2 case

進(jìn)一步利用WindSim 的定向模擬結(jié)果來(lái)分析激光雷達(dá)點(diǎn)位的風(fēng)廓線形成機(jī)理。根據(jù)參考文獻(xiàn)［17］的介紹，激光雷達(dá)觀測(cè)點(diǎn)位周?chē)匦纬尸F(xiàn)北高南低的特征，其西北方向和東南方向有低矮丘陵分布。如圖5所示，選擇激光雷達(dá)點(diǎn)位南-北方向的地形剖面進(jìn)行分析，在時(shí)段A 內(nèi)，受西北偏北氣流控制，當(dāng)氣流經(jīng)過(guò)雷達(dá)北側(cè)2 km 附近的丘陵地形時(shí)，由于受到山體阻擋、遲滯，翻山后氣流下沉，風(fēng)速降低，形成回流區(qū)并擴(kuò)展至激光雷達(dá)點(diǎn)位上空，導(dǎo)致激光雷達(dá)V'（z）更為傾斜，W'（z）出現(xiàn)負(fù)值。在時(shí)段C 內(nèi)，受偏南氣流控制，氣流受地形強(qiáng)迫抬升，導(dǎo)致激光雷達(dá)W'（z）出現(xiàn)正值?？梢钥吹?，上述現(xiàn)象與特定風(fēng)向內(nèi)10～100 m 空間分辨率量級(jí)的局地地形效應(yīng)有關(guān)，在T2 試驗(yàn)中隨著空間分辨率的提升，這種局地地形效應(yīng)才得以體現(xiàn)，而在空間分辨率有限的T1試驗(yàn)中，無(wú)法被中尺度氣象模式所反映。

圖5 T2試驗(yàn)中的局地風(fēng)速剖面對(duì)比（南-北方向）Fig.5 Comparison of local wind speed profiles h in T2 case（from north to sout）

值得注意的是，由于臺(tái)風(fēng)眼壁強(qiáng)風(fēng)區(qū)內(nèi)湍流活動(dòng)旺盛，時(shí)段A、C內(nèi)的大氣穩(wěn)定度可能偏離中性，因此基于中性大氣穩(wěn)定度的T2 試驗(yàn)在垂直速度模擬上仍存在偏差。另一方面，時(shí)段A、C 內(nèi)風(fēng)向會(huì)發(fā)生一定幅度擺動(dòng)；而隨著風(fēng)向的不同，氣流途徑地形的坡度、坡向發(fā)生變化，地形效應(yīng)也出現(xiàn)差異。因此基于固定邊界條件的T2 試驗(yàn)無(wú)法反映臺(tái)風(fēng)環(huán)流的實(shí)際特征，存在一定不足。

3.3 T3試驗(yàn)

鑒于T2 試驗(yàn)的邊界條件與大氣穩(wěn)定度局限，通過(guò)WRF +WindSim 的中、微尺度嵌套模擬，響應(yīng)臺(tái)風(fēng)眼壁強(qiáng)風(fēng)區(qū)逐10 min的風(fēng)速、風(fēng)向、氣壓、氣溫變化，并將高空間分辨率的CFD模擬擴(kuò)展至整個(gè)時(shí)段A和時(shí)段C。

同樣針對(duì)激光雷達(dá)點(diǎn)位進(jìn)行對(duì)比，將時(shí)段A、C 全時(shí)段內(nèi)的水平風(fēng)速廓線、垂直速度廓線進(jìn)行平均，如圖6 所示。WRF +WindSim 的中、微尺度嵌套模擬既修正了單獨(dú)WRF 模擬的地形效應(yīng)偏差，又改善了單獨(dú)WindSim 模擬的邊界條件和大氣穩(wěn)定度，垂直速度模擬精度更高?？梢钥吹?，中、微尺度嵌套模擬方式結(jié)合了中尺度氣象模式、微尺度CFD 工具各自的模擬優(yōu)勢(shì)，對(duì)比圖3—5 的模擬值，與實(shí)測(cè)值的偏差更小。

圖6 T3試驗(yàn)中激光雷達(dá)點(diǎn)位的水平風(fēng)速廓線對(duì)比Fig.6 Comparison of horizontal wind speed profiles at LiDAR location in T3 case

筆者在時(shí)段C 的激光雷達(dá)觀測(cè)過(guò)程中，還采集了激光雷達(dá)水平方位掃描（PPI）的徑向速度數(shù)據(jù)，在此按照2.3節(jié)的表達(dá)式，將T3試驗(yàn)中時(shí)段C 的三維風(fēng)場(chǎng)模擬結(jié)果合成出徑向速度Vr，與激光雷達(dá)實(shí)際測(cè)量的Vr進(jìn)行對(duì)比。如圖7（a）所示，WRF+WindSim 的中、微尺度嵌套模擬得到的徑向速度場(chǎng)中的零速度線方位、正值區(qū)和負(fù)值區(qū)的方位與覆蓋范圍、徑向速度的具體大小等均與圖7（b）的實(shí)測(cè)結(jié)果具有一致性。由于徑向速度是風(fēng)矢量在某一特定方位和傾角下的投影，綜合體現(xiàn)了局地風(fēng)速、風(fēng)向的實(shí)際特征，因而徑向速度模擬的精確程度在很大程度上可以反映三維風(fēng)場(chǎng)模擬的準(zhǔn)確性。因此認(rèn)為WRF +WindSim 的中、微尺度嵌套模擬具備臺(tái)風(fēng)影響期間復(fù)雜地形下的三維風(fēng)場(chǎng)模擬能力。

圖7 T3試驗(yàn)中風(fēng)場(chǎng)空間分布對(duì)比Fig.7 Comparison of wind field spatial distributions in T3 case

3.4 討論

結(jié)合3.1至3.3節(jié)的結(jié)果分析，T3試驗(yàn)中的中、微尺度嵌套模擬方式綜合了中尺度氣象模式與微尺度CFD工具各自的優(yōu)勢(shì)，在臺(tái)風(fēng)“妮妲”登陸期間沿海局地風(fēng)場(chǎng)的數(shù)值模擬上具有更高精度，能夠?yàn)轱L(fēng)災(zāi)事故的原因分析和情景再現(xiàn)提供更為豐富的風(fēng)場(chǎng)資料；下一步還需要收集更多臺(tái)風(fēng)個(gè)例進(jìn)行模擬結(jié)果檢驗(yàn)，確定中、微尺度嵌套模擬方式的可靠性。

另外，值得注意的是，T3 試驗(yàn)中的中、微尺度嵌套模擬方式在提高模擬精度的同時(shí)，逐10 min的CFD 模擬也需要消耗較大量的計(jì)算資源與時(shí)間，從而無(wú)法直接在沿海輸電線路風(fēng)災(zāi)預(yù)警中應(yīng)用。因此在下一步的研究中還需要發(fā)展新的中、微尺度嵌套預(yù)測(cè)方法，在提高空間分辨率的同時(shí)保證相應(yīng)的預(yù)報(bào)時(shí)效。

4 結(jié)論

本文針對(duì)1614 號(hào)強(qiáng)臺(tái)風(fēng)“妮妲”登陸期間的沿海輸電線路局地風(fēng)場(chǎng)進(jìn)行數(shù)值模擬研究，通過(guò)單獨(dú)中尺度模擬、單獨(dú)微尺度模擬以及中、微尺度嵌套模擬等三組試驗(yàn)，并與臺(tái)風(fēng)期間的實(shí)測(cè)資料進(jìn)行對(duì)比，取得如下結(jié)論。

1）基于WRF 的中尺度模擬較好地再現(xiàn)出“妮妲”登陸期間的氣象要素時(shí)間變化過(guò)程，WRF 模式對(duì)氣壓場(chǎng)強(qiáng)度的模擬具有高精度，對(duì)風(fēng)場(chǎng)的模擬精度不及氣壓，眼壁強(qiáng)風(fēng)區(qū)影響時(shí)段內(nèi)的風(fēng)廓線與實(shí)測(cè)值相比出現(xiàn)偏離，原因與局地風(fēng)場(chǎng)的地形阻塞效應(yīng)在WRF模式中被低估有關(guān)。

2）基于WindSim 的微尺度模擬較好地再現(xiàn)特定風(fēng)向內(nèi)10～100 m 空間分辨率量級(jí)的風(fēng)場(chǎng)地形效應(yīng)；當(dāng)空間分辨率提高為100 m 甚至更高的30 m 時(shí)，會(huì)出現(xiàn)某個(gè)風(fēng)向或某兩個(gè)臨近風(fēng)向的風(fēng)廓線迅速逼近實(shí)測(cè)值的現(xiàn)象，揭示了臺(tái)風(fēng)環(huán)流背景下局地風(fēng)廓線的形成機(jī)理。

3）WRF +WindSim 的中、微尺度嵌套模擬既能修正單獨(dú)WRF 模擬的地形效應(yīng)偏差，又能改善單獨(dú)WindSim 模擬的邊界條件和大氣穩(wěn)定度，具備臺(tái)風(fēng)影響期間復(fù)雜地形下的三維風(fēng)場(chǎng)模擬能力，模擬結(jié)果精度高于單獨(dú)中、微尺度模擬方式。