趙飛達(dá)，董睿，吳琳，宋積文，姜靜，張大勇

1.中海油信息科技有限公司，北京 100027

2.大連理工大學(xué) 海洋科學(xué)與技術(shù)學(xué)院，遼寧 盤錦 124221

渤海具有豐富的油氣資源，我國(guó)渤海的遼東灣由于其地理、氣象和水文條件，歷年都存在海域覆冰現(xiàn)象[1]。典型的氣象條件如寒潮、過(guò)冷空氣等均會(huì)對(duì)海域覆冰有顯著影響。由于冬季海域結(jié)冰問(wèn)題，海洋平臺(tái)勘探、海上航運(yùn)等工作也會(huì)受到不同程度的影響[2]。歷史上因輕視海冰威脅而發(fā)生較多安全事故，如1969 年的渤海大冰封事件[3]，海冰漂移產(chǎn)生巨大推力導(dǎo)致平臺(tái)坍塌，對(duì)平臺(tái)上部設(shè)施有嚴(yán)重?fù)p壞，同時(shí)對(duì)作業(yè)人員造成安全隱患；2000 年遼東灣某油氣平臺(tái)因海冰作用造成上部管線斷裂，引發(fā)安全事故。嚴(yán)重的冰情對(duì)各類海上作業(yè)都存在一定程度的威脅[4-5]。

為規(guī)范歷年冰情研究，中國(guó)海洋石油總公司（China National Offshore Oil Corporation，CNOOC）在20 世紀(jì)初制定發(fā)布了《中國(guó)海海冰條件及應(yīng)用規(guī)定》[6]（后文簡(jiǎn)稱《規(guī)定》）。然而，近年來(lái)渤海的海冰條件發(fā)生了變化，傳統(tǒng)的標(biāo)準(zhǔn)和經(jīng)驗(yàn)無(wú)法應(yīng)對(duì)新的海冰形勢(shì)。大連理工大學(xué)海冰管理項(xiàng)目組開(kāi)發(fā)了一套完整的海冰監(jiān)測(cè)技術(shù)，于20 世紀(jì)末投入應(yīng)用，為適應(yīng)海冰新形勢(shì)與時(shí)代新的技術(shù)，不斷更新升級(jí)[7-8]。本文利用JZ20-2 點(diǎn)式平臺(tái)監(jiān)測(cè)的海冰數(shù)據(jù)，分析了2012—2022 年的海冰狀況，為更新和修訂渤海新的海冰標(biāo)準(zhǔn)提供依據(jù)。

1 海冰定點(diǎn)監(jiān)測(cè)系統(tǒng)

每年入冬，渤海的部分海域會(huì)出現(xiàn)結(jié)冰現(xiàn)象，隨著冷空氣不斷入侵，這些海域結(jié)冰情況愈發(fā)嚴(yán)重，海冰覆蓋面積明顯增大，冰厚不斷加厚，次年春季時(shí)海冰開(kāi)始融化，直至消失。海冰的生消包括初冰期、盛冰期與融冰期3 個(gè)時(shí)期，在海冰監(jiān)測(cè)中現(xiàn)場(chǎng)監(jiān)測(cè)是一種傳統(tǒng)且關(guān)鍵的監(jiān)測(cè)手段，但在冰厚、冰速、冰密集度、冰類型等海冰要素的現(xiàn)場(chǎng)監(jiān)測(cè)方面仍舊存在很多困難[9-10]。由于海冰參數(shù)測(cè)量結(jié)果的缺失，早期的冰荷載原型測(cè)量實(shí)驗(yàn)數(shù)據(jù)不能得到有效的解釋，對(duì)冰荷載的研究帶來(lái)巨大的挑戰(zhàn)，所以在研究過(guò)程中需使用適合海域環(huán)境且能準(zhǔn)確獲得海冰要素的方法[11]。同時(shí)，冰情變化也受環(huán)境條件的影響，故對(duì)現(xiàn)場(chǎng)的溫度、風(fēng)速等氣象數(shù)據(jù)進(jìn)行監(jiān)測(cè)，氣象數(shù)據(jù)由成熟完備的小型氣象站準(zhǔn)確監(jiān)測(cè)?；谀壳暗难芯壳闆r，以JZ20-2 海域的海洋平臺(tái)為研究目標(biāo)，海冰管理項(xiàng)目組設(shè)計(jì)海冰定點(diǎn)監(jiān)測(cè)系統(tǒng)，研發(fā)海冰要素監(jiān)測(cè)技術(shù)。

