TFG Rodrigo Sanchez Vera
TFG Rodrigo Sanchez Vera
TFG Rodrigo Sanchez Vera
El agua es el recurso natural más valioso sobre la Tierra. Es omnipresente en todos los
ámbitos de la vida y transversal. En la industria química no iba a ser menos, empleándola
como materia prima, como disolvente o para el transporte de calor entre otros. Por ello
es esencial controlar ciertos parámetros para la correcta utilización. Uno de esos
parámetros es la dureza del agua. Es necesario controlarla ya que si no puede provocar
incrustaciones en equipos tales como intercambiadores, calderas, bombas o tuberías y
reducir su rendimiento o dañarlos seriamente.
La dureza del agua es la presencia de iones metálicos, principalmente Ca2+ y Mg2+, en el
agua bajo forma de cloruros y sulfatos (responsables de la dureza permanente) y
bicarbonatos (responsables de la dureza temporal). Estos precipitan y sedimentan formando
costras de carbonatos, óxidos y sulfatos que es la que provoca daños en los equipos como
se ha explicado anteriormente.
Se han desarrollado un gran número de métodos de análisis para determinar la dureza del
agua en las muestras. El método que se ha empleado en este TFG es el de la valoración
complexométrica con AEDT. De esta manera el ión metálicos forma un complejo con el
ácido etilendiaminotetraacético o AEDT tal y como se ve en la figura 1. Se lleva a cabo la
valoración a pH 10 con Negro de Eriocromo T o NET como indicador. En presencia de iones
metálicos la disolución presentará un color rosado y tornará a azul cuando ya todos los
iones hayan reaccionado con el AEDT.
A raíz de los problemas que provoca la dureza del agua se han ideado diversos procesos de
ablandamiento o eliminación de los iones calcio y magnesio. El que se estudia en este TFG
es el del uso de una resina de intercambio iónico. En estas resinas, formada por esferas de
polímeros compuestos de vinilbencenos y divinilbencenos, se introducen grupos activos que
luego llevan a cabo intercambios iónicos con los cationes de la disolución que se quiere
tratar. En este caso la resina contiene H+ y los intercambia por los iones Ca2+ presentes en
el agua tal y como se ve en la reacción siguiente:
2 R – H+ + Ca2+ R2 – Ca2+ + 2H+
A medida que se lleva a cabo el intercambio la resina sufre un proceso llamado saturación,
que significa que no puede intercambiar más iones porque están ocupados todos los puntos
activos. Cuando la columna se satura se debe llevar a cabo una regeneración. Esta
Una vez puesto a punto todo el sistema se ha procedido a optimizarlo. Se han realizado
diversos ensayos midiendo el pH automáticamente y tomando muestras de la dureza del
agua a la salida para valorarlas manualmente y comprobar si el sistema funcionaba. Se ha
Los costes totales del este trabajo se estiman en 7.822,13 €, siendo la partida de personal el
coste más elevado asociado a este proyecto.
Se considera que los resultados obtenidos son muy satisfactorios y cumplen con los
objetivos planteados inicialmente, aunque se proponen también ciertas líneas de mejora,
como puede ser el control de la alimentación o la inserción del proceso de regeneración en
la automatización.
Palabras clave:
Códigos UNESCO
230331 Química del agua
330308 Desionización
330806 Regeneración del agua
I) INTRODUCCIÓN Y OBJETIVOS
“El agua es la fuerza de la naturaleza”, dijo Leonardo Da Vinci. Presente en todos los
aspectos de la vida, el agua posee un papel protagonista en nuestra sociedad. Ya sea
en el ámbito de la salud, de la agricultura, de la ganadería o de muchos otros, el agua
tiene un impacto indirecto en nuestro día a día mucho más allá del simple consumo
Por supuesto el sector de la industria química no se queda atrás, llegando a emplear el
11% del agua dulce en Europa en 2016 [1] ya sea como materia prima, como disolvente
o para el transporte de calor entre otros. Sin embargo, hay que controlar ciertos
parámetros específicos en este sector para asegurar el buen funcionamiento de los
procesos. Entre estos parámetros se encuentra la dureza del agua, que provoca
calcificación en los equipos tales como las bombas, los calentadores, las calderas, los
intercambiadores, o las mismas tuberías.
Si bien no hay un límite legal de dureza del agua en España para el agua de uso
industrial, se llevan a cabo diversos procesos de ablandamiento del agua para
conseguir valores recomendables u óptimos. Un método de eliminación de la dureza
contrastado es el intercambio iónico. Este proceso se basa en el uso de una resina de
polímeros con grupos activos, en el caso de este trabajo fin de grado grupos activos
sulfónicos que presentan protones que se intercambian por los iones Calcio o
Magnesio, que se quedan fijados en la resina, pasando esos H + al efluente y
permitiendo medir el pH,.
El reto que presenta este proceso es la saturación de la columna. Esto significa que
cuando la resina ha intercambiado la gran mayoría de los iones de sus grupos activos,
el rendimiento del intercambio decrece de forma drástica. Por ello es necesario saber
en qué momento se llega a ese punto para controlar de forma óptima el proceso de
intercambio.
El objetivo de este trabajo fin de grado es el de diseñar un sistema de control para el
intercambio iónico con un analizador de bajo coste. Para ello hay diversos objetivos
parciales:
- Demostrar la relación entre las curvas de ruptura del pH y la de la dureza
- Diseñar un equipo donde introducir una sonda que esté sumergida en todo
momento y que permita medir en continuo valores del pH a través de un
analizador que funcione con Arduino
- Diseñar una transcripción de las medidas al ordenador
- Minimizar el coste total del sistema y su complejidad.
1) Definición
2) Métodos de análisis
A lo largo de los años se han desarrollado un gran número de métodos para determinar la
dureza del agua. Estos métodos se basan en las propiedades físicas y químicas de los
cationes Ca2+ y Mg2+. Los principales métodos conocidos se pueden dividir en dos grandes
categorías según lo que se mida. Por un lado están los métodos físicos (conductividad o
espectroscopia) y los métodos químicos (gravimétrico, valoración complexométrica con
AEDT y métodos de Hehner y Blacher).
El método de la conductividad se basa en la conducción eléctrica a través de la muestra
entre dos puntos de diferente potencial y sigue la ley de Ohm. Se puede suponer una
equivalencia de 2 µS/cm con 1 mg de CaCO3. Sin embargo esta equivalencia depende de
diversos parámetros, como la temperatura, que a su vez condiciona la viscosidad, de la
distancia entre electrodos y de la presencia de otros iones en la muestra. También hay que
tener en cuenta que la relación no es proporcional en todo momento ya que a muy bajas y
muy altas concentraciones hay interacciones entre los iones que afectan a la conductividad.
Además no es un método selectivo [3].
(1)
El método de Hehner es una valoración con un ácido fuerte, en concreto HCl, y permite
calcular la dureza permanente como la temporal basándose en la alcalinidad de la muestra y
el hecho de que la dureza temporal se pueda eliminar mediante ebullición. Se valora por un
lado la muestra con HCl (generalmente 0,1 M) y en presencia de naranja de metilo y se
obtiene un volumen Va. Posteriormente se toma una muestra de mismo volumen, se lleva a
ebullición y una vez evaporado el líquido se diluye en agua destilada y luego se filtra. Se
repite la valoración obteniendo un volumen Vb. La diferencia de volúmenes entre las dos
valoraciones se asocia a la cantidad de carbonatos de calcio y de magnesio, o lo que es lo
mismo, a la dureza temporal. El mismo proceso se repite para la dureza permanente pero
añadiéndole carbonato sódico. El método de Blacher es una volumetría y se basa en la
reacción del jabón con los elementos responsables de la dureza. Se emplea el palmitato
potásico (un tipo de jabón de fórmula CH3(CH2)14COOH que provoca la precipitación de los
cationes. Cuando todos han precipitado, el añadir palmitato provoca una variación muy
marcada de pH que se puede detectar con fenolftaleína. Este método se basa en el antiguo
método de Clark, ya obsoleto, que consistía en añadir jabón y agitar hasta que se cree una
espuma permanente, lo que indica que se ha llegado al punto de equivalencia. [3]
3) Métodos de ablandamiento
1) Intercambio iónico
Las resinas se pueden dividir en cuatro grupos en función de sus grupos activos y de los
iones que intercambian: las resinas débilmente básicas (WBA), las resinas débilmente
ácidas (SBA), las resinas fuertemente básicas (WAC) y las resinas fuertemente ácidas
(SAC). Son estas últimas las que se usan para el tratamiento del agua y la que se ha usado
para este trabajo [6].
La estructura de una resina es polimérica. Son pequeñas esferas compuestas por un
polímero cuyo monómero principal es el vinilbenceno. También se le añaden monómeros de
divinilbenceno para entrecruzar las cadenas y conseguir una estructura tridimensional. En la
figura 2 se puede ver el polímero compuesto de vinilbencenos y, en naranja, un enlace entre
las cadenas provocado por la inserción de un monómero de divinilbenceno. La
polimerización se lleva a cabo generalmente en un reactor agitado [7].
C) Arduino.
Por un lado se ha preparado el lecho de resina. Este lecho se ha colocado en una columna
de vidrio para intercambio iónico que previamente ha sido limpiada y se ha verificado su
buen funcionamiento, en especial el de la llave de salida. En esta columna hay una lámina
de frita porosa de vidrio cerca de la llave que permitirá que la resina no sea arrastrada por el
agua.
La resina es de tipo ácida o catiónica. Cabe destacar que cuando se adquirió esta resina fue
bajo forma sódica y hubo que tratarla, en el marco de otro proyecto, con HCl para obtener la
resina de la que hoy disponemos. Puesto que se desconocía el último uso que se dio a la
resina, se decidió llevar a cabo una regeneración y un lavado de ésta para dejarla preparada
para su correcto uso.
La regeneración de la columna se ha llevado a cabo vertiendo en ella 60 mL de la disolución
de HCl 3 M. Al regenerar la columna después de cada ensayo, se ha comprobado cómo
aumenta el volumen del lecho y recupera su color amarillo claro en vez del naranja que
había adquirido. Una vez regenerada hay que lavar la columna con agua destilada para
eliminar los cloruros que aparecen al regenerarla. Es importante eliminarlos para que no
afecte al pH y, por ende, a los resultados. Es en la fase de lavado en la que la columna
retoma su volumen inicial, que es el más grande de todo el proceso.
Para comprobar que se han eliminado los cloruros, se recoge una muestra de agua a la
salida de la columna y se le añaden unas gotas de la disolución de nitrato de plata 0,1 M. En
presencia de cloruros la muestra torna a un color blanquecino, como se puede ver en la
imagen 1 debido a la aparición de un precipitado de cloruro de plata, y en ausencia de ellos
la muestra se mantiene incolora, como se ve en la imagen 2. Una vez se ha obtenido una
muestra sin cloruros se ha de dejar la columna siempre con una cantidad de agua por
encima de la resina para que no se deshidrate. El proceso de regeneración se ha de repetir
después de cada ciclo de intercambio o ensayo.
Para ello se ha hecho un primer ensayo con una alimentación que contiene 5 g de
CaCl22H2O disueltos en un litro de agua, lo que equivale a una dureza de 3.390 mg de
CaCO3/L. Se ha llevado a cabo el ensayo manteniendo el caudal fijo en 65 mL/min dentro de
lo posible y se han recogido en viales de vidrio muestras de 15 mL cada minuto.
Posteriormente se han medido el pH y la dureza de cada una de las muestras por duplicado.
La medición de pH se ha realizado empleando un pH-metro del laboratorio previamente
calibrado, modelo ORP PCE-228. La dureza del agua se ha medido mediante valoración
con AEDT. Para la valoración se han pipeteado 5 mL de la muestra y se le han añadido 2
mL de la disolución tampón NH4Cl/NH3, para mantener el pH estable, y unas gotas de NET.
Este proceso se ha llevado a cabo dos veces por cada muestra y, en caso de haber poca
concordancia entre los dos resultados, se ha hecho una tercera vez, ya que el volumen de la
muestra, 15 mL, se ha elegido de manera que permita un tercer ensayo de comprobación.
En el primer ensayo de intercambio llevado a cabo se obtuvieron resultados poco
clarificadores, como se ve en la gráfica 1, ya que el punto de ruptura aparece demasiado
pronto y no permite comprobar bien la relación entre la curva del pH y la de la dureza. Esto
es debido a que la columna se ha saturado demasiado rápido por la alta concentración de
Ca2+. En esta gráfica se representan el pH de la muestra y la relación de concentraciones,
es decir, la concentración de Ca2+ de la muestra (C) dividida entre la concentración de Ca2+
de la alimentación (Co).
Gráfica 1: Ensayo 1
Sin embargo, se ha podido observar la evolución del lecho a lo largo del proceso. Si bien el
lecho ha comenzado siendo de un color amarillo claro, poco a poco se ha ido compactando
y se ha podido comprobar cómo el color ha ido cambiando a naranja verticalmente de
manera homogénea. Este fenómeno se puede observar en la imagen 3. Esta zona
compacta y anaranjada corresponde a la parte de la resina que está saturada.
Se ha hecho un segundo ensayo de intercambio con una dureza del agua inferior, más
específicamente de 870 mg de CaCO3/L. El proceso de medida de pH y de la dureza ha sido
el mismo que para el ensayo 1. Los resultados se pueden ver en la gráfica 2 y se observa
que el punto de ruptura del pH y de la dureza del agua ocurren a la vez, lo cual corrobora la
hipótesis de que midiendo la variación del pH durante el proceso de ablandamiento
podemos indirectamente saber la evolución de la dureza.
Se han realizado dos ensayos con la misma dureza del ensayo 2 pero con diferentes
caudales de salida, como se ha explicado anteriormente. El ensayo 2.A se ha realizado con
un caudal de 65 mL/min y el ensayo 2.B con un caudal de 47 mL/min. Conociendo los
caudales de salida y el volumen total, se ha calculado el intervalo necesario entre cada
muestra para conseguir un número idéntico de muestras (15) y así poder comparar si son
proporcionales, es decir, que para la muestra 8 de ambos ensayos habrá atravesado la
columna el mismo volumen de agua, 500 mL. Los resultados se pueden ver en la gráfica 3 y
se observa que la evolución de la relación de durezas es prácticamente idéntica, por lo que
el caudal de salida no afecta a la capacidad de intercambio, al menos en los rangos
estudiados.
Como se ha explicado anteriormente, el objetivo final del TFG es diseñar un equipo con el
que poder controlar la saturación de la columna en tiempo real y maximizando el
rendimiento. Para ello se han estudiado opciones muy diversas que se exponen en este
apartado.
Uno de los grandes retos del diseño ha sido asegurarse que el sensor de pH estuviese en
todo momento sumergido y que pudiese hacer medidas en continuo. El primer parámetro es
una exigencia propia del sensor para el buen funcionamiento propio. El segundo tiene por
objetivo mejorar la toma de medidas. En TFGs anteriores se diseñaron instalaciones con un
depósito al que caía el agua, se esperaba a que se llenase lo suficiente para poder tomar la
medida y se vaciaba mediante una válvula en cada toma de medida. El problema de este
proceso es que las medidas en realidad no son en continuo y, por lo tanto, la precisión con
la que se detecta el momento en el que se satura la columna no es óptima.
De cara a solucionar este problema una primera opción fue cambiar la disposición de la
columna de vertical a horizontal. Esta solución solventaría el problema de la inmersión de la
sonda a la par que haría disminuir el trabajo de la bomba, ya que sólo habría que impulsar el
agua y no habría una diferencia de altura que superar como hay en el diseño vertical.
También se pensó que al estar en vertical, la regeneración a contracorriente sería más
sencilla.
Sin embargo, esta solución presenta diferentes problemas. El principal es la disposición del
lecho. De manera similar al montaje inicial, se dispondría de una lámina de frita porosa de
vidrio a un lado para que no se arrastrase la resina con el agua y otra lámina fritada al lado
opuesto para que no se arrastrase con el regenerador. Esta segunda lámina, al ser fija,
habría que posicionarla de manera que cupiera el lecho cuando su volumen es máximo
(después del lavado). El problema surge cuando, durante la operación, el volumen del lecho
se reduce. Al estar en posición horizontal, el lecho se quedaría en la parte inferior de la
columna dejando todo un espacio arriba por donde el agua podría pasar sin ser tratada
(figura 9). Además, se aumentaría la pérdida de carga al añadir una segunda frita porosa.
Una solución a este problema consistiría en poner una lámina móvil que se adaptase al
volumen del lecho. Esto es problemático por un lado porque al estar en horizontal haría falta
algún tipo de muelle para la adaptación (con su consiguiente pérdida de carga), y por otro
por como asegurar el deslizamiento de la lámina sin que haya escapes ni se descompacte el
lecho. Por estas razones se desechó la opción de la horizontalidad.
La siguiente opción que se ha analizado ha sido la de añadir un sifón en la salida de la
columna. Este tipo de dispositivos, muy comunes en fontanería, permiten mantener un
volumen de agua en todo momento mediante vasos comunicantes y diferencia de altura. De
esta manera se podría introducir la sonda de pH y medir en continuo.
Analizando los pros y los contras del sifón surgió la duda de si aparecían volúmenes
muertos en los codos y/o caminos preferentes en la curva interior como se ilustra en la figura
10. Sin embargo no se ha encontrado ninguna bibliografía ni información acerca de ello.
Una vez se tienen los dos sifones hay que decidir dónde se introduce la sonda de pH. El
parámetro que se ha considerado más importante es el de que se pueda retirar fácilmente,
ya que la sonda hay que calibrarla cada día. Sin embargo, el que sea fácilmente removible
dificulta la estanqueidad del equipo. Por ello se ha decidido que se va a perforar en el codo
3 para el sifón negro y su equivalente en el sifón blanco. De esta manera la perforación no
está sumergida pero la sonda se puede introducir hasta el fondo del sifón.
Esto representa sin embargo un problema puesto que se está midiendo el pH a la salida del
sifón y no a la salida de la columna. Debido al problema de la estanqueidad no se ha podido
introducir la sonda antes y por lo tanto hay cierto volumen de retardo. Se tendrá en cuenta
en las líneas futuras.
D) Diseño analizador.
La escritura del código Arduino presenta tres secciones diferentes. La primera sección
permite iniciar variables y crear funciones propias. La segunda sección, llamada void setup,
se ejecuta una sola vez y la tercera, llamada void loop, se ejecuta en bucle. En Arduino se
escribe el código en la sección void loop en general por cuestiones prácticas y es lo que se
ha hecho también en este proyecto.
El primer código que se ha escrito es el de calibración del sensor de pH, aspecto
especialmente importante para obtener medidas fiables. Este sensor, que va conectado a la
placa Arduino por medio del módulo de pH, lee el voltaje de la disolución gracias a un
electrodo. La relación entre este voltaje y el pH es lineal como se puede ver en la ecuación
siguiente:
Una vez obtenidos estos valores se pueden conseguir transformar los voltajes en pH. Si bien
se podría hacer con un código con el Arduino, se ha decidido que este cálculo se va a llevar
a cabo en el Excel para así aligerar el código del Arduino.
Se ha escrito un código para el proceso completo. El código completo aparece en la figura
15 y el esquema explicativo en la figura 14.
Se ha puesto también una alarma para que suene al alcanzar el punto de ruptura. La idea
inicial era poner un valor alarma, pero se descartó por dos razones. Por un lado el valor
variaba en función de la calibración y por otro los valores leídos al principio del proceso
provocaban que sonase antes de lo deseado. Se ha decidido entonces que la alarma suene
cuando haya un cambio de tendencia.
Para ello se han introducido unas funciones para comparar los valores leídos con la media
de los últimos 100 valores. Para que funcione bien se han empezado a guardar los valores
en el momento en el que el agua entra en contacto con la sonda. Si se detecta una
diferencia del valor con respecto a la media se tiene en cuenta en un contador. Este
contador se pone a cero si el valor siguiente no presenta esa diferencia con la media, y de
esa manera se evita que suene la alarma con errores de lectura. Al llegar a 3 el contador
confirma que es un cambio de tendencia, ya sea por ruptura de la curva o por fallo en el
proceso, y suena la alarma.
En paralelo se escriben todos los valores leídos por la sonda en el ordenador. Este código
es válido para este TFG cuyo objetivo es comprobar que el proceso funciona correctamente.
Para una aplicación industrial se puede hacer que el Arduino deje de enviar información
cuando suene la alarma.
La toma de muestras durante todo el proceso es continua y se registra una gran cantidad de
información. Aprovechando que el sistema está automatizado y que Arduino está conectado
al ordenador sería muy positivo poder transcribir esas medidas puesto que el registro de
éstas es muy interesante a muchos niveles, permitiendo analizar el funcionamiento del
proceso y abrir vías de mejora. En otros trabajos similares de fin de grado se consiguió
comunicar el Arduino con MATLAB y así obtener estos datos en tiempo real.
Sin embargo se pensó que sería interesante no solo tener los datos de los pH o de los
tiempos del ensayo sino también poder guardarlos y comparar entre las diferentes
repeticiones, permitiendo así analizar tendencias, detectar errores, mejorar la
automatización e incluso extrapolar estos datos para otros experimentos.
Para llevar a cabo esta idea se consideró la opción de utilizar Microsoft Excel. Esta
aplicación ofrece la opción de tratar una gran cantidad de datos con una interfaz operativa
muy completa y a la vez accesible a todo el mundo. Además cuenta con un servicio de
ayuda muy completo y foros de usuarios muy activos. Una de las grandes ventajas del
programa de cara a este TFG es que se pueden guardar los datos de cada ensayo en hojas
diferentes de un mismo libro y luego se pueden relacionar estos datos entre ellos,
permitiendo así un análisis más profundo como se explica en el párrafo anterior.
El primer paso ha sido el de conectar Excel con el Arduino. Para ello hay que programarlo
para que abra el puerto serie y lea los datos que recibe el ordenador por ahí. Y para
conseguir todo este proceso es necesario usar una macro de Excel o VBA.
Una macro o VBA es un conjunto de acciones que se pueden automatizar en Excel de
manera que se puedan repetir tantas veces como sea necesario. Las macros también
permiten acceder a ciertas funcionalidades del Excel más complejas mediante la
programación de ciertas funciones y rutinas. Estas funciones se recogen en lo que son
conocido como librerías y todo el mundo puede crear su propia librería. Existen también
librerías públicas.
Para este TFG se ha buscado en internet y se ha encontrado una librería que permite
conectar el Excel con el puerto serie. Esta librería fue creada por David M. Hitchner y es de
libre acceso en internet. Si bien ofrece una gran cantidad de funciones muy útiles y bien
desarrolladas, debido al nivel de conocimiento de VBA del autor de este TFG se ha decidido
limitarse a una programación básica y al uso de las funciones más simples.
El primer paso para la transcripción de las medidas ha sido el de plantear cómo se quieren
visualizar y que elementos han de aparecer.
Por un lado es importante que aparezca una tabla con los valores de la calibración de la
sonda del pH ya que, como se ha explicado en el apartado anterior, la lectura que recibe el
Excel es la del voltaje y ha de traducirse a pH mediante una fórmula con los valores de la
calibración. Esta tabla se pone al principio y consta de cuatro valores como se puede ver en
la tabla 2. Al introducir los valores de voltaje correspondientes a pH 7 y pH 4, se ha
programado el Excel de manera que calcule directamente la pendiente y la ordenada en el
origen según (2) y (3).
pH 7
pH 4
Pendiente
Ordenada
en el
origen
(2)
(3)
En lo respectivo a los datos en sí, se ha decidido tomar tres medidas en una misma abscisa.
Es decir, que si el intervalo de lectura del voltaje es de i ms, se considerará que se mide
cada 3 x i ms y en cada intervalo habrá tres lecturas. Esto se traduce en el Excel en una
tabla de seis columnas para cada intervalo de tiempo como se ve en la tabla 3.
De esta manera en los voltajes de entrada se van escribiendo las lecturas de la sonda de pH
y el Excel calcula automáticamente los pH de entrada.
Para la escritura automática de los voltajes se han empleado tres macros. En un primer
lugar una macro para abrir el puerto serie, es decir, que sea Excel el que tiene acceso a este
puerto lo que desinhibe a otros programas (por ejemplo el monitor serie del Arduino). Esta
macro, llamada Abrir Puerto, se le asigna a un botón de Excel y su esquema se puede ver
en la figura 16. El número del puerto hay que introducirlo manualmente en el código. En
caso de haber algún tipo de error se ha hecho que aparezca un mensaje indicándolo. Este
error puede ser debido a que ya esté abierto el puerto o a que otro programa lo esté usando.
Finalmente se escribe una última macro llamada Cerrar Puerto y cuya función es la de
liberar el puerto. También funciona con botón y se puede ver el código en la figura 18.
Gráfica 4: Ensayo 3
Otro fenómeno que se observa es que el pH comienza siendo muy alto hasta el minuto 1
aproximadamente y luego desciende bruscamente. El hecho de que la sonda no esté
sumergida al principio provoca la medida del pH de 4,25. Y en cuanto el agua llega a la
sonda ésta empieza a medir el pH del líquido. Esta evolución viene ilustrada en la figura 20.
Sin embargo, se puede ver en el gráfico 4 que la curva de ruptura del pH va con retraso con
respecto a la curva de ruptura de la relación de concentraciones. Esto conlleva que la
dureza aumenta en el agua de salida antes siquiera de que varíe el pH, lo cual imposibilita
cualquier control por medio del análisis del pH, al menos en esta disposición.
Uno de los fenómenos que podría provocar este desajuste es el de los volúmenes muertos
en el sifón. Como se explica en el apartado III-C y se ve en la figura X, es posible que dentro
del sifón se creen volúmenes muertos y caminos preferentes y que uno de esos volúmenes
muertos sea en el que está la sonda de pH. Se ha pensado que se puede solucionar
colocando la sonda más arriba, a la altura del nivel del agua cuando ya está saliendo. En la
figura 21 se describen las dos posiciones del pHímetro, siendo la posición A en la que se ha
hecho el ensayo 3 y la posición B en la que se ha hecho el ensayo 4.
Gráfica 6: Ensayo 4
En la gráfica 7 se puede ver la calibración de la bomba peristática entre 20 y 120 rpm, que
son los valores límites en los que funciona la bomba. Sin embargo se ha observado que hay
una condición que no se cumple por debajo de las 120 rpm y es que no se mantiene una
capa de agua por encima de la resina, aspecto importante para el correcto funcionamiento
de la resina. Esto es debido a que el caudal de salida de la columna únicamente y con una
mínima columna de agua sobre el lecho es superior a la alcanzada por la bomba peristática.
Esto provoca problemas a nivel del funcionamiento de la columna puesto que se crean
caminos preferentes. El hecho de tener una capa de líquido provoca que el agua atraviese
todo el lecho y de manera homogénea ya que hay una misma presión sobre toda la sección
circular superior. Sin embargo, si esa capa de líquido no existe, entonces el chorro de la
alimentación cae sobre un punto fijo y se crea un camino preferente a través de la resina
que, consecuentemente, dará como resultado que se saturará antes, cuando en realidad
una parte no ha sido utilizada.
Esta misma calibración se repitió variando la altura que debía superar la bomba para
alcanzar la entrada de la columna con el objetivo de ver si afectaba a los caudales. Se
observó que no había ninguna variación ya que se obtuvieron exactamente los mismos
resultados.
Una vez tenemos hecha la calibración de la bomba, se hace el montaje completo del
sistema. Como se ve en la figura 22, se conectan el tubo de entrada de la bomba peristática
(1) directamente a la alimentación y el de salida se conecta a la entrada de la columna de
intercambio iónico (2). El agua saliente de la columna va directamente al sifón (3), a la salida
del cual se recoge el agua tratada, al igual que las muestras. En el sifón se encuentra la
Una vez hecho el montaje final, se ha llevado a cabo el ensayo 5, utilizando la bomba
peristática a 120 rpm y con una alimentación de 430 mg de CaCO3/L. Los resultados se
pueden observar en la gráfica 8.
Gráfica 8: Ensayo 5
Durante este ensayo se ha producido un problema a nivel del lecho. Al empezar a saturarse,
el lecho reduce su volumen, lo que a su vez reduce la pérdida de carga del lecho y provoca
que esa capa de agua que ha de mantenerse por encima de la resina no se produzca. Esto
provoca caminos preferentes y que la columna se sature de forma desigual como se puede
ver en la imagen 9, dónde la parte clara del lecho es resina sin saturar y la parte más oscura
es la resina saturada.
Las consecuencias de la creación de estos caminos preferentes se pueden ver en la gráfica
X. Hay subidas seguidas de pequeños valles en los minutos 18 y 23. Esto significa que se
han saturado ciertas partes de la resina mientras que otras aún no lo están. Se ha decidido
no realizar la valoración de las muestras puesto que se ha considerado este ensayo como
fallido.
Para solucionar este problema hay que actuar sobre los caudales del sistema. En lo relativo
a la bomba peristática no se puede aumentar el caudal más de lo que ya está, como se
indicó anteriormente. El único caudal que queda sobre el que se puede actuar es sobre el de
salida de la columna. En efecto la columna consta de una llave manual para la salida de
líquidos. En todos los ensayos se ha operado con la llave abierta al 100% puesto que, al ser
manual, es la única manera de tener un caudal replicable. Se ha decidido entonces cerrar
ligeramente la llave de salida para disminuir el caudal de salida de la columna. En todos los
ensayos que siguen se ha ajustado el caudal de la bomba de tal manera que quede siempre
esa lámina de fluido sobre la resina. Eso implica que, durante la operación, cuando la resina
se va saturando y por ende disminuyendo su volumen, se ha de aumentar manualmente el
caudal de la bomba para mantener siempre el mismo nivel de agua y manteniendo así
constante el caudal de salida de la columna.
Se ha realizado el ensayo 6 con una alimentación de 340 mg de CaCO3/l y con las
condiciones descritas anteriormente. Los datos se representan en la gráfica 9.
Gráfica 9: Ensayo 6
Podemos ver que en este ensayo los resultados son muy positivos ya que ambas curvas
rompen en el mismo punto y siguen una progresión casi idéntica. También podemos
observar que la duración del ensayo es mucho mayor debido a la disminución del caudal.
Por esta razón la frecuencia de recogida de muestras ha aumentado de un minuto a tres
minutos. Puesto que los resultados han sido tan positivos se ha decidido hacer un ensayo 7
en las mismas condiciones para confirmar y evaluar la reproducibilidad de la instalación y
del sensor. Se realizó con una dureza de 320 mg de CaCO3/L.
Como se puede ver en la gráfica 10, los resultados del ensayo 7 son igual de positivos, lo
cual confirma que nuestro sistema de control de la dureza es válido y replicable. Se puede
observar también como, al disminuir la dureza, la curva se desplaza a la derecha y tarda un
poco más en saturarse.
Vistos los buenos resultados obtenidos, se ha decidido probar el mismo sistema con el sifón
blanco. Se ha realizado un ensayo 8 con una alimentación de 490 mg de CaCO3/L.
En la gráfica 11 se pueden ver como los resultados son bastante malos. Probablemente sea
debido al volumen muerto presente justo en la parte central inferior del sifón, que también
provoca variaciones en la mecánica del fluido en el interior de éste. Se ve claramente en la
curva de relación de concentraciones como el final de la ruptura no es tan claro como en los
otros ensayos, lo cual indica volúmenes muertos. Se ha concluido que la forma del sifón
blanco no permite un buen funcionamiento del sistema de control y, por lo tanto, no se harán
más ensayos con este sifón.
Dureza de la
Intervalo de Tiempo Momento pH del
Caudal alimentación
Ensayo Fecha Bomba Sifón Exitoso tiempo entre total del del punto punto de Caudal (l/h)
(ml/min) (mg de
datos (ms) ensayo de rotura rotura
CaCO3/l)
3 18-07-19 NO Negro SI 500 20,50 11,975 1,49 48,78 2,927 460,000
4 22-07-19 NO Negro SI 500 20,65 8,20 1,27 48,43 2,906 430,000
5 24-07-19 SI Negro NO 500 29,80 17,33 2,05 33,56 2,013 430,000
6 25-07-19 SI Negro SI 500 96,00 65 1,192 10,417 0,625 340,000
7 25-07-19 SI Negro SI 500 95,38 61,5 1,135 10,485 0,629 320,000
8 26-07-19 SI Blanco NO 500 58,35 27,73 1,646 17,138 1,028 490,000
V. PLANIFICACIÓN TEMPORAL
VI. COSTES
En un primer lugar se han calculado los costes del material empleado como se puede
ver en la tabla 6. Los costes ascienden a 404,41 €.
Equipo Suministrador Unidades Coste por unidad (€) Coste total (€)
Base Labbox 3 13,41 40,23
Placa Arduino UNO Ebay 1 7,85 7,85
Protoboard Ebay 1 7,89 7,89
Cables Arduino Ebay 1 1,00 1,00
Electrodo + pH Shield DFRobot 1 25,98 25,98
Zumbador Array Electrónica Profesional 1 1,35 1,35
Columna de vidrio Alamo 1 100,00 100,00
Aspirador manual para pipetas Selecta 2 5,06 10,12
Espátula cuchara plana Labbox 1 1,15 1,15
Marcador para laboratorio Edding 400 1 2,10 2,10
Matraz aforado 100 ml Labbox 5 2,71 13,55
Matraz aforado 1000 ml Labbox 1 7,81 7,81
Matraz aforado 200 ml Labbox 2 3,60 7,20
Matraz aforado 25 ml Labbox 32 2,38 76,16
Pipeta graduada de vidrio 2 ml Labbox 1 1,15 1,15
Pipeta graduada de vidrio 5 ml Alamo 1 2,97 2,97
Probeta graduada 100 ml Alamo 1 13,26 13,26
Sifón de fontanería Fontanería Venecia 1 5,00 5,00
Tubo de PVC de 25 mm, 2,5 m Leroy Merlin 1 1,65 1,65
Codo liso 25 mm con rosca 3/4' hembra Leroy Merlin 4 1,10 4,40
Pe macho 3/4' Leroy Merlin 2 1,70 3,40
IVA 21% 70,19
TOTAL 404,41
Por otro lado hay que contabilizar el coste de los reactivos. El desglose se ve en la tabla
7 y el coste total asciende a 48,02€.
También hay que tener en cuenta el coste de los equipos. Sin embargo para los equipos hay
que considerar la amortización y el tiempo que se han usado para obtener el coste real. Se
ha supuesto para todos los equipos una disponibilidad de 8 horas al día, 229 días al año, lo
que equivale a 1832 horas al año. Para calcular el coste real del equipo en el proyecto hay
que calcular sus horas de vida en función de sus años de vida útil y ver qué proporción de
esa vida útil se ha empleado en este proyecto. Los datos se presentan en la tabla 8 y el
coste total asciende a 16,11 €.
Equipo Modelo Precio (€) Vida útil (años) Utilización (h) Coste (€)
Balanza analítica Precisa XB 220A 1402 3 1,5 0,38
Bomba peristática Heidolph Pumpdrive 5201 3146 5 15 5,15
Cabezal bomba peristática Heidolph Pumpdrive 5201 508,2 3 15 1,39
Placa calefactora e imantada Janke & Kunel KS 501 D 590 8 18 0,72
Ordenador portátil Toshiba 6HHTDEE 248 4 250 8,46
TOTAL 16,11
Además hay que sumar el coste de personal, tanto del alumno como de los tutores. Los
datos se recogen en la tabla 9 y ascienden a 6642,5 €.
Para terminar hay que sumar todos los costes y añadirle los costes indirectos como son el
agua y la electricidad que se han considerado un 10% de los costes asociados al no tener
acceso al detalle de los gastos del laboratorio. El coste total del proyecto asciende entonces
a 7.822,13 €.
Como se ve la partida que más eleva el coste asociado a este proyecto es la de personal.
VII. CONCLUSIONES
1) Conclusones generales
A lo largo de este trabajo de fin de grado se han llevado a cabo diversos ensayos con
diferentes parámetros e innovando y arreglando problemas hasta que en el ensayo 6 y
el 7 se consigue un montaje con unos resultados muy correctos. Se ha conseguido
encontrar un espectro de durezas en el que trabajar y ver resultados con la cantidad de
resina disponible, se ha encontrado un método por el cual la sonda de pH está
sumergida en todo momento, se ha implementado un sistema por el cual hay
mediciones en continuo y las lecturas se llevan en tiempo real al ordenador dónde son
analizadas y comparadas con otros ensayos, todo con un montaje sencillo, eficiente y
económico.
VIII. BIBLIOGRAFÍA
[3] Neira Gutiérrez, Marco Antonio. “Dureza en aguas de consumo humano y uso
industrial, impactos y medidas de mitigación. Estudio de caso: Chile”,
Departamento de Ingeniería Civil, Facultad de Ciencias Físicas y Matemáticas.
Universidad de Chile, 2006.
[6] Cortijo Herrera, Dora. “Desalcalinización del agua mediante intercambio iónico”,
Universidad de Lima, Perú.
[8] Nomura Micro. “Cation detector for detecting cation leakage in water treatment
plants compromises measuring column filled with ion exchange resin and pH
meter to mmeasure pH of water entering and leaving measuring column”.
JP2001056323.
[11] Ruiz Gutiérrez, Jose Manuel. “Arduino: Manual de programación”. Disponible en:
https://arduinobot.pbworks.com/f/Manual+Programacion+Arduino.pdf
A) Índice de tablas
B) Índice de figuras
C) Índice de gráficas
- Gráfica 1: Ensayo 1
- Gráfica 2: Ensayo 2.A
- Gráfica 3: Comparación entre los ensayos 2.A y 2.B
- Gráfica 4: Ensayo 3
- Gráfica 5: Detalle del ensayo 3 entre el minuto 2 y el minuto 4
- Gráfica 6: Ensayo 4
- Gráfica 7: Calibración de la bomba peristática
- Gráfica 8: Ensayo 5
- Gráfica 9: Ensayo 6
- Gráfica 10: Ensayo7
- Gráfica 11: Ensayo 8
D) Índice de imágenes
X. ANEXOS
En este apartado se muestran los datos de las diferentes valoraciones de las muestras de todos
los ensayos, así como datos técnicos de cada ensayo automatizado. Se ha decidido no incluir
los datos leídos por Arduino ya que son en algunos casos más de 11.500 datos por ensayo y se
ha considerado que no añadirían más que hojas a este TFG, teniendo en cuenta que esos datos
se pueden leer y analizar con las gráficas del apartado IV.
Ensayo 1:
Volúmen (mL) Moles
Muestra mg CaCO3/L C/Co pH
1 2 medio AEDT/Ca
1 4,3 3,8 4,05 0,0000405 810 0,23893805 0
2 16,9 16,6 16,75 0,0001675 3350 0,98820059 2,11
3 16,7 16,8 16,75 0,0001675 3350 0,98820059 4,61
4 16,6 16,6 16,6 0,000166 3320 0,97935103 5,11
5 16,6 16,7 16,65 0,0001665 3330 0,98230088 5,46
6 16,7 16,7 16,7 0,000167 3340 0,98525074 6,02
7 16,7 16,9 16,8 0,000168 3360 0,99115044 6,09
8 16,8 17 16,9 0,000169 3380 0,99705015 6,16
9 17 16,9 16,95 0,0001695 3390 1 5,98
10 16,9 16,8 16,85 0,0001685 3370 0,99410029 6,05
11 16,8 16,9 16,85 0,0001685 3370 0,99410029 6,11
12 16,7 17,2 16,95 0,0001695 3390 1 6,11
13 16,8 16,9 16,85 0,0001685 3370 0,99410029 6,09
14 16,9 16,9 16,9 0,000169 3380 0,99705015 6,07
15 16,7 16,8 16,75 0,0001675 3350 0,98820059 5,96
0 16,9 17 16,95 0,0001695 3390
Ensayo 2.A
Volúmen (mL) Moles
Muestra mg CaCO3/L C/Co pH
1 2 medio AEDT/Ca
1 0 0 0 0 0 0 1,6
2 0 0 0 0 0 0 1,6
3 0 0 0 0 0 0 1,6
4 0 0 0 0 0 0 1,6
5 0,3 0,2 0,25 0,0000025 50 0,05747126 1,6
6 0,5 0,4 0,45 0,0000045 90 0,10344828 1,64
7 1,3 1,1 1,2 0,000012 240 0,27586207 1,79
8 4,2 4,3 4,25 0,0000425 850 0,97701149 3,46
9 4,3 4,3 4,3 0,000043 860 0,98850575 3,89
10 4,2 4,3 4,25 0,0000425 850 0,97701149 4,09
11 4,3 4,3 4,3 0,000043 860 0,98850575 4,27
12 4,3 4,2 4,25 0,0000425 850 0,97701149 4,43
13 4,3 4,3 4,3 0,000043 860 0,98850575 4,56
14 4,3 4,2 4,25 0,0000425 850 0,97701149 4,77
15 4,3 4,3 4,3 0,000043 860 0,98850575 4,9
0 4,4 4,3 4,35 0,0000435 870
Ensayo 2.B
Volúmen (mL) Moles
Muestra mg CaCO3/L C/Co pH
1 2 medio AEDT/Ca
1 0 0 0 0 0 0 1,6
2 0 0 0 0 0 0 1,6
3 0 0 0 0 0 0 1,6
4 0 0 0 0 0 0 1,6
5 0,2 0,2 0,2 0,000002 40 0,04597701 1,6
6 0,5 0,6 0,55 0,0000055 110 0,12643678 1,64
7 1,4 1,3 1,35 0,0000135 270 0,31034483 1,79
8 4,2 4,2 4,2 0,000042 840 0,96551724 3,46
9 4,3 4,2 4,25 0,0000425 850 0,97701149 3,89
10 4,2 4,3 4,25 0,0000425 850 0,97701149 4,09
11 4,3 4,3 4,3 0,000043 860 0,98850575 4,27
12 4,3 4,2 4,25 0,0000425 850 0,97701149 4,43
13 4,3 4,3 4,3 0,000043 860 0,98850575 4,56
14 4,3 4,2 4,25 0,0000425 850 0,97701149 4,77
15 4,3 4,3 4,3 0,000043 860 0,98850575 4,9
0 4,4 4,3 4,35 0,0000435 870
Ensayo 3:
Intervalo de tiempo entre datos (ms) 500
Tiempo total del ensayo 20,5
Momento del punto de rotura 11,975
pH del punto de rotura 1,49
Caudal (ml/min) 48,78
Caudal (l/h) 2,927
Ensayo 4:
Intervalo de tiempo entre datos (ms) 500
Tiempo total del ensayo 20,65
Momento del punto de rotura 8,2
pH del punto de rotura 1,27
Caudal (ml/min) 48,43
Caudal (l/h) 2,906
Ensayo 5:
Intervalo de tiempo entre datos (ms) 500
Tiempo total del ensayo 29,8
Momento del punto de rotura 17,325
pH del punto de rotura 2,05
Caudal (ml/min) 33,56
Caudal (l/h) 2,013
Ensayo 6:
Ensayo 7:
Ensayo 8: