Se desarrolla el modelo acoplado de flujo de fluidos y estrés con efecto dependiente del tiempo. En este modelo, la masa rocosa se simula como un continuo de doble porosidad y doble permeabilidad homogéneo e isotrópico. La ley de Darcy se utiliza para describir el flujo de fluidos en un medio poroso, y la ley cúbica se utiliza para describir el flujo en fracturas. Se introduce el método de elementos finitos para resolver el sistema, y se consideran numéricamente los efectos de las fracturas y el espaciamiento de las fracturas. Los resultados numéricos indican que las fracturas tienen un impacto significativo en el desplazamiento de las rocas y la presión de poro. También muestra un aumento en la deformación plástica con la disminución del espaciamiento de las fracturas, al igual que la deformación por fluencia. El proceso acoplado es altamente sensible al espaciamiento de las fracturas. Se realizan una serie de simulaciones numéricas para comprender mejor los complejos procesos acoplados, lo que lleva a una mejor comprensión de diferentes aspectos de los reservorios naturalmente fracturados y puede impactar en el diseño experimental para explorar estos atributos en una situación de reservorio real.