A teoria da poroelasticidade surge do acoplamento entre o balanço de forças na estrutura sólida de um meio poroso com a mecânica do fluido escoando dentro deste meio. ela possui uma vasta gama de aplicações pelas quais desperta um crescente interesse em diversos pesquisadores. o primeiro modelo matemático para a solução de problemas de poroelasticidade foi proposto por biot em 1941, o qual ainda hoje é utilizado como referência para a elaboração de novos modelos. muitos trabalhos que utilizam o modelo de biot simplificam as variações de porosidade e densidade e desconsideram por completo as variações espaciais destas propriedades. as simplificações desta abordagem reduzem drasticamente a complexidade do problema e os custos computacionais associados à sua solução numérica, no entanto, elas também limitam a variedade de situações nas quais este modelo pode ser aplicado. neste trabalho uma abordagem diferente será proposta, na qual as simplificações previamente mencionadas não serão aplicadas, obtendo assim um modelo mais geral o qual não é limitado pelas hipóteses simplificativas. para a correção dos valores do campo de densidades será utilizado um modelo em função da pressão do fluido, enquanto que para a porosidade será utilizada uma expressão baseada no modelo de biot. ademais, para garantir uma solução numericamente conservativa será utilizado o método de volumes finitos baseado em elementos (ebfvm). os resultados mostram que as hipóteses do modelo linear nem sempre são válidas, sugerindo portanto que utilizar o modelo não linear é uma abordagem mais consistente. adicionalmente, o uso de uma técnica de estabilização para resolver um importante problema de benchmark é também demonstrado neste trabalho.