A field study (Massachusetts, USA) of the factors controlling the depth of groundwater flow systems in crystalline fractured-rock terrain

Groundwater movement and availability in crystalline and metamorphosed rocks is dominated by the secondary porosity generated through fracturing. The distributions of fractures and fracture zones determine permeable pathways and the productivity of these rocks. Controls on how these distributions vary with depth in the shallow subsurface (<300 m) and their resulting influence on groundwater flow is not well understood. The results of a subsurface study in the Nashoba and Avalon terranes of eastern Massachusetts (USA), which is a region experiencing expanded use of the fractured bedrock as a potable-supply aquifer, are presented. The study logged the distribution of fractures in 17 boreholes, identified flowing fractures, and hydraulically characterized the rock mass intersecting the boreholes. Of all fractures encountered, 2.5% are hydraulically active. Boreholes show decreasing fracture frequency up to 300 m depth, with hydraulically active fractures showing a similar trend; this restricts topographically driven flow. Borehole temperature profiles corroborate this, with minimal hydrologically altered flow observed in the profiles below 100 m. Results from this study suggest that active flow systems in these geologic settings are shallow and that fracture permeability outside of the influence of large-scale structures will follow a decreasing trend with depth.     

Mechanical and hydraulic properties of rocks related to induced seismicity

Witherspoon, P.A. and Gale, J.E., 1977. Mechanical and hydraulic properties of rocks related to induced seismicity. Eng. Geol., 11(1): 23–55.The mechanical and hydraulic properties of fractured rocks are considered with regard to the role they play in induced seismicity. In many cases, the mechanical properties of fractures determine the stability of a rock mass. The problems of sampling and testing these rock discontinuities and interpreting their non-linear behavior are reviewed. Stick slip has been proposed as the failure mechanism in earthquake events. Because of the complex interactions that are inherent in the mechanical behavior of fractured rocks, there seems to be no simple way to combine the deformation characteristics of several sets of fractures when there are significant perturbations of existing conditions. Thus, the more important fractures must be treated as individual components in the rock mass.In considering the hydraulic properties, it has been customary to treat a fracture as a parallel-plate conduit and a number of mathematical models of fracture systems have adopted this approach. Non-steady flow in fractured systems has usually been based on a two-porosity model, which assumes the primary (intergranular) porosity contributes only to storage and the secondary (fracture) porosity contributes only to the overall conductivity. Using such a model, it has been found that the time required to achieve quasi-steady state flow in a fractured reservoir is one or two orders of magnitude greater than it is in a homogeneous system. In essentially all of this work, the assumption has generally been made that the fractures are rigid.However, it is clear from a review of the mechanical and hydraulic properties that not only are fractures easily deformed but they constitute the main flow paths in many rock masses. This means that one must consider the interaction of mechanical and hydraulic effects. A considerable amount of laboratory and field data is now available that clearly demonstrates this stress-flow behavior. Two approaches have been used in attempting to numerically model such behavior: (1) continuum models, and (2) discrete models. The continuum approach only needs information as to average values of fracture spacing and material properties. But because of the inherent complexity of fractured rock masses and the corresponding decrease in symmetry, it is difficult to develop an equivalent continuum that will simulate the behavior of the entire system. The discrete approach, on the other hand, requires details of the fracture geometry and material properties of both fractures and rock matrix. The difficulty in obtaining such information has been considered a serious limitation of discrete models, but improved borehole techniques can enable one to obtain the necessary data, at least in shallow systems. The possibility of extending these methods to deeper fracture systems needs more investigation. Such data must be considered when deciding whether to use a continuum or discrete model to represent the interaction of rock and fluid forces in a fractured rock system, especially with regard to the problem of induced seismicity. When one is attempting to alter the pressure distribution in a fault zone by injection or withdrawal of fluids, the extent to which this can be achieved will be controlled in large measure by the behavior of the fractures that communicate with the borehole. Since this is essentially a point phenomenon, i.e., the changes will propagate from a relatively small region around the borehole, the use of a discrete model would appear to be preferable.     

Theoretical and numerical study of the steady‐state flow through finite fractured porous media

