基于改進(jìn)自適應(yīng)遺傳算法的海洋平臺(tái)管路優(yōu)化設(shè)計(jì)方法＊

伊 鵬 劉衍聰 石永軍 任紅偉 鄭煥鵬

(中國(guó)石油大學(xué)(華東)機(jī)電工程學(xué)院)

基于改進(jìn)自適應(yīng)遺傳算法的海洋平臺(tái)管路優(yōu)化設(shè)計(jì)方法＊

伊 鵬 劉衍聰 石永軍 任紅偉 鄭煥鵬

(中國(guó)石油大學(xué)(華東)機(jī)電工程學(xué)院)

1 系統(tǒng)模型

海洋平臺(tái)的設(shè)計(jì)進(jìn)行到管道設(shè)計(jì)階段,平臺(tái)上的主要設(shè)備及構(gòu)件已基本確定,其三維坐標(biāo)已知。為方便數(shù)學(xué)描述,可以對(duì)平臺(tái)上的三維空間進(jìn)行分解,對(duì)其上設(shè)備根據(jù)外形進(jìn)行立方體規(guī)整并采用正交布管。合理簡(jiǎn)化系統(tǒng)模型,可提高計(jì)算效率,同時(shí)也符合三維空間管道設(shè)計(jì)的實(shí)際要求。

1.1 平臺(tái)空間描述

海洋平臺(tái)的管路空間主要由管路、管端、彎頭和障礙物等組成。考慮正交布管,該三維空間中某單元設(shè)置如圖1所示。

將圖中空間單元立方體的1～8個(gè)頂點(diǎn)作為管路節(jié)點(diǎn),其中1號(hào)節(jié)點(diǎn)和7號(hào)節(jié)點(diǎn)分別為管路的起

圖1 空間單元示意圖

針對(duì)海洋平臺(tái)管路系統(tǒng)設(shè)計(jì)周期長(zhǎng)、效率低等問題,以平臺(tái)管路優(yōu)化布置為研究對(duì)象,以管段長(zhǎng)度、彎頭數(shù)和架設(shè)高度為控制目標(biāo)建立數(shù)學(xué)模型,對(duì)通用自適應(yīng)遺傳算法進(jìn)行改進(jìn),設(shè)計(jì)了獎(jiǎng)勵(lì)函數(shù)處理約束條件。模擬計(jì)算結(jié)果表明,改進(jìn)的自適應(yīng)遺傳算法進(jìn)化時(shí)間大幅度減少,穩(wěn)定性和進(jìn)化效率明顯增強(qiáng);基于該算法的海洋平臺(tái)管路優(yōu)化設(shè)計(jì)方法,可用于進(jìn)行多約束條件下的平臺(tái)管路設(shè)計(jì),可以提高設(shè)計(jì)效率、降低管路建設(shè)成本,并為海洋平臺(tái)管路三維優(yōu)化設(shè)計(jì)軟件的開發(fā)提供理論支持。

海洋平臺(tái) 管路優(yōu)化設(shè)計(jì) 遺傳算法 獎(jiǎng)勵(lì)函數(shù)

