孟丁丁 陳曉東 季順迎

(大連理工大學(xué)工業(yè)裝備結(jié)構(gòu)分析國家重點實驗室，遼寧大連 116024)

海冰力學(xué)性質(zhì)對海洋結(jié)構(gòu)物的設(shè)計和海冰動力學(xué)分析有重要的影響。對海冰力學(xué)性質(zhì)的研究包括海冰的彎曲強(qiáng)度、壓縮強(qiáng)度、剪切強(qiáng)度和斷裂韌度等[1-3]，其中海冰彎曲強(qiáng)度是海洋平臺、斜面護(hù)坡結(jié)構(gòu)、破冰船與海冰作用產(chǎn)生冰載荷的主要影響因素[4-7]。因此，對海冰彎曲強(qiáng)度的研究有重要的意義。

海冰彎曲強(qiáng)度受冰晶結(jié)構(gòu)、溫度、鹽度、海冰密度等物理因素的影響而不同，一般通過三點彎曲試驗、四點彎曲試驗和現(xiàn)場懸臂梁試驗測試確定[8]。冰晶結(jié)構(gòu)對海冰彎曲強(qiáng)度的影響主要體現(xiàn)在冰晶形狀、海冰加載時的冰晶方向等方面。關(guān)于海冰溫度對彎曲強(qiáng)度的研究表明，一定溫度范圍內(nèi)，溫度越高彎曲強(qiáng)度值越低[9-11]。對溫度和鹽度的綜合研究，則以鹵水體積的形式呈現(xiàn)[12]。大量的彎曲試驗數(shù)據(jù)表明鹵水體積的平方根和彎曲強(qiáng)度呈指數(shù)關(guān)系[13-14]。國內(nèi)外研究學(xué)者對影響海冰彎曲強(qiáng)度的因素作進(jìn)一步探討[10-11,13]，發(fā)現(xiàn)對海冰的加載速率也對其彎曲強(qiáng)度有不同影響。然而，以上研究均是基于單因素或者兩個因素對海冰彎曲強(qiáng)度的影響進(jìn)行的，海冰受多種因素的共同影響的彎曲強(qiáng)度還未有研究。

近年來機(jī)器學(xué)習(xí)的發(fā)展為多因素影響下輸出結(jié)果的預(yù)測研究提供了可能，其中較為典型的BP神經(jīng)網(wǎng)絡(luò)[15]早在20世紀(jì)80年代被提出，為了更好地解決具體領(lǐng)域的特定問題衍生出循環(huán)神經(jīng)網(wǎng)絡(luò)（recurrent neural network，RNN）[16]，LSTM[17]、CNN[18]、RBF[19]等神經(jīng)網(wǎng)絡(luò)模型。其中RNN模型克服了傳統(tǒng)神經(jīng)網(wǎng)絡(luò)對輸入和輸出的許多限制，對輸出結(jié)果的預(yù)測不僅受到當(dāng)前參數(shù)的影響，而且也參考前一步的預(yù)測結(jié)果。這也在不同工程研究領(lǐng)域得到了廣泛應(yīng)用，李澤龍等[20]利用LSTM-RNN模型對多參數(shù)影響下鐵水硅含量進(jìn)行預(yù)測得到了預(yù)測誤差穩(wěn)定、預(yù)測精度高的訓(xùn)練模型。聶少軍等[21]在風(fēng)洞試驗中采用RNN對激波風(fēng)洞測力系統(tǒng)振動特性進(jìn)行了高效建模，基本消除了天平輸出信號中的大幅振動干擾信號，取得了較為理想的效果。此外，RNN模型在自然語言處理[22]、文本識別[23]、圖像音頻[24]等領(lǐng)域均取得了驚人的成績。機(jī)器學(xué)習(xí)在關(guān)于海冰的物理力學(xué)性質(zhì)的研究中也得到了應(yīng)用。Kellner等[25]利用主成分分析法及決策樹方法建立了小型冰試驗數(shù)據(jù)庫，通過冰參數(shù)預(yù)測冰的峰值應(yīng)力及行為。然而RNN方法在海冰力學(xué)性質(zhì)的預(yù)測研究中還未得到應(yīng)用。

