Voltage sags and power interruptions are important power quality problems that affect sensitive customers, mainly because they cause annual massive economical losses to the industrial sector as a result of unexpected production process disruptions. In this sense, to propose corrective and preventive measures and improve the power quality of distribution systems, stochastic methodologies have been proposed in the literature to estimate annual voltage sags and power interruptions. However, these methodologies generally use typical cumulative distribution functions of voltage sag duration (PSgD), which can not reflect the real estate of the network under study. To solve this constraint, this paper proposes a novel methodology to estimate a proper PSgD considering the information of the distribution network (i.e., topology and coordination schemes of the protection system) and the stochastic behaviors of short-circuits that can affect the distribution system. Moreover, the proposed methodology allows estimating permanent failure rates and average repair time considering known or expected values of reliability indicators. The results show that the proposed methodology is capable to adapt from an initial PSgD curve to another one with fidelity to achieve real values of expected annual power interruptions.