For a symbolic distribution dist, the r factorial moment is given by Expectation[FactorialPower[x, r], xdist].
For a multivariate symbolic distribution dist, the factorial moment is given by Expectation[FactorialPower[x1, r1]FactorialPower[x2, r2]..., {x1, x2, ...}dist].