PERAMALAN IPR UNTUK TEKANAN RESERVOIR DI ATAS TEKANAN BUBBLE POINT PADA RESERVOIR BERTENAGA DORONG GAS TERLARUT Oleh: Dody Irawan Z Pembimbing: Dr. Pudjo Sukarno Amega Yasutra, MS Sari Selama reservoir memproduksi fluida reservoir maka tekanan reservoir akan mengalami penurunan, demikian pula halnya tekanan alir di dasar sumur juga akan cenderung menurun. Jika pada mulanya tekanan reservoir lebih besar dari tekanan bubble point, maka dengan berjalannya waktu produksi, tekanan dasar sumur akan mengalami penurunan dan dapat mencapai tekanan yang lebih rendah dari tekanan bubble point. Pada keadaan ini, gas yang terlarut dalam minyak akan terbebaskan sebagai fasa gas dan jika saturasi gas telah melampaui harga saturasi gas kritik, maka gas akan mengalir bersama minyak. Pada keadaan ini di dalam reservoir mengalir fasa gas dan minyak. Namun jika sebagian dari reservoir terutama yang jauh dari lubang sumur masih mempunyai tekanan diatas tekanan bubble point, maka pada bagian reservoir ini masih terjadi aliran minyak satu fasa. Keadaan aliran ini berbeda dengan model pengembangan Kurva IPR tak berdimensi yang dilakukan oleh Vogel, Sukarno, Wiggins, dan beberapa peneliti yang lain yang memodelkan sumur pada tekanan reservoir di bawah tekanan bubble point. Untuk menggambarkan kemampuan reservoir memproduksikan fluida pada kondisi tekanan reservoir diatas tekanan bubble point, maka dikembangkan kurva IPR (Inflow Performance Relationship) tak berdimensi yang baru. Ada beberapa persamaan IPR yang telah dipublikasikan dan diterima secara umum penggunaannya seperti Gilbert4, Vogel9, Sukarno10, dll. Begitu juga dengan persamaan peramalan kurva IPR yang telah dipublikasikan oleh Standing8, Eckmier5, dan Sukarno10. Persamaan-persamaan tersebut dikembangkan berdasarkan anggapan bahwa tekanan reservoir lebih rendah dari tekanan bubble point. Makalah ini bertujuan untuk mempelajari adanya pengaruh tekanan reservoir yang berbeda dan menghasilkan usulan persamaan peramalan IPR tak berdimensi yang sesuai dengan kondisi tersebut. Keadaan ini dimodelkan menggunakan simulator CMG dengan menggunakan data batuan dan fluida reservoir seperti yang digunakan Vogel9 dalam melakukan penelitiannya. Kata kunci : Produktvitas Sumur, kurva IPR dua fasa, peramalan kurva IPR
Abstract The reservoir pressure and the flowing well bore pressure will decline continuously during the production of reservoir fluids. In the case of reservoir pressure above the bubble point pressure, when the decreasing well bore pressure has passed below the bubble point pressure, then the dissolved gas in oil will be liberated as gas phase in the form of bubbles and flow with oil. This condition occurs when the saturation value of gas has passed its critical value, at this time there are gas phase and oil phase flow in the reservoir. When the pressure away from the well bore is still higher than the bubble point pressure, then in the reservoir there would be two section of flow, which are two phase flow in the section that close to the well bore, and single phase flow in the other section up to the reservoir boundary. This condition were modeled by Vogel, and other researchers, where most of them assuming that the initial reservoir pressure equal to or less than the bubble point pressure. Therefore, an IPR curve is needed to depict reservoir ability in producing fluid, in this condition. There are some IPR equations that have been published and accepted in general, like Gilbert4,Vogel9, Sukarno10, etc. There are also IPR equation for prediction purposes that has been published by Standing8, Eckmier5, and Sukarno10. These equations were developed with assumption of below bubble point pressure. This Paper investigates the influence of different bubble point pressure, the study also expected to give the new IPR curve forecasting equation referred to the given condition. The study was done by CMG simulator, the study also using the same rock data and reservoir fluid as those which used by Vogel9 in conducting its research. Keywords ,Well Productivity, Two phase IPR, future IPR prediction *
Mahasiswa Program Studi Teknik Perminyakan ITB.
Dody Irawan Z , 12203051, Sem 2 2007/2008 1
1.
Pendahuluan
Kurva Inflow Performance Relationship (IPR) merupakan plot antara laju alir fluida reservoir dengan tekanan dasar sumur, yang menunjukkan produktivitas suatu sumur. Untuk keperluan perencanaan produksi dan analisa ekonomi dimasa yang akan datang diperlukan persamaan peramalan kurva IPR. Sudah banyak persamaan peramalan IPR yang dikembangkan oleh para ahli, seperti peramalan Eckmeir5, Sukarno10, dll. Semua peramalan tersebut dikembangkan untuk tekanan reservoir di bawah tekanan bubble point dan pada harga tekanan bubble point konstan. Sedangkan peramalan kurva IPR untuk kondisi reservoir yang mempunyai tekanan reservoir diatas tekanan bubble point belum dilakukan penelitiannya. Namun Riri1 melakukan penilitian bahwa untuk solution gas drive, dengan berjalannya waktu produksi, terjadi perubahan komposisi minyak sebagai akibat terbebaskannya gas dari minyak, maka akan terjadi perubahan tekanan bubble point. Dari penelitian itu Riri1 mengusulkan persamaan kurva IPR dua fasa dengan memperhatikan tekanan bubble point yang berubah. Adanya perubahan tekanan bubble point disebabkan oleh perubahan komposisi minyak sebagai akibat terbebaskannya gas dari minyak. Akibat tekanan bubble point yang berubah tersebut, maka perlu dilakukan peramalan kurva IPR untuk tekanan reservoir diatas tekanan bubble point, namun tekanan alir dasar sumur dibawah tekanan bubble point.
2.
Studi Pustaka
Solution gas drive merupakan salah satu mekanisme pendorong alami untuk mendorong minyak dari reservoir ke sumur produksi. Menurut Cole(1969)4, reservoir dengan energi pendorong solution gas drive memiliki ciri-ciri sebagai berikut: 1. Volume konstan, tidak ada perubahan pada ukuran awal reservoir, dan tidak ada aliran air dalam reservoir. 2. Penurunan tekanan yang terjadi sangat cepat dan kontinyu akibat tidak terdapatnya fasa air dan gas bebas yang mengisi kekosongan dari reservoir setelah diproduksikan. 3. Pada awal produksi Gas Oil Ratio (GOR) di permukaan rendah kemudian mengalami peningkatan yang cepat dan turun kembali. 4. Gas keluar sebagai larutan bersama minyak tetapi tidak membentuk gas cap. 5. Faktor perolehan 5 % - 30 %. Pada profil produksi solution gas drive Gambar 1, terlihat bahwa saat tekanan dasar sumur di atas
tekanan bubble point, fasa yang mengalir adalah fasa minyak. Selama operasi produksi, tekanan dasa, gas yang terlarut dalam minyak mulai terbebaskan membentuk gelembung-gelembung gas bebas. Saat saturasi gas di dalam batuan berpori itu melampaui harga saturasi gas kritik, maka akan terjadi aliran dua fasa.r sumur terus menurun dan pada saat tekanan dasar sumur turun sampai dibawah tekanan bubble point
Gambar 2. Profil Solution Gas Drive2 Untuk menggambarkan kemampuan produksi suatu sumur pada reservoir solution gas drive, Vogel9 mengembangkan kelakuan produksi sumur tersebut dengan asumsi-asumsi sebagai berikut: 1. Reservoir berbentuk lingkaran, terbatas, sumur terletak di tengah reservoir, dan seluruh tebal reservoir dalam kondisi terbuka. 2. Reservoir bersifat homogen dan isotropik dengan saturasi air konstan pada suatu harga saturasi connate water tertentu. 3. Pengaruh gravitasi di abaikan. 4. Kompreisbilitas air dan batuan diabaikan. 5. Pada kondisi awal tekanan reservoir berada pada bubble point pressure . 6. Tidak ada pengaruh skin. Berdasarkan hasil penelitiannya, Vogel9 menghasilkan persamaan kurva IPR tak berdimensi sebagai berikut:
Pwf Qo = 1 − 0.2 Qo max Pr
P − 0.8 wf Pr
2
..............(1)
Sukarno10 kemudian melanjutkan studi Vogel dengan menggunakan asumsi yang sama tetapi kali ini
Dody Irawan Z , 12203051, Sem 2 2007/2008 2
mengikutsertakan pengaruh skin dalam pengembangan kurva IPR-nya. Persamaan IPR tak berdimensi yang diusulkan oleh Sukarno10 kemudian dimodifikasi oleh Fadli dkk11 sehingga dihasilkan persamaan sebagai berikut:
Pwf Qo = 1 − 0.0446 Q o max Pr
P − 1 .0446 wf Pr
2
....(2)
Sedangkan kurva IPR untuk tekanan reservoir diatas tekanan bubble point (satu fasa) merupakan garis lurus yang dinyatakan sebagi berikut:
Pwf Qo = 1 − Q o max Pr
.......................(3)
Menurut Riri1, pada reservoir solution gas drive dimana terjadi perubahan tekanan bubble point yang didasari oleh perubahan solution gas-oil ratio (Rso) yang terjadi selama waktu produksi. Untuk itu Riri1 mengembangkan IPR dua fasa dengan memasukkan unsur tekanan bubble point yang selalu berubah. Persamaan kurva IPR tak berdimensi yang diusulkan Riri1 tanpa memasukkan pengaruh skin adalah sebagai berikut:
Pwf Qo = 1 − 0.27 Qo max Pr
P − 0 .73 wf Pr
2
...........(4)
Untuk meramalkan produktivitas sumur di masa yang akan datang, maka diperlukan pengembangan persamaan peramalan kurva IPR. Prosedur untuk mendapatkan persamaan peramalan kurva IPR yang diperkenalkan pertama kali oleh Standing8 dengan memodifikasi persamaan Vogel9, sehingga didapatkan persamaan peramalan sebagai berikut:
Q o max f Qo max p
=
m (Prf
) …….….........(5) m (P ) rp
Dimana m(Pr) adalah pseudo pressure function pada harga tekanan reservoir tertentu. Selain Standing8 terdapat beberapa persamaan peramalan kurva IPR lain, yaitu persamaan peramalan IPR Eckmier5 dan persamaan peramalan kurva IPR Sukarn10. Persamaan peramalan Eckmier5 merupakan gabungan antara persamaan Vogel9 dan Fetkovich 8, dengan anggapan bahwa harga n (flow exponent) sama dengan satu, yang diperoleh persamaan peramalan sebagai berikut:
Qo max f Qo max p
Prf = Prp
3
.....................(6)
Sukarno10 mengembangkan persamaan peramalan kurva IPR untuk sumur yang berproduksi dari reservoir bertenaga dorong gas terlarut, dengan anggapan bahwa tidak terjadi perubahan jari-jari pengurasan dan faktor skin, dalam bentuk perbandingan antara laju produksi maksimum pada dua waktu produksi yang berbeda yang dinyatakan dalam persamaan berikut: Untuk API < 40, Q o max f Prf ……(5) = 0.033210 exp 3.429922 Qo max p P rp Untuk API > 40, Q o max f Prf ……(6) = 0.015215 exp 4.152343 Qo max p Prp Karena persamaan peramalan IPR untuk tekanan reservoir di atas tekaan bubble point belum tersedia, maka persamaan usulan yang diperoleh dibandingkan dengan ketiga metode peramalan diatas, yaitu metode Standing8, Eckmier5 dan Sukarno10 sebagai pembanding.
3.
Pengembangan Model
Reservoir dengan solution gas drive dimodelkan menggunakan simulator CMG berbentuk silinder dengan pembagian grid block menggunakan fungsi logaritmik, yaitu grid block makin membesar semakin jauh dari lubang sumur. Ketebalan seragam, homogen dan isotropik. Pengaruh tekanan kapiler dan gravitasi di reservoir diabaikan. Di dalam reservoir pada keadaan awal terdapat minyak dan connate water dengan tekanan reservoir di atas tekanan bubble point 3.1 Model Reservoir Reservoir radial dengan sumur vertikal terletak ditengah reservoir menembus seluruh lapisan produktif. Penelitian ini menggunakan model sumur tunggal. Dengan model ini dapat diketahui secara jelas mengenai gerakan aliran minyak dan gas menuju sumur produksi. Model reservoir yang digunakan pada penelitian ini dapat dilihat pada Gambar 2. 3.2 Data Reservoir Data reservoir yang digunakan dalam makalah ini adalah sama dengan data reservoir dari Vogel9. Kondisi reservoir pada keadaan undersarurated (Pr ≥ Pb) dan bertenaga dorong gas terlarut dimana reservoir bersifat homogen, isotropik dan sistem reservoir terbatas. Data
Dody Irawan Z , 12203051, Sem 2 2007/2008 3
model dan batuan reservoir dapat dilihat pada Tabel 1, di Lampiran A.
akan terdapat hubungan antara plot PD vs t D pada kertas grafik semi-log yang menghasilkan kemiringan garis 1.151. Untuk aliran minyak(satu fasa),
PD =
h (Pr − Pwf ) k o 141 .2Qo µ o Bo
................... (8)
tD =
0 .0002637 2 φrw
ko µ oCt
.................... (9)
Dari data tekanan dan waktu hasil perhitungan simulator dibuat plot PD vs tD yang hasilnya ditunjukkan pada Gambar 3. Pd vs t d
Gambar 2. Model Reservoir
16 14
Kurva permeabilitas relatif dari batuan reservoir dan sifat fisik dari fluida minyak dan gas dapat dilihat pada Lampiran A. Data sifat fisik dari fluda reservoir yang digunakan pada penelitian ini dapat dilihat pada Tabel 2.
12
Pd
10 8 6 4
4.
y = 1. 151Ln(x) - 1. 415 R2 = 1
2
Initialisasi dan Validasi Model
0
Untuk memastikan bahwa model akan memberikan hasil yang benar, perlu dilakukan initialisasi dan validasi. Setelah itu barulah dilakukan perhitungan parameter yang dibutuhkan dengan menggunakan model (simulator) tersebut.
10000
100000
1000000
td
Gambar 3. Chart pD vs tD
4.1 Initialisasi Model
Dari plot tersebut dihasilkan harga slope sebesar 1.15 yang harganya cukup mendekati harga 1.151 sehingga model dapat dikatakan valid untuk penelitian ini.
Initialisasi model dilakukan dengan membandingkan OOIP hasil similator terhadap volumetrik. Dari hasil perhitungan dengan menggunakan simulator diperoleh hasil OOIP sebesar 7.914 MMRBBL. Perhitungan cadangan dengan metoda volumetrik dengan persamaan,
5.
OOIP =
Ah φ (1 − S wc ) ........................(7) Boi
Dengan parameter-parameter pada Tabel 1 dan 2, OOIP dengan persamaan di atas bernilai 7.829 MMRBBL. Perbedaan OOIP dari simulator dengan perhitungan volumetrik cukup kecil yaitu sekitar 1.07 %. 4.2 Validasi Slope 1.151 Pada awal produksi dimana aliran masih satu fasa (tekanan di atas tekanan bubble point) untuk model silinder homogen dan isotropis, pada periode transien
Metode Penelitian
Penelitian tentang pengaruh adanya perubahan tekanan bubble point saat tekanan reservoir berada di atas tekanan bubble point disimulasikan dengan memasukkan harga tekanan bubble point awal dalam simulator dan sumur diproduksikan dengan laju produksi konstan. Saat tekanan dasar sumur turun di bawah tekanan bubble point awal dimana respon tekanan belum sampai ke batas reservoir sumur ditutup. Ketika sumur diproduksikan kembali, tekanan bubble point di sekitar sumur yaitu sampai batas respon tekanan lebih tinggi dari tekanan buble point awal. Kondisi tekanan reservoir saat sumur ditutup dan dibuka kembali dapat dilihat pada Gambar 4 dan Gambar 5 dengan tekanan bubble point awal 1000 psi.
Teknik Perminyakan-FTTM-ITB, Sem 2 2007/2008 4
Peramalan kurva IPR untuk tekanan reservoir di atas tekanan bubble point memerlukan tekanan reservoir dan laju alir maksimum minyak yang diperoleh dari kurva IPR dimana tekanan reservoir masih diatas tekanan bubble point awal. radial oil Pressure (psi) 1901-10-30 -10
K layer: 1
0
File: 300.irf User: dody irawan z Date: 4/1/2006
10
Scale: 1:45.429162 Y/X: 1.00:1 Axis Units: ft
1,153 1,127
0
0
NewWell
1,102 1,076 1,050 1,024 999
-10
-10
973 947 0.00
2.50 0.00
5.00 feet
0.50
921
1.00 meters
896 -10
0
10
Gambar 4. Distribusi Tekanan Saat Sumur Ditutup radial oil Bubble Point Pressure (psi) 1901-11-30 -10
0
K layer: 1 File: 300.irf User: dody irawan z Date: 4/1/2006
10
Scale: 1:45.429162 Y/X: 1.00:1 Axis Units: ft
0
0
1,124
NewWell
mengalami penurunan. Untuk reservoir yang tidak berbentuk silinder maka periode late transient akan ditemui, pada periode ini belum semua respons tekanan mencapai batas reservoir. Ketika penurunan tekanan telah mencapai batas reservoir, maka tekanan pada batas reservoir akan berkelakuan seperti penurunan tekanan pada volume konstan ini berarti tekanan reservoir akan mengalami penurunan dengan laju penurunan yang tetap. Kondisi inilah yang di sebut kondisi pseudo steady-state, oleh sebab itu tekanan pada kondisi ini dapat dianggap mewakili tekanan reservoir.
1,101
Dalam melakukan peramalan diperlukan beberapa harga harga tekanan reservoir dan laju alir maksimum minyak pada masa yang akan datang. Untuk mendapatkan data tersebut, mula-mula sumur diproduksikan dengan laju produksi konstan kemudian sumur ditutup sampai tekanan di reservoir konstan. Tekanan reservoir yang konstan tersebut merupakan tekanan reservoir yang akan datang. Selanjutnya sumur diproduksi kembali, untuk mendapatkan IPR selanjutnya. Hal ini dilakukan berulang-ulang sampai tekanan reservoir turun dibawah tekanan bubble point. Metode untuk mendapatkan IPR terlampir pada Lampiran B. Sedangkan laju alir maksimum minyak diperoleh dengan persamaan Riri1. Persamaan IPR usulan diperoleh dengan membuat plot antara Qomaxf/Qomaxp vs Prf/Prp. Selanjutnya Persamaan peramalan IPR usulan ini dibandingkan dengan persamaan peramalan Standing8, Eckmier5, dan Sukarno10.
1,078 1,055 1,032 1,010 987 964 -10
-10
0.00
2.50 0.00
0.50
5.00 feet
941 919
1.00 meters
896 -10
0
10
Gambar 5. Distribusi Tekanan Saat Sumur Dibuka Kurva IPR tersebut diperoleh dengan memplot tekanan dasar sumur terhadap laju alir minyak dengan berbagai sensitivitas sifat fisik batuan dan fluida. Penentuan tekanan alir dasar sumur dilakukan dengan memproduksi sumur dengan laju alir konstan hingga mencapai kondisi pseudo-steady state. Ketika reservoir diproduksi dengan laju alir konstan, maka di dalam reservoir tekanan akan turun melalui tiga periode yaitu transient, late transient dan pseudo-steady state. Pada periode transient penurunan tekanan berkelakuan seperti penurunan tekanan di reservoir yang tak berbatas dikarenakan tekanan di batas luar reservoir belum
6. Studi Sensitivitas Studi sensitivitas perlu dilakukan untuk mengetahui persamaan peramalan kurva IPR usulan berlaku pada rentang sifat fisik batuan dan fluida. Penentuan parameter sensitivitas untuk penelitian ini didasarkan pada parameter-parameter yang dianggap dapat mempengaruhi pelepasan gas dari minyak. Sifat fisik batuan dan fluida yang diuji sebagai berikut: 6.1 Sensitivitas API Gravity Sensitivitas API Gravity mewakili perubahan kelakuan produksi untuk sifat fisik fluida yang berbeda. Semakin besar API Gravity maka membuat fraksi ringan dalam minyak semakin banyak sehingga viskositas ikut turun dan laju alir minyak meningkat. Gambar 6 menunjukkan pengaruh dari perubahan API terhadap kurva IPR. Dari gambar tersebut terlihat bahwa perubahan harga API minyak sangat kecil pengaruhnya terhadap kurva IPR untuk rentang API Gravity 30-60.
Dody Irawan Z , 12203051, Sem 2 2007/2008 5
mudahnya gas keluar dari minyak dengan tingginya temperatur.
6.2 Sensitivitas Specific Gravity Gas Sensitivitas terhadap Specific Gravity gas menunjukkan kandungan fasa ringan dalam minyak. Semakin besar Specific Gravity gas maka kandungan fasa ringan yang terkandung dalam minyak semakin banyak sehingga minyak semakin mudah mengalir. Gambar 7 menunjukkan pengaruh dari perubahan Specific Gravity terhadap kurva IPR. Dari gambar tersebut terlihat bahwa dengan meningkatnya Specific Gravity dari 0.7 sampai 0.9 akan meningkatkatkan sedikit laju alir minyak.
6.4 Sensitivitas Saturasi Gas Kritik Dari hasil run simulator, semakin besar saturasi gas kritik semakin lama waktu yang diperlukan gas untuk ikut mengalir bersama minyak. Pengaruh sensitivitas saturasi gas kritik terhadap kurva IPR dapat dilihat pada Gambar 9. IPR
2500
IPR
3500
2000 Pwf (psi)
3000 Pwf (psi)
2500 2000 1500
1500 1000 500
1000
0
500
0
0 0
2000 A PI 30
4000 6000 Qo (bbl/day) A PI 35
A PI 45
8000
1000
2000
10000
3000 4000 Qo (bbl/d)
T 150
A PI 60
T 200
5000
6000
T 250
Gambar 8. Kurva IPR Sensitivitas Temperatur (Pb=1000 psi, Sgc=0.1, SG=0.9, API=35)
Gambar 6. Kurva IPR Sensitivitas API Grafity (Pb=1500 psi, Sgc=0.1, T=200 F, SG=0.9)
IPR
3000
IPR
2500
2500 Pwf (psi)
Pwf (psi)
2000 1500
2000 1500 1000
1000
500
500
0
0 0
1000
2000
SG 0.7
3000 4000 Qo (bbl/d) SG 0.8
5000
6000
SG 0. 9
Gambar 7. Kurva IPR Sensitivitas Specific Grafity Gas (Pb=1000 psi, Sgc=0.1,T=200 F, API=35) 6.3 Sensitivitas Temperatur Pengaruh sensitivitas temperatur terhadap kurva IPR dapat dilihat pada Gambar 8. Pada gambar tersebut menunjukkan laju alir minyak meningkat dengan naiknya temperatur. Hal ini disebabkan oleh semakin
0
1000 Sgc 0. 05
2000
3000 4000 Qo (bbl/d) Sgc 0. 1
5000
6000
Sgc 0. 15
Gambar 9. Kurva IPR Sensitivitas Saturasi Gas Kritik (Pb=1000 psi, T=200 F, SG=0.9, API=35) 6.5 Sensitivitas Tekanan Bubble Point Sensitivitas terhadap tekanan bubble point menunjukkan bahwa laju alir minyak akan lebih besar bila tekanan buble ponit jauh lebih kecil dari tekanan reservoir. Hal ini disebabkan oleh ikut terproduksikannya gas bersama dengan minyak. Pengaruh sensitivitas tekanan bubble point terhadap IPR dapat dilihat pada Gambar 10. Teknik Perminyakan-FTTM-ITB, Sem 2 2007/2008 6
pada model adalah black oil. Untuk itu data sensitivitas API 60 tidak digunakan dalam persamaan peramalan kurva IPR usulan pada Gambar 12.
IPR
3500 3000 2000
1
1500
0.9
1000
0.8 Qo max2/Qo max1
Pwf (psi)
2500
500 0 0
2000
4000 6000 Qo (bbl/d)
Pb 2000
8000
10000
0.7 0.6 0.5 0.4
Pb 1500
Pb 1000
y = x1.04 R2 = 0.99
0.3
Gambar 10. Kurva IPR Sensitivitas Tekanan Bubble Point (T=200, psi, Sgc=0.1, SG=0.9, API=35)
0.2 0.2
0.4
0.6 Pr2/Pr1
0.8
1
Gambar 12. Peramalan IPR Usulan 7.
Dengan melakukan power regression dari data model maka didapatkan persamaan:
Peramalan Kurva IPR
Qomf/Qomp
Data hasil uji sensitivitas kurva IPR digunakan untuk membuat persamaan peramalan IPR. Hasil uji sensitivitas kurva IPR terhadap faktor perolehan minyak pada tiap kasus diperoleh pasangan data tekanan reservoir (Pr) dan laju alir maksimum (Qmaks). Data tersebut terlampir pada Lampiran C. Tiap pasang data ini kemudian dibandingkan yang satu dengan yang lain. Hasil plot dapat dilihat pada Gambar 11. 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.2 BA SE Sgc=0. 05
0. 4 SG 0. 7 Sgc 0.15
Gambar 11. Plot
0. 6 Prf/Prp SG 0.8 A PI 30
T 150 A PI 45
Pr 2
Q o maks 2
Pr 1
Vs
0.8
Q o max f Qo max p
Prf = Prp
1.04
................... (10)
Dengan adanya ekponen sebesar 1.04 pada persamaan peramalan usulan, menunjukkan bahwa dengan adanya tekanan bubble point yang berbeda pada solution gas drive membuat laju alir pada tekanan reservoir diatas tekanan bubble point tidak sepenuhnya satu fasa. Hal ini disebabkan oleh terbebasnya gas hingga melewati saturasi gas kritis disekitar lubang sumur. Peristiwa ini dapat diamati pada saturasi fluida di sekitar sumur Gambar 13.
1 T 250 A PI 60
Q o maks 1
Dari Gambar 11 terlihat bahwa penyebaran data tidak baik untuk harga API 60. Berdasarkan teori, API Gravity 60 merupakan range dari sifat fluida resorvoir jenis retrogade gas. Sedangkan fluida yang digunakan
Gambar 13. Saturasi Fluida Di Sekitar Sumur (Pb=1000 psi, Sgc=0,1, T=100 F, API=45)
Dody Irawan Z , 12203051, Sem 2 2007/2008 7
Peramalan kurva IPR usulan ini dibandingkan dengan peramalan metode Standing8, Eckmeir5, dan Sukarno10. Grafik perbandingan dapat dilihat pada Gambar 14. 1
Qomf/Qomp
0.8 0.6 0.4 0.2
Peramalan IPR untuk tekanan reservoir di bawah tekanan bubble point mengunakan persamaan Eckmeir5 Dapat dilihat pada Gambar 16 sedangkan persamaan Sukarno10 pada Gambar 17. Dari gambar tersebut dapat dilihat bahwa peramalan IPR untuk tekanan reservoir di atas tekanan bubble point yang menggunakan persamaan peramalan IPR usulan tidak terdapat perbedaan yang signifikan dengan IPR hasil simulator. Sedangkan untuk peramalan IPR untuk tekanan reservoir di bawah tekanan bubble point terdapat sedikit perbedaan. Hal ini disebabkan persamaan IPR yang digunakan adalah persamaan IPR Riri1. Data IPR simulator, peramalan IPR Eckmeir5, dan peramalan IPR Sukarno10 Dapat dilihat pada Lampiran D.
0 0.2
0.4
0.6
0.8
1
EIKHMEIR
PS -A PI<40
PS -A PI>40
3500 3000
Gambar 14. Perbandingan Peramalan IPR Dari gambar tersebut terlihat bahwa terdapat perbedaan yang nyata antara peramalan kurva IPR di atas tekanan bubble point dengan peramalan dibawah tekanan bubble point. Persamaan peramalan model hampir berimpit dengan persamaan peramalan Standing8. Hal ini disebabkan oleh pseudo pressure function yang digunakan pada persamaan Standing8 merupakan fungsi dari tekanan reservoir diatas tekanan bubble point. Dengan menggunakan persamaan yang didapatkan dari model, kita akan mendapatkan harga Qo maks future yang akan menjadi input kembali ke persamaan Riri1 untuk mendapatkan kurva IPR saat tekanan reservoir diatas tekanan bubble point pada solution gas drive. Untuk membuktikan hal tersebut, maka dilakukan percobaan pada tekanan bubble point 2000 psi. Gambar 15 memperlihatkan hasil percobaan dengan simulator
Pwf(psi)
1500 1000 500 0 0 IPR#1
2000 IPR#2
IPR#3
4000
6000 8000 10000 Qo (bbl/d) IPR#4 IPR#5 IPR#6 IPR#7
12000 IPR#8
Gambar 16. Peramalan IPR Eikhmeir IPR
4000 3500 3000 2500
500 0
3000 2500 2000 1500
0 IPR#1
1000 500 0
2000 IPR#2
IPR#3
4000
6000 8000 10000 12000 Qo (bbl/d) IPR#4 IPR#5 IPR#6 IPR#7 IPR#8
Gambar 16. Peramalan IPR Sukarno 0
IPR#1
2000
2000 1500 1000
IPR
4000 3500
2500
Pwf (psi)
STA NDING
Pwf (psi)
MODEL
IPR
4000
Prf/Prp
2000 IPR#2
IPR#3
4000
6000 Qo(bbl/d)
IPR#4
IPR#5
8000
10000
12000
8. IPR#6
Gambar 15. IPR Hasil Simulator
IPR#7
Kesimpulan
IPR#8
Berdasarkan penelitian yang telah dilakukan maka dapat diambil kesimpulan sebagai berikut: Teknik Perminyakan-FTTM-ITB, Sem 2 2007/2008 8
1.
2.
Dengan berubahnya tekanan bubble point disekitar sumur akan membuat saturasi gas didaerah sekitar sumur berubah, sehingga bila saturasi gas bertambah dan melewati saturasi gas kritik aliran yang terjadi disektar sumur menjadi 2 fasa. Makalah ini mengusulkan persamaan peramalan IPR di masa depan pada tekanan reservoir berada diatas tekanan bubble point yang didapatkan secara empiris menggunakan data simulasi. Persamaan yang diusulkan adalah:
Q o max f Qo max i 3.
9.
Prf = Pri
2.
= waktu dimensionless = waktu, jam = waktu periode pseudosteady- state, jam = porositas, fraksi = viskositas minyak, cp = tekanan reservoir sebagai base, psi = tekanan reservoir pada waktu tertentu, psi = tekanan dasar sumur, psi = tekanan reservoir saat ini, psi = tekanan reservoir yang akan datang, psi = Laju alir maksimum di masa depan (bbl/d) = Laju alir maksimum saat ini (bbl/d).
1 .04
Selama Peramalan IPR masa depan dibawah tekanan bubble point untuk reservoir solution gas drive dengan tekanan bubble point berubah belum ada, maka dapat menggunakan metode yang telah ada seperti Eckmeir5, dan Sukarno10.
Saran 1.
tD t tc φ µo Pr 1 Pr 2 Pwf prf prp Qomaksf Qomaksp
Perlu dilakukan peramalan IPR di atas tekanan bubble point yang ikut memperhatikan faktor skin pada solution gas drive reservoir. Perlu dilakukan peramalan IPR untuk tekanan reservoir di bawah tekanan bubble point Untuk mendapatkan peramalan IPR yang tepat.
10. Ucapan Terima Kasih Syukur Alhamdulillah penulis ucapkan atas kehadirat Allah SWT dan segala rahmat serta karunia-Nya sehingga penulis dapat menyelesaikan tugas akhir ini. Kepada keluarga, penulis mengucapkan terima kasih kepada orang tua dan adik-adik yang selalu memberikan dukungan, doa, dan kasih sayangnya. Dengan segenap ketulusan hati penulis mengucapkan terima kasih kepada Dr. Pudjo Sukarno dan Amega Yasutra MS, sebagai pembimbing tugas akhir ini yang telah meluangkan waktu, tenaga, dan pikiran dalam membantu penulis menyelesaikan permasalahan-permasalahan yang di hadapi dalam penyelesaian tugas akhir ini. Terakhir penulis ingin mengucapkan terima kasih kepada perusahaan software CMG yang telah mengizinkan penulis untuk menggunakan software CMG dalam studi penelitian ini.
11. Daftar Pustaka 10. Daftar Simbol 1. Bo Bg cf GOR h k kr m(Prf) m(Prp) p pb pD pr Rso rw re So Sgc Sw
= faktor volume formasi minyak, rb/stb = faktor volume formasi gas, rb/Mscf = kompresibilitas batuan, 1/psi = Gas Oil Ratio, scf/stb = ketebalan vertikal, ft = permeabilitas, mD = permeabilitas relatif, fraksi = Pseudo pressure function masa depan, psi = Pseudo pressure function sekarang, psi = tekanan, psi = tekanan bubble point, psi = tekanan dimensionless = tekanan reservoir, psi = solution gas-oil ratio, scf/stb = jari-jari sumur, ft = jari-jari pengurasan, ft = saturasi minyak, fraksi = saturasi gas kritik, fraksi = saturasi air, fraksi
2.
3.
4. 5.
6.
Riri, Armelia Suska, “Pengembangan IPR 2 Fasa pada Reservoir dengan Solution Gas Drive Pada Kondisi Tekanan Reservoir Di Atas Tekanan Bubble Point”,Tesis, TM-ITB, 2006. Smith, C.R., G.W., dan Farrar, R.L., “Applied Reservoir Engineering Vol.2”, OGCI Publications, Tulsa, Oklahoma, 1992. Tarek,Ahmed., “Reservoir Engineering Handbook”, 2nd edition, Gulf Publishing Company, Houston, Texas,2000. Tarek,Ahmed., “Advance Reservoir Engineering”, Gulf Publishing Company, Houston, Texas,2000. Brown, K.E.: Technology of Artificial Lift Methods, Volume IV, Penn Well Publishing Company, Tulsa, Oklahoma, 1980. McCain, W.D.: The Properties of Petroleum Fluids, second Edition, Penn Well Publishing Company, Tulsa, Oklahoma, 1990.
Dody Irawan Z , 12203051, Sem 2 2007/2008 9
7.
Earlougher, R.C., “Advances in Well Test Analysis”, 2nd printing, Society of Petroleum Engineering of AIME, New York, Dallas, 1977. 8. Beggs, Dale.H., “Production Optimization”, OGCI Publications, Tulsa, Oklahoma, 1991. 9. Vogel, J.V.: “Inflow Performance Relationship Curve for Solution Gas Drive Wells,” JPT, Januari 1968. 10. Sukarno, Pudjo :”Inflow Performance Relationship In Two and Three Phase Conditions”, PhD
Dissertation, The University of Tulsa, Oklahoma, 1986. 11. Rahman, F., Yahya, M., dan Leonardus, S.: ”Modifikasi Persamaan Sukarno untuk Menentukan IPR Sumur Pada Reservoir dengan Tekanan di bawah Tekanan Saturasi,” Tugas Mata Kuliah Teknik Produksi, Jurusan Teknik Perminyakan, Institut Teknologi Bandung, 2006.
Teknik Perminyakan-FTTM-ITB, Sem 2 2007/2008 10
Lampiran A Tabel 1. Data Model dan Batuan Reservoir
Tabel 2. Data sifat fisik fluida reservoir
Dody Irawan Z , 12203051, Sem 2 2007/2008 11
Data PVT dan Permeabilitas relatif Base Case
Teknik Perminyakan-FTTM-ITB, Sem 2 2007/2008 12
Lampiran B Metode Penelitian Pb = 2000 Psi radial oil Well-1 10 Layer (Pb=2000).irf 3,500
1,200
1,000
800
2,500
600
400
Oil Rate SC - Daily (bbl/day)
Well Bottom-hole Pressure (psi)
3,000
2,000 200
1,500
0 1901-4
1901-7
1901-10 Time (Date)
1902-1
1902-4
Well Bottom-hole Pressure Oil Rate SC - Daily
IPR 4000
Pwf(psi)
3500 3000 2500 2000 1500 1000 500 0 0
2000 IPR#1
4000
6000 Q(bbl/d)
IPR#2
8000 IPR#3
10000
12000
IPR#4
Dody Irawan Z , 12203051, Sem 2 2007/2008 13
Pb = 1500 Psi radial oil Well-1 10 layer (Pb=1500).irf
Lampiran B Kurva Relatif Permeability untuk Base Case
3,500
1,200
1,000
600
2,000
Oil Rate SC (bbl/day)
800 2,500
400
1,500 200
1,000
0 1901-4
1901-7
1901-10 Time (Date)
1902-1
1902-4
Well Bot tom -hole Pressure Oil R ate SC
IPR 4000 3500 3000 Pwf(psi)
Well Bottom-hole Pressure (psi)
3,000
2500 2000 1500 1000 500 0 0
2000 IPR#1
4000 6000 Q(bbl/d) IPR#2 IPR#3 IPR#4
8000 IPR#5
10000 IPR#6
Teknik Perminyakan-FTTM-ITB, Sem 2 2007/2008 14
Pb = 1000 Psi
IPR
4000 3500
Pwf(psi)
3000 2500 2000 1500 1000 500 0 0 IPR#1
2000 IPR#2
IPR#3
4000 6000 Q(bbl/d) IPR#4
IPR#5
8000 IPR#6
IPR#7
10000 IPR#8
Dody Irawan Z , 12203051, Sem 2 2007/2008 15
Lampiran C Data Tekanan Reservoir dan Laju Alir Minyak Maksimum Pada Sensitivitas Pb, API, SG, T, dan Sgc
Teknik Perminyakan-FTTM-ITB, Sem 2 2007/2008 16
Dody Irawan Z , 12203051, Sem 2 2007/2008 17
Lampiran D Data IPR simulator, peramalan IPR Eikhmeir9, dan peramalan IPR Sukarno10
IPR Hasil Simulator
Peramalan IPR Eikhmeir
Peramalan IPR Sukarno
Teknik Perminyakan-FTTM-ITB, Sem 2 2007/2008 18