Polynomial fitting can increase signal-to-noise ratio without damaging resolution. In this paper, a 3D curved-surface fitting and stacking method is introduced. It can separate and compound linearity fitting, hyperbola fitting and flexure fitting using decomposability of 3D compound orthogonal polynomial. A 3D seismic dataset is processed in offset domain by this method. High quality section is obtained. Processing results prove that the method can improve signal-to-noise ratio and resolution effectively.