This study systematically and experimentally evaluates data integration methods for the isotopic analysis of Pb at ultratrace levels using thermal ionization mass spectrometry (TIMS) with a continuous heating method. The evaluation utilized a certified reference material of Pb (SRM 981). The experimental evaluations encompass different data calculation methods (methods I, II, and III) and integration ranges (full, over 1%, 25%, and 75%). Method I, in which isotope ratios were calculated based on summed ion signal intensities compensating for mass fractionation, was consistent with the certified values for 10 and 1 ng standard samples across all integration ranges. For 100 pg samples, full range calculations failed for specific isotope ratios, but reduced ranges (over 1%, 25%, and 75%) yielded values overlapping with certified ones. Method II, in which isotope ratios were calculated by averaging the precalculated isotope ratios, exhibited inferior performance compared to method I. Method III, using weighted averaging to reduce anomalous values, showed results consistent with those of method I but was recommended only for single measurements. An integration range of over 1% or 25% is preferred to exclude anomalies while compensating for mass fractionation. The optimized method was validated by comparing two different instruments used for the isotopic analysis of the reference material. The enhanced accuracy and precision provide valuable insights for researchers working in ultratrace-level Pb isotopic analysis using TIMS.