09.03.2014 Views

L - Transportu

L - Transportu

L - Transportu

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

PRACE NAUKOWE POLITECHNIKI WARSZAWSKIEJ<br />

z. 78 Transport 2011<br />

Marek Szczotka<br />

Wydzia Zarzdzania i Informatyki<br />

Akademia Techniczno-Humanistyczna<br />

MODELOWANIE I ANALIZA UKADU<br />

DO OBSUGI MODUÓW INSTALOWANYCH<br />

NA DNIE MORZA<br />

Rkopis dostarczono, marzec 2011<br />

Streszczenie: W pracy przedstawiono model ukadu przeznaczonego do instalacji offshore moduów<br />

niezbdnych do transportowania ropy i gazu z odwiertów w dnie morskim. System jest montowany na<br />

specjalnych jednostkach pywajcych, wyposaonych w ukady dynamicznego pozycjonowania.<br />

Jednym z istotniejszych elementów ukadu jest system AHC (Active Heave Compensation), który<br />

kompensuje ruchy statku spowodowane falowaniem. Zastosowano uproszczony ukad pneumatycznohydrauliczny,<br />

wraz z odpowiednimi algorytmami sterowania prac wcigarek. adunek jest<br />

opuszczany na dno za pomoc wcigarki z kompensacj. Zamodelowano równie dodatkowe liny<br />

prowadzce adunek. Przedstawiono wyniki przykadowych oblicze numerycznych, dla ukadu<br />

pracujcego w zadanych warunkach, przy wczonym oraz wyczonym ukadzie prowadzenia<br />

moduu w trakcie instalacji.<br />

Sowa kluczowe: urzdzenia offshore, instalacja moduów, modelowanie<br />

1. SYSTEMY DO TRANSPORTU MODUÓW (MHS)<br />

Nieustajcy rozwój gospodarek wielu krajów wiata generuje zwikszone<br />

zapotrzebowanie na surowce energetyczne. Dlatego ju od kilku dekad znaczca cz<br />

energii pochodzi ze róde znajdujcych si pod dnami mórz i oceanów. Ich wydobycie<br />

staje si coraz bardziej kosztowne, poniewa istniejce zasoby w akwenach o maych<br />

gbokociach (do kilkuset metrów) s ju eksploatowane od duszego czasu. Nowe zoa<br />

ropy i gazu ziemnego, dla których powstaje obecnie wikszo infrastruktury zwizanej<br />

z ich eksploatacj, charakteryzuj si wikszymi wymaganiami w stosunku do systemów<br />

budowanych w minionych dekadach: du gbokoci (2000m i wicej), wysokim stanem<br />

mórz i oceanów, przy których prace wydobywcze musz by prowadzone w sposób cigy,<br />

a take silnymi prdami morskimi i niskimi temperaturami (obszary arktyczne). Coraz<br />

bardziej ekstremalne warunki wymagaj coraz wikszych nakadów przy jednoczesnym<br />

skróceniu czasu eksploatacji pól (s one czsto coraz mniejsze i coraz bardziej<br />

rozproszone) oraz obsugujcych ich urzdze. Czynnik ten bdzie w znacznym stopniu


86 Marek Szczotka<br />

wpywa midzy innymi na ceny paliw w nastpnych latach. Z drugiej strony, te trudne<br />

warunki pracy wymuszaj szybki rozwój technologii rodków transportu i maszyn,<br />

niezbdnych do funkcjonowania brany wydobywczej. Bezawaryjna praca tych urzdze,<br />

niejednokrotnie w trudnych warunkach pogodowych (silne falowanie, prdy morskie oraz<br />

wiatr), wymaga specjalnego podejcia do zagadnie niezawodnoci i bezpieczestwa [1-3].<br />

Czstym rozwizaniem przy rónorodnych pracach przeadunkowych jest stosowanie<br />

urawi offshore, montowanych na statkach i platformach. Maszyny te nale do<br />

najbardziej efektywnych i uniwersalnych, charakteryzuj si licznymi zaletami. Jednak nie<br />

zawsze tego typu urzdzenie jest najodpowiedniejsze do wykonania pewnych prac, poza<br />

tym ich cena jest bardzo wysoka. Ponadto urawie czsto w tym samym czasie wykonuj<br />

inne zadania, wic zachodzi potrzeba uycia drugiego dwigu. W artykule przedstawiono<br />

dedykowany system do opuszczania i podnoszenia moduów uywanych przy eksploatacji<br />

morskich rurocigów na polach naftowych i gazowych. Moe on stanowi alternatyw dla<br />

urawi, dziki dostosowaniu konstrukcji do specyfiki prac instalacyjnych. W szczególnoci<br />

moe umoliwia prowadzenie adunku za pomoc ukadów dodatkowych lin i wcigarek,<br />

co jest trudne stosujc klasyczny uraw.<br />

W pracy krótko omówiono model systemu MHS (Multi-purpose Handling System),<br />

