See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/350939694

Estimation model of the Manning roughness coefficient as a function of

granulometry in the Santa River, Recuay-Carhuaz sector, Áncash, Perú

Article · April 2021

CITATIONS READS

0 100

4 authors, including:

Abelardo Diaz Edilberto Guevara-Perez

Universidad Nacional Santiago Antúnez de Mayolo, Huaraz Universidad de Carabobo, UC
4 PUBLICATIONS   0 CITATIONS    28 PUBLICATIONS   51 CITATIONS   

SEE PROFILE SEE PROFILE

Some of the authors of this publication are also working on these related projects:

Doctoral Thesis View project

All content following this page was uploaded by Abelardo Diaz on 17 April 2021.

The user has requested enhancement of the downloaded file.

UC Universidad
de Carabobo Díaz-Salas et al. / Revista Ingeniería UC, Vol. 27, No 3, Diciembre, 2020 328 – 342

Estimation model of the Manning roughness coefficient as a function of

granulometry in the Santa River, Recuay–Carhuaz sector, Áncash, Perú
Abelardo Díaz-Salas ∗ iD , Edilberto Guevara-Pérez iD , Juan Rosales-Cueva iD ,
Escuela Profesional de Ingeniería Agrícola. Universidad Nacional Santiago Antúnez de Mayolo. Huaraz, Perú

Abstract.- Given that the Manning equation is one of the most used application methods to determine the flow in ungauged
rivers, an investigation was carried out to establish a model to determine the Roughness coefficcient (n), in a mountainous sector
of the river Santa in the department of Ancash, Peru. The novelty of the model is its simplicity since it only uses the effective
diameter (di), of the river material obtained from the granulometric analysis as an input variable. Three estimation procedures
were used, speed measurement, Manning’s equation and Cowan’s model, and for the analysis, linear and non-linear correlation
and regression models; goodness of fit was performed by the Student test, t. The results were validated with field measurements
of the hydraulic characteristics of the section at eight sampling points and granulometric analysis of the bed material samples.
The best fit model for estimating the Manning roughness coefficient is the potential model, being its mathematical expression
n = 0,0429(d10 )0,1371 with a determination coefficient of R2 = 0,956.

Keywords: roughness coefficient; effective diameter; granulometry; Manning’s equation; fluvial hydraulics.

Modelo de estimación del coeficiente de rugosidad de Manning en

función de la granulometría en el río Santa, sector Recuay – Carhuaz,
Áncash, Perú
Resumen.- Dado a que la ecuación de Manning es uno de los métodos de aplicación más utilizados para determinar el
caudal en ríos no aforados, se llevó a cabo una investigación para establecer un modelo para determinar el coeficiente de
rugosidad (n), en un sector montañoso del río Santa en el departamento de Ancash, Perú. Lo novedoso del modelo está en su
simplicidad ya que solo usa como variable de entrada el diámetro efectivo (di), del material de río obtenido del análisis de
granulometría. Se utilizaron tres procedimientos de estimación, medición de velocidades, ecuación de Manning y modelo de
Cowan y para el análisis, modelos de correlación y regresión lineal y no lineal; la bondad de ajuste se efectuó con la prueba t
de Student. Los resultados se validaron con mediciones de campo de las características hidráulicas del tramo en ocho puntos
de muestreo y análisis granulométrico de las muestras de material de fondo. El modelo de mejor ajuste para la estimación del
coeficiente de rugosidad de Manning es el potencial, siendo su expresión matemática n = 0,0429(d10 )0,1371 con un coeficiente
de determinación de R2 = 0,956.

Palabras clave: coeficiente de rugosidad; diámetro efectivo; granulometría; Ecuación de Manning; hidráulica fluvial.

Recibido: 26 de agosto, 2020. disponen de registros de caudales; sin embargo,

Aceptado: 27 de octubre, 2020. esa información es necesaria para el estudio del
régimen hidrométrico y el diseño de las estructuras
hidráulicas para el aprovechamiento del agua.
1. Introducción
Bajo las condiciones expuestas se suelen usar
En los países en vías de desarrollo, como es el procedimientos indirectos desarrollados para zonas
caso del Perú, las estaciones de monitoreo de las carentes de información. La ecuación de Manning
variables hidrometeorológicas están en un estado permite determinar el caudal en ríos no aforados en
incipiente de desarrollo por lo que no siempre se base a los parámetros hidráulicos radio hidráulico,
pendiente hidráulica y área hidráulica; siendo
∗ Autor para correspondencia: la dificultad la determinación del coeficiente de
Correo-e:[email protected] (A. Díaz-Salas) rugosidad n de Manning, el mismo que expresa la
328 Revista Ingeniería UC, ISSN: 1316–6832, Online ISSN: 2610-8240.
UC Universidad
de Carabobo Díaz-Salas et al. / Revista Ingeniería UC, Vol. 27, No 3, Diciembre, 2020 328 – 342

resistencia del cauce al escurrimiento del flujo. se comparan de manera espacial y temporal. En
En el Perú no existen investigaciones para todos ellos los valores obtenidos se adaptan a
la estimación directa de n en los ríos; para las condiciones en las que se llevó a cabo la
los estudios hidráulicos se suele adoptar valores investigación, y al criterio del investigador.
establecidos en condiciones similares en otros Para el estudio de la estimación del coeficiente
países con la desventaja que las similitudes de rugosidad de Manning cuyos resultados se
son relativas, arrojando resultados generales de presentan aquí, se usaron cuatro procedimientos
rangos muy imprecisos, especialmente para ríos de hidráulicos generales: medición de la velocidad;
características muy particulares, como es el caso de medición de la rugosidad; procedimiento de la
los ríos andinos del país, típicos ríos de montaña fórmula de Manning; y Cowan; con el objeto de
con presencia de sedimentos gruesos y pendientes entender mejor el proceso seguido, a continuación,
elevadas. se describe la base teórica utilizada en la
Por lo expuesto, se ha llevado a cabo la presente investigación.
investigación con el objeto de establecer un
modelo de estimación del coeficiente de rugosidad 2.1. Distribución de velocidades en flujo turbu-
de Manning en función de la granulometría en lento
sitios no aforados del río Santa en el sector
Para la aplicación de los métodos de medición
montañoso Recuay–Carhuaz. Lo novedoso del
de velocidad y medición de rugosidad, Osío [1]
modelo está en su sencillez ya que facilita la
indica que para Prandtl, la capa límite turbulenta
determinación de dicho coeficiente solo en función
en el sentido de que la longitud de mezcla (l),
de la granulometría del cauce.
es proporcional a la profundidad del flujo (y), la
fuerza cortante (τ), es constante e igual al esfuerzo
2. Revisión de la literatura cortante unitario (τ0 ), según se muestran en las
ecuaciones (1) y (2).
Existen muchos estudios realizados para obtener
el coeficiente de rugosidad que aplican una serie r
τ0
Vf = gRS =
p
de métodos y/o fórmulas empíricas. (1)
ρ
Investigaciones como la de: Osío [1] indica
que se estimó el coeficiente de rugosidad de y
 
