Accurate horizon recognition within poststack seismic sections is important in seismic interpretation. Horizon picking techniques influence diverse seismic structural analyses and inversion methodologies. Despite encountering challenges such as computational demands and time constraints, the past few years have witnessed the development of numerous 2D and 3D methods. With the progress of modern computing, more robust and efficient techniques, harnessing artificial intelligence, have come to the fore. In this context, we develop an innovative cost-effective algorithm grounded in a global optimization framework. This algorithm combines the very fast simulated annealing method with a set of stability parameters and coherence measurements between neighboring seismic traces, which is used as the objective function. The search is executed on continuous and sequential clusters of seismic traces, with a focus on maximizing coherence amidst the presented events. We evaluate the algorithm’s efficacy by applying it to two distinct input data sets — enveloped and nonenveloped seismic traces. The initial application is to the time-migrated Marmousi data set, a synthetic marine data set originated from a complex geologic sedimentary basin situated in Angola, Africa. Subsequently, we apply the algorithm to depth-migrated land data extracted from a past survey conducted in the Tacutu Basin in northern Brazil. Outcomes arising from the application to 2D synthetic and field data sets underscore the method’s viability as a compelling alternative for seismic horizon picking within the time or depth domain.