Reservoir Grid and Geometric Models

Ringkasan

Modul ini membahas bagaimana reservoir dipetakan menjadi grid (kumpulan blok-blok kecil) untuk keperluan simulasi numerik, mulai dari pertimbangan memilih tipe grid, berbagai jenis grid modern (corner point, Voronoi/PEBI, hybrid, dan lainnya), sampai efek pilihan grid terhadap akurasi hasil simulasi. Selain itu, modul ini juga membahas dua topik pendukung yang erat kaitannya dengan grid: bagaimana sumur dimodelkan di dalam satu grid block (Peaceman well index), dan bagaimana persamaan governing yang dibahas di modul-modul sebelumnya benar-benar diselesaikan lewat skema explicit, implicit, atau IMPES.

Signifikansi

Grid adalah fondasi dari seluruh proses reservoir simulation: sebelum satu pun persamaan aliran fluida dapat dihitung, reservoir harus dulu dipecah menjadi blok-blok diskrit tempat persamaan tersebut diterapkan. Pilihan tipe grid, ukuran blok, dan orientasinya terhadap sumur bukan sekadar detail teknis; pilihan yang buruk dapat membuat hasil simulasi menyimpang jauh dari kondisi fisik yang sesungguhnya, misalnya lewat fenomena grid orientation effect dan numerical dispersion yang dibahas di modul ini.

Modul ini juga menyatukan banyak benang dari modul-modul sebelumnya (konsep grid block dan mass balance dari Modul 5, serta metode solusi seperti IMPES dari modul sebelumnya, semuanya muncul kembali di sini dalam konteks yang lebih luas) termasuk bagaimana sumur (yang berukuran jauh lebih kecil dari satu grid block) dimodelkan secara konsisten dengan persamaan-persamaan tersebut.

Konsep Kunci

• Block-centered grid: grid di mana titik hitung (node) diletakkan di tengah setiap blok; pendekatan yang paling umum dan familiar bagi engineer.
• Point-distributed grid: grid di mana titik-titik node ditentukan lebih dulu, lalu batas-batas blok dibangun mengelilingi titik tersebut.
• Corner point grid: grid dengan bentuk blok yang benar-benar tidak beraturan (bukan kotak), tetapi tetap “logically rectangular” karena setiap blok didefinisikan oleh 8 titik sudut; digunakan untuk menangkap geometri reservoir yang kompleks.
• Voronoi/PEBI grid: grid di mana setiap blok adalah wilayah yang lebih dekat ke satu titik grid dibanding titik grid lain mana pun; generalisasi dari point-distributed grid.
• Grid orientation effect: fenomena di mana orientasi garis grid terhadap posisi sumur memengaruhi hasil simulasi (misalnya waktu breakthrough air), meski ukuran grid sama.
• Numerical dispersion (smearing): front saturasi yang tajam secara fisik jadi tampak “melebar”/menghalus di hasil simulasi, karena satu grid block hanya dapat memiliki satu nilai saturasi rata-rata.
• Peaceman well index / equivalent radius ($r_o$): cara menghubungkan tekanan hasil hitungan simulator di satu grid block dengan tekanan aliran nyata di dinding sumur (wellbore), lewat radius ekuivalen $r_o$.
• Skema explicit vs implicit: dua cara mengevaluasi suku aliran antar blok pada waktu perhitungan (waktu lama vs waktu baru), yang menentukan tingkat akurasi dan kestabilan numerik.

Pembahasan

Pertimbangan dalam Memilih Tipe Grid

Menurut Aziz (SPE 25233), ada banyak faktor yang menentukan tipe grid yang tepat untuk suatu studi simulasi: geologi dan ukuran reservoir, jenis proses displacement (misalnya waterflood atau injeksi gas), lokasi dan jenis sumur, tingkat akurasi numerik yang diinginkan, software yang tersedia, tujuan studi simulasi, kompetensi tim engineer, hingga sumber daya komputer yang tersedia. Tidak ada satu tipe grid yang selalu terbaik: pemilihannya selalu berupa trade-off antara akurasi, kompleksitas, dan biaya komputasi.

Simulator modern umumnya menyediakan beberapa opsi: local grid refinement, hybrid grid, curvilinear/streamline grid, Voronoi/PEBI grid, corner point grid, dynamic grid, dan automatic grid generation. Masing-masing memiliki kelebihan dan kelemahan yang dibahas satu per satu di bawah.

Efek Orientasi Grid dan Numerical Dispersion

Salah satu temuan penting (Mattax and Dalton) adalah bahwa orientasi garis grid terhadap posisi sumur dapat memengaruhi hasil simulasi. Fluida cenderung bergerak lebih mudah mengikuti arah garis grid. Contohnya, jika ada sumur injeksi air, waktu breakthrough air di sumur produksi dapat berbeda antara konfigurasi grid yang “sejajar” dan yang “diagonal” terhadap sumur: meskipun ukuran gridnya sama. Ini adalah artefak numerik, bukan fenomena fisik yang sesungguhnya terjadi di reservoir.

