Estimaciones de precipitaciones - Centro de Estudios Hidrográficos ...
Estimaciones de precipitaciones - Centro de Estudios Hidrográficos ...
Estimaciones de precipitaciones - Centro de Estudios Hidrográficos ...
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