03.06.2013 Views

Estimaciones de precipitaciones - Centro de Estudios Hidrográficos ...

Estimaciones de precipitaciones - Centro de Estudios Hidrográficos ...

Estimaciones de precipitaciones - Centro de Estudios Hidrográficos ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Combinación <strong>de</strong> fuentes <strong>de</strong> información para la<br />

estimación <strong>de</strong> <strong>precipitaciones</strong><br />

Director <strong>de</strong> programa<br />

<strong>Centro</strong> <strong>de</strong> <strong>Estudios</strong> <strong>Hidrográficos</strong> <strong>de</strong>l CEDEX<br />

Angel Luis Aldana Valver<strong>de</strong><br />

Dr. Ingeniero <strong>de</strong> Caminos, Canales y Puertos<br />

Estimación <strong>de</strong> <strong>precipitaciones</strong><br />

Coordinador <strong>de</strong> la red temática<br />

PROHIMET<br />

<strong>de</strong>l programa CYTED


<strong>Estimaciones</strong> <strong>de</strong> <strong>precipitaciones</strong><br />

• A partir <strong>de</strong> pluviómetros<br />

– Inversa a la distancia elevado a un exponente<br />

• Con restricciones<br />

• Con reparto volumétrico<br />

• A partir <strong>de</strong> radar<br />

• Combinación<br />

• Evaluaciones <strong>de</strong> coherencia<br />

– Por intervalo<br />

• Rechazo <strong>de</strong> valor en pluviómetro<br />

– Acumulado<br />

Estimación <strong>de</strong> <strong>precipitaciones</strong>


Interpolación <strong>de</strong> <strong>precipitaciones</strong><br />

• Interpolación para obtener el valor <strong>de</strong> la precipitación P en el punto (x, y):<br />

P<br />

( x,<br />

y)<br />

∑<br />

l = 1<br />

f l (x,y): funciones base<br />

• Fórmula general <strong>de</strong> interpolación en función <strong>de</strong> los valores puntuales:<br />

=<br />

p<br />

a<br />

Estimación <strong>de</strong> <strong>precipitaciones</strong><br />

l<br />

*<br />

f<br />

l<br />

( x,<br />

y)<br />

n<br />

n ⎛<br />

P( x,<br />

y)<br />

= ∑ λ<br />

* ; ⎜ ∑<br />

l Pl<br />

λ l<br />

l = 1 ⎝ l = 1<br />

=<br />

⎞<br />

1⎟<br />


Pon<strong>de</strong>ración inversa a la distancia<br />

• Expresión general <strong>de</strong> los coeficientes <strong>de</strong> pon<strong>de</strong>ración, cumpliéndose<br />

que su sumatorio es la unidad:<br />

λ<br />

l<br />

=<br />

ϖ<br />

l<br />

n<br />

∑ ϖk<br />

k = 1<br />

• Expresión <strong>de</strong> funciones <strong>de</strong> pon<strong>de</strong>ración <strong>de</strong>l método <strong>de</strong> inversa a la<br />

distancia h l (x,y) al pluviómetro l <strong>de</strong> coor<strong>de</strong>nadas (x l , y l ) elevado a un<br />

exponente θ :<br />

ϖ<br />

l<br />

=<br />

θ<br />

1<br />

1<br />

=<br />

h ( x,<br />

y)<br />

−<br />

l<br />

( ) θ<br />

2<br />

2<br />

( x − x ) + ( y y )<br />

Estimación <strong>de</strong> <strong>precipitaciones</strong><br />

l<br />

l


Restricciones <strong>de</strong> distancia<br />

• Distancia máxima (d max ):<br />

– Se or<strong>de</strong>nan los pluviómetros según el valor <strong>de</strong> la distancia h l y se<br />

seleccionan aquellos m pluviómetros que distan menos que d max <strong>de</strong>l punto<br />

en el que obtener el valor interpolado.<br />

• Número máximo <strong>de</strong> pluviómetros (n max ):<br />

– Si m> n max se seleccionan los primeros n max.<br />

Estimación <strong>de</strong> <strong>precipitaciones</strong>


Método <strong>de</strong>l reparto volumétrico<br />

• Método basado en una fórmula <strong>de</strong>l tipo general pero sobre valores <strong>de</strong><br />

precipitación en pluviómetro modificados, es <strong>de</strong>cir, suponiendo que tienen<br />

asociado un cierto error:<br />

*<br />

Pl<br />

n<br />

n<br />

* ⎛<br />

P( x,<br />

y)<br />

= ∑ λl<br />

* Pl<br />

; ⎜∑<br />

λl<br />

l=<br />

1 ⎝ l=<br />

1<br />

: valor corregido en cada pluviómetro<br />

Estimación <strong>de</strong> <strong>precipitaciones</strong><br />

⎞<br />

= 1⎟<br />


Expresión <strong>de</strong> valores corregidos<br />

• Los valores pluviométricos se corregirán mediante<br />

P<br />

*<br />

l<br />

=<br />

n<br />

∑<br />

l=<br />

1<br />

con coeficientes que tengan una expresión <strong>de</strong>l tipo:<br />

c<br />

c<br />

k<br />

l<br />

k<br />

l<br />

=<br />

=<br />

a<br />

l<br />

(C<br />

1<br />

k<br />

( h )<br />

l<br />

auto<br />

a<br />

*D<br />

c<br />

k<br />

l<br />

k<br />

*<br />

P<br />

Estimación <strong>de</strong> <strong>precipitaciones</strong><br />

θ<br />

l<br />

si<br />

min<br />

)<br />

θ<br />

≠<br />

l<br />

si<br />

k<br />

k<br />

=<br />

l


Coeficiente <strong>de</strong> autopeso<br />

• Don<strong>de</strong> se ha empleado la distancia mínima entre pluviómetros:<br />

y un parámetro C auto que se <strong>de</strong>nomina coeficiente <strong>de</strong> autopeso.<br />

• Como:<br />

D<br />

min<br />

n<br />

∑<br />

k = 1<br />

=<br />

P<br />

*<br />

k<br />

min<br />

{ k<br />

h ; k = 1,...,<br />

n;<br />

l = 1,...,<br />

n}<br />

l<br />

k<br />

∑ ∑ cl<br />

*<br />

Pl<br />

= ∑ ⎜∑<br />

= n<br />

n<br />

k = 1 l = 1<br />

l = 1 k = 1<br />

Estimación <strong>de</strong> <strong>precipitaciones</strong><br />

n<br />

⎛<br />

⎝<br />

n<br />

c<br />

k<br />

l<br />

⎞<br />

⎟<br />

⎠<br />

P<br />

l


• La igualdad <strong>de</strong> sumas:<br />

se cumple si:<br />

lo que implica que:<br />

Igualdad <strong>de</strong> sumas<br />

a<br />

l<br />

= n<br />

∑<br />

k = 1;<br />

k ≠ l<br />

n<br />

∑<br />

k = 1<br />

P<br />

n<br />

∑<br />

k = 1<br />

*<br />

k<br />

k<br />

cl<br />

1<br />

( h )<br />

k<br />

l<br />

θ<br />

∑<br />

= n<br />

=<br />

+<br />

l=<br />

1<br />

1<br />

1<br />

(C<br />

P<br />

Estimación <strong>de</strong> <strong>precipitaciones</strong><br />

l<br />

auto<br />

1<br />

*D<br />

min<br />

)<br />

θ


<strong>Estimaciones</strong> basadas en radar<br />

meteorológico<br />

Estimación <strong>de</strong> <strong>precipitaciones</strong>


Combinación <strong>de</strong> fuentes <strong>de</strong><br />

información<br />

• Parámetros <strong>de</strong> usuario<br />

– Fiabilidad<br />

• relaciona la información radar y la información pluviométrica en magnitud, <strong>de</strong><br />

– Forma<br />

tal forma que cuanto mayor es este parámetro, mayor será el peso <strong>de</strong> la<br />

información pluviométrica en perjuicio <strong>de</strong> la información radar.<br />

• relaciona la información radar y la pluviométrica atendiendo a la forma,<br />

siendo el resultado <strong>de</strong> la malla <strong>de</strong> combinación más parecida a la malla <strong>de</strong><br />

interpolación basada en pluviómetros cuanto mayor sea este parámetro.<br />

Estimación <strong>de</strong> <strong>precipitaciones</strong>


Método fiabilidad-forma <strong>de</strong> combinación <strong>de</strong> funciones.<br />

Fiabilidad = 0.5<br />

Forma =0.5<br />

Estimación <strong>de</strong> <strong>precipitaciones</strong>


Método fiabilidad-forma <strong>de</strong> combinación <strong>de</strong> funciones.<br />

Fiabilidad =<br />

0.1<br />

Forma =0.1<br />

Estimación <strong>de</strong> <strong>precipitaciones</strong>


Método fiabilidad-forma <strong>de</strong> combinación <strong>de</strong> funciones.<br />

Fiabilidad = 0.9<br />

Forma =0.9<br />

Estimación <strong>de</strong> <strong>precipitaciones</strong>


Evaluación <strong>de</strong> coherencia<br />

• Índice <strong>de</strong> Coherencia:<br />

coherencia<br />

=<br />

• Rechazo <strong>de</strong> pluviómetro por umbrales <strong>de</strong>:<br />

– Precipitación mínima<br />

– Coherencia mínima<br />

⎛ estimado − real ⎞<br />

⎜ 1− ⎟<br />

⎝ estimado + real ⎠<br />

Estimación <strong>de</strong> <strong>precipitaciones</strong><br />

* 100


EDIMACHI_CEnP<br />

Principales características<br />

• Configuración <strong>de</strong> directorios <strong>de</strong> datos y<br />

aplicaciones<br />

• Consulta y edición <strong>de</strong> datos <strong>de</strong>scriptivos<br />

básicos (curvas <strong>de</strong> gasto, <strong>de</strong> embalses,<br />

conexiones entre elementos)<br />

• Consulta y edición <strong>de</strong> datos temporales<br />

(selección, filtrado, relleno <strong>de</strong> huecos)<br />

• Cálculos básicos, tales como caudales a partir<br />

<strong>de</strong> curvas, volúmenes y balances en embalses.<br />

• Cálculo <strong>de</strong> <strong>precipitaciones</strong>, basados en datos<br />

pluviométricos que, opcionalmente, pue<strong>de</strong>n<br />

combinarse con mallas <strong>de</strong> valores estimados,<br />

como las generadas a partir <strong>de</strong> los datos <strong>de</strong><br />

radares meteorológicos.<br />

• Consulta <strong>de</strong> resultados <strong>de</strong> cálculo realizados<br />

con esta aplicación y con las periféricas.<br />

Estimación <strong>de</strong> <strong>precipitaciones</strong>


EDIMACHI-CEnP.<br />

Edición <strong>de</strong> topología<br />

Estimación <strong>de</strong> <strong>precipitaciones</strong>


EDIMACHI-CEnP.<br />

Consulta <strong>de</strong> datos temporales<br />

Estimación <strong>de</strong> <strong>precipitaciones</strong>


EDIMACHI-CEnP.<br />

Relleno, filtrado, coherencia<br />

Estimación <strong>de</strong> <strong>precipitaciones</strong>


EDIMACHI-CEnP.<br />

Consulta topológica<br />

Estimación <strong>de</strong> <strong>precipitaciones</strong>


Presentación <strong>de</strong> <strong>precipitaciones</strong><br />

Estimación <strong>de</strong> <strong>precipitaciones</strong>


Salida <strong>de</strong> resultados<br />

Estimación <strong>de</strong> <strong>precipitaciones</strong>


EDIMACHI-CEnP. Radar<br />

Zona 1 Zona 2<br />

Estimación <strong>de</strong> <strong>precipitaciones</strong><br />

Zona 3

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!