本文結(jié)合渤海海冰試驗得到的數(shù)據(jù)，對海冰彎曲強(qiáng)度的影響因素進(jìn)行分析，研究基于RNN預(yù)測海冰彎曲強(qiáng)度的模型。該模型實現(xiàn)了多參數(shù)影響下海冰彎曲強(qiáng)度的預(yù)測，僅需輸入海冰的溫度、鹽度、密度及加載速率即可快速得到結(jié)果。該結(jié)果的預(yù)測精度相較于普遍采用的鹵水體積的平方根和彎曲強(qiáng)度關(guān)系擬合得到的公式預(yù)測結(jié)果有了較大提高，可為海洋斜面結(jié)構(gòu)的設(shè)計載荷提供參考。

1 海冰彎曲強(qiáng)度試驗及影響因素分析

海冰彎曲強(qiáng)度的試驗研究在海冰和結(jié)構(gòu)物相互作用過程中有重要地位，主要有簡支梁和懸臂梁等試驗方法。Timco等[9]、Schwarz等[26]和李志軍[27]均對海冰彎曲試驗進(jìn)行了研究，分析了不同試驗條件和物理因素對海冰彎曲強(qiáng)度的影響。本文依據(jù)海冰試驗規(guī)范進(jìn)行試驗，對環(huán)渤海海冰的彎曲強(qiáng)度試驗進(jìn)行總結(jié)，給出了試驗方法及相關(guān)影響因素分析。

1.1 海冰彎曲強(qiáng)度的試驗方法及數(shù)據(jù)采集

為研究渤海海域的海冰物理力學(xué)性質(zhì)，2009?2010年冬季在環(huán)渤海海岸進(jìn)行了海冰彎曲強(qiáng)度試驗[5]。為得到較全面的渤海海冰數(shù)據(jù)，選擇采冰地點分布如圖1所示。在9個測點進(jìn)行了145組海冰彎曲強(qiáng)度測試。

圖1 渤海海冰彎曲強(qiáng)度測試位置分布Fig. 1 Location distribution of sea ice bending strength test in Bohai Sea

試驗采用三點彎曲的試驗方法，將海冰切割成75 mm×75 mm×700 mm的長方體，為探究海冰物理性質(zhì)對彎曲強(qiáng)度的影響，分別測量了海冰的鹽度、溫度和密度。由于在與船舶、海洋平臺等海洋結(jié)構(gòu)物作用時海冰發(fā)生彎曲破壞的方向多為海冰生成方向，為了更接近工程實際，試驗時選擇沿此方向進(jìn)行加載。圖2為海冰三點彎曲示意圖，其中P為載荷，L為試樣長度，L0為加載跨距，b為試樣寬度，h為試樣高度。由簡支梁相關(guān)理論，當(dāng)加載點位于跨中時，可求解得到海冰試樣的彎曲強(qiáng)度，即

圖2 海冰三點彎曲試驗加載示意圖Fig. 2 Loading diagram of sea ice three-point bending test

式中，σf為海冰彎曲強(qiáng)度，Pmax為海冰彎曲破壞時的最大載荷。海冰三點彎曲試驗中記錄了加載過程中的受力時程。典型的載荷歷程曲線如圖3所示。該時程曲線很好地體現(xiàn)了海冰脆性破壞特性。

圖3 海冰彎曲試驗的載荷時程曲線Fig. 3 Load time history curve of sea ice bending test

1.2 海冰彎曲強(qiáng)度影響因素分析

海冰彎曲強(qiáng)度受冰晶結(jié)構(gòu)、溫度、鹽度、密度以及加載速率的影響。早期研究中主要針對單個因素影響下的海冰彎曲強(qiáng)度進(jìn)行分析，得到各個影響因素和彎曲強(qiáng)度的函數(shù)關(guān)系[5,8-10,13]。為了確定渤海海冰各參數(shù)對海冰彎曲強(qiáng)度的影響，本文依據(jù)環(huán)渤海海冰彎曲強(qiáng)度的測試結(jié)果對海冰彎曲強(qiáng)度的影響因素進(jìn)行逐個分析。

1.2.1 溫度對海冰彎曲強(qiáng)度的影響

溫度是海冰冰晶結(jié)構(gòu)形成的重要影響因素，也是鹵水體積的主要參數(shù)，對海冰彎曲強(qiáng)度有顯著影響[9]。國內(nèi)外研究表明，海冰彎曲強(qiáng)度σf同溫度T呈線性關(guān)系[10,28]。本文對2009?2010年在渤海海冰彎曲試驗的數(shù)據(jù)進(jìn)行分析，在4 ppt～6.7 ppt鹽度范圍內(nèi)，得到溫度和彎曲強(qiáng)度的關(guān)系如圖4所示。對此進(jìn)行線性擬合得到

圖4 海冰溫度與彎曲強(qiáng)度的對應(yīng)關(guān)系Fig. 4 The relation between sea ice temperature and bending strength

1.2.2 鹽度對海冰彎曲強(qiáng)度的影響

海冰的鹽度是海冰融化后海水的鹽度。在空氣溫度較低時，海水凍結(jié)成海冰的速度較快，海水內(nèi)鹽分來不及析出即凍結(jié)在海冰中，因此形成的海冰是冰晶固體和鹵汁的混合物[29]。在–5 ℃ ～–15 ℃溫度范圍內(nèi)，基于海冰彎曲強(qiáng)度的實測數(shù)據(jù)擬合后得到彎曲強(qiáng)度與鹽度的線性關(guān)系如圖5所示。鹽度越大，彎曲強(qiáng)度越小。這是由于鹽度越高越容易結(jié)冰，從而導(dǎo)致海冰結(jié)構(gòu)中含有較多的鹵水及雜質(zhì)，致使海冰彎曲強(qiáng)度較小。

圖5 海冰鹽度與彎曲強(qiáng)度的對應(yīng)關(guān)系Fig. 5 The relation between sea ice salinity and bending strength

1.2.3 鹵水體積對海冰彎曲強(qiáng)度的影響

海冰的溫度和鹽度可綜合表征為鹵水體積，其關(guān)系公式為[30]

式中，vb為海冰鹵水體積， 1 0?3；S為海冰鹽度，10?3；T為海冰溫度，℃。

Timco等[9]總結(jié)了鹵水體積的平方根與彎曲強(qiáng)度呈指數(shù)關(guān)系，統(tǒng)計分析大量的試驗數(shù)據(jù)得到

式中，σf為海冰彎曲強(qiáng)度，MPa。

由于不同海域的海冰生成環(huán)境不同，導(dǎo)致海冰的物理力學(xué)性質(zhì)也不盡相同。部分學(xué)者得到σf與呈線性關(guān)系[28]。本文對渤海沿岸試驗得到的145個試樣進(jìn)行統(tǒng)計分析，結(jié)果如圖6所示。從中可以發(fā)現(xiàn)，海冰彎曲強(qiáng)度與鹵水體積的平方根與彎曲強(qiáng)度呈指數(shù)關(guān)系，其可寫作

圖6 海冰鹵水體積對彎曲強(qiáng)度的影響Fig. 6 Effect of sea ice brine volume on bending strength

1.2.4 加載速率對彎曲強(qiáng)度的影響

加載速率對海冰彎曲強(qiáng)度的影響研究較少，但也引起部分學(xué)者的重視。Timco等[31]分析大量的海冰數(shù)據(jù)后發(fā)現(xiàn)海冰彎曲強(qiáng)度和加載速率并無明顯的函數(shù)關(guān)系，另有其他研究學(xué)者也對彎曲強(qiáng)度和加載速率的關(guān)系進(jìn)行了研究，Kermani等[32]發(fā)現(xiàn)淡水冰的彎曲強(qiáng)度和加載速率呈線性關(guān)系，隋吉學(xué)等[33]得到在不同溫度條件下，隨應(yīng)力率的增加海冰彎曲強(qiáng)度呈先減小后增大再急劇減小的趨勢。隋吉學(xué)等[13]的研究表明不同加載速率下的海冰彎曲強(qiáng)度呈現(xiàn)一定的韌脆轉(zhuǎn)化關(guān)系，對河冰的研究則表明隨應(yīng)力速率的增加彎曲強(qiáng)度呈馬鞍型分布[34]。而在單軸壓縮試驗中，李志軍等[35]研究團(tuán)隊得出海冰壓縮強(qiáng)度隨加載速率的增大先增加后減小，加載速率對壓縮強(qiáng)度的影響存在明顯的韌脆轉(zhuǎn)化關(guān)系。本文得到加載速率與海冰彎曲強(qiáng)度的關(guān)系如圖7所示，整體關(guān)系如圖7(a)所示。由于本文海冰取自不同地點，呈現(xiàn)的海冰冰晶結(jié)構(gòu)、冰溫和鹽度均不同，因而綜合所有數(shù)據(jù)的加載速率與海冰彎曲強(qiáng)度無明顯的函數(shù)關(guān)系。大量的試驗數(shù)據(jù)分析表明，鹵水體積平方根與彎曲強(qiáng)度呈現(xiàn)較為明顯的指數(shù)關(guān)系，為降低鹵水體積的影響，將鹵水體積平方根分為0～0.2、0.2～0.3、0.3～0.6三部分。圖7(b)中隨加載速率的演化海冰彎曲強(qiáng)度呈現(xiàn)先增大后減小的趨勢；而圖7(c)和(d)呈現(xiàn)線性關(guān)系，在圖7(d)中表現(xiàn)更為明顯。

