鄧林青， 朱耀文， 王宏偉， 張勇青， 李彤濱

(1. 中海石油深海開(kāi)發(fā)有限公司，深圳 518067；2. 哈爾濱工程大學(xué) 船舶工程學(xué)院，哈爾濱 150001；3. 海洋石油工程股份有限公司，天津 300451)

0 引 言

浮式生產(chǎn)儲(chǔ)卸油平臺(tái)FPSO(floating production storage and offloading)在海洋石油勘探開(kāi)發(fā)中被廣泛應(yīng)用，一般采用系泊系統(tǒng)限制其在海面的運(yùn)動(dòng)。由于人類對(duì)資源需求的日益增長(zhǎng)，海洋油氣的開(kāi)發(fā)不斷走向深海，從而使浮式平臺(tái)及其附屬系統(tǒng)面臨更為惡劣的海況[1]。為了保證各類平臺(tái)的安全作業(yè)，浮式結(jié)構(gòu)的運(yùn)動(dòng)響應(yīng)范圍以及系泊纜的張力極值必須在規(guī)定的范圍之內(nèi)。因此，在設(shè)計(jì)初期，需要借助復(fù)雜的時(shí)域耦合分析軟件，對(duì)復(fù)雜的外部環(huán)境荷載進(jìn)行相應(yīng)的計(jì)算與分析。在外部風(fēng)、浪、流等載荷的聯(lián)合作用下，浮式結(jié)構(gòu)運(yùn)動(dòng)與系泊及立管系統(tǒng)之間的耦合作用不能忽視，深水細(xì)長(zhǎng)結(jié)構(gòu)物顯著的動(dòng)態(tài)特性也需要準(zhǔn)確地模擬，這些因素將導(dǎo)致整個(gè)系統(tǒng)在時(shí)域仿真求解時(shí)需要耗費(fèi)較多的計(jì)算時(shí)間。

當(dāng)前已有大量的研究投入到尋求一個(gè)可靠的、耗時(shí)較少的浮式結(jié)構(gòu)系泊系統(tǒng)耦合分析方法[2-4]。人工神經(jīng)網(wǎng)絡(luò)是一種基于數(shù)據(jù)驅(qū)動(dòng)的方法，在輸入與輸出之間起到了類似傳遞函數(shù)的作用，適用于尋找隨機(jī)輸入與相應(yīng)輸出之間的映射關(guān)系，其在船舶與海洋工程領(lǐng)域已經(jīng)得到了廣泛的應(yīng)用，包括浮式結(jié)構(gòu)及其系泊系統(tǒng)設(shè)計(jì)參數(shù)的優(yōu)化，波浪有義波高、環(huán)境荷載、船舶運(yùn)動(dòng)的預(yù)報(bào)，纜線與結(jié)構(gòu)的疲勞分析等[5-8]。Guarize等人[9]提出了一種基于人工神經(jīng)網(wǎng)絡(luò)與有限元模型的混合方法，用于細(xì)長(zhǎng)海洋結(jié)構(gòu)物的時(shí)域動(dòng)態(tài)分析，前期使用仿真軟件得出短時(shí)間內(nèi)的結(jié)果，作為神經(jīng)網(wǎng)絡(luò)的訓(xùn)練數(shù)據(jù)，將已知的浮體運(yùn)動(dòng)輸入訓(xùn)練好的神經(jīng)網(wǎng)絡(luò)中進(jìn)行張力預(yù)報(bào)，得到一個(gè)長(zhǎng)期的時(shí)域響應(yīng)。孫麗萍等人[10]將BP神經(jīng)網(wǎng)絡(luò)和遺傳算法相結(jié)合，對(duì)FPSO系泊系統(tǒng)的長(zhǎng)度進(jìn)行優(yōu)化，縮短了優(yōu)化所需的時(shí)間。張隆輝等人[11-12]在構(gòu)建系泊系統(tǒng)主動(dòng)式截?cái)喾抡孢^(guò)程中，使用NARX神經(jīng)網(wǎng)絡(luò)模型代替了傳統(tǒng)纜索數(shù)值模型，用于預(yù)報(bào)系泊線截?cái)帱c(diǎn)的運(yùn)動(dòng)響應(yīng)。該模型能夠加快計(jì)算求解的速度，有助于主動(dòng)式混合模型試驗(yàn)的實(shí)現(xiàn)。Pina[13-15]使用NARX神經(jīng)網(wǎng)絡(luò)進(jìn)行了系泊線和立管的力學(xué)分析。以浮體的運(yùn)動(dòng)為輸入?yún)?shù)，系泊和立管的頂端張力為輸出參數(shù)，通過(guò)對(duì)比神經(jīng)網(wǎng)絡(luò)模型預(yù)報(bào)結(jié)果與有限元計(jì)算結(jié)果，驗(yàn)證了該模型能夠準(zhǔn)確地預(yù)報(bào)頂端張力，并具有較高的計(jì)算效率。隨后，將小波函數(shù)作為激活函數(shù)，使用小波神經(jīng)網(wǎng)絡(luò)代替?zhèn)鹘y(tǒng)的有限元仿真，進(jìn)行海洋細(xì)長(zhǎng)結(jié)構(gòu)物的分析。上述方法均是基于預(yù)先規(guī)定的浮體運(yùn)動(dòng)，沒(méi)有考慮浮體與系泊系統(tǒng)之間的耦合作用。為此，本研究將波浪的時(shí)歷升高作為神經(jīng)網(wǎng)絡(luò)模型的輸入?yún)?shù)，使用神經(jīng)網(wǎng)絡(luò)模型進(jìn)行浮式生產(chǎn)系統(tǒng)的時(shí)域耦合分析。

