The decline in the exploration of new oil sites necessitates the development of efficient strategies to maximize recovery from existing reservoirs. This study employs a molecular dynamics (MD) approach to investigate oil detachment from silica surfaces of varying hydrophobicity using a combination of bis-cationic gemini surfactants (GS) and functionalized silica nanoparticles (SNPs). Density profiles and radial distribution function (rdf) plots revealed a multilayered oil adsorption model. A reduction in oil-silica interaction energy was observed with an increase in surface hydrophobicity, highlighting the importance of polar interactions. Standard waterflooding studies, involving oil detachment solely with water, were conducted to assess baseline recovery efficiency. All the GS-SNP combinations outperformed standard waterflooding methods. SNPs significantly mitigated GS adsorption on reservoir beds, as evidenced by center-of-mass measurements. However, the effectiveness of the added injectants (GS-SNP) went downhill with increasing surface hydrophobicity, further validating the existence of a potential barrier for oil detachment, as known previously. Finally, supervised machine learning (ML) models were generated to predict the GS-SNP combination for a given silica surface, with MD generated descriptors. In most cases, boosting models,