一种基于特征识别的空调系统仿真模拟方法与流程

技术特征：

1.一种基于特征识别的空调系统仿真模拟方法，其特征在于，该方法包括以下步骤：

(1)根据既有空调系统的实测运行数据，采用最小二乘法分别求解得到以下几种模型的模型参数：空调系统中冷水机组性能预测模型、表冷器性能预测模型、冷却塔性能预测模型、水泵模型、流体输配管路阻力模型，将所求得的模型参数作为表征各部件结构特性的特征参数；

(2)按照所仿真模拟空调系统的实际连接关系，将所述空调系统中冷水机组性能预测模型、表冷器性能预测模型、冷却塔性能预测模型、水泵模型、流体输配管路阻力模型进行模拟连接；

(3)设定以下系统运行指标：冷水机组蒸发温度、冷水机组冷凝温度、表冷器进口冷冻水温和冷却塔进口冷却水温，将所述系统运行指标与模拟工况参数一并输入已在所述步骤(1)中确定了特征参数的各模型，计算得到空调系统各运行状态参数。

2.根据权利要求1所述的一种基于特征识别的空调系统仿真模拟方法，其特征在于，所述步骤(1)中，冷水机组性能预测模型包括蒸发器模型、冷凝器模型、压缩机模型和节流阀模型：

a.所述蒸发器模型在变水量工况下为：

$Q_{e1}=\frac{{\mathrm{\Δt}}_{e1}}{{A^{*}}_e{m_{w,e}}^{-x_1}+{B^{*}}_{1,e}{m_r}^{-y_1}+{C^{*}}_e}---\left(1\right)$

$Q_{e2}=\frac{{\mathrm{\Δt}}_{e2}}{{A^{*}}_e{m_{w,e}}^{-x_1}+{B^{*}}_{2,e}{m_r}^{-y_2}+{C^{*}}_e}---\left(2\right)$

Q_e＝m_w,ec_p,w(t_wi,e-t_wo,e)＝m_r(h_eo-h_ei) (3)

${\mathrm{\Δt}}_{e1}=\frac{(t_{w1,e}-t_{wo,e})}{ln(t_{w1,e}-t_e)-\ln (t_{wo,e}-t_e)}---\left(4\right)$

${\mathrm{\Δt}}_{e2}=\frac{(t_{wi,e}-t_{eo})-(t_{w1,e}-t_e)}{ln(t_{wi,e}-t_{eo})-ln(t_{w1,e}-t_e)}---\left(5\right)$

所述蒸发器模型在变水温工况下为：

$Q_{e1}=\frac{{\mathrm{\Δt}}_{e1}}{{B^{*}}_{1,e}{m_r}^{-y_1}+{C^{*}}_e}---\left(6\right)$

$Q_{e2}=\frac{{\mathrm{\Δt}}_{e2}}{{B^{*}}_{2,e}{m_r}^{-y_2}+{C^{*}}_e}---\left(7\right)$

Q_e＝m_w,ec_p,w(t_wi,e-t_wo,e)＝m_r(h_eo-h_ei) (8)

${\mathrm{\Δt}}_{e1}=\frac{(t_{w1,e}-t_{wo,e})}{ln(t_{w1,e}-t_e)-\ln (t_{wo,e}-t_e)}---\left(9\right)$

${\mathrm{\Δt}}_{e2}=\frac{(t_{wi,e}-t_{eo})-(t_{w1,e}-t_e)}{ln(t_{wi,e}-t_{eo})-ln(t_{w1,e}-t_e)}---\left(10\right)$

