MODEL

KUALITAS AIR PERMUKAAN

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

• Model Mass Balance / Neraca Masa

• Model Streeter – Phelps
• Model Komputasi
Neraca Massa/Mass Balance
Input
V2;C2

Input
V1;C1
Output
Sistem Aliran
VT;CT

Input Rumus umum pengenceran:

V1C1 + V2C2 + ...... + VnCn = VTCT
Vn;Cn V = Volume
C = Konsentrasi
Model Mass Balance dalam sistem sungai
A

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)

Jarak atau Waktu

Titik masuknya pencemar
Dasar Perhitungan Streeter-Phelps
Persamaan Oxygen Sag: K d Lo
Dt  (e kd t  e kat )  Do e kat
Ka  Kd
di mana:
Dt = defisit oksigen = DOsaturasi – DO
Kd = koefisien deoksigenasi
Ka = koefisien reaerasi
Do = defisit oksigen mula-mula
Lo = BOD ultimate mula-mula
t = time travel
Keterbatasan metoda Streeter-Phelps:
Hanya untuk sumber pencemar tunggal
Terbatas pada parameter BOD dan DO
 Prinsip dalam Penggunaan
Model Komputasi Qual2KW
• Segmentasi Badan Air
• Mass and Flow Balance
• Karakteristik Hidraulik
– Weir, Waterfall, Manning Equation, Rating Curve
– Time Travel
– Longitudinal Dispersion
• Model Temperatur
– Surface heat flux
– Solar radiasi, wave radiasi, konduksi/konveksi
 Segmentasi Badan Air
Headwater boundary
1
Point source
2
Point abstraction
Point abstraction
3

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

bottom algae sediments

 Flow Balance dalam segmentasi
Qin,i Qab,i Qnpt

Qi1 Qi
25%Qnpt 25%Qnpt 50%Qnpt
i1 i i+1

Qi  Qi 1  Qin,i  Qab,i
start end

• Point Source dan non point source

 Karakteristik Hidraulik - WEIR
(a) Side (b) Cross-section

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

1/ 2 5/3 Q = flow [m3/s], S0 = Slope [m/m],

S A
Q 0 c
2/3
n = koefisien kekasaran manning,
n P Ac = luas penampang[m2], and
Q
U P = keliling basah [m]
Ac
 Koefisien Kekasaran Manning
MATERIAL n
Man-made channels
Concrete 0.012
Gravel bottom with sides:
concrete 0.020
mortared stone 0.023
riprap 0.033
Natural stream channels
Clean, straight 0.025-0.04
Clean, winding and some weeds 0.03-0.05
Weeds and pools, winding 0.05
Mountain streams with boulders 0.04-0.10
Heavy brush, timber 0.05-0.20
Reaerasi sebagai fungsi kedalaman
dan kecepatan aliran
O’Connor
Dobbins
10
0.05
0.1
D e p th (m )

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

solar atmospheric water conduction evaporation

shortwave longwave longwave and and
radiation radiation radiation convection condensation

net absorbed radiation water-dependent terms

 Kinetika dan Transfer Massa
Kinetika Reaksi
Pelarutan/dissolution (ds)
Reaksi hydrolysis (h)
Rekasi oxidation (x)
Reaksi nitrifikasi (n), denitrification (dn)
Reaksi photosynthesis (p), death (d), and respiration/excretion (r).
Mass transfer processes
reaeration (re),
settling (s),
sediment oxygen demand (SOD),
sediment exchange (se), and sediment inorganic carbon flux (cf). Note that
the subscript x for the stoichiometric conversions stands for chlorophyll a
(a) and dry weight (d) for phytoplankton and bottom algae, respectively
Model Kinetika dan Transfer Massa
re
rdx
ds cT o
h
mo rcd cs cf cf sod
s se x

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

Total Inorganic Carbon

S cT  rcco FastCOxid  rcca PhytoResp  rccd BotAlgResp

 rcca PhytoPhoto  rccd BotAlgPhoto  CO2Reaer

WATER Reaksi pada sedimen (SOD)
Jpom cf o na o nn pi
AEROBIC

CH4 CO2 NH4p NH4d NO3 N2 PO4p PO4d

CH4(gas)
POC
ANAEROBIC

PON NH4p NH4d NO3 N2

POP PO4p PO4d

DIAGENESIS METHANE AMMONIUM NITRATE PHOSPHATE

 Perhitungan dalam Model