Manning utilizando las funciones de distribución V = 2,5V f ln (2)
y0
de velocidades mediante aforo y que el valor de
n puede determinarse en función del diámetro Donde:
medio de las partículas del material de cauce, V f : velocidad de fricción, [ m/s ]
Fernández [2] menciona que se calcula el n de g: gravedad, [m/s2 ]
Manning por diferentes métodos plasmados en la R: radio hidráulico, [ m ]
literatura y se realiza un análisis de los mismos, S: pendiente hidráulica, [ m/m ]
Burgos [3] indica que se logró obtener el n de τ0 : fricción unitaria en la capa límite, [ kg·m/s ]
Manning mediante una ecuación potencial que está y0 : altura de rugosidad igual a k/30 para superficies
en función del diámetro representativo del lecho rugosas, [ m ]
de río y Pastora [4] desarrolla la evaluación de k: altura de rugosidad, [ m ]
la ecuación de Manning en un flujo de corriente ρ: densidad del fluido, [ kg/m3 ]
natural mediante registro de aforos y determina un V: velocidad del agua, [ m/s ]
coeficiente promedio para época seca y lluviosa. υ: viscosidad cinemática, [ m2 /s]
Existen otros métodos de estimación, Colmenárez y: profundidad de sección [ m ]
[5], usa modelos digitales de rugosidad topográfica,
obtenidos a partir de modelos LIDAR y GPS, De acuerdo con Chow [7], la ecuación (2) se
y Guzmán [6] utiliza valores de caudales que conoce como la Ley Universal de Distribución de
Revista Ingeniería UC, ISSN: 1316–6832, Online ISSN: 2610-8240. 329
UC Universidad
de Carabobo Díaz-Salas et al. / Revista Ingeniería UC, Vol. 27, No 3, Diciembre, 2020 328 – 342

Velocidades de Prandtl – Von Karman, que viene a ecuaciones (7) y (8), se obtiene el n de Manning
ser una función logarítmica. descrita por la ecuación (9).
Para superficies lisas, la constante y0 queda  
R 1/6
expresada por la ecuación (3). n=φ k , (9)
k
υ
y0 = (3) donde:
9V f
(R/k)1/6
 
De donde se obtiene la ecuación (4). R
φ =   (10)
k 12,2 · R

y9V f
 18 log
V = 5, 75 · V f log (4) k
υ
De las ecuaciones (9) y (10), se deduce que el
En superficies rugosas, se encuentra que la valor de k está dado por la ecuación (11).
constante y0 depende de la altura de la rugosidad,
12,2 · R
ésta se indica en la ecuación (5). k= (1/6)
(11)
©R
y0 = mk
ª
(5) ­
18n
®
10« ¬

donde, la constante m es igual aproximadamente donde, k es conocido como la altura efectiva de

a 1/30. Sustituyendo y simplificando para superfi- las irregularidades formados por los elementos de
cies rugosas, se obtiene la ecuación (6). rugosidad.
  Las investigaciones realizadas sobre el tema
30y
V = 5, 75 · V f log (6) indican que para un amplio rango de R/k, la
k variación en φ (R/k) es pequeña, y ésta, puede
dada como la distribución de velocidades en suponerse constante e igual a un valor promedio.
flujos turbulentos sobre superficies rugosas. La Es así que, la ecuación (9) queda determinado en
altura de rugosidad, k, corresponde al diámetro función de la potencia de la altura de rugosidad (k),
medio de los granos de arena utilizado por [7].
Nikuradse, Chow [7]. Según la distribución teórica de velocidades en
Mediante el desarrollo de ecuaciones teóricas canales rugosos, se ha desarrollado los métodos de
para flujo uniforme French [8], al comparar la medición de velocidades y medición de rugosidad,
ecuación de Chezy con la ecuación de Manning [7]. Asimismo, se pueden usar dos métodos
se obtiene la ecuación (7). adicionales que se describen a continuación para
obtener n de Manning.
R1/6
C= (7) 2.2. Método de medición de la velocidad
n
De la ecuación (6), el coeficiente n de Manning
Además, Rocha [9], indica que para conductos
puede determinarse por medio de mediciones de
rugosos se tiene la ecuación (8).
velocidades a 2/10 de y, desde el espejo de agua
12,2R (V0,2 ), dado por la ecuación (12).
C = 18 log (8)
k 
24y

V0,2 = 5,75 · V f · log (12)
donde, C es el coeficiente de Chezy m1/2 /s. k
Siguiendo el análisis de Chow [7] para la
Asimismo, a 8/10 de y, (V0,8 ), expresado por la
interpretación teórica del coeficiente de rugosidad,
ecuación (13).
la relación entre el n de Manning y el valor
 
teórico de la rugosidad del canal para superficies 6y
rugosas, al eliminar el valor de C de Chezy de las V0,8 = 5,75 · V f · log (13)
k
330 Revista Ingeniería UC, ISSN: 1316–6832, Online ISSN: 2610-8240.
UC Universidad
de Carabobo Díaz-Salas et al. / Revista Ingeniería UC, Vol. 27, No 3, Diciembre, 2020 328 – 342

Eliminando V f en las ecuaciones (12) y (13), La granulometría permite calcular mediante

V0,2 el análisis granulométrico por tamizado los
sustituyendo x = , simplificando y resolviendo
V0,8 porcentajes de granos en función de sus tamaños.
se obtiene la ecuación (14). Se designa por dn que es el tamaño, tal que el
n % del peso del material es menor que dicho
(x − 1)y 1/6
n= (14) parámetro.
5,57(x + 0, 95)
2.7. Correlación y regresión
Donde, y representa la profundidad en la sección
2.7.1. Coeficiente de correlación (r)
considerada en metros.
Villón [12], indica que el coeficiente de
2.3. Método de la medición de la rugosidad correlación es el estadístico que permite medir el
grado de asociación de dos variables linealmente
En este método se supone que la ecuación (10)
relacionadas. Sus valores varían entre -1 y 1. Si X
para la función φ(R/k) es aceptable. Luego, el valor
e Y son independientes, r = 0.
de n puede determinarse mediante la ecuación (9)
conocido la altura de rugosidad, k. 2.7.2. Coeficiente de determinación (r 2 )
De acuerdo con Villón [12], r 2 , indica que la
2.4. Método de la fórmula de Manning
proporción o porcentaje de la variación total de la
Conocidos los parámetros involucrados en dicha variable dependiente Y , es explicada por la variable
fórmula, en el punto de interés, se despeja la independiente X; varía en el rango de 0 ≤ r 2 ≤ 1, o
fórmula de Manning para obtener el valor de n, de 0 a 100 %. Cuanto mayor sea el r 2 , mayor será la
que viene a estar dado por la ecuación (15). confiabilidad de las predicciones y menores serán
los errores predictivos. Para la regresión lineal
AR2/3 S 1/2
n= (15) simple se tiene que, r 2 = R2 y pueden ser denotados
Q con letras mayúsculas, significando lo mismo.
Donde, A, es el área hidráulica, R, radio
2.7.3. Regresión lineal simple
hidráulico, S, pendiente media y Q, es el caudal
Según Villón [12], la regresión lineal simple
en la sección transversal de interés.
permite estimar en forma lineal una variable
2.5. Método de Cowan dependiente Y , en función de una variable
independiente X; se suele usar para completar datos
El modelo de Cowan está definido por la
o extender un registro. El modelo de regresión
ecuación (16).
lineal está dado por la ecuación (17).
n = (n0 + n1 + n2 + n3 + n4 ) · m5 (16) y = a+b· x (17)
Los valores se obtienen de la Tabla 1 los cuales Donde, x, es la variable independiente o
están en función del perímetro, irregularidades, conocida, y, es la variable dependiente o que
variación de la sección transversal, efecto de las se trata de predecir, a y b son parámetros que
obstrucciones, vegetación y meandros. son determinados por el método de mínimos
cuadrados.
2.6. Granulometría de material de río
La toma de muestras (García [10]), y la 2.7.4. Regresión no lineal
distribución granulométrica del lecho de río son La regresión no lineal indica que existen una
de gran utilidad, los resultados obtenidos de la asociación no lineal entre la o las variables
granulometría son utilizados en diversos procesos dependientes e independientes. Entre los modelos
fluviales como la estimación de parámetros de no lineales más conocidos se tiene: exponencial,
rugosidad, Crozzoli [11]. potencial y logarítmico. A continuación, se indican
Revista Ingeniería UC, ISSN: 1316–6832, Online ISSN: 2610-8240. 331
UC Universidad
de Carabobo Díaz-Salas et al. / Revista Ingeniería UC, Vol. 27, No 3, Diciembre, 2020 328 – 342

Tabla 1: Criterios y valores numéricos para la ecuación de Cowan [9]

Condición del cauce Valores

Tierra 0,020
Corte en roca 0,025
Material considerado n0
Grava fina 0,024
Grava gruesa 0,028
Suave 0,000
Menor 0,005
Grado de Irregularidad n1
Moderado 0,010
Severo 0,020
Gradual 0,000
Variación de la Sección
Ocasionalmente alternante n2 0,005
Transversal
Frecuentemente alternante 0,010–0,015
Insignificante 0,000
Menor 0,010–0,015
Efecto de las Obstrucciones n3
Apreciable 0,020–0,030
Severo 0,040 - 0,060
Baja 0,005–0,010
Media 0,010–0,025
Vegetación n4
Alta 0,025–0,050
Muy alta 0,050–1,000
Menor 1,000
Cantidad de Meandros Apreciable m5 1,150
Severa 1,300

las ecuaciones del modelo de regresión no lineal: Ambiental, Facultad de Ciencias del Ambiente
Exponencial, está dado por la ecuación (18). de la Universidad Nacional Santiago Antúnez
de Mayolo - UNASAM. Luego se llevó a cabo
y = abx (18) un análisis estadístico relacionando los diámetros
efectivos de granulometría por cada punto de aforo
Potencial, está dado por la ecuación (19). con los valores del coeficiente de rugosidad para
los tres métodos seleccionados.
y = ax b (19)
Descripción del área de estudio
Logarítmica, está dado por la ecuación (20). El área de estudio, se ubica en la región
Ancash, provincias de Recuay, Huaraz y Carhuaz,
y = a ln x + b (20)
y comprende los distritos de Catac, Ticapampa,
Estas ecuaciones pueden ser resueltas por el Recuay, Olleros, Huaraz, Independencia, Jangas,
método de los mínimos cuadrados con un artificio Tarica, Pariahuanca y Anta, ubicada entre las
adecuado pudiéndose reducir a ecuaciones lineales coordenadas UTM 8915258.00 m S, 233116 m
y aplicándose una regresión lineal simple. E aguas abajo del puente Catac – Recuay, y,
8963881.28 m S, 215570.84 m E aguas arriba del
aeropuerto Anta – Carhuaz, entre una altitud que
3. Metodología
varía de 3545 msnm y 2743 msnm respectivamente,
Para obtener el modelo, en primer lugar, se han perteneciente a la Cuenca del Río Santa (Figura 1).
obtenido los valores del coeficiente de rugosidad
3.1. Información básica
por los métodos: distribución de velocidades,
A continuación, se describen resumidamente las
fórmula de Manning y Cowan, con datos obtenidos
fuentes y la información usada en la investigación:
mediante aforos que es una actividad importante
para conocer la disponibilidad del recurso hídrico, 1. Información cartográfica: Cartas Nacionales
Díaz [13] realizados por el Laboratorio de Calidad (19-h, 19-i, 20-h y 20-i), con escala 1:100000,
332 Revista Ingeniería UC, ISSN: 1316–6832, Online ISSN: 2610-8240.
UC Universidad
de Carabobo Díaz-Salas et al. / Revista Ingeniería UC, Vol. 27, No 3, Diciembre, 2020 328 – 342

Figura 1: Ubicación de puntos de aforo, sector Recuay–Carhuaz

y el modelo de elevación digital (DEM) con Ambiental – UNASAM, en periodo de estiaje.

resolución de 30 metros. Los puntos de aforo se muestran en la Figura 1.

2. Información topográfica: Obtenida mediante 3.2. Fases de estudio

