Автоматизация контроля качества гомогенизации молока методами ультразвуковой спектроскопии Текст научной статьи по специальности «Математика»

Информационные технологии, моделирование и управление

УДК: 637.073:534.2

Профессор В.К. Битюков, профессор А.А. Хвостов, доцент Д.И. Ребриков, аспирант В.Е. Мерзликин

(Воронеж. гос. ун-т. инж. технол.) кафедра информационных и управляющих систем. тел. (473) 255-38-75

E-mail: [email protected], [email protected]

Professor V.K. Bityukov, professor A.A. Khvostov, associate Professor D.I. Rebrikov, graduate V.E. Merzlikin

(Voronezh state university of engineering technologies) Department of information and control systems. phone (473) 255-38-75

E-mail: [email protected], [email protected]

Автоматизация контроля качества гомогенизации молока методами

ультразвуковой спектроскопии

Automation of quality control of milk homogenization by ultrasonic spectroscopy methods

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

Summary. The paper deals with the possibility of determining homogenization degree of milk and dairy products using ultrasonic vibrations absorption spectra. Advantages of this method application in automated manufacturing systems were examined. Theoretical background of the method, as well as the possibility of determining the distribution of the fat globules in milk, depending on their sizes were substantiated. We derived mathematical equations, showing the relationship between the homogenization degree of dairy products and the acoustic properties of the medium, such as the propagation velocity of ultrasonic vibrations in the medium tested and the absorption coefficient in the medium at a specific frequency of the ultrasonic impact and medium temperature. It was shown that the measurement of the propagation of the milk acoustic properties in frequency, evaluation of their density function of the relaxation times spectrum, and then from the relaxation times to the particles masses allow operational control of the fat globules masses distribution in fractions. Theoretical studies were confirmed by experimental research carrying out whose results demonstrate clearly the dependence of absorption degree of ultrasonic vibrations on the degree of milk homogenization. The dependence between the estimates of the relaxation spectra and the first two moments of the statistical distribution of milk fat globules, which demonstrated their relationship was studied. Possible improvements the method under consideration to increase the reliability of the results obtained were proposed.

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

Keywords: dairy products, dairy products homogenization, acoustic properties, mass distribution, homogenization degree, fat globules,

ultrasonic vibrations absorption. _

© Битюков В.К., Хвостов А. А., Ребриков Д.И., Мерзликин В.Е., 2015

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

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

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

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

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

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

Наряду с прямыми методами контроля, использующими визуальное наблюдение и определение степени гомогенизации по специальным методикам, существует ряд косвенных методов, использующих связь поддающихся измерению физических величин, коррелирующих со степенью гомогенизации [2]. Одним из удобных для практического применения способов косвенного контроля степени гомогенизации и других свойств молочных продуктов являются ультразвуковые (УЗ) методы контроля. На базе этих методов созданы лабораторные анализаторы, использующиеся в молочной отрасли [3, 15, 16, 17].

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

Измерение степени гомогенизации с использованием УЗ основано на связи акустических свойств среды (скорости звука с, коэффициента поглощения а при определённой частоте УЗ воздействия и температуре), его вязко-упругих характеристик (компонентов комплексного модуля упругости ¿(о)) и структурно-молекулярных свойств. Непрерывное измерение акустических свойств реализуется с помощью УЗ излучателя и приемника, позволяющих излучать и принимать прошедшие через исследуемую среду акустические импульсы. По параметрам излученных и принятых импульсов рассчитывают акустические свойства среды [4]:

h 1 А

с = —, а = — ■ 1п- л

t к АР

где к - расстояние между поверхностями излучателя и приемника, м; / - время прохождения импульса, с; Аизл - амплитуда сигнала, подаваемого на излучатель; Апр - амплитуда сигнала

на приемнике, В.

Существующие УЗ измерители степени гомогенизации обычно основаны на интеграль-

ной оценке этого показателя по косвенным измерениям (вязкости, плотности, скорости осаждения и т.д.) и расчету степени гомогенизации по эмпирическим формулам [6].

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

Одной из областей применения методов акустической спектрометрии является оценка размеров частиц и их распределения по размерам [5, 6]. Акустические свойства среды а, с позволяют определять компоненты комплексного модуля упругости и их частотные характеристики [7]: ¿(ю) = Е' + гЕ" = Е' + гаг! =

1 -

= рс

ас а

1

^ а

■ + г ■ 2рс2

ас а

,.(1)

1+а

^ а

где а-частота, Гц; Е',Е" - действительная и мнимая части комплексного модуля упругости, Па; Г - динамическая вязкость, Па ■ с ; р - плотность, кг/м3.

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

Ё(со) =

ыг Е,а2т2 ] . у 1 1

у л 2 2

1=11 + а т2 1 1 /

(2)

+га

