Satellite image sequences visualise the ocean surface and allow assessing its dynamics. Processing these data is then of major interest to get a better understanding of the observed processes. As demonstrated by state-of-the-art, image assimilation permits to retrieve surface motion, based on assumptions on the dynamics. In this paper, we demonstrate that a simple heuristics, such as the Lagrangian constancy of velocity, can be used and successfully replaces the complex physical properties described by the Navier-Stokes equations for assessing surface circulation from satellite images. A data assimilation method is proposed that adds an acceleration term a(t) to this Lagrangian constancy equation, which summarises all physical processes other than advection. A cost function is designed that quantifies discrepancy between satellite data and model values. This cost function is minimised by the BFGS solver with a dual method of data assimilation. The result is the initial motion field and the acceleration terms a(t) on the whole temporal interval. These values a(t) model the forces, other than advection, that contribute to surface circulation. Our approach was tested on synthetic data and with Sea Surface Temperature images acquired on Black Sea. Results are quantified and compared to those of state-of-the-art methods.