王淑蘭， 葛于席

(1.江陰市江港堤閘管理處， 江蘇 無(wú)錫 214431； 2.江蘇全方測(cè)繪地理信息有限公司， 江蘇 無(wú)錫 214431)

數(shù)字高程模型，簡(jiǎn)稱(chēng)DEM，是通過(guò)有限的地形高程數(shù)據(jù)實(shí)現(xiàn)對(duì)地面地形的數(shù)字化模擬，它是用有序數(shù)值陣列形式表示地面高程的實(shí)體地面模型，是數(shù)字地形模型的分支。DEM可以描述包括高程在內(nèi)的各種地貌因子，如坡度、坡向、坡度變化等因子在內(nèi)的線性和非線性組合的空間分布，其中DEM是零階單純的單項(xiàng)數(shù)字地貌模型，其他如坡度、坡向及坡度變化率等地貌特性可在DEM的基礎(chǔ)上派生[1-2]。

本文主要分析DEM在水下地形中的應(yīng)用，以長(zhǎng)江下游江陰丁壩為闡述對(duì)象，通過(guò)對(duì)長(zhǎng)江河道水下地形測(cè)繪成果的對(duì)比分析，掌握水下地形的變化規(guī)律與趨勢(shì)，做出相應(yīng)的應(yīng)對(duì)措施，為江堤及節(jié)點(diǎn)控制工程丁壩的維護(hù)和安全提供數(shù)據(jù)參考，同時(shí)以點(diǎn)帶面地詮釋DEM在描述水下地形演變中的重要作用。

1 丁壩地理區(qū)位

長(zhǎng)江江陰段西起與常州交界的老桃花港，東至與張家港交界的陸家莊，全長(zhǎng)35 km。由于江陰長(zhǎng)江南岸臨江山體基巖露頭，無(wú)山體段地表以下土質(zhì)堅(jiān)硬，在20世紀(jì)80年代后，人工護(hù)岸工程的加強(qiáng)，河道主流趨于穩(wěn)定，長(zhǎng)期保持少變的趨勢(shì)。尤其是江陰桃花港段長(zhǎng)江水道主流偏南，深泓緊靠南岸一側(cè)，岸線微向江內(nèi)凸起，水面在此收窄、水流變急，河床易遭洪水的沖刷，此處建有人工節(jié)點(diǎn)丁壩控制調(diào)整水流，和下游鵝鼻嘴節(jié)點(diǎn)一起有效控制了長(zhǎng)江主弘的擺動(dòng)和走向，使江陰段長(zhǎng)江河勢(shì)基本穩(wěn)定。

2018年丁壩處桃花港1.5 km長(zhǎng)江河道被定為Ⅲ級(jí)崩岸險(xiǎn)情監(jiān)測(cè)預(yù)警段，每年須組織開(kāi)展江陰桃花港1.5 km長(zhǎng)江河道Ⅲ級(jí)崩岸險(xiǎn)情預(yù)警監(jiān)測(cè)工作。長(zhǎng)期水流的沖涮和往來(lái)船只的影響，因水下地理?xiàng)l件的不同，造成的河沙沉積及對(duì)河道堤防的沖刷狀況也不盡相同。

通過(guò)2013年、2016年2個(gè)年度的數(shù)據(jù)的DEM模型計(jì)算與分析，闡明水下地形監(jiān)測(cè)的必要性與丁壩對(duì)水下地形的影響。

2 數(shù)據(jù)與分析

2.1 數(shù)據(jù)格式說(shuō)明

水下數(shù)據(jù)采用DGPS+測(cè)深儀+潮汐改正，組合成數(shù)據(jù)采集工具，平面坐標(biāo)系為1954北京坐標(biāo)系，高程基準(zhǔn)為1985黃海高程，為了數(shù)據(jù)分析的準(zhǔn)確性，所測(cè)得數(shù)據(jù)不進(jìn)行任何內(nèi)插值。原始數(shù)據(jù)記錄格式通過(guò)ArcGIS軟件中的Tin模塊，實(shí)現(xiàn)水下地形的數(shù)字化模擬。

數(shù)據(jù)提取2013年10月和2016年12月水下地形測(cè)量成果進(jìn)行對(duì)比分析，數(shù)據(jù)橫斷面點(diǎn)間距為20 m，斷面間距為40 m。建立DEM分析該段水域的河床泥沙沖刷與淤積、坡度、河床粗糙度等方面的變化情況，建立丁壩區(qū)域水下傾斜柵格圖，可見(jiàn)河底“深槽”。

2.2 泥沙沖刷與淤積變化

可以通過(guò)2期DEM數(shù)據(jù)的柵格運(yùn)算來(lái)分析泥沙的沖刷與淤積，判斷河床未來(lái)演變方向以及上下游防洪安全，分析河床的演變態(tài)勢(shì)。

由于丁壩的存在改變了原有的水動(dòng)力外部條件，至2016年12月溝脊長(zhǎng)度超200 m，溝底最寬處約200 m。通過(guò)2次DEM柵格計(jì)算，在丁壩壩頭長(zhǎng)江側(cè)河底地形變化較小，2個(gè)年度高程值相近，趨于穩(wěn)定，溝脊處變化較大，沖刷位置向上游移動(dòng)?！吧畈邸碧幙v斷面顯示縱向沖刷線平均向上游移動(dòng)約8.3 m，縱軸同一數(shù)值明顯向左移動(dòng)(上游)；橫斷面數(shù)據(jù)顯示“深槽”向江岸側(cè)移動(dòng)大約2.0 m以上，溝逐漸變窄、變深。

2.3 坡度變化分析

坡度是用來(lái)描述河床地形變化起伏的一系列數(shù)值，與河沙的淤積、流失相關(guān)聯(lián)。坡度數(shù)值越大，數(shù)值之間變化越大，說(shuō)明河床越陡峭，水流動(dòng)力等各方面因數(shù)復(fù)雜；反之則河床平緩穩(wěn)定，水流也相對(duì)平緩，對(duì)河床、堤岸等影響較小。

通過(guò)2年數(shù)據(jù)對(duì)比如表1，該區(qū)域0°～12°為主要地勢(shì)，合計(jì)占比80% 以上，地形連續(xù)變化較大，對(duì)比數(shù)值主要變化在2°～ 5°和5°～ 12°，變化主要發(fā)生在“橢圓形”區(qū)域內(nèi)(深槽附近)?！吧畈邸苯斗较颍炊螀^(qū)影響域內(nèi)，坡度為0°～ 2°，該地形對(duì)長(zhǎng)江堤岸、河床等周邊環(huán)境影響較小；2°～ 12°變化量最大，正是河床土壤分離和淤積的主要區(qū)域；丁壩盡頭長(zhǎng)江方向坡度在22°～ 27°區(qū)域內(nèi)，變化數(shù)值不大，但需注意后期的觀測(cè)數(shù)據(jù)。

表1 坡度占比范圍 單位：%

2.4 河床粗糙度變化

地面粗糙度是指在某個(gè)區(qū)域內(nèi)， 柵格表面積與其投影面積之比，是反映地表形態(tài)的一個(gè)宏觀指標(biāo)，是水下地形表達(dá)的參數(shù)之一，也是河床穩(wěn)定性參考因子之一，粗糙度越大意味地貌表面越粗糙[3]。

粗糙度通過(guò)地形的柵格坡度值可以計(jì)算出，公式如下:

M=1/cos(β3.14159/180)

(1)

式中，M為粗糙度，β為柵格坡度值,在 ArcGIS中，通過(guò)提取坡度得到的值是角度，而計(jì)算默認(rèn)的角度值是弧度值，所以在計(jì)算時(shí)必須把角度轉(zhuǎn)為弧度[4]。

使用ArcGIS中柵格計(jì)算器，輸入1/cos(“坡度柵格”×3.14159/180)，得到2013年和2016年測(cè)量區(qū)域柵格粗糙度數(shù)據(jù)，最大值分別為1.16597、1.17419，最小值均為1。由數(shù)值分析結(jié)合坡度圖可知，水下地形粗糙度在緩慢增加，柵格的粗糙度可能會(huì)隨水下地形采樣點(diǎn)的密度不同而不同，可以看出柵格表面的粗糙程度緩慢變化，和DEM柵格運(yùn)算得出沖刷和流失區(qū)域相符。

2.5 DEM在長(zhǎng)江水下地形中的分析結(jié)果

沖淤變化：用丁壩2期的水下地形數(shù)據(jù)作對(duì)比分析，對(duì)河床泥沙的流失與淤積的統(tǒng)計(jì)、計(jì)算、分析，可以看出沿丁壩方向側(cè)水下地形有緩慢的變化，水下沖刷槽緩慢向上游移動(dòng)。

坡度分析：可以看出水下地形坡度變化主要在2°～5°、5°～ 12°這兩個(gè)較平坦的區(qū)域，其余坡度范圍變化程度不大。平坦區(qū)域占比較大，說(shuō)明該區(qū)域在丁壩的作用下水下泥沙交換較為平衡，變化相對(duì)較為平穩(wěn)。

粗糙度分析：通過(guò)粗糙度分析該段河床的粗糙度，該區(qū)域主要為平滑區(qū)域，不平滑區(qū)域主要在“深槽”和丁壩區(qū)域。變化數(shù)值較大的主要在丁壩附近，主要是泥沙的淤積，對(duì)岸線和江堤的影響較小，但需注意丁壩附近的沖刷槽變化情況。

3 結(jié) 語(yǔ)

觀測(cè)數(shù)據(jù)處理采用DEM技術(shù)，可以及時(shí)、直觀發(fā)現(xiàn)水下地形的變化，進(jìn)而可以采取科學(xué)手段早干預(yù)、早預(yù)防，確保江堤、水上及水下安全，這對(duì)長(zhǎng)江防洪、水上運(yùn)輸，航行路線，都有著重大的實(shí)際意義。針對(duì)該技術(shù)在應(yīng)用分析過(guò)程中還存在影響精度的一些因素，比如采集的數(shù)據(jù)點(diǎn)間距，長(zhǎng)江的水況、環(huán)境等，在以后的實(shí)施過(guò)程中可以通過(guò)減小斷面間距和采集間隔，換用多波束測(cè)深系統(tǒng)等來(lái)提高數(shù)據(jù)的質(zhì)量，以達(dá)到準(zhǔn)確分析的預(yù)期，為保障江堤、岸線的安全提供直觀可靠的數(shù)據(jù)支持。