非能动安全壳空气冷却系统运行特性分析方法

1.本发明属于核反应堆安全壳事故现象分析技术领域，具体涉及一种非能动安全壳空气冷却系统运行特性分析方法。
技术背景
2.安全壳是核电厂纵深防御的最后一道屏障，对防止放射性物质进入环境具有重要作用。多用途模块式小型反应堆采用了非能动安全壳空气冷却系统，该系统内层为钢制安全壳，外层为混凝土屏蔽构筑物，中间环形夹层是空气冷却通道。当核反应堆一回路发生破口事故时，非能动安全壳空气冷却系统通过辐射换热、气体冷却、液膜冷却等方式将安全壳内部热量导出至大气环境，以保持安全壳的完整性。
3.目前已有三维流体力学计算软件可以对非能动安全壳空气冷却系统进行模拟计算，但是需要消耗大量的计算资源且耗时较长。目前缺少能够快速、高效、精确分析非能动安全壳空气冷却系统运行特性的计算程序。


技术实现要素：

4.为填补上述研究空白，本发明提供了一种非能动安全壳空气冷却系统运行特性分析方法，可以快速计算非能动安全壳空气冷却系统运行特性，对优化空气冷却系统设计、提升核反应堆安全性具有重要意义。
5.为实现上述目的，本发明采用的技术方案是：
6.一种非能动安全壳空气冷却系统运行特性分析方法，包括如下步骤：
7.步骤1：初始化钢制安全壳内部和空气冷却系统各隔间温度、压力、气体成分以及钢制安全壳壁面温度、液膜的温度和厚度；
8.步骤2：基于钢制安全壳内部热工水力状态判断空气冷却系统是否启动，具体包括如下内容：
9.1)当钢制安全壳内部压力未达到启动限值，判定空气冷却系统不启动；
10.2)当钢制安全壳内部压力达到启动限值，判定空气冷却系统启动，开始工作，并记录压力达到限值时间点tn，设定时间步长δτ；
11.步骤3：利用多节点安全壳分析程序计算钢制安全壳内部隔间与钢制安全壳内壁面换热，具体包括如下内容：
12.1)基于tn时刻钢制安全壳内部温度、压力和气体成分以及tn时刻钢制安全壳壁面温度计算δτ时间步长内两者之间的换热量q
in
；
13.2)基于tn时刻钢制安全壳内的热工水力状态、质量源项和δτ时间步长内钢制安全壳内部和钢制安全壳内壁面的换热量，计算该时间步长结束t
n 1
时刻钢制安全壳内部的温度、压力、气体成分热工水力状态；
14.步骤4：计算空气冷却系统隔间与钢制安全壳外壁面换热及其热工水力状态，具体包括如下内容：
15.1)计算空气冷却系统隔间与钢制安全壳外壁面换热
16.i.基于tn时刻空气冷却系统隔间气体温度、钢制安全壳壁面温度和液膜厚度计算δτ时间步长内空气冷却系统隔间气体和钢制安全壳外壁面之间的辐射换热量
[0017][0018]
式中：
[0019]qrad
——δτ时间步长内空气冷却系统隔间气体和钢制安全壳外壁面之间的辐射换热量；
[0020]ffilm
——液膜修正系数，若钢制安全壳外壁面未被液膜覆盖，则f
film
＝1，否则f
film
《1；
[0021]
a——钢制安全壳壁面面积；
[0022]
σ——玻尔兹曼常数；
[0023]
tg——tn时刻空气冷却系统隔间气体温度；
[0024]
t
wall
——tn时刻钢制安全壳壁面温度；
[0025]eg
——空气冷却系统隔间气体发射率；
[0026]ewall
——钢制安全壳壁面发射率；
[0027]
ii.i若钢制安全壳外壁面被液膜覆盖，则空气冷却系统隔间与钢制安全壳外壁面的对流换热量q
conv
＝0；
[0028]
基于tn时刻空气冷却系统隔间气体温度、气体成分和液膜温度及物性计算δτ时间步长内空气冷却系统隔间和液膜的传热传质：
[0029]
δτ时间步长内空气冷却系统隔间和液膜之间的传质即液膜蒸发量采用下式计算：
[0030][0031]
式中：
[0032]mevap
——δτ时间内液膜蒸发量；
[0033]hcd
——tn时刻液膜蒸发传热系数；
[0034]
t
film
——tn时刻液膜温度；
[0035]
δτ——时间步长；
[0036]
δh——汽化潜热；
[0037]
蒸发传热系数h
cd
采用uchida关系式进行计算：
[0038][0039]
式中：
[0040]
ρ
steam
——tn时刻空气冷系统隔间蒸汽密度；
[0041]
ρ
ng
——tn时刻空气冷系统隔间不凝性气体密度；
[0042]
δτ时间步长内液膜和空气冷却系统隔间的传热量采用下式计算：
[0043]qfilm_rb
＝m
evaphevap
ah
conv
(t
film-tg)δτ
[0044]
式中：
[0045]qfilm_rb
——δτ时间内液膜和空气冷却系统隔间传热量；
[0046]hevap
——tn时刻蒸汽比焓；
[0047]hconv
——液膜-隔间气体自然对流换热系数；
[0048]
因为液膜流速小，故研究液膜-隔间气体自然对流换热问题时将液膜视为固体，h
conv
采用如下方法计算：
[0049][0050]
nu——努塞尔数；
[0051]
l——钢制安全壳壁面特征尺寸；
[0052]
k——空气冷却系统隔间气体导热系数；
[0053]
努塞尔数nu的计算关系式与流动类型相关，流动类型通过雷诺数re2和格拉晓夫数gr的关系判断：
[0054][0055]
ρg——空气冷却系统隔间气体密度；
[0056]
v——空气冷却系统隔间气体流速；
[0057]
μ——空气冷却系统隔间气体动力粘度；
[0058][0059]
g——重力加速度；
[0060]
①
若re2＜1.0gr，为自然对流，努塞尔数nu采用以下关系式：
[0061][0062]
ra＝gr
·
pr
[0063]
pr——空气冷却系统隔间气体普朗特数；
[0064]
ra——空气冷却系统隔间气体瑞利数；
[0065][0066]
ν——空气冷却系统隔间气体运动粘度；
[0067]
α——空气冷却系统隔间气体热扩散系数；
[0068]cp
——空气冷却系统隔间气体定压比热容；
[0069]
②
若re2＞10gr，为强迫对流，努塞尔数nu采用以下关系式：
[0070][0071]
nu
forced
＝0.023re
0.8
pr
0.3
ꢀꢀꢀꢀꢀꢀꢀꢀ
re＞6000
[0072]
lnnu
forced
＝a1lnre b1ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
2000≤re≤6000
[0073]dh
——空气冷却系统隔间水力学直径；
[0074]
③
若1.0gr≤re2≤10gr，为混合对流，努塞尔数nu采用以下关系式：
[0075]
nu
mixed
＝[(re2/gr-1)/9][nu
forced-nu
natural
] nu
natural
[0076]
ii.ii若钢制安全壳外壁面无液膜覆盖，则m
evap
＝0，q
film_rb
＝0；
[0077]
δτ时间步长内空气冷却系统隔间与钢制安全壳外壁面的对流换热量基于空气冷却系统隔间气体温度、气体成分和钢制安全壳壁面温度计算：
[0078]qwall_rb
＝ah
conv
(t
wall-tg)δτ
[0079]
式中的对流换热系数h
conv
计算方法同步骤ii.i所述；
[0080]
2)基于tn时刻空气冷却系统隔间压力、温度、气体成分、δτ时间步长内空气冷却系统隔间与钢制安全壳外壁面的辐射换热量和对流换热量、空气冷却系统隔间与液膜之间的传热传质计算t
n 1
时刻空气冷却系统隔间压力、温度、气体成分；
[0081]
i.更新t
n 1
时刻空气冷却系统隔间蒸汽质量
[0082]mst_n 1
＝m
st
m
evap
[0083]
式中：
[0084]mst_n 1
——t
n 1
时刻空气冷却系统隔间蒸汽质量；
[0085]mst
——tn时刻空气冷却系统隔间蒸汽质量；
[0086]
ii.更新t
n 1
时刻空气冷却系统隔间温度
[0087][0088]
式中：
[0089]
t
g_n 1
——t
n 1
时刻空气冷却系统隔间气体温度；
[0090]cg
——tn时刻空气冷却系统隔间气体比热容；
[0091]
m——tn时刻空气冷却系统隔间气体质量；
[0092]
iii.更新t
n 1
时刻空气冷却系统隔间压力
[0093][0094]
p
n 1
——t
n 1
时刻气体压力；
[0095]
p
st_n 1
——t
n 1
时刻蒸汽分压，由蒸汽质量、气相体积和气体温度计算得到；
[0096]nmn_n 1
——t
n 1
时刻不凝气体摩尔数；
[0097]
rg——气体常数；
[0098]vg
——空气冷却系统隔间体积；
[0099]
步骤5：判断是否计算到最后一个隔间，若是，则到步骤7；否则，到步骤6；
[0100]
步骤6：计算钢制安全壳外壁面液膜传递过程，并对下一空气冷却系统隔间进行步骤3计算，具体包括如下内容：
[0101]
1)更新第i号钢制安全壳外壁面步长末t
n 1
时刻液膜质量
[0102]mfilm_n 1
＝m
film
m
in-m
out
m
sp-m
evap
[0103]mfilm_n 1
——t
n 1
时刻钢制安全壳外壁面液膜质量；
[0104]mfilm
——tn时刻钢制安全壳外壁面液膜质量；
[0105]min
——δτ时间步长内从i-1号钢制安全壳外壁面流入的液膜质量，若第i号钢制安全壳壁面上部无壁面，则m
in
＝0，否则m
in
为第i-1号壁面的m
out
值；
[0106]mout
——δτ时间内流入i 1号钢制安全壳外壁面的液膜质量，采用下式计算：
[0107][0108]
时刻钢制安全壳外壁面液膜平均流速，采用下式计算：
[0109][0110]
ρ
film
——tn时刻钢制安全壳外壁面液膜密度；
[0111]
μ
film
——tn时刻钢制安全壳外壁面液膜动力粘度；
[0112]
δ——tn时刻钢制安全壳外壁面液膜厚度，采用下式计算：
[0113][0114]
2)更新第i号钢制安全壳外壁面步长末t
n 1
时刻液膜温度
[0115]
i.基于tn时刻液膜温度和钢壳壁面温度计算δτ时间步长内液膜-钢制安全壳外壁面导热量：
[0116][0117]
式中：
[0118]qfilm-wall
——δτ时间步长内液膜和钢制安全壳外壁面导热量；
[0119]
δ
film
——tn时刻钢制安全壳外壁面液膜厚度；
[0120]
δ
wall
——钢制安全壳厚度；
[0121]kfilm
——tn时刻钢制安全壳外壁面液膜热导率；
[0122]kwall
——钢制安全壳热导率；
[0123]
ii.计算t
n 1
时刻液膜温度
[0124][0125]hfilm
——tn时刻钢制安全壳外壁面液膜比焓；
[0126]hfilm_n 1
——t
n 1
时刻钢制安全壳外壁面液膜比焓；
[0127]
已知更新后t
n 1
时刻的空气冷却系统隔间内压力和液膜比焓则利用插值法求出更新后钢制安全壳外壁面液膜温度；
[0128]
步骤7：计算钢制安全壳壁面温度，具体包括如下内容：
[0129]
1)计算得到各部分安全壳壁面的净吸热量：
[0130]
δq＝q
in-q
out
[0131]
式中：
[0132]
δq——壁面的净吸热量；
[0133]qin
——钢制安全壳内壁面从钢制安全壳内部隔间吸收的热量；
[0134]qout
——钢制安全壳外壁面向壳空气冷却系统散失的热量；
[0135]qout
＝q
film_wall
q
film_rb
[0136]
2)更新t
n 1
时刻各部分钢制安全壳壁面的温度：
[0137][0138]
t
wall_n 1
——t
n 1
时刻钢制安全壳壁面温度；
[0139]mwall
——钢制安全壳壁面质量；
[0140]cwall
——钢制安全壳壁面的定压比热容；
[0141]
步骤8：将t
n 1
时刻作为新的tn时刻，重复步骤3至7，直到达到指定计算时间。
[0142]
与现有技术相比，本发明具有如下优点：
[0143]
1.基于多节点安全壳分析程序开发，计算速度快、精度高；
[0144]
2.可以模拟非能动安全壳空气冷却系统的气体冷却、液膜冷却和辐射换热以及钢制安全壳外壁面液膜流动；
[0145]
3.可以分析空气冷却系统的热工水力状态，包括气体成分、温度和压力。
附图说明
[0146]
图1为非能动安全壳空气冷却系统运行特性分析计算流程图；
[0147]
图2为非能动安全壳空气冷却系统；
[0148]
图3为隔间划分示意图；
[0149]
图4为钢制安全壳外壁面-空气冷却系统隔间换热示意图。
具体实施方式
[0150]
下面通过结合附图对本发明进行详细描述。
[0151]
本发明安全壳内多组非能动热量导出系统运行特性分析方法，具体流程如图1所示，包括如下步骤：
[0152]
步骤1：如图2所示，非能动安全壳空气冷却系统内层是钢制安全壳，外层是混凝土屏蔽构筑物，中间环形通道为空气冷却系统隔间；初始化钢制安全壳内部和空气冷却系统各隔间温度、压力、气体成分以及钢制安全壳壁面温度、液膜的温度和厚度；
[0153]
步骤2：基于钢制安全壳内部热工水力状态判断空气冷却系统是否启动，具体包括如下内容：
[0154]
1)当钢制安全壳内部压力未达到启动限值，判定空气冷却系统不启动；
[0155]
2)当钢制安全壳内部压力达到启动限值，判定空气冷却系统启动，开始工作，并记录压力达到限值时间点tn，设定时间步长δτ；
[0156]
步骤3：利用多节点安全壳分析程序计算钢制安全壳内部隔间与钢制安全壳内壁面换热，具体包括如下内容：
[0157]
1)基于tn时刻钢制安全壳内部温度、压力和气体成分以及tn时刻钢制安全壳壁面温度计算δτ时间步长内两者之间的换热量q
in
；
[0158]
2)基于tn时刻钢制安全壳内的热工水力状态、质量源项和δτ时间步长内钢制安
全壳内部和钢制安全壳内壁面的换热量，计算该时间步长结束t
n 1
时刻钢制安全壳内部的温度、压力、气体成分热工水力状态；
[0159]
步骤4：图3中15-28号隔间为空气冷却系统隔间，计算某个空气冷却系统隔间与钢制安全壳外壁面换热及其热工水力状态，具体包括如下内容：
[0160]
1)计算空气冷却系统隔间与钢制安全壳外壁面换热
[0161]
i.基于tn时刻空气冷却系统隔间气体温度、钢制安全壳壁面温度和液膜厚度计算δτ时间步长内空气冷却系统隔间气体和钢制安全壳外壁面之间的辐射换热量
[0162][0163]
式中：
[0164]qrad
——δτ时间步长内空气冷却系统隔间气体和钢制安全壳外壁面之间的辐射换热量；
[0165]ffilm
——液膜修正系数，若钢制安全壳外壁面未被液膜覆盖，则f
film
＝1，否则f
film
《1；
[0166]
a——钢制安全壳壁面面积；
[0167]
σ——玻尔兹曼常数；
[0168]
tg——tn时刻空气冷却系统隔间气体温度；
[0169]
t
wall
——tn时刻钢制安全壳壁面温度；
[0170]eg
——空气冷却系统隔间气体发射率；
[0171]ewall
——钢制安全壳壁面发射率；
[0172]
ii.i如图4所示，部分钢制安全壳外壁面被液膜覆盖，空气冷却系统隔间与钢制安全壳外壁面的对流换热量q
conv
＝0；
[0173]
液膜自上而下流动时被钢制安全壳加热，伴随液膜蒸发现象，空气冷却系统隔间和液膜之间发生传热传质；
[0174]
基于tn时刻空气冷却系统隔间气体温度、气体成分和液膜温度及物性计算δτ时间步长内空气冷却系统隔间和液膜的传热传质：
[0175]
δτ时间步长内空气冷却系统隔间和液膜之间的传质即液膜蒸发量采用下式计算：
[0176][0177]
式中：
[0178]mevap
——δτ时间内液膜蒸发量；
[0179]hcd
——tn时刻液膜蒸发传热系数；
[0180]
t
film
——tn时刻液膜温度；
[0181]
δτ——时间步长；
[0182]
δh——汽化潜热；
[0183]
蒸发传热系数h
cd
采用uchida关系式进行计算：
[0184][0185]
式中：
[0186]
ρ
steam
——tn时刻空气冷系统隔间蒸汽密度；
[0187]
ρ
ng
——tn时刻空气冷系统隔间不凝性气体密度；
[0188]
δτ时间步长内液膜和空气冷却系统隔间的传热量采用下式计算：
[0189]qfilm_rb
＝m
evaphevap
ah
conv
(t
film-tg)δτ
[0190]
式中：
[0191]qfilm_rb
——δτ时间内液膜和空气冷却系统隔间传热量；
[0192]hevap
——tn时刻蒸汽比焓；
[0193]hconv
——液膜-隔间气体自然对流换热系数；
[0194]
因为液膜流速小，故研究液膜-隔间气体自然对流换热问题时将液膜视为固体，h
conv
采用如下方法计算：
[0195][0196]
nu——努塞尔数；
[0197]
l——钢制安全壳壁面特征尺寸；
[0198]
k——空气冷却系统隔间气体导热系数；
[0199]
努塞尔数nu的计算关系式与流动类型相关，流动类型通过雷诺数re2和格拉晓夫数gr的关系判断：
[0200][0201]
ρg——空气冷却系统隔间气体密度；
[0202]
v——空气冷却系统隔间气体流速；
[0203]
μ——空气冷却系统隔间气体动力粘度；
[0204][0205]
g——重力加速度；
[0206]
①
若re2＜1.0gr，为自然对流，努塞尔数nu采用以下关系式：
[0207]
ra＝gr
·
pr
[0208]
pr——空气冷却系统隔间气体普朗特数；
[0209]
ra——空气冷却系统隔间气体瑞利数；
[0210][0211]
ν——空气冷却系统隔间气体运动粘度；
[0212]
α——空气冷却系统隔间气体热扩散系数；
[0213]cp
——空气冷却系统隔间气体定压比热容；
[0214]
②
若re2＞10gr，为强迫对流，努塞尔数nu采用以下关系式：
[0215][0216]
nu
forced
＝0.023re
0.8
pr
0.3
ꢀꢀꢀꢀꢀꢀꢀꢀ
re＞6000
[0217]
lnnu
forced
＝a1lnre b1ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
2000≤re≤6000
[0218]dh
——空气冷却系统隔间水力学直径；
[0219]
③
若1.0gr≤re2≤10gr，为混合对流，努塞尔数nu采用以下关系式：
[0220]
nu
mixed
＝[(re2/gr-1)/9][nu
forced-nu
natural
] nu
natural
[0221]
ii.ii如图4所示，部分钢制安全壳外壁面未被液膜覆盖，则m
evap
＝0，q
film_rb
＝0；空气自下而上流动时，空气冷却系统隔间和钢制安全壳外壁面发生对流换热；
[0222]
δτ时间步长内空气冷却系统隔间与钢制安全壳外壁面的对流换热量基于空气冷却系统隔间气体温度、气体成分和钢制安全壳壁面温度计算：
[0223]qwall_rb
＝ah
conv
(t
wall-tg)δτ
[0224]
式中的对流换热系数h
conv
计算方法同步骤ii.i所述；
[0225]
2)基于tn时刻空气冷却系统隔间压力、温度、气体成分、δτ时间步长内空气冷却系统隔间与钢制安全壳外壁面的辐射换热量和对流换热量、空气冷却系统隔间与液膜之间的传热传质计算t
n 1
时刻空气冷却系统隔间压力、温度、气体成分；
[0226]
i.更新t
n 1
时刻空气冷却系统隔间蒸汽质量
[0227]mst_n 1
＝m
st
m
evap
[0228]
式中：
[0229]mst_n 1
——t
n 1
时刻空气冷却系统隔间蒸汽质量；
[0230]mst
——tn时刻空气冷却系统隔间蒸汽质量；
[0231]
ii.更新t
n 1
时刻空气冷却系统隔间温度
[0232][0233]
式中：
[0234]
t
g_n 1
——t
n 1
时刻空气冷却系统隔间气体温度；
[0235]cg
——tn时刻空气冷却系统隔间气体比热容；
[0236]
m——tn时刻空气冷却系统隔间气体质量；
[0237]
iii.更新t
n 1
时刻空气冷却系统隔间压力
[0238][0239]
p
n 1
——t
n 1
时刻气体压力；
[0240]
p
st_n 1
——t
n 1
时刻蒸汽分压，由蒸汽质量、气相体积和气体温度计算得到；
[0241]nmn_n 1
——t
n 1
时刻不凝气体摩尔数；
[0242]
rg——气体常数；
[0243]vg
——空气冷却系统隔间体积；
[0244]
步骤5：判断是否计算到最后一个隔间，若是，则到步骤7；否则，到步骤6；
[0245]
步骤6：计算钢制安全壳外壁面液膜传递过程，并对下一空气冷却系统隔间进行步骤3计算，具体包括如下内容：
[0246]
1)更新第i号钢制安全壳外壁面步长末t
n 1
时刻液膜质量
[0247]mfilm_n 1
＝m
film
m
in-m
out
m
sp-m
evap
[0248]mfilm_n 1
——t
n 1
时刻钢制安全壳外壁面液膜质量；
[0249]mfilm
——tn时刻钢制安全壳外壁面液膜质量；
[0250]min
——δτ时间步长内从i-1号钢制安全壳外壁面流入的液膜质量，若第i号钢制安全壳壁面上部无壁面，则m
in
＝0，否则m
in
为第i-1号壁面的m
out
值；
[0251]mout
——δτ时间内流入i 1号钢制安全壳外壁面的液膜质量，采用下式计算：
[0252][0253]
时刻钢制安全壳外壁面液膜平均流速，采用下式计算：
[0254][0255]
ρ
film
——tn时刻钢制安全壳外壁面液膜密度；
[0256]
μ
film
——tn时刻钢制安全壳外壁面液膜动力粘度；
[0257]
δ——tn时刻钢制安全壳外壁面液膜厚度，采用下式计算：
[0258][0259]
2)更新第i号钢制安全壳外壁面步长末t
n 1
时刻液膜温度
[0260]
i.基于tn时刻液膜温度和钢壳壁面温度计算δτ时间步长内液膜-钢制安全壳外壁面导热量：
[0261][0262]
式中：
[0263]qfilm-wall
——δτ时间步长内液膜和钢制安全壳外壁面导热量；
[0264]
δ
film
——tn时刻钢制安全壳外壁面液膜厚度；
[0265]
δ
wall
——钢制安全壳厚度；
[0266]kfilm
——tn时刻钢制安全壳外壁面液膜热导率；
[0267]kwall
——钢制安全壳热导率；
[0268]
ii.计算t
n 1
时刻液膜温度
[0269][0270]hfilm
——tn时刻钢制安全壳外壁面液膜比焓；
[0271]hfilm_n 1
——t
n 1
时刻钢制安全壳外壁面液膜比焓；
[0272]
已知更新后t
n 1
时刻的空气冷却系统隔间内压力和液膜比焓则利用插值法求出更新后钢制安全壳外壁面液膜温度；
[0273]
对下一隔间进行步骤4计算，隔间计算顺序为从上往下计算；对于图3中空气冷却系统隔间，计算顺序为由15号隔间、16号隔间依次计算至28号隔间；
[0274]
步骤7：计算钢制安全壳壁面温度，具体包括如下内容：
[0275]
1)计算得到各部分安全壳壁面的净吸热量：
[0276]
δq＝q
in-q
out
[0277]
式中：
[0278]
δq——壁面的净吸热量；
[0279]qin
——钢制安全壳内壁面从钢制安全壳内部隔间吸收的热量；
[0280]qout
——钢制安全壳外壁面向壳空气冷却系统散失的热量；
[0281]qout
＝q
film_wall
q
film_rb
[0282]
2)更新t
n 1
时刻各部分钢制安全壳壁面的温度：
[0283][0284]
t
wall_n 1
——t
n 1
时刻钢制安全壳壁面温度；
[0285]mwall
——钢制安全壳壁面质量；
[0286]cwall
——钢制安全壳壁面的定压比热容；
[0287]
步骤8：将t
n 1
时刻作为新的tn时刻，重复步骤3至7，直到达到指定计算时间。
[0288]
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施方式仅限于此，对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单的推演或替换，都应当视为属于本发明由所提交的权利要求书确定专利保护范围。