Научная статья на тему 'ОПТИМИЗАЦИЯ ТЕХНОЛОГИЧЕСКИХ ПАРАМЕТРОВ ПРОЦЕССА СУХОГО ОХМЕЛЕНИЯ С РАЗРАБОТКОЙ МЕТОДИКИ ИНТЕГРАЛЬНОЙ ОЦЕНКИ КАЧЕСТВА ПИВА'

ОПТИМИЗАЦИЯ ТЕХНОЛОГИЧЕСКИХ ПАРАМЕТРОВ ПРОЦЕССА СУХОГО ОХМЕЛЕНИЯ С РАЗРАБОТКОЙ МЕТОДИКИ ИНТЕГРАЛЬНОЙ ОЦЕНКИ КАЧЕСТВА ПИВА Текст научной статьи по специальности «Прочие технологии»

CC BY
35
3
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ХМЕЛЬ / НЕЙРОННЫЕ СЕТИ / SAPHIRE / CHINOOK / КРИТЕРИЙ ОПТИМИЗАЦИИ / СУХОЕ ОХМЕЛЕНИЕ

Аннотация научной статьи по прочим технологиям, автор научной работы — Новикова Инна Владимировна, Рукавицын Павел Владимирович, Муравьев Александр Сергеевич

Развитие крафтовых сортов пива со сложными фруктовыми, цветочными и цитрусовыми тонами, поиск оригинального вкуса и возрождение старинных рецептур служат предпосылками для применения сухого охмеления предприятиями пивоваренной промышленности. Процесс сухого охмеления проводили динамическим способом на экспериментальной установке, включающей аппарат для охмеления, сетчатый фильтр, емкость брожения, насос. Использовали хмели Saphir и Chinook. Исследование включало определение физико-химических показатели пива: концентрация гумулинонов, изо-а-кислот и а-кислот; общее содержание полифенолов; pH; объемная плотность напитка. Показатели сенсорной оценки включали: интенсивность и качество аромата хмеля; полнота вкуса; качество горечи; общее предпочтение. Для аппроксимации результатов исследования использовали радиальные нейронные сети. Определены оптимальные технологические параметры процесса: хмель Saphire, длительность процесса 12 ч, температура 3,0 °C. При сухом охмелении хмелем Chinook оптимальная длительность процесса составляет 10,5 ч, оптимальная температура составляет 3,0 °C. При температурах процесса сухого охмеления ниже 3,0 °C критерий оптимизации получается незначительно ниже, чем для 3,0 °C. При температурах 1,5 и 2,0 °C оптимальная длительность процесса составляет 8,15 и 11,25 суток соответственно. Использование комбинации хмелей Saphire + Chinook не дает улучшение общего критерия оптимальности по сравнению с использованием только хмеля Saphire. При длительности процесса 4 ч удается достигнуть значений критерия оптимизации более 0,60, однако только при длительности порядка 12 ч удается достигнуть критерия оптимизации более 0,65.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по прочим технологиям , автор научной работы — Новикова Инна Владимировна, Рукавицын Павел Владимирович, Муравьев Александр Сергеевич

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

OPTIMIZATION OF TECHNOLOGICAL PARAMETERS OF THE DRY HOPPING PROCESS WITH THE DEVELOPMENT OF A METHODOLOGY FOR THE INTEGRAL ASSESSMENT OF THE QUALITY OF BEER

The development of craft beers with complex fruit, floral and citrus tones, the search for an original taste and the revival of old recipes are the prerequisites for the use of dry-hoping in the brewing industry. The dry-hopping process was carried out dynamically in an experimental setup with a hop apparatus, screen filter, fermentation tank and pump. Saphir and Chinook hops were used. The study included determination of physical and chemical parameters of beer: concentration of humulinones, iso-a-acids and a-acids; total content of polyphenols; pH; volume density. Sensory evaluation indicators included: hop aroma intensity and quality; flavor fullness; bitterness quality; and overall preference. Radial neural networks were used to approximate the results of the study. Optimal process parameters were determined: Saphir hops, process duration 12 hours, temperature 3,0 °C. The optimum process duration is 10,5 h and the optimum temperature is 3,0 °C. Processing temperatures below 3.0 °C result in a slightly lower optimization criterion than the 3.0 °C processing temperature. At temperatures of 1.5 and 2.0 °C the optimum process duration is 8.15 and 11.25 days respectively. The use of the combination of Saphir & Chinook hops offers no improvement in the overall optimality criterion compared with the use of Saphir alone. The process duration of 4 hours achieves an optimization criterion of more than 0.60, but only with durations of around 12 hours does an optimization criterion of more than 0.65 become attainable.

Текст научной работы на тему «ОПТИМИЗАЦИЯ ТЕХНОЛОГИЧЕСКИХ ПАРАМЕТРОВ ПРОЦЕССА СУХОГО ОХМЕЛЕНИЯ С РАЗРАБОТКОЙ МЕТОДИКИ ИНТЕГРАЛЬНОЙ ОЦЕНКИ КАЧЕСТВА ПИВА»

УДК 633.791:663.4 https://doi.org/10.36107/spfp.2021.234

Оптимизация технологических параметров процесса сухого охмеления с разработкой методики интегральной оценки качества пива

Новикова Инна Владимировна

Воронежский государственный университет инженерных технологий Адрес: 394036, Россия, г. Воронеж, пр-т Революции, д. 19

E-mail: noviv@list.ru

Рукавицын Павел Владимирович

Воронежский государственный университет инженерных технологий Адрес: 394036, Россия, г. Воронеж, пр-т Революции, д. 19

E-mail: kafedra-tbisp@mail.ru

Муравьев Александр Сергеевич

Воронежский государственный университет инженерных технологий, Адрес: 394036, Россия, г. Воронеж, пр-т Революции, д. 19

E-mail: hntrun@mail.ru

Развитие крафтовых сортов пива со сложными фруктовыми, цветочными и цитрусовыми тонами, поиск оригинального вкуса и возрождение старинных рецептур служат предпосылками для применения сухого охмеления предприятиями пивоваренной промышленности. Процесс сухого охмеления проводили динамическим способом на экспериментальной установке, включающей аппарат для охмеления, сетчатый фильтр, емкость брожения, насос. Использовали хмели Saphir и Chinook. Исследование включало определение физико-химических показатели пива: концентрация гумулинонов, изо-а-кислот и а-кислот; общее содержание полифенолов; pH; объемная плотность напитка. Показатели сенсорной оценки включали: интенсивность и качество аромата хмеля; полнота вкуса; качество горечи; общее предпочтение. Для аппроксимации результатов исследования использовали радиальные нейронные сети. Определены оптимальные технологические параметры процесса: хмель Saphire, длительность процесса 12 ч, температура 3,0 °С. При сухом охмелении хмелем Chinook оптимальная длительность процесса составляет 10,5 ч, оптимальная температура составляет 3,0 °С. При температурах процесса сухого охмеления ниже 3,0 °С критерий оптимизации получается незначительно ниже, чем для 3,0 °С. При температурах 1,5 и 2,0 °С оптимальная длительность процесса составляет 8,15 и 11,25 суток соответственно. Использование комбинации хмелей Saphire + Chinook не дает улучшение общего критерия оптимальности по сравнению с использованием только хмеля Saphire. При длительности процесса 4 ч удается достигнуть значений критерия оптимизации более 0,60, однако только при длительности порядка 12 ч удается достигнуть критерия оптимизации более 0,65.

Ключевые слова: хмель, нейронные сети, Saphire, Chinook, критерий оптимизации, сухое охмеление

Введение

Для распространенных сортов пива сейчас характерно высокое содержание ароматических компонентов, экстрагируемых из хмеля, наряду с нормируемым количеством горьких веществ (Leitao et al., 2011; Takoi et al., 2016; Walsh et al., 2018; Меледина, 2003). Развитие крафтовых сортов пива со сложными фруктовыми, цветочными и цитрусовыми тонами, поиск оригинального вкуса и возрождение старинных рецептур служат

предпосылками для применения сухого охмеления различными предприятиями данного сегмента пищевой промышленности (Dresel et al., 2015; Schönberger & Kostelecky, 2011; Гернет, 2017; Ko-ростелев и др., 2020; Матвеева & Титов, 2014; Фе-доренко & Трэнкл, 2017). Сухое охмеление может стать источником ул учшения сенсорных характеристик пива с вероятностью потенциального увеличения содержания в пивоваренной продукции соединений, оказывающих благоприятное воздействие на организм человека, позволя-

ет в большей степени раскрыть потенциал хмеля (Donaldson, Bamforth, & Heymann, 2012; Barry, Muggah, McSweeney, & Walker, 2018; Kippensberger, Hanke, Biendl, Stettner, & Lagemann, 2014; Zanoli & Zavatti, 2008).

Известны данные исследователей о специальных технологиях внесения хмеля и соответствующем оборудовании (Бородулин, Сафонова, Просин, & Миленький, 2017; Матвеева, & Титов, 2014; Матвеева, & Титов, 2015; van Opstaele, de Rouck, de Clippeleer, de Cooman, 2010). Однако недостаточно раскрыты вопросы по оптимизации условий и продолжительности динамического процесса охмеления с помощью математического аппарата.

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

Материалы и методы исследования Материалы

Объектом исследования являлось светлое пиво верхового брожения, плотность начального сусла 11 %; содержание этанола 5,2 %; рН 3,82.

Применяли два сорта хмеля - Saphir и Chinook урожая 2020 года с содержанием компонентов: Saphir (а-кислота 3,5 %); Chinook (а-кисло-та 15,0 %).

Оборудование

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

Методы

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

и их относительные концентрации в хмеле и пиве определяли с помощью метода высокоэффективной жидкостной хроматографии (ВЭЖХ) (Нао, et al., 2018; Haseleu, Intelmann, Hofmann, 2009; Чес-нокова, 2011). Соединения выделяли из пива следующим образом: образец пива в количестве 5 см3 подкисляли ортофосфорной кислотой (100 см3) и экстрагировали изооктаном (10 см3) в течение 30 мин. Экстракт переносили в стеклянную пробирку и выпаривали. К остатку после выпаривания добавляли ацетонитрил (2 см3) и анализировали методом ВЭЖХ в трех повторностях. Для кислот хмеля и их производных были получены градуировочные графики на основе стандартов в диапазоне 1, 5, 10, 20,40, 60 мг/дм3. Общее содержание полифенолов определяли спектрофо-тометрическим методом (Чеснокова, 2011).

Оценки состава и свойств исследуемых объектов осуществляли в соответствии с нормативной документацией: методы отбора проб - по ГОСТ 12786-801; определение спирта - по ГОСТ 12787-812; определение экстрактивности начального сусла - по ГОСТ 12787-813; определение кислотности - по ГОСТ 12788-874; определение рН - по ГОСТ 31764-20125; определение цвета - по ГОСТ 12789-876; определение органолептических показателей - по ГОСТ 30060-937.

Процедура исследований

В аппарат для сухого охмеления вносили гранулы хмеля (Saphir или Chinook, расход 100 г./гл), заполняли его пивом, после чего с помощью насоса проводили рециркуляцию суспензии. Продолжительность процесса до 8-10 часов. Во избежание вспенивания, температуру пива не поднимали выше 3 °С; а разность давления между аппаратом и емкостью не превышала 0,2 атм (Новикова, Рукавицын, Муравьев, 2019).

Анализ данных

Для статистической оценки исследований применяли специальное ПО - Design Expert, v. 11.1.0, Stat-Ease Inc. и функции MS Excel. Экспериментальные данные обрабатывали вероятностными методами математической ста-

1 ГОСТ 12786-80. (2013). Пиво. Правила приемки и методы отбора проб. М.: Стандартинформ.

