Universidad Catlica de la Santsima Concepcin

Facultad de Ingeniera
Ingeniera Civil

Desarrollo de herramientas computacionales para el post-proceso

de datos registrados por un Perfilador de Corrientes Doppler
Acstico

Roberto Felipe Cheuque Aroca

Informe de proyecto de titulo para optar al ttulo de

Ingeniero Civil

Profesor gua:

Dr. Diego Caamao Avendao

Profesor Informante: Dr. Edwin Behrens Rincn

Concepcin, diciembre de 2012

Resumen
El presente proyecto tiene como fin desarrollar una herramienta computacional capaz de
post-procesar los datos obtenidos en terreno mediante equipos ADCP (Acoustic Doppler Current
Profiler). Los equipos ADCP presentan una gran ventaja al medir velocidades y caudales en ros,
son ms precisos, ms rpidos, y ms seguros para el operador, que los equipos convencionales.
Adems se reconoce que el correcto post-proceso de los datos recopilados entrega informacin
relevante para la estimacin del transporte de sedimentos en ros.
El objetivo del proyecto es implementar una herramienta computacional, a travs de la
suite de desarrollo MATLAB, que post-procese la informacin recopilada por el ADCP y
automticamente entregue los parmetros necesarios para la estimacin del transporte de
sedimentos de fondo y cuantifique la capacidad terica de transporte. El programa se disea sobre
la base de mdulos con el fin de desarrollar una suite de clculo abierta para la incorporacin de
futuros mdulos de clculo basados en los mismos datos. Cabe notar que las aplicaciones del
software sern verificadas con datos recopilados mediante ADCP en el ro Biobo.
Los resultados esperados corresponden a la plataforma informtica, que a travs de una
interface grfica ser capaz de clasificar parmetros hidrulicos y de transporte de slidos, as
como tambin de determinar la capacidad de transporte de fondo en un ro. Adems, y mediante
la herramienta desarrollada, se estimar la capacidad de transporte de sedimentos de fondo en el
ro Biobo.

Abstract
This project aims at developing a computational tool capable of post-process the field data
by ADCP equipment (Acoustic Doppler Current Profiler). The ADCP equipment have a huge
advantage when measuring velocities and flows in rivers, are more accurate, faster, and safer for
the operator than conventional equipment.
The project objective is to implement a computational tool, through the MATLAB
development suite, which post-process the data collected by the ADCP and automatically deliver
the necessary parameters for estimating the background sediment transport and quantify the
theoretical capacity of transport. The program is designed based on modules in order to develop a
calculation suite open for future incorporation of calculation modules based on the same data.
Worth noting, that software applications will be verified with ADCP data collected through the
Biobo River.
The expected results correspond to the computer platform, which through a graphical
interface will be able to classify and transport hydraulic parameters of solids, as well as
determining transmission capacity of a river bottom. Furthermore, the tool developed by be
estimated transport capacity of bottom sediments in the river Biobo.

ii

A mi madre Betty y a mi familia

iii

Agradecimientos
Al cuerpo acadmico y paradocente de la UCSC, en especial a los profesores de la carrera
por su invaluable apoyo durante toda mi estancia como estudiante.
Especial mencin merecen los Profesores Dr. Diego Caamao y Dr. Edwin Behrens por su
apoyo, ayuda y dedicacin en el desarrollo de este trabajo de titulo.
Al proyecto Fondecyt Sustainability of morphology features in alluvial rivers and coastal
systems in Chile: The Biobo river as study case y a INNOVA BIO BIO por el financiamiento del
presente trabajo.
Finalmente agradecer a mi familia y mi madre por toda la ayuda, apoyo y fuerza que me
han prestado durante toda mi vida.

iv

Tabla de contenidos
1

Introduccin ................................................................................................................................ 1
1.1

Motivacin .......................................................................................................................... 2

1.2

Marco .................................................................................................................................. 2

1.3

Objetivos ............................................................................................................................. 3

Mtodos ...................................................................................................................................... 4
2.1

Descripcin de la zona de estudio ....................................................................................... 4

2.1.1

Descripcin general de Ro Biobo............................................................................... 4

2.1.2

Clima ............................................................................................................................ 4

2.1.3

Suelos .......................................................................................................................... 5

2.2

Descripcin de un ADCP ...................................................................................................... 6

2.2.1

Descripcin general ..................................................................................................... 6

2.2.2

Limitaciones del ADCP ................................................................................................. 8

2.3

Obtencin de datos en terreno ........................................................................................... 9

2.3.1

Equipo de terreno ....................................................................................................... 9

2.4

Formato de datos .............................................................................................................. 10

2.5

Pre-procesamiento de datos ............................................................................................. 12

2.5.1

Lectura ....................................................................................................................... 12

2.5.2

Correccin ................................................................................................................. 18

2.6

Post proceso de datos ....................................................................................................... 20

2.6.1

Mdulo principal ....................................................................................................... 21

2.6.2

Mdulo de transporte de fondo (ecuaciones empricas) .......................................... 27

2.6.3

Mdulo de anlisis esttico ....................................................................................... 30

2.7

Formulaciones empricas para la estimacin del transporte de fondo ............................ 35

2.7.1

Mtodo de implementacin ..................................................................................... 35

2.7.2

Meyer-Peter and Muller (1948) ................................................................................ 35

2.7.3

Correcciones MPM .................................................................................................... 37

2.7.4

Shields (1936) ............................................................................................................ 38

2.7.5

Van Rijn (1984) .......................................................................................................... 40

2.7.6

Wilcock and Kenworthy (2001) ................................................................................. 42

2.8

Desarrollo de la interfaz grfica ........................................................................................ 44

v

2.8.1

Descripcin ................................................................................................................ 44

2.8.2

Variables globales...................................................................................................... 45

2.8.3

Interfaz primaria........................................................................................................ 46

2.8.4

Cuadrcula de carga ................................................................................................... 48

2.8.5

Interfaz secundaria .................................................................................................... 49

Resultados ................................................................................................................................. 52
3.1

Estimacin del esfuerzo de corte sobre el lecho por distintos mtodos .......................... 52

3.1.1

Ubicacin de las mediciones utilizadas en el clculo ................................................ 52

3.1.2

Perfiles utilizados para el clculo del esfuerzo de corte ........................................... 53

3.2

Estimacin del transporte de sedimento de fondo........................................................... 59

3.2.1

Estimacin y validacin de las ecuaciones empricas implementadas en el software

59

3.2.2
Comparacin del esfuerzo de corte sobre el lecho obtenido en el perfil transversal y
en las mediciones estticas. ...................................................................................................... 64
4

Discusin y conclusiones ........................................................................................................... 66

Bibliografa ................................................................................................................................ 68

Anexos ....................................................................................................................................... 70
6.1

Anexo 1, Tablas de ordenamiento de datos por ensamble. ............................................. 70

6.2

Anexo 2, Perfiles y datos utilizados para el clculo del esfuerzo de corte. ...................... 73

6.2.1

Correspondencia entre los archivos del ADCP y los perfiles medios obtenidos ....... 73

6.2.2

Perfiles estticos; Velocidad y error asociado........................................................... 74

6.2.3

Perfiles medios de velocidad ..................................................................................... 95

6.2.4

Perfiles de esfuerzo de corte calculados con el TKE. .............................................. 100

vi

ndice de figuras
Figura 2-1. rea de estudio de este proyecto de titulo. ..................................................................... 5
Figura 2-2. Mtodo de medicin del equipo....................................................................................... 6
Figura 2-3. Workhorse Ro Grande ADCP. ........................................................................................... 7
Figura 2-4. ADCP obteniendo datos en el ro Biobo, montado de manera lateral sobre una
embarcacin del tipo Zodiac. ...................................................................................................... 7
Figura 2-5. Perfil indicando las zonas ciegas de medicin del equipo (Simpson, 2001). .................... 8
Figura 2-6. Perfil transversal (transecto) tipo obtenido desde un ADCP, delimitando las zonas
ciegas del mismo. ........................................................................................................................ 8
Figura 2-7. Muestreador de sedimento de fondo tipo Helley Smith utilizado en este proyecto de
titulo. ........................................................................................................................................... 9
Figura 2-8. Descripcin grafica de un ensamble. ............................................................................. 10
Figura 2-9. Flujo de datos en la funcin lectura. ............................................................................... 13
Figura 2-10. Estructura de datos horizontal del programa. .............................................................. 15
Figura 2-11. Estructura vertical de datos del Fixed leader data (ejemplo). ...................................... 15
Figura 2-12. Flujo de datos resumido para la funcin Ensamble. .................................................. 17
Figura 2-13. Ejemplo del trabajo realizado por la funcin corte....................................................... 19
Figura 2-14. Operatoria bsica sobre las matrices realizada por la funcin Corte. .......................... 20
Figura 2-15. Estructura del programa contenedor de los mdulos. ................................................. 21
Figura 2-16. Correccin de la direccin por falta de datos. .............................................................. 23
Figura 2-17. Aplicacin de la media mvil sobre un vector .............................................................. 24
Figura 2-18. Ejemplo de aplicacin del mtodo de media mvil sobre una matriz. ......................... 24
Figura 2-19. Representacin grafica del clculo del ancho de celda o ensamble............................. 25
Figura 2-20. Zonas en las que se estimaran las velocidades ............................................................. 26
Figura 2-21. Representacin del clculo de las dos velocidades utilizadas por ensamble para el
clculo del caudal de fondo....................................................................................................... 26
Figura 2-22. Estimacin de velocidades superiores y laterales para el clculo del caudal. .............. 27
Figura 2-23. Sistema coordenado utilizado por el ADCP para la toma de datos. ............................. 28
Figura 2-24. Ajuste del MPM y su correccin sobre los datos registrados (Wong & Parker, 2006). 38
Figura 2-25. Curva de Shields (Vanoni, 2006). .................................................................................. 39
Figura 2-26. Relacin entre los esfuerzos de corte de referencia respecto al contenido de arena del
lecho (Wilcock & Kenworthy, 2001). ......................................................................................... 43
Figura 2-27. Bosquejo de la interfaz primaria del programa desarrollado. ...................................... 44
Figura 2-28. Bosquejo de las distintas interfaces secundarias del programa. .................................. 45
Figura 2-29. Ejemplo de las secciones graficas que manejan las variables globales descritas. ........ 46
Figura 2-30. Ejemplo grafico del funcionamiento de la grilla de carga. ............................................ 48
Figura 2-31. Plantilla del Mdulo principal. ...................................................................................... 49
Figura 2-32. Plantilla del Mdulo de ecuaciones empricas. ............................................................. 50
Figura 2-33. Plantilla del Mdulo de anlisis esttico. ...................................................................... 51
Figura 3-1. Datos obtenidos frente a la estacin de la DGA, sector Pedro de Valdivia. ................... 52
vii

Figura 3-2. Datos obtenidos en puntos cercanos al puente Juan Pablo Segundo. ........................... 52
Figura 3-3, Puntos estticos en el sector de la desembocadura. ...................................................... 53
Figura 3-4. Perfiles de velocidad media para los puntos 1 y 9 (de izquierda a derecha). ................. 54
Figura 3-5. Perfiles de variacin del esfuerzo de corte calculado a partir del TKE. .......................... 54
Figura 3-6. Datos del punto 25. ......................................................................................................... 55
Figura 3-7. Datos del punto 5. ........................................................................................................... 56
Figura 3-8. Histograma comparativo de los esfuerzos de corte calculados a partir del mtodo del
perfil promedio (logaritmo), mejor ajuste logartmico y TKE mximo (media = 9.75 N/m2). .. 58
Figura 3-9. Histograma de frecuencias incluyendo todos los mtodos de estimacin del esfuerzo de
corte (media = 6.7 N/m2) .......................................................................................................... 58
Figura 3-10. Ubicacin de los transectos XS1 y XS2 entre los puentes Juan Pablo Segundo y
LLacoln. .................................................................................................................................... 59
Figura 3-11. Puntos que componen al transecto XS1. ...................................................................... 59
Figura 3-12. Ubicacin especifica del transecto XS2. ........................................................................ 60
Figura 3-13. Perfiles de velocidad utilizados para la estimacin del transporte de sedimento en el
transecto XS1 (Eje x= Ensamble, Eje y=Profundidad (m), Velocidades en m/s)........................ 62
Figura 3-14. Perfiles de velocidad utilizados para la estimacin del transporte de sedimento en el
transecto XS2 (Eje x= Ensamble, Eje y=Profundidad (m), Velocidades en m/s)........................ 63
Figura 3-15. Mediciones que componen al transecto XS1. ............................................................... 64

viii

1 Introduccin
La medicin directa del transporte de sedimentos en los ros es un proceso largo,
engorroso y caro. Dependiente de una serie de condiciones ambientales e instrumentales que
adems lo hacen poco preciso en la estimacin de la cantidad real de sedimento que est
transportando un ro en un determinado momento. Es por ello que contar con herramientas
capaces de estimar el flujo de sedimentos de manera rpida y con una precisin aceptable a bajo
costo operacional es primordial para el desarrollo de actividades econmicas que utilizan los
sistemas fluviales.
La presente tesis aborda el problema de la estimacin del transporte de sedimento de
fondo a travs de su estimacin indirecta por medio de un Perfilador de Corrientes Doppler
Acstico (ADCP por sus siglas en ingles) y el desarrollo de herramientas computacionales capaces
de procesar los datos entregados por el equipo de manera rpida y fiable.
Dado que los datos de medicin tienden a ser consistentes pero no ordenados, el conjunto
de las herramientas computacionales desarrolladas debe ser capaz de leer los datos, corregirlos,
procesarlos y entregar resultados que se ajusten a la realidad.
Los equipos ADCP tienen una serie de ventajas sobre otros equipos para los mismos fines.
Entre estas ventajas se encuentran la gran cantidad de datos que es capaz de medir el equipo,
adems de la nula intervencin a la que es sometido el ro durante el proceso de medicin, sin
embargo, tambin cuenta con algunas desventajas que deben ser subsanadas antes del
procesamiento de datos. Entre stas encontramos las zonas ciegas del equipo en que no se tienen
datos y se realizan interpolaciones para rellenar la informacin faltante. Pero la ms relevante
hace mencin a las limitantes que presenta el software entregado por los fabricantes para realizar
estudios ms detallados como por ejemplo la estimacin del volumen de sedimento transportado.
Por lo tanto existe una tremenda oportunidad para potenciar estos equipos mediante el desarrollo
de herramientas que permitan post-procesar los datos.

1.1 Motivacin
La gran cantidad de estudios realizados sobre el transporte de sedimentos hace necesario
el desarrollo de herramientas que permitan poner en prctica estos conocimientos generados de
manera sencilla y amigable para su aplicacin en las distintas reas de la ingeniera.
Es por ello que el desarrollo de una suite de software altamente especializado que
permita potenciar el desarrollo terico con mediciones de alta calidad obtenidas a travs del ADCP
se hace tremendamente necesario.
La suite de software incluir distintas aplicaciones que se enfocaran segn su uso. En lo
que se refiere a este trabajo, el desarrollo se centrar en una plataforma de carcter modular
compuesta por un Mdulo de estimacin de transporte de fondo mediante formulaciones
empricas probadas cientficamente.

1.2 Marco
Este proyecto de ttulo se encuentra bajo el marco del proyecto Fondecyt denominado
Sustainability of morphology features in alluvial rivers and coastal systems in Chile: The Biobo
river as study case, del Dr. Diego Caamao A., el que tiene por objetivo general evaluar la
sostenibilidad de las caractersticas morfolgicas en el ro Biobo.
Para lograr este objetivo se requerirn tanto datos de terreno como una gran cantidad de
utilizacin de herramientas computacionales que permitan describir los procesos morfolgicos a
los que est sometido el ro.
Cabe destacar que el ro Biobo y su cuenca es uno de los principales focos de desarrollo
del pas el que durante los ltimos aos sufri grandes perturbaciones debido a la utilizacin de
sus aguas para fines productivos sumado adems al efecto que tuvo el reciente terremoto del ao
2010 sobre la cuenca del ro.
Este proyecto estima el cambio en el flujo de sedimentos y las distribuciones de llanuras
aluviales como lo son las zonas de inundaciones en la ribera del ro antes y despus de la entrada
en operacin de los proyectos hidroelctricos del Alto Biobo. Luego de este primer paso se
utilizarn los datos recabados como condiciones de contorno para implementar un modelo
morfodinmico que ser calibrado y validado a partir de distintos mtodos.

1.3 Objetivos
Los objetivos generales de este proyecto son:
- Generar algoritmos que permitan leer directamente los datos entregados por el
equipo obviando el programa del fabricante.
- Desarrollar una interfaz grfica de usuario que permita una fluida interaccin del
usuario con los datos registrados por el equipo y posteriores clculos que se realicen
con ellos.
Los objetivos especficos del proyecto son:
- Implementar ciclos de lectura directos, que permitan la importacin de los datos
registrados por el ADCP de manera directa sin pasar por el software propietario.
- Desarrollar funciones y scripts especficos para determinar algunas de las propiedades
hidrulicas del ro.
- Generar scripts que permitan hacer una estimacin del transporte de sedimento de
fondo en un perfil transversal de un ro.
- Implementar de manera modular programas que permitan un fcil acceso a los datos
calculados y registrados.

2 Mtodos
2.1 Descripcin de la zona de estudio
La zona de estudio del presente proyecto comprende la parte baja del Ro Biobo,
especficamente las zonas aledaas a la ciudad de Concepcin en donde se tomarn los datos y
muestras que permitan el clculo y la validacin de los resultados entregados por el software
desarrollado.

2.1.1

Descripcin general de Ro Biobo

La cuenca del ro Biobo cuenta con un rea de 24.260

, est ubicada entre los 36
45 y los 39 de latitud sur. El ro tiene su origen en los lagos Galletu e Icalma ubicados en la
provincia de Malleco, Cordillera de los Andes, novena regin. Su desembocadura est ubicada en
la comuna de San Pedro de la Paz, sobre el Golfo de Arauco, punto en el cual el caudal vara entre
los 250 y 2.000

(caudales medios mensuales).

En su avance el ro Biobo atraviesa tres geoformas caractersticas del pas, La Cordillera de

los Andes formada principalmente por rocas gneas producto de la gran actividad volcnica de la
zona, El valle central o depresin intermedia que se compone esencialmente de sedimentos
aluviales y fluvioglaciares junto a depsitos volcnicos; y por ltimo se encuentra la cordillera de la
costa constituida principalmente por rocas granticas y metamrficas con aportes de sedimentos
marinos terciarios y fluviomarinos del cuaternario, estos ltimos encontrados en la plataforma
costera adyacente (Cisternas, 1993).

2.1.2

Clima

El clima de la cuenca presenta caractersticas de transicin entre los climas Mediterrneo y

Templado hmedo, predominante en la zona sur.
En el sector de estudio (costero) se presenta un clima Templado Hmedo con
precipitaciones que fluctan entre los 1.200 y 2.000
anuales de norte a sur (MOP-DGA,
2004).

2.1.3

Suelos
Se distinguen tres sectores:
-

Seccin superior de la cuenca: Sobre los 3.500

los suelos son esquelticos
sobre afloramiento rocoso. Existe gran alteracin por el proceso de remocin de masa
y erosin.

Seccin media de la cuenca: Entre los 3500 y 800

la zona es
precordillerana con un valle hacia el poniente. Los materiales parentales dominantes
son las cenizas volcnicas. Se distinguen 3 sectores: en el norte suelos con
caractersticas de trumaos, de alta fertilidad y regados. En el centro se encuentran
extensos arenales de baja fertilidad y con problemas de drenaje y en el sur, suelos
limosos bien drenados y de alta fertilidad.

Seccin baja de la cuenca: Comprende la Cordillera de la Costa y desembocadura. Son