który jest wyposaony w ukad wcigarek z systemami AHC, a take zawiera wcigarki<br />

pomocnicze, stosowane do prowadzenia adunku. Poniej na Rys. 1 przedstawiono widok<br />

urzdzenia.<br />

wcigarki<br />

kabinasterowania<br />

moonpool<br />

wózektransportowy<br />

Rys. 1. Widok systemu MHS od strony torowiska umoliwiajcego poziomy transport moduów<br />

(urzdzenie zaprojektowane i wyprodukowane przez AXTech AS)<br />

Struktura nona („wiea”) stanowi rodzaj konstrukcji ramowej, na której zamocowano<br />

wcigarki dla lin prowadzcych, gówn wcigark z ukadem AHC, kabin sterownicz<br />

oraz wiele innych komponentów dodatkowych (krki, szyny mechanizmu jazdy ramy<br />

prowadzcej adunek, cigi komunikacyjne, oprzyrzdowanie sterujce). MHS zosta


Modelowanie i analiza ukadu do obsugi moduów instalowanych na dnie morza 87<br />

zamontowany na statku uczestniczcym w pracach instalacyjnych i obsudze wydobycia<br />

ropy i gazu, przy czym do stosowania omawianego systemu niezbdny jest specjalny kana<br />

w kadubie statku (moonpool). Moduy s opuszczane i podnoszone z dna przez ten kana,<br />

posiadajcy midzy innymi waciwoci eliminujce efekty dynamiczne pojawiajce si<br />

podczas przechodzenia adunku przez lustro wody. Wykonanie podobnej operacji za<br />

pomoc klasycznego urawia offshore (a wic w przypadku kiedy adunek jest opuszczany<br />

przez lustro wody poza obrysem kaduba), wymaga znacznie wikszego zapasu udwigu<br />

urawia, poniewa operacje takie cechuj si znacznie wikszymi siami dynamicznymi.<br />

Zatem, przy tych samych warunkach pogodowych, prezentowany MHS wymaga<br />

mniejszych nonoci ni wymagaby uraw pokadowy. Moduy s transportowane po<br />

pokadzie statku z hangaru (znajdujcego si w innej czci pokadu) do wiey,<br />

wykorzystujc specjalne wózki i system szyn.<br />

wcigarkaAHC35t<br />

HPU<br />

Rys. 2. MHS od strony generatora i wcigarki z ukadem automatycznej kompensacji<br />

Przedstawiony w pracy ukad jest przeznaczony do opuszczania moduów o wymiarach<br />

6.5m x 6.5m x 8m, na gboko do 1500m. Maksymalna zdolno kompensacji falowania<br />

zostaje osignita przy wysokoci fal 5.0m. Cakowita masa wasna urzdzenia to 120t<br />

(wraz z osprztem i linami). Centraln czci systemu jest specjalnie zaprojektowany<br />

ukad napdowy dla wcigarek (Rys. 2): gównej z ukadem AHC o udwigu nominalnym<br />

35t, dwóch specjalnych dla utrzymania staego nacigu lin prowadzcych modu w trakcie<br />

podnoszenia, obsugujcej ram prowadzc modu przez kadub statku (zabezpieczenie<br />

przez uderzeniami o ciany boczne) oraz pomocnicze dla obsugi pomostu zamykajcego<br />

moonpool i kurtyn bocznych. Podstawow funkcj napdu wcigarki z systemem AHC jest<br />

zapobieganie przed uderzeniem adunku o dno w czasie jego podnoszenia i opuszczania.<br />

Wcigarka kompensuje ruchy nurzania i koysania statku spowodowane falowaniem,<br />

poprzez odpowiedni dobór dugoci liny. Klasyczny ukad AHC wcigarki z napdem<br />

elektrycznym opisano midzy innymi w pracy [4].<br />

W niniejszej pracy przedstawiono model matematyczny umoliwiajcy symulacj<br />

jednej z moliwych operacji wykonywanych za pomoc MHS. Skupiono si na operacji<br />

podnoszenia i opuszczania adunku przy dodatkowym istnieniu lin prowadzcych.


88 Marek Szczotka<br />

W analizie uwzgldniono efekty dynamiczne wywoane zarówno ruchem unoszenia statku<br />

w wyniku jego ruchu podczas falowania morza, ale take oddziaywania lin. Wymuszenia<br />

pochodz take od si hydrodynamicznych spowodowanych prdkociami czsteczek<br />

wody (od falowania i prdów morskich).<br />

2. PODSTAWOWE ZAOENIA I OPIS MODELU UKADU<br />

Model matematyczny, umoliwiajcy symulacj pracy ukadu i ocen wpywu<br />

wybranych parametrów na zachowanie si adunku i powstajce w ukadzie siy<br />

dynamiczne, zbudowano uwzgldniajc ponisze zaoenia:<br />

