Fourth SPE Comparative Solution Project: Comparison of Steam Injection Simulators

Ringkasan

Paper ini adalah studi banding keempat dalam seri SPE Comparative Solution Project, kali ini berfokus pada simulator injeksi uap (steam injection): teknik yang digunakan untuk memproduksi minyak berat (heavy oil) dengan cara memanaskan reservoir. Enam perusahaan diminta menjalankan simulator masing-masing pada tiga soal yang sama persis (data reservoir, fluida, dan skenario operasi sudah ditentukan), lalu hasilnya dibandingkan berdampingan. Tujuannya sederhana: apabila enam simulator berbeda diberi soal identik, seberapa mirip jawaban yang mereka hasilkan? Hasilnya secara umum cukup sepakat untuk besaran-besaran yang penting secara praktis, meski ada juga perbedaan yang cukup mencolok pada beberapa aspek. Studi ini penting sebagai salah satu tonggak awal validasi simulator reservoir thermal (panas).

Signifikansi Studi

Simulator reservoir sulit divalidasi karena analytical solution (rumus matematis pasti) untuk proses yang rumit seperti injeksi uap nyaris tidak ada. Satu-satunya cara praktis untuk menguji apakah sebuah simulator “masuk akal” adalah membandingkannya dengan simulator lain yang independen, dikembangkan oleh tim berbeda, pada soal yang sama persis. Apabila semua sepakat, itu memberi rasa percaya diri (meski tidak menjamin semuanya benar); apabila tidak sepakat, itu jadi sinyal ada sesuatu yang perlu digali: entah bug, perbedaan interpretasi soal, atau perbedaan numerical method.

Studi ini merupakan kelanjutan dari tiga studi banding SPE sebelumnya (black-oil oleh Odeh, coning tiga fase oleh Chappelear & Nolen, dan gas cycling kondensat oleh Kenyon & Behie), dan yang pertama yang menyasar khusus proses thermal: proses yang jauh lebih kompleks karena melibatkan perpindahan panas, penguapan/kondensasi, dan perubahan viskositas minyak yang drastis terhadap suhu.

Latar Belakang Studi

• Jenis reservoir/masalah: reservoir minyak berat (heavy oil, gravitasi 14° API) yang diproduksi dengan injeksi uap, ditambah satu kasus dengan minyak yang sebagian dapat menguap (distillable oil, 22° API).
• Diorganisir oleh K. Aziz (Stanford University), A.B. Ramesh (Scientific Software-Intercomp), dan P.T. Woo (Chevron Oil Field Research), atas permintaan panitia program SPE Reservoir Simulation Symposium 1985 di Dallas.
• Enam perusahaan berpartisipasi: Arco, Chevron, Computer Modelling Group (CMG), Mobil, Societe Nationale Elf Aquitaine (Elf/SNEA), dan Scientific Software-Intercomp (SSI). Hanya tiga di antaranya (Chevron, CMG, SSI) yang mengirim hasil untuk kasus ketiga (compositional).
• Tiga masalah uji yang dibandingkan:
  1. Injeksi uap siklik (cyclic steam / “huff and puff”) pada grid radial 2D, tiga siklus injeksi-rendam-produksi.
  2. Steamflood (pendorongan uap terus-menerus) pada pola sumur inverted nine-spot, disimulasikan sebagai satu-perdelapan elemen simetri.
  3. Sama seperti soal kedua, tetapi minyaknya bersifat “distillable” (dapat menguap sebagian) sehingga butuh model compositional, bukan black-oil.

Konsep Kunci

• Cyclic steam injection (huff and puff): teknik memanaskan reservoir di sekitar satu sumur dengan menginjeksi uap, mendiamkannya (soak), lalu memproduksi minyak yang telah menjadi lebih encer dari sumur yang sama.
• Steamflood / steam displacement: uap diinjeksikan terus-menerus dari sumur injeksi untuk mendorong minyak menuju sumur produksi di sekitarnya, biasanya dalam pola sumur tertentu (misalnya nine-spot).
• Inverted nine-spot pattern: susunan sumur di mana satu sumur injeksi dikelilingi delapan sumur produksi; karena simetri, cukup satu-perdelapan atau satu-perempat bagian pola yang perlu disimulasikan.
• Grid orientation effect: kesalahan numerik pada simulasi grid persegi (five-point) yang membuat hasil aliran fluida sangat bergantung pada arah grid, bukan pada fisika sesungguhnya; salah satu alasan beberapa simulator menggunakan skema nine-point.
• Skema nine-point vs five-point: cara menghitung aliran antar-blok grid; nine-point memperhitungkan blok-blok diagonal tambahan untuk mengurangi efek arah grid tadi.
• Model Stone (Stone’s second model): metode standar untuk menghitung relative permeability minyak pada sistem tiga fase (minyak-air-gas), digunakan semua peserta studi ini dengan variasi kecil.
• Cold water equivalent (CWE): cara menyatakan laju injeksi uap dalam satuan setara volume air dingin, memudahkan perbandingan antar simulator.
• Oil-steam ratio (OSR): rasio minyak yang diproduksi dibanding uap yang diinjeksikan, ukuran efisiensi ekonomi steamflood.

