The 3D common-reflection-surface (CRS) stack operator depends on eight kinematic wavefield attributes that must be extracted from the prestack data. These attributes are obtained by an efficient optimization strategy based on the maximization of the coherence measure of the seismic reflection events included by the CRS stacking operator. The main application of these kinematic attributes is to simulate zero-offset stacked data; however, they can also be used for regularization of the prestack data, prestack migration, and velocity model determination. The initial implementations of the 3D CRS stack used grid-search techniques to determine the attributes in several steps with the drawback that accumulated errors can deteriorate the final result. In this work, the global optimization very fast simulated annealing algorithm is used to search for the kinematic attributes by applying three optimization strategies for implementing CRS stacking: (1) simultaneous global search of five kinematic attributes of the 3D common-diffraction-surface stacking operator, (2) two-step global optimization strategy to first search for three attributes and then five attributes of the CRS stacking operator, and (3) simultaneous global search of eight kinematic attributes of the CRS operator. The proposed CRS stacking algorithms are applied to land data of the Potiguar Basin, Brazil. It is demonstrated that the one-step optimization strategy of the eight parameters produces the best results, however, with a higher computational cost.