КАРТОЧКА ПРОЕКТА ФУНДАМЕНТАЛЬНЫХ И ПОИСКОВЫХ НАУЧНЫХ ИССЛЕДОВАНИЙ,
ПОДДЕРЖАННОГО РОССИЙСКИМ НАУЧНЫМ ФОНДОМ
Информация подготовлена на основании данных из Информационно-аналитической системы РНФ, содержательная часть представлена в авторской редакции. Все права принадлежат авторам, использование или перепечатка материалов допустима только с предварительного согласия авторов.
ОБЩИЕ СВЕДЕНИЯ
Номер проекта 25-11-20026
НазваниеОбратное моделирование в задачах анализа баланса парниковых газов и теплового режима почв Северных регионов
Руководитель Пятков Сергей Григорьевич, Доктор физико-математических наук
Организация финансирования, регион федеральное государственное бюджетное образовательное учреждение высшего образования "Югорский государственный университет" , Ханты-Мансийский АО - Югра
Конкурс №101 - Конкурс 2025 года «Проведение фундаментальных научных исследований и поисковых научных исследований отдельными научными группами» (региональный конкурс)
Область знания, основной код классификатора 01 - Математика, информатика и науки о системах; 01-219 - Математическое моделирование в науках о Земле и проблемах окружающей среды
Ключевые слова идентификация параметров тепломассопереноса, обратная задача, система конвекции-диффузии, тепловые режимы почв, распространение загрязнений (нефть, рассолы) почв в криолитозоне, диффузионный и конвективный перенос, парниковые газы, слоистая среда, экосистема, круговорот углерода, нелинейная динамика, физико-химические процессы в почвах, фазовый портрет, устойчивость стационарных решений, предельные циклы,,
Код ГРНТИ27.35.45, 27.31.17, 27.35.25
ИНФОРМАЦИЯ ИЗ ЗАЯВКИ
Ожидаемые результаты
Ожидаемые результаты:
а) Новые алгоритмы и численные методы решения обратных задач возникающих при описании потоков парниковых газов потоков на границе по интегральным данным, баланса и процессов поглощения и рождения парниковых газов в почвах и водоемах. Теоретические обоснования предложенных методов, включая в том числе описания классов решений и их свойств. Результаты численных экспериментов.
б) Новые алгоритмы построения решений задач восстановления теплофизических и массообменных характеристик по точечным и интегральным данным (коэффициенты диффузии, теплоемкость, и пр.) в задачах описания тепловых режимов грунтов и описания процессов тепломассопереноса при техногенном загрязнении почв, влияния паводковой воды и нефтяного загрязнения на тепломассообменный режим грунта. Описание качественных свойств решений и результаты численных экспериментов. Теоретическое обоснование результатов в том числе и в случае слоистых сред с различными типами условий сопряжения.
г) Создание и регистрация комплексов программ, посвященных в первую очередь задачам описания балансов и потоков парниковых газов и тепловых режимов почв в условиях антропогенного воздействия.
е) Обоснование новых математических моделей циклов углерода, основанных на нелинейной динамике;
ж) Исследование устойчивости стационарных решений, построение фазовых портретов математических моделей круговорота углерода
з) Численный расчет динамики круговорота углерода, в целях прогноза его концентрации в экосистеме
Предполагаемые теоретические результаты предполагают развитие новых подходов и методов для теоретического исследования широких классов обратных задач, основанных на теории дифференциальных уравнений и систем, что позволит моделировать динамику, прогнозировать концентрации углерода в экосистемах на основе новых знаний и измерениях циклов углерода в экосистемах «растения-почва», а также, изучать последствия антропогенного воздействия на экосистему, что чрезвычайно важно для нефтедобывающего региона.
Практическая значимость результатов довольно высокая: на основе полученных результатов возможно построение новых численных алгоритмов и создание комплексов программ для решения конкретных экологических задач, создание основы современной автоматизированной сети мониторинга потоков углерода и инновационной системы обработки и анализа данных; дальнейшая работа в этом направлении может послужить основой для создания сервисов по формирования оценок и прогнозирования углеродной ёмкости территории.
ОТЧЁТНЫЕ МАТЕРИАЛЫ
Аннотация результатов, полученных в 2025 году
1. Исследованы вопросы корректности обратных задач об определении потока на границе области с условиями переопределения в виде интегралов от решения с весом по границе области или ее части. Поток (коэффициент теплопередачи) ищется в виде конечного отрезка ряда Фурье с неизвестными коэффициентами зависящими от времени. Были исследованы вопросы существования и единственности решений и наличие оценок устойчивости, а также разработаны численные алгоритмы построения решения по определению потока. Найденные условия корректности задачи локальные по времени получены в виде отличия от нуля некоторого определителя, построенного по данным задачи. Разработаны алгоритмы численного решения задачи об определении потока как в случае размерности по пространству 2 так и в случае размерности 3, проведены численные эксперименты. В качестве пространственной области, в которой рассматривается уравнение, была взята цилиндрическая области. Области такого типа наиболее часто возникают в приложениях и при численном моделировании. Программы предполагается использовать при численном моделировании в задачах определения удельных потоков парниковых газов из почвы методом обратной задачи.
Алгоритмы основаны на полученных теоретических результатах а также на методе конечных элементов и методе конечных разностей.
2. В качестве модели для расчета концентраций метана в почке рассматирвалось параболическое уравнение второго порядка. Теоретически исследована задача определения концентрации метана в почве с учетом поглощения метана метонотрофами по данным замеров концентраций в точках области. Задача сводится в решения обратной задачи по определению коэффициента поглощения. Это функция является младшим коэффициентом в параболическом уравнении, который ищется в виде конечного отрезка ряда с неизвестными коэффициентами Фурье. В общем случае получены теоремы существования и единственности (в некотором шаре) решений обратной задачи для малых данных. С целью определения средней величины поглощения была рассмотрена задача об определении числового коэффициента перед коэффициентом поглощения. В этом частном случае получена глобальная по времени теорема существования и единственности решений задачи об определении этой постоянной и распределения концентрации газа.
Для всех численных экспериментов был использован метод глобальной оптимизации -дифференциальная эволюция, обеспечивающий устойчивость при поиске решений в задачах с многомерными и негладкими пространствами параметров. Исследование проводилось для одномерного стационарного, одномерного, двумерного и трехмерного нестационарного случаев. Получены и описаны результаты серии численных экспериментов, направленных на решение обратной задачи восстановления параметров функции скорости переноса газа в пористой среде. В одномерном стационарном случае и в одномерном и двумерном нестационарном случае обратная задача решается с высокой точностью. Метод дифференциальной эволюции показывает устойчивость и способность находить глобальный минимум, несмотря на скачкообразный характер оптимизации. В четырёхмерной постановке (время + пространственная переменная) точность заметно снижается, что указывает на необходимость увеличения числа измерений концентрации, усложнения параметризации функции $g(x,t)$ и применения гибридных схем.
3. Исследованы вопросы коррректности обратных задач определения коэффициента теплопередачи на границе раздела двух сред. Коэффициент теплопередачи ищется в виде конечного отрезка ряда Фурье с неизвестными коэффициентами зависящими от времени. В качестве условий сопряжения используются условия неидеального контакта. В качестве условий переопределения используются набор средних – интегралов от решения с весом по пространственной области. Условия переопределения в виде задания набора некоторых средних величин в литературе используются редко. Кроме всего прочего, такие данные могут быть использованы как приближения точечных условий переопределения, для построения приближенного решения обратных задач с точечными условиями переопределения. Получены оптимальные условия существования и единственности решений а также наличие оценок устойчивости. Все рассмотрения проводятся в пространствах Соболева. Доказательство основано на получаемых априорных оценках и теореме о неподвижной точке. Разработан алгоритмы построение решения таких задач основанный на методе конечных элементов и конечных разностей (неявная схема). Проведены численные эксперименты. Описаны результаты численных экспериментов. Они показывают хорошую сходимость алгоритма, его точность и устойчивость при случайном возмущении данных.
4. Для квазилинейного параболического уравнения исследвана рассмотрена задача определения коэффициентов старшей части уравнения, которая нелинейным образом зависит от решения. Задача в частности, возникает при описании теплдовых режимов почв северных территорий. Получены теоремы существования и единственности в случае когда неизвестные функции входящие в нелинейный коэффициент теплопроводности зависят о времени. Работа сейчас в стадии оформления. Рассмотрены интегральные условия переопределения.
5. Разработана новая динамическая математическая модель, учитывающая процессы конвекции и диффузии газов в почвах, поглощение газов растениями, а также влияние температурных режимов на интенсивность этих процессов. Модель основана на физически обоснованном разделении процессов трансформации органического вещества и явно учитывает влияние уровня грунтовых вод на переход между аэробными и анаэробными условиями в почвенном профиле. Модель описывает углеродные потоки в болотных экосистемах с использованием двух основных пулов органического углерода: Live (относительный запас углерода в живой растительной биомассе) и Mort (относительный запас углерода в отмерших органических материалах). Модель включает расширенный метановый цикл, детализирующий процессы метаногенеза в анаэробных зонах и последующее окисление метана (метанотрофию) как в аэробном поверхностном слое почвы, так и в ризосфере растений. Особое внимание уделено механизмам переноса газов, включая молекулярную диффузию, конвективный перенос с почвенными растворами и критически важный растение-опосредованный транспорт метана через аэренхиму, который является доминирующим путем эмиссии CH4 в атмосферу. Ключевым элементом модели является функция контроля уровня грунтовых вод, которая нелинейно регулирует соотношение объемов аэробной и анаэробной зон и тем самым перераспределяет потоки углерода между процессами аэробного дыхания и метаногенеза.
Публикации
1.
Пятков С.Г., Сафонов Е.И., Петров А.А.
INVERSE PROBLEM OF RECOVERING FLUXES FROM INTEGRAL DATA
Вестник Южно-Уральского государственного уни верситета. Серия Математическое моделирование и программирование, vol. 18, no. 4, pp. 32–44 (год публикации - 2025)
10.14529/mmp250404
2.
Тукмачева Ю.А., Пятков С.Г.
Определение параметра поглощения в уравнении конвекции-диффузии
Челябинский Физико-Математический Журнал, Т. 10, № 4. С. 701–714 (год публикации - 2025)
10.47475/2500-0101-2025-10-4-701-714
3.
Тукмачева Ю.А., Пятков С.Г.
Определение коэффициента поглощения по точечным данным
Вестник Югорского государственного университета, Том 21, № 3, С. 54-59, (год публикации - 2025)
10.18822/byusu20250354-59
4.
Солдатов О.А.
ИДЕНТИФИКАЦИЯ КОЭФФИЦИЕНТА ТЕПЛООБМЕНА ПО ГРАНИЧНЫМ ИНТЕГРАЛЬНЫМ ДАННЫМ
Математические заметки СВФУ, Т. 32. № 3. С. 82-94 (год публикации - 2025)
10.25587/2411-9326-2025-3-82-94
5. Пятков С.Г., Сафонов Е.И. Программа численного определения потоков на границе области по интегральным измерениям в двумерном случае Свидетельство о государственной регистрации программ № 2025692351, Заявка № 2025691392 Дата поступления 10 ноября 2025 г. Дата государственной регистрации в Реестре программ для ЭВМ 20 ноября 2025 (год публикации - 2025)