Modelo QUAL-UFMG

Dados de entrada

Temperatura

Vazão do rio

Oxigênio dissolvido de saturação

Oxigênio dissolvido do rio

DBO5 do rio

Nitrogênio orgânico do rio

Nitrogênio amoniacal do rio

Nitrogênio de nitrito do rio

Nitrogênio de nitrato do rio

Fósforo orgânico do rio

Fósforo inorgânico do rio

Coliformes no rio (NMP/100ml)

Extensão do rio

Coeficiente de desoxigenação (K1)

Coeficiente de decomposição de DBO (Kd)

Coeficiente de sedimentação de DBO (Ks)

Coeficiente de reaeração (K2)

Coeficiente de sedimentação de nitrogênio orgânico (Kso)

Coeficiente de conversão de nitrogênio orgânico-amoniacal (Koa)

Coeficiente de conversão de nitrogênio amoniacal-nitrito (Kan)

Coeficiente de conversão de nitrito-nitrato (Knn)

Liberação de nitrogênio amoniacal pelo sedimento de fundo (Snamon)

Oxigênio para conversão amônia-nitrito (O2Namon)

Oxigênio para conversão nitrito-nitrato (O2Nnitri)

Coeficiente de inibição da nitrificação por baixo oxigênio dissolvido (knitr)

pH

Coeficiente de sedimentação do fósforo orgânico (Kspo)

Coeficiente de conversão fósforo orgânico-inorgânico (Koi)

Taxa de decaimento dos coliformes (Kb)

Comprimento da base menor do rio

Inclinação do rio

Coeficiente de Manning

Fator de secção (z) do rio

Dados de saída

Vazão

Área da seção transversal

Profundidade

DBO5 a jusante

Fósforo orgânico a jusante

Fósforo inorgânico a jusante

Nitrato a jusante

Nitrito a jusante

Nitrogênio amoniacal a jusante

Nitrogênio orgânico a jusante

Oxigênio dissolvido a jusante

Coliformes a jusante (NMP/100ml)

Curvas

Vazão

Modelo QUAL-UFMG

Descrição

O modelo QUAL-UFMG é uma ferramenta abrangente para simulação da qualidade da água em rios, desenvolvido para avaliar a distribuição espacial de poluentes e parâmetros de qualidade da água ao longo de cursos d'água. O modelo considera múltiplos processos biogeoquímicos e permite análise detalhada da evolução de diversos constituintes da qualidade da água.

 

O modelo oferece 3 configurações geométricas de seção transversal:

 

  • Seção trapezoidal: Para canais naturais com taludes laterais
  • Seção retangular: Para canais artificiais com paredes verticais
  • Seção triangular: Para canais em forma de V

 

Processos modelados

O modelo QUAL-UFMG simula os seguintes processos:

  • Ciclo do oxigênio: Consumo por decomposição da matéria orgânica, nitrificação e reaeração atmosférica
  • Ciclo do nitrogênio: Conversão de nitrogênio orgânico → amoniacal → nitrito → nitrato
  • Ciclo do fósforo: Transformação entre fósforo orgânico e inorgânico
  • Decaimento de coliformes: Morte natural de organismos indicadores
  • Sedimentação: Remoção de material particulado por deposição

Equações fundamentais

O modelo utiliza as seguintes equações diferenciais para descrever a variação temporal das concentrações:

Oxigênio dissolvido:

dODdt=K2(ODsatOD)KdDBOKTO2NamonNNH3KanfnitrO2NnitriNNO2Knnfnitr\frac{dOD}{dt} = K_2(OD_{sat} - OD) - K_d \cdot DBO \cdot K_T - O2Namon \cdot N_{NH_3} \cdot K_{an} \cdot f_{nitr} - O2Nnitri \cdot N_{NO_2} \cdot K_{nn} \cdot f_{nitr}

DBO:

dDBOdt=KrDBO\frac{dDBO}{dt} = -K_r \cdot DBO

onde Kr=Kd+KsK_r = K_d + K_s

Nitrogênio orgânico:

dNorgdt=(Kso+Koa)Norg\frac{dN_{org}}{dt} = -(K_{so} + K_{oa}) \cdot N_{org}

Nitrogênio amoniacal:

dNNH3dt=KoaNorgKanNNH3fnitr+SNH3h\frac{dN_{NH_3}}{dt} = K_{oa} \cdot N_{org} - K_{an} \cdot N_{NH_3} \cdot f_{nitr} + \frac{S_{NH_3}}{h}

Fósforo orgânico:

dPorgdt=(Kspo+Koi)Porg\frac{dP_{org}}{dt} = -(K_{spo} + K_{oi}) \cdot P_{org}

Correção de temperatura

Todos os coeficientes cinéticos são corrigidos para a temperatura usando:

KT=K20θ(T20)K_T = K_{20} \cdot \theta^{(T-20)}

onde θ\theta varia entre 1,024 e 1,080 dependendo do processo.

Hidráulica do canal

Para seções trapezoidais:

  • Área molhada: A=h(b+zh)A = h(b + zh)
  • Perímetro molhado: P=b+2h1+z2P = b + 2h\sqrt{1 + z^2}

Para seções retangulares:

  • Área molhada: A=bhA = bh
  • Perímetro molhado: P=b+2hP = b + 2h

Para seções triangulares:

  • Área molhada: A=zh2A = zh^2
  • Perímetro molhado: P=2h1+z2P = 2h\sqrt{1 + z^2}

A profundidade é calculada usando a equação de Manning:

Q=AnRh2/3S1/2Q = \frac{A}{n} \cdot R_h^{2/3} \cdot S^{1/2}

Dados de entrada

Parâmetro
Unidades padrão
Descrição
Temperatura
°C
Temperatura da água para correção dos coeficientes cinéticos
Vazão do rio
m³/s
Vazão volumétrica do curso d'água
Oxigênio dissolvido de saturação
mg/L
Concentração de saturação do OD na temperatura especificada
Oxigênio dissolvido do rio
mg/L
Concentração inicial de OD no rio
DBO5 do rio
mg/L
Demanda bioquímica de oxigênio em 5 dias inicial
Nitrogênio orgânico do rio
mg/L
Concentração inicial de nitrogênio orgânico
Nitrogênio amoniacal do rio
mg/L
Concentração inicial de nitrogênio amoniacal
Nitrogênio de nitrito do rio
mg/L
Concentração inicial de nitrito
Nitrogênio de nitrato do rio
mg/L
Concentração inicial de nitrato
Fósforo orgânico do rio
mg/L
Concentração inicial de fósforo orgânico
Fósforo inorgânico do rio
mg/L
Concentração inicial de fósforo inorgânico
Coliformes no rio
NMP/100ml
Número mais provável de coliformes por 100 ml
Extensão do rio
km
Comprimento do trecho de rio a ser simulado
Coeficiente de desoxigenação (K1)
d⁻¹
Taxa de remoção de OD por decomposição da DBO
Coeficiente de decomposição de DBO (Kd)
d⁻¹
Taxa de degradação da matéria orgânica
Coeficiente de sedimentação de DBO (Ks)
d⁻¹
Taxa de remoção de DBO por sedimentação
Coeficiente de reaeração (K2)
d⁻¹
Taxa de reoxigenação atmosférica
Coeficiente de sedimentação de nitrogênio orgânico (Kso)
d⁻¹
Taxa de sedimentação do nitrogênio orgânico
Coeficiente de conversão de nitrogênio orgânico-amoniacal (Koa)
d⁻¹
Taxa de conversão Norg → NH₃
Coeficiente de conversão de nitrogênio amoniacal-nitrito (Kan)
d⁻¹
Taxa de nitrificação NH₃ → NO₂
Coeficiente de conversão de nitrito-nitrato (Knn)
d⁻¹
Taxa de oxidação NO₂ → NO₃
Liberação de nitrogênio amoniacal pelo sedimento de fundo (Snamon)
g/m²/d
Fluxo de NH₃ do sedimento para a coluna d'água
Oxigênio para conversão amônia-nitrito (O2Namon)
-
Relação estequiométrica O₂:NH₃ na nitrificação
Oxigênio para conversão nitrito-nitrato (O2Nnitri)
-
Relação estequiométrica O₂:NO₂ na nitratação
Coeficiente de inibição da nitrificação por baixo oxigênio dissolvido (knitr)
d⁻¹
Parâmetro de inibição da nitrificação em baixas concentrações de OD
pH
-
Potencial hidrogeniônico da água
Coeficiente de sedimentação do fósforo orgânico (Kspo)
d⁻¹
Taxa de sedimentação do fósforo orgânico
Coeficiente de conversão fósforo orgânico-inorgânico (Koi)
d⁻¹
Taxa de mineralização Porg → Pinorg
Taxa de decaimento dos coliformes (Kb)
d⁻¹
Taxa de morte natural dos coliformes
Comprimento da base menor do rio
m
Largura do fundo do canal (seções trapezoidal/retangular)
Inclinação do rio
-
Declividade longitudinal do canal
Coeficiente de Manning
-
Coeficiente de rugosidade para cálculo hidráulico
Fator de secção (z) do rio
-
Inclinação das margens em seções trapezoidais/triangulares

Dados de saída

Parâmetro
Unidades padrão
Descrição
Vazão
m³/s
Vazão volumétrica confirmada
Área da seção transversal
Área molhada da seção transversal
Profundidade
m
Profundidade hidráulica calculada
DBO5 a jusante
kg/m³
Concentração final de DBO5 no final do trecho
Fósforo orgânico a jusante
mg/L
Concentração final de fósforo orgânico
Fósforo inorgânico a jusante
mg/L
Concentração final de fósforo inorgânico
Nitrato a jusante
mg/L
Concentração final de nitrato
Nitrito a jusante
mg/L
Concentração final de nitrito
Nitrogênio amoniacal a jusante
mg/L
Concentração final de nitrogênio amoniacal
Nitrogênio orgânico a jusante
mg/L
Concentração final de nitrogênio orgânico
Oxigênio dissolvido a jusante
mg/L
Concentração final de oxigênio dissolvido
Coliformes a jusante
NMP/100ml
Concentração final de coliformes
Curvas
-
Gráficos da evolução espacial de todos os parâmetros

Referências

  • Von Sperling, M. (2007 e 2014). Estudos e modelagem da qualidade da água de rios. DESA-UFMG e Editora UFMG.