國(guó) 堯

（天津市眼科醫(yī)院 天津 300020）

1 引言

中央空調(diào)的優(yōu)化節(jié)能控制不僅對(duì)資源匱乏的我國(guó)能源戰(zhàn)略方面起到至關(guān)重要的作用，而且在使用過(guò)程中為室內(nèi)的溫度實(shí)時(shí)調(diào)節(jié)提供了新的研究領(lǐng)域［1～3］。中央空調(diào)令人們的生產(chǎn)生活變得更加舒適的同時(shí)，高耗能的壓縮機(jī)通過(guò)冷凍水系統(tǒng)對(duì)空氣進(jìn)行壓縮抽取，利用能量轉(zhuǎn)換實(shí)現(xiàn)了整個(gè)大環(huán)境內(nèi)的溫度的降低［4］。醫(yī)院作為大型的建筑物載體，中央空調(diào)的使用是必不可少的，而面對(duì)著不同房間因面積大小與環(huán)境的不同，如何提高中央空調(diào)的節(jié)能率，盡可能地減少電能消耗成為了目前研究領(lǐng)域的重點(diǎn)難題［5～6］。本研究針對(duì)醫(yī)院類型的中央空調(diào)系統(tǒng)節(jié)能優(yōu)化控制，設(shè)計(jì)了一套從房間內(nèi)的數(shù)據(jù)信息采集到溫控調(diào)節(jié)的優(yōu)化模型，并利用退火算法對(duì)設(shè)定參數(shù)預(yù)測(cè)模型的預(yù)測(cè)值與實(shí)際使用的實(shí)測(cè)值之間實(shí)時(shí)的計(jì)算偏差，將調(diào)控?cái)?shù)據(jù)以控制信號(hào)的形式通過(guò)總線傳輸?shù)紻DC控制器中，以此完成對(duì)中央空調(diào)各個(gè)設(shè)備的實(shí)時(shí)控制。

2 優(yōu)化建模

2.1 基本思路

針對(duì)醫(yī)院重要空調(diào)在實(shí)際部署中收到工況環(huán)境的條件約束，建立其優(yōu)化控制模型以減少中央空調(diào)電力消耗和管理人員工時(shí)操作。本研究以AI/DI傳感器到DDC控制器之間設(shè)計(jì)了一套面向醫(yī)院中央空調(diào)的優(yōu)化控制算法模型。具體的優(yōu)化控制建?；舅悸啡鐖D1所示。

圖1 中央空調(diào)優(yōu)化控制流程圖

由圖1所示，醫(yī)院中央空調(diào)優(yōu)化控制分為三個(gè)階段：采集過(guò)程、優(yōu)化過(guò)程和控制過(guò)程。其中，采集過(guò)程所用設(shè)置在醫(yī)院各大廳和具體房間內(nèi)部的溫控傳感器對(duì)控制度量進(jìn)行收集。通過(guò)總線傳輸?shù)侥茉垂芾砜刂葡到y(tǒng)（EMCS）中［7］，整合傳感器數(shù)據(jù)轉(zhuǎn)化為數(shù)據(jù)信號(hào)進(jìn)入系統(tǒng)管理層進(jìn)行識(shí)別，結(jié)合PID調(diào)節(jié)器對(duì)數(shù)據(jù)參數(shù)完成動(dòng)態(tài)調(diào)控［8］；優(yōu)化過(guò)程通過(guò)建立工藝參數(shù)（溫度、濕度和氣壓等）的目標(biāo)函數(shù)計(jì)算中央空調(diào)的能耗函數(shù)［9］，根據(jù)操作人員對(duì)控制變量組合運(yùn)用退火算法（SA）構(gòu)建約束條件計(jì)算最有控制過(guò)程，根據(jù)最優(yōu)設(shè)定值的操作指令通過(guò)系統(tǒng)傳輸反饋調(diào)節(jié)（PID調(diào)節(jié)器）輸出控制信號(hào)，在EMCS中依據(jù)總線傳輸分散到各個(gè)DDC控制器以調(diào)節(jié)各大廳和房間的中央空調(diào)運(yùn)轉(zhuǎn)。

2.2 數(shù)學(xué)建模

