Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Geotecnia Nivel Freatico

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 142

UNIVERSIDAD NACIONAL DE SAN AGUSTÍN DE AREQUIPA

ESCUELA DE POSGRADO
UNIDAD DE POSGRADO DE LA FACULTAD DE
GEOLOGÍA, GEOFÍSICA Y MINAS

“MODELAMIENTO NUMÉRICO Y DISEÑO DEL SISTEMA DE DRENAJE PARA


EL ABATIMIENTO DEL NIVEL FREÁTICO EN LA VÍA VARIANTE DE
UCHUMAYO TRAMO II-AREQUIPA”

Tesis presentada por la Bachiller:


MARISOL MILY LLAMOCCA ALFEREZ
Para optar el Grado Académico de Maestra en Ciencias
de la Tierra, con mención en Geotecnia

Asesor:
MSC. LIZARDO GUILLERMO HUAMANI
HUAYHUA

AREQUIPA – PERÚ
2020
Dedicatoria

A Dios, mi amado Padre, quién me dio la


oportunidad de emprender este reto y en quién
encontré la fuerza y fe para culminarlo.

A mi familia, a mi papá Vicente por su apoyo


incondicional y por enseñarme con su vida a seguir
adelante a pesar de las circunstancias; a mi mamá Isabel,
una mujer admirable por ser la mejor maestra y por sus
sabios consejos, a mis hermanos Sandra y Milhuar por
estar conmigo en las buenas y en las malas y a
Emmanuel, la alegría de la casa.

A todas las personas que me ayudaron a


concluir este reto importante en mi vida profesional.

“En las alturas abriré ríos, y fuentes en medio de los


valles; abriré en el desierto estanques de aguas, y
manantiales de aguas en la tierra seca.” La Biblia.

i
Agradecimientos

Mis sinceros agradecimientos a todos los que hicieron posible que culminara este proyecto,

gracias por su apoyo, colaboración y por compartir sus conocimientos y experiencia

profesional.

Un agradecimiento especial a mi asesor de tesis el Msc. Lizardo Guillermo Huamaní

Huayhua por dedicar horas de su tiempo para absolver mis dudas, la disposición en compartir

su experiencia profesional y su constante motivación a lo largo de toda la tesis.

Al Ing. Eduardo Bautista Chambi y la empresa Rock&Soil S.A.C. por haber podido ser parte

de los estudios en la vía y por la información brindada.

Al Msc. Pablo Meza por sus aportes oportunos y su disposición para la revisión de la Tesis.

A los docentes de la Maestría de Ciencias de la Tierra de la Facultad de Geología, Geofísica

y Minas gracias por sus enseñanzas y apoyo.

ii
Resumen

El presente estudio fue realizado en el la vía Variante de Uchumayo Tramo II, en el

departamento de Arequipa. Tiene por objeto modelar el comportamiento del medio

hidrogeológico de la zona y diseñar la configuración o posición de pozos profundos para el

drenaje. El modelamiento se realizó en el software de elementos finitos FEFLOW 6.2, en el

que se simuló el comportamiento del medio en régimen permanente y en régimen variable.

La caracterización geomorfológica, geológica e hidrogeológica de la zona de estudio

permitió elaborar el modelo conceptual en el que se describe al acuífero como un medio poroso

no consolidado, que tiene el comportamiento de un acuífero libre; en el que se diferencian tres

Unidades Hidrogeológicas (Detrítico, Volcánico e Intrusivo).

El área de estudio fue de 34.1 𝑘𝑚2 este se discretizó espacialmente en una malla triangular

y bidimensional de elementos finitos. Luego, se realizó la calibración del modelo en flujo

estacionario variando los valores de permeabilidad y tomando como puntos de control las cotas

piezométricas de pozos y las obtenidas a partir de las perforaciones. En el modelo de flujo de

simulaciones transitorias se analizaron tres escenarios de drenaje, pozos laterales, pozos

centrales y en tresbolillo. A partir del análisis se ha optado como la solución óptima el arreglo

de tresbolillo, que consiste en 30 pozos de configuración triangular, construidos de 20 a 25 m

de profundidad y con capacidad de bombeo entre 5 a 10 l/s, con los que se logra cumplir los

objetivos de rebajamiento del nivel freático.

Palabras claves: Variante de Uchumayo, pozos, acuífero, modelamiento, Feflow, sistema

de drenado.

iii
Abstract

The present study was carried out on the “Variante de Uchumayo Tramo II” road, in the

department of Arequipa. Its objective is to model the behavior of the hydrogeological

environment in the area and to design the configuration or position of the deep wells. The

modeling was carried out in the finite element software FEFLOW 6.2, in which the behavior

of the medium in steady state and in variable regime was simulated.

The geomorphological, geological and hydrogeological characterization of the study area

allowed to elaborate the conceptual model in which the aquifer is described as an

unconsolidated porous medium, which has the behavior of a free aquifer; in which three

Hydrogeological Units are differentiated (Detritic, Volcanic and Intrusive).

The study area was 34.1 𝑘𝑚2 , it was spatially discretized in a two-dimensional triangular

finite element mesh. Then, the steady flow model was calibrated by varying the permeability

values and taking as control points the piezometric heights of wells and those obtained from

boreholes. In the variable flow model, three drainage scenarios, lateral wells, central wells and

staggered wells were analyzed. From the analysis, the optimal solution has been chosen the

staggered arrangement, which consists of 30 wells of triangular configuration, built from 20 to

25 m deep and with a pumping capacity between 5 to 10 l / s, with which manages to meet the

objectives of lowering the water table.

Keywords: Variante de Uchumayo, wells, aquifer, modeling, Feflow, drainage system.

iv
ÍNDICE

Dedicatoria ........................................................................................................................ i

Agradecimientos ...............................................................................................................ii

Resumen .......................................................................................................................... iii

Abstract ............................................................................................................................ iv

Capítulo I ......................................................................................................................... 1

Introducción..................................................................................................................... 1

1.1 Planteamiento del problema .................................................................................... 3

1.1.1 Situación problemática ..................................................................................... 3

1.1.2 Formulación del problema ............................................................................... 3

1.1.3 Justificación ...................................................................................................... 3

1.1.4 Objetivos .......................................................................................................... 4

1.2 Variables de la Investigación ................................................................................... 4

1.2.1 Variable Independiente .................................................................................... 4

1.2.2 Variable Dependiente ....................................................................................... 4

1.3 Hipótesis .................................................................................................................. 4

1.4 Metodología de la investigación: ............................................................................. 5

1.5 Tipo y diseño de la investigación ............................................................................ 5

1.6 Cronograma de trabajo de la investigación: ............................................................ 5

1.7 Presupuesto: ............................................................................................................. 6

Capítulo II ........................................................................................................................ 7

Marco Teórico ................................................................................................................. 7

v
2.1 Conceptos previos.................................................................................................... 7

2.1.1 Acuífero ............................................................................................................ 7

2.1.2 Tipos de acuífero .............................................................................................. 7

2.2 Formaciones geológicas como acuíferos ................................................................. 9

2.3 Parámetros hidrogeológicos .................................................................................. 10

2.3.1 Porosidad total ................................................................................................ 10

2.3.2 Porosidad eficaz y retención específica ......................................................... 11

2.3.3 Permeabilidad ................................................................................................. 13

2.3.4 Transmisividad ............................................................................................... 16

2.3.5 Coeficiente de almacenamiento ..................................................................... 17

2.4 Mapas de isopiezas: ............................................................................................... 18

2.4.1 Construcción de isopiezas .............................................................................. 19

2.5 Ecuaciones que gobiernan el flujo subterráneo ..................................................... 20

2.5.1 Ley de Darcy .................................................................................................. 20

2.5.2 Balance de masa de agua ................................................................................ 20

2.5.3 Ecuación de flujo ............................................................................................ 24

2.5.4 Régimen permanente y régimen no permanente: ........................................... 26

2.6 Hidráulica de pozos ............................................................................................... 27

2.6.1 Flujo radial a un pozo de bombeo en acuífero confinado en régimen

permanente ....................................................................................................................... 27

2.6.2 Flujo radial a un pozo de bombeo en acuífero libre en régimen permanente 28

vi
2.7 Metodología analítica para la obtención de caudales y número de pozos de bombeo

29

2.8 Modelo conceptual hidrogeológico ....................................................................... 30

2.8.1 Delimitación, descripción y geomorfología ................................................... 31

2.8.2 Geología y litoestratigrafía ............................................................................. 31

2.8.3 Propiedades hidráulicas .................................................................................. 31

2.8.4 Flujo subterráneo ............................................................................................ 32

2.8.5 Aspectos mínimos para el modelamiento conceptual .................................... 32

2.9 Modelo numérico de flujo subterráneo .................................................................. 33

2.9.1 Métodos de elementos finitos......................................................................... 34

2.9.2 Etapas del modelamiento ............................................................................... 35

2.9.3 Discretización: ................................................................................................ 35

2.9.4 Condiciones de contorno ................................................................................ 36

2.9.5 Calibración y validación de un modelo numérico .......................................... 37

2.10 Sistemas de desaguado en construcción de obras civiles................................... 37

2.10.1 Bombeos superficiales ................................................................................ 38

2.10.2 Wellpoint .................................................................................................... 39

2.10.3 Bombeo con pozos profundos .................................................................... 41

2.10.4 Sistema mixto pozos de bombeo y wellpoint ............................................. 43

Capítulo III .................................................................................................................... 44

Metodología ................................................................................................................... 44

3.1 Alcances................................................................................................................. 44

vii
3.2 Procedimiento ........................................................................................................ 44

3.3 Inventario de fuentes de agua ................................................................................ 45

3.4 Investigación de campo ......................................................................................... 47

3.4.1 Perforación diamantina .................................................................................. 47

3.4.2 Estudios geofísicos ......................................................................................... 48

3.5 Estudios de parámetros geológicos e hidrogeológicos .......................................... 56

3.5.1 Permeabilidad ................................................................................................. 57

3.6 Parámetros hidráulicos .......................................................................................... 58

3.6.1 Nivel freático .................................................................................................. 58

Capítulo IV .................................................................................................................... 59

Caracterización de la zona de estudio ......................................................................... 59

4.1 Situación ................................................................................................................ 59

4.2 Caracterización geomorfológica ............................................................................ 61

4.2.1 Geomorfología Regional: ............................................................................... 61

4.2.2 Geomorfología local:...................................................................................... 65

4.3 Caracterización geológica...................................................................................... 67

4.3.1 Geología regional ........................................................................................... 67

4.3.2 Geología local ................................................................................................ 76

4.3.3 Estratigrafía local ........................................................................................... 79

4.4 Caracterización geotécnica .................................................................................... 83

4.4.1 Para el Puente Juan Santos Atahualpa............................................................ 83

viii
4.4.2 Para el intercambio a desnivel Taboada – Rotonda ....................................... 84

4.5 Caracterización hidrometereológica ...................................................................... 85

4.5.1 Análisis hidrológico de precipitaciones pluviales .......................................... 85

4.6 Caracterización hidrogeológica ............................................................................. 87

4.6.1 UH Detrítica ................................................................................................... 87

4.6.2 UH Volcánico ................................................................................................. 89

4.6.3 UH Intrusiva ................................................................................................... 90

4.7 Modelo conceptual ................................................................................................ 93

Capítulo V ...................................................................................................................... 95

Modelamiento numérico y diseño ................................................................................ 95

5.1 Construcción de modelo numérico ........................................................................ 95

5.1.1 Geometría ....................................................................................................... 95

5.1.2 Parámetros hidrogeológicos ........................................................................... 97

5.1.3 Condiciones de contorno ................................................................................ 98

5.1.4 Condiciones iniciales y objetivos de calibración ........................................... 99

5.1.5 Calibración del modelo para la condición inicial ......................................... 100

5.2 Selección del sistema de desaguado .................................................................... 103

5.3 Diseño del sistema de desaguado ........................................................................ 104

5.3.1 Configuración de pozos laterales ................................................................. 105

5.3.2 Configuración de pozos centrales ................................................................ 106

5.3.3 Configuración de pozos en tresbolillo .......................................................... 106

ix
Capítulo VI .................................................................................................................. 113

Análisis y resultados .................................................................................................... 113

6.1 Modelo conceptual hidrogeológico ..................................................................... 113

6.2 Modelo numérico hidrogeológico ....................................................................... 114

6.2.1 Modelo hidrogeológico en estado estacionario ............................................ 114

6.2.2 Modelo hidrogeológico con pozos de desaguado ........................................ 116

Capítulo VII ................................................................................................................. 118

Conclusiones y Recomendaciones .............................................................................. 118

7.1 Conclusiones ........................................................................................................ 118

7.2 Recomendaciones ................................................................................................ 121

Referencias bibliográficas ............................................................................................. 123

x
LISTA DE TABLAS

Tabla 1. Cronograma de trabajo ............................................................................................ 5

Tabla 2. Presupuesto ............................................................................................................. 6

Tabla 3. Clasificación Meinzer ........................................................................................... 11

Tabla 4. Valores de porosidad eficaz .................................................................................. 12

Tabla 5. Métodos de laboratorio para determinar la porosidad total y eficaz ..................... 12

Tabla 6. Métodos de campo para determinar la porosidad eficaz ....................................... 12

Tabla 7. Clasificación de rocas sedimentarias no consolidadas en relación con su

permeabilidad ........................................................................................................................... 13

Tabla 8. Métodos de campo para la conductividad hidráulica ............................................ 15

Tabla 9. Diferencia de los coeficientes de almacenamiento por tipo de acuíferos ............. 17

Tabla 10. Diferencia entre métodos de diferencias finitas y elementos finitos................... 34

Tabla 11. Recomendaciones del diámetro de tubería de entubación y de perforación ....... 42

Tabla 12. Rendimiento de los métodos de perforación a diferentes formaciones............... 42

Tabla 13. Pozos dentro del área de estudio ......................................................................... 46

Tabla 14. Cotas piezométricas registradas de perforaciones............................................... 47

Tabla 15. Resumen de las perforaciones ............................................................................. 47

Tabla 16. Identificación de materiales encontrados en la línea RS-01. .............................. 49

Tabla 17. Identificación de materiales encontrados en la línea RS-02 ............................... 50

Tabla 18. Identificación de materiales encontrados en la línea RS-03. ............................. 51

Tabla 19. Identificación de materiales encontrados en la línea RS-04. ............................. 52

Tabla 20. Descripción de los estratos determinados de las líneas de MASW .................... 53

Tabla 21. Permeabilidades del sondaje DH-01 ................................................................... 57

Tabla 22. Permeabilidades del sondaje DHA-01 ................................................................ 57

Tabla 23. Parámetros geotécnicos M-01 en puente Juan Santos Atahualpa ....................... 83

xi
Tabla 24. Estados de compacidad y consistencia................................................................ 84

Tabla 25. Resumen de los N corregidos en el estribo del puente Juan Santos Atahualpa. . 84

Tabla 26. Ensayos en muestras a diferentes profundidades. ............................................... 85

Tabla 27. Ubicación de las estaciones pluviométricas ........................................................ 86

Tabla 28. Análisis de la precipitación total mensual estación la Pampilla. ........................ 86

Tabla 29. Análisis de la precipitación total mensual estación Huasacache. ....................... 86

Tabla 30. Ensayo de Bombeo. Urb. Angélica Palacios Sachaca......................................... 88

Tabla 31. Coeficiente de almacenamiento .......................................................................... 88

Tabla 32. Permeabilidad de la Formación Añashuayco ...................................................... 89

Tabla 33. Parámetros hidráulicos finales ............................................................................ 98

Tabla 34. Objetivos de Calibración ..................................................................................... 99

Tabla 35. Resultados estadísticos de la calibración del modelo ....................................... 101

Tabla 36. Ventajas y limitaciones de los sistemas de desaguado ..................................... 104

Tabla 37. Rango de valores de conductividad hidráulica asignados. ................................ 115

Tabla 38. Comparación según configuración de pozos..................................................... 116

Tabla 39.Descensos máximos de rebajamiento ................................................................ 117

xii
LISTA DE FIGURAS

Figura 1. Acuífero libre ......................................................................................................... 8

Figura 2. Acuífero confinado ................................................................................................ 9

Figura 3. Agua en los poros del suelo ................................................................................. 10

Figura 4. Conductividad equivalente .................................................................................. 14

Figura 5. Esquema de Transmisividad ................................................................................ 17

Figura 6. Esquema de coeficiente de almacenamiento ....................................................... 18

Figura 7. Representación del flujo horizontal ..................................................................... 19

Figura 8. Construcción de isopiezas.................................................................................... 19

Figura 9. Flujo en medios porosos ...................................................................................... 20

Figura 10. Balance de masa ................................................................................................ 21

Figura 11. Flujo radial en acuífero confinado ..................................................................... 28

Figura 12. Flujo radial en acuífero libre .............................................................................. 28

Figura 13. Radio equivalente .............................................................................................. 29

Figura 14. Modelo hidrogeológico conceptual ................................................................... 32

Figura 15. Esquema genérico para bombeo superficial ...................................................... 38

Figura 16. Bombeo superficial. ........................................................................................... 39

Figura 17. Esquema genérico del Sistema Wellpoint ......................................................... 40

Figura 18. Abatimiento del nivel freático con pozos profundos ......................................... 41

Figura 19. Combinación de sistemas de vacíos con pozos ................................................. 43

Figura 20. Líneas de Refracción sísmica y ensayos MASW. ............................................. 48

Figura 21, Línea de Refracción sísmica RS-01. .................................................................. 49

Figura 22, Línea de Refracción sísmica RS-02. .................................................................. 50

Figura 23, Línea de Refracción sísmica RS-03. Fuente: Rock&Soil, 2018 ........................ 51

xiii
Figura 24. Línea de Refracción sísmica RS-04. .................................................................. 52

Figura 25. Ubicación de las líneas de tomografía eléctrica................................................. 54

Figura 26. Modelo de Sección de Resistividad de la TE_01 a TE_06. ............................... 55

Figura 27. Modelo de Sección de Resistividad de la TE_07 a TE_13. ............................... 55

Figura 28. Modelo de Sección de Resistividad de la TE_01_02_03_04. ........................... 56

Figura 29. Ubicación del proyecto ...................................................................................... 59

Figura 30. Sección típica de la vía ...................................................................................... 61

Figura 31. Terraza aluvial en la zona de excavaciones para la vía principal. ..................... 65

Figura 32. Columna estratigráfica ....................................................................................... 74

Figura 33. Material puzolánico en zona de excavaciones ................................................... 76

Figura 34. Material aluvial compuesto por arenas .............................................................. 77

Figura 35. Estratigrafía Puente Juan Santos Atahualpa (Perforación DHA-01) ................. 80

Figura 36. Estratigrafía Intercambio a desnivel Taboada (Perforación DH-01) ................. 82

Figura 37. Estaciones pluviométricas de la zona de estudio ............................................... 85

Figura 38. Muestra el mallado de elementos finitos y el refinamiento ............................... 97

Figura 39. Condiciones de Contorno .................................................................................. 98

Figura 40. Objetivos de Calibración ................................................................................. 100

Figura 41. Cargas hidráulicas calculadas vs observadas. .................................................. 102

Figura 42. Rangos de aplicación de sistemas de desaguado ............................................. 103

Figura 43. Esquema de necesidades de rebajamiento. ...................................................... 105

Figura 44. Simulación con distribución de pozos laterales – Parte 1. ............................... 107

Figura 45. Simulación con distribución de pozos laterales – Parte 2. ............................... 108

Figura 46. Simulación con distribución de pozos centrales – Parte 1. .............................. 109

Figura 47. Simulación con distribución de pozos centrales – Parte 2. .............................. 110

Figura 48. Simulación con distribución de pozos en tresbolillo – Parte 1. ....................... 111

xiv
Figura 49. Simulación con distribución de pozos en tresbolillo – Parte 2. ....................... 112

Figura 50. Mapa de isopiezas en condiciones estacionarias ............................................. 116

xv
LISTA DE PLANOS

Plano N°1. Plano de Ubicación ........................................................................................... 60

Plano N°2. Plano Geomorfológico Regional ...................................................................... 64

Plano N°3. Plano Geomorfológico local ............................................................................. 66

Plano N°4. Plano Geológico regional ................................................................................. 75

Plano N°5. Plano Geológico local ....................................................................................... 78

Plano N°6. Plano Geológico del área del modelo ............................................................... 91

Plano N°7. Plano Hidrogeológico ....................................................................................... 92

Plano N°8. Plano del Modelo Conceptual........................................................................... 94

Plano N°9. Plano Límites del Modelo................................................................................. 96

xvi
Capítulo I

Introducción

La presente tesis plantea una alternativa de diseño del sistema de drenaje de aguas

subterráneas para el abatimiento del nivel freático durante las excavaciones, mediante pozos

profundos, a partir de un modelo numérico utilizando el programa FEFLOW.

En el Capítulo I se desarrolla el planteamiento del problema; en el capítulo II se dan los

fundamentos teóricos a seguir; en el Capítulo III se describe la metodología del trabajo que

inicia con la recolección de datos tales como los niveles piezométricos obtenidos del inventario

de pozos realizados por el ANA en el año 2017, además de los parámetros hidrogeológicos

obtenidos de la investigación de campo que consistió en 2 perforaciones diamantinas, 26 líneas

de tomografía eléctrica, 4 líneas de refracción sísmica y 4 ensayos de MASW; a partir de los

cuales se pudo realizar la caracterización de la zona de estudio descrito en el Capítulo IV y

asimismo elaborar los mapas morfológico, geológico, hidrogeológico y el modelo conceptual.

Se determina el comportamiento del acuífero como libre y conformado por tres Unidades

Hidrogeológicas (UH Detrítica, UH Volcánica y UH Intrusiva).

En el Capítulo V se desarrolla el modelamiento numérico y diseño. El modelamiento fue

realizado utilizando el software de elementos finitos FEFLOW; de donde se conoce que el área

de la región simulada es de 34.1 𝑘𝑚2 , discretizado en una malla de elementos finitos de 5046

nodos; los parámetros hidrogeológicos utilizados son la permeabilidad y el almacenamiento

específico. Inicialmente el modelo fue realizado en condiciones de flujo estacionario; el cual

se calibró variando las permeabilidades hasta lograr cargas hidráulicas semejantes a las

establecidas o medidas en los puntos de observación. A partir del modelo calibrado se realizó

el modelamiento en condición transitoria donde se analizaron tres arreglos diferentes de pozos

(centrales, laterales y tresbolillo).

1
En el Capítulo VI se trata el análisis y resultados. Al realizar la comparación entre las tres

configuraciones de pozos, se optó por la configuración de tres bolillos en el que se distribuyen

30 pozos de capacidad de 5 a 10 l/s y entre 20 a 25 m de profundidad con los que se alcanzan

