Modelos de análise hidráulica de redes de distribuição de água considerando demanda dirigida pela pressão e vazamentos

Water distribution networks analysis considering leakage and pressure driven demand modeling

ARTIGO TÉCNICO

Modelos de análise hidráulica de redes de distribuição de água considerando demanda dirigida pela pressão e vazamentos

Water distribution networks analysis considering leakage and pressure driven demand modeling

Klebber Teodomiro Martins Formiga^I; Fazal Hussain Chaudhry^II

^IEngenheiro Civil –UFCG. Mestre em Engenharia Civil - UFCG. Doutor em Hidráulica e Saneamento EESC/USP. Professor Adjunto da Escola de Engenharia Civil UFG

^IIEngenheiro Civil - Punjab University. Mestre em Engenahria Hidráulica - Asian Institute Of Technology. Doutor em Engenharia Civil - Colorado State University. Professor Titular do Departamento de Hidráulica e Saneamento EESC/USP

^{Endereço para correspondência}

RESUMO

As ferramentas para análise hidráulica são componentes importantes na avaliação do funcionamento das redes de distribuição de água para abastecimento. Existem diversos métodos que podem ser utilizados para essa análise, no entanto, os modelos que procuram resolver o sistema de equações correspondente através do método Newton-Raphson ou por meio de linearizações sucessivas são os mais eficientes. Quatro formulações baseadas nestes esquemas são avaliadas neste trabalho. O objetivo deste trabalho é fazer uma comparação dos métodos Newton-Raphson, Teoria Linear, Híbrido e Gradiente para a análise de redes de distribuição de água em regime permanente, considerando a demanda dirigida pela pressão e os Vazamentos. Para tanto, foi utilizado um layout de rede frequentemente empregado na literatura dotado de válvulas. O método do Gradiente foi o que convergiu em um número menor de iterações para redes mais simples, o Método Híbrido foi o que mais se adaptou para sistemas mais complexos.

Palavras-chave: Modelagem, redes de distribuição de água, vazamentos, demanda dirigida pela pressão.

ABSTRACT

The hydraulic analysis tools are important in the performance evaluation of water distribution networks. Various methods are available for such analysis. However, the hydraulic models that solve the system of equations describing the flow problem through Newton-Raphson or through its successive linearizations are the most efficient. It is the purpose of this paper to compare Newton-Raphson, Linear Theory, Hybrid and Gradient methods for steady-state hydraulic network analysis, considering leakage and pressure driven demand modeling. The network layout with hydraulic components frequently used in the literature was employed for this analysis. The Gradient method was found to produce best results in terms of number of iterations for the more simple networks, whereas the Hybrid method was better for more complex networks.

Keywords: Modeling, water distribution networks, leakage, pressure driven demand.

INTRODUÇÃO

A análise de redes de distribuição de água consiste na forma mais eficiente de previsão do comportamento do sistema – vazões e perdas de carga nos trechos, pressões e cargas hidráulicas nos nós – quando ocorrem modificações devido a incrustações, mudança da rugosidade em virtude do envelhecimento dos tubos, vazamentos e quebras, entre outras situações comuns às redes. O problema hidráulico correspondente pode ser posto como o de determinação das vazões nos diversos trechos da rede, de maneira a satisfazer as condições expressas através das equações da conservação de massa nos nós e de energia nos trechos da rede.

Como as equações do problema não são lineares, este só pode ser resolvido diretamente em apenas algumas situações especiais. No entanto, na maioria dos casos, é preciso empregar métodos iterativos assumindo uma solução aproximada a cada passo até que o problema tenha sido resolvido, tais métodos podem ser classificados da seguinte forma: a) gradiente local; b) Newton-Raphson; c) linearização; d) minimização numérica; e e) métodos híbridos.

Os métodos pertencentes às três primeiras classes utilizam ferramentas de resolução de sistemas de equações não lineares para achar a resposta do problema, enquanto que o quarto método transforma o problema na minimização de uma função não linear convexa, considerando restrições de igualdade e desigualdade lineares. Os métodos híbridos são resultados da união de dois ou mais métodos anteriores.

