1.一种基于绿证协商交易的综合能源服务商合作运行方法，其特征在于，包括以下步骤：（1）采集市场中可再生能源发电成本和绿证单价数据，建立绿证协商交易模型；所述的绿证协商交易模型包括绿证交易约束模型和绿证成本模型；

（2）基于绿证协商交易模型，建立用于计算由绿证获得免费碳配额的绿证‑碳交互模型；所述绿证交易约束模型如下：；

其中，gmax为绿证最高价格； gs为单位可再生能源发电政府补贴金额；gmin为最低绿证价格；S、Q为各服务商获得和需求的绿证数量；g0为绿证基准价格； 为供需调整系数； 为成本调整系数；y为服务商系统内部的可再生能源发电成本， 为可再生能源的平均发电成本，g为综合能源服务商的交易价格；

所述绿证成本模型如下：

；

其中，Cgreen为卖方服务商的绿证成本；Psum(t)为t时刻服务商系统内化石能源设备的总出力，所述化石能源设备包括燃气轮机、燃气锅炉和热电联供装置；Pnew(t)为t时刻服务商系统内可再生能源的总出力；为各服务商的可再生能源责任配额系数； 为可再生能源发电转化为绿证的量化系数；

所述的绿证‑碳交互模型包括由绿证获得免费的碳配额计算模型和碳排放成本计算模型；

所述碳配额计算模型如下：

；

其中，Bst碳排放额；Bstg为综合能源系统需承担的的碳排放配额， 为单位功率产生的碳配额与分配碳配额的差值；Psum(t)为t时刻服务商系统内化石能源设备的总出力，所述化石能源设备包括燃气轮机、燃气锅炉和热电联供装置；Bgreen为可再生能源发电转化成绿证所获得的碳配额， 为转化系数；

所述碳排放成本计算模型如下：

；

其中，CET为碳排放成本，为单位碳配额成本；n为阶梯长度；为增长率；

（3）基于设备运行成本、绿证成本和碳排放成本，建立综合能源服务商日运行成本纳什谈判模型；

所述的综合能源服务商日运行成本纳什谈判模型如下：

；

其中，Ux为服务商x的日运行成本； 为谈判破裂点，即服务商没有参与合作运行的情况下各自的综合成本，N为参与合作运行的服务商个数；Cop为设备运行成本；Cgreen为卖方服务商的绿证成本；CET为碳排放成本；Pnew(t)为t时刻服务商系统内可再生能源的总出力；

为各服务商之间的功率交互成本；Pi(t)和Hi(t)分别表示在t时刻的化石能源供电设备和化石能源供热设备的输出功率， 、 分别为电热化石能源设备的单位功率成本；

、 分别为电热交互功率的单位成本；Pimax和Himax分别表示化石能源供电设备和化石能源供热设备的最大输出功率，n和m分别表示化石能源供电设备和化石能源供热设备的数量， 和 分别表示服务商在t时刻的电热交互传输功率，HL(t)、PL(t)为t时刻系统的热电负荷大小；

（4）基于ADMM算法原理对纳什谈判合作运行模型进行转换，等效为联盟成本最小化和收益分配两个子问题；

（5）采用灰狼搜索算法对联盟成本最小化和收益分配两个子问题进行求解，确定各服务商之间的交互功率、绿证交成本，并根据求解结果对服务商内的化石能源和可再生能源设备的出力进行合理调度。

2.根据权利要求1所述的一种基于绿证协商交易的综合能源服务商合作运行方法，其特征在于，所述步骤（4）具体为：（4.1）基于ADMM算法原理对纳什谈判合作运行模型进行转换，联盟成本最小化子问题具体为：根据均值不等式知：

；

Ux为服务商x的日运行成本；为谈判破裂点，即服务商没有参与合作运行的情况下各自的日运行成本，N为参与合作运行的服务商个数；当服务商x的成本降低量 相等时，纳什议价合作运行模型取得最值，在此条件下，能够将变量相乘的非凸优化问题转化成变量相加的凸优化问题；

各服务商进行合作运行时，绿证成本和收益相互抵消，将模型简化为联盟成本最小化子问题，该问题主要求解的是各服务商之间的共享交互功率大小和自身系统内的化石能源和可在生能源设备的出力，具体子问题如下：；

其中，Cop,x、CET,x、Cinter,x分别为服务商的设备运行成本、碳排放成本和能源交互成本；

（4.2）基于ADMM算法原理对纳什谈判合作运行模型进行转换：求解联盟成本最小化子问题后，确定服务商之间的交互功率、化石能源和可再生能源设备运行功率，在此基础上，将绿证协商中的绿证成本简化为收益分配子问题，将联盟成本最小化子问题的最优解 带入纳什谈判合作运行模型中，该问题主要求解的是各服务商之间的绿证成本，具体收益分配子问题如下：；

其中，Cgreen,x为服务商x的绿证成本/收益。

3.一种基于绿证协商交易的综合能源服务商合作运行系统，其特征在于，包括以下模块：绿证协商交易模型建立模块：用于采集市场中可再生能源发电成本和绿证单价数据，建立绿证协商交易模型；所述的绿证协商交易模型包括绿证交易约束模型和绿证成本模型；

绿证‑碳交互模型建立模块：基于绿证协商交易模型，建立用于计算由绿证获得免费碳配额的绿证‑碳交互模型；

综合能源服务商日运行成本纳什谈判模型建立模块：基于设备运行成本、绿证成本和碳排放成本，建立综合能源服务商日运行成本纳什谈判模型；

模型转换模块：基于ADMM算法原理对纳什谈判合作运行模型进行转换，等效为联盟成本最小化和收益分配两个子问题；

求解调度模块：采用灰狼搜索算法对联盟成本最小化和收益分配两个子问题进行求解，确定各服务商之间的交互功率、绿证交成本，并根据求解结果对服务商内的化石能源和可再生能源设备的出力进行合理调度。

4.一种计算机存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现如权利要求1‑2中任一项所述的一种基于绿证协商交易的综合能源服务商合作运行方法。

5.一种计算机设备，包括储存器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1‑2中任一项所述的一种基于绿证协商交易的综合能源服务商合作运行方法。