水力发电机组甩负荷过渡过程的参数敏感性分析

谢卓健; 闫懂林; 肖汉; 麦先春; 王峻峰; 陈启卷; Zhuo-jian XIE; Dong-lin YAN; Han XIAO; Xian-chun MAI; Jun-feng WANG; Qi-juan CHEN

中国农村水利水电 ›› 2021 ›› (8) : 195-199,206.

水电建设

水力发电机组甩负荷过渡过程的参数敏感性分析

作者信息 +

Parameter Sensitivity Analysis of Transition Process for Hydraulic Generating Set

Author information +

稿件信息 +

摘要

以国内某水电站为对象，在不确定框架下构建了水力发电机组甩负荷过渡过程的数学模型。并以此模型为基础，引入扩展傅里叶幅度检验法量化输入参数对系统响应的影响，获得了机组在甩负荷过渡过程中输入参数相对于转速、水击压力以及调压室涌浪的敏感性特性。可以得到飞轮力矩对转速最大上升值的影响最为显著。同时，也可以发现各参数的不确定度会明显影响对水击压力的敏感性，而不会影响对调压室涌浪的敏感性。

Abstract

Taking a hydropower station in China as an example， a mathematical model of load rejection process of hydroelectric generating units is constructed under the uncertain framework. Based on this model， an Extended Fourier Amplitude Sensitivity Test is introduced to quantify influences of input parameters on the system response， and sensitivity characteristics of input parameters relative to the rotational speed， water hammer pressure and surge in the surge-chamber during the load rejection transition are obtained. It can be concluded that the flywheel torque has the most significant effect on the maximum speed appreciation. At the same time， it can be found that the uncertainty of each parameter will obviously affect the sensitivity to the water hammer pressure， but will not affect the sensitivity to the surge in the surge chamber.

关键词

水力发电机组 / 甩负荷 / 全局敏感性 / 扩展傅里叶幅度检验法

Key words

hydroelectric generating set / load shedding / global sensitivity / Extended Fourier Amplitude Sensitivity Test

基金

国家自然科学基金项目(52009096)

大唐集团公司科研项目“大型混流式水电机组有功振荡控制优化关键技术及实践”

引用本文

EndNote

Ris (Procite)

Bibtex

导出引用

谢卓健 , 闫懂林 , 肖汉 , 麦先春 , 王峻峰 , 陈启卷. 水力发电机组甩负荷过渡过程的参数敏感性分析[J].中国农村水利水电, 2021(8): 195-199,206

Zhuo-jian XIE , Dong-lin YAN , Han XIAO , Xian-chun MAI , Jun-feng WANG , Qi-juan CHEN. Parameter Sensitivity Analysis of Transition Process for Hydraulic Generating Set[J].China Rural Water and Hydropower, 2021(8): 195-199,206

0 引言

在过去的二十年里，我国水电的开发和利用得到了快速的发展，在西南地区，各种大流量、高水头、包含复杂水力系统的水电站不断兴起。到2019年底，我国的水电总装机容量已经达到3.56 亿kW^［1］。对于水电站来说，由于电网故障等原因，运行中随时可能会经历机组甩负荷过程，而一些水电站的重大事故，例如飞逸、压力钢管爆裂等，往往就是由于这种较为极端的工况引起^［2］。如何有效地保证甩负荷过渡过程中机组的安全可靠，是实际水电工程中所面临的重要问题。