對(duì)于海上石油鉆井,由于海洋平臺(tái)允許的可變載荷有一定的限制范圍且平臺(tái)上空間有限,要求管匯系統(tǒng)必須設(shè)計(jì)得結(jié)構(gòu)緊湊、質(zhì)量輕、效率高;又由于平臺(tái)吊機(jī)能力的限制,還要求整個(gè)管匯系統(tǒng)模塊化程度高,安裝連接方便且節(jié)省時(shí)間和空間。海洋平臺(tái)管匯系統(tǒng)的上述特點(diǎn),使得高自動(dòng)化程度優(yōu)化布管設(shè)計(jì)手段和程序的研究開發(fā)更為迫切。管路布置設(shè)計(jì)是管系詳細(xì)設(shè)計(jì)階段的重要組成部分,目前對(duì)管路布置設(shè)計(jì)的研究已經(jīng)取得了一定成果,提出了迷宮法、逃逸法、遺傳算法等方法[1-3],其中的遺傳算法以目標(biāo)函數(shù)值為搜索依據(jù),通過群體優(yōu)化搜索和執(zhí)行遺傳運(yùn)算實(shí)現(xiàn)種群的不斷進(jìn)化,對(duì)解決此類非線性的組合優(yōu)化問題具有一定的優(yōu)勢(shì),但傳統(tǒng)的遺傳算法在實(shí)際應(yīng)用過程中存在不穩(wěn)定或局部收斂的問題,設(shè)計(jì)效率較低。本文將改進(jìn)的自適應(yīng)遺傳算法應(yīng)用于海洋平臺(tái)三維管匯布置的優(yōu)化設(shè)計(jì),在保證管匯系統(tǒng)有效運(yùn)行的前提下,尋求管道建設(shè)和運(yùn)行成本最小的方案,以管段長(zhǎng)度和彎頭數(shù)量為優(yōu)化設(shè)計(jì)變量構(gòu)造目標(biāo)函數(shù),從一組隨機(jī)產(chǎn)生的初始管路方案出發(fā),計(jì)算可行管匯布置方案的目標(biāo)函數(shù)值,通過不斷的搜索評(píng)價(jià)和遺傳操作,進(jìn)化、形成目標(biāo)函數(shù)最小的布置方案,實(shí)現(xiàn)管匯布置的優(yōu)化設(shè)計(jì)。點(diǎn)Si(a,b,c)和終點(diǎn)Ei(x,y,z),取由Si點(diǎn)指向Ei點(diǎn)的方向?yàn)檎?寫出該空間單元的節(jié)點(diǎn)管路連接路徑距離矩陣如下。

該矩陣為一反對(duì)稱矩陣,它反映了三維空間管路節(jié)點(diǎn)的正交連接性和路徑距離。通過該矩陣可進(jìn)行正交驗(yàn)證,在優(yōu)化迭代過程中由算法生成管路組群,并對(duì)每條管路的單元Di值進(jìn)行連乘,若乘積為“0”,則認(rèn)為是非可行布管,在算法的懲罰函數(shù)中進(jìn)行體現(xiàn)并將此路徑淘汰。

將平臺(tái)上已知設(shè)備據(jù)形狀虛擬為大小不等的立方體,可知每個(gè)立方體的頂點(diǎn)坐標(biāo),用包容盒描述立方體邊界,平臺(tái)設(shè)備簡(jiǎn)化及管道布置投影如圖2所示。

圖2 平臺(tái)布管投影示意圖

1.2 數(shù)學(xué)模型

進(jìn)行海洋平臺(tái)管匯布置時(shí),在考慮實(shí)現(xiàn)其功能的同時(shí)要重點(diǎn)考慮其經(jīng)濟(jì)因素,即管路越長(zhǎng)、彎頭越多,直接成本就越高,另外要注意到彎管會(huì)導(dǎo)致進(jìn)水壓及背壓的增加,水頭損失會(huì)隨之增大,因此管匯布置優(yōu)化目的可簡(jiǎn)化為,在控制管路長(zhǎng)度的同時(shí)盡量減少?gòu)濐^的使用量,即設(shè)計(jì)目標(biāo)變量為各管段長(zhǎng)度di和彎頭數(shù)量n,構(gòu)造優(yōu)化過程目標(biāo)函數(shù)如下:式(1)中:w1、w2、w3分別為管段長(zhǎng)度、彎頭數(shù)量和管段架設(shè)高度的權(quán)重系數(shù);DSE為管道起點(diǎn)至終點(diǎn)的直線距離;di為管段長(zhǎng)度;n為彎頭數(shù)量;hi為管段架設(shè)高度;r為懲罰項(xiàng)。

由于搜索空間不連續(xù),目標(biāo)函數(shù)值波動(dòng)較大,設(shè)計(jì)了選擇函數(shù)以提高進(jìn)化的穩(wěn)定性,構(gòu)造個(gè)體適應(yīng)度函數(shù)fi為