Permasalahan kedua adalah numerical dispersion atau smearing. Pada waterflood 1-D, displacement front yang tajam secara fisis bergerak maju dari satu blok ke blok berikutnya dengan saturasi air bernilai 1 di belakang front dan bernilai awal di depannya. Karena satu grid block hanya merepresentasikan satu nilai saturasi rata-rata, simulator tidak dapat menangkap distribusi saturasi subgrid. Akibatnya, front tampak semakin menyebar seiring bertambahnya jarak rambat dibandingkan dengan solusi acuan.

Block-Centered vs Point-Distributed Grid

Ada dua filosofi dasar dalam menata grid. Pada block-centered grid, reservoir dibagi menjadi blok-blok dulu, lalu titik hitung (node) diletakkan di tengah tiap blok: ini adalah pendekatan yang paling umum digunakan dan paling familiar bagi engineer. Pada point-distributed grid, urutannya dibalik: titik-titik node ditentukan lebih dulu, lalu batas blok dibangun mengelilingi tiap titik.

Perbedaan filosofi ini berdampak pada akurasi: block-centered grid memberikan perhitungan suku akumulasi yang lebih rigorous, sementara point-distributed grid lebih akurat untuk suku aliran (flow term). Namun jika grid yang digunakan reguler (semua blok berukuran sama) atau mendekati reguler, perbedaan antara kedua metode ini menjadi tidak signifikan.

Grid Tidak Beraturan dan Local Refinement

Untuk problem simulasi besar, grid yang halus sesungguhnya hanya dibutuhkan di bagian reservoir tempat tekanan atau saturasi berubah cepat (misalnya di sekitar sumur atau front saturasi). Menggunakan grid ireguler standar di seluruh reservoir cenderung menghasilkan blok-blok kecil yang tidak perlu di bagian-bagian lain, sehingga memboroskan sumber daya komputasi.

Local grid refinement mengatasi ini dengan membagi satu sel (atau sekelompok sel) menjadi banyak sel yang lebih kecil hanya di area yang diperlukan. Masalah utamanya adalah akurasi perhitungan aliran antara blok kasar dan blok halus di batas antara kedua area tersebut. Hybrid grid mengambil pendekatan berbeda: menggabungkan dua atau lebih tipe grid berbeda menjadi satu model. Ini membantu menangkap profil saturasi tajam di dekat sumur dan meningkatkan akurasi perhitungan water-oil ratio (WOR) dan gas-oil ratio (GOR), tetapi konsekuensinya struktur matriks pita (banded) dari sistem persamaan menjadi hilang, sehingga proses penyelesaian matriks menjadi lebih rumit dan kurang efisien.

Corner Point, Voronoi/PEBI, Streamline, dan Dynamic Grid

Corner point grid memungkinkan bentuk blok yang benar-benar tidak beraturan (bukan kotak), namun tetap “logically rectangular” karena setiap blok didefinisikan oleh 8 titik sudut. Grid tipe ini digunakan untuk menangkap geometri reservoir yang kompleks (misalnya mengikuti struktur patahan atau lapisan yang melengkung). Konsekuensinya, aliran lintas satu muka blok dapat bergantung pada lebih dari dua tekanan di kedua sisi muka tersebut: dan ketika grid sangat skewed (miring), hubungan antar blok tidak lagi orthogonal terhadap muka bloknya, yang menambah kompleksitas perhitungan.

Grid Voronoi (juga disebut PEBI, “perpendicular bisector”) adalah konsep lama yang pertama didefinisikan tahun 1908, dan penggunaannya dalam reservoir simulation dipelopori oleh Dr. Heinemann di University of Leoben. Satu blok grid Voronoi didefinisikan sebagai “wilayah ruang yang lebih dekat ke titik gridnya sendiri dibanding ke titik grid lain mana pun”: garis yang menghubungkan dua titik grid tetangga selalu tegak lurus terhadap batas blok di antara mereka. Grid ini adalah generalisasi dari point-distributed grid.

Pada streamline grid, urutan pembuatannya justru dibalik: lintasan streamline dihitung dulu (dengan menyelesaikan persamaan tekanan IMPES), lalu corner point grid dibangun mengikuti lintasan streamline tersebut. Karena metode finite difference standar tetap digunakan untuk menyelesaikan persamaan aliran, fluida masih dapat bergerak melintasi streamline (tidak terkunci mengikutinya secara sempurna). Terakhir, dynamic grid adalah pendekatan di mana grid diubah selama simulasi berjalan: sel-sel diperhalus atau dikasarkan tergantung kecepatan aliran fluida di titik dan waktu tersebut.

Persamaan Governing dan Skema Explicit vs Implicit

Setelah grid ditentukan, langkah berikutnya adalah menuliskan persamaan mass balance di tiap grid block (pola yang sama dengan yang sudah dibahas di Modul 5) suku akumulasi (perubahan $V_p \cdot S/B$ antara waktu baru dan lama, dibagi $\Delta t$) harus sama dengan jumlah suku aliran (net flow rate) dari semua blok tetangga.