This paper investigates the two‐dimensional flow problem through an anisotropic porous medium containing several intersecting curved fractures. First, the governing equations of steady‐state fluid flow in a fractured porous body are summarized. The flow follows Darcy's law in matrix and Poiseuille's law in fractures. An infinite transversal permeability is considered for the fractures. A multi‐region boundary element method is used to derive a general pressure solution as a function of discharge through the fractures and the pressure and the normal flux on the domain boundary. The obtained solution fully accounts for the interaction and the intersection between fractures. A numerical procedure based on collocation method is presented to compute the unknowns on the boundaries and on the fractures. The numerical solution is validated by comparing with finite element solution or the results obtained for an infinite matrix. Pressure fields in the matrix are illustrated for domains containing several interconnected fractures, and mass balance at the intersection points is also checked. Copyright © 2013 John Wiley & Sons, Ltd.     

Effects of mechanical layering on hydraulic fracturing in shale gas reservoirs based on numerical models

Generally, induced hydraulic fractures are generated by fluid overpressure and are used to increase reservoir permeability through forming interconnected fracture systems. However, in heterogeneous and anisotropic rocks, many hydraulic fractures may become arrested or offset at layer contacts under certain conditions and do not form vertically connected fracture networks. Mechanical layering is an important factor causing anisotropy in sedimentary layers. Hence, in this study, with a shale gas reservoir case study in the Longmaxi Formation in the southeastern Chongqing region, Sichuan Basin, we present results from several numerical models to gain quantitative insights into the effects of mechanical layering on hydraulic fracturing. Results showed that the fractured area caused by hydraulic fracturing indicated a linear relationship with the neighboring layer’s Young’s modulus. An increase of the neighboring layer’s Young’s modulus resulted in better hydraulic fracturing effects. In addition, the contact between two neighboring layers is regarded as a zone with thickness and mechanical properties, which also influences the effects of hydraulic fracturing in reservoirs. The initial hydraulic fracture was unable to propagate into neighboring layers under a relatively low contact’s Young’s modulus. When associated local tensile stresses exceeded the rock strength, hydraulic fractures propagated into neighboring layers. Moreover, with the contact’s Young’s modulus becoming higher, the fractured area increased rapidly first, then slowly and finally became stable.     

Finite-difference modelling of anisotropic wave scattering in discrete fracture model

The presence of fractures in reservoir rocks causes scattering of seismic wave energy. In this paper, we utilize the finite-difference modelling technique to study these scattering effects to gain more insights into the effects and assess the validity of using anisotropic wave scattering energy as a diagnostic tool to characterize fractured hydrocarbon reservoirs. We use a simplified fractured reservoir model with four horizontal layers with a fractured-layer as the third layer. The fractures are represented by grid cells containing equivalent anisotropic medium by the use of the linear slip equivalent model. Our results show that the scattered energy, quantified through estimates of the seismic quality factor (Q) is anisotropic, exhibiting a characteristic elliptical (\(\cos 2\theta \)) variations relative to the survey azimuth angle \(\theta \). The fracture normal is inferred from the minor axis of the Q ellipse. This direction correlates with the direction of maximum wave scattering. Minimum wave scattering occurs in the fracture strike direction inferred from the major axis of the Q ellipse. These results provide more complete insights into anisotropic wave scattering characteristics in fractured media and thus, validate the practical utility of using anisotropic attenuation attribute as an additional diagnostic tool for delineation of fracture properties from seismic data.     

基于零厚度内聚力单元的水力裂缝随机扩展方法研究

为有效模拟裂缝性页岩储层中水力裂缝随机扩展过程,基于单元节点的拓扑数据结构,利用网格节点分裂方式,建立了一种基于有限元网格嵌入零厚度内聚力单元的水力裂缝随机扩展新方法。利用KGD模型解析解和2种室内试验,验证了新方法的准确性和有效性。同时,通过数值算例研究了水平地应力差和储层非均质性对水力裂缝随机扩展过程的影响。研究表明:(1)该方法弥补了ABAQUS平台内置的内聚力单元无法有效模拟水力裂缝随机扩展的不足;(2)在较高水平地应力差下页岩储层非均质性越强,与水力裂缝相交的高角度天然裂缝越容易开启。所建方法能准确地描述复杂水力裂缝的随机扩展行为,可为裂缝性页岩储层的数值模拟提供新手段。     

离散裂隙渗流方法与裂隙化渗透介质建模

