Preferencje help
Widoczny [Schowaj] Abstrakt
Liczba wyników

Znaleziono wyników: 6

Liczba wyników na stronie
first rewind previous Strona / 1 next fast forward last
Wyniki wyszukiwania
Wyszukiwano:
w słowach kluczowych:  flow with mass and momentum exchange
help Sortuj według:

help Ogranicz wyniki do:
first rewind previous Strona / 1 next fast forward last
PL
Zjawisko kolmatacji występuje w przyrodzie wszędzie tam, gdzie dochodzi w ośrodkach porowatych do przepływów cieczy niosących zawieszone cząstki stałe. Nawet „najczystsza” woda dopływająca do studni po pewnym czasie spowoduje jej zakolmatowanie, a tym samym spadek jej wydajności, co jest zjawiskiem negatywnym. Badania prowadzone w naszym ośrodku od lat 60-tych ubiegłego wieku [1,4-14] doprowadziły do opracowania matematycznego opisu zjawiska kolmatacji [4-8,13,14] oraz przeprowadzenia szeregu eksperymentów ją weryfikujących [9-14]. Uzyskane wyniki wykorzystano również podczas prób uszczelnienia górotworu wokół wyrobiska górniczego [12]. W niniejszym artykule podjęto próbę określenia obszaru К w przypadku przebiegu zjawiska kolmatacji zachodzącego zgodnie z kinetyką pierwszą, oraz podania zależności opisujących rozkład ciśnienia h(x,t) dla przepływu bez kolmatacji i z kolmatacją bez linearyzacji wyrażenia ε(x,t)-3 w otoczeniu εo gdzie ε określa porowatość ośrodka w funkcji położenia i czasu. Określenie obszaru К pozwala nam na jednoznaczne wyprowadzenie dokładnego rozkładu ciśnienia h(x,t) podczas przepływu z kolmatacją przez ośrodek porowaty bez linearyzacji, a następnie porównania rozwiązań układu równań kolmatacji metodą linearyzacji i metodą dokładną przy wykorzystaniu bezwymiarowej postaci funkcji ξ. W trakcie prowadzenia badań eksperymentalnych próba dopasowania rzeczywistego zjawiska do opracowanego modelu matematycznego obarczona była dużą niepewnością wynikającą prawdopodobnie z zastosowania linearyzacji członu ε(x,t)-3. W artykule autorzy wyjaśniają co generuje odstępstwo metody przybliżonej od dokładnego rozwiązania oraz zwracają uwagę, że dokładne rozwiązanie bardziej oddaje sens fizyczny matematycznego modelu opisu zjawiska oraz zdefiniowanych współczynników kolmatacji a w szczególności parametru αo
EN
The phenomenon of colmatage occurs in nature wherever there is a flow of fluid carrying suspended solid particles through porous media. Even the "cleanest" water flowing into the well after some period of time will become clogged and therefore its efficiency will decrease, which is a negative phenomenon. Research conducted in our center, since the 1960s [1,4-14], has led to: a theoretical description of the phenomenon of colmatage [4-8,13,14] and a number of experiments verifying it [9-14]. The obtained results were used during tests to seal the rock mass around a mining excavation [12]. This article attempts to determine the area К in the case of the colmatage phenomenon occurring in accordance with the first kinetics, and to identify/formulate relationships describing the pressure distribution h(x,t) for the flow without colmatage and with colmatage without linearization of the expression ε(x,t)-3 in the surroundings εo , where ε determines the porosity of the medium as a function of position and time. Determining the area К allows us to clearly derive the exact pressure distribution h(x,t) during flow with colmatage through a porous medium without linearization, and then compare the solutions of the system of colmatage equations using the linearization method and the exact method using the dimensionless form of the ξ function. During the experimental research, the attempt to match the actual phenomenon with the developed mathematical model was burdened with high uncertainty, probably resulting from the use of linearization of the ε(x,t)-3 term. In the article, the authors explain what generates the deviation of the approximate method from the exact solution and point out that the exact solution better reflects the physical meaning of the mathematical model for describing the phenomenon and the defined colmatation coefficients, the parameter εo , in particular.
2
Content available remote Flows with substance exchange and change of its volume in radial geometry
EN
The subject of this paper is a certain model of the process of eolmatage in which the swelling of colmatant settled in a porous medium proceeds. An axi-symmctric flow is considered. The paper is a continuation of the article (Trzaska, Sobowska 2000) in which a similar problem for the one-dimensional case was examined. It is assumed that the swelling occurs when colmatant contacts a certain substance. In the accepted model it is assumed that suspension carrying solid particles of the colmatant is initially forced into the porous medium. The particles settle in the pores. This part of the process is conventionally called its first stage. The stage is described by a system of balance-transport (1) and kinetics (2) equations. The obtaining of a certain distribution of colmatant P0(r) settled in the porous medium (16) is the result of stage one of the process. During stage two, substance which can cause the swelling of the colmatant settled earlier is forced into the medium. Part of this substance stays in the medium pores reacting on the colmatant and leading to the increase of its volume. This stage of the process is described by equations (17), (18), (24). The distribution of the colmatant in the medium pores P(r,t) is obtained in the form (25). In this paper computations arc made which help to determine the distribution of pressure in the medium during stage two of the process. The equation of motion (26) was used. The flow at known discharge q(t) was examined. The distribution of pressures h(r,t) is expressed then by formula (29).
PL
Niniejsza publikacja jest kontynuacją pracy Trzaski i Sobowskicj (2000), w której podobny problem rozpatrywany był dla przypadku przepływu jednowymiarowego wzdłuż tworzących równoległych do osi przepływu. Tutaj rozważany jest model osiowo-symetryczny. Obydwa wymienione artykuły należą do serii prac poświęconych kolmatacji i będących efektem badań tego zjawiska prowadzonych przez ten zespół autorów. Przypomnijmy, że przez kolmatację rozumiemy zjawisko zachodzące w trakcie przepływu zawiesiny przez ośrodek porowaty i polegające na wymianie cząstek stałych z ośrodka ciekłego do stałego. W efekcie dochodzi do zmian własności fizycznych ośrodka porowatego, zmienia się w szczególności jego porowatość i w konsekwencji współczynnik przepuszczalności. To z kolei powoduje zmianę ciśnienia oraz przy pewnych uwarunkowaniach wpływa na wielkość wydatku przepływu. W dotychczasowych opracowaińach rozpatrywano procesy, w których własności fizyczne kol-matanta po jego osadzeniu się w porach ośrodka nic ulegają zmianie. W niniejszym artykule, podobnie jak w pracy Trzaski i Sobowskicj (2000) rozważany jest model proeesu, w którym dochodzi do zmian objętości cząstek kolmatanta osadzonego w ośrodku porowatym. Przyjęto założenie, że zmiany te (pęcznienie, namnażanie) następują w wyniku kontaktu kolmatanta z pewną substancją. W pracy badany jest model procesu składającego się z dwóch etapów. Pierwszy etap to zatłaczanie do ośrodka zawiesiny niosącej cząstki stałe. W trakcie przepływu zawiesiny dochodzi do osadzania się tych cząstek w przestrzeni porowej, zachodzi zjawisko kolmatacji. Ten etap opisano układem równań bilansu-transportu (1) i kinetyki (2). W wyniku rozwiązania układu (1), (2) z warunkami początkowo--brzegowymi (3), (5) uzyskano funkcję rozkładu kolmatanta w ośrodku porowatym w dowolnej chwili t trwania procesu, a w szczególności w chwili tt zakończenia tego etapu (16). Uzyskane rozwiązanie stanowi warunek początkowy dla drugiego etapu procesu. Teraz do ośrodka zatłaczana jest substancja powodująca pęcznienie lub namnażanie kolmatanta osadzonego wcześniej. Do opisu tego etapu użyto układu równań (17), (18), (24), który rozpatrywano z warunkami początkowo-brzegowymi (19), (20), (23). Uzyskano funkcje koncentracji omawianej substancji w cieczy płynącej przez ośrodek oraz zatrzymanej w przestrzeni porowej i reagującej z kolmatantem. Następnie wyznaczono funkcję rozkładu pęczniejącego kolmatanta w ośrodku porowatym (25). Kolejnym etapem pracy było wyznaczenie rozkładu ciśnienia panującego w ośrodku porowatym w trakcie drugiego etapu procesu. Wykorzystano równanie ruchu postaci (26) i uzyskano funkcję rozkładu ciśnienia dla przepływu przy zadanej, na przykład stałej wartości wydatku przepływu (29).
3
Content available remote Process of colmatage with transient boundary condition
EN
The subject of this publication is a certain model of the process of colmatage in a porous medium with a closed circulation of suspension. The process is investigated in which a liquid flowing out of the medium flows into the container filled with suspension. It is mixed there with the rest of the liquid and forced back into the medium. The investigations are carried out on the basis of the system of balance-transport (2) and kinetics (3} equations and on that of the balance of forced suspension (9). Initial-boundary conditions are accepted in the form (4), (5). Function of time n(t) of forced suspension (13) have been obtained followed by the distribution of concentration of suspension N(x,t) (14) flowing through the medium, and the distribution of the medium porosity E[epsilon](x,t) (16). Basing on the equation of motion (17) the distribution of pressure in the porous medium (19) has been determined.
PL
Tematem niniejszej publikacji jest pewien model przebiegu procesu kolmatacji w ośrodku porowatym przy zamkniętym obiegu zawiesiny. W trakcie takiego procesu koncentracja zatłaczanej do ośrodka zawiesiny, która w chwili t = 0 posiada wartość n0, zmienia się wskutek osadzania w ośrodku transportowanych przez ciecz cząstek. Wartość koncentracji na wlocie nie ulega zmianie przez okres wyznaczony dojściem czoła fali z punktu x = 0 do punktu x = L. W tym momencie, który w pracy oznaczamy jako t = t1, wypływającą z ośrodka zawiesinę nawracamy do zbiornika, z którego jak poprzednio po dokładnym, permanentnym wymieszaniu jest zatłaczana do ośrodka porowatego. Poczynając od chwili t = t1 koncentracja zawiesiny na wejściu do ośrodka staje się funkcją czasu n(t). Na jej wartość wpływa koncentracja zawiesiny w zbiorniku w chwili t = 0, oraz charakter przebiegu zjawiska, w wyniku którego koncentracja na wyjściu z ośrodka może przyjmować różne wartości 0 < N (L,t) < n0. Gdy N (L,t) = n0 nie zachodzi proces kolmatacji. Wtedy przepływająca przez ośrodek zawiesina nie poddawana jest wymianie masy z otoczeniem: z ośrodka ciekłego do porowatego. W przypadku przepływu z wymianą masy zachodzi warunek 0 < N (L,t) < n0. Z punktu widzenia matematycznego opisu obiektem naszego zainteresowania są przepisy funkcyjne takich wielkości, jak rozkład koncentracji unoszonych i zatrzymanych w ośrodku porowatym cząstek kolmatanta, opis rozkładu porowatości ośrodka, a co za tym idzie i rozkład ciśnień do jakiego dochodzi w wyniku przebiegu omawianego procesu na drodze x i w czasie t jego trwania. Wymieniony opis teoretyczny podajemy w oparciu o układ równań bilansu-transportu i kinetyki procesu kolmatacji, który to układ ze względu na przyjęty model przebiegu zjawiska ma postać daną wzorami [2) i (3) z warunkami początkowo-brzegowymi (4), (5). W uzyskanym rozwiązaniu powyżej omawianych równań dostajemy funkcję określającą rozkład przepływających cząstek kolmatanta N(x,t) w postaci (8) z niewiadomą funkcją n(t), która występuje w nieustalonym warunku brzegowym (5) opisywanego zjawiska. Przepis na wymienioną funkcję otrzymujemy rozwiązując liniowe równanie różniczkowe, do którego dochodzimy dokonując bilansu cząstek stałych znajdujących się w zbiorniku w różnych czasach t. Rozwiązanie uzyskanego równania (9) otrzymujemy stosując przekształcenie Laplace'a. Sposób rozwiązania przedstawiono w apendyksie. Uzyskaną w wyniku funkcję n((t) (12) wykorzystujemy podstawiając ją do wzoru [7) i otrzymując w ten sposób przepis informujący o rozkładzie przepływających cząstek N(x,t) (14). Drugą szukaną funkcję e[epsilon](x,t) dostajemy podstawiając wzór (14) do równania (3) i całkując je z warunkiem początkowym (4). Postać tej funkcji dana jest wzorem (16). Kolejna funkcja, której poszukujemy i podajemy, informuje o rozkładzie ciśnienia i jego zmienności w trakcie przebiegu zjawiska kolmatacji, przy czym zakładamy, że przepływ zachodzi ze stałą prędkością filtracji q. Rozkład ciśnienia uzyskujemy całkując równanie ruchu (17) z warunkiem (18) po podstawieniu do niego wyliczonej poprzednio funkcji e(x,t) danej wzorem (16). Otrzymana funkcja h[x,t) (19) jest ostatnim przepisem opisującym całokształt zagadnień związanych z przebiegiem zjawiska kolmatacji zachodzącego w omawianym ośrodku.
4
Content available remote Process of colmatage in porous medium with closed circulation of suspension
EN
In this paper a theoretical description of the phenomenon of colatage observed in a porous medium with a finite lenght "L" and a closed circulation has been presented. Such circulation allows one to used the same suspension many times because after flowing through the medium, it can be forced again into it. Theoretical considerations have been presented on the basis of a systems of balance-transport equation (1) and those of the kinetics of the colmatage process (2) with initial-boundary conditions (3), (4), (5). Functions obtained determine the distribution of the concentration of flowing suspension "N" (x,t) (17) and the porosity epsilon (x, t) (18) in space x and time "t" of the proceeding phenomenon. Based on equation of motion (20) the distribution of pressure in porous medium were determined (23).
PL
W niniejszej publikacji został przedstawiony teoretyczny opis zjawiska kolmatacji obserwowanego w ośrodku porowatym o skończonej długości L, w zamkniętym obiegu zawiesiny. Taki obieg pozwala wielokrotnie wykorzystać tę samą zawiesinę, którą po przepływie przez ośrodek powtórnie zatłacza się do tego ośrodka. Rozważania teoretyczne zostały przedstawione w oparciu o układ równań bilansu--transportu (1) i kinetyki procesu kolmatacji (2) z warunkami początkowo-brzegowymi (3),(4), (5). Uzyskane funkcje określają rozkład - w przestrzeni x i czasie t trwającego zjawiska - koncentracji przepływającej zawiesiny N(x,t) (17) i porowatości [epsilon](x,t) (18). Opierając się na równaniu ruchu postaci (20) wyznaczono rozkład ciśnienia w ośrodku porowatym (23).
EN
In this paper we assume that the flow of suspension is accompanied by the exchange of solid particles from a liquid medium into a solid one. We take it as true that both the porous medium and the suspension are characterised by the homogeneity of colmatage properties. We also assume that the flow proceeds along the generating lines parallel to the [chi] axis. The phenomenon under discussion will be described with a suitable system of partial differential equations which consists of the balance-transport equation (1), that of the process kinetics (2) and the equation of motion (6) with boundary conditions (3). Solution of this system results in determining the distribution function of colmatage mass arrested P ([chi], t), that of porosity [epsilon] ([chi], t), that of pressure h ([chi], t) and the function of unitary flow discharge q (t). Basing on the theory, on computation procedures suitably chosen and on an experiment determining the distribution of pressure h ([chi], t) and the discharge of flow q (t) we find out all colmatage parameters. Their knowledge using theoretical description helps to determine every funtion above mentioned and, additionally, a function describing the changeability of medium permeability k([chi],t) when the exchange of mass from a liquid medium into a solid one proceeds during the flows, i.e. when the flows take place with colmatage in the whole space and time of the phenomenon duration. This is obtained using formula (16) for flows without colmatage, or (17) if the flow with colmatage occurs at constant difference of pressures at points [chi] = 0 and [chi] = L.
PL
Niniejsza publikacja dotyczy pewnych aplikacyjnych zagadnień związanych z przepływem zawiesiny przez ośrodki porowate. Zakładamy przy tym, że zawiesina i ośrodek porowaty są jednorodne, a także, że przepływająca przez taki ośrodek porowaty zawiesina osadza się w tym ośrodku. W efekcie tego osadzania dochodzi do ubytku przestrzeni porowej ośrodka, a co za tym idzie do zmniejszenia jego porowatości, a w konsekwencji do przepuszczalności oraz uzmiennienia rozkładu ciśnienia i wydatku przepływu. Jak już wielokrotnie pisaliśmy, taki przebieg zjawiska, z takimi jego konsekwencjami, nazywamy kolmatacją. Zauważmy, że przepływy z kolmatacją mogą być realizowane dwojako: przy stałym wydatku przepływu (i z samoregulującym się rozkładem ciśnienia) lub przy stałej różnicy ciśnień (i zmieniejącym się wydatku). Ze względu na powyższy fakt dzielimy teorię kolmatacji na te dwie podstawowe grupy, tzn. na opis z przepływem wymuszonym charakteryzującym się tym, że q (r) = const i taki, gdzie q = q (t). Ponieważ zjawisku kolmatacji w obu przypadkach towarzyszy ubytek przestrzeni porowej ze wszystkimi podanymi wyżej konsekwencjami, to dla eksperymentatora badającego takie przepływy istotną rzeczą jest określenie rozkładu zatrzymanej masy P([chi],t), rozkładu porowatości [epsilon]([chi],t), jak również rozkładu ciśnienia h([chi],t), zmienności przepuszczalności ośrodka k([chi],t), a zatem i zmienności wydatku przepływu q(t) — w każdej chwili trwającego procesu. Eksperymentator jest w stanie dokonać tego (tak jak wielokrotnie robił to nasz zespół badawczy) poprzez zebranie informacji o rozkładach ciśnienia, poczynając od chwili zerowej trwającego procesu. Następnie poprzez rozebranie kolumny kolmatacyjnej na poszczególne segmenty i wypłukanie osadzonego w przestrzeni porowej kolmatanta, ma możliwość wyznaczenia zatrzymanej w przestrzeni porowej masy tego kolmatanta i wyliczenia jego objętości. Dzięki tym czynnościom mamy informację o rozkładzie ciśnienia h([chi], ti) oraz wyłapanej masy w przestrzeni porowej ośrodka P ([chi], ti). Powtarzając tę czynność dla różnych czasów trwania zjawiska kolmatacji, otrzymuje się ostatecznie rozkład zatrzymanej masy kolmatanta P([chi],t). Ponieważ związek między masą zatrzymaną w ośrodku P([chi],t) a porowatością tego ośrodka jest znany, łatwo określić tym samym i rozkład porowatości w przestrzeni i czasie opisywanego zjawiska: [epsilon]([chi],t). Opisany tok postępowania jest wyjątkowo czasochłonny i pracochłonny. Wymaga wielokrotnego powtarzania eksperymentu z precyzyjnym zachowaniem tych samych warunków. Zauważmy, że jest to trudne lub szczególnie w warunkach naturalnych niemożliwe do zrealizowania. Ze względu na ten fakt i wymogi jakie stoją przed eksperymentatorem opisującym przepływ z kolmatacją pokażemy jak z jednego eksperymentu przebiegającego w określonej przestrzeni i czasie — w którym dokonujemy prostego pomiaru h([chi],t) i q(t) — można kolejno określić wyżej wymienione wielkości fizyczne, a więc: [epsilon]([chi],t). P([chi],t) i k([chi],t). Znajomość tych wielkości jest istotna w wielu dziedzinach gospodarki narodowej. Dokonujemy tego w oparciu o eksperymenty i wielokrotnie zweryfikowaną, zmodyfikowaną i dopracowaną teorię kolmatacji. Polega to na tym, że w oparciu o istniejącą teorię kolmatacji i eksperyment oraz opracowane procedury obliczeniowe dopasowujemy krzywą teoretyczną do posiadanych przebiegów eksperymentalnych. W efekcie tego dopasowania uzyskujemy informację o wszystkich wielkościach charakteryzujących proces kolmatacji. W oparciu o te wielkości i teorię możemy już tym razem wyznaczyć to, co nas interesuje w danym przypadku, a więc: rozkład ciśnienia h([chi],t) w całym obszarze i czasie trwania zjawiska, rozkład porowatości [epsilon]([chi],t), rozkład przepuszczalności k([chi],t) i zmienność wydatku przepływu q(t). Dobre dopasowanie krzywej teoretycznej do eksperymentalnej pozwala więc uniknąć pracochłonnych eksperymentów, wiernie zastępując je opisem teoretycznym. Teoretyczny opis omawianego zjawiska dobrany został w oparciu o układ równań różniczkowych cząstkowych (1), (2) i (6) z warunkami brzeżno-początkowymi (3). Uzyskane z rozwiązania funkcje P([chi],t) i [epsilon]([chi],t) dane są odpowiednio wzorami (4) i (5), a ich przebiegi ilustrują rysunki 3, 4, 7 i 8. Funkcje rozkładu ciśnienia uzyskane z równania ruchu z kolmatacją (7) opisują kolejno te rozkłady dla przepływu realizowanego ze stałym wydatkiem (9) i stałej różnicy ciśnień (10). Funkcja (8) opisuje rozkład ciśnienia przy przepływie cz>stych cieczy nie niosących ze sobą zawiesiny, a więc wtedy, gdy n = 0 lub gdy przepływająca zawiesina nic podlega wymianie z ośrodka ciekłego do stałego, tzn. gdy przepływ jest realizowany bez kolmatacji lub w początkowej chwili procesu z kolmatacją. Ten związek dla jednorodnych ośrodków porowatych lub inaczej mówiąc ta charakterystyka dana jest w postaci „prostej o ujemnym nachyleniu". Należy tu dodać, że ta prosta ilustrująca rozkład ciśnienia dla przepływu bez kolmatacji jest równocześnie charakterystyką określającą rozkład ciśnienia dla przepływu z kolmatacją. ale w chwili t = 0. Kontynuacja zjawiska kolmatacji doprowadza do odstępstwa od tej prostej polegającego na wzroście ciśnienia w każdym punkcie 0 [mniejszy-równy][chi]< L dla przepływu realizowanego ze stałym wydatkiem, i do spadku ciśnienia w każdym punkcie 0 < [chi] < L w przypadku przepływu realizowanego ze stałą różnicą ciśnień. Na rys. 2 i rys. 6 wyraźnie widać, na czym polega to odstępstwo, a opis graficzny tego „odstępstwa" od „prostej o ujemnym nachyleniu" podaje funkcja (9) i (10). Rysunki te ilustrują równocześnie stopień dopasowania przebiegów teoretycznych (linie ciągłe! do danych eksperymentalnych (krzyżyki). Stopień tego dopasowania świadczy o tym, że teoretyczne opisy zjawiska kolmatacji zostały poprawnie sformułowane w swych założeniach. Daje to możliwość określenia z jednego eksperymentu — stosując wymienione procedury — wszystkich interesujących nas wielkości kolmatacyjnych. Poza wymienionymi wyżej funkcjami P([chi],t), [epsilon]([chi], t) i h([chi],t), także — bardzo istotną z aplikacyjnego punktu widzenia — funkcję określającą zmienność w przestrzeni i czasie pod wpływem przebiegającego zjawiska kolmatacji przepuszczalności ośrodka porowatego, jak również zmienność wydatku przepływu. Określenie funkcji przepuszczalności uzyskujemy w oparciu o związek (16) dla przepływów bez kolmatacji albo (17) dla przepływów z kolmatacją realizowanych ze stałym wydatkiem lub (18) dla przepływów z kolmatacją realizowanych przy założeniu stałej różnicy ciśnień w punkcie [chi] = 0 i [chi] = L. Badania teoretyczne dotyczące zjawiska kolmatacji zostały zapoczątkowane w naszym ośrodku badawczym w latach 60. przez Profesora Jerzego Litwiniszyna pracami (Litwiniszyn i inni, 1961 - 1969), których już niestety nie będzie kontynuował. Niezależnie od tego, co byśmy napisali na temat tych prac, nigdy nie jesteśmy w stanie przecenić ich wagi i wartości. Były one pionierskie w skali światowej i otworzyły nowy rozdział w badaniach przepływu zawiesin przez ośrodki porowate. Prowadząc dalej te badania, czujemy się Jego uczniami i jesteśmy świadomi kontynuacji Jego dokonań jako prekursora i nauczyciela.
PL
W publikacji przedstawiono wyniki badań eksperymentalnych rozkładu ciśnienia i zmienności wydatku przepływu w ośrodku porowatym w odniesieniu do zjawiska kolmatacji zachodzącego przy przepływie ze swobodnym zwierciadłem dla różnych: zawiesin i ich koncentracji oraz różnicy ciśnień w punktach x = 0 i x = L
EN
This paper presents the results of experimental research on pressure distribution and variability of flow discharge in porous media with reference to the process of colmatage occuring during flow with a free surface of liquid for different suspensions, and their concentrations as well as for the differences in pressure at points x = 0 and x = L
first rewind previous Strona / 1 next fast forward last
JavaScript jest wyłączony w Twojej przeglądarce internetowej. Włącz go, a następnie odśwież stronę, aby móc w pełni z niej korzystać.