到目前为止，针对水电机组甩负荷过渡过程进行了大量研究。部分学者对实际电站发生的甩负荷事故进行了分析和讨论。黄菊珠^［3］分析了山美水电站因110 kV线路开关跳闸，造成3台机组满负荷高频切机，引发机组过速停机的事故，并提出了相关的预防措施。曲东^［4］分析了平班水电站在甩负荷试验时调速器调节失灵，导叶不能回关导致机组过速停机的事故，提出了相应的处理方法，并对调速器重要参数进行了优化整定。张伟^［5］针对西溪电站甩负荷试验时出现的机组过速停机及导叶快关时间不稳定等问题进行了分析，并提出了故障处理方法。除此之外，也有大量学者对水力发电机组甩负荷过渡过程进行了数值模拟研究。付亮^［6］针对有调压井长引水隧洞水电站，利用数值模拟的方法研究了调压井水位波动与流道压力的变化过程。王威^［7］等以伊朗卡仑水电站为例，分析了具有上游调压室、一管两机布置的水电站机组甩负荷过渡过程，并讨论了导叶关闭规律对转速上升和水击压力的影响，提出了合理的导叶关闭规律。叶文波等^［8］针对南方某大型水电站，进行了机组甩负荷大波动过渡过程的计算，分析了引水系统最高压力升高和机组转速的最大升高特性。陈胜等^［9］研究了水电站中机组相继甩负荷情况下调压室涌浪的叠加特性。张洋^［10］研究了设置气垫式调压室的超长引水隧洞水电站的甩负荷过渡过程，提出了气垫式调压室对系统水击压力的影响效应。在这些分析和讨论中，都是在确定性框架下研究水电站机组的甩负荷过渡过程，而并没有考虑到参数不确定性对甩负荷过渡过程的影响。实际上，水电站过渡过程分析往往参数众多且包含复杂的非线性动态过程，参数尺度多样化及制造误差等多种不确定性因素均会导致计算结果的不确定性。所以，在不确定框架下研究参数对甩负荷过渡过程的影响效应就显得更有必要。

参数敏感性分析是量化系统参数在设计或设定的可能范围内变化时对系统输出响应影响的一种方法^［11］。通过敏感性分析，可以获得对系统输出影响较为显著的输入参数，无论在设计、运行还是维护过程中，只需要重点关注这些参数而忽略一些不敏感的参数，就可以大大简化工作量。通常的参数敏感性分析包括局部敏感性分析和全局敏感性分析，局部敏感性分析只能反映单一参数变化对系统响应的影响，而全局敏感性分析则可以有效地反映多个参数同时变化对系统响应的影响^［12］。在全局敏感性分析方法中，主要包括：回归分析法^［13］、筛选法^［14］，Sobol指数法^［15］和扩展傅里叶幅度检验法（extended Fourier amplitude sensitivity test，EFAST）^［16］等。其中EFAST方法可以获得输入参数之间的交互效应，已经在很多领域得到了广泛应用^{［17，18］}。因此，本研究中将利用EFAST全局敏感性分析方法来研究水力发电机组在甩负荷过渡过程中输入参数对转速、水击压力和调压室涌浪的敏感性特性。

相比于之前的研究，本文的主要创新在于考虑参数的不确定性，在不确定框架下建立了水力发电机组的甩负荷过渡过程的模型，并以此模型为基础进行了数值研究，量化分析了输入参数对机组甩负荷过程中机组转速、水击压力和调压室涌浪的敏感性特性。相关的理论方法和研究结论可以应用于水电站的设计及维护工作。

1 水力发电机组甩负荷过渡过程建模

水电站的引水系统在空间分布上主要包括上游水库，引水隧洞，上游引水调压室，压力钢管，蜗壳，尾水管及尾水明渠段。本节将首先利用特征线方法建立该电站引水系统的数学模型，同时结合水轮机特性曲线的插值模型，调速器数学模型，发电机及负载数学模型，建立完整的水电机组甩负荷过渡过程的数学模型。

1.1 引水系统模型

1.1.1 有压过流管道特征线方程

利用特征线法将有压输水系统非恒定流的运动方程和连续性方程转化为在特征线dx/dt = V±a上的常微分方程组［19］，可以得到：

