This invention provides a method for determining the integrity support data of BeiDou PPP-B2b service. It utilizes long-term statistical data from BeiDou PPP-B2b messages, broadcast
ephemeris, and precise
ephemeris to calculate instantaneous
ranging error and worst-case spatial
signal ranging error. A two-step over-boundary
envelope method is used to enclose the worst-case spatial
signal ranging error distribution to obtain a two-step
Gaussian over-boundary envelope distribution, subsequently determining the sum of the values in the ISD parameters. Based on the sum, the upper
limit value UB of the worst-case spatial ranging error is obtained. The instantaneous spatial
signal ranging error within the PPP-B2b service area is calculated using a grid
simulation method. Then, using the upper
limit value UB as an
anomaly detection threshold to remove outliers, the standard deviation is calculated to obtain the sum of the values in the ISD parameters. Finally, using the upper
limit value UB as an
anomaly detection threshold, the number of faults and the
average duration of faults for individual satellites and constellations are statistically analyzed to calculate the sum of the values in the ISD parameters. This provides accurate prior input parameters for ARAIM, improving the reliability of ARAIM fault monitoring,
elimination, and protection level calculation. This method can also be used as an evaluation method for the long-term service performance of PPP-B2b.