Os métodos Hardy Cross (MHC) baseiam-se no trabalho de Cross (1936) em que o sistema de equações original é resolvido empregando o conceito de gradiente local. Embora sejam os mais fáceis de calcular manualmente, os MHC e suas derivações podem apresentar um comportamento instável, não convergindo para redes com muitas malhas, sendo muito dependente da escolha dos anéis para se efetuar o balanço de energia. Além disso, o método perde eficiência numérica quando aplicado às redes de grande porte, pois apresenta uma taxa de convergência muito baixa (Lemieux, 1972). Outra desvantagem desta formulação está na dificuldade de se acrescentar dispositivos na rede como bombas e válvulas redutoras de pressão.

A segunda classe de métodos resolve o sistema de equações não lineares por meio do método Newton-Raphson – MNR (Shamir e Howard, 1968). Este método, embora mais eficiente que o MHC, apresenta deficiências inerentes à formulação que depende principalmente dos valores iniciais para que apresente uma convergência rápida. A falta desses valores pode ser crítica em sistemas mais complexos.

Os métodos de linearização das equações, com destaque para o Método da Teoria Linear – MTL (Wood e Charles, 1972; Isaacs e Mills, 1980), empregam aproximações lineares sucessivas da equação de perda de carga nos trechos, produzindo bons resultados. A principal vantagem deste método são suas propriedades de convergência, que diferentemente do MNR, não requer boas soluções iniciais. No entanto a taxa de convergência das iterações finais deste método é inferior ao MNR (Nielsen, 1989).

O emprego de técnicas de otimização na análise de redes é feito de modo a minimizar uma função objetivo não linear e convexa, sujeita a um conjunto de equações lineares de restrições que representam as equações de conservação de massa nos nós e energia nos trechos. A principal vantagem deste método está na formulação, uma vez que a convexidade da função objetivo, combinada com a linearidade de equações garante a existência e a unicidade da solução. Infelizmente, a solução numérica para este problema requer a utilização de um método de otimização que possa ser eficiente até mesmo quando aplicado a redes de grande porte. Embora as técnicas de otimização não lineares estejam avançadas atualmente, estas ainda enfrentam problemas de dimensionalidade, ou seja, para redes de grande porte, a capacidade e o tempo de processamento necessário são muito elevados. Exemplos de métodos que utilizam otimização para a análise hidráulica de redes são os trabalhos de Collins, Cooper, Helgason et al (1978) e Gomes (2004) que empregam respectivamente o Convex Simplex Method e o Generalized Reduced Gradient (GRG2), na resolução deste problema.

Os métodos híbridos incorporam vantagens de dois, ou mais métodos acima citados de modo a reduzir as limitações dos mesmos. O Método do Gradiente (Todini e Pilati, 1987) efetua a união de técnicas de minimização numérica com o método Newton-Raphson. Esta união resultou em um algoritmo muito eficiente e que implicou na sua utilização em quase todos os pacotes de simuladores hidráulicos atualmente disponíveis como: EPANET, WATERCAD, SARA, MIKENET, PIPENET. Outro método híbrido foi o proposto por Nielsen (1989) que uniu as melhores características do MNR e MTL, resultando em um método robusto, capaz de convergir, em um número muito menor de iterações, quando comparado com os resultados obtidos com essas técnicas isoladas (Souza, 1994).

Embora a adoção destes modelos de análise tenha avançado nos últimos anos, eles ainda carecem de um melhor refinamento que possibilite simular de modo mais realista as condições reais do sistema (Todini, 2003). Estes requisitos são ainda mais importantes em estudos mais específicos da rede, como análise de perdas por vazamentos e avaliação da confiabilidade do sistema. Para tanto, é necessário efetuar uma modificação dos métodos desenvolvidos até então. Esta é feita incorporando modelos de vazamentos e de demanda dirigida pela pressão (DDP) à formulação hidráulica tradicional (Formiga, 2005).

Considerando que o vazamento é consequência direta de defeitos na rede, pode-se caracterizá-lo baseado na equação de descarga em orifícios (Jowitt e Xu, 1990). Desta forma, para avaliar as perdas o método de análise hidráulica deve considerar que ocorrerão fugas de água na rede e que estas serão proporcionais à pressão. Alguns métodos de avaliar estas perdas foram propostos por Germanopoulos e Jowitt (1989), Martínez, Conejos e Vercher (1999) e Tucciarelli, Criminisi e Termini (1999).

