Institut Teknologi Bandung, Jl. Ganesha No. 10, Bandung Abstrak. Abstract

Prosiding Seminar Nasional Fisika (E-Journal) SNF2016 http://snf-unj.ac.id/kumpulan-prosiding/snf2016/

VOLUME V, OKTOBER 2016

p-ISSN: 2339-0654 e-ISSN: 2476-9398

DOI: doi.org/10.21009/0305020124

SOLUSI PERSAMAAN DIFUSIVITAS ALIRAN FLUIDA MINYAK PADA RESERVOIR DAN MODIFIKASINYA UNTUK UJI ALIR (PRESSURE DRAWDOWN TEST) DAN UJI TUTUP (PRESSURE BUILDUP TEST) Hardiyanto1,a), Lilik Hendrajaya 2,b) 1,2

Institut Teknologi Bandung, Jl. Ganesha No. 10 , Bandung 40132

Email: a)[email protected], b)[email protected] Abstrak Tekanan transien fluida minyak pada reservoir dapat diperoleh dari solusi persamaan difusivitas. Untuk mencari solusi persamaan difusivitas ini diperlukan dua syarat batas sebagai syarat dalam dimensi ruang dan satu syarat keadaan awal sebagai syarat dimensi waktu. Solusi tersebut ditentukan dengan menggunakan Transformasi Laplace. Solusi tersebut selanjutnya digunakan sebagai analisis uji penurunan tekanan pada Uji alir (Pressure Drawdown Test) dan Uji Tutup (Pressure Buildup Test). Hasil plot dari solusi diperoleh nilai parameter reservoir diantaranya permeabilitas, faktor “skin” (kerusakan sekitar sumur), volume reservoir, Dietz Shape faktor dan tekanan rata-rata. Kata-kata kunci: persamaan difusivitas, tekanan transien,drawdown-Buildup.

Abstract Pressure transient of oil in reservoir can be obtained from the solution diffusivity equation. To find a solution of this diffusivity equation takes two boundary conditions as space dimension and one dimension of the initial state as a condition for the time dimension. The solution is determined by using the Laplace transform and then used as a test analysis of pressure drop in a flow test (Pressure Drawdown Test) and Pressure buildup Test. The plot of the solution used to prediction reservoir parameters such as permeability, skin faktor , volume of the reservoir, Dietz Shape faktor and the average pressure. Keywords: diffusion equation, the pressure transient, the drawdown-buildup.

1. Pendahuluan Suatu karakterisitik reservoir migas yang perlu dianalisa untuk mengetahui suatu kapasitas atau kemampuan produksi suatu reservoir diantaranya adalah besarnya porositas dan permeabilitas. Sehingga dikembangkan suatu model matematik yang mendeskripsikan aliran dalam reservoir sebagai media aliran berpori yang dikenal dengan persamaan difusivitas. Solusi persamaan difuivitas ini bisa digunakan sebagai dasar analisis aliran dalam reservoir disebut dengan test tekanan transient. Test tekanan transient yang sering diterapkan pada sumur produksi umumnya terdiri dari uji alir (Drawdown test) dan uji tutup (build up test). Kedua uji tersebut digunakan untuk mengetahui nilai permeabilitas maupun faktor-faktor lain yang ikut mempengaruhi prilaku aliran fluida dalam media berpori atau reservoir seperti faktor kerusakan sekitar lubang sumur atau biasa disebut dengan faktor skin (skin faktor), bentuk reservoir, volume pengurasan, dan tekanan rata rata reservoir. Dalam penelitian ini dilakukan penurunan kembali persamaan difusivitas menggunakan transformasi

laplace yang selanjutnya dimodfikasi sehingga bisa digunakan untuk memprediksi beberapa parameter reservoir

2. Metode Penelitian Penyelesaian secara analitik umumnya dilakukan menggunakan metode transformasi. Metode transformasi yang dapat digunakan adalah transformasi laplace. skema penggunaan transformasi dapat dilihat dari bagan berikut.

Gambar 1. Skema Transformsi Laplace Invers transformasi laplace dapat menggunakan cara analtik maupun numerik. Solusi analitik terhadap

Seminar Nasional Fisika 2016 Prodi Pendidikan Fisika dan Fisika, Fakultas MIPA, Universitas Negeri Jakarta

SNF2016-CIP-125


syarat awal untuk aliran satu fasa didalam media berpori bisa di peroleh dalam dua bentuk yaitu Solusi eksak, yaitu solusi dalam bentuk transformasi laplace Solusi pendekatan (approximation solution). Persamaan umum transformasi laplace dan inversnya dapat dituliskan sebagai berikut[4]. 

    f s    f t e st dt

L f t

(1)

0

dan inversnya dapat dituliskan: 1

L

f s   f t 

(2)

dengan f(t) suatu fungsi. Transformsi laplace maupun inversnya ini sudah dipermudah dari beberapa sumber buku [4] dalam bentuk tabel sehinggga bisa digunakan dengan mudah. Hanya ada beberapa invers transformasi memerlukan metode komputasi untuk menyelesaikannya. Invers trasnformasi laplace yang tidak bisa langsung dilihat dari tabel bisa digunakan metode penyederhanan yang terbagi menjadi beberapa methode diantarnya yaitu: a. Metode pecahan parsial. b. Metode Deret c. Metode Persamaan diferensial Dalam penelitian ini terdapat solusi persamaan difusi ini dalam bentuk modifikasi fungsi bessel sehingga diperlukan teknik atau beberapa metode untuk menyelesaikan invers transformasi laplace.

konstan dan tekanan sumur produksi konstan sebagai berikut: pD 

tD

kh( pi  p )

2,634 x104 kt , 

Persamaan diferensial dasar untuk aliran radial dalam media berpori ini akan menunjukkan bagaimana aliran fluida di daerah sekitar lubang bor seperti terlihat pada gambar 2. untuk mendapatkan persamaan diferesnisal diasumsikan aliran bergerak secara radial menuju lubang bor, reservoir dianggap homogen, reservoir memiliki permeabilitas isotropik, tekanan konstan dan formasi sepenuhnya disaturasi 1 fasa. Dengan mengaplikasikan hukum kekekalan aliran massa, persamaan keadaan dan hukum darcy didapat persamaan sebagai berikut : (3)

atau dapat dituliskan dalam satuan lapangan sbagai berikut: 2 d p 1 dp c dp   2  4 dt r dr dr 2,634 x10 k

(4)

Untuk mendapaatkan solusi analitik persamaan difusivitas, terlebih dahulu persamaan tersebut ditransformasikan kedalam bentuk tak berdimensi dengan mendefinisikan variabel tak berdimensi berdasarkan keadaan keadaan produksi di sumur produksi. Keadaaan sumur produksi dapat dibagi menjadi dua keadaan yaitu laju produksi sumur

rDe 

2

crw

r rw

(5)

re rw

Untuk mendapaatkan solusi analitik persamaan difusivitas, terlebih dahulu persamaan tersebut ditransformasikan kedalam bentuk tak berdimensi dengan mendefinisikan variabel tak berdimensi berdasarkan keadaan keadaan produksi di sumur produksi. Keadaaan sumur produksi dapat dibagi menjadi dua keadaan yaitu laju produksi sumur konstan dan tekanan sumur produksi konstan. Variabel tak berdimensi untuk tekanan sumur produksi konstan dapat dituliskan sebagai berikut: qD 

pD 

141, 2 qB



kh pi  p wf



(10)

pi  p pi  p wf

Sehingga persamaan difusi (persamaan 4) di transformasikan kedalam variabel tak berdimensi menjadi : 2 d p D

2



1 dpD r dr D

D

3. Hasil dan Pembahasan

rD 

,

141, 2 qB

dr

2 d p 1 dp c dp   2 r dr k dt dr

p-ISSN: 2339-0654 e-ISSN: 2476-9398




dp

(11)

D

dt

D

D

Persamaan diatas selanjutnya ditrasnformasikan kedalam ruang laplace untuk menyederhanakan persamaan difusi tersebut sehingga diperoleh persamaan berikut: 

dengan 

2  2 d p 1 dp 2   p  0 2  d d

 rD s

(12)

.

Persamaan 12 merupakan bentuk umum persamaan diferensial yang memiliki solusi umum dalam bentuk fungsi bessel termodisikasi. Solusi persamaan tersebut dapat dituliskan sebagai berikut:  (13) p  AI 0 rD s  BK0 rD s









Untuk mencari solusi khusus tertentu digunakan dua syarat batas diatas. Adapun dua syarat batas itu adalah syarat batas pada sumur dan syarat batas luar reservoir yang dalam penelitian ini penyelesaian solusi dibagi dalam beberapa kasus sebagai berikut: a. Kasus 1 Sumur diproduksi dengan laju produksi konstan dan reservoir tak terbatas (infinite acting). syarat batas dalam pada sumur : p 141, 2 qB lim r  rw0 r kh

Syarat batas luar reservoir:


SNF2016-CIP-126


lim p( r , t )  pi re e

Syarat batas tersebut di transformasikan ke ruang laplace sebelum digunakan untukmencari solusi khusus persamaan 13. Sehingga diperoleh solusi khusus sebagai berikut:



 1 p  K 0 rD s s



Persamaan tersebut di inverskan dari ruang laplace[4] sehingga diperoleh: pD

 r2    Ei   D   2   4t D  1




 

70,6 qB kh

  cr 2    0,00105kt   

(15)

Kasus 2 Sumur diproduksi dengan tekanan konstan dan jari-jari reservoir tak terbatas (infinite acting) Syarat batas dalam pada sumur :



pD 

 rD2  2 rDe  1   4 1

Syarat batas luar reservoir: lim p( r , t )  pi re e








  

K rD s  p rD , t D  0 sK 0 s





 (1  e  0

c.

2



u 2tD



) J 0 (u )Y0 (urD )  Y0 (u ) J 0 (urD ) 2 2 2 u J 0 (u )  Y0 (u )







(16)

Kasus 3 Sumur diproduksi dengan laju konstan dan reservoir tertutup (tidak ada aliran dibatas reservoir)

syarat batas dalam pada sumur : p 141, 2 qB lim r  rw0 r kh

atau jika ditulis dalam satuan lapangan dan mengganti variabel tak berdimensi menjadi:



p ( rw , t )  pi  141, 2

qB  2t D kh

 r 2  De



3

 ln r  De 4  



  (18)   2t   n D J 2  r   2e qB   1  De n     141, 2         kh n  1 2 2    J  r    J 2       n  1  De n  1  n     d.

Kasus 4 Sumur diproduksi dengan laju konstan dan tekanan konstan pada batas reservoir

p 141, 2 qB lim r  rw0 r kh

syarat batas luar reservoir:

   pi

p re , t

Syarat batas tersebut di transformasikan keruang laplace sebelum digunakan untukmencari solusi khusus persamaan 13. Sehingga diperoleh solusi khusus sebagai berikut:



 





     

 I rDe s K 0 rD s  K 0 rDe s I 0 rD s p 0 3 s 2 K 0 rDe s I1 s  I 0 rDe s K1 s

 

  

Menurut Matthews dan Russell ( Pressure Buildup and Flow Test in Wells, Society of Petroleum Engineers of AIME, 1967) yang mengikuti carslaw dan jaegar (Conduction of heat in solid, 1959)[3], solusi untuk kasus ini adalah :





p rw , t  pi  141, 2

syarat batas luar reservoir: p

(17)

syarat batas dalam pada sumur :

Invers untuk transformasi laplace ini menurut Van Everdingen dan Hurst (The Application of the Laplace Transformation to Flow Problems in Reservoirs, Petroleum Transactions AIME, 305-324, 1949) dapat dituliskan menjadi: p rD , t D 

2  rDe   tD   r 2  1 ln rD  De

4 4 2 2 3rDe  4 rDe ln rDe  2 rDe  rDe  1  2 2 4( rDe  1)

b.

p ( r  rw , t )  pw

           

Menurut Matthews dan Russell ( Pressure Buildup and Flow Test in Wells, Society of Petroleum Engineers of AIME, 1967)[3] sebagai berikut.

(14)

Ei 

 

 K1 rDe s I 0 rD s  I1 rDe s K 0 rD s p 3 s 2 I1 rDe s K1 s  K1 rDe s I1 s

atau jika ditulis dalam satuan lapangan dan mengganti variabel tak berdimensi menjadi: p r , t  pi 

p-ISSN: 2339-0654 e-ISSN: 2476-9398


qB kh



ln rDe



  n2 t D 2 e J 0 rDe  n  2  n 1  2 J 2   J 2 r  n 1 n 0 De n

qB

lim r  141, 2 0 r rDe r kh

  


SNF2016-CIP-127



(20)




e.

Kasus 5 Sumur diproduksi dengan tekanan konstan dan tekanan pada batas reservoir juga konstan

syarat batas dalam pada sumur :



p r  rw , t

  pw

syarat batas luar reservoir:

kondisi tidak ada aliran (no-flow outer boundary condition) dan sumur berproduksi dengan laju alir konstan. Sehingga digunakan persamaan 18 untuk menggambarkan aliran. Untuk priode aliran pseudosteady-state terjadi pada masa produksi yang sudah lama (pada harga t yang besar) maka solusi pendekatan dan asumsi rDe>>1 dan rD=1 dapat dapat dutuliskan:



lim p( r , t )  pi re e

                

 K rDe s I 0 rD s  I 0 rDe s K 0 rD s p 0 s K 0 rDe s I 0 s  K 0 s I 0 rDe s

