Comparison of Solutions to a Three-Dimensional Black-Oil Reservoir Simulation Problem

Ringkasan

Paper ini adalah studi banding pertama dari seri SPE Comparative Solution Project. Penulisnya, Aziz S. Odeh dari Mobil Research and Development Corp., memberikan satu soal reservoir tiga dimensi yang persis sama kepada tujuh perusahaan berbeda, lalu meminta masing-masing menjalankannya dengan simulator black-oil milik sendiri. Hasilnya dikumpulkan dan dibandingkan dalam bentuk grafik untuk melihat seberapa dekat (atau seberapa jauh) simulator-simulator yang berbeda pendekatan numeriknya dapat sepakat terhadap kasus yang identik. Studi ini menjadi rujukan klasik karena menunjukkan secara terbuka bahwa perbedaan numerical method antar simulator dapat menghasilkan perbedaan hasil yang nyata, bukan sekadar isu teoretis.

Signifikansi Studi

Sebelum era software reservoir simulation yang matang dan tervalidasi seperti sekarang, insinyur reservoir di berbagai perusahaan mengembangkan simulatornya sendiri-sendiri, masing-masing dengan asumsi numerik dan skema solusi yang berbeda. Tanpa ada kasus uji bersama, sulit menjawab pertanyaan mendasar: apakah simulator-simulator ini benar-benar menghasilkan jawaban yang sama untuk masalah fisik yang sama? Studi banding seperti ini penting sebagai bentuk validasi silang (cross-validation): bukan dibandingkan dengan analytical solution yang “benar”, melainkan dibandingkan satu sama lain untuk melihat konsistensi.

Kasus uji yang dipilih di sini sengaja dibuat sederhana secara geometris (grid 10x10x3, hanya dua sumur) tetapi mengandung fenomena fisik yang menantang secara numerik: injeksi gas ke reservoir minyak undersaturated yang kemudian melepaskan gas bebas seiring tekanan turun di bawah titik gelembung (bubble point). Fenomena “gas keluar dari larutan” ini terkenal sulit ditangani secara numerik karena melibatkan pergantian jumlah fase secara dinamis di setiap sel grid.

Latar Belakang Studi

• Jenis reservoir/masalah: reservoir minyak 3D, undersaturated pada awalnya, dengan satu sumur injeksi gas dan satu sumur produksi minyak (masalah dua fase gas-minyak, saturasi air dianggap tetap/tidak bergerak).
• Diorganisir oleh Aziz S. Odeh (Mobil Research and Development Corp.), dipublikasikan sebagai SPE 9723 di Journal of Petroleum Technology, Januari 1981.
• Tujuh perusahaan/tim berpartisipasi: Amoco Production Co., Computer Modelling Group of Calgary (CMG), Exxon Production Research Co., Intercomp Resource Development and Engineering Inc., Mobil Research and Development Corp., Shell Development Co., dan Scientific Software Corp. (SSC).
• Dua skenario utama diuji (Case 1 (saturation pressure/bubble point dianggap konstan) dan Case 2 (saturation pressure bervariasi mengikuti saturasi gas) kasus bubble point variabel yang lebih sulit secara numerik).
• Hasil yang dibandingkan: laju produksi minyak vs waktu, GOR (gas-oil ratio) vs waktu, tekanan di sel injektor dan produser, serta saturasi gas di beberapa titik grid, sampai 8-10 tahun simulasi.

Konsep Kunci

• Black-oil model: model reservoir simulation yang menyederhanakan fluida menjadi tiga komponen (minyak, gas, air) dengan gas yang dapat larut dalam minyak, cocok untuk sebagian besar kasus produksi minyak konvensional.
• Bubble point (saturation pressure): tekanan di mana gas mulai keluar dari larutan minyak membentuk fase gas bebas; di bawah tekanan ini reservoir menjadi dua fase (minyak + gas).
• IMPES (Implicit Pressure, Explicit Saturation): numerical scheme yang menghitung tekanan secara implicit namun saturasi secara explicit; lebih cepat dihitung namun butuh langkah waktu (time step) lebih kecil agar stabil.
• Metode fully implicit: semua variabel (tekanan dan saturasi) dihitung implicit secara bersamaan; lebih stabil untuk langkah waktu besar namun perhitungan per langkah lebih berat.
• Adaptive Implicit Method (AIM): pendekatan yang dikembangkan SSC yang memilih tingkat keimplisitan berbeda untuk tiap sel grid secara dinamis, hanya menerapkan kalkulasi berat pada sel yang benar-benar membutuhkannya.
• Single-point vs two-point upstream weighting: cara menentukan nilai mobility fluida di batas antar sel grid; sebagian besar peserta menggunakan satu titik upstream, Shell menggunakan dua titik upstream.
• GOR (Gas-Oil Ratio): rasio volume gas terhadap volume minyak yang diproduksikan, indikator penting kapan gas bebas mulai ikut terproduksi di sumur.

Pembahasan

Deskripsi Kasus

