Akademia Górniczo-Hutnicza
Im. Stanisława Staszica w Krakowie
Wydział Wiertnictwa Nafty i Gazu
Obliczanie zasobów gazu ziemnego
metodą objętościową i spadku ciśnienia
dla złoża F25
Budowa geologiczna:
Złoże typu warstwowego w asymetrycznej antyklinie o nachyleniu skrzydła SW 60˚, skrzydła NE 25˚. Poziom zbiornikowy stanowią piaskowce średnioziarniste wieku Pg2, podścielone czerwonymi łupkami tego samego wieku. Profil antykliny(ponad poziomem zbiornikowym): 100 m czarnych łupków, 200 m drobnorytmicznego fliszu łupkowo-piaskowcowego, 300 m łupków menelitowych - Pg3, 200 m gruboławicowych piaskowców krośnieńskich, ok. 400 m piaskowcowo-łupkowych warstw krośnieńskich(Pg3). Jądro antykliny na powierzchni budują łupki menelitowe, występujące w dolinie potoku(300 m n.p.m.), a wzgórza po obu stronach doliny(ok. 450-500 m n.p.m.) tworzą piaskowce warstw krośnieńskich.
Obliczanie zasobów gazu ziemnego metodą objętościową dla złoża F25
Qg= F*hef*Φef(1-Sw)/Bg
F - powierzchnia złoża=1 350 000 m2
hef - średnia efektywna miąższość poziomu gazonośnego lub też strefy
nasycenia=40m
Φef - średnia porowatość efektywna poziomu gazonośnego=10%
Sw - nasycenie wodą związaną dla K=330 mD
Log330=2,52
Sw=23,6
Bg - parametr termobaryczny
Bg= P0/ T0+ ( Pzł * z)/ Tzł
P0 - ciśnienie normalne, czyli standardowe = 1 atm
T0 - temperatura normalna, czyli standardowa = 14,5˚C=287,5 K
Tzł - temperatura złożowa, czyli temperatura panująca w złożu, zależy ona od
obszaru=55,025˚C=328,025 K
Hk - głębokość konturu złoża, liczona od powierzchni terenu=1425 m
Pzł - ciśnienie złożowa, czyli ciśnienie panująca w złożu = 15377175 Pa = =151,7727173 atm
γw=1,1g/cm3
z - współczynnik ściśliwości gazu ziemnego zależny od składu chemicznego gazu
z= zdw*ydw+(1-ydw)*zg
zg - współczynnik ściśliwości węglowodorów
Składnik Gazu |
Ciężar Cząsteczk. |
Konc. mol. |
Ciężar mol |
Tkr [K] |
Pkr [ata] |
konc.mol. x Tkr |
konc.mol x Pkr |
CH4 |
16.043 |
0.89345 |
1.,334 |
191.0 |
46.0 |
170.649 |
41.099 |
C2H6 |
30.070 |
0.02525 |
0.759 |
305.0 |
50.0 |
7.701 |
1.263 |
C3H8 |
44.097 |
0.01946 |
0.858 |
369.0 |
42.5 |
7.181 |
0.827 |
nC4H10 |
58.124 |
0.01734 |
1.009 |
426.0 |
37.0 |
7.387 |
0.642 |
iC4H10 |
58.124 |
0.01282 |
0.745 |
407.0 |
38.0 |
5.218 |
0.487 |
nC5H12 |
72.151 |
0.00845 |
0.610 |
470.0 |
33.0 |
3.972 |
0.279 |
iC25H12 |
72.151 |
0.00247 |
0.178 |
461.0 |
34.0 |
1.139 |
0.084 |
C6H14 |
86.178 |
0.00299 |
0.258 |
508.0 |
29.5 |
1.519 |
0.088 |
CO2 |
44.010 |
0.01777 |
0.782 |
304.3 |
75.2 |
- |
- |
Σ=19.533 Σ=204.466 Σ=44.769
Tpkr= konc.mol.xTnkr
Ppkr= konc.mol.xPnkr
Tr= Tzł / Tpkr =1,60
Pr= Pzł / Ppkr =3,39
zg=0,82
zdw - współczynnik ściśliwości CO2 w warunkach Pzł i Tz=0,4
ydw - zawartość CO2 w gazie=0.01777
z= 0,4*0.01777+(1-0.01777)*0,82=0,811366125
Bg=0,00610
Qg= 1 350 000 *40*Φef(1-0,236)/ 0,00610=695 mln m3
Obliczanie zasobów gazu ziemnego metodą spadku ciśnienia dla złoża F25
Moment eksploatacji |
Pds |
Pgł |
zn |
Pds / zn |
Zasoby gazu (wydobyte) |
Zasoby pozostałe w złożu - Qg |
*1t0 |
15.38 |
15.37 |
z1=0,811 |
18.964 |
Qg= - |
695 |
*2t50% |
7.60 |
6.71 |
z2=0,811 |
8.274 |
Q50%=342.5 |
342.5 |
*3tk1 |
0.85 |
0,75 |
z3=1 |
0.85 |
Qw=655 |
40 |
z3 obliczamy analogicznie do z1 wykorzystując do obliczeń ciśnienie Pds=0.85 więc:
Pr= Pds / Ppkr =0,02
zg=1
zdw=1
z3= 1*0.01777+(1-0.01777)*1=1
1