Modul 4: Numerical Formulation — Continuous Model

Ringkasan

Modul ini menjelaskan dua hal dasar yang jadi fondasi numerical method dalam reservoir simulation: pertama, tiga gaya alami yang bekerja pada fluida di reservoir (viskos, gravitasi, kapiler); kedua, bagaimana persamaan diferensial yang biasanya sulit diselesaikan secara eksak dapat didekati secara numerik menggunakan Taylor series dan skema finite difference. Modul ini juga menunjukkan mengapa pendekatan numerik dibutuhkan: reservoir nyata umumnya tidak homogen, bentuknya tidak beraturan, dan sifat fluidanya berubah terhadap tekanan (nonlinier), sehingga analytical solution jarang dapat digunakan.

Signifikansi

Modul 3 menurunkan persamaan aliran fluida secara analitik untuk kasus-kasus yang disederhanakan: geometri sederhana, sifat fluida dan batuan yang dianggap seragam. Di dunia nyata, asumsi-asumsi itu jarang berlaku: reservoir memiliki bentuk tidak beraturan, sifat batuan yang bervariasi antar lokasi, dan sifat fluida yang berubah seiring tekanan turun. Modul 4 ini adalah jembatan menuju metode yang digunakan simulator reservoir sungguhan: mengubah persamaan diferensial kontinu menjadi persamaan aljabar yang dapat dihitung komputer secara berulang (iteratif), menggunakan pendekatan finite difference.

Memahami Taylor series dan turunannya (forward, backward, central difference) penting karena pilihan skema ini menentukan akurasi dan stabilitas hasil simulasi: topik yang akan terus muncul di modul-modul simulasi numerik berikutnya.

Konsep Kunci

• Gaya viskos: gaya akibat gesekan/hambatan aliran fluida melalui pori batuan, sebanding dengan viskositas fluida dan laju aliran, berbanding terbalik dengan permeabilitas.
• Gaya gravitasi: gaya akibat perbedaan densitas antar fluida (misalnya minyak vs air) yang mendorong pemisahan/segregasi fluida secara vertikal.
• Gaya kapiler: gaya akibat tegangan permukaan antar fluida, terkait dengan gradien saturasi fluida di suatu titik.
• Bilangan gravitasi (Gravity Number): rasio antara gaya gravitasi dan gaya viskos, digunakan untuk menilai fenomena mana yang lebih dominan dalam suatu aliran.
• Taylor Series: cara menuliskan nilai fungsi di suatu titik sebagai jumlah tak hingga dari turunan-turunannya di titik lain; dasar matematis semua skema finite difference.
• Truncation error: kesalahan yang muncul karena Taylor series yang sesungguhnya tak hingga dipotong hanya sampai beberapa suku saja.
• Skema finite difference (forward, backward, central difference): cara mendekati turunan suatu fungsi menggunakan selisih nilai fungsi di titik-titik diskret yang berdekatan.
• Orde akurasi, $O(\Delta x)$: notasi yang menunjukkan seberapa cepat galat mengecil ketika jarak antar titik grid ($\Delta x$) diperkecil.

Pembahasan

Tiga Gaya yang Bekerja di Reservoir

Sebelum masuk ke numerical method, modul ini mengingatkan bahwa aliran fluida di reservoir dipengaruhi tiga gaya:

$\text{Gaya Viskos} = \frac{\mu}{0.00633\,k}\,u \quad [\text{psi/ft}]$

$\text{Gaya Gravitasi} = \frac{\Delta\rho}{144} \quad [\text{psi/ft}]$

$\text{Gaya Kapiler} = \nabla P_c = \frac{dP_c}{dS}\nabla S \quad [\text{psi/ft}]$

Gaya viskos berasal dari “gesekan” fluida saat mengalir melalui pori-pori batuan: semakin kental fluidanya ($\mu$ besar) atau semakin kecil permeabilitasnya ($k$ kecil), semakin besar gaya ini. Gaya gravitasi muncul dari perbedaan densitas antar fluida ($\Delta\rho$): inilah yang menyebabkan minyak, air, dan gas cenderung memisah secara vertikal seiring waktu. Gaya kapiler muncul dari tegangan permukaan antar fluida yang tidak saling campur, dan besarnya tergantung seberapa tajam perubahan saturasi fluida ($\nabla S$) di suatu titik.

Modul menekankan tiga kondisi khusus: apabila tidak ada aliran, gaya viskos = 0; apabila tidak ada perbedaan densitas (fluida satu fasa), gaya gravitasi = 0; apabila tidak ada gradien saturasi (fluida satu fasa), gaya kapiler = 0. Rasio gaya gravitasi terhadap gaya viskos disebut Gravity Number, yang berguna untuk menilai kapan efek gravitasi (misalnya segregasi gas-minyak) menjadi signifikan dibanding dorongan tekanan.

Kebutuhan akan Pendekatan Numerik

