Área com cana e parte do etanol no que os postos vendem, nos municípios de São Paulo

Português · English

Prancha Nº 111, a primeira de cinco da coleção Bioenergia cruzada; a seguinte, a prancha 112.

A coleção cruza a bioenergia paulista com o que se vende nos postos, com o campo e com as prefeituras, nos 645 municípios de São Paulo. Os pares saíram de uma triagem de correlações entre o potencial de biometano do CP2B, as vendas de combustíveis do Anuário da SEMIL, a silvicultura, o rebanho, o Censo Agropecuário, a frota e a MUNIC. Cada par passou por correlação de postos bruta e parcial (descontadas urbanização, densidade e renda), amostra efetiva corrigida pela autocorrelação espacial, valores q de Benjamini–Hochberg e LISA bivariado. A última prancha é descritiva, sem teste, e o texto diz isso.

As chamadas do mapa apontam Ibaté, onde o etanol é 69 % da energia vendida no ciclo Otto (a maior parte entre os municípios de 10 mil habitantes ou mais) e a cana cobre 51 % da área; São Luiz do Paraitinga, sem cana e com 16 % (a menor); e São Paulo, com 39 %, pouco abaixo do estado.

O que vem nos arquivos

Camada (GPKG)	Geometria	Conteúdo
municipios	polígono	645 municípios com cana e etanol, biometano e diesel, frota, mata nativa e madeira para celulose, bovinos e adubo orgânico, distâncias ao gasoduto, gás canalizado e as respostas da MUNIC sobre resíduos; classes b111, b113, b114, arranjo115 e aglomerados lisa111, lisa113, lisa114
gasodutos	linha	gasodutos de transporte e rede de distribuição (tipo)
regioes_intermediarias, limite_estadual	polígono	referência

O GPKG é o mesmo nas cinco pranchas e traz um estilo do QGIS para cada uma. O CSV tem a tabela dos municípios, sem geometria. Campos vazios são falta de dado, não zero. Os aglomerados LISA vêm como AA (alto-alto), BB (baixo-baixo), AB, BA e ns (não significativo); a primeira letra é a variável do município, a segunda a dos vizinhos.

Método

1. pct_area_cana: área com cana-de-açúcar (MapBiomas, coleção 10, 2024) sobre a área do município, em %.
2. etanol_l e gasolina_l: litros de etanol hidratado e de gasolina C vendidos no município em 2024 (SEMIL/SP, Anuário de Energéticos por Município 2025). pct_etanol_otto = etanol × 21,3 MJ/L ÷ (etanol × 21,3 + gasolina C × 32,2 MJ/L), em %; a gasolina C já traz etanol anidro misturado, que fica na parte da gasolina.
3. Terços de cana: até 3,4 %, de 3,4 % a 35,6 % e mais de 35,6 %; de etanol: até 39,8 %, de 39,8 % a 52,5 % e mais de 52,5 %. b111: classe bivariada, dezena = terço de cana (1 a 3), unidade = terço de etanol (1 a 3). De baixo para cima e da esquerda para a direita, as nove casas têm 169, 39, 6; 38, 102, 74; 7, 73, 135 municípios. Anhumas e Redenção da Serra não têm venda de derivados no Anuário e ficam de fora (643 municípios).
4. Cana × etanol no ciclo Otto: correlação de postos de Spearman 0,70 (643 municípios), parcial 0,59, amostra efetiva 310, p < 0,001, q < 0,001. Sem a capital, 0,59; fora da região intermediária de São Paulo, parcial 0,56.
5. lisa111: aglomerados LISA bivariados (cana no município × etanol nos vizinhos): 179 alto-alto, 176 baixo-baixo, 2 alto-baixo e 8 baixo-alto; I de Moran bivariado 0,68.
6. Correlação parcial: correlação entre os resíduos dos postos de cada variável regredidos nos postos dos controles, que são urbanização (parte de moradores em situação urbana), log da densidade demográfica e log da renda por habitante (Censo 2022). Amostra efetiva sob autocorrelação espacial: n·(1 − Iₓ·Iᵧ)/(1 + Iₓ·Iᵧ), com Iₓ e Iᵧ o I de Moran dos resíduos numa vizinhança dos 6 vizinhos mais próximos; valores q de Benjamini–Hochberg sobre todos os pares da triagem. LISA bivariado: valor do município contra a média dos 6 vizinhos na outra variável, 999 permutações, semente 20261019, p < 0,05.

Limitações

• Venda não é consumo: o Anuário da SEMIL conta o que postos e distribuidoras vendem no município, inclusive a quem só passa pela estrada ou abastece frotas de fora. Municípios com rodovia movimentada podem vender bem mais diesel do que gastam.
• O preço relativo na bomba, a frota flex e a distância da usina mais próxima podem estar por trás da relação; o mapa não separa esses efeitos. A parte do etanol é medida em energia, não em litros: em litros, ela seria maior.
• Leitura ecológica: o mapa compara municípios, não pessoas nem fazendas. Uma correlação entre municípios não diz que as mesmas propriedades ou os mesmos motoristas têm as duas características, nem que uma causa a outra.