Numerical approximation of poroelasticity with random coefficients using Polynomial Chaos and Hybrid High-Order methods · Study