Gasoline is one of the most valuable products in an oil refinery and can account for as much as 60-70% of total profit. Optimal integrated scheduling of gasoline blending and order delivery operations can significantly increase profit by avoiding ship demurrage, improving customer satisfaction, minimizing quality give-aways, reducing costly transitions and slop generation, exploiting low-quality cuts, and reducing inventory costs. In this article, we first introduce a new unit-specific event-based continuous-time formulation for the integrated treatment of recipes, blending, and scheduling of gasoline blending and order delivery operations. Many operational features are included such as nonidentical parallel blenders, constant blending rate, minimum blend length and amount, blender transition times, multipurpose product tanks, changeovers, and piecewise constant profiles for blend component qualities and feed rates. To address the nonconvexities arising from forcing constant blending rates during a run, we propose a hybrid global optimization approach incorporating a schedule adjustment procedure, iteratively via a mixed-integer programming and nonlinear programming scheme, and a rigorous deterministic global optimization approach. The computational results demonstrate that our proposed formulation does improve the mixed-integer linear programming relaxation of Li and Karimi, Ind. Eng. Chem. Res., 2011, 50, 9156-9174. All examples are solved to be 1%-global optimality with modest computational effort.Before we develop our mathematical formulation for this scheduling problem, we need to select one time representation. Although several time representations have been proposed in the literature, 23-26 the advantages of an approach based on unit-specific event-based model [27][28][29][30][31][32][33][34][35] are well established in the literature. It can significantly reduce bilinear terms, which determine the complexity of the nonconvex model. 36 The unitspecific event-based time representation is illustrated in Figure 2. The details about this time representation can be referred to Figure 2. Event points definition for each unit.We drop all nonconvex bilinear constraints (i.e., Eqs. 22-25) from the original MINLP models (i.e.LK2011: The model of Li and Karimi (2011). RMILP is the relaxed MILP value from the MPM-NRD or MPM of LK2011. a The entire LMPM-NRD model. b LMPM-NRD model without tightening constraints 118-120.