圖7 海冰加載速率對彎曲強(qiáng)度的影響Fig. 7 Effect of sea ice loading rate on bending strength

1.2.5 密度及孔隙率對海冰彎曲強(qiáng)度的影響

海冰密度表征了海冰的鹽度、溫度以及氣泡等的綜合情況，反映了海冰孔隙率的大小。海冰中含有的鹵水、氣泡以及其他雜質(zhì)，隨著時間的推移，鹵水沿海冰微裂縫逐漸滲出，因此海冰的冰齡越大其密度越小[36]。本文應(yīng)用的海冰數(shù)據(jù)取自渤海中的當(dāng)年冰，其密度均值為898.9 kg/m3。單從海冰密度和彎曲強(qiáng)度的關(guān)系分析，如圖8所示，總體上海冰密度和彎曲強(qiáng)度并無一定的函數(shù)關(guān)系，在同一個海冰密度下海冰呈現(xiàn)不同的彎曲強(qiáng)度值，這是由于不同地點的海冰形成條件不同，導(dǎo)致的海冰冰晶結(jié)構(gòu)有差異。由于海冰密度對海冰彎曲強(qiáng)度影響的不確定性，因此對孔隙率和海冰彎曲強(qiáng)度的影響進(jìn)行分析，計算孔隙率的公式為[37-38]

圖8 海冰密度與彎曲強(qiáng)度的關(guān)系Fig. 8 Relationship between sea ice density and bending strength

式中，v為孔隙率，‰；va為氣泡體積比，‰；vb為鹵水體積比，‰；ρ為海冰密度，kg/m3；ρi為純冰密度，kg/m3；Si為海冰鹽度，‰；T為海冰溫度，℃。F1(T) 和F2(T) 為關(guān)于T的三次多項式，其系數(shù)列于表1。

表1 F1(T)和F2(T)多項式系數(shù)Table 1 Polynomial coefficients of F1(T) and F2(T)

本文得到海冰孔隙率與彎曲強(qiáng)度的關(guān)系如圖9所示，呈現(xiàn)指數(shù)關(guān)系，即

圖9 海冰孔隙率與彎曲強(qiáng)度的關(guān)系Fig. 9 Relationship between sea ice porosity and bending strength

式中，v為海冰內(nèi)孔隙率，‰。

以上研究均基于同一個加載方向?qū)ζ秸M(jìn)行的試驗，從以上各單因素分析可知，海冰彎曲強(qiáng)度與溫度、鹽度、密度及加載速率存在著復(fù)雜的非線性關(guān)系。

2 海冰彎曲強(qiáng)度預(yù)測的RNN模型

RNN模型在傳統(tǒng)的神經(jīng)網(wǎng)絡(luò)模型的基礎(chǔ)上發(fā)展起來，因而也具有很強(qiáng)的非線性數(shù)據(jù)處理能力。本文根據(jù)這一功能，并結(jié)合RNN在前一個輸入和后一個輸入有關(guān)系的基礎(chǔ)上，搭建了可應(yīng)用于海冰彎曲強(qiáng)度預(yù)測的模型。其中海冰彎曲強(qiáng)度預(yù)測的RNN的結(jié)構(gòu)圖如圖10所示，X表示海冰參數(shù)(溫度、鹽度、密度、加載速率)的輸入，O表示海冰彎曲強(qiáng)度的輸出，U表示輸入層到輸出層之間的鏈接權(quán)重，V表示隱藏層之間的鏈接權(quán)重，W表示相鄰隱藏層之間的鏈接權(quán)重，S表示隱藏層之間的值，圖中左側(cè)的結(jié)構(gòu)可以展開成右邊形式。每一個時間節(jié)點只與相鄰的時間節(jié)點進(jìn)行關(guān)聯(lián)，同時連接權(quán)重在不同時刻共享參數(shù)，即不同時刻的連接權(quán)重都相同，同時在模型更新的時候也一起更新。

圖10 海冰彎曲強(qiáng)度預(yù)測的RNN模型結(jié)構(gòu)圖Fig. 10 RNN model structure diagram for prediction of sea ice bending strength

海冰彎曲強(qiáng)度的RNN模型訓(xùn)練過程分為前向傳播和反向傳播，前向傳播過程為

式中，S表示隱藏層之間的值,f(·) 和g(·) 為激活函數(shù)，常用的激活函數(shù)有tanh，sigmoid，Relu等。

前向傳播完成后將產(chǎn)生海冰彎曲強(qiáng)度輸出值Ot，Ot和實際根據(jù)海冰參數(shù)及相應(yīng)公式計算得到的海冰彎曲強(qiáng)度值比較產(chǎn)生誤差et，利用反向傳播算法訓(xùn)練整個網(wǎng)絡(luò)，通過梯度下降法降低誤差完成訓(xùn)練。

為實現(xiàn)快速收斂，本文采用隨機(jī)梯度下降法進(jìn)行訓(xùn)練，對每一個參數(shù)都計算損失函數(shù)，然后求梯度更新參數(shù)。計算時將海冰數(shù)據(jù)的訓(xùn)練集分成多個小批量數(shù)據(jù)進(jìn)行訓(xùn)練，每個批次數(shù)據(jù)量表示為batchsize。與全數(shù)據(jù)集訓(xùn)練相比，在尋找最優(yōu)時增加了一些隨機(jī)噪聲，用以修正由局部數(shù)據(jù)得到的梯度，盡量避免因batchsize過大陷入局部最優(yōu)。其他參數(shù)一定時，batchsize越小，一個批次選擇數(shù)據(jù)的隨機(jī)性越大，訓(xùn)練結(jié)果越不易收斂。但batchsize越小計算速度越快，訓(xùn)練得到的權(quán)值更新越頻繁。由于采用隨機(jī)選擇batch的方法，對于非凸損失函數(shù)更便于尋找全局最優(yōu)，因而選擇隨機(jī)梯度下降法收斂更快更容易達(dá)到全局最優(yōu)。

訓(xùn)練過程中激勵函數(shù)分別采用tanh，sigmoid，Relu進(jìn)行訓(xùn)練，發(fā)現(xiàn)用sigmoid函數(shù)進(jìn)行訓(xùn)練的結(jié)果最好，這是因為sigmoid函數(shù)將結(jié)果映射為(0, 1)之間，和歸一化后的數(shù)據(jù)結(jié)果有較強(qiáng)的相關(guān)性。另外，在反向傳播過程中，將均方根誤差(RMSE)作為損失函數(shù)。在隨機(jī)梯度下降法中選擇Adam優(yōu)化算法，此算法僅需較小的內(nèi)存即可完成計算，而且使用動量和自適應(yīng)學(xué)習(xí)率可加快收斂速度，能夠糾正學(xué)習(xí)率消失或高方差參數(shù)更新導(dǎo)致?lián)p失函數(shù)波動較大的問題。

3 海冰彎曲強(qiáng)度的循環(huán)神經(jīng)網(wǎng)絡(luò)預(yù)測

3.1 海冰彎曲強(qiáng)度預(yù)測的RNN模型訓(xùn)練結(jié)果

