高 云 云1,，魯 程 鵬1，王 茂 枚，姜 果，趙 鋼

(1.河海大學(xué) 水文水資源學(xué)院，江蘇 南京 210098； 2.江蘇省水利科學(xué)研究院，江蘇 南京 210017)

水下拋石效果評(píng)價(jià)是拋石護(hù)岸工程驗(yàn)收工作的重要部分，對(duì)水下拋石工程質(zhì)量評(píng)價(jià)有重要意義。此外，水下拋石屬于隱蔽工程，護(hù)岸段及其岸坡主要為第四系松散堆積物，極易沖刷流失[1-2]，應(yīng)更重視水下拋石的施工質(zhì)量檢驗(yàn)。目前評(píng)價(jià)的重要指標(biāo)是測(cè)點(diǎn)增厚值合格率和斷面增厚值合格率，按照順?biāo)鞣较?0～50 m一個(gè)斷面，每個(gè)斷面5～10 m測(cè)一個(gè)點(diǎn)進(jìn)行合格率計(jì)算[3-4]。對(duì)具有隱蔽特性的水下拋石工程而言，第三方的監(jiān)理評(píng)價(jià)決定著工程是否能夠被驗(yàn)收。在實(shí)際工作中，為減輕工作量，常采用40 m一個(gè)斷面進(jìn)行評(píng)價(jià)，斷面的選取對(duì)最終評(píng)價(jià)結(jié)果有著顯著的影響，從而導(dǎo)致工程的監(jiān)理作用大打折扣，不利于發(fā)揮工程的實(shí)際效益[5]。隨著現(xiàn)代電子、計(jì)算機(jī)、信號(hào)處理技術(shù)的發(fā)展，多波束測(cè)深技術(shù)發(fā)生了質(zhì)的飛躍，具有高精度、高密度、高效率和全覆蓋的特點(diǎn)，逐漸被用于河道水下地形測(cè)量，精度可達(dá)0.1 m[6-7]。此時(shí)如仍沿用以往的驗(yàn)收標(biāo)準(zhǔn)，既不能對(duì)實(shí)際拋投狀況進(jìn)行客觀全面的評(píng)價(jià)，還浪費(fèi)了檢驗(yàn)工作所耗費(fèi)的物力和財(cái)力。但多波束測(cè)探技術(shù)所獲得的數(shù)據(jù)點(diǎn)云密度一般過(guò)小[8]，實(shí)際工作中對(duì)精度的要求不高，因此有必要探討一種適合拋石效果評(píng)價(jià)的點(diǎn)云密度，為工程質(zhì)量檢驗(yàn)評(píng)價(jià)提供依據(jù)。

本文以張家港市老海壩河段整治工程中的拋石護(hù)岸工程區(qū)為研究區(qū)，首先，在CARIS軟件上根據(jù)施工期間獲得的多波束水深點(diǎn)云數(shù)據(jù)分別導(dǎo)出不同密度的點(diǎn)云數(shù)據(jù)；然后，對(duì)不同施工區(qū)域前后兩期數(shù)據(jù)進(jìn)行疊加處理獲取拋石增厚值，并分析不同點(diǎn)云密度對(duì)拋石效果評(píng)價(jià)的影響；最后，引入測(cè)量不確定度的概念對(duì)拋石效果評(píng)價(jià)指標(biāo)進(jìn)行不確定性分析。

1 資料與方法

1.1 研究區(qū)概況

研究區(qū)位于張家港市老海壩河段，是長(zhǎng)江下游澄通河段重要的河勢(shì)控制節(jié)點(diǎn)，也是重要的深水航道碼頭集中區(qū)。近幾十年來(lái)，受如皋中汊和瀏海沙水道匯流頂沖影響，河勢(shì)變化劇烈，坍岸嚴(yán)重，經(jīng)拋石治理，現(xiàn)河勢(shì)已基本穩(wěn)定。但近岸處深槽不斷沖刷，出現(xiàn)多個(gè)水深50～70 m的沖刷坑，離岸距離僅150 m，嚴(yán)重威脅河勢(shì)的穩(wěn)定和近岸工礦企業(yè)的防洪安全。為維護(hù)老海壩河段河勢(shì)和航道穩(wěn)定，保障沿線碼頭運(yùn)行安全，自2014年開(kāi)始進(jìn)行節(jié)點(diǎn)整治工程，于2016年完工并驗(yàn)收。驗(yàn)收評(píng)價(jià)時(shí)關(guān)于拋石增厚值的重要指標(biāo)是固定斷面增厚率，由于該指標(biāo)是每隔20～40 m選擇一個(gè)斷面進(jìn)行采樣，故斷面的選取情況對(duì)評(píng)價(jià)結(jié)果的影響較大[5]。

1.2 數(shù)據(jù)來(lái)源

數(shù)據(jù)來(lái)源于節(jié)點(diǎn)整治工程二期3個(gè)標(biāo)段在施工期間(非汛期)的多波束水深數(shù)據(jù)，在施工期間共進(jìn)行7次測(cè)量，從相鄰測(cè)次間施工區(qū)域中分散性地選取12塊50 m×100 m的區(qū)域，然后分別按0.5，1，2，3，4，5，6，7，8，9 m和10 m的密度輸出施工前和施工后對(duì)應(yīng)日期測(cè)次的多波束點(diǎn)云數(shù)據(jù)，并對(duì)每塊區(qū)域施工前后的兩期數(shù)據(jù)進(jìn)行疊加分析(利用ArcGIS平臺(tái))，得到每個(gè)區(qū)域在施工前后不同點(diǎn)云密度下的水深變化值，即測(cè)點(diǎn)增厚值。

1.3 分析方法

1.3.1增厚值的統(tǒng)計(jì)分析

根據(jù)DB32/T2334.2-2013《水利工程施工質(zhì)量檢驗(yàn)與評(píng)定標(biāo)準(zhǔn)》[4]，測(cè)點(diǎn)增厚值大于設(shè)計(jì)值75%的點(diǎn)為合格點(diǎn)，水深流急區(qū)可放寬至大于70%。為提高工作效率，避免選取固定斷面使高密度點(diǎn)云數(shù)據(jù)在不同軟件中切換，根據(jù)相關(guān)文獻(xiàn)[9]，將合格點(diǎn)個(gè)數(shù)占總個(gè)數(shù)的百分比定義為測(cè)點(diǎn)合格率，所有測(cè)點(diǎn)增厚值的平均值定義為平均增厚值，所有測(cè)點(diǎn)增厚值的分布形態(tài)采用偏度和峰度來(lái)體現(xiàn)。利用統(tǒng)計(jì)分析軟件(SPSS)計(jì)算所有區(qū)域不同點(diǎn)云密度下增厚值的統(tǒng)計(jì)特征值和評(píng)價(jià)結(jié)果，如偏度、峰度、測(cè)點(diǎn)合格率、平均增厚值等。

1.3.2誤差分析

根據(jù)施工設(shè)計(jì)，研究區(qū)內(nèi)拋石石塊的粒徑范圍在0.16～0.65 m之間，認(rèn)為0.5 m的點(diǎn)云密度下評(píng)價(jià)結(jié)果接近或等于真實(shí)值，這也是多波束測(cè)量系統(tǒng)可獲取的點(diǎn)云密度。故將該點(diǎn)云密度下的測(cè)點(diǎn)合格率和平均增厚值作為參考值，用于分析其他點(diǎn)云密度評(píng)價(jià)結(jié)果的優(yōu)劣程度。為統(tǒng)一分析所有區(qū)域不同點(diǎn)云密度下評(píng)價(jià)結(jié)果的優(yōu)劣，取測(cè)量學(xué)中相對(duì)誤差的概念[10-11],計(jì)算公式如下：

(1)

式中，δi定義為不同點(diǎn)云密度下的測(cè)點(diǎn)合格率相對(duì)誤差(為了得出誤差正負(fù)，不取絕對(duì)值)；b為0.5 m的點(diǎn)云密度下測(cè)點(diǎn)合格率；bi為其他點(diǎn)云密度下測(cè)點(diǎn)合格率；i取1，2，3，4，5，6，7，8，9，10(平均增厚值相對(duì)誤差的計(jì)算同理)。

1.3.3不確定性分析

測(cè)量不確定度在1993年發(fā)布的《測(cè)量不確定度表示指南》中被首次提出[12]，它是一個(gè)與測(cè)量結(jié)果相關(guān)聯(lián)的參數(shù)，表征合理賦予的被測(cè)量值的分散性，并反映了不同概率水平下可能的誤差分布范圍。測(cè)量不確定度分為2類：標(biāo)準(zhǔn)不確定度和擴(kuò)展不確定度[13-14]。當(dāng)誤差分布為正態(tài)分布時(shí)，對(duì)于某一次測(cè)量i的標(biāo)準(zhǔn)不確定度為uc(即標(biāo)準(zhǔn)偏差)，則其真值Zi落在(Zi-uc,Zi+uc)的可能性為68.26%，即置信概率為68.26%時(shí)測(cè)量結(jié)果的置信范圍為(Zi-uc，Zi+uc)。同理，其擴(kuò)展不確定度為2uc或3uc時(shí)，置信概率則分別為95.40%和99.73%。