Reservoirnya berupa grid areal 10x10 dengan 3 lapisan (layer) vertikal, masing-masing memiliki ketebalan, porositas, dan permeabilitas berbeda. Sumur injeksi gas ditempatkan di titik grid (1,1) dan sumur produksi minyak di titik grid (10,10): posisi diagonal berlawanan sehingga aliran fluida harus melewati seluruh reservoir. Reservoir mula-mula undersaturated (tekanan awal di atas tekanan bubble point) pada 4.800 psia. Data PVT lengkap (formation volume factor, viskositas, densitas, solution GOR) untuk minyak, gas, dan air disediakan dalam tabel, begitu juga data relative permeability minyak-gas. Batasan operasi juga diberikan: laju injeksi gas tetap 100 MMscf/D, laju produksi minyak maksimum 20.000 STB/D, tekanan alir dasar sumur minimum 1.000 psi, dan run dihentikan pada akhir 10 tahun atau ketika GOR mencapai 20.000 scf/STB.

Dua skenario dijalankan (Case 1 menganggap tekanan bubble point konstan sepanjang simulasi (lebih sederhana), sedangkan Case 2 membiarkan tekanan bubble point berubah mengikuti saturasi gas lokal) pendekatan yang jauh lebih realistis namun lebih sulit ditangani secara numerik karena simulator harus melacak titik gelembung yang bergerak di setiap sel.

Perbedaan Pendekatan Antar Simulator

Setiap peserta menggunakan strategi numerik yang berbeda-beda: Amoco menggunakan IMPES dengan penyesuaian semi-implicit pada laju sumur; CMG menjalankan model fully implicit tiga fase tiga dimensi dengan formulasi pseudo solution-gas yang unik untuk menangani hilangnya fase gas; Exxon menggunakan metode sequential implicit (GPSIM) yang memisahkan solusi tekanan dan saturasi secara berurutan; Intercomp (BETA II) menawarkan beberapa formulasi mulai dari IMPES sampai implicit penuh; Mobil (ALPURS) menggunakan metode fully implicit yang sepenuhnya coupled dengan Newton–Raphson iteration; SSC memperkenalkan Adaptive Implicit Method (AIM) yang baru dikembangkan saat itu; dan Shell menggunakan mode IMPES atau implicit dengan pembobotan upstream dua titik (berbeda dari peserta lain yang menggunakan satu titik).

Hasil dan Titik Perbedaan

Untuk saturation pressure (bubble point), ketujuh perusahaan sepakat dalam rentang yang sangat sempit (dalam 20 psi satu sama lain): jadi pada aspek ini semua metode konsisten. Namun untuk kurva laju produksi minyak, GOR, tekanan sel, dan saturasi gas terhadap waktu, hasilnya menunjukkan pola “kesepakatan sebagian”: pasangan-pasangan simulator tertentu cocok pada rentang waktu tertentu lalu menyimpang di rentang lain (misalnya Mobil cocok dengan CMG sampai tahun ke-7 lalu menyimpang, sementara Intercomp mengambil alih kecocokan itu di tahun-tahun berikutnya). Pola ini lebih terlihat jelas pada Case 2 (bubble point variabel) dibanding Case 1, konsisten dengan dugaan bahwa kasus yang secara numerik lebih menantang menghasilkan sebaran hasil yang lebih lebar antar simulator. Paper tidak secara explicit membuktikan simulator mana yang “paling benar” karena tidak ada analytical solution pembanding: perbandingan murni antar simulator.

Perlu diperhatikan bahwa artikel ini juga memuat tanggapan (discussion) dari William Hurst yang mempertanyakan nilai dari perbandingan semacam ini karena tidak dibandingkan terhadap analytical solution (misalnya solusi Muskat), dan menyarankan agar studi berikutnya memasukkan pembanding analitik tersebut.

Temuan dan Implikasi

• Simulator reservoir dengan numerical scheme berbeda (IMPES, sequential implicit, fully implicit, adaptive implicit) dapat menghasilkan jawaban yang berbeda untuk kasus fisik yang identik: perbedaan numerical method bukan hal sepele.
• Kasus dengan bubble point variabel (gas keluar dari larutan secara dinamis) secara konsisten lebih sulit dan menghasilkan sebaran hasil yang lebih besar dibanding kasus bubble point konstan.
• Variabel “mudah” seperti saturation pressure dapat sangat konsisten antar simulator, sementara variabel turunan seperti laju produksi dan GOR jauh lebih sensitif terhadap pilihan numerical scheme.
• Tidak ada satu numerical method yang secara jelas “menang” dalam studi ini: pilihan IMPES vs implicit vs adaptif lebih soal trade-off stabilitas dan biaya komputasi daripada akurasi mutlak.
• Studi banding semacam ini menjadi cikal bakal budaya benchmarking di industri reservoir simulation, yang kemudian dilanjutkan dengan studi kedua dan ketiga dengan kasus yang lebih kompleks.
• Kritik dari pembaca (Hurst) mengingatkan bahwa perbandingan antar simulator saja tidak cukup untuk memvalidasi “kebenaran”: perbandingan terhadap analytical solution tetap penting jika tersedia.
• Deskripsi kasus uji yang detail dan terbuka (grid, PVT, batasan operasi) menjadikan paper ini dapat direproduksi oleh siapa pun yang ingin menguji simulatornya sendiri terhadap kasus klasik ini.