本文將NARX神經(jīng)網(wǎng)絡(luò)應(yīng)用于某一分布式FPSO系泊系統(tǒng)時(shí)域耦合分析中，該模型包括兩種訓(xùn)練模式，序列并行模式(NARX-SP)與并行模式(NARX-P)。前者需要將當(dāng)前時(shí)間步之前的真實(shí)輸出目標(biāo)數(shù)據(jù)進(jìn)行反饋，而后者需要將預(yù)報(bào)出的數(shù)據(jù)反饋回輸入層，由于輸出的真實(shí)目標(biāo)數(shù)據(jù)無(wú)法提前得知，因此選取并行模式(NARX-P)進(jìn)行訓(xùn)練與預(yù)報(bào)。首先使用OrcaFlex軟件計(jì)算得到短時(shí)的時(shí)域耦合結(jié)果，以已知波浪的時(shí)歷升高數(shù)據(jù)作為輸入，系泊纜張力響應(yīng)為目標(biāo)輸出，使該神經(jīng)網(wǎng)絡(luò)具有預(yù)報(bào)余下時(shí)間歷程張力響應(yīng)的能力。經(jīng)比較發(fā)現(xiàn)，預(yù)報(bào)結(jié)果與時(shí)域耦合分析計(jì)算結(jié)果之間的最大相對(duì)誤差為7.61%，而且該替代模型能夠大幅減小數(shù)值仿真的時(shí)間。

1 NARX神經(jīng)網(wǎng)絡(luò)與實(shí)施方案

NARX神經(jīng)網(wǎng)絡(luò)是一種動(dòng)態(tài)的循環(huán)神經(jīng)網(wǎng)絡(luò)，以誤差反向傳播的多層前饋神經(jīng)網(wǎng)絡(luò)(BPNN)為基礎(chǔ)，引入了適當(dāng)?shù)难訒r(shí)參數(shù)與反饋?zhàn)饔?。此神?jīng)網(wǎng)絡(luò)的輸出不僅取決于當(dāng)前時(shí)間步的輸入，還受到過(guò)去幾個(gè)時(shí)間步神經(jīng)網(wǎng)絡(luò)輸入、輸出值的影響，在一定程度上能夠緩解時(shí)域數(shù)據(jù)間存在的長(zhǎng)期依賴問(wèn)題，改善神經(jīng)網(wǎng)絡(luò)的性能，適用于時(shí)域非線性動(dòng)態(tài)系統(tǒng)的預(yù)報(bào)，其作為一種強(qiáng)大的時(shí)域預(yù)報(bào)模型已被廣泛地應(yīng)用[16-17]。

1.1 神經(jīng)元與激活函數(shù)

圖1 McCulloch-Pitts模型Fig.1 McCulloch-Pitts model

神經(jīng)元是神經(jīng)網(wǎng)絡(luò)的基礎(chǔ)，也是神經(jīng)網(wǎng)絡(luò)的基本單元，其有興奮和抑制兩種狀態(tài)。在一般情況下，大多數(shù)的神經(jīng)元處于抑制狀態(tài)，但是一旦某個(gè)神經(jīng)元受到刺激，將導(dǎo)致它超過(guò)某個(gè)閾值，這個(gè)神經(jīng)元就會(huì)被激活，處于“興奮”的狀態(tài)，進(jìn)而向其他的神經(jīng)元傳遞信息。1943年，心理學(xué)家McCulloch和數(shù)學(xué)家Pitts參考了生物神經(jīng)元的結(jié)構(gòu)，發(fā)表了抽象的神經(jīng)元模型——MP模型，如圖1所示。MP模型接受已知特征數(shù)據(jù)的輸入xi，i=1，N，每個(gè)神經(jīng)元使用權(quán)重wi、偏置b與輸入數(shù)據(jù)進(jìn)行線性組合后得到權(quán)值輸出z，為了引入非線性特性，加權(quán)輸出還需經(jīng)過(guò)激活函數(shù)f(z)的轉(zhuǎn)換得到輸出y。MP模型的建立為神經(jīng)網(wǎng)絡(luò)的發(fā)展奠定了基礎(chǔ)，但是由于其預(yù)先定義的權(quán)值無(wú)法進(jìn)行學(xué)習(xí)，后來(lái)逐步提出了單層感知器、多層感知器以及誤差的反向傳播算法。當(dāng)前各種形式的神經(jīng)網(wǎng)絡(luò)模型被不斷地提出，廣泛應(yīng)用于各個(gè)領(lǐng)域。

神經(jīng)網(wǎng)絡(luò)已有多種激活函數(shù)可供選擇，例如常用于分類問(wèn)題的邏輯函數(shù)，以及用于回歸問(wèn)題的線性整流函數(shù)。本文使用的是如圖2所示的雙曲正切激活函數(shù)，它是一個(gè)值域具有上下界的連續(xù)函數(shù)：

(1)

式中：y為輸出結(jié)果；z為加權(quán)輸出。

圖2 雙曲正切激活函數(shù)Fig.2 Hyperbolic tangent activation function

1.2 NARX神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)

基于圖3所示的雙層BP神經(jīng)網(wǎng)絡(luò)(由于輸入層未進(jìn)行任何變換，沒(méi)有作為單獨(dú)的一層考慮)只包含了一個(gè)隱藏層，圖4展示了NARX神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)。當(dāng)前的神經(jīng)網(wǎng)絡(luò)輸出不僅取決于當(dāng)前和過(guò)去的輸入，還取決于過(guò)去的輸出，其輸入集和最終輸出結(jié)果可以分別表示為

X(t)=(x(t),x(t-Δt), …,x(t-pΔt),
y(t-Δt), …,y(t-qΔt))

(2)

yt=f(w，X(t))

(3)

式中：x(t)為神經(jīng)網(wǎng)絡(luò)的輸入；yt為輸出結(jié)果；w為權(quán)重參數(shù)；f(·)代表了NARX神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)建立的映射關(guān)系模型；Δt為時(shí)域耦合分析所用的時(shí)間步長(zhǎng)；p和q分別為輸入與輸出反饋的延遲階數(shù)；pΔt和qΔt為延時(shí)的時(shí)間長(zhǎng)度。輸入與輸出的延遲效應(yīng)使得NARX具有了一定的記憶能力。

