Second Comparative Solution Project: A Three-Phase Coning Study

Ringkasan

Ini adalah studi banding kedua dari seri SPE Comparative Solution Project, kali ini fokus pada masalah “coning”: fenomena di mana gas atau air ditarik naik/turun membentuk kerucut menuju lubang sumur produksi akibat penurunan tekanan lokal di sekitar sumur. Sebelas perusahaan dan lembaga riset menjalankan simulator black-oil tiga fase mereka masing-masing pada satu penampang radial (radial cross section) dengan satu sumur produksi di tengah. Kasusnya sengaja dibuat menantang secara numerik: densitas minyak dan air yang hampir sama membuat zona transisi kapiler minyak-air meluas tinggi ke kolom minyak, ditambah variasi laju produksi yang tajam dan GOR larutan yang tidak lazim tinggi. Hasilnya, secara umum, simulator-simulator ini sepakat dengan cukup baik meskipun numerical method yang digunakan sangat beragam.

Signifikansi Studi

Setelah studi banding pertama (masalah reservoir 3D skala lapangan) mendapat sambutan sangat positif di industri (bahkan penulisnya, Aziz Odeh, diminta mempresentasikan ulang paper itu di hadapan SPE London Petroleum Section) muncul dorongan kuat untuk melanjutkan tradisi ini dengan kasus yang berbeda karakter. Studi kedua ini sengaja memilih masalah coning skala sumur tunggal (bukan skala lapangan) karena coning dikenal sebagai salah satu tantangan numerik tersulit dalam reservoir simulation: gradien saturasi dan tekanan yang tajam di sekitar sumur menuntut simulator memiliki skema time-stepping dan linierisasi yang benar-benar stabil.

Nilai penting studi semacam ini bagi industri adalah sebagai bentuk validasi silang antar simulator komersial dan riset yang beredar saat itu, sekaligus sebagai “kasus uji” standar yang dapat digunakan siapa pun untuk menguji simulator baru: bahkan penulis paper ini secara explicit menyarankan beberapa modifikasi pada data input yang dapat membuat soal ini lebih sulit lagi, untuk keperluan pengujian simulator di masa depan.

Latar Belakang Studi

• Jenis reservoir/masalah: coning tiga fase (minyak, air, gas) pada penampang radial (radial cross section) dengan satu sumur produksi di pusat, melibatkan gas coning, water coning, dan repressurisasi gas (gas repressuring).
• Diorganisir menyusul rekomendasi Aziz Odeh kepada SPE, dengan ketua program Khalid Aziz pada SPE Symposium on Reservoir Simulation 1982; dipublikasikan sebagai SPE 10489 di Journal of Petroleum Technology, Maret 1986.
• Sebelas perusahaan/lembaga berpartisipasi: ARCO Oil and Gas Co., Chevron Oil Field Research Co., D&S Research and Development Ltd., Franlab Consultant S.A., Gulf Research and Development Co., Harwell, Intercomp, McCord-Lewis Energy Services, J.S. Nolen & Assocs., Scientific Software Corp. (SSC), dan Shell Development Co.
• Model reservoir: grid radial dengan 10 blok radial dan 15 lapisan vertikal, sumur diselesaikan (completed) di dua blok bagian atas lapisan.
• Skenario produksi: laju produksi minyak sengaja divariasikan tajam dalam beberapa periode waktu (1.000 STB/D turun ke 100 STB/D lalu naik lagi ke 1.000 STB/D, dst.) selama total 900 hari, untuk menguji respons transien simulator.
• Yang dibandingkan: jumlah fluida awal di tempat (initial fluids in place), laju produksi minyak, water cut, GOR, tekanan alir dasar sumur (BHP), pressure drawdown, serta data performa komputasi (jumlah time step, jumlah iterasi nonlinier, waktu CPU, jenis komputer yang digunakan).

Konsep Kunci

• Coning: fenomena penarikan air dan/atau gas membentuk profil seperti kerucut menuju perforasi sumur akibat drawdown tekanan yang kuat di sekitar sumur, dapat menyebabkan air atau gas ikut terproduksi lebih awal dari perkiraan.
• Radial cross section: model reservoir dua dimensi berbentuk irisan radial (jari-jari dan kedalaman) di sekitar satu sumur, cocok untuk studi perilaku dekat-sumur seperti coning tanpa perlu memodelkan seluruh lapangan.
• Zona transisi kapiler: rentang kedalaman di mana saturasi fluida berubah bertahap akibat efek kapiler, bukan berubah tiba-tiba di kontak fluida; pada studi ini zona ini sengaja dibuat tebal karena densitas minyak dan air yang hampir sama.
• Gas repressuring: kondisi ketika tekanan di sekitar sumur naik kembali (misalnya karena laju produksi diturunkan) sehingga gas bebas yang sempat terbentuk larut kembali ke dalam minyak.
• Point-centered vs block-centered grid: dua cara berbeda menempatkan titik hitung dalam sel grid; keduanya dapat memberi hasil yang tampak berbeda kecuali dikoreksi secara matematis terhadap posisi titik hitung.
• Adaptive Implicit Method (AIM): metode dari Scientific Software Corp. yang menerapkan tingkat keimplisitan berbeda per sel grid tergantung kebutuhan lokal, sehingga tidak semua sel “dipaksa” dihitung penuh implicit.
• Stone’s method (metode Stone): salah satu model matematis populer untuk menghitung relative permeability tiga fase (minyak-air-gas) dari data dua fase yang lebih mudah diukur; digunakan oleh mayoritas peserta studi ini dengan variasi masing-masing.