suelos generados a expensas del material formativos de la cordillera de la costa, hoy
intensamente erosionados por la accin antrpica de cultivos de trigo. Los suelos de
las Planicies Costeras de Arauco, poseen erosin leve y estn sometidos a una fuerte
presin por el sector forestal.

En la desembocadura mostrada en la Figura 2-1, el material que compone los suelos es

arena proveniente del sistema del Antuco, cuya permeabilidad es alta, escorrenta
lenta y baja capacidad de erosin fluvial (MOP-DGA, 2004).

Figura 2-1. rea de estudio de este proyecto de titulo.

2.2 Descripcin de un ADCP

2.2.1

Descripcin general

Los ADCP (Acoustic Doppler Current profiler) son equipos desarrollados en la dcada de
los 80. Su tecnologa y mtodos de medicin se basan en el efecto Doppler. Los ADCP aprovechan
este efecto enviando una seal de sonido con una frecuencia y amplitud normalizada hacia el
fondo del cauce, esta seal rebota sobre los slidos en suspensin que tericamente tienen la
misma velocidad del medio fluido en que estn inmersos, producto que las partculas tienen
velocidad, la onda que podemos llamar eco tiene una frecuencia distinta a la seal de sonido
normalizada. Son estas dos seales las que el equipo compara para estimar la velocidad del fluido
producto de la existencia del efecto Doppler. Se muestra un ejemplo del proceso en la Figura 2-2.

Figura 2-2. Mtodo de medicin del equipo.

Adems de la velocidad, el equipo es capaz de capturar otras variables del entorno como
la temperatura superficial del agua, la velocidad del sonido de la onda que se propaga, la salinidad,
presin y eco son algunas de ellas. Especial nfasis tiene la amplitud del eco que recibe el
instrumento, ya que sta se puede asociar directamente a la concentracin de sedimento
producto del principio de funcionamiento del equipo explicado anteriormente.
El equipo a utilizar en el desarrollo de esta tesis es el Workhorse Ro Grande ADCP
fabricado por RDI-Teledyne mostrado en la Figura 2-3. Este equipo est diseado para montarlo
sobre una embarcacin, con la cual se pueden trazar transectos1 sobre el ro para as obtener la
descarga del mismo. Este equipo cuenta adems con la capacidad de capturar numerosas
variables del entorno como las que se mencionan en el prrafo anterior, en el caso particular de
esta tesis, ste se utilizar para la obtencin de los transectos sobre el ro Biobo y de mediciones
estacionarias para determinar el movimiento del transporte de fondo del mismo ro.

Trayecto a travs del cual se realizan observaciones o toma de datos.

Figura 2-3. Workhorse Ro Grande ADCP.

Debido a la gran variabilidad de ambientes en los cuales pueden ser operados, estos
equipos cuentan con modos de funcionamiento, los que sirven esencialmente para configurar el
equipo de manera que ste adapte su funcionamiento interno para obtener mediciones de buena
calidad bajo diferentes condiciones ambientales (ver Figura 2-4). Adems cuenta con 4 beams2 por
lo que es posible cuantificar el error cometido en las mediciones de velocidad que realiza el equipo
(para el clculo de la velocidad son solamente necesarios 3 beams).

Figura 2-4. ADCP obteniendo datos en el ro Biobo, montado de manera lateral sobre una embarcacin del tipo
Zodiac.
2

Transductor cermico que permite la emisin y deteccin de seales de sonido.

2.2.2

Limitaciones del ADCP

Los inconvenientes que presentan este tipo de equipos estn principalmente asociados a
las reas ciegas de medicin, donde el equipo es incapaz de medir fsicamente lo que pasa en un
determinado punto del cauce o el equipo mide, pero esta medicin es incorrecta debido a la
interferencia causada por el ruido del side-lobe (onda de sonido generada en el borde del beam
del instrumento de manera parasitaria) y que imposibilita la medicin de datos en profundidades
cercanas al lecho del ro (ver Figura 2-5).
La distancia en la que el side-lobe acta es directamente proporcional del ngulo de
orientacin respecto a la horizontal que tiene el beam y de la profundidad del cauce en el punto
de medicin.

Figura 2-5. Perfil indicando las zonas ciegas de medicin del equipo (Simpson, 2001).

Dados estos inconvenientes, un ADCP puede obtener mediciones de perfiles transversales

similares a las mostradas en la Figura 2-6 donde quedan en evidencia las zonas donde es posible
registrar datos. Es por estos motivos que estos equipos actualmente cuentan con modos de
operacin, siendo cada uno de stos diseados para determinadas circunstancias.

.
Figura 2-6. Perfil transversal (transecto) tipo obtenido desde un ADCP, delimitando las zonas ciegas del mismo.

2.3 Obtencin de datos en terreno

El proceso de toma de datos en terreno comienza con una serie de consideraciones a
tener en cuenta para que el proceso mismo se desarrolle de manera ordenada.
-

2.3.1

Establecer una poltica para el manejo y almacenamiento de datos.

Verificar que el hardware y el software del ADCP funcionan correctamente.
Identificar otros equipos que pueden necesitarse en el proceso de toma de datos.
Comprobar que el equipo es capaz de registrar los datos bajo las condiciones que el
terreno impone.
Juntar y chequear los equipos externos acoplados al ADCP, verificar su correcto
funcionamiento junto a este.

Equipo de terreno
El equipo a utilizarse durante la toma de datos de validacin ser, adems del ADCP.
1) Equipo de posicionamiento global (GPS): Este equipo servir como referencia de los
puntos de toma de datos adems proveer la informacin necesaria sobre la
trayectoria real del ADCP durante el test de fondo que se realizar.
2) Muestreadores de sedimento: Se utilizaron muestreadores de sedimento del tipo
Helley Smith, de 6.4
, 15
de peso y con una malla de 250
(ver Figura
2-7).
3) Embarcacin: Bote a motor del tipo Zodiac, donde se montar el ADCP y que se
utilizara para construir los transectos sobre el ro Biobo.
4) Equipo de seguridad personal: Chalecos salvavidas, botas de pierna completa, zapatos
aptos para la tarea y dems equipamiento de seguridad menor.
5) Equipos secundarios utilizados por el ADCP: Computador porttil que permita recoger
en tiempo real los datos de descarga provenientes de las mediciones del ADCP.

Figura 2-7. Muestreador de sedimento de fondo tipo Helley Smith utilizado en este proyecto de titulo.

2.4 Formato de datos

Los datos registrados por el ADCP Ro Grande fabricado por Teledyne son almacenados
bajo la extensin .PD0. Este tipo de archivo tiene la capacidad de guardar datos en formato
binario y de manera secuencial segn las mediciones registradas.
La estructura bsica de los archivos .PD0 se compone de ensambles como muestra el
ejemplo grafico de la Figura 2-8, los que a su vez se subdividen en grupos, que corresponden a
identificadores, datos de configuracin del equipo (fijos) y datos fsicos medidos por el equipo
(variables).

Figura 2-8. Descripcin grafica de un ensamble.

Los ensambles se componen de sub- grupos de datos, segn el libro de comandos del
equipo (Teledyne, RD instruments, 2007) estos subgrupos de datos se definen en la Tabla 2-1.

Tabla 2-1. Estructuracin de un archivo .PD0 ordenado de manera secuencial por ensambles.
Ensamble

Tipos de datos

Header

2
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.

Header

Fixed
Variable
leader data leader Data
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
Fixed
Variable
leader data leader Data

Velocity
.
.
.
.
.
.
.
.
Velocity

Correlation
magnitude

Echo
Intensity

Percent
good

Bottom
track data

.
.
.
.
.
.
.
.
Correlation
magnitude

.
.
.
.
.
.
.
.
Echo
Intensity

.
.
.
.
.
.
.
.
Percent
good

.
.
.
.
.
.
.
.
Bottom
track data

Reserved

Checksum

.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.

Reserved

Checksum

10

La funcin de cada uno de estos subgrupos de datos se define como se indica a

continuacin:
Header: Corresponde al primer grupo de datos almacenado por el ensamble, entre los
datos que se registran aqu se tienen el identificador Hexadecimal del ensamble, el tamao del
mismo y la ubicacin relativa al inicio del ensamble de los dems grupos de datos.
Fixed leader data: Aqu se almacenan la variables que no cambian durante la realizacin
de la medicin, estas pueden ser, la configuracin del equipo, el firmware del equipo, en que se
empieza a medir, etc.
Variable leader data: Este conjunto de datos almacena parmetros fsicos que no varan
en el ensamble de medicin pero que si varan a lo largo de toda la medicin, entre estos datos se
encuentran la temperatura del agua, la velocidad del sonido en el agua, la salinidad, el nmero de
celdas de medicin, etc.
Velocity: Tal como su nombre lo indica en esta parte del ensamble quedan registradas las
velocidades absolutas de los sedimentos transportados por el flujo. Esta velocidad puede variar su
formato dependiendo de la configuracin del equipo que puede ser configurado para que registre
las velocidades respecto al norte magntico, respecto al bote al que esta adosado el ADCP o
respecto a los Beams del mismo equipo.
Correlation magnitude: Este coeficiente relaciona la intensidad del eco con sus respectivas
frecuencias, con ello se obtiene un valor normalizado representativo de cuan slido es el objetivo
sobre el cual rebot la onda. Su valor va desde 0 a 255 donde el primer valor representa la
ausencia de un medio slido y el segundo valor, representa que la totalidad de la celda medida es
slida.
Echo intensity: Intensidad del sonido emitido que rebota sobre las partculas y que capta
el ADCP. Este valor es normalizado y su escala va desde 0 a los 114.75 dB.
Percent good: Indica la calidad de la informacin recibida en cada celda de medicin, tiene
distintos formatos de almacenamiento dependiendo de la forma en que este configurada la
medicin de velocidades.
Bottom track data: Conjunto de datos que permiten calcular el posicionamiento del
instrumento en el instante de medicin, esto se realiza a travs del clculo de velocidades
respecto al fondo del ro y al registro de direccionamiento del equipo.
Reserved: Datos que coloca el fabricante del equipo y de los cuales no hay descripcin.

11

Checksum: Bits que se encargan de cerrar el ensamble y que verifiquen la integridad de

ste (se registra finalmente si el ensamble es apto para ser interpretado como un dato fsico real).

2.5 Pre-procesamiento de datos

La etapa de pre-procesamiento de datos se divide en las sub-etapas de Lectura y
Correccin de las mediciones.

2.5.1

Lectura
La lectura de datos se realiza siguiendo un esquema secuencial por ensamble, de manera
general el procedimiento es el siguiente:
1) Lectura del archivo de datos en forma de un vector binario (vbin).
2) Ubicacin del tamao del ensamble
3) Con el punto 2 se determinan las posiciones iniciales de cada ensamble relativas al
comienzo del vector vbin.
4) A partir de las posiciones iniciales, se crea un ciclo de lectura que permita extraer los
datos de cada ensamble, esto se logra creando matrices referidas a las posiciones de
cada nmero en formato binario que debe ser transformado a formato decimal.
5) Finalmente los datos de cada ensamble se organizan en matrices con el nombre del
grupo al que pertenecen, donde las filas de la matriz corresponden al ensamble, y las
columnas a la variable registrada por el equipo.
Las funciones o scripts desarrollados para este proceso son 3, los detalles de cada uno se
muestran a continuacin.

2.5.1.1 Lectura (a,b)

Esta funcin es la encargada en primer lugar de seguir la ruta del archivo, o conjunto de
archivos, a travs de las variables a y b que sern ledas, luego cada archivo es descompuesto en
un vector de valores binarios positivos para ello Matlab incorpora funciones que leen la raz
binaria de estos archivos, el vector creado como se mencion anteriormente se denomina vbin.
Este vector puede ser nico en el caso de leer slo una medicin o puede ser un conjunto de
vectores en el caso de un conjunto de mediciones, el tamao del vector es variable y depende del
tamao del archivo original que se est leyendo.
Esta funcin otorga los datos necesarios para una segunda funcin de lectura, denominada
Ensamble, su diagrama de flujo se muestra en la Figura 2-9.

12

[ne,ntd,offs,vld,fld,veld,comd,ecod,pgod,btrd,name,date,bytes]=Lectura(a,b)

No

Si
a=0

Dat=dir([b,a])

Dat=dir([d,\PD0]

Nd=largo(dat)

Para i=1 hasta Nd Hacer

Name(i)=Dat(i).name
Date(i)=Dat(i).date
Bytes=Dat(i).bytes

F=abrir([d,Dat(i).name],r)
Vbin(i)=fread(F,binario)
[ne(i),ntd(i),offs(i),vld(i),fld(i),veld(i),comd(i),ecod(i),pgod(i),btrd(i)]=ensamble(vbin(I))

FUNCIONES DE
CORRECCION
DE DATOS

Figura 2-9. Flujo de datos en la funcin lectura.

A travs de la funcin Lectura se obtienen las variables que sern utilizadas por toda la
suite del software, stas estn ordenadas por sus respectivos nombres y la descripcin de cada
una de ellas es:
ne: Matriz en la que se almacenan los tamaos de cada ensamble por archivo ledo,
variable solamente referencial obtenida del proceso de lectura de los ensambles
.
ntd: Matriz que contiene el nmero de tipos de datos registrados en cada ensamble por
archivo ledo, supone una referencia para saber qu datos se registraron en dicho
ensamble(Header data, Fixed data, Velocity data, etc.).

13

offs: Matriz de datos que contiene las ubicaciones relativas del inicio de cada tipo de datos
registrados en el ensamble, estn guardadas por archivo ledo.
vld: Matriz de datos en la que se registran los valores del conjunto Variable leader data
por archivo ledo, el orden en que se almacenan los datos es del tipo matricial, donde el ndice de
las filas representa el ensamble correspondiente, en tanto el ndice de las columnas representa el
dato que se est obteniendo referenciado segn las tablas del Anexo 1, Tablas de ordenamiento
de datos por ensamble.
fld: Matriz de datos en la que se registran los valores del Fixed leader data, tiene la misma
estructura que la Matriz vld, sin embargo, los ndices de las columnas representan datos distintos
definidos en el Anexo 1, Tablas de ordenamiento de datos por ensamble.
veld: Matriz en la que se registran las velocidades del flujo medidas por el instrumento, su
ordenamiento es similar a las matrices anteriores en cuanto a las filas que son las representantes
del ensamble que se est leyendo, sin embargo en el ordenamiento de las columnas difieren ya
que la estructura es distinta, incluyendo ahora las variables de celdas y Beams del equipo segn
la configuracin del mismo como se muestra en el Anexo 1, Tablas de ordenamiento de datos por
ensamble.
comd: Matriz donde se registran los datos del Correlation Magnitude data, para su
ordenamiento interno sigue la misma estructura de la matriz veld.
pgod: Matriz en la que se registran los datos del Percent Good data, su ordenamiento
interno sigue la misma estructura de la matriz veld.
ecod: Matriz en la que se registran los valores del Echo data, su estructura interna es
igual a la de la matriz veld.
btrd: Matriz en la que se guardan los registros del Bottom Track data, esta matriz tiene
una estructura interna similar a la matriz vld pero donde los ndices de las filas representan los
datos que muestra en el Anexo 1, Tablas de ordenamiento de datos por ensamble.
La estructura horizontal de organizacin de datos se describe en la Figura 2-10, la cual
separa los datos de cada archivo ledo segn su tipo y los almacena en cada uno de los Tipos de
datos.

14

Fixed leader
data
Variable leader
data
Velocity data

Tipo de datos

Percent good
data
Echo data
Correlation
magnitude data
Bottom track
data

Figura 2-10. Estructura de datos horizontal del programa.

Fixed leader data

Archivo 1

Archivo ...

Archivo M

Identificador
Ensamble 1 hasta N.

Identificador
Ensamble 1 hasta N.

Identificador
Ensamble 1 hasta N.

Cpu FW rev.

Cpu FW rev.

Cpu FW rev.

Ensamble 1 hasta N.

Ensamble 1 hasta N.

Ensamble 1 hasta N.

.......

.......

.......

Blank after transmit

Ensamble 1 hasta N.

Blank after transmit

Ensamble 1 hasta N.

Blank after transmit

Ensamble 1 hasta N.

.......

.......

.......

Beam angle
Ensamble 1 hasta N.

Beam angle
Ensamble 1 hasta N.

Beam angle
Ensamble 1 hasta N.

Figura 2-11. Estructura vertical de datos del Fixed leader data (ejemplo).

15

En tanto la estructura vertical de datos mostrada en la Figura 2-11, ordena cada una de las
matrices obtenidas anteriormente por archivo ledo, as la matriz del tipo Fixed leader data del
Archivo 1 queda almacenada en la primera posicin del grupo Fixed leader data. Cabe destacar
que la matriz almacena a su vez datos organizados en filas y columnas.
Para mayor informacin acerca de los datos mencionados en las tablas del Anexo 1, se
debe consultar el Workhorse Command and output data (Teledyne, RD instruments, 2007).

2.5.1.2 Ensamble (vbin)

Esta funcin es la encargada de la transformacin del archivo binario en variables
contenidas en matrices ordenadas, su proceso de funcionamiento es el siguiente (descrito en la
Figura 2-12).
1) En primer lugar la funcin extrae desde el archivo binario vbin los tamaos de cada
ensamble registrando adems el inicio de cada uno de ellos relativos al comienzo del
archivo vbin.
2) A partir de las posiciones iniciales de cada uno de los ensambles relativas al inicio del
archivo, se suman a ellas las posiciones de cada uno de los tipos de datos
guardndolos en cada ensamble, estas posiciones son entregadas en el Header data de
cada ensamble.
3) Con las posiciones determinadas de cada uno de los valores del ensamble, es posible
extraer el ensamble complemente en formato binario desde el vector vbin, la
conversin de estos datos est a cargo de la funcin Bitodec (binario).
4) Finalmente los datos son ordenados segn su tipo, asignados a las matrices
anteriormente descritas, as por ejemplo los datos de velocidad sern guardados en la
matriz veld, siendo separados los datos segn el ensamble y archivo al que
corresponden.

16

Inicio

Lectura Archivo
PD0

Entrega de las
posiciones relativas
del los grupos de
datos
Of=E+6+n de datos

Transformacin del
archivo PD0 a
binario
Vbin

i=1
E(i)=4

Datos de ejemplo
Fixed Leader Data
Fld=[]
J=0

E(i+1)=4+Vbin(E(i))

J=0 hasta largo(E)

Si
E(i+1)<Largo(Vbi
n)

Fld(j)=vbin(j)+Of(j)

Fin
E=E-4

Figura 2-12. Flujo de datos resumido para la funcin Ensamble.

2.5.1.3 Bitodec (binario)

La funcin bitodec se encarga de la transformacin de los datos desde el formato binario
en que estn guardados al formato decimal de fcil lectura, cabe mencionar que algunos datos
deben ser pasados del binario al hexadecimal, sin embargo para el alcance de esta tesis todos los
valores utilizados estn en formato decimal.
Los datos sobre los que se realiza la operacin aprovecha las ventajas de programacin
que tiene Matlab, as la formula de transformacin desde nmeros binarios a decimales queda
establecida por la ecuacin 2.1.

2.1

Donde:

Numero de columnas de
Matriz de nmeros binarios donde el subndice indica la fila en la que est el
nmero a convertir.
Matriz de nmeros decimales obtenidos a partir de la conversin.
17

La formula anterior demostr ser la ms eficiente para convertir una gran cantidad de
datos binarios a decimales, esta se desarroll con el fin de evitar la utilizacin de las formulas
incluidas en Matlab debido a su extrema lentitud en el proceso de conversin y al tiempo que
conlleva este proceso cuando se deben procesar una gran cantidad de datos binarios.

