Основные модели роста численности популяций.

Оглавление

Введение. 4

Основные модели роста численности популяций. 4

Демографические таблицы популяций. 18

Расчет компонент сбалансированного техноценоза. 23

Зонирование территории по степени поврежденности экосистем. 27

Расчет рассеивания выбросов в атмосфере по методике ОНД-86. Определение норматива ПДВ и необходимой степени очистки выбросов. 34

Расчет платы за загрязнение атмосферы.. 41

Расчет платы за загрязнение водных объектов. 46

Организация системы мониторинга заданного объекта. 49

Расчет материального баланса веществ при сжигании твердого топлива. 56

Расчет разбавления вредных примесей по методу В.А. Фролова - И.Д. Родзиллера. Определение норматива НДС и необходимой эффективности работы очистных сооружений. 60

Расчет кислородного баланса водоема. Модель Стритера-Фелпса. 64

Список литературы.. 68

Введение

 

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

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

 

 

Занятие 1

КОНТРОЛЬНОЕ ЗАДАНИЕ.

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

Таблица 1.4

ЧИСЛЕННОСТЬ МИРОВОГО НАРОДОНАСЕЛЕНИЯ

Год
Млн. чел

 

Год
Млн. чел

 

2. Определить массу зерна, которую необходимо положить на последнюю клетку шахматной доски, если, начиная с одного зернышка на первой клетке на каждую следующую кладется в два раза большее число зерен. Вес одного зерна принять равным 0.1 г.

 

3. Определить время, за которое потомство одной бактерии, размножающейся делением и весом 10^-12 г. достигнет веса Земного шара - 10²⁸ г., если одно деление происходит с интервалом в 20 минут.

 

4. Экспоненциально растущая популяция эксплуатируется в стационарном режиме, необходимо увеличить снимаемый с нее урожай в «к» раз за время «t». Определить, при каких «t» возможно увеличение урожая в «к» раз.

 

5. Популяция эксплуатируется в режиме на уничтожение.Определить время, за которое численность популяции сократится до 10% от первоначальной, см. табл.1.4.

 

Таблица 1.4

Исходные данные для выполнения контрольного задания по вариантам

Варианты	n0	r	k
		1.0	1.3
		0.8	1.0
		1.2	5.6
		0.3	1.2
		0.6	1.5
		1.3	2.6
		0.1	0.3
		0.7	1.6
		0.08	0.1
		0.7	1.4

 

6. Площадь лесов планеты составляет приблизительно 3.85 млрд. га. С учетом процессов искусственного и естественного возобновления скорость сведения лесов составляет 20 га/мин. Определить время, за которое площадь лесов сократится до 10% от первоначальной, предполагая, что их сведение происходит по экспоненте.

 

7. По данным таблиц хода роста древостоев основных лесообразующих пород начиная с возраста 60 лет и старше определить их популяционные параметры приняв в качестве закона роста J - образную кривую. cм.табл.1.5.

 

8. По данным таблиц хода роста основных лесообразующих пород определить их популяционные параметры приняв в качестве закона роста модифицированный логистический. Использовать данные о ходе роста по запасу. cм.табл.1.5.

 

Таблица 1.5.

ХОД РОСТА ОСНОВНЫХ ЛЕСООБРАЗУЮЩИХ ПОРОД, м³/га.

Возраст,	Номера вариантов и данные для определения популяционных параметров
лет

 

Занятие 2

Занятие 3.

Занятие 4.

Указания к решению задачи

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

По формуле 4.3 рассчитать плотности поврежденной растительности на каждом из расстоянийr_i по двум исходным направлениям

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

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

По формуле 4.10 рассчитать параметр f для каждого из двух исходных градиентов.

По формуле 4.11 рассчитать параметр f для оставшихся шести градиентов.Расчет следует производить на основе ближайших градиентов. Например, для первого варианта расчет северо-восточного градиента правильней производить на основе данных по северному градиенту, а не по юго-западному.

На основе уравнения 4.8 определить значение параметра d. Значения параметра d по искомым градиентам может немного различаться в зависимости от того, какой из двух исходных градиентов был выбран ранее за основу для расчета параметра f по тому или иному градиенту.

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

