MODELOS MATEMATICOS

EN
HIDROLOGIA

JULIO SANCHEZ ORDOEZ

BOCOTA , COLOMBiA
DICIEMBRE. 1974

l
MODELOS MATEMATICOS EN HIDROLOGI.A

JULIO SANCHEZ ORDOEZ

Ingen iero Civil

BOGOTA, COLOMBIA
Diciembre 1974
 I N D ICE

MODELOS MATEMATlCOS EN HIDROLOGIA

pgina
INTRODUCCION .............. .... . 1

CAPITULO I - CONCEPTOS GENERALES

1.1 SISTEMAS HIDROLOGICOS .. . ........ 3

1. 2 MODELOS DE SIMULACION .... , .... 5

1.2.1 Simulacin - Generalidades................................ 5
1.2.2 Modelos Hidrolgicos de simulacin ......................... 8
1.2.3 Elementos de un modelo de simulacin............ ........... 10
1.2.4 Aplicacin de un modelo de simulacin...................... 12

CAPITULO II - CLASIFICACION DE MODELOS HIDROLOGICOS

2.1 DETERMINISTlCOS y ESTOCASTICOS .... ....... ..... 15

2.2 ANALlTlCOS y EMPIRICOS ........... ..... ................ 18
2.3 MODELOS A EMBALSES Y CON OPERADOR PLWIOHIDROLO
GICO ............................... ~ ~................ o" 19
2.4 LINEALES Y NO LINEALES ...................... .......... 19
2.5 DESCRIPTIVOS Y CONCEPTUALES ................... ....... 19

CAPITULO 111 - MODELOS DETERMINISTlCOS

3.1 ANTECEDENTES .......... " .... ... 23

3.2 SIMULACJON DE ESCURRIMIENTO EN ZONAS INUNDABLES
(Modelo SOGREAH - UN ESC O .. . .. .. .. .. .. ... . .. .. . .. .. 25

3.2.1 Caracterfsticas del problema.. .. ........ .......... .. .... .... 25

A. Generalidades......................................... 25
 pgina

B. Conguracin de lo cuenca .... ".. ""."" ...... """"""" .. "" ........... .. 26

3.2.2 Desc:ripcin del modelo ....... .. " 26 ~ " .......................... " ..

A. Leyes entre embalses fluviales " " . " " .. " .... " " " .. " " " .. " " 26
B. Leyes entre embalses no fluviales ...... ".... "" .. """ ...... ".... ".. ".. " 29
3.2.3 Formulacin matemtica del modelo ,. " " "" 31
.. " .... "" ...... " " .. " " " " " . " ....
A. Hiptesis previ 0$ .. "" .... " .. " .... " .... " .. "" .. " " " .. " ........ " .. " " ............ " 31
B. Ecuacin de continuidad . " . ".. ".. ".. """"" .... """" .. """"" .. """ . " 31
C. Ecuaciortes de descargas ".... "........ ".... "" ........ """" """" """"".... 32
Unin de tipo fluvial"" .... "" ............ " s." ...... " "" .... " .. """" .... " 32
lJrlin de tipo vertedero .... " .. "" .. "" .. ".. " .. "" .. " " " .. " "",, " "...... "" 35

3.2.4 Desa.rrollo del m:odelo .... "....... ". "" ........ "................ " 36
A. D,strl'buc .Ion
' topo "oglca
.. "" ...... " .......... " .. " .... " ...... " .... " " " .... " 36
B. Conformacin de las motrices . " .. " .......... " " ........... " .............. " .. 38
C. Condiciones Irmites . ". ".... o. . . . . o.o.o.o. . . . . . . . o. .. o. ..... o. . . . . . . . . . . . . . 42
Condiciones aguas arriba ....................................... 42
Condiciones aguas abajo ............ # ................... . 44
D. Tcnica de solucin ............ o. . . . . . . . . . . . . . . . . . . 'o. . . . . . 45

3.2.5 Esquema de I programa de computador ......................... . 47

3.3 MODELO DE LA UNIVERSIDAD DE TEXAS (Modelo Stonford IV

Mocfificodo) ................................................. . 48
3.3.} Antecedentes ..................... . ................................... 48
3.3.2 Caracterrsticas de I modelo ................................ . 49
A. Origen ............ o. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . o. . . . . . . . . . . . . o. . . . . . . 49
B. Almacenamiento de intercepcin ................................ . 52
C. Almacenamiento superficial ..... 'o. . . . . . . . . . . . . . . . . . . 'o. .. . 52
D. Infiltracin ... ... , ........................................... . ....... 55
E. Fluio o travs de la zona no saturada ......................... .. 59
F. Flujo sllbsuPf!mcial ......... o. . . . . . . . . . . . . . . .- . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 60
G. F luio subterrneo ......................................................... .. 62
H. 1:.._ .. "
I;;vuporaclon ..
y evapotransprraclon ..'
..................................... . 64-
I . Propaga.cir, de descargas .......................................... .. 66

11
 pgina

3.3.3 Aplicaciones del modelo ............ ""............ 66

3.4 MODELO SSARR (Streamflow Synthesis ond Reservoir Regulotion).. 70

3.4. I General idodes ... . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70

3.4.2 Evolucin histrica del modelo . .. . ...... ..... 70
3.4.3 Caractersticas del programa SSARR.......................... 72
A. Conceptos hidrolgicos incluidos ......... 72
B. Subdivisin de la cuenca............ . . . . . . . . . . . 72
C. Intervalos de Nempo . . . . . . . . . . . . . . 72
D. Datos de en trada ...... . . . . . . . . . . . . . . . . .. . . . 73
E. Usos del modelo................................. .. 73

3.4.4 Tipos bsicos de simulacin empleados........................ 73

3.4.5 Transformacin IIuvia-descarga.............................. 74
A. Determinacin de lo precipitacin media en la cuenca . 76
B. Distribucin de la lluvia ........................................ ....... .. 77
C. C!eulo de la escorrenta total ......... 77
D. Clculo de la humedad del suelo ....................... .. 78
E. Determinacin del flujo bose ........................... .. 78
F. Separacin de los flujos superficial y suhsuperficial .... 80
G. Propagacin de las componentes ............. 81
H. Definicin de las caractersticas de la cuenca ........... .. 81

3.4.6 Propagacin en el cauce ................................... . 83

A. Ecuaciones bsicos . ...................................... lO ' . . . .. 83
.' geornetnca
.. . ................................... .
B .. Interpretaclon 87
C. Esquema de clculo .................................. .. 87
D. Propagacin con efecto de remanso ................ 87

3.4.7 Aplicaciones prcticas del modelo ................. 91

CAPITULO IV - MODELOS ESTOCASTlCOS

4.1 ALGUNOS ASPECTOS DE LOS EVENTOS HIDROLOGICOS .. 93

4.1.1 Procesos hidrolgicos - Defin iciones . . . . . . . . .. . . . ... 93
4.1.2 Caractersticas de los datos hidrolgicos ...... '" . ... 93

111
 pgina
4.1.3 Naturaleza de las magnitudes hidrolgicas.................... 96

4.2 APLICACIONES DE ESTADISTICA EN HIDROLOGIA , . 96

4 . 2. 1 Nociones generales ....................................................... ~ .................. .. 97
4 .... I rr uc rOI'l8S emp.' rlcas .................. ..............................
2 2 D'st'b lO . . . . . . . . . . . . . . . . .. 98
4.2.3 Parmetros estadrsticos ....................................................................... . 99
A. Medidas de tendencia central ....... 99
8:.. ~didas de d;spersin . o......................................................... 100
C. Medidas de asimetrfa . . . . . . . . . . . .... .. 100

- y regreslon
2 4 Corre loClon
4 .... - ..................................................................... 101
4.2.5 Distribuciones de probabilidad .................. o..... 104
A. Distribucin normal ......... ............................... o . . . . . . . . . . . . . . . . . . . . . . .. 104
8. Distribucin lognonnal .................................................' .............. .. 105
C .. D.strl .- gcamma .................................................................. ..
'bUClorl 106
D. Esquema para seleccin de la distribucin ..... 106

4.3 SELECCION DE UN MODELO PARA SINTETIZAR CAUDALES... 108

4.3.1 El modelo de Maricov......................................... 109
A. Caudales con distribucin normal......................... 110
B. Caudales con distribucin lognormal ... " " 112
C. Caudales con distribucin gamma ....................... 112

4.3.2 Modelas estacionales de Markov .,... .. . .. 112

A. Caudales con distribucin normal y lognormal ...... 113
B. Caudales can distribucin gamma .............. 113

REFERENClAS . ....... ...... . ... ... ........ .... ................... 115

IV
 INTRODUCCION

En los lHmos aos se ha incrementado considerablemente el uso de computadores digi-

tales en el anlisis de datos hidrolgicos paro diversidad de fines. En un principio los
computadores se usaran bsicamente paro labores rutinarios toles como conversin de do
tos, propagacin de hidrogramas, clculo de remansos etc.; los avances logrados en la
tecnolaga de computadores han tomado prcticos lo programacin lineal y dinmico, por
ejemplo, y avanzados tcnicas estadsticos de anlisis; el anlisis numrico y los mode-
los matemticos han aumentado su grado de sofisticacin como resultado de la disponibi
lidad de computadores cado vez ms grandes y rpidos. -

Hasta antes del advenimiento de estos tcnicos, los parmetras derivados analticamente
pora lo mayora de los modelos matemticos provenan de regresiones y correlaciones;
estos pormetros estaban, par lo tonto, sujetos a condiciones de independencia, normali
dad y linealidad, condiciones a los cuales no pueden restringirse los parmetros de lo
mayora de los sistemas hidrolgicos, o tal punto que lospormetros tipo coeficiente de
madelos tan restringidos pueden llegar o no tener ningn significado fsico.

La simulacin puede llevarse a cabo con mecan ismos analgicos o con computadores di
gitales; en el primer coso lo simulacin es posible usando la similitud que existe entre-
los expresiones matemticos que describen el flujo elctrico yel flujo del aguo; en for-
ma digital, el computador se usa paro resolver numricamente toles ecuaciones. Cada
uno de estos tipos tiene ventajas y desventajas, osi como campos de aplicacin espec-
ficos, pero la computacin digital es extremadamente til en hidrologa de superficie,
por ejemplo poro resolver situaciones probabilsticos en los cuales las variables (preci-
pitaciones y escurrimientos, por ejemplo), varan aleatoriamente dentro de lmites proba
bies. -

El presente trabajo hoce referencia o 'o simulacin con computadores digitales, e nclu
ye lo descripcin de algunos modelos de amplio uso, exponiendo sus fundamentos y con-:
ceptos bsicos.

I
 CAPITULO

CONCEPTOS GENERALES

1.1 SISTEMAS HIDROLOGICOS

A travs de la historia de la humanidad, el hambre siempre na estado lleno de

curiosidad y se ha interrogado acerco del media ambiente donde vive. Con es-
te deseo y curiosidad de conocer la naturaleza, lo nocin de un cc lo hidrol-
gicose desarroll desde pocas muy tempranos. Hacia fines del siglo I A.C.,
Marcus Vitruvius Pollio observ y describi lo siguiente: (ref 14) "El sol calien-
ta el agua de los arroyos, ras, lagos y mares, formando vapores que se elevan
paro formar nubes. Las nubes son transportadas por el vi en to y, por lo perturba
cin que sufren cuando ascienden montaas y forman tempestades, se rompen y
disperson sobre la tierra. Las vapores, nubes y exhalaciones que se elevan de
la tierra parecen flotar gracias al intenso calor, los fuertes vientos, el aire h-
medo y su gran proporcin de agua. As, cuando de la frialdad de la noche asis
tida por lo oscuridad, los vientos se levantan y se forman las nubes en los lugO
res hmedos, el sol en su ascenso golpea la tierra con gran poder, caliento er
aire, y eleva las vapores y el roco al mismo tiempo". Esta descripcin es po-
tica y algo mstica, pero sin embargo refleja el primer concepto cientfico racia
nal sobre el ciclo hidrolgico. -

Los hidrlogos madernos consideran el ciclo hidrolgico como un gran sistema

hidrolgico. Con el objeto de alcanzar propositos prcticos, la hidrologa cien
tmca debe ir ms all de hacer una descripcin terica pura de naturaleza cuo
litativa de sus problemas, puesto que el hidrloga prctico reclama cantidades-
y no est interesado en cmo funciona el ciclo hidrolgico ni en las expresio-
nes matemticos que rigen sus partes, sino en conocer las cantidades de aguo
involucradas en l. Con el concepto de sistema, es posible dar un enfoque ms
cuantitativo y racional o lo hidrologa.

El concepto de sistema no es, por lo tanto, nuevo pora nosotros; sistemas de to-
das clases nos son comunes. Un sistema conlleva un grupo de elementos interde
pendientes o interactuantes que conforman un tado un ificado. Por ejemplo, er
sistemo solar, compuesto de una variedad de cuerpos celestes que operan como
una unidad; el ecosistema, que relaciona los seres vivientes y el ambiente fsico
donde habitan; los sistemas del cuerpo humano (circulatorio, digestivo, etc. l,son

3
 todos excelentes ejemplos de componentes inter-actuantes relacionados pera cons
tltuir un sistema unificado.

Una caracterstica importante de cualquier sisteme, es que sus componentes estn

relae:ionados entre s de tal manera que un cambio en uno de ellos afectar a
los dems, e implicar por lo tanto un reajuste total hacia un nuevo equilibrio.
El anlisis de un sistema implico, entonces, lo comprensin de las caractersti-
cos de coda elemento y el conocimiento de cmo se relaciona con el sistema to
tal; tal conocimiento nos llevar a predecir los principoles efectos resultantes -
de cualquier cambio porticular en alguna componente especfica del sistema.

El anlisis de sistemes podra ser considerado como "anl isis de un cuerpo libre",
puesto que siempre es necesario establecer las fronteros del sistema, cuya aeer-
toda delimitacin es una clave importante pora el xito del anlisis. Esto puede
vel$e muy claramente en un ejemplo que se relaciono con la estabilidad de ta-
ludes ocurrido recientemente en lo costa de California. En estos reas muchos
de los lboles de races profundos fueron remplazados par pastos; e I pasto prote
ge o lo superficie del suelo de lo erosin, pero las races profundas de los rbo-:
les tenan dos influencias adicionales sobre el sistema: 1) extraan aguo desde -
profundidades considerables y contribuan as a mantener un e levado grodo de hu
medad en los taludes, y 2) las races eran un elemento estabilizante del suelo -
profundo. As, el retiro de los rboles produjo un reajuste en el sistema, que se
tradujo en deslizamientos hasta alcanzar pendientes menores; es muy posible que
un hidrlogo experimentodo, 01 mirar el sistema como un toda, hubiese podido
prever estas reajustes. Aquellos que retiraron los rboles consideraron solo una
porte de I si stema .

Un sistema se puede definir como un conjunto de elementos o procesos fsicos,

qumicos o biolgicos unidos o travs de alguna forma de interdependencia, que
acto sobre un grupo de variables de entrado para convertirlos en las de sali-
do.

Los sistemas hidrolgicos naturales son tan complejos que an no se han deserro
lIado leyes exactas que puedan explicar completo y exactamente los fenmenos-
hidro/gc:os naturales.

Codo uno de los elementos o procesos integrantes de un sistema hidrolgico es

el resultado de complicadas interrelaciones de muchos factores de gran variabili
dad espacial y temparal, cuyos caractersticos fsicas difcilmente pueden medir
se o computarse directamente. -

Un sistema hidrolgico es fsico, puesto que pertenece 01 mundo real, secuencial,

porque consto de una entrado, solida y Un medio de trabajo (materia, energa o
informacin) conocido como "throughput" que pasa o travs del sistema; es ade-
ms, dinmico porque recibe ciertos entrados cuantitativas y acta de manera -
controlado bajo restricciones definidas pora producir uno serie de solidos cuanti-
..
 tativas.

la relacin explcita entrada -salida de un sistema se puede reDresentar matem

ticamente por una ecuacin del tipo:

y (t} -== 3 x (ti

donde Y (t) es la salida

X (t) es la entrado
% es el operador de tronsformocin o funcin de trarsferenca, que re
presenta la accin intemo del sistema sobre la entrado X t) para cO'1vertirla -
en la salida Y (t). Un hidrograma unitario es un .. emplo de op.. rodor de trons
formacin paro pasar de lluvias a descargas.

1.2 MODELOS DE SIMULACION

1.2.1 Simulacin - generalidades

la simulacin matemtica es una tcnica para producir una serie de eventos a

travs de expresiones matemticas que representan las situaciOnes del mundo fr-
sico bajo estudio; los simulaciones generalmente son efectuados por computado-
res que aplican iterativamente esos expresiones, o partir de ciertas condiciones
iniciales, y progresan en el tiempo, repitiendo las relaciones funcionales entre
los diferentes elementos del sistema.

La simulacin no es una tcnico nueva; por este trmino se entiende el antiguo

arte de la construccin de modelos destinados tanto a la apreciacin esttica -
(pintura, escultura) como a los estudios cientficos: modelos reducidos o analfti
COS desde el punto de vista cientfico, el objetivo fundamental de la simulacin
es el conocimiento del futuro.

Los mtodos deductivos, por lo tanto, son clasificados dentro de la filosofa es-
peculativa y han sido estudiados desde los filsofos griegos has!" el siglo XVII
con Sir Francis Bacon.

Aparecieron entonces sistemos filosficos que intentaron construir una fsica en

lo cual la explicacin cientfica era remplazodo por simples "nalogas con expe
riencias de la vida diaria. En 1.620, Bacon reconoci las limitaciones de 10-
filosofa especulativo como metodologa de la prediccin del futuro: "Lo razn
sola no tiene ninguna capacidad de previsin; la consigue solamente en combi-
nacin COn lo observacin. Los mtodos de prediccin de la razn estn conte-
nidos en las operaciones lgicas par medio de los cuales construimos un orden
dentro de lo moterio observada y sacamos conclusiones, a trovs de inducciones
lgicos". IReichenbach, 1. 95])

5
 Es as( como la filosaffa cienttr.ca o mtodo cientmco, puede desarrollarse en
cuatro fases :

Primera fase: Identificacin y observacin del sistema f(sico de inters.

Segunda fase: Formulacin de una hiptesis (o de un modelo) que intente expl.!.

car los observaciones del sistema.

Ten:era fase: Prediccin de I comportamiento del sistema en base a esta hipte-

sis, mediante deducciones matemticas o lgicas (obtener soluc~
nes pora el modelo).

Cuarta fase: Realizacin de experiencias pora comprobar la validez de la hi-

ptesis (o modelo).

Sin embargo, la experiencia demuestra que no siempre es posible la real izadn

de estas cuatro fases en todos los sistemas, y es allr entonces donde la simula-
cin puede ser un substituto til para la fase o foses que estn causando dificul
tades. Como ejemplo de dificultades que pueden tomar impracticable en hidro=
logfa alguna de las fases del mtodo cientffico, pueden citarse :

En la primera fase: En general existen imposibilidades f(sicas o econmicas pora

Observar en el mundo real los procesos constitutivos de un sistema hidrolgico o
de reCUr505 hfdricos.

.En la segunda fase: Estos sistemas son muchas 'teces tan complejos que es im-
. '
posible, al menos por ahora, describirlo exactamente con ecuaciones matematr- .
caso

En la ten:era Fase: En el caso de ser posible la representacin matemtica de

un proceso, podrTa ser impracticable, sin emborgo, obtener una solucin poro las
ecuaciones a travs de tcnicas directas.

En la cuarta fase: Puede ser diffcil o sumamente costoso realizar experiencias

vlidas en los modelos que describen el sistema.

las observaciones de la entrada y la salida de I sistema deben usarse coma in-

formacin adicional para determinar la estructura del mismo y los pormetros des
conocidos; tambin es esencial para ello un conocimiento de la formo generar
de las relaciones entre las cantidades individuales y sus valores ff,cos permisi-
bles. ~(t la constl'\lCcin de un modelo significa lo determinacin de un siste-
ma ptimo de ecuaciones que correspondan de lo mejor manera posible a las
funciones de entrado y salida observadas y o un conjunto dada de limitaciones.
las leyes fundamentales de lo f(sica (ecuaciones de energfa, balances hfdricos,
ecuaciones del movimiento, etc. \ y los resultados experimentales, sirven como
base poro la identificacin del sistema hidrolgico.
6
la simulacin numrica de los procesos hidrolgicos que ocurren en uno cuenco
ha surgido como un mtodo potec,cialm"nte til paro lo previsin de eventos hi-
drolgicos, y como una herramienta poro mejorar notablemente la comprensin de
tales procesos. Desde el punto de visto de aplicacin a lo ingeniera, la hidro
loga generolmente proporciona informacion"s como magnitud y frecuencia de-
avenidos, caudales de diseo de obras hidrulicos, Mc., mientros que desde el
punto de visto de lo hidrologa cientfico, el hombre est interesado en i"cre-
mentar sus conocimientos sobre la estructuro intema de coda uno de los proc~sos
fsicos que ocurren en los distintos fos .. s del ciclo hidrolgico.

los tcnicos tradicionales de anlisis hidrolgico han sido en g"'ne ro I fragmenta-

rios, con mtodos diferentes desarrollados poro cada problema especfico. Los-
mtodos referentes a los problemas de avenidas, por ejemplo, se basan en el -
anlisis de series histricas, variaciones de lo frmula "racional", variaciones del
mtodo del hidrograma unitario, anlisis de correlaciones y varias frmulas espe-
cficos. En "stas mtodos donde las caractersticas de las descargas se estiman o
portir de lluvias, una de las mayores inexactitudes o incertezos es la determina
cin de la lluvia neta (o "fectivo) y la manero como es afectada por los condr-
ciones anteriores de la cuenca.

La simulacin se ha definido como "el desarrollo y aplicacin de modelos mate-

mticos para representar la vcriacin temporol de la interaccin de los procesos
fsicos". Lo simulacin numrico ofrece un importante avance en la solucin de
problemas hidrolgicos segn se ve en las siguientes proposiciones:

o- En muchos lugares existen mejores y ms confiables registros de lluvias que

de descargas.

b - En muchos lugares urbanos o rurales han ocurrido cambios en los cuencos en

el tronscurso del tiempo; la simulacin numrica basada en caractersticas fr
sicas reales, ofrece la mejor posibilidad de evaluar los cambios en las desear
gos originados por modificaciones de lo cuenca conocidas o anticipadas, na-
turoles o artificiales.

La simulacin con computadores tiene las siguientes importantes ventajas:

a - El sistema puede probarse sin destruirlo.

b - Pueden probarse modificaciones propuestas de los sistemas existentes.

c - Pueden probarse muchos posibilidades en cortos perradas de tiempo.

d - Pueden ensayarse diseos hipot';ticos para estudios de factibilidad o compa-

racin con sistemas alternativos.
e - Se incrementa el conocimiento interno del sistema estudiado.

7
 1.2.2 Modelos hidrolgicos de simulacin

En el pasado se han usado ampliamente modelos hidrulicos para simular procesos

hidrolgicos, pero esto tal vez seo debido ms o tradicin que o uno justificacin
tcnico o econmico, puesto que muchos de estos procesos no pueden representar
se con modelos hidrulicos sin grandes distorsiones; adems tanto paro la reproduc
cin del proceso con la ayuda del modelo hidrulico como para lo simulacin con-
computadores electrnicos, es necesario un modelo matemtico del proceso
lo importancia de los modelos de simulacin para la investigacin cientrfiea fue
establecido por Rosenblieth y Wiener (1.945): "Ninguno parte substancial del uni
verso es tan simple que puedo ser comprendida y controlada sin abstraccin. Lo
abstraccin es lo substitucin de la porte del universo en estudio por un modelo
semejante, aunque de estructura ms simple. Las modelos constituyen por lo ton-
to, uno necesidad primordial de cualquier procedimiento cient(fico.

Un modelo cient(fico puede definirse como uno abstraccin de un sistema real

que se presto o propsitos de prediccin y cantrol. A travs de l, un analis-
ta tendrfo condiciones para determinar en qu proporcin las modificaciones de
determinados aspectos del sistema afectaran otros aspectos o el sistema en conjun
too En hidrologa, el uso de modelos 5010 ha sido posible gracias o la octual -
disponibilidad de computadores analgicos y digitales, la cual ha permitido desa
rrollar, por ejemplo, complejos modelos de lo fose precipitacin -escorrenta de!
ciclo hidrolgico.

Estos modelos basados en las leyes fsicas que controlan los procesos constitutivas
del sistema, son parcialmente subjetivos, puesto que prcticamente en ningn ca-
so existe un mtodo generalmente aceptado pora describir tales leyes fsicas. Asr,
las aproximaciones empricas util izadas afectan no solamente lo exactl tud de la
prediccin logrado con el mode lo, sino lo respuesta aparente de otros componen
tes del sistema hidrolgico. Adems, el criterio de exactitud o bondad del ojwii
incide directamente en el nmero y clase de los parmetros "ptimos que descri
ben las caractersticas fsicas de una cuenca dado. Por lo tanto los cualidades-
deseables en un mode lo contrastan entre s(:

a - El modelo deber ser uno aproximacin razonablemente precisa del sistema

real y contener el mayor nmero posible de los aspectos importantes del mis-
mo.

b - El modelo no debero ser tan complejo que se tome difcil de comprender o

manipular.

El uso potencial ms importante de un modelo es lo prediccin paro propsitos

de diseo. La bondad de un mode lo no est re (acionada solamente al criterio
mismo del ajuste obtenido, sino o la sensibilidad de la prediccin, cama canse-
8
cuencia de modificaciones en los parmetros usados. Estos dos medios de juzgar
un modelo deben ser tenidos en cuenta tanto por el que construye un modelo co-
mo par el que lo uso. f"ingn modelo es perfecto y por eso, codo uno est en
un estado contrnuo de evolucin, y cada cambio meiora el esquema concebido
originalmente.

Por lo tonto, deben tenerse muy ~n cuenta las limitaciones de la simulacin nu

mrica, tales como lo disponibilidad de datos de entrada adecuados; el procese
de simulacin en si, depende de las relaciones funcionales escogidas para repre
sentar los procesos frsicos que ocurren en la cuenco y de los parmetros usodos-
para describir las caractersticas de lo misma. El arte de desarrollar un modelo
de simulacin satisfactorio estribo en alcanzar un grado apropiado de fidelidad
en la reproduccin de tos distintos procesos individuales y su interrelacin, omi-
tiendo los detalles insignificantes e innecesarios que slo agregan complicacin
al modelo.

Los procesos hidrolgicos, por tanto, deben considerarse como sistemas dinmicos
con estructuras complejas y parcialmente desconocidas; as, el problema de simu
larlos se puede esquematizar en tres etapas:

a - Delimitacin del sistema dE' inters y establecimiento de los caroctersticas

que lo identifican;

b - Expresin del fenmeno con relaciones matemticos que reproduzcan los pro-
c"sos fsicos en formo aproximado; puede ocurrir qu~ el fenmeno no tenga
ecuaciones que lo representen o que tenga E'cuaciones diferenciales capaces
de simularlo con precisin; salvo este ltimo caso, la descripcin del fen-
meno es forzosamente incompleta;

c - Determinacin de los parmetros del modelo, los cuales vendrn a suplir esa
deficiencia en la expresin matemtica del proceso.

En resumen, el d,.sarrollo de un modelo matemtico operativo requiere dos pasos.

El primero es la creacin dE' un modelo conceptual que represente los elementos
y sistemas del mundo real;E'sta conceptualizacin se baso en lo informacin cono
cida y en hiptesis relatvas o los componQntes del sistema y sus interrelacioneS;
en general SE' formulo 'm trminos de los datos disponibles.

El segundo paso es lo transicin d 4 1 modelo conceptual al modelo de computador

propiamente dicho. Durante esto etapa se trotan de pxpresar en forma motemti
ca los procesos y relacion .. , identificados por el modelo conceptual. Es decir, -
se trata de la convE'rsin de los conceptos relativos al mundo real, en trminos
que puedan programarse en un computador. Este paso generalmente implico ms
simplificaciones y por consiguiente prdida de informacin que viene a sumarse
o la yo ocurrido en el paso del mundo real al modelo conceptual.

9
 La adquisici~ de datos adicionales usualmente signiflca uno mejora en el made
lo conceptual, lo cual provee uno base paltl el perfeccionamiento del modelo de
computador. La salida del modelo puede, por supuesto, compaltl!Se con las sal
das observadas en el mundo real y, si existen discrepancias entre ellas, se req4t
rirn ajustes tanto en el modelo conceptual como en el de computador. -

1.2.3 Elementos de un modelo de simulacin

En el campo de la hidrologa, los elementos de un modelo matemtico pueden

clasificarse en " cate garras (ref. 3) :

a - Componentes
b - Variables
c - Relaciones funcionales
d - Parmetros

0- Componentes: Son los elementos escogidos para integrar el modelo en

funcln d los Objetivos del mismo y de la importancia que se presume que estos
elementos tengan sobre el sistema en conjunto; entre las componentes ms comu-
nes de modelos de simulacin hidrolgica encontramos:

- Lluvia
- Escurrimiento supe rficial
- Evapotranspiracin
- Demando de agua para consumo urbano
- Demando de agua para irrigacin
- Reservas de agua en embalses
Demanda de agua para otros fines (generacin de energa, etc)

b - Va r i a b I e s : Son una medida o representacin cuantitativa de las compo-

nentes del sistema hdrica, que presentan variacin espac ial o temporal; entre ellas
podemos citar:

- Precipitacin horaria dentro de un rea especfica

- Escurrimiento medio mensual en una seccin
- Descargo instantnea en uno seccin
- Consumo horario de agua potable
- Demando mensual de agua paro irrigacin
- Demanda instantnea de una planta hidroelctrica

Pueden clasificarse estas variables dentro de tres grupos :

10
Va r i o b I e s e x gen o S : son los variables independien.tes o de entrada al mo
delo y se suponen conocidas previamente; son independientes del sistemareprese
todo por el modelo, es decir actuantes en el sistema pero no influidas por l. Es-
tas variables se pueden representar de dos moneros: como pormetros, dados por-
las condiciones del problema y previamente determinados y lerdos por el compu-
tador cama dotas de entrada, o como variables estadsticos que pueden ser gene
radas internamente par el computador en el desarrollo del modelo.

Voriab les endgenos o de sal ida: son los variables dependientes del
sistema, generadas por lo interaccion entre las variables exgenos y las variables
de estado, de acuerdo con los caractersticos operacionales del sistema.

Variables de estado: son las que describen el estado de un sistema o de

uno de sus componentes a travs del tiempo (puede ser al principia, 01 fin de un
intervalo de tiempo dado, o an durante el transcurso de un corto perodo). Es-
tas variables interocton con los dos anteriores siguiendo los relaciones funciono
les establecidas previamente. Por lo tanto, dependen no solo de los volares di
una OMS variables exgenas relativas o un corto perodo de tiempo anterior,
sino tambin de ciertas variables de salida del sistema en intervelos de tiempo
anteriores.

En la simulacin de lo propagacin de descargas en un trecho de r(a, por ejem-

plo, lo descargo de entrado 01 trecho sera una variable exgena, el almacena-
miento en e I mismo se ro una variable de estado y la descarga de sal ida uno va
riable endgeno.

c - Relaciones funcionales: se pueden presentar baja dos formas: 10$

identidQdes y las caracterfstkas operacionales.

Las identidades son definiciones o proposiciones relativas a las componentes del

modelo; por ejemplo, se define como descarga afluente a un embolse en un pe-
rodo determinado, al volumen medio de agua por unidad de tiempo que capt
el embolse en ese perodo.

Las caroctersticos operacionales corresponden o hiptesis generalmente expresa

das o travs de ecuaciones matemticos que relacionan las variables exgenas y
de estado de un sistema con sus variables endgenas; en los procesos estocsti-
cos, los caractersticos operacionales toman la formo de funciones de densidad
de probabilidad.

d - Po rme t ros son elementos cuyo funcin es expresor (o travs de los

relaciones funcionales) cuantitativamente el proceso estudiado; pueden ser deter
minados por tentativas, como en el coso del ajuste de modelos matemticos plu=
vio-hidrolgicos o uno regin determinado, o calculados por inferencia estadrsti
ca en el caso de modelos de transformacin pluvio-hidromtrica bosodos en una
11
 regresin mltiple, o en el caso de la estimativa del valor esperado de una va-
riable. &1 general, se consideran como parmetros las series histricas utiliza-
das como variables exgenas.

1.2.4 Aplicacin de un modelo de simulacin

La aplicacin de un modelo hidrolgico de simulacin o una cuenca particular

se hace a travs de un procedimiento de verificacin donde se establecen los
valores de ciertos parmetros paro un sistema tpico particular. La verificacin
de un madelo se hace en dos etapas: la calibracin o fijacin de los valores
de los porometros y la pruebo del madelo; para los dos se requieren datos del
sistema prototipo.

La calibracin significa la realizacin de una serie de modificaciones en los po

rometros hasta obtener un ojuste o aproximacin adecuados entre las funciones -
de salida observada y calculada; de aqu se deduce, evidentemente, que la exoc:
titud del modelo no puede exceder lo de los dotos histricos. .

la evaluacin de los parmetros del madelo puede hacerse mediante cualquier

procedimiento deseado. El procedimiento ms racional (Hill, 1972) consiste en
asignar a cado coeficiente o variable un valar inicial, valores extremos y un
nmero de incrementos que cubran eso variacin. Se hace variar la primera va
riable a travs del intervalo especificado, mientras se mantienen las dems e
sus valores iniciales. Se escriben.los valores de la funcin objetivo o de salida
y se almacena el valor que produzca el mfnimo. Despus de completar el ci-
elo para lo primero variable, se toma la segundo y se repite el procedimiento.
Despus de haber variada tados los coeficientes se escoge el con junto de valores
que produjo cada mnimo local y con l se aplica el programo de simulacin;
lo funcin as obtenida se compara con lo mnimo obtenida en las pruebas pasa
das. El vector que produce la m(nima funcin objetivo se selecciona como vec=-
tor inicial pora la prximo fase, y el proceso se repite hasta encontrar un vec
tor que praduzca una razonable carrespondencia"entre las solidos observada y-
calculada.

la calibrocin del modelo se alcanza con un proceso de ajuste que establece

los pormetros del madelo pora un coniunto particular de datos, de una unidad
hidrolgica dada. La prueba del madela requiere el uso de un segundo conjunto
de datos de lo misma unidad hidrolgica, con el cual se aplica el modelo y se
determina el grado de ajuste entre las datos observadas y calculados.

Es muy posible que estas pruebas indiquen la necesidad de efectuar ajustes, ya

seo en los datos de entrodo o en lo estructuro mismo del modelo. Cuando se ha
ya realizado una verificacin adecuado, el modelo est listo para ser operado. -

12
El anlisis de sensibilidad se llevo a coba mediante la modificacin del valor de
una voriable, mientras se mantienen constantes las dems, para observor los cam
bias en las funciones de salida. Si cambios pequeos en un parmetro produce
grondes cambios en la respuesto del modelo, se dice que el sistema es altamente
sensib le a ese parmetro.

En esta forma es posible determinar lo importancia relativa de los parmetros y

los funciones de entrada con respecto o la respuesta del sistema.

13
 CAPITULO 11

CLASIFICACION DE MODELOS HIDROLOGiCOS

2.1 MODELOS DETERMINISTICOS y ESTOCASTlCOS

Cuando los dos palabras estocstico y determinfstico se uson en hidrologfa y en

recu/'5O$ h(dricos, y como ellas normalmente separan los dos enfoques ms impor
tantes de la hidrologfa, es conveniente una reflexin sobre la larga controver
sia histrico en ffsica y filosoffa entre determinismo y probobili$ll1a, que com"en
za desde la antiguo Grecia o an antes de esa civilizacin.

Es suficiente el estudio de la historia de lo ffsica y la astroffsica para hacer

evidente eso controversia. los representates de ambas escuelas, la detenninfsfi
ca y la estocst ica, han estado batallando entre sf desde hace cerca de dos 0-:
tres mil aos; ambas escuelas existfan en la ciencia y lo fi losoffo griegas. Uno
era la del mundo predeterminado, con tomos concebidos como pequei'las partrc:u
las muy bien definidos; la otra consideraba que el caos aparente demando que -
los fenmenos naturales estn sujetos a las leyes del azar y estima que el deter-
mini$ll1O es slo un aspecto de lo qLe sucede en lo Naturaleza. Desde las9l'ie
gos estas dos escue las se han remplazado y suplementado mutuamente. -

El concepto original de la ffsica de pequeas part(culas discretas hasta el des-

cubrimiento de la molcula, enfatiz el punto de vista determinfstico para 10-
explicacin de la estructura de la materia. Cuando se descubri el movimien-
to Browniano, se demostr la existencia de un movimiento catico de las mol-
culos de Ifquidos y gases para pequeas escalas de tiempo y espacio; entonces
la escala de espacio se redujo an ms hasta el tamao del tomo. El concep
to original de ffsica de partfculas determinfsticos se restableci y mantuvo has":
to el descubrimiento del quantum, el cual revivi nuevamente el concepto de lo
teorra estocstica sobre lo composicin de lo materia. Cuando se penetr al in
terior del tomo, se volvi a tener en cuenta el concepto determinfstic:o para ::
explicar $U estructura y la de la materia. Pero cuando los aceleradores de par
t(culos fraccionaron los ncelos atmicos, se encontr un movimiento catico -
dentro de ellos y surgi una vez ms un tipo de planteamiento probabilrstico
para lo explicacin de la estructura de la materia.

15
 Estos contrnuos cambios tambin se extendieron a la mecnica de flurdos; al mi
rar los procesos del flujo laminar a la luz del movimiento Browniano, las Irneas-
de flujo san consideradas camo procesos determinfstic:os, propiedad que vino a
convertirse en la base paro el determinfsmo de la mecnica de flurdos; sin em-
bargo, si se amplfa considerablemente lo escala del tiempa, el movimiento Brow
niano mostrar los efectos aleatorios de las Ifneas de corriente. Cuando la velo
cidad del agua se aumenta suficientemente, el flujo se convierte. en turbulento,-
como un proceso estocstico gobernado por las leyes del acaso. Cambiando las
escalas de tiempo y espacio, ocurren virojes alternativos de la concepcin deter
mlnrstica a la estocstica, pora explicar los procesos de la mecnica de' flufdos:-

De manera similar, si vamos de la escala del tamao de la tierra a una escala

del orden de magnitud del sistema solar, se observa que los componentes de este
sistema se mueven de una manera predeterminada; as( fu creado y mantenido el
punto de vista determinfstica de la astronomfa, en el cual se basan mudtos prin-
cipios ffsicas y filosficos resultantes del movimiento en el sistema salar. Estos
principios pueden haber dada origen a la creencia de que todo est predetermina
do en la Naturaleza, o que sus leyes son siempre determinrsticos. Pero cuando-
se penetr en el espacio profundo y se estudi el nmero de estrellas contenido
en un espacio dado, a uno esc"ala muchrsimo mayar, fu muy diffcil -si no mpo
sible- describir determinfsticamente lo distribucin de los cuerpos celestes. Se re-- '
curri entances o los principios estocsticos poro establecer la probabilidad de
encontrar cuntas estrellas pueden hallarse en un espocio dado y en un tiempo
dada.

La anterior exposicin sirve pora demostror cmo es posible encontrar fenmenos

naturales en hldrologra que pueden ser estudiados y descritos determinfsticamente,
siempre que puedon simplificarse, y cmo hay fenmenos que no pueden serla;-
en consecuencia, deben usarse mtodos estocsticos y determinrsticos de investi-
gacin, anlisis y descripcin, paro un moyor entendimiento de los procesos hi-
drolgicos naturales.

Existen, pues, dos criterios bsicos para simular sistemos hidrolgicos: el primero
a travs del uso de un "modelo determinfstico" cuya respuesta seo equivalente a
la del sistema frsico. Un modelo es determinrstico cuando cualquiera que sea el
valor de la variable tiempo, la respuesta a una entrada dado es siempre la mis
ma, pora un mismo estado inicial del sistema.

&! un modelo determinrstico el problema se reduce bsicamente a la determina-

cin y aiuste de 'los parmetros que describen el sistema; por tal motivo estos mo
delos son llamados tambin "paramtricos". Debe notarse sin embargo, que en -
real dad no es posible Formular un sistema hidrolgico natural en trminos estrlc
tameryte determinrsticos pa cuanto:

a - Existe variabilidad de los sistemas hidrolgicos en el tiempo, debido o los -

cambios introducidos por el hambre directa o indirectamente, y o los proce-
16
 $OSnaturales de erosin, cambios climticos y otros que constituyen la evolu
cin geomorfolgica de la tierra;

b - Existe incerteza con respecto a las magnitudes y distribucin espacial y tem

poral de las entradas y salidas de los sistemas hidrolgicos y con respecto (j
los estados y propiedades de sus elementos interiores;

c - Existen dificultades en la formulacin matemtica de los complejos procesos

no lineales de transferencia de maso y energra que constituyen el ciclo hi-
drolgico. La variabilidad temporal, o no estacionalidad, es de importancia
en el estudio de sistemas hidrolgicos puesto que un sistema en proceso de
cambio no puede producir relaciones invariables entrada-salida, debido a que
tales sistemas tienen "memoria", si se entiende por este trmino el efectoque
las entradas anteriores tienen sobre la operacin presente del proceso; est
asociada con todas las etapas intermedios cuya influencia es siempre trans-
portada hacia adelante en el tiempo.

Por lo tanto, mientras estos modelos operan de un modo determinrstico sobre los
datos de entrada que se les suministran, podemos a lo sumo hacer afirmaciones
probabilfsticas respecto a las sal idas correspondientes. En otras palabras, el he-
cho de que el modela sea determinrstico no significa que el comportamiento del
sistema hidrolgico natural por l representado lo sea.

El segundo gran camino de simulacin es la determinacin de los pan:imetros esta

drsticas que describen la respuesta de I sistema, para ser usados posteriormente en-
la generacin de series de datos estadrsticamente indistinguibles de las series ob
servadas; esta es la llamada simulacin estocstica.

Un sistema se puede llamar probabilrstico cuando, a igualdad de estado inicial,

la respuesta ante una misma entrada es aleatorio; por lo tanto, para una misma
entrada pueden obtenerse salidas diferentes, siguiendo una determinada distribu-
cin de probabilidades. El sistema ser llamado estocstico si, adems, existe una
relacin secuencial (proceso en cadena) entre las OCUrrencias del fenmeno.

Un modelo estocstico es, pues, un tipo de modelo probabi/(stico.

En los modelos determinrstico se pretende simular de manera contrnuo en el tiem

po uno sucesin de eventos hidrolgicos, y lo comprobacin de la validez del mo
delo se hace por comparacin con la sucesin de eventos observados; con los mo
delos estocsticos no se busca uno simulacin contrnua de eventos; su objetivo es,
en general, la generacin de series no observadas de igual probabilidad de ocu-
rrencia que las observados.

Cada uno de los criterios expuestos tiene ventaias y limitaciones que los hacen
aconsejables paro distintos tipos de problema. Los modelos paromtricos general
mente requieren datos de entrada con intervalos de tiempo cortos y sintetizan :'
17
 bien respuestas para incrementos del mismo orden. (datOs horarios o diarios), ro
zn par la cual se usan ampliamente para simular eventos hidrolgicos con inter
valos cortos de tiempo. La simulacin estocstica se usa fundamentalmente paro
predicciones can intervalos largos de tiempo, donde lo que interesa no es la res
puesta instantnea del sistema sino los valares caracterfsticos; como las modelos
estocsticos trabajan con informacin estadfstica, es muy difi'ci I modelar con ellos
fenmenos instantneos a de intervalos cortos de tiempo; par eso, la simulacin
estocstica es usada principalmente para prOpsitos de planeamiento, para gene-
rar series "igualmente probables" de caudales medios mensuales, por ejemplo.

2.2 MODELOS ANALlTlCOS y EMPIRICOS

Un modelo hidrolgico es analrtico cuando est basada en las ecuaciones bsicas

de la hidrodinmica; tambin es llamado "de sistema distribuido". &, estos mo-
delos el sistema se trota como un con junto de reas distribuidas en el espacia, y
se simula el comportamiento de las diferentes partes; son modelos en los cuales
se analiza y se trata de seguir el proceso interno del sistema. El hecho de can
siderar las coordenadas espaciales adems del tiempo, implica que en estos mo':-
delos se trabaja con ms de una variable independiente. Eh el caso de que los
relaciones den origen a diferenciales totales funciones de varias variables simul
tneamente, se expresan o travs de un sistema de ecuaciones diferenciales a efe
rivadas parciales.

Las modelos empfricos o "de sistema bloque", estn detArminados, en trminos ge

nerales, slo por los datos de entrada y salida disponibles (lluvias y descargas e
una cuenca, par ejemplo); en ellos las coordenadas de posicin carecen de impar
tancia y todas las partes del sistema simulado se tratan como si estuvieran locali-
zados en un mismo punto del espacio. Tambin suelen llamarse "de caja negra"-
porque no se siguen paso a paso los mecanismos que involucra el proceso interno
del sistema. El modelo est constituido por un "bloque" que acepta datos de en
trada y proparciona datos de salida, sin que interese conocer lo que sucede de
tro del bloque; s/n embarga, el grado de empirismo puede conducir a distintos ti
pos de "bloques"; en efecto, el bloque est constituido par el coniunto de algo
ritmos que permiten sintetizar las datos de salida en funcin de 10$ datos de en-:
trada; estos algoritmos pueden ser de cualquier fndale y tienen, en general, como
nica limitante la ecuacin de continuidad. Por eiemplo, un modelo pluviohidro
mtrico que tiene como datos de entrada lluvias y cuyo "caja" est constituida-
por un algoritmo que incluya coeficientes de prdida, sin preocuparse par espe-
cificar el proceso que conduce a las prdidas por infiltracin, evaporacin desde
el suelo, evaporacin desde los almacenamientos superficiales transitarios, etc.,
podra llamarse con propiedad "de caja negra"; no asf si los algoritmos permiten
individualizar coda uno de los procesos.

En Hidrologra, a causa de que muchos procesos no san an suficientemente co-

nocidos, todos los modelos tienen algo de empirismo; incluso en las analfticos o
18
 hidrodinmicos se puede llegar o un punto en donde es preciso acudir 01 empiris
mo y originar as, en lo prctico, modelos denominados semi-empricos.

2.3 MODELOS A EMBALSES Y CON OPERADOR PLI'VIOHIDROLOGICO OPH)

Se clasifican as segn el mecanismo de transformacin de escorrenta en desear

90. En los primeros, eso transformacin se hace a travs de uno serie de fases,-
cada una de las cuales simula el paso del agua a travs de un embolse, y est
especificado por dos ecuaciones: lo de continuidad y uno ecuacin de almacena
mienta. En estos modelos lo escorrenta puede dividirse hasta en cuatro compo":
nentes; sin embargo, lo subdivisin en ms de dos componentes no siempre est
justificado. En los modelos con OPH la transformacin de escorrenta en desear
90 de salido se hace a travs de un operador de transformacin pluviohidromtrl
ca tal como un hidragroma unitario. En ellos generalmente se consideran slo-
dos componentes del flujo, con un operador de transformacin para cado fose.

2.4 MODElOS LINEALES Y NO LINEALES

Un modelo ser lineal o no segn el carcter de la funcin de transformacin

usada; en trminos generales, un operador es lineal cuando, actuando sobre una
determinada entrada, produce una salida proporcional a ello; si un operador de
transFormacin es independiente de lo entrada y la solido podemos decir que es
lineal; el hidrograma unitario, por ejemplo, es un operador lineal. En otras po
labras, un modelo lineal es aquel representado. por ecuaciones lineales; si se tia
ta de modelos de optimizacin de recursos hfdricos, todas las restricciones, o por
te de ellos, y la funcin objetivo, son lineales; por ejemplo: -

ecuaciones lineales: Y ~ 5X l + 6X2 + 7X3

ecuaciones no lineales Y ~ 5X~ + 6X 03

y = Log XI

2.5 MODELOS DESCRIPTIVOS Y CONCEPTUALES

Es otro manera de enfocar lo clasificacin de modelos matemticos. El modelo

descriptivo usualm~nte se disea para reproducir fenmp.nos observados mientros
que el modelo conceptual se construye coro dilucidar la esencia de lo teora
que puede interpretor el fenmeno. En la prctica, la distincin entre los dos
clases na es muy definida y a veces, en efecto, suelen combinarse. Sin embar
go, los extremos de los dos clases de modelos son claromente reconocibles. -
En un modelo descriptivo d~ un sistemo hidrolgico, la base ms comn pora su

19
 formulacin es la ecuacin de continuidad o balance de masas. El modelo con-
ceptual se fundamenta en un cierto concepto o teora; por ejemplo, hoy modelos
estocsticos que estn basados en la teora de los procesos estocsticos.

Si denominamos por x (t) las variables de entrado y sal ida de un sistema en el

tiempo t, un modelo matemtico podra definirse por
d~ Jy.. 'O'y
f [ 1<0, "'t-(O''dT' j" 'i{T
'1)

donde g, 7 G, san pormetros estimados a portir de los datos.

Cuando el sistema en estudio es muy campleio, puede ser ventajoso adoptar una
forma ms simple f * (.) de lo funcin f l.) de lo ecuacin 11), y expresar la
"falta de ajuste" del modelo mediante un residuo o error; entonces, lo ecuacin
(1), remplazando las derivados porcioles por diferencias finitas, se convierte en
(2),
=0
2)

donde Et es el residuo en ~I tiempo. La acertada escogencia de la funcin.f* (.)

es lo ms importante en lo definicin del modelo.

Clorke (ref. 19) ha propuesto ogrupor los diferentes modelos presentados en la li-
teratura hidrolgica en cuatro grandes categoras:

Estocstico Conceptual lEC)

Estocstico Emprico (EE)
Determinstico - Conceptual (OC)
Determinstico - Emprico 'DE)

Para modelar series hidrolgicas es necesario usar modelos motemticos que es-
tn de acuerdo can los caractersticas del problema que va
ser resuelto. EI-
escurrimiento en cursos de agua es un proceso estocstico, lo cual implica una
total dependencia temporal de las valores aleatorios (descargas\. Si este porme
tro (tiempo) en una etapo dada asume todos los valores reoles, tendremos un prO
ceso estocstico contnuo (descarga instantnea); pero si asume valores integrares
que son mltiplos de alguna duracin de tiempo (ao, mes, etc.), entonces esta
remos ante una secuencia aleatoria o un proceso aleatorio discreto (descargo me::
dio anual, mensual, etc.). Por tonto, si alguna de las variables 1-~, 'ff.., E. en
la ecuacin (2) es aleatoria, con distribucin de probabilidades, el' modelo ser
Estocstico, me jor que estadrstico, poro enfatizar la dependencia temporal de los
variables. Si todas las variables en la ecuacin 2) estn libres de variaciones
aleatorias, o sea que ninguno tiene distribucin de probabilidades, el modelo-
ser considerado Detenninrstico. En otras palabras, si se toma en cuenta la po
sibilidad de ocurrencia de los variables y el concepto de probabilidad se introd

20
ce en la formulacin del modelo, el modelo ser estocstico o probabilstico.

El modelo, por otra porte, ser conceptual o emprico, segn que la formo de
la funcin f* (.) en la ecuacin (2) sea, o na, sugerida de acuerdo con los pro
cesos fsicos que acton sobre las variables de entrada para converti rlas en las
de solida.

Describir un modelo como conceptual requiere mucho menOs justificacin que ca

lificarlo de -estocstico; sin embargo, los propiedades estocsticas tanto de la lI
vio (variable de entradal, como de la descorga (salida), han sido ampliamente-
recanocidas por varios autores: Sutcliffe 1] 966) I Hutchinson (1970,1972), Dickinson
(1967).

Como ejemplos de modelos dentro de cada uno de los grupos, podemos citar:

Modelos EC - a) Modelo de Dawdy - Donnell (1965)

b) Modelo de Nash-Suteliffe (1969)
e) Modelo SSARR (1968)

Modelos EE - a) Mode los de regresin descritos por Guillot (1973)

b) Modelo'de Thomas - Fiering
cl Modelo de Mandelbrat (1971)
d) Donnel (J 973)
e) Hidrograma Unitario Instantneo, con ardenodas
estimadas por mnimos cuadrodos (Snyder, 1955,
Donnell, 1966, Jenkins y Watts, 1968)

Modelos De - al Modelo de Freeze (J 971),

b) Modelo hidrulico de Wooding (1966),
el Ecuacin de Laplace para flujo variado en acuiTe-
ros no conFinados (Eagleson, 1970)

Modelos DE - a) Hidrograma Unitario Instantneo, estimado par ex-

pansiones de series armnicas (Donnefl, 1961)
b) Amorocho y Orl ob (1961)
el Kulandaiswamy (1973.

21
 CAPITULO ;[1

\\ODELOS DET ER/'v'IN IST[C es

3.1 ANTECEDENTES

Los primeros hidrlogos q..'e comenzaron o obsprvor los fenmenos naturales a fines
del siglo XIX y comienzos del XX, descubrieron inmediatamente en las observacio
nes de precipitacin, evaporacin y escurrimiento, un significante "ruido" o "es:.
tocasticidod" en las series temporales observados; sin embargo, esto aleatoriedad
de los fenmenos se asimil a algn componente sistemtico de la naturaleza en
la forma de periodicidades en las series; '!ste aspecto di un gran impulso al an
lisis de las series hidrolgicas desde el punto de vista determinstico, en busca-
de regularidades peridicas.

Una segundo ero de determinismo est ~n desarrollo actualmente, en la forma-

de investigaciones de funciones de respuesta determinstico entre variables hidro
lgicas aleatorias, se observo uno variable aleatoria y se busco una relacin fun
donal matemtica determinstico o experimental entre esa variable y otro u oh-as
variables. Un eiemplo tpico os lo relacin precipitacin -escurrimiento en la -
cual lo precipitacin de entroda o una cuenco se convierte en lo salido (escurri
miento) o trovs de funciones de respuesto fija a modelos determinsticos emprri
COSi este tipo de respuesto hidrolgico o reaccin del sistemo, es la base de la
actual era de desarrollo de los procesos hidrolgicos desde este punto de visto.

Sin embargo, algunos conceptos determinsticos corrientemente usados influyen de

manero significativo en el anlisis de los procesos estocsticos naturoles y ocosio
non dificultades de varios tipos.

En la literatura hidrolgica se encuentro a menudo el concepto de que existe uno

especie de rastro fsico en los eventos hidrolgicos (porte determinsticoJ, as co-
mo tambin algo de "ruido" '.porte estocstico). Con el tiempo, lo relacin en-
tre el rastro explicado por los [ey"s ffsicas, al "ruido" no explicado, tiende a -
aumentari en otros polabras, el contnuo progreso en el conocimiento fsico de los
procesos hace que la porte desconocido vaya disminuyendo, "" lo medido en que
se profundizo en el mecanismo fsico que los controla; esto significa qu~ la expli
cocin determinstico de las relaciones causa -efecto, ser el resultado final de-
23
 las investigaciones, mientras que la explicacin estocstico de los fenmenos oleo
torios notun:lles desaparecer con el tiempo como uno propiedad indeseable y tem-
poral de estos fenmenos. -

Deben analizarse con especial cuidado los consideraciones emitidas en el prrafo

anterior, pues stas san, probablemente, uno de las equivocaciones bsicas que
se encuentran en las aproximaciones filosficas a muchos problemas ffsicas en ge
neral o geofsicos e hidrolgicos en particular. -

Por ejemplo, en la determinacin de la precipitacin efectiva como eSCUrrImIen-

to superficial en lo aplicacin de un hidrograma unitario, se omite la aleatorie-
dad de la evaporacin, infiltrocin, parmetros de la cuenca y otros factores que
intervienen; para aplicar uno aproximacin determinrstica, se trota generalmente
de esquivar el carcter aleatorio de los procesas naturales. Para aclarar esto, es
suficiente contemplar el coso de la rugosidad de las superficies. Es un hechosim
pie que la rugosidad de un terreno, o de una cuenca, es un factor decisivo en -
la respuesta de una cuenca. Si se conservan iguales los dems factores geom-
tricos de la cuenca, es fcil demostrar que la rugosidad tiene un efecto muy gron
de sobre el hidrograma unitario. Entre ms lisa sea la superficie y menor seo 10-
profundidad de la capa de flujo superficial, menor ser la retencin; entre menor
sea la rugosidad del lecho Fluvial y moyor la velocidad del flujo, menor ser la
profundidad y por lo tanto menor el almacenamiento y ms rpida 'a salida, y
el tiempo de concentrocin y los efectos del almacenamiento sern pequellos. -
Cuando hay vegetacin (pastos) desarrollada en el rea, mayor ser la resisten-
cia al flujo y mayor la capa de escurrimiento superficial. Esto cambiar el hi-
drograma unitario. La rugosidad de la cuenca se puede modificar en forma pe-
ridico y tambin aleatoriamente a causa de los contfnuos cambios de la vege-
tacin debidos al clima y a los muchos factores aleatorias involucrodos. El con
cepto del hidrograma unitario constante para uno cuenca es, pues, una concep
ci" artificial para aplicar una aproximacin determinrstico a un procesa hidro":
lgica bsicamente muy complejo; como se mostr oqur, la rugosidad cambiante
presenta dificultades en la aplicacin del concepto de uno respuesta unitaria -
constante bajo condiciones naturales.

No obstante las premisas anotadas anteriormente, los modelos determinrsticos -

orientados dentro de los principios establecidos, y aplicados racionalmente a los
eventos aleatorios hidrolgicos, son uno poderosa arma de amplia utilizacin en
la actualidad que ha permitido resolver numerosos problemas prcticos y profundi
zar notablemente no slo en el conocimiento de los procesos ffsicos, sino tambi"en
en su descripcin matemtica. A manero de ejemplo, se presentan 10$ fundamen
tos de tres modelO$ que han sido probados y usados extensamente con resultados-
satisfactoriO$

24
3.2 SIMULACION DE ESCURRIMIENT O .0. TRAVES DE UNA ZONA DE INUNDA
CION

(Modelo SOGREAH-UNESCO, usado en el delta del Mekong)

3.2.1 Caractersticas del problema

A. Genera I i dades

Los consideraciones siguientes se refieren al coso especrFico de lo simulacin

del escurrimiento en zonas inundables, donde 10$ modelos del tipo propaga-
cin sufren innumerables compl icaciones adicionales, ya que entran en jue-
go fenmenos fsicos e hidrulicos que aumentan el nmero de parmetros que
deben tenerse en cuenta en el made lo .

la teora fu desarrollada en el modelo SOGREAH-UNESCO paro ser aplica

da en el delta del ro Mekong, y una adaptacin se us con buenos resulta
dos en el Centro de Hidrologa Aplicada de Parto Alegre, Brasil, para l0si
mulacin del ro CaL

Se presento aqu la solucin matemticJ dada 01 problema, en vista de lo

similitud de condiciones presentado por e I ro Magdalena en su zona inunda
ble, como uno sugerencia del enfoque que podra darse para la formulacin-
de un modelo matemtico capaz de simular aceptablemenle el escurrimiento
en esa zona, donde la abundancia de cinagas y lagunas profusamente inter
conectadas, con canales de flujo reversible segn el estado del rfo, dificuf=
to grandemente los onlisis tradicionales.

Los siguientes puntos resumen los caractersticos que son comunes o este ti-
po de problema :

o) Escurrimiento multidimensional (o multidirecciono l), causado como ya se

anot, por lo zona logunorio que acta como amortiguador de las crecien
tes.

b) Pequeas velocidades, ya que el flujo en su mayora se efecto por los

canales de drenaje principales, comportndose lo zona inundada como -
una serie de embalses de almacenamiento.

cl Tres variables que entran en juega: dos de espacio Ilos coordenadas X

e y de un punto cualquiero de lo cuenca), y e l tiempo t

d) Condiciones lmites conocidos aguas arriba yaguas abajo de la zona si

mulada.
25
 e) Efecto posible de marea, ya que las zonas inundables normalmente (no
siempre) se encuentran cer ca de lo desembocadura. Esto impl ica que
los trmin05 de inercia en los ecuaciones diferenciales que rigen los
fenmenos, no pueden despreciarse.

f) Puede presentarse otra dificultad adicional: en ocasiones, el efecto de

marea se comunica de manera compleja o travs de una red de cana-
les, o seo que lo zona de inundacin que une el rfo al mar no se pre
senta como una superficie contnua. -

B. Configuracin de la cuenco

En la figura 1 se muestra un esquema de lo cuenca que va a ser simulada,

lo cual se subdivide en una serie de subreas con base en lo topografa y
en la configuracin gener~1 del escurrimiento (canales, ros, carreteras, zo-
nas de vegetacin que pueden impedir la propagacin del escurrimiento, di-
ques y defensas yo construidas, etc.) (figuro 2) .

Cado una de estas subreas ser considerado como un embalse que intercam
bia agua con las adyacentes, y sigue lo ecuacin de continuidad y las leyes
de descarga entre estos embalses. Esto solucin permite analizar el escurri-
miento como unidimensional solamente, entre cada por de embalses (aunque
son dos sentidos posibles).

3.2.2 Descripcin del modelo

Evidentemente, el modelo estaro fundamentado en los leyes de intercambio entre

cada par de embalses adyacentes, de los cuales reconoCf'mos lo existencia de dos
clases: de tipa flwial y no fluvial.

A. Leyes entre embolses fluviales

Seron tratadas en esta forma las uniones entre embalses comunicados directo
mente, sin ninguno barrera fsica de importancia, tales como trechos de ros
y canales, unin de canal y ro, etc., ocasiones ~n los que rigen los leyes
del movimiento uniforme en canales, obedeciendo por tonto lo ecuacin de
Strickler

Qk = kstr Ak Rk 2/3 }/2 (1)

donde
Qik = caudal que pasa del embolse a! k

26
\

"
~
"

- ---- -- ..........

,,
.
... .
~
Q
e
/
......~/:::::==::----. - ... "
e
;;
e
-o
~

E
CI< u i
::>
c:> " -
~
~

\
,
... ",'"
:;:~
e
O e
\ U .-
\
\
\
\
\
\
\ /
\ /
\
\
\

27
 kstr = coeficiente de Strickler

Aik = rea mojado de lo seccin de poso entre dos embolses

Rik = radio hidrulico correspondiente a la seccin Aik

J = pendiente de lo Irneo de aguo (= (Zk -Z i ) /L ik)

El valor de k str est defin ido por

21. 1

'lf;'
siendo e lo rugosidad absoluta de las pared.. s del canal; existen tablas que
dan el valor de k str en funcin de la naturaleza del material constitutivo
de las paredes :

Canales con revestimiento de concreto K = 53 a 57

Conales con revestimiento de concreto liso K '" 80 a 90
Galerras excavados en roco K = 250 40
Canales antiguos con vegetacin K = 43 a 52
Canales de tierra K = 33 a 40
Canales de arena K = 50 Q 90

R(os y arroyos :
de fondo rocoso, rugoso . K = 20
de fondo medianamente rugoso K=20a28

Sin embargo, dada la dificultad para estimar su valor exactamente a causa

del desconocimiento de los carocter(sticas frsicas de las secciones de paso,
asr como de sus variaciones, kstr es tratado como un parmetro de aiuste
del modelo.

En esta formo, para cada seccin de paso (de tipo fluvial) deben conocerse
los siguientes datos:

a) Un volar de "str
b) Un valor de Lik (distancia entre los puntos centrales de dos subreas
o embalses adyacentes)
c) Caracterrsticas geomtricos de la seccin de paso, con informacin sobre:
Z altura del aguo
Aik - rea mojada
28
 Rik - radia hidrulico correspondiente

d) Nmero de subdivisiones en cada seccin de ooso, er. el caso "" que

sta puede subdividirse en otras ms sencillas ver figura 3;.

B. Leyes entre embalses no fluviales

Se incluir aqu lombi"n ,,1 caso de intercamio entre un emoolse fluvial y

uno no fluvial; estomas anTe el caso de flua sobre un vertedero, por locuol
podremos denominar los secciones como de tipo v"rtedero. ver figura fi

Para oolicar las ecuaciones del vertedero, distinguimos dos coses: v"rtedera
ahogado y no ahogado, con las siguientes ~otaciones;

Seon f A Y 4D dos variables qU'l dependen de las caractersticas ffsicas

del vertedera, para caso ahogado y no ahogado respectivamente, tales que

iD = )"1 e "9
h = "'2 D' 29

donde b longitud til del vertedero

.1' = coeficiente de descargo del vertedero

Se ha estimado que ,r-o

= :J, 386 ~A aproximadamente, ya que el coefi-
ciente de descargo vorf" ,pgn la condicin hidrulico de la secci" de ~
so.
En lo tabla siguiente se resumen los ecuaciones pora coda uno de los dos
cosos mencionados, con el siguiente significado:

'-k = altura de la lmina en el embalse que cede agua

ZI = altura de la lmina er- el embalse que recibe agua

Zs = altura de la solera

REGIMEt'-' Ii CONDICION FORMULA DE CALCULO

3/2
-,. - -s- +D ~Zk -Zs)
~ ~

No ahogado
.'
" 1Zk - -s) Qik =
3
Ahogado
I
-i
-
- - -s ~
- - --
! :'k '51 I
Qik = ~A
,Z _ L
SJ Zk -Li

29
 4

SUBDIVISION DE UNA SECCION DE PASO

FIGURA 3

EIftbCII i
T
I 1
Zk I i
lo
Zi

I,

UNIONES TIPO VERTEDERO

FIGURA 4

Ahi

_1st i

l. Si ..1

FIGURA 5

30
 As, los parmetros necesarios 01 modelo, poro este clculo, seron

o - Nmero de vertederos en cado seccin de paso

b - Paro cado vertedero: valores de ;; y b;./ ser tomado como parme-
tro de ajuste.

3.2.3 Formulacin motemtica del modelo

A. Hiptesis previos

Para lo formulacin de las ecuaciones que definen el comportamiento hidr~

!ico del sistema, aceptamos dos hiptesis previas:

a - La superficie del aguo se mantiene horizontal; suposicin razonable pues

to que la velocidad del flujo en las zonas inundados es muy reducida.

b - Relacin biunvoca entre caudal y alturo; condicin derivada de la an-

terior.

El problema podr ser resuelto con el uso de dos ecuaciones: la de continui

dad y la de descargos, ecuaciones que conformarn un sistema parablico;"
puesto que no se tendrn en cuenta los trminos de inercia (no se considera
efecto de mareo).

B. Ecuacin de continuidad
Haciendo referencia o la figura 5, donde se ha representado un embolse i,
si .1 hi es el aumento de altura en el interv<:rlo lH, tendremos que

s, llhi '= ,
p. + (2)

donde
es lo precipitacin cada en lo superficie del embolse i,

es la descargo media intercambiada con todos los embolses ad

yacentes en el perodo t; t; entonces,
t t + 4 t
Q ik = (1 - IX. ) Qik + oQik 13)

donde los ndices t y t + .i t sealan que Qik es calculado en los

instantes t y t + t respectivamente, yo( es un coeficiente odimensional
de ponderacin que depende de lo manera como varan los descargos duran-
te el intervalo {j t; as, o( = _1_ significo que tomamos la descarga media
2
31
 entre la del comienzo y la del final del intervalo,

&ttonces, ....,
() Qik
+ ( a hi ) t fJ. hi +

~Qik t
+ ( 'J hk ) ~ hk ] (4)

Tomando las derivadas de primer orden (desarrollo en 5eries de potencias),

Al combinar (2) y (4) Y dividir por (\ t obtendremos

Si Pi
.1hi - - = 2 Q'k
At LIt k '
remplazando (4) Y aceptando r/. igual poro todos los embalses,

~ Q'k +Z:- a Qk
- Si h',- Pi --
o{ 6 t ",,ot
1
o<. k ' k
" ~i
Ll hi +

d Qik
+~ .1 hk
k 'a hk

es decir

--,,--
c( <1 t
~ d Qik hk = - -1 [ --'
Si ] Cl hi + o:::-
k
P' +z"Q'k
a hk oc At k '
1 (5)

La ecuacin (5) ser una de las qu-. conforman el sistema parablico; la

segunda, ser la expresin de descarga del vertedero, (convencin de siS
nos: Qik = - Qki ),

C. EcuacionRs de descargas y sus derivadas

Unin de tipo fluvial

Como se dijo atrs, se opl ica la ecuacin de Strickler en lo formo

Q'k
, = E I'k t * IoQ ni
k sr "-, - hk /I '1/2 (6)

32
donde
hi , hk son las alturas de los planos de aguo en lo porte cenlral
de cada sub rea i y k o

:. ik es el signo de (hi - hk): ser positivo si hoy flujo de k

hocio i y negativo en el coso contrario.

FQ es uno funcin que depO~dQ de la configuracin geomtrica

de los pores de embalses, dada por

R 2/3
FQ =
donde
Lik es lo distancio entre los puntos centrales de los sub reas
i y k.

k str es el coeficiente de Strickler

La tuncin FQ se puede obtener por interpolacin a portir de

h + hk o

a - una tabla de valores hik = --,-'-.,-_ _

2

b - lo tabla correspondiente de FQ en funcin de hik; si tenemos

Q'ik =
Qik
entonces
() Q ik
= ksIr aQ'ik
kslr el hi 'O hi

Segn lo direccin del flujo, son posibles das casos:

Caso lo Cuando hi >hk, entonces 0k = - 1 Y tendremos, derivando

(6) : . I

OQ'ik l FQ hk/ 1/2 __

1 FG /hi _ h / - 1/2
=- "'b hi
'hi -
2
k
" hi
(7)

'O Q'ik (1 ~C 1/2 ~Q ./L hk.'' -1/2

= + , ni
o _

O hk "O hk

Coso 20 Cuando \ > hi , por lo cual E ik = 1, tendremos:

33
OQ'lk d FQ /h -h /1/2 _ -2 FQ /hk -hi /-1/2
k
'a hi 1} h

(8)

/h
k -hi /1/2 +.2.. FQ / hk - h / -112
2

hi + hk
teniendo en cuenta que FQ = f (h ik ) y h
ik = ~-2-';'':'';'

vemos que

o FQ
=
OFa
-* =
1 dFQ (pues ~ ~~ =-+) (9)
2 dI. k

y 1 dFQ
= =
2
(10)
d hik

remplazando (9) Y (1 O} en (7) y (8) tendremos:

o Q'ik
'O hi
= --2 1 dFQ
dhik
11.. _
. I
h /1/2
k -2' FQfh-h~-1/2
1

(7' \
Q Q'k
'O hk
:: .. - dFQ
/ hi - hk, ;1/2 +2"
1 FQ JI,. - h ;-1/2
'1 k
2 dh k

OQ'ik
= -2I dFQ
1\ - h / j/2 _ .2.. FQ fhk - 1..1-
1/ 2
~ hI dhik '2 I

(8' )
"O Q'ilo: dFQ ; 1/2
o hk
= 1
2"dIiik
/ hk - hi + -21 FQfhk - ni 1 -1/2
Podemos ver que los rndi ces y k pueden intercambiarse en /h i - hl/
si llamamos

a la expresin FQ hk - hi - 1 12
I
2
1 dFQ I hk - hi /1/2
y 11 a la expresin
2 dhik

lograremos uno solo expresin, cualquiera que sea el caso:

-o Q'ik
'H
~
- 1 + E
ik
1I
I

(10\
O Q'ik
~
I + t ik * I1
"O hk

Unin de tipo vertedero

Con referencia a la notocin indicada en la figura 6, distinguimos 3 casos:

Caso 1. Cuando hm - h <

DHTOL, donde DHTOL, es un valor de tole-
rancia fijado adecuadamente, o bien hm < hs simplemente; entonces,

o (ausencia de flujo)

~ ~ O

Coso 2. Cuando (h - hsl < 32 (h m - hs ), tendremos reg men no aho-

godo, por lo cual,
:..
'k rD . (lJl

de acuerdo con lo expuesto antes, llamando A o lo derivada de Qik

en (11), tendremos que

A ~ 3 3/2
2

Si es el embalse oguas arribo,

35
 '" -A y '" o

mientros que si es el embalse aguas abajo,

= O y = A

2
Caso 3. Cuondo (hj - hs ) ~ T (hm - hs ) estamos en el caso de rgi-
men ahogodo y de acuerdo a lo expuesto en el numerol correspondiente,

y al igual que en el caso anterior, las derivadas san:

A =
2
1/2
(hm - h) -"21 1.
'fN (h - h5 ) (h
1 m-
h.) - 1/2
I

Resumiendo los considerociones anteriores, si es el emba lse oguos arri-

bo,

= - A y = - B

mientros que si es el embolse aguas abajo

= B y = A

3.2.4 Desarrollo del modelo

A. Distribucin topolgica

Paro plantear la solucin al problema, se dispone de sub reas en una can-

figurocin esquemtica de lneas (o filos) sucesivas, de tal manero que cual-
36
 NOTACION USADA
FIGURA 6

ETAPAS

61

EJEMPl.O DE DISTRIBUCION TOPOl.OGICA

FIGUR A 1

37
 quier sub rea de una lnea pueda intercambiar agua solomente.con otros sub
reas de las dos lneas adyacentes o de la misma lnea. Puede verse que
existen varios configuraciones posibles, pero como el nmero de sub reas
de coda fila da el toma'lo de uno motriz que deber invertirse, debe bus-
carse un nmero reducido de sub reas en cada lnea, de acuerdo con el ti
po de computador dspon ib le.

Tomamos como ejemplo lo distribucin topolgica mostrada en lo figura 7,

suponiendo que corresponde a las etapas de clculo 9, 10 y JI; cada cosi-
lla represento uno sub rea que intercambio aguo con sus vecinas, como in-
dican los conexiones.

B. Conformacin de los motrices

la formo general de lo ecuacin de continuidad puede escribirse como:

[z
k
- ~ 1~h
ri. LI t i
+ ~ d G ik llh _
k ( el hk k) - -
1 [Pi + Z
T"Kt k
0.1,k

donde Si represento el rea de lo sub rea i; multiplicando por r:i. y reor-

denando, obtenemos

S'1
lit Llt
l
- _' L1 h = - -P'' +. L<;;"" O'k ]
k '
(J 2\

Aplicando ahora lo ecuacin (12) o coda sub reo, obtendremos tantas ecua
cianes cuantas sub reas tengomos; as, lo 56 por eemplo, len.; comunicacion
con lo 51, 52, 57, 58 Y 59.

Sub rea 56 = 56
k = 5 1, 52, 57,58 y 59

La ecuacin (J 2\ produce :

eL ) 0 56 ,51 20 56,57
; h /j hS7 +
51 hS7

+ ri ) 056,58
OhS8 11 hS9 +J_6_-
... hS6
(0 56,51 +056,52 +

38
 Q Q
P56
+ 56,57 + 56,59) - +
t

En formo similar se puede obtener uno ecuacin para cada sub rea; el sis-
tema resultante se puede expresar en formo matricial asr:

,
h
51
1 Ih59
hS2 I rh56 h
60
,I ' i
h53
I + 11 I h57 - 111 I h61 =

h54 I 1hS3 j
lj
h2

h55

-l P56
t
+
k
Q56, k 1 (k = (51,52,57,58, 59)

,-l
!
P57
t
+2:
k
Q57,k) (k = 53,56,60,61)
I

l- [ ;:' +~
k
QS8, k 1 Ik "53,55,56,62)

Donde cada una de los motrices 1, /, I i: son los siguientes

39
 1 :

Q 56,52
d. OQ56,51 r:f.
Oh51 ~ h52 O O O

O O r:f. ~ Q57, 53 O O
. CJ h53

aQ dQ
O O IX 58,53 O o( 58,55
o hS3 e hS5

11 :
!
t<.;: dQ56,k 556
I ~Q5,57 ~Q56,58
oh56
- dI
o{
o(
" h57 h58
k = 51,52,57,58,59

;, Q57,56
C(. dQ S7,k - II
-tl5S7 O
o( ; h57 t
el Fi S6 k = 53,56, 60, 61
~Q58 k 558
Q(
21 QS8 , 56 a(~ ' - -
O J hS8 t
o h56 k = 53, 55, 56, 62

1II

\ DI
a Q56,59
o o I o
" h59 1 i
)Q
57,60 I ~Q57 61 :
O Q(
c< ' O
\
,
'2> h60
I ti h6 1
I
I

1
; dQ
1
I oc 58,62
I
\
O
I O
i O
' el h62 I
\
40
 Para simplificar, denominaremos las matrices de la siguiente manera:

1: Matriz precedente P (corresponde a la etapa precedente a la ac-
tual en desarrollo)

11 : Matriz central C (etapa en actual desarrollo)

111: Motriz siguiente S (correspondiente a la prxima etapa)

Vector V: Matriz de los trminos independientes.

En general, si tenemos que:

ni es el nmero de sub reas de la etapa precedente

n' es el nmero de sub reas de la etapa central, y
I
nk es el nmero de sub reas de la etapa siguiente, podemos escribir la
ecuacin matricial general como

,
Ah1, 1 r .1h.1, 1 ~ hk , 1

Ah1, 2 d h1, 2 6 hk, 2

[pd x ......... + [C] x ....... I

+ 1Si x
....... '" [ Vi]
. ... . . . .......... . .......
11 h.1, nI. Ll h. hk , nk
1, ni

y las dimensiones de las matrices son

MATRIZ No. DE FILAS No. DE COLUMNAS

I
P. n.
! nI
I I I
C. n n
I
I I

n n
Si I k
V. n.
I I
41
 C. Condiciones Irmites

Estn representodos por las condiciones aguas arribo yaguas abajo, tenie~
do en cuento que las incgnitas del sistema son los h:

= VI

= V
2

~"~' .. " " ' . " . " " " " " " " .. ' . ' .. " " " ' ... '

P., /;. 1 ,- + ei D. i+l = V.I

" " .. , , , w , " " " w " " " " , , ,

= vn

Condiciones aguas arribo

Tomando como ejemplo el modelo aplicado en el rro Car, mencionado en

A, tenemos como condiciones aguas arribo 3 entradas, representadas par las
descargas que entran en (2), (3) y (4) (ver diagrama de la figura 8); pora
facilitar lo comprensin del esquema, suponemos que esos descargos provienen
de 3 sub reas ficticios (O, -1 Y -2), La ecuacin de continuidad aplicada
a la sub rea (4) proporciona:

51 L! h. , = dt ~ Q
k ik + P.1

o seo

54 ,h 4 = .1t [ 04 , O + ~,3 + Q4,8 \ ..L P,

como
j Qk d 0ik
Q. k = Qik + o( ;; h,
6 h, ... o( t. hk (4)
"
':3 hk

entonces

54 6 h4 , :1 Q4 , 3 .,.
= .1 t 1 Q4 O Q4 3 ... :1 Ll h4
. t
I I
d 1'14

42
 , ......, ,-, ,.,
t- f } ,, 0\I
.....,,, '
\ -2 '
''T'
,,
OO'

)----{S

CONDICIONES AGUAS ARRIBA

FIGURA 8

70 71 I
,
I

,
I

,,
I
I
I I
, .. 1, ,.1.,
(75) 74 176)
.... ' '.'
,
f
I
I
I
I
f I
I f
I f
,.1., , ...1."
(771 17S'
'.'
, I
,-

CONDICIONES AGUAS ABAJO

FIGURA 11

43
 ""

?Q
+ Il(
?
Q4,3
h3
.1 h3 + Q4 , 8+ ot

4,ajh
h4
4

;, Q4 , a P4
+ o(.
ha
.1 h8) +
t. t

qgrvpando y haciendo las siguientes convenciones:

" 04 a )4 -,
~ ,
el = o(
75
Q4,3
h3
tl h3 + [~ h4 --:6tJ
.1 h4

= ()(
Q4,a
SI .6 ha
\ d ha
1J
i
P4
VI = [ t
+ Q4,O + 4,3 + Q4 , 8

\ ,
podemos expresar la ecuacloo en la formo
.

quedando incluidas en VIlos condiciones de entrada al sistema (datos de

descargas de entr<:lda a la zona de inundacin l.

Condiciones aguas abajo

Al igual que para las condiciones aguas arriba, aquf se supone la existen-
cia de sub reas ficticias (o casillas), donde se conocen los niveles.

Siguiendo con el ejemplo referido antes, en la figur<:l 9 se muestr<:l el final

del diagr<:lma topolgico, para las etapas 18 y 19. los niveles en las casi-
llas 75, 76, 77 Y 78, que son conocidos, se incluyen en la ltima matriz,
V9. Las matrices correspondientes a las dos ltimas etapas serfan :

44
 IXIJ Q6~, ~r I
ah., I I
! I !
:
"
: I, I
I
i I I

T
c( ;'Q92.14 i!l.' OQ'12.7~ .1h
, I I 0",2 ~ -ri "h
o ",
.." - Q'2 76
J

eX. I
C/u 17
D~'3
,~
0(*

~Q9<.'8
oh,.

El resto de tnninos en cada matriz se pncuentra normalmente; aqur slo se

han mostrado los que incluyen las sub reas ficticias.

D. Tcnico de solucin
El problema se ha reducido ahora a la solucin del sistema matricial

C1 11 1 + 51 !J 2 = VI

P2 ~1 + C2 2 + S" lJ. 3
~
= Y2

............. ~ .... ~ ........ ....... ......... .

~ ~
(13)

Pi t:" i-l + Ci fJ. i + Si IJ. i+l y.I

Si hacemos .. ( l; 2) en la primero ecuacin del sistema (13),

tendremos
45
 -1 e -1
III = - el SI .1 2 + 1

o sea
Al = MI !J. 2 + W

si MI = -C1- 1 SI y

remplazando en la segunda ecuacin de (13), obtenemos

o sea

de donde

si hacemos

M2 = (P2 MI + C2 ) -1 (-S2 )

obtendremos finalmente

En las dems ecuaciones se contina remplazando sucesivamente :

Vi' I = VI

tJ i = Mi Ll i+l + Wi con M = p Mi_1 + ei

W = (Pi Mi-l + en -1 (Vi - Pi Wi-l)

 Si tenemos por ejemplo, 19 etapas de clculo, lo ltimo ecuacin ser

con w19 = (P19 M 18 +

e 19) -1 (V19 -P19 W18)

En esto formo, se procede desde lJ 19 hacia atrs hasta LI 1 ' Y se obtie

ne un valar paro codo matriz cada vez que asignamos volares o los parme
tros de ajuste; tendremos, entonces, 01 final de codo uno de estos procesos;-
los niveles alcanzadas por el aguo en codo sub rea.

3.2.5 Esquema del programo de computador

El esquema general adjunto muestro lo secuencio de clculos o travs de los di-

ferentes etapas :

Plt'EL I
Lee los parmetros en tarjetas y los grabo en el disco;
se hoce un programo paro codo cree ido.

Lee los dotas de codo crecido (alturas de aguo en co-

da sub rea) y algunos parmetros que varran con h. Si
se tienen dotas cado 12 h., se descomponen en n, I:J. t
probando el modelo con los datos medidos:

t total de lo crecido !J t = 12
N = y
12 horas n Ll t

Hoce las interpolaciones necesarios paro calculas de

1/2 en 1/2 hora, los volares de los parmetros.

114.." ----.
I

S
______ J
Modifico y grabo nuevamente alguno$ parmetros que
yo estn grabados y deben ser cambiados se requiere
lo llave porque no siempre es necesario esto operacin.

Calculo, paro coda seccin de paso, todas los motri-

ces y encuentro los valores de h en el instante
t + t los entrego o travs de SALID.
'--~".~.~.6:::_:-:-~- - 1
IJlI""':, . . . ".. 1
______ J
Grabo en el disco los valores de los olturas mximas
de ogua en coda caja para comparar varias crecidas.

Es el programo que pruebo el resultado del modelo.

47
3.3 MODELO DE LA UNIVERSIDAD DE TEXAS

3.3.1 Antecedentes

Aunque el modelo Stanford IV es capaz de reproducir el escurrimiento adecuado

mente, se pens en la Univel'$idad de Texas que algunos de los procesos constit
tivos del modelo eran susceptibles de mejorane paro hacerlos ms compatibles co;,
los proc:esos ffsicos reales y utilizar parmetros con significado fsico definido. En
el modelo Stanford IV hay tres parmetros cuyos valores deben ser suministradas
por el usuario, y son los valores del almacenamiento en las zonas superior e infe
rior del suelo y el ndice de infiltracin, parmetros que estn inter-relacionadOS,
lo cual agrega considerable dificultod a lo determinacin de sus valores.

Para mejorar el modelo se consider que 105 conceptos de almacenamiento de de-

presin, infiltracin, contenido de humedad del suelo y movimiento del agua den-
tro del suelo deberan ser modificados y adaptados ms adecuadamente o los pro-
cesos conocidos de lo ffsica del suelo. Con esto en mente se efectu un anlisis
intensivo del flujo en la zona no saturado del suelo, para desarrollar un programo
de simulacin modificado.

Entre las fuerzas que acton para producir el movimiento del agua dentro del sue
lo pueden tenene en cuenta las siguientes: (Remsan, 1962)

o - lo fuerzo gravitacional que siempre acto en direccin vertical.

b - las fuerzas capilares
c - las fuerzas qumicas y osmticas
d - Fuerzas debidas a gradientres de presin po rcia I que acton sobre el vapor
de agua.

las fuerzas qumicos y osmticas generalmente originan movimientos muy lentos,

por lo cual no se tienen en cuenta en el anlisis. El volumen de agua que se
mueve en estado de vapor es iguolmente reducido comparado con el volumen l-
quido, por lo cual el estudio se limita a la consideracin de las dos primeras
fuanas mencionadas.

Se dedic mucho tiempo 01 estudio del modelo de diferencias finitas de una ca

lumna vertical de suelo, besado en el usa de una forma de la ecuacin de Oo';y
para flujo en la zona no saturada:

o(lf+y)
v = -K (1)
() y

y la ecuacin de continuidad
48
 = (2)
" t

donde
V es la velocidad en la direccin y, es decir, la rata de flujo
a travs de la unidad de rea;

K es el coeficiente de permeabilidad del suelo en unidades L

T
'f es el potencial copilar o tensin copilar en unidades de longi-
tud;

s es el contenido de humedad del suelo, y

y es la distancia vertical desde un plano arbitrario de referencia,

considerado positivo hocio arriba.

La formo de lo ecuacin resultante era oplicoble o lo zona no saturada pe-

ro no o lo zona saturado; sto condujo o lo necesidad de determinar lo loco
lizocin de la interfase entre los dos zonas como porte de lo solucin, con-
una forma apropiada de la ley de Darcy poro cada zona. Sin embargo, no
se encontraron mtodos satisfactorios para manejar este problema y evaluar
las derivados por diferencias finitos en la regin de la interfase, par locual
fue abandonada esto I(neo de accin.

Estas dificultades podran evitarle utilizando un mtodo que no use formas di

ferentes de la ley de Doreyen cado zona. Rubin (1966) ha hecho posible-
sto con una combinacin adecuada de permeabilidad, humedad del suelo y
tensin capilar; con esto tcnico es posible desarrollar un modelo de infiltra
cin ms adecuado que el descrito en los pginas siguientes, lo cual puede-
ser objeto de futuras investigaciones.

3.3.2 Caractersticos del modelo

A. Origen

Poro desarrollar el modelo en lo Universidad de Texas se hizo uno utiliza-

cin extensivo de los procesos del modelo Stonford IV, particularmente en
lo relacionado COn la escorrento superficial y lo propagacin. El estudio
del flujo en lo zona no saturada, deserHo atrs, se us como guo para el
proceso de clculo del mencionado flujo. En las figuras 10 Y 11 se muestran
esquemas del modelo Stanford IV y del modelo descrito aqu, cada una de
cuyos fases de simulacin ser objeto de un comentario.
49
 PRECIPITACION

EVA"ORACION

AGUA SlPERFICIAl.

Al.MACENA-
NTO SUPER
FICIAL

D
ALMACENAMlE/j Al.MACENA
TO SUBTE- TO EN l.A
RRANEO ZONA INFERIO

tHHt
"OTENCIAL
\ CORRIENTE

DIAGRAMA ESOUEMATlCO DEL MODELO STANFORO IV

FIGURA 10

50
 AREA AREA
PERMEABLE IMPERMEABLE

ESCOR~NTIA
ALMACENAMIENTO
DE INTERCEPCION

'k\~
~ ~..
",Vd/;

ALMACE-! ALMACE-
NAMIENTO NAMIENTO
DE DEPRE- DE FLUJO 1--. ESCORRENTlA
SIDN SUPERFICIAL

V
I
ALMACENAMIENTO
EN LA CAPA
SUPERIOR

L ~L....-----'
PARTE SUPERIOR
DE LA CAPA
INTERMEDIA
.------,
'---------"'ALMACENAMIENTO

SUBSLf'ERFICI AL

PARTE INFERIOR
DE LA CAPA
INTERMEDIA

ALMACENAMIENTO
1--. ESCORRENTI A
SUBTERRANEO

~
FLUJO BASE

DIAGRAMA ESQUEMATICO
DEL MODELO DE LA UNIVERSIDAD DE TEXAS

FIGURA 11

51
--~- _.
 El programo de computador (Clabom y Moore, 1970) fue escrito para aco-
modarlo o lluvia medido en intervalos irregulares o an en formo acumula-
da como en ciertos tipos de pluvimetros.

Estos datos se convierten 01 perrodo de tiempo especificado, que puede va-

riar desde un minuto hasta un dfa. Puede prepararse un registro similar con
los datos de descargas cuando los hay disponibles.

B. Almacenamiento de intercepcin

Todo la precipitacin que cae en el rea permeable se .prapaga para simu-

-lar el almacenamiento de intercepcin; el volumen de agua interceptado es
t controlado por un parmetro (VINSM) que vara mensualmente. Este pa-=
rmetro puede estimarse como una funcin del tipo y densidad del material
constitutivo de la capa superficial de suelo; en terrenos agrcolas, el alma-
cenamiento de intercepcin mximo puede ser pequeo cuando se hoce lo -
siembra, y crece cuando se desarrollan las plantas; por otra parte, en reas
cubiertas can bosques perennes o en zonas urbanizadas, es sensiblemente -
constante durante el ao.

c. Almacenamiento superficial

El aguo que excede el mxima almacenamiento de intercepcin llega hasta

la superficie y puede caer en uno de das tipos de almacenamiento: el de -
depresin y el de flujo superficial; el almacenamiento de depresin, es el
volumen retenido en las depresiones de todo tamaa que tengo la cuenco.

Dentro del rea permeable de la cuenco se pueden delimitar tres tipos de zo

nos, como se muestra en lo figura 12. Estas zonas estn caracterizados ca=-
mo (A), reas donde el flujo superficial alimenta directamente la corriente;
(B), reas donde el flujo superficial va o uno depresin; y (el, reas de 01
macenamiento de depresin; stos ltimos varan desde surcos causados por-
pisodas hasta estanques y otras estructuras que retienen el flujo.

El rea superficial del aguo almacenado en depresiones se considero como

una funcin lineal del volumen del almacenamiento (figura 13). Lo lluvia que
alcanza la superficie del suelo se divide en almacenamiento de flujo superfi
cial (LS) y almacenamiento de depresin (DS) en base 01 rea superficial de-
codo tipo de almacenamiento. L S se divide posteriormente en una parte
que va directamente o la corriente y otra que va a los depresiones; la in-
filtracin y lo evaporacin, tambin fenmenos superficiales, se basan en el
rea superficial de cada tipo de almacenamiento.

52
 TRES TIPOS DE AREAS EN LA CUENCA
Al Fluj-o superficial que d.,corvo
directam.nte a la corriente
8) Flujo superficial que vo o un
almacenamiento de deprllin
e) Almacenamienlo de dl"'"n

FIGURA 12

Volumen del almacenamiento de depresin

Volumen mx.imo del almacenamiento de depr.,on

RELACION ASUMIDA ENTRE AREA Y VOLUMEN

DEL ALMACENAMIENTO DE OEPRESION

FIGU RA 13

tl.,~ ..........."...... 'o H,.,.lot'~ Wail' Rn.e,ea 'ubllc.tl-., 'orl Col eo ........ '171

53
 El modelo de flujo superficial es tomado del modelo Stanford IV, constituida
por la ecuacin de IZZARD (formo para flujo cuasiturbulenta) :

OLF ~
64200 s
nL
( ~) 5/3 (1+0.6 (-&) 3\ 5/3 (3)

donde :

OLF es la rata de escorrenta en pulQ/hora por unidad de rea.

n es la rugosidad de Manning para flujo superficial (ROUGH).

L es la longitud del trecho donde hay escurrimiento superficial, en

pies (VLENGH).

s es la pendiente del trecho anterior (SLOPE)

D es el almacenamiento en ft 3/ft, definido como el promedio al pr~

cipio y 01 final del perodo (OFSAVG).

De es el almacenamiento en ft 3/ft para uno condicin de equilibrio

(OFSEQU), definido por:

~
0.000818 .0.6I n
0.6
OFSEQU (4)

donde i es la intensidad de precipitacin, en pulgadas por hora. Este

modelo se basa en las siguientes relaciones empricas entre flujo, profun-
didad y almacenamiento:

Y ~
D
T [ 1+0.6 ~) 3] (5)

donde Y es lo profundidod,en pies, en el extremo inferior del plano del

flujo.

Estas ecuaciones fueron desarrollodos DOr Lnsley y Crawford (1966) y don

buenos ajustes con los resultados obtenidos con la solucin de la ecuacin
diferencial parcial del flujo variado por diferencias finitas.

El flujo superficia!' generado en la parte de la cuenca clasificada como

contribucin directo a la corriente, se acumula en la variable RO; la par-
te originada en las reas que constituyen el almacenamiento de depresin
se suma a ste. Las depresiones pequei'las se llenan rpidamente 01 princi-
pio de lo lluvia y comienzan o producir escurrimiento, lo que no sucede
54

I
 con depresiones grandes, donde el escurrimiento comienza mucho tiempo des
pus. En la figura 15 se muestra una relacin tfpico entre el rea que pro
duce escurrimienta a travs del almacenamiento de depresin y el volumen-
de ese almacenamiento; la curva se asume como una parbola con vrtice
en el origen, y su ecuacin puede escribirse:

= (.2.
'c
Ds
DsMAX
) 0.5

donde Xa es la fraccin del rea originalmente designada como depresin

o su rea contribuyente pero que est ya produciendo escurrimiento.

DS es el volumen del almacenamiento de depresin, en pulgadas de pro-

fundidad.

DsMAX es el mximo OS, y

C es una constante, norma Imente con un va lar de 1.

La forma de esta relacin varra de una cuenca a otra; sin embargo, parece
que la parbola es adecuada para la mayorfa de las cuencas, segn inves-
tigacianes efectuadas por Clabom y Moore (1970)

D. Infiltracin

La infiltracin depende de: a) el agua disponible, es decir, la intensi-

dad de precipitacin, b) la capacidad del suelo para conducir el agua a
travs de la zona no saturada, y c) profundidad de la lmina de aguo que
cubre el suelo, el espesor de la zona satura y su permeabilidad.

La ecuacin (1) ?uede escribirse como

vo =
(7)

donde Va es la tasa de infiltrocin; esto ecuacin es apl icable en las con

diciones a) y b) dadas atrs, y s/<)y debe considerarse como la va-
riable dependiente.

En relacin con la figura 15, la tasa de infiltracin de la lluvia encharca

da se puede calcular as:

55
--1---

r PONO - 'f OC

l SATOPH

donde PON O y SATOPH tienen el significado mostrado en la figura, y

es el volar de la tensin en la interfase suelo-oire.
oc

Holtan (1967) propuso una ecuacin de infiltracin

n
f = a F + f (9)
p c

donde f es la rata de infiltracin cuando el abastecimiento no es el fac-

tor Iimitante.

fc es la rata de infiltracin constante, o permeabilidad saturada del suelo.

Fp es el volumen total que puede inA Itrarse antes de alcanzar una [ata de
infiltracin constante,

a y n son constantes.

Holtan recomienda tomar para a el valor 0.62 y para n 1.387. Fp es

una medida de los vacios que quedan en la columna de suelo en un tiem-
po dado, por lo cual es tambin una medida del agua en almacenamiento
dentro de la columna. Como la ecuacin de Holtan se aplica cuando exis-
te exceso de lluvia, entonces

= - f (10)
c

donde el signo menos es una convencin; si hacemos

K
sat
(PONO - If-' oc
SATOPH
1= o ( 11)

la ecuacin de Holtan se aproximo ecuacin (8). Holtan encontr a

trovs de datos experimentales que puede expresarse como

= KSo (12)

donde el subscrito cero indico evaluacin en el tiempo cero,

56
 - I AOI.FSF

o Almacenamiento de depresin
Mximo alrnocblmiento de depr.aiQn

REL.ACION ASUMIDA ENTRE VOL.UMEN DEL. AL.MACENAMIENTO

DE OEPRESION Y ESCURRIMIENTO. SUPERFICIAL. A PARTIR DE ESE AL.MACENAMlENTO
FIGURA '4

PONO
,
SATOPH '"
~

----------- '"
~
..e
a. Perfil d.1 suelo

o Humedad del suelo S

b. Perfil de humedad

INFIL.TRACION DE LA L.L.UVIA RETENIDA EN DEPRESIONES

FIGURA 15

K,
S Zona de f~lrocin
,/

1- L
-1

MODEL.O IDEAL.IZADO DE FL.UJO SUBSUPERFICIAL.

FIGURA 16

57
 K es el factor vegetativo dado en la tabla 1, Y
S es el espacio de poros disponible hasta una profundidad de 21 pulgadas.

TABLA

FACTOR VEGETATIVO DE HOLTAN

Cubierta Valor del factor K

,
gramaneas 1.00
hierba y olfalfa 0.70
forraies 0.45
alfalfa 0.35
maleza 0.30

La ecuacin de Holtan se aplica particularmente para modelos del tipo de

balance contrnuo en una cuenca, puesto que la variable independiente es
el espacio poroso no saturado en la capo superior de suelo. Definida esta
zona como la zona superior, la ecuacin (9) se puede escribir

f = a (VZST - UZS) n + SATPRM

(13)
= e 1 (VEST UZS)C2 + SATPRM

donde

UZST es el espacia total de poros en la zona superior, dada en pulg 3/

pulg 2 (porosidad por espesor)

UZS es el volumen normal de agua en lo zona superior, en pulg 3pulg2

SATPRM es la permeabilidad saturada de la zona superior, en pulgadas por

hora.

C1 y C2 son constantes de entrada, correspondientes a las de Holtan

01 K Y n. Holtan no hizo ninguno concesin respecto al aguo que aban-
dona la zona superior, es decir, el volumen de espacio poroso remanente

58
 en un instante t despus 'de comenzado e I aguacero, se ca Icu I como el
espacio poroso disponible de to al instante t. El modelo de la cuenca re
flejar ms exactamente el volumen real de aguo almacenada dentro de 1';-
zona superior. En consecuencia, esta zona puede tener menor espesor que
las 21 pulgadas establecidas por Holtan. En tierras agrcolas la profundidad
de cultivo alcanza siempre un valor razonable, mientras que en reas no
cultivados el espesor puede variar de medio pulgada, paro suelos muy com-
pactos, hasta varias pulgadas en suelos arenosos. La presenc ia de una zo-
na menos permeable interrumpe la zona superior. La ecuacin (13) d lo
tasa de infiltracin poro el perodo en que el abastecimiento no es el fac-
tor limitante, es decir, es lo infiltracin potencial. Cuando lo intensidad
de precipitacin es menar que la infiltracin potencial, toda la lluvia se
infiltra; sto disminuye el espacio porosa disponible, y en consecuencia dis
minuye la infiltracin potencial para el prximo perodo de tiempo. Rubin-
(1966) estableci que, para un abastecimiento uniforme menor que lo tasa
de infiltracin potencial, la columna de suelo eventualmente alcanza un es
todo de equilibrio en sus reas superiores donde el gradiente -

l K lJI +y)
oy
es justamente el adecuado para transferir el abastecimiento. Esto puede -
tambin ocurrir en el modelo, cuando el flujo que entra en la zona supe-
rior es exactanente igual 01 que sale. La ecuacin de Holtan, por cuanto
no es exacta, no ofrece una bueno aproximacin del proceso de infiltra-
cin, tanto en lo situacin de la infiltracin limitada por el abastecimien-
to, como en la de la infiltracin limitada por las condiciones del suelo.

E. Flujo o travs de la zona no saturada

Idealmente, se podra describir el flujo o travs de la zona no saturado,

con alguno formo de lo ecuacin (1). Sin embargo, en un modelo de cuen
ca sto debe hacerse en formo aproximado. En los modelos de balance coo
trnuo, se debe obtener el volumen de aguo paro codo intervalo de profun=-
didad; no obstante, cuando el nmero de estos intervalos aumenta, el tiem-
po necesario para real izar los clculos tambin aumenta y rpidamente se
encuentra un lmite econmico; por eso debe limitarse el nmero de tales
zonas. Las 'variables independientes en la ecuacin (1) se pueden e'scribir
como funciones de lo humedad del suelo, S.

ji = f :S' (14)

K = 9 :S) ( 15)

59
 la solucin de la ecuacin (1), debe dar

v = h (5) (16)

Se han propuesto varias formas para las funciones f (5) Y Y (5), Y han sido
discutidas por Clabom y Moore (1970). El conjunto de funciones escogidas
poro su usa en el programa de simulacin es

K = Al (5+8 ) C1 + DI (17)
1
_r = A2 (S+B2) C2 + D2 (18)

En la tabla 2 se muestran los valores de los coeficientes poro diferentes

suelos ensayados, y se puede usar como gura pora determinar esos coeficien-
tes.

El programa transfiere la humedad de una zona a la otro usando la ecua-

cin (1) escrita en la forma

(19)

El valor de la permeabilidad es el promedio ponderado de la permeabilidad

en cada zona como se determin en la ecuacin (17). El valor de '3of /ay
se obtiene a parti r de

= (20)

donde los subscritos 1 y 2 hacen referencia a las dos zonas. El valor de

11' se encuentra con la ecuacin (18); el valor de S poro cada zona se en
cuentra dividiendo la humedad total en la zona por su espesor.

F. Flujo subsuperficial (o hipodrmico)

El fenmeno del flujo subsuperficial, aunque observado frecuentemente en la

Naturaleza, no ha sido estudiado en detalle; en la figura 16 se muestra un
modelo simplificado del flujo subsuperficiol; la aparicin de agua en la zo-
na de filtracin depende de la pendiente de la zona menos permeable, S,
lo permeabilidad de esta zona en relacin con la que est encima, la pro-
fundidad de la zona menos permeable, y el rea contribuyente al proceso de
flujo subsuperficial, es decir, la distancia L en la figura 16. Estas variables
y tal vez otras se combinan poro convertir una porcin del hidrograma de in-
60
 TABLA 2 COEFICIENTES PARA APROXIMACIONES DE LA PERMEABILIDAD Y EL POTENCIAL CAPILAR

P E R M E A B I L I D A D

DESCRIPCION y FUENTE Promedio

A B C D
mox. R2 ,
R. E. R. E.
I
Limo
Banks y Bawers (l962) 797,7 -O, 114 7,75 2,9. 10- 9 0,400 0,139 0,949

Arcilla
Phillip (1957)
Gardner (1958\ 3,37 -0,103 6,01 -3, O 10- 7 0,372 I O, 147 0,936
Watson !l967) 206,8 0,002 4,06 -6,9.10-4 0,231 , 0,091 0,919
I

Greda !

Hanks y Bowers (1962) 1677,2 0,044 8,70 -1,6'10- 6 0,652 0,284 0,936

Arena seca
Liokopoulos (1966) 1,689 0,026 2,45 2,4'10- 3 cO,233 0,079 0,899
,

Arena hmeda
Liakopoulos 11966 1,983 -0,00585 2,66 3 2'10- 5 0,312 0,188 0,915
Vachaud (1966) 117,32 -0,1185 5,69 7', l' 10- 4 0,187 0,068 0,954
Remson (1965) 1,6104 0,087 25,16 I , 0'10- 9 .~!790 . 1, 179 O_l 883
- ------- --- -- - --~---
.. _ . _
p O T E N C I A L C A P I l A R
c -----
Limo 0,061 rO,063 r_7,, r
-8,11 1 0,014 10,004 1 0,963
Arcilla Da t os d a dos en fa r ma de ecuaClon .-
Watsan -47,6 -0,047 -0,033 33,74 0,008 0,001 0,923
Greda 26,20 0,7745 - 23, 91 0,965 2,68 0,753 0,898
Arena seca 225,3 1, O - 11, 35 3,64 0,59 0,239 0,900
Vachaud 9,937 0,6555 - 21,28 1,207 1,765 0,478 0,894
o- Remson _. - 249,76 0,816 - 30, 89 -2,653 0,480 0,167 0,909
- -----_.- ._~. - -- -- - .. --_._-

,
~,
 filtracin en el hidrograma de flujo subsuperficial. La formulacin matem-
tiCG ms simple del procesa serra entonces una que deteJminara el valumen
del fluio subsuperficial y produjera el hidragrama subsuperficial propaganda
este volumen; ste es el mtodo usado en este modelo. El volumen del flu-
jo subsuperficlal del nivel superior hacia la zona intermedia se determina con

YINFlO = C 16 (TlZS - CID VINST) b.t (21 )

donde
YINFLO es el volumen de agua adicionada al procesa subsuperf..!,
col durante el perodo lit.

TIZS es el volumen de agua almacenada en la mitad superior

de la zona intermedia.

VINST es el volumen total de la zona intermedia .

san parmetros de entrada que deben reflejar las diferen-

cias en la permeabil idad de las dos zonas.

Cuando la cantidad (TIZS - C 10 . VINST) es menor o igual a cero,

VINFLO se hace cero. C 10 parece ser un parmetro que establece el l-
mite inferior del volumen de agua que debe estar presente para que el flu-
io subsuperficial pueda iniciarse. La forma del hidrograma del fluio hipo-
drmico est determinada por una funcin de retraso.

(22)

donde
X es el volumen de aguo descargada a la corriente desde el
almacenamiento hipodrmico en el tiempo A. t

DEUNG es el volumen de agua en almacenamiento subsuperfic!al.

C14, C 15 son parmetros de entrada.

La ecuacin (22) no tiene explicacin fsica; es probable que C 15 pueda

tomarse como cero; por analoga con las ecuaciones de propagacin en cau
cas, la salida deber relacionarse al valumen de almacenamiento. -

G. Flujo subterrneo
Mucho de lo que se dijo sobre las dificultades de la modelacin matemtica
62
del proceso subsuperficial se aplica al flujo del agua subterrnea. El mode-
lo considera la posibilidad de flujo subterrneo directamente sobre la carrien
te y tambin fluia hacia afuera de la cuenca. Este ltimo Se detennina co=-
mo una funcin del agua almacenada en la porcin saturada del perfil del su~
lo.

23)

donde
y es el volumen de agua que sale de la cuenca en el pero-
do l> t

GWS es el volumen de agua en almacenamiento subterrneo, y

C 13 es un parmetro de entrada.

El flujo que va directamente a la corriente se determina de manera similar.

x = (GWS - C 11 ) . (24)

donde
X es el volumen de agua que fluye del sistema subterrneo a
la corriente durante t.

es el volumen de agua en almacenamiento subterrneo par

debajo del cual na hay flujo hacia la corriente.

C 12 es un parmetro de entrado.

Para el perado de flujo base del hidrogramo, la ecuacin !24) es anloga

a la ecuacin comunmente usada, descompuesta

G.1 = k. G.1- 1 = k i r-
V
o (25)

donde
G.1 es el flujo para el perodo de tiempo j- smo.

G.1- 1 es el flujo poro el perodo de tiempo Ii- n

Go es el flu jo al comienzo del primer perodo de tiempo, y

k es la constan te de recesin.

63
 Por la ecuac in (24)

K =
GWS + C11 (26)

GWS - Cll
Se ha notado que K no es constante, pero se acerca a cero cuando el al-
macenamiento subterrneo se aproxima a C 11' El efecto de variar la con s
tante de recesin es incrementar o reducir la recesin dada por la ecuaci
(25); sto se acenta con la remocin de un vol umen de agua subterrnea
que sale de la cuenca.

H. Evaporacin y evapotranspiracin

El proc~ de evaporacin depende en grado sumo de la energra reclbida en

forma de radiacin solar; la cantidad de energra recibida afecta tanto la ew
poracin total como su distribucin temporal. Otros factores climticos impor
tantes en el proceso de evaporacin son los vientos y la humedad re lativa; eT
efecto combinado de estos factores se reflejar en la cantidad de aguo que
se evapora desde un tanque. En consecuencia, el modelo hace uso de la -
evaporacin mensual de tanque para establecer la evaporacin potencial men
sual de la cuenca; de aquf se extrae un media diaria, y finalmente una tasa
instantnea. Se asume que lo evaporacin potencial comienza media hora -
despus de la salida del sol y va hasta una hora antes de la puesta, con su
mximo a los tres cuartos de ese tiempo. Para estimar la tasa instantnea de
evaporacin potencial se usa una curva en forma de campana asimtrica a;us
tada a estos tres puntos. La distribucin temporal de la evaporacin poten--
clal concuerda bien con los datos suministrados por Bavel (1966), como se -
muestra en la figura 17.

La evaporacin tiene lugar primero en el almacenamiento de intercepcin.

Si no hay suficiente agua almacenada allr para satisfacer la evaporacin
potencial de este perrada, sern retiradas cantidades adicionales de agua,
en orden, del alll\Qcenamiento de flujo superficial, almacenamiento de de-
presin y del almacenamiento en la zona superior del suelo. Se asume que
la evaporacin se real iza a ritmo potencial desde el almacenamiento super-
ficial y el de depresin, asf como en la zona superior, en la cual disminu
ye a medida que la zona decrece, por la dificultad del suelo para mantener
la humedad de la superficie. Esta situacin se aproxima en el modelo, usan
do -
UZS UZMIN) OC
E = EVAPOT - ( -:U~Z:-:::S=-T--~U~Z-:-M~IN'-'-- (27)

64
1.0 Superficie de D9UQ

E vaporacin
medido

Aproximocion /
.5
de una curva
sinusoidal '1
/
:
0600 1200 1800 1400

Hora del dio

Suelo hmldo
1.0

EvaporoctOn
/medidO

0600 1200 1800 2400

Horo del dia

APROXIMACIOH DE UNA CURVA SINUSOIDAL

PARA MEDIDAS DE EVAPORAC/ON (VAN BAVEL, 1966)

FIGURA 11

65
 donde
E es lo toso de evaporacin real de lo superficie del suelo.

EVAPOT es la evoporaciefl patencial

UES es el volumen corriente de agua almacenoda en la zona supe-

rior.

UZMIN es el volumen mnimo permitido en la zona superior.

UZST es el volumen total permitido en lo ZOna superior, y

oc 1',un parmetro que describe lo capacidad del suelo paro trans-

portar aguo hacia lo superficie.

lo transpiracin del aguo del suelo paro los necesidades de los planlcs pa-
rece ser independiente del contenido de humedad del suelo, cuando est por
encima del punto de marchitez; en consecuencia, en los meses donde hay
transpiracin, uno cantidad adicional de aguo va o ser extrado de lo zo-
na radicular (RTZONE); lo extracciefl est limitado por el punto de maTchi
tez, asumido igual al volumen mnimo para cado zona. El espesor de la:ro
na radicular no necesariamente coincide con lo zona superior, pues para-
rboles de races profundas es mucho mayor. As, la evaporacin mensual
medida en el tanque se aumenta con el uso consuntivo mensual (TRANPO).

l. Propagac in de descargas

El escurrimiento para cada perodo de tiempo, calculado por el modelo en

la forma descrita someramente, debe propagarse o travs de lo cuenca pa-
ro producir el hidrograma de salida; paro sto, se uso la tcnico del grfi-
co de distribucin empleado en el modelo Stanford IV. El flujo generado du
rante un perodo se distribuye de acuerdo con los parmetros de entrada e
intervalos futuros de escurrimiento, "S decir, el flujo que llega o lo corrien
te durante el intervalo actual de tiempo, aparecer en lo seccin de con=-
trol distribudo en el tiempo.

El programo tiene capacidad para simular cuencas compuestos de varias sub

cuencas o segmentos, en sentido aguas arribo -aguas abajo, por lo cual 10-
cinto magntica que contiene los dotas de precipitacin debe prepararse en
ese sentidc, tambin.

3.3.3 Aplicaciones del modelo

Aunque los oportunidades pora apl icor el programo han sido Iimitados, se ha -
66
usado en dos cuencas de caractersticos bien diferentes y con distintas objetivos.

Lo primera aplicacin se hizo a una cuenco pequeo (0.1 km2 cerca de

Edwardsville, IHinois, U. S.A., para la cual se dispona de informacin suficien-
te para la estimacin de los parmetros. Se persegua un propsito muy especi'fi
ca: seleccionar los parmetros doe la simulacin con anticipacin, can base en las
caractersticas conocidas de la cu"nca. Se trota en realidad de una prueba muy
severa del programa y de la habilidad para interpretar las caractersticas conoci
das de lo cuenca en relacin con los procesas fsicos que se estn simulando. -

Inicialmente se escogieron algunos hidrogramas tpicos de los diferentes condicio

nes encontradas en el modelo; en lo figuro 18 soe muestro uno comparacin entre-
los hidrogramos observado y calculado en uno situacin tpica.

Aunque los resultados podran mejorarse, son aceptables dentro de los condicio-
nes en que fu usado el modelo; anlisis posteriores mostraron que los descargos
en ste perodo fueron afectadas por fenmenos de deshielo, que no estn con-
templados en e I modelo.

La segunda aplicacin se hizo en el arroyo Mukewater en Texas, U.S.A. (183

km 2), 01 cual haba sido aplicado anteriormente el modelo Stonford IV, del cual
se tomaron los parmetros comunes. Los resultados de la simulacin para agua-
ceros individuales se muestron en las figuras 19 y 20 Y las descargas medias dio
rias en la figuro 21.

Se puede observar que lo correspondencia entre descargas observadas y simula-

das puede mejorarse si se tiene en cuenta la distribucin espacial de las Iluvia,.
hecho que, sin embargo, no ha sido contemplado ni en ste modelo, ni en el mo
delo Stonford IV.

67
 .,..------------,
ti

...

.
'.


o
o
.--y-r "10 11 O .

COMPARACION DE DESCARGAS MEDIDA

Y SIMULADA PARA EL AGUACERO DE o
JUNIO 26, 1942
FIGURA 18 ..
o

1000

l.

o "--- AGUACERO DE MAYO 17 - 18, 1957

FIGURA 20

"
AGUACERO DE' MAYO 11-13, 1957

'.eUItA l '

68
 1.6 r
~ 14 r
o
'O
o
'" 1.2
"5 r-
"-
1.0
e
-<>
0.8
~

"
u 0.6
;-
r -
"- 04 -
0.2
r
10
,
20
r-l
30
r

" , , .,.,
- -
- ----
Medido
Simulada
1000

1400
~

"
$
1200

."..,o ,....
~2 1000 e

..
.2
'O
BOO
-
E
o
:s'"u BOO

O r-,
400 '-
-..,
'- rt-, ,

--
200 I

h Q "A. , 6..
20
,-
30
10

PRECIPITACION y ESCURRIMIENTO PARA MAYO 1957

FIGURA 21

... "....... c. "5.,.......... rou. . . H,drOIOey'; Wa.tef R g ..... c . . P"bllcCllion" Forl COllml ColOl'odo, 1971
69
3.4 MODELO S. S.A. R. R.
(Streamflaw Synthesis and Reservair Regulation)

3.4.1 Generalidades
Se trata de un modelo pluviohidrolgico con propagacin par embalses que per-
mite transformar precipitacin en descarga, trabajando con tres componentes de
la escorrentfa: flujo superficial, subsuperficial y base (ver figuro 22). El modelo
completo fu desarrollado poro sintetizar la respuesta de un complejo sistema flu
vial a la entrada de precipitacin y temperatura, adems del efecto de regula--
cin de estructuras construidas par el hombre.

3.4.2 Evolucin histrica del modelo

El modelo SSARR, de caracterfsticas muy generales y de gran flexibilidad, fu

concebido par David M. Rockwood* para sintetizar descargas y ser usado en pla
neamienta y previsin operacional diaria de sistemas hfdricos. La versin origi-:
nal del programa fu desarrollada para un computador IBM 650 en 1957, para sin
tetizar descargas en la cuenca del rfo Columbia, diseado especfficamente para-
las condiciones de esta cuenca. Inicialmente fu formulado para simular escurr
miento de origen nival, y fu usado operacianalmente para previsin de escurri-
miento en la mencionada cuenca durante el perfodo de aguas altas, asf comopa
ra determinar avenidas de diseo para proyectos de embalses en el rea. -

Con la experiencia asf adquirida, en 1960 se iniciaron trabajos para mejorar el

programa lgico, tanto en lo que se refiere a tcnicas hidrolgicas como a lo
utilizacin del computador. Esta "segunda generacin" diseada para un IBM
1920, se describi cama un modelo general para sfntesis de descargas que pue-
de ser aplicado a cualquier cuenca sin limitaciones en cuanto a configuracin,
rgimen hidrolgico a sistema de embalses. El programa se us extensivamente
paro previsin de descargas, regulacin de embalses, estudias de caudales de di-
seo y regulacin de inundaciones en EE. UU., Alaska y Canad; se aplic en
reas de rgimen combinado (pluvial-nival) y, en conexin con el programa de
entrenamiento en anlisis de sistemas del Comit del rfo Mekong, se utiliz en
los tributarios de este rfo.

El programa actual es la "tercero generacin", y es diseado como un modelo

hidrolgico y de regulacin de embalses, con aplicacin o un ranga mucho ma-
yor de condiciones que la versin anterior, y con notable perfeccionamiento pa
ro simular procesos hidrolgicos, tcnicas de propagacin en canales y mtodos-
de regulacin de embalses. Es completamente flexible en su capacidad para re-

*Jefe de la Seccin de Hidrologfa e Hidrulico, North Pacific Division, U.S. Corps

of Engenee rs.
70
 LLUVIA

Evaporacin

Flujo Total

Flujo Flujo Flujo

Bose Subsuperficial Superficial

Prop iOn PropaqaciCn

OB OSS Os

Trocho 2 I

ESQUEMA DEL MODELO

,
FIGURA 22

71
 presentar relaciones frsicas de varias variables utilizando tablas previamente es-
tructuradas; estas relaciones pueden ser lineales o n y derivadas, ya seo de me
diciones directas o de funciones generales, a travs de coeficientes empfricos;
para sintetizar las condiciones naturales del escurrimiento.

3.4.3 c.aracterrsticas del Programa SSARR

A.Conceptos hidrolgicos inclurdos

Una cuenco de drenaje se puede separar en unidades hidrolgicas homogneas
de un tamai'lo y carcter tales que pueden ser usadas como una representacin
lgica de una cuenca moyor y sus subreas elementales; el almacenamiento en
el cauce puede determinarse para trechos de rro y representar asr el retardo
natural al escurrimiento que se encuentra en los sistemas naturales; en igual
forma, el almacenamiento en lagos o embalses puede evaluarse de acuerdo con
las condiciones especificadas de capacidad de los mismos. En esta forma, el
escurrimiento puede ser determinado tanto para el cauce principal coma para
105 principales tributarios.

El pragrama elaborado para este madelo combina la evaluacin de varias fun-

ciones hidrometeorolgicas que representan el praceso completo de simulacin
del escurrimiento; puede ser tan general que el escurrimiento puede sintetizar-
se ajustondo las diferentes parmetros de una cuenca particular que afecten la
evaluacin del volumen de agua de entrada (lluvia o nieve), la humedad del
suelo, la evapatranspiracin, el almacenamiento de depresin, almacenamien-
to superficial, subsuperficial o subterrneo, la infiltracin o los acui'feros,tiem
po de retardo de cada tipo de escurrimiento y efectos del almacenamiento en-
el cauce; cuando estos parmetras se han evaluado para una cuenca dada, re-
presentan la calibracin de esa rea particular o unidad homognea.

B. Subdivisin de lo cuenca
El escurrimiento se puede sintetizar digitalmente a partir de una serie de sub-
divisiones de la cuenca, cada una de las cuales puede estor individualmente
representada por su conjunto particular de coeficientes y funciones, y no hoy
Ifmite para el nmero que de toles subdivisiones se quiera hocer. La descar
90 calculada de tadas estas pequeas subdivisiones representa la entrado 01-
sistema del canal fluvial donde es propagada y combinado para constituir el -
flujo total en cualquier punto de inters.

C. Intervalos de tiempo
Los ecuaciones diferenciales que definen el flujo se resuelven por diferencias
finitos; el incremento de tiempo poro propagacin es determinado en codo coso
particular (puede variar de tres horas a un dro) y se puede modificar en el
72
 transcurso de una misma simulacin _ Los volares de entrado de las variables
hidrometeorolgicas se darn segn el intervalo de tiempo en el cual se requie
ren los datos de sal ido y el programo sintetizar el escurrimiento pora cualquier
duracin total de tiempo, dependiendo del tipo de estudio. Por ejemplo, se pue
den sintetizar descargos paro estudiar sus variaciones en perodos de un ao 0-
ms en estudios hidrolgicos de cuencas; poro sintetizar descargas de diseo a
portir de lluvias o para previsin, pueden usarse perodos de unos pocos das -
hasta varios meses de durac in.

D. Datos de entrada
Los datos hidrometeorolgicos que usa el programa se suministran en un gran
nmero de formas poro permitir flexibilidod en el manejo de la gran cantidad
de datos que generalmente requiere la simulacin paramtrica. Estos datos son,
precipitacin diaria, temperatura media diaria del aire, a valores de variables
hidrometeoro lgicas que afectan la fusin de la nieve, cuando es el coso.

Para operar el programa, es necesario proveer las condiciones iniciales, que

representan las condiciones de todas las variables funciones del tiempo en el
instante t o , y que son suministrados en tarjetas especiales; en esta fanma, la
sntesis puede comenzarse en cualquier instante deseado sin tener en cuenta
las condiciones hidrolgicas. Adems, existe la posibilidad de ajustar los con
diciones iniciales de valores previamente calculados, a condiciones observadas
en un punto particular en el tiempo, propiedad muy adecuado en la operacin
del modelo poro previsin de descargas diarios.

E. Usos del Modelo

Como se dijo antes, el programa fu desarrollado para su utilizacin en cual
quier estudio que incluyo el cmputo de descargos en uno comple jo cuenca-
fluvial, o partir de lluvias (y fusin de nieve).

En ingeniera se puede usar paro clculo de avenidas de diseo, estudios de

regulacin de embolses, previsin de crecidas y detenminacin de normas de
operacin de embalses. Tambin puede usarse poro extender perodos cortos
de descargos en base a precipitaciones y en estudios poro detenminar relacio
nes precipitacin - escurrimiento y caracter(sticas del almacenamiento en 10-
cuenco o travs de procesos de tentativas.

3.4.4 Tipos bsicos de simulacin empleados en el modelo

En el esquema simplificado que se expone en el presente trabajo, se describen

los dos tipos bsicos de simulacin que se incluyen en el modelo; la transforma-
cin lluvia - descarga (denominado tambin propagacin en lo cuenco o propaga
cin lluvia - descarga), y la propagacin descargo - descargo o propagacin e
el couce.

73
 la primera incluye todos los mecanismos de conversin de Ifuvia en descarga a
travs de cada uno de los procesos que constituyen el ciclo hidrolgico y sinteti
za el hidrograma generado por uno lluvia (no se considerar el problema de fu=
sin de nieve) en un punto dado de la cuenca. La segunda toma como entrada
ese hidrograma y lo propaga hasta otro punto del cauce, haciendo uso de las ca-
racterrsticas hidrulicas y geomtricas del mismo (ver figura 23).

3.4.5 Transformacin lluvia - descarga

El modelo hidrolgico utilizado en el programa SSARR considera los siguientesele

mentas hidrolgicos bsicos que deben ser evaluados individualmente para cada sub
rea en que fu dividida la cuenca:

o - Precipitacin neta en la cuenca; se determina ponderando los datos observa

dos.

b - Incremento de humedad del suelo; se determino como un (ndice (variable

temporal) que representa las variaciones de la humedad del suelo.

c - Evapotranspiracin; se determina como evapotranspiracin potencial expre-

sada en fundon de los valores medios mensuales en la cuenca y ajus/adame
dionte una reduccin en los dfas de lluvia en funcin de su intensidad. -

d - Exceso de escurrimiento; se determina a partir de la precipitacin neta del

suelo, conformando asf la cantidad de agua que entra a cada una de los tres
zonas de almacenamiento temporal.

e - Almacenamiento de retencin; representa el retardo temparal causada por

las pequei'las irregularidades de la superficie del suelo que na tienen dre-
naje libre.

f - Almacenamiento de superficie; representa el efecto de almacenamiento de

las capas superiores del suelo, causante de retardo en la componente super-
ficial.

g - Almacenamiento subsuperficial ; representa el efecto de almacenamiento de

los acufferos intermedios, causantes de retardo en la componente subsuperfi-
cial.

h - Almacenamiento subterrneo; representa el efecto de'la percolacin profun-

da hacia los acurFeros inferiores y su retardo en el flua base.

i - Tiempos de propagacin; para calcular los efectos del retardo en cada una
de las componentes del flujo.

74
 Propa9Qcin en la cuenca

Hidroqramo SintetiZado
, __-}---- IL__::::==- en A

ARa CUbierla PropoQOCn dlscaroa destorga

l ~rQmQ
por D lIuyiQ

_godo
_ _-... hasla B

TIPOS DE SIMULACION

FI GURA 23

ROP

" .... "" " , \

"" Diferentes funciones

, , '"
/
posibles

,
~~~~ ________________________--'SMI

FUNCION PARAMETRO ROP - SMI

FIGURA 24

75
 j - Relaciones para separacin del flujo; permiten calcular la porcin del exc!
so de agua que entra en cadO zona de almacenamiento.

Can el uso de estos fndices y relaciones, el modelo es una adecuada repre-

sentacin terica de los procesos ffsicos del escurrimiento, que lo hace utili
zable en grandes segmentos de cuenca teniendo en cuenta la disponibilidad
de datos y otras limitaciones prcticas.

A. Determinacin de la precipitacin en la cuenca

La precipitacin en la cuenca se calcula a partir de los datos diarios de un
cierto nmero de pluvimetros, 05(:

Pl Wl + P2*W2 .... +P n * Wn )
Pd =
s
donde
p.J = precipitacin diaria en las n estaciones

W = factor de ponderacin (en %) de cada pluvimetro

S = nme ro tata I de estac iones

Pd = precipitacin diaria ponderada de la cuenca

El factor de ponderacin W depende del rea de influencia de cada pluvi

metro, de la precipitacin normal anual en esa parte de la subcuenca y efe
la precipitacin normal anual en cada estacin, incluyendo asr la tendencia
que pueda tener la precipitacin. El factor se calcula asr:

An (NAPa) * N
Wn =
(NAP n )

donde
Wn = factor de ponderacin para la estacin n

An = rea de influencia de la estacin n, en % del total

NAPa = precipitacin normal anual para el rea An

NAPn = precipitacin normal anual en la estacin n

N = nmero total de estaciones usadas para calcular la preci- '

pitacin en la sub-rea.
76
B. Distribucin de lo 11 uvia

Codo dfo se subdivide en perradas, no siempre iguales. Cuando no hoy plu

vigrofos que permitan hacer el cmputo exacto, pueden hacerse estimocioes
basados en lo experiencia. Lo precipitacin ponderado poro codo perrodo es
entonces
WP = Pd * DIST
donde
W = precipitacin ponderado poro un perrada dado
p
DIST = distribucin (en%) de lo lluvia diario poro ese periodo;
puede calcularse as;

DIST = x TOO

donde
precipitacin registrado por un pluvigrofo durante el pe-
rodo
precipitacin total en el da

C. Clculo de lo escorrento total

Lo porte del volumen total de agua de entrodo, disponible poro producir es

currimiento superficial, subsuperficiol y base, se encuentro o partir de relo--
ciones derivados empricamente entre el Porcentaje de Escurrimiento ROP
(Runoff Percent) y el Indice de Humedad del Suelo SMI (Soil Moisture Index).
Poro uno cuenco particular, sto relacin se entrego 01 computador en for-
mo de tabla. El flujo total generado por perodo (RGP) est dodo por

RGP = ROP * WP
donde
RGP = flujo total por perodo

WP = precipitacin ponderado por perrada

ROP = % que indico qu porte de lo lluvia WP se convirti en

escorrentro. Es uno funcin del ndice de humedad del sue
lo, 5MI; ROP es uno funcin parmetro que puede tener-
la formo indicado en la figuro 24; la formo de lo CUrvO
porte de lo aceptacin de gue el mnimo ROP se obtiene
poro el 5MI mnimo y consecuentemente poro el mximo.
El mximo volar de SMI es el que debe o;ustarse en el ma
dela.

77
 D. Clculo de la humedad del suelo
El ndice de humedad del suelo SMI es una medida de la humedad del suela
usada para determinar el escurrimiento. SMI se considera nula cuando la -
cantidad de 'agua contenida en e l suelo alcanzo el punto de marchitez per-
manente; cuando la lluvia recarga la humedad del suelo, el valor de SMI
aumenta hasta alcanzar un mximo valor considerado como la capacidad de
compo. El valor de SMI al final de cado perrada se calcula mediante la
ecuacin
PH * KE * ElI
SMI2 = SMII + (WP-RGP) -
24
donde
SMI
1 = ndice de humedad del suelo al comienzo del perrada
SMI = rndice de humedad del suelo al final del perado

PH = duracin de I perrodo en horas

ElI = ndice de evapotranspiracin diaria expresada en base a la

media mensual

KE = factor de reduccin de la evaporacin en d(as de lluvia,

especificado al computador en forma de tabla de KE contra
intensidad de precipitacin. En la figura 25 puede apreciar
se que KE es igual a 1 cuando no hay lluvia, par lo cuar
SMI se reducir durante el perrodo por el factor constante
(PH/24) * ElI; KE es, entonces, otra funcin parmetra, de
creciente con la intensidad de lluvia.

E. Determinacin del flujo base

El volumen de agua que contribuye a la formacin del flujo base es una
funcin del Indice de Infiltracin del Flujo Base, BII (base flow infiltration
index), valor sinttico representativo de la influencia ponderada del escurri
miento precedente; ste dispositivo simula la elevada infiltracin al princi:
pio de una lluvia, con una disminucin paulatina a medida que sta progre
so. La relacin se suministra al computador como una funcin contnua en-
forma de tabla; el rndice de infiltracin del flujo base se calcula para cada
perodo con la ecuacin

BIII + (24* RG - BIII ) *

PH
TSBII + PH
donde
2
BII = ndice de infiltracin del flujo base (en cm/24 h) al comien
zo de I perodo

78
 KE

Varios funciones
1.0
/ , posibles

0.5

QO~------------------------~------' mm/h

FUNCION PARAMETRO KE - i

FIGURA 25

BFP

Vorias funciones
-..................... ~ posibles

>-< "- ....

......
"-,,-
"- ....
....
......
BU

FUNCION PARAMETRO BFP - BIl

FIGURA 26

79
 = rndice de infiltracin del flujo base (en cm/24 h) 01 final
del perrodo
RG = RGP/PH = escurrimiento en cm/hara

TSB 11 = tiempo de retardo, Ts (tiempo de almacenamiento) pora el

clculo del cambio de BII.

Puede asimilarse BII o un rndice de almacenamiento de depresin que man

tiene el escurrimiento disponible poro percolacin profundo. El flujo base
(RB)se calculo entonces asr:

RB = RG * BFP

donde BFP es uno funcin del rndice de infiltracin, BII, como puede ver-
se en lo figuro 26; BFP es, entonces, el % del fuja total que se convierte
en flujo base. Por lo tonto lo entrado poro el flujo superficial y subsuperfi
cial RGS o flujo directo se calculo con: -

RGS = RG * (1. O - BFP) = RG-RB

F. Separacin de los flujos superficial y subsuperficial

Para uno cuenco dado, lo relacin del escurrimiento superficial RS con el
flujo directo RGS se entrego 01 computador mediante uno tabla en lo cual
se puede asumir cualquier tipo de relacin, pero en general lo separacin
se baso en los siguientes consideraciones:

1. El mrnimo volar de lo componente superficial RS es el 10% del total RGS.

2. Lo componente subsuperficial alcanzo un mximo KSS y permanece cons

tante pora volares de RGS mayores que el 200% de KSS.

3. Los volores usados comunmente se calculan as(:

RS = [ O. 1 + 0.2 (RGS/KSS) ] RGS
y :
si RS < KSS, RSS = RG - RS
si RS > KSS, RSS = KSS y RS = RGS - RSS
donde
RS = componente superficial

RGS = fl ujo di recto

KSS = mximo valor de lo componente subsuperficial

80
 En la figura 27 se aprecia una interpretacin geomtrica de estos ndices y
la figura 28 es un esquema representativo de lo separacin de los diferentes
volmenes de agua.

G. Propagacin de las componentes

Una vez calculodo el volumen de cada omponente, se convierte en m 3/$
-
aSI !
OB =
RBxAxl0 2
36
RSS x A x 10 2
055 = A en Km 2
36

os =
RS x A x 10 2
36

Hecho sto, se efecta la propagaci-, de cada componente o travs de una

serie de embalses ficticios como se explicar en la parte correspondiente o
la propagacin descargo - descargo; en esta forma, como se muestra en la
figura 29, se ca' vi ,-te un hidrograma rectangular, representativo del volu-
men sintetizado para cada fase, en hidrogromos muy semejantes a los natura
les; la modificacin de la forma de esos hidrogramas est determinada por el
nmero de incrementos considerados y el tiempo de almacenamiento por in-
cremento.

H. Definicin de los caractersticas de la cuenco

El uso del modelo depende de la determinacin de varios coeficientes y re

lociones especrficas para una cuenco dada; algunas de esas relaciones so
generales y oplicables a otras subcuencas dentro de la misma rea y otras
se determinan especiticamente para una sola cuenca.

El proceso bsica para definir esas caractersticas es el uso del programa de

computador en tentativas sucesivas hasta obtener el mejor ajuste a los datos
histricos; ste proceso puede considerarse anlogo al ajuste de un modelo
reducido,modificando la rugosidad del canal, por ejemplo.

Los rangos de variacin de algunas variables son conocidos con base en ex-
periencias y anlisis ya efectuados en los diferentes campos de lo hidrologa;
as, las tasas de evopotranspiracin y la capacidad de infiltracin de los sue
los estn dentro de ciertos lmites conocidos; las relaciones entre lo precipi=-

] tocin puntual (en un pluvimetrol y la precipitacin media en uno cuenca,

se conocen relativamente bien cuando se dispone de datos adecuados; sin
embargo, como en la mayora de las cuencas grandes no existen datos sufi-
cientes' en el modelo SSARR el trmino usodo es ms bien un ndice de lo
SI
 RS,RGS

w.==~-i~-------L--------~------ __________-'RGS
INOICES OE SEPARACION OE FLU.JOS

FIGURA 27

RGS
Lluwla RG P
RS

SEPARACION DE FL.U.JOS

FIGURA 28

82
."', ..LL L",."
 precipitacin media en la cuenca. En un sentido similar se usa lo separa-
cin del Flujo en tres componentes, ya que es un concepto difcil de definir
fsicamente con exactitud y que, o pesar de ser hecho a travs de funciones
arbitrarias, se puede demostrar que mejora notablemente el comportamiento
del modelo. Las tiempos de almacenamiento (Ts) son tambin forzosamente
empricos y se determinan, igualmente, por tentativas.

En resumen, los siguientes son las elementos y funciones que requieren ser
definidos para representar el escurrimiento en el modelo de uno cuenca:

a - Estaciones de precipitacin y foctores de panderaciPn de codo uno pa-

ra calcular la lluvia media.

b - Relocin del ndice de humedad d.. 1 suelo (5MI) con el porcentaje de

escurrimiento (ROP), entregada en forma de tabla.

c - Indice de evapotranspirocin (ETI), expresado como evapotronspiracin

potencial media mensual.

d - Fadar de reduccin de lo evaporacin en das de lluvia expresado en

forma de tabla, como uno funcin de la intensidad de precipitacin.

e - Relacin del ndice de infiltracin del flujo base (BII) con el parcenta-
je del flujo total que se convierte en flujo base (BFP), en formo de ta-
bla.

f - Tiempo de retardo TSB 1I expresado como un solo volar en horas.

g - Relaciones de separacin de flujos superficial y subsuperficial expresa-

dos en una tab la.

h - Coeficientes poro propagaci n de los componentes, expresados como:

o) nmero de incrementos (o embalses) y b) tiempo de almacenamiento
por incremento, en horas, paro cada componente.

i - Area de la cuenco en km 2 .

3.4.6 Propagacin en el cauce

A. Ecuaciones bsicas

Como todo modelo o embalses, utiliza la ecuacin de continuidad y uno

ley de descargo de los embolses:

83
 o . Q

HlcIrOgramo
inlciol

/"

TRANSFORMACION DE HIDROGRAMAS

FIGURA 29

Ta

Tli f--~

1
()j o

I
VARIACION DE Ta

FI!ilURA 30
Ecuacin de continuidad:

11 + '2 01 + =---
2 liS
2 2 AT
o sea:
1 + 2 l!.S
ds (1)
1m -
2 AT
o = "'"'dt

Donde It y 0t son lo entrada y sal ida respectivamente y ds/dl la rato

de cambio del almacenamiento en el tiempo

Paro casos donde el almacenamiento es slo funcin de lo salida como en lo

gas o trechos cortos de ro donde lo cuo de almacenamiento es desprecio::
ble en comparacin con el prisma de almacenamiento, lo ecuacin de desear
go de los embolses es (ver figuro 30) :

s = Ts * O S = Ts * a (2)

donde Ts es el factor de proporcionalidad entre el almacenamiento y lo des

carga de solido.

En el programo SSARR, la variacin de Ts con la descarga se calcula con

la ecuacin
KTs
Ts = - - -
on (3)

donde:
Ts = tiempo de almacenamiento poro un incremento dado, en ho-
ros
KTs = coefic ient,e que representa el tiempo de retardo de I trecho
en cuestin (constante)
a '" descarga al comienzo del perodo
n = coeficiente (O <: n <: 1 ) relativo o la variacin del tiem
po de almacenamiento en funcin de la descarga.

Diferenciando (3) tendremos:

nTs On-1 da + OndTs = O

divi d iendo por O n-1

nTs da + OdTs = O

85
 o sea dTs '" _ ( Ts ) n (4)
dO e Q

Por otra parte, derivando (2) tendremos :

dS dO dTs dO
= Ts * + O 15)
dt dt dO dt

Remplazando (4) en (5)

dS
'" Ts
dQ
+Q (-~) dQ
'" (l-n) Ts
dO (6)
dt dt O dt dt

S se toma entonces un valor de KTs que involucre la condicin (1-n) de

la ecuacin (6), se puede reformular la ecuacin como

dS dQ
--
dt
= Ts
dt
Y llevarla o la de continuidad (1) para dar :
01 + 02 ~O
1m - 2
= Ts
6t

1m At -6t-
2
1 - -r
lit 02
= Ts2 - Ts0 1

2 (Ts + lit ) = (Ts - ~) 01 + Im .6. t

-r 2

6t
Ts - --
2 1m
2 =
1 +
<1t
Ts +...M...... Ts + .M.-
2 2

como llt
Ts -
2 M
Ts +.M-
= 1 - Ts+~
i
2 2
tl.t
\ tendremos 2 = 1 + ( 1m - 1 ) (7)
Ts At
+
--z-
86
\,
\

1,j,i""', e' eh
 que es la ecuacin bsica para hallar la descarga de salida 02 de un tre-
cha cualquiera en un intervalo de tiempo t.t, conocida la de entrada 01

B. Interpretacin geomtrica

En la figura 31 se representa un esquema de los hidragromas de entrada y

salida y se puede visualizar fcilmente el significado de los trminos de la
ecuacin (7), la cual se usa en el programa SSARR en incrementos finitos de
tiempo a travs del nmero predeterminado de embalses poro obtener final-
mente el hidrograma de salida, como se muestra en la figura 32.

C. Esquema de clculo
El esquema de la figura 33 muestra e I proceso iterativo de clculo usado en
el programa de simulacin; conocidas las condiciones iniciales representadas
por la primera fila, se hace la propagacin a travs de un embalse durante
un incremento de tiempo; al tratrse de un rra, ste se descompondrra en n
trechos, cada uno de los cuales se asimila a un embalse; los pasos del pro-
ceso son los siguientes: (figura 34) '.

a - Clculo de :
1m =
11 + 12
2

b - Clculo de : 1m - 1
c - Con 01 se encuentra Ts en la tabla que define esta relacin.

~t
d -
60 = [
, Ts + tH 1( 1m - 1 )
J
2

e - 2 = 01 + t::.o

Las valores de 01 y 02 se convierten entonces en los 11 e 12 para el si-

guiente embalse, y as puede continuarse el clculo del cuadra por filas ca
mo se hace normalmente, poro producir un perfil instantneo del ro ya que
cada Irnea se refiere a un instante de tiempo determinado. El clculo se
puede hacer tambin por columnas para producir el hidragrama en una sec-
cin dada del ro.

D. Propagacin con efecto de remanso

En casos normales, la propagacin efectuada con el modelo SSARR se hace
en circunstancias donde no hay efectos apreciables de las condiciones aguas
87
 1,0

-"1IIo ..... OIliO di onlroda

/
1, l. Htdr_ di ootido
I
----
..................
,.,,--.!Q!. ... ....
dI
... ....
.... ....
dO,
dI

t,o
--- T,

f-111/. --l-

t,

1
lit
11' Im- 01 , r, .Al.
M 1

+

} - - - 111 -----i

INTERPAETACIOH HOMETRICA DE !.AS ECUACIONES

FIGURA 31

88
 Q

Hidro9rama
de salida

----
OBTENCION DEL HIDROGRAMA DE SAL.IDA

FIGURA 32

~ O 121 n 1 n

Condiciones
O Qo,o QI,O Qz,o Qn-I,o Qn,o ..-.-..- ni e i oles

1 I Q 0,1 Q 1,1 Q2,1

2 QO,2 Qt,2 j Qz,z I

I !
!

1ft Qo,m Qn,m

l ConcflCionts lmUes (Hidro9romo de entrado)

ESQUEMA DE CALCULO
FIGURA 33

89
 x
t o

o ro 00

Il 01

IN ICIACION DE LAS ITERACIONES

fiGURA 34

HIDROGRAMAS DE LA CUENCA DEL RIO SE BANG HIEHG,LAOS

FIGURA
, 35
90
 abajo, es decir donde hoy uno relacin biunvoca entre nivel y descarga en
cualquier lugar Q lo Jorgo del canal.

Hoy cosos, sin embargo, donde la descargo es afectada en diversos grados

por los niveles en diferentes lugares, como en los trechos aguas arriba de
embalses o de la desembocadura de afluentes importantes y en estuarios de
ros donde hoy efecto de marea.

En estos casos el modelo SSARR trobaja con una relacin de tres variables
entre nivel y descargo paro lo propagacin, que puede expresarse en las
formas:
1 = 'E 1 , E2
o = i ,El, 02
2
donde
1 = descarga en el lugar

E1 = nivel del aguo en el lugar

E2 = nivel del aguo en el lugar 2

= descarga en el lugar 2
2

Estas relaciones especficas pueden derivarse de clculos previos a

partir de datos observados que relacionen descargos a niveles en -
ciertos puntos especiales de observacin.

3.4.7 Aplicaciones prcticas del modelo SSARR en sntesis de des-

cargos

Para i1ustror la aplicacin del modelo completo a lo cuenca del ro Mekong, en

lo figuro 35 se muestro el grado de ajuste obtenido en una porte de lo simula-
cin del ro Bang Hieng en Loas, can una cuenca de 3.600 km 2 de rea, ajusta
do para un perodo de 4 aos y cuyos resultados han sido considerados aceptables
si se tiene en cuento lo gran deficiencia de informacin pluviomtrica en eso -
zona.

El mecanismo de propagacin descarga - descarga del modelo SSARR fue usado

igualmente can buen xito en un modelo de simulacin de descargas e<'I el ro
Bogot (rea 1.418 km 2) (ref.12) y fu tambin involucrado coma elemento en
un modelo desarrollado para sintetizar descargas a partir de lluvias en el ro
Cauco (rea 13.000 km 2) (ref.6.

91
 El modelo completo se aplic paro lo reconstruccin de hidrogromos en lo cuen
ca del r(o Mc Kenzie (reo 4.800 Km2) con resultados satisfactorios en compara:"
cin con otros modelos existentes y demostr ser adecuado paro fines de previsin
de descargos, donde el tiempo de prediccin es esencial y los dotas sOn general-
mente insuficientes. Se debe hacer notar que el mejoramiento de los coeficien-
tes a parmetros en uno cuenco dodo es un proceso cont(nuo, y uno de las ven-
tajas del modelo es pennitir el perfeccionamiento de los tcnicas sin alterar su
estructura; el grado de refinamiento que puede obtenerse depende de la dsponib
I idad de datos bsicos, osr como de lo iniciativo de las personas en la aplicaciOn
de los d~tos al problema de s(ntesis de descorgos.

92
 CAPITULO IV

MODELOS ESTOCASTICOS

4.1 ALGUNOS ASPKTOS DE LOS EVENTOS r'IDROLOGICOS

4.1.1 Procesos hidrolgicos - definiciones

Las regiones de lo tierra donde se presentan fenmenos hidrolgicos particulares,

pueden denominarse ambientes hidrolgicos, cualquiera de los cuales puede de-
cirse que consiste en entrados de agua, respuestas del ambiente y almacenamien
to de aguo ante esas entradas, y solidos de aguo.

Los procesos hidrolgicos naturales son de naturaleza estocstico, o uno combi-

nacin determinstico - estocstico de procesos.

Los procesos determinsticos de lo hidrologa SOn el resultado de leyes fsicas,qu

micos y biolgicos, pero primordialmente de las leyes y regularidades de la me=-
cnica de fludos y lo termodinmico.

Lo mayora de los fenmenos hidrolgicos en lo Naturaleza son procesos estocs

ticas, o procesos gooemados por los leyes del ocaso. Estrictamente hablando, no
hoy procesos determinsticos puros en el mundo natural; por eso, lo hidrologa no
puede ser ntegramente comprendido, descrito y aplicado aJecuadamente, sin el
uso extensivo de mtodos basados en lo teora de prOOabili:lades. La precipita-
cin anual, en un punto o sobre un rea definido, lo descarga anual en uno es-
tacin, y otros procesos hidrolgicos simi lores de volares anuales, son procesos
estocsticos, siempre que no hoyo actividad humano artificial que influyo en -
ellas.

El estudia de lo evolucin de lo configuracin de un sistema fsico en el tiem-

po, presupone dos principios bsicos:

a - Dado un sistema fsico S, los efectos generados sobre el mismo Dar accin de
lo influencio de cuerpos exteriores en intervalos de tiempo detenninados, se
acumulan y constituyen las causes de ias variaciones de la configuracin de
S.
b - Dados dos configuraciones de un sistema Fsica S, existe uno probabilidad
asociado al poso de lo primera modificacin poro la segundo.
93
 A portir de estos principios viene la definicin de proceso aleatorio:
Se llama proceso aleatorio aquel poro el cual, o intervalos aleatorios o de-
terminados de tiempo, corresponden cambios de estado del sistema frsico S
asociados a leyes de probabilidades.

Proceso estocstico es el proceso aleatorio ms general, en el cual lo distri-

bucin que caracteriza o cado instante lo variable aleatoria del pro<:eso es
funcin del tiempo, o sea, varra de instante a instante. Paro su completo
cono<:imiento, por tonto, es necesario que se muestre de qu manero la dis
tribucin que lo caracteriza es funcin del tiempo.

Procesos markavianos son procesos aleatorios independientes del tiempo, don

de el valor de la variable aleatorio en el instante CerO determina estadrsti-:
comente el valor de la variable en el instante 1, que o su vez lo determi-
na en el instante 2 y asr sucesivamente. Este proceso es de gran importancia
en el estudio de procesos hidrolgicos, dado el grado de autocorreb::in exis-
tente en algunas series temporoles hidrolgicas.

Procesos probab i Irsticos son procesos a leatorios independientes de I tiempo, don

de la secuencia de ocurrencia de las variables involucradas en un proceso -
es ignorada y lo probabilidad de su ocurrencia se considera segn uno distri-
bucin de probabilidades determinado, en lo cual las variables se toman co-
mo puramente aleatorios (no autocorrelacionadas).

Procesos ergdicos son procesos donde lo distribucin de lo variable aleatoria

es independiente del tiempo y, adems de eso, dados n procesos iguales, el
muestreo hecho en un nico instante es estadrsticamente equivalente a un -
muestreo hecho sobre un nico proceso en n instantes distintas de tiempo.
Son, por tanto, pro<:esos sin historia, donde las estructuras fundametales son
independientes del tiempo.

El mtodo de la estacin-ao es bosodo en la ergodicidod del proceso de deter

minados tipos de precipitacin. Dadas n regiones climticamente hamogneas,-
y si las causos de la precipitacin en coda regin ro dependen de las causas en
los dems, el muestreo de los precipitaciones en 1 ao en las n regiones 01
mismo tiempo, equivalen a n muestras (o n aos) de precipitacin en una ni-
ca regin, si el proceso de las precipitaciones fuero ergdico.

4.1.2 Caracterrsticas de los datos hidrolgicos

Poro investigar las propiedodes de las variobles aleatorias con la distribucin es
pocial y temporal de lo contidad y cal idad del agua, hay cuotro tipos de datos-
disponibles:

o - Datos histricos o cronolgicos, u observaciones de procesos en e I tiempo,

de las cuales resultan series hidrolgicas temporoles discretos. Lo mayorro
94
 de los datos hidrolgicos actuales pertenecen a ste tipo;

b - Observaciones de campo a lo largo de I(neas o reas determinadas, como la

determinacin de profundidades de agua subterrneo en estratos porosos. Uno
parte sustancial de los datos actuales puede ser de este tipo;

e - Datos de laboratorio o experimentos de campo, usados principalmente en in-

vestigacin.

d - Medidas simultneas de dos a ms variables aleatorias para establecer una

relacin entre ellas, principalmente para transferir informacin estadstica en
tre ellas.

Frecuentemente es til hacer una distincin entre el valor verdadero, original y

observado de una variable hidrolgica. El valor verdadero de una observacin
nunca se conoce a causa de que el data obtenido, por inevitables errores en la
observacin, no es un valor exacto. El valor original es el producido por con-
diciones no modificadas de un ambiente; el valor observado es el resultado de
varias investigaciones, registros o experimentos; es el valor usualmente publicado
por tos servicios hidrolgicos.

Los datos observados tienen errores en cado una de los cuatro fases de deteccin,
transmisin, registro y procesamiento; los valores calculados siempre contienen -
errores ya sea aleatorios o sistemticos; en ste ltimo caso se dice que son in-
consistentes; los datos afectados por cambios ambientales, naturales o artificiales,
son no homogneos.

Los tres conceptos de errores aleatorios, sistemticos y de no homogeneidad, son

necescrios y tiles en las aplicaciones prcticas hidrolgicas. Errores aleatorios
hoy siempre en los datos; a menuda se asume que estos errores se distribuyen sim
tricamente alrededor del valor verdadero, siguiendo una distribucin normal; la .-
magnitud de estos errores se mide entonces o travs de su desviacin estndar.

La inconsistencia es causada por errores sistemticos en los datos hidrolgicos; se

define como lo diferencia entre los datos observados, con el error sistemtico in-
herente, y los valores verdaderos.

La no homogeneidad es la diferencia entre los valores verdaderos y los originales

de una variable hidrolgica; es causada por cambios introducidos en las condi-
ciones naturales.

El anlisis de estos errores y sus efectos en los propiedades de las series, es in-
dispensable cuando se requiere conocer el grado de confabilidad de la informa-
cin estadstica hidrolgica.

95
4.1.3 Naturaleza de los magnitudes hidrolgicos

Es conveniente hacer uno distincin entre los trminos "fenmeno hidrolgico"

(precipitacin, evaporacin, caudal, transporte de sedimentos, etc.), "variable
hidrolgico" que describe esos fenmenos (descargo instantneo, mensual, anual,
etc.), y"parmetros hidrolgicos estadsticos", que describen los funciones de pro
babilidod y procesos estocsticos de esos variables (medio, mediano, varianza, cae
ficiente de variacin, coeficiente de asimetra, coeficientes de correlacin de :
serie etc.). En lo mismo formo debe distinguirse entre uno varJoble bsica (cau
dal, intensidad de lluvia, etc.), y una variable derivado (valores diarios, mensuo
les, anuales, picos, mnimos, volmenes totales, etc.).

Uno de los objetivos de los numerosos anlisis efectuados con lo aplicacin de la

teora de probabilidades, matemtico estadstico y procesos estocsticos, es el -
desorrollo de procedimientos prcticos paro determinar los caractersticas de las
variables aleatorios derivadas de las variables aleatorios bsicos.

4.2 APLICACIONES DE ESTADISTICA EN HIDROLOGIA

Se incluyen en este captulo algunos conceptos elementales de estadstico, con

el nico objeto de unificar criterios y aclarar las proposiciones que se haron ms
adelante; es obvio que para mayores detalles se deber recurrir o algunos de los
numerosos textos existentes sobre lo materia.

La mayor parte de los datos hidrolgicos disponibles actualmente conforman uno

serie de muestros estadsticos de tamao limitado; por tanto, lo estadstico mote
mtico es lo principal herramienta poro extraer de esos muestras lo mayor centT'
dad de inFormacin posible, y deducir conclusiones sobre las caractersticas de-
las poblaciones hidrolgicos de donde provienen esas muestros.

Los principios, posibilidades o limitociones de 105 mtodos estadsticos en el com

po de la hidrologa, son ms importantes que las complejos tcnicas para opli--
carlos o, en otras palabras, no son tan importantes los tcnicas utilizadas, como
las condiciones bajo las cuales ellas son vlidas.

Los principales objetivos que se persiguen con lo aplicacin de los principios e~

todsticos o la hidrologa se pueden resumir as,

o - La interpretacin de los observaciones;

b - Lo bsqueda de regularidades probabilsticas en los datos hidrolgicos;
c - Lo extraccin del mximo de informacin de los datos hidrolgicos;
d - Lo presentacin de la informocin hidrolgica en form~ condensada, median
te nmeros, ecuaciones, grficos y parmetros de grcrl aplicoc in sobre tOOo-
96
 en planeacin, conservacin, desarrollo y control de recursos hdricos.

El uso indebido de la estadstica en hidrologa puede causar interpretacin err-

nea de los datos en los siguientes casos:

a - Esperar mucho mas de lo que los datos y mtodos estadsticos pueden produ-
cir;

b - Creer que los mbdos estadsticos pueden de alguna manera remplazar la es-
casez y falta de confiabilidad de los datos, lo cual es impasible, y

c - Intentar remplazar la necesaria interpretacin fsica de las fenmenos por

proposiciones o clculos estadsticos.

4.2.1 Nociones generales

Una reunin de individuos constituye una poblacin; debe entenderse por indivi
duo uno entidad que posee o no existencia propia; un hombre, un vegetal, uno-
presin atmosfrica, una temperatura, un caudal, una estrella, san ejemplos de in
dividuas; el conjunto de ltis estaturas de los hombres, las dimensiones de las es-=
tre 1105, las presiones en Bogot a las 06.00 horas, por e jemplo, son las poblacio
nes. El estudio de estas poblaciones no se puede hacer analizando todos los i
dividuos, puesto que el nmero de stos es usualmente infinito \cama lOS'Couda=-
les de una estacin dado), o fi"ito pero muy grande (cama el nmero de habitan
tes de un pas), o an indeterminada (coma el nmero de estrellas); es imperioso,
por tanto, estudiar un grupo d= individuas relativamente restringida, que constituye
uno muestra de la pablocin, la cual debe ser extrada de tal manera que sea-
una imagen lo ms fiel posible de la poblacin de donde proviene; para que es-
to sea as, el nmero de individuos debe ser suficientemente grande y, por otra
parte, los miembros deben ser escogidos al azar entre la poblacin. De aqu se
deduce que una extraccin premeditada de ciertos individuos, tender a "ampli
ficar" los caractersticas de la poblacin comunes a esos individuos, mientras qe
dejar de lado otras aspectos que poseen los miembros excludos.

Los principales caroctersticas de una poblacin estn condensadas en una cier-

fa cantidad, usualmente pequea, de nmeros denominadas" parmetros", las cua-
les dan una descripcin esquemtica de una poblacin, los parmetros de una
poblacin na son directamente accesibles por no serlo la poblacin misma; dben,
por tonto, encontrarse a travs de una muestra de la poblacin.

Una muestra est descrito, igualmente, por sus parmetros media, desviacin tr-
pico, etc.), los cuales son calculables puesto que la muestra es totalmente cono-
cido, pero Jo que interesa conocer son los parmetros de la poblacin. En con-
secuencia, ya que la muestra es una imagen de la poblacin, se dice que lospa
rmelros de la muestra son estimaciones de los correspondientes a la poblacin-
falta sober, entonces, con qu precisin ellos representan los valores reales de

97
 los parmetros de lo poblacin, lo cual se hace trovs del estudio de los le-
yes de distribucin de los estimadores.

4.2.2 Distribuciones empricas (Distribuciones de frecuencia)

Los observaciones de uno variable pueden estor en formo de uno muestro de va-
lores discretos o de un registro contnuo.

Poro una muestro de N volares Xi de uno variable discreto X, hoy n valores

de un evento K. Estos ni nmeros son los frecuencias absolutas de X eventos;
los relaciones F = n/N son sus frecuencias relativos. El grfico de frecuencias
muestro cmo los frecuencias (absolutos o relativos) se distribuyen sobre estos even
tos discretos, y se representon usualmente en diagromas de barros. lo sumo de las
frecuencias absolutos, ni' es N, y lo sumo de los frecuencias relativos es l. La
frecuencia acumulado absoluta, N, de un evento discreto,est dado por

N = ~ ni
i =1

y lo frecuencia ocumuloda relativo, F (X ), por

i
F (X i) = P (X ~ Xi) =

1 L ni
i = 11

El rango de una muestra es lo diferencio entre el mayor y menor volar observa-

dos; veces es til clasificar los volares Xi en subgrupos, denominadas interva
los de clase; el nmero de valores que caen dentro del intervalo es su frecue
cia absoluto de clase, ni; usualmente el volar caracterstica de codo clase es su-
punta media.

La frecuencia relativo, f, del intervalo' es lo relacin de su frecuencia abso

luto n 01 tamaa de lo muestra, o seo

f = n/N

lo frecuencia relativa acumulada para el volar Xi es F (X) = N/N.

La grfico a curva de frecuencias describe lo manero coma los frecuencias se
distribuyen sobre sus intervalos de clase; paro evitar confusin, es necesario ocio
ror siempre si se troto de frecuencias absolutas a relativos (densidodes de frecuer.
cia). Si los frecuencias de cada intervalo se representan par barros horizontaleS,
el grfico se I lamo histograma de frecuencias; si los puntas medios de cada inter-
valo se unen con lneas rectas, se obtiene el polgono de frecuencias.

98

I
-~ ..
 Si lo longitud de los intervalos tiende o cero, el polgono y el histograma de
frecuencias tienden o uno curvo lmite ideal, denominado Curvo de frecuencias.

4.2.3 Parmetros estadsticos

Como se dijo antes, los parmetros permiten dar uno descripcin cuantitativo es-
quemtico de lo poblacin estudiado; de uno muestro, sin emoarga, podemos o la
sumo deducir estimaciones de esos parmetros. El conocimiento de la ley de re-
particin de esas estimaciones permite determinar el intervalo en el cual estn
comprendidos los parmetros de lo poblacin.

Lo Curva de frecuenc ias asoc iodo a una muestra puede carac terizarse por:

o - Su posicin; determinado por los denominados promedios 'medio, mediano,

modo, medio geomtrico, medio armnico), o medidos de tendencia centrol;

b - Lo dispersin de los diferentes observaciones 01 rededor del volar medio; los

medidos de dispersin ms generolmente empleados son lo desviacin tpico,
el coeficiente de variacin y lo desviacin inter-cuartil.

c - Su asimetra o sesgo.

A. Medidos de tendencia central

Se llaman as porque tienden a situarse en el centro del conjunto de datos

ordenados segn su magnitud.

La medio aritmtica es un valor alrededor del cual los valores observados

varran de uno manero aleatorio; se represento por X y se define como
n

L X
I
=
x =--'----
n

donde n es e I nmero de datos

Se llamo Mediana Me el valor de la variable tal que el nmero de ob-

servaciones poro las cuales el volar de X es inferior o "Ae es igual 01 nu
mero poro el cual el volar de X es superior a Me.

El Modo es el valor de la variable correspondiente 01 mximo de la cur-

vo Clel"recuencia; es el volar que se presento ms frecuentemente.

99
 Lo Medio Geomtrico G de uno serie de n volares Xl, X2 ,X n ,
est definido por
l/n
G = (XIx X2 x X3 ---x Xn )

o, definido de otra manero, el logaritmo de lo medio geomtrico de uno

serie de volares es el volar medio de lo sumo de sus logaritmos.

Lo Medio Armnico H est definido por

1
"R"" = n1 '"
L- (+)
n

donde ni es lo frecuencia del volar Xi

B. Medidos de dispersin
Los medidos de dispersin no slo ayudan o describir uno distribucin, sino
que tambin praporcionon informacin sobre lo precisin con que lo medio
de uno muestra caracteriza lo media de la poblacin.

La desviacin tfpica (o desviacin estndar) se define como la rafz cua-

drada del cuadrada medio de las desviaciones:

Se denomina varianza (02) al cuadrado de la desviacin estndar.

El coeficiente de variacin es un parmetro adimensional definido por la

relacion entre la desviacion estndar y la media

cJ
=
X

C. Medidas de asimetrfa
El coeficiente de asimetrfa (o de sesgo) es un fndice de la falta de sime
Ida de una distribucion; se define por

L (X - Xl 3
a
~ = ______________
n i =1 __
n
21I 1,5
~ (X -)()
100 i = 1 J
4.2.4 Correlacin y regresin

Se define correlacin como la asociacin de dos o ms variables aleatorios que

explica parcialmente la variacin total de uno variable por la variacin de otros
variables aleatorias involucradas en esa asociacin. Lo correlacin es, pues, un
proceso que define el grado de asociacin entre dos o ms variables.

Esta asociacin tiene dos extremos: el caso en que hay correlacin total, equiva
lente a una relacin funcional, y el caso en que hoy ausencia de cualquier co:
rrelacin.

La Regresin represento uno ecuaclon matemtica, que expresa la correlacin en-

tre dos o ms variables aleatorias. La ecuacin de regresin puede ser cualquier
funcin factible de aiustarse a una nube de puntos.

Si dos variables dadas como una serie de volares concurrentes IXi, vi; muestran
una concentracin alrededor de una curva imaginaria al graficarlos ver figura 36),
siempre habr una distribucin de Y valores para un valor dado Xi. El valor me
dio Yo' de todos los valores Y para este intervalo X alrededor de Xi , es el va
lar esperado de Y para el valor dado de X- Xi. Una curva aiustada a todos 105-
valores medios Yo, se denomina lnea o curva de regresin de Yen X. Por otra
parte, la curva aiustada a todos los valores esperados X, para un dado Y=Yi ,de
fine la regresin de X en Y. Estas dos lneas no coinciden y tienen parmetros
diferentes.

En hidrologfa son comunes dos tipos de regresin

a - Relaciones causa-efecto, donde una variable Y se relaciona correlativamen

te o factores causales Xi, los cuales producen o afectan el resultado Y; --
e iemplos tpicos son las relaciones precipitacin - escurrimiento, puesto que
la lluvia es la causo bsica del escurrimiento, con factores como la geome
tro de la cuenca, suelo, humedad, clima, etc., que afectan la relacin be
sica causa-efecto;

b - Relaciones de variables aleatorias que tienen por causa los mismos factores,
tales como la asociacin correlativa de las descargas de ros adyacentes.

El coeficiente de correlacin es una definicin matemtica de la asociacin

entre las variables, sin que esta definicin implique, como ya se puntuali-
z, relacin de causa a efecto. Asumamos que se tienen n observaciones
simultneas de las variables aleatorias X e Y, y denotamos por (Xi, Vi) el
i-sima par. El coeficiente de correlacin r est dado por
lo n

L (X-X)
I
(Yi -YI Xi Yi-n XY
i=l i-l
r = =
n x 0y n x uy

101
 y
(l i t yi.)

y=I(.)

i i
'i yi

I I
I

x
CONCEPTO DE REGRESION
FIGURA 36

1(.)
0.8

a- =0.5

a- = 1.0
/
03

a- = 2.0
0.1

-3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.5 1.0 1.5 2.0 2.5 3.0
x
CURVAS DE DENSIDAD NORMAL DE PROBABILIDAD
FIGURA 37

102
donde Ux y ay son las desviaciones trpicas de X e y,

Es conveniente resaltar lo diferencia entre 105 trminos correlacin y regre-

sin, yo que lo primera est representada por un coeficiente, mientras que
lo segunda do una ecuacin que permite calcular los valores de una varia-
ble en trminos de la otra.

Los datos usados en un onlisis de correlacin deben ser razonablemente horno

gneos; si la nube de puntos muestra un agrupamiento de puntos en dos o mas
zonas, puede resultar una alta con:elocin aparente, pero debida solamente a
la heterogel'\eidad de 105 datos; tales correlaciones se llaman falsas o apare~
tes.

Cuando se tiene una serie temporal hidrolgica y se desea saber el grado de

asociacin que existe entre pares de volares separados k unidades de tiem-
po, se usa el coeficiente de correlacin de serie o coeficiente de autoco-
rrelacin, calculado en la misma forma indicada atrs, con los valores de X
e y siguientes:
X Y

Xl Xl+ k
X2 X 2+ k
X3 X3 + k
......
donde k es un nmero entero llamado perrada o retardo, el cual debe espe-
cificarse siempre que se habla de este tipo de correlacin.

En el caso de caudales, ste coeficiente indica el grado de persistencia, es

decir, la dependencia que presenta un caudal, de otro situado k unidades
de tiempo atrs. Un ejemplo de persistencia es el hecho observado de que
a caudales altos siguen generalmente caudales altos y a caudales bajos si-
guen tambin caudales bajos.

El valor del coeficiente de autocorre locin de orden 1 puede tomarse como

un indicador de la cantidad de informacin agregado al estudio de un proc~
so, por coda valor de lo variable obtenida en observaciones.

Esta cantidad de informacin varra proporcionalmente al inverso del coeficien

te de outocorrelacin de orden 1. Asr, para un valor nulo, la cantidad de-
informacin es mximo, y nula para un valor unitario, en cada nuevo valar
anexado a las observaciones. Asr, una serie temporal aleatoria (coeficiente
de autocorrelacin nulo) de una determinada extensin, presenta mayores in
formaciones que una serie temporal no aleatoria (coeficiente de autocorrela
cin no nulo) de igual extensin.

103
4.2.5 Distribuciones de probabilidad
Casi todas las distribuciones de variables aleatorias hidrolgicos son de naturale
za emprica, obtenidas a partir de datos observados muy limitados, y su mejor-
descripcin se logra con el apropiado ajuste de funciones de distribucin de pro
babilidad. -

Las funciones de distribucin de probabilidad de las variables hidrolgicas bsi-

cas no pueden, en general, deducirse tericamente, y las funciones verdaderas
nunca se conocen; lo que se hace en lo mayora de los casos es seleccionar una
funcin y estimar sus parmetros o partir de lo muestra. Si se encuentro un ajus-
te aceptable, se asume que la misma funcin se apl ica a la poblacin de esa va
riable.

Las funciones de distribucin de probabilidades se pueden clasificar desde diferen

tes puntos de vista, segn se considere el aspecto estadstico o el fsico; sin em- -
bargo, se pueden dividir, en general, en funciones tericas, semi-tericas y em-
" .
plflcas.

En hidrologra es posible encontrar ejemplos de aplicacin de todas las funciones

de distribucin de probabilidad, tales como la hipergeometrica,binomial, de Poi-
55on, y otras. La experiencia ha mostrado, no obstante, que slo en limitado n
mero de ellas se puede aplicar adecuadamente a distribuciones hidrolgicas empi
ricas, y entre ellas estn la normal (de Gauss), lognormal (de Galton), Gamma-
(particularmente en su forma de dos parmetros y la Pearson tipo 111).

Si se tiene en cuenta el hecho de que la mayora de los registros de caudales,

especi'ficamente, son sumamente cortos, las pruehas estadsticas rigurosas que exi
gen gran cantidad de datos para estimar la bondad del ajuste no constituyen e
este caso una ayuda significativa para seleccionar la distribucin ms adecuada.
Se denomina bondad del ajuste al grado de precisin con que los valores de una
variable aleatoria se ajustan o una distribucin determinada. En la mayor parte
de los casos la seleccin de la distribucin se basa en la experiencia; a menudo
se sabe que ciertos tipos de distribucin se ajustan bien a determinadas variables
hidrolgicas; es corriente, tambin, que dos o tres funciones se ajusten igualmen
te bien. -

Las distribuciones que, de acuerdo con la experiencia acumulada en el anlisis

de series de caudales, se ajustan mejor a estos datos, son la normal, lognormal
con 2 parmetros y gamma con 2 parmetros, por lo cual se dar una rpida des
cripcin de coda una. -

A. Distribucin normal

La funcin de densidad normal de probabilidad de Gauss", oe Gauss-Lapl~

ce) de una variable X es 'ver figura 37
104
 f (X) = e - (X"1!) 2/2 ('2
6VZTT
donde Jl es el valor esperado (o medio) y ()2 lo varianza de lo variable
X; lo fo""o dado en lo ecuacin anterior se simbolizo por N (.)J. (J ),
con pormetros estimados o porti r de lo muestra disponible.

Con lo variable normalizado o reducido t = (X -./" ) / () , la ecuacin

es
f () =
V2Tr
y representa una distribucin con .JI' = O Y r:r 2 = 1, simbolizado por N
(O, 1).

Se ha encontrado (ref. 24) que lo distribuc in norma I se a justa bien a distn

bueiones empfricas de precipitaciones anuales sin embargo, la principal di
ficultad terica para el ajuste de esta distribucin estriba en el hecho de -
que una variable no""ol puede asumir cualquier valor entre - c:t) y + ct:> ,
mientras que lo moyorro de los variables hidrolgicas asumen solo valore$
patitvos; sin embargo, si la medio es mayor que 3 \ o 46"", lo probabili-
dad de que el valor X sea negativo es muy pequea.

Los mtodos ms sencillos para dete""inar si una serie de caudales presen-

to o n una distribucin normol son el papel de probabilidades y el alculo
del coeficiente de asimetrra.

Con el papel de probabilidad, se groficon los valores de los caudales con-

tra el porcentaje de valores iguales o mayores a ellos, y se ajusto manual
mente una curva; si sto se aproxima o uno recto, se puede concluir que-
los dotas se pueden representar por una distribucin no""ol.

Adems, si consideramos que el coeficiente de osmetrra de uno distribucin

normal es cero, mientras ms se acerque o este valor el coeficiente de los
caudales histricos, ms se aproximo su distribucin a lo no""ol.

B. Di s tri b uc in lag - no r mo I

Cuando los logaritmos, In X, de uno variable X estn normalmente distri-

buidas, se dice que lo variable X sigue uno distribucin de probabilidades
lag-normal; est representado por lo ecuacin

f (X) = exp [ _~ (ln X ~n

X<nl/ 27r' (n

f (X) = O para X <:. O

105
 donde el ndice n en./-'- n y o n indico que son deducidas de los 1090-
ritmos naturales, y/A n es la media de 1n X as como 0n es su desvia
cin estndar.

Se ha encontrado a travs de uno amplia experimentacin, que la funcin

lognormal se ajusta bien a un gran nmero de distribuciones en hidrologa,
a causa de ser asimtrico y por otras de sus propiedades (ver figuro 38), en
tre las cuales se destaco la de no producir caudales negativos por ejemplo"';
por cuanto los logaritmos no estn definidos poro nmeros negativos; sto lo
hoce muy recomendable cuando se hacen estudios en el ranga de caudales
muy be jos (sequas).

Para determinar si una serie de caudales histricos presento uno distribucin

lognormol, se usan los mismos principios recomendados poro lo normal, pero
trabajando con los logaritmos de lo variable.

c. Distribucin Gamma

Lo integral
1 e -X d
x
converge si c( > O Y se llamo la funcin gamma de o{ , y se puede calcular
por expansin de series e integracin numrico.

La variable continuo X, cuyo funcin de densidad de probabilidad es

1 ri- 1 e- x/p
f (X) = X
/3 r (e(' )
se dice que sigue uno distribucin gamma de dos parmetros, donde ex es
el parmetro de forma (ver figura 39). Su uso es tan comn en hidrologa
como el de la 109Oormal.

Mancovic (ref. 24) demostr que prcticamente no hay diferencio en el oju,!

te de distribuciones lognormal y gamma a distribuciones empricas de preci-
pitaciones y caudales anuales.

D. Esquema para seleccin de lo distribucin

Fiering y Jacksan (" Synthetic Streamflow". 1971) establecieron la secuencia
de pasos que deben seguirse en lo seleccin de la distribucin terica que
ms se ajuste a una distribucin emprico de caudales;
a - Examen de los registros histricos, para detectar posibles cambios origina
dos por la actividad del hombre, o si los caudales han sido regulados e
el pasado.
106
 .. 2
f (1 )

- - - -(J"' 0.5
---1\'0

0.4

-- ..... _-
2 , 1 4

Seis curvos de probabilidad lognormol , '" n = 0.05, Y LOO, poro IT ~ = 0.5,

Y IT *= 0.1, 0.5 Y 2, poro ",. = O.

FIGURA 38

, ( 1 )

1.00
<1'2
0.75

0.50 a, I
~, 2
025
X

0.5 1.0 2.0 3.0 4.0 5.0 6.0

Curvos de densidad de probabilidad gamma 2 poro ~= 2,

Y tres valores de f3 (O.50, 1, Y 2)

FIGIJ'flA 39

. . . . . . . . . d. ",,"-"'UI ... SteU,t... le H 'l Vu . . le" Fort CoUiltl Colltt'.Oo 1972

107
 b - Clculo de lo medio X de los caudales poro un intervalo de tiempo
considerado importante, (meses, aos, etc).
c- Clculo de lo desviacin estndar de los dotos histricos.

d- Clculo del coeficiente de asimetra tx de lo muestro.

e - Clculo de los logaritmos (Yi,) de los datos, as como su media, des
viocin tpico y el coeficiente de asimetra.
f - Graficar los valores de Xi en el pope I de probabilidades poro verifi-
car si se puede ajustar una distribucin normal; el clculo de x ayu -t
da a sto verificacin.

g - Si del poso anterior se deduce que los caudales no estn distribuidos

normalmente, se repite el paso f poro los volares de y i; si y i po-
rece normal, es adecuado la distribucin lognormol.

h - Si ninguna de los suposiciones anteriores se confirmo, son posibles va-

rias oltemotivos. Si el coeficiente de asimetra es pronunciado y se
deseo que ste seo reflejado en los caudales sintetizados con base en
uno distribucin, puede usarse la gamma. De otra forma, deben usar-
se tcnicas tericas de decisin pora elegir la distribucin ms apropio
do. -

i- Clculo del coeficiente de autocorrelocin de orden 1 que puede ser

calculado por

n-) n-) n
X.I Xi+) ( 2 X;) (2" x)
n-)
i=1 i=) i=2
rl =

[~ O,5[~2
r
n- ) n
2
X , I ('Z'
n:r .,-) Xi) 2]
2
XI -~(
n ,-. 2
x)2
i

donde Xi son los caudales de la distribucin emprica.

4.3 SELECCION DE UN MODELO PARA SINTETIZAR CAUDALES IreL 23)

Lo esencia de lo simulacin es la reproduccin del comportamiento de un sistema

dado en codo uno de sus aspectos importantes. As pues, en los programas de si
mulocin de sistemas de recursos hidrulicos de propsitos mltiples, se incluyen-
variables tales como el dficit mensuol de energo y de irrigacin y las dencien
cias acumulados a lo larga del perrado de estudio, los oven idas mximos, los cos
tos de irrigacin, control de avenidas, generacin de energa, costos de capital,-
mantenimiento y operacin, beneficios netos, etc.
108
 Pero, si bien es cierto que la simulacin ofrece la ventaja de poder analizar el
comportamiento de un sistema bajo muchas condiciones de operocin, por otro
porte presento ciertos limitaciones, ya mencionados en el capftulo 1:

o - No lleva directamente o uno solucin ptimo, aunque despus de probor va

rias combinaciones de variables, es posible escoger entre estas la que pro-:
duce mayores beneficias, con mayores ventajas que los mtodos convencio-
nales;
b - Existe falta de flexibilidad en los procedimientos de operacin del sistema,
por cuanto paro cado forma de operacin, es necesario un progroma dife-
rente;

c - En general faltan registros de caudal lo suficientemente extensos paro que

puedan considerarse como uno representacin de la posible configuracin de
los caudales.

Todas los limitaciones tienen soluciones parciales, y aqu se consideraron slo

las correspondientes a registros de caudales.

Puede darse el coso de registros que, por su corto extensin, no incluyan pero
dos cr(tcos de caudales muy altas o muy bajos. Sin embargo, o menos que eT
registro sea sumamente corto,generalmente es suficiente para determinar parme-
tros estad(sticos confiables que puedan usarse dentro de un modelo generador de
caudales.

As se puede obtener una serie de caudales de lo extensin deseada, que repro

duce las principales caracter(sticos de lo muestra, y que puede ser usada en ro
simulacin del sistema completo.

los caudales son una de los principales variables del modelo simulador; si se
generan, par ejemplo, 500 aos de datos, y el perodo de diseo es de 50 aos,
se tienen 10 posib les muestras paro usor.

4.3.1 El mode I o de Ma rkov

Entre los cosos extremos de los modelos puramente determinsticos y los purarnen
te estocsticos, el modelo de Markov ocupo uno posicin intermedio; en stos T
timos, el estado de un sistema es completamente independiente de estados ante=
riores, y obedece nicamente o los leyes de lo probabilidad; tiene mayor ven-
taja que los modelos determinfsticos en aquellos fenmenos donde los procesos
bsicos son complejos y sujetos a influencias que no pueden ser evaluadas exac-
tamente.

El modelo de Markovocupa una posicin de de~ndencia p?rcial entre los dos

extremos de dependencia completa y completa independenCia; lo formo general
109
 perrexlo siguiente seo alto, ir en su mayor parte o 105 ocuileros y se manten
dr lo tendencia de 105 caudales bajos.

B. Caudales con distribucin lognormal

Cuando los caudales histricos tienen sta distribucin, se tendr uno secuen
cia de valores de logaritmos hi, h2, h3 ... hn y la ecuacin de los caudales
generados se r
qi = exp (h i)

Debe recordarse que este procedimiento reproduce los parmetros de los lo-
garitmos de los caudales, pero no preserva el coeficiente de asimetrra y el
de correlacin de serie de 105 caudales.

c. Caudales con distribucin ,gamma

lo formo general de la ecuacin

puede adoptarse a caudales que presenten distribucin gamma, si se introdu

cen

-tx como coeficiente de asimetrra de lo muestro

y "lE. como el coeficiente de asimetrra de la variable aleatorio.

~ se define como

~= 1_;3 ?t-x
(1-/2) 1.5
Si ti tiene distribucin N (0,1), se puede usar lo transformacin:

i = 2 nt?: ~2 3
- (1+ --)
~E 6 36
y dor valores aleatorios con distribucin gamma con media cero, varianza
unitaria y coeficiente de asimetrra '1-",.

4.3.2 Modelos estacionales de Markov

los modelos anteriores se recomiendan cuando se necesiton d

les, pero cuando es necesario obtener caudales estacional cau a es anuo-
112 es o mensuales ,

,L_ 'pg
rn
el mismo modelo de Markov que se seguir usando se toma ms complejo.

Es lgico que cuando se usan caudales estacionales, no pueden usarse una media
y varianza nicas puesto que estos parmetros difieren de una estacin o otra; y
sto hace necesario el uso de dos subndices, el primero de los cuales indica el
afia de ocurrencia del evento y el segundo la estacin dentro del afio dado; el se-
gundo vara de 1 a n y se repite dclicamente, siendo n el nmero de estacio-
nes de cada ao.

As, }-'- indica la media de lo estacin j-sima a lo larga de todos los m aflos
de registro, es decir, la suma de los caudales de la estacin y, dividida entre
m lo misma sucede con la desviacin estndar () i Y con JI', que es el cae
ficiente de correlacin entre pares de volares correspondientes o los estaciones-
(j) y (-l).

A. Caudales con distribucin normal y lognormol

En el primer caso el modelo generador se convierte en

~'()' 2
qi,i=/j + '~j_: (qi,i- 1 -~i-l)+ ti,jCfJ(1- )!
y do el caudal en el ao i-simo y en la estacin j-sima, compuesta par:

- lo media de los caudales en lo estacin j-simo,

- lo diferencio entre el caudal anterior de lo estacin j-1 y su media, mul-
tiplicada por el coeficiente de correlacin y par (J /cr;-L
_ t.1, j, lo variable aleatoria nonnolmente distribuido.

Cuando el ndice j tiene el valor 1, el volar j-l debe sustituirse por n

(12, en modelos mensuales), y se debe disminuir el ndice i en 1.

Para caudales con distribucin lognonnal se sigue exactamente el mismo pro-

cedimiento, con sus logaritmos.

8. Caudales con distribucin gamma

Se tienen los parmetros ./'" j, <Jj,!'j definidos anterionnente se calcu

lo 1'" j, coeficiente de asimetra de los caudales de lo estacin j, Se de
fine ~ t, j' como
[;tj _/0-1)3 Yj-l]
=
( 1 -ti 2 ] 1.5

113
 Y se usa la transformacin t1,
r, " a partir de ti' (normal) :

2 t t ti . (/2t, i )3 2
ti, .... I i = ( 1+ " " 36 - tt,;
J t, i 6

Asr ti 't- es una variable con distribucin gamma oon media cera, des
'." ~ 'do r y cae f"IClente de asnnetna
vlaclon estan . ~ J o.
4' t, '1
SI se usa ti,'" " en 1-
u-
gar de 'a\,
i en el modelo generador,
p. 6'"'
q' . "',u. + / ,
'ti , / , cJ-l
, (q, ;-1 - / H ) -t-I, ~, i o (1- i 2) ~

producir caudales con distribucin gamma.

114