白丝美女被狂躁免费视频网站,500av导航大全精品,yw.193.cnc爆乳尤物未满,97se亚洲综合色区,аⅴ天堂中文在线网官网

首頁 / 專利庫(kù) / 地球科學(xué) / 地表徑流 / 一種城區(qū)地表徑流二維數(shù)值模擬方法

一種城區(qū)地表徑流二維數(shù)值模擬方法

閱讀:927發(fā)布:2020-05-13

專利匯可以提供一種城區(qū)地表徑流二維數(shù)值模擬方法專利檢索,專利查詢,專利分析的服務(wù)。并且本 發(fā)明 提供了城區(qū) 地表徑流 二維數(shù)值模擬方法。獲取城區(qū)地形、 地貌 數(shù)據(jù)后,采用三 角 形非結(jié)構(gòu)網(wǎng)格離散城區(qū)計(jì)算區(qū)域, 建筑物 所在區(qū)域不參與剖分,其輪廓線作為網(wǎng)格剖分的控制線,每根建筑物雨 水 管通過匯流連接邊與地面相連。先計(jì)算每根雨水管的匯流面積,再計(jì)算每根雨水管的出流過程。二維網(wǎng)格單元設(shè)置初始計(jì)算條件后啟動(dòng)二維Godunov模型進(jìn)行模擬計(jì)算,雨水管出流過程作為與連接邊對(duì)應(yīng)單元的流量邊界條件,地面網(wǎng)格單元內(nèi)的降雨在模型中當(dāng)作源項(xiàng)處理。本發(fā)明充分考慮城區(qū)內(nèi)建筑物匯流的特點(diǎn),既能保證城區(qū)洪澇計(jì)算時(shí)整體的 質(zhì)量 守恒特性,同時(shí)又能考慮雨水管出口水流動(dòng)量對(duì)地面徑流的影響,克服現(xiàn)有城市洪澇模型無法考慮建筑物匯流過程的缺點(diǎn),提高計(jì)算 精度 。,下面是一種城區(qū)地表徑流二維數(shù)值模擬方法專利的具體信息內(nèi)容。