為保證布管方案實(shí)際可行,需使之滿足約束條件的限制,即約束范圍內(nèi)的個(gè)體被認(rèn)為是優(yōu)良個(gè)體,約束適應(yīng)度反映了個(gè)體對(duì)約束條件的滿足程度。據(jù)此,本文主要針對(duì)管道與設(shè)備的碰撞干涉約束,采用區(qū)域搜索法以管道與包容盒的相交干涉為判據(jù)對(duì)碰撞干涉進(jìn)行檢查,其步驟為:先對(duì)管路與設(shè)備的碰撞進(jìn)行檢查,然后進(jìn)行管路間的碰撞檢查,最后對(duì)檢查出碰撞干涉的管路進(jìn)行處理。該處理過程通過實(shí)時(shí)修改管路距離矩陣,將矩陣中的相關(guān)距離修改為一個(gè)極大值,使其成為不可行路徑來實(shí)現(xiàn)。

2 遺傳算法的改進(jìn)

管匯系統(tǒng)優(yōu)化設(shè)計(jì)的核心問題就是求解約束條件下的目標(biāo)最優(yōu)化問題。遺傳算法對(duì)所求解問題的數(shù)學(xué)模型要求不高,能夠進(jìn)行自適應(yīng)概率下的全局搜索,但其搜索效率及精度受到多方面條件的限制,因此需要對(duì)通用自適應(yīng)遺傳算法進(jìn)行改進(jìn),使之能較好地完成設(shè)計(jì)任務(wù)。

2.1 約束條件的處理

常用的各種遺傳算法均采用懲罰函數(shù)法對(duì)約束條件進(jìn)行處理,通過削弱不良染色體的競(jìng)爭(zhēng)力降低其遺傳概率[4-5],然而該逆向懲罰法極易造成相當(dāng)數(shù)量個(gè)體的適應(yīng)度值劇烈波動(dòng),特別是在進(jìn)化的中后期,超出約束范圍的個(gè)體仍有較大的規(guī)模,基于概率的選擇機(jī)制不能迅速將其淘汰,同時(shí)由于個(gè)體中可能存在優(yōu)良的基因片段,因此也有保留進(jìn)行交叉操作的必要;而逆向懲罰后,部分染色體個(gè)體適應(yīng)度值迅速背離,并在一段時(shí)間進(jìn)化后回歸,多次的交叉、變異操作將再次導(dǎo)致大量此類染色體的產(chǎn)生,體現(xiàn)為種群平均適應(yīng)度值波動(dòng)劇烈,算法的搜索和收斂效率低。

因此,本文通過設(shè)計(jì)獎(jiǎng)勵(lì)函數(shù)實(shí)現(xiàn)對(duì)約束條件的處理,即構(gòu)造動(dòng)態(tài)獎(jiǎng)勵(lì)函數(shù)對(duì)優(yōu)良染色體進(jìn)行操作,且獎(jiǎng)勵(lì)幅度隨進(jìn)化過程適當(dāng)增大,使約束范圍內(nèi)個(gè)體適應(yīng)度隨遺傳進(jìn)程逐代減小,體現(xiàn)其競(jìng)爭(zhēng)遺傳的優(yōu)勢(shì),控制不良個(gè)體的適應(yīng)度波動(dòng),提高進(jìn)化收斂效率。本獎(jiǎng)勵(lì)函數(shù)法中的選擇算子Ci可表示為

式(3)中:g為迭代次數(shù);r∈(0,1)為獎(jiǎng)勵(lì)指數(shù);p為管系中超出約束的管段數(shù)目。

2.2 交叉、變異操作

遺傳算法依賴的概率參數(shù)主要有交叉概率Pc和變異概率Pm。Pc決定著算法的探測(cè)能力,Pc越大越容易產(chǎn)生新個(gè)體,但過大的Pc將導(dǎo)致遺傳模式的混亂,使具有高適應(yīng)度值的個(gè)體結(jié)構(gòu)很快被破壞,而過小的Pc則會(huì)使搜索過程緩慢;變異操作可保持種群的多樣性,Pm太小難以產(chǎn)生新的基因塊,太大則會(huì)使遺傳算法變成隨機(jī)搜索。隨著種群的進(jìn)化,個(gè)體適應(yīng)度開始分化,交叉、變異的概率值也需隨進(jìn)化過程而變化,初始化兩概率值易造成搜索停滯不前或陷入局部最優(yōu)狀況,理論上尚沒有成熟的方法用于選取算法的最優(yōu)參數(shù)。

