A comercialização de energia oferece grandes oportunidades de ganho, mas também a exposição ao risco devido à volatilidade dos preços. Um modelo de preço aderente ao mercado é um elemento-chave da análise de risco por empresas atuantes no mercado de eletricidade e representa uma forte ferramenta para fornecer um melhor entendimento da dinâmica do mercado. Este trabalho aplica Simulação de Monte-Carlo aliada ao Movimento Browniano Geométrico e a decomposição de Cholesky, procedimentos estocásticos amplamente utilizados nos mercados financeiro e de ações, para simular curvas estocásticas de preços futuros para o mercado brasileiro de energia, com foco no submercado Sudeste/Centro-Oeste, que corresponde à maior parte da comercialização de energia no país. Os resultados deste trabalho representam convenientemente a incerteza nos preços do mercado brasileiro de energia, demonstrando resultados robustos mesmo quando aplicados para os recentes períodos de alta volatilidade que o mercado brasileiro vem passando devido às crises hídricas e econômicas desde o ano 2012 até a recente crise com a pandemia de 2020. Ao fim, um caso prático de otimização de portfólio consolida a utilização das ferramentas estocásticas propostas neste trabalho como dados de entrada para a solução do problema de uma comercializadora de energia.