1.一種城區(qū)地表徑流二維數(shù)值模擬方法,其特征在于:
包括以下步驟:
(1)獲取城區(qū)地形數(shù)據(jù)和樓房、街道地貌幾何特征數(shù)據(jù),樓房地貌幾何特征數(shù)據(jù)包含雨管數(shù)據(jù);
(2)計(jì)算區(qū)域離散:采用三形非結(jié)構(gòu)網(wǎng)格離散城區(qū)計(jì)算區(qū)域,建筑物所在區(qū)域不參與剖分,建筑物輪廓線作為網(wǎng)格剖分的控制邊界線;剖分完成后,建筑物輪廓線被離散作為多個(gè)新生成的三角形單元的單元邊,建筑物雨水管在水平面的投影點(diǎn)落在的單元邊,被定義為連接建筑物和普通城市地面的匯流連接邊;
(3)計(jì)算每根雨水管的樓頂集水面積;
(4)計(jì)算每條匯流連接邊的流量過流:建筑物頂部的降雨徑流通過雨水管對(duì)應(yīng)的匯流連接邊進(jìn)入網(wǎng)格單元,每根雨水管末端的出流流量過程如式(1)所示:
Qk(t)=I(t-τ)*Ak??(1)
其中,Qk(t)為第k根雨水管出口流量過程;t為時(shí)間變量;τ為建筑物頂部到地面的平均匯流時(shí)間,表述為τ=τ1+τ2,τ1為雨水管對(duì)應(yīng)集水區(qū)的平面匯流時(shí)間,τ2為雨水管垂向的匯流時(shí)間;I為降雨強(qiáng)度;Ak為第k根雨水管對(duì)應(yīng)的集水面積;
(5)設(shè)置計(jì)算初始條件:采用地形數(shù)據(jù)對(duì)三角形網(wǎng)格節(jié)點(diǎn)和網(wǎng)格型心進(jìn)行高程插值,獲取網(wǎng)格高程數(shù)據(jù);給每個(gè)網(wǎng)格設(shè)置糙率和初始條件,初始流速和水深均設(shè)置為0;
(6)城區(qū)二維地表徑流計(jì)算:采用基于完整二維淺水方程組的Godunov模型模擬城區(qū)二維地表徑流過程,各條匯流連接邊處的Qk(t)作為各對(duì)應(yīng)單元流量邊界條件處理;普通地面網(wǎng)格單元內(nèi)的降雨作為源項(xiàng)處理,地表匯流過程直接采用二維Godunov模型進(jìn)行模擬。
2.根據(jù)權(quán)利要求1所述的城區(qū)地表徑流二維數(shù)值模擬方法,其特征在于:
步驟(1)中獲取城區(qū)地形數(shù)據(jù)和樓房、街道地貌幾何特征數(shù)據(jù),比例尺為1:2000或1:
2000以上。
3.根據(jù)權(quán)利要求1所述的城區(qū)地表徑流二維數(shù)值模擬方法,其特征在于:步驟(2)計(jì)算區(qū)域離散中,計(jì)算區(qū)域建筑物輪廓線的剖分尺度不超過5m。
4.根據(jù)權(quán)利要求1所述的城區(qū)地表徑流二維數(shù)值模擬方法,其特征在于:如果有建筑物頂部高程數(shù)據(jù),則根據(jù)D8算法計(jì)算每根雨水管的集水面積;如果無樓頂高程數(shù)據(jù),則將建筑物頂部的集水面積平均分配到每根雨水管。
5.根據(jù)權(quán)利要求1所述的城區(qū)地表徑流二維數(shù)值模擬方法,其特征在于:步驟(6)城市二維地表徑流計(jì)算中,所采用的模型為基于完整二維淺水方程組的Godunov模型,該模型所采用的控制方程組的守恒形式如式(2)所示:
其中:
h為水深,u,v分別為x,y方向的流速,t為時(shí)間, 分別為x,y方向的單位向量;
分別為x,y方向的坡度,Zb為地面高程,g為重加速度;
分別為x,y方向的摩阻項(xiàng),其中n為Manning糙
率系數(shù),r為凈雨源項(xiàng),由降雨數(shù)據(jù)經(jīng)過產(chǎn)流計(jì)算得到。
6.根據(jù)權(quán)利要求1所述的城區(qū)地表徑流二維數(shù)值模擬方法,其特征在于:步驟(6)城市二維地表徑流計(jì)算中,在處理連接邊處的流量過程Qk(t)時(shí),直接采用Roe格式的近似Riemann解求解通過匯流連接邊處的數(shù)值通量,公式如式(3)所示:
其中,F(xiàn)*·n為通過匯流連接邊法線方向的數(shù)值通量;h*,u*,v*分別為該連接邊中點(diǎn)處的水深、x方向流速和y方向流速;nx,ny分別為匯流連接邊外法線單位向量在x方向和y方向的分量;g為重力加速度

說明書全文

一種城區(qū)地表徑流二維數(shù)值模擬方法

技術(shù)領(lǐng)域

[0001] 本發(fā)明屬于利工程技術(shù)領(lǐng)域,特別涉及一種考慮建筑物匯流特點(diǎn)的城區(qū)地表徑流二維數(shù)值模擬方法。

背景技術(shù)