Modul lalu beralih ke topik utama: pendekatan finite difference secara kontinu. Alasan pendekatan numerik dibutuhkan dijelaskan lewat tiga poin: pertama, banyak masalah aliran reservoir tidak memiliki analytical solution; kedua, reservoir nyata heterogen dan bentuknya tidak beraturan sehingga tidak cocok dengan asumsi geometri sederhana pada analytical solution; ketiga, persamaan alirannya sendiri bersifat nonlinier karena sifat fluida (viskositas, faktor compressibility) berubah seiring tekanan.

Taylor Series: Fondasi Semua Skema Finite Difference

Taylor Series menyatakan nilai fungsi $p$ di titik $(x+\Delta x)$ sebagai jumlah nilai $p$ di titik $x$ ditambah kontribusi dari turunan-turunannya:

$p(x+\Delta x) = p(x) + \Delta x\,p'(x) + \frac{\Delta x^2}{2!}p''(x) + \frac{\Delta x^3}{3!}p'''(x) + \cdots$

Secara fisis, Taylor series mengestimasi nilai suatu besaran, misalnya tekanan, pada lokasi yang berdekatan berdasarkan nilai fungsi dan turunannya di suatu titik acuan. Deret tersebut bersifat tak hingga dan eksak secara teoritis, tetapi implementasi numerik hanya mempertahankan sejumlah suku pertama. Suku yang diabaikan menghasilkan truncation error; penambahan jumlah suku meningkatkan akurasi sekaligus kompleksitas komputasi.

Ekspansi maju pada titik $i+1$ dan ekspansi mundur pada titik $i-1$ menghasilkan dua persamaan yang dapat dikombinasikan untuk mengaproksimasi turunan tekanan. Formulasi ini menjadi dasar finite-difference scheme.

Skema Finite Difference: Forward, Backward, Central

Dari Taylor series yang dipotong hanya sampai suku orde pertama, didapat tiga cara mendekati turunan pertama $\partial p/\partial x$:

• Forward difference ($\dfrac{\partial p}{\partial x} \approx \dfrac{p_{i+1}-p_i}{\Delta x}$) menggunakan nilai pada titik acuan dan titik berikutnya.
• Backward difference ($\dfrac{\partial p}{\partial x} \approx \dfrac{p_i - p_{i-1}}{\Delta x}$) menggunakan nilai pada titik acuan dan titik sebelumnya.
• Central difference untuk turunan kedua ($\dfrac{\partial^2 p}{\partial x^2} \approx \dfrac{p_{i-1}-2p_i+p_{i+1}}{(\Delta x)^2}$) menggunakan nilai pada kedua sisi titik acuan.

Forward dan backward difference memiliki galat berorde $O(\Delta x)$, artinya apabila jarak antar titik ($\Delta x$) diperkecil setengahnya, galatnya juga mengecil kira-kira setengahnya (turun linier). Central difference untuk turunan kedua memiliki galat berorde $O(\Delta x^2)$, jauh lebih akurat karena apabila $\Delta x$ diperkecil setengahnya, galatnya mengecil menjadi seperempatnya.

Skema Orde Tinggi (High Order Finite Difference)

Modul ini juga menunjukkan bahwa dengan melibatkan lebih banyak titik grid di sekitar titik yang dihitung (misalnya $p_{i-2}, p_{i-1}, p_i, p_{i+1}, p_{i+2}$), skema forward dan backward difference dapat ditingkatkan akurasinya dari $O(\Delta x)$ menjadi $O(\Delta x^2)$, dan skema central difference untuk turunan kedua dapat ditingkatkan dari $O(\Delta x^2)$ menjadi $O(\Delta x^4)$. Konsekuensinya adalah rumus yang lebih rumit (melibatkan lebih banyak titik dan koefisien) namun hasilnya jauh lebih presisi untuk jumlah grid yang sama. Ini adalah trade-off klasik dalam numerical method antara kompleksitas komputasi dan akurasi hasil.

Kesimpulan Utama

• Tiga gaya alami di reservoir: viskos (akibat gesekan aliran), gravitasi (akibat perbedaan densitas), dan kapiler (akibat gradien saturasi).
• Gravity Number adalah rasio gaya gravitasi terhadap gaya viskos.
• Pendekatan numerik dibutuhkan karena reservoir nyata heterogen, bentuknya tidak beraturan, dan persamaannya nonlinier.
• Taylor Series adalah dasar matematis dari semua skema finite difference; memotongnya menimbulkan truncation error.
• Forward dan backward difference berakurasi $O(\Delta x)$; central difference untuk turunan kedua berakurasi $O(\Delta x^2)$.
• Melibatkan lebih banyak titik grid meningkatkan orde akurasi, namun menambah kompleksitas perhitungan.
• Konsep-konsep ini adalah fondasi menuju penyusunan persamaan aljabar (matriks) yang digunakan simulator reservoir untuk menghitung tekanan dan saturasi tiap sel grid.