1.1 海冰要素監(jiān)測(cè)技術(shù)

海冰要素監(jiān)測(cè)主要為冰厚、冰速和冰密集度的監(jiān)測(cè)。海冰管理項(xiàng)目組依據(jù)圖像測(cè)距原理，設(shè)計(jì)開(kāi)發(fā)了一套基于視頻圖像技術(shù)的海冰要素監(jiān)測(cè)系統(tǒng)，并已應(yīng)用于渤海的海冰定點(diǎn)監(jiān)測(cè)中。

1.1.1 測(cè)量基本原理

1）冰厚、冰速測(cè)量

海冰冰厚測(cè)量方法包括目測(cè)法、實(shí)地測(cè)量法、遙感測(cè)量法等。實(shí)地測(cè)量可以通過(guò)人工取樣、圖像測(cè)量等方式開(kāi)展，遙感測(cè)量可以通過(guò)雷達(dá)、聲吶、衛(wèi)星等方式開(kāi)展。綜合考慮設(shè)備的測(cè)量精度、現(xiàn)場(chǎng)設(shè)備安裝可靠性等問(wèn)題，研究采用圖像測(cè)量法測(cè)量冰厚，該方法的原理為

式中：L為物體的實(shí)際長(zhǎng)度，l為物體的成像長(zhǎng)度，H為攝像機(jī)的透鏡與物體的距離，h為攝像機(jī)的透鏡與成像的距離。

在現(xiàn)場(chǎng)圖像測(cè)量應(yīng)用中，在攝像機(jī)焦距不變的情況下，首先標(biāo)定攝像頭，通過(guò)式（1）計(jì)算出成像比例n，由此推導(dǎo)被測(cè)物體的實(shí)際長(zhǎng)度S與其成像長(zhǎng)度s的關(guān)系為

則被測(cè)物體的實(shí)際長(zhǎng)度S為

海冰與結(jié)構(gòu)發(fā)生相互作用后出現(xiàn)破碎，攝像機(jī)對(duì)掉落至海中的碎冰塊進(jìn)行圖像采集，碎冰的斷面參數(shù)被捕獲，前期根據(jù)焦距不變對(duì)攝像機(jī)進(jìn)行標(biāo)定后，即可通過(guò)標(biāo)定物的尺寸與海冰斷面像素點(diǎn)尺寸來(lái)計(jì)算冰厚，以上為圖像法測(cè)量冰厚的基本原理。海冰翻轉(zhuǎn)現(xiàn)場(chǎng)圖片如圖1 所示，以錐體的防浪孔為標(biāo)定（如圖1（a）中AB），該孔的孔徑HD為20 cm，圖像上像素長(zhǎng)度為Dh，海冰斷面像素尺寸為T（如圖1（b）），則實(shí)際冰厚Tice為

圖1 定點(diǎn)監(jiān)測(cè)時(shí)測(cè)量冰厚

此方法可以用來(lái)測(cè)量冰厚、冰速等海冰要素[12]。

2）海冰密集度測(cè)量

海冰的覆蓋率表明結(jié)冰區(qū)海冰的密集度，通過(guò)圖像分割方法確定海水和海冰的占比，來(lái)完成海冰密集度的圖像識(shí)別。首先預(yù)處理圖像；其次采用射影方法對(duì)圖像進(jìn)行校正；最后通過(guò)最大類別方差法實(shí)現(xiàn)海水和海冰的分割，以確定海冰的占比面積。

其中，在海洋平臺(tái)上采集海冰圖像時(shí)，會(huì)因?yàn)榕臄z角度出現(xiàn)射影失真，所以校正工作非常重要。與此同時(shí)，海冰圖像識(shí)別時(shí)容易受海平面的干擾，如圖2 所示，為避免誤差，在識(shí)別時(shí)添加海平面提取步驟，將圖像中天空部分刪除，以此實(shí)現(xiàn)對(duì)海冰圖像的單獨(dú)處理，從而有效避免將天空判定為海冰的情況。

圖2 海冰密集度圖像識(shí)別

1.1.2 海冰參數(shù)測(cè)量軟件

