3.1.

METODOLOGÍA EMPLEADA La determinación del radio de influencia (R0) se realiza a partir de la fórmula de Sichardt, que estima el radio de
influencia de un pozo genérico (Rp):
Para abordar los cálculos de caudal se han utilizado varios métodos analíticos de cálculo que, en esencia, asumen
que la obra prevista se va a comportar como una captación de agua subterránea, resolviendo las ecuaciones de
R p =3000· d· √ K
flujo correspondientes.
… así como del Radio del pozo equivalente (r). En Jiménez Salas (1980) se asimila directamente como la anchura
Los dos primeros corresponden a la Fórmula del Pozo Equivalente y la Fórmula de la Zanja Equivalente. Los de la excavación. Sin embargo, aquí se ha optado por definirla asumiendo que el área de influencia de dicho pozo
desarrollos matemáticos y consideraciones correspondientes vienen recogidas en tres publicaciones: “Geotecnia es igual al área del agotamiento (Sa):
y Cimientos”. Vol. III, parte 2 (Jiménez Salas, 1980), “Manual de Hidrología Subterránea”. Vol. I y II (Custodio, E. y
Llamas, R., 1983), “Mecánica de Suelos” (Berry, P. L. y Reid, D., 1993) e “Hidrogeología aplicada a la obra civil”
r =√ Sa/ π
(Garrido Ruíz, M. S., 2012).
La relación entre estos dos parámetros, finalmente, define el Radio de influencia del pozo equivalente según la
siguiente expresión:
El tercer método, recogido en

R0 =√ R p +r
2 2
Asimismo, y en contraste con los métodos analíticos, se realiza un modelo matemático de las secciones XXXX…

Ambas expresiones asumen que la excavación penetra totalmente en el acuífero, lo cual, por regla general, no se En el caso de tratarse de excavaciones de solares de grandes dimensiones, donde el radio equivalente (r) puede
ajusta a la realidad; por ello se suele considerar cierta reducción de los valores obtenidos, si bien resulta superar la tercera parte del radio de influencia del pozo (Rp), resulta preferible la utilización de la expresión
aconsejable no rebajar el caudal por estos motivos más de un 30%. siguiente, derivada de la fórmula original:

3.1.1.

[Link].
MÉTODOS ANALÍTICOS

Fórmula del Pozo Equivalente

Q=π·K· ( 2 H−d ) · d ·
( 2r
RP
+ 0,25
)
En este sentido, la Fórmula denominada del Pozo Equivalente (según Custodio y Llamas, 1982 y Jimenez Salas,
1980) considera que el caudal total susceptible de filtrarse al interior de la excavación es similar al que captaría [Link]. Fórmula de la Zanja Equivalente (según Custodio y Llamas, 1982)
un pozo circular cuya superficie fuera igual a la del solar. Esta fórmula es de aplicación en excavaciones de
geometría aproximadamente cuadrada, cuando la anchura del solar es superior a la cuarta parte de su largo. El El segundo método de cálculo utilizado corresponde a la aplicación de la Fórmula denominada de la Zanja
caudal viene determinado a partir de la Fórmula de Dupuit, para acuíferos libres, introduciendo en ella el radio equivalente (según Custodio y Llamas, 1982), de aplicación cuando el agotamiento tiene una geometría alargada
del pozo equivalente: (la relación de aspecto entre el largo y ancho del solar supera el valor 4: L/a > 4). El cálculo asume cada lateral de
la excavación como una zanja en un acuífero libre, alimentada por sólo por un lateral:
( 2 H−d ) · d d
Q=π·K· ≈ 2· π·H·K·
K
( )
R
ln 0
r ( )
R
ln 0
r
Q=L·
2 RZ
· ( 2 H−d ) · d

Donde: Donde:

Q es el Caudal captado, en m3/s. Q es el Caudal captado, en m3/s.

K es el Coeficiente de permeabilidad, en m/s. K es el Coeficiente de permeabilidad, en m/s.
H es el Espesor estático de acuífero saturado, en m. H es el Espesor estático de acuífero saturado, en m.
d es el Descenso producido, en m. d es el Descenso producido, en m.
R0 es el Radio de influencia del pozo equivalente, en m. RZ es la Distancia de influencia de la captación, en m.
r es el Radio del pozo equivalente, en m. L es la Longitud de la zanja, en m.

La aplicación de esta fórmula ha de realizarse, lógicamente, a cada uno de los laterales de la excavación, puesto
que el método de la Zanja equivalente, asume que el caudal proviene sólo de un lateral de la zanja. En el caso
concreto que nos ocupa, al tratarse de un recinto rectangular, se tienen dos zanjas de 120 m y otras dos de 20 m. El perfil a modelizar consta de dos capas horizontales de terreno natural que, de techo a muro, son:
Por ello, el cálculo para un solo lateral se multiplica, en ambos casos, x 2.
- CAPA 1: Capa correspondiente al Tosco arenoso.
La estimación de la Distancia de influencia (RZ), se aplica nuevamente la fórmula de Sichardt, si bien, en esta - CAPA 2: Capa correspondiente a la Peñuela.
ocasión, cambia el coeficiente a utilizar: - CAPA 3: Capa correspondiente a la pantalla de pilotes en tramo de Tosco.
- CAPA 4: Capa correspondiente a la pantalla de pilotes en tramo de Peñuela.
R Z =15 00 ·d· √ K
En la siguiente figura se muestra la malla final y los materiales asignados a cada elemento. Se ha exagerado la
Donde: escala vertical al doble de la horizontal, con objeto de que los detalles en la vertical de las sucesivas simulaciones
sean apreciables. En el apéndice a la presente memoria se incluye la colección completa de perfiles obtenidos en
K es el Coeficiente de permeabilidad, en m/s. los distintos pasos de la simulación.
d es el Descenso producido, en m.
RZ es la Distancia de influencia de la captación, en m. CAPA 3
CAPA 1