[0002] 近幾年,城市暴雨洪澇災(zāi)害呈常態(tài)化現(xiàn)象,成為干擾城市正常運(yùn)行、威脅城市公共安全的重大問題。在城市洪澇問題的研究工作中,城市洪澇模型成為解決城市規(guī)劃、城市洪澇預(yù)警等領(lǐng)域的重要科學(xué)工具。
[0003] 目前,城市洪澇過程的模擬大致來講主要有三種模式:水文模型、簡(jiǎn)化模型和水動(dòng)模型。相比于水文模型和簡(jiǎn)化模型,基于物理機(jī)理的水動(dòng)力模型在模擬具有復(fù)雜城市下墊面條件的水流運(yùn)動(dòng)時(shí)具有更好的準(zhǔn)確性,能夠提供更加翔實(shí)的水力要素信息。這類模型通常把計(jì)算區(qū)域離散成三形(四邊形)非結(jié)構(gòu)網(wǎng)格或正方形網(wǎng)格,通過淺水方程或簡(jiǎn)化的動(dòng)力波、擴(kuò)散波模型來模擬地表的降雨徑流過程?,F(xiàn)有的城市洪澇水動(dòng)力模型方面的分析,較多的是關(guān)注一個(gè)較大區(qū)域的城市洪水淹沒趨勢(shì),因此在網(wǎng)格離散方面選取了相對(duì)較大的尺度,對(duì)城市地貌進(jìn)行了較大的概化處理。然而,由于城市復(fù)雜下墊面條件下的水流運(yùn)動(dòng)極易受到微地形的影響,因此,精細(xì)化是城市洪澇模型發(fā)展的必然趨勢(shì)。
[0004] 城市地表分布有密集的建筑物,這點(diǎn)與普通的自然流域有顯著不同,建筑物的匯流過程與普通城市地面匯流過程有明顯不同,建筑物一般是通過雨水管將房頂?shù)挠晁占蠹信诺降孛?,然后再參與到城市地表的匯流運(yùn)動(dòng)中來。當(dāng)前,尚無城市洪澇模型考慮城區(qū)建筑物與普通城市地面匯流過程的差異,勢(shì)必對(duì)城市洪澇過程的模擬精度造成影響。

發(fā)明內(nèi)容