本文選擇影響海冰彎曲強(qiáng)度的溫度、鹽度、密度以及加載速率作為輸入?yún)?shù)，與之對應(yīng)的海冰彎曲強(qiáng)度值作為輸出。為了使訓(xùn)練的收斂效果更好，首先采用隨機(jī)分配數(shù)據(jù)的方法將數(shù)據(jù)的順序打亂，并采用歸一化方法，將數(shù)據(jù)變?yōu)?0, 1)之間的小數(shù)，其歸一化公式為

本文共得到了145組海冰參數(shù)及對應(yīng)的彎曲強(qiáng)度數(shù)據(jù)，在訓(xùn)練時按照8∶1∶1的比例將數(shù)據(jù)分為訓(xùn)練集，驗證集和測試集。訓(xùn)練集適用于模型學(xué)習(xí)的樣本，驗證集為訓(xùn)練過程中留出的樣本，可以為模型調(diào)參和模型評估提供參考，而測試集不參與任何的訓(xùn)練和測試，用來測試樣本外的數(shù)據(jù)模型評估。由于數(shù)據(jù)集的限制，本文選擇batchsize = 2訓(xùn)練得到的結(jié)果最好。RNN模型通過不斷訓(xùn)練調(diào)整權(quán)重和閾值，訓(xùn)練次數(shù)過少會導(dǎo)致精度不夠，而訓(xùn)練次數(shù)過多將導(dǎo)致過擬合，使得訓(xùn)練不能收斂。隱藏層層數(shù)和節(jié)點數(shù)也要隨著訓(xùn)練進(jìn)行調(diào)整，通過多次測試，最終調(diào)整隱藏層層數(shù)為2，每一個隱藏層的節(jié)點數(shù)目為30，將不同參數(shù)的訓(xùn)練結(jié)果列于表2。

表2 海冰彎曲強(qiáng)度值預(yù)測的RNN模型不同參數(shù)值的訓(xùn)練結(jié)果Table 2 Training results of different parameter values of RNN model for prediction of sea ice bending strength

從表中看出，訓(xùn)練次數(shù)為500，學(xué)習(xí)率選取0.001時得到的均方根誤差RMSE值最小，其訓(xùn)練結(jié)果如圖11所示。訓(xùn)練評估指標(biāo)包括均方根誤差、擬合優(yōu)度、平均絕對值誤差等，可對參數(shù)的調(diào)整提供指導(dǎo)。隨訓(xùn)練次數(shù)的增加，訓(xùn)練集、驗證集和測試集的均方根誤差均較為一致趨于收斂，最優(yōu)訓(xùn)練結(jié)果的訓(xùn)練集均方根誤差為0.13。

圖11 海冰彎曲強(qiáng)度值預(yù)測的RNN模型訓(xùn)練結(jié)果Fig. 11 Training results of RNN model for prediction of bending strength of sea ice

3.2 彎曲強(qiáng)度的預(yù)測結(jié)果

在海冰物理力學(xué)性質(zhì)的研究中，鹵水體積（溫度和鹽度）與彎曲強(qiáng)度的函數(shù)關(guān)系最為明顯。因此，本文利用鹵水體積的平方根和海冰彎曲強(qiáng)度的關(guān)系圖（圖6），擬合得到公式（式5），通過鹵水體積完成對海冰彎曲強(qiáng)度的計算。本節(jié)利用RNN模型對本文得到的145組數(shù)據(jù)進(jìn)行訓(xùn)練，將海冰溫度、鹽度、密度和加載速率作為參數(shù)輸入，海冰彎曲強(qiáng)度作為輸出，對測試集以及所有數(shù)據(jù)進(jìn)行預(yù)測，并采用經(jīng)驗公式計算的結(jié)果和預(yù)測結(jié)果進(jìn)行對比。

圖12為測試集14個數(shù)據(jù)的RNN預(yù)測值、經(jīng)驗公式計算值與實測計算得到的彎曲強(qiáng)度值對比結(jié)果。從圖中可以看出，RNN預(yù)測的結(jié)果和基準(zhǔn)線的擬合度更好，利用RNN預(yù)測的結(jié)果要明顯好于經(jīng)驗公式預(yù)測結(jié)果。與實際測量值對比，RNN預(yù)測結(jié)果的平均誤差為12.01%，均方根誤差為11.46%，有7個樣本的預(yù)測誤差在10%以下，5個樣本的誤差在10%～20%，預(yù)測樣本的最大誤差為31.49%。 而經(jīng)驗公式預(yù)測結(jié)果的平均誤差為37.95%，均方根誤差為35.68%，最大誤差為148.60%。僅有2個樣本的誤差處于10%以下，誤差大于50%的多達(dá)8個，將誤差對比結(jié)果列于表3，可以明顯看出RNN預(yù)測結(jié)果遠(yuǎn)比經(jīng)驗公式預(yù)測結(jié)果好。