2 ГОСТ 12787-80. (2013). Пиво. Методы определения спирта, действительного экстракта и расчет сухих веществ в начальном сусле. М.: Стандартинформ.

3 Там же.

4 ГОСТ 12788-87. (2015). Пиво. Методы определения кислотности. М.: Стандартинформ.

5 ГОСТ 31764-2012. (2013). Пиво. Метод определения рН. М.: Стандартинформ.

6 ГОСТ 12789-87. (2011). Пиво. Методы определения цвета. М.: Стандартинформ.

7 ГОСТ 30060-93. (2013). Пиво. Методы определения органолептических показателей и объема продукции. М.: Стандартинформ.

тистики - дисперсионным и корреляционным анализом (Грачев & Плаксин, 2005).

Для аппроксимации зависимостей показателей от времени использовали радиальные нейронные сети (Bishop, 1995; Шелехова, Скворцова, Полтавская, 2021).

Результаты и их обсуждение

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

Факторами оптимизации являются: t - длительность процесса охмеления; Т - температура процесса охмеления; с - концентрация хмеля Chinook в комбинации хмелей Saphir и Chinook (представляет собой вещественное число от 0 до 1).

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

позволит получить плавные функции при изменении времени, температуры или состава хмелей.

Разработанная и изложенная ниже методика прогноза позволяет определить следующие физико-химические показатели пива: сг- концентрация гумулинонов; сиак - концентрация изо-а-кислот; сак - концентрация а-кислот; спф -общее содержание полифенолов; рН - водородный показатель (показатель кислотности); р - объемная плотность напитка.

Для заданной комбинации хмелей Saphir и Chinook можно определить показатели сенсорной оценки в приближении линейной комбинации с учетом концентрации хмелей: Иах - интенсивность аромата хмеля; Ках - качество аромата хмеля; Пв - полнота вкуса; Кг - качество горечи; Псз - предпочтение (среднее значение).

На основе спрогнозированных для данного набора t, Т, с физико-химических показателей и показателей сенсорной оценки определяются частные критерии оптимизации: Кфхп - критерий оптимизации физико-химических показателей; ^ - критерий оптимизации сенсорной оценки (Новикова и др., 2018).

На основе частных критериев оптимизации определяется общий критерий оптимизации Коп, который представляет собой интегральную оценку качества пива. Таким образом общий критерий оптимизации является функцией трех технологических параметров: K^t, Т, с).

Оптимизация заключается в определении значений или диапазонов (что лучше для практических

Показатели качесгва

Гех но леи л-je ские параметры

Пр 0ГН03 влияния ил свойств л га ¡в я

Расчет сенсорном оценки в линемном пр п&шжши II

■-ia^iHbi? Kpirreprni оптимизации

Общий Kprrrepiti"! опгпмкацш!

Рисунок 1. Постановка задачи оптимизации технологических параметров процесса сухого охмеления

целей) параметров t, Т, с, при которых функция принимает максимальное значение K^t, Т, с).

Коп {t>T,c) —> max => topt,Topt,copt,

где topt, Т0pt, copt - оптимальные значения технологических параметров.

Собранные в процессе экспериментального исследования данные позволяют получить функцию Kon(t, Т, с) в аналитическом виде на основе линейной и нейросетевой аппроксимаций. По известной функции Kon(t, Т, с) численным методом производится поиск максимума и определяются оптимальные параметры.

Составление общего критерия оптимизации из частных критериев

В соответствии с часто используемым подходом общий критерий оптимизации К„п представляет собой линейную комбинацию двух частных критериев оптимизации К^ и К,,0:

Коп = wéxn -Кфхп + Wcc

где Нфи, и 1УС0 - весовые коэффициенты, которые могут принимать вещественные значения от О до 1.

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

оценке качестве напитка. Весовые коэффициенты определены методом экспертных оценок и составляют W^ = 0,5, Wm = 0,5.

Частные критерии рассчитываются по шести (Кфхп) и по пяти (Ксо) показателям путем линейной комбинации показателей с весовыми коэффициентами и нормировкой на характерные значения показателей:

мР

Кч =

i=1

1

Pil-Pü

Pi-Pib di = 1;

Pa-Pi, dt=-1,

ще К, - частный критерий оптимизации (К^ или Кео); г - номер показателя; Л/р - количество показателей; м^ - весовой коэффициент г-го показателя; Р1 - показатели, на основе которых рассчитывается частный критерий оптимизации (сг, сиак, и т. д.). Рп и Ра - минимальное и максимальное значения показателя.^ в экспериментальном исследовании; с- коэффициент благоприятной тенденции (равен 1, если необходимо добиться максимума показателя и -1, если необходимо добиться минимума показателя).

Весовые коэффициенты щ, определенные экспертным путем и отражающие относительную важность показателей, приведены в Таблицах 1 и 2 (Коростелев и др., 2020). Сумма весовых коэффициентов в каждой группе равна 1.

Аппроксимация экспериментальных данных

Для определения оптимальных технологических параметров используются функции К^, Т, с) - зависимость физико-химических показателей Щ

Таблица 1

Весовые коэффициенты физико-химических показателей

Показатель Обозначение Весовой коэффициент w, (относительная важность) Оптимальные значения Коэффициент благоприятной тенденции dt

Концентрация гу-мулинонов Cr 0,20 max 1

Концентрация изо-а-кислот ^иок 0,40 min -1

Концентрация а-кислот ^ах 0,15 2,5 ± 0,5 0

Общее содержание полифенолов Спф 0,15 max 1

Водородный показатель pH 0,05 4,7 ±0,1 0

Плотность - 0,05 1,0045-1,0055 0

Таблица 2

Весовые коэффициенты показателей сенсорной оценки

Показатель Обозначение Весовой коэффициент Wj (относительная важность) Оптималь- Коэффициент благо-ные значения приятной тенденции di

Интенсивность аромата хмеля и» 0,10 1

Качество аромата хмеля Кцх 0,25 1

Полнота вкуса П, 0,15 max 1

Качество горечи к„ 0,25 1

Предпочтение (ср. знач.) Общая оценка Псз 0,25 1

(в частности, сг, сиак, сак, спф, pH, р) от длительности t процесса сухого охмеления, температуры Т, концентрации с хмеля Chinook в комбинации хмелей Saphire + Chinook.

Ранее в работе (Новикова, Рукавицын, Муравьев, 2018) представлены экспериментальные данные по зависимостям K^t, Т, с). Для каждой функции получены 24 экспериментальные точки: для шести значений длительности процесса t (0...12 ч), двух значений температуры Г (1,5 и 3,0 °С) и двух значений концентрации с (0 и 1, то есть только хмель Saphire и только Chinook).

Так как 24 экспериментальные точки расположены сравнительно редко в пространстве факторов оптимизации (£, Т, с), было принято решение использовать аналитическую аппроксимацию для того, чтобы частные критерии оптимизации Kfi, Т, с) были непрерывными функциями в пространстве (t, Т, с).

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

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

При расчете показателя K^t, Т,с)в первую очередь производили линейную интерполяцию по переменной с, то есть по составу комбинированного хмеля:

Kt (t,T) = К{ (t,T,0) + с[К( (t,T,l) Ц К, (*,Г,0)],

где Ki(t, Т) - показатель i для заданной концентрации с в комбинации хмелей Saphir и Chinook; K^t,

T, 0) и Ki(t, T, 1) - показатель i для хмелей Saphir и Chinook.

Далее производили линейную интерполяцию данных по температуре процесса Т:

Kt (?) = К, (/,1.5) + (',3.0) Ц К, (,,1.5)],

где 1,5 и 3,0 - значения температуры, при которых проводили эксперименты по сухому охмелению.

Для дальнейшего получения аналитической зависимости для функции Kt(t) (при заданном составе комбинации хмелей с и температуре Т) использовали радиальную нейронную сеть.

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

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

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

программа «Программа для прогнозирования результатов процесса сухого охмеления» на языке Object Pascal в среде программирования Delphi 7.

Пример аппроксимации показателя cHaK(t) для сухого охмеления хмелем Chinook (с = 1) при 1,5 °С (Т= 1,5) показан на Рисунке 2. График аппроксимирующей функции практически проходит близко к шести экспериментальным точкам сиа1д(^), при этом кривая является плавной, начало и конец кривой не стремятся к бесконечности (как обычно происходит при полиномиальной аппроксимации).

С использованием формул и получены аппроксимирующие функции не только для зависимости показателя от времени Щф, но и от концентрации с в комбинации хмелей и температуры процесса Т, то есть полноценные функции K^t, Т, с). Для примера на Рисунке 3 приведена серия графиков функций cr(t, Т= 1,5 °С, с).

Таким образом, разработанный математический аппарат - комбинация линейной интерполяции и аппроксимации радиальной нейронной сетью -позволяют получить плавные функции при изме-

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

Решение общей задачи оптимизации

Разработанный выше математический аппарат позволяет определить общий критерий оптимизации Kon(t, Т, с) для различных комбинаций параметров £, Т, с. Основной задачей оптимизации является поиск такой комбинации параметров t, Т, с, при которой критерий Коп принимает наибольшее значение. Для поиска максимума Коп в пространстве (t, Т, с) разработана компьютерная программа на языке Object Pascal в среде программирования Delphi 7. В программе производится поиск максимума Кои численным методом - методом упорядоченного перебора с равным шагом по трем переменным. Время поиска составляет около 2 с машинного времени. По каждой переменной исследуемый диапазон разбили на 101 уровень (переменная цикла изменяется от 0 до 100). В ходе перебора переменная t изменяется от 0 до 12 часов с шагом 0,5 ч.; переменная Г изменяется от 1,5 до 3,0 °С с шагом 0,15 °С; переменная с изменяется от 0,00 до 1,00 с шагом 0,01.

» —»

V.

О 2 4 6 8 10 12 г. ч

Рисунок 2. Пример аппроксимации экспериментальных данных радиальной нейронной сетью: зависимость концентрации изо-а-кислот в процессе сухого охмеления хмелем Chinook при 1,5 °С

Ст.-

иг/дм3 60-

40

20

0

О 2 4 б 8 10 12 Г. ч

Рисунок 3. Зависимость концентрации гумулинонов сг от длительности процесса сухого охмеления при температуре 1,5 °С при использовании различных комбинаций хмелей 1-е Saphire + с Chinook

Оптимальная точка обычно может располагаться либо внутри исследуемой области факторного пространства, либо на границе. В данном случае точка с максимальным общим критерием Коп = = 0,679 расположена на границе исследуемой области: t= 12 ч; Т= 3,0 °С; с = 0. Таким образом, наилучшие физико-химические показатели и показатели сенсорной оценки достигаются при сухом охмелении хмелем Saphire в течении 12 ч при температуре 3,0 °С.

Для более детального анализа поверхности отклика построены три серии картограмм (Рисунок 4). На картограммах отмечена точка оптимума обозначением «opt».

Использование комбинации хмелей Saphire + Chinook не дает улучшение общего критерия оптимизации по сравнению с использованием только хмеля Saphire (исчезает красная и оранжевая области на Рисунок 4, б). Однако, если необходимо использовать только хмель Chinook, оптимальная длительность процесса составляет 10,5 ч (точка лежит уже не на границе диапазона); оптимальная температура составляет 3,0 °С, однако при любых температурах критерий оптимизации практически одинаков и составляет 0,514 (Рисунок 4, в).

При температурах ниже 3,0 °С критерий оптимизации получается незначительно ниже, чем для

3,0 °С (Рисунок 5). Так, при температурах 1,5 и 2,0 °С удается достигнуть максимального критерия оптимизации 0,623 и 0,638 соответственно при длительности процесса 8,15 (Рисунок 5, а) и 11,25 (Рисунок 5, б) часов соответственно.

Аналогичные исследования других авторов показали важность контроля температуры пива при проведении сухого охмеления (Titus, et al., 2021; Матвеева & Титов, 2015).

Даже при небольшой длительности процесса 4 часов удается достигнуть значений критерия оптимизации более 0,60 (оранжевая область на Рисунке 6, а). Однако только при длительности процесса порядка 12 часов удается достигнуть

критерия оптимизации более 0,65 (красная область на Рисунке 6, в).

С точки зрения экономической эффективности желательно, чтобы длительность процесса сухого охмеления была минимальной при высоких показателях качества (Cibaka, et al., 2018; Maye, Smith, Leker, 2016; Матвеева & Титов, 2015). Поэтому более подробно изучено изменение с течением времени t критерия оптимизации Коп, а именно, максимального значения KonM(t) среди различных температур (Г от 0 до 20 °С) и различных концентраций комбинации хмелей (с от 0 до 1).

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

Рисунок 4. Картограммы оптимизации Kon(t, Т) при разных концентрациях с хмеля Chinook в комбинации хмелей (-lc) Saphire + с Chinook

Рисунок 5. Картограммы оптимизации Коп(1, с) при разных температурах Г процесса сухого охмеления

Т.. с) в..К^СП:Гс)