假設(shè)某醫(yī)院中央空調(diào)系統(tǒng)分為n個(gè)分散的空調(diào)設(shè)備，若對(duì)于第i個(gè)空調(diào)設(shè)備，其運(yùn)行過(guò)程中費(fèi)用開(kāi)銷為Ji，輸入與輸出變量分別為xi，yi，對(duì)于不可控、連續(xù)可控和離散可控變量設(shè)定為 fi，ui和Mi，則中央空調(diào)的費(fèi)用開(kāi)銷目標(biāo)函數(shù)為

其中，對(duì)于等式與不等式約束條件分別為|g(f，m，u)|=0 和 | h(f，m，u)|≥0。整個(gè)醫(yī)院的中央空調(diào)系統(tǒng)中的每一個(gè)設(shè)備作為分系統(tǒng)收到控制以減少電能消耗，優(yōu)化的基本原理如圖2所示。

圖2 優(yōu)化原理

本研究以構(gòu)建中央空調(diào)電能消耗為目標(biāo)函數(shù)，以運(yùn)行過(guò)程中盡可能地降低運(yùn)行成本為約束條件，綜合中央空調(diào)的能源分配方式，假設(shè)空調(diào)設(shè)備的能耗分為外掛機(jī)組P1、壓縮機(jī)組P2、制冷機(jī)組P3、冷卻機(jī)組P4、風(fēng)扇機(jī)組P5。因此，系統(tǒng)總能耗P優(yōu)化目標(biāo)函數(shù)為

2.3 優(yōu)化預(yù)測(cè)

中央空調(diào)設(shè)備節(jié)能的優(yōu)化控制過(guò)程依賴參數(shù)設(shè)定，則在優(yōu)化控制前需要預(yù)測(cè)判斷結(jié)果是否適應(yīng)條件函數(shù)的約束［10］。因此，我們首先根據(jù)中央空調(diào)使用的環(huán)境，設(shè)定參數(shù)估計(jì)模型并對(duì)各分系統(tǒng)的空調(diào)設(shè)備進(jìn)行濾定參數(shù)，再由系統(tǒng)建模結(jié)果對(duì)各設(shè)備進(jìn)行電能消耗的計(jì)算，最終采用本研究的優(yōu)化模型對(duì)醫(yī)院中央空調(diào)運(yùn)行系統(tǒng)進(jìn)行優(yōu)化控制。預(yù)測(cè)判斷流程如圖3所示。

圖3 預(yù)測(cè)判斷

由圖3的預(yù)測(cè)判斷流程，通過(guò)建立自適應(yīng)函數(shù)的方式對(duì)空調(diào)設(shè)備進(jìn)行系統(tǒng)識(shí)別，其具體過(guò)程分為兩個(gè)方面：

1）根據(jù)傳感器與控制器在每個(gè)時(shí)刻采集數(shù)據(jù)樣本自適應(yīng)得到模型參數(shù)在系統(tǒng)中參數(shù)濾定，為下一時(shí)刻的預(yù)測(cè)判斷做準(zhǔn)備；

2）根據(jù)AI/DI傳感器在當(dāng)前時(shí)刻探測(cè)到的各房間采集樣本測(cè)量值與預(yù)測(cè)判斷得出的預(yù)測(cè)值進(jìn)行偏差計(jì)算，并通過(guò)系統(tǒng)傳輸?shù)竭^(guò)濾器中以判斷下一時(shí)刻控制條件的校正參數(shù)。

將數(shù)據(jù)送入優(yōu)化器，進(jìn)行設(shè)備模型的系統(tǒng)辨識(shí)和預(yù)測(cè)校正，并綜合形成過(guò)程模型的輸入數(shù)據(jù)文件［11］。根據(jù)當(dāng)前的過(guò)程約束和控制策略，計(jì)算優(yōu)化操作工況。對(duì)優(yōu)化結(jié)果進(jìn)行合理性判斷，并檢查是否位于操作可行域內(nèi)。在該段時(shí)間內(nèi)，操作人員可以采納全部或部分指導(dǎo)信息，也可以完全拒絕。采納時(shí)則手動(dòng)修改控制回路設(shè)定值。優(yōu)化結(jié)果不經(jīng)過(guò)人的參與或干預(yù)，直接送至到下層作為控制回路的設(shè)定值，并將約束變量作為控制約束條件，操作人員監(jiān)視過(guò)程的變化。

3 退火算法（SA）

3.1 狀態(tài)函數(shù)

