Model Kualitas Air Sungai
Model Kualitas Air Sungai
Model Kualitas Air Sungai
ALI MASDUQI
MK PEMODELAN LINGKUNGAN
PROGRAM STUDI ILMU DAN TEKNOLOGI LINGKUNGAN
UNIVERSITAS AIRLANGGA
2013
DASAR PERHITUNGAN MODEL
KUALITAS BADAN AIR
PEMODELAN KUALITAS AIR
Beberapa Model Kualitas Badan Air
Input
V1;C1
Output
Sistem Aliran
VT;CT
B
2
Berdasarkan konsep mass balance, maka dapat dirumuskan sebagai berikut:
Q1C1 + QACA + QBCB = Q2C2 atau Q C QAC A QB C B
C2 1 1
di mana: Q2
Q2 = Q1 + QA + QB
Model Perhitungan Streeter-Phelps
Titik masuknya pencemar
DOsat
Dissolve
Oksigen (DO)
Point source 4
5 Non-point
6 abstraction
Non-point
7
source
8 Point source
Downstream boundary
Segmentasi
Proses alamiah (Fisik/Kimia/Biologi) sangat dinamis dan
kompleks diperlukan pembagian segmen untuk
mempermudah pendefinisian sifat-sifat proses
Satu segmen sebisa mungkin memiliki karakter yang
mirip/hampir sama diasumsikan sama
Setiap segmen seharusnya cukup representatif diwakili
oleh nilai rata-rata data segmen tersebut, misalnya kedalaman
rata-rata, lebar rata-rata, slope rata-rata (beda elevasi), debit
Setiap segmen jelas point source dan point abstraction,
non point source dapat didekati secara matematis dengan
flow balance dan sampling kualitas
Mass Balance
atmospheric
transfer
mass load mass abstraction
inflow outflow
i
dispersion dispersion
Qi1 Qi
25%Qnpt 25%Qnpt 50%Qnpt
i1 i i+1
Qi Qi 1 Qin,i Qab,i
start end
Bi
Hh
Hi Hd
Hi
Hw Hw
Hi+1
elev2i elev2i
elev1i+1 elev1i+1
Qi 1.83Bi H h3 / 2
Kecepatan aliran di Weir
Qi
Ac ,i Bi H i Ui
Ac,i
Manning Equation
S0
1 H 1
ss1 ss2
Q, U
B0
C h u r c h ill
0.2
0.5
1
1 2
5
10
20
50
100 Owens
Gibbs
0.1
0.1 1
Velocity (mps)
Model Temperatur
atmospheric
transfer
heat load heat abstraction
inflow outflow
i
dispersion dispersion
sediment-water
transfer
sediment
Surface Heat Flux
radiation terms non-radiation terms
air-water
interface
cT o x s
h n dn s
rnd no na nn Alk
s se se mi
s
d
rnx
rcn
r
h p ab ap
rpd po pi rcp s
s s se
rpx cT o
Reaksi-reaksi kimia fundamental
• Fotosintesis dan Respirasi P
106CO2 16NH 4 HPO 24 108H 2 O C106 H 263O110 N16 P1 107O 2 14H
R
P
106CO 2 16 NO 3 HPO 24 122H 2 O + 18H +
C106 H 263 O110 N16 P1 138O 2
R
• Nitrifikasi
NH 4 2O 2 NO 3 H 2 O 2H
• Denitrifikasi
5CH 2 O 4NO 3 4H 5CO 2 2N 2 7H 2 O
• Reaksi degradasi organik aerobik dan anaerobik,
hidrolisis,dll
Beberapa persamaan yang digunakan dalam
perhitungan model (mass balance)
Dissolve Oxygen
S o roa PhytoPhoto ro d BotAlgPhot o roc FastCOxid roc SlowCOxid ron NH4Nitr
roa PhytoResp rod BotAlgResp OxReaer CODoxid
CBOD
S cs rod DetrDiss SlowCHydr SlowCOxid
Organic Nitrogen
S no rna PhytoDeath q N BotAlgDeath ONHydr ONSettl
Organic Phospor
S po rpa PhytoDeath q P BotAlgDeath OPHydr OPSettl
CH4(gas)
POC
ANAEROBIC
Debit Izin
Lokasi
No Nama Perusahaan Pembuangan
Outlet
Limbah (m3/hr)
1 ASX, PT Km 40 18000
2 KSX, PT Km 35 45
3 KD, PT Km 30 150
4 MX, PT Km 28 10000
5 XX, PT Km 15 50
6 SAX, PT Km 10 30
7 ASX, PT Km 5 8500
8 YYY, PT Km 4 15000
Data Meteorologi
• Angin rata-rata (arah dan kecepatan)
• Kelembaban udara
• Lama penyinaran dan intensitas matahari
• Dewpoint temperature (temperatur
pengembunan)
• Prosentase mendung (clode cover)
• Zona waktu (GMT)
Data-data kinetika kimia (opsional)
• Data laju reaksi dan koefisien reaksi tiap segmen
– sLaju aerasi dan reaerasi
– Laju nitrifikasi dan denitrifikasi
– Laju oksidasi organik
– Laju hidrolisis
– Respirasi dan fotosintesis, dll
• Data nilai parameter kualitas badan air
– pH, warna, SS, organik (COD/BOD), N-org, N-anorg,
P-org, DO
Sheet Excel untuk Entri Data (Qual2kw)
Head water sheet
Reach sheet
Reach rates sheet Digunakan
Dew point temperatur sheet untuk
Air temperature sheet
membangun
Wind speed
Model
Cloude cover
Shade
Light and heat
Point source
Diffuse Source
Sheet Excel untuk Entri Data (Qual2kw)
Data 5 tahun di titik-titik
pemantauan
• Data hidrolik
• Data temperatur Digunakan
untuk
• Water Quality (WQ) mencocokkan
Data Model dengan
• WQ Data min hasil monitoring
• WQ Data max
• Diel Data
Contoh Entri Data
Pada Sheet Excel
ENTRI DATA (QUAL2KW)
System ID:
River name Kali Simulasi
Saved file name Sim-simulasi-1
Directory where the input/output files are saved G:\PROJECT\Pelatihan\Simulasi
Month 6
Day 23
Year 2010
Time zone Pacific
Daylight savings time No
Simulation and output options:
Calculation step 6minutes
Number of days 2days
Solution method (integration) Runge-Kutta
Solution method (pH) Newton-Raphson
Simulate hyporheic exchange and pore water quality Level 1
Display dynamic diel output No
State variables for simulation All
Simulate sediment diagenesis No
Program determined calc step 5,625minutes
Time elapsed during last model run 1,25minutes
Time of sunrise 5:24 AM
Time of solar noon 11:26 AM
Time of sunset 5:28 PM
Photoperiod 12,06hours
ENTRI DATA (KUALITAS HEADWATER)
Headwater Flow 21,880m3/s
Prescribed downstream boundary? No
Headwater Water Quality Units 12:00 AM 1:00 AM
Temperature C 28,95 28,95
Conductivity umhos 400,00 400,00
Inorganic Solids mgD/L 0,00 0,00
Dissolved Oxygen mg/L 6,22 6,22
CBODslow mgO2/L 1,00 1,00
CBODfast mgO2/L 2,00 2,00
Organic Nitrogen ugN/L 0,00 0,00
NH4-Nitrogen ugN/L 100,00 100,00
NO3-Nitrogen ugN/L 0,00 0,00
Organic Phosphorus ugP/L 0,00 0,00
Inorganic Phosphorus (SRP) ugP/L 50,00 50,00
Phytoplankton ugA/L 0,00 0,00
Detritus (POM) mgD/L 0,00 0,00
Pathogen cfu/100 mL 10,00 10,00
Generic constituent user defined 4,00 4,00
Alkalinity mgCaCO3/L 100,00 100,00
pH s.u. 7,85 7,85
REACH
Reach for
diel plot: 5
Downstrea
m Elevation Downstream
Upstrea Downstrea
Reach Downstream location m m Latitude Longitude
Degree Minute Second Degree Minute Second
Label end of reach label (km) (m) (m) s s s s s s
8/10/18 52
Hasil Simulasi skenario 2
8/10/18 53
Hasil Simulasi skenario 3
8/10/18 54
Hasil Simulasi skenario 4
8/10/18 55
Penentuan Daya Tampung
• Penentuan daya tampung ini memanfaatkan hasil
simulasi tanpa ada sumber pencemar (misal skenario
3) dan simulasi air sungai memenuhi baku mutu (misal
skenario 4) dari worksheet ‘Sources Summary’.
• Besarnya daya tampung dihitung dari selisih inflow
beban pencemaran hasil simulasi tanpa ada sumber
pencemar dan hasil simulasi air sungai memenuhi
baku mutu (Tabel 1).
8/10/18 56
Tabel 1:
Contoh Hasil Penentuan Daya Tampung
BOD COD
Reach KM
kg/hari kg/hari
Reach 1 (Mlirip) 42,3-36,9 -972,86 -4993,92
Reach 2 36,9-31,6 -3284,93 -16467,84
Reach 3 31,6-21,7 197,86 660,96
Reach 4 21,7-12,0 1477,44 4665,60
Reach 5 12,0-5,6 1280,88 3548,88
Reach 6 5-6-2,6 548,64 483,84
Reach 7 (Jagir) 2,6-0,0 6480,00 9590,40
Jumlah 5727,02 -2512,08
Contoh pembagian reach
Pemanfaatan Hasil Penentuan Daya Tampung (1)
Nilai daya tampung (Tabel 1) setiap reach dibandingkan
dengan nilai beban existing (Tabel 2) .
Bila nilai daya tampung lebih besar daripada beban
existing, maka sungai masih bisa menerima beban
pencemaran tambahan masih boleh membuang
limbah ke sungai
Bila nilai daya tampung lebih kecil daripada beban
existing, maka sungai tidak bisa menerima beban
pencemaran tambahan tidak boleh membuang limbah
ke sungai dan beban pencemaran harus diturunkan.
Pemanfaatan Hasil Penentuan Daya Tampung (2)
• Besarnya persen penurunan beban pencemaran
dapat dihitung dari persentase selisih inflow
beban pencemaran hasil simulasi skenario 1
dan hasil simulasi skenario 4. Simulasi skenario
4 adalah kondisi air Surabaya tidak tercemar.
Contoh hasil perhitungan besarnya persen
penurunan beban pencemaran yang harus
diturunkan agar air sungai tidak tercemar dapat
dilihat pada Tabel 3.
Tabel 2:
Contoh Beban Existing
BOD COD
Reach KM
kg/hari kg/hari
Reach 1 (Mlirip) 42,3-36,9 4371,93 12301,80
Reach 2 36,9-31,6 58299,28 145228,51
Reach 3 31,6-21,7 536,16 1994,90
Reach 4 21,7-12,0 1017,79 2905,25
Reach 5 12,0-5,6 14096,14 45267,43
Reach 6 5-6-2,6 1662,18 3192,08
Reach 7 (Jagir) 2,6-0,0 10849,75 20205,21
Jumlah 90833,24 231095,18
Tabel 3:
Contoh Besarnya Penurunan Beban Pencemaran
8/10/18 63