Your example does not work ( on python 2.7.9), so I only sketch my idea:
calculate sp
take the derivative via sp.derivative and evaluate it at the relevant times (probably the same times at which you measured your data)
Set the relevant points to zero (e.g. the value at t=0)
Calculate another spline from the derivative values.
Integrate your spline function. I guess you will have to do this numerically, but that should not be a problem. Do not forget to add a constant, to get your original function.