There is a method for finding a best distribution of source elements that form a vibratory source array. The method includes inputting plural constraints for the source elements; generating plural distributions of the source elements that fulfill the plural constraints; calculating for each distribution a performance index characterizing the source array; and selecting the best distribution from the plural distributions based on a value of the performance index.