Dari berbagai aspek prinsip (mass balance, flow balance,
karakteristik hidraulik, kinetika reaksi, dan transfer
massa) maka akan diperoleh nilai perubahan konstituen
(parameter/zat) kualitas badan air  didapatkan nilai r
(laju perubahan=dC/dt) dan k (konstanta laju reaksi)
Nilai r dan k  digunakan untuk memodelkan perubahan
kualitas badan air, berdasarkan travel time (waktu reaksi)
Perubahan polutan  sangat dipengaruhi dengan
perubahan oksigen terlarut (DO) dan kinetika reaksi
seluruh konstituen yang berhubungan  kompleks/saling
terkait antar parameter pencemar
 Beberapa nilai r dan k model (contoh)
Parameter Value Units Symbol
Stoichiometry:      
Carbon 49,726 gC gC
Nitrogen 7,2 gN gN
Phosphorus 1 gP gP
Dry weight 100 gD gD
Chlorophyll 1 gA gA
Inorganic suspended solids:      
Settling velocity 0,31902 m/d vi
Oxygen:      
Reaeration model Internal
Temp correction 1,024 qa
Reaeration wind effect None
O2 for carbon oxidation 2,69 gO2/gC roc
O2 for NH4 nitrification 4,57 gO2/gN ron
Oxygen inhib model CBOD oxidation Half saturation
Oxygen inhib parameter CBOD
oxidation 0,60 mgO2/L Ksocf
Oxygen inhib model nitrification Half saturation
Parameter VALUE UNIT simbol
Oxygen enhance model denitrification Half saturation
Oxygen enhance parameter
denitrification 0,60mgO2/L Ksodn
Oxygen inhib model phyto resp Half saturation

Oxygen inhib parameter phyto resp 0,60mgO2/L Ksop

Oxygen enhance model bot alg resp Half saturation
Oxygen enhance parameter bot alg
resp 0,60mgO2/L Ksob
Slow CBOD:      

Hydrolysis rate 0,3/d khc

Temp correction 1,047 q hc

Oxidation rate 0,2/d kdcs

Temp correction 1,047 q dcs

Fast CBOD:      

Oxidation rate 4/d kdc

Temp correction 1,047 q dc

Parameter Value Unit Simbol
Organic N:      
Hydrolysis 0/d khn
Temp correction 1,07 q hn
Settling velocity 0m/d von
Ammonium:      
Nitrification 2/d kna
Temp correction 1,07 q na
Nitrate:      
Denitrification 4/d kdn
Temp correction 1,07 q dn
Sed denitrification transfer coeff 0,36422m/d vdi
Temp correction 1,07 q di
Organic P:      
Hydrolysis 0/d khp
Temp correction 1,07 q hp
Settling velocity 0m/d vop
Inorganic P:      
Settling velocity 6m/d vip
Sed P oxygen attenuation half sat constant 1,40878mgO2/L kspi
Phytoplankton:      
Max Growth rate 1,5/d kgp
Temp correction 1,07 q gp
Respiration rate 0,1/d krp
Temp correction 1,07 q rp
Death rate 0/d kdp
Temp correction 1 q dp
 Dasar Perhitungan Model
• Nilai – nilai r dan k, dapat diubah dan disesuaikan
sehingga model kualitas yang diperoleh menghasilkan
nilai kualitas dan kuantitas air yang mirip (simpangan
kecil) dengan kualitas dan kuantitas di titik-titik
pemantauan pada masing-masing segmen (dari
pemodelan data 4 atau 5 tahun)
• Jika nilai r dan k telah sesuai  dilakukan validasi
dengan data 1 tahun terakhir, pada head water, dan
setiap segment, akan dilihat seberapa besar
simpangan yang terjadi.
– Jika hasil validasi sangat baik  model bisa digunakan
untuk menghitung daya tampung badan air
– Jika hasil validasi tidak baik  model perlu di evaluasi
kembali
ENTRI DATA YANG DIPERLUKAN DALAM
MODEL PERHITUNGAN DAYA TAMPUNG
(QUAL2KW)
Review Model Qual2kw
 Data-data yang diperlukan
Peta Sungai, koordinat posisi, elevasi sungai, profil hidrolis
Data kedalaman dan lebar sungai
Debit head water /input (5 tahun)
Data pemantauan kualitas head water (5 tahun)
Data weir, waterfall, koeffisien manning, slope
Koefisien laju reaksi tiap reach
Data meteorologi: dew point, temperatur, solar radiasi, wind
speed
Point source; posisi, debit, kualitas
Diffuse source/non point source; range posisi, debit, kualitas
Data-data monitoring; data hidrolik (debit, kedalaman, kecepatan,
waktu alir), data temperatur (rata-rata, min, max), data kualitas
(rata-rata, min, max), diel data
 DATA DAN ANALISIS
