CN116822252A - 抽水蓄能电站水能参数的快速确定方法 - Google Patents
抽水蓄能电站水能参数的快速确定方法 Download PDFInfo
- Publication number
- CN116822252A CN116822252A CN202311089919.8A CN202311089919A CN116822252A CN 116822252 A CN116822252 A CN 116822252A CN 202311089919 A CN202311089919 A CN 202311089919A CN 116822252 A CN116822252 A CN 116822252A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- warehouse
- dam
- capacity
- equation
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 121
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000005086 pumping Methods 0.000 claims abstract description 27
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 238000009412 basement excavation Methods 0.000 claims description 63
- 239000000463 material Substances 0.000 claims description 20
- 238000005056 compaction Methods 0.000 claims description 12
- 238000007599 discharging Methods 0.000 claims description 9
- 238000004146 energy storage Methods 0.000 claims description 9
- 238000010276 construction Methods 0.000 claims description 7
- 238000012876 topography Methods 0.000 claims description 7
- 238000010248 power generation Methods 0.000 claims description 4
- 230000037237 body shape Effects 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 description 7
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- -1 geology Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000007620 mathematical function Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Geometry (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Data Mining & Analysis (AREA)
- Structural Engineering (AREA)
- Civil Engineering (AREA)
- Architecture (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Control Of Non-Electrical Variables (AREA)
Abstract
一种抽水蓄能电站水能参数的快速确定方法，涉及抽水蓄能电站工程领域。本发明是基于抽发比、库容和挖填三个平衡条件下，通过选定轴线，建立库容、坝体工程量等相关函数并通过联立五元二次方程方程直接锁定最优解的思路，求得主要水能参数后，即可分别进行上下库的库坝设计，后继的工作与常规设计方法完全相同。通过联立五元二次方程方程直接锁定最优解的思路和方法，可最大程度上减轻设计人员的设计计算工作量，解决了传统设计方法试算工作量巨大且无法追踪到最优解只能得到相对优解的痛难点。
Description
技术领域
本发明涉及抽水蓄能电站工程领域，尤其是一种基于抽发比、库容和挖填三个平衡条件下的抽水蓄能电站水能参数的快速确定方法。
背景技术
随着抽水蓄能电站项目迎来大发展，是对效率、质量和数字化应用的挑战，传统的生产力与生产关系难以适应当前局面。当前，我国水电工程建设已基本形成机械化、信息化、数字化为主体的建设格局，随着云计算、大数据、物联网、人工智能等现代信息技术、数据技术与传感技术的快速发展，正朝着智能化方向迈进。我国抽水蓄能电站工程建设中，智能化设计技术正得到越来越广泛的应用。
目前，抽水蓄能电站相比较于传统水电站主要有以下特点：
（1）抽水蓄能电站的功能决定了必须同时具备上、下两个水库，且上、下水库通过一系列水能参数动态关联，同时与输水发电建筑物、机电设备的设计也密切相关；
（2）抽水蓄能电站的上、下库一般为中小型水库，库容需求不大；且上库多为在山顶洼地、河沟修建，天然库容较小，一般都需要开挖扩库，在增加库容的同时，开挖料还可用做坝体填筑料，可最大程度降低工程投资。
上述特点决定了上下库的设计工作相对于常规水电站，要复杂、困难得多。主要体现在：
（1）抽水蓄能电站的主要水能设计参数数量较多，且动态关联，一般情况下，需要同时满足抽发比限制条件、调节库容和挖填平衡等三个基本条件，才是最为经济合理的。但三个基本条件涉及电力规划、水文泥沙、地质、水工、施工等多个主要专业，专业配合协调难度极大；
（2）找到一组能够同时满足抽发比限制、调节库容需求、挖填基本平衡的水能参数，一般采用试算方法，试算工作量较大，且很难找到准确结果；
（3）抽水蓄能电站库坝址、坝线选择范围较大，一般情况下可选择库坝址、坝线、库形开挖方案的数量均较多，且主要的水能参数均需要经济技术比选后确定，每一次比较的组合方案的数量都是上、下库和库形的乘积，设计和计算工程量巨大，方案筛选十分困难。
发明内容
本发明所要解决的问题就是抽水蓄能电站设计的上述困难，提供一种基于抽发比、库容和挖填三个平衡条件下的抽水蓄能电站水能参数的快速确定方法。
本发明的抽水蓄能电站水能参数的快速确定方法，基于抽发比、库容和挖填量三个平衡条件，其特征在于该确定方法，包括以下步骤：
S1，根据选定的上库和下库的库址以及坝轴线的地形地质情况，依据设计工程师经验和抽水蓄能电站运行要求，确定坝体断面体型和库形开挖体型参数；
S2，设上库天然库容为𝑢𝑅(𝑍 𝑢 )，下库的天然库容为d𝑅(𝑍 d )，上库正常蓄水位记为Z u1 ，上库死水位记为Z u2 ，下库正常蓄水位记为Z d1 ，下库死水位记为Z d2 ，则𝑢𝑅(𝑍 𝑢1 )为上库正常蓄水位对应的天然库容，𝑢𝑅(𝑍 𝑢2 )为上库死水位对应的天然库容，d𝑅(𝑍 d1 )为下库正常蓄水位对应的天然库容，d𝑅(𝑍 d2 )为下库死水位对应的天然库容；
S3，设上库大坝填筑工程量为𝑢𝐵𝑉(𝑍 𝑢1 +Δ𝐻) ，下库大坝填筑工程量为d𝐵𝑉 (𝑍 d1 +Δ𝐻)，上库大坝侵占库容量为𝑢𝑀𝑅(𝑍 𝑢2 ,𝑍 𝑢1 +Δ𝐻），下库大坝侵占库容量为d𝑀 𝑅(𝑍 d2 ,𝑍 d1 +Δ𝐻）；
S4，设上库开挖工程量为𝑢𝐸𝑉(𝑍 𝑢2 )，下库开挖工程量为d𝐸𝑉(𝑍 d2 )，上库开挖扩库量为𝑢𝐴𝑅(𝑍 𝑢1 ,𝑍 𝑢2 )，下库开挖扩库量为d𝐴𝑅(𝑍 d1 ,𝑍 d2 )；
S5，上库和下库所需的的调节库容量相等，设调节库容量为reqR(Z u1 ，Z u2 ，Z d1 ，Z d2 )；
S6，通过抽水蓄能电站抽发比限制条件、库容和挖填平衡条件，分别建立以下方程式：
方程1：抽发比限制条件：抽发比=允许值
；
式中，ℎ 𝑐 为电站抽水时的水头损失，ℎ 𝑓 为电站发电时的水头损失，μ为抽发比的允许值；
方程2：上库调节库容：上库兴利库容+上库开挖扩库量-上库大坝侵占库容量=调节库容
𝑢𝑅(𝑍𝑢1)−𝑢𝑅(𝑍𝑢2)+𝑢𝐴𝑅(𝑍𝑢1,𝑍𝑢2)−𝑢𝑀𝑅(𝑍𝑢2,𝑍𝑢1+Δ𝐻）=𝑟𝑒𝑞𝑅(𝑍𝑢1,𝑍𝑢2,𝑍d1,𝑍d2)；
上库兴利库容为上库正常蓄水位与上库死水位之间的库容差值；
方程3：下库调节库容：下库兴利库容+下库开挖扩库量-下库大坝侵占库容量=调节库容
d𝑅(𝑍d1)−d𝑅(𝑍d2)+d𝐴𝑅(𝑍d1,𝑍d2)−d𝑀𝑅(𝑍d2,𝑍d1+Δ𝐻）=𝑟𝑒𝑞𝑅(𝑍𝑢1,𝑍𝑢2,𝑍d1,𝑍d2)；
下库兴利库容为下库正常蓄水位与下库死水位之间的库容差值；
方程4：上库大坝填筑方量=上库有用料开挖量+上库洞挖填筑方量
𝑢𝐵𝑉(𝑍𝑢1+Δ𝐻)×ξ𝑢=𝑢𝐸𝑉(𝑍𝑢2)×η𝑢+𝑢𝑇𝑉；
式中，𝑢𝑇𝑉为上库洞挖填筑方量，ξ𝑢为上库坝体填筑压实系数，上库大坝填筑工程量与上库坝体填筑压实系数的乘积为上库大坝填筑方量，η 𝑢 为上库有用料开挖量系数，上库开挖工程量与上库有用料开挖量系数的乘积为上库有用料开挖量，上库大坝填筑方量为自然方；
方程5：下库大坝填筑方量=下库有用料开挖量+下库洞挖填筑方量
d𝐵𝑉(𝑍d1+Δ𝐻)×ξd=d𝐸𝑉(𝑍d2)×ηd+d𝑇𝑉；
式中，dTV为下库洞挖填筑方量，ζd为下库大坝填筑压实系数，下库大坝填筑工程量与下库坝体填筑压实系数的乘积为下库大坝填筑方量，η d 为下库有用料开挖量系数，下库开挖工程量与下库有用料开挖量系数的乘积为下库有用料开挖量，下库大坝填筑方量为自然方；
S7，上述方程组，求解过程如下：
S7-1：根据大坝建设要求，设定上库正常蓄水位的初始值，根据坝顶高程、坝体形状以及地形参数，建立大坝填筑三维模型，得到上库大坝填筑方量；
S7-2：通过方程4和上库大坝填筑方量，求出上库死水位；
S7-3：根据设定得上库正常蓄水位初始值和上库死水位，通过方程2，得到上库的调节库容；
S7-4：根据设定的上库正常蓄水位初始值，S7-2得到的上库死水位，通过方程1得到下库正常蓄水位与下库死水位的关系式；
S7-5：根据S7-4得到的下库正常蓄水位和死水位关系式，结合方程5，求得下库正常蓄水位和死水位；
S7-6：带入S7-2至S7-5求得的上库死水位、下库正常蓄水位和下库死水位以及设定的上库正常蓄水位初始值，当数值满足方程3和方程4，则上库和下库特征水位即为所求；当不满足方程3和方程4，进入下一过程；
S7-7：当S7-6中不满足方程3和方程4，采用二分法，设定上库的正常蓄水位的第二个初始值，之后重复S7-1至S7-6，若依然不满足，则设定第三个上库正常蓄水位初始值，进行迭代计算，直至求得满足方程3和方程4的解；
S8，当方程组迭代不收敛，无解时，则可蓄能量超出，调整大坝体型参数和坝轴线，再按S7-1至S7-7步骤重新计算；
S9，根据求得的特征水位参数，进一步复核库坝设计、库容、抽发比以及满足方程4和方程5，既能同时确定抽水蓄能电站的动能参数、坝形参数、工程量参数以及方案可比投资的设计参数。
本发明的抽水蓄能电站水能参数的快速确定方法，通过设定合理的上库大坝蓄水初始值，即可通过方程组求得上下库主要的水能参数，随后即可分别进行上下库的库坝设计，后继的工作与常规设计方法完全相同；通过联立五元二次方程方程快速求解和验证参数的合理性，直接锁定最优解的思路和方法，通过建立多个参数之间的关联性，利用方程组迭代求解的方式，求得多个参数最优匹配方案，得到一组能够同时满足抽发比限制、调节库容需求、挖填基本平衡的水能参数，大大减轻了设计人员的设计计算工作量，解决了传统设计方法试算工作量巨大且无法追踪到最优解只能得到相对优解的痛难点。
附图说明
图1本发明的抽水蓄能电站水能参数的确定方法流程图。
图2抽水蓄能电站挖填平衡示意图。
具体实施方式
实施例1：一种抽水蓄能电站水能参数的快速确定方法，基于抽发比、库容和挖填量三个平衡条件，该确定方法，包括以下步骤：
S1，根据选定的上库和下库的库址以及坝轴线的地形地质情况，依据设计工程师经验和抽水蓄能电站运行要求，确定坝体断面体型和库形开挖体型参数；
S2，设上库天然库容为𝑢𝑅(𝑍 𝑢 )，下库的天然库容为d𝑅(𝑍 d )，上库正常蓄水位记为Z u1 ，上库死水位记为Z u2 ，下库正常蓄水位记为Z d1 ，下库死水位记为Z d2 ，则𝑢𝑅(𝑍 𝑢1 )为上库正常蓄水位对应的天然库容，𝑢𝑅(𝑍 𝑢2 )为上库死水位对应的天然库容，d𝑅(𝑍 d1 )为下库正常蓄水位对应的天然库容，d𝑅(𝑍 d2 )为下库死水位对应的天然库容；
S3，设上库大坝填筑工程量为𝑢𝐵𝑉(𝑍 𝑢1 +Δ𝐻) ，下库大坝填筑工程量为d𝐵𝑉 (𝑍 d1 +Δ𝐻)，上库大坝侵占库容量为𝑢𝑀𝑅(𝑍 𝑢2 ,𝑍 𝑢1 +Δ𝐻），下库大坝侵占库容量为d𝑀 𝑅(𝑍 d2 ,𝑍 d1 +Δ𝐻）；
S4，设上库开挖工程量为𝑢𝐸𝑉(𝑍 𝑢2 )，下库开挖工程量为d𝐸𝑉(𝑍 d2 )，上库开挖扩库量为𝑢𝐴𝑅(𝑍 𝑢1 ,𝑍 𝑢2 )，下库开挖扩库量为d𝐴𝑅(𝑍 d1 ,𝑍 d2 )；
S5，上库和下库所需的的调节库容量相等，设调节库容量为reqR(Z u1 ，Z u2 ，Z d1 ，Z d2 )；
S6，通过抽水蓄能电站抽发比限制条件、库容和挖填平衡条件，分别建立以下方程式：
方程1：抽发比限制条件：抽发比=允许值
；
式中，ℎ 𝑐 为电站抽水时的水头损失，ℎ 𝑓 为电站发电时的水头损失，μ为抽发比的允许值；
方程2：上库调节库容：上库兴利库容+上库开挖扩库量-上库大坝侵占库容量=调节库容
𝑢𝑅(𝑍𝑢1)−𝑢𝑅(𝑍𝑢2)+𝑢𝐴𝑅(𝑍𝑢1,𝑍𝑢2)−𝑢𝑀𝑅(𝑍𝑢2,𝑍𝑢1+Δ𝐻）=𝑟𝑒𝑞𝑅(𝑍𝑢1,𝑍𝑢2,𝑍d1,𝑍d2)；
上库兴利库容为上库正常蓄水位与上库死水位之间的库容差值；
方程3：下库调节库容：下库兴利库容+下库开挖扩库量-下库大坝侵占库容量=调节库容
d𝑅(𝑍d1)−d𝑅(𝑍d2)+d𝐴𝑅(𝑍d1,𝑍d2)−d𝑀𝑅(𝑍d2,𝑍d1+Δ𝐻）=𝑟𝑒𝑞𝑅(𝑍𝑢1,𝑍𝑢2,𝑍d1,𝑍d2)；
下库兴利库容为下库正常蓄水位与下库死水位之间的库容差值；
方程4：上库大坝填筑方量=上库有用料开挖量+上库洞挖填筑方量
𝑢𝐵𝑉(𝑍𝑢1+Δ𝐻)×ξ𝑢=𝑢𝐸𝑉(𝑍𝑢2)×η𝑢+𝑢𝑇𝑉；
式中，𝑢𝑇𝑉为上库洞挖填筑方量，ξ𝑢为上库坝体填筑压实系数，上库大坝填筑工程量与上库坝体填筑压实系数的乘积为上库大坝填筑方量，η 𝑢 为上库有用料开挖量系数，上库开挖工程量与上库有用料开挖量系数的乘积为上库有用料开挖量，上库大坝填筑方量为自然方；
方程5：下库大坝填筑方量=下库有用料开挖量+下库洞挖填筑方量
d𝐵𝑉(𝑍d1+Δ𝐻)×ξd=d𝐸𝑉(𝑍d2)×ηd+d𝑇𝑉；
式中，dTV为下库洞挖填筑方量，ζd为下库大坝填筑压实系数，下库大坝填筑工程量与下库坝体填筑压实系数的乘积为下库大坝填筑方量，η d 为下库有用料开挖量系数，下库开挖工程量与下库有用料开挖量系数的乘积为下库有用料开挖量，下库大坝填筑方量为自然方；
S7，上述方程组，求解过程如下：
S7-1：根据大坝建设要求，设定上库正常蓄水位的初始值，根据坝顶高程、坝体形状以及地形参数，建立大坝填筑三维模型，得到上库大坝填筑方量；
S7-2：通过方程4和上库大坝填筑方量，求出上库死水位；
S7-3：根据设定得上库正常蓄水位初始值和上库死水位，通过方程2，得到上库的调节库容；
S7-4：根据设定的上库正常蓄水位初始值，S7-2得到的上库死水位，通过方程1得到下库正常蓄水位与下库死水位的关系式；
S7-5：根据S7-4得到的下库正常蓄水位和死水位关系式，结合方程5，求得下库正常蓄水位和死水位；
S7-6：带入S7-2至S7-5求得的上库死水位、下库正常蓄水位和下库死水位以及设定的上库正常蓄水位初始值，当数值满足方程3和方程4，则上库和下库特征水位即为所求；当不满足方程3和方程4，进入下一过程；
S7-7：当S7-6中不满足方程3和方程4，采用二分法，设定上库的正常蓄水位的第二个初始值，之后重复S7-1至S7-6，若依然不满足，则设定第三个上库正常蓄水位初始值，进行迭代计算，直至求得满足方程3和方程4的解；
S8，当方程组迭代不收敛，无解时，则可蓄能量超出，调整大坝体型参数和坝轴线，再按S7-1至S7-7步骤重新计算；
S9，根据求得的特征水位参数，进一步复核库坝设计、库容、抽发比以及满足方程4和方程5，既能同时确定抽水蓄能电站的动能参数、坝形参数、工程量参数以及方案可比投资的设计参数。
上述方法中，设定的参数为：
（1）上、下库天然库容参数：𝑢𝑅(𝑍𝑢)和d𝑅(𝑍d)；
（2）调节库容：reqR(Z u1 ,Z u2 ,Z d1 ,Z d2 )；
（3）大坝填筑工程量：上库𝑢𝐵𝑉(𝑍𝑢1+Δ𝐻)和下库d𝐵𝑉(𝑍d1+Δ𝐻)；
（4）库区开挖工程量：上库𝑢𝐸𝑉(𝑍𝑢2)和下库d𝐸𝑉(𝑍d2)；
（5）开挖扩库量：上库𝑢𝐴𝑅(𝑍𝑢1,𝑍𝑢2)和下库d𝐴𝑅(𝑍d1,𝑍d2)；
（6）坝体侵占库容量: 上库𝑢𝑀𝑅(𝑍 𝑢2 ,𝑍 𝑢1 +Δ𝐻）和下库d𝑀𝑅(𝑍 d2 ,𝑍 d1 +Δ 𝐻）；
共计十一个函数，函数特征和运用规律如下：
（1）天然库容参数对于选定的库址和坝线，是固定不变的，相当于一个常量；
（2）调节库容是一个数学函数，大小取决于蓄能量和额定水头，在装机容量和满发小时数给定的情况下，蓄能量是定值，而额定水头则随水位的变化而变化，是因变量；
（3）大坝填筑工程量、坝体侵占库容量，对选定坝轴线和坝体断面，曲线也不会再变化，也相当于一个常量；
（4）库区开挖工程量以及开挖扩库量，对于不同的库区开挖方案有不同的工程量曲线和库容曲线，也就对应不同的特征水位，特征水位和库形开挖方案具有一一对应的关系；
（5）库区开挖方案则变化多端，主要影响因素为库区地形条件、地质情况，开挖的范围、大小、边坡高度以及边坡稳定性等；因此，水位的直接影响因素就是库形开挖方案，水位的比选就是库区开挖方案的比选；
（6）抽发比决定着调节库容的体型，瘦高还是扁平；库容平衡决定着调节库容的大小；挖填平衡决定着调节库容的位置；
（7）上述十一个函数均是单调函数，因此方程组在有解的情况下，是唯一确定的解，无解的情况则是可蓄能量偏大，需要调整可蓄能量数值；
（8）调节库容，也还是可蓄能量的函数；当可蓄能量较小时，显然方程组各式都容易满足，但可蓄能量较大时，就不一定能找到满足各式的解，因此必定存在着一个可蓄能量的最大值，这个最大值即集中反映了十条曲线确定成库条件的优劣，最大值越大，则成库条件越好；在装机容量给定的情况下，就对应存在一个满发小时数的最大值，这个小时数最大值也就集中反映了十条曲线决定的成库条件；当可蓄能量极小时，说明该处并不适合建立抽水蓄能电站，在选址过程中，就已淘汰，因此可以不考虑可蓄能量极小的情况；
（9）对于不同的上、下库成库型式，根据成库特点，灵活调整上述参数和方程组结构，上述计算方法均能适用。
Claims (1)
1.一种抽水蓄能电站水能参数的快速确定方法，基于抽发比、库容和挖填量三个平衡条件，其特征在于该确定方法，包括以下步骤：
S1，根据选定的上库和下库的库址以及坝轴线的地形地质情况，依据设计工程师经验和抽水蓄能电站运行要求，确定坝体断面体型和库形开挖体型参数；
S2，设上库天然库容为𝑢𝑅(𝑍 𝑢 )，下库的天然库容为d𝑅(𝑍 d )，上库正常蓄水位记为Z u1 ，上库死水位记为Z u2 ，下库正常蓄水位记为Z d1 ，下库死水位记为Z d2 ，则𝑢𝑅(𝑍 𝑢1 )为上库正常蓄水位对应的天然库容，𝑢𝑅(𝑍 𝑢2 )为上库死水位对应的天然库容，d𝑅(𝑍 d1 )为下库正常蓄水位对应的天然库容，d𝑅(𝑍 d2 )为下库死水位对应的天然库容；
S3，设上库大坝填筑工程量为𝑢𝐵𝑉(𝑍 𝑢1 +Δ𝐻) ，下库大坝填筑工程量为d𝐵𝑉(𝑍 d1 +Δ𝐻)，上库大坝侵占库容量为𝑢𝑀𝑅(𝑍 𝑢2 ,𝑍 𝑢1 +Δ𝐻），下库大坝侵占库容量为d𝑀𝑅 (𝑍 d2 ,𝑍 d1 +Δ𝐻）；
S4，设上库开挖工程量为𝑢𝐸𝑉(𝑍 𝑢2 )，下库开挖工程量为d𝐸𝑉(𝑍 d2 )，上库开挖扩库量为𝑢𝐴𝑅(𝑍 𝑢1 ,𝑍 𝑢2 )，下库开挖扩库量为d𝐴𝑅(𝑍 d1 ,𝑍 d2 )；
S5，上库和下库所需的的调节库容量相等，设调节库容量为reqR(Z u1 ，Z u2 ，Z d1 ，Z d2 )；
S6，通过抽水蓄能电站抽发比限制条件、库容和挖填平衡条件，分别建立以下方程式：
方程1：抽发比限制条件：抽发比=允许值
；
式中，ℎ 𝑐 为电站抽水时的水头损失，ℎ 𝑓 为电站发电时的水头损失，μ为抽发比的允许值；
方程2：上库调节库容：上库兴利库容+上库开挖扩库量-上库大坝侵占库容量=调节库容
𝑢𝑅(𝑍𝑢1)−𝑢𝑅(𝑍𝑢2)+𝑢𝐴𝑅(𝑍𝑢1,𝑍𝑢2)−𝑢𝑀𝑅(𝑍𝑢2,𝑍𝑢1+Δ𝐻）=𝑟𝑒𝑞𝑅(𝑍𝑢1,𝑍𝑢2,𝑍d1,𝑍d2)；
上库兴利库容为上库正常蓄水位与上库死水位之间的库容差值；
方程3：下库调节库容：下库兴利库容+下库开挖扩库量-下库大坝侵占库容量=调节库容
d𝑅(𝑍d1)−d𝑅(𝑍d2)+d𝐴𝑅(𝑍d1,𝑍d2)−d𝑀𝑅(𝑍d2,𝑍d1+Δ𝐻）=𝑟𝑒𝑞𝑅(𝑍𝑢1,𝑍𝑢2,𝑍d1,𝑍d2)；
方程4：上库大坝填筑方量=上库有用料开挖量+上库洞挖填筑方量
𝑢𝐵𝑉(𝑍𝑢1+Δ𝐻)×ξ𝑢=𝑢𝐸𝑉(𝑍𝑢2)×η𝑢+𝑢𝑇𝑉；
式中，𝑢𝑇𝑉为上库洞挖填筑方量，ξ𝑢为上库坝体填筑压实系数，上库大坝填筑工程量与上库坝体填筑压实系数的乘积为上库大坝填筑方量，η 𝑢 为上库有用料开挖量系数，上库开挖工程量与上库有用料开挖量系数的乘积为上库有用料开挖量，上库大坝填筑方量为自然方；
方程5：下库大坝填筑方量=下库有用料开挖量+下库洞挖填筑方量
d𝐵𝑉(𝑍d1+Δ𝐻)×ξd=d𝐸𝑉(𝑍d2)×ηd+d𝑇𝑉；
式中，dTV为下库洞挖填筑方量，ζd为下库大坝填筑压实系数，下库大坝填筑工程量与下库坝体填筑压实系数的乘积为下库大坝填筑方量，η d 为下库有用料开挖量系数，下库开挖工程量与下库有用料开挖量系数的乘积为下库有用料开挖量，下库大坝填筑方量为自然方；
S7，上述方程组，求解过程如下：
S7-1：根据大坝建设要求，设定上库正常蓄水位的初始值，根据坝顶高程、坝体形状以及地形参数，建立大坝填筑三维模型，得到上库大坝填筑方量；
S7-2：通过方程4和上库大坝填筑方量，求出上库死水位；
S7-3：根据设定得上库正常蓄水位初始值和上库死水位，通过方程2，得到上库的调节库容；
S7-4：根据设定的上库正常蓄水位初始值，S7-2得到的上库死水位，通过方程1得到下库正常蓄水位与下库死水位的关系式；
S7-5：根据S7-4得到的下库正常蓄水位和死水位关系式，结合方程5，求得下库正常蓄水位和死水位；
S7-6：带入S7-2至S7-5求得的上库死水位、下库正常蓄水位和下库死水位以及设定的上库正常蓄水位初始值，当数值满足方程3和方程4，则上库和下库特征水位即为所求；当不满足方程3和方程4，进入下一过程；
S7-7：当S7-6中不满足方程3和方程4，采用二分法，设定上库的正常蓄水位的第二个初始值，之后重复S7-1至S7-6，若依然不满足，则设定第三个上库正常蓄水位初始值，进行迭代计算，直至求得满足方程3和方程4的解；
S8，当方程组迭代不收敛，无解时，则可蓄能量超出，调整大坝体型参数和坝轴线，再按S7-1至S7-7步骤重新计算；
S9，根据求得的特征水位参数，进一步复核库坝设计、库容、抽发比以及满足方程4和方程5，既能同时确定抽水蓄能电站的动能参数、坝形参数、工程量参数以及方案可比投资的设计参数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311089919.8A CN116822252B (zh) | 2023-08-28 | 2023-08-28 | 抽水蓄能电站水能参数的快速确定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311089919.8A CN116822252B (zh) | 2023-08-28 | 2023-08-28 | 抽水蓄能电站水能参数的快速确定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116822252A true CN116822252A (zh) | 2023-09-29 |
| CN116822252B CN116822252B (zh) | 2023-11-10 |
Family
ID=88116920
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311089919.8A Active CN116822252B (zh) | 2023-08-28 | 2023-08-28 | 抽水蓄能电站水能参数的快速确定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116822252B (zh) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN204626355U (zh) * | 2015-01-21 | 2015-09-09 | 深圳中广核工程设计有限公司 | 排水蓄能电站 |
| CN107067339A (zh) * | 2017-01-10 | 2017-08-18 | 中国电建集团北京勘测设计研究院有限公司 | 一种寒冷及严寒地区抽水蓄能电站水库冰冻库容的计算方法 |
| CN108060655A (zh) * | 2017-12-15 | 2018-05-22 | 国家电网公司 | 一种三面环山一面冲沟抽水蓄能电站上水库的构造及其施工方法 |
| CN108425349A (zh) * | 2018-03-13 | 2018-08-21 | 陈德华 | 一种高坝大库抽水蓄能水电站的建造方法 |
| CN112232986A (zh) * | 2020-10-19 | 2021-01-15 | 北京许继电气有限公司 | 抽水蓄能电站发电和抽水能力数据处理方法和*** |
| CN113158286A (zh) * | 2021-01-15 | 2021-07-23 | 中国电建集团华东勘测设计研究院有限公司 | 一种基于最大规模准则的抽水蓄能电站水能参数计算方法 |
| DE102021004099A1 (de) * | 2020-04-30 | 2021-11-11 | Gerhard Luther | Pumpspeicherkraftwerk mit einem von einer Ringstaumauer umschlossenen Speicherbecken |
| CN114908718A (zh) * | 2022-06-22 | 2022-08-16 | 中国电建集团成都勘测设计研究院有限公司 | 利用既有水库的抽水蓄能电站下水库布置结构 |
| CN116151421A (zh) * | 2022-12-01 | 2023-05-23 | 中国电建集团北京勘测设计研究院有限公司 | 抽水蓄能电站选址方法、数据处理终端及可读存储介质 |
-
2023
- 2023-08-28 CN CN202311089919.8A patent/CN116822252B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN204626355U (zh) * | 2015-01-21 | 2015-09-09 | 深圳中广核工程设计有限公司 | 排水蓄能电站 |
| CN107067339A (zh) * | 2017-01-10 | 2017-08-18 | 中国电建集团北京勘测设计研究院有限公司 | 一种寒冷及严寒地区抽水蓄能电站水库冰冻库容的计算方法 |
| CN108060655A (zh) * | 2017-12-15 | 2018-05-22 | 国家电网公司 | 一种三面环山一面冲沟抽水蓄能电站上水库的构造及其施工方法 |
| CN108425349A (zh) * | 2018-03-13 | 2018-08-21 | 陈德华 | 一种高坝大库抽水蓄能水电站的建造方法 |
| DE102021004099A1 (de) * | 2020-04-30 | 2021-11-11 | Gerhard Luther | Pumpspeicherkraftwerk mit einem von einer Ringstaumauer umschlossenen Speicherbecken |
| CN112232986A (zh) * | 2020-10-19 | 2021-01-15 | 北京许继电气有限公司 | 抽水蓄能电站发电和抽水能力数据处理方法和*** |
| CN113158286A (zh) * | 2021-01-15 | 2021-07-23 | 中国电建集团华东勘测设计研究院有限公司 | 一种基于最大规模准则的抽水蓄能电站水能参数计算方法 |
| CN114908718A (zh) * | 2022-06-22 | 2022-08-16 | 中国电建集团成都勘测设计研究院有限公司 | 利用既有水库的抽水蓄能电站下水库布置结构 |
| CN116151421A (zh) * | 2022-12-01 | 2023-05-23 | 中国电建集团北京勘测设计研究院有限公司 | 抽水蓄能电站选址方法、数据处理终端及可读存储介质 |
Non-Patent Citations (1)
| Title |
|---|
| 唐文华等: "潘家口抽水蓄能电站水能参数的设计", 水力发电, no. 8, pages 40 - 44 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116822252B (zh) | 2023-11-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110046400B (zh) | 一种基于h-p型有限元法对混凝土重力坝进行静力特性分析的高精度数值模拟计算方法 | |
| CN103235853B (zh) | 用于采煤沉陷地边采边复的基于土方平衡的基塘布局方法 | |
| CN101709645A (zh) | 一种矿区厚煤层开采沉陷区就地取土复垦方法 | |
| CN116822252B (zh) | 抽水蓄能电站水能参数的快速确定方法 | |
| CN110821453B (zh) | 基于三维地质数值模型的注气油藏开发方案设计方法 | |
| CN104851131A (zh) | 一种局部精度约束的三维地形生成方法 | |
| Xu et al. | Developments and characteristics of pumped storage power station in China | |
| CN110990923A (zh) | 水利水电工程双曲拱坝的快速工程计算方法、***及介质 | |
| CN118013628A (zh) | 一种挖填平衡水库成库设计方法 | |
| CN117113785A (zh) | 一种孔隙水压力测量方法、装置、电子设备及存储介质 | |
| CN104747154B (zh) | 一种利用驱油效率比提高蒸汽驱剩余油研究精度的方法 | |
| CN115099092B (zh) | 基于三维建模的尾矿库渗流计算方法 | |
| CN114021240B (zh) | 基于3dexperience平台的拱坝建基面开挖边坡模型建模方法 | |
| CN112536914B (zh) | 土心墙坝心墙土料掺砾方法 | |
| CN113158597B (zh) | 一种基于catia计算的水闸应力稳定分析方法 | |
| CN108805981A (zh) | 一种分支河道型三角洲前缘训练图像建立方法 | |
| CN1862912A (zh) | 梯级串联电站联合运行控制调节器 | |
| CN110069850B (zh) | 一种基于bim创建堤坝横断面的方法 | |
| CN110826131B (zh) | 大江截流龙口流量、上游水位参数矩阵计算方法、***及介质 | |
| Qian et al. | Theoretical analysis on the influence of the slope angle of plate for the failure mechanism of the concrete expanded-plates pile applied for oceanographic engineering under the horizontal force | |
| CN205134355U (zh) | 用于浅覆盖层风力发电机组的基础结构 | |
| CN111191315B (zh) | 一种基于bim技术的土石堤坝智能化设计方法 | |
| CN111931273B (zh) | 边坡安全系数理论值的快速计算方法 | |
| CN113139298A (zh) | 冰水堆积细粒防渗土料超径块碎石含量的评价方法 | |
| CN205653760U (zh) | 尾水水位调节结构 |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |