Reported in #918.
new_scip = Model(sourceModel=orig_scip) loses 226 parameters including dual presolvers, some heuristics. Solving new_scip is usually slow than solving orig_scip with same random seed. The new copy function can full copy orig_scip problem except solution storage.
Also need to find out why it loses some of the original model's parameters.