式中，Q_e1、Δt_e1分别为两相区的换热量和换热温差；Q_e2、Δt_e2分别为过热区的换热量和换热温差；m_w,e为冷冻水流量；c_p,w为水的比热；m_r为制冷剂流量；A^＊_e为与蒸发器结构和污垢热阻有关的模型参数；x₁为根据不同结构形式蒸发器的水侧换热系数经验公式得到的常系数；B^＊_1,e为与蒸发器结构和污垢热阻有关的模型参数；y₁为根据不同结构形式蒸发器的制冷剂侧换热系数经验公式得到的常系数；C^＊_e为与蒸发器结构和污垢热阻有关的模型参数；B^＊_2,e为与蒸发器结构和污垢热阻有关的模型参数；y₂为根据不同结构形式蒸发器的制冷剂侧换热系数经验公式得到的常系数；Q_e为蒸发器总换热量；t_wi,e，t_wo,e分别为冷冻水进口温度和出口温度；t_w1,e为两相区冷冻水入口温度；h_ei，h_eo分别为蒸发器入口焓值和出口的焓值；t_e为冷水机组蒸发温度；t_eo为压缩机吸气温度，即蒸发器出口的制冷剂温度；

b.所述冷凝器模型在变水量工况下为：

Q_c＝m_w,cc_p,w(t_wo,c-t_wi,c)＝m_r(h_ci-h_co) (11)

$\begin{array}{c}{Q}_c=\frac{{\mathrm{\Δt}}_{c1}}{{A^{*}}_c{m_{w,c}}^{-x_2}+{B^{*}}_{1,c}{m_r}^{-y_3}+{C^{*}}_c}+\frac{{\mathrm{\Δt}}_{c2}}{{A^{*}}_c{m_{w,c}}^{-x_2}+{B^{*}}_{2,c}{m_r}^{-y_4}+{C^{*}}_c}\\ +\frac{{\mathrm{\Δt}}_{c3}}{{A^{*}}_c{m_{w,c}}^{-x_2}+{B^{*}}_{3,c}{m_r}^{-y_5}+{C^{*}}_c}\end{array}---\left(12\right)$

${\mathrm{\Δt}}_{c1}=\frac{(t_c-t_{w1,c})-(t_{co}-t_{wi,c})}{ln\[(t_c-t_{w1,c})/(t_{co}-t_{wi,c})\]}---\left(13\right)$

${\mathrm{\Δt}}_{c2}=\frac{(t_{w1,c}-t_{w2,c})}{ln\[(t_c-t_{w2,c})/(t_c-t_{w1,c})\]}---\left(14\right)$

${\mathrm{\Δt}}_{c3}=\frac{(t_{ci}-t_{wo,c})-(t_c-t_{w2,c})}{ln\[(t_{ci}-t_{wo,c})/(t_c-t_{w2,c})\]}---\left(15\right)$

所述冷凝器模型在变水温工况下为：

Q_c＝m_w,cc_p,w(t_wo,c-t_wi,c)＝m_r(h_ci-h_co) (16)

$Q_c=\frac{{\mathrm{\Δt}}_{c1}}{{B^{*}}_{1,c}{m_r}^{-y_3}+{C^{*}}_c}+\frac{{\mathrm{\Δt}}_{c2}}{{B^{*}}_{2,c}{m_r}^{-y_4}+{C^{*}}_c}+\frac{{\mathrm{\Δt}}_{c3}}{{B^{*}}_{3,c}{m_r}^{-y_5}+{C^{*}}_c}---\left(17\right)$

${\mathrm{\Δt}}_{c1}=\frac{(t_c-t_{w1,c})-(t_{co}-t_{wi,c})}{ln\[(t_c-t_{w1,c})/(t_{co}-t_{wi,c})\]}---\left(18\right)$

${\mathrm{\Δt}}_{c2}=\frac{(t_{w1,c}-t_{w2,c})}{ln\[(t_c-t_{w2,c})/(t_c-t_{w1,c})\]}---\left(19\right)$

${\mathrm{\Δt}}_{c3}=\frac{(t_{ci}-t_{wo,c})-(t_c-t_{w2,c})}{ln\[(t_{ci}-t_{wo,c})/(t_c-t_{w2,c})\]}---\left(20\right)$

