O cenário atual da indústria de gás natural brasileira é caracterizado por baixa maturidade e dinamismo de mercado. O comportamento estocástico da demanda por gás, somado volatilidade do preço de mercado do GNL, motiva a utilização de estocagem subterrânea como forma de inserir flexibilidade no suprimento, além de promover proteção contra flutuação no preço. No entanto, a literatura existente carece de uma uma ferramenta analítica mais robusta para apoiar uma análise quantitativa dos benefícios que a atividade UNGS poderia proporcionar à indústria de gás natural. Nesta tese, propomos um modelo de programação dinâmica estocástica para planejamento de longo/médio prazo, a fim de determinar a política ótima de fornecimento juntamente com a possibilidade de armazenamento de gás. Um modelo markoviano caracteriza a demanda termoelétrica, enquanto o preço de GNL é representado por um processo estocástico temporalmente independente. O modelo proposto é eficientemente resolvido usando o algoritmo de programação dinâmica dual estocástica para o estudo de caso brasileiro, considerando dados dos setores de gás e setor elétrico. Para uma escolha exógena, mas significativa, da localização e tamanho do armazenamento subterrâneo, observamos os benefícios operacionais e econômicos da flexibilidade que esta atividade poderia proporcionar. Além disso, comparando os custos de OPEX e CAPEX de investimentos em infraestrutura de armazenamento em campos depletados e cavernas de sal com as economias proporcionadas pelo armazenamento na operação de fornecimento, é possível observar o benefício econômico da atividade de estocagem. A estrutura proposta fornece suporte quantitativo importante para discussões sobre precificação de infraestrutura e modelo de negócios para Armazenamento Subterrâneo de Gás Natural.