This paper focuses on the model reduction problem for a special class of linear parameter-varying systems. This kind of systems can be reformulated as bilinear dynamical systems. Based on the bilinear system theory, we give a definition of the H2 norm in the generalized frequency domain. Then a model reduction method is proposed based on the gradient descent on the Grassmann manifold. The merit of the method is that by utilizing the gradient flow analysis, the algorithm is guaranteed to converge, and further speedup of the convergence rate can be achieved as well. Two numerical examples are tested to demonstrate the proposed method.