在海冰參數(shù)測(cè)量與數(shù)據(jù)處理的過(guò)程中，基于圖像測(cè)量原理，如果人工參與測(cè)量過(guò)程較多，會(huì)導(dǎo)致測(cè)量結(jié)果存在誤差，且工作量過(guò)大。因此，項(xiàng)目組設(shè)計(jì)開(kāi)發(fā)了海冰參數(shù)自動(dòng)測(cè)量軟件，用于實(shí)現(xiàn)海冰參數(shù)測(cè)量的半自動(dòng)化，保證精確性，減少人力消耗。通過(guò)圖像射影矯正、圖像識(shí)別與匹配等方法對(duì)現(xiàn)場(chǎng)視頻進(jìn)行分析與計(jì)算處理，計(jì)算所得結(jié)果與測(cè)量數(shù)據(jù)將同步存儲(chǔ)。歷年冰情分析依托于該軟件所測(cè)得的海冰參數(shù)數(shù)據(jù)，軟件界面如圖3 所示，今后將通過(guò)現(xiàn)場(chǎng)實(shí)際應(yīng)用不斷完善。

圖3 海冰參數(shù)測(cè)量軟件界面

2 JZ20-2 海域歷年冰情統(tǒng)計(jì)分析

JZ20-2 隸屬的渤海海域位于北緯37°～41°，是全世界緯度最低的結(jié)冰海域，每年冬季都會(huì)出現(xiàn)覆冰現(xiàn)象，歷年冰情受氣候影響存在顯著不同。冰情嚴(yán)重年份對(duì)海上航運(yùn)、石油開(kāi)采、捕魚(yú)養(yǎng)殖等相關(guān)活動(dòng)具有很大影響，甚至直接影響整個(gè)環(huán)渤海經(jīng)濟(jì)發(fā)展[13-14]。海冰的監(jiān)測(cè)、預(yù)報(bào)及研究工作能在一定程度上有益于海上作業(yè)，本文以JZ20-2 海域監(jiān)測(cè)系統(tǒng)的實(shí)測(cè)數(shù)據(jù)為基礎(chǔ)，參考?xì)鉁亍⒈?、冰厚、海冰覆蓋率等資料數(shù)據(jù)，統(tǒng)計(jì)近10 a 中JZ20-2 海域的冰情變化趨勢(shì)。

2.1 冰期數(shù)據(jù)統(tǒng)計(jì)分析

JZ20-2 海域一年中僅冬季出現(xiàn)覆冰現(xiàn)象，每年最早出現(xiàn)海冰的日期為初冰日，翌年海冰消融的日期為終冰日。從每年的初冰日到終冰日這段時(shí)期稱為該年的總冰期，總冰期內(nèi)海冰會(huì)經(jīng)歷從生成到發(fā)展再到相對(duì)穩(wěn)定，直至最后消融的過(guò)程。隨海冰的發(fā)展，總冰期又可分為初冰期、盛冰期以及終冰期3 個(gè)階段，其中海冰進(jìn)入盛冰期的日期為盛冰日，進(jìn)入終冰期的日期為融冰日。將近10 a 中JZ20-2 海域的冰情情況進(jìn)行統(tǒng)計(jì)，如表1 所示，分析結(jié)果如圖4 所示。受環(huán)境影響，會(huì)出現(xiàn)總冰期內(nèi)海域不結(jié)冰的現(xiàn)象。據(jù)統(tǒng)計(jì)，無(wú)冰日約占總冰期的1/10～1/5。因此，有效冰期的統(tǒng)計(jì)能更加精確有效地體現(xiàn)冰情。

表1 近10 a 中JZ20-2 海域冰期

據(jù)統(tǒng)計(jì)，JZ20-2 海域初冰日通常在初年12 月中下旬或翌年1 月上旬，于翌年2 月或3 月終冰。受冰情變化影響，歷年冰期不等。由圖4 可以看出，近10 a 中，JZ20-2 海域總冰期在40～80 d內(nèi)變動(dòng)。總冰期于2012—2013 年際達(dá)到最長(zhǎng)，為76 d，最短為2016—2017 年際，僅有41 d；有效冰期最長(zhǎng)為72 d，最短為33 d，分別為2013—2014年際與2016—2017 年際。受歷年氣溫具體變化的差異影響，有效冰期小于總冰期2～14 d 不等。2012—2017 年冰期整體呈縮短趨勢(shì)，2017—2019 年冰期受氣溫影響有所波動(dòng)，近4 a 逐漸產(chǎn)生回升趨勢(shì)。

