CN110222362B - 一种通过隧洞和调压井的多机形式建立一管多机微分方程模型的方法 - Google Patents

一种通过隧洞和调压井的多机形式建立一管多机微分方程模型的方法 Download PDF

Info

Publication number
CN110222362B
CN110222362B CN201910343377.XA CN201910343377A CN110222362B CN 110222362 B CN110222362 B CN 110222362B CN 201910343377 A CN201910343377 A CN 201910343377A CN 110222362 B CN110222362 B CN 110222362B
Authority
CN
China
Prior art keywords
pressure regulating
ith
flow
tunnel
machine
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
CN201910343377.XA
Other languages
English (en)
Other versions
CN110222362A (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.)
Kunming University of Science and Technology
Original Assignee
Kunming University of Science and Technology
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 Kunming University of Science and Technology filed Critical Kunming University of Science and Technology
Priority to CN201910343377.XA priority Critical patent/CN110222362B/zh
Publication of CN110222362A publication Critical patent/CN110222362A/zh
Application granted granted Critical
Publication of CN110222362B publication Critical patent/CN110222362B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation
    • 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
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/20—Hydro energy

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Control Of Water Turbines (AREA)

Abstract

本发明涉及一种通过隧洞和调压井的多机形式建立一管多机微分方程模型的方法，属于水轮机及水力机组稳定性分析和控制技术领域。在调压井之后有n路钢管，每路钢管末端有m根分岔管，机组总的台数为n×m台；水力系统动态包括隧洞、调压井、共用管道和分岔管四部分；以水轮机额定流量Qr和额定水头Hr为基值，Qr单位米3/秒，Hr单位米，具体步骤如下：步骤一：将隧洞和调压井水力动态分解为多机形式；步骤二：构建一管多机带调压井的微分代数多机模型。本方法为研究复杂水力系统和水电站水机电多场耦合条件下的稳定性，以及水电站多机协同非线性控制设计提供一种便捷的计算方法和手段。

Description

一种通过隧洞和调压井的多机形式建立一管多机微分方程模型的方法
技术领域
本发明涉及一种通过隧洞和调压井的多机形式建立一管多机微分方程模型的方法，属于水轮机及水力机组稳定性分析和控制技术领域。
背景技术
在具有长引水系统的水电站，通常采用隧洞引水引水至靠近厂房的调压井，调压井之后采用钢管引水至厂房前，然后分叉出几路管道至水轮机，即一管多机带调压井的复杂引水系统。这种复杂引水系统共用管道是水电站多机之间的水力耦合渠道，其水力动态的计算是研究水电机组稳定性和控制设计的主要困难之一。
尽管经典的水击特征线方法对这种复杂水力系统进行计算，但是，由于水击特征线方法耗时太长，在水电机组的控制设计和稳定性分析中很少使用。传统的应用中主要采用水力暂态的传递函数形式。随着非线性理论的发展，非线性分析和控制理论中一般采用一阶微分方程形式，这种传递函数形式的水力暂态应用很不方便。目前，在设计非线性理论的应用中，水力系统动态一般采用最简单的一阶微分方程形式，即刚性水击的简单情况。近年来，一些学者已建立了考虑弹性水击的水力系统微分方程模型，而且也逐步扩展到带有共用管的一管多机水力系统的微分方程形式。随着研究的深入，考虑水电站水机电多场耦合条件下，水电机组的稳定性分析和协同控制设计成为技术发展的必然。为满足这种应用需求，复杂水力系统必须分解为多机形式，才能将水轮机、发电机联合构成多机形式的模型进行分析研究。
发明内容
针对上述现有技术存在的问题及不足，本发明提供一种通过隧洞和调压井的多机形式建立一管多机微分方程模型的方法。本方法为研究复杂水力系统和水电站水机电多场耦合条件下的稳定性，以及水电站多机协同非线性控制设计提供一种便捷的计算方法和手段。本发明通过以下技术方案实现。
一种通过隧洞和调压井的多机形式建立一管多机微分方程模型的方法，如图1所示，在调压井之后有n路钢管，每路钢管末端有m根分岔管，机组总的台数为n×m台；水力系统动态包括隧洞、调压井、共用管道和分岔管四部分；以水轮机额定流量Qr和额定水头Hr为基值，Qr单位米3/秒，Hr单位米，具体步骤如下：
步骤一：将隧洞和调压井水力动态分解为多机形式；
以水头为和流量的构成为核心，经过分析推导后得到以下两个微分方程：
根据刚性水击传递函数形式，可以直接写出隧洞水力动态方程如下：
Figure BDA0002041455910000021
Figure BDA0002041455910000022
式中，hs是调压井处的水头相对值，即hs＝Hs/Hr，Hs是调压井水头(米)，h0是水轮机静水头(上下游水位差)相对值，即h0＝H0/Hr，H0是水轮机静水头(米)，fpT是隧洞内的擦损失水头系数，qT是隧洞流量的相对值，即qT＝QT/Qr，QT是隧洞流量(米3/秒)，ΔqT是隧洞流量增量相对值，即ΔqT＝ΔQT/Qr，ΔQT是隧洞流量增量(米3/秒)，ΔhT是隧洞流量变化在隧洞出口断面引起的水头动态变化相对值，即ΔhT＝ΔHT/Hr，ΔHT是隧洞水头增量(米)，TwT是隧洞段的水流惯性时间常数。上述t为时间，秒；
这里只考虑最简单的直筒式调压井，忽略调压井摩擦损失，其动态描述方程为：
Figure BDA0002041455910000023
式中，Δhs是调压井水头增量相对值，即Δhs＝ΔHs/Hr，ΔHs是调压井水头增量(米)，Δqs是调压井流量增量相对值，即Δqs＝ΔQs/Qr，ΔQs是调压井流量增量(米3/秒)，Cs是调压井的储能常数(秒)，定义为：Cs＝AsHr/Qr，As是调压井断面面积(米2)。
在调压井处的流量平衡方程为：
Figure BDA0002041455910000024
令：
Figure BDA0002041455910000025
式中，qs是调压井流量相对值，即qs＝Qs/Qr，Qs是调压井流量(米3/秒)，q(i)是连接水轮机的第i路分岔管流量相对值，即q(i)＝Q(i)/Qr，Q(i)是第i路分叉管增量(米3/秒)，ΔqT(i)表示第i台机组流量变化在隧洞中引起的流量变化相对值，即ΔqT(i)＝ΔQT(i)/Qr，ΔQT(i)是第i台机组流量变化在隧洞中引起的流量变化(米3/秒)。
式(1)改写为：
hs＝hs0+Δhs   (6)
即：
Figure BDA0002041455910000031
Δhs＝-fpT(2qT0+ΔqT)ΔqT+ΔhT   (8)
式中，hs0是调压井的稳态水头相对值，即hs0＝Hs0/Hr，Hs0是调压井稳态水头(米)，qT0是隧洞的稳态流量相对值，即qT0＝QT0/Qr，QT0是隧洞稳态流量(米3/秒)，Δhs是调压井水头增量相对值，即Δhs＝ΔHs/Hr，ΔHs是调压井水头增量(米)。
式(7)实际上是式(1)的稳态形式，稳态工况下暂态水头ΔhT＝0。
忽略该式(8)中流量增量的二次项，上式进一步改写为：
Figure BDA0002041455910000032
利用上式(9)、式(5)、和式(2)，有：
Figure BDA0002041455910000033
在调压井处的流量平衡方程式(4.6-6)改写为：
Figure BDA0002041455910000034
若令：
Figure BDA0002041455910000035
则有：
Δqs(i)＝ΔqT(i)-Δq(i)   (12)
式中，Δqs(i)是第i路分岔管流量变化在调压井中引起的流量变化相对值，即Δqs(i)＝ΔQs(i)/Qr，ΔQs(i)是第i路分岔管流量变化在调压井中引起的流量变化(米3/秒)，Δhs(i)是第i路分岔管流量变化在调压井中引起的水头变化相对值，即Δhs(i)＝ΔHs(i)/Hr，ΔHs(i)是第i路分岔管流量变化在调压井中引起的水头变化(米)。
上式表明，单一机组流量变化在调压井断面处也满足流量连续性方程。
假设暂态过程中，调压井水头的增量变化相对值可视为各分岔管流量变化在调压井中引起的水头变化相对值之和，即
Figure BDA0002041455910000041
代入式(3)，对比得到：
Figure BDA0002041455910000042
Δq(i)为第i台水轮机流量增量相对值，即Δq(i)＝ΔQ(i)/Qr，ΔQ(i)是第i台机组流量增量，米3/秒；上述式(10)和式(13)已变为多机形式。
步骤二：构建一管多机带调压井的微分代数多机模型；
将步骤一得到一管多机带调压井的微分代数多机模型与一管多机微分方程模型进行连接，选取状态变量x1(i)，x2(i)，x3(i)，x4(i)＝q(i)，x5(i)＝y(i)，x6(i)＝ΔqT(i),，x7＝Δhs(i)，得到统一的微分代数模型如下：
微分方程模型：
Figure BDA0002041455910000043
Figure BDA0002041455910000044
代数方程模型：
Figure BDA0002041455910000045
式中，x1(i)，x2(i)，x3(i)是与水轮机相连的第i路分叉管动态中间状态变量，无确定的物理意义，下标“(i)”表示与水轮机相连的第i路分叉管的参数，Te(i)是第i路分叉管的弹性时间，秒，TeC是共用管的弹性时间，秒，Zn(i)是第i路分叉管的水力涌浪阻抗的规格化值，ZnC是共用管的水力涌浪阻抗的规格化值，fpC是连接共用管的摩擦损失系数，fp(i)是第i路分叉管的摩擦损失系数，u(i)是第i台水轮机的输入控制，y(i)是第i台水轮机主接力器位移相对值，即y(i)＝Y(i)/Ymax，Y(i)是第i台机主接力器位移，毫米，Ymax是主接力器位移最大值，毫米；y0(i)是第i台机主接力器位移初值相对值，即y0(i)＝Y0(i)/Ymax，Y0(i)是第i台机主接力器位移初值，毫米，yr(i)是第i台机额定工况下主接力器位移相对值，即yr(i)＝Yr(i)/Ymax，Yr(i)是第i台机额定工况下的导叶开度，毫米；Ty(i)是第i台机组主接力器时间常数，秒，pt(i)＝Pt(i)/Pr(i)是第i台水轮机的功率相对值，Pt(i)是第i台水轮机的功率，千瓦，Pr(i)是第i台水轮机的额定功率，千瓦，At(i)是第i台水轮机的增益系数，qnl(i)＝Qnl(i)/Qr是第i台水轮机的空载流量相对值，Qnl(i)是第i台水轮机的空载流量，米3/秒，ht(i)＝Ht(i)/Hr是第i台水轮机的水头相对值，Ht(i)是第i台水轮机的水头，米。
上述公式中相同标号代表相同的含义。或者为本领域技术人员公知的含义。
本发明的有益效果是：
1、将隧洞和调压井水力动态解耦为单机形式，解决了水电站复杂水力系统的多机建模问题。借助这种水力系统的解耦形式，可建立具有复杂引水系统水电站的多机微分代数模型，是水电站局部多机系统多控制器协同设计和稳定性分析的理论基础。
2、隧洞和调压井水力动态的多机分解形式，揭示了水力系统满足叠加原理，属于线性系统的本质特性。
3、本发明所提出隧洞和调压井水力暂态多机形式，也是一种水力耦合的解耦方法，对复杂水力系统的研究具有重要的参考价值。
附图说明
图1是本发明一管多机带调压井的水力系统示意图；
图2是本发明实施例1号机组负荷调节暂态过程中调压井水位和4台水轮机水头的暂态变化。
具体实施方式
下面结合附图和具体实施方式，对本发明作进一步说明。
实施例1
某电站水力系统布置及基本参数为：水为进水口至调压井为引水隧道，其直径8米、长度932米；调压井高65.8米、内径13米；调压井后两根直径4.6米的压力钢管，第一根的长度为517米、第二根的长度490米，其末端又各自分为两根直径由进口处的3.2米过渡到2.2米的岔管连接4台水轮机。仿真中编号1号和2号机组共用一根共管，3号、4号机组共用另一根共管。
水轮机参数:Hr＝312m，Qr＝53.5m3/s，Pr＝150MW，nr＝333.3rpm，隧道水流惯性时间常数Tw＝3.24376(s)，调压井参数：Cs＝774.06(s)。
水轮机调速器采用典型的并联PID结构，控制参数为：Kp＝5.0，KD＝1.5，KI＝2.5，bp＝0.04，控制器执行周期40ms。励磁控制器采用无功的PI控制，KPI＝1.0，KII＝1.5，控制器执行周期20ms。4台机组具有相同的特征及参数。
仿真工况：电站4台机组带额定负荷运行。设定在t＝1s时刻，1号机给定负荷pt(1)＝0.5pr，其余机组负荷保持不变。
该通过隧洞和调压井的多机形式建立一管多机微分方程模型的方法，具体步骤如下：
步骤一：将隧洞和调压井水力动态分解为多机形式；
以水头为和流量的构成为核心，经过分析推导后得到以下两个微分方程：
以第i台水轮机流量相对值表示的隧洞流量增量动态方程为：
Figure BDA0002041455910000061
其中ΔqT是隧洞流量增量相对值，即ΔqT＝ΔQT/Qr，ΔQT是隧洞流量增量，米3/秒；t为时间，秒；TwT是隧洞段的水流惯性时间常数；fpT是隧洞内的擦损失水头系数；qT0是隧洞的稳态流量相对值，即qT0＝QT0/Qr，QT0是隧洞稳态流量，米3/秒；ΔqT(i)表示第i台机组流量变化在隧洞中引起的流量变化相对值，即ΔqT(i)＝ΔQT(i)/Qr，ΔQT(i)是第i台机组流量变化在隧洞中引起的流量变化，米3/秒；Δhs(i)是第i路分岔管流量变化在调压井中引起的水头变化相对值，即Δhs(i)＝ΔHs(i)/Hr，ΔHs(i)是第i路分岔管流量变化在调压井中引起的水头变化，米；
以第i台水轮机流量相对值表示的调压井水位增量动态方程为：
Figure BDA0002041455910000062
其中Δhs(i)是第i路分岔管流量变化在调压井中引起的水头变化相对值，即Δhs(i)＝ΔHs(i)/Hr，ΔHs(i)是第i路分岔管流量变化在调压井中引起的水头变化，米；t为时间，秒；Cs是调压井的储能常数(秒)，定义为：Cs＝AsHr/Qr，As是调压井断面面积，米2；ΔqT(i)表示第i台机组流量变化在隧洞中引起的流量变化相对值，即ΔqT(i)＝ΔQT(i)/Qr，ΔQT(i)是第i台机组流量变化在隧洞中引起的流量变化，米3/秒；Δq(i)为第i台水轮机流量增量相对值，即Δq(i)＝ΔQ(i)/Qr，ΔQ(i)是第i台机组流量增量，米3/秒；
步骤二：构建一管多机带调压井的微分代数多机模型；
将步骤一得到一管多机带调压井的微分代数多机模型与一管多机微分方程模型进行连接，选取状态变量x1(i)，x2(i)，x3(i)，x4(i)＝q(i)，x5(i)＝y(i)，x6(i)＝ΔqT(i),，x7＝Δhs(i)，得到统一的微分代数模型如下：
微分方程模型：
Figure BDA0002041455910000071
Figure BDA0002041455910000072
代数方程模型：
Figure BDA0002041455910000073
式中，x1(i)，x2(i)，x3(i)是与水轮机相连的第i路分叉管动态中间状态变量，无确定的物理意义，下标“(i)”表示与水轮机相连的第i路分叉管的参数，Te(i)是第i路分叉管的弹性时间，秒，TeC是共用管的弹性时间，秒，Zn(i)是第i路分叉管的水力涌浪阻抗的规格化值，ZnC是共用管的水力涌浪阻抗的规格化值，fpC是连接共用管的摩擦损失系数，fp(i)是第i路分叉管的摩擦损失系数，u(i)是第i台水轮机的输入控制，y(i)是第i台水轮机主接力器位移相对值，即y(i)＝Y(i)/Ymax，Y(i)是第i台机主接力器位移，毫米，Ymax是主接力器位移最大值，毫米；y0(i)是第i台机主接力器位移初值相对值，即y0(i)＝Y0(i)/Ymax，Y0(i)是第i台机主接力器位移初值，毫米，yr(i)是第i台机额定工况下主接力器位移相对值，即yr(i)＝Yr(i)/Ymax，Yr(i)是第i台机额定工况下的导叶开度，毫米；Ty(i)是第i台机组主接力器时间常数，秒，pt(i)＝Pt(i)/Pr(i)是第i台水轮机的功率相对值，Pt(i)是第i台水轮机的功率，千瓦，Pr(i)是第i台水轮机的额定功率，千瓦，At(i)是第i台水轮机的增益系数，qnl(i)＝Qnl(i)/Qr是第i台水轮机的空载流量相对值，Qnl(i)是第i台水轮机的空载流量，米3/秒，ht(i)＝Ht(i)/Hr是第i台水轮机的水头相对值，Ht(i)是第i台水轮机的水头，米。
通过本实施例构建一管多机带调压井的微分代数多机模型，获得暂态过程中调压井水位、四台机组水轮机水头变化图，如图2所示。从图2中可以看出，粗线是调压井水头变化相对值；由于1号机组和2号机组有共用管一，1号机组的调节对2号机组水轮机水头的影响较大，使得1号机和2号机水轮机水头具有显示的变化趋势；3号机和4号机拥有相同的共用管，受1号机调节的影响很小，其水头变化主要受调压井水位变化的影响。仿真结果与工程实践是一致的。
以上结合附图对本发明的具体实施方式作了详细说明，但是本发明并不限于上述实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下作出各种变化。

Claims (1)

1.一种通过隧洞和调压井的多机形式建立一管多机微分方程模型的方法，其特征在于：在调压井之后有n路钢管，每路钢管末端有m根分岔管，机组总的台数为n×m台；水力系统动态包括隧洞、调压井、共用管道和分岔管四部分；以水轮机额定流量Qr和额定水头Hr为基值，Qr单位米3/秒，Hr单位米，具体步骤如下：
步骤一：将隧洞和调压井水力动态分解为多机形式；
以水头为和流量的构成为核心，经过分析推导后得到以下两个微分方程：
以第i台水轮机流量相对值表示的隧洞流量增量动态方程为：
其中ΔqT是隧洞流量增量相对值，即ΔqT＝ΔQT/Qr，ΔQT是隧洞流量增量，米3/秒；t为时间，秒；TwT是隧洞段的水流惯性时间常数；fpT是隧洞内的擦损失水头系数；qT0是隧洞的稳态流量相对值，即qT0＝QT0/Qr，QT0是隧洞稳态流量，米3/秒；ΔqT(i)表示第i台机组流量变化在隧洞中引起的流量变化相对值，即ΔqT(i)＝ΔQT(i)/Qr，ΔQT(i)是第i台机组流量变化在隧洞中引起的流量变化，米3/秒；Δhs(i)是第i路分岔管流量变化在调压井中引起的水头变化相对值，即Δhs(i)＝ΔHs(i)/Hr，ΔHs(i)是第i路分岔管流量变化在调压井中引起的水头变化，米；
以第i台水轮机流量相对值表示的调压井水位增量动态方程为：
其中Δhs(i)是第i路分岔管流量变化在调压井中引起的水头变化相对值，即Δhs(i)＝ΔHs(i)/Hr，ΔHs(i)是第i路分岔管流量变化在调压井中引起的水头变化，米；t为时间，秒；Cs是调压井的储能常数,秒，定义为：Cs＝AsHr/Qr，As是调压井断面面积，米2；ΔqT(i)表示第i台机组流量变化在隧洞中引起的流量变化相对值，即ΔqT(i)＝ΔQT(i)/Qr，ΔQT(i)是第i台机组流量变化在隧洞中引起的流量变化，米3/秒；Δq(i)为第i台水轮机流量增量相对值，即Δq(i)＝ΔQ(i)/Qr，ΔQ(i)是第i台机组流量增量，米3/秒；
步骤二：构建一管多机带调压井的微分代数多机模型；
将步骤一得到一管多机带调压井的微分代数多机模型与一管多机微分方程模型进行连接，选取状态变量x1(i)，x2(i)，x3(i)，x4(i)＝q(i)，x5(i)＝y(i)，x6(i)＝ΔqT(i)，x7＝Δhs(i)，得到统一的微分代数模型如下：
微分方程模型：
代数方程模型：
式中，x1(i)，x2(i)，x3(i)是与水轮机相连的第i路分叉管动态中间状态变量，无确定的物理意义，下标“(i)”表示与水轮机相连的第i路分叉管的参数，Te(i)是第i路分叉管的弹性时间，秒，TeC是共用管的弹性时间，秒，Zn(i)是第i路分叉管的水力涌浪阻抗的规格化值，ZnC是共用管的水力涌浪阻抗的规格化值，fpC是连接共用管的摩擦损失系数，fp(i)是第i路分叉管的摩擦损失系数，u(i)是第i台水轮机的输入控制，y(i)是第i台水轮机主接力器位移相对值，即y(i)＝Y(i)/Ymax，Y(i)是第i台机主接力器位移，毫米，Ymax是主接力器位移最大值，毫米；y0(i)是第i台机主接力器位移初值相对值，即y0(i)＝Y0(i)/Ymax，Y0(i)是第i台机主接力器位移初值，毫米，yr(i)是第i台机额定工况下主接力器位移相对值，即yr(i)＝Yr(i)/Ymax，Yr(i)是第i台机额定工况下的导叶开度，毫米；Ty(i)是第i台机组主接力器时间常数，秒，pt(i)＝Pt(i)/Pr(i)是第i台水轮机的功率相对值，Pt(i)是第i台水轮机的功率，千瓦，Pr(i)是第i台水轮机的额定功率，千瓦，At(i)是第i台水轮机的增益系数，qnl(i)＝Qnl(i)/Qr是第i台水轮机的空载流量相对值，Qnl(i)是第i台水轮机的空载流量，米3/秒，ht(i)＝Ht(i)/Hr是第i台水轮机的水头相对值，Ht(i)是第i台水轮机的水头，米，hs是调压井处的水头相对值，hs0是调压井的稳态水头相对值，qT是隧洞流量相对值。
CN201910343377.XA 2019-04-26 2019-04-26 一种通过隧洞和调压井的多机形式建立一管多机微分方程模型的方法 Active CN110222362B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910343377.XA CN110222362B (zh) 2019-04-26 2019-04-26 一种通过隧洞和调压井的多机形式建立一管多机微分方程模型的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910343377.XA CN110222362B (zh) 2019-04-26 2019-04-26 一种通过隧洞和调压井的多机形式建立一管多机微分方程模型的方法

Publications (2)

Publication Number Publication Date
CN110222362A CN110222362A (zh) 2019-09-10
CN110222362B true CN110222362B (zh) 2023-05-02

Family

ID=67820070

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910343377.XA Active CN110222362B (zh) 2019-04-26 2019-04-26 一种通过隧洞和调压井的多机形式建立一管多机微分方程模型的方法

Country Status (1)

Country Link
CN (1) CN110222362B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111143954A (zh) * 2019-12-05 2020-05-12 三峡大学 可有效降低局部水头损失的分叉输水管的制作方法
CN112651180B (zh) * 2020-12-31 2022-06-14 昆明理工大学 一种一管多机水电机组调节系统微分方程计算方法
CN115776256B (zh) * 2022-11-26 2023-07-07 昆明理工大学 一管多机水电机组励磁系统和调速系统联合控制的方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103593508A (zh) * 2013-10-23 2014-02-19 广东电网公司电力科学研究院 大型抽水蓄能电站通用仿真平台
CN106436660A (zh) * 2016-09-26 2017-02-22 昆明理工大学 一种不同位置的引水式电站发电厂房布置

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102663241A (zh) * 2012-03-29 2012-09-12 昆明理工大学 一种弹性水击下水轮机暂态功率的模拟计算方法
CN104504224B (zh) * 2015-01-16 2017-12-01 大连理工大学 一种水轮机组暂态模型的仿真方法
GB201601094D0 (en) * 2016-01-20 2016-03-02 Soliton Holdings Corp Delaware Corp Method for computational fluid dynamics and apparatuses for jet-effect use
CN105956350B (zh) * 2016-07-05 2017-04-19 华中科技大学 一种抽水蓄能机组过水系统的建模方法
CN106372296B (zh) * 2016-08-29 2019-11-08 昆明理工大学 一种具有共用管道的水轮机多机微分方程计算方法
CN106485064B (zh) * 2016-09-28 2017-09-12 华中科技大学 一种抽水蓄能机组水轮机工况智能开机方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103593508A (zh) * 2013-10-23 2014-02-19 广东电网公司电力科学研究院 大型抽水蓄能电站通用仿真平台
CN106436660A (zh) * 2016-09-26 2017-02-22 昆明理工大学 一种不同位置的引水式电站发电厂房布置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
T. Xu, L. Zhang, Y. Zeng and J. Qian.Hamiltonian Model of Hydro Turbine with Sharing Common Conduit.2012 Asia-Pacific Power and Energy Engineering Conference.2012,1-5. *

