КАРТОЧКА ПРОЕКТА ФУНДАМЕНТАЛЬНЫХ И ПОИСКОВЫХ НАУЧНЫХ ИССЛЕДОВАНИЙ,
ПОДДЕРЖАННОГО РОССИЙСКИМ НАУЧНЫМ ФОНДОМ

Информация подготовлена на основании данных из Информационно-аналитической системы РНФ, содержательная часть представлена в авторской редакции. Все права принадлежат авторам, использование или перепечатка материалов допустима только с предварительного согласия авторов.

 

ОБЩИЕ СВЕДЕНИЯ


Номер 22-11-20031

НазваниеМатематическое модели и методы обработки данных мониторинга потоков парниковых газов болотных экосистем севера Западной Сибири

РуководительПятков Сергей Григорьевич, Доктор физико-математических наук

Организация финансирования, регион федеральное государственное бюджетное образовательное учреждение высшего образования "Югорский государственный университет", Ханты-Мансийский АО - Югра

Период выполнения при поддержке РНФ 2022 г. - 2024 г. 

Конкурс№66 - Конкурс 2022 года «Проведение фундаментальных научных исследований и поисковых научных исследований отдельными научными группами» (региональный конкурс).

Область знания, основной код классификатора 01 - Математика, информатика и науки о системах, 01-219 - Математическое моделирование в науках о Земле и проблемах окружающей среды

Ключевые словаидентификация параметров тепломассопереноса, обратная задача, система конвекции-диффузии, определение потоков, конвективный перенос, парниковые газы, слоистая среда, экосистема, круговорот углерода, нелинейная динамика, фазовый портрет, устойчивость стационарных решений, предельные циклы.

Код ГРНТИ27.31.17, 27.35.45, 27.35.25


 

ИНФОРМАЦИЯ ИЗ ЗАЯВКИ


Аннотация
Проблемы глобального изменения климата вызывают большой интерес в силу их огромной практической значимостью. Круговорот углерода, как один из наиболее интенсивных процессов в биосфере в значительной степени определяет климат планеты. Углерод в форме углекислого газа поглощается растениями и выделяется при дыхании мертвого органического вещества. Глобальный круговорот углерода изучался в международном проекте GCP (Global Carbon Project), опубликовано большое количество статьей и монографий, посвященные этому вопросу. Вместе с тем, представляет интерес исследование круговорота углерода в локальных экосистемах, где значительным элементом экосистемы является болота, а также велико антропогенное воздействие. Последнее обстоятельство, в частности, связано с процессами нефтедобычи, проявляется в сжигании попутного газа, а также применением углекислого газа для интенсификации нефтеотдачи пластов. Указанные обстоятельства являются характерными для региона ХМАО-Югра. Важную роль играет создание и анализ математических моделей, описывающих протекающие в биосфере динамические процессы. Предполагается исследование динамики углерода в болотных экосистемах, в том числе: 1) Аналитическое и численное исследование задач восстановления теплофизических, массообменных характеристик, граничных режимов и потоков парниковых газов в системе почва-атмосфера на основе решения обратных задач. 2) Построение новых численных алгоритмов и их обоснование для систем конвекции-диффузиии, проведение численных экспериментов и описание процессов конвективного и диффузионного переноса газов (метан, углекислый газ), процесса поглощения и также температурных режимов в болотистых почвах с учетом региональных особенностей (Западная Сибирь). 3) Аналитическое и численное исследование математических моделей нелинейной динамики экосистем; нахождение особых точек, исследование устойчивости стационарных решений, построение фазовых портретов. Применительно к тематике, исследование прежде всего связано с неопределенностью значительного количества параметров модели системы (граничных условий, теплоемкости, теплопроводности, коэффициента диффузии и т.д.). Традиционный подход с использованием значений характеристик, полученных из эксперимента, часто приводит к неверным результатам. Поэтому задача определения параметров модели является актуальной. Трудности при решении коэффициентной обратной задачи тепломассообмена или диффузии связаны с тем, что в почвах происходит сложный взаимосвязанный термодинамический процесс: идет миграция воды (порового раствора), осуществляется как конвективный так и диффузионный перенос газов а также их поглощение за счет бактерий (метанотрофы) в почве. Исследования газообмена болотных экосистем особенно актуальны в связи с необходимостью инвентаризации основных наземных источников СО2 и метана, являющихся основными парниковыми газами. При описании процессов фильтрации возникают также и постановки в слоистых средах с различными типами условий сопряжения, что усложняет теоретические исследования и методы численного построения решений. Теоретических результатов о разрешимости и единственности решений подобных обратных задач особенно в нелинейном случае крайне мало, хотя имеется очень большое количество работ, посвященных численным методам определении решения. В целом, можно сказать, что основные результаты связаны с методами численного решения подобных задач, причем многие из них далеко не всегда обоснованы. Почти всегда численные методы основаны на сведении задачи к задаче оптимального управления и минимизации соответствующего функционала, что как правило, требует больших вычислительных возможностей и не всегда приводит к нужному результату. Поэтому актуальными являются вопросы разработки новых численных методов решения подобных задач, их теоретическому обоснованию, а также использование полученных результатов в задачах описания процессов конвективного и диффузионного переноса газов (метан, углекислый газ) в почвах, вычислению потоков процесса парниковых газов на поверхности земли, описанию процесса поглощения газов и также температурных режимов в болотистых почвах с учетом региональных особенностей (Западная Сибирь). Предполагается искать неизвестные функции в виде конечных отрезков ряда по некторому базису в отличие от некоторых известных постановок.

