Batch-sequential design and heteroskedastic surrogate modeling for delta smelt conservation