проксгіи]ашкя гарантирует для компонентов скорости третий поря-
док 
[51].
ля временнбй аппроксимации в программе РЕЕГО\ возмо-
жен выбор неявной схемы (у==!) или схемы Кранка-—-Николсона
(}==0,5) с порядком аппроксимации А/ или А.
Так как при у==! вводится численная дисперсия с погрешно-
стью АР == А0%/2ло [см. уравнение (7.5)], это соответственно_ли-
митирует величину А. Поэтому в области времени, в которой не
вводятся скачкообразные изменения на границах, используется_не-
стабильная, но более точная схема Кранка--Николсона, причем
для стабилизации рекомендуется предварительно задавать внут-
ренние граничные условия второго рода для конструкции [51].
После пространственной и временнбй дискретизации составля-
ется система уравнений
[А] (Х (Е+- а6) = 18] (Х () + 1Е (0
© субвектором Х 
Х(Н,, на,_дуль о) в точке 1 [ 
] 
{ 
} —
символы матрицы и вектора, Эта система несимметрична и реша-
ется_посредством так называемого фронтального и 
профильного
алгорнтма [50]. Величины, входящне в матричные и векторные
элементы А, В и Ё, принимаются в РЕЕГ.О\® поэлементно постоян-
ными, Вычисление двойных интегралов 
проводится 
посредством
Зх3-квадратной схемы Гаусса. Программа записана применитель-
но к машине БЭСМ-6.
Посредством программы РЕЕГО® решался пример для харак-
терной области береговой фильтрации к водозабору (рис. 98) при
дискретизации:поля_с. 96.четырехугольными: элементами, 347 узло-
выми точками и общим числом 1168 нензвестных М, Фхль Фр И б
Перенос в слонстых пластах
Как уже отмечалось (см. гл. 4), слонстость пластов весьма су-
щественно влияет на характер и’количественные закономерности
процессов переноса. Моделирование переноса в слоистых пластах
отличается рядом методических особенностей, требующих специ-
ального рассмотрения применительно к различной структуре по-
тока.
Продольный_конвективный перенос с 
поперечной 
диффузией
(дисперсией). В пласте, состоящем из слоев 
© резко различной
проницаемостью 
(например, песчаных и глинистых), продольный
перенос осуществляется конвективным Путем 
в хорошопроницае-
мых слоях путем поперечной диффузии 
и 
в слабопроницаемых.
При этом распределение концентрации мигранта 
в 
каждом ГМ
слое при гомогенном строении пород описывается дифференциаль-
ным уравнением
Щ-Ё—і‹-і-’!},ікд—— 
(7.16)
132