Reservoir Grid and Geometric Models

Halaman pembuka modul Grid Model — modul terpanjang dan terakhir dalam rangkaian TM-6011 Advanced Reservoir Simulation.

Modul ini membahas satu pertanyaan besar yang belum disentuh modul-modul sebelumnya: bagaimana reservoir sesungguhnya dipetakan menjadi grid (jaringan gridblock), pilihan-pilihan geometri apa yang tersedia, konsekuensi numeriknya masing-masing (dispersi numerik, orientasi grid), bagaimana sumur dimodelkan di dalam satu gridblock yang jauh lebih besar dari radius sumur sesungguhnya (model Peaceman), dan akhirnya bagaimana semua ini disatukan dalam skema solusi Fully Implicit dan IMPES — menutup seluruh alur derivasi persamaan aliran yang dimulai dari Modul 2.

Aziz (SPE 25233) mendaftar pertimbangan-pertimbangan dalam memilih pendekatan gridding:

• Geologi dan ukuran reservoir
• Tipe proses perpindahan fluida (displacement process)
• Lokasi dan tipe sumur
• Akurasi numerik yang diinginkan
• Software yang tersedia
• Tujuan studi simulasi
• Kompetensi insinyur/tim simulasi
• Sumber daya komputer yang tersedia

Daftar ini menjadi kerangka pikir untuk seluruh sisa modul: setiap pilihan gridding yang dibahas berikutnya (block-centered, point-distributed, corner point, Voronoi, dsb.) adalah jawaban berbeda terhadap kombinasi pertimbangan-pertimbangan ini.

Simulator modern menawarkan beragam pilihan gridding:

• Local grid refinement — penghalusan grid lokal
• Hybrid grids — kombinasi dua atau lebih tipe grid
• Curvilinear grid (streamline grid) — grid mengikuti garis alir
• Voronoi/PEBI grid
• Corner Point Grid
• Dynamic Grid
• Automatic Grid Generation

Halaman ini adalah daftar isi visual untuk sisa modul: setiap opsi di atas akan dibahas satu per satu di halaman-halaman berikutnya, lengkap dengan kelebihan dan kekurangannya masing-masing.

Orientasi grid relatif terhadap posisi sumur dapat memengaruhi hasil simulasi.

Fluida bergerak secara preferensial mengikuti garis-garis grid. Ilustrasi membandingkan dua tata letak grid berukuran sama: grid persegi biasa (kiri) versus grid diagonal (kanan), masing-masing dengan sumur injeksi (titik merah) dan sumur produksi (titik hijau) di posisi yang sama.

Jika titik merah adalah sumur injeksi air, air akan breakthrough (tembus) paling cepat pada grid diagonal (kanan) — meskipun dimensi grid sama persis. Ini adalah artefak numerik murni, bukan fenomena fisik: arah garis grid yang kebetulan sejajar dengan garis lurus antar-sumur mempercepat transmisi numerik di sepanjang arah tersebut.

Melanjutkan halaman sebelumnya, ilustrasi dari Mattax and Dalton (Fig. 5.20) membandingkan dua cara sebuah gridblock “melihat” tetangganya:

(a) Formulasi five-point — aliran hanya dihitung ke 4 tetangga langsung (utara, selatan, timur, barat), digambarkan sebagai 4 anak panah tegak lurus dari titik pusat.

(b) Formulasi nine-point — aliran juga dihitung ke 4 tetangga diagonal, sehingga total 8 arah keluar dari titik pusat (mirip papan catur, gerakan raja).

Formulasi nine-point mengurangi efek grid orientation dari halaman sebelumnya (karena aliran diagonal ikut diperhitungkan secara explicit, tidak hanya “bocor” melalui zig-zag sepanjang garis grid horizontal/vertikal), tetapi menambah kompleksitas komputasi karena setiap gridblock kini terhubung ke 8 tetangga, bukan 4.

Masalah numerik kedua setelah grid orientation: dispersi numerik (numerical dispersion) — penyamaran (smearing) front fluida yang seharusnya tajam.

• Secara fisik, pada waterflood 1D, front air seharusnya maju sebagai batas tajam yang bergerak satu sel per timestep.
• Masalah muncul karena simulator tidak dapat merepresentasikan distribusi saturasi di dalam satu gridblock — setiap gridblock hanya menyimpan satu nilai saturasi rata-rata, bukan profil kontinu.

