Several problems of interest to Process Systems Engineering can be formulated as nonconvex quadratically constrained programs (QCPs). These are challenging to solve to global optimality due to the existence of many local optimal solutions and the difficulty of computing a tight dual bound. Although recent research has shown that mixed-integer programming (MIP) relaxations of the bilinear terms in the constraints can play an important role addressing this challenge, because they take considerably more time to solve than their linear programming counterparts, state-of-the art commercial solvers have found limited use for them. We now propose a global optimization algorithm relying extensively on MIP relaxations, for reducing the variable domain through optimality-based bound tightening and improving the dual bound. MIP relaxations are derived from the multiparametric disaggregation technique, using a logarithmic partitioning scheme to ensure tractability and base-2 to minimize the jump in complexity when increasing the number of intervals in the partition for a subset of the bilinear variables. Major improvements are achieved doing this sequentially, with refinements occurring within a spatial branch and bound procedure. Using a set of 54 benchmark instances from the literature, we show that the new algorithm significantly outperforms GUROBI 9.5.2 and BARON 22.7.23, in terms of optimality gap at termination (reduction up to four orders of magnitude), and number of problems solved.