uso de equipos topográficos en cada uno de 3.2.1. Fase de campo
los puntos de aforo y sus respectivos tramos. En la fase de campo, se procedió a la recolección
de la información y datos de la zona de estudio
3. Información hidráulica: Obtenida mediante tales como: parámetros hidráulicos de secciones
aforos, Villón [14] por el Servicio de transversales de sitios específicos o puntos de aforo
monitoreo de agua del Laboratorio de Calidad (PA), pendiente, material de lecho y velocidades
Revista Ingeniería UC, ISSN: 1316–6832, Online ISSN: 2610-8240. 333
UC Universidad
de Carabobo Díaz-Salas et al. / Revista Ingeniería UC, Vol. 27, No 3, Diciembre, 2020 328 – 342

de flujo mediante aforo con correntómetro. Así 5. Conocido el A, R, S y Q, mediante la

mismo, fotografías e imágenes satelitales. También ecuación (15) se determinó el coeficiente n
se realizó la verificación y planificación en campo, de Manning.
como el reconocimiento del entorno y su ubicación
de puntos de aforo (PA), seleccionándose ocho 6. De la Tabla 1 y la ecuación (16) se
(08) puntos con su respectivo tramo de estudio. determinaron los valores de n de Manning.
La pendiente se realizó por medición directa
haciendo uso del nivel de ingeniero y la extracción 7. El análisis granulométrico del lecho del río
de material de río se realizó en cada punto de por punto de aforo (PA) se hizo mediante
muestreo. el tamizado siguiendo la norma ASTM D–
422, Gutiérrez [17], obteniéndose así las
curvas granulométricas con sus respectivos
3.2.2. Fase de gabinete diámetros efectivos d10 , d30 , d50 , d60 , d84 , y
En la fase de gabinete se procesó la información, d90 .
estimando el coeficiente n de Manning, los
diámetros efectivos de partículas, di , y estable- 8. Mediante análisis de regresión se estableció
ciendo las relaciones entre ambos parámetros el modelo de estimación del coeficiente
y los modelos de estimación correspondientes. de rugosidad. Se probaron funciones de
Finalmente se validó el modelo establecido de relación, lineal y no lineal aplicando las
modo tal de aceptarlo como representativo para ecuaciones (17), (18), (19) y (20), Sánchez
el área de estudio. A continuación, se resume el [18]. Mediante la prueba t de Student. El t
procedimiento: calculado, tc esta dado por la ecuación (21).
√
r n1 − 2
1. Identificación de los puntos de aforo (PA) y tc = √ (21)
procesamiento de la información de la sección 1 − r2
transversal.
Donde, n1 , es número de pares de valores. Así
mismo, el t tabular, tt se obtiene de las tablas
2. División del ancho de la sección transversal,
estadísticas de la distribución t preparadas
en un número de subsecciones con un mínimo
para este fin con un nivel de significación α o
de un 10 % por cada una de éstas, Villón [14].
una Probabilidad de (1 − α) y con Grado de
Libertad (GL = n1 − 2). Es decir, α/2 = 0,025
3. De las curvas de distribución de velocidades se y GL = n1 − 2. Por último, se aplicó el criterio
obtuvo el valor correspondiente a la velocidad de decisión: si |tc | ≤ tt , entonces r = 0
a dos y a ocho decimos de la profundidad V0,2 por tanto no hay correlación significativa; en
y V0,8 ; es decir, a 0,8 y 0,2 respectivamente cambio, si |tc | > tt , entonces r , 0 por tanto
del fondo del canal, ICC [15] y Tapia [16]. El es significativo y existe correlación entre las
coeficiente n de Manning se calcula en cada variables.
sub sección, en las secciones transversales
de los PA aplicándose la ecuación (14), 9. Finalmente se contrasta y se valida el
conocido el tirante hidráulico y, y la relación modelo con mayor valor de coeficiente de
de velocidades x, el valor de n para la sección correlación propuesto en la estimación del
total se estimó como la media de los valores coeficiente de rugosidad en función de la
de n de las subsecciones. granulometría, con el método de la fórmula de
Manning usando los parámetros hidráulicos
4. El valor de n se calcula en función de la altura medidas directamente en campo, los cuales
de rugosidad k, mediante las ecuaciones (9), se consideran los menos subjetivos de los
(10) y (11). métodos de estimación del n de Manning.
334 Revista Ingeniería UC, ISSN: 1316–6832, Online ISSN: 2610-8240.
UC Universidad
de Carabobo Díaz-Salas et al. / Revista Ingeniería UC, Vol. 27, No 3, Diciembre, 2020 328 – 342

4. Análisis y discusión de resultados

4.1. Descripción de los puntos de aforo (PA)

En la Tabla 2 se presentan la ubicación de
los ocho (8) puntos de aforo seleccionados en
función de una combinación de criterios, tramos
vulnerables y facilidades de medición de los
parámetros hidráulicos.

4.2. Parámetros hidráulicos de secciones trans-

versales
Figura 2: Coeficiente de rugosidad, n, y altura de
En la Tabla 3 se presentan los parámetros rugosidad, k
hidráulicos de las ocho secciones de aforo
seleccionadas (A, área hidráulica, P, perímetro
mojado, R, radio hidráulico, S, pendiente media Dónde: n, desviación estándar de n, R, radio
y Q, caudal) que son aplicados al cálculo del n de hidráulico y k, altura de rugosidades han sido
Manning. obtenidos según la metodología establecida en la
sección 3.2.2.
4.3. Coeficiente de rugosidad de Manning
En la Tabla 4 se presentan los resultados
obtenidos para los n de Manning para cada uno
de los puntos de muestreo y para cada uno de los
procedimientos hidráulicos seleccionados.

4.4. Coeficiente de rugosidad mediante medición

de rugosidad
En la Tabla 5 se presenta los resultados obtenidos
de las alturas de rugosidades, k.
Con dichos resultados se ha realizado un análisis
gráfico entre las variables del coeficiente de Figura 3: Relación de φ(R/k) y R/k.
rugosidad, n, (n calculado mediante medición de
velocidades) y la altura de rugosidad, k, que se La Figura 3, muestra que para un amplio
muestra en la Figura 2. La línea de tendencia a rango de R/k, la variación de φ(R/k) es baja;
la que mejor se ajusta es a una regresión lineal cuyo coeficiente de variación es de 22,90 %,
simple, con un coeficiente de determinación R2 = por consiguiente la relación podría aproximarse
0,9879; lo que indica que el 98,79 % de la variación mediante el valor medio; 0,047 (media de valores
del coeficiente de rugosidad, n, es explicada por de φ(R/k)).
la altura de rugosidad, k, y el 1,21 % restante Combinando la Figura 3 y la ecuación (9) se
es debido a los errores y a otras variables no obtiene n en función de la altura de rugosidad, k,
consideradas. dado por la ecuación (22).
De igual forma, mediante el análisis estadístico,
el tc es igual a 22,23 y el tt es igual a 2,447 lo n = 0,047 · k 1/6 (22)
que significa que el tc es mayor que el tt , y se
concluye que existe correlación significativa entre 4.5. Diámetros representativos del lecho de río
los coeficientes de rugosidad, n y las alturas de De análisis granulométrico se han obtenido
rugosidad, k. los diámetros efectivos di , para las muestras de
Revista Ingeniería UC, ISSN: 1316–6832, Online ISSN: 2610-8240. 335
UC Universidad
de Carabobo Díaz-Salas et al. / Revista Ingeniería UC, Vol. 27, No 3, Diciembre, 2020 328 – 342

Tabla 2: Puntos de aforo (PA) en el río Santa, sector Recuay–Carhuaz

Puntos de aforo Localidad Coordenadas UTM Altitud Progresiva

E S (msnm) (km)
PA-01 Distrito de Catac 233116,00 8915258,00 3545,00 0+000,00
PA-02 Distrito de Ticapampa 232155,00 8920136,00 3460,00 5+166,05
PA-03 Distrito de Recuay 231340,00 8922889,00 3419,00 8+188,40
PA-04 Distrito de Olleros 226914,00 8932961,00 3256,00 20+347,90
PA-05 Distrito de Huaraz 221892,00 8944496,00 3042,00 34+460,25
PA-06 Distrito de Independencia 221386,00 8950581,00 2945,00 41+113,30
PA-07 Distrito de Jangas 219370,00 8957413,00 2850,00 49+158,10
PA-08 Distrito de Pariahuanca 215570,00 8964181,00 2743,00 58+024,35

Tabla 3: Parámetros hidráulicos en cada sección transversal

Puntos de aforo A (m2 ) P (m) R (m) S (m/m) Q (m3 /s)

PA-01 52,420 158,360 0,3310 0,0045 27,724
PA-02 51,460 176,390 0,2917 0,0035 31,179
PA-03 36,950 131,010 0,2820 0,0075 31,018
PA-04 98,450 263,570 0,3735 0,0019 54,290
PA-05 91,300 242,720 0,3762 0,0043 71,346
PA-06 111,650 238,030 0,4691 0,0037 101,387
PA-07 159,280 364,750 0,4367 0,0018 102,944
PA-08 115,100 235,200 0,4894 0,0032 110,651

Tabla 4: Coeficiente de rugosidad, n (medición de velocidades, fórmula de Manning y Cowan)

Método
Puntos de aforo Medición de velocidades Fórmula de Manning Método Cowan
(n) (n) (n)
PA-01 0,0680 0,0607 0,0530
PA-02 0,0421 0,0428 0,0430
PA-03 0,0324 0,0444 0,0480
PA-04 0,0315 0,0410 0,0380
PA-05 0,0345 0,0435 0,0380
PA-06 0,0276 0,0403 0,0480
PA-07 0,0280 0,0377 0,0330
PA-08 0,0237 0,0364 0,0330
Promedio 0,0360 0,0433 0,0418
Desv estándar 0,0141 0,0075 0,0074
Coef, Variación 0,3905 0,1740 0,1782

sedimentos en puntos de aforo, los que se presentan 4.6. Relación de los coeficientes de rugosidad y
en la Tabla 6. los diámetros efectivos

Para analizar las relaciones entre las variables

Los diámetros efectivos d10 , d30 , d50 , d60 , d84 y n de Manning y los diámetros efectivos, di , se
d90 corresponden a 10 %, 30 %, 50 %, 60 %, 84 % generaron los modelos estadísticos descritos en
y 90 % de los porcentajes acumulados que pasan la la metodología, dados por modelos lineales y no
curva granulométrica. lineales (exponencial, potencial y logarítmico).
336 Revista Ingeniería UC, ISSN: 1316–6832, Online ISSN: 2610-8240.
UC Universidad
de Carabobo Díaz-Salas et al. / Revista Ingeniería UC, Vol. 27, No 3, Diciembre, 2020 328 – 342

Tabla 5: Altura de rugosidades k y función φ(R/k)

Puntos de aforo n Des. Est. de n R (m) k (m) R/k φ(R/k)

PA-01 0,0680 0,0292 0,3310 0,8447 0,3919 0,0699
PA-02 0,0421 0,0319 0,2917 0,3003 0,9714 0,0515
PA-03 0,0324 0,0282 0,2820 0,1411 19,991 0,0449
PA-04 0,0315 0,0201 0,3735 0,1447 25,808 0,0434
PA-05 0,0345 0,0278 0,3762 0,1958 19,207 0,0452
PA-06 0,0276 0,0272 0,4691 0,0958 48,963 0,0408
PA-07 0,0280 0,0202 0,4367 0,1002 43,587 0,0411
PA-08 0,0237 0,0245 0,4894 0,0500 97,971 0,0391
Promedio 0,0360 0,0470
Desv. Est. 0,0141 0,0100
Coef. Var. 0,3905 0,2290

Tabla 6: Diámetros efectivos de material de río

Resumen de los diámetros de las partículas (mm)

Puntos de aforo d10 d30 d50 d60 d84 d90
PA - 01 9,651 18,344 25,934 29,681 41,422 54,464
PA - 02 1,060 5,418 13,656 18,762 35,293 44,316
PA - 03 1,586 10,879 17,868 22,615 35,543 42,344
PA - 04 0,659 3,545 8,835 13,759 30,529 35,301
PA - 05 0,904 5,720 15,727 22,436 48,377 58,811
PA - 06 0,939 3,145 6,372 8,433 41,340 54,412
PA - 07 0,491 2,634 8,551 13,137 33,420 44,442
PA - 08 0,202 0,734 3,473 8,300 24,170 29,394

Tabla 7: Análisis estadístico del Modelo Lineal entre el coeficiente de rugosidad de Manning y los
diámetros efectivos d10 , d30 , d50 , d60 , d84 y d90

dd Método
Medición de velocidad (n) Fórmula de Manning (n)
R2 r tc r sig. R2 r tc r sig.
d10 0,887 0,942 6,873 Si 0,931 0,965 8,984 Si
d30 0,764 0,874 4,402 Si 0,921 0,960 8,370 Si
d50 0,762 0,873 4,377 Si 0,841 0,917 5,642 Si
d60 0,680 0,825 3,570 Si 0,739 0,860 4,117 Si
d84 0,165 0,406 1,088 No 0,235 0,485 1,356 No
d90 0,203 0,451 1,237 No 0,259 0,509 1,447 No
r sig.: correlación significativa, tt = 2,447