los objetivos de rebajamiento de nivel freático. Y finalmente en el Capítulo VII se mencionan

las conclusiones y recomendaciones.

2
1.1 Planteamiento del problema

1.1.1 Situación problemática

Durante la etapa de ejecución de la obra “Carretera Variante de Uchumayo Tramo II” se

observó la presencia de filtraciones de agua entre las progresivas KM 1+560.00 al KM 2+520.

Ante este evento se planteó como solución la instalación del drenaje subterráneo que permita

rebajar el nivel freático.

En marzo del 2018 a solicitud del Gobierno Regional de Arequipa se realizó la investigación

geotécnica y estudio hidrogeológico en los puentes proyectados de Juan Santos Atahualpa y

Taboada por la empresa Rock & Soil Ingeniería SAC. la cual concluye con la presencia de agua

subterránea a una profundidad de 5.30 m en el puente Atahualpa y 5.00 m en el intercambio

vial de Taboada a partir del nivel de terreno natural.

Para la construcción de la vía principal intermedia será necesario la remoción de material,

pero debido a la presencia del nivel freático elevado es que se debe asegurar un adecuado

drenaje de las aguas subterráneas que permitan realizar las excavaciones.

1.1.2 Formulación del problema

No se cuenta con un sistema de drenaje adecuado para el abatimiento del nivel freático en

las excavaciones de la Vía Variante de Uchumayo en el Tramo II.

1.1.3 Justificación

1.1.3.1 Justificación teórica

Es necesario conocer el comportamiento del medio hidrogeológico de la zona de estudio

siguiendo los aspectos mínimos del modelamiento conceptual y numérico. Asimismo diseñar

un adecuado sistema de drenaje considerando las características del proyecto y los niveles de

rebajamiento del nivel freático que permita realizar los trabajos de excavación de manera

eficiente.

3
1.1.3.2 Justificación práctica

A través de la presente investigación se pretende desarrollar un modelo numérico de flujo

en FEFLOW que permita diseñar el sistema de drenaje para el abatimiento del nivel freático

en las excavaciones a través del bombeo con pozos profundos analizando tres configuraciones

de pozos.

1.1.4 Objetivos

1.1.4.1 Objetivo General

Diseñar el sistema de drenaje para el abatimiento del nivel freático en las excavaciones en

la zona de la Vía Variante de Uchumayo Tramo II, a partir del desarrollo de un modelo

numérico de aguas subterráneas.

1.1.4.2 Objetivos Específicos.

 Analizar la información geológica, geotécnica e hidrogeológica de la zona de estudio.

 Construir un modelo numérico de flujo subterráneo en el programa FEFLOW 6.2.

 Analizar alternativas de dispositivos de desaguado combinando el número y la

ubicación de pozos de bombeo.

1.2 Variables de la Investigación

1.2.1 Variable Independiente

La variable independiente es la posición de los pozos de bombeo.

1.2.2 Variable Dependiente

Las variables dependientes son el caudal de drenaje, el número de pozos y la eficiencia de

los conos de bombeo.

1.3 Hipótesis

A través del modelamiento hidrogeológico del área de estudio y de la simulación de pozos

de bombeos en el modelo calibrado es posible diseñar la distribución de pozos óptima que

permita alcanzar los niveles de rebajamiento del nivel freático.

4
1.4 Metodología de la investigación:

La primera etapa consiste en la recolección y análisis de la información para la construcción

del modelo conceptual. A partir del modelo conceptual se elabora el modelo numérico

hidrogeológico en el programa FEFLOW primero en estado estacionario para lograr su

calibración. Con el modelo calibrado se procede con la etapa de diseño que consiste en la

selección del sistema de desaguado y selección de la configuración óptima de pozos entre tres

alternativas: “pozos laterales”, “centrales” y “tres bolillos” a partir de simulaciones realizadas

en el modelo numérico.

1.5 Tipo y diseño de la investigación

Según el tipo de diseño la investigación es experimental ya que se manipula la variable

independiente, para ver su efecto en la variable dependiente.

El experimento consta de tres escenarios donde la variable independiente es la posición de

los pozos de bombeo que afecta las variables dependientes como es el caso del caudal de

drenaje, la cantidad de los pozos y la eficiencia de los conos de bombeo.

1.6 Cronograma de trabajo de la investigación:

Se muestra el cronograma en la Tabla 1.

Tabla 1. Cronograma de trabajo

Mes Mes Mes Mes Mes Mes Mes Mes


N° Descripción
1 2 3 4 5 6 7 8
1 Permisos y recolección de datos x
2 Análisis de información preliminar x
3 Construcción del modelo de flujo x x
Calibración y validación del
4 x x
modelo
5 Escenarios de simulación x x x
Análisis de resultados,
6 x x x
conclusiones y recomendaciones
7 Revisión final x x x

Fuente: Propia

5
1.7 Presupuesto:

El presupuesto de la presente investigación es de 2,500 soles.

Tabla 2. Presupuesto

Descripción Costo (S/)


Trámites 600
Información bibliográfica 500
Transporte 300
Impresiones 600
Otros 500
Total 2,500

Fuente: Propia

6
Capítulo II

Marco Teórico

2.1 Conceptos previos

2.1.1 Acuífero

El acuífero es una formación geológica capaz de almacenar y transmitir agua subterránea.

El agua se transmite a través de los poros o grietas (en el caso de formaciones rocosas). Los

acuíferos son fuentes de abastecimiento aprovechables para diversos usos. (Rojas, 2008)

2.1.2 Tipos de acuífero

Según Tupa (1989) se puede agrupar en tres clases de acuíferos; de acuerdo a la capacidad

de almacenamiento de agua, a la presión hidrostática del agua contenida y a las condiciones

litológicas. A continuación se desarrolla cada una de ellas.

De acuerdo a la movilidad del agua pueden clasificarse en:

 Los acuíferos, que almacenan y transmiten agua subterránea con facilidad. Por ejemplo:

las arenas, gravas, rocas compactas con fracturas apreciables y las calizas con cierto grado

de karsticidad.

 Los acuitardos, que son buenos almacenes, pero malos transmisores de agua subterránea

porque la transmiten lentamente. Por ejemplo: los limos, arenas arcillosas, areniscas y

rocas compactas con fracturas moderadas.

 Los acuicludos, que son buenos almacenes, pero nulos transmisores. El agua no puede

circular en condiciones naturales. Por ejemplo: las arcillas, los limos

 Los acuífugos que no almacenan ni transmiten agua. Por ejemplo: granitos y cuarcitas

inalteradas o no fisuradas.

De acuerdo a la presión hidrostática del agua y condiciones limitantes:

 Acuíferos libres: Aquel en que el agua subterránea del acuífero está en contacto con la

presión atmosférica a través de los poros. Al efectuar una perforación hasta su superficie

7
freática el agua permanecerá en ese mismo nivel y si la superficie freática se encuentra por

encima del terreno se generará lagunas o humedales.

Figura 1. Acuífero libre

Fuente: Adaptado de Sánchez (2017)

 Acuíferos confinados: Aquel en que el agua subterránea está a mayor presión que la

atmosférica. Al efectuar una perforación el agua ascenderá por encima del estrato, lugar

donde se equilibra la presión hidráulica con la atmosférica. Los pozos perforados en este

acuífero se llaman pozos artesianos. Si la superficie piezométrica se encuentra por encima

del terreno se generará un área surgente.

8
Figura 2. Acuífero confinado

Fuente: Adaptado de Sánchez (2017)

 Acuíferos semiconfinados: Aquel cuyo límite superior está constituido por un estrato

semipermeable que permite la filtración vertical (acuitardo) y el límite interior por un

estrato impermeable.

De acuerdo a sus características litológicas:

 Medio poroso, donde el agua fluye a través de los poros.

 Medio fracturado: el agua fluye a través de las grietas.

2.2 Formaciones geológicas como acuíferos

Las formaciones geológicas son capaces de transmitir y almacenar agua de manera distinta.

Según Custodio y Llamas (1996) los acuíferos más frecuentes son los formados por depósitos

no consolidados de materiales sueltos como gravas, arenas o mezcla de ambos estos pueden

tener un origen fluvial por ejemplo los materiales aluviales de los ríos o sus terrazas; deltaico

si se acumulan en las desembocaduras de los ríos; sedimentario por la acumulación de

9
partículas transportadas por gravedad, viento o de glaciares, son los más adecuados para la

explotación de agua subterránea. Además menciona que las propiedades hidrogeológicas de

los acuíferos que forman las rocas volcánicas dependen de las características físico químico de

las propias rocas, su grado de alteración, edad entre otros y que las rocas ígneas y metamórficas

como los granitos, dioritas, gabros, pizarras y esquistos son los de peores condiciones para el

almacenamiento de agua a menos que se encuentren en una superficie muy alterada o con

apreciable presencia de fallas y diaclasas.

2.3 Parámetros hidrogeológicos

Los parámetros hidrogeológicos nos permiten definir el funcionamiento o respuesta del

acuífero frente a determinadas acciones externas.

2.3.1 Porosidad total

El suelo está compuesto por una parte sólida y una parte hueca. La parte hueca puede estar

parcialmente ocupada por un fluido (Figura 3) o por aire. Una de las relaciones que se utilizan

para relacionar estas fases del suelo es la porosidad total.

Figura 3. Agua en los poros del suelo

Fuente: Adaptado de FCIHS (2009)

La porosidad total relaciona el volumen de vacíos que puede ser ocupado por fluidos o gases,

respecto del volumen total (Custodio y Llamas, 1996).

10
𝑉𝑜𝑙𝑢𝑚𝑒𝑛 𝑡𝑜𝑡𝑎𝑙 𝑑𝑒 𝑝𝑜𝑟𝑜𝑠 𝑐𝑜𝑛𝑒𝑐𝑡𝑎𝑑𝑜𝑠 𝑒𝑛𝑡𝑟𝑒 𝑠𝑖 (1)
𝑚=
𝑉𝑜𝑙𝑢𝑚𝑒𝑛 𝑡𝑜𝑡𝑎𝑙 𝑑𝑒 𝑙𝑎 𝑚𝑢𝑒𝑠𝑡𝑟𝑎

En la clasificación de Meinzer mostrada en la Tabla 3 vemos la relación de la estructura de

la roca con su porosidad, observamos que los depósitos sedimentarios de granulometría

homogénea presentan la porosidad más alta.

Tabla 3. Clasificación Meinzer

Tipo de estructura Porosidad


Depósito sedimentario de
Alta
granulometría homogénea
Depósito sedimentario de Media-
granulometría heterogénea Baja
Depósito sedimentario de
granulometría homogénea con Baja
minerales cementantes
Roca con huecos originados por
Baja
disolución
Roca con huecos originados por
Baja
fracturación.

Fuente: Adaptado de Custodio y Llamas (1996)

Según Custodio y Llamas (1996) en la porosidad influyen factores como la forma,

disposición y tamaño de los granos y otros factores geológicos como la compactación,

presencia de arcilla y demás.

2.3.2 Porosidad eficaz y retención específica

La porosidad eficaz relaciona el volumen de poros interconectados donde el agua es drenada

por gravedad respecto al volumen total. La porosidad total 𝑚 es la suma de la porosidad eficaz

𝑚𝑒 más la porosidad de retención o retención específica 𝑚𝑟 . (FCIHS, 2009).

𝑚 = 𝑚𝑒 + 𝑚𝑟 (2)

En la Tabla 4 se muestran algunos valores referenciales de porosidad eficaz.

11
Tabla 4. Valores de porosidad eficaz

Porosidad eficaz (%)


Material Máxima Mínima Media
Arcilla arenosa 18 3 7
Arena fina 28 10 21
Arena media 32 15 26
Arena gruesa 35 20 27
Arena con grava 35 20 25
Gravas finas 35 21 25
Gravas medias 26 13 23
Arena gruesa 26 12 22

Fuente: Adaptado de Custodio y Llamas (1996)

Existen diferentes métodos de laboratorio para determinar los valores de porosidad total y

eficaz (Tabla 5). Sin embargo los métodos de campo son los que nos dan los valores más fiables

de porosidad (Tabla 6).

Tabla 5. Métodos de laboratorio para determinar la porosidad total y eficaz

Parámetro Método
Saturación y drenado de muestras
𝑚 Métodos nucleares
Métodos volumétricos
Saturación y drenado de muestras
Correlaciones granulométricas
𝑚𝑒 Drenado por centrifugación
(humedad equivalente centrífuga)
Inyección de mercurio

Fuente: FCIHS (2009)

Tabla 6. Métodos de campo para determinar la porosidad eficaz

Parámetro Método
Bombeos controlados con
pozos de observación (Largos
periodos de tiempo)
𝑚𝑒
Balance hídrico en periodos
cortos
Ensayos con trazadores

Fuente: FCIHS (2009)

12
2.3.3 Permeabilidad

La permeabilidad es la facilidad que un medio ofrece para transmitir un fluido, en este caso

el agua, bajo un gradiente potencial hidráulico. Las propiedades que definen la permeabilidad

dependen de las propiedades del fluido y del medio poroso en conjunto (Fonseca, 2019).

2.3.3.1 Conductividad hidráulica o coeficiente de permeabilidad

La conductividad hidráulica es el parámetro que evalúa la resistencia que ofrece el suelo al

flujo del agua y es variable en los diversos medios geológicos. Los rangos van desde centenares

de m/d en gravas gruesas a valores mínimos entre 10−4 a 10−6 m/d en limos y arcillas

(Custodio y Llamas, 1996). La Tabla 7 puede servir como referencia para sus valores.

Tabla 7. Clasificación de rocas sedimentarias no consolidadas en relación con su


permeabilidad

k (m/d) 104 103 102 10 1 10−1 10−2 10−3 10−4 10−5

Gravas y Arenas finas o


Clasificación Gravas Limos o Arcillas
arenas, o arenas con mezcla de
geológica limpias arcillas arenosas compactas
gruesas limos

Clasificación Prácticamente
Buen acuífero Acuífero pobre Acuitardo
hidrogeológica impermeable

Fuente: Custodio E. - Llamas, M. R. (1996)

En un medio heterogéneo con estratos de diferentes conductividades hidráulicas, la

conductividad hidráulica equivalente es un valor global considerando al conjunto de estratos

como una unidad (Sánchez, 2017).

13
Figura 4. Conductividad equivalente

Fuente: Adaptado de Sánchez (2017)

(∑ 𝐾𝑖 . 𝑏𝑖 ) (3)
𝐾ℎ =
𝐵

Donde:

𝐾ℎ : Conductividad hidráulica horizontal equivalente

𝐾𝑖 : Conductividad hidráulica de una de las capas

𝑏𝑖 : Espesor de una de las capas

B: Espesor total, suma de todos los espesores

2.3.3.2 Métodos de laboratorio para determinar la conductividad hidráulica

En laboratorio podemos obtener la conductividad con los siguientes métodos (FCIHS,

2009):

a) Uso de permeámetros de carga fija, variable o diferencial sobre las muestras.

b) Métodos intrínsecos mediante fórmulas empíricas, granulometrías u otros, como los

siguientes:

 Hazen:

14
𝑘0 = 𝐶 ∗ (𝑑𝑒𝑓𝑖𝑐𝑎𝑧 )2 (4)

El rango de C está entre 90 y 120

 Slichter:
2
𝑘0 = 𝐶 ∗ (𝑑𝑒𝑓𝑖𝑐𝑎𝑧 ) ∗ 𝑚𝑡 𝑛 (5)

Donde n=3.3

 Kozeny, Fair y Hatch:

𝑉 2 𝑚3 (6)
𝑘0 = (𝐶 ∗ ( ) ) ∗ ( )
𝑆 (1 − 𝑚)2

𝑉 2
Donde V/S es la superficie específica de los granos y el valor de 𝐶 ∗ ( 𝑆 ) puede

1
admitirse como 200 según Custodio y Llamas (1996).

c) Gráficos de Breddin

2.3.3.3 Métodos de campo para determinar la conductividad hidráulica

Los valores obtenidos en campo serán los más confiables al igual que en los valores de

porosidad.

Tabla 8. Métodos de campo para la conductividad hidráulica

Parámetro Método
Ensayos de bombeo controlado.
Medidas puntuales. Ensayo Lefranc y
k Lugeon.
Ensayo de trazadores naturales o
artificiales.

Fuente: Adaptado de FCIHS (2009)

 Pruebas Lefranc

Estas pruebas permiten determinar la permeabilidad de los suelos arenosos, limosos o rocas

muy fracturadas por debajo del nivel freático aunque también pueden utilizarse en gravas o

15
suelos arcillosos. Existen dos tipos de ensayo Lefranc, el de carga constante y el de carga

variable. (Sánchez, 2011)

El ensayo Lefranc de carga constante consiste en medir el caudal necesario para mantener

el nivel de agua en la boca del sondaje durante 10 minutos, se va tomando nota del caudal de

agua cada minuto, en este tipo de ensayo de nivel o carga constante, se realiza en suelos cuyos

valores de permeabilidad son altos, tales como las arenas o gravas, en donde la extracción de

muestras inalteradas para realizar ensayos en laboratorio se torna imposible o muy dificultosa.

(López, 2018)

El ensayo Lefranc de carga variable es similar al anterior, pero en lugar de medir los

caudales necesarios para mantener el nivel de agua constante se mide la velocidad de descenso

del nivel de agua en el interior del sondeo, estas mediciones se realizan con una sonda (water

level meter) cada minuto o menos si el descenso es muy rápido. El método de carga variable

debe ser usado en condiciones que la permeabilidad del suelo sea lo suficientemente baja, para

permitir una determinación del nivel de agua certero, por ejemplo, limos o arcillas. (Sánchez,

2011)

2.3.4 Transmisividad

Según Custodio y Llamas (1996) el concepto de transmisividad fue dado por Theis en 1935

como el caudal que se filtra a través de una franja vertical de terreno, de ancho unitario y de

altura igual a la del manto permeable saturado bajo un gradiente hidráulico.

La transmisividad es el parámetro que resulta de la combinación de la conductividad

hidráulica y el espesor del acuífero.

𝑇 = 𝐾. 𝑏 (7)

Donde:

K: Conductividad hidráulica

b: Espesor del acuífero

16
Figura 5. Esquema de Transmisividad

Fuente: Adaptado de Sánchez (2017)

2.3.5 Coeficiente de almacenamiento

Según lo mencionado por la FCIHS el coeficiente de almacenamiento es el volumen de agua

que se obtiene de un prisma de acuífero, que tenga por base la unidad de área y por altura, la

altura saturada del acuífero, cuando se deprime o rebate el nivel piezométrico en una unidad

(2006).

En la Tabla 9 se muestra la diferencia del coeficiente de almacenamiento de acuerdo al tipo

de acuífero.

Tabla 9. Diferencia de los coeficientes de almacenamiento por tipo de acuíferos

Acuífero libre Acuífero cautivo


Efectos elásticos del agua y del acuífero.
Problemas de subsidencia o hundimiento
del terreno progresivos por explotación
Simple drenaje físico excesiva del acuífero cautivo.
Relación con los cambios de presión
atmosférica (eficiencia barométrica) o a
las mareas (eficiencia a las mareas)
S = γ · b · (m · β + α)
γ: Peso específico del agua
b: Espesor del acuífero
m: Porosidad del acuífero
S≈me
β: Coeficiente de compresibilidad del
agua
α: Coeficiente de compresibilidad del
terreno

Fuente: FCIHS (2009)

17
En el acuífero libre el coeficiente de almacenamiento es igual a la porosidad eficaz que es

el volumen de agua gravífica que contienen los poros y que se extrae por vaciado físico de una

unidad de acuífero saturado. El acuífero confinado aporta agua por descompresión y depende

del espesor del acuífero (Sánchez, 2017).

𝑉𝑜𝑙𝑢𝑚𝑒𝑛 𝑑𝑒 𝑎𝑔𝑢𝑎 𝑙𝑖𝑏𝑒𝑟𝑎𝑑𝑜 ( 8)


𝑆=
𝑉𝑜𝑙𝑢𝑚𝑒𝑛 𝑡𝑜𝑡𝑎𝑙 𝑞𝑢𝑒 ℎ𝑎 𝑏𝑎𝑗𝑎𝑑𝑜 𝑙𝑎 𝑠𝑢𝑝𝑒𝑟𝑓𝑖𝑐𝑖𝑒 𝑝𝑖𝑒𝑧𝑜𝑚é𝑡𝑟𝑖𝑐𝑎

Figura 6. Esquema de coeficiente de almacenamiento

Fuente: Adaptado de Sánchez (2017)

2.4 Mapas de isopiezas:

Al igual que un mapa topográfico representa la superficie tridimensional en un plano, el

mapa de isopiezas representa el flujo tridimensional sobre un plano horizontal mediante líneas

de flujo y líneas equipotenciales (FCIHS, 2009).

18
Figura 7. Representación del flujo horizontal

Fuente: FCIHS (2009)

2.4.1 Construcción de isopiezas

Conociendo la ubicación y cota piezométrica de varios puntos, se trazan líneas auxiliares

que una los puntos, Luego se dibujan las isopiezas utilizando el método de la triangulación

(FCIHS, 2009).

Figura 8. Construcción de isopiezas

Fuente: Adaptado de FCIHS (2009)

19
2.5 Ecuaciones que gobiernan el flujo subterráneo

2.5.1 Ley de Darcy

La ley de Darcy establece la relación entre el gradiente de potencial y el movimiento del

agua. La expresión de la ley de Darcy nos indica que la velocidad de Darcy (o flujo específico)

es proporcional al gradiente (Sánchez, 2017):

𝑑ℎ
𝑞 = −𝑘 ( )
𝑑𝑙 (9)

Donde q es el caudal que circula en 1 𝑚2 de sección, k es la conductividad hidráulica, 𝜕ℎ/ 𝜕𝑙

es el gradiente hidráulico.

Figura 9. Flujo en medios porosos

Fuente: Adaptado de Sánchez (2017)

2.5.2 Balance de masa de agua

El principio de la conservación de la masa establece que, en un volumen determinado de

medio poroso saturado, la masa del fluido que entra en un intervalo de tiempo es igual a la que

sale en el mismo intervalo de tiempo (FCIHS, 2009).

𝐸𝑛𝑡𝑟𝑎𝑑𝑎𝑠 − 𝑆𝑎𝑙𝑖𝑑𝑎𝑠 = 𝑉𝑎𝑟𝑖𝑎𝑐𝑖ó𝑛 𝑑𝑒 𝑎𝑙𝑚𝑎𝑐𝑒𝑛𝑎𝑚𝑖𝑒𝑛𝑡𝑜

A partir del volumen de control mostrado en la Figura 10 de sistema de coordenadas X, Y,

Z y de lados ∆x, ∆y, ∆z analizaremos el balance de masas en los siguientes ítems. (FCIHS,

2009)

