Índice

Introdução

O uso de excesso de calor para aquecimento urbano.

Entradas e saídas

Camadas e parâmetros de entrada

Fornecido por Toolbox

  • Áreas de aquecimento urbano (por enquanto fornecidas diretamente pelo potencial de aquecimento urbano CM)
  • Banco de dados industrial (por padrão, fornecido pela caixa de ferramentas)

Fornecido pelo usuário

  • Min. demanda de calor em hectare. Consulte DH Potencial CM
  • Min. demanda de calor em uma área de DH. Consulte DH Potencial CM
  • Raio de pesquisa em km
  • Vida útil do equipamento em anos
  • Taxa de desconto em%
  • Fator de custo
  • Custos operacionais em%
  • Valor limite para linhas de transmissão em ct / kWh

Parâmetros de desempenho

  • Resolução de tempo (hora, dia, semana, mês, ano)
  • Resolução espacial em km

Camadas e indicadores de saída

  • Linhas de transmissão
  • Excesso total de calor na área selecionada em GWh
  • Excesso de calor conectado em GWh
  • Excesso de calor usado em GWh
  • Custo da rede em €
  • Custos anuais da rede em € / ano
  • Custos nivelados do suprimento de calor em ct / kWh
  • Gráfico mostrando o potencial DH, excesso de calor total, excesso de calor conectado e excesso de calor usado
  • Gráfico mostrando a demanda mensal de calor e excesso
  • Gráfico mostrando a demanda diária média de calor e excesso

Método

visão global

O elemento chave do módulo de excesso de calor é o modelo de coletor de origem usado. Constrói uma rede de transmissão de comprimento mínimo e calcula o fluxo para cada hora do ano com base em perfis de carga de aquecimento residencial com resolução Nuts2 e perfis de carga do setor com resolução Nuts0. Com base nos fluxos de pico médios ao longo do ano, é possível calcular os custos de cada linha de transmissão e trocador de calor no lado da fonte e do dissipador.

Detalhes

Modelagem de fontes

Com base no ID Nuts0 e no setor industrial, um perfil de carga resolvido a cada ano é atribuído a todas as fontes.

Modelagem de pias

Com base no módulo de cálculo do potencial de aquecimento urbano, são criados pontos de entrada equidistantes nas áreas coerentes. Dependendo do ID Nuts2 dos pontos de entrada, um perfil de carga é atribuído.

Pesquisa de raio fixo

Dentro de um raio definido, é verificado quais fontes estão dentro do alcance uma da outra, quais pias estão dentro do alcance uma da outra e quais pias estão dentro do alcance das fontes. Isso pode ser representado por um gráfico com fontes e pias formando os vértices e os vértices no intervalo sendo conectados por uma aresta.

Redução ao comprimento mínimo da rede

Uma árvore de abrangência mínima é calculada com a distância das arestas como pesos. Isso resulta em um gráfico mantendo sua conectividade enquanto possui um comprimento total mínimo de arestas. Observe que os pontos de entrada de áreas coerentes são conectados internamente gratuitamente, pois formam sua própria rede de distribuição.

Cálculo de fluxo

O fluxo máximo das fontes para os sumidouros é calculado para cada hora do ano.

Determinação de custos

O pico de fluxo do ano, calculado em média por 3 horas, determina a capacidade necessária para as linhas de transmissão e trocadores de calor. Os custos das linhas de transmissão dependem do comprimento e da capacidade, enquanto os custos dos trocadores de calor são influenciados apenas pela capacidade. No lado da fonte, é assumido um trocador de calor ar / líquido com bomba integrada para a linha de transmissão e, no lado da pia, um trocador de calor líquido / líquido.

Variação da rede

Como o custo e o fluxo de cada linha de transmissão são conhecidos, as linhas com a maior relação custo / fluxo podem ser removidas e o fluxo recalculado até que seja alcançado o custo por fluxo desejado.

Implementação

Pesquisa de raio fixo

Para o cálculo da distância entre dois pontos, é utilizada uma pequena aproximação do ângulo do comprimento do loxódromo. Embora exista também uma implementação precisa da distância do ortodromo, a maior precisão não tem nenhum benefício real, devido às pequenas distâncias em sua maioria inferiores a 20 km e à incerteza do comprimento real da linha de transmissão devido a muitos fatores, como topologia. Se dois pontos estiverem dentro do alcance do raio, ele será armazenado em uma lista de adjacências. A criação dessas listas de adjacência é realizada entre fontes e fontes, sumidouros e sumidouros e fontes e sumidouros. O motivo da separação está na flexibilidade de adicionar certos requisitos de temperatura para fontes ou sumidouros.


Exemplo de uma pesquisa de raio fixo. Os vértices vermelhos representam fontes e os azuis afundam. Os números representam a distância entre os pontos. Desenho não está em escala.

