Background <p>Parametric reduced-order model (PROM) is constructed by interpolating the substructures in the parameter space. Since the substructures at interpolation points are constructed using the classic Craig-Bampton method, these models may encounter a risk of insufficient accuracy due to the omission of residual flexibility. It leads to a decrease in the accuracy of the PROM, especially for high-order vibration analysis.</p> Methods <p>This paper presents an enhanced high-order PROM modeling method that parameterizes the residual flexibility of substructures at interpolation points based on the direct Taylor expansion, thereby automatically incorporating the effect of the high-order residual flexibility into the PROM.</p> Results <p>Two numerical examples with geometric mistuning are presented to demonstrate the effectiveness of the proposed method, comparing the eigenvalues, modal assurance criterion, frequency response functions, and response of the enhanced PROM with that of the previous PROM and the full-order model.</p> Conclusion <p>Results indicate that this method effectively improves the accuracy of predicting high-order vibration responses.</p>