圖12 不同方法對測試集數(shù)據(jù)預(yù)測的海冰彎曲強(qiáng)度值與實際測量計算值的對比結(jié)果Fig. 12 The results of the comparison between the sea ice bending strength predicted by different methods and the actual measured values are obtained

表3 不同方法對測試集數(shù)據(jù)預(yù)測的海冰彎曲強(qiáng)度準(zhǔn)確度Table 3 Accuracy of sea ice bending strength predicted by different methods on test set data

對于所有數(shù)據(jù)的彎曲強(qiáng)度對比結(jié)果，如圖13所示，與經(jīng)驗公式預(yù)測的海冰彎曲強(qiáng)度結(jié)果對比，RNN預(yù)測結(jié)果較為集中地分布在基準(zhǔn)線兩側(cè)，RNN預(yù)測值和基準(zhǔn)線擬合效果更好。

圖13 不同方法對所有數(shù)據(jù)預(yù)測的海冰彎曲強(qiáng)度值與實際計算值的對比結(jié)果Fig. 13 The comparison results between the predicted and calculated values of sea ice bending strength by different methods

為了使對比數(shù)據(jù)更加清晰，繪制了圖14所示不同方法對所有數(shù)據(jù)的鹵水體積平方根和彎曲強(qiáng)度的預(yù)測結(jié)果，并和實際計算的彎曲強(qiáng)度值對比，誤差分析結(jié)果列于表4中。RNN預(yù)測的結(jié)果平均誤差為16.47%，而經(jīng)驗公式的平均誤差值超過RNN預(yù)測結(jié)果的兩倍，為37.32%；均方根誤差中RNN的11.82%也明顯優(yōu)于經(jīng)驗公式的32.91%；誤差值分別在10%以下和20%以下預(yù)測結(jié)果中，RNN的預(yù)測結(jié)果也明顯優(yōu)于經(jīng)驗公式，表明RNN在預(yù)測驗證集之外的數(shù)據(jù)時也有較好的準(zhǔn)確性，和真實結(jié)果更為接近，體現(xiàn)出RNN模型具有較強(qiáng)的泛化能力。

圖14 不同方法對鹵水體積平方根和彎曲強(qiáng)度的關(guān)系預(yù)測結(jié)果Fig. 14 Prediction results of the relationship between square root of brine volume and bending strength by different methods

表4 不同方法對所有數(shù)據(jù)預(yù)測的海冰彎曲強(qiáng)度準(zhǔn)確度Table 4 Accuracy of sea ice bending strength predicted by different methods for all data

4 結(jié)論

本文介紹了海冰三點彎曲試驗的測量方法并對渤海海冰彎曲強(qiáng)度的試驗結(jié)果進(jìn)行了總結(jié)，分析了海冰溫度、鹽度、加載速率和密度等因素影響下海冰彎曲強(qiáng)度的演變規(guī)律，得到各因素和海冰彎曲強(qiáng)度間復(fù)雜的非線性關(guān)系。由于單個影響因素對彎曲強(qiáng)度的影響呈現(xiàn)不同的非線性關(guān)系的特點，不能綜合各個因素完成彎曲強(qiáng)度的精確預(yù)測。由此本文構(gòu)建了泛化能力較強(qiáng)的RNN預(yù)測模型。通過采用隨機(jī)梯度下降法優(yōu)化模型，對數(shù)據(jù)歸一化處理后進(jìn)行訓(xùn)練。調(diào)參優(yōu)化后完成了海冰溫度、鹽度、密度和加載速率影響下的海冰彎曲強(qiáng)度的預(yù)測。RNN預(yù)測的海冰彎曲強(qiáng)度同經(jīng)驗公式對比表明，該RNN模型對海冰彎曲強(qiáng)度有很好的準(zhǔn)確性。