式中，Q_c为冷凝器换热量；m_w,c为冷却水流量；t_wi,c，t_wo,c分别为冷凝器的冷却水进口温度和出口温度；h_ci，h_co分别为冷凝器入口焓值和出口的焓值；Δt_c1，Δt_c2，Δt_c3分别为冷凝器过冷区、两相区和过热区的换热温差；A^＊_c为需要确定的冷凝器模型参数；x₂为根据不同结构形式冷凝器的水侧换热系数经验公式得到的常系数；B^＊_1,c为需要确定的冷凝器模型参数；y₃为根据不同结构形式冷凝器的制冷剂侧换热系数经验公式得到的常系数；C^＊_c为需要确定的冷凝器模型参数；B^＊_2,c为需要确定的冷凝器模型参数；y₄为根据不同结构形式冷凝器的制冷剂侧换热系数经验公式得到的常系数；B^＊_3,c为需要确定的冷凝器模型参数；y₅为根据不同结构形式冷凝器的制冷剂侧换热系数经验公式得到的常系数；t_c为冷水机组冷凝温度；t_co为制冷剂在冷凝器出口的温度；t_w1,c，t_w2,c分别为两相区冷却水进口温度和出口温度；t_ci为压缩机排气温度；

c.所述压缩机模型为：

$m_r=\mathrm{\λ}\frac{V_{th}}{v_1}---\left(21\right)$

V_th＝ψV_th0 (22)

$T_{ci}=T_{eo}{\left(\frac{p_c}{p_e}\right)}^{(k-1)/k}---\left(23\right)$

$P_{th}={\mathrm{\λV}}_{th}\frac{p_e\·k}{k-1}\[{\left(\frac{p_c}{p_e}\right)}^{(k-1)/k}-1\]---\left(24\right)$

$P_{in}=\frac{P_{th}}{{\mathrm{\η}}_e}---\left(25\right)$

式中，m_r为制冷剂的质量流量；λ为输气系数；V_th为压缩机理论输气量；v₁为压缩机吸气比容；ψ为冷水机组负荷率；V_th0为压缩机额定工况下理论输气量；T_eo为吸气温度；T_ci为排气温度；p_c为冷凝压力；p_e为蒸发压力；k为压缩过程多变指数；P_th为压缩机理论功率；P_in为压缩机实际功率；η_e为压缩机的电能效率；

d.节流阀模型为：

h_co＝h_ei (26)

t_eo＝t_e+△t_e (27)

t_co＝t_c-△t_c (28)

式中，h_co为节流阀入口焓值；h_ei为节流阀出口焓值；Δt_e为过热度；Δt_c为过冷度；

所述表冷器性能预测模型包括表冷器换热量模型、表冷器传热效能模型、表冷器接触系数模型、表冷器传热单元数模型、热容比模型、出风参数模型；

所述表冷器换热量模型为：

Q_b＝m_a,b(h_ai,b-h_ao,b)＝m_w,bc_p,w(t_wo,b-t_wi,b) (29)

所述表冷器传热效能模型为：

${\mathrm{\ϵ}}_{1,b}=\frac{t_{gi,b}-t_{go,b}}{t_{gi,b}-t_{wi,b}}=\frac{1-e^{-NTU(1-\mathrm{\γ})}}{1-{\mathrm{\γe}}^{-NTU(1-\mathrm{\γ})}}---\left(30\right)$

所述表冷器接触系数模型为：

${\mathrm{\ϵ}}_{2,b}=1-\frac{t_{go,b}-t_{so,b}}{t_{gi,b}-t_{si,b}}=1-e^{\frac{-{\mathrm{\α}}_{o,b}{F}_b}{m_{a,b}{c}_{p,a}}}---\left(31\right)$

所述表冷器传热单元数模型为：

$NTU=\frac{K_b{F}_b}{m_{a,b}{c}_{p,a}}---\left(32\right)$

所述热容比模型为：

$\mathrm{\γ}=\frac{m_{a,b}{c}_{p,a}}{m_{w,b}{c}_{p,w}}---\left(33\right)$

所述出风参数模型，在干工况时为：

t_go,b＝t_gi,b-ε_1,b(t_gi,b-t_wi,b) (34)

$t_{so,b}=t_{Li}+\frac{(t_{go,b}-t_{Li})(t_{si,b}-t_{Li})}{t_{gi,b}-t_{Li}}---\left(35\right)$

所述出风参数模型，在湿工况时为：

t_go,b＝t_gi,b-ε_2,b(t_gi,b-t_b) (36)

t_so,b＝t_go,b-(1-ε_2,b)(t_gi,b-t_si,b) (37)

所述出风参数模型，在临界工况时，采用上述干工况时、湿工况时的出风参数模型均可；

式中，Q_b为表冷器换热量；m_a,b为空气质量流量；h_ai,b为空气入口焓值；h_ao,b为空气出口焓值；m_w,b为水流量；t_wi,b，t_wo,b分别为表冷器进口水温和出口水温；ε_1,b为表冷器传热效能；t_gi，t_go分别为空气进口干球温度和出口干球温度；γ为热容比；NTU为表冷器传热单元数；ε_2,b为表冷器接触系数；t_so,b，t_si,b分别为空气出口湿球温度和进口湿球温度；a_o,b为空气侧换热系数；F_b为表冷器总换热面积；c_p,a为空气比热；K_b为表冷器总传热系数；t_Li为空气进口露点温度；

所述冷却塔性能预测模型，在冷却塔内空气在非饱和状态下时为：

$\frac{{\mathrm{dm}}_{w,t}}{dz}={\mathrm{\β}}_t{\mathrm{aF}}_z(X_{s,w}-X)---\left(38\right)$

$\frac{dX}{dz}=\frac{{\mathrm{\β}}_t{\mathrm{aF}}_z(X_{s,w}-X)}{m_{a,t}}---\left(39\right)$

