A system and method for controlling multiple drilling tools inside wellbores makes use of Bayesian optimization with range constraints. A computing device samples observed values for controllable drilling parameters such as weight-on-bit, mud flow rate and drill bit rotational speed in RPM and evaluates a selected drilling parameter such a rate-of-penetration and hydraulic mechanical specific energy for the observed values using an objective function. Range constraints including the physical drilling environment and the total power available to all drilling tools within the drilling environment can be continuously learned by the computing device as the range constraints change. A Bayesian optimization, subject to the range constraints and the observed values, can produce an optimized value for the controllable drilling parameters to achieve a predicted value for the drilling parameters. The system can then control the drilling tool using the optimized value to achieve the predicted value for the selected drilling parameter.