CN113282874B - 一种海水淡化负荷参与双边交易的价格区间优化方法 - Google Patents

一种海水淡化负荷参与双边交易的价格区间优化方法 Download PDF

Info

Publication number
CN113282874B
CN113282874B CN202110546354.6A CN202110546354A CN113282874B CN 113282874 B CN113282874 B CN 113282874B CN 202110546354 A CN202110546354 A CN 202110546354A CN 113282874 B CN113282874 B CN 113282874B
Authority
CN
China
Prior art keywords
sea water
desalination plant
wind
water desalination
power
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110546354.6A
Other languages
English (en)
Other versions
CN113282874A (zh
Inventor
楚帅
张诗钽
葛维春
刘闯
李音璇
崔岱
董辉
唐婧怡
蔡国伟
葛延峰
王顺江
罗桓桓
杨俊友
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Tianjin Electric Power Co Ltd
Shenyang University of Technology
State Grid Liaoning Electric Power Co Ltd
Northeast Electric Power University
Binhai Power Supply Co of State Grid Tianjin Electric Power Co Ltd
Original Assignee
State Grid Tianjin Electric Power Co Ltd
Shenyang University of Technology
Northeast Dianli University
State Grid Liaoning Electric Power Co Ltd
Binhai Power Supply Co of State Grid Tianjin Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Tianjin Electric Power Co Ltd, Shenyang University of Technology, Northeast Dianli University, State Grid Liaoning Electric Power Co Ltd, Binhai Power Supply Co of State Grid Tianjin Electric Power Co Ltd filed Critical State Grid Tianjin Electric Power Co Ltd
Priority to CN202110546354.6A priority Critical patent/CN113282874B/zh
Publication of CN113282874A publication Critical patent/CN113282874A/zh
Application granted granted Critical
Publication of CN113282874B publication Critical patent/CN113282874B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q30/00—Commerce
    • G06Q30/02—Marketing; Price estimation or determination; Fundraising
    • G06Q30/0283—Price estimation or determination
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S50/00—Market activities related to the operation of systems integrating technologies related to power network operation or related to communication or information technologies
    • Y04S50/14—Marketing, i.e. market research and analysis, surveying, promotions, advertising, buyer profiling, customer management or rewards

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Development Economics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Business, Economics & Management (AREA)
  • Marketing (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Accounting & Taxation (AREA)
  • Finance (AREA)
  • Databases & Information Systems (AREA)
  • Water Supply & Treatment (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Game Theory and Decision Science (AREA)
  • Public Health (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Primary Health Care (AREA)
  • Tourism & Hospitality (AREA)
  • Separation Using Semi-Permeable Membranes (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种海水淡化负荷参与双边交易的价格区间优化方法，包括：建立海水淡化厂运行能耗模型；建立双边交易参与主体的收益模型，参与主体包括海水淡化厂和风电场；将海水淡化厂每日收益最大作为目标函数，建立海水淡化厂运行优化模型，求解不同弃风电价下，双边交易参与主体的收益及消纳弃风电量；分别向双边交易参与主体的收益、消纳弃风电量中引入市场中的利益分配参数、消纳弃风预期范围参数，得到弃风电价优化区间。在满足供水需求和海水淡化负荷运行边界等关键参数的前提下，充分尊重双边交易主体的利润分配关系，以消纳清洁能源为目标，优化双边交易的签约价格，得到合理的双边交易的签约价格。

Description

一种海水淡化负荷参与双边交易的价格区间优化方法
技术领域
本发明属于清洁能源消纳方式技术领域，涉及一种海水淡化负荷参与双边交易的价格区间优化方法。
背景技术
为响应中国政府“碳达峰”和“碳中和”目标，中国风电装机容量必然会持续迅猛增加，因调峰困难导致的弃风电量也会随之增长。未来，清洁能源并网规模持续增长的背景下，向负荷侧挖掘调峰资源是一种经济可行的方式。在建设新一代电力系统的大背景下，负荷侧参与电力市场的程度不断加深，负荷不再仅作为电能的消费者，而是向产消者过度。负荷侧多元化的发展趋势为负荷参与主动配电网消纳清洁能源提供了可能性。未来应高度重视负荷侧需求响应，依托市场化手段，引导用户主动提供调峰调频辅助服务、优化用电负荷特性，全面提升系统运行的整体效率。
海水淡化负荷具有可控、可时移、高耗能的特点，可全年运行，不受季节影响，在参与电网调峰、消纳清洁能源方面展现出巨大潜力。随着电力市场制度的逐渐完善，在我国东北地区，可中断负荷已逐步可以参与到调峰辅助服务市场中。考虑中长期电力需求，在电力调度机构的审核与监管下，海水淡化负荷与风电企业签订双边交易合同。在电力市场建设初期，双边交易的签约价格处于探索阶段，只有保证交易价格的合理与公正，才能促进双边交易主体参与辅助服务市场的积极性，保证电力市场建设的长期稳定发展。本发明旨在充分考虑两个交易主体利润分配的前提下，在满足消纳弃风电量的约束下，制定合理的双边交易价格区间。
发明内容
本发明的目的是提供一种海水淡化负荷参与双边交易的价格区间优化方法，能得到合理的双边交易的签约价格。
本发明所采用的技术方案是，一种海水淡化负荷参与双边交易的价格区间优化方法，包括以下步骤：
步骤1、建立海水淡化厂运行能耗模型；
步骤2、建立双边交易参与主体的收益模型，参与主体包括海水淡化厂和风电场；
步骤3、将海水淡化厂每日收益最大作为目标函数，建立海水淡化厂运行优化模型，求解不同弃风电价下，双边交易参与主体的收益及消纳弃风电量；
步骤4、分别向双边交易参与主体的收益、消纳弃风电量中引入市场中的利益分配参数、消纳弃风预期范围参数，得到弃风电价优化区间。
本发明的特点还在于：
步骤1具体包括以下步骤：
步骤1.1、计算在海水淡化过程中维持反渗透过程的最小压力渗透压Δπ：
Δπ＝fosC0/(1-R) (1)；
上式中，fos是海水的渗透压系数，C0是海水的初始盐浓度，R为海水淡化厂的利用率，
R＝V/V0 (2)；
上式中，V为海水淡化厂每时刻的产水量，V0为海水淡化厂每时刻的进水量；
步骤1.2、设海水平均温度为25℃，渗透压系数为73.9Pa/(mg/L)，按照如下方式计算维持最小压力渗透压所需要的能量DSE1：
步骤1.3、计算生产每吨淡水消耗的电能DSE2：
DSE2＝DSE1+2.78×10-7Rmj (4)；
上式中，j为水的通量，Rm是渗透膜的阻抗；
步骤1.4、根据生产每吨淡水消耗的电能DSE2和高压泵的效率，计算海水淡化厂的实际能耗W1：
W1＝DSE2/ηhp (5)；
上式中，ηhp为高压泵运行效率；
步骤1.5、根据实际能耗W1、海水淡化厂在第t个时段的产水量V[t]得到海水淡化厂生产淡水的功率P1：
上式中，ΔT为每个时段的时长。
步骤2具体包括以下步骤：
步骤2.1、海水淡化厂的售水收入R1的计算方式为：
上式中，Sfw为售水单价，Vneed[t]为用户在时段t内的用水需求量，T为供水时间；
步骤2.2、根据海水淡化厂的每日额定产水量Vn计算海水淡化厂的建设成本Cn，并将初始投资成本折算到规划年，得到规划年初始投资的分摊成本C1：
Cn＝6.7×Exp(0.81×lnVn+8.96) (8)；
上式中，CRF为资本回收率，γ为海水淡化厂资金的年利率，n为海水淡化厂生命周期；
步骤2.3、海水淡化厂生产淡水电耗成本Cdian计算过程如下：
上式中，pwc为弃风电价，Pwi[t]为时刻t利用的弃风功率，Set为电网购电分时电价，Pgrid[t]为时刻t的上网电功率，其计算过程如下：
Pgrid[t]＝P1[t]-Pwi[t] (11)；
步骤2.4、设化学添加剂的预处理成本约为0.0225元/m3，化学添加剂的费用OCCH计算过程如下：
OCCH＝0.0225V0 (12)；
步骤2.5、反渗透膜更换每日成本OCME的计算过程如下：
OCME＝PriME×MOD×ζre/365 (13)；
上式中，PriME为反渗透膜的单价，MOD为反渗透系统中包含的反渗透膜数量，ζre为反渗透膜年更换率；
步骤2.6、设备保养成本OCMN的计算过程如下：
OCMN＝OCnom×CoeMN (14)；
上式中，OCnom为常规操作费用，CoeMN为设备保养成本在常规操作费用中占比；
步骤2.7、根据工人的工资、人数计算人工成本OCLB：
OCLB＝OCnom×CoeLB (15)；
上式中，CoeLB为人工成本在常规操作费用中占比；
步骤2.8、结合公式(10)-(15)，计算海水反渗透的运行费用C2：
C2＝Cdian+OCCH+OCME+OCMN+OCLB (16)；
步骤2.9、海水淡化厂每日的最大收益f1计算过程如下：
f1＝R1-C1-C2 (17)；
步骤2.10、风电场通过供给海水淡化厂弃风电的收益为f2计算过程如下：
f2＝(pwc-cw)Qwc (18)；
上式中，cw为风电场的边际成本，Qwc为海水淡化厂消纳弃风电量(kWh)。
步骤3具体包括以下步骤：
步骤3.1、以海水淡化厂每日的收益f1最大为目标，确定目标函数为：
maxf1(V,H,Pwi,P1,Pgrid) (19)；
上式中，V为海水淡化厂的产水量，为储水罐水位高度，Pwi为利用的弃风功率，P1为海水淡化厂生产淡水的功率，Pgrid为上网电功率；
步骤3.2、目标函数的不等式约束条件如下：
0≤Pwi[t]≤min{P1[t],Pwimax[t]} (20)；
Vmin≤V[t]≤Vmax (21)；
Hmin≤H[t]≤Hmax (22)；
0≤P1[t]≤PN (23)；
Vneed[t]≤H[t]At (24)；
上式中，Pwimax[t]为t时刻最大可调度弃风功率，Vmin为海水淡化厂每时刻最小产水量，Vmax为海水淡化厂每时刻最大产水量，Hmin为储水罐最低水位，Hmax为储水罐最高水位，H[t]为t时刻储水罐水位高度，PN为海水淡化厂的额定运行功率，At为储水罐底面积；
步骤3.3、目标函数的等式约束条件如下：
H[0]＝H[24] (25)；
上式中，H(0)第0小时的储水罐水位，H(24)第24小时的储水罐水位，H[t-1]为t-1时刻储水罐水位高度；
步骤3.4、根据公式(19)-(27)，建立海水淡化厂运行优化模型如下，并求解不同弃风电价下，双边交易参与主体的收益及消纳弃风电量：
上式中，gi(xi)≥0为公式(20)-公式(24)的不等式约束集合，hj(xj)＝0为公式(25)-公式(27)的等式约束集合。
步骤4具体包括以下步骤：
步骤4.1、根据步骤3得到海水淡化厂运行优化模型得到不同弃风电价下，双边交易参与主体的收益；向双边交易参与主体的收益中引入市场中的利益分配参数α，对海水淡化厂和风电场的收益进行合理分配，需满足约束条件为：
上式中，[α,1-α]表示市场中利益分配的范围，通过考虑式(29)的约束条件，得出电价区间为[pwcα1,pwcα2]；
步骤4.2、向步骤3得到的消纳弃风电量中引入消纳弃风预期范围参数为β1、β2，则市场中预期弃风消纳范围满足约束条件为：
上式中，β1、β2为弃风消纳预期范围的上下限；
得出预期弃风消纳范围下的电价区间为[pwcβ1,pwcβ2]；
步骤4.3、综合满足步骤4.1-4.2的约束条件，得到弃风电价优化区间为：
[pwc1,pwc2]＝[pwcα1,pwcα2]∩[pwcβ1,pwcβ2] (32)。
本发明的有益效果是：
本发明一种海水淡化负荷参与双边交易的价格区间优化方法，充分考虑反渗透海水淡化技术的非线性运行特性，通过对反渗透海水淡化的能耗与海水淡化厂利用率之间的关系进行建模，在满足供水需求和海水淡化负荷运行边界等关键参数的前提下，充分尊重双边交易主体的利润分配关系，以消纳清洁能源为目标，优化双边交易的签约价格，得到合理的双边交易的签约价格；有效适应中长期风电消纳需求情况下，海水淡化厂与风电企业交易的公正性、持久性，具有较高的实际应用价值。
附图说明
图1是本发明一种海水淡化负荷参与双边交易的价格区间优化方法的流程图；
图2是本发明一种海水淡化负荷参与双边交易的价格区间优化方法的实施例中某日弃风电量分布图；
图3是本发明一种海水淡化负荷参与双边交易的价格区间优化方法实施例中某日淡水需求分布图；
图4是本发明一种海水淡化负荷参与双边交易的价格区间优化方法实施例中电网分时电价图；
图5a是本发明一种海水淡化负荷参与双边交易的价格区间优化方法实施例中弃风电价0.1元/kWh下海水淡化厂用电量与弃风利用情况图；
图5b是本发明一种海水淡化负荷参与双边交易的价格区间优化方法实施例中弃风电价0.2元/kWh下海水淡化厂用电量与弃风利用情况图；
图5c是本发明一种海水淡化负荷参与双边交易的价格区间优化方法实施例中弃风电价0.28元/kWh下海水淡化厂用电量与弃风利用情况图；
图5d是本发明一种海水淡化负荷参与双边交易的价格区间优化方法实施例中弃风电价0.3元/kWh下海水淡化厂用电量与弃风利用情况图；
图6是本发明一种海水淡化负荷参与双边交易的价格区间优化方法实施例中典型弃风电价下各时刻水位高度图；
图7a是本发明一种海水淡化负荷参与双边交易的价格区间优化方法实施例中典型弃风电价下各时刻产水量图；
图7b是本发明一种海水淡化负荷参与双边交易的价格区间优化方法实施例中典型弃风电价下各时刻利用率图；
图8是本发明一种海水淡化负荷参与双边交易的价格区间优化方法实施例中海水淡化厂能耗与弃风电价的关系图；
图9是本发明一种海水淡化负荷参与双边交易的价格区间优化方法实施例中海水淡化厂收益成本及风电场收益分布图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
一种海水淡化负荷参与双边交易的价格区间优化方法，包括以下步骤：
步骤1、根据反渗透海水淡化技术的非线性运行特性，建立海水淡化厂运行能耗模型；
步骤1.1、计算在海水淡化过程中维持反渗透过程的最小压力渗透压Δπ：
Δπ＝fosC0/(1-R) (1)；
上式中，fos是海水的渗透压系数，C0是海水的初始盐浓度，R为海水淡化厂的利用率；
R＝V/V0 (2)；
上式中，V为海水淡化厂每时刻的产水量，V0为海水淡化厂每时刻的进水量；
步骤1.2、设海水平均温度为25℃，渗透压系数为73.9Pa/(mg/L)，按照如下方式计算维持最小压力渗透压所需要的能量DSE1：
步骤1.3、海水淡化厂的实际运行中，除需要保持淡水与海水之间平衡所需的能量(即维持最小压力渗透压所需要的能量DSE1)以外，还需要额外的能量驱动非零的水通量，则通过下式计算生产每吨淡水消耗的电能DSE2：
DSE2＝DSE1+2.78×10-7Rmj (4)；
上式中，j为水的通量，Rm是渗透膜的阻抗；
步骤1.4、在海水淡化过程中，高压泵将电能转化为机械能，推动海水通过渗透膜。所以，在计算海水淡化厂的实际能耗W1时需考虑高压泵的效率。则根据生产每吨淡水消耗的电能DSE2和高压泵的效率，计算海水淡化厂的实际能耗W1：
W1＝DSE2/ηhp (5)；
上式中，ηhp为高压泵运行效率；
步骤1.5、根据实际能耗W1、海水淡化厂在第t个时段的产水量V[t]得到海水淡化厂生产淡水的功率P1：
上式中，ΔT为每个时段的时长。
步骤2、建立双边交易参与主体的收益模型，参与主体包括海水淡化厂和风电场；
步骤2.1、海水淡化厂主要通过出售淡水获得收益，海水淡化厂的售水收入R1的计算方式为：
上式中，Sfw为售水单价，Vneed[t]为用户在时段t内的用水需求量，T为供水时间；
步骤2.2、根据海水淡化厂的每日额定产水量Vn计算海水淡化厂的建设成本Cn，并将初始投资成本折算到规划年，得到规划年初始投资的分摊成本C1：
Cn＝6.7×Exp(0.81×lnVn+8.96) (8)；
上式中，CRF为资本回收率，γ为海水淡化厂资金的年利率，n为海水淡化厂生命周期；
步骤2.3、海水淡化厂生产淡水电耗成本Cdian与海水淡化厂每时刻生产淡水功率P1和电网购电分时电价有关。海水淡化的能耗成本Cdian由购买弃风电量成本与电网购电成本组成，则海水淡化厂生产淡水电耗成本Cdian计算过程如下：
上式中，pwc为弃风电价，Pwi[t]为时刻t利用的弃风功率，Set为电网购电分时电价，Pgrid[t]为时刻t的上网电功率，其计算过程如下：
Pgrid[t]＝P1[t]-Pwi[t] (11)；
步骤2.4、设化学添加剂的预处理成本约为0.0225元/m3，化学添加剂的费用OCCH计算过程如下：
OCCH＝0.0225V0 (12)；
步骤2.5、反渗透膜更换每日成本OCME的计算过程如下：
OCME＝PriME×MOD×ζre/365 (13)；
上式中，PriME为反渗透膜的单价，MOD为反渗透系统中包含的反渗透膜数量，ζre为反渗透膜年更换率；
步骤2.6、设备保养成本OCMN的计算过程如下：
OCMN＝OCnom×CoeMN (14)；
上式中，OCnom为常规操作费用，CoeMN为设备保养成本在常规操作费用中占比；
步骤2.7、根据工人的工资、人数计算人工成本OCLB：
OCLB＝OCnom×CoeLB (15)；
上式中，CoeLB为人工成本在常规操作费用中占比；
步骤2.8、结合公式(10)-(15)，计算海水反渗透的运行费用C2：
C2＝Cdian+OCCH+OCME+OCMN+OCLB (16)；
步骤2.9、海水淡化厂每日的最大收益f1计算过程如下：
f1＝R1-C1-C2 (17)；
步骤2.10、风电场通过供给海水淡化厂弃风电的收益为f2计算过程如下：
f2＝(pwc-cw)Qwc (18)；
上式中，cw为风电场的边际成本，Qwc为海水淡化厂消纳弃风电量(kWh)。
步骤3、将海水淡化厂每日收益最大作为目标函数，建立海水淡化厂运行优化模型，求解不同弃风电价下，双边交易参与主体的收益及消纳弃风电量；
步骤3.1、以海水淡化厂每日的收益f1最大为目标，确定目标函数为：
maxf1(V,H,Pwi,P1,Pgrid) (19)；
上式中，V为海水淡化厂的产水量，为储水罐水位高度，Pwi为利用的弃风功率，P1为海水淡化厂生产淡水的功率，Pgrid为上网电功率；
步骤3.2、为了约束目标函数的变量并满足供水需求要求，该模型应包含海水淡化厂利用弃风约束、产水量约束、储水罐水位高度约束、海水淡化厂实际运行功率约束与供水需求约束。目标函数的不等式约束条件如下：
0≤Pwi[t]≤min{P1[t],Pwimax[t]} (20)；
Vmin≤V[t]≤Vmax (21)；
Hmin≤H[t]≤Hmax (22)；
0≤P1[t]≤PN (23)；
Vneed[t]≤H[t]At (24)；
上式中，Pwimax[t]为t时刻最大可调度弃风功率，Vmin为海水淡化厂每时刻最小产水量，Vmax为海水淡化厂每时刻最大产水量，Hmin为储水罐最低水位，Hmax为储水罐最高水位，H[t]为t时刻储水罐水位高度，PN为海水淡化厂的额定运行功率，At为储水罐底面积；
步骤3.3、为保证海水淡化系统的稳定运行，海水淡化厂需满足储水罐每日的运行约束、产水量与储水罐水位高度关系约束、海水淡化厂消纳弃风需求约束，则目标函数的等式约束条件如下：
H[0]＝H[24] (25)；
上式中，H(0)第0小时的储水罐水位，H(24)第24小时的储水罐水位，H[t-1]为t-1时刻储水罐水位高度；
步骤3.4、根据公式(19)-(27)，建立海水淡化厂运行优化模型如下，并求解不同弃风电价下，双边交易参与主体的收益及消纳弃风电量：
上式中，gi(xi)≥0为公式(20)-公式(24)的不等式约束集合，hj(xj)＝0为公式(25)-公式(27)的等式约束集合。
步骤4、分别向双边交易参与主体的收益、消纳弃风电量中引入市场中的利益分配参数、消纳弃风预期范围参数，得到弃风电价优化区间。
步骤4.1、根据步骤3得到海水淡化厂运行优化模型得到不同弃风电价下，双边交易参与主体的收益；向双边交易参与主体的收益中引入市场中的利益分配参数α，在海水淡化厂和风电场都具有收益的情况下，对海水淡化厂和风电场的收益进行合理分配，需满足约束条件为：
上式中，[α,1-α]表示市场中利益分配的范围，通过考虑式(29)的约束条件，得出电价区间为[pwcα1,pwcα2]；
步骤4.2、为了让海水淡化厂在提供调峰辅助服务时消纳更多的弃风，向步骤3得到的消纳弃风电量中引入消纳弃风预期范围参数为β1、β2，则市场中预期弃风消纳范围满足约束条件为：
上式中，β1、β2为弃风消纳预期范围的上下限；
得出预期弃风消纳范围下的电价区间为[pwcβ1,pwcβ2]；
步骤4.3、综合满足步骤4.1-4.2的约束条件，得到弃风电价优化区间为：
[pwc1,pwc2]＝[pwcα1,pwcα2]∩[pwcβ1,pwcβ2] (32)。
应用时，如图1所示，将弃风电价、分时弃风电量、电网购电电价、淡水需求输入海水淡化厂运行优化模型中，得到海水淡化厂与风电场的收益、消纳弃风电量，根据海水淡化厂与风电场在市场总利益中的占比，引入不同范围参数，对收益进行分配，得到最优弃风电价交易区间。
通过以上方式，一种海水淡化负荷参与双边交易的价格区间优化方法，充分考虑反渗透海水淡化技术的非线性运行特性，通过对反渗透海水淡化的能耗与海水淡化厂利用率之间的关系进行建模，在满足供水需求和海水淡化负荷运行边界等关键参数的前提下，充分尊重双边交易主体的利润分配关系，以消纳清洁能源为目标，优化双边交易的签约价格，得到合理的双边交易的签约价格；有效适应中长期风电消纳需求情况下，海水淡化厂与风电企业交易的公正性、持久性，具有较高的实际应用价值。
实施例
针对某反渗透海水淡化系统，研究本发明的优化策略效果。该省某日的弃风情况如图2所示。海水淡化厂需满足每时刻居民的用水需求。某日淡水需求情况如图3所示。该海水淡化厂的主要参数在表1中给出。电网售电单价如图4所示，用电高峰时段的电价为低谷时段的3倍，海水淡化厂的能耗成本具有显著优化空间。
表1海水淡化厂的主要参数
弃风电价通常偏低，本实施例首先在电网购电最低电价0.2643元/kWh附近，选取[0.1,0.2,0.28,0.3]元/kWh作为典型弃风电价对海水淡化厂进行运行优化。然后，对0.1-0.3元/kWh区间的连续弃风电价进行运行能耗与利润分析，引入不同范围参数对价格区间进行优化。
图5a-5d显示典型弃风电价下，各时段海水淡化厂的用电情况和弃风消纳情况。当弃风电价低于0.2643元/kWh时，海水淡化厂优先使用弃风电量满足产水需求(图5a,图5b)。当弃风电价高于0.2643元/kWh时，海水淡化厂在用电低谷时段全部从电网购电满足产水需求，此时被消纳的弃风电量为0kWh，导致弃风消纳量急剧减少(图5c,图5d)。
图6显示了储水罐水位高度的变化情况。图7a-7b显示了海水淡化厂在一天中的产水量(图7a)和海水淡化厂利用率(图7b)变化情况。海水淡化厂在用电高峰期降低产水量，在平价时期适当提升产水量，优先利用弃风电量，减少电网电量使用，利用储水罐中储备淡水平衡供水需求，获得更大利润(图9)。
图8给出了海水淡化厂在不同弃风电价下消纳弃风电量和从电网购电情况。图中海水淡化厂的产水总能耗随弃风电价的升高缓慢下降。弃风电价较低时，总能耗升高的原因是由于弃风阶段海水淡化厂的利用率较高导致的。尽管较高的海水淡化厂利用率会增加弃风阶段的产水能耗，但是因为弃风电价非常低，一天中的总能耗成本和净利润都是可接受的(图9)。
图9给出了海水淡化厂和风电场的净利润以及海水淡化厂的电费成本。随着弃风电价的升高，海水淡化厂的电费成本提高，海水淡化厂利润减少，海水淡化厂利用弃风电量的经济性降低。当弃风电价低于电网低谷售电电价(0.2643元/kWh)时，风电场净收益随着弃风电价的升高而增长；当弃风电价高于电网低谷售电电价(0.2643元/kWh)时，风电场收益出现骤降现象，这是由于在此阶段下，海水淡化厂为了保证收益最大化，只会在用电高峰和平谷时期消纳弃风电量，风电场出售给海水淡化厂的弃风电量显著降低，从而导致风电场收益骤降。
海水淡化厂与风电场在市场总利益中的占比与弃风利用率情况如表2所示。引入不同范围参数，对弃风电价区间进行优化，结果如表3所示。从表中可知价格区间范围宽度会影响双边交易中交易双方的利益占比，交易价格区间越小，交易双方在市场中利益占比分配的可能性越少，双方利益分配差距缩小。当市场分配利益不变时，弃风利用率不断增加，价格区间会缩小。因此若希望增加海水淡化厂的弃风利用率时，弃风电价应适当减少。
表2海水淡化厂与风电场在市场总利益中的占比
表3不同价格区间优化
/>

Claims (2)

1.一种海水淡化负荷参与双边交易的价格区间优化方法，其特征在于，包括以下步骤：
步骤1、建立海水淡化厂运行能耗模型；
步骤2、建立双边交易参与主体的收益模型，所述参与主体包括海水淡化厂和风电场；
步骤3、将海水淡化厂每日收益最大作为目标函数，建立海水淡化厂运行优化模型，求解不同弃风电价下，双边交易参与主体的收益及消纳弃风电量；
步骤4、分别向所述双边交易参与主体的收益、消纳弃风电量中引入市场中的利益分配参数、消纳弃风预期范围参数，得到弃风电价优化区间；
步骤1具体包括以下步骤：
步骤1.1、计算在海水淡化过程中维持反渗透过程的最小压力渗透压Δπ：
Δπ＝fosC0/(1-R) (1)；
上式中，fos是海水的渗透压系数，C0是海水的初始盐浓度，R为海水淡化厂的利用率，
R＝V/V0 (2)；
上式中，V为海水淡化厂每时刻的产水量，V0为海水淡化厂每时刻的进水量；
步骤1.2、设海水平均温度为25℃，渗透压系数为73.9Pa/(mg/L)，按照如下方式计算维持最小压力渗透压所需要的能量DSE1：
步骤1.3、计算生产每吨淡水消耗的电能DSE2：
DSE2＝DSE1+2.78×10-7Rmj (4)；
上式中，j为水的通量，Rm是渗透膜的阻抗；
步骤1.4、根据生产每吨淡水消耗的电能DSE2和高压泵的效率，计算海水淡化厂的实际能耗W1：
W1＝DSE2/ηhp (5)；
上式中，ηhp为高压泵运行效率；
步骤1.5、根据实际能耗W1、海水淡化厂在第t个时段的产水量V[t]得到海水淡化厂生产淡水的功率P1：
上式中，ΔT为每个时段的时长；
步骤2具体包括以下步骤：
步骤2.1、海水淡化厂的售水收入R1的计算方式为：
上式中，Sfw为售水单价，Vneed[t]为用户在时段t内的用水需求量，T为供水时间；
步骤2.2、根据海水淡化厂的每日额定产水量Vn计算海水淡化厂的建设成本Cn，并将初始投资成本折算到规划年，得到规划年初始投资的分摊成本C1：
Cn＝6.7×Exp(0.81×lnVn+8.96) (8)；
上式中，CRF为资本回收率，γ为海水淡化厂资金的年利率，n为海水淡化厂生命周期；
步骤2.3、海水淡化厂生产淡水电耗成本Cdian计算过程如下：
上式中，pwc为弃风电价，Pwi[t]为时刻t利用的弃风功率，Set为电网购电分时电价，Pgrid[t]为时刻t的上网电功率，其计算过程如下：
Pgrid[t]＝P1[t]-Pwi[t] (11)；
步骤2.4、设化学添加剂的预处理成本约为0.0225元/m3，化学添加剂的费用OCCH计算过程如下：
OCCH＝0.0225V0 (12)；
步骤2.5、反渗透膜更换每日成本OCME的计算过程如下：
OCME＝PriME×MOD×ζre/365 (13)；
上式中，PriME为反渗透膜的单价，MOD为反渗透系统中包含的反渗透膜数量，ζre为反渗透膜年更换率；
步骤2.6、设备保养成本OCMN的计算过程如下：
OCMN＝OCnom×CoeMN (14)；
上式中，OCnom为常规操作费用，CoeMN为设备保养成本在常规操作费用中占比；
步骤2.7、根据工人的工资、人数计算人工成本OCLB：
OCLB＝OCnom×CoeLB (15)；
上式中，CoeLB为人工成本在常规操作费用中占比；
步骤2.8、结合公式(10)-(15)，计算海水反渗透的运行费用C2：
C2＝Cdian+OCCH+OCME+OCMN+OCLB (16)；
步骤2.9、海水淡化厂每日的最大收益f1计算过程如下：
f1＝R1-C1-C2 (17)；
步骤2.10、风电场通过供给海水淡化厂弃风电的收益为f2计算过程如下：
f2＝(pwc-cw)Qwc (18)；
上式中，cw为风电场的边际成本，Qwc为海水淡化厂消纳弃风电量(kWh)；
步骤3具体包括以下步骤：
步骤3.1、以海水淡化厂每日的收益f1最大为目标，确定目标函数为：
maxf1(V,H,Pwi,P1,Pgrid)(19)；
上式中，V为海水淡化厂的产水量，为储水罐水位高度，Pwi为利用的弃风功率，P1为海水淡化厂生产淡水的功率，Pgrid为上网电功率；
步骤3.2、所述目标函数的不等式约束条件如下：
0≤Pwi[t]≤min{P1[t],Pwimax[t]} (20)；
Vmin≤V[t]≤Vmax (21)；
Hmin≤H[t]≤Hmax (22)；
0≤P1[t]≤PN (23)；
Vneed[t]≤H[t]At (24)；
上式中，Pwimax[t]为t时刻最大可调度弃风功率，Vmin为海水淡化厂每时刻最小产水量，Vmax为海水淡化厂每时刻最大产水量，Hmin为储水罐最低水位，Hmax为储水罐最高水位，H[t]为t时刻储水罐水位高度，PN为海水淡化厂的额定运行功率，At为储水罐底面积；
步骤3.3、所述目标函数的等式约束条件如下：
H[0]＝H[24] (25)；
上式中，H(0)第0小时的储水罐水位，H(24)第24小时的储水罐水位，H[t-1]为t-1时刻储水罐水位高度；
步骤3.4、根据公式(19)-(27)，建立海水淡化厂运行优化模型如下，并求解不同弃风电价下，双边交易参与主体的收益及消纳弃风电量：
上式中，gi(xi)≥0为公式(20)-公式(24)的不等式约束集合，hj(xj)＝0为公式(25)-公式(27)的等式约束集合。
2.根据权利要求1所述的一种海水淡化负荷参与双边交易的价格区间优化方法，其特征在于，步骤4具体包括以下步骤：
步骤4.1、根据步骤3得到海水淡化厂运行优化模型得到不同弃风电价下，双边交易参与主体的收益；向双边交易参与主体的收益中引入市场中的利益分配参数α，对海水淡化厂和风电场的收益进行合理分配，需满足约束条件为：
上式中，[α,1-α]表示市场中利益分配的范围，通过考虑式(29)的约束条件，得出电价区间为[pwcα1,pwcα2]；
步骤4.2、向步骤3得到的所述消纳弃风电量中引入消纳弃风预期范围参数为β1、β2，则市场中预期弃风消纳范围满足约束条件为：
上式中，β1、β2为弃风消纳预期范围的上下限；
得出预期弃风消纳范围下的电价区间为[pwcβ1,pwcβ2]；
步骤4.3、综合满足步骤4.1-4.2的约束条件，得到弃风电价优化区间为：
[pwc1,pwc2]＝[pwcα1,pwcα2]∩[pwcβ1,pwcβ2] (32)。
CN202110546354.6A 2021-05-19 2021-05-19 一种海水淡化负荷参与双边交易的价格区间优化方法 Active CN113282874B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110546354.6A CN113282874B (zh) 2021-05-19 2021-05-19 一种海水淡化负荷参与双边交易的价格区间优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110546354.6A CN113282874B (zh) 2021-05-19 2021-05-19 一种海水淡化负荷参与双边交易的价格区间优化方法

Publications (2)

Publication Number Publication Date
CN113282874A CN113282874A (zh) 2021-08-20
CN113282874B true CN113282874B (zh) 2024-04-12

Family

ID=77279940

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110546354.6A Active CN113282874B (zh) 2021-05-19 2021-05-19 一种海水淡化负荷参与双边交易的价格区间优化方法

Country Status (1)

Country Link
CN (1) CN113282874B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101704690A (zh) * 2009-10-27 2010-05-12 窦观一 利用风电场和非粮田能源植物生产的肥料/饲料
CN108460482A (zh) * 2018-01-31 2018-08-28 南京理工大学 考虑可中断负荷和日前小时电价的风电消纳方法
CN109242350A (zh) * 2018-10-17 2019-01-18 燕山大学 计及可平移负荷的冷热电联供系统容量优化配置方法
CN109447492A (zh) * 2018-11-05 2019-03-08 湘潭大学 一种热电厂参与风电供暖的博弈机制构建方法
CN111382901A (zh) * 2020-02-25 2020-07-07 沈阳工业大学 一种反渗透海水淡化厂的建模方法
CN111882131A (zh) * 2020-07-31 2020-11-03 南方电网调峰调频发电有限公司 一种基于双工况最优效率跟踪的变速海水抽蓄电站容量优化方法
CN111969655A (zh) * 2020-08-17 2020-11-20 东北大学 考虑负荷响应不确定性的多源多荷协调调度方法及系统

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101704690A (zh) * 2009-10-27 2010-05-12 窦观一 利用风电场和非粮田能源植物生产的肥料/饲料
CN108460482A (zh) * 2018-01-31 2018-08-28 南京理工大学 考虑可中断负荷和日前小时电价的风电消纳方法
CN109242350A (zh) * 2018-10-17 2019-01-18 燕山大学 计及可平移负荷的冷热电联供系统容量优化配置方法
CN109447492A (zh) * 2018-11-05 2019-03-08 湘潭大学 一种热电厂参与风电供暖的博弈机制构建方法
CN111382901A (zh) * 2020-02-25 2020-07-07 沈阳工业大学 一种反渗透海水淡化厂的建模方法
CN111882131A (zh) * 2020-07-31 2020-11-03 南方电网调峰调频发电有限公司 一种基于双工况最优效率跟踪的变速海水抽蓄电站容量优化方法
CN111969655A (zh) * 2020-08-17 2020-11-20 东北大学 考虑负荷响应不确定性的多源多荷协调调度方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
大规模电储热参与调峰辅助双边交易容量及价格区间优化模型;孙婧卓等;《电网技术》;第第43卷卷(第第6期期);正文第1995-2000页 *

Also Published As

Publication number Publication date
CN113282874A (zh) 2021-08-20

Similar Documents

Publication Publication Date Title
Spyrou et al. Design study of a stand-alone desalination system powered by renewable energy sources and a pumped storage unit
Divakar et al. Optimal allocation of bulk water supplies to competing use sectors based on economic criterion–An application to the Chao Phraya River Basin, Thailand
Kaldellis et al. The water shortage problem in the Aegean archipelago islands: cost-effective desalination prospects
WO2019006733A1 (zh) 一种跨省互联水电站群长期联合调峰调度方法
CN108229865A (zh) 一种基于碳交易的电热气综合能源系统低碳经济调度方法
Yuansheng et al. A two-stage planning and optimization model for water-hydrogen integrated energy system with isolated grid
CN107895221A (zh) 市场环境下梯级水电站中长期调度与检修计划优化方法
CN104915790A (zh) 一种促进风电消纳的峰谷电价优化方法
CN102254237A (zh) 区域水资源供需态势预测及动态调控方法
CN107153975A (zh) 一种基于联盟博弈的多控制性水库梯级电站补偿效益分摊方法
Ma et al. Optimal capacity configuration of hydro-wind-PV hybrid system and its coordinative operation rules considering the UHV transmission and reservoir operation requirements
CN107133840A (zh) 一种面向电供热取暖促进风电就地消纳的多风场竞价方法
Desfiandi et al. Building an energy consumption model and sustainable economic growth in emerging countries
Howe The coming conflicts over water
CN112365108B (zh) 一种园区综合能源系统多目标优化协同运行方法
CN113282874B (zh) 一种海水淡化负荷参与双边交易的价格区间优化方法
CN105809281A (zh) 一种考虑多业主新增效益分配的水库群调度方法
CN110570033B (zh) 一种基于合作博弈法的水库多目标优化调度方法
Pulido-Calvo et al. Optimal design of pumping stations of inland intensive fishfarms
CN114548609B (zh) 一种面向生态发电供水的水资源协同调配方法
Srinivasan The Indian solar photovoltaic industry: a life cycle analysis
Roumasset et al. Ordering renewable resources: Groundwater, recycling, and desalination
Mu et al. A food security approach to analyse irrigation efficiency improvement demands at the country level
Roumasset et al. The economics of groundwater
Meier Mini hydropower for rural development: a new market-oriented approach to maximize electrification benefits with special focus on Indonesia

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20240313

Address after: No. 394, Yingkou Road, Tanggu, Binhai New Area, Tianjin 300450

Applicant after: STATE GRID TIANJIN ELECTRIC POWER COMPANY BINHAI POWER SUPPLY BRANCH

Country or region after: China

Applicant after: STATE GRID TIANJIN ELECTRIC POWER Co.

Applicant after: SHENYANG University OF TECHNOLOGY

Applicant after: STATE GRID LIAONING ELECTRIC POWER SUPPLY Co.,Ltd.

Applicant after: NORTHEAST DIANLI University

Address before: 110870 No. 111 Shenyang West Road, Shenyang economic and Technological Development Zone, Liaoning

Applicant before: SHENYANG University OF TECHNOLOGY

Country or region before: China

Applicant before: STATE GRID LIAONING ELECTRIC POWER SUPPLY Co.,Ltd.

Applicant before: NORTHEAST DIANLI University

GR01 Patent grant
GR01 Patent grant