$\begin{array}{c}\frac{{\mathrm{dT}}_{w,t}}{dz}=\frac{{\mathrm{\β}}_t{\mathrm{aF}}_z}{m_{w,t}{c}_{p,w}}\[Le(T_{w,t}-T_{a,t})(c_{p,a}+c_{p,v}X)+(r_0+c_{p,v}{T}_{w,t}-c_{p,w}{T}_{w,t})\\ \·(X_{s,w}-X)\end{array}---\left(40\right)$

$\begin{array}{c}\frac{{\mathrm{dT}}_{a,t}}{dz}=\frac{{\mathrm{\β}}_t{\mathrm{aF}}_z}{m_{a,t}(c_{p,a}+c_{p,v}X)}\[Le(T_{w,t}-T_{a,t})(c_{p,a}+c_{p,v}X)\\ +(c_{p,v}{T}_{w,t}-c_{p,v}{T}_{a,t})(X_{s,w}-X)\]\end{array}---\left(41\right)$

$Le={0.866}^{2/3}(\frac{X_{s,w}+0.622}{X+0.622}-1){\[ln\frac{X_{s,w}+0.622}{X+0.622}\]}^{-1}---\left(42\right)$

${\left({\mathrm{\β}}_t{\mathrm{aF}}_z\right)}_j=Me\frac{m_{wi,t}}{H_t}---\left(43\right)$

所述冷却塔性能预测模型，在冷却塔内空气在饱和状态下时为：

$\frac{{\mathrm{dm}}_{w,t}}{dz}={\mathrm{\β}}_t{\mathrm{aF}}_z(X_{s,w}-X_{s,a})---\left(44\right)$

$\frac{dX}{dz}=\frac{{\mathrm{\β}}_t{\mathrm{aF}}_z(X_{s,w}-X_{s,a})}{m_{a,t}}---\left(45\right)$

$Le={0.866}^{\frac{2}{3}}(\frac{X_{s,w}+0.622}{X_{s,a}+0.622}-1){\[ln\frac{X_{s,w}+0.622}{X_{s,a}+0.622}\]}^{-1}---\left(46\right)$

$\begin{array}{c}\frac{{\mathrm{dT}}_{w,t}}{dz}=\frac{{\mathrm{\β}}_t{\mathrm{aF}}_z}{m_{w,t}{c}_w}\[(r_0+c_{p,v}{T}_{w,t}-c_{p,w}{T}_{w,t})(X_{s,w}-X_{s,a})\\ +Le(T_{w,t}-T_{a,t})(c_{p,a}+c_{p,w}(X-X_{s,a})+Le(T_{w,t}-T_{a,t})\\ \·(c_{p,a}+c_{p,w}(X-X_{s,a})+c_{p,v}{X}_{s,a})\]\end{array}---\left(47\right)$

$\begin{array}{c}\frac{{\mathrm{dT}}_{a,t}}{dz}=-\frac{{\mathrm{\β}}_t{\mathrm{aA}}_z}{m_{a,t}}\[c_{p,a}Le(T_{a,t}-T_{w,t})-X_{s,w}(r_0+c_{p,v}{T}_{w,t})+c_w\left(Le\right(T_{a,t}-T_{w,t})\\ \·\left(X-X_{s,a}\right)+T_{a,t}\left(X_{s,w}-X_{s,a}\right))+X_{s,a}(r_0+c_{p,v}Le(T_{a,t}-T_{w,t})+c_{p,v}{T}_{w,t})\]\\ /\[c_{p,a}+c_{p,w}X+\frac{{\mathrm{dX}}_{s,a}}{{\mathrm{dT}}_{a,t}}(r_0+c_{p,v}{T}_{a,t}-c_{p,w}{T}_{a,t})+X_{s,a}(c_{p,v}-c_{p,w})\]\end{array}---\left(48\right)$

${\left({\mathrm{\β}}_t{\mathrm{aF}}_z\right)}_j=Me\frac{m_{wi,t}}{H_t}---\left(49\right)$

式中，m_w,t为冷却塔水质量流量；dz为冷却塔填料竖直方向微元长度；β_t为冷却塔传质系数；a为填料比表面积；F_z为冷却塔填料横截面积；X_s,w为水温对应的饱和空气含湿量；X为冷却塔内空气水蒸气质量组分；T_w,t为水温；Le为刘易斯数；c_p,v为水蒸气定压比热；r₀为水的汽化潜热；T_a,t为冷却塔内空气温度；m_a,t为冷却塔空气质量流量；X_s,_a为饱和空气含湿量；(β_taF_Z)_j为换热过程中的等效值，Me为模型参数，m_wi,t表示入口水流量，H_t表示填料高度；

所述水泵模型为：

H_p0＝a₀+a₁V₀+a₂V₀² (50)

P_in,p0＝b₀+b₁V₀+b₂V₀² (51)

$\frac{V_1}{V_0}=\frac{n_1}{n_0}---\left(52\right)$

$\frac{H_{p1}}{H_{p0}}={\left(\frac{n_1}{n_0}\right)}^2={\left(\frac{V_1}{V_0}\right)}^2---\left(53\right)$

$\frac{P_{in,p1}}{P_{in,p0}}={\left(\frac{n_1}{n_0}\right)}^3={\left(\frac{V_1}{V_0}\right)}^3---\left(54\right)$

式中，V₀、H_p0、P_in,p0分别为泵在转速n₀下的流量、扬程与功率，V₁、H_p1、P_in,p1分别为泵在转速n₁下的流量、扬程与功率；a₀、a₁、a₂、b₀、b₁、b₂为模型参数；

所述流体输配管路阻力模型为：

式中，H_f为设备或管路阻力，χ_i为第i段管路沿程阻力系数，l_i为第i段管路长度，d为管路水力直径，ζ_j为第j个管道局部阻力件阻力系数，V为管道体积流量；为模型参数。

3.根据权利要求1所述的一种基于特征识别的空调系统仿真模拟方法，其特征在于，所述步骤(1)中，既有空调系统的实测运行数据具体为：求解所述冷水机组性能预测模型的模型参数，即求解冷凝器的模型参数、蒸发器的模型参数和压缩机的模型参数时，需要的实测数据有压缩机吸气温度t_eo、压缩机吸气压力p_e、压缩机排气温度t_ci、压缩机排气压力p_c、冷水机组负荷率ψ、压缩机功率P_in、冷却水流量m_w,c、冷凝器进口水温t_wi,c、冷凝器出口水温t_wo,c、冷冻水流量m_w,e、蒸发器进口水温t_wi,e、蒸发器出口水温t_wo,e；

求解所述表冷器性能预测模型的模型参数时，需要的实测数据为：冷冻水流量m_w,e、表冷器进口水温t_wi,b、表冷器出口水温t_wo,b、表冷器进口干球温度t_gi,b、表冷器进口湿球温度t_si,b、表冷器出口干球温度t_go,b、表冷器出口湿球温度t_so,b；

求解所述冷却塔性能预测模型的模型参数时，需要的实测数据为：冷却水流量m_w,c、冷却塔进口水温t_wi,t、冷却塔出口水温t_wo,t、冷却塔进口干球温度t_gi,t、冷却塔进口湿球温度t_si,t、冷却塔出口干球温度t_go,t、冷却塔出口湿球温度t_so,t；

求解所述水泵模型的模型参数时，需要的实测数据为：水泵在转速n₀下的流量、扬程和功率；

求解所述流体输配管路阻力模型的模型参数，需要实测的数据为：冷冻水流量m_w,e、冷却水流量m_w,c、水泵扬程H₁。

4.根据权利要求1所述的一种基于特征识别的空调系统仿真模拟方法，其特征在于，所述步骤(1)中，采用最小二乘法分别求解空调系统中冷水机组性能预测模型、表冷器性能预测模型、冷却塔性能预测模型、水泵模型和流体输配管路阻力模型的特征参数时，使用克莱姆法则保证线性方程组解的存在性和唯一性，最终求得特征参数的唯一解。

5.根据权利要求1、2、3或4所述的一种基于特征识别的空调系统仿真模拟方法，其特征在于，所述步骤(2)中模拟连接的具体方法为：按照工作介质流动方向，将各模型入口端的工作介质状态，作为与该入口端连接的模型出口端的工作介质状态，所述工作介质包括制冷剂、冷冻水和冷却水；将冷冻水管路阻力和设备阻力之和作为冷冻水泵扬程，将冷却水管路阻力和设备阻力之和作为冷却水泵扬程。

6.根据权利要求1、2、3或4所述的一种基于特征识别的空调系统仿真模拟方法，其特征在于，所述步骤(3)中，模拟工况参数具体为：冷水机组模型输入参数包括冷水机组负荷率、冷冻水流量、冷却水流量；表冷器性能预测模型输入参数包括：冷冻水流量、表冷器进口风量、表冷器进口干球温度和湿球温度；冷却塔性能预测模型输入参数包括：冷却水流量、冷却塔进口风量、冷却塔进口干球温度和湿球温度。

7.根据权利要求6所述的一种基于特征识别的空调系统仿真模拟方法，其特征在于，所述步骤(3)的具体流程为：

1)设定表冷器进口冷冻水温，将其与冷冻水流量、表冷器进口风量、表冷器进口干湿球温度一同输入表冷器性能预测模型，对表冷器出口冷冻水温进行模拟计算，所述表冷器出口冷冻水温即为蒸发器进口冷冻水温；

设定冷却塔进口冷却水温，将其与冷却水流量、冷却塔进口风量、冷却塔进口干湿球温度一同输入冷却塔性能预测模型，对冷却塔出口冷却水温进行模拟计算，所述冷却塔出口冷却水温即为冷凝器进口冷却水温；

2)设定冷水机组蒸发温度和冷凝温度，将其与冷水机组负荷率一同输入压缩机模型，根据常规压缩多变过程计算得到制冷剂流量、排气温度和压缩机电功率；

