Este trabalho desenvolveu um modelo matemático numérico baseado na conservação da massa. Este modelo desenvolvido foi concebido incorporando e investigando o impacto do parâmetro governante - permeabilidade aparente do gás, difusão de Knudsen, fator Forchheimer e isoterma de Langmuir - para definir a adsorção de gás dependente de pressão necessária para a transição de fluxo da matriz para o influxo da matriz de fratura e, em seguida, a fratura para o wellbore. O modelo desenvolvido também foi utilizado para analisar o efeito da contribuição do mecanismo de fluxo não-darcy da matriz na produção de gás de xisto. O modelo desenvolvido é uma ferramenta sensível que descreve a complexidade da produção dos reservatórios de gás xistoso e identifica regimes de fluxo para várias condições de produção. As variáveis sem dimensão ou o seu produto definem a produção de gás de xisto e a contribuição do domínio da matriz ou da fractura.