利用海圖數(shù)據(jù)與衛(wèi)星影像計算海灣納潮量——以膠州灣為例＊

李君益，紀(jì)育強(qiáng)，鄭全安，吳永亭，陳義蘭，劉自力，辛海英

（國家海洋局 第一海洋研究所，山東 青島266061）

利用海圖數(shù)據(jù)與衛(wèi)星影像計算海灣納潮量
——以膠州灣為例＊

李君益，紀(jì)育強(qiáng)，鄭全安，吳永亭，陳義蘭，劉自力，辛海英

（國家海洋局 第一海洋研究所，山東 青島266061）

以膠州灣為例，嘗試獲取膠州灣的三維地形，并利用海灣高、低潮所容納的海水體積之差計算海灣納潮量。計算結(jié)果表明：膠州灣1992年納潮量為9.8×108m3，2005年膠州灣納潮量為8.9×108m3，膠州灣納潮能力下降0.9×108m3。

納潮量；Landsat衛(wèi)星影像；海水體積；膠州灣

（張 騫 編輯）

海灣納潮量定義為高潮水量與低潮水量之差，其數(shù)值取決于海灣高潮與低潮水域面積變化、潮差和潮間帶地形。研究人員使用了不同方法對納潮量的計算做了大量工作。鄭全安等［1］利用式（1）計算膠州灣納潮量：

式中，S1，S2分別為高、低潮時刻水域面積；H為其潮差；W為海灣納潮量。

喬貫宇等［2］利用膠州灣口走航ADCP測得的海水流速，對觀測的流量值進(jìn)行時間積分獲得了膠州灣的納潮量。楊世倫等［3］針對膠州灣這種半封閉海灣，潮間帶被部分圍墾后，對納潮量的計算進(jìn)行了探討。陳紅霞等［4］對不同方式納潮量的計算給出了比較。可以看出，上述研究使用的計算方法是對納潮量進(jìn)行直接求解，本研究將利用文章開頭給出的納潮量定義，通過求解海灣高、低潮水量來計算海灣納潮量。

1 納潮量計算方法

理論上，如果已知海灣三維模型，可以計算任意給定潮位海灣所容納的海水體積，進(jìn)而通過計算高、低潮海灣海水體積差來得到海灣納潮量。其中，海灣三維模型的獲取是海灣納潮量計算的關(guān)鍵。

一般情況下，海灣三維模型由實測海灣地形數(shù)據(jù)進(jìn)行插值，并輔以岸界的方式給出。由于進(jìn)行現(xiàn)場實測海灣地形的代價過于昂貴，因此可以采用數(shù)字化海圖的方式獲得0m等深線以下的地形數(shù)據(jù)，但是大部分海圖沒有0m等深線至岸線的潮間帶高程信息。Mason等［5］將水邊線假設(shè)為等高線，疊加多種潮汐條件下的水邊線，生成潮灘數(shù)字高程模型。因此，把提取的多幅海灣衛(wèi)星影像水邊線信息與海圖數(shù)字化后的等深線信息轉(zhuǎn)換到相同坐標(biāo)系與基準(zhǔn)面下，加入岸界信息，可以計算出海灣的三維模型。其中，三維模型的岸界可以由衛(wèi)星影像給出。利用海圖數(shù)據(jù)與衛(wèi)星影像計算海灣納潮量的流程如圖1所示。

本研究以膠州灣為研究實例。由于膠州灣1992年海圖含有潮間帶地形數(shù)據(jù)，因此僅利用海圖數(shù)據(jù)就可以計算該年海灣納潮量；而膠州灣2005年海圖不含有潮間帶數(shù)據(jù)，因此要利用該年份海圖及衛(wèi)星影像計算海灣納潮量。

圖1 利用海圖、衛(wèi)星影像計算納潮量的流程圖Fig.1 A flow sheet for tidal prism computation by using charts and satellite images

2 數(shù)據(jù)來源

1）中國航海圖書出版社出版的膠州灣2005年海圖。該圖采用WGS－84坐標(biāo)系，墨卡托投影，深度基準(zhǔn)為理論深度基準(zhǔn)面。

2）國家海洋局第一海洋研究所測量的1992年膠州灣海圖。該數(shù)據(jù)采用54北京大地坐標(biāo)系，中央經(jīng)線為120°E的3°分帶高斯投影，深度基準(zhǔn)為理論深度基準(zhǔn)面。

