Aplicación y limitaciones de un modelo 2D de transporte de calor para estimar la recarga del río Guadalfeo al acuífero Motril-Salobreña

Abstract

Un modelo 2D de transporte de calor fue aplicado para determinar la interacción entre en acuífero Motril-Salobreña y el río Guadalfeo. La infiltración del río al acuífero se estimó a partir de este método. La infiltración máxima que se simuló fue un 18 % del caudal del río al inicio de los periodos de recarga (Marzo y Noviembre) y cercano a 0 relacionado con los momentos de nivel piezométrico elevado (Abril y Mayo). Las diferencias entre los resultados del modelo y los perfiles de temperaturas medidas en pozos cercanos al río fueron analizados para determinar las limitaciones del método de simulación de transporte de calor. Las principales dificultades para ajustarlo a las observaciones se localizaron a 20-30 metros de profundidad y cerca de la superficie topográfica. Tanto en las zonas profundas como en las cercanas a la superficie sería de gran utilidad incluir un flujo perpendicular a la sección simulada para incrementar la velocidad de los cambios y evitar la sobreestimación de la infiltración. Sin embargo esto no fue posible al utilizar un modelo en dos dimensiones

A 2D heat transport model was applied to determine the interaction between the Motril-Salobreña aquifer and the Guadalfeo River. The river infiltration to the aquifer was estimated from this method. The maximum simulated infiltration was 18 % of the total river discharge at the beginning of the recharge period (March and November), and close to 0 related to the timing of the increase of the water table in the aquifer (April, May). The differences between the model results and the observed temperatures logs in wells located in the vicinity of the river were analyzed looking for the limitations of the heat transport simulation method. The main difficulties matching observations were located at 20-30 meters depth and close to the surface. In these shallow and deeper zones, it would be useful to include perpendicular flow (underflow) to the simulated section for increasing the velocity of the changes and avoid the infiltration overestimation. However, this was not possible with a 2D modeling approach