據(jù)此，對(duì)相對(duì)誤差引入測(cè)量不確定度的概念，并利用SPSS軟件進(jìn)行檢驗(yàn)，若通過(guò)正態(tài)分布假設(shè)，則可計(jì)算其在不同概率水平下相對(duì)誤差可能的誤差分布范圍[15]。最后利用公式(1)反推不同概率水平下評(píng)價(jià)結(jié)果的置信范圍，具體見(jiàn)公式(2)。

(2)

式中，δmax，δmin分別為相對(duì)誤差的最大值、最小值；bmax，bmin分別為測(cè)量結(jié)果的可能最大值和最小值；bi為評(píng)價(jià)采用點(diǎn)云密度下的評(píng)價(jià)結(jié)果。

2 結(jié)果與分析

2.1 點(diǎn)云密度對(duì)增厚值空間分布的影響

為分析不同點(diǎn)云密度對(duì)增厚值空間分布的影響，選取增厚值分布不均勻的區(qū)域并分別對(duì)0.5，1，5 m和10 m的點(diǎn)云密度下的增厚值進(jìn)行分析。由于不同點(diǎn)云密度下點(diǎn)的密度不一，不能直觀觀察到空間分布的區(qū)別，故對(duì)數(shù)據(jù)點(diǎn)云進(jìn)行克里克金插值處理[16-17]，對(duì)比分析不同點(diǎn)云密度下的增厚值分布情況，結(jié)果見(jiàn)圖1，對(duì)所有區(qū)域不同點(diǎn)云密度下增厚值的偏度和峰度作箱式圖，見(jiàn)圖2～3。在箱圖中，最上方、最下方的粗線分別代表最大值、最小值，箱子的上邊緣、下邊緣和中間線分別代表上四分位數(shù)、下四分位數(shù)和中位數(shù)，最上方和最下方的星號(hào)表示樣本數(shù)據(jù)中的極端值[18]。

結(jié)果初步表明：① 當(dāng)點(diǎn)云密度為1 m時(shí)，拋石增厚值的分布和0.5 m點(diǎn)云密度下的分布并無(wú)明顯差異，細(xì)節(jié)部分也能精確展現(xiàn)；當(dāng)點(diǎn)云密度為5 m時(shí)，增厚值的分布趨勢(shì)和0.5 m點(diǎn)云密度下相同，但局部地方的分布形態(tài)發(fā)生變化；當(dāng)點(diǎn)云密度達(dá)到10 m時(shí)，增厚值的分布已基本不能體現(xiàn)，大部分區(qū)域的增厚值發(fā)生變化，說(shuō)明此時(shí)增厚值的分布形態(tài)已發(fā)生明顯變化，已不能作為最終評(píng)價(jià)的基礎(chǔ)數(shù)據(jù)。② 點(diǎn)云密度的變化對(duì)反映增厚值分布形態(tài)系數(shù)(偏度和峰度)的影響并不顯著，符合工作實(shí)際情況和概率論相關(guān)理論。

圖1 不同點(diǎn)云密度下增厚值空間分布Fig.1 Spatial distribution of thickening values at different point cloud densities

圖2 不同點(diǎn)云密度下偏度分布箱式圖Fig.2 Box diagram of deviation distribution under different cloud densities

圖3 不同點(diǎn)云密度下峰度分布箱式圖Fig.3 Box diagram of kurtosis distribution under different cloud densities

2.2 點(diǎn)云密度對(duì)評(píng)價(jià)結(jié)果的影響

根據(jù)公式(1)，計(jì)算不同點(diǎn)云密度下測(cè)點(diǎn)合格率和平均增厚值的相對(duì)誤差，并對(duì)相對(duì)誤差作箱式圖，得出不同點(diǎn)云密度下相對(duì)誤差的分布情況，見(jiàn)圖4～5。結(jié)果表明：① 當(dāng)點(diǎn)云密度為1 m時(shí)，兩者的相對(duì)誤差均分布在±2%之間，可作為精確評(píng)價(jià)分析時(shí)的點(diǎn)云密度；② 當(dāng)點(diǎn)云密度小于5 m時(shí)，測(cè)點(diǎn)合格率和平均增厚值均逐漸增加；③ 點(diǎn)云密度大于5 m時(shí)，測(cè)點(diǎn)合格率和平均增厚值的相對(duì)誤差分布范圍逐漸穩(wěn)定，大部分?jǐn)?shù)據(jù)均分布在±8%之間，但其離散趨勢(shì)有所增加；④ 該圖證實(shí)了點(diǎn)云密度越小結(jié)果越精確的普遍規(guī)律，結(jié)果具備普適性。