3）美國陸地資源衛(wèi)星5號（Landsat－5）TM、美國陸地資源衛(wèi)星7號（Landsat－7）ETM＋影像。影像采用WGS－84坐標(biāo)系，UTM 投影，其中，Landsat－5 4，5，7波段象元大小為30m×30m，Landsat－7 4，5，7波段象元大小為25m×25m。選用1∶50 000地形圖作為理想基礎(chǔ)底圖進(jìn)行幾何精確校正和輻射校正，校正后的幾何精度達(dá)一個象元，其參數(shù)見表1。

表1 Landsat遙感資料參數(shù)表Table 1 Parameters of Landsat images

3 數(shù)據(jù)處理

3.1 海圖矢量化

利用Mapinfo軟件，把膠州灣1992年和2005年掃描海圖的等深線、岸線、水深注記信息由像素形式轉(zhuǎn)化為點和線形式，即矢量化。由于膠州灣沿岸的不斷開發(fā)，1992年海圖矢量化的結(jié)果很難驗證其精度，因此這里只給出2005年矢量化的精度驗證。選用自2005年以來，膠州灣沿岸尚未發(fā)生變化的兩處（分別為紅島參池與黃島碼頭），利用現(xiàn)場GPS定點連續(xù)觀測1h，獲取這2處的空間坐標(biāo)，轉(zhuǎn)換為與海圖投影一致的平面墨卡托x，y坐標(biāo)，其中，x軸對應(yīng)地球赤道，y軸對應(yīng)本初子午線，單位為m；把紅島參池與黃島碼頭的經(jīng)投影轉(zhuǎn)換后的GPS實測坐標(biāo)與矢量化的墨卡托x，y坐標(biāo)相比較（表2），并給出矢量化坐標(biāo)與相應(yīng)GPS實測坐標(biāo)之間的距離差，即矢量化誤差。按照《1∶500、1∶1000、1∶2000地形圖數(shù)字化規(guī)范》［6］和《數(shù)字測繪產(chǎn)品檢查驗收規(guī)定和質(zhì)量評定》［7］中規(guī)定，相對于工作底圖，點狀要素平面位移中誤差不超過圖上的±0.25mm，線狀、面狀要素平面位移中誤差不超過圖上的±0.3mm，該矢量化結(jié)果可信。

表2 GPS紅島、黃島觀測點與海圖上相應(yīng)位置處的矢量化結(jié)果對照表Table 2 Comparison of the GPS data measured at the observation points in Hongdao and Huangdao with the coordinates vectorized from the corresponding positions on the digitalized charts

3.2 水邊線及岸線提取

選用第4，5，7波段進(jìn)行假彩色合成，利用表2中2005年Landsat影像提取水邊線，由潮汐表查得各個影像的過境潮位。研究中的膠州灣水邊線提取標(biāo)準(zhǔn)遵循以下規(guī)定［1］：基巖型海岸線直接用圖像解譯的結(jié)果進(jìn)行對比分析提??；潮汐河口處以河口兩側(cè)的海岸線與河口入?？陂T兩岸交點的連線定義為其海岸線；人工堤壩與周圍地形形成封閉區(qū)的堤壩，以向海的一邊作為海岸線的邊界。提取的海岸線的坐標(biāo)采用WGS－84坐標(biāo)系，6°分帶的 UTM 投影，中央經(jīng)線123°E。膠州灣口取團(tuán)島（120°16′58.188″E，36°02′38.796″N）—薛家島（120°17′14.532″E，36°00′57.528″N）連線（圖2中的膠州灣口處的黑粗虛線部分）。圖2中的1號線為膠州灣岸線（最外側(cè)線段），由膠州灣2005年海圖中提取。提取的2005年水邊線如圖2所示，分別為：2號線為07－28衛(wèi)星過境潮位；3號線為03－06高潮跡線；4號線為05－09衛(wèi)星過境潮位線；5號線為03－06衛(wèi)星過境潮位線；6號線為01－09衛(wèi)星過境潮位線。

圖2 由Landsat衛(wèi)星影像中提取出2005年膠州灣的水邊線Fig.2 The edge water lines of the Jiaozhou Bay extracted from the Landsat images in 2005

