本发明涉及管理用于泵送源自可能含有杂质的自然环境的水的设施，所述设施旨在利用至少一个旋转式过滤器来净化所泵送的水，而所述杂质中的至少一些杂质至少部分地堵塞所述旋转式过滤器。具体地，规定了至少基于以下来估计由堵塞所述旋转式过滤器的杂质引起的压力损失随时间推移的演变：与所述自然环境相关的数据、所述旋转式过滤器的尺寸，以及与所述过滤器上游的至少一个水位Namont和所述过滤器下游吸入的水的流速Qaspiré有关的局部测量结果。对所述压力损失随时间推移的所述演变的这种估计使得能够预期所述过滤器(FIL)下游的泵送系统(POM)供水不足的可能风险，或甚至使得能够设计出适应于特定泵送设施的需求的过滤器(FIL)。
对旋转式过滤器涉及堵塞风险的优化
[0001] 本发明涉及从自然场地或室外场地过滤水的领域。例如，场地可以是河流、湖泊、大海或海洋。因此，水中可能会掺入杂质(藻类、海洋动物、树叶或枯枝和其它废物，尤其是
有机废物)。然后，这些元素可能会堵塞旋转式过滤器(例如鼓式过滤器)，所述旋转式过滤
器旨在在使用水之前对其进行过滤，例如但不限于电力生产设施中的蒸汽产生回路的上
游。
[0002] 这些元素(在下文中称为“堵塞剂”(植物或活生物))的确可以被吸入任何用于供应水的设施的泵站中。这可能是发电厂、饮用水生产厂或其它设施的泵站。在以上任何类型
的水泵应用中，存在堵塞设置在泵站中的旋转式过滤器(鼓式过滤器或链式过滤器)的风
险。
[0003] 在以上类型的某些应用中，不允许旋转式过滤器下游流量不足或为零的风险。并且，期望一种以下解决方案：根据过滤器下游的最小期望流速、或至少过滤器下游的最大容
许压力损失或至少此压力损失随时间推移的可容许的演变来优化过滤器的设计。
[0004] 在设计水泵站的移动过滤系统时，不存在对这些过滤器的设计参数(过滤器大小、转速、洗涤系统)进行定量定义的，顾及到这些堵塞剂带来的影响的精确的定量技术。
[0005] 本发明改善了此情况。
[0006] 为此，本发明提出了一种用于管理用于泵送源自可能含有杂质的自然环境的水的设施的方法，所述设施旨在利用至少一个旋转式过滤器来净化所泵送的水，而所述杂质中
的至少一些杂质至少部分地堵塞所述旋转式过滤器。
[0007] 具体地，所述方法包括至少基于以下来估计由杂质对所述旋转式过滤器的堵塞引起的压力损失随时间推移的演变：
[0008] ‑与所述自然环境相关的数据，
[0009] ‑所述旋转式过滤器的尺寸，以及
[0010] ‑与所述过滤器上游的至少一个水位Namont以及所述过滤器下游吸入的水的流速Qaspiré有关的局部测量结果。
[0011] 然后，这种实施例使得能够预期随时间推移的所述压力损失，以优化所述泵送设施的操作条件，或甚至设计出针对此泵送设施的要求而优化的过滤器。
[0012] 在一个实施例中，所述旋转式过滤器包括至少一个圆筒：
[0013] ‑所述至少一个圆筒的给定半径为RF，
[0014] ‑所述至少一个圆筒相对于给定参考绕高度为Naxe的给定轴线旋转，
[0015] ‑并且所述至少一个圆筒的给定宽度为LF，所述给定宽度被限定为平行于所述至少一个圆筒的旋转轴线，
[0016] 并且，所述压力损失由所述过滤器上游的水位Namont与所述过滤器下游的水位Naval之差限定，所述压力损失随时间推移的所述演变是根据由下式给出的所述过滤器下游的水
位随时间推移的演变推断的：
[0017]
[0018] 其中：
[0019] ‑Naval是根据所述过滤器下游吸入的水的流速Qaspiré的测量结果来测量或推断的所述过滤器下游的水位；
[0020] ‑Naxe是相对于与所述过滤器下游的水位Naval相同的参考而确定的所述过滤器的旋转轴线的高度水平，
[0021] ‑Qfiltre是穿过所述过滤器的水的流速。
[0022] 在此实施例中，穿过所述过滤器的水的流速Qfiltre可以至少根据所述过滤器上游的水位Namont来计算，所述过滤器上游的水位是相对于与所述过滤器下游的水位Naval相同的
参考而确定的，对于：
[0023]
[0024] 如下：
[0025]
[0026] ，其中：
[0027] ‑g是万有引力常数，
[0028] ‑f1是部分堵塞的过滤器的孔隙率并且与所述过滤器中的水的空隙容积与所述过滤器的总容积的比率相对应，
[0029] ‑Sem与未浸没过滤器区段相对应，并且对于鼓式过滤器由下式表示：
[0030] Sem＝LFRF(αNamont‑αNaval)，其中
[0031] 并且
[0032] ‑对于Re＜400，kRe＝1+0.7e‑0.0106Re并且对于Re＞400，kRe＝1，Re是特定于所述过滤器的雷诺数。
[0033] 其中 的这种情况下(与所述压力损失相比问题更少)给出了针对以下详细描述呈现的Qfiltre的略有不同的等式。
[0034] 在以上表达式中，所述比率f1由f1＝f0(1‑Tc1)给出，f0是与干净过滤器的孔隙率有关的常数，并且Tc1是当前时刻所述过滤器的堵塞程度。
[0035] 在所述过滤器包括交叉金属丝网筛的实施例中，以上提及的雷诺数由下式给出：
[0036] 其中
[0037] 其中：
[0038] ‑δ0是干净无杂质的过滤器金属丝的直径，
[0039] ‑a0是过滤网的大小，
[0040] ‑Vam是所述过滤器上游的水的速度，并且
[0041] ‑Tc1是当前时刻所述过滤器的堵塞程度，
[0042] ‑v是在25℃的温度下的值为0.000001m2/s的水的运动粘度。
[0043] 在此实施例中，所述过滤器上游的水的速度是通过下式计算的：
[0044] 其中ξ是所述过滤器的压力损失系数。
[0045] 在一个实施例中，ξ是所述过滤器的浸没区段的所述压力损失系数，由适应于细滤网的Idel′chik相关性给出，并且由下式表示：
[0046]
[0047] 在一个实施例中，所述压力损失随时间推移的所述演变是当前时刻所述过滤器的至少一个堵塞程度和此堵塞程度的时间导数的函数，并且其中所述方法包括求解与此堵塞
程度相关联的至少一个微分方程以及利用与所述自然环境相关的数据的步骤。
[0048] 在所述旋转式过滤器包括一个或多个喷嘴的实施例中，所述一个或多个喷嘴用于洗涤所述过滤器的内部外围上的水射流冲击的至少一个给定点，所述旋转式过滤器的堵塞
程度是针对三种不同类型的过滤器部分来确定的：
[0049] ‑所述过滤器的第一部分，所述第一部分浸没在水中并且与对应于所述过滤器的整个浸没外围的接触表区域S0上的杂质接触，所述第一部分具有第一堵塞程度Tc1，
[0050] ‑所述过滤器的第二部分，所述第二部分未被浸没并且在外围表面区域之上延伸，所述外围表面区域最多延伸到来自所述洗涤喷嘴的所述冲击，所述第二部分表示为Sbal‑
S0，其中Sbal≤Sαlav，Sαlav是所述过滤器的延伸到来自所述喷嘴的所述冲击的所述未浸没部
分的外围表面区域，此第二部分与第二堵塞程度Tc2相对应，
[0051] ‑所述过滤器的第三部分，所述第三部分未被浸没并且与所述第二部分互补，所述第三部分由Smax‑Sbal限定，Smax是所述过滤器的总外围表面区域，此第三部分与第三堵塞程
度Tc3相对应。
[0052] 因此，来自所述喷嘴的所述射流在所述过滤器的内部外围壁上的所述冲击的位置在压力损失演变计算中起作用，如以下参考图2C进一步详细解释的。
[0053] 在一个实施例中，所述第一堵塞程度、所述第二堵塞程度和所述第三堵塞程度通过以下微分方程关联：
[0054]
[0055]
[0056]
[0057] 其中：
[0058] ‑C是在所述自然环境中测量的堵塞杂质的浓度，
[0059] ‑P是特定于所述堵塞杂质的类型的堵塞能力，
[0060] ‑LF是被限定成平行于所述过滤器的旋转轴线的所述过滤器的宽度，
[0061] ‑Qfiltre是穿过所述过滤器的水的流速，
[0062] ‑VF是作为压力损失阈值的函数的所述过滤器的转速，
[0063] ‑rlav是与测量数据或统计数据相对应的由所述喷嘴实现的过滤器洗涤效率。
[0064] 所述过滤器的所述转速VF通常是由所述过滤器根据所述测量的压力损失阈值进行操作而施加的。
[0065] 在此实施例中，对于鼓式旋转式过滤器：所述鼓式旋转式过滤器的给定半径为RF，所述鼓式旋转式过滤器相对于给定参考绕高度为Naxe的给定轴线旋转，并且所述鼓式旋转
式过滤器的给定宽度为LF，所述给定宽度被限定为平行于所述鼓式旋转式过滤器的旋转轴
线，
[0066] ‑浸没表面区域S0可以由下式给出：
[0067] Sam＝RF×LF×αNamont，其中
[0068] ‑所述总外围表面区域Smax＝2×π×RF×LF
[0069] ‑延伸到来自所述洗涤喷嘴的所述冲击的所述未浸没部分的所述外围表面区域可以由下式给出：
[0070] 其中：
[0071] ‑αlav是所述过滤器固有的角度值并且是来自所述喷嘴的所述冲击的位置的函数，
[0072] ‑Namont是根据所述过滤器下游吸入的水的流速Qaspiré的测量结果来测量或推断的所述过滤器上游的水位；
[0073] ‑Naxe是相对于与所述过滤器上游的水位Namont相同的参考而确定的所述过滤器的旋转轴线的高度水平。
[0074] 在包括通过链条连接的上部旋转圆筒和下部旋转圆筒的链式过滤器类型的所述旋转式过滤器的替代性实施例中：
[0075] ‑所述浸没表面区域S0可以由下式给出： 条件是
[0076] Namont＞Naxe，
[0077] 并且如果Namont＜Naxe，则由 给出
[0078] ‑并且在不考虑来自任何洗涤喷嘴的冲击的情况下，所述未浸没过滤器区段Smax‑S0可以由下式给出：
[0079] 其中Naxe＜Naval并且Namont＜NFsup，
[0080] 其中：
[0081] ‑RFsup是所述上部圆筒的半径，
[0082] ‑NFsup是上部旋转轴线的水平，并且
[0083] ‑αF是过滤表面相对于竖直线的倾斜度，
[0084] ‑Namont是根据所述过滤器下游吸入的水的流速Qaspiré的测量结果来测量或推断的所述过滤器上游的水位；
[0085] ‑Naval是根据所述过滤器下游吸入的水的流速Qaspiré的测量结果来测量或推断的所述过滤器下游的水位；
[0086] ‑Naxe是相对于与所述过滤器下游的水位Naval相同的参考而确定的所述过滤器的旋转轴线的高度水平。
[0087] 在一个实施例中，穿过所述过滤器的水的流速Qfiltre可以根据与以上提及的所述过滤器的所述第一部分中的所述第一堵塞程度Tcl相对应的所述过滤器的堵塞程度来计算。
[0088] 因此，在一个可能的实施例中，所述过滤器的尺寸可以根据所述设施在吸入水的流速方面的要求并且在预期来自所述自然环境的杂质堵塞所述过滤器的情况下进行选择。
[0089] 在另外的或可替代的实施例中，所述旋转式过滤器然后可以以作为估计的预期压力损失的函数的可变速度旋转，所述估计的预期压力损失是通过计算给定性质和尺寸的过
滤器的所述压力损失随时间推移的所述演变给出的。
[0090] 在另外的或可替代的实施例中，对于所述泵送设施的要求，所述过滤器下游吸入的水的流速可以根据给定性质和尺寸的过滤器的所述压力损失随时间推移的所述演变的
计算来预先估计。
[0091] 本发明还涉及一种包括用于实施根据本发明的方法的处理电路的装置。
[0092] 本发明还涉及一种计算机程序，所述计算机程序包括用于在处理器(和/或可能以非暂时性方式存储这种计算机程序的数据的存储介质)执行此程序时实施根据本发明的方
法的指令。
[0093] 通过阅读一些示例性实施例的以下详细描述并且通过检查附图，本发明的其它特征和优点将变得显而易见，其中：
[0094] ‑图1示意性地展示了包括旋转式过滤器FIL的过滤设施，
[0095] ‑图2A示意性地展示了鼓式过滤器类型的旋转式过滤器，
[0096] ‑图2B示意性地展示了链式过滤器类型的旋转式过滤器，
[0097] ‑图2C展示了由来自洗涤喷嘴的射流在旋转式过滤器的内部外围壁上冲击而形成的角度，所述旋转式过滤器在所展示的实例中为鼓式过滤器类型，
[0098] ‑图3A、3B和3C展示了压力损失演变计算中涉及的旋转式过滤器的外围的相应表面区域，
[0099] ‑图4展示了根据本发明的一个实施例的方法的主要步骤，
[0100] ‑图5展示了用于实施根据本发明的方法的装置。
[0101] 在以下示例性实施例中呈现的本发明提出对旋转式过滤器处的压力损失的演变(或旋转式过滤器下游的水位的演变)进行建模，此压力损失由堵塞剂的流入产生。
[0102] 以下提出了用于过滤器操作的不同配置以及针对每种配置的堵塞动力学。
[0103] 在下文中，假设在泵站上游的过滤设施中，旋转式过滤器如以下参考图1所呈现的来操作。过滤设施可以包括格栅GR，以保留存在于天然水道中的一些堵塞剂COL(如藻类、树
叶或其它植物碎屑或更通常的是其它碎屑)，旋转式过滤器FIL的上游包括带有开口的外壁
(例如金属丝网)以在允许水通过内部的同时将碎屑COL保持在外部。取决于穿过过滤器的
压力损失，过滤器以更高或更低的速度旋转并且由上游的水位高度Namont减去下游的水位高
度Naval给出。通常，上游水位是恒定的(至少在此处考虑的时间段内)，而下游水位则根据过
滤器的堵塞而变化(下降)。可以通过压力传感器等来测量下游水位。如果注意到下游水位
低于预定阈值，则旋转过滤器FIL(通常像洗衣机滚筒一样)以释放尽可能多的堵塞剂，根据
测量的压力损失来调整过滤器的转速。另外，洗涤喷嘴BL通常设置在过滤器内部的固定位
置处(并且不像过滤器那样旋转)，使得其在过滤器内部外围上的冲击点与水平线形成恒定
角度，如以下参考图2C将进一步看到的。然后，来自喷嘴的喷射使得能够从内部清洁过滤器
的外周表面，以尽可能多地去除堵塞剂。
[0104] 因此，可以如图1所展示的收集过滤器内部的经过过滤的水，然后由例如发电厂的冷却回路等泵POM泵送。
[0105] 另外，旋转式过滤器的两种类型区别如下：
[0106] ‑如图2A所展示的第一类型的圆柱形(圆形横区段)鼓式过滤器类型，并且
[0107] ‑如图2B所展示的第二类型称为链类型，并且包括由链条连接的两个滚筒(下部和上部)。
[0108] 在下表中解释了两个图2A和2B中出现的参数。
[0109]
[0110] 表1：过滤器的几何参数
[0111] 以下参考图2C解释了以上表1中与洗涤喷嘴的位置有关的最后一个参数。洗涤喷嘴位于滚筒内部的固定位置处。所述洗涤喷嘴的水射流在过滤器内部外围壁上的冲击点相
对于定向的水平轴线(图的右侧)形成角度αlav(以度为单位)。图2C展示了在鼓式过滤器的
情况下的角度αlav的定义。在链式过滤器的情况下，洗涤喷嘴在上部圆筒上。
[0112] 以下表2中呈现了过滤器操作参数。
[0113]参数 单位 符号
上游水位 m Namont
下游水位 m Naval
3
由泵站吸入的流速 m/s Qasp
针对每个压力损失阈值的过滤器的转速 m/s VF
用于触发转速的压力损失阈值 m Dp
洗涤效率 ‑ rlav
[0114] 表2：过滤器操作参数
[0115] 进一步限定表征堵塞剂的参数：
[0116] ‑C，堵塞剂的浓度(kg/m3)
[0117] ‑P，所研究的堵塞剂类型的堵塞能力。此参数P被限定为堵塞1m2的表面所需的质2 2
量。其也可以被看作是lm均匀分布的堵塞剂的重量。其以kg/m表示。
[0118] 堵塞剂浓度数据与当地水生生物学知识相关。可以例如基于在先前堵塞事件期间在泵站去除的生物量来对其进行估计。这些数据可以从观察中得出并且也可以测量得出。
因此，其可以表征以下等式中所呈现的参数。
[0119] 以下呈现了在堵塞剂流入的情况下计算压力损失演变的原理。
[0120] 旋转式过滤器下游的水位变化具体地取决于以下进行计算：
[0121] ‑过滤器的类型及其操作，
[0122] ‑上游水位，
[0123] ‑堵塞剂的浓度，
[0124] ‑以及下游吸入的流速。
[0125] 另外，采用了其中堵塞剂的浓度均匀地分布在水柱中的假设。此浓度可以随时间推移而恒定或可变。
[0126] 还采用了到达过滤器的所有堵塞剂都保留在其上的假设。
[0127] 接下来，考虑瞬变现象以正确地表示堵塞动态是有利的：穿过过滤器的流速不一定等于堵塞期间的吸入流速。
[0128] 在初始时间，假设清洁过滤器和稳定状态，使得最初穿过过滤器的流速与吸入流速相对应。然后，用耦合方程组计算三个序列中时间t与t+dt之间的主要变量的变化，如下：
[0129] ‑首先计算堵塞程度的演变，
[0130] ‑由此根据堵塞程度和压力损失来确定穿过过滤器的流速，
[0131] ‑然后使用因此计算出的过滤器的流速和吸入流速来求解下游水位的演变。此变化取决于过滤器的类型(鼓式过滤器或链式过滤器)。
[0132] 下表中总结了计算符号。
[0133]
[0134] 表3：计算变量
[0135] 如以下呈现的计算堵塞程度的演变。
[0136] 根据定义，堵塞程度Tc是堵塞的流体体积与初始流体体积的比率。换句话说，通过将堵塞的过滤器的孔隙率表示为f，将所述孔隙率表示为流体的空隙体积与所考虑的过滤
器总体积(更多固体流体)的比率，得到：
[0137] f0是清洁过滤器的孔隙率。
[0138] 固体体积与清洁过滤器体积和累积堵塞体积的总和相对应，如下：
[0139]
[0140] 对于“二维”过滤器(格栅、鼓式过滤器或链式过滤器类型的组合过滤器)，这些体积比率与区段比率相对应。
[0141] 因此，以下使用与设计过滤器的物理尺寸相对应的区段来描述数值模型。
[0142] 具体地，试图比较过滤器的堵塞表面区域与流体表面区域。
[0143] 流体区段由下式给出：
[0144] Sfluide＝f0(Stot‑Scolm)
[0145] 总堵塞区段与堵塞程度相关：
[0146]
[0147] 此外，堵塞剂的特征在于由P表示的“堵塞能力”并且使得：
[0148] Mc＝PScolm＝TcPStot，Mc是过滤器上堵塞剂的质量(以kg为单位)并且与堵塞表面相关联。
[0149] 根据图3A、3B和3C的三种情况来限定与堵塞剂接触的鼓式过滤器或链式过滤器底部的表面区域。
[0150] 在图3A中，仅浸没在水中的滚筒部分与堵塞剂接触，因此限定了外部接触表面区域S0。在图3B中，堵塞剂的量使得堵塞剂的一部分还保持与可变表面区域Sbal之上的未浸没
滚筒的一部分的接触：
[0151] Sbal＝S0+LFVFt，因此与大于S0但小于滚筒Smax的总外表面面积的扫掠表面区域相对应。在示出了最不利的情况的图3C中，滚筒的整个表面区域Smax完全被堵塞剂覆盖。
[0152] 可以更具体地考虑相对于旋转式过滤器的第一一次或多次旋转(通常在若干次旋转之后建立稳定状态之前的过滤器的第一一次或多次旋转)中的扫掠表面区域Sbal的三个
相邻且互补的区：
[0153] ‑在堵塞S0区Zl中，与堵塞程度Tc1相对应，
[0154] ‑在与洗涤喷嘴之前的Sbal‑S0相对应的中性区Z2中，其中Sbal≤Sαlav，与堵塞程度Tc2相对应(Sαlav与延伸到洗涤喷嘴的冲击并且由角度αlav限定的外围表面区域相对应)，
[0155] ‑在互补区Z3中，Smax‑Sbal(因此包含在针对Sbal＞Sαlav的由喷嘴清洁的区域内)，与堵塞程度Tc3相对应。
[0156] 限定了在时间间隔dt内在区之间交换的质量：
[0157] ‑dMl，从区Z3朝着区Zl进入的质量为：
[0158] dM1＝PTc3LFVFdt
[0159] ‑dM2，离开区Zl朝着区Z2的质量为：
[0160] dM2＝PTc1LFVFdt
[0161] ‑dM3，离开区Z2朝着区Z3的质量为：
[0162] dM3＝PTc2LFVFdt
[0163] 在区Zl(S0区段)中，将随时间推移的质量平衡记为：
[0164] dMentrant+dM1‑dM2＝PdScolm1
[0165] ‑进入过滤器的质量为：
[0166] dMentrant＝CQfiltredt
[0167] ‑将区Z1中的堵塞区段及其导数记为：
[0168] Scolm1＝Tc1S0
[0169] dScolm1＝S0dTc1
[0170] 获得第一关系式：
[0171]
[0172] 在区Z2(Sbal‑S0区段)中，将随时间推移的质量平衡记为：
[0173] dM2‑dM3＝PdScolm2
[0174] ‑在此区中，将堵塞部分及其导数记为：
[0175] Scolm2＝Tc2(Sbal‑S0)
[0176]
[0177] 获得第二关系式：
[0178]
[0179] 更准确地说，根据来自喷嘴的射流的冲击位置来更详细地限定区Z2，并且只要Sbal＜Sαlav，以上关系式就有效。
[0180] 另一方面，对于Sbal之Sαlav，得到：
[0181] Sbal＝Sαlav
[0182]
[0183] 在区Z3(Smax‑Sbal区段)中，对于Sbal≤Sαlav，堵塞程度为零并且记为：
[0184]
[0185]
[0186] ‑对于Sαlav＜Sbal≤Smax，将随时间推移的质量平衡记为：
[0187] M3‑dM1‑dMsortt＝PdScolm3
[0188] ○由于对区Z3进行清洁而离开过滤器的质量(dMsortt)为：
[0189] dMsortt＝PTc2dSnet＝rlavPTc2LFVFdt
[0190] ○将区Z3中的堵塞区段及其导数记为：
[0191] Scolm3＝Tc3(Sbal‑Sαlav)
[0192]
[0193] 获得第三关系式：
[0194]
[0195] ‑并且，对于Sbal≥Smax：
[0196] Sbal＝Smax
[0197]
[0198] 接下来，穿过过滤器的流速可以计算如下。其由浸没流速和未浸没流速构成，第一个由Idel′chik等式描述，第二个由通过根据使用已知孔隙率的多孔板在通道中进行的测
试而建立的定律描述(根据申请人的内部研究)：
[0199] Qfiltre(t)＝Qim(t)+Qem(t)
[0200] 可以根据压力损失来通过Idel′chik定律计算淹没流速Qim(t)。
[0201] 在堵塞程度沿浸没上游表面变化的一般情况下，浸没流速为：
[0202]
[0203] 根据过滤器的压力损失来通过下式计算区段dS上游的局部速度：
[0204]
[0205] 系数ξ取决于区段dS中的局部堵塞程度。如果浸没区段之上的堵塞程度是均匀的，则系数ξ是均匀的并且获得：
[0206] 此术语Sam最终仅代表先前所计算的区段S0。
[0207] 对于堵塞的过滤网，系数ξ由适应于细滤网的Idel′chik相关性给出：
[0208] 此处具体地，f＝f1＝f0(1‑Tc1)，f0是与清洁过滤器的孔隙率有关的常数，并且Tc1是第一区Zl中过滤器的堵塞程度。
[0209] 系数kRe取决于由下式计算的雷诺数Re：
[0210] 其中v是水的运动粘度(在25℃下恒定并且等于0.000001m2/s)并且δ是堵塞的滤网的金属丝的直径，使得：
[0211]
[0212] 对于Re＜400：kRe＝1+0.7e‑00106 Re
[0213] 对于Re＞400：kRe＝1
[0214] 此处Tc＝Tc1。
[0215] 对于鼓式过滤器，上游区段为：
[0216]
[0217] 对于链式过滤器，上游区段为：
[0218] ‑如果Namont＞Naxe
[0219]
[0220] ‑如果Namont＜Naxe
[0221]
[0222] 然后，通过根据使用已知孔隙率的多孔板在通道中进行的测试建立的相关性来计算未浸没流速(基于申请人的内部结果)。根据这些测试，通过下式将未浸没流速计算作为
流体(fSem)通过的横区段积和压力损失的函数：
[0223]
[0224] 其中流量系数Cd取决于：
[0225] ‑对于x＜0.65：Cd＝‑0.8x3+2.48x2‑2.60x+1.51
[0226] ‑对于x＞0.65：Cd＝0.6
[0227] 因此，将Cd的这些值计算作为取决于压力损失(Namont‑Naval)的参数x的函数。这些公式是通过由申请人进行的在通道中进行的试验建立的并且因此已经证明取决于x＞0.65
还是x＜0.65考虑两个公式是优选的。
[0228] 对于鼓式过滤器，未浸没过滤器的区段为：
[0229] Sem＝LFRF(αNamont‑αNaval)
[0230] 无论对于Namont或N＝Naval
[0231] 对于链式过滤器，在最常见的Naxe＜Naval并且Namont＜NFsup的情况下，未浸没过滤器的区段为：
[0232]
[0233] 在链式过滤器中，洗涤喷嘴在水平NFsup上方的上部圆筒上并且如鼓式过滤器在计算中使用射流冲击的位置，但其中公式略有不同。
[0234] 然后通过求解下游体积的质量平衡方程来确定下游水位的演变：
[0235]
[0236] ‑其中在每个瞬间均根据堵塞程度和压力损失来计算过滤器的流速Qfiltre，
[0237] ‑并且随时间推移反复地测量吸入流速Qaspiré。
[0238] 将与水平N相关联的过滤器下游的容积表示为Vf(N)。
[0239] 对于鼓式过滤器，此体积等于：
[0240]
[0241] 其中：
[0242] 因此：
[0243] 对于链式过滤器，此体积取决于HN和αF：
[0244] HN＝N‑Naxe
[0245]
[0246] 对于Namont＞Naxe，此体积为 并且得到：
[0247]
[0248] 对于Namont＜Naxe，下游体积及其导数的计算方法与鼓式过滤器相同。
[0249] 总结如下，如图4所展示的，在第一步骤S0中，考虑过滤器(鼓式或链式)固有的尺寸和功能参数：
[0250] ‑RF，过滤器的半径(被认为是链式过滤器的下部圆筒)，
[0251] ‑LF，过滤器的宽度(平行于其旋转轴线)，
[0252] ‑Naxe，过滤器的旋转轴线的水平(相对于泵送设施的参考水平的高度，以米为单位)，
[0253] ‑a0，清洁过滤器的滤网大小(在没有堵塞剂的情况下)，以米为单位，
[0254] ‑δ0，清洁过滤器的金属丝直径(以m为单位)，
[0255] ‑f0，清洁过滤器的滤网孔隙率
[0256] ‑αlav，以相对于过滤器的旋转轴线的弧度表示的过滤器中洗涤喷嘴的位置，
[0257] ‑rlav，喷嘴的过滤器洗涤效率。
[0258] 在步骤S1中，还考虑了特定于泵送场地的特性：
[0259] ‑C，可以在泵送场地(当地水生生物学)测量的水中的堵塞剂的浓度(kg/m3)，
[0260] ‑P，所研究的堵塞剂类型的堵塞能力。
[0261] 此外，在步骤S2中，还考虑了以下变量，所述变量易于随时间推移而变化并且其构成了其后进行的用于计算的输入参数：
[0262] ‑VF(t)，针对每个压力损失阈值的过滤器的转速，
[0263] ‑Qasp(t)，以m3/s为单位的由泵站吸入的流速，
[0264] ‑Namont(t)，相对于泵送设施中的参考的上游水位(以米为单位)。
[0265] 实际上，由用于水平Namont和Naval的传感器在不同时间提供压力损失的局部测量结果。基于此测量结果，鼓式过滤器取决于周期性地测量的压力损失以更低或更高的速度VF
(t)旋转。
[0266] 实际上，过滤器(鼓式或链式)的正常操作是根据由传感器测量(测量上游水位和下游水位)的压力损失进行旋转。更具体地，在获得以上提及的压力损失测量结果之后，转
速VF取决于超过压力损失阈值的值。因此，转速可以作为依次到达的相应的压力损失阈值
Dp的函数依次触发。
[0267] 在S3步骤中，根据此推导出用于计算的以下变量并且所述变量特定于如现场安装的过滤器的数据(鼓式的情况下)：
[0268] 表面区域Smax＝2×π×RF×LF
[0269] 表面区域S0＝SNamont＝RF×LF×αNamont
[0270] 表面区域
[0271] αNamont以弧度为单位并且αlav以度为单位，
[0272] 其中
[0273] 接下来，在初始化步骤S4中，在时间t＝0处确定初始条件，在所述时间假设过滤器是清洁的并且流速处于稳定状态，如下：
[0274] Tc1(0)＝Tc2(0)＝Tc3(0)＝0
[0275] Qfiltre(0)＝Qasp(0)
[0276]
[0277] 然后，在步骤S5中，在每个时间增量(dt)处：
[0278] ‑扫掠表面区域Sbal计算如下：
[0279] 其中Sbal(t)＝S0+LFVFt并且dSbal(t)＝LFVFdt对于Sbal＜Smax
[0280] 并且另外
[0281] 然后在步骤S6中，可以假设具有三个随时间推移的微分方程的方程组，所述方程组与堵塞程度的演变有关，如以上所看到的：
[0282]
[0283] 对于Sbal≤Sαlav，
[0284] dTc3＝0
[0285] 对于Sbal＞Sαlav
[0286]
[0287] 对于Sαlav＜Sbal＜Smax，
[0288]
[0289] 对于Sbal＞Smax
[0290]
[0291] 在此方程组中，穿过过滤器的流速Qfiltre最初由Qfiltre(0)＝Qasp(0)给出并且此值Qasp(0)是其最初在时间t＝0+dt处用来计算dTc1的值，
[0292] 以及用来计算dTc1、dTc2和dTc3，并且然后根据此推导出Tc1、Tc2和Tc3的新值。
[0293] 接下来，在步骤S7中，根据堵塞率Tc1和压力损失推导穿过过滤器的新流速，如下，此处举例来说，在最常见的实施例中，其中过滤器是旋转滚筒，使得：
[0294]
[0295] 其中：f1＝f0(1‑Tc1)，浸没过滤器的第一区的孔隙率并且因此根据Tc1推导所述孔隙率，并且：
[0296] ·Sem＝LFRF(αNamont‑αNaval)其中
[0297] ·对于Re＜400kRe＝1+0.7e‑0.0106Re；对于Re＞400kRe＝1，Re是特定于过滤器的雷诺数，使得：
[0298] 其中 并且
[0299] 然后可以将以此方式发现的穿过过滤器的流速Qfiltre的值再注入第一方程dt以确定下一时间t+dt的dTc1，以及dTc1、dTc2和dTc3、，
并且从那里再次确定Tc1、Tc2和Tc3。
[0300] 最后，也可以在步骤S8中通过下式给出根据过滤器的流速和吸入流速的下游水位随时间推移的演变：
[0301]
[0302] 本发明可以用于：
[0303] ‑优化旋转过滤系统的操作，例如通过选择不同的压力损失阈值和对应的速度，以预测未来的压力损失并且立即确定最合适的速度；
[0304] ‑模拟过去的堵塞事件，以完善以上方程，具体地但不仅限于特定于场地水生生物学的参数(以上步骤S1)；
[0305] ‑确定用于需要给定最大下游流量的新泵送设施的旋转式过滤器的理想尺寸，并且考虑泵送场地的水生生物学。
[0306] 因此，图4的方法可以用于优化旋转式过滤器的操作(微调用于触发转速的阈值、用于降低泵送流速的阈值或其它值)或用于确定新泵站的旋转式过滤器的大小(扩建现有
泵站，翻新或创建新泵站)。
[0307] 此方法的优点如下：
[0308] ·对压力损失演变的估计进行定量，
[0309] ·其取决于特定场地，
[0310] ·其在资源和计算时间方面是物美价廉的，
[0311] ·其可以应用于在电力生产中为蒸汽发生器泵送水，但是更通常地应用于涉及用于从自然环境(大海、湖泊或河流)中泵送水的泵站的任何工业(化学工业或水处理)，或更
通常地应用于旋转式过滤器设计工业。
[0312] 此方法以及以上相关联的计算可以通过计算机程序来实施，其中通用算法遵循以上参考图4所呈现的流程图。因此，本发明还涉及这种计算机程序，以及如图5所展示的计算
装置并且包含处理电路DIS，所述处理电路DIS包括：
[0313] ‑输入接口IN，所述输入接口IN用于接收压力损失测量结果(或至少下游水位)或用于设计旋转式过滤器的这种测量结果的模拟值，
[0314] ‑存储器MEM，所述存储器MEM用于存储本发明意义内的计算机程序的指令(以及可能的过滤器参数数据和/或中间计算数据)，
[0315] ‑处理器PROC，所述处理器PROC与存储器MEM协作并且执行以上方法，以及
[0316] ‑输出接口OUT，所述输出接口OUT用于检索压力损失演变的计算，并且可能地警报信号AL，所述警报信号AL用于通知由于基于此演变的估计的时间长度之后的预期压力损失
过高而引起的泵送系统POM的供水不足，或推荐数据RECO，所述推荐数据RECO在模拟输入IN
的测量数据MES的情况下用于优化过滤器设计。
