Modul 3: Reservoir Fluid-Flow Models

Ringkasan

Modul ini membahas bagaimana persamaan matematis untuk aliran fluida di dalam batuan reservoir diturunkan. Dasarnya adalah Darcy’s Law, yang dipadukan dengan hukum kekekalan massa dan persamaan keadaan (equation of state) fluida. Dari kombinasi ini, dibahas beberapa bentuk aliran: aliran linier (garis lurus) dan aliran radial (menuju satu titik, seperti ke arah sumur), masing-masing untuk incompressible fluid, sedikit termampatkan, dan gas (sangat termampatkan). Modul juga memperkenalkan konsep tekanan reservoir rata-rata dan aliran multifasa horizontal, saat minyak, air, dan gas mengalir bersamaan.

Signifikansi

Simulasi reservoir memerlukan formulasi persamaan fisik yang merepresentasikan mekanisme aliran secara tepat. Modul 3 menyajikan fondasi matematis bagi reservoir simulation; numerical model pada modul-modul berikutnya, termasuk numerical formulation pada Modul 4, merupakan pendekatan untuk menyelesaikan persamaan tersebut ketika kompleksitas kondisi lapangan tidak memungkinkan penyelesaian analitik.

Memahami perbedaan antara aliran linier dan radial, serta bagaimana sifat fluida (incompressible, slightly compressible, gas) mengubah bentuk persamaan, membantu menjelaskan mengapa simulator reservoir butuh iterasi dan linearisasi ketika fluidanya adalah gas atau minyak dengan sifat compressibility yang signifikan.

Konsep Kunci

• Darcy’s Law: hubungan linier antara laju aliran fluida dan gradien tekanan dalam media berpori, dengan permeabilitas sebagai konstanta proporsionalitas.
• Potensial aliran ($\Phi$): versi “tekanan yang disesuaikan” yang memperhitungkan efek gravitasi/ketinggian, digunakan untuk aliran yang tidak horizontal.
• Fluida incompressible, slightly compressible, dan compressible: tiga kategori fluida berdasarkan seberapa besar volumenya berubah terhadap tekanan; menentukan bentuk persamaan alirannya.
• Aliran linier vs radial: aliran linier terjadi searah garis lurus (misalnya di core sample atau antar sumur berjajar), sedangkan aliran radial mengalir menuju/dari satu titik (khas aliran menuju sumur).
• Fungsi pseudopressure $m(p)$: transformasi matematis yang digunakan untuk aliran gas supaya persamaan alirannya dapat “dilinierkan” walau viskositas dan faktor z gas berubah terhadap tekanan.
• Tekanan reservoir rata-rata: tekanan efektif reservoir, yang secara praktis diwakili oleh tekanan pada radius sekitar $0.61\,r_e$ dari sumur.

Pembahasan

Dasar Penurunan Persamaan Aliran

Materi dibuka dengan penegasan bahwa setiap persamaan aliran fluida dalam media berpori adalah gabungan tiga hal: hukum kekekalan massa, hukum aliran dalam media berpori (Darcy’s Law), dan persamaan keadaan (equation of state) yang menghubungkan densitas fluida dengan tekanan. Bentuk akhir persamaannya tergantung pada jenis fluida dan rezim alirannya.

Untuk menyederhanakan penurunan awal, digunakan empat asumsi: aliran laminar (viskos, bukan turbulen), steady-state (tekanan di tiap titik reservoir tidak berubah terhadap waktu), fluida incompressible, dan formasi homogen.

Aliran Linier Incompressible Fluid

Untuk aliran satu arah (horizontal) dari fluida yang volumenya dianggap tidak berubah oleh tekanan, Darcy’s law diintegralkan sepanjang jarak L:

$\frac{q}{A}\int_0^L dx = -\frac{k}{\mu}\int_{p_1}^{p_2} dp \implies q = \frac{kA(p_1-p_2)}{\mu L}$

Di sini $q$ adalah laju alir, $k$ permeabilitas absolut, $A$ luas penampang aliran, $\mu$ viskositas fluida, dan $L$ jarak tempuh. Tanda minus pada Darcy’s law asli hilang karena arah aliran didefinisikan dari tekanan tinggi ($p_1$) ke tekanan rendah ($p_2$). Dalam satuan lapangan (bbl/day, md, psi, cp, ft), muncul konstanta konversi 0.001127 di depan rumus.

Apabila alirannya tidak horizontal, tekanan biasa ($p$) diganti dengan potensial $\Phi$, yang memasukkan efek berat kolom fluida:

$\Phi_i = p_i + \left(\frac{\rho}{144}\right)\Delta z_i$

Fluida akan mengalir dari titik dengan potensial lebih tinggi ke titik dengan potensial lebih rendah, bukan sekadar dari tekanan tinggi ke rendah, karena gravitasi ikut “mendorong” atau “menahan” aliran tergantung arah ketinggiannya.

Aliran Fluida yang Sedikit Termampatkan

Fluida seperti minyak sesungguhnya sedikit berubah volumenya terhadap tekanan, dinyatakan lewat compressibility c:

$V = V_{ref}\left[1+c\,(p_{ref}-p)\right]$

Karena laju alir q sebanding dengan volume yang mengalir, q juga mengikuti bentuk serupa. Ketika ini disubstitusikan ke Darcy’s law dan diintegralkan, hasilnya adalah persamaan logaritmik untuk laju alir referensi:

$q_{ref} = \left[\frac{0.001127\,kA}{\mu c L}\right] \ln\left[\frac{1+c(p_{ref}-p_2)}{1+c(p_{ref}-p_1)}\right]$

Bentuk logaritma natural (ln) ini muncul karena compressibility membuat densitas (dan karenanya laju alir volumetrik) berubah sepanjang lintasan aliran, sehingga integrasinya tidak sesederhana kasus incompressible.

Aliran Gas: Dasar Penggunaan p²

Gas adalah fluida paling kompresibel, sehingga persamaan keadaannya (hukum gas nyata, $pV = znRT$) harus ikut dimasukkan. Setelah laju alir gas di kondisi reservoir dikonversi ke kondisi standar (scf/day) dan disubstitusikan ke Darcy’s law, hasil integrasinya memunculkan selisih kuadrat tekanan:

$Q_{sc} = \frac{0.003164\,T_{sc}\,A\,k\,(p_1^2 - p_2^2)}{p_{sc}\,T\,L\,\bar{z}\mu_g}$

Munculnya $p^2$ (bukan $p$ biasa) adalah ciri khas aliran gas, konsekuensi langsung dari sifat gas yang densitasnya berbanding lurus dengan tekanan. Ini penting diingat karena artinya hubungan laju alir dan tekanan pada aliran gas pada dasarnya tidak linier seperti pada minyak.

Aliran Radial: Menuju Sumur

Berbeda dari aliran linier, aliran menuju sumur berbentuk radial: fluida datang dari segala arah menuju satu titik pusat. Luas penampang aliran di sini bukan konstan, melainkan semakin kecil saat mendekati sumur ($A_r = 2\pi r h$), sehingga kecepatan aliran semakin tinggi di dekat sumur. Untuk fluida incompressible, hasil integrasinya adalah persamaan produktivitas sumur yang sangat dikenal di teknik reservoir:

$Q_o = \frac{0.00708\,k\,h\,(p_e - p_w)}{\mu_o B_o \ln(r_e/r_w)}$

dengan $r_e$ jari-jari pengurasan (drainage radius), $r_w$ jari-jari sumur, dan $B_o$ formation volume factor minyak (mengoreksi perbedaan volume minyak di reservoir vs di permukaan). Bentuk $\ln(r_e/r_w)$ ini konsisten dengan pola aliran radial (gradien tekanan paling tajam terjadi tepat di sekitar sumur, lalu melandai seiring menjauh) sehingga profil tekanannya secara visual berbentuk kurva curam dekat sumur dan datar dekat batas reservoir.

Konsep terkait yang muncul di sini adalah tekanan reservoir rata-rata, yaitu tekanan efektif yang mewakili seluruh area pengurasan sumur. Secara praktis, tekanan pada radius sekitar $0.61\,r_e$ kurang lebih sama dengan tekanan rata-rata volumetrik seluruh area itu, sehingga rumus produktivitas dapat ditulis ulang menggunakan $\ln(r_e/r_w) - 0.5$.

Untuk gas, aliran radial menggunakan transformasi pseudopressure:

$m(p) = \psi = \int_0^p \frac{2p}{\mu_g z}\,dp$

Transformasi ini “melinierkan” persamaan aliran gas walau viskositas dan faktor compressibility gas ($z$) berubah-ubah terhadap tekanan, sehingga bentuk akhirnya dapat ditulis mirip kasus incompressible namun menggunakan $\psi$, bukan $p$.

Aliran Multifasa

Ketika minyak, air, dan gas mengalir bersamaan (kondisi paling umum di lapangan produksi) masing-masing fasa memiliki permeabilitas efektifnya sendiri ($k_o$, $k_w$, $k_g$) yang nilainya tergantung saturasi fasa tersebut. Persamaan Darcy ditulis terpisah untuk tiap fasa, dengan bentuk yang mirip namun permeabilitasnya bukan lagi permeabilitas absolut, melainkan permeabilitas efektif fasa itu.

Kesimpulan Utama

• Persamaan aliran fluida = kekekalan massa + Darcy’s Law + persamaan keadaan fluida.
• Asumsi dasar Darcy: aliran laminar, steady-state, fluida incompressible, formasi homogen.
• Aliran non-horizontal menggunakan potensial $\Phi$ (tekanan + koreksi gravitasi), bukan tekanan biasa.
• Fluida slightly compressible menghasilkan persamaan berbentuk logaritmik ($\ln$) karena densitas berubah sepanjang lintasan aliran.
• Aliran gas menghasilkan persamaan berbentuk kuadrat tekanan ($p_1^2 - p_2^2$) karena densitas gas sebanding dengan tekanan.
• Aliran radial (menuju sumur) menghasilkan persamaan produktivitas berbentuk $\ln(r_e/r_w)$, karena luas penampang aliran mengecil mendekati sumur.
• Tekanan reservoir rata-rata secara praktis diwakili oleh tekanan di radius sekitar $0.61\,r_e$.
• Aliran multifasa menggunakan permeabilitas efektif tiap fasa ($k_o$, $k_w$, $k_g$), bukan permeabilitas absolut tunggal.