2.5.2

Correccin

La correccin de los datos se centra en la preparacin de stos para ser utilizados en los
siguientes procesos, existen tres grandes tems de correccin de datos.
1) Correccin de los nmeros negativos: Esta correccin se realiza debido a que los
nmeros negativos son almacenados por el equipo en formato inverso, esto quiere
decir que los nmeros negativos se guardan como un nmero positivo muy cercano a
la precisin dada por la profundidad de bits asignada por el equipo para este dato en
particular. La correccin se logra encontrando en primer lugar esos nmeros para
despus restarles la precisin con que fueron almacenados (generalmente 8 o 16 bits
de precisin).
2) Correccin de datos errneos registrados: En este caso el equipo asigna un nico
nmero entero para los datos que son registrados con errores, por ejemplo en la
medicin de la velocidad. La correccin se realiza identificando el valor entero
asignado por el equipo para este tipo de datos (generalmente de 16 bits de precisin)
y colocando en su lugar un cero. La eleccin de colocar un cero es para que no
intervenga en operaciones aritmticas posteriores entre distintas matrices de datos.
3) Correccin de datos por porcentajes bien medidos: Por ultimo este tipo de error est
asociado al equipo, este entrega una matriz llamada Percent good data donde queda
registrada la fiabilidad de los datos medidos, en el caso del software, cualquier dato
medido con una fiabilidad menor al 50% es eliminado y reemplazado con un cero
numrico.
4) Corte de los datos: Dado que el ADCP se configura antes de cada medicin con
parmetros preestablecidos, todos los ensambles tienen la misma cantidad de celdas
por tanto la profundidad mxima que mide el instrumento es constante, gracias a la
recoleccin de la profundidad real del ro registrados por el equipo y almacenados en
el conjunto de datos denominados Bottom track data, es posible cortar las matrices
de acuerdo al perfil de fondo del ro. Esto se realiza a travs de la creacin de una
matriz de profundidades de cada celda, esta matriz es comparada con los datos de
profundidad del entregados por el Bottom track data, recortando las profundidades
mayores al valor registrado por el instrumento (Figura 2-13).

18

Figura 2-13. Ejemplo del trabajo realizado por la funcin corte.

2.5.2.1 Corte (Matriz a recortar, btrd, fld)

Esta funcin requiere de los argumentos matriz a recortar, la que por ejemplo puede ser
la primera imagen de la Figura 2-13. Esta matriz ingresa en estado crudo, cada celda que la
compone se le asigna una determinada profundidad, esta profundidad es comparada con la
profundidad del lecho medida por el instrumento, si la profundidad de la celda es menor que la del
instrumento, la celda toma el valor de 1, si es mayor su valor es cero. Luego se multiplica esta
matriz consistente de 0 y 1 que llamaremos rectificadora por la matriz de velocidad, generando
esta multiplicacin termino a trmino de las matrices, la versin recortada de la matriz a
recortar, este proceso generalmente deja en blanco una gran cantidad de celdas debajo del perfil
del fondo medido por el instrumento. Se hace necesario entonces ajustar el tamao de la matriz
para que no ocupe recursos adicionales durante las operaciones que se realizan con ellas.
El corte de las filas sobrantes se lleva a cabo determinando la suma de los valores
ordenados en columnas de la matriz rectificadora, dado que esta matriz se compone de 0 y 1
adems que se sabe que la matriz por encima del perfil del fondo solo toma el valor de 1, la suma
de los valores ordenados en las columnas entrega la dimensin de las mismas, si se obtiene el
mximo valor de estos valores se obtendr la ubicacin de la fila ms profunda que tiene
mediciones validas, con este dato se cortan las filas sobrantes generando una matriz ms pequea
y liviana de operar.
El procedimiento expuesto anteriormente queda explicado de manera grafica en la Figura
2-14.

19

1
2
3
4
5
6
7
8
9
10
11

1
2
3
4
5
6
7
8
9
10
11

1
2
3
4
5
6
7
8
9
10
11

1
2
3
4

1
2
3
4

1
2
3
4

1
2
3
4

1
2
3
4

1
2
3
4

Matriz de
profundidades
8
9
10
11

8
9
10
11

8
9
10
11

8
9
10
11

8
9
10
11

8
9
10
11

1
2
3
4
5
6
7
8
9
10
11

1
2
3
4
5
6
7
8
9
10
11

1
2
3
4
5
6
7
8
9
10
11

Profundidad del
lecho

1
0
0
0
0
0
0
0
0
0
0

1
1
0
0
0
0
0
0
0
0
0

1
1
1
1
1
1
1
0
0
0
0

v1
0
0
0
0
0
0
0
0
0
0

v1
v2
0
0
0
0
0
0
0
0
0

v1
v2
v3
v4
v5
v6
v7
0
0
0
0

1
1
1
1

1
1
1
1

1
1
1
1

1
1
1
1

1
1
1
1

1
1
1
1

Matriz rectificadora
0
0
0
0

1
0
0
0

1
0
0
0

0
0
0
0

0
0
0
0

0
0
0
0

v1
v2
v3
v4

v1
v2
v3
v4

v1
v2
v3
v4

v1
v2
v3
v4

v1
v2
v3
v4

v1
v2
v3
v4

1
1
0
0
0
0
0
0
0
0
0

1
1
1
1
1
1
1
0
0
0
0

1
1
1
1

1
1
1
1

1
1
1
1

1
1
1
1

1
1
1
1

1
1
1
1

1
1
1
1
1
0
0
0
0
0
0

1
1
1
0
0
0
0
0
0
0
0

1
1
0
0
0
0
0
0
0
0
0

v1
v2
v3
v4
v5
0
0
0
0
0
0

v1
v2
v3
0
0
0
0
0
0
0
0

v1
v2
0
0
0
0
0
0
0
0
0

v1
v2
v3
v4
v5
Matriz de datos de
0
salida
0
0 v8 v8 0 0 0 0

v1
v2
v3
0
0
0
0
0

v1
v2
0
0
0
0
0
0

Matriz rectificadora
0
0
0
0

1
0
0
0

1
0
0
0

0
0
0
0

0
0
0
0

0
0
0
0

v1
v2
v3
v4

v1
v2
v3
v4

v1
v2
v3
v4

v1
v2
v3
v4

v1
v2
v3
v4

1
1
1
1
1
0
0
0
0
0
0

1
1
1
0
0
0
0
0
0
0
0

1
1
0
0
0
0
0
0
0
0
0

v1
0
0
0
0
0
0
0
0
0
0

v1
v2
0
0
0
0
0
0
0
0
0

v1
v2
v3
v4
v5
v6
v7
0
0
0
0

v1
v2
v3
v4

0
0
0
0

v8
0
0
0

v8
0
0
0

0
0
0
0

0
0
0
0

0
0
0
0

v1
v2
v3
v4
v5
0
0
0
0
0
0

v1
v2
v3
0
0
0
0
0
0
0
0

v1
v2
0
0
0
0
0
0
0
0
0

v1
0
0
0
0
0
0
0

v1
v2
0
0
0
0
0
0

v1
v2
v3
v4
v5
v6
v7
0

v1
v2
v3
v4

v1
v2
v3
v4

v1
v2
v3
v4

v1
v2
v3
v4

v1
v2
v3
v4

v1
v2
v3
v4

Matriz recortada

Fila

v1 v1 v1 v1 v1 v1 v1 v1 v1 v1 v1 v1
v2 v2 v2 v2 v2 v2 v2 v2 v2 v2 v2 v2
v3 v3 v3 v3 v3 v3 v3 v3 v3 v3 v3 v3
v4 v4 v4 v4 v4 v4 v4 v4 v4 v4 v4 v4
v5 v5 v5
v5 v5 v5
v6 v6 v6 Matriz a recortar v6 v6 v6
v7 v7 v7
v7 v7 v7
v8 v8 v8 v8 v8 v8 v8 v8 v8 v8 v8 v8
v9 v9 v9 v9 v9 v9 v9 v9 v9 v9 v9 v9
v10v10v10v10v10v10v10v10v10v10v10v10
v11v11v11v11v11v11v11v11v11v11v11v11

1
0
0
0
0
0
0
0
0
0
0

1
0
0
0
0
0
0
0
0
0
0

1
1
0
0
0
0
0
0
0
0
0

1
1
1
1
1
1
1
0
0
0
0

1
1
1
1

1
1
1
1

1
1
1
1

1
1
1
1

1
1
1
1

1
1
1
1

Matriz rectificadora
0
0
0
0

1
0
0
0

1
0
0
0

0
0
0
0

0
0
0
0

0
0
0
0

1
1
1
1
1
0
0
0
0
0
0

1
1
1
0
0
0
0
0
0
0
0

1
1
0
0
0
0
0
0
0
0
0

1
2
3
4
5
6
7
8

Matriz recortada
0
0
0
0

v8
0
0
0

v8
0
0
0

0
0
0
0

0
0
0
0

0
0
0
0

Suma 1 2 7 7 8 8 7 7 7 5 3 2
Maximo = 8

Figura 2-14. Operatoria bsica sobre las matrices realizada por la funcin Corte.

Las variables btrd y fld aportan la profundidad del lecho, el tamao de la celda, la
profundidad de la primera celda de datos (desde la superficie del agua) y la cantidad de celdas por
ensamble.
Las matrices que se recortan en el programa son las de velocidad, eco, percent good data y
la magnitud de correlacin, dado que son estos datos los que son registrados en forma de perfil
vertical (ensamble) sobre todo un transecto de datos.

2.6 Post proceso de datos

El post proceso de datos es donde los datos obtenidos del proceso anterior son tratados a
travs de distintos mtodos segn su uso final (visualizacin o clculo).
Dada la organizacin modular del programa mostrada en la Figura 2-15, se analizar a
travs de ella los scripts de programacin desarrollados para su funcionamiento.

20

Datos de navegacin

Mdulo Principal

Datos de Velocidad y
eco

Datos tabulados

ADCP

Seleleccin de
archivo
Mdulo de
transporte de fondo

Mtodo de
estimacin del
transporte de fondo

Mdulo de
transporte en
suspensin

Salida de variables
para su posterior
implementacin

Por formula
emprica
Por movimiento del
fondo con lecturas
estacionarias

Figura 2-15. Estructura del programa contenedor de los mdulos.

2.6.1

Mdulo principal

Este mdulo tiene como propsito la visualizacin y entrega de datos proporcionados por
el equipo. Estos datos son guardados y trabajados de manera tabular, mejorando as la velocidad
de todo el proceso.
Los scripts utilizados en este mdulo son los siguientes:

2.6.1.1 Magnitud (velocidad)

Esta funcin tiene por objetivo separar y calcular las componentes de los datos de
velocidad, eco, coeficiente de correlacin y fiabilidad de la medicin (percent good data). Funciona
aislando las componentes de cada uno de estos datos ordenados de igual manera sobre una
matriz (Tabla 6-3) segn su estructura de organizacin, adems para el caso de la velocidad ajusta
esta en caso de ser una medicin absoluta (velocidad de la embarcacin + velocidad de las
partculas) a slo la velocidad de las partculas.

2.6.1.2 Tracking (btrd,fld,vld)

Script que permite calcular la trayectoria de la embarcacin a partir de los datos
entregados por el ADCP, se utilizan variables como el ngulo de la brjula, las componentes de
velocidad del bote medidas respecto al lecho del ro y el tiempo entre cada ensamble, con esto es
posible el clculo de cada elemento diferencial de distancia que sumados permiten obtener la
trayectoria total respecto al norte magntico y en coordenadas N-E.
21

Formulas utilizadas para la transformacin de las coordenadas

2.2

2.3

Donde
Diferenciales de direccin norte y este respectivamente, se calculan por
ensamble, por ello son incrementos en la direccin respecto al punto anterior.
Tiempo de duracin de la medicin del ensamble, es obtenido como la diferencia
del tiempo acumulado registrado por el reloj del equipo, se utiliza para el clculo de la distancia
junto con la velocidad de la embarcacin.
Velocidades de la embarcacin registradas por el equipo respecto al
fondo del ro, estas estn ordenadas como la velocidad en direccin estribor-babor y proa-popa
respectivamente.
Producto de que el equipo registra las velocidades en milmetros, es necesaria su
conversin a metros, por ello se agrega a las formulas el divisor 1000, la precisin de los resultados
no se ve afectada dado que se est trabajando con una profundidad de datos en coma flotante
suficiente para evitar errores por truncamiento de datos.
El proceso anteriormente descrito utiliza datos corregidos para calcular la trayectoria,
dado que los datos ms vulnerables a errores son las velocidades respecto al fondo, estas deben
corregir a las otras variables utilizadas en el clculo, por ejemplo cuando una velocidad es
registrada de manera errnea en un ensamble, esta toma el valor de 0, automticamente el valor
correspondiente al ngulo de la brjula de ese ensamble toma el valor anterior, adems la
velocidad tambin es corregida y se le asigna el valor del ltimo dato correcto, Un ejemplo grfico
se muestra en la Figura 2-16.

22

Figura 2-16. Correccin de la direccin por falta de datos.

Finalmente para el clculo de la trayectoria total del barco solo hace falta hacer la suma
acumulativa de los incrementos en direccin norte y este calculados a travs de las ecuaciones 2.2
y 2.3 respectivamente.

2.6.1.3 Media mvil central

Este script se utiliza para la suavizacin de los datos errneos registrados por el equipo, se
divide en dos tipos de media mvil, en una dimensin (1d) para su aplicacin sobre vectores y en 2
dimensiones (2d) para su aplicacin sobre matrices, el algoritmo es parecido en los dos casos.
La media mvil utilizada corresponde al tipo central donde los datos a suavizar
corresponden al promedio de los datos circundantes al mismo.
La frmula utilizada para el caso de vectores (una dimensin)

2.4

23

Donde:
Valor medio calculado para el valor en la posicin .
Valor almacenado en el vector a suavizar en la posicin .
Valor indicador de cuantos valores
sern utilizados para el clculo de .
Nmero de valores vlidos para el clculo de la media (valores vlidos se definen
como todos aquellos distintos de cero.
Un ejemplo de aplicacin de la ecuacin 2.4 con un valor de
Figura 2-17.

es el mostrado en la

Valores a promediar

Vector promedio

2,5

2,5

2,5

Figura 2-17. Aplicacin de la media mvil sobre un vector

En el caso de las matrices, el procedimiento es similar, en este caso la ecuacin 2.4 cambia
su forma a la siguiente:
2.5

La definicin de las variables que componen la ecuacin es igual a la anterior, salvo la

diferencia de los trminos matriciales de posicin ( ) y los trminos de la cantidad de datos a
tomar para realizar el ajuste ( ).
Un ejemplo de la aplicacin de la ecuacin 2.5 utilizando valores de
frmula se muestra en la Figura 2-18.

Matriz a
promediar
Parte de la
matriz
promediada

1
5
0

0
6
4

de esta

4
4
3

6 .
.
.
3,85 .
.
.
4,25

Figura 2-18. Ejemplo de aplicacin del mtodo de media mvil sobre una matriz.

24

El mtodo de media mvil se utiliz en la realizacin de este programa para evitar el

ruido generado por la medicin de la turbulencia en las velocidades. La turbulencia no se puede
eliminar de los datos, sin embargo se puede suavizar a partir de este algoritmo. Se probaron
mtodos basados en Wavelets, sin embargo este tipo de mtodos es aplicable a ruido ms bien
pequeo dentro de una onda por lo que su aplicacin sobre la magnitud de la velocidad no
produjo mejoras importantes.

2.6.1.4 Clculo de caudal proyectado

Dado los datos de magnitud de velocidad, ancho del ensamble y alto del ensamble es
posible calcular la descarga tomada por la medicin, el dato de magnitud de la velocidad es
calculado de manera simple a travs de la funcin anterior llamada Magnitud, en tanto el dato de
la profundidad de la celda es entregado por el equipo, por lo tanto el nico dato faltante es el del
ancho de la celda.
El dato de ancho de celda es calculado realizando la proyeccin de los vectores (ver Figura
2-19) que forman los incrementos
y
con un vector unitario formado por dos puntos
elegidos sobre el transecto de medicin (por lo tanto su direccin es conocida).

Figura 2-19. Representacin grafica del clculo del ancho de celda o ensamble.

Luego cada una de las magnitudes de velocidad vlidas es multiplicada por el rea formada
por la altura de la celda y su respectivo ancho teniendo as el clculo del caudal por celda, si se
suman las celdas queda finalmente el caudal total que se midi.
Como se explic anteriormente el equipo cuenta con zonas ciegas producto de
limitaciones fsicas o tecnolgicas del mismo, para estimar el caudal de estas zonas es necesaria la
aplicacin de regresiones sobre los datos medidos para obtener con ellos los datos faltantes.
25

Figura 2-20. Zonas en las que se estimaran las velocidades

1) Estimacin del caudal de fondo (ver Figura 2-20, rea 1)

Para la estimacin de este caudal se ajusta sobre cada uno de los ensambles la regresin
logartmica de von Krmn-Prandtl, con la cual es posible estimar los parmetros del perfil
y
con los que se calculan las velocidades bajo el lmite que tiene el equipo para medir, lo anterior
queda representado grficamente en la Figura 2-21.
La distancia desde el fondo hasta la primera medicin vlida del equipo, es dividida en dos
celdas con rea y altura conocidas para tener una mejor estimacin del caudal que se transporta
en el fondo del ro.

Figura 2-21. Representacin del clculo de las dos velocidades utilizadas por ensamble para el clculo del caudal de
fondo.

2) Estimacin del caudal sobre la medicin (ver Figura 2-20, rea 2)

La estimacin del caudal sobre la medicin se realiza con las velocidades de la parte
superior del perfil logartmico ajustado, con esta velocidad conocida se estima el caudal de agua
cercano a la superficie que fluye por el ro como se muestra en la Figura 2-20 y en la Figura 2-22.

26

3) Estimacin de los caudales laterales

Una vez realizada la estimacin del caudal sobre la medicin, se procede a estimar los
caudales laterales del ro, para ello simplemente se ingresan las distancias desde la ribera del ro
hasta el punto donde el equipo comienza a medir, esta distancia es utilizada para crear un
tringulo en conjunto con las coordenadas del fondo, el rea del tringulo es de fcil clculo por
los tanto el dato geomtrico es conocido, la velocidad asignada a este tringulo corresponde a la
velocidad estimada ms cercana dividida en 2 para asegurar una distribucin pareja sobre el
tringulo de estimacin (ver Figura 2-22).

Figura 2-22. Estimacin de velocidades superiores y laterales para el clculo del caudal.

2.6.2

Mdulo de transporte de fondo (ecuaciones empricas)

Este mdulo tiene por objetivo estimar la cantidad de sedimento que lleva el fondo del ro
a travs de distintas formulaciones empricas. La base de trabajo son regresiones logartmicas
realizadas a los perfiles de velocidad entregados por el ADCP con los que es posible estimar
parmetros hidrulicos, tales como la velocidad de corte
y la rugosidad absoluta del fondo .
Este mdulo se construye sobre la base de los siguientes scripts:

2.6.2.1 Magnitud de velocidades

Script encargado de calcular la magnitud de las velocidades medidas por el equipo. En el
caso particular de que las coordenadas de orientacin y toma de datos del equipo estn
configuradas respecto al movimiento de la embarcacin como se muestra en la Figura 2-23, se
tiene la frmula generalizada 2.6.
27

2.6

Donde:
Magnitud de la velocidad medida por el equipo
Velocidad longitudinal registrada en el sentido de proa a popa de la embarcacin.
Velocidad trasversal registrada en el sentido de estribor a babor de la
embarcacin.
Velocidad vertical registrada por el equipo.
ndice correspondiente al ensamble o distancia que se est analizando.
ndice correspondiente a la profundidad o celda que se est analizando.