3.3 水邊線及岸線數(shù)據(jù)坐標(biāo)的轉(zhuǎn)換

由于1992年海圖數(shù)據(jù)采用北京54坐標(biāo)系，而海圖及衛(wèi)星圖片數(shù)據(jù)采用WGS－84坐標(biāo)系，為了計算結(jié)果有可比性，要把北京54坐標(biāo)系下的數(shù)據(jù)轉(zhuǎn)換到WGS－84坐標(biāo)系。利用3個公共點在WGS－84和北京54坐標(biāo)系下的坐標(biāo)（x54，y54）T、（x84，y84）T，計算出相應(yīng)坐標(biāo)系下的空間直角坐標(biāo)（X54，Y54，Z54）T、（X84，Y84，Z84）T；然后用式（2）坐標(biāo)轉(zhuǎn)換模型［8］求得該地區(qū) WGS－84坐標(biāo)系與北京54坐標(biāo)系的轉(zhuǎn)換參數(shù)：（X0，Y0，Z0）T，m，ω1，ω2，ω3。

式中，（X0，Y0，Z0）T為坐標(biāo)平移量；m為縮放尺度；ω1，ω2，ω3為坐標(biāo)系旋轉(zhuǎn)角。

應(yīng)用上述轉(zhuǎn)換參數(shù)把矢量化后的1992年膠州灣水深及岸線數(shù)據(jù)轉(zhuǎn)換成WGS84大地坐標(biāo)系，墨卡托投影。

3.4 水深數(shù)據(jù)處理

水深數(shù)據(jù)點的數(shù)量有限，而且只能代表相應(yīng)坐標(biāo)處的水深，為了得到膠州灣的水深空間分布，本研究使用空間插值法，求得特定坐標(biāo)處的水深?？臻g插值法有反距離加權(quán)插值法、多項式插值法、Kriging，Natural－Neighbor插值法等。Kriging插值常用來解決空間連續(xù)性變化的屬性非常不規(guī)則的問題［9］，因此這里使用Surfer8.0的Kriging插值方法對膠州灣的實測水深數(shù)據(jù)及海圖矢量化后水深數(shù)據(jù)的進(jìn)行空間插值［10］，插值的同時得到水深網(wǎng)格數(shù)據(jù)，便于計算灣內(nèi)海水體積。為與影像像元尺寸匹配，選取插值網(wǎng)格間距ΔL為30m。

采用交叉檢驗法（Cross－validation）對水深數(shù)據(jù)空間插值結(jié)果進(jìn)行精度檢驗，即假設(shè)某一位置處的水深未知，用其它所有水深數(shù)據(jù)來估算該位置處的水深值。這里進(jìn)行體積計算，因此采用平均誤差來反映插值效果。隨機(jī)抽取800組數(shù)據(jù)進(jìn)行交叉檢驗，水深數(shù)據(jù)插值后的帶來的誤差小于《海洋工程地形測量規(guī)范》［11］規(guī)定測深誤差，因此插值結(jié)果滿足計算要求［12］。

4 膠州灣納潮量的計算

將坐標(biāo)轉(zhuǎn)換后的衛(wèi)星水邊線數(shù)據(jù)及矢量化后的海圖水深數(shù)據(jù)，經(jīng)插值得到膠州灣的三維地形模型，分別得到高、低潮膠州灣海水的空間分布。應(yīng)用式（3）計算膠州灣海水體積V，應(yīng)用式（4）計算膠州灣納潮量P。

式中，V為海水體積；ΔS為水深網(wǎng)格面元，ΔS＝ΔL2；hi為網(wǎng)格i處的水深值；Δh為潮位；V1為高潮時所容納的海水體積；V2為低潮時容納的海水體積。

由潮汐表計算可知，膠州灣1992年平均高潮潮位為3.80m；平均低潮潮位為1.00m；2005年平均高潮潮位為3.79m；平均低潮潮位為1.00m。經(jīng)計算，膠州灣1992年和2005年0m等深線以下海水體積V0、高潮時所容納的海水體積V1、低潮時容納的海水體積V2及納潮量P見表3。

表3 膠州灣納潮量（×108 m3）Table 3 Tidal prism in the Jiaozhou Bay（×108 m3）

5 討 論

