Geologic storage of carbon dioxide (CO2) is one of the potential technological options to mitigate human-induced climate change. Depleted hydrocarbon reservoirs are promising candidates for storing CO2. However, these reservoirs contain residual hydrocarbons that may migrate beyond primary reservoirs during CO2 storage operations. Therefore, it is imperative to quantify the leakage risks of residual hydrocarbons (primarily methane) and CO2 to shallow aquifers. We focus on developing models that quantify leakage risks from wellbores that provide potential leakage pathways. To offset the computational intensity of high-fidelity reservoir simulations, we develop a Reduced-order model (ROM) enabling fast and accurate CO2 and hydrocarbon leakage predictions through wellbores. The ROM is generated using a dataset of 1000 high-fidelity, compositional reservoir simulations of CO2 injection into generic hydrocarbon reservoirs. The input parameters include reservoir depth, permeability, net to the gross ratio (NTG), reservoir pressure multiplier, wellbore permeability, average water saturation, oil compositions, and time. We analyze the performance of various ROM development techniques, including multivariate adaptive regression splines, gradient boosting, and neural networks. While all ROM development techniques yield excellent agreement with R2 higher than 0.95, the neural network model performs best compared to the other two methods. Furthermore, we explore developing a ROM as a collection of multiple sub-ROMs to improve accuracy across a wide range of predictions. We observe that using a set of sub-ROMs outperforms the prediction accuracy of a single ROM. Our work enables fast and reliable risk-assessment tools for CO2 geologic storage in depleted hydrocarbon fields.