The method presented here is a modeling technique for surface zero-offset multiples in 3D. It does not require any picking, and the multiples are modeled from the data themselves by using post-stack migration tools and a demigration with a delayed exploding reflector condition (DDERC). This leads to a kinematically exact technique for modeling the zero-offset double travel-path surface multiples produced by any generator with any dip, such as the sea-floor first-order multiple, for example. The 3D multiple model thus found can then be subtracted from the zero-offset input sections by using waveform and amplitude adaptive techniques. The cost-effectiveness of this modeling method resides in the different levels of approximations that can be accepted, either in the migration or demigration stage, or the in the kind of multiples that we intend to model. This 3D multiple modeling method has been implemented in the f-k domain by using the phase-shift technique. It is thus in principle absolutely exact for 1D velocity backgrounds only. Nevertheless, the robustness of the method has been verified in laterally variant velocity media in deep seafloor contexts. An approximation of the kinematics of the demigration with the delayed exploding reflector condition (DDERC), based on Taylor's expansions of the traveltimes curves corresponding to the DDERC operator, leads us to a cheaper approach for the constant velocity case based on Stolt's migration and demigration techniques by modifying the velocity value. At last, we can show that some additional approximations allow us not only to model accurately the double travel path multiples, but also the peg-legs and the diffracted multiples.