PENGGUNAAN DATA PLT DAN DATA PRODUKSI DALAM MENENTUKAN POLA ALIRAN FLUIDA PADA SUMUR X
Tugas Akhir Oleh: BRAVO MAHENDRA 12206050
Diajukan sebagai syarat untuk memperoleh gelar SARJANA TEKNIK PERMINYAKAN Institut Teknologi Bandung
PROGRAM STUDI TEKNIK PERMINYAKAN FAKULTAS TEKNOLOGI PERTAMBANGAN DAN PERMINYAKAN INSTITUT TEKNOLOGI BANDUNG 2011
PENGGUNAAN DATA PLT DAN DATA PRODUKSI DALAM MENENTUKAN POLA ALIRAN FLUIDA PADA SUMUR X
Tugas Akhir Oleh: BRAVO MAHENDRA 12206050
Diajukan sebagai syarat untuk memperoleh gelar SARJANA TEKNIK PERMINYAKAN Institut Teknologi Bandung
Disetujui oleh: Dosen Pembimbing,
Penulis,
Ir. Hernansjah
Bravo Mahendra, ST
NIP. 131695140
NIM. 12206050
PENGGUNAAN DATA PLT DAN DATA PRODUKSI DALAM MENENTUKAN POLA ALIRAN FLUIDA PADA SUMUR X Oleh: Bravo* dan Hernansjah**, Teknik Perminyakan – Institut Teknologi Bandung, 2010. Sari Salah satu kegunaan Production Logging Tool (PLT) adalah untuk mengetahui profil aliran fluida dalam sumur. Asumsi pola aliran yang salah dapat mengakibatkan interpretasi profil aliran menjadi kurang akurat dan berujung pada pengambilan keputusan workover yang salah pula. Karena profil aliran dipengaruhi oleh kecepatan alir, holdup, kecepatan slip dan pola aliran, maka diperlukan data PLT dan data produksi yang dapat digunakan untuk menentukan hal-hal yang mempengaruhi profil aliran yg terjadi dalam lubang sumur tersebut. Untuk itu perlu dilakukan beberapa kali rangkaian tes PLT pada berbagai kondisi debit aliran sumur. Kata kunci: profil aliran, pola aliran, holdup, data produksi, logging produksi, kecepatan slip Abstract Production Logging Tool (PLT) has been used to identify flow profile in the well. Flow profile can be altered by velocity, holdup, slippage velocity and flow regime. Wrong taken assumption on flow regime can cause flow profile interpretation goes inaccurate and would lead to a wrong workover decision. PLT response and production data can be used to determine the flow regime inside the wellbore. Production data can also be used to determine the slippage velocity with additional data of holdup data from PLT. By knowing the precise flow regime and slippage velocity, one can deliver a more accurate flow profile interpretation. Keyword: flow profile, flow regime, holdup, production data, production logging, slippage velocity
Pendahuluan Fluida yang mengalir dalam lubang sumur seringkali terdiri dari dua fasa atau lebih. Perbedaan densitas antar fasa tersebut akan membentuk suatu pola aliran tertentu. Pemahaman pola aliran fluida multifasa dalam lubang sumur yang baik akan menghasilkan analisa performa reservoir dan pengambilan keputusan workover yang tepat. Production Logging Tool (PLT, rangkaian peralatan logging produksi) dapat digunakan untuk mengetahui karakteristik aliran fluida multifasa dalam sumur, seperti: 1. Total kecepatan alir fluida, oleh spinner atau flowmeter, dan 2. Besarnya fraksi tiap fasa fluida yang mengalir menggunakan densitometer ataupun water-holdup meter. Tekanan dan temperatur pun bisa didapatkan dengan menambahkan pressure dan temperature gauge ke dalam rangkaian PLT. Respon PLT yang sama namun memiliki pola aliran yang berbeda dapat menghasilkan hasil interpretasi yang berbeda pula8. Pola aliran tersebut dapat ditentukan dengan bantuan data produksi (debit aliran tiap fasa yang terekam di permukaan) ketika proses pengambilan data PLT dilakukan diturunkan. Dengan mengetahui pola aliran yang tepat, hasil analisa PLT yang akurat pun dapat diraih.
*) Mahasiswa Teknik Perminyakan – Institut Teknologi Bandung **) Dosen Teknik Perminyakan – Institut Teknologi Bandung
Teori Dasar Fraksi atau Holdup Holdup atau fraksi berarti besarnya volume pipa yang ditempati oleh suatu fasa fluida8 (digunakan sebagai prinsip densitometer). Dengan asumsi pola aliran yang seragam sepanjang pipa yang diukur, holdup dapat berarti besarnya bagian yang ditempati oleh suatu fasa fluida dalam suatu luas penampang pipa3 (digunakan sebagai prinsip water-holdup meter, seperti FloView#). Pada aliran fluida fasa gas dan air10: (1) (2)
Gambar 1. Pola Aliran Mist
#) hak cipta Schlumberger
Reynolds Number yang dihasilkan oleh aliran fluida tersebut. (4) Hendy (2010)4 telah menyusun sebuah chart yang memperlihatkan hubungan antara VPCF dengan Reynolds Number. Gambar 2. Pola Aliran Annular-Mist
Terlihat dari Gambar 1 dan 2, dengan nilai fraksi yang sama dapat dihasilkan pola aliran yang berbeda. Holdup selalu lebih besar atau sama (pada aliran mist) dengan rasio fasa fluida di permukaan (GOR dan Water-Cut). Hal ini disebabkan oleh adanya perbedaan kecepatan aliran tiap fasa9. Gambar 4. NRe vs VPCF4
Kecepatan Alir Fluida Dari kalibrasi log spinner1, dapat dihasilkan kecepatan alir fluida total sesuai dengan hubungan: (3)
b
Kecepatan alir fluida total merupakan jumlah seluruh kecepatan superfisial tiap fasanya6. Pada aliran fluida fasa gas dan air: (5) Kecepatan superfisial bukanlah kecepatan alir aktual dalam aliran multifasa, karena kecepatan yang terekam tersebut merupakan kecepatan alir rata-rata akibat fasa fluida tidak sepenuhnya menempati luas penampang pipa. Pada kenyataannya, fasa fluida mengalir lebih cepat dari kecepatan superfisial, dan diartikan sebagai kecepatan fasa fluida ketika bernilai fraksi satu6. Pada aliran fluida fasa gas dan air:
Gambar 3. Flow Response Line1
(6) RPS o adalah respon log spinner ketika peralatan PLT berada dalam kondisi statis. b adalah gradien respon log spinner yang dihasilkan dari plot data RPS dengan kecepatan kabel. V t adalah kecepatan minimal yang dibutuhkan oleh spinner untuk mulai berputar. Nilai b dan V t dipengaruhi oleh viskositas fluida yang dilaluinya8. Karena viskositas fluida dipengaruhi oleh kecepatan alir dan komposisi campuran fluida, Flow Response Line yang paling akurat akan dihasilkan ketika spinner melalui fluida yang terdiri dari satu fasa dan berada dalam kondisi statis (sumur shut-in). Nilai kecepatan alir fluida yang dihasilkan dari Flow Response Line merupakan nilai kecepatan alir fluida pada bagian tengah penampang10. Untuk mendapatkan nilai kecepatan alir fluida rata-rata, maka nilai kecepatan alir fluida maksimum tersebut harus dikalikan dengan VPCF (Velocity Performance Correction Factor)4. VPCF dipengaruhi oleh besarnya nilai
Maka, hubungan antara debit aliran fluida dengan fraksi dan kecepatan aktual fluida tersebut adalah9: (7)
Gambar 5. Ilustrasi Sederhana Aliran Dua Fasa9
Dari persamaan 7 pun ditunjukkan hubungan antara debit aliran fluida yang terekam di dalam sumur dengan debit aliran fluida yang terekam di permukaan (adanya formation volume factor, dipengaruhi oleh tekanan dan temperatur). Kecepatan Slip Pada aliran multifasa, fasa fluida dengan densitas yang lebih ringan akan mengalir lebih cepat. Perbedaan kecepatan fasa fluida ini disebut dengan kecepatan slip. Pada aliran fasa gas dan air: (8) Dengan mengkombinasikan persamaan 3, 5, 6 dan 8 (penurunan pada lampiran 1) didapatkan: (9) Kecepatan slip dapat ditentukan menggunakan korelasi yang dikembangkan oleh Smollen (1996)10: • Untuk ρ g < 0.5 gr/cc dan Y w < 0.35: (aliran mist) •
Untuk ρ g < 0.5 gr/cc dan Y w > 0.35:
(10)
Gambar 6. Ilustrasi Pola Aliran Fluida
2
Studi Kasus Sumur X adalah sumur penghasil gas bertipe deviated-wellbore, dengan rata-rata kemiringan sumur sebesar 19 derajat. Interval perforasi berada pada kedalaman 3533-3534m, 35783587m dan 3841-3843m. Pada studi kasus kali ini, sumur X telah mengalami tiga kali rangkaian tes PLT, masing-masing pada kondisi debit aliran maksimum (Qmax, 9,73 MMSCF/day), setengah dari maksimum (Q1/2max, 4,87 MMSCF/day), dan saat sumur tidak mengalir (shut-in, 0 MMSCF/day). Untuk tiap rangkaian, dilakukan pula pembacaan dengan kondisi peralatan statis pada stasiun 1-3, masing-masing pada kedalaman 3515, 3560 dan 3860 m.
Metode lain untuk menentukan kecepatan slip adalah dengan menggunakan persamaan 7 dan 8, di mana debit fasa fluida dalam lubang sumur didapat dengan mengkonversikan data produksi ke dalam kondisi kedalaman (menggunakan data temperatur dan tekanan pada kedalaman tersebut). Pola Aliran Fluida Ketika dua fasa fluida mengalir dalam suatu pipa vertikal, maka perbedaan kecepatan superfisial tiap fasanya akan menghasilkan jenis pola aliran yang berbeda-beda.2 Dengan melakukan analisa terhadap respon PLT, pola aliran fluida multifasa dalam sumur dapat diketahui.
Gambar 7. Diagram Konfigurasi Sumur X
Interpretasi Pola Aliran Diambil data PLT statis untuk sumur X pada stasiun paling atas ketika debit aliran maksimum (Qmax).
Tabel 1. Data PLT Stasiun Teratas (Qmax) Y-w = 0,19022 RPSo = 116,7 P = 1424,08 T= 689,63 A = 0,083811
fraction RPS psia R ft2
Dari hasil kalibrasi spinner (lampiran 2), didapatkan nilai gradien sebesar 0,4 RPS/(m/min). Faktor koreksi diambil nilai satu, dengan mempertimbangkan nilai viskositas fluida campuran yang kecil (gas) dan kecepatan alir fluida yang tinggi akan menghasilkan Reynolds Number yang besar (validasi asumsi pada lampiran 5). Didapatkan nilai kecepatan fluida total sebesar 296,18 m/min. Data produksi digunakan sebagai indikator keakuratan hasil interpretasi PLT: Tabel 2. Data Produksi (Qmax)
Smolen (1996) tidak tepat untuk diterapkan pada kasus ini. Ketidakakuratan interpretasi ini disebabkan oleh salahnya interpretasi pola aliran. Data produksi (Tabel 2) dan data formation volume factor (Tabel 3) digunakan untuk mengetahui kecepatan slip dengan memodifikasi persamaan 7: (12) Dengan menggunakan persamaan 11 dan 7, didapatkan nilai kecepatan slip sebesar 321,93 m/min dan hasil interpretasi data PLT: Tabel 5. Interpretasi Data PLT (Qmax), Setelah Dikoreksi Vslip = fasa water gas
Q-w = 437,67 STB/day Q-g = 9,48 MMSCF/d γ-w = 1,005642 SG γ-g = 0,797 SG
Formation Volume Factor pada kondisi temperatur dan tekanan stasiun teratas, dihitung dengan menggunakan korelasi yang dikembangkan oleh McCain (1988) 5(lampiran 3), didapatkan: Tabel 3. Formation Volume Factor Bg = 0,012047 ft3/scf Bw = 1,047896 bbl/STB
Digunakan metode perhitungan kecepatan slip oleh Smolen (1996)10 pada persamaan 10. Kecepatan superfisial didapatkan oleh persamaan 5 dan 9. Menggunakan data pada Tabel 1, 2 dan 3, dapat dibandingkan besarnya debit aliran hasil interpretasi (Q-spinner) dengan data produksi yang terekam (Q-prod):
321,93 m/min u Q-spinner (m/min) 6,75 454,21 289,43 9,51
Q-prod
error
437,67 9,48
0,04 0,00
Didapatkan hasil interpretasi yang lebih akurat (tingkat error lebih kecil dari 10%). Hal ini menunjukkan bahwa terdapat kecepatan slip pada fraksi fasa berat lebih kecil dari 0,35. Data log densitas dan fraksi berat pada lampiran 8 menunjukkan respon yang cenderung stabil pada nilai kecil. Adanya kecepatan slip dan karakteristik respon log seperti ini menunjukkan pola aliran fluida adalah annular mist (untuk pola aliran slug/churn, respon log tidak akan stabil). Interpretasi Kontribusi Tiap Interval Perforasi Data respon PLT pada kondisi aliran Qmax: Tabel 6. Data Tiap Stasiun station 3 2 1
depth (m) 3515 3560 3860
RPSo
Yw
116,72 104,18 0,00
0,19 0,20 0,98
Sesuai dengan persamaan 9 dan 6, Tabel 4. Interpretasi Data PLT (Qmax), Sebelum Dikoreksi Vslip = fasa water gas
(11)
0 m/min u (m/min) 56,34 239,84
Q-spinner
Q-prod
error
3791,48 7,88
437,67 9,48
7,66 0,17
Sesuai dengan pernyataan Smolen (1996)10, ketika didapatkan nilai fraksi fasa berat (air) yang lebih kecil dari 0,35 maka aliran berpola mist flow. Oleh karenanya, kecepatan slip sama dengan nol. Ternyata, hal ini menghasilkan interpretasi debit aliran dengan tingkat error lebih besar dari 10% (tidak akurat). Maka dapat diambil kesimpulan sementara bahwa metode
dengan nilai kecepatan total fluida dan kecepatan aliran fasa air yang sama, untuk fraksi air yang semakin besar, maka semakin besar pula kecepatan slip yang dihasilkan. Pada pola aliran bubble dihasilkan kecepatan slip yang terbesar dan pada mist dihasilkan kecepatan slip terkecil (sesuai dengan persamaan 10). Namun, hubungan ini tidak berlaku untuk aliran annular mist, akibat nilai kecepatan alir fasa air yang jauh lebih kecil 7. Untuk interpretasi kontribusi tiap laju fasa untuk tiap interval perforasi, digunakan nilai kecepatan
slip yang sama (321,93 m/min) untuk fraksi holdup air lebih kecil dari 0,35. Untuk fraksi holdup air lebih besar dari 0,35, digunakan persamaan 10, sesuai dengan pernyataan Smollen (1996)10. Ketika spinner tidak berputar (bernilai 0 RPS), kecepatan fluida tidak dapat terdeteksi. Pada kondisi tersebut, diasumsikan kecepatan fluida adalah sebesar kecepatan threshold.
1.
2.
3.
4. Tabel 7. Interpretasi PLT Tiap Stasiun (Qmax) station 3 2 1
v-slip (m/min) 321,93 321,93 38,18
v-f (m/min) 296,18 264,42 0,72
Q-w (STB/D) 454,21 81,98 5,25
Q-g (mmscfd) 9,51 8,65 0,02
Tabel 8. Interpretasi Kontribusi Tiap Interval Perforasi (Qmax) perfo depth Q-w (cont) Q-g (cont) (m) (STB/D) (mmscfd) 3533-3534 372,23 0,86 3578-3587 76,73 8,63 3841-3843 5,25 0,02
Dengan menggunakan metode yang sama, telah dilakukan interpretasi data PLT pada sumur X dengan kondisi aliran Q1/2max (lampiran 4). Hasil yang akurat yang dihasilkan (tingkat error lebih kecil dari 10%) menunjukkan bahwa metode interpretasi pada tulisan ini dapat diaplikasikan pada kasus lain. Kesimpulan Dari studi kasus ini, dapat diambil kesimpulan: 1. Logging produksi dan data produksi (debit aliran fluida tiap fasa) dapat digunakan untuk mengidentifikasi pola aliran fluida di dalam sumur. 2. Korelasi pada buku Smollen (1996)10 tidak tepat untuk diaplikasikan pada pola aliran annular mist. 3. Kecepatan slip pada aliran annular mist dapat dihitung dengan menggunakan bantuan data produksi (debit aliran fluida tiap fasa). 4. Pola aliran pada stasiun paling atas sumur X adalah annular mist. 5. Dengan mengetahui pola aliran fluida dalam lubang sumur yang tepat, maka kontribusi aliran tiap fasa untuk tiap lapisan yang akurat dapat diraih. Hal ini dapat digunakan untuk menghasilkan keputusan workover yang tepat. Saran Beberapa hal yang dapat dikembangkan dengan mengetahui pola aliran fluida dalam sumur dari studi kasus ini adalah:
Menutup (plug) interval perforasi pengkontribusi air terbesar (35333534m). Menentukan tekanan kepala sumur maksimal agar tidak terjadi liquid loading (lampiran 6 mencantumkan metode yang dapat digunakan). Mencari korelasi untuk menentukan kecepatan slip pada pola aliran annular mist, dengan adanya tambahan data. Membuktikan bahwa data PLT dapat digunakan untuk menentukan pola aliran yang lain (selain mist dan annular mist), dengan adanya tambahan data.
Daftar Notasi = luas area penampang pipa (ft2) = gradient kalibrasi (RPS/(m/min)) = formation volume factor gas(cuft/scf) = formation volume factor air(bbl/STB) = diameter penampang sumur (in) = specific gravity gas = specific gravity air = viskositas fluida (cp) = Reynolds Number = tekanan (psia) = debit gas dalam sumur (cuft/d) = debit air dalam sumur (bbl/d) = debit produksi gas (MMSCF/d) = debit produksi air (STB/d) = densitas gas (gr/cc) = densitas campuran (gr/cc) = densitas air (gr/cc) = banyaknya putaran = banyaknya putaran ketika Vkabel statis (bernilai nol) = gaya tegangan antar permukaan (dynes/cm) θ = kemiringan sumur (derajat) = temperatur (˚R) = kecepatan superfisial gas (m/min) = kecepatan superfisial air (m/min) = kecepatan aliran fluida total (m/min) = kecepatan aliran gas (m/min) = kecepatan slip (m/min) = kecepatan threshold (m/min) = kecepatan aliran air (m/min) = kecepatan aliran gas kritikal (m/min) = fraksi gas = fraksi air = faktor kompresibilitas gas Referensi 1. Atlas Wireline Services: “Interpretive Methods for Production Well Logs”, Third Edition, 1988. 2. Brill, J. P., and Mukherjee, H.: “Multiphase Flow in Wells”, Henry L.
Doherty Memorial Fund of AIME, Society of Petroleum Engineers Inc., Richardson, Texas, 1999. 3. Ding, Z. X., Ullah, K., and Huang, Y.: “A Comparison of Predictive Oil/Water Holdup Models for Production Log Interpretation in Vertical and Deviated Wellbores”, paper presented at SPWLA 35th Annual Logging Symposium, 1994. 4. Hendy, and Hernansjah: “Perhitungan Laju Alir Tiga Fasa Pada Sumur Comingle Dengan Data Logging Produksi”, Tugas Akhir Mahasiswa Teknik Perminyakan ITB, 2010 5. McCain, William D. Jr.: “The Properties of Petroleum Fluids”, PennWell Books, Tulsa Oklahoma, 1990. 6. Nicolas, Y., and Witterholt, E. J.: “Measurements of Multiphase Fluid Flow”, SPE paper 4023, 1972. 7. Ouyang, L.B: “Field Case Histories Demonstrating Critical Role of PLT Flow Model Selection for Improved Water Shut-off Results in Offshore Thailand”, IPTC paper 12628, 2008 8. Roesner, R. E., LeBlanc, A. J., and Davarzani, M. J.: “Effects of Flow Regimes on Production Logging Instruments’ Response”, SPE paper 18206, 1988. 9. Schlumberger: “Production Log Interpretation”, 1973. 10. Smolen, J. J.: “Cased Hole and Production Log Evaluation”, PennWell Books, Tulsa, Oklahoma, 1996. 11. Turner, R.G., Hubbard, M.G., and Dukler, A.E.: “Analysis and Prediction of Minimum Flow Rate for the Continuous Removal of Liquid From Gas Wells”, paper SPE 2198, 1969.
Lampiran 3. 1.
Penurunan rumus kecepatan superfisial:
Formation Volume Factor5: Untuk dry gas: (11) Untuk air formasi: (12) di mana,
2.
Kalibrasi log spinner (Flow Response Line) disusun menggunakan data respon spinner pada saat kondisi sumur tidak mengalir (shut-in). Kalibrasi dilakukan dengan cara mengubah-ubah respon spinner hingga didapatkan nilai R2 sebesar 1 pada saat peralatan naik dan turun. Respon spinner yang mengakibatkan nilai R2 lebih kecil dari 1 tersebut adalah akibat adanya kesalahan pembacaan ataupun kesalahan alat.
4.
Data-data yang digunakan untuk analisa pada debit aliran Q1/2max sebesar 4,87 MMSCF/day adalah: Tabel 10a. Data PLT Stasiun Teratas (Q1/2max) Y-w = 0,13071 RPSo = 42,9 P = 3163,53 T= 709,55 Bg = 0,00558 Bw = 1,051907
Tabel 9. Data Respon Spinner Kedalaman 3515m, Pada Kondisi Shut-In Vcable (m/min) 30 20 10 0 -10 -20 -30
RPS RPS (sebelum) (setelah) 15,11 12,09 9,39 8,13 4,70 4,19 0,00 0,00 -4,65 -3,17 -8,36 -7,07 -12,95 -11,07
Tabel 10b. Data Produksi (Q1/2max) Q-w = Q-g =
296,30 STB/day 4,91 MMSCF/d
Tabel 11. Data Tiap Stasiun station 3 2 1
Gambar 9. Kalibrasi Spinner Pada Kedalaman 3515m, Sebelum Dikoreksi
fraction RPS psia R ft3/scf bbl/STB
depth (m) 3515 3560 3860
RPSo
Yw
42,92 36,18 0,00
0,13 0,12 1,00
Dari nilai fraksi fasa berat yang lebih kecil dari 0,35 dan respon log yang statis menunjukkan bahwa pola aliran pada stasiun teratas adalah annular mist. Kecepatan slip pada pola aliran ini didapatkan dengan menggunakan metode yang sudah dibahas sebelumnya (menggunakan data produksi), sebesar 45,79 m/min. VPCF pada kasus ini diasumsikan sebesar 0,663 (validasi asumsi pada lampiran 5). Tabel 12. Interpretasi Data PLT (Q1/2max) Vslip = fasa water gas
Gambar 10. Kalibrasi Spinner Pada Kedalaman 3515m, Setelah Dikoreksi
45,79 m/min u Q-spinner (m/min) 4,37 293,19 68,89 4,89
Q-prod
error
296,30 4,91
0,01 0,00
Sebagai bahan validasi untuk menunjukkan ketepatan asumsi, dilakukan metode yang sama untuk debit aliran yang lain. Data-data yang dibutuhkan didapat dari hasil regresi debit gas.
Tabel 13. Interpretasi PLT Tiap Stasiun (Q1/2max) station 3 2 1
v-slip (m/min) 45,79 45,79 38,18
v-f (m/min) 73,27 61,83 0,72
Q-w (STB/D) 293,19 181,92 45,34
Q-g (mmscfd) 4,89 4,20 0,00
Tabel 19. Validasi Matching Debit Gas Hasil Interpretasi Dengan Hasil Produksi (Sudah Cocok)
Tabel 14. Interpretasi Kontribusi Tiap Interval Perforasi (Q1/2max)
kondisi
Q-g (prod) VPCF Q-g (int) error Q-g (MMSCF/d) (mmscfd) Q3/4max 7,194975 0,709 7,21 0,00 Q3/5max 5,82462 0,666 5,89 0,01 Q4/5max 7,65176 0,74 7,65 0,00
perfo depth Q-w (cont) Q-g (cont) (STB/D) (mmscfd) 3533-3534 111,26 0,69 3578-3587 136,58 4,19 3841-3843 45,34 0,00
5.
Tabel 20. Validasi Matching Debit Air Hasil Interpretasi Dengan Hasil Produksi (Sudah Cocok)
Nilai VPCF dipengaruhi oleh Reynolds Numbers10 (gambar 4). Walaupun nilai kecepatan alir fluida dan diameter lubang sumur pada persamaan Reynold Number (persamaan 4) dapat diketahui, namun nilai densitas campuran dan viskositas campuran fluida berpola aliran annular mist dalam lubang sumur tidak bisa diketahui secara akurat oleh karena peralatan PLT tidak dapat mendeteksi karakteristik fluida pada lapisan liquid-film di sepanjang dinding sumur. Metode lain untuk mengetahui nilai VPCF adalah dengan mengubah-ubah nilai VPCF hingga didapatkan nilai interpretasi debit yang cocok dengan nilai debit data produksi (matching).
kondisi
Q-w (prod) VPCF (STB/d) Q3/4max 366,9875 0,709 Q3/5max 324,578 0,666 Q4/5max 381,124 0,74
Q-w (int) error Q-w (STB/D) 371,5981 0,01 335,8022 0,03 381,111 0,00
Tabel 21. Hubungan Antara Debit Gas Dengan VPCF kondisi Qmax Q4/5max Q3/4max Q3/5max Q1/2max
Q-g (prod) (MMSCF/d) 9,48 7,65 7,19 5,82 4,91
VPCF 1 0,74 0,709 0,666 0,663
Tabel 15. Matching Debit Gas Hasil Interpretasi Dengan Hasil Produksi (Belum Cocok) kondisi
Q-g (prod) VPCF (MMSCF/d) Qmax 9,48 Q1/2max 4,91
Q-g (int) error Q-g (mmscfd) 1 9,51 0,00 1 7,12 0,45
Tabel 16. Matching Debit Air Hasil Interpretasi Dengan Hasil Produksi (Belum Cocok) kondisi
Q-w (prod) VPCF (STB/d) Qmax 437,67 Q1/2max 296,30
Q-w (int) error Q-w (STB/D) 1 454,21 0,04 1 609,40 1,06
Gambar 11. Plot Validasi VPCF
Terlihat pada gambar 11, terdapat hubungan yang sebanding antara debit gas dengan VPCF. Hal ini membuktikan bahwa asumsi VPCF yang digunakan adalah valid.
Tabel 17. Matching Debit Gas Hasil Interpretasi Dengan Hasil Produksi (Sudah Cocok) kondisi
Q-g (prod) VPCF Q-g (int) error Q-g (MMSCF/d) (mmscfd) Qmax 9,48 1 9,51 0,00 Q1/2max 4,91 0,663 4,84 0,01
Tabel 18. Matching Debit Air Hasil Interpretasi Dengan Hasil Produksi (Sudah Cocok) kondisi
Q-w (prod) VPCF Q-w (int) error Q-w (STB/d) (STB/D) Qmax 437,67 1 454,21 0,04 Q1/2max 296,30 0,663 286,48 0,03
6.
Fenomena liquid-loading terjadi ketika kecepatan aliran gas lebih kecil dari kecepatan kritikalnya untuk dapat membawa butiran air secara sempurna11. Butiran air yang tidak terangkat akan mengendap dan terakumulasi di dasar sumur. Ketika tekanan hidrostatis akumulasi butiran air ini lebih besar daripada tekanan reservoir, maka sumur akan mati. (13)
Dengan mengetahui kecepatan aliran gas, tekanan dan temperatur dasar sumur pada dua debit alir yang berbeda, dapat dilakukan iterasi pada persamaan 13 untuk mencari nilai kecepatan kritikal gas. Selanjutnya, dengan memplot kecepatan kritikal gas tersebut ke dalam kurva tekanan kepala sumur terhadap kecepatan aliran gas, dapat didapatkan tekanan kepala sumur maksimum agar tidak terjadi liquid-loading.
7.
Data Produksi Sumur X
8.
Data PLT sumur X, kondisi aliran Qmax (bagian 1)
9.
Data PLT sumur X, kondisi aliran Qmax (bagian 2)
10. Data PLT sumur X, kondisi aliran Q1/2max
11. Data PLT sumur X, kondisi shut-in