Шпаргалки к экзамену по динамике подземных вод
СОДЕРЖАНИЕ: 1. Уравнение неразрывности массы жидкости. Для вывода ур-я неразрывности выделим в напорном водоносном пласте кубик dx, dy, dz. Ч-з заднюю грань вытекает масса жидкости = V1. Уравнение неразрывности массы жидкости.
Для вывода ур-я неразрывности выделим в напорном водоносном пласте кубик dx, dy, dz. Ч-з заднюю грань вытекает масса жидкости = Vx dydzdt = M.
V дв-я (массовая) выражается: d(Vx )/dx. На уч-ке dx массовая V получ приращение [d(Vx )/dx*dx]dydxdt. dx – расстояние.
M = Vx dydzdt + [d(Vx )/dx*dx]dydxdt = [Vx + d(Vx )/dx*dx] dydxdt.
Разница между массами жидк, втек и вытек ч/з заданые границы:
M – M = Vx dydzdt - [Vx + d(Vx )/dx*dx] dydxdt = - d(Vx )/dx*dxdydxdt
Проводя аналог вычесл по yz, получим М — разница массы жидк, вошедшей и вышедш со всех сторон =
М = (- d(Vx )/dx - d(Vy )/dy - d(Vz )/dz)dxdydxdt;
М либо накапл в кубике, либо получается за счет уменьш запасов воды в нем. Упругие запасы = М = ndxdydzdt, а скорость их изм во t опр-ся частной произв-й d(М)/dt = d(n)/dt*dxdydxdt. Приравнивая
М = М = (- d(Vx )/dx - d(Vy )/dy – d(Vz )/dz)dxdydxdt = d(n)/dt*dxdydxdt + d(Vx )/dx + d(Vy )/dy + d(Vz )/dz – d(n)/dt = 0 – конечное уравнение неразр массы жидкости. Перечисл выводы могут быть повторены для расчетного эл-та изолир-го напорного пласта мощн m. При этом xy лежат в пл-ти напластования. Расходы жидк ч/з верх и нижн грань z = 0 (изолир пл-т). Ур-е неразрывности:
d(mVx )/dx + d(mVy )/dy + d(mn)/dt = 0
Осн закон-ти (ур состояния, ур движ Дарси, ур неразр) образ систему опред-х ур-й, из к-х можно получить результируещее диф ур-е фильтрации, отраж всю инф о процессе и содерж единст неизв — напор Н, завис от xyzt.
2. Емкостные св-ва г.п.
Водообильность — не общий V воды в порах, а вода кя может быть извлечена водозабором. Способность породы отдавать воду связана с емкостными св-ми. Удаление воды происх 2мя путями: свободн стекание под дейст Fтяж — гравит емкость, отжатие воды под воздейст доп нагрузки — упругая емкость _тк пр-с сжатия носит упр хар-р). Емк св-ва хар-ся и обратным пр-сом — когда порода приним доп V воды.
Гравит емкость — H – изменение уровня воды соотв V вытекшей воды. Для однородного уплотн песка V/wH = V/H = ; w — S попер сечения. — коэф гравит емкости (водоотдачи), отн V вытекшей воды к Vосуш породы. Понятие «осушенности» - условно. После стекания воды в породе ост иммобилизац вода и вода на стыках пор, к-я явл гравит. Аналог коэф гравит водоотдачи при подьеме при подьеме уровня п.в. и заполнения пор водой явл коэф недостатка водонасыщ. Но его числ зн-е меньше, чем , тк в породе сохр защемл воздух. Крупноз пески и гравий — 0,25 — 0,3; мелео и ср пески — 0,15 — 0,2; супеси — 0,05 — 0,1; суглинки и глины — 0,001 — 0,01.
Упругая емкость — пески ост полностью водонасыщ, те на порный водоносн компл отдают воду не осушаясь. 2 пути: за счет расширения воды за счет уменьш в ней гидростат P. За счет сжатия породы.Оба мех-ма опред упругую емк пласта и проявл снижением напоров во время откачки. При откач гидростат Р пониж, уменьш взвеш дейст воды, в породе увел эффектив напряж-я (от веса самой породы), происх уменьш V породы за счет уменьш пор и трещин. Показатель уменьш V пор по мере нагруж породы служит коэф сжим = измен коэф пористости/изм эффект напряж. Ас = /е . Величина, хар-ся изм V жидк в ед-це V породы при единичном изм гидрост напора над коэф упругой емк породы V0 /Vn = *H; * - коэф упр емкости. [*] = [1/м]. При сниж напора к-во жидк в пласте уменьш, те кажд ед V породы отд некот V воды = *Н. * = *m. * - коэф упругой водоотдачи. V0 = *Н. Пески (0,5 — 1)10-4 ; трещин породы 10-5 ; для трещин порист-ти 10-4 .
3. Гидростатический и гидродинамич напор.
Гидростат — механ енерг нек-го V покой жидк-ти с пост опред ее потенц составл и измер-ся работой, кот нужно совершить, чтобы преод дейст 2х полей — гравитации (поднять V жидк на высоту z относ пл-ти сравнения); силы гидростат Р, дейст в V жидк. Н — гидрост напор, hр — пьезометр высота; z — высота. Н = hр + z; Н = Р/sg +z;
hр хар-я долю потенц Е в жидк-ти, связ с действ гидростат давл-я. Z — геом высота полож-я рассматр V жидк относ выбран пл-ти ср-я — уровень М0 . Н — м абс высоты.
Гидродинамич — полная мех Е движ V жидк-ти включ потенц и кинетич Е и опред ур-м Бернулли: Н = Р/sg +z + u2 /2g; u2 /2g – скоростн напор, хар подьем воды, обусл v движ потока жидк. При медл дв-ии пв, скоростной составл пренебрегают, тк при этом доп подьем воды сост не 0,01 мм, тогда полная Е движ пв будет опред-ся гидростат напором.
4. Основные гидродинамические представления ДПВ. Геофильтрация пв. Фильтрационный поток. Область фильтрации. Краевые условия области фильтрации.
Фундаментальные закономерности:
уравнения состояния, отражающие возможный характер изменения физ св-в изучаемой среды во время фильтр пр-са. 1)з-н Гука, отраж зависимость плотности от гидростат давл. 2) компрес ур-е, отр связь пористости с эффект давлением. ас = d/dэ .
Уравнение движения — связь между потерей Е потока и работой сил сопрот — з-н Дарси. v = - kdH/dl.
Усл сохр-е массы жидк — ур неразрывности.
d(mVx )/dx + d(mVy )/dy + d(mn)/dt = 0
Фильтрация — движение вода в насыщ порами пр-ве обусл наличием гидравлич градиента (перепад напоров)
Стр-ра потока пв опред взаимн располож линий тока и линий равного напора. По простр стр-ре выдел 1мерн (между 2мя || реками, дв воды к соверш скв), 2мерн (между не || реками), 3х мерными.
5. З-н Дарси. Усл применимости.
Предпосылки вывода: связь между силами сопрот и изм Е потока (fтр = -m *dH/dl; dH/dl – гидровл град, m = m g.); если движ жидк ламинарное, то завис между v гидр град линейное; из ф-лы Гагено-Пуазейля для сист || капилляр радиусом r. Q = -1/8*ж gr4 ж /ж *dH/dl*N; N = n/r2 к — число капил, - вязкость ж, Q - сум расход; ч/з 1нич капил Q = vd n; vd = -r2/8*ж /* dH/dl;
Дарси проводил опыты по трубкам и определил, что получается прямолин график завис v от I. V = Q/F=kI; F – S попер сеч, Q - расход ж, I - перепад напоров к длине, k - коэф пропорц, завис от грунта и жидк; vф — отн расхода к попер сечен.
З-н Дарси - v = - kdH/dl. Можно мпроект на оси, если св-ва среды в разл напр различны.
Общий вид — v = - k0 ж g/ж * dH/dl. k0 - коэф прониц, зависит только от св-в среды. k = k0 ж g/ж ;
Применим только для ламин дв-я, ограничения:
верх предел — переход ламин в турбулент, сильнозакарст породы, широкие трещины, вход воды из гп в скв, число Рейнальдса Re = 7-9ед.
Ниж предел — слабопрон глин породы, где движ своб воды может происходить только при наличии высокого напора.
6. Разделение потоков пв по усл залегания и гидравл сост. Расход фильтр потока. Скорость фильтрации. Гидростат напор фильтр потока.
Стр-ра потока пв опред взаимн располож линий тока и линий равного напора. По простр стр-ре выдел 1мерн (между 2мя || реками, дв воды к соверш скв), 2мерн (между не || реками), 3х мерными.
Расход Q = kIF. V фильтр = V = Q/F=kI = - kdH/dl.
Гидростат — механ энерг нек-го V покой жидк-ти с пост опред ее потенц составл и измер-ся работой, кот нужно совершить, чтобы преод дейст 2х полей — гравитации (поднять V жидк на высоту z относ пл-ти сравнения); силы гидростат Р, дейст в V жидк. Н — гидрост напор, hр — пьезометр высота; z — высота. Н = hр + z; Н = Р/sg +z;
hр хар-я долю потенц Е в жидк-ти, связ с действ гидростат давл-я. Z — геом высота полож-я рассматр V жидк относ выбран пл-ти ср-я — уровень М0 . Н — м абс высоты.
7. Диф ур упругого режима фильтрации
Более общая модель, к-я учитывает сжимаемость пласта и воды: dV/V = 1/E*dp (коэф пористости); ас (коэф компрессии) = -dt/d (изм эфект напряж).
Счит что из-за пониж v фильтр и малой простр изменч-ти плотности будем иметь: d(vx )/dx = dvx /dx + [vxd/dx]можно пренебр из-за малой изм воды ~ dvx /dx;
Ур-е сост находит отражение в зависимости, включ коэф упругой емкости пласта: dV0/Vn = *dH [30]; V0 = nVn; Vn = V0/n; dV0/(V0/n) = ndV0/V0; [31]; Масса воды в единичн V породы: M0 = n; Выразим 31 ч/з массу: ndM0/M0 = ndn/n = dn/ [32]; Из 30: dn/ = *dH [33]; * - упр емкость [1/м]. Из 33 dn = *dH [34]; Если 34/t: dn/dt = *dH/dt [35];
Ур-е неразрывности: dVx /dx + dVy /dy + dVz /dz + d(n)/dt = 0; 35 в ур-е неразр и /:
dVx /dx + dVy /dy + dVz /dz + dH/dt = 0; С учетом з-на Дарси: -d/dx*(KxdH/dx) - d/dy*(KydH/dy) - d/dz*(KzdH/dz) + dH/dt = 0;
Для однор изотропной среды (Kx=Ky=Kz): K(d2 H/dx2 + d2 H/dy2 + d2 H/dz2 ) = dH/dt; d2H = /K*dHdt [36]; d2H = 1/a*dH/dt; [37] – ур-е Фурье; a* = K/ – коэф пьезопров [м2 /сут];
В отл от жесткого режима ур-е упругого режима опред-ся 3мя простр коорд и временем, и движ не стационарное.
Физически — по мере уменьш напоров во t в водоносн пласте постеп срабатыв его упругие запасы. Высвоб V воды вклад в общий баланс фильтр потока. Реакция от возмущ напоров на гр-це или в какой-то обл водоносн пл-та распростр постеп, и v распр тем выше, чем выше прониц и ниже упругоемкость. Отсюда коэф пьезопров-ти а* явл показат v изм напора (гидростат давл) в пласте.
8. Диф ур жесткого режима фильтрации.
Основываясь на ур неразрывности делаем допущение что порода и вода не сжимаемы: = А, n = В. А и В = const не завис от t. Тогда в ур неразр пропадает временная проиизводная d(n)/dt = 0; Можно считать что интенсивность пространственной изменчивости плотности очень мала: vx d/dx dvx /dx. Тогда из ур неразравн получим: dvx /dx + dvy /dy + dvz /dz = 0;
Подст з-н Дарси: d/dx*(Kx dH/dx) + d/dy*(Ky dH/dy) + d/dz*(Kz dH/dz)=0; Kx Ky Kz – зн-е коэф фильтр по напр осн осей анизотропии. Для однор изотропн пласта: d2 H/(dx)2 + d2 H/(dy)2 + d2 H/(dz)2 = 2 H = 0; Kx =Ky =Kz ;
2 H – сумма 2х призв-х, оператор Лапласса для ф-ии напора. Для 2мерного дв-я: 2 H = d2 H/(dx)2 + d2 H/(dy)2 ;[1]; Для радиального: 1/*H2 = 1/r*d/dr*(rdH/dr);
Ур-е жесткого режима [1] не содержит t в явном виде = при пеизменных напорах на гр-х обл фильтр-ии дв-е должно быть стационарным. H = f(xyz) – не завис от t = водоносн пласт не отдает воду и не принимает. Вся вода проходит транзитом от обл питания к обл разгрузки. При этом р-я на возмущ в любой обл пласта мгновенно распр по всему пласту как в абс жесткой физ системе. Данная модель применима когда транз поток значит превыш упругие запасы пласта.
9. Упругая водоотдача.
* = *m. * - коэф упругой водоотдачи. V0 = *Н. Пески (0,5 — 1)10-4 ; трещин породы 10-5 ; для трещин порист-ти 10-4 .
Величина, хар-ся изм V жидк в ед-це V породы при единичном изм гидрост напора над коэф упругой емк породы V0 /Vn = *H; * - коэф упр емкости. [*] = [1/м]. При сниж напора к-во жидк в пласте уменьш, те кажд ед V породы отд некот V воды = *Н.