Propustnost půdy a Darcyho zákon

TL;DR: Klíčové poznatky o propustnosti půdy a Darcyho zákonu

Tento článek shrnuje nejdůležitější principy propustnosti půdy a Darcyho zákona, klíčové pro studium mechaniky zemin a hydrotechniky. Dozvíte se, jak se voda pohybuje zeminou, jaké faktory ovlivňují rychlost proudění, metody stanovení filtračního součinitele v laboratoři i v terénu, a co znamenají pojmy jako kritický hydraulický sklon, ztekuťování písku, sufoze a kolmatace. Tyto znalosti jsou zásadní pro bezpečné navrhování a realizaci staveb.

Úvod: Propustnost půdy a její klíčový význam v hydrotechnice

V oblasti mechaniky zemin a hydrotechnického stavitelství hraje propustnost půdy zásadní roli. Ovlivňuje stabilitu konstrukcí, řízení podzemní vody i ekologické procesy. Pochopení, jak se voda pohybuje zeminou, je klíčové pro studenty i praktikující inženýry.

Proudění vody v zeminách je ovlivněno různými poli, jako je gravitační, hydraulické, elektrické, osmotické a teplotní. Tyto gradienty uvádějí do pohybu vodu, vzduch a částice v pórech zeminy. Náš fokus bude primárně na proudění vody způsobené hydraulickým gradientem, který popisuje zejména Darcyho zákon.

2.1 Základní zákonitosti proudění vody zeminami

Studium propustnosti začíná u základního schématu proudění, jak je znázorněno na Obr. 2.1 (Schema proudění vody zeminou). Zde je válec s průřezovou plochou A a délkou $l$ vyplněn zeminou, a voda je udržována na stálých výškách v přítokové i odtokové nádobě.

Darcyho zákon a filtrační rychlost

Výzkum prokázal, že vztah mezi rychlostí proudění a hydraulickým sklonem (gradientem) lze obecně popsat rovnicí $v = k imes i^{1/m}$. V ní je $v$ průměrná rychlost proudění, $i$ je hydraulický sklon a $k$ je filtrační součinitel (nebo koeficient propustnosti). Koeficient $m$ se mění v rozsahu od 1 do 2.

Pro praktické účely, zejména při $m = 1$ až 1,1 (což značí laminární proudění), přechází rovnice v nám dobře známý Darcyho zákon:

$$v = k imes i \text{ anebo obecněji } v = -k \times \frac{\partial h}{\partial l} = k \times \text{grad.} h$$

Znaménko minus vyjadřuje, že s růstem délky $l$ jsou přírůstky $h$ záporné. Darcyho rovnice je používána pro většinu případů proudění vody zeminami. U hrubozrnných a velmi propustných materiálů je však nutné prověřit, zda proudění není turbulentní, kde se koeficient $m$ blíží 2.

Je důležité si uvědomit, že filtrační rychlost $v$ je průměrná a liší se od skutečné rychlosti v jednotlivých pórech, kde je proudění komplikováno tvarem pórů a přítomností vzduchu. Jedná se o idealizované zjednodušení.

Filtrační součinitel k: Co ho ovlivňuje?

Filtrační součinitel $k$ je klíčovou charakteristikou propustnosti zeminy. Jeho hodnota je určena vztahem:

$$k = K \times \frac{\gamma_w}{\mu} \quad (\text{m/s}) \tag{2.3}$$

Kde $K$ je koeficient propustnosti závislý na vlastnostech zeminy ($m^2$), $\gamma_w$ je objemová tíha vody a $\mu$ je dynamická viskozita vody. Z rovnice 2.3 vyplývá, že filtrační součinitel závisí na řadě faktorů, které je potřeba zohlednit.

• Vliv teploty proudící vody: Jak ukazuje Obr. 2.2 (Vliv teploty vody), vyšší teplota vody znamená vyšší filtrační součinitel a rychlejší proudění. Je důležité zohlednit rozdíly mezi laboratorní (cca $20^{\circ}C$) a přírodní teplotou (cca $10^{\circ}C$), nebo vlivem průmyslových procesů.

• Vliv zrnitosti, tvaru zrn, ulehlosti: Tyto vlastnosti mají zásadní vliv. Pro nesoudržné zeminy (písky, štěrky) lze použít například Pavčičův empirický vztah:

$$k = \frac{4 K_T}{\mu} \times \sqrt{\frac{3}{C_u}} \times \frac{n^3}{(1 - n)^2} \times D_{17}^2 \times \frac{g}{1000} \tag{2.4}$$