流体渗流模拟的连续介质方法通常适用于多孔地质体，并不一定适用于裂隙岩体，由于裂隙分布及其特征与孔隙差异较大。若流体渗流主要受裂隙的控制，对于一定尺寸的裂隙岩体，多孔介质假设则较难刻划裂隙岩体的渗流特征。离散裂隙渗流方法不但可直接用于模拟裂隙岩体非均质性和各向异性等渗流特征，而且可用其确定所研究的裂隙岩体典型单元体及其水力传导（渗透）张量大小。主要讨论了以下问题：（1）饱和裂隙介质中一般的离散流体渗流模拟；（2）裂隙岩体中的REV（典型单元体）及其水力传导（渗透）张量的确定；（3）利用离散裂隙网络流体渗流模型研究裂隙方向几何参数对水力传导系数和REV的影响；（4）在二维和三维离散裂隙流体渗流模型中对区域大裂隙和局部小裂隙的处理方法。调查结果显示离散裂隙流体渗流数学模型可用来评价不同尺度上的裂隙岩体的水力特征，以及裂隙方向对裂隙化岩体的水力特征有着不可忽视的影响。同时，局部小裂隙、区域大裂隙应当区别对待，以便据其所起的作用及水力特征，建立裂隙化岩体相应的流体渗流模型。     

裂隙岩体水力特征数值模拟试验的初步结果

本文利用文献中的数据进行数值模拟试验,结果表明:流量等效和渗透作用力变形效果等效互不包容;不同级别结构面体系难以用连续介质方法处理,同级别结构面体系在工程作用力变化范围内,可由其静态水力特征推演动态特征;工程岩体结构面间充分连通的假设过于粗糙,常规试验手段获取的量值倾向于局部特征的反映。     

渗流应力耦合下含裂缝页岩力学特性和破坏模式研究

为研究在渗流应力耦合作用下含裂缝页岩的力学性质以及破坏模式,利用岩石破坏过程分析系统RFPA2D-Flow,建立7组含不同角度预制裂缝的页岩模型,并对其破裂过程进行了数值模拟。研究表明：页岩的破裂过程分为线性变形、屈服变形以及完全破坏三个阶段,源于裂缝的存在,页岩的弹性模量和抗压强度均表现出明显的各向异性的特点,且试样的最终破裂模式可分为斜线型、X型、λ型与崩坏型四种。对于不同层理倾角下的裂纹分布,渗流—应力耦合作用下页岩试样的破坏主要是拉伸破坏,并伴随着剪切破坏,其中试验表明,当α=90°时,试样所受的剪切破坏最严重。该研究结果对页岩开采具有重大的意义。     

Two-dimensional patterns of metamorphic fluid flow and isotopic resetting in layered and fractured rocks'pa

One-dimensional advection-dispersion models predict that characteristic δ¹⁸O vs. distance and δ¹⁸O vs. δ¹³C profiles should be produced during isothermal metamorphic fluid flow under equilibrium conditions. However, the patterns of isotopic resetting in rocks that have experienced fluid flow are often different from the predictions. Two-dimensional advection-dispersion simulations in systems with simple geometries suggest that such differences may be as a result of fluid channelling and need not indicate disequilibrium, high dispersivities, or polythermal flow. The patterns of isotopic resetting are a function of: (1) the permeability contrast between more permeable layers ('channels') and less permeable layers ('matrix'); (2) the width and spacing of the channels; (3) the width and spacing of discrete fractures; and (4) the orientation of the pressure gradient with respect to layering. In fractured systems, the efficiency of isotopic transport depends on the fracture aperture and the permeability of the surrounding rock. Resetting initially occurs along and immediately adjacent to the fractures, but with time isotopic resetting because of flow through the rock as a whole increases in importance. Application of the one-dimensional advection-dispersion equations to metamorphic fluid flow systems may yield incorrect estimates of fluid fluxes, intrinsic permeabilities, dispersivities, and permeability contrasts unless fluid flow occurred through zones of high permeability that were separated by relatively impermeable layers.     

Development of a full 3D numerical model to investigate the hydraulic fracture propagation under the impact of orthogonal natural fractures

