In this work, we study three higher-dimensional Virasoro integrable models, namely the (3+1)-dimensional Nizhnik-Novikov-Veselov equation, the (3+1)-dimensional breaking soliton equation, and a (3+1)-dimensional extended breaking soliton equation. The three equations are among the Virasoro integrable models. We use the simplified form of the Hirota’s method to establish multiple soliton solutions for each equation. We determine the constraint conditions between the coefficients of the spatial variables to guarantee the existence of the multiple soliton solutions for each model.