由各設(shè)備優(yōu)化控制過(guò)程中預(yù)測(cè)判斷的各個(gè)時(shí)刻產(chǎn)生狀態(tài)函數(shù)，考慮到中央空調(diào)在使用過(guò)程中隨著外界環(huán)境的不斷變化［12］，狀態(tài)函數(shù)的解空間是隨機(jī)產(chǎn)生［13］。候選解的計(jì)算結(jié)果由AI/DI傳感器探測(cè)到的采集樣本決定，通常當(dāng)前溫度條件緊隨上一時(shí)刻的結(jié)果而產(chǎn)生，以不同的接受概率條件給出參數(shù)濾定結(jié)果。具體遵循以下兩個(gè)原則：

1）當(dāng)醫(yī)院內(nèi)部溫度在目標(biāo)函數(shù)下降時(shí)作為可接受的概率條件時(shí)，目標(biāo)函數(shù)值將候選預(yù)測(cè)值增加以加大與現(xiàn)實(shí)系統(tǒng)測(cè)量值之間的誤差，通過(guò)過(guò)濾器識(shí)別參數(shù)為負(fù)值并由總線傳輸?shù)絻?yōu)化控制系統(tǒng)中，系統(tǒng)利用DDC控制器完成升溫的調(diào)節(jié)；

2）當(dāng)醫(yī)院內(nèi)部溫度在目標(biāo)函數(shù)上升時(shí)作為可接受的概率條件時(shí)，目標(biāo)函數(shù)值將候選預(yù)測(cè)值增加以加大與現(xiàn)實(shí)系統(tǒng)測(cè)量值之間的誤差，通過(guò)過(guò)濾器識(shí)別參數(shù)為正值并由總線傳輸?shù)絻?yōu)化控制系統(tǒng)中，系統(tǒng)利用DDC控制器完成降溫的調(diào)節(jié)；

考慮到狀態(tài)函數(shù)的引入是算法實(shí)現(xiàn)全局搜索的最關(guān)鍵因素［14］，因此試驗(yàn)證明狀態(tài)接受函數(shù)的具體形式對(duì)算法性能的影響不敏感，算法中通常采用作為狀態(tài)接受函數(shù)。

3.2 溫控調(diào)節(jié)

假設(shè)醫(yī)院室內(nèi)的初始溫度為t0，初溫越大，獲得高質(zhì)量解的幾率也越大，但花費(fèi)的計(jì)算時(shí)間將增加。因此，通過(guò)設(shè)計(jì)的參數(shù)以調(diào)節(jié)初溫改變優(yōu)化質(zhì)量和效率。具體的算法的溫控調(diào)節(jié)步驟如下：

步驟一：均勻抽樣并產(chǎn)生一組狀態(tài)，以各狀態(tài)目標(biāo)值的方差為初溫；

步驟二：隨機(jī)產(chǎn)生一組狀態(tài)，確定組內(nèi)兩兩狀態(tài)之間的最大目標(biāo)值 ||Δmax，并依據(jù)差值條件；

步驟三：在外循環(huán)中修改溫度值以更新函數(shù)方式，即下一時(shí)刻與上一時(shí)刻呈現(xiàn)線性的變化，tk+1=λtt，其中 λ∈(0，1)可以不斷調(diào)整變化。

3.3 算法準(zhǔn)則

通過(guò)AI/DI傳感器對(duì)初溫的探測(cè)傳入系統(tǒng)到計(jì)算預(yù)測(cè)模型預(yù)測(cè)值與現(xiàn)實(shí)系統(tǒng)實(shí)測(cè)值之間的誤差，退火算法對(duì)誤差值的計(jì)算存在內(nèi)外兩個(gè)循環(huán)，具體遵循的準(zhǔn)則分別如下：

1）內(nèi)循環(huán)計(jì)算目標(biāo)函數(shù)的均值是否已經(jīng)達(dá)到穩(wěn)定，若連續(xù)的目標(biāo)值變化范圍在可接受的概率條件下，按照隨時(shí)刻樣本抽樣的方式判斷中央空調(diào)最終溫度是否穩(wěn)定；

2）外循環(huán)終止準(zhǔn)則也就是算法終止準(zhǔn)則，用于決定算法何時(shí)結(jié)束。雖然算法的收斂性理論中要求趨于零，但實(shí)際實(shí)際中并不這樣執(zhí)行。首先設(shè)置終值溫度閥值和外循環(huán)的迭代次數(shù)，其次，算法搜索到的最優(yōu)值是否已經(jīng)連續(xù)若干步保持不變，最后檢驗(yàn)系統(tǒng)熵是否已經(jīng)穩(wěn)定。

