1. A method of determining a horizon-picking solution within a geological volume of interest, the method comprising: obtaining volume information representing a geological volume of interest; obtaining an initial horizon-picking solution for the geological volume of interest derived from the volume information, wherein the initial horizon-picking solution includes an initial set of horizons through the geological volume of interest; and determining a final horizon-picking solution for the geological volume of interest to obtain a final set of horizons through the geological volume of interest by perturbing the horizons in the initial set of horizons simultaneously to follow local features of the geological volume of interest represented in the volume information using the following equation: E =  [ [ ɛ ⁢ ⁢ W s ⁡ ( τ ) ⁢ S W d ⁡ ( τ ) W d ⁡ ( τ ) ] ⁡ [ I ▽ x ▽ y ] ⁢ τ - [ τ 0 p x ⁡ ( τ ) p y ⁡ ( τ ) ] ]  2 where S is a smoothing operator, ε is a smoothing parameter, W d is a diagonal weight applied to a dip fitting goal, W s is a diagonal weight applied to a stratal slicing goal, I is an identity matrix, τ is the final horizon-picking solution, and τ 0 is the initial horizon-picking solution.