本发明公开了一种考虑多重不确定性的天然气概率‑模糊能流分析方法，主要包括以下步骤：1)建立天然气系统模型。2)建立所述非平衡节点负荷的概率模型。3)利用拉丁超立方法对所述非平衡节点负荷的概率模型抽样，从而得到N×NL维负荷概率样本矩阵。4)建立管道综合参数的模糊模型。5)利用α‑截集法对所述管道综合参数的模糊模型进行抽样。6)利用牛顿法计算所述天然气系统中的管道能流和管道非平衡节点承受压力。7)计算得到输出变量的可能性测度Pos和必然性测度Nec，从而得到输出变量的信任度函数Bel和似真度函数Pl。本发明结合概率和模糊两类不确定因素的影响，可以有效、准确的计算出天然气系统输出变量的变化范围。