Nanoparticle (NP) transport is increasingly relevant to subsurface engineering applications such as aquifer characterization, fracture electromagnetic imaging and environmental remediation. An efficient field-scale simulation framework is critical for predicting NP performance and designing subsurface applications. In this work, for the first time, a streamline-based model is presented to simulate NP transport in field-scale subsurface systems. It considers a series of behaviors exhibited by engineered nanoparticles (NPs), including time-triggered encapsulation, retention, formation damage effects and variable nanofluid viscosity. The key methods employed by the algorithm are streamline-based simulation (SLS) and an operator-splitting (OS) technique for modeling NP transport. SLS has proven to be efficient for solving transport in large and heterogeneous systems, where the pressure and velocity fields are firstly solved on underlying grids using finite-difference (FD) methods. After tracing streamlines, one-dimensional (1D) NP transport is solved independently along each streamline. The adoption of OS enhances flexibility for the entire solution procedure by allowing different numerical schemes to solve different governing equations efficiently and accurately. For the NP transport model, an explicit FD scheme is used to solve the advection term, an implicit FD scheme is used for the diffusion term and an adaptive numerical integration is used to solve the retention terms.