3.4 算法步驟

針對(duì)本研究提出的優(yōu)化控制模擬退火算法的步驟描述如下：

步驟一：給定初溫t=t0，隨機(jī)產(chǎn)生初始狀態(tài)s=s0，令 k=0；

步驟二：Repeat：

Repeat：

產(chǎn)生新?tīng)顟B(tài)sj=Genete(s)；

if min{1，exp[(-C(sj)-

C(s))/tk]}≥random[0，1]s=sj；

until抽樣滿足算法穩(wěn)定準(zhǔn)則；

退溫tk+1=update(tk)并令k=k+1；

步驟三：until滿足算法終止準(zhǔn)則；

步驟四：輸出算法搜索結(jié)果。

4 實(shí)驗(yàn)?zāi)M

針對(duì)類似醫(yī)院這類大型的中央空調(diào)使用過(guò)程中的優(yōu)化控制，本研究利用模擬退火算法對(duì)所構(gòu)建的預(yù)測(cè)模型的預(yù)測(cè)值與現(xiàn)實(shí)系統(tǒng)的測(cè)量值之間構(gòu)建了實(shí)時(shí)調(diào)控的誤差項(xiàng)，用于傳輸?shù)竭^(guò)濾器對(duì)DDC控制器進(jìn)行動(dòng)態(tài)溫控的調(diào)整，進(jìn)而計(jì)算中央空調(diào)的費(fèi)用開(kāi)銷。

4.1 實(shí)驗(yàn)環(huán)境

測(cè)試地點(diǎn)為某醫(yī)院中央空調(diào)系統(tǒng)，測(cè)試時(shí)間是2016年8月15日，實(shí)驗(yàn)對(duì)象是中央空調(diào)系統(tǒng)中的外掛機(jī)、壓縮機(jī)、制冷機(jī)、冷卻機(jī)和風(fēng)扇機(jī)各1組，測(cè)試參數(shù)包含室內(nèi)溫度、室內(nèi)相對(duì)濕度、外掛機(jī)、壓縮機(jī)、制冷機(jī)、冷卻機(jī)和風(fēng)扇機(jī)能耗等參數(shù)。依據(jù)本研究所設(shè)計(jì)的中央空調(diào)優(yōu)化控制，AI/DI傳感器將采集到的信息通過(guò)數(shù)字信號(hào)傳輸?shù)较到y(tǒng)進(jìn)行識(shí)別。設(shè)定6h內(nèi)以12min為周期的傳感器現(xiàn)場(chǎng)采集操作，從初始時(shí)刻到結(jié)束時(shí)刻共對(duì)采集的參數(shù)信息31次。假設(shè)傳感器儀器與醫(yī)院環(huán)境干擾不足以令實(shí)驗(yàn)測(cè)試出現(xiàn)較大誤差項(xiàng)。

4.2 優(yōu)化控制

利用本研究設(shè)計(jì)的中央空調(diào)優(yōu)化系統(tǒng)，對(duì)預(yù)測(cè)值與測(cè)量值之間的偏差進(jìn)行退火算法的模擬計(jì)算，以外掛機(jī)設(shè)備在電能消耗過(guò)程中的溫度為空調(diào)設(shè)備的測(cè)試對(duì)象，系統(tǒng)優(yōu)化前后的采樣數(shù)據(jù)如表1所示。

表1 系統(tǒng)優(yōu)化前后采樣數(shù)據(jù)對(duì)比

由表1可見(jiàn)，中央空調(diào)系統(tǒng)在優(yōu)化前后的外掛機(jī)設(shè)備溫度差異較大，這主要是由于過(guò)多的電能消耗為熱量表現(xiàn)。運(yùn)用退火算法可對(duì)中央空調(diào)系統(tǒng)進(jìn)行實(shí)時(shí)的數(shù)據(jù)采集以及時(shí)刻的模擬計(jì)算調(diào)節(jié)系統(tǒng)的控制信號(hào)反饋給DDC控制器。在以整個(gè)中央空調(diào)系統(tǒng)的能耗與節(jié)約用電角度分析，選擇8：00～10：00h內(nèi)的11次優(yōu)化前后用電能耗進(jìn)行實(shí)測(cè)，具體的優(yōu)化前后能耗對(duì)比與節(jié)能率如圖4所示。

