СПОСОБ ОПРЕДЕЛЕНИЯ СОДЕРЖАНИЯ ТВЕРДОГО ЖИРА ПО ДАННЫМ ЯМР-РЕЛАКСАЦИИ, ПРЯМОЙ МЕТОД
METHOD TO DETERMINE CONTENT OF SOLID FAT ACCORDING TO DATA OF NMR-RELAXATION, DIRECT METHOD

Изобретение относится к области пищевой (масложировой, кондитерской), микробиологической, фармацевтической промышленности, в частности к методам контроля качества сырья и готовой продукции, для быстрого определения процентного содержания твердой фазы в двухфазных жировых системах.
Содержание твердого жира (solid fat content, SFC) - важный параметр, влияющий на органолептические свойства готового продукта, зависящий от компонентного состава и температуры. В некоторых случаях он может служить индикатором замены натурального жира на более дешевый аналог. Производство высококачественной продукции, как правило, включает контроль этого параметра на той или иной стадии приготовления.
Известен способ определения содержания твердого жира ГОСТ Р 53158-2008 "Масла растительные, жиры животные и продукты их переработки. Определение содержания твердого жира методом импульсного ядерно-магнитного резонанса" ("прямой метод", без использования дополнительных манипуляций с образцом). Способ основан на измерении уровней сигнала магнитной индукции через 11 и 70 микросекунд после поляризующего магнитного импульса и расчете содержания твердого жира по формуле 1 с использованием фактора коррекции f.
Формула 1
где SFC - содержание твердого жира в процентах,
1
2
S и S - измеренные уровни сигнала магнитной индукции через 11 мкс и 70 мкс соответственно.
Недостатком данного способа является то, что формула 1 не учитывает сложный характер спада магнитной индукции, и при ее использовании в широком диапазоне значений SFC и температур в вычислениях присутствует неустранимая ошибка. В частности фактор коррекции, подобранный при одной температуре, не подходит при другой температуре, что затрудняет получение зависимостей содержания твердого жира в образце от температуры ("кривых плавления"). В то же время именно кривые плавления в наилучшей степени характеризуют контролируемые свойства образца.
Техническим результатом изобретения является устранение ошибки и повышение точности измерений до уровня, полностью использующего аппаратурные возможности современных ЯМР-анализаторов.
Технический результат достигается путем применения следующего алгоритма: при помощи ЯМР-анализатора, работающего в комплексе с персональным компьютером, получают полную кривую спада магнитной индукции после поляризующего импульса с заданным числом накоплений, которая в виде файла, содержащего пары чисел (время и соответствующее значение амплитуды в каждой точке) записывается на ПК. Затем осуществляют автоматизированный подбор параметров математической модели, описываемой формулой
до наилучшего совпадения с формой полной кривой спада (минимальной суммы средних квадратичных отклонений), и затем рассчитывают содержание твердого жира с использованием подобранных таким образом параметров по формуле
. Кривая, соответствующая математической модели с автоматически подобранными параметрами, показана на фиг.1 (гладкая кривая).
Изменение (спад) магнитной индукции в образце после поляризующего импульса представляет собой зависимость сложной формы. Регистрация начинается после истечения так называемого «мертвого времени» прибора, составляющего около 10 мкс. Современные ЯМР-анализаторы позволяют регистрировать по несколько тысяч точек спадов с минимальным расстоянием между точками 0,5 мкс, которые затем можно обрабатывать на ПК с использованием различных программ.
Теоретически абсолютно точно описать форму спада невозможно, так как вклад в форму вносят не только свойства самого образца (который не вполне однороден), но и неоднородность магнитного поля в данном конкретном приборе, а также многие другие факторы. Тем не менее, во многих случаях форма спада может быть с высокой точностью описана при помощи математических моделей с не слишком большим количеством параметров.
Для автоматизированного анализа двухфазных жировых систем нами была выбрана модель в виде суммы трех функций:
Формула 2:
где S(t) -наблюдаемая интенсивность сигнала;
t - время;
Ж
Г
A
А, А, A - амплитуды жидкой части, функций Гаусса и Абрагама;
2Ж
2Г
2A
2
Т, Т и T - соответствующие времена релаксации T;
2H
T - временной коэффициент для учета неоднородности поля;
b - характеристика осцилляции сигнала.
Ж
2Ж
2H
Ж
2Ж
2H
Ж
Г
2Г
2
2
2
Первый член суммы А×exp(-t/Т)×exp(-t/T) описывает релаксацию жидкой части образца. Для полностью жидких образцов характерны спады, по форме близкие к экспоненте А×exp(-t/Т). Дополнительный множитель exp(-t/T) позволяет учесть влияние на форму спада неоднородности поля [Е. Trezza, A.M. Haiduc, G.J.W. Goudappel and J.P.M. van Duynhoven, Rapid phase-compositional assessment of lipid-based food products by time domain NMR, Magnetic Resonance in Chemistry 2006; 44:1023-1030]. В двухфазном образце амплитуда А соответствует содержанию жидкой компоненты. Твердые аморфные образцы жира имеют форму спада, хорошо описывающуюся функцией Гаусса А×exp(-t/Т).
Форма спада у твердых образцов, имеющих признаки квазикристалличности структуры, часто имеет характерную волну, для описания которой удобно использовать функцию Абрагама [Abragam A. The Principles of Nuclear Magnetism, Oxford: London, 1961]
.
Г
А
Содержанию твердого жира в твердом аморфном и твердом квазикристаллическом состояниях соответствует сумма амплитуд А+А. Содержание твердого жира в процентах определяют по формуле 3:
Согласно предлагаемому способу определение содержания твердого жира производится следующим образом:
- образец помещают в ампулу для ЯМР измерений, а ампулу - в датчик ЯМР-анализатора;
- при помощи ЯМР-анализатора, работающего в комплексе с персональным компьютером, получают полную кривую спада магнитной индукции с заданным числом накоплений, которая в виде файла, содержащего пары чисел (время и соответствующее значение амплитуды в каждой точке) записывается на ПК;
- в параллельный столбец заносятся значения, рассчитанные по математической модели;
- параметры модели при помощи программы "поиск решения" в Microsoft Excel, либо при помощи программы Micorocal Origin, либо любой другой аналогичной по функциям программы - подбираются такими, чтобы сумма квадратов отклонений в каждой точке полной кривой спада и ее математической модели была минимальной, т.е. кривые совпадали по форме в максимальной степени. Подобранные таким образом значения амплитуд можно считать отражающими содержание компонент в образце наиболее точно.
1
2
На фиг.1 показаны полная кривая спада («зашумленная» кривая), на которой отмечены точки измерения S и S по прототипу.
На фиг.2 показан пример экспериментально полученного спада для образца, приготовленного из равных частей триолеина и тристеарина (кривая с небольшим приборным «шумом»), и кривая, соответствующая описанию формы спада при помощи математической модели с автоматически подобранными параметрами (гладкая кривая). Математическая модель точно описывает форму спада, в том числе и в области «волны» в районе 20-30 микросекунд от начала спада.
На фиг.3 показан пример экспериментально полученного спада для образца масла какао и его математическая модель. Математическая модель точно описывает форму спада, в том числе и в области «волны».
На фиг.4 показан пример экспериментально полученного спада для образца пальмового жира и его математическая модель. В данном случае «волна» в области 20-30 микросекунд от начала спада отсутствует. Математическая модель точно описывает форму спада.
Точная математическая модель спада позволяет указать на источник ошибок, возникающих при использовании метода по ГОСТ Р 53158-2008. Формула 1 является довольно грубым упрощением. Подобрать фактор коррекции f таким образом, чтобы формула 1 позволяла получать точные значения SFC в широком диапазоне и при разных температурах в принципе невозможно. Соответственно при смене условий измерения каждый раз требуется повторная калибровка.
Предлагаемый способ лишен этого недостатка, и сам не требует калибровки. В то же время с его помощью можно изготовить калибровочные образцы с известным содержанием твердого жира для любой нужной температуры. Ниже приведены примеры иллюстрирующие применение предлагаемого способа определения содержания твердого жира, в том числе и в стандартных образцах со строго заданным составом, предназначенных для калибровки различных ЯМР приборов.
Пример 1:
Образец, приготовленный из 50% тристеарина (твердого при 40°C) и 50% триолеина (жидкого при 40°C) путем точного взвешивания компонент, смешивания их при температуре выше точки плавления тристеарина, помещают в ампулу для ЯМР измерений и затем охлаждают. Погрешность приготовления составляет не более ±0,3%. Ампулу помещают в термостат, выдерживают в течение 15 минут при температуре 40°C, затем помещают в датчик ЯМР-анализатора «ХРОМАТЭК-ПРОТОН 20М» и проводят измерение содержания твердого жира: получают полную кривую спада магнитной индукции, осуществляют автоматизированный подбор параметров математической модели, описываемой формулой
е
до наилучшего совпадения с формой полной кривой спада, и затем рассчитывают содержание тврдого жира с использованием подобранных таким образом параметров поформуле
.
Экспериментально полученная полная кривая спада («зашумленная» кривая) и ее математическая модель (гладкая кривая) показаны на фиг.2. Рассчитанное содержание твердого жира при 40°C составляет 50±0,25%. Приготовленный таким образом образец может использоваться в качестве стандартного образца с содержанием твердого жира 50% для градуировки и поверки различных ЯМР-анализаторов.
Пример 2:
Образец, приготовленный из 25% тристеарина (твердого при 40°C) и 75% триолеина (жидкого при 40°C) путем точного взвешивания компонент, смешивания их при температуре выше точки плавления тристеарина, помещают в ампулу для ЯМР измерений и затем охлаждают. Погрешность приготовления составляет не более ±0,3%. Ампулу помещают в термостат, выдерживают в течение 15 минут при температуре 40°C, затем помещают в датчик ЯМР-анализатора «ХРОМАТЭК-ПРОТОН 20М» и проводят измерение и расчет содержания твердого жира при как в примере 1.
Рассчитанное содержание твердого жира при 40°C составляет 25±0,2%
Приготовленный таким образом образец может использоваться в качестве стандартного образца с содержанием твердого жира 25% для градуировки и поверки различных ЯМР-анализаторов.
Пример 3:
Образец масла какао помещают в ампулу для ЯМР измерений, выдерживают в термостате при температуре 25°C в течение 15 минут, затем быстро перемещают в датчик ЯМР-анализатора и проводят измерение и расчет содержания твердого жира при как в примере 1.
Экспериментально полученная полная кривая спада показана на фиг.3 («зашумленная» кривая) вместе с ее математической моделью (гладкая кривая). Рассчитанное содержание твердого жира в масле какао при 25°C составляет 65,3±0,3%.
Пример 4:
Образец пальмового жира помещают в ампулу для ЯМР измерений, выдерживают в термостате при температуре 20°C в течение 15 минут, затем быстро перемещают в датчик ЯМР-анализатора и проводят измерение и расчет содержания твердого жира при как в примере 1.
Экспериментально полученная полная кривая спада («зашумленная» кривая) и ее математическая модель (гладкая кривая) показаны на фиг.4. Рассчитанное содержание твердого жира в пальмовом жире при 20°C составляет 23,8±0,3%.