CN113282874A - 一种海水淡化负荷参与双边交易的价格区间优化方法 - Google Patents
一种海水淡化负荷参与双边交易的价格区间优化方法 Download PDFInfo
- Publication number
- CN113282874A CN113282874A CN202110546354.6A CN202110546354A CN113282874A CN 113282874 A CN113282874 A CN 113282874A CN 202110546354 A CN202110546354 A CN 202110546354A CN 113282874 A CN113282874 A CN 113282874A
- Authority
- CN
- China
- Prior art keywords
- seawater desalination
- desalination plant
- wind
- water
- price
- 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.)
- Granted
Links
- 238000010612 desalination reaction Methods 0.000 title claims abstract description 180
- 239000013535 sea water Substances 0.000 title claims abstract description 179
- 230000002146 bilateral effect Effects 0.000 title claims abstract description 63
- 238000000034 method Methods 0.000 title claims abstract description 57
- 238000005457 optimization Methods 0.000 title claims abstract description 43
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 108
- 230000005611 electricity Effects 0.000 claims abstract description 45
- 238000005265 energy consumption Methods 0.000 claims abstract description 23
- 238000009826 distribution Methods 0.000 claims abstract description 20
- 230000008901 benefit Effects 0.000 claims abstract description 18
- 238000003860 storage Methods 0.000 claims description 29
- 238000001223 reverse osmosis Methods 0.000 claims description 25
- 238000004364 calculation method Methods 0.000 claims description 24
- 230000008569 process Effects 0.000 claims description 24
- 239000013505 freshwater Substances 0.000 claims description 22
- 239000012528 membrane Substances 0.000 claims description 16
- 230000003204 osmotic effect Effects 0.000 claims description 13
- 239000000126 substance Substances 0.000 claims description 7
- 239000000654 additive Substances 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 5
- 230000004907 flux Effects 0.000 claims description 4
- 238000012423 maintenance Methods 0.000 claims description 4
- 238000010521 absorption reaction Methods 0.000 claims description 3
- 230000000996 additive effect Effects 0.000 claims description 3
- 238000011084 recovery Methods 0.000 claims description 3
- 150000003839 salts Chemical class 0.000 claims description 3
- 239000002699 waste material Substances 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 11
- 238000004519 manufacturing process Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 238000009412 basement excavation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
Images
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 (5)
1.一种海水淡化负荷参与双边交易的价格区间优化方法，其特征在于，包括以下步骤：
步骤1、建立海水淡化厂运行能耗模型；
步骤2、建立双边交易参与主体的收益模型，所述参与主体包括海水淡化厂和风电场；
步骤3、将海水淡化厂每日收益最大作为目标函数，建立海水淡化厂运行优化模型，求解不同弃风电价下，双边交易参与主体的收益及消纳弃风电量；
步骤4、分别向所述双边交易参与主体的收益、消纳弃风电量中引入市场中的利益分配参数、消纳弃风预期范围参数，得到弃风电价优化区间。
2.根据权利要求1所述的一种海水淡化负荷参与双边交易的价格区间优化方法，其特征在于，步骤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为每个时段的时长。
3.根据权利要求1所述的一种海水淡化负荷参与双边交易的价格区间优化方法，其特征在于，步骤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利用的弃风功率，
为电网购电分时电价，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)。
4.根据权利要求1所述的一种海水淡化负荷参与双边交易的价格区间优化方法，其特征在于，步骤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)，建立海水淡化厂运行优化模型如下，并求解不同弃风电价下，双边交易参与主体的收益及消纳弃风电量：
maxf1(V,H,Pwi,P1.Pgrid)
上式中，gi(xi)≥0为公式(20)-公式(24)的不等式约束集合，hj(xj)＝0为公式(25)-公式(27)的等式约束集合。
5.根据权利要求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)。
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 true CN113282874A (zh) | 2021-08-20 |
| CN113282874B 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)
| 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 | 东北大学 | 考虑负荷响应不确定性的多源多荷协调调度方法及系统 |
-
2021
- 2021-05-19 CN CN202110546354.6A patent/CN113282874B/zh active Active
Patent Citations (7)
| 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)
| Title |
|---|
| 孙婧卓等: "大规模电储热参与调峰辅助双边交易容量及价格区间优化模型", 《电网技术》, vol. 43, no. 6, pages 1995 - 2000 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113282874B (zh) | 2024-04-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114362168B (zh) | 一种能源互联系统的设备选型方法 | |
| Spyrou et al. | Design study of a stand-alone desalination system powered by renewable energy sources and a pumped storage unit | |
| CN106981888B (zh) | 基于多源互补的风蓄水火电力系统的多目标动态调度方法 | |
| CN106682810B (zh) | 巨型水电站动态投产下跨流域梯级水电站群长期运行方法 | |
| Yuansheng et al. | A two-stage planning and optimization model for water-hydrogen integrated energy system with isolated grid | |
| CN108229865A (zh) | 一种基于碳交易的电热气综合能源系统低碳经济调度方法 | |
| CN114707403B (zh) | 基于抽水蓄能调节的地区配电网多能源协调优化调度方法 | |
| CN107153975A (zh) | 一种基于联盟博弈的多控制性水库梯级电站补偿效益分摊方法 | |
| CN107895221A (zh) | 市场环境下梯级水电站中长期调度与检修计划优化方法 | |
| CN116169682B (zh) | 一种考虑碳排放流及风光消纳的综合能源系统优化调度策略 | |
| CN106126923B (zh) | 流域梯级龙头水电站效益补偿方法及系统 | |
| CN114243694A (zh) | 考虑阶梯碳交易和需求响应的并网型微电网优化配置方法 | |
| CN102254237A (zh) | 区域水资源供需态势预测及动态调控方法 | |
| CN112365108B (zh) | 一种园区综合能源系统多目标优化协同运行方法 | |
| CN111668878A (zh) | 一种可再生微能源网的优化配置方法和系统 | |
| Desfiandi et al. | Building an energy consumption model and sustainable economic growth in emerging countries | |
| Han et al. | Multi-stage distributionally robust optimization for hybrid energy storage in regional integrated energy system considering robustness and nonanticipativity | |
| CN114862457B (zh) | 梯级水电站群参与日前现货市场的竞价电量边界确定方法 | |
| CN110570033B (zh) | 一种基于合作博弈法的水库多目标优化调度方法 | |
| Yujie et al. | Optimal operation of new coastal power systems with seawater desalination based on grey wolf optimization | |
| CN117526446A (zh) | 梯级水风光多能互补发电系统风光容量双层优化配置方法 | |
| CN109829559A (zh) | 一种促进清洁能源消纳的调峰权交易方法 | |
| CN113282874B (zh) | 一种海水淡化负荷参与双边交易的价格区间优化方法 | |
| CN117094847A (zh) | 一种基于用水户水权交易的灌区水资源优化配置方法 | |
| Pulido-Calvo et al. | Optimal design of pumping stations of inland intensive fishfarms |
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 |