Fractal geostatistics can be used to generate distributions of such reservoir properties as porosity and permeability. We compare three different methods that use these distributions in reservoir simulations. The first method is a hybrid finite-difference/stream-tube method that starts with 2D vertical cross sections with relatively fine grids, the second uses pseudo-relative permeabilities and also is based on 2D cross sections, and the third uses 3D distributions. A mature waterflood in a carbonate reservoir serves as a field test case.