Srinvivas等提出一種自適應(yīng)遺傳算法,其特點(diǎn)主要體現(xiàn)在其變化的交叉、變異概率設(shè)計(jì)上[6],但在優(yōu)化的中后期進(jìn)化趨于停滯時(shí),對(duì)兩概率的處理也不能有效地使算法擺脫困境。因此本文對(duì)其進(jìn)行了改進(jìn),當(dāng)種群適應(yīng)度值趨于一致時(shí)適當(dāng)?shù)卦黾覲c、Pm值,當(dāng)種群適應(yīng)度值較分散時(shí)減小Pc、Pm值,同時(shí)對(duì)適應(yīng)度值低于平均值的個(gè)體減小Pc以避免不良基因片段進(jìn)入其他個(gè)體,并增大Pm使其得以更新,增加種群的多樣性。另外,當(dāng)進(jìn)化趨于停滯時(shí),種群多樣性依靠適應(yīng)度值低于平均值的個(gè)體的較大變異概率來調(diào)整,可跳出局部最優(yōu)狀況,改進(jìn)后的兩概率值公式為

式(4)、(5)中:fmax為適應(yīng)度最大值;favg為群體適應(yīng)度平均值;fi為變異個(gè)體的適應(yīng)度值;Pc0為初始交叉參數(shù);Pm0為初始變異參數(shù)。

2.3 評(píng)價(jià)準(zhǔn)則

算法改進(jìn)后的效果可引入評(píng)價(jià)機(jī)制進(jìn)行評(píng)估, DeJong曾提出兩個(gè)測(cè)度用于定量分析遺傳算法從而對(duì)算法進(jìn)行評(píng)估[7],其中離線性能測(cè)度表達(dá)式為

式(6)表明,其收斂特性的離線性能測(cè)試是特定時(shí)刻最佳性能的累積平均。然而由于進(jìn)行交叉、變異操作,適應(yīng)度最優(yōu)值隨進(jìn)化波動(dòng)的情況普遍存在,用離線性能常常不能準(zhǔn)確判別算法的優(yōu)劣,而適應(yīng)度最優(yōu)值的波動(dòng)程度及算法對(duì)波動(dòng)的抑制體現(xiàn)著算法的選擇進(jìn)化能力,其波動(dòng)幅值統(tǒng)計(jì)量可用于對(duì)性能進(jìn)行評(píng)價(jià),因此構(gòu)造了波動(dòng)幅值函數(shù)作為其離線性能并表達(dá)為

進(jìn)化的終止條件可以種群目標(biāo)函數(shù)平均值與最小值之差為參考,表達(dá)式為

式(8)中:e取實(shí)驗(yàn)值,當(dāng)進(jìn)化至某代種群目標(biāo)函數(shù)值滿足式(8)時(shí)即可跳出循環(huán)取得尋優(yōu)結(jié)果。

基于遺傳算法的管路布置優(yōu)化計(jì)算流程如圖3所示,在M atlab環(huán)境下編制應(yīng)用程序包以完成管路布置的優(yōu)化任務(wù)。

圖3 基于遺傳算法的管路布置優(yōu)化流程

3 模擬實(shí)驗(yàn)分析

3.1 參數(shù)設(shè)置

某平臺(tái)管系設(shè)計(jì)模型空間如圖4所示,該部分平臺(tái)上布置有7臺(tái)設(shè)備,管路起點(diǎn)、終點(diǎn)坐標(biāo)分別為(0,0,0)、(100,100,100),根據(jù)包容盒規(guī)則對(duì)障礙物所占據(jù)空間進(jìn)行簡(jiǎn)化,形成圖中模型空間,靠近管路起點(diǎn)、終點(diǎn)的各包容盒對(duì)角點(diǎn)坐標(biāo)見表1。

圖4 平臺(tái)空間模型

表1 包容盒對(duì)角點(diǎn)坐標(biāo)(x,y,z)

