время к исключительно широкому и эффективному их
внедрению в практику гидрогеологических исследований
(примеры такого рода приведены в разделе 8.3); методи-
ческие аспекты моделирования более подробно рассмот-
рены, например, в работах [7, 14] . Вместе стем, в послед-
ние десятилетия в гидрогеологии наиболее широко стало
использоваться математическое моделирование на ЭВМ
[19, 20, 40, 48]. 
`
4.3.3. Исходные представления о схемах численного
моделирования нестационарной фильтрации 
На ЭВМ
Рассмотрим простейшее уравнение одномерной филь-
трации (4.1) в конечно-разностной форме:
НН *НЁ+1 -2 Ніц Ні,‹"*'1
=а 
_
А 
(Ах)? 
(4.75)
Пусть для исходного уравнения заданы граничные ус-
ловия: Н (0, Р = Н, = сопы!, Н 1 , В = Н, = соп$! и
начальное условие Н (х, 0) = И(х) (рис. 4.8). Йа выбран-
ной конечно-разностной сётке краевые условия запишут-
ся в виде 
`
Нё=Н, Н:=Н, Н?=ЙН, 
(для#>0),
где нижний индекс { = 0 отвечает левой границе, а 7-М —
правой границе (М = // А х); верхний индекс & = 0 отве-
чает начальному моменту времени.
Перепишем уравнение (4.75) в виде
+ 
к 
@" АЕ [ оу 
Кр 
оруй
НЁ 1 :.-:Ні +*(*°Б*.“х“5*2 (Ні+1 — 2 Ні +Ні_1) , (4_75&)
где в правой части стоят лишь значения напоров на &-ом
временном слое. Положим & =0 и 7= 1, тогда из уравнения
(4.75) получаем формулу для расчета напора Н 11 в первой
узловой точке (т.е. для х = 1 -Ах) на первом временном
слое (т.е. на момент #= 1-Ар:
239