SIMULACIÓN DE INTRUSIÓN MARINA EN ACUÍFEROS CONFINADOS SIMULACIÓN DE INTRUSIÓN MARINA EN...
-
Upload
junipero-sosa -
Category
Documents
-
view
9 -
download
3
Transcript of SIMULACIÓN DE INTRUSIÓN MARINA EN ACUÍFEROS CONFINADOS SIMULACIÓN DE INTRUSIÓN MARINA EN...
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
SIMULACIÓN DE INTRUSIÓN MARINAEN ACUÍFEROS CONFINADOS
María Pool Ramírez
Grupo de Hidrología Subterránea
Universidad Politécnica de Cataluña
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
Simulación de intrusión marina en acuíferos confinados
1-Introducción.
2- Objetivos.
3- Ecuaciones Básicas.
3.1.Ecuación de Flujo.3.2.Ecuación de Transporte.
4- Metodología.
5- Resultados.
6- Conclusiones.
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
Se llama intrusión marina al movimiento permanente o temporal del agua salada tierra adentro, desplazando al agua dulce. Como estos dos fluidos (agua dulce-agua salada) son miscibles, no existe una interfaz brusca, sino que se pasa de un fluido a otro a través de una zona de mezcla, cuya anchura depende de la difusividad y dispersividad del medio y de las características del movimiento (Hidrología Subterránea. E.Custodio & R.M.Llamas)
1- INTRODUCCIÓN
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
1- INTRODUCCIÓN
Dr. Kelly Li
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
Los modelos son una herramienta muy útil, y cada día más frecuente, a la hora de gestionar acuíferos costeros con problemas de intrusión marina. Es necesario evaluar en qué casos se puede simplificar tomando la densidad constante, problema desacoplado, y modelar en dos dimensiones para evitar problemas numéricos y un alto coste.
1-Análisis de los procesos de intrusión marina en acuíferos costeros confinados estudiando el efecto de la densidad variable
2- Evaluar los posibles casos en los cuales se puede simplificar el problema a dos dimensiones.
2- OBJETIVOS
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
3- ECUACIONES BÁSICAS
0 0( ) ( )C C CC
Densidad de fluido para C=C0=0 (agua pura),
Fracción másica de soluto
ConstanteVariación de la densidad con respecto a la fracción másica.Para mezclas de agua dulce-agua salada a 20ºC, tomando C0=0 y ρ0 =998.2 (Kg/m3), el factor es aproximadamente 700 (kg/m3).
C
0 fM V
sC M M
r
w
kkv p g
S
2
2
( , (, ), ) /
( , (, ))
( , (, ), ) 1
/
r
v x y z t L s
k x y z L
k x y z t
g L s
LEY DE DARCY
flujo de agua
matriz de permeabilidades
permeabilidad relativa
vector gravedad, definido como (ELEVACIÓN).g g
Densidad Variable
Densidad Constante ρ = ρ0
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
BALANCE DE MASA DE FLUIDO
( )( )w
w p
SS v Q
t
ε porosidad
Sw grado de saturación de agua
Qp (M/(L3s) fuente de masa de fluido
γ (M/(L3 s) fuente de masa de solutos
*( )( )w
w w m w p
S Cf S vC S D I D C S C Q C
t
ECUACIÓN DE BALANCE DE SOLUTOS
Adsorción Advección Difusión-Dispersión Degradación Fuente-Sumidero
3- ECUACIONES BÁSICAS
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
4- METODOLOGÍA
5000 m
5000 m
Espesor 50 mPorosidad 0.25Permeabilidad Intrínseca 1.25e-11 m2Dispersividad longitudinal máxima 20mDispersividad longitudinal media 20mDispersividad longitudinal mínima 2mDispersividad transversal máxima, media y mínima 2m
Caudal de agua dulce preescrito Q = 1.18 hm3/año
Presión preescrita P = ρsg*(h-z)
PENDI ENTE 5%
Malla 54571 nodos48000 elementos
α
Código SUTRA.Caso teórico.
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
Se han simulado los siguientes casos:
Modelos 3D:-Problema desacoplado (densidad constante).-Problema acoplado (densidad variable).
Modelos 2D:-Problema con densidad constante.-Problema con densidad variable.
1- Análisis de sensibilidad al caudal de agua dulce.2- Análisis de sensibilidad a las dispersividades.3- Análisis de sensibilidad a las pendientes.4- Influencia de un bombeo.
4- METODOLOGÍA
-
Pozo de extracción
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
5- RESULTADOS
DIFERENCIAS
-La longitud del pie, así como el ancho de la zona de mezcla.
-La velocidades
MODELOS 3Ddensidad variable (rojo)densidad constante (azul)
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
R eference Vectors
2 .8 3 6 E -0 0 81 .9 1 9 E -0 0 6
R eference Vectors
2 .8 3 6 E -0 0 81 .9 1 9 E -0 0 6
q
q
y
DENSIDAD CONSTANTE
DENSIDAD VARIABLE
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000-300
-100
0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000-300
-100
M o d e lo 3 D co n d en s id ad C o n s tan te
M o d e lo 3 D co n d en s id ad V ariab le
5- RESULTADOS
Existen diferencias en la zona de descarga.
-Modelo con densidad constante la componente vertical es muy importante.-Modelo con densidad variable la componente principal es la transversal.
Plano XZ (Y = 0)
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
5- RESULTADOS
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
R efe ren ce V ec to rs
7 .4 0 9 E -0 0 92 .2 3 2 E -0 0 6
R efe ren ce V ec to rs
7 .4 0 9 E -0 0 92 .2 3 2 E -0 0 6
M o d e lo 2 D d en sid ad V ariab le M o d e lo 2 D d en sid ad C o n stan te
BALANCES!!!
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS Análisis de sensibilidad al caudal de agua dulce.
Caudales de 0.6 hm3/año, 1.18 hm3/año y 2.36 hm3/año. La línea marca el 90% de concentración.
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
3 D d en sid ad C o n stan te
2 D d en s id ad C o n stan te
2 D d en s id ad V ariab le
3 D d en sid ad V ariab le
5- RESULTADOS
-Al aumentar el gradiente las diferencias con respecto a la interfaz disminuyen notablemente
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
R eference Vectors
5.057638568E-0082.198217204E-006
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
Reference Vectors
05.9916084E-006
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
-En los modelos 3D al aumentar el gradiente con densidad variable las velocidades aumentan, mientras que en el caso 3D con densidad constante disminuyen.
5- RESULTADOS
Análisis de sensibilidad al caudal de agua dulce.
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
5- RESULTADOS
Análisis de sensibilidad al caudal de agua dulce.Presión en el eje Y (X = 0, Z = -300)
Pre
sión
Y (m)
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
Detalle de modelos 2D
El gradiente de concentración es elevado(desciende mucho en pequeño espacio), por eso el cambio brusco de pendiente
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
5- RESULTADOS
Análisis de sensibilidad al caudal de agua dulce.
Pre
sión
Y (m)
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
Presión en el eje Y(X = 0, Z = -300)
En los modelos 3D el cambio de pendiente coincide aproximadamente con la isolínea de 70% de concentración.
2.98E+06
3.00E+06
3.02E+06
3.04E+06
3.06E+06
3.08E+06
3.10E+06
3.12E+06
0.00E+00 5.00E+02 1.00E+03 1.50E+03 2.00E+03 2.50E+03 3.00E+03 3.50E+03 4.00E+03 4.50E+03 5.00E+03
2D constante0.5q
2D constanteq
2D constante2q
ISO constanteq
ISO constante2q
ISO constante0.5q
5- RESULTADOS
Análisis de sensibilidad al caudal de agua dulce.P
resi
ón
Y (m)
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
BALANCES
En los modelos 3D con densidad variable al aumentar el caudal de entrada al acuífero, aumenta el caudal de entrada de mar. En el resto de simulaciones (3D con densidad constante y modelos 2D) disminuye.
5- RESULTADOS
Análisis de sensibilidad al caudal de agua dulce.
0.00
10.00
20.00
30.00
40.00
50.00
60.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
Qd (kg/sg)
Qs
(kg
/sg
)285.00
290.00
295.00
300.00
305.00
310.00
315.00
320.00
3d densidad Variable 2d densidad Variable
2d densidad Constante 3d densidad Constante
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
ρ Variable ρ Constante ρ Variable ρ Constante ρ Variable ρ Constanteflujo salida (kg/sg) 18.75 18.75 37.50 37.50 74.75 75.00
flujo entra por contorno de Presión (kg/sg) 0.57 8.61 0.375 4.445 0.18 1.15
flujo sale por contorno de Presión (kg/sg) -19.31 -27.36 -37.874 -41.940 -75.18 -76.14
concentración (kg soluto/sg) 0.02 0.33 0.01 0.17 0.01 0.04
MODELOS 3D ρ Variable ρ Constante ρ Variable ρ Constante ρ Variable ρ Constante
flujo salida (kg/sg) 18.75 18.75 37.50 37.50 74.75 74.75
flujo entra por contorno de Presión (kg/sg) 21.62 315.88 35.688 306.994 49.32 290.00
flujo sale por contorno de Presión (kg/sg) -40.36 -334.58 -73.186 -344.494 -124.06 -364.75
concentración (kg soluto/sg) 0.83 12.19 1.38 11.85 1.90 11.19
Q= 2.36 hm3/añoMODELOS 2D
Q= 0.6 hm3/año Q= 1.18 hm3/año
5- RESULTADOS
Análisis de sensibilidad al caudal de agua dulce.
Tomando la densidad variable.↑Qd ↑Qs -Dispersión.-Presiones.
Con densidad constante↑Qd ↓Qs
mar
QdQs
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
5- RESULTADOS
Análisis de sensibilidad a las dispersividades.
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
120-100180-60
120-80
250-80
200-80
280-90
300-120
3D P R O M E D IO
600-200
2D
250-120 280-60
300-80 1000-150
20-10 40-10 80-40 80-80
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
Los valores de las dispersividades con las que se consigue un mejor ajuste son: Dispersividad longitudinal 280m.Dispersividad transversal 60m.
2D densidad variable (naranja)3D densidad variable promedio (azul)
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
BALANCES: modelo 3D con densidad variable Qs = 35.6 Kg/s.modelo 2D con densidad variable aumentando las dispersividades Qs = 1.29 kg/sg (valor bastante mayor que tomando las dispersividades iniciales 0.37 kg/sg).
perfil YZ
0.00000
0.00500
0.01000
0.01500
0.02000
0.02500
0.03000
0.03500
0.04000
0.04500
0 500 1000 1500 2000 2500 3000
eje Y (m) X = 0
20-2
20-10
40-10
80-40
80-80
120-80
120-100
180-60
200-80
250-80
250-120
280-60
280-90
300-80
300-120
600-200
1000-150
3d
3d promedio
5- RESULTADOS
Análisis de sensibilidad a las dispersividades.
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
5- RESULTADOSAnálisis de sensibilidada la pendiente.
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
5 0 0
1 0 0 0
1 5 0 0
2 0 0 0
2 5 0 0
3 0 0 0
3 5 0 0
4 0 0 0
4 5 0 0
5 0 0 0
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
5 0 0
1 0 0 0
1 5 0 0
2 0 0 0
2 5 0 0
3 0 0 0
3 5 0 0
4 0 0 0
4 5 0 0
5 0 0 0
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
5 0 0
1 0 0 0
1 5 0 0
2 0 0 0
2 5 0 0
3 0 0 0
3 5 0 0
4 0 0 0
4 5 0 0
5 0 0 0
0 .5 Q
0 .2 5Q
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
5 0 0
1 0 0 0
1 5 0 0
2 0 0 0
2 5 0 0
3 0 0 0
3 5 0 0
4 0 0 0
4 5 0 0
5 0 0 0
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
5 0 0
1 0 0 0
1 5 0 0
2 0 0 0
2 5 0 0
3 0 0 0
3 5 0 0
4 0 0 0
4 5 0 0
5 0 0 0
0 .1 25 Q
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
5 0 0
1 0 0 0
1 5 0 0
2 0 0 0
2 5 0 0
3 0 0 0
3 5 0 0
4 0 0 0
4 5 0 0
5 0 0 0
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
5 0 0
1 0 0 0
1 5 0 0
2 0 0 0
2 5 0 0
3 0 0 0
3 5 0 0
4 0 0 0
4 5 0 0
5 0 0 0
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
5 0 0
1 0 0 0
1 5 0 0
2 0 0 0
2 5 0 0
3 0 0 0
3 5 0 0
4 0 0 0
4 5 0 0
5 0 0 0
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
5 0 0
1 0 0 0
1 5 0 0
2 0 0 0
2 5 0 0
3 0 0 0
3 5 0 0
4 0 0 0
4 5 0 0
5 0 0 0
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
5 0 0
1 0 0 0
1 5 0 0
2 0 0 0
2 5 0 0
3 0 0 0
3 5 0 0
4 0 0 0
4 5 0 0
5 0 0 0
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
5 0 0
1 0 0 0
1 5 0 0
2 0 0 0
2 5 0 0
3 0 0 0
3 5 0 0
4 0 0 0
4 5 0 0
5 0 0 0
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
5 0 0
1 0 0 0
1 5 0 0
2 0 0 0
2 5 0 0
3 0 0 0
3 5 0 0
4 0 0 0
4 5 0 0
5 0 0 0
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
5 0 0
1 0 0 0
1 5 0 0
2 0 0 0
2 5 0 0
3 0 0 0
3 5 0 0
4 0 0 0
4 5 0 0
5 0 0 0
0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
5 0 0
1 0 0 0
1 5 0 0
2 0 0 0
2 5 0 0
3 0 0 0
3 5 0 0
4 0 0 0
4 5 0 0
5 0 0 0
Q
M O D E L O S 3 D D E N S ID A D C O N S T A N T E M O D E L O S 3 D D E N S ID A D V A R IA B L E
M O D E L O S 2 D D E N S ID A D C O N S T A N T E M O D E L O S 2 D D E N S ID A D V A R IA B L E
PENDIENTE 1%
Las diferencias existentes entre los modelos 3D y 2D con densidad constante son menores al disminuir la pendiente.Con respecto a los modelos con densidad variable, al disminuir la pendiente aumentan las diferencias.
3D densidad Cte 3D densidad Variable 2D densidad Cte 2D densidad Variable
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
Con respecto a los balances de masas ocurre lo mismo que para los modelos con pendiente del 5%.
5- RESULTADOS
Análisis de sensibilidada la pendiente.
PENDIENTE 1%
0
10
20
30
40
50
60
70
0 50 100 150 200 250
Q dulce (kg/sg)
Q s
alin
o (
kg/s
g)
0
50
100
150
200
250
300
350
2d variable 2d constante 3d variable 3d constante
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
2
0 5 10 15 20 25 30 35 40
Qf
Qs
2d variable 2d constante
Detalle de los modelos 2D
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000 0 500 1000 1500 2000 2500 3000 3500 4000 4500 50000
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000
P E N D IE N T E 1 0 % P E N D IE N T E 5 % P E N D IE N T E 1 %
5- RESULTADOS
Análisis de sensibilidada la pendiente.
Al aumentar la pendiente las diferencias aumentan al comparar los modelos con densidad constante y disminuyen al comparar los modelos con densidad variable.
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
5- RESULTADOS
Análisis de sensibilidada la pendiente.
PENDIENTE 0%Modelo 3D densidad variable-Modelo 2D sección densidad variable
0.00
20.00
40.00
60.00
80.00
100.00
120.00
140.00
160.00
0 500 1000 1500 2000 2500 3000 3500 4000
Q f (kg/s)
Qs
(kg
/s)
3d var malla ref inada 2d 0% var 2D espesor 300m 2D espesor 100m 2D espesor 200m
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
5- RESULTADOS
Influencia de un bombeo.
Pozo totalmente penetrante:
-Caudal de extracción de 4 hm3/año-Período de tiempo de 92 años.
El descenso calculado es muy similar en todos los casos:
2D-densidad constante d = 26.23m,2D- densidad variable d = 25.47m,3D-densidad constante d = 25.95m,3D densidad variable d = 25.22m.
0 5 0 0 1 0 0 0 1 5 0 0 2 0 0 0 2 5 0 0 3 0 0 0 3 5 0 0 4 0 0 0 4 5 0 0 5 0 0 00
5 0 0
1 0 0 0
1 5 0 0
2 0 0 0
2 5 0 0
3 0 0 0
3 5 0 0
4 0 0 0
4 5 0 0
5 0 0 0
0 5 0 0 1 0 0 0 1 5 0 0 2 0 0 0 2 5 0 0 3 0 0 0 3 5 0 0 4 0 0 0 4 5 0 0 5 0 0 00
5 0 0
1 0 0 0
1 5 0 0
2 0 0 0
2 5 0 0
3 0 0 0
3 5 0 0
4 0 0 0
4 5 0 0
5 0 0 0
0 5 0 0 1 0 0 0 1 5 0 0 2 0 0 0 2 5 0 0 3 0 0 0 3 5 0 0 4 0 0 0 4 5 0 0 5 0 0 00
5 0 0
1 0 0 0
1 5 0 0
2 0 0 0
2 5 0 0
3 0 0 0
3 5 0 0
4 0 0 0
4 5 0 0
5 0 0 0
0 5 0 0 1 0 0 0 1 5 0 0 2 0 0 0 2 5 0 0 3 0 0 0 3 5 0 0 4 0 0 0 4 5 0 0 5 0 0 00
5 0 0
1 0 0 0
1 5 0 0
2 0 0 0
2 5 0 0
3 0 0 0
3 5 0 0
4 0 0 0
4 5 0 0
5 0 0 0M O D E L O 3 D D E N S ID A D C O N S T A N T E
M O D E L O 2 D D E N S ID A D C O N S T A N T E M O D E L O 2 D D E N S ID A D V A R IA B L E
M O D E L O 3 D D E N S ID A D V A R IA B L E
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
concentración pozo
-2.00E-03
0.00E+00
2.00E-03
4.00E-03
6.00E-03
8.00E-03
1.00E-02
1.20E-02
1.40E-02
1.60E-02
1.80E-02
0.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 80.0 90.0 100.0
tiempo (años)
co
nc
en
tra
ció
n
3d constante 3d variable 2d constante 2d variable promedio C 3d Cons promedio C 3d var
5- RESULTADOS
Influencia de un bombeo.
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
V.I. C.I. Variable ConstanteBombeo (kg/sg) -63.42 -63.42 -63.42 -63.42
flujo salida (kg/sg) 37.50 37.50 37.50 37.50flujo entra por contorno de Presión (kg/sg) 41.51 338.69 33.77 37.13flujo sale por contorno de Presión (kg/sg) -15.32 -313.01 -7.58 -11.25
concentración (kg soluto/sg) 1.60 13.07 1.30 1.43Concentración alcanzada en el pozo 1.67E-02 1.37E-02 1.44E-02 1.36E-02
92 años de simulaciónPOZO Q 3D POZO Q 2D
pie en equilibrio (m) pozo distancia (m) tiempo (años) velocidad (m/año)
2d constante 875 3100 2225 26.7 83.332d variable 2025 3100 1075 8.29 129.67
3d constante 837.5 3100 2262.5 27.16 83.303d variable 1750 3100 1350 14.27 94.60
5- RESULTADOS
Influencia de un bombeo.
La máxima velocidad se alcanza en el modelo 2D con densidad variable, y la mínima en el modelo 3D con densidad constante.
Los balances de masas son muy parecidos en todos los casos excepto en el caso del modelo 3D con densidad constante.
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
5- RESULTADOS
Influencia de un bombeo.
CUIDADO!!!!!
NO SE HA LLEGADO AL ESTACIONARIO
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
1- Longitud del pie es mayor en los modelos con densidad variable que con densidad constante.
2- La forma de la interfaz es tendida en los modelos 3D con densidad variable, siendo muy verticalizada para el resto de los casos.
3- Balance. En los modelos 3D con densidad constante, para un Qd fijo, entra mayor cantidad de agua salada que en los modelos 3D con densidad variable, porque el gradiente de presión es mayor. En los modelos 2D la entrada de agua salina es mucho menor que en los modelos 3D.
6- CONCLUSIONES
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
4- En los modelos 3D con densidad constante la componente vertical es importante, en los modelos 3D con densidad variable la componente principal es la transversal.
5- Al aumentar el caudal de agua dulce llega un momento que el gradiente es tan elevado que la forma de la interfaz y el ancho de zona de mezcla son similares en todos los casos, sin embargo las velocidades de entrada, y los balances difieren mucho.
6- Al aumentar el caudal de agua dulce de entrada al acuífero aumentan las velocidades en los modelos 3D con densidad variable, en los modelos 3D con densidad constante sucede lo contrario. En los modelos 2D al aumentar el caudal de entrada las velocidades prácticamente no varían.
6- CONCLUSIONES
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
7- Al aumentar el caudal de agua dulce, la cuña de intrusión marina se retira pero entra mayor cantidad de agua para los modelos 3D con densidad variable. Esto en los modelos 2D y en el modelo 3D con densidad constante no queda reflejado.
8- Al aumentar los valores de las dispersividades en los modelos 2D con densidad variable se ajusta la forma de la interfaz (deja de ser verticalizada) y el ancho de la zona de mezcla, la entrada de agua salina al acuífero aumenta, pero no lo suficiente como para alcanzar los valores resultantes del modelo 3D con densidad variable.
9- Las diferencias existentes entre los modelos 3D y 2D con densidad constante aumentan al aumentar la pendiente. Con respecto a los modelos 2D y 3D con densidad variable, al aumentar la pendiente disminuyen las diferencias.
6- CONCLUSIONES
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
10- Al introducir un pozo de extracción tanto en los modelos 3D como en los modelos 2D la forma de la interfaz, el descenso, las concentraciones alcanzadas en el pozo y los balances de masas son muy similares en todos los casos, excepto en el caso del modelo 3D con densidad constante.
11- Cuando existe un gradiente importante hacia el continente, el modelo más coherente con respecto al modelo 3D con densidad variable es el modelos 2D con densidad constante. Sin embargo estos resultados han de tomarse con cautela debido a que en ninguno de los casos se ha alcanzado el estado estacionario.
6- CONCLUSIONES
SIM
UL
AC
IÓN
DE
IN
TR
US
IÓN
MA
RIN
A E
N A
CU
ÍFE
RO
S C
ON
FIN
AD
OS
SIMULACIÓN DE INTRUSIÓN MARINA EN ACUÍFEROS CONFINADOS
¡GRACIAS!