“Año del Fortalecimiento de la Soberanía Nacional”
UNIVERSIDAD DE HUANUCO
Facultad de Ingeniería
ESCUELA ACADEMICO PROFESIONAL DE INGENIERIA AMBIENTAL
MODELAMIENTO DE TRANSPORTE DE HIDROCARBUROS
EN SUELOS
DOCENTE: Vasquez Baca, Yasser
CURSO: Modelamiento y Simulación
INTEGRANTES:
Alvarez Velásquez, Jhomira Sheyla
Arias Bonilla, Brimson
Cabello Salas, Jordy Anibal
Camacho Armillón, Meier
Lozano Santiago, Manuel
Huánuco, 2023
MODELAMIENTO Y SIMULACIÓN ESPACIAL DE LA CONDUCTIVIDAD
O HIDRAULICA POR CONTAMINACIÓNDE HIDROCARBURS EN SUELOS DE
TRUJILLO MEDIANTE MÉTODO DIRECTO
B GENERAL ESPECÍFICOS
J
Realizar el Modelamiento y simulación ▪ Medir los abatimientos y tiempos de
espacial de la Conductividad Hidráulica filtración del petróleo y gasolina en los
por contaminación de Hidrocarburos en suelos de la zona de estudio mediante el
E suelos de Trujillo mediante Método método directo.
Directo. ▪ Determinar las propiedades mecánicas de
T
los suelos permeables adyacentes.
▪ Desarrollar un programa interactivo Veper-
I
H1, aplicando el lenguaje de programación
de Visual Basic para la simulación
hidráulica.
V ▪ Realizar el modelamiento y simulación
espacial de la conductividad hidráulica de
O
hidrocarburos en suelos permeables para un
tiempo de respuesta.
S
▪ Validar los resultados obtenidos de la
simulación mediante los datos reales
obtenidos por medición directa.
MARCO TEÓRICO
NIVEL FREÁTICO
Quintanilla (2020) define
como el nivel debajo del suelo y
rocas que estén permanentemente
saturados; su profundidad bajo la
superficie del suelo varía.
ABATIMIENTO
Según Manilla, Garnica y VEPER-H1
Martínez (2001) corresponde al
Acorde a Narváez (2021)
intervalo drenado por la
es un sistema computacional que
declinación del nivel del agua, al
emplea lenguaje de programación
multiplicarlo por el área, se
Visual Basic, este programa tiene
obtiene el volumen del suelo en el
como objeto la simulación
que se llevarán a cabo los
hidráulica de la permeabilidad.
cambios significantes en el campo
de esfuerzo efectivo.
METODOLOGÍA
FASE DE CAMPO
Se realizó mediciones directas en el terreno de
prueba, obteniendo datos de abatimiento (Z), tiempo
FASE DE GABINETE
de filtración (t) y dimensiones de la perforación. Método directo por la Interpretación de la
En esta fase de realizó gráfica de Abatimientos
perforaciones en el Método Teórico por interpretación de la
suelo de 10 cm de velocidad de Filtración
diámetro por 12 cm de
profundidad, en cada Método de la Ecuación Experimental del
perforación se colocó Cuerpo de Ingenieros de U.S. Army
un combustible para
Programa VEPER-H1
luego tomar lectura los
abatimientos respecto Mapa de Profundidad del nivel freático
al tiempo de filtración.
R
DATOS DE CAMPO Y
ABATIMIENTO
DATOS DEL PROCESO DEL MÉTODO DIRECTO
POR INTERPRETACIÓN DE GRÁFICA DE
ABATIMIENTOS E
S
U
L
T
DATOS DEL PROCESO DEL MÉTODO
DE LA ECUACIÓN EXPERIMENTAL
A
D
O
S
DATOS DEL MÉTODO TEÓRICO POR DATOS DE LA SIMULACIÓN DE LA
INTERPRETACIÓN DE LA VELOCIDAD DE PROFUNDIDAD DE TRÁNSITO DEL LÍQUIDO
FILTRACIÓN
Se hizo uso del programa VEPER-H1, determinando rangos
aceptables del comportamiento experimental de la
permeabilidad, obtenido de 27 ensayos en 9 puntos
aplicando el método directo, considerando una efectividad
del modelo como herramienta de garantía para la toma de
decisiones al reportar la simulación un % de error menor del
10% debido a la granulometría, temperatura, humedad y
variación del estrato del suelo.
Del análisis de los resultados se concluyó que, en
suelos con arena limosa y arena arcillosa
permeable a la carretera Panamericana Norte, la
conductividad hidráulica promedio de la gasolina Se logro establecer desde el punto de vista
es de 2.11x10-2 cm/s siendo mayor a la hidráulico la profundidad de filtración de la gasolina
conductividad hidráulica promedio del petróleo de y del petróleo para un mismo tiempo de derrame en
1.68x10-2cm. suelos permeables adyacentes a la carretera
Panamericana Norte: de 3 a 4 horas alcanza el nivel
freático de 2 a 3m en los cinco primeros puntos y de
Se concluyó que en los suelos permeables a la 10 horas a más en los puntos de ensayo 6 a 9 para
carretera Panamericana Norte para un mismo un nivel freático de 8 a 14 m.
tiempo de derrame la gasolina tiene mayor
profundidad de filtración que el petróleo, siendo
mas evidente cuando el tiempo de derrame se
prolonga.
ESTIMACIÓN DE PARÁMETROS EN MODELOS DE TRANSPORTE DE
O AGUA Y NITRÓGENO EN EL SUELO
GENERAL ESPECÍFICOS
B
J
Se consideraran dos modelos de 1. Analizar la sensibilidad de los parámetros en
simulación de la dinámica del agua los modelos LEACHN y EURotate_N.
y el nitrógeno: 2. Estudiar el funcionamiento de distintos
E
- LEACHM algoritmos de búsqueda, tanto local como
- EU-ROTATE global, para la determinación de los
T
parámetros de ambos modelos.
3. Aplicar distintos índices para la
determinación de la calidad de los
I
parámetros obtenidos
4. Determinar los valores de los parámetros en
los modelos LEACHN y EURotate_N en
V
distintos experimentos.
5. Utilizar los modelos LEACHN y EU-Rotate_N
una vez calibrados como sistemas de
O recomendación de abonado nitrogenado y de
evaluación de escenarios relacionados con
S
la lixiviación de nitrato
MARCO TEÓRICO
La base teórica del método consiste
básicamente en explorar el espacio de
parámetros del modelo utilizando una serie
de curvas de búsqueda y evaluar la
contribución principal de cada parámetro
sobre la salida del modelo.
METODOLOGÍA
El primer paso necesario para los modelos de transporte que
utilizan agua y nitrógeno El suelo tiene su calibración, lo que
significa que es necesario conocer los valores óptimos de los
parámetros de los que depende el modelo matemático que
describe el sistema. Alguno Estos parámetros se pueden
determinar experimentalmente o i en la literatura, pero en
otros casos deben ser determinadas por mediciones
experimentales de diversas variables en el modelo de
intervención.
Es por eso Para este propósito, es necesario desarrollar
métodos útiles de estimación de parámetros. Tipo de modelo.
Como se mencionó anteriormente, el modelo general de
transporte Agua y nitrógeno en el suelo como se introduce en
el Codex LEACHN y EU-Rotate_N requieren muchos
parámetros Decidir. Por otra parte, las puntuaciones de
distintas variables en estas previsiones Los modelos de
desarrollo comercial no son fáciles de implementar y tienen
costos Por esta razón, los métodos propuestos para estimar
estos parámetros del modelo deben poder adaptarse a un
número limitado eventos Espanol Seleccione el nivel de
paráfrasis
R
E
S
U
L
Muestran que, al igual que ocurría con el modelo
LEACHM, los dos análisis de sensibilidad
atribuyen una importancia diferente a los
T
distintos parámetros del módulo del nitrógeno de
EU-Rotate_N, independientemente de la
conguración asociada a los parámetros
A
hidráulicos ajustados que se utilice
D
O
S
C La calibración de los modelos de transporte de agua y nitrógeno
O estudiados en este trabajo, como los modelos LEACHM y
EURotate_N, requiere la estimación de un gran número de
N
parámetros a partir de los datos experimentales disponibles, en lugar
de muchos de ellos. En este caso surge el problema de la
C
reparametrización del modelo, es necesario reducir de alguna manera
el número de parámetros a calibrar, por lo que debe ser un criterio
por el cual elegir qué parámetros utilizar en el proceso de
L optimización relacionado con la calibración del modelo. .
U
S
I
O
N
E
S
O MODELAMIENTO DE LA PROPAGACIÓN DE LOS
B
FRENTESDE CONTAMINACION GENERADOS POR
EL TRANSPORTE DE COMBUSTIBLE POR
POLIDUCTOS
J
GENERAL
E Proponer un modelo matemático para predecir el crecimiento
T
de la zona de mezcla o frente de contaminación que se genera
entre combustibles derivados del petróleo como
consecuencia de su transporte secuencias por poliductos.
I
V
O
S
MARCO TEÓRICO
Poliducto
METODOLOGÍA
Ecuación de transporte de masa
ρ es la densidad de la mezcla, wA la fracción másica del
fluido desplazante, AT n es el "flux" total de masa de la
especie A, t es el tiempo y z la dirección de flujo. El término
AT n se puede descomponer como la suma de tres
componentes.
n es el "flux" de masa debido a la difusión molecular pura,
AC n es el "flux" de masa debido al movimiento total de la
mezcla "bulk" y ADC n es el "flux" de masa debido al aporte
difusivo a la convección.
La forma en que se propaga el frente de
contaminación en un tramo de poliducto de longitud L
Movimiento de la zona de mezcla a lo largo de un
tramo de tubería
R El diámetro de la línea afecta en menor escala, aunque
E
puede llegar a ser una variable determinante, pues a medida
que el diámetro crece, el área transversal disponible para la
formación de los remolinos es mayor favoreciéndose el
S
mezclado.
Los cambios más fuertes en el volumen de contaminación
U se dan por causa de la longitud que recorre el lote de
producto, no así por la variación de la tasa de bombeo.
L
T
A
D
O
S
C
O • Se logró implementar un modelo matemático simplificado que
N
permite estimar el volumen de los frentes de contaminación generados
por el transporte de productos del petróleo por baches en un
poliducto. Dicho modelo construido a partir de principios básicos de
C transferencia de masa se puede resolver en forma analítica o
numérica.
L
• La mayor dificultad para aplicar el modelo es el cálculo del
coeficiente efectivo de transferencia de masa. Una aproximación
analítica para el cálculo de dicho coeficiente sería sumamente
U complejo, motivo por el cual se han desarrollado diferentes
correlaciones empíricas para su cálculo.
S
• Antes de aplicar el modelo a un tramo de poliducto dado, éste debe
ajustarse usando datos de campo. La forma de ajuste propuesta en
este trabajo es seleccionando la correlación para el coeficiente
I
efectivo de transferencia que mejor simule una medida de campo del
volumen del frente de contaminación.
• El modelo desarrollado permite analizar el efecto de variables como
O la tasa de flujo y las dimensiones del poliducto, diámetro y longitud, en
el crecimiento del frente de contaminación. De las variables
N
mencionadas, la longitud recorrida por los baches de productos es la
que afecta en mayor medida el volumen del frente de contaminación.
E
S
O
“MODELO UNIDIMENSIONAL PARA TRANSPORTE DE
ELEMENTOS CONTAMINANTES EN DIFERENTES TIPOS DE
B SUELOS”
J GENERAL
E
Generar un modelo unidimensional para el transporte
de elementos contaminantes en diferentes tipos de
suelo.
T ESPECÍFICOS
I
• Analizar las propiedades físicas, y químicas e
hidráulicas de los diferentes tipos de suelos.
• Generar de forma experimental el transporte de
V
contaminantes en diferentes tipos de suelos
utilizando columnas de experimentación.
• Identificar las principales variables de flujo y
O
transporte para elementos contaminantes mediante
un análisis acoplado con las plataformas SEEP/W y
CTRAN/W del código Geostudio fundamentado en el
S
método de elementos finitos para suelos saturados
MARCO TEÓRICO
Para el modelamiento se utiliza el software GEOSTUDIO versión
2012, El paquete de programas Geostudio 2012 está compuesto
de varias herramientas con distintos usos y funcionalidades.
METODOLOGÍA
Se determinó usar el Sistema Munsell que es uno de los sistemas más
influyentes, fue inventado por Alberto Henry Munsell, deseó crear una forma
de describir el color.
Se determina mediante el análisis granulométrico de los suelos, es uno de
los ensayos más básicos para caracterizar un suelo, junto con los límites de
Atterberg y los ensayos de estado permiten identificar claramente un suelo
Procesamiento y análisis de datos - Software GEOSTUDIO versión 2012, extensión SEEP/W y CTRAM/W.
RESULTADOS
CONCLUSIONES
El color del estrato 1 presenta las características de un suelo con alto
contenido orgánico y de acuerdo a la Tabla de Munsell, aunque los
dos puntos difirieron levemente, se determinó que es un suelo media
fertilidad, mientras en el Estrato 2, el suelo es de media baja fertilidad
y en el Estrato 3 la fertilidad es nula, los colores son determinantes
para reconocer la capacidad del suelo y su posible constitución.
De acuerdo con el Cuadro de
clasificación USCS se
determinó que el Estrato 1
corresponde a suelo de tipo
arcilloso, el estrato 2 es un
suelo de tipo arenoso y, por
último, el Estrato 3 presenta las
características de un suelo
gravoso.
O
SIMULACIÓN DE MODELOS PROPUESTOS POR MICHAELIS-MENTEN
Y MONOD, PARA BIODEGRADACIÓN DE BENCENO, TOLUENO Y
XILENO EN SUELOS ARENOSOS, UTILIZANDO SOFTWARE
B GENERAL
J Utilizar los modelos propuestos por Michaelis-Menten y Monod, para simular los
procesos de movimiento y biodegradación de aguas contaminadas por BTX, utilizando
el software FENLAB 31 en suelos arenosos.
E ESPECÍFICOS
Determinar la viabilidad de utilizar los Determinar si esta herramienta permitirá
T modelos de Michaelis-Menten y Monod,
para simular proceso de movimiento de
aguas contaminadas con BTX, en suelos
evaluar posibles impactos ambientales
generados por contaminación de EDI, en
lechos freáticos.
I arenosos.
V
O
S
MARCO TEÓRICO
❑ Biodegradación: La biodegradabilidad es la capacidad que tienen las
sustancias y los materiales orgánicos de descomponerse en sustancias
más simples mediante la actividad (enzimática) de microorganismos.
❑ Crecimiento microbiano: incremento en el número de células
microbianas de una población. La velocidad de crecimiento es el
incremento en el número de células o en la masa celular por unidad de
tiempo.
❑ Napa freática: superficie que toma los puntos donde la presión del agua y
la presión atmosférica son iguales.
❑ Cinética: energía potencial; se encarga de analizar movimientos bajo la
influencia de fuerzas, es decir, la cinética toma en consideración también
las causas del movimiento.
❑ Régimen de flujo estacionario: el patrón de líneas de corriente
permanece inalterado en el transcurso del tiempo.
❑ Concentración microbiana: Normalmente, las bacterias crecen a
la concentración de CO2 atmosférico (0.03%). De tal forma que, la
composición del medio de cultivo puede influir sobre diferentes aspectos
en la fisiología del microorganismo y por ende en la producción de
metabolitos primarios y secundarios
❑ Reacciones catalizadas: La catálisis es el proceso por el cual cambia la
velocidad de una reacción química a causa de una sustancia llamada
catalizador.
METODOLOGÍA
MODELO CINÉTICO DE
MONOD
Monod (1950) propuso que la tasa de
crecimiento depende no sólo de la
MDELO CINÉTICO MICHAELIS- concentración de microorganismos,
MENTEN sino también de la concentración de
sustrato. Describió esta relación con
Para explicar La relación entre la velocidad inicial N y una función hiperbólica similar a la
la concentración inicial de sustrato ([S],) Michaelis y propuesta por Michaelis-Menten para la
Menten propusieron que las reacciones catalizadas interacción enzima-sustrato (Seghezzo y
enzimáticamente ocurren en dos etapas: En la Carden, 2001):
primera etapa se forma el complejo enzima-sustrato
y en la segunda, el complejo enzima-sustrato da
lugar a la formación del producto, liberando enzima
libre (Martínez, 2004)
RESULTADOS
Los resultados obtenidos en la utilización de los modelos cinéticos de Monod y
Michaelis-Menten, utilizando el FEMLAB 3.1, para su modelamiento, consideramos
satisfactorios ya que al compararse con los :saltados propuestos por Mobammed y
Allayla (1996), siguen la tendencia de las curvas obtenidas, como puede verse en las
Figuras N°06 y 07.
CONCLUSIONES
❑ Es posible utilizar los modelos de Michaelis-Menten y
Monod, para simular procesos de movimiento de aguas
contaminadas con BTX en suelos arenosos
❑ l software utilizado en este estudio fue el FENLAB [versión
1111] que resultó adecuado para la simulación de los
modelos de Michaelis-Menten y Monod, para la
biodegradación de BTX en suelos arenosos
❑ También se puede indicar que la herramienta numérica
utilizada, presenta resultados confiables y rápidos, desde
el punto de vista de evaluación del impacto ambiental, del
problema de derramamiento de BTX, por ser una causa
importante de contaminación de lechos freáticos.
REFERENCIAS BIBLIOGRÁFICAS
Narváez, R. (2021). Modelamiento y simulación espacial de la conductividad hidráulica por
contaminación de hidrocarburos en suelos de Trujillo mediante método directo [Tesis de Grado,
Universidad Nacional de Trujillo].
[Link]
[Link]?sequence=1&isAllowed=y
Manilla, A., Garnica, P. y Martínez G. (2001). Metodología para determinar la influencia de la
extracción de agua subterránea en la Infraestructura del Transporte [Archivo PDF].
[Link]
Quintanilla, S. (2020). Dinámica del nivel freático de aguas subterráneas por explotación
artesanal con fines de consumo – Comunidad de Azapampa distrito de Chilca 2019 [Tesis de Pregrado,
Universidad Continental].
[Link]
[Link]
Ibarcena W. (2014) simulación de modelos propuestos por Michaelis-Menten y Monod, para
biodegradación de benceno, Tolueno y Xileno en suelos arenosos, utilizando software. Ciencia &
Desarrollo [Link]
REFERENCIAS BIBLIOGRÁFICAS
Manuel S. O. (2016) Estimación de Parámetros en Modelos de Transporte de Agua y
Nitrógeno en el Suelo. Universidad Politécnica de Valenciafile [Archivo PDF].
[Link]
%20Estimaci%F3n%20de%20par%E1metros%20en%20modelos%20de%20transporte%20de%20ag
ua%20y%20nitr%F3geno%20en%20el%20suelo..pdf?sequence=1
Demoya, A. C., Cruz, H. Q., & Hernández, P. N. B. (2007). Modelamiento de la propagación de los
frentes de contaminación generados por el transporte de combustibles por poliducto. Dyna,
74(152), 89-96. [Link]
Mamani, I. C. & Vargas, G. V. (2019 ). Modelo Unidimensional para Transporte de
Elementos Contaminantes en Diferentes Tipos de Suelos.
[Link]