Pilotes y Pantallas Termoactivas. Proceso Constructivo y Caracterización Térmica
Pilotes y Pantallas Termoactivas. Proceso Constructivo y Caracterización Térmica
Proceso
constructivo y caracterización térmica
ALBERTO MAZARIEGOS (*), LUIS CARLOS ANTÓN (**), RODRIGO VALLE (**), JAVIER URCHUEGUIA (***),
SALVADOR QUILIS (****), SALVADOR MARTÍNEZ (****) y TERESA MAGRANER (****)
RESUMEN La Cimentación Termoactiva es una tecnología aplicable a los elementos de las estructuras de hormigón
armado de las cimentaciones especiales, “pilotes”, “pantallas” y “losas”, para la obtención de energía para climatización de
edificios a partir del subsuelo. Este tipo de cimentación se basa en el aprovechamiento de la temperatura del terreno para
aumentar el rendimiento de las bombas de calor geotérmicas.
En estas estructuras de hormigón armado, el intercambio geotérmico se puede realizar por medio de un circuito cerrado
instalado en las armaduras de la cimentación. Este circuito cerrado esta formado por tubos de PEAD “sondas geotérmi-
cas”, a través de los cuales circula agua o agua con anticongelante, produciéndose un intercambio de calor entre este fluido
y el terreno. El fluido es conducido a una bomba de calor geotérmica “BCG”, generando la energía suficiente para la com-
pleta climatización de un edificio.
ENERGESIS, empresa puntera en el campo de la climatización de edificios a través de la energía geotérmica y RODIO
KRONSA, empresa puntera en el campo de la ingeniería del suelo, con la colaboración de la Universidad Politécnica de Va-
lencia (UPV) y la Universidad Politécnica de Madrid (UPM), han desarrollado una línea de actuación para proyectar y eje-
cutar instalaciones de ahorro y aprovechamiento energético mediante el uso de la energía geotérmica para edificaciones,
asociadas a las cimentaciones.
La necesidad de disponer de un procedimiento de actuación para la ejecución de pilotes y pantallas termoactivas y su ca-
racterización térmica, es el motivo para realizar una serie de pruebas en obra, con la finalidad de verificar el comporta-
miento térmico de estos dos tipos de cimentación termoactiva.
THERMOACTIVE PILES AND DIAPHRAGM WALLS BUILDING PROCESS AND THERMAL RESPONSE TEST
ABSTRACT Thermoactive foundation technology uses building foundation, “piles”, “diaphragm wall” and “slabs” as
ground heat exchangers to provide buildings with heating and cooling. This type of foundation is based on the use of the
ground temperature to enhance the heat pump coefficient of performance (COP).
In this kind of structures, the heat exchange can be done through a closed circuit installed on the reinforced of the founda-
tion. This closed circuit is made up of HDPE pipes “geothermal probes”, through which water or brine circulates producing
a heat exchange between this fluid and the ground. This fluid is led to a geothermal heat pump “GHP”, generating the
energy required for air-conditioning a building.
ENERGESIS, leading company in the air-conditioning field using geothermal energy and RODIO KRONSA, leading com-
pany in the engineering field of soil and subsoil, with the collaboration of the Universidad Politécnica de Valencia (UPV)
and the Escuela Universitaria de Ingeniería Técnica de Obras Públicas (UPM), have developed an action line to design and
implement savings energy facilities by using geothermal energy associated with the foundations.
A series of tests at work have been done in order to have a procedure to implement the thermal response of thermoactive pi-
les and diaphragm walls.
Palabras clave: Pilote termoactivo, Pantalla termoactiva, Caracterización térmica de pilotes y pantallas,
Bombas de calor geotérmico.
Keywords: Energy piles, Thermoactive diaphragm wall, Thermoactive foundations thermal response test,
Ground coupled heat pump system.
(*) Departamento de Ingeniería Civil: Tecnología de la Construcción. E.U.I.T.O.P. Universidad Politécnica de Madrid.
(**) RODIO KRONSA.
(***) Departamento de Física Aplicada, Universidad Politécnica de Valencia.
(****) Energesis Ingeniería.
Nº de sondas ML de armadura
Zona Nº de pilotes Nº de sondas ML de sondas
por pilote por pilote
Hangar 2 25 2 50 14 700
Área 3 22 2 44 14 616
TABLA 2. Distribución de
pilotes y sondas geotérmicas. TOTAL 243 486 6.804
2.2. PILOTE TERMOACTIVO PREFABRICADO HINCADO. La capa superficial de la parcela, está ocupada por un paquete
PROCESO CONSTRUCTIVO de rellenos antrópicos heterogéneos compuestos por una trama
granosoportada de gravas, bolos y bloques poligénicos de caliza,
La prueba con el pilote termoactivo prefabricado hincado ofita y escorias de fundición con matriz arenolimosa. Infraya-
han tenido dos fases diferenciadas. La primera, la fabrica- cente al nivel de rellenos se describen depósitos aluviales del
ción de 8 pilotes de sección cuadrada, CK 300, en la factoría Cuaternario, constituidos por depósitos finos, arcillas y limos are-
de KRONSA en Algete y la segunda, la hinca de cuatro de nosos, y gruesos, gravas con arenas y limos, del río Urumea, Ta-
estos pilotes en la parcela SU - 24.2, propiedad de AMR Re- bla 3. El sustrato rocoso (Cretácico Superior), está constituido por
fractarios Hernani S.A., en el polígono industrial Ibaiondo, margas y lutitas calcáreas limosas que alternan en bancos centi-
Hernani (Guipúzcoa), donde KRONSA estaba realizando la métricos a métricos con calizas arenosas y areniscas calcáreas.
cimentación de varias naves con este tipo de pilotes.
2.2.2. Características del pilote termoactivo prefabricado
El pilote termoactivo prefabricado, donde se han realizado las
2.2.1. Naturaleza y características del subsuelo medidas de conductividad térmica, es un “pilote prefabricado
Los materiales naturales presentes en la zona están repre- de sección cuadrada CK-300”, habitualmente utilizado por la
sentados por depósitos aluviales del Cuaternario en contacto empresa KRONSA, especialista en la realización de este tipo
discordante sobre una alternancia de margas y calizas are- de cimentaciones con pilotes prefabricados. Para la prueba
nosas del Flysch detrítico-calcáreo del Cretácico Superior, de prevista se fabricaron 8 pilotes (2 pilotes de cada tipo) en la
la Unidad Estructural San Sebastián, en la zona noreste de factoría de KRONSA en Algete (Madrid). Las características
la Cuenca Vasco-Cantábrica. de estos pilotes se reflejan en la Tabla 4.
Nº de pilotes fabricados 2 2 2 2
TABLA 4. Características del pilote termoactivo prefabricado. *Fecha de fabricación: 4 de febrero de 2008. *Fecha de hincado 3 Marzo de 2008.
FIGURA 4. Detalle de la fabricación de los pilotes y los pilotes termoactivos almacenados fuera de los moldes.
La fabricación de cuatro tipos de pilote, diferenciando úni- La hinca de cuatro (4) pilotes termoactivos (uno de cada
camente la distancia del pie de sonda al pie del pilote “azu- tipo fabricado, creando un encepado con una distancia entre
che”, se ha realizado con el fin de comprobar si durante la centros de pilote de 1 m), se llevó a cabo con un equipo JUNT-
hinca del pilote el pie de la sonda geotérmica sufre alguna al- TAN PM 25 HLC 1450, con martillo de 9 Tn. Finalizada la
teración. hinca se procedió al descabezado de los pilotes y verificación de
La fabricación se ha realizado en la Factoría de Pilotes que su estanqueidad, Figura 5.
la empresa KRONSA tiene en Algete (Madrid). Los pilotes
prefabricados CK cuentan con los sellos CE (UNE EN 2.3. PILOTE TERMOACTIVO. CARACTERIZACIÓN DE LA
12794:2006) y N certificado por AENOR que garantiza su fa- RESPUESTA TÉRMICA DEL SUELO
bricación. La caracterización de la respuesta térmica del suelo mediante
Las dimensiones del pilote propuesto, 300 x 300 mm, per- la realización de medidas “in situ” se introdujo en Suecia y Es-
mite la instalación de una sola sonda geotérmica ya que más tados Unidos en 1995, y se ha desarrollado desde entonces en
de una haría que se produjeran interferencias en su comporta- otros países.
miento térmico. Existen diferentes modelos analíticos y numéricos para
En la Figura 4 se observan los detalles de la armadura, ple- evaluar la respuesta térmica del suelo. En [Link] y Canadá
tina de protección y sonda geotérmica en los moldes durante el es frecuente emplear modelos cilíndricos conjuntamente con
proceso de hormigonado, y los pilotes termoactivos almacena- técnicas de estimación de parámetros, mientras que en Eu-
dos, listos para su transporte a obra. ropa el método más empleado, y el que se ha seguido en este
FIGURA 5. Hinca de los pilotes termoactivos y detalle de un pilote descabezado con las sondas geotérmicas.
trabajo, es el de la línea infinita debido a su simplicidad y a y una vez obtenida la conductividad térmica (λ) se puede cal-
la rapidez de su aplicación. cular la resistencia térmica del bucle enterrado (Rb) mediante
El modelo de línea infinita se basa en la aproximación de la expresión (1).
que el intercambiador vertical es una fuente lineal infinita, La caracterización de ambos tipos de pilotes termoactivos
considerando el suelo como un medio homogéneo a una tempe- la ha realizado ENERGESIS mediante el ensayo conocido
ratura inicial constante Tg. El modelo inicial fue desarrollado como “Thermal Response Test” consistente en someter al
por Lord Kelvin; Ingersoll and Plass (1948) aplicaron el mo- terreno, durante un periodo predeterminado, nunca inferior a
delo a los intercambiadores de calor enterrados, mientras que 60 horas, a una inyección o extracción de calor constantes me-
investigaciones posteriores han validado el uso de este modelo diante ensayos de inyección térmica a potencia constante. Este
para la estimación de las propiedades térmicas del suelo (Mo- ensayo simula el comportamiento al que se vería afectado el
gensen, 1983; Kavanaugh, 1984; Austin, 1998; Gehlin, 1998). pilote si se estuviera trabajando en modo refrigeración, mien-
El modelo aproxima el proceso de extracción o inyección de tras que el ensayo de extracción equivale al funcionamiento
calor mediante la fórmula: del sistema en calefacción.
A lo largo del proceso se van registrando telemétricamente
−(qvρc(Tout − Tin )) ⎡ ⎡4at⎤ ⎤ qvρc(Tout − Tin )) Rb las lecturas de temperatura a la entrada y salida del bucle geo-
Tf = ⎢ln⎢ 2 ⎥− γ ⎥+ + Tg (1) térmico. La evolución temporal de dichas temperaturas permite,
4πλH ⎣ ⎣ r0 ⎦ ⎦ H
mediante ajuste al modelo de la línea infinita, obtener paráme-
donde: tros como la conductividad térmica efectiva y resistividad tér-
Tf : temperatura media del bucle enterrado. mica del pilote. Dichos parámetros resultan críticos para conse-
qv: caudal circulante. guir diseños geotérmicos con relaciones óptimas entre la longitud
ρ: densidad fluido circulante. del intercambiador y la energía inyectada o extraída del terreno.
c: calor específico fluido circulante. De esta manera se procede a inyectar un flujo de calor
Tout: temperatura de retorno del bucle enterrado. constante, controlando el caudal de agua inyectado y las tem-
Tin: temperatura ida al bucle enterrado. peraturas a la entrada y salida de las sondas geotérmicas, Fi-
Tg: temperatura estable del suelo. gura 6. Finalmente se obtiene una gráfica que muestra la evo-
lución de la temperatura media entre la entrada y la salida en
λ: conductividad térmica del suelo.
función del tiempo.
H: longitud tubería enterrada.
Los resultados de los ensayos de caracterización térmica de
Rb: resistencia térmica del bucle enterrado.
ambos pilotes se muestran en los siguientes apartados.
γ: constante de Euler.
t: tiempo. 2.3.1. Pilote termoactivo hormigonado “in situ”
ro: diámetro perforación.
k: coeficiente de regresión Tf con ln (t). La caracterización del pilote termoactivo realizado se ha lle-
vado a cabo mediante ensayo de inyección térmica a potencia
a: difusividad térmica.
constante, Tabla 5.
Esta expresión puede emplearse bajo la siguiente condición: El ensayo de inyección ha consistido en someter al pilote
5r02 durante un tiempo determinado a una potencia constante, este
t≥ ensayo simula el comportamiento al que se vería afectado el
a pilote si estuviéramos trabajando en modo refrigeración. De
Por lo que las propiedades térmicas se pueden estimar con esta manera se procede a inyectar un flujo de agua constante a
los datos medidos en el ensayo como: una temperatura determinada por la entrada de las sondas
−(qvρc(Tout − Tin )) geotérmicas (conectadas en serie) registrando el valor de dicha
λ= (2) temperatura y la que se obtiene a la salida de las sondas geo-
4πHk térmicas (ver Figura 6). Así se va registrando la evolución de
Calor Flujo
Tcaliente Tfría
Th Tc
Intercambiador
enterrado
FIGURA 6. Esquema del ensayo de caracterización térmica del suelo y prototipo móvil para realizar medidas de conductividad térmica.
Ensayo Ensayo 1
estas temperaturas a lo largo del tiempo. Finalmente se obten- La pendiente de la gráfica representa el coeficiente de re-
drá una gráfica que muestra la evolución de la temperatura gresión que nos permitirá hallar la conductividad térmica del
media entre la entrada y la salida en función del tiempo. suelo. Se observa que la conductividad del suelo es de:
En la Figura 7, se puede observar la potencia total del pi- K = 1, 02 ± 0,174 W / mK
lote y la potencia por metro de perforación térmica efectiva a
la que se ha sometido el pilote en el ensayo. Se ha inyectado La resistencia térmica asociada al pilote calculada a partir
una potencia de 62,13 W/m. En la Figura 8, se puede ver la de la conductividad obtenida en el terreno es de:
evolución de todos los parámetros de control (temperatura de
Rb = 0,14 K / (W/m)
entrada al pilote, temperatura de salida, temperatura am-
biente, caudal) que se han medido a lo largo del ensayo.
En el análisis de datos se ha despreciado la toma de datos 2.3.2. Pilote termoactivo prefabricado hincado
durante las 6 primeras horas, por considerarse un período no Los ensayos se han realizado considerando que en la tipología
estabilizado del sistema. En la Figura 8, se observan pequeñas de pilotes prefabricados, como tienen poca profundidad efec-
oscilaciones que se deben tanto a interferencias producidas por tiva (unos 12 m.), el intercambio de calor se ve condicionado
fluctuaciones diarias de la temperatura ambiente, como a una por la evolución de la temperatura ambiente (la temperatura
inyección de potencia del terreno no estrictamente constante, del terreno es estable a profundidades mayores de 5 metros,
que presenta una desviación típica de un 7,2 %. en los primeros metros la evolución de la temperatura del te-
En la Figura 9 se representa la temperatura media del rreno queda amortiguada y desfasada respecto a la ambiente).
fluido (Tf) en función del ln del tiempo. La pendiente de la Se han llevado a cabo dos ensayos, uno de inyección y otro
gráfica representa el coeficiente de regresión que nos permi- de extracción, con valores de potencia constante inyectada, Ta-
tirá hallar la conductividad térmica del suelo. bla 6.
Potencia inyectada
940
840
740
640
Potencia
540
340 Potencia/longitud
240
140
40
0 20 40 60 80
FIGURA 7. Potencia inyectada
/longitud del pilote. Horas
40.00 100
Ta entrada
32.00 80
Ta salida
24.00 60
Temperaturas
Ta ambiente
16.00 40
8.00 20
Caudal
0.00 0
0.00 6.25 12.50 19.75 26.00 31.25 37.50 43.75 50.00 56.25 62.50 68.75 76.00 FIGURA 8. Evolución de los
Simulation Time = 75.00 [hr] parámetros de control.
40
35
y = 1,4664x + 11,604
30
Temperatura (°C)
25
20
15
Average Fluid T
(°C)
10
5
FIGURA 9. Temperatura media 10 10,5 11 11,5 12
del fluido/Ln del tiempo. LN (Ts)
30
Ta terreno
25
20
Temperatura (°C)
15
Ta ambiente
10
0
10 20 30 40 50 60 70
FIGURA 10. Ensayo 1. Inyección.
30
25
20
Temperatura (°C)
Ta ambiente
15
10
Ta terreno
0
10 20 30 40 50 60 70
FIGURA 11. Ensayo 2. Extracción . Tiempo (h)
les, a 5 plantas de sótano, e integración, dentro del área de La dimensión de la armadura máxima, 2,85 x 32,50 m, per-
aparcamiento de un espacio para la práctica deportiva de per- mite la instalación del colector geotérmico acoplado a esta en el
sonas discapacitadas. lado terreno. El colector geotérmico instalado en esta armadura
tiene una longitud de 180 metros, formado por 3 bucles de 60 m
3.1.1. Naturaleza y características del terreno de longitud cada bucle. La sonda se colocó a lo largo de los 32,50
Los materiales naturales presentes en la zona están represen- m de longitud de la armadura de la pantalla, Figura 12.
tados por depósitos aluviales de origen marino del Cuaternario El total de paneles con colectores geotérmicos y el total de
en contacto discordante sobre una alternancia de margas y ca- metros lineales instalados se refleja en la Tabla 8. Los paneles
lizas arenosas del Flysch detrítico-calcáreo del Cretácico Supe- o bataches están formados, en general por dos columnas de ar-
rior, de la Unidad Estructural San Sebastián, en la zona no- madura.
reste de la Cuenca Vasco-Cantábrica. La perforación de las pantallas se realizó con una hidro-
El aparcamiento de La Concha se ubica sobre unos rellenos fresa, Evolution III, de Rodio. Se trata de una grúa LIEB-
antrópicos que descansan sobre una terraza aluvial de origen HERR HS 875HD modificada ligeramente a la que se le incor-
marino de naturaleza arenosa. El sustrato rocoso (Cretácico pora, el conjunto de la fresa.
Superior) aparece a 30 m de profundidad junto al hotel Lon- La perforación con hidrofresa, al contrario que el resto de
dres y a 45 m junto al parque Alderdi Eder y está constituido sistemas convencionales, es una actividad continua y perfecta-
por calizas arcillosas y margocalizas, Tabla 7. mente estructurada, que da lugar a un ciclo de trabajo for-
mado por cuatro fases interrelacionadas entre sí:
3.1.2. Características de la pantalla continua 1. Fabricación y almacenamiento del lodo de perforación.
termoactiva 2. Impulsión del lodo hacia el panel de trabajo de la hidro-
La UTE Ferrovial-Moyua contrató a RODIOKRONSA la ejecu- fresa.
ción de la pantalla continua termoactiva del aparcamiento subte- 3. Perforación continúa con hidrofresa e impulsión de re-
rráneo de La Concha en la Plaza de Cervantes, en San Sebastián. torno hacia la planta del material excavado mezclado
Este tipo de pantalla continua ha sido realizada con hidro- con el lodo de perforación.
fresa, de 800 mm de espesor y 40 m de profundidad media, 4. Separación del material excavado y del lodo de perfora-
con una longitud máxima de armadura de 32,50 m. La panta- ción, lodo que a su vez se incorpora a la impulsión hacia
lla se ha realizado de acuerdo con el PEC-24, de RODIO. el panel de trabajo de la hidrofresa.
Temperaturas Temperaturas
–– T_Inlet –– _flujo
–– T_Outlet
–– T_Inletreal
–– T_Outletreal
–– T_ambiente
60.0 1.500
Temperatura salida
24.0 900
0.0 500
0.0 17.9 35.8 53.8 71.7 89.6 107.5’ 125.4 143.3 161.3 179.2 197.1 215.0
Simulation Time = 215.00 [hr]
3.2.1. Validación del modelo realizado dencia entre estas temperaturas y la ambiente, esto es debido
En la gráfica obtenida, Figura 15, se aprecia la evolución de a la longitud de tubería existente entre el laboratorio móvil y
las temperaturas de entrada y salida del colector y el caudal la pantalla termoactiva. Dicha tubería tiene una longitud de
inyectado en la pantalla. En esta se observa una fuerte depen- 25 metros tanto la de ida como la de retorno entre el test y la
Temperaturas Temperaturas
–– T_Inlet –– _flujo
–– T_Outlet
–– T_Inletreal
–– T_Outletreal
–– T_ambiente
50.00
40.00
30.00
Temperaturas
20.00
10.00
–– Potencia
50.00 80
41.00 60
32.00 40
23.00 20
14.00 0
5.00 –20
FIGURA 17. Evolución de la potencia 0.0 17.9 35.8 53.8 71.7 89.6 107.5’ 125.4 143.3 161.3 179.2 197.1 215.0
inyectada a la pantalla. Simulation Time = 215.00 [hr]
pantalla. Para minimizar los efectos de aporte de calor, estas ensayo realizado y se comparan los resultados, los reales y los
tuberías fueron aisladas y protegidas. simulados, Figura 16. De esta manera se van ajustando las va-
A partir de los datos obtenidos se verifica el modelo de riables desconocidas hasta que se obtiene un error entre los re-
transferencia de calor específico desarrollado para este sultados por debajo de un valor preestablecido.
caso, frente el comportamiento real de la pantalla termoactiva En esta gráfica se observa que se produce un buen acople
para dar validez al modelo. entre la curva del modelo y los datos reales obtenidos, con un
El modelo implementado ha tenido en cuenta los aportes de error máximo entre las dos curvas de 0,5°C.
calor generados por la red de tuberías de conexiones entre el Debido al aporte de potencia de la red de tuberías de cone-
laboratorio y las conexiones de la pantalla. Para validar el mo- xión ya comentada, el ensayo realizado ha sido a potencia va-
delo se realiza una simulación con los parámetros reales del riable a lo largo del tiempo. La potencia inyectada a la panta-
50.00 50.00
40.00 40.00
Heat transfer rates
30.00 30.00
Temperaturas
20.00 20.00
10.00 10.00
0.00 0.00
0.0 16.7 33.3 50.0 66.7 83.3 100.0 116.7 133.3 150.0 166.7 183.3 200.0 FIGURA 18. Evolución de la
Simulation Time = 200.00 [hr] temperatura P = 10 W/m.
50.00 50.00
40.00 40.00
20.00 20.00
10.00 10.00
0.00 0.00
FIGURA 19. Evolución de la 0.0 16.7 33.3 50.0 66.7 83.3 100.0 116.7 133.3 150.0 166.7 183.3 200.0
temperatura P = 15 W/m. Simulation Time = 200.00 [hr]
lla, definida como la potencia total del sistema (constante) me- Las simulaciones se llevaron a cabo inyectando un ratio de
nos la potencia de los ramales, se representa en la Figura 17. potencia por metro lineal de colector de 10W/m, 15W/m y
20W/m, durante un periodo de 200 horas. Las evoluciones de
3.2.2. Capacidad térmica de la pantalla las temperaturas de entrada y salida para los ratios de poten-
Una vez validado el modelo, se realizaron tres simulaciones cia inyectada mencionados se representan en las Figuras 18,
para diferentes niveles de potencia. 19 y 20.
50.00 50.00
40.00 40.00
Heat transfer rates
30.00 30.00
Temperaturas
20.00 20.00
10.00 10.00
0.00 0.00
0.0 16.7 33.3 50.0 66.7 83.3 100.0 116.7 133.3 150.0 166.7 183.3 200.0 FIGURA 20. Evolución de la
Simulation Time = 200.00 [hr] temperatura P = 20 W/m.