Analise de Sensibilidade e Calibração

[brunodinato]

Olá, estou tentando fazer a Analise de sensibilidade. Pelo ArcSwat, eu tenho o manual que esta em anexo, porém quando coloco a tabela com dados observados, ele não roda a Analise de sensibildade, se alguem tiver uma tabela pronta para eu poder verificar a minha, ficarei muito agradecido, tenho também o SWAT-CUP e seu manual, porém não sei utiliza-lo, se alguém tiver um passo-a-passo ou para o SWAT-CUP ou para ArcSwat, agradeço muito.



Bruno I. Dinato

[Wellington Nunes]

Eu utilizei esse arquivo .Dat e a minha análise de sensibilidade rodou beleza. No caso os observados são de vazão diária. PS Não estou conseguindo anexar o arquivo de exemplo aqui... deixo aqui meu contato caso você ou mais alguém queira receber o arquivo. welling...@gmail.com

[Camila Guimaraes]

Boa noite!

Wellington, teria como me enviar o seu modelo do arquivo dos dados observados fazendo favor?

Obrigada!

[José Monteiro]

Olá,

Eu aconselho o uso de SWAT-CUP, pois além da análise de sensibilidade, tb é possível fazer a calibração e análise de incerteza. Assumindo que vcs precisem calibrar apenas para vazão em uma única estação, aqui vão umas dicas:

Os arquivos que devem ser modificados depois do projeto ser criado são
observed.txt
observed_rch.txt
var_file_name.txt
var_file_rch.txt
(todos na pasta SUFI2.IN)
SUFI2_extract_rch.def

Tanto em observed.txt como observed_rch.txt as observações precisam ser organizadas em 3 colunas, na primeira uma data genérica começando no primeiro mes (ou dia, no caso de simulações diárias) de simulação. A segunda coluna é uma referencia a data gregoriana, normalmente eu abrevio a variavel, estaçao e data (e.g. q_18_1998-01-00 para vazao na sub-bacia 18, janeiro de 1998), e na terceira o valor de vazao para essa data. Mais exemplos sobre esses arquivos podem ser encontrados dentro da pasta do programa SWAT-CUP em exemplos de modelos.

Em var_file_name e var_file_rch deve constar o nome que vc quer que os arquivo gerado receba (para esta estação). Por exemplo: q_18.txt

Segue um exemplo de SUFI2_extract_rch.def e comentários (marcados com #)

output.rch     : swat output file name
1              : number of variables to get # 1 pq calibramos apenas para uma variável, ou seja, vazão
7              : variable column number(s) in the swat output file (as many as the above number) # a coluna da vazão é 7

39              : total number of subbasins in the project # 39 sub-bacias no projeto

1              : number of subbasins to get for the first variable # extraindo para apenas uma sub-bacia
18     : subbasin numbers for the first variable # sub-bacia onde se localiza a estação com os dados observados

1973           : beginning year of simulation not including the warm up period
2006           : end year of simulation # inicio e fim da simulação, descontando-se os anos de warm up

2              : time step (1=daily, 2=monthly, 3=yearly) # 2 para simulações mensais

Se tiverem algum problema, podem escrever.

Um abraço e boa sorte,
Zé

-- 
José A. F. Monteiro (PhD)
Visiting Researcher at Applied Limnology Laboratory, UFSJ, Brazil
Visiting Researcher at Spatial Sciences Laboratory, Texas A&M University, USA

E-Mail: jose.m...@gmx.ch
ResearchGate: https://www.researchgate.net/profile/Jose_Alberto_Monteiro

[Rafael Almeida]

Bom dia Jose,

Estou entrando em contato pois vi que voce sempre responde e estou em apuros...

Nao estou conseguindo rodar o SWAT-CUP, mesmo apos seguir o manual e as suas explicacoes...

Ele roda e no final fala que nao existe arquivo 95 ppu-plot e nenhum outro na aba Calibration Outputs - segue o print com a mensagem que ele da apos rodar, e alguns outros com informacoes que possam ajudar...

Algumas informacoes importantes: Estou simulando apenas vazoes para uma bacia no leste de MG, sendo que o periodo escolhido eh de 01-09-2006 ate 31-12-2014.... 2 anos de warm-up (2006 e 2007), 3 anos para calibracao (2008, 2009 e 2010) e os demais para a validacao (2011, 2012, 2013 e 2014). tenho dados de 6 estacoes, porem quero calibrar para uma e validar nas demais...

Desculpe a falta de acentos, mas o meu teclado esta desconfigurado...

 

Espero que possa me ajudar...
Atenciosamente,

Rafael Alvarenga Almeida

[José Monteiro]

Olá Rafael,

Só com os sceenshots eu não posso avaliar o que está acontecendo, mas duas coisas me chamaram a atenção. Seu threshold para a função objetivo é muito alto, seria melhor vc usar um valor mais baixo, tipo 0.3 (não tenho muita experiência em usar NS, geralmente eu uso bR2=0.3).

Pelo número de pontos de dados eu assumo que vc esteja fazendo calibração diária. Eu recomendaria vc começar pela calibração mensal para ajustar o balanço hídrico. Além disso, a calibração mensal provavelmente seria suficiente para um projeto com o objetivo de verificar o balanço hídrico com sazolalidade, independentemente de eventos específicos de precipitação.

Se vc não conseguir encontrar o problema, me mande em privado as pastas SUFI2.IN, SUFI2.OUT e os arquivos SUFI2_extract_rch.def e file.cio.

Um abraço,
José

[Rafael Almeida]

Bom dia José, 

Primeiramente muito obrigado pela atenção!

Já sai deste problema e cai em outro... hahaha

A função objetivo não interfere no caso, pois ela nunca se ajusta, pois no meu novo problema, a calibração do SWAT-Cup me retorna uma linha reta sempre com o mesmo valor correspondente a soma das vazões no período.

Estou usando uma calibração diária, pois os objetivos do meu trabalho estão relacionados a esta. Não entendi como uma simulação em base mensal pode ajudar, se realizada primeiro...

Tem alguma ideia do que pode estar acontecendo para o SWAT-Cup retornar uma simulação inteiramente linear?

[Anaís Ostroski]

Olá Rafael, 

Eu obtive o mesmo erro, ele roda mas os arquivos do Calibration Output são inexistentes.

Pode me falar como arrumou?

Agradeço

[José Monteiro]

Olá Rafael,

Eu não disse que ajustando a fnção objetivo resolveria o problema, mas o valor dela está alto demais. A calibração mensal é mais simples pois o modelo roda mais rápido e vc pode ter resultados prévios mais rápidos e re-ajustar os parâmetros. Além disso, a variabilidade temporal e espacial da precipitação, que tem alta influência nos dados diários, tem efeito bem menor nos dados mensais: https://www.researchgate.net/publication/283014354_Accuracy_of_grid_precipitation_data_for_Brazil_Application_in_river_discharge_modelling_of_the_Tocantins_catchment

Da uma olhada se vc não mandou imprimir os dados de vazão em escala logarítmica (print log flow).

Um abraço,
José

[Rafael Almeida]

Bom dia José,

Não marquei a opção print log flow...
Estou bem perdido com os resultados, ou a falta destes... =(

[José Monteiro]

Outra opçõ que pode gerar uma diferença enorme entre observação e simulação é um posicionamento incorreto da estação. Se a estação se localiza próxima a uma confluência, verifica se ela esta corretamente na sub-bacia do tributário ou do canal principal. O posicionamento automático é apenas uma sugestão.

[Rafael Almeida]

Boa tarde José...

Também não é o caso, pois ela está situada corretamente, porém apareceu uma coisa nova e interessante...

Tive observando que quando abro o arquivo de saída do SWAT (output) - observo que a coluna que está a vazão(FLOW_OUTcms) é a 8 e não a 7...
Por coincidência variei a coluna no CUP, e ele apresentou um erro. Porém gerou um arquivo de saída a qual existem simulações...

Agora fiquei mais confuso ainda...

Saberia ajudar?

[Rafael Almeida]

Seguem os prints em anexo, e outro fato interesante que acabei de observar, a constante que o CUP retorna corresponde a coluna 6 do arquivo de saída do SWAT, que é a AREAkm2...

[José Monteiro]

Olá Rafael,

A coluna para vazão é a 7. Que coluna vc está usando?

[Rafael Almeida]

tive pesquisando e parece que nem sempre é a 7, instruções de Karim Abbaspour no fórum específico do SWAT CUP... no meu caso é a 9, olhei no arquivo do SWAT (output.rch)...

Agora acontece um novo problema, começa a corrida e essa nunca termina... já deixei por mais de 25 horas rodando, e ele nunca passa para a simulação 2... Mais uma vez me encontro perdido...

[Magda Melo]

Eu não tenho dados observados de vazão, como faço para observar o fósforo. Alguém sabe me dizer, pelo que vi, tenho que mudar a coluna, mas onde identifico a coluna certa e o que mudar na programação?

Muito obrigada.

Em quinta-feira, 24 de setembro de 2015 11:16:17 UTC-3, José Monteiro escreveu:

[Magda Melo]

Bruno eu segui esse tutorial, o meu trava na primeira simulação.

SWAT-CUP | SWAT Model Calibration | SUFI-2 Algoritham | Auto-Calibration

 agora esse video ensina a montar a tabela

Preparing Data Input For SWAT-CUP | Auto-Calibration

todos no you tube

[viniciusraposo]

Olá, Magda. Não sou nenhum especialista, somente alguém que já pelejou muito no SWAT-CUP. Mas assim, há uma recomendação no manual do SWAT-CUP para primeiramente calibrar a vazão: "se você deseja calibrar usando, por exemplo, vazão, nitrato e fosfato, vocês deve calibrar primeiro a vazão, por ser a principal variável de controle". Você pode enviar um print do erro? Muitas vezes é em alguma configuração. Há um outro grupo no Google chamado SWAT-CUP, onde, inclusive, participa o Dr. Abbaspour, o qual também costuma ajudar os demais usuários. Abraços.

[washington luiz]

Estou com esse problema, sou iniciante e estou tentando resolver os problemas, Pelo que entendi, não gera o SUFI.Out, mas organizo todo os dados, reviso e mesmo assim o erro persiste.

Alguém poderia me dar uma ajuda?