Nonclassical correlations provide a resource for many applications in quantum technology as well as providing strong evidence that a system is indeed operating in the quantum regime. Optomechanical systems can be arranged to generate quantum entanglement between the mechanics and a mode of travelling light. Here we propose automated optimisation of the production of quantum correlations in such a system, beyond what can be achieved through analytical methods, by applying Bayesian optimisation to the control parameters. Two-mode optomechanical squeezing experiment is simulated using a detailed theoretical model of the system, while the Bayesian optimisation process modifies the controllable parameters in order to maximise the non-classical two-mode squeezing and its detection, independently of the inner workings of the model. The Bayesian optimisation treats the simulations or the experiments as a black box. This we refer to as emph{theory-blind} optimisation, and the optimisation process is designed to be unaware of whether it is working with a simulation or the actual experimental setup. We find that in the experimentally relevant thermal regimes, the ability to vary and optimise a broad array of control parameters provides access to large values of two-mode squeezing that would otherwise be difficult or intractable to discover. In particular we observe that modulation of the driving frequency around the resonant sideband, when added to the set of control parameters, produces strong nonclassical correlations greater on average than the maximum achieved by optimising over the remaining parameters. We also find that using our optimisation approach raises the upper limit to the thermal regime in which squeezing can be achieved. This extends the range of experimental setups in which non-classical correlations could be generated beyond the region of high quantum cooperativity.