In this study, lattice Boltzmann method (LBM) is utilised for three-dimensional simulation of fluid flow through two porous structures, consisting of grains with the same diameter: (i) a homogeneous porous domain, in which the grains are placed with a simple cubic packing configuration, and (ii) a randomly-packed porous domain. An ultra-fine mesh size is considered to perform the simulations in three orders of magnitude of Reynolds number (), covering laminar to turbulent flow regimes, and capture different flow signatures. Pore velocity fields are derived, and their sample probability density functions (PDF) are analyzed versus time to investigate the dynamics of the flow. The analysis of the PDFs clearly shows that stagnant zones play a significant role in the formation of the pore flow fields, manifested by multimodal PDFs, and the distribution of the velocities in porous media at various cannot be characterized by a single PDF model regardless of the pore structure. While the velocities at the stagnant regions and in the vicinity of the solid boundaries are primarily affected by the viscous forces and exhibit a power-law PDF at different , the velocities in the main (preferential) flow pathways away from the boundaries are shown to be influenced by the inertial forces, hence having an exponential PDF when is low. At high, however, depending on the tortuosity of the porous structure, the velocities may exhibit an exponential or even Laplace PDF.