Seismic petrophysics is an important link between seismic elastic properties and reservoir physical properties. Based on the petrological and microstructure characteristics of shale in the Qingshankou formation of Sanzhao sag in the north of Songliao Basin, this paper presents an anisotropy petrophysical model with complex pore structure suitable for organic shale constructed with the use of the Voigt-Reuss-Hill average model, an anisotropy self-consistent approximation+differential effective medium model, and the layering of clay and kerogen is simulated by using the Voigt-Reuss-Hill average and bond transform to achieve the simulation of shale anisotropy. Based on the proposed model, the effects of the organic volume fraction, porosity, and pore aspect ratio on rock elastic properties are discussed. The result shows that with the increase of matrix porosity, all elastic parameters show a decreasing trend; with the increase of the organic volume fraction, except shear modulus, other elastic parameters show an increasing trend. Through comparative analysis, the elastic parameters (Lamé impedance and Shear impedance) sensitive to the organic volume fraction and porosity are optimized; the seismic petrophysical cross-plot template with core calibration is constructed. The application shows that the predicted S-wave velocity based on the proposed model is in good agreement with the S-wave velocity derived from dipole source logging. Combined with the high-precision prestack elastic parameter inversion, the “sweet spot” characteristics can be well described, and the research could contribute to a better “sweet spot” description and provide a better support for shale exploration in Sanzhao sag.