Although hydraulic fracturing has been massively studied and applied as a key technique to enhance the gas production from tight formations, some problems and uncertainties exist to accurately predict and analyze the fracture behavior in complex reservoirs, especially in the naturally fractured reservoirs like shale reservoirs. This paper presents a full 3D numerical model (FLAC3D) to study hydraulic fracturing behavior under the impact of preexisting orthogonal natural fractures. In this numerical model, the hydraulic fracture propagation direction is assumed perpendicular to the minimum principal stress and activated only by tensile failure, whereas the preexisting natural fractures can be activated by tensile or shear failure or a combination of them, and only tensile failure can open the natural fracture as well. The newly developed model was used to study the impact of preexisting orthogonal natural fractures on hydraulic fracturing behavior, based on a multistage hydraulic fracturing operation in a naturally fractured reservoir from the Barnett Shale formation, northwest of Texas in USA. In this multistage operation, two more representative stages, i.e., stage 1 with a relatively large horizontal stress anisotropy of 3.3 MPa and stage 4 with a comparatively small one of 1.3 MPa, were selected to conduct the simulation. Based on the numerical results, one can observe that the interaction between hydraulic and natural fracture is driven mainly by induced stress around fracture tip. Besides, the horizontal stress anisotropy plays a key role in opening the natural fracture. Thus, no significant opened fracture is activated on natural fracture in stage 1, while in stage 4 an opened fracture invades to about 90 m into the first natural fracture. Conversely, the hydraulic fracture length in stage 1 is much longer than in stage 4, as some fluid volume is stored in the opened natural fracture in stage 4. In this work, the shear failure on natural fractures is treated as the main factor for inducing the seismic events. And the simulated seismic events, i.e., shear failure on natural fractures, are very comparable with the measured seismic events.

     

Detection and characterization of hydraulically active fractures in a carbonate aquifer: results from self-potential, temperature and fluid electrical conductivity logging in the Combioula hydrothermal system in the southwestern Swiss Alps

A geophysical and geochemical study has been conducted in a fractured carbonate aquifer located at Combioula in the southwestern Swiss Alps with the objective to detect and characterize hydraulically active fractures along a 260-m-deep borehole. Hydrochemical analyses, borehole diameter, temperature and fluid electrical conductivity logging data were integrated in order to relate electrokinetic self-potential signals to groundwater flow inside the fracture network. The results show a generally good, albeit locally variable correlation of variations of the self-potential signals with variations in temperature, fluid electrical conductivity and borehole diameter. Together with the hydrochemical evidence, which was found to be critical for the interpretation of the self-potential data, these measurements not only made it possible to detect the hydraulically active fractures but also to characterize them as zones of fluid gain or fluid loss. The results complement the available information from the corresponding litholog and illustrate the potential of electrokinetic self-potential signals in conjunction with temperature, fluid electrical conductivity and hydrochemical analyses for the characterization of fractured aquifers, and thus may offer a perspective for an effective quantitative characterization of this increasingly important class of aquifers and geothermal reservoirs.     

Modelling of steady‐state fluid flow in 3D fractured isotropic porous media: application to effective permeability calculation

In this paper, 3D steady‐state fluid flow in a porous medium with a large number of intersecting fractures is derived numerically by using collocation method. Fluid flow in the matrix and fractures is described by Darcy's law and Poiseuille's law, respectively. The recent theoretical development presented a general potential solution to model the steady‐state flow in fractured porous media under a far‐field condition. This solution is a hypersingular integral equation with pressure field in the fracture surfaces as the main unknown. The numerical procedure can resolve the problem for any form of fractures and also takes into account the interactions and the intersection between fractures. Once the pressure field and then the flux field in fractures have been determined, the pressure field in the porous matrix is computed completely. The basic problem of a single fracture is investigated, and a semi‐analytical solution is presented. Using the solution obtained for a single fracture, Mori‐Tanaka and self‐consistent schemes are employed for upscaling the effective permeability of 3D fractured porous media. Copyright © 2012 John Wiley & Sons, Ltd.     

蒲石河抽水蓄能电站上水库裂隙岩体渗透张量的确定

本文提出了一种既能反映裂隙岩体的渗透特性,又相对准确的确定裂隙岩体渗透张量的方法。首先通过裂隙在空间展布状况的测量,用统计学方法初步确定裂隙岩体的渗透张量,获得渗透主值及主方向,然后根据野外压水试验得到的岩体透水率,利用巴布什金公式计算各试段岩体的渗透系数,求出修正系数,从而得到研究区裂隙岩体的修正渗透张量。并运用上述方法对蒲石河抽水蓄能电站上水库坝址区裂隙岩体的渗透张量进行了计算。结果表明,该方法能较好地刻画裂隙岩体渗透性的各向异性特征,可为岩体渗透性分区及防渗帷幕的优化提供科学依据。     

