Abstract: A seismic data analysis method including: loading at least one seismic data volume to a database, performing unbiased sampling on a portion of the seismic data volume to obtain a plurality of unbiased sampling traces, selecting a plurality of trace positions from the unbiased sampling traces, computing a plurality of samples, identifying a first set of wavelets using the samples, computing at least one genotype, creating at least two subpopulations by applying the genotype to the seismic data volume, analyzing the subpopulations for fitness to identify at least one fit population, forming a merged subpopulation by mapping similarities in the first set of wavelets, computing a second sample by identifying a second set of wavelets, creating a mutated subpopulation by applying the genotype to the second sample, and inserting the mutated subpopulation into the seismic data volume to create a geopopulation that defines geological parameters of a subsurface formation.