本发明的一种基于投影积分的气电耦合园区综合能源系统动态仿真方法，适用于解决具有刚性特征的气电耦合园区综合能源系统仿真问题。其具体算法包括：先进行k次内部积分器的积分运算，再根据第k‑2次、第k‑1次和第k次求解得到的系统状态变量进行一次外部积分器的运算。在内部积分器中，电气系统模型采用前向欧拉法差分化，微分方程和代数方程交替求解，天然气系统模型采用后向欧拉法差分化，偏微分方程和代数方程联立求解；外部积分器采用外推进行求解。本发明方法既能满足仿真所需的数值精度与数值稳定性，又可以在一定程度上提高计算的速度，减小由于计算占用的空间和资源。