20
Figura 10. Balance de masa

Fuente: Adaptado de FCIHS (2009)

2.5.2.1 Entradas menos salidas

Analizamos el primer término del balance que son las entradas menos las salidas.

La masa del fluido que entra por la cara 1:

𝜌𝑞𝑥 |𝑥 ∆𝑦∆𝑧 ( 10 )

Donde 𝑞𝑥 |𝑥 es el caudal o flujo específico en la dirección x en el punto de coordenadas x, 𝜌

es la densidad y ∆𝑦∆𝑧 el área transversal.

La masa del fluido que sale por la cara 2:

𝜌𝑞𝑥 |𝑥+∆𝑥 ∆𝑦∆𝑧 ( 11 )

Donde 𝑞𝑥 |𝑥+∆𝑥 es el flujo específico en la coordenada 𝑥 + ∆𝑥.

Esta última expresión por series de Taylor también lo podemos expresar de la siguiente

forma:

𝜕 ( 12 )
𝜌𝑞𝑥 |𝑥+∆𝑥 ∆𝑦∆𝑧 = 𝜌𝑞𝑥 |𝑥 ∆𝑦∆𝑧 + ∆𝑥 (𝜌𝑞𝑥 )∆𝑦∆𝑧
𝜕𝑥

21
Por lo que el balance de masas entre las caras 1 y 2 que es la diferencia de lo que entra

menos lo que sale resulta:

𝜕
𝜌𝑞𝑥 |𝑥 ∆𝑦∆𝑧 − 𝜌𝑞𝑥 |𝑥+∆𝑥 ∆𝑦∆𝑧 = − (𝜌𝑞𝑥 )∆𝑥∆𝑦∆𝑧 ( 13 )
𝜕𝑥

De la misma forma para los ejes Y y Z las entradas menos las salidas para las seis caras

sería:

𝜕𝑞𝑥 𝜕𝑞𝑦 𝜕𝑞𝑧 ( 14 )


−( + + ) 𝜌∆𝑥∆𝑦∆𝑧
𝜕𝑥 𝜕𝑦 𝜕𝑧

2.5.2.2 Variación de almacenamiento

Analizamos el otro término del balance que es la variación del almacenamiento.

La masa de agua en el volumen de control es el producto del volumen por la densidad, el

volumen de agua es el volumen total por la porosidad y el volumen del material sólido es la

diferencia del volumen total menos el volumen de agua:

𝑉𝑤 𝜌 = 𝑉𝑇 𝑚𝜌 ( 15 )

𝑀 = 𝑚𝜌∆𝑥∆y∆z ( 16 )

𝑉𝑆 = (1 − 𝑚)∆𝑥∆𝑦∆𝑧
( 17 )

Donde M es la masa de agua, m es la porosidad, 𝜌 es la densidad, 𝑉𝑤 el volumen de agua,

𝑉𝑇 es el volumen total y 𝑉𝑆 es el volumen del material sólido.

La variación temporal de la masa M respecto al tiempo sería:

𝜕𝑀 𝜕 𝜕𝑚 𝜕𝜌 𝜕∆𝑧 ( 18 )
= (𝑚𝜌∆𝑥∆𝑦∆𝑧) = 𝜌∆𝑥∆𝑦∆𝑧 + 𝑚∆𝑥∆𝑦∆𝑧 + 𝑚𝜌∆𝑥∆𝑦
𝜕𝑡 𝜕𝑡 𝜕𝑡 𝜕𝑡 𝜕𝑡

Los parámetros que pueden provocar variaciones de masa en el volumen de control son la

compresibilidad del agua que depende de la densidad de agua y compresibilidad del terreno

22
que depende de la porosidad. Ambos parámetros varían al variar la presión y sus expresiones

son respectivamente:

1 𝑑𝜌 ( 19 )
𝛽=
𝜌 𝑑𝑃

1 𝑑(∆𝑧) ( 20 )
𝛼=
∆𝑧 𝑑𝑃

Donde P es presión, 𝛽 compresibilidad del agua y 𝛼 compresibilidad del terreno.

La compresibilidad del terreno se considera despreciable es decir el volumen permanece

constante (𝑑𝑉𝑆 = 0):

𝑑𝑉𝑆 = 𝑑[(1 − 𝑚)∆𝑥∆𝑦∆𝑧] = 0 ( 21 )

Donde la diferencial del producto se reduce a dos términos:

−𝑑𝑚 ∆𝑥∆𝑦∆𝑧 + (1 − 𝑚)∆𝑥∆𝑦𝑑(∆𝑧) = 0 ( 22 )

(1 − 𝑚)𝑑(∆𝑧)
𝑑𝑚 = ( 23 )
∆𝑧

Para obtener la expresión de la variación de la densidad con la presión, sabemos que:

𝑃 = 𝑃𝑜 + 𝜌𝑔ℎ ( 24 )

Donde h es la profundidad del punto evaluado, g es la gravedad, 𝑃𝑜 es la presión en el nivel

de referencia y P la presión en el volumen de control.

Luego:

𝑑𝑃 = 𝜌𝑔 𝑑ℎ ( 25 )

Sabemos de la compresibilidad del agua que:

𝑑𝜌 = 𝜌 (𝜌𝑔 𝑑ℎ) 𝛽 ( 26 )

Lo mismo con la compresibilidad del terreno:

𝑑(∆𝑧) = 𝛼 ∆𝑧 𝑑𝑃 ( 27 )

23
𝑑(∆𝑧) = 𝛼 ∆𝑧 (𝜌𝑔 𝑑ℎ)
( 28 )

Finalmente para la porosidad tenemos:

(1 − 𝑚)𝑑(∆𝑧) ( 29 )
𝑑𝑚 = = (1 − 𝑚)𝛼𝜌 𝑑(∆ℎ)
∆𝑧

Reemplazamos las expresiones halladas en la ecuación de la variación de la masa:

𝜕𝑀 𝜕 𝜕𝑚 𝜕𝜌 𝜕∆𝑧 ( 30 )
= (𝑚𝜌∆𝑥∆𝑦∆𝑧) = 𝜌∆𝑥∆𝑦∆𝑧 + 𝑚∆𝑥∆𝑦∆𝑧 + 𝑚𝜌∆𝑥∆𝑦
𝜕𝑡 𝜕𝑡 𝜕𝑡 𝜕𝑡 𝜕𝑡

𝜕𝑀 𝜕∆𝑧 𝜕𝑚 𝜕𝜌
= (𝜌𝑚 + 𝜌∆𝑧 + 𝑚∆𝑧 ) ∆𝑥∆𝑦 ( 31 )
𝜕𝑡 𝜕𝑡 𝜕𝑡 𝜕𝑡

La ecuación de variación de almacenamiento quedaría:

𝜕𝑀 𝜕ℎ
= (𝛼𝜌𝑔 + 𝑚𝛽𝜌𝑔)𝜌∆𝑥∆𝑦∆𝑧 ( 32 )
𝜕𝑡 𝜕𝑡

La expresión del balance general se obtiene igualando la expresión resultante de las entradas

menos las salidas con la ecuación de variación de almacenamiento:

𝜕𝑞𝑥 𝜕𝑞𝑦 𝜕𝑞𝑧 𝜕ℎ


−( + + ) 𝜌∆𝑥∆𝑦∆𝑧 = (𝛼𝜌𝑔 + 𝑚𝛽𝜌𝑔)𝜌∆𝑥∆𝑦∆𝑧 ( 33 )
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑡

Simplificando:

𝜕𝑞𝑥 𝜕𝑞𝑦 𝜕𝑞𝑧 𝜕ℎ


−( + + ) = (𝛼𝜌𝑔 + 𝑚𝛽𝜌𝑔) ( 34 )
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑡

2.5.3 Ecuación de flujo

A partir de la ley de Darcy y el balance de masas procedemos a derivar la ecuación general

del flujo subterráneo (FCIHS, 2009).

Tenemos la ley de Darcy para las direcciones X, Y, Z:

24
𝜕ℎ 𝜕ℎ 𝜕ℎ ( 35 )
𝑞𝑥 = −𝑘𝑥 ; 𝑞𝑦 = −𝑘𝑦 ; 𝑞𝑧 = −𝑘𝑧
𝜕𝑥 𝜕𝑦 𝜕𝑧

Se sustituye en la Ecuación 35:

𝜕 𝜕ℎ 𝜕 𝜕ℎ 𝜕 𝜕ℎ 𝜕ℎ ( 36 )
(𝑘𝑥 ) + (𝑘𝑦 ) + (𝑘𝑧 ) = 𝜌𝑔(𝛼 + 𝛽𝑚)
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑧 𝜕𝑧 𝜕𝑡

𝑆𝑆 = 𝜌𝑔(𝛼 + 𝛽𝑚)
( 37 )

Donde 𝑆𝑆 es el coeficiente de almacenamiento.

Finalmente la ecuación general del flujo de agua subterránea es:

𝜕 𝜕ℎ 𝜕 𝜕ℎ 𝜕 𝜕ℎ 𝜕ℎ ( 38 )
(𝑘𝑥 ) + (𝑘𝑦 ) + (𝑘𝑧 ) = 𝑆𝑆
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑧 𝜕𝑧 𝜕𝑡

Tenemos algunos casos particulares de la aplicación de la ecuación general del flujo

subterráneo (FCIHS, 2009).

2.5.3.1 Estado estacionario

Aquí el flujo no cambia a lo largo del tiempo es decir los niveles no varían con respecto al

tiempo.

𝜕 𝜕ℎ 𝜕 𝜕ℎ 𝜕 𝜕ℎ ( 39 )
(𝑘𝑥 ) + (𝑘𝑦 ) + (𝑘𝑧 ) = 0
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑧 𝜕𝑧

2.5.3.2 Flujo horizontal

No depende de la coordenada Z.

𝜕 𝜕ℎ 𝜕 𝜕ℎ 𝜕ℎ ( 40 )
(𝑘𝑥 ) + (𝑘𝑦 ) = 𝑆𝑆
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑡

Si la conductividad hidráulica es la misma en ambas direcciones:

𝜕 𝜕ℎ 𝜕 𝜕ℎ 𝜕ℎ ( 41 )
(𝑘 ) + (𝑘 ) = 𝑆𝑆
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑡

25
𝜕 2ℎ 𝜕 2ℎ 𝜕ℎ ( 42 )
𝑘 2 + 𝑘 2 = 𝑆𝑆
𝜕𝑥 𝜕𝑦 𝜕𝑡

Si el espesor del acuífero b se considera constante y sabemos que:

𝑇 = 𝑘𝑏; 𝑆 = 𝑆𝑆 𝑏 ( 43 )

Donde T es la transmisividad y S coeficiente de almacenamiento.

Podemos escribir la ecuación de flujo de la siguiente manera:

𝜕 2ℎ 𝜕 2ℎ 𝜕ℎ ( 44 )
𝑇 2 +𝑇 2 =𝑆
𝜕𝑥 𝜕𝑦 𝜕𝑡

𝜕 2 ℎ 𝜕 2 ℎ 𝑆 𝜕ℎ
+ = ( 45 )
𝜕𝑥 2 𝜕𝑦 2 𝑇 𝜕𝑡

2.5.3.3 Flujo en acuífero libre

Haciendo balance de masas en un flujo en acuífero libre tenemos la siguiente ecuación:

𝜕 2 ℎ2 𝜕 2 ℎ2 𝑚 𝜕ℎ ( 46 )
+ = 2
𝜕𝑥 2 𝜕𝑦 2 𝑘 𝜕𝑡

Si además se trata de régimen permanente:

𝜕 2 ℎ2 𝜕 2 ℎ2 ( 47 )
+ =0
𝜕𝑥 2 𝜕𝑦 2

2.5.4 Régimen permanente y régimen no permanente:

El régimen del sistema es permanente, estacionario o en equilibrio cuando el flujo no varía

en el tiempo. Si el flujo varía el régimen es no permanente o variable.

En el régimen permanente, la recarga sobre un área es constante y la salida es constante e

igual a la recarga. Este sistema se da en zonas húmedas donde la precipitación es casi constante

en largos periodos. En el régimen no permanente, la recarga sobre un área no es constante al

igual que la descarga. Ocurre en zonas donde existen precipitaciones fluctuantes en periodos

cortos y zonas de altas intensidades de lluvia. (Fonseca, 2019)

26
Cuando bombeamos un pozo a caudal constante poco a poco va descendiendo el cono de

influencia. El periodo en que ocurren los descensos se llama régimen no permanente o régimen

variable (Custodio y Llamas, 1996).

2.6 Hidráulica de pozos

A partir de las siguientes condiciones analizaremos las expresiones matemáticas del flujo

en pozos de bombeo (Fonseca, 2019).

 Acuífero infinito, isótropo y homogéneo

 Pozo único

 Caudal de bombeo constante

 Radio del pozo de bombeo despreciable

 Pozo totalmente penetrante

 Espesor del acuífero constante

 No existen pérdidas de carga en el pozo

2.6.1 Flujo radial a un pozo de bombeo en acuífero confinado en régimen

permanente

El flujo de agua subterránea hacia un pozo, en un acuífero confinado y régimen permanente

fue determinado por Thiem (Custodio y Llamas, 1996).

ℎ 𝑟
𝑄 𝑑𝑟 ( 48 )
∫ 𝑑ℎ = ∫
ℎ𝑤 2𝜋𝑇 𝑟𝑤 𝑟

Integrando obtenemos:

2𝜋𝑇(𝐻 − ℎ𝑤 )
𝑄=
𝑅 ( 49 )
𝑙𝑛 𝑟
𝑤

Donde Q es el caudal constante de bombeo, 𝑟𝑤 es el radio del pozo, ℎ𝑤 es la altura hidráulica

en el pozo, H es la altura de agua inicial y R el radio de influencia.

27
Figura 11. Flujo radial en acuífero confinado

Fuente: Ferrer (2010)

2.6.2 Flujo radial a un pozo de bombeo en acuífero libre en régimen permanente

El flujo en un acuífero libre no es radial y la altura del nivel freático varía con el caudal de

bombeo dentro del radio de influencia. Aplicamos la hipótesis de Dupuit asumiendo que la

base del acuífero es horizontal y tomando la altura hidráulica con respecto a esta base (Ferrer,

2010).

Figura 12. Flujo radial en acuífero libre

Fuente: Ferrer (2010)

ℎ 𝑟
𝑄 𝑑𝑟
∫ ℎ𝑑ℎ = ∫ ( 50 )
ℎ𝑤 2𝜋𝐾 𝑟𝑤 𝑟

28
(𝐻 2 − ℎ𝑤 2 )
𝑄 = 𝜋𝐾 ( 51 )
𝑅
𝑙𝑛 𝑟
𝑤

2.7 Metodología analítica para la obtención de caudales y número de pozos de bombeo

La obtención analítica de caudales y número de pozos se puede realizar con la siguiente

metodología. Ferrer (2010).

 Cálculo del caudal del “pozo equivalente” de las dimensiones de la obra, es decir, como

si toda la obra fuera un solo pozo, posteriormente se calcula el número de pozos reales,

metodología aplicable para recintos de obras civiles menores a 200 m de longitud.

 Cálculo del radio equivalente.

Figura 13. Radio equivalente

Fuente: Adaptado de Ferrer (2010)

𝑎. 𝑏 ( 52 )
𝑟𝑒𝑞 = √
𝜋

 Determinación de la conductividad hidráulica equivalente.

 Estimación del radio de influencia.

𝑅 = 𝐶. ℎ. √𝑘 ( 53 )

29
Donde C es una constante de valor de 3000, k es el valor de la conductividad hidráulica

en m/s y h el descenso del pozo.

 Conociendo todos los parámetros necesarios aplicamos la ecuación de Thiem para

acuíferos confinados y las simplificaciones de Dupuit si se trata de un acuífero libre.

 Estimamos un diámetro de pozo para determinar el rendimiento específico máximo en

función de su diámetro y luego obtener el número de metros mojados de perforación

necesarios para bombear el caudal total.

 Con el espesor del acuífero, podremos obtener el “número de metros mojados”

disponibles para cada pozo y el potencial máximo en cada pozo. Para obtener el caudal

máximo de cada pozo multiplicamos el “número de metros mojados” por el caudal por

metro mojado exigible (caudal específico). Por tanto, conocido el caudal y el caudal

por pozo, podemos determinar el número teórico de pozos. Se recomienda elevar el

número de pozos entre 15 a 20%.

 Verificar que el diámetro del pozo no limite la instalación del equipo de bombeo

sumergible.

2.8 Modelo conceptual hidrogeológico

El modelo conceptual es la representación cualitativa y cuantitativa del sistema de agua

subterránea. Correlaciona las características geológicas e hidrológicas. Asimismo, define las

fronteras, conceptualiza el flujo, componentes de balance de agua y las propiedades hidráulicas

del acuífero subterráneo. El modelo conceptual simplifica la situación real pero necesita ser

calibrada. (Ingol y Castañeda, 2017).

Los modelos son una simplificación de la realidad partiendo de un número limitado de datos

y consideraciones que se describen a continuación.

30
2.8.1 Delimitación, descripción y geomorfología

El área de influencia física se define de acuerdo a los objetivos del modelo, es decir, se debe

estudiar la extensión del terreno donde la forma o relieve influyan directamente en el

comportamiento del modelo. Por ejemplo identificar las zonas de recarga con mayores

probabilidad de escorrentía, infiltración, recargas y descargas, zonas que son favorables a la

presencia y movilidad de agua subterránea entre otros. (Ingol y Castañeda, 2017)

2.8.2 Geología y litoestratigrafía

La caracterización geológica es la base de los estudios hidrogeológicos, en esta etapa

identificamos las rocas o estructuras geológicas y podemos analizar sus propiedades frente a la

circulación de agua o almacenamiento. (Martínez, 2018)

Es importante reconocer las unidades litoestratigráficas de acuerdo a sus condiciones

genéticas, mineralógicas y estructurales. Las unidades litológicas pueden ser Formaciones,

Unidades, Súper Unidades, entre otros. Esta descripción puede provenir de fuentes

bibliográficas, mapeos superficiales y registros de perforaciones que pueden ser

complementadas con prospección geofísica. (Ingol y Castañeda, 2017)

Los parámetros geológicos para un modelo hidrogeológico se orientan hacia el

almacenamiento y transporte del agua.

2.8.3 Propiedades hidráulicas

Se describen las diferentes características hidráulicas del medio poroso saturado, los

parámetros obtenidos pueden zonificarse para agruparlos a diferentes regiones. El

comportamiento hidráulico del medio puede estar en función de la granulometría, geología

estructural, tipo de acuífero y similares. Los valores pueden ser estimados en campo, en

laboratorio o referirse a fuentes bibliográficas (en caso donde no exista información). (Ingol y

Catañeda, 2017)

31
Figura 14. Modelo hidrogeológico conceptual

Fuente: Adaptado de FCIHS (2009)

2.8.4 Flujo subterráneo

A partir de los datos observados y las mediciones en los puntos de observación se puede

conocer la dinámica del flujo subterráneo, indicando la velocidad de flujo, zonas de

confinamiento, el análisis piezométrico y otros. A partir de esta información se pueden elaborar

mapas indicando el sentido del flujo. (Ingol y Catañeda, 2017)

2.8.5 Aspectos mínimos para el modelamiento conceptual

Según Wels (2012) mencionado en las directrices para el modelamiento de agua subterránea

en el Perú (Ingol y Castañeda, 2017) los aspectos mínimos que se deben considerar al realizar

un modelamiento conceptual son:

 Recopilación de estudios e información básica. Formaciones geológicas, dirección

del flujo, fronteras hidrológicas, parámetros hidráulicos, extracción o inyección de

pozos, cargas observadas y otros.

32
 Descripción general del sistema de agua subterránea. En base a los estudios

mencionados se describe el sistema del acuífero bajo estudio. Describir la

naturaleza de las fronteras.

 Conceptualización del flujo de agua subterránea. Con la información piezométrica

evaluar la configuración y sentido del flujo de agua subterránea. Poner énfasis en

que las cotas del terreno sean congruentes con los niveles piezométricos.

 Descripción del tipo de acuífero. Se describe y evalúa que tipo de acuífero (su

geometría, interconexión, geología o geomorfología) fueron considerados para el

modelo numérico.

 Componentes y balance hídrico preliminar. Describir y evaluar los mecanismos de

recarga y descarga.

 Propiedades hidráulicas. Describir y evaluar las propiedades de transmisión y

almacenamiento del acuífero (conductividad hidráulica, coeficiente de

almacenamiento y porosidad); para cada grupo o zona de acuerdo al marco

hidrogeológico. Para cada unidad hidrogeológica establecer un rango de valores

basado en las pruebas de bombeo desarrolladas y de la literatura existente.

 Visualización del modelo en 2D, 3D

2.9 Modelo numérico de flujo subterráneo

La naturaleza tiene sistemas complejos que dependen de numerosos parámetros que a su

vez dependen de otras magnitudes físicas, que difícilmente se puede conocer su

comportamiento mediante los métodos analíticos y es preciso el uso de métodos numéricos.

Mediante los modelos numéricos podemos reducir las ecuaciones de flujo a un sistema de

ecuaciones lineales, cuyas variables y solución son estimadas en un número finito de puntos

discretos en el espacio y en el tiempo. (Rojas, 2008)

33
La ecuación diferencial de flujo puede resolverse a través de soluciones numéricas, por los

métodos de diferencias finitas o elementos finitos. En la Tabla 10 se observamos las ventajas

y limitaciones de los métodos de diferencias finitas y elementos finitos.

Tabla 10. Diferencia entre métodos de diferencias finitas y elementos finitos

Método Ventajas Limitaciones

Diferencias finitas Conserva la masa Red rectangular


Intuitivo, simple, conocido Difícil adaptación a
y comprobado, fácil de sistemas muy complejos
programar. con capas deformes.
Pozos extraen agua de toda
la celda.

Elementos finitos Discretización más No tan intuitivos, teoría


eficiente en zonas de más compleja y abstracta
interés solamente (malla Conserva masa solo
triangular) globalmente
Pozos pueden ser Solución más compleja,
representados como nudos discretización más
discretos elaborada, requiere
Bien adaptado a geología software especializados.
compleja.

Fuente: Rojas (2008)

2.9.1 Métodos de elementos finitos

El método de elementos finitos fue usado inicialmente en el estudio de los campos

potenciales en aplicados a problemas de flujo de fluidos luego fue aplicado al flujo de medios

porosos. (Rojas, 2008)

En el método de los elementos finitos, se discretiza el espacio en elementos triangulares y

se obtiene una ecuación algebraica para cada elemento, en función de la altura piezométrica de

34
cada una de ellas y de los contiguos, y de los parámetros físicos del acuífero. La resolución de

ecuaciones lineales algebraicas, en función de la altura piezométrica en cada elemento, da una

solución aproximada de la ecuación diferencial en una serie de puntos discretos (Custodio y

Llamas, 1996).

2.9.2 Etapas del modelamiento

Después de elaborar el modelo conceptual hidrogeológico, de forma general Fonseca (2019)

menciona la siguiente secuencia de pasos:

 Diseño del modelo: Se selecciona las condiciones de borde, condiciones iniciales, rango

de valores y distribución de parámetros, definir puntos y metas de calibración y definir

los escenarios de simulación.

 Construcción del modelo: Preparar la malla, ingresar las condiciones de borde, ingresar

las propiedades, observaciones, fuentes y sumideros, definir los periodos de estrés y

pasos de tiempo (discretización temporal).

 Calibrar y validar el modelo: Ajustar los parámetros del modelo inicial hasta llegar a

las metas de calibración; usar el modelo en otra situación de estrés para comprobar su

calidad en reproducir datos.

 Evaluar la sensibilidad del modelo a cambios de parámetros dentro de un rango

aceptable.

 Simular escenarios futuros.

 Evaluar la incertidumbre de las predicciones

2.9.3 Discretización:

La discretización consiste en dividir el sistema en elementos, se ubican una serie de puntos

con coordenadas conocidas y se unen mediante segmentos para formar triángulos en caso de

un modelamiento en 2D. (Fonseca, 2019)

35
Según la FCIHS (2009), para realizar la discretización se debe tener en cuenta las siguientes

consideraciones:

 Ajustarse a contornos y contacto entre materiales.

 Evitar ángulos muy agudos u obtusos.

 Procurar que los nudos coincidan con los pozos.

 Que la malla se adapte a red de flujo anisotropía.

 El tamaño de elementos debe ser pequeño en zonas de alto gradiente y tener cambio de

tamaño suave en lo posible.

2.9.4 Condiciones de contorno

Las condiciones de contorno son los valores que toman los campos en las fronteras de los

dominios del modelo. Desde un marco físico las condiciones de contorno en el modelamiento

de aguas subterráneas son de tres tipos (FCIHS, 2009):

 Nivel fijo (Condición de tipo Dirichlet o de primer orden):

En esta condición la carga hidráulica no resulta de la simulación, en estos nodos la carga

es el valor de condición de frontera. Las condiciones se aplican en casos donde se conoce

el potencial hidráulico, cuando se presente contacto entre el acuífero y el mar o un río y

donde se produce recarga o descarga subterránea del acuífero.

Se debe asignar este límite fijando el nivel a los nodos donde no habrá necesidad de

resolver la ecuación de flujo subterráneo.

 Caudal conocido (Condición de Neumann o de segundo orden):

Esta condición aplica un flujo predefinido a nodos a lo largo de una línea (modelo 2D) Las

condiciones de límite de flujo se aplican en los casos en que la velocidad del gradiente,

entrada o salida se conoce de antemano. (Fonseca, 2019)

Se utiliza cuando el acuífero se encuentra en contacto con formaciones geológicas

donde no se produzca entradas o salidas de agua (zonas impermeables)

36
 Mixto (Condición de Cauchy o de tercer orden)

Esta condición de borde se utiliza cuando debido al tamaño del acuífero no se puede

modelar por completo y se asume gradiente constante, en contornos permeables donde

el gasto de entrada de salida es conocido. Se debe asignar el gradiente deseado a los

nudos del borde.

Puede ser utilizado en los casos de ríos o lagos con una conexión limitada al agua

subterránea (el nivel del agua superficial es igual al nivel del agua subterránea) o para

desagües parcialmente obstruidos. (Fonseca, 2019)

2.9.5 Calibración y validación de un modelo numérico

La calibración es la etapa del modelamiento en la cual se ajusta, verifica calibra o valida el

modelo para lograr una buena aproximación entre los valores del registro conocido y los valores

calculados. Se puede variar los valores de partida (coeficiente de permeabilidad, coeficiente de

almacenamiento, transmisividad) hasta conseguir que la diferencia entre los valores de los

niveles conocidos y los niveles calculados estén dentro de la precisión necesaria. (Rojas, 2008)

En el proceso de validación se utiliza el modelo en otra situación de estrés para comprobar

su calidad en reproducir los datos. (FCIHS, 2009)

2.10 Sistemas de desaguado en construcción de obras civiles

Los sistemas de desaguado subterráneo o control de nivel freático en obras civiles ha

evolucionado en las últimas décadas lográndose obtener sistemas más eficientes. (Ferrer,

2010).

Estos sistemas pueden ser barreras físicas que impidan el acceso de agua a la excavación,

bombas que depriman el nivel freático o la combinación de ambos. Los pozos de bombeos

solos o combinados son una solución económica a este problema, su utilización dependerá de

la naturaleza del terreno y la profundidad de la excavación. (Cassiraga y Ferrer, 2008)

Los principales métodos basados en pozos de bombeos son (Ferrer, 2010):

37
 Bombeos superficiales

 Wellpoint

 Bombeo con pozos profundos

 Sistemas mixtos pozos de bombeo y wellpoint

2.10.1 Bombeos superficiales

Los bombeos superficiales drenan el agua subterránea procedente del fondo de una

excavación ubicado bajo el nivel freático, pueden conducirse a través de zanjas perimetrales

hasta su captación. Generalmente el área de excavación no rebaja el nivel freático hasta tener

el área seca, sino que conduce el agua de escorrentía hasta los puntos de captación. Para evitar

los arrastres se debe minimizar la velocidad de flujo. (Ferrer, 2010)

Figura 15. Esquema genérico para bombeo superficial

Fuente: Adaptado Ferrer (2010)

a) Aspectos constructivos:

De acuerdo al reporte del Department of the Environment and the North West Water

Authority (2015) se debe considerar lo siguiente:

 El bombeo superficial debe ubicarse en una cota de profundidad suficiente para recibir

el agua captada por medio de un flujo por gravedad en las zanjas, dentro de esta se

38
coloca la tubería perforada rodeada de material drenante, si es necesario se utilizará una

membrana de geotextil para evitar el paso de materiales finos.

 En excavaciones de pequeñas dimensiones, las bombas se distribuyen perimetralmente

y al exterior de la excavación.

 Dentro del sumidero se instala una tubería perforada, el tubo es colocado dentro de una

calicata y es rodeado por material prefiltro. El tamaño de la ranura y la granulometría

dependerán del tipo de suelo.

Figura 16. Bombeo superficial.

Fuente: Ferrer (2010)

2.10.2 Wellpoint

El sistema wellpoint combina el bombeo de agua y bombeo de aire que causan la depresión

del nivel freático por medio de filtros que se introducen en el suelo. Los componentes

principales de este sistema son los filtros wellpoint, los cuales se encuentran ubicados

perimetralmente a la excavación, la tubería ascendente, los latiguillos flexibles, el colector

general y la bomba de aspiración que descarga el agua fuera del área de trabajo. (Ferrer, 2010)

39
Figura 17. Esquema genérico del Sistema Wellpoint

Adaptado: Ferrer (2010)

a) Aspectos constructivos:

Según Ferrer (2010) los aspectos a considerar son los siguientes:

 En terreno granular sin finos, se hincan directamente en el terreno un número

considerable de lanzas, tubos de 50 mm de diámetro, por medio de inyección de agua a

presión entre 20 a 25 lt/s, estos en su tramo más profundo adopta una forma de filtro.

 En caso de terrenos con alto contenido de partículas finas se instalan los tubos dentro

de pozos drenantes, el espacio anular resultante se rellena con arenas gruesas y grava a

excepción de los 1.50 m que se rellenan de limos arcillosos.

 Los filtros wellpoint deben ubicarse entre 2 a 6 m debajo del nivel freático y el colector

principal en la cota del nivel freático.

40
 La altura de agotamiento puede alcanzar los 7m, en caso de mayor altura, puede

colocarse un sistema escalonado de varias líneas de Wellpoint (Catálogo Wellpoint de

Ischebeck, 2019).

2.10.3 Bombeo con pozos profundos

Este sistema consta de un conjunto o red de pozos profundos con bombas sumergibles, estos

generan el descenso del nivel freático en una gran extensión.

Para el diseño de este sistema se debe conocer el comportamiento del acuífero y la influencia

de la obra a proyectar en su entorno, esto se logra a través de modelos analíticos o matemáticos

de acuerdo a la complejidad. Si podemos conocer los caudales de bombeo esperados con una

adecuada exactitud podremos elegir el número y tipo de equipos de bombeo necesarios.

(D.W.W.A., 2015)

Figura 18. Abatimiento del nivel freático con pozos profundos

Fuente: Propia

a) Aspectos constructivos:

Según Ferrer (2010) los aspectos a considerar son los siguientes:

 Se realiza el entubado con una tubería de PVC con un tramo con ranuras en la zona de

filtro que permita el paso del agua. En el caso de terreno con alta presencia de finos

(limos y arcillas), la perforación será mayor al diámetro del tubo y en el espacio vacío

anular se coloca material granular drenante.

41
 El equipamiento del material a utilizar en este sistema por ser de uso temporal suele ser

más económicos en comparación con los utilizados para suministro de agua.

 Se debe garantizar la estabilidad del fondo de excavación por lo que se dispone de un

margen de seguridad de 2 m por debajo de la excavación, en el presente estudio se

considera mantener el nivel de agua a 3 m del fondo de la excavación.

 El diámetro de la tubería se elige de acuerdo a la dimensión de la bomba que se va a

instalar. Asimismo, el método de perforación a realizar dependerá de las condiciones

del terreno y las limitaciones que tiene el área del proyecto.

Tabla 11. Recomendaciones del diámetro de tubería de entubación y de perforación

Caudal Diámetro Diámetro


máximo de mínimo mínimo
bombas recomendado recomendado
sumergibles de tubería filtro de perforación
(l/s) (mm) (mm)
5 125-152 250-275
10 152-203 300-325
15 165-250 300-375
20 180-250 300-375
25 203-300 325-425
44 250-350 375-475

Fuente: Preene y Cashman (2001)

Tabla 12. Rendimiento de los métodos de perforación a diferentes formaciones

Percusión Rotación Dual o OD


Formación
(a cable) (con lodos) (con martillo)
Arenas y gravas Medio Rápido Muy rápido
Arcillas y limos Medio Rápido Medio
Arenisca
cementada Lento Medio Medio
No
Basaltos granitos Lento Lento recomendado
Basaltos No
fracturados Lento Imposible recomendado
Areniscas con Lento a No
cavernas Lento imposible recomendado

Fuente: Ferrer (2010)

42
2.10.4 Sistema mixto pozos de bombeo y wellpoint

Cada sistema de bombeo presenta limitaciones en su uso, por ejemplo el sistema de

wellpoint en arenas y gravas tiene una altura de succión de 5 m a 6 m pero puede verse limitada

de 3.5 m a 4.5 m en suelos más finos o casos similares (Ferrer, 2010). Por lo tanto, en casos

donde se necesite drenar grandes profundidades y con áreas apropiadas para el uso de

Wellpoints se puede optar por la combinación de sistemas.

Figura 19. Combinación de sistemas de vacíos con pozos

Fuente: Adaptado de Ferrer (2010)

43
Capítulo III

Metodología

3.1 Alcances

La elaboración del presente trabajo se inicia con la recolección de información. A partir del

inventario de los pozos de Arequipa elaborado por la Autoridad Nacional del Agua (2017) se

obtuvieron datos de pozos y niveles piezométricos para control. La profundidad del nivel

freático fue analizada a partir del “Informe del estudio hidrogeológico del proyecto:

Mejoramiento de la carretera Variante de Uchumayo, Tramo II” y el expediente de “Drenaje

subterráneo del proyecto: Mejoramiento de la carretera Variante de Uchumayo, Tramo II”

proporcionados por la empresa Rock&Soil S.A.C. (2018). Las variables para la caracterización

geológica, hidrogeológica y geotécnica del “Informe de investigaciones geotécnicas mediante

perforación diamantina - Puente Atahualpa y Taboada” realizado por la empresa PERFOCON

S.A.C. (2018). Estas fueron las fuentes principales de información que se utilizaron.

Para la etapa de modelamiento y simulación se utilizó el software FEFLOW 6.2. El modelo

numérico simula el flujo en dos dimensiones mediante aproximaciones de elementos finitos.

La calibración se realizó con las cargas hidráulicas de los puntos de control. La etapa del diseño

consistió en optimizar la configuración de los pozos profundos, para ello se simuló tres

escenarios; pozos laterales, pozos centrales y tres bolillos, eligiendo la opción de “tres bolillo”

con la que se logró los objetivos de desaguado con un número adecuado de pozos y con la

ventaja que estos se encuentren ubicados fuera del área de excavaciones. En el Capítulo IV se

explicará a detalle el desarrollo de la etapa de modelamiento y diseño.

3.2 Procedimiento

El procedimiento realizado para el estudio se describe a continuación:

1. Análisis de la información obtenida.

2. Delimitación del área de influencia (Plano N°9) en escala 1/25000.

44
3. Sobre la base de la investigación de campo y los estudios geológicos del INGEMMET

se elaboran los planos geomorfológicos (Plano N°2 y N°3), geológicos (Plano N°4, N°5

y N°6) e hidrogeológico (Plano N°7).

4. Elaboración del plano conceptual (Plano N° 8) de acuerdo al análisis de las cargas

hidráulicas de los pozos y los niveles de descenso del nivel observados en el área del

proyecto (Puente Taboada e Intercambio a desnivel de Atahualpa) además de los

parámetros hidrogeológicos estudiados.

5. Elaboración de la malla triangular, se asigna la permeabilidad y la transmisividad.

6. Calibración del modelo en base a la piezometría calculada.

7. Simulación del modelo y diseño (configuración) de los pozos.

3.3 Inventario de fuentes de agua

El INRENA realizó el estudio hidrogeológico para la evaluación y ordenamiento en la

cuenca del río Chili en el año 2003 donde muestra el inventario de 645 pozos de agua

subterránea, de los cuales 378 se encontraban operativos. En el año 2014 la Autoridad Nacional

del Agua (ANA) monitoreó los niveles estáticos de 70 pozos. Y nuevamente en el año 2015 y

2017 se realizaron campañas de los mismos pozos. A partir del inventario realizado el año 2017

y de acuerdo al área de influencia del estudio es que se seleccionaron 11 pozos según se muestra

en la Tabla 13.

En la zona de la vía principal se obtuvieron 3 puntos de control a partir de perforaciones y

observaciones en campo (Tabla 14).

45
Tabla 13. Pozos dentro del área de estudio

Coordenadas Niveles del agua y caudal

Cota del p.r. suelo N. estático caudal


N° Nombre /Empresa Distrito
terreno
Este Norte
PROF.
(m) m.s.n.m. l/s
(m)
FIORELA CARMEN CELMIRA
P-22 Yanahuara 228168 8185975 2206 0 8 2312.00 1.02
UGARTE POSTIGO
P-41 INCALPACA TPX S.A. Sachaca 227137 8184683 2315 0 30.64 2284.36 5.7
UNIÓN CERVECERA BACKUS Y
P-42 Sachaca 227288 8184425 2285.38 0.5 25.28 2265.10 70
JOHNSTON S.A.A. BACKUS
UNIÓN CERVECERA BACKUS Y
P-44 Sachaca 226 988 8184426 2287.49 0.45 21.38 2275.00 30
JOHNSTON S.A.A. BACKUS
EMBOTELLADORA SAN MIGUEL
P-45 Sachaca 224640 8183980 2279 0.00 35.00 2275.00 3.00
DEL SUR SAC
MUNICIPALIDAD DISTRITAL DE
P-50 Sachaca 223171 8183732 2274 0 14.8 2253.00 16.42
SACHACA
Cerro
P-405 CAMAL FRIGORÍFICO DON GOYO 224798 8187247 2376.17 0.56 27.06 2355.00 10
Colorado
Cerro
P-416 EMPRESA CEPCO AVÍCOLA E.I.R.L. 225271 8185200 2320 0 25 2310.0.00 4
Colorado
KRISTEL STEPHANIE ABARCA Cerro
P-420 224290 8185884 2341 0 23.7 2325.00 3.1
LINARES Colorado
GRANJA LA RINCONADA DEL SUR Cerro
P-421 225024 8184760 2311 0 10 2299.00 1.8
S.A. Colorado
Cerro 81850
P-427 PITATA SAC. 224152 2310 0 11.24 2305.76 2
Colorado 15

Fuente: Adaptado de Inventario de pozos del ANA (2017)

46
Tabla 14. Cotas piezométricas registradas de perforaciones

Coordenadas
Cota del Cota
IRHS
Este Norte terreno Piezométrica

CP-1 225478.44 8184597.44 2303 2298


CP-2 225717.64 8184611.9 2300 2295
CP-3 226048.24 8184476.29 2297 2282

Fuente: Propia
3.4 Investigación de campo

3.4.1 Perforación diamantina

Se realizó un total de 02 perforaciones diamantinas. Estas perforaciones fueron realizadas

en febrero del 2018. La Tabla 15 muestra un resumen de las perforaciones realizadas.

Tabla 15. Resumen de las perforaciones

Coordenadas UTM (WGS 84) Cota


N° Códigos Ubicación
Este(m) Norte(m) msnm
Puente J. S.
1 DHA-01 225478.44 8184597.44 2303 Atahualpa (Prog.
1+750)
Estribo derecho
2 DH-01 225717.64 8184611.90 2300 Puente Taboada
(Prog. 2+290)

Fuente: Propia

El método de perforación diamantina (DDH) utilizado fue con recuperación de muestras en

diámetro HQ (96 mm) y con tuberías de revestimiento PQ (114.3 mm). Este método de

perforación implicó el uso de brocas de corona diamantada que, mediante procesos de rotación,

corta al material y por medio de un sistema de izaje la muestra es extraída hacia la superficie

dentro del tubo interior. Como fluido de perforación únicamente se usó agua (sin aditivos), y

el agua de retorno de la perforación fue almacenada en pozas debidamente impermeabilizadas.

(Rock&Soil, 2018). En el Capítulo III se describirá la estratigrafía de cada perforación.

47
3.4.2 Estudios geofísicos

Se realizaron 26 líneas de Tomografía Eléctrica, 4 líneas de Refracción Sísmica y 4 ensayos

MASW que sirvieron para definir los diferentes estratos y la ubicación de materiales saturados,

a continuación se realiza la descripción de estos de acuerdo al informe geofísico elaborado por

la empresa Rock&Soil (2018) .

3.4.2.1 Refracción sísmica y MASW

El estudio sísmico se realizó en tres sectores diferentes del tramo II de la variante

Uchumayo, el primero en la zona del puente Circunvalación, el segundo en el puente Juan

Santos Atahualpa y el tercero en el intercambio a desnivel Taboada-Rotonda.

Figura 20. Líneas de Refracción sísmica y ensayos MASW.

Fuente: Rock&Soil (2018)

48
 Línea RS-01

Se ubica en la zona del puente de Circunvalación. Los estratos de la sección sísmica se

describen en la Tabla 16.

Tabla 16. Identificación de materiales encontrados en la línea RS-01.

Vp
Vp Espesor
Línea Estrato Promedio Descripción
(m/s) (m/s) (m)

Material suelto de fácil


ripabilidad y de baja
A 400 - 555 4.5 – 14.2 compacidad asociado a suelo
810 de cobertura correspondiente
a arenas limosas.
Material no ripable
RS-01 1330 - posiblemente asociado a
B 1855 12.5 – 22.3
2380 gravas con alto grado de
saturación.
Material no ripable de alta
3100 - compacidad posiblemente
C 3800 3450 -
asociado a roca ligeramente
fracturada.

Fuente: Rock&Soil (2018)

Figura 21, Línea de Refracción sísmica RS-01.

Fuente: Rock&Soil, 2018

49
 Línea RS-02

Se ubica en la zona del puente Juan Santos Atahualpa, presentando la siguiente sección.

Tabla 17. Identificación de materiales encontrados en la línea RS-02

Vp Vp
Espesor
Línea Estrato Promedio Descripción
(m/s) (m/s) (m)

Material suelto de fácil


ripabilidad y de baja
A 400 - 900 650 7.5 – 14.0 compacidad asociado a
suelo intemperizado con
presencia de arenas limosas.
Material no ripable
RS-02 2300 -2800 posiblemente asociado a
B 2550 5.0 – 10.8
gravas con contenido de
humedad.
Material no ripable de alta
3300 -3800 compacidad posiblemente
C 3550 -
asociado a roca ligeramente
fracturada.

Fuente: Rock&Soil (2018)

Figura 22, Línea de Refracción sísmica RS-02.

Fuente: Rock&Soil (2018)

50
 Línea RS-03

Se ubica en la zona del intercambio a desnivel Taboada-Rotonda, los estratos se

describen a continuación.

Tabla 18. Identificación de materiales encontrados en la línea RS-03.

Vp Vp
Espesor
Línea Estrato Promedio Descripción
(m/s) (m/s) (m)

Material de ripado fácil no muy


bien consolidado e
A 400 - 1300 intemperizado, asociado a
850 3.6 – 5.2
depósitos superficiales
correspondientes a suelos
blandos como arenas sueltas.
RS-03 Material no ripable
B 2300- 2800 10.0 –
2550 posiblemente asociado a suelo
2.8
denso o aluvión muy compacto.
Material no ripable de alta
C 3000- 3800 compacidad posiblemente
3400 -
asociado a roca ligeramente
fracturada.

Fuente: Rock&Soil (2018)

Figura 23, Línea de Refracción sísmica RS-03. Fuente: Rock&Soil, 2018

Fuente: Rock&Soil (2018)

51
 Línea RS-04

La línea se desarrolló en la zona perpendicular al intercambio a desnivel Taboada-

Rotonda, los estratos se describen a continuación.

Tabla 19. Identificación de materiales encontrados en la línea RS-04.

Vp Vp
Espesor
Línea Estrato Promedio Descripción
(m/s) (m/s) (m)

Material de ripado fácil, asociado


400- a depósitos superficiales como
A 1100 750 3.0 – 7.5 suelos de cobertura, suelos
blandos, arenas sueltas con
presencia de material gravoso.
2100 - Material no ripable posiblemente
RS-04 2400
B 2700 15.2– 21.1 asociado a suelo muy denso con
contenido de grava húmeda.
Material no ripable de alta
3300 - compacidad posiblemente
C 3800 3550 -
asociado a roca ligeramente
fracturada.

Fuente: Rock&Soil (2018)

Figura 24. Línea de Refracción sísmica RS-04.

Fuente: Rock&Soil (2018)

52
 Ensayos MASW

En los ensayos de MASW se diferenciaron dos tipos de suelos. En el puente de

Circunvalación se observa un suelo de tipo C asociado a suelo muy denso o roca blanda y

para la zona del puente Juan Santos Atahualpa y el intercambio a desnivel Taboada–

Rotonda se asocia a un suelo rígido. (Rock&Soil, 2018)

Tabla 20. Descripción de los estratos determinados de las líneas de MASW

Vs Vs
Vs30 Espesor
Línea Estrato Promedio Descripción
(m/s)
(m/s) (m/s) (m)

100 - Material
400 posiblemente
A 200 5.2
asociado a suelo muy
denso o roca blanda.
MASW 01 458.3 Material asociado
480 - posiblemente a suelo
710 muy denso o roca
B 595 -
blanda con presencia
de una delgada capa
de suelo rígido.
Material
150 - posiblemente
A 190 4.7
230 asociado a suelo
rígido.
MASW 02 215.8
Material asociado
230 - posiblemente a suelo
B 390 -
550 muy denso o roca
blanda.
Material asociado a
170 – depósitos
A 185 5.1
200 superficiales o suelo
rígido
MASW 03 222.8 Material
posiblemente
asociado a suelo
B 200 - 325 -
rígido con presencia
450
de suelo de roca
blanda.

Fuente: Adaptado de Rock&Soil (2018)

53
3.4.2.2 Tomografía eléctrica

Las líneas de tomografía eléctrica se ubican a lo largo de la vía, dentro de las excavaciones

y en secciones transversales, como se muestra a continuación.

Figura 25. Ubicación de las líneas de tomografía eléctrica

Fuente: Rock&Soil (2018)

A continuación se realiza la descripción de las líneas más representativas de la zona del

proyecto.

3.4.2.3 LINEA TE_01 a TE_06 (Progresivas 1+780 a 2+500)

Los valores de la resistividad son variables con los cuales se sectorizo en tres tipos de

litologías (Rock&Soil, 2018):

Zona A: Los valores de resistividad oscilan entre 20 – 80 ohm (tonalidad azul a turquesa),

materiales asociados a zonas con saturación, con alto contenido de agua.

Zona B: Los valores de resistividad oscilan entre 80 – 200 ohm (tonalidad verde), materiales

asociados a arenas – limosas con mediano contenido de humedad.

54
Zona C: Los valores de resistividad oscilan entre 200 – 500ohm (tonalidad roja), asociados

a materiales puzolánicos, provenientes de la formación Sencca con bajo contenido de humedad.

Figura 26. Modelo de Sección de Resistividad de la TE_01 a TE_06.


Fuente: Rock&Soil (2018)

3.4.2.4 LINEA TE_07 a TE_13 (Progresivas 1+090 a 1+780)

La Zona A se asocia a materiales con filtraciones de agua continua y la zona B a materiales

como arenas limosas húmedas (Rock&Soil, 2018).

Figura 27. Modelo de Sección de Resistividad de la TE_07 a TE_13.

Fuente: Rock&Soil (2018)

3.4.2.5 LINEA TE_01_02_03_04

La línea está ubicada dentro de la excavación, se observan dos tipos de materiales

predominantes. El tipo A se asocia a materiales saturados (puzolana) con filtraciones de agua

y el tipo B asociado a arenas limosas húmedas (Rock&Soil, 2018).


55
Figura 28. Modelo de Sección de Resistividad de la TE_01_02_03_04.

Fuente: Rock&Soil (2018)

En síntesis, los resultados de los estudios geofísicos (métodos indirectos) determinaron tres

horizontes de materiales con valores de resistividades diferentes que varían entre 20 - 80 ohm.m

para un primer horizonte (arenas limosas con cierto grado de saturación del volcánico Sencca),

150 ohm.m para el segundo horizonte (materiales con grava), 350 - 500 ohm.m para el tercer

horizonte (suelo con grado de compactación elevado). El hecho de encontrar el primer

horizonte con cierto grado de saturación es coherente la presencia del nivel freático medido en

las perforaciones geotécnicas, con un valor alrededor de 5 m de profundidad.

3.5 Estudios de parámetros geológicos e hidrogeológicos

Para conocer mejor el comportamiento del flujo de agua subterránea necesitamos conocer

algunos parámetros geológicos e hidrogeológicos como la capacidad de almacenamiento,

porosidad, litología y permeabilidad.

Los parámetros como la capacidad de almacenamiento y porosidad, fueron tomados de

estudios geológicos previos; la litología y permeabilidad se determinaron en base a la

investigación en campo. La descripción de cada uno de estos parámetros se en el ítem 4.6 del

Capítulo IV.

Los valores de permeabilidad tomados en campo se describen a continuación.

56
3.5.1 Permeabilidad

Se realizó un total de 11 ensayos Lefranc de carga variable. En el sondeo DHA-01 se

realizaron 5 ensayos y en el sondeo DH-01, 6 ensayos.

Tabla 21. Permeabilidades del sondaje DH-01

Código de Permeabilidad k
Desde (m) Hasta (m) Litología
perforación (cm/s)
DHA-01 5.5 5.6 Suelo arenoso 1.40E-04
DHA-01 6.9 7 Suelo arenoso 1.30E-04
DHA-01 10.4 11 Suelo arenoso 1.13E-04
DHA-01 16 17 Suelo arenoso 7.38E-04
Suelo tobáceo
DHA-01
20.1 21 de elevada 2.78E-05
plasticidad

DHA-01 24.1 25 Suelo arenoso 4.85E-04

Fuente: Perfocon S.A.C. (2018)

Tabla 22. Permeabilidades del sondaje DHA-01

Código de Permeabilidad
Desde (m) Hasta (m) Litología
perforación k(cm/s)
Suelo
DHA-01 6.8 7.5 3.14E-03
puzolánico
Suelo
DHA-01 12.2 13 puzolánico 6.61E-04
Suelo
DHA-01 17.4 18 puzolánico 8.51E-03
Suelo
DHA-01 21.5 22 arenoso 3.67E-03
Suelo
DHA-01 24.3 25 arenoso 5.61E-03

Fuente: Perfocon S.A.C. (2018)

57
3.6 Parámetros hidráulicos

3.6.1 Nivel freático

Sobre la base de las perforaciones realizadas para la caracterización geotécnica descritas en

el siguiente capítulo se conoce que la profundidad de agua se sitúa a 5.30 m en el puente de

Atahualpa y 5.00 m en el intercambio a desnivel de Taboada con respecto del nivel terreno

natural.

El Plano N° 8 muestra el mapa conceptual muestra la piezometría del área del modelo, este

fue construido a partir de los niveles piezométricos medidos en las perforaciones realizadas y

de la información disponible de los pozos de agua (ANA, 2017). Un total de 14 puntos de

control piezométrico fueron considerados para la construcción de las isopiezas.

El mapa piezométrico indica que las cargas hidráulicas varían entre 2170 msnm y 2360

msnm, con una dirección del flujo subterráneo principalmente de norte a sur. Las Tablas 13 y

14 muestran la ubicación de los pozos y puntos de control considerados.

La mayor profundidad de excavación en terreno saturado se encuentra en la progresiva

2+020 m cuya cota de rasante es de 2298.6 m.s.n.m. y que tiene la cota de excavación de 2286.2

m.s.n.m. es decir se requiere una excavación de 12.4 m. Debido a que los trabajos de

excavación se sitúan por debajo del nivel freático será necesario la instalación de pozos de

bombeo que hagan posible la excavación en seco y que asegure la estabilidad y seguridad de

la zona de trabajo.

58
Capítulo IV

Caracterización de la zona de estudio

4.1 Situación

El área a drenar está ubicada en el distrito de Sachaca, corresponde al tramo de la Vía

Variante de Uchumayo Tramo II entre las progresivas KM 1+560.00 al KM 2+520. El área

ocupa una superficie de 28337 m2 aproximadamente (ver Plano N°1).

Figura 29. Ubicación del proyecto

Fuente: Rock&Soil (2018)

59
AYACUCHO

JUAN SANTO
LOS

S ATAHUALPA
ALAM
OS

O
VARIANTE UCHUMAY

MADRE DE
DIOS
VA
RIA
NTE
UC
HU
MAY
O

ADA
BO
TA
LE
L
CA
RESERVORIO
AMAZO
NAS

PLANO PROVINCIAL
20km.E
5 20km.E
6 20km.E
7 820km.E 20km.E
9

00km.N

00km.N
84

84
CUSCO

ICA AYACUCHO

Cotahuasi
Orcopampa Caylloma

Sibayo

00km.N

00km.N
Andagua
PUNO

83

83
Machahuay CAYLLOMA
CASTILLA Chivay
Caraveli CONDESUYOS
Tipan
Chuquibamba

Aplao

AREQUIPA

00km.N

00km.N
82

82
CAMANA Arequipa

Camana YARABAMBA

MOQUEGUA
Mollendo
ISLAY

00km.N

00km.N
81

81
20km.E
5 20km.E
6 20km.E
7 820km.E 20km.E
9

PLANO DISTRITAL

CERRO COLORADO YANAHUARA

Huasamayo
SACHACA
VITOR

R
TE
UBICACION DEL

UN
TIABAYA

J.H
UCHUMAYO
PROYECTO El Ramal

SOCABAYA

UNIVERSIDAD NACIONAL DE SAN AGUSTIN


POSTGRADO DE LA FACULTAD DE
La sección tipo de la vía principal considerada en toda la longitud consta de una berma

exterior de 3.0 m y dos carriles de 3.60 m para cada sentido de circulación. El separador central

es de 3.0 m conformado por una barrera de concreto y bermas interiores de 1.20 m.

Figura 30. Sección típica de la vía

Fuente: Rock&Soil (2018)

4.2 Caracterización geomorfológica

4.2.1 Geomorfología Regional:

La región de Arequipa está ubicada a 90 Km. en línea recta del Océano Pacífico, a 16º de

latitud sur, mostrando fuertes contrastes de altitud (desde los más de 6000 m.s.n.m. de la

cumbre del Chachani).

La geomorfología representa las formas topográficas del terreno, originadas por la presencia

de uno o varios agentes como son la erosión, velocidad del viento, gravedad, agua y

sedimentación, o las generadas por la acción de fuerzas interiores como son los volcanes o

sismos.

En el Plano N°2 Plano Geomorfológico Regional podemos distinguir tres unidades

geomorfológicas, las unidades presentadas están dentro de la geomorfología de la ciudad de

Arequipa:

61
 Arco volcánico del Barroso: Esta unidad geomorfológica, comprende los flujos de

lodo y lava, que están asociadas a las erupciones volcánicas y a depósitos por la

actividad post-volcánica de los tres volcanes que son: El Chachani, Misti y Pichu Pichu

(Koasaka y Macedo, 2000). Comprenda la siguiente subunidad:

- Estrato Volcán (Es-v): Se caracteriza por presentar pendientes moderadas (entre

15° y 40°), y un drenaje radial, con quebradas de 10 a 50 m de profundidad.

Formado por flujos de lava del volcán Misti, emplazados entre 50 y 120 mil

años (INGEMMET, 2016).

 Cordillera de Ladera: Es una cadena montañosa con una elongación de Noroeste a

Sureste. Está constituido por roca plutónica, en la parte Sureste, predominan las rocas

metamórficas del Complejo Basal y en la parte Sureste, dentro de la masa, intrusiva, se

encuentran algunos remanentes de rocas sedimentarias. (Vargas, 1970). Comprende las

siguientes subunidades.

- Colina y lomada en roca sedimentaria (RCL-rs): Compuesta por rocas

sedimentarias. Están representadas por colinas y lomadas con pendiente

moderada. Ubicadas al Noroeste de la ciudad de Arequipa.

- Colina y lomada en roca intrusiva (RCL-ri): Compuesta por rocas intrusivas.

Están representadas por colinas y lomadas con pendiente moderada. Ubicadas

al Suroeste de la ciudad de Arequipa.

- Montañas y colinas en roca intrusiva (RMC-ri): Compuesta por rocas intrusivas.

Están representadas por colinas y lomadas con pendiente fuerte. Ubicadas al Sur

de la ciudad de Arequipa.

 Penillanura de Arequipa: Superficie suavemente ondulada de forma groseramente

triangular, comprende entre las localidades de Arequipa, Yura y la confluencia de los

62
ríos Chili y Yura. A su vez está rodeada de cerros altos que forman parte de la Cordillera

de Laderas. (Vargas, 1970). Compuesto por las siguientes subunidades:

- Superficie de flujo piroclásticos (Sfp): Formado por tufos del Volcánico Sencca,

cubren una gran superficie de la ciudad y están en el noroeste en el distrito de

Cerro Colorado.

- Vertiente Glacio-Fluvial (V-gfl): Formado por diversos tipos de materiales, se

encuentra cubriendo partes de la superficie de flujo piroclásticos.

- Terraza Aluvial (T-al): Formado por materiales aluviales recientes.

63
216 000 E 220 000 E 224 000 E 228 000 E 232 000 E
8 192 000 N

8 192 000 N
Es-v
Sfp Es-v

Sfp

RCL-rs 8204
Coordendas UTM

Tramo del proyecto


8 188 000 N

8 188 000 N
Red Vial Nacional

V-gfl

RCL-rs

UNIDADES Y SUBUNIDADES
Area del
RCL-rs proyecto
Sfp UNIDAD
SISTEMA
LITOESTRAT.
8 184 000 N

8 184 000 N
Es-v
del Barroso

Sfp Superficie de flujo

Plenillanura de
Glacial Arequipa V-gfl Vertiente glacio-fluvial

Fluvio-aluvial T-al Terraza aluvial


Ki-gb
RCL-ri

RCL-rs Colina y lomada en roca


sedimentaria
RCL-ri RCL-ri Cordillera de RCL-ri Colina y lomada en roca
8 180 000 N

8 180 000 N
Laderas intrusiva
RCL-ri

Sfp
RMC-ri
V-glf
RCL-ri
RCL-ri roca intrusiva

RCL-ri

RCL-ri
RCL-ri

V-gfl

RCL-ri
RCL-ri
8 176 000 N

8 176 000 N
RMC-ri
Sfp
V-glf
RCL-ri

216 000 E 220 000 E 224 000 E 228 000 E 232 000 E

Escala
PLANO
Metros
TESISTA: MARISOL LLAMOCCA A. UNIVERSIDAD NACIONAL DE SAN AGUSTIN
0 250 500 1000 1500
ASESOR: MSc. LIZARDO HUAMANI H. POSTGRADO DE LA FACULTAD DE
Escala 1/50000
FECHA: SETIEMBRE 2020
FUENTE: ADAPTADO DE INGEMMET

DATUM: WGS84
TESIS:
02
ZONA: 19S
4.2.2 Geomorfología local:

En el Plano N° 03 Plano Geomorfológico Local observamos las siguientes subunidades:

 Superficie de flujo piroclástico (Sfp): Superficie de relieve ondulado ubicado alrededor

de la zona del proyecto.

 Terraza aluvial (T-al): Formado por acumulaciones de material aluvial (bancos de suelo

que son una mezcla grava, arenas y limos dispuestos en estratos) que se encuentran a

lo largo de la vía.

Figura 31. Terraza aluvial en la zona de excavaciones para la vía principal.

Fuente: Propia

65
PLANO GEOMORFOLOGICO LOCAL

224 400 E 224 600 E 224 800 E 225 000 E 225 200 E 225 400 E 225 600 E 225 800 E 226 000 E 226 200 E 226 400 E 226 600 E

8 184 800 N 8 184 800 N

PTE. JUAN SANTOS


ATAHUALPA

T-al 1+700
1+800

8 184 600 N 1+90


0
8 184 600 N
T-al
2+
00
0

2+
100

Sfp
2+
20
0

INTERCAMBIO
A
T-al 2+
30
0 DESNIVEL
8 184 400 N T-al TABOADA 8 184 400 N

T-al 2+
4 00

Sfp 2+
500

Sfp 2+
60
0 FIN DEL
TRAMO II

2+
65
2.8
58

8 184 200 N 8 184 200 N

8 184 000 N 8 184 000 N


224 400 E 224 600 E 224 800 E 225 000 E 225 200 E 225 400 E 225 600 E 225 800 E 226 000 E 226 200 E 226 400 E 226 600 E

0 50 100 200
Metros

Escala 1/750

Coordendas UTM 8204

UNIDAD
SISTEMA
LITOESTRAT. Grillado
Progresiva 1+500
Sfp Superficie de flujo
Curvas de Nivel
Plenillanura de
Fluvio-aluvial Arequipa T-al Terraza aluvial

PLANO
TESISTA: MARISOL LLAMOCCA A. UNIVERSIDAD NACIONAL DE SAN AGUSTIN
ASESOR: MSc. LIZARDO HUAMANI H. POSTGRADO DE LA FACULTAD DE
FECHA: SETIEMBRE 2020
ESCALA: INDICADA

DATUM: WGS84
TESIS:
03
ZONA: 19S
4.3 Caracterización geológica

4.3.1 Geología regional

Se observan rocas y suelos de origen sedimentario, ígneo y metamórfico con edades que

comprenden desde la era precámbrica hasta la era reciente; estos forman las diferentes unidades

litoestratigráficas. En orden de antigüedad tenemos primeramente al complejo basal de la costa

que es un conjunto de rocas metamórficas que se encuentran en la parte noroeste y son de la

era del Precámbrico. En segundo lugar la formación Chocolate, así como la formación Socosani

formada por sedimentos calcáreos y también el grupo Yura formado por cinco miembros:

Puente, Chachíos, Labra, Gramadal y Hualhuani. Ellos en su mayoría pertenecen al Jurásico.

En tercer lugar el material en capas rojas pertenecientes a la formación Huanca en la región

noreste y la formación Sotillo en la zona suroeste. Ambos del Paleógeno. En cuarto lugar los

depósitos de la formación Millo y el volcánico Sencca del Neógeno. Ellos cubren un área

extensa que tras un intenso vulcanismo dio origen a los volcánicos Chila y Barroso. Finalmente

los depósitos recientes de material aluvial y eólico y las rocas intrusivas que están en la

cordillera de Laderas del Cuaternario. Podemos observar estas unidades en el Plano N°4 Plano

Geológico Regional.

Para más alcance se describe cada una de ellas según lo mencionado en el Boletín N°24 del

Cuadrángulo de Arequipa:

a) Del Pre-Cámbrico:
 Complejo Basal de la Costa (Pe-gn, Pre-Cámbrico): Bellido y Narváez (1960) lo

definen como un conjunto de rocas metamórficas de características similares

presentes en la parte NO del Cuadrángulo de Arequipa, formando una faja ancha y

en el límite Sur estas rocas se exponen como una prolongación del afloramiento de

Cerro Verde.

b) Del Mesozoico

67
 Formación Chocolate (Ji-cho, Jurásico Inferior): Vargas (1970) menciona que el

nombre fue dado por W. Jenks (1948) a causa de la cantera de Chocolate situada a

20 km. al NO de Arequipa. Las relaciones estratigráficas de la formación Chocolate

con las rocas del Complejo Basal no están definidas, e infrayace con discordancia

erosional a la Formación Socosani. A 1.5 km de esta cantera se halla otro

afloramiento, separados entre sí, por los tufos pliocénicos. Aquí la secuencia es

predominantemente de origen volcánico, alternan derrames, brechas y tufos, con

capas de calizas marrones ferruginosas y pobremente fosilíferas.

 Formación Socosani (Jm-so, Jurásico Medio): V. Benavides (1962) diferencia tres

miembros dentro de esta formación. La formación Socosani descansa en

discordancia erosional sobre los volcánicos Chocolate y también subyace en igual

condición al miembro Puente del grupo Yura. En la carretera Variante de

Uchumayo, la formación Socosani está de forma alargada de NO a SE, en las rocas

plutónicas, mostrando un metamorfismo de contacto más o menos intenso y algunas

perturbaciones estructurales, como pliegues, fallas y fracturas que han originado el

desprendimiento de bloques, (Vargas, 1970)

 Grupo Yura: V. Benavides (1962) indica que la sección típica mide 2,142 m.,

dividiéndola en cinco miembros en el siguiente orden, de la base al tope, Puente,

Cachíos, Labra, Gramadal y Hualhuani; luego, Guevara (1969) reconoce los

mismos miembros con un grosor de 2,900 m.

- Formación Puente (Js-pu, Jurásico Superior): Se caracteriza porque predominan

las areniscas cuarcíticas de grano muy fino, de colores amarillentos, pardos y

verdosos, que intemperizan a colores pardo y rojo grisáceo, generalmente en

estratos medianos y en ocasiones gruesos, impuros, estratificados con

68
cantidades menores de lutitas negras, carbonáceas. Se le da 600 m de espesor.

(Vargas, 1970)

- Formación Cachíos (Js-ca, Jurásico Superior): V. Benavides (1962) le da el

nombre por llamarse así la quebrada en la que se midió la sección típica; indica

que está constituido por lutitas negras y gris oscuras, tufáceas y carbonáceas;

con intercalaciones menores de lutitas y siltitas grises a beiges, de grano fino,

en capas medianas, fuertemente lenticulares y con un grosor total de 603 m.

- Formación Labra (Ki-la, Cretácico Inferior): Caracterizado por una

secuencia litológica bastante variable, predominando las areniscas cuarcíticas y

cuarcitas sobre las lutitas y limolitas, que a su vez, por ser más duras, sobresalen

en la topografía. (Vargas, 1970)

- Formación Gramadal (Ki-gr, Cretácico Inferior): Esta integrado

esencialmente por bancos de 1 a 3.50 m. de calizas arrecifales de color marrón

a gris oscuro, con gran contenido fosilífero. Se intercalan con las calizas, en

menor proporción, hay estratos de lutitas y limolitas amarillo verdosas hasta

marrón violáceas. En la parte media, se aprecia la presencia de bancos gruesos

de arenisca cuarcítica de color blanquecino. (Vargas, 1970)

- Formación Hualhuani (Ki-hu, Cretácico Inferior): Estas rocas son duras y

compactas formando una topografía escarpada. En localidad de Hualhuani este

Miembro se caracteriza por estar constituido de areniscas cuarcíticas y cuarcitas

de grano fino, de color blanco que por intermperismo se torna rojo amarilleno

claro. (Vargas, 1970)

 Formación Murco (Ki-mu, Cretácico Inferior): Se encuentra formada por lutitas

abigarradas con estratificación laminar y fácilmente deleznables, las cuales se

69
intercalan con capas de areniscas y algunos estratos calcáreos, esto en la zona del

río Pichirigma en Huanca. (Vargas, 1970)

 Formación Acurquina (Kms-ar, Cretácico Medio a Superior): la base de la

Formación encontramos fangolitas y areniscas marrón rojizas que se intercalan con

capas gruesas de calizas, sobre esta zona se presenta margas, conglomerados

calcáreos y capas de areniscas verdosas intercaladas con calizas. (Vargas, 1970)

 Formación Chilcane (Ks-chi, Cretácico Superior): Constituido por yeso es de color