圖3 雙層BPNN結(jié)構(gòu)Fig.3 Double-layer BPNN structure

圖4 NARX神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)Fig.4 NARX neural network structure

神經(jīng)網(wǎng)絡(luò)模型輸入?yún)?shù)的選擇對(duì)于最終的結(jié)果有很大的影響。在浮式結(jié)構(gòu)系泊系統(tǒng)時(shí)域耦合分析過(guò)程中，浮式結(jié)構(gòu)水動(dòng)力性能(AQWA軟件預(yù)先計(jì)算)與外部的環(huán)境荷載是已知的條件，但這些都是不變的常數(shù)項(xiàng)，為了能夠預(yù)報(bào)時(shí)域的輸出，需要有相應(yīng)的時(shí)域輸入數(shù)據(jù)。文中選用波高的時(shí)歷曲線作為輸入，F(xiàn)PSO系泊系統(tǒng)的系泊線頂部張力值作為輸出。

需要注意的是，NARX神經(jīng)網(wǎng)絡(luò)具有兩種訓(xùn)練模式。第一種模式是序列并行模式(NARS-SP)，這種方式的輸入層由實(shí)際的目標(biāo)值d(t)組成，即反饋回去的是實(shí)際的軟件仿真結(jié)果，可表示為

X(t)=(x(t),x(t-Δt), …,
x(t-pΔt),d(t-Δt), …,d(t-qΔt))

(4)

另一種模式為并行模式(NARX-P)，不同之處在于其將預(yù)報(bào)的輸出反饋到輸入層，如式(2)所示。在實(shí)際進(jìn)行預(yù)報(bào)時(shí)，只有在訓(xùn)練集與驗(yàn)證集中的d(t)是已知的，可以通過(guò)OrcaFlex軟件進(jìn)行短時(shí)的時(shí)域耦合分析得到，而測(cè)試集的結(jié)果則使用訓(xùn)練完成的神經(jīng)網(wǎng)絡(luò)模型進(jìn)行預(yù)報(bào)。

本文使用神經(jīng)網(wǎng)絡(luò)模型預(yù)報(bào)系泊線的頂端張力響應(yīng)，以提高計(jì)算效率，減小時(shí)域耦合計(jì)算時(shí)間。由于無(wú)法得知當(dāng)前時(shí)刻實(shí)際的張力值，因此使用NARX-P模式進(jìn)行訓(xùn)練和預(yù)報(bào)。后續(xù)提到的NARX神經(jīng)網(wǎng)絡(luò)均表示使用NARP-P模式進(jìn)行訓(xùn)練與預(yù)報(bào)的模型，該模型使用MATLAB提供的神經(jīng)網(wǎng)絡(luò)工具箱來(lái)完成[18]。

1.3 應(yīng)用方案

如NARX此類的監(jiān)督神經(jīng)網(wǎng)絡(luò)均需要一定數(shù)量的訓(xùn)練集進(jìn)行訓(xùn)練，再通過(guò)驗(yàn)證集對(duì)模型的參數(shù)等進(jìn)行適當(dāng)?shù)恼{(diào)整，最終運(yùn)用到測(cè)試集上去預(yù)報(bào)未知的數(shù)據(jù)結(jié)果，這需要一個(gè)短時(shí)的時(shí)域耦合計(jì)算結(jié)果。同時(shí)，神經(jīng)網(wǎng)絡(luò)模型對(duì)訓(xùn)練集與驗(yàn)證集上的輸入與輸出數(shù)據(jù)必須是已知的，由于FPSO及其系泊系統(tǒng)在時(shí)域耦合分析過(guò)程中無(wú)法預(yù)先得知船舶運(yùn)動(dòng)時(shí)歷數(shù)據(jù)，即無(wú)法使用船舶運(yùn)動(dòng)作為神經(jīng)網(wǎng)絡(luò)的輸入?yún)?shù)。本文選擇波高時(shí)歷數(shù)據(jù)作為輸入，系泊線頂端張力作為輸出。

浮式結(jié)構(gòu)系泊系統(tǒng)時(shí)域耦合運(yùn)動(dòng)模擬時(shí)間設(shè)置為3 h(10 800 s)，計(jì)算步長(zhǎng)為0.1 s。首先使用OrcaFlex軟件計(jì)算800 s內(nèi)的浮體系泊系統(tǒng)時(shí)域結(jié)果，其中600 s用于訓(xùn)練，200 s用于驗(yàn)證，每一個(gè)訓(xùn)練周期后，該模型都會(huì)預(yù)報(bào)驗(yàn)證集上的系泊張力，并與實(shí)際3 h數(shù)值模擬計(jì)算結(jié)果進(jìn)行對(duì)比。如果在6個(gè)周期中驗(yàn)證集上的誤差均呈現(xiàn)增長(zhǎng)，則表明該模型無(wú)法繼續(xù)優(yōu)化，即完成訓(xùn)練，余下時(shí)間的張力響應(yīng)可由建立的NARX神經(jīng)網(wǎng)絡(luò)模型預(yù)報(bào)得出。為了驗(yàn)證最終在測(cè)試集上的預(yù)報(bào)結(jié)果，文中進(jìn)行了3 h的時(shí)域計(jì)算，實(shí)際使用時(shí)只需要計(jì)算短時(shí)目標(biāo)數(shù)據(jù)(800 s)即可。實(shí)施方案如圖5所示，其中的誤差指的是均方誤差：

(5)

式中：n為輸入數(shù)據(jù)的數(shù)量；y(t)和d(t)分別為預(yù)報(bào)值和實(shí)際目標(biāo)值。

圖5 NARX神經(jīng)網(wǎng)絡(luò)應(yīng)用方案Fig.5 The application scheme of NARX neural network

2 數(shù)值計(jì)算分析

2.1 FPSO系泊系統(tǒng)模型與環(huán)境參數(shù)

上述方法被應(yīng)用于某一分布式FPSO系泊系統(tǒng)，作業(yè)水深為1 200 m，平臺(tái)與系泊纜的具體參數(shù)如表1和表2所示，共有12根系泊纜，分為4組，每組3根。每根系泊纜均鏈條-鋼纜-鏈條三段材料構(gòu)成，系泊編號(hào)及模型如圖6所示。系統(tǒng)受到210°方向的風(fēng)、浪、流的聯(lián)合作用，表面流速為1.53 m/s，平均風(fēng)速為25.6 m/s，波浪有義波高與峰值周期分別為7.25 m和19.78 s。在整個(gè)仿真過(guò)程中，流速與風(fēng)速的大小和方向都是定常的，而波高是隨時(shí)間變化的。為了預(yù)報(bào)系泊纜張力的時(shí)歷輸出，使用波高時(shí)歷數(shù)據(jù)作為輸入。波浪采用JONSWAP譜，可定義為

(6)

式中:α為能量尺度參數(shù)；γ為譜峰升高因子；σ為峰形參數(shù)；ωm為譜峰頻率。

表1 FPSO平臺(tái)參數(shù)

表2 系泊線參數(shù)

波浪譜確定后，即可以得到相對(duì)于空間某點(diǎn)的波面升高時(shí)歷曲線。在OrcaFlex軟件中可以得到整個(gè)時(shí)域過(guò)程的波高時(shí)域數(shù)據(jù)，這些數(shù)據(jù)在計(jì)算之前即可得知，這是前提條件。在實(shí)際應(yīng)用時(shí)會(huì)發(fā)現(xiàn)，輸入多點(diǎn)的波面升高更有利于NARX神經(jīng)網(wǎng)絡(luò)發(fā)現(xiàn)數(shù)據(jù)間的映射關(guān)系。因此，基于FPSO系泊系統(tǒng)靜態(tài)平衡位置(-48.5, -55.6)處的波面升高時(shí)歷曲線，再選擇(-28.5, -35.6)、 (-38.5, -45.6)、 (-58.5, -65.6)、 (-68.5, -75.6)，共五個(gè)點(diǎn)的波浪升高作為輸入?yún)?shù)。圖7給出了(-48.5, -55.6)位置處的波面升高時(shí)歷曲線，上述靜態(tài)平衡位置的計(jì)算也可以避免時(shí)域耦合分析時(shí)浮式平臺(tái)的瞬態(tài)效應(yīng)，防止不合理數(shù)據(jù)的生成。