• Orientasi Lapangan: penyusuran sungai, penentuan koordinat,
elevasi, pointsource, non point source, pengambilan data primer (debit dan
kualitas) dsb
• Kondisi Hidrolis Sungai: data pengukuran penampang sungai
• Debit dan Kualitas Air Sungai: data dari monitoring rutin selama 5
tahun
• Debit dan Kualitas Air Sumber Pencemar: data limbah
industri hasil pemantauan 5 tahun, data anak sungai, data
kualitas limbah domestik, pertanian, dsb
• Penggunaan Air Sungai: data dari instansi terkait, ijin
pemakaian air, intake PDAM, irigasi (anak sungai), dsb
• Klimatologi: data dari BMG
• Peta Sungai: Bakosurtanal dan Google earth
Peta Sungai dan Pembagian Reach/Segmen
DATA HIDROLIK SEGMEN / REACH
Kedalaman Kecepatan aliran
Lebar Sungai
KM Sungai Rata-rata Rata-rata
Rata-rata (m)
(m) (m/detik)
42,3 - 36,9 3,54 32,22 0,36
36,9 - 31,6 3,05 43,82 0,39
31,6 - 21,7 3,19 39,32 0,41
21,7 – 11,9 2,96 42,22 0,41
11,9 – 5,6 4,31 47,14 0,22
5,6 - 2,6 3,95 51,18 0,18
2,6 - 0 3,66 52,73 0,20

Debit air Debit air

Debit air
Lokasi pada musim pada musim
Rata-rata
pengukuran hujan kemarau
(m3/det)
(m3/det) (m3/det)
Dam 1 29,2 24,18 21,2
Jemb 1 63,2 48.34 33,5
Dam 2 50,2 32,69 16,2
 Data-data Profil Sungai
• Peta Sungai dan Pembagian Reach/Segmen
• Koordinat Posisi dan elevasi
• Data debit dan kualitas headwater
• Koordinat point-source atau letak point source
(dalam km)
• Range letak non-point source (koordinat atau km)
• Data posisi titik pemantauan/monitoring (koordinat
atau km)
• Data dam, weir, waterfall, dll
 Stasiun Pemantauan
Stasiun Pemantauan Kualitas Air
Kode KM Lokasi
A KM 41,2 Pemantauan 1
B KM 35,3 Pemantauan 2
C KM 24,15 Pemantauan 3
D KM 16,65 Pemantauan 4
E KM 11,75 Pemantauan 5
F KM 9,25 Pemantauan 6
G KM 8,1 Pemantauan 7
H KM 2,65 Pemantauan 8
I KM 0,0 Pemantauan 9
 Point Source, Industri

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

Headwater 42,300 19,738 -7,00 -26 -43 112,00 27 25

Reach 1 RC1 36,900 19,738 12,212 -7,00 -24 -42 112,00 28 59

Reach 2 RC2 31,600 12,212 9,416 -7,00 -23 -50 112,00 31 33

Reach 3 RC3 21,700 9,416 6,512 -7,00 -22 -20 112,00 35 18

Reach 4 RC4 12,000 6,512 4,474 -7,00 -21 -6 112,00 39 44

Reach 5 RC5 5,600 4,474 4,375 -7,00 -19 -40 112,00 42 38

Reach 6 RC6 2,600 4,375 3,882 -7,00 -18 -29 112,00 43 12

Reach 7 RC7 0,000 3,882 1,798 -7,00 -18 -2 112,00 44 27

Interpretasi Hasil Pemodelan
Daya Tampung
 Gambar
apakah ini?
 Simulasi
• Simulasi dilakukan menggunakan model
kualitas air dengan software Qual2Kw ver 5.1
untuk mengestimasi kualitas air dengan
skenario yang dikehendaki.
• Parameter kualitas air yang disimulasi
disesuaikan dengan kemampuan Qual2Kw dan
ketersediaan data, misal: temperatur, pH, DO,
BOD, COD, TSS, N-organik, NO3-N, NH4-N, P-
organik, dan PO4-P.
 Contoh skenario simulasi
• Skenario 1: Kondisi sumber pencemar sesuai
kondisi eksisting pada debit minimum
• Skenario 2: Kondisi sumber pencemar memenuhi
baku mutu pada debit minimum
• Skenario 3: Kondisi tidak ada sumber
pencemaran pada debit minimum
• Skenario 4: Kondisi sumber pencemar di-trial-
error hingga kualitas air sepanjang sungai
memenuhi baku mutu sungai
Contoh Hasil Simulasi
 Hasil Simulasi skenario 1

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

Reach KM BOD COD

Reach 1 (Mlirip) 42,3-36,9 66,5% 41,6%
Reach 2 36,9-31,6 99,9% 99,8%
Reach 3 31,6-21,7 60,0% 62,7%
Reach 4 21,7-12,0 - -
Reach 5 12,0-5,6 82,7% 81,8%
Reach 6 5-6-2,6 63,9% 76,7%
Reach 7 (Jagir) 2,6-0,0 35,3% 39,3%
 TERIMA KASIH

8/10/18 63