Из уравнения 4.7 определить статистическую сумму распределения Z для всех градиентов

Из уравнения 4.7 определить плотность поврежденной растительности на каждом выбранном ранее расстоянии по каждому из шести искомых градиентов.

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

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

 

Рис. 4.1. Пример определения размеров площадей зон повреждения

 

 

Занятие 5

Занятие 6.

Работа 7

Работа 8

Занятие 9

Задание 10

Занятие 11

КОНТРОЛЬНОЕ ЗАДАНИЕ

1. При заданных параметрах модели определить X_max, D_max, t_max.

Таблица 11.3.

Исходные данные для выполнения контрольного задания

Вариант	L0, мг/л	D0, мг/л	k1, 1/сутки	k2, 1/сутки	V, м/сутки
1.	5,2	0,50	0,15	0,45
2.	5,6	1,00	0,11	0,42
3.	2,8	0,07	0,13	0,45
4.	1,4	0,12	0,21	0,53
5.	3,5	0,16	0,12	0,51
6.	5,7	0,70	0,18	0,52
7.	2,3	0,02	0,13	0,50
8.	3,7	0,50	0,32	0,52
9.	4,3	0,15	0,17	0,45
10.	5,2	0,10	0,10	0,30

 

2. При заданных постоянных реаэрации (k₁) и скорости окисления (k₂) при­месей (см. таблицу) определить максимально допустимую вели­чину БПК в зоне загрязнения (L_o) такую, чтобы минимальная концентрация кислорода в любом месте реки составляла с=6 мг/л. Положить равновесную концентрацию кислорода в воде при отсутствии загрязнений - c_o=9.8 мг/л.

 

Список литературы

1. Федеральный закон «Об охране атмосферного воздуха» №96-ФЗ. 1999г.

2. Постановление Правительства РФ №461. 2000 г. «Правила разработки и утверждения нормативов образования отходов и лимитов на их размещение».

3. Агаханянц П.Ф. Средства и методы контроля состояния окружающей природной среды: Учебное пособие/СПбГИЭУ.-СПб.:СПбГИЭУ,2011.-167 с. ПК

4. СанПиН 2.1.4.1074-01 - Гигиенические нормативы содержания вредных веществ в питьевой воде.

5. Закон «Об охране окружающей среды» М. 2002г

6. Закон «Об экологической экспертизе», М. 1995г.

7. Алексеев А.С. Мониторинг лесных экосистем. 2-изд. СПб., 2003. 116 с.

8. Методика расчета концентраций в атмосферном воздухе вредных веществ, содержащихся в выбросах предприятий. ОНД – 86. – Л.: Гидрометеоиздат, 1987. – 94 с.

9. Озерова Е.М.Оценка и прогнозирование состояния окружающей среды:Конспект лекций для специальности 080502/4 Экономика и управление на предприятии природопользования. СПб:СПбГИЭУ, 2011 – 83 с.

10. Алексеев А.С. Математические модели и методы в лесном хозяйстве. Л.: ЛТА, 1988, 88 с.

Оглавление

Введение. 4

Основные модели роста численности популяций. 4

Демографические таблицы популяций. 18

Расчет компонент сбалансированного техноценоза. 23

Зонирование территории по степени поврежденности экосистем. 27

Расчет рассеивания выбросов в атмосфере по методике ОНД-86. Определение норматива ПДВ и необходимой степени очистки выбросов. 34

Расчет платы за загрязнение атмосферы.. 41

Расчет платы за загрязнение водных объектов. 46

Организация системы мониторинга заданного объекта. 49

Расчет материального баланса веществ при сжигании твердого топлива. 56

Расчет разбавления вредных примесей по методу В.А. Фролова - И.Д. Родзиллера. Определение норматива НДС и необходимой эффективности работы очистных сооружений. 60

Расчет кислородного баланса водоема. Модель Стритера-Фелпса. 64

Список литературы.. 68

Введение

 

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

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

 

 

Занятие 1

Основные модели роста численности популяций.

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

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

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