本研究計算結(jié)果表明，從1992—2005年，14a間膠州灣0m等深線以下海水體積下降0.6×108m3；高潮時容納的海水體積下降1.8×108m3；低潮容納海水體積下降0.9×108m3，可以看出，膠州灣容納海水的能力總體呈減少的趨勢；而0m線至岸線部分的萎縮是造成膠州灣容納海水能力下降的主要原因；0m線至低潮線部分的減少與低潮線至高潮線部分的減少（即納潮量的減少）對膠州灣容納海水能力的減少均有貢獻(xiàn)。

本研究計算得出的膠州灣1992年納潮量為9.8×108m3、2005年納潮量為8.9×108m3，與鄭全安等［13］計算的膠州灣納潮量相差0.32×108m3（1988年9.48×108m3）；與楊世倫等［2］的計算結(jié)果相差0.918×108m3（1992年10.718×108m3）；與喬貫宇等［3］的計算結(jié)果相差1.217×108m3（2006年春季納潮量為9.522×108m3、秋季納潮量為10.713×108m3，平均納潮量為10.117×108m3）；與胡澤建等①胡澤建，邊淑華，趙可光，等.半封閉海灣淤積災(zāi)害預(yù)測關(guān)鍵技術(shù)研究——以膠州灣為例.國家海洋局第一海洋研究所.2000.的計算結(jié)果相差0.51×108m3（1992年9.29×108m3）。考慮到納潮量的計算時間，本研究結(jié)果與鄭全安等與胡澤建等的結(jié)果相近，原因是本研究建立起的海灣三維模型是對式（1）的微分近似，不同之處是模型注重海灣的三維結(jié)構(gòu)，而鄭全安等與胡澤建等的計算注重對不同潮位海灣面積的測量。

6 結(jié) 論

采用膠州灣衛(wèi)星影像建立起的潮間帶地形數(shù)據(jù)，與矢量化海圖數(shù)據(jù)結(jié)合可以建立膠州灣三維模型，利用該模型可以計算岸線以下任意潮位海灣容納的海水體積，進(jìn)而計算出納潮量。計算結(jié)果表明：1992—2005年膠州灣納潮量減少0.9×108m3，通過與其它納潮量計算方法相比，該計算方便可靠。海灣三維模型對實際潮間帶的模擬精度，取決于衛(wèi)星影像的數(shù)量。

（References）：

［1］ZHENG Q A，WU L Y，DAI M Y，et al.Remote sensing research in the Jiaozhou Bay－1.Measurements of the water area and coastline length［J］.Acta Oceanologica Sinica，1991，22（3）：193－199.鄭全安，吳隆業(yè)，戴懋瑛，等.膠州灣遙感研究－1.總水域面積和總岸線長度量算［J］.海洋與湖沼，1991，22（3）：193－199.

［2］QIAO G Y，HUA F，F(xiàn)AN B，et al.Tidal water capacity computed form ADCP Measurement at Bay Mouth［J］.Advances in Marine Scince，2008，26（3）：282－291.喬貫宇，華鋒，范斌，等.基于 ADCP灣口測流的納潮量計算［J］.海洋科學(xué)進(jìn)展，2008，26（3）：282－291.

［3］YANG S L，CHEN Q M，ZHU J，et a1.Discussion of the tidal prism calculation in a semi－closed bay with the reclaimed fields in the tide—land of Jiaozhou Bay as an example［J］.Ocean Science，2003，27（8）：43－47.楊世倫，陳啟明，朱俊，等.半封閉海灣潮間帶部分圍墾后納潮量計算的商榷—以膠州灣為例［J］.海洋科學(xué)，2003，27（8）：43－47.

［4］CHEN H X，HUA F，LIU N，et al.Comparison among different methods for tidal prism calculation－neap tide of Jiaozhou Bay in autumn 2006as an expounded Example［J］.Advances in Marine Scince，2009，27（1）：11－15.陳紅霞，華鋒，劉娜，等.不同方式的納潮量計算比較—以膠州灣2006年秋季小潮為例［J］.海洋科學(xué)進(jìn)展，2009，27（1）：11－15.

［5］MASON D C，MIN M，DAVENPORT I J，et al.Measurement of recent intertidal sediment transport in Morecambe Bay using the waterline method［J］.Estuarine，Coastal and Shelf Seience，1999，49：427－456.

