王露　靳麗萍　李濤

摘要：為了解福建福清市東張水庫淤泥的密度及分布情況，利用基于SILAS系統(tǒng)的淤泥密度測量技術(shù)進行了現(xiàn)場測量和分析。Densitune密度計在3條SILAS測線上均勻采集密度點數(shù)據(jù)。通過選擇不同數(shù)量和位置的密度點數(shù)據(jù)參與SILAS走航數(shù)據(jù)的密度分層計算，來分析最優(yōu)的密度點采集位置和點位數(shù)量。利用SILAS系統(tǒng)和密度計開展淤泥密度的測量和分析工作，可高效獲取測區(qū)的淤泥厚度和分布情況。分析數(shù)據(jù)表明：納入密度分層計算的密度點位淤泥越厚，點位分布越均勻，點位數(shù)量越多，淤泥的密度分層精度越高。研究成果可為水庫及河湖等水域的淤泥研究提供數(shù)據(jù)采集和分析方法。

關(guān)鍵詞：水庫淤泥監(jiān)測;淤泥密度測量;SILAS系統(tǒng);采集與匹配;東張水庫;福建

中圖法分類號：TV698.1

文獻標志碼：A

DOI：10.15974/j .cnki.slsdkb.2019. 10.009

1 研究背景

福建省福清市東張水庫于1958年12月建成蓄水。壩址位于石竹山下的龍江中游，是一座以防洪、農(nóng)業(yè)灌溉、源水供應為主，結(jié)合發(fā)電、旅游等綜合開發(fā)利用的大型水庫。水庫東側(cè)主要為斜陡坡，有大量大孤石;水庫北側(cè)、西側(cè)及南側(cè)西半部分地勢較平緩，植被較茂盛，多以雜草為主，靠近水源部分有沼澤;水庫南側(cè)東半部主要為斜陡坡，由大量小碎石組成。鑒于東張水庫復雜的自然地理條件，對其淤泥密度和厚度的研究成果可推廣應用至同類型水庫中。

為了解水庫多年淤積情況，認識底泥及底質(zhì)的現(xiàn)狀，并將傳統(tǒng)的淤泥厚度量測轉(zhuǎn)化為對淤泥密度的空間分析，加之密度是地形法和輸沙量法計算水庫淤積的重要轉(zhuǎn)換參數(shù)，因此開展庫區(qū)淤泥密度觀測工作非常必要。而淤泥密度觀測對水庫及河湖的淤泥研究及生態(tài)清淤的標準確定也有著重要的意義，水庫底泥是水庫生態(tài)系統(tǒng)的重要組成部分，是入庫物質(zhì)如有機質(zhì)、污染物等的蓄積庫，也是水土界面物質(zhì)進行物理、化學、生物性交換的地帶，因此開展水庫底泥密度的研究，對提升庫區(qū)的有效保護、科學管理起著積極作用[l-3]。S的管子內(nèi)介質(zhì)密度為r）相當于電路中的電感;聲容類似于電容;聲壓類似于電壓;體積速度類似于電流等。基于聲學和電學存在的基本對應關(guān)系，進而可用電路的物理量，通過A/D轉(zhuǎn)換，量化和量測聲波的各個物理量[4-5]。SILAS系統(tǒng)利用普通回聲測深儀向水底發(fā)射低頻聲學信號，聲波到達水底后，部分聲波被反射而另一部分聲波將穿透水底。反射強度和透射強度取決于水底沉積層不同介質(zhì)的密度變化。密度變化對應的是介質(zhì)的特性阻抗變化，根據(jù)上面的內(nèi)容可知，反射強度R、入射角和傳播速度都可以通過電路直接或間接量測得到。因此，可獲得聲波反射強度的比值，即與密度梯度之間的對應關(guān)系，進而可對密度梯度進行定量化處理[6-7]。經(jīng)過量化的密度梯度結(jié)合Densitune密度計獲得的代表點（標定點）絕對密度，建立起絕對密度與密度梯度的之間的對應關(guān)系[8-9]，即獲得了絕對密度與聲波反射強度的對應關(guān)系，從而提取航線上相應密度值曲線。根據(jù)相關(guān)行業(yè)規(guī)范或規(guī)定，確定底泥厚度。SILAS系統(tǒng)工作結(jié)構(gòu)圖如圖1所示。

3 應用實例

3.1 密度儀的標定

現(xiàn)場采集數(shù)據(jù)時直接調(diào)用廠方校正好的密度校正文件是Densitune密度計使用過程中常出現(xiàn)的誤區(qū)，然而因為不同區(qū)域底質(zhì)成分構(gòu)成及流變特性存在差異，調(diào)入不適用的校正文件會導致采集的密度特征點數(shù)值產(chǎn)生誤差。因此，密度儀的標定是SI-LAS系統(tǒng)按密度劃分淤泥的一個關(guān)鍵環(huán)節(jié)。在測量區(qū)域內(nèi)挖取底泥，使用稱重法進行現(xiàn)場淤泥密度標定。標定完成后，采用同樣的方法測量稱重法調(diào)制的固定密度的泥漿，統(tǒng)計精度如表1所示。

由表1可知，Densitune泥漿密度儀的測量中誤差為+5.32 g/L，其相對誤差不大于0.5%的樣本有15個，占比75%，0.5% - 1%的樣本有5個，占比25%，沒有大于1%的樣本。最大誤差百分比為-0.85%，證明儀器具有良好的精度和穩(wěn)定性。

3.2 現(xiàn)場數(shù)據(jù)采集

區(qū)域淤泥密度數(shù)據(jù)的獲取分為兩個步驟，第一步是SILAS系統(tǒng)采集淺地層的簡諧縱波的反射信號，并形成瀑布圖，第二步是采集特征點的密度數(shù)據(jù)[10-13]。

3.2.1 SILAS系統(tǒng)采集淺地層數(shù)據(jù)

以固定架設的GNSS基站為參考站，載有SILAS系統(tǒng)的測量船為流動站，沿著計劃線行進，行進中利用差分定位系統(tǒng)和SILAS系統(tǒng)，連續(xù)獲取測量點的平面位置和高程，并記錄航線上水下地層信息[14]。

3.2.2 Densitune密度計采集特征點密度值

（1）由于可直接影響到淤泥表層及分層數(shù)據(jù)的相對基面深度數(shù)值，密度計采集數(shù)據(jù)前需要進行深度和溫度校正。

（2）密度計特征點位的選擇原則是在SILAS數(shù)據(jù)采集以后，根據(jù)SILAS采集的瀑布圖選擇合適的點位。

（3）密度點位的數(shù)量原則是越多越好，數(shù)量越多，SILAS劃分的密度層越準確。

（4）在具有代表性的庫區(qū)中選擇中間3條線，每條線均勻采集10個點，采用分類對比的方法進行密度值誤差統(tǒng)計分析。

（5）庫區(qū)密度采樣點點位示意圖詳見圖2。

4 成果分析

4.1 庫區(qū)淤泥分布特點

從表2中可以看出，東張水庫淤泥厚度呈現(xiàn)區(qū)域性分布。庫區(qū)北部沼澤及南部地勢平緩區(qū)域淤泥層較厚，最大厚度為0.85 m，平均厚度為0.36 m，總體呈現(xiàn)出南、北區(qū)域比中間淤泥量大的趨勢，也正好符合庫區(qū)南北地勢平緩、沼澤多的特點。

4.2 密度采樣點選擇對密度分層精度的影響

選擇庫區(qū)中間3條SILAS走航測線，每條測線均勻采集10個密度采樣點，選取不同位置和數(shù)量的密度采樣點參與SILAS走航數(shù)據(jù)密度分層計算，未參與分層的密度點作為分層數(shù)據(jù)的外符合精度比對點。SILAS走航數(shù)據(jù)密度分層計算是將采集的密度點的各特征密度（1050，1250 g/cm和1 300 g/cm）分別標定出來，然后將各個密度點的特征密度導人SILAS淺地層數(shù)據(jù)中參與整條測線的淤泥密度分層計算，并獲取測線特征密度（1 050，1 250 g/cm和1 300 g/cm）的相應分層數(shù)據(jù)。本文以國內(nèi)普遍適用的適航水深密度1 300 g/cm3的數(shù)據(jù)為參考數(shù)據(jù)[15]，1 050 g/cm為淤泥表層密度，得到各點1 300 g/cm和1 050 g/cm密度對應的深度差值所代表的淤泥厚度，即為該點的1 300 g/cm密度分層厚度。