Рассмотрим основные, наиболее часто используемые для описания роста популяций закономерности, каковыми являются гиперболическая, экспоненциальная, J - образная и логистическая (S - образная).

 

1. Гиперболический рост.

Графически изменение численности популяции по гиперболе описывается графиком, изображенным на рис. 1.1.

n

Время, t.

n0

t*

Рис. 1.1. График гиперболического роста популяции.

Динамика роста популяции по гиперболе описывается следующим уравнением:

(1.1)

 

 

его решение имеет вид:

 

 

 

 

 

 

(1.2)

 

где n - численность популяции, t* - момент времени, при котором численность популяции становится равной бесконечно большой величине, а – константа.

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

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

 

Пусть

 

Тогда исходная функция примет следующий вид:

 

Введем функцию

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

 

 

(1.3)

 

где N – число замеров численности особей в популяции. Решая полученную систему уравнений любым доступным методом, можно вычислить коэффициенты "а" и "b"

При решении системы уравнений полезно составить следующую таблицу:

Таблица 1.1.

Форма таблицы для решения системы уравнений (для примера по первому варианту, описанного в задании)

ni	ti	ni^'=1/ni	ti^2	ni^'∙ti
		5*10^-3
		3,3*10^-3	10^6	3,3
		2,9*10^-3	1,44*10^6	3,43
		1,5*10^-4	4*10^6	0,3
	∑ti=	∑ni^'=	∑(ti^2)=	∑(ni^'∙ti)=

 

2. Экспоненциальный рост.

Уравнение, описывающее изменение численности популяции имеет следующий вид:

(1.4)

 

его решение имеет следующий вид:

(1.5)

Здесь: r – биотический потенциал популяции, n₀ - начальная численность. Биотический потенциал - величина, характеризующая на какую долю от общего количества изменяется численность популяции в единицу времени.

Вид кривой, описывающий экспоненциальный рост, представлен на рис. 1.2.

n

Рис.1.2. Рост популяций по экспоненте

Время, t.

n0

Для экспоненциально растущей популяции важной характеристикой является время удвоения ее численности – T.

 

при

 

t = Т

 

(1.6)

 

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

 

(1.7)

 

откуда получаем выражение (1.8)

 

(1.8)

 

с учетом (1.8) выражение (1.5) принимает следующий вид (1.9):

 

(1.9)

 

где t/T – число удвоений численности популяции, произошедшее за время t.

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

Многие популяции, растущие по экспоненте, представляют хозяйственную ценность и могут активно эксплуатироваться. Процесс эксплуатации экспоненциально растущей популяции описывается с помощью коэффициента промысловой смертности – k, характеризующего интенсивность изъятия особей из популяции.

 

 

(1.10)

 

 

t

n

n0

Щадящий режим

r>k

t

n

n0

Стационарный режим

r=k

t

n

n0

Режим на уничтожение

r<k

Рис.1.3. Режимы эксплуатации экспоненциально растущей популяции

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

 

3. Рост с ограничением или J - образный.

Уравнение, описывающее рост популяции по J-образной кривой, имеет следующий вид

 

(1.11)

 

его решение:

 

(1.12)

 

Где: К – максимально возможная численность популяции, т.е. емкость среды обитания.

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

 

 

 

Подставив последние выражения в (1.12) можно получить следующее линейное рекуррентное соотношение между соседними значениями численности популяции, отстоящими друг от друга во времени на величину dt:

 

(1.13)

 

где и - численности популяции в два смежных, отстоящих друг от друга на величину времени dt момента времени.

Для определения констант А и В применим к полученному линейному выражению метод наименьших квадратов, описанный выше:

 

(1.14)

 

Таблица 1.2.

Форма таблицы для решения системы уравненийJ-образного роста (на примере первого варианта)

 

n1	n2	n1∙n2	n1^2
∑ n1	∑ n2	∑ (n1∙n2)	∑ n1^2

 

Решая систему (1.14) получаем значения А и В. После чего не составляет труда определить параметры J - образной кривой:

 

(1.15)

Величина dt=10 для всех вариантов.

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

 