4.6.1. Regresión lineal 4.6.2. Regresión no lineal

Las estadísticas del modelo lineal entre el n En la Tabla 8, Tabla 9 y Tabla 10, se muestran
de Manning (método de medición de velocidades, los resultados del análisis estadístico entre el n
fórmula de Manning y Cowan) y los diámetros de Manning (método medición de velocidad y
efectivos, di , de la granulometría se muestran en fórmula de Manning) de los modelos: exponencial,
la Tabla 7 y la Tabla 11. potencial y logarítmico respectivamente, con los
diámetros efectivos, di . Asimismo, en la Tabla 11
y Tabla 12 se muestran los resultados del análisis
estadístico entre el n de Manning (método Cowan)

Revista Ingeniería UC, ISSN: 1316–6832, Online ISSN: 2610-8240. 337

UC Universidad
de Carabobo Díaz-Salas et al. / Revista Ingeniería UC, Vol. 27, No 3, Diciembre, 2020 328 – 342

Tabla 8: Análisis estadístico del Modelo Exponencial entre el coeficiente de rugosidad de Manning y los
diámetros efectivos d10 , d30 , d50 , d60 , d84 y d90

dd Método
Medición de velocidad (n) Fórmula de Manning (n)
R2 r tc r sig. R2 r tc r sig.
d10 0,788 0,888 4,717 Si 0,886 0,941 6,842 Si
d30 0,736 0,858 4,093 Si 0,929 0,964 8,840 Si
d50 0,810 0,900 5,053 Si 0,875 0,935 6,484 Si
d60 0,745 0,836 4,188 Si 0,777 0,939 4,567 Si
d84 0,213 0,462 1,272 No 0,272 0,522 1,498 No
d90 0,236 0,486 1,362 No 0,286 0,535 1,548 No
r sig.: correlación significativa, tt = 2,447

Tabla 9: Análisis estadístico del Modelo Potencial entre el coeficiente de rugosidad de Manning y los
diámetros efectivos d10 , d30 , d50 , d60 , d84 y d90

dd Método
Medición de velocidad (n) Fórmula de Manning (n)
R2 r tc r sig. R2 r tc r sig.
d10 0,838 0,915 5,577 Si 0,941 0,970 9,809 Si
d30 0,667 0,817 3,469 Si 0,751 0,867 4,259 Si
d50 0,693 0,832 3,676 Si 0,696 0,834 3,710 Si
d60 0,662 0,814 3,427 Si 0,643 0,802 3,284 Si
d84 0,242 0,492 1,385 No 0,296 0,544 1,588 No
d90 0,253 0,503 1,427 No 0,292 0,540 1,573 No
r sig.: correlación significativa, tt = 2,447

Tabla 10: Análisis estadístico del Modelo Logarítmico entre el coeficiente de rugosidad de Manning y los
diámetros efectivos d10 , d30 , d50 , d60 , d84 y d90

dd Método
Medición de velocidad (n) Fórmula de Manning (n)
R2 r tc r sig. R2 r tc r sig.
d10 0,826 0,909 5,330 Si 0,919 0,959 8,267 Si
d30 0,649 0,806 3,327 Si 0,691 0,831 3,660 Si
d50 0,585 0,765 2,908 Si 0,638 0,799 3,250 Si
d60 0,566 0,752 2,795 Si 0,595 0,771 2,968 Si
d84 0,186 0,431 1,169 No 0,254 0,504 1,429 No
d90 0,212 0,460 1,269 No 0,262 0,512 1,459 No
r sig.: correlación significativa, tt = 2,447

modelo exponencial, y, potencial y logarítmico para la estimación del coeficiente de rugosidad

respectivamente, con los diámetros efectivos di . de Manning, n, por cada método desarrollado
(método de medición de velocidades, ecuación de
4.7. Modelo de estimación del coeficiente de Manning y Cowan). Los resultados se presentan
rugosidad en la Tabla 13 resaltados con un sombreado.
Se seleccionaron los modelos (lineal y no Observando los valores de dicha tabla, se encuentra
lineal) que muestran mejor grado de asociación que para estimar el n de Manning por el método
y, por tanto, mayor coeficiente de correlación, r, de la medición de velocidades, el modelo lineal

338 Revista Ingeniería UC, ISSN: 1316–6832, Online ISSN: 2610-8240.

UC Universidad
de Carabobo Díaz-Salas et al. / Revista Ingeniería UC, Vol. 27, No 3, Diciembre, 2020 328 – 342

Tabla 11: Análisis estadístico del Modelo Lineal y Exponencial entre el coeficiente de rugosidad de
Manning y los diámetros efectivos d10 , d30 , d50 , d60 , d84 y d90

dd Modelo
Lineal. Método de Cowan (n) Exponencial. Método de Cowan (n)
R2 r tc r sig. R2 r tc r sig.
d10 0,483 0,695 2,366 Si 0,434 0,659 2,147 Si
d30 0,630 0,794 3,199 Si 0,598 0,773 2,985 Si
d50 0,484 0,696 2,371 No 0,471 0,686 2,310 No
d60 0,336 0,580 1,744 No 0,330 0,574 1,720 No
d84 0,242 0,492 1,385 No 0,257 0,507 1,442 No
d90 0,272 0,522 1,495 No 0,277 0,526 1,516 No
r sig.: correlación significativa, tt = 2,447

Tabla 12: Análisis estadístico del Modelo Potencial y Logarítmica entre el coeficiente de rugosidad de
Manning y los diámetros efectivos d10 , d30 , d50 , d60 , d84 y d90

dd Modelo
Potencial. Método de Cowan (n) Logarítmico. Método de Cowan (n)
R2 r tc r sig. R2 r tc r sig.
d10 0,735 0,857 4,076 Si 0,756 0,869 4,315 Si
d30 0,629 0,793 3,187 Si 0,621 0,788 3,138 Si
d50 0,415 0,644 2,062 No 0,408 0,639 2,031 No
d60 0,247 0,497 1,403 No 0,245 0,495 1,394 No
d84 0,309 0,556 1,638 No 0,292 0,540 1,572 No
d90 0,305 0,552 1,621 No 0,297 0,545 1,593 No
r sig.: correlación significativa, tt = 2,447

arroja el mejor ajuste, con un coeficiente de Y, por último, para el método de la medición de
correlación r = 0,942 entre las variables n y el la rugosidad se ha obtenido un modelo Potencial
d10 ; estableciéndose así el modelo lineal (MLin) (MPot-01) dado en la ecuación (26)
dada en la ecuación (23).
n = 0,047(k)1/6 (26)
n = 0,0042(d10 ) + 0,0278 (23)
Cabe señalar que el modelo MPot-01 está en
Para el método de fórmula de Manning, el función de la altura de rugosidad k, en metros,
modelo Potencial posee mejor ajuste entre las pero relacionándolo mediante análisis estadístico
variables n y el d10 con un coeficiente de con los diámetros efectivos di de la granulometría
correlación, r = 0,970. De esta forma se y contrastando con n obtenidos por fórmula de
obtiene el modelo potencial (MPot-02) dada en la Manning, se ha encontrado mejor ajuste con el d10
ecuación (24). (mm), es así que se puede reemplazar el k, por el
d10 obteniendo así el modelo Potencial (MPot-01)
n = 0,0429(d10 )0,1371 (24)
dado en la ecuación (27).
Para el método de Cowan, el modelo Logarít-
mico presenta mejor ajuste con un coeficiente de n = 0,047(d10 )1/6 (27)
correlación,r = 0,869 entre las variables n y el d10 .
4.8. Validación del modelo general
De esta forma se obtiene el modelo logarítmico
(MLog) dada por la ecuación (25). Para validar el modelo general de la inves-
tigación, se estiman los n de Manning por los
n = 0,0058 ln(d10 ) + 0,0418 (25) cuatro procedimientos desarrollados y expresados
Revista Ingeniería UC, ISSN: 1316–6832, Online ISSN: 2610-8240. 339
UC Universidad
de Carabobo Díaz-Salas et al. / Revista Ingeniería UC, Vol. 27, No 3, Diciembre, 2020 328 – 342

Tabla 13: Modelos de regresión para cada método de estimación desarrollado

Modelo
Med. de Vel. (n) Manning (n) Cowan (n)
Regresión R2 di r r sig. R2 di r r sig. R2 di r r sig.
Lineal 0,887 d10 0,942 Si 0,930 d10 0,964 Si 0,630 d30 0,794 Si
Exponencial 0,810 d50 0,900 Si 0,928 d30 0,963 Si 0,598 d30 0,773 Si
Potencial 0,838 d10 0,915 Si 0,941 d10 0,970 Si 0,735 d10 0,857 Si
Logarítmica 0,826 d10 0,909 Si 0,919 d10 0,959 Si 0,756 d10 0,869 Si

mediante las ecuaciones (23), (24), (25) y (27) Manning se encuentra dentro del rango establecido
y se contrasta con el n obtenido por el método de 0,030 a 0,050 para corrientes de montañas
de la fórmula de Manning (n observados). En la con fondo: grava, piedras y pocos cantos rodados
Tabla 14, se presentan los resultados numéricos propuesto por Chow [7].
comparativos de este análisis y en la Figura 3 se
muestra los resultados gráficos de los coeficientes La función obtenida para estimar n de Manning
de correlación, r, de n de Manning observados por el método de la medición de la rugosidad, está
y estimados por los modelos encontrados. Así dado por: n = 0,047(k)(1/6) que viene a ser mayor
mismo de la Tabla 14, MLin, MPot-01, MPot-02 y a la ecuación de n = 0,0336(k)(1/6) obtenido en
MLog son abreviaturas del: modelo lineal, modelo los grandes ríos de Venezuela [1], y análogo a
potencial 01, modelo potencial 02 y modelo la ecuación n = 0,044(k)(1/6) obtenido en el río
logarítmico respectivamente. Chonta, distrito de Baños del Inca – Cajamarca
[3].
4.9. Comparativo entre los coeficientes de rugo-
sidad de Manning Se observa que en la aplicación del análisis de
regresión lineal y no lineal entre el coeficiente de
La Tabla 15 muestra una comparación de los
rugosidad n de Manning y los diámetros efectivos
valores de n de Manning estimados por los métodos
di , en la Tabla 7, Tabla 8, Tabla 9, Tabla 10, Tabla 11
de estimación (medición de velocidades, medición
y Tabla 12, para los diámetros menores o iguales al
de rugosidades, fórmula de Manning y Cowan).
diámetro medio d50 se obtienen mejor coeficiente
Donde se observa que el promedio de n de
de correlación, r, a diferencia de los diámetros
los cuatro métodos es 0,042 y el coeficiente de
mayores que no guardan relación entre las dos
variación es 11,37 %, por esta razón se asume el
variables, López [19].
valor de n como representativo para el río Santa en
la zona de estudio. Del contraste de los modelos, se observa que
son dos los modelos de mejor ajuste: MPot-02,
4.10. Discusiones
con r = 0,9776 y MPot-01 con r = 0,9804. La
A partir de los resultados obtenidos, de la estimación de n de Manning mediante el modelo
Tabla 15, se ha determinado un nmín = 0,024, MPot-02 da como resultado un npromedio = 0,043
nmáx = 0,069 y npromedio = 0,042 (considerado y la estimación mediante el modelo MPot-01 da
como la media de los métodos desarrollados: como resultado un npromedio = 0,048 (Tabla 14).
medición de velocidades, medición de rugosidad,
fórmula de Manning y Cowan). Se puede decir que el n medio del MPot-02
El método de medición de la rugosidad no se aproxima al n estimado como la media de los
establece directamente el valor de n, sino que métodos de medición de velocidades, fórmula de
establece una función para calcular n en función Manning y Cowan, y el n medio de MPot-01
de la altura de rugosidad, k. Es así que, entre diverge del n estimado como media de los métodos
los métodos antes mencionados, el valor de n de desarrollados.
340 Revista Ingeniería UC, ISSN: 1316–6832, Online ISSN: 2610-8240.
UC Universidad
de Carabobo Díaz-Salas et al. / Revista Ingeniería UC, Vol. 27, No 3, Diciembre, 2020 328 – 342

Tabla 14: Valores del coeficiente de rugosidad de Manning (n) observados y estimados

PA Observados Estimados
n Manning n (MLin) n (MPot-01) n (MPot-02) n (MLog)
PA-01 0,0607 0,0683 0,0686 0,0585 0,0549
PA-02 0,0428 0,0322 0,0475 0,0432 0,0421
PA-03 0,0444 0,0345 0,0508 0,0457 0,0445
PA-04 0,0410 0,0306 0,0438 0,0405 0,0394
PA-05 0,0435 0,0316 0,0462 0,0423 0,0412
PA-06 0,0403 0,0317 0,0465 0,0425 0,0414
PA-07 0,0377 0,0299 0,0418 0,0389 0,0377
PA-08 0,0364 0,0286 0,0360 0,0345 0,0325
Promedio 0,0433 0,0359 0,0476 0,0433 0,0417
Desv. Est. 0,0075 0,0132 0,0095 0,0070 0,0064
Coef Var. 0,1740 0,3676 0,2004 0,1617 0,1540

Tabla 15: Valores del coeficiente de rugosidad n de Manning

Puntos de aforo Medic. de Vel. Manning Cowan Medic. de Rug.

PA-01 0,0680 0,0607 0,0530 0,0686
PA-02 0,0421 0,0428 0,0430 0,0475
PA-03 0,0324 0,0444 0,0480 0,0508
PA-04 0,0315 0,0410 0,0380 0,0438
PA-05 0,0345 0,0435 0,0380 0,0462
PA-06 0,0276 0,0403 0,0480 0,0465
PA-07 0,0280 0,0377 0,0330 0,0418
PA-08 0,0237 0,0364 0,0330 0,0360
Promedio 0,0360 0,0433 0,0418 0,0476
Desv. Est. 0,0141 0,0075 0,0074 0,0095
Coef. Var. 0,3904 0,1740 0,1780 0,2004

5. Conclusiones 6. Recomendaciones

Realizar estudios en épocas de lluvias y en más

El coeficiente de rugosidad de Manning (n) secciones del río Santa y en otros ríos del país.
estimado por el método de medición de velocidades Extraer por cada punto de aforo cinco muestras del
es lineal en función del d10 . material de lecho de río.
El coeficiente de rugosidad de Manning estima- Se recomienda sistematizar los cálculos de
do por el método de la fórmula de Manning y por el estimación del n de Manning para su mejor
método de la medición de la rugosidad es potencial aplicación, así como implementar redes de
en función del d10 . monitoreo hidrometeorológicos e instalación de
estaciones.
El coeficiente de rugosidad de Manning (n)
estimado por el método de Cowan es logarítmico
en función del d10 . 7. Referencias
El modelo que presenta mayor coeficiente de [1] M. Osío, F. Valencia, E. Guevara, y H. Cartaya di Lena,
correlación es el modelo potencial y el que presenta “Cálculo del coeficiente de rugosidad “n” de manning
menor coeficiente de correlación es el modelo en los grandes ríos de Venezuela,” Revista Ingeniería
logarítmico. UC, vol. 7, no. 2, pp. 1–2, 2000.
[2] C. Fernández, A. León, y Y. Rodríguez, “Revista
El coeficiente de rugosidad de Manning adecua- INGENIERÍA HIDRÁULICA Y AMBIENTAL,”
do para la zona de estudio es de 0,042. Influencia del método de estimación en el coeficiente

Revista Ingeniería UC, ISSN: 1316–6832, Online ISSN: 2610-8240. 341

 UC Universidad
de Carabobo Díaz-Salas et al. / Revista Ingeniería UC, Vol. 27, No 3, Diciembre, 2020 328 – 342

de Manning para cauces naturales, vol. 39, no. 1, pp. Salle, 2009.
17–31, 2018. [18] R. Sánchez, “t-Student. Usos y abusos,” Revista
[3] N. Burgos, “Estimación del coeficiente de rugosidad Mexicana de Cardiología, vol. 26, no. 1, pp. 59–61,
de Manning mediante mediciones de velocidad y 2015.
profundidad, empleando un molinete hidrométrico, en [19] R. López, J. Barragán, y A. Colomer, “Predicción de
el río Chonta, Cajamarca, 2016,” Tesis de Maestría, la resistencia al flujo en ríos de montaña,” Revista
Universidad Nacional de Cajamarca, 2017. Ingeniería hidráulica en México, vol. 23, no. 3, pp.
[4] D. Pastora, “Evaluación de la fórmula de Manning en 65–76, 2008.
el río Ostua,” Tesis de Maestría, Universidad de San
Carlos de Guatemala, 2010.
[5] G. Colmenárez and J. Pardo, “Estudio de la relación
de la rugosidad topográfica obtenida a partir de datos
LIDAR y GPS con el coeficiente de rugosidad n de
Manning,” Revista C& G, vol. 24, no. 1–2, pp. 143–
159, 2010.
[6] G. Guzmán y T. Vera, “Estudio del intercambio de agua
superficial y agua subterranea en el río Cumbe,” Tesis
de Maestría, Universidad de Cuenca, 2013.
[7] V. Chow, Hidráulica de Canales Abiertos. MacGraw
Hill Latinoamericana, 2004.
[8] R. French, Hidráulica de Canales Abiertos. McGraw-
Hill Latinoamericana, 1988.
[9] A. Rocha, Tuberías y Canalizaciones. Perú:
Universidad Nacional de Ingeniería, 2007.
[10] C. García y J. Martín, “Caracterización granulométrica
del lecho móvil de un río de gravas efímero: Aplicación
a un tramo de la riera de Les Arenes,” Revista Acta
Geológica Hispánica, vol. 36, no. 1–2, pp. 137–147,
2001.
[11] L. Crozzoli y R. Batalla, “Aplicación de la fotografía
al análisis granulométrico de ríos con lecho de gravas,”
Revista C& G, vol. 17, no. 3–4, pp. 29–39, 2003.
[12] M. Villón, Hidrología Estadística. Instituto
Tecnológico de Costa Rica. Escuela de Ingeniería
Agrícola, 2005.
[13] C. Díaz, J. Segovia, M. Garduño, y S. Tejeda, “Medición
de caudales mediante la implementación de un
vehículo acuático teleoperado,” Revista Internacional
de Contaminación Ambiental, vol. 28, no. 1, pp. 73–91,
2012.
[14] M. Villón, Hidrología. Editorial Tecnológica de Costa
Rica, 2004.
[15] Instituto Privado de Investigación sobre el cambio
climático, Manual de medición de caudales. Instituto
Privado de Investigación sobre el cambio climático,
2017.
[16] G. Tapia, J. Molina, G. Pérez, y A. Torres, “Metodología
para la medición de la velocidad de flujo en un río en el
diagnóstico de la socavación en pilas de un puente,
utilizando un dispositivo electrónico,” Universidad
Michoacana de San Nicolás de Hidalgo, Reporte
Técnico, 2012.
[17] J. Gutiérrez, “Cálculo del coeficiente de rugosidad de
manning utilizando gravilla, arena, piedra pegada y
tierra como fondo mediante un canal a escala como
modelo físico,” Tesis de Maestría, Universidad de la

342 Revista Ingeniería UC, ISSN: 1316–6832, Online ISSN: 2610-8240.

View publication stats