Classe NetworkGraph

Com base na biblioteca igraph, uma classe NetworkGraph é implementada com todas as funcionalidades necessárias para o módulo de cálculo. Embora o igraph esteja mal documentado, oferece um desempenho muito melhor do que os módulos python puros como o NetworkX e um suporte mais amplo à plataforma além do Linux, diferentemente da ferramenta gráfica. A classe NetworkGraph descreve apenas uma rede na superfície, mas contém 3 gráficos diferentes. Primeiramente, o gráfico que descreve a rede como é definido pelas três listas de adjacência. Em segundo lugar, o gráfico de correspondência que conecta internamente pias da mesma área coerente e dura o gráfico de fluxo máximo usado para o cálculo do fluxo máximo.

Gráfico

Contém apenas as fontes reais e sumidouros como vértices.


Exemplo de um gráfico. Os vértices vermelhos representam fontes e os azuis afundam.
Gráfico de correspondência

Cada coletor precisa de um ID de correspondência, que indique se está conectado internamente por uma rede já existente, como em áreas coerentes. Pias com o mesmo ID de correspondência são conectadas a um novo vértice com arestas com pesos zero. Isso é crucial para o cálculo de uma árvore de abrangência mínima e a razão pela qual o gráfico de correspondência é usado para ela. Esse recurso também é implementado para fontes, mas não usado.


Exemplo de um gráfico de correspondência. Os vértices vermelhos representam fontes e os azuis afundam. As três pias à direita são coerentes conectadas por um vértice maior adicional
Gráfico de fluxo máximo

Como o igraph não suporta múltiplas fontes e afunda em sua função de fluxo máximo, é necessário um gráfico auxiliar. Introduz uma fonte infinita e um vértice de coletor. Toda fonte real é conectada à fonte infinita e todo coletor real é conectado ao coletor infinito por uma borda. Observe que se um coletor estiver conectado a um vértice de correspondência, esse vértice será conectado em vez do próprio coletor.


Exemplo de um gráfico de fluxo máximo.
Cálculo mínimo da árvore de expansão

Com base no gráfico de correspondência, a árvore de abrangência mínima é calculada. As arestas que conectam os dissipadores coerentes sempre têm o peso 0, portanto sempre permanecerão parte da árvore de abrangência mínima.


Exemplo de um gráfico de correspondência com os pesos de cada aresta e sua árvore de abrangência mínima.
Cálculo do fluxo máximo

O fluxo através das bordas que conectam as fontes reais ou sumidouros à fonte ou coletor infinito, respectivamente, é limitado à capacidade real de cada fonte ou coletor. Por razões numéricas, as capacidades são normalizadas, de modo que a maior capacidade seja 1. O fluxo através do subconjunto de arestas contido no gráfico de correspondência é limitado a 1000, o que deve, para todos os propósitos intensos, oferecer fluxo irrestrito. Em seguida, o fluxo máximo da fonte infinita para o coletor infinito é calculado e o fluxo redimensionado para seu tamanho original. Como os coletores coerentes não estão diretamente conectados ao vértice de coletor infinito, mas pelo vértice de correspondência, o fluxo através dele é limitado à soma de todos os coletores coerentes.


Exemplo de um gráfico de fluxo máximo e as capacidades de cada fonte e coletor. O gráfico à direita mostra o fluxo máximo permitido através de cada borda após a normalização. Observe que o fluxo máximo permitido pelas bordas com o símbolo de infinito é realmente limitado a 1000 na implementação.

A implementação da função de fluxo máximo do igraph usa o algoritmo Push-re-label. Esse tipo de algoritmo não é sensível ao custo e nem sempre pode encontrar a maneira mais curta de rotear o fluxo. Um algoritmo sensível ao custo não está disponível no igraph e o desempenho provavelmente será baixo para poder resolver um fluxo com base horária ao longo do ano. Porém, devido à redução prévia a uma extensão mínima, os casos em que uma solução não ideal é escolhida são muito limitados e improváveis. O algoritmo Push-re-label também tem tendência para rotear o fluxo pela menor quantidade de arestas. A implementação do igraph parece ser determinística na ordem de alocação do fluxo se os gráficos são pelo menos automorfismos, o que é importante para o cálculo de fluxo baseado em uma hora, uma vez que qualquer oscilação de fluxo introduzida artificialmente entre as bordas é indesejável.


Fluxo calculado pelo algoritmo de fluxo máximo e redimensionamento para o tamanho original.

Fontes de calor

As fontes de calor são retiradas do banco de dados industrial. Com base no excesso de calor, o Nuts0 ID e o setor industrial criam um perfil de carga que cobre cada hora do ano para cada local. A adição personalizada de sites está planejada.

Dissipadores de calor

Os dissipadores de calor são baseados em áreas coerentes com uma demanda de calor conhecida. As áreas coerentes formam uma máscara para uma grade na qual pontos equidistantes são colocados como pontos de entrada. Dependendo do ID Nuts2 selecionado, um perfil de aquecimento residencial é atribuído às pias. A adição personalizada de pontos de entrada e pias é planejada.


Exemplo de uma área coerente e seus pontos de entrada gerados.

Carregar perfis

Os perfis de carga mencionados consistem em 8760 pontos que representam a carga para cada hora dos 365 dias. Mais informações sobre os perfis de carga podem ser encontradas aqui.

Cálculo de custos

Como os sistemas de aquecimento urbano têm uma grande capacidade de aquecimento, um pico no fluxo não significa que as linhas de transmissão precisam fornecer instantaneamente esse curto pico de calor. Portanto, as capacidades exigidas das linhas de transmissão e trocadores de calor são determinadas pelo pico de carga médio. Especificamente, a função de convolução numpy é usada para calcular a média do fluxo nas últimas três horas, convoluindo com uma função constante. Dependendo desse valor, uma linha de transmissão da tabela a seguir é escolhida.

Custos específicos das linhas de transmissão utilizadas

| Potência em MW | Custos em € / m | Temperatura em ° C | | ------------- |: -------------: | -----: | | 0,2 | 195 <150 | 0,3 206 <150 | 0,6 | 220 <150 | 1.2 240 <150 | 1,9 261 <150 | 3,6 | 288 <150 | 6.1 323 <150 | 9,8 | 357 <150 | 20 426 <150 | 45 564 <150 | 75 701 <150 | 125 839 <150 | 190 976 <150 | > 190 976 <150

Os custos do trocador de calor no lado da fonte, que é assumido como ar para líquido, são computados com

C HSource (P) = P pico * 15.000 € / MW.

Os custos do trocador de calor líquido para líquido no lado da pia são determinados com

C HSink (P) = P pico * 265.000 € / MW se P pico <1MW ou

C HSink (P) = P pico * 100.000 € / MW mais.

Os custos da bomba seguem

Bomba C (P) = P pico * 240.000 € / MW se P pico <1MW ou

Bomba C (P) = P pico * 90.000 € / MW mais.

Remoção de linhas de transmissão

Com um limite de custo para fluxo para linhas de transmissão, eles podem ser removidos se excederem para melhorar a relação fluxo / custo. Após a remoção das arestas, o fluxo deve ser recalculado, pois a continuidade do fluxo no gráfico não é mais garantida. A relação custo / fluxo também pode aumentar para outras arestas agora, portanto, esse processo é repetido até que a soma de todos os fluxos não seja mais alterada.

Descrição da rotina completa

Primeiro, as fontes de calor e os dissipadores são carregados com seus perfis de carga. Em seguida, a busca por raio fixo é realizada e a rede inicializada. Posteriormente, a rede é reduzida à sua árvore de abrangência mínima e o fluxo máximo é calculado para cada hora do ano. Com base no fluxo, são calculados os custos para cada trocador de calor, bomba e linha de transmissão. Se uma relação de custo-limite for definida, o procedimento de remoção da linha de transmissão é executado. No final, o custo total e o fluxo total da rede e o layout da rede são retornados.

Amostra de execução

Amostra executada em Aalborg.

Amostra executada em Aalborg. As aeras azuis representam o aquecimento urbano. O laranja aponta a fonte de calor e o amarelo aponta os pontos de entrada para a rede de aquecimento urbano.

Os custos totais são de 13,7 M € e o fluxo anual total é de 185 GWh, o que resulta em 0,74 ct / kWh por um período de investimento de 10 anos.

Autores e revisores

Esta página foi escrita por Ali Aydemir * e David Schilling *

  • [] Esta página foi revisada por Tobias Fleiter *.

* Fraunhofer ISI Fraunhofer ISI, Breslauer Str. 48, 76139 Karlsruhe

Licença

Copyright © 2016-2018: Ali Aydemir, David Schilling

Licença Internacional Creative Commons Attribution 4.0 Este trabalho está licenciado sob uma Licença Internacional Creative Commons CC BY 4.0.

Identificador de licença SPDX: CC-BY-4.0

Texto da licença: https://spdx.org/licenses/CC-BY-4.0.html

Reconhecimento

Gostaríamos de expressar nossa profunda gratidão ao Projeto Horizonte Hotmaps 2020 (Contrato de Doação número 723677), que forneceu o financiamento para a realização da presente investigação.

To Top

This page was automatically translated. View in another language:

English (original) Bulgarian* Croatian* Czech* Danish* Dutch* Estonian* Finnish* French* German* Greek* Hungarian* Irish* Italian* Latvian* Lithuanian* Maltese* Polish* Romanian* Slovak* Slovenian* Spanish* Swedish*

*: machine translated