Ini adalah konsekuensi langsung dari diskretisasi (Modul 5): setiap gridblock adalah satu titik data, sehingga transisi tajam antar dua kondisi fisik yang berbeda dalam satu blok “tersebar” menjadi nilai rata-rata tunggal — diilustrasikan secara kuantitatif di dua halaman berikutnya.

Ilustrasi kuantitatif pertama dari dispersi numerik, pada waktu $t = t_1$: lima gridblock berjajar, front air (batas antara $S_w=1.0$ dan $S_w=0$) berada di tengah blok ke-4.

True Solution (solusi sejati): front tajam — tiga blok pertama penuh air ($S_w=1.0$), blok keempat terbagi tajam menjadi sebagian air dan sebagian minyak, blok kelima penuh minyak ($S_w=0$).

Simulator: karena tiap gridblock hanya dapat menyimpan satu nilai saturasi, blok keempat yang sesungguhnya terbagi tajam direpresentasikan sebagai satu nilai rata-rata $S_w=0.4$ — front yang sejatinya tajam sudah mulai “kabur” menjadi nilai tunggal yang lebih rendah dari kondisi sesungguhnya di sebagian blok tersebut.

Melanjutkan halaman sebelumnya, satu timestep kemudian ($t = t_1 + \Delta t$):

True Solution: front tajam sudah bergerak maju — tiga blok pertama masih penuh air ($S_w=1.0$), blok keempat kini penuh minyak ($S_w=0$), blok kelima juga minyak ($S_w=0$). Front berpindah rapi satu blok ke depan tanpa “kabur”.

Simulator: karena blok keempat pada langkah sebelumnya sudah terlanjur “kabur” menjadi $S_w=0.4$ (bukan tajam), hasil simulator pada langkah ini justru semakin melebar: blok ketiga tetap $S_w=1.0$, blok keempat menjadi $S_w=0.7$, blok kelima $S_w=0.2$ — front yang seharusnya satu batas tajam kini tersebar ke dua blok sekaligus, keduanya bernilai antara.

Ini menunjukkan sifat kumulatif dispersi numerik: kesalahan smearing dari satu timestep terbawa dan semakin melebar di timestep berikutnya — semakin banyak timestep dijalankan, semakin lebar zona transisi buatan ini dibandingkan front sejati yang tajam.

Setelah membahas dua masalah numerik (grid orientation, dispersi numerik), modul beralih ke jenis-jenis geometri grid itu sendiri, dimulai dari yang paling umum.

Block-Centered Grid: reservoir dibagi menjadi blok-blok, dan node (titik representasi) diletakkan di pusat setiap blok (SPE 25233). Diagram menunjukkan tata letak blok $(i,j)$ dengan node ditandai titik di tengah masing-masing kotak, diberi label sesuai indeks $(i-1,j)$, $(i,j+1)$, $(i+1,j)$, $(i,j-1)$, dst.

Batas antar-blok berada tepat di tengah jarak antar dua node bertetangga — pendekatan yang konsisten dengan seluruh penurunan mass balance per gridblock di Modul 5.

Dua poin singkat menutup pembahasan block-centered grid:

• Block-centered grid adalah yang paling umum digunakan di industri.
• Grid ini merupakan pendekatan yang “familiar” bagi para insinyur — sesuai intuisi fisik langsung (satu blok = satu titik data pusat, sama seperti yang dibayangkan pada Modul 5).

Kelaziman ini adalah alasan mengapa seluruh contoh numerik di Modul 5 (Oil Material Balance) secara implicit mengasumsikan block-centered grid.

Alternatif dari block-centered grid: Point-Distributed Grid.

Pada pendekatan ini, urutannya dibalik: node diletakkan lebih dahulu, baru kemudian gridblock dikonstruksi di sekitarnya (SPE 25233). Diagram menunjukkan node-node ($i-1,j$; $i,j+1$; $i+1,j$; $i,j-1$; dst.) sebagai titik-titik yang sudah ditentukan posisinya, dengan garis-garis batas blok (garis tebal) yang mengikuti letak node tersebut, bukan sebaliknya.

Perbedaan urutan konstruksi ini (node dulu vs. blok dulu) berdampak pada akurasi relatif suku akumulasi vs. suku aliran, sebagaimana dibahas di halaman perbandingan berikutnya.