圖4 不同點(diǎn)云密度下測(cè)點(diǎn)合格率的相對(duì)誤差分布箱式圖Fig.4 Relative error distribution box plot of qualified points under different point cloud densities

圖5 不同點(diǎn)云密度下平均增厚值的相對(duì)誤差分布箱式圖Fig.5 Relative error distribution box plot of average thickening values under different point cloud densities

2.3 點(diǎn)云密度對(duì)評(píng)價(jià)結(jié)果的不確定性影響

2.3.1假設(shè)檢驗(yàn)

由于實(shí)際工作以合格率作為驗(yàn)收標(biāo)準(zhǔn)，現(xiàn)對(duì)不同點(diǎn)云密度下的測(cè)點(diǎn)合格率的相對(duì)誤差引入測(cè)量不確定度的概念，利用統(tǒng)計(jì)分析軟件(SPSS)對(duì)其進(jìn)行正態(tài)分布檢驗(yàn)，結(jié)果見(jiàn)表1。由表1可知：除個(gè)別情況外，在顯著性水平為0.05時(shí)，不同點(diǎn)云密度下測(cè)點(diǎn)合格率的相對(duì)誤差均呈正態(tài)分布，符合自然界誤差呈正態(tài)分布的規(guī)律；隨著點(diǎn)云密度的增加，標(biāo)準(zhǔn)偏差呈現(xiàn)明顯的增加趨勢(shì)，平均值則先減小后增大。

2.3.2不確定性分析

由上節(jié)可知，增厚值合格率的相對(duì)誤差呈正態(tài)分布，引入測(cè)量不確定度的概念，依據(jù)1.3.3節(jié)的方法計(jì)算不同概率水平下合格率的可能范圍，結(jié)果見(jiàn)圖6。

圖6 不同置信概率下測(cè)點(diǎn)合格率可能分布范圍Fig.6 The possible distribution range of qualified rate under different confidence probabilities

由圖6可知：① 當(dāng)置信概率相同時(shí)，點(diǎn)云密度為5 m時(shí)的相對(duì)誤差范圍始終小于點(diǎn)云密度為10 m時(shí)的相對(duì)誤差范圍；② 在相同置信概率下，隨著測(cè)量合格率的增加，合格率可能分布的范圍也逐漸增加；③ 點(diǎn)云密度為5 m時(shí)，當(dāng)測(cè)量得到的合格率為75%時(shí)，在置信概率分別為95.44%和68.26%條件下，合格率真值的可能分布范圍分別為65%～83%和70%～78%；點(diǎn)云密度為10m時(shí)，當(dāng)測(cè)量得到的合格率為75%時(shí)，在置信概率分別為95.44%和68.26%條件下，合格率真值的可能分布范圍分別為60%～87%和67%～80%，顯然點(diǎn)云密度為5m時(shí)區(qū)間較窄，評(píng)價(jià)結(jié)果精度更高，結(jié)果的可靠性也較高；④ 在進(jìn)行拋石效果分析時(shí)，如遇重要施工區(qū)，如河流水勢(shì)急劇變化區(qū)域，應(yīng)選擇精度較高的點(diǎn)云數(shù)據(jù)(如5 m)，在河勢(shì)較平穩(wěn)區(qū)域及水下地形沖淤變化不明顯的區(qū)域，可選擇點(diǎn)云密度為10 m的數(shù)據(jù)進(jìn)行拋石效果評(píng)價(jià)。

3 結(jié) 論

本文通過(guò)對(duì)選取區(qū)域的拋石增厚值作誤差分析并引入測(cè)量不確定度理論進(jìn)行研究，可得到如下結(jié)論。

(1) 點(diǎn)云密度對(duì)河勢(shì)急劇變化區(qū)的拋石效果評(píng)價(jià)結(jié)果的影響較為顯著，評(píng)價(jià)結(jié)果的誤差隨點(diǎn)云密度的增加呈先增加后逐漸穩(wěn)定的變化趨勢(shì)。

(2) 在置信概率相同的條件下，點(diǎn)云密度越低，其拋石效果評(píng)價(jià)結(jié)果的不確定性越小，相應(yīng)的可靠程度越大。

研究不僅證明了評(píng)價(jià)結(jié)果的誤差呈正態(tài)分布的普遍規(guī)律，還給出了不同點(diǎn)云密度下評(píng)價(jià)結(jié)果的不確定性，為河勢(shì)急劇變化區(qū)水下拋石評(píng)價(jià)選取合適的點(diǎn)云密度提供了借鑒。