圖4 能耗對(duì)比與節(jié)能率

從圖4可見(jiàn)，（a）中給出了優(yōu)化后的中央空調(diào)各設(shè)備的總能耗明顯低于優(yōu)化前的情況，即采用退火算法對(duì)中央空調(diào)預(yù)測(cè)值與測(cè)量值之間搭建的偏差數(shù)值的實(shí)時(shí)計(jì)算可及時(shí)地調(diào)節(jié)中央空調(diào)系統(tǒng)的電能消耗，降低整個(gè)系統(tǒng)的能耗且提高綜合利用率。由（b）所示，整個(gè)中央空調(diào)系統(tǒng)的各個(gè)設(shè)備可節(jié)省能耗最大值為86kW，節(jié)能率為16.24%，2h內(nèi)的平均節(jié)能率為9.12%。若全天24h使用中央空調(diào)可節(jié)省電能1257.6kW/h，按照每度電價(jià)格0.56元計(jì)算，共可節(jié)省電費(fèi)704.26元。通過(guò)中央空調(diào)的整體優(yōu)化調(diào)控可有效地降低電能的浪費(fèi)，最終達(dá)到節(jié)能的目的。

5 結(jié)語(yǔ)

本文根據(jù)目前醫(yī)院中央空調(diào)系統(tǒng)在使用過(guò)程中的存在能量消耗過(guò)大的問(wèn)題，從AI/DI傳感器到DDC控制器之間設(shè)計(jì)了一套面向醫(yī)院中央空調(diào)的優(yōu)化控制系統(tǒng)，該系統(tǒng)中首先將采集到的參數(shù)信息以數(shù)字信號(hào)的形式傳輸?shù)较到y(tǒng)中已備識(shí)別，通過(guò)建立各空調(diào)設(shè)備費(fèi)用開(kāi)銷與總能耗的目標(biāo)函數(shù)，利用退火算法對(duì)模型的樣本測(cè)量值與預(yù)測(cè)值之間進(jìn)行實(shí)時(shí)的偏差計(jì)算以調(diào)節(jié)最優(yōu)設(shè)定值操作，解決了大型的中央空調(diào)設(shè)備能耗浪費(fèi)的問(wèn)題，實(shí)現(xiàn)了對(duì)建筑物室內(nèi)問(wèn)題的統(tǒng)一化管理。通過(guò)實(shí)驗(yàn)測(cè)試與分析結(jié)果表明，退火算法對(duì)采集信息的快速計(jì)算為中央空調(diào)系統(tǒng)的實(shí)時(shí)優(yōu)化調(diào)控提供了最優(yōu)參數(shù)保證。當(dāng)然該優(yōu)化控制還存在一些不足，在今后的工作中還將繼續(xù)研究改進(jìn)。最后，希望本研究對(duì)醫(yī)院中央空調(diào)的優(yōu)化控制工作提供一些借鑒。

［1］吳俊峰，張秀平，張朝暉.中央空調(diào)系統(tǒng)節(jié)能改造方案的節(jié)能量測(cè)量及驗(yàn)證方法初探［J］.制冷與空調(diào)，2013，13（9）：78-83.

WU Junfeng，ZHANG Xiuping，ZHANG Chaohui.Energy Saving Measurement and Verification Method of Energy Saving Reconstruction Scheme for Central Air-condition?ing System［J］.Refrigeration and Air Conditioning，2013，13（9）：78-83.

［2］戚野白，王丹，賈宏杰.基于局部終端溫度調(diào)節(jié)的中央空調(diào)需求響應(yīng)控制策略［J］.電力系統(tǒng)自動(dòng)化，2015（17）：82-88.

QI Yebai，WANG Dan，JIA Hongjie.Central air condition?ing demand response control strategy based on local termi?nal temperature regulation［J］.Automation of Electric Pow?er Systems，2015（17）：82-88.

［3］王蓓蓓，朱峰，嵇文路.中央空調(diào)降負(fù)荷潛力建模及影響因素分析［J］.電力系統(tǒng)自動(dòng)化，2016，40（19）：44-52.

WANG Peipei，ZHU Feng，JI Wenlu.Modeling of Central Air Conditioning Load Drop Potential and Analysis of In?fluencing Factors［J］.Power System Automatic.2016，40（19）：44-52.

