MOD DEL SEBA ARAN TU UMPAHA AN MINYA AK DI AL LUR PELAY YARAN PE ELABUHAN TANJJUNG INT TAN CILA ACAP, WA TENG GAH JAW
Oleh: R Rizka Safitrii C C64104026
PRO OGRAM ST TUDI ILMU U DAN TEK KNOLOGI KELAUTA AN FAKULTA AS PERIKA ANAN DAN N ILMU KEL LAUTAN IN NSTITUT PERTANIA P AN BOGOR R 2009
PERNYATAAN MENGENAI SKRIPSI DAN SUMBER INFORMASI Dengan ini saya menyatakan bahwa Skripsi yang berjudul:
MODEL SEBARAN TUMPAHAN MINYAK DI ALUR PELAYARAN PELABUHAN TANJUNG INTAN CILACAP, JAWA TENGAH adalah benar merupakan hasil karya sendiri dan belum diajukan dalam bentuk apa pun kepada perguruan tinggi mana pun. Semua sumber data dan informasi yang berasal atau dikutip dari karya yang diterbitkan maupun tidak diterbitkan dari penulis lain telah disebutkan dalam teks dan dicantumkan dalam Daftar Pustaka dibagian akhir Skripsi ini. Bogor, Agustus 2009
RIZKA SAFITRI C64104026
RINGKASAN RIZKA SAFITRI. Model Sebaran Tumpahan Minyak di Alur Pelayaran Pelabuhan Tanjung Intan Cilacap, Jawa Tengah. Dibimbing oleh I WAYAN NURJAYA dan ANDRI PURWADANI. Ramainya alur pelayaran Cilacap setelah didirikannya kilang minyak tahun 1983 menyebabkan risiko pencemaran minyak akibat aktivitas pelayaran semakin meningkat. Berdasarkan kondisi tersebut maka perlu dilakukan tindakan penanggulangan jika terjadi tumpahan minyak dengan membuat studi yang dituangkan dalam sebuah model. Penelitian dengan topik pemodelan sebaran tumpahan minyak ini bertujuan untuk memodelkan pola sebaran minyak pada beberapa daerah titik rawan kecelakaan maupun rawan kebocoran di alur pelayaran Cilacap.pada tahun 2007 Penelitian dilaksanakan di Laboratorium P-TISDA, BPPT Jakarta pada bulan April 2008 – Februari 2009. Skenario pemodelan dijalankan pada bulan Februari 2007 (musim barat) dan bulan Agustus 2007 (musim timur). Pemodelan dibuat dengan asumsi terjadi kebocoran yang berasal dari kapal tanker maupun dari limbah pelabuhan dengan menggunakan DHI Software Mike 21. Data masukan pemodelan hidrodinamika berupa data arah dan kecepatan angin serta data pasang surut perairan Cilacap tahun 2007. Minyak yang diasumsikan tumpah antara lain diesel, avtur, minyak mentah, dan aspal. Kondisi pemodelan yang diamati yaitu hidrodinamika perairan, proses pelapukan masing-masing jenis minyak.dan sebaran total lapisan minyak pada saat perairan berada dalam posisi pasang tertinggi, surut terendah, menjelang pasang dan menjelang surut pada saat muka laut berada dalam posisi MSL. Data angin masukan model dan data hasil pengukuran lapang memiliki pola bertiup yang cenderung sama. Meskipun demikian, kecepatan angin yang digunakan sebagai masukan pemodelan cenderung lebih besar dari hasil pengukuran lapang sehingga memperbesar pengaruh angin dalam pemodelan. Data pasang surut masukan model dan data insitu memiliki fase yang sama sehingga pola sebaran lapisan minyak yang keluar masuk domain pemodelan sesuai dengan kondisi pasang surut di perairan Cilacap yang sebenarnya. Hasil pemodelan menunjukkan pola sebaran lapisan minyak di perairan Cilacap sangat dipengaruhi oleh resultan antara gaya yang diberikan oleh arus pasang surut dan angin permukaan. Perairan di sekitar pesisir Pulau Nusakambangan menjadi daerah yang paling rawan terkena dampak pencemaran pada musim barat karena pergerakan arus dan arah angin pada musim tersebut mengarah ke bagian tenggara dan timur domain. Tingkat kerawanan pencemaran minyak pada perairan Cilacap bersifat sementara serta high recovery dikarenakan sebagian besar lapisan minyak cenderung menyebar meninggalkan domain menuju Samudera Hindia. Pada musim timur, daerah pantai timur Cilacap serta di sepanjang aliran kanal utama dan Kali Donan memiliki risiko tertinggi terhadap pencemaran minyak karena arah arus dan angin permukaan yang bertiup dominan mengarah ke bagian barat dan barat laut domain. Tingkat kerawanan pencemaran minyak di Cilacap pada musim timur lebih tinggi dan lebih persistent dibandingkan pada musim barat.
Proses evaporasi lapisan avtur dan minyak mentah cenderung lebih tinggi pada musim timur, sedangkan lapisan diesel mengalami proses evaporasi tertinggi pada musim barat. Seluruh jenis minyak mengalami proses disolusi tertinggi pada musim barat, kecuali pada minyak mentah. Proses emulsifikasi lapisan avtur dan minyak mentah memiliki nilai yang lebih tinggi pada musim timur. Pada musim yang sama, aspal lebih cepat keluar dari dalam domain karena terdorong oleh arus dan angin permukaan sehingga proses emulsifikasi pada musim tersebut lebih rendah dibandingkan dengan musim barat. Proses emulsifikasi pada diesel terjadi pada seluruh lapisan minyak yang mengalami akumulasi. Seluruh jenis minyak mengalami proses dispersi vertikal tertinggi pada musim barat dibandingkan pada musim timur. Nilai exceedance frequency lapisan minyak (kecuali aspal) pada musim timur lebih tinggi jika dibandingkan dengan musim barat dikarenakan minyak yang tumpah pada musim tersebut banyak terakumulasi di dalam domain perairan. Nilai time exposition lapisan minyak lebih pendek pada musim timur jika dibandingkan dengan musim barat dikarenakan pada musim tersebut lapisan minyak lebih cepat hilang dari domain perairan akibat terbawa oleh arus dan terdorong oleh angin permukaan yang cukup besar. Nilai time exposition paling besar terjadi pada lapisan diesel dikarenakan diesel diasumsikan memasuki perairan secara konstan sehingga tetap berada dalam domain. Kesimpulan yang di dapat dari penelitian ini yaitu kondisi hidrodinamika dan tumpahan minyak yang dimodelkan dalam studi ini dapat diterima dengan baik karena data masukan pemodelan dan data hasil pengukuran insitu memiliki kemiripan tinggi. Pada musim timur, lapisan minyak cenderung menyebar ke arah barat domain perairan, sementara pada musim barat lapisan tersebut menyebar ke arah sebaliknya. Proses evaporasi pada lapisan minyak yang lebih tebal semakin meningkat seiring dengan peningkatan luas permukaan akibat dari turbulensi, sedangkan proses evaporasi pada lapisan minyak yang lebih tipis lebih banyak mendapat pengaruh dari suhu lingkungan sekitar. Proses disolusi pada musim barat lebih tinggi dibandingkan dengan musim timur diduga karena pada musim timur fraksi minyak dengan berat molekul rendah (aromatic) cenderung lebih dahulu terevaporasi daripada terdisolusi. Proses emulsifikasi seluruh jenis minyak memiliki nilai tertinggi pada musim timur dibandingkan dengan musim barat diduga karena pada musim timur lapisan minyak tersebut lebih banyak mengalami turbulensi. Proses dispersi vertikal pada musim barat justru lebih tinggi dibandingkan pada musim timur diduga karena pada musim timur lapisan minyak telah mengalami peningkatan viskositas yang cukup besar akibat proses emulsifikasi. Laju perubahan konsentrasi fraksi minyak bernilai tinggi di sekitar sumber tumpahan dan di bagian pusat lapisan karena lapisan minyak memiliki ketebalan yang lebih besar pada bagian lapisan tersebut.
MODEL SEBARAN TUMPAHAN MINYAK DI ALUR PELAYARAN PELABUHAN TANJUNG INTAN CILACAP, JAWA TENGAH
Skripsi Sebagai salah satu syarat untuk memperoleh gelar Sarjana Perikanan pada Fakultas Perikanan dan Ilmu Kelautan Institut Pertanian Bogor
Oleh: Rizka Safitri C64104026
PROGRAM STUDI ILMU DAN TEKNOLOGI KELAUTAN FAKULTAS PERIKANAN DAN ILMU KELAUTAN INSTITUT PERTANIAN BOGOR 2009
© Hak cipta milik Rizka Safitri, tahun 2009 Hak cipta dilindungi Dilarang mengutip dan memperbanyak tanpa izin tertulis dari Institut Pertanian Bogor, sebagian atau seluruhnya dalam bentuk apapun, baik cetak, fotocopy, microfilm, dan sebagainya
Judul
: MODEL SEBARAN TUMPAHAN MINYAK DI ALUR PELAYARAN PELABUHAN TANJUNG INTAN CILACAP, JAWA TENGAH Nama : Rizka Safitri NRP : C64104026
Disetujui, Pembimbing I
Pembimbing II
Dr. Ir. I Wayan Nurjaya, M.Sc NIP. 19640801 198903 1 001
Ir. Andri Purwandani NIP. 680 003 898
Mengetahui, Dekan Fakultas Perikanan dan Ilmu Kelautan
Prof. Dr. Ir. Indra Jaya, M.Sc NIP. 19610410 198601 1 002
Tanggal Lulus: 5 Agustus 2009
KATA PENGANTAR Segala puji dan syukur penulis panjatkan kepada Allah SWT atas seluruh rahmat dan karunia yang terlimpah bagi hambanya hingga saat ini. Shalawat dan salam tak lupa penulis haturkan pula bagi Rasul tercinta Nabi Muhammad SAW yang telah menjadi panutan dan tauladan bagi umatnya. Skripsi yang berjudul Model Sebaran Tumpahan Minyak di Alur Pelayaran Pelabuhan Tanjung Intan Cilacap, Jawa Tengah telah diselesaikan oleh penulis sebagai salah satu syarat kelulusan dalam menyelesaikan program studi S1 di Departemen Ilmu dan Teknologi Kelautan, Fakultas Perikanan dan Ilmu Kelautan Institut Pertanian Bogor. Dalam pembuatan skripsi ini, penulis menyampaikan terima kasih yang sebesar-besarnya atas waktu dan bimbingan dari Dr. Ir. I. Wayan Nurjaya, M.Sc selaku pembimbing utama dan kepada Ir. Andri Purwandani dari P3TISDA-BPPT selaku pembimbing anggota. Ucapan terima kasih penulis sampaikan pula kepada Prof. Dr. Ir. Harpasis S. Sanusi, M.Sc selaku dosen penguji dan kepada Dr. Ir. Henry M. Manik, MT selaku ketua komisi pendidikan di Departemen Ilmu dan Teknologi Kelautan. Besar harapan penulis, semoga hasil penelitian yang telah dilakukan dapat memberi manfaat bagi mahasiswa di Departemen Ilmu dan Teknologi Kelautan pada khususnya, serta masyarakat luas pada umumnya.
Bogor, Agustus 2009
Penulis
viii
UCAPAN TERIMA KASIH Atas terselesaikannya skripsi ini penulis mengucapkan terima kasih kepada: 1. Allah S.W.T atas rahmat dan Kasih-Nya, sehingga penulis dapat menghadapi segala permasalahan yang dihadapi. 2. Ibu, Adik dan keluarga di Solo atas kasih sayang, dukungan, dan doanya. 3. Dosen pembimbing skripsi, Dr. Ir. I Wayan Nurjaya, M.Sc. atas perhatian, bantuan dan bimbingan sehingga skripsi ini dapat terselesaikan. 4. Dosen pembimbing skripsi, Ir. Andri Purwandani beserta keluarga atas bantuan, bimbingan, saran, dan kritik selama proses penelitian. 5. Dosen Oseanografi Fisika Dr. Ir. Mulia Purba, M.Sc dan Tri Hartanto, S.Pi atas izinnya untuk menggunakan fasilitas di laboratorium Oseanografi Fisika. 6. Anugerah Trihatmojo atas semua waktu, tenaga, semangat, doa, hiburan, teknologi, pemikiran, pengetahuan, harapan, kasih sayang dan kepercayaan kepada penulis. 7. Seluruh teman-teman di Departemen Ilmu dan Teknologi Kelautan Institut Pertanian Bogor angkatan 2004 atas dukungan dan bantuan kepada penulis.
ix
DAFTAR ISI Halaman DAFTAR TABEL .........................................................................................
xiii
DAFTAR GAMBAR ......................................................................................
xiv
DAFTAR LAMPIRAN ................................................................................
xviii
1. PENDAHULUAN.....................................................................................
1
1.1. Latar Belakang ..................................................................................
1
1.2. Tujuan ...............................................................................................
2
2. TINJAUAN PUSTAKA ...........................................................................
3
2.1. Kondisi Umum Alur Pelayaran Cilacap ............................................
3
2.2. Arus ...................................................................................................
3
2.3. Angin .................................................................................................
5
2.4. Pasang Surut ......................................................................................
6
2.5. Tumpahan Minyak ............................................................................
7
2.5.1. Karakteristik Minyak ...........................................................
7
2.5.2. Sumber Pencemaran Minyak ...............................................
9
2.5.3. Perilaku Minyak di Laut ......................................................
11
2.6. Pemodelan Tumpahan Minyak .........................................................
17
3. METODOLOGI .......................................................................................
18
3.1. Waktu dan Lokasi .............................................................................
18
3.2. Sumber Data ......................................................................................
19
3.2.1. Data Masukan Model...........................................................
19
3.2.2. Data Verifikasi .....................................................................
20
3.3. Peralatan yang Digunakan ................................................................
20
3.4. Desain Skenario Model .....................................................................
20
3.4.1. Lokasi Pemodelan................................................................
22
3.4.2. Syarat Batas .........................................................................
23
3.4.3. Waktu Pemodelan ................................................................
24
3.4.4. Skenario Tumpahan Minyak ...............................................
24
x
3.5. Parameter Pemodelan ........................................................................
27
3.5.1. Parameter Hidrodinamika ....................................................
27
3.5.2. Parameter Spill Analysis ......................................................
35
3.6. Persamaan Utama ..............................................................................
40
3.7. Parameter Oil Spill ............................................................................
41
3.7.1. Spreading .............................................................................
41
3.7.2. Evaporation .........................................................................
42
3.7.3. Vertical Dispersion ..............................................................
44
3.7.4. Dissolution ...........................................................................
45
3.7.5. Emulsification ......................................................................
46
3.7.6. Heat transport......................................................................
47
3.7.7. Sifat Fisik dan Kimia Minyak .............................................
52
4. HASIL DAN PEMBAHASAN ................................................................
55
4.1. Verifikasi Hasil Pemodelan...............................................................
55
4.1.1.
Verifikasi Angin .................................................................
55
4.1.1.1. Musim Barat .........................................................
55
4.1.1.2. Musim Timur ........................................................
58
Verifikasi Pasang Surut ......................................................
61
4.1.2.1. Musim Barat .........................................................
61
4.1.2.2. Musim Timur ........................................................
62
4.2. Hasil Pemodelan Hidrodinamika ......................................................
64
4.1.2.
4.2.1.
Musim Barat ........................................................................
64
4.2.1.1. Menjelang Pasang (Muka Laut pada Kondisi MSL) 64
4.2.2.
4.2.1.2. Pasang ...................................................................
66
4.2.1.3. Menjelang Surut (Muka Laut pada Kondisi MSL)
67
4.2.1.4. Surut......................................................................
68
Musim Timur .......................................................................
69
4.2.2.1. Menjelang Surut (Muka Laut pada Kondisi MSL)
69
4.2.2.2. Surut......................................................................
71
4.2.2.3. Menjelang Pasang (Muka Laut pada Kondisi MSL) 72 4.2.2.4. Pasang ...................................................................
73
4.3. Hasil Pemodelan Pola Sebaran Total Minyak...................................
75
4.3.1.
Musim Barat ........................................................................
76
4.3.1.1. Kondisi Awal ........................................................
76
xi
4.3.1.2. Menjelang Pasang (Muka Laut pada Kondisi MSL) 77
4.3.2.
4.3.1.3. Pasang ...................................................................
78
4.3.1.4. Menjelang Surut (Muka Laut pada Kondisi MSL)
80
4.3.1.5. Surut......................................................................
81
Musim Timur .......................................................................
83
4.3.2.1. Kondisi Awal ........................................................
83
4.3.2.2. Menjelang Surut (Muka Laut pada Kondisi MSL)
84
4.3.2.3. Surut......................................................................
86
4.3.2.4. Menjelang Pasang (Muka Laut pada Kondisi MSL) 87 4.3.2.5. Pasang ...................................................................
88
4.4. Pembahasan Pola Sebaran Tumpahan Minyak .................................
90
4.5. Hasil Pemodelan Proses Pelapukan Tumpahan Minyak di Laut ......
93
4.5.1.
4.5.2.
Musim Barat ........................................................................
93
4.5.1.1. Jam ke-1 ................................................................
93
4.5.1.2. Jam ke-12 ..............................................................
95
4.5.1.3. Jam ke-24 ..............................................................
97
4.5.1.4. Jam ke-96 ..............................................................
98
Musim Timur ....................................................................... 100 4.5.2.1. Jam ke-1 ................................................................ 100 4.5.2.2. Jam ke-12 .............................................................. 102 4.5.2.3. Jam ke-24 .............................................................. 104
4.6. Pembahasan Proses Pelapukan Tumpahan Minyak di Laut............. 106 5. KESIMPULAN DAN SARAN ................................................................ 117 6. DAFTAR PUSTAKA ............................................................................... 119 7. LAMPIRAN .............................................................................................. 121 8. RIWAYAT HIDUP .................................................................................. 137
xii
DAFTAR TABEL
Halaman 1. Fraksi-Fraksi Minyak Bumi Berdasarkan Titik Didih ...............................
8
2. Informasi Lokasi, Jumlah Tumpahan dan Waktu Pengeluaran Skenario Model Tumpahan Minyak di Perairan Cilacap ..........................................
26
3. Informasi Spasial, Jumlah dan Waktu Tumpahan Masing-Masing Jenis Minyak yang di Skenariokan Tumpah di Perairan Cilacap .......................
35
4. Informasi Nilai Konstanta Transfer Bahang Minyak .................................
37
5. Informasi Nilai Konstanta Emulsifikasi Masing-Masing Minyak Dalam Skenario Model Tumpahan Minyak ..........................................................
38
6. Volume Fraksi Masing-Masing Minyak yang Diasumsikan Tumpah di Perairan Cilacap .........................................................................................
39
7. Fraksi Minyak Berdasarkan Struktur Kimia ..............................................
53
8. Perbandingan Pola Sebaran Total Lapisan Diesel, Avtur, Minyak Mentah, dan Aspal pada Berbagai Kondisi Muka Laut saat Musim Barat dan Musim Timur di Perairan Cilacap Tahun 2007 ..................................
91
9. Perbandingan Proses Evaporasi dan Disolusi Seluruh Jenis Minyak yang Dimodelkan Tumpah di Perairan Cilacap pada Jam Ke-12 pada Musim Barat dan Timur Tahun 2007 ..................................................................... 109 10. Perbandingan Proses Emulsifikasi dan Dispersi Vertikal Seluruh Jenis Minyak yang Dimodelkan Tumpah di Perairan Cilacap pada Jam Ke-12 pada Musim Barat dan Timur Tahun 2007 ................................................ 112 11. Perbandingan Exceedance Frequency dan Time Exposure Seluruh Jenis Minyak yang Dimodelkan Tumpah di Perairan Cilacap pada Jam Ke-12 pada Musim Barat dan Timur Tahun 2007 ................................................ 116
xiii
DAFTAR GAMBAR
Halaman 1. Proses Pelapukan Lapisan Minyak yang Tumpah di Permukaan Laut ......
11
2. Tingkat Evaporasi Berbagai Jenis Minyak Pada Suhu 15oC .....................
14
3. Proses Emulsifikasi pada Lapisan Minyak yang Membentuk "Chocolate Mousse" ......................................................................................................
15
4. Peta Lokasi Penelitian Model Sebaran Tumpahan Minyak di Perairan Cilacap, Jawa Tengah (Sumber: Google Earth, 2008) ...............................
18
5. Diagram Alir Pemodelan Sebaran Tumpahan Minyak di Perairan Cilacap dengan Menggunakan DHI Software Mike 21 ..............................
22
6. Domain Dasar Pemodelan Tumpahan Minyak di Perairan Cilacap dengan Menggunakan Program Mike21 ....................................................
23
7. Lokasi Skenario Sumber Tumpahan Minyak di Domain Perairan Cilacap .......................................................................................................
25
8. Batimetri Perairan Cilacap Hasil Survey Sounding Dasar Laut (Sumber: JANHIDROS, 2007) ..................................................................
28
9. Syarat Batas Terbuka pada Domain Model Hidrodinamika di Perairan Cilacap .......................................................................................................
29
10. Tinggi Muka Air Laut pada Seluruh Batas Terbuka Domain Perairan Cilacap pada Musim Barat Tahun 2007 .....................................................
30
11. Tinggi Muka Air Laut pada Seluruh Batas Terbuka Domain Perairan Cilacap pada Musim Timur Tahun 2007 ...................................................
30
12. Lokasi Pengamatan Data Pasang Surut Hasil Pengukuran Lapang dengan Data Masukan Model di Cilacap Tahun 2007 ...............................
31
13. Pola Nilai Tahanan Dasar (Manning Number) dalam Domain Model Perairan Cilacap .........................................................................................
33
14. Lokasi Pengamatan Data Angin Hasil Insitu dan Data Angin Masukan Model di Cilacap Tahun 2007 ....................................................................
34
15. Transfer Bahang Antara Udara, Lapisan Minyak, dan Air Laut ................
47
xiv
16. Arah [°] dan Kecepatan Angin [m/s] Masukan Model pada Musim Barat 2007 di Perairan Cilacap .................................................................. 17. Arah [°] dan Kecepatan Angin [m/s] Insitu pada Musim Barat 2007 di Perairan Cilacap .........................................................................................
55 56
18. Windrose Arah dan Kecepatan Angin Masukan Model dan Insitu pada Musim Barat 2007 di Perairan Cilacap ......................................................
57
19. Arah dan Kecepatan Angin Masukan Model pada Musim Timur 2007 di Perairan Cilacap .........................................................................................
58
20. Arah dan Kecepatan Angin Insitu pada Musim Timur 2007 di Perairan Cilacap .......................................................................................................
59
21. Windrose Arah dan Kecepatan Angin Masukan Model dan Insitu pada Musim Timur 2007 di Perairan Cilacap .....................................................
59
22. Pola Scattering Data Angin Masukan Model dan Insitu di Perairan Cilacap pada Musim Timur 2007 ...............................................................
61
23. Grafik Tinggi Muka Air Laut Hasil Prediksi Model pada Musim Barat 2007 di Perairan Cilacap ............................................................................
62
24. Grafik Tinggi Muka Air Laut Hasil Prediksi Model pada Musim Timur 2007 di Perairan Cilacap ............................................................................
62
25. Grafik Tinggi Muka Air Laut Hasil Pengukuran Insitu pada Musim Timur 2007 di Perairan Cilacap .................................................................
63
26. Perbandingan Fluktuasi Tinggi Muka Air Laut Hasil Pemodelan dan Tinggi Muka Air Laut Hasil Pengukuran Insitu di Perairan Cilacap pada Bulan Agustus 2007 ..........................................................................
63
27. Kondisi Hidrodinamika Perairan Cilacap Menjelang Pasang (Muka Laut Berada pada Kondisi MSL) pada Bulan Februari 2007 .............................
65
28. Kondisi Hidrodinamika Perairan Cilacap saat Pasang pada Bulan Februari 2007 .............................................................................................
66
29. Kondisi Hidrodinamika Perairan Cilacap Menjelang Surut (Muka Laut Berada pada Kondisi MSL) pada Bulan Februari 2007 .............................
68
30. Kondisi Hidrodinamika Perairan Cilacap saat Surut pada Bulan Februari 2007 ............................................................................................................
69
31. Kondisi Hidrodinamika Perairan Cilacap Menjelang Surut (Muka Laut Berada pada Kondisi MSL) pada Bulan Agustus 2007 .............................
70
xv
32. Kondisi Hidrodinamika Perairan Cilacap saat Surut pada Bulan Agustus 2007 ............................................................................................................ 33. Kondisi Hidrodinamika Perairan Cilacap Menjelang Pasang (Muka Laut Berada pada Kondisi MSL) pada Bulan Agustus 2007 .............................
71 72
34. Kondisi Hidrodinamika Perairan Cilacap saat Pasang pada Bulan Agustus 2007..............................................................................................
74
35. Pola Sebaran Total Lapisan Avtur saat Kondisi Awal di Perairan Cilacap pada Bulan Februari 2007 .............................................................
76
36. Pola Sebaran Total Lapisan Avtur Menjelang Pasang (Muka Laut pada Kondisi MSL) di Perairan Cilacap pada Bulan Februari 2007 ..................
78
37. Pola Sebaran Total Lapisan Avtur saat Pasang di Perairan Cilacap pada Bulan Februari 2007...................................................................................
79
38. Pola Sebaran Total Lapisan Avtur Menjelang Surut (Muka Laut pada Kondisi MSL) di Perairan Cilacap pada Bulan Februari 2007 ..................
81
39. Pola Sebaran Total Lapisan Avtur saat Surut di Perairan Cilacap pada Bulan Februari 2007...................................................................................
82
40. Pola Sebaran Total Lapisan Avtur saat Kondisi Awal di Perairan Cilacap pada Bulan Agustus 2007 ..........................................................................
84
41. Pola Sebaran Total Lapisan Avtur Menjelang Surut (Muka Laut pada Kondisi MSL) di Perairan Cilacap pada Bulan Agustus 2007...................
85
42. Pola Sebaran Total Lapisan Avtur saat Surut di Perairan Cilacap pada Bulan Agustus 2007 ...................................................................................
87
43. Pola Sebaran Total Lapisan Avtur Menjelang Pasang (Muka Laut pada Kondisi MSL) di Perairan Cilacap pada Bulan Agustus 2007...................
88
44. Pola Sebaran Total Lapisan avtur saat Pasang di Perairan Cilacap pada Bulan Agustus 2007 ...................................................................................
89
45. Sebaran serta Proses Pelapukan Lapisan Avtur Jam ke-1 di Perairan Cilacap pada Bulan Februari 2007 .............................................................
94
46. Sebaran serta Proses Pelapukan Lapisan Avtur Jam ke-12 di Perairan Cilacap pada Bulan Februari 2007 .............................................................
96
47. Sebaran serta Proses Pelapukan Lapisan Avtur Jam ke-24 di Perairan Cilacap pada Bulan Februari 2007 ............................................................... 98
xvi
48. Sebaran serta Proses Pelapukan Lapisan Avtur Jam ke-96 di Perairan Cilacap pada Bulan Februari 2007 ............................................................... 99 49. Sebaran serta Proses Pelapukan Lapisan Avtur Jam ke-1 di Perairan Cilacap pada Bulan Agustus 2007 ............................................................. 102 50. Sebaran serta Proses Pelapukan Lapisan Avtur Jam ke-12 di Perairan Cilacap pada Bulan Agustus 2007 ............................................................. 103 51. Sebaran serta Proses Pelapukan Lapisan Avtur Jam ke-24 di Perairan Cilacap pada Bulan Agustus 2007 ............................................................. 105
xvii
DAFTAR LAMPIRAN
Halaman 1. Contoh Laporan Hasil Pemodelan Hidrodinamika Pada Musim Barat dengan Menggunakan DHI Software Mike 21Hydrodynamic Modul ........ 122 2. Contoh Laporan Hasil Pemodelan Tumpahan Minyak Menggunakan DHI Software Mike 21 Spill Analysis Modul ..................................................... 129 3. Contoh Sumber Data Mentah Minyak Jenis Crude Oil ............................. 136
xviii
1. PENDAHULUAN
1.1.
Latar Belakang Cilacap merupakan satu-satunya daerah administratif di Propinsi Jawa
Tengah yang memiliki pelabuhan laut di pantai selatan Pulau Jawa. Sebagai daerah yang memiliki fasilitas pelabuhan ekspor dan impor, perairan Cilacap juga berfungsi sebagai jalur pelayaran/lalu lintas berbagai kapal pengangkut komoditi perdagangan. Ramainya alur pelayaran pelabuhan Cilacap menjadikan perairan tersebut sangat berpotensi mengalami pencemaran, khususnya pencemaran oleh minyak. Sejak peresmian Perluasan Kilang Minyak Cilacap tahun 1983, kegiatan di Pelabuhan Cilacap terus meningkat, terutama lalu lintas kapal-kapal tanker. Berbagai jenis minyak baik yang mentah maupun yang telah diolah diangkut dengan menggunakan kapal tanker. Seiring dengan ramainya lalu lintas kapal tersebut, berbagai macam kasus pencemaran akibat tumpahan minyak pernah terjadi di sekitar alur pelayaran. Hal tersebut membuktikan bahwa di beberapa tempat di sepanjang alur pelayaran ini terdapat area yang rawan terjadi kecelakaan. Kecelakaan yang dialami oleh kapal, khususnya kapal tanker, dapat menyebabkan terjadinya tumpahan minyak (oil spill) yang akan membawa dampak negatif terhadap lingkungan di sekitarnya. Maka dari itu, perlu dilakukan suatu tindakan untuk mencegah berulangnya kembali kasus serupa. Saat ini sudah menjadi suatu keharusan bagi perusahaan yang dalam kegiatan operasionalnya berpotensi mencemari lingkungan untuk membuat suatu kajian resiko (Risk Assessment). Tindakan dasar yang dilakukan dalam membuat kajian resiko adalah
1
2 meneliti daerah-daerah yang rawan kecelakaan serta mempelajari arah penyebaran minyak jika suatu saat terjadi kecelakaan kembali. Penelitian ini dibuat untuk memodelkan penyebaran tumpahan minyak dengan asumsi terjadi kecelakaan kapal tanker yang menyebabkan kebocoran minyak di beberapa tempat di alur pelayaran yang dianggap rawan kecelakaan. Selain itu dimodelkan pula kebocoran minyak yang berasal dari sumber-sumber lain di sekitar alur pelayaran Cilacap yang berpotensi mencemari lingkungan. Sifat dari sebagian jenis minyak cepat sekali menyebar ketika memasuki lautan. Untuk itu, jika tidak dilakukan penanganan yang tepat dikhawatirkan penyebaran lapisan minyak akan lebih meluas dan dampak yang ditimbulkan akan semakin besar. Diharapkan dari model ini, dihasilkan pemodelan pola sebaran tumpahan minyak di laut yang dapat mewakili kondisi sebenarnya untuk membantu proses penanganan pencemaran minyak secara cepat dan tepat.
1.2.
Tujuan Penelitian ini bertujuan untuk membuat skenario model pola sebaran
minyak pada beberapa daerah rawan kecelakaan maupun rawan kebocoran di alur pelayaran Cilacap. Skenario model dibuat dengan asumsi terjadi kebocoran yang berasal dari kapal tanker maupun dari kapal-kapal lainnya yang melewati alur pelayaran tersebut dengan menggunakan DHI Software Mike 21 modul Hydrodynamic dan modul Spill Analysis. Skenario model dimodelkan pada bulan Februari 2007 sebagai representatif pada musim barat dan pada bulan Agustus 2007 sebagai representatif pada musim timur.
2. TINJAUAN PUSTAKA
2.1.
Kondisi Umum Alur Pelayaran Cilacap Pelabuhan Tanjung Intan Cilacap merupakan satu-satunya pelabuhan di
pantai selatan Pulau Jawa yang merupakan pintu gerbang perekonomian bagi daerah Jawa Tengah bagian selatan untuk perdagangan ekspor dan impor maupun pasar antar Pulau. Selain memiliki dermaga umum, terdapat beberapa perusahaan besar di Cilacap yang memiliki pelabuhan khusus tersendiri di luar pelabuhan tersebut, seperti Pelabuhan Minyak Pertamina UP IV dan pelabuhan semen milik Holcim (Wikipedia, 2007). Alur pelayaran di sekitar pelabuhan mempunyai kedalaman rata-rata -11 m s/d -12 m LWS. 2.2.
Arus Arus laut yaitu proses pergerakan massa air laut menuju kesetimbangan
yang menyebabkan perpindahan horizontal dan vertikal massa air laut secara terus menerus. Berdasarkan asal penyebabnya, terdapat dua gaya yang berhubungan dengan arus yaitu gaya eksternal dan gaya internal. Gaya eksternal terdiri dari angin, perbedaan tekanan udara, gaya gravitasi, gaya tektonik serta gaya tarik benda-benda angkasa yang dipengaruhi oleh tekanan dasar laut. Gaya internal arus antara lain perbedaan densitas air laut, gradien tekanan mendatar dan gesekan lapisan air (Gross, 1972). Kecepatan arus permukaan sangat bergantung dari kecepatan dan lamanya angin bertiup. Kecepatan arus permukaan besarnya kurang dari 2% dari kecepatan angin (Gross, 1972). Arah pergerakan arus permukaan ini tidak searah
3
4 dengan arah pergerakan angin dikarenakan oleh adanya gaya Coriolis yang menyebabkan timbulnya perubahan arah arus sebesar ± 45° dari arah angin. Dinamika pasut akan menimbulkan perbedaan tekanan hidrostatik pada beberapa tempat sehingga dapat terjadi arus yang dikenal sebagai arus pasang surut (Gross, 1972). Arus pasang surut dengan tipe rotary dominan terdapat di laut terbuka dan di perairan dekat pantai, sedangkan arus tipe bolak-balik (reversing current) umum terjadi di perairan yang berbentuk terusan, selat dan alur pelayaran yang relatif sempit. Penyebaran lapisan minyak yang berada di permukaan laut sangat dipengaruhi oleh arus permukaan. Jika lapisan minyak dekat dengan daratan dimana kecepatan angin kurang dari 10 km/jam, maka lapisan tersebut 100% menyebar mengikuti arus permukaan. Pengaruh angin pada lapisan minyak dalam kondisi tersebut tidak lebih dari 3% (CRC, 2000). Melalui hasil pengukuran arus di alur pelayaran Cilacap 1992 oleh Dishidros TNI-AL pada tiga stasiun diperoleh hasil bahwa di perairan tersebut arus pasut lebih dominan dan arus nonpasut relatif lebih rendah (Dishidros,1992 in Firdaus, 1997). Pada saat pasang, massa air laut akan mengalir dari Samudera Hindia masuk melalui perairan antara Pulau Jawa dan Nusakambangan melalui pintu terusan timur (pantai Pulau Jawa) dan pintu terusan barat (Nusakambangan). Tetapi massa air yang masuk melalui pintu timur lebih dominan daripada yang masuk dari pintu barat. Pada saat surut, massa air mengalir kembali menuju Samudera Hindia melalui jalan yang sama. Periode aliran massa air keluar lebih panjang dibandingkan dengan periode aliran massa air yang masuk. Kecepatan arus permukaan di perairan Cilacap berkisar antara 3.5 knot sampai 4.0 knot
5 (Pertamina UP IV, 1992 in Harimurthy, 2001). Arus di sekitar alur pelayaran yang berasal dari hulu perairan Kali Donan sangat kecil, karena pada dasarnya perairan tersebut bukan perairan sungai (Ilham, 2002).
2.3.
Angin Angin didefinisikan sebagai gerakan udara mendatar (horizontal) yang
disebabkan oleh perbedaan tekanan udara antara dua tempat. Angin yang berhembus di permukaan perairan akan menimbulkan wind wave, yaitu gelombang yang ditimbulkan oleh angin. Peristiwa ini merupakan pemindahan tenaga angin menjadi tenaga gelombang di permukaan air dan gelombang itu sendiri meneruskan tenaganya kepada peristiwa lainnya, diantaranya molekul air. Selain menimbulkan gelombang di permukaan air, angin juga dapat menyebabkan terjadinya arus. Penyebaran lapisan minyak yang berada di permukaan laut dipengaruhi oleh angin permukaan. Jika kecepatan angin bertiup lebih besar dari 20 km/jam yang tentu saja terjadi pada laut terbuka, maka penyebaran lapisan minyak ditentukan oleh kondisi angin setempat. Hal tersebut tidak berlaku jika kecepatan angin kurang dari 10 km/jam dimana angin tidak memainkan peranan penting dalam proses penyebaran minyak (CRC, 2000). Pola angin yang sangat berperan di Indonesia adalah angin muson. Angin muson bertiup secara mantap ke arah tertentu pada satu masa sedangkan pada masa lainnya angin bertiup secara mantap pula dalan arah yang berlawanan. Bulan Desember, Januari dan Februari adalah musim dingin di belahan bumi utara dan musim panas di belahan bumi selatan. Pada saat itu terbentuklah pusat tekanan udara tinggi di atas daratan Asia dan pusat tekanan rendah di atas daratan
6 Australia. Keadaan ini menyebabkan angin berhembus dari Asia menuju Australia, yang di Indonesia umumnya dikenal sebagai angin muson barat. Sebaliknya, pada bulan Juli –Agustus terjadilah pusat tekanan tinggi di atas daratan Australia dan pusat tekanan udara rendah di atas daratan Asia sehingga mengakibatkan berhembusnya angin muson timur di Indonesia. Dua kali dalam setahun angin muson berganti arah. Kecepatan angin rata-rata bulanan di Cilacap pada umumnya bervariasi antara dua hingga enam knot [mil/jam]. Arah angin yang paling dominan atau yang paling sering adalah ke arah tenggara (Firdaus, 1997).
2.4.
Pasang Surut Pasang-surut atau pasut adalah proses naik turunnya paras laut (sea level)
secara berkala yang ditimbulkan oleh adanya kombinasi gaya sentrifugal dan gaya tarik dari benda-benda angkasa, terutama matahari dan bulan, terhadap massa air di bumi (Pariwono, 1989). Gerakan pasut menyebabkan permukaan air laut senantiasa berubah setiap saat. Periode selama permukaan air laut naik disebut air pasang (flood tide), sedangkan kedudukan saat permukaan air laut mencapai puncaknya disebut air tinggi (high water). Saat permukaan air laut menurun akibat gaya pasut disebut air surut (ebb tide) dan kedudukan permukaan air laut rendah disebut air rendah (low water). (Gross, 1972). Pasang surut di perairan Cilacap adalah penjalaran langsung pasang surut di Samudera Hindia yang bertipe campuran, dimana komponen setengah hariannya masih dominan (semi diurnal) (Pariwono, 1989). Dengan sifat tersebut, maka terjadi dua kali pasang dan dua kali surut dalam 24 jam yang tidak teratur
7 dengan beda pasang surut (tidal range) antara 1.5 sampai 2.0 meter (Pertamina UP IV, 1992 in Harimurthy, 2001).
2.5. Tumpahan Minyak 2.5.1. Karakteristik Minyak Minyak mentah (crude oil) adalah campuran kompleks hidrokarbon dengan jumlah atom karbon antara 4-26 atom dalam satu molekul. Susunan atom karbon dapat membentuk rantai lurus dan rantai cabang (alifatik), rantai siklik (alisiklik), dan rantai aromatik (Clark, 1986). Komponen hidrokarbon aromatik jumlahnya relatif kecil jika dibandingkan dengan komponen hidrokarbon alifatik dan alisiklik. Namun demikian, komponen aromatik justru lebih beracun, mudah berubah menjadi gas, dan menguap. Secara umum toksisitas minyak mentah meningkat dengan memanjangnya rantai karbon (Mukhtasor, 2007). Komposisi senyawa hidrokarbon dari minyak mentah berbeda-beda antar sumur minyak yang satu dengan yang lain, tergantung pada sumber penghasil minyak tersebut. Agar dapat digunakan, minyak mentah terlebih dahulu harus melewati proses penyulingan. Penyulingan (refining) adalah proses destilasi minyak mentah untuk memutuskan ikatan rantai karbon yang berbeda titik didihnya menjadi beberapa fraksi (Clark, 1986). Hasil pengelompokan fraksi minyak mentah berdasarkan titik didihnya disajikan dalam Tabel 1. Seluruh komponen dari minyak mentah dapat diuraikan oleh bakteri dengan tingkat kecepatan yang bervariasi. Minyak dengan komposisi rantai karbon yang sederhana, lurus, maupun bercabang dapat terurai dengan cepat. Minyak dengan komposisi molekuler yang rumit, berupa ter, atau minyak yang membentuk gumpalan membutuhkan waktu yang sangat lama untuk terurai
8 disebabkan luas permukaan minyak dalam bentuk tersebut lebih kecil jika dihubungkan dengan ukuran volumenya (Clark, 1986). Tabel 1. Fraksi-Fraksi Minyak Bumi Berdasarkan Titik Didih Fraksi
Titik Didih Ukuran [°C] Molekuler < 25 C3 - C4 Refinery gases Gasolin 40 – 150 C4 - C10 Naptha 150 – 200 C10 - C12 Kerosin 200 – 250 C12 - C18 Minyak gas 250 – 300 C18 - C20 Minyak pelumas 300 – 400 C20 - C25 Minyak sisa > 400 > C25 Sumber : Bishop (1983) dalam Mukhtasor (2007)
Volume [%] 2 25 6 10 15 17 25
Diesel merupakan pencampuran kompleks dari minyak hasil penyulingan. Diesel mudah menyala dengan titik didih antara 150 – 380 0C. Diesel biasa digunakan sebagai bahan bakar kendaraan bermotor. Bahan bakar mesin diesel sebagian besar terdiri dari senyawa hidrokarbon dan senyawa nonhidrokarbon. Senyawa hidrokarbon yang dapat ditemukan dalam bahan bakar diesel antara lain parafinik, naftenik, olefin dan aromatik. Sedangkan untuk senyawa nonhidrokarbon terdiri dari senyawa yang mengandung unsur non logam, yaitu S, N, O dan unsur logam seperti vanadium, nikel dan besi. Avtur adalah campuran minyak tanah dengan hidrokarbon cair. Avtur digunakan sebagai bahan bakar untuk pesawat terbang jet yang terdiri atas hidrokarbon sedang dengan karakteristik distilasi dan titik nyala seperti minyak tanah dan kandungan aromatik maksimum 25% terhadap volume. Kekentalan avtur di bawah 8 cST pada temperatur -20 derajat C dan titik beku di bawah -47
9 derajat C. Beberapa jenis komposisi avtur antara lain Paraffin, Olefin, Naptha, dan Aromatic (Pertamina, 2006).
2.5.2. Sumber Pencemaran Minyak Minyak masuk ke lingkungan perairan laut dengan beberapa cara, yaitu: a.
Eksplorasi Lepas Pantai Sumber pencemaran minyak yang berasal dari eksplorasi lepas pantai
cenderung kecil jika dibandingkan dengan jumlah total minyak yang masuk ke lingkungan laut. Namun, jika terjadi kecelakaan tertentu seperti semburan sumur minyak (well blow-out), kerusakan struktur platform, maupun kerusakan perlatan, maka sejumlah besar minyak dipastikan akan mencemari laut (Mukhtasor, 2007). b.
Transportasi Laut Polutan yang berasal dari transportasi laut dapat berasal dari
pengoperasian kapal dan tanker maupun kecelakaan kapal dan tanker. Dari beberapa sumber tersebut, input polutan terbesar berasal dari pengoperasian kapal tanker dalam proses deballasting (sistem kestabilan kapal menggunakan mekanisme bongkar-muat air). Air ballast adalah air laut yang diisikan ke dalam tanki sebuah tanker yang kosong dimana tanki tersebut sebelumnya merupakan wadah minyak mentah. Untuk mengisi kembali tanki tersebut dengan minyak, maka air ballast yang terdapat di dalamnya harus dibuang ke laut dengan membawa sisa-sisa minyak yang terdapat di dinding tanki (Clark, 1986). Jika dibandingkan dengan proses deballasting, polutan dari kecelakaan tanker hanya berkontribusi sangat kecil dari keseluruhan minyak yang masuk ke laut. Namun kecelakaan tanker tetap menjadi masalah yang besar karena menghasilkan buangan minyak yang volumenya relatif besar pada suatu lokasi.
10 Semakin besar ukuran tanker, maka diperkirakan input polutan minyak ke laut juga semakin besar. Konsentrasi polutan dalam jumlah besar tentunya dapat menyebabkan kerusakan lingkungan pada area tersebut (Mukhtasor, 2007). Sumber lapisan minyak lainnya yang berasal dari transportasi laut antara lain docking atau perawatan kapal. Dalam proses tersebut, semua sisa bahan bakar yang ada dalam tangki harus dikosongkan untuk mencegah terjadinya ledakan dan kebakaran. Selain itu proses bongkar muat tanker yang dilakukan di tengah laut juga banyak menimbulkan resiko tumpahan minyak akibat seperti pipa yang pecah, bocor maupun kecelakaan karena kesalahan manusia. Proses scrapping kapal (pemotongan badan kapal untuk menjadi besi tua) dapat pula menyebabkan banyak kandungan metal dan lainnya termasuk kandungan minyak terbuang ke laut. c.
Sumber dari Darat Input polutan yang berasal dari darat bersumber dari berbagai aktivitas
manusia, seperti pemakaian minyak untuk keperluan industri, limbah rumah tangga, kegiatan perbengkelan, kilang minyak, run off dari daerah perkotaan, maupun hasil pembakaran hidrokarbon di atmosfer yang terbawa oleh hujan. Limbah minyak tersebut terbawa oleh sistem saluran air yang menuju ke sungai dan bermuara ke laut. Apabila diakumulasi, jumlah limpasan minyak yang berasal dari darat menjadi sumber utama polutan minyak yang masuk ke kawasan pesisir dan laut (Clark, 1986). d.
Sumber Alami Laut merupakan tempat dimana minyak bumi secara alami akan
menyembur ke permukaan bumi di dasar laut dan merembes masuk ke perairan.
11 Sumber polutan dalam kasus ini merupakan suatu fenomena alami, meskipun total masukan polutan yang berasal dari rembesan tersebut kemungkinan jumlahnya dua kali lebih besar dari pada masukan polutan dari kecelakaan tanker (Clark, 1986).
2.5.3. Perilaku Minyak di Laut Minyak yang masuk ke dalam lingkungan laut akan mengalami berbagai proses, baik secara fisika maupun secara kimia (Gambar 1). Proses-proses tersebut antara lain membentuk lapisan (slick formation), menyebar (dissolution), menguap (evaporation), emulsifikasi (emulsification), minyak dalam air (oil in water emulsions), fotooksidasi (photooxidation), biodegradasi mikroba (microbial biodegradation), sedimentasi (sedimentation), dicerna oleh plankton (plankton ingestion), dan bentukan gumpalan ter (tur lump formation). Semua proses tersebut secara kolektif disebut dengan weathering of oil (Mukhtasor, 2007).
Gambar 1. Proses Pelapukan Lapisan Minyak yang Tumpah di Permukaan Laut (Sumber : ITOPF, 2007)
12 Penyebaran, penguapan, dispersi, emulsifikasi, dan pelarutan adalah proses-proses penting selama tahap awal tumpahan. Sementara oksidasi, sedimentasi, dan biodegradasi adalah proses weathering jangka panjang yang akan membantu proses penguraian minyak. Menurut Krough (1980) dalam Firdaus (1997), berdasarkan kekekalannya (persistent) tumpahan minyak dapat dibedakan menjadi dua bagian, yaitu tumpahan minyak yang tidak kekal (nonpersistent) dan tumpahan minyak yang kekal (persistent). Tumpahan minyak non-persistent akan berangsur-angsur menghilang dari permukaan laut akibat adanya proses fisika-kimia, sedangkan tumpahan minyak yang kekal (persistent) akan menyebar secara perlahan sehingga mencemari lingkungan laut. a.
Penyebaran (Spreading) Minyak yang keluar di permukaan air akan dengan segera bertambah luas
permukaannya. Mekanisme spreading dipengaruhi oleh karakteristik minyak itu sendiri antara lain perbedaan densitas minyak dan air laut, dan tegangan permukaan. Semakin rendah nilai viskositas minyak, maka minyak akan menyebar semakin cepat. Kecepatan dari penyebaran minyak serta ketebalan lapisannya tergantung dari suhu perairan dan jenis minyak yang tumpah (Clark, 1986). Proses penyebaran tumpahan minyak juga dipengaruhi oleh arus air, pola pasang surut, kecepatan angin, dan kekasaran muka laut ( Fay, 1971 in Mukhtasor, 2007 ). Angin dan arus pasang surut memindahkan unsur-unsur dari lapisan minyak secara relatif satu sama lain dan mempercepat proses penyebaran. Ketika lapisan membentuk gumpalan dengan luas permukaan yang stabil, hanya dispersi horizontal yang memindahkan unsur-unsur minyak menjauh dari pusat
13 massa. Untuk perairan yang tertutup dan estuari, pergerakan lapisan minyak lebih banyak mendapat pengaruh dari arus dan pasang surut setempat (DHI, 2006b). Dalam mekanisme spreading, minyak dapat menyebar secara horizontal meskipun tanpa angin. Proses penyebaran minyak disebabkan oleh gaya gravitasi dan tegangan permukaan antara minyak dan air. Gaya-gaya tersebut berlawanan dengan gaya yang diberikan oleh pengaruh viskositas minyak. b.
Penguapan (Evaporation) Proses penyebaran minyak akan menyebabkan lapisan menjadi lebih tipis
dan proses penguapan meningkat. Proses penguapan pada tumpahan minyak dipengaruhi oleh komposisi minyak, suhu udara dan air laut, area tumpahan, kecepatan angin, radiasi matahari, dan ketebalan lapisan minyak (DHI, 2006b). Secara umum, komponen dalam minyak dengan berat molekul rendah atau minyak yang memiliki titik didih di bawah 200°C cenderung mengalami penguapan dalam waktu 24 jam. Kekasaran muka laut, kecepatan angin, dan temperatur yang tinggi akan meningkatkan penguapan ( ITOPF, 2007). Selain itu luas permukaan minyak juga sangat berperan dalam proses ini. Sifat minyak dapat berubah secara signifikan seiring terjadinya proses penguapan. Hilangnya sebagian material yang bersifat mudah menguap mengakibatkan berat jenis minyak menjadi lebih berat. Berikut ditampilkan tingkat evaporasi berbagai jenis minyak (Gambar 2)
14
G Gambar 2. Tingkat T Evaaporasi Berbaagai Jenis Minyak M Pada Suhu 15oC (Sumber: CR RC, 2000) c c.
Disppersi vertikall Disppersi merupakkan proses mekanik. m Tuurbulensi airr memecahkaan lapisan
m minyak men njadi butiran dan memasuukannya ke dalam kolom m perairan, emulsi minyyak-dalam-aair. Butiran yang m membentuk y berukuuran sangat kecil k (<20 μ relatif lebih stabil bercampur μm) b deengan air lauut di kolom pperairan mennjadi b bentuk tersuuspensi. Butiran yang beerukuran bessar (>100 μm m) cenderung g akan k kembali naik k dan bergabbung ke lapissan minyak di permukaaan laut (ITOP PF, 2007). D Dispersi verrtikal bergantung pada siifat minyak dan d jumlah energi e laut. Kecepatan d dispersi minnyak akan sem makin tinggi jika viskossitas minyak rendah dan nilai k kekasaran m muka laut bessar (ITOPF, 2007). Dalaam cuaca yaang buruk, mekanisme m d dispersi lebih dominan disebabkan d ooleh pecah gelombang. S Sebaliknya pada p cuaca
15 yang tenang, mekanisme dispersi yang paling signifikan terjadi karena stretching compression dari lapisan, yang menyebabkan terbentuknya droplet (DHI, 2006b). d.
Emulsifikasi Salah satu proses penting yang menyebabkan bertahan lamanya minyak di
permukaan laut yaitu dengan membentuk emulsi air-dalam-minyak, yang mengubah minyak menjadi campuran yang sangat kental. Emulsi terbentuk jika terdapat dua cairan (liquid) yang bercampur, dimana salah satu dari cairan tersebut tersuspensi dalam cairan lainnya (Clark, 1986). Emulsi tersebut dapat menyerap hingga 80% air. Kestabilan dari bentuk ini sangat tergantung pada jenis minyak dan kondisi lingkungan. Kestabilan dari emulsi sangat berhubungan dengan jumlah kehadiran surfactant (resin dan aspal) dalam minyak, sedangkan tingkat pengambilan air sangat berhubungan dengan kondisi laut setempat seperti gelombang dan turbulensi air (DHI, 2006b). Dalam beberapa kondisi, emulsi akan membentuk lapisan tebal di permukaan laut dan berwujud kental yang disebut sebagai ”chocolate mousse” (Clark, 1986). Emulsi dapat terpisah kembali menjadi minyak dan air jika dipanaskan oleh sinar matahari pada kondisi permukaan laut yang tenang atau saat terdampar di pantai (ITOPF, 2007). Wujud dari emulsi minyak ditampilkan dalam Gambar 3.
Gambar 3. Proses Emulsifikasi pada Lapisan Minyak yang Membentuk "Chocolate Mousse"
16 Proses emulsifikasi akan mempengaruhi volume lapisan minyak serta meningkatkan viskositas minyak. Jika emulsifikasi minyak terdampar di pantai maka akan mengganggu kehidupan ekosistem di daerah tersebut. e.
Disolusi Komponen dari minyak yang dapat larut dalam air akan terlarut dalam
kolom perairan. Proses pelarutan tersebut akan cepat terjadi pada minyak yang telah lebih dulu terdispersi dalam air. Minyak umumnya hanya sedikit mengandung komponen yang dapat larut dalam air. Salah satu komponen yang paling cepat terlarut dalam air adalah hidrokarbon aromatik dengan berat jenis rendah dan komponen polar resin. Komponen-komponen yang dapat terdisolusi tersebut umumnya beracun. Meskipun demikian, senyawa aromatik biasanya akan lebih dulu menguap dibandingkan terlarut dikarenakan proses penguapan terjadi 10-100 kali lebih cepat dibandingkan proses melarut (ITOPF, 2007). Batas kadar minyak yang diizinkan berada di kolom perairan yaitu 0.01 ppm. f.
Sedimentasi Sedimentasi merupakan proses dimana minyak terdeposisi ke dasar laut.
Sedimentasi terjadi ketika butir minyak mencapai densitas tinggi dibandingkan dengan densitas air setelah berinteraksi dengan mineral tersuspensi di dalam kolom perairan. Minyak juga bereaksi terhadap oksigen dan menghasilkan bentuk persistent yang disebut ter (tars) akibat adanya proses oksidasi minyak dengan viskositas tinggi. Proses ini terjadi dengan sangat lambat pada lapisan minyak yang terekspos sinar matahari. Beberapa hasil dari oksidasi minyak ini memiliki densitas yang besar dan dapat tenggelam di air payau atau di perairan dangkal (ITOPF, 2007).
17 g.
Biodegradasi Dalam kolom perairan, terdapat beberapa jenis mikroorganisme yang bisa
menguraikan sebagian atau seluruh komponen minyak. Hasil penguraian tersebut dapat berupa komponen yang dapat terlarut dalam air atau terkadang berupa karbondioksida dan air. Unsur utama yang berpengaruh terhadap efisiensi proses biodegradasi yaitu nutrien (nitrogen dan fosfor), temperatur, dan oksigen terlarut. Minyak yang berbentuk butiran atau partikel lebih mudah mengalami proses biodegradasi karena menyediakan luas permukaan yang lebih besar dibandingkan minyak yang berbentuk lapisan tebal atau gumpalan (ITOPF, 2007). Tingkat biodegradasi tinggi pada rantai jenuh (12-20 atom karbon) dan lambat pada rantai aromatik dan aspal.
2.6.
Pemodelan Tumpahan Minyak Rau dan Woten (1980) dalam Firdaus (1997) menyatakan bahwa model
merupakan penampakan dari sistem yang sebenarnya. Perilaku dan konsentrasi polutan di laut dapat diperkirakan atau diestimasi menggunakan pemodelan dengan bantuan komputer. Karena umumnya perilaku maupun konsentrasi polutan di alam memiliki proses yang kompleks, maka pemodelan dapat dimanfaatkan untuk menyederhanakan proses tersebut. Pemodelan hidrodinamika dapat menjadi salah satu cara untuk mengetahui proses penyebaran polutan. Pemodelan yang akurat membutuhkan representasi yang baik mengenai parameter, proses, dan kondisi batas pemodelan. Secara umum, pemodelan perilaku dan penyebaran polutan terdiri dari dua komponen pokok, yaitu model hidrodinamika serta model perilaku dan penyebaran dari polutan itu sendiri (Mukhtasor, 2007).
3. METODOLOGI 3.1.
Waktu dan Lokasi Penelitian berupa pemodelan sebaran tumpahan minyak di Perairan
Cilacap, Jawa Tengah dilakukan pada bulan April 2008 – Februari 2009 menggunakan DHI software Mike 21 dengan modul Hydrodynamic dan Spill Analysis (Gambar 4). Pemodelan dilaksanakan dengan menggunakan perangkat komputer Laboratorium Pusat Teknologi & Inventarisasi Sumberdaya Alam (PTISDA) bertempat di Badan Pengkajian dan Penerapan Teknologi (BPPT) Jakarta.
Gambar 4. Peta Lokasi Penelitian Model Sebaran Tumpahan Minyak di Perairan Cilacap, Jawa Tengah (Sumber: Google Earth, 2008)
18
19 3.2.
Sumber Data Sumber data yang diperlukan dalam penelitian ini meliputi dua hal, yaitu
data untuk masukan model serta data untuk kepentingan verifikasi.
3.2.1. Data Masukan Model Untuk membangun skenario model, diperlukan beberapa data masukan yang didapat dari berbagai sumber, antara lain: a. Data kedalaman (batimetri) perairan Cilacap, yaitu: 1) peta batimetri hasil pemetaan Jawatan Hidro-Oseanografi (JANHIDROS) TNI-AL tahun 2007 Nomor 108 dengan skala 1 : 15.000; 2) peta batimetri hasil pemetaan PT. PERTAMINA (PERSERO) Unit Pengolahan IV Cilacap tahun 2006 Nomor Gambar CS 05/X/06 dan CS 07/IX/06 dengan skala 1 : 2000; 3) peta batimetri hasil survey sounding Kolam Pelabuhan Tanjung Intan Cilacap PT. (PERSERO) Pelabuhan Indonesia III Surabaya tahun 2006 Nomor Gambar DL 427/2006 dengan skala 1 : 250; b. Data arah dan kecepatan angin di Cilacap bulan Februari dan Agustus tahun 2007 dengan interval data per enam jam dan bersumber dari QuickScat & Seawind (IFREMER); c. Data pasang surut perairan Cilacap bulan Februari dan Agustus tahun 2007 dengan interval data per 15 menit bersumber dari Topex Poesidon & Jason; d. Data lalu lintas perkapalan dan rute/alur pelayaran tahun 2007 diperoleh dari PT. (PERSERO) Pelabuhan Indonesia III cabang Tanjung Intan Cilacap, Jawa Tengah;
20 e. Data oil properties dari Lembaga Minyak & Gas (LEMIGAS) Jakarta dan dari berbagai sumber (Lampiran 3).
3.2.2. Data Verifikasi Data yang diperlukan untuk verifikasi masukan skenario model antara lain: a. Data arah dan kecepatan angin di Cilacap bulan Februari dan Agustus tahun 2007 hasil pengukuran Badan Meteorologi dan Geofisika (BMG) Cilacap yang direkam setiap jam selama 28 hari; b. Data pasang surut perairan Cilacap bulan Februari dan Agustus tahun 2007 hasil pengukuran Badan Koordinasi Survey dan Pertanahan Nasional (BAKOSURTANAL) Cibinong dengan interval pengukuran data per 15 menit.
3.3.
Peralatan yang Digunakan Sistem perangkat keras yang dipakai dalam pemodelan maupun
pengolahan data masukan (input) yaitu menggunakan sistem perangkat komputer di BPPT. Pembuatan skenario pemodelan sebaran tumpahan minyak diproses dengan menggunakan berbagai modul, antara lain Mike Zero Bathymetries, Mike Zero Time Series, Mike Zero Profile Series, Mike Zero Data Extraction, Mike Zero Toolbox, dan Mike 21 Flow Model. Untuk skenario analisis tumpahan minyak, digunakan modul Hydrodynamic Modul dan Spill Analysis Modul.
3.4.
Desain Skenario Model Model diawali dengan pengolahan data masukkan untuk menyimulasikan
modul hidrodinamika pada program Mike 21. Data masukkan yang diolah antara lain pembuatan domain model dengan menggunakan data kedalaman perairan,
21 pengolahan data arah maupun kecepatan angin dari IFREMER yang dihitung tiaptiap grid dan berubah terhadap ruang dan waktu, serta data prediksi pasang surut yang dihasilkan dari satelit Topex Poseidon dan Jason. Data tersebut kemudian diverifikasi dengan menggunakan data hasil pengukuran lapang. Proses selanjutnya adalah membuat skenario pemodelan hidrodinamika dengan memasukkan data input angin dan pasang surut yang telah diverifikasi serta melengkapi data-data parameter pendukung dalam modul hidrodinamika tersebut. Modul hidrodinamika yang telah lengkap kemudian dimodelkan dan menghasilkan keluaran berupa dua buah model hidrodinamika. Bagian hidrodinamika pertama digunakan untuk melihat kondisi hidrodinamika di perairan Cilacap atara lain berupa arah dan kecepatan arus (U dan V) serta perubahan tinggi muka air laut (surface elevation) terhadap Mean Sea Level (MSL). Bagian hidrodinamika kedua memiliki keluaran berupa debit perairan/fluks dalam arah u dan v serta perubahan kedalaman perairan terhadap waktu (water level). Keluaran hidrodinamika bagian kedua tersebut bersamasama dengan data karakteristik minyak digunakan kembali sebagai masukkan untuk menjalankan modul Spill Analysis berikutnya. Keluaran yang dihasilkan dari pemodelan modul Spill Analysis tersebut selanjutnya menjadi hasil akhir dari seluruh proses pemodelan. Diagram alir dari seluruh proses pemodelan disajikan pada Gambar 5. Kondisi pemodelan yang dilakukan berupa pemodelan pola sebaran tumpahan minyak dengan pengaruh angin (timur dan barat) dan pasang surut setempat. Sedangkan kondisi pemodelan yang diamati yaitu pada saat muka air
22 laut berada pada posisi tertinggi (flood tide), posisi terendah (ebb tide), menjelang pasang dan menjelang surut pada kondisi pertengahan (Mean Sea Level).
Gambar 5. Diagram Alir Pemodelan Sebaran Tumpahan Minyak di Perairan Cilacap dengan Menggunakan DHI Software Mike 21
3.4.1. Lokasi Pemodelan Dalam memutuskan area yang tercakup dalam model, harus pula dipertimbangkan lingkup area, posisi dan tipe dari batas model hidrodinamika yang akan digunakan. Model sebaran tumpahan minyak dibangun dengan skenario di lokasi yang memungkinkan terdapat sumber buangan atau tumpahan minyak masuk ke dalam perairan Cilacap. Desain domain pemodelan berbentuk empat persegi panjang dengan posisi geografis terletak pada 7°46’23” LS 7°41’20” LS dan 108°59’01” BT - 109°03’51” BT ditunjukkan pada Gambar 6.
23 Daerah perairan yang dimodelkan meliputi aliran Kali Donan, Muara Sungai Serayu, alur pelayaran Pelabuhan Tanjung Intan, dan Teluk Penyu. Dalam domain ini digunakan proyeksi WGS 1984 UTM Zone 49S. Domain dibagi ke dalam grid 8850 x 9350 sel dengan lebar ∆x = ∆y = 10 meter.
Gambar 6. Domain Dasar Pemodelan Tumpahan Minyak di Perairan Cilacap dengan Menggunakan Program Mike21
3.4.2. Syarat Batas Syarat batas area pemodelan ditentukan oleh variasi tinggi muka laut yang terdiri dari dua bagian yaitu, syarat batas tertutup dan syarat batas terbuka.
3.4.2.1. Syarat Batas Tertutup Syarat batas tertutup pada area model yaitu berupa garis pantai dimana massa air tidak memungkinkan untuk melewatinya. Berikut ini merupakan lokasi dari syarat batas tertutup pada area model :
24 a) Bagian utara : garis pantai pesisir Cilacap dan Daerah Aliran Sungai (DAS) Donan. b) Bagian selatan : garis pantai pesisir Pulau Nusakambangan. c) Bagian barat : Daerah Aliran Sungai (DAS) Donan.
3.4.2.2. Syarat Batas Terbuka Syarat batas terbuka adalah batas daerah pada model yang berbatasan dengan laut terbuka. Pada area model ini, syarat batas terbuka yaitu antara lain: a) Bagian selatan : garis lurus yang ditarik sejajar dengan Pulau Nusakambangan b) Bagian barat : garis lurus yang memotong aliran Sungai Serayu c) Bagian utara : garis lurus yang memotong aliran Kali Donan d) Bagian timur : garis lurus yang memotong perairan Teluk Penyu
3.4.3. Waktu Pemodelan Waktu pemodelan hidrodinamika terdiri dari dua musim, yaitu musim barat dan musim timur. Pemodelan hidrodinamika pada musim barat dimodelkan pada bulan Februari 2007, sedangkan pemodelan pada musim timur dimodelkan pada bulan Agustus 2007. Waktu pemodelan untuk musim barat yaitu tanggal 1 Februari 2007 hingga 28 Februari 2007. Sedangkan waktu pemodelan untuk musim timur yaitu tanggal 1 Agustus 2007 hingga 28 Agustus 2007.
3.4.4. Skenario Tumpahan Minyak Dalam pemodelan ini terdapat beberapa skenario sumber tumpahan minyak yang berpotensi mencemari perairan Cilacap. Minyak yang akan dimodelkan tumpah dan mencemari perairan Cilacap antara lain avtur, solar
25 (diesel), minyak mentah (crude oil) dan aspal. Sumber tumpahan minyak diskenariokan mengeluarkan minyak dalam jenis, jumlah flux, dan waktu tertentu. Skenario yang disajikan dalam Tabel 2 telah disesuaikan dengan kondisi tumpahan yang memungkinkan terjadi berdasarkan dari data perkapalan setempat. Lokasi terjadinya tumpahan masing-masing minyak ditampilkan pada Gambar 7 di bawah ini .
Gambar 7. Lokasi Skenario Sumber Tumpahan Minyak di Domain Perairan Cilacap
Tabel 2. Informasi Lokasi, Jumlah Tumpahan dan Waktu Pengeluaran Skenario Model Tumpahan Minyak di Perairan Cilacap
Sumber
Bujur (BT)
Lintang (LS)
Lokasi
Potensi
Jenis Minyak
1 2 3 4 5 6 7 8 9 10 11
108°59'24" 108°59'10" 108°59'16" 108°59'24" 108°59'10" 108°59'16" 108°59'05" 108°59'12" 108°59'06" 108°59'16" 108°59'05"
07°46'17" 07°46'15" 07°46'19" 07°46'17" 07°46'15" 07°46'19" 07°46'09" 07°46'15" 07°46'07" 07°46'07" 07°45'59"
Teluk Penyu Jetty Area 70 Jetty CIB Teluk Penyu Jetty Area 70 Jetty CIB Dermaga umum Jetty Area 70 Dermaga umum PPSC Jetty Area 60
Tanker karam Tabrakan tanker Kebocoran loading Tanker karam Tabrakan tanker Kebocoran loading Tabrakan tongkang Limbah dermaga Limbah dermaga Limbah kapal nelayan Tabrakan tanker
Avtur Avtur Avtur Crude Oil Crude Oil Crude Oil Diesel Diesel Diesel Diesel Asphalt
Jumlah Discharge Waktu Tumpahan [m3/s] [menit] 3 [m ] 1800 1.5 25 1800 3 10 300 0.5 10 1800 1.5 25 1800 3 10 300 0.5 10 900 1.5 10 688.84 0.0003 konstan 65.3184 0.000027 konstan 18.6624 0.000008 konstan 600 0.5 20
26
27 3.5. Parameter Pemodelan 3.5.1. Parameter Hidrodinamika Parameter hidrodinamika diawali dengan membuat batimetri pada program Mike 21 sebagai domain model. Perairan Cilacap memiliki nilai batimetri yang bervariasi dengan kisaran kedalaman laut berada di antara nol hingga 25 meter di bawah permukaan laut. Posisi batas selatan dan timur domain berbatasan langsung dengan Samudera Hindia. Kontur batimetri menunjukkan nilai tertinggi pada perairan di sekitar kedua batas tersebut yang ditunjukkan dengan warna ungu. Warna tersebut menunjukkan kisaran kedalaman antara 24 – 25 meter di bawah permukaan laut. Nilai kedalaman semakin mengalami penurunan saat perairan mendekati garis pantai. Perairan pada batas barat maupun utara domain masing-masing berbatasan langsung dengan aliran Sungai Serayu dan Kali Donan. Kedalaman perairan di kedua batas domain tersebut memiliki nilai yang rendah yang ditunjukkan dengan warna kontur hijau dan jingga. Kontur batimetri di perairan Cilacap disajikan pada Gambar 8. Kontur kedalaman laut di perairan Teluk Penyu terlihat semakin merapat saat mendekati garis pantai. Perairan Kali Donan memiliki kontur kedalaman yang rapat dengan kisaran kedalaman bernilai antara 0.88 – 10.56 meter di bawah permukaan laut. Kedalaman perairan di bagian tengah aliran Kali Donan serta di sekitar kolam dermaga/pelabuhan dibuat lebih besar hingga mencapai -11.44 meter. Alur pelayaran Tanjung Intan di sepanjang kanal utama memiliki morfologi dasar laut yang lebih curam dengan kontur kedalaman yang lebih rapat. Kedalaman laut di sepanjang alur pelayaran tersebut berkisar antara 1.76 – 20.73 m di bawah permukaan laut dan terletak memanjang hingga ke perairan Teluk
28 Penyu. Morfologi dasar laut pada alur pelayaran Tanjung intan merupakan morfologi buatan yang dibuat dan dipertahankan untuk kepentingan pelayaran. Terdapat beberapa daerah perairan dangkal di sekitar pantai Cilacap dan Pulau Nusakambangan, yaitu di sepanjang aliran Kali Donan dan di muara Sungai Kaliyasa. Daerah perairan dangkal terdapat pula di sekitar muara Sungai Serayu yang berada di batas barat domain, di mulut alur pelayaran Tanjung Intan dan di sekitar pesisir Pulau Nusakambangan.
Gambar 8. Batimetri Perairan Cilacap Hasil Survey Sounding Dasar Laut (Sumber: JANHIDROS, 2007) Waktu pemodelan hidrodinamika dibagi ke dalam dua musim, yaitu musim timur dan musim barat. Skenario hidrodinamika musim barat dimodelkan pada tanggal 1 Februari 2007 pukul 12:00 AM hingga 28 Februari 2007 pukul 12:00 AM. Skenario hidrodinamika musim timur dimodelkan pada tanggal 1
29 Agustus 2007 pukul 12:00 AM hingga 28 Agustus 2007 pukul 12:00 AM. Langkah waktu masing-masing pemodelan ditentukan sebesar 10 detik disesuaikan dengan syarat kestabilan domain (Courant Number). Courant Number menunjukkan banyaknya grid yang memproses hasil selama pemodelan berjalan dalam satu satuan waktu. Domain area pada skenario pemodelan menggunakan variasi pasang surut air laut pada keempat batas terbuka yaitu, batas utara, batas selatan, batas timur, dan batas barat (Gambar 9).
Gambar 9. Syarat Batas Terbuka pada Domain Model Hidrodinamika di Perairan Cilacap Masing-masing variasi pasang surut pada keempat batas terbuka domain perairan Cilacap yang dimodelkan untuk musim barat disajikan dalam Gambar 10, sedangkan pada musim timur ditampilkan pada Gambar 11
30
Gambar 10. Tinggi Muka Air Laut pada Seluruh Batas Terbuka Domain Perairan Cilacap pada Musim Barat Tahun 2007
Gambar 11. Tinggi Muka Air Laut pada Seluruh Batas Terbuka Domain Perairan Cilacap pada Musim Timur Tahun 2007 Data pasang surut hasil pemodelan bersumber dari data prediksi pasang surut yang didapat dari Jason dan Topex Poseidon. Data tersebut diverifikasi dengan data pasang surut hasil pengukuran insitu yang bersumber dari Bakosurtanal. Masing-masing data pasang surut diukur setiap 15 menit selama 27 hari. Data pasang surut yang diambil pada tanggal 1 - 28 Februari 2007 mewakili kondisi pasang surut pada musim barat, sedangkan data pasang surut yang diambil
31 pada tanggal 1 - 28 Agustus 2007 mewakili kondisi pasang surut pada musim timur. Pengamatan kedua data pasang surut tersebut dilakukan pada posisi 07° 34’ LS - 108° 59’ BT (Gambar 12).
Gambar 12. Lokasi Pengamatan Data Pasang Surut Hasil Pengukuran Lapang dengan Data Masukan Model di Cilacap Tahun 2007 Domain model perairan Cilacap sangat dipengaruhi oleh kondisi pasang surut setempat sehingga perlu ditentukan nilai Drying depth dan Flooding depth. Nilai Drying depth ditentukan dengan memasukan nilai kedalaman minimum yaitu 0.2 dan nilai kedalaman maksimum untuk Flooding depth sebesar 0.3. Nilai masukan parameter tersebut menandakan bahwa perhitungan pemodelan pada masing-masing grid tidak akan dihitung pada kedalaman di atas 0.3 m maupun pada kedalaman di bawah 0.2 m dari Mean Sea Level.
32 Initial surface merupakan nilai awal tinggi muka laut domain saat memulai pemodelan dalam satuan meter. Parameter Initial surface ditentukan dengan memasukkan nilai awal tinggi muka laut yang didapat dari rata-rata tinggi muka laut pada seluruh syarat batas terbuka. Nilai Initial surface pada musim barat ditentukan sebesar 0.12 m dan pada musim timur sebesar 0.4 m. Parameter Source & Sink digunakan untuk menentukan adanya titik sumber masukan dan keluaran air dalam domain. Pada skenario pemodelan hidrodinamika ini, nilai Source & Sink tidak ditentukan karena pada domain tidak diskenariokan terdapat sumber masukan maupun keluaran air. Parameter Eddy Viscosity berhubungan dengan gaya gesek antara molekulmolekul fluida yang bergerak dengan kecepatan berbeda dan menghasilkan gerak turbulen (Alonso dan Finn, 1992). Dalam pemodelan hidrodinamika ini parameter tersebut ditentukan dengan menggunakan formula Smagorinsky. Tipe formula Smagorinsky dihitung berdasarkan kecepatan mengalir fluida dengan nilai konstan sebesar 0.5. Nilai tahanan dasar (bed resistance) pada domain model diberikan dalam parameter Resistance. Nilai tahanan dasar berhubungan dengan kekasaran dasar laut dan gaya gesek antara dasar laut dengan air (DHI, 2007). Konstanta tahanan dasar dalam pemodelan ini menggunakan nilai Manning Number [m1/3/s] dimana pada laut terbuka bernilai 32, sedangkan pada laut dangkal menggunakan nilai tahanan dasar 27 (Gambar 13).
33
Gambar 13. Pola Nilai Tahanan Dasar (Manning Number) dalam Domain Model Perairan Cilacap Data angin yang digunakan untuk masukan model didapat dari IFREMER. Data angin tersebut merupakan data hasil pengamatan satelit yang diukur setiap enam jam. Data angin masukan model kemudian diverifikasi dengan menggunakan data hasil pengukuran insitu yang dilakukan oleh Badan Meteorologi Klimatologi dan Geofisika (BMKG) Cilacap. Data angin insitu merupakan data yang direkam setiap jam selama 28 hari. Pada musim barat, data angin diambil dari tanggal 01 - 28 Februari 2007. Sedangkan pada musim timur, data angin diambil dari tanggal 01 - 28 Agustus 2007. Pengamatan kedua data angin tersebut dilakukan pada posisi 07° 44’ LS - 109° 01’ BT (Gambar 14).
34
Gambar 14. Lokasi Pengamatan Data Angin Hasil Insitu dan Data Angin Masukan Model di Cilacap Tahun 2007 Nilai tekanan yang diberikan oleh angin terhadap permukaan laut diskenariokan bervariasi terhadap ruang dan waktu. Nilai friksi angin pada pemodelan ini diskenariokan bervariasi terhadap kecepatan angin dimana pada saat kecepatan angin bernilai nol, maka besar friksinya 0.0016. Nilai tersebut bervariasi linier dimana pada saat kecepatan angin 16 m/s maka nilai friksinya sebesar 0.0026. Hasil keluaran dari pemodelan hidrodinamika tersebut kemudian dibagi menjadi dua bagian. Bagian hidrodinamika pertama memiliki output berupa surface elevation, U-velocity, dan V-velocity. Sedangkan bagian hidrodinamika kedua memiliki output berupa water level, P flux, Q flux. Contoh hasil laporan pemodelan modul hidrodinamika pada musim barat terdapat pada Lampiran 1.
35
3.5.2. Parameter Spill Analysis Pemodelan pada modul Spill Analysis dibagi menjadi dua bagian, yaitu Basic parameter dan Oil Spill parameter. a.
Basic Parameter Pemodelan basic parameter diawali dengan menentukan Starting
Condition berupa Oil Spill Analysis. Masing-masing minyak yang diasumsikan tumpah memiliki waktu terjadinya tumpahan yang berbeda-beda. Lapisan diesel dan aspal dimodelkan selama 10 hari, lapisan avtur dimodelkan selama 14 hari dan lapisan minyak mentah dimodelkan selama 21 hari baik pada musim barat maupun pada musim timur. Parameter Hydrodynamic Data diisi dengan menggunakan hasil keluaran modul hidrodinamika bagian kedua. Pada output tersebut, arus diberikan dengan variasi terhadap ruang dan waktu. Informasi mengenai lokasi tumpahan (dalam grid), jumlah tumpahan, serta waktu keluaran tumpahan minyak dalam parameter Source disajikan pada Tabel 3 berikut ini. Tabel 3. Informasi Spasial, Jumlah dan Waktu Tumpahan Masing-Masing Jenis Minyak yang di Skenariokan Tumpah di Perairan Cilacap
X 323 93 Diesel 131 447 697 Avtur 449 269 697 Crude 449 269 Asphalt 89
Grid Y 247 450 499 518 201 139 266 201 139 266 748
Discharge [m3/detik] 0.0003 1.5 0.000027 0.000008 1.5 0.5 3 1.5 0.5 3 0.5
Time Release [Time step] Konstan 60 Konstan Konstan 150 60 60 150 60 60 120
36 Koefisien dispersi dalam parameter Dispersion diskenariokan memiliki nilai yang besarnya proporsional terhadap arus. Nilai Longitudinal direction memiliki faktor proporsional sebesar satu, begitu juga dengan nilai Transversal direction. Sementara itu, nilai Vertical direction memiliki faktor proporsional terhadap arus sebesar 0.01. Vertical direction bernilai kecil karena proses dispersi pada lapisan minyak diasumsikan lebih banyak dipengaruhi oleh gerak arus horizontal dibandingkan gerak arus vertikal. Profil arus secara horizontal dipengaruhi oleh gesekan terhadap permukaan dasar laut yang ditentukan dalam parameter Eddy & Logarithmic Velocity Profile. Tipe Velocity profile yang digunakan dalam pemodelan ini yaitu logarithmic velocity profile, sementara nilai Bottom roughness ditentukan sebesar 0.1 m. Parabolic eddy profile disertakan karena berpengaruh dalam penyesuaian proses dispersi vertikal berdasarkan pendekatan gradien. Informasi mengenai suhu dan salinitas air laut di perairan Cilacap pada musim barat dan musim timur ditentukan dalam parameter Water Properties. Suhu dan salinitas air laut tersebut diskenariokan bernilai konstan sepanjang pemodelan. Pada musim barat, salinitas permukaan laut diasumsikan bernilai 33.5 dengan suhu permukaan 290C. Pada musim timur, salinitas permukaan laut diasumsikan bernilai 34 dengan suhu permukaan 250C. Kondisi angin yang diberikan dalam parameter Wind Condition sama dengan data yang diberikan pada modul hidrodinamika. Kedua data tersebut memiliki nilai yang bervariasi terhadap ruang dan waktu. Exceeding Concentration merupakan laju perubahan konsentrasi fraksi minyak. Nilai batas tertinggi pengeluaran konsentrasi minyak pada parameter
37 Exceeding Concentration [%] diskenariokan bernilai 100 mm. Time Exposition merupakan parameter yang digunakan untuk merekam waktu perjalanan lapisan minyak saat mencapai suatu area. Dalam pemodelan ini, time exposition disertakan untuk melihat resident time lapisan minyak dalam domain model. Parameter Line Discharge berfungsi untuk menghitung volume materi yang melewati suatu transek. Dalam pemodelan ini, parameter tersebut tidak digunakan.
b. Oil Spill Parameter Informasi perawanan dan suhu udara di Cilacap pada musim barat dan musim timur disajikan ditentukan dalam parameter Air Properties. Sumber data untuk masukan kedua parameter ini didapat dari BMKG Cilacap. Pada musim barat, nilai Cloudiness diskenariokan konstan sebesar 0.58 dengan temperatur udara 27.565°C. Sementara pada musim timur, nilai Cloudiness diskenariokan sebesar 0.13 dengan temperatur udara 27.527°C Parameter Heat transport digunakan untuk menghitung pertukaran bahang antara minyak dengan air laut dan minyak dengan udara. Nilai konstanta yang digunakan pada perhitungan proses evaporasi dan transfer bahang dalam pemodelan (Tabel 4) menggunakan konstanta yang telah tersedia dalam program Mike 21 (default). Tabel 4. Informasi Nilai Konstanta Transfer Bahang Minyak Heat Balance Albedo Emissivity of Oil Emissivity of Water Emissivity of Air
0.14 0.82 0.95 0.82
Evaporation Evaporation
0.029
38 Proses emulsifikasi pada lapisan minyak ditentukan oleh kehadiran surfactant yaitu kandungan aspal dan wax. Informasi nilai konstanta yang digunakan dalam proses emulsifikasi masing-masing minyak disajikan dalam Tabel 5. Nilai K1 dan K2 berasal dari nilai default yang telah tersedia dalam program Mike 21. Tabel 5. Informasi Nilai Konstanta Emulsifikasi Masing-Masing Minyak Dalam Skenario Model Tumpahan Minyak
Diesel Avtur Crude Asphalt
Max Water Content
Asphaltens Content [wt%]
Wax Content [wt%]
0.85 0.85 0.85 0.85
1 1 0.05 50
2 2 7.04 10
K1 due to water uptake [kg/m3]
5 x 10-7 5 x 10-7 5 x 10-7 5 x 10-7
K2 due to water release [kg/s2]
1.2 x 10-5 1.2 x 10-5 1.2 x 10-5 1.2 x 10-5
Parameter Dissolution & Entrainment disertakan untuk menghitung volume lapisan minyak yang meninggalkan lapisan tersebut karena proses disolusi. Untuk itu, perlu ditentukan nilai Mass transfer coefficient yaitu sebesar 2.36 x 10-6 ks dan nilai Oil in water interfacial tension sebesar 47.2 dyne/cm. Kedua nilai tersebut merupakan nilai default dan telah tersedia dalam program Mike 21. Proses pelapukan pada lapisan minyak ditentukan oleh komponen kimia dari masing-masing minyak. Nilai volume fraksi masing-masing minyak yang diasumsikan tumpah di perairan Cilacap dan digunakan dalam input parameter Oil Properties disajikan dalam Tabel 6 sedangkan laporan hasil pemodelan pada modul Spill Analysis disertakan pada Lampiran 2.
39 Tabel 6. Volume Fraksi Masing-Masing Minyak yang Diasumsikan Tumpah di Perairan Cilacap Diesel Avtur Crude Oil Asphalt [% v/v] [% v/v] [% v/v] [% v/v] Oil Properties C6-C12 (Paraffin) 14.7 5.1 C13-C25 (Paraffin) 32.7 3.78 5 C6-C12 (Cycloparaffin) 34.2 C13-C23 (Cycloparaffin) 43.2 16.2 5 C6-C11 (Aromatic) 9.1 1.8 C12-C18 (Aromatic) 24.1 10 C9-C25 (Naphtean) 42.4 4.1 Residual 69.02 80 Temperatur Reference Temperature 20 20 40 23.5 Viscositas at Reference 6.94 8 4.05 800 Temperature [cs] Oil Temperature Constant [deg C] 25 25 25 25 Selanjutnya, hasil pemodelan modul Spill Analysis dengan menggunakan DHI Software Mike 21 yaitu antara lain: a) Instantaneous oil slick thickness [mm]: Ketebalan total lapisan minyak setelah mengalami proses pelapukan. b) Instantaneous emulsification rate [mm]: Ketebalan lapisan minyak yang mengalami proses emulsifikasi. c) Instant oil evaporation [mm]: Ketebalan lapisan minyak yang mengalami proses evaporasi. d) Instant oil dissolution [mm]: Ketebalan lapisan minyak yang mengalami proses disolusi. e) Instant vertical dispersion [mm]: Ketebalan lapisan minyak yang mengalami proses dispersi vertikal. f) Exceeding concentration [-]: Laju perubahan konsentrasi fraksi minyak.
40 g) Time exposition [second]: Waktu yang dibutuhkan oleh lapisan minyak untuk berada dalam suatu grid.
3.6.
Persamaan Utama Model hirodinamika dalam MIKE 21 HD merupakan sistem model
numerik umum untuk pemodelan permukaan air dan arus. MIKE 21 HD memodelkan arus dua dimensi dalam satu lapis fluida yang diasumsikan homogen secara vertikal. Persamaan berikut merupakan konservasi dari massa dan momentum yang terintegrasi secara vertikal, serta menggambarkan variasi arus dan tinggi muka air:
∂ζ ∂p ∂q ∂d + + = ∂t ∂x ∂y ∂t ∂p ∂t
+
+ gh
(1)
∂ ⎛ p2 ⎞ ∂ ⎛ pq ⎞ ⎜ ⎟+ ⎜ ⎟ ∂x ⎜⎝ h ⎟⎠ ∂y ⎜⎝ h ⎟⎠ gp p 2 + q 2 1 ⎡ ∂ ∂ζ (h τ xx ) + ∂ (h τ xy )⎤⎥ − Ω q + ⎢ 2 2 ρ w ⎣ ∂x ∂x ∂y C .h ⎦
− fVV x +
h ∂
ρ w ∂x
( pa ) = o
(2)
∂q ∂ ⎛ q2 ⎞ ∂ ⎛ pq ⎞ ⎜ ⎟+ + ⎜ ⎟ ∂t ∂ y ⎜⎝ h ⎟⎠ ∂ x ⎝ h ⎠ + gh
gp p 2 + q 2 1 ⎡ ∂ ∂ζ (h τ + ∂y ρ w ⎢⎣ ∂ y C 2 .h 2
− fVVy +
yy
)+
h ∂ (p ) = o ρ w ∂xy a
∂ (h τ xy )⎤⎥ + Ω p ∂x ⎦
Keterangan: h(x,y,t)
= kedalaman perairan [= ζ –d, m]
d(x,y,t)
= kedalaman perairan bervariasi terhadap waktu [m]
ζ(x,y,t)
= elevasi permukaan [m]
(3)
41 p,q(x,y,t)
= densitas flux dalam arah x- dan y- [m3/s/m] = (uh,vh); (u,v) = kecepatan rata2 kedalaman dalam arah x- dan y-
C(x,y)
= Chezy resistance [m½/s]
G
= percepatan gravitasi [m/s2]
f(V)
= faktor gesekan angin
V, Vx, Vy(x,y,t) = kecepatan angin dan komponen dalam arah x- dan y- [m/s] Ω(x,y)
= parameter Coriolis, tergantung latitude [s-1]
Pa(x,y,t)
= tekanan atmosfir [kg/m/s2]
ρw
= densitas air [kg/m3]
x, y
= jarak koordinat [m]
t
= waktu [s]
τxx, τxy, τyy
= komponen shear stress
3.7.
Parameter Oil Spill
3.7.1. Spreading Fay (1969) dalam DHI (2006b) telah membangun teori tiga fase spreading dari lapisan minyak, yaitu : 1. Fase primer, hanya gravitasi (spreading) dan inersia (perlambatan); 2. Fase intermediate, gravitasi dan viskositas (perlambatan); 3. Fase final, tegangan permukaan (spreading) equilibrium dengan viskositas. Mackay et al. (1980) dalam DHI (2006b) kemudian membangun modifikasi formula viskositas-gravitasi dari teori Fay untuk perluasan area lapisan minyak berdasarkan asumsi berikut ini: 1. Minyak dapat dianggap sebagai massa yang homogen;
42 2. Lapisan minyak diasumsikan menyebar sebagai lapisan tipis dan kontinu dalam bentuk melingkar; 3. Diasumsikan tidak ada massa yang hilang dari lapisan. Berdasarkan asumsi tersebut, perubahan area lapisan minyak (Aoil) terhadap waktu dapat digambarkan dalam persamaan berikut:
⎛ dAoil ⎞ 13⎛V ⎜ ⎟ = K a Aoil ⎜⎜ oil ⎝ dt ⎠ ⎝ Aoil
⎞ ⎟⎟ ⎠
43
(4)
Keterangan : Ka
= konstanta [detik-1]
t
= waktu [detik]
Aoil
= πRoil2 [m2]
Volume lapisan minyak didapat dengan menggunakan persamaan berikut: 2 Voil = Roil π ⋅ hs
(5)
Untuk mengetahui ketebalan awal lapisan minyak, dapat diestimasi dengan cara: hs
= 10 cm pada t = 0
Beberapa waktu setelah terjadinya tumpahan di laut, minyak akan berhenti menyebar hingga titik tuang dari partikel-partikel minyak tersebut melewati suhu air laut.
3.7.2.
Evaporation Untuk menghitung tingkat penguapan minyak, diberikan beberapa asumsi
sebagai berikut:
43 1. Tidak terdapat batas difusi dalam lapisan minyak. Hal ini secara umum merupakan asumsi pada temperatur minyak di atas 0°C dan ketebalan lapisan minyak di bawah 5-10 cm. 2. Minyak tercampur sempurna (ideal). 3. Komponen tekanan parsial di udara dapat diabaikan jika dibandingkan dengan tekanan uap. Dengan asumsi tersebut maka tingkat evaporasi dapat digambarkan sebagai berikut:
Nie = kei ⋅ Pi SAT / RT ⋅
Mi
ρi
[
⋅ X i m3 / m2 s
] (6)
Keterangan : Ne
= tingkat penguapan
ke
= koefisien transpor massa
PSAT
= tekanan uap
R
= konstanta gas
T
= suhu
M
= berat molekul
X
= fraksi mol
ρ
= densitas dari fraksi minyak
i
= jenis fraksi minyak ke-i
Perkiraan nilai kei dapat dihitung berdasarkan pada Mackay et al. (1980) dalam DHI (2006b) dan didefinisikan sebagai berikut: 0.045 k ei = k ⋅ Aoil ⋅ S Ci
−2 3
⋅ U w0.78
[m / s]
(7)
44 Keterangan : k
= konstanta (dapat diestimasi)
Aoil
= luas lapisan minyak [m2]
SCi
= konstanta penguapan Schmidts untuk komponen i
Uw
= kecepatan angin [m/detik]
3.7.3.
Vertical Dispersion Fraksi dispersi lapisan minyak di permukaan laut yang masuk ke kolom
perairan per unit waktu dihitung sebagai fraksi yang hilang di permukaan laut dengan kondisi non-wave breaking, dan dapat dihitung dengan persamaan:
D = Da Db
(8)
dimana Da merupakan fraksi dari dispersi minyak di permukaan laut per detik, sedangkan Db merupakan fraksi dari dispersi minyak yang tidak kembali ke lapisan minyak yang dapat digambarkan dengan persamaan:
0.11(1+Uw ) 3600
2
Da =
(9)
dimana : Uw
= kecepatan angin
dan
Db =
1 1+ 50μoilhsγ ow
dimana : μoil
= viskositas minyak [cP]
hs
= ketebalan lapisan minyak [cm]
γow
= tegangan permukaan minyak-air [dyne cm-1]
(10)
45 Tingkat naiknya kembali dispersi butiran minyak ke permukaan dapat dihitung dengan persamaan:
dVoil = Da (1 − Db ) dt
3.7.4.
(11)
Dissolution Dengan asumsi bahwa konsentrasi dari hidrokarbon dapat diabaikan jika
dibandingkan dengan solubility, maka tingkat pelarutan (disolusi) dapat digambarkan sebagai berikut:
dVdsi M = KsiCisat X moli i Aoil dt ρi
(12)
Keterangan :
Cisat
= daya larut fraksi minyak ke-i [mg/kg air]
Xmol
= molar fraksi dari fraksi minyak ke-i
M
= berat molar dari fraksi minyak ke-i [kg/mol]
ρ
= densitas fraksi i [kg/m3]
Aoil
= area tumpahan minyak [m2]
Koefisien transfer massa untuk proses disolusi pada persamaan diatas dapat dihitung sebagai berikut: K Si = 2.36 ⋅10 −6 ei
dimana : 1,4 untuk alkanes ei =
2,2 untuk aromatics 1,8 untuk oilfines
(13)
46 3.7.5.
Emulsification Proses emulsifikasi dapat diketahui dengan mengasumsikan reaksi yang
terjadi sebagai reaksi yang setimbang.
Oil + water
→ ←
Water − in − oil emulsion
Perubahan kandungan air terhadap waktu
dy w dt
dapat dijelaskan sebagai berikut:
dyw = R1 − R2 dt
(14)
Dimana R1 adalah tingkat pengambilan air. Nilai R1 tersebut akan bertambah seiring dengan bertambahnya suhu dan kecepatan angin. Nilai R1 dapat dijabarkan dalam persamaan berikut:
(1 + U w )2
R1 = K1
μ oil
( ywmax − yw ) (15)
Keterangan : Uw
= kecepatan angin
μ oil
= viskositas minyak
y wmax
= kandungan air maksimum (masukan)
yw
= kandungan air dalam minyak
K1
= koefisien yang harus diestimasi (masukan)
R2 adalah tingkat pelepasan air. Nilai R2 berkurang seiring dengan peningkatan kandungan aspal, wax dan surfactan dalam minyak, dan dengan penambahan viskositas minyak. Hal tersebut dapat dijelaskan sebagai berikut:
R2 = K2
1 ⋅ yw As⋅Wax⋅ μoil
Keterangan :
As
= kandungan aspal dalam minyak [wt%]
(16)
1 K 10 ⋅= 5
47 Wax
= kandungan lilin dalam minyak [wt%]
K2
= koefisien yang diestimasi
Selanjutnya konstanta emulsifikasi K1 dan K2 dapat diestimasi dengan: [kg/m3] [kg(wt%)/s] Koefisien K1 dan K2 dapat diestimasi menggunakan data eksperimen dari pengendalian oil spill (Haltenbanken, 1984 in DHI, 2006b).
3.7.6.
Heat transport Tekanan uap dan viskositas lapisan minyak sangat dipengaruhi oleh suhu.
Suhu pada lapisan minyak dapat menjadi lebih hangat daripada udara dan laut sekitar. Oleh karena itu dibangun model untuk menghitung suhu lapisan minyak. Gambar 15 memperlihatkan transfer bahang dari lapisan minyak ke udara dan air laut.
Gambar 15. Transfer Bahang Antara Udara, Lapisan Minyak, dan Air Laut Keterangan: 1 = transfer bahang antara lapisan minyak dan udara, 2 = lapisan minyak menerima dan memancarkan radiasi dari dan ke udara, 3 = lapisan minyak menerima radiasi matahari, 4 = bahang hilang dari lapisan minyak akibat evaporasi,
48 5 = transfer bahang antara lapisan minyak dan air laut, 6 = lapisan minyak menerima dan memancarkan radiasi dari dan ke laut.
1) Transfer bahang antara minyak dan udara Transfer bahang antara lapisan minyak dan atmosfir dapat dijelaskan dalam persamaan berikut:
H Toil − air = Aoil k Hoil − air (Tair − Toil )
(17)
dimana : 0.67
k
oil−air H
⎛S ⎞ = k m ρa C pa ⎜⎜ c ⎟⎟ ⎝ Pr ⎠ air
(18)
Keterangan Sc
= Schmidt’s number
Toil
= suhu minyak [Kelvin]
Tair
= suhu udara [Kelvin]
ρa
= densitas udara [kg/m3]
Cpa
= kapasitas bahang udara [j/kg/°C]
Sedangkan bilangan Prandtl’s dihitung sebagai :
Pr =
C pa ρ a
0.0241(0.18055+ 0.003Tair )
(19)
Dimana kapasitas bahang udara diberikan dalam persamaan (40). Jika tidak terdapat evaporasi, maka kHoil-air dapat dengan mudah dihitung oleh Duffie dan Beckmann (1974) dalam DHI (2006b) sebagai berikut:
k Hoil − air = 5,7 + 3,8U w
(20)
49 2) Transfer bahang antara minyak dan air Transfer bahang antara lapisan minyak dan air dijelaskan sebagai berikut:
H Toil − water = Aoil k Hoil − water (Twater − Toil )
(21)
dimana kHoil-water adalah koefisien transfer bahang yang dihitung oleh Bird et al. (1960) dalam DHI (2006b): −2
k Hoil − water = 0.332 + ρ wC pw Re −0.5 Prw
3
(22)
Kapasitas bahang dari air diberikan dalam persamaan (42). Bilangan Prandtl dari air didasarkan pada persamaan berikut (Duffie dan Beckman, 1974 in DHI, 2006b). ⎛ ⎞ 1 ⎟⎟ Pr w = C pw v w ρ w ⎜⎜ ⎝ 0 . 330 + 0 . 000848 (T w − 273 . 15 ) ⎠
(23)
Re merupakan bilangan Reynolds untuk menghitung koefisien transfer bahang
minyak-air yang dijelaskan dalam persamaan berikut:
Re =
vrel
4 Aoil
ηw
π
(24)
dimana vrel merupakan viskositas kinematik dari lapisan minyak. 3) Solar Radiation
Radiasi matahari yang diterima oleh lapisan minyak tergantung pada beberapa parameter, seperti lokasi terjadinya tumpahan minyak, hari dan waktu penyinaran, perawanan, kandungan air, debu dan ozon di udara. Variasi radiasi matahari dalam satu hari diasumsikan menjadi sinusoidal:
50 ⎧ ⎛ t − t sunrise ⎞ sunrise max 〈t 〈t sunset ⎪ Kt , H o ⋅ sin ⎜⎜ π sunset sunrise ⎟⎟, t t t − ⎠ ⎝ ⎪⎪ H (t ) = ⎨ ⎪ 0, otherwise ⎪ ⎪⎩
(25)
Penyinaran dimulai pada tsunrise dan berakhir pada tsunset. tsunrise = waktu matahari terbit tsunset
= waktu matahari terbenam
tsunset
bisa dihitung dengan menambahkan panjang hari (Td) pada tsunrise
tsunset
= tsunrise + Td
[S]
(26)
Panjang hari dihitung dengan persamaan berikut: Td
= acos(tan ø tan ς)
(27)
dimana ø
= lintang (utara positif)
ς
= deklinasi (posisi angular matahari pada tata surya )
⎛ ⎝
ς ≅ 23.45 ⋅ sin⎜ 360 ⋅
284 + n ⎞ ⎟ 365 ⎠
(28)
Homax merupakan radiasi pada siang hari, yang dihitung oleh Duffie dan
Beckmann (1974) in DHI (2006b).
H omax =
12 ⋅ K t ⎛ ⎛ 360 ⋅ n ⎞ ⎞ ⋅ I sc ⋅ ⎜⎜1 + 0.033 ⋅ cos⎜ ⎟ ⎟⎟ sunrise t −t ⎝ 365 ⎠ ⎠ ⎝ sunset
⋅ (cos(φ ) ⋅ cos(ζ ) ⋅ sin(ω s ) + ω s ⋅ sin(φ ) ⋅ sin(ζ ) ) Keterangan: Isc
= konstanta matahari = 1.353 [W/m]
n
= jumlah hari dalam satu tahun
(29)
51 ωs
= sudut matahari terbit, matahari siang dianggap nol, dan setiap jam sama dengan 15° longitude
dan
Kt =
H radiasi matahari yang mencapai permukaan = Ho radiasi keseluruhan
(30)
Jika langit tidak berawan, Kt ≅ 0.75; Jika tidak Kt akan berkurang seiring pertambahan perawanan. Fraksi besar, a (albedo), dari radiasi matahari yang mencapai tanah akan dipantulkan. Maka masukan bahang bersih dari radiasi matahari dihitung menjadi: (1 − a ) ⋅ H ( t ) [W / m 2 ]
(31)
4) Memancarkan dan menerima radiasi Lapisan minyak akan kehilangan dan menerima bahang dari panjang gelombang radiasi yang dipancarkan. Jumlah bahang yang diterima dan hilang karena radiasi dengan mudah dihitung dengan menggunakan hukum StefanBoltzman. Jumlah bahang bersih yang diterima oleh lapisan minyak dihitung dengan persamaan:
[
rad 4 4 H total = σ (ludara ⋅ Tudara + lair ⋅ Tair4 − 2 ⋅ lmin yak ⋅ Tmin W / m2 yak )
]
Keterangan : σ
= konstanta Boltzman = 5,72·108 [W/(m2K)]
ludara, lair, lminyak
= emisivitas udara, air dan minyak
Tudara,Tair, Tminyak
= temperatur udara, air dan minyak
(32)
52 5) Bahang hilang akibat evaporasi Pendinginan lapisan minyak akibat evaporasi akan menyebabkan lapisan minyak tersebut kehilangan bahang. numberof component
H
vapour
=
∑ i
[
Ni ⋅ ΔH vi W / m2
]
(33)
dimana : ∆Hvi
= bahang penguapan dari komponen i [J/mol]
Keseimbangan bahang dinamis untuk lapisan minyak diberikan dalam persamaan berikut ini:
[
(
dToil 1 4 (1 − a )⋅ H + lairTair4 + lwater Twater = − 2 ⋅ loil Toil4 ξ Cp h dt
)]
+ how (Twater − Toil ) + hoa (Tair − Toil ) − ∑ N i ΔH vi dV ⎛ dV ⎞ + ⎜ water ⋅ ξ w ⋅ C pw + oil ⋅ ξ oil ⋅ C poil ⎟ ⋅ (Twater − Toil ) ⋅ Aoil dt ⎝ dt ⎠
(34)
Keterangan : dWwater = tingkat pengambilan air [m3/s] dt
dVoil dt
= jumlah butiran air terdispersi yang muncul ke permukaan [m3/s]
Cpo
= kapasitas bahang minyak [J/kg°C]
Cpw
= kapasitas bahang air [J/kg°C]
3.7.7.
Sifat Fisik dan Kimia Minyak Sifat dari minyak secara menyeluruh tergantung pada sifat dari unsur-
unsur penyusunnya. Unsur-unsur penyusun minyak tersebut mengalami pelapukan pada tingkat yang berbeda, maka sifat dari lapisan minyak akan berubah terhadap waktu. Sifat dari minyak dijelaskan dengan membagi minyak
53 ke dalam delapan fraksi, yang ditentukan berdasarkan sifat distilasi dan struktur kimianya (alkana atau aromatic). Tabel 7 menjelaskan mengenai fraksi tersebut. Tabel 7. Fraksi Minyak Berdasarkan Struktur Kimia Fraksi Deskripsi 1 C6-C12 (Paraffin) 2 C13-C25 (Paraffin) 3 C6-C12 (Cycloparaffin) 4 C13-C23 (Cycloparaffin) 5 C6-C11 (Aromatic) 6 C12-C18 (Aromatic) 7 C9-C25 (Naphteno-aromatic) 8 Residual (incl. heterocycles) Sumber: DHI, 2006b
Boiling Range 69-230°C 230-405°C 70-230°C 230-405°C 80-240°C 240-400°C 180-400°C >400°C
Viskositas minyak akan bertambah selama proses pelapukan, terutama akibat proses evaporasi dan emulsifikasi. Viskositas sangat tergantung pada temperatur lapisan. Viskositas minyak dapat dihitung dalam tiga langkah. Pertama, menghitung viskositas lapisan minyak tanpa masukan air pada Tref=100°F, menggunakan persamaan Kendall-Monroe: ⎛ 8 ⎞ vToilREF = ⎜ ∑ X i vi1 3 ⎟ ⎝ i =1 ⎠
3
(35)
dimana : Xi
= fraksi model dari i
Kedua, menghitung viskositas lapisan minyak pada temperatur aktual: 0.7
log log
log log
0.7
Keterangan : T
= temperatur [K]
ν
= viskositas kinematik pada suhu T [Cs]
B
= 3.98
log
(36)
54 Ketiga, menghitung viskositas lapisan pada suhu aktual dan kandungan air, menggunakan persamaan Hossain dan Mackay (1980) dalam DHI (2006b). exp
(37)
Keterangan : C4
= kandungan dimensionless dalam minyak [wt%]
Fe
= fraksi minyak yang menguap Efek kombinasi dari emulsifikasi dan evaporasi dapat dihitung dengan
menggunakan persamaan berikut, dimana penjumlahan dari dua efek tersebut dalam bentuk diferensial yaitu .
(38)
Tegangan permukaan dari minyak dapat dengan mudah dihitung dengan: ∑
(39) Kapasitas bahang dari udara, minyak dan air diberikan dalam persamaan
berikut ini dengan suhu dalam Kelvin: 998.73 1684.74 4.3684
. ·
0.133 .
0.00061
.
10
(40) (41) (42)
Untuk minyak tanpa kandungan air , perhitungan titik tuang dapat menggunakan pendekatan berikut ini: ,
(43)
Titik tuang bertambah untuk emulsifikasi dan dihitung dengan persamaan berikut:
Pp oil−water = Pp oil + Pp oil K p 2 y
(44)
4. HASIL DAN PEMBAHASAN
4.1. Verifikasi Hasil Pemodelan 4.1.1. Verifikasi Angin 4.1.1.1. Musim Barat Kecepatan angin masukan model memiliki nilai maksimum pada bulan Februari 2007 sebesar 4.2 meter/detik dengan arah menuju timur laut dan nilai minimum sebesar 0.25 meter/detik dengan arah menuju timur. Rata-rata kecepatan angin masukan model pada bulan tersebut adalah 1.90 meter/detik. Grafik nilai kecepatan serta arah angin masukan model selama bulan Februari 2007 di Cilacap ditunjukkan pada Gambar 16. Kecepatan angin hasil pengukuran lapang di Stasiun Meteorologi Cilacap pada bulan yang sama memiliki nilai maksimum sebesar 6.18 meter/detik dengan arah menuju tenggara dan nilai minimum sebesar nol meter/detik. Nilai rata-rata kecepatan angin hasil pengukuran insitu yaitu sebesar 1.19 meter/detik. Nilai kecepatan serta arah angin insitu selama bulan Februari 2007 di Cilacap ditunjukkan pada Gambar 17.
Gambar 16. Arah [°] dan Kecepatan Angin [m/s] Masukan Model pada Musim Barat 2007 di Perairan Cilacap
55
56
Gambar 17. Arah [°] dan Kecepatan Angin [m/s] Insitu pada Musim Barat 2007 di Perairan Cilacap Pada bulan Februari 2007, angin masukan model maupun hasil pengukuran insitu bertiup dengan kecepatan dan arah yang bervariasi (Gambar 18). Angin masukan model yang bertiup ke arah timur memiliki frekuensi tertinggi yaitu lebih besar dari 20%. Frekuensi angin yang mengarah ke tenggara sekitar 20%, sedangkan frekuensi angin yang bertiup ke arah selatan kurang dari 15%. Kecepatan angin tertinggi, yaitu pada kisaran empat hingga lima meter/detik terutama terjadi saat angin sedang bertiup ke arah timur dan timur laut. Angin hasil pengukuran insitu dominan bertiup ke arah selatan dengan frekuensi bertiup lebih dari 50%. Kecepatan angin yang bertiup ke arah tersebut sebagian besar berada pada kisaran nol sampai satu meter/detik. Sementara itu angin yang bertiup ke arah lainnya pada musim yang sama umumnya memiliki intensitas masing-masing sekitar lebih kurang 10%. Kecepatan angin tertinggi yaitu di atas lima meter/detik terjadi pada saat arah angin sedang bertiup ke tenggara.
57 INPUT MODEL
INSITU
Gambar 18. Windrose Arah dan Kecepatan Angin Masukan Model dan Insitu pada Musim Barat 2007 di Perairan Cilacap Dari penjelasan di atas, dapat dilihat bahwa data angin yang menjadi masukan model memiliki pola yang cukup berbeda dengan data angin hasil pengukuran lapang BMKG. Data angin yang digunakan untuk masukan model memiliki nilai rata-rata yang lebih besar dari data angin insitu. Salah satu faktor penyebab perbedaan nilai kedua data tersebut yaitu terdapat perbedaan metode pengukuran arah dan kecepatan angin antara IFERMER dan BMKG. Selain itu, perbedaan interval pengukuran antara IFREMER dan BMKG akan mempengaruhi data angin yang dihasilkan. Interval pengukuran yang lebih rapat akan memperbesar keakuratan data angin yang dihasilkan. Nilai data masukan model yang lebih besar dapat memperbesar data hasil keluaran. Data angin yang digunakan untuk model memiliki frekuensi arah bertiup hampir merata ke segala arah terutama ke arah timur dan tenggara, sedangkan data
58 angin hasil pengukuran lapang memiliki nilai intensitas yang cenderung dominan ke arah selatan.
4.1.1.2. Musim Timur Pada musim timur, kecepatan angin maksimum hasil masukan model yaitu sebesar 7.75 meter/detik dengan arah bertiup menuju barat laut. Sedangkan kecepatan angin minimum dari hasil masukan model tersebut adalah 3.35 meter/detik dengan arah bertiup juga menuju barat laut. Rata-rata kecepatan angin hasil masukan model pada musim timur ini yaitu sebesar 5.27 meter/detik. Grafik kecepatan serta arah angin hasil model selama bulan Agustus 2007 di Cilacap disajikan pada Gambar 19. Kecepatan angin hasil pengukuran lapang pada musim timur 2007 memiliki nilai maksimum sebesar 7.725 meter/detik dan bertiup ke arah barat. Sedangkan nilai minimum kecepatan angin pada musim tersebut yaitu sebesar nol meter/detik. Nilai rata-rata kecepatan angin insitu pada musim timur adalah 2.80 meter/detik. Grafik kecepatan angin insitu selama bulan Agustus 2007 di Cilacap ditampilkan pada Gambar 20.
Gambar 19. Arah dan Kecepatan Angin Masukan Model pada Musim Timur 2007 di Perairan Cilacap
59
Gambar 20. Arah dan Kecepatan Angin Insitu pada Musim Timur 2007 di Perairan Cilacap Gambar 21 menampilkan grafik kecepatan dan pola arah angin masukan model dan insitu pada musim timur 2007 di Cilacap. Angin masukan model pada musim timur dominan bertiup ke arah barat laut dengan frekuensi bertiup sebesar 70% . Sementara itu, hanya 30% dari arah keseluruhan angin yang bertiup mengarah ke barat. Kecepatan angin terbesar yaitu lebih dari tujuh meter/detik terjadi pada saat angin bertiup menuju timur laut. MASUKAN MODEL
INSITU
Gambar 21. Windrose Arah dan Kecepatan Angin Masukan Model dan Insitu pada Musim Timur 2007 di Perairan Cilacap
60 Angin insitu pada musim timur dominan bertiup ke arah barat dengan frekuensi mencapai 65% dan berkecepatan lebih dari tujuh meter/detik (Gambar 21). Sebanyak 15% dari total keseluruhan angin yang bertiup pada bulan Agustus 2007 menuju ke arah barat laut dan 15% sisanya bertiup ke selatan, sedangkan kurang dari 5% angin yang bertiup ke arah barat daya. Dari penjelasan di atas dapat disimpulkan bahwa pada bulan Agustus 2007 rata-rata kecepatan angin yang digunakan dalam masukan model memiliki nilai yang jauh lebih besar dibandingkan dengan rata-rata kecepatan angin hasil pengukuran lapang. Perbedaan nilai tersebut dipengaruhi oleh perbedaan metode pengukuran antara IFREMER dengan BMKG. Selain itu, interval pengukuran yang dilakukan BMKG lebih rapat jika dibandingkan dengan IFREMER. Hal tersebut akan mempengaruhi keakuratan data, dimana interval pengukuran yang lebih rapat akan semakin mendekati kondisi angin yang sebenarnya. Besarnya kecepatan angin masukan model akan menyebabkan pengaruh angin pada model sebaran lapisan minyak di permukaan laut Cilacap pada musim timur menjadi lebih besar dari kondisi sebenarnya. Namun masing-masing dari data angin tersebut menunjukkan frekuensi arah bertiup yang cenderung sama yaitu dominan menuju arah barat dan barat laut. Perbandingan sebaran data angin insitu dan data angin yang digunakan untuk pemodelan dalam bentuk vektor U dan vektor V pada bulan Agustus 2007 disajikan dalam Gambar 22. Masing-masing data insitu maupun data angin masukan pemodelan dibuat dengan interval waktu yang sama. Persamaan interval waktu dilakukan untuk memudahkan dalam membandingkan kedua data tersebut. Sebagian besar dari kedua data tersebut, baik yang digunakan untuk pemodelan
61 maupun hasil pengukuran lapang memiliki pola sebaran yang hampir serupa. Kedua data tersebut memiliki sebaran yang seragam dan secara dominan berada pada kuadran IV. Vektor angin yang digunakan dalam masukan model memiliki sebaran yang lebih rapat dengan nilai Vektor U berada di bawah -2 radian. Sementara nilai vektor angin yang didapat dari pengukuran lapang memiliki sebaran yang lebih luas dengan nilai Vektor U dimulai dari nol.
Gambar 22. Pola Scattering Data Angin Masukan Model dan Insitu di Perairan Cilacap pada Musim Timur 2007
4.1.2. Verfikasi Pasang Surut 4.1.2.1. Musim barat Gambar 23 menyajikan perubahan tinggi muka air laut masukan model pada bulan Februari 2007 di Cilacap. Perairan Cilacap memiliki pola pasang surut campuran dominasi ganda. Tinggi muka air laut pada saat pasang tertinggi mencapai 0.86 meter di atas permukaan laut. Sedangkan tinggi muka air laut pada saat surut terendah mencapai 0.75 meter di bawah Mean Sea Level (MSL).
62
Gambar 23. Grafik Tinggi Muka Air Laut Hasil Prediksi Model pada Musim Barat 2007 di Perairan Cilacap
4.1.2.2. Musim timur Perubahan tinggi muka air laut hasil masukan model pada bulan Agustus 2007 di perairan Cilacap mewakili kondisi pasang surut pada musim timur dan disajikan pada Gambar 24. Kenaikan muka air laut tertinggi yaitu mencapai 0.87 meter di atas permukaan laut, sedangkan muka air laut terendah pada grafik yaitu 0.75 meter di bawah Mean Sea Level (MSL). Hasil pengukuran lapang menunjukkan bahwa nilai muka air laut tertinggi terjadi saat air laut pasang yaitu 1.03 meter di atas permukaan laut (Gambar 25), sedangkan nilai muka air laut terendah saat perairan mengalami surut yaitu 0.95 di bawah Mean Sea Level.
Gambar 24. Grafik Tinggi Muka Air Laut Hasil Prediksi Model pada Musim Timur 2007 di Perairan Cilacap
63
Gambar 25. Grafik Tinggi Muka Air Laut Hasil Pengukuran Insitu pada Musim Timur 2007 di Perairan Cilacap Dari penjelasan di atas dapat dilihat bahwa data pasang surut hasil masukan model memiliki nilai tinggi muka laut maksimum dan minimum yang lebih rendah dari data pasang surut insitu. Gambar 26 menampilkan perbandingan fluktuasi tinggi muka air laut hasil pemodelan maupun hasil pengukuran lapang di perairan Cilacap dari tanggal 13 Agustus 2007 12:00 AM hingga 16 Agustus 2007 12:00 AM.
Gambar 26. Perbandingan Fluktuasi Tinggi Muka Air Laut Hasil Pemodelan dan Tinggi Muka Air Laut Hasil Pengukuran Insitu di Perairan Cilacap pada Bulan Agustus 2007
64 Dari gambar tersebut terlihat bahwa tidak terdapat beda fase antara pasang surut hasil masukan model dan pasang surut insitu. Namun terdapat perbedaan nilai amplitudo pada kedua data pasang surut tersebut. Umumnya, tinggi muka laut hasil pengukuran insitu memiliki nilai amplitudo yang lebih besar dari data hasil masukan model. Perbedaan amplitudo pada kedua grafik pasang surut tersebut mencapai 0.15 meter. Perbedaan nilai tersebut cukup kecil dan tidak banyak berpengaruh pada sebaran lapisan minyak di Perairan Cilacap.
4.2.
Hasil Pemodelan Hidrodinamika Kondisi hidrodinamika yang diamati setiap musimnya mengacu pada
kondisi pasang surut perairan setempat, antara lain: kondisi pasang tertinggi, surut terendah, menjelang pasang dan menjelang surut saat muka laut berada pada posisi Mean Sea Level (MSL). Penentuan kondisi hirodinamika berdasarkan posisi tinggi muka laut ini bertujuan untuk membandingkan pola pergerakan arus di setiap kondisi tersebut yang akan mempengaruhi pola sebaran lapisan minyak di permukaan laut.
4.2.1. Musim Barat 4.2.1.1. Menjelang Pasang (Muka Laut pada Kondisi MSL) Gambar 27 menampilkan kondisi hidrodinamika perairan Cilacap pada bulan Februari 2007. Kondisi hidrodinamika ditinjau saat air laut di titik P dalam kondisi menjelang pasang dimana muka laut berada dalam posisi Mean Sea Level (MSL). Tinggi muka air laut pada saat MSL seluruhnya berada pada kisaran nol
hingga 0.07 meter di atas permukaan laut. Dalam kondisi tersebut tidak terjadi perbedaan gradien tinggi muka air laut di seluruh perairan dalam domain model.
65 Kondisi angin di titik P (Gambar 27) terlihat mengarah ke tenggara dengan kecepatan angin sebesar 1.3 meter/detik, sedangkan kondisi arus di titik P bergerak menuju timur laut dengan besar kecepatan arus mencapai 0.0124 meter/detik. Arus yang berada pada batas barat domain mengalir di sepanjang kanal utama hingga keluar menuju muara kanal. Di sepanjang Kali Donan terlihat bahwa arus yang mengalir di dalamnya bergerak menuju utara dengan kecepatan yang sangat kecil. Arus di seluruh perairan Teluk Penyu dengan kecepatan rendah bergerak cenderung menuju utara. Sebagian dari arus yang mengalir di perairan Teluk Penyu tersebut mendapat pengaruh dari pembelokan arus yang berasal dari kanal utama.
Gambar 27. Kondisi Hidrodinamika Perairan Cilacap Menjelang Pasang (Muka Laut Berada pada Kondisi MSL) pada Bulan Februari 2007 Arus yang terbentuk di perairan secara dominan masih dipengaruhi oleh kondisi surut pada fase sebelumnya. Hal ini terlihat dari arah arus yang mengarah keluar dari kanal utama menuju Teluk Penyu yang berbatasan dengan Samudera Hindia.
66 4.2.1.2. Pasang Kondisi hidrodinamika hasil pemodelan di perairan Cilacap saat perairan mengalami pasang pada bulan Februari 2007 tersaji dalam Gambar 28. Tinggi muka laut pada saat pasang di titik P yaitu 0.8 meter di atas permukaan laut, sedangkan kondisi tinggi muka air laut secara keseluruhan di perairan Cilacap berkisar antara 0.76 - 0.84 meter di atas permukaan laut. Kondisi angin pada titik P memiliki kecepatan sebesar dua meter/detik dengan arah bertiup menuju tenggara, sedangkan kondisi arus pada titik yang sama memiliki kecepatan sebesar 0.0088 meter/detik dengan arah mengalir menuju barat laut.
Gambar 28. Kondisi Hidrodinamika Perairan Cilacap saat Pasang pada Bulan Februari 2007 Arus pada batas timur domain bergerak langsung menuju ke dalam perairan Cilacap dan mengalir menuju utara, sesuai dengan kondisi perairan yang sedang mengalami pasang. Semakin mendekati pantai, kecepatan arus yang dihasilkan semakin berkurang. Kondisi tersebut disebabkan oleh pengaruh gesekan dasar yang semakin besar akibat perubahan kedalaman di wilayah pantai
67 yang relatif lebih dangkal. Arus pada batas barat domain justru mengarah ke timur atau keluar dari perairan Cilacap. Hal tersebut disebabkan sebagian perairan Cilacap masih dipengaruhi oleh fase surut yang terjadi sebelumnya (beda fase). Arus tersebut kemudian bergabung dengan arus yang berasal dari kanal utama dan bergerak membelok menuju aliran Kali Donan.
4.2.1.3. Menjelang Surut (Muka Laut pada Kondisi MSL) Gambar 29 menampilkan kondisi hidrodinamika perairan Cilacap pada bulan Februari 2007. Kondisi hidrodinamika ditinjau saat air laut di titik P dalam kondisi menjelang surut dimana muka laut berada dalam posisi Mean Sea Level (MSL). Sebagian besar perairan Cilacap memiliki tinggi muka laut antara nol
hingga 0.07 meter di atas permukaan laut, sedangkan disekitar batas barat domain, bagian barat kanal utama, dan disekitar muara Kali Donan memiliki tinggi muka laut berkisar antara -0.08 hingga nol meter di atas permukaan laut. Kondisi angin di titik P memiliki kecepatan sebesar 2.2 meter/detik dengan arah bertiup menuju tenggara. Sedangkan kondisi arus pada titik P memiliki kecepatan sebesar 0.0108 meter/detik serta mengarah ke barat laut. Pola arus pada batas timur domain bergerak masuk menuju Teluk Penyu dengan kecepatan kurang dari 0.25 meter/detik. Arus tersebut bergerak menyusuri pantai Cilacap, kemudian keluar menuju batas timur domain bagian utara. Sebagian arus yang berasal dari batas timur domain berbelok menuju kanal utama kemudian mengalir menuju barat domain dengan kecepatan yang semakin besar. Semakin besarnya kecepatan arus pada daerah kanal tersebut dikarenakan kondisi perairan masih mendapat pengaruh dari fase pasang sebelumnya. Selain
68 itu kondisi geografi perairan yang menyempit dan berbentuk kanal menyebabkan arus bergerak lebih cepat.
Gambar 29. Kondisi Hidrodinamika Perairan Cilacap Menjelang Surut (Muka Laut Berada pada Kondisi MSL) pada Bulan Februari 2007
4.2.1.4. Surut Gambar 30 menampilkan kondisi perairan Cilacap pada saat terjadi surut pada bulan Februari 2007. Seluruh perairan Cilacap memiliki tinggi muka laut yang merata yaitu antara 0.75 – 0.66 meter di bawah permukaan laut. Pada saat surut, angin pada titik P bergerak ke arah tenggara dengan kecepatan bertiup mencapai 2.4 meter/detik, sedangkan arus yang dimodelkan pada titik yang sama bergerak ke arah barat laut dengan kecepatan mengalir mencapai 0.0144 meter/detik. Arus di seluruh perairan Teluk Penyu Cilacap bergerak masuk dari batas timur domain bagian utara, menyusuri pantai Cilacap kemudian berbelok keluar domain melalui batas timur bagian selatan. Pergerakan arus tersebut sesuai dengan fase pasang (flood tide), yaitu bergerak keluar dari domain perairan.
69 Sebagian kecil dari arus tersebut bergerak membelok ke arah kanal utama dengan kecepatan yang semakin kecil. Arus pada kanal utama cenderung mengarah ke barat dan bertemu dengan arus yang berasal dari Kali Donan kemudian keluar dari domain melewati batas barat domain. Dari pola arus yang terbentuk di sepanjang kanal utama, terlihat bahwa sebagian arus bergerak masuk ke alur pelayaran pada saat kondisi laut mengalami surut dikarenakan perairan di kanal utama masih dipengaruhi oleh fase pasang yang terjadi sebelumnya. Arus tersebut kemudian mengalami transisi menuju kondisi surut jika dilihat dari pola arus balik di batas utara domain dan di aliran Kali Donan yang mengarah keluar menuju muara Kali Donan.
Gambar 30. Kondisi Hidrodinamika Perairan Cilacap saat Surut pada Bulan Februari 2007
4.2.2. Musim Timur 4.2.2.1. Menjelang Surut (Muka Laut pada Kondisi MSL) Kondisi hidrodinamika perairan Cilacap menjelang surut dalam posisi muka laut berada pada Mean Sea Level di bulan Agustus 2007 disajikan dalam Gambar 31. Tinggi muka air laut di seluruh perairan Cilacap berkisar antara 0.08
70 hingga nol meter di bawah permukaan laut. Arah angin di titik P bertiup menuju barat dengan kecepatan sebesar 5.9 meter/detik, sedangkan kondisi arus pada titik yang sama memiliki kecepatan sebesar 0.002 meter/detik menuju barat daya. Arus yang masuk dari batas timur domain bagian selatan bergerak menuju utara menyusuri pantai Teluk Penyu dan memutar keluar di batas timur domain bagian utara. Arus yang berasal dari batas timur domain bagian selatan sebagian mengalami pembelokkan menuju kanal utama. Kondisi kanal yang menyempit menyebabkan arus yang mengalir di sepanjang kanal membesar dan bergerak menuju batas barat domain. Sebagian kecil dari arus tersebut membelok ke perairan Kali Donan menuju utara. Arus di sepanjang kali Donan kembali membesar karena aliran sungai yang semakin menyempit di sekitar dermaga tanker.
Gambar 31. Kondisi Hidrodinamika Perairan Cilacap Menjelang Surut (Muka Laut Berada pada Kondisi MSL) pada Bulan Agustus 2007 Keseluruhan pola arus saat menjelang surut terlihat masih mengarah memasuki perairan Cilacap. Kondisi perairan tersebut memiliki arah yang
71 berkebalikan dikarenakan masih dipengaruhi oleh fase pasang yang terjadi sebelumnya.
4.2.2.2. Surut Gambar 32 menyajikan kondisi hidrodinamika di perairan Cilacap pada bulan Agustus 2007 saat muka laut berada pada kondisi surut. Seluruh perairan Cilacap memiliki tinggi muka air laut antara 0.75 – 0.66 meter di bawah permukaan laut, terkecuali pada perairan di sekitar mulut kanal hingga batas timur domain bagian selatan yang memiliki tinggi muka laut lebih rendah yaitu antara 0.83 – 0.75 meter di bawah permukaan laut. Kondisi angin di titik P mengarah ke barat dengan kecepatan bertiup mencapai 5.4 meter/detik, sedangkan kecepatan arus di titik yang sama berkisar antara 0.016 meter/detik dengan arah mengalir menuju ke barat daya.
Gambar 32. Kondisi Hidrodinamika Perairan Cilacap saat Surut pada Bulan Agustus 2007 Keseluruhan pola arus di perairan Cilacap sesuai dengan fase surut yang sedang terjadi dimana sebagian besar arus mengarah keluar dari perairan Cilacap.
72 Arus di sekitar Teluk Penyu bergerak masuk dari batas timur domain bagian utara setelah bergerak menyusuri pantai. Arus tersebut kemudian memutar keluar melewati batas timur domain bagian selatan. Sebagian kecil dari arus tersebut bergerak membelok dan menyusuri kanal utama menuju batas barat domain dengan kecepatan rendah. Kecepatan arus pada aliran Kali Donan juga rendah dan cenderung bergerak ke selatan menuju muara Kali Donan.
4.2.2.3. Menjelang Pasang (Muka Laut pada Kondisi MSL) Gambar 33 menampilkan kondisi hidrodinamika perairan Cilacap menjelang pasang pada bulan Agustus 2007 saat muka laut berada dalam posisi MSL. Dari gambar tersebut terlihat bahwa kondisi perairan seluruhnya memiliki tinggi muka laut yang seragam yaitu berada pada kisaran nol hingga 0.07 meter di atas permukaan laut.
Gambar 33. Kondisi Hidrodinamika Perairan Cilacap Menjelang Pasang (Muka Laut Berada pada Kondisi MSL) pada Bulan Agustus 2007 Kondisi angin di titik P mengarah ke barat dengan kecepatan bertiup mencapai 4.1 meter/detik, sedangkan kondisi arus di titik yang sama mengarah ke barat laut dengan kecepatan mengalir sebesar 0.008 meter/detik.
73 Kondisi arus di perairan Cilacap saat menjelang pasang masih dipengaruhi oleh kondisi arus saat terjadi fase surut sebelumnya. Hal ini terlihat dari pola arus perairan yang masih bergerak keluar dari perairan Cilacap. Di sekitar perairan Teluk Penyu, arus bergerak masuk dari batas timur domain bagian utara kemudian keluar melalui batas timur domain bagian selatan. Sebagian kecil arus di perairan Teluk Penyu bergerak membelok ke dalam kanal utama. Perairan kanal utama juga menerima arus yang berasal dari Kali Donan dan Sungai Serayu. Pertemuan dua arus dengan arah yang berlawanan menyebabkan arah arus pada kanal menjadi tidak beraturan.
4.2.2.4. Pasang Kondisi hidrodinamika perairan Cilacap pada bulan Agustus 2007 saat perairan sedang mengalami pasang diperlihatkan pada Gambar 34. Tinggi muka laut saat terjadi pasang di titik P mencapai 0.82 meter di atas permukaan laut, sedangkan tinggi muka laut di seluruh perairan Cilacap berkisar antara 0.76 - 0.84 meter di atas permukaan laut. Kondisi angin di titik P memiliki arah menuju barat dengan kecepatan bertiup mencapai lima meter/detik. Kecepatan arus di titik P mencapai 0.008 meter/detik dengan arah mengalir menuju barat laut. Kondisi keseluruhan arus saat terjadi pasang umumnya mengarah ke dalam perairan Cilacap. Hal ini membuktikan bahwa kondisi perairan tidak lagi mendapat pengaruh dari fase surut yang terjadi sebelumnya. Arus di sekitar perairan Teluk Penyu masuk dari batas timur domain bagian selatan kemudian bergerak menyusuri pantai Cilacap. Arus tersebut juga membelok ke kanal utama menuju batas barat domain serta berbelok ke Kali Donan. Arus yang dihasilkan di
74 sepanjang kanal utama semakin membesar seiring menyempitnya aliran sungai tersebut.
Gambar 34. Kondisi Hidrodinamika Perairan Cilacap saat Pasang pada Bulan Agustus 2007 Dari penjelasan di atas dapat diketahui bahwa arus yang mengalir di domain perairan Cilacap sangat dipengaruhi oleh pasang surut perairan setempat. Berdasarkan pengamatan pada Titik P, arus yang dihasilkan oleh model saat menjelang pasang maupun saat menjelang surut pada kedua musim masih mendapat pengaruh dari fase sebelumnya yaitu fase surut maupun fase pasang. Hal tersebut menyebabkan arus yang dihasilkan saat kondisi perairan menjelang pasang ataupun menjelang surut saat muka laut berada pada kondisi MSL memiliki arah mengalir yang berkebalikan (beda fase). Kecepatan arus di Titik P pada saat surut relatif memiliki nilai yang lebih besar jika dibandingkan dengan kecepatan arus pada saat pasang. Hal tersebut disebabkan karena arus yang melewati Titik P berasal dari dalam kanal yang lebih
75 sempit. Semakin sempit luas penampang zat cair, maka kecepatan mengalirnya akan semakin besar. Kondisi arus, khususnya arus permukaan di perairan Cilacap juga mendapat pengaruh dari angin yang bertiup di atasnya. Pada saat surut, kecepatan arus pada musim barat relatif lebih besar dibandingkan pada musim timur. Hal tersebut disebabkan pada musim barat, arus saat surut searah dengan arah bertiup angin sehingga resultan keduanya semakin menguatkan. Pada musim timur, arus saat surut dan angin memiliki arah yang berkebalikan sehingga resultan keduanya akan saling melemahkan. Kondisi serupa terjadi pada saat pasang, dimana kecepatan arus pada musim timur relatif lebih besar dibandingkan pada musim barat.
4.3.
Hasil Pemodelan Pola Sebaran Total Minyak Jenis minyak yang dimodelkan dalam skenario model tumpahan minyak di
Peraran Cilacap, Jawa Tengah antara lain diesel, avtur, minyak mentah, dan aspal. Dalam sub bab Hasil Pemodelan Pola Sebaran Total Minyak ini, hanya akan ditampilkan salah satu hasil pemodelan tumpahan minyak yang berasal dari jenis avtur dimana seluruh sumber tumpahannya berasal dari kapal tanker. Sementara pembahasan hasil pemodelan tumpahan minyak lainnya secara keseluruhan akan dibahas pada sub bab Pembahasan Pola Sebaran Total Minyak. Visual hasil pemodelan tumpahan minyak yang disertakan dalam penulisan ini oleh penulis hanya dapat ditampilkan dalam ukuran minimalis. Untuk melihat hasil pemodelan tumpahan minyak tersebut secara utuh dan jelas, dapat dilihat dalam DVD Hasil Pemodelan Tumpahan Minyak (terlampir).
76 44.3.1. Musim Barat 4 4.3.1.1. Kon ndisi Awal Gam mbar 35 menaampilkan koondisi awal teerjadinya tum mpahan minnyak jenis a avtur di peraairan Cilacapp pada bulann Februari 20007. Terdappat tiga sumbber m masukan minyak ke perm mukaan lautt. Sumber tuumpahan perrtama ditand dai dengan l lingkaran meerah yang diiasumsikan m masuk ke linngkungan lauut disebabkaan oleh k kecelakaan k kapal tanker pengangkutt avtur. Jum mlah minyak yang diskennariokan t tumpah men ncapai 1800 m3dengan waktu w keluaraan selama 100 menit.
G Gambar 35. Pola Sebaraan Total Lapisan Avtur saat s Kondisi Awal di Perrairan Cilacap padda Bulan Febbruari 2007 Tum mpahan minyak yang diassumsikan terrjadi akibat kebocoran k peengisian m muatan avtuur ke dalam kapal k tanker ditandai denngan lingkarran berwarnaa kuning. V Volume minnyak yang diiskenariokann tumpah berrjumlah berjuumlah 300 m3 juga
77 dengan durasi tumpahan 10 menit. Lingkaran berwarna hijau menandakan sumber tumpahan minyak yang terjadi akibat karamnya kapal tanker yang bermuatan avtur. Jumlah total minyak yang diskenariokan tumpah yaitu 1800 m3 dengan durasi tumpahan 25 menit. Kondisi awal perairan saat tejadi tumpahan avtur yaitu menjelang surut. Pada kondisi tersebut, lapisan minyak belum menyebar jauh dan masih berada di sekitar lokasi titik sumber dengan ketebalan masing-masing melebihi 144 mm.
4.3.1.2. Menjelang Pasang (Muka Laut pada Kondisi MSL) Pola sebaran tumpahan lapisan avtur di perairan Cilacap menjelang pasang pada bulan Februari 2007 disajikan dalam Gambar 36. Lapisan minyak yang bersumber dari tabrakan tanker telah menyebar menjauhi titik sumber hingga melewati Transek E1 - E2. Ketebalan pada bagian tengah lapisan minyak mencapai lebih dari 144 mm dan semakin tipis saat menjauhi pusat lapisan. Arus yang mengalir pada kanal utama memiliki kecepatan cukup besar, sehingga dapat dengan mudah membawa lapisan minyak keluar dari kanal utama. Lapisan minyak kedua yang bersumber dari kebocoran pengisian muatan tanker menyebar mendekati pantai di arah tenggara dikarenakan terpengaruh oleh angin permukaan dan terbawa oleh arus menyusur pantai. Lapisan minyak tersebut memiliki ketebalan antara 54 - 60 mm di bagian pusat lapisan. Lapisan minyak yang berasal dari kapal tanker karam terlihat menyebar ke utara sesuai dengan gerak arus disekitarnya. Ketebalan lapisan minyak tersebut juga mencapai lebih dari 144 mm dibagian tengahnya.
78
G Gambar 36. Pola Sebaraan Total Lapisan Avtur Menjelang M Pasang (Muka Laut pada Konddisi MSL) di Perairan Cillacap pada B Bulan Februaari 2007
4 4.3.1.3. Pasang Padaa Gambar 377 disajikan poola sebaran tumpahan t laapisan avtur di d perairan C Cilacap saatt mengalami pasang padaa bulan Febrruari 2007. L Lapisan min nyak yang b bersumber dari d tabrakann kapal tankeer telah menyyebar ke perrairan Teluk Penyu m menuju ke arah a tenggaraa sesuai denggan arah perrgerakan arus. Pada konndisi p pasang, aruss di mulut kaanal mengalaami transisi dari d kondisi surut yang terjadi t s sebelumnya . Arus yangg berbalik terrsebut mengaarah ke dalaam perairan Cilacap C s sehingga meenyebabkan lapisan minyyak tertahann di mulut kaanal. Arus yang b berasal dari batas timur domain (Traansek T1 - T2) T bergerak mengarah ke k utara d menyebabkan ujungg lapisan minnyak tersebuut menyebar mengikuti arah dan a arus k utara sehiingga mempperluas perm ke mukaan minyak. Keteballan lapisan minyak m
79 t tersebut telaah jauh berkuurang, yaitu antara 114 - 120 mm di bagian tengaahnya d dikarenakan n telah mengaalami prosess pelapukan.
G Gambar 37. Pola Sebaraan Total Lapisan Avtur saat s Pasang di d Perairan Cilacap C pada Bulann Februari 20007 Minyyak yang berrasal dari kebocoran penngisian muattan tanker tellah m menyebar menjauhi m sum mber tumpahan dan menddekati pantaii utara Pulauu N Nusakamban ngan. Lapisan minyak yang y berasal dari kapal taanker karam m menyebar k arah utaraa kemudian condong ke timur laut mendekati ke m Trransek T1 - T2. T lapisan miny P Penyebaran yak ini mem miliki lintasan n yang palinng jauh dikarrenakan a yang masuk arus m disekittar batas timuur domain berkecepatan tinggi. Lap pisan m minyak terseebut juga meengalami perrluasan perm mukaan lapissan dengan ketebalan k l lebih besar dari d 144 mm m.
80 4.3.1.4. Menjelang Surut (Muka Laut pada Kondisi MSL) Pola sebaran tumpahan minyak jenis avtur pada bulan Februari 2007 saat perairan Cilacap menjelang surut disajikan pada Gambar 38. Minyak yang berasal dari tabrakan kapal tanker menyebar dengan arah berbalik mendekati mulut kanal. Pola sebaran ini sesuai dengan pola sebaran arus di perairan sepanjang kanal yang masih dipengaruhi oleh kondisi pasang sebelumnya. Pusat lapisan minyak masih berketebalan lebih dari 144 mm, namun sebagian besar lapisan memiliki ketebalan di bawah 108 mm. Perubahan ketebalan lapisan tersebut disebabkan oleh adanya proses pelapukan. Lapisan minyak yang bersumber dari kebocoran pengisian muatan tanker menyebar sesuai dengan arah pergerakan arus yaitu menuju ke barat. Lapisan tersebut tetap berada di sekitar garis pantai utara Pulau Nusakambangan dan tidak meyebar jauh dikarenakan kecepatan arus menyusuri pantai di lokasi tersebut juga tidak terlalu besar. Ketebalan maupun luas permukaan lapisan minyak tersebut telah jauh berkurang disebabkan adanya proses pelapukan minyak. Minyak yang berasal dari kebocoran kapal tanker karam, tidak lagi terdapat dalam domain. Lapisan minyak tersebut sebelumnya telah keluar dari domain melewati Transek T1 - T2.
81
G Gambar 38. Pola Sebaraan Total Lapisan Avtur Menjelang M Surut (Muka Laut pada Kondisi MSL) di Perairan Cilacap pada Bulan Februari 200 07
4 4.3.1.5. Surrut Gam mbar 39 menyyajikan pola sebaran tum mpahan lapissan avtur di perairan p C Cilacap saatt mengalami surut pada bbulan Februaari 2007. Laapisan minyaak yang b berasal dari tabrakan kappal tanker, menyebar m meenjauhi muluut kanal dan c cenderung b bergerak menndekati pantaai Pulau Nussakambangaan. Lapisan minyak m t tersebut kem mudian terpissah menjadi dua bagian. Lapisan peertama meruppakan b bagian dari pusat p lapisann sebelumnyya dan terlihaat mengalam mi penurunann ketebalan l lapisan. Lap pisan kedua terbentuk daari hasil akum mulasi lapisaan minyak awal a dan m memiliki kettebalan yang g lebih tingggi pada bagiaan tengahnyaa. Lapisan minyak m y yang berasall dari kebocooran pengisian tanker tettap menyebaar di sekitar garis
82 p pantai Pulau u Nusakambaangan. Keteebalan lapisaan minyak teersebut umum mnya b berada di baawah enam milimeter. m
G Gambar 39. Pola Sebaraan Total Lapisan Avtur saat s Surut di Perairan Cillacap pada Bulan Febrruari 2007 Padaa gambar di atas a terlihat bbahwa lapisan minyak yang y menyebbar tepat m melewati titiik monitor hanya h terdapaat pada Titikk Monitor E.. Lapisan minyak m y yang melewati titik monnitor tersebutt berasal darri peristiwa tabrakan kap pal tanker d dalam kannal utama. Lapisan di L avtuur yang terdaapat di dalam m domain moodel u umumnya m menyebar meelewati pingggir transek. Keteebalan minyaak yang meliintasi masingg-masing traansek tidak terlihat t j jelas pada gaambar. Hal ini dikarenaakan lapisan minyak tersebut melintaasi transek m memiliki diaameter perm mukaan yang kecil. Namuun jika dilihhat dalam Viddeo Total O Avtur Baarat (terlamppir), maka daari seluruh trransek tersebbut hanya Trransek T1 Oil
83 – T2 dan Transek E1 – E2 yang dilintasi oleh lapisan minyak. Transek T1 – T2 terletak pada batas timur domain, sedangkan Transek E1 – E2 terletak pada mulut kanal utama. Ketebalan lapisan minyak yang melewati Transek E1 – E2 mencapai lebih dari 144 mm dengan diameter mencapai 200 meter, sedangkan ketebalan lapisan minyak yang melewati Transek T1 – T2 mencapai lebih dari 144 mm dengan diameter mencapai 400 meter. Pergerakan lapisan minyak yang hanya melewati kedua transek tersebut disebabkan oleh pengaruh kondisi musim barat, dimana angin dominan bertiup dari arah barat.
4.3.2. Musim Timur 4.3.2.1. Kondisi awal Gambar 40 menyajikan kondisi awal pola sebaran tumpahan lapisan avtur di perairan Cilacap pada bulan Agustus 2007. Tumpahan minyak pada musim timur diasumsikan memiliki sumber yang sama dengan musim barat. Kondisi perairan saat awal tumpahan yaitu sedang mengalami pasang. Lapisan minyak yang bersumber dari tabrakan kapal tanker telah menyebar jauh ke dalam kanal utama hingga berada di tepi utara kanal dengan ketebalan lapisan lebih dari 14 mm pada bagian tengahnya. Sebaran minyak tersebut bergerak ke arah barat laut sesuai dengan vektor arus dan angin. Vektor arus dan angin pada musim timur sangat berpengaruh terhadap sebaran lapisan minyak dikarenakan kecepatan angin yang cukup besar.
84
G Gambar 40. Pola Sebaraan Total Lapisan Avtur saat s Kondisi Awal di Perrairan Cilacap padda Bulan Aggustus 2007 Lapisan minyak yang berasaal dari kebocoran saat penngisian tankker m menyebar seearah dengann arah angin hingga men ncapai mulutt kanal. Bag gian tengah l lapisan avturr memiliki ketebalan k hinngga 60 - 66 mm sementtara minyak yang m melewati paada Titik Monitor E mem miliki keteballan lebih darri 100 mm. Lapisan m minyak yang g bersumberr dari tanker karam meny yebar menujuu ke arah baarat laut p perairan Teluk Penyu. Ketebalan K laapisan tersebuut mencapaii 144 m.
4 4.3.2.2. Men njelang Surrut (Muka L Laut pada Kondisi K MSL L) Pola sebaran tum mpahan minyyak di perairan Cilacap m menjelang suurut pada b bulan Agusttus 2007 disaajikan pada Gambar G 41. Lapisan miinyak yang bersumber b d tabrakann kapal tankker telah mennyebar menu dari uju ke arah barat b laut mengikuti a arah pergeraakan arus meenyusur panttai di sepanjaang kanal. P Pergerakan lapisan
85 m minyak terseebut kemudiian terhalangg oleh dermaaga yang terlletak di sebeelah barat d daratan Cilaacap, sehinggga minyak teerjebak dan terakumulasi t i. Luas perm mukaan l lapisan miny yak tersebut menjadi berrkurang nam mun ketebalannnya bertam mbah di s seluruh bagiian lapisan hingga h mencaapai lebih daari 144 mm.
G Gambar 41. Pola Sebaraan Total Lapisan Avtur Menjelang M Surut (Muka Laut pada Kondisi MSL) di Perairan Cilacap pada Bulan Agustus 20007 Lapisan minyak yang berasaal dari kebocoran pengisiian muatan tanker t m menyebar jaauh ke dalam m kanal sesuaai dengan peergerakan aruus di sepanjaang kanal d menyebar di bagian selatan daraatan Cilacap. Luas perm dan mukaan lapisaan tersebut s semakin mem mbesar, nam mun ketebalaannya semakkin berkurangg yaitu beradda pada k kisaran antarra 54 - 60 mm. m Lapisann minyak yanng bersumbeer dari kapal tanker k karam telah menyebar ke k utara peraiiran Teluk Penyu. P Luas permukaan lapisan
86 tersebut semakin membesar dan ketebalan lapisannya mencapai lebih dari 144 mm. Sementara di bagian tepi lapisan, ketebalannya hanya mencapai 12 - 18 mm.
4.3.2.3. Surut Gambar 42 menyajikan pola sebaran tumpahan lapisan avtur di perairan Cilacap saat mengalami surut pada bulan Agustus 2007. Arah arus yang bergerak ke barat laut membuat lapisan minyak yang bersumber dari tabrakan kapal tanker tetap terperangkap di sekitar dermaga pelabuhan. Keadaan tersebut yang disertai dengan proses pelapukan membuat lapisan minyak mengalami pengurangan luasan permukaan maupun ketebalan lapisan minyak. Lapisan minyak yang berasal dari kebocoran pengisian muatan tanker menyebar mengikuti arah arus menyusuri pantai menuju ke batas barat domain (Transek B1 - B2). Lapisan minyak tersebut mengalami penurunan luas permukaan serta pengurangan ketebalan lapisan hingga mencapai 12 - 18 mm. Lapisan minyak yang berasal dari kapal tanker karam menyebar mendekati pantai Cilacap. Pada saat terjadi pasang, lapisan minyak bergerak mendekati pantai menuju utara. Namun saat terjadi surut, arus berbalik kembali ke selatan dan diikuti oleh lapisan minyak. Lapisan minyak tetap berada di sekitar pantai akibat dipengaruhi oleh arus menyusur pantai serta angin yang bertiup. Luas permukaan lapisan minyak semakin membesar dengan ketebalan lapisan tetap lebih besar dari 144 mm pada bagian pusatnya.
87
G Gambar 42. Pola Sebaraan Total Lapisan Avtur saat s Surut di Perairan Cillacap pada Bulan Agustus 2007
4 4.3.2.4. Men njelang Passang (Muka Laut pada Kondisi MS SL) Pola sebaran totaal lapisan avvtur pada di perairan p Cilaacap saat meenjelang p pasang padaa bulan Agusstus 2007 dissajikan dalam m Gambar 43. Lapisan minyak m y yang berasall dari tabrakan kapal tannker maupunn yang bersum mber dari keebocoran p pengisian m muatan tankerr sudah tidakk tampak lag gi pada gambbar. Lapisann tersebut t telah mengaalami pelapuk kan serta terrbawa arus keluar k dari doomain modeel. M Meskipun beegitu, masih h terdapat miinyak yang teerakumulasii dan terdam mpar di s sekitar beberrapa dermagga dengan keetebalan lebiih dari 144 m mm. Semeentara itu, laapisan minyaak yang beraasal dari kapal tanker karram m menyebar keembali ke paantai Cilacapp sesuai denggan pergerakkan arus mennyusur p pantai setelaah sebelumnyya keluar daari domain model m melew wati Transek U1 - U2
88 d dengan diam meter mencap pai 200 meteer. Terdapatt lapisan minnyak yang mengarah m m masuk ke daalam mulut dermaga d nelaayan. Keteb balan lapisann avtur meninngkat d dimana ham mpir seluruh permukaann p nya memilikii ketebalan leebih dari 1444 mm.
G Gambar 43. Pola Sebaraan Total Lapisan Avtur Menjelang M Pasang (Muka Laut pada Konddisi MSL) di Perairan Cillacap pada B Bulan Agustuus 2007
4 4.3.2.5. Pasang Gam mbar 44 menyyajikan pola sebaran tum mpahan lapissan avtur di perairan p C Cilacap saatt pasang pada bulan Aguustus 2007. Lapisan L minnyak yang masih b berada dalam m domain model m berasall dari kapal karam k di perrairan Teluk Penyu. A Adanya geraak arus yangg memasuki pperairan Cilaacap saat terj rjadi pasang m menyebabka an lapisan minyak m tersebbut menyebaar memasuki dermaga neelayan. K Ketebalan laapisan minyaak tersebut bberkisar antaara 24 mm hiingga mencaapai lebih d 144 mm dari m.
89
G Gambar 44. Pola Sebaraan Total Lapisan avtur saaat Pasang ddi Perairan Cilacap C pada Bulann Agustus 20007 Padaa gambar di atas a terlihat bbahwa hanyya Titik Monnitor E yang dilintasi d o lapisan avtur dengaan ketebalan lebih dari 100 mm. Polla sebaran laapisan oleh m minyak umuumnya melew wati bagian pinggir p dari masing-massing transek.. Lapisan m minyak yang g berasal darri kebocorann pengisian muatan m tankeer menyebarr melewati T Titik Monito or E dan mellewati Transsek E1 – E2 dengan diam meter mencappai 100 m meter. Lapisan avtur yaang dimodelkkan tumpah pada p musim m timur samaa sekali t tidak menyeebar melewatti Transek T1 T – T2. Hall tersebut terj rjadi karena sebaran l lapisan miny yak di permuukaan laut diipengaruhi oleh o kondisi angin pada musim m t timur yang dominan d berrgerak dari tiimur menujuu barat dan bbarat laut.
90 Kecilnya diameter permukaan minyak saat melintasi transek menyebabkan sebaran lapisan avtur tidak tampak jelas. Transek U1 – U2 dilintasi lapisan minyak dengan posisi lintasan cenderung berada di sekitar Titik U1 yang berbatasan dengan garis pantai Cilacap. Pola sebaran minyak yang demikian disebabkan oleh adanya gerak arus menyusuri pantai serta pengaruh kondisi angin yang bertiup menuju barat laut dengan kecepatan yang cukup besar. Lapisan minyak yang melewati Transek U1 – U2 memiliki ketebalan lebih dari 150 mm dan diameter mencapai 200 meter. Transek W1 – W2 dan Transek B1 – B2 juga dilintasi oleh lapisan minyak yang berasal dari dalam kanal utama. Lapisan minyak yang melewati kedua transek tersebut memiliki luas permukaan yang kecil namun ketebalan lapisannya cukup besar. Lapisan minyak yang melewati Transek B1 – B2 memiliki ketebalan mencapai 120 mm dan diameter mencapai 100 meter, sedangkan lapisan minyak yang melintasi Transek W1 – W2 memiliki ketebalan lebih dari 150 mm.
4.4.
Pembahasan Pola Sebaran Tumpahan Minyak Seluruh lapisan minyak yang diasumsikan tumpah di perairan Cilacap
dalam pemodelan ini, mengalami proses pelapukan seperti evaporasi, disolusi, emulsifikasi dan dispersi vertikal. Total ketebalan minyak dari berbagai proses tersebut selama mengalami pelapukan disebut sebagai total minyak (total oil). Total ketebalan lapisan pada masing-masing jenis minyak memiliki nilai yang bervariasi (Tabel 8).
Pasaang
ASPAL
CRUDE
AVTUR
DIESEL
Tabel 8. Perbandinga P an Pola Sebaaran Total Lapisan L Dieesel, Avtur, Minyak Meentah, dan Aspal A pada Berbagai B K Kondisi Muk ka Laut saat Musim Barrat dan Mussim Timur di Perairan Cilacap Taahun 2007 MUSIIM TIMUR MUS SIM BARA AT R Meenjelang Pa asang Pa asang Meenjelang Su urut Su urut Meenjelang Su urut Surrut Men njelang Pasang
na Skala Warn Diesel
Avtu ur
Crude
Aspaal
91
9 91
92 Total lapisan minyak yang tumpah pada musim barat dan musim timur memiliki pola sebaran yang berbeda. Perbedaan pola sebaran tersebut disebabkan oleh beberapa faktor, antara lain perbedaan kondisi awal pemodelan pada masingmasing musim serta perbedaan kondisi angin pada kedua musim. Pada musim barat, pemodelan tumpahan minyak dimulai saat kondisi perairan sedang mengalami pasang sehingga arus di sekitar perairan bergerak ke luar domain. Kondisi angin pada musim barat yang bertiup ke arah barat dan tenggara menyebabkan resultan antara arus dan angin saling menguatkan sehingga lapisan minyak ikut menyebar jauh menuju Teluk Penyu. Daerah perairan yang rawan terkena pencemaran tumpahan minyak pada terutama di sekitar pesisir utara Pulau Nusakambangan dan perairan Teluk Penyu serta tepi timur aliran Kali Donan. Pada musim barat, tingkat kerawanan pencemaran minyak pada perairan Cilacap bersifat sementara serta high recovery dikarenakan sebagian besar lapisan minyak cenderung menyebar meninggalkan domain menuju Samudera Hindia. Pada musim timur, lapisan minyak dimodelkan tumpah saat perairan sedang mengalami kondisi surut sehingga arus laut di sekitarnya mengarah masuk ke dalam perairan Cilacap. Resultan arus tersebut semakin diperkuat oleh kondisi angin pada musim timur yang bertiup kencang menuju timur dan barat laut sehingga menyebabkan lapisan minyak yang tumpah di permukaan laut tersebar cukup jauh ke dalam perairan Cilacap. Daerah perairan Cilacap yang sangat rawan terhadap pencemaran minyak yaitu meliputi aliran kanal utama, tepi barat aliran Kali Donan, dan daerah sekitar pesisir Pantai Cilacap. Lebih dari itu, tumpahan minyak pada musim timur dapat memasuki daratan melalui aliran sungai Kaliyasa yang mengalir membelah daratan Cilacap. Pola sebaran lapisan
93 minyak yang bergerak menuju ke dalam perairan Cilacap menyebabkan risiko kerawanan pencemaran minyak terhadap perairan Cilacap pada musim timur lebih tinggi dan lebih persistent dibandingkan pada musim barat. Pada keseluruhan proses pelapukan yang terjadi pada semua jenis minyak, lapisan aspal memiliki total ketebalan lapisan tertinggi yaitu lebih dari 192 mm. Lapisan minyak avtur dan minyak mentah memiliki ketebalan mencapai lebih dari 144 mm, sementara diesel memiliki ketebalan lapisan terkecil yaitu lebih dari 0.0288 mm. Lapisan minyak tersebut memiliki ketebalan yang berbeda-beda dikarenakan mengalami proses pelapukan dengan tingkat berbeda pula.
4.5.
Hasil Pemodelan Proses Pelapukan Tumpahan Minyak di Laut Minyak yang masuk ke lingkungan perairan laut akan membentuk lapisan
tipis di atas permukaan laut. Lapisan tersebut kemudian akan mengalami proses pelapukan yang disebabkan oleh pengaruh cuaca. Proses pelapukan yang terjadi pada lapisan minyak tergantung dari masing-masing jenis minyak yang tumpah. Salah satu contoh hasil pemodelan yang akan dijelaskan dalam sub bab ini yaitu proses pelapukan minyak dari jenis avtur. Tingkat pelapukan minyak tersebut dibahas per satu jam, 12 jam, 24 jam dan menjelang akhir waktu pemodelan. Pembahasan mengenai proses pelapukan masing-masing jenis minyak lainnya akan dibahas pada sub bab Pembahasan Proses Pelapukan Tumpahan Minyak di Laut.
4.5.1. Musim Barat 4.5.1.1. Jam ke-1 Minyak jenis avtur yang tumpah di perairan Cilacap diasumsikan tumpah dari tiga titik sumber utama. Lapisan minyak pertama bersumber dari tabrakan
94 k kapal tankerr di sekitar dermaga tankker. Lapisann minyak keddua bersumb ber dari k kebocoran p pengisian muuatan tanker,, sedangkan lapisan ketigga berasal daari kapal t tanker karam m di alur pelaayaran. Kettiga sumber tumpahan t avvtur dimodellkan mulai t tumpah padaa waktu yangg bersamaann namun mem miliki total w waktu tumpaah yang b berbeda. Gaambar 45 meerupakan polla sebaran seerta proses pelapukan p lappisan avtur y yang tumpah h di perairann Cilacap settelah satu jam m pada bulann Februari 2007.
G Gambar 45. Sebaran serta Proses Peelapukan Lappisan Avtur JJam ke-1 di Perairan Cilacap padda Bulan Febbruari 2007 Padaa gambar terssebut terlihaat bahwa satuu jam setelahh terjadinya tumpahan, l lapisan miny yak tersebut masih menyyebar di sekiitar titik sum mber. Luas permukaan p ing tumpahaan relatif keccil jika diban m masing-mas ndingkan denngan volumee avtur y yang tumpah h ke laut. Kondisi K lapisaan minyak yang belum m menyebar serrta kondisi
95 luas permukaan yang sempit membuat penyebaran lapisan minyak tersebut masih mudah untuk ditangani.
4.5.1.2. Jam ke-12 Sebaran serta proses pelapukan lapisan avtur yang tumpah di perairan Cilacap setelah jam ke-12 pada bulan Februari 2007 dapat dilihat pada Gambar 46. Dari gambar tersebut terlihat bahwa hanya terdapat dua lapisan avtur yang berada dalam domain perairan Cilacap. Lapisan minyak pertama yang bersumber dari tabrakan kapal tanker menyebar di mulut kanal dan memiliki luas permukaan yang cukup besar dengan panjang lapisan mencapai 1.2 km dan lebar lapisan 0.4 km. Lapisan avtur kedua yang bersumber dari kebocoran pengisian muatan menyebar di dekat pantai utara Pulau Nusakambangan dan membentuk garis dengan panjang lapisan mencapai 0.3 km. Lapisan minyak pertama menyebar keluar dari kanal utama dan telah melakukan perjalanan sepanjang 2.2 km, sedangkan lapisan minyak kedua menyebar menuju pantai Pulau Nusakambangan sejauh 1.2 km dari titik sumbernya. Lapisan avtur ketiga yang bersumber dari kapal tanker karam menyebar dan bergerak menuju utara sejauh 2.8 km dengan lebar lintasan mencapai 0.3 km. Lapisan tersebut kemudian keluar meninggalkan perairan Cilacap melewati batas timur domain pada jam ke-12 setelah terjadinya tumpahan. Proses evaporasi lapisan minyak yang bersumber dari tabrakan kapal tanker memiliki tingkat evaporasi antara nol hingga 0.0012 mm. Sedangkan pada lapisan minyak yang bersumber dari kebocoran pengisian muatan memiliki tingkat evaporasi antara 0.0050 - 0.0056 mm. Proses disolusi lapisan minyak
96 y yang pertam ma memiliki tingkat t disollusi antara 5 x 10-8 – 1.2 x 10-7 mm, sedangkan p pada lapisan n kedua mem miliki tingkatt disolusi anttara 3.8 x 100-7 – 4.2 x 100-7 mm.
G Gambar 46. Sebaran serta Proses Peelapukan Lappisan Avtur JJam ke-12 di d Perairan Cilacap padda Bulan Febbruari 2007 Ting gkat dispersi vertikal padda lapisan miinyak yang pertama p sebaagian besar b berada pada kisaran 2.3 x 10-7 – 7.1 x 10-7 mm, sedangkan ppada lapisan minyak k kedua memiiliki ketebalaan di atas 1.66 x 10-6 mm.. Ketebalan emulsifikasi pada l lapisan avturr pertama yaaitu mencapaai 12.50 mm m, sedangkann pada lapisaan avtur y yang kedua memiliki m ketebalan antaara 56.25 - 622.50 mm. Exceeedance frequency meruppakan laju perubahan koonsentrasi fraaksi t tumpahan minyak m saat berada b di perrmukaan lautt. Lapisan avtur a yang peertama, m memiliki nillai exceedance frequencyy tertinggi paada lapisan yang y terletak k di sekitar s sumber tumppahan yaitu antara 3.7 - 4.1 %. Lajuu perubahan konsentrasi fraksi
97 tumpahan minyak tersebut semakin berkurang saat menjauhi sumber tumpahan. Exceedance frequency pada pusat lapisan minyak memiliki nilai yang lebih tinggi
yaitu di atas 4.5 %, sedangkan pada perairan yang baru saja dilintasi oleh lapisan minyak, laju perubahan konsentrasi fraksi tumpahan minyak tersebut memiliki nilai yang rendah, yaitu antara nol hingga 0.4 %.
4.5.1.3. Jam ke-24 Sebaran serta proses pelapukan lapisan avtur yang tumpah di perairan Cilacap setelah jam ke-24 pada bulan Februari 2007 disajikan pada Gambar 47. Pada gambar tersebut terlihat bahwa kedua lapisan minyak yang terdapat di perairan Cilacap telah terakumulasi menjadi satu lapisan. Lapisan minyak pertama yang berasal dari tabrakan kapal tanker bergerak menuju pantai pulau Nusakambangan dan menyatu dengan lapisan minyak kedua. Lapisan minyak kedua yang bersumber dari kebocoran pengisian muatan sebelumnya telah lebih dahulu menyebar hingga ke pantai Pulau Nusakambangan. Panjang pantai Nusakambangan yang telah dilewati oleh lapisan avtur tersebut mencapai dua km. Di daerah pesisir Pulau Nusakambangan yang membentuk cekungan, lapisan minyak terkonsetrasi dan memiliki ketebalan evaporasi antara 0.0005 mm hingga mencapai 0.0065 mm. Sementara itu, proses disolusi pada lapisan minyak tersebut yaitu antara 2 x 10-7 mm hingga lebih dari 4.6 x 10-7 mm. Proses dispersi vertikal pada lapisan minyak di sepanjang pantai Pulau Nusakambangan memiliki ketebalan antara 7.5 x 10-8 mm hingga 3.9 x 10-7 mm. Sementara proses emulsifikasi lapisan minyak tersebut bervariasi antara 12.5068.75 mm.
98 Laju perubahan konsentrasi k ffraksi minyaak cenderungg semakin keecil pada s seluruh areaa yang telah dilewati d olehh lapisan minnyak. Laju pperubahan konsentrasi k f fraksi tetap tinggi t pada lapisan l minyyak yang meemiliki ketebbalan tinggi.
G Gambar 47. Sebaran serta Proses Peelapukan Lappisan Avtur JJam ke-24 di d Perairan Cilacap padda Bulan Febbruari 2007
4 4.5.1.4. Jam m ke-96 Sebaaran serta prooses pelapukkan lapisan avtur a yang tuumpah di perrairan C Cilacap setelah jam ke-996 pada bulaan Februari 2007 2 disajikaan dalam Gaambar 48. D Dalam wakttu empat harii sejak terjaddinya tumpaahan, lapisann minyak tidaak lagi m menyebar dii dalam dom main perairann Cilacap meelainkan telahh hilang akib bat proses p pelapukan seerta terbawaa arus menujuu ke luar dom main perairaan Cilacap. Lapisan m minyak telah h menyebar terutama di sekitar muluut kanal hinggga mencapaai daerah di s sekitar batass timur domaain. Jarak linntasan minyyak yang mem mbujur dari barat ke
99 t timur telah mencapai m em mpat kilometter, sedangkaan jarak darii utara ke sellatan juga m mencapai em mpat kilometter.
G Gambar 48. Sebaran serta Proses Peelapukan Lappisan Avtur JJam ke-96 di d Perairan Cilacap padda Bulan Febbruari 2007 Dalaam Gambar 48 4 tersebut tiidak lagi terllihat adanyaa lapisan minnyak yang m mengalami proses p evapo orasi, disolussi, dipersi veertikal, mauppun emulsifiikasi. Laju p perubahan k konsentrasi fraksi fr lapisann minyak di permukaan p l setelah jam laut j ke-96 b besar beradaa pada kisaraan nol hingga 0.4%. Maasih terdapat nilai laju peerubahan k konsentrasi fraksi yang agak a tinggi di d sekitar panntai Pulau N Nusakambanggan, yaitu a antara 1.6 - 2.0%. 2 Lapissan minyak ttersebut mem miliki nilai eexceedance frequency f l lebih tinggi dari sekitarn nya dikarenaakan masih teerdapat akum mulasi minyak yang m membuat lap pisan minyak k di area terssebut lebih tebal t dari sekkitarnya.
100 Dari penjelasan diatas, dapat diketahui pola proses pelapukan lapisan minyak dari masing-masing parameter yaitu nilai evaporasi, disolusi, dispersi vertikal dan emulsifikasi. Pada satu jam pertama setelah terjadinya tumpahan, lapisan avtur terlihat belum menyebar luas dari lokasi sumber. Luas permukaan lapisan minyak yang kecil menyebabkan proses pelapukan tidak dapat terlihat jelas. Pada jam ke-12, telah terlihat nilai pelapukan pada masing-masing parameter. Nilai tersebut semakin meningkat pada time step jam ke-24, kecuali pada proses dispersi vertikal. Selain itu, nilai exceedance frequency lapisan minyak pada time step tersebut juga mengalami penurunan. Setelah Jam ke-96, lapisan avtur tersebut tidak lagi mengalami proses pelapukan disebabkan seluruh lapisan minyak telah terbawa arus keluar dari domain model perairan Cilacap. Nilai yang terlihat pada time step tersebut hanya exceedance frequency yang kini memiliki nilai sangat kecil. Semakin mengecilnya laju perubahan konsentrasi fraksi lapisan minyak disebabkan karena komponen kimia yang terkandung di dalam lapisan minyak sebagian besar telah mengalami pelapukan.
4.5.2. Musim Timur 4.5.2.1.Jam ke-1 Gambar 49 merupakan sebaran serta proses pelapukan lapisan avtur yang tumpah di perairan Cilacap setelah jam ke-1 pada bulan Agustus 2007. Terdapat tiga sumber tumpahan minyak yang diasumsikan tumpah ke laut dimana jenis minyak, lokasi dan jumlah tumpahannya sama seperti pada musim barat. Lapisan minyak pertama berada di pantai selatan Cilacap dengan panjang lapisan mencapai 600 meter. Lapisan minyak kedua berada di sekitar mulut kanal utama dengan panjang lapisan mencapai 400 meter. Sedangkan lapisan minyak
101 ketiga masih berada di tengah-tengah perairan Teluk Penyu dengan panjang lapisan mencapai 400 meter. Proses penguapan yang terjadi pada bagian tengah lapisan avtur yang pertama mencapai 0.0062 mm, sedangkan ketebalan evaporasi pada lapisan kedua berada di bawah 0.0006 mm dan pada lapisan ketiga mencapai 0.0019 mm. Proses evaporasi pada lapisan pertama lebih besar dibandingkan dengan lapisan lainnya. Proses disolusi pada lapisan avtur pertama memiliki ketebalan antara 1.9 x 10-7 mm hingga 2.1 x 10-7 mm. Lapisan avtur kedua dan ketiga masing-masing memiliki ketebalan di bawah 3 x 10-8 mm dan 7 x 10-8 mm dalam proses disolusi. Hal ini berbanding lurus dengan tingkat evaporasi serta proses dispersi vertikal yang terjadi pada masing-masing lapisan minyak. Proses dispersi vertikal pada lapisan minyak yang pertama sebagian besar memiliki ketebalan di atas 6.9 x 10-5 mm. Lapisan minyak kedua mengalami proses dispersi vertikal sebanyak 5.7 x 10-6 – 6.3 x 10-6 mm, dan lapisan minyak ketiga mencapai lebih dari 6.9 x 10-6 mm pada bagian tengah lapisan. Tingkat emulsifikasi pada bagian tengah lapisan minyak pertama mencapai 225 - 250 mm. Sementara itu, lapisan minyak kedua mengalami proses emulsifikasi dengan ketebalan di bawah 25 mm. Lapisan minyak ketiga mengalami proses emulsifikasi dengan ketebalan antara 75-100 mm. Lapisan minyak pertama mengalami proses emulsifikasi yang lebih tinggi jika dibandingkan dengan proses emulsifikasi yang terjadi pada lapisan minyak kedua dan ketiga.
102 Laju perubahan konsentrasi k ffraksi lapisaan minyak jennis avtur padda satu j pertamaa menunjukk jam kan tingkat yyang tinggi, yaitu y lebih besar b dari 18% % di s seluruh lapissan yang tum mpah. Hal inni disebabkaan pada awall mula tumpaahan, l lapisan miny yak mengalaami proses peelapukan seh hingga konsentrasi lapisan minyak t tersebut teru us mengalam mi perubahann.
G Gambar 49. Sebaran serta Proses Peelapukan Lappisan Avtur JJam ke-1 di Perairan Cilacap padda Bulan Aguustus 2007
4 4.5.2.2.Jam ke-12 Sebaaran serta prooses pelapukkan lapisan avtur a yang tuumpah di perrairan C Cilacap setelah jam ke-112 pada bulaan Agustus 2007 2 disajikaan dalam Gaambar 50. L Lapisan avtu ur kedua yan ng berasal daari kebocoran n pengisian muatan tankker telah t tergabung deengan lapisaan avtur yangg pertama. Gabungan G daari kedua lappisan t tersebut men nyebar menuuju ke dalam m kanal utam ma. Panjang llintasan gabuungan
103 l lapisan minyyak tersebut mencapai em mpat km dann menyebar di tepi utaraa alian k kanal utamaa. Lapisan minyak m ketigaa telah menccapai pantai Cilacap denngan m menempuh perjalanan p seepanjang dua km dari tim mur ke baratt. Prosees penguapaan terjadi teruutama pada lapisan minyyak yang berrada di p pantai timurr Cilacap. Laapisan minyyak yang mem miliki panjanng hampir 6000 meter t tersebut men ngalami evapporasi dengaan ketebalann antara 0.00662 - 0.0068 mm, s sedangkan laapisan minyak lainnya mengalami m evvaporasi di bbawah 0.00112 mm.
G Gambar 50. Sebaran serta Proses Peelapukan Lappisan Avtur JJam ke-12 di d Perairan Cilacap padda Bulan Aggustus 2007 Prosees disolusi memiliki m keteebalan tertin nggi terutamaa pada lapisaan minyak y yang berasall dari kapal tanker t karam m yaitu di ataas 2.3 x 10-7 mm. Pada lapisan l a avtur lainnyaa hampir tiddak terlihat adanya proses disolusi diikarenakan luuas
104 permukaannya telah mengecil atau memiliki ketebalan disolusi di bawah 5 x 10-8 mm. Lapisan minyak di timur pantai Cilacap masih mengalami dispersi vertikal dengan ketebalan di bawah 1 x 10-6 mm. Proses emulsifikasi tertinggi terjadi pada lapisan minyak di sisi timur pantai Cilacap dengan ketebalan lebih dari 275 mm. Lapisan minyak lainnya mengalami emulsifikasi dengan ketebalan di bawah 50 mm. Laju perubahan konsentrasi fraksi tertinggi terjadi pada lapisan minyak di pantai timur dan selatan Cilacap hingga mencapai 15%.
4.5.2.3.Jam ke-24 Sebaran serta proses pelapukan lapisan avtur yang tumpah di perairan Cilacap setelah jam ke-24 pada bulan Agustus 2007 disajikan dalam Gambar 51. Lapisan minyak yang menyebar memasuki Sungai Kaliyasa melalui dermaga nelayan. Proses penguapan pada lapisan minyak tidak lagi terlihat, terkecuali pada lapisan minyak yang berada di dalam aliran Sungai Kaliyasa yang berada pada kisaran 0.0062 - 0.0068 mm. Hal yang sama juga tejadi pada proses disolusi dimana ketebalannya mencapai lebih dari 2.3 x 10-7 mm. Proses dispersi vertikal tidak lagi terlihat nilainya, namun proses emulsifikasi masih berlangsung dengan ketebalan antara 250 - 275 mm pada lapisan minyak di aliran Sungai Kaliyasa. Diperkirakan proses dispersi vertikal tersebut masih berlangsung namun memiliki ketebalan di bawah nilai kisaran terendah skala. Laju perubahan konsentrasi fraksi lapisan minyak di sekitar aliran kanal berada di bawah 1%, sedangkan lapisan minyak yang berada di Teluk Penyu
105 m masih mengalami laju perubahan koonsentrasi fraaksi minyak sebesar 5%. S Selanjutnya galami perubbahan lagi diikarenakan sebaran lapiisan avtur inni tidak meng l lapisan avturr tersebut tellah hilang daari kolom peerairan Cilaccap karena mengalami m p pelapukan.
G Gambar 51. Sebaran serta Proses Peelapukan Lappisan Avtur JJam ke-24 di d Perairan Cilacap padda Bulan Aggustus 2007 Dari penjelasan di d atas dapatt disimpulkaan bahwa padda musim tim mur l lapisan avturr lebih cepatt hilang atauu meninggalkkan permukaaan laut padaa domain m model diban ndingkan denngan musim barat. Keceepatan anginn yang bertiuup lebih b besar pada musim m timurr memperceppat proses peelapukan lapisan minyakk. Pada J ke-1, laapisan minyaak pertama m Jam memiliki ting gkat pelapukkan paling tinnggi jika d dibandingka an dengan tin ngkat pelapuukan pada lappisan minyaak kedua dann ketiga m meskipun juumlah tumpaahan lapisan minyak perttama memiliiki nilai yangg sama
106 dengan jumlah tumpahan lapisan ketiga. Tingginya nilai pelapukan pada lapisan minyak pertama dikarenakan lapisan minyak tersebut mendapat pengaruh dari posisinya yang dekat dengan daratan. Perairan yang dekat dengan daratan umumnya lebih dangkal dan lebih cepat menerima bahang. Selain itu perairan yang berada dekat dengan daratan juga lebih banyak mengalami turbulensi akibat gesekan dasar. Turbulensi tersebut dapat meningkatkan luas permukaan lapisan minyak sehingga mempercepat proses pelapukan. Tingginya tingkat pelapukan yang terjadi serta kondisi kecepatan arus pada kanal yang cukup besar menyebabkan volume lapisan minyak tersebut cepat berkurang. Kondisi yang sama juga terjadi pada lapisan minyak kedua dan terlihat pada time step Jam ke-12. Lapisan minyak ketiga pada time step tersebut justru memiliki nilai pelapukan yang paling tinggi. Pada time step tersebut, lapisan minyak ketiga baru saja sampai ke tepi pantai Cilacap. Selanjutnya pada jam ke24, hanya lapisan minyak ketiga saja yang masih berada dalam domain model dan mengalami proses pelapukan. Proses pelapukan yang terjadi pada lapisan minyak memiliki nilai yang sama dengan jam ke-12 sebelumnya, kecuali pada nilai exceedance frequency yang semakin menurun. Setelah jam ke-24 tersebut,
lapisan avtur tidak lagi terlihat berada di dalam kolom perairan.
4.6.
Pembahasan Proses Pelapukan Tumpahan Minyak di Laut Masing-masing jenis minyak yang dimodelkan tumpah di perairan Cilacap
memiliki tingkat pelapukan yang bervariasi sesuai dengan karakteristik masingmasing minyak. Pada proses pelapukan kali ini, akan dibahas kondisi ketebalan lapisan masing-masing jenis minyak setelah tumpah di perairan selama 12 jam.
107 Perbandingan proses evaporasi dan disolusi masing-masing jenis minyak pada musim barat dan musim timur ditampilkan pada Tabel 9. Proses evaporasi lapisan avtur dan minyak mentah cenderung lebih tinggi pada musim timur. Minyak mentah memiliki ketebalan evaporasi tertinggi yaitu 1.8333 mm pada musim barat dan 3.6 mm musim timur, sementara avtur memiliki ketebalan evaporasi sebesar 0.0056 mm pada musim barat dan 0.0068 mm pada musim timur. Jika dilihat dari pola penyebarannya, kedua lapisan minyak tersebut cenderung menyebar di dekat daratan pada musim timur. Lebih tingginya tingkat turbulensi lapisan minyak pada perairan di sekitar daratan menyebabkan peningkatan luas permukaan lapisan minyak yang akan meningkatkan ketebalan evaporasi pada awal terjadinya tumpahan. Selain itu posisi awal sumber tumpahan kedua jenis minyak yang cenderung berada di perairan terbuka mendapatkan pengaruh angin lebih besar pada musim timur sehingga meningkatkan proses evaporasi. Aspal memiliki ketebalan evaporasi sebesar 0.0183 mm pada musim barat dan 9.2 x 10-4 mm pada musim timur, sementara diesel dengan sumber masukan konstan memiliki ketebalan evaporasi mencapai 0.0053 mm pada musim barat dan 0.0041 pada musim timur. Suhu permukaan laut yang lebih tinggi pada musim barat lebih banyak mempengaruhi proses evaporasi lapisan minyak tipis seperti diesel sehingga cenderung lebih tinggi dibandingkan pada musim timur. Sementara pada musim timur, aspal lebih cepat keluar dari dalam domain sehingga proses evaporasi pada musim tersebut terlihat lebih rendah dibandingkan dengan musim barat.
108 Proses disolusi tertinggi terjadi pada lapisan minyak mentah dengan ketebalan lapisan mencapai 3.4 x 10-6 mm pada musim barat dan 4.6 x 10-6 mm pada musim timur. Berat molekul minyak mentah yang lebih besar dibandingkan jenis minyak lainnya memudahkan larutnya partikel minyak mentah ke dalam kolom perairan. Sementara lapisan minyak avtur dan diesel memiliki densitas yang lebih ringan sehingga cenderung menyebar rata di permukaan laut. Kondisi tersebut memungkinkan proses evaporasi berjalan lebih cepat dibandingkan dengan proses disolusi dimana keduanya sama-sama dialami oleh fraksi minyak aromatik. Lapisan avtur memiliki ketebalan disolusi sebesar 4.2 x 10-7 mm pada musim barat dan 2.3 x 10-7 pada musim timur, sementara lapisan aspal memiliki ketebalan disolusi 9.1 x 10-7 pada musim barat dan 2.3 x 10-8 pada musim timur. Lapisan diesel memiliki ketebalan disolusi sebesar 5 x 10-8 pada musim barat dan 1.3 x 10-8 pada musim timur. Seluruh jenis minyak mengalami proses disolusi tertinggi pada musim barat, kecuali pada minyak mentah. Turbulensi yang lebih tinggi pada kolom perairan menyebabkan lapisan minyak mentah mengalami peningkatan luas permukaan lapisan minyak yang dapat meningkatkan ketebalan disolusi pada musim timur. Namun pada lapisan avtur, peningkatan luas permukaan lapisan pada musim timur justru mempercepat terjadinya proses evaporasi dibandingkan dengan disolusi. Lapisan aspal pada musim barat lebih lama berada dalam domain sehingga nilai disolusi terlihat lebih tinggi.
Aspal
Crude
Avtur
Diesel
Tabel 9. Perb bandingan Prosees Evaporasi daan Disolusi Selurruh Jenis Minyaak yang Dimodeelkan Tumpah ddi Perairan Cilaccap pada Jam Kee12 pada Musim Baarat dan Timur Tahun T 2007 Musim Baarat Musim Timur Evaporassi Disolussi Evaporasi Diso olusi
109
109
110 Perbandingan proses emulsifikasi dan dispersi vertikal seluruh jenis minyak pada musim barat dan musim timur disajikan dalam Tabel 10. Proses emulsifikasi lapisan avtur dan minyak mentah memiliki nilai yang lebih tinggi pada musim timur. Minyak mentah memiliki ketebalan emulsifikasi tertinggi yaitu 75 mm pada musim barat dan 275 mm musim timur, sementara avtur memiliki ketebalan evaporasi sebesar 62.5 mm pada musim barat dan 275 mm pada musim timur. Kedua lapisan minyak tersebut cenderung menyebar menuju permukaan laut yang berada di dekat daratan pada musim timur. Lebih tingginya tingkat turbulensi lapisan minyak pada perairan di dekat daratan menyebabkan peningkatan intrusi air laut ke dalam lapisan minyak (water-uptake) sehingga meningkatkan proses emulsifikasi dan meningkatkan ketebalan lapisan minyak seluruhnya. Pengaruh angin yang lebih besar pada musim timur terutama pada perairan terbuka juga meningkatkan turbulensi permukaan laut sehingga meningkatkan proses emulsifikasi. Aspal memiliki ketebalan emulsifikasi sebesar 275 mm pada musim barat dan 25 mm pada musim timur, sementara diesel dengan sumber masukan konstan memiliki ketebalan evaporasi mencapai 0.037 mm baik pada musim barat maupun pada musim timur. Pada musim timur, aspal lebih cepat keluar dari dalam domain karena terdorong oleh arus dan angin permukaan sehingga proses emulsifikasi pada musim tersebut terlihat lebih rendah dibandingkan dengan musim barat. Proses emulsifikasi pada aspal banyak terjadi pada lapisan yang terdampar di sekitar tepi aliran Kali Donan. Sedangkan proses emulsifikasi pada diesel terjadi pada seluruh lapisan yang mengalami akumulasi. Nilai emulsifikasi yang rendah
111 pada lapisan diesel disebabkan karena diesel memiliki jumlah kandungan surfactant (aspal dan wax) dalam jumlah sedikit.
Proses dispersi vertikal tertinggi terjadi pada lapisan minyak mentah dengan ketebalan lapisan mencapai 1.8 x 10-6 mm pada musim barat dan 1 x 10-6 mm pada musim timur. Minyak mentah memiliki berat molekul yang lebih besar dibandingkan jenis minyak lainnya sehingga memudahkan masuknya partikel minyak mentah ke dalam kolom perairan. Selain itu, densitas partikel minyak mentah yang terdispersi memiliki nilai yang dapat mendekati densitas air laut di sekitarnya sehingga dapat mempertahankan kedudukan partikel untuk tetap berada dalam kolom perairan. Sementara itu, lapisan avtur dan diesel memiliki densitas yang lebih ringan sehingga cenderung kembali lagi ke lapisan minyak yang berada di permukaan laut. Hal yang sama terjadi pada lapisan avtur dimana memiliki ketebalan dispersi vertikal sebesar 1.6 x 10-6 mm pada musim barat dan 1 x 10-6 pada musim timur. Sementara itu, lapisan aspal memiliki ketebalan dispersi vertikal sebesar 9.2 x 10-7 pada musim barat. Pada musim timur, lapisan aspal telah terdispersi pada awal tumpahan sebelum akhirnya terbawa keluar dari domain perairan akibat pengaruh dari pergerakan arus dan angin permukaan. Lapisan diesel memiliki ketebalan dispersi vertikal sebesar 4.6 x 10-9 pada musim barat dan 4.5 x 10-9 pada musim timur. Proses dispersi vertikal pada lapisan diesel memiliki nilai yang sangat kecil dikarenakan masukan minyak tersebut ke dalam permukaan laut juga bernilai sangat kecil. Seluruh jenis minyak mengalami proses dispersi vertikal tertinggi pada musim barat dikarenakan pada musim timur lapisan minyak lebih banyak mengalami emulsifikasi sehingga menyebabkan viskositas lapisan minyak meningkat.
Aspal
Crude
Avtur
Diesel
Tabel 10. Perrbandingan Prosses Emulsifikasii dan Dispersi V Vertikal Seluruh Jenis Minyak yyang Dimodelkaan Tumpah di Peerairan Cilacap padda Jam Ke-12 paada Musim Barat dan Timur Taahun 2007 Musim Timur Musim Baarat Emulsifikaasi Dispersi Veertikal Emulsiffikasi Dispersi Vertikal
112
112
113 Laju perubahan konsentrasi fraksi minyak (exceedance frequency) dan time exposure seluruh jenis minyak pada musim barat dan musim timur
ditampilkan pada Tabel 11. Baik pada musim barat maupun pada musim timur, lapisan diesel tidak terlihat memiliki nilai exceedance frequency. Hal tersebut dikarenakan laju perubahan fraksi minyak yang tumpah setelah 12 jam tersebut masih sangat kecil. Lapisan diesel diasumsikan memasuki perairan dengan volume yang sangat kecil namun konstan sehingga menyebabkan ketebalan lapisan minyak setelah tumpah selama 12 jam belum terlalu signifikan dibandingkan dengan total diesel yang diskenariokan tumpah di perairan Cilacap. Avtur memiliki nilai exceedance frequency sebesar 4.1% pada musim barat dan 15% pada musim timur, sedangkan minyak mentah memiliki nilai exceedance frequency sebesar 4.5% pada musim barat dan 13% pada musim timur. Nilai exceedance frequency pada musim timur lebih tinggi jika dibandingkan dengan
musim barat dikarenakan minyak yang tumpah pada musim tersebut banyak terakumulasi di dalam domain perairan. Akumulasi minyak tersebut tidak terlalu besar jumlahnya dibandingkan dengan total minyak avtur maupun minyak mentah yang dimodelkan tumpah di perairan Cilacap. Nilai exceedance frequency pada aspal bernilai sebesar 9% pada musim barat dan 2.5-9% pada musim timur. Nilai exceedance frequency pada musim timur lebih besar jika dibandingkan dengan
musim barat dikarenakan pada musim timur, sebagian besar lapisan aspal telah terbawa ke luar domain akibat pergerakan arus dan angin yang cukup besar. Konsentrasi fraksi minyak yang tertinggal dalam domain bernilai kecil jika dibandingkan dengan total volume aspal yang diasumsikan tumpah di perairan Cilacap dikarenakan aspal juga terlah mengalami proses pelapukan. Keseluruhan
114 nilai exceedance frequency umumnya tinggi di sekitar sumber tumpahan dan di bagian pusat lapisan. Hal tersebut disebabkan karena lapisan minyak memiliki ketebalan yang lebih besar di sekitar sumber tumpahan maupun pada bagian tengah lapisan. Lapisan minyak yang terakumulasi juga memiliki nilai exceedance frequency yang lebih tinggi dari lapisan minyak di sekitarnya. Time exposition merupakan lamanya suatu grid terkena atau terpapar oleh
minyak. Setelah 12 jam, tumpahan diesel terlihat telah bergerak ke luar kanal utama pada musim barat, sedangkan pada musim timur lapisan tersebut bergerak pula ke dalam kanal. Pada musim barat, lapisan avtur terlihat bergerak menuju ke timur domain yang berbatasan langsung dengan Samudera Hindia, sedangkan pada musim timur lapisan tersebut bergerak ke dalam domain. Hal yang sama terjadi pada tumpahan minyak mentah. Pada musim barat pesisir Pulau Nusakambangan menjadi daerah yang paling rawan terkena risiko pencemaran minyak, sedangkan pada musim timur lapisan minyak menyebar di sekitar pantai timur Cilacap, disepanjang aliran kanal utama dan aliran Kali Donan. Lapisan aspal pada musim barat terlihat menyebar disekitar tepi timur aliran Kali Donan, sedangkan pada musim timur lapisan tersebut menyebar di tepi sungai bagian barat. Keseluruhan nilai time exposition yang dimodelkan pada tumpahan diesel mencapai 240 jam pada musim barat dan musim timur. Lapisan avtur dan minyak mentah memiliki total time exposition hingga 96 jam pada musim barat dan 24 jam pada musim timur, sementara lapisan aspal memiliki nilai time exposition hingga 120 jam pada musim barat dan 48 jam pada musim timur. Nilai time exposition lebih pendek pada musim timur jika dibandingkan dengan musim barat
dikarenakan pada musim tersebut lapisan minyak lebih cepat hilang dari domain
115 perairan akibat terbawa oleh arus dan terdorong oleh angin permukaan yang cukup besar. Nilai time exposition paling besar terjadi pada lapisan diesel dikarenakan diesel diasumsikan memasuki perairan secara konstan sehingga tetap berada dalam domain. Lapisan avtur, minyak mentah dan aspal sebagian besar akan terbawa keluar dari domain perairan namun tetap menyisakan sejumlah minyak yang terdampar di daratan dalam domain model.
Aspal
Crude
Avtur
Diesel
Tabel 11. Perrbandingan Exceeedance Frequeency dan Time Exposure E Seluruuh Jenis Minyak yang Dimodelk kan Tumpah di Perairan P Cilacapp padda Jam Ke-12 paada Musim Barat dan Timur Taahun 2007 Musim Timur Musim Baarat E Exceedance Freequency Time Expoosure Exceedance F Frequency Time Exxposure
116
116
5. KESIMPULAN DAN SARAN 5.1.
Kesimpulan Model hidrodinamika dan tumpahan minyak yang dimodelkan dalam studi
ini dapat diterima dengan baik karena data masukan pemodelan dan data hasil pengukuran insitu memiliki kemiripan tinggi. Perairan di sekitar pesisir Pulau Nusakambangan menjadi daerah yang paling rawan terkena dampak pencemaran pada musim barat karena pergerakan arus dan arah angin pada musim tersebut mengarah ke bagian tenggara dan timur domain. Pada musim timur, daerah pantai timur Cilacap serta di sepanjang aliran kanal utama dan Kali Donan memiliki risiko tertinggi terhadap pencemaran minyak karena arah arus dan angin dominan mengarah ke bagian barat dan barat laut domain. Proses evaporasi pada lapisan minyak yang lebih tebal (minyak mentah dan avtur) semakin meningkat seiring dengan peningkatan luas permukaan akibat dari turbulensi, sedangkan proses evaporasi pada lapisan minyak yang lebih tipis (diesel) lebih banyak mendapat pengaruh dari suhu lingkungan sekitar. Proses disolusi pada musim barat lebih tinggi dibandingkan dengan musim timur diduga karena pada musim timur fraksi minyak dengan berat molekul rendah (aromatic) cenderung lebih dahulu terevaporasi daripada terdisolusi. Proses emulsifikasi seluruh jenis minyak memiliki nilai tertinggi pada musim timur dibandingkan dengan musim barat diduga karena pada musim timur lapisan minyak tersebut lebih banyak mengalami turbulensi. Proses dispersi vertikal pada musim barat justru lebih tinggi dibandingkan pada musim timur
117
118 diduga karena pada musim timur lapisan minyak telah mengalami peningkatan viskositas yang cukup besar akibat proses emulsifikasi. Laju perubahan konsentrasi fraksi minyak bernilai tinggi di sekitar sumber tumpahan dan di bagian pusat lapisan karena lapisan minyak memiliki ketebalan yang lebih besar pada bagian lapisan tersebut. Kecilnya nilai laju perubahan fraksi minyak pada semua jenis minyak menandakan lapisan minyak telah menyebar, kecuali pada lapisan minyak yang terdampar di sekitar domain. Resident time pada musim timur lebih pendek jika dibandingkan dengan
musim barat karena pada musim tersebut lapisan minyak lebih cepat hilang dari domain perairan akibat terbawa oleh arus dan terdorong oleh angin permukaan yang cukup besar.
5.2.
Saran Saran yang dapat diberikan dalam penelitian ini yaitu perlu dilakukan
penelitian serupa dengan memasukkan faktor gelombang ke dalam model hidrodinamika. Faktor gelombang sangat berperan terhadap proses pelapukan minyak seperti emulsifikasi dan dispersi vertikal. Selain itu dalam penelitian selanjutnya dapat dilakukan verifikasi hasil model baik dari segi hidrodinamika maupun dari segi oil spill. Pemodelan dengan menggunakan input jenis minyak yang lain sangat disarankan untuk mendapatkan hasil perbandingan yang lebih beragam, sedangkan untuk mendapatkan hasil perbandingan tingkat pelapukan pada masing-masing jenis minyak secara nyata dapat dilakukan dengan menyamakan jumlah tumpahan minyak, waktu tumpahan, serta lokasi sumber tumpahan minyak.
DAFTAR PUSTAKA Alonso, M dan E.J. Finn. 1992. Dasar-dasar fisika universitas. Penerbit Erlangga. Jakarta. Clark, R. B. 1986. Marine Pollution. Clarendon Press-Oxford. UK. CRC. 2000. Chapter 4: Behaviour of Oil in The Environment. CRC Press LCC. DHI. 2006. Mike 21 HD Scientific Documentation. DHI Water & Environment. Denmark. DHI. 2006. Mike 21 & Mike 3 PA/SA: Particle Analysis and Oil Spill Analysis Module. DHI Water & Environment. Denmark. DHI. 2007. Mike 21 Flow Model FM: Hydrodynamic Module. DHI Water & Environment. Denmark. Firdaus. 1997. Pemodelan dan simulasi komputer pola arus dan trayektori tumpahan minyak di Perairan Cilacap dengan metode beda hingga eksplisit. Fakultas Perikanan dan lmu Kelautan. Skripsi: Tidak dipublikasikan. Institut Pertanian Bogor. Gross, M.G. 1972. Oceanography a View The Earth. Prentice Hall International Inc. London, UK. Harimurthy, S. 2002. Tipologi komunitas makrozoobentos sebagai bioindikator pencemaran perairan di muara Sungai Donan, Cilacap, Jawa Tengah. Fakultas Perikanan dan Ilmu Kelautan. Skripsi: Tidak dipublikasikan. Institut Pertanian Bogor. IFREMER. 2007. Sumber Data Angin. ftp://ftp.ifremer.fr/pub/ifremer/cersat/products/gridded/mwf-blended/data/ [1 Oktober 2008] Ilham. 2002. Identifikasi permasalahan dan langkah penanganan mitigasinya. Jurusan Studi Ilmu Lingkungan Program Pasca Sarjana. Universitas Gajah Mada. ITOPF. 2007. About Marine Spills – Weathering Process. http://www.itopf.com/marine%2Dspills/fate/weathering%2Dprocess/ [27 Maret 2008] JANHIDROS. 2007. Peta batimetri Jawatan Hidro-Oseanografi TNI-AL Nomor 108. JANHIDROS, Jakarta. Mukhtasor. 2007. Pencemaran pesisir dan laut. PT Pradnya Paramita, Jakarta.
119
120 Pariwono, J. I. 1989. Kondisi pasang surut di Indonesia. in Pasang Surut. Penyunting: O. S. R. Ongkosongko dan Suyarso. Pusat Penelitian dan Pengembangan Oseanologi Lembaga Ilmu Pengetahuan. Pertamina, 2006. Warta Pertamina: Avtur sei pakning potensial lebih efisien. http://www.pertamina.com/index.php?option=com_content&task=view&i d=1312&Itemid=507 [27 Juli 2009] Wikipedia. 2007. Kabupaten Cilacap. http://id.wikipedia.org/wiki/Kabupaten_Cilacap [27 Maret 2008]
LAMPIRAN
121
122 Lampiran 1. Contoh Laporan Hasil Pemodelan Hidrodinamika Pada Musim Barat dengan Menggunakan DHI Software Mike 21Hydrodynamic Modul 20080925 184906 ========================================================== ============Target: MIKE21_NESTED_MODEL Section: NESTED_MODEL_GLOBALS -------------------------------------Section: MODULE_SELECTION ---------------------------------------Simulation_Type : 0 Include_DB_computations : False Section: HYDRODYNAMIC_MODULE --------------------------------------Section: OPTION_PARAMETERS --------------------------------------Print_Misc_Output : 0 0 0 0 0 0 0 0 0 0 (default) Section: BATHYMETRY_SELECTION -----------------------------------No_of_Areas : 1 Hot_Start : 0 Section: BATHYMETRY_SELECTION -----------------------------------Section: AREA_1 -----------------------------------------------Section: DATA_FILE ------------------------------------------File_Name : G:\Rizka Cilacap\SC01\Bathy\Bathy.dfs2 Item_Numbers : 1 (default) Section: OPTION_PARAMETERS --------------------------------------dzdt_File_Name : (default) Section: BATHYMETRY_SELECTION -----------------------------------Projection : PROJCS["WGS_1984_UTM_Zone_49S",GEOGCS["GCS_WGS_1984",DATU M["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIM EM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["T ransverse_Mercator"],PARAMETER["False_Easting",500000],PARAMETER["F alse_Northing",10000000],PARAMETER["Central_Meridian",111],PARAMETE R["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0],UNIT["Meter ",1]] CoriolisForce : True Land_Slides : False Section: SIMULATION_PERIOD --------------------------------------Number_of_Timesteps : 466560 Time_Step_Interval : 5 Warm_Up_Period : 7200 Start_Time : 2007 2 1 0 0 0 Section: SURFACE_ELEVATION --------------------------------------Section: AREA_1 -----------------------------------------------Type : 1 Constant_Value : 0.12 Section: FLOOD_AND_DRY ------------------------------------------Enable_Flood_and_Dry_Checking : 1 Drying_Depth : 0.2
123 Flooding_Depth : 0.3 Section: OPTION_PARAMETERS --------------------------------------Flood_Dry_Time_Step_Interval : 1 (default) Section: OPTION_PARAMETERS --------------------------------------fldhw : 0.6 (default) Section: OPTION_PARAMETERS --------------------------------------epsf : 0.000199 (default) Section: OPTION_PARAMETERS --------------------------------------epsh : 0.000099 (default) Section: OPTION_PARAMETERS --------------------------------------RainOnDryLand : False (default) Section: RESISTANCE ---------------------------------------------Formulation : 0 Section: AREA_1 -----------------------------------------------Data_Format : 2 Section: DATA_FILE ------------------------------------------File_Name : G:\Rizka Cilacap\SC01\Ress\Manning.dfs2 Item_Numbers : 1 Include_Pier : False --------------------------- Data File Summary --------------------------Data origin: Grid editor Title : G:\@-Cil_Oil\Pros\bathy_ascii.dfs2 Delete : -1.000E-30 Axis Sets Interval Axis origin Unit 1 885 10 0 meter 2 935 10 0 meter Item name Item type Minimum Maximum Missing Unit Ressistance Manning's M 8 32 0 m^(1/3)/s ---------------------------------------------------------------------Section: OPTION_PARAMETERS --------------------------------------Reference_Point_in_Area_1 : -99 0 (default) Section: EDDY_VISCOSITY -----------------------------------------Data_Format : 3 Formulation : 1 Section: EDDY_VISCOSITY -----------------------------------------Section: AREA_1 -----------------------------------------------Smagorinsky_Format : 0.5 Section: OPTION_PARAMETERS --------------------------------------Smag_File_Name_Area_1 : (default) Section: OPTION_PARAMETERS --------------------------------------MinimumEddy_AREA_1 : 0 (default) Section: OPTION_PARAMETERS --------------------------------------MaximumEddy_AREA_1 : 200 (default) Section: OPTION_PARAMETERS --------------------------------------MinEddyDepth : 50 (default) Section: WAVE_RADIATION -----------------------------------------Section: AREA_1 -----------------------------------------------Included : False
124 Section: OPEN_BOUNDARY ------------------------------------------Number_Of_Boundaries : 4 Section: OPEN_BOUNDARY ------------------------------------------Section: BOUNDARY_1 -------------------------------------------First_Point : 0 357 Last_Point : 0 495 Section: OPEN_BOUNDARY ------------------------------------------Section: BOUNDARY_1 -------------------------------------------Formulation : 1 Data_Format : 4 FAB_Type : 12 Tilt_Type : 1 Tilting_Point : 473 Section: USER_DIRECTIONS ------------------------------------Enable_User_Directions : False Section: DATA_FILE ------------------------------------------File_Name : G:\Rizka Cilacap\SC01\Pasut\Barat-Kiri.dfs1 --------------------------- Data File Summary --------------------------Data origin: //apv Title : G:\@-Cil_Oil\SC02\Pasut\Barat-Kiri.dfs1 Delete : -1.000E-30 Axis Sets Interval Axis origin Unit time 8065 300 2007-02-01 00:00:00 1 139 10 0 meter Item name Item type Minimum Maximum Missing Unit Predicted Water Level -0.775733 0.853062 0 meter ---------------------------------------------------------------------Section: OPEN_BOUNDARY ------------------------------------------Section: BOUNDARY_2 -------------------------------------------First_Point : 117 934 Last_Point : 169 934 Section: OPEN_BOUNDARY ------------------------------------------Section: BOUNDARY_2 -------------------------------------------Formulation : 1 Data_Format : 4 FAB_Type : 12 Tilt_Type : 1 Tilting_Point : 143 Section: USER_DIRECTIONS ------------------------------------Enable_User_Directions : False Section: DATA_FILE ------------------------------------------File_Name : G:\Rizka Cilacap\SC01\Pasut\Barat-Atas.dfs1 --------------------------- Data File Summary --------------------------Data origin: //apv Title : G:\@-Cil_Oil\SC02\Pasut\Barat-Atas.dfs1 Delete : -1.000E-30 Axis Sets Interval Axis origin Unit time 8065 300 2007-02-01 00:00:00
125 1 53 10 0 meter Item name Item type Minimum Maximum Missing Unit Predicted Water Level -0.749254 0.833902 0 meter ---------------------------------------------------------------------Section: OPEN_BOUNDARY ------------------------------------------Section: BOUNDARY_3 -------------------------------------------First_Point : 884 0 Last_Point : 884 871 Section: OPEN_BOUNDARY ------------------------------------------Section: BOUNDARY_3 -------------------------------------------Formulation : 1 Data_Format : 4 FAB_Type : 12 Tilt_Type : 1 Tilting_Point : 191 Section: USER_DIRECTIONS ------------------------------------Enable_User_Directions : False Section: DATA_FILE ------------------------------------------File_Name : G:\Rizka Cilacap\SC01\Pasut\Barat-Kanan.dfs1 --------------------------- Data File Summary --------------------------Data origin: //apv Title : G:\@-Cil_Oil\SC02\Pasut\Barat-Kanan.dfs1 Delete : -1.000E-30 Axis Sets Interval Axis origin Unit time 8065 300 2007-02-01 00:00:00 1 872 10 0 meter Item name Item type Minimum Maximum Missing Unit Predicted Water Level -0.793907 0.869539 0 meter ---------------------------------------------------------------------Section: OPEN_BOUNDARY ------------------------------------------Section: BOUNDARY_4 -------------------------------------------First_Point : 729 0 Last_Point : 884 0 Section: OPEN_BOUNDARY -----------------------------------------Section: BOUNDARY_4 -------------------------------------------Formulation : 1 Data_Format : 4 FAB_Type : 12 Tilt_Type : 1 Tilting_Point : 858 Section: USER_DIRECTIONS ------------------------------------Enable_User_Directions : False Section: DATA_FILE ------------------------------------------File_Name : G:\Rizka Cilacap\SC01\Pasut\Barat-Bawah.dfs1 --------------------------- Data File Summary --------------------------Data origin: //apv Title : G:\@-Cil_Oil\SC02\Pasut\Barat-Bawah.dfs1 Delete : -1.000E-30
126 Axis Sets Interval Axis origin Unit time 8065 300 2007-02-01 00:00:00 1 156 10 0 meter Item name Item type Minimum Maximum Missing Unit Predicted Water Level -0.793929 0.869562 0 meter ---------------------------------------------------------------------Section: OPTION_PARAMETERS --------------------------------------Bingham_Fluid : False (default) Section: OPTION_PARAMETERS --------------------------------------Fixedresimap : False (default) Section: OPTION_PARAMETERS --------------------------------------gID_lower : 0.05 (default) Section: OPTION_PARAMETERS --------------------------------------gID_upper : 0.059999 (default) Section: OPTION_PARAMETERS --------------------------------------Chezymin : 5 (default) Section: OPTION_PARAMETERS --------------------------------------Chezymax : 30 (default) Section: OPTION_PARAMETERS --------------------------------------Include_Structures : False (default) Section: OPTION_PARAMETERS --------------------------------------Read_decoupled_files : False (default) Section: WIND_CONDITIONS ----------------------------------------Data_Format : 3 Neutral_Pressure : 1013 Section: DATA_FILE --------------------------------------------File_Name : G:\Rizka Cilacap\SC01\Angin\Angin-01-02-07-Barat-FN.dfs2 IncludeAirPressureCorrections : False (default) --------------------------- Data File Summary --------------------------Data origin: Grid editor Title : File Title Delete : -1.000E-30 Axis Sets Interval Axis origin Unit time 112 21600 2007-02-01 00:00:00 1 885 10 0 meter 2 935 10 0 meter Item name Item type Minimum Maximum Missing Unit SL. Pressure Air Pressure 1009.90002 1012.59997 0 hectoPascal U-Wind Wind speed -2.6106 5.9242 0 m/s V-Wind Wind speed -3.4601 2.771 0 m/s ---------------------------------------------------------------------Section: WIND_CONDITIONS ----------------------------------------Type_of_Wind_Friction : 1 Linear_Friction_1 : 0.0016 Linear_Friction_2 : 0.0026 Linear_Speed_1 : 0 Linear_Speed_2 : 16 Section: SOURCE_AND_SINK -----------------------------------------
127 Number_Of_Sources : 0 Section: SOURCE_AND_SINK ----------------------------------------PrecipitationType : -1 Precipitation : 0 NetPrecipitation : False RainOnDryLand : False EvaporationType : -1 Evaporation : 0 Section: MASS_BUDGET --------------------------------------------NoOfPolygons : 0 Section: OUTPUT_SPECIFICATION -----------------------------------Number_of_Output_Areas : 1 Generate_Hotstart : False Section: OUTPUT_SPECIFICATION -----------------------------------Section: OUTPUT_AREA_1 ----------------------------------------File_Name : G:\Rizka Cilacap\SC01\HD\HD-Barat-PL.dfs2 Type : 2 X_Range_And_Interval : 0 884 1 Y_Range_And_Interval : 0 934 1 Area_No : 1 First_Time_Step : 1440 Last_Time_Step : 466560 Time_Step_Interval : 180 Output_Item : 3 ITEM_NUMBERS : 4 5 6 Title : HD-Barat Section: OUTPUT_SPECIFICATION -----------------------------------Section: INUNDATION_STATISTICS --------------------------------IncludeInundationStatistics : False InundationFileName : InundationFileStoreStep : 0 Section: KALMAN_FILTER_MODULE -------------------------------------Section: OPTION_PARAMETERS --------------------------------------bdel21 : True (default) Section: OPTION_PARAMETERS --------------------------------------Include_DB_computations : False (default) Section: OPTION_PARAMETERS --------------------------------------Include_Upwinding : True (default) Section: OPTION_PARAMETERS --------------------------------------Fr1 : 0.25 (default) Section: OPTION_PARAMETERS --------------------------------------Fr2 : 1 (default) Section: OPTION_PARAMETERS --------------------------------------alpha0 : 0 (default) Section: OPTION_PARAMETERS --------------------------------------Include_Forward_Centering_Mass_Eq : False (default) Section: OPTION_PARAMETERS --------------------------------------Include_Dissipative_Interface : False (default)
128 Section: OPTION_PARAMETERS --------------------------------------Dissipation_coefficient : 0.5 (default) Section: OPTION_PARAMETERS --------------------------------------No_Border_FLD : False (default) Section: OPTION_PARAMETERS --------------------------------------IncludeInundationStatistics : False (default) ======================== COMPUTATION STARTED ========================= Section: HYDRODYNAMIC_MODULE --------------------------------------Section: OPTION_PARAMETERS --------------------------------------Atm_Pressure_Correction : 0 (default) Section: OPTION_PARAMETERS --------------------------------------Logical_Bamc : True (default) Section: OPTION_PARAMETERS --------------------------------------Logical_Bamc : True (default) ========================= COMPUTATION ENDED ========================== ========================= Output Statistics ========================== --------------------------- Data File Summary --------------------------Data origin: Title : HD-Barat Delete : -1.000E-30 Axis Sets Interval Axis origin Unit time 2585 900 2007-02-01 02:00:00 1 885 10 0 meter 2 935 10 0 meter Item name Item type Minimum Maximum Missing Unit Surface elevation Surface Elevation -0.800179 0.868936******** meter U velocity u-velocity component -0.956005 0.692321******** m/s V velocity v-velocity component -0.798391 1.547274******** m/s ---------------------------------------------------------------------========================================================== ============ =========================== Run Statistics =========================== Total number of computational points : 176553115648. CPU time : 269280 [s] Points/sec : 655648.8 ========================================================== ============ Normal run completion
129 Lampiran 2. Contoh Laporan Hasil Pemodelan Tumpahan Minyak Menggunakan DHI Software Mike 21 Spill Analysis Modul 20090129 013603 =============================================================== ====== Target: MIKE321_PA_SA Section: BASIC_PARAMETERS -----------------------------------------Section: OPTION_PARAMETERS --------------------------------------Print_Misc_Output : 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 (default) Specify_Neutral_Dispersion : 0 (default) Include_Vert_Disp_From_Wind : False (default) Include_mud_consolidation : 1 (default) Include_Misc_Timeserie : False (default) Section: STARTING_CONDITIONS ------------------------------------Module_Selection : 1 Number_of_Computational_Areas : 1 Section: SIMULATION_PERIOD --------------------------------------Number_Of_Time_Steps : 103680 Time_Interval : 10 Start_Time : 2007 2 1 2 0 0 Section: HYDRODYNAMIC_DATA --------------------------------------Data_Selection : 0 Section: CURRENT_CONDITIONS -----------------------------------HD_Selection : 2 Section: CURRENT_CONDITIONS -----------------------------------Section: AREA_1 ---------------------------------------------Section: DATA_FILE ----------------------------------------File_Name : G:\Rizka Cilacap\SC01\HD\HD-Barat-Oil.dfs2 Working on: G:\Rizka Cilacap\SC01\HD\HD-Barat-Oil.dfs2 --------------------------- Data File Summary --------------------------Data origin: C:\Program Files\DHI\MIKE Title : HD-Barat Delete : -1.000E-30 Axis Sets Interval Axis origin Unit time 2585 900 2007-02-01 02:00:00 1 885 10 0 meter 2 935 10 0 meter Item name Item type Minimum Maximum Missing Unit H Water Depth m Water Level 0.925467 25.663668******** meter P Flux m^3/s/m Flow Flux -21.332455 14.248095******** m^3/s/m Q Flux m^3/s/m Flow Flux -12.573159 37.359356******** m^3/s/m ---------------------------------------------------------------------Projection : PROJCS["WGS_1984_UTM_Zone_49S",GEOGCS["GCS_WGS_1984",DATUM["D_ WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwic h",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"], PARAMETER["False_Easting",500000],PARAMETER["False_Northing",10000000],P ARAMETER["Central_Meridian",111],PARAMETER["Scale_Factor",0.9996],PARAM ETER["Latitude_Of_Origin",0],UNIT["Meter",1]]
130 Section: SOURCES ------------------------------------------------Number_of_Sources : 5 Number_of_Particles_per_Timestep : 10 Include_Max_Age : 0 Section: SOURCE_NO_1 ------------------------------------------Area_No : 1 Section: SOURCE_COORDINATES ---------------------------------Data_Type : 0 Source_Coordinates : 323 247 Section: SOURCE_DEPTH ---------------------------------------Data_Type : 0 Source_Depth_Below_Datum : 0 Section: SOURCE_FLUX ----------------------------------------Data_Type : 0 Mass_Flux : 0.0003 Section: SOURCE_NO_2 ------------------------------------------Area_No : 1 Section: SOURCE_COORDINATES ---------------------------------Data_Type : 0 Source_Coordinates : 93 450 Section: SOURCE_DEPTH ---------------------------------------Data_Type : 0 Source_Depth_Below_Datum : 0 Section: SOURCE_FLUX ----------------------------------------Data_Type : 1 Section: DATA_FILE ----------------------------------------File_Name : G:\Rizka Cilacap\SC01\flux\flux7b.dfs0 Item_Numbers : 1 Working on: G:\Rizka Cilacap\SC01\flux\flux7b.dfs0 --------------------------- Data File Summary --------------------------Data origin: MIKE Zero Title : Flux Barat Delete : -1.000E-30 Axis Sets Interval Axis origin Unit time232561 10 2007-02-01 02:00:00 Item name Item type Minimum Maximum Missing Unit Flux Volume Flux 0 1.5 0 m^3/s ---------------------------------------------------------------------Section: SOURCE_NO_3 ------------------------------------------Area_No : 1 Section: SOURCE_COORDINATES ---------------------------------Data_Type : 0 Source_Coordinates : 131 499 Section: SOURCE_DEPTH ---------------------------------------Data_Type : 0 Source_Depth_Below_Datum : 0 Section: SOURCE_FLUX ----------------------------------------Data_Type : 0 Mass_Flux : 0.000027 Section: SOURCE_NO_4 ------------------------------------------Area_No : 1 Section: SOURCE_COORDINATES ----------------------------------
131 Data_Type : 0 Source_Coordinates : 447 518 Section: SOURCE_DEPTH ---------------------------------------Data_Type : 0 Source_Depth_Below_Datum : 0 Section: SOURCE_FLUX ----------------------------------------Data_Type : 0 Mass_Flux : 0.000007 Section: SOURCE_NO_5 ------------------------------------------Area_No : 1 Section: SOURCE_COORDINATES ---------------------------------Data_Type : 0 Source_Coordinates : 731 660 Section: SOURCE_DEPTH ---------------------------------------Data_Type : 0 Source_Depth_Below_Datum : 0 Section: SOURCE_FLUX ----------------------------------------Data_Type : 1 Section: DATA_FILE ----------------------------------------File_Name : G:\Rizka Cilacap\SC01\flux\flux12b.dfs0 Item_Numbers : 1 Working on: G:\Rizka Cilacap\SC01\flux\flux12b.dfs0 --------------------------- Data File Summary --------------------------Data origin: MIKE Zero Title : Flux Barat Delete : -1.000E-30 Axis Sets Interval Axis origin Unit time232561 10 2007-02-01 02:00:00 Item name Item type Minimum Maximum Missing Unit Flux Volume Flux 0 0.079999 0 m^3/s ---------------------------------------------------------------------Section: DISPERSION ---------------------------------------------Type_of_Dispersion : 1 Section: LONGITUDINAL_DISPERSION ------------------------------Section: AREA_1 ---------------------------------------------Longitudinal_Proportionality_Factor : 1 Max_Longitudinal_Dispersion_Coefficient : 1 Min_Longitudinal_Dispersion_Coefficient : 0 Section: TRANSVERSAL_DISPERSION -------------------------------Section: AREA_1 ---------------------------------------------Transversal_Proportionality_Factor : 1 Max_Transversal_Dispersion_Coefficient : 1 Min_Transversal_Dispersion_Coefficient : 0 Section: VERTICAL_DISPERSION ----------------------------------Section: AREA_1 ---------------------------------------------Vertical_Proportionality_Factor : 0.009999 Max_Vertical_Dispersion_Coefficient : 0.009999 Min_Vertical_Dispersion_Coefficient : 0 Section: WIND_CONDITIONS ----------------------------------------Type_of_Wind : 3 Section: WIND_CONDITIONS ----------------------------------------Section: DATA_FILE ---------------------------------------------
132 File_Name : G:\Rizka Cilacap\SC01\Angin\Angin-01-02-07-Barat-FN.dfs2 Working on: G:\Rizka Cilacap\SC01\Angin\Angin-01-02-07-Barat-FN.dfs2 --------------------------- Data File Summary --------------------------Data origin: Grid editor Title : File Title Delete : -1.000E-30 Axis Sets Interval Axis origin Unit time 112 21600 2007-02-01 00:00:00 1 885 10 0 meter 2 935 10 0 meter Item name Item type Minimum Maximum Missing Unit SL. Pressure Air Pressure 1009.90002 1012.59997 0 hectoPascal U-Wind Wind speed -2.6106 5.9242 0 m/s V-Wind Wind speed -3.4601 2.771 0 m/s ---------------------------------------------------------------------Section: WIND_CONDITIONS ----------------------------------------Type_of_Wind_Friction : 1 Lower_Limit_Speed : 0 Lower_Limit_Friction : 0.015 Upper_Limit_Speed : 24 Upper_Limit_Friction : 0.02 Depth_of_Wind_Influence : 1 Deflection_Angle : 15 Depth_of_Longshore_Current_Zone : 5 Section: VELOCITY_AND_EDDY_PROFILE ------------------------------Include_Logarithmic_Velocity_Profile : 0 Include_Parabolic_Eddy_Profile : 1 Section: EXCEEDING_CONCENTRATION --------------------------------Include_Exceeding_Concentration : 1 Exceeding_Concentration_Limit : 100 Section: TIME_EXPOSITION ----------------------------------------Include_Time_Exposition : 1 Section: WATER_PROPERTIES ---------------------------------------Section: WATER_TEMPERATURE ------------------------------------Data_Type : 0 Section: AREA_1 ---------------------------------------------Constant_Water_Temperature : 29 Section: WATER_SALINITY ---------------------------------------Data_Type : 0 Section: AREA_1 ---------------------------------------------Constant_Water_Salinity : 33.5 Section: OIL_SPILL_PARAMETERS -------------------------------------Section: AIR_PROPERTIES -----------------------------------------Section: TEMPERATURE ------------------------------------------Data_Type : 0 Constant_Air_Temperature : 27.565 Section: CLOUDINESS -------------------------------------------Data_Type : 0 Constant_Cloudiness : 0.58 Section: HEAT_BALANCE -------------------------------------------Include_Evaporation : 1 Include_Heat_Balance : 1
133 Albedo_Value : 0.14 Emissivity_of_Oil : 0.82 Emissivity_of_Water : 0.95 Emissivity_of_Air : 0.82 Evaporation_Constant : 0.029 Section: EMULSIFICATION -----------------------------------------Include_Emulsification : 1 Maximum_Water_Content : 0.85 Asphaltenes_Content : 1 Wax_Content : 2 K1_Due_to_Water_Uptake : 5.E-07 K2_Due_to_Water_Release : 0.000012 Section: DISSOLUTION --------------------------------------------Include_Dissolution : 1 Mass_Transfer_Coefficient : 0.000002 Include_Entrainment : 1 Oil_Water_Tension : 29.9 Section: OIL_PROPERTIES -----------------------------------------Section: SOURCE_NO_1 ------------------------------------------C6_C12_Paraffin : 14.7 C13_C25 : 0 C6_C12_Cycloparaffin : 34.2 C13_C23 : 0 C6_C11_Aromatic : 9.1 C12_C18_Aromatic : 0 C9_C25_Naphtheon : 42.4 Residual : 0 Reference_Temperature : 20 Oil_Viscosity_at_Reference_Temperature : 6.94 Temp_Type : 0 Constant_Oil_Temperature : 25 Section: OIL_PROPERTIES -----------------------------------------Section: SOURCE_NO_2 ------------------------------------------C6_C12_Paraffin : 14.7 C13_C25 : 0 C6_C12_Cycloparaffin : 34.2 C13_C23 : 0 C6_C11_Aromatic : 9.1 C12_C18_Aromatic : 0 C9_C25_Naphtheon : 42.4 Residual : 0 Reference_Temperature : 20 Oil_Viscosity_at_Reference_Temperature : 6.94 Temp_Type : 0 Constant_Oil_Temperature : 25 Section: OIL_PROPERTIES -----------------------------------------Section: SOURCE_NO_3 ------------------------------------------C6_C12_Paraffin : 14.7 C13_C25 : 0 C6_C12_Cycloparaffin : 34.2 C13_C23 : 0 C6_C11_Aromatic : 9.1 C12_C18_Aromatic : 0
134 C9_C25_Naphtheon : 42.4 Residual : 0 Reference_Temperature : 20 Oil_Viscosity_at_Reference_Temperature : 6.94 Temp_Type : 0 Constant_Oil_Temperature : 25 Section: OIL_PROPERTIES -----------------------------------------Section: SOURCE_NO_4 ------------------------------------------C6_C12_Paraffin : 14.7 C13_C25 : 0 C6_C12_Cycloparaffin : 34.2 C13_C23 : 0 C6_C11_Aromatic : 9.1 C12_C18_Aromatic : 0 C9_C25_Naphtheon : 42.4 Residual : 0 Reference_Temperature : 20 Oil_Viscosity_at_Reference_Temperature : 6.94 Temp_Type : 0 Constant_Oil_Temperature : 25 Section: OIL_PROPERTIES -----------------------------------------Section: SOURCE_NO_5 ------------------------------------------C6_C12_Paraffin : 14.7 C13_C25 : 0 C6_C12_Cycloparaffin : 34.2 C13_C23 : 0 C6_C11_Aromatic : 9.1 C12_C18_Aromatic : 0 C9_C25_Naphtheon : 42.4 Residual : 0 Reference_Temperature : 20 Oil_Viscosity_at_Reference_Temperature : 6.94 Temp_Type : 0 Constant_Oil_Temperature : 25 Section: OUTPUT_SPECIFICATION -----------------------------------Number_Of_Output_Areas : 1 Section: OUTPUT_AREA_1 ----------------------------------------File_Name : G:\Rizka Cilacap\SC01\oil\diesel barat_I.dfs2 Title : Diesel Barat I Reference_Area : 1 XRange : 0 884 1 YRange : 40 800 1 TRange : 720 103680 60 Section: OUTPUT_INCLUDE -------------------------------------Output_As : 0 Total_Concentration : True Oil_Fraction_No_1 : False Oil_Fraction_No_2 : False Oil_Fraction_No_3 : False Oil_Fraction_No_4 : False Oil_Fraction_No_5 : False Oil_Fraction_No_6 : False Oil_Fraction_No_7 : False
135 Oil_Fraction_No_8 : False Emulsification_Rate : True Oil_Evaporation : True Dissolution : True Vertical_Dispersed_Oil : True Exceeding_Concentration : True Time_Exposition : True U_Current_Velocity : True V_Current_Velocity : True Section: BASIC_PARAMETERS -----------------------------------------Section: DISCHARGE_LINES ----------------------------------------Number_Of_Discharge_lines : 0 Working on: G:\Rizka Cilacap\SC01\HD\HD-Barat-Oil.dfs2 --------------------------- Data File Summary --------------------------Data origin: C:\Program Files\DHI\MIKE Title : HD-Barat Delete : -1.000E-30 Axis Sets Interval Axis origin Unit time 2585 900 2007-02-01 02:00:00 1 885 10 0 meter 2 935 10 0 meter Item name Item type Minimum Maximum Missing Unit H Water Depth m Water Level 0.925467 25.663668******** meter P Flux m^3/s/m Flow Flux -21.332455 14.248095******** m^3/s/m Q Flux m^3/s/m Flow Flux -12.573159 37.359356******** m^3/s/m ---------------------------------------------------------------------Timestep: 0 ( 50 ) Timestep: 1 ( 100 ) Timestep: 2 ( 150 ) Timestep: 3 ( 200 ) Timestep: 4 ( 250 ) Timestep: 5 ( 300 )
dst. hingga Timestep:
86082 (
1783056 )
136 Lampiran 3. Contoh Sumber Data Mentah Minyak Jenis Crude Oil
RIWAYAT HIDUP Penulis dilahirkan di Jakarta, 16 April 1986 dari pasangan Bapak (Alm.) Hamidi dan Ibu Nani Hendani. Lulus dari SMA Negeri 1 Depok pada tahun 2004, penulis langsung melanjutkan studi di Institut Pertanian Bogor pada Program Studi Ilmu Kelautan, Departemen Ilmu dan Teknologi Kelautan, Fakultas Perikanan dan Ilmu Kelautan melalui jalur USMI. Selama kuliah di Institut Pertanian Bogor, penulis aktif dalam organisasi kemahasiswaan yaitu HIMITEKA sebagai Staf Divisi Kewirausahaan periode 2005-2006. Selain itu, penulis juga aktif menjadi Panitia Musyawarah Nasional HIMITEKINDO pada tahun 2006. Untuk menyelesaikan studi dan memperoleh gelar Sarjana Perikanan di Fakultas Perikanan dan Ilmu Kelautan, penulis menyusun skripsi dengan judul “Model Sebaran Tumpahan Minyak di Alur Pelayaran Pelabuhan Tanjung Intan Cilacap, Jawa Tengah”.
137