Figura 2-23. Sistema coordenado utilizado por el ADCP para la toma de datos.

Con los datos de la magnitud de la velocidad, se pueden calcular las velocidades de corte
generadas en el fondo, esto se realiza a travs de la formulacin de von Krmn-Prandtl que
asigna un perfil logartmico a la distribucin de velocidades en el eje vertical y que, relaciona la
profundidad con la velocidad y la velocidad de corte en el fondo.

28

2.6.2.2 Ajuste logartmico

Este script tiene por funcin proporcionar las variables de velocidad de corte y rugosidad
cerca del lecho. Funciona a partir de la relacin de von Krmn-Prandtl expresada en la ecuacin
2.7 donde el perfil de velocidades tiene una variacin logartmica con la profundidad, el clculo del
perfil se realiza ensamble por ensamble, de este ajuste estadstico se obtienen tres variables:
1) La velocidad de corte cerca del lecho
2) La rugosidad efectiva cerca del lecho
3) El coeficiente de correlacin lineal
A partir de los parmetros obtenidos es posible calcular el esfuerzo de corte sobre el lecho
del ro y as estimar el esfuerzo sobre las partculas que lo componen.
Adems, dado que los datos de velocidad se ven afectados por la turbulencia, es necesario
verificar la calidad del ajuste logartmico que se est realizando. Esto se lleva a cabo mediante la
clasificacin por el coeficiente de correlacin lineal obtenido en el ajuste y que indica la calidad del
mismo. Este parmetro estadstico vara entre -1 y 1 siendo los valores -1 y 1 las mejores
correlaciones de datos, valores cercanos a 0 indican una nula relacin lineal entre los datos. Su
formulacin es la planteada por la ecuacin 2.8.

2.7

Donde:
Velocidad transversal del perfil a una altura .
Velocidad de corte cerca del fondo del ro.
Profundidad a la que se est calculando la velocidad.
Profundidad a la que la velocidad se hace cero (Teledyne, RD instruments, 2007).
Producto de que el ajuste logartmico se debe realizar con datos de buena calidad, es
necesario eliminar los valores nulos propios de la medicin y que distorsionan los resultados del
perfil logartmico, adems se debe tener en cuenta que el perfil logartmico se desarrolla dentro
de una rango de profundidades, que generalmente van de 0 veces la profundidad hasta el 38% de
la profundidad del perfil, as se asegura una correcta estimacin de los parmetros anteriores.

2.8

29

Donde:
Coeficiente de correlacin lineal.
Covarianza de las variables aleatorias
Desviacin estndar de la variable .
Desviacin estndar de la variable .
2.6.3

Mdulo de anlisis esttico

Este Mdulo tiene como finalidad el clculo de la energa cinemtica turbulenta o TKE por
sus siglas en ingls. Este parmetro fsico es posible utilizarlo dado que para su clculo se
requieren las 3 componentes turbulentas de la velocidad. Las componentes turbulentas son
obtenidas a partir de las velocidades instantneas medidas sobre un punto dentro de cierto
perodo de tiempo.
Dado que el equipo entrega perfiles de velocidad instantnea, es posible calcular el TKE
para cada profundidad donde se ubica un grupo de celdas. Una vez determinado el TKE es posible
asociarlo al esfuerzo de corte sobre cada profundidad a la que fue calculado, con estos datos y si
se elige el punto ms cercano al lecho del ro se tiene una estimacin del esfuerzo de corte
actuando sobre ste. Finalmente este esfuerzo de corte puede ser comparado con los obtenidos a
partir de la ley logartmica de von Krmn-Prandtl para as realizar una calibracin entre los
esfuerzos obtenidos en un transecto y los puntos estticos posicionados en determinados puntos
del transecto.
Los scripts utilizados para el clculo de los parmetros anteriormente mencionados son:
2.6.3.1 Corte
Esta funcin implementada para el caso esttico tiene por objetivo truncar los datos
crudos de velocidad a una profundidad promedio. Su funcionamiento es muy parecido a la funcin
indicada en el punto 2.5.2.1. Se diferencian solamente en que esta nueva funcin de corte
promedia la profundidad a la que se truncar la matriz de velocidades dado que se asume que el
equipo est esttico registrando los datos y por lo tanto no hay variacin en la batimetra.
2.6.3.2 Turbulencia
La turbulencia es obtenida a partir de las velocidades instantneas medidas por el equipo,
la forma de obtenerla es relativamente sencilla, los pasos son:
1) Promediar respecto a una profundidad fija la componente de velocidad de la que se
requiere obtener la turbulencia.
2) A la velocidad instantnea medida restarle su propio promedio.
3) Repetir este procedimiento para cada una de las componentes y profundidades.
30

Con los pasos anteriores se obtienen las componentes turbulentas en cada direccin
separadas por velocidad y profundidad, con estos datos es posible calcular el TKE y a partir del TKE
se calcula el esfuerzo de corte promedio en la columna de agua y en el lecho.

2.6.3.3 Proyeccin de las componentes de velocidad

Dado que algunos esfuerzos de corte necesitan entre sus variables velocidades en la
direccin del flujo, es necesario definir transformaciones geomtricas de las velocidades obtenidas
con el equipo para poder obtener estas nuevas componentes.
La transformacin se realiza a partir de la matriz 2.9 que permite proyectar las magnitudes
de velocidad sobre la direccin deseada, los datos requeridos son:
-

Magnitudes de velocidad en el plano con direccin definida

Direccin de las proyecciones a calcular
La matriz de proyeccin utilizada se define como:
2.9

Donde
Velocidades referidas a la orientacin del equipo.
Velocidades referidas a la orientacin del transecto de proyeccin definido.
Angulo que hay entre las velocidades
y el perfil de proyeccin.

2.6.3.4 Esfuerzos de corte

Los esfuerzos de corte estimados en este mdulo tienen como base, el anlisis de las
componentes turbulentas de las velocidades, el desarrollo de dichas expresiones comienza por la
formulacin diferencial bsica que define el comportamiento del campo de velocidades en un
fluido.
Las ecuaciones que gobiernan el fenmeno escritas en forma tensorial son las siguientes:
Conservacin de la masa: ecuacin de continuidad

2.10

31

Conservacin del momentum: ecuaciones de Navier-Stokes

2.11

Conservacin de la energa trmica o de una sustancia disuelta en el fluido

2.12

Donde:
Componente de la velocidad instantnea en la direccin .
Presin esttica instantnea (no incluyendo la presin hidrosttica calculada con la
densidad de referencia ).
Unidad escalar que se puede reemplazar por la temperatura o la concentracin de
una sustancia.
Fuente volumtrica que puede expresar la generacin de calor (temperatura) o
reacciones qumicas o biolgicas (sustancia) en el fluido.
Viscosidad cinemtica del fluido.
Difusividad del fluido.
Aceleracin de gravedad en la direccin .

Las ecuaciones anteriores no se pueden resolver actualmente para fluidos turbulentos que
tengan relevancia prctica, por lo que se hace una aproximacin estadstica de las variables a
travs de las siguientes expresiones.
,

2.13

Donde
corresponden a los promedios de las variables vistas anteriormente
(velocidad, presin, temperatura o concentracin de alguna sustancia). Las variables
son las fluctuaciones de las variables respecto a sus valores medios, as se descomponen cada una
de las variables que conforman las ecuaciones 2.10, 2.11 y 2.12 en una componente representado
el promedio de la variable a la que se le suma la fluctuacin debido a la turbulencia. Esto se realiza
para poder sacar el trmino de la turbulencia de las variables.

32

Aplicando las ecuaciones 2.13 a las ecuaciones 2.10, 2.11 y 2.12 quedan reescritas de la
siguiente forma:
Ecuacin de continuidad

2.14

Conservacin del momentum

2.15

Conservacin de la energa trmica o de la concentracin de alguna sustancia

2.16

El proceso de promediar las variables introduce las correlaciones desconocidas

,
asociada a las fluctuaciones de velocidades y
, que es la fluctuacin de la cantidad escalar .
Fsicamente estas correlaciones multiplicadas por la densidad representan el transporte de
momentum y el transporte de calor o masa que estn descritas por las expresiones 2.17 y 2.18
respectivamente
2.17
2.18

La expresin 2.17, es conocida como el esfuerzo de corte de Reynolds y fsicamente es el

transporte de momentum sobre la direccin . Esta expresin sera la primera a utilizar para la
estimacin del esfuerzo de corte sobre el lecho del ro con la informacin recopilada por el ADCP.
2.19

33

Continuando con el desarrollo de las expresiones, (desde el punto de vista terico) se debe
encontrar un parmetro que relacione el trmino
con algn parmetro fsico del fluido (o
sobre el dominio) sobre el que se desarrollara la turbulencia, esta relacin se muestra en la
expresin 2.20.

2.20

La expresin 2.20 introduce el trmino de viscosidad turbulenta , que es una propiedad

de la turbulencia y no del fluido, este parmetro tiene una fuerte dependencia con el estado de
desarrollo de la turbulencia y puede variar de manera significativa desde un punto a otro.
El trmino corresponde a la energa cintica del movimiento turbulento, est definida en
trminos de la ecuacin 2.20, en donde se evalan los siguientes trminos vlidos producto del
operador de Kronecker .
As se llega a las siguientes expresiones

2.21

Reemplazando y despejando k desde la ecuacin 2.20.

2.22

La energa cintica turbulenta o TKE por sus siglas en ingls es utilizada para la estimacin
del esfuerzo de corte sobre el lecho del ro, la relacin se define como.
2.23

Donde la constante
toma el valor de 0.20 (Soulsby & Dyer, 1981) aunque el valor ms
extendido de esta constate sea 0.19 (Stapleton & Huntley, 1995) .

34

Una formulacin alternativa a la expresin es.

2.24

Las velocidades instantneas ortogonales ,

y
se les debe asignar una direccin
dado que al estar escritas en forma tensorial, estas no cuentan con direccin definida. Una forma
simple de hacer esto es aplicar un sistema de referencia sobre el cual
es la velocidad principal
de escurrimiento del fluido,
es la velocidad transversal al escurrimiento y
se puede definir
como la velocidad vertical al escurrimiento. Se debe tener en cuenta que se est asignado un
sistema coordenado a las velocidades, este sistema por lo tanto debe ser fijo y conocido.

2.7 Formulaciones empricas para la estimacin del transporte de fondo

2.7.1

Mtodo de implementacin

A partir de las velocidades de corte sobre el lecho y de la rugosidad total del lecho es
posible determinar los diferentes parmetros de las ecuaciones que sern utilizadas para el
transporte de sedimento, la aplicacin de las mismas se realiz por ensambles (perfiles verticales)
de velocidad.
Todas las ecuaciones tienen en comn la utilizacin de parmetros como el esfuerzo de
corte sobre el lecho obtenido a partir de la relacin descrita en la ecuacin 2.25.
2.25

En general las ecuaciones utilizadas entregan los resultados en trminos de volumen de

sedimento transportado, en cambio, en el programa desarrollado, el gasto de slidos es entregado
en unidades de masa, para as tener una mejor idea de la cantidad de material de fondo que se
est transportando en el cauce.

2.7.2

Meyer-Peter and Muller (1948)

El modelo desarrollado por estos autores es uno de los ms conocidos y aplicados en el

estudio del transporte de sedimento de fondo, este fue desarrollado a partir de experimentos de
laboratorio sobre canales con una ancho entre los 0.15 a 2
con pendientes que varan entre
los 0.0004 a 0.02 y profundidades del agua desde 1 a 120
.

35

Los sedimentos utilizados tenan una gravedad especifica ( ) de 1.25, los dimetros
medios oscilan entre los 0.4 a los 30

La ecuacin utilizada para el clculo de la descarga es la siguiente

2.26

Donde
es la tasa de transporte de sedimento adimensional, es el factor de forma del
lecho (correccin para eliminar el efecto de las formas del lecho) y es el esfuerzo de corte
adimensional. Cada uno de estos parmetros se calcula con las siguientes frmulas.

2.27

2.28

2.29

2.30

2.31

Donde
Transporte volumtrico de sedimento
Densidad relativa
Aceleracin de gravedad
Dimetro medio de las partculas
Profundidad del cauce
Rugosidad efectiva del lecho
Densidad del sedimento
Densidad del agua
obtenido de la ecuacin 2.25.

36

Esta frmula fue implementada por ensambles de esta manera se obtiene mayor certeza
sobre el clculo realizado ya que no se recurre a promedios de velocidad u otros parmetros que
distorsionen las mediciones de terreno.

2.7.3

Correcciones MPM

Distintos autores han propuesto correcciones a la ecuacin MPM, en este caso la

correccin que ser implementada en el programa es la propuesta por Wong & Parker, 2006.
Esta correccin surge del re-anlisis de los datos utilizados en la elaboracin de la ecuacin
MPM. La hiptesis planteada para el re-anlisis de los datos es que el factor de correccin del
esfuerzo de corte sub-estima el esfuerzo sobre las partculas en el caso de que el lecho sea plano y
no desarrolle rugosidades asociadas al cambio de geometra del mismo (formas del lecho como
dunas, antidunas, etc.).
Las ecuaciones presentadas utilizan los mismos parmetros de la formulacin original
(MPM), estas son:
2.32

2.33

Donde:
2.34

Esfuerzo de corte corregido sobre el lecho del ro (la correccin utilizada es la misma
que para el caso de MPM).
El esfuerzo de corte actuante sobre las partculas es el mismo que para el caso del mtodo
(MPM), dado que a partir de los datos del equipo es posible calcular la rugosidad efectiva del lecho
del ro (que incluye las formas de fondo) adems de la velocidad de corte (que incluye todos los
efectos de las rugosidades presentes en el lecho) no se hace necesaria la correccin sobre el
esfuerzo de corte por lecho plano dado que esto se realiza de manera automtica segn los
valores calculados.
Sin embargo cabe destacar que los dos ajustes anteriores no sobreestiman el transporte
de sedimentos tal como lo hace la formulacin (MPM) esto puede ser observado en la Figura 2-24.
37

Figura 2-24. Ajuste del MPM y su correccin sobre los datos registrados (Wong & Parker, 2006).

La ecuacin 2 de la Figura 2-24 corresponde a la formulacin original de MPM, en tanto la

ecuacin 22 corresponde a la primera correccin escrita anteriormente (ecuacin 2.32), los puntos
observables son los datos que se utilizaron para la estimacin de las dos relaciones.

2.7.4

Shields (1936)

La ecuacin de descarga de Shields fue desarrollada en 1936. Esta se obtuvo a partir de

experimentos realizados en dos canales de 40 y 80
de ancho respectivamente. La gravedad
especfica de los materiales utilizados como sedimento en los experimentos estaba en un rango de
1.06 a 4.2 y el dimetro medio de las partculas utilizadas vara desde los 1.56 a 2.5
.
La formulacin propuesta por Shields es:

2.35

Donde
Caudal del ro
Pendiente del lecho
Peso especfico del sedimento
Peso especfico del agua
Esfuerzo cortante sobre el lecho del ro
Esfuerzo crtico de corte
38

Esta ecuacin se basa en la utilizacin del parmetro que se interpreta como el esfuerzo
lmite sobre el que comienza el movimiento de las partculas sobre el lecho del ro, este parmetro
depende del fluido en que est inmerso el sedimento, del peso especfico del sedimento y del
tamao de las partculas. La relacin entre estas variables encontrada por Shields (ver ecuacin
2.36).

2.36

El tamao de partcula
se asume igual al
de la distribucin granulomtrica del
sedimento,
corresponde a la velocidad de corte crtica que puede ser escrita en trminos del
esfuerzo de corte crtico

, siendo

la densidad del agua. Por ltimo

es la viscosidad

cinemtica del agua.

El trmino del lado izquierdo de la ecuacin 2.36 se denomina esfuerzo de corte
adimensional, en tanto el trmino del lado derecho contenido en la funcin lleva por nombre
Nmero de Reynolds limite.
Shields y otros autores construyeron una curva mostrada en la Figura 2-25 partir de
experimentos desarrollados en laboratorio con la que es posible determinar la relacin entre las
dos partes que conforman la ecuacin 2.36.

Figura 2-25. Curva de Shields (Vanoni, 2006).

39

La curva de la Figura 2-25 muestra el esfuerzo de corte al que debe estar sometida una
partcula para iniciar el movimiento y por lo tanto el transporte, esta curva no tiene una ecuacin
explicita que pueda ser evaluada para encontrar por lo tanto lo que se hizo para encontrar este
parmetro esta explicado a continuacin:
1) Tomar distintos puntos de la curva de Shields.
2) Ajustar sobre estos puntos (Numero de Reynolds lmite) un cambio de variable a la
que una vez evaluada tuviera la forma de una curva polinmica.
3) Ajustar una regresin polinmica sobre la curva anterior para determinar los valores
de los coeficientes que acompaan al polinomio.
4) Escribir en trminos de la ecuacin 2.36.
5) Encontrar las races de la anterior a partir de soluciones numricas.
6) Esperar un resultado vlido de las iteraciones.
7) Utilizar el en las ecuaciones de transporte de sedimento que lo requieran.

2.7.5

Van Rijn (1984)

Van Rijn desarroll la formulacin a partir del trabajo de (Bagnold, 1954). Este se bas en
que el movimiento de las partculas de fondo era dominado por el fenmeno de saltacin de las
mismas y que las principales fuerzas actuantes sobre ellas eran las del fluido y la gravedad. El autor
trabaj sobre 130 canales distintos, con partculas que tenan dimetros medios entre los 200 a los
2000
, con profundidades del cauce mayores a 0.1
y nmeros de Froude menores que
0.9.
La ecuacin desarrollada tiene la siguiente forma.
2.37

La ecuacin anterior tiene un lmite de aplicabilidad para valores del parmetro

Para el resto de los casos la expresin utilizada es la siguiente:

2.38

Con los parmetros

calculados de la siguiente manera

2.39

40

2.40

Donde
Dimetro del porcentaje de partculas que pasa en la distribucin granulomtrica.
Aceleracin de gravedad
Viscosidad cinemtica del agua
Densidad del suelo relativa a la del agua (
)
Esfuerzo crtico obtenido desde la curva de Shields

El esfuerzo crtico efectivo

se calcula a partir del esfuerzo de corte del lecho
factor de correccin por la forma del lecho .

y el

2.41

En este caso

se define como

2.42

Donde

2.43

2.44

Profundidad del cauce

Rugosidad efectiva del lecho

41

2.7.6

Wilcock and Kenworthy (2001)

Esta formulacin para el clculo del transporte de sedimento tiene en consideracin el

efecto que tiene la arena sobre la estabilidad global del lecho de grava y su movimiento. A partir
de diversos experimentos los autores de esta formulacin encontraron una relacin directa entre
la cantidad de arena que estaba disponible sobre el lecho de grava y la cantidad de grava que era
transportada.
La ecuacin utilizada para el clculo es:

2.45

Donde
Tasa de transporte adimensional del material (arena o grava).
Peso especfico relativo (
).
Aceleracin de gravedad.
Descarga de slidos de la fraccin por unidad de ancho.
Fraccin del material .
Velocidad de corte al cubo.

La tasa de transporte adimensional se define como:

2.46

