=…2 , “%,“*=2ель, № 4, 2015 Гелдиев Х.А., кандидат физикоматематических наук, докторант Международного университета нефти и газа (Туркменистан) ЧИСЛЕННОЕ РЕШЕНИЕ ЗАДАЧИ ИДЕНТИФИКАЦИИ С УЧЁТОМ ПРОНИЦАЕМОСТИ Для нелинейного дифференциального уравнения в частных производных параболического типа , описивающего плоскорадиальное движение несжимаемой жидкости к скважине с учетом изменения проницаемости пласта , приводятся алгоритм численного решения методом прямых, что позволяет использовать в дальнейшем методы теории оптимальных процессов. <...> Для приближенного решения уравнений параболического типа обычно используется как метод сеток так и метод прямых [1]. <...> Последний позволяет непосредственно использовать при решении задач оптимизации ( идентификации) систем, описываемых уравнениями в частных производных, методы, разработанным в теории оптимизации [2] для систем, описываемых обыкновенными дифференциальными уравнениями. <...> Предполагается, что имеется горизонтальный пласт постоянной толщины h* и неограниченной (ограниченной ) протяженности. <...> В пласте пробурена одна скважина, вскрывшая его на всю толщину и имеющая открытый забой. <...> Задача о плоскорадиальной фильтрации жидкости к скважине радиуса rc плуатацию в момент t = 0 с постоянным дебетом Q0 , пущенной в экс, с учетом зависимости проницаемости пласта k от давления p сводится к решению уравнения нелинейного уравнения в частных производных параболического типа [3] = (k0 + αL) 2 r R, L, где L(r,t) – функция Лейбензона двух переменных, t - время, rc 0 t при следующих краевых условиях : L( r , 0) = 0 L( R , t) = 0 = Здесь - толщина пласта. проницаемости при p = pk Уравнение (1) получено в предположении: вязкость =const, коэффициент пористости = const. <...> = An где An = ( k0 +αun ), Bn = Краевое условие на правой границе: uN (t) = 0 На левой границе задается производная: = μ = const, которая заменяется соотношением = . <...> Теперь рассмотрим задачу идентификации, заключающейся в уточнении параметра α в уравнении (1) в процессе <...>