stringtranslate.com

Transporte de sedimentos

Polvo que sopla desde el desierto del Sahara sobre el Océano Atlántico hacia las Islas Canarias

El transporte de sedimentos es el movimiento de partículas sólidas ( sedimento ), generalmente debido a una combinación de la gravedad que actúa sobre el sedimento y el movimiento del fluido en el que se arrastra el sedimento. El transporte de sedimentos ocurre en sistemas naturales donde las partículas son rocas clásticas ( arena , grava , cantos rodados , etc.), lodo o arcilla ; el fluido es aire, agua o hielo; y la fuerza de la gravedad actúa para mover las partículas a lo largo de la superficie inclinada sobre la que descansan. El transporte de sedimentos debido al movimiento de fluidos ocurre en ríos , océanos , lagos , mares y otros cuerpos de agua debido a corrientes y mareas . El transporte también es causado por los glaciares a medida que fluyen y en superficies terrestres bajo la influencia del viento . El transporte de sedimentos debido solo a la gravedad puede ocurrir en superficies inclinadas en general, incluidas laderas de colinas , escarpes , acantilados y la plataforma continental —límite del talud continental.

El transporte de sedimentos es importante en los campos de la geología sedimentaria , la geomorfología , la ingeniería civil , la ingeniería hidráulica y la ingeniería ambiental (ver aplicaciones, a continuación). El conocimiento del transporte de sedimentos se utiliza con mayor frecuencia para determinar si se producirá erosión o deposición , la magnitud de esta erosión o deposición y el tiempo y la distancia en los que se producirá.

Entornos

eólico

Arena que se desprende de una cresta en las dunas de Kelso del desierto de Mojave , California

El término eólico o eólico (según el análisis de æ ) se utiliza para designar el transporte de sedimentos por el viento . Este proceso da lugar a la formación de ondulaciones y dunas de arena . Normalmente, el tamaño del sedimento transportado es arena fina (<1 mm) y menor, porque el aire es un fluido con baja densidad y viscosidad y, por lo tanto, no puede ejercer mucha fuerza de cizallamiento sobre su lecho.

Las formaciones de lecho se generan por el transporte de sedimentos eólicos en el entorno terrestre cercano a la superficie. Las ondulaciones [1] y las dunas [2] se forman como una respuesta natural de autoorganización al transporte de sedimentos.

El transporte de sedimentos eólicos es común en las playas y en las regiones áridas del mundo, porque es en estos ambientes donde la vegetación no impide la presencia y el movimiento de campos de arena.

El polvo de grano muy fino arrastrado por el viento es capaz de entrar en la atmósfera superior y moverse por todo el globo. El polvo del Sahara se deposita en las Islas Canarias y en las islas del Caribe [3] , y el polvo del desierto de Gobi se ha depositado en el oeste de los Estados Unidos [4] . Este sedimento es importante para el presupuesto del suelo y la ecología de varias islas.

Los depósitos de sedimentos glaciares de grano fino arrastrados por el viento se denominan loess .

Fluvial

Depósitos glaciofluviales profundos y erosionados a lo largo del río Matanuska , Alaska

En geografía y geología , los procesos de sedimentación fluvial o el transporte de sedimentos fluviales están asociados con ríos y arroyos y los depósitos y accidentes geográficos creados por sedimentos . Puede dar lugar a la formación de ondulaciones y dunas , a patrones de erosión en forma de fractales , a patrones complejos de sistemas fluviales naturales y al desarrollo de llanuras aluviales y la aparición de inundaciones repentinas . Los sedimentos transportados por el agua pueden ser más grandes que los sedimentos transportados por el aire porque el agua tiene una mayor densidad y viscosidad . En los ríos típicos, el sedimento transportado más grande es del tamaño de arena y grava , pero las inundaciones más grandes pueden arrastrar guijarros e incluso rocas .

Cuando el arroyo o los ríos están asociados con glaciares , capas de hielo o casquetes polares , se utiliza el término glaciofluvial o fluvioglacial , como en los flujos periglaciales y las inundaciones por desbordamiento de lagos glaciares . [5] [6] Los procesos de sedimentación fluvial incluyen el movimiento de sedimentos y la erosión o deposición en el lecho del río . [7] [8]

Costero

Ondulaciones en la arena , playa Laysan, Hawái . El transporte de sedimentos costeros produce estas ondulaciones uniformemente espaciadas a lo largo de la orilla. Foca monje como referencia.

El transporte de sedimentos costeros se produce en entornos cercanos a la costa debido al movimiento de las olas y las corrientes. En las desembocaduras de los ríos, los procesos de transporte de sedimentos costeros y fluviales se combinan para crear deltas fluviales .

El transporte de sedimentos costeros da lugar a la formación de accidentes geográficos costeros característicos, como playas , islas barrera y cabos. [9]

Un glaciar se une al glaciar Gorner , Zermatt, Suiza . Estos glaciares transportan sedimentos y dejan tras de sí morrenas laterales .

Glacial

A medida que los glaciares se desplazan sobre sus lechos, arrastran y mueven material de todos los tamaños. Los glaciares pueden transportar los sedimentos más grandes, y las áreas de deposición glaciar a menudo contienen una gran cantidad de rocas erráticas glaciares , muchas de las cuales tienen varios metros de diámetro. Los glaciares también pulverizan la roca en " harina glaciar ", que es tan fina que a menudo es arrastrada por los vientos para crear depósitos de loess a miles de kilómetros de distancia. El sedimento arrastrado en los glaciares a menudo se mueve aproximadamente a lo largo de las líneas de flujo glaciar , lo que hace que aparezca en la superficie en la zona de ablación .

Ladera

En el transporte de sedimentos en laderas, una variedad de procesos desplazan el regolito pendiente abajo. Entre ellos se incluyen:

Estos procesos generalmente se combinan para dar a la ladera un perfil que parece una solución a la ecuación de difusión , donde la difusividad es un parámetro que se relaciona con la facilidad de transporte de sedimentos en la ladera en particular. Por esta razón, las cimas de las colinas generalmente tienen un perfil parabólico cóncavo hacia arriba, que se va transformando en un perfil convexo hacia arriba alrededor de los valles.

Sin embargo, a medida que las laderas se hacen más empinadas, se vuelven más propensas a deslizamientos episódicos y otros eventos de erosión masiva . Por lo tanto, los procesos en laderas se describen mejor mediante una ecuación de difusión no lineal en la que la difusión clásica domina para las laderas poco profundas y las tasas de erosión tienden al infinito a medida que la ladera alcanza un ángulo crítico de reposo . [10]