Transformsi laplace dari solusi diatas tidak bisa diselesaikan secara manual. Namun Zimmerman(Flow in Porous media ,2003-2004) melakukan penurunan dengan cara yang berbeda mendapat kan solusi untuk kasus ini sebagai berikut :



ln rD / rDe



ln rDe n2tD

    

 U n (  n rD ) e J 0  n J 0 rDe  n   2 2 n1 J 0  n  J 0 rDe  n

 

dengan

  



  



(21)



Selanjutnya solusi solusi persamaan difusi dyang dikemukakan diatas dikoreksi dengan mengambil aproksimasi agar bisa digunakan untuk menngambarkan aliran saat proses uji alir dan uji tutup berlangsung. Aproksimasi ini berdasarkan pembagian waktu aliran pada saat gangguan tekanan diberikan. Pembagian dikelompokan ke dalam tiga periode yaitu periode transient, periode pseudosteady state dan steady state. Periode transient terjadi pada saat awal produksi ketika efek batas luar reservoir belum terasa di sumur sehingga reservoir berperilaku seperti halnya tidak ada batas (reservoir bersifat infinite-acting). Sehingga untuk keadaan ini digunakan persamaan 15 untuk menggambarkan keadaan aliran. Untuk harga argumen x yang kecil (yaitu x < 0.01) maka Ei(x) dapat didekati oleh harga logaritmik [1] sehingga persamaan 15 dapat dituliskan menjadi:

 

162,6 qB kh

2t D 3  ln rDe  2 4 rDe

atau jika ditulis dalam bentuk lain sebagai berikut: p wf  p i  141, 2

re qB  0,000527 kt 3  ln   (23)  kh  c re2 rw 4   

Steady-state terjadi pada selang waktu yang sangat besar atau lama ( sumur sudah diproduksi dengan sangat lama) pada suatu sistem reservoir dengan kondisi batas luar reservoir berupa tekanan konstan atau dengan kata lain tidak ada perubahan tekanan terhadap waktu. Sehingga untuk mendapatkan solusi yang sesuai dengan keadaan diatas maka dalam kasus diambil syarat batas constant pressure outary boundary dan constant-rate production. Solusi untuk periode ini dapat diambil dari persamaan yang solusinya merupakan persamaan darcy dalam arah radial. p wf  pi  141, 2

U n (  n rD )  Y0  n J 0 rD  n  J 0  n Y0 rD  n

p r , t  pi 



p D 1, t D 

Syarat batas tersebut di transformasikan keruang laplace sebelum digunakan untukmencari solusi khusus persamaan 13. Sehingga diperoleh solusi khusus sebagai berikut:

pD  

p-ISSN: 2339-0654 e-ISSN: 2476-9398


qB kh

ln

re

(24)

rw

Untuk membuktikan bahwa persamaan-persamaan diatas bisa digunakan untuk memprediksi suatu aliran pada suatu reservoir maka diambil data produksi dan reservoir sumur X[2]. q = 250 STBPD k = 246 md  = 0,8 cp h = 69 ft B =1,55 bbl/STB f = 0,039 Pi =4150 psi c = 0,000017 1/psi re =1489 ft rw = 0,198 ft Setelah data itu dimasukkan ke persamaan sesuai dengan periode periode yang telah ditentukan diatas maka diperoleh grafik sebagai berikut:

  kt     3, 23 (22) log  2   cr  

Pada periode pseudosteady-state terjadi ketika semua batas reservoir pada closed reservoir system sudah “terasa” yaitu gangguan akibat aktivitas produksi sudah sampai di batas reservoir. Kondisi pseudosteady state ini terkait dengan keadaan reservoir terbatas (finite-bounded) yang mempunyai Seminar Nasional Fisika 2016 Prodi Pendidikan Fisika dan Fisika, Fakultas MIPA, Universitas Negeri Jakarta

SNF2016-CIP-128


Gamba 2. Hasil plot distribusi tekanan transient Dari grafik diatas terlihat bahwa saat sumur diproduksi dengan laju konstan maka terjadi penurunan tekanan. Hal ini analogi dengan keadaan suatu reservoir saat diawal mulai diproduksi dan penurunan tekanan dibatas reservoir belum terasa. Selanjutnya dengan data yang sama diterapkan pada persamaan periode pseudosteady-state maka diperoleh grafik sebagai berikut:

p-ISSN: 2339-0654 e-ISSN: 2476-9398


