Математическая модель развития и накопления биомассы деревьев на примере сосны обыкновенной

Размещено на http://www.allbest.ru/

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ РАЗВИТИЯ И НАКОПЛЕНИЯ БИОМАССЫ ДЕРЕВЬЕВ НА ПРИМЕРЕ СОСНЫ ОБЫКНОВЕННОЙ

Математическое моделирование все глубже проникает не только в технику, но и во все сферы человеческой деятельности. Под математическим моделированием в технике понимают адекватную замену исследуемого технического устройства или процесса соответствующей математической моделью и ее последующее изучение методами вычислительной математики с привлечением средств вычислительной техники.[1]

В настоящее время значительный интерес для исследователей представляет моделирование различных биологических процессов, осуществляемое с целью управления ростом и развитием биологических объектов. В области лесного хозяйства разработка математических моделей ведется в рамках популяций или экосистем деревьев, в рамках которых отдельно дерево является элементарным объектом (неделимым, но обладающим набором свойств) [2],[3] и в рамках отдельно стоящих деревьев или их малых групп, в раках которых рассматривается развитие элементов отдельного дерева[4],[5],[6],[7].

При моделировании нестационарных процессов роста биологических объектов необходимо учитывать особенности образования и накопления биомассы, стимулирующие и ограничивающие эти процессы факторы, влияние на эти процессы совокупности внешних факторов, значительно изменяющихся в течение года. Для этого важным является выбор адекватного и оптимального принципа моделирования.

Принципы моделирования

Для моделирования был выбран метод, представленный на рисунке 1, который представляет собой изменённый метод конечных разностей. Сначала осуществляется ввод входных величин, которые представляют собой параметры объекта моделирования и внешние параметры моделирования (или совокупность внешних условий, присущих району, типу леса и т.д. на разных этапах моделирования). Важным моментом является подстановка параметров в расчетную часть, в рамках которой определяется характер влияния (функциональный, случайный) и условия изменения параметров с течением времени.

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

Затем осуществляются запись значений для шага моделирования в таблицу результатов. Производится изменение параметров объекта с учетом протекания процессов в течение шага моделирования и внешних условий в зависимости от их циклических колебаний.

Измененные в течение шага входные параметры подставляются в расчетную часть на следующем шаге моделирования.

Рисунок 1 - Общий принцип моделирования динамических процессов роста

Рассмотрим основные этапы моделирования применительно к модели роста дерева, реализованной в среде Microsoft Office Excel, более подробно.

Описание модели

Ввод входных величин

В начале моделирования вводятся следующие входные величины:

1) период моделирования в годах - продолжительность последовательного изменения состояния объекта;

2) начальные параметры саженца - количество живой и неживой биомассы в кг, расположенной в хвое, сучьях, стволе и корнях саженца;

3) постоянные в рамках моделирования параметры почвы - среднее содержание в ней основных элементов в усвояемой форме (азота, калия, фосфора и серы);

Рисунок 2 - Форма для введения начальных параметров саженца и постоянных в рамках моделирования величин

Рисунок 3 - Таблица для циклически изменяющихся в процессе моделирования внешних факторов

4) циклически изменяющиеся в процессе моделирования внешние факторы - типичные для климатической зоны и географической широты изменения в течение года температуры, дневной инсоляции, продолжительности дня, влажности воздуха и почвы. Таблица параметров составляется для каждого дня года с учетом случайного характера величин (кроме продолжительности дня).

Подстановка условий в расчетную часть

Затем осуществляется подстановка входных величин в зависимости от их характера и условий изменения в течение моделирования.

Значение периода моделирования непосредственно используется в рамках разработанной в среде Visual Basic for Application программы в качестве величины, характеризующей количество циклических изменений состояния объекта.

Начальные параметры саженца подставляются в расчетную схему в качестве текущих значений, характеризующих состояние объекта в начале цикла моделирования, с помощью следующего программного кода:

For j = 1 To 4

Sheets("Рабочий лист").Cells(2 + j, 7).Value = Sheets("Рабочий лист").Cells(2 + j, 3).Value 'присваиваем текущему значению начальные условия живой биомассы'