Рисунок 6. Картограммы оптимизации Кои(Т, с) при разной длительности Г процесса сухого охмеления

(Рисунок 7). Точка А - это точка наименьшей длительности процесса 3,8 ч, при которой достигается высокие показатели качества пива - критерий оптимизации равен 0,636. Благодаря такой малой длительности процесса может быть достигнут наибольший экономический эффект.

Вторая характерная точка В представляет собой локальную точку максимума и соответствует примерно вдвое большей длительности процесса 6,2 ч, при этом достигается больший критерий оптимизации 0,650. С точки зрения экономического эффекта в данном случае необходимо оценивать баланс между снижением производительности в 1,75 раза процесса и повышением качества на 2,2%. Данный вариант не позволяет получить пиво заметно высокого качества, при высокой длительности, поэтому не рекомендуется для реального производства.

Третья характерная точка С соответствует глобальному максимуму Коп" = 0,679 и требует значительной длительности процесса 13,1 ч. В таком технологическом процессе получается наиболее качественное пиво (на 6,8 % качественнее, чем для

К,Л/) 0,7 0.6 0.5 0.4

0 2 4 б 8 10 12 Л ч

Рисунок 7. Зависимость от времени С критерия оптимизации Коим, максимального в диапазонах температур Гот 1,5 до 3,00 и концентраций с 0-1

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

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

Выводы

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

Определены оптимальные технологические параметры процесса сухого охмеления: хмель Saphire, длительность процесса 12 ч, температура 3,0 °С. Использование комбинации хмелей Saphire + Chinook не дает улучшение общего критерия оптимальности по сравнению с использованием только хмеля Saphire. При сухом охмелении хмелем Chinook оптимальная длительность процесса составляет 10,5 ч, оптимальная температура составляет 3,0 °С.

При температурах процесса сухого охмеления ниже 3,0 °С критерий оптимизации получает-

ся незначительно ниже, чем для 3,0 °С. При температурах 1,5 и 2,0 °С оптимальная длительность процесса составляет 8,15 и 11,25 суток соответственно. Уже при небольшой длительности процесса 4 ч удается достигнуть значений критерия оптимизации более 0,60, однако только при длительности порядка 12 ч удается достигнуть критерия оптимизации более 0,65.

С точки зрения экономической эффективности возможны два варианта процесса сухого охмеления. Для получения пива обычного класса с достаточно высоким показателем качества 0,636 длительность процесса должна составлять 3,8 ч. Для получения пива премиального класса с показателем качества 0,679 (на 6,8 % качественнее, чем обычного класса) длительность процесса должна составлять 13,1 ч.

Литература

Бородулин, Д. М., Сафонова, Е. А., Просин, М. В., & Миленький, И. О. (2017). Исследование процесса охмеления пивного сусла с применением современного оборудования. Современные материалы, техника и технологии, 3 (11), 16-21. Гернет, М. В. (2017). Перспективы расширения ассортимента напитков брожения для пивоваренных заводов малой мощности. Пиво и напитки, 3,14-17. Грачев, Ю. П., & Плаксин, Ю. М. (2005). Математические методы планирования экспериментов. М.: ДеЛи Принт. Коростелев, А. В., Рукавицын, П. В., Новикова, И. В., Кучменко, Т. А., Умарханов, Р. У., & Муравьев, А. С. (2020). Исследование ароматобразую-щих компонентов хмеля с применением химических сенсоров. Известия Вузов. Прикладная Химия и Биотехнология, 10(3), 479-486. ЬИрв:// doi.org/10.21285/2227-2925-2020-10-3-479-486 Матвеева, Н. А., & Титов, А. А. (2014). Выбор сорта хмеля для технологии сухого охмеления. Научный журнал НИУ ИТМО. Процессы и аппараты пищевых производств, 4,120-125. Матвеева, Н. А., & Титов, А. А. (2015). Применение технологии сухого охмеления в пивоварении. Научный журнал НИУ ИТМО. Процессы и аппараты пищевых производств, 1, 111-118. Меледина, Т. В. (2003). Сырье и вспомогательные

материалы в пивоварении. СПб.: Профессия. Новикова, И. В., Агафонов, Г. В., Коротких, Е. А., Коростелев, А. В., Рукавицын, П. В., & Новиков, С. В. (2018). Разработка методики проектирования безалкогольных напитков с учетом сбалансированности состава. Пиво и напитки, 2,26-30.

Новикова, И. В., Рукавицын, П. В., & Муравьев, А. С.

(2018). Влияние основных параметров процесса сухого охмеления на физико-химические показатели пива. Технология и товароведение инновационных пищевых продуктов, 6,9-17.

