Based on a surrogate model of the performance function with an adaptive learning strategy, the meta-model-based importance sampling method (Meta-IS) can approximate the optimal importance sampling probability density function (IS-PDF) for estimating failure probabilities, making it an efficient approach for reliability analysis. However, when dealing with extremely small failure probabilities, estimating the normalization factor in the IS-PDF becomes computationally expensive for Meta-IS. To mitigate the computational burden, this paper proposes a meta-model-based double importance sampling method (Meta-IS2) for estimating extremely small failure probabilities. In the proposed method, a hierarchical weighted clustering strategy is designed to construct an IS-PDF for estimating the normalization factor, thereby improving the computational efficiency of the reliability analysis. Under equivalent accuracy, the computational efficiency of the proposed method significantly outperforms that of the existing Meta-IS in the cases with extremely small failure probabilities, as demonstrated through three numerical and one engineering case studies.