Structure-preserving $H^2$ optimal model reduction based on Riemannian trust-region method