To evaluate pore structures of the Bakken Shale, which is one of the most important factors that affect petrophysical properties, high-pressure mercury intrusion was employed in this study. Pore structures such as pore-throat size, pore-throat ratio, and fractal attributes are investigated in this major shale play. Pore-throat size from 3.6 to 200um is widely distributed in these shale samples. Accordingly, pore-throat size distributions demonstrate the multimodal behavior within the samples. The whole pore-throat network can be divided into four clusters: one set of large pores, two transitional/intermediate pore groups, and one set of smaller pores. The fractal analysis revealed that fractal dimensions decrease as the pore-throat size decreases. The multifractal analysis demonstrated that as the maturity of the shale samples increases, pore-throat size distributions would become more uniform and pore structures tend to become more homogeneous. The results are compared to our previous results obtained from nitrogen gas adsorption for further verifications of fractal behavior. Finally, although fractal analysis of mercury intrusion and nitrogen gas adsorption were comparable, the results of multifractal analysis from these two methods were not identical.