Tenth SPE Comparative Solution Project: A Comparison of Upscaling Techniques

Ringkasan

Paper ini adalah studi perbandingan kesepuluh dalam seri SPE Comparative Solution Project, ditulis oleh M.A. Christie (Heriot-Watt University) dan M.J. Blunt (Imperial College), dipublikasikan pada 2001. Topiknya adalah upscaling: proses menyederhanakan model geologi yang sangat detail (jutaan sel kecil) menjadi grid simulasi yang jauh lebih kasar (ribuan sel) agar dapat dihitung dalam waktu wajar, tanpa kehilangan terlalu banyak akurasi. Sembilan tim/perusahaan menguji berbagai metode upscaling pada dua kasus: satu model kecil sederhana yang mudah dihitung penuh, dan satu model besar (1,1 juta sel) yang nyaris mustahil dihitung penuh tanpa penyederhanaan.

Signifikansi Studi

Model geologi reservoir modern sering dibangun dari data geostatistik dengan jutaan sel kecil untuk menangkap detail heterogenitas batuan. Namun menjalankan simulasi aliran fluida langsung pada grid sehalus itu dapat sangat mahal secara komputasi, bahkan hampir tidak mungkin untuk model penuh lapangan. Karena itu, industri sangat bergantung pada teknik upscaling untuk memampatkan model tersebut menjadi grid simulasi yang lebih kasar. Studi ini penting karena menjadi salah satu perbandingan independen pertama yang menguji berbagai pendekatan upscaling (mulai dari yang sangat sederhana (rata-rata permeabilitas) sampai yang memanfaatkan hasil simulasi grid halus untuk menyesuaikan kurva relative permeability) pada kasus uji yang sama, sehingga terlihat jelas metode mana yang lebih dapat diandalkan dan dalam kondisi seperti apa.

Latar Belakang Studi

• Dua model diuji: Model 1 adalah penampang vertikal 2D sederhana (100x1x20 sel, 2.000 sel total) untuk kasus injeksi gas, cukup kecil sehingga solusi grid halusnya mudah dihitung sebagai pembanding. Model 2 adalah waterflood 3D dari model geostatistik besar (60x220x85 = sekitar 1,122 juta sel), mewakili formasi Tarbert (35 lapisan atas) dan Upper Ness (50 lapisan bawah) dari deret Brent, yang aslinya dibuat untuk proyek PUNQ.
• Diorganisir oleh M.A. Christie dan M.J. Blunt; sembilan organisasi berpartisipasi: Chevron (CHEARS), Coats Engineering (SENSOR), GeoQuest (ECLIPSE 100 dan simulator streamline FRONTSIM), Landmark (VIP), Phillips Petroleum (SENSOR), Roxar (Nextwell/RMSsimgrid), Streamsim (simulator streamline 3DSL), TotalFinaElf (ECLIPSE dan 3DSL), dan University of New South Wales (CMG IMEX, hanya Model 1).
• Skenario utama yang diuji: upscaling ke grid sangat kasar tetap yang ditentukan (misalnya 5x1x5 untuk Model 1), upscaling ke grid pilihan bebas peserta (maksimum 100 sel untuk Model 1), serta perbandingan berbagai metode upscaling permeabilitas (arithmetic-harmonic, harmonic-arithmetic, power-law, flow-based dengan kondisi batas no-flow atau linear-pressure) pada grid kasar yang sama untuk Model 2.

Konsep Kunci

• Upscaling: proses menggabungkan banyak sel geologi kecil menjadi satu sel simulasi yang lebih besar, dengan menghitung properti “efektif” (permeabilitas, porositas) yang mewakili gabungan sel-sel aslinya.
• Pseudoization / pseudo relative permeability: teknik menyesuaikan kurva relative permeability pada grid kasar (bukan hanya permeabilitas absolut) agar perilaku aliran multifasenya meniru hasil grid halus, biasanya dengan bantuan hasil simulasi fine-grid sebagai acuan regresi.
• Upscaling berbasis aliran (flow-based upscaling): metode menghitung permeabilitas efektif dengan mensimulasikan aliran (biasanya aliran satu fasa sederhana) melalui blok kasar dan menyamakan fluksnya dengan hasil pada grid halus.
• Kondisi batas (boundary condition) pada upscaling: asumsi tentang bagaimana aliran diperlakukan di tepi blok saat menghitung permeabilitas efektif, misalnya “no-flow” (tidak ada aliran tegak lurus tepi) atau “linear pressure” (gradien tekanan linear); pilihan ini dapat sangat memengaruhi hasil upscaling.
• Simulasi streamline: pendekatan simulasi alternatif (digunakan FRONTSIM dan 3DSL di studi ini) yang melacak jalur aliran fluida, sering digunakan untuk memverifikasi hasil upscaling secara cepat pada model besar.
• Productivity Index (PI) sumur: angka yang menghubungkan laju alir sumur dengan penurunan tekanan di sekitarnya; sangat sensitif terhadap ukuran grid dan lokasi sumur relatif terhadap sel grid, dan berpengaruh besar terhadap prediksi tekanan rata-rata lapangan.
• Grid halus (fine grid) vs grid kasar (coarse grid): istilah untuk model geologi asli beresolusi tinggi versus model simulasi yang sudah disederhanakan; perbandingan keduanya adalah inti dari validasi metode upscaling.

Pembahasan

Deskripsi Kasus

Model 1 adalah kasus injeksi gas dua fasa (oil-gas) yang sangat sederhana secara geometri (hanya penampang vertikal tanpa kemiringan atau sesar) tetapi memiliki medan permeabilitas terkorelasi secara geostatistik. Karena ukurannya kecil (2.000 sel), solusi grid halusnya mudah dihitung semua peserta sebagai acuan, sehingga fokus pengujian ada pada seberapa dekat hasil upscaling ke grid 5x1x5 (dan grid pilihan bebas hingga 100 sel) terhadap solusi acuan tersebut.

Model 2 jauh lebih menantang (model waterflood 3D dengan lebih dari satu juta sel, meniru bagian formasi Tarbert (representasi lingkungan pengendapan near-shore yang mengalami progradasi) di atas Upper Ness (lingkungan fluvial)) kombinasi dua tipe formasi geologi berbeda yang membuat pola heterogenitas permeabilitasnya cukup kompleks, termasuk rasio permeabilitas vertikal-horizontal yang berbeda antara zona channel dan latar belakangnya. Dengan ukuran sebesar itu, menjalankan simulasi penuh pada grid asli nyaris tidak praktis: sehingga model ini sengaja dirancang untuk memaksa peserta benar-benar menggunakan teknik upscaling, bukan sekadar menjalankan grid halus secara langsung (meski beberapa peserta dengan simulator paralel dan streamline tetap berhasil menjalankannya untuk keperluan validasi).

Metode Upscaling yang Dibandingkan

Ada dua kelompok besar pendekatan yang muncul dari sembilan peserta:

1. Pendekatan berbasis regresi/pseudoisasi (Landmark, Phillips, Coats Engineering): peserta menjalankan (atau memanfaatkan) hasil grid halus terlebih dahulu, lalu menyesuaikan kurva relative permeability pada grid kasar melalui regresi agar hasilnya cocok dengan grid halus. Pendekatan ini secara konsisten memberi hasil yang sangat dekat dengan solusi grid halus.
2. Pendekatan upscaling satu fasa murni (Chevron, TotalFinaElf, GeoQuest, Streamsim, Roxar): permeabilitas absolut di-upscale dengan berbagai metode (arithmetic-harmonic, harmonic-arithmetic, power-law, flow-based dengan kondisi batas berbeda), tanpa mengubah kurva relative permeability asli. Pendekatan ini menghasilkan sebaran hasil yang jauh lebih lebar, terutama untuk prediksi water cut dan laju sumur individual.

