Politechnika Wrocławska Wrocław, dn. 5. 12. 2008

Wydział Geoinżynierii

Górnictwa i Geologii

Studia Stacjonarne Inżynierskie

Rok II, rok akad. 2008/2009

Grupa: PTN godz. 8:15-9:45

Hydrogeologia

LABORATORIUM

Parametry filtracji nielaminarnej

Ćwiczenie wykonał:

Bartłomiej Bilski

Tomasz Chmiel

Damian Mrozek

Iwona Skoczylas

Paweł Strzałkowski

Adam Śpiewak

Konstanty Zdrojewski

I. Wstęp teoretyczny.

Filtracja jest bardzo powolne przesączanie się wody przez system kanalików i porów. Ruch ten można podzielić na dwa typy: filtrację laminarną czyli warstwową, lub ruch burzliwy inaczej też zwany turbulentny. Ruch laminarny polega na równoległym poruszaniu się cząsteczek wody względem siebie i osi przewodu. Występuje on przy niewielkich prędkościach przepływu cieczy, zaś w miarę wzrostu prędkości i przekroczeniu pewnych granicznych wartości, takich jak graniczny spadek hydrauliczny i prędkość graniczna, pojawia się ruch burzliwy. W ruchu tym oprócz ruchu postępowego możemy obserwować także ruch wirowy cząsteczek wody i poprzeczny do kierunku przepływu. Dla przepływu laminarnego obowiązuje prawo Darcy zgodnie z którym, prędkość filtracji układa się w zależności od spadku hydraulicznego.

, [m/s], gdzie:

v - prędkość filtracji [m/s],

k - współczynnik filtracji [m/s],

J - spadek hydrauliczny [-].

Prawo to jednak obowiązuje tylko dla przepływu laminarnego i po przekroczeniu Jgr i vgr zależność prędkości filtracji od spadku hydraulicznego przestaje być liniowa, a wykładnicza. Przy zwiększaniu spadków hydraulicznych zwiększa się przewaga ruchu turbulentnego nad ruchem laminarnym. Zjawisko ruchu turbulentnego w hydrogeologii nazywano fluacją lub filtracją turbulentną. Zależność natomiast między prędkością filtracji a spadkiem hydraulicznym w tym ruchu opisuje między innymi wzór Chezy-Krasnopolskiego:

A f •
, [m/s], gdzie:

Af - współczynnik fluacji [m/s],

n - wykładnik potęgowy [-].

Istnieję jeszcze ruch mieszany. Jest to ruch, w którym mamy do czynienia z ruchem turbulentnym i ruchem laminarnym naraz. W części kanalików o najmniejszym przekroju występuję filtracja laminarna, natomiast w kanalikach o dużych przekrojach filtracja turbulentna. Do przedstawienia tego ruchu istnieje wiele wzorów. Jednym z nich jest formuła Smrekera-Missbacha:

, [m/s], gdzie:

A - współczynnik filtracji nielaminarnej [m/s],

n - wykładnik potęgowy [-].

Wykładnik potęgowy zawarty się między 1 a 2. Jeżeli przyjmuję on wartość 1 to mamy przypadek filtracji laminarnej i współczynnik A w formule Smrekera-Missbacha jest zastępowany współczynnikiem k, który obowiązuję w prawie Darcy. Jeżeli natomiast wykładnik potęgowy jest równy 2 to formuła Smrekera-Missbacha zastępowana jest wzorem Chezy-Krasnopolskiego, który opisuje filtrację burzliwą. Wówczas współczynnik A = Af. Jeżeli natomiast wykładnik znajduję się pomiędzy 1 a 2 to ruch podlega filtracji mieszanej. Zatem istnieję taka granica dla danego gruntu, przy której prędkość filtracji jest opisana dwoma prawami, prawem Darcy i prawem przepływu nielaminarnego. Wówczas dla tego punktu są spełnione następujące równania:

vgr  k •Jgr

i

vgr = A•
gr

Porównując oba równania otrzymujemy:

