The coupled fluid-flow and stress model with time-dependent effect is developed. In this model, rock mass is simulated as a homogeneous and isotropic dual-porosity, dual-permeability continuum. Darcy's law is used to describe the fluid flow in a porous medium, and cubic law is used to describe the flow in fractures. The finite element method is introduced to solve the system, and the effects of fractures and fracture spacing are considered numerically. The numerical results indicate that fractures have a significant impact on the rocks' displacement and pore pressure. It also shows an increase in plastic strain with decreasing fracture spacing, as does the creep strain. The coupled process is highly sensitive to the fracture spacing. A series of numerical simulations are conducted to better understand the complex coupled processes, which leads to an improved understanding of different aspects of naturally fractured reservoirs and may impact on experimental design to explore these attributes in a real reservoir situation. on modeling fractured formations with the dual-porosity, dual-permeability approach [5,[15][16][17][18]. Some researchers followed the conventional fluid flow modeling coupled with geomechanics through stress dependent rock properties and during the development process [19][20][21]. Considering pore deformation, Valliappan and Khalili-Naghadeh [20] derived a set of coupled differential equations governing the behaviour of fissured porous media. Chen and Teufel [19] focused on theoretical developments with emphasis on the identification of the critical coupling and physical interpretations of the parameters involved. Taking advantage of the joint-mechanics theory, Bagheri and Settari [21] developed general, rigorous coupling between the fluid-flow equation and deformation of fractured media. The need for coupled modeling of hydro-mechanical model has been addressed and investigated in various studies whether the analytical and numerical modelling [22][23][24][25]. For example, Zhang et al. [22] used the finite element method to investigate the effects of fracture spaces on the displacement, stress and pressure in dual-porosity media. However, they neglected the matrix-fracture interaction transfer function. Rutqvist [23] proposed a linked multicontinuum and crack tensor approach for modelling of coupled geomechanics, fluid flow and transport in fractured rock. Hu et al.[25] used a numerical manifold method for analyzing fully coupled hydro-mechanical processes in porous masses with discrete fractures. Some even considered the chemo-hydro-mechanical model for the coupled multiphysics of carbon dioxide sequestration [10,11].Although all of these works provide better extension to the existence theory and models, the long-term creep effect on the matrix with coupled fields is seldom considered. In practice, the time-dependent effect has a great impact on the characteristics of rock material and the stability of rock engineering, such as deep underground rock engineering, high slope engineering, dam ba...