n

t

n0

K

Рис.1.4. Рост с ограничением

 

4. Логистический или S - образный рост. В этом случае уравнение, описывающее рост популяции имеет вид:

 

(1.16)

 

(1.17)

 

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

 

 

введем новые обозначения:

 

 

тогда:

 

(1.18)

 

(1.19)

(1.20)

 

Для решения системы уравнений 1.20 полезно использовать таблицу 1.3

 

Таблица 1.3.

Форма таблицы для решения системы уравнений (1.20) логистического роста. (на примере первого варианта)

n1	n2
		∑	∑	∑

 

С помощью выражения (1.20) по данным о численности популяции в смежные равноотстоящие моменты времени могут быть определены ее биотический потенциал и емкость среды обитания. Принцип решения системы уравнений (1.20) аналогичен рассмотренному выше для гиперболического типа роста.

 

1.21

 

Величина dt=10 для всех вариантов.

График, описывающий логистический рост имеет следующий вид (рис. 1.5):

n

t

n0

I

II

III

Рис.1.5. Логистический рост популяций

K

Как видно из приводимого графика, на S образной кривой выделяют три участка:

I - начальная стадия. Отличается низкой скоростью роста численности популяции. Это обусловлено низкой плотностью популяции, вследствие чего особи противоположного пола встречаются достаточно редко.

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

III - стадия замедления. Отличается замедлением скорости роста числа особей в популяции. Это связано с тем, что ресурсов среды обитания становится недостаточно для обеспечения всех особей.

Оценка параметров логистического закона в этом случае может быть осуществлена с помощью соотношения (1.20) или другим специальным методом.

 

КОНТРОЛЬНОЕ ЗАДАНИЕ.

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

Таблица 1.4

ЧИСЛЕННОСТЬ МИРОВОГО НАРОДОНАСЕЛЕНИЯ

Год
Млн. чел

 

Год
Млн. чел

 

2. Определить массу зерна, которую необходимо положить на последнюю клетку шахматной доски, если, начиная с одного зернышка на первой клетке на каждую следующую кладется в два раза большее число зерен. Вес одного зерна принять равным 0.1 г.

 

3. Определить время, за которое потомство одной бактерии, размножающейся делением и весом 10^-12 г. достигнет веса Земного шара - 10²⁸ г., если одно деление происходит с интервалом в 20 минут.

 

4. Экспоненциально растущая популяция эксплуатируется в стационарном режиме, необходимо увеличить снимаемый с нее урожай в «к» раз за время «t». Определить, при каких «t» возможно увеличение урожая в «к» раз.

 

5. Популяция эксплуатируется в режиме на уничтожение.Определить время, за которое численность популяции сократится до 10% от первоначальной, см. табл.1.4.

 

Таблица 1.4

Исходные данные для выполнения контрольного задания по вариантам

Варианты	n0	r	k
		1.0	1.3
		0.8	1.0
		1.2	5.6
		0.3	1.2
		0.6	1.5
		1.3	2.6
		0.1	0.3
		0.7	1.6
		0.08	0.1
		0.7	1.4

 

6. Площадь лесов планеты составляет приблизительно 3.85 млрд. га. С учетом процессов искусственного и естественного возобновления скорость сведения лесов составляет 20 га/мин. Определить время, за которое площадь лесов сократится до 10% от первоначальной, предполагая, что их сведение происходит по экспоненте.

 

7. По данным таблиц хода роста древостоев основных лесообразующих пород начиная с возраста 60 лет и старше определить их популяционные параметры приняв в качестве закона роста J - образную кривую. cм.табл.1.5.

 

8. По данным таблиц хода роста основных лесообразующих пород определить их популяционные параметры приняв в качестве закона роста модифицированный логистический. Использовать данные о ходе роста по запасу. cм.табл.1.5.

 

Таблица 1.5.

ХОД РОСТА ОСНОВНЫХ ЛЕСООБРАЗУЮЩИХ ПОРОД, м³/га.

Возраст,	Номера вариантов и данные для определения популяционных параметров
лет

 

Занятие 2