Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 147 vistas

    Calibración - ModeloGR4 J (4)

    Cargado por

    Mateo López Marín
    Este documento describe el proceso de calibración del modelo hidrológico GR4J para la estación La Boca. Se descargó el modelo, se ingresaron los datos de precipitación, evapotranspiración y caudales, y se filtraron los datos faltantes. Luego, se ajustaron manualmente y automáticamente los parámetros del modelo para maximizar la eficiencia de Nash, obteniendo valores finales de 56,25% para la calibración y 34,7% para la validación. El comportamiento irregular de la relación precipitación-caudal se atribuye a

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd

    Calibración - ModeloGR4 J (4)

    Cargado por

    Mateo López Marín
    147 vistas13 páginas
    Este documento describe el proceso de calibración del modelo hidrológico GR4J para la estación La Boca. Se descargó el modelo, se ingresaron los datos de precipitación, evapotranspiración y caudales, y se filtraron los datos faltantes. Luego, se ajustaron manualmente y automáticamente los parámetros del modelo para maximizar la eficiencia de Nash, obteniendo valores finales de 56,25% para la calibración y 34,7% para la validación. El comportamiento irregular de la relación precipitación-caudal se atribuye a

    Descripción original:

    Ejemplo de calibración model GR4-J. Calcula y calibra modelos

    Título original

    Calibración_ModeloGR4-J-{4}

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Este documento describe el proceso de calibración del modelo hidrológico GR4J para la estación La Boca. Se descargó el modelo, se ingresaron los datos de precipitación, evapotranspiración y caudales, y se filtraron los datos faltantes. Luego, se ajustaron manualmente y automáticamente los parámetros del modelo para maximizar la eficiencia de Nash, obteniendo valores finales de 56,25% para la calibración y 34,7% para la validación. El comportamiento irregular de la relación precipitación-caudal se atribuye a

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    147 vistas13 páginas

    Calibración - ModeloGR4 J (4)

    Cargado por

    Mateo López Marín
    Este documento describe el proceso de calibración del modelo hidrológico GR4J para la estación La Boca. Se descargó el modelo, se ingresaron los datos de precipitación, evapotranspiración y caudales, y se filtraron los datos faltantes. Luego, se ajustaron manualmente y automáticamente los parámetros del modelo para maximizar la eficiencia de Nash, obteniendo valores finales de 56,25% para la calibración y 34,7% para la validación. El comportamiento irregular de la relación precipitación-caudal se atribuye a

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    de 13

    Taller # 3.

    Calibración modelo GR4J

    Grupo 4: Mateo López Marín - Cristian Aristizábal Grajales

    El primer paso, es descargar el modelo GR4J en formato xlsx de la siguiente página Daily
    hydrological model GR4J – Équipe Hydrologie des bassins versants (inrae.fr)

    Una vez descargado el modelo, se traen todos los datos diarios calculados e interpolados para la
    precipitación (Pluie mm), el valor de la evapotranspiración potencial (ETP mm) y los caudales
    registrados en la estación limnimétrica La Boca (Débit m3/s), hallados previamente en el Taller N°2.
    Una vez plasmada la información en el modelo, se identifican los datos donde la estación limnimétrica
    no registró caudales y se elimina completamente la información asociada a ese día sin registro.

    Con los registros filtrados, procedemos a asignarle al modelo el nombre de la estación de estudio y
    el área de la cuenca en km2 relacionada a dicha estación. También, cuantificamos los períodos;
    asignamos 1 año (365 días) para el período de calentamiento, 17 años aproximadamente (6103 días)
    para el período de calibración y se deja un último año para el período de validación.

    Con la información entrante suministrada, se observa el valor del Nash promedio de todos los
    criterios de eficiencia para comenzar a calibrar el modelo.

    Se deben variar y ajustar en un rango especifico los “Parámetros del modelo” x1, x2, x3 y x4 hasta
    encontrar la mayor eficiencia promedio.

    Para graduar cada parámetro se define el rango de ajuste, el cual iniciará para todos los parámetros
    desde 0 hasta el doble de las condiciones iniciales suministradas por el modelo.

    x1: [0 hasta 11.54] x2: [0 hasta 3.24]

    x3 [0 hasta 8.48] x4: [-0.24 hasta 0]


    Una vez definido el rango, se construye la curva de calibración de cada parámetro. Para esto fue
    necesario asignarle un incremento o paso, en el caso de x1, x2 y x3 fue de 0.1 y x4 de 0.01. Luego,
    se debe evaluar en el modelo el valor de cada parámetro y calcular su Nash promedio
    correspondiente.

    Dado que el volumen de valores a evaluar se vuelve difícil e ineficiente de administrar manualmente,
    se automatiza el proceso de evaluación de los parámetros para una calibración manual más rápida.

    Esta automatización se realizó con códigos de programación (VBA) que permiten calcular el Nash
    promedio para cada paso en el rango del parámetro.

    Iniciamos, creando un botón de ejecución.


    Posteriormente, se crea un formulario con 4 botones para parámetro

    Desarrollamos las líneas de código para automatizar el proceso


    Se comienza con la calibración de cada parámetro

    Finalmente, se elaboran las curvas de calibración y se selecciona el valor del parámetro que arroja
    una mayor eficiencia de Nash.

    Calibración x1
    Calibración x2
    Calibración x3
    Calibración x4
    Por últimos, se realiza una calibración automática por medio de la herramienta “Solver” de Excel.

    Seleccionamos como objetivo maximizar el valor de Nash promedio variando los parámetros x1, x2,
    x3 y x4 del modelo.
    Luego de realiza la calibración, la gráfica de caudal observado versus el caudal simulado arrojó los
    siguientes resultados:

    Para la validación del modelo, se usó el año restante entre el 01/01/2014 y el 31/12/2014. Donde se
    obtiene que, para los parámetros de x1 al x4 calibrados anteriormente, la eficiencia del modelo según
    el Nash promedio es de 56,25.

    En la siguiente imagen, se ilustran los resultados correspondientes a la validación:


    Modélisation pluie-débit à l'aide du modèle journalier GR4J (Version Perrin et al., 2003)

    (Rq: seules les cellules en jaune sont à modifier, les autres sont issues de calculs; les formules des cellules G41 à BQ41
    sont à recopier sur les lignes suivantes pour les calculs des pas de temps suivants)

    Nom du bassin La Boca

    Superficie du bassin (km²) 47.63

    Paramètres modèle Transf. Réels


    x1: Capacité rés. production (mm) 6.48 649.22
    Calibración Manual
    x2: Paramètre d'échange (mm) 2.44 5.71
    x3: Capacité rés. routage (mm) 4.07 58.85
    x4: Délai (jours) -0.20 1.32

    Valeurs initiales
    Taux de remplissage initial S0/x1 0.60
    Taux de remplissage initial R0/x3 0.70

    Période
    Longueur de la période de mise en route (j) 6469
    Durée de la période test (j) 6832
    Date de départ 1/01/2014
    Date de fin 31/12/2014

    Moyenne des pluies observées (mm/j) 8.77


    Moyenne des ETP observées (mm/j) 3.02
    Moyenne des débits observés (mm/j) 12.31
    Moyenne des racines des débits observés 3.07
    Moyenne des log des débits observés 2.07

    Critères d'efficacité (%) Nash Promedio


    Nash(Q) 34.7 56.25
    Nash(VQ) 61.0
    Nash(ln(Q)) 73.0
    Bilan 89.6

    Conclusiones:

    • Por medio de la calibración manual se logró identificar que, para este modelo, el parámetro
    x2 es el que más variabilidad genera en la eficiencia de este. Moviendo el Nash en un rango
    entre -13 y 21.
    • Por su parte, el parámetro x3 tiene un comportamiento similar al anterior, variando la
    eficiencia del modelo entre magnitudes de -3 y 21.
    • x1 genera menor variabilidad que los anteriores, pero sigue siendo determinante en el
    calculo de la eficiencia del modelo. Variando la magnitud del Nash entre -2 y 16.
    • Por último, el parámetro x4 es el que muestra menor relevancia en la definición en la
    eficiencia del modelo, ya que la magnitud del Nash varia en un rango de 0.15.
    • Con el fin de evaluar la poca eficiencia que se observa en el modelo, se hizo un análisis
    exploratorio de los datos de precipitación promedio respecto al caudal observado y, se
    encontró un comportamiento irregular en la relación de estos donde, para valores altos de
    precipitación se obtienen valores medios y altos de caudal, mientras que; en algunas
    temporadas se registraron valores picos de caudal aun cuando se presentaron
    precipitaciones de baja intensidad o simplemente no se presentaron.
    • Ese comportamiento irregular se puede explicar en que la estación limnimétrica La Boca se
    encuentra sobre el cauce del río Batatas, el cual es uno de los 3 afluentes que alimentan el
    embalse del Guavio (Rio Guavio, Rio Batatas y rio Chivor). Tanto el rio Batatas como el
    Chivor alimentan el embalse por medio un trasvase, el cual, para el caso del Rio Batatas se
    da aguas arriba de la estación de medición. Esta condición externa de toma de caudal no
    permite que el modelo represente fielmente el comportamiento de la cuenca, mucho menos
    cuando dicha condición externa no se está considerando, ya que no se tiene en cuenta el
    caudal que se le está arrebatando a esta.

    También podría gustarte