ruch statku powstay w wyniku falowania morskiego traktuje si jako znany<br />

i opisany dowoln funkcj zalen od czasu i parametrów stanu morza<br />

i jednostki pywajcej,<br />

nie uwzgldnia si podatnoci struktury nonej wiey oraz pokadu statku,<br />

napd wcigarek lin prowadzcych zapewnia wzgldnie stay nacig lin,<br />

realizowany przez dostarczenie staego cinienia oleju do silników,<br />

wcigarka gówna posiada hybrydowy ukad napdowy (hydro-pneumatyczny),<br />

natomiast prdko obrotowa bbna jest wymuszana kinematycznie,<br />

uwzgldnia si podatno lin prowadzcych oraz podatno gównej liny nonej,<br />

poprzez jej dyskretyzacj metod sztywnych elementów skoczonych (SES),<br />

adunek traktuje si jako bry sztywn o szeciu stopniach swobody,<br />

prowadnice adunku mog by traktowane jako elementy sprysto-tumice<br />

o dowolnie definiowalnych charakterystykach.<br />

Ruch statku, na którym zamontowane jest urzdzenie, jest okrelony jeli znane s<br />

skadowe wektora q D , bdce wspórzdnymi uogólnionymi statku:<br />

D<br />

<br />

<br />

s<br />

z<br />

<br />

q f t ,RAO H , T<br />

(1)<br />

gdzie:<br />

H S , T z - charakterystyczna wysoko i okres fali,<br />

q <br />

T<br />

D xD<br />

yD<br />

zD<br />

D D D<br />

- wspórzdne okrelajce ruch statku,<br />

xD<br />

xD<br />

t ,..., D<br />

D<br />

t , gdzie x D, yD,<br />

zD<br />

s przemieszczeniami jednostki<br />

wzdu osi ukadu bezwadnociowego, natomiast kty D , D , D<br />

s<br />

ktami obrotu Eulera ZYX wzgldem tego ukadu,<br />

RAO - funkcje przejcia statku (Response Amplitude Operators).<br />

Schemat ukadu przedstawiono na rysunku 3.


Modelowanie i analiza ukadu do obsugi moduów instalowanych na dnie morza 89<br />

wiea<br />

kreklinygównej<br />

profilfali<br />

{D}<br />

a <br />

y<br />

moonpool<br />

v <br />

ruchczstekwody<br />

{L}<br />

linagówna<br />

linyprowadzce<br />

modu<br />

z {0}<br />

x<br />

dno<br />

Rys. 3. Schemat ukadu MHS<br />

Poniewa zaoono, e konstrukcja nona („wiea”) jest nieodksztacalna, trajektorie<br />

ruchu krków mona wyznaczy stosujc przeksztacenia jednorodne. Zatem ruch<br />

elementu liny stycznego do krka zamocowanego na wiey mona okreli za pomoc<br />

operacji [5]:<br />

(0)<br />

r T <br />

<br />

(2)<br />

s, l D t rs,<br />

l<br />

(0)<br />

gdzie: r<br />

s,l<br />

- wspórzdne punktu zejcia liny z krka w ukadzie globalnym { 0}<br />

,<br />

r s , l const - znane wspórzdne w ukadzie statku {D}<br />

,<br />

T D - macierz transformacji jednorodnej wspórzdnych z ukadu {D}<br />

do { 0}<br />

.<br />

Macierz transformacji jednorodnych jest zdefiniowana zgodnie z konwencj któw<br />

ZYX Eulera, w sposób podany w [5].<br />

Liny prowadzce adunek, jak równie gówn lin non, dyskretyzowano za pomoc<br />

metody sztywnych elementów skoczonych, która szczegóowo opisana jest w monografii<br />

(l)<br />

[6]. Pojedynczy sztywny element skoczony ses i jest opisany przez wspórzdne<br />

uogólnione (Rys. 4):<br />

i<br />

( l)<br />

i<br />

i<br />

( l)<br />

i<br />

i<br />

( l)<br />

i<br />

~( l)<br />

qi<br />

<br />

<br />

x<br />

( l)<br />

i<br />

y<br />

( l)<br />

i<br />

z<br />

( l)<br />

i<br />

<br />

( l)<br />

i<br />

<br />

( l)<br />

i<br />

<br />

( l)<br />

T<br />

i<br />

( l)<br />

( l)<br />

( l)<br />

gdzie: x , y , z - przemieszczenia elementu i w ukadzie {D}<br />

,<br />

, , - kty obrotu Eulera elementu wzgldem ukadu {D}<br />

,<br />

(3)


90 Marek Szczotka<br />

l - numer liny prowadzcej adunek (w pracy zakadano l 1, 2 ).<br />

a)<br />

( l)<br />

sesi<br />

1<br />

( )<br />

y ˆ l<br />

i<br />

( )<br />

(l)<br />

z ˆ l<br />

i est<br />

( ) i<br />

x ˆ l<br />

i<br />