3)将设定的冷水机组蒸发温度和冷凝温度，与冷冻水流量、所述步骤2)中压缩机模型计算得到的制冷剂流量、所述步骤1)中得到的蒸发器进口冷冻水温，一同输入蒸发器模型，得到蒸发器进口制冷剂焓值和出口制冷剂焓值，结合制冷剂流量计算得到蒸发器内制冷剂侧换热量，根据热量守恒方程，利用蒸发器内制冷剂侧换热量等于冷冻水侧换热量计算得到蒸发器出口冷冻水温；根据传热基本方程，利用蒸发器总传热量和换热系数计算得到蒸发器中制冷剂和冷冻水的传热温差，所述蒸发器总传热量等于蒸发器内制冷剂侧换热量；

将设定的冷水机组冷凝温度，与冷却水流量、所述步骤2)中压缩机模型计算得到的排气温度、制冷剂流量、所述步骤1)中得到的冷凝器进口冷却水温，一同输入冷凝器模型，得到冷凝器进口制冷剂焓值和出口制冷剂焓值，结合制冷剂流量计算得到冷凝器内制冷剂侧换热量，根据热量守恒方程，利用冷凝器内制冷剂侧换热量等于冷却水侧换热量计算得到冷凝器出口冷却水温，根据传热基本方程，利用冷凝器总传热量和换热系数计算得到冷凝器中制冷剂和冷却水的传热温差，所述冷凝器总传热量等于冷凝器内制冷剂侧换热量；

