劉 非,李陪然,陳 垚,袁紹春,劉 臻,曾梓健,任春嬌,劉 俊,張 茅,朱俊成
(1. 重慶交通大學(xué)河海學(xué)院,重慶 400074; 2. 重慶交通大學(xué)環(huán)境水利工程重慶市工程實(shí)驗(yàn)室,重慶 400074;3. 中建二局重慶建設(shè)發(fā)展有限公司,重慶 401329; 4. 中國建筑第二工程局有限公司,北京 100160)
氣候變暖和人類活動(dòng)改變了自然界的水文循環(huán),增加了極端水文事件發(fā)生的概率,突發(fā)性暴雨引發(fā)的城市內(nèi)澇問題日益嚴(yán)重,給我國人民的生命財(cái)產(chǎn)安全帶來嚴(yán)重威脅。
隨著計(jì)算機(jī)技術(shù)和模型理論的不斷發(fā)展,數(shù)值模型在城市內(nèi)澇治理和人員防災(zāi)中扮演越來越重要的角色?;谖锢磉^程的雨洪模型可以較好再現(xiàn)和預(yù)測(cè)研究區(qū)內(nèi)澇發(fā)生規(guī)律,提供較完整的積水范圍、水深、浸水時(shí)間等洪水時(shí)空要素信息,對(duì)防洪管理和應(yīng)急減災(zāi)意義重大。雨洪模型的準(zhǔn)確性和運(yùn)行效率通常受到諸多因素影響,比如輸入數(shù)據(jù)的完整性和精確性、計(jì)算區(qū)域大小和排水系統(tǒng)規(guī)模,以及計(jì)算設(shè)備和模型算法的先進(jìn)性等。目前,國內(nèi)外使用較廣泛的一維管網(wǎng)模型有美國環(huán)保局開發(fā)的暴雨徑流管理模型(Storm Water Management Model,SWMM)[1],二維水文水動(dòng)力模型軟件有InfoWorks ICM、MIKE Urban 等。在眾多開源或半開源二維水動(dòng)力模型中,LISFLOOD-FP 模型具有學(xué)習(xí)成本低、計(jì)算效率高的優(yōu)點(diǎn),但在國內(nèi)推廣使用不多[2,3]。無論采用哪一款模型,管網(wǎng)資料通常都是進(jìn)行城市內(nèi)澇模擬分析的重要基礎(chǔ)數(shù)據(jù),因其特殊性和涉密性,往往不易獲得。因而,無管勘資料背景下的城市內(nèi)澇模擬技術(shù)是當(dāng)下雨洪管理領(lǐng)域的一個(gè)熱點(diǎn)和難點(diǎn)。
近年來,已有部分學(xué)者對(duì)此進(jìn)行了探索和嘗試,并總結(jié)出了目前較為常用的管網(wǎng)替代方法,針對(duì)研究區(qū)完全沒有管網(wǎng)數(shù)據(jù)的情形,主要以下滲等效法和降雨率折現(xiàn)法[4,5]為主,前者包括雨水井等效排水法、道路等效排水法、全域等效排水法等;對(duì)于研究區(qū)域尚有部分管網(wǎng)資料,但不完整的情況,除了下滲等效法外,還有管網(wǎng)概化法等。楊東等[6]以陜西省西咸新區(qū)為研究區(qū),對(duì)比了不同下滲等效法和降雨率折現(xiàn)法,并指出雨水井等效排水法的實(shí)際效果最佳。杜佳鍇等[7]以粵港澳大灣區(qū)某流域?yàn)檠芯繀^(qū),探討了降雨率折現(xiàn)法、下滲增加法和管網(wǎng)概化法的適用性。以上方法中,下滲等效法為代表的管網(wǎng)替代方法對(duì)數(shù)據(jù)依賴程度最低,且可通過調(diào)整下滲量對(duì)模型準(zhǔn)確性進(jìn)行調(diào)整,在無管網(wǎng)資料地區(qū)具有良好適用性;降雨率折現(xiàn)法由于是將降雨量降低達(dá)到替代管網(wǎng)排水作用,極易導(dǎo)致模擬的內(nèi)澇程度較實(shí)際偏低。在下滲等效法中,雨水井等效法在模擬積水面積準(zhǔn)確度上最優(yōu)[6]。但既有的研究較少有提供高效準(zhǔn)確的雨水檢查井?dāng)?shù)據(jù)獲取方法,并將其應(yīng)用到開源或半開源模型中(例如LISFLOOD-FP 模型),形成從基礎(chǔ)數(shù)據(jù)獲取到構(gòu)建模型的完整技術(shù)方法流程。若采取人工實(shí)地調(diào)查獲取檢查井位置的方式,當(dāng)研究區(qū)域范圍較大時(shí),會(huì)有成本過高和效率較低的缺點(diǎn)。因此,亟需一種基于較便捷的檢查井?dāng)?shù)據(jù)獲取方式和開源(或半開源)模型的無管勘資料地區(qū)雨洪模擬技術(shù)。
為解決城市內(nèi)澇模擬研究中管網(wǎng)資料缺失的問題,本文以山地城市重慶某排水分區(qū)為研究區(qū)域,利用百度地圖開放平臺(tái)提供的全景靜態(tài)API 功能,對(duì)街道檢查井進(jìn)行識(shí)別和定位。在此基礎(chǔ)上,發(fā)揮LISFLOOD-FP 模型計(jì)算高效的優(yōu)勢(shì),采用雨水井等效排水法進(jìn)行雨洪建模,對(duì)不同暴雨情景下的研究區(qū)內(nèi)澇積水面積、積水量、積水分布等進(jìn)行分析。模型采用直接降雨法進(jìn)行全水動(dòng)力地表徑流計(jì)算,并與SWMM 耦合考察橋面徑流對(duì)局部?jī)?nèi)澇的影響。所提出的模擬計(jì)算方法可為類似無管勘資料山地城市區(qū)域利用LISFLOOD-FP 模型進(jìn)行暴雨內(nèi)澇風(fēng)險(xiǎn)高效評(píng)估提供參考。
研究區(qū)域位于重慶市主城南岸區(qū)某排水分區(qū),面積約4.98 km2,如圖1(a)所示,北面通過菜園壩大橋和重慶長江大橋與渝中區(qū)相連,在南北交通上具有十分突出的地位。根據(jù)《南岸區(qū)城市內(nèi)澇治理系統(tǒng)化實(shí)施方案》和《中心城區(qū)排水防澇專項(xiàng)規(guī)劃(修編)》,該片區(qū)存在多處易澇點(diǎn),對(duì)城市交通運(yùn)行構(gòu)成不利影響。由于缺少高精度的地形資料,研究區(qū)地形數(shù)據(jù)從地理空間數(shù)據(jù)云平臺(tái)免費(fèi)獲?。╤ttps://www.gscloud.cn/home),空間分辨率為30 m。
圖1 研究區(qū)域范圍和地形Fig.1 The study area and topography
采用國內(nèi)雨洪模擬研究中適用性良好的芝加哥降雨過程線[8]進(jìn)行暴雨情景設(shè)計(jì)。考慮到國內(nèi)洪澇規(guī)劃標(biāo)準(zhǔn)中對(duì)于重現(xiàn)期的取值一般不超過百年一遇,而諸如2021年“7·21”鄭州極端暴雨事件發(fā)生的幾率在氣候變化背景下呈逐漸上升趨勢(shì),因此,本研究中增加了500年一遇降雨情景作為極端情況,設(shè)計(jì)降雨重現(xiàn)期包括1、5、20、50、100、500 年等6 種情況,雨峰系數(shù)設(shè)為0.17[9],降雨歷時(shí)為120 min。參考《南岸區(qū)城市內(nèi)澇治理系統(tǒng)化實(shí)施方案》,重慶市南岸區(qū)暴雨強(qiáng)度公式為式(1),設(shè)計(jì)降雨過程線如圖2所示。
圖2 設(shè)計(jì)降雨過程線Fig.2 Designed rainfall hyetograph
圖3 研究區(qū)(藍(lán)色方框)雨水檢查井分布Fig.3 The distribution of extracted manholes in the study area
圖4 基于SWMM的橋面徑流模型Fig.4 Bridge deck runoff model based on SWMM
圖5 模擬計(jì)算框架Fig.5 The modeling framework
式中:q為設(shè)計(jì)暴雨強(qiáng)度,L/(s·hm2);P為設(shè)計(jì)暴雨重現(xiàn)期,a;t為降雨歷時(shí),min。
利用百度地圖開放平臺(tái)提供的全景靜態(tài)API 功能,開發(fā)基于Python 的道路全景批量自動(dòng)下載和處理程序,對(duì)研究區(qū)內(nèi)的雨水檢查井進(jìn)行識(shí)別和定位。主要步驟包括:首先,在百度地圖開放平臺(tái)提供的坐標(biāo)拾取界面(https://api.map.baidu.com/lbsapi/getpoint/index.html)對(duì)研究區(qū)所有道路坐標(biāo)點(diǎn)進(jìn)行提?。蝗缓?,將提取的所有坐標(biāo)作為輸入,運(yùn)行所編寫的基于全景靜態(tài)API 接口的Python 程序,批量下載指定坐標(biāo)位置的街道全景照片;最后,通過對(duì)所下載的照片進(jìn)行篩選,從而確定檢查井的分布。
LISFLOOD-FP 模型[10-13]是由英國布里斯托大學(xué)開發(fā)的二維洪水動(dòng)力模型,以柵格化網(wǎng)格為計(jì)算單元,可實(shí)現(xiàn)基于淺水方程簡(jiǎn)化形式的一維河道水力和二維水動(dòng)力過程模擬。LISFLOOD-FP 模型針對(duì)不同應(yīng)用場(chǎng)景提供了多種模型求解器以及邊界條件輸入方式。構(gòu)建城市區(qū)域洪水漫流模型,所需準(zhǔn)備的數(shù)據(jù)文件主要包括地形(*.asc)、邊界條件(*.bci和*.bdy)和模型運(yùn)行文件(*.par)。在本研究中,對(duì)于二維地表漫流計(jì)算,采用基于完整淺水方程求解的FV1求解器。
在模型中采用雨水井等效排水法[6,14],將管網(wǎng)排水概化為通過雨水檢查井所在點(diǎn)排出,根據(jù)地區(qū)管網(wǎng)設(shè)計(jì)標(biāo)準(zhǔn)結(jié)合堰流公式對(duì)檢查井入流量進(jìn)行確定,而后將其作為點(diǎn)源邊界條件在模型中進(jìn)行設(shè)置[15]。
隨著城市立體交通的日趨成熟,橋梁作為其重要組成部分,往往建設(shè)在區(qū)域連接的關(guān)鍵節(jié)點(diǎn),一旦因降雨集中或橋梁排水系統(tǒng)運(yùn)行不暢導(dǎo)致橋面積水,不僅會(huì)直接影響橋面行車安全,甚至?xí)渣c(diǎn)帶面加重局部?jī)?nèi)澇程度。而在以往的城市內(nèi)澇研究中,對(duì)橋梁排水系統(tǒng)和橋面徑流關(guān)注較少,例如,劉非等[16]利用InfoWorks ICM 對(duì)武漢長江隧道內(nèi)暴雨積水過程進(jìn)行了模擬分析,但未考慮隧道口附近沙湖大橋橋面雨水徑流的影響。
為考察研究區(qū)北面兩座橋梁橋面徑流對(duì)局部?jī)?nèi)澇的影響,利用SWMM 模型進(jìn)行橋面徑流計(jì)算,橋面徑流量模擬結(jié)果作為邊界條件輸入LISFLOOD-FP 模型中進(jìn)行二維洪水地表漫流計(jì)算。
菜園壩大橋和長江大橋是連接渝中區(qū)與南岸區(qū)的過江通道,位于長江水道之上。菜園壩大橋線路全長7 km,主橋長約800 m,橋面為雙向六車道快速路(寬度約30.5 m)。長江大橋全長1.1 km,橋面為雙向八車道(寬度約40 m)。橋面橫向泄水管按照20 m間距布設(shè)考慮,橋面徑流系數(shù)取為0.95[17,18]。
暴雨作為內(nèi)澇產(chǎn)生的根本原因,其對(duì)內(nèi)澇積水的影響因素除雨型、降雨強(qiáng)度、降雨歷時(shí)外,還受到地形及用地類型等的作用。在傳統(tǒng)的半分布耦合模式[19]中,降雨作為輸入條件被關(guān)聯(lián)到每個(gè)匯水區(qū),分別計(jì)算產(chǎn)流,而后進(jìn)入到管網(wǎng)系統(tǒng),若管網(wǎng)排水不暢產(chǎn)生節(jié)點(diǎn)溢流,則通過檢查井與地表進(jìn)行水量交換。這種概化方式雖有一定水文學(xué)理論依據(jù),但忽略了徑流在進(jìn)入管網(wǎng)系統(tǒng)前的地表漫流過程,與實(shí)際匯流仍有差別。在Info-Works ICM 等軟件中,模型提供了與實(shí)際過程更符合的直接降雨法,但該方法缺點(diǎn)是計(jì)算較為耗時(shí),需采用較為高效的模型算法[20],但總體而言相關(guān)研究較少。
研究區(qū)采用直接降雨法進(jìn)行徑流模擬,即基于柵格的徑流模擬方法,將整個(gè)研究區(qū)按地形柵格劃分結(jié)構(gòu)計(jì)算單元,每個(gè)計(jì)算單元包含高程、粗糙系數(shù)、滲透參數(shù)、不透水區(qū)百分比和洼地深度等信息,地表產(chǎn)流由每個(gè)單元獨(dú)立產(chǎn)生[21]。每個(gè)計(jì)算單元按用地類型或地表覆蓋物種類劃分為透水區(qū)和不透水區(qū),分別進(jìn)行產(chǎn)流計(jì)算。對(duì)不透水區(qū),徑流量等于降雨強(qiáng)度,而在透水區(qū),徑流量等于降雨強(qiáng)度減去入滲損失。最終,所有計(jì)算單元所得產(chǎn)流均作為邊界條件輸入LISFLOOD-FP 模型進(jìn)行二維洪水地表漫流計(jì)算。
式中:Qm為第m個(gè)計(jì)算單元的產(chǎn)流量,m3/s;i為降雨強(qiáng)度,mm/min;f為下滲量,mm/h;Sm為第m個(gè)計(jì)算單元的面積,m2。
采用Horton模型對(duì)透水區(qū)入滲率進(jìn)行估算?;痉匠虨椋?/p>
式中:ft為t時(shí)刻的下滲率,mm/h;fc為穩(wěn)定下滲率,mm/h;f0為土壤初始下滲率,mm/h;k為衰減系數(shù),與土壤的物理性質(zhì)有關(guān),1/h;t為時(shí)間,h。參數(shù)取值參考SWMM模型手冊(cè)。
為考察將直接降雨法所得徑流量結(jié)果作為邊界條件驅(qū)動(dòng)LISFLOOD-FP 模型進(jìn)行洪水計(jì)算的可行性和準(zhǔn)確性,通過對(duì)理想算例進(jìn)行水量平衡分析作為驗(yàn)證。算例為長寬均為100 m的方形區(qū)域,以10 m 網(wǎng)格進(jìn)行劃分。建立同區(qū)域的SWMM 模型,同樣降雨作為輸入,對(duì)比驗(yàn)證公式法計(jì)算的徑流量精確值、SWMM模型模擬的徑流量、LISFLOOD-FP模型模擬所得洪水量三者的吻合程度。模型時(shí)間為4 h。降雨輸入和結(jié)果對(duì)比如圖6 所示。圖6(b)中LISFLOOD-FP 模型模擬的地表洪水量與公式法計(jì)算所得徑流量精確值、SWMM 模型模擬所得徑流量均吻合良好。計(jì)算區(qū)域初始和最終時(shí)刻的水位分布如圖7所示。
圖6 理想算例降雨輸入和不同徑流量模擬方式所得結(jié)果對(duì)比Fig.6 The rainfall input and comparison of different runoff calculation methods
圖7 初始和最終時(shí)刻區(qū)域水位分布Fig.7 Water surface elevation at the initial and final time step using LISFLOOD-FP model
利用2022 年7 月30 日實(shí)測(cè)場(chǎng)次降雨對(duì)所提出的模擬方法進(jìn)行驗(yàn)證,降雨數(shù)據(jù)來源為重慶交通大學(xué)山地城市降雨和內(nèi)澇風(fēng)險(xiǎn)監(jiān)測(cè)系統(tǒng),降雨測(cè)量裝置安裝在校園內(nèi)港航樓樓頂[如圖8(a)]。圖9 為模擬所得內(nèi)澇點(diǎn)與文獻(xiàn)[22,23]中所提到研究區(qū)典型易澇點(diǎn)對(duì)比,與實(shí)際監(jiān)測(cè)記錄位置基本一致。表1 和表2為該場(chǎng)降雨期間易澇點(diǎn)位置的峰值積水面積和積水深度模擬值與實(shí)測(cè)值對(duì)比,可見所提出的模擬方法在確定積水面積上準(zhǔn)確度較高,而對(duì)積水深度的高精度模擬存在不足,以及易澇點(diǎn)D的積水范圍和積水深度的模擬值與實(shí)測(cè)值相對(duì)誤差較大,可能是由于地形精度較低導(dǎo)致。
表1 峰值積水范圍模擬值與實(shí)測(cè)值對(duì)比Tab.1 Comparison of the peak flood area between the simulated and measured
表2 峰值積水深度模擬值與實(shí)測(cè)值對(duì)比Tab.2 Comparison of the peak flood depth between the simulated and measured
圖8 雨量計(jì)與實(shí)測(cè)降雨Fig.8 Rain gauge and the measured rain process
圖9 模擬所得內(nèi)澇點(diǎn)與實(shí)際易澇點(diǎn)對(duì)比(易澇點(diǎn)數(shù)據(jù)來源于文獻(xiàn)[22,23])Fig.9 Comparison of simulated flood-prone points and the reported ones
將不同重現(xiàn)期的降雨過程線應(yīng)用直接降雨法所得徑流量作為邊界條件輸入LISFLOOD-FP 模型,得到研究區(qū)內(nèi)澇積水變化過程。積水量可以表征內(nèi)澇區(qū)域的總體積水情況,對(duì)各重現(xiàn)期降雨條件下研究區(qū)的積水量和積水面積隨時(shí)間變化進(jìn)行統(tǒng)計(jì),如圖10 所示。總體而言,積水面積變化規(guī)律與對(duì)應(yīng)的降雨過程線變化較為相似,重現(xiàn)期越大,則同時(shí)刻的積水面積越大。積水面積峰值出現(xiàn)的時(shí)間與降雨重現(xiàn)期大小無關(guān),但相比雨峰有滯后;積水面積在降雨結(jié)束(7 200 s)后出現(xiàn)較快下降,并最終保持穩(wěn)定。隨著重現(xiàn)期的增加,積水量也在相應(yīng)的增加。通過對(duì)比積水面積變化,區(qū)域積水量增長速率最快的時(shí)刻與積水面積達(dá)到峰值的時(shí)刻基本一致,均在約1 800 s 時(shí),晚于雨峰出現(xiàn)時(shí)間。在積水面積達(dá)到峰值并開始逐漸下降階段,積水量仍在增加,這表明積水范圍雖在快速縮小,但陡坡路段可能面臨較大洪水流速,對(duì)行人,尤其是兒童和老年人仍存在較大威脅,需注意避開水流較急的危險(xiǎn)路段。
圖10 積水面積和積水量隨時(shí)間變化Fig.10 Time-series variation of the flood area and flood volume
內(nèi)澇積水分布圖可用于排查易澇點(diǎn)和分析積水原因,繪制出研究區(qū)不同降雨重現(xiàn)期下最大積水時(shí)刻的積水分布圖,如圖11所示。隨著降雨重現(xiàn)期的增加,積水點(diǎn)的數(shù)量和范圍均在增加,但易澇點(diǎn)(降雨重現(xiàn)期較小時(shí)即發(fā)生較嚴(yán)重積水的地段)的空間位置變化不大,在研究區(qū)最北端、中部、東南部、西南部均存在,只是當(dāng)重現(xiàn)期增大后,原先零散的易澇點(diǎn)有互相連接成片的趨勢(shì)。通過查看該區(qū)域地形[圖1(b)],可知內(nèi)澇點(diǎn)基本分布在地形較為低洼處,比如研究區(qū)最北端連接長江大橋的地段,該區(qū)域的積水是由與之相鄰區(qū)域的洪水匯集而成。由于所采用地形數(shù)據(jù)精度有限,不能十分精準(zhǔn)地刻畫該區(qū)域城市建筑、道路等,因此該積水分布可作為內(nèi)澇點(diǎn)快速分析的參考。在未來將通過無人機(jī)探測(cè)等手段獲取高分辨率的數(shù)字地表模型(DSM)來構(gòu)建更高精度的內(nèi)澇計(jì)算模型,提供更加準(zhǔn)確的積水風(fēng)險(xiǎn)信息。
圖11 研究區(qū)域不同重現(xiàn)期下最大積水分布(t=6 120 s時(shí)刻積水分布)Fig.11 Peak flood distribution of the study area under different rainfall scenarios
由于研究區(qū)北面長江大橋連接處地勢(shì)較低,易形成較高積水風(fēng)險(xiǎn),將SWMM 計(jì)算所得橋面徑流量(P=0.2%)作為邊界條件輸入LISFLOOD-FP 模型,所得最大積水時(shí)刻的積水分布如圖12 所示。與未加入橋面情況[圖11(f)]相比,此時(shí)研究區(qū)北部的積水范圍有明顯增加,且主要集中在與菜園壩大橋連接處地段,與長江大橋連接處的積水程度也略有增加。因此,在針對(duì)地形變化復(fù)雜的山地城市內(nèi)澇風(fēng)險(xiǎn)排查分析中,應(yīng)對(duì)橋梁排水系統(tǒng)的運(yùn)行情況以及橋面徑流的影響納入考慮。
圖12 最大積水時(shí)刻研究區(qū)積水分布(P=0.2%)Fig.12 Peak flood distribution of the study area
山地城市由于地形起伏變化大,降雨時(shí)空異質(zhì)性明顯,內(nèi)澇的發(fā)生和發(fā)展往往較為迅速,對(duì)較大范圍的雨洪風(fēng)險(xiǎn)點(diǎn)進(jìn)行高效計(jì)算和識(shí)別是預(yù)防和減少災(zāi)害損失的重要前提。研究區(qū)LISFLOOD-FP 網(wǎng)格規(guī)模為92×99個(gè),網(wǎng)格尺寸為30 m×30 m,模型時(shí)間為5 h,模擬時(shí)間步長設(shè)置為1 s。在普通個(gè)人電腦(配置:8 核;Intel Core i7-1065G7 CPU、主頻1.30 GHz)上,耗時(shí)約36 s,而在高性能計(jì)算機(jī)(配置:48 核;Intel Xeon 8275、主頻3.0 GHz)上僅需不到3 s 即可完成計(jì)算,耗時(shí)均為同工況測(cè)試5 次平均結(jié)果??梢?,應(yīng)用該模擬方法可實(shí)現(xiàn)對(duì)大區(qū)域內(nèi)澇風(fēng)險(xiǎn)進(jìn)行快速診斷識(shí)別。
以山地城市重慶某城區(qū)為研究區(qū)域,提出了無管勘資料條件下內(nèi)澇高效模擬計(jì)算方法,采用理想算例和歷史降雨進(jìn)行了驗(yàn)證,并對(duì)區(qū)域內(nèi)不同暴雨情景下的內(nèi)澇積水進(jìn)行了模擬分析。在無管勘資料區(qū)域,利用百度地圖開放平臺(tái)全景靜態(tài)API功能可對(duì)區(qū)域內(nèi)的道路檢查井進(jìn)行有效確定?;贚ISFLOOD-FP,采用雨水井等效排水法和直接降雨法,可實(shí)現(xiàn)對(duì)較大區(qū)域范圍雨洪的高效計(jì)算和對(duì)易澇點(diǎn)進(jìn)行識(shí)別。
積水面積和積水量的變化均與降雨變化較為相似,重現(xiàn)期越大,則同時(shí)刻的積水面積和積水量越大;積水量增長速率最快的時(shí)刻與積水面積達(dá)到峰值的時(shí)刻基本一致,晚于雨峰出現(xiàn)時(shí)間。在積水面積下降階段,積水量仍在增加,陡坡路段可能面臨較大洪水流速,對(duì)行動(dòng)不便的群體仍存在較大威脅,需注意防范。對(duì)于地形起伏變化大的山地城市區(qū)域,橋面徑流可能會(huì)加重局部?jī)?nèi)澇的程度,應(yīng)予以關(guān)注。
所提出的模擬計(jì)算方法可為類似無管勘資料區(qū)域內(nèi)澇風(fēng)險(xiǎn)高效評(píng)估提供技術(shù)參考,對(duì)LISFLOOD-FP 模型在國內(nèi)的使用提供借鑒。