(l)<br />

ses i<br />

( )<br />

yˆ l<br />

i<br />

( )<br />

{} l<br />

( )<br />

zˆ l i<br />

i<br />

( )<br />

xˆ l<br />

i<br />

( l)<br />

( )<br />

esti<br />

1 yˆ l<br />

i<br />

1<br />

( l)<br />

( )<br />

zˆ l sesi<br />

1<br />

i 1<br />

( )<br />

xˆ l<br />

i<br />

1<br />

b)<br />

( )<br />

zˆ l<br />

i<br />

( )<br />

z ˆ l<br />

i<br />

~ ( l)<br />

r<br />

L,<br />

i<br />

( )<br />

xˆ l<br />

i<br />

( l)<br />

sesi<br />

1<br />

~ ( l)<br />

r i<br />

x<br />

i<br />

( )<br />

{} l<br />

i<br />

( )<br />

yˆ l<br />

i<br />

(l)<br />

ses i<br />

( )<br />

zˆ l i 1 ( )<br />

{} l<br />

i1<br />

( )<br />

xˆ l1 i <br />

( )<br />

ˆ l<br />

~ ( l)<br />

r<br />

R,<br />

i1<br />

( )<br />

y ˆ l<br />

i<br />

~ ( l )<br />

( l)<br />

i<br />

x<br />

z<br />

i<br />

, y<br />

( l)<br />

i<br />

y<br />

{0}<br />

( )<br />

yˆ l1 i <br />

, z<br />

( l)<br />

i<br />

x<br />

Rys. 4. Dyskretyzacja metod sztywnych elementów skoczonych<br />

(l)<br />

a) sztywny element skoczony ses i fragmentu liny l (przed deformacj)<br />

( )<br />

b) ukad po deformacji: odksztacenie liniowe ~ l<br />

~ ( )<br />

i ktowe l<br />

elementu<br />

(l)<br />

Deformacja elementu sprysto-tumicego est i umoliwia, po przyjciu<br />

odpowiedniego modelu materiau, wyznaczenie si i momentów si wystpujcych<br />

( l)<br />

(l)<br />

w poczeniu elementów sesi<br />

1 i ses i liny. Uzewntrznienie si i momentów si<br />

wynikajcych z odksztace elementów, w odrónieniu od obliczania skadników równa<br />

ruchu pochodzcych od energii potencjalnej odksztace sprystych, pozwala na<br />

modelowanie dowolnych, nieliniowych charakterystyk sprystych i sprystoplastycznych<br />

materiau [7].<br />

Transformacji wspórzdnych z ukadu lokalnego, zwizanego ze sztywnym elementem<br />

(l)<br />

skoczonym ses , dokonuje si stosujc podobnie przeksztacenie:<br />

gdzie: x y z 1 T<br />

0 <br />

0<br />

i<br />

r -<br />

( l)<br />

i, loc<br />

<br />

( l<br />

<br />

i,<br />

)<br />

<br />

0<br />

)<br />

loc<br />

0<br />

( l)<br />

i,<br />

loc<br />

r i<br />

i,<br />

loc<br />

i,<br />

loc<br />

i<br />

(l)<br />

est i<br />

( l)<br />

( ) ( ) ( ) ( )<br />

0 ( ) ~<br />

l l l l<br />

r TD qD<br />

t Ti<br />

qi<br />

r<br />

Bi<br />

r<br />

(4)<br />

( l)<br />

i,<br />

loc<br />

T<br />

wektor wspórzdnych w ukadzie<br />

bezwadnociowym,<br />

r x y z 1 - wektor wspórzdnych w ukadzie lokalnym,<br />

<br />

~ ( l)<br />