Suku aliran ini dapat dievaluasi dengan dua skema: skema explicit menghitung transmissibility dan mobility pada waktu lama ($n$), sementara skema implicit menghitungnya pada waktu baru ($n+1$). Mobility dievaluasi menggunakan upstream weighting (nilai relative permeability diambil dari blok yang menjadi sumber aliran), dan formation volume factor serta viskositas antar blok dihitung sebagai rata-rata dari kedua blok yang berdekatan.

Pemodelan Sumur: Peaceman Well Index

Sumur menimbulkan tantangan khusus karena radiusnya jauh lebih kecil dibanding ukuran satu grid block: akibatnya, tekanan yang dihitung simulator untuk blok tempat sumur berada bukanlah tekanan aliran sesungguhnya di dinding sumur (wellbore). Selisih antara keduanya dihubungkan lewat sebuah koefisien yang disebut Productivity/Injectivity Index.

Peaceman menemukan bahwa tekanan yang dihitung untuk suatu wellblock ternyata sama dengan tekanan aliran pada suatu radius tertentu di sekitar sumur, yang disebut radius ekuivalen ($r_o$): didefinisikan sebagai “radius di mana tekanan aliran steady-state dari sumur yang sesungguhnya sama dengan tekanan hasil hitungan numerik untuk blok sumur tersebut”. Untuk reservoir yang seragam dan isotropik dengan grid persegi, $r_o = 0.28 \Delta x$; untuk grid non-uniform isotropik, $r_o = 0.14\sqrt{\Delta x^2 + \Delta y^2}$; dan untuk reservoir anisotropik ada rumus tambahan yang melibatkan rasio permeabilitas arah x dan y.

Dari radius ekuivalen ini disusun Well Index ($WI$), yang menghubungkan laju alir sumur (untuk fase minyak, air, atau gas) dengan selisih antara tekanan grid block dan tekanan dasar sumur (BHP):

$WI = \frac{2\pi k h}{\ln(r_o/r_w) + s}$

di mana $r_w$ adalah radius sumur sesungguhnya dan $s$ adalah faktor skin. Sumur produksi dapat dioperasikan dengan beberapa kendala (constraint): laju alir konstan, BHP konstan, atau drawdown konstan; sementara sumur injeksi dapat dioperasikan dengan laju injeksi konstan atau tekanan injeksi konstan.

Kaitan dengan Metode Solusi: Fully Implicit dan IMPES

Modul ini juga menyinggung bagaimana persamaan-persamaan governing di atas benar-benar diselesaikan. Salah satu pendekatan adalah menyusunnya sebagai fungsi residual (selisih antara nilai kiri dan kanan mass balance untuk minyak, gas, air) yang harus disamakan dengan nol: persoalan ini setara dengan persoalan root-finding, yang dapat diselesaikan lewat metode iteratif seperti Newton-Raphson atau Quasi-Newton, dengan menghitung matriks Jacobian dari sistem persamaan tersebut di setiap iterasi. Pendekatan inilah yang mendasari metode fully implicit.

Sebagai alternatif yang lebih efisien secara komputasi, ada metode IMPES (Implicit Pressure, Explicit Saturation) yang sudah dibahas secara khusus di modul sebelumnya: ketiga persamaan fase digabung menjadi satu persamaan tekanan yang diselesaikan implicit, lalu saturasi dihitung explicit setelahnya. Munculnya kembali IMPES di modul ini menegaskan bahwa pilihan grid, pemodelan sumur, dan metode solusi numerik semuanya adalah bagian dari satu sistem yang saling terkait dalam reservoir simulation.

Kesimpulan Utama

• Pemilihan tipe grid selalu berupa trade-off antara akurasi, kompleksitas, dan biaya komputasi, dan dipengaruhi banyak faktor (geologi, jenis displacement, lokasi sumur, dst).
• Grid orientation effect menunjukkan bahwa orientasi garis grid terhadap sumur dapat memengaruhi hasil simulasi meski ukuran grid sama.
• Numerical dispersion membuat front saturasi yang tajam secara fisik tampak melebar di hasil simulasi karena keterbatasan representasi per grid block.
• Block-centered grid lebih umum digunakan; point-distributed grid lebih akurat untuk suku aliran namun kurang rigorous untuk suku akumulasi.
• Corner point grid dan Voronoi/PEBI grid digunakan untuk menangkap geometri reservoir yang kompleks, dengan trade-off kompleksitas perhitungan aliran antar blok.
• Local refinement dan hybrid grid memperhalus grid di area yang butuh akurasi tinggi, namun dapat merusak struktur matriks yang efisien.
• Peaceman well index menghubungkan tekanan grid block dengan tekanan dasar sumur lewat konsep radius ekuivalen ($r_o$).
• Persamaan governing dapat diselesaikan lewat metode fully implicit (Newton-Raphson) yang lebih akurat namun berat, atau IMPES yang lebih ringan namun memiliki keterbatasan tersendiri.