[0005] 本發(fā)明的目的是提供一種考慮建筑物匯流特點(diǎn)的城區(qū)地表徑流二維數(shù)值模擬方法,該方法區(qū)別于傳統(tǒng)的城市洪澇模型方法,詳細(xì)考慮了建筑物匯流過程的特點(diǎn),首先定義建筑物與普通地面的單元連接邊,再計(jì)算每根雨水管對(duì)應(yīng)的集水面積,然后計(jì)算每根雨水管出口的流量過程,該過程作為連接邊對(duì)應(yīng)三角形單元的入流過程進(jìn)行處理。該方法在保證整體水量守恒特性的同時(shí),還能夠考慮雨水管出口水流動(dòng)量的影響,從而達(dá)到提高城市洪澇模型計(jì)算精度的目的。
[0006] 本發(fā)明的目的是通過以下方案實(shí)現(xiàn)的:
[0007] 本發(fā)明為一種考慮建筑物匯流特點(diǎn)的城區(qū)地表徑流二維數(shù)值模擬方法,在城市精細(xì)化地形、地貌數(shù)據(jù)的基礎(chǔ)上,定義城市建筑物與普通地面的連接邊,確定每根雨水管對(duì)應(yīng)的建筑物頂部的集水面積,計(jì)算每根雨水管出口處的流量過程,該過程作為單元連接邊的流量邊界條件加入到計(jì)算中去。通過該方法能夠很精細(xì)的將建筑物的匯流過程和普通城市地面的匯流過程區(qū)分和聯(lián)系起來,提高城市洪澇模型的計(jì)算精度。該方法的具體步驟如下:
[0008] (1)獲取城區(qū)詳細(xì)的地形數(shù)據(jù)和樓房(含雨水管數(shù)據(jù))、街道等地貌幾何特征數(shù)據(jù)。
[0009] (2)計(jì)算區(qū)域離散:采用三角形非結(jié)構(gòu)網(wǎng)格離散城區(qū)計(jì)算區(qū)域,建筑物所在區(qū)域不參與剖分,其輪廓線作為網(wǎng)格剖分的控制線。剖分完成后,建筑物輪廓線被離散作為多個(gè)新生成的三角形單元的單元邊,建筑物雨水管在水平面的投影點(diǎn)落在哪條單元邊上,哪條邊則被定義為連接建筑物和普通城市地面的匯流連接邊。
[0010] (3)計(jì)算每根雨水管的樓頂集水面積。
[0011] (4)計(jì)算每條匯流連接邊的流量過程:建筑物頂部的降雨徑流通過雨水管對(duì)應(yīng)的匯流連接邊進(jìn)入網(wǎng)格單元,每根雨水管末端的出流流量過程如式(1)所示:
[0012] Qk(t)=I(t-τ)*Ak??(1)
[0013] 其中,Qk(t)為第k根雨水管出口流量過程;t為時(shí)間變量;τ為建筑物頂部到地面的平均匯流時(shí)間,表述為τ=τ1+τ2,τ1為雨水管對(duì)應(yīng)集水區(qū)的平面匯流時(shí)間,τ2為雨水管垂向的匯流時(shí)間;I為降雨強(qiáng)度;Ak為第k根雨水管的集水面積。
[0014] (5)設(shè)置計(jì)算初始條件:采用地形數(shù)據(jù)對(duì)三角形網(wǎng)格節(jié)點(diǎn)和網(wǎng)格型心進(jìn)行高程插值,獲取網(wǎng)格高程數(shù)據(jù);給每個(gè)網(wǎng)格設(shè)置糙率和初始條件,初始流速和水深均設(shè)置為0。
[0015] (6)城區(qū)二維地表徑流計(jì)算:采用基于完整二維淺水方程組的Godunov模型模擬城區(qū)二維地表徑流過程,各條匯流連接邊處的Qk(t)作為各對(duì)應(yīng)單元的流量邊界條件處理,這樣既能保證建筑物所在區(qū)域的徑流質(zhì)量守恒特性,又能考慮雨水管末端水流自身的動(dòng)量影響。普通地面網(wǎng)格單元內(nèi)的降雨作為源項(xiàng)處理,地表匯流過程直接采用二維Godunov模型進(jìn)行模擬。
[0016] 進(jìn)一步優(yōu)選的方案為:步驟(1)中獲取城區(qū)地形數(shù)據(jù)和樓房、街道地貌幾何特征數(shù)據(jù),比例尺為1:2000或1:2000以上。
[0017] 進(jìn)一步的,步驟(2)計(jì)算區(qū)域離散中,為保證數(shù)值計(jì)算精度,計(jì)算區(qū)域建筑物輪廓線的剖分尺度不宜超過5m。
[0018] 進(jìn)一步的,步驟(3)中,如果有詳細(xì)的建筑物頂部高程數(shù)據(jù),則可以根據(jù)D8算法計(jì)算每根雨水管的集水面積;如果無樓頂高程數(shù)據(jù),可將建筑物頂部的集水面積平均分配到每根雨水管。
[0019] 進(jìn)一步的,步驟(6)城市二維地表徑流計(jì)算中,所采用的模型為基于完整二維淺水方程組的Godunov模型,該模型所采用的控制方程的守恒形式如式(2)所示:
[0020]
[0021] 其中:
[0022]
[0023] h為水深,u,v分別為x,y方向的流速,t為時(shí)間, 分別為x,y方向的單位向量;
[0024] 分別為x,y方向的坡度,Zb為地面高程,g為重力加速度;
[0025] 分別為x,y方向的摩阻項(xiàng),其中n為Manning糙率系數(shù),r為凈雨源項(xiàng),由降雨數(shù)據(jù)經(jīng)過產(chǎn)流計(jì)算得到;
[0026] 在模型中,采用Roe格式計(jì)算界面通量,底坡項(xiàng)直接積分處理,摩阻項(xiàng)采用半隱式的處理方法,該模型能夠捕捉城市地表徑流運(yùn)動(dòng)的細(xì)節(jié),如十字路口水流交互發(fā)生的流態(tài)變化以及坡度較大道路處發(fā)生的滾波等。
[0027] 進(jìn)一步的,步驟(6)城市二維地表徑流計(jì)算中,在處理連接邊處的流量過程Qk(t)時(shí),采用Roe格式的近似Riemann解直接求解通過匯流連接邊處的數(shù)值通量,公式如式(3)所示:
[0028]
[0029] 其中,F(xiàn)*·n為通過匯流連接邊外法線方向的數(shù)值通量;h*,u*,v*分別為該連接邊中點(diǎn)處的水深、x方向流速和y方向流速;nx,ny分別為匯流連接邊外法線單位向量在x方向和y方向的分量;g為重力加速度
[0030] 本發(fā)明的有益效果:
[0031] 在本發(fā)明中,充分考慮城區(qū)建筑物匯流過程與城市普通地面匯流過程的不同特點(diǎn),在網(wǎng)格剖分時(shí)即將建筑物與普通地面區(qū)分開來,建筑物所在區(qū)域不參與網(wǎng)格剖分,將雨水管所對(duì)應(yīng)的網(wǎng)格單元邊定義為連接邊(連接建筑物與普通地面),雨水管出口的流量過程作為連接邊對(duì)應(yīng)單元的流量邊界處理,普通網(wǎng)格內(nèi)的降雨作為源項(xiàng)處理,網(wǎng)格內(nèi)的水流運(yùn)動(dòng)采用二維Godunov模型進(jìn)行模擬,這種做法可以充分保證城區(qū)洪澇過程計(jì)算時(shí)整體的水量質(zhì)量守恒特性,同時(shí)又能很好的考慮雨水管出口水流動(dòng)量對(duì)地面徑流的影響。在本發(fā)明中,有區(qū)別的考慮建筑物匯流和城市普通地面匯流過程的不同,可以更好的模擬城區(qū)精細(xì)化的地表匯流過程,提高現(xiàn)有城市洪澇過程模擬的精度。附圖說明
[0032] 圖1為本發(fā)明的城區(qū)地表徑流二維數(shù)值模擬方法流程圖;
[0033] 圖2為本發(fā)明的城區(qū)地表徑流數(shù)二維值模擬方法示意圖。

