Сравнение подходов решения систем нелинейных уравнений в моделях гидроразрыва пласта

Семинар: Информационно-вычислительные технологии в задачах поддержки принятия решений
Начало заседания: 11:00

Дата выступления: 19 Апрель 2022

Организация: НГУ

Авторы: Кочарина Алена Романовна

Гидроразрыв пласта (ГРП) в настоящее время – один из наиболее эффективных методов интенсификации добычи нефти и газа. Его проведение на российских месторождениях невозможно без поддержки математического моделирования. Большинство моделей ГРП включают в себя уравнения движения жидкости, деформации породы и ее разрушения, что приводит к необходимости решения систем нелинейных уравнений (СНУ), получающихся после их аппроксимации, на каждом шаге распространения трещины. 

Поэтому в области математического моделирования и численных методов является актуальной задача решения СНУ таким итерационным методом, который дает результат за наименьшее время.

В работе на примере модели радиальной трещины проведено сравнение методов решения такой СНУ. Движение жидкости в этой модели осесимметрично и описывается в рамках приближения модели ньютоновской жидкости уравнениями неразрывности и количества движения, которые можно записать в виде одного дифференциального уравнения второго порядка для давления. Это уравнение после аппроксимации запишется в виде СЛАУ с трехдиагональной матрицей с зависящими от раскрытия коэффициентами. СЛАУ, полученная при аппроксимации уравнений упругого равновесия будет, как и в общем случае, полностью заполненной. Из двух полученных систем можно составить одну общую систему уравнений для раскрытия и давления, но она уже будет нелинейной.

Для ее решения были рассмотрены три метода: метод Ньютона, метод релаксации и метод Левенберга-Марквардта. В ходе работы реализованы все три метода, проведено их сравнение по точности и скорости.

Метод Левенберга-Марквардта имеет более широкую область сходимости по сравнению с методом Ньютона, что и обусловливало его применение, однако в рассмотренной задаче всегда оказалось легко выбрать начальное приближение, обеспечивающее сходимость метода  Ньютона. 

Научный руководитель – Лапин В.Н.

 

Семинар будет проведен в смешанном формате. Очное заседание пройдет в  конференц-зале ФИЦ ИВТ (к.513), дистанционно к семинару можно будет подключится по ссылке https://vcs-6.ict.nsc.ru/b/grz-ayy-7ne

Запись семинара (выступление Алены Романовны с 19:00): https://vcs-6.ict.nsc.ru/playback/presentation/2.3/caeaa151fe7743faa3d160ccb904deddc5b7d1cb-1650340299531