Разработка алгоритма численного решения задач фильтрации подземных вод в многослойной пористой среде методом конечных элементов. Характерные значения напоров и водопроводимости. Движение грунтовых и напорных вод в первых двух пластах от поверхности земли.
Аннотация к работе
Применение метода конечных элементов к решению задач установившейся фильтрации в многослойных пластахРассмотрим установившееся движение подземных вод в многослойном пласте, состоящем из основного хорошо проницаемого напорного горизонта, покрытого малопроницаемой покровной толщей и подстилаемого снизу слабопроницаемой прослойкой, через которую происходит связь с нижележащим водоносным горизонтом в жестком режиме (рис. При расчетах фильтрации в слоистых водоносных системах обычно используются общие предпосылки перетекания, в которых предполагается, что движение через раздельные относительно малопроницаемые слои происходит только по вертикали, а в хорошо проницаемых слоях - только по горизонтали. Движение подземных вод в многослойной среде с учетом указанных предпосылок описывается следующей системой дифференциальных уравнений [1-3]: где h(x, y), H(x, y) и Z(x, y) - отметки уровня грунтовых вод (УГВ) в верхнем слое и напоров в основном и нижележащем напорных пластах соответственно; Т(x, y)=к(x, y) m(x,y) - водопроводимость основного водоносного горизонта; кв(x, y), к(x, y) и кн= const - коэффициенты фильтрации верхнего, основного и слабопроницаемого пластов соответственно; в(x,y) - поверхность раздела покровного и основного напорного слоев; m(x,y) и мн(x,y) - мощности напорного и слабопроницаемого пластов; f(x, y) - функция инфильтрации; W(x,y) - функция, учитывающая работу эксплуатационных скважин, отбирающих воды основного водоносного горизонта. Образуя начальные приближения h(0) (х, у) и Н(0) (х, у), подставим их в формулы (5) и (10) вместо функций h и Н и решаем уравнения (8) и (9) совместно с краевыми условиями (3) и (4) соответственно.