专利名称：：一种检测热设备内部温度的方法
技术领域：
：本发明涉及一种监测热设备内部温度及其分布的方法。
背景技术：
：工业生产中的很多设备都与热过程有关，这些设备通常被称为热设备。热设备内部的温度及其分布包含设备运行状态的重要信息。由于在许多情况下，热设备的内部温度及其分布不便或难以直接测量，如何准确获得热设备内部温度信息一直是热工过程测试及状态监测
技术领域：
的重要课题之一。采用反演方法确定热设备内部温度，是解决该问题的一种较为有效的方法。该方法将热设备内部温度检测问题归结为一类传热学反问题，根据能够通过直接测量获得的相关温度(如热设备外表面等处的温度)，结合传热反问题方法来反演热设备的内部温度及其分布。其中，建立有效的传热学反问题求解方法，是采用反演方法检测热设备内部温度的关键。优化方法是目前求解各类传热学反问题的基本求解方法。其中，共轭梯度法(CGM)、Levenberg-Marquardt(L-M)法以及最速下降法(SDM)等已经在传热学反问题求解中得到了广泛应用。采用各类基于梯度的优化算法求解传热学反问题，容易陷入局部极值，其解对初始猜测值也有严重的依赖性。传热学反问题通常是Hadamard意义下的不适定(Ill-Posed)问题。这种不适定性主要表现在两个方面一方面，反问题中的输入信息往往是欠定的或者是超定的，进而导致反问题的解非唯一或不存在严格意义上的解；另一方面，传热学反问题的解对输入信息往往不具有连续依赖性，输入数据的观测误差极有可能在反演过程中被显著地放大，从而导致反演过程的不稳定。此外，从本质上讲，传热学反问题又是一类推理问题，即由传热系统的部分观测结果到系统定解条件的一种反向推理问题，推理的直接证据来源于对实际传热系统的观测。至少，从实际系统中得到的观测信息不可避免地包含着干扰噪声，并具有一定的观测误差和一定的不完整性，从而导致了推理证据的某种不确定性。所以，实际的传热学反问题是一类不确定性推理问题。目前所开展的传热学反问题研究方法均可归结为精确推理方法范畴，对于处理具有不适定性和不确定性的推理问题所固有的局限性十分明显。
发明内容本发明的目的是，提供一种具有良好抗不适定性、能够有效地克服内部温度反演结果对初始猜测值的严重依赖性、同时还能够有效地抑制内部温度反演过程中存在的误差放大现象的检测热设备内部温度的方法。实现所述发明目的技术方案是这样一种检测热设备内部温度的方法，该检测方法包括如下步骤；其中，与现有技术相同步骤有(一)对热设备内部温度检测问题进行初始化处理，通过初始化处理，设定反演过程迭代次数标识A:的初值为^:二0，根据热设备待检测温度区间和可检测温度区间的位置，确定传热问题的求解区域/2及其边界，该边界具体包括给定边界条件的第一边界r;和第二边界尸2、该求解区域的待检测边界/1，以及该求解区域的温度的可测量边界/;,建立该求解区域的稳态导热过程模型，给定待检测边界/l上的待检测温度《的初始猜测值7和可测量边界/;上各测点7'处的温度的实际测量值^/;其中，！'=7,2,w;_/=7,2,";"为可测量边界/;上设置的温度测点数目(二)根据已知的该待检测温度《的第A:次猜测值7气，构造待检测边界厂3的边界条件，采用数值解法求解传热学的正问题，获得可测量边界/;温度的计算结果f/7';本发明改进的步骤是(三)利用获得的可测量边界/;温度的计算结果t"，产生可测量边界/;上各测点/处的温度误差^，若该温度误差^满足迭代停止条件，则反演迭代过程结束，并以前述的待检测温度的《窠A:次猜测值7、作为热设备内部温度的检测结果；反之，则进入步骤(四)；所述迭代停止条件为te/^S，其中，S为预定的温度误差阈值，且5>0;(四)对于可测量边界/;上的各测点/，设置与温度误差^/对应的温度调整分量Zl^，通过选择模糊子集、模糊子集的隶属度函数A和模糊推理规则，对可测量边界/;上各测点7处的温度误差~和对应的温度调整分量分别进行模糊化处理,建立与可测量边界厂4上各测点7'对应的一维模糊推理单元组(F^，F/2，...，Fi")，经过模糊推理产生与可测量边界/;上各测点y对应的温度调整分量d^的值；(五)确定温度调整分量^^对于各待检测点Z'处温度补偿量d《的权值W,j，通过对该一维模糊推理单元组(，Fi2，...，输出的温度调整分量dWy进行加权综合，获得各待检测点/的温度补偿量d《，并利用该温度补偿量J《对待检测温度f的第A:次猜测值7~进行修正，产生待检测温度f的第^:+1次猜测值？^+~;(六)以产生的待检测温度6的第t+l次猜测值f;替代待检测温度6的第A次猜测值7~，并作为求解区域/3待检测边界/^的边界条件，同时对迭代次数标识A:进行更新，转入步骤(二)。从方案中可以看出，与现有的基于精确推理的热设备内部温度反演方法不同，本发明将热设备内部温度的检测问题归结为由热设备外部温度观测信息到内部温度信息的模糊推理问题，通过模糊反演来获得热设备内部温度及其分布。据此，本发明建立了一种基于模糊推理的传热学反问题求解方案，在此基础上提供了一种基于模糊反演的检测热设备内部温度的方法。由于模糊推理有对输入信息具有明显的抗干扰能力，推理过程具有良好的鲁棒性和容错能力；能够有效利用不精确、不确定和不完备信息进行推理和决策的特点，所以，与现有的基于精确推理的热设备内部温度反演方法相比较，本发明方法不但具有良好的抗不适定性，能够有效地克服内部温度反演结果对初始猜测值的严重依赖性，同时还能够有效地抑制内部温度反演过程中的存在的误差放大现象。下面结合附图对本发明作进一步的说明。图1为本发明方法的流程图图2为传热问题求解区域示意图图3为模糊子集的隶属度函数图图4为本发明方法在不同初始猜测值下的反演结果图5为L-M方法在不同初始猜测值下的反演结果图6为标准差0"=0.01时两种方法反演结果比较图7为标准差a=0.03时两种方法反演结果比较具体实施例方式一种检测热设备内部温度的方法，该方法包括如下步骤(参考图l):(一)对热设备内部温度检测问题进行初始化处理，通过初始化处理，设定反演过程迭代次数标识A:的初值为A:-O，根据热设备待检测温度区间和可检测温度区间的位置，确定传热问题的求解区域/2及其边界(参考图2)，该边界具体包括给定边界条件的第一边界/]和第二边界尸2、该求解区域的待检测边界^，以及该求解区域的温度的可测量边界/;，建立该求解区域的稳态导热过程模型，给定待检测边界/l上的待检测温度《的初始猜测值r和可测量边界厂4上各测点_/处的温度的实际测量值^/，其中，/=门入…，"；_/=/,2,"(下同)；w为可测量边界/;上设置的温度测点数目；(二)根据巳知的该待检测温度《的第A次猜测值7、'，构造待检测边界厂j的边界条件，采用数值解法求解传热学的正问题，获得可测量边界/;温度的计算结果^^;(三)利用获得的可测量边界/;温度的计算结果f/;'，产生可测量边界/;上各测点y处的温度误差^.，即^=一6，若该温度误差^满足迭代停止条件，则反演迭代过程结束，并以前述的待检测温度《的第A:次猜测值？^作为热设备内部温度的检测结果；反之，则进入步骤(四)；所述迭代停止条件为1>/SS式(1)式(1)中，S为预定的温度误差阈值，且￡>0;(四)对于可测量边界/;上的各测点,设置与温度误差^对应的温度调整分量zi.，通过选择模糊子集、模糊子集的隶属度函数A和模糊推理规则，对可测量边界〔上各测点处的温度误差~和对应的温度调整分量分别进行模糊化处理，建立与可测量边界厂4上各测点_/对应的一维模糊推理单元组(F^，Fi^，...，Fi)，经过模糊推理产生与可测量边界/;上各测点y对应的温度调整分量^",的值；(五)确定温度调整分量^w,对于各待检测点f处温度补偿量zl《的权值Wu，通过对该一维模糊推理单元组(Fi^，F72，...，输出的温度调整分量zl^.进行加权综合，获得各待检测点/的温度补偿量^《，并利用该温度补偿量」《对待检测温度《的第A:次猜测值7~进行修正，产生待检测温度《的第&+1次猜测值7"~;t;7(六)以产生的待检测温度《的第*+1次猜测值纟/替代待检测温度《的第A:次猜测值7~，并作为求解区域/2待检测边界r;的边界条件，同时对迭代次数标识A:进行更新，即取A:-A:+1，转入步骤(二)。进一步讲，本发明在步骤(一)中所说的根据热设备待检测温度区间和可检测温度区间的位置，确定传热问题的求解区域/2及其边界，需要保证待检测的温度点Z处于求解区域/3的待检测边界/;之上；同时，还要保证求解区域/2的可测量边界/;上的各测点y处的温度可以直接测量；本发明在步骤(一)中所说的求解区域"的稳态导热过程模型，包括求解区域/2内温度的控制方程和第一边界厂7和第二边界/^上的边界条件，该求解区域"内温度的控制方程为式(2)中的X和y为该二维系统的坐标，"JC，y)为求解区域^内点空间点(x，;;)处的温度。进一步讲，在本发明的(四)中所述的一维模糊推理单元组(尸^，尸及2，...，F/)，包括w个分别与可测量边界/;上各测点y对应的一维模糊推理单元Fi^;该一维模糊推理单元F/^.的输入为温度误差^，输出为温度调整分量zlw"将温度误差^和温度调整分量^",分成7级，即温度误差~和温度调整分量经过量化后的论域X为55少2式(2)X={-3,-2,-/，0,+7，+2,+3}式(3)温度误差^的模糊子集￡,和温度调整分量^",.的模糊子集U,分别取为:￡={肌層，肌za尸s,亂尸5}式(4)={肌層，肌尸S,肌尸5}式(5)其中，A为"负大"，iVM为"负中"，iVS为"负小"，ZO为"零"，尸S为"正小"，尸M为"正中"，尸5为"正大"；各模糊子集(a^,iVM,ivs,zaM,P仏/^}的隶属函数//选用三角形函数(参考图3)，一维模糊推理单元Fi,采用if-then规则进行模糊推理，对应的模糊推理规则见表1;表lF」R,的模糊推理规则表<table>tableseeoriginaldocumentpage10</column></row><table>然后，根据一维模糊推理单元F^产生模糊推理的结果，采用重心法确定与可测量边界/;上各测点y对应的温度调整分量d^。进一步讲，在本发明的步骤(五)中所述的温度调整分量^^对于各待检测点/处温度补偿量^《的权值W,j，按式(6)确定式(6)其中/=7,2，…，my=7,2,w(下同)；式(6)中的、.代表待检测边界/l上各待检测点/到可测量边界/;上各测点/之间的空间距离；在步骤(五)中所述的各待检测点；处的温度补偿量^1《，通过対温度调整分量d^进行加权综合获得，艮P:<formula>formulaseeoriginaldocumentpage10</formula>式(7)乂w根据各待检测点/处的温度补偿量」f，按式(8)对待检测温度《的第A:次猜测值？^'进行修正，产生待检测温度《的第A+1次猜测值^"^:<formula>formulaseeoriginaldocumentpage10</formula>下面给出一个采用本发明方法检测热设备内部温度的试验验证实例，并用对比的方法说明本发明的优越性。在验证实例中，传热问题的求解区域/2为一常物性、无内热源的矩形求解区域，该求解区域的在直角坐标系中的范围为0Sx《fl及0S^《6;其中，a=0.3m，6=0.12/。该求解区域的边界的具体位置分别为给定边界条件的第一边界/]:;c-O及0^y^6给定边界条件的第二边界厂2待检测边界/^:可测量边界/;:M""及y"在验证实例中，给定边界条件的第一边界A和可测量边界/;采用绝热边界条件，给定边界条件的第二边界厂2采用定热流边界条件。该求解区域"内温度的控制方程和边界条件为<formula>formulaseeoriginaldocumentpage11</formula>昧力=0,"65a:式(9)式(10)式(11)式(12)在以上诸式中，C为环境温度，A为材料的导热系数，"为边界上的对流传热系数；在本验证实例中，分别取^-20。C，；i=50『/m.oC，a=60^//2.。C。在本验证实例中，可测量边界/;上各测点_/处的温度的实际测量值^/采用仿真实验方法获得，即首先按照式(13)设置待检测边界厂3上均匀分布的11个待检测点/处的温度实际值z,<formula>formulaseeoriginaldocumentpage11</formula>式(13)其中，A为待检测边界尸j上待检测点/处的义坐标值；！'=1，2，…，11;_/=1，2，…，11(下同)。由上述温度实际值？^确定待检测边界/^的边界条件，求解区域/2内温度的控制方程，获得可测量边界/;上测点_/处的温度^，并以该温度^作为可测量边界/;上各测点/处的温度的实际测量值6，然后根据该实际测量值^反演待检测边界/^上的待检测温度《。在本验证实例中，通过本发明方法与常规的基于优化技术的内部温度反演方法的比较，说明本发明优越性；这里的基于优化技术的内部温度反演方法中，采用L-M方法求解对应的传热学反问题；在求解传热学反问题过程中，取预定的温度误差阈值￡=/0—J。首先比较采用本发明方法和L-M方法得到的热设备内部温度的反演结果对于初始猜测值的依赖性。分别取待检测温度的初始猜测值r,:2(TC和？^—200cC，采用本发明方法和L-M方法确定待检测边界7^上的待检测温度《。通过比较表明，本发明所提供的基于模糊推理技术的检测热设备内部温度的方法，对于待检测温度的初始猜测值的设置不敏感(参考图4)，反演结果具有较高的精度；基于L-M算法的检测热设备内部温度的方法则对于待检测温度的初始猜测值的设置十分敏感，检测结果严重依赖初始猜测值(参考图5)。进一步比较测量误差对于热设备内部温度的反演结果的影响。取待检测温度的初始猜测值7=200。C，当可测量边界分别存在标准差o"=0.01和标准差a=0.03的测量误差时，采用本发明方法和L-M方法确定待检测边界/^上的待检测温度f。当0"=0.01时，本发明方法的反演结果的最大误差为3.79%，L-M方法的反演结果的最大误差为7.86%(参考图6);当0"=0.03时，本发明方法的反演结果的最大误差为3.61%，L-M方法的反演结果的最大误差为10.89%(参考图7)。可见，随着标准差的增大，L-M方法的反演结果明显变差，本发明方法的反演结果具有良好的抗不适定性。权利要求1.一种检测热设备内部温度的方法，该方法包括对热设备内部温度检测问题进行初始化处理，通过初始化处理，设定反演过程迭代次数标识k的初值为k＝0，根据热设备待检测温度区间和可检测温度区间的位置，确定传热问题的求解区域(Ω)及其边界，该边界具体包括给定边界条件的第一边界(Γ1)和第二边界(Γ2)、该求解区域的待检测边界(Γ3)，以及该求解区域的温度的可测量边界(Γ4)，建立该求解区域的稳态导热过程模型，给定待检测边界(Γ3)上的待检测温度id="icf0001"file="A2009101909040002C1.tif"wi="9"he="6"top="83"left="84"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>的初始猜测值id="icf0002"file="A2009101909040002C2.tif"wi="11"he="4"top="84"left="124"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>和可测量边界(Γ4)上各测点j处的温度的实际测量值(tj)的步骤(一)；根据已知的待检测温度id="icf0003"file="A2009101909040002C3.tif"wi="9"he="6"top="94"left="158"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>的第k次猜测值id="icf0004"file="A2009101909040002C4.tif"wi="14"he="5"top="106"left="30"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>构造待检测边界(Γ3)的边界条件，采用数值解法求解传热学的正问题，获得可测量边界(Γ4)温度的计算结果(<overscore>t</overscore>ckj)的步骤(二)；其中j＝1，2，...，n；i＝1，2，...，n；n为可测量边界(Γ4)上设置的温度测点数目；其特征在于，在步骤(一)和步骤(二)的基础上，还包括如下步骤(三)利用获得的可测量边界(Γ4)温度的计算结果(<overscore>t</overscore>ckj)，产生可测量边界(Γ4)上各测点j处的温度误差(ej)，若该温度误差(ej)满足迭代停止条件，则反演迭代过程结束，并以前述的待检测温度id="icf0005"file="A2009101909040002C5.tif"wi="9"he="6"top="179"left="78"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>的第k次猜测值id="icf0006"file="A2009101909040002C6.tif"wi="12"he="5"top="179"left="121"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>作为热设备内部温度的检测结果；反之，则进入步骤(四)；所述迭代停止条件为id="icf0007"file="A2009101909040002C7.tif"wi="20"he="11"top="201"left="68"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>其中，ε为预定的温度误差阈值，且ε＞0；(四)对于可测量边界(Γ4)上的各测点j，设置与温度误差(ej)对应的温度调整分量(Δuj)，通过选择模糊子集、模糊子集的隶属度函数(μ)和模糊推理规则，对可测量边界(Γ4)上各测点j处的温度误差(ej)和对应的温度调整分量(Δuj)分别进行模糊化处理，建立与可测量边界(Γ4)上各测点j对应的一维模糊推理单元组(FR1，FR2，...，FRn)，经过模糊推理产生与可测量边界(Γ4)上各测点j对应的温度调整分量(Δuj)的值；(五)确定温度调整分量(Δuj)对于各待检测点i处温度补偿量id="icf0008"file="A2009101909040003C1.tif"wi="11"he="5"top="39"left="145"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>的权值(wij)，通过对该一维模糊推理单元组(FR1，FR2，...，FRn)输出的温度调整分量(Δuj)进行加权综合，获得各待检测点i的温度补偿量id="icf0009"file="A2009101909040003C2.tif"wi="12"he="5"top="59"left="100"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>并利用该温度补偿量id="icf0010"file="A2009101909040003C3.tif"wi="11"he="5"top="59"left="155"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>对待检测温度id="icf0011"file="A2009101909040003C4.tif"wi="8"he="6"top="69"left="32"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>的第k次猜测值id="icf0012"file="A2009101909040003C5.tif"wi="12"he="5"top="69"left="74"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>进行修正，产生待检测温度id="icf0013"file="A2009101909040003C6.tif"wi="8"he="6"top="69"left="139"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>的第k+1次猜测值(六)以产生的待检测温度id="icf0015"file="A2009101909040003C8.tif"wi="8"he="6"top="90"left="81"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>的第k+1次猜测值id="icf0016"file="A2009101909040003C9.tif"wi="15"he="5"top="90"left="128"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>替代待检测温度id="icf0017"file="A2009101909040003C10.tif"wi="8"he="6"top="89"left="176"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>的第k次猜测值id="icf0018"file="A2009101909040003C11.tif"wi="13"he="5"top="101"left="52"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>并作为求解区域(Ω)待检测边界(Γ3)的边界条件，同时对迭代次数标识k进行更新，转入步骤(二)。2.根据权利要求1所述检测热设备内部温度的方法，其特征在于，步骤(四)中所述的一维模糊推理单元组(f^，尸/2，...，/^)，包括"个分别与可测量边界(/;)上各测点/对应的一维模糊推理单元(F及p;该一维模糊推理单元(F/p的输入为温度误差(~)，输出为温度调整分量(」",)；将温度误差(~)和温度调整分量(J)分成7级，即温度误差(。)和温度调整分量(」"7.)经过量化后的论域(义)为<formula>formulaseeoriginaldocumentpage3</formula>温度误差(。)的模糊子集和温度调整分量(J~)的模糊子集(f/p分别取为<formula>formulaseeoriginaldocumentpage3</formula>其中，iV5为"负大"，iVM为"负中"，JVS为"负小"，Z6>为"零"，户5"为"正小"，尸Af为"正中"，尸5为"正大"；各模糊子集(JV5,M仏JVS,ZO,尸S,尸B)的隶属函数(A)选用三角形函数，一维模糊推理单元(F尺,)采用if-then规则进行模糊推理，对应的模糊推理规则见表l;<table>tableseeoriginaldocumentpage4</column></row><table>然后，根据一维模糊推理单元(Fip产生模糊推理的结果，采用重心法确定与可测量边界(厂4)上各测点j'对应的温度调整分量(zi"y)。3.根据权利要求1或2所述检测热设备内部温度的方法，其特征在于，在步骤(五)中所述的温度调整分量()对于各待检测点/处温度补偿量(^《)的权值(),按下式确定其中的C代表待检测边界(厂j)上各待检测点/到可测量边界(厂》上各测点y之间的空间距离；在步骤(五)中所述的各待检测点/处的温度补偿量(^《)，通过对温度调整分量()进行加权综合获得，即df=￡(W,;^"》；根据各待检测点/处的温度补偿量(J《)，按下式对待检测温度(f)的第&次猜测值t■u(fO进行修正，产生待检测温度(《)的第A+1次猜测值(^"、)严,=P,+Zl《。全文摘要一种检测热设备内部温度的方法。该方法是一种采用模糊推理技术反演热设备内部温度的方法，该方法的基本特征在于，将热设备内部温度的检测问题归结为由设备外部温度观测信息到内部温度信息的模糊推理问题，通过模糊反演获得热设备内部温度及其分布。与现有的基于反演方法的热设备内部温度检测方法相比，该方法具有良好的抗不适定性，能够有效地克服内部温度反演结果对初始猜测值的严重依赖性，同时还能够有效地抑制内部温度反演过程中的存在的误差放大现象。文档编号G01K3/00GK101660951SQ20091019090公开日2010年3月3日申请日期2009年9月21日优先权日2009年9月21日发明者朱丽娜,王广军,红陈申请人:重庆大学