Donde
Relacin entre el esfuerzo de corte del lecho y el esfuerzo de corte de referencia
(
Parmetros de ajuste de la curva que describe la relacin entre la descarga de
sedimento y la relacin del esfuerzo de corte .
Los valores recomendados por los autores para los parmetros ,
1.27, respectivamente.

son 115, 0.923 y

42

El trmino
denomindado esfuerzo de corte de referencia, es calculado a partir de las
siguientes expresiones.
2.47

La ecuacin 2.47 es el esfuerzo de corte adimensional calculado por Shields, sin embargo
el parmetro
debe estimarse de otra manera, utilizando las siguientes frmulas propuestas en
el estudio.
2.48

Figura 2-26. Relacin entre los esfuerzos de corte de referencia respecto al contenido de arena del lecho (Wilcock &
Kenworthy, 2001).

La frmula anterior no es ms que un ajuste de los datos obtenidos por el autor mostrados
en la Figura 2-26. Esta frmula tiene lmites tanto para la fraccin de arena (sand o s) como la de
grava (g).

es el dimetro medio de la fraccin de arena (s) o grava (g),

del sedimento del lecho del ro.

es la fraccin de arena

43

Con estos datos es posible calcular la descarga total de sedimento sumando la descarga
producto de la fraccin de arena y la descarga de la grava como muestra la ecuacin 2.49.
2.49

La descarga estimada con este mtodo est en unidades de

, se debe multiplicar por el

ancho del cauce y el peso especfico del sedimento para obtener unidades de

2.8 Desarrollo de la interfaz grfica

2.8.1

Descripcin

La interfaz grfica de usuario se divide en mdulos de operacin, que estn clasificados

segn el archivo de datos que se est trabajando. Este proceso se realiza a travs de paneles
desplegables cuya cantidad vara segn el nmero de archivos que se estn consultando
simultneamente.
La interfaz primaria (ver Figura 2-27) se compone de los siguientes elementos:
1) Ventana de carga de archivos y mdulos
2) Ventana de seleccin de archivos
3) Men principal
La funcin principal de esta interfaz es servir de plataforma de las interfaces secundarias
que automticamente se generaran al llevar a cabo la lectura de los archivos de datos. Adems, en
esta plataforma se organiza el almacenamiento de los datos requeridos para el funcionamiento de
los diferentes mdulos de clculo.

Figura 2-27. Bosquejo de la interfaz primaria del programa desarrollado.

44

La interfaz secundaria corresponde a los archivos que se despliegan junto con el Mdulo
principal y que se encarga de visualizar los datos entregados por el instrumento (ver Figura 2-28).

Figura 2-28. Bosquejo de las distintas interfaces secundarias del programa.

Se compone de los siguientes mdulos para el post-proceso de datos.

1)
2)
3)
4)

2.8.2

Mdulo principal
Mdulo de ecuaciones empricas
Mdulo de anlisis esttico
Mdulo de Transporte en suspensin

Variables globales

El ADCP genera una gran cantidad de datos que deben organizarse y almacenarse en
diversas variables de modo que faciliten las operaciones y clculos posteriores. Se utiliza una
organizacin matricial que permiten reconocer el tipo de archivo y datos que se han almacenado
(ver Figura 2-29).
Son cinco las variables globales propuestas para el funcionamiento del programa:
Data: Es la primera variable global. Corresponde a una matriz de celdas en la que los datos
de los archivos ledos por la funcin lectura son ordenados en filas y columnas, las primeras se
utilizan para identificar el archivo y las segundas para reconocer el tipo de datos que est
guardado en cada posicin.
Tp: Corresponde a la variable que controla los paneles generados por archivo abierto, de
esta variable el dato ms importante que se puede obtener es el panel seleccionado actualmente,
a partir de este dato se pueden leer los datos requeridos por los mdulos desde la matriz data.
Bp: Controla y almacena los paneles incrustados en Tp, que corresponden a los mdulos
disponibles por archivos. Adems, una de sus funciones es servir de puntero de posicin, tal como
Tp pero solo de los mdulos activos visualizados en pantalla.

45

Mgh: Es una matriz de celdas similar a data, pero almacena los datos calculados en los
mdulos de cada archivo activo, sus filas son controladas por la seleccin de Tp y sus columnas por
Bp.
Hg: Matriz de estructura similar a Mgh pero que almacena el identificador para cada
modulo y archivo de cada interfaz grfica desplegada. Esta informacin se utiliza para activar y
desactivar la visualizacin de los objetos que no se estn utilizando. De esta manera prevalece una
estructura limpia y se libera memoria de la mquina.

Figura 2-29. Ejemplo de las secciones graficas que manejan las variables globales descritas.

2.8.3

Interfaz primaria

2.8.3.1 Ventana de carga de mdulos y archivos

Esta interfaz es la plataforma de todos los elementos grficos disponibles en el programa.
Sus despliegue grfico est acelerado bajo el estndar OpenGL que mejora significativamente la
velocidad de las visualizaciones de los elementos desplegados.
La interfaz es siempre ejecutada en el monitor primario ocupando aproximadamente un
80% de la superficie visible de la pantalla. Cabe notar que es posible maximizar y minimizar la
ventana de visualizacin haciendo ms cmodo su manejo cuando se trabaja con mltiples
aplicaciones al mismo tiempo.

2.8.3.2 Ventana de seleccin de archivos

La ventana de seleccin de archivos tiene por funcin proveer una lista de archivos
disponibles para su visualizacin y procesamiento. Adems, seleccionando alguno de los archivos
listados se despliega automticamente un mdulo asociado a ste. En el mdulo recin activado
es posible ver los datos registrados por el ADCP y los datos son asignados segn el orden con que
se han desplegado en pantalla archivos de datos anteriores, esto es:
46

a = Nmero de archivo activo creado (si es el primer archivo a desplegar, el numero

asignado es 1.
n = Nmero identificador del archivo desplegado en la lista visualizada.
La asignacin se realiza de la siguiente manera:
Data(a,:)=dat(n,:)
La variable Data crea una fila donde se guardan todos los datos contenidos en la variable
dat que proviene de la lectura del archivo y que est ordenada segn la lista desplegada.
Este procedimiento se realiza de manera de relacionar la posicin de los datos con el
identificador de las ventanas desplegadas, de esta manera es mucho ms fcil llamar los datos que
se necesitan, adems si se tienen archivos abiertos y se requiere cargar ms archivos, los datos
registrados por el equipo no se pierden dado que si estn desplegados en pantalla estarn
guardados de manera redundante en la variable Data.

2.8.3.3 Men principal

Este es el nico men disponible en el programa a travs de l, se accede al men de
archivos donde estn las opciones para la carga de los datos registrados por el ADCP en formato
PD0. Luego se tiene el men de exportar donde los datos registrados por el equipo es posible
exportarlos sin ningn tipo de tratamiento.
Otros tems del men principal son:
Mdulos: Men donde se selecciona el o los mdulos que se cargaran para un archivo
desplegado en pantalla.
Utilidades: Men que an no tiene una funcin definida pero que se le agregarn
funcionalidades en el futuro.
Ayuda: tem en el cual se puede acceder a la ayuda por contenidos (no disponible al
tiempo de la elaboracin del presente trabajo pero prevista en el futuro desarrollo de esta
herramienta) y un acceso al Acerca de, del programa donde se detalla la versin y los
desarrolladores que han trabajo sobre el programa.

47

Al estar compuesto de botones, el men principal llama funciones callbacks3 con la

finalidad de reconocer las acciones indicadas por el usuario. Estas funciones estn registradas
sobre el archivo .m del programa, y en ellas se registran las llamadas a las funciones que crean las
grillas de los elementos grficos.

2.8.4

Cuadrcula de carga

Cada mdulo tiene una funcin de creacin llamada cuadrcula (rejilla o grilla), su trabajo
es juntar las GUI4 independientes creadas para cada Mdulo y que cumplen una tarea especfica
para despus incrustarlas sobre los diferentes paneles creados para el despliegue de los datos.
Las cuadrculas de carga tienen una estructura definida y que debe ser seguida de manera
rigurosa si se quiere agregar un mdulo de operacin sobre el programa, su estructura de
funcionamiento es la siguiente y se explica de manera grfica en la Figura 2-30.
1)
2)
3)
4)
5)
6)

Declaracin de la funcin
Generacin del elemento grfico contenedor
Generacin de la cuadrcula de elementos
Carga de las distintas GUI sobre la cuadrcula de elementos
Visualizacin de los elementos cargados
Cierre de la funcin

Figura 2-30. Ejemplo grafico del funcionamiento de la grilla de carga.

Callbacks: Funciones que son llamadas a ejecutarse cuando se cumple cierta condicin sobre un elemento
grafico un ejemplo es hacer clic sobre un botn.
4
GUI: Graphics User Interface (Interfaz grafica de usuario), permite al usuario interactuar con los datos
procesados de manera grafica y simple, omitiendo las llamadas de datos a travs de cdigo y funciones.

48

2.8.5

Interfaz secundaria

2.8.5.1 Mdulo principal

El mdulo principal mostrado en la Figura 2-31, cuenta con tres GUI desarrolladas de
manera independiente que estn programadas para interactuar entre ellas. La primera de estas
GUI corresponde a un visualizador que muestra la trayectoria del ADCP en coordenadas N-E. Estas
coordenadas son obtenidas a travs de la funcin tracking y los datos entregados por esta funcin
son almacenados en la matriz Mgh(1,archivo)(1).
La segunda GUI es una grfica de superficie que indica 1) el perfil transversal de
velocidades, 2) eco, 3) el coeficiente de correlacin y 4) el porcentaje de mediciones con calidad
superior a la definida como lmite (i.e. mediciones buenas). Adems sobre esta grfica se trabaja el
algoritmo de correccin media mvil, que es activado al seleccionar Velocidades->Magnitud
corregida. Los datos procesados aqu son almacenados en la matriz Mgh(1,archivo)(3).
Por ltimo al costado izquierdo de la grfica se muestra el perfil de datos en la columna de
agua referido a un ensamble determinado. Este perfil se utiliza para observar la variacin de los
datos respecto a la profundidad de ro y tambin permite evaluar visualmente si los datos se
ajustan a un perfil logartmico.
En la parte inferior de cada una de las ventanas se encuentran cuadros de texto los que
tienen por funcin mostrar datos relevantes como la temperatura del agua, velocidades mximas,
coordenadas Norte-Este, etc.
Este mdulo es cargado automticamente al seleccionar una archivo desde la lista
desplegable. La carga del mdulo se realiza a travs de una grilla de carga donde el identificador
de esta grilla es almacenado en la matriz Hg, esto se realiza para hacer los llamados de datos una
vez cargada la grilla sobre la ventana principal.

Figura 2-31. Plantilla del Mdulo principal.

49

2.8.5.2 Mdulo de ecuaciones empricas (transporte de fondo)

El mdulo de ecuaciones empricas (ver Figura 2-32) es activado manualmente desde el
men mdulos. Aqu interactan cuatro GUIs, la primera es slo una ventana donde se
seleccionan las ecuaciones a utilizar en la estimacin del transporte de sedimento. La segunda
GUI, muestra los parmetros necesarios para las ecuaciones seleccionadas en color rojo y permite
ingresar estos datos en forma tabular.
La tercera GUI tiene como propsito seleccionar la distancia sobre el transecto sobre la
cual se calcular la capacidad de transporte de slidos de fondo. Esta ventana tambin ofrece la
visualizacin y las opciones de ajuste para 1) el perfil de velocidad medido, 2) la velocidad de corte
calculada a partir de la regresin logartmica, y 3) la rugosidad efectiva del fondo que tambin
proviene de la regresin logartmica.
La cuarta GUI cuenta con un grfico y una tabla de modo tal que se puedan comparar los
distintos resultados obtenidos a partir de las ecuaciones seleccionadas.

Figura 2-32. Plantilla del Mdulo de ecuaciones empricas.

2.8.5.3 Mdulo de anlisis esttico

Este mdulo (ver Figura 2-33) es activado manualmente desde el men de mdulos, est
compuesto por dos GUIs. La primera selecciona el ngulo respecto al cual se calcularn las
velocidades y muestra la velocidad calculada del fondo (distancia recorrida/ tiempo de medicin),
En tanto la segunda GUI est encargada de mostrar las mediciones de velocidad, las turbulencias
para cada una de las componentes de velocidad y calcular con ellas el TKE que es visualizado como
el esfuerzo de corte a distintas profundidades.

50

Figura 2-33. Plantilla del Mdulo de anlisis esttico.

51

3 Resultados
3.1 Estimacin del esfuerzo de corte sobre el lecho por distintos
mtodos
En esta seccin se estima el esfuerzo de corte sobre el lecho del ro a partir de los distintos
mtodos sealados anteriormente, con ello se obtiene un mximo de cinco valores de esfuerzos
de corte, los que se compararan para establecer cul es el valor que tiene mayor validez para cada
caso analizado.

3.1.1

Ubicacin de las mediciones utilizadas en el clculo

Figura 3-1. Datos obtenidos frente a la estacin de la DGA, sector Pedro de Valdivia.

Figura 3-2. Datos obtenidos en puntos cercanos al puente Juan Pablo Segundo.

52

Figura 3-3, Puntos estticos en el sector de la desembocadura.

No fueron ingresados todos los puntos mostrados en la Tabla 6-5 del Anexo 2, Perfiles y
datos utilizados para el clculo del esfuerzo de corte. Esto debido a que muchos puntos
registrados son dobles y tienen la misma ubicacin, solo difiere en la duracin del registro de
datos.

3.1.2

Perfiles utilizados para el clculo del esfuerzo de corte

Para el clculo del esfuerzo de corte a partir del ajuste logartmico es necesario el perfil
promedio de velocidades sobre un punto. Para esto el promedio de cada perfil fue realizado con el
software desarrollado, aplicando sobre los datos de velocidad el algoritmo de la media mvil
correspondiente a las matrices.
La densidad del agua se asumi como 1000

, para todos los clculos de los esfuerzos

aqu presentados.
Se promediaron las matrices de velocidad por fila para cada perodo de toma de datos (i.e.
1, 2 o 3 minutos segn sea el caso). De esta manera se obtienen los perfiles medios de velocidad.
Ms detalles de esta operacin se encuentran en el Anexo 2, Perfiles y datos utilizados para el
clculo del esfuerzo de corte. Cabe notar que las profundidades son consideradas en forma
creciente desde la superficie hacia el fondo (i.e. indicado por una lnea horizontal).

53

-0.2

-0.8

-0.4

-1

-0.6

Profundidad (m)

Profundidad (m)

-0.6

-1.2
-1.4

-0.8
-1
-1.2

-1.6

-1.4

-1.8

-1.6

-2
0

0.5

1.5

0.5

Velocidad (m/s)

Velocidad (m/s)

Figura 3-4. Perfiles de velocidad media para los puntos 1 y 9 (de izquierda a derecha).

Se incluyen tambin otros esfuerzos de corte obtenidos del anlisis de la turbulencia. Un

ejemplo de ello son los perfiles de esfuerzo de corte calculados a partir del TKE, los que estn
adjuntos en el Anexo 2, Perfiles y datos utilizados para el clculo del esfuerzo de corte.

-0.6

-0.8

Profundidad (m)

Profundidad (m)

-0.6

-1
-1.2
-1.4
-1.6

-0.8
-1
-1.2

10

11

Esfuerzo de corte (N/m2)

10

12

14
2

Esfuerzo de corte (N/m )

Figura 3-5. Perfiles de variacin del esfuerzo de corte calculado a partir del TKE.

Un resumen de los datos calculados es mostrado en la Tabla 3-1, en ella se registraron los
datos de esfuerzo de corte calculados a partir de diversos mtodos. Cabe mencionar que la
velocidad promedio en la zona estudiada del ro es cercana a 1

, esto analizando los perfiles

mostrados en la Figura 3-4 y en el Anexo 2 donde se adjuntan los perfiles de velocidad para los
dems puntos. El esfuerzo de corte presenta un distribucin espacial dispar ya que analizando
puntos cercanos (por ejemplo los puntos 5 y 25) se tiene una variacin desde los 8
40

a ms de

, este fenmeno de variacin se explica debido a que entre los puntos registrados existen

datos de velocidades con errores asociados bastante grandes, claro ejemplo de ello es el mismo
punto 25 (ver Figura 3-6), el que tiene una velocidad promedio cercana al 1
mximos de hasta 1.5

pero con errores

, lo que influye de manera drstica en el clculo de la turbulencia dado

que debido a este nivel de error, se incluye mas ruido sobre la medicin creando de manera
artificial turbulencia que finalmente es procesada y entrega valores de corte sobre el lecho que
54

son muy distintos a los calculados con un conjunto de datos con pocos errores (por ejemplo el
punto 5, Figura 3-7) .
a) Promedio de velocidades (en

Profundidad (m)

-0.5

-1

-1.5

-2

-2.5
0

0.5

Velocidad (m/s)

b) Registro de velocidades (en

Profundidad (m)

0
-0.5

-1

-1.5
-2
50

100

150

200

250

300

350

400

Ensamble

c) Error (en

Profundidad (m)

-0.5
-1

-1.5

-1

-2
50

100

150

200

250

300

350

400

Ensamble

Figura 3-6. Datos del punto 25. El punto a) representa el promedio de todas las velocidades instantneas obtenidas
durante la toma de datos para el punto, b) presenta el registro de velocidades instantneas y c) presenta el error
asociado a cada uno de los registros de las velocidades instantneas.

55

a) Promedio de velocidades (en

-0.5

Profundidad (m)

-1

-1.5

-2

-2.5
0

0.5

Velocidad (m/s)

b) Velocidad (en

Profundidad (m)

1.5

-0.5

1
-1

0.5

-1.5
-2
50

100

150

200

250

Ensamble

c) Error (en

Profundidad (m)

0.5

-0.5
-1

-1.5

-0.5

-2
50

100

150

200

250

Ensamble

Figura 3-7. Datos del punto 5. Tal como en la Figura 3-6, el punto a) representa el promedio de todas las velocidades
instantneas obtenidas durante la toma de datos para el punto, b) presenta el registro de velocidades instantneas y
c) presenta el error asociados a cada uno de los registros de las velocidades instantneas.

56

Tabla 3-1. Esfuerzos de corte sobre el lecho del ro calculado por distintos mtodos (en N/m2).

Punto

Fecha

Promedio
Velocidad
Logaritmo

Esfuerzo de corte
de Reynolds

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41

Mircoles 25-7
Mircoles 25-7
Mircoles 25-7
Mircoles 25-7
Mircoles 25-7
Mircoles 25-7
Mircoles 25-7
Mircoles 25-7
Jueves 26-7
Jueves 26-7
Jueves 26-7
Jueves 26-7
Jueves 26-7
Jueves 26-7
Jueves 26-7
Jueves 26-7
Jueves 26-7
Jueves 26-7
Jueves 26-7
Jueves 26-7
Viernes 27-7
Viernes 27-7
Viernes 27-7
Viernes 27-7
Viernes 27-7
Viernes 27-7
Viernes 27-7
Viernes 27-7
Viernes 27-7
Viernes 27-7
viernes 16-11
viernes 16-11
jueves 15-11
jueves 15-11
jueves 15-11
viernes 16-11
viernes 16-11
viernes 16-11
viernes 16-11
viernes 16-11
viernes 16-11

8,84
10,00
25,60
28,87
1,90
3,50
2,40
2,10
2,90
10,00
14,40
14,40
12,00
6,20
0,00
0,00
4,90
3,40
3,14
3,90
0,00
1,00
0,00
0,22
0,15
0,00
0,00
0,00
0,00
0,12
8,40
7,20
2,70
7,90
8,00
7,70
2,50
0,00
0,41
0,00
1,40

2,80
3,20
1,70
1,20
3,70
4,90
1,60
1,50
2,40
3,10
1,60
1,80
5,10
0,24
3,20
0,73
1,70
0,27
6,40
1,50
6,20
14,00
52,00
13,00
18,00
20,00
26,00
62,00
0,13
15,00
11,00
7,20
40,00
23,00
20,00
0,12
8,80
1,30
13,00
1,00
6,10

Esfuerzo de corte Aproximacin del

