A 3D basin model based on a detailed stratigraphic database was built on the Paris basin. In this paper, the model ability to reproduce the permeability heterogeneity at the regional scale is tested. In a ''genetic'' way, the model reproduces processes such as sedimentation, erosion and compaction. In a geostatistical way, a facies model, as a major input for the basin simulator is generated by means of kriging and conditional simulations. A precise facies model is attributed to the Callovo-Oxfordian semi-permeable and to the Dogger and Lusitanian aquifers; a ''discrete'' facies model is attributed to the other layers. Petrophysical laws relating intrinsic permeability to porosity are tested to describe the permeability within the three layers of interest: Archie's law, fractal model and Kozeny-Carman relationship for the siliciclastic sediments and petrographic relationships for the carbonates. The available permeability data for the two aquifers are upscaled and used to calibrate the models, providing a description of the heterogeneity associated with facies and porosity in agreement with those shown by the data. Conditional simulations of facies distribution allows us to estimate the dispersion on the results.