Also Published As

Publication number Publication date
CN110222362A (zh) 2019-09-10

Similar Documents

Publication Publication Date Title
CN110222362B (zh) 一种通过隧洞和调压井的多机形式建立一管多机微分方程模型的方法
CN107291990B (zh) 基于电-气互联综合能源系统暂态模型的能量流仿真方法
Xu et al. A review of dynamic models and stability analysis for a hydro-turbine governing system
CN112966394B (zh) 一种水力耦合条件水电机群动态特性的仿真方法及系统
Yu et al. Stability analysis of hydro-turbine governing system including surge tanks under interconnected operation during small load disturbance
CN112651180B (zh) 一种一管多机水电机组调节系统微分方程计算方法
CN112417662B (zh) 一种实现集中供热管网系统动态水力优化的方法
Mao et al. Investigation on optimization of self-adaptive closure law for load rejection to a reversible pump turbine based on CFD
CN102663241A (zh) 一种弹性水击下水轮机暂态功率的模拟计算方法
CN106874595A (zh) 基于节点参数化技术的输水管网计算方法
Tang et al. Dynamic regulation characteristics of pumped-storage plants with two generating units sharing common conduits and busbar for balancing variable renewable energy
Martínez-Lucas et al. Risk of penstock fatigue in pumped-storage power plants operating with variable speed in pumping mode
CN106372296B (zh) 一种具有共用管道的水轮机多机微分方程计算方法
Zeng et al. Instability analysis of pumped-storage stations under no-load conditions using a parameter-varying model
Lu et al. Influence of water diversion system topologies and operation scenarios on the damping characteristics of hydropower units under ultra-low frequency oscillations
Guo et al. Stability control of dynamic system of hydropower station with two turbine units sharing a super long headrace tunnel
CN112947076B (zh) 一种一管多机水电机组协同控制器设计方法
CN113806972B (zh) 一种考虑供需双侧灵活性的综合能源系统可靠性分析方法
CN101178051B (zh) 洞群差动式调压室
CN110941868A (zh) 一种水电站调压室临界稳定断面的计算方法
CN114329835A (zh) 一种多机共管道水轮机调节系统稳定域确定方法及系统
Lu et al. Comprehensive stability analysis of complex hydropower system under flexible operating conditions based on a fast stability domain solving method
CN100413180C (zh) 梯级串联电站联合运行控制调节器
Wang et al. Multi-objective intelligent start-up strategy for small hydropower plants based on adaptive variable parameters of water head
CN114936440A (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