Quimica Incertidumbre
Quimica Incertidumbre
Quimica Incertidumbre
[4]
siendo:
s la desviacin estndar de repeticiones actuales del ensayo o correspondiente a valores
histricos
n la cantidad de valores de repeticiones a partir de la cual se hace la estimacin del valor del
mensurando (y NO de la cantidad de valores con los cuales se calcula la desviacin estndar)
3 Se establece el modelo matemtico:
y = f(x1, x2, ... , xn) [5]
en donde:
y: valor verdadero del mensurando
xk: valor verdadero de la k-sima magnitud de entrada
4 Se calcula la incertidumbre estndar combinada mediante us mediante la frmula de propagacin:
( )
2 2 2 2 2 2
1 1 2 2 12 12 13 13
2
... ...
cov ,
s n n
k
k
ij
i j
ij i j
u c u c u c u c u c u
Y
c
X
Y
c
X X
u X X
+ + + + + +
[6]
siendo:
uk:la componente sistemtica de la incertidumbre estndar combinada correspondiente a la k-
sima magnitud de entrada
uij: la covarianza entre la i-sima y la j-sima magnitudes de entrada; si son independientes, uij= 0
ck: coeficiente de sensibilidad respecto a la k-sima magnitud de entrada
5 Para calcular el nmero efectivo de grados de libertad, ef, se utiliza la siguiente frmula de
propagacin (a menos que se desee usar la aproximacin k0.95 2 mencionada anteriormente):
4
f 4 4 4
a
1
1
e
n
k k
k k
u
c u u
r
[7]
siendo:
r: nmero de repeticiones del ensayo (necesariamente r > 1)
k: nmero efectivo de grados de libertad correspondiente a la k-sima magnitud de entrada
Los dems smbolos tienen el significado establecido anteriormente
6 Para una medicin directa de una magnitud x se debe considerar el siguiente modelo de error:
x = + 1 + ... + m [8]
siendo:
: valor verdadero del mensurando
1, ... , m: componentes de error (independientes)
Se tiene:
5
Page of 10
XVII Jornadas IRAM-Universidades
UNLM, Octubre 2002
2 2 2
1 2
4
f 4
1
...
m
e
m
k
k k
u u u u
u
u
+ + +
[9]
aqu:
u: componente sistemtica de incertidumbre correspondiente a x
uk: componente de incertidumbre debida a k
7 Para una componente de error, se tiene:
a) Si el error puede tomar cualquier valor slo en un intervalo [-a; a] (donde a puede ser
conocido o desconocido), al mismo se le atribuye una distribucin uniforme o triangular.
b) La distribucin triangular es la ms conservadora, a los efectos del clculo de incertidumbre,
en el sentido de que conduce a valores mayores. La distribucin triangular es utilizada
cuando se considera que los errores estn concentrados cerca del 0.
c) Clculos:
( )
( )
( )
2
Uni , , evaluacin tipo B
3
Triang , , evaluacin tipo B
6
si la distribucin del error es conocida
si la distribucin del error no es conocida
( es el estimador de la incert 2
a
si a a
u
a
si a a
u
V
'
u u
:
:
idumbre)
'
[10]
EL CASO DE LOS ENSAYOS QUIMICO ANALITICOS
El estudio de las incertidumbres en los ensayos qumico analticos presenta dificultades particulares.
La posible complejidad de un procedimiento tcnico que puede comenzar con el muestreo, la apertura de la
muestra, continuar con el posible efecto de matriz antes de llegar a la obtencin de un resultado y a la
estimacin de la incertidumbre hace parecer muy complicado este enfoque. Tambin, lleva a preguntarse
cuan significativo puede ser el resultado numrico obtenido.
La estimacin de las incertidumbres en los ensayos qumicos se ve ayudada por la publicacin de una gua
especfica de Eurachem (5) basada en la GUM (4).
Antes de incursionarse en los clculos, es muy importante tener una visin cualitativa del problema lo que
ayuda sobremanera al momento de tener que identificar cuales son las componentes no cuantificables o si
hay algunas de ellas que sean despreciables.
Nuestra experiencia
En el Departamento de Anlisis del Instituto de Ingeniera Qumica en la Facultad de Ingeniera se han
validado varios procedimientos, hacindose la estimacin de la incertidumbre del mtodo. En todos los
casos, fue un trabajo interdisciplinario con el Departamento de Matemtica de la Facultad de Qumica. A
estos efectos, se estudian todas las fuentes de incertidumbre que puedan aportar al resultado final. Ese
estudio se realiza tanto observando cuales son los parmetros que aparecen en la expresin del clculo del
resultado como siguiendo los pasos experimentales, buscando qu otras fuentes de incertidumbre pueden
influir en ese resultado aunque no aparezcan en la frmula de clculo (por ejemplo variaciones de
temperatura o de humedad). Luego de hacer ese despliegue, se busca que una misma fuente de incertidumbre
no se repita, se ve cuales se pueden estimar, cuales pueden ser desestimadas.
Algunos de los casos estudiados ya han sido publicados (6), (7).
Enfoque cualitativo
Deben buscarse todas las fuentes de incertidumbre relevantes para el ensayo en estudio; esto es til tanto para
comprender el proceso como para evitar la sobrecobertura al hacer el clculo.
6
Page of 10
Incertidumbre
Derrgibus. Fuentes
Una manera de realizarlo es utilizando un diagrama causa-efecto o de Ishikawa o espina-pez. Para
establecerlo se debe:
- escribir la ecuacin completa y con sus parmetros formar las ramas principales del diagrama (a veces es
necesario agregar una rama principal correspondiente a la correccin general por el sesgo); esto permite
identificar los efectos sobre un resultado
- estudiar los pasos del mtodo y agregar otros factores tales como condiciones ambientales o efecto de
matriz
- para cada uno de los parmetros, agregar ramas secundarias con los factores que tienen influencia sobre
ella hasta considerar que los efectos son despreciables
- resolver duplicaciones, reordenar, agrupar causas relacionadas entre s
- afinar, refinar
Esto puede conducir a ms de un tipo de arreglo posible.
De la observacin del diagrama, puede volverse a la ecuacin para agregar algn parmetro si fuera
necesario.
Esto permite tambin integrar parmetros que luego no aparecen en la expresin de la cuantificacin de las
incertidumbres, pero que deben estar bajo control para que el mtodo lo est. Es el caso, por ejemplo, del
control de la temperatura ambiente o el control de la temperatura de un digestor de muestras.
Algunos datos
Se presentan los resultados del estudio de algunos procedimientos tcnicos. En todos los casos por nosotros
estudiados se aplic una metodologa similar y, dentro de lo posible, se trabaj con el mismo equipamiento..
Elegido el procedimiento analtico, se consideran las siguientes fuentes que pueden afectar la incertidumbre
del resultado:
- pesada en la balanza analtica
- pureza de diferentes reactivos
- utilizacin de matraces aforados de diferentes volmenes
- utilizacin de pipetas aforadas y graduadas de diferentes volmenes
- escaln de la bureta
- otras particulares del mtodo en estudio
Se cuantifican para establecer cuales son aquellos componentes cuyas incertidumbres sean despreciables y
para cuales otros va a ser necesario evaluarlos y, de ser posible, minimizarlos.
Balanza Se trabaja con una balanza de precisin a la dcima de miligramo calibrada Mettler AB 204-S
Precisin: = 1 10
-4
g
Mximo error de linealidad estimado: 2 10
-4
g
Mximo error de excentricidad de carga, estimado: 2 10
-4
g
Incertidumbre (considerando error sistemtico de ajuste despreciable):
( )
2 2
4 4
2
4 4
210 210
110 g 1.6310 g
3 2 3
u
| ` | `
+ +
. , . ,
Pureza de los reactivos Es una componente sistemtica que se evala con la informacin proporcionada por
el fabricante.
Material volumtrico Se trabaja con material volumtrico de clase A. Se presentan a modo de ejemplo,
algunos de ellos con las tolerancias establecidas por sus fabricantes (Cuadro N 1).
Cuadro N 1- Caractersticas del material volumtrico
Material Tipo Volumen (ml) Tolerancia (ml) Error mximo
relativo permitido
Incertidumbre
Matraz 50
t 0,05
0,001 0,029
7
Page of 10
XVII Jornadas IRAM-Universidades
UNLM, Octubre 2002
200
t 0,10
0,0005 0,058
1000
t 0,40
0,0004 0,023
Pipeta Graduada 5
t 0,05
0,01 0,029
Aforada 10
t 0,02
0,002 0,012
Aforada 20
t 0,03
0,0015 0,017
Aforada 25
t 0,03
0,0012 0,017
Las incertidumbres fueron estimadas mediante una evaluacin tipo B considerando una distribucin
rectangular.
Los datos presentados permiten concluir que, a priori, es aceptable utilizar los valores de tolerancia al
momento de estimar las incertidumbres respectivas. Esto debe ser controlado luego de tener cuantificados
todos los componentes de la incertidumbre, bsicamente la aleatoria.
La excepcin en ese cuadro, podra ser la pipeta graduada de 5 ml.
Caso 1
Se presenta el clculo de las incertidumbres realizado en el proceso de validacin del procedimiento de
anlisis por el mtodo del reflujo cerrado y volumetra de la Demanda Qumica de Oxgeno (DQO) total en
efluentes.
La importancia del error relativo de la pipeta graduada (tal como se indica en el Cuadro N 1) que se utiliza
tanto en la medida del volumen de muestra como en la del reactivo hace necesaria su calibracin. Como
todas las medidas se hacen en el rango 0 y 2,5 ml, es el que se calibra tal como se indica en el Cuadro N 2;
los resultados presentados corresponden a 3 rplicas pesadas 3 veces cada una.
En los casos en que se calibra material de vidrio, se trabaja sobre la base de la Norma ISO correspondiente
(8).
Cuadro N 2- Calibracin de una pipeta de 5 ml graduada en 0,1 (rango: 0- 2,5 ml)
V nominal V calibracin
|Diferencia|
Error relativo
2,5 2,522389 0,022389 0,00888
De donde se concluye que el error estimado por calibracin es considerablemente menor que la tolerancia
indicada de fbrica. De esta manera, se mejora la incertidumbre sobre el resultado al disminuir el error
relativo en un instrumento cuya medida influye en el clculo final.
En el caso de las medidas de los gastos de la solucin de titulacin se utiliza una bureta de 10 ml graduada en
0,05 ml la que fue calibrada segn se indica en el Cuadro N 3.
Cuadro N 3- Calibracin de una bureta de 10 ml graduada en 0,05 ml
V (ml) Diferencia* Desviacin estndar
0-2 -0.0011 1.78E-02
2-4 0.0124 1.10E-02
4-6 -0.0289 4.19E-03
6-8 0.0187 9.26E-03
8-10 0.0077 1.28E-03
* Diferencia entre el valor nominal y el medido
De la observacin de la segunda columna se puede ver que las diferencias obtenidas por calibracin son
considerablemente menores que el escaln de la bureta. Los errores relativos son despreciables por lo que no
se hace correccin de volumen en los gastos de la titulacin a excepcin del intervalo 4 - 6 ml.
Estimacin de la componente de incertidumbre por no corregir el resultado:
cal
0.0289
ml 0.017ml
3
u
Componente aleatoria de la incertidumbre (caso ms desfavorable):
2
aleat
10 ml u
8
Page of 10
Incertidumbre
Derrgibus. Fuentes
La componente aleatoria de la incertidumbre se obtiene midiendo soluciones de biftalato de potasio (KHP)
ppa a diferentes concentraciones. Se hacen 5 medidas para cada concentracin por dos analistas diferentes.
En el Cuadro N 4 se presentan los resultados
Cuadro N 4- Volumetras de soluciones de KHP
Conc.
Nominal
(mgO2/l)
Conc. Experimental
(mgO2/l)
Desv. Est. Incertidumbre Coef. Var.
100 87,7 13,8 6,2 0,070
200 188,7 15,4 6,9 0,036
300 281,1 8,9 4,0 0,014
Para el caso en estudio, los principales componentes de la incertidumbre en el resultado son:
1) la precisin de las medidas
2) el error sistemtico en la pipeta graduada
3) el error sistemtico en la bureta
La pipeta graduada fue calibrada y se muestra que su error sistemtico es mucho menos importante que el
indicado por la tolerancia, es decir que es irrelevante.
Caso 2
Se trata tambin de una volumetra por retorno pero aplicada a la determinacin de cloruros en hormign por
el mtodo de Volhard.
Se trabaj, dentro de lo posible, con el mismo equipamiento que en el caso anterior y en las mismas
condiciones experimentales para la determinacin de la precisin.
Se realizaron otras calibraciones en los casos que se encontraron que los valores de la incertidumbre
calculada con los datos de las tolerancias era elevada.
Calibracin de un matraz aforado de 100 ml
Desviacin estndar: 4.9 10
-2
g
Incertidumbre al considerar correccin de volumen (de acuerdo a la calibracin del matraz):
2
corr
4.910 ml u
Incertidumbre al no considerar correccin de volumen (es decir, trabajando con matraz verificado):
2 mat
mat
0.08
ml 4.610 ml
3 3
T
u
en donde umat es la incertidumbre en el contenido del matraz al no aplicar la correccin, y Tmat es la tolerancia
de clase; la incidencia de umat en la incertidumbre del mtodo es despreciable, por lo cual no se efectuaron
correciones al valor nominal.
Calibracin de pipeta aforada de 5 ml
Cantidad de medidas de calibracin: 12
Promedio de las medidas: 4.9996 ml
Desviacin estndar de calibracin: 1.6310
-4
ml
Estimacin de la componente sistemtica de la incertidumbre de la pipeta (al utilizarla sin correccin):
4
sist
2.310 ml
3
E
u
siendo E=-0.0004 ml la estimacin del error sistemtico de la pipeta
Por lo tanto, esta componente de incertidumbre resulta tambin despreciable.
9
Page of 10
XVII Jornadas IRAM-Universidades
UNLM, Octubre 2002
Medida de cloruros Se realiza la apertura cida de la muestra de hormign en 5 rplicas y se miden los
cloruros en las soluciones resultantes. Los resultados se presentan en el Cuadro N 5.
Cuadro N 5- % Cloruro en hormign
m V1 %Cl
1.0031 4.00 0.121517
1.0028 4.05 0.112203
0.9993 4.03 0.117288
0.9985 4.05 0.112686
0.9976 4.05 0.112788
Media: 0.115
: Desv. Est. 4.04E-03
ualeat = 4.04 10
-3
g
utot =4.6 10
-3
g
Caso 3
En un trabajo de Ramachandran y Rashmi (9) sobre la estimacin de incertidumbres para un mtodo
espectrofotomtrico, se puede tambin observar que la componente principal en el valor de la incertidumbre
total es la que proviene de las medidas experimentales de la absorbancia., es decir, que la componente
aleatoria es, una vez ms, la ms importante.
CONCLUSIONES
Para el clculo de la incertidumbre de un ensayo, se desglosa el mismo en etapas elementales, se realizan
las estimaciones de las incertidumbres para cada una de ellas y luego se introducen todos estos trminos en la
frmula de propagacin. Es all donde se puede determinar la incidencia relativa sobre el resultado final de
cada una de las componentes de incertidumbre involucradas, lo que en general lleva a despreciar la mayora
de ellas, siendo muy pocas las dominantes (en la mayora de los casos 1, 2 o 3 componentes).
En los casos citados, las componentes predominantes de incertidumbre son las de Tipo A (componentes
aleatorias), por lo cual la incertidumbre estndar combinada se reduce a la desviacin estndar del promedio
de los resultados experimentales.
BIBLIOGRAFA/REFERENCIAS
1. UNIT-ISO-IEC Norma 17025 (2000) Requisitos generales para la competencia de laboratorios de ensayo y
calibracin UNIT-Uruguay
2. Centro Espaol de Metrologa (1998) Gua para la expresin de la Incertidumbre de Medida 1 Edicin,
Ministerio de Fomento, Espaa
3. Pezet Sandoval, F.- Mendoza Illescas,J. (1996) Vocabulario Internacional de trminos fundamentales y generales
de Metrologa Centro Nacional de Metrologa, Mxico
4. ISO/TAG4/WG3 (1992) Guide to the expression of uncertainty in
measurement (GUM) ISO-IEC-BIPM Ginebra.
Suiza
5. Eurachem-Citac (2000) Quantifying uncertainty in analytical measurement EURACHEM. UK.
6. Derrgibus, M. Fuentes, J.- Hermida, S.-Vazquez, M. (2000) Acreditacin de laboratorios. Validacin de
un procedimiento analtico II Encuentro de Ingeniera Qumica. AIQ, Uruguay
7. Derrgibus, M.T.- Fuentes, J. (2001) Cloruros en hormigones: Ensayos Interlaboratorios IV Simposio AOAC.
Uruguay
8. ISO 4787 (1984) Laboratory glassware- Volumetric glassware- Methods for use and testing of capacity. ISO.
Suiza.
9. Ramachandran, R- Rashmi (1999) Uncertainty of measurements in spectrophotometric analysis: a case study
Analyst, 124: 1099-1103
10