圖4中一條直線管路路徑對(duì)應(yīng)種群中的一個(gè)染色體,其坐標(biāo)由起點(diǎn)、終點(diǎn)的彎頭坐標(biāo)組成,采用可變長(zhǎng)度的實(shí)數(shù)編碼。染色體的基因相互關(guān)聯(lián),一個(gè)節(jié)點(diǎn)坐標(biāo)分別同時(shí)存在于前、后兩個(gè)染色體的首、尾基因片段中,隨機(jī)生成初始種群。算法參數(shù)設(shè)置見表2,一般管路布置應(yīng)盡量處于低重力勢(shì)能區(qū)域,但由于對(duì)不同的管路架設(shè)有不同的設(shè)計(jì)要求,這里取A、B兩種優(yōu)化策略進(jìn)行對(duì)比分析,且B策略下較A策略下增大了彎頭數(shù)權(quán)重系數(shù)而減小了架設(shè)高度權(quán)重系數(shù),即放松了對(duì)架設(shè)高度的限制,以滿足不同管路的設(shè)計(jì)要求。

表2 不同優(yōu)化策略下管路布置算法參數(shù)設(shè)置

3.2 模擬實(shí)驗(yàn)結(jié)果分析

對(duì)兩種優(yōu)化策略程序分別運(yùn)算10次以減小偶然性的影響,優(yōu)化結(jié)果如圖5所示。路徑1、路徑2為A優(yōu)化策略下得出的兩種結(jié)果,路徑3為B優(yōu)化策略下得出的結(jié)果。可以看出,由本文優(yōu)化算法得出的不同策略下的各條路徑均可合理地避開障礙物且靠近空間邊界,本文算法取得了架設(shè)長(zhǎng)度較短和彎頭數(shù)量較少的布管結(jié)果。在B優(yōu)化策略下,由于減小了架設(shè)高度的權(quán)重系數(shù),路徑3較路徑1和2的管路平均架設(shè)高度有所增大但彎頭數(shù)量明顯減少,因此路徑3屬于合理布管路徑,且適用于需在高處架設(shè)的管路,而路徑1和2也各具特點(diǎn),設(shè)計(jì)人員可根據(jù)現(xiàn)場(chǎng)具體情況進(jìn)行選用。優(yōu)化結(jié)果表明,本文利用基于改進(jìn)后算法的路徑優(yōu)化方法較好地完成了一定約束條件下的海洋平臺(tái)管匯設(shè)計(jì)任務(wù)。

圖5 管路布置優(yōu)化結(jié)果

3.3 算法改進(jìn)效果討論

對(duì)于算法的改進(jìn)效果,提取改進(jìn)前、后的算法評(píng)估值繪制波動(dòng)幅度比較曲線如圖6所示。可以看出,改進(jìn)后算法的進(jìn)化時(shí)間較改進(jìn)前大幅度減少;同時(shí)可以看出,改進(jìn)后算法對(duì)算法波動(dòng)的抑制能力具有明顯優(yōu)勢(shì),體現(xiàn)為提高了種群中劣勢(shì)個(gè)體的淘汰效率,增大了優(yōu)勢(shì)基因片段在個(gè)體基因段中所占比例,因此進(jìn)化初期適應(yīng)度波動(dòng)幅值迅速衰減,在40代前、后種群進(jìn)化方向趨于一致,而改進(jìn)前算法需進(jìn)化至130代才能達(dá)到該抑制水平,可見改進(jìn)后算法的進(jìn)化穩(wěn)定性明顯增強(qiáng)且效率更高。

圖6 改進(jìn)前、后算法波動(dòng)幅值比較

改進(jìn)前、后算法的收斂過程曲線分別如圖7、8所示。對(duì)于管路優(yōu)化設(shè)計(jì)問題,改進(jìn)前算法收斂緩慢,頻繁出現(xiàn)劇烈跳動(dòng);改進(jìn)后算法收斂迅速平穩(wěn),跳動(dòng)的頻率及幅度較小,50代前、后目標(biāo)函數(shù)值趨于穩(wěn)定。由以上分析可見,本文算法的改進(jìn)效果顯著,效率高且收斂迅速、穩(wěn)定。

