Solution 1:

Nobody answered. My teacher gave me some explanation. I trusted him. I let the problem fall. Now I'm back on it to prepare the presentation for my thesis. I thought about it, and here is the explanation I've come up with. With any field $X$ on a manifold $M$, it makes sense to define a two-index family of diffeomorphisms $\phi_s^t$, where $\phi_s^t(x)$ is where you end up at time $t$ if you are in $x$ at time $s$. We can then define:

$$\mathcal L_{X_s}\omega=\frac{d}{dt}(\phi_s^t)^\ast\omega\Big|_{t=s}.$$

This is the Lie derivative along a field that evolves like $X$, only reparametrizing time in such away that what $X$ sees as time $s$ is seen by this field as time 0. So I can call this field $X_s$. Naturally, if $X$ is constant, all these derivatives are identical, and $X_s\equiv X$. What this double index is useful for is when $X$ actually evolves. It is easy to see that, for any $t,s,u$:

$$\phi_s^t\circ\phi_u^s=\phi_u^t.$$

The LHS means you start in $x$ at $u$, and move till time $s$, then you call $y$ the result, and see where you end up at time $u$ if you start at $y$ at time $s$. The RHS does the two steps in one: start in $x$ at $u$ and see where you are in $t$. Of course, when the field is time-independent, this reduces to $\phi_0^t\circ\phi_0^s=\phi_0^{t+s}$, for the LHS's first function is $\phi_s^{t+s}$. That is, if $X$ does not evolve in time, $\phi_s^t=\phi_{u+s}^{u+t}$ for all $u\in\mathbb{R}$. Again, this changes in the time-dependent case.

With that in mind, we move to solving the problem that sparked this question. We have the family $\omega_t=\omega_0+t(\omega-\omega_0)$. We need to find $\phi^t:(\phi^t)^\ast\omega_t\equiv\omega_0$. We try to build a time-dependent $X$ such that:

$$(\phi_0^t)^\ast\omega_t\equiv\omega_0.$$

Deriving this at any time $s$ we have to get 0. But then:

\begin{align*} \frac{d}{dt}\Big|_{t=s}(\phi_0^t)^\ast\omega_t={}&\lim_{h\to0}\frac{(\phi_0^{s+h})^\ast\omega_{s+h}-(\phi_0^s)^\ast\omega_s}{h}={} \\ {}={}&\lim_{h\to0}\frac{(\phi_0^{s+h})^\ast\omega_{s+h}-(\phi_0^s)^\ast\omega_{s+h}}{h}+\lim_{h\to0}\frac{(\phi_0^s)^\ast\omega_{s+h}-(\phi_0^s)^\ast\omega_s}{h}={} \\ {}={}&\lim_{h\to0}\frac{(\phi_s^{s+h}\circ\phi_0^s)^\ast\omega_{s+h}-(\phi_0^s)^\ast\omega_{s+h}}{h}+(\phi_0^s)^\ast\frac{d}{dt}\Big|_{t=s}\omega_t={} \\ {}={}&(\phi_0^s)^\ast\lim_{h\to0}\frac{((\phi_s^{s+h})^\ast-(\phi_s^s)^\ast)\omega_{s+h}}{h}={} \\ {}={}&(\phi_0^s)^\ast\left\{\mathcal{L}_{X_s}\omega_s+\frac{d}{dt}\omega_t\Big|_{t=s}\right\}. \end{align*}

And this, except for a couple more indices, is precisely the desired formula. The proof I have of Cartan's magic formula readily generalizes to this case, so I am finally done on this.

Update

No. Cartan doesn't work as I expected, but rather I managed to show $\mathcal{L}_{X_s}$ as above define doesn't depend on $s$! It always gives $d\iota_{X_0}+\iota_{X_0}d$, where $X_0$ is the field at time 0. I'm really giving up on this. There is NO FUCKING WAY I can set that proof of Darboux straight, because Hofer's way I can justify neither the use of Cartan nor the passage in which $\mathcal{L}_{X_s}$ pops up, nor indeed how I define that derivative. The above way has no Cartan. The field family way (trying to construct a family of fields $X_t$ with flows with two indices, so $X_s$ would have flow $\phi_s^t$) gives a damn $\mathcal{L}_{X_s}\omega_u+\frac{d}{dt}\omega_t\Big|_{t=u}$, so having me conclude all fields are equal. I'm leaving this answer as an attempt to solve the problem, but it didn't solve it AT ALL. It just provided a third non-solution.

Update 2

Turns out all I had to do was define $\mathcal{L}_{X_s}$ in an arbitrary way so as to get the formula, and them prove a "Cartan formula" for $\mathcal{L}_{X_s}$ so defined, to get to the end without too much trouble. And the proof of this "Cartan formula" was identical to that for the real Cartan formula, the one with no time dependency. So I set this proof straight at last.