blanco tintes rojizos y verdosos; con una ligera estratificación paralela a la de la

Formación infrayacente, intercalándas con algunas capas delgadas de lodolita roja

y lutitas verdes. (Vargas, 1970)

 Rocas Ígneas Intrusivas (Cretácico Tardío-Terciario Inferior): Forma parte del

Batolito Costanero que forman una cadena montañosa que se extiende desde el

límite Sur, hasta cerca del río Siguas en el Norte.

- Grupo Gabro-diorita: Se encuentra afloramientos en los cerros Gloria y San

Ignacio, al norte del río Vitor, en Huasamayo, en Quishuarani y Cerro Gloria.

En la ciudad de Arequipa se hallan en los alrededores de Tiabaya y Hunter. La

roca es de color gris oscuro y gris verdoso por alteración y algo rojizo por

intemperismo; el grano es mediano a grueso, con cristales de feldespato y

hornblenda, y en algunos casos biotita y cuarzo.

- Tonalita de laderas: Los afloramientos presentan relieves suaves o algo ásperos,

y la roca es de grano grueso, de color blanco rosáceo y gris oscuro por

contaminación con las del grupo gabrodiorita. Los afloramientos presentan

relieves suaves o algo ásperos, y la roca es de grano grueso, de color blanco

rosáceo y gris oscuro.

70
- Granodiorita: Se encuentran en la parte meridional del cuadrángulo, en los

cerros Calderas y San Ignacio. Presenta grandes bloques rectangulares y por

intemperismo una ligera disyunción catafilar, lo que ha originado bloques

desprendidos en el primer caso, y esferas imperfectas, en el segundo.

c) Del Cenozoico

 Formación Huanca (Ti-h, Terciario Inferior): Esta Formación en la parte inferior de su

afloramiento en la quebrada del río Pichirigma, afluente del río Siguas, está formada

por areniscas arcósicas marrón rojizas de grano mediano a grueso, elementos de cuarzo

angulares a subangulares, tiene poca compacidad. La mayor parte del afloramiento está

conformado por conglomerados, constituidos principalmente por elementos de

areniscas y cuarcitas (derivados posiblemente del grupo Yura y de la formación Murco),

intrusivos muy alterados y en menor proporción calizas y cherts que tal vez provienen

de la formación Arcurquina.

 Formación Sotillo (Ti-s, Pre- Oligoceno): Expuesta en el sector Suroeste del

cuadrángulo, con afloramientos en ambas laderas del valle de Vítor, a la altura del

pueblo de Sotillo, y también en el flanco occidental del batolito. Las relaciones

estratigráficas son las siguientes: descansa sobre una superficie de erosión en las rocas

intrusivas y está cubierta en discordancia erosional por la Formación Millo. Tiene dos

miembros Superior e Inferior.

 Grupo Tacaza (Tm-t, Terciario Medio): Conjunto litológico, predominantemente

