## MODEL DISPERSI BAHANG HASIL BUANGAN AIR PROSES PENDINGINAN PLTGU CILEGON CCPP KE PERAIRAN PANTAI MARGASARI DI SISI BARAT TELUK BANTEN

## THERMAL DISPERSION MODEL OF WATER COOLING PLTGU CILEGON CCPP DISCHARGE INTO MARGASARI COASTAL WATERS AT THE WESTERN COAST OF BANTEN BAY
I Wayan Nurjaya P 1) P dan Heron Surbakti P 2) P P 1) P Bagian Oseanografi, Dept. Ilmu dan Teknologi Kelautan, FPIK-IPB
P 2) P Jurusan Ilmu Kelautan, FMIPA-UNRI

## ABSTRACT
Thermal dispersion model based on the hydrodynamics model was applied on PLTGU Cilegon (electric power industry based on gas and steam) at the coast of Margasari, Pulo Ampel District, Serang-Banten.  This PLTGU used around 60.000 m P 3 P /hour of seawater as cooling water system. Therefore, it produced water with high temperature of about 5 P o P C higher than the sourounding of seawater temperature. This high water temperature was flowed out into the coastal waters. This study tried to predict their distribution according to southeast and northeast monsoon. Model verification was conducted both to hydrodynamics component (tide and current) and water temperature. The verification results show good enough patterns between the model results and field measurement.
Key words: Bay of Banten, cooling water system, thermal dispersion, hydrodynamics model

## ABSTRAK
Model dispersi bahang yang didasari dengan model hidrodinamika diaplikasikan di perairan pantai timur Desa Margasari, Kecamatan Pulo Ampel, Kabupaten Serang, Provinsi Banten. Kegiatan PLTGU Cilegon menggunakan air laut dengan debit 60.000 m P 3 P /jam untuk mendinginkan mesinnya ( cooling water system ). Hasil proses pendinginan ini menghasilkan air buangan bersuhu tinggi, selisih nilai suhu air buangan adalah 5 P o P C lebih tinggi dari suhu lingkungan yang disalurkan ke perairan pantai melalui outlet . Air buangan bersuhu tinggi tersebut secara kontinyu dialirkan ke pantai dengan debit 60.000 m P 3 P /jam. Dalam model sebaran bahang dilakukan verifikasi hasil model dengan hasil pengukuran di lapangan, baik verifikasi terhadap komponen hidrodinamika (pasang surut dan arus), maupun verifikasi terhadap komponen suhu. Hasil verifikasi menunjukkan ada kesamaan pola yang cukup baik antara hasil model dan pengukuran di lapangan.
Kata K unci : Teluk Banten, cooling water system , dispersi bahang, model hidrodinamika

