This paper addresses the issue of Naturally Fractured Reservoirs (NFR), where fracturation patterns are assumed to be closely related to stress history. Introducing uncertainties on geologic data, we propose to compute a new stochastic parameter the Failure Probability based on geomechanics. This parameter can then be used to simulate 3D fracture density and orientation analogue.