Pembahasan

Deskripsi Kasus

Ketiga soal menggunakan data reservoir yang serupa: kedalaman sekitar 1.500 ft, empat lapisan dengan permeabilitas berbeda, porositas 0,3, dan suhu awal reservoir 125°F yang dipanaskan oleh uap hingga 450°F (kualitas uap 0,7). Soal pertama menggunakan grid radial (r-z) dengan 13 titik grid ke arah radial, mensimulasikan tiga siklus steam-huff-and-puff selama total sekitar tiga tahun. Soal kedua dan ketiga menggunakan grid 9x5x4 yang mewakili satu-perdelapan pola inverted nine-spot seluas 2,5 acre, disimulasikan hingga 10 tahun. Soal ketiga memiliki komposisi minyak tiga pseudo-komponen (dua yang dapat menguap dan satu yang tidak), lengkap dengan data K-value (kesetimbangan fase) sebagai fungsi tekanan dan suhu.

Apa yang Dibandingkan

Setiap peserta menjelaskan model thermal mereka secara singkat (metode diferensial, cara menangani sumur, metode penyelesaian persamaan linear, dsb.), lalu mengirim hasil numerik untuk dibandingkan pada grafik yang sama: produksi minyak dan air kumulatif, laju produksi, kehilangan panas (heat loss) kumulatif, tekanan dasar sumur (BHP), waktu hingga steam breakthrough di setiap sumur produksi, serta beban komputasi (waktu CPU, jumlah langkah waktu, jumlah iterasi).

Hasil dan Temuan Utama

Untuk soal pertama (cyclic steam), kehilangan panas kumulatif di akhir siklus pertama bervariasi dari 465 hingga 666 MMBtu: variasi yang cukup besar meski pola kurva produksi antar peserta relatif mirip. Untuk soal kedua (steamflood), waktu yang dibutuhkan setiap sumur produksi untuk mencapai laju uap 10 STB/D bervariasi cukup jauh antar peserta (misalnya sumur produksi jauh berkisar dari 1.398 sampai 2.245 hari): perbedaan besar ini sebagian disebabkan pemilihan skema grid dan cara menghitung indeks produktivitas sumur yang berbeda-beda. CMG secara explicit melaporkan bahwa skema five-point pada soal ini menghasilkan efek grid-orientation yang parah sehingga uap “salah” terobos ke sumur jauh terlebih dahulu: mereka menambahkan sedikit tekanan kapiler untuk meredam gejala numerik ini. Untuk soal ketiga (compositional), hanya tiga peserta yang berpartisipasi dan hasilnya relatif konsisten.

Paper ini secara explicit mendaftar sembilan kemungkinan penyebab perbedaan hasil antar simulator: cara menangani suku antar-blok, perlakuan sumur, toleransi konvergensi iterasi, prosedur pemilihan langkah waktu, prosedur perhitungan kehilangan panas, penggunaan skema nine-point, pemilihan parameter kontrol program, kemungkinan bug, dan kesalahan input data. Karena persamaan aliran dasar dan sifat fluida sudah disamakan di awal, penulis menyimpulkan bahwa perbedaan terbesar kemungkinan berasal dari perlakuan sumur dan ketidakstabilan numerik akibat konvergensi yang kurang baik.

Temuan dan Implikasi

• Bahkan ketika data input (reservoir, fluida, kondisi operasi) disamakan persis, simulator berbeda tetap dapat menghasilkan angka yang cukup berbeda: terutama untuk besaran yang sensitif terhadap detail numerik seperti waktu breakthrough.
• Efek grid orientation adalah masalah numerik nyata pada steamflood dengan pola sumur, bukan sekadar isu teoretis: salah satu peserta studi ini melaporkannya secara langsung memengaruhi urutan breakthrough antar sumur.
• Cara menangani sumur (well treatment) dan konvergensi numerik disebut sebagai penyebab perbedaan hasil yang lebih besar dibanding perbedaan formulasi fisika dasar.
• Skema nine-point finite-difference dikembangkan justru untuk mengatasi bias arah grid pada simulasi pola sumur seperti nine-spot.
• Kesepakatan hasil antar simulator tidak membuktikan hasil itu benar: namun ketidaksepakatan adalah alasan kuat untuk menyelidiki lebih lanjut, sebagaimana ditegaskan penulis mengutip semangat studi banding SPE sebelumnya.
• Perbandingan seperti ini berbeda dari studi akademis yang membangun ulang model dari deskripsi publikasi; di sini simulator komersial asli langsung digunakan, sehingga perbedaan cenderung mencerminkan perbedaan implementasi nyata, bukan salah tafsir prosedur.
• Beban komputasi (waktu CPU, jumlah iterasi) juga sangat bervariasi antar simulator untuk soal yang identik, menunjukkan trade-off antara efisiensi numerik dan detail pemodelan.