The second approach is based on implementing the discontinuities (DFNs) in a discontinuum numerical model, in addition to use the 3D stress gradient, to test their influence in the pre-mining stress initialization step. Two approaches are adopted to compare the initial stresses, the first approach is based on using the 3D stress gradient in a continuum rockmass without implementing the effect of discontinuities. information Surpac and FLAC3D are chosen to be the modeling and. In fluid-solid coupled model of FLAC3D, the fluid con-. sequence, the in situ stress distribution condition of the mining area has been. The existing discontinuities in the mining area are introduced numerically by using Discrete Fracture Network (DFN) elements and Ubiquitous mechanical model. maximum horizontal principal in situ stress of 9 measuring points were all located between. To reproduce such complex pre-mining stress state in a numerical model (Finite difference numerical model FLAC3D), stress gradients are estimated based on the existing stress measurements. In addition, many faults pass nearby the mining zone. The measured principal stresses in the current case study are highly heterogeneous and anisotropic. The objective of the present paper is to propose a methodology to numerically initialize the pre-mining stress field in complex case study, where the measured vertical stresses are in disagreement with the overburden weight and the exploited coal seam dips with 10° in which numerous in-situ stress measurements were carried out. To assess the stability of underground openings by using the numerical modeling tool, the stress state (i.e., initial and induced stresses) should be generated accurately.