Streeter-Phelps er en modelleringsmetode for iltindholdet i åer og floder. Den tager højde for ændringen af iltindholdet, som følge af forurening med organisk materiale. Forureningen medfører en biologisk nedbrydning, hvor ilt forbruges. Streeter-Phelps standardformel tager kun højde for den biologiske nedbrydning af organisk materiale, mens den modificerede formel giver mulighed for at tilføje flere processer, der ligeledes påvirker iltkoncentrationen i åer og floder.
Streeter-Phelps standardformel:
Løsningen af differentialligningen:
Hvor D er givet ved mg O2/L.
Historie
H. W. Streeter og Earle B. Phelps udgav i 1925 en rapport over undersøgelser af Ohio River. Her publicerede de for første gang deres matematiske model for iltudviklingen i forurenede floder. Ud fra målinger udført i årene 1914-1915 af forureninger samt den medfølgende biologiske nedbrydning, kunne Streeter og Phelps udlede en matematisk sammenhæng, der stadig anvendes verden over.
Phelps var ansat ved U.S. Public Health Service i Cincinnati, Ohio, hvor han var leder af afdelingen for kemi og hygiejne. I 1912 forlangte den Amerikanske Kongres en undersøgelse af den skadelige effekt på mennesker ved udledning af kloakvand. Phelps udviklede i samarbejde med Streeter en matematisk 1. ordens nedbrydningsreaktion for iltkoncentrationen. Ud fra målinger af carbonholdige forbindelser i Ohio River, der kan nedbrydes biologisk (biokemisk iltforbrug – BOD), og iltkoncentrationen i floden (dissolved oxygen – DO) kunne en nedbrydningskonstant (kBOD) findes. Ved at benytte Streeter-Phelps formel kan det udregnes hvor langt nedstrøms fra forureningskilden, at den mindste iltkoncentration forekommer.EEJournalV13P1.indd
Efter publiceringen af standardformlen for iltnedbrydningen, hvor kun forureningens iltoptag var medtaget, blev der tilføjet flere processer, der påvirker iltkoncentrationen. Dette medførte en mere avanceret og modificeret model, hvor forskellige processer kunne medtages i forhold til relevans for den enkelte flod. De forskellige processer, der kan medtages, er vist under afsnit 2.2 Modificeret model.
Den modificerede formel:
Løsningen af differentialeligningen:
Hvor D er givet ved mg O2/L.
Udledning I udledningen af Streeter-Phelps formlen vil der ikke blive gjort forskel på det kvælstofholdige biokemiske iltforbrug (NBOD) og det kulstofholdige biokemisk iltforbrug (CBOD). I andre kilder såsom i Schnoor (1996)Schnoor, J. L., Environmental Modeling, Fate and Transport of Pollutants in Air, Water and Soil, vil denne forskel blive gjort. Først udledes den oprindelige Streeter-Phelps formel. Udledningen er baseret på en Plug Flow Rector (PFR) som beskriver kemiske reaktioner mens de transporteres i en rørledning, i dette tilfælde repræsentere rørledningen en flod.
Oxideringen af BOD ved en konstant strømningshastighed og ligevægtstilstand kan beskrives som en 1. ordens reaktion:
hvor BOD = koncentrationen af BOD
kBOD = 1. ordens nedbrydningsrate
Løsningen til differentialligningen fås ved:
Herefter integreres begge sider:
Hvorefter løsningen er:
Det er dernæst muligt at opstille en simpel massebalance, som beskriver udviklingen af iltindholdet som funktion af tiden:
hvor DO = Opløst ilt koncentration
DOsat = Mættet opløst ilt koncentration
kO2 = 1. ordens geniltningsrate
De to led på højre side af ligningen beskriver således iltforbruget fra BOD og graden af udveksling af ilt med atmosfæren. Denne ligning kan dog med fordel beskrives i form af underskuddet af opløst oxygen i vandet, dvs. D = DOsat – DO, hvor D er opløst oxygen deficit. Det medfører at:
Det ses, at fortegnet for de to led på højre side af ligningen ændres grundet det omvendte fortegn for D i forhold til DO. D beskriver underskuddet af DO, hvorfor et input af oxygen vil have en negativ effekt på D og vice versa. Den ordinære differentialligning (2) kan efterfølgende løses. Løsningen (1) indsættes i (2) og alle led indeholdende D samles på venstre side af lighedstegnet.
Ligning (3) kan skrives ved hjælp af Integration Factor Method på formen:
Ved at sammenholde ligning (3) og (4) fås følgende definition hvor:
Metoden for løsningen af ligning (4) er baseret på en integrations faktor, M:
Integrationsfaktoren ganges på begge sider af ligning (4):
På grund af den valgte integrations faktor kan M×p(t) erstattes af dM/dt hvilket simplificerer ligningen til:
Ved at benytte produktreglen kan venstre side omskrives til:
Ved at integrere begge sider fås:
Løsningen for differentialligningen (4) vil være på formen:
Konstanten C er lig startkoncentrationen i deficit, D(0). Løsningen til differentialligningen (3) bestemmes ud fra den generelle løsning (5):
Det bestemte integrale er:
Det kan indsættes i (6) og løsningen til differentialligningen er:
Ligning (7) beskriver udviklingen i opløst oxygen deficit, D, som funktion af tiden, t, udtrykt ved det oprindelige opløste oxygen deficit, D(0) og den oprindelige BOD koncentration, BOD(0), samt reaktionskonstanterne kBOD og kO2. Hvilket er den oprindelige Streeter-Phelps ligning. Koncentrationen af opløst oxygen, DO, i vandet kan slutteligt findes ved relationen mellem D og DO således:
Kritiske distance, tid og deficit Udviklingen af D vil typisk forløbe således, at D stiger til et maksimalt kritisk deficit, for derefter at falde igen, hvorved den opløste oxygen koncentration stiger og går mod det oprindelige niveau igen. Det kritiske opløste oxygen deficit, Dcrit, og den dertilhørende kritiske tid, tcrit, kan forholdsvist simpelt beregnes. Følgende gør sig gældende for t = tcrit:
Ved at indsætte dette i (3), og løse den for Dcrit fås:
Et udtryk for tcrit fås ved at substituere (8) ind i (7) og isolere tcrit.
Den kritiske distance findes ved:
Modificeret model Den modificeret model tager højde for flere processer, herunder ses eksempler på iltproducerende og iltforbrugende processer, der kan inkluderes i modelleringen. Tabel 1 - Processor relateret til oxygen Dannelse Forbrug Fotosyntese Respiration Geniltning Forbrug via sediment Baggrunds BOD BOD
Processerne kan løses enkeltvist ved først at opstille en massebalance og derefter benytte superpositionsprincippet til at finde det samlede udtryk. Resultatet bliver en langt mere omfattende Streeter-Phelps formel.
Baggrunds BOD Baggrunds BOD indikerer et non-point input til den respektive flod, å, bæk osv., hvilket kan skyldes afstrømning fra landbrug og motorvej samt regnvands udledning. Baggrunds BOD forårsager et generelt underskud af opløst oxygen under ligevægtstilstand.
Ændringen i BOD over tid inklusiv baggrunds BOD er:
For en konstant BOD gælder:
Massebalancen for opløst oxygen deficit:
Under ligevægtstilstand (ændringen i deficit over tid er lig nul) og med konstant baggrunds BOD, er opløst oxygen deficit en konstant:
Det ovenstående er et udtryk for baggrunds koncentration af opløst oxygen deficit.
Sedimentering af BOD Yderligere kan der tages højde for mulig sedimentering af BOD. Total suspenderet faststof (TSS) indeholder organisk materiale, som muligvis er iltforbrugende. Streeter-Phelps formel (den oprindelig) modificeres således, at den nu tager højde for mulig sedimentering af BOD. Det er nødvendigt med en kombineret rate bestående af en sedimenteringsrate for BOD, ks samt kBOD, som medtages i en fælles rate, kr:
kr = ks + kBOD
Ændringen i BOD over tid:
Massebalancen for opløst oxygen deficit, D:
Løsningen til differentialligningen findes ved hjælp af Integration Factor Method og er:
Sedimentets iltforbrug Sedimentet forbruger ilt til nedbrydning af bundfældet organisk materiale og tabet antages for værende konstant:
Massebalance for opløst oxygen deficit:
Løsningen til differentialligningen findes ved hjælp af Integration Factor Method og er:
Fotosyntese og respiration Respiration, R, vil være dominerende om natten og bidrage til opløst oxygen deficit hvorimod om dagen dominerer fotosyntesen, P, og opløst oxygen er dannet.
Massebalance for opløst oxygen deficit:
Løsningen til differentialligningen er:
Det er således muligt at opstille den samlede massebalance for DO deficit for en plug flow flod:
Løsningen til differentialligningen er summen af de enkelte løsninger som konsekvens af superpositionsprincippet:
Parameter Antallet af parametre i Streeter-Phelps formel afhænger af, om den er udvidet til at indeholde mere end forureningens iltforbrug og udvekslingen af ilt med luften, dvs. om hvorvidt der er tale om den modificeret model eller ej.
Tabel 2 - Oversigt over parametre til Streeter-Phelps Symbol Symbolforklaring Forhold Typiske værdier Endhed kO2 GeniltningsratenDybe åer med lavt flow 0,1 - 0,4 Typiske forhold 0,4 – 1,5 Dybe floder, hurtigt flow 1,5 – 4,0 Lavvandede floder, hurtigt flow 4,0 – 10 DOpløst oxygen deficitDOKoncentarion af opløst oxygenDOsatMættede koncentration af opløst iltMellem 10 °C – 20 °C10,92 – 8,84kBOD BOD fjernelsesraten Temperaturafhængig 0,05 – 1,7BODBiokemisk iltforbrugSSedimentets iltforbrugIngen til lav forurening0,1 – 10Moderat til meget forurening5 – 10HGennemsnitlige dybdekr0,5 – 5ks Sedimenteringsraten for BODRRespiration1 – 10PFotosynteseMaximal daglig produktion2 – 20(P-R)Netto primær iltproduktionDagligt gennemsnit0,5 – 10
k02 k02 afhænger i høj grad af vandmassens fysiske udformning, dvs. om der er tale om en å, flod, sø eller åbent hav. Ydermere er temperaturen en vigtig faktor, da den afgør mætningskoncentrationen af ilt i vandet. k02 beskriver raten af geniltning fra atmosfæriske luft til vandmassen, og sammenhængen er beskrevet ved Fick’s første lov om diffusionTrapp, S. & Matthies, M., Chemodynanics and Environmental modeling – An introduction, .
Generelt kan raten af geniltningen findes ud fra følgende formel:
Hvor a og b afhænger af vandmassens fysiske udformning. Flowhastigheden (v) opgives i m/s og middeldybden (H) i meter, derved er enheden for kO2 per døgn (d-1). En typisk værdi ved 20 °C ville ligge på mellem 0,4-1,5 d-1, men den kan ligge mellem 0,1 d-1 og helt op til 10 d-1. Den situationspecifikke værdi kan findes ved hjælp af følgende formler:
Tabel 3 - Geniltningsrate ved forskellige strømningsforhold Forhold Formel RegressionLangsomt strømmende i åer Langbein-Durum Store floder med H>3,5 m og v>0,5 m/s Southworth/Churchill Mix af forskellige floder IHP/ Jolánkai Åbent hav Liss og Slater Søer, ingen vind Mackay og Yeun
Temperaturkorrigering kan foretages ved brug af følgende formel:
Deficit,D Tabel 4 - Iltmætning ved forskellige temperatureVDG, 1989, Water quality determination standard protocol, Germany [°C] [mg O2/L] [°C] [mg O2/L] [°C] [mg O2/L] [°C] [mg O2/L] 0 14,16 1 13,77 11 10,67 21 8,68 31 7,42 2 13,40 12 10,43 22 8,53 32 7,32 3 13,05 13 10,20 23 838 33 7,22 4 12,70 14 9,98 24 8,25 34 7,13 5 12,37 15 9,76 25 8,11 35 7,04 6 12,06 16 9,56 26 7,99 36 6,94 7 11,76 17 9,37 27 7,86 37 6,86 8 11,47 18 9,18 28 7,75 38 6,76 9 11,19 19 9,01 29 7,64 39 6,68 10 10,92 20 8,84 30 7,53 40 6,59
D er oxygen deficit, dvs. underskuddet i oxygen koncentrationen i forhold til den fuldstændig mættet koncentration. Da den mættede koncentration afhænger af temperaturen, variere D ligeledes med temperaturen. Deficit kan beregnes ved:
Hvor DO er opløst oxygen, og DOsat er den mættede koncentration af ilt i vandet. DOsat kan aflæses i tabellen til højre.
DOsat for specifikke temperature kan udregnes ved brug af følgende formelJolánkai, G., Basic river water quality models, computing aided learning (CAL) programme on water quality modeling (WQMCAL version 1.1):
kBOD
kBOD er fjernelsesraten for BOD, og den kan bestemmes eksperimentel ud fra BOD målingerne. Den har en typisk værdi på mellem 0,1-1,7 d-1.
Den er temperatur afhængig og kan korrigeres ud fra følgende formel:
BOD Denne parameter måles for udløbet til åen, desuden måles BOD i åen før udløbet, og derefter beregnes en vægtet værdi for BOD for åen efter udløbet. Denne BOD-værdi efter et udløb kan udregnes ved brug af følgende formel: