TFM Gonzalez Vondaniken 2019
TFM Gonzalez Vondaniken 2019
TFM Gonzalez Vondaniken 2019
Presentado por:
D. Boris González von Daniken
Dirigido por:
Dra. IRENE DE BUSTAMANTE
0
I Índice
Índice de Figuras 3
Indice de Tablas 3
II Resumen 1
1. INTRODUCCIÓN 2
2 Objetivo general 4
3 Objetivos específicos 4
4. Marco teórico 5
4.1 ADCP 5
4.2 Caudal medido con ADCP en modalidad dinámica 5
4.3 Incertidumbre 6
4.3.1 Incertidumbre típica combinada uc(y) 6
4.3.2 Incertidumbre expandida U 6
4.3.3 Factor de cobertura k 6
4.3.4 Evaluación Tipo A (errores aleatorios) 7
4.3.5 Evaluación Tipo B (errores sistemáticos) 7
4.3.6 Incertidumbre Combinada 8
4.4 Estimación de la incertidumbre de las mediciones de caudal realizadas con
molinetes (ISO 5168, 2005), (ISO 748, 2007) 8
4.4.1 Incertidumbre debida a errores aleatorios en mediciones con molinete 9
4.4.2 Incertidumbre sistemática en mediciones con molinetes 9
4.5 Incertidumbre en mediciones de Caudal con ADCP en modalidad estática 10
4.6. Ley de propagación de incertidumbres 11
5 Plan de Trabajo y Métodos 14
7 Conclusiones 32
8 Referencias Bibliográficas 34
ANEJOS 36
Índice de Figuras
Figura 8 Perfil de velocidades por verticales del aforo con ADCP en modalidad
estática. .........................................................................................................................31
Figura 9 Río Llobregat en la sección medida, a la izquierda durante la creciente con
gran transporte de sedimentos, a la derecha en la recesión, se aprecia a simple vista
como reduce la turbidez. ...............................................................................................31
Indice de Tablas
resultados del aforo en modo dinámico del día 9 de mayo del 2019 _____________ 26
1
1. INTRODUCCIÓN
En la generación de información hidrométrica de una red de puntos de control, se
realizan mediciones de caudal, para posteriormente calibrar curvas de relación nivel
caudal. Al ser la medición del caudal un parámetro inicial del proceso de producción de
información hidrológica, el resto de los parámetros obtenidos a partir de éste
involucrarán los errores e incertidumbres que hayan sido cometidos en el procedimiento
del aforo.
2
por usar el ADCP en modo estático, separando la sección en varias estaciones por
verticales, para así aplicar el cálculo de forma análoga a los equipos convencionales de
medición puntual. Sin embargo, durante avenidas, las condiciones del flujo dificultan el
mantener estático el instrumento, por lo que se considera importante avanzar en el
estudio de posibles metodologías alternativas que permitan obtener una estimación de
la incertidumbre de medida para la medición de caudal con ADCP en modo dinámico.
3
2 Objetivo general
El objetivo general de este trabajo, consiste en realizar una comparación de la calidad
de las mediciones de caudal con ADCP en modalidades dinámica y estática, enfocado
en las fuentes de incertidumbre de medida.
3 Objetivos específicos
Y como objetivos específicos, se persigue:
4
4. Marco teórico
A continuación se relacionan los conceptos teóricos utilizados para el desarrollo del
trabajo.
4.1 ADCP
&'(
𝑉"#$ = (1)
)'*
Donde,
𝑉"#$ es la velocidad relativa entre el agua con respecto al flotador,
𝐹- La frecuencia de la onda emitida,
𝐹" La frecuencia de la onda de retorno,
𝐶 La velocidad del sonido en el medio.
En el caso de aforos dinámicos, se mide la velocidad relativa del agua con respecto al
ADCP, y usando la función de rastreo de fondo, la velocidad del ADCP respecto al lecho,
con ambas velocidades se puede obtener la velocidad del agua respecto al lecho, la
velocidad es medida en celdas de dimensión constante y definida por el operario, dentro
de las limitaciones de diseño del instrumento. Resultando así el caudal medido por el
ADCP según la ecuación 2.
5
1
𝑄= 𝑉0 𝐴0 (2)
023
Donde,
𝑉0 es la velocidad del agua respecto al lecho para cada celda i; y 𝐴0 es el área de la celda.
4.3 Incertidumbre
“Parámetro, asociado con los resultados de una medición, que caracteriza la dispersión
de los valores que razonablemente podrían ser atribuidos al mensurando.
Nota: Las incertidumbres se expresan como un valor absoluto y no toman un signo
positivo o negativo” (UNE-ISO-5168, 2006).
Incertidumbre típica del resultado de una medición cuando ese resultado se obtiene a
partir de los valores de un número de otras magnitudes, igual a la raíz cuadrada positiva
de una suma de términos, que son las varianzas o covarianzas de estas otras
magnitudes ponderadas según como varía el resultado de la medición con cambios en
estas magnitudes.
Magnitud que define un intervalo alrededor del resultado de una medición que puede
esperarse abarque una gran fracción de la distribución de valores que podrían
razonablemente atribuirse al mensurando.
6
para obtener una incertidumbre expandida. Depende de la distribución de probabilidad
con la que se evalúa la muestra, suele emplearse la distribución normal en presencia de
errores aleatorios y con muestras representativas, cuando las muestras son pequeñas
se suele usar la distribución T de Student.
7
Figura 1 Representación gráfica de las componentes de errores de la incertidumbre
Una vez determinados todos los factores que introducen errores en la medición, las
incertidumbres típicas de las magnitudes de entrada y sus coeficientes de sensibilidad
asociados 𝑐0 , puede determinarse la incertidumbre total de la magnitud de salida de
acuerdo con la ecuación:
= )
𝑢: 𝑦 = 023 𝑐0 𝑢(𝑥0 ) (3)
?
𝑄 = 𝐹? 023 𝑏0 𝑑0 𝑣0 (4)
8
4.4.1 Incertidumbre debida a errores aleatorios en mediciones con molinete
) ) 3 ) ) )
𝐸C "
= 𝑒'? + 𝐸F " + 𝐸- " + 𝐸G " (5)
?
) ) ) )
𝐸GI " = 𝐸J + 𝐸: " + 𝐸# " (6)
"
En donde:
(Ep)r es la incertidumbre debida al número limitado de puntos de la vertical;
(Ec)r es la incertidumbre de la medida del molinete (se obtiene en la calibración del
equipo);
(Ee)r es la incertidumbre debido a las fluctuaciones del flujo.
) ) ) )
𝐸C K
= 𝐸F K + 𝐸- K + 𝐸: K (7)
En donde:
(Eb)s es el porcentaje de incertidumbre sistemática del instrumento que mide la anchura;
(Ed)s es el porcentaje de incertidumbre sistemática del instrumento que mide la
profundidad;
(Ec)s es el porcentaje de incertidumbre sistemática del tanque de calibración del
9
molinete.
) )
𝐸C = 𝐸C "
+ 𝐸C K
(8)
Las mediciones con ADCP en modalidad estática, guardan similitudes con la medición
de caudal con molinetes, la medición del ancho y el posicionamiento del equipo
separando la sección perpendicular al flujo en m verticales, se realiza con un
instrumento de medida de distancia, cinta métrica para secciones no muy anchas (<50m
por lo general), pudiéndose usar otros dispositivos para medidas cuando el ancho es
mayor, por ejemplo, teodolitos; la profundidad en cada vertical puede medirse bien por
el haz acústico del ADCP o con una barra de vadeo, cable con contra peso o similares.
Por lo tanto, y es lo que muchos hacen en la práctica, se puede usar la misma ecuación
para el cálculo de incertidumbre que en mediciones con molinetes, solo hay que
considerar algunos términos ya que el ADCP mide la velocidad con un principio
diferente:
La velocidad es medida en múltiples celdas en la vertical, teniendo solo áreas sin medir
cerca de la superficie del agua y cerca del fondo. La extrapolación en estas áreas suele
hacerse por el método potencial, que coincide con el perfil teórico en el que se
fundamentan los métodos de medición puntual contenidos en la norma ISO 748:2007.
10
En un estudio realizado (Lee k. et al.-2014) se halló que la desviación estándar de las
velocidades medidas por los 4 haces de un equipo StreamPro, tiende a reducirse al cabo
del tiempo. En la práctica, si se mantiene el equipo estable, este efecto puede
apreciarse, el software muestra la desviación estándar de la velocidad y del ángulo con
respecto al norte de la brújula. En caso de que la primera supere el 10% o el segundo
20º, el software advierte y recomienda extender la medición unos segundos más.
MN O ? 3 ?U3
𝑢C = 𝑢? ) +𝑢:L$ ) + 023 𝜎Q0
) + 2 ) )
023 𝜎Q0 𝜎Q0R3 𝑟0,TR3 (9)
CO CO
11
Si se tiene una variable funcional que depende de varias magnitudes medidas, en
nuestro caso el caudal.
Pueden distinguirse los dos tipos de incertidumbres: Tipo A y Tipo B. Las de Tipo A
pueden
Calcularse a partir de la desviación estándar de una muestra de mediciones
independientes entre sí.
3 1
𝑠 𝑋0 = ?23(𝑋0,? − 𝑋I )) (12)
(1U3)
K(`a )
𝑆 𝑋I = (13)
1
En general para una función dependiente de parámetros medidos, cada uno con sus
incertidumbres asociadas, siempre que sean independientes de otros parámetros puede
establecerse la siguiente relación para la obtención de la incertidumbre:
1 cC
𝛿𝑄 = 3 𝛿𝑋0 c` (14)
a
1 ? cC cC
023 T23 d` d` 𝛿𝑋0 𝛿𝑋T (15)
a e
12
Por último, para que la ley de propagación sea correcta, todos sus términos de
incertidumbres deben tener el mismo nivel de confianza o grado de cobertura.
1 cC )
𝛿𝑄 ) = 3 𝛿𝑋0 (16)
c`a
)
1 cC cC
𝛿𝑄 ) = 3 d` d` 𝛿𝑋0 𝛿𝑋T (17)
a e
13
5 Plan de Trabajo y Métodos
5.1 Equipos Utilizados
Para este trabajo se utilizó un ADCP marca RDI, modelo StreamPro, el cual posee las
siguientes especificaciones y limitaciones:
Especificación Descripción
Perfilador de Caudal:
Rastreo de fondo:
Transductor:
Frecuencia 2MHz
Configuración Janus de 4 haces con 20° de amplitud de haz
14
Figura 1 Sonda del ADCP Stream Pro (Simpson M., 2001)
5.1.2 Flotador
Se empleó los software WinRiver II, versión 2.20.0.0 (2016), para modalidad dinámica y
SxS PRO, versión 1.15 (2015), para modalidad estática, ambos desarrollados por
Teledyne RD Instruments.
15
Todos los ADCP poseen unas franjas de la sección perfilada en las que no es posible
medir o bien no se desea medir por considerarse que poseen perturbaciones por el
propio instrumento y su flotador. En esas franjas el Caudal se estima por extrapolación
de las velocidades medidas por el equipo.
Los programas informáticos que controlan las mediciones con ADCP realizan
estimaciones de la velocidad en estas zonas ciegas, a partir de las velocidades medias
de los conjuntos de celdas más próximos a estas zonas que den un valor aceptable. El
usuario puede elegir la forma del margen (cuadrado o triangular por lo general).
Se puede intuir, que en la mayor parte de los casos, el flujo que discurre en las cercanías
de los márgenes y el fondo es poco representativo para el caudal total, ya que las
velocidades son sustancialmente inferiores en estas zonas, no ocurre lo mismo para el
flujo en la zona en que está sumergido el transductor, en contacto con la superficie libre
del agua, en donde las velocidades suelen ser superiores a la media, debido a esto, en
16
la práctica, es recomendable que esta distancia sea la menor posible para asegurar que
el sensor no salga del agua durante los movimientos del flotador, y debe medirse, cuanto
se ha sumergido el transductor, con aceptable precisión. En las mediciones realizadas
para este trabajo, se calculó esta profundidad con la diferencia entre la profundidad total
medida utilizando barras de vadeo convencionales calibradas y la profundidad efectiva
que se obtiene con el ADCP, promediando en tres puntos diferentes de la sección, dos
cercanos a los márgenes y uno en el centro.
La velocidad del flujo medido por el ADCP es un vector por lo que tiene magnitud y
dirección. El equipo ADCP modelo StreamPro, mide la velocidad en la dirección de cada
17
uno de sus cuatro haces, para obtener un vector tridimensional a partir de estas
componentes de velocidad, es necesario utilizar 3 de las velocidades medidas en la
dirección de los haces, el valor medido por el cuarto haz es redundante y se usa para
aplicar una corrección a los errores asociados a movimientos de cabeceo y rotación del
flotador. Las ecuaciones para calcular las componentes de la velocidad relativa del flujo
con respecto al flotador son las siguientes (USGS, 2009):
𝐵3 − 𝐵) 𝐵n − 𝐵o 𝐵3 + 𝐵) + 𝐵o +𝐵n
𝑉f = ; 𝑉m = ; 𝑉p =
2 sin 𝜃 2 sin 𝜃 4 cos 𝜃
(18, (19) y (20).
El cuarto haz redundante sirve para calcular el error de velocidad según la siguiente
expresión y es utilizado para aplicar corrección por el software interno:
18
propone la necesidad de inferir más en este apartado en futuros trabajos, que con
suficientes mediciones e introduciendo variabilidad de cabeceo y rotación controlada,
orienten a ampliar el conocimiento del error producido por estos movimientos.
(1) Presencia de fondo móvil. Este trabajo limita el estudio a lechos estáticos.
Los sedimentos en forma de transporte de fondo dan errores en el trazado
de la trayectoria del ADCP, porque la consideración para el cálculo de la
velocidad del flotador con respecto al lecho supone que el lecho es estático,
existen formas de detectar el lecho móvil y aplicar correcciones, el software
WinRiver II incluye un test de lecho móvil pero eso introduce otras fuentes de
incertidumbre que complican aún más el proceso de determinación de la
incertidumbre de medida.
Los errores de medición del ángulo de la velocidad pueden ser inducidos por dos
fuentes (Lee et al, 2013):
19
entre 5% y 10% introduce una sobrestimación de caudal entre 0.38% y 1.31%
en el caudal total.
Durante las mediciones realizadas para este trabajo se tomaron fuertes medidas en
cuanto a calibración de la brújula (criterio de aceptación: inferior a 0.5 grados de error
en el test de calibración) y perpendicularidad de la trayectoria del instrumento con
respecto a la dirección del flujo.
&'(
𝑉} = (22)
)'*
20
Donde,
𝑉} es la velocidad del flotador,
𝐹- La frecuencia de la onda emitida,
𝐹" La frecuencia de la onda de retorno,
𝐶 La velocidad del sonido en el medio.
𝐷 = 𝑉} 𝑡 (23)
El StreamPro puede medir la profundidad usando el tiempo que tarda el sonido emitido
por sus transductores en reflejarse en el lecho y retornar a la fuente de origen, la
profundidad en cada posición será estimada a partir del promedio de los cuatro haces.
&•
𝑍= (24)
)
Donde,
𝐶 es la velocidad del sonido en el medio y t el tiempo entre el pulso emitido y el recibido,
el factor igual a 2 es porque el eco de la onda recorre dicha distancia dos veces, ida y
vuelta.
21
5.4. Ley de propagación de incertidumbres en mediciones con ADCP en modo dinámico.
Análisis de los factores determinantes en la estimación de incertidumbre en estas
mediciones
La resolución de este problema no es nada sencillo, las mediciones con ADCP contienen
muchos factores de incertidumbre que no solo dependen del equipo sino de las
condiciones del lugar donde se realiza la medida y que además son condiciones propias
del instante de tiempo en que se dan, además las incertidumbres de medida de las
diferentes variables que definen al caudal (velocidad del agua, profundidad y posición
en el ancho de la sección) no son independientes unas de otras, la velocidad del agua
se calcula referente a la velocidad del flotador, y de la velocidad del flotador, a su vez,
depende la coordenada de posición que al final se convierte en la medida del ancho.
Por lo que aparecerán términos del tipo de la ecuación (17).
La ecuación teórica general para calcular caudal a partir de mediciones por efecto
doppler, en modalidad dinámica, separando una sección en celdas de medidas de
velocidad (Simpson, 2001):
• -
𝑄= „ „
((𝑉‚ × 𝑉F ). 𝑘)𝑑𝑧𝑑𝑡 (25)
Donde Q es el caudal, 𝑉‚ 𝑦 𝑉F los vectores de las velocidades del agua y del flotador, 𝑘
el vector unitario en dirección perpendicular al plano de los vectores de velocidad (en
condiciones ideales es vertical), t el tiempo total de medición, d la profundidad en cada
vertical de celdas o conjunto (ensamble) de zumbidos (pings), dz es la diferencial de
profundidad y dt la diferencial de tiempo.
22
u1 Resolución espacial (el equipo asume la velocidad en una celda centrada entre 4
haces acústicos con ángulos divergentes, a medida que avanza en profundidad, más
diferencia hay entre el volumen medido y el asumido).
u3 ambigüedad de velocidad (el ADCP mide con ángulos desfasados entre diferentes
pulsos)
u6 velocidad del sonido (el ADCP asume la temperatura del agua constante e igual a la
de su sensor que está cercano a la superficie, asume además la salinidad contante que
es el principal parámetro que afecta a la velocidad del sonido, por lo general los ríos se
asumen con salinidad cero, el río Llobregat puede ser un caso especial por tener
formaciones de sal en su cuenca)
u8 velocidad del flotador (mucha velocidad del bote involucra mayor error)
u9 tiempo de muestreo (no es tan crítico para el caudal pero sí para la velocidad media)
u10 cerca del transductor (el periodo de espera del zumbido de retorno causa un tiempo
ciego que genera errores en las celdas superiores ( próximas al transductor)
23
u13 posicionamiento de las celdas (coloca la primera celda en que logra medir y el resto
sucesivas a esta según la longitud que se especifique, que no necesariamente coincide
con la real)
u15 tiempo (necesario para establecer la velocidad y la posición del flotador y debe
esperar que retorne antes de emitir un sonido nuevo)
u16 edge (distancia de la primera medición que puede realizar cerca de los márgenes,
esa área la extrapolará, con los errores que esto acarrea)
u18 modelo de caudal (modelos de velocidad área son usados en el algoritmo que
calcula)
u19 sumas finitas (como se hace en aforos con molinete, debe hallarse una relación
similar al número de verticales)
u20 elección del lugar y operatividad (un parámetro muy subjetivo, elegir lugares sin
convergencia o divergencia de flujo cercano, sin corrientes secundarias, con
condiciones estables, poner atención a la posibilidad de fondo móvil, de haberlo hay que
realizar test de fondo móvil y aplicar las correcciones pertinentes, aunque por lo general
en estos casos, la modalidad estática es recomendable ya que no usa rastreo de fondo
para el posicionamiento a través de la sección del flujo)
Se han propuesto varios métodos para estimar la incertidumbre del caudal medido con
ADCP en modo dinámico, todos suelen requerir un gran número de transectas medidos
y de condición de flujo permanente entre los transectas. Una solución teórica presentada
por (Huang, 2008), se basa en un gran número de suposiciones ideales, la sección
elegida debe ser perpendicular al flujo, la velocidad del flotador debe ser menor a la
velocidad media del agua y permanecer constante a lo largo del recorrido, no considera
fondo móvil, considera que la distribución de los caudales medidos en sucesivos
24
transectas se distribuyen de forma normal.
) ˆ‰
𝐸C = (26)
‡ C
25
6. Análisis de Resultados
En la tabla 1 se puede apreciar los datos correspondientes a la ubicación y
características generales de la estación donde se realizaron las mediciones.
Z= 490 Msnm
TIPO DE CURSO DE
AGUA: Río
MD: Hormigón; Roca y cantos de gran
MATERIAL MÁRGENES: MI: Roca MATERIAL LLIT / SOLERA: tamaño
ANCHO DE LA SECCIÓN: 20,20 m PROFUNDIDAD MEDIA: 1,43 m
LECTURA ESCALA (i): 1,06 m LECTURA ESCALA (f): 1,02 m
3 3
CAUDAL: 38,349 m /s ± m /s
TEMPERATURA MEDIA DEL
1,246 m/s
VELOCIDAD MEDIA DE LA SECCIÓN: AGUA 12,6 ºC
Tabla 1 carátula de un aforo con la información de la estación aforada y resumen de resultados del aforo
en modo dinámico del día 9 de mayo del 2019
El día 9 de mayo del 2019 se realizaron maniobras del embalse La Baells, ubicado cerca
de la ciudad de Berga, Cataluña, para generar una crecida con fines de mantenimiento
de la calidad ecológica del río Llobregat en sus tramos medio y bajo. El informe
proporcionado por el equipo de gestión de embalses de la Agencia Catalana del Agua,
donde se aprecian los hidrogramas esperados en distintas estaciones ubicadas aguas
abajo de la presa tras la apertura de las compuertas se muestra a continuación en la
figura 5.
26
70
60
50
40
30
Cabal (m3/s)
20
10
0
109
121
133
145
157
169
181
193
205
217
229
241
253
265
277
289
301
313
325
337
349
361
373
385
397
409
421
433
445
457
1
13
25
37
49
61
73
85
97
La estación EA066 Llobregat en Berga (Olvan), es una de las más cercanas a la presa,
el tiempo de tránsito esperado hasta este punto es de una hora, el hidrograma esperado
tiene forma escalonada ya que por su cercanía al embalse se esperaba un
comportamiento de respuesta similar a las aperturas escalonadas de las compuertas,
comportamiento que se suaviza a medida que el hidrograma transita por el río, este
comportamiento escalonado fue determinante para elegir la estación para realizar las
mediciones, ya que se esperaba contar con períodos de estabilidad (flujo casi
permanente) suficientemente largos para realizar los aforos, cada aforo se inició cuando
transcurridos unos 5 minutos, no se apreciaba cambios en el nivel en la escala,
adicionalmente durante el aforo se hacían sucesivas lecturas de la escala para
determinar si se producían cambios. La realidad no coincidió con el comportamiento
esperado, aunque sí hubo un cierto comportamiento escalonado, el nivel no permaneció
estable durante los tiempos como que se preveía en la planificación, presentando
cambios significativos durante los aforos de mayor duración. En los aforos en modalidad
estática, se decide reducir el tiempo de muestreo de velocidades en cada vertical a 40s
(inicialmente se propuso tener tiempos de muestreo de 60s), y reducir el número de
verticales al mínimo aceptable para agilizar el proceso de medición (aunque esto no
cumpla las recomendaciones de la normativa ISO 748,2007), teniendo un promedio de
duración para cada aforo en esta modalidad de 22 minutos, se realizan un total de 8
aforos, 4 en modalidad dinámica y 4 en modalidad estática, cuyos resultados se
exponen en la tabla 2.
27
Caudal
Hora Caudal por curva de Método de
Nivel (m) aforado
UTC gasto VISCAT (m3/s) aforo
(m3/s)
Tabla 2 Resumen de resultados de caudales medidos en las dos modalidades con ADCP
En esta tabla, se puede ver el resultado de los aforos y las lecturas de escala,
contrastados con el caudal calculado por curva de gasto, además se detalla la hora en
que fue realizado cada aforo y la modalidad en que fue realizado.
28
permanente).
Caudal por
curva de
Caudal Método de Desviación
Hora UTC Nivel (m) gasto
Medio (m3/s) aforo Standard
M 2/ 𝑀 EQ%
VISCAT
(m3/s)
9:15 0,9055 30,486 25,923 Estático N/A 11 N/A 7,5%
9:25 0,90875 30,734 30,77 Dinámico 4,36 3 1,15 14,3%
10:00 0,99275 38,596 37,004 Dinámico 3,55 3 1,15 11,1%
12:00 1,16215 55,723 52,748 Estático N/A 18 N/A 6,1%
12:20 1,21381 60,295 56,946 Estático N/A 20 N/A 5,6%
12:50 1,20914 59,86 60,345 Dinámico 0,36 4 1,00 36,3%
14:20 1,01502 40,962 38,921 Dinámico 0,10 6 0,82 7,8%
14:30 1,01223 40,687 38,349 Estático N/A 13 N/A 7,1%
Figura 5. Perfil de velocidades del ADCP en modo estático durante la crecida, con alto transporte de
sedimentos.
29
Figura 6 Perfil de velocidades para un aforo en modalidad dinámica, durante la recesión de la creciente,
cuando la concentración de sedimentos había disminuido
30
Figura 7 Perfil de velocidades por verticales del aforo con ADCP en modalidad estática.
Figura 8 Río Llobregat en la sección medida, a la izquierda durante la creciente con gran transporte de sedimentos,
a la derecha en la recesión, se aprecia a simple vista como reduce la turbidez.
31
7 Conclusiones
Desde que se han comercializado los equipos ADCP para medidas de caudal, se ha
intentado resolver el problema de estimar el grado de aceptación que tienen las
mediciones realizadas con estos equipos, las complejidades físicas del funcionamiento
del equipo han representado una serie de obstáculos por que complican la ley de
propagación de incertidumbre al utilizar los parámetros que realmente mide (frecuencias
e intensidades de ondas de sonido) para obtener, por medio de un algoritmo, lleno de
complicados cálculos y suposiciones, las diferentes variables que definen el caudal, la
velocidad del flujo y el área de la sección transversal al flujo.
Se ha realizado una investigación entre diferentes fuentes que han estudiado este
problema. De ellas se han recopilado y expuesto los factores que están involucrados en
la incertidumbre de medidas en modalidad dinámica. Comprendiendo bien los factores
pueden diseñarse nuevos experimentos en campo, buscando eliminar o disminuir
algunos, y al mismo tiempo, variar otros para intentar obtener tablas de calibración que
puedan dirigir hacía soluciones empíricas análogas a las tablas propuestas por la ISO
748:2007 para medidas de caudal con molinetes.
Basado en los resultados obtenidos en las mediciones en campo, se pudo apreciar que
en presencia de flujo con importante transporte de sedimentos y durante eventos de
crecidas, el método en modalidad estática puede ser la alternativa que dé resultados
más aceptables, siempre que no haya grandes cambios en el nivel durante la medición.
En el momento de la recesión, cuando el transporte de sedimentos disminuye, la
modalidad dinámica produce mejores resultados y es preferible ya que el tiempo de
medición se reduce, disminuyendo los errores asociados a la no permanencia del flujo.
32
ancho de sección, medidas de forma automatizada por el instrumento, generando
incertidumbres complejas de estimar y que requieren mayor profundidad en el estudio
de los factores involucrados, así como la ejecución de más ensayos de campo con el fin
de obtener resultados empíricos que logren mejorar el estudio de esta área.
33
8 Referencias Bibliográficas
González-Castro Juan A., Kevin Oberg, James J. Duncker (2000), Effect of Temporal
Resolution on the Accuracy of ADCP Measurements, 9p.
Hening Huang, R.D. Instruments, Inc., (2008), Estimating Precision of Moving Boat
ADCP Discharge Measurement, 12p.
Herschy, R. (2009). Streamflow Measurement (3rd ed., p. 507). Abingdon, Oxon: Taylor
& Francis.
Kyutae Lee, Hao-Che Ho, Muste Marian, Chun-Hung Wua, (2014), Uncertainty in open
channel discharge measurements acquired with StreamPro ADCP, 15p.
R.D. Instruments, Inc., (2015),. sxs pro software user’s guide, Vol. P/N 957-6260-00.
R.D. Instruments, Inc., (2011), Acoustic Doppler Current Profiler Principles of Operation
A Practical Primer, Vol. P/N 951-6069-00, 62p.
34
Doppler Current Profiler,121p.
Tokyay T., Constantinescu C., Gonzalez-Castro Juan A., (2009), Investigation of Two
Elemental Error Sources in Boat-Mounted Acoustic Doppler Current Profiler
Measurements by Large Eddy Simulations, 13p.
US Geological Survey, Water Supply Paper number 2175, online available at:
http://water.usgs.gov/pubs/wsp/wsp2175/, (1982).
35
ANEJOS
Ilustración A.1 Ficha del último aforo realizado ese día, puede apreciarse en la gráfica como caen los aforos
en la curva de gastos, se eligió solo uno de los aforos en modalidad dinámica para ser presentados, el resto
no dieron resultados satisfactorios.
Ilustración A.2 Cálculo de incertidumbre según ISO 748 automatizado en la Hoja de aforos.
36
Relació del nivell de l'escala - aforaments a EA066
80,0
70,0
60,0
50,0
m3/s
40,0
30,0
20,0
10,0
0,0
0,0 0,2 0,4 0,6 0,8 1,0 1,2 1,4
m
La creu vermella representa l'aforament del present informe. Les creus blaves són aforaments més antics. Les creus negres
corresponen als aforaments fets per OFITECO durant el present contracte.
Gráfico A.1 Curva de gastos de la estación con los aforos históricos y en negro los 8 aforos realizados el
9/05/2019, en los que se basó este trabajo.
37