“…This fracture modeling strategy is valid if fractures are more permeable than the matrix, that is, it is limited to fractures being channels for flow drainage (Segura & Carol, 2008b,a; Jha & Juanes, 2014; Cueto‐Felgueroso et al., 2017, 2018; Andrés et al., 2021). This approach agrees with proposals of other authors to deal with the fracture network modeling (McClure & Horne, 2013; Garipov & Hui, 2019; Berre et al., 2020), including the full simulation of matrix and fractures with poroelastic coupling. It also allows to vary parameters such as fracture aperture or porosity depending on the stress or strain level, although this barely affects the results (see Appendix ) (Bandis et al., 1983; Chin et al., 2000; Hsiung et al., 2005; Rutqvist et al., 2005; Cappa & Rutqvist, 2011).…”