## I. PENDAHULUAN
Kebutuhan akan energi listrik belakangan ini meningkat. Peningkatan tersebut disebakan karena permintaan energi listrik untuk pemukiman dan isndustri meningkat. Kegiatan industri pembangkit listrik tenaga gas dan uap (PLTGU) di Indonesia lebih banyak dibangun di dekat pantai misalnya  PLTU
Suralaya di Merak Banten, PLTU Paiton di Paiton Jawa Timur, PLTU Muara Karang di Muara Karang Jakarta. Alasan utamanya adalah lebih mudah memperoleh air yang digunakan untuk proses pendinginan mesin atau cooling water system .
Permasalahan utama yang selalu muncul dari kegiatan buangan air hasil
proses cooling water system adalah suhu air buangan dari sistem tersebut jauh lebih tinggi dari suhu lingkungan di sekitarnya, umunya suhu air buangan tersebut dapat mencapai 40 P o P C. Namun, karena ada aturan bahwa nilai maksimum perbedaan suhu air buangan dengan suhu alami adalah 5 P o P C maka umumnya air buangan tadi dikelola dulu dengan cara menurunkan suhunya dari 40 P o P C menjadi 34 P o P C. PT.Gunanusa Utama Fabricators adalah salah satu perusahan galangan kapal dan pembuat offshore platform atau anjungan minyak lepas pantai yang beroperasi di atas lahan pantai seluas 18 hektar di depan areal PLTGU Cilegon. Adanya peningkatan permintaan atau pemesanan fasilitas penunjang produksi minyak dan gas bumi, maka PT Gunanusa memerlukan lahan yang lebih luas menjadi 24 hektar.  Salah satu upaya yang dilakukannya adalah melalui reklamasi pantai. Oleh karena kegiatan reklamasi pantai ini lokasinya berdampingan dengan lokasi intake dan outlet milik PLTGU, maka diperkirakan akan mengganggu kebutuhan kuantitas dan kualitas air pendingin akibat adanya peningkatan kekeruhan dan pendangkalan di sekitar intake dan outlet . Proses fisik yang paling mendasar dalam transpor bahang ( heat transport ) adalah ketika limbah bahang masuk ke dalam badan air, hal ini menyebabkan suhu air meningkat sampai terjadi kehilangan keseimbangan bahang di permukaan. Untuk keperluan komputasi, masuknya bahang ke dalam badan air dikelompokkan menjadi dua zona, yakni badan air yang dekat dengan sumber buangan ( near-field ) dan jauh dari sumber buangan ( far-field ). Pada zona
pertama, buangan bahang ( heated discharge ) diencerkan oleh adanya turbulen ( discharge-induced turbulence ). Pada zona kedua yang luas dan berdekatan dengan zona pertama, distribusi bahang diatur oleh proses konveksi dengan adanya arus ( convection by ambient currents ), difusi karena adanya turbulen dan adanya pertukaran bahang melalui permukaan laut.
Penelitian  ini bertujuan untuk mengetahui seberapa jauh pengaruh kegiatan reklamasi terhadap lingkungan sekitar, khususnya disperse bahang hasil buangan proses cooling water sehingga hasilnya dapat dipakai sebagai pertimbangan dalam rencana reklamasi ke depan. Karena kegiatan reklamasi belum dilakukan, maka dalam penelitian ini dilakukan beberapa skenario (sebelum reklamasi dan setelah reklamasi) dengan menggunakan model hidrodinamika dan transpor polutan (bahang).

## II. METODE PENELITIAN
Penelitian ini dilakukan sejak awal Januari 2009 sampai awal Februari 2009. Data yang digunakan dalam penelitian ini adalah data primer hasil pengukuran di lapangan (arus, pasut dan suhu air laut), data  sekunder baik yang diperoleh dari internet maupun dari beberapa publikasi ilmiah dan hasil penelitian serta data hasil monitoring di perairan Teluk Banten.
Lokasi penelitian terletak di sisi barat perairan Teluk Banten, persisnya di pantai Desa Margasari,  Kecamatan Pulo Ampel, Kabupaten Serang, Propinsi Banten. Secara geografis lokasi penelitian digambarkan pada Gambar 1.
Gambar 1. Lokasi Penelitian

## 2.1.  Pengumpulan Data Primer
Pengukuran arus dan pasang surut dilakukan selama tiga hari tiga malam pada posisi 5 P o P 55,4934’ L.S dan 106 P o P 07,4975’ B.T. Alat ukur pasang-surut dan arus dipasang di bagan pada kedalaman laut 14 m dengan system Mooring. Tide gauge dengan SBE Type dan Aandera
Current meter tipe RCM 7 diset dengan interval waktu perekaman setiap 10 menit selama selang waktu pengukuran.
Alat yang digunakan mengukur pasang-surut, arus (komponen kecepatan dan arah) serta posisi di tengah laut berturut-turut adalah tide gauge , current meter dan GPS = Global Positioning System (Gambar 2).
(a) (b) (c)
Gambar 2. Peralatan yang digunakan dalam survei lapang (a) GPS Garmin, (b) Tide
Gauge SBE Type , dan (c) Aandera Current Meter tipe RCM 7

## 2.2.  Model

## 2.2.1.  Model Hidrodinakia (Pola Arus)
Model pola arus dibangun dengan menggunakan Piranti lunak MIKE 21 Versi 2005 yang dikembangkan oleh Danish Hydraulics Institute ( DHI ) Water and Environment , Denmark. Persamaan yang digunakan pada model ini adalah persamaan kontinuitas dan persamaan momentum dengan perata-rataan kedalaman.  Model ini menggunakan pendekatan metode beda hingga ( finite difference method ) untuk menyelesaikan persamaan yang digunakan.  Berikut uraian persamaan yang digunakan dalam model diatas (Abbott et al ., 1978):
2.2.2.  Persamaan kontinuitas:
2.2.3.  Persamaan momentum U pada sumbu x ………… (1) U pada sumbu y …………(2)
Dimana: :   kedalaman air ( )
:   kedalaman  air bervariasi terhadap waktu (m)
:  elevasi muka laut (m)
:  kerapatan densitas arah
sumbu x dan y (m P 3 P /s/m)
= ; =  kecepatan rerata terhadap kedalaman pada sumbu x dan y : Chezy resistance (m P 1/2 P /s) :  percepatan gravitasi (m P 2 P /s)
:  faktor gesekan angin
:  kecepatan angin pada sumbu x dan y (m/s) :  parameter Coriolis (s P -1 P )
:  tekanan atmosfer (kg/m/s P 2 P ) :  densitas air laut (kg/m P 3 P )
: komponen shear stress efektif
Bed shear stress pada arah x dan y dapat dihitung dengan (Melor, 1998):
……….. (3)
……….  (4) dimana c R f R adalah koefisien gesekan dan dapat dihitung sebagai berikut
…………….. (5)
dimana C : Koefisien Chezy n : Koefisien Manning
λ : 1,486 untuk unit Inggris dan 1,0 untuk SI
Kedalaman rerata gesekan turbulen dapat dihitung menggunakan konsep viskositas Eddy dari Boussinesq, yaitu :
….………… (6)
1 2 2 2 2 2 1 b b bx f z z c U U V x y τ ρ     ∂ ∂   = + + +       ∂ ∂         1 2 2 2 2 2 1 b b by f z z c V U V x y τ ρ     ∂ ∂   = + + +       ∂ ∂         1 3 2 2 2 f g gn c C H λ = = xx xx U U v x x τ ρ ∂ ∂   = +   ∂ ∂  
…………. (7) ………….. (8) dimana, v adalah 0,3 ± 0,6 U*H V : Kecepatan arus pada arah y U : Kecepatan arus pada arah x

