IMPES Method

Ringkasan

IMPES (Implicit Pressure, Explicit Saturation) adalah salah satu metode solusi numerik paling klasik dan banyak digunakan dalam reservoir simulation. Idenya: dari persamaan mass balance minyak, air, dan gas yang sudah diturunkan sebelumnya (lihat Modul 5), ketiga persamaan itu digabung menjadi satu persamaan tekanan tunggal. Persamaan tekanan ini diselesaikan secara implicit (lewat sistem matriks) untuk mendapatkan tekanan di semua grid block pada waktu baru. Setelah tekanan diketahui, saturasi minyak, air, dan gas dihitung secara explicit: artinya langsung dihitung dari tekanan yang baru saja diperoleh, tanpa perlu menyelesaikan sistem matriks lagi. Nama “IMPES” persis menggambarkan urutan kerja ini: tekanan implicit, saturasi explicit.

Signifikansi

IMPES relevan karena menjadi salah satu metode solusi numerik yang paling banyak digunakan dalam sejarah reservoir simulation, terutama karena lebih efisien secara komputasi dibanding metode fully implicit (yang menyelesaikan tekanan dan saturasi sekaligus lewat Newton–Raphson iteration). Dengan memisahkan perhitungan tekanan (implicit) dari saturasi (explicit), ukuran sistem matriks yang harus diselesaikan jauh lebih kecil: hanya satu unknown (tekanan) per grid block, bukan empat unknown (tekanan + tiga saturasi) sekaligus.

Modul ini juga adalah kelanjutan langsung dari Modul 5. Persamaan mass balance minyak, air, dan gas yang diturunkan di modul sebelumnya adalah dasar formulasi yang di sini diolah menjadi satu persamaan tekanan yang dapat diselesaikan. Memahami IMPES berarti memahami bagaimana teori mass balance benar-benar diubah menjadi algoritma yang dapat dijalankan simulator, langkah demi langkah.

Konsep Kunci

• IMPES (Implicit Pressure, Explicit Saturation): metode yang menghitung tekanan secara implicit lewat sistem matriks, lalu saturasi secara explicit dari tekanan yang sudah diperoleh.
• 4 unknown per blok: tekanan ($p$), saturasi minyak ($S_o$), saturasi air ($S_w$), dan saturasi gas ($S_g$); IMPES butuh cara untuk memecah keempatnya menjadi tahap-tahap yang lebih sederhana.
• Syarat saturasi ($S_o + S_w + S_g = 1$): persamaan tambahan yang memastikan volume fluida di setiap blok selalu pas mengisi pore volume yang tersedia.
• Compressibility batuan ($c_f$): mengukur seberapa besar pore volume berubah akibat perubahan tekanan; digunakan untuk mengaproksimasi $V_p$ pada waktu baru.
• Compressibility fluida ($c_m$): mengukur seberapa besar formation volume factor ($B$) tiap fase berubah akibat perubahan tekanan.
• Compressibility total ($c_t$): gabungan compressibility batuan dan fluida (dibobot saturasi tiap fase), muncul sebagai koefisien utama dalam persamaan tekanan akhir.
• Persamaan tekanan (pressure equation): persamaan tunggal hasil penggabungan tiga persamaan fase (minyak, air, gas) yang diselesaikan lewat sistem matriks untuk mendapatkan tekanan di semua blok sekaligus.

Pembahasan

Formulasi Awal: Tiga Persamaan Fase

Seperti dibahas di Modul 5, mass balance untuk minyak, air, dan gas masing-masing membentuk persamaan sendiri-sendiri yang melibatkan koefisien antar-blok ($\alpha$) dan tekanan pada waktu baru ($p^{n+1}$), sama dengan suku perubahan volume fluida tersimpan ($V_p \cdot S/B$) antara waktu baru dan waktu lama, dibagi $\Delta t$. Masalahnya, tiap grid block memiliki empat unknown yang belum diketahui: tekanan $p$, dan tiga saturasi $S_o$, $S_w$, $S_g$. Tiga persamaan mass balance saja tidak cukup untuk menyelesaikan empat unknown: dibutuhkan satu persamaan tambahan, yaitu syarat bahwa ketiga saturasi harus berjumlah 1 (memastikan seluruh pore volume terisi penuh oleh ketiga fase).

Tiga Langkah Inti Metode IMPES

Metode IMPES memecah masalah ini menjadi tiga langkah berurutan:

