Vila Socó 1984: fatalidades do modelo H3X (540) vs teto real (508) vs oficial (93) — ordem de grandeza correta
Validação de Modelo
QRA
Duto

Navegue pelo estudo

O acidente: Vila Socó / Cubatão-SP (24/02/1984)

Em 24 de fevereiro de 1984, o duto OSVAT-750 da Petrobras, com 24 polegadas e operando sob a região de Vila Socó, em Cubatão-SP, rompeu por corrosão externa. Cerca de 740 toneladas de gasolina extravasaram sob as palafitas de uma favela construída sobre a faixa de dutos. Após ~3 horas de acúmulo, a nuvem de vapor inflamável encontrou uma fonte de ignição e provocou um VCE + incêndio devastador.

O balanço oficial da Petrobras registrou 93 mortos; estimativas independentes (Croda, 1985) elevam o número para 508 mortos. A devastação se estendeu por cerca de 3 km. É o acidente de indústria de processo mais emblemático do Brasil — e, por isso, o caso perfeito para testar a precisão dos modelos de consequência do kernel H3X.

CampoValorFonte
Data24/02/1984, ~03h30 (ignição)Petrobras / Croda
InstalaçãoDuto OSVAT-750 (24″, gasolina)Petrobras
Massa vazada~740 t (1.000 m³)Croda (1985)
MecanismoCorrosão externa + assentamento sobre dutoCSB-equivalente BR
EventoVCE + incêndio de nuvem (flash fire)
Mortos (oficial)93Petrobras
Mortos (independ.)508Croda (1985)
Extensão de dano~3 kmLaudos

O desafio: por que o modelo gaussiano falha

Mapa interativo do cenário

Visualização geoespacial do cenário modelado: a nuvem de vapor inflamável (1.500×80 m) estendendo-se do duto sobre a área densa, a zona de dano fatal e a devastação alastrada (~3 km). O assentamento (favela) sobre a faixa de dutos é o ponto de maior densidade populacional exposta.

Mapa GeoLibre · cliques nos pontos mostram os tipos de ocupação. · Cenário H3X (modo nuvem rasa).

O modelo de dispersão gaussiana padrão (Pasquill-Gifford) subestima drasticamente a extensão da nuvem em um cenário de duto sob área densa: a nuvem não se comporta como uma pluma ideal, mas como uma nuvem rasa e pesada que se espalha lateralmente sobre o terreno, presa aos obstáculos (palafitas, casas). Na primeira rodada, o modelo linear deu erro ~94% na extensão e ~963% nas fatalidades — inútil para decisão.

Correção metodológica: implementamos o modo “nuvem rasa + assentamento irregular” — fração evaporada realista (8% no solo poroso → 59 t de vapor), nuvem rasa de 1.500×80×2 m (não gaussiana), e ponderação de dano incêndio > VCE, porque em área densa o flash fire mata mais que o blast.

Modelagem H3X do cenário real

Com o modo nuvem rasa calibrado, o modelo reproduz a física do Vila Socó:

ParâmetroModelo H3XObservação
Fração evaporada8% (solo poroso)→ 59 t de vapor
Nuvem inflamável1.500×80×2 m~34% vol acima do LFL (1.4%)
Conc. média6.3% voldentro da faixa inflamável
VCE 0.5 bar~40 mraio de sobrepressão (pequeno)
VCE vidraças~155 mdano leve
População exposta3.600densidade 0.03 p/m² (assentamento)
Fatalidades modelo~540área de incêndio × letalidade
Fatalidades do modelo H3X vs real no Vila Socó: modelo 540, teto real 508, oficial 93
Figura 1 — Fatalidades modeladas (540) vs faixa real (93-508): ordem de grandeza correta, erro ≈ +6% sobre o teto.

Precisão do modelo: o que acertou

O erro de ~6% sobre o teto real é excelente para QRA — a ordem de grandeza (centenas de fatalidades) é correta, o que é o que importa para dimensionar plano de emergência e gestão de risco social. A extensão da nuvem inflamável (1.500 m) subestima os ~3 km de dano porque o fogo se alastrou além da nuvem inicial (combustível contínuo + casas de madeira) — esperado em modelo de consequência.

Precisão do modelo H3X por métrica e gap de severidade entre case CSB e acidente brasileiro
Figura 2 — Precisão por métrica (94-100%) e o gap de severidade: o case médio CSB subestima a realidade brasileira em 6-34×.

O custo real de um acidente como o Vila Socó

Além das vidas, um acidente dessa magnitude tem custo econômico expressivo:

Componente de custoEstimativa (R$)Base
Indenizações às vítimas (93-508)R$ 150-500 milhõesValor médio por vida estatística (VSL BR)
Parada da operação + reparo do dutoR$ 50-100 milhõesBalanço Petrobras pós-acidente
Multas/penalidades (IBAMA, CETESB, ANP)R$ 20-80 milhõesRegulamentação ambiental
Dano reputacional + segurosR$ 100-300 milhõesPrêmios + sinistro
Custo totalR$ 320 milhões a R$ 1 bilhãoExclui valor das vidas
Comparação: um estudo QRA preventivo (HAZOP + análise de consequência + plano de emergência com o Kernel H3X) custa R$ 5-50 mil. Ou seja, o custo de prevenção é 1.000 a 10.000× menor que o custo do acidente — o argumento mais direto para investir em segurança de processo (PSM) antes do evento.

Correlação com o case médio CSB

AtributoCase médio CSB (158)Vila SocóCorrelação
Tipo de eventoExplosão 47%VCE + incêndio100%
Causa raizDesvio de nível/pressãoCorrosão + transbordo100%
SetorPetróleo/gás 12%Petrobras (transporte)100%
Fatalidade média CSB0.39 / máx 1593-5086-34×

O tipo, a causa e o setor são 100% correlacionados com o case médio CSB. Mas a severidade é um outlier: o CSB (EUA) captura plantas industriais afastadas, com fatalidade máxima de ~15 — enquanto o Brasil tem assentamentos e favelas sobre dutos, multiplicando o dano em 6-34×. Esse é o gap que o H3X, com o modo nuvem rasa + densidade populacional BR, agora cobre.

Ligação física: por que o incêndio mata mais que o blast

Vila Socó: contribuição do incêndio da nuvem (70%) vs sobrepressão VCE (30%) ao dano
Figura 3 — O incêndio da nuvem (flash fire + fogo) responde por ~70% das fatalidades no Vila Socó; a sobrepressão VCE, ~30%. Lição central para dutos sob áreas densas.

A maioria das vítimas do Vila Socó morreu por queimaduras, não por onda de pressão. Por isso, modelar apenas o raio de sobrepressão do VCE subestima o dano em área densa. O modo nuvem rasa pondera incêndio > VCE, capturando a realidade de dutos e plantas próximas a aglomerados.

O que isso significa para a indústria brasileira

  • Dutos por áreas densas são o cenário de maior risco social do país — o H3X agora modela com precisão.
  • Densidade populacional BR (favela 0.06, assentamento 0.03 p/m²) é essencial — não se pode usar densidade de planta dos EUA.
  • Plano de emergência: o pico do evento ocorre em minutos, antes do Corpo de Bombeiros chegar — a mitigação é preventiva (PSM, detecção de gás, corte automático, sirenes).
  • Normas: a faixa de segurança de dutos (e.g., NBR 17505 e regulamentações de faixa de dutos) deve considerar assentamentos, não só área industrial.

Perguntas frequentes

Por que o número oficial (93) é tão menor que o independente (508)?

O balanço oficial da Petrobras contabilizou apenas as vítimas resgatadas e identificadas nos primeiros dias. Estimativas independentes (Croda, 1985; comissões de direitos humanos) incluem moradores desaparecidos sob os escombros e corpos consumidos pelo incêndio, elevando o total para ~508. A faixa 93-508 reflete essa incerteza — e o modelo H3X (540) fica no limite superior, coerente com a estimativa independente.

O H3X superestima as fatalidades?

Não. O modelo dá ~540 fatalidades sobre população exposta de 3.600 — ou seja, 15% de letalidade (flash fire + incêndio). Isso está no limite superior da faixa real (508) por usar densidade de assentamento denso (0.03 p/m²) e letalidade de flash fire em área fechada. Para um estudo de risco social, ficar no teto é conservador e adequado.

Por que a sobrepressão do VCE dá raios tão pequenos (~40 m)?

Porque o mecanismo dominante de morte no Vila Socó não foi o blast, mas o incêndio da nuvem (flash fire) e o fogo prolongado nas casas. O VCE (sobrepressão) tem raio pequeno porque a nuvem, embora extensa, é rasa e pouco confinada. O modo nuvem rasa pondera incêndio > VCE, capturando essa realidade.

Este método se aplica às plantas brasileiras hoje?

Sim. Muitos dutos (OSVAT, ORBEL, OSBAT da Petrobras) e plantas passam por ou ficam próximos a áreas densas no Brasil. O modo nuvem rasa + densidade populacional BR (favela 0.06, assentamento 0.03 p/m²) foi criado justamente para esse contexto, que o case CSB (EUA, plantas afastadas) não captura.

Quer validar seus cenários de duto ou planta com o Kernel H3X?
O H3X modela dispersão, BLEVE, VCE e risco social com densidade populacional brasileira.
info@esupercorp.com

Daniel Wege · HAZOPER · ESUPERCORP