CANARINA:

Principal

Aire

Agua

Sonido

Electromagnético

¿Quiénes somos?

Testimonios

Preguntas

Contacto

DEMOS

Precio

Prensa

Mundo

 

DESCAR:

DESCAR

Aplicaciones

Ventajas

Imágenes

Algoritmos I

Algoritmos II

Algoritmos III

Algoritmos IV

Algoritmos V

Datos I

Datos II

Promedio

Funciones

Picnoclina I

Picnoclina II

Gráficas I

Gráficas II

Contaminantes I

Contaminantes II

Contaminantes III

Contaminantes IV

El mar I

El mar II

Mapas Google

3D y 2D

SIG

 

APLICACIONES:

1 emisarios submarinos y contaminación marina

2 contaminación DBO y BOD 

3 vertidos urbanos e industriales

4 emisario submarino

5 efluentes

6 vertido de residuos

7 vertido de aguas fecales

8 vertidos al mar

9 salud ambiental

10 diagnóstico ambiental

11 administración ambiental

12 información ambiental

13 contaminación hídrica

14 desagües de vertidos

15 contaminación marina

16 aguas residuales

17 plantas de tratamiento

18 aguas residualesles

19 contaminación del mar

20 estudio de impacto ambiental

21 contaminación del agua

22 proyectos ambientales

23 ingeniería medioambiental

24 contaminación del mar

25 depuradoras y aguas residuales

26 gerencia ambiental

 

 

 

 

 

 

 

 

Algoritmos V · software · información ambiental

                   

 

2.4 Cálculo de las concentraciones

 

La concentración en un punto de la pluma determinado por sus coordenadas X,Y y Z vendrá dada por la expresión:

 

  C(X,Y,Z)=(C0/S) F0(t)F1(t)F2(Y,t)F3(Z,t)

 

C0=Concentración de contaminante en el efluente

S=Dilución inicial de la capa de mezcla

siendo t una variable auxiliar que se calcula mediante t=X/Ua. F0(t) tiene en cuenta la autodepuración de la pluma no conservativa y tiene la forma:

   F0(t)=10-t/T90    (**)

El resto de las funciones tienen forma diferente dependiendo del punto de aplicación:

 

(a)  Zona próxima al punto de surgencia:

Aquí el espesor de la capa de mezcla es inferior a la profundidad y, por tanto, existe dispersión vertical. El perfil de concentraciones tiene forma de meseta.

 

F1(t)=1

F2(Y,t)=(1/2)[erf[(B/2+Y)/(σy21/2)]+ erf[(B/2-Y)/(σy21/2)]]

F3(Z,t)=(1/2)[erf[(e+Z)/(σz21/2)]+ erf[(e-Z)/(σz21/2)]]

 

Siendo σy=(2Kyt)1/2 y σz=(2Kzt)1/2. Las funciones error erf las evalúa numéricamente el programa mediante integración. El método de integración que usaremos es mediante la descomposición de trapecios de la función área de la curva, el más sencillo que hay. El número de rectángulos a considerar lo podremos variar mediante el parámetro N_int de la función PARAMETROS DEL CÁLCULO del programa. Aumentando su valor, aumentamos la convergencia numérica de las integrales pero ralentizamos el programa. Se deberá tomar un valor lo óptimo (si lo aumentamos más no notamos cambios en los resultados) para dicho parámetro.

 

(b)  Zona alejada al punto de surgencia:

Aquí se supone que la pluma se homogeneizó en vertical cuando la profundidad era Hh que es la profundidad en el punto donde el espesor de la pluma empieza a ocupar toda la capa de agua. Tenga en cuenta que el programa no permite introducir relieves topográficos en el fondo del mar por lo que Hh coincide con la altura a la superficie del agua del emisario submarino más profundo que tenemos. En caso de que quisiera considerar una situación de relieve en el que la altura del punto de surgencia Hh fuese mayor podría introducir otro emisario, con emisión nula, con una profundidad igual a Hh en cuyo caso el programa podría realizar el cálculo tomando como Hh la altura a la superficie del agua del emisario submarino más profundo que tenemos (que hemos introducido artificialmente). Aproximaremos las integrales de la función error erf de tal manera que

 

F1(t)=(2π)-1/2B σy-1/2

F2(Y,t)=exp[(-Y2/2σy2)]

F3(Z,t)=e/Hh

 

Siendo ahora σy=(B2/16+2Kyt)1/2. El cálculo hidráulico del emisario en el caso de que quiera tener en cuenta pérdidas por fricción se podrá realizar por los métodos de tanteo que sugiere la orden ministerial.

 

 

Notas sobre la modelización (errores y comentarios del modelo):

 

(*) Hemos encontrado lo que creemos, a nuestro juicio, errores tipográficos en la Orden del 13 de Julio de 1993 del Ministerio de Obras Públicas y Transportes del Reino de ESPAÑA, B.O.E. Martes 27 de Julio de 1993, página 22861, que hemos subsanado atendiendo a criterios de consistencia matemática. El software asume las presentes correcciones en el cálculo.

 

NOTA PARA LAS VERSIONES DESCAR 3.0 e INFERIORES:

 

(**) En el modelo oficial, T90 se expresa en horas, tal y como asume el programa. Sin embargo, la ecuación F0(t)=10-t/T90  del modelo oficial(1) viene con el tiempo en segundos. Atendiendo a criterios de consistencia matemática y de coherencia de resultados tendría que venir T90 expresada en segundos (multiplicar por 3600 segundos que tiene una hora). Dejamos al usuario la posibilidad de usar el modelo oficial tal y como está, o bien introducir la rectificación correspondiente. Por ejemplo, para un T90=2 horas el usuario debe introducir:

(a) Atendiendo al modelo oficial:

1/T90=0,5 horas-1 como dato de entrada del programa. Es decir, introducir 0,5 en la casilla del programa para un T90=2 horas.

(b) Atendiendo a criterios de consistencia matemática y de coherencia de resultados:

1/T90=1/(2 x 3600)=0,000278 como dato de entrada del programa. Es decir, introducir 0,000278 en la casilla del programa para un T90=2 horas.

(1) Orden del 13 de Julio de 1993 del Ministerio de Obras Públicas y Transportes del Reino de ESPAÑA, B.O.E. Martes 27 de Julio de 1993, página 22861.

(***) En el modelo oficial, no encontramos la relación que hay entre Qb y Q, ya que debería tener un coeficiente n (número de bocas). En la ventana de parámetros del modelo hemos añadido esa posibilidad.

 

Referencias Científicas:

Orden del 13 de Julio de 1993 del Ministerio de Obras Públicas y Transportes del Reino de ESPAÑA, B.O.E. Martes 27 de Julio de 1993, página 22861, I. Disposiciones generales. Proyecto de conducciones de vertidos desde tierra al mar.

IKSR 2000 : M. Braun, „The Pathways for the most important hazardous substances in the rhine basin (during floods)“, International Commission for the Protection of the Rhine, Koblenz, Germany, in Int. Symposium on River Flood Defence, Kassel,, Kassel Reports of Hydraulic Engineering No. 9/2000 

 

 

 

 

 

información ambiental

 

  información ambiental

 

DESCAR

DE14  DE15  DE16  DE17  DE18  DE19  DE20  DE21  DE22  DE23  DE24  DE25  DE26 

 

APLICACIONES DESCAR

emisarios submarinos y contaminación marina  contaminación DBO y BOD  vertidos urbanos e industriales  emisario submarino  efluentes  vertido de residuos  vertido de aguas fecales  vertidos al mar  salud ambiental  diagnóstico ambiental  administración ambiental  información ambiental  contaminación hídrica