GeoQuest secara khusus menguji enam metode upscaling permeabilitas berbeda pada grid kasar yang sama (15x55x17), memungkinkan perbandingan “apel-ke-apel” murni dari pilihan metode. Hasilnya: kondisi batas linear pressure memberikan prediksi laju alir yang paling buruk (terlalu tinggi secara signifikan), sedangkan kondisi batas no-flow dan metode arithmetic-harmonic memberikan prediksi laju dan water cut yang cukup baik, dan metode harmonic-arithmetic memberikan prediksi tekanan rata-rata lapangan paling akurat.

Hasil dan Temuan Utama

Solusi grid halus dari berbagai simulator (baik solver konvensional maupun streamline) sangat sepakat satu sama lain untuk kedua model: ini penting karena berarti “acuan” yang digunakan untuk menilai kualitas upscaling memang solid. Untuk Model 2, kelompok pseudo-based secara konsisten paling akurat, terutama untuk memprediksi laju produksi dan water cut sumur individual (Producer 1). Kelompok upscaling satu fasa murni menunjukkan sebaran hasil yang jauh lebih besar, khususnya pada waktu terobosan air (water breakthrough) yang dapat bervariasi antara sekitar 200 hingga lebih dari 400 hari tergantung metode.

Satu temuan menarik: akurasi prediksi tekanan rata-rata lapangan ternyata tidak berkorelasi langsung dengan akurasi prediksi laju minyak/air. Karena compressibility sistem ini rendah, pembagian produksi antar sumur lebih ditentukan oleh rasio productivity index (PI) antar sumur, sementara nilai absolut tekanan lapangan justru sangat sensitif terhadap perhitungan PI absolut masing-masing sumur: dua hal yang ternyata dapat “salah bersama-sama” secara konsisten tanpa terlihat dari akurasi laju alir saja.

Uji tambahan dari Heriot-Watt University (menjalankan grid dari 30x55x17 hingga sekasar 5x11x6 hanya dengan upscaling satu fasa) menunjukkan bahwa penjarangan grid yang berlebihan tanpa pseudoisasi menghasilkan kesalahan besar: menegaskan bahwa ada batas seberapa kasar grid boleh dibuat sebelum akurasi benar-benar runtuh.

Temuan dan Implikasi

• Solusi grid halus dari simulator dan pendekatan yang berbeda-beda (finite difference maupun streamline) dapat sangat sepakat satu sama lain, memberikan keyakinan bahwa “kebenaran acuan” untuk menilai upscaling memang dapat dipercaya.
• Metode upscaling yang memanfaatkan informasi aliran dari grid halus (pseudoisasi/regresi relative permeability) secara konsisten memberi hasil paling akurat dibanding metode upscaling satu fasa murni.
• Di antara metode upscaling satu fasa “buta” (tanpa menggunakan info grid halus), pendekatan flow-based dengan kondisi batas no-flow terbukti paling andal untuk kasus di studi ini.
• Kondisi batas linear-pressure yang digunakan dalam upscaling ternyata dapat menjadi pilihan yang buruk untuk model tertentu (studi ini), meski di studi geologi lain kondisi batas ini pernah terbukti baik: menunjukkan bahwa tidak ada satu metode upscaling yang universal cocok untuk semua model.
• Prediksi laju alir sumur individual jauh lebih sensitif terhadap pilihan metode upscaling dibanding prediksi total produksi lapangan: variasi hasil antar-metode dapat sebesar variasi antar-peserta yang berbeda simulator sekalipun.
• Akurasi tekanan rata-rata lapangan bukan indikator yang baik untuk menilai kualitas prediksi laju minyak/air pada kasus compressibility rendah seperti ini, karena tekanan lebih dipengaruhi oleh estimasi productivity index sumur daripada kualitas upscaling permeabilitas itu sendiri.
• Penjarangan grid yang terlalu agresif tanpa penyesuaian kurva relative permeability (pseudoisasi) berisiko menghasilkan kesalahan prediksi yang besar, terutama saat ukuran blok mendekati skala pola sumur pada model lapangan penuh.
• Karena upscaling selalu melibatkan trade-off antara ukuran grid dan akurasi, pemilihan metode upscaling sebaiknya disesuaikan dengan apakah data hasil grid halus tersedia (untuk pseudoisasi) atau tidak (sehingga harus mengandalkan upscaling satu fasa saja).