[0048] 下面结合附图以及具体实施例对本发明作进一步的说明，但本发明的保护范围并不限于此。
[0049] 通过瞬态压力传感器监测静压值。发生空化时，核主泵进口的压力脉动频谱会发生显著的变化：在正转核主泵工况下，核主泵进口主频即轴频的幅值会出现先稍有增大后快速降低的趋势，而在低频区则会出现快速增强的低频信号；在正转逆流工况下，核主泵进口的压力脉动主频幅值呈递减趋势，低频信号幅值则持续增大。由于空化状态会使上述频率的幅值发生特定地变化，因此可以利用这些压力脉动信号的幅值构造函数来对核主泵内的空化状态进行判断。
[0050] 如图1所示，本发明所述的核主泵空化状态监控方法，包括如下步骤：
[0051] 压力数据采集：采集核主泵进口的压力脉动数据；选取低频目标频率F1＝0.42Hz，F2＝2.92Hz，F3＝17.92Hz，F4＝19.17Hz，主频信号fn为48.3Hz。
[0052] 确定低频信号变量 在压力脉动数据中，选取其中幅值变化明显的低频信号作为低频目标频率Fj；对所述低频信号进行量化，构造低频目标频率为Fj时的低频信号变量[0053]
[0054]
[0055] 其中：
[0056] Fj为选取的第j个低频目标频率，1≤j＜N，N为自然数；
[0057] 为当低频目标频率Fj时，第i次试验所获取的频域幅值，1≤i＜M，M为自然数；
[0058] 为当低频目标频率Fj时的低频信号变量；
[0059] 为未发生空化时，低频目标频率Fj采集n次的平均幅值；n为采集次数，n不小于200次。
[0060] 确定主频信号变量Bi：在压力脉动数据中，选取其中主频信号fn信号量化，构造主频信号变量Bi：
[0061]
[0062]
[0063] 其中：
[0064] bi为主频信号fn在第i次试验所获取的频域幅值；
[0065] Bi为主频信号变量；
[0066] B0为未发生空化时，主频信号fn采集n次的平均幅值；n为采集次数，n不小于200次。
[0067] 构造判断空化状态函数σi：利用低频信号变量 与主频信号变量Bi构造判断空化状态函数σi为：
[0068]
[0069] 如图2所示，为低频目标频率F1＝0.42Hz、F2＝2.92Hz、F3＝17.92Hz、F4＝19.17Hz下的 曲线，主频信号fn为48.3Hz下的Bi曲线以及构造函数σi的数值曲线，其中， 曲线与Bi曲线共用左边的纵坐标，构造函数σi的数值曲线使用右边的纵坐标。
[0070] 构造阙值函数τi：根据判断空化状态函数σi构造阙值函数τi为[0071] 其中： 为未空化时σi的平均值，此值可以通过数值模拟的方式获得；
[0072] 求解临界空化点：通过插值法求得临界空化点，当检测出的σi值大于临界空化点的值时，判定当前为空化状态。
[0073] 所述 和bi的值均采用量纲均一化处理，即：
[0074]
[0075]
[0076] 其中：
[0077] pa为低频目标频率Fj时的瞬态静压值，KPa；
[0078] 为低频目标频率Fj时的平均静压值，KPa；
[0079] pb为主频信号fn时的瞬态静压值，KPa；
[0080] 为主频信号fn时的平均静压值，KPa；
[0081] ρ为流体密度，kg/m3；
[0082] u为叶轮转速，rad/s。
[0083] 最终得到如图2所示的阙值函数曲线图，通过插值法求得τi＝2.5％时σi的值，此数值即临界空化点。求得此时构造函数值σi＝-0.0284，利用插值法求得此时进口压力约为-69kpa，这与试验测得的临界空化点(P＝-70kpa)相差不大，证明了本发明是具有可实施性的。
[0084] 目标频率选取的越多，最终的预测结果越接近实际情况。
[0085] 本发明所述的核主泵空化状态监控系统，包括瞬态压力传感器和控制系统，所述核主泵的进口位置安装瞬态压力传感器，所述控制系统接收瞬态压力传感器输出的信号；所述控制系统根据本发明所述的核主泵空化状态监控方法，判断核主泵是否出现空化状态。所述控制系统可以内置数据分析模块和判定模块，这个可以根据控制系统自行设计。
[0086] 所述实施例为本发明的优选的实施方式，但本发明并不限于上述实施方式，在不背离本发明的实质内容的情况下，本领域技术人员能够做出的任何显而易见的改进、替换或变型均属于本发明的保护范围。