Sheets("Рабочий лист").Cells(2 + j, 8).Value = Sheets("Рабочий лист").Cells(2 + j, 4).Value 'присваиваем текущему значению начальные условия общей биомассы'

Next j

В расчетной схеме используются прямые ссылки на постоянные в рамках моделирования параметры почвы.

Значения циклически изменяющихся в процессе моделирования внешних факторов подставляются в расчетную схему согласно составленной таблицы в соответствии с текущим значением счетчика циклов моделирования. Один цикл моделирования соответствует одному дню роста, то есть значения изменяющихся внешних факторов повторяются через 365 циклов. Большинство изменяющихся внешних факторов являются случайными величинами, которые распределяются относительно средних значений с некоторым отклонением. Для учета случайного характера воздействия величин используется функция возвращения обратного нормального распределения для расчета текущих значений случайных величин, использующихся для расчетов в рамках цикла моделирования.

Расчет в рамках шага моделирования

При проведении расчетов в рамках шага (цикла) моделирования определяются параметры, характеризующие изменение состояние объекта в течение шага.

В рамках ассимиляции рассматриваются процессы получения энергии деревом и расходования энергии на дыхание. Для этого по формуле (1) определяется площадь кроны дерева ():

,

где - Отношение радиуса к высоте кроны, м/м; - текущая масса хвои, кг;

- плотность конуса хвои (конус без малого конуса), кг/м3; - уменьшение линейных размеров малого конуса в %.

С учетом коэффициента эффективности хвои (с учетом проникания света вглубь кроны), м2/м2 () и коэффициента рабочей зоны хвои, м2/м2 () определяется эффективная площадь кроны () по формуле:

.

По формуле (3) определяется получаемая деревом солнечная энергия () :

,

где - инсоляция ср.дневная, Вт/м2; - продолжительность дня, ч; - коэффициент использования света,%.

При этом при понижении температуры окружающего воздуха () ниже температуры начала жизнедеятельности () значение получаемой деревом солнечной энергии принимается равным 0.

По формуле (4) определяется энергия, необходимая дереву на дыхание () :

,

где - коэффициент влияния температуры на дыхание (кратность увеличения при изменении Т на 10 оС); - коэффициент влияния массы дерева на дыхание, кДж/кг; - текущая масса дерева, кг; - коэффициент влияния вл. воздуха на дыхание, кДж/%; - влажность воздуха, % от макс абсолютной.

Таким образом, разность между получаемой энергией () и затрачиваемой на дыхание () представляет собой энергию, доступную для синтеза биомассы ():

.

Энергия, затрачиваемая деревом на дыхание () определяет количество усваиваемого кислорода ():

,

где - энергия, затрачиваемая при дыхании на кг кислорода, кДж/кг.

В рамках процесса ассимиляции синтез биомассы () лимитируется доступной энергией () и количеством усваиваемого кислорода () в соответствии с зависимостями:

,

,

где - масса образующихся органических соединений, кг; - энергия, необходимая на 1 моль органических соединений, кДж/моль; - кислород, необходимый на 1 моль органических соединений, кг/моль.

Модуль расчета процесса ассимиляции представлен на рисунке 4.

Рисунок 4 - Модуль расчета параметров процесса ассимиляции

При описании процесса корневого питания растения рассматриваются процессы усвоения воды и основных минеральных веществ из почвы.

Для этого по формуле (8) определяется объем питательного кома почвы ():

,

где - текущая масса корней, кг; - удельный объем, занимаемый корнями в почве, м3/кг.

Так же определяется площадь кома (конуса), занимаемого корнями, () по формуле:

,

где - отношение радиуса питания к глубине корней, м/м.

По формуле (10) определяются количества доступных в почве воды и минеральных веществ () :

,

где - содержание в почве воды и минеральных веществ, кг/м3.

А по формуле (11) определяются максимальные коэффициенты усвоения воды и минеральных веществ из почвы ():

.

При расчете количества доступных для питания воды и веществ () производится сравнение коэффициентов усвоения воды и минеральных веществ, полученных по содержанию веществ в почве (), и известных значений коэффициентов для дерева (). При недостаточном содержании веществ в почве расчет проводится по формуле:

,

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

При достаточном содержании веществ в почве расчет проводится по формуле:

.