Г

Ыг

- у

¿—'л 2 2

;=11 + а т

1 / модуль

упругости;

где Е0 - равновесный Е1 - модуль упругости 1-го релаксационного механизма; г - нерелаксирующая вязкость; Г - вязкость 1-го релаксационного механизма; Т = Е^г - время релаксации 1-го релаксационного механизма.

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

Ё{со)=\ Н{т)

2 2 а т

1 + а2т2

d 1пт +

(3)

+г

I Н (т) -

-да 1

ат

2 2 + а т

^ 1пт

где Н(т) - функция плотности спектра времен

релаксации (релаксационный спектр).

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

Получение в аналитическом виде функции Н(т) возможно при аналитическом задании функций Е'(а), Е"(а) в интервале [-да,да].

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

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

+да

2

В другом источнике [7] приведен прямой метод расчета релаксационного спектра. Метод состоит в определении пар величин Ei-Hi. Пары Ei-Hi варьируются, и рассчитанные значения Е 'аос (а) и Е "са1с (а) сопоставляются с экспериментальными данными Е'ехр(а) и Е" (0).

Расхождение между экспериментом и результатами расчетов (получаемыми при варьировании пар Е-Н) оценивается функционалом ошибок. Эта величина может рассчитываться различным образом. Известен также [8] метод линеаризации, который основан на представлении зависимостей Е '(а) и Е "(а) или релаксационной функции, экспериментально измеренных в некоторых пределах в виде степенных рядов по переменной ш. В противоположность описанному выше методу, в данном случае минимизация функционала ошибок становится линейной задачей, и поэтому определение коэффициентов имеет однозначное решение.

Также в [7] приведен метод регуляризации, который обсуждался во многих работах [9, 10]. Этот метод представляет собой частный случай решения так называемых интегральных уравнений Фредгольма первого рода. Все обратные задачи определения времен релаксации относятся к задачам этого типа. Один из наиболее общепринятых методов решения таких задач основывается на методе регуляризации Тихонова [11]. Метод заключается в минимизации критерия У(Л):

(

Е -

г (<.0 ^

//

(4)

N 1

V (Л) = 1—

г =1 СТ

+Л\\Ш\\2,те 1Т

где X - параметр регуляризации; оператор L выбирается обычно из условий (Ш=Н) или (1И= Н"); а - ошибка; 1т - множество значений времен релаксации т.

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

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

иметь информацию об одном из компонентов модуля. Для жидкостей удобнее использовать модуль потерь Е". Одним из способов приближенного аналитического описания спектра времен релаксации по известному спектру Е "(а)

показан в [7].

По определению, плотность спектра времен релаксации Н(т), умноженная на От имеет смысл суммарной вязкости всех максвеллов-ских механизмов, времена релаксации которых заключены в интервале от т до т+ёт.

Способ заключается в использовании до-

ат

пущения о том, что выражение -

1 + ст

аппроксимируется ¿-функцией. В силу свойств ¿-функции:

I 5(х-х0) f (х)°х = f (х0) ,

(5)

можно записать мнимую часть комплексного модуля упругости как:

Е "(с) = аг](а) = а?]х+]

Н(т) атё(ст)

I л 2 2

ат 1 + 0т

Н| -а

(6)

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

Н(т) = о3(71(о3)-71х)\св=1. (7)

т

Пренебрегая г]х , в соответствии с [7]: Н (т) = ат?(а)1л = Е "(а)^. (8)

т т

Таким образом, для определения дискретного спектра времен релаксации, т.е. множества значений Н(т), необходимо измерить множество значений Е "(а>г). Значение модуля потерь может быть рассчитано с использованием формул (1), (2) по измеряемым акустическим параметрам среды: коэффициенту затухания и скорости распространения ультразвука.

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

2

В работах [12, 13] использованы подходы к оценке масс, вовлекаемых в колебательное движение по временам их релаксации. Одной из форм связи массы структурного элемента и его времени релаксации является вид:

т = k ■ тр , (9)

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

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

Для проведения исследований были проведены следующие действия (рисунок 1) [18, 19]:

Рисунок 1. Методика проведения экспериментальных исследований

1. Осуществление микросъемки образцов молока с разной глубиной обработки на гомогенизаторе и получение множества

R = ,г = 1,М,j = 1,N},где г - номер пикселя по горизонтали, 1 - номер пикселя по вертикали.

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

Ьр = ^рк,k = 1К: УЬрк = {Хгк,Тгк,Dк }},

где k,Хгк,Угк,Dk - порядковый номер, координаты и диаметр шарика.

3. Построение гистограмм, описывающих распределение жировой фазы по массам в исследуемом объеме.

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

5. Ультразвуковой анализ этих же образцов молока в доступном частотном диапазоне и построение частотных характеристик модуля потерь Е"(а) и функции плотности спектра

времен релаксации Н (т).

7. Изучение корреляции плотности спектра времен релаксации и аппроксимирующих гистограммы функций.

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

ЗЕ-12

2Е-12

2Е-12

2 2Е-12

о.

(11

I 1Е-12

го

а.

5Е-13 0Е-01 -5Е-13

< <

< с

[ 1 к < < :

Е 1 [ а \

< < 1 < 1 < I 1

30 60 90 120 150 180

давление: МПа

п - среднее значение; □- среднее значение ± стандартная ошибка; I среднее значение ±2* стандартное отклонение; ° - выбросы. Рисунок 2 Распределение жировых шариков по размерам, для образцов в ходе процесса гомогенизации.

Предварительная экспериментальная проверка сделанных предположений осуществлялась на образцах в разной степенью гомогенизации, которая достигалась за счет изменения давления в гомогенизаторе (4, 8, 12 МПа). План проведения исследований представлен на рисунке 1. Распределение жировых фракций представлено на рисунке 2. Для определения частотного распределения модуля потерь использовалось «прозвучивание» образцов молока с использованием сигнала «белый шум» (схема экспериментальной установки на ри-

сунке 3) и расчет релаксационного спектра в соответствии с (10). Параллельно проводилось 4 измерения при температуре 18 0С.

Рисунок 3 Экспериментальная установка проведения исследований

Для формализации связи между моментами распределения жировых шариков по диаметру и релаксационного спектра использована зависимость вида: Р1

/ (х) =

/ (х) =

- для первого момента;

х + q1 р1- х2 + р2 ■ х + р3

- для второго момента,

х + q1

параметры которых представлены в таблице.

На рисунке 4 и 5 представлена связь первых и вторых (соответственно) моментов распределений жировых шариков по диаметру и релаксационного спектра.

Т а б л и ц а 1

Параметры математической модели для описания моментов

Параметр модели М1 М2

Р1 0.001533 4.986

р2 - -2.365

р3 - 0.2806

д! -18.15 -0.225

СКО 0.00063 0.000056

Коэф. парной корреляции 0.9924 0.9936

Скорректированный коэф. парной корреляции 0.9908 0.9872

Среднеквадратичная ошибка 0.01126 0.0043

X

Рисунок 4. Связь первых моментов распределений жировых шариков по диаметру и релаксационного спектра

• y1 vs. х1

-untitled fit n

\

0.23 0.232 0.234 0.236 0.238 х1

Рисунок 5. Связь вторых моментов распределений жировых шариков по диаметру и релаксационного спектра

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

ЛИТЕРАТУРА

диаметру и релаксационных спектра (0,98 для первого момента и 0,9 для второго). Кроме того, ввиду зависимости времени релаксации от температуры представляется целесообразным для повышения метрологических характеристик предложенного метода провести исследования при различных значениях температур.

1 Крусь Г.Н. Технология молока и молочных продуктов. М.: Колосс, 2003. 315 с.

2 Kress-Rogers E., Christopher J. B. Instrumentation and sensors for the food industry. Abing-ton: Woodhead Publishing Limited, 2000.

3 Научно-производственное предприятие «Биомер» [Электронный ресурс]. Режим доступа: http://www.biomer.ru (11 марта 2015 г.)

4 Бражников Н.И., Шумиловский Н.Н. Ультразвуковые методы. М.-Л.: Энергия, 1965. 248 с.

5 Riebel U., Loffler F. The Fundamentals of particle size analysis by means of Ultrasonic Spectrometry. Part. 1. // Syst. Charact. 1989. № 6. P. 135-143.

6 Mc Clements D.J. Principles of ultrasonic droplet size determination // Langmuir. 1996. .№12. P. 3454-3461.

7 Михайлов И.Г., Соловьев В.А., Сырников Ю.П. Основы молекулярной акустики. М.: Наука, 1964. 516 с.

8 Ферри Дж. Вязкоупругие свойства полимеров. М.: Изд-во иностр. лит-ры, 1963. 536 с.

9 Honerkamp J., Weese J. Nonlinear regu-larization method for the calculation of relaxation spectra // Rheologica acta. 1993. №32. P. 65-73.

10 Honerkamp J. et al. Determination of the relaxation time spectrum from dynamic moduli using an edge preserving regularization method // Rheologica acta. 2000. № 39. P. 163-173.

11 Honerkamp J., Weese J. Tikhonovs regularization method for ill-posed problems Continuum // Mechanics and Thermodynamics. 1990. V. 2. № 1. P. 17-30.

12 Wasserman H.S., Graessley W.W.J. Rheologica acta. 1992. V. 36. P. 543.

13 Малкин А.Я., Исаев А.И. Реология: концепции, методы, приложения. СПб.: Профессия, 2007. 560 с.

14 Roths T. // Rheologica acta. 2000. V. 39. P. 163.

15 Протопопов И.И. Автоматизация процессов производства сыра. М.: Легкая и пищевая промышленность, 1982. 136 c.

16 Ekomilk BOND ultrasonic milk analyzer // Eon trading inc. [Электронный ресурс]. Режим доступа: http://www.eonbg.com/ (11 марта 2015 г.)

17 Ultrasonic milk analyzers // Milkotronic company. [Электронный ресурс]. Режим доступа: http://www.lactoscan.com/ (11 марта 2015 г.)

18 Хвостов А.А., Ребриков Д.И., Мерзли-кин В.Е. Алгоритм обработки микрофотографий молочных продуктов с использованием ImageJ // Научный вестник ВГАСУ. 2014. №2(4). C. 72-76.

19 Битюков В.К. Автоматизация обработки микрофотографий молочных продуктов с использованием ImageJ и Statistica // Вестник ВГУИТ. 2014. №4(62). С. 58-63.

BeemnvKjBTy^T, №1, 205_

REFERENCES

1 Krus' G.N., Khramtsov A.G., Volokitina Z.V., Karpychev S.V. Tehnologiya moloka i molochnykh produktov [Technology of milk and dairy products]. Moscow, Koloss, 2003, 315 p. (In Russ.).

2 Kress-Rogers E., Christopher J. B. Instrumentation and sensors for the food industry. Abing-ton: Woodhead Publishing Limited, 2000.

3 Nauchno-proizvodstvennoe predpriyatie «Biomer» [Scientific-production enterprise «Biomer»]. Available at: URL: http://www. http://bio-mer.ru. (In Russ.).

4 Brazhnikov N. I., Shunmilovskii N.N. Ul-trazvukovye metody [Ultrasonic methods]. Moscow-Leningrad: Energiya, 1965. 248 p. (In Russ.)

5 Riebel U., Loffler F. The Fundamentals of particle size analysis by means of Ultrasonic Spectrom-etry. Part. 1. [Syst. Charact], 1989, no. 6, pp. 135-143.

6 Mc Clements D.J. Principles of ultrasonic droplet size determination [Langmuir], 1996, no.12, pp. 3454-3461.

7. Mikhailov I.G., Solov'ev V.A., Syrnikov Yu.P. Osnovy molekulyarnoy akustiki [Fundamentals of molecular acoustics]. Moscow: Nauka, 1964. 516 p. (In Russ.)

8 Ferri Dzh. Vyazkouprugie svoistva po-limerov [Viscoelastic properties of polymers]. Moscow, Izdatel'stvo inostrannoi literatury, 1963. 536 p. (In Russ.)

9 Honerkamp J., Weese J. Nonlinear regu-larization method for the calculation of relaxation spectra [Rheologica acta], 1993, no. 32, pp. 65-73.

10 Honerkamp J. et al. Determination of the relaxation time spectrum from dynamic moduli

using an edge preserving regularization method [Rheologica acta], 2000, no. 39, pp. 163-173.

11 Honerkamp J., Weese J. Tikhonovs regularization method for ill-posed problems Continuum [Mechanics and Thermodynamics], 1990, vol. 2, no. 1, pp. 17-30.

12 Wasserman H.S., Graessley W.W.J. [Rheologica acta], 1992, vol. 36, pp. 543.

13 Malkin A. Ya., Isaev A. I. Reologiya: kontseptsii, metody, prilozheniya [Rheology: concepts, methods, and applications] Saint Petersburg: Professiya, 2007. 560 p. (In Russ.).

14 Roths T. [Rheologica acta], 2000, vol. 39, pp.163.

15 Protopopov I.I. Avtomatizatsiya protsessov proizvodstva syra [Automate the production of cheese]. Moscow: Legkaya i pishchevaya promyshlennost', 1982. 136 p. (In Russ.).

16 Ekomilk BOND ultrasonic milk analyzer [Eon trading inc]. Available at: http://www.eonbg.com/ (Accessed 11 March 2015)

17 Ultrasonic milk analyzers [Milkotronic company]. Available at: http://www.lactoscan.com/ (Accessed 11 March 2015)

18 Khvostov A.A., Rebrikov D.I., Merzlikin V.E. The processing algorithm microphotographs of dairy products using ImageJ. Nauchnyi vestnik VGASU [Scientific Bulletin of VSUACE], 2014, no. 2(4), pp. 72-76. (In Russ.).

19 Bityukov V.K., Khvostov A.A., Rebrikov D.I., Merzlikin V.E. Automated processing of micrographs of dairy products using ImageJ and Statis-tica. Vestnik VGUIT. [Bulletin of VSUET], 2014, no. 4(62), pp. 58-63. (In Russ.).