4 結(jié)論

(1)本文提出的基于改進(jìn)自適應(yīng)遺傳算法的海洋平臺(tái)管路優(yōu)化設(shè)計(jì)方法,可用于進(jìn)行多約束條件下的平臺(tái)管路設(shè)計(jì),可以提高設(shè)計(jì)效率,并有效控制管段長(zhǎng)度和彎頭數(shù)量,從而降低建設(shè)成本。

(2)引入評(píng)估函數(shù)對(duì)改進(jìn)前、后算法的進(jìn)化過程和優(yōu)化結(jié)果進(jìn)行對(duì)比分析表明,改進(jìn)的遺傳算法收斂迅速且穩(wěn)定,群體適應(yīng)度波動(dòng)小、頻率低,遺傳進(jìn)化效率顯著提高,性能優(yōu)于改進(jìn)前算法。

[1] KN IATA.Optimization of three-dimensional piperouting[J].Ship Technology Research,2000,47:111-114.

[2] KANGSSS,SEHYUNM,HANS SH.A design expert system for auto-routing of ship pipes[J].J of Ship Production,1999,15(1): 1-9.

[3]SANDURKAR S,CHENW.GARPUS-genetic algorithms based pipe routing using tessellated objects[J].Computers in Industry,1999, 38:209-223.

[4] 關(guān)曉晶,魏立新,楊建軍.基于混合遺傳算法的油田注水系統(tǒng)運(yùn)行方案優(yōu)化模型[J].石油學(xué)報(bào),2005,26(3):114-117.

[5] 王躍宣,劉連臣,牟盛靜,等.處理帶約束的多目標(biāo)優(yōu)化進(jìn)化算法[J].清華大學(xué)學(xué)報(bào):自然科學(xué)版,2005,45(1):103-106.

[6] 王小平,曹立明.遺傳算法理論、應(yīng)用與軟件實(shí)現(xiàn)[M].西安:西安交通大學(xué)出版社,2003.

[7] 高強(qiáng),呂文芝,杜小山.遺傳算法優(yōu)化性能評(píng)價(jià)準(zhǔn)則研究[J].西安交通大學(xué)學(xué)報(bào),2006,40(7):803-806.

(編輯:張金棣)

Optimized piping design method for offshore platform based on improved adaptive genetic algorithm

Yi Peng Liu Yancong Shi Yongjun Ren Hongwei Zheng Huanpeng
(College of Mechanical and Electronic Engineering, China University of Petroleum,Shandong,257061)

For solving the problems of low efficiency and long period of design in offshore platform piping design process,the manifold optimizing layout on platform was taken as research object and the mathematical model was established with the controlled parameters such as pipe length,number of elbow s and erection height.The common adaptive genetic algorithm was imp roved and the encouraging function was designed to deal with constraints.The simulation results have indicated that,the evolution time of the optimized algorithm reduced significantly,besides,the evolution stability and efficiency were imp roved obviously;the optimized method based on the imp roved adaptive genetic algorithm can accomplish the design of piping on offshore platform with constraints,so that the design efficiency raised and construction costs deduced.This piping design method rovided the theoretical supports for development of three-dimensional optimized piping design software for offshore platform.

offshore platform;optimized piping design;genetic algorithm;encouraging function

＊國(guó)家高技術(shù)研究發(fā)展計(jì)劃(863計(jì)劃)經(jīng)費(fèi)資助項(xiàng)目(2006AA 09A 104)及中國(guó)石油大學(xué)(華東)研究生創(chuàng)新基金資助項(xiàng)目(S2008-13)部分研究成果。

伊鵬,男,中國(guó)石油大學(xué)(華東)機(jī)電工程學(xué)院在讀博士研究生,主要研究方向?yàn)楹Ｑ筱@井平臺(tái)管匯系統(tǒng)設(shè)計(jì)。地址:山東省東營(yíng)市中國(guó)石油大學(xué)(郵編:257061)。電話:0546-7878345。E-mail:yipupc@yahoo.cn。

2009-05-19 改回日期:2009-10-26