INTRODUCCIÓN
Los modelos analíticos de calibración basados en espectroscopia ultravioleta-visible (UV-Vis) para la determinación de los hidrocarburos totales del petróleo (HTP) en suelos y de hidrocarburos aromáticos (HA) en aceites vegetales, alimentos asados y agua, son una alternativa rápida para el análisis químico cuantitativo de muestras problema y han demostrado ser adecuados para este fin. La construcción de modelos de calibración lineal simple (MCLS) se lleva a cabo mediante una calibración en la que se emplean estándares del analito de interés y permite garantizar una relación fiable entre la cantidad de éste y una señal dada por el instrumento (Valcárcel y Cárdenas 2000). No debe confundirse con la calibración instrumental, la cual se realiza para asegurar el correcto funcionamiento del equipo empleando una sustancia diferente al analito de interés.
Con relación el método para la cuantificación de HTP en suelo, Wang et al. (2011) validaron un modelo lineal basándose en el coeficiente de correlación. A su vez, Dost e Ideli (2012) obtuvieron sus modelos de calibración mediante cromatografía líquida de alta resolución (HPLC) equipada con un detector UV-Vis. La linealidad de dichos modelos, desarrollados para la cuantificación de HA, fue evaluada por el método de mínimos cuadrados.
Si bien la mayoría de las técnicas convencionales de cuantificación de HTP se basan en el método gravimétrico por su sencillez, relativa rapidez y bajo costo, éste sólo es recomendable para suelos que contienen niveles superiores a 50 000 mg/kg de HTP o de asfaltenos. En suelos con bajas concentraciones de hidrocarburos disminuye la precisión del método debido al error que implica el manejo de pesos pequeños (Fernández et al. 2006).
En estudios cuyo objetivo es el análisis cuantitativo o identificación de compuestos específicos en mezclas complejas de hidrocarburos como los aromáticos policíclicos, comúnmente se aplican técnicas como la cromatografía de gases (GC) o la HPLC, entre otras (Dost e Ideli 2012), pero estos métodos requieren largos tiempos de análisis por la necesidad de extracción de compuestos específicos, alto consumo de reactivos y elevado costo de los instrumentos.
Los modelos de calibración obtenidos por espectrofotometría son resultado de una calibración indirecta, ya que la cantidad del analito presente en las soluciones de referencia se determina en función de la señal dada por el instrumento analítico (Cuadros-Rodríguez et al. 2001).
En los métodos de calibración por UV-Vis reportados, las relaciones entre las concentraciones de referencia de hidrocarburos y la respuesta del espectrofotómetro son descritas por modelos lineales simples. Éstos, además de la estadística descriptiva, consideran como parámetros de validación de sus modelos de calibración el coeficiente de determinación (R2), el coeficiente de correlación (r), los límites teóricos de detección y de cuantificación, y su precisión en condiciones de reproducibilidad (Wang et al. 2011, Dost e Ideli 2012).
Barwick (2003) establece paso a paso el proceso para la construcción de calibraciones, los parámetros estadísticos y la validación interna de los modelos de calibración, utilizando los siguientes indicadores: 1) el cálculo y representación gráfica de los residuales y 2) los estadísticos de la regresión, como r, , R2, la desviación estándar de los residuales o error estándar, el análisis de varianza (ANDEVA) y los coeficientes de regresión m y a (pendiente y ordenada al origen), los cuales no son considerados en su totalidad en los estudios realizados por Vogt et al. (2000), Wang et al. (2011) y Dost e Ideli (2012).
Este estudio presenta una propuesta para la validación interna y externa de modelos de calibración lineal simple obtenidos por espectrofotometría UV-Vis, para determinar los HTP en suelos o en agua extraídos con diferentes disolventes y en niveles que no superan los límites máximos permitidos (3000 mg/kg base seca de hidrocarburo de la fracción pesada, en suelos agrícolas) establecidos por la NOM-138-SEMARNAT/SSA1-2012 (SEMARNAT 2013).
MATERIALES Y MÉTODOS
Petróleos crudos y disolventes orgánicos
Se utilizaron dos tipos de petróleo crudo, ligero (PCL) y pesado (PCP), proporcionados por Petróleo Mexicanos, cuyas densidades son similares a las reportadas por Santos et al. (2014): 0.83-0.87 y 0.92-1.0 g/cm3, respectivamente. Las densidades correspondientes por grados del American Petroleum Institute (API) son de 31.1 a 39.0 y de 10.2 a 22.3 ºAPI, respectivamente (IMP 2014).
Los disolventes empleados en la preparación de las soluciones de referencia o estándares de los petróleos crudos fueron hexano ultrapuro y éter etílico grado ACS. El primero tiene una polaridad más baja que el segundo con base en sus constantes dieléctricas: 1.9 y 4.3, respectivamente (Masschelein-Kleiner 2004).
Espectrofotómetro UV-Vis
Para las determinaciones se utilizó un espectrofotómetro UV-Vis marca Perkin Elmer, modelo Lambda 25 de doble haz, con un intervalo de escaneo de 200 a 1100 nm. Los registros de los espectros fueron realizados cada 2 nm con una velocidad de escaneo de 240 nm/min.
Análisis preliminar de un calibrado de petróleo crudo ligero mediante UV-Vis
En un análisis preliminar se obtuvo el espectro completo de la región UV-Vis (200 a 1100 nm) de soluciones de referencia de concentraciones de 0 a 1000 mg/L. Esto se realizó con la finalidad de encontrar una longitud de onda (λ) en el espectro UV-Vis en la que se observara una respuesta lineal entre la concentración y la absorbancia.
Una vez que se obtuvo la λ esperada, se procedió a desarrollar la metodología, aplicando dicha λ para la obtención de los modelos con las combinaciones restantes de éter etílico y el PCP. Para esta etapa preliminar se prepararon 17 concentraciones (número mayor de siete, que es el mínimo indicado por Barwick [2003]) de 0, 1, 2, 4, 6, 8, 10, 20, 40, 60, 80, 100, 200, 400, 600, 800 y 1000 mg/L de PCL. Las soluciones se prepararon en matraces volumétricos clase A de 25 mL.
Las concentraciones de 0 a 10, de 20 a 100 y de 200 a 1000 mg/L se prepararon a partir de soluciones stock de 10, 100 y 1000 de PCL, respectivamente, todas a un volumen final de 25 mL diluido con hexano (el disolvente).
Preparación y análisis por UV-Vis de los calibrados de petróleo crudo en dos diferentes disolventes
Al incluir los dos disolventes y los dos tipos de petróleo crudo en esta segunda etapa de la experimentación, se obtuvieron las siguientes combinaciones: 1) hexano-petróleo crudo ligero, 2) hexano-PCP, 3) éter-PCL y 4) éter-PCP. A partir de éstas se obtuvieron cuatro calibrados.
Cada calibrado se preparó con ocho concentraciones de referencia: 0, 4, 8, 20, 40, 60, 80 y 100 mg/L (el mínimo indicado por Barwick [2003] es siete), con cinco repeticiones de cada una para tener un total de 40 soluciones por calibrado.
La absorbancia de todas las soluciones (160 de los cuatro calibrados) se midió en un espectrómetro UV-Vis a 254 nm. En este valor de λ se encontró la mejor correlación lineal entre la concentración de petróleo y la absorbancia, definida en la etapa preliminar. El ajuste del cero de absorbancia se realizó con agua desionizada.
Validación de los modelos de calibración
Para llevar a cabo la validación primero se dividieron en dos grupos los resultados de las absorbancias de las 40 soluciones de referencia de cada combinación. Para el primer grupo, denominado grupo de entrenamiento, sólo se tomaron los resultados de 35 soluciones, equivalentes al 90 %, y se utilizaron para la construcción del modelo mediante una regresión lineal simple de los puntos graficados con objeto de realizar la calibración o validación interna.
Para el segundo grupo, denominado grupo de validación, se tomaron las absorbancias de las cinco soluciones de referencia (equivalente al 10 % de las soluciones) omitidas en el grupo de entrenamiento. Este grupo se utilizó para la validación externa del modelo.
Validación interna
La validación interna de los modelos de calibración se llevó a cabo aplicando los parámetros que se describen a continuación.
Validación cruzada
Se dividieron las absorbancias de las 35 soluciones de referencia del grupo de entrenamiento en 35 grupos, de los cuales 34 se tomaron para establecer una primera calibración mediante regresión lineal simple. El otro grupo se tomó para la predicción o validación.
El proceso se repitió iterativamente hasta que cada una de las concentraciones de referencia fue predicha con el resto de los grupos de calibración.
Coeficiente de determinación (R2)
Se obtuvieron los R2 de los cuatro modelos de regresión de los resultados de las 35 soluciones de referencia correspondientes al grupo de entrenamiento, no así los de la validación cruzada.
ANDEVA de la regresión
Se realizó el ANDEVA de la regresión de las soluciones del grupo de entrenamiento. Se estableció la hipótesis nula de que la absorbancia (y) no es una función de la concentración (x), frente a la hipótesis alternativa de que y sí está en función de las variaciones de x, aceptándose esta última si F > F crítica. Al aceptarse la hipótesis alternativa se indica que el modelo lineal es significativo (Weisberg 2005).
Residuales y error absoluto medio (EAM)
Se presentan los gráficos de los residuales y los valores del EAM de los cuatro modelos.
Coeficiente de correlación (r)
Se propone utilizar r para explicar el grado de asociación entre las concentraciones estimadas por dichos modelos y el conjunto de concentraciones de referencia correspondientes, para medir su precisión.
Prueba t de Student
Se propone realizar una prueba t de Student para datos apareados con el objeto de definir si existen o no diferencias estadísticas entre los mismos valores que se usaron para r.
La prueba t se realizó con una probabilidad del 5 %, con la hipótesis nula de que no existen diferencias estadísticas entre las concentraciones de referencia y las concentraciones estimadas por los modelos, frente a la hipótesis alternativa de que sí existe diferencia entre los valores de referencia y los estimados, aceptándose la primera si p-valor > 0.05.
Los estadísticos empleados para la validación de los modelos se describen a continuación (Marqués 2010).
Coeficiente de determinación (R 2 ) del modelo de regresión lineal
donde r es el coeficiente de correlación para las variables de las curvas de calibración
donde x i es el valor observado en el eje x, y i el valor observado en el eje y, la media de los valores x i ; la media de los valores y i , y n el número de puntos de calibración.
Análisis de varianza (ANDEVA) de la regresión de las curvas de calibración
Suma de cuadrados de los residuales:
Suma de cuadrados de la regresión:
Pendiente (m) de la línea de regresión:
Ordenada al origen (a) de la línea de regresión
donde = valores de y estimados por la ecuación de regresión para un valor dado de x.
Prueba t de Student para datos apareados entre concentraciones estimadas por los modelos y las concentraciones de referencia
donde x̅ d es la media de las diferencias de los valores de las variables, s d el número total de diferencias y la desviación estándar de las diferencias de los valores de las variables.
Validación externa
Durante la validación externa se estimaron las cinco concentraciones de referencia del grupo de validación utilizando el modelo construido con las 35 concentraciones del grupo de entrenamiento durante la validación interna. De esta manera se tuvo un conjunto de concentraciones de referencia y otro de concentraciones estimadas, para cada modelo de calibración.
Límites de detección y cuantificación de los cuatro modelos de calibración
Se calcularon los límites de detección y de cuantificación de cada modelo. Primero se obtuvo la media del conjunto de concentraciones del blanco (concentración 0 mg/L) estimadas por los modelos. El límite de detección se obtuvo sumando a la media de las concentraciones del blanco, tres veces su desviación estándar.
El límite de cuantificación se calculó sumando a la media de las concentraciones del blanco 10 veces su desviación estándar (Shrivastava y Gupta 2011).
RESULTADOS Y DISCUSIÓN
Análisis preliminar de un calibrado de petróleo crudo mediante espectrofotometría
La obtención del espectro completo de la región UV-Vis (Fig. 1) de las soluciones de referencia (0 a 1000 mg/L) de PCL en hexano, permitió identificar que en la λ de 254 nm las soluciones presentan un pico de absorción con respuesta lineal entre la concentración de PCL y la absorbancia. Este resultado fue ratificado repitiendo la medición de dichas soluciones, pero únicamente a una λ de 254 nm, también en el intervalo de concentración de 0 a 1000 mg/L de PCL. A partir de las absorbancias obtenidas del segundo análisis se construyó el gráfico con la recta ajustada por regresión, la cual se muestra en la figura 2.
La λ de 254 nm seleccionada en el presente estudio coincide con la reportada por Dost e Ideli (2012), quienes observaron que es la mejor respuesta lineal para nueve hidrocarburos aromáticos policíclicos (HAP) en la región UV. Por su parte, Wang et al. (2011) observaron una absorción estable a un valor de λ de 230 y 200 nm para una mezcla de hidrocarburos disueltos en diclorometano y hexano, respectivamente.
Preparación y análisis por UV-Vis de los calibrados de petróleo crudo en diferentes disolventes
En la figura 3 se presentan los gráficos de las cuatro curvas de calibración, con cinco repeticiones para cada concentración de los petróleos crudos. La representación gráfica de los puntos de calibración (absorbancia-concentración de los petróleos crudos) es una opción rápida para observar el comportamiento de los datos en cuanto a linealidad, y constituye una forma cualitativa de evaluarla. También se observa el efecto de la combinación de disolventes-petróleos crudos en las diferentes pendientes (m) de cada caso.
Las absorbancias de los calibrados de PCL resultaron menores debido a la baja densidad de éste respecto del PCP, cuyas absorbancias fueron mayores. Además, el PCL contiene hidrocarburos de menor peso molecular, a diferencia del PCP. Todo esto influyó también en las variaciones observadas entre las m de los modelos ambos petróleos.
Validación de los modelos
Coeficiente de determinación (R2) del modelo de regresión lineal
Como resultado del proceso de validación interna se construyeron cuatro modelos lineales entre la absorbancia a 254 nm y las concentraciones de petróleos crudos (Fig. 4) con valores de R2 de 0.9952 a 0.9996. Esto significa que los modelos explican en promedio el 99.74 % de la variabilidad de los resultados. Valores de R2 similares fueron obtenidos por Dost e Ideli (2012) en la validación de un método de detección de hidrocarburos aromáticos policíclicos (HAP) por HPLC/UV-Vis; y por Wang et al. (2011) en mezclas de hidrocarburos de petróleo.
En la figura 4 se observa también que para los modelos de calibración con PCL y PCP, el disolvente usado influye mediante la modificación de los valores de a en las ecuaciones de regresión, siendo en ambos casos mayores con éter etílico que con hexano. Las variaciones en los valores de a pueden atribuirse a la densidad de los disolventes, ya que el éter etílico es más denso que el hexano (713.4 kg/m3 y 654.8 kg/m3, respectivamente).
Por otra parte, los valores de m son influenciados por el tipo de hidrocarburo, siendo mayores con el PCP en ambos disolventes. En el caso de los modelos hechos con PCL, m conserva su valor indistintamente del disolvente usado (0.004) y en los modelos hechos con PCP, m mostró un ligero aumento con el hexano respecto al éter etílico (0.0169 y 0.0156, respectivamente). Aun así, son valores muy cercanos entre sí.
Los cambios en los valores de m en las curvas debidos al tipo de petróleo se pueden explicar por las características físicas de cada uno, específicamente por su densidad. El PCP es más denso que el PCL, por lo que las absorbancias fueron mayores para los calibrados con PCP, llegando hasta 1.7 a 100 mg/L, mientras que los calibrados con PCL mostraron una absorbancia de 0.5 a la misma concentración.
Al tratarse de cuatro combinaciones de disolvente-petróleo crudo preparados en condiciones e intervalos de concentración estándar iguales, las variaciones observadas en las absorbancias evidencian la existencia de efectos matriciales por la composición de las soluciones, mismas que afectan directamente los valores de los coeficientes de los modelos. El efecto matricial en los modelos de calibración se descarta al ajustar el cero de absorbancia empleando el blanco de reactivos de los dos disolventes usados.
Danzer y Currie (1998) afirman que el efecto matricial de las muestras problema puede contrarrestarse utilizando el método de adición de estándar, el cual consiste en adicionar soluciones estándar a la muestra para crear un comportamiento similar del conjunto de calibración y de la muestra siempre que el analito se añada en la forma de la misma especie.
En la figura 4 se observa que el intervalo de confianza (IC) al 99 % (región enmarcada por las bandas punteadas a ambos lados de la recta ajustada) es muy estrecho para las cuatro líneas ajustadas, especialmente para los calibrados en que se usó hexano como disolvente, lo cual coincide con los bajos valores de los residuales. Estos resultados son congruentes con los conceptos de IC reportados por Newcombe (1998) y Barwick (2003).
ANDEVA de la regresión
El análisis de varianza es considerado un método no gráfico para evaluar la linealidad de un modelo de calibración (Raposo 2016). El ANDEVA de los modelos de regresión de los cuatro calibrados (Cuadro I) muestra que las sumas de cuadrados de las regresiones son superiores a la suma de cuadrados de los residuales en todos los casos. Para magnificar estas diferencias se tiene el valor del estadístico F, que resultó significativamente mayor que su valor crítico, lo cual es indicativo de que los cuatro modelos son significativamente lineales.
Hexano-PCL | |||||
Grados de libertad | Suma de cuadrados | Promedio de los cuadrados | F | Valor crítico de F | |
Regresión | 1 | 0.71846936 | 0.71846936 | 86096.4487 | 5.8204E-58 |
Residuos | 33 | 0.00027538 | 8.3449E-06 | ||
Total | 34 | 0.71874474 | |||
Hexano-PCP | |||||
Regresión | 1 | 12.0796147 | 12.0796147 | 42363.2316 | 6.9845E-53 |
Residuos | 33 | 0.00940975 | 0.00028514 | ||
Total | 34 | 12.0890244 | |||
Éter-PCL | |||||
Regresión | 1 | 0.73864744 | 0.73864744 | 6878.49859 | 6.9648E-40 |
Residuos | 33 | 0.0035437 | 0.00010738 | ||
Total | 34 | 0.74219114 | |||
Éter-PCP | |||||
Regresión | 1 | 11.2470541 | 11.2470541 | 15714.2794 | 8.738E-46 |
Residuos | 33 | 0.02361882 | 0.00071572 | ||
Total | 34 | 11.270673 |
PCL: petróleo crudo ligero, PCP: petróleo crudo pesado
El modelo del calibrado hexano-PCL es capaz de describir el 99.96 % de la variabilidad de los datos, el del hexano-PCP el 99.92 %, el del éter-PCL el 99.52 % y el del éter-PCP el 99.79 %.
Residuales y EAM
De acuerdo con la figura 5, la distribución de los residuales indica un comportamiento más aleatorio para el calibrado hexano-PCL que para el resto, destacando además que el intervalo de dichos valores fue el más cercano a cero (0.0122 a -0.0033).
Los intervalos correspondientes a los residuales de los calibrados éter-PCL, hexano-PCP y éter-PCP presentaron una dispersión respecto al cero del orden siguiente: 0.0209 a -0.0215, 0.0450 a -0.0359 y 0.0557 a -0.0663, respectivamente.
El EAM o promedio de los valores absolutos de los residuales se cataloga como una medida de desviación (Mayer y Butler 1993) aplicada en la validación de modelos de calibración. Los valores para los calibrados hexano-PCL, hexano-PCP, éter-PCL y éter-PCP fueron 0.0017, 0.0133, 0.0081 y 0.0214, respectivamente.
Prueba t de Student para datos apareados entre concentraciones estimadas y concentraciones de referencia
Los resultados de la prueba t de Student para datos apareados o relacionados mostraron que no existió diferencia significativa (p > 0.05) entre las concentraciones estimadas por los modelos y las concentraciones de referencia de las muestras del grupo de entrenamiento, en los cuatro modelos analizados.
Coeficiente de correlación (r) entre concentraciones estimadas y concentraciones de referencia
Los valores de r entre las concentraciones estimadas y las de referencia (Fig. 6) obtenidos en la validación interna, fueron superiores a 0.9996 para los calibrados hexano-PCL y hexano-PCP, y en el caso de los calibrados éter-PCL y éter-PCP fueron superiores a 0.9976. En todos los casos estos resultados indican una alta asociación entre dichas concentraciones como lo indica el valor cercano a 1. El uso de r como medida de la linealidad del modelo es ampliamente reconocido (Barwick 2003, Raposo 2016). La correlación entre las variables analizadas puede ser o no lineal (van Loco et al. 2002).
Validación externa de los modelos de calibración
Prueba t de Student para datos apareados
Con la prueba t de Student se encontró que no existió diferencia estadística (p > 0.05) entre las medias de las concentraciones estimadas por los modelos y las concentraciones de referencia del conjunto del 10 % de las mediciones que se omitieron en la construcción y calibración de los modelos.
Coeficiente de correlación (r) entre concentraciones estimadas y concentraciones de referencia
Los valores de r obtenidos en la validación externa (Fig. 7) para las concentraciones estimadas con los modelos y las concentraciones de referencia del grupo de validación, fueron superiores a 0.9993 para los calibrados hexano-PCL, hexano-PCP y éter-PCP, e igual a 0.9983 para el calibrado éter-PCL. En todos los casos existe un alto grado de asociación entre ambos conjuntos de concentraciones por los valores cercanos a 1.
Límites de detección y cuantificación de los calibrados
Una vez validados los modelos se calcularon los límites de detección y los límites de cuantificación a partir de los datos obtenidos con los calibrados, para cada uno de los modelos de calibración. Dichos límites se muestran en el cuadro II, que también incluye las medias y desviaciones estándar de los blancos.
Hexano-PCL | Hexano-PCP | Éter-PCL | Éter-PCP | |
Límite de detección (mg/L) | 0.82 | 1.03 | 3.70 | 2.81 |
Límite de cuantificación (mg/L) | 3.72 | 1.21 | 6.98 | 4.7 |
Media de los valores estimados del blanco | -0.425 | 0.9556 | 2.3 | 1.9743 |
DE de los valores estimados del blanco | 0.4145 | 0.0256 | 0.4677 | 0.2775 |
PCL: petróleo crudo ligero, PCP: petróleo crudo pesado, DE: desviación estándar
Los valores más bajos en el límite de detección y de cuantificación correspondieron a ambos tipos de petróleo crudo con hexano, mientras que para los calibrados con éter fueron mayores. Este resultado coincide con los valores de las ordenadas al origen en los modelos mostrados en la figura 4.
En contraste, los valores encontrados para modelos generados para HAP utilizando concentraciones estándares de HA (Dost e Ideli 2012) y los límites de detección para HA (Vogt et al. 2000) fueron menores a los de este estudio. La diferencia en dichos valores respecto a los encontrados en este trabajo puede deberse a que los autores emplearon compuestos puros para sus análisis y en este trabajo se emplearon petróleos crudos cuya composición es más compleja.
CONCLUSIÓN
La evaluación interna y externa de los MCLS propuesta para calificar las calibraciones por espectrofotometría UV-Vis indican que las estimaciones para cuantificar HTP en suelos o en agua, extraído con diferentes disolventes, son adecuadas para su aplicación en concentraciones menores a las establecidas como límites máximos permitidos en la normativa mexicana. La linealidad de los modelos se demostró cuantitativamente en la validación interna mediante un valor del coeficiente de determinación muy cercano de la unidad; asimismo, por medio del ANDEVA, en el que el modelo lineal resultó altamente significativo, y la representación gráfica de los datos experimentales, que es una demostración cualitativa de la linealidad. Los residuales presentaron una distribución aleatoria cercana al cero.
La capacidad para estimar la concentración de HTP de los modelos se evaluó en la validación externa mediante las pruebas t de Student y de correlación. Con la primera se demostró que no hay diferencias significativas entre los datos estimados con los modelos y los valores de referencia que se omitieron en la construcción y calibración de los modelos; y la segunda evidenció que hay un alto grado de asociación entre ambos conjuntos de datos. Además, los límites de detección y de cuantificación variaron en función del tipo de disolvente empleado (efecto matricial), siendo ambos valores más altos en los calibrados con éter etílico que con hexano, lo cual coincide con las ordenadas al origen en los modelos. Todos los indicadores estadísticos usados demuestran la robustez de los modelos propuestos y son la base para la evaluación de calibraciones lineales.