Colloidal contaminants infiltrate and can be attached onto grain surfaces of soils and aquifers, where they may persist. In this study, Lagrangian particle tracking is used to investigate particle trajectories and attachment in pore and fracture spaces modeled as three‐dimensional constricted tubes with diverse geometries and orientations relative to gravity. A comprehensive force balance arising from hydrodynamic drag and lift, gravitational settling, Brownian motion, and attractive DLVO interactions is simulated. Results show that the collection efficiency η is primarily governed by the dimensionless settling number 𝑆, representing the relative dominance of gravitational over hydrodynamic forces experienced by the particles. High‐𝑆 scenarios have larger η and are more sensitive to pore orientation, while low‐𝑆 scenarios are more sensitive to pore geometry. For all scenarios but especially low‐S scenarios, the majority of colloid attachment occurs near pore extremities, where fluid velocities are low, such that mechanical remobilization of particles attached is improbable. In low‐𝑆 scenarios, particles may spread and become immobilized at greater distances from the contamination source owing to lower η, are harder to mechanically remobilize as they attach more disproportionately at pore extremities, and have trajectories more sensitive to minor forces, rendering their environmental fates complex. Characterizing the collection efficiency and deposition morphology for various pore space geometries and orientations is crucial in understanding particle fate and developing continuum‐scale models of colloid transport in real soils, where pore spaces are heterogeneous and advection paths are tortuous.