Volume armazenado no reservatório subterrâneo

Prezados,

Gostaria de saber se é possível acessar dados dos volumes armazenados no reservatório subterrâneo? de acordo com na fórmula abaixo:


Se este não for possível, o escoamento subterrâneo (DBas) o é? ou a variação do armazenamento subterrâneo ("Delta" S)?

Atenciosamente,

Laíne

Olá Laíne, tudo bem?

Complementando a resposta da sua pergunta anterior: https://mgb-iph.forumeiros.com/t206-mgb-subterraneo
Existe um arquivo de saída na pasta Output chamado "QBASTUDO_bin.MGB" que contém para todas as mini bacias e todos os intervalos de tempo os valores da vazão de saída do reservatório subterrâneo "QBAS".
Trata-se de um arquivo binário, você pode ler esses dados em um programa como MATLAB ou OCTAVE usando o script a seguir:
OBS: preencher as informações de mini bacias e intervalos de tempo de acordo com os valores disponíveis no arquivo LOG_MGB

clear all

% Insira o número de minibacias
nMini = 0;

% Insira o número de intervalos de tempo
NT = 0;

% Abrir arquivo de vazão QBASTUDO_bin
fid = fopen('QBASTUDO_bin.MGB','r');

% Gravar dados do arquivo QBAS em uma variável matricial
QBAS = fread(fid,[nMini,NT],'single');

% Fechar arquivo
fclose(fid);


Atenciosamente,
Gustavo.

Prezado Gustavo,

Muito obrigada! Gostaria de saber qual a definição desses arquivos de saída selecionados?

Olá Laíne, seguem abaixo algumas definições dos arquivos encontrados na pasta output do MGB:

YTUDO.MGB: arquivo binário que contém para cada minibacia e para cada intervalo de tempo da simulação o nível de água do trecho de rio em metros.
NOSOLO: esse arquivo já se encontra desativado nas versões mais recentes do MGB (tratava-se de um arquivo de teste do balanço de água no solo para URHs específicas em estudo). Logo, sugiro que desconsidere os resultados encontrados nesse arquivo de saída.
WATERDEPTH.QCL: arquivo com os níveis de água para cada intervalo de tempo para as minibacias selecionadas com postos fluviométricos que foram utilizados na calibração dos hidrogramas simulados. De forma análoga, o arquivo VAZAO.QCL informa a vazão observada nos postos fluviométricos selecionados. Caso você tenha realizado uma simulação inercial, você pode observar que os valores encontrados no arquivo WATERDEPTH.QCL são os mesmos descritos nos arquivos SIM_INERC_Hfl_x para cada minibacia "x".
QPROP.PRP: também é um arquivo de testes mas que ainda está presente nas versões atuais do MGB, esse arquivo serve para informar a origem do escoamento em uma minibacia (nesse caso, o número da minibacia escolhida deve ser alterado diretamente no código fonte do MGB). O arquivo separa o escoamento de base, sub-superficial ou superficial fazendo um acumulado das vazões de todas as minibacias que estão a montante da minibacia desejada (OBS: esse arquivo só apresenta resultados utilizando o método de propagação de vazões de Muskingum Cunge).

Atenciosamente,
Gustavo.

Prezado Gustavo,

Há possibilidade de ler o arquivo QPROP.PRP no matlab? e se possível, consigo ler os dados de escoamento de base, sub-superficial ou superficial de todas? E qual a diferença dele com o QITUDO.binMGB? o arquivo QITUDO informa a vazão incremental, certo? a vazão incremental não é a vazão originária na minibacia (sem acumulações)?

Outra dúvida, mas agora é a respeito da fórmula MUSKINGUM-CUNGE utilizada no modelo. O tempo t+1 e t, da fórmula abaixo, é diário? ou seja a vazão no dia 1 e no dia 2?

Muito Obrigada

Olá Laíne,

Você pode ler o arquivo QPROP.PRP no matlab sim, mas isso não é necessário pois ele não se encontra em formato binário. Você pode acessar seus dados diretamente abrindo-o com um editor de texto qualquer. Nesse arquivo, estão separadas em colunas para cada intervalo de tempo as vazões acumuladas do escoamento superficial, sub-superficial e de base para a mini bacia configurada no código fonte.
Em relação ao arquivo QITUDO.binMGB sua definição está certa, trata-se da vazão incremental gerada em cada mini bacia do modelo.