具體實(shí)施方式

[0034] 下面結(jié)合附圖1和附圖2對(duì)本發(fā)明作進(jìn)一步詳細(xì)說明。
[0035] 本發(fā)明為一種考慮建筑物匯流特點(diǎn)的城區(qū)地表徑流二維數(shù)值模擬方法,充分考慮城區(qū)內(nèi)建筑物匯流與普通地面匯流的區(qū)別與聯(lián)系(即建筑物頂部的降雨徑流過程一般是通過雨水管道流到建筑物周圍的地面上后,再參與到城區(qū)地面匯流的過程中來)來進(jìn)行更精確的城區(qū)地表徑流數(shù)值過程模擬。該方法主要包括以下步驟:
[0036] (1)獲取城區(qū)詳細(xì)的地形數(shù)據(jù)和樓房(含雨水管數(shù)據(jù))、街道等地貌幾何特征數(shù)據(jù),比例尺1:2000以上(含1:2000)。
[0037] (2)計(jì)算區(qū)域離散:采用三角形非結(jié)構(gòu)網(wǎng)格離散城區(qū)計(jì)算區(qū)域,建筑物所在區(qū)域不參與剖分,其輪廓線作為網(wǎng)格剖分的控制邊界線。剖分完成后,建筑物輪廓線被離散作為多個(gè)新生成的三角形單元的單元邊,建筑物雨水管在水平面的投影點(diǎn)落在哪條單元邊上,哪條邊則被定義為連接建筑物和普通城市地面的匯流連接邊。為保證數(shù)值計(jì)算精度,建筑物計(jì)算區(qū)域輪廓線的剖分尺度不宜超過5m。
[0038] (3)計(jì)算每根雨水管的樓頂集水面積:如果有詳細(xì)的建筑物頂部高程數(shù)據(jù),則可以根據(jù)D8算法計(jì)算每根雨水管的集水面積;如果無樓頂高程數(shù)據(jù),可將建筑物頂部的集水面積平均分配到每根雨水管。
[0039] (4)計(jì)算每條匯流連接邊的流量過程:建筑物頂部的降雨徑流通過雨水管對(duì)應(yīng)的匯流連接邊進(jìn)入網(wǎng)格單元,每根雨水管末端的出流流量過程如式(1)所示:
[0040] Qk(t)=I(t-τ)*Ak??(1)
[0041] 其中,Qk(t)為第k根雨水管出口流量過程;t為時(shí)間變量;τ為建筑物頂部到地面的平均匯流時(shí)間,表述為τ=τ1+τ2,τ1為雨水管對(duì)應(yīng)集水區(qū)的平面匯流時(shí)間,τ2為雨水管垂向的匯流時(shí)間;I為降雨強(qiáng)度;Ak為第k根雨水管的集水面積。
[0042] (5)設(shè)置計(jì)算初始條件:采用地形數(shù)據(jù)對(duì)三角形網(wǎng)格節(jié)點(diǎn)和網(wǎng)格型心進(jìn)行高程插值,獲取網(wǎng)格高程數(shù)據(jù);給每個(gè)網(wǎng)格設(shè)置糙率和初始條件,初始流速和水深均設(shè)置為0。
[0043] (6)城區(qū)二維地表徑流計(jì)算:采用基于完整二維淺水方程組的Godunov模型模擬城區(qū)二維地表徑流過程,各條匯流連接邊處的Qk(t)作為各對(duì)應(yīng)單元的流量邊界條件處理,這樣既能保證建筑物所在區(qū)域的徑流質(zhì)量守恒特性,又能考慮雨水管末端水流自身的動(dòng)量影響。普通地面網(wǎng)格單元內(nèi)的降雨作為源項(xiàng)處理,地表匯流過程直接采用二維Godunov模型進(jìn)行模擬。
[0044] 基于完整二維淺水方程組的Godunov模型,該模型所采用的控制方程的守恒形式如式(2)所示:
[0045]
[0046] 其中:
[0047]
[0048] h為水深,u,v分別為x,y方向的流速,t為時(shí)間, 分別為x,y方向的單位向量;
[0049] 分別為x,y方向的坡度,Zb為地面高程,g為重力加速度;
[0050] 分別為x,y方向的摩阻項(xiàng),其中n為Manning糙率系數(shù),r為凈雨源項(xiàng),由降雨數(shù)據(jù)經(jīng)過產(chǎn)流計(jì)算得到;
[0051] 在模型中,采用Roe格式計(jì)算界面通量,底坡項(xiàng)直接積分處理,摩阻項(xiàng)采用半隱式的處理方法,該模型能夠捕捉城市地表徑流運(yùn)動(dòng)的細(xì)節(jié),如十字路口水流交互發(fā)生的流態(tài)變化以及坡度較大道路處發(fā)生的滾波等。具體的模型構(gòu)建過程可以參見如下文獻(xiàn)(張大偉,權(quán)錦,建明,向立.基于Godunov格式的流域地表徑流二維數(shù)值模擬.水利學(xué)報(bào),2018,49(7):787-794.)
[0052] 在處理連接邊處的流量過程Qk(t)時(shí),采用Roe格式的近似Riemann解直接求解通過匯流連接邊處的數(shù)值通量,公式如式(3)所示:
[0053]
[0054] 其中,F(xiàn)*·n為通過匯流連接邊外法線方向的數(shù)值通量;h*,u*,v*分別為該連接邊中點(diǎn)處的水深、x方向流速和y方向流速;nx,ny分別為匯流連接邊外法線單位向量在x方向和y方向的分量;g為重力加速度。
[0055] 上述的實(shí)施例僅是本發(fā)明的部分體現(xiàn),并不能涵蓋本發(fā)明的全部,在上述實(shí)施例以及附圖的基礎(chǔ)上,本領(lǐng)域技術(shù)人員在不付出創(chuàng)造性勞動(dòng)的前提下可獲得更多的實(shí)施方式,因此這些不付出創(chuàng)造性勞動(dòng)的前提下獲得的實(shí)施方式均應(yīng)包含在本發(fā)明的保護(hù)范圍內(nèi)。
高效檢索全球?qū)@?/div>

專利匯是專利免費(fèi)檢索,專利查詢,專利分析-國(guó)家發(fā)明專利查詢檢索分析平臺(tái),是提供專利分析,專利查詢,專利檢索等數(shù)據(jù)服務(wù)功能的知識(shí)產(chǎn)權(quán)數(shù)據(jù)服務(wù)商。

我們的產(chǎn)品包含105個(gè)國(guó)家的1.26億組數(shù)據(jù),免費(fèi)查、免費(fèi)專利分析。

申請(qǐng)?jiān)囉?/a>

QQ群二維碼
意見反饋