Na produção de gás natural, é de grande importância que todo e qualquer líquido presente no poço seja carregado pelo gás até o separador. o acúmulo de líquido no poço devido à perda de capacidade do gás em arrastá-lo decorre do próprio envelhecimento do reservatório, mas também da sensibilidade das propriedades dos hidrocarbonetos a variações de pressão e temperatura durante transientes de produção. desta forma, a redução da quantidade de movimento do gás ascendente gera um escoamento oscilatório do líquido que, ao se acumular no poço, aumenta sua queda de pressão, reduzindo assim, ou até mesmo interrompendo, a produção. a fim de abordar o fenômeno acima descrito, denominado carga de líquido (liquid loading), é necessário compreender a interação entre o poço e a região adjacente do reservatório. adicionalmente, é necessário possuir ferramentas capazes de resolver o escoamento bifásico no poço em função tanto do espaço quanto do tempo. embora a literatura seja rica em trabalhos sobre escoamentos gás-líquido em canais verticais, estes resultam, em sua maioria, em correlações empíricas ou modelos simplificados para o cálculo da queda de pressão e das frações volumétricas das fases em condições de regime permanente. mesmo que tais relações sirvam de critério para a ocorrência do fenômeno de liquid loading, elas não possibilitam a descrição da sequência de eventos que levam à transição entre os regimes de escoamento unidirecional ascendente e oscilatório. este trabalho apresenta um modelo diferencial unidimensional para o cálculo do escoamento gás-líquido transiente em tubos verticais com altas frações de gás (padrões anular e churn). equações de conservação de natureza hiperbólica para a massa, quantidade de movimento e energia são propostas para as fases gasosa e líquida, a qual é dividida entre os campos de filme contínuo e gotículas arrastadas pelo núcleo de gás. relações de fechamento para o cálculo do atrito na parede e na interface do filme, além de expressões para as taxas de entranhamento e deposição de gotículas, foram obtidas da literatura. um algoritmo de solução por diferenças finitas baseado no método da separação da matriz de coeficientes foi implementado a fim de melhor lidar com variações acentuadas no domínio espaço-temporal, tais como ondas de pressão e de retenção de fases. os resultados do modelo foram comparados com dados experimentais de regime permanente obtidos de oito referências da literatura, somando mais de 1300 pontos para o gradiente de pressão, a vazão do filme de líquido e a retenção gasosa. para estas variáveis, a concordância entre o modelo e os dados foi em margens inferiores a ± 20%. o modelo foi avaliado também contra dados experimentais para escoamentos gás-líquido transientes em um tubo vertical de 42 m de comprimento e diâmetro de 0,049 m (waltrich, 2012). transientes induzidos por variações de pressão e de vazão de líquido na seção de testes foram simulados, conferindo níveis de concordância entre dados experimentais e modelo matemático também inferiores a ± 20%.