Table summarizes the simulation setup for each oil case. For all cases, we draw 2500 samples from the posterior distributions and used 25 Monte Carlo steps in each CATMIP stage, which are draws and n steps CATMIP arguments in PyMC3, respectively. We used two heavier components characterization for heavy oils: based on Soave-Redlich-Kwong 17 (SRK) and Peng-Robinson 35 (PR) cubic equations of state. The characterized and lumped mixture for each oil case is provided at the Supporting Information file (see A.1). The posterior distributions for CSP tuning parameters and standard deviation (hyper-parameter) are available at A.2. For the Elementary Effects method, we used p = 4 grid levels and 40 trajectories. We considered a region of experimentation construct by a range of ±25% for each CSP tuning parameter, i.e., 0.75 to 1.25. The results for each case is provided below.