Ninth SPE Comparative Solution Project: A Reexamination of Black-Oil Simulation

Ringkasan

Paper ini adalah studi perbandingan kesembilan dalam seri SPE Comparative Solution Project, ditulis oleh J.E. Killough dari University of Houston dan dipresentasikan pada 1995. Setelah delapan studi sebelumnya banyak membahas topik-topik khusus (kompositional, dual porosity, thermal, gridding), studi ini “kembali ke dasar” dengan menguji ulang simulasi black-oil biasa, tetapi kali ini pada model berukuran sedang (9.000 sel) dengan tingkat heterogenitas tinggi karena medan permeabilitasnya dibangun secara geostatistik (acak namun terkorelasi secara statistik, meniru kondisi reservoir nyata). Sembilan perusahaan/tim menjalankan kasus yang identik dan membandingkan hasilnya.

Signifikansi Studi

Sejak dua studi comparative pertama di awal 1980-an, kemampuan simulator black-oil sudah berkembang jauh: metode solver linear lebih canggih, penanganan sumur lebih sofistikasi, dan penggunaan solusi fully implicit semakin umum untuk kasus-kasus sulit. Studi ini penting karena memberikan “titik acuan terbaru” (up-to-date benchmark) untuk memvalidasi generasi simulator black-oil yang lebih baru, sekaligus secara khusus menguji bagaimana simulator-simulator itu menangani tantangan heterogenitas permeabilitas gaya geostatistik: sesuatu yang semakin umum digunakan di industri untuk memodelkan reservoir nyata, tetapi dapat menimbulkan kesulitan numerik tersendiri (misalnya pada konvergensi Newton-Raphson).

Latar Belakang Studi

• Kasus yang disimulasikan: reservoir miring (dip 10 derajat) berlapis 15 layer dengan grid 24x25x15 (9.000 sel), memiliki medan permeabilitas heterogen yang dibangun secara geostatistik dengan panjang korelasi sekitar 1.800 kaki (enam blok grid) di satu arah horizontal saja.
• Diorganisir dan dilaporkan oleh J.E. Killough (University of Houston); sembilan peserta memberikan data: AEA Technology (TechSIM), ARCO (simulator in-house), CMG (IMEX), Intera (ECLIPSE 100), Fina, SENSOR, SSI (SIMBEST II), TIGRESS, dan Western Atlas Software (VIP/DESKTOP-VIP).
• Skenario uji: 25 sumur produksi minyak yang tersebar cukup acak, dengan satu sumur injeksi air; laju produksi maksimum diubah tiga kali sepanjang simulasi 900 hari untuk memicu perubahan tekanan dan pembentukan tudung gas sekunder.
• Data PVT dan kurva relative permeability mengambil dari studi comparative kedua (Second CSP), termasuk kurva tekanan kapiler air-minyak yang memiliki diskontinuitas tidak biasa pada saturasi air sekitar 35%.

Konsep Kunci

• Black-oil simulation: pemodelan reservoir dengan asumsi tiga fase (minyak, gas, air) yang sifatnya hanya bergantung pada tekanan, tanpa perlu menghitung komposisi kimia detail seperti pada simulasi kompositional.
• Medan permeabilitas geostatistik: permeabilitas antar sel dibuat “acak” tetapi tetap saling berkorelasi secara statistik (diukur lewat semi-variogram) sehingga meniru pola heterogenitas batuan sungguhan, bukan angka seragam.
• IMPES vs fully implicit: dua strategi penyelesaian numerik; IMPES (Implicit Pressure, Explicit Saturation) menghitung tekanan secara implicit namun saturasi secara explicit sehingga langkah waktunya harus sangat kecil pada kasus yang tidak stabil, sedangkan fully implicit menghitung semuanya sekaligus sehingga jauh lebih stabil meski satu langkahnya memerlukan sumber daya komputasi yang lebih besar.
• Newton Iteration: proses iteratif yang digunakan simulator fully implicit untuk menyelesaikan persamaan aliran nonlinear pada setiap langkah waktu; semakin sulit konvergensinya, semakin banyak iterasi dibutuhkan.
• Well index / productivity index (PI): angka yang menghubungkan laju alir sumur dengan perbedaan tekanan antara reservoir dan lubang sumur; cara menghitungnya dapat berbeda antar simulator dan berdampak besar pada laju produksi sumur individual.
• Stone’s Model I vs Model II: dua pendekatan berbeda untuk menghitung relative permeability tiga fase (minyak-air-gas sekaligus); pemilihan model ini dapat mengubah signifikan laju produksi air yang diprediksi.
• Gas Oil Ratio (GOR): indikator munculnya gas bebas di reservoir; di studi ini digunakan untuk menunjukkan seberapa sensitif hasil terhadap ukuran langkah waktu (time step).

Pembahasan

Deskripsi Kasus

Reservoirnya berupa lapisan miring dengan 15 layer yang porositas dan ketebalannya bervariasi per lapisan (dari 8% hingga 17% porositas, ketebalan 8 hingga 100 kaki). Data PVT minyak dan gas, serta fungsi saturasi gas-minyak, diambil ulang dari studi comparative kedua. Kondisi awal reservoir memiliki kontak minyak-air pada 9.950 kaki subsea, dan karena kurva tekanan kapiler tidak mencapai saturasi air 100%, seluruh model memiliki sedikit saturasi minyak sisa bahkan di zona air. Yang membuat kasus ini unik dibanding studi black-oil sebelumnya adalah medan permeabilitas yang dihasilkan secara geostatistik: didistribusikan sel per sel kepada peserta, lengkap dengan semi-variogramnya, sehingga semua peserta memodelkan heterogenitas yang identik.