Ожидаемые результаты
а) Новые алгоритмы и численные методы решения обратных задач об определении потоков на границе по точечным или интегральным данным, в первом случае будет рассмотрено два варианта: значения концентраций переносимых веществ заданы на границе области и во втором случае точки расположены внутри области (в этом случае задача становится некорректной по Адамару). Теоретические обоснования предложенных методов, включая в том числе описания классов решений и их свойств. Результаты численных экспериментов. б) Новые алгоритмы и численные методы решения обратных задач об определении граничных режимов (определение коэффициента теплопередачи, термического сопротивления, и т.д. по точечным или интегральным данным). Теоретическое обоснования предложенных методов, включая в том числе описания классов решений и их свойств. Результаты численных экспериментов. в) Новые алгоритмы построения решений задач, восстановления теплофизических и массообменных характеристик по точечным и интегральным данным (коэффициенты диффузии, теплоемкость, и пр.) описания их качественных свойств в задачах результаты численных экспериментов. г) Новые алгоритмы построения решений задач, диффузии и переноса в слоистых средах (в том числе задач по определению граничных режимов), обоснования и результаты численных экспериментов; д) Создание и регистрация комплексов программ, посвященных в первую очередь задачам определения потоков парниковых газов. е) Обоснование новых математических моделей циклов углерода, основанных на нелинейной динамике; ж) Исследование устойчивости стационарных решений, построение фазовых портретов математических моделей круговорота углерода з) Численный расчет динамики круговорота углерода, в целях прогноза его концентрации в экосистеме Предполагаемые теоретические результаты предполагают развитие новых подходов и методов для теоретического исследования широких классов обратных задач, основанных на теории дифференциальных уравнений и систем, что позволит моделировать динамику, прогнозировать концентрации углерода в экосистемах на основе новых знаний и измерениях циклов углерода в экосистемах «растения-почва», а также, изучать последствия антропогенного воздействия на экосистему, что чрезвычайно важно для нефтедобывающего региона. Практическая значимость результатов довольно высокая: на основе полученных результатов возможно построение новых численных алгоритмов и создание комплексов программ для решения конкретных экологических задач, создание основы современной автоматизированной сети мониторинга потоков углерода и инновационной системы обработки и анализа данных; дальнейшая работа в этом направлении может послужить основой для создания сервисов по формирования оценок и прогнозирования углеродной ёмкости территории.


 

ОТЧЁТНЫЕ МАТЕРИАЛЫ


Аннотация результатов, полученных в 2022 году
1. Была рассмотрена обратная задача об определении потока через границу для систем уравнений тепломассопереноса. Неизвестная функция, входящая в граничное условие, ищется в виде линейной комбинации некоторых базисных функций с неизвестными коэффициентами зависящими от времени. Определению вместе с решением подлежат эти неизвестные функции от времени. В качестве данных были рассмотрены стандартные начально краевые условия и данные переопределения, которыми являлись значения решения в некотором наборе точек, лежащих на границе рассматриваемой пространственной области. Была получена теорема существования и единственности решений в пространствах Соболева. При этом Предполагаются выполненными условия гладкости на коэффициенты: требуется непрерывность старших коэффициентов в уравнении и принадлежность младших коэффициентов некоторому пространству $L_q(Q)$. Условия на коэффициенты граничных операторов стандартные. Тоже самое относится к данным функциям, входящим в начальные и граничные условия. Предполагаются выполненными условия согласования данных. Кроме стандартных условий требуется выполнение дополнительных условий гладкости на все данные задачи в некоторой окрестности точек замеров. Условия гладкости и согласования дополняются условием корректности (некоторый определитель, построенный по данным задачи отличен от 0). Был построен новый метод численного решения задачи и проведены численные эксперименты. Теоретические результаты опубликованы: https://cyberleninka.ru/article/n/identifikatsiya-granichnogo-usloviya-v-zadachah-teplomassoperenosa; https://www.mathnet.ru/php/archive.phtml?wshow=paper&jrnid=chfmj&paperid=283&option_lang=rus. Были подготовлены данные результатов замеров парниковых газов (углекислый газ и метан) и описаний метеоусловий в районе болотного массива Мухрино (Ханты-Мансийский автономный округ) для последующего использования. 2. Была рассмотрена обратная задача об коэффициента теплопередачи для параболического уравнения второго порядка. В качестве данных были рассмотрены стандартные начально-краевые условия и данные переопределения, которыми являлись значения решения в некотором наборе точек, лежащих на границе рассматриваемой пространственной области. В пространствах Соболева мы получаем теорему существования и единственности решений. Предполагаются выполненными условия согласования данных и условия гладкости на коэффициенты: требуется непрерывность старших коэффициентов в уравнении и принадлежность младших коэффициентов некоторому пространству $L_q(Q)$. Условия на коэффициенты граничных операторов стандартные. Тоже самое относится к данным функциям, входящим в начальные и граничные условия. Кроме стандартных условий требуется выполнение дополнительных условий гладкости на все данные задачи в некоторой окрестности точек замеров. Условия гладкости и согласования дополняются условием корректности, которое представляет из себя отличие от нуля некоторого определителя, построенного по данным задачи. На этой основе можно построить и новый метод численного решения задачи, который и был построен и проведены численные эксперименты. Теоретические результаты оформлены и приняты печати в журнале "Математические заметки" (т.113, вып. 1, 2023). Работа с результатами численных экспериментов пока только готовится к печати. Результаты работы по определению коэффициента теплопередачи будут использованы в более сложных моделях в задачах определения температурных режимов почв. 3. Была рассмотрена обратная задача об потока через границу для параболического уравнения с независящими от времени коэффициентами. В качестве данных были рассмотрены стандартные начально-краевые условия и данные переопределения (значения решения в некотором наборе точек, лежащих внутри рассматриваемой пространственной области). Такая задача некорректна по Адамару. Используя преобразование Лапласа и асимптотику функций Грина эллиптических задач с параметром мы получили условия на данные, гарантирующие существование и единственность решений. Последнее позволяет обосновать некоторые из численных методов решения этой задачи. Условия на данные в целом являются стандартными, однако, для данных переопределения в приведенной обратной задаче, (здесь мы привели задачу к задаче с однородными начально-краевыми условиями) требуются уже условия принадлежности значений решения в фиксированных точках, где мы задаем данные переопределения, некоторому классу бесконечно дифференцируемых функций, определяемого при помощи преобразования Лапласа. Полученные результаты опубликованы: https://www.researchgate.net/publication/360422150_Existence_and_Uniqueness_Theorems_in_the_Inverse_Problem_of_Recovering_Surface_Fluxes_from_Pointwise_Measurements; https://www.mdpi.com/2227-7390/10/9/1549. К сожалению, методика сведения подобных задач к уравнению такого типа в случае коэффициентов, зависящих от всех переменных, пока не разработана. Этот факт существенен, поскольку в практически х задачах вычисления потоков парниковых газов, коэффициенты турбулентной диффузии зависят от времени и других переменных. Для решения подобных задач численно ранее обычно использовался метод сведения задачи к задаче управления и затем минимизация соответствующего функционала качества. Работ, посвященных численному решению задачи достаточно много в силу большого количества приложений данной математической модели. Но нельзя сказать, что имеются достаточно хорошие методы для ее решения. Нами был построен численный метод решения этой задачи в двумерном случае, основанный на методе конечных разностей и методе конечных элементов. В частности, при этом были использованы условия существования решений, из теоретических результатов. Эта программа будет далее оформлена и зарегистрирована. 4. Мы получили теоретические результаты (теоремы существования и единственности) по определению источников атмосферного загрязнения по данным замеров концентраций в стационарном случае. Результаты приняты к печати в Lobachevskii Journal of Mathematics. 5. Построена модель динамики круговорота углерода в экосистеме «растение-почва» с учетом антропогенного воздействия в регионах интенсивной нефтедобычи. Предполагается, что содержание углекислого газа в атмосфере и океане велико, и не изменяется в результате круговорота углерода в наземной растительности. Комбинация указанные потоков составляют динамические уравнения круговорота углерода, построенные по схеме модели Холлинга-Тэннера. Произведен анализ разработанных на текущий момент математических моделей в области динамики углерода, произведено их сравнение применительно к задачам исследования. Описаны условия устойчивости стационарных решений, и фазовых портретов для математической модели. Построен численный алгоритм нахождения решений и написана соответствующая программа. Теоретические результаты частично опубликованы http://izvestiya.asu.ru/article/view/%282022%294-16 и частично приняты к печати: Lobachevskii Journal of Mathematics, Вестник Югорского государственного университета. Результаты численного исследования модели опубликованы: http://konf.asu.ru/math-altai/?page=sborniki-mak, http://izvestiya.asu.ru/article/view/%282022%294-16

 