［6］GB/T17160－1997 Specification for digitsing of 1：500，1：1000and 1：2000topographic maps［S］.Beijing：China Standard Press，1998.GB/T17160－1997 1：500、1：1000、1：2000地形圖數(shù)字化規(guī)范［S］.北京：中國標(biāo)準(zhǔn)出版社，1998.

［7］GB/T18316－2001 Specifications for inspection，acceptance and quality assessment of digital surveying and mapping products［S］.Beijing：China Standard Press，2001.GB/T18316－2001 數(shù)字測繪產(chǎn)品檢查驗收規(guī)定和質(zhì)量評定［S］.北京：中國標(biāo)準(zhǔn)出版社，2001.

［8］KONG X Y，GUO J M，LIU Z Q.Foundation of geodesy［M］.Wuhan：Wuhan University Press，2006.孔祥元，郭際明，劉宗泉.大地測量學(xué)基礎(chǔ)［M］.武漢：武漢大學(xué)出版社，2006.

［9］CHEN H H，LI X，DING W X.Twelve kinds of gridding methods of Surfer8.0in isoline drawing［J］.Chinese Journal of Engineering Geophysics，2007，4（1）：52－57.陳歡歡，李星，丁文秀.Surfer8.0等值線繪制中的十二種插值方法［J］.工程地球物理學(xué)報，2007，4（1）：52－57.

［10］BIAN S H，HU Z J，F(xiàn)ENG A P，et al.An exploration on the natural morphology and the siltation and erosion development of the Jiaozhou Bay in recent 130years［J］.Journal of Oceanography of Huanghai＆Bohai Seas，2001，19（3）：46－53.邊淑華，胡澤建，豐愛平，等.近130年膠州灣自然形態(tài)和沖淤演變探討［J］.黃渤海海洋，2001，19（3）：46－53.

［11］GB17501－1998 Specification for marine engineering topographic surveying［S］.Beijing：China Standard Press，1999.GB17501－1998海洋工程地形測量規(guī)范［S］.北京：中國標(biāo)準(zhǔn)出版社，1999.

［12］YE X M ，JI Y Q ，LIU Z L，et al.Topographic mapping method for Jiaozhou Bay physical model［J］.Coastal Engineering，2008，27（4）：77－82.葉小敏，紀(jì)育強(qiáng)，劉自力，等.膠州灣物理模型地形圖制作方法［J］.海岸工程，2008，27（4）：77－82.

［13］ZHENG Q A，WU L Y，DAI M Y，et al.A remote sensing study of Jiaozhou Bay－Ⅱ.Calculation of dynamical parameters［J］.Acta Oceanologica Sinica，1992，23（1）：1－6.鄭全安，吳隆業(yè)，戴懋瑛，等.膠州灣遙感研究－Ⅱ動力參數(shù)計算［J］.海洋與湖沼，1992，23（1）：1－6.

Calculation of the Tidal Prism in the Jiaozhou Bay by Using Chart Data and Landsat－TM5Images

LI Jun－yi，JI Yu－qiang，ZHENG Quan－an，WU Yong－ting，CHEN Yi－lan，LIU Zi－li，XIN Hai－ying
（FirstInstituteofOceanography，SOA，Qingdao 266061，China）

It is physically sound to compute the tidal prism of a bay by using volume difference between high tide and low tide if a 3Dmodel of the bay is obtained.Herein we take the Jiaozhou Bay as an example and calculate the tidal prism in the bay by using this method.The results show that the tidal prism in the Jiaozhou Bay was 9.8×108m3in 1992and 8.9×108m3in 2005，indicating that a decrease of 0.9×108m3in the tidal prism occurred in the Jiaozhou Bay from 1992to 2005.

tidal prism；Landsat image；Jiaozhou Bay；seawater volume

January 14，2011

P731.24

A

1671－6647（2012）02－0205－07

2011－01－14

海洋公益性行業(yè)科研專項經(jīng)費項目——膠州灣海灣水動力物理模擬實驗研究（200705012）

李君益（1984－），男，遼寧本溪人，碩士研究生，主要從事區(qū)域性海洋學(xué)方面研究.E－mail：ljyi＠fio.org.cn

致謝：感謝張建對GPS數(shù)據(jù)的處理工作。