Jgr = (
)

i

vgr = k•(
)

2. Opis badań

Celem ćwiczenia było wyznaczenie dla wybranego gruntu granic stosowalności prawa

Darcy tj. granicznego spadku hydraulicznego Jgr, granicznej prędkości vgr oraz parametrów przepływu nielaminarnego A i n. Pomiary przeprowadzono na stanowisku, które składa się z cylindra pomiarowego. W dolnej części ma on perforowane dno i siatkę filtracyjną. Dolna część cylindra łączy się z komorą mającą połączenie z piezometrem i zaworem, przez który następuje wypływ przefiltrowanej wody. Górna część cylindra połączona jest natomiast z górną komorą, do której podłączony jest drugi piezometr, przez który doprowadza się wodę do urządzenia. Próbkę badanego gruntu umieszczono w cylindrze pomiarowym. W tym celu wsypano grunt małymi porcjami i zagęszczano. Następnie przez próbkę przepuszcza się wodę w celu usunięcia z porów powietrza, dlatego wodę doprowadzono bardzo wolno (ok. 20 minut) w kierunku od spodu. Potem próbkę obciążono, następnie podłączono pizometry i odpowietrzono cały układ hydrauliczny. Przed wykonywania ćwiczenia zmierzono średnicę d i wysokość l próbki. Badanie polegało na wytworzeniu różnicy wysokości hydraulicznych ΔH w komorach wodnych i pomiarze wydatku Q przepływającej przez próbkę przy danej różnicy wysokości hydraulicznej. Wydatek wody wyznacza się za pomocą wzoru, który określa się mierząc objętość wypływającej wody V w czasie t.

[
]

Dla każdej różnicy wysokości przeprowadzono 2 pomiary wydatku, jeśli zaś pomiary te znacznie się różniły od siebie, przeprowadzono dodatkowy pomiar. Badanie przeprowadzono przy różnych spadkach hydraulicznych J. Spadek hydrauliczny ustalono za pomocą odpowiedniej regulacji dolnego zaworu. Różnicę ΔH odczytywano z podziałki, połączonej do obu kolumn, spadki zaś obliczono za pomocą wzoru:

[-], gdzie:

ΔH - różnica wysokości hydraulicznych [m],

H1 - wysokość hydrauliczna na powierzchni wpływu wody do próbki [m],

H2 - wysokość hydrauliczna na powierzchni wypływu wody z próbki [m],

l - długość drogi filtracji [m].

Prędkość filtracji wyraża się wzorem:

gdzie:

V - objętość przepływającej przez próbkę wody w czasie t [cm3],

F - powierzchnia przekroju próbki prostopadłego do kierunku filtracji [cm2],

d - średnica próbki [cm],

t - czas [s].

Wartością spadków hydraulicznych i prędkości filtracji towarzyszą błędy wynikające z dokładności instrumentów pomiarowych i specyfiki przepływu burzliwego. Niepewności pomiarowe ΔJ i Δv określono metodą różniczki zupełnej za pomocą wzorów:

, gdzie:

ΔV , Δl; Δd; Δt oznaczają błędy wynikające z dokładności mierzonych wielkości.

3. Wyniki badań.

Granice stosowania prawa Darcy i parametry przepływu nielaminarnego określono dla piasku gruboziarnistego. Badanie przeprowadzono na próbce o średnicy 70,0 mm i wysokości37 mm w temperaturze 22 o C. Dla małych spadków hydraulicznych przyjęto czas pomiaru 60 s i zmniejszano go dla coraz większych spadków do 10 s. Czas mierzono stoperem z dokładnością 0,1 s. Największy uzyskany spadek hydrauliczny przy całkowicie otwartym zaworze wyniósł 10. Przy określaniu spadku hydraulicznego jako wartości H1 i H2 przyjmowano poziomy wody w piezometrach odczytywane na skali milimetrowej. W przypadku wahań zwierciadła wody występujących po przekroczeniu granicy przepływu laminarnego przyjmowano średnią wartość uwzględniając amplitudę wahań. Jeżeli wahania nie były widoczne jako błąd oznaczania ΔH1 i ΔH2 przyjmowano 1 mm. W przypadku wahania poziomu wody w piezometrze jako błąd oznaczenia przyjmowano amplitudę wahań. Pomiar objętości wykonywano menzurkami o różnej pojemności 50, 100 i 250cm³. Dokładność podziałek dla menzurek o pojemności 50 i 100 cm³ wynosiła 1 cm³, a dla menzurki o pojemności 250 cm³ - 2 cm³. Jako błąd oznaczenia objętości ΔV przyjmowano maksymalną różnicę objętości przepływającej wody przy tym samym spadku hydraulicznym nie mniejszą jednak niż objętość wynikająca z podziałki. Jest to związane z pulsacyjnym przepływem przy ruchu nielaminarnym. Pozostałe błędy oznaczeń wynikały z dokładności używanych instrumentów pomiarowych tj. Δl - 1,0 mm; Δd - 0,1 mm; Δt - 0,1 s.

Wyniki badań zestawiono w tabeli nr 1.

h1 [cm]	h2 [cm]	Δh [cm]	J [-]	ΔJ [-]	V [cm3]	Vśr [cm3]	t [s]	tśr [s]	Q [cm3/s]	Qśr [cm3/s]	v	Δv
166,9	166,1	0,8	0,2162	0,0599	28	29	60	60	0,466667	0,483333	0,0126	1,5591
															30		60		0,5
					166,7		164,8		1,9				0,5135	0,0679	87	88	60	60	1,45	1,466667	0,0381	1,5592
															89		60		1,483333
					165,7		162,7		3				0,8108	0,076	77	76,5	30	30	2,566667	2,55	0,0663	0,7797
															76		30		2,533333
					164,9		161		3,9				1,0541	0,0825	106	107	30	30	3,533333	3,566667	0,0927	1,5593
															108		30		3,6
					163,7		158,5		5,2				1,4054	0,092	92	92	20	20	4,6	4,6	0,1195	1,0397
															92		20		4,6
					164		158,2		5,8				1,5676	0,0964	112	112	20	20	5,6	5,6	0,1455	1,0398
															112		20		5,6
					162,5		154,5		8				2,1622	0,1125	118	118	15	15	7,866667	7,866667	0,2044	0,7801
															118		15		7,866667
					160		150,4		9,6				2,5946	0,1242	142	142	15	15	9,466667	9,466667	0,246	0,7803
															142		15		9,466667
					157		145		12				3,2432	0,1417	126	126	10	10	12,6	12,6	0,3274	0,5208
															126		10		12,6
					154,8		141,2		13,6				3,6757	0,1534	140	141	10	10	14	14,1	0,3664	0,5209
															142		10		14,2
					151,9		136,2		15,7				4,2432	0,1687	170	170	10	10	17	17	0,4417	0,5212
															170		10		17
					147,7		129,9		17,8				4,8108	0,1841	192	191	10	10	19,2	19,1	0,4963	0,5214
															190		10		19
					142,5		122,5		20				5,4054	0,2001	216	216	10	10	21,6	21,6	0,5613	0,5216
															216		10		21,6
					134		116,2		17,8				4,8108	0,1841	260	260	10	10	26	26	0,6756	0,522
															260		10		26
					125		98		27				7,2973	0,2513	280	280	10	10	28	28	0,7276	0,5221
															280		10		28
					118		89		29				7,8378	0,2659	315	313,5	10	10	31,5	31,35	0,8146	0,5224
															312		10		31,2
					98		62		36				9,7297	0,317	380	380	10	10	38	38	0,9874	0,523
															380		10		38
					86		47		39				10,541	0,3389	415	415	10	10	41,5	41,5	1,0784	0,5233
															415		10		41,5
					77		35,5		41,5				11,216	0,3572	435	433,5	10	10	43,5	43,35	1,1264	0,5235
															432		10		43,2

Tabela nr 1. Wyniki badań.

---