［4］沈?qū)m新，李云明，劉其和.中央空調(diào)仿真培訓(xùn)系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)［J］.計(jì)算機(jī)測(cè)量與控制，2014，22（5）：1623-1625.

SHEN Gongxin，LI Yunming，LIU Qihe.Design and Real?ization of Simulation Training System for Central Air Con?ditioning［J］.Computer Measurement&Control.2014，22（5）：1623-1625.

［5］王文謙.論物聯(lián)網(wǎng)與云計(jì)算在中央空調(diào)行業(yè)方面的應(yīng)用［J］.電子技術(shù)與軟件工程，2015（13）：12-13.

WANG Wenqian.Application of Internet of Things and Cloud Computing in Central Air-conditioning Industry［J］.Electronic Technology and Software Engineering.2015（13）：12-13.

［6］張過(guò)有.基于物聯(lián)網(wǎng)技術(shù)的建筑電氣節(jié)能應(yīng)用研究［J］.通訊世界，2016（16）：20-21.

ZHANG Guoyou.Application of Building Electrical Ener?gy-saving Based on Internet of Things［J］.Communica?tion World，2016（16）：20-21.

［7］陶維成.面向上下文圖形可視化挖掘網(wǎng)絡(luò)行為與管理［J］.計(jì)算機(jī)與現(xiàn)代化，2012（3）：143-146.

TAO Weicheng.Context-oriented graph visualization min?ing network behavior and management［J］.Computer and Modernization，2012（3）：143-146.

［8］劉曉東，韋燚.基于PID調(diào)節(jié)器的Buck電路的系統(tǒng)設(shè)計(jì)與仿真［J］.新型工業(yè)化，2015（11）：27-33.

LIU Xiaodong，WEI Yi.System Design and Simulation of Buck Circuit Based on PID Regulator［J］.New Type of In?dustrialization.2015（11）：27-33.

［9］羅亮，吳文峻，張飛.面向云計(jì)算數(shù)據(jù)中心的能耗建模方法［J］.軟件學(xué)報(bào)，2014（7）：1371-1387.

LUO Liang，WU Wenjun，ZHANG Fei.Modeling Method for Energy Consumption in Cloud Computing Data Center［J］.Journal of Software，2014（7）：1371-1387.

［10］張法.網(wǎng)絡(luò)能耗系統(tǒng)模型及能效算法［J］.計(jì)算機(jī)學(xué)報(bào)，2012，35（3）：603-615.

ZHANG Fa.twork energy consumption system model and energy efficiency［J］.Chinese Journal of Computers，2012，35（3）：603-615.

［11］陳友明，張楠.基于能耗目標(biāo)函數(shù)的空調(diào)新風(fēng)量尋優(yōu)控制研究［J］.湖南大學(xué)學(xué)報(bào)（自然科學(xué)版），2012，39（2）：25-29.

CHEN Youming，ZHANG Nan.Study on Optimal Con?trol of Fresh Air Volume Based on Energy Consumption Objective Function［J］.Journal of Hunan University（Natural Sciences），2012，39（2）：25-29.

［12］朱美強(qiáng)，程玉虎，李明.一類基于譜方法的強(qiáng)化學(xué)習(xí)混合 遷 移 算 法［J］.自 動(dòng) 化 學(xué) 報(bào) ，2012，38（11）：1765-1776.

ZHU Meiqiang，CHENG Yuhu，LI Ming.A Hybrid Learn?ing Algorithm Based on Spectral Method［J］.Journal of Automation，2012，38（11）：1765-1776.

［13］賴雄鳴，張勇，王成.基于OELM重構(gòu)極限狀態(tài)函數(shù)的可靠度計(jì)算方法［J］.浙江大學(xué)學(xué)報(bào)（工學(xué)版），2015，49（4）：692-698.

LAI Xiongming，ZHANG Yong，WANG Cheng.Reliabil?ity calculation method based on oelm reconstruction limit state function［J］.Journal of Zhejiang University（Engi?neering Science），2015，49（4）：692-698.

［14］鄭毅，鄭蘋(píng).基于模糊熵和模擬退火算法的雙閾值圖像分割［J］.電子測(cè)量與儀器學(xué)報(bào)，2014，28（4）：360-367.

ZHENG Yi，ZHENG Ping.Double threshold image seg?mentation based on fuzzy entropy and simulated anneal?ing algorithm［J］.Journal of Electronic Measurement and Instrument，2014，28（4）：360-367.