Abstract: The shape of rock and soil particles is distinct and their orientation may prefer a certain direction, and determining the influence of particle shape and orientation on the characteristics of granular column collapse is of great significance to geological disasters prevention. The collapse process of granular columns is simulated with the superquadric discrete element method, and the influence of particle shape (aspect ratio A and blockiness B) and initial orientation angle θ of long axis on the collapse characteristics of granular column is analyzed, and further the underlying physical mechanism is explored. The results show that, compared with the case A=1 (isotropic), the collapse of particles with larger A (elongated) or smaller A (platy) results in shorter final runout distance and increases the deposit height. The collapse of particles with larger B (more angular) results in shorter runout distance and higher deposit height. The column with θ ≈ 0° (near horizontal) particles produces an obviously smaller final runout distance, and moves in a way that some particles being first squeezed out from the lateral free surface followed by the gradual slide of particles in the upper parts. For short column, the flow regime is slip-dominant, and the collapse of particles with θ ≈ 135° (near the orientation of subsequent accumulation slope) produces the longest final runout distance. For tall column, the early falling effect is obvious, and the collapse of particles with θ ≈ 90° (nearly upright) has the longest final runout distance. The influence of different particle shapes and orientation angles on the characteristics of granular column collapse is closely related to the transformation process of energy and the meso-structure.