Spektroskopia w bliskiej podczerwieni (NIR) zyskała powszechną akceptację jako szybka, wydajna, nieniszcząca i bezemisyjna nowoczesna technologia analityczna; metoda ta była stosowana w ciągu ostatnich kilku lat do wykrywania i analizy jakości produktów oraz pomiaru składników chemicznych w procesach przemysłowych. Najistotniejszą cechą tej metody jest zdolność do rejestrowania widm próbek stałych i ciekłych bez żadnego przygotowania wstępnego, co sprawia, że NIRS jest szczególnie odpowiednia do bezpośredniego i szybkiego wykrywania oraz analizy produktów naturalnych i syntetycznych1,2. W przeciwieństwie do tradycyjnych czujników, które mierzą zmienne procesowe (np. temperaturę, ciśnienie, poziom cieczy itp.) w skali makroskopowej i nieuchronnie ulegają zakłóceniom zewnętrznym i tłu, NIRS wykrywa informacje strukturalne o składzie chemicznym w skali mikroskopowej i molekularnej. Dzięki temu kluczowe informacje mogą być mierzone dokładniej i skuteczniej niż za pomocą innych metod3,4.
Eter polifenylu, jako jeden z tworzyw inżynieryjnych, jest szeroko stosowany ze względu na swoją odporność termiczną, właściwości trudnopalne, izolacyjność, właściwości elektryczne, stabilność wymiarową, odporność na udary, odporność na pełzanie, wytrzymałość mechaniczną i inne właściwości5. Co ważniejsze, jest on nietoksyczny i nieszkodliwy w porównaniu z innymi tworzywami inżynieryjnymi. Obecnie 2,6-ksylenol jest jednym z podstawowych surowców do syntezy eteru polifenylenu, a zazwyczaj jest otrzymywany metodą katalitycznego alkilowania fenolu metanolem6. Głównymi produktami tej metody otrzymywania są o-krezol i 2,6-ksylenol. Po serii etapów separacji i ekstrakcji 2,6-ksylenol jest wykorzystywany do produkcji eteru polifenylenu. Jednak w 2,6-ksylenolu pozostają śladowe ilości o-krezolu. O-krezol nie bierze udziału w syntezie eteru polifenylenu i pozostaje w produkcie końcowym, co prowadzi do obniżenia jakości produktu lub nawet do niespełnienia norm jakościowych. Obecnie większość firm nadal analizuje skład złożonych mieszanin organicznych, takich jak produkty eteru polifenylu w fazie ciekłej zawierające zanieczyszczenia (np. o-krezol), za pomocą fizycznej lub chemicznej analizy separacyjnej, takiej jak chromatografia7,8. Zasadą separacji w chromatografii jest wykorzystanie niewielkich różnic w rozpuszczalności, analizie, adsorpcji, desorpcji lub innym powinowactwie składników mieszaniny względem fazy stacjonarnej i fazy ruchomej. Gdy obie fazy poruszają się względem siebie, składniki są wielokrotnie rozdzielane w obu fazach poprzez opisane procesy. W zależności od badanego obiektu, operacja separacji złożonego materiału trwa zazwyczaj od kilku do kilkudziesięciu minut. Można zatem zauważyć, że wydajność pomiaru jest niska.
Obecnie pomiar jakości produktu oraz zaawansowana technologia sterowania oparta na tej analizie dla nowoczesnego przemysłu chemicznego materiałów precyzyjnych stanowią kluczowy kierunek dalszej poprawy jakości produktów. W przemyśle procesowym produkcji polieterów polifenylowych pomiar w czasie rzeczywistym zawartości o-krezolu w produkcie polieteru polifenylenu ma ogromne znaczenie rozwojowe. Analiza chromatograficzna ewidentnie nie spełnia wymagań zaawansowanej technologii sterowania w zakresie pomiaru substancji i sprzężenia zwrotnego sygnału w czasie rzeczywistym. Dlatego proponujemy metodę regresji cząstkowych najmniejszych kwadratów (PLSR) w celu ustanowienia modelu liniowego między danymi NIRS a stężeniem o-krezolu, co pozwala na pomiar online zawartości o-krezolu w ciekłym produkcie polieteru polifenylenu na wylocie.
Wstępne przetwarzanie danych w spektroskopii NIRS odgrywa kluczową rolę przed zastosowaniem wielowymiarowego modelowania statystycznego. Liczby falowe NIRS w widmie bliskiej podczerwieni oraz rozmiary cząsteczek próbek biologicznych są porównywalne, co prowadzi do występowania nieoczekiwanych efektów rozproszenia wpływających na zarejestrowane widma próbek. Dzięki zastosowaniu odpowiednich metod wstępnego przetwarzania efekty te można w znacznym stopniu wyeliminować9. Najczęściej stosowane techniki wstępnego przetwarzania w NIRS klasyfikuje się jako metody korekcji rozproszenia oraz metody pochodnych widmowych. Pierwsza grupa metod obejmuje multiplikatywną korekcję rozproszenia (MSC), detrending, transformację standard normal variate (SNV) oraz normalizację. Metody pochodnych widmowych obejmują zastosowanie pierwszej i drugiej pochodnej.
Przed opracowaniem ilościowego modelu regresji ważne jest usunięcie niesystematycznych wariacji rozproszonych z danych NIRS, ponieważ mają one istotny wpływ na dokładność modelu predykcyjnego, jego złożoność oraz oszczędność. Wybór odpowiedniej metody wstępnego przetwarzania powinien zawsze zależeć od kolejnego etapu modelowania. W tym przypadku, jeśli zbiór danych widmowych NIR nie spełnia prawa Lamberta-Beera, inne czynniki mają tendencję do kompensowania nieidealnego zachowania predykcji dla przewidywanych komponentów. Wadą istnienia takich zbędnych czynników jest zwiększenie złożoności modelu, a najprawdopodobniej również zmniejszenie jego odporności. Zatem zastosowanie pochodnych widma oraz konwencjonalnej normalizacji danych widmowych stanowi niezbędną część tej metody.
Po wstępnym przetwarzaniu widm uzyskuje się dane NIRS o wysokim stosunku sygnału do szumu i niskich zakłóceniach tła. Nowoczesna analiza NIRS umożliwia szybką akwizycję dużych ilości danych o absorbancji w odpowiednim zakresie widmowym. Skład chemiczny próbki jest następnie przewidywany poprzez wyodrębnienie istotnych zmiennych z informacji zawartych w krzywej widmowej. Zazwyczaj NIRS łączy się z technikami analizy wielowymiarowej w celu przeprowadzenia analiz jakościowych lub ilościowych10. Wielowymiarowa regresja liniowa (MLR) jest powszechnie stosowana do opracowywania i badania zależności matematycznych między danymi a składnikami w procesach przemysłowych i znajduje szerokie zastosowanie w analizie NIRS.
Istnieją jednak dwa fundamentalne problemy podczas wdrażania regresji wieloliniowej (MLR) dla wstępnie przetworzonych danych NIRS. Pierwszym problemem jest redundancja zmiennych. Wysoka wymiarowość danych NIRS często sprawia, że prognozowanie zmiennej zależnej staje się niewiarygodne, ponieważ do modelu włączane są zmienne, które nie wykazują korelacji z komponentami. Te nadmiarowe zmienne obniżają efektywność informacyjną danych widmowych i wpływają na dokładność modelu. Aby wyeliminować redundancję zmiennych, niezbędne jest opracowanie i zmaksymalizowanie korelacji między danymi NIRS a prognozowanymi komponentami.
Kolejnym problemem jest kwestia współliniowości w danych NIRS. Jednym z ważnych założeń modeli regresji liniowej wielorakiej jest brak liniowej zależności między zmiennymi objaśniającymi modelu regresji. Jeśli taka zależność liniowa występuje, oznacza to obecność współliniowości w modelu regresji liniowej i naruszenie tego założenia. W regresjach liniowych wielorakich, takich jak regresja metodą najmniejszych kwadratów (OLSR), wielokrotne korelacje między zmiennymi wpływają na szacowanie parametrów, zwiększają błąd modelu oraz wpływają na jego stabilność. Aby wyeliminować korelację wieloliniową pomiędzy danymi widmowymi NIR, stosujemy metody selekcji zmiennych, które maksymalizują inherentną zmienność próbek.
W niniejszej pracy proponujemy zastosowanie metody PLSR, która stanowi uogólnienie wielokrotnej regresji liniowej i jest szeroko wykorzystywana w dziedzinie NIRS11,12. Metoda PLSR integruje podstawowe funkcje MLR, analizy korelacji kanonicznej (CCA) oraz analizy głównych składowych (PCA), łącząc analizę prognostyczną z analizą konotacyjną danych bezmodelowych. PLSR można podzielić na dwie części. W pierwszej części, za pomocą analizy częściowych najmniejszych kwadratów (PLS), wybierane są składowe zmiennych charakterystycznych oraz składowe przewidywane. PLS maksymalizuje wewnętrzną zmienność głównych składowych poprzez dążenie do uzyskania jak największej kowariancji między głównymi składowymi a składowymi przewidywanymi podczas ekstrakcji głównych składowych. Następnie dla wybranych głównych składowych tworzony jest model OLSR dla stężenia o-kresolu. PLSR nadaje się do analizy danych zaszumionych z licznymi zmiennymi niezależnymi, które wykazują silną kolinearność i wysoką korelację, a także do jednoczesnego modelowania kilku zmiennych odpowiedzi. Ponadto PLSR pozwala na ekstrakcję istotnych informacji ze widm próbek, rozwiązuje problem współliniowości oraz charakteryzuje się wysoką stabilnością i dużą dokładnością przewidywań13,14.
Poniższy protokół opisuje proces wykorzystania modelu PLSR do pomiaru stężenia o-krezolu z wykorzystaniem danych spektralnych NIR. Niezawodność i dokładność modelu są oceniane ilościowo przy użyciu współczynnika determinacji (
), współczynnika korelacji predykcji (
) oraz średnio kwadratowego błędu predykcji walidacji krzyżowej (MSPECV). Ponadto, aby intuicyjnie przedstawić zalety metody PLSR, wskaźniki oceny zostały zwizualizowane na kilku wykresach w ramach analizy jakościowej. Na koniec, wskaźniki oceny eksperymentu zostają przedstawione w formie tabeli, aby ilościowo zilustrować niezawodność i precyzję modelu PLSR.