Flujo de escombros

Grandes masas de material se desplazan en corrientes de detritos , mezclas hiperconcentradas de lodo, clastos que pueden alcanzar el tamaño de una roca y agua. Las corrientes de detritos se desplazan como flujos granulares por valles y cauces montañosos empinados. Debido a que transportan sedimentos como una mezcla granular, sus mecanismos y capacidades de transporte varían de manera diferente a los de los sistemas fluviales.

Aplicaciones

Sedimentos suspendidos de un arroyo que desemboca en un fiordo ( Isfjorden , Svalbard, Noruega)

El transporte de sedimentos se aplica para resolver muchos problemas ambientales, geotécnicos y geológicos. Por lo tanto, la medición o cuantificación del transporte de sedimentos o la erosión es importante para la ingeniería costera . Se han diseñado varios dispositivos de erosión de sedimentos para cuantificar la erosión de sedimentos (por ejemplo, el Simulador de Erosión de Partículas (PES)). Uno de estos dispositivos, también conocido como BEAST (Herramienta de Evaluación Ambiental Bentónica de Sedimentos), se ha calibrado para cuantificar las tasas de erosión de sedimentos. [11]

El movimiento de sedimentos es importante para proporcionar hábitat a los peces y otros organismos en los ríos. Por lo tanto, a los administradores de ríos muy regulados, que a menudo carecen de sedimentos debido a las represas, se les suele aconsejar que realicen breves inundaciones para renovar el material del lecho y reconstruir las barras. Esto también es importante, por ejemplo, en el Gran Cañón del río Colorado , para reconstruir los hábitats costeros que también se utilizan como lugares de acampada.

La descarga de sedimentos en un embalse formado por una presa forma un delta del embalse . Este delta llenará la cuenca y, con el tiempo, será necesario dragar el embalse o retirar la presa. El conocimiento del transporte de sedimentos puede utilizarse para planificar adecuadamente la prolongación de la vida útil de una presa.

Los geólogos pueden utilizar soluciones inversas de relaciones de transporte para comprender la profundidad, la velocidad y la dirección del flujo de rocas sedimentarias y depósitos jóvenes de materiales aluviales.

El flujo en alcantarillas, sobre represas y alrededor de pilares de puentes puede provocar erosión del lecho. Esta erosión puede dañar el medio ambiente y exponer o desestabilizar los cimientos de la estructura. Por lo tanto, es importante que los ingenieros civiles e hidráulicos tengan un buen conocimiento de la mecánica del transporte de sedimentos en un entorno construido.

Cuando el transporte de sedimentos en suspensión aumenta debido a las actividades humanas, causando problemas ambientales, incluido el llenado de canales, se denomina sedimentación debido a la fracción de tamaño de grano que domina el proceso.

Iniciación del movimiento

Mecanismos de movimiento en el transporte de sedimentos. (a) Rodar: la partícula de sedimento gira por el esfuerzo cortante a lo largo del suelo. (b) Elevación: la partícula de sedimento se eleva por el esfuerzo cortante hacia el interior del volumen. (c) Arrancar: la partícula de sedimento se extrae de las grietas del suelo por el esfuerzo cortante.

Equilibrio del estrés

Para que un fluido comience a transportar sedimentos que se encuentran en reposo sobre una superficie, la tensión de corte límite (o de lecho) ejercida por el fluido debe superar la tensión de corte crítica para el inicio del movimiento de los granos en el lecho. Este criterio básico para el inicio del movimiento se puede escribir como:

.

Esto se representa típicamente mediante una comparación entre una tensión cortante adimensional y una tensión cortante crítica adimensional . La adimensionalización se realiza para comparar las fuerzas impulsoras del movimiento de partículas (tensión cortante) con las fuerzas de resistencia que lo harían estacionario (densidad y tamaño de partículas). Esta tensión cortante adimensional, , se denomina parámetro de Shields y se define como: [12]

.

Y la nueva ecuación a resolver queda:

.

Las ecuaciones incluidas aquí describen el transporte de sedimentos clásticos o granulares . No funcionan para arcillas y lodos porque estos tipos de sedimentos floculares no se ajustan a las simplificaciones geométricas de estas ecuaciones y también interactúan a través de fuerzas electrostáticas . Las ecuaciones también fueron diseñadas para el transporte de sedimentos fluviales de partículas transportadas en un flujo de líquido, como el de un río, canal u otro canal abierto.

En esta ecuación se considera un solo tamaño de partícula. Sin embargo, los lechos de los ríos suelen estar formados por una mezcla de sedimentos de varios tamaños. En caso de movimiento parcial, en el que solo se mueve una parte de la mezcla de sedimentos, el lecho del río se enriquece con grava grande a medida que los sedimentos más pequeños son arrastrados. Los sedimentos más pequeños presentes debajo de esta capa de grava grande tienen una menor posibilidad de movimiento y el transporte total de sedimentos disminuye. Esto se denomina efecto de blindaje. [13] Otras formas de blindaje de sedimentos o tasas decrecientes de erosión de sedimentos pueden ser causadas por alfombras de esteras microbianas, en condiciones de alta carga orgánica. [14]

Esfuerzo cortante crítico

Diagrama original de Shields, 1936

El diagrama de Shields muestra empíricamente cómo la tensión cortante crítica adimensional (es decir, la tensión cortante adimensional requerida para el inicio del movimiento) es una función de una forma particular del número de Reynolds de la partícula , o del número de Reynolds relacionado con la partícula. Esto permite reescribir el criterio para el inicio del movimiento en términos de una solución para una versión específica del número de Reynolds de la partícula, denominada .

Esto puede resolverse utilizando la curva de Shields derivada empíricamente para encontrar como función de una forma específica del número de Reynolds de la partícula llamado número de Reynolds límite. La solución matemática de la ecuación fue dada por Dey . [15]

Número de Reynolds de la partícula

En general, un número de Reynolds de partícula tiene la forma:

Donde es una velocidad de partícula característica, es el diámetro de grano (un tamaño de partícula característico), y es la viscosidad cinemática, que viene dada por la viscosidad dinámica, , dividida por la densidad del fluido, .

