A method of automatic optimization of ROP. The method obtains a plurality of drilling surface parameters for a field of interest, and determines an UCS data and a MSE data for a targeted formation based on well logs. The method further trains a ML model using the drilling surface parameters as inputs, and outputs a plurality of weights for drilling parameters in a ROP equation and in a Teale's MSE equation for the field of interest. The method further combines the ML ROP equation with the Teale's MSE equation to determine a plurality of optimum drilling parameters by simultaneously solving the set of ML ROP equation and the Teale's MSE equation. Furthermore, the method generates a work order to adjust the drilling parameters and cause display of the work order and the determined optimum drilling parameters in a user interface of a client device.