Nuclear magnetic resonance (NMR) relaxation experimentation is an effective technique for probing the dynamics of proton spins in porous media but interpretation requires the application of appropriate spin diffusion models. Molecular dynamics (MD) simulations of porous silicate-based systems containing a quasi-two-dimensional water-filled pore are presented. The MD simulations suggest that the residency time of the water on the pore surface is in the range 0.03-12 ns, typically 2-5 orders of magnitude less than values determined from fits to experimental NMR measurements using the established surface-layer (SL) diffusion models of Korb and co-workers [Phys. Rev. E 56, 1934Rev. E 56, , (1997]. Instead, MD identifies four distinct water layers in a tobermorite-based pore containing surface Ca 2+ ions. Three highly-structured water layers exist within 1 nm of the surface and the central region of the pore contains a homogeneous region of bulk-like water. These regions are referred to as layer 1 and 2 (L1, L2), transition layer (TL) and bulk (B), respectively. Guided by the MD simulations, a two-layer (2L) spin-diffusion NMR relaxation model is proposed comprising two two-dimensional layers of slow-and fast-moving water associated with L2 and layers TL+B respectively. The 2L model provides an improved fit to NMR relaxation times obtained from cementitious material compared to the SL model, yields diffusion correlation times in the range 18-75 ns and 28-40 ps in good agreement with MD, and resolves the surface residency time discrepancy. The 2L model, coupled with NMR relaxation experimentation, provides a simple yet powerful method of characterising the dynamical properties of proton-bearing porous silicate-based systems such as porous glasses, cementitious materials and oil-bearing rocks.