Abstract: In the unresolved energy region, due to the dense distribution of resonance peaks, the precise resonance parameter data cannot be obtained due to the resolution limitation of the nuclear measuring instrument. For this energy region, the Evaluation Nuclear Data Files only provides the average values and distribution functions of the resonance parameters. In order to consider the resonance self-shielding effect in the unresolved resonance energy region, it is necessary to calculate the effective self-shielding cross-section based on these parameters with probability distribution properties. For the treatment of resonance self-shielding effect in the unresolved resonance energy region, a probability table module PUnresXS and an effective self-shielding cross-section calculation module UnresXS were developed in the advanced nuclear cross-section processing program AXSP based on the "Ladder Sampling" and integral statistics principle method. The algorithms in the probability table calculation module were also improved. By comparing with the reference solution obtained from NJOY2016 calculations, the PUnresXS and UnresXS modules were independently validated. The results showed that the efficiency of the probability table module was improved by more than 60%, the k_eff and neutron spectrum were in good agreement with the reference solution. The developed resonance self-shielding cross-section calculation module exhibited similar accuracy to NJOY2016, demonstrating the correctness of these two modules developed in the AXSP program.