Esfuerzo calculado a
Esfuerzo de corte a
calculado a partir esfuerzo de corte partir del promedio de
partir del TKE en el
del TKE promedio con la componente
las velocidades de
lecho
(0,19*TKE*p)
vertical (0,9*p*w2) corte mejor ajustadas
5,10
7,20
5,60
5,50
8,30
8,70
4,80
6,20
9,90
10,00
5,10
8,80
6,10
6,60
5,80
5,20
8,10
6,40
5,30
5,10
60,00
56,00
60,00
56,00
43,00
48,00
100,00
87,00
75,00
78,00
4,40
6,10
51,00
47,00
52,00
9,00
33,00
6,80
3,80
8,20
2,60

1,50
1,80
0,62
0,52
3,20
3,50
1,70
1,90
2,70
3,30
1,90
2,40
1,90
1,60
2,20
2,30
2,90
2,20
5,70
2,10
8,80
19,00
44,00
24,00
19,00
16,00
59,00
49,00
33,00
28,00
1,70
0,87
25,00
13,00
9,20
1,20
2,20
1,20
1,20
2,20
1,80

11,45
12,32
26,80
28,73
13,36
13,69
10,53
10,00
29,24
31,83
12,32
11,88
16,33
15,93
10,00
7,06
7,40
9,22
6,39
8,46
88,51
99,48
220,90
106,15
193,60
48,40
144,40
159,76
136,90
126,31
10,63
9,01
100,30
62,50
87,08
9,70
3,14
7,94
5,63
15,50
7,74

5,40
6,00
3,00
3,00
11,00
11,60
6,50
7,20
16,00
16,00
6,00
9,60
5,50
6,30
8,00
7,00
6,70
7,20
6,90
5,30
80,00
62,00
170,00
75,00
60,00
90,00
170,00
135,00
96,00
125,00
6,80
6,10
68,00
58,00
58,00
7,00
33,00
7,20
4,00
9,00
4,50

Desde esta tabla se pueden hacer algunos anlisis con respecto a la relacin entre los
distintos esfuerzos de corte, en general las mejores relaciones se dan entre el esfuerzo de corte
calculado a partir de los perfiles de velocidad que mejor se ajustan a un perfil logartmico y el
esfuerzo calculado a partir de la energa cinemtica turbulenta. Teniendo en consideracin los
esfuerzos de corte logartmico y calculado a partir del TKE adems de obviando los puntos con
errores en el registro de velocidades demasiado excesivos se tiene un valor acotado con una
media de 9.75

(ver Figura 3-8).

57

Histograma de frecuencias
25

N de datos registrados

20

15

10

10

15
20
Esfuerzo de corte (N/m2)

25

30

35

Figura 3-8. Histograma comparativo de los esfuerzos de corte calculados a partir del mtodo del perfil promedio
(logaritmo), mejor ajuste logartmico y TKE mximo (media = 9.75 N/m2).

No se incluy el esfuerzo de Reynolds debido a que en general su valor siempre es ms

bajo que los otros esfuerzos calculados, por lo que creara un sesgo hacia la izquierda
promediando un valor menor al obtenido. En general y tal como se menciona anteriormente los
valores de esfuerzo de corte calculados se mantienen en rangos acotados, esto se debe
principalmente a que la velocidad media del ro es baja y casi constante para el sector analizado,
adems la configuracin del lecho bajo estas condiciones no produce altos niveles de turbulencia
haciendo el clculo del esfuerzo de corte menos certero por esta va.
Histograma de frecuencias
60

50

N de datos registrados

40

30

20

10

10

15
20
Esfuerzo de corte (N/m2)

25

30

35

Figura 3-9. Histograma de frecuencias incluyendo todos los mtodos de estimacin del esfuerzo de corte (media = 6.7
N/m2)

58

3.2 Estimacin del transporte de sedimento de fondo

Aqu se analiza el comportamiento predictivo del software para la estimacin del
transporte de sedimento en una porcin de un transecto del ro Biobo. Se comparan los
resultados de las distintas formulaciones empricas implementadas con resultados de campo
obtenidos a travs de la medicin directa del gasto de fondo.
Se procesan tambin los perfiles estticos realizados sobre el transecto. Con stos se
obtiene un patrn de comportamiento del esfuerzo de corte que acta sobre el lecho del ro.

3.2.1

Estimacin y validacin de las ecuaciones empricas implementadas en el

software

3.2.1.1 Ubicacin del transecto y puntos

Figura 3-10. Ubicacin de los transectos XS1 y XS2 entre los puentes Juan Pablo Segundo y LLacoln.

Figura 3-11. Puntos que componen al transecto XS1.

59

Figura 3-12. Ubicacin especifica del transecto XS2.

3.2.1.2 Comparacin de resultados

El sedimento transportado por el ro corresponde a arenas medias a gruesas, que son
caractersticas de la zona y se denominan arenas Biobo. Se reconocen siguientes propiedades
fsicas.

Tal como en el clculo de los esfuerzos de corte, se utiliz una densidad del agua con valor
de 1000

Los resultados de las mediciones con el Helley Smith se muestran en la Tabla 3-2, con
mediciones de tiempo constante iguales a 1 minuto y un factor del instrumento de 0.567 debido a
la relacin entre el rea de salida y el rea de entrada del mismo (Emmett, William W., 2010).

Tabla 3-2. Resultados de toma de muestras con el Helley Smith para el transecto XS1.

Identificador
pto1b
pto2b
pto3b
pto4b
pto0

Pto
42
36
38
40
0

Masa (gr)
1541,2
687,6
2543,5
571,3
343,8

Tiempo (s)
60
60
60
60
60

Distancia Wi (m)
4,4665
5,0155
4,7085
2,533
2,2405

Factor del instrumento

0,567
0,567
0,567
0,567
0,567

60

El transporte de sedimento de fondo estimado para el perfil XS1 a partir de los datos de la
Tabla 3-2 es:

Tabla 3-3. Resultados de toma de muestras con Helley Smith para el transecto XS2.

Identificador
p1
p2
p3
p4
p00

Masa (gr)
378,2
429,9
23,6
212,4
106,2

Tiempo (s)
60
60
60
60
60

Distancia Wi (m)
4,5
9
9
9
4,5

Factor del instrumento

0,567
0,567
0,567
0,567
0,567

El transporte de sedimento de fondo calculado en este caso es:

Cabe destacar que en el caso del segundo transecto (XS2) la distribucin de los datos
obtenidos con el Helley Smith se realiz de manera aproximada, ya que slo se contaba con la
posicin geogrfica del punto p1 por lo tanto estos datos no deben tomarse como definitivos ni
decisorios debido al error intrnseco en la distancia calculada para distribuir la masa capturada por
el Helley Smith.
Las expresiones que se utilizaron para el clculo del transporte de sedimento fueron 3, la
formulacin propuesta por Meyer-Peter and Muller, la formulacin de VanRijn y por ltimo las
correcciones realizadas por Parker y Wong a la formulacin MPM.
La estimacin del transporte de sedimento utilizando el ADCP es la siguiente:
Transecto XS1

Formulacin
MPM
Van Rijn
Correc. 1 MPM
Correc. 2 MPM
Promedio [Kg/s]

Material transportado [Kg/s]

Transecto 002
Transecto 004
Transecto 006
24,285
17,874
24,552
11,580
8,812
11,282
13,867
9,910
14,078
11,948
8,772
12,081
14,087
61

a) Transecto 002, suavizado del corte 0.06

2
1

0
-1

0.5

-2
-3
10

20

30

40

50

60

70

80

90

100

110

b) Transecto 004, suavizado del corte 0.08

3
2

1
0

0.5

-1
-2
-3
10

20

30

40

50

60

70

80

90

100

110

c) Transecto 006, suavizado del corte 0.08

0
-1

0.5

-2
-3
10

20

30

40

50

60

70

80

90

Figura 3-13. Perfiles de velocidad utilizados para la estimacin del transporte de sedimento en el transecto XS1 (Eje x=
Ensamble, Eje y=Profundidad (m), Velocidades en m/s).

Transecto XS2
Las estimaciones del transporte de sedimentos obtenidas a partir del ADCP son:

Formulacin
MPM
VanRijn
Correc. 1 MPM
Correc. 2 MPM
Promedio [Kg/s]

Material transportado [Kg/s]

Transecto 001
Transecto 002
Transecto 003
19,165
15,611
18,743
9,710
8,597
10,093
10,409
8,379
10,247
9,375
7,621
9,165
11,420

62

a) Transecto 001, suavizado del corte 0.1

1

1
0

0.5
-1
-2
20

40

60

80

100

120

b) Transecto 002, suavizado del corte 0.13

1

1
0

0.5

-1
-2
20

40

60

80

100

120

c) Transecto 003, suavizado del corte 0.15

1

1
0

0.5

-1
-2
20

40

60

80

100

120

140

Figura 3-14. Perfiles de velocidad utilizados para la estimacin del transporte de sedimento en el transecto XS2 (Eje x=
Ensamble, Eje y=Profundidad (m), Velocidades en m/s).

En general los resultados obtenidos con cada una de las formulaciones para estimar el
transporte de sedimento tiene un comportamiento similar en las seis mediciones realizadas para
el transecto XS1 y el transecto XS2.
En particular el transecto XS1 tiene el mejor nivel de prediccin sobre los datos medidos
por el Helley Smith, este transecto es el que tiene la mejor calidad de mediciones. Para el
transecto XS2 el comportamiento disminuye en precisin que se atribuye a la falta de puntos de
GPS con la ubicacin del Helley Smith. Cabe notar, que a pesar de ello, las formulaciones de Parker
and Wong presentan valores bastante cercanos a los reales, haciendo ver que el error est en la
falta de datos de coordenadas para estimar el punto exacto de extraccin de muestras.
El nivel de suavizado de corte se refiere al parmetro que permite calcular la media mvil
sobre un vector y que determina el largo del vector donde se calcular la media, este parmetro se
utiliza para rellenar y suavizar los datos faltantes de velocidad de corte y en el transecto de
medicin.

63

Para todos los clculos que involucren la utilizacin del mtodo del mejor perfil
logartmico ajustado, el coeficiente de correlacin lineal se ajusta a un valor de 0.8. Esto implica
que se asegura cierta calidad sobre los datos procesados ya que se debe cumplir que el ajuste
logartmico de velocidades tenga una coeficiente de correlacin R mayor a 0.8.

3.2.2

Comparacin del esfuerzo de corte sobre el lecho obtenido en el perfil

transversal y en las mediciones estticas.

Se valid solo el transecto XS1 debido a que cuenta con todos los datos georeferenciados.
Los tres transectos que los componen y los puntos estticos de medicin son mostrados en la
Figura 3-15.

Transectos XS1 referidos a un punto en la ribera

-10

-8

Este (m)
-4

-6

-2

2
5

Transecto 002

Transecto 004
Transecto 006
-5

Norte (m)

pto1b
pto2b
pto3b

-10

pto4b
Pto0

-15

-20

Figura 3-15. Mediciones que componen al transecto XS1.

Se calcul el esfuerzo de corte de los datos estticos a partir de los mtodos mencionados
en el captulo de mtodos, en tanto se eligi el transecto 006 con un nivel de suavizado de 0.08,
para estimar el esfuerzo de corte comparable al perfil por ser el ms cercano a los puntos
estticos.
Tabla 3-4. Comparacin del esfuerzo de corte calculado en las mediciones estticas y en el transecto total.

Archivo

Este (m)

Norte (m)

Seccin 1 008 sbt

Seccin 1 P2_0_000_SBT.PDO
Seccin 1 P2_0_002_SBT.PDO
Seccin 1 P2_0_004_SBT.PDO

670966,891
670972,023
670966,864
670964,885

5923720,957
5923728,759
5923725,821
5923716,304

Punto Identificador

42
36
38
40

pto1b
pto2b
pto3b
pto4b

Logaritmo
promedio
velocidades
0
7,7
0
0

Esfuerzo de corte [N/m^2 ]

Estacionario
Esfuerzo de
Esfuerzo de
corte a
Esfuerzo de
corte a
Mejor ajuste
partir de la
tke fondo
Reynolds
partir del
logartmico
turbulencia
TKE
vertical
0
0
0
0
0
0,1
9,0
1,2
9,7
7,0
1,3
6,8
1,2
7,9
7,2
1,0
8,2
2,2
15,5
9,0

Transecto 006
Mejor ajuste
logartmico
0
5,3
12,1
25,6

64

Los valores entregados en la Tabla 3-4 son consistentes, en general se aprecia que los 3
puntos de los que se tienen datos los valores de esfuerzo de corte (Por logaritmo y TKE) son
menores a los valores obtenidos por el transecto 006, sin embargo la diferencia entre los valores
no es exagerada y puede deberse al suavizado aplicado o a que el transecto no est precisamente
sobre los puntos estticos medidos lo que puede generar cierta variabilidad en los valores
obtenidos.

65

4 Discusin y conclusiones
Con las mediciones realizadas en terreno se pone a prueba la capacidad del software
diseado y se entregan resultados correspondientes a esfuerzos de corte promediados en el
espacio (i.e. transecto) y en el tiempo (i.e. estacionario). Futuras mediciones en terreno estarn
enfocadas en verificar el software con mediciones de velocidad de diversa magnitud, y bajo la
presencia de diferentes sustratos de fondo.
Los esfuerzos de corte en el sector bajo del ro Biobo son similares en magnitud, y se
hipotetiza se debe a que las velocidades registradas por el equipo son similares en todas las
mediciones realizadas, combinado con la ausencia de sedimento grueso en el lecho del ro
Las ecuaciones utilizadas en la estimacin del gasto de slidos de fondo fueron las de
Meller-Peter and Muller, VanRijn y las correcciones al MPM desarrolladas por Parker y Wong.
Estas tres formulaciones tienen condiciones de desarrollo que permiten aplicarlas al caso de la
arena Biobo, ya que por ejemplo la granulometra de la arena Biobo est en el rango de las
granulometras utilizadas para el desarrollo de las expresiones antes nombradas. Caso aparte es la
ecuacin desarrollada por Shields ya que a pesar de cumplir con los rangos de granulometra, se
omitieron los resultados de est por existir diferencias en los rdenes de magnitud de los valores
obtenidos. Una explicacin a este fenmeno es que a todas las expresiones evaluadas los autores
le incorporan un factor de correccin del esfuerzo de corte por las formas de fondo del ro (dunas,
risos, etc.), en cambio la formulacin de Shields aplica directamente el esfuerzo de corte total
calculado en las partculas del lecho, sobrestimando el transporte total de sedimento.
Futuros trabajos apuntarn a evaluar la formulacin de Wilcock and Kenworthy aplicable a
sedimentos gruesos (Grava) con fracciones de arena.
La tcnica utilizada para repartir el esfuerzo de corte a partir de los ajustes al perfil
logartmico de velocidades demostr ser bastante til y con una precisin aceptable. Sin embargo,
siempre es conveniente revisar los datos sobre los que se est trabajando e intentar estimar al
menos un valor de comparacin por otros medios o tcnicas como por ejemplo la estimacin por
turbulencia.
El mdulo de anlisis esttico se utiliz en la aplicacin de las mediciones estticas, donde
se analiza la turbulencia del ro. En este caso queda de manifiesto a partir de los perfiles de corte
calculados (a partir del TKE) que la turbulencia es generada en el lecho del ro, es all donde en la
gran mayora de los perfiles obtenidos alcanza su mximo, disminuyendo y tendiendo a un valor
constante al acercarse a la superficie del agua.

66

El diseo del software aqu presentado es de carcter experimental por lo que faltan
muchas horas de validacin de resultados. Adems, cabe destacar que por ejemplo el clculo del
gasto slido de fondo no siempre debiese ser igual al medido en terreno ya que a pesar de que
este parmetro es dependiente de la fuerza de arrastre del agua, tambin lo es de la cantidad de
sedimento disponible por lo que lo que se calcula finalmente es una tasa potencial de transporte
de sedimento, y no la tasa real.
Finalmente el software fue pensado desde un principio para que pudiera mejorar a partir
de la inclusin de nuevos mdulos de operacin, el proceso de creacin de los mdulos es
explicado en el Manual del desarrollador, en general es un proceso algo lento pero no complejo si
se tiene claro lo que se pretende hacer y los resultados esperables, adems las caractersticas de
funcionamiento modular generan una total independencia entre cada modulo de operacin lo que
conlleva una serie de ventajas desde el punto de vista del desarrollo de nuevas herramientas.

67

5 Bibliografa
Bagnold, R. A. (1954). Experiments on a Gravity-free Dispersion of Large Solid Spheres in a
Newtonian Fluid Under Shear. Proc. Royal Soc. , Vol225A.
Cisternas, M. (1993). Descripcion sedimentologica (granulometrica) de segmentos del curso
inferior y superior del rio Biobio. En EULA-Chile, Monografias cientificas Volumen 12. Santiago:
Grafica Andes.
Emmett, William W. (2010). U.S. geological Survey Scientific Investigations Report , 159-170.
Harvard center for geographic analysis. (s.f.). Chile Earthquake Data Portal. Obtenido de
http://cegrp.cga.harvard.edu/chile/?q=resources_data
Kim, S. C., Friedrichs, C. T., Maa, J. P., & Wright, L. D. (2000). Estimating Bottom Stress in Tydal
Boundary Layer From Acustic Doppler Velocimeter Data. Journal of Hydarulics Engineering , 399406.
Meyer-Peter, E., & Mller, R. Formulas for Bedload Transport. 2nd Congress IAHR. Stockholm.
MOP-DGA. (2004). Cuenca del rio Biobio. Santiago: Sinia.
R., M. D. (2009). Measuring discharge with acoustic Doppler profiler from a moving boat. Reston,
Virginia: U.S. Geological Survey.
Shields, A. (1936). Anwendung der Aenlichkeitsmechanik und der Turbulenzforshung auf die
Geschiebebewegung. Mitteilungen der Preussischen Versuchsanstalt fur Wasserbau und Shiffbau .
Simpson, M. R. (2001). Discharge Measurements Using a Broad-Band Acoustic Dopplet Current
Profiler. Sacramento, California: U.S. Geological Survey.
Soulsby, R. L., & Dyer, K. R. (1981). The formof the near-bed velocity profile in a tidally accelerating
flow. Journal Geophys , 8067-8074.
Stapleton, K. R., & Huntley, D. A. (1995). Seabed stress determination using inertial dissipation
method and turbulent kinetic energy method. Earth Surface Processes and Landforms 20 , 807815.
Teledyne, RD instruments. (2007). Workhorse, Command and output data format.
van Rijn, L. C. (1984). Sediment Transport, Part I: Bed Load Trasnport. Journal of Hydraulic
Engineering , 1431-1456.
Vanoni, V. A. (2006). Sedimentation Engineering. Reston, Virginia: ASCE Manual.

68

Wilcock, P. R., & Kenworthy, S. T. (2001). A two-fraction model for the transport of sand/gravel
mixtures. Water Resource Res. 38(10) , 1194.
Wolfang, R. (1984). Turbulence models and thei application in hydraulics. Rotterdam: IAHR.
Wong, M., & Parker, G. (2006). Re-analysis and correction of bedload relation of Meyer-Peter and
Mller using their own database. Journal of Hydraulics Engineering, Vol 132, No. 11 , 1159-1168.

69

6 Anexos
6.1 Anexo 1, Tablas de ordenamiento de datos por ensamble.
Tabla 6-1. Ordenamiento de datos del variable leader data (Posicin = ndice de la columna).

Posicion
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23

N de bits
16
16
8
8
8
8
8
8
8
8
16
16
16
16
16
16
16
16
8
8
8
8
8

Variable leader data