Pembahasan

Deskripsi Kasus

Reservoir dimodelkan sebagai penampang radial dengan 10 blok radial (dari jari-jari sumur 0,25 ft hingga radius luar 2.050 ft) dan 15 lapisan vertikal dengan ketebalan, permeabilitas horizontal/vertikal, dan porositas yang berbeda-beda tiap lapisan. Reservoir mula-mula berada dalam kesetimbangan kapiler-gravitasi dengan kontak gas-minyak di kedalaman 9.035 ft dan kontak minyak-air di 9.209 ft. Karena densitas minyak dan air yang berdekatan, zona transisi kapiler minyak-air menjadi sangat tebal: inilah salah satu sumber kesulitan numerik utama. Sumur produksi ditempatkan di pusat radial, diselesaikan pada dua blok teratas. Jadwal produksi dibuat berubah-ubah tajam (naik-turun antara 100 dan 1.000, lalu 3.000 STB/D) sehingga simulator harus menangani transien tekanan dan saturasi yang signifikan, bukan kondisi kuasi-statis yang lebih mudah dikonvergensikan.

Keragaman Pendekatan Numerik

Paper ini mendeskripsikan sebelas simulator dengan pendekatan yang sangat bervariasi: dari yang fully implicit penuh (Chevron, D&S, Gulf) dengan solusi Newton-Raphson, sequential implicit yang menjamin mass balance (Harwell PORES), formulasi semi-implicit dengan bobot upstream (ARCO, McCord-Lewis), sampai Adaptive Implicit Method dari SSC yang secara dinamis memilih tingkat keimplisitan per sel. Metode penghitungan relative permeability tiga fase juga bervariasi: sebagian besar menggunakan variasi metode Stone (metode pertama atau kedua), sementara Chevron menggunakan interpolasi diagram ternary sendiri. Solver aljabar linier juga beragam: eliminasi Gauss langsung dengan D4 ordering, metode iteratif, hingga metode conjugate gradient terpotong (truncated conjugate gradient).

Hasil dan Titik Perbedaan

Secara umum, penulis melaporkan tingkat kesesuaian yang tinggi meskipun discretization scheme dan metode solusi yang digunakan beragam. Kurva laju produksi minyak, water cut, GOR, dan BHP dari sebelas peserta hampir berimpit selama 900 hari simulasi. Penyimpangan utama terdapat pada kurva pressure drawdown milik Gulf yang menggunakan point-centered grid, sedangkan sebagian besar peserta menggunakan block-centered grid. Perbedaan posisi titik hitung pertama, yaitu 2,0 ft dibandingkan 0,84 ft dari sumur, menghasilkan drawdown yang tampak lebih tinggi. Setelah dilakukan koreksi logaritmik terhadap posisi titik hitung, hasil Gulf menjadi konsisten dengan hasil peserta lain. Hal ini menegaskan bahwa perbedaan konvensi numerik dapat menghasilkan deviasi tanpa menunjukkan perbedaan pada model fisik.

Paper juga mencatat data performa komputasi (jumlah time step, iterasi nonlinier, waktu CPU) dari tiap peserta yang dijalankan di komputer berbeda-beda (IBM 3033, Cray-1S, VAX 11/780, CDC, Univac, dll.), meskipun penulis mengakui data ini sulit ditarik kesimpulan umum karena ukuran soal yang relatif kecil dan perbedaan mesin.

Catatan untuk Perancang Uji Simulator

Bagian menarik lain dari paper ini adalah daftar saran explicit tentang cara membuat soal ini lebih sulit lagi untuk keperluan pengujian simulator baru: misalnya mengubah rasio $k_v/k_h$ untuk memperbesar efek coning, menaikkan laju produksi untuk memperbesar transien tekanan, atau mengubah $kh$ sumur pada lapisan ketat (tight layer) untuk membuat alokasi produksi antar lapisan lebih bergantung pada BHFP.

Temuan dan Implikasi

• Coning adalah salah satu masalah paling menantang secara numerik dalam reservoir simulation karena melibatkan gradien tajam di dekat sumur: cocok dijadikan kasus uji stabilitas dan konvergensi.
• Meskipun numerical method sangat beragam (fully implicit, sequential implicit, semi-implicit, adaptive implicit), hasil sebelas simulator pada studi ini secara umum sepakat dengan baik.
• Perbedaan hasil tidak selalu berarti ada yang salah secara fisik: kasus Gulf menunjukkan bahwa konvensi grid (point-centered vs block-centered) dapat menghasilkan angka yang tampak berbeda padahal secara matematis konsisten setelah dikoreksi.
• Zona transisi kapiler yang tebal (akibat densitas minyak-air yang berdekatan) dan variasi laju produksi yang tajam adalah dua “pendekatan” desain yang efektif untuk memaksa simulator bekerja pada kondisi sulit.
• Data performa komputasi (waktu CPU, jumlah time step) sulit dibandingkan secara adil ketika peserta menggunakan mesin yang sangat berbeda: perbandingan performa lintas simulator butuh kehati-hatian ekstra.
• Studi banding yang baik tidak berhenti pada “siapa yang benar”, namun juga menghasilkan rekomendasi konkret (di sini: cara mempersulit kasus uji) yang berguna bagi pengembang simulator berikutnya.
• Rangkaian SPE Comparative Solution Project terbukti memiliki nilai berkelanjutan bagi industri, terbukti dari permintaan presentasi ulang dan dorongan aktif dari komunitas untuk melanjutkan studi serupa dengan kasus yang berbeda.