Correct modeling of the geometric architecture of the sand/shale geometry within shaly layers of the Wilmington field is critical to predict performance recovery accurately. Stochastic images of that clastic sequence are obtained with sequential indicator simulation with non-Cartesian coordinates to ensure stratigraphic conformity.