Публикации

1. Дюкарев Е.А., Семенов С.П. Численное моделирование биогеохимических циклов углерода в болотных экосистемах Известия АлтГУ. Математика и механика, № 4 (126), с.104-109 (год публикации - 2022) https://doi.org/10.14258/izvasu(2022)4-16

2. М.В. Липатов, С.П. Семенов Численное исследование динамики углерода в болотных экосистемах Cборник: «Математики – Алтайскому краю». [Электронный ресурс]: сборник трудов всероссийской конференции по математике с международным участием. – Барнаул: Алтайский государственный университет, 2022. – 389 с., c.113-116 (год публикации - 2022)

3. Пятков С.Г. IDENTIFICATION OF THERMOPHYSICAL PARAMETERS IN MATHEMATICAL MODELS OF HEAT AND MASS TRANSFER Journal of Computational and Engineering Mathematics, v. 9, no. 2, p. 52-66 (год публикации - 2022) https://doi.org/10.14529/jcem220205

4. Пятков С.Г., Баранчук В.А. ИДЕНТИФИКАЦИЯ ГРАНИЧНОГО УСЛОВИЯ В ЗАДАЧАХ ТЕПЛОМАССОПЕРЕНОСА Челябинский физико-математический журнал, Т. 7, вып. 2. с. 234–253 (год публикации - 2022) https://doi.org/10.47475/2500-0101-2022-17206

5. Пятков С.Г., Шиленков Д.В. Existence and Uniqueness Theorems in the Inverse Problem of Recovering Surface Fluxes from Pointwise Measurements Mathematics, v. 10, 1549 (год публикации - 2022) https://doi.org/10.3390/math10091549

6. Пятков С.Г. ON INVERSE PROBLEMS WITH POINTWISE OVERDETERMINATION FOR MATHEMATICAL MODELS OF HEAT AND MASS TRANSFER Вестник ЮУрГУ. Серия ≪Математическое моделирование и программирование≫, Т. 15, № 3. С. 34–50 (год публикации - 2022) https://doi.org/10.14529/mmp220303


