In efforts to further improve final migrated images we have developed a new methodology for post-migration multiple removal in the migration depth domain. The typical input for the prediction phase of this process is a 3D depth migration volume and the corresponding velocity field. A post-stack Wavefield Extrapolation (WFE) based multiple prediction is used to identify/confirm the multiple events in the migration depth domain. Once multiple events are identified, an effective and efficient demultiple procedure called interactive demultiple is applied to remove the residual multiple from the final migration. The key ingredient of this new interactive demultiple methodology is the attribute-based subtraction. We will describe the main steps of this methodology, and demonstrate its effectiveness by showing some field data applications.