Profilowanie konopii na podstawie składu pierwiastkowego Część II walidacja metody
Marzena Kuras Marek Wachowicz Profilowanie konopi na podstawie skÅ‚adu pierwiastkowego cz. II (walidacja metody) Na Å‚amach poprzedniego numeru analitycznÄ…. Należą do nich: dokÅ‚ad- kich wymienionych parametrów. Na Problemów Kryminalistyki [1] zosta- noSć, precyzja, selektywnoSć, specy- poczÄ…tku konieczne jest zastanowie- Å‚y opisane metody oznaczania pier- ficznoSć, zakres prostoliniowoSci nie siÄ™, które parametry najlepiej opi- wiastków technikÄ… ICP-OES ze krzywych kalibracji oraz zakres robo- sujÄ… metodÄ™. Oczywiste jest, że np. szczególnym uwzglÄ™dnieniem korek- czy, granica wykrywalnoSci i ozna- przy walidacji metody oznaczania cji efektów matrycowych. We wcze- czalnoSci. Należy podkreSlić, że nie wapnia w próbkach roSlinnych meto- Sniejszych publikacjach omówiono jest niezbÄ™dne wyznaczanie wszyst- dÄ… ICP-OES nie jest potrzebne wy- podstawy metody oraz badania znaczanie granicy wstÄ™pne [2, 3, 4]. W niniejszym wykrywal noSci , Tabela 1 opracowaniu zaS autorzy przedsta- gdyż zawartoSć te- Warunki operacyjne pracy spektrometru wiÄ… zagadnienia zwiÄ…zane z wali- go pierwiastka w ro- Spectrometer set-up parameters dacjÄ… prezentowanej wczeSniej Slinach wynosi kilka metody takie jak: dokÅ‚adnoSć i pre- procent w przelicze- Warunki W a r u n k i Parametr P a r a m e t r odporne o d p o r n e cyzja metody, granica wykrywalno- niu na masÄ™ próbki. Sci, granica oznaczalnoSci, zakres Ważne bÄ™dzie wy- PrzepÅ‚yw gazu plazmowego [l/min] 15 P r z e p Å‚ y w g a z u p l a z m o w e g o [ l / m i n ] prostoliniowy i roboczy krzywych znaczenie zakresu PrzepÅ‚yw gazu pomocniczego [l/min] 0,5 P r z e p Å‚ y w g a z u p o m o c n i c z e g o [ l / m i n ] kalibracyjnych, selektywnoSć me- prostoliniowoSci P r z e p Å‚ y w g a z u p r z e z r o z p y l a c z [ l / m i n ] tody, szacowanie niepewnoSci me- PrzepÅ‚yw gazu przez rozpylacz [l/min] 0,5 krzywych kalibracji. tody. Należy zatem za- Moc plazmy [W] 1450 M o c p l a z m y [ W ] Dla przypomnienia należy po- planować proces WysokoSć obserwacji plazmy [mm] 15 W y s o k o S ć o b s e r w a c j i p l a z m y [ m m ] dać, że w badaniach wykorzysty- walidacji tak, by jak wany jest spektrometr ICP-OES PrzepÅ‚yw próbki [ml/min]* 1,5 najmniejsza liczba P r z e p Å‚ y w p r ó b k i [ m l / m i n ] * Optima 3100XL firmy Perkin El- wyznaczonych pa- Czas opóxnienia [s] 60 C z a s o p ó x n i e n i a [ s ] mer. Na podstawie przeprowadzo- rametrów caÅ‚kowi- nych dotychczas badaÅ„ dalsze cie charakteryzo- * W metodzie wzorca wewnÄ™trznego przepÅ‚yw próbki wynosi 0,65 ml/min, a czas opóxnienia 90 s prace prowadzone bÄ™dÄ… w warun- waÅ‚a stworzonÄ… kach odpornych (ang. robust) pra- cy spektrometru, których parame- odporność try zamieszczono w tabeli 1. specyficzność Próbki ziela konopi przygotowy- liniowość wane sÄ… do badaÅ„ na drodze minera- granica lizacji na mokro z wykorzystaniem oznaczalnoÅ›ci energii mikrofalowej w ukÅ‚adzie za- granica wykrywalnoÅ›ci mkniÄ™tym za pomocÄ… systemu Multi- odtwarzalność wave firmy Anton Paar (Perkin El- powtarzalność mer) z użyciem mieszaniny kwasu dokÅ‚adność azotowego i wody utlenionej. Metoda JakoÅ›ciowa (identyfikacja, potwie rdzenie) Walidacja metody IloÅ›ciowa - wysoki poziom stężeÅ„ IloÅ›ciowa - niski poziom stężeÅ„ Każda nowo opracowana metoda - wymagane zawsze - nie ma potrzeby wyznaczać analityczna powinna być poddana walidacji. Celem walidacji jest wyzna- czenie i dokumentacja parametrów Ryc. 1. Elementy walidacji wymagane przy okreSlonych metodach analitycznych [5] charakteryzujÄ…cych danÄ… metodÄ™ Fig. 1. Elements of validation required in certain analytical methods [5] PROBLEMY KRYMINALISTYKI 253/06 15 metodÄ™. Parametry niezbÄ™dne do wy- Najpowszechniej stosowanÄ… jest Swiadczalnego z wartoSciÄ… referen- znaczenia w procesie walidacji meto- analiza materiałów referencyjnych. cyjnÄ…. dy jakoSciowej i iloSciowej zamiesz- SÄ… to próbki o znanym stężeniu Analiza kilku materiałów certyfiko- czono na rycinie 1. skÅ‚adników w niej zawartych. Ich wanych o różnej zawartoSci skÅ‚adni- Istnieje wiele dokumentów, publi- analiza pozwala na okreSlenie, jaki ka oznaczanego pozwala na oszaco- kacji i przewodników dotyczÄ…cych jest rozrzut wyników miÄ™dzy warto- wanie dokÅ‚adnoSci metody na róż- walidacji i szacowania niepewnoSci Sciami certyfikowanymi a wyznaczo- nym poziomie stężeÅ„. Realizuje siÄ™ metod analitycznych, a najbardziej nymi doSwiadczalnie. Jedynym man- to poprzez wykreSlenie zależnoSci znaczÄ…ce zostaÅ‚y opublikowane kamentem jest fakt, że dla wielu pró- wartoSci certyfikowanej od wyzna- przez AOAC International, Internatio- bek nie sÄ… dostÄ™pne materiaÅ‚y refe- czonej doSwiadczalnie. Otrzymany nal Conference and Harmonization rencyjne, np. próbki konopi z certyfi- wykres powinien przedstawiać liniÄ™ (ICH) i Eurachem [6, 7, 8, 9, 10, 11]. kowanymi zawartoSciami pierwiast- prostÄ…. IdealnÄ… dokÅ‚adnoSć metody Ze wzglÄ™du na znacznÄ… iloSć publika- ków. DokÅ‚adnoSć zwykle wyraża siÄ™ w caÅ‚ym badanym zakresie stężeÅ„ cji zwiÄ…zanych z zagadnieniem wali- jako procentowy odzysk (%Rc): uzyskuje siÄ™ przy współczynniku ko- dacji, jej parametry zostanÄ… opisane relacji równym 1. w skrócie. Innym sposobem wyznaczania do- kÅ‚adnoSci jest porównanie wartoSci X *100 exp DokÅ‚adnoSć (1), wyznaczonych nowo opracowanÄ… me- %Rc = X todÄ… analitycznÄ… i metodÄ… sprawdzo- cert Celem pracy każdego analityka nÄ…, która jest dokÅ‚adna, oraz przez po- jest oznaczenie prawdziwego stęże- równanie miÄ™dzylaboratoryjne. nia substancji. W rzeczywistoSci jed- gdzie: Jak już zostaÅ‚o to podkreSlone nak możliwe jest tylko pewne przybli- Xexp wartoSć wyznaczona doSwiad- w niniejszym artykule, na rynku nie żenie tej wartoSci. MetodÄ™ analitycz- czalnie, sÄ… dostÄ™pne próbki konopi z certyfi- nÄ… uznaje siÄ™ za dokÅ‚adnÄ…, gdy war- Xcert wartoSć certyfikowana. kowanymi zawartoSciami pierwiast- toSć mierzona jest równa wartoSci ków. Dlatego też do okreSlenia do- prawdziwej. Istnieje wiele sposobów WartoSć %Rc równa 100% Swiad- kÅ‚adnoSci metody wybrano cztery ro- na sprawdzenie dokÅ‚adnoSci metody. czy o idealnej zgodnoSci wyniku do- Slinne materiaÅ‚y certyfikowane. Ze Tabela 2 DokÅ‚adnoSć metody przedstawiona na podstawie wyników analizy materiałów certyfikowanych Method accuracy presented on the basis of results of certified material analysis INCT TL1 INCT MPH2 CTA VTL2 CTA OTL1 I N C T T L 1 I N C T M P H 2 C T A V T L 2 C T A O T L 1 ZawartoSć Z a w a r t o S ć pierwiastka p i e r w i a s t k a WartoSć WartoSć WartoSć WartoSć W a r t o S ć W a r t o S ć W a r t o S ć W a r t o S ć Certyfikat Certyfikat Certyfikat Certyfikat C e r t y f i k a t C e r t y f i k a t C e r t y f i k a t C e r t y f i k a t otrzymana otrzymana otrzymana otrzymana o t r z y m a n a o t r z y m a n a o t r z y m a n a o t r z y m a n a B [µg/g] 26* 29 B [ µ g / g ] Ba [µg/g] 43,2 Ä… 3,9 42,7 Ä… 1,5 32,5 Ä… 2,5 31,8 Ä… 0,8 42,7 Ä… 6,6 38,3 Ä… 1,6 84,2 Ä… 11,5 80,9 Ä… 0,8 B a [ µ g / g ] 0,582 Ä… 0,578 Ä… 1,08 Ä… 3,58 Ä… 3,25 Ä… Ca [%] 1,08 Ä… 0,07 3,60 Ä… 0,15 3,17 Ä… 0,12 C a [ % ] 0,052 0,034 0,04 0,15 0,05 7,70 Ä… Cu [µg/g] 20,4 Ä… 1,5 21,1 Ä… 0,9 7,77 Ä… 0,53 18,2 Ä… 0,9 17,2 Ä… 0,7 14,1 Ä… 0,5 13,8 Ä… 0,2 C u [ µ g / g ] 0,35 Fe [µg/g] 432* 448 Ä… 15 460* 459 Ä… 19 1083 Ä… 33 1075 Ä… 51 989* 1004 Ä… 24 F e [ µ g / g ] 0,224 Ä… 0,220 Ä… 0,292 Ä… 0,296 Ä… 0,510 Ä… 0,519 Ä… 0,447 Ä… 0,439 Ä… Mg [%] M g [ % ] 0,017 0,009 0,018 0,010 0,023 0,014 0,021 0,064 1570 Ä… Mn [µg/g] 1559 Ä… 54 191 Ä… 12 195 Ä… 3 79,7 Ä… 2,6 78,1 Ä… 2,5 412 Ä… 14 398 Ä… 7 M n [ µ g / g ] 110 Sr [µg/g] 20,8 Ä… 1,7 22,0 Ä… 1,2 37,6 Ä… 1,1 38,2 Ä… 1,1 110 Ä… 12 119 Ä… 5 201 Ä… 20 184 Ä… 4 S r [ µ g / g ] Zn [µg/g] 34,7 Ä… 2,7 32,2 Ä… 1,4 33,5 Ä… 2,1 32,7 Ä… 1,4 43,3 Ä… 2,1 42,1 Ä… 1,9 49,9 Ä… 2,4 47,8 Ä… 0,8 Z n [ µ g / g ] * WartoSć informacyjna 16 PROBLEMY KRYMINALISTYKI 253/06 wzglÄ™du na podobnÄ… matrycÄ™ próbek operacyjnych, to samo laboratorium Na powtarzalnoSć metody skÅ‚ada zaÅ‚ożono, że próbki te bÄ™dÄ… dobrze oraz powtarzanie badaÅ„ w krótkim siÄ™ wiele elementów. Można wyzna- odzwierciedlać dokÅ‚adnoSć opraco- odstÄ™pie czasu. czać np.: wywanej metody. OdtwarzalnoSć Sr można powtarzalnoSć aparatu odczy- DokÅ‚adnoSć metody wyznaczono oszacować w jednym laboratorium tu wyników podczas analizy po- na podstawie wyników analiz mate- lub poprzez badania miÄ™dzylaborato- jedynczej próbki, riałów certyfikowanych w warunkach ryjne. Jest to precyzja oznaczeÅ„ da- powtarzalnoSć obejmujÄ…cÄ… pre- odtwarzalnych. Wyniki przedstawio- nego skÅ‚adnika w zmiennych warun- cyzjÄ™ mineralizacji próbki oraz ne w tabeli 2 to Srednie z oSmiu po- kach. Można zmieniać nastÄ™pujÄ…ce aparatu podczas jej analizy, miarów jednej próbki. warunki: metodÄ™ pomiarowÄ…, zasadÄ™ powtarzalnoSć caÅ‚ej metody od MetodÄ™ należy uznać za dokÅ‚ad- pomiaru, analityka, aparat, odczynni- pobrania próbki poprzez homo- nÄ…, gdyż wyniki analiz czterech mate- ki, laboratorium i czas pomiaru. genizacjÄ™, mineralizacjÄ™ do ana- riałów certyfikowanych nie różniÄ… siÄ™ PrecyzjÄ™ metody zwykle wyraża lizy. istotnie od wartoSci referencyjnych. siÄ™ przez odchylenie standardowe Precyzja analizy w nowoczesnych Należy zwrócić uwagÄ™, że zawartoSci (SD) lub wzglÄ™dne procentowe od- aparatach analitycznych jest bardzo (stężenia) pierwiastków w poszcze- chylenie standardowe (%RSD), na- duża i zwykle nie ma znaczÄ…cego gólnych materiaÅ‚ach zmieniajÄ… siÄ™ zywane coraz czÄ™Sciej współczynni- udziaÅ‚u w powtarzalnoSci caÅ‚ej meto- w szerokim zakresie, np. dla wapnia kiem zmiennoSci (ang. coefficient of dy. PowtarzalnoSć analizy metodÄ… jest to zakres od 0,582% do 3,60%, variation) o skrócie CV. Wzory po- ICP-OES, którÄ… zamieszczono w ta- a w przypadku strontu od 20,8 mg/g zwalajÄ…ce wyznaczyć SD i CV wyglÄ…- beli 3, okreSlono na podstawie 8 po- do 201 mg/g, zatem opracowana me- dajÄ… nastÄ™pujÄ…co: wtórzeÅ„ analizy próbki CTA VTL2 toda jest dokÅ‚adna w szerokim zakre- (materiaÅ‚ certyfikowany liScie tyto- sie stężeÅ„ dla każdego oznaczanego niu z Wirginii). PowtarzalnoSć wyra- pierwiastka. (2) żono za pomocÄ… współczynnika (xi - x )2 " SD = zmiennoSci (CV). n - 1 Precyzja PowtarzalnoSć analizy dla wiÄ™k- szoSci pierwiastków, wyrażona jako SD * 100 Na precyzjÄ™ metody skÅ‚adajÄ… siÄ™ (3), CV, jest mniejsza od 1%. Tylko dla CV = dwa elementy: powtarzalnoSć i od- x manganu i strontu wartoSci CV sÄ… twarzalnoSć wyników pomiarów. wiÄ™ksze i wynoszÄ… odpowiednio PowtarzalnoSć SR jest to pre- 1,1% i 3,8%. NastÄ™pnie wyznaczono cyzja metody wyznaczona w powta- gdzie: CV na podstawie wyników analizy rzalnych warunkach. Obejmuje ona xi wynik analityczny, próbek materiaÅ‚u certyfikowanego tÄ™ samÄ… procedurÄ™ pomiaru, tego sa- x Srednia z wyników analitycznych, CTA VTL2 przygotowanego w powta- mego analityka, ten sam sprzÄ™t po- n liczba wyników analitycznych. rzalnych warunkach. Wyniki przed- miarowy o takich samych warunkach stawiono w tabeli 4. Tabela 4 Tabela 3 PowtarzalnoSć wyrażona jako CV wyznaczona na podsta- PowtarzalnoSć analizy wykonywanej metodÄ… ICP-OES wie analizy materiaÅ‚u certyfikowanego CTA VTL2 Repeatability of analysis performed by ICP-OES Repeatability expressed as CV determined as the result of analysis of CTA VTL2 certified material Åšrednia Åšrednia Pierwiastek SD CV Pierwiastek SD CV zawartość zawartość B [µg/g] 22 0,2 0,7 B [µg/g] 23 1 4,3 Ba [µg/g] 38,3 0,1 0,2 Ba [µg/g] 38,8 1,6 4,2 Ca [%] 3,58 0,03 0,9 Ca [%] 3,60 0,15 4,3 Cu [µg/g] 17,2 0,1 0,5 Cu [µg/g] 17,3 0,7 3,8 Fe [µg/g] 1075 3 0,3 Fe [µg/g] 1100 51 4,6 Mg [%] 0,517 0,01 1,1 Mg [%] 0,512 0,014 2,6 Mn [µg/g] 78,1 0,2 0,3 Mn [µg/g] 78,6 2,5 3,1 Sr [µg/g] 119 5 3,8 Sr [µg/g] 112 5 4,3 Zn [µg/g] 42,1 0,1 0,3 Zn [µg/g] 42,5 1,9 4,5 PROBLEMY KRYMINALISTYKI 253/06 17 Wszystkie wartoSci CV mieszczÄ… przez jej ujednorodnienie i minerali- W celu wyznaczenia odtwarzalno- siÄ™ w granicach 2,6÷4,6%. zacjÄ™ aż do analizy. Wyznaczone Sci wyrażonej za pomocÄ… CV anali- Należy zauważyć, że materiaÅ‚ cer- wartoSci CV na podstawie wyników zie poddano również 8 próbek mate- tyfikowany to próbka idealnie jedno- analizy próbki konopi zamieszczono riaÅ‚u certyfikowanego CTA VTL2. rodna. Z tego wzglÄ™du przy wyzna- w tabeli 5. Próbki byÅ‚y pobierane, mineralizowa- czaniu powtarzalnoSci metody na Jak wynika z danych zamieszczo- ne i analizowane w różnych dniach podstawie jego analizy pomijane sÄ… nych w tabeli 3, dla baru i manganu analitycznych, w tych dniach przygo- takie etapy jak pobranie próbki i jej współczynniki zmiennoSci wyniosÅ‚y towywano również wzorce i wykony- homogenizacja, które, jak można odpowiednio 6,2 i 5,9%. Dla pozosta- wano kalibracjÄ™. NastÄ™pnie spraw- przypuszczać, majÄ… istotny wkÅ‚ad Å‚ych pierwiastków nie przekroczyÅ‚ on dzano kalibracjÄ™ za pomocÄ… próbki w powtarzalnoSć caÅ‚ej metody. Po- 5%. kontrolnej o certyfikowanych zawar- branie próbki reprezentatywnej roSli- Bardzo istotnÄ… skÅ‚adowÄ… powta- toSciach oznaczanych pierwiastków. ny konopi jest niezwykle trudnym eta- rzalnoSci metody jest powtarzalnoSć KalibracjÄ™ akceptowano, gdy odchy- pem. Otóż do analizy trafiajÄ… zwykle procesu kalibracji. By jÄ… wyznaczyć, lenie od stężenia certyfikowanego bÄ…dx caÅ‚e roSliny konopi, bÄ…dx próbki przygotowano trzy niezależne krzywe nie przekraczaÅ‚o 5%. W przeciwnym Tabela 5 Tabela 6 CV wyznaczone na podstawie wyników analizy próbki konopi PowtarzalnoSć procesu kalibracji metody wyrażona jako CV CV determined basing on results of hemp sample analysis Repeatability of process of method calibration expressed as CV Åšrednie Pierwiastek SD CV Åšrednia nachylenie Pierwiastek SD CV zawartość B 77000 2000 2,6 B [µg/g] 77 3 3,4 Ba 152000 1000 0,7 Ba [µg/g] 27,3 1,7 6,2 Ca 95000 2000 2,1 Ca [%] 3,94 0,14 3,5 Cu 153000 2000 1,3 Cu [µg/g] 21,4 0,2 1,0 Ca 97800 900 0,9 Fe [µg/g] 900 43 4,8 Fe 181000 6000 3,3 0,01 Mg [%] 0,584 2 2,1 Mg 382000 5000 1,3 Mn [µg/g] 128 8 5,9 Mn 599000 1000 0,2 Sr [µg/g] 87 4 4,3 Sr 4870 10 0,2 Zn [µg/g] 70,8 2,2 3,1 Zn 131000 4000 3,1 Tabela 7 konopi wstÄ™pnie rozdrobnione. W ce- kalibracyjne. JednÄ… OdtwarzalnoSć wyrażona jako CV wyznaczona lu wykonania analiz porównawczych z wielkoSci charakte- na podstawie analizy materiaÅ‚u certyfikowanego niezbÄ™dne jest opracowanie proce- ryzujÄ…cych krzywÄ… Reproducibility expressed as CV determined as the result dury poboru próbki, która umożliwia- kalibracji jest jej na- of analysis of certified material Å‚aby wiarygodne porównywanie wyni- chylenie. Porównano ków uzyskanych dla różnych próbek. nachylenia trzech Åšrednia Pierwiastek SD CV Zbadanie zawartoSci pierwiastków wyznaczonych krzy- zawartość w poszczególnych czÄ™Sciach roSliny wych. Uzyskane wy- B [µg/g] 20,5 1,2 5,6 uÅ‚atwiÅ‚o opracowanie procedury po- niki zamieszczono Ba [µg/g] 31,8 0,8 2,5 boru próbek konopi. w tabeli 6. W tym celu z próbki konopi pobra- Wyniki zamiesz- Ca [%] 1,08 0,04 3,7 no 8 próbek analitycznych. Poddano czone w tab. 4 wska- Cu [µg/g] 7,70 0,35 4,5 je ujednorodnieniu w oddzielnych po- zujÄ…, że kalibracja to jemnikach do homogenizacji mÅ‚ynka proces charakteryzu- Fe [µg/g] 459 19 4,2 planetarnego. NastÄ™pnie wykonano jÄ…cy siÄ™ dużą powta- Mg [%] 0,296 0,010 3,3 mineralizacjÄ™ i analizowano. Taki pro- rzalnoSciÄ…. Dla wiÄ™k- Mn [µg/g] 195 3 1,3 ces wyznaczania precyzji w warun- szoSci pierwiastków kach powtarzalnych obejmuje wartoSci współczyn- Sr [µg/g] 38,2 1,1 2,8 wszystkie etapy postÄ™powania anali- nika zmiennoSci nie Zn [µg/g] 32,7 1,4 4,2 tycznego od pobrania próbki, po- przekraczaÅ‚y 2%. 18 PROBLEMY KRYMINALISTYKI 253/06 razie przygotowywano nowÄ… krzywÄ… SDb odchylenie standardowe Slepej Kolejnym sposobem wyznaczania kalibracji. Wyniki przedstawiono próby. granicy wykrywalnoSci jest metoda w tabeli 7. W drugim przypadku wyrażenie to wizualna. Polega ona na ocenie wi- Podobnie jak w przypadku powta- przyjmuje postać: zualnej sygnałów uzyskanych dla rzalnoSci, wyznaczono nastÄ™pnie od- różnych stężeÅ„ oznaczanej substan- kc0RSDb twarzalnoSć metody, poddajÄ…c anali- (5), cji. Stężenia te powinny być tak do- GW = zie próbkÄ™ konopi pobranÄ…, homoge- brane, by obejmowaÅ‚y zakres niewy- Ra-b nizowanÄ…, mineralizowanÄ… i analizo- krywany i wykrywany przez aparat. wanÄ… w różnych dniach analitycz- Na podstawie oceny wizualnej wybie- nych. Uzyskane wyniki zamieszczo- gdzie: ramy sygnaÅ‚, który jest wyraxnie wy- no w tabeli 8. k staÅ‚a, Tabela 8 Tabela 9 OdtwarzalnoSć wyznaczona na podstawie analizy WartoSci granicy wykrywalnoSci i oznaczalnoSci metody próbki konopi wyznaczone na podstawie wyników analizy Slepych prób Reproducibility determined basing on analysis Values of detection and determination limits obtained basing of hemp sample on results of analysis of blind samples Åšrednia Pierwiastek SD CV Åšrednie stężenie SD GW GO zawartość Pierwiastek [mg/l] [mg/l] [mg/l] [mg/l] B [µg/g] 26 1 5,7 B 0,05 0,02 0,10 0,21 Ba [µg/g] 43,1 1,3 3,0 Ba 0,01 0,00 0,02 0,03 Ca [%] 1,25 0,05 3,9 Ca 0,82 0,22 1,47 2,99 Cu [µg/g] 10,1 0,4 4,3 Cu 0,01 0,00 0,03 0,06 Fe [µg/g] 512 23 4,5 Fe 0,06 0,04 0,17 0,44 Mg [%] 0,273 0,011 4,0 Mg 0,02 0,01 0,04 0,07 Mn [µg/g] 205 4 2,0 Mn 0,002 0,001 0,006 0,015 Sr [µg/g] 58,1 2,0 3,5 Sr 0,37 0,04 0,49 0,76 Zn [µg/g] 50,8 2,4 4,8 Zn 0,02 0,01 0,03 0,07 Granica wykrywalnoSci c0 stężenie odpowiadajÄ…ce sygnaÅ‚o- odrÄ™bniony z tÅ‚a, a stężenie mu odpo- i oznaczalnoSci wi analitu ya, wiadajÄ…ce to granica wykrywalnoSci. RSDb wzglÄ™dne odchylenie standar- Granica oznaczalnoSci to naj- GranicÄ™ wykrywalnoSci, zgodnie dowe Slepej próby, mniejsze stężenie, jakie można wy- z zaleceniem organizacji IUPAC, de- Ra-b stosunek sygnałów tÅ‚a i analitu kryć danÄ… metodÄ… z odpowiedniÄ… pre- finiuje siÄ™ jako najmniejsze stężenie, o stężeniu c0. cyzjÄ…. By wyznaczyć granicÄ™ ozna- któremu odpowiada sygnaÅ‚ różniÄ…cy czalnoSci (GO), należy we wzorach siÄ™ statystycznie (istotnie) od sygnaÅ‚u W kolejnej metodzie do wyznacza- 4÷6 przyjąć wartoSć staÅ‚ej k=6. Slepej próby. NajczÄ™Sciej stosowane nia granicy wykrywalnoSci stosowany Istnieje wiele metod wyznaczania sÄ… podejScia oparte albo na stosunku jest nastÄ™pujÄ…cy wzór: granicy wykrywalnoSci i oznaczalno- sygnaÅ‚u do szumu, albo na stosunku Sci. W realizowanych badaniach po- sygnaÅ‚u do tÅ‚a i wzglÄ™dnym odchyle- stanowiono przedstawić wartoSci GW = X + kSDb blank niu standardowym tÅ‚a. W pierwszym (6), granicy wykrywalnoSci i oznaczalno- przypadku wyrażenie na granicÄ™ wy- Sci uzyskane różnymi metodami. Ja- krywalnoSci GW można przedstawić ko pierwszÄ… wybrano wyznaczanie jako [12]: gdzie: tych parametrów na podstawie analiz k staÅ‚a, Slepej próby. Wyniki dla 10 Slepych (4), Xblank wartoSć Slepej próby, prób uzyskano w warunkach odtwa- kc0SDb GW = SDb odchylenie standardowe Slepej rzalnych. Na ich podstawie wyzna- ya próby. czono wartoSć granicy wykrywalnoSci (GW) i oznaczalnoSci (GO) wedÅ‚ug gdzie: Należy zaznaczyć, że dla granicy wzorów 4÷6. Uzyskane wyniki przed- k staÅ‚a, wykrywalnoSci wartoSć staÅ‚ej k wyno- stawiono w tabeli 9. c0 stężenie odpowiadajÄ…ce sygnaÅ‚o- si 3. KolejnÄ… metodÄ… wyznaczania GW wi analitu ya, i GO jest metoda wizualna. Jak PROBLEMY KRYMINALISTYKI 253/06 19 wspomniano, polega ona na zareje- ny być dobrane tak, by obejmowaÅ‚y azotowy. Widma pierwiastków o wzra- strowaniu i analizie widm pierwiast- poziom, którego aparat nie jest w sta- stajÄ…cym stężeniu zamieszczono na ków w próbkach wzorców o różnym nie wykryć. W badaniach przygotowa- rycinach 2÷10. Należy zaznaczyć, że stężeniu. Przygotowywana jest seria no serie wzorców pierwiastków o stÄ™- zamieszczono widma tylko dla tych roztworów o wzrastajÄ…cym stężeniu żeniach: 0; 0,01; 0,02; 0,05; 0,1; 0,2; stężeÅ„ pierwiastków, które sÄ… nie- oznaczanego pierwiastka. Należy za- 0,5; 1 i 2 mg/l. MatrycÄ… wszystkich zbÄ™dne do wyodrÄ™bnienia granicy wy- znaczyć, że stężenia wzorców powin- wzorców jest 20-procentowy kwas krywalnoSci i oznaczalnoSci. B Ba 0 mg/l 0 mg/l 18000 47000 0,01 mg/l 17000 0,01 mg/l 44000 16000 0,02 mg/l 41000 15000 0,02 mg/l 0,05 mg/l 38000 14000 35000 13000 0,1 mg/l 0,05 mg/l 12000 32000 0,2 mg/l 11000 29000 0,1 mg/l 10000 0,5 mg/l 26000 233.505 233.515 233.525 233.535 233.545 249.65 249.66 249.67 249.68 249.69 249.7 249.71 DÅ‚ugość fali [nm] DÅ‚ugość fali [nm] Ryc. 2. Widma uzyskane dla boru o różnych stężeniach Ryc. 3. Widma uzyskane dla baru o różnych stężeniach Fig. 2. Spectra of boron in various concentrations Fig. 3. Spectra of barium in various concentrations Cu 0 mg/l Ca 0 mg/l 135000 0,01 mg/l 63000 0,01 mg/l 130000 61000 0,02 mg/l 125000 0,02 mg/l 59000 0,05 mg/l 120000 57000 0,05 mg/l 115000 55000 0,1 mg/l 110000 53000 0,1 mg/l 0,2 mg/l 105000 51000 0,2 mg/l 0,5 mg/l 49000 100000 315.86 315.87 315.88 315.89 315.9 315.91 315.92 327.36 327.38 327.4 327.42 DÅ‚ugość fali [nm] DÅ‚ugość fali [nm] Ryc. 4. Widma uzyskane dla wapnia o różnych stężeniach Ryc. 5. Widma uzyskane dla miedzi o różnych stężeniach Fig. 4. Spectra of calcium in various concentrations Fig. 5. Spectra of copper in various concentrations Fe Mg 0 mg/l 0 mg/l 30000 76000 0,01 mg/l 0,01 mg/l 28000 71000 66000 26000 0,02 mg/l 0,02 mg/l 24000 61000 22000 56000 0,05 mg/l 0,05 mg/l 20000 51000 0,1 mg/l 0,1 mg/l 18000 46000 285.17 285.19 285.21 285.23 285.25 238.18 238.19 238.2 238.21 238.22 238.23 DÅ‚ugość fali [nm] DÅ‚ugość fali [nm] Ryc. 6. Widma uzyskane dla żelaza o różnych stężeniach Ryc. 7. Widma uzyskane dla magnezu o różnych stężeniach Fig. 6. Spectra of iron in various concentrations Fig. 7. Spectra of magnesium in various concentrations 20 PROBLEMY KRYMINALISTYKI 253/06 Intensywność [cps] Intensywność [cps] Intensywność [cps] Intensywność [cps] Intensywność [cps] Intensywność [cps] 0 mg/l Mn Sr 0 mg/l 0,01 mg/l 24500 21500 0,02 mg/l 23500 0,05 mg/l 20500 22500 0,1 mg/l 21500 0,01 mg/l 19500 20500 0,2 mg/l 18500 19500 0,5 mg/l 18500 17500 1 mg/l 17500 0,02 mg/l 16500 16500 2 mg/l 257.57 257.59 257.61 257.63 257.65 232.21 232.22 232.23 232.24 232.25 232.26 DÅ‚ugość fali [nm] DÅ‚ugość fali [nm] Ryc. 8. Widma uzyskane dla manganu o różnych stężeniach Ryc. 9. Widma uzyskane dla strontu o różnych stężeniach Fig. 8. Spectra of manganese in various concentrations Fig. 9. Spectra of strontium in various concentrations Tabela 10 Zn 0 mg/l WartoSci granicy wykrywalnoSci i oznaczalnoSci metody 10500 wyznaczone na podstawie wizualnej oceny widm 10000 Values of detection and determination limits obtained basing 9500 0,01 mg/l on visual assessment of spectra 9000 8500 8000 G W G O 0,02 mg/l GW GO Pierwiastek P i e r w i a s t e k 7500 [mg/l] [mg/l] [ m g / l ] [ m g / l ] 7000 6500 B 0,2 0,5 B 0,05 mg/l 6000 206.176 206.186 206.196 206.206 206.216 Ba 0,05 0,1 B a DÅ‚ugość fali [nm] Ca 0,05 0,2 C a Cu 0,2 0,5 C u Ryc. 10. Widma uzyskane dla cynku o różnych stężeniach Fig. 10. Spectra of zinc in various concentrations Fe 0,05 0,1 F e Mg 0,05 0,1 M g W tabeli 10 podano wartoSci grani- Mn 0,01 0,02 M n cy wykrywalnoSci i oznaczalnoSci wy- Sr 1 2 S r znaczone na podstawie oceny wizu- Zn Z n alnej widm. 0,02 0,05 Zakres prostoliniowy i roboczy W analizie chemicznej oczekiwane n krzywych kalibracyjnych (x - x)(y - y) sÄ… korelacje lepsze od 0,99. " i i i =1 (7), By wyznaczyć zakres prostolinio- r = n n LiniowoSć to niezwykle istotny pa- 2 2 wy metody, należy poddać analizie (x - x) (y - y) "" i i rametr metody analitycznej. OkreSla seriÄ™ co najmniej piÄ™ciu roztworów i =1 i =1 ona zakres prostoliniowej zależnoSci analitu o różnym stężeniu, a nastÄ™p- pomiÄ™dzy stężeniem oznaczanego nie obliczyć współczynnik korelacji. skÅ‚adnika i uzyskanym sygnaÅ‚em gdzie: Zakres roboczy krzywej kalibracji analitycznym. Ze wzglÄ™dów praktycz- xi, yi współrzÄ™dne kolejnych punk- to zakres stężeÅ„, z jakiego korzysta nych pożądane jest, by zakres pro- tów w analizowanym zbiorze, siÄ™ w danej metodzie analitycznej. stoliniowoSci byÅ‚ jak najwiÄ™kszy. n liczba punktów, Musi charakteryzować siÄ™ dobrym Umożliwia to oznaczanie maÅ‚ych x, y odpowiednio wartoSci Srednie współczynnikiem korelacji. i dużych stężeÅ„ skÅ‚adnika bez ko- współrzÄ™dnych. RóżnicÄ™ miÄ™dzy zakresem prosto- niecznoSci rozcieÅ„czania próbki. Ko- liniowym a roboczym metody przed- relacjÄ™ pomiÄ™dzy stężeniem a sygna- ZależnoSć miÄ™dzy y (zmienna za- stawiono na rycinie 11. Å‚em można opisać za pomocÄ… współ- leżna) a x (zmienna niezależna) bÄ™- Metoda ICP-OES ma duży zakres czynnika korelacji r. dzie idealnie liniowa, gdy współczyn- prostoliniowoSci krzywych kalibracyj- nik korelacji wyniesie 1 albo 1. nych, zwykle siÄ™gajÄ…cy kilku rzÄ™dów PROBLEMY KRYMINALISTYKI 253/06 21 Intensywność [cps] Intensywność [cps] Intensywność [cps] nachylenie krzywej kalibracji za- mieszczono wczeSniej w tabeli 6. Tabela 11 Zakres roboczy krzywych kalibracji Working range of calibration curves Zakres roboczy Pierwiastek [mg/l] B 0÷10 Ba 0÷10 Ca 0÷500 Cu 0÷10 Fe 0÷100 Mg 0÷200 Ryc. 11. Zakres prostoliniowy i roboczy metody analitycznej Fig. 11. Linear and working range of analytical method Mn 0÷50 Sr 0÷10 wielkoSci. Jednak każda krzywa kali- Powyżej stężenia 200 mg/l nastÄ™- bracji odznacza siÄ™ ograniczonym puje przeÅ‚adowanie linii Mg 285,213. Zn 0÷10 zakresem. Zatem oznaczanie wyższych stężeÅ„ W przypadku analizy próbek kono- za pomocÄ… tej linii nie jest możliwe. SelektywnoSć pi najwiÄ™ksze stężenia oznaczane sÄ… WartoSć współczynnika (0,9995) dla wapnia i magnezu do odpo- wskazuje na dużą prostoliniowoSć SelektywnoSć to zdolnoSć metody wiednio 2100 mg/l i 1270 mg/l. Jedy- krzywej kalibracji. do odróżniania oznaczanego analitu nie w przypadku tych dwóch pier- W przypadku wapnia krzywa jest od innych substancji. Badanie selek- wiastków jest uzasadnione wyzna- prostoliniowa w zakresie 0÷700 mg/l. tywnoSci zwykle przeprowadza siÄ™ czanie prostoliniowego zakresu krzy- Współczynnik korelacji wynosi poprzez analizÄ™ próbek analitów, do wych kalibracji. W tym celu przygoto- 0,9998. Powyżej 700 mg/l nastÄ™puje których dodano potencjalne interfe- wano i analizowano serie roztworów przeÅ‚adowanie linii wapnia, co unie- renty, i obserwacjÄ™ uzyskiwanych sy- kalibracyjnych zawierajÄ…cych: możliwia oznaczanie wyższych stÄ™- gnałów. dla wapnia 0, 100, 200, 400, żeÅ„. Metoda ICP-OES nie odznacza 600, 700, 800, 900, 1000 mg/l, W tabeli 11 zamieszczono zakresy siÄ™ niestety dobrÄ… selektywnoSciÄ… ani dla magnezu 0, 50, 100, 150, robocze krzywych kalibracji dla specyficznoSciÄ…. MożliwoSć jedno- 175 i 200 mg/l. wszystkich oznaczanych pierwiast- czesnego oznaczania zdecydowanej Uzyskane wyniki zamieszczono ków. CzuÅ‚oSć metody wyrażona jako wiÄ™kszoSci pierwiastków ukÅ‚adu na rycinach 12 i 13. okresowego jest dużą zaletÄ…, lecz Mg 285.213 Ca 315.877 80000000 80000000 70000000 70000000 60000000 60000000 50000000 50000000 40000000 40000000 30000000 30000000 20000000 20000000 2 10000000 10000000 2 R = 0.9995 R = 0.9998 0 0 0 100 200 300 0 100 200 300 400 500 600 700 800 stężenie Mg [mg/l] stężenie Ca [mg/l] Ryc. 12. Zakres prostoliniowy krzywej kalibracji wyznaczonej dla magnezu Ryc. 13. Zakres prostoliniowy krzywej kalibracji wyznaczonej dla wapnia Fig. 12. Rectilinear range of calibration curve for magnesium Fig. 13. Rectilinear range of calibration curve for calcium 22 PROBLEMY KRYMINALISTYKI 253/06 intensywność [cps] intensywność [cps] może również być wadÄ…. WystÄ™po- mocÄ… którego wykonywano analizy, przygotowano seriÄ™ czterech roztwo- wanie sygnałów pierwiastka przy wyposażony jest w detektor w zakre- rów wzorcowych oraz SlepÄ… próbÄ™. wielu różnych dÅ‚ugoSciach fali może sie UV. W tym zakresie stront emituje W każdym roztworze stężenie strontu prowadzić do ich nakÅ‚adania siÄ™ z sy- tylko jednÄ… falÄ™ o stosunkowo maÅ‚ej byÅ‚o staÅ‚e i wynosiÅ‚o 2,5 mg/l. Stęże- gnaÅ‚ami innych pierwiastków. Naj- czuÅ‚oSci. Interferencja spektralna nie żelaza wynosiÅ‚o odpowiednio 0, wiÄ™ksze bÅ‚Ä™dy zwiÄ…zane z nakÅ‚ada- mogÅ‚aby uniemożliwić dokÅ‚adne 10, 20 i 50 mg/l. Należy zaznaczyć, niem siÄ™ sygnałów obserwowane sÄ…, oznaczanie tego pierwiastka w prób- że stężenia strontu i żelaza dobrano gdy obok siebie wystÄ™pujÄ… silne linie kach konopi. Jedynym pierwiastkiem, tak, by odzwierciedlaÅ‚y poziom wy- pierwiastka oznaczanego i interferen- który mógÅ‚by zakłócać sygnaÅ‚ Sr, jest stÄ™pujÄ…cy w próbkach konopi. Wyniki ta. Emitowanie wielu fal daje możli- żelazo. Linia Fe232,233 jest bardzo analizy przedstawiono na rycinach woSć prostej eliminacji interferencji sÅ‚aba, lecz ze wzglÄ™du na idealne na- 14, 15 oraz w tabeli 13. spektralnych poprzez wybranie innej Å‚ożenie z liniÄ… Sr postanowiono nieobarczonej interferencjami linii sprawdzić jej analitycznej. Jednak czÄ™sto takie linie wpÅ‚yw. Dodatko- 24500 Sr sÄ… mniej czuÅ‚e i wykrycie stężeÅ„ Sla- wym czynnikiem dowych za ich pomocÄ… jest utrudnio- przemawiajÄ…cym 22000 ne. za sprawdzeniem 19500 W trakcie walidacji metody ozna- wystÄ™powania tej czania pierwiastków w próbkach ko- interferencji jest 17000 nopi sprawdzono, jakie potencjalne fakt, że stront interferenty emitujÄ… sygnaÅ‚y w pobli- w próbkach kono- 14500 żu linii analitycznych oznaczanych pi wystÄ™puje na 12000 pierwiastków. Zestawienie tych infor- poziomie stężeÅ„ 232.22 232.23 232.24 232.25 232.26 macji zamieszczono w tabeli 12. Bra- kilkakrotnie niż- dÅ‚ugość fali [nm] no pod uwagÄ™ linie interferentów szym od żelaza. Sr 0 mg/l + Fe 0 mg/l Sr 2.5 mg/l znajdujÄ…ce siÄ™ w odlegÅ‚oSci mniejszej Zatem żelazo mo- Sr 2.5 mg/l + Fe 10 mg/l Sr 2.5 mg/l + Fe 20 mg/l niż 0,005 nm od linii oznaczanego że mieć faktyczny Sr 2.5 mg/l + Fe 50 mg/l pierwiastka. wpÅ‚yw na wyniki W wyznaczaniu interferencji spek- oznaczeÅ„. W celu Ryc. 14. Widma uzyskane dla strontu w matrycy żelaza o różnym stężeniu tralnych najwiÄ™kszÄ… uwagÄ™ zwrócono sprawdzenia, czy Fig. 14. Spectra of strontium in ferrous matrix in various concentrations na liniÄ™ strontu. Spektrometr, za po- tak jest w istocie, Tabela 12 Analiza widm zamieszczonych na Potencjalne interferenty [13] rycinie 14 wskazuje, że obecnoSć że- Possible interferents laza powoduje przede wszystkim wzrost poziomu tÅ‚a, który jest propor- DÅ‚ugość fali cjonalny do stężenia żelaza w prób- Pierwiastek Potencjalny interferent/dÅ‚ugość fali [nm] [nm] ce. Na podstawie tych widm trudno okreSlić, czy wystÄ™pujÄ… interferencje Tc 249,677; Sn 249,677; Hg 249,678; Cr 249,681 B 249,677 spektralne, czyli czy nastÄ™puje Ba 233,527 Nb 233,531; V 233,533 wzrost intensywnoSci sygnaÅ‚u strontu w obecnoSci coraz wiÄ™kszych stężeÅ„ Be 234,861 Mo 234,858; Ta 234,859; Zr 234,859 żelaza. Be 313,107 Th 313,107; Zr 313,111 Na rycinie 15 zamieszczono wid- Ca 315,887 Ce 315,888 ma takie same jak na rycinie 14, lecz z przesuniÄ™tym poziomem tÅ‚a. Wy- Nb 327,389; U 327,390; Th 327,392; Ce 327,393; stÄ™powanie tÅ‚a na podobnym pozio- Cu 327,393 Co 327,393; Ta 327,396; Mo 327,396; Ce 327,396; mie uÅ‚atwia interpretacjÄ™ wizualnÄ… Sb 327,397 widm. Okazuje siÄ™, że wysokoSć sy- Fe 238,204 V 238,203; Hg 238,206 gnaÅ‚u Sr roSnie wraz ze wzrostem za- Mg 285,213 U 285,209; W 285,210; Ce 285,212; Mo 285,213; wartoSci żelaza w próbce. Rwiadczy to o wystÄ™powaniu interferencji spek- Mn 257,610 Ru 257,609; Co 257,610; Zr 257,610 tralnych. Sr 232,235 Fe 232,233 Dane zestawione w tabeli 13 po- Zn 206,200 Nb 206,297; In 206,200; V 206,200 zwalajÄ… na oszacowanie wielkoSci tej interferencji. Okazuje siÄ™, że stężenie Y 371,029 Ce 371,025; W 371,029; U 371,031 żelaza w wysokoSci 10 mg/l powodu- PROBLEMY KRYMINALISTYKI 253/06 23 intensywność [cps] Tabela 13 Sr 22500 IntensywnoSci emisji uzyskane dla strontu w obecnoSci 21000 matrycy żelaza o różnym stężeniu 19500 Emission intensities obtained for strontium in presence of ferrous matrix 18000 16500 Stężenie Stężenie Intensywność Intensywność wzglÄ™dna 15000 Sr Fe dla Sr dla Sr 13500 [mg/l] [mg/l] [cps] [%] 12000 2,5 0 10537 100 232.22 232.23 232.24 232.25 232.26 dÅ‚ugość fali [nm] 2,5 10 11330 108 Sr 0 mg/l + Fe 0 mg/l Sr 2.5 mg/l 2,5 20 12037 114 Sr 2.5 mg/l + Fe 10 mg/l Sr 2.5 mg/l + Fe 20 mg/l Sr 2.5 mg/l + Fe 50 mg/l 2,5 50 12805 122 Ryc. 15. Wystandaryzowane widma uzyskane dla strontu w matrycy żelaza o różnym stężeniu Fig. 15. Standardised spectra of strontium in ferrous matrix in various concentrations je wzrost intensywnoSci strontu o 8%, niepewnoSć wzorca wyjSciowe- kÄ™ i wzorce robocze, objÄ™toSć pipeto- a w wysokoSci 50 mg/l aż o 22% go, z którego w wyniku rozcieÅ„czenia wanego wzorca oraz masa próbki. w stosunku do sygnaÅ‚u uzyskanego przygotowywane sÄ… wzorce robocze, ProstokÄ…tnemu rozkÅ‚adowi praw- bez dodatku interferenta. niepewnoSć objÄ™toSci pipetowa- dopodobieÅ„stwa, dla którego niepew- Wyniki uzyskane dla 105 próbek nego wzorca (niepewnoSć pipety), noSć standardowa jest równa konopi wskazujÄ…, że oznaczane stÄ™- niepewnoSć kolb, w których Sd żenie żelaza nie przekracza 40 mg/l. przygotowywane sÄ… wzorce robocze. U= 3 Należy zaznaczyć, że dla 72% pró- SkÅ‚adniki niepewnoSci zwiÄ…zane bek oznaczone stężenie tego pier- z analizÄ… wzorców roboczych metodÄ… podlegajÄ… objÄ™toSci kolb, w któ- wiastka nie przekraczaÅ‚o 10 mg/l, ICP-OES: rych przygotowywano próbkÄ™ i wzor- a tylko w 5% próbek wyniosÅ‚o wiÄ™cej niepewnoSć wyznaczenia inten- ce robocze, a także objÄ™toSć pipeto- niż 20 mg/l. Te informacje wskazujÄ…, sywnoSci sygnaÅ‚u dla pierwiastków wanego wzorca. Natomiast masa że efekt interferencyjny dla strontu we wzorcach roboczych, próbki, czyli niepewnoSć wagi, bÄ™dzie jest istotny dla 28% badanych pró- skÅ‚adniki niepewnoSci zwiÄ…zane podlegać trójkÄ…tnemu rozkÅ‚adowi bek. z analizÄ… próbki konopi, prawdopodobieÅ„stwa, dla którego niepewnoSć wyznaczenia inten- Sd Szacowanie niepewnoSci metody sywnoSci sygnaÅ‚u dla pierwiastków U= 6 w próbce konopi, Na niepewnoSć metody oznacza- niepewnoSć wyznaczenia inten- Aby opracować model szacowa- nia pierwiastków w próbce konopi sywnoSci sygnaÅ‚u dla pierwiastków nia niepewnoSci, należy również technikÄ… ICP-OES wpÅ‚ywa wiele w Slepej próbie. znać model matematyczny, wedÅ‚ug czynników. SkÅ‚adniki niepewnoSci należy po- którego spektrometr oblicza stężenie NiepewnoSć wynikajÄ…ca z procesu grupować w zależnoSci od tego, ja- analitu na podstawie zmierzonej in- przygotowania próbki obejmuje na- kim typom szacowania niepewnoSci tensywnoSci sygnaÅ‚u. stÄ™pujÄ…ce skÅ‚adniki: podlegajÄ… (A czy B). Równanie krzywej wzorcowej niepewnoSć wyznaczania masy W tym przypadku typ A szacowa- można przedstawić nastÄ™pujÄ…co: próbki (niepewnoSć wagi analitycz- nia niepewnoSci obejmuje skÅ‚adniki, nej), które zostaÅ‚y wyznaczone doSwiad- y = ax + b (8), niepewnoSć kolby, w której przy- czalnie, tj. intensywnoSci pierwiast- gotowywana jest próbka. ków dla próbki, intensywnoSci pier- gdzie: NiepewnoSć wyznaczania kalibra- wiastków dla Slepej próby i intensyw- y intensywnoSć sygnaÅ‚u analityczne- cji i oznaczania pierwiastka w próbce noSci pierwiastków dla wzorców. Dla go, konopi obejmuje cztery rodzaje tych skÅ‚adników niepewnoSć standar- x stężenie, skÅ‚adników: dowa (U) jest równa odchyleniu stan- a nachylenie krzywej, skÅ‚adniki niepewnoSci zwiÄ…zane dardowemu. b punkt przeciÄ™cia krzywej z osiÄ… y. z przygotowywaniem roztworów SkÅ‚adniki podlegajÄ…ce typowi B wzorcowych, szacowania niepewnoSci to: objÄ™toSć WartoSć a obliczana jest ze wzoru 9. kolb, w których przygotowywano prób- WartoSć b obliczana jest ze wzoru 10: 24 PROBLEMY KRYMINALISTYKI 253/06 intensywność [cps] (x1 - x)(y1 - y)+ (x2 - x)(y2 - y)+ (x3 - x)(y3 - y)+ (x4 - x)(y4 - y)matycznego (wzór 13). Ponieważ to jeszcze bardzie komplikuje model ma- a = 2 2 2 2 (x1 - x) + (x2 - x) + (x3 - x) + (x4 - x) tematyczny, zrezygnowano z wsta- (9) wiania tego wzoru w tym miejscu pra- cy. PeÅ‚ny model matematyczny jest gdzie: rozcieÅ„czanie z odpowiedniego uwzglÄ™dniony w skoroszycie Microsoft xi wartoSć stężenia wzorca, wzorca wyjSciowego. Stężenie wzor- Excel, który byÅ‚ wykorzystywany w ob- yi wartoSć sygnaÅ‚u analitycznego. liczeniach. y1 + y2 + y3 + y4 (x1 - x)(y1 - y)+ (x2 - x)(y2 - y)+ (x3 - x)(y3 - y)+ (x4 - x)(y4 - y)" x1 + x2 + x3 + x4 b = - 2 2 2 2 4 4 (x1 - x) + (x2 - x) + (x3 - x) + (x4 - x) (10) Stężenie nieznanej próbki w mg/l ca do kalibracji można wyrazić wzo- W tabeli 14 podano wzglÄ™dny obliczane jest z zależnoSci 11. rem 14. udziaÅ‚ poszczególnych elementów w niepewnoSci zÅ‚ożonej. Dla pierwiastków B, Ba, Cu, Sr i Zn (yox - yosp )- b (11), najwiÄ™kszy udziaÅ‚ w niepewnoSci me- xo = (14), tody ma stężenie wzorca głównego. Cwz "Vi a xi = NiepewnoSć tego wzorca stanowi od Vk gdzie: 39,9% caÅ‚kowitej niepewnoSci dla x0 wartoSć stężenia [mg/l], miedzi do 61,0% dla cynku. Otóż jest y0x wartoSć Srednia sygnaÅ‚u anali- gdzie: to zwiÄ…zane z faktem, że roztwór tego tycznego dla nieznanej próbki [mg/l], Cwz stężenie wzorca wyjSciowego wzorca przygotowano poprzez roz- y0sp wartoSć Srednia sygnaÅ‚u anali- [mg/l], cieÅ„czenie wzorca wyjSciowego tycznego dla Slepej próby [mg/l]. Vi objÄ™toSć odpipetowanego wzorca o stężeniu każdego z pierwiastków [ml], 1000 mg/l. Oczywiste jest, że wraz Stężenie analitu wyrażone w mg/l Vk objÄ™toSć kolby, w której przygoto- z rozcieÅ„czaniem próbki (wzorca) ro- przeksztaÅ‚cane jest na mg/kg wedÅ‚ug wywano wzorzec [ml]. Snie niepewnoSć metody. By zmniej- wzoru (12) z uwzglÄ™dnieniem masy szyć zatem niepewnoSć oznaczania próbki oraz objÄ™toSci kolby, w której TÄ… zależnoSć należy wstawić tych pierwiastków, należy zastosować jÄ… przygotowano: w miejsce x1, x2, x3 do modelu mate- roztwór wzorca o stężeniu 100 mg/l. Tabela 14 x0 "Vp (12), WzglÄ™dny udziaÅ‚ poszczególnych elementów metody w niepewnoSci zÅ‚ożonej C = m Relative contribution of individual method components in complex uncertainty gdzie: Element metody B Ba Ca Cu Fe Mg Mn Sr Zn C wartoSć stężenia [mg/kg], Intensywność próbki [cps] 26,2 23,7 16,3 14,57 0,2 18,8 1,8 30,6 17,8 Intensywność Å›lepej próbki [cps] <0,1 <0,1 <0,1 <0,1 <0,1 <0,1 <0,1 <0,1 <0,1 x0 wartoSć stężenia [mg/l], Intensywność wzorca 1 [cps] <0,1 <0,1 <0,1 <0,1 <0,1 <0,1 <0,1 3,2 <0,1 Vp objÄ™toSć kolby, w której przygoto- Intensywność wzorca 2 [cps] 0,9 0,4 <0,1 2,42 <0,1 0,5 0,11 <0,1 0,1 wano próbkÄ™ [ml], Intensywność wzorca 3 [cps] 5,6 4,5 0,7 10,4 3,6 1,6 2,13 6,9 7,1 m masa próbki [g]. Intensywność wzorca 4 [cps] 0,03 0,5 7,2 2,3 4,0 6,2 5,3 1,2 0,7 ObjÄ™tość próbki [ml] 0,3 0,3 1,3 0,20 2,6 1,0 1,50 0,2 0,3 Masa próbki [g] <0,1 <0,1 0,2 <0,1 0,2 0,1 0,1 <0,1 <0,1 Po poÅ‚Ä…czeniu tych wzorów 8 12 Stężenie wzorca głównego [mg/l] 55,0 58,3 0,2 39,9 0,4 0,4 0,2 49,1 61,0 otrzymujemy model matematyczny ObjÄ™tość wzorca 1 [ml] <0,1 <0,1 <0,1 <0,1 <0,1 <0,1 <0,1 <0,1 <0,1 opisany wzorem 13. ObjÄ™tość wzorca 2 [ml] 3,6 3,5 <0,1 9,7 15,2 <0,1 0,1 1,7 3,7 ObjÄ™tość wzorca 3 [ml] 7,3 7,5 2,3 16,6 37,8 11,0 28,1 4,1 8,0 Jednak należy dodatkowo ObjÄ™tość wzorca 4 [ml] 0,7 1,0 70,4 3,6 33,4 59,5 59,1 2,7 1,0 uwzglÄ™dnić fakt, że wzorce kalibra- ObjÄ™tość roztworu wzorca [ml] 0,3 0,3 1,3 0,2 2,6 1,0 1,5 0,2 0,3 cyjne zostaÅ‚y przygotowane przez ëÅ‚ y1 + y2 + y3 + y4 (x1 - x)(y1 - y)+ (x2 - x)(y2 - y)+ (x3 - x)(y3 - y)+ (x4 - x)(y4 - y)" x1 + x2 + x3 + x4 öÅ‚ ìÅ‚ ÷Å‚ (13) Vp " (yox - y0sp)- - 2 2 2 2 ìÅ‚ ÷Å‚ 4 îÅ‚mg Å‚Å‚ (x1 - x) + (x2 - x) + (x3 - x) + (x4 - x) 4 íÅ‚ Å‚Å‚ CïÅ‚ śł = (x1 - x)(y1 - y)+ (x2 - x)(y2 - y)+ (x3 - x)(y3 - y)+ (x4 - x)(y4 - y) kg ðÅ‚ ûÅ‚ m " 2 2 2 2 (x1 - x) + (x2 - x) + (x3 - x) + (x4 - x) PROBLEMY KRYMINALISTYKI 253/06 25 Dla pozostaÅ‚ych pierwiastków Podsumowanie BIBLIOGRAFIA głównym elementem determinujÄ…cym niepewnoSć metody jest dozowanie W opublikowanych już artykuÅ‚ach 1. M. Kuras, M. Wachowicz: Profilo- objÄ™toSci roztworu wzorca głównego, zwiÄ…zanych z opracowywaniem me- wanie konopi na podstawie skÅ‚adu pier- czyli by zmniejszyć niepewnoSć me- tody profilowania konopi na podsta- wiastkowego cz. I (efekty matrycowe), tody, należy kupić pipety o wiÄ™kszej wie skÅ‚adu pierwiastkowego przed- Problemy Kryminalistyki 2006, nr 252, precyzji dozowania objÄ™toSci. stawiono: badania wstÄ™pne zwiÄ…zane s. 21 30. W tabeli 15 podano wzglÄ™dne pro- z analizÄ… materiaÅ‚u roSlinnego, stabil- 2. M. Kuras, praca magisterska Ana- centowe wartoSci niepewnoSci zÅ‚ożo- noSciÄ… stosowanych roztworów, efek- liza elementarna wybranych narkotyków nej (U) metody oznaczania każdego tami matrycowymi i walidacjÄ… meto- oraz półproduktu i produktu syntezy siar- pierwiastka. Zamieszczono w niej dy. DoprowadziÅ‚y one do opracowa- czanu 4-etoksyamfetaminy , WydziaÅ‚ również wartoSć stężenia, gdyż waż- nia wiarygodnej i pewnej metody wy- Chemii UW, Warszawa 2002. ny jest poziom stężeÅ„, dla którego znaczania skÅ‚adu pierwiastkowego 3. M. Wachowicz, M. Kuras: WstÄ™p wyznaczono niepewnoSć. Podana konopi. Na podstawie wykonanej pra- do profilowania konopi na podstawie skÅ‚a- wartoSć Ur to wartoSć wzglÄ™dna nie- cy okreSlone zostaÅ‚y zasady i para- du pierwiastkowego, Problemy Krymina- pewnoSci rozszerzonej. Obliczono jÄ…, metry mineralizacji konopi oraz wa- listyki 2003, nr 240, s. 10 19. mnożąc wartoSć U przez 2. Niepew- runki analizy technikÄ… ICP-OES. 4. M. Wachowicz, M. Kuras: Minerali- noSć rozszerzona oddaje peÅ‚niej war- Wyniki opisanych etapów pracy zacja mikrofalowa jako jedna z technik toSć niepewnoSci metody. Jest to autorzy wykorzystali w nastÄ™pnych przygotowania próbek do badaÅ„ porów- spowodowane tym, że przy szacowa- badaniach zwiÄ…zanych już z wÅ‚aSci- nawczych, Problemy Kryminalistyki niu niepewnoSci czÄ™sto nie można wym, bezpoSrednim badaniem ziela 2002, nr 238, s. 8 22. przewidzieć i okreSlić wszystkich konopi pod kÄ…tem rozkÅ‚adu pierwiast- 5. P. van Zoonen, R. Hoogerbrugge, czynników, które na niÄ… wpÅ‚ywajÄ…. ków w roSlinie (korzeniach, Å‚odygach, S.M. Gort, H.J. van de Wiel, H.A. van`t Jak można siÄ™ byÅ‚o spodziewać, liSciach, kwiatostanach oraz nasio- Klooster: Some practical examples of najwiÄ™kszÄ… wzglÄ™dnÄ… niepewnoSć nach). Ponadto przeprowadzono method validation in the analytical labora- analizy próbek kono- tory, Trends in Analytical Chemistry , 18 Tabela 15 pi włóknistych pobra- (1999), s. 584 593. WartoSci niepewnoSci standardowej (U) i rozszerzonej nych z plantacji usy- 6. Eurachem/CITAC Guide: Quantify- (Ur) metody oznaczania pierwiastków technikÄ… ICP OES tuowanych w róż- ing uncertainty in analytical measure- Values of standard (U) and extended (Ur) uncertainty nych rejonach Polski, ment, Second edition, 2000. for element determination method by ICP OES a także próbek kono- 7. B.N. Taylor, C.E. Kuyatt: Guideli- pi narkotycznych , nes for evaluating and expressing the un- które byÅ‚y jednocze- certainty of NIST measurement results, Stężenie U Ur Pierwiastek Snie przedmiotem NIST Technical Note 1297, 1994 Edition. [mg/kg] [%] [%] badaÅ„ w eksperty- 8. M. Thompson, S.L.R. Ellison, B 29,6 5,5 11,0 zach opracowywa- R. Wood: Harmonized guidelines for sin- nych w CLK KGP. gle laboratory validation of methods of Ba 31,4 5,3 10,6 Otrzymane wyniki analysis (IUPAC Technical Report), Pure Ca 11428 2,5 5,1 zostaÅ‚y poddane za- Apel. Chem. , vol. 74, nr 5 (2002), s. 835 855. Cu 7,27 6,4 12,9 awansowanej anali- 9. Eurachem/CITAC Guide: Traceabili- zie statystycznej ty in chemical measurement, a guide to Fe 423 2,8 5,6 w celu oceny możli- achieving comparable results in chemical Mg 2952 2,8 5,6 woSci grupowania measurement, 2003. Mn 184 2,9 5,8 próbek konopi. 10. Analytical Methods Committee: Opracowana me- Uncertainty of measurement: implications Sr 39,1 5,8 11,6 toda jest również of its use in analytical science, Analyst Zn 31,4 5,2 10,4 wykorzyst ywana 1995, vol. 120, s. 2303 2308. w bieżącej pracy 11. R.J.N. Bettencourt da Silva, M. Fi- metody uzyskano dla pierwiastków, WydziaÅ‚u Chemii CLK KGP przy bez- lomena, G.F.C. Camões, J. Seabra e które oznaczano na poziomie Slado- poSrednim porównywaniu konopi do- Barros: Validation and quality control wym, tj. B, Ba, Cu, Sr i Zn, i mieSciÅ‚a wodowych (użytkownik) i porównaw- schemies based on the expression of re- siÄ™ ona w zakresie 10,4÷12,9%. Dla czych (hurtownik) pod kÄ…tem ich sults with uncertainty, Analytica Chimica pierwiastków, których stężenia pier- wspólnego pochodzenia. Wyniki Acta 1999, 393, s. 167 175. wiastków w próbkach konopi sÄ… z tych etapów pracy autorzy przed- 12. W. Hyk, Z. Stojek: Analiza statystycz- znacznie wyższe, niepewnoSć roz- stawiÄ… w nastÄ™pnym numerze Pro- na w laboratorium analitycznym, Komitet szerzona wynosi od 5,1% dla wapnia blemów Kryminalistyki . Chemii Analitycznej PAN, Warszawa 2000. 13. ICP WinLab Bonus Pack, Perkin Elmer. do 5,8% dla manganu. 26 PROBLEMY KRYMINALISTYKI 253/06