Dua puluh lima sumur produksi diberi batas laju maksimum 1.500 STB/hari yang diturunkan menjadi 100 STB/hari pada hari ke-300, lalu dinaikkan kembali ke 1.500 STB/hari pada hari ke-360 hingga akhir simulasi (900 hari): skenario ini sengaja dirancang untuk memicu penurunan tekanan signifikan dan pembentukan tudung gas sekunder di puncak reservoir yang sangat dipengaruhi oleh pola heterogenitas permeabilitas.

Metode dan Parameter yang Dibandingkan

Analisis awal menunjukkan sensitivitas yang tinggi terhadap pemilihan numerical scheme. Pada simulator yang sama, skema IMPES memerlukan sekitar 2.180 time step, sedangkan skema fully implicit hanya memerlukan 27 time step dengan total 109 Newton iteration. Perbedaan tersebut berkaitan dengan ketidakstabilan IMPES ketika gas mulai mengalami percolation menuju bagian atas reservoir. Perbandingan time step 60 hari, 10 hari, dan skema IMPES menghasilkan variasi GOR lapangan hampir 10%. Dengan demikian, sebagian variasi antarpeserta berasal dari konfigurasi numerik, bukan dari perbedaan physical model.

Hasil dan Temuan Utama

Secara keseluruhan, tingkat kesepakatan antar sembilan peserta cukup baik pada level lapangan (field-wide): laju produksi minyak total bervariasi dalam rentang sekitar 9% dari nilai rata-rata, laju gas sekitar 11%, tekanan reservoir hanya berbeda sekitar 4%, dan saturasi gas di tudung gas sekunder berbeda sekitar 5%. Namun ada beberapa titik di mana hasil menyebar cukup lebar:

• Laju air: variasi mencapai sekitar 20%, kemungkinan besar karena perbedaan pendekatan relative permeability tiga fase (Stone’s Model I vs Model II) dan asumsi inisialisasi saturasi air di zona akuifer (sebagian peserta menggunakan saturasi air 100% di area akuifer, yang cenderung menaikkan produksi air).
• Saturasi air di dekat baris sumur produksi: variasi sekitar 25% menjelang akhir simulasi, sekali lagi terkait perlakuan kurva relative permeability dan tekanan kapiler, serta laju injeksi air yang berbeda-beda akibat batasan tekanan dasar sumur injeksi.
• Laju produksi minyak sumur 21 (well individu) (dua peserta menunjukkan laju yang jauh lebih tinggi dari peserta lainnya, dan penulis menduga ini disebabkan perbedaan cara menghitung well index (productivity index)) peserta lainnya justru sepakat dalam rentang sekitar 4% saja.

Dari sisi performa komputasi, seluruh peserta menyelesaikan simulasi dalam kurang dari 60 langkah waktu dengan rata-rata 4-5 Newton iteration per langkah, dan pada workstation kelas atas saat itu (IBM RS/6000), simulasi ini dapat diselesaikan hanya dalam hitungan menit: menegaskan bahwa masalah berukuran sedang seperti ini sudah cukup ringan untuk komputasi rutin di pertengahan 1990-an.

Temuan dan Implikasi

• Untuk kasus reservoir heterogen berukuran sedang seperti ini, simulator black-oil modern (saat itu) dapat mencapai kesepakatan yang cukup kuat pada level total produksi lapangan: variasi umumnya di bawah 10%.
• Variasi hasil yang jauh lebih besar justru muncul pada level sumur individual dan pada parameter-parameter yang sensitif terhadap asumsi relative permeability/tekanan kapiler tiga fase, seperti laju air dan saturasi air.
• Perbedaan treatment well index/productivity index antar simulator dapat menghasilkan perbedaan laju produksi sumur individual yang sangat besar (hingga 30% menurut abstrak), meskipun hasil keseluruhan lapangan tetap sepakat.
• Skema time-stepping (IMPES vs fully implicit) bukan sekadar soal kecepatan komputasi: pada kasus dengan gas percolation dan heterogenitas tinggi, IMPES dapat menjadi sangat tidak stabil dan butuh langkah waktu jauh lebih kecil.
• Sebagian variasi hasil antar-peserta ternyata bukan berasal dari perbedaan physical model sama sekali, melainkan murni dari pilihan pengaturan numerik (ukuran langkah waktu): penting diingat saat membandingkan hasil simulator yang berbeda.
• Asumsi inisialisasi kondisi awal (misalnya saturasi air di akuifer) dapat menghasilkan perbedaan hasil produksi air yang signifikan meski model geologinya identik.
• Medan permeabilitas heterogen bergaya geostatistik menghasilkan pola aliran fluida yang jauh lebih rumit dan tidak intuitif (misalnya distribusi saturasi gas yang “berliku”) dibanding model homogen sederhana.
• Studi ini menjadi acuan validasi yang relevan untuk simulator generasi baru karena menyediakan kasus black-oil heterogen dengan data lengkap dan hasil pembanding dari sembilan simulator komersial/riset berbeda.