ВведениеВ товарной аквакультуре скорость роста навески рыбы является одним из важнейших показателей экономической и производственной эффективности. Значения показателей ожидаемого прироста рыбы лежат в основе производственного плана и определяют планы потребления кормов, логистики кормов, сортировки, вылова и отгрузки рыбы. На скорость роста рыбы влияет большое количество параметров: качество рыбопосадочного материала, возраст рыбы, используемые корма и режимы кормления, качество и химический состав воды и ряд других факторов. При выращивании рыбы в установках замкнутого водоснабжения (УЗВ) число неуправляемых параметров значительно снижается по сравнению с садковым или прудовым выращиванием.Для анализа прироста, как правило, используются показатели абсолютной, относительной и удельной скорости роста рыбы, которые относятся к описательным моделям роста [1]. Популярность этих математических моделей связана с простотой использования, хорошей сопоставимостью результатов и понятной биологической интерпретацией. Однако по построению такие модели являются апостериорными и не могут дать прогнозных значений. К описательным моделям также относится температурный коэффициент роста (thermal growth coefficient, TGC), который может быть использован для прогнозирования размеров рыбы. Однако в исследованиях отмечается, что TGC не обладает хорошей точностью, поэтому предпринимаются попытки улучшения модели (см., например, [2]).Помимо простых описательных моделей существует целый ряд нелинейных математических моделей роста рыбы, обзоры которых можно найти, например, в работах [1, 3–5]. Наиболее известной моделью роста является уравнение фон Берталанфи, основанное на предположении о том, что изменение массы тела во времени является результатом различия между процессами анаболизма и катаболизма [6]. На основе базовой модели роста фон Берталанфи были разработаны несколько ее разновидностей, например линейная и сезонная модели Берталанфи [7, 8]. Как отмечают авторы работы [3], к наиболее часто используемым моделям в рыбохозяйственных исследованиях относятся трехпараметрическая модель фон Берталанфи, трехпараметрическая модель Гомперца и трехпараметрическая логистическая модель.В общем случае рост рыбы, в том числе и аквакультурной радужной форели, описывается нелинейным процессом, состоящим из 3-х фаз: экспоненциальной (молодь), линейной (молодая рыба) и ограниченной (зрелая рыба), поэтому соответствующая математическая модель роста должна описывать некоторую S-образную кривую. При исследовании процесса роста следует использовать мультимодельный подход, при котором на каждом этапе жизни для моделирования используется отдельная функция. Так, например, в работе [9] делается вывод о том, что для моделирования роста молоди лучше использовать степенную функцию роста, а для роста взрослой рыбы лучше подходит модель фон Берталанфи.Как правило, математические модели роста рыбы не учитывают внешних условий и режимов кормления, поскольку были разработаны в рамках исследований по рыболовству. В аквакультуре для оценки роста рыбы часто используются таблицы кормления (параметрические модели), описывающие суточные нормы кормления и коэффициенты прироста в зависимости от навески рыбы и внешних параметров, таких как температура воды и количество растворенного кислорода. Цель представленной работы заключается в оценке адекватности и подборе оптимальных параметров восьми нелинейных математических моделей роста рыбы в приложении к прогнозированию роста аквакультурной форели в УЗВ. Для выполнения оценки и сравнения были рассмотрены 8 нелинейных математических моделей роста рыбы: двух- и трехпараметрические модели фон Берталанфи, двух- и трехпараметрические модели Гомперца, двух- и трехпараметрические логистические модели, экспоненциальная и степенная модели. В качестве эталона взяты табличные (параметрические) модели кормов Skretting [10] Nutra HP и Optiline RC для выращивания форели в УЗВ. С помощью стандартных инструментов библиотеки Scipy языка Python, основанных на нелинейных методах наименьших квадратов, был выполнен подбор оптимальных параметров, приближающих каждую из математических моделей к эталонной. Сравнение результатов проводилось при помощи оценки среднеквадратичной ошибки, информационного критерия Акаике с поправкой на небольшой объем выборки, байесовского информационного критерия.В итоге представлены параметры, обеспечивающие наилучшее приближение рассматриваемых нелинейных математических моделей роста рыбы к эталонной табличной функции. Результаты оценки погрешностей показали, что трехпараметрическая модель фон Берталанфи демонстрирует лучшую точность, при этом другие модели (за исключением экспоненциальной) также показали точность, достаточную для использования на практике.Практическую значимость работы представляют выводы об адекватности использования нелинейных математических моделей для моделирования роста форели в УЗВ, а также представленные параметры моделей, обеспечивающие наилучшее приближение к эталонной табличной функции.Материалы и методыВ качестве модельных данных была использована ранее упомянутая табличная модель роста радужной форели в УЗВ, представленная компанией Skretting. Табличная модель устанавливает суточные нормы кормления и кормовые коэффициенты для рыбы различной навески при различных температурах воды, определяющие траекторию роста.Для сравнения были выбраны 8 нелинейных математических моделей роста: двух- и трехпараметрические модели фон Берталанфи (2-VBGM, 3-VBGM), двух- и трехпараметрические модели Гомперца (2-GGM, 3-GGM), двух- и трехпараметрические логистические модели (2-LGM, 3-LGM), экспоненциальная модель (EGM) и степенная модель (PGM). Формулы расчета для моделей и набор параметров, некоторые из которых несут биологическую смысловую нагрузку, представлены в табл. 1, в которой используются следующие обозначения: W(t) – прогнозируемая навеска рыб (вес) в момент времени t; W∞ – предельная навеска рыбы; k – экспоненциальная скорость роста; t0 – возраст рыбы при нулевом весе; T – расположение точки перегиба; b, A0 – некоторые параметры. Таблица 1Table 1Формулы расчета для моделей ростаCalculation formulas for growth models МодельФормулаНабор параметровИсточник2-VBGM  W(t) = W∞ (1 – e–kt)b{W∞, k, b}[11]3-VBGMW(t)= W ∞ ⋅ 1- e -k(t- t 0 ) b {W∞, k, t0, b}[12]2-GGMW(t)= W ∞ ⋅ e -be -kt P= W ∞ ,k, b [13]3-GGMW(t)= W ∞ ⋅ e -be -k(t-T) P= W ∞ ,k, b, T [14]2-LGMW(t)= W ∞ /(1+ be -kt ) P= W ∞ ,k, b [13]3-LGMW t = W ∞ /( 1+be -k(t-T) ) P= W ∞ ,k, b, T [15]EGMW(t)= A 0 e kt P= A 0 ,k [16]PGMW(t)= A 0 t k P= A 0 ,k [17]  Для оценки параметров каждой модели использовался алгоритм Левенберга – Марквардта, который является одним из нелинейных методов наименьших квадратов [18, 19], реализованныйс помощью стандартных инструментов библиотеки Scipy языка Python. Точность подобранных моделей роста была оценена с помощью следующих метрик эффективности:– среднеквадратичная ошибка где n – объем данных; yi – исходные данные;  – прогнозируемые данные; – информационный критерий Акаике (AICc) с поправкой на небольшой объем выборки ,где K – количество параметров в модели [20, 21]; – байесовский информационный критерий [22] BIC = K ln(n) + n ln(MSE).  Модель с наименьшими значениями среднеквадратической ошибки и информационных критериев считается наилучшей для моделирования роста радужной форели. Результаты и обсуждениеРезультаты подбора параметров моделей, обеспечивающих наилучшее приближение к эталонной табличной функции, представлены в табл. 2.  Таблица 2Table 2Оценка набора параметров восьми моделей ростаEvaluation of a set of parameters for eight growth modelsМодельПараметрW∞kbt0TA02-VBGM2,6938E+127,04E-072,60–––3-VBGM205 530,560,000593,32–42,54––2-GGM8 222,110,00377,25–––3-GGM8 222,110,00376,06–48,54–2-LGM2 350,490,0138131,24–––3-LGM2 350,490,01389,85–187,86–EGM–0,00935–––42,2PGM–2,5985–––0,00028  Отдельные параметры исходных математических моделей имеют биологическое наполнение. Так, параметр W∞ определяет предельную навеску рыбы, однако в результате подбора оптимальное значение варьирует в диапазоне от 2 350,49 г до 2,6938E+12 г, где максимальные значения получены для моделей фон Берталанфи, а наименьшие – для логистических моделей. Учитывая, что вес наиболее типичной особи радужной форели составляет 2 000–3 000 г, а максимальный зафиксированный вес форели достигает 25 400 г [23], наиболее адекватные значения параметра W∞ дают логистические модели и модели Гомперца. Параметр экспоненциальной скорости k находится в диапазоне от 2,49E-07 до 2,5985, причем его наименьшие значения получены для двухпараметрической модели фон Берталанфи. Максимальные значения параметра k имеет степенная модель. Параметр t0, указывающий на возраст молоди с нулевым весом в трехпараметрической модели фон Берталанфи, имеет отрицательное значение. Такое значение параметра объясняется тем, что в таблице кормления отсчет возраста рыб с массой 2,1 г начинается с 1. Потеря биологического содержания подобранных параметров является следствием ограниченности практической максимальной навески форели, выращиваемой в УЗВ.О степени близости нелинейных моделей к эталонной можно судить по показателям эффективности, представленным в табл. 3.Можно заметить, что модели показывают невысокую погрешность, а значит, каждую из них (кроме экспоненциальной модели) можно использовать на практике.В последнем столбце табл. 3 представлен рейтинг моделей относительно их приближенности к эталонной. Согласно рейтингу, трехпараметрическая модель фон Берталанфи демонстрирует наилучшую точность аппроксимации.  Таблица 3Table 3Оценки качества аппроксимации моделейEstimates of the quality of model approximation МодельКритерийРейтингMSEAICcBIC2-VBGM  19,441 160,352 631,7543-VBGM7,61797,53813,2812-GGM36,571 406,161 417,9933-GGM36,571 408,211 423,2142-LGM277,462 194,452 206,2853-LGM277,462 196,492 212,246EGM2 061,362 972,542 980,437PGM19,421 157,861 165,752  ЗаключениеВ товарной аквакультуре скорость роста навески рыбы лежит в основе производственного плана и определяет планы потребления кормов, логистики кормов, сортировки, вылова и отгрузки рыбы. Для контроля скорости роста рыбы на предприятиях аквакультуры в основном применяются табличные (параметрические) модели, разрабатываемые производителями кормов для своей продукции. Наряду с табличными моделями существуют различные математические модели, используемые для прогнозирования и оценки темпов роста рыбы. Среди них наиболее часто используются нелинейные модели, такие как двух- и трехпараметрические модели фон Берталанфи, двух- и трехпараметрические модели Гомперца, двух- и трехпараметрические логистические модели, экспоненциальная и степенная модели.В представленной работе выполнена оценка адекватности и подбор параметров восьми нелинейных математических моделей в приложении к прогнозированию роста аквакультурной форели в УЗВ в сравнении с принятой за эталонную табличной моделью. Такой подход имеет некоторые особенности. Во-первых, навеска рыбы по практическим причинам ограничивается 1,5 кг и двумя стадиями роста (молодь и молодая рыба), не образующими S-образной кривой. Во-вторых, в силу постоянства внешних условий в УЗВ, полученные модели имеют прогностическое свойство.В ходе исследования получены параметры, наилучшим образом приближающие рассматриваемые модели к эталонной. Результат оценки погрешностей показал, что трехпараметрическая модель фон Берталанфи демонстрирует лучшую точность. При этом другие модели (за исключением экспоненциальной) также показали точность, достаточную для использования на практике. Необходимо отметить, что в силу ограниченной навески рыбы параметры моделей перестают нести исходное заложенное биологическое наполнение (в частности, W∞ – предельная навеска рыбы; t0 – возраст рыбы при нулевом весе).Результаты исследования могут использоваться для моделирования роста форели на действующих предприятиях.