A finalidade básica dos sistemas de distribuição de água é transportar água de um lugar para o outro, garantindo que esta chegue em quantidade e qualidade desejadas ao usuário final. Assim, ao disponibilizar água ao consumidor, é preciso que exista ainda uma sobra de energia que será usada para superar as perdas de carga dentro dos imóveis. Se a pressão não for adequada, a quantidade de água suprida será menor do que aquela requerida, mas, não será nula, havendo então um atendimento parcial das demandas. No entanto, se a pressão estiver abaixo de determinado valor, que depende de cada caso, não haverá suprimento de nenhuma demanda (Wagner, Shamir e Marks, 1988).

Desta forma, ao se tentar modelar este comportamento, deve-se levar em conta a diminuição das demandas atendidas, o que modifica em parte o cenário que se está simulando. A forma de se avaliar esse fenômeno é considerar a demanda também como uma variável, que teria o seu valor calculado em função da pressão nos nós, ou seja, tem-se a demanda dirigida pela pressão (DDP). Algumas formulações de DDP foram propostas por Wagner, Shamir et al (1988), Reddy e Elango (1989) e Fujiwara e Ganesharajah (1993). Gupta e Bhave (1996) fizeram um comparativo de diferentes técnicas que calculam a demanda como dependente da pressão, concluindo que os métodos parabólicos foram os que apresentaram melhores resultados.

A introdução dos modelos de vazamento e demanda dirigida pela pressão (DDP), no algoritmo de análise hidráulica foi proposta teoricamente por Salgado, Rojo e Zepeda (1993) para o método gradiente. No entanto, a adoção direta deste provoca instabilidade no modelo (Ackley, Tanyimboh, Tahar et al, 2001), podendo inviabilizar a sua utilização em algumas situações. Para se contornar este problema, foram desenvolvidas alguma rotinas para o cálculo da DDP e dos vazamentos (Martínez, Conejos et al, 1999; Todini, 2003, Soares, Cheung, Reis et al, 2004), no entanto, nestes trabalhos, o cálculo dos vazamentos era efetuado fora do algoritmo de análise hidráulica, devendo os valores das demandas serem recorrentemente atualizados. Embora este procedimento possibilite a convergência dos resultados, a sua eficiência é limitada, visto que o número de iterações necessário ao cálculo é elevado.

Este trabalho tem por finalidade efetuar a inclusão da DDP e dos vazamentos explicitamente na formulação dos métodos de análise hidráulica, verificando a eficiência e estabilidade numérica dos mesmos.

MÉTODOS DE ANÁLISE HIDRÀULICA DE REDES

A seguir será feita uma explanação sobre a formulação dos métodos avaliados neste trabalho. No desenvolvimento do modelo de análise hidráulica será considerada a seguinte convenção: para um trecho i que conecte dois nós j_i,1 e j_i,2 Seja uma rede qualquer composta por NN nós interiores, NR nós fontes (reservatórios ou tanques), e NT trechos que conectam os nós. A perda de carga entre dois nós, j_j,1 e j_j,2, será expressa por:

em que h é a carga hidráulica no nó j, Ki é o coeficiente de resistência da tubulação no trecho i; a é o expoente da fórmula de perda de carga.

Método Newton-Raphson (MNR)

O método Newton-Raphson (Shamir e Howard, 1968) para a resolução de sistemas de equações não lineares é desenvolvido a partir da seguinte equação:

em que h é o vetor de carga hidráulica nos nós da rede na iteração k, e Dh é o gradiente da função de balanço de massa nos nós, obtido a partir da matriz Jacobiana conforme procedimento demonstrado por Nielsen (1989).

Desenvolvendo este sistema (Souza, 1994) obtêm-se a equação matricial expressa por:

em que A é a matriz de incidência da rede, a é o expoente da fórmula de perda de carga, Q é o vetor de demanda nos nós, h_r é a carga disponível nos reservatórios do sistema, e J_k é a matriz Jacobiana expressa por :

J_k = A^TCA

em que C é a matriz diagonal:

O sistema expresso em (3) é resolvido sucessivamente até que a solução encontrada atenda a uma precisão pré-determinada, ou seja, a diferença das variáveis entre duas iterações seja pequena.

Método da Teoria Linear (MTL)

O método da Teoria Linear (MTL) foi proposto inicialmente por Wood e Charles (1972) para formulação dos anéis, e modificado posteriormente por Isaacs e Mills (1980) para a formulação dos nós. Este método lineariza (1), resultando:

O vetor da carga hidráulica que será avaliado iterativamente por:

Método Híbrido (MH)

Os métodos MNR e MTL apresentam vantagens para se trabalhar com determinados tipos de problemas. De um modo geral, o MNR converge mais rápido (quadraticamente) quando comparado com o MTL, porém necessita de boas estimativas iniciais das cargas nos nós para que isto ocorra. Caso a aproximação inicial não esteja nas proximidades da solução final, o MNR pode consumir mais iterações para convergir, podendo até não fazê-lo (Nielsen, 1989). Essa dificuldade é contornada implementando-se sub-rotinas para a determinação de melhores valores da solução inicial (Souza, 1994).

O MTL apresenta uma taxa de convergência inicial melhor do que o MNR, não precisando de valores iniciais tão bons. No entanto, o MTL tende a oscilar em torno da solução final, apresentando dificuldades de convergência para sistemas complexos (Hansen, Madsen e B, 1991).

O método proposto por Nielsen (1989), doravante chamado de híbrido (MH), faz inicialmente algumas iterações, com o MTL, obtendo-se uma melhor aproximação das vazões e cargas hidráulicas nos nós. A aproximação assim obtida, será a solução inicial do MNR, que tende a convergir em um número menor de iterações. Assim, o MTL serviria como uma rotina para determinação dos valores de partida para o método de Newton-Raphson propriamente dito.

Como as formulações do MTL e MNR são semelhantes (3) e (6), a implementação do método Híbrido é realizada de acordo com a expressão (7):

em que g = 1 quando k < ii e g = a quando k > ii, onde ii é o número de iterações iniciais com o MTL.

Método Gradiente (MG)

As equações de conservação de energia podem ser escritas matricialmente da seguinte forma (Todini e Pilati, 1987):

enquanto que a conservação de massa é expressa por:

O sistema de equações resultante é dado por:

Este sistema (10) é não linear, podendo ser resolvido através do esquema Newton-Raphson, para isso, é necessário efetuar uma diferenciação do sistema com relação à q e h, que resulta em:

Após alguns procedimentos algébricos aplicados à expressão (11), encontram-se as seguintes equações:

MODIFICAÇÃO DOS MÉTODOS DE RESOLUÇÃO

A inclusão dos eventuais vazamentos na rede e da utilização da demanda como uma variável dependente da pressão nos nós fará com que o vetor demanda (Q) seja função também da altura hidráulica (h). Desta forma, é preciso reformular o sistema de equações original para atender a essa nova situação. As modificações no método de análise hidráulica são apresentadas a seguir.

A formulação empregada para vazamento neste trabalho foi proposta por Tucciarelli, Criminisi et al (1999), em que os coeficientes de perdas são considerados constantes por unidade de superfície dos tubos pertencentes a determinadas áreas da rede. Essa suposição é importante no processo de calibração de redes, em que se procura trabalhar com o menor número de parâmetros possível, além de ser melhor ajustada à formulação hidráulica de nós, que é empregada neste trabalho. Os autores consideraram que o vazamento pode ser expresso por:

sendo z_j a cota do nó j, _ij é o coeficiente de perdas por unidade de superfície do tubo entre os nós i e j; q é o expoente de vazamento tomado como sendo 1,18 (Jowitt e Xu, 1990), e d_ij é o diâmetro da tubulação.

A expressão empregada para DDP, conhecida como parabólica, foi proposta por Wagner, Shamir e Marks (1988) de acordo com o seguinte equacionamento:

em que Q^dem é a demanda atendida; Q^req é a demanda potencial; h^min é a carga hidráulica mínima necessária haja suprimento de água, ou seja a cota mínima em que há consumo; h^req é a carga hidráulica requerida para que a demanda seja completamente suprida e n é o expoente tomado como 2.

A adoção da vazão como sendo dependente da carga no nó torna necessário introduzir uma mudança na formulação original dos métodos para a análise hidráulica da rede. Nos métodos MNR, MTL E MH esta mudança ocorrerá na definição da matriz Jacobiana que será:

em que CVaz será avaliado pela formulação de Tucciarelli, Criminisi et al, 1999:

Esta expressão representa a derivada de 1ª ordem do termo de vazamento com relação à carga hidráulica h_j de (14).

A matriz C_Dem, que representa a derivada da equação de demanda dependente da pressão, é expressa por:

em que:

Salgado, Todini e O’Connell (1993) propuseram uma modificação no MG de modo a considerar um componente da demanda como sendo uma função da carga hidráulica no nó.

em que:

Desenvolvendo de modo semelhante ao feito no MG, obtêm-se as seguintes expressões:

em que:

Válvulas

Neste trabalho será avaliada uma rede com diferentes configurações. Estas configurações dizem respeito à consideração dos vazamentos, DDP e válvulas. As válvulas empregadas neste estudo são: válvula de retenção (VR) e válvula redutora de pressão, que têm o seu equacionamento apresentado a seguir.

Válvula de Retenção (VR): têm por finalidade fazer com que a água flua em apenas uma direção, fechando-se quando o fluxo começa a reverter. Seja um trecho i, com resistência igual a K_i, que possui uma válvula de retenção instalada. Esse dispositivo hidráulico pode operar de duas formas diferentes: a) fechada - se a altura de jusante for maior que a de montante a válvula será ativada fazendo com que: K_i =+¥ q_i = 0. b) aberta - a VR estará inativa permanecendo inalterada a resistência do trecho.

Válvula Redutora de Pressão (VRP): servem para limitar a pressão à jusante, mantendo-a abaixo de valores que podem ser prejudiciais à rede. A VRP pode estar em três situações: a) inativa - quando a carga a montante da válvula for inferior à pressão de serviço, nesta situação as condições do sistema permanecem inalteradas; b) fechada - quando a carga de jusante é maior que a de montante, neste caso a VRP se comportará como se fosse uma VR atuando uma resistência infinita c) ativa - em que permanece parcialmente fechada reduzindo a pressão à jusante da válvula, neste caso a resistência do trecho será:

em que K^valv é a resistência da válvula expressa por:

Termo de relaxação para equação de demanda dirigida pela pressão

Uma das dificuldades apontadas quando se trabalha com DDP diz respeito à descontinuidade da formulação (15), o que tem inviabilizado a adoção do equacionamento explícito para o método. Para exemplificar isto, considere um sistema adutor simples conforme a Figura 1. Para a formulação original proposta (15) tem-se que na primeira iteração, a carga hidráulica no nó será menor do que a mínima requerida (Figura 1a). Assim na iteração seguinte, a demanda no nó será nula, e a pressão será a mesma do reservatório (Figura 1b). Este processo se repetirá indefinidamente fazendo a rotina entrar em um "loop infinito", impedindo assim a convergência para uma solução.

Se em uma situação relativamente simples, esta ocorrência implica na falta de convergência do método, para redes maiores, que contenham dispositivos como válvulas e bombas, tal fato pode inviabilizar completamente este tipo de abordagem. Assim, para se superar tal empecilho, este trabalho propõe uma relaxação no cálculo da DDP de modo a permitir a convergência, em que a demanda na iteração i seria calculada de seguinte forma:

Esta equação seria válida até que o método estabilizasse, pequeno erro entre duas iterações (neste trabalho empregou-se 0,1), a partir do qual voltaria a ser empregada a formulação original até atingir o critério de convergência.

REDE EMPREGADA NO ESTUDO

Para avaliar a eficiência e o comportamento dos modelos hidráulicos abordados neste trabalho foi escolhida uma rede exemplo usada por Souza (1994), que tem o seu layout apresentado na Figura 2. Esta rede é formada por 27 trechos conectados por 21 nós de consumo e 1 reservatório. A principal característica dessa configuração está no fato dela possuir os componentes hidráulicos que serão também avaliados neste trabalho: válvulas de retenção e redutora de pressão, sendo, ideal para a avaliação completa do comportamento de modelos de análise hidráulica. Os dados referentes aos nós da rede – cota, demanda, cargas hidráulicas desejada e mínima e coeficiente de vazamento – são apresentados na Tabela 1, enquanto que as informações dos trechos – comprimento, diâmetros e coeficiente linear de vazamento – são encontradas na Tabela 2.

O foco das simulações efetuadas neste trabalho é verificar o funcionamento das metodologias de análise hidráulica da rede considerando a demanda dependente da pressão e os vazamentos da rede. Por este motivo, foram feitas pequenas modificações nos dados originais. Quanto aos componentes da rede, a válvula redutora de pressão (VRP) funcionará quando a pressão à montante for superior a 90 e 25 mH₂O_, nas simulações com e sem DDP respectivamente. A rugosidade dos tubos foi tomada como sendo igual a 1 mm, uma vez que a equação empregada será a de Darcy-Weissbach para perda de carga, implicando em um valor de 2 para o expoente a de (3). Foram adotados dois níveis de água nos reservatórios, 50 m nas simulações em que a DDP foi empregada e 100 m nas restantes. O critério de convergência escolhido foi o de que a maior diferença entre as cargas hidráulicas nos nós obtidas em de duas iterações consecutivas, fosse menor que uma tolerância (t) pré-estabelecida, neste trabalho utilizou-se t = 10^-5. A formulação apresentada para os métodos descritos neste trabalho foi implementada em linguagem de programação C++, utilizando o compilador Borland C++ Builder 6.

Os métodos de balanço hidráulico de redes foram testados para um total de 16 configurações de componentes e diferentes comportamentos da rede quanto à demanda dependente da pressão e vazamentos (Tabela 3)

RESULTADOS

Os métodos comparados neste trabalho apresentam formulações semelhantes, sendo diferentes em poucos detalhes. Do ponto de vista de eficiência computacional, o tempo gasto por esses métodos a cada iteração é muito similar, uma vez que os sistemas de equações lineares, que consomem maior tempo computacional, resolvidos em cada passo são do mesmo tamanho (nxn) e possuem as mesmas características matriciais. Dessa forma, aferiu-se a eficiência das metodologias para análise hidráulica de rede a partir do número de iterações necessárias para que o método convergisse, considerando a tolerância pré-definida.

Os números de iterações necessários para a convergência dos métodos para cada configuração são mostrados na Tabela 4. Os valores das cargas hidráulicas nos nós e das vazões nos trechos são mostrados nas Tabelas 5 e 6.

A partir dos resultados apresentados na Tabela 4 verifica-se que o MG foi o mais eficiente para as configurações mais simples que não consideram a DDP e os vazamentos (Configurações 1 a 4), bem como nas situações em que a DDP e vazamentos foram consideradas sem a presença de componentes hidráulicos (Configurações 5, 9 e 13). No entanto, para as outras redes o desempenho do método foi inferior ao MH e em alguns casos (Configurações 6, 8, 14 e 16), não apresentou convergência. Verifica-se assim uma grande capacidade do MG para redes mais simples e uma dificuldade de encontrar soluções para sistemas mais complexos, onde ocorrem descontinuidades no sistema, geradas principalmente pela presença de VRPs.

Para configurações mais complexas (Configurações 3, 4, 7, 11, 12, 13 e 15), o MH foi o que apresentou melhor convergência (Tabela 4). Isto decorre especialmente pela estabilidade numérica proporcionada pelo MTL, que é empregado nas primeiras iterações daquele método. Embora o MH não tenha sido o método que convergiu em menos iterações nas outras configurações, este não ficou muito distante dos métodos mais rápidos, conseqüência da sua capacidade em lidar com sistemas simples ou complexos.

O MTL foi o mais lento na análise de configurações menos complexas, chagando a gastar 3 vezes mais tempo do que o MG (Configurações 1, 5 e 9). Entretanto, o MTL mostrou-se o mais rápido nas situações mais complexas analisadas. Isto se deve principalmente à estabilidade do método que permite uma aproximação gradual, porém estável, da solução final e permite ao método lidar bem com as descontinuidades das funções resultante das configurações mais complexas. No caso dos métodos mais rápidos, como o Gradiente e Newton-Raphson, pode-se incorrer em dificuldades tais como instabilidade numérica e não convergência. Mesmo que os equacionamentos do MNR (3) e do MTL (6) sejam muito semelhantes, o MTL convergiu para todas as situações enquanto que o MNR apresentou dificuldade de convergência para os casos mais complexos.

A relaxação introduzida melhorou a convergência do MG nas Simulações 3 e 4 e permitiu a convergência do MG e do MNR para a situação de DDP (Simulação 5 e 13). No entanto, a relaxação da função DDP não alterou a convergência dos métodos MTL e MH. Isto se deve ao fato do MTL já ter incorporado em sua formulação o principio de cálculo das vazões de cada interação levando-se em conta a vazão da iteração anterior.

Das configurações analisadas, aquelas que possuem válvulas redutoras de pressão foram as que requereram maior esforço computacional, fato motivado pelo tipo de condição de contorno que este componente apresenta, com três possíveis estados e descontinuidade nas derivadas de primeira ordem. As condições de contorno das válvulas de retenção, embora mais simples, são descontínuas, operando em apenas dois tipos de condições: aberta ou fechada, provocando também um aumento no número de iterações.

CONCLUSÕES

Os modelos de análise hidráulica de redes de distribuição são ferramentas essenciais nos estudos modernos de projeto e operação dos sistemas de abastecimento, uma vez que eles permitem avaliar diferentes cenários de funcionamento e auxiliam o processo de tomada de decisão. Para que o uso destes modelos seja ampliado é preciso que os mesmos sejam adaptados para descrever a situações mais realísticas que os sistemas, principalmente os brasileiros, enfrentam atualmente, como pressões insuficientes em alguns pontos da rede e vazamentos.

Esse trabalho procurou comparar as principais formulações de análise hidráulica de redes em regime permanente desenvolvidas ao longo das últimas décadas para situações de demanda dirigida pela pressão (DDP) e vazamentos na rede. As formulações MNR (Mét. Newton-Raphson), MTL MTL (Mét.Teoria Linear), MH (Mét.Híbrido) e MG (Mét.Gradiente) foram modificadas de modo a considerar os vazamentos e a DDP dentro do modelo e implementadas em linguagem de programação C++.

O MG mostrou ser o mais robusto, uma vez que trabalha diretamente com as duas variáveis envolvidas no problema, carga hidráulica e vazão. Isto se verifica pela sua melhor adaptação a configurações menos complexas, como nas situações em que a DDP e os vazamentos não são considerados, ou quando estes ocorrem individualmente sem a presença de válvulas. Esta eficiência do MG o faz estar presente na maior parte dos softwares disponíveis no mercado. Ao se avaliar cenários mais complexos, o MG apresentou alguns problemas de convergência, havendo casos em que o mesmo não convergiu, mesmo após a incorporação da rotina de relaxação.

O MH mostrou-se o método que melhor se adapta aos diferentes cenários de complexidade, pois alia as características positivas do MTL, de estabilidade para as iterações iniciais, e do MNR, relativas à convergência para a solução final. O MTL apresentou uma fraca eficiência para redes mais simples, sendo até três vezes mais lento que o MG, porém em sistemas mais complexos, foi o método que melhor se comportou. O MNR revelou-se um algoritmo muito instável para analisar redes mais complexas, produzindo o pior resultado dentre as metodologias avaliadas, na maioria dos casos.

Para melhorar a eficiência dos métodos na análise de redes com DDP foi proposta uma forma de cálculo da vazão do sistema, considerando a estimativa da iteração anterior. Esta relaxação melhorou a taxa de convergência do método MG e MNR na maioria das configurações, porém, para o MTL e MH não houve melhora. Isto ocorreu porque a linearização do MTL emprega um esquema de atualização de vazões semelhante à relaxação proposta neste trabalho.

AGRADECIMENTOS

Este trabalho teve apoio da FAPESP que subsidiou os estudos do primeiro autor por meio de bolsa de Doutorado, e do CNPq, Processo 478810/2006-7, cujo de auxílio financeiro possibilitou o prosseguimento das investigações. Agradecemos também aos revisores que deram importantes contribuições à versão final deste artigo.

Recebido: 30/04/07

Aceito: 08/04/08

• ACKLEY, J. R. L., TANYIMBOH, T. T., TAHAR, B. Head Driven Analysis of Water Distribution Systems In: B. Ulanicki, B. Boulbeck, et al (Ed.). Water Software Systems: Theory and Applications. London: Research Studies Press, v.1, p.183-192, 2001.
• COLLINS, M., COOPER, L., HELGASON, R. Solving the Pipe Network Analysis Problem Using Optimization Techniques. Management Science, v.24, n.7, p.747-760. 1978.
• CROSS, H. Analysis of Flow in Networks of Conduits or Conductors Enginnering Experimental Station: University of Illinois. (Bulletin Nş 286), 1936.
• FORMIGA, K. T. M. Otimização multiobjetivo de projetos de redes de distribuição de água (Tese de Doutorado). Hidráulica e Saneamento, Universidade de São Paulo, São Carlos, 276 p. 2005.
• FUJIWARA, O. E GANESHARAJAH, T. Reliability assessment of water supply systems with storage and distributions networks. Water Resources Research, v.29, n.8, p.2917-2924. 1993.
• GERMANOPOULOS, G. E JOWITT, P. W. Leakage reduction by excessive pressure minimization in a water supply network. Proceedings of the Institution of Civil Engineers, v.87, n.2, p.195-214. 1989.
• GOMES, H. P. Sistema de Abastecimento de Água - Dimensionamento Econômico e Operação de Redes e Elevatórias João Pessoa: Editora Universitária - UFPB. 242 p. 2004.
• GUPTA, R. E BHAVE, P. R. Comparison of Methods for Predicting Deficient-Network Performance Journal of Water Resources Planning and Management, ASCE, v.122, n.3, p.214-217. 1996.
• HANSEN, C. T., MADSEN, K. e B, N. H. Optimization of pipe networks Mathematical Programming, v.52, n.1, p.45-58. 1991.
• ISAACS, L. T. E MILLS, K. G. Linear Theory Method for Pipe Network Analysis. Journal of Hydraulic Engineering, v.106, n.7, p.1191-2001. 1980.
• JOWITT, P. W. e Xu, C. Optimal valve control in water-distribution networks. Journal of Water Resources Planning and Management, ASCE, v.116, n.4, p.455-472. 1990.
• LEMIEUX, P. F. Efficient Algorithm for Distribution Networks. Journal of the Hydraulics Division, v.98, n.11, p.1911-1920. 1972.
• MARTíNEZ, F., CONEJOS, P. E VERCHER, J. Developing an Integrated Model for Water Distribution Systems Considering both Distributed Leakage and Pressure-Dependent Demands. 26th ASCE Water Resources Planning and Management Division Conference. Tempe, AZ: ASCE, 1999.
• NIELSEN, H. B. Methods for Analyzing Pipe Networks. Journal of Hydraulic Engineering, ASCE, v.125, n.2, p.139-157 1989.
• REDDY, L. E ELANGO, K. Analysis of Water Distribution Networks with Head Dependent Outlets. Civil Engineering Systems, v.6, n.3, p.102-110. 1989.
• SALGADO, R., ROJO, J. E ZEPEDA, S. Extended gradient method for fully non-linear head and flow analysis in pipe networks. In: B. Coulbeck (Ed.). Integrated Computer Applications of Water Supply - Methods and Procedures for Systems Simulation and Control. New York: John Wiley & Sons, Inc., v.1, p.49-60 1993.
• SHAMIR, U. E HOWARD, C. D. D. Water Distribution Systems Analysis. Journal of Hydraulics Division, ASCE, v.94, n.1, p.219-234. 1968.
• SOARES, A. K., CHEUNG, P. B., REIS, L. F. R. Avaliação das perdas físicas de um setor da rede de abastecimento de água de Campo Grande-MS via modelo inverso. Revista Brasileira de Engenharia Sanitária e Ambiental, v.9, n.4, p.312-321. 2004.
• SOUZA, R. S. Aspectos computacionais da analise de redes de distribuição de água com componentes hidráulicos em regime permanente (Dissertação de Mestrado). Hidráulica e Saneamento, Universidede de São Paulo, São Carlos, 236 p. 1994.
• TODINI, E. A more realistic to the "extended period simulation" of water distribution networks. In: C. Maksimovic, D. Butler, et al (Ed.). Advances in Water Supply Management. Lisse A. A. Balkema Publishers, p.732, 2003.
• TODINI, E. E PILATI, S. A Gradient Algorithm for ht Analysis of Pipe Networks In: (Ed.) International Conference on Computer Applications for Water Supply and Distribution: Leicester Polytechnic, p.1-20, 1987.
• TUCCIARELLI, T., CRIMINISI, A. E TERMINI, D. Leak Analysis in Pipeline Systems by Means of Optimal Valve Regulation Journal of Hydraulic Engineering, ASCE, v.125, n.3, p.277-285. 1999.
• WAGNER, J. M., SHAMIR, U. E MARKS, D. H. Water Distribution Reliability: Analytical Methods. Journal of Water Resources Planning and Management, ASCE, v.114, n.3, p.276-292. 1988.
• WOOD, D. J. E CHARLES, C. O. A. Hydraulic Network Analysis Using Linear Theory Journal of Hydraulics Division, ASCE, v.98, n.7, p.1157-1170. 1972.

Datas de Publicação

• Publicação nesta coleção
  30 Jun 2008
• Data do Fascículo
  Jun 2008

Histórico

• Aceito
  08 Abr 2008
• Recebido
  30 Abr 2007