1. Hitung koefisien persamaan tekanan: menyusun ulang ketiga persamaan fase menjadi satu persamaan yang hanya memuat tekanan sebagai unknown.
2. Selesaikan persamaan tekanan secara implicit: persamaan tekanan disusun sebagai sistem matriks (karena tekanan di satu blok terkait dengan tekanan di blok-blok tetangganya) dan diselesaikan sekaligus untuk seluruh grid, menghasilkan $p^{n+1}$ di semua blok.
3. Selesaikan persamaan saturasi secara explicit: dengan tekanan $p^{n+1}$ yang sudah diketahui, saturasi $S_o^{n+1}$, $S_w^{n+1}$, dan $S_g^{n+1}$ dihitung langsung dari persamaan fase masing-masing tanpa perlu menyelesaikan matriks lagi.

Dari Tiga Persamaan Menjadi Satu: Menyusun Persamaan Tekanan

Triknya adalah menjumlahkan ketiga persamaan mass balance (minyak, air, gas) sekaligus. Karena ketiganya dikalikan dengan formation volume factor $B$ masing-masing, penjumlahan ini menghasilkan satu persamaan gabungan yang di ruas kanannya memuat suku $1/V_p^{n+1}$ dikurangi jumlah semua ($B \cdot V_p S/B$) pada waktu lama: dan yang penting, syarat $S_o + S_w + S_g = 1$ sudah otomatis terpakai dalam penjumlahan ini sehingga saturasi tidak lagi menjadi unknown terpisah pada tahap ini.

Masalah berikutnya: ruas kanan persamaan gabungan masih mengandung $V_p^{n+1}$ dan $B^{n+1}$ (pore volume dan formation volume factor pada waktu baru), yang keduanya sesungguhnya fungsi dari tekanan baru yang justru sedang dicari. Di sinilah compressibility masuk sebagai pendekatan (aproksimasi).

Menggunakan Compressibility untuk Mengeliminasi $V_p$ dan $B$ pada Waktu Baru

$V_p^{n+1}$ didekati menggunakan compressibility batuan $c_f$, lewat hubungan chord slope (kemiringan rata-rata) antara $V_p$ dan tekanan:

$V_p^{n+1} \approx V_p^{n} \left(1 + c_f \left(p^{n+1} - p^{n}\right)\right)$

Artinya, perubahan pore volume akibat kenaikan/penurunan tekanan diaproksimasi sebagai fungsi linear sederhana dari selisih tekanan.

Dengan cara serupa, formation volume factor tiap fase $B_m^{n+1}$ didekati menggunakan compressibility fluida $c_m$:

$B_m^{n+1} \approx B_m^{n} \left(1 - c_m \left(p^{n+1} - p^{n}\right)\right)$

Dengan kedua aproksimasi ini, semua suku yang tadinya bergantung pada $V_p^{n+1}$ dan $B^{n+1}$ sekarang dapat ditulis ulang murni sebagai fungsi dari $p^{n+1}$ dan $p^{n}$: sehingga persamaan gabungan benar-benar menjadi persamaan tekanan tunggal yang siap diselesaikan.

Bentuk Akhir Persamaan Tekanan

Persamaan tekanan akhirnya berbentuk: jumlah koefisien antar-blok ($\alpha$) dikalikan tekanan di setiap blok, sama dengan ($V_p^{n} c_t / \Delta t$) dikalikan selisih tekanan ($p^{n+1} - p^{n}$):

$\sum \alpha \, p = \frac{V_p^{n} c_t}{\Delta t}\left(p^{n+1} - p^{n}\right)$

Di sini $c_t$ adalah compressibility total sistem, yaitu penjumlahan compressibility batuan ($c_f$) dengan compressibility tiap fase fluida ($c_m$) yang dibobot oleh saturasi fase tersebut ($S_m$). Persamaan inilah (satu per grid block, saling terkait lewat koefisien $\alpha$ antar blok tetangga) yang disusun menjadi sistem matriks dan diselesaikan implicit untuk mendapatkan medan tekanan baru di seluruh reservoir.

Kesimpulan Utama

• IMPES = Implicit Pressure, Explicit Saturation: tekanan dihitung implicit lewat matriks, saturasi dihitung explicit setelahnya.
• Setiap grid block memiliki 4 unknown ($p, S_o, S_w, S_g$); syarat $S_o + S_w + S_g = 1$ dibutuhkan untuk menutup sistem persamaan.
• Tiga persamaan mass balance (minyak, air, gas) dijumlahkan menjadi satu persamaan tekanan tunggal.
• Compressibility batuan ($c_f$) digunakan untuk mengaproksimasi pore volume pada waktu baru.
• Compressibility fluida ($c_m$) digunakan untuk mengaproksimasi formation volume factor pada waktu baru.
• Compressibility total ($c_t$) adalah gabungan $c_f$ dan $c_m$ (dibobot saturasi), menjadi koefisien utama persamaan tekanan.
• Setelah tekanan baru diperoleh, saturasi dihitung langsung (explicit) tanpa perlu menyelesaikan matriks lagi.
• IMPES populer karena lebih efisien secara komputasi dibanding metode fully implicit yang menyelesaikan semua unknown sekaligus.