Sobre a equação de Muskingum-Cunge, o intervalo de tempo da fórmula não é o mesmo intervalo de tempo dos dados de vazão gerados pelo modelo. O intervalo de tempo equação de MC é calculado internamente e apresenta valores diferentes para cada mini bacia gerada.

Atenciosamente,
Gustavo.

Olá Gustavo,

Abri o arquivo QPROP.PRP (rodei o MGB versão mapwindow 2017) em .txt e apresentou como na figura abaixo. Não entendi o valor grifado em vermelho e amarelo, o verde são os dias. Pensei que havia colunas separadas, para cada intervalo de tempo, as vazões acumuladas do escoamento superficial, sub-superficial e base, em todas mini bacias.

Oi Laíne,

O arquivo QPROP.PRP não apresenta dados do escoamento para todas as minibacias, ele apresenta apenas para uma minibacia específica de estudo que você pode selecionar alterando a numeração diretamente no código fonte do MGB.
Sobre a separação das colunas do arquivo, pode ser algo relacionado a versão do MGB que você está usando. Tente usar a versão mais atual disponível no site com o QGIS, acredito que assim o arquivo deverá ficar com colunas separadas indicando os intervalos de tempo, o escoamento superficial, sub-superficial e de base (nesta ordem) igual a figura abaixo:



Atenciosamente,

Gustavo.

Boa Tarde Gustavo,

Eu tentei rodar na última versão mas está dando problemas ao poligonizar o arquivo catchment e também a ferramenta de uso de climas diários. Assim, estou tentando usar a versão MGB 4.2 QGIS 2.18, mas ao interpolar dados de chuva está dando essa mensagem abaixo. O que há de errado?




Consulte o final desta mensagem para obter detalhes sobre como chamar a
depuração just-in-time (JIT) em vez desta caixa de diálogo.

************** Texto de Exceção **************
System.UnauthorizedAccessException: O acesso ao caminho 'C:\Program Files\QGIS 2.18\apps\qgis-ltr\python\plugins\MGB\bin\proj.dirmgb' foi negado.
em Microsoft.VisualBasic.FileSystem.FileOpen(Int32 FileNumber, String FileName, OpenMode Mode, OpenAccess Access, OpenShare Share, Int32 RecordLength)
em Form_plu.Form_plu.Form_plu_Load(Object sender, EventArgs e)
em System.EventHandler.Invoke(Object sender, EventArgs e)
em System.Windows.Forms.Form.OnLoad(EventArgs e)
em System.Windows.Forms.Control.CreateControl(Boolean fIgnoreVisible)
em System.Windows.Forms.Control.CreateControl()
em System.Windows.Forms.Control.WmShowWindow(Message& m)
em System.Windows.Forms.Control.WndProc(Message& m)
em System.Windows.Forms.Control.ControlNativeWindow.WndProc(Message& m)
em System.Windows.Forms.NativeWindow.Callback(IntPtr hWnd, Int32 msg, IntPtr wparam, IntPtr lparam)