以未參與分層密度點為真值的各點誤差曲線見圖3-5。

通過表3的淤泥采樣點中誤差統(tǒng)計表和3條斷面的誤差分布圖可以看出：

（1）參與分層計算的密度點數(shù)量越多，分層計算的密度值越接近于真值（密度儀實測數(shù)值）。

（2）參與分層計算的密度點淤泥厚度越大，分層計算的密度值精度越高。

（3）參與計算的密度點數(shù)量達到一定數(shù)量時（本項目以4個為例），增加參與分層計算的密度點數(shù)量對整體測線密度分層的精度提高效果不明顯。

（4）參與分層計算的密度點分布越均勻，分層計算的密度值精度越高。

由上述的統(tǒng)計結(jié)果得出3種提高密度分層精度的方法：①增加參與分層計算密度點的數(shù)量;②選擇淤泥厚度大的參與分層的密度點;③選擇均勻分布的參與分層的密度點。

5 結(jié)論

本文通過SILAS走航式淤泥測量系統(tǒng)在東張水庫的實際應用研究，得出以下結(jié)論。

（1）利用SILAS淤泥數(shù)據(jù)采集系統(tǒng)開展庫區(qū)的淤泥觀測工作，可以高效地獲取測區(qū)的淤泥厚度和分布情況。實測數(shù)據(jù)表明，東張水庫平均淤泥厚度為0.36 m。

（2）庫區(qū)、河湖等靜態(tài)水域的淤泥分布與地形地勢關(guān)系密切，沼澤、平坦地區(qū)較陡峭地勢區(qū)域淤泥厚度大。

（3）在準確確定淤泥分布的情況下，分析了密度采樣點的選擇對分層精度的影響，提出了一種提高分層精度的可靠的密度采樣點選擇方法，有助于提高走航式SILAS密度數(shù)據(jù)的分層精度。

參考文獻：

[1]謝津平.應用Silas技術(shù)測量太湖湖底淤泥[D].天津：天津大學，2008.

[2]龔聲勝，龔彥，熊偉.SILAS適航水深系統(tǒng)在連云港港口測量中的應用[J].科技信息，2010，12（3）：23-25.

[3]亢文強，謝津平，魯志文.SILAS系統(tǒng)在太湖湖底淤泥調(diào)查中的應用[J"地質(zhì)找礦論叢，2005，20（4）：17-18.

[4]許寶華，王真祥，象山港進港航道外干門淺段試挖槽浮泥觀測研究[J].人民長江，2009，40（22）：29-31.

[5]劉大偉，王真祥，薛劍鋒.基于SILAS適航水深測量系統(tǒng)的試挖航道淤泥實測研究[J].測繪通報，2012（S）：700-703.

[6] 牛桂芝，裴文斌.三爪砣測量適航水深技術(shù)分析與對策[J].水道港口，2006，27（4）：265-268.

[7]牛桂芝，沈小明，裴文斌.SILAS適航水深測量系統(tǒng)測試研究[J].海洋測繪，2003（5）：24-27.

[8]程和森，魏明成，王順中，等，用γ-射線方法觀測甌江口浮泥容重分布和泥沙沉降固結(jié)[R].南京：南京水利科學研究院，2004.

[9]王志東.Silas和Densitune系統(tǒng)在航道工程中的應用[J].港口技術(shù)，2011. 48（1）：10-11.

[10] 錢平生，天津港適航水深下界面淤泥容重值確定與適航水深測量[J].海洋測繪，2002，22（3）：35-39.

[11] 時鐘.杭州灣深水航道浮泥的觀測[J].海洋通報，2001，20（6）：40-50.

[12]陳鈞，萬軍，施衛(wèi)星.雙頻測深儀測深研究[J].海洋測繪，2008，18（2）：32-33.

[13]張俊，顧亞平，查雨，等.雙頻測深儀對淤泥層測定的研究[J]儀器儀表學報，2002，14（2）：16-18.

[14]牛桂芝，裴文斌，三爪砣測量適航水深技術(shù)分析與對策[J].水道港口，2006，27（4）：265-268.

[15]JTJ 319-99疏浚工程技術(shù)規(guī)范[Sl.

（編輯：李曉瀠）