Parametrization including frequency, dip and spacial coordinates allows to map particular smoothly continuous events such as wavefronts or reflectors into separate or in practice almost separate sets of coefficients, which makes curvelet domain favorable in regard to seismic data processing. Even if the subsurface structure is complicated with conflicting dips, intersecting curved events, faults, etc. the curvelet transform is still able to decompose the data into linear sum of needle-shaped atom functions. Our work-flow based on curvelet denosing already proved its robustness with pre-stack and 2D/3D post-stack seismic data.