Nombre
Posicion N de bits
Variable leader id
24
8
Ensemble number
25
8
Rtc year (TS)
26
8
Rtc month (TS)
27
8
Rtc day (TS)
28
8
Rtc hour (TS)
29
8
Rtc minute (TS)
30
8
Rtc second (TS)
31
8
Rtc hundredths (TS)
32
8
Ensemble number msb
33
32
Bit Result
34
16
Speed of sound (EC)
35
32
Depth of Transducer (ED)
36
32
Heading (EH)
37
8
Pitch (EP)
38
8
Roll (EP)
39
8
Salinity (ES)
40
8
Temperature (ET)
41
8
Mpt minutes
42
8
Mpt seconds
43
8
Mpt hundredths
44
8
Hdg std dev
45
8
Pitch std dev

Nombre
Roll std dev
Adc channel 0
Adc channel 1
Adc channel 2
Adc channel 3
Adc channel 4
Adc channel 5
Adc channel 6
Adc channel 7
Error status word (ESW) (CY)
Spare
Pressure
Pressure sensir variance
Spare
Rtc century
Rtc year
Rtc month
Rtc day
Rtc hour
Rtc minute
Rtc second
Rtc hundredths

70

Tabla 6-2. Ordenamiento de las variables del Fixed leader data.

Posicion
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19

N de bits
16
8
8
16
8
8
8
8
16
16
16
8
8
8
8
16
8
8
8

Fixed leader data

Nombre
Posicion N de bits
Nombre
Fixed leader id
20
8
Coordinate transform (EX)
Cpu FW ver
21
16
Heading Alignment (EA)
Cpu FW rev
22
16
Heading bias (EB)
Sistem configuration
23
8
Sensor source (EZ)
Real/Sim flag
24
8
Sensor avaible
Lag length
25
16
BIN 1 distance
Number of beams
26
16
XMIT pulse length based on (WT)
Number of cells (WN)
27
8
WP ref layer average (starting cell)
Pings per ensemble (WP)
28
8
WP ref layer average (ending cell)
Depth cell lentgh (WS)
29
8
False target tresh (WA)
Blank after transmit (WF)
30
8
Spare
Profiling Mode (WM)
31
16
Transmit lag distance
Low corr thresh (WC)
32
64
Cpu board serial number
No. Code reps
33
16
Sistem bandwidth (WB)
%GD minimum (WG)
34
8
Sistem power (CQ)
Error velocity Maximum (WE)
35
8
Spare
Tpp minutes
36
32
Instrument serial number
Tpp seconds
37
8
Beam angle
Tpp Hundredths (TP)

Filas

Tabla 6-3. Estructura de la matriz de datos de velocidad registrada por el equipo.

N de Identif Vel celda 1

ensam icador Beam 1
Vel celda 1
1
.
Beam 1
Vel celda 1
2
.
Beam 1
Vel celda 1
.
.
Beam 1
Vel celda 1
.
.
Beam 1
Vel celda 1
.
.
Beam 1
Vel celda 1
.
.
Beam 1
Vel celda 1
N
N
Beam 1

Vel celda 1
Beam 2
Vel celda 1
Beam 2
Vel celda 1
Beam 2
Vel celda 1
Beam 2
Vel celda 1
Beam 2
Vel celda 1
Beam 2
Vel celda 1
Beam 2
Vel celda 1
Beam 2

Vel celda 1
Beam 3
Vel celda 1
Beam 3
Vel celda 1
Beam 3
Vel celda 1
Beam 3
Vel celda 1
Beam 3
Vel celda 1
Beam 3
Vel celda 1
Beam 3
Vel celda 1
Beam 3

Columnas
Vel celda 1 Vel celda 2
Beam 4
Beam 1
Vel celda 1 Vel celda 2
Beam 4
Beam 1
Vel celda 1 Vel celda 2
Beam 4
Beam 1
Vel celda 1 Vel celda 2
Beam 4
Beam 1
Vel celda 1 Vel celda 2
Beam 4
Beam 1
Vel celda 1 Vel celda 2
Beam 4
Beam 1
Vel celda 1 Vel celda 2
Beam 4
Beam 1
Vel celda 1 Vel celda 2
Beam 4
Beam 1

Vel celda M
Beam 4
Vel celda M
Beam 4
Vel celda M
Beam 4
Vel celda M
Beam 4
Vel celda M
Beam 4
Vel celda M
Beam 4
Vel celda M
Beam 4
Vel celda M
Beam 4

71

Tabla 6-4. Estructura de ordenamiento de las columnas para la matriz Bottom track data.

Posicion
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30

Bottom-track data
N de bits
Nombre
Posicion N de bits
16
Bottom track id
31
16
16
Bt pings per ensemble (BP)
32
16
16
Bt delay beforre re-acquire (BD)
33
16
8
Bt corr mag min (BC)
34
16
8
Bt eval amp min (BA)
35
16
8
Bt percent good min (BG)
36
16
8
Bt mode (BM)
37
8
16
Bt err vel max (BE)
38
8
32
Reserved
39
8
16
Beam #1 range
40
8
16
Beam #2 range
41
8
16
Beam #3 range
42
8
16
Beam #4 range
43
8
16
Beam #1 bt vel
44
8
16
Beam #2 bt vel
45
8
16
Beam #3 bt vel
46
8
16
Beam #4 bt vel
47
8
8
Beam #1 bt corr
48
8
8
Beam #2 bt corr
49
16
8
Beam #3 bt corr
50
8
8
Beam #4 bt corr
51
8
8
Beam #1 bt eval amp
52
8
8
Beam #2 bt eval amp
53
8
8
Beam #3 bt eval amp
54
8
8
Beam #4 bt eval amp
55
8
8
Beam #1 bt %good
56
8
8
Beam #2 bt %good
57
8
8
Beam #3 bt %good
58
8
8
Beam #4 bt %good
59
32
16
Ref layer min (BL)

Nombre
Ref layer near (BL)
Ref layer far (BL)
Beam #1 ref layer vel
Beam #2 ref layer vel
Beam #3 ref layer vel
Beam #4 ref layer vel
Beam #1 ref corr
Beam #2 ref corr
Beam #3 ref corr
Beam #4 ref corr
Beam #1 ref int
Beam #2 ref int
Beam #3 ref int
Beam #4 ref int
Beam #1 ref %good
Beam #2 ref %good
Beam #3 ref %good
Beam #4 ref %good
Bt max depth (BX)
Bm#1 rssi amp
Bm#2 rssi amp
Bm#3 rssi amp
Bm#4 rssi amp
Gain
Beam #1 range msb
Beam #2 range msb
Beam #3 range msb
Beam #4 range msb
Reserved

72

6.2 Anexo 2, Perfiles y datos utilizados para el clculo del esfuerzo de

corte.
6.2.1

Correspondencia entre los archivos del ADCP y los perfiles medios obtenidos

Tabla 6-5. Correspondencia de los archivos del ADCP con los puntos registrados por el GPS

Archivo
BB P1 0_000.PDO
BB P1 0_000_SBT.PDO
BB P1 0_001.PDO
BB P1 0_001_SBT.PDO
STATIONBB_25_TARDE_0_000.PDO
STATIONBB_25_TARDE_0_000_SBT.PDO
STATIONBB_25_TARDE_0_001.PDO
STATIONBB_25_TARDE_0_001_SBT.PDO
Station_desembocadura_0_000.PDO
Station_desembocadura_0_000_SBT.PDO
Station_desembocadura_0_001.PDO
Station_desembocadura_0_001_SBT.PDO
Station_desembocadura_0_002.PDO
Station_desembocadura_0_002_SBT.PDO
Station_desembocadura_0_003.PDO
Station_desembocadura_0_003_SBT.PDO
Station_desembocadura_0_004.PDO
Station_desembocadura_0_004_SBT.PDO
Station_desembocadura_0_005.PDO
Station_desembocadura_0_005_SBT.PDO
StationBBDGA_0_000.PDO
StationBBDGA_0_000_SBT.PDO
StationBBDGA_0_001.PDO
StationBBDGA_0_001_SBT.PDO
StationBBDGA_0_002.PDO
StationBBDGA_0_002_SBT.PDO
StationBBDGA_0_003.PDO
StationBBDGA_0_003_SBT.PDO
StationBBDGA_0_004.PDO
StationBBDGA_0_004_SBT.PDO
Seccion 1_0_007_SBT.PDO
Seccion 1_0_009_SBT.PDO
Station_0_004_SBT.PDO
Station_0_005_SBT.PDO
Station_0_006_SBT.PDO
Seccion 1 P2_0_000_SBT.PDO
Seccion 1 P2_0_001_SBT.PDO
Seccion 1 P2_0_002_SBT.PDO
Seccion 1 P2_0_003_SBT.PDO
Seccion 1 P2_0_004_SBT.PDO
sseccion 2_0_004_SBT.PDO

Pto
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41

Fecha
Miercoles 25-7
Miercoles 25-7
Miercoles 25-7
Miercoles 25-7
Miercoles 25-7
Miercoles 25-7
Miercoles 25-7
Miercoles 25-7
Jueves 26-7
Jueves 26-7
Jueves 26-7
Jueves 26-7
Jueves 26-7
Jueves 26-7
Jueves 26-7
Jueves 26-7
Jueves 26-7
Jueves 26-7
Jueves 26-7
Jueves 26-7
Viernes 27-7
Viernes 27-7
Viernes 27-7
Viernes 27-7
Viernes 27-7
Viernes 27-7
Viernes 27-7
Viernes 27-7
Viernes 27-7
Viernes 27-7
viernes 16-11
viernes 16-11
jueves 15-11
jueves 15-11
jueves 15-11
viernes 16-11
viernes 16-11
viernes 16-11
viernes 16-11
viernes 16-11
viernes 16-11

73

6.2.2

Perfiles estticos; Velocidad y error asociado

Los siguientes grficos corresponden a los perfiles de velocidad

y su error asociado

(mismas unidades que la velocidad).

a) Velocidad

Profundidad (m)

-0.5

-1

0.5
-1.5

20

40

60

80

100

120

140

160

180

200

220

Ensamble

Profundidad (m)

b) Error
0

0.4

-0.5

0.2

-1

-1.5

-0.2
20

40

60

80

100

120

140

160

180

200

220

Ensamble

Figura 6-1. Punto esttico 1.

a) Velocidad

Profundidad (m)

-0.5
-1

0.5
-1.5

100

200

300

400

500

600

700

800

900

Ensamble

Profundidad (m)

b) Error
0

0.4

-0.5

0.2
0

-1

-0.2
-1.5

-0.4
100

200

300

400

500

600

700

800

900

Ensamble

Figura 6-2. Punto esttico 2.

74

a) Velocidad

Profundidad (m)

-0.5

-1

0.5
-1.5

50

100

150

200

250

Ensamble

b) Error

Profundidad (m)

0.6
0.4
0.2
0
-0.2
-0.4

-0.5

-1

-1.5

50

100

150

200

250

Ensamble

Figura 6-3. Punto esttico 3.

a) Velocidad

Profundidad (m)

1.5

-0.5

1
-1

0.5
-1.5

100

200

300

400

500

600

700

800

900

Ensamble

b) Error

Profundidad (m)

0.6
0.4
0.2
0
-0.2
-0.4

-0.5

-1

-1.5

100

200

300

400

500

600

700

800

900

Ensamble

Figura 6-4. Punto esttico 4.

75

a) Velocidad

Profundidad (m)

1.5

-0.5

1
-1

0.5

-1.5
-2
50

100

150

200

250

Ensamble

b) Error

Profundidad (m)

0.5

-0.5
-1

-1.5

-0.5

-2
50

100

150

200

250

Ensamble

Figura 6-5.Punto esttico 5.

a) Velocidad

Profundidad (m)

1.5

-0.5

1
-1

0.5

-1.5
-2
100

200

300

400

500

600

700

800

900

1000

Ensamble

b) Error

Profundidad (m)

0.5

-0.5
-1

-1.5

-0.5

-2
100

200

300

400

500

600

700

800

900

1000

Ensamble

Figura 6-6.Punto esttico 6.

76

a) Velocidad

Profundidad (m)

-0.5
-1

0.5
-1.5
-2
50

100

150

200

250

300

Ensamble

Profundidad (m)

b) Error
0

0.4

-0.5

0.2
0

-1

-0.2

-1.5

-0.4
-2
50

100

150

200

250

300

Ensamble

Figura 6-7. Punto esttico 7.

a) Velocidad

Profundidad (m)

-0.5
-1

0.5
-1.5
-2
100

200

300

400

500

600

700

Ensamble

b) Error
0.5

Profundidad (m)

0
-0.5

-1
-1.5

-0.5

-2
100

200

300

400

500

600

700

Ensamble

Figura 6-8. Punto esttico 8.

77

a) Velocidad

Profundidad (m)

1.5

-0.5

-1

0.5

-1.5
50

100

150

200

250

300

Ensamble

b) Error
0

Profundidad (m)

0.5
-0.5

0
-1

-0.5

-1.5
50

100

150

200

250

300

Ensamble

Figura 6-9.Punto esttico 9.

a) Velocidad

Profundidad (m)

1.5
-0.5

1
-1

0.5
-1.5
100

200

300

400

500

600

700

800

900

Ensamble

b) Error

Profundidad (m)

0.5

-0.5

-1

-0.5

-1.5
100

200

300

400

500

600

700

800

900

Ensamble

Figura 6-10. Punto esttico 10.

78

a) Velocidad
0

Profundidad (m)

-0.2
-0.4

0.5

-0.6
-0.8
-1
50

100

150

200

250

300

Ensamble

b) Error
0

0.4

Profundidad (m)

-0.2
-0.4

0.2

-0.6

-0.8

-0.2

-1

-0.4
50

100

150

200

250

300

Ensamble

Figura 6-11. Punto esttico 11.

a) Velocidad
0

Profundidad (m)

-0.2
-0.4

0.5

-0.6
-0.8
-1
100

200

300

400

500

600

700

800

900

Ensamble

b) Error
0

0.4

Profundidad (m)

-0.2

0.2

-0.4
-0.6

-0.8

-0.2

-1

-0.4
100

200

300

400

500

600

700

800

900

Ensamble

Figura 6-12. Punto esttico 12.

79

a) Velocidad

Profundidad (m)

-0.5

-1

0.5

-1.5

20

40

60

80

100

120

140

160

180

200

Ensamble

b) Error

Profundidad (m)

0.4

-0.5

0.2
0

-1

-0.2

-1.5

-0.4
20

40

60

80

100

120

140

160

180

200

Ensamble

Figura 6-13. Punto esttico 13.

a) Velocidad

Profundidad (m)

1.2
1
0.8
0.6
0.4
0.2

-0.5

-1

-1.5

100

200

300

400

500

600

700

800

900

1000

Ensamble

b) Error
0

Profundidad (m)

0.4
-0.5

0.2
0

-1

-0.2

-1.5

-0.4
100

200

300

400

500

600

700

800

900

1000

Ensamble

Figura 6-14. Puntos esttico 14.

80

a) Velocidad
0

Profundidad (m)

-0.2
-0.4

0.5

-0.6
-0.8
-1
50

100

150

200

250

300

Ensamble

b) Error

Profundidad (m)

0
-0.2

0.2

-0.4

-0.6

-0.2

-0.8

-0.4

-1
50

100

150

200

250

300

Ensamble

Figura 6-15. Punto esttico 15.

a) Velocidad
0

Profundidad (m)

-0.2
-0.4

0.5

-0.6
-0.8
-1
100

200

300

400

500

600

700

800

Ensamble

b) Error
0

0.4

Profundidad (m)

-0.2

0.2

-0.4
-0.6

-0.8

-0.2

-1
100

200

300

400

500

600

700

800

-0.4

Ensamble

Figura 6-16. Punto esttico 16.

81

a) Velocidad

Profundidad (m)

1
0.8
0.6
0.4
0.2

-0.5

-1

-1.5
20

40

60

80

100

120

140

160

180

200

220

Ensamble

b) Error

Profundidad (m)

0.2

-0.5

0
-1

-0.2
-1.5

-0.4
20

40

60

80

100

120

140

160

180

200

220

Ensamble

Figura 6-17. Punto esttico 17.

a) Velocidad
0

Profundidad (m)

1
-0.5

0.5

-1

-1.5
100

200

300

400

500

600

700

800

Ensamble

Profundidad (m)

b) Error
0

0.4

-0.5

0.2
0

-1

-0.2
-1.5

-0.4
100

200

300

400

500

600

700

800

Ensamble

Figura 6-18. Punto esttico 18.

82

a) Velocidad

Profundidad (m)

-0.5

0.5
-1

50

100

150

200

250

Ensamble

b) Error
0

Profundidad (m)

0.4
0.2

-0.5

0
-1

-0.2
50

100

150

200

250

Ensamble

Figura 6-19. Punto esttico 19.

a) Velocidad

Profundidad (m)

-0.5

0.5

-1

-1.5
50

100

150

200

250

300

350

Ensamble

b) Error

Profundidad (m)

0.2

-0.5

0
-1

-0.2
-1.5

-0.4
50

100

150

200

250

300

350

Ensamble

Figura 6-20. Punto esttico 20.

83

a) Velocidad

Profundidad (m)

-0.5

-1
-1.5

-2
50

100

150

200

250

300

350

400

450

Ensamble

b) Error

Profundidad (m)

-0.5
-1

-1.5

-1

-2
50

100

150

200

250

300

350

400

450

Ensamble

Figura 6-21. Punto esttico 21.

a) Velocidad

Profundidad (m)

-0.5
-1

-1.5
-2
200

400

600

800

1000

1200

1400

1600

Ensamble

b) Error

Profundidad (m)

-0.5
-1

-1.5

-1

-2
200

400

600

800

1000

1200

1400

1600

Ensamble

Figura 6-22. Punto esttico 22.

84

a) Velocidad

Profundidad (m)

-0.5

-1

2
1

-1.5
50

100

150

200

250

300

350

400

Ensamble

b) Error
0

Profundidad (m)

1
-0.5

-1

-1

-1.5
50

100

150

200

250

300

350

400

Ensamble

Figura 6-23. Punto esttico 23.

a) Velocidad

Profundidad (m)

-0.5

-1

-1.5

200

400

600

800

1000

1200

1400

1600

Ensamble

b) Error
0

Profundidad (m)

1
-0.5

-1

-1.5

-1
200

400

600

800

1000

1200

1400

1600

Ensamble

Figura 6-24. Punto esttico 24.

85

a) Velocidad

Profundidad (m)

0
-0.5

-1

-1.5
-2
50

100

150

200

250

300

350

400

Ensamble

b) Error

Profundidad (m)

-0.5
-1

-1.5

-1

-2
50

100

150

200

250

300

350

400

Ensamble

Figura 6-25. Punto esttico 25.

a) Velocidad

Profundidad (m)

0
-0.5

-1

-1.5

-2
200

400

600

800

1000

1200

1400

Ensamble

b) Error

Profundidad (m)

-0.5
-1

-1.5

-1

-2
200

400

600

800

1000

1200

1400

Ensamble

Figura 6-26. Punto esttico 26.

86

a) Velocidad
0

Profundidad (m)

-0.2
-0.4

-0.6

-0.8
-1
50

100

150

200

250

300

350

400

Ensamble

b) Error
0

Profundidad (m)

-0.2
-0.4

-0.6
-0.8

-1

-1
50

100

150

200

250

300

350

400

Ensamble

Figura 6-27. Punto esttico 27.

a) Velocidad
0

Profundidad (m)

-0.2
-0.4

-0.6

-0.8
-1
200

400

600

800

1000

1200

1400

1600

1800

Ensamble

b) Error
0

Profundidad (m)

-0.2
-0.4

-0.6
-0.8

-1

-1
200

400

600

800

1000

1200

1400

1600

1800

Ensamble

Figura 6-28. Punto esttico 28.

87

a) Velocidad
0

Profundidad (m)

