本发明提供了一种计及可靠性的多能耦合系统分布式协同优化方法。本发明方法步骤如下：步骤S1，可靠性指标和其他指标组成共a个指标集，在由a个指标构成的a维指标空间中进行节点的聚类；步骤S2，生成子系统整体可靠性指标和其他指标；步骤S3，利用抽样方法获取扰动场景；步骤S4，建立分布式协同优化模型，利用优化算法对多子系统不同优化目标的优化调度模型进行求解。本发明能够为多能耦合系统提供有效的子系统可靠性分析方法，实现计及生产生活差异的多层级、多主体协同优化，为其在生产生活过程扰动下满足不同主体的不同优化需求并实现系统效率和可靠性协同提供指导。