[Link]. Fórmula de la Zanja Equivalente

3.1.2. MODELO MATEMÁTICO DE FLUJO EN 2D CAPA 4

CAPA 2
Adicionalmente, y a modo de comprobación, se ha procedido a modelizar una sección en dos dimensiones
mediante el módulo SEEP/W del paquete de software geotécnico GEO-SLOPE 2012, que permite la construcción
de secciones mediante elementos finitos y la modelización de flujos de agua en ellas, una vez definidos los
parámetros hidrogeológicos (permeabilidad, coeficiente de almacenamiento, etc.) y condiciones de contorno
para los distintos materiales. Para ello, se ha construido un perfil tipo transversal a la estructura E-3, a partir de
las testificaciones de los sondeos más próximos a la estructura (SO-7 y SO-8). Se trata de un perfil idealizado que
busca representar las características del citado tramo y que permite estimar los caudales de filtración a través del
suelo de la excavación.
Figura 1. Malla del perfil transversal objeto de modelización (elaboración propia).
Se ha procedido a modelizar una sección tipo del recinto, considerando los siguientes parámetros:

Ancho Cota
Cota Cota Prof. Cota Asimismo, se ha considerado el conjunto de las pantallas de pilotes como un tercer elemento del modelo. Para la
Sondeo entre fondo Cota NF
Sección sondeo terreno empotramiento contacto
Ref. pantallas excavac. (msnm) estimación de sus correspondientes parámetros geotécnicos se ha tenido en cuenta que el diámetro de los
(msnm) (msnm) pantallas (m) (msnm)
(m) (msnm)
pilotes es de aprox. 0,80 m, que están separados aprox. 1,20 m (con respecto a sus respectivos ejes) y que, por
SO-07 613,3 tanto, éstos no son secantes. Ello implica una sección hormigonada del 66,67%, dejando una superficie libre de
Sección 1 12,4 612,76 605,0 601,0 607,0 610,0
SO-08 612,7
terreno natural en la pantalla de pilotes del 33,33%. Esta estimación de la permeabilidad de las pantallas se
Tabla 1. Consideraciones geométricas de la sección modelizada (elaboración propia). parametriza en términos hidrogeológicos de la siguiente forma.

Cabe recordar que, de forma general, el programa asume que la sección modelizada presenta fondo de 1 m, por Parámetros hidrogeológicos
lo que los valores de caudal que se obtengan como resultado de cada simulación, han de interpretarse como
caudales unitarios por sección. Para la obtención del caudal total correspondiente al recinto completo, habrá que Según las estimaciones efectuadas en el apdo. [Link].1., se han tomado los siguientes valores iniciales de
multiplicar cada dato directo por la longitud de la estructura según su eje. referencia para definir cada uno de los materiales del modelo.

Construcción de la malla
 POTENCIA PERMEABILIDAD RELACIÓN COEF.
CAPA MATERIAL
(m) HORIZ. KH (m/s) KV/KH ALMACEN.
Arenas finas –
CAPA 1 2,76. 1,6 x 10-05 0,5 0,32 (32%)
Terreno medias, arcillosas
natural
CAPA 2 > 30 m Arcillas 1,03 x 10-07 0,5 0,02 (2%)

CAPA 3 Hormigonado Tosco:

(permeabilización 5,33 x 10-6 (Tosco) 0,5
(PILOTES) Material 0,11 (10,66%)
1 m. del 33,3% de la
CAPA 4 artificial Peñuela:
sección según el 3,43 x 10-8 (Peñuela) 0,5
(PILOTES) terreno). 0,0067 0,67%)
Tabla 2. Valores iniciales de referencia para los parámetros hidrogeológicos del modelo

Figura 2. Curvas de variación del contenido volumétrico de agua en saturación, en función de la presión de
La permeabilidad asignada al Hormigón, 1,0 x 10e-12, se ha tomado de la ponencia “Hormigón: Propiedades del
poro para los distintos materiales del modelo.
hormigón endurecido” (UNIVERSIDADE DA CORUÑA E.T.S.I. Caminos, Canales y Puertos).
Por su parte, el comportamiento de la permeabilidad en saturación, en función de la presión de poros, se ha
El coeficiente de almacenamiento para cada uno de los materiales considerados en el modelo se parametriza a estimado a partir del contenido volumétrico de agua y del valor de referencia de la permeabilidad horizontal (KH)
través del Contenido Volumétrico de Agua (Volumetric Water Content), que es una función de la presión de en saturación (Pression = 0), para cada tipo de material considerado. En la figura que se muestra a continuación
poros. Valores positivos de la presión de poros corresponden a medios saturados, mientras que valores negativos se presentan las curvas de ajuste para cada tipo de material del terreno.
corresponden a medios insaturados en agua. En la siguiente figura se presentan las curvas de ajuste para cada
Asimismo, para cada material se ha asumido una permeabilidad vertical (KV) un orden de magnitud inferior a la
tipo de material del terreno.
horizontal, de forma que la relación KH / KV = 0,5.
 Figura 3. Curvas de variación de la permeabilidad saturada en función de la presión de poro para los
distintos materiales del modelo

Se ha considerado un nivel freático constante a una cota de 607,0 m.s.n.m. (coherentemente con las
consideraciones establecidas para las estimaciones mediante la aplicación de las Fórmulas de la Zanja
Equivalente y del Pozo Equivalente).

Asimismo, a partir de las testificaciones de los sondeos más próximos se ha establecido una cota de contacto
entre el tramo multicapa de las granulometrías más gruesas Tosco arenoso (TTA), Arena tosquiza (TAT) y Tosco (TT)
y el tramo inferior dominado por las facies Peñuela en torno a la cota 610,0 m.s.n.m.

Por último, se ha estimado un peso específico del agua de 9,807 kN/m3 y un coeficiente de compresibilidad del
agua constante (Mv) igual a 1 x 10-05 kPa.

Condiciones de contorno

En los extremos de la primera capa del modelo se ha establecido la condición de potencial hidráulico constante
(Head - H) igual a la cota 607,0 m.s.n.m., correspondiente al nivel freático. Ello implica que el flujo se establece
exclusivamente hacia el interior del recinto objeto de modelización.

Asimismo, se ha establecido la condición de flujo nulo (Total Flux = 0) a través de las paredes de la pantalla de
pilotes, lo que condiciona que la entrada de agua se produzca exclusivamente a través del fondo de la
excavación.

Por último, a efectos de asegurar el flujo hacia el interior de la excavación, rodeando el pie de las pantallas, se ha
establecido la condición de Presión de poros (Press. Head) = 0.

Asimismo, se ha dispuesto de una sección de flujo (Flux section) a lo largo de toda la excavación, a fin de evaluar
el caudal filtrado.
3.2. CALCULO DE LOS CAUDALES DE INFILTRACIÓN

3.2.1. CONDICIONES DE CONTORNO

Para abordar las estimaciones de los caudales de infiltración se ha asumido una situación de máximos, a fin de
estar de lado de la seguridad. En este sentido, a continuación se especifica la parametrización realizada con
respecto a la piezometría de referencia y al dimensionamiento del solar.

- En un primer escenario, se asume que la totalidad que los caudales susceptibles de producirse en obra
van a provenir de la primera capa terciaria considerada (que constituye el acuífero multicapa), integrada
por un conjunto de lentejones interdigitados de Tosco arenoso (TTA), Arena tosquiza (TAT) y Tosco (TT). A
efectos de cálculo, se considera como un conjunto homogéneo al que se le asignará la permeabilidad
determinada en el apartado [Link].1.1. (1,60E-05 m/s).

- Se analiza un segundo escenario en el que la totalidad que los caudales susceptibles de producirse en
obra van a provenir de los materiales de la Peñuela (TP). A efectos de cálculo, se considera como un
conjunto homogéneo al que se le asignará la permeabilidad determinada en el apartado [Link].1.2.
(1,03E-07 m/s).

- Se ha tomado como referencia el nivel freático más elevado registrado en las mediciones del 24 de enero
de 2019, que corresponde al sondeo SO-7 (cota freática de 605,96 m.s.n.m.). A fin de simular un
escenario de aguas altas razonable, se ha optado por incrementar dicha cota 1 metro (redondeando a
607 m.s.n.m.), suponiendo así que ésta sería la altura máxima que alcanzaría el nivel freático con una
recarga máxima en aguas altas.

- Esto permite estimar un espesor saturado máximo de unos 6,85 m en la zona de la excavación.

- Por lo que se refiere a la cota de agotamiento, se ha tomado la cota mínimo de excavación en el punto
más bajo del trazado (604,774 m.s.n.m.), que se sitúa en el entorno del sondeo S-3. Con respecto a ésta
se ha asumido necesaria una depresión del nivel del agua en esta zona, aproximadamente 1 metro por
debajo (hasta la cota 603,5 m.s.n.m.).

Figura 4. Perfil geológico – geotécnico y Planta estructural del entorno de la excavación. En el primero, se
indica la longitud máxima considerada para el mismo.

- Esto arroja un abatimiento freático máximo de 3,5 metros aprox. La zona donde debería centrarse este
agotamiento corresponde al área entre los sondeos S-3 y SE-45, que es donde la cota de excavación
resulta mínima.

- Se considera la excavación en terreno natural, sin tener en cuenta las pantallas de pilotes.
 - Con respecto a las dimensiones de la excavación, se ha asumido una longitud de 120 m, lo cual (como se Ratio r/(R/3) (*) r/(R/3) 1.43 ---
puede apreciar en la figura anterior), no abarca la totalidad del tramo soterrado. Esto es así debido a que (*) En grandes excavaciones donde el radio equivalente (r) puede superar la 3ª parte del radio de influencia del pozo (Rp =
el contacto entre las facies T TA y TP se encuentra a una cota más elevada en la mitad occidental del mismo r/(Ro/3) > 1), el cálculo del caudal se realiza según la expresión para grandes solares.
(en el entorno del S-9) de forma que apenas existe zona saturada en esta zona, por lo que no es objeto
de agotamiento. En este sentido, el rebaje del nivel freático sólo afectaría a la zona saturada presente en Tabla 4. Radios de influencia del agotamiento, según la Fórmula del Pozo equivalente.
el área de los sondeos SE-45, S-2 y S-3.

Con respecto a la anchura del solar, como se aprecia en la figura anterior, ésta se cifra netamente en 10,4
m; no obstante, en los cálculos se ha asumido un valor de 20 m, de forma conservativa. Caudal estimado Abrev. Valor Unidades

Caudal captado (m3/s) Expresión general: Q 0.003 m3/s

Caudal captado (m3/día) ( 2 H−d ) · d Q 259.10 m3/día
Q=π·K·
( )
3.2.2. ESTIMACIÓN DE CAUDALES SEGÚN EL DRENAJE EN MATERIALES DE TOSCO ARENOSO (TTA) Caudal captado (L/s) R0 Q 3.00 L/s
ln
Caudal captado (L/día) r Q 2.59E+05 L/día
[Link]. Aplicación de Fórmulas analíticas Expresión aproximada:
Caudal captado (m3/s) Q 0.00 m3/s

Cálculo mediante la aplicación de la fórmula del Pozo Equivalente Caudal captado (m3/día) d Q 348.01 m3/día
Q ≈ 2· π·H·K·
R0
( )
Caudal captado (L/s) Q 4.03 L/s
En las tablas siguientes se resumen los parámetros utilizados y los resultados obtenidos: ln
Caudal captado (L/día) r Q 3.48E+05 L/día
Caudal captado (m3/s) Expresión para grandes solares: Q 0.003 m3/s
Parámetro Abrev. Valor Unidades
Caudal captado (m3/día) Q=π·K· ( 2 H−d ) · Q 242.79 m3/día
Espesor saturado estático del acuífero (m) Ho 6.85 m

( 2r
)
Caudal captado (L/s) Q 2.81 L/s
Nivel piezométrico a conseguir (m) H 3.35 m · d· +0,25
Caudal captado (L/día) RP Q 2.43E+05 L/día
Descenso producido (m) d 3.5 m
Tabla 5. Estimación de caudales a drenar en el recinto de la excavación, según la aplicación de la Fórmula
Coef. permeabilidad (m/s) K 1.60E-05 m/s
del Pozo equivalente.
Coef. permeabilidad (m/día) K 1.38E+00 m/día

Transmisividad (m2/s) T 1.10E-04 m2/s

Transmisividad (m2/día) T 9.47E+00 m2/día

Como puede apreciarse en la tabla anterior, según se aplique una aproximación u otra a la fórmula de Dupuit, se
Longitud del solar (m) L 120.00 m obtiene un caudal total a drenar de entre 2,81 y 4,03 L/s. Asimismo, se estima un radio de influencia de la
Anchura del solar (m) a 20.00 m captación equivalente en torno a 50 m. Teniendo en cuenta que el punto de agua inventariado más próximo a la
Área del solar (m2) - Área del Pozo equivalente Sa 2400.00 m2 estructura (1922-4-0114) se encuentra a una distancia de 38 m, según la perpendicular a la misma, a priori, éste
se vería algo afectado por el área de influencia de un posible drenaje de aguas subterráneas. Sin embargo, esto
Relación aspecto L/a (*) L/a 6.00 ---
no ocurriría en el caso de la cabecera del Ayo. de Rejas, situada a más de 100 m de la estructura.
(*) Cuando L/a > 4, ha de calcularse como una Zanja Drenante.

Tabla 3. Parámetros utilizados para el cálculo de caudales en el recinto de la excavación, según la También es preciso mencionar que, como puede observarse en la primera de las tablas, la relación de aspecto de
aplicación de la Fórmula del Pozo equivalente. la excavación aconseja utilizar la expresión correspondiente a la Zanja Equivalente.

Cálculo mediante la aplicación de la fórmula de la Zanja Equivalente

Radio de influencia Abrev. Valor Unidades

Radio Influencia captación genérica (Formula Sichart) Rp 42.00 m De la misma forma, en las siguientes tablas se resumen los parámetros utilizados y los resultados obtenidos de la
Radio equivalente del solar r 27.64 m aplicación de la Fórmula de la Zanja equivalente a cada uno de los laterales de 120 m. Se ofrece también el
Radio de influencia del Pozo Equivalente Ro 50.28 m
caudal por metro lineal de ambos laterales, en conjunto:
 Parámetro Abrev. Valor Unidades Parámetro Abrev. Valor Unidades
Espesor saturado estático del acuífero (m) Ho 6.85 m Transmisividad (m2/s) T 1.10E-04 m2/s
Nivel piezométrico a conseguir (m) H 3.35 m Transmisividad (m2/día) T 9.47E+00 m2/día
Descenso producido (m) d 3.5 m Longitud de la zanja (m) L 20.00 m
Coef. permeabilidad (m/s) K 1.60E-05 m/s Tabla 9. Parámetros utilizados para el cálculo de caudales a drenar en el recinto de la excavación, según
Coef. permeabilidad (m/día) K 1.38E+00 m/día la aplicación de la Fórmula de la Zanja equivalente (laterales de 20 m).
Transmisividad (m2/s) T 1.10E-04 m2/s

Transmisividad (m2/día) T 9.47E+00 m2/día

Longitud de la zanja (m) L 120.00 m Radio de influencia Abrev. Valor Unidades

Tabla 6. Parámetros utilizados para el cálculo de caudales a drenar en el recinto de la excavación, según Distancia de Influencia captación (Formula Sichart) RZ 21.00 m
la aplicación de la Fórmula de la Zanja equivalente (laterales de 120 m). Tabla 10. Radios de influencia del agotamiento, según la Fórmula de la Zanja equivalente.

Radio de influencia Abrev. Valor Unidades Valor Valor

Caudal estimado Abrev. Unidades
Distancia de Influencia captación (Formula Sichart) RZ 21.00 m 1 sólo lateral 2 laterales (x2)

Tabla 7. Radios de influencia del agotamiento, según la aplicación de la Fórmula de la Zanja equivalente. Caudal captado (m3/s) Q 0.0003 0.001 m3/s
Expresión general:
Caudal captado (L/s) Q 0.272 0.544 L/s
K
3
Caudal específico (m /s x m.l.) Q=L· · ( 2 H−d ) · d q 0.00001 0.00003 m3/s x m.l.
2 RZ
Valor Valor Caudal específico (L/s x m.l.) q 0.014 0.027 L/s x m.l.
Caudal estimado para L = 120 m Abrev. Unidades
1 sólo lateral 2 laterales (x2) Tabla 11. Estimación de caudales a drenar en el recinto de la excavación, según la aplicación de la Fórmula
3
Caudal captado (m /s) Q 0.002 0.003 m /s3 de la Zanja equivalente (laterales de 20 m).
Expresión general:
Caudal captado (L/s) Q 1.632 3.264 L/s
K
Caudal específico (m3/s x m.l.) Q=L· · ( 2 H−d ) · d q 0.00001 0.00003 m3/s x m.l.
2 RZ
Caudal específico (L/s x m.l.) q 0.014 0.027 L/s x m.l. Como puede comprobarse, entre los cuatro laterales se obtiene un caudal total de 3,80 L/s para todo el recinto
Tabla 8. Estimación de caudales a drenar en el recinto de la excavación, según la aplicación de la Fórmula de la excavación. Por su parte, se ha estimado un caudal específico para cada metro de zanja, de 0.014 L/s por m.
de la Zanja equivalente (laterales de 120 m). l. y un radio de influencia máximo perimetral de la excavación de aproximadamente 21 m. No obstante, éste
quedará acotado por la limitación en la extensión de los lentejones más permeables, por lo que posiblemente
sea inferior. En cualquier caso, el punto de agua inventariado más próximo a la estructura (1922-4-0114) se
encuentra a una distancia de 38 m, según la perpendicular a la misma, por lo que, a priori, quedaría fuera del
Por su parte, en las siguientes tablas se resumen los parámetros utilizados y los resultados obtenidos de la área de influencia de un posible drenaje de aguas subterráneas. Tampoco se vería afectada la cabecera del Ayo.
aplicación de la Fórmula de la Zanja equivalente a cada uno de los laterales de 20 m. Se ofrece también el caudal de Rejas, situada a más de 100 m de la estructura.
por metro lineal de ambos laterales, en conjunto:

Parámetro Abrev. Valor Unidades

[Link]. Cálculo mediante modelo matemático con SEEP/W
Espesor saturado estático del acuífero (m) Ho 6.85 m

Nivel piezométrico a conseguir (m) H 3.35 m De forma paralela a los anteriores cálculos analíticos, se ha procedido a realziar una estimación de caudales de
Descenso producido (m) d 3.5 m infiltración en terreno natural, con las características del Tosco arenoso (TTA), según la parametrización expuesta
en el apdo. [Link].3. En este primer escenario, y de forma coherente con la asunción de los cálculos anteriores,
Coef. permeabilidad (m/s) K 1.60E-05 m/s
se ha asumido que la totalidad que los caudales susceptibles de producirse en obra van a provenir de la primera
Coef. permeabilidad (m/día) K 1.38E+00 m/día
capa terciaria considerada (que constituye el acuífero multicapa), integrada por un conjunto de lentejones
interdigitados de Tosco arenoso (TTA), Arena tosquiza (TAT) y Tosco (TT).

De esta forma, queda una malla constituida por una única capa equivalente a la CAPA 1 del modelo genérico
construido:

POTENCIA PERMEABILIDAD RELACIÓN COEF.

CAPA MATERIAL
(m) HORIZ. KH (m/s) KV/KH ALMACEN.
Terreno Arenas finas –
CAPA 1 2,76. 1,6 x 10-05 0,5 0,32 (32%)
natural medias, arcillosas
Tabla 12. Valores iniciales de referencia para los parámetros hidrogeológicos del modelo

En las siguientes figuras se ilustran tanto la constitución de la malla inicial como el resultado de la modelización.
Éstas han sido incorporadas también al Apéndice 4.

Figura 6. Resultados de la modelización de flujo en el perfil transversal objeto de modelización en la

hipótesis simple de Tosco arenoso (elaboración propia).

Como puede apreciarse en la última de ellas, la sección de flujo arroja un valor de 2,9664E-05 m 3/s y por metro
lineal de sección transversal, que equivalen a 2,9664E-02 L/s (0,029664 L/s). Considerando la longitud total de la
estructura asumida en los cálculos analíticos (120 m), el caudal obtenido según la modelización resulta ser de
3,55 L/s. Este resultado es coherente con los rangos obtenidos mediante la aplicación de las fórmulas del Pozo
equivalente (entre 2,81 y 4,03 L/s) y la Zanja equivalente (en torno a 3,80 L/s). No obstante, puede considerarse
la estimación mediante modelización matemática de flujo más precisa puesto que incluye en el cálculo otros
parámetros hidrogeológicos (como el coeficiente de almacenamiento, la succión, etc.) que no se contemplan en
las fórmulas analíticas.

3.2.3. ESTIMACIÓN DE CAUDALES SEGÚN EL DRENAJE EN MATERIALES DE PEÑUELA (TP)

[Link]. Aplicación de Fórmulas analíticas

Figura 5. Malla del perfil transversal objeto de modelización en la hipótesis simple de Tosco arenoso
(elaboración propia).
Cálculo mediante la aplicación de la fórmula del Pozo Equivalente

En las tablas siguientes se resumen los parámetros utilizados y los resultados obtenidos:

Parámetro Abrev. Valor Unidades

Espesor saturado estático del acuífero (m) Ho 6.85 m

Nivel piezométrico a conseguir (m) H 3.35 m

Descenso producido (m) d 3.5 m

Coef. permeabilidad (m/s) K 1.03E-07 m/s

 Parámetro Abrev. Valor Unidades Tabla 3. Estimación de caudales a drenar en el recinto de la excavación, según la aplicación de la Fórmula
Coef. permeabilidad (m/día) K 8.90E-03 m/día
del Pozo equivalente.

Transmisividad (m2/s) T 7.06E-07 m2/s Como puede apreciarse en la tabla anterior, según se aplique una aproximación u otra a la fórmula de Dupuit, se
Transmisividad (m2/día) T 6.10E-02 m2/día obtiene un caudal total a drenar de entre 0,19 y 2,1 L/s. Por lo que respecta al radio de influencia de la captación
Longitud del solar (m) L 120.00 m equivalente, en este escenario de drenaje a través de las facies Peñuela, se cifra en 27,84 m, por lo que el punto
de agua 1922-4-0114, a una distancia de 38 m según la perpendicular a la misma, quedaría fuera del área de
Anchura del solar (m) a 20.00 m
influencia. Tampoco se vería afectada la cabecera del Ayo. de Rejas, situada a más de 100 m de la estructura.
Área del solar (m2) - Área del Pozo equivalente Sa 2400.00 m2
Relación aspecto L/a (*) L/a 6.00 --- Igualmente al escenario anterior, y como puede observarse en la primera de las tablas, la relación de aspecto de
(*) Cuando L/a > 4, ha de calcularse como una Zanja Drenante. la excavación aconseja utilizar la expresión correspondiente a la Zanja Equivalente.

Tabla 1. Parámetros utilizados para el cálculo de caudales en el recinto de la excavación, según la Cálculo mediante la aplicación de la fórmula de la Zanja Equivalente
aplicación de la Fórmula del Pozo equivalente.
De la misma forma, en las siguientes tablas se resumen los parámetros utilizados y los resultados obtenidos de la
aplicación de la Fórmula de la Zanja equivalente a cada uno de los laterales de 120 m. Se ofrece también el
caudal por metro lineal de ambos laterales, en conjunto:
Radio de influencia Abrev. Valor Unidades

Radio Influencia captación genérica (Formula Sichart) Rp 3.37 m Parámetro Abrev. Valor Unidades

Radio equivalente del solar r 27.64 m Espesor saturado estático del acuífero (m) Ho 6.85 m

Radio de influencia del Pozo Equivalente Ro 27.84 m Nivel piezométrico a conseguir (m) H 3.35 m

Ratio r/(R/3) (*) r/(R/3) 17.80 --- Descenso producido (m) d 3.5 m
(*) En grandes excavaciones donde el radio equivalente (r) puede superar la 3ª parte del radio de influencia del pozo (Rp = Coef. permeabilidad (m/s) K 1.03E-07 m/s
r/(Ro/3) > 1), el cálculo del caudal se realiza según la expresión para grandes solares.
Coef. permeabilidad (m/día) K 8.90E-03 m/día
Tabla 2. Radios de influencia del agotamiento, según la Fórmula del Pozo equivalente. Transmisividad (m2/s) T 7.06E-07 m2/s

Transmisividad (m2/día) T 6.10E-02 m2/día

Longitud de la zanja (m) L 120.00 m

Caudal estimado Abrev. Valor Unidades Tabla 4. Parámetros utilizados para el cálculo de caudales a drenar en el recinto de la excavación, según
Expresión general: la aplicación de la Fórmula de la Zanja equivalente (laterales de 120 m).
Caudal captado (m3/s) Q 0.002 m3/s
Caudal captado (m3/día) Q 135.29 m3/día
( 2 H−d ) · d Radio de influencia Abrev. Valor Unidades
Q=π·K·
( )
Caudal captado (L/s) R0 Q 1.57 L/s
Distancia de Influencia captación (Formula Sichart) RZ 1,69 m
ln
Caudal captado (L/día) r Q 1.35E+05 L/día Tabla 5. Radios de influencia del agotamiento, según la aplicación de la Fórmula de la Zanja equivalente.
Caudal captado (m3/s) Expresión aproximada: Q 0.002 m3/s
Valor Valor
Caudal captado (m3/día) d Q 181.71 m3/día Caudal estimado para L = 120 m Abrev. Unidades
Q ≈ 2· π·H·K· 1 sólo lateral 2 laterales (x2)

( )
R
Caudal captado (L/s) Q 2.10 L/s
Caudal captado (m3/s) Q 0.0001 0.0003 m3/s
ln 0 Expresión general:
Caudal captado (L/día) r Q 1.82E+05 L/día
Caudal captado (L/s) Q 0.131 0.262 L/s
3
Caudal captado (m /s) Expresión para grandes solares: Q 0.0002 m3/s K
3
Caudal específico (m /s x m.l.) Q=L· · ( 2 H−d ) · d q 0.000001 0.000002 m3/s x m.l.
3
Caudal captado (m /día) Q 16.63 3
m /día 2 RZ
Q=π·K· ( 2 H−d ) · Caudal específico (L/s x m.l.) q 0.001 0.002 L/s x m.l.

(2r
)
Caudal captado (L/s) Q 0.19 L/s
· d· +0,25 Tabla 6. Estimación de caudales a drenar en el recinto de la excavación, según la aplicación de la Fórmula
Caudal captado (L/día) RP Q 1.66E+04 L/día de la Zanja equivalente (laterales de 120 m).
Por su parte, en las siguientes tablas se resumen los parámetros utilizados y los resultados obtenidos de la [Link]. Cálculo mediante modelo matemático con SEEP/W
aplicación de la Fórmula de la Zanja equivalente a cada uno de los laterales de 20 m. Se ofrece también el caudal
por metro lineal de ambos laterales, en conjunto: De forma paralela a los anteriores cálculos analíticos, se ha procedido a realziar una estimación de caudales de
infiltración en terreno natural, con las características de la Peñuela (T P), según la parametrización expuesta en el
Parámetro Abrev. Valor Unidades apdo. [Link].3. En este primer escenario, y de forma coherente con la asunción de los cálculos anteriores, se ha
asumido que la totalidad que los caudales susceptibles de producirse en obra van a provenir de la segunda capa
Espesor saturado estático del acuífero (m) Ho 6.85 m
terciaria considerada, integrada, como se comentaba, por las facies de la Peñuela (TP).
Nivel piezométrico a conseguir (m) H 3.35 m

Descenso producido (m) d 3.5 m De esta forma, queda una malla constituida por una única capa equivalente a la CAPA 1 del modelo genérico
construido:
Coef. permeabilidad (m/s) K 1.03E-07 m/s

Coef. permeabilidad (m/día) K 8.90E-03 m/día POTENCIA PERMEABILIDAD RELACIÓN COEF.

CAPA MATERIAL
2 2 (m) HORIZ. KH (m/s) KV/KH ALMACEN.
Transmisividad (m /s) T 7.06E-07 m /s
Terreno
Transmisividad (m2/día) T 6.10E-02 m2/día CAPA 2 > 30 m Arcillas 1,03 x 10-07 0,5 0,02 (2%)
natural
Longitud de la zanja (m) L 20.00 m Tabla 13. Valores iniciales de referencia para los parámetros hidrogeológicos del modelo
Tabla 7. Parámetros utilizados para el cálculo de caudales a drenar en el recinto de la excavación, según
la aplicación de la Fórmula de la Zanja equivalente (laterales de 20 m).

En las siguientes figuras se ilustran tanto la constitución de la malla inicial como el resultado de la modelización.
Radio de influencia Abrev. Valor Unidades
Éstas han sido incorporadas también al Apéndice 4.
Distancia de Influencia captación (Formula Sichart) RZ 1,69 m
Tabla 8. Radios de influencia del agotamiento, según la Fórmula de la Zanja equivalente.

Valor Valor
Caudal estimado Abrev. Unidades
1 sólo lateral 2 laterales (x2)
Caudal captado (m3/s) Q 0.00002 0.00004 m3/s
Expresión general:
Caudal captado (L/s) Q 0.022 0.044 L/s
K
Caudal específico (m3/s x m.l.) Q=L· · ( 2 H−d ) · d q 0.000001 0.000002 m3/s x m.l.
2 RZ
Caudal específico (L/s x m.l.) q 0.001 0.002 L/s x m.l.
Tabla 9. Estimación de caudales a drenar en el recinto de la excavación, según la aplicación de la Fórmula
de la Zanja equivalente (laterales de 20 m).

Como puede comprobarse, entre los cuatro laterales se obtiene un caudal total de 0,31 L/s para todo el recinto
de la excavación. Por su parte, se ha estimado un caudal específico para cada metro de zanja, de 0.001 L/s por m.
l. y un radio de influencia máximo perimetral de la excavación de aproximadamente 1,44 m, por lo que el punto
de agua 1922-4-0114 quedaría igualmente fuera del área de influencia del posible drenaje. Tampoco se vería
afectada la cabecera del Ayo. de Rejas, situada a más de 100 m de la estructura.
Figura 7. Malla del perfil transversal objeto de modelización en la hipótesis simple de Peñuela (elaboración
propia).
 CAPA 3

CAPA 1

CAPA 4
CAPA 2

Figura 8. Resultados de la modelización de flujo en el perfil transversal objeto de modelización en la Figura 9. Malla del perfil transversal objeto de modelización (elaboración propia).
hipótesis simple de Peñuela (elaboración propia).
POTENCIA PERMEABILIDAD RELACIÓN COEF.
CAPA MATERIAL
(m) HORIZ. KH (m/s) KV/KH ALMACEN.
Arenas finas –
CAPA 1 2,76. 1,6 x 10-05 0,5 0,32 (32%)
3 Terreno medias, arcillosas
Como puede apreciarse en la última de ellas, la sección de flujo arroja un valor de 2.0177E-7 m /s y por metro
natural
lineal de sección transversal, que equivalen a 2.0177E-04 L/s (0.00020177 L/s). Considerando la longitud total de CAPA 2 > 30 m Arcillas 1,03 x 10-07 0,5 0,02 (2%)

la estructura asumida en los cálculos analíticos (120 m), el caudal obtenido según la modelización resulta ser de CAPA 3 Hormigonado
5,33 x 10-6 (Tosco) 0,5 0,11 (10,66%)
(PILOTES) Material (permeabilización
0,0242 L/s. Este resultado es aproximadamente 1 orden de magnitud inferior que los caudales más bajos 1 m.
CAPA 4 artificial 33,3% de la
obtenidos según métodos analíticos (entre 0,16 y 2,10 L/s), especialmente al obtenido por medio de la aplicación 3,43 x 10-8 (Peñuela) 0,5 0,0067 0,67%)
(PILOTES) sección).
de la Fórmula de la Zanja equivalente (0,26 L/s) y de la expresión para grandes excavaciones del Pozo equivalente Tabla 14. Valores iniciales de referencia para los parámetros hidrogeológicos del modelo
(0,16 L/s). No obstante, puede considerarse la estimación mediante modelización matemática de flujo más
precisa puesto que incluye en el cálculo otros parámetros hidrogeológicos (como el coeficiente de
almacenamiento, la succión, etc.) que no se contemplan en las fórmulas analíticas.

3.2.4. ESTIMACIÓN CONJUNTA DE MATERIALES DE TOSCO ARENOSO (TTA) Y PEÑUELA (TP)

Por último, se ha abordado la simulación de un modelo más completo de la sección de excavación, según la
descripción expuesta en el apdo. [Link].3., con dos capas de terreno y dos estimaciones de pantallas de pilotes
con cierta permeabilidad, según el tipo de terreno.

En la tabla siguiente se recuerdan la parametrización de las características hidrogeológicas establecidas para

cada una de las 4 capas. Asimismo, en las siguientes figuras se ilustran la constitución de la malla inicial como el
resultado de la modelización. Éstas han sido incorporadas también al Apéndice 4.

Figura 10. Resultados de la modelización de flujo en el perfil transversal objeto de modelización en la

hipótesis simple de Peñuela (elaboración propia).
Como se puede ver en los resultados de la simulación, el caudal de infiltración obtenido es de 1,654E-7 m 3/s, que
equivale a 1,654E-4 L/s (0.0001654 L/s). En el conjunto de los 120 m estimados para la estructura, se obtienen
aproximadamente 0.0198 L/s (≈ 0,02 L/s). Como resulta lógico, este resultado se corresponde en orden de
magnitud con el estimado para la hipótesis simple del drenaje a través de las facies Peñuela, dado que, según las
cotas obtenidas para el nivel freático, es precisamente a través de esta unidad por donde se efectuará
mayoritariamente el drenaje subterráneo.

3.2.5. DISCUSIÓN DE LOS RESULTADOS

Las expresiones utilizadas para la estimación de los caudales susceptibles de producirse durante la excavación,
resultan extremadamente sensibles a dos parámetros: el valor de la permeabilidad utilizado y las dimensiones
del solar. Cabe tener en cuenta que, ninguna de ellas tiene en consideración otros parámetros hidrogeológicos
aparte de la permeabilidad.

En el primero de ellos, la variación de un orden de magnitud en la permeabilidad se traduce una variación similar
en los caudales. Así, por ejemplo, en el escenario de drenaje a través de las facies de Tosco arenoso (T TA), y
empleando la permeabilidad máxima de K = 3,38E-04 m/s (obtenida según la aplicación de la Fórmula de Hazen
para el SE-24), la estimación de caudales totales según las fórmulas del Pozo equivalente oscilaría entre 20 y 30
L/s. En el caso de la Zanja equivalente, ésta arrojaría valores en torno a 17,5 L/s. En este último caso, se obtienen
además caudales específicos de 0,062 L/s x m. l. (multiplicados por 280 m perimetrales, se obtiene un caudal
total de 17,36 L/s.).

El grado de incertidumbre de las estimaciones analíticas, supone ya un orden de magnitud completo.

Por lo que respecta al segundo parámetro, y considerando la aplicación del Fórmula de la Zanja equivalente, la
asunción de un recinto cuadrado de 120 x 120 m (algo implícito en la aproximación según el Pozo equivalente),
en lugar del rectangular considerado (120 x 20 m), aumenta los caudales infiltrados en torno a medio orden de
magnitud (sobre 6,6 L/s).

En cuanto a las estimaciones a partir del modelo matemático bidimensional de flujo, cabe recordar lo ya
mencionado anteriormente: éste incluye en el cálculo otros parámetros hidrogeológicos (como el coeficiente de
almacenamiento, la succión, etc.) que no se contemplan en las fórmulas analíticas, por lo que (respaldado por la
convergencia con los órdenes de magnitud de las fórmulas analíticas) puede considerarse algo más precisa y
representativa.

En el caso de utilizar la permeabilidad máxima estimada (3,38E-4 m/s) para el escenario de drenaje a través del
Tosco arenoso (TTA), el caudal unitario por sección resultaría de 0,00062666 m 3/s (0.62666 L/s) que, para toda la
estructura, se traduce en 75 L/s.