~ ( l<br />

i qi<br />

T -<br />

(l)<br />

macierz transformacji z ukadu lokalnego ses i<br />

do ukadu statku {D}<br />

.


Modelowanie i analiza ukadu do obsugi moduów instalowanych na dnie morza 91<br />

Równania ruchu liny prowadzcej l wyprowadzono z równa Lagrange’a II rodzaju:<br />

V<br />

D<br />

(5)<br />

( l)<br />

( l)<br />

k<br />

Q<br />

( l)<br />

( l)<br />

k<br />

qk<br />

q<br />

k<br />

gdzie:<br />

( l)<br />

d T<br />

T<br />

<br />

k<br />

<br />

( l)<br />

( l)<br />

dt q<br />

k<br />

qk<br />

- operator Lagrange’a,<br />

E , V - energia kinetyczna i potencjalna,<br />

D - funkcja dyssypacji energii,<br />

Q - niepotencjalna sia uogólniona,<br />

<br />

<br />

(l)<br />

k<br />

k 0,...,n l<br />

, - n l jest liczb elementów skoczonych.<br />

(l)<br />

Równania ruchu elementu ses mona zapisa w postaci:<br />

gdzie:<br />

<br />

( l)<br />

k<br />

<br />

l<br />

l l l l<br />

q<br />

~ ( )<br />

k<br />

F q<br />

~ ( )<br />

k<br />

q<br />

~ ( )<br />

, ,<br />

k<br />

, q<br />

~ ( )<br />

k<br />

, q<br />

~ ( )<br />

<br />

1<br />

k 1<br />

A<br />

k<br />

t , q<br />

<br />

( l)<br />

<br />

( l)<br />

Mk<br />

0<br />

A<br />

k<br />

<br />

( l)<br />

<br />

- macierz mas,<br />

<br />

0 Ak<br />

<br />

<br />

<br />

( l )<br />

( l ) ( l<br />

,<br />

) ,<br />

( l<br />

M )<br />

k<br />

diag mk<br />

mk<br />

mk<br />

- podmacierz 3x3 o staych elementach,<br />

(l) ( l)<br />

( l)<br />

( l)<br />

(l)<br />

m<br />

k<br />

, I<br />

x, k<br />

, I<br />

y,<br />

k<br />

, Iz,<br />

k<br />

- masa i momenty bezwadnoci ses<br />

k<br />

,<br />

( l)<br />

( l)<br />

( l)<br />

( l)<br />

( l)<br />

podmacierz 3x3 o zmiennych elementach<br />

Ak<br />

Ak<br />

<br />

k<br />

, k<br />

, k<br />

-<br />

(l)<br />

(l)<br />

(l)<br />

(zalenych od , , ),<br />

(l)<br />

Fk<br />

<br />

k<br />

D<br />

k<br />

<br />

wektor prawych stron zawierajcy skadowe<br />

od si uogólnionych, ruchu unoszenia, si<br />

cikoci itp.<br />

<br />

(l)<br />

Warunki pracy urzdzenia okrelaj, e element ses<br />

n l<br />

jest zamocowany na dnie morza.<br />

W modelu realizuje si to poprzez uwzgldnienie równania:<br />

<br />

( l)<br />

( l)<br />

( l)<br />

r B const . (7)<br />

n , 0 n<br />

rn<br />

, loc<br />

l<br />

l<br />

<br />

(l)<br />

Równanie (7) reprezentuje poczenie kuliste, czce ses<br />

n l<br />

z dnem w punkcie<br />

( )<br />

o wspórzdnych lokalnych okrelonych w r l<br />

. n l , loc<br />

( l)<br />

Element ses 0 jest obciony si w linie prowadzcej. Ruch bbna wcigarki<br />

(l)<br />

prowadzcej okrelony jest wspórzdn uogólnion GW , któr wyznacza si z równania:<br />

<br />

( l)<br />

( l)<br />

( l)<br />

( l)<br />

( l)<br />

I T F r<br />

(8)<br />

GW<br />

GW<br />

l<br />

const<br />

l<br />

GW<br />

k<br />

(6)<br />

gdzie:<br />

(l)<br />

IGW<br />

- masowy moment bezwadnoci bbna wcigarki,


92 Marek Szczotka<br />

(l)<br />

T const - stay moment napdowy,<br />

r - promie podziaowy,<br />

(l)<br />

GW<br />

(l)<br />

Fl<br />

- sia w linie l .<br />

<br />

Zakada si, i adunek jest bry sztywn o szeciu stopniach swobody. Jego ruch jest<br />

okrelony w ukadzie bezwadnociowym przez równania ruchu wyprowadzone równie<br />

z równa Lagrange’a II rodzaju:<br />

<br />

A Lq L FL<br />

t,<br />

q L , q<br />

L <br />

(9)<br />

( L)<br />

A , ( L)<br />

( L)<br />

( L)<br />

T<br />

a <br />

<br />

k , j<br />

trB k<br />

H B j ,<br />

<br />

<br />

(L)<br />

H jest macierz pseudo-masow,<br />

q x y z ,<br />

( L)<br />

gdzie: L a k<br />

, j k,<br />

j 1,..., 6<br />

L<br />

<br />

L<br />

L<br />

( L)<br />

( L)<br />

k<br />

L<br />

L<br />

L<br />

T<br />

L<br />

( L)<br />

<br />

k<br />

<br />

B (L)<br />

B , B - macierz transformacji z {L}<br />

do { 0}<br />

,<br />

q<br />

F L - wektor prawych stron równa ruchu.<br />

<br />

adunek opuszczany jest za pomoc liny gównej. Aby ograniczy wpyw prdów<br />

morskich, stosuje si prowadzenie adunku z wykorzystaniem dodatkowych lin, Rys. 5.<br />

<br />

prdy morskie,<br />

fale<br />

lina<br />

prowadzca<br />

(1)<br />

F<br />

l<br />

lina nona<br />

{ 0}<br />

modu<br />

ẑ<br />

L<br />

ŷ<br />

L<br />

przeguby kuliste<br />

xˆ<br />

L<br />

(2)<br />

F<br />

l<br />

prowadnica<br />

lina<br />

prowadzca<br />

<br />

<br />

Rys. 5. adunek (modu) opuszczany przy wykorzystaniu lin prowadzcych<br />

<br />

W celu okrelenia si dziaajcych na adunek oraz lin l , zaoono, e prowadnice<br />

posiadaj sztywno i tumienie w kierunkach osi xˆ L oraz ẑ L . Obliczajc wspórzdne<br />

(l)<br />

pocztku i koca ses i bdcego w kontakcie z prowadnic:<br />

<br />

( l)<br />

( l)<br />

( l)<br />

r <br />

T<br />

a, i Bi<br />

ri<br />

, a,<br />

loc<br />

xa<br />

ya<br />

za<br />

1<br />

(10.1)<br />

( l)<br />

( l)<br />

( l)<br />

r B r x y z 1 T<br />

(10.2)<br />

b, i i i,<br />

b,<br />

loc<br />

<br />

( l)<br />

l<br />

gdzie: r <br />

0 0 T ( l)<br />

l<br />

, 0 0 T<br />

i,<br />

a,<br />

loc<br />

( l )<br />

i<br />

2<br />

1<br />

i,<br />

b,<br />

loc<br />

b<br />

b<br />

( l )<br />

i<br />

2<br />

1<br />

r ,<br />

b


Modelowanie i analiza ukadu do obsugi moduów instalowanych na dnie morza 93<br />

<br />

oraz zakadajc wspórzdn y L, 0 okrelon przez:<br />

y<br />

L,0<br />

( L)<br />

( L)<br />

rk<br />

gdzie: 0<br />

1 0 0,<br />

(L)<br />

r<br />

k<br />

- wektor wspórzdnych prowadnicy w ukadzie }<br />

( L)<br />

(L)<br />

k 1,..., , n - liczba prowadnic,<br />

<br />

<br />

n p<br />

p<br />

B<br />

(11)<br />

{L ,<br />

mona wyznaczy siy wynikajce z odksztace prowadnicy z zalenoci:<br />

gdzie:<br />

l)<br />

xc<br />

mxb<br />

xa<br />

<br />

xa<br />

l)<br />

zc<br />

mzb<br />

za<br />

<br />

za<br />

( ,<br />

( ,<br />

yL,0<br />

ya<br />

m ,<br />

y y<br />

b<br />

a<br />

( l)<br />

c,<br />

x,<br />

k<br />

F<br />

( l)<br />

c,<br />

z,<br />

k<br />

F<br />

( l)<br />

k,<br />

z<br />

<br />

<br />

L,0<br />

( l)<br />

c<br />

( l)<br />

c<br />

<br />

<br />

c x x<br />

(12.1)<br />

( l)<br />

k,<br />

z<br />

c z z<br />

(12.2)<br />

x L, 0, zL,0<br />

- wspórzdne otrzymane z zalenoci (11) zakadajc odpowiedni wersor<br />

w operatorze zwajcym ,<br />

( l)<br />

( l)<br />

c c - wspóczynniki sztywnoci.<br />

k, x,<br />

k,<br />

z<br />

<br />

W pracy zaoono, i siy hydrodynamiczne mog by opisane zgodnie z równaniem<br />

Morisona [8]. Zaoenie to jest suszne zarówno dla lin jak i samego adunku, jeli wymiar<br />

charakterystyczny jest co najmniej piciokrotnie mniejszy od dugoci fali [9].<br />

Równania ruchu caego ukadu, wraz z równaniami wizów, zapisano w formie ukadu<br />

równa róniczkowych:<br />

<br />

Aq<br />

DR Ht,<br />

q,<br />

q<br />

<br />

T<br />

(13)<br />

D q<br />

G t,<br />

q,<br />

q<br />

gdzie: A diag A (1)<br />

, A (2)<br />

, I (1)<br />

GW , I (2)<br />

GW<br />

, A ( l)<br />

~<br />

L, <br />

( l ) ~ ( l<br />

diag A<br />

)<br />

0<br />

,..., A<br />

<br />

T<br />

T<br />

(1) (2) (1)<br />

q q GW<br />

T<br />

T T<br />

( l)<br />

l<br />

l<br />

q<br />

~<br />

... q<br />

~ <br />

0<br />

n<br />

,<br />

l<br />

(2)<br />

GW<br />

T<br />

T<br />

L<br />

<br />

L,0<br />

A ,<br />

q q - wektor wspórzdnych uogólnionych,<br />

q<br />

<br />

<br />

H - wektor prawych stron równa ruchu,<br />

<br />

(2)<br />

T<br />

z<br />

(1) (1) (1) (2) (2)<br />

R Rx<br />

Ry<br />

Rz<br />

Rx<br />

Ry<br />

R - wektor nieznanych reakcji wizów,<br />

D, G - macierze wspóczynników równa wizów.<br />

<br />

Równania (13) cakowano metod Rungego-Kutty IV rzdu ze staym krokiem<br />

cakowania. Warunki pocztkowe zagadnienia dynamiki wyznaczano metod Newtona<br />

rozwizujc ukad równa nieliniowych w zadaniu statyki.<br />

<br />

nl


94 Marek Szczotka<br />

3. PRZYADOWE SYMULACJE NUMERYCZNE<br />

W oparciu o wyprowadzony model matematyczny ukadu zbudowano odpowiedni<br />

program komputerowy wykorzystujc rodowisko programistyczne Microsoft Visual C++.<br />

Podstawowe dane zaoone do oblicze przedstawiono w Tabeli 1.<br />

Przyjte wymuszenie statku w kierunku pionowym i wzdunym (ruchy nurzania<br />

i oscylacji wzdunych) jest opisane prostymi funkcjami harmonicznymi, co odpowiada<br />

przypadkowi fali regularnej (Rys. 6 a). Kty obrotu statku przyjto jako równe zeru.<br />

Zaoono równie brak wymuszenia ruchu bbna wcigarki gównej (dugo liny<br />

odwinitej z bbna nie zmienia si).<br />

Na dynamik adunku oddziauj, oprócz wymusze spowodowanych ruchem statku<br />

i/lub nawijaniem liny na bben gównej wcigarki, siy hydrodynamiczne. Przyjto, e<br />

prdko czstek wody spowodowana pywami wynosi: v c na powierzchni akwenu<br />

i maleje liniowo do zera przy powierzchni dna. Dodatkowe zmiany prdkoci wody na<br />

danej gbokoci wynikaj z ruchu orbitalnego czsteczek, wywoanego falowaniem.<br />

Uwzgldniono take siy hydrodynamiczne dziaajce na liny prowadzce i non.<br />

Tablica. 1<br />

Podstawowe parametry modelu przyjte w obliczeniach<br />

Parametr Warto Parametr Warto<br />

( l)<br />

masa adunku m L<br />

4000kg nacig lin F , l 1, 2<br />

16.2kN<br />

l,<br />

nom<br />

gboko d 400m prdko prdu morza v c<br />

wysoko fali<br />

wspóczynnik oporu<br />

hydrodynamicznego adunku<br />

H s<br />

1.0m okres fali T z<br />

C D , L<br />

promie bbna wcigarki liny<br />

prowadzcej<br />

2.2<br />

wspóczynnik oporu<br />

hydrodynamicznego liny<br />

C D , l<br />

0.0 - 0.5 s<br />

m<br />

9 s<br />

1.4<br />

0.15m rednica nominalna lin prowadzcych 24mm<br />

a )<br />

b)<br />

przemieszczenia[m]<br />

vL[m/s]<br />

czas[s]<br />

czas[s]<br />

Rys. 6. a) Przebiegi wymuszenia ruchu jednostki pywajcej,<br />

b) prdko adunku wzdu osi y


Modelowanie i analiza ukadu do obsugi moduów instalowanych na dnie morza 95<br />

Dla adunku opuszczonego na gboko 300m obliczono przemieszczenia wywoane<br />

falowaniem i parciem prdu, a wyniki przedstawiono na Rys. 7. W przypadku braku<br />

systemu lin prowadzcych, odchylenie adunku w kierunku prdkoci prdu dochodzi do<br />

4m. Zastosowanie lin prowadzcych, nawet przy stosunkowo niskim nacigu, prowadzi do<br />

ograniczenia tego przemieszczenia. W analizowanym przykadzie odchylenie adunku<br />

wynosi ok. 0.5m. Przyjty tryb pracy wcigarek lin prowadzcych zapewnia wzgldnie<br />

stay nacig lin, dlatego w przypadku pojawienia si znacznych si bocznych (wywoanych<br />

na przykad silnym prdem), lina zostanie odwinita z bbna i adunek podlega<br />

zwikszonym przemieszczeniom poprzecznym. Po ustaniu si poprzecznych, stay moment<br />

napdowy powoduje nawijanie si lin, i adunek jest sprowadzany do pozycji ustalonej<br />

przez napicie lin.<br />

a )<br />

b)<br />

xL[m]<br />

Fliny[kN]<br />

czas[s]<br />

czas[s]<br />

Rys. 7. a) Wpyw lin prowadzcych na wspórzdn x adunku,<br />

b) siy nacigu liny prowadzcej przy bbnie oraz reakcja w miejscu mocowania liny w dnie<br />

Warto bezwzgldn reakcji pionowej obliczonej dla punktu zamocowania liny<br />

prowadzcej w dnie przedstawia rysunek 7 b). Bezwadno ukadu powoduje, e<br />

sterowanie prac wcigarki (zastosowano stay moment napdowy bez informacji o stanie<br />

ukadu) powinno uwzgldnia równie prdko zmiany pooenia statku. Niska warto<br />

