Дискретно-непрерывная модель в задаче анализа критического уровня эксплуатации биоресурсов

Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2014. № 1 (34). С. 145—155

УДК 517.958:57; 517.938

ДИСКРЕТНО-НЕПРЕРЫВНАЯ МОДЕЛЬ В ЗАДАЧЕ АНАЛИЗА КРИТИЧЕСКОГО УРОВНЯ ЭКСПЛУАТАЦИИ БИОРЕСУРСОВ

А. Ю. Переварюха

Санкт-Петербургский институт информатики и автоматизации РАН,

Россия, 199178, Санкт-Петербург, 14 линия, 39.

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

Ключевые слова: динамика популяционных процессов, моделирование эксплуатации биоресурсов, сценарные модели.

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

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

ISSN: 2310-7081 (online), 1991-8615 (print); doi: http://dx.doi.org/10.14498/vsgtu1116 © 2014 Самарский государственный технический университет.

Образец цитирования: А. Ю. Переварюха, “Дискретно-непрерывная модель в задаче анализа критического уровня эксплуатации биоресурсов” // Вестн. Сам. гос. техн. унта. Сер. Физ.-мат. науки, 2014. № 1(34). С. 145-155. doi: 10.14498/vsgtu1116.

Сведения об авторе: Андрей Юрьевич Переварюха (к.т.н., с.н.с.), старший научный сотрудник, лаб. прикладной информатики.

E-mail address: madelf@pisem.net

145

А. Ю. Переварюха

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

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

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

У. Рикер [2] предложил часто применяемую функцию R = f (S) связи «запасы-пополнение», отражающую негативное действие на выживаемость поколения повышенной плотности запаса. Будем рассматривать эту модель как дискретную динамическую систему в виде итераций {Aj)}jyo, где Ro, Ri, R2, ... — последовательность точек, описывающих эволюцию системы и определённых условием

Rj+1 = aRj exp(- bRj).

Её свойства как функции и недостатки как популяционной модели обсуждались неоднократно, например в [3].

Биологический смысл модели становится яснее, если перейти к дифференциальному уравнению убыли текущей численности N(t):

dN

_ = —

Константы заданного на промежутке времени обыкновенного дифференциального уравнения (ОДУ) соотносятся с константами формулы Рикера:

a = Л exp(—в0), b = a0,

где Л — средняя плодовитость популяции. Согласно теории, повышенная плотность на нерестилищах становится негативным фактором, увеличивающим смертность. Предполагается, что на смертность влияет именно начальная

146

Дискретно-непрерывная модель в задаче анализа критического уровня эксплуатации биоресурсов

численность поколения. Большое значение имеет длина интервала $ — периода уязвимости поколения по отношению к основным факторам смертности. Мы предлагаем на практике для вариативности в оценках параметров использовать модель именно в виде численного решения ОДУ (1) с начальным условием N(0) = ЛS.

Динамическая система V j > 0 обладает аттрактором A таким, что

VR0 lim {^>(j)} = A.

Для простого случая аттрактор глобальный и является устойчивым состоянием равновесия:

A = R* = ф(Я*).

Первый метаморфоз A определён условием нарушения критерия устойчивости неподвижной точки

\ф'Д*)\ < 1,

следующего из теоремы, независимо доказанной Д. Гробманом и Ф. Хартманом в [4,5]. При значении производной ф1(Е¥) = —1, когда a = e2, происходит бифуркация удвоения периода. У функциональной итерации возникают две новые циклические точки

yn(R*) = yn+2(R*),

являющиеся неподвижными точками второй итерации ^2(R). При дальнейшем увеличении параметра две неподвижные точки ф2(Д) аналогично потеряют устойчивость и появятся n = 4 точек ^4(R). Любая ранее устойчивая R* становится репеллером и делит область притяжения образовавшегося цикла на субинтервалы Yi, Т2. Для R0 « R0, R0 е Yi, R0 е Т2, отличается порядок обхода цикла

lim ^(R0)n = R\ при n = 2г

П^Ж

и, соответственно,

lim ф(К0)п = R\ при n = 2г.

