Simultaneously modelling steam flow and heat transfer in vertical and horizontal wellbores during SAGD cyclic pre-heating stage