Новикова, И. В., Рукавицын, П. В., & Муравьев, А. С.

(2019). К вопросу перехода ароматических соединений хмеля в пиво при реализации сухого охмеления. Пищевая промышленность, 1,69-73.

Осовский, С. (2004). Нейронные сети для обработки

информации. М.: Финансы и статистика. Федоренко, Б. Н., & Трэнкл, Ф. (2017). Инновационное технологическое оборудование для крафто-вого пивоварения. Пиво и напитки, 1,22-24. Чеснокова, А. Н. (2011). Фенольные соединения хмеля обыкновенного (Humulus Lupulus L.) и хмеле-продуктов [Кандидатская диссертация, Иркутский государственный технический университет]. Иркутск, Россия. Шелехова, Н. В., Шелехова, Т. М., Скворцова, Л. И., & Полтавская, Н. В. (2021). Идентификация спиртных напитков с применением искусственных нейронных сетей. Пищевая промышленность, 3, 43-48. https://doi.org/10.24412/0235-2486-2021-3-0028 Barry, S., Muggah, Е. М., McSweeney, М. В., & Walker, S. (2018). A preliminary investigation into differences in hops' aroma attributes. International Journal of Food Science & Technology, 55(3), 804-811. https ://doi.org/10.1111/ijfs. 13656 Bishop, С. M. (1995). Neural networks for pattern

recognition. Clarendon Press. Cibaka, M.-L. K., Ferreira, C. S., Decourriere, L., Lorenzo-Alonso, C.-J., Bodart, E., & Collin, S. (2018). Dry hopping with the dual-purpose varieties amarillo, citra, hallertau blanc, mosaic, and sorachi ace: Minor contribution of hop terpenol glucosides to beer flavors. Journal of the American Society of Brewing Chemists, 75(2), 122-129. https://doi.org/10.1094/ASBCJ-2017-2257-01 Donaldson, B. R., Bamforth, C. W., & Heymann, H. (2012). Sensory descriptive analysis and free-choice profiling of thirteen hop varieties as whole cones and after dry hopping of beer. Journal of the American Society of Brewing Chemists, 70(3), 176-181. https://doi.org/10.1094/ASBCJ-2012-0710-01 Dresel, M., Praet, Т., van Opstaele, F., van Holle, A., Naudts, D., Keukeleire, D. de, Cooman, L. de, & Aerts, G. (2015). Comparison of the analytical profiles of volátiles in single-hopped worts and beers as a function of the hop variety. BrewingScience, 68(1-2), 8-28. Hao, J., Dong, J., Yin, H., Yan, P., Ting, P. L., Li, Q., Tao, X., Yu, J., Chen, H., & Li, M. (2018). Optimum method of analyzing hop derived aroma compounds in beer by headspace solid-phase microextraction (SPME) with GC/MS and their evolutions during chínese

lager brewing process. Journal of the American Society of Brewing Chemists, 72(4), 261-270. https://doi. org/10.1094/ASBCJ-2014-1021-01 Haseleu, G., Intelmann, D., & Hofmann, T. (2009). Identification and RP-HPLC-ESI-MS/MS quantitation of bitter-tasting beta-acid transformation products in beer. Journal of Agricultural and Food Chemistry, 57(16), 7480-7489. https://doi. org/10.1021/jf901759y Kippensberger, M., Hanke, S., Biendl, M., Stettner, G., & Lagemann, A. (2014). Transfer of nitrate and various pesticides into beer during dry hopping. BrewingScience, 67,1-9. Leitao, C., Marchioni, E., Bergaentzle, M., Zhao, M., Didierjean, L., Taidi, B., & Ennahar, S. (2011). Effects of processing steps on the phenolic content and antioxidant activity of beer. Journal of Agricultural and Food Chemistry, 59(4), 1249-1255. https://doi.org/10.1021/jfl04094c Maye, J. P., Smith, R., & Leker, J. (2016). Humulinone formation in hops and hop pellets and its implications for dry hopped beers. Technical Quarterly, 55(1), 23-27. https://doi.org/10.1094/TQ-53-1-0227-01

Schönberger, C., & Kostelecky, T. (2011). 125th Anniversary review. The role of hops in brewing. Journal of the Institute of Brewing, 117(5), 259-267. https://doi.Org/10.1002/j.2050-0416.2011.tb00471.x

Takoi, K., Tokita, K., Sanekata, A., Usami, Y., Itoga, Y., Koie, K., Matsumoto, I., & Nakayama, Y. (2016). Varietal difference of hop-derived flavour compounds in late-hopped/dry-hopped beers. BrewingScience, 69,1-7. Titus, B. M., Lerno, L. A., Beaver, J. W., Byrnes, N. K., Heymann, H., & Oberholster, A. (2021). Impact of dry hopping on beer flavor stability. Foods, 10(6), 1264. https://doi.org/10.3390/foodsl0061264 Van Opstaele, F., de Rouck, G., de Clippeleer, ]., de Cooman, L. (2010). Analytical and Sensory Assessment of Hoppy Aroma and Bitterness of Conventionally Hopped and Advanced Hopped Pilsner Beers. Journal of the Institute of Brewing, 116(4), 445-458. https://doi. org/10.1002/j.2050-0416.2010.tb00796.x Walsh, D. B., O'Neal, S. D., George, A. E., Groe-nendale, D. P., Henderson, R. E., Groe-nendale, G. M., & Hengel, M.J. (2018). Evaluation of pesticide residues from conventional, organic, and nontreated hops on conventionally hopped, late-hopped, and wet-hopped beers. Journal of the American Society of Brewing Chemists, 74(1), 53-56. https://doi.org/10.1094/ASBCJ-2016-1115-02 Zanoli, P., & Zavatti, M. (2008). Pharmacognostic and pharmacological profile of Humulus lupulus L. Journal of Ethnopharmacology, 116(3), 383-396. https://doi.Org/10.1016/j.jep.2008.01.011