Can tensile stress develop in fractured multilayers under compressive strain conditions?

To contribute to the understanding of how opening-mode fractures (joints) form and open or close at depth in layered rocks, we present a 2D numerical study aiming to determine whether tensile stress can develop in pre-fractured elastic multilayers submitted to biaxial compressive strain conditions.First, we investigate the role of the elastic and geometrical properties of the layers on the development of tensile stress in models with five bonded layers and containing one open fracture in the central layer. Our results indicate that, in absence of elastic contrast (in Young's modulus) between the layers, no tensile stress develops in the models. However, when the fractured layer is stiffer than the two adjacent layers directly above and below, a lobe of horizontal tensile stress develops centered on the pre-existing fracture. The creation of this tensile stress is contingent upon the partial closing of the fracture. The levels of tensile stress and the thickness of the lobe of tensile stress increase logarithmically with an increase in the elastic contrast and are systematically larger for a larger Soft/Stiff ratio (ratio of the total thickness of the soft layers with the total thickness of the stiff layers).Second, we investigate the role of fracture interaction in the development of tensile stress in models containing a pair of open fractures. We observe that the levels of tensile stress in the region between the fractures are systematically higher than those observed in identical models containing a single fracture. This increase in tensile stress is very large for small elastic contrasts between the layers but diminishes when the elastic contrast increases. Furthermore, the spacing between the pre-existing fractures plays an important role in the stress distribution in the region between them. When the fracture spacing is equal to or lower than 1.15 times the height of the fractured layer for the experimental conditions chosen, the lobes of tensile stress centered on the fractures coalesce. This results in the formation of vast areas of tensile stress in models under remote compressive loading conditions. Such tensile areas are likely to allow the initiation and propagation of subsequent opening-mode fractures.The results obtained provide new insights into the formation of joints in layered rocks in compressive environments, with important consequences on fluid flow.     

单轴压缩条件下岩溶化裂隙岩体损伤破坏特征研究

岩溶化裂隙岩体是普遍发育于自然界中具有初始损伤的岩体。为了研究岩溶化裂隙岩体损伤破坏特征,本文以贵州某地赋存的溶蚀岩体为研究对象,运用损伤力学理论构建岩溶化裂隙岩体在单轴压缩条件下的损伤演化模型,并建立岩溶化裂隙岩体损伤演化方程。采用颗粒流数值软件进行单轴压缩数值试验,进一步研究岩溶化裂隙岩体试件在单轴压缩条件下的损伤演化特征,分析岩溶化裂隙岩体的微观损伤特征。结果表明：岩溶化裂隙岩体的初始损伤主要包括溶蚀损伤和裂隙损伤。岩溶化裂隙岩体的初始损伤随着溶蚀率的增加而增加,最终增加速率趋于平缓;岩溶化裂隙岩体的损伤演化曲线均呈“S”型变化,先缓慢增加,再迅速增加,最后缓慢增加至损伤值1;岩溶化裂隙岩体存在异构特征,导致破坏裂隙起源于具有初始损伤的溶蚀孔洞和裂隙处,随后裂隙经历萌发、扩张和剪切作用、数量和长度增加以及裂隙贯通4个阶段后发生宏观破坏。     

基于模拟退火粒子群优化算法的裂缝型储层各向异性参数地震反演

发育垂直定向排列裂缝的地下岩石可等效为具有水平对称轴的横向各向同性（horizontal transverse isotropic,HTI）介质。针对HTI介质模型,本文研究了裂缝型储层的各向异性参数地震振幅随方位角变化（amplitude variations with azimuth,AVAZ）的反演方法。首先,在地震AVAZ反演流程中,提出采用模拟退火粒子群优化算法实现裂缝型储层各向异性参数反演。之后,通过理论模型测试,验证了基于模拟退火粒子群优化算法的地震AVAZ反演的有效性。最后,将反演方法应用于四川盆地龙马溪组页岩气储层实际方位地震数据;在反演之前先利用傅里叶级数方法估计裂缝方位并对实际数据进行方位校正,以提供更准确的输入数据;通过计算得到的P波、S波各向异性参数可用于评价裂缝发育程度。反演结果表明,研究区域构造顶部裂缝较发育,与地质基本理论一致,验证了反演方法的合理性。     

Steeply-dipping extension fractures in the Newark basin, New Jersey

Late Triassic and Early Jurassic bedrock in the Newark basin is pervasively fractured as a result of Mesozoic rifting of the east–central North American continental margin. Tectonic rifting imparted systematic sets of steeply-dipping, en échelon, Mode I, extension fractures in basin strata including ordinary joints and veins. These fractures are arranged in transitional-tensional arrays resembling normal dip-slip shear zones. They contributed to crustal stretching, sagging, and eventual faulting of basin rift deposits. Extension fractures display progressive linkage and spatial clustering that probably controlled incipient fault growth. They cluster into three prominent strike groups correlated to early, intermediate, and late-stage tectonic events reflecting about 50– 60° of counterclockwise rotation of incremental stretching directions. Finite strain analyses show that extension fractures allowed the stretching of basin strata by a few percent, and these fractures impart stratigraphic dips up to a few degrees in directions opposing fracture dips. Fracture groups display three-dimensional spatial variability but consistent geometric relations. Younger fractures locally cut across and terminate against older fractures having more complex vein-cement morphologies and bed-normal folds from stratigraphic compaction. A fourth, youngest group of extension fractures occur sporadically and strike about E–W in obliquely inverted crustal blocks. A geometric analysis of overlapping fracture sets shows how fracture groups result from incremental rotation of an extending tectonic plate, and that old fractures can reactivate with oblique slip components in the contemporary, compressive stress regime.     

各向异性低渗透油藏菱形反九点井网压裂裂缝优化

在低渗透油藏开发中普遍采用菱形反九点井网,但采油井等缝长压裂将导致非均匀驱替,当油藏物性各向异性时,非均匀驱替现象更加明显。以往针对菱形反九点井网的研究主要基于等缝长压裂,未考虑利用不等缝长设计改善开发效果。根据某各向异性低渗透油藏基本参数建立数值模型,进行不等缝长设计,分析油藏渗透性各向异性及油水井压裂情况对菱形反九点井网水驱效果的影响。随着储层Kx∶Ky的增加,角井裂缝优化穿透比减小,而边井裂缝优化穿透比增大。在相等Kx∶Ky条件下,注水井是否压裂对角井裂缝优化穿透比影响较小,对边井裂缝优化穿透比影响较大。数值模拟结果表明,在各向异性低渗透油藏中进行不等缝长优化设计,能有效改善菱形反九点井网的开发效果,优化结果对于各向异性低渗透油藏的水力压裂设计以及井网布置具有一定的指导意义。     

Fully coupled hydro-mechanical numerical manifold modeling of porous rock with dominant fractures

Coupled hydro-mechanical (HM) processes are significant in geological engineering such as oil and gas extraction, geothermal energy, nuclear waste disposal and for the safety assessment of dam foundations and rock slopes, where the geological media usually consist of fractured rock masses. In this study, we developed a model for the analysis of coupled hydro-mechanical processes in porous rock containing dominant fractures, by using the numerical manifold method (NMM). In the current model, the fractures are regarded as different material domains from surrounding rock, i.e., finite-thickness fracture zones as porous media. Compared with the rock matrix, these fractured porous media are characterized with nonlinear behavior of hydraulic and mechanical properties, involving not only direct (poroelastic) coupling but also indirect (property change) coupling. By combining the potential energy associated with mechanical responses, fluid flow and solid–fluid interactions, a new formulation for direct HM coupling in porous media is established. For indirect coupling associated with fracture opening/closure, we developed a new approach implicitly considering the nonlinear properties by directly assembling the corresponding strain energy. Compared with traditional methods with approximation of the nonlinear constitutive equations, this new formulation achieves a more accurate representation of the nonlinear behavior. We implemented the new model for coupled HM analysis in NMM, which has fixed mathematical grid and accurate integration, and developed a new computer code. We tested the code for direct coupling on two classical poroelastic problems with coarse mesh and compared the results with the analytical solutions, achieving excellent agreement, respectively. Finally, we tested for indirect coupling on models with a single dominant fracture and obtained reasonable results. The current poroelastic NNM model with a continuous finite-thickness fracture zone will be further developed considering thin fractures in a discontinuous approach for a comprehensive model for HM analysis in fractured porous rock masses.