si nacigu lin prowadzcych oraz znaczna masa adunku, powoduje due wartoci<br />

przemieszcze ktowych (i prdkoci) bbnów wcigarek, dlatego efekty dynamiczne<br />

ujawniaj si w postaci waha siy napicia lin.<br />

4. PODSUMOWANIE<br />

W pracy przedstawiono ukad systemu do obsugi moduów stosowanych w instalacjach<br />

offshore przy wydobyciu ropy i gazu. Uproszczony model matematyczny umoliwia<br />

wstpn analiz si powstajcych w trakcie pracy urzdzenia w warunkach falowania<br />

morskiego. Wyniki symulacji mog by przydatne w fazie projektowania urzdzenia,<br />

definiowania wymaga i specyfikacji poszczególnych komponentów ukadu, jak równie


96 Marek Szczotka<br />

do okrelenia warunków granicznych, przy których moliwa jest jeszcze bezpieczna praca<br />

urzdzenia. Ze wzgldu na znaczne koszty poszczególnych moduów instalowanych na<br />

dnie morza, operacje zwizane z ich przemieszczaniem musz by wykonane w sposób<br />

kontrolowany. Opracowany model i program wzgldnie szybko pozwalaj na ocen<br />

ustawie podstawowych parametrów eksploatacyjnych, geometrycznych i tych zwizanych<br />

z bezwadnoci, dla danego rodzaju wymuszenia (warunków pogodowych).<br />

Innym moliwym zastosowaniem przedstawionego modelu jest sterowanie prac<br />

wcigarek AHC oraz pomocniczych, do uzyskiwania danego nacigu lin prowadzcych.<br />

Moliwe jest poczenie modelu dynamiki z algorytmami sterowania, oraz czujnikami<br />

ruchu jednostki. Ukad taki zostay zbudowany oraz zainstalowany w prezentowanym<br />

w pracy urzdzeniu. Jako kompensacji waha adunku, przy masach rzdu 20-30 ton<br />

wynosi 95% przy charakterystycznej wysokoci fal 5.0m. Wynik ten cakowicie wystarcza<br />

na bezpieczn instalacj infrastruktury w warunkach bardzo intensywnego falowania.<br />

Autor dzikuje firmie AXTech AS (Molde, Norwegia), producentowi systemów MHS, za zgod na<br />

opublikowanie zdj i szereg danych do oblicze.<br />

Bibliografia<br />

1. Hann M.: Komputerowa analiza niezawodnoci i bezpieczestwa maszyn i konstrukcji okrtowych<br />

poddanych koysaniom. Okrtownictwo i egluga, 2001.<br />

2. Hann M.: Statics and dynamics of multi-cable systems for semi-submersibles. Marine Structures 1995,<br />

Vol. 8, 1, s. 555-583.<br />

3. Szelangiewicz T.: Wpyw kotwicznego systemu utrzymania pozycji na koysania statku w obecnoci<br />

wiatru prdu i falowania. Prace Naukowe Politechniki Szczeciskiej, Nr 523, 1995.<br />

4. Szczotka M.: Simulation of an AHC system during offshore installation. Logistyka 2010, 2, s. 2311-2320.<br />

5. Craig J.J.: Wprowadzenie do robotyki. WNT, Warszawa 1995.<br />

6. Wittbrodt E., Adamiec-Wójcik I., Wojciech S.: Dynamics of flexible multibody systems. The rigid finite<br />

element method. Springer 2006.<br />

7. Szczotka M.: Pipe laying simulation with an active reel drive. Ocean Eng. 2010, Vol. 37, 2, s. 539-548.<br />

8. Morison J.R., O'Brien M.P., Johnson J.W., Schaaf S.A.: The force exerted by surface waves on piles.<br />

Petroleum Transactions, 1950, Vol. 189, s. 149-154.<br />

9. Chakrabarti S.K.: Hydrodynamics of offshore structures. WIT Press 1987.<br />

MODELLING AND SIMULATION OF AN SUBSEA MODULE HANDLING SYSTEM<br />

Summary: Dedicated, multi-purpose module handling system (MHS) is presented in the paper. Typically<br />

such equipment is assembled on a special offshore vessel, performing various construction works related to<br />

subsea field developments. Special modules (for example pumps, compressors, special tanks) are lowered<br />

towards the seabed and need to be positioned precisely on defined foundation. This is possible thanks to the<br />

application of the guide wires and an AHC (Active Heave Compensation) system. The module handling<br />

system is useful tool in the case of strong sea currents. When the module has to be lowered directly into a<br />

precise location on the seabed, general, offshore crane cannot be used easily. In this case, the MHS is far<br />

better solution. A simplified mathematical model of the system is shortly presented in the article, together<br />

with some example simulation results. The influence of the guide wires on the behavior of a load is presented<br />

through the results.<br />

Keywords: MHPS, simulation of module installation, offshore equipment<br />

Recenzent: Lech Kobyliski

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!