Optimization of Technological Parameters of the Dry Hopping Process With the Development of a Methodology for the Integral Assessment of the Quality of Beer

Inna V. Novikova

Voronezh State University ofEngineer.ng Technologies 19, Revolution Avenue, Voronezh, 394036, Russian Federation

E-mail: noviv@list.ru

Pavel V. Rukavitsyn

Voronezh State University of Engineering Technologies 19, Revolution Avenue, Voronezh, 394036, Russian Federation

E-mail: kafedra-tbisp@mail.ru

Alexander S. Muravyov

Voronezh State University of Engineering Technologies, 19, Revolution Avenue, Voronezh, 394036, Russian Federation

E-mail: hntrun@mail.ru

The development of craft beers with complex fruit, floral and citrus tones, the search for an original taste and the revival of old recipes are the prerequisites for the use of dry-hoping in the brewing industry. The dry-hopping process was carried out dynamically in an experimental setup with a hop apparatus, screen filter, fermentation tank and pump. Saphir and Chinook hops were used. The study included determination of physical and chemical parameters of beer: concentration of humulinones, iso-a-acids and a-acids; total content of polyphenols; pH; volume density. Sensory evaluation indicators included: hop aroma intensity and quality; flavor fullness; bitterness quality; and overall preference. Radial neural networks were used to approximate the results of the study. Optimal process parameters were determined: Saphir hops, process duration 12 hours, temperature 3,0 °C. The optimum process duration is 10,5 h and the optimum temperature is 3,0 °C. Processing temperatures below 3.0 °C result in a slightly lower optimization criterion than the 3.0 °C processing temperature. At temperatures of 1.5 and 2.0 °C the optimum process duration is 8.15 and 11.25 days respectively. The use of the combination of Saphir & Chinook hops offers no improvement in the overall optimality criterion compared with the use of Saphir alone. The process duration of 4 hours achieves an optimization criterion of more than 0.60, but only with durations of around 12 hours does an optimization criterion of more than 0.65 become attainable.

Keywords: hops, neural networks, Saphire, Chinook, optimization criterion, dry hopping

References

Borodulin, D. M., Safonova, E. A., Prosin, M. V., & Milen'kii, I. O. (2017). Issledovanie protses-sa okhmeleniya pivnogo susla s primeneniem sovremennogo oborudovaniya [Investigation of the process of hopping beer wort with the use of modern equipment]. Sovremennye materialy, tekh-nika i tekhnologii [Modern Materials, Equipment and Technologies], 5,16-21.

Chesnokova, A. N. (2011). Fenol'nye soedineniya khme-lya obyknovennogo (Humulus Lupulus L.) i khmele-produktov [Candidate Dissertation, Irkutskii go-

sudarstvennyi tekhnicheskii universitet]. Irkutsk, Russia.