П^Ж

При изменении управляющего параметра в диапазоне значений e2 < a < a реализуется каскад бифуркаций удвоения периода, происходящих при ф2п '(R*) = —1.

Когда период цикла становится бесконечным, образуется странный аттрактор, гомеоморфный канторовскому множеству. Скорость увеличения периода цикла описывается универсальными константами М. Фейгенбаума [6], в честь которого называют данный сценарий хаотизации. Из-за бесконечного количества малых субинтервалов наблюдается чувствительная зависимость от начальных условий. Практические цели предметной области диктуют необходимость биологической интерпретации сценария хаотизации. Бифуркационный параметр a в модели Рикера трактуется как мера репродуктивного потенциала популяции, что приведёт исследователя к очевидным заключениям. Мы отметим ряд проблем для поиска практических решений при моделировании. Во-первых, сценарий хаотизации связан с возникновением

147

А. Ю. Переварюха

ряда других нелинейных эффектов. Метаморфозы происходят как с аттрактором, так и с границами его области притяжения, и подобным явлениям трудно дать популяционную интерпретацию. Во-вторых, существуют альтернативные модели вида «запас - пополнение», в которых хаотизация происходит при увеличении параметров, трактуемых иначе. В-третьих, наблюдается огромная, маловероятная для реальных популяционных процессов, амплитуда хаотических колебаний. Сценарий хаотизации является свойством определенного класса дискретных динамических систем SU-отображений, значит, необходимо вести речь о биологической интерпретации свойств отображений, удовлетворяющих критериям теоремы Д. Сингера [7]. Задачи промысловой ихтиологии пересекаются с фундаментальными и ещё не полностью разрешёнными проблемами нелинейных эффектов, т. к. в реальных моделях может присутствовать несколько бифуркационных параметров.

2. Модель регулируемого плотностью воспроизводства. Характер зависимости решения (1) между запасом и пополнением с резко выраженным единственным экстремумом далеко не всегда согласуется с данными наблюдений ихтиологов и с точки зрения математики означает возможность появления циклов всевозможных периодов в соответствии с порядком теоремы А. Н. Шарковского [8,9]. Опубликованные результаты наблюдений [10] свидетельствуют о возможности появления двух максимумов Smi, Sm2 на графике зависимости нерестового запаса и пополнения. Очевидный и простой способ описания сложной формы кривой — воспользоваться второй итерацией формулы Рикера

f (f (S)) = f 2(S).

Подобные попытки моделирования известны. Однако трудно представить, что выполнится

f (Smi) = f (Sm2),

более того, данный способ некорректен из-за свойств функции R = f (S):

df 2(S) dS

df (S) dS

ae bS(1 — bS),

lim

S^Q

df (S) dS

a;

a2 exp (—bS — abSe bS)(1 — abSe bS — bS + ab2S2e bS),

lim

S^Q

df 2(S) dS

2

a

5

подразумевающих, что при деградации запаса скорость воспроизводства стре-

22 мится к предельному значению a, но в f2 возрастает уже до a2.

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

Rn+i = f (f (Rn)), f (R) = aRe~bR.

После потери устойчивости возникает не цикл, а аттракторы Ai, A2, с которыми происходят удвоения периода. Задачу описания разнообразных сложных форм зависимости мы предлагаем решать при помощи введения функционалов с подобранной ограниченной областью значений E(0) = [1,1), I > 0 в правую часть (1).

148

Дискретно-непрерывная модель в задаче анализа критического уровня эксплуатации биоресурсов

Анализ опубликованных данных о воспроизводстве находящейся на грани выживании популяции волжского осетра позволил сделать вывод о сильном действии эффекта Олли [11]. Это описанное в 30-годы XX века зоологами нелинейное явление выражается в том, что если плотность популяции в большом ареале обитания становится ниже оптимальной, то непропорционально резко падает эффективность воспроизводства. Оценка смещения максимума распределения веса молоди в прудах при увеличенной сверх норматива плотности позволил нам предложить модель в виде системы ОДУ с функционалом 0(S), описывающую убыль поколения на интервале уязвимости t £ [0, T]:

