A method and system includes determining a subsurface fluid seismic attribute comprising recording a signal at a wellhead related to pumping fracture stimulation fluid to obtain a pressure pulse pump signal. A deconvolution operator is determined from the obtained pressure pulse pump signal. Seismic data are acquired from a plurality of sensors. Travel time differences are computed for the seismic data between the plurality of sensor locations and the subsurface position. Seismic data are deconvolved with the deconvolution operator to obtain a plurality of deconvolution coefficients associated with the subsurface position. The computed travel time differences are used to sum the plurality of deconvolution coefficients associated with the subsurface position to obtain a subsurface fluid seismic attribute at the subsurface position.