C^{+} : \{\begin{array}{l} \frac{d H}{d t} + \frac{a}{g A} \frac{d Q}{d t} + \frac{a f}{2 g D A^{2}} Q | Q | = 0 \\ \frac{d x}{d t} = a \end{array}

（1）

C^{-} : \{\begin{array}{l} \frac{d H}{d t} - \frac{a}{g A} \frac{d Q}{d t} - \frac{a f}{2 g D A^{2}} Q | Q | = 0 \\ \frac{d x}{d t} = - a \end{array}

（2）

式中：H表示测压管水头；a表示水击波速；f为达西-维斯巴哈摩擦系数；D表示管道直径；A为流道截面面积；Q为流道中的流量。在这里，V和H是与时间和空间相关的函数，其中t为时间，x表示管道长度。

在利用特征线法离散化求解微分方程时，需要对有压管道按照等距离进行划分，得到长度为Δx的N段管，并用序号i进行标记，定义起始端为i=1，则末端为i=N+1。设定仿真步长为ΔT=Δx/a。利用特征线法进行计算机数值求解的原理如图1所示，横坐标表示管道的空间分布，纵坐标表示采样时刻分布。从中可以看出，根据上一时刻的节点i-1和节点i+1的值，可以计算出下一时刻节点i的水头、流量值。

图1 特征线法原理示意

Fig.1 Schematic diagram of characteristics method

Full size|PPT slide

对公式（1）和（2）的特征线方程积分并离散化，可以得到t+ΔT时刻节点i处的水头、流量关系方程式为：

\{\begin{array}{l} H_{i, t + Δ T} = C_{P i, t + Δ T} - B_{P i, t + Δ T} \cdot Q_{i, t + Δ T} \\ H_{i, t + Δ T} = C_{M i, t + Δ T} + B_{M i, t + Δ T} \cdot Q_{i, t + Δ T} \end{array}

（3）

式中：B = a/Ag； R = fΔx/（2gDA ²）； C_Pi _， _t _+Δ _T = H_i _-1， _t +BQ_i _-1， _t；B_Pi _， _t _+Δ _T = B+R|Q_i _-1， _t|；C_Mi _， _t _+Δ _T = H_i _+1， _t -BQ_i _+1， _t；B_Mi _， _t _+Δ _T = B+R|Q_i _+1， _t|。

1.1.2 水轮机边界

水轮机的瞬态工作水头定义为水轮机蜗壳进口和尾水管出口的水头差，如图2所示，当水轮机上、下游节点为1和2时，则水轮机工作水头为H_T = H ₁-H ₂，瞬态的过机流量为Q_T =Q ₁=Q ₂。

图2 水轮机简化结构图

Fig.2 Structure chart of hydro-turbine

Full size|PPT slide

根据水轮机上下游进出口节点的水头平衡方程，可以得到：

\begin{array}{l} Q_{T, t + Δ T} = - C_{V} B_{T, t + Δ T} \\ + \sqrt[]{{(C_{V} B_{T, t + Δ T})}^{2} + 2 C_{V} (C_{P 1, t + Δ T} - C_{M 2, t + Δ T})} \end{array}

（4）

式中：

C_{V} = (Q_{11} D_{1}^{2})^{2} / 2

；

B_{T, t + Δ T} = B_{P 1, t + Δ T} + B_{M 2, t + Δ T}

；D ₁为转轮直径；Q ₁₁为单位流量，可以通过特性曲线得到。

水库边界，上游调压室边界，蜗壳及尾水管当量管的计算可以参考文献［19］。

1.2 水轮机特性曲线插值模型

依据文献［20］的方法，对原始的水轮机综合特性曲线进行扩展，可以得到Q ₁₁=f ₁（n ₁₁，y）和M ₁₁=f ₂（n ₁₁，y）的关系，用于甩负荷过渡过程的计算。在这两个关系式中，n ₁₁表示单位转速，M ₁₁表示单位力矩，y表示导叶开度。以当前研究的水电站为例，得到扩展后的Q ₁₁-n ₁₁-y和M ₁₁-n ₁₁-y关系曲线如图3所示。

图3 单位流量和单位力矩图

Fig.3 Diagrams of unit discharge and unit moment

Full size|PPT slide

1.3 发电机及负载数学模型

水力发电机组的力矩平衡方程^［19］可以表示为：

J \frac{d ω}{d t} = M - M_{g}

（5）

式中：J为机组转动部分的转动惯量，J=GD ²/4，t·m²；GD ²为飞轮力矩，t·m²；ω表示机组旋转角速度，rad/s；M和M_g 分别是水轮机的主动力矩和阻力矩，kN·m。

假定水轮机力矩在ΔT时刻里线性变化，则可以利用梯形积分公式对上式积分。同时，将旋转角速度（ω）转换为转速（n），可以得到：

n_{t + Δ T} = n_{t} + \frac{120}{π G D^{2}} (\frac{M_{t} + M_{t + Δ T}}{2} - M_{g}) Δ T

（6）

至此，水力发电机组甩负荷过渡过程的确定性模型被建立，在这个确定性模型中将确定性的参数用不确定的分布或者范围进行替代，即可以获得不确定性的模型，以此为基础，相关的参数敏感性可以利用方便的非侵入式EFAST方法进行讨论。

2 全局敏感性分析

2.1 全局敏感性分析方法

本文采用EFAST方法来对水电机组甩负荷过渡过程的参数敏感性进行分析。EFAST法的原理是在输入参数的多维空间中挑选合适的搜索曲线，并给所有输入参数分配非线性相关的整数频率值，从而引入搜索函数，使模型作为独立参数的周期函数，则可以把多维积分问题转换为一维积分问题。然后通过把目标函数变为傅里叶级数，计算出各频率的傅里叶频谱曲线，由频谱曲线获得由输入参数x_i 引起的系统输入方差和总方差，其比值则被定义为该参数的敏感度指标。文献［16-18］中给出了具体的算法推演流程，这里仅给出敏感度计算的公式。

总方差可以写成各个子方差和的形式：

V = \sum_{i} V_{i} + \sum_{i \neq j} V_{i j} + \dots + \sum_{i \neq j \neq \dots \neq k} V_{i j \dots k}

（7）

式中：V为模型总方差；V_i 表示单一输入变化引起的模型方差；V_ij 表示通过x_i 和x_j 交互作用引起的方差；V_ij _… _k 表示x_i，x_j，…，x_k 交互作用引起的耦合方差。

主敏感性（一阶敏感性）指标表示单个参数独立作用时对系统输出响应的贡献，可以表示为：

S_{i} = \frac{V_{i}}{V}

（8）

总敏感性指标为：

S_{T i} = \frac{V - V_{- i}}{V}

（9）

式中：V _- _i 表示不包括参数x_i 的所有参数方差和。

2.2 仿真设计与参数选择

以国内某水电站为例，进行机组甩100%负荷时过渡过程计算的参数敏感性分析，涉及的主要参数如表1所示。敏感性分析的输入参数均值选取设计值，方差分别取变异系数（=标准差/均值）为0.01和0.05两种情况式的值，代表不同的不确定度。选择Monte Carlo方法进行参数的随机采样，采样样本为5 000，满足EFAST法对样本个数的要求。另外，由于其他水工建筑物的影响，该电站不同高程下的调压室面积会有所有差别，在高程H_s ₀<1 458.7 m时，面积A_s ₂为801.975 m²；当1 458.7 m<H_s ₀< 1 476.5时，面积A_s ₃为2 280.975 m²；当H_s ₀>1 476.5 m时，面积A_s ₄为2 562 m²。

表1 仿真参数

Tab.1 Simulation parameters

参数	符号	值
机组飞轮力矩/（t·m²）	GD ²	65 000
调压室阻抗孔面积/m²	A_s ₁	58.205
调压室第一段面积/m²	A_s ₂	801.975
调压室第二段面积/m²	A_s ₃	2 280.975
调压室第三段面积/m²	A_s ₄	2 562
调压室流入系数	r ₁	3.564e-5
调压室流出系数	r ₂	3.564e-5
引水隧洞波速/（m·s^-1）	a ₁	1 100
引水隧洞摩擦系数	f ₁	0.03
引水隧洞面积/m²	A ₁	162.163 9
调压井至钢管下平段首波速/（m·s^-1）	a ₂	1 166.667
调压井至钢管下平段首摩擦系数	f ₂	0.02
调压井至钢管下平段首面积/m²	A ₂	74.365 2
钢管下平段首至机组段波速/（m·s^-1）	a ₃	1 065
钢管下平段首至机组段摩擦系数	f ₃	0.005
钢管下平段首至机组段面积/m²	A ₃	52.738 9
尾水管段波速/（m·s^-1）	a ₄	1 025
尾水管段摩擦系数	f ₄	0.01
尾水管段面积/m²	A ₄	62.841 1
尾水管延伸段波速/（m·s^-1）	a ₅	1 107.33
尾水管延伸段摩擦系数	f ₅	0.044
尾水管延伸段面积/m²	A ₅	141.761 8

3 分析与讨论

在机组的甩负荷过程中，最高转速、最大水击压力、调压室最大涌浪和最小涌浪需要被重点关注，也是评估机组甩负荷过程安全性的重要指标，所以，在本节中，参数相对于这几个输出变量的敏感性将被重点分析。

3.1 转速升高

机组甩负荷过程中最高转速的全局敏感性分析结果如图4所示。

图4 最高转速的全局敏感性

Fig.4 Global sensitivity of maximum speed

Full size|PPT slide

以甩负荷过程中的机组最高转速作为输出时，无论在变异系数为0.01还是0.05的情况下，机组旋转部分的飞轮力矩都是对其最为敏感的参数，其他参数的影响微乎其微。另外，对比飞轮力矩的主效应和总效应指标，可以发现在不同的不确定度下其差值都很小，即飞轮力矩对甩负荷过程中转速升高的影响也主要是通过自身的作用实现，与其他参数的交互作用对转速升高的影响可以忽略。最后，对比在变异系数为0.01和0.05情况下对最高转速的参数敏感性特性，发现变异系数为0.05时上游引水隧洞面积对最高转速的敏感性要高于变异系数为0.01时的敏感性。其主要原因是，在系统变异系数从0.01同尺度放大到0.05时，引水隧洞面积对最高转速影响的放大尺度要高于其他参数对最高转速影响的放大尺度，从而导致不同不确定度下会出现这种差异。当然从整体的角度来说，飞轮力矩的绝对主导地位仍然不会受到威胁。总体来说，飞轮力矩对于甩负荷过程中的机组转速升高的影响最为显著。

3.2 水击压力

除了转速的升高，水击压力也是甩负荷过程中影响机组安全的关键指标。如图5所示，给出了机组甩负荷过程中最高水击压力的全局敏感性分析结果。

图5 最高水击压力的全局敏感性

Fig.5 Global sensitivity of maximum water hammer pressure

Full size|PPT slide

当变异系数为0.01时，输入参数对水击压力敏感性的主效应指标排序为：A ₂>a ₁>A_s ₄>r ₁>A_s ₃>a ₂>a ₄>a ₃>A ₅>f ₁>f ₂>A_s ₁>A ₄>A ₃>GD ²>f ₅>f ₃>A ₁>a ₅>A_s ₂>f ₄>r ₂；输入参数对水击压力敏感性的总效应指标排序为：A ₂>a ₃>A_s ₄>A_s ₃>A ₅>a ₁>f ₁>a ₂>f ₅>r ₁>a ₄>f ₂>a ₅>GD ² r ₂>A_s ₁>A ₃>A ₁>A_s ₂>f ₄>A ₄>f ₃。无论是主效应还是总效应，压力钢管面积对水击压力的敏感性最为显著，且从中可以看出A ₂独立的影响占46.6%，与其他参数的交互效应占53.4%，也就是说A ₂对水击压力的独立作用和与其他参数的交互作用效果基本持平。此外，蜗壳水击波速a ₃，调压室面积A_s ₄和A_s ₃的总效应指标都超过了0.4，它们对最高水击压力的影响也仅次于A ₂。与A ₂的作用方式不同的是，a ₃、A_s ₄和A_s ₃都是主要通过与其他参数的交互作用来实现对水击压力的影响。

当变异系数为0.05时，输入参数对水击压力敏感性的主效应指标排序为：a ₃ > A ₂ > r ₁ > A ₄ > A_s ₃ > A ₁ > A ₃ > a ₄ > f ₁ > GD ² > A ₅ > A_s ₁ > A_s ₄ > A_s ₂> a ₅ > a ₂ > a ₁ > f ₂ > f ₃ > f ₅ > f ₄ > r ₂；输入参数对水击压力敏感性的总效应指标排序为：a ₃ > A ₂ > A_s ₃ > r ₁ > A ₄ > f ₅ > A ₅ > f ₁ > A_s ₄ > GD ² > f ₂ > a ₁ > f ₄ > A_s ₂ > a ₅ > a ₄ > a ₂ > A ₃ > r ₂ > A_s ₁ > A ₁ > f ₃。在这种情况下，蜗壳水击波速a ₃会超越压力钢管面积A ₂成为最敏感的参数。同时，a ₃独立作用对水击压力影响的占比也从变异系数为0.01时的11.51%增加到49.54%，这也就意味着随着系统不确定度的增加，a ₃独立的影响会被放大。其主要原因是系统不确定度的增加会使得参数的变化范围也相应地扩大，在这个过程中，a ₃在偏离标准值越远时其单参数项对水击压力影响的变化程度更加显著，从而使得其独立影响的占比在不确定度增加后提升而交互效应的占比下降。此外，相比于变异系数为0.01时的敏感性指标，A_s ₄和A_s ₃的影响在变异系数为0.05时会被削弱。值得注意的是，在变异系数为0.05时，调压室流入系数r ₁和蜗壳截面积A ₄对最大水击压力的影响会显著提升。

3.3 调压室涌浪

图6描述了机组甩负荷过程中调压室最高涌浪的全局敏感性分析结果。在变异系数为0.01时，调压室截面积（A_s ₂、A_s ₃和A_s ₄），调压室流入系数（r ₁）以及引水隧洞的面积（A ₁）对最高涌浪较为敏感，它们之间的总敏感性排序为：A_s ₃ > A ₁ > A_s ₄ > r ₁ > A_s ₂ > 其他。同时，可以注意到参数A_s ₃，A ₁，A_s ₄和r ₁主要是通过它们独立的作用直接影响最高涌浪高度，而参数A_s ₂则是通过与其他参数的交互过程形成对调压室最高涌浪高度的影响。剩余的参数对调压室最高涌浪的影响则不明显。在变异系数为0.05时，输入参数对调压室最高涌浪的总敏感性排序仍然为：A_s ₃ > A ₁ > A_s ₄ > A_s ₂ > r ₁ > 其他。另外，与变异系数为0.01时的参数敏感性指标进行对比，可以发现它们之间的差异并不明显。也就是说，参数不确定性的变化对调压室最高涌浪的影响比较微弱，可以忽略不计。

图6 调压室最高涌浪的全局敏感性

Fig.6 Global sensitivity of highest surge

Full size|PPT slide

图7给出了机组甩负荷过程中调压室最低涌浪的全局敏感性分析结果。无论变异系数为0.01还是0.05，调压室面积A_s ₃和压力隧洞面积A ₁是对调压室最低涌浪影响最为显著的两个参数，而其他参数的影响则相对较小。与调压室最高涌浪的敏感性结果对比可以发现，参数A ₁和A_s ₄的影响会明显减小，但A_s ₃对最高涌浪和最低涌浪的影响无论在哪个不确定度下都十分显著。

图7 调压室最低涌浪的全局敏感性

Fig.7 Global sensitivity of minimum surge

Full size|PPT slide

4 结论

本文以国内某水电站为例，建立了水力发电机组甩负荷过渡过程的数学模型，并考虑到参数的不确定性，引入EFAST方法对输入参数在不同不确定度下进行了敏感性分析，可以得到以下结论。

（1）对于机组甩负荷过程中转速的最大上升值，机组的飞轮力矩对其影响占据了绝对的主导地位，并且这种影响效应是飞轮力矩独立作用的效果而不是与其他参数交互影响的结果，且参数的不确定度对参数敏感性的影响不明显。

（2）在较小的不确定度下，压力钢管面积A ₂，蜗壳波速a ₃，调压室面积A_s ₄和A_s ₃，以及尾水管面积A ₅是对水击压力影响最为显著的5个参数。而当参数的不确定度增加后，蜗壳波速a ₃对水击压力的影响会提升并成为影响最显著的参数，而调压室面积A_s ₄和A_s ₃的影响在不确定度增大后会被削弱。

（3）对于调压室涌浪，调压室截面积A_s ₃及引水隧洞的面积A ₁对调压室最高最低涌浪都较为敏感，而调压室面积A_s ₂和A_s ₄以及调压室流入系数r ₁仅对最高涌浪的影响较为明显。同时，参数的不确定度并不会影响对调压室涌浪的敏感性特性。

考虑到本文的主题是讨论参数不确定框架下系统参数对输出响应的敏感性特性，以及考虑到当前关于参数集引入数学模型展开敏感性分析的实际困难，这里没有对特性曲线引起的模型不确定性进行分析，仅针对参数的不确定性进行了研究，在后续的工作中可以继续考虑改进分析方法来展开模型不确定性的相关研究。

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序

1	周建平，杜效鹄，周兴波.全球水电开发现状及未来趋势［J］.中国电业，2020（7）：26-29. 本文引用 [1]

2	张良君，李小阳.水电站压力钢管水压试验安全质量事故分析［J］.中国科技信息，2014（8）：251-253. 本文引用 [1]

3	黄菊珠.山美水电站机组甩负荷事件分析［J］.广西水利水电，2020（4）：87-89. 本文引用 [1]

4	曲东.平班水电站甩负荷试验异常分析及处理［J］.红水河，2015，34（3）：46-49. 本文引用 [1]

5	张伟.西溪电站机组甩负荷试验故障分析［J］.水电站机电技术，2006（2）：44-45. 本文引用 [1]

6	付亮，王义国.有调压井水电站甩负荷试验与仿真分析［J］.水电能源科学，2012，30（6）：154-157. 本文引用 [1]

7	王威，王知行，钟师胜.伊朗卡仑I水电站机组甩负荷过渡过程的数值模拟计算［J］.大电机技术，2005（3）：44-48，52. 本文引用 [1]

8	叶文波，王煜，姜锋，等.水电站甩负荷大波动过渡过程模拟［J］.水电自动化与大坝监测，2015，39（1）：4-6，37. 本文引用 [1]

9	陈胜，张健，俞晓东.水电站输水系统相继甩负荷下调压室涌浪叠加研究［J］.水利学报，2015，46（11）：1 321-1 328. 本文引用 [1]

10	张洋，杨建东，郭文成.设气垫式调压室的超长引水隧洞水电站大波动过渡过程探讨［J］.大电机技术，2016（4）：42-48，64. 本文引用 [1]

11	何亮，侯英雨，赵刚，等.基于全局敏感性分析和贝叶斯方法的WOFOST作物模型参数优化［J］.农业工程学报，2016，32（2）：169-179. 本文引用 [1]

12	姜书鹏，乔颖，徐飞，等.风储联合发电系统容量优化配置模型及敏感性分析［J］.电力系统自动化，2013，37（20）：16-21. 本文引用 [1]

13	IOOSS B， LEMATRE P. A review on global sensitivity analysis methods［J］. Operations Research/ Computer ence Interfaces Series， 2014，59. 本文引用 [1]

14	WANG C， PENG M， XIA G， Sensitivity analysis based on Morris method of passive system performance under ocean conditions［J］， Ann. Nucl. Energy， 2020，137：107 067. 本文引用 [1]

15	SALTELLI A， ANNONI P， AZZINI I， et al. Variance based sensitivity analysis of model output. Design and estimator for the total sensitivity index［J］. Computer Physics Communications， 2010，181（2）：259-270. 本文引用 [1]

16	JIN X， LI Z， NIE C， alet， Parameter sensitivity analysis of the AquaCrop model based on extended fourier amplitude sensitivity under different agro-meteorological conditions and application ［J］， F. Crop. Res.， 2018，226：1-15. 本文引用 [2]

17	邢会敏，相诗尧，徐新刚，等.基于EFAST方法的AquaCrop作物模型参数全局敏感性分析［J］.中国农业科学，2017，50（1）：64-76. 本文引用 [1]

18	李艳，黄春林，卢玲.基于EFAST方法的SEBS模型参数全局敏感性分析［J］.遥感技术与应用，2014，29（5）：719-726. 本文引用 [2]

19	陈家远.水力过渡过程的数学模拟及控制［M］.成都：四川大学出版社，2008. 本文引用 [2]

20	叶文波. XW水电站水力机组大波动过渡过程研究［D］.2015. 本文引用 [1]

PDF(2662 KB)

2513

访问

引用

详细情况

段落导航

收稿日期	出版日期
2020-10-27	2021-08-15
发布日期
2021-09-24

选择文件类型/文献管理软件名称

选择包含的内容