Methods and devices according to various embodiments perform full-wave inversion jointly for datasets acquired at different times over the same underground formation using a time-lag cost function with target regularization terms. This approach improves the 4D signal within reservoirs and suppresses 4D noise outside.