volcánico. Presenta tres unidades: la unidad inferior sedimentaria (conglomerados

intercalados con arenisca tufácea, la unidad media tufácea (tufos brechoides y

compactos) y unidad superior lávica (derrames gris oscuro).

 Formación Millo (Ts-mi, Plioceno Inferior a Medio): Conjunto estratigráfico que

consiste de conglomerados, moderadamente consolidados; areniscas tufáceas, grisáceas

71
de grano medio a fino y lentes de tufos que forman capas lenticulares con lapilli y

pomez.

 Volcánico Sencca (Tp-vse, Plioceno Medio a Superior): En la región de Arequipa cubre

grandes áreas desde el límite Norte del cuadrángulo (Pampa de Arrieros), en el Sur una

faja angosta entre el río Yura y las lavas de grupo Barroso, hasta llegar al denominado

“triángulo de Arequipa” (Yura, Arequipa y Quishuarani, W. Jenks, 1948).

Compuesto por tufos son de composición dacítica o riolítica, con granos de cuarzo,

feldespatos y lamelas de biotita muy visibles. Además contienen fragmentos de pómez

y lavas, que pueden ser redondeados o angulosos y de tamaño variable.

De consistencia compacta, con cohesión aunque también hay de tipo poco consistente

y fácilmente desmenuzable. Se encuentran dando lugar a columnas cortadas por planos

horizontales. Se distinguen dos niveles según el color, uno superior, rosado a marrón

rojizo, y otro inferior, gris claro a blanco. También alternados con los tufos se encuentran

capas de arenas, gravas y tufos retrabajados con poca consistencia.

En la zona del proyecto se presenta como un tufo rosado de poca consistencia e

inestabilidad debido a la presencia de agua y que se encuentra cubierto por depósitos

aluviales en varios sectores.

 Grupo Barroso: En Arequipa se encuentran dos de sus unidades:

- Volcánico Chila (Q-vchi, Plio- Pleistocénica): Representa la unidad inferior del

grupo. En el cañón del río Chili (Charcani) las coladas se prolongan cerca del

distrito de Miraflores (falda occidental del Misti).

- Volcánico Barroso (Q-vba, Plio- Pleistocénica): Está distribuida desde la esquina

Noreste del cuadrángulo, en las faldas del cerro Las Minas o Suni, hasta cerca del

aeropuerto de Arequipa, esta unidad forma la mayor parte de la cadena del Nevado

Chachani, incluyendo a los aparatos volcánicos más recientes de los cerros

72
Cortaderas y El Tolar. La litología consiste de lavas andesíticas en bancos gruesos,

gris oscuras y gris claras, que por intemperismo toman colores rojizos y marrones.

 Flujos de Barro (Q-fb, Neo- pleistocénica): Los depósitos de flujos de barro ocupan la

esquina Sureste del cuadrángulo (Socayaba, Urbanización Hunter, y más al Norte, en

la falda occidental del Misti). Estos, al tener poca compactación, han sido erosionados

dando lugar a algunas quebradas profundas. En general, consisten de una masa formada

por fragmentos y bloques angulares de rocas de origen volcánico (tufos y lavas), dentro

de una matríz areno-tufácea.

 Depósitos piroclásticos recientes (Holoceno): Una parte se encuentra en las faldas

orientales del nevado Chachani (cuadrángulo Characato), se observa que estos

depósitos están sobre los fluvioglaciares y debajo de gravas aluviales recientes. Esta

unidad está formada principalmente por capas de lapilli de color amarillo y

blanquecino, que contienen pequeños fragmentos de pómez, lavas y algunas escorias y

bombas, cementados por ceniza en parte arenosa y poco consolidada. Alternando con

dichas capas de hallan otras más delgadas de arena de grano grueso y de color gris

oscuro, que parece haber sido acarreada por el agua.

 Depósitos aluviales recientes (Qh-al, Holoceno): Constituyen el suelo de la mayor parte

de las llanuras y depresiones. Podemos verlos ocupando el lecho y laderas de los

diferentes riachuelos, quebradas y a lo largo de los valles principales, como el Vítor,

Chili y Yura. El material consiste principalmente de gravas o conglomerados poco

consolidados, arenas y limos.

73
Figura 32. Columna estratigráfica

Fuente: Adaptado de INGEMMET

74
216 000 E 220 000 E 224 000 E 228 000 E 232 000 E
8 192 000 N

8 192 000 N
Qh-al

Qpl-mm/ct
Qpl-ma/an

Qpl-mm/ct

Qh-al
Qpl-cu Coordendas UTM 8204

Qpl-pi
Tramo del proyecto
8 188 000 N

8 188 000 N
Red Vial Nacional

Qpl-pi

Qh-al
Qpl-cu

Qh-al

Area del
proyecto UNIDAD
Qpl-cu Qpl-cu ERA SISTEMA SERIE
LITOESTRAT.
8 184 000 N

8 184 000 N
PE-34C
Holocena Qh-al Bloques, gravas y
Aluviales arenas
Qh-al
Flujos Qpl-pi

Pleistocena
Qpl-mm/ct

Cenozoica
Cuaternario Superior Lapilli y cenizas

Qpl-cu Arena y conglomerados


Ki-gb Ki-gb
Ki-gb
Qpl-pi/bxa Bloques, arenas y
Media avalancha del
Qpl-pi/bxa cenizas
Pichu
Disc. Erosional
Qpl-pi/bxa Qh-al Tobas soldadas y no
Ki-gb Ki-gb Ki-gb Pliocena Tp-vse
8 180 000 N

8 180 000 N
Ki-gb
soldadas
Discontinuidad
Ki-gb Superunidad ks-Ti/gb Granodiorita. Segmento

Mesozoica
Qh-al Tiabaya Toquepala
Ki-to/di Superior
Qpl-pi/bxa Qpl-pi/bxa
Superunidad Ilo Ki-to/di
Tonalita-diorita.
Ki-to/di Segmento Toquepala
Ki-Gb
Gabrodiorita. Segmento
Ki-gb tempranas
Ki-gb
Arequipa
Qh-al
Ki-to/di
Ki-gb
8 176 000 N

8 176 000 N
Ks-ti/gd
Qpl-pi/bxa
Qh-al
Ki-gb

216 000 E 220 000 E 224 000 E 228 000 E 232 000 E

Escala
PLANO
Metros
TESISTA: MARISOL LLAMOCCA A. UNIVERSIDAD NACIONAL DE SAN AGUSTIN
0 250 500 1000 1500
ASESOR: MSc. LIZARDO HUAMANI H. POSTGRADO DE LA FACULTAD DE
Escala 1/50000
FECHA: SETIEMBRE 2020
FUENTE: ADAPTADO DE INGEMMET

DATUM: WGS84
TESIS:
04
ZONA: 19S
4.3.2 Geología local

La zona del proyecto está sobre la Unidad del Volcánico Sencca cubiertos por depósitos

aluviales pertenecientes al Aluvial Reciente como se muestra en el Plano N°5 Plano Geológico

Local.

 Volcánico Sencca (Tp-vse)

En el puente Juan Santos Atahualpa se presenta un tufo areno-limoso de color rosado

con textura fina, cuya compacidad aumenta con la profundidad, desmoronable a la presión de

los dedos pero que en presencia de humedad y a mayores profundidades es más denso y

consolidado. Esporádicamente está intercalado por gravillas de 5 mm. No presenta plasticidad

ni límites de consistencia.

Según Aguilar (1991) el material de esta zona está compuesto por tufos estratificados

que siguen la topografía del terreno, adopta un color rosado a marrón rojizo, tiene biotita y

fragmentos de pómez y lava.

Figura 33. Material puzolánico en zona de excavaciones

Fuente: Propia

76
 Aluvial reciente (Qr-a):

En la zona del Intercambio de Taboada se encontraron sectores formados por arena gruesa

a media de color gris oscuro y lentes de material areno limoso con cierta plasticidad. A más

profundidad el material cambia a arena media a fina de color gris oscuro y trazas de limos de

color gris claro y para luego nuevamente aparecer una mezcla de arenas gruesas a media con

gravillas y gravas esporádicas. A los 18 metros cambia a un suelo tobáceo de color beige.

Figura 34. Material aluvial compuesto por arenas

Fuente: Propia

77
ESCALA: 1/750

224 400 E 224 600 E 224 800 E 225 000 E 225 200 E 225 400 E 225 600 E 225 800 E 226 000 E 226 200 E 226 400 E 226 600 E

8 184 800 N 8 184 800 N

Tp-vse

Tp-vse
PTE. JUAN SANTOS
Qh-al ATAHUALPA
1+800
1+700

8 184 600 N 1+90 8 184 600 N


0

Qh-al 2+
00
0

Qh-al
2+
100

2+
20
Tp-vse ST_01
0

ST_03

INTERCAMBIO
2+
30
A
0 DESNIVEL
8 184 400 N Qh-al Qh-al DH_01
TABOADA 8 184 400 N
Qh-al
Qh-al 2+
40
0

2+
50
0

2+
600 FIN DEL
TRAMO II

2+
652
.85
8

8 184 200 N 8 184 200 N

Tp-vse

Metros
8 184 000 N
0 50 100 200 300 8 184 000 N

224 400 E 224 600 E 224 800 E 225 000 E 225 200 E 225 400 E 225 600 E 225 800 E 226 000 E 226 200 E 226 400 E 226 600 E

NW
PUENTE DE SAN JUAN
ESCALA: 1/500
2,320 ATAHUALPA INTERCAMBIO A DESNIVEL
Qh-al DH_01
2,310 TABOADA - ROTONDA
Qh-al Qh-al Tp-vse
2,300
2,290
2,280 Tp-vse
Tp-vse
1+550.00

1+570.00

1+590.00

1+610.00

1+630.00

1+650.00

1+670.00

1+690.00

1+710.00

1+730.00

1+750.00

1+770.00

1+790.00

1+810.00

1+830.00

1+850.00

1+870.00

1+890.00

1+910.00

1+930.00

1+950.00

1+970.00

1+990.00

2+010.00

2+030.00

2+050.00

2+070.00

2+090.00

2+110.00

2+130.00

2+150.00

2+170.00

2+190.00

2+210.00

2+230.00

2+250.00

2+270.00

2+290.00

2+310.00

2+330.00

2+350.00

2+370.00

2+390.00

2+410.00

2+430.00

2+450.00

2+470.00

2+490.00
ERA SISTEMA SERIE

Coordendas UTM 8204


CUATERNARIO HOLOGENA Qh-al Bloques, Gravas y arenas.
CENOZOICO

PLANO
Grillado
Discont. Erosional
TESISTA: MARISOL LLAMOCCA A. UNIVERSIDAD NACIONAL DE SAN AGUSTIN
Progresiva 1+500 PLIOCENO
Discontinuidad
Tp-vse ASESOR: MSc. LIZARDO HUAMANI H. POSTGRADO DE LA FACULTAD DE
Curvas de Nivel
Perforaciones DH_01
FECHA: SETIEMBRE 2020
FUENTE: ADAPTADO DE ROCK&SOIL

DATUM: WGS84
TESIS:
05
ZONA: 19S
4.3.3 Estratigrafía local

Según los estudios realizados se describe a continuación los perfiles litológicos

interceptados con las perforaciones.

 Para el Puente Juan Santos Atahualpa

- 0.0 – 2.30 m. (SM) Niveles de suelos residuales conformado por arenas limosas

que son acompañadas por gravas de compacidad medianamente densa que se

presentan de forma esporádica pobremente graduadas.

- 2.30 – 18 m. (SP-SM) Presenta niveles de suelo areno-limoso de color rosado

con una textura fina compuesto de arenas finas mal graduadas no presenta

plasticidad. Su estimación visual 5% de arenas gruesas, 85% de arenas medias

y finas y un 10% de finos pertenecientes al volcánico Sencca conocida como la

Puzolana, en el tramo de 16.05 – 16.50 se muestra una textura más gruesa

predominando las arenas gruesas. Se registró que el nivel freático

correspondiente al acuífero libre, se sitúa a una profundidad de 5.30 m

aproximadamente desde el nivel del terreno natural.

- 18.00 – 21.35 m. (SP) Se encuentra representado por niveles de suelos arenosos

de color gris con tonalidad rojiza, arenas pobremente graduada, no presenta

cohesión; su estimación visual es de 95% de arenas y 5% de finos. Suelo

puzolánico del volcánico Sencca.

- 21.35 m. (SM) Suelo areno limoso de color rosado, de textura fina sin cohesión,

presenta arenas gruesas de forma esporádica; tiene una estimación visual de

70% de arenas medias a finas y un 30% de suelos finos del volcánico Sencca.

79
Figura 35. Estratigrafía Puente Juan Santos Atahualpa (Perforación DHA-01)

Fuente: Propia

80
 Para el intercambio a desnivel Taboada

- 0 - 2.45 m. (SP-SM) se aprecian niveles de arenas limosas de compacidad

medianamente densa a suelta, color pardo amarillento con clastos esporádicos

de 0.5 - 2 cm.

- 2.45 – 9.00 m. (SP) se tiene niveles de suelo arenoso con presencia de gravas,

las arenas predominantes son de granos fino a grueso, también presenta lentes

de arena pobremente graduadas con presencia de limos. La arena presenta

gravillas heterogéneas de formas sub-redondeadas a sub-angulares con TM= 1

½”. Se registró que el nivel freático se encuentra a una profundidad de 5.00 m

aproximadamente con respecto del nivel del terreno natural.

- 9.00 – 18.00 m. (SP-SM) Se evidencian suelos limo-arenosos de color gris a

pardo oscuro, del mismo modo se encuentran gravillas de manera esporádica de

pómez alteradas, el suelo arenoso está compuesto por arenas de grano fino a

grueso. Entre los 11.00 a 11.95 m aparece un suelo limo-arenoso de color beige

probablemente procedente de la alteración de roca Tobácea del volcánico

Sencca.

- 18.00 – 25.00 m. (SM) Suelo tobáceo de color beige, de alta plasticidad

conformado mayormente por limos y en baja proporción arenas finas y medias.

81
Figura 36. Estratigrafía Intercambio a desnivel Taboada (Perforación DH-01)

Fuente: Propia

82
4.4 Caracterización geotécnica

Según los estudios e investigación geotécnica y geofísica se realizó la caracterización

geotécnica:

4.4.1 Para el Puente Juan Santos Atahualpa

El material presente en la zona del puente Atahualpa es granular con alto porcentaje de

arenas limosas de naturaleza volcánica de textura media a fina y gravas en menor

porcentaje.

El material predominante areno-limoso tiene una clasificación SUCS de SM, de color

rosado a gris rojizo, compacto, no presenta plasticidad y conforme aumenta la profundidad

la textura del material es más fina.

Se muestran a continuación algunos parámetros geotécnicos de una muestra tomada de

11.30 a 13.00 m de profundidad:

Tabla 23. Parámetros geotécnicos M-01 en puente Juan Santos Atahualpa

Grav Arenas Finos SUCS Humedad LL(%) LP(%) ∅ (°) C


as(%) (%) (%) (%) (kg/cm2)
2.4 88.4 9.3 SP- 30.26 NP NP 28.38 4.12
SM

Fuente: Perfocon S.A.C. (2018)

4.4.1.1 Ensayo de SPT

El ensayo de penetración estándar se aplica normalmente en suelos granulares y nos

determina de manera indirecta la resistencia del terreno. A partir del número de golpes podemos

conocer la consistencia del material.

En el estribo derecho del puente Atahualpa se realizaron ensayos de SPT a diferentes

profundidades como se muestra a continuación:

83
Tabla 24. Estados de compacidad y consistencia.

Índice de
resistencia a la Designación
penetración N
≤4 Suelto
5a8 Poco compacto
Medianamente
9 a 18
compacto
19 a 40 Compacto
>40 Muy compacto

Fuente: Perfocon S.A.C. (2018)

Tabla 25. Resumen de los N corregidos en el estribo del puente Juan Santos Atahualpa.

(N1)60
Profundidad N campo
(CORR)
12.20 a 12.65m 11 10
13.70 a 14.15m 11 10
15.20 a 15.65m 30 21
16.70 a 17.15m 53 30
19.20 a 19.65m 7 6
20.70 a 21.15m 9 7
23.00 a 23.27m Rechazo Rechazo

Fuente: Perfocon S.A.C. (2018)

4.4.2 Para el intercambio a desnivel Taboada – Rotonda

Predomina el suelo areno-limoso con una clasificación SUCS de SP-SM hasta una

profundidad de 21.35 m. Este no presenta límites de plasticidad y tiene un contenido de finos

entre 5 y 10%, el suelo se presenta estratificado y en estado compacto.

A mayor profundidad el suelo se clasifica como SP y se torna de color rojizo (suelo

puzolánico).

84
En la Tabla 26 se muestran los valores obtenidos en muestras a diferentes profundidades:

Tabla 26. Ensayos en muestras a diferentes profundidades.

Grava Arena Finos Humedad


Muestra Profundidad SUCS LL LP
(%) (%) (%) (%)
M-1 6.00-6.70m 7.1 87 5.3 SP-SM NP NP -
M-2 7.00-7.90m 1 89 10 SP-SM NP NP 17.66
M-3 11.00-11.95m 0.2 99.1 0.7 SP NP NP 31.65
M-4 14.30-15.60 2.2 91 6.8 SP-SM NP NP 6.13
M-5 18.55-21.00 8.5 26.7 64.8 MH 62.13 48 64.44

Fuente: Perfocon S.A.C. (2018)

4.5 Caracterización hidrometereológica

4.5.1 Análisis hidrológico de precipitaciones pluviales

Existen dos estaciones cercanas al área de estudio, la estación de la Pampilla y la estación

de Huasacache.

Figura 37. Estaciones pluviométricas de la zona de estudio

Fuente: Adaptado de SENAMHI (2020)

85
Tabla 27. Ubicación de las estaciones pluviométricas

Estación la Pampilla
Latitud 16°24´18.42”
Longitud 71°31´24.02”
Altitud 2365 msnm
Estación Huasacache
Latitud 16°27´27.95”
Longitud 71°33´58.96”
Altitud 2207 msnm

Fuente: Adaptado de Marín (2017)

A continuación se muestra las precipitaciones totales mensuales para las estaciones antes

indicadas.

Tabla 28. Análisis de la precipitación total mensual estación la Pampilla.

Precipitación total mensual entre 1997 al 2015 (mm)


MES ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC
MAX 106.41 86.51 46.3 19.31 2.8 0.8 6.02 12.4 4.5 0.6 0.1 17.1
PROMEDIO 29.321 7.869 4.095 1.732 0.206 0.206 0.638 0.691 0.29 0.086 0.009 2.812
MIN 1.9 0 0 0 0 0 0 0 0 0 0 0

Fuente: SENAMHI. Citado en Marín (2017)

Tabla 29. Análisis de la precipitación total mensual estación Huasacache.

Precipitación total mensual entre 2012 al 2015 (mm)


MES ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC
MAX 19.9 68.2 23.7 5.9 0.3 0.1 0 0 0 0 0 0.2
PROMEDIO 6.78 20.16 7.92 1.18 0.06 0.02 0 0 0 0 0 0.04
MIN 0 0 0 0 0 0 0 0 0 0 0 0

Fuente: SENAMHI. Citado en Marín (2017)

En el modelamiento no se considera recargas debido principalmente al periodo corto en que

se analiza el comportamiento del acuífero y también, a que las excavaciones se realizarían entre

los meses de mayo, junio y julio meses correspondientes a la época de estiaje.

86
4.6 Caracterización hidrogeológica

A partir de la geología regional y la geología local se desarrolló la geología de la zona de

estudio mostrado en el Plano N°6 en la que se pueden diferenciar diferentes tipos de depósitos

dentro del área del modelo.

Los depósitos encontrados en el área de estudio se han agrupado en tres Unidades

Hidrogeológicas (UH) como se muestra en el Plano N°7 Plano Hidrogeológico. La definición

de las UH se realizó sobre la base de la litología y de los parámetros geológicos e hidráulicos

obtenidos a partir de las pruebas realizadas en la zona de estudio y datos registrados en otros

estudios con litologías similares. Las UHs se describen a continuación.

4.6.1 UH Detrítica

4.6.1.1 Características geológicas

El acuífero detrítico está conformado por los depósitos del Aluvial Reciente (Q-al) que son

una mezcla de material limoarenoso con grava. El acuífero que forma este material se clasifica

como un acuífero poroso sedimentario de material no consolidado. (ANA, 2018)

4.6.1.2 Parámetros hidrogeológicos

 Permeabilidad

Para determinar la permeabilidad del acuífero se analizaron los resultados obtenidos de los

ensayos de permeabilidad en la perforación DH-01 mostrado en la Tabla 30 donde observamos

que los valores de permeabilidad van desde 0.02 m/d a 0.63 m/d.

Aguilar (1991) en su estudio de microzonificación sísmica de Arequipa menciona que la

conductividad hidráulica del Aluvial ubicado en Acequia Alta está en el rango de 0.1 a 15 m/d.

También en el ensayo de bombeo realizado en la Urb. Angélica Palacios en Sachaca (Tabla 30)

se obtuvo una conductividad hidráulica de 7.01 m/d (8.12𝑥10−3 cm/s).

87
Tabla 30. Ensayo de Bombeo. Urb. Angélica Palacios Sachaca

Ubicación Tiempo de Conductividad


del pozo recuperación (s) hidráulica (cm/s)

Sachaca 19630 8.12𝑥10−3

Fuente: Tupa (1989)

Al comparar las permeabilidades del ensayo in situ con las obtenidas de otras fuentes

notamos una diferencia esto debido a que la unidad litológica de la zona del proyecto

muestra un mayor contenido de materiales finos en su matriz, que conlleva a una menor

permeabilidad.

En este sentido, se define que la permeabilidad de esta UH varía entre 0.1 m/d y 7.01

m/d. Particularmente para la zona del proyecto, tras la calibración del modelo numérico,

para esta unidad hidrogeológica, se asignó una conductividad hidráulica de 0.25 m/d.

 Almacenamiento específico

Al tratarse de un acuífero libre, un coeficiente de almacenamiento aceptable es de 0.05,

que considerando un espesor promedio de materiales granulares entre 80 -100 m, brinda un

almacenamiento específico de 5x 10-4 m-1. Los valores referenciales de coeficiente de

almacenamiento se muestran en la Tabla 31.

Tabla 31. Coeficiente de almacenamiento

Material Permeable Acuífero S (Valores medios)


Libre 2x10-2
Kárstico: Calizas y dolomías jurásicas Semiconfinado 5x10-2
Confinado 5x10-5
Libre 2x10-2 a 6x10-2
Calizas y dolomías cretácicas y terciarias Semiconfinado 10-3 a 5x10-4
Confinado 10-4 a 5x10-5
Libre 5x10-2 a 15x10-2
Material intergranular gravas y arenas Semiconfinado 10-3
Confinado 10-4

Fuente: Tomado de Villanueva e Iglesias (1984). Mencionado por el ANA (2018)

88
4.6.2 UH Volcánico

4.6.2.1 Características geológicas

El acuífero es de origen volcánico y está conformado por los materiales del Volcánico

Sencca (Tp-vse), formado por tufo limoarenoso poco compacto con intercalaciones de grava.

El acuífero que forma este material se denomina como un acuífero de depósitos fufo arenosos.

(ANA, 2018)

4.6.2.2 Parámetros hidrogeológicos

 Permeabilidad

Los valores de conductividad hidráulica obtenidos en la zona del proyecto están en el rango

de 0.51 m/d (5.65𝑥10−4 cm/s) a 7.35 m/d (8.51𝑥10−3 cm/s).

El ANA (2018) mediante ensayos de bombeo en este tipo de material obtuvo los siguientes

valores de permeabilidad:

Tabla 32. Permeabilidad de la Formación Añashuayco

Permeabilidad (KX10-4)
Pozo IRHS Geología Ubicación Descripción Descenso Recuperación
(m/d) (m/d)
Fm.
04/01/17-21 Sachaca Tufo arenoso brechoides 14.58 7.86
Añashuayco

Fuente: ANA (2018)

En ensayos de bombeos realizados en Semi-Rural Pachacútec se obtuvo una conductividad

hidráulica de 0.48 m/d (5.55𝑥10−4 cm/s) y en Mariscal Castilla, de 0.51 m/d (5.90𝑥10−4

cm/s). (Tupa, 1989).

Acorde a los datos descrito anteriormente, para esta UH se define una permeabilidad que

varía entre 0.5 m/d y 14.6 m/d.

Tras la calibrac1ión del modelo numérico se asignó como conductividad hidráulica

representativa para esta UH el valor de 0.85 m/d.

89
 Almacenamiento específico

Similar a la UH Detrítica, en el proceso de calibración del modelo se definió un valor de 5.0

x 10-4 m-1.

4.6.3 UH Intrusiva

4.6.3.1 Características geológicas

La unidad hidrogeológica intrusiva está formada por la Granodiorita de Tiabaya (Kti-gb),

compuesto por granodioritas moderadamente resistentes que tienen el comportamiento de un

acuífugo (no almacenan ni transmiten agua) formando una barrera hidráulica impermeable.

(ANA, 2018). En los primeros metros se encuentra fisurada por lo que para el modelo se

considerará como un acuitardo.

4.6.3.2 Parámetros hidrogeológicos

 Permeabilidad

Se tiene permeabilidades bajas entre 0.01 m/d (1.16𝑥10−5 cm/s) y 1𝑥10−5 m/d (1.16𝑥10−8

cm/s) de acuerdo al grado de fracturación de la roca (Custodio y Llamas, 1996). Dada la poca

representatividad de la UH en el proyecto se definió como valor representativo de 0.08 m/d

para la conductividad hidráulica.

 Almacenamiento específico

Debido a que los primeros metros de esta unidad hidrogeológica se muestra fracturada

debido a los procesos de meteorización se asignó un valor de 1x10-4 m-1, y se validó en el

proceso de calibración del modelo numérico.

90
223 000 E

224 000 E

225 000 E

226 000 E

227 000 E

228 000 E

229 000 E
Qpl-pi

Tp-vs

Qh-al Coordendas UTM 230000

Tramo del proyecto


Red Vial Nacional

Tp-vs
UNIDAD
LITOESTRAT.
Flujos
Qpl-pi

Disc. Erosional
Qh-al
Aluviales
Disc. Erosional
Tp-vse
Tp-vs
Ki-Gb
tempranas

PE-34C

Qh-al

Ki-gb

PLANO
TESISTA: MARISOL LLAMOCCA A. UNIVERSIDAD NACIONAL DE SAN AGUSTIN
ASESOR: MSc. LIZARDO HUAMANI H. POSTGRADO DE LA FACULTAD DE
0 250 500
Metros

Escala 1/25000
1000
FECHA: SETIEMBRE 2020
ESCALA: INDICADA

DATUM: WGS84
TESIS:
06
ZONA: 19S
223 000 E

224 000 E

225 000 E

226 000 E

227 000 E

228 000 E

229 000 E
Qr-e Q-aaa

2340

LEYENDA:
Tipo:

Qr-a Isopiezas

2320 Curvas de Nivel

Tp-vs Grillado

U. H. Intrusivo

2280

Qr-au 2260

Kti-gb

Qr-a
2220 Qr-au
U. H. Intrusivo

Qr-e

2180

TESISTA: MARISOL LLAMOCCA A. UNIVERSIDAD NACIONAL DE SAN AGUSTIN


ASESOR: MSC.LIZARDO HUAMANI H. POSTGRADO DE LA FACULTAD DE
0 250 500
Metros

Escala 1/25000
1000
FECHA: SETIEMBRE 2020
ESCALA: INDICADA
TESIS:
07
DATUM: WGS84 ZONA: 19S
4.7 Modelo conceptual

A partir de la caracterización de la zona de estudio se definió el modelo hidrogeológico

conceptual. El acuífero está formado por un medio poroso no consolidado que tiene el

comportamiento de un acuífero libre. Está compuesto por tres unidades hidrogeológicas. La

UH de origen detrítico que está compuesto por el acuífero de depósitos aluviales. La UH de

origen volcánico que es parte del acuífero de la Formación Añashuayco. Y la UH intrusiva que

es un acuitardo derivado del acuífugo de la Super Unidad Tiabaya debido su fracturamiento.

Las conductividades hidráulicas se obtuvieron inicialmente de ensayos de permeabilidad y

referencias bibliográficas de ensayos de bombeo.

El nivel de cota piezométrica estática medido en los pozos varía entre 2190 a 2360 msnm

con la dirección del flujo de NW a SW en el lado occidental y de NE a SE en el lado oriental.

El área de estudio se encuentra a una altitud entre los 2200 a 2400 msnm con una

precipitación aproximadamente de 200 mm/año (ANA, 2018). Esta zona no se presenta notable

escorrentía y además por el tipo de modelamiento que es de periodo corto es que no se

considerarán los efectos de recarga producto de las precipitaciones.

En el Plano N° 8 Modelo Conceptual observamos lo descrito.

93
UH

Aluvial Reciente (Q-al) que son


una mezcla de material
limoarenoso con grava.
2340
Formado por los materiales del

formado por tufo limoarenoso


poco compacto con intercalaciones
de grava.
Formada por la Granodiorita de
Tiabaya (Kti-gb), compuesto por
Intrusivo granodioritas moderadamente
2320 resistentes.

2320 PERMEABILIDADES:

UH K (m/d)

0.1 - 7.01

2280 0.5-14.6

Intrusivo 0.00001-0.01
2300

2260

2280
LEYENDA:
Tipo:
2260
2220 Curvas de Nivel

2240
Sentido del flujo
2220 Curva de nivel
Pozos

2200
UH Intrusiva
2180 Red vial nacional

2180

0 250 500 1000


Metros

Escala 1/25000

TESISTA: ING. MARISOL LLAMOCCA A.


UNIVERSIDAD NACIONAL DE SAN AGUSTIN
ASESOR: MSC.LIZARDO HUAMANI H. POSTGRADO DE LA FACULTAD DE
FECHA: SETIEMBRE 2020

ESCALA: INDICADA

DATUM: WGS84 ZONA: 19S


TESIS:

DE UCHUMAYO TRAMO II-AREQUIPA

MODELO CONCEPTUAL
08
Capítulo V

Modelamiento numérico y diseño

5.1 Construcción de modelo numérico

La modelación numérica representa matemáticamente el comportamiento del sistema del

modelo de funcionamiento hidrodinámico, es decir, del modelo conceptual establecido. En este

caso se ha utilizado el código de elementos finitos en tres dimensiones FEFLOW 6.2 (DHI-

WASY GMBH, 2017), el cual permite una adaptación a la forma del dominio, a los distintos

materiales geológicos y a las zonas de mayor interés.

Los modelos numéricos de flujo, generalmente tienen 04 etapas; (1) modelación, (2)

calibración, (3) verificación y (4) validación. En proyectos de obras civiles con interacción del

agua subterránea, los modelos numéricos son utilizados como herramientas que brindan

soporte para definición del diseño del sistema de desaguado, por lo que, en esta primera fase,

se realiza las etapas (1) y (2) de modelación y calibración. Las etapas de verificación y

validación (3) y (4) se realizan una vez se ponga en operación el sistema de drenaje definido,

con la data real de los pozos de bombeo (caudales, eficiencia, radio de influencia) y los

descensos medidos en los piezómetros de control.

El modelo numérico calibrado se constituye como la herramienta fundamental sobre la cual

se soporta el análisis de alternativas de desaguado para optimizar número de pozos, caudales

de bombeo a drenar, así como las estimaciones de la interacción del agua con la infraestructura.

5.1.1 Geometría

El dominio del modelo ha sido adaptado a la piezometría conceptual definida. Los bordes

este y oeste se han seleccionado de manera que funcionen como contornos impermeables

(perpendiculares a la dirección de flujo subterráneo). A diferencia los contornos norte y sur

coinciden con bordes de carga hidráulica impuesta, siendo estas 2360 msnm y 2170 msnm

respectivamente. Los límites del dominio se pueden observar en el Plano N° 9.

95
LEYENDA:
Tipo:
Curvas de Nivel

Grillado

0 250 500 1000


Metros

Escala 1/25000

TESISTA: ING. MARISOL LLAMOCCA A.


UNIVERSIDAD NACIONAL DE SAN AGUSTIN
ASESOR: MSC.LIZARDO HUAMANI H. POSTGRADO DE LA FACULTAD DE
FECHA: SETIEMBRE 2020

ESCALA: INDICADA

DATUM: WGS84 ZONA: 19S


TESIS:

DE UCHUMAYO TRAMO II-AREQUIPA


09
El dominio ha sido discretizado en una malla de elementos finitos, la cual consta de 9,980

elementos con un total de 5,046 nodos. La malla ha sido adaptada y refinada en aquellos

entornos para los cuáles se esperan mayores gradientes hidráulicos y se requiere de mayor

precisión en la solución, tales como los pozos de bombeo o la huella del carril central de la vía.

El tamaño de los elementos en las zonas más refinadas es del orden de 0.2 m. mientras en las

partes más alejadas de las zonas de interés, las dimensiones son aproximadamente 250 m.

La Figura 38 muestra el mallado de elementos finitos y el refinamiento hecho en las zonas

de interés.

Figura 38. Muestra el mallado de elementos finitos y el refinamiento

Fuente: Propia

5.1.2 Parámetros hidrogeológicos

La variación espacial de los parámetros hidrogeológicos iniciales se ha representado

mediante zonas de materiales los cuales se han definido tomando como base los ensayos

realizados durante la campaña de Geotecnia y los valores típicos descritos en la bibliografía

disponible. Durante el proceso de calibración estos valores iniciales han sido ajustados de tal

97
manera el modelo sea capaz de reproducir los niveles piezométricos medidos en campo (ver

Tabla 33).

Tabla 33. Parámetros hidráulicos finales

Conductividad
Unidad Almacenamiento
Descripción Hidráulica
Hidrogeológica específico (m-1)
(m/d)
Arenas y limos con presencia
UH Detrítico 0.25 5.0 x10-4
de gravas
UH Volcánico Puzolana 0.85 5.0 x10-4

UH Intrusiva Roca intrusiva 0.08 1.0 x10-4

Fuente: Propia

5.1.3 Condiciones de contorno

El contorno norte y sur se modelaron con una condición de contorno fija (Dirichlet) de nivel

prescrito igual a la piezometría calculada con los niveles de agua medidos en pozos de

observación. En el contorno norte se ha utilizado un valor de 2360 msnm y en el contorno sur

un valor de 2170 msnm.

Condición de flujo nulo (Neumann) fue definido para el contorno este y oeste. Esta

condición de contorno indica que no hay transferencia de flujo a través de dichos contornos.

La Figura 39 Muestra las Condiciones de contorno implementadas al modelo numérico.

Figura 39. Condiciones de Contorno

Fuente: Propia

98
5.1.4 Condiciones iniciales y objetivos de calibración

Los puntos de control (piezómetros y pozos) son la base principal para la validación del

funcionamiento del modelo numérico. Con el fin de contar con información del nivel

piezométrico en el área de estudio para la calibración del modelo, se ha utilizado los registros

de niveles de agua de los 14 puntos de monitoreo descritos en el apartado 3.5 (ver Figura 40).

La Tabla 34 muestra las cargas hidráulicas utilizadas para la calibración.

Tabla 34. Objetivos de Calibración

Carga
ID Este (m) Norte (m) Hidráulica
(msnm)
P-405 224798.00 8187247.00 2355.00
P-416 225271.00 8185200.00 2310.00
P-420 224290.00 8185884.00 2325.00
P-421 225033.00 8184624.00 2299.00
P-427 224152.00 8185015.00 2305.76
P-41 227137.00 8184683.00 2284.36
P-42 227288.00 8184425.00 2265.10
P-44 226988.00 8184426.00 2275.00
P-45 224640.00 8183980.00 2275.00
P-50 223171.00 8183732.00 2273.00
PY22 228168.00 8185975.00 2312.00
PC1 225478.44 8184597.44 2298.00
PC2 225717.64 8184611.90 2295.00
PC3 226048.24 8184476.29 2292.00

Fuente: Adaptado del ANA (2017)

99
Figura 40. Objetivos de Calibración

Fuente: Propia

5.1.5 Calibración del modelo para la condición inicial

La etapa de calibración del modelo puede realizarse de dos métodos, (1) calibración manual

y (2) calibración automática. En el presente trabajo se ha optado por el método de calibración

manual, realizando ajustes de valores de parámetros hidráulicos (conductividades hidráulicas)

para conseguir los niveles piezométricos medidos o definidos para este fin.

La calibración manual de los parámetros hidráulicos se ha realizado sobre la base de las

cargas hidráulicas observadas. El proceso de calibración consistió en minimizar los valores

estadísticos mostrados en la Tabla 35 a partir de la variación de los parámetros hidráulicos

(principalmente permeabilidad), entendiendo que valores estadísticos más bajos representan un

100
mejor ajuste entre los valores observados de carga hidráulica (niveles de agua medidos en

piezómetros) y los calculados por el modelo.

La Figura 41 presenta los niveles de agua simulados en condición estacionaria en

comparación con los niveles de agua medidos en los puntos de monitoreo.

Todos los estadísticos, están basados en el residual, o diferencia entre los valores observados

(Xobs) y los calculados por el modelo (Xcal) el cual se muestra en la ecuación (14).

𝑅𝑒𝑠𝑖𝑑𝑢𝑎𝑙 = 𝑋𝑐𝑎𝑙 − 𝑋𝑜𝑏𝑠 ( 54 )

La desviación estándar (s) es calculada a partir de la siguiente ecuación:

𝑁
( 55 )
1
𝑠=√ ∑(𝑋𝑜𝑏𝑠,𝑖 − 𝑋𝑐𝑎𝑙,𝑖 )2
𝑁−1
𝑖=1

Asimismo, el Error Residual Medio Cuadrático (RMSE, por sus siglas en inglés) es

calculado con la ecuación (16):

( 56 )
∑𝑁 (𝑋𝑜𝑏𝑠 − 𝑋𝑐𝑎𝑙 )2
𝑅𝑀𝑆𝐸 = √ 𝑖=1
𝑁

Tabla 35. Resultados estadísticos de la calibración del modelo

Estadístico Unidad Valor

Residual Medio (m) m -1.26


Promedio del Valor Absoluto de Residuales (m) m 3.52
Desviación Estándar de los Residuales (m) m 4.34
Suma de Residuales al Cuadrado, RSS m2 267.46
Error medio Cuadrático, RMS (m) m 4.37
Residual Mínimo (m) m -10.04
Residual Máximo (m) m 6.31
Numero de Observaciones - 14
Rango de Observaciones (m) m 89.90
Desviación Estándar Escalado Porcentaje 4.8%
RMS escalado Porcentaje 5%

Fuente: Propia

101
Figura 41. Cargas hidráulicas calculadas vs observadas.

Fuente: Propia

Tal y como se indica en Anderson y Woessner (2015): El criterio de cuando el modelo se

encuentra lo suficientemente ajustado a la realidad es subjetivo. Una regla común de aceptación

es que, en un modelo bien calibrado, el RMS debería ser menor que el 10% del rango de

mediciones en todo el dominio (RMS escalado). El residual promedio obtenido de la

calibración del modelo fue ligeramente negativo (-1.26 m) indicando una ligera tendencia

general de que los niveles de agua se encuentran por debajo de los observados. Los resultados

de la calibración indican un modelo con ajuste bastante aceptable a esta escala regional con un

promedio del valor absoluto de los residuales (ARM) de 3.52 m y un RMS escalado de 5 %.

102
5.2 Selección del sistema de desaguado

Las características hidrogeológicas de la zona del proyecto y la altura a drenar condicionan

preliminarmente la selección del sistema de desaguado. Asimismo se debe considerar el

proceso constructivo, presupuesto y cronograma de la obra.

Las permeabilidades del área de estudio varían entre 0.1 m/d ( 10−6 m/s) a 14.6 m/d

(1.6𝑥10−4 m/s) y 9 m de rebajamiento del nivel freático en la zona más crítica.

Figura 42. Rangos de aplicación de sistemas de desaguado

Fuente: Cassiraga y Ferrer (2008)

En el Capítulo II se realizó la descripción de los diferentes sistemas de desaguado. En la

Tabla 36 adicionalmente se muestra las limitaciones y ventajas de estos sistemas.

De acuerdo a las recomendaciones vistas según el tipo de suelo, la profundidad de

abatimiento del nivel y tipo de proyecto es que se eligió el sistema de bombeo mediante pozos

profundos.

103
Tabla 36. Ventajas y limitaciones de los sistemas de desaguado

Bombeo con pozos


Descripción Zanjas de drenaje Sistema de Wellpoints
profundos
Mejor rendimiento en suelos de
Arenas bien graduadas, Suelos de grano grueso
Suelos matriz arenosa sin limos, no
gravas limpias (grava y arenas).
plástico
- Los descensos conseguidos
- Rápida instalación (hasta 1000 están sólo limitados por la
metros por día en buenas potencia y profundidad
condiciones) con de la bomba y por la
- Suponen un menor coste máquinas especiales para estratificación del terreno.
presentando un resultado excavar zanjas. - Los pozos pueden ser
Ventajas satisfactorio como sistema - No se requiere circulación de emplazados fuera del área de
cuando las condiciones son agua para su construcción. trabajo.
apropiadas - Una vez que la tubería ha sido - Los pozos son construidos
instalada, la puesta en marcha y a gran distancia unos de los
desmantelamiento del sistema otros lo cual reduce
es simple y rápido. los problemas de acceso a la
obra.
- La altura de succión está
- Su menor impacto limitada a los 5 o 6 metros. Si
- Los costos de instalación
económico hace que en los suelos son blandos, el peso
por cada pozo son altos por
ocasiones, este método sea de las máquinas excavadoras
lo que requiere optimizar
utilizado en condiciones puede
su número.
inapropiadas, con resultados hacer inviable su utilización.
- Se requiere información
negativos, demoras en la - Las máquinas excavadoras
cualitativa acerca del sitio
Limitaciones ejecución, accidentes y pueden sufrir graves daños si el
donde se van a instalar
ocasionalmente grandes suelo a excavar
los pozos.
daños personales y contiene adoquines y rocas.
- Las bombas utilizadas
materiales. - La movilización de las
son eléctricas por lo que se
-Requiere de una mayor máquinas excavadoras es caro
requiere del correspondiente
experiencia, sentido común pudiendo condicionar
suministro de energía.
del experto. la instalación del sistema en
caso de obras pequeñas.

Fuente. Adaptado de Cassiraga y Ferrer (2008)

5.3 Diseño del sistema de desaguado

El objetivo de las simulaciones realizadas es presentar un plan de desaguado que permita

abatir el nivel del agua a lo largo de la traza de la obra, a fin de asegurar la seguridad y

estabilidad de los trabajos de excavación.

Tras la calibración del modelo numérico en condiciones estacionarias, se implementó

diferentes configuraciones de pozos de bombeo para evaluar los efectos que producía en el área

del Proyecto, y comprobar que efectivamente el nivel freático quedará por debajo de la cota de

104
máxima excavación y 1 metro adicional como factor de seguridad. La Figura 43 muestra un

esquema con las necesidades de rebajamiento por sectores.

Figura 43. Esquema de necesidades de rebajamiento.

Fuente: Propia

Acorde a lo descrito en el párrafo anterior, para la medida de los descensos del nivel freático

en el eje de la traza se han probado diferentes configuraciones de dispositivos de desaguado,

es decir, diferentes ubicaciones de los pozos con el fin de encontrar el más óptimo desde el

punto de vista del cumplimiento del rebajamiento y económicamente viable. A continuación,

se describe los resultados de los diferentes escenarios simulados.

5.3.1 Configuración de pozos laterales

Este arreglo consistió en ubicar un total de 38 pozos a ambos lados el Carril Central, es

decir, 19 pozos de bombeo por cada lado y distanciados a 40 - 60 m (Ver Figuras 44 y 45). Los

pozos de bombeo propuestos en esta configuración tendrían las siguientes características y

presentarían los siguientes resultados:

105
 38 pozos de bombeo, construidos de 20 a 25 m de profundidad.

 Diámetro de entubación recomendable de 8 pulgadas.

 Capacidad de bombeo de 10 a 5 l/s.

 Al final de la etapa de bombeo el caudal total se atenúa a 180 l/s.

 Los objetivos de rebajamiento del nivel freático se cumplen para todas las zonas.

5.3.2 Configuración de pozos centrales

Este arreglo consistió en ubicar un total de 21 pozos alineados al eje central del Carril

Central (Ver Figuras 46 y 47). Los pozos de bombeo propuestos en esta configuración tendrían

las siguientes características y presentarían los siguientes resultados:

 21 pozos de bombeo, construidos de 24 a 28 m de profundidad.

 Diámetro recomendable de 8 pulgadas.

 Capacidad de bombeo de 12 a 8 l/s.

 Al final de la etapa de bombeo el caudal total se atenúa a 210 l/s.

 Los objetivos de rebajamiento del nivel freático se cumplen para todas las zonas.

5.3.3 Configuración de pozos en tresbolillo

Este arreglo consistió en ubicar un total de 30 pozos, de los cuales 16 ubicados al norte del

carril y 14 al hacia el sur (Ver Figuras 48 y 49). Los pozos de bombeo propuestos en esta

configuración tendrían las siguientes características y presentarían los siguientes resultados:

 30 pozos de bombeo, construidos de 20 a 25 m de profundidad.

 Diámetro de entubación recomendable de 8 pulgadas.

 Capacidad de bombeo de 10 a 5 l/s.

 Al final de la etapa de bombeo el caudal total se atenúa a 175 l/s.

 Los objetivos de rebajamiento del nivel freático se cumplen para todas las zonas.

106
Figura 44. Simulación con distribución de pozos laterales – Parte 1.
Fuente: Propia

107
Figura 45. Simulación con distribución de pozos laterales – Parte 2.
Fuente: Propia

108
Figura 46. Simulación con distribución de pozos centrales – Parte 1.
Fuente: Propia

109
Figura 47. Simulación con distribución de pozos centrales – Parte 2.
Fuente: Propia

110
Figura 48. Simulación con distribución de pozos en tresbolillo – Parte 1.
Fuente: Propia

111
Figura 49. Simulación con distribución de pozos en tresbolillo – Parte 2.
Fuente: Propia

112
Capítulo VI

Análisis y resultados

6.1 Modelo conceptual hidrogeológico

Para realizar el modelo conceptual hidrogeológico se han seguido las directrices

mencionados en el capítulo III.

A partir de la investigación de campo que incluyen dos perforaciones diamantinas realizadas

por Perfocon y los estudios geofísicos que incluyen 26 líneas de Tomografía eléctrica, 4 de

Refracción Sísmica y 4 ensayos de MASW realizados por Rock&Soil S.A.C. se obtuvo la

geología local en la que se identificó la presencia de dos unidades litológicas los Depósitos

Aluviales (Qh-al) y el Volcánico Sencca (Tp-vsc), debido a que no se contó con ensayos fuera

del área de la zona de trabajos de la carretera variante de Uchumayo-Tramo II que permitan

definir a mayor precisión la geología del área de influencia del modelo numérico se

complementó con la geología regional.

Los estudios mencionados también permitieron conocer el nivel freático a precisión en la

zona de la vía, adicionalmente para abarcar toda el área de influencia del modelo se utilizaron

los datos piezométricos de pozos inventariados por la ANA en el año 2017, actualmente la

ANA se encuentra realizando la actualización de esta información pero aún la información no

es de acceso público. Con estos datos se logró realizar el mapa piezométrico del área de

influencia siguiendo lo mencionado en el ítem 2.4.1.

La caracterización de los depósitos encontrados en la zona permitió definir las Unidades

Hidrogeológicas Detrítica, Volcánica e Intrusiva y obtener sus parámetros hidrogeológicos

estos fueron obtenidos a partir de los estudios de campo, para el presente estudio fue necesario

tomar datos de otros estudios para poder tener un rango de valores asociados a cada unidad

como se menciona en el ítem 3.12.

113
A partir de lo mencionado se ha logrado definir el modelo hidrogeológico conceptual en el

que se describe a un acuífero libre, con cotas piezométricas entre 2190 a 2360 msnm, cuyo

flujo va de NW a SW y que debido al periodo corto (carácter temporal) de implementación del

sistema de desaguado subterráneo y considerando en este periodo las precipitaciones promedio

van de 0.2 a 0.6 mm/mes según la estación la Pampilla es que no se considerarán recargas.

6.2 Modelo numérico hidrogeológico

Se elaboró el modelo numérico en FEFLOW a partir de los del modelo conceptual.

Las simulaciones se realizaron en estado estacionario y transitorio con escenarios

hidrogeológicos de corto plazo. FEFLOW nos permite observar mediante gráficos la

información de los niveles de agua subterránea o líneas piezométricas en los que podremos

verificar la dirección del flujo y observar los niveles de rebajamiento con diferentes arreglos

de pozos y en diferentes periodos de tiempo.

6.2.1 Modelo hidrogeológico en estado estacionario

Utilizando la calibración manual se ha realizado una serie de iteraciones en los que se varia

la conductividad hidráulica para ajustar los objetivos de calibración (cotas piezométricas)

observando como variaba el RMS hasta lograr el valor recomendado, se logró calibrar el

modelo con un RMS escalado (5%) menor que el 10% del rango de mediciones en todo el

dominio que nos da un indicativo de un nivel de ajuste aceptable.

En la Tabla 34 se mostraron los objetivos de calibración (cotas piezométricas) que se

utilizaron. También se obtuvieron las permeabilidades de las tres Unidades Hidrogeológicas

las cuales se encontraron dentro del rango establecido en el modelo conceptual.

114
Tabla 37. Rango de valores de conductividad hidráulica asignados.

k (m/d)
Unidad Hidrogeológica
Rango de valores iniciales Valor calibrado

UH Detrítica 0.1 a 7.01 0.25

UH Volcánica 0.5 a 14.6 0.85

UH Intrusiva 10-5 a 0.01 0.08

Fuente: Propia

Los valores de conductividad hidráulica de esta UH obtenidos a partir de los ensayos de

permeabilidad en la zona van de 0.02 m/d a 0.63 m/d, el valor calibrado se ubica dentro de este

rango guardando correlación con las propiedades del material aluvial ubicado en la zona de la

vía.

Los valores de conductividad de la UH Volcánica tienen un rango de valores en la zona de

la vía entre 0.51 m/d y 7.35 m/d, el valor calibrado se encuentra dentro de este rango.

La UH Intrusiva no tiene representatividad en el área del modelo, su conductividad

hidráulica es baja en comparación a las otras UH y cercano al valor máximo del rango

establecido para esta, aceptando el valor calibrado de 0.08.

A partir del modelo hidrogeológico calibrado se pudo interpretar hidrogeológicamente el

acuífero en condiciones estáticas y nos permite obtener el mapa de cargas hidráulicas (Figura

50) donde se muestra correlación en el sentido del flujo y las cotas piezométricas del modelo

conceptual.

115
Figura 50. Mapa de isopiezas en condiciones estacionarias

Fuente: Propia

6.2.2 Modelo hidrogeológico con pozos de desaguado

El modelamiento se realizó en estado transitorio. Se realizaron tres configuraciones de pozos

laterales, centrales y en tres bolillos (triangular). De los resultados obtenidos en los

modelamientos se tiene el siguiente cuadro comparativo. La Tabla 38 muestra las

características de las bombas a utilizar y el caudal final.

Tabla 38. Comparación según configuración de pozos

Descripción Laterales Centrales Tres bolillo


N° de pozos 38 21 30
Profundidad 20-25 m 24-28 m 20-25 m
Diámetro 8" 8" 8"
Capacidad 10-5 l/s 12-8 l/s 10-5 l/s
Caudal final 180 l/s 210 l/s 175 l/s
Niveles de rebajamiento Cumple Cumple Cumple

Fuente: Propia

116
Se observa que el número de pozos del arreglo del tipo “pozos centrales” es menor en

comparación de los otros arreglos y también, que el caudal bombeado del arreglo en “tres

bolillo”, con el que se cumple los niveles de rebajamiento, es el menor de todos.

En la Tabla 39 se muestra los niveles máximos de rebajamiento a lo largo de diferentes

periodos de tiempo notando que los tres arreglos cumplen los niveles de rebajamiento

establecidos en la Figura 43.

Tabla 39.Descensos máximos de rebajamiento

Configuración de Pozos Configuración de pozos Configuración de pozos


Periodos
Laterales Centrales Tres bolillo
30 días 10 7 9
60 días 11 8 10
90 días 11 9 10
120 días 11 11 10
150 días 11 11 10
180 días 11 12 10
210 días 11 12 11
240 días 11 13 11

Fuente: Propia

El arreglo de tipo “pozos centrales” tiene el menor número de pozos y alcanza niveles de

rebajamiento mayores a los requeridos, pero debido al proceso constructivo de la vía es que se

elige el arreglo de “tres bolillos”; además se observa que el arreglo de “tres bolillo” a los 30

días ya sobrepasa los 9m de rebajamiento, que es el mayor descenso de nivel freático requerido.

Según el cronograma de obra, las excavaciones no deberían prolongarse en más de 80 días;

tiempo en que el nivel de rebajamiento para el arreglo “tres bolillos” es de 10 m.

117
Capítulo VII

Conclusiones y Recomendaciones

En el presente apartado se describe las principales conclusiones y recomendaciones

como resultados del estudio realizado.

7.1 Conclusiones

 A partir de los estudios geológico, geotécnico e hidrogeológico se obtuvo la

caracterización de la zona de estudio, y se concluye:

o Según la caracterización geológica, la zona del proyecto está ubicada en la

formación Volcánica Sencca (Tp-vs); formada por material tobáceo

limoarenoso (SM-SP); arenas de textura gruesa, media y fina; medianamente

denso a denso y sin plasticidad. Sobre esta encontramos sectores con material

del Aluvial Reciente (Q-al) de espesor variable; formada por arena, limos y

gravas heterogéneas (SM-SP); con intercalaciones de gravillas de piedra pómez

alterada y de consistencia medianamente denso a denso.

o Según la caracterización geotécnica en la zona del puente Atahualpa, el índice

de resistencia a la penetración obtenido del ensayo SPT fue de 10 a los 14 m de

profundidad y de 30 a los 17 m; pasando de consistencia medianamente

compacto a compacto. Además, de acuerdo a los ensayos de MASW los valores

de Vs están entre 185 m/s a 595 m/s. Por lo que, según la norma E030, se puede

clasificar al suelo como suelo intermedio (S2); es decir, un suelo medianamente

rígido; típico de arenas densas gruesas a densas medias. Y se concluye que el

suelo del proyecto tiene buenas propiedades mecánicas y que mejoran conforme

aumenta su profundidad.

o Según la caracterización hidrogeológica, se ha definido tres (03) unidades

hidrogeológicas (UH): La primera, la UH Detrítica formado por los depósitos

118
del Aluvial Reciente (Q-al); constituida por arenas y limos con presencia de

gravas. La segunda, la UH Volcánica formado por los materiales del Volcánico

Sencca (Tp-vse), constituido por el tufo limoarenoso. Y la tercera, la UH

Intrusiva representada por las rocas intrusivas que afloran en el sector sureste

del área del modelo hidrogeológico. La UH Detrítica y Volcánica forman parte

del acuífero libre, mientras que la UH Intrusiva forma parte del acuitardo.

 En el modelo conceptual se identificó que la permeabilidad de la UH Detrítica varía

entre 0.1 m/d a 7 m/d; de la UH Volcánica, varía entre 0.5 m/d a 14.6 m/d y de la UH

Intrusiva, varía entre 0.01 m/d a 1𝑥10−5 m/d. Asimismo, a partir de las mediciones de

nivel de agua realizadas en las perforaciones geotécnicas y la información disponible

de los niveles piezométricos de pozos; se pudo definir el mapa piezométrico y la

dirección del flujo subterráneo en la zona de estudio; donde se observa que las cargas

hidráulicas varían entre 2170 msnm a 2360 msnm y la dirección del flujo subterráneo

es principalmente de norte a sur.

 Se logró construir el modelo numérico de flujo subterráneo con el programa FEFLOW

6.2 a partir del modelo conceptual. Y tras la calibración del modelo numérico; se han

determinado permeabilidades de 0.25 m/d para la UH Detrítica, 0.85 para la UH

Volcánica y 0.08 para la UH Intrusiva; valores que se encuentra dentro del rango de

valores típicos para estos tipos de materiales.

 Se logró calibrar el modelo numérico hidrogeológico comparando los valores de nivel

freático medidos en campo; con los valores calculados por el modelo numérico en

estado estacionario, el valor obtenido con el método del Error de la raíz promedio

cuadrática (RMS) fue de 5% menor al 10% recomendado.

 En el modelo hidrogeológico en estado transitorio, como parte del diseño de desaguado;

se plantearon tres (03) escenarios de modelamiento para evaluar; el caudal de

119
extracción, la distribución de pozos de bombeo, el diseño de los pozos y el rebajamiento

del nivel freático. Los resultados obtenidos se describen a continuación:

o Esquema de desaguado con pozos laterales

- 38 pozos de bombeo de 20 a 25 m de profundidad y de 8 pulgadas de

diámetro, que permita la instalación de una bomba sumergible de 6

pulgadas de diámetro y capacidad bombeo entre 5 a 10 l/s.

- El caudal total bombeado es de 180 l/s.

- Se alcanzará los objetivos de rebajamiento tanto en el carril central, así

como en los tramos del puente Atahualpa y el intercambio a desnivel

Taboada.

o Esquema de desaguado con pozos centrales

- 21 pozos de bombeo de 24 a 28 m de profundidad y de 8 pulgadas de

diámetro, que permita la instalación de una bomba sumergible de 6

pulgadas de diámetro y capacidad bombeo entre 8 a 12 l/s.

- El caudal total bombeado es de 210 l/s.

- Se alcanzará los objetivos de rebajamiento tanto en el carril central, así

como en los tramos del puente Atahualpa y el intercambio a desnivel

Taboada.

o Esquema de desaguado con pozos en tresbolillo

- 30 pozos de bombeo de 20 a 25 m de profundidad y de 8 pulgadas de

diámetro, que permita la instalación de una bomba sumergible de 6

pulgadas de diámetro y capacidad bombeo entre 5 a 10 l/s.

- El caudal total bombeado es de 175 l/s.

120
- Se alcanzará los objetivos de rebajamiento tanto en el carril central, así

como en los tramos del puente Atahualpa y el intercambio a desnivel

Taboada.

 Dada la particularidad de las propiedades hidrogeológicas, la ingeniería del proyecto y

las solicitaciones en el descenso del nivel freático, como parte del presente estudio se

sugiere optar por el esquema de desaguado con pozos distribuidos en tresbolillo; en el

que se distribuyen la cantidad de 30 pozos y con el que se alcanzan los objetivos de

rebajamiento del nivel freático a los 30 días de implementado el sistema; resultando ser

la configuración más óptima.

 A partir del modelamiento numérico del flujo subterráneo en el programa FEFLOW 6.2

se logró evaluar los niveles de descenso del nivel freático y simular los escenarios de

extracción de corto y largo plazo, lo que permite fundamentar una correcta toma de

decisiones y lograr el diseño del sistema de desaguado de manera óptima.

7.2 Recomendaciones

 La construcción de los pozos debe realizarse con perforación tipo circulación reversa,

con recuperación de detritos del material perforado.

 Se recomienda la realización de pruebas airlift durante la perforación de los pozos, a

tramos de cada 6 m una vez interceptado el nivel freático.

 El desarrollo de los pozos mediante la inyección de aire comprimido o pistoneo debe

considerar un periodo de 6 a 8 horas como mínimo.

 Se recomienda la realización de pruebas de bombeo en algunos pozos, una vez

construido el pozo, para validar la capacidad de producción de los pozos y la

eficiencia de su construcción.

 Para el monitoreo de los niveles piezométricos durante la construcción de la obra y

garantizar la excavación en seco se recomienda la construcción de piezómetros tipo

121
Casagrande de 2” de diámetro con profundidades similares al de los pozos. Debe

considerarse piezómetros tanto como fuera y dentro del carril central.

 Durante la operación del bombeo se debe colocar un flujómetro en cada brocal de pozo

para contabilizar la cantidad de agua extraída.

 Tomar muestras de agua para análisis de manera mensual para controlar la

hidroquímica del agua extraída.

 Realizar pruebas de bombeo en los pozos y pruebas slug en los piezómetros durante la

implementación de los pozos y antes de la puesta en operación. Las pruebas de bombeo

deben considerar la realización de pruebas de bombeo escalonado y a caudal constante,

con el fin de verificar la eficiencia de la construcción de los pozos y el radio de

influencia. Las pruebas slug permitirán obtener información adicional de permeabilidad

de los materiales con el fin de realizar el proceso de verificación del modelo numérico

calibrado.

 Con el fin de realizar el monitoreo de los niveles de agua antes, durante y al cese de la

operación del sistema de desaguado, se recomienda implementar la construcción de

piezómetros de monitoreo tipo Casagrande de 2” de diámetro, situados a una separación

de 50 m, alcanzando un total de 20 piezómetros. Las profundidades de los piezómetros

deben ser no mayor a de 15 m.

 Se recomienda validar el modelo 30 días después de ser implementado a partir de los

caudales y niveles de agua registrados en campo.

 Al finalizar el bombeo se recomienda clausurar los pozos con concreto, para evitar

posibles accidentes y realizar la entrega de los piezómetros a la entidad correspondiente.

 Se recomienda contar con suministros alternos de energía como grupos electrógenos.

122
Referencias bibliográficas

Aguilar Bardales, Z. (1991). Microzonificación sísmica de la ciudad de Arequipa. Universidad

Nacional de Ingeniería. Lima, Perú.

Alva Hurtado, J., Aguilar Bardales, Z. (2002). III Curso internacional microzonificación y su

aplicación en la mitigación de desastres. Recuperado de

www.jorgealvahurtado.com/files/redacis19_a.pdf.

Anderson, M., Woessner, W. (2015). Applied groundwater modelling. Elsevier. Londres.

Cashman, P., Preene, M. (2001). Groundwater lowering in contruction. Taylor & Francis.

Londres.

Autoridad Nacional del Agua (2018). Estudio hidrogeológico del acuífero del río Chili.

Recuperado de http://repositorio.ana.gob.pe/handle/20.500.12543/4361.

Comisión Docente Curso Internacional de Hidrología subterránea. (2009). Hidrogeología.

FCIHS. Barcelona, España.

Custodio, E., Llamas M. (1996). Hidrología subterránea, Tomo I y II. Omega. Barcelona,

España.

Department of the Environment and the North West Water Authority (2015). Control of

groundwater for temporary works. CIRIA. Londres.

DHI-WASY GmbH. (2018). Modelación de agua subterránea con FEFLOW 6.2.

Ferrer Granell, A. (2010). Control de las aguas subterráneas en la ingeniería civil. Interacción

entre la obra y el medio hidrogeológico, síntesis de métodos de control y aplicación de modelos

matemáticos. Recuperado de https://riunet.upv.es/handle/10251/14198.

123
Gobierno Regional de Arequipa. (2018). Informe de investigaciones geotécnicas mediante

perforación diamantina. Puente Atahualpa y Taboada. Arequipa. Perú.

INGEMMET. (2016). Geología y mapas de peligros del volcán Misti. Lima. Perú.

ISCHEBECK. (2019). Catálogo de sistema de agotamiento Wellpoint. Recuperado de

http://www.ischebeck.es/

Koasaka, R., Macedo, L. (2000). Mapa de peligro potencial del volcán Misti. Universidad

Nacional de San Agustín. Arequipa. Perú.

Marín Torres, Alberto Raúl. (2017). Estudio de aguas subterráneas del río Yura, su

aprovechamiento hídrico para su parque industrial y zonas urbanas. Universidad Nacional de

San Agustín. Arequipa, Perú.

NBR 6484. Ensayos simples de reconocimiento con SPT. Asociación Brasilera de Normas

Técnicas. Rio de Janeiro. Brasil. Febrero del 2001.

Norma Técnica E030. Diseño Sismoresistente del Reglamento Nacional de Edificaciones.

Ministerio de Vivienda Saneamiento y Construcción. Lima. Perú. Febrero del 2019.

Rojas Rubio, H. (2008). Modelamiento para el análisis de la variación hidrogeológica espacial

del acuífero la Yarada – Tacna. Universidad Nacional de Ingeniería. Lima, Perú.

Rock&Soil Ingeniería S.A.C. (2018). Expediente Técnico del drenaje subterráneo del proyecto:

“Mejoramiento de la carretera Variante Uchumayo, entre el puente San Isidro y la Vía de

Evitamiento distrito de Sachaca, Yanahuara y Cerro Colorado tramo II”. Arequipa. Perú.

Sánchez San Román, F. (2017). Hidrología superficial y subterránea. CreateSpace Independent

Publishing Platform. Salamanca, España.

124
Tupa Fernandez, N. (1989), La hidráulica subterránea de la ciudad de Arequipa. Universidad

Nacional de San Agustín. Arequipa, Perú.

Vargas V., L. (1970). Geología del Cuadrángulo de Arequipa 33-s. INGEMMET. Lima. Perú.

125

También podría gustarte