Optimisation of hydrokinetic turbine array layouts via surrogate modelling