RESUMO
Esta pesquisa objetivou projetar um sistema de controle adaptativo para o modelamento e controle, em tempo real, de um sistema de distribuição de água. Os controladores atuaram no controle da pressão por meio da variação da velocidade de rotação do sistema de bombeamento e do ângulo de abertura de válvulas de controle. O sistema de controle é do tipo Variância Mínima Generalizado (GMV) autoajustável, cujos parâmetros foram estimados em tempo real pelo método dos mínimos quadrados recursivo (MQR). Os ensaios realizados na bancada experimental apresentaram resultados satisfatórios, com um erro máximo no controle de pressão de 2,12% (0,42 mca) e uma redução de 19,2% no consumo de energia elétrica para o sistema atuando na condição mais desfavorável.
Palavras-chave:
sistema de bombeamento; eficiência energética; sistema de controle.
ABSTRACT
This research aimed to design an adaptive control system for modeling and real-time control of a water distribution system. Air traffic controllers acted on pressure control by varying the rotational speed of the pumping system and opening angle control valve. The system of control is of type Generalized Minimum Variance self-tuning, whose parameters were estimated in real time by the recursive least squares method. The tests performed on experimental workbench presented satisfactory results, with a maximum error on pressure control of 2.12% (0.42 mca) and a decrease of 19.2% in the consumption of electrical energy to the system in the most unfavorable condition.
Keywords:
pumping system; energy efficiency; control system.
INTRODUÇÃO
Água e energia são dois recursos essenciais e possuem uma histórica interdependência. Há uma forte sinergia existente entre os recursos hídricos e energéticos. Dentre os consumidores desses recursos, o setor de saneamento se destaca entre os usuários com maiores demandas. O setor de saneamento consome de 2 a 10% da energia total de um país (PELLI & HITZ, 2000PELLI, T. & HITZ, H.U. Energy indicators and savings in water supply. (2000) Journal AWWA , v. 92, n. 6, p. 55-62.), sendo que a maior parte é consumida pelos sistemas de bombeamento (WANG et al., 2009WANGA, J.Y.; CHANG, T.P.; CHEN, J.S. (2009) An enhanced genetic algorithm for bi-objective pump scheduling in water supply., Expert Systems with Applications v. 36, p. 249-258.). Nos Estados Unidos, segundo Sanders e Webber (2012SANDERS, K.T. & WEBBER, M.E. (2012) Evaluating the energy consumed for water use in the United States. Environmental Research Letters , v. 7, p. 1-11.), o setor consumiu cerca de 8,4% da energia gerada em 2010. No Brasil, estima-se que o setor de saneamento consome aproximadamente 2,5% do total de energia elétrica (ELETROBRAS, 2008ELETROBRÁS - Centrais Elétricas Brasileiras S.A. (2008). Disponível em: <http://www.eletrobras.com>. 2008. Acesso em: 28/10/2008.
http://www.eletrobras.com...
).
O crescimento da população dos centros urbanos tem impulsionado o aumento da demanda por água e energia. Cada vez mais estações de bombeamento e boosters têm sido empregados em larga escala para manter o sistema de distribuição de água dentro de níveis de pressão e vazão compatíveis às necessidades. O uso eficiente de água e energia elétrica deve ser prioridade nas companhias. Globalmente, a energia está entre os três principais itens de custo para os prestadores de serviços de água e esgoto, e em muitos casos representa a segunda despesa, perdendo apenas para os custos com pessoal. Nos países em desenvolvimento, a energia é geralmente o custo mais elevado associado com a fonte de água.
Diante do exposto, torna-se primordial o desenvolvimento de metodologias e técnicas que proporcionem o aumento da eficiência energética dos sistemas de bombeamento e a redução das perdas de água. Esta pesquisa objetivou projetar um sistema de controle adaptativo para o controle, em tempo real, da pressão de sistemas de distribuição de água por meio da atuação na rotação do sistema de bombeamento e na abertura de válvulas.
MODELAGEM MATEMÁTICA DO SISTEMA
As redes de distribuição de água com injeção direta estão sujeitas a variação da demanda ao longo do tempo, fenômeno que ocasiona alterações significativas nas condições de operação dos sistemas de controle. Os métodos clássicos de controle usualmente empregados no setor de saneamento apresentam limitações quando aplicados a processos com características complexas (não linearidades, instabilidade, tempo morto elevado, comportamento não estacionário e acoplamento entre as variáveis, etc.). Desse modo, as técnicas adaptativas foram adotadas porque o tipo de sistema de distribuição de água estudado requer o controle simultâneo de pressão por meio de equipamentos com características e comportamentos diferentes. A literatura apresenta inúmeras referências que evidenciam que os algoritmos adaptativos são uma solução adequada para o controle de sistemas dinâmicos complexos (PURWAR et al., 2008PURWARA, S.; KARB, I.N.; JHAB, A.N. (2008) Adaptive output feedback tracking control of robot manipulators using position measurements only. Expert Systems with Applications , v. 34, n. 4, p. 2789-2798.; LIUZZO & TOMEI, 2008LIUZZO, S. & TOMEI, P.A. (2008) A global adaptive learning control for robotic manipulators. Automatica , v. 44, n. 5, p. 1379-1384.).
O princípio básico do funcionamento de um sistema de controle adaptativo é que os parâmetros dos controladores se alteram a cada mudança dos parâmetros da planta. Foram adotados nesta pesquisa controladores de Variância Mínima Generalizada (GMV) autoajustáveis, cujos parâmetros foram estimados em tempo real pelo algoritmo dos mínimos quadrados recursivo (MQR), visando à atualização do sistema de controle diante das não linearidades apresentadas na operação da rede hidráulica.
O estimador dos mínimos quadrados manipula as medidas de entrada/saída nas formas não iterativa e iterativa por meio de um algoritmo recursivo. Este considera, a priori , que a ordem do modelo é conhecida e que as amostras entrada/saída estão disponíveis a cada período de amostragem no universo da experimentação. Diversos modelos de entrada/saída estão disponíveis na literatura, sendo escolhido o modelo linear Autorregressivo com Entradas Exógenas - ARX (Autoregressive with External Input ) porque este apresenta a vantagem de identificar um sistema MIMO (múltiplas entradas e múltiplas saídas) tão facilmente como um sistema SISO (uma entrada e uma saída) (HUUSOMA et al., 2012HUUSOMA, J.K.; POULSENB, N.K.; JØRGENSENA, S.B.; JØRGENSENB, J.B. (2012) Tuning SISO offset-free Model Predictive Control based on ARX models. Journal of Process Control , v. 22, n. 10, p. 1997-2007.).
O modelo ARX inclui em sua estrutura o erro e por isso é normalmente classificado como pertencendo à classe de modelos de erro na equação. O "erro" é interpretado como ruído e engloba todos os efeitos que não forem explicados nos dados por parte do modelo, tais como não linearidades, imprecisão de medições, incertezas, presenças de perturbações. Para um sistema SISO e com uma perturbação, a equação do modelo ARX é expressa pela Equação 1. A melhor previsão da saída (valor estimado), ŷ(t), ocorre quando o erro E(t) for zero. Logo, o estimador foi obtido pela minimização do erro de previsão (Equação 2).
onde y(t) é a saída do sistema (observadas e/ou medidas); ŷ(t) são as saídas estimadas do modelo; φm(t) é o vetor de medidas; θ(t) é o vetor de parâmetros do sistema; E(t) é o erro de modelagem.
A minimização da Equação 2 resulta nas Equações 3 e 4, que representam o estimador dos mínimos quadrados não recursivo. A Equação 4 é a matriz de covariância, P(t), da qual depende a precisão das estimativas.
O algoritmo estimador dos parâmetros é iterativo, com o modelo do sistema sendo atualizado para cada período de amostragem. Dessa forma, o estimador dos mínimos quadrados não recursivo foi adaptado para resultar no algoritmo MQR, equação 5. No MQR, as estimativas anteriormente calculadas, ao instante de tempo t, são calculas para obtenção dos parâmetros atualizados θ(t) (CARVALHO, 2009CARVALHO, J.S.F. (2009) Controle adaptativo aplicado em um robô manipulador de dois graus de liberdade planar . Tese (Doutorado) - Universidade Federal da Paraíba, João Pessoa.).
A qualidade do modelo estimado depende, em parte, da natureza do sinal de entrada aplicado no sistema. Um sinal Pseudo Random Binary Sequence (PRBS) foi aplicado como entrada porque este excita os modos do sistema que correspondem ao conteúdo espectral, prevenindo a ocorrência de um condicionamento numérico ruim e a formação de colunas iguais na matriz de covariância.
Para processos variantes no tempo, é necessário fornecer ao MQR capacidade de adaptação mínima, o que foi considerado por meio da regulagem de um elemento chamado de fator λ, que deu maior importância às novas medidas. Com isso, o MQR assumiu a forma da Equação 6, onde k(t) é o ganho do estimador (Equação 7). O MQR com o fator λ evitou que os elementos da matriz de covariância tendessem para zero, mantendo o estimador em alerta para rastrear dinâmicas variantes.
O tipo de sistema de distribuição de água avaliado possui múltiplas entradas e múltiplas saídas. Reescrevendo a Equação 1 para um sistema MIMO, esta assume a forma dada pelos itens da Equação 8.
onde YM(t) é a saída do sistema (observadas e/ou medidas); U(t) são as variáveis de entrada; EM(t) é o erro de modelagem; Θ(t) é a matriz de parâmetros referente ao sistema MIMO; A e B são conjuntos de matrizes de polinômios com parâmetros ai e bj, respectivamente; Ψ(t) é o vetor de medidas de entradas e saídas referentes ao sistema MIMO.
O modelo matemático da bancada experimental foi identificado, em tempo real, a partir dos dados de excitação do inversor de frequência, u(t), e as respostas obtidas nos transdutores de pressão, y1(t) e y2(t). Definiu-se, a priori , que a estrutura do modelo dos ramais do sistema experimental tinha três polos, dois zeros e não possuía atraso, o que resultou nos itens da Equação 9 para o modelo, sendo ai e bj os parâmetros do sistema experimental de distribuição de água.
De forma análoga a um sistema SISO, a função custo do estimador MIMO, JMQ, foi estimada por meio da equação 10, resultando na matriz de parâmetros estimados Θ(t)(MQR) e no ganho do estimador k(t), ambos dados pelas Equações 11 e 12, respectivamente. Com a solução destas equações, obtêm-se os parâmetros do sistema experimental de distribuição de água ai e bj para cada período de amostragem. As respostas estimadas, ŶM(t), e os erros de previsão, ΞM(t), foram obtidos por meio das equações 13 e 14, respectivamente.
Para a identificação da bancada experimental em estudo, em tempo real, o MQR foi executado com o auxílio do programa computacional Matlab(r), com o tempo de amostragem de 5 ms. Nas excitações impulsivas aplicadas aos dois pontos de medição (um em cada ramal do sistema experimental) foram adicionados sinais PRBS com amplitudes da ordem de 0,5% das respectivas excitações.
CONTROLADOR DE VARIÂNCIA MÍNIMA GENERALIZADA AUTOAJUSTÁVEL
A técnica de controle GMV autoajustável foi desenvolvida na década de 1970 por Clarke e Gawthrop (1975CLARKE, D.W. & GAWTHROP, P.J. (1975) Self-Tuning controller. Proceedings of the Institution of Electrical Engineers , v. 122, n. 9, p. 929-934.), a partir da generalização do regulador de variância mínima proposto por Åström e Wittenmark (1973ÅSTRÖM, K.J. & WITTENMARK, B. (1973) On self-tuning regulators. Automatica , v. 9, n. 2, p. 185-199.). Baseia-se na minimização de uma função custo, a qual é sintetizada em função do processo e da dinâmica desejada para a malha de controle por meio de uma saída generalizada.
O sinal de controle é otimizado de maneira a determinar, pelos parâmetros de projeto, a dinâmica transitória, reduzindo a sobre-elevação e eliminando ou atenuando o erro em regime permanente. Neste caso, optou-se pela liberdade de incluir parâmetros para atender às especificações de projeto, conferindo uma maior flexibilidade à estrutura do controlador. Considerando tal conceito, a função custo a ser minimizada foi composta para possibilitar a ponderação dos sinais de controle, referência e saída do sistema.
O diagrama esquemático do sistema de controle adaptativo supervisionado é mostrado na Figura 1 e a equação que rege tal sistema é escrita de acordo com a Equação 15. Observa-se que esta equação é composta por informações disponíveis no instante k e por informações futuras referentes ao instante k+d. Sendo assim, consideraram-se a identidade polinomial sugerida por Marchi et al. (1998MARCHI, P.A.; ROCHA, C.R.; VENTURO, C.J.; BOUZON, G.; SIMAS, H. (1998) Aspectos em identificação, projeto e experimentação de algoritmos de controle adaptativo e predito indireto e direto. Congresso Brasileiro de Automática, 12, Anais... Uberlândia: CBA.) e os processos de cálculos estabelecidos pelo modelo ARX, que resultou nas Equações 16 e 17. A Equação 17 é composta por duas parcelas. A primeira (entre colchetes) é composta por informações disponíveis até o instante k, e a segunda, por informações futuras, no instante k+d.
onde Φ(k+d) é o sistema generalizado para o instante d passos à frente do instante t; P, Q e R são polinômios de ponderação sobre os respectivos sinais de saída, controle e referência; A, B, F e G são polinômios discretos de parâmetros; w(k) é a referência para a saída do sistema original no instante de tempo k; e(k+d) é o erro no instante (d+k); z-1 é o operador de atraso, y(k)z-1=y(k-1). Os polinômios P, Q e R foram ajustados de acordo com os critérios de desempenho estabelecidos.
O mínimo da função custo é obtido igualando o termo entre colchetes da Equação 17 a zero. Assim, a lei de controle de GVM é dada pela equação 18. Observa-se que a variável de controle é função dos valores de entrada e saída e da referência a ser seguida, obtidos no instante k, além de ser dependente dos parâmetros do sistema original, os quais são identificados recursivamente (VAZ, 1999VAZ, G.L.F. (1999) Controladores adaptativos de Variância Mínima e Dahlin: uma revisão e novas concepções de projeto. Dissertação (Mestrado) - Universidade Federal de Santa Catarina, Florianópolis.). Estruturando o GMV e substituindo as equações encontradas tem-se o diagrama mostrado na Figura 2 para a dinâmica do sistema em malha fechada. Este é descrito pela Equação 19.
Os parâmetros P, Q e R, presentes na Equação 19, devem ser selecionados pelo projetista. Neste caso, adotou-se R(z-1) = P(1) e Q(z-1) = q0(1-z-1), de maneira a garantir um comportamento transitório e de regime permanente adequado ao sistema em malha fechada. O procedimento de projeto do controlador consistiu em calcular o valor da variável de controle, supondo os parâmetros do sistema conhecidos e, em seguida, fazendo a sua substituição pelos valores dos parâmetros estimados, de forma recursiva.
Os polinômios A1(z-1), B1(z-1), A2(z-1) e B2(z-1) compõem as funções de transferência dos ramais e são dados pelos itens da Equação 20. Os polinômios são compostos pelos parâmetros que representam a dinâmica "própria" de cada ramal. Os polinômios F e G são representados pelas Equações 21 e 22, respectivamente.
onde an, bn, fn e gn são parâmetros discretos dos polinômios.
Realizando as devidas manipulações nas equações e considerando os parâmetros ai e bj estimados para estes modelos, obtêm-se as Equações 23 e 24 representando as leis de controle dos ramais do sistema experimental. Tais equações foram adotadas no controlador GMV autoajustável implementado no ambiente LabVIEW.
onde yi=wi(k)-Ei(k); wi(k) é a referência preestabelecida para o ramal i; ei(k) é o sinal de erro de posicionamento referente ao ramal i.
DESCRIÇÃO DA BANCADA DE TESTES E DOS ENSAIOS
O sistema de controle adaptativo foi implementado e testado em um sistema experimental de distribuição de água (Figura 3) localizado no Laboratório de Eficiência Energética e Hidráulica em Saneamento (LENHS) da Universidade Federal da Paraíba (UFPB), e é constituído de um reservatório circular de 7,64 m3, um conjunto motor-bomba (CMB) de 3 cv, um inversor de frequência, duas válvulas de controle (CVc e CVs), dois transdutores de pressão (PT1 e PT2), dois medidores de vazão, fontes de alimentação, manômetros tipo bourdon , além das conexões e tubos de PVC (DN 50) da rede hidráulica.
A rede hidráulica contém dois ramais que emulam a alimentação de dois setores em sistemas de distribuição de água reais. A CVs, a jusante do PT1, tem a função de alterar as condições de operação do sistema, proporcionando a variação da vazão demandada. À medida que a CVs fecha, a vazão diminui e, consequentemente, a pressão a montante aumenta. O inversor de frequência será responsável por manter a pressão no ponto mais desfavorável igual à pressão desejada (pressão de referência), enquanto a CVc elimina o excesso de pressão no ponto de medição do ramal 1.
Realizaram-se três experimentos para atestar a eficiência do controlador adaptativo e o seu impacto no consumo de energia elétrica e água, cuja diferença entre estes é a abertura da CVs. Os experimentos 1, 2 e 3 foram realizados com a CVs nas posições 22,5°, 45,0° e 67,5°, respectivamente.
Os ensaios experimentais foram realizados com o sistema em malha fechada para uma entrada tipo degrau igual a 20 mca para os dois pontos de monitoramento de pressão. Estes iniciaram com a frequência da tensão de alimentação do motor em 30 Hz e a CVc totalmente fechada. O tempo de duração dos experimentos foi de 5 minutos e a aquisição de dados foi feita a cada 2 segundos, gerando uma massa de 150 amostras por ensaio.
RESULTADOS E DISCUSSÃO
Foram realizados diversos testes, em que foi constatado que as respostas seguem, em cada experimento, um mesmo padrão. O inversor de frequência, atuando de maneira isolada, não tem capacidade de otimizar a pressão em mais de um ponto de um sistema de distribuição de água. Desse modo, o inversor foi responsável por manter a pressão no ponto de monitoramento mais desfavorável (maior déficit de pressão) no valor desejado, enquanto a válvula de controle eliminou o excesso de pressão a sua jusante.
O comportamento dos experimentos 1, 2 e 3 é representado nas Figuras 4 a 6, respectivamente. Observa-se que enquanto as pressões do sistema são inferiores à pressão de referência, o controlador atua aumentando a velocidade de rotação do CMB e o ângulo de abertura da válvula CVc, até que a pressão nos ramais atinja o valor da referência (20 mca). Na Tabela 1 é apresentada a sinopse dos resultados dos três experimentos. Os valores dos parâmetros demonstram a aplicabilidade do controlador adaptativo GMV autoajustável no controle de sistemas de distribuição de água.
Na prática, variações de pressão inferiores a 2 mca são plenamente aceitas na operação de redes hidráulicas. Nas Figuras 7 e 8 são mostradas, respectivamente, as curvas de resposta da CVc e do inversor de frequência. Nas Figuras 9 a 11 é mostrado o comportamento do erro nos pontos de medição dos experimentos. O erro máximo no regime estacionário dos ensaios foi de 2,12%, que equivale a 0,42 mca.
Após a verificação da aplicabilidade do sistema de controle, analisou-se o comportamento da potência elétrica do motor para a CVs nas posições 0°, 22,5°, 45°, 67,5° e 90° (Figura 12).
Comparando a energia consumida pelo sistema experimental operando com e sem a atuação do sistema de controle tem-se uma redução de 19,2% na energia consumida e de 16,7% no consumo específico de energia elétrica para a condição mais desfavorável de ganho de energia (vazão máxima). À medida que a CVs é fechada, a vazão no setor diminui e, consequentemente, há uma exigência menor da bomba, resultando na diminuição significativa do consumo de energia elétrica à medida que o ângulo de abertura da CVs diminui.
O controle da pressão é a principal ação para o controle de perdas reais (físicas) em sistemas de distribuição de água (BEZERRA et al., 2012BEZERRA, S.T.M.; SILVA, S.A.; GOMES, H.P. (2012) Operational optimisation of water supply networks using a fuzzy system. Water SA , v. 38, n. 4, p. 565-572.; BABIĆ et al., 2014BABIĆ, B.; DUKIĆ, A.; STANIĆ, M. (2014) Managing water pressure for water savings in developing countries. Water SA , v. 40, n. 2, p. 221-232.). A manutenção da pressão próxima a valores predefinidos em redes reais proporciona, também, a redução no volume de água perdido por vazamentos. Segundo Babić et al. (2014BABIĆ, B.; DUKIĆ, A.; STANIĆ, M. (2014) Managing water pressure for water savings in developing countries. Water SA , v. 40, n. 2, p. 221-232.), uma redução de 10% na pressão do sistema resulta, geralmente, em uma diminuição no volume perdido em vazamentos entre 5 e 15%, podendo atingir, em alguns casos, uma redução de 23%.
CONCLUSÃO
Conclui-se que o controlador adaptativo GMV autoajustável projetado é estável e robusto, tendo apresentado um bom desempenho no controle da pressão e uma redução satisfatória no consumo de energia elétrica na operação da bancada experimental apresentada. Extrapolando a aplicação para sistemas reais e considerando que uma das principais ações de controle de perdas reais de água é o controle de pressão, a manutenção da pressão próxima a valores predefinidos em redes reais resultaria, também, em uma diminuição significativa no volume de água perdido em vazamentos.
Devido às características adaptativas do sistema desenvolvido, espera-se que este se ajuste satisfatoriamente às condições operacionais impostas em redes hidráulicas reais, que, comumente, possuem múltiplos ramos e derivações, muitas vezes sem setorização. De forma geral, o sistema de controle garante que os valores de pressão nos pontos de medição permaneçam próximos aos valores de referência, porém ressalta-se que, em aplicações reais, o desempenho global do sistema está diretamente relacionado com a definição desses pontos pelo gestor da operação.
REFERÊNCIAS
- ÅSTRÖM, K.J. & WITTENMARK, B. (1973) On self-tuning regulators. Automatica , v. 9, n. 2, p. 185-199.
- BABIĆ, B.; DUKIĆ, A.; STANIĆ, M. (2014) Managing water pressure for water savings in developing countries. Water SA , v. 40, n. 2, p. 221-232.
- BEZERRA, S.T.M.; SILVA, S.A.; GOMES, H.P. (2012) Operational optimisation of water supply networks using a fuzzy system. Water SA , v. 38, n. 4, p. 565-572.
- CARVALHO, J.S.F. (2009) Controle adaptativo aplicado em um robô manipulador de dois graus de liberdade planar . Tese (Doutorado) - Universidade Federal da Paraíba, João Pessoa.
- CLARKE, D.W. & GAWTHROP, P.J. (1975) Self-Tuning controller. Proceedings of the Institution of Electrical Engineers , v. 122, n. 9, p. 929-934.
- ELETROBRÁS - Centrais Elétricas Brasileiras S.A. (2008). Disponível em: <http://www.eletrobras.com>. 2008. Acesso em: 28/10/2008.
» http://www.eletrobras.com - HUUSOMA, J.K.; POULSENB, N.K.; JØRGENSENA, S.B.; JØRGENSENB, J.B. (2012) Tuning SISO offset-free Model Predictive Control based on ARX models. Journal of Process Control , v. 22, n. 10, p. 1997-2007.
- LIUZZO, S. & TOMEI, P.A. (2008) A global adaptive learning control for robotic manipulators. Automatica , v. 44, n. 5, p. 1379-1384.
- MARCHI, P.A.; ROCHA, C.R.; VENTURO, C.J.; BOUZON, G.; SIMAS, H. (1998) Aspectos em identificação, projeto e experimentação de algoritmos de controle adaptativo e predito indireto e direto. Congresso Brasileiro de Automática, 12, Anais... Uberlândia: CBA.
- PELLI, T. & HITZ, H.U. Energy indicators and savings in water supply. (2000) Journal AWWA , v. 92, n. 6, p. 55-62.
- PURWARA, S.; KARB, I.N.; JHAB, A.N. (2008) Adaptive output feedback tracking control of robot manipulators using position measurements only. Expert Systems with Applications , v. 34, n. 4, p. 2789-2798.
- SANDERS, K.T. & WEBBER, M.E. (2012) Evaluating the energy consumed for water use in the United States. Environmental Research Letters , v. 7, p. 1-11.
- VAZ, G.L.F. (1999) Controladores adaptativos de Variância Mínima e Dahlin: uma revisão e novas concepções de projeto. Dissertação (Mestrado) - Universidade Federal de Santa Catarina, Florianópolis.
- WANGA, J.Y.; CHANG, T.P.; CHEN, J.S. (2009) An enhanced genetic algorithm for bi-objective pump scheduling in water supply., Expert Systems with Applications v. 36, p. 249-258.
-
1
Trabalho realizado no Laboratório de Eficiência Energética e Hidráulica em Saneamento do Centro de Tecnologia da Universidade Federal da Paraíba (UFPB) - João Pessoa (PB), Brasil.
-
Reg. ABES:
107226
Datas de Publicação
-
Publicação nesta coleção
Set 2015
Histórico
-
Recebido
27 Nov 2012 -
Aceito
01 Dez 2014