圖6 FPSO系泊系統(tǒng)模型

圖7 波浪升高時(shí)歷曲線Fig.7 Time history of wave elevations

2.2 超參數(shù)設(shè)定

首先使用OrcaFlex計(jì)算800s內(nèi)的時(shí)域耦合數(shù)據(jù)，作為NARX神經(jīng)網(wǎng)絡(luò)的訓(xùn)練集和驗(yàn)證集。隨后，需要確定該模型的輸入與輸出延遲數(shù)，以及隱藏層的神經(jīng)元數(shù)量，這些超參數(shù)對(duì)神經(jīng)網(wǎng)絡(luò)模型的性能有很大影響。采用試湊法對(duì)隱藏層數(shù)和延遲數(shù)分別為5～20和5～30，間隔為5的區(qū)間范圍內(nèi)進(jìn)行搜索，尋找結(jié)果較優(yōu)的超參數(shù)，輸入與輸出使用相同的延遲數(shù)。選用系泊纜Line3，由于系泊纜張力量級(jí)較大，為了清晰地表示真實(shí)值與預(yù)報(bào)值之間的誤差，表3展示的為均方根的誤差。此處只給出了部分結(jié)果，當(dāng)隱藏層神經(jīng)元數(shù)為20時(shí)，輸入與輸出的延遲數(shù)為15時(shí)，均方根誤差最小，為87.92 kN。在此過(guò)程中，NARX神經(jīng)網(wǎng)絡(luò)的訓(xùn)練也已經(jīng)完成，當(dāng)相應(yīng)超參數(shù)確定后，即可用此參數(shù)下的NARX模型對(duì)測(cè)試集上的系泊纜張力進(jìn)行預(yù)報(bào)。

一般神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)能力能夠隨隱藏層數(shù)與延遲數(shù)的增大而加強(qiáng)，但是太過(guò)復(fù)雜的神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)會(huì)導(dǎo)致在訓(xùn)練過(guò)程中的過(guò)擬合問(wèn)題，使其在驗(yàn)證集上的性能反而不好。網(wǎng)絡(luò)層數(shù)太少又會(huì)導(dǎo)致訓(xùn)練時(shí)的欠擬合，使得該神經(jīng)網(wǎng)絡(luò)無(wú)法學(xué)習(xí)到輸入和輸出之間的關(guān)系。因此，對(duì)于超參數(shù)的正確選擇是神經(jīng)網(wǎng)絡(luò)模型能否有效的關(guān)鍵因素之一，此處采用試湊法確定了一個(gè)較優(yōu)的神經(jīng)網(wǎng)絡(luò)模型結(jié)構(gòu)。

2.3 結(jié)果與分析