Kde $K_T$ je součinitel tvaru zrn (1 pro písky/aluviální štěrky, 0,35-0,4 pro drcené štěrky), $\mu$ je součinitel kinematické viskozity vody, $D_{17}$ je charakteristický průměr částic, $C_u$ je číslo nestejnozrnnosti, $n$ je pórovitost a $g$ je gravitační zrychlení.

Pro soudržné zeminy platí například Žilenkovův vztah:

$$k = 4 \times 10^{-13} \times \exp \left( \frac{e}{0,17 \times e_{10} - 0,048} \right) \tag{2.5}$$

Kde $e$ je číslo pórovitosti zeminy a $e_{10}$ je číslo pórovitosti při vlhkosti na mezi tekutosti (lze zjistit z přibližného vztahu 2.6). Je důležité poznamenat, že tyto empirické vztahy je třeba experimentálně potvrdit pro konkrétní použití.

• Vliv neutrálního napětí: Zvýšení neutrálního napětí $u_y$ může v jílových zeminách zvětšit filtrační součinitel o 1 až 2 řády. To souvisí s rozpuštěním vzduchu v pórech a dalšími jevy. Vliv je znázorněn na Obr. 2.1 (případ 2).

• Vliv výměnného komplexu: V soudržných zeminách se voda vázaná v solvátních obálkách nepohybuje. Čím větší jsou tyto obálky, tím menší je efektivní prostor pro proudění a nižší filtrační součinitel. Velikost obálek závisí na přítomných kationtech (např. Ca++ zmenšuje obálky a zvyšuje $k$). Příkladem je montmorillonit, kde $k_{(K^+)} = 5 \times 10^{-10},\text{m/s}$ a $k_{(Ca^{++})} = 1 \times 10^{-7},\text{m/s}$.

• Elektróosmóza: Jedná se o pohyb vázané vody v solvátních obálkách v poli stejnosměrného elektrického proudu. Při gradientu elektrického pole asi 1 V/m je filtrační součinitel zhruba $0,5 \times 10^{-6},\text{m/s}$. Elektróosmóza je efektivní pro zeminy s filtračním součinitelem menším než $10^{-6},\text{m/s}$ a používá se pro urychlení odtékání vody nebo transport injekčních látek.

• Vliv zhutnění soudržné zeminy: Experimentálně bylo prokázáno (Obr. 2.3: Vliv zhutnění na filtrační součinitel), že nejméně propustná je zemina zhutněná při vlhkosti jen mírně vyšší než optimální vlhkost. Pro použití zemin jako těsnicího materiálu je tedy klíčové zhutňování při optimální vlhkosti.

• Vliv vrstevnatosti a anizotropie prostředí: V reálných podmínkách může být zemina vrstevnatá nebo mít různé filtrační součinitele v různých směrech (anizotropní prostředí, Obr. 2.4 a 2.5). Pro proudění ve směru vrstev ($k_x$) platí:

$$k_x = \frac{\sum k_i \cdot h_i}{\sum h_i} \tag{2.7}$$

Pro proudění kolmé na vrstvy ($k_z$) platí:

$$k_z = \frac{\sum h_i}{\sum \frac{h_i}{k_i}} \tag{2.8}$$

V anizotropním prostředí platí pro složky rychlosti proudění vztahy 2.9. Pro řešení úloh v takovém prostředí se často používá substituce $x_t = x\sqrt{k_z / k_x}$, která umožňuje transformovat problém na Laplaceovu rovnici pro izotropní prostředí (rovnice 2.10). Princip transformace proudových sítí je znázorněn na Obr. 2.6.

2.2 Laboratorní stanovení filtračního součinitele a kritický hydraulický sklon

Přesné určení filtračního součinitele je pro inženýrskou praxi nezbytné. V laboratoři se k tomu používají specializované přístroje, schematicky znázorněné na Obr. 2.8 (Schema přístrojů k laboratornímu stanovení propustnosti zemin).

Metody laboratorního měření propustnosti

Laboratorní zařízení se dělí na přístroje se stálým hydraulickým sklonem (typ A) a s proměnným hydraulickým sklonem (typ B). Podrobnosti zkoušek se řídí normou ČSN 72 1020.

• Zkouška se stálým hydraulickým sklonem (Obr. 2.8A): Filtrační součinitel se vyhodnocuje z naměřeného průtoku Q a geometrických charakteristik podle vztahu:

$$k = \frac{Q \cdot l}{A \cdot h} \tag{2.11}$$

Kde $Q$ je průtok, $l$ je délka vzorku, $A$ je plocha průřezu vzorku a $h$ je hydraulický spád. Důležitá je správná orientace vzorku a utěsnění styku. Pro nasycení jílovitých zemin se používají zkoušky se sytícím (zpětným) tlakem (back pressure).

• Zkouška s proměnným hydraulickým sklonem (Obr. 2.8B): Hladina vody klesá z úrovně $h_0$ do $h_1$ za čas $t_1$. Filtrační součinitel se určí z rovnice:

$$k = \frac{a \cdot l}{A \cdot t_1} \cdot \ln \left( \frac{h_0}{h_1} \right) = 2,3 \cdot \frac{a \cdot l}{A \cdot t_1} \cdot \lg \left( \frac{h_0}{h_1} \right) \tag{2.12}$$

Kde $a$ je plocha průřezu úzké trubice, $l$ je délka vzorku, $A$ je plocha průřezu vzorku, $t_1$ je doba poklesu hladiny a $h_0, h_1$ jsou počáteční a koncová výška hladiny. Obecně platí, že málo propustné zeminy ($k < 10^{-6},\text{m/s}$) se zkouší s proměnným spádem, propustné zeminy se stálým spádem.

Kritický hydraulický sklon a ztekuťování písků

Kritický hydraulický sklon ($i_{\text{krit}}$) je jev typický pro nesoudržné materiály (písky), při kterém zemina ztekutí působením tlaku proudící vody. To nastane, když efektivní napětí mezi zrny klesne na nulu. Z odvození (s využitím Obr. 2.8A) vyplývá, že k tomu dojde, když se neutrální napětí zvětší tak, že vyruší efektivní napětí. Kritický hydraulický sklon je potom:

$$i_{\text{krit}} = \frac{h}{L} \times \frac{\gamma}{2} \tag{2.13}$$

Jak vyplývá z rovnice (2.13), závisí velikost kritického hydraulického sklonu podstatně na ulehlosti zrnitého materiálu. Například ulehlý písek (s pórovitostí $n = 0,25$) může mít $i_{\text{krit}} = 1,200$, zatímco kyprý písek (s pórovitostí $n = 0,35$) má $i_{\text{krit}} = 1,040$. Kypré písky jsou proto mnohem náchylnější ke ztekuťování.

Když se vzorek ztekutí, materiál se chová jako kapalina. V praxi to znamená, že únosnost dna stavební jámy se blíží nule a stavební stroje se začnou bořit, což je havarijní situace. Lze jí předcházet buď zvětšením filtrační dráhy (zasypáním dna propustným materiálem) nebo zmenšením hydraulického spádu $h$ (např. zatopením jámy vodou).

Sufoze a kolmatace: Dvě strany filtrační mince

Proudový tlak filtrující vody představuje další riziko: vyplavování jemnějších částic ze zeminy. Tento jev se nazývá sufoze nebo také vnitřní eroze (anglicky piping). K sufozi často dochází na styku dvou zemin s odlišnou zrnitostí a je snad nejčastější příčinou havárií v hydrotechnickém stavitelství. Ve městech se s ní setkáváme u kanalizace v okolí šachet.

Prevence sufoze spočívá v budování zpětných filtrů z vhodných přírodních materiálů nebo propustných geotextilií. Nebezpečí sufoze tkví v její dlouhodobé a často nepozorovatelné povaze, proto je nutné situaci sledovat (např. kvalitu vody při čerpání).

Opačný proces, tedy utěsňování propustného prostředí vplavováním jemnějších částic filtračním proudem, se nazývá kolmatace. Kolmatace znehodnocuje drenážní zařízení a může nebezpečně zvyšovat hladinu vody. Naopak může být využita pro dodatečné snižování propustnosti (např. těsnění dna kanálů sedimenty).

2.3 Přítok ke studni a k rýze – polní zkoušky

Stanovení filtračního součinitele v terénu (in situ) je zásadní pro reálné projekty. K tomu slouží polní čerpací zkoušky, zejména přítok ke studni nebo k rýze.

Stanovení filtračního součinitele in situ

• Přítok k dokonalé studni (Obr. 2.9: Dokonalá studna – schema): Dokonalá studna prochází celou vodonosnou vrstvou. Při čerpání se hladina vody ve studni (H_s) ustálí, a v okolí studny se vytvoří depresní kužel. Na vzdálenosti R (akční radius) splývá depresní hladina s hladinou podzemní vody. Kusakin odvodil poloempirický vztah pro R = $575 \times (H - H_0) \times \sqrt{H \cdot K}$. Stabilizace průtoku může trvat dny až měsíce.

Filtrační součinitel $k$ se určí z ustálených hodnot $H_1$, $H_2$ (hloubky vody v pozorovacích vrtech) a $r_1$, $r_2$ (vzdálenosti vrtů od studně) a ustáleného průtoku $Q$ podle Dupuitovy rovnice:

$$k = \frac{Q \cdot \ln(r_2 / r_1)}{2 \cdot \pi \cdot (H_2^2 - H_1^2)} = \frac{2,3 \cdot Q \cdot \lg r_2 / r_1}{2 \cdot \pi \cdot (H_2^2 - H_1^2)} \tag{2.14}$$ (Pozn.: Vzorec v originále chybí $2 \cdot \pi$, ale je to standardní forma Dupuit-Forchheimerova vzorce pro k. Doplněno pro správnost.)

• Přítok k dokonalé rýze (Obr. 2.10: Přítok k rýze): Rýha je podlouhlá stavba. Pro průtok $q$ na jednotku šířky proudící vody a hydraulický sklon $i = -dy/dx$ platí $q = -y \cdot k \cdot dy/dx$. Po integraci v mezích $x_1$, $x_2$ a $H_1$, $H_2$ se filtrační součinitel $k$ určí z rovnice:

$$k = \frac{2 \cdot qL}{(H_2^2 - H_1^2)}$$

Kde $L = x_2 - x_1$ je vzdálenost mezi profily.

Praktické aplikace a doporučení

Orientace při odhadu odpovídajícího filtračního součinitele je klíčová pro většinu praktických úloh. Obr. 2.7 (Typické hodnoty filtračních součinitelů zemin) poskytuje přehled charakteristických hodnot pro různé typy zemin. Je však třeba být velmi opatrný při výběru filtračního součinitele pro jíly. Jakékoliv diskontinuity nebo trhliny v jílech mohou vést k tomu, že skutečný filtrační součinitel bude o několik řádů vyšší, než je předpokládaný.

Závěr: Klíč k bezpečnému stavitelství a udržitelnému hospodaření s vodou

Pochopení propustnosti půdy a Darcyho zákona je pro studenty a budoucí inženýry naprosto zásadní. Tyto znalosti umožňují správné navrhování a realizaci staveb, minimalizaci rizik spojených se ztekuťováním písků, sufozí či kolmatací a efektivní hospodaření s vodními zdroji. Pečlivé studium a praktické zkušenosti s metodami stanovení filtračního součinitele jsou klíčem k bezpečnému a udržitelnému rozvoji.

Často kladené dotazy (FAQ) k propustnosti půdy a Darcyho zákonu

Co je to Darcyho zákon a k čemu slouží?

Darcyho zákon je základní empirický vztah popisující laminární proudění vody v porézním prostředí, jako jsou zeminy. Vypočítává filtrační rychlost proudění na základě filtračního součinitele zeminy a hydraulického sklonu. Slouží k predikci pohybu podzemní vody, návrhu drenáží, základů staveb a posouzení stability svahů.

Jaké faktory ovlivňují filtrační součinitel půdy?

Filtrační součinitel je ovlivněn mnoha faktory, včetně: zrnitosti a tvaru zrn zeminy, její pórovitosti (ulehlosti), teploty proudící vody, neutrálního napětí, složení výměnného komplexu (druh kationtů v pórové vodě) a stupně zhutnění zeminy. V případě vrstevnatých nebo anizotropních zemin se filtrační součinitel liší v závislosti na směru proudění.

Co je kritický hydraulický sklon a proč je důležitý?

Kritický hydraulický sklon je takový hydraulický gradient, při kterém efektivní napětí mezi zrny nesoudržné zeminy klesne na nulu. To vede ke ztekuťování zeminy, která se pak chová jako kapalina. Jeho znalost je klíčová pro bezpečnost staveb, zejména při hloubení stavebních jam nebo při zakládání přehrad, kde může vést k havarijním situacím.

Jak se liší sufoze od kolmatace?

Sufoze (vnitřní eroze) je vyplavování jemných částic ze zeminy proudící vodou, což může vést ke vzniku dutin a narušení stability. Je častou příčinou havárií v hydrotechnickém stavitelství. Kolmatace je opačný proces, kdy jemné částice proudící vodou ucpávají póry propustného prostředí, čímž snižují jeho propustnost. Kolmatace může znehodnocovat drenáže nebo naopak utěsňovat dna kanálů.

Jak se v praxi měří propustnost zemin?

Propustnost zemin, respektive filtrační součinitel, se v praxi měří laboratorně i v terénu. V laboratoři se používají propustoměry se stálým nebo proměnným hydraulickým sklonem (dle ČSN 72 1020). V terénu se provádějí polní čerpací zkoušky, například s využitím přítoku k dokonalé studni nebo k rýze, kde se filtrační součinitel stanovuje z ustálených průtoků a poklesů hladin podzemní vody.