|
Sa se proiecteze o conducta pentru transportul gazelor naturale, cu urmatoarele date initiale, tinand cont de numarul de ordine j =15.
Debitul de gaze
Qg j [Stm³/zi] ;
Qg [Stm³/zi]
Lungimea conductei
L = 120+ 0,1 · j [km] ;
L [km]
Presiunea finala
p2 = 6 + 0,1·j [bar] ;
p2 = 6+0.1·15 =7,5[bar] ;
Diametrul conductei
D = 12 in ;
D = 304,8 [mm] = 0,3048 [m]
Densitatea relativa a gazelor
= 0,554 ;
Vascozitatea relativa a gazelor
μ = 0,012 cP ; μ = 0,012 ·10 ‾ 3 [N/m2]
Densitatea aerului
ρaer = 1,293 [kg/m3] ;
Rugozitatea conductei
ks = 0,03 [cm] ;
Factorul de neidealitate
Z = 0,9 + 0,001·j
Z = 0,9 + 0,001·15 = 0,915;
Calculam Z cu relatia lui Istomin:
Recalculam pe Z pentru CH4 cu relatia lui Dranchuk si Abou-Kassem:
in care:
ρr
Tr = temperature pseudoredusa;
Pr = presiunea pseudoredusa;
A1 = 0,3265; A2 = -1,0700; A3 = -0,5339; A4 = 0,01569; A5 = -0,05165; A6 = 0,5475;
A7 = -0,7361; A8 = 0,1844; A9 = 0,1056; A10 = 0,6134; A11 = 0,7210.
Tcr = 190 [K]
Tabs= 15°C
Pcr = 4,54[MPa]
Pabs= 7,9 [bar]
=
Z=0,9673
Pentru stabilirea regimul de curgere trebuie calculat numarul lui Reynolds.
Luand drept criteriu numarul lui Reynolds, in functie de marimea acestuia se pot distinge 3 regimuri fundamentale de curgere: laminar, instabil si turbulent.
Calculul debitului
Q = [Stm3/s]
Q = = 1,331 [Stm3/s]
Calculul vitezei medii
vm = [m/s]
vm = [m/s]
Determinarea densitatii gazelor transportate
ρg = δ · ρaer [kg/m3]
ρg [kg/m3]
Calculul numarului lui Reynolds
Re ==331719.936 > 2300
In urma calcularii numarului Reynolds se observa ca regimul de curgere este turbulent.
I.2 Calculul coeficientului de rezistenta
D = [m]
ks = [m]
Verificarea cu formula lui Weymouth:
v
v
RcrI
RecrII
RecrII = 1389301,04
RecrII > Re> RecrI
Regim turbulent mixt.
I.3 Calculul presiunii de intrare in conducta - relatia lui Weymouth
p1 = ;
unde: Q - debitul [Stm3/h] ;
L - lungimea conductei [km] ;
D - diametrul conductei [cm] ;
Z - factor de neidealitate; - presiune [bar].
p1 = = 8,2956[bar]
Pentru S.I.:
p1=
p1 =7.5
Conform Normativului departamental 3915-94 pentru proiectarea si construirea conductelor colectoare si de transport gaze naturale emis in anul 1995.
p1 =
P0 = Pa→P0=1bar=105 Pa
T0 = K→T0=273,15 K
L = m→L = 121,5 Km=121,5·103 m
T = K →T =288,15 K
Ra=J/KgK→Ra=278,041 J/KgK
Q = Stm3/h→Q= Stm3/h
D= cm→ D=30,48 cm
P1 = 7,5 bar
Pentru Z 0,9673→P1 7,5 bar
p1 =
T = 288.15 K
L = 121,5 ·103 m
Q= Stm3/h
D=30.48 cm
p1 =
p1 = 750008,63 Pa.
Pentru determinarea coeficientului λ normativului departamental 3915-94 face trimitere la formula Colebrook- White.
Q- Stm3/s
μ - N/m
D ,48 cm = 304,8 mm
I.4 Calculul pantei hidraulice
i = = 0,006548; ;
I.5 Calculul presiunii medii
pm = [bar]
I.6 Calculul caderii de presiune pe tronsoane
xi=1
Impartim conducta in zece bucati adica L/10:
px =
px1 = = 8.2195 [bar]
px2 = = 8.1427 [bar]
px3 = = 8.0584 [bar]
px4 = = 7.9868[bar]
px5 = =7.9078 [bar]
px6= =7.8279[bar]
px7 ==7.7472[bar]
px8 = =7.6657 [bar]
px9 = =7.5833[bar]
px10 = =7,5 [bar]
Tronsonul xi [km] |
Presiunea Pxi [bar] |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
I.7 Calculul volumului de gaze stocat in conducta
Conducta se poate folosi si pentru stocarea temporara a gazelor, ea comportandu-se ca un rezervor orizontal cu un diametru mai redus.
[SI]
unde: D - diametrul interior al conductei [m] ;
V - cantitatea de gaze stocata in conducta [Stm3] ;
p0 - presiunea normala ; p0 = 101325 Pa ;
T0 - temperatura normala ; T0 = 273,15 [K] ;
T - tempertura gazelor din conducta ; T = 15 °C ;
p1,2 - presiunea initiala, respectiv finala ;
L - lungimea conductei [m] ;
= 64282.46727 [Stm3]
[MKFS]
II.1 Calculul grosimii peretelui tevii conductei de transport gaze
unde:
e-grosimea nominala a peretelui
Dext-diametrul exterior al conductei;
Dext 323,8mm;
Pr max presiunea maxima de regim a gazelor [bar]
φ - coeficientul de calitate ale imbinarilor sudate;
σc - rezistenta minima la curgere a otelului materialului tubular [N/m2];
σc 3870 N/m2;
S - coeficient de siguranta care se ia functie de clasa de locatie a traseului conductei;
S
a - adaos la grosimea de perete a conductei pentru coroziune,abraziune,toleranta de fabricatie;
a = a1+a2+a3
a = 0.
Se alege din STAS A.P.I. specificatia 5L:
e = 7.6mm
II.2 Calculul cantitatii de metal necesar pentru o conducta de transport gaze
G = Vmet·103·2,467 [t/km] ;
unde: G - cantitatea de metal [t/km] ;
Dint - diametrul interior al conductei [cm] ;
Dext- diametrul exterior al conductei [cm] ;
G [t/km] .
II.3 Calculul temperaturii in conducta de transport gaze
Spre deosebire de gazele ideale, gazele reale se pot lichefia prin comprimare si/sau racire. Lichefierea are loc numai la anumite presiuni si temperaturi, specifice pentru fiecare gaz in parte.
Temperatura la care apare prima picatura de condens se numeste temperatura punctului de roua.
;
unde: Tx - temperatura in conducta de transport gaze la distanta x de punctul initial al conductei [K];
Ts - temperatura solului [K];
T0 - temperatura intiala a gazelor [k]; T0 = 288,15[k];
k - coeficientul de transfer al caldurii [kcal/m2hK]; k = 1,3 [kcal/m2hk];
Dext - diametrul exterior al conductei [m]; Dext = 0,3238 m
Q - debitul de gaze [stm3/h] ;Q = 4183[stm3/h]
- densitatea gazelor [kg/m3] ; ,7112[kg/m3]
cp - caldura specifica la p0 [kcal/kgk] ;
cp = 0,36 + 0,599 10-3T [kcal/kgk] ;
cp 10-3 288,15 = 0,533 [kcal/kgk] ;
Determinarea distantei corespunzatoare aparitiei punctului de roua:
unde x este distanta de la punctul initial pana la punctul in care se formeaza punctul de roua.
- temperatura punctului de roua:Tx = 285,15 [k]
temperatura solului pentru iarna: Ts = 275,15 [k] ;
temperatura solului pentru vara: Ts = 284,15 [k] ;
Pentru iarna :
Pentru vara :
Calculul debitului:
debitul de gaze [Stm3/s]
T0 - temperatura normala ; T0 = 273,15 [K] ;
p1,2 - presiunea initiala, respectiv finala [Pa]
λ - coeficient de frecare hidraulica;
L = 121,5 Km = 121,5·103 m;
D = 0,3048 m;
Z = 0,9673;
Ra 278,041[J/kgK].
II.4 Calculul conductelor de gaze cu ramificatii care primesc sau livreaza gaze pe traseu
Conform schemei de mai jos, sa se evalueze care ar fi presiunile in punctele: B, C, D, E ale conductei de transport gaze daca:
[m]
q1 =
q2 = [m3·h-1]
q3 = [m3·h-1]
QA=Q QA=4791.66 [m3 ora
QB=Q - q1 QB=2395.83 [m3 ora
QC=QB - q2 QC=1197.9134 [m3 ora
QD=QB + q3 QD=5989.58 [m3 ora
QE=QD QE=5989.58 [m3 ora
Q = 1,331 [Stm3/s]=1,331·3600 [Stm3/h]
L = 121,5 Km = 121,5·103 m
D = 0,3048 m = 30,48 cm
Z = 0,9673
[SI]
Pentru tronsonul L4
Pentru tronsonul L3
Pentru tronsonul L2:
Pentru tronsonul L1
Pentru sistemul tolerant de unitatii de masura
Q = 1,331 [Stm3/s]=1,331·3600 [Stm3/h]
L = 121,5 Km
D = 0,3048 m = 30,48 cm
Z = 0,9673
Pentru tronsonul L4
Pentru tronsonul L3
Pentru tronsonul L2:
Pentru tronsonul L1:
III. Comprimarea gazelor
Presiunea necesara penru transportul gazelor naturale prin conducte se realizeaza in statii de comprimare.
Cand presiunea gazelor nu mai este suficienta pentru a asigura debitul prevazut se instaleaza o statie de compresoare in punctul initial al conductei si o serie de statii intermediare pe traseul conductei.
Procesul de comprimare al gazelor se realizeaza cu ajutorul compresoarelor care pot fi de mai multe tipuri. Se instaleaza o statie de comprimare a gazelor cu piston in punctul final al conductei.
Se cunosc:
Coeficientul transformarii politrope: n= 1,28
Spatiul mort: m = 0,15
Raportul de comprimare: r = 2,3
Cursa pistonului: s = 0,38[m]
Diametrul pistonului: DC = 0,15[m]
Turatia compresorului: nC = 310[rot/min]
Diametrul tijei: dt = 0,025[m]
Numarul de cilindrii: NC
Volumul de aspiratie: VA = 1[m3]
Coeficientul transformarii adiabate: k
III.1 Calculul lucrului mecanic
III.1.1 Lucrul mecanic in destinderea izoterma
Liz = paValn
Unde: Liz - lucrul mecanic in destinderea izoterma;
pa - presiunea de aspiratie:
III.1.2 Lucrul mecanic in destinderea adiabata
Lad =
III.1.3 Lucrul mecanic in destinderea politropa
Lpol =
III.2 Randamentul volumetric al compresorului
Raportul dintre volumul de gaze aspirat Va si plasamentul pistonului Vc reprezinta randamentul volumetric al pistonului.
III.2.1 Randamentul compresorului in detenta izoterma
ηvizt = 1 - m (r - 1) = 1 - 0,15 (2,3 -1) = 0,805 = 80,5 % ;
III.2.2 Randamentul compresorului in detenta adiabata
ηvad = 1 - m ( - 1) = 1 - 0,15 ( -1) = 0,859 =85,9 % ;
III.2.3 Randamentul compresorului in detenta politropa
ηvpol = 1-m ( - 1) = 1 - 0,15 ( -1) = 0,862 = 86,2 % ;
III.2.4 Randamentul efectiv al compresorului
include mai multe pierderi:
Pierderi cauzate de neetanseitatea supapelor, segmentilor
Se alege η1 =0,97
Pierderi la aspiratie si refulare
Se alege η2 =0,97
Pierderi de debit datorita pierderii gazului aspirat
Se alege η3=0,98
Pierderi de gaze datorita umiditatii gazului aspirat
sau
ηe = η1· η2· η3· η4
ηe = 1-m ( - 1) - 0,01r = 1-0,15 ( -1) -0,01
sau:
ηe = 0,97-m ( - 1) 0.97-0,15 ( -1)=0.832=83.2%
III.3 Temperatura gazelor la intarea in compresor
unde :
T0 este temperatura de referinta T0 = 288,15 [ K ] ;
Ts - temperatura solului [K];
k - coeficientul de transfer al caldurii [kcal/m2hK]; k = 1,26 [kcal/m2hk];
Dext - diametrul exterior al conductei [m]; Dext = 0,3238 m
Q - debitul de gaze [stm3/h] ;Q = 1,331[stm3/s]
cp - caldura specifica la p0 [kcal/kgk] ; cp= [kcal/kgk];
III.3.1 Temperatura gazelor la intarea in compresor vara
Temperatura solului: Ts=284,15[K]
III.3.2 Temperatura gazelor la intarea in compresor iarna
Temperatura solului: Ts=275,15[K]
Calculul randamentului efectiv al compresorului calculat cu:
Sau:
III.4 Debitul masic al compresorului
G=Q g [kg/s] ;
III.5.1 Puterea pentru comprimarea izoterma
Niz = Liz ·G = 624681.8422[J/kg] · 0,9533[kg/s] = 595509.2002[W] = 595.5092 [kW]
III.5.2 Puterea pentru comprimarea adiabata
Nad = Lad ·G = 681399.2947 [J/kg] · 0,9533[kg/s] = 649577.9476[W] =649.5779 [kW]
III.5.3 Puterea pentru comprimarea politropa
Npol = Lpol ·G = 685209.5508[J/kg] · 0,9533[kg/s] = 653210.2648 [W] = 653.2102[kW]
III.6 Calculul debitul compresorului monoetajat
cu simplu efect
vara
iarna
cu dublu efect
vara
iarna
III.7 Calculul debitul total al compresorului monoetajat
cu simplu efect
vara
q = q'·Nc=14.438·5=72.19
iarna
q = q'·Nc=14.8987·5=74.4935
cu dublu efect
vara
q = qd·Nc=28.9451·5=144.7255
iarna
q = qd·Nc=29.7974·5=148.987
III.8. Calculul puterii compresorului monoetajat
III.8.1 Calculul puterii compresorului monoetajat cu simplu efect
1. Neluand in calcul influenta factorului Z
MKfs
Vc - volumul geometric al cilindrilor descris de piston
Daca se considera ca acest compresor are 5 cilindri atunci:
V5 = 5·Vc = 5·2,08 = 10,4 [m3/min]
Pa = 7,5·1,02 = 7.65 at
Pa = 7.65 +1 = 7.75 ata
N = 193.665 [cP]
N=193.665·735,5= 142440.6281 [W]
S.I.
Luand in calcul influenta factorului Z
vara:
iarna:
III.8.2 Calculul puterii compresorului monoetajat cu dublu efect
1)Luand in considerare influenta factorului Z:
MKfs:
S.I.:
N= 429173.942[W]
2)Neluand in considerare influenta factorului Z:
Vara:
Iarna:
IV. Calculul de racire al compresorului
Vom considera ca transformarea gazelor in compresor este politropa si se realizeaza cu un compresor cu piston.
Sa se calculeze cantitatea de caldura ce trebuie evacuata din cilindrul compresorului si cantitatea de apa necesara pentru racirea compresorului.
Temperatura gazului la intarea in compresor T1=Tr
Temperatura gazului la iesirea din compresor T2
IV.1. Calculul de racire al compresorului cu piston
IV.1.1. Cantitatea de caldura ce se degaja la racirea cilindrilor
Vara:
Qc=96319.0233 [KJ/h]
Iarna:
Qc=92060.7312 [KJ/h]
IV.1.2. Cantitatea de caldura ce trebuie evacuata din cilindru intermediar
Vara:
Qi=1765501.049 [KJ/h]
Iarna:
Qi=1687447.733 [KJ/h]
IV.1.3. Cantitatea totala de caldura ce trebuie evacuata
Unde x este nr de trepte al compresorului.
1)Pentru compresorul cu piston intro singura treapta:
Vara:
=96319.0233 [KJ/h]
Iarna:
=92060.7312 [KJ/h]
2)Pentru compresorul cu piston cu doua treapte:
Vara:
[KJ/h]
Iarna:
[KJ/h]
IV.1.4. Cantitatea de apa necesara pentru racirea compresorului
1)Pentru compresorul cu piston intro singura treapta:
Vara:
Iarna:
2)Pentru compresorul cu piston cu doua trepte:
Vara:
Iarna:
IV.1.5. Cantitatea de apa necesara pentru racirea motorului
Pentru simplu efect:
Pentru dublu efect:
Racirea se realizeza numai vara deoarece iarna temperatura este mentinuta sub 60 C si motorul nu are nevoie de racire suplimentara.
IV.2. Calculul numarului de agregate separate pentru fiecare treapta
C
c-coeficient.
IV.3. Calculul puterii absorbite de un compresor centrifugal
Puterea absorbita de un compresor:
-cresterea entalpiei gazelor
Vara:
Iarna:
IV.4. Calculul puterii absorbite totala
Vara:
Iarna:
IV.5. Calculul randamentului termic al turbocompresorului
Hi=8500 KJ/kg;
Vara:
Iarna:
Continutul de apa in gazele naturale este in functie de presiunea si temperatura lor, parametrii care variaza pe parcurs. Acest continut este maxim in zacamant, apoi o parte din vaporii de apa condenseaza si se separa, iar gazul ramane saturat.
Extragerea vaporilor de apa din gazele naturale se poate realiza prin:
comprimare, racire si separare;
absorbtie cu o solutie higroscopica (glicoli);
adsorbtie cu un desicant solid (silicagel).
Deshidratarea prin adsorbtie:
In procesul de adsorbtie retinerea vaporilor de apa din gazele naturale se face prin concentrarea lor pe suprafata adsorbantului care este foarte mare pe unitatea de greutate.
Cand suprafetele s-au saturat cu vapori de apa si se ajunge la un echilibru, desicantul trebuie regenerat, ceea ce se realizeaza prin scaderea presiunii gazului care trece prin el sau prin cresterea temperaturii. Acest procedeu este raspandit in toate instalatiile industriale.
Utilizarea desicantilor solizi ofera urmatoarele avantaje:
o se pot obtine depresiuni mari ale punctului de roua;
o se pot exploata cu parametrii variabili;
o simplicitate in proiectare si exploatare;
o nu pun probleme de coroziune, spumare etc.;
o se pot realiza instalatii pentru debite foarte reduse.
Principalele dezavantaje constau in suscebilitatea de contaminare si pierderea capacitatii de adsorbtie, consumul de caldura mai ridicat decat in cazul instalatiilor de absorbtie si in sfarsit, caderile de presiune mai mari la trecerea gazelor prin masa adsorbantului, decat prin coloana cu talere.
Se considera o instalatie de deshidratare prin adsorbtie cu site moleculare avand urmatoarele caracteristici:
Debitul masic de gaze: Q = 1,15 105 [m3/zi]
Temperatura de intrare a gazelor: Tgi = 288,15 [K]
Punctul de roua: Tr = 285,15 [K]
Presiunea in coloana de adsorbtie: pab = 80 [bar]
Durata unui ciclu: tad = 6 [ore]
Capacitatea reala de adsorbtie: a = 18 %
Continutul de umiditate la intrare: u1 = 1,7 [g/m3N]
Continutul de umiditate la iesire: u2 = 0,055 [g/m3N]
Diametrul coloanei de adsorbtie: Dad = 1,2 [m]
Densitatea sitelor moleculare : [kg/m3]
V.1 Debitul maxim
V.2 Volumul ocupat de adsorbant in coloana
V.3 Inaltimea coloanei de adsorbtie
CONCLUZII
Proiectarea unui sistem de transport gaze si de comprimare a gazelor este un calcul complex care are in vedere realizarea unei scheme tehnice de transport astfel incat cheltuielile efectuate pentru acestea sa fie cat mai reduse, dar cu un randament maxim.
Unii dintre parametri care intervin in calcul, depind de proprietatile gazului transportat, deci proiectarea sistemului de transport tine seama de acesti parametri, deci de respectivele proprietati(densitate, vascozitate, presiune, temperatura etc.) si este valabila numai pentru gazul transportat pentru care s-a facut calculul.
Tratarea gazelor pentru transport se face prin procesul de adsorbtie a apei din gaze pe site moleculare. Fenomenul de adsorbtie se bazeaza pe proprietatea pe care o au unele corpuri solide de a concentra pe suprafata lor anumite substante in stare de vapori sau lichid. Tratarea gazelor se face pentru a preveni formarea criohidratilor care blocheaza conducta.
Oroveanu, T., David, V., Stan, Al., Trifan, C., - Colectarea, transportul, distributia si depozitarea produselor petroliere si gazelor - Editura Didactica, Pedagogica, Bucuresti, 1983;
Soare, Al., - Transportul si depozitarea fluidelor - vol.II - Editura Unviresitatii din Ploiesti
Puscoiu, N. - Extractia gazelor naturale. Aplicatii de calcul. Editura Tehnica Bucuresti 1989
Grigoras, D. - Note de curs
*** - STAS 715/2-82 ;