使用訓(xùn)練好的神經(jīng)網(wǎng)絡(luò)模型去預(yù)報(bào)余下10 000 s系泊纜Line3的時(shí)域張力結(jié)果，在此基礎(chǔ)上，可以使用相同的神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)進(jìn)行系泊纜Line4的訓(xùn)練和預(yù)報(bào)。分別使用OrcaFlex軟件和神經(jīng)網(wǎng)絡(luò)模型進(jìn)行計(jì)算，圖8～圖11分別展示了2 000～2 400 s和10 400～10 800 s的Line3和Line4張力時(shí)歷曲線對(duì)比結(jié)果?？梢?jiàn)，預(yù)報(bào)結(jié)果與目標(biāo)值(OrcaFlex軟件計(jì)算結(jié)果)之間的時(shí)歷曲線趨勢(shì)幾乎完全一致，結(jié)果吻合良好，但在每個(gè)時(shí)間段的峰值上略有差異。兩種計(jì)算方式在時(shí)域結(jié)果上的相對(duì)誤差如表4和表5所示，其中系泊纜頂端張力平均值的相對(duì)誤差較小，系泊纜Line3和Line4張力最小值的相對(duì)誤差分別為7.10%和7.61%，張力最大值的相對(duì)誤差分別為5.33%和5.81%。整體上，相對(duì)誤差的最大值為7.61%，在可接受的范圍內(nèi)。

圖8 Line3頂部張力時(shí)域曲線2 000～2 400 sFig.8 Line3 top tension time series 2 000～2 400 s

圖9 Line3頂部張力時(shí)域曲線10 400～10 800 sFig.9 Line3 top tension time series 10 400～10 800 s

圖10 Line4頂部張力時(shí)域曲線2 000～2 400 sFig.10 Line4 top tension time series 2 000～2 400 s

圖11 Line4頂部張力時(shí)域曲線10 400～10 800 sFig.11 Line4 top tension time series 10 400～10 800 s

表4 Line3頂部張力時(shí)域數(shù)據(jù)分析

表5 Line4頂部張力時(shí)域數(shù)據(jù)分析

此種方法只需計(jì)算短時(shí)(800 s)時(shí)域耦合分析數(shù)據(jù)，使用神經(jīng)網(wǎng)絡(luò)模型預(yù)報(bào)其后10 000 s的系泊纜頂部張力，可以大幅減小OrcaFlex軟件進(jìn)行浮式結(jié)構(gòu)系泊系統(tǒng)時(shí)域耦合分析的時(shí)間，計(jì)算效率大約為常規(guī)方法的13倍。這里的計(jì)算效率沒(méi)有考慮神經(jīng)網(wǎng)絡(luò)構(gòu)建以及訓(xùn)練、預(yù)報(bào)的時(shí)間，其中神經(jīng)網(wǎng)絡(luò)的構(gòu)建，包括文中提到的超參數(shù)的選擇會(huì)耗費(fèi)相對(duì)較多的時(shí)間。但是，對(duì)于同一類的問(wèn)題，使用相同或相近的網(wǎng)絡(luò)結(jié)構(gòu)，也可以得到較好的結(jié)果。本文中的神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)與相關(guān)文獻(xiàn)所用結(jié)構(gòu)十分接近[13]，使用了具有20個(gè)神經(jīng)元的單隱藏層、延遲數(shù)為15的神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)。

3 結(jié) 語(yǔ)

以某一分布式FPSO系泊系統(tǒng)為例，基于NARX神經(jīng)網(wǎng)絡(luò)模型進(jìn)行FPSO系泊系統(tǒng)耦合分析中系泊纜頂部張力預(yù)報(bào)，對(duì)比OrcaFlex軟件計(jì)算結(jié)果與NARX神經(jīng)網(wǎng)絡(luò)模型預(yù)報(bào)結(jié)果，可知兩者的最大相對(duì)誤差為7.61%。此方法是將OrcaFlex軟件仿真計(jì)算與NARX神經(jīng)網(wǎng)絡(luò)模型相結(jié)合的一種方法，可以大幅減小軟件時(shí)域耦合分析的時(shí)間，計(jì)算效率大約為常規(guī)方法的13倍。但是，文中使用的NARX神經(jīng)網(wǎng)絡(luò)模型只適用于上述浮式平臺(tái)系泊系統(tǒng)的張力預(yù)報(bào)，如果需要更換浮式平臺(tái)以及系泊系統(tǒng)的參數(shù)，則需要重新對(duì)該神經(jīng)網(wǎng)絡(luò)模型進(jìn)行訓(xùn)練。