Modelling binary black hole waveforms using Gaussian process regression