dN

dt

dw

dt

— (aw(t)N (t) + 0(S)e)N (t), g ---------- l

Nk + (’

0(S ) =

1 — exp(—cS) ’

где S — величина нерестового запаса; w(t) — функция, отражающая условный уровень размерного развития поколения, влияющий на увеличение пищевых потребностей; g — параметр, учитывающий ограниченность количества доступных кормовых объектов; 0(S) — некоторая убывающая функция (0(S) ^ 1 при S ^ те), не влияющая на вычисление N(T); T = 45 сут., если численность запаса достаточно велика.

Введение в систему ОДУ (2) быстро убывающей функции 0(S) отражает снижение эффективности воспроизводства при деградации популяции, связанное с уменьшением вероятности встречи особей в местах размножения. Влияние данного эффекта можно оценить по имеющимся данным о состоянии осетровых [12]. Далее в (2): ( — параметр, учитывающий ограничение темпов развития, не зависящее от численности популяции; c — параметр, характеризующий степень выраженности эффекта Олли; а — мгновенный коэффициент компенсационной смертности; в — мгновенный коэффициент декомпенсационной смертности. Начальные условия для системы (2) определяются следующим образом1:

N

Типичный график моделируемой зависимости запаса и пополнения R = = f (S) = N(T) (кривая 1 на рис. 1), полученный при численном решении задачи Коши (2), (3) в инструментальной среде AnyLogic, является унимодальной кривой с двумя нетривиальными пересечения R|, R| с биссектрисой N(T) = S (прямая 2 на рис. 1). Кривая 1 отличается пологой ниспадающей правой ветвью с ненулевой горизонтальной асимптотой и крутым наклоном восходящей ветви только после первого пересечения с биссектрисой координатного угла R| < R < Rmax.

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

Ni+1(0) = Nj(T )A, w(0) = const,

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

149

А. Ю. Пер ев ар юх а

Рис. 1. Типичный график R = f

[Figure 1. Typical graph of R = f (S) (curve number 1) obtained from the numerical solution of the Cauchy problem (2), (3) with t = T, and image of coordinate angle bisector (line number 2)]

и задающих дискретную траекторию. При организации модельных экспериментов необходимо вводить соответствие условного модельного времени отрезков t £ {[0,T]i} календарному. Использование событийно-управляемого подхода позволяет выделить некое особое, приводящие к переменам состояние, что в сочетании с непрерывным описанием популяционного процесса перспективно при рассмотрении последствий изменяющегося по некоторым условиям внешнего воздействия.

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

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

Поведение дискретно-непрерывной системы с оператором эволюции на основе численного решения (2) качественно отличается от итераций функций Рикера и Шепарда разделением фазового пространства на две области притяжения. Граница областей dQi представляет простой случай для мультистабильных систем, являясь неустойчивой точкой равновесия дQi = R|, или «репеллером». Анализ устойчивости неподвижных точек динамической системы, реализованной в инструментальной среде моделирования, можно проводить с использованием свойства второй итерации ф2(х). Необходимым

150

Дискретно-непрерывная модель в задаче анализа критического уровня эксплуатации биоресурсов

и достаточным условием устойчивости неподвижной точки х* одномерного отображения являются неравенства

^2(х) > х при х < х* и ^2(х) < х при х > х*.

Устойчивым является тривиальное состояние равновесия, таким образом, модель подразумевает наличие критически низкой численности для существования биологического вида. Интересным представляется поведение траектории в окрестности точки репеллера R*±e, когда под внешним воздействием, характеризующим некоторой неопределённостью, траектория медленно покидает эту окрестность.

Рассмотрим сценарий, когда после длительного периода получения стабильных уловов возникает предположение, что запас популяции Sst недоиспользуются промыслом. Пример подобных заключений экспертов о промысле осетровых в Волге — статья [13]. Отметим, что учёт запасов специалистами производится по количеству пропущенных на нерест особей, тогда как статистика промысла ведётся по массе уловов, что может вносить искажения в оценки имеющихся данных. Разработанная модель пополнения оперирует единицами численности. Пусть первоначально установленная доля изъятия Но в момент 1 увеличивается:

Н1 = Н0 + АН, АН и 0.15.

В дискретно-непрерывной системе начальные условия теперь сопрягаются так:

NmLo = (1 — Н1)^и.

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

В сцепленных сценариях будем моделировать популяции с различной введённой дополнительной промысловой нагрузкой и последующим её изменением в некоторый момент времени, получая на графике две временные диаграммы величины уловов. На диаграммы на ось t, где отложена шкала времени инструментальной среды, так же выводится пересчёт модельного времени в нумерованные календарные сезоны Y (одному сезону соответствуют 8 тактов модельного времени AnyLogic5). Для удобства наблюдения отсчёт запускается от момента 1, когда был изменен режим промысла. Далее принятие решений и развитие процесса может происходить разными путями, отражёнными в предикативно заданных условиях двух модельных экспериментов с одинаковой исходной численностью запасов при Но = 0.4, Л = 227000.2

После кратковременного увеличения уловы снижаются до прежних значений и продолжают медленно уменьшаться со скоростью, зависящей от АН. Эксперты предполагают стабилизацию популяции, не вносят изменения в режим промысла и далее происходит резкое падение уловов, после которого остановка лова не спасает от деградации (кривая 1 на рис. 2):

Ri < Sst.

2Взята оценка среднемноголетней плодовитости каспийской севрюги.

151

А. Ю. Пер ев ар юх а

Время, календарные сезоны [Time (Calendar Seasons)]

Такты модельного времени [Simulation Time Steps]

Рис. 2. Сценарии деградации при различном уровне перелова [Figure 2. Degradation scenarios at various level of overfishing]

Длительный сценарий деградации растянулся на 33 модельных сезона. Предположим, для второй популяции Hi значительно превышает оптимальный уровень эксплуатации. Объёмы вылова падают быстрее, и эксперты, осмыслив негативные последствия, решают через 21 сезон вернуться к уровню Но, когда состояние популяции находится в е окрестности репеллера R1 (псевдостабильно). Некоторое время уловы постоянны, но незначительные изменения смертности неожиданно резко выводят популяцию из нестабильного равновесия, и происходит деградация (кривая 2 на рис. 2). Поведение популяции низкой численности перестаёт удовлетворительно описываться детерминированной динамикой средних величин, и мы получаем очень мало сведений о реальном состоянии потерявших промысловое значение биоресурсов.

В случае, когда эксперты после получения данных о сокращении уловов за предыдущие 3 сезона своевременно переводят промысловый режим Н2 = = Но — АН при Rj > Rl, обе имитированные в сценариях популяции за 39 модельных сезонов восстанавливают свои промысловые запасы (рис. 3).

Длительный сценарий деградации (подобный кривой 1 на рис. 2) реализо-

Время, календарные сезоны [Time (Calendar Seasons)]

О 6 12 18 24 30 36 42

Такты модельного времени [Simulation Time Steps]

152

Рис. 3. Сценарии восстановления после прекращения перелова [Figure 3. Restore scenarios because of a stop overfishing]

Дискретно-непрерывная модель в задаче анализа критического уровня эксплуатации биоресурсов

вался для популяций осетровых Каспийского моря. В 1979 г. уловы достигли рекордных 27 тыс. тонн. В 1988-89 гг. произошло резкое падение уловов, промысел не был приостановлен вовремя, эксперты рассчитывали на успех технологии искусственного воспроизводства. Только в ноябре 2010 г. по инициативе Казахстана главы пяти государств Каспийского бассейна одобрили введение пятилетнего моратория на промышленный лов осетровых. Три вида включены в список Конвенции по международной торговле видами дикой фауны и флоры, находящимися под угрозой уничтожения (CITES). Однако остановить незаконный промысел, несмотря на декларируемые усилия, не удаётся, т. к. браконьерство является частью системной социальной проблемы.

За счёт гибкой системы регулирования удалось предотвратить деградацию запасов горбуши тихоокеанского побережья Канады, где установился баланс «промысел - пополнение». Падение уловов привело к уменьшению (из-за низкой рентабельности промысла) количества купивших лицензию на лов рыболовецких судов.

Заключение. Данные о динамике в системе «запас - пополнение - промысел» для популяций рыб почти всегда неполны, приблизительны и труднодоступны для стороннего исследователя. В отечественной практике редко можно говорить о наличии достоверных данных об изменении одновременно и запаса, и промыслового усилия. В ситуации с осетровыми Каспийского моря данные о промысловом усилии неизвестны. Сведения о современных размерах нелегального промысла крайне противоречивы. В работе [14] были опубликованы данные о динамике вылова, количестве орудий лова и состоянии запасов основных промысловых популяций Горьковского водохранилища. Сведения подтверждают динамику развития «оптимистичного» сценария. После начавшегося перелова численность популяции резко снижается и далее «псевдостабилизируется». В случае снижения промыслового усилия происходит быстрый рост численности видов рыб, т. к. эффективность воспроизводства из-за снижения конкуренции на ограниченных пригодных нерестилищах в водохранилищах достаточна. Промысел может использовать столь продуктивное состояние популяции, однако оно близко к Щ, и из-за кажущего незначительным превышения вылова может вновь произойти истощение ресурсов популяции на длительный период.

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

Работа выполнена при поддержке РФФИ (проект № 14-07-00066).

This work is supported by RFBR, project no. 14-07-00066.

153

А. Ю. Переварюха

СПИСОК ЛИТЕРАТУРЫ/ REFERENCES

1. R. M. May, “Simple mathematical models with very complicated dynamics”, Nature, 1976, vol. 261, pp. 459-467. doi: 10.1038/261459a0.

2. W. E. Ricker, “Stock and Recruitment”, J. Fish. Res. Bd. Canada, 1954, vol. 11, no. 5, pp. 559-623. doi: 10.1139/f54-039.

3. А. Ю. Переварюха, “Анализ опыта моделирования биоресурсов с точки зрения современной теории динамических систем”// Экологические системы и приборы, 2011. №11. С. 17-21. [A. Yu. Perevarykha, “Analysis of the simulation practices of bioresources in terms of modern theory of dynamical systems”, Ekologicheskiye sistemy i pribory [Ecological Systems and Devices], 2011, no. 11, pp. 17-21. (In Russian)].

4. Д. М. Гробман, “Гомеоморфизм систем дифференциальных уравнений” // ДАН СССР, 1959. Т. 128, № 5. С. 880-881. [D. M. Grobman, “Homeomorphisms of systems of differential equations”, Dokl. Akad. Nauk SSSR, 1959, no. 5, pp. 880-881. (In Russian)].

5. Ph. Hartman, “On the local linearization of differential equations”, Proc. Amer. Math. Soc., 1963, vol. 14, no. 4, pp. 568-573. doi: 10.1090/s0002-9939-1963-0152718-3.

6. M. J. Feigenbaum, “Universal behavior in nonlinear systems”, Physica D: Nonlinear-Phenomena, 1983, vol. 7, no. 1-3, pp. 16—39. doi: 10.1016/0167-2789(83)90112-4.

7. D. Singer, “Stable Orbits and Bifurcation of Maps of the Interval”, SIAM J. Appl. Math., 1978, vol. 35, no. 2, pp. 206-267. doi: 10.1137/0135020.

8. А. Н. Шарковский, “Сосуществование циклов непрерывного преобразования прямой в себя”// Укр. мат. ж., 1964. Т. 16, №1. С. 61-71. [A. N. Sharkovskii, “Coexistence of cycles of a continuous transformation of a line into itself”, Ukrain. Mat. Zh., 1964, vol. 16, no. 1, pp. 61-71. (In Russian)].

9. А. Н. Шарковский, С. Ф. Коляда, А. Г. Сивак, В. В. Федоренко, Динамика одномерных отображений, Киев: Наукова думка, 1989. 216 с.; A. N. Sharkovsky, S. F. Kolyada, A. G. Sivak, V. V. Fedorenko, Dynamics of One-Dimensional Maps, Mathematics and Its Applications, vol. 407, Berlin, Springer, 1997, ix+262 pp.

10. M. P. Shepard, F. C. Withler, “Spawning Stock Size and Resultant Production for Skeena Sockeye”, J. Fish. Res. Bd. Canada, 1958, vol. 15, no. 5, pp. 1007—1025. doi: 10.1139/ f58-055.

11. W. C. Allee, Edith S. Bowen, “Studies in animal aggregations: Mass protection against colloidal silver among goldfishes”, Journal of Experimental Zoology, 1932, vol. 61, no. 2, pp. 185-207. doi: 10.1002/jez.1400610202.

12. П. В. Вещев, Г. И. Гутенева, “Эффективность естественного воспроизводства осетровых в низовьях Волги в современных условиях”// Экология, 2012. №2. С. 123-128; P. V. Veshchev, G. I. Guteneva, R. S. Mukhanova, “Efficiency of natural reproduction of sturgeons in the lower Volga under current conditions”, Russian Journal of Ecology, 43, vol. 2012, no. 2, pp. 142-147. doi: 10.1134/S1067413612020154.

13. А. П. Сливка, А. В. Павлов, “Биологические основы изменения режима промысла осетровых (Acipenseridae) в дельте Волги”// Вопросы ихтиологии, 1982. №5. С. 738745. [A. P. Slivka, A. V. Pavlov, “Biological basics of regime change sturgeon fishery (Acipenseridae) in the Volga delta”, Voprosy ikhtiologii [Journal of Ichthyology], 1982, no. 5, pp. 738-745. (In Russian)].

14. Д. И. Постнов, А. Е. Минин, А. А. Клевакин, “Динамика рыбных запасов и возможности их освоения на Горьковском и Чебоксарском водохранилище” // Рыбное хозяйство, 2012. №1. С. 60-62. [A. E. Minin, A. A. Klevakin, D. I. Postnov, “Fish stocks dynamics and possibilities of their exploitation in Gorky and Cheboksary Reservoirs”, Rybnoye khozyaystvo, 2012, no. 1, pp. 60-62. (In Russian)].

Поступила в редакцию 18/IX/2012; в окончательном варианте — 12/XI/2013; принята в печать — 17/I/2014.

154

Discrete-continuous Model for the Problem of Analysis Critical Level of Exploitation of Bioresources

MSC: 37D45, 34A38; 92D25, 93A30

DISCRETE-CONTINUOUS MODEL FOR THE PROBLEM OF ANALYSIS CRITICAL LEVEL OF EXPLOITATION OF BIORESOURCES

A. Yu. Perevaryukha

St. Petersburg Institute for Informatics and Automation of RAS,

39, 14-ya liniya V. O., St. Petersburg, 199178, Russian Federation.

The article deals with the model developed in the framework of the research on the problem of replenishment of fish populations. A new approach takes into account the peculiarities of describing the changes in the number of generations by the system of ordinary differential equations. The method of models implementation in computing environment enables the simulation of various types of human impact, including the case of introduction of the new species to the environment. The author analyzes the behavior of scenarios of a population model with different levels of fishing impact and concludes the most dangerous case for the systematic management of bioresources. Examples and compared data on the dynamics of some commercial fish populations of the Volga basin are discussed. It is noted that a dangerous scenario for managing the bioresources is realized with continued slight excess of quotas.

Keywords: dynamics of population processes, modeling of exploitation of bioresources, use case model.

Received 18/IX/2012;

received in revised form 12/XI/2013;

accepted 17/I/2014.

ISSN: 2310-7081 (online), 1991-8615 (print); doi: http://dx.doi.org/10.14498/vsgtu1116 © 2014 Samara State Technical University.

Citation: A. Yu. Perevaryukha, “Discrete-continuous Model for the Problem of Analysis Critical Level of Exploitation of Bioresources”, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. & Math. Sci.], 2014, no. 1

Author Details: Andrey Yu. Perevaryukha (Cand. Technical (Computer) Sci.), Senior Researcher, Lab. of Applied Informatics.

E-mail address: madelf@pisem.net