Among all properties of reservoir rocks, permeability is an important parameter which is typically measured from core samples in the laboratory. Due to limitations of core drilling all over a reservoir, simulation of rock porous media is demanded to explore more scenarios not seen in the available data. One of the most accurate methods is cross correlation based simulation (CCSIM) which recently has broadly applied in geoscience and porous media. The purpose of this study is producing realizations with the same permeability trend to a real sample. Berea sandstone sample is selected for this aim. Permeability results, extracted from smaller sub-samples of the original sample, showed that classic Kozeny-Carman permeability trend is not suitable for this sample. One reason can be due to lack of including geometrical and fractal properties of pore-space distribution in this equation. Thus, a general trend based on fractal dimensions of pore-space and tortuosity of the Berea sample is applied in this paper. Results show that direct 3D stochastic modeling of porous media preserves porous structure and fractal behavior of rock. On the other hand, using only 2D images for constructing the 3D pore structures does not reproduce the measured experimental permeability. For this aim, a hierarchical sampling is implemented in two and three steps using both 2D and 3D stochastic modeling. Results showed that two-step sampling is not suitable enough, while the utilized three-step sampling occurs to be show excellent performance by which different models of porous media with the same permeability trend as the Berea sandstone sample can be generated.