圖4 近10 a 中JZ20-2 海域冰期統(tǒng)計(jì)

2.2 冰厚數(shù)據(jù)統(tǒng)計(jì)分析

國(guó)家海洋局根據(jù)中國(guó)海的結(jié)冰特點(diǎn)，依據(jù)結(jié)冰范圍和冰厚資料于1973 年研究制定了《中國(guó)海冰情預(yù)報(bào)等級(jí)》，將中國(guó)海冰情共劃分5 級(jí)，并規(guī)定這5 個(gè)等級(jí)作為海冰預(yù)報(bào)標(biāo)準(zhǔn)，如表2[15]。

表2 冰情預(yù)報(bào)等級(jí)

冰厚是歷年冰情分析的關(guān)鍵參數(shù)，現(xiàn)將近10 a 中JZ20-2 海域覆冰后的冰厚進(jìn)行統(tǒng)計(jì)，數(shù)據(jù)以單層冰為主，如表3 所示，分析結(jié)果如圖5 所示。

表3 近10 a 中JZ20-2 海域冰厚 cm

分析圖5 可知，近10 a 中JZ20-2 海域冰厚變化整體呈下降趨勢(shì)。2012—2013 年際最大冰厚為30 cm，此后年際最大冰厚均未超過(guò)20 cm。歷年平均冰厚在3～11 cm 不等，最大時(shí)出現(xiàn)在2012—2013 年際，為10.7 cm；最小時(shí)出現(xiàn)在2018—2019年際，為3.3 cm。根據(jù)冰情預(yù)報(bào)等級(jí)，2012—2013年際冰厚較大，屬于偏輕冰年；其余各年際最大冰厚均小于20 cm，平均冰厚均小于10 cm，屬于輕冰年。

圖5 近10 a 中JZ20-2 海域冰厚統(tǒng)計(jì)

2.3 海冰覆蓋率數(shù)據(jù)統(tǒng)計(jì)分析

海冰的覆蓋率表明結(jié)冰區(qū)域海冰的密集度，根據(jù)不同年份JZ20-2 海域覆冰面積程度的不同將冰量劃分為10 個(gè)等級(jí)?，F(xiàn)將近10 a 中JZ20-2海域的海冰覆蓋率進(jìn)行統(tǒng)計(jì)，如表4 所示，分析結(jié)果如圖6 所示。

表4 近10 a 中JZ20-2 海域海冰覆蓋率 %

圖6 近10 a 中JZ20-2 海域海冰覆蓋率統(tǒng)計(jì)

經(jīng)統(tǒng)計(jì)，近10 a 中JZ20-2 海域最大覆冰率多達(dá)八成及以上，平均最大覆冰率為84%，有時(shí)可實(shí)現(xiàn)海域全面覆冰。平均覆冰量為整個(gè)冰期每日覆冰量的平均值，可反映某冰年整體覆冰情況。近10 a 中JZ20-2 海域平均覆冰率最大可達(dá)88.1%，最低為42.1%，分別為2012—2013 年際與2020—2021年際。歷年平均覆冰率多浮動(dòng)于60%～70%內(nèi)，近2 a 最大覆冰率與平均覆冰量均呈下滑狀態(tài)，海域整體覆冰量有下滑趨勢(shì)。

2.4 氣溫?cái)?shù)據(jù)統(tǒng)計(jì)分析

JZ20-2 隸屬的渤海海域?yàn)榈湫偷募竟?jié)性結(jié)冰海域，受西伯利亞高壓影響，盛行寒冷干燥的西北風(fēng)。每年冬季，渤海海域在過(guò)冷空氣不斷入侵下開(kāi)始結(jié)冰，有時(shí)爆發(fā)寒潮會(huì)導(dǎo)致覆冰面積增大，海冰加厚。氣候變化直接影響著該海域海冰的生消與冰情強(qiáng)弱。據(jù)統(tǒng)計(jì)，渤海海域在冬季的平均氣溫在-3 ℃左右，1 月是最為嚴(yán)寒的月份，平均氣溫僅-7 ℃?，F(xiàn)將近10 a 中JZ20-2 海域平均氣溫的變化情況進(jìn)行統(tǒng)計(jì)，分析結(jié)果如圖7。分析圖7 得出，JZ20-2 海域近10 a 中冬季平均氣溫變化范圍為-5.5～-2.0 ℃。平均氣溫最低為2012 —2013年際，低至-5.5 ℃；2013 —2014、2014—2015、2016—2017 及2017—2018 這4 年際平均氣溫均不低于渤海海域冬季平均氣溫-3℃，分別為-2.5、-2、-3 及-2.8 ℃。結(jié)合冰期、冰厚與覆冰率數(shù)據(jù)，發(fā)現(xiàn)氣溫對(duì)海域結(jié)冰有明顯影響。受全球變暖影響，近10 a 中JZ20-2 海域冬季氣溫整體呈回暖現(xiàn)象，又因歷年不同過(guò)冷空氣與寒潮影響有所波動(dòng)。自2014—2015 年際始，平均氣溫有所下降且伴隨著小范圍回升。