Аннотация результатов, полученных в 2023 году
1. Оформлена и зарегистрирована программа по определению потоков парниковых газов по данным замеров концентраций в двумерном случае. Поток газа находился в виде конечного отрезка ряда с неизвестными коэффициентами, зависящими от времени. Использовалась математическая модель тепломассопереноса (параболическое уравнение второго порядка). Алгоритм основан на методе конечных элементов и методе конечных разностей. Задача является некорректной по Адамару, поэтому использовалась регуляризация по Тихонову. Проведены численные эксперименты, которые показали в том числе и устойчивость алгоритма по отношению к малым случайным возмущениям данных. 2. В квазистационарном случае были получены теоремы существования, единственности и оценки устойчивости решений обратной задачи об определении потока, построен и обоснован численный алгоритм, основанный на построении специальных решений сопряженного уравнения и сведении задачи к системе алгебраических уравнений с ненулевым определителем, использовался метод конечных элементов. Задача в этом случае является корректной по Адамару. Проведены численные эксперименты. В модельных случаях трехмерная квазистационарная задача сведена к двумерной задаче об определении источников, записанных в виде в виде конечного отрезка ряда Фурье. Для этой модели получены условия гарантирующие существование решений, их единственность и устойчивость. Результаты получены в общем случае для эллиптической системы уравений с произвольным заданием точек замеров. 3. Рассмотрены вопросы аналитического (в необходимых случаях) и численного исследования задач по определению потоков по интегральным данным. Данные это есть интегралы от решения по области с весом. Теоретические результаты в первом случае, когда данные есть интегралы по оласти с весом, уже известны. Проведены численные эксперименты. 4. Проведено аналитическое и численное исследование задачи об определения коэффициента теплопередачи как в случае одной среды, так и в случае слоистых сред. В результате теоретических исследований (они были начаты еще в 2022 г.) было показано, что задачи с точечным переопределением (т.е. концентрации переносимого вещества или температура дополнительно задаются в точках пространственной области) задача является корректной по Адамару в случае если точки замеров лежат на границе рассматриваемой области. В противном случае задача корректной не является. В первом случае были получены теоремы существования и единственности решений, а также получены оценки устойчивости решений. Была офрмлена программа по численному определению коэффициента теплопередачи в случае двух пространственных переменных была оформлена, проведены численные эксперименты. Полученный алгоритм довольно быстро сходится и устойчив в случайным возмущениям данных. Создана программа по определению коэффициента теплопередачи, входящего в условие сопряжения типа неидеального контакта в случае двух пространственных переменных в слоистых средах также создана, проведен ряд численных экспериментов. 5. Проведен теоретический анализ обратных задач по определению теплофизических характеристик как с смысле теории так и с точки зрения численных методов и предложены и обоснованы в теоретическом плане некоторые новые подходы, прежде всего связанные с определением неизвестных коэффициентов в виде конечных отрезков рядов по известному базису в слоистой среде. Получены теоремы о разрешимости, единственности решений и оценки устойчивости решений в задаче об определении теплофизический характеристик (теплоемкость, теплопроводность) в случае слоистой среды с условиями сопряжения типа неидеального контакта как в случае точечных так и в случае интегральных условий переопределения. 6. Были проанализированы имеющиеся в литературе результаты по определению скорости поглощения метана в среде. С использованием имеющися алгоритмов (работающих в случае большого количества замеров концентраций) проведены численные эксперименты с использованием различных численных методов (такими как метод внутренней точки (численный алгоритм для решения задач оптимизации с ограничениями), алгоритм последовательного квадратичного программирования и генетический алгоритм). Произведен сравнительный анализ результатов. В случае если число замеров недостаточно, т.е. концентрация метана замерена в небольшом количестве точек, коэффициент поглощения строился в виде конечного отрезка ряда (в нашем случае в виде разложения по базисным кусочно-линейным функциям метода конечных элементов). Была предложена другая на наш взгляд более совершенная схема решения задачи, которая сводилась к некоторому нелинейному уравнению и затем решалась численно. 7. В результате проведенных исследований разработана оригинальная математическая модель динамики углерода в болотных экосистемах, учитывающая влияние таких важных факторов, как температура, интенсивность фотосинтеза, скорости деструкции органического вещества. Модель позволяет количественно оценить пространственно-временную динамику потоков углерода между растительным покровом, торфяной залежью и атмосферой. Разработанный инструментарий может быть эффективно применен для расчетно-прогностических оценок эмиссии парниковых газов в зависимости от климатических условий. Получены количественные характеристики потоков метана и углекислого газа имеют в целях уточнения вклада болотных экосистем Западной Сибири в глобальный углеродный бюджет и понимания механизмов регуляции климатических процессов. Разработанная модель предоставляет эффективный инструментарий для исследования чувствительности углеродного обмена к изменениям климатических характеристик, таких как температура воздуха и влагообеспеченность. В разработанной математической модели динамики углерода особое внимание уделено описанию температурных зависимостей скоростей ключевых процессов биотрансформации органического вещества, таких как фотосинтез, автотрофное и гетеротрофное дыхание растительности и микроорганизмов. Использована обобщенная нелинейная параметризация, позволяющая учесть наличие минимальной, оптимальной и максимальной температур для каждого из процессов. Количественные параметры функций определены на основе литературных данных для типичных представителей растительных сообществ болот Западной Сибири. Реализация температурной зависимости скоростей биохимических превращений позволяет моделировать наблюдаемую сезонную и многолетнюю динамику интенсивности фотосинтеза, дыхания и чистой первичной продукции болотных экосистем в зависимости от динамики температурных условий. Проведенный анализ показал существенное влияние динамики температуры на баланс углерода в моделируемых болотных экосистемах Западной Сибири. Реализованный в разработанной модели подход к описанию температурной зависимости биогеохимических процессов является одним из ключевых компонентов, определяющих прогностические возможности модели для оценки отклика углеродного цикла на изменяющиеся климатические условия. Описаны результаты применения полученного программного продукта к известным реальным данным.

 

Публикации

1. Бананчук В.А., Пятков С.Г. On Some Inverse Problems of Recovering Sources in Stationary Convection-Diffusion Models Lobachevskii Journal of Mathematics, Vol. 44, No. 3, pp. 1111–1118 (год публикации - 2023) https://doi.org/10.1134/S1995080223030071

2. Пятков С.Г., Баранчук В.А. Определение коэффициента теплопередачи в математических моделях тепломассопереноса Математические заметки, Т. 113, в. 1. c. 90-108 (год публикации - 2023) https://doi.org/10.4213/mzm13573

3. Пятков С.Г., Баранчук В.А. Determination of the Heat Transfer Coefficient in Mathematical Models of Heat and Mass Transfer Mathematical Notes, Vol. 113, No. 1, pp. 93–108 (год публикации - 2023) https://doi.org/10.1134/S0001434623010108

4. Семенов С.П., Дюкарев Е.А., Ташкин А.О. Biogeochemical Carbon Cycles Numerical Modeling in Wetland Ecosystems Lobachevskii Journal of Mathematics, Vol. 44, No. 3, pp. 1223–1228 (год публикации - 2023) https://doi.org/10.1134/S1995080223030307

5. Семенов С.П., Ташкин А.О. Анализ моделей углеродного цикла применительно к исследованию болотных экосистем Западной Сибири Вестник Югорского государственного университета, Т. 18, №4. С.145-152 (год публикации - 2023) https://doi.org/10.18822/byusu202204145-152

6. Шергин С.Н., Пятков С.Г. Recovering of the heat transfer coefficient from the temperature measurements Вестник ЮУрГУ. Серия «Математическое моделирование и программирование», Т. 16, № 3. С. 51–64 (год публикации - 2023) https://doi.org/10.14529/mmp230304

7. Пятков С.Г., Сафонов Е.И. Программа численного определения поверхностных потоков на границе области по точечным измерениям ФЕДЕРАЛЬНАЯ СЛУЖБА ПО ИНТЕЛЛЕКТУАЛЬНОЙ СОБСТВЕННОСТИ (РОСПАТЕНТ), Свидетельство о регистрации программ для ЭВМ № 2023680110 от 26 сентября 2023 г. (год публикации - 2023)

8. Пятков С.Г., Сафонов Е.И. Программа численного определения поверхностных потоков на границе области по точечным измерениям -, Свидетельство о регистрации программ для ЭВМ № 2023680110 от 26 сентября 2023 г (год публикации - )