A method for combustion tuning, comprises collecting exhaust parameters indicating combustion status of a boiler by a sensor array; determining whether the exhaust parameters of the boiler match a preset optimization target; and optimizing combustion, if the exhaust parameters do not match the preset optimization target by selecting a model from a model repository based on a current boiler condition, wherein the model corresponds to a relationship between model input variables and the exhaust parameters; determining at least one optimized model input variable of the boiler for realizing the optimization target, based on the selected model; and adjusting actuators of the boiler according to the optimized model input variable.