圖7 近10 a 中JZ20-2 海域平均氣溫變化

2.5 實(shí)測(cè)數(shù)據(jù)與標(biāo)準(zhǔn)中數(shù)據(jù)的對(duì)比分析

《規(guī)定》中對(duì)于冰日、冰期、冰厚、覆蓋率以及氣溫等數(shù)據(jù)給出規(guī)范性建議[6]，本文將近10 a中JZ20-2 海域的冰情數(shù)據(jù)分析結(jié)果與《規(guī)定》進(jìn)行對(duì)比。

對(duì)比結(jié)果發(fā)現(xiàn)，在冰日方面，初冰日、盛冰日、融冰日推遲，平均融冰日、終冰日提前，對(duì)比結(jié)果如表5 所示；在冰期方面，初冰期增長(zhǎng)，盛冰期、融冰期縮短，總冰期縮短，有效冰期增長(zhǎng)，對(duì)比結(jié)果如表6 所示；在冰厚方面，單層冰厚度降低，海冰危害降低，對(duì)比結(jié)果如表7 所示；在盛冰期海冰覆蓋率方面，最大覆蓋率提高，最小覆蓋率降低，平均覆蓋率小幅降低，冰情嚴(yán)重情況不好確定，對(duì)比結(jié)果如表8 所示；在氣溫方面，冰期平均氣溫整體升高，冰情情況與氣溫呈正相關(guān)，對(duì)比結(jié)果如表9 所示。

表5 冰日對(duì)比表

表6 冰期對(duì)比 d

表7 單層冰厚對(duì)比 cm

表8 盛冰期海冰覆蓋率對(duì)比 %

表9 冰期氣溫對(duì)比 ℃

3 結(jié)論

本文基于海冰監(jiān)測(cè)系統(tǒng)將近10 a 中JZ20-2 海域的冰期、冰厚、覆冰率、氣溫?cái)?shù)據(jù)進(jìn)行統(tǒng)計(jì)，結(jié)合《規(guī)定》，分析近10 a 的冰情情況，得到結(jié)論如下：

1）統(tǒng)計(jì)近10 a 中JZ20-2 海域冰期數(shù)據(jù)發(fā)現(xiàn)，2012—2017 年冰期整體呈縮短趨勢(shì)，2017—2019年冰期受氣溫影響有所波動(dòng)，近4 a 逐漸產(chǎn)生回升趨勢(shì)。

2）統(tǒng)計(jì)近10 a 中JZ20-2 海域冰厚數(shù)據(jù)發(fā)現(xiàn)，冰厚整體變化情況呈下降趨勢(shì)，2012—2013 年際冰厚較大，屬于偏輕冰年；其余各年際最大冰厚均小于20 cm，平均冰厚均小于10 cm，屬于輕冰年。

3）統(tǒng)計(jì)近10 a 中JZ20-2 海域海冰覆蓋率數(shù)據(jù)發(fā)現(xiàn)，覆冰率整體呈下降趨勢(shì)，最大覆冰量平均在6 成以上。

4）統(tǒng)計(jì)近10 a 中JZ20-2 海域氣溫?cái)?shù)據(jù)發(fā)現(xiàn)，平均氣溫整體呈上升趨勢(shì)，低于-3 ℃的年際冰情嚴(yán)重，冰情受氣溫影響較大。

本文對(duì)JZ20-2 海域近10 a 的冰情數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，總結(jié)變化規(guī)律，為渤海海冰相關(guān)標(biāo)準(zhǔn)的更新與完善提供數(shù)據(jù)支持。