В рамках процесса корневого питания синтез биомассы лимитируется количеством доступных для питания воды и минеральных веществ () в соответствии с зависимостью:

,

где - вода и питательные вещества, необходимые на 1 моль органических соединений, кг/моль.

Модуль расчета процесса корневого питания представлен на рисунке 5.

В рамках ограничения процесса синтеза биомассы в соответствии с зависимостями (7), (14) определяется лимитирующий фактор, по которому в дальнейшем рассчитывается получаемое количество органических соединений () и биомассы ( ) (Рис.6):



где - количество лимитирующего фактора; - количество лимитирующего фактора, необходимое для синтеза 1 моля органических соединений, 1/моль; - количество молей органических соединений, необходимое для синтеза 1 кг биомассы, кг/моль.

Рисунок 5 - Модуль расчета параметров процесса корневого питания

Распределение биомассы между элементами дерева (хвоя, сучья, ствол, корни) осуществляется в соответствии с выбранным распределением для шага моделирования.



Рисунок 6 - Расчет количества биомассы и ее распределения между элементами дерева

При этом предполагается наличие четырех возможных видов распределения биомассы, при которых вся биомасса расходуется на рост одного элемента дерева (Рис. 7). Выбор сценария роста осуществляется на основе анализа лимитирующих факторов и механических свойств элементов дерева.

Рисунок 7 - Рассматриваемые виды распределения биомассы

На основе анализа лимитирующих факторов определяется процесс (ассимиляционный или корневого питания), который ограничивает развитие дерева. Если развитие ограничивается процессом корневого питания, то для текущего шага моделирования принимается первый вид распределения, при котором вся биомасса идет на развитие корневой системы дерева.

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

Рисунок 8 - Механические ограничения при развитии ассимиляционного аппарата

При этом в модели реализован процесс отмирания элементов дерева через определенные для каждого элемента интервалы времени. Реализованный в модели процесс отмирания элементов заключается в вычитании из текущих значений биомассы значений ее прироста, произошедшего ранее (ровно на период отмирания). При этом так же ведется учет общей биомассы, в которой не учитывается отмирание элементов (кроме хвои). сосна моделирование дерево биомасса

Запись текущих результатов в таблицу

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

For I = 1 To Sheets("Рабочий лист").Cells(22, 7).Value 'присваиваем текущим значениям динамических параемтров значения следующего шага'

For m = 1 To 4

Sheets("Выходные данные").Cells(3 + I, 1 + m).Value = Sheets("Рабочий лист").Cells(2 + m, 7).Value 'записываем текущие значениям в таблицу результатов живой биомассы'

Sheets("Выходные данные").Cells(3 + I, 5 + m).Value = Sheets("Рабочий лист").Cells(2 + m, 8).Value 'записываем текущие значениям в таблицу результатов общей биомассы'

Next m

Следующие операции шага моделирования

Next I

Получаемая в процессе моделирования таблица выходных данных представлена на рисунке 9.

Рисунок 9 - Таблица выходных данных моделирования

Изменение параметров объекта и внешних условий моделирования

В течение шага изменяется текущее состояние объекта - осуществляется прирост биомассы дерева. Для этого в расчетной схеме определяются значения для следующего шага моделирования, а затем они подставляются с помощью приведенного ниже программного кода в ячейки для текущего состояния объекта и используются на следующем шаге моделирования (Рис. 10).

For I = 1 To Sheets("Рабочий лист").Cells(22, 7).Value 'присваиваем текущим значениям динамических параметров значения следующего шага'

Предыдущие операции шага моделирования

For k = 1 To 4

A(k) = Sheets("Рабочий лист").Cells(17 + k, 7).Value 'присваиваем значения следующего шага для Живой биомассы переменным'

B(k) = Sheets("Рабочий лист").Cells(17 + k, 8).Value 'присваиваем значения следующего шага для Общей биомассы переменным'

Next k

For n = 1 To 4

Sheets("Рабочий лист").Cells(2 + n, 7).Value = A(n) 'присваиваем текущим значениям Живой биомассы значения переменных'

Sheets("Рабочий лист").Cells(2 + n, 8).Value = B(n) 'присваиваем текущим значениям Общей биомассы значения переменных'

Next n

Sheets("Рабочий лист").Cells(23, 7).Value = I 'счетчик Циклов'

