Struktura szeregu czasowego natężenia dopływu ścieków do oczyszczalni
MONIżÿ CHUCHRO * Struktura szeregu czasowego natężenia dopÅ‚ywu Å›cieków do oczyszczalni SÅ‚owa kluczowe natężenie dopÅ‚ywu Å›cieki analiza statystyczna szeregi czasowe Streszczenie Analizowano statystycznie siedmioletnie szeregi czasowe przepÅ‚ywu o rozdzielczoÅ›ci dobowej. Badano ich strukturÄ™ czasowÄ…, skÅ‚adowe periodyczne oraz przyrodniczÄ… korelacje pogodowÄ…, a także korelacje z tygodniowÄ… aktywnoÅ›ciÄ… przemysÅ‚owÄ… i komunalnÄ…. Wysoka rozdzielczość i dÅ‚ugość szeregów pozwala poszukiwać stosunkowo sÅ‚abo odzwierciedlajÄ…cych siÄ™ zależnoÅ›ci. Zas- tosowano standardowe metody analizy spektralnej oraz autokorelacji. Do badania cyklicznoÅ›ci tygodniowej zastosowano me- tody analizy falkowej. Analizowane dane pochodziÅ‚y z trzech komunalnych oczyszczalni Å›cieków w: Krakowie oczyszczalnia Kujawy, Sandomierzu oraz w Warszawie oczyszczalnia Czajka. Podstawowym badanym parametrem byÅ‚ dopÅ‚yw Å›cieków do oczyszczalni. DziÄ™ki analizie falkowej udaÅ‚o siÄ™ wykonać tygodniowy model zmiennoÅ›ci dopÅ‚ywu Å›cieków do oczyszczalni Å›cieków. Dodatkowo analiza spektralna potwierdziÅ‚a prawidÅ‚owość utworzonej falki oraz wystÄ™pujÄ…ce w dopÅ‚ywie Å›cieków okresowoÅ›ci. Wprowadzenie Proces oczyszczania Å›cieków zależy od wielu czynników. Jednym z najważniejszych z nich jest natężenie dopÅ‚ywu Å›cieków do oczyszczalni. Zapisywane w regularnych odcinkach czasu wartoÅ›ci parametrów, w tym także natężenia dopÅ‚ywu, stanowiÄ… typowe szeregi czasowe. Analizowane szeregi czasowe pochodzÄ… z trzech komunalnych oczyszczalni Å›cieków, znajdujÄ…cych siÄ™ w dorzeczu WisÅ‚y. Pierwszym z obiektów jest oczyszczalnia Å›cieków w Krakowie. Jest to duża i nowoczesna oczyszczalnia dziaÅ‚ajÄ…ca od 1999 roku, przyjmujÄ…ca Å›cieki z dzielnicy Nowa Huta, o Å›rednim dobowym natężeniu dopÅ‚ywu wynoszÄ…cym 52,0 tys. m3. Sandomierska oczyszczalnia jest Å›redniej wielkoÅ›ci obiektem, jednoczeÅ›nie najmniejszÄ… oczyszczalniÄ… Å›cieków spoÅ›ród analizowanych. Åšrednie natężenie dopÅ‚ywu wynosi w niej 2,9 tys. m3/dobÄ™. NajwiÄ™kszÄ… badanÄ… oczyszczalniÄ…, a jednoczeÅ›nie najwiÄ™kszÄ… w Polsce, jest warszawska oczyszczalnia Å›cieków. Åšrednie dobowe natężenie * Akademia Górniczo-Hutnicza, Kraków e-mail: monika-chuchro@o2.pl 209 Chuchro M.: Struktura szeregu czasowego natężenia dopÅ‚ywu Å›cieków do oczyszczalni dopÅ‚ywu Å›cieków do tego obiektu wynosi 179,0 m3. Dla dokÅ‚adnoÅ›ci przeprowadzanych analiz wybrano szeregi cza- sowe o rozdzielczoÅ›ci dobowej. We wszystkich oczyszczalniach wybrany okres czasu (od 01.01.2000 do 31.12.2007) charakteryzowaÅ‚ siÄ™ najmniejszÄ… iloÅ›ciÄ… brakujÄ…cych danych. Analiza szeregów czasowych pozwala na okreÅ›lenie krótko- oraz dÅ‚ugookresowych zmian, wystÄ™pujÄ…cych w dopÅ‚y- wach do badanych oczyszczalni Å›cieków. Znajomość zmian zachodzÄ…cych w objÄ™toÅ›ci Å›cieków dopÅ‚ywajÄ…cych do oczyszczalni umożliwia przeprowadzenie procesu oczyszczania z wyższÄ… skutecznoÅ›ciÄ…. Modele zmiennoÅ›ci natężenia przepÅ‚ywu uÅ‚atwiajÄ… także zaprojektowanie nowych obiektów o wysokiej skutecznoÅ›ci oczyszczania Å›cieków. Głównym celem badaÅ„ byÅ‚o okreÅ›lenie struktury i korelacji szeregów czasowych pochodzÄ…cych z trzech oczyszczalni Å›cieków. 1. Podstawowe informacje o badanych szeregach czasowych Podstawowe informacje o natężeniu dopÅ‚ywu uzyskano obliczajÄ…c Å›redniÄ… arytmetycznÄ…, odchylenie standardowe, skoÅ›ność oraz współczynniki zmiennoÅ›ci i korelacji. Najbardziej zróżnicowane pomiary zaobserwowano w oczysz- czalni Å›cieków obsÅ‚ugujÄ…cej miasto Sandomierz, odchylenie standardowe dla tych danych wynosi 0,63, a współczynnik zmiennoÅ›ci 22%. W pozostaÅ‚ych oczyszczalniach wartość współczynnika zmiennoÅ›ci wynosi ok. 17%. Wartość tego współczynnika wskazuje na to, że znaczna ilość pomiarów miaÅ‚a wartość zbliżonÄ… do Å›redniej arytmetycznej (rys. 1), jednakże wystÄ™pujÄ… pomiary o wartoÅ›ci dużo wyższej i dużo niższej od wartoÅ›ci Å›redniej. Szeregi czasowe z Krakowa i Sandomierza wykazujÄ… silnÄ… prawostronnÄ… skoÅ›ność, której wartość przekracza 3. W danych pochodzÄ…cych z tych 2 oczyszczalni zaobserwowano niewiele pomiarów bardzo niskiego natężenia dopÅ‚ywu. Wysokie natężenie dopÅ‚ywu byÅ‚o natomiast zjawiskiem dość powszechnym (rys. 1). Dane z Warszawy wykazujÄ… skoÅ›ność, jednak jest ona dużo sÅ‚absza niż w pozostaÅ‚ych badanych oczyszczalniach Å›cieków i wynosi 0,9. Różnice w skoÅ›noÅ›ci poszczególnych szeregów czasowych widoczne sÄ… na rysunku 1. SkoÅ›ność danych zostaÅ‚a usuniÄ™ta przy użyciu przeksztaÅ‚cenia Boxa-Coxa. Rysunek 1. Histogramy szeregów czasowych: a) Kraków-Kujawy b) Sandomierz c) Warszawa-Czajka Figure 1. Time series histograms: a) Kraków-Kujawy b) Sandomierz c) Warszawa-Czajka Analiza danych obejmowaÅ‚a także wykonanie macierzy korelacyjnej. Najsilniejsza korelacja wystÄ™puje miÄ™dzy danymi z oczyszczalni sandomierskiej i krakowskiej i wynosi ona 0,4. NajsÅ‚absza korelacja wystÄ™puje natomiast miÄ™dzy danymi pochodzÄ…cymi z Krakowa oraz z Warszawy (0,29). Może być to spowodowane wiÄ™kszÄ… odlegÅ‚oÅ›ciÄ… miÄ™dzy oczyszczalniami, która powoduje wiÄ™ksze różnice w warunkach atmosferycznych. 2. Model natężenia dopÅ‚ywu Å›cieków Do utworzenia modelu natężenia dopÅ‚ywu wykorzystano analizÄ™ czasowo-skalowÄ… (falki). Falki pozwalajÄ… lepiej analizować procesy zawierajÄ…ce trendy deterministyczne i stochastyczne, a także sezonowość zmiennÄ…, zaÅ‚amania strukturalne oraz obserwacje nietypowe. Dane z oczyszczalni Å›cieków zostaÅ‚y zestandaryzowane osobno w obrÄ™bie 210 Krakowska Konferencja MÅ‚odych Uczonych 2008 każdego tygodnia, zaczynajÄ…c od czwartku. Zmiana pierwszego dnia rozpoczynajÄ…cego tydzieÅ„ obliczeniowy miaÅ‚a na celu wyrazniejsze zaznaczenie różnicy w natężeniu dopÅ‚ywu, jaka nastÄ™puje podczas weekendu. Kolejnym kro- kiem do utworzenia modelu falkowego byÅ‚o obliczenie Å›redniej wielkoÅ›ci natężenia dopÅ‚ywu dla poszczególnych dni tygodnia. NastÄ™pnie wykonano tygodniowy model dopÅ‚ywu Å›cieków do oczyszczalni (rys. 2). Nietypowość utwo- rzonej falki zwiÄ…zana jest z strukturÄ… badanych szeregów czasowych. Wielkość natężenia dopÅ‚ywu do oczyszczalni Å›cieków jako szereg czasowy wykazuje zmienność wewnÄ…trztygodniowÄ…, która jest zwiÄ…zana z dziaÅ‚alnoÅ›ciÄ… czÅ‚owieka. W badanych szeregach czasowych zaobserwowano silne nieregularnoÅ›ci wywoÅ‚ane zmianami pogodowymi. Rysunek 2. Tygodniowy model natężenia dopÅ‚ywu Å›cieków do oczyszczalni Å›cieków Figure 2. Weekly inżÿow model in the wastewater treatment plant W modelach natężenia dopÅ‚ywu Å›cieków wystÄ™pujÄ… zmiany zachodzÄ…ce w ciÄ…gu tygodnia, jak widać na rysunku 2 szczególnie widoczne w okresie czasu od piÄ…tku do niedzieli. Aby je uwidocznić, model tygodniowy zostaÅ‚ powtó- rzony. We wszystkich trzech falkach powtarza siÄ™ ten sam schemat: najniższe wartoÅ›ci natężenia dopÅ‚ywu wystÄ™pujÄ… w niedziele, maksimum natężenia dopÅ‚ywu w sobotÄ™, Å›rodowy wzrost wielkoÅ›ci natężenia dopÅ‚ywu. Zmiany te zwiÄ…zane sÄ… z dziaÅ‚alnoÅ›ciÄ… czÅ‚owieka. Maksimum sobotnie można powiÄ…zać ze wzmożonÄ… dziaÅ‚alnoÅ›ciÄ… w zakÅ‚adach przemysÅ‚owych, punktach usÅ‚ugowych oraz w gospodarstwach domowych. W niedzielÄ™ wiÄ™kszość zakÅ‚adów przemysÅ‚owych oraz punktów usÅ‚ugowych jest nieczynna. DziÄ™ki temu minimum wystÄ™pujÄ…ce w niedzielÄ™ można Å‚atwo wytÅ‚umaczyć. Od poniedziaÅ‚ku wysokość natężenia dopÅ‚ywu powoli wzrasta, by osiÄ…gnąć maksimum lokalne w Å›rodÄ™, po którym nastÄ™puje niewielkie obniżenie wielkoÅ›ci natężenia dopÅ‚ywu, trwajÄ…ce w czwartek i piÄ…tek. W przypadku Warszawy Å›rodowe i sobotnie maksima sÄ… sobie równe. WystÄ™powanie Å›rodowego maksimum jest ciekawym zjawiskiem, które być może możemy powiÄ…zać z wzmożonÄ… dziaÅ‚alnoÅ›ciÄ… w gospodarst- wach domowych oraz dziaÅ‚alnoÅ›ciÄ… maÅ‚ych przedsiÄ™biorstw odprowadzajÄ…cych Å›cieki do kanalizacji miejskiej. NastÄ™pnym krokiem analizy byÅ‚o porównanie modelu tygodniowego z danymi dotyczÄ…cymi natężenia dopÅ‚ywu. W tym celu wytypowano kilka okresów kilkutygodniowych. Okresy czasu wykorzystane w dalszej analizie wyznaczono przy użyciu korelacji rangowej. Do analizy wybrano przedziaÅ‚y czasowe, w których natężenie dopÅ‚ywu we wszystkich oczyszczalniach wykazujÄ… kilkutygodniowÄ… stabilność przy ujemnych lub dodatnich wartoÅ›ciach korelacji rangowej. PrzykÅ‚adem takiej stabilnoÅ›ci jest okres czasu od 160 do 166 tygodnia Å‚Ä…cznie, czyli od 18.01.2003 do 07.03.2003 roku. Okres ten wyróżniÅ‚ siÄ™ 7 tygodniowÄ… stabilizacjÄ… dodatniej korelacji rangowej dla badanych oczyszczalni Å›cieków. Wykre- sy dopÅ‚ywu Å›cieków znajdujÄ… siÄ™ poniżej. PorównujÄ…c przebieg wykresów, można zaobserwować podobieÅ„stwo miÄ™dzy wykresem pochodzÄ…cym z oczyszczalni z Krakowa i Sandomierza w drugiej poÅ‚owie wybranego okresu czasu. Niedzielne minimum modelu wystÄ™puje także w zestandaryzowanych szeregach natężenia dopÅ‚ywu Å›cieków do oczy- szczalni. Podobnie maksimum sobotnie można także zaobserwować w danych zestandaryzowanych (rys. 3, rys. 4, rys. 5). Dodatkowo w celu wyjaÅ›nienia wysokich i bardzo niskich wartoÅ›ci dopÅ‚ywu dane z wykresów porównano z danymi pogodowymi. Dane pogodowe dotyczÄ…ce nasÅ‚onecznienia i opadów pochodzÄ… z obszaru Krakowa (tab. 1). Przez prawie 211 Chuchro M.: Struktura szeregu czasowego natężenia dopÅ‚ywu Å›cieków do oczyszczalni caÅ‚y analizowany okres panowaÅ‚y dobre warunki atmosferyczne. W ciÄ…gu 49 dni tylko siedmiokrotnie zanotowano opady, zaobserwowano także jednodniowÄ… odwilż. Reszta okresu czasu byÅ‚a sÅ‚oneczna. Opady wystÄ™pujÄ…ce od 23.01.2003 do 28.01.2003 widoczne sÄ… na wykresie w postaci dość dużej wartoÅ›ci natężenia dopÅ‚ywu Å›cieków do oczyszczalni. Rysunek 3. Porównanie modelu natężenia dopÅ‚ywu z danymi z oczyszczalni Å›cieków z Krakowa z tygodni od 160 do166 Figure 3. Comparison of inżÿow model with the wastewater treatment plant data from Kraków, weeks from 160 to 166 Rysunek 4. Porównanie modelu natężenia dopÅ‚ywu z danymi z oczyszczalni Å›cieków z Sandomierza z tygodni od 160 do166 Figure 4. Comparison of inżÿow model with the wastewater treatment plant data from Sandomierz, weeks from 160 to 166 Tabela 1. Dane pogodowe dla miasta Krakowa za okres 2003-01-18 do 2003-03-07 (160 do 166 tydzieÅ„) Table 1. Weather data for Kraków, from 2003-01-18 to 2003-03-07 (weeks from 160 to 166) Data Pogoda Data Pogoda Data Pogoda Data Pogoda Data Pogoda 2003-01-18 S 2003-01-28 S/OP 2003-02-07 S 2003-02-17 S 2003-02-27 S 2003-01-19 S 2003-01-29 S 2003-02-08 S 2003-02-18 S 2003-02-28 ODW 2003-01-20 S 2003-01-30 S 2003-02-09 S 2003-02-19 S 2003-03-01 S 2003-01-21 S 2003-01-31 S 2003-02-10 S 2003-02-20 S 2003-03-02 S/OP 2003-01-22 S 2003-02-01 S 2003-02-11 S 2003-02-21 S 2003-03-03 OP 2003-01-23 OP 2003-02-02 S 2003-02-12 S 2003-02-22 S 2003-03-04 S 2003-01-24 OP 2003-02-03 S 2003-02-13 S 2003-02-23 S 2003-03-05 S 2003-01-25 OP 2003-02-04 S 2003-02-14 S 2003-02-24 S 2003-03-06 S 2003-01-26 S 2003-02-05 S 2003-02-15 S 2003-02-25 S 2003-03-07 S 2003-01-27 S/OP 2003-02-06 S 2003-02-16 S 2003-02-26 S S dzieÅ„ soneczny, OP opady, ODW odwilż 212 Krakowska Konferencja MÅ‚odych Uczonych 2008 Rysunek 5. Porównanie modelu natężenia dopÅ‚ywu z danymi z oczyszczalni Å›cieków z Warszawy z tygodni od 160 do166 Figure 5. Comparison of żÿow model with the wastewater treatment plant data from Kraków, weeks from 160 to 166 3. Analiza widmowa szeregów czasowych DopÅ‚ywy Å›cieków do oczyszczalni poddano dodatkowo Szybkiej Transformacie Fouriera. WartoÅ›ci periodogramu zostaÅ‚y wygÅ‚adzone przy pomocy ważonej Å›redniej ruchomej okna Hamminga. WartoÅ›ci wag Hamminga widoczne sÄ… na rysunku 6. Na wszystkich trzech wykresach analizy widmowej zaobserwowano zwiÄ™kszonÄ… gÄ™stość widmowÄ… w tych samych punktach. Na wykresie analizy widmowej (rys. 6) dotyczÄ…cym szeregu czasowego z oczyszczalni Å›cieków w Sandomierzu można zaobserwować piki przy czÄ™stotliwoÅ›ciach 0,0180; 0,0324; 0,074; 0,1482; 0,284; 0,48, a także odpowiadajÄ…cÄ… im gÄ™stość widmowÄ…. CzÄ™stotliwość 0,284 odpowiada okresowi 3,5 dnia. Analiza widmowa potwierdziÅ‚a wiÄ™c prawidÅ‚owość utworzonego modelu dopÅ‚ywu Å›cieków. PozostaÅ‚e czÄ™stotliwoÅ›ci odpowiadajÄ… kolejno okresom: 2 miesiÄ…ce, 1 miesiÄ…c, 13,5 dnia, 6,7dnia oraz 2,1 dnia. Na wykresie analizy widmowej widoczna jest czÄ™stotliwość odpowiadajÄ…ca 3,5 dnia oraz jej harmoniczne, a także czÄ™stotliwość odpowiadajÄ…ca w przybliżeniu zmiennoÅ›ci ty- godniowej oraz jej harmoniczne. Rysunek 6. Analiza widmowa natężenia dopÅ‚ywu do oczyszczalni Å›cieków w Sandomierzu Figure 6. Inżÿow spectrum analysis for Sandomierz wastewater treatment plant Podsumowanie RozkÅ‚ad analizowanych szeregów czasowych jest niesymetryczny i wykazuje prawostronnÄ… skoÅ›ność; 213 Chuchro M.: Struktura szeregu czasowego natężenia dopÅ‚ywu Å›cieków do oczyszczalni MiÄ™dzy danymi pochodzÄ…cymi z różnych oczyszczalni Å›cieków wystÄ™puje dodatnia korelacja wynoszÄ…ca od 0,29 do 0,4; Utworzona falka prawidÅ‚owo przedstawia zmienność dopÅ‚ywu Å›cieków do oczyszczalni w ciÄ…gu tygod- nia. Potwierdza to porównanie modelu dopÅ‚ywu tygodniowego z danymi zestandaryzowanymi. Okresy niezgodnoÅ›ci modelu z danymi mogÄ… być zwiÄ…zane z warunkami pogodowymi, co można zaobserwować porównujÄ…c okresy niezgodnoÅ›ci z danymi dotyczÄ…cymi pogody; Modele przepÅ‚ywu utworzone dla oczyszczalni w Sandomierzu i w Krakowie majÄ… podobny przebieg i war- toÅ›ci dla poszczególnych dni; W analizowanych szeregach czasowych zaobserwowano okresowość 3,5-; 7-; 14-; 30-; 60-dobowÄ…; Możliwe jest utworzenie dokÅ‚adniejszego modelu dotyczÄ…cego dopÅ‚ywu Å›cieków do oczyszczalni Å›cieków, przy uwzglÄ™dnieniu zmiennoÅ›ci okresowej oraz w arunków pogodowych; Stworzenie modelu pozwalajÄ…cego z wysokÄ… skutecznoÅ›ciÄ… przewidywać natężenie dopÅ‚ywu pozwoliÅ‚oby na usprawnienie procesu oczyszczania Å›cieków; Model natężenia dopÅ‚ywu Å›cieków może być przydatnym elementem podczas projektowania nowych obiek- tów, a także modernizacji już istniejÄ…cych oczyszczalni Å›cieków. Literatura [1] Hamilton J.D.: Time Series Analysis, Princeton University Press, New Jersey 1994. [2] Last M., Kandel A., Bunke H.: Data Mining in Time Series Databases, [in:] Series in Machine Perception and Artiżÿcial Intelligence, World Scientiżÿc 2004. [3] www.data-mining.pl/textbook/sżÿimser.html. MONIżÿ CHUCHRO Time Series Structure of Wastewater Inżÿow into Sewage Treatment Plant Keywords inżÿow wastewater statistical analysis time series Abstract żÿe overall goals of investigation were analysis of inżÿow time series from wastewater treatment plant (WWTP). We have daily data of seven years time period from three WWTP: Kraków-Kujawy, Sandomierz and Warszawa-Czajka. We investigated time structure, periodic components and weather correlation of time series. We also checked correlation of inżÿow time series with weekly industry and municipal data activity. Time series high resolution and length helps in searching weakly visible relationships. For investigation we used standard spectrum analysis and autocorrelation methods. We used wavelet methods for checking weekly periodic. żÿe wavelet methods helped to made weekly model of inżÿow changeability for data from WWTP. Spectrum analysis corroborated correctness of inżÿow model, and weekly periodic found in inżÿow data. In the model we observed changes during a week. On Saturday inżÿow had the biggest value, from the all week. żÿe second high value we observed on Wednesday, the lowest values of inżÿow were on Sunday. Comparison of inżÿow data and model showed similarities, between graphs. 214