The invention discloses a fatigue crack initiation prediction method considering surface integrity, and belongs to the field of engineering material fatigue failure. According to the method, a fatigue crack initiation micro-nano model is provided based on the crystal plasticity constitutive structure and in combination with dislocation strain energy, the fatigue crack initiation position, direction and cycle can be accurately predicted by establishing the fatigue crack initiation prediction model and in combination with crystal plasticity finite element simulation, the influence of the integrity of the machined surface, including the surface roughness, residual stress and surface layer texture, can be considered at the same time, and the method plays an important role in anti-fatigue manufacturing and failure analysis of important parts; a multilevel microstructure modeling method is adopted, martensite microstructure modeling is achieved, meanwhile, the relation between martensite block and martensite beam crystal orientation is combined, the martensite crystal orientation is accurately established, and compared with a Thiessen polygon modeling method, the method has the advantage of being high in precision, and the method plays an important role in mechanical property analysis and failure analysis of martensite materials.