3
-0.5

2
1

-1

50

100

150

200

250

300

350

400

450

Ensamble

b) Error
0

Profundidad (m)

1
-0.5

0
-1

-1
50

100

150

200

250

300

350

400

450

Ensamble

Figura 6-29. Punto esttico 29.

a) Velocidad
4

Profundidad (m)

-0.5

2
-1

200

400

600

800

1000

1200

1400

1600

Ensamble

b) Error
0

Profundidad (m)

1
-0.5

0
-1

-1
200

400

600

800

1000

1200

1400

1600

Ensamble

Figura 6-30. Punto esttico 30.

88

a) Velocidad

Profundidad (m)

-0.5
-1

0.5

-1.5
-2
20

40

60

80

100

120

Ensamble

b) Error

Profundidad (m)

0.2

-0.5
-1

-1.5

-0.2

-2
20

40

60

80

100

120

Ensamble

Figura 6-31. Punto esttico 31.

a) Velocidad
0

Profundidad (m)

-0.5
-1

0.5

-1.5
-2

10

20

30

40

50

60

70

80

Ensamble

b) Error
0

0.2

Profundidad (m)

-0.5
-1

-1.5
-2

-0.2
10

20

30

40

50

60

70

80

Ensamble

Figura 6-32. Punto esttico 32.

89

a) Velocidad
0

Profundidad (m)

-0.5

-1
-1.5

-2
-2.5
10

20

30

40

50

60

70

80

90

100

Ensamble

b) Error
0

Profundidad (m)

-0.5
-1

-1.5
-2

-1

-2.5
10

20

30

40

50

60

70

80

90

100

Ensamble

Figura 6-33. Punto esttico 33.

Profundidad (m)

a) Velocidad
0

-1

-2

-3
10

20

30

40

50

60

Ensamble

b) Error
0

Profundidad (m)

1
-1

0
-2

-1
-3
10

20

30

40

50

60

Ensamble

Figura 6-34. Punto esttico 34.

90

a) Velocidad
0

Profundidad (m)

3
-1

2
-2

-3
20

40

60

80

100

120

140

Ensamble

b) Error
0

Profundidad (m)

1
-1

0
-2

-1
-3
20

40

60

80

100

120

140

Ensamble

Figura 6-35. Punto esttico 35.

a) Velocidad

Profundidad (m)

1.2
1
0.8
0.6
0.4
0.2

-1

-2

-3
10

20

30

40

50

60

70

Ensamble

b) Error

Profundidad (m)

0.2
-1

-2

-3
10

20

30

40

50

60

70

-0.2

Ensamble

Figura 6-36. Punto esttico 36.

91

a) Velocidad
0

Profundidad (m)

-0.5

-1
-1.5

0.5

-2
-2.5
5

10

15

20

25

30

35

40

Ensamble

b) Error
0

0.2

Profundidad (m)

-0.5
-1

-1.5
-2

-0.2

-2.5
5

10

15

20

25

30

35

40

Ensamble

Figura 6-37. Punto esttico 37.

a) Velocidad
0

1.2
1
0.8
0.6
0.4
0.2

Profundidad (m)

-0.5
-1
-1.5
-2
-2.5
20

40

60

80

100

120

Ensamble

b) Error
0

Profundidad (m)

-0.5

0.2

-1

-1.5
-2

-0.2

-2.5
20

40

60

80

100

120

Ensamble

Figura 6-38. Punto esttico 38.

92

a) Velocidad
0

Profundidad (m)

-0.5

-1

0.5

-1.5
-2
10

20

30

40

50

60

70

80

90

Ensamble

b) Error
0

0.2

Profundidad (m)

-0.5
-1

-1.5
-2

-0.2
10

20

30

40

50

60

70

80

90

Ensamble

Figura 6-39. Punto esttico 39.

a) Velocidad

Profundidad (m)

1
0.8

-0.5

0.6
0.4
0.2

-1
-1.5
-2
10

20

30

40

50

60

70

80

90

Ensamble

b) Error

Profundidad (m)

0.4

-0.5

0.2

-1

0
-1.5

-0.2

-2
10

20

30

40

50

60

70

80

90

Ensamble

Figura 6-40. Punto esttico 40.

93

a) Velocidad

Profundidad (m)

-0.5
-1

0.5
-1.5
-2
10

20

30

40

50

60

70

80

Ensamble

b) Error

Profundidad (m)

0.2

-0.5

-1
-1.5

-0.2

-2
10

20

30

40

50

60

70

80

Ensamble

Figura 6-41. Punto esttico 41.

94

Perfiles medios de velocidad

-0.6

-0.6

-0.8

-0.8

-0.8

-1

-1

-1

-1.2
-1.4

Profundidad (m)

-0.6

Profundidad (m)

Profundidad (m)

6.2.3

-1.2
-1.4

-1.2
-1.4

-1.6

-1.6

-1.6

-1.8

-1.8

-1.8

-2

-2

-2
0

0.5

1.5

0.5

1.5

0.5

1.5

Velocidad (m/s)

Velocidad (m/s)

Velocidad (m/s)

Figura 6-42. Punto esttico 1.

Figura 6-43. Punto esttico 2.

Figura 6-44. Punto esttico 3.

-0.4

-0.5

-0.5

-1

-1

-1
-1.2
-1.4

Profundidad (m)

-0.8

Profundidad (m)

Profundidad (m)

-0.6

-1.5

-2

-1.5

-2

-1.6
-1.8

-2.5
0

0.5

1.5

-2.5
0

0.5

0.5

Velocidad (m/s)

Velocidad (m/s)

Velocidad (m/s)

Figura 6-45. Punto esttico 4.

Figura 6-46. Punto esttico 5.

Figura 6-47. Punto esttico 6.

-0.5

-0.5

-0.2

-1.5

-1

Profundidad (m)

-1

Profundidad (m)

Profundidad (m)

-0.4

-1.5

-0.6
-0.8
-1
-1.2
-1.4

-2

-2
0

0.5

-1.6
0

0.5

0.5

Velocidad (m/s)

Velocidad (m/s)

Velocidad (m/s)

Figura 6-48. Punto esttico 7.

Figura 6-49. Punto esttico 8.

Figura 6-50. Punto esttico 9.

95

-0.4

-0.4

-0.5

-0.5

-0.6

-0.6

-0.6

-0.8
-1

Profundidad (m)

-0.4

Profundidad (m)

Profundidad (m)

-0.2

-0.7
-0.8

-0.7
-0.8

-1.2

-0.9

-0.9

-1.4

-1

-1

-1.6

-1.1
0

0.5

-1.1
0

0.2

0.4

0.6

0.8

0.2

0.4

0.6

0.8

Velocidad (m/s)

Velocidad (m/s)

Velocidad (m/s)

Figura 6-51. Punto esttico 10.

Figura 6-52. Punto esttico 11.

Figura 6-53. Punto esttico 12.

-0.5

-0.5

-1

-1

-0.5

-1.5

Profundidad (m)

Profundidad (m)

Profundidad (m)

-0.6

-1.5

-0.7

-0.8

-0.9

-2

-2
0

0.5

-1
0

0.5

0.2

0.4

0.6

0.8

Velocidad (m/s)

Velocidad (m/s)

Velocidad (m/s)

Figura 6-54. Punto esttico 13.

Figura 6-55. Punto esttico 14.

Figura 6-56. Punto esttico 15.

-0.5

-0.4

-0.4

-0.6

-0.6

-0.8

-0.8

-0.7

-0.8

Profundidad (m)

Profundidad (m)

Profundidad (m)

-0.6

-1
-1.2

-1
-1.2

-1.4

-1.4

-1.6

-1.6

-0.9

-1

-1.8
0

0.5

1.5

-1.8
0

0.2

0.4

0.6

0.8

0.2

0.4

0.6

0.8

Velocidad (m/s)

Velocidad (m/s)

Velocidad (m/s)

Figura 6-57. Punto esttico 16.

Figura 6-58. Punto esttico 17.

Figura 6-59. Punto esttico 18.

96

-0.5

-0.4

-0.6

-0.6
-0.5

-0.8
-0.9
-1

-0.8

Profundidad (m)

Profundidad (m)

-0.7

Profundidad (m)

-1
-1.2

-1

-1.5

-1.4

-1.1

-2
-1.6

-1.2
-1.3

-1.8
0

0.2

0.4

0.6

0.8

-2.5
0

0.2

0.4

0.6

0.8

Velocidad (m/s)

Velocidad (m/s)

Velocidad (m/s)

Figura 6-60. Punto esttico 19.

Figura 6-61. Punto esttico 20.

Figura 6-62. Punto esttico 21.

-0.4

-0.4

-0.6

-0.6

-0.8

-0.8

-1

-1.5

Profundidad (m)

Profundidad (m)

Profundidad (m)

-0.5

-1
-1.2

-1
-1.2

-1.4

-1.4

-1.6

-1.6

-2

-2.5

-1.8
0

0.5

-1.8
0

0.5

1.5

0.5

1.5

Velocidad (m/s)

Velocidad (m/s)

Velocidad (m/s)

Figura 6-63. Punto esttico 22.

Figura 6-64. Punto esttico 23.

Figura 6-65. Punto esttico 24.

-0.5

-0.5

-0.4

-1

-1.5

Profundidad (m)

Profundidad (m)

Profundidad (m)

-0.5

-1

-1.5

-0.6
-0.7
-0.8
-0.9

-2

-2
-1

-2.5

-2.5
0

0.5

-1.1
0

0.5

0.5

1.5

Velocidad (m/s)

Velocidad (m/s)

Velocidad (m/s)

Figura 6-66. Punto esttico 25.

Figura 6-67. Punto esttico 26.

Figura 6-68. Punto esttico 27.

97

-0.4

-0.4

-0.4

-0.6

-0.6

-0.7
-0.8

Profundidad (m)

-0.6

Profundidad (m)

Profundidad (m)

-0.5

-0.8

-1

-0.8

-1

-0.9
-1.2

-1.2

-1
-1.1

-1.4
0

-1.4
0

Velocidad (m/s)

-0.5

-1

-1

-1.5

-0.5

-1.5

0.5

-1.5

-2.5

-1

-2

-2.5

-2.5

0.5

Figura 6-73. Punto esttico 32.

-0.5

-0.5

-0.5

-1

Profundidad (m)

-2
-2.5

Profundidad (m)

-1.5

-1
-1.5
-2

-3.5
1

Velocidad (m/s)

Figura 6-75. Punto esttico 34.

-1.5
-2
-2.5

-3.5

-3
0

1.5

-3

-2.5

-3

Figura 6-74. Punto esttico 33.

-1

0.5

Velocidad (m/s)

Velocidad (m/s)

Velocidad (m/s)

Figura 6-72. Punto esttico 31.

1.5

-2

-2

Figura 6-71. Punto esttico 30.

Profundidad (m)

-0.5

0.5

Velocidad (m/s)

Figura 6-70. Punto esttico 29.

Profundidad (m)

Profundidad (m)

Velocidad (m/s)

Figura 6-69. Punto esttico 28.

Profundidad (m)

Velocidad (m/s)

Figura 6-76. Punto esttico 35.

0.5

1.5

Velocidad (m/s)

Figura 6-77. Punto esttico 36.

98

-0.5

-0.5
-1

-1.5

-2

-1

Profundidad (m)

Profundidad (m)

-1

Profundidad (m)

-0.5

-1.5
-2

-2

-2.5

-2.5

-3
0

0.5

1.5

-2.5
0

Velocidad (m/s)

0.5

1.5

Velocidad (m/s)

Figura 6-78. Punto esttico 37.

0.5

1.5

Velocidad (m/s)

Figura 6-79. Punto esttico 38.

-0.4

-1.5

Figura 6-80. Punto esttico 39.

-0.5

-1

-0.8

Profundidad (m)

Profundidad (m)

-0.6

-1
-1.2
-1.4

-1.5

-2

-1.6
-1.8

-2.5
0

0.5

1.5

Velocidad (m/s)

Figura 6-81. Punto esttico 40.

0.5

1.5

Velocidad (m/s)

Figura 6-82. Punto esttico 41.

99

Perfiles de esfuerzo de corte calculados con el TKE.

-0.8

Profundidad (m)

Profundidad (m)

-0.8
-1
-1.2
-1.4

-0.8

Profundidad (m)

6.2.4

-1
-1.2
-1.4

-1.6
4.8

5.1

5.2

-1.6

5.3

6.5

Esfuerzo de corte (N/m2)

7.5

Esfuerzo de corte (N/m2)

Figura 6-84. Punto esttico 2.

Profundidad (m)

-0.8
-1
-1.2
-1.4

Esfuerzo de corte (N/m2)

Figura 6-85. Punto esttico 3.

-0.6

-0.6

-0.8

-0.8

Profundidad (m)

Figura 6-83. Punto esttico 1.

Profundidad (m)

-1.2
-1.4

-1.6
4.9

-1

-1
-1.2
-1.4
-1.6

-1
-1.2
-1.4
-1.6

-1.6
3.5

4.5

5.5

Esfuerzo de corte (N/m2)

10

11

-0.6

-0.8

-0.8

-1
-1.2
-1.4
-1.6
5

5.5

-1
-1.2
-1.4

-1.2

6.5

10

10

12

14
2

Figura 6-92. Punto esttico 10.

14

16

Figura 6-91. Punto esttico 9.

-0.5

-0.5

-0.55

-0.55
Profundidad (m)

-0.6
-0.65
-0.7
-0.75
-0.8

Esfuerzo de corte (N/m )

12

Esfuerzo de corte (N/m )

Figura 6-90. Punto esttico 8.

Profundidad (m)

-1

-1

Esfuerzo de corte (N/m2)

-0.6
-0.8

-0.8

-1.2

Esfuerzo de corte (N/m2)

Figura 6-89. Punto esttico 7.

11

-0.6

-1.6
4.5

10

Figura 6-88. Punto esttico 6.

Profundidad (m)

-0.6

Esfuerzo de corte (N/m2)

Figura 6-87. Punto esttico 5.

Profundidad (m)

Profundidad (m)

Figura 6-86. Punto esttico 4.

Profundidad (m)

Esfuerzo de corte (N/m2)

-0.6
-0.65
-0.7
-0.75
-0.8

4.8

5.2

5.4

5.6

5.8
2

Esfuerzo de corte (N/m )

Figura 6-93. Punto esttico 11.

8.5

9.5
2

Esfuerzo de corte (N/m )

Figura 6-94. Punto esttico 12.

100

-0.5

-0.5

-0.6

-1
-1.2
-1.4

Profundidad (m)

-0.8

Profundidad (m)

Profundidad (m)

-0.55

-1

-0.6
-0.65
-0.7
-0.75
-0.8

-1.6

-1.5
5

10

6.4

6.5

6.6

6.7

Esfuerzo de corte (N/m )

Figura 6-95. Punto esttico 13.

6.5

7.5

Esfuerzo de corte (N/m )

Figura 6-96. Punto esttico 14.

-0.5

5.5

Esfuerzo de corte (N/m )

Figura 6-97. Punto esttico 15.

-0.5

-0.5

-0.65
-0.7
-0.75

Profundidad (m)

-0.6

Profundidad (m)

Profundidad (m)

-0.55

-1

-1

-0.8
-1.5
4.5

5.5

6.5

-1.5
6.5

7.5

8.5

Esfuerzo de corte (N/m )

Figura 6-98. Punto esttico 16.

Figura 6-100. Punto esttico 18.

-0.5
-0.6

-0.7
-0.8

Profundidad (m)

Profundidad (m)

Profundidad (m)

-0.6

-1

-0.9

-0.8
-1
-1.2
-1.4

-1

-1.5
5

5.5

6.5

-1.6
4.6

Figura 6-101. Punto esttico 19.

4.8

5.2

5.4

55

Esfuerzo de corte (N/m )

60

65

70

75

80

Esfuerzo de corte (N/m )

Esfuerzo de corte (N/m )

Figura 6-102. Punto esttico 20.

Figura 6-103. Punto esttico 21.

-0.5

-0.6

-1
-1.2
-1.4

Profundidad (m)

-0.6

-0.8

Profundidad (m)

Profundidad (m)

Esfuerzo de corte (N/m )

Figura 6-99. Punto esttico 17.

-0.5

6.5
2

Esfuerzo de corte (N/m )

-1

-0.8
-1
-1.2

-1.6
-1.4

-1.5
50

55

60
2

Esfuerzo de corte (N/m )

Figura 6-104. Punto esttico 22.

60

80

100

120

140

160

Esfuerzo de corte (N/m )

Figura 6-105. Punto esttico 23.

55

60

65

70
2

Esfuerzo de corte (N/m )

Figura 6-106. Punto esttico 24

101

-0.8

-1
-1.2
-1.4
-1.6
-1.8

-0.45
-0.5
Profundidad (m)

-0.6

-0.8

Profundidad (m)

Profundidad (m)

-0.6

-1
-1.2
-1.4
-1.6

45

50

-0.6
-0.65
-0.7

-1.8
40

-0.55

-0.75

55

50

60

70

80

90

100

Esfuerzo de corte (N/m )

Figura 6-107. Punto esttico 25.

150

200

250
2

Esfuerzo de corte (N/m )

Esfuerzo de corte (N/m )

Figura 6-108. Punto esttico 26.

Figura 6-109. Punto esttico 27.

-0.5

-0.5

-0.6

-0.6

-0.6

-0.7

Profundidad (m)

Profundidad (m)

Profundidad (m)

-0.5
-0.7
-0.8
-0.9
-1

-0.8
60

80

100

120

140

160

-0.7
-0.8
-0.9
-1

60

80

100

120

70

Esfuerzo de corte (N/m )

Figura 6-110. Punto esttico 28.

80

90

100

110

120

Esfuerzo de corte (N/m )

Esfuerzo de corte (N/m )

Figura 6-111. Punto esttico 29.

Figura 6-112. Punto esttico 30.

-0.5

-1
-1.2
-1.4
-1.6

-1

Profundidad (m)

Profundidad (m)

Profundidad (m)

-0.8

-1.5

Esfuerzo de corte (N/m2)

-0.5

-0.5

-1

-1

-1.5
-2

-1.5
-2

60

70

Esfuerzo de corte (N/m2)

Figura 6-116. Punto esttico 34.

45

50

55

60

65

Esfuerzo de corte (N/m2)

Figura 6-115. Punto esttico 33.

-1
-1.5
-2
-2.5

-2.5

-2.5
50

40

Profundidad (m)

Figura 6-114. Punto esttico 32.

Profundidad (m)

Profundidad (m)

Esfuerzo de corte (N/m2)

Figura 6-113. Punto esttico 31.

40

-1.5

-2

-2

-1.8

-1

45

50

55

60

Esfuerzo de corte (N/m2)

Figura 6-117. Punto esttico 35.

10

Esfuerzo de corte (N/m2)

Figura 6-118. Punto esttico 36.

102

-1.5

-2

-1

Profundidad (m)

Profundidad (m)

Profundidad (m)

-0.8
-1

-1.5

-2
29

30

31

32

33

-1.4
-1.6
-1.8

34

6.6

Esfuerzo de corte (N/m2)

6.7

6.8

6.9

7.1

Esfuerzo de corte (N/m2)

Figura 6-119. Punto esttico 37.

-1
-1.2

Esfuerzo de corte (N/m2)

Figura 6-120. Punto esttico 38.

Figura 6-121. Punto esttico 39.

-0.8

-0.8

Profundidad (m)

Profundidad (m)

-0.6

-1
-1.2

-1
-1.2
-1.4
-1.6
-1.8

10

11

Esfuerzo de corte (N/m2)

Figura 6-122. Punto esttico 40.

2.5

3.5

Esfuerzo de corte (N/m2)

Figura 6-123. Punto esttico 41.

103