楊 斌 ,趙 睿 ,
(1. 湖北工業(yè)大學(xué),湖北 武漢 430068;2. 長江中游水文水資源勘測局,湖北 武漢 430010)
基于SWMM模型的城市暴雨洪水模擬與分析
楊 斌1,趙 睿1,2
(1. 湖北工業(yè)大學(xué),湖北 武漢 430068;2. 長江中游水文水資源勘測局,湖北 武漢 430010)
為了進(jìn)行城市暴雨洪水模擬分析,以湖北省恩施市中心城區(qū)為研究區(qū)域,在河流水系分布和洪水來源分析的基礎(chǔ)上,提出洪水分析計算方案,構(gòu)建一、二維耦合 SWMM 模型計算典型頻率(10,20,50,100 年一遇)河道洪水演進(jìn)和漫堤淹沒過程,獲得較為準(zhǔn)確的洪水淹沒范圍、水深、淹沒歷時和洪峰到達(dá)時間等信息,并結(jié)合 GIS 進(jìn)行洪水影響評價。結(jié)果表明,清江干流主要河段基本滿足 20 年一遇的防洪標(biāo)準(zhǔn),清江干流和龍洞河部分河段現(xiàn)狀防洪能力較低。SWMM 模型是一種成熟的城市雨洪徑流模型,適宜于城市暴雨洪水模擬與分析應(yīng)用。
暴雨洪水管理模型;一、二維耦合;暴雨洪水模擬;洪水影響分析
伴隨著我國城市化進(jìn)程的不斷加快,城市人口急劇增加導(dǎo)致城市下墊面環(huán)境狀況發(fā)生顯著變化。在全球變化的背景下,城市化進(jìn)程直接影響了流域產(chǎn)匯流過程,以及區(qū)域水資源供、用、耗、排水條件,使得城市小尺度氣候發(fā)生變化。受到下墊面環(huán)境變化和氣候變化等多因素的共同作用,我國城市洪澇災(zāi)害發(fā)生次數(shù)不斷增加,危害程度不斷提升[1]。2012 年 7 月 21 日,北京市出現(xiàn)特大暴雨,全市受災(zāi)人口 160.2 萬人,緊急轉(zhuǎn)移 9.7 萬人,損毀公共設(shè)施448 處,因災(zāi)直接經(jīng)濟(jì)損失高達(dá) 116.4 億元。因此,在城市重點區(qū)域內(nèi),針對可能發(fā)生的暴雨洪水淹沒水深、淹沒范圍、洪峰到達(dá)時間等洪水演進(jìn)過程要素進(jìn)行模擬計算,并結(jié)合 GIS 進(jìn)行洪水影響分析,以標(biāo)示城市范圍內(nèi)受洪水災(zāi)害的危險程度,對于城市防洪除澇減災(zāi)、人口轉(zhuǎn)移疏散、防洪預(yù)案編制等具有重要的意義[2]。
為了解決城市內(nèi)澇和排水等問題,世界各國研究了城市雨洪徑流模型,其中比較典型的是暴雨洪水管理模型(SWMM)、蓄水、處理和溢流模型(STORM)和沃林福德模型(Wallinford)[3-4]。SWMM 模型是在美國環(huán)境保護(hù)署資助下,由美國水資源有限公司、梅特卡夫有限公司和佛羅里達(dá)大學(xué)共同研發(fā)的,提供了地表產(chǎn)流、河道徑流、管網(wǎng)匯流和水污染負(fù)荷的模擬分析,其中地表徑流采用非線性水庫模型,管網(wǎng)匯流采用圣維南方程,水污染負(fù)荷采用累積—沖刷模型[5]。STORM 模型是由美國陸軍工程兵團(tuán)開發(fā)的,主要用于降雨過程中水量、水質(zhì)的模擬和城市蓄水、地表漫流等計算[6]。Wallinford 模型是由英國環(huán)境部開發(fā)的,可進(jìn)行城市水循環(huán)過程的模擬分析、城市雨水、污水排水系統(tǒng)能力模擬[7]。由于 SWMM 模型是提出最早、應(yīng)用最廣泛的城市雨洪徑流模型,本研究選用了該模型進(jìn)行城市暴雨洪水模擬計算。
國內(nèi)外學(xué)術(shù)界對 SWMM 模型參數(shù)敏感性分析、模型適用范圍、模型擴(kuò)展應(yīng)用等相關(guān)研究較為深入,特別是基于 GIS 的 SWMM 模型應(yīng)用成為一大熱點。Zaghloul 等[8]運用神經(jīng)網(wǎng)絡(luò)算法對徑流參數(shù)、匯水區(qū)域參數(shù)等進(jìn)行分析,初步提出了 SWMM 模型參數(shù)適用范圍。Smith 等[9]將 GIS 和 SWMM 相結(jié)合,研究了土地利用類型對降雨徑流的影響。董欣等[10]在運用不確定分析的 HSY 算法和 Monte Carlo采樣法,在實測數(shù)據(jù)的基礎(chǔ)上,對 SWMM 模型參數(shù)進(jìn)行驗證,為模型參數(shù)的確定提供依據(jù)。Krebs 等[11]利用基因遺傳算法研究了加入低影響開發(fā)(LID)措施的現(xiàn)代城市區(qū)域,提出了加入 LID 措施的模型參數(shù)適用范圍。周力寧[12]提出暴雨內(nèi)澇風(fēng)險識別技術(shù),進(jìn)行了不同暴雨強(qiáng)度下內(nèi)澇風(fēng)險等級分析,并基于 SWMM 模型計算內(nèi)澇發(fā)生時積水區(qū)域的積水深度和內(nèi)澇風(fēng)險等級。申邵洪等[13]選擇湖北省保康縣為研究區(qū)域,利用 TUFLOW 軟件構(gòu)建暴雨洪水計算和淹沒分析模型,對城區(qū)受淹地區(qū)、淹沒水深、淹沒面積等進(jìn)行了深入分析。
恩施市位于高山河流沖刷形成的峽谷盆地中,中心城區(qū)地勢低,境內(nèi)河流眾多,水系形態(tài)復(fù)雜,涉河工程多,下游出口斷面窄,城市排水不暢,暴雨山洪引發(fā)的城市內(nèi)澇積水頻發(fā),在山區(qū)城市暴雨洪水分析方面具有一定的代表性。因此,本文選擇恩施市中心城區(qū)為研究區(qū)域,在河流水系分布和洪水來源分析的基礎(chǔ)上,提出洪水分析計算方案,構(gòu)建了一、二維耦合 SWMM 模型,計算典型頻率(10,20,50,100 年一遇)洪水演進(jìn)過程中淹沒范圍、水深、淹沒歷時和洪峰到達(dá)時間等信息,并結(jié)合 GIS 技術(shù)進(jìn)行洪水影響評價,為城市暴雨洪水內(nèi)澇分析和防洪預(yù)案編制提供范例。鑒于恩施市排水管網(wǎng)資料匱乏,為了簡化暴雨洪水過程分析,主要使用 SWMM 模型的地表產(chǎn)流、河道徑流模塊。
恩施市是恩施土家族苗族自治州政府所在地,位于湖北省西南部,清江中上游。全市總面積 3 967 km2,轄 3 個街道、4 個鎮(zhèn)、9 個鄉(xiāng),1 個辦事處,172 個行政村和 34 個社區(qū)。據(jù) 2013 年統(tǒng)計,全市戶籍總?cè)丝?81.43 萬人,其中城鎮(zhèn)人口(含鄉(xiāng)鎮(zhèn)居委會)36.76 萬人,鄉(xiāng)村人口 44.67 萬人。
恩施市城區(qū)主要河流是清江干流,發(fā)源于利川西北部齊躍山東麓的龍洞溝,自北向南從流經(jīng)利川市城區(qū)、騰龍洞、雪照河、恩施市龍鳳鎮(zhèn)、紅廟開發(fā)區(qū)、小渡船、舞陽壩、六角亭進(jìn)入峽口。清江流域面積 2 928.0 km2,在恩施境內(nèi)流長 127.0 km,其中城區(qū)段 18.0 km,主河道坡降 7.0‰。清江干流城區(qū)河段中下游左岸有支流龍洞河、帶水河、蔡家河匯入,下游段峽口右岸有支流高橋河、巴公溪匯入。其中,清江干流進(jìn)入恩施市區(qū)時,流經(jīng)了一個大型水利樞紐工程(大龍?zhí)端畮欤?/p>
龍洞河是清江左岸一級支流,發(fā)源于恩施市大樹埡鄉(xiāng)魚仙洞,流至恩施市舞陽壩官坡與清江交匯,流域面積 38.5 km2,河長 15.3.0 km,主河道坡降 1.8‰。研究區(qū)主要河流水系分布如圖 1 所示。
圖 1 恩施市城區(qū)主要河流水系分布
恩施市所在的清江流域?qū)儆趤啛釒Ъ撅L(fēng)性氣候,降雨充沛。多年平均降水量 1 470.1 mm,最大年降水量 1 918.4 mm。降雨年內(nèi)分配主要集中在 4—10 月,占全年水量的 70% 以上,暴雨多集中在 7—9 月,最大 24 h 降水 326.9 mm。
在恩施市防汛抗旱指揮部辦公室的支持下,系統(tǒng)性搜集了城區(qū)地形、高分辨率遙感影像、水文氣象、社會經(jīng)濟(jì)等數(shù)據(jù),并進(jìn)行了河道斷面測量工作。其中,城區(qū)地形和水文數(shù)據(jù)來源于恩施州水文局,高分辨率遙感影像數(shù)據(jù)來源于 Google Earth 下載的 Worldview 影像,氣象數(shù)據(jù)來源于恩施州氣象局,社會經(jīng)濟(jì)數(shù)據(jù)來源于恩施市統(tǒng)計年鑒(2014 版)。為了獲得更為精細(xì)的河道斷面觀測資料,在湖北省連續(xù)運行衛(wèi)星定位服務(wù)綜合系統(tǒng)(CORS)的支持下,使用 GPS RTK 設(shè)備,在研究區(qū)范圍內(nèi)清江干流和龍洞河河段進(jìn)行了控制斷面測量工作。
根據(jù)恩施市河流水系分析,城區(qū)主要洪水來源為清江干流。在城區(qū)中東部(舞陽壩,恩施州人民政府所在地)附近,有支流龍洞河匯入清江。因此,恩施市主要洪水來源為:1)清江干流洪水;2)龍洞河洪水。
由于恩施市城市防洪標(biāo)準(zhǔn)定為 20 年一遇,因此主要進(jìn)行 10,20,50 ,100 年一遇 4 種洪水量級的洪水分析計算,主要考慮清江干流、龍洞河及相應(yīng)洪水組合,共設(shè)計了 8 種計算方案,如表 1 所示。
表 1 恩施市洪水分析計算方案
采用水動力學(xué)方法進(jìn)行暴雨洪水分析計算。根據(jù)恩施城區(qū)地形、水文氣象及河道斷面觀測資料,采用一維水動力學(xué)模型進(jìn)行河道洪水演進(jìn)分析,采用二維水動力學(xué)模型進(jìn)行河道外城區(qū)洪水淹沒分析。其中,一維水動力學(xué)方法選用圣維南方程,二維水動力學(xué)模型選用深度平均二維淺水波方程。構(gòu)建一、二維耦合模型開展暴雨洪水分析能夠計算出暴雨洪水淹沒范圍和水深、淹沒歷時及洪鋒到達(dá)時間等信息,滿足洪水影響評價的需求。
1)模型范圍。一、二維耦合模型的范圍如圖 1所示,其中藍(lán)色線條代表河道一維水動力學(xué)模型范圍,紅色線條范圍內(nèi)的面狀區(qū)域是二維模型區(qū)域。一維河道在橫向的寬度由沿著左右河岸劃定的線條確定,河岸線內(nèi)是一維河道,河岸線外是二維網(wǎng)格區(qū)域。
模型邊界:清江干流的上邊界范圍為大龍?zhí)端畮煲绾榈莱隹?,下邊界范圍為恩施連珠大橋下游約1 000 m 處;龍洞河的上邊界范圍為恩施東陽酒店上游約 700 m 處,下邊界為龍洞河匯入清江的入口處。
2)一、二維耦合模型。按照 SWMM 模型的設(shè)置要求,進(jìn)行下墊面的概化。將研究區(qū)域下墊面進(jìn)行二維網(wǎng)格劃分,構(gòu)建一、二維耦合 SWMM 模型。其中,一維河道屬于二維網(wǎng)格的非活動區(qū)域,即位于一維河道的網(wǎng)格不參與二維運算,河道結(jié)點處水位由一維河道斷面控制;二維區(qū)域通過網(wǎng)格進(jìn)行運算,模型可將地形(.txt 格式)、初始水位、土地利用等數(shù)據(jù)賦在網(wǎng)格上。
一、二維連接線連接河道結(jié)點和一、二維交界面上,兩控制斷面間河道水位由斷面處水位線性內(nèi)插獲得。當(dāng)河道內(nèi)水位超過一、二維交界面處所在網(wǎng)格的地形高程,則該處網(wǎng)格水位等于河道水位,并以該水位作為二維運算的邊界條件;同時,將二維運算得到的流量作為一維河道的邊界條件。
3)模型邊界條件。對于清江干流洪水分析,上游邊界條件為大龍?zhí)端畮旄黝l率設(shè)計下泄洪水過程線,中間加入龍洞河相應(yīng)洪水匯入。在清江最下游斷面給定實測的河道坡降,以該斷面處的水位流量關(guān)系曲線作為下邊界條件。
對于龍洞河洪水分析,龍洞河上邊界條件為根據(jù)《湖北省暴雨徑流查算圖表》(2012 年)計算的設(shè)計洪水過程線,下邊界條件為清江 20 年一遇洪水水位頂托。
1)模型糙率。通過實地調(diào)研過程中拍攝的照片獲取地面植被情況,參考水力學(xué)計算手冊[14]選取洪水分析模型范圍的糙率系數(shù),其參考糙率值如表 2所示。
2)模型計算參數(shù)。水力學(xué)模型計算時,清江干流各方案計算時間 72 h,二維模型的計算時間步長為 3.0 s,一維模型的計算時間步長為 1.5 s。龍洞河各方案計算時間為 40 h,洪峰值出現(xiàn)在第 10~12 h,二維模型的計算時間步長為 3.0 s,一維模型的計算時間步長為 1.5 s。
表 2 恩施市水力學(xué)計算參考糙率
1)清江干流設(shè)計洪水過程。在一維水力學(xué)模型起點斷面以上的清江干流上,有大龍?zhí)端畮?。因此,進(jìn)入研究區(qū)的清江干流洪水是經(jīng)過大龍?zhí)端畮煺{(diào)蓄后的下泄洪水。根據(jù)《清江大龍?zhí)端麡屑~初步設(shè)計報告》,不同設(shè)計洪水削峰流量如表 3 所示。
表 3 大龍?zhí)端畮煜鞣辶髁勘?/p>
因此,經(jīng)過調(diào)洪演算之后,清江干流設(shè)計洪水過程如圖 2 所示。
圖 2 經(jīng)過大龍?zhí)端畮煺{(diào)蓄后清江干流設(shè)計洪水過程線
2)龍洞河設(shè)計暴雨洪水過程。龍洞河流域內(nèi)無實測洪水資料,且龍洞河流域面積小于 1 000 km2,適合采用《湖北省暴雨徑流查算圖表》(2002 年)推求龍洞河城區(qū)河段控制斷面10,20,50,100 年一遇的設(shè)計洪水過程線。設(shè)計洪水過程結(jié)果如圖 3所示。
圖 3 龍洞河設(shè)計暴雨洪水過程線
以 100 年一遇為例。在清江入口處給定 100 年一遇洪水經(jīng)過大龍?zhí)端畮煺{(diào)蓄后的下泄洪水,龍洞河入口處給定 20 年一遇設(shè)計洪水過程。下游邊界條件給定實測的河道坡降,換算成下游出口的水位流量關(guān)系曲線。
1)洪水淹沒水深、歷時。如圖 4 所示,在 100 年一遇洪水條件下,清江干流洪水淹沒范圍主要集中在恩施市小渡船、黃泥壩、施州大道、東門及南門,平均淹沒水深約為 0.8~1.5 m,平均淹沒歷時約為 3~6 h;龍洞河洪水淹沒范圍主要集中在官坡、棲鳳橋及土橋壩,平均淹沒水深約為 0.5~0.8 m,平均淹沒歷時約為 2~4 h。
圖 4 100 年一遇洪水淹沒計算結(jié)果
2)洪峰到達(dá)時間。以洪水通過恩施城區(qū)的典型控制斷面(恩施市水文站,位于施州大橋上游 100 m處)為基礎(chǔ),計算洪峰到達(dá)時間。
如圖 5 所示,在 100 年一遇洪水條件下,清江干流洪峰到達(dá)時間約為 6~12 h,龍洞河洪峰到達(dá)時間約為 3~6 h。
圖 5 100 年一遇洪峰到達(dá)時間計算結(jié)果
1)洪水淹沒分析。如圖 6 所示,利用 ArcGIS軟件,結(jié)合高分辨率遙感影像數(shù)據(jù)和社會經(jīng)濟(jì)數(shù)據(jù),進(jìn)行恩施市洪水淹沒分析。在 10~100 年一遇洪水條件下,恩施市主要受淹區(qū)域集中在小渡船、黃泥壩、施州大道、東門、南門、官坡、棲鳳橋及土橋壩一帶,受淹區(qū)域總面積約為 2.69 km2,詳細(xì)計算結(jié)果如表 4 所示。
2)現(xiàn)狀防洪能力評價。如圖 7 所示,利用ArcGIS 軟件結(jié)合高分辨率遙感影像和社會經(jīng)濟(jì)數(shù)據(jù),進(jìn)行恩施市現(xiàn)狀防洪能力評價。根據(jù)《山洪災(zāi)害分析評價技術(shù)要求(試行)》文件[15],將現(xiàn)狀防洪能力低于 10 a 的河段定義為極高危險區(qū),現(xiàn)狀防洪能力在 10~20 a 之間的河段為高危險區(qū),現(xiàn)狀防洪能力大于 20 a 的河段為危險區(qū)。清江干流紅旗大橋以上基本上是裸露河床,無河道堤防工程,實際現(xiàn)狀防洪能力低于 20 a(高危險區(qū));紅旗大橋至黃石大橋全部河段均修筑漿砌石堤防,堤防高度在 4~6 m 左右,受到龍洞河洪水頂托作用,實際現(xiàn)狀防洪能力在 20~50 a 之間(危險區(qū));黃石大橋至連珠大橋部分河段修筑漿砌石堤防,堤防高度在 4~6 m 左右,但受到巴公溪、高橋河等洪水頂托作用,實際現(xiàn)狀防洪能力低于 10 a(極高危險區(qū))。龍洞河只有部分下游河段修筑漿砌石堤防,實際現(xiàn)狀防洪能力在10~20 a 之間(高危險區(qū))。
圖 6 恩施市城區(qū)主要洪水淹沒區(qū)
圖 7 恩施市城區(qū)現(xiàn)狀防洪能力評價
1)SWMM 模型是一種成熟的城市雨洪徑流分析模型,能夠較為方便地構(gòu)建一、二維耦合水力學(xué)模型,快速計算河道洪水過程和洪水淹沒分布,適宜于城市暴雨洪水模擬與分析。
2)根據(jù)恩施市城區(qū)河流水系、洪水來源、地形地貌等情況,設(shè)計了一套洪水分析計算方案,并基于 SWMM 模型進(jìn)行雨洪過程模擬分析,取得了較為準(zhǔn)確的洪水淹沒水深、淹沒歷時、洪峰到達(dá)時間等結(jié)果,為城市雨洪淹沒分析和現(xiàn)狀防洪能力評價提供依據(jù)。結(jié)果表明,清江干流主要河段基本滿足了 20 年一遇的防洪標(biāo)準(zhǔn),清江干流和龍洞河部分河段現(xiàn)狀防洪能力較低。這一結(jié)果與恩施市防汛抗旱指揮部辦公室發(fā)布的《恩施市城區(qū)防洪應(yīng)急預(yù)案》[16]基本吻合。
表 4 恩施市雨洪淹沒面積統(tǒng)計表
3)下一步準(zhǔn)備通過重點單位走訪和實地調(diào)研,盡可能地搜集完善恩施市排水管網(wǎng)空間分布、管徑大小等數(shù)據(jù)資料,利用 SWMM 排水管網(wǎng)模塊進(jìn)一步分析城市洪水內(nèi)澇消退過程。
[1] 張冬冬,嚴(yán)登華,王義成,等. 城市內(nèi)澇災(zāi)害風(fēng)險評估及綜合應(yīng)對研究進(jìn)展[J]. 災(zāi)害學(xué),2014 (1): 144-149.
[2] 戴晶晶,劉增賢,陸沈釣. 基于數(shù)值模擬的城市內(nèi)澇風(fēng)險評估研究—以蘇州市城市中心區(qū)為例[J]. 中國水利,2015 (13): 20-23.
[3] 胡偉賢,何文華,黃國如,等. 城市雨洪模擬技術(shù)研究進(jìn)展[J]. 水科學(xué)進(jìn)展,2010 (1): 137-144.
[4] 秦語涵,王紅武,張一龍. 城市雨洪徑流模型研究進(jìn)展[J]. 環(huán)境科學(xué)與技術(shù),2016 (1): 13-19.
[5] CHIRISTOPHER Z. Review of urban storm watermodels[J].Enviornmental Modeling & Software, 2001 (16):195-231.
[6] GORDON S I. Stormwater Runoff Models [M]. U. S. :Springer, 1985: 145-162.
[7] HUANG G, WU S Y. Simulation study on effect of rainwater utilization measures on urban stormwater based on infoworks CS[J]. Water Resources and Power, 2013 (5): 65-71.
[8] ZAGHLOUL N A, M A Abu Kiefa. Neural network solution of inverse parameters used in the sensitivity-calibration analysis of the SWMM model simulations[J]. Advances in Engineering Software, 2001, 32 (7): 87-595.
[9] SMITH D, LI J,BANTING D. A PCSWMM/GIS-based water balance model for the reesor Creek watershed[J].Atmospheric Research, 2005, 77 (1): 338-406.
[10] 董欣,杜鵬飛,李志一,等. SWMM 模型在城市不透水區(qū)地表徑流模擬中的參數(shù)識別與驗證[J]. 環(huán)境科學(xué),2008 (6): 1495-1501.
[11] KREBS G, KOKKONEN T, VALTANEN M, et al. A high resolution application of a storm water management model(SWMM) using genetic parameter optimization[J]. Urban Water Journal, 2013,10 (6): 394-410.
[12] 周力寧. 基于 SWMM 的城市內(nèi)澇風(fēng)險識別研究[M]. 成都:西南交通大學(xué),2016: 35-44.
[13] 申邵洪,江炎生,蔡斌,等. 城鎮(zhèn)區(qū)域洪水風(fēng)險分析研究-以湖北省??悼h為例[J]. 長江科學(xué)院院報,2016,33 (11): 109-115.
[14] 李煒. 水力計算手冊 [M]. 2 版. 北京:中國水利水電出版社,2006: 53-65.
[15] 全國山洪災(zāi)害防治組. 山洪災(zāi)害分析評價技術(shù)要求(試行)[M]. 北京:國家防汛抗旱總指揮部辦公室,2014:79-112.
[16] 恩施州水文水資源勘測局. 恩施市城區(qū)防洪應(yīng)急預(yù)案[M].恩施:恩施市防汛抗旱指揮部辦公室,2009: 33-54.
Urban storm flood simulation and analysis based on storm water management model
YANG Bin1, ZHAO Rui1,2
(1. Hubei University of Technology, Wuhan 430068, China;2. Middler Changjiang River Bureau of Hydrological and Water Resources Survey, Wuhan 430010, China)
In order to conduct urban storm flood simulation and analysis, the paper starts with river system distribution and flood source analysis in case study of Enshi City, and a flood calculation and analysis program are proposed. On the basis of storm water management model (SWMM), an one dimension and two dimension coupled hydraulics model is constructed, which is applied to calculate river flood evolution and diffuse embankment process under the designed flood conditions of 0.01, 0.02, 0.05 and 0.10 (probability). The comprehensive information of flood area, water depth, duration and peak arrival time is also retrieved, which is used to calculate flood effect analysis.Results reveal that the main parts of Qingjiang River basically meet the flood prevention standard of 0.05 (probability)in Enshi City, while a small section of Qingjiang River and Longdong River has lower flood prevention ability than its flood prevention standard. And results also suggest that SWMM is a reliable urban storm water model, which is suitable for urban storm flood simulation and analysis applications.
storm water management model (SWMM); one dimension and two dimension couple; storm flood monitoring; flood effect analysis
TV223
A
1674-9405(2017)05-0056-07
10.19364/j.1674-9405.2017.05.011
2017-06-14
楊 斌(1969-),男,湖南常德人,教授,主要從事管理科學(xué)與工程、智慧城市等方面研究。