Pseudosteady state artinya tekanan di setiap titik di reservoir menurun terhadap waktu dengan laju penurunan konstan Dari grafik tersebut terlihat bahwa penurunan tekanan di lubang sumur produksi munurun secara konstan selama rentang wajktu tertentu. Hal tersebut bersesuai dengan definisi Pseudosteady state artinya tekanan di setiap titik di reservoir menurun terhadap waktu dengan laju penurunan konstan Dengan memodifikasi solusi untuk periode yang dijelaskan diatas dengan memasukan faktor kerusakan disekitar lubang sumur (skin faktor), dan faktor volume resesrvoir maka hasilnya dapat digunakan untuk menganalisa beberapa parameter reservoir dalam uji alir yang disajikan pada tabel 1 berikut ini.

Gambar 3. Hasil plot distribusi tekanan pseudosteady-state Tabel 1. Pengelompokan solusi dan peneranpan persamaan difusivitas pada reservoir Priode Parameter Yang di No. Hasil modifikasi persamaan Prediksi 162,6qB permeabilitas k mh

1.

Transient

Faktor kerusakan disekitar lubang sumur (skin faktor)

Bentuk reservoir (Dietz Shape faktor (CA))

 pi

s  1,151



CA

Tekanan rata-rata reservoir

  2  rw exp ln 4 A    

p  p wf  141, 2

Laju Produksi

qB kh





2 p  p wf 141, 2

qB kh

  re  ln      rw 

k j h p  p wf



   2s      3  s 4 



  re  3   r   4    w 

141, 2 B ln 

Pseudosteady-state PI 

Productivity Index

kh

1  4A 141, 2 B  ln  2   2  C A rw FE 

Efesiensi aliran

Dalam analisis uji tutup (pressure drawdown test) di turunkan dengan prinsip superposisi (superposotion in time). Misalkan suatu sumur dialirkan dengan laju produksi konstan sebesar q.

  kt    3,23  cr 2    w 

 log 

m

q

2

 p wf

    s   

p  p wf  p s p  p wf

Pada watu t = tp, sumur kedua yang berlokasi sama dengan sumur pertama dialirkan dengan laju konstan sebesar –q (diinjeksikan) sementara sumur


SNF2016-CIP-129


pertama dibiarkan tetep mengalir dengan laju alir q. waktu pengaliran sumur kedua dinyatakan t. Ketika pengaruh kedua sumur dijumlahkan sebagai aplikasi dari superposisi, hasilnya adalah model yang menggambarkan sebuah sumur yang diproduksi pada laju q selama tp dan kemudian ditutup selama t. Oleh karena itu, untuk menganalisis data pressure buildup test digunakan aproksimasi logaritmik persamaan 22 untuk masing-masing sumur sehingga diperoleh : p ws  p i  162,6

qB kh

 t p  t   t   

log 

(25)

Dari hasil plot test buildup ini bisa digunakan untuk memprediksi permeabilitas dari kemiringan hrner plot: 162,6qB k

mh

Dengan cara yang sama seperti pada analisis data hasil pressure drawdowwn test, untuk menentukan faktor skin dapat diperoleh:

 p ws  p wf   1,696crw 2       log     m k  t     s  1,151   t p  t       log  t    p    4. Simpulan


p-ISSN: 2339-0654 e-ISSN: 2476-9398

Ucapan Terimakasih Terimakasih kepada Institut Teknologi bandung yang menjadi wadah sekaligus tempat menimba ilmu, terimakasih juga kepada para dosen dan rekan-rekan yang membantu penyelesaian penelitian ini.

Daftar Acuan Jurnal [1] Odeh, A.S.: “Pressure Drawdown Analysis, Variabel-Rate Case” Journal Of Petroleum Technology, vol 960. 08-1995 Buku [2] John lee. Well testing. First Printing. New York, Macmillan (1982), [3] Mattews, C.S. and Russell, D.G.Pressure Buildup and Flow Tests in Wells SPE Monograph vol.1, Henry L. Doherty series SPE of AIME, N.Y. 1967 [4] Abramowitz, M. & Stegun, I. A. (eds.) Handbook of Mathematical Functions Dover Publications, inc., NY,

Dari pembahasan diatas maka dapat disimpulkan bahwa dari solusi persamaan difusi dapat digunakan untuk memeperediksi nilai parameter reservoir diantaranya permeabilitas, faktor “skin” (kerusakan sekitar sumur), volume reservoir, Dietz Shape faktor dan tekanan rata-rata dan beberapa faktor produksi lain seperti productivity index dan efesiensi aliran .


SNF2016-CIP-130

Institut Teknologi Bandung, Jl. Ganesha No. 10, Bandung Abstrak. Abstract

Recommend Documents