Fedorenko, B. N., & Trenkl, F. (2017). Innovatsionnoe tekhnologicheskoe oborudovanie dlya kraftovogo pivovareniya [Innovative process equipment for craft brewing]. Pivo i napitki [Beer and drinkj, 1,22-24.

Gernet, M. V. (2017). Perspektivy rasshireniya assor-timenta napitkov brozheniya dlya pivovarenny-kh zavodov maloi moshchnosti [Prospects for expanding the range of fermentation drinks for small capacity breweries]. Pivo i napitki [Beer and Beverages ], 3,14-17.

Grachev, Yu. P., & Plaksin, Yu. M. (2005). Matema-ticheskie metody planirovaniya eksperimentov

[Mathematical methods for planning experiments]. Moscow: DeLi print.

Korostelev, A. V., Rukavitsyn, P. V., Novikova, I. V., Kuchmenko, T. A., Umarkhanov, R. U., & Murav'-ev, A. S. (2020). Issledovanie aromatobrazu-yushchikh komponentov khmelya s primeneniem khimicheskikh sensorov [Study of aroma-forming components of hops using chemical sensors]. Izvestiya Vuzov. Prikladnaya Khimiya i Biotekhnologiya [Proceedings of Universities. Applied Chemistry & Biotechnology], 10(3), 479-486. https://doi.org/10.21285/2227-2925-2020-10-3-479-486

Matveeva, N. A., & Titov, A. A. (2014). Vybor sorta khmelya dlya tekhnologii sukhogo okhmeleniya [Selection of hop variety for dry hop technology]. Nauchnyi zhurnal NIU ITMO. Protsessy i apparaty pishchevykh proizvodstv [Scientific journal of ITMO University. Processes and apparatus for food production], 4,120-125.

Matveeva, N. A., & Titov, A. A. (2015). Primenenie tekhnologii sukhogo okhmeleniya v pivovarenii [Application of dry hops technology in brewing]. Nauchnyi zhumal NIU ITMO. Protsessy i apparaty pishchevykh proizvodstv [Scientific journal of ITMO University. Processes and apparatus for food production], 1,111-118.

Meledina, T. V. (2003). Syr'e i vspomogatel'nye materi-aly v pivovarenii [Raw materials and auxiliary materials in brewing]. S-Petersburg: Professiya.

Novikova, I. V., Agafonov, G. V., Korotkikh, E. A., Korostelev, A. V., Rukavitsyn, P. V., & Novikov, S. V. (2018). Razrabotka metodiki proektirovaniya be-zalkogol'nykh napitkov s uchetom sbalansirovan-nosti sostava [Development of design methodology for non-alcoholic beverages taking into account the balance of composition]. Pivo i napitki [Beer and Beverages], 2,26-30.

Novikova, I. V., Rukavitsyn, P. V., & Murav'ev, A. S.

(2018). Vliyanie osnovnykh parametrov protsessa sukhogo okhmeleniya na fiziko-khimicheskie po-kazateli piva [Influence of the main parameters of the dry hops process on physico-chemical parameters of beer]. Tekhnologiya i tovarovedenie inno-vatsionnykh pishchevykh produktov [Technology and Merchandising of Innovative Foods], 6,9-17.

Novikova, I. V., Rukavitsyn, P. V., & Murav'ev, A. S.

(2019). K voprosu perekhoda aromaticheskikh soedinenii khmelya v pivo pri realizatsii sukhogo okhmeleniya [On the transition of hop aromatic compounds into beer during the implementation of dry hops]. Pishchevaya promyshlennost' [Food Industry], 1,69-73.

Osovskii, S. (2004). Neironnye seti dlya obrabotki infor-matsii [Neural networks for information processing]. Moscow: Finansy i statistika.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Shelekhova, N. V., Shelekhova, T. M., Skvortsova, L. I., & Poltavskaya, N. V. (2021). Identifikatsiya spirt-nykh napitkov s primeneniem iskusstvennykh nei-ronnykh setei [Identification of alcoholic beverages using artificial neural networks]. Pishchevaya promyshlennost' [Food Industry], 3, 43-48. https:// doi.org/10.24412/0235-2486-2021-3-0028 Barry, S., Muggah, E. M., McSweeney, M. B., & Walker, S. (2018). A preliminary investigation into differences in hops' aroma attributes. International Journal of Food Science & Technology, 53(3), 804-811. https ://doi.org/10.1111/ijfs. 13656 Bishop, C. M. (1995). Neural networks for pattern recognition. Clarendon Press. Cibaka, M.-L. K., Ferreira, C. S., Decourriére, L., Lorenzo-Alonso, C.-J., Bodart, E., & Collin, S. (2018). Dry hopping with the dual-purpose varieties amarillo, citra, hallertau blanc, mosaic, and sorachi ace: Minor contribution of hop terpenol glucosides to beer flavors. Journal of the American Society of Brewing Chemists, 75(2), 122-129. https:// doi.org/10.1094/ASBCJ-2017-2257-01 Donaldson, B. R., Bamforth, C. W., & Heymann, H. (2012). Sensory descriptive analysis and free-choice profiling of thirteen hop varieties as whole cones and after dry hopping of beer. Journal of the American Society of Brewing Chemists, 70(3), 176-181. https://doi.org/10.1094/ASBCJ-2012-0710-01 Dresel, M., Praet, T., van Opstaele, F., van Holle, A., Naudts, D., Keukeleire, D. de, Cooman, L. de, & Aerts, G. (2015). Comparison of the analytical profiles of volátiles in single-hopped worts and beers as a function of the hop variety. BrewingSdence, 68(1-2), 8-28. Hao, J., Dong, J., Yin, H., Yan, P., Ting, P. L., Li, Q., Tao, X., Yu, J., Chen, H., & Li, M. (2018). Optimum method of analyzing hop derived aroma compounds in beer by headspace solid-phase mi-croextraction (SPME) with GC/MS and their evolutions during Chinese lager brewing process. Journal of the American Society of Brewing Chemists, 72(4), 261-270. https://doi.org/10.1094/ ASBCJ-2014-1021-01 Haseleu, G., Intelmann, D., & Hofmann, T. (2009). Identification and RP-HPLC-ESI-MS/MS quantitation of bitter-tasting beta-acid transformation products in beer. Journal of Agricultural and Food Chemistry, 57(16), 7480-7489. https://doi. org/10.102 l/jf901759y Kippensberger, M., Hanke, S., Biendl, M., Stettner, G., & Lagemann, A. (2014). Transfer of nitrate and various pesticides into beer during dry hopping. BrewingSdence, 67,1-9. Leitao, C., Marchioni, E., Bergaentzlé, M., Zhao, M., Didierjean, L., Taidi, B., & Ennahar, S. (2011). Effects of processing steps on the phenolic content and antioxidant activity of beer. Journal of

Agricultural and Food Chemistry, 59(4), 1249-1255. https://doi.org/10.1021/jfl04094c Maye, J. P., Smith, R., & Leker, J. (2016). Humulinone formation in hops and hop pellets and its implications for dry hopped beers. Technical Quarterly, 55(1), 23-27. https://doi.org/10.1094/TQ-53-l-0227-01 Schönberger, C., & Kostelecky, T. (2011). 125th Anniversary review. The role of hops in brewing. Journal of the Institute of Brewing, 117(5), 259-267. https://doi.Org/10.1002/j.2050-0416.2011.tb00471.x Takoi, K., Tokita, K., Sanekata, A., Usami, Y., Itoga, Y., Koie, K., Matsumoto, I., & Nakayama, Y. (2016). Varietal difference of hop-derived flavour compounds in late-hopped/dry-hopped beers. Bre-wingScience, 69,1-7. Titus, B. M., Lerno, L. A., Beaver, J. W., Byrnes, N. K., Heymann, H., & Oberholster, A. (2021). Impact of dry hopping on beer flavor stability. Foods, 10(6), 1264. https ://doi.org/10.3390/foods 10061264

Van Opstaele, F., de Rouck, G., de Clippeleer, J., de Cooman, L. (2010). Analytical and Sensory Assessment of Hoppy Aroma and Bitterness of Conventionally Hopped and Advanced Hopped Pilsner Beers. Journal of the Institute of Brewing, 116(4), 445-458. https://doi. org/10.1002/j.2050-0416.2010.tb00796.x Walsh, D. B., O'Neal, S. D., George, A. E., Groen-endale, D. P., Henderson, R. E., Groenendale, G. M., & Hengel, M. J. (2018). Evaluation of pesticide residues from conventional, organic, and nontreated hops on conventionally hopped, late-hopped, and wet-hopped beers. Journal of the American Society of Brewing Chemists, 74(1), 53-56. https://doi. org/10.1094/ASBCJ-2016-1115-02 Zanoli, P., & Zavatti, M. (2008). Pharmacognostic and pharmacological profile of Humulus lupulus L. Journal of Ethnopharmacology, 116(3), 383-396. https://doi.Org/10.1016/j.jep.2008.01.011

i Надоели баннеры? Вы всегда можете отключить рекламу.