4)根据所述步骤1)中得到的蒸发器进口冷冻水温、步骤3)中计算得到的蒸发器出口冷冻水温、蒸发器中制冷剂和冷冻水的传热温差，计算得到冷水机组蒸发器的蒸发温度；

根据所述步骤1)中得到的冷凝器进口冷却水温、步骤3)中计算得到的冷凝器出口冷却水温、冷凝器中制冷剂和冷却水的传热温差，计算得到冷水机组冷凝器的冷凝温度；

5)当以下四组条件均成立时，则将蒸发温度、冷凝温度、表冷器进口冷冻水温、表冷器出口冷冻水温、冷凝器进口冷却水温、冷凝器出口冷却水温、压缩机电功率、蒸发器总传热量和冷凝器总传热量作为系统模拟运行参数输出，否则用每组条件中的计算值更新与之比较的设定值，返回步骤1)：

所述步骤1)中设定的表冷器进口冷冻水温与所述步骤3)中计算得到的蒸发器出口冷冻水温相等，所述步骤1)中设定的冷却塔进口冷却水温与步骤3)中计算得到的冷凝器出口冷却水温相等，所述步骤3)中设定的冷水机组蒸发温度与步骤4)中计算得到的蒸发温度相等，所述步骤3)设定的冷水机组冷凝温度与步骤4)中计算得到的冷凝温度相等。

8.根据权利要求1、2、3或4所述的一种基于特征识别的空调系统仿真模拟方法，其特征在于，所述冷水机组性能预测模型、表冷器性能预测模型、冷却塔性能预测模型、水泵模型和流体输配管路阻力模型均是基于集总参数法，将既有空调系统中各部件的未知结构参数进行集总建立得到。

9.根据权利要求1、2、3或4所述的一种基于特征识别的空调系统仿真模拟方法，其特征在于，所述空调系统为单台冷水机组中央空调系统或由多台冷水机组、多台水泵和多套表冷器末端/混合末端所构成的空调系统。