************** Conjuntos de Módulos (Assemblies) Carregados **************
mscorlib
Versão do Conjunto de Módulos (Assembly): 2.0.0.0
Versão do Win32: 2.0.50727.8806 (QFE.050727-8800)
Base de Código: file:///C:/Windows/Microsoft.NET/Framework64/v2.0.50727/mscorlib.dll
----------------------------------------
Form_plu
Versão do Conjunto de Módulos (Assembly): 1.0.0.0
Versão do Win32: 1.0.0.0
Base de Código: file:///C:/Program%20Files/QGIS%202.18/apps/qgis-ltr/python/plugins/MGB/bin/Form_plu.exe
----------------------------------------
Microsoft.VisualBasic
Versão do Conjunto de Módulos (Assembly): 8.0.0.0
Versão do Win32: 8.0.50727.5483 (Win7SP1GDR.050727-5400)
Base de Código: file:///C:/Windows/assembly/GAC_MSIL/Microsoft.VisualBasic/8.0.0.0__b03f5f7f11d50a3a/Microsoft.VisualBasic.dll
----------------------------------------
System
Versão do Conjunto de Módulos (Assembly): 2.0.0.0
Versão do Win32: 2.0.50727.8806 (QFE.050727-8800)
Base de Código: file:///C:/Windows/assembly/GAC_MSIL/System/2.0.0.0__b77a5c561934e089/System.dll
----------------------------------------
System.Windows.Forms
Versão do Conjunto de Módulos (Assembly): 2.0.0.0
Versão do Win32: 2.0.50727.5491 (Win7SP1GDR.050727-5400)
Base de Código: file:///C:/Windows/assembly/GAC_MSIL/System.Windows.Forms/2.0.0.0__b77a5c561934e089/System.Windows.Forms.dll
----------------------------------------
System.Drawing
Versão do Conjunto de Módulos (Assembly): 2.0.0.0
Versão do Win32: 2.0.50727.5495 (Win7SP1GDR.050727-5400)
Base de Código: file:///C:/Windows/assembly/GAC_MSIL/System.Drawing/2.0.0.0__b03f5f7f11d50a3a/System.Drawing.dll
----------------------------------------
System.Runtime.Remoting
Versão do Conjunto de Módulos (Assembly): 2.0.0.0
Versão do Win32: 2.0.50727.8771 (QFE.050727-8700)
Base de Código: file:///C:/Windows/assembly/GAC_MSIL/System.Runtime.Remoting/2.0.0.0__b77a5c561934e089/System.Runtime.Remoting.dll
----------------------------------------
System.Data
Versão do Conjunto de Módulos (Assembly): 2.0.0.0
Versão do Win32: 2.0.50727.8762 (QFE.050727-8700)
Base de Código: file:///C:/Windows/assembly/GAC_64/System.Data/2.0.0.0__b77a5c561934e089/System.Data.dll
----------------------------------------
System.Xml
Versão do Conjunto de Módulos (Assembly): 2.0.0.0
Versão do Win32: 2.0.50727.8773 (QFE.050727-8700)
Base de Código: file:///C:/Windows/assembly/GAC_MSIL/System.Xml/2.0.0.0__b77a5c561934e089/System.Xml.dll
----------------------------------------
System.Windows.Forms.resources
Versão do Conjunto de Módulos (Assembly): 2.0.0.0
Versão do Win32: 2.0.50727.5420 (Win7SP1.050727-5400)
Base de Código: file:///C:/Windows/assembly/GAC_MSIL/System.Windows.Forms.resources/2.0.0.0_pt-BR_b77a5c561934e089/System.Windows.Forms.resources.dll
----------------------------------------
mscorlib.resources
Versão do Conjunto de Módulos (Assembly): 2.0.0.0
Versão do Win32: 2.0.50727.8806 (QFE.050727-8800)
Base de Código: file:///C:/Windows/Microsoft.NET/Framework64/v2.0.50727/mscorlib.dll
----------------------------------------

************** Depuração JIT **************
Para habilitar a depuração just-in-time (JIT), o arquivo .config deste
aplicativo ou computador (machine.config) deve ter o valor
jitDebugging definido na seção system.windows.forms.
O aplicativo também deve ser compilado com a depuração
habilitada.

Por exemplo:

<configuration>
<system.windows.forms jitDebugging="true" />
</configuration>

Quando a depuração JIT está habilitada, qualquer exceção não tratada
será enviada ao depurador JIT registrado no computador,
em vez de ser tratada nesta caixa de diálogo.

Oi Laíne,

Apenas para confirmar, você encontrou um erro ao tentar poligonizar o arquivo de mini bacias utilizando qual versão do QGIS?
É importante que seja utilizada a versão (3.4.14) do QGIS com a última versão do MGB (4.5). Porém, eu pessoalmente já testei a versão 3.10 do QGIS com a última versão do MGB e não tive problemas para rodar nenhuma ferramenta.

Quanto ao problema com a interpolação de chuva, você deve garantir que no período escolhido da simulação não haja nenhum ano em que não exista ao menos uma estação com dados completos de chuva. Você pode fazer essa análise clicando na opção "time availability". Acredito que possa ser esse o erro que você está encontrando.

Atenciosamente,
Gustavo.

Olá Gustavo,

Eu tentei a última versão do MGB com o qgis 3.4.14 e não consegui poligonizar, já com a versão 3.10 consegui, entretanto na simulação dá mensagem de erro e só tem 3 arquivos no output. Quanto a interpolação dos climas diários, a ferramenta nem abre, tanto no qgis 3.4.14 quanto na versão 3.10.

Gustavo Gabbardo dos Reis escreveu:Oi Laíne,

Apenas para confirmar, você encontrou um erro ao tentar poligonizar o arquivo de mini bacias utilizando qual versão do QGIS?
É importante que seja utilizada a versão (3.4.14) do QGIS com a última versão do MGB (4.5). Porém, eu pessoalmente já testei a versão 3.10 do QGIS com a última versão do MGB e não tive problemas para rodar nenhuma ferramenta.

Quanto ao problema com a interpolação de chuva, você deve garantir que no período escolhido da simulação não haja nenhum ano em que não exista ao menos uma estação com dados completos de chuva. Você pode fazer essa análise clicando na opção "time availability". Acredito que possa ser esse o erro que você está encontrando.

Atenciosamente,
Gustavo.

Prezado Gustavo,

Consegui rodar o MGB 4.5 no Qgis 3.10.9, no entanto o problema no arquivo QPROP e a ferramenta Using Daily data continua. Vou relatar abaixo alguns problemas que tive.

1) Conferi se a foz estava formada corretamente, comparando com as cartas topográficas do IBGE. Há possibilidade de união/cruzamento de linhas de drenagem, e a minha não estava certa, então utilizei a ferramenta Stream-Burning do ArcHydro, inserido no Arcgis. No entanto ao usar o Qgis tentei utilizar a ferramenta QGIS ASCII-IPHHT não deu certo, mas consegui utilizar o arquivo sem precisar utilizar essa ferramenta.

2) Ao rodar o MGB com intuito de criar o arquivo watershed verifiquei como estava o valor min e máx. Cliquei com o botão direito em cima do arquivo watershed e fui em propriedades, na parte de simbologia, mais especificamente em renderização da banda, verifiquei se o min=1 e máx=1, e não estava, então corrigi.

3) Ao criar o arquivo catchment e urh, também verifiquei se ele está iniciando a partir do número 1, e não estava, então fiz a mesma correção do item anterior.

4) Quando fui realizar a interpolação de chuvas com os dados que eu já tinha (formatado com o super manejo de dados) percebi que KB estava pequeno, então utilizei a ferramenta Ana data acquisition e deu certo.

5) Ao interpolar a vazão não estava dando certo, percebi que os valores ID das minibacias do discharge (utilizei a ferramenta automatically sugest catchment e update unit-catchments polygon) e run não coincidiram. Então comparei os dados dos postos fluviométricos da Hidroweb/DAEE-SP com os da internal database, e não estavam coincidindo, assim, corrigi-os e deu certo.

6) Também percebi que apesar da simulação ter ocorrido normalmente, a minha pasta do input não há o arquivo dos parâmetros fixos.

Atenciosamente,

Laíne

Ola Laíne, tudo bem?

- O IPH HydroTools aceita os formatos ASCII gerados pelo QGIS, então realmente não seria necessário utilizar a ferramenta de conversão no QGIS nesse caso.
- É estranho que você tenha feito uma simulação sem carregar o arquivo de parâmetros fixos do solo. Como é um arquivo gerado em uma etapa anterior, acredito que ele não esteja na pasta de input pois não deve ter sido devidamente carregado na etapa posterior com a ferramenta "Create/Edit Simulation", não tenho certeza se é esse o caso.
- Como no seu caso você conseguiu realizar sem problemas o pré-processamento dos dados, eu vou te enviar o link de um vídeo no qual está gravado apenas alguns detalhes de como rodar o modelo MGB. Tente identificar se você está fazendo todos os passos da mesma forma que é feito no vídeo e se não seria uma alternativa para você a utilização de dados de normais climatológicas para o Brasil (INMET) ao invés de utilizar uma base de dados diária, que aparentemente está dando erro na sua simulação.

https://www.youtube.com/watch?v=KrS-eSE4MLQ

Atenciosamente,
Gustavo.

Prezado Gustavo,

Eu já havia assistido o vídeo. Acredito fortemente que o problema esta ocorrendo porque eu fiz a correção com Stream-Burning do ArcHydro. Depois que te mandei uma mensagem fui rodar o MGB para um posto fluviométrico perto da foz (minha bacia tem 2 postos com dados em períodos diferentes), e percebi que de novo meu catchment não começava a partir de 1. Não conseguia interpolar chuvas e não consegui simular. Justamente porque a área de drenagem do posto contém a minibacia 1, onde estava o erro das drenagens para a formação da foz. Se você quiser te mando arquivos do MDE sem correção e o MDE com correção do Stream-Burning. Muito Obrigada!

Atenciosamente,

Laíne