本发明公开了一种基于蒙特卡罗方法与特征线方法耦合的计算辐射屏蔽的方法，该方法通过耦合这两种方法，在几何复杂的辐射源区域采用蒙特卡罗方法模拟，而针对几何简单的屏蔽层区域则采用非均匀修正的特征线方法模拟来获取屏蔽层中的辐射粒子通量分布结果，其中，蒙特卡罗方法优点在于能够精确地模拟多维复杂几何下的粒子输运过程，而特征线方法优点在于其对计算模型几何限制小，计算速度快且不存在射线效应，可以获得高精度的计算结果。相对传统的屏蔽计算方法而言，本发明可以对深穿透屏蔽区域获得更好的几何适应性、更高的计算效率和更精确的计算精度。