As the retrieval of unconventional oil and gas resources extends to the deep and ultra-deep domains, the issue of cement sheath failure in shale oil wellbores seriously endangers wellbore safety, making it imperative to uncover the relevant damage mechanism and develop effective assessment approaches. In response to the limitations of conventional finite element methods in representing mesoscopic damage, in this study, we determined the mesoscopic parameters of cement paste via laboratory calibrations; constructed a 3D casing–cement sheath–formation composite model using the discrete element method; addressed the restriction of the continuum assumption; and numerically simulated the microcrack initiation, propagation, and interface debonding behaviors of cement paste from a mesomechanical viewpoint. The model’s reliability was validated using a full-scale cement sheath sealing integrity assessment apparatus, while the influences of fracturing location, stage count, and internal casing pressure on cement sheath damage were analyzed systematically. Our findings indicate that the DEM model can precisely capture the dynamic evolution features of microcracks under cyclic loading, and the results agree well with the results of the cement sheath sealing integrity evaluation. During the first internal casing pressure loading phase, the microcracks generated account for 84% of the total microcracks formed during the entire loading process. The primary interface (casing–cement sheath interface) is fully debonded after the second internal pressure loading, demonstrating that the initial stage of cyclic internal casing pressure exerts a decisive impact on cement sheath integrity. The cement sheath in the horizontal well section is subjected to high internal casing pressure and high formation stress, resulting in more frequent microcrack coalescence and a rapid rise in the interface debonding rate, whereas the damage progression in the vertical well section is relatively slow.
Wang et al. (2026) studied this question.