AbstractGiven the growing importance of organic-rich shale as unconventional reservoirs, a thorough understanding of the elastic and anisotropic behavior of shales is of great concern. However, for lacustrine shales, the complex lithofacies assemblage with geological deposition makes it challenging. Four lithofacies (argillaceous, mixed, siliceous, and calcareous) are recognized for 40 lacustrine shale samples from Jurassic formation in Sichuan Basin on the basis of their mineral compositions. We perform ultrasonic velocity measurements on 40 pairs of shale plugs at varied confining pressures, attempting to uncover the controls on the anisotropic properties of different lithofacies. The experimental results reveal that the total porosity, clay, and organic matter would positively contribute to velocity anisotropy of Jurassic shales. Combined with microstructure and pressure-dependent velocity analysis, the preferred orientations of platy clay particles and lenticular kerogen, the development of clay pores along clay fabric, and the sub-parallel micro-cracks induced by hydrocarbon expulsion are treated to be the controlling mechanisms. We add the total porosity, clay content, and kerogen volume together, intending to distinguish the elastic and anisotropic properties of four lithofacies. Generally, argillaceous shales, the dominant lithofacies in Jurassic formation, could be characterized by the highest clay and total organic content (TOC), the lowest bedding-normal velocities, and the strongest velocity anisotropy. Finally, with the laboratory data, two rock-physics-driven exponential relationships are proposed to predict the P- and S-wave velocity anisotropy with the bedding-normal velocities.