## 2.2.4.  Model Dispersi Panas
Model sebaran panas dibuat dengan menggunakan persamaan Adveksi/Dispersi.   Model adveksi- dispersi diselesaikan dengan persamaan adveksi-dispersi dua dimensi.  Persamaan adveksi-dipersi tersebut dapat dituliskan sebagai berikut (Ramming  and  Kowalik, 1980):
H p C J T T y D A x D
A y T v x T u t T ρ + ∂ ∂ ∂ ∂ + ∂ ∂ ∂ ∂ + ∂ ∆ ∂ − ∂ ∆ ∂ − = ∂ ∆ ∂             ∆ ∆ ) ( ) ( ) ( ) ( ) (
.....................(9) dimana:
:  beda suhu air terhadap suhu ambien
:  kecepatan horizontal pada arah
x dan y (m/s)
: koefisien difusi pada arah x dan y (m P 2 P /s) t :  waktu
J :  Total fluks bahang di permukaan
ρ :  densitas air laut
C R P R : Specific heat H :  kedalaman penyebaran limbah
panas efektif di permukaan laut
Nilai u , v dan h untuk tiap time step diperoleh dari hasil simulasi pola arus. Total fluks bahang di permukaan dalam kasus ini adalah (Chapra, 1997):
J = J R sn R + J R an R – J R br R – J R c R – J R e R ……… (10) J R an R = σ (T+ ∆ T + 273) P 4 P (0,7 + (0,031 Eud P 0.5 P ))0,97 J R br R = ε σ (T+ ∆ T + 273) P 4
J R c R = c R 1 R f(Uw) (T+ ∆ T – Tud) J R e R = f(Uw) (Es – Eud) 133.32
dimana:
J R sn R =  radiasi gelombang pendek (dari matahari) J R an R =  radiasi gelombang panjang (dari atmosfer)
J R br R =  radiasi balik gelombang panjang (dari air)
J R c R =  konduksi/konveksi di permukaan (lapisan udara-air)
J R e R
=  evaporasi/kondensasi (lapisan udara-air)
T =  suhu ambien air laut, ε =  emisifitas air,
σ =  konstanta Stefan-Boltzmann, c R 1 R
=  koefisien Bowen,
f(Uw) =  pengaruh angin dipermukaan,
E R ud R =  tekanan uap air udara,
E R s R =  tekanan permukaan air dan T R ud R =  suhu udara.
Pengaruh angin dipermukaan dinyatakan dalam rumus empirik (Chapra, 1997): f(Uw) = 19 + 0,95 Uw P 2 P ………..  (11) dimana U R w R = kecepatan angin di permukaan.

## 2.2.5.  Desain Model dan Syarat Batas:
Model ini diselesaikan dengan pendekatan numerik yaitu metode beda hingga.  Simulasi model dilakukan untuk kasus arus yang dibangkitkan oleh pasang-surut (pasut) dengan memasukkan juga pengaruh angin pada daerah penelitian.
Syarat batas untuk model hidrodinamika (Gambar 1) adalah:
- syarat batas tertutup (daratan)
( ) 0 , , = η V U ……………..  (12)
- syarat batas terbuka atau open boundary (laut)
Elevasi pasut yang diberikan di kondisi batas adalah berdasarkan hasil peramalan pasut NAO 99 selama 15 hari. Selanjutnya untuk kecepatan arus di batas dan elevasi yang belum diberikan,
xy yx xy U V v y x τ τ ρ   ∂ ∂ = = +   ∂ ∂  
yy yy V V v y y τ ρ   ∂ ∂ = +   ∂ ∂  
menggunakan syarat batas radiasi Orlanski (Kowalik dan Murty, 1993):
0 = ∂ ∂ ± ∂ ∂ x F C t F …………………….  (13)
dimana:
F :  kecepatan arus rata-rata atau elevasi pasut, dan C :  ( gH ) P 0.5 P , yaitu kecepatan gelombang panjang.
Syarat batas untuk model dispersi bahang  di laut mengikuti syarat batas hidrodinamika, dengan mengganti kecepatan dan elevasi menjadi nilai suhu. Pada titik buangan air panas ( outlet ) diberikan sumber  air panas secara kontinu dengan nilai ∆T =5 °C.

## 2.3.  Data Input Model
2.3.1.  Pasang Surut
Data pasang surut yang digunakan dengan menggunakan data peramalan NAO.99b ( National Astronomical Observatory , Jepang).  NAO Tide merupakan suatu model peramalan pasang surut global dengan resolusi 1/2 P o P x 1/2 P o P merupakan data asimilasi dari TOPEX/Poseidon selama 5 tahun.  Data pasut digunakan sebagai batas terluar model ( open boundary condition ).  Pada model digunakan empat batas terluar. Gambaran pasut pada masing-masing batas terluar adalah Gambar 3 untuk musim barat dan Gambar 4 untuk kondisi musim timur.
Data pasang surut dari NAO tide diinterpolasi menjadi tiap 5 detik (sesuai dengan langkah waktu/time step simulasi) dengan menggunakan cubic spline untuk mendapatkan stabilitas model.
Gambar  3.  Elevasi Permukaan Laut pada Batas Terbuka saat Musim Barat (a) Kiri atas Daerah Model (b) Kanan Atas Daerah Model (c) Kiri Bawah Daerah Model (d) Kanan Bawah Daerah Model
Gambar  4.  Elevasi Permukaan Laut pada Batas Terbuka saat Musim Timur (a) Kiri atas Daerah Model (b) Kanan Atas Daerah Model (c) Kiri Bawah Daerah Model (d) Kanan Bawah Daerah Model

## 2.3.2.  Angin
Data angin yang digunakan dengan menggunakan data prediksi CERSAT dengan resolusi 15’.  Data angin yang digunakan sebagai masukan pada model adalah data bulan Januari (musim Barat) dan data bulan Juli (musim Timur).  Pada Lampiran Gambar 1 terlihat angin dominan terjadi pada musim barat adalah angin dari arah Barat Laut dan Utara, sedangkan pada musim  timur, angin dominan bertiup dari arah tenggara.

## 2.3.3.  Batimetri
Sumber peta batimetri yang digunakan berasal dari Dishidros TNI AL
tahun 2008 skala 1 : 200.000 dan peta Lingkungan Pantai Indonesia untuk daerah Teluk Banten dari Bakosurtanal dengan skala 1:50000 (Lembar Peta 1110-09) serta hasil sounding di lokasi kegiatan.  Keseluruhan data kedalaman tersebut diatas digambarkan pada Gambar 1.

## 2.4.  Diskritisasi Daerah Model
Gambaran mengenai model yang digunakan disajikan dalam Tabel 1. Grid yang dibangun pada area model ditunjukkan pada Gambar Lampiran 2.
Tabel 1. Jumlah grid, dimensi grid, time step yang digunakan dalam model
Parameter Value Unit Number Of X (I) Cells 150 - Number Of Y (J) Cells 150 - X Grid Size 20 m Y Grid Size 20 m Time Step Of  Model 5 second Max. Number of Time Steps 518400 second

## III. HASIL DAN PEMBAHASAN

## 3.1.  Verifikasi Model
Verifikasi model arus dan sebaran suhu (bahan) telah dilakukan untuk mengetahui apakah model yang dibangun sesuai atau tidak dengan kondisi di lapangan.  Hasil verifikasi antara data hasil pemodelan dan data hasil pengukuran di lapangan digambarkan pada grafik pencar, fluktuasi muka laut dan suhu berturut-turut pada Gambar 5, 6, dan 7.
Pada Gambar  5 terlihat bahwa pola pergerakan arus hasil pengukuran (tanda dot warna biru) adalah arah barat laut- tenggara.  Pola pergerakan arah arus dari
model hidrodinamika juga memperlihat arah yang sama, namun besarannya antara kedua hasil tersebut terlihat masih ada perbedaan. Data arus hasil pengukuran di ukur pada lapisan kedalaman 5 m, sedangkan data arus hasi model merupakan hasil perata-rataan kedalaman 14 m, sehingga kecepatan arus hasil pengukuran cenderung lebih besar.
Fluktuasi muka laut hasil pengukuran (garis berwarna hitam) selama 3 hari dan tiga malam diflotkan bersama dengan hasil running model (garis berwarna orange) pada saat yang sama menghasilkan Gambar 6.
1.
= Data Hasil Pengukuran                           = Data Hasil Model
Gambar 5.  Hasil verifikasi data arus hasil pengukuran dan model simulasi
Gambar 6.    Hasil verifikasi data pasut hasil pengukuran tinggi paras laut di lapangan dan hasil model simulasi
Pada gambar tersebut diatas terlihat bahwa fluktuasi muka laut dari hasil model dan hasil pengukuran di lapangan terdapat  sedikit perbedaan, meskipun pola naik turunnya muka laut sudah mendekati sama. Perbedaan ketinggian muka laut yang besar terjadi saat pasang dan surut, perbedaannya dapat mecapai 0,3 m. Faktor penyebabnya diperikirakan karena data pasut dari NAO tide yang dipergunakan sebagai kondisi batas pada batas terbuka pada domain model kurang baik diaplikasikan pada perairan dangkal di daerah teluk atau selat dimana konfigurasi pantainya yang rumit.

## 3.2.  Verifikasi Hasil Sebaran Bahang
Verifikasi terhadap data suhu juga dilakukan, dengan cara melakukan pengukuran suhu air laut  di 7 titik sampling yang berbeda, mulai dari titik sumber ( outlet ) kemudian menjauh ke
tengah laut dengan interval jarak 0,5 km. Kedua data suhu tersebut baik dari hasil pengukuran dan hasil model diplotkan pada Gambar 7 untuk melihat seberapa jauh kesamaan dari dua data suhu tersebut.
Data suhu air laut hasil pengukuran di lapangan selalu lebih besar dari data suhu air laut hasil model. Perbedaan ini terutama diprediksi karena data suhu hasil pengukuran hanya mewakili 1 lapisan kedalaman (suhu permukaan laut), sedangkan data suhu dari model adalah hasil dari perata-rataan kedalaman sehingga selalu lebih rendah dari hasil pengukuran. Akan tetapi, pola penurunan suhu hasil buangan cooling water di perairan pantai Desa Margasari, Kecamatan Pulo Ampel, Kabupaten Serang, Propinsi Banten mempunyai kemiripan yang tinggi antara data suhu dari pengukuran di lapangan dan hasil dari model.
Gambar 7. Verifikasi hasil model disperse bahang dan hasil pengukuran pada 7 titik sampling dari sumber buangan (outlet) dan dengan interval jarak 0,5 km menjauhi titik outlet

## 3.3.  Hasil Simulasi Pola Arus
Arus merupakan massa air yang bergerak dapat dibangkitkan oleh berbagai gaya baik secara ekstrenal maupun internal. Pergerakan massa air atau arus ini merupakan media yang mampu memindahkan bahang dari sumbernya ke tempat lain. Untuk mengetahui kecepatan, arah dan pola gerakan arus dapat dilakukan dengan berbagai cara, diantaranya dengan melakukan pengukuran langsung di laut (metode euler atau langrangian), membuat model simulasi menggunakan persamaan matematika dengan bantuan komputer melalui pemahaman karakteristik fisika air laut dan faktor- faktor pembangkitnya.
Pola arus yang terjadi di perairan Teluk Banten dipengaruhi oleh fenomena pasang-surut (pasut) dan angin musim (monsoon) yang terjadi di Indonesia khususnya di Laut Jawa (Wyrtki, 1961; Hosoya and Muchari, 1986; Qu et al ., 2005).  Dalam model simulasi arus ini sangat mempertimbangan kondisi tersebut  agar mendapatkan kondisi yang sesungguhnya di perairan pantai Desa Margasari, Pulo Ampel, Serang-Banten.
Gambar 8 sampai dengan Gambar 11 merupakan pola arus hasil simulasi
pada musim barat pada empat posisi muka air yang berbeda baik pada saat pra dan pasca reklamasi, yakni saat muka laut menuju surut minimum, surut minimum, muka laut laut menunju pasang maksimum dan saat pasang maksimum.  Perbedaan pola arus saat menuju surut minimum dan pasang maksimum sangat jelas. Saat muka laut menuju surut minimum massa air bergerak meninggalkan perairan pantai Desa Margasari.  Sebaliknya, pada saat muka laut menuju pasang maksimum massa air bergerak mendekati pantai, sebagian bergerak kearah barat laut sejajar pantai dan sebagian bergerak kearah selatan dengan kecepatan arus berkisar antara 0,1 -0,9 m/det.
Pada musim timur juga dilakukan simulasi yang sama, pada empat posisi muka laut yang berbeda. Gambar 12 sampai dengan Gambar 15 merupakan hasil simulasi arus pada musim timur, dimana pola arus secara umum terlihat sama. Pada saat air laut menuju surut maksimum massa air terlihat bergerak menjauhi pantai, sebaliknya pada saat muka laut menuju pasang maksimum massa air bergerak mendekati pantai. Ketika mendekati pantai massa air bergerak ke barat laut sejajar pantai, sebagian lainnya bergerak ke selatan.
Gambar  8.   Pola arus hasil simulasi saat MSL menuju surut minimum pada saat musim barat (Panel kiri atas sebelum reklamasi, Panel Kanan atas sesudah reklamasi, panel kiri bawah posisi muka laut, dan panel kanan  bawah adalah kecepatan dan arah angin)
Gambar 9.   Pola arus hasil simulasi saat surut minimum pada saat musim barat (Panel kiri atas sebelum reklamasi, Panel Kanan atas sesudah reklamasi, panel kiri bawah posisi muka laut, dan panel kanan bawah adalah kecepatan dan arah angin)
Gambar 10.    Pola arus hasil simulasi saat MSL menuju pasang maksimum pada saat musim barat (Panel kiri atas sebelum reklamasi, Panel Kanan atas sesudah reklamasi, panel kiri bawah posisi muka laut, dan panel kanan bawah adalah kecepatan dan arah angin)
Gambar 11.    Pola arus hasil simulasi saat pasang maksimum pada saat musim barat (Panel kiri atas sebelum reklamasi, Panel Kanan atas sesudah reklamasi, panel kiri bawah posisi muka laut, dan panel kanan bawah adalah kecepatan dan arah angin
Gambar 12.  Pola arus hasilsimulasi saat MSL menuju surut minimum pada saat musim timur (Panel kiri atas sebelum reklamasi, Panel Kanan atas sesudah reklamasi, panel kiri bawah posisi muka laut, dan panel kanan bawah adalah kecepatan dan arah angin)
Gambar 13.   Pola arus hasil simulasi saat surut minimum pada saat musim timur (Panel kiri atas sebelum reklamasi, Panel Kanan atas sesudah reklamasi, panel kiri bawah posisi muka laut, dan panel kanan bawah adalah kecepatan dan arah angin)
Gambar 14.   Pola arus hasil simulasi saat MSL menuju pasang maksimum pada saat musim timur (Panel kiri atas sebelum reklamasi, Panel Kanan atas sesudah reklamasi, panel kiri bawah posisi muka laut, dan panel kanan bawah adalah kecepatan dan arah angin)
Gambar 15.  Pola arus hasil simulasi saat pasang maksimum pada saat musim timur (Panel kiri atas sebelum reklamasi, Panel Kanan atas sesudah reklamasi, panel kiri bawah posisi muka laut, dan panel kanan bawah adalah kecepatan dan arah angin)
Tabel  2. Kecepatan arus hasil simulasi untuk musim barat dan timur pada empat posisi muka laut berbeda baik sebelum dan setelah reklamasi MUSIM Posisi Muka Laut Kecepatan Arus (m/s) Sebelum Reklamasi Setelah Reklamasi BARAT  MSL ST Maksimum 0.958 0.961 Rata-rata 0.364 0.366 ST Maksimum 0.238 0.238 Rata-rata 0.067 0.068 MSL PT Maksimum 0.437 0.437 Rata-rata 0.211 0.211 PT Maksimum 0.288 0.288 Rata-rata 0.122 0.122 TIMUR  MSL ST Maksimum 1.003 1.006 Rata-rata 0.371 0.374 ST Maksimum 0.243 0.243 Rata-rata 0.073 0.073 MSL PT Maksimum 0.466 0.466 Rata-rata 0.224 0.225 PT Maksimum 0.331 0.331 Rata-rata 0.144 0.145
U Keterangan U :
MSL ST= muka laut di Mean Sea Level menuju titik surut minimum
ST = muka laut di titik surut minimum MSL PT= muka laut di Mean Sea Level menuju titik pasang maksimum
PT = muka laut di titik pasang maksimum
Kecepatan arus pada saat posisi muka laut di titik Mean Sea Level menuju Surut terendah (MSL ST), di titik surut terendah (ST), Mean Sea Level menuju pasang tertinggi (MSL PT) dan pasang tertinggi (PT) baik musim barat maupun musim timur serta kondisi sebelum reklamasi dan stelah reklamasi ditabulasikan pada Tabel 2 diatas. Secara umum terlihat bahwa adanya kegiatan reklamasi terutama saat muka laut menuju titik surut minimum pada musim timur dapat meningkatkan kecepatan arus sebesar 0,003 m/det. Pada posisi muka laut lainnya di kedua musim tidak terjadi perubahan kecepatan arus.

## 3.4.  Hasil Simulasi Dispersi Bahang
Gambar  16 dan Gambar 17 merupakan beberapa contoh hasil dispersi bahang setelah di- running selama 15 hari, baik pada musim barat dan musim timur serta sebelum dan setelah reklamasi. Setiap Gambar disajikan dengan 4 panel, panel kiri atas menunjukkan hasil dispersi bahang sebelum dilakukan reklamasi pantai, panel kanan atas menggambarkan dispersi bahang setelah reklamasi. Pada panel kiri bawah menggambarkan pola arus sebelum ada reklamasi dan panel kanan bawah setelah reklamasi. Skala warna yang ada disisi kanan pada kedua panel diatas  menunjukkan skala suhu,
nilai suhu meningkat kearah warna merah, sebaliknya suhu menurun bila warna biru. Skala warna yang ada disisi kanan pada panel-panel dibawah menunjukkan kedalaman laut yang berubah-ubah sesuai dengan tinggi muka laut  diakibatkan oleh fenomena pasang surut. Kedalaman laut bertambah tentu pada periode pasang ditandai dengan warna biru, sebaliknya kedalaman laut berkurang pada periode surut, ditandai dengan warna biru muda.
Dispersi bahang yang keluar pada outlet pembuangan air hasil proses cooling water system ditandai dengan warna merah kira-kira 150 m sebelah selatan batas reklamasi, baik pada
kondisi pra dan pasca-reklamasi sebarannya terlihat bergerak kearah selatan menyusur garis  pantai. Garis isothermal 0,5 P o P C terlihat memanjang dan bergerak ke selatan sejauh 600 meter pada periode tidak ada reklamasi. Bertambah ke selatan kira-kira sampai 750 m dari titik outlet pada periode pasca reklamasi. Sebaran thermal dari garis pantai hamper sama antara sebelum dan setelah reklamasi.
Kondisi yang hampir sama terlihat pada skenario musim timur (Gambar 17),yang berbeda adalah jarak maksimum sebaran isothermal 0,5 P o P C lebih pendek bila dibandigkan dengan musim barat.
Gambar 16.  Sebaran panas (delta suhu 5 P o P C) terjauh hasil simulasi pada saat musim barat sebelum dan setelah reklamasi.
Gambar 17.  Sebaran panas (delta suhu 5 P o P C) terjauh hasil simulasi pada saat musim timur sebelum dan setelah reklamasi.

## IV.  KESIMPULAN
Hasil model simulasi arus mempunyai pola yang sama yakni bergerak dengan sumbu utama barat laut dan tenggara, namun besaran kecepatan arus dari model lebih kecil bila dibandingkan dengan hasil pengukuran arus di lapangan.
Pola dispersi bahang hasil buangan dari proses pendinginan turbin PLTGU yang di- running selama 15 hari baik pada saat musim timur dan musim barat hampir sama, yaitu 600-750 ke arah selatan dengan beda suhu 0,5 ⁰ C.
Data pasut dari NAO tide kurang baik digunakan di perairan pantai yang dangkal dan memiliki konfigurasi yang kompleks, seperti di perairan sisi barat Teluk Banten.

## SARAN
Untuk verifikasi model, data arus dan suhu hasil pengukuran di lapangan sebaiknya diukur pada setiap lapisan kedalaman, kemudian dirata-ratakan berdasarkan kedalaman. Ketika diverifikasi dengan data dari model hasilnya lebih baik.
Untuk input data pada batas terbuka di perairan dangkal dan mempunyai konfigurasi pantai yang komplek, sebaiknya tidak menggunakan data pasut dari NAO.

## Ucapan Terima Kasih
Penulis mengucapkan terima kasih kepada PT. Gunanusa dan LPPM IPB atas kesempatan yang diberikan dalam studi kajian sebaran bahang hasil
buangan air pendingin PLTGU Cilegon di Serang, Provinsi Banten.

## Daftar Pustaka
Abbott, M.B., H.M. Petersen and O. Skovgard. 1978. On the Numerical Modelling of Short Waves in Shallow Water. Journal of Hydraulic Research, 16(3).
Chapra, 1997.  Surface water quality modelling. McGraw-Hill. Singapore.
Hydrodynamics Module of MIKE 21 Flow Model, 2007.   Scientific Documentation, DHI Water Environmental. Denmark. 58 p. Hosoya, M. and M. Muchari. 1986.
Survey on the oceanographical conditions in Banten Bay October 1979 to September 1980. In Scientific Report of Mariculture Research and Development Project (ATA-192) in Indonesia March 1986. JICA Jakarta. 289- 321 p.
Kowalik, Z. dan T.S. Murty. 1993. Numerical modelling of ocean dynamics. Advance series on ocean engineering, vol. 5. World Scientific Publishing Co. Pte. Ltd. Singapore Melor, G. L. 1998. User guide for a three dimensional primitive equation,
numerical ocean model. Princeton
University.
MIKE. 2005Coastal Hydraulics and Oceanography, Hydrodynamic
Module, Scientific
Documentation, DHI Software.
MIKE. 2005.  Environmental Hydraulics Advection, Mud Transport Module, User Guide, DHI Software.
MIKE. 205.  Environmental Hydraulics Advection, Mud Transport Module, Scientific Background,
DHI Software.
MIKE 21 Environmental Hydraulics Advection, Dispersion Module, Reference Manual, DHI Software.
MIKE. 2005.   Environmental Hydraulics Advection,  Dispersion Module, scientific Documentation, DHI Software.
Qu, T., Yan D., J. Strachan, G. Meyers and J. Slingo. 2005. Sea  surface temperature and its variability in the Indonesian region. Oceanography , 18(4):50-61.
Ramming, H-G. dan Z. Kowalik.  1980. Numerical Modelling of Marine Hydrodynamics; Applications to Dynamic Physical Processes. Elsevier Scientific Publishing Company. Amsterdam Wyrtki, K. 1961. Physical Oceanography of the Southeast Asian Waters. Naga Report. Vol. 2. The University of California Scripps Institution of Oceanography, La Jolla, California.
Lampiran Gambar 1. Windrose diolah dari data angin selama 10 tahun terakhir
Lampiran Gambar 2. Grid daerah model sebelum dan setelah kegiatan reklamasi