Summary Neural network predictive models are popular for production forecasting in unconventional reservoirs due to their ability to learn complex relationships between well properties and production responses from extensive field data. The intricate flow behavior in hydraulically fractured unconventional reservoirs, which remains poorly understood, makes these statistical models particularly useful. Various neural network variants have been developed for production prediction in these reservoirs, each offering predictive capability of varying levels of granularity, accuracy, and robustness against noisy and incomplete data. Neural network predictive models that integrate physical principles are especially useful for subsurface systems, as they provide predictions that adhere to physical laws. This work introduces a new dynamic physics-guided deep learning (DPGDL) model that incorporates physical functions into neural networks and employs residual learning to compensate for the imperfect description of the physics, under variable data support. The new formulation allows for dynamic residual correction, avoids unintended bias due to less-than-ideal input data, and provides robust long-term predictions. The DPGDL model improves upon a static formulation by utilizing a masked loss function to enable learning from wells with varying production lengths and by improving the results when partially-observed timesteps are present. In addition, a sequence-to-sequence residual model has been developed to correct additional biases in the long-term predictions from the physics-constrained neural networks. Several synthetic data sets with increasing complexity as well as a field data set from the Bakken are used to demonstrate the performance of the new DPGDL model.