下面是小编为大家整理的基于相似抛物线和插值法水泵转速确定方法研究（2022年）,供大家参考。
基于相似抛物线和插值法的水泵转速确定方法研究
我国水资源时空分布不均，与人口、生态、经济产业布局不相匹配，随着经济社会的高速发展，水资源供需矛盾日益突出，进入 21世纪以来，各地建设大型调水工程的速度明显加快。珠江三角洲水资源配置工程是国务院要求加快建设的全国 172 项节水供水重大水利工程之一，其首级取水泵站——鲤鱼洲泵站流量与扬程变化幅度极大，因此在设计中采用机组并联加变频调速的运行方式［1］。目前，我国已建的南水北调中线工程中的惠南庄泵站、牛栏江-滇池补水工程干河泵站、山西万家寨引黄工程三级提水泵站、引汉济渭工程黄金峡泵站以及正在建设的滇中引水工程石鼓水源泵站等，均存在水泵运行工况变化复杂的特点，因此泵站采用了变速调节的运行方式［2-6］。大容量变频调速技术的日益成熟，使得单机大功率且运行工况复杂的水泵变频调速成为可能，这既方便了机组的启停，又可使机组在各种复杂工况下高效运行。在泵站设计和水泵选型阶段，主要关注水泵的最佳转速［7］；而在泵站建成投入运行后，则要求水泵在满足用户需水量要求的情况下运行，此时所需流量已知，且进、出水池水位和管路阻力特性已经确定，需要确定的是水泵的实际转速。
1 水泵最佳转速的确定
对于扬程较高的泵站，当进、出水池水位变幅较大（即扬程变化较大），
且多年平均扬程与水泵的设计扬程相差较大时，或当用户的需水量变化较大时，常常采用变速调节的方法，使得水泵以最高的效率运行。水泵在不同工况下工作时，对应最高效率点的转速称为最佳转速nopt。如图 1 所示，过水泵额定转速 n0 的流量-扬程（Q-H）曲线上对应最高效率工况点 A0(QA0,HA0)作相似工况抛物线，该抛物线上的点所对应的转速即为水泵的最佳转速，相似工况抛物线方程如下式所示：
H=HA0Q2A0Q2
(1) 式中：H 为与流量相对应的扬程；HA0 为最高效率工况点所对应的扬程；QA0 为最高效率工况点所对应的流量；Q 为流量。
图 1 最佳转速确定 Fig. 1Determination of the optimum speed
根据进、出水池水位变幅得到装置静扬程 Hst，然后作装置需要扬程曲线，其与相似工况抛物线的交点为 A1、A2 等，交点所对应的转速即不同静扬程下的最佳转速，装置需要扬程方程如下式所示：
H=Hst+SQ2
(2) 式中：Hst 为装置静扬程；S 为管路系统阻力系数。
联立式（1）和（2）可得
Q=QA0HstHA0− SQ2A0− − − − − − − − − − √
(3)
H=HA0HstHA0− SQ2A0
(4) 将所得的Q和H值代入比例律公式即可用解析法求得水泵的最佳转速：
nopt=n0A0=n0HstHA0− SQ2A0− − − − − − − − − − √
(5) 式中：nopt 为最佳转速；n0 为额定转速。
水泵最佳转速的确定对泵站设计以及水泵选型具有一定的参考价值。而对于建成投入运行泵站，则需要根据用户的需水量（即流量）来进行调速。本文提出了基于相似抛物线和变斜率线性插值的方法，对调速运行中的水泵进行转速确定。当水泵调速范围不大时、水泵变速后
扬程和效率不需要修正［8，9］，可根据比例律和相似抛物线借助计算机程序来求解水泵实际运行转速（简称相似抛物线法）；当水泵调速范围较大、水泵变速后用比例律公式换算产生的误差较大时，扬程和效率均需要修正［10，11］，则借助计算机的插值程序来求解水泵转速（简称插值法）比较合适。
2 基于相似工况抛物线的水泵转速确定
当调速范围较小时，可以认为变速前后水泵的效率不变，相似工况点的流量、扬程等参数满足比例律且相似工况点位于同一条相似抛物线上。为了便于计算，往往需要根据水泵性能试验实测的数据求出水泵拟合性能曲线的数学表达式［12，13］。
2.1 基本性能曲线的数学表达式
水泵的基本性能曲线可用 1 组 n 次多项式来表示，以 Q-H 性能曲线为例，对于水泵性能试验实测的数据点(Qi,Hi)，1≤i≤N（N 为工况点最大个数），可用下面的 n 次多项式进行拟合：
H=f(Q)=a0+a1Q+a2Q2+⋯ +anQn=∑k=0nakQk
(6) 式中：i 为试验实测工况点编号；Qi、Hi 分别为第 i 个工况点水泵流量和扬程；ak（k=0,1,⋯ ,n）为多项式的系数。
为了使拟合出的近似曲线尽可能地反映所给数据的变化趋势，要求在所有数据点上的残差|δi|=|f(Qi)− Hi|为较小值。为此可以令上述偏差的平方和∑i=1N(δi)2=∑i=1N[f(Qi)− Hi]2 为最小值。这种方法称为最小二乘原则，利用这一原则确定拟合多项式 f(x)的方法即为最小二乘法多项式拟合。
确定上述多项式的过程也就是确定 f(Q)中的系数 ak（k=0,1,⋯ ,n）的过程，根据最小二乘原则，偏差平方和应为系数 ak 的函数，即
S(a0,a1,⋯ ,an)=∑i=1N(δi)2=∑i=1N[f(Qi)− Hi]2
(7) 为了使上式取值最小，则其关于 ak （k=0,1,⋯ ,n）的 1 阶导数应为 0，即
⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧∂ S∂ a0=∑i=1N2[f(Qi)− Hi]=0⇒∑i=1Nf(Qi)=∑i=1NHi∂ S∂ a1=∑i=1N2Qi[f(Qi)− Hi]=0⇒∑i=1NQif(Qi)=∑i=1NQiHi⋯ ⋯ ∂S∂ ak=∑i=1N2kQki[f(Qi)−Hi]=0⇒∑i=1NQkif(Qi)=∑i=1NQkiHi⋯ ⋯ ∂S∂ an=∑i=1N2nQni[f(Qi)− Hi]=0⇒∑i=1NQnif(Qi)=∑i=1NQniHi
（8）
将上面各等式写成方程组的形式，则有
⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧∑i=1Nf(Qi)=∑i=1NHi⇒a0N+a1∑i=1NQi+a2∑i=1NQ2i+⋯ +an∑i=1NQni=∑i=1NHi∑i=1NQif(Qi)=∑i=1NQiHi⇒a0∑i=1NQi+a1∑i=1NQ2i+a2∑i=1NQ3i+⋯ +an∑i=1NQn+1i=∑i=1NQiHi⋯ ⋯ ∑i=1NQkif(Qi)=∑i=1NQkiHi⇒a0∑i=1NQki+a1∑i=1NQk+1i+a2∑i=1NQk+2i+⋯ +an∑i=1NQn+ki=∑i=1NQkiHi⋯ ⋯ ∑i=1NQnif(Qi)=∑i=1NQniHi⇒a0∑i=1NQni+a1∑i=1NQn+1i+a2∑i=1NQn+2i+⋯ +an∑i=1NQ2ni=∑i=1NQniHi
（9）
写成矩阵形式，则有
⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧N∑i=1NQi⋮ ∑i=1NQki⋮ ∑i=1NQni∑i=1NQi∑i=1NQ2i⋮ ∑i=1NQk+1i⋮∑i=1NQn+1i⋯ ⋯ ⋮ ⋯ ⋮ ⋯ ∑i=1NQki∑i=1NQk+1i⋮ ∑i=1NQ2ki⋮ ∑i=1NQn+ki⋯ ⋯ ⋮ ⋯ ⋮ ⋯ ∑i=1NQni∑i=1NQn+1i⋮ ∑i=1NQn+ki⋮ ∑i=1NQ2ni⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⋅ ⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧a0a1⋮ ak⋮ an⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧=⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧∑i=1NHi∑i=1NQiHi⋮ ∑i=1NQkiHi⋮ ∑i=1NQniHi⎧⎧⎧⎧
⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧
(10) 解此方程组，即可求得 Q-H 曲线数学表达式的系数 ak(k=0,1,⋯ ,n)。
为了提高计算速度和精度，还应考虑相关系数 R 和标准残差 s，其表达式分别为
R=∑i=1N(Qi− Q¯ ¯ ¯ )(Hi− H¯ ¯ ¯ )∑i=1N(Qi− Q¯ ¯ ¯ )2(Hi−H¯ ¯ ¯ )2− − − − − − − − − − − − − − − − − ⎧
(11)
s=1N∑i=1N[∑k=0nakQki− Hi]2− − − − − − − − − − − − − − −− − ⎧
(12) 式中：Q¯ ¯ ¯ =1N∑i=1NQi，H¯ ¯ ¯ =1N∑i=1NHi。
相关系数表示 1 组数据的离散程度，通过计算 RR 值，可确定拟合多项式的方次，以达到减少不必要的计算，提高计算速度的目的。RR值越小，说明该组数据的离散程度越大。标准残差 ss 表示曲线方程对各点拟合的程度，ss 值越小，说明拟合程度越高。
在 实 际 计 算 过 程 中 ， 为 了 防 止 运 算 溢 出 ， QiQi 常 用
Qi− Q¯ ¯ ¯ Qi-Q¯ 代替，则拟合多项式为
H=a0+a1(Q− Q¯ ¯ ¯ )+a2(Q− Q¯ ¯ ¯ )2+⋯ +an(Q− Q¯ ¯ ¯ )n=∑k=0nak(Q− Q¯ ¯ ¯ )k
(13) 利用上述方法即可求出适当精度的 Q-H、Q-P（轴功率）、Q-NPSH（空化余量）和 Q-η（水泵效率）曲线的数学表达式。
2.2 满足需水量要求的水泵实际运行转速
水泵的性能曲线与水泵的转速有关。改变水泵的转速，可使水泵的性能曲线发生改变，从而达到调节水泵工作点的目的。
变速调节只改变了水泵的性能曲线，装置需要扬程曲线并未改变，如图 2 所示。图中曲线 0、1 分别为水泵转速为 n0n0 和 n1n1 时的 QQ-HH曲线，方程 Hr=Hst+SQ2Hr=Hst+SQ2 对应的曲线为装置需要扬程（HrHr）曲线。由图 2 可知，当水泵在转速 n0n0 下运行时，其工作点为 A0A0，对应的流量为 QA0QA0；转速减小为 n1n1 后，工作点变为 A1A1，流量、功率和效率也分别从原来的QA0QA0、PA0PA0和ηA0ηA0变为QA1QA1、PA1PA1 和 ηA1ηA1。
图 2 变速调节 Fig. 2Speed regulation
通常遇到的情况如图 3 所示，已知水泵在额定转速为 n0n0 的 Q-H 曲线和装置需要扬程曲线 Hr=Hst+SQ2Hr=Hst+SQ2，此时水泵的工作点为 A0(QA0,HA0)A0(QA0,HA0)，对应的流量为 QA0QA0，但实际用户所需流量为 QA1QA1，由进、出水池水位和管路阻力特性计算出需要扬程为 HA1HA1，对应的水泵工作点 A1(QA1,HA1)A1(QA1,HA1)不在转速为 n0n0 时的 Q-H 曲线上。为了使水泵能在新的工作点 A1A1 工作，需要求出转速 n1n1。
图 3 变速调节的转速确定 Fig. 3Speed determination of variable speed regulation
假定变速后的转速 n1n1 和水泵额定转速 n0n0 相差不大，变速前后各效率可以认为是相等的，由比例律可知，水泵的流量、扬程与转速的关系为
QA2QA1=n0n1
(14)
HA2HA1=n20n21
(15) 只有相似工况抛物线上的工况点才能应用比例律来进行性能的换算，为此，将所需流量 QA1QA1 值与扬程 HA1HA1 值代入相似工况抛物线方程 H=KQ2=Q2⋅ HA1/Q2A1H=KQ2=Q2⋅ HA1/QA12 ， 求 出K=HA1/Q2A1K=HA1/QA12 的值，相似工况抛物线通过点 A1A1 并与原有的 QQ-HH 曲线相交于 A2A2 点，由此求出 A2A2 点的坐标（QA2QA2，HA2HA2），如图 3 所示，用下式可求出 n1n1：
n1=n0QA1QA2
(16) 由以上可得，采用相似抛物线法数值求解所需流量 QA1QA1 所对应水泵转速 n1n1 的计算步骤如下：
1 ）
将 用 户 实 际 需 流 量 QA1QA1 代 入 装 置 需 要 扬 程 方 程Hr=Hst+SQ2Hr=Hst+SQ2，计算得到变速后工作点 A1A1 的扬程 HA1HA1。
2）计算过点A1A1的相似工况抛物线的KK值，K=HA1/Q2A1K=HA1/QA12。
3）联立求解相似工况抛物线方程和水泵在转速为 n0n0 时的 QQ-HH 曲线方程，如下式所示：
⎧⎧⎧⎧⎧⎧⎧H=KQ2H=a0+a1(Q− Q¯ ¯ ¯ )+a2(Q− Q¯ ¯ ¯ )2+⋯ +a
n(Q− Q¯ ¯ ¯ )n
(17) 求解得到的交点坐标（QA2QA2，HA2HA2）即为工况点 A2A2 的流量QA2QA2 和扬程 HA2HA2。
4）采用比例律公式，计算水泵对应于工况点 A1A1 时的转速 n1n1：
n1=n0QA1QA2 或 n1=n0HA1HA2− − − − √
(18)
3 基于插值法的水泵转速确定
当水泵调速范围较大时，变速前后水泵的流动效率不相等，虽然相似工况点的流量、扬程等参数仍满足相似律，相似工况点仍位于同一条相似曲线上，但该曲线不再为抛物线，无法如前所述推导出相似抛物线方程，基于相似工况抛物线的方法无法实现对水泵转速的确定，针对这种情况，提出了插值法。
3.1 变速通用性能曲线
基本性能曲线反映了水泵在额定转速下的性能。然而，在实际运行中
常常采用改变转速的方法来适应工况变化对水泵性能的要求。为此，需要知道水泵在各种不同转速时的性能。同一直角坐标中绘制的一台水泵在不同转速时的一簇性能曲线称为变速通用性能曲线。通过对水泵的性能试验，可以得到某一转速下水泵的基本性能曲线，如果在试验中改变水泵转速并重复试验，可以得到不同转速下水泵相应的基本性能曲线。但该曲线使用不方便，且转速变化档次越多，线条越密集。为便于使用，通常将性能曲线变换为仅在图中反映不同转速下的QQ-HH 曲线，且用等效率线代替不同转速时的效率曲线，如图 4 所示。
图 4 不同转速下的通用性能曲线 Fig. 4General performance curves at different rotating speeds
水泵在不同转速下的基本性能曲线还可以通过理论公式换算得到。按照叶片泵的相似理论，同一台水泵的转速变化时，其变化前后的性能是相似的。因此，可以利用相似律换算变速后的性能，即可根据已知转速的性能曲线换算出另一转速的性能曲线。
但相似律只适用于 2 个相似的工况，即只有同一条相似工况曲线上的点才能应用相似律公式进行相似性能换算。水泵变转速下的流量、扬程和流动效率按下式换算：
⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧Q0Q1=n0n1⋅ ηq0ηq1H0H1=n20n2
1⋅ ηh0ηh11− ηh1− ηh0=(n0n1)β
(19) 式中：ηqηq、ηhηh 分别为容积效率、水力效率；下标 0 代表额定转速 n0n0 下的相应参数；下标 1 代表其他任意转速 n1n1 下的相应参数；ββ 由试验确定，不同的文献所给出的 ββ 值有所不同，ββ=0.16~0.5 ［14-19］，本文按 IEC 60193 的规定取 ββ=0.16 ［9］；一般认为，水泵改变转速时其机械效率 ηmηm 和容积效率 ηqηq 不变［9］，即有 ηq0=ηq1ηq0=ηq1。
可以看出，通过实测得到的不同转速下通用性能曲线来发现不同转速下相似工况点之间的关系是困难的，并且即便使用相似换算，但变速前后效率不相等时，无法如前所述导出相似工况抛物线方程进而求出水泵变速运行时的转速，此时，采用插值法确定水泵变速调节运行时的转速更为合适。
3.2 插值法确定水泵变速调节运行时的转速
泵站建成投入运行后，水泵机组一般按用户需水量要求来运行。假定根据用户需要的流量确定水泵流量为 QAQA，由于水泵的工作点须落在装置需要扬程曲线上，根据装置需要扬程方程确定水泵扬程为HAHA。
如图 5 所示，装置需要扬程 ...