Next I



Рисунок 10 - Использование значений для следующего состояния в качестве текущих на следующем шаге

При выполнении цикла моделирования изменяется значение счетчика, который используется для определения текущего состояния моделирования и для циклического изменения внешние параметры по описанной выше схеме.

Измененные в течение шага входные параметры подставляются в расчетную часть на следующем шаге моделирования. Затем производится расчет уже с новыми входными величинами. Циклы повторяются, пока счетчик не достигнет требуемого периода моделирования.

Рисунок 11 - Результаты моделирования с помощью разработанной модели

Выводы

1. При моделировании нестационарных процессов роста биологических объектов необходимо учитывать особенности образования и накопления биомассы, стимулирующие и ограничивающие эти процессы факторы, влияние на эти процессы совокупности внешних факторов, значительно изменяющихся в течение года. Для этого важным является выбор адекватного и оптимального принципа моделирования.

2. Разработана математическая модель развития и накопления биомассы дерева, позволяющая установить зависимость между условиями развития дерева и количеством накапливаемое в процессе роста биомассы в различных элементах дерева.

3. Таким образом, общий вид график соответствует S-образной кривой роста [8], что позволяет сделать вывод об адекватности разработанной математической модели и возможности использования ее результатов для прогнозирования роста деревьев.

Работа выполнена в рамках реализации ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009-2013 годы (Соглашение № 14.B37.21.0301).

Библиографический список

1. Зарубин В.С. Математическое моделирование в технике - М.:Изд-во МГТУ им. Н.Э. Баумана, 2003.

2. Карев Г. П. Структурные модели лесных экосистем / Г. П. Карев // Сибирский экологический журнал. - 1999. - №4. - с. 381-396.

3. Карев, Г. П. Структурные модели и динамика древесных популяций: дис... докт. физ.-мат. наук : 03.00.02 / Карев Георгий Петрович - Москва, 1993. - 224 c.

4. Медведев С. С. Физиология растений. - СПб.:Изд-во С.-Петерб. ун-та, 2004.

5. Колобов А.Н. Численно-аналитическое исследование модели роста дерева в условиях конкуренции за свет / А.Н. Колобов // Математическая биология и биоинформатика. - 2012. - Т.7. - №1. - с. 125-138.

6. Омелько А. М. Модель роста деревьев темнохвойных пород на основе L-систем/ А. М. Омелько // Сибирский экологический журнал. - 2006. - №2. - с. 181-188.

7. Омелько А. М. Математическая модель роста дерева в древостое темнохвойных пород/ А. М. Омелько // Биологические исследования на горнотаежной станции. Сб. научн. трудов - 2006. - Вып.10. - с. 86-98.

Аннотация

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ РАЗВИТИЯ И НАКОПЛЕНИЯ БИОМАССЫ ДЕРЕВЬЕВ НА ПРИМЕРЕ СОСНЫ ОБЫКНОВЕННОЙ

Медяков Андрей Андреевич к.т.н.

Онучин Евгений Михайлович к.т.н., доцент

Перетягин Павел Андреевич студент

Поволжский государственный технологический университет, Йошкар-Ола, Россия

В статье приведено описание модели нестационарного роста биологических объектов на примере Сосны обыкновенной. Рассмотрены особенности моделирования роста биологических объектов

Ключевые слова: БИОЛОГИЧЕСКИЕ ОБЪЕКТЫ, АССИМИЛЯЦИОННЫЙ АППАРАТ, ЗОНА ПИТАНИЯ, СОСНА ОБЫКНОВЕННАЯ

MATHEMATICAL MODEL OF THE DEVELOPMENT AND ACCUMULATION OF TREE BIOMASS

Medyakov Andrey Andreevich Cand.Tech.Sci.

Onychin Evgeny Mihailovich Cand.Tech.Sci., assistant professor

Peretyagin Pavel Andreevich student

Volga State University of Technology, Ioshkar-Ola, Russia

The article describes the model of the unsteady growth of biological objects on the example of pine. The features of the growth modeling of biological objects are presented

Keywords: BIOLOGICAL OBJECTS, ASSIMILATION APPARATUS, FEED ZONE, PINE

Размещено на Allbest.ru