El número de Reynolds de la partícula específica de interés se denomina número de Reynolds límite y se forma reemplazando el término de velocidad en el número de Reynolds de la partícula por la velocidad de corte , , que es una forma de reescribir la tensión de corte en términos de velocidad.

donde es la tensión cortante del lecho (descrita a continuación), y es la constante de von Kármán , donde

.

Por lo tanto, el número de Reynolds de la partícula viene dado por:

Esfuerzo cortante del lecho

El número de Reynolds límite se puede utilizar con el diagrama de Shields para resolver empíricamente la ecuación.

,

que resuelve el lado derecho de la ecuación

.

Para resolver el lado izquierdo, expandido como

,

Es necesario encontrar la tensión cortante del lecho . Hay varias formas de calcular la tensión cortante del lecho. El enfoque más simple es suponer que el flujo es constante y uniforme, utilizando la profundidad y la pendiente promediadas en el tramo. Debido a que es difícil medir la tensión cortante in situ , este método también es uno de los más utilizados. El método se conoce como el producto profundidad-pendiente .

Producto profundidad-pendiente

Para un río que experimenta un flujo de equilibrio uniforme y aproximadamente constante, de profundidad h y ángulo de pendiente θ aproximadamente constantes en el tramo de interés, y cuyo ancho es mucho mayor que su profundidad, la tensión cortante del lecho está dada por algunas consideraciones de momento que establecen que el componente de fuerza de gravedad en la dirección del flujo es exactamente igual a la fuerza de fricción. [16] Para un canal ancho, se obtiene:

Para los ángulos de pendiente poco profundos, que se encuentran en casi todos los arroyos naturales de tierras bajas, la fórmula de ángulo pequeño muestra que es aproximadamente igual a , que viene dada por , la pendiente. Reescrito de esta manera:

Velocidad de corte, velocidad y factor de fricción

Para el caso estable, extrapolando el producto profundidad-pendiente y la ecuación para la velocidad de corte:

,

El producto profundidad-pendiente se puede reescribir como:

.

está relacionada con la velocidad media del flujo, , a través del factor de fricción generalizado de Darcy-Weisbach , , que es igual al factor de fricción de Darcy-Weisbach dividido por 8 (por conveniencia matemática). [17] Insertando este factor de fricción,

.

Flujo inestable

Para todos los flujos que no pueden simplificarse como un canal infinito de pendiente única (como en el producto profundidad-pendiente , arriba), la tensión cortante del lecho se puede encontrar localmente aplicando las ecuaciones de Saint-Venant para continuidad , que consideran aceleraciones dentro del flujo.

Ejemplo

Configuración

El criterio de iniciación del movimiento, establecido anteriormente, establece que

.

En esta ecuación,

, y por lo tanto
.
es una función del número de Reynolds límite, un tipo específico de número de Reynolds de partícula.
.

Para una partícula particular, el número de Reynolds será una constante empírica dada por la curva de Shields o por otro conjunto de datos empíricos (dependiendo de si el tamaño del grano es uniforme o no).

Por lo tanto, la ecuación final a resolver es:

.

Solución

Algunas suposiciones permiten la solución de la ecuación anterior.

La primera suposición es que una buena aproximación de la tensión cortante promediada en el tramo se obtiene mediante el producto profundidad-pendiente. La ecuación puede entonces reescribirse como:

.

Mover y volver a combinar los términos produce:

donde R es la gravedad específica sumergida del sedimento.

El segundo supuesto es que el número de Reynolds de las partículas es alto. Esto se aplica típicamente a partículas del tamaño de la grava o más grandes en un arroyo, y significa que la tensión de corte crítica es constante. La curva de Shields muestra que para un lecho con un tamaño de grano uniforme,

.

Investigadores posteriores [18] han demostrado que este valor está más cerca de

para lechos clasificados de manera más uniforme. Por lo tanto, el reemplazo

se utiliza para insertar ambos valores al final.

La ecuación ahora se lee:

Esta expresión final muestra que el producto de la profundidad del canal y la pendiente es igual al criterio de Shield multiplicado por la gravedad específica sumergida de las partículas multiplicada por el diámetro de las partículas.

Para una situación típica, como un sedimento rico en cuarzo en el agua , la gravedad específica sumergida es igual a 1,65.

Conectando esto a la ecuación anterior,

.

Para el criterio de Shield de 0,06 * 1,65 = 0,099, que se encuentra dentro de los márgenes de error estándar de 0,1. Por lo tanto, para un lecho uniforme,

.

Para estas situaciones, el producto de la profundidad por la pendiente del flujo debe ser el 10% del diámetro del grano medio.

El valor del lecho de tamaño de grano mixto es , lo que está respaldado por investigaciones más recientes como de aplicación más amplia porque la mayoría de los arroyos naturales tienen tamaños de grano mixtos. [18] Si se utiliza este valor y D se cambia a D_50 ("50" para el percentil 50, o el tamaño de grano medio, como un valor apropiado para un lecho de tamaño de grano mixto), la ecuación se convierte en:

Lo que significa que la profundidad multiplicada por la pendiente debe ser aproximadamente el 5% del diámetro de grano medio en el caso de un lecho de tamaño de grano mixto.

Modos de arrastre

Los sedimentos arrastrados en un flujo pueden transportarse a lo largo del lecho como carga de fondo en forma de granos deslizantes y rodantes, o en suspensión como carga suspendida transportada por el flujo principal. [16] Algunos materiales de sedimentos también pueden provenir de los tramos superiores y ser transportados río abajo en forma de carga de lavado .

Número de Rouse

La ubicación en el flujo en la que se arrastra una partícula está determinada por el número de Rouse , que está determinado por la densidad ρ s y el diámetro d de la partícula de sedimento, y la densidad ρ y la viscosidad cinemática ν del fluido determinan en qué parte del flujo será arrastrada la partícula de sedimento. [19]

Aquí, el número de Rouse se da por P . El término en el numerador es la velocidad de sedimentación del sedimento (hacia abajo) w s , que se analiza a continuación. La velocidad hacia arriba en el grano se da como un producto de la constante de von Kármán , κ  = 0,4, y la velocidad de corte , u .

La siguiente tabla proporciona los números de Rouse requeridos aproximados para el transporte como carga de fondo , carga suspendida y carga de lavado . [19] [20]

Velocidad de asentamiento

Líneas de corriente alrededor de una esfera que cae a través de un fluido. Esta ilustración es precisa para el flujo laminar , en el que el número de Reynolds de la partícula es pequeño. Esto es típico para partículas pequeñas que caen a través de un fluido viscoso; las partículas más grandes darían lugar a la creación de una estela turbulenta .
Para partículas con una pequeña velocidad de sedimentación, la difusión aumentará la complejidad de la trayectoria de la partícula hasta el fondo y el tiempo que tarda en sedimentarse en comparación con partículas con altas velocidades de sedimentación.

La velocidad de sedimentación (también llamada "velocidad de caída" o " velocidad terminal ") es una función del número de Reynolds de la partícula . Generalmente, para partículas pequeñas (aproximación laminar), se puede calcular con la Ley de Stokes . Para partículas más grandes (números de Reynolds de partículas turbulentas), la velocidad de caída se calcula con la ley de arrastre turbulento . Dietrich (1982) recopiló una gran cantidad de datos publicados a los que ajustó empíricamente curvas de velocidad de sedimentación. [21] Ferguson y Church (2006) combinaron analíticamente las expresiones para el flujo de Stokes y una ley de arrastre turbulento en una sola ecuación que funciona para todos los tamaños de sedimento, y la probaron con éxito contra los datos de Dietrich. [22] Su ecuación es

.

En esta ecuación, w s es la velocidad de sedimentación, g es la aceleración debida a la gravedad y D es el diámetro medio del sedimento. es la viscosidad cinemática del agua , que es aproximadamente 1,0 x 10 −6 m 2 /s para agua a 20 °C.

y son constantes relacionadas con la forma y suavidad de los granos.

La expresión para la velocidad de caída se puede simplificar de modo que se pueda resolver solo en términos de D . Usamos los diámetros de tamiz para granos naturales, , y los valores dados anteriormente para y . A partir de estos parámetros, la velocidad de caída viene dada por la expresión:

Alternativamente, la velocidad de sedimentación de una partícula de sedimento se puede derivar utilizando la Ley de Stokes suponiendo que el fluido está quieto (o quieto) en estado estable . La formulación resultante para la velocidad de sedimentación es,

Gráfico que muestra la relación entre el diámetro de las partículas de sedimento y la velocidad de sedimentación de Stokes.

donde  es la constante gravitacional ; es la densidad del sedimento;  es la densidad del agua ;  es el diámetro de la partícula del sedimento (comúnmente asumido como el diámetro medio de la partícula, a menudo mencionado  en estudios de campo); y  es la viscosidad molecular del agua. La velocidad de sedimentación de Stokes puede considerarse como la velocidad terminal resultante de equilibrar la fuerza de flotación de una partícula (proporcional al área de la sección transversal) con la fuerza gravitacional (proporcional a la masa). Las partículas pequeñas tendrán una velocidad de sedimentación más lenta que las partículas más pesadas, como se ve en la figura. Esto tiene implicaciones para muchos aspectos del transporte de sedimentos, por ejemplo, qué tan lejos río abajo puede ser transportada una partícula en un río.

Diagrama de Hjulström-Sundborg

La curva logarítmica de Hjulström

En 1935, Filip Hjulström creó la curva de Hjulström , un gráfico que muestra la relación entre el tamaño del sedimento y la velocidad necesaria para erosionarlo (levantarlo), transportarlo o depositarlo. [23] El gráfico es logarítmico .

Posteriormente, Åke Sundborg modificó la curva de Hjulström para mostrar curvas separadas para el umbral de movimiento correspondiente a varias profundidades de agua, como es necesario si se utiliza la velocidad del flujo en lugar de la tensión cortante límite (como en el diagrama de Shields) para la resistencia del flujo. [24]

Esta curva no tiene más que un valor histórico en la actualidad, aunque su simplicidad sigue siendo atractiva. Entre los inconvenientes de esta curva están que no tiene en cuenta la profundidad del agua y, lo que es más importante, que no muestra que la sedimentación es causada por la desaceleración de la velocidad del flujo y la erosión es causada por la aceleración del flujo . El diagrama adimensional de Shields es ahora aceptado unánimemente para el inicio del movimiento de sedimentos en los ríos.

Tarifa de transporte

Diagrama esquemático de los diferentes tipos de carga de sedimentos que se transportan en la corriente. La carga disuelta no es sedimento: está compuesta de iones disociados que se desplazan junto con la corriente. Sin embargo, puede constituir una proporción significativa (a menudo varios porcentajes, pero en ocasiones más de la mitad) de la cantidad total de material transportado por la corriente.

Existen fórmulas para calcular la tasa de transporte de sedimentos para los sedimentos que se desplazan en distintas partes del flujo. Estas fórmulas suelen dividirse en carga de lecho , carga suspendida y carga de lavado . A veces también pueden dividirse en carga de material de lecho y carga de lavado.

Carga de cama

La carga de fondo se mueve rodando, deslizándose y saltando (o saltando ) sobre el fondo, y se mueve a una pequeña fracción de la velocidad del flujo del fluido. En general, se cree que la carga de fondo constituye el 5-10% de la carga total de sedimentos en un arroyo, lo que la hace menos importante en términos de balance de masa. Sin embargo, la carga de material del fondo (la carga del fondo más la parte de la carga suspendida que comprende material derivado del fondo) a menudo está dominada por la carga de fondo, especialmente en ríos con lecho de grava. Esta carga de material del fondo es la única parte de la carga de sedimentos que interactúa activamente con el fondo. Como la carga de fondo es un componente importante de eso, juega un papel importante en el control de la morfología del canal.

Las tasas de transporte de carga de fondo se expresan generalmente en relación con el exceso de esfuerzo cortante adimensional elevado a cierta potencia. El exceso de esfuerzo cortante adimensional es una medida adimensional del esfuerzo cortante de fondo en torno al umbral de movimiento.

,

Las tasas de transporte de carga de fondo también pueden expresarse mediante una relación entre la tensión de corte del fondo y la tensión de corte crítica, que es equivalente tanto en el caso dimensional como en el no dimensional. Esta relación se denomina "etapa de transporte" y es importante porque muestra la tensión de corte del fondo como un múltiplo del valor del criterio para el inicio del movimiento.

Cuando se utiliza para fórmulas de transporte de sedimentos, esta relación normalmente se eleva a una potencia.

La mayoría de las relaciones publicadas para el transporte de sedimentos de fondo se expresan en peso de sedimento seco por unidad de ancho de canal (" ancho "):

.

Debido a la dificultad de estimar las tasas de transporte de carga de fondo, estas ecuaciones normalmente sólo son adecuadas para las situaciones para las que fueron diseñadas.

Fórmulas de transporte de carga de fondo notables

Meyer-Peter Müller y derivados

La fórmula de transporte de Meyer-Peter y Müller, desarrollada originalmente en 1948, [25] fue diseñada para grava fina bien clasificada en una etapa de transporte de aproximadamente 8. [19] La fórmula utiliza la no dimensionalización anterior para el esfuerzo cortante, [19]

,

y la no dimensionalización de Hans Einstein para la descarga volumétrica de sedimentos por unidad de ancho [19]

.

Su fórmula dice:

. [19]

Su valor determinado experimentalmente es 0,047, y es el tercer valor comúnmente utilizado para esto (además del 0,03 de Parker y el 0,06 de Shields).

Debido a su amplio uso, a lo largo de los años se han realizado algunas revisiones de la fórmula que muestran que el coeficiente de la izquierda ("8" arriba) es una función de la etapa de transporte: [19] [26] [27] [28]

[26]
[27] [28]

Las variaciones del coeficiente se generalizaron posteriormente como función de la tensión cortante adimensional: [19] [29]

[29]
Wilcock y Crowe

En 2003, Peter Wilcock y Joanna Crowe (ahora Joanna Curran) publicaron una fórmula de transporte de sedimentos que funciona con múltiples tamaños de grano en el rango de arena y grava. [30] Su fórmula funciona con distribuciones de tamaño de grano de superficie, a diferencia de los modelos más antiguos que utilizan distribuciones de tamaño de grano del subsuelo (y por lo tanto infieren implícitamente una clasificación de grano de superficie ).

Su expresión es más complicada que las reglas básicas de transporte de sedimentos (como las de Meyer-Peter y Müller) porque tiene en cuenta múltiples tamaños de grano: esto requiere considerar tensiones de corte de referencia para cada tamaño de grano, la fracción del suministro total de sedimentos que cae en cada clase de tamaño de grano y una "función de ocultación".

La "función de ocultamiento" tiene en cuenta el hecho de que, si bien los granos pequeños son inherentemente más móviles que los granos grandes, en un lecho de granos de distintos tamaños pueden quedar atrapados en bolsas profundas entre granos grandes. De la misma manera, un grano grande en un lecho de partículas pequeñas quedará atrapado en una bolsa mucho más pequeña que si estuviera en un lecho de granos del mismo tamaño. En los ríos con lecho de grava, esto puede causar una "movilidad igual", en la que los granos pequeños pueden moverse con la misma facilidad que los grandes. [31] A medida que se agrega arena al sistema, se aleja de la parte de "movilidad igual" de la función de ocultamiento hacia una en la que el tamaño del grano nuevamente importa. [30]

Su modelo se basa en la etapa de transporte, o la relación entre la tensión de corte del lecho y la tensión de corte crítica para el inicio del movimiento del grano. Debido a que su fórmula funciona con varios tamaños de grano simultáneamente, definen la tensión de corte crítica para cada clase de tamaño de grano, , como igual a una "tensión de corte de referencia", . [30]

Expresan sus ecuaciones en términos de un parámetro de transporte adimensional (donde " " indica adimensionalidad y " " indica que es una función del tamaño del grano):

es la tasa de transporte de carga volumétrica del lecho de la clase de tamaño por unidad de ancho del canal . es la proporción de la clase de tamaño que está presente en el lecho.

Se les ocurrieron dos ecuaciones, dependiendo de la etapa de transporte, para :

y para :

.

Esta ecuación alcanza asintóticamente un valor constante de a medida que se hace grande.

Wilcock y Kenworthy

En 2002, Peter Wilcock y TA Kenworthy, siguiendo a Peter Wilcock (1998), [32] publicaron una fórmula de transporte de carga de lecho de sedimentos que trabaja con sólo dos fracciones de sedimentos, es decir, fracciones de arena y grava. [33] Un modelo de transporte de carga de lecho de sedimentos de tamaño mixto que utiliza sólo dos fracciones ofrece ventajas prácticas en términos de modelado tanto computacional como conceptual al tener en cuenta los efectos no lineales de la presencia de arena en lechos de grava sobre la tasa de transporte de carga de lecho de ambas fracciones. De hecho, en la fórmula de carga de lecho de dos fracciones aparece un nuevo ingrediente con respecto a la de Meyer-Peter y Müller que es la proporción de fracción en la superficie del lecho donde el subíndice representa la fracción de arena (s) o grava (g). La proporción , en función del contenido de arena , representa físicamente la influencia relativa de los mecanismos que controlan el transporte de arena y grava, asociada con el cambio de un lecho de grava soportado por clastos a uno soportado por matriz. Además, dado que los intervalos están entre 0 y 1, los fenómenos que varían con incluyen los efectos de tamaño relativo que producen "ocultamiento" de granos finos y "exposición" de granos gruesos. El efecto de "ocultamiento" tiene en cuenta el hecho de que, si bien los granos pequeños son inherentemente más móviles que los granos grandes, en un lecho de tamaños de grano mixtos, pueden quedar atrapados en bolsas profundas entre granos grandes. Del mismo modo, un grano grande en un lecho de partículas pequeñas quedará atrapado en una bolsa mucho más pequeña que si estuviera en un lecho de granos del mismo tamaño, a lo que se refiere la fórmula de Meyer-Peter y Müller. En los ríos con lecho de grava, esto puede causar "movilidad igual", en la que los granos pequeños pueden moverse con la misma facilidad que los grandes. [31] A medida que se agrega arena al sistema, se aleja de la parte de "movilidad igual" de la función de ocultamiento a una en la que el tamaño del grano nuevamente importa. [33]

Su modelo se basa en la etapa de transporte, es decir , o la relación entre el esfuerzo cortante del lecho y el esfuerzo cortante crítico para el inicio del movimiento del grano. Debido a que su fórmula funciona solo con dos fracciones simultáneamente, definen el esfuerzo cortante crítico para cada una de las dos clases de tamaño de grano, , donde representa la fracción de arena (s) o grava (g). El esfuerzo cortante crítico que representa el movimiento incipiente para cada una de las dos fracciones es consistente con los valores establecidos en el límite de lechos de arena y grava puros y muestra un cambio brusco con el aumento del contenido de arena durante la transición de un lecho soportado por clastos a uno soportado por matriz. [33]

Expresan sus ecuaciones en términos de un parámetro de transporte adimensional (donde " " indica adimensionalidad y " " indica que es una función del tamaño del grano):

es la tasa de transporte de carga volumétrica del lecho de la clase de tamaño por unidad de ancho del canal . es la proporción de la clase de tamaño que está presente en el lecho.

Se les ocurrieron dos ecuaciones, dependiendo de la etapa de transporte, para :

y para :

.

Esta ecuación alcanza asintóticamente un valor constante de a medida que se hace grande y los símbolos tienen los siguientes valores:

Para aplicar la formulación anterior es necesario especificar los tamaños de grano característicos para la porción de arena y para la porción de grava de la capa superficial, las fracciones de arena y grava, respectivamente, en la capa superficial, la gravedad específica sumergida del sedimento R y la velocidad de corte asociada a la fricción superficial .

Kuhnley otros.

Para el caso en el que la fracción de arena es transportada por la corriente sobre y a través de un lecho de grava inmóvil, Kuhnle et al. (2013), [34] siguiendo el análisis teórico realizado por Pellachini (2011), [35] proporciona una nueva relación para el transporte de carga de lecho de la fracción de arena cuando las partículas de grava permanecen en reposo. Vale la pena mencionar que Kuhnle et al. (2013) [34] aplicaron la fórmula de Wilcock y Kenworthy (2002) [33] a sus datos experimentales y descubrieron que las tasas de carga de lecho predichas de la fracción de arena eran aproximadamente 10 veces mayores que las medidas y se acercaban a 1 a medida que la elevación de la arena se acercaba a la parte superior de la capa de grava. [34] También plantearon la hipótesis de que el desajuste entre las tasas de carga de lecho de arena predichas y medidas se debe al hecho de que la tensión de corte del lecho utilizada para la fórmula de Wilcock y Kenworthy (2002) [33] era mayor que la disponible para el transporte dentro del lecho de grava debido al efecto de protección de las partículas de grava. [34] Para superar este desajuste, siguiendo a Pellachini (2011), [35] asumieron que la variabilidad de la tensión de corte del lecho disponible para que la arena sea transportada por la corriente sería una función de la denominada "Función Geométrica de Rugosidad" (RGF), [36] que representa la distribución de las elevaciones del lecho de grava. Por lo tanto, la fórmula de carga del lecho de arena es la siguiente: [34]

dónde

El subíndice se refiere a la fracción de arena, s representa la relación donde es la densidad de la fracción de arena, es el RGF en función del nivel de arena dentro del lecho de grava, es el esfuerzo cortante del lecho disponible para el transporte de arena y es el esfuerzo cortante crítico para el movimiento incipiente de la fracción de arena, que se calculó gráficamente utilizando la relación de tipo Shields actualizada de Miller et al. (1977) . [37]

Carga suspendida

La carga suspendida se transporta en las partes inferiores y medias del flujo y se mueve a una gran fracción de la velocidad media del flujo en la corriente.

Una caracterización común de la concentración de sedimentos en suspensión en un flujo se da mediante el perfil de Rouse. Esta caracterización funciona para la situación en la que se puede cuantificar la concentración de sedimentos a una determinada altura por encima del lecho . Se da mediante la expresión:

Aquí, es la elevación sobre el lecho, es la concentración de sedimento suspendido a esa elevación, es la profundidad del flujo, es el número de Rouse y relaciona la viscosidad de remolino para el momento con la difusividad de remolino para el sedimento, que es aproximadamente igual a uno. [38]

Los trabajos experimentales han demostrado que los rangos van desde 0,93 a 1,10 para arenas y limos. [39]

El perfil de Rouse caracteriza las concentraciones de sedimentos porque el número de Rouse incluye tanto la mezcla turbulenta como la sedimentación bajo el peso de las partículas. La mezcla turbulenta da como resultado el movimiento neto de partículas desde regiones de altas concentraciones a concentraciones bajas. Debido a que las partículas se sedimentan hacia abajo, para todos los casos en los que las partículas no son neutralmente flotantes o lo suficientemente livianas como para que esta velocidad de sedimentación sea insignificante, existe un gradiente de concentración negativo neto a medida que se asciende en el flujo. Por lo tanto, el perfil de Rouse proporciona el perfil de concentración que proporciona un equilibrio entre la mezcla turbulenta (neta hacia arriba) de sedimentos y la velocidad de sedimentación hacia abajo de cada partícula.

Carga del material de la cama

La carga del material del lecho comprende la carga del lecho y la parte de la carga suspendida que proviene del lecho.

Tres relaciones comunes de transporte de material de lecho son las fórmulas de "Ackers-White", [40] "Engelund-Hansen" y "Yang". La primera es para arena a grava de tamaño granular , y la segunda y la tercera son para arena [41] aunque Yang luego amplió su fórmula para incluir grava fina. El hecho de que todas estas fórmulas cubran el rango de tamaño de arena y dos de ellas sean exclusivamente para arena se debe a que el sedimento en los ríos con lecho de arena se mueve comúnmente de manera simultánea como carga de lecho y carga suspendida.

Engelund-Hansen

La fórmula de carga de material del lecho de Engelund y Hansen es la única que no incluye ningún tipo de valor crítico para el inicio del transporte de sedimentos. Dice así:

donde es la no dimensionalización de Einstein para la descarga volumétrica de sedimentos por unidad de ancho, es un factor de fricción y es la tensión de Shields. La fórmula de Engelund-Hansen es una de las pocas fórmulas de transporte de sedimentos en la que no existe un umbral de "tensión cortante crítica".

Carga de lavado

La carga de lavado se transporta dentro de la columna de agua como parte del flujo y, por lo tanto, se mueve con la velocidad media de la corriente principal. Las concentraciones de carga de lavado son aproximadamente uniformes en la columna de agua. Esto se describe mediante el caso de endmember en el que el número de Rouse es igual a 0 (es decir, la velocidad de sedimentación es mucho menor que la velocidad de mezcla turbulenta), lo que conduce a una predicción de un perfil de concentración vertical perfectamente uniforme del material.

Carga total

Algunos autores han intentado formulaciones para la carga total de sedimentos transportados en el agua. [42] [43] Estas fórmulas están diseñadas principalmente para arena, ya que (dependiendo de las condiciones del flujo) la arena a menudo puede transportarse como carga de fondo y carga suspendida en la misma corriente o costa.

Mitigación de sedimentos por carga de fondo en estructuras de captación

Las estructuras de toma de agua de los ríos que se utilizan para el suministro de agua , la desviación de canales y el enfriamiento del agua pueden experimentar el arrastre de sedimentos de carga de fondo (del tamaño de la arena). Estos sedimentos arrastrados producen múltiples efectos nocivos, como la reducción o el bloqueo de la capacidad de toma, daños o vibraciones en el impulsor de la bomba de agua de alimentación y dan lugar a la deposición de sedimentos en tuberías y canales aguas abajo. Las estructuras que modifican las corrientes secundarias locales de campo cercano son útiles para mitigar estos efectos y limitar o prevenir la entrada de sedimentos de carga de fondo. [44]

Véase también

Referencias

  1. ^ Anderson, R (1990). "Ondulaciones eólicas como ejemplos de autoorganización en sistemas geomorfológicos". Earth-Science Reviews . 29 (1–4): 77. doi :10.1016/0012-8252(0)90029-U.
  2. ^ Kocurek, Gary; Ewing, Ryan C. (2005). "Autoorganización de los campos de dunas eólicos: implicaciones para la formación de patrones de campos de dunas simples y complejos". Geomorfología . 72 (1–4): 94. Bibcode :2005Geomo..72...94K. doi :10.1016/j.geomorph.2005.05.005.
  3. ^ Goudie, A; Middleton, NJ (2001). "Tormentas de polvo saharianas: naturaleza y consecuencias". Earth-Science Reviews . 56 (1–4): 179. Bibcode :2001ESRv...56..179G. doi :10.1016/S0012-8252(01)00067-8.
  4. ^ "Una tormenta de polvo se extiende desde el desierto de Gobi". Earthobservatory.nasa.gov. 13 de abril de 2006. Consultado el 8 de mayo de 2022 .
  5. ^ Neuendorf, Klaus KE; Mehl, James P. Jr; Jackson, Julia A., eds. (2011). Glosario de geología (quinta edición revisada). Alexandria, Virginia: American Geological Institute. pág. 800. ISBN 978-3-642-06621-4.OCLC 751527782  .
  6. ^ Wilson, WE y Moore, JE 2003. Glosario de hidrología, Instituto Geológico Americano, Springer, 248pp.
  7. ^ Charlton, Ro (2008). Fundamentos de la geomorfología fluvial . Londres: Rutledge. p. 234. ISBN. 978-0-415-33454-9.
  8. ^ Wohl, Ellen (2014). Ríos en el paisaje: ciencia y gestión . Wiley-Blackwell. pág. 330. ISBN 978-1-118-41489-7.
  9. ^ Ashton, Andrew; Murray, A. Brad; Arnault, Olivier (2001). "Formación de características costeras por inestabilidades a gran escala inducidas por olas de gran ángulo". Nature . 414 (6861): 296–300. Bibcode :2001Natur.414..296A. doi :10.1038/35104541. PMID  11713526. S2CID  205023325.
  10. ^ Roering, Joshua J.; Kirchner, James W.; Dietrich, William E. (1999). "Evidencia de transporte difusivo no lineal de sedimentos en laderas de colinas e implicaciones para la morfología del paisaje". Water Resources Research . 35 (3): 853. Bibcode :1999WRR....35..853R. doi : 10.1029/1998WR900090 .
  11. ^ Grant, J.; Walker, TR; Hill, PS; Lintern, DG (2013). "BEAST: un dispositivo portátil para la cuantificación de la erosión en núcleos de sedimentos intactos". Métodos en Oceanografía . 5 : 39–55. doi :10.1016/j.mio.2013.03.001.
  12. ^ Shields, A. (1936) Anwendung der Ähnlichkeitsmechanik und der Turbulenzforschung auf die Geschiebebewegung; En Mitteilungen der Preussischen Versuchsanstalt für Wasserbau und Schiffbau, Heft 26 (en línea archivado el 18 de julio de 2011 en Wayback Machine ; PDF; 3,8 MB)
  13. ^ Sharmeen, Saniya; Willgoose, Garry R. (2006). "La interacción entre el blindaje y la meteorización de partículas para la erosión de los paisajes". Procesos y formas del terreno en la superficie de la Tierra . 31 (10): 1195–1210. Bibcode :2006ESPL...31.1195S. doi : 10.1002/esp.1397 . S2CID  91175516.
  14. ^ Walker, TR; Grant, J. (2009). "Cuantificación de las tasas de erosión y estabilidad de los sedimentos del fondo en los sitios de acuicultura de mejillones en la Isla del Príncipe Eduardo, Canadá". Journal of Marine Systems . 75 (1–2): 46–55. Bibcode :2009JMS....75...46W. doi :10.1016/j.jmarsys.2008.07.009.
  15. ^ Dey S. (1999) Umbral de sedimento. Modelado matemático aplicado , Elsevier, vol. 23, núm. 5, 399-417.
  16. ^ de Hubert Chanson (2004). La hidráulica del flujo en canales abiertos: una introducción . Butterworth-Heinemann, 2.ª edición, Oxford, Reino Unido, 630 páginas. ISBN 978-0-7506-5978-9.
  17. ^ Whipple, Kelin (2004). "Rugosidad hidráulica" (PDF) . 12.163: Procesos superficiales y evolución del paisaje . MIT OCW . Consultado el 27 de marzo de 2009 .
  18. ^ ab Parker, G (1990). "Relación de transporte de carga superficial para ríos de grava". Revista de investigación hidráulica . 28 (4): 417–436. Código Bibliográfico :1990JHydR..28..417P. doi :10.1080/00221689009499058.
  19. ^ abcdefgh Whipple, Kelin (septiembre de 2004). "IV. Fundamentos del transporte de sedimentos" (PDF) . 12.163/12.463 Procesos superficiales y evolución del paisaje: notas del curso . MIT OpenCourseWare . Consultado el 11 de octubre de 2009 .
  20. ^ Moore, Andrew. "Conferencia 20: algunos cabos sueltos" (PDF) . Notas de la conferencia: Transporte de sedimentos fluviales . Universidad Estatal de Kent . Consultado el 23 de diciembre de 2009 .
  21. ^ Dietrich, WE (1982). "Velocidad de sedimentación de partículas naturales" (PDF) . Investigación de recursos hídricos . 18 (6): 1615–1626. Bibcode :1982WRR....18.1615D. doi :10.1029/WR018i006p01615. S2CID  128625889.
  22. ^ Ferguson, RI; Church, M. (2006). "Una ecuación universal simple para la velocidad de sedimentación del grano". Revista de investigación sedimentaria . 74 (6): 933–937. doi :10.1306/051204740933.
  23. ^ El perfil largo: procesos cambiantes: tipos de erosión, transporte y deposición, tipos de carga; la curva de Hjulström. coolgeography.co.uk. Último acceso: 26 de diciembre de 2011.
  24. ^ Temas especiales: Introducción a los movimientos de fluidos, transporte de sedimentos y estructuras sedimentarias generadas por corrientes; como se enseñó en: otoño de 2006. Instituto Tecnológico de Massachusetts . 2006. Último acceso: 26 de diciembre de 2011.
  25. ^ Meyer-Peter, E; Müller, R. (1948). Fórmulas para el transporte de cargas de fondo . Actas de la 2.ª reunión de la Asociación Internacional para la Investigación de Estructuras Hidráulicas. págs. 39–64.
  26. ^ ab Fernández-Luque, R; van Beek, R (1976). "Erosión y transporte de sedimentos de fondo". J. hidrol. Res . 14 (2): 127-144. Código Bib : 1976JHydR..14..127F. doi : 10.1080/00221687609499677.
  27. ^ ab Cheng, Nian-Sheng (2002). "Fórmula exponencial para el transporte de carga de fondo". Journal of Hydraulic Engineering . 128 (10): 942. doi :10.1061/(ASCE)0733-9429(2002)128:10(942). hdl : 10356/83917 .
  28. ^ ab Wilson, KC (1966). "Transporte de carga de fondo a alta tensión de corte". J. Hydraul. Div . 92 (6). ASCE: 49–59. doi :10.1061/JYCEAJ.0001562.
  29. ^ ab Wiberg, Patricia L. ; Dungan Smith, J. (1989). "Modelo para calcular el transporte de sedimentos por carga de fondo". Journal of Hydraulic Engineering . 115 : 101. doi :10.1061/(ASCE)0733-9429(1989)115:1(101).
  30. ^ abc Wilcock, Peter R.; Crowe, Joanna C. (2003). "Modelo de transporte basado en la superficie para sedimentos de tamaño mixto". Journal of Hydraulic Engineering . 129 (2): 120. doi :10.1061/(ASCE)0733-9429(2003)129:2(120). S2CID  40651180.
  31. ^ ab Parker, G.; Klingeman, PC; McLean, DG (1982). "Carga de fondo y distribución de tamaño en arroyos con lecho de grava pavimentado". Revista de la División de Hidráulica . 108 (4). ASCE : 544–571. doi :10.1061/JYCEAJ.0005854.
  32. ^ Wilcock, PR (1998). "Modelo de dos fracciones del movimiento inicial de sedimentos en ríos con lecho de grava". Science . 280 (5362): 410–412. Bibcode :1998Sci...280..410W. doi :10.1126/science.280.5362.410. PMID  9545213.
  33. ^ abcde Wilcock, Peter R.; Kenworthy, T. (2002). "Un modelo de dos fracciones para el transporte de mezclas de arena y grava". Water Res . 38 (10): 1194. Bibcode :2002WRR....38.1194W. doi : 10.1029/2001WR000684 .
  34. ^ ABCDE Kuhnle, RA; Reyezuelo, DG; Langendoen, EJ; Rigby, JR (2013). "Transporte de arena sobre un sustrato de grava inmóvil". Revista de Ingeniería Hidráulica . 139 (2): 167-176. doi :10.1061/(ASCE)HY.1943-7900.0000615.
  35. ^ ab Pellachini, Corrado (2011). Modelado del transporte de sedimentos finos sobre un lecho de grava inmóvil (tesis doctoral). Trento: Unitn-eprints.
  36. ^ Nikora, V; Goring, D; McEwan, I; Griffiths, G (2001). "Flujo en canal abierto promediado espacialmente sobre lecho rugoso". J. Hydraul. Eng . 127 (2): 123–133. doi :10.1061/(ASCE)0733-9429(2001)127:2(123).
  37. ^ Miller, MC; McCave, IN; Komar, PD (1977). "Umbral de movimiento de sedimentos bajo corrientes unidireccionales". Sedimentología . 24 (4): 507–527. Bibcode :1977Sedim..24..507M. doi :10.1111/j.1365-3091.1977.tb00136.x.
  38. ^ Harris, Courtney K. (18 de marzo de 2003). «Lecture 9: Suspended Sediment Transport II» (PDF) . Procesos de transporte de sedimentos en entornos costeros . Instituto de Ciencias Marinas de Virginia . Archivado desde el original (PDF) el 28 de mayo de 2010. Consultado el 23 de diciembre de 2009 .
  39. ^ Moore, Andrew. "Conferencia 21: Transporte de sedimentos suspendidos" (PDF) . Notas de la conferencia: Transporte de sedimentos fluviales . Universidad Estatal de Kent . Consultado el 25 de diciembre de 2009 .
  40. ^ Ackers, P.; White, WR (1973). "Transporte de sedimentos: nuevo enfoque y análisis". Revista de la División de Hidráulica . 99 (11). ASCE : 2041–2060. doi :10.1061/JYCEAJ.0003791.
  41. ^ Ariffin, J.; AA Ghani; NA Zakaira; AH Yahya (14–16 de octubre de 2002). "Evaluación de ecuaciones sobre la carga total de material del lecho" (PDF) . Conferencia internacional sobre hidrología urbana para el siglo XXI . Kuala Lumpur .
  42. ^ Yang, C (1979). "Ecuaciones de potencia de corriente unitaria para carga total". Revista de hidrología . 40 (1–2): 123. Código Bibliográfico :1979JHyd...40..123Y. doi :10.1016/0022-1694(79)90092-1.
  43. ^ Bailard, James A. (1981). "Un modelo energético de transporte de sedimentos de carga total para una playa plana en pendiente". Revista de investigación geofísica . 86 (C11): 10938. Código Bibliográfico :1981JGR....8610938B. doi :10.1029/JC086iC11p10938.
  44. ^ Natato, T.; Ogden, FL (1998). "Control de sedimentos en tomas de agua a lo largo de ríos con lechos arenosos". Journal of Hydraulic Engineering . 126 (6): 589–596. doi :10.1061/(ASCE)0733-9429(1998)124:6(589).

Enlaces externos