(中核集團三門(mén)核電有限公司，浙江 三門(mén) 317112)
魏 來(lái)，陳 森
                                  
    魏來(lái)(1981-)男，安徽省淮南市人，工程師，華中科技大學(xué)工學(xué)學(xué)士，研究方向為全范圍模擬機開(kāi)發(fā)及維護。

    摘要：三門(mén)核電站采用AP1000全范圍模擬機，所有電站模型按照系統類(lèi)型和建模工具的不同可以劃分為堆芯系統、一回路熱工水力系統、氣液兩相流體系統、氣體或者液體單相流體系統、儀控系統以及電氣系統。本文主要介紹AP1000電站流體系統的建模原理和仿真過(guò)程，以AP1000乏燃料池冷卻系統為例，使用GSE公司的圖形化建模工具JTOPMERET對建模和仿真過(guò)程進(jìn)行說(shuō)明。

    關(guān)鍵詞：建模；仿真；流體系統

    Abstract: The AP1000 plant adopts full scope real time simulator, all the plant models divide into reactor core models , primary hydraulic fluid system, gas and liquid two phase fluid system, gas or liquid single phase fluid system,instrumentation and control system and electric system. This article mainly introduces the philosophy of the AP1000 plant’s fluid system modeling and process of the simulation, takes the spend fuel cooling system as an example,using the GSE company’s graphic modeling tools Jtopmeret to explain the process of the modeling and simulation.

    Key words: Modeling; Simulation; Fluid

    1 流體系統建模的主要原理以及組成

    1.1 電站模型分類(lèi)以及JTOPMERET簡(jiǎn)介

    AP1000電站全范圍模擬機由美國WEC公司總包，美國GSE公司分包，CNPO作為GSE的子承包商負責一部分核島、常規島以及公用系統建模。GSE公司和CNPO公司負責AP1000電站除儀控系統外的全部系統模型的仿真模型建立，AP1000電站模型根據系統以及建模工具的不同可以劃分為堆芯模型、一回路熱工水力模型、兩相流系統模型、電氣系統模型以及儀控系統模型，其中流體系統建模使用JTOPMERET圖形化建模軟件。

    JTOPMERET軟件是GSE公司于1990年開(kāi)始開(kāi)發(fā)，最初的著(zhù)眼點(diǎn)是解決兩相流的壓力矩陣，隨后又逐漸完善了熱傳遞方程，并增加了傳感器。1993年開(kāi)始圖形化編輯界面的開(kāi)發(fā)，使用核心方程計算反應率，放射性衰變熱等。經(jīng)過(guò)近十幾年的發(fā)展，JTOPMERET已經(jīng)逐漸成為流體系統建模的首要解決方案，其廣泛應用于核電、火電、石化等大型行業(yè)，以運行穩定、仿真精度高、界面友好易掌握獲得大量客戶(hù)的好評。

    1.2 流體系統建模遵循的主要方程以及計算的主要變量

    流體系統建模遵循動(dòng)量平衡、質(zhì)量平衡以及能量守恒三大定律，對于氣相以及液相均獨立計算動(dòng)量、質(zhì)量以及能量方程。

    動(dòng)量方程

        (1)

    質(zhì)量方程

      (2)

    能量方程

      (3)

    流體系統主要計算壓力、流量、焓、溫度、熱傳遞、氣相濃度、沸騰和凝結率、濃度、反應性、電導率、可溶化學(xué)濃度等重要參數。

    1.3 JTOPMERET的程序結構

    首先了解以下幾個(gè)概念，有助了解JTOPMERET的程序結構，如圖1所示，MST以及SST處于程序結構的最頂端，驅動(dòng)RTEXEC以及IEXEC，均處于服務(wù)器層。工程師層有控制模塊，控制模塊控制各程序段，各程序段控制模型內各組件，組件又由某些特定編碼的子程序構成：

    (1)MST(Master Synchronization Task)主計算機同步任務(wù)管理，控制主計算機實(shí)時(shí)系統執行；

    (2)SST(Slave Synchronization Task)從計算機同步任務(wù)管理，控制從計算機實(shí)時(shí)任務(wù)執行；

    (3)RTEXEC(Real Time Executive)實(shí)時(shí)執行，由MST或者SST驅動(dòng)，用于模擬機模型間的集成執行；

    (4)IEXEC(Interactive Executive)內部執行，用于實(shí)時(shí)多模型共同測試；

    (5)Control Modual控制模塊，其下層為程序段(segment)；

    (6)Segment 程序段，其下層為組件(component)，如泵、閥等；

    (7)Subroutine 子程序，通常為具有某些特定功能的子程序。
                               
                                       圖1  JTOPMERET程序結構 
    1.4 JTOPMERET流體網(wǎng)絡(luò )的主要組成

    JTOPMERET流體網(wǎng)絡(luò )主要由結點(diǎn)(n o d a l)、連接(link)、壓力邊界(pressure boundary)、流體邊界(flow boundary)、熱邊界(heat boundary)、泵、風(fēng)機、閥門(mén)、熱交換器等組成，下面簡(jiǎn)要介紹節點(diǎn)以及流體通道的類(lèi)型以及計算方程。

    1.4.1 結點(diǎn)

    JTOPMERET共有9種結點(diǎn)，依次為常規型、汽輪機型、容器型、冷卻塔型、熱交換器殼側型、熱交換器管側型、管道型、氣體分離器、液體分離器型，常規型結點(diǎn)約占全部結點(diǎn)使用率的90%，一般型結點(diǎn)主要用于將系統管道流程圖(P&ID)中的所有管道按照走向、高度以及尺寸劃分為流體網(wǎng)絡(luò )，將管道以及箱體轉化為結點(diǎn)。其中氣體分離器型結點(diǎn)可以精確的實(shí)現對汽輪機除濕再熱器的模擬。

    1.4.2 結點(diǎn)的主要計算方程

    流導(admittance)：
      (4)  
    

      

    水力直徑(Hydraulic Diameter)：

    A代筆面積，Pw代表周長(cháng)

    hd=4A/pw

    正方形區域面積周長(cháng)計算：

    
   
  
          
            

    三角型區域面積周長(cháng)計算：

      (5)  

      (6)  
         
            

    1.4.3 Link

    JTOPMERET共有3種link，依次為質(zhì)量流量連接，熱流量連接以及簡(jiǎn)單連接三種。其中質(zhì)量流量連接用于連接結點(diǎn)，熱流量連接用于連接熱邊界，簡(jiǎn)單連接用于連接各種測量?jì)x表，如壓力變送器，溫度變送器等。連接是沒(méi)有任何體積和質(zhì)量的單元，其中質(zhì)量流量連接用于連接兩個(gè)結點(diǎn)，閥門(mén)放置在連接之上。

    1.4.4 邊界

    邊界(boundary)有四種，分別對應壓力邊界、流量邊界以及熱流量邊界以及熱結點(diǎn)邊界。流量邊界與壓力邊界總是成對出現，通常在相互接口的兩個(gè)系統中一個(gè)系統為壓力邊界，則相對應的一個(gè)系統就為流量邊界。熱邊界代表熱源，如乏燃料池中乏燃料釋放的衰變熱就可以用熱邊界來(lái)模擬，只需輸入不同工況下的對應熱量釋放即可。敞口容器需要加上壓力邊界，用于連接系統所在區域環(huán)境，如乏燃料池處于核輔助廠(chǎng)房VAS范圍，那么就使用該壓力邊界連接SFS系統與VAS系統，與之類(lèi)似換料腔處于安全殼CNS系統內。

    2 乏燃料池冷卻系統建模
                     
                                  圖2  乏燃料池冷卻系統過(guò)程流程圖
                
                                圖3  乏燃料池冷卻系統JTOPMERET 網(wǎng)絡(luò )結點(diǎn)圖
  2.1 JTOPMERET網(wǎng)絡(luò )圖繪制一般原則

    JTOPMERET 網(wǎng)絡(luò )圖的繪制遵循從上到下，從左到右的原則。建立模型前首先要仔細閱讀系統建模要求說(shuō)明書(shū)，模擬范圍必需與系統劃定模擬范圍的P&ID圖保持一致，節點(diǎn)和部件擺放相對位置考慮P&ID圖以及過(guò)程流程圖參照真實(shí)位置高度，節點(diǎn)的編號以及連接的次序遵循與工藝流程一致的原則，泵、熱交換器以及過(guò)濾器的選擇以電廠(chǎng)設計參數為參考進(jìn)行選擇。邊界的選擇遵循體積大的，流量對其影響不大的結點(diǎn)選擇流量邊界，對應的其他系統的邊界選擇壓力邊界。

    一張JTOPMERET結點(diǎn)圖規定的最大結點(diǎn)數量為200，結點(diǎn)數量超過(guò)200必需考慮分為兩張圖。JTOPMERET規定了泵、閥門(mén)以及熱交換器等部件的命名規則。以泵為例，命名規則為‘系統簡(jiǎn)稱(chēng)+部件代號+編號’，以SFS系統01A泵為例其編號為SFSMP01A。

    2.2 乏燃料池冷卻系統中非標準部件的建模實(shí)現

    JTOPMERET有一個(gè)標準組件庫，一般的泵、風(fēng)機、閥門(mén)、熱交換器等部件均有，但是也并不能完全涵蓋AP1000電廠(chǎng)的所有設備，比如爆破閥、過(guò)濾器等。這些非標準設備就需要重新建立部件庫，重新編碼生成，也可使用JTOPMERET自帶的部件庫中的部件進(jìn)行組合來(lái)生成新的滿(mǎn)足功能的部件。以過(guò)濾器為例，可以將其根據體積以及高度差分為上下兩個(gè)結點(diǎn)，壓降平均分配，中間使用連接連接，在連接上放置閥門(mén)，以閥門(mén)的關(guān)閉來(lái)模擬過(guò)濾器的堵塞故障。

    2.3 故障的模擬

    故障是模型必需重點(diǎn)考慮的組成部分，是操縱員培訓必不可少的情景培訓。故障的種類(lèi)和數量要真實(shí)反映電廠(chǎng)設計和實(shí)際運行情況，常見(jiàn)的流體系統故障有管道小破口、大破口、泵堵轉卡殼、泵失去電源、泵超速、過(guò)濾器堵塞、閥門(mén)泄露誤動(dòng)作等通用故障。在乏燃料池冷卻系統中，根據系統模型設計說(shuō)明書(shū)，主要考慮設置了泵出口管道泄漏、熱交換器入口管道泄漏、除鹽器堵塞過(guò)濾器堵塞這三個(gè)類(lèi)型的故障。

    2.4 乏燃料池冷卻系統結點(diǎn)圖

    如前2.1所述，按照乏燃料池冷卻系統P&ID以及過(guò)程流程圖按照JTOPMERET繪制規則完成乏燃料池冷卻系統結點(diǎn)圖的繪制。

    2.5 模型輸入數據的準備

    模型輸入數據主要包括結點(diǎn)數據輸入、link數據輸入、邊界數據輸入、泵、風(fēng)機、熱交換器等部件的數據輸入，下面將介紹結點(diǎn)、link以及泵的數據輸入。

    2.5.1 結點(diǎn)的輸入數據

    主要包括結點(diǎn)類(lèi)型、雙相流或者單相流、是否過(guò)濾器選項、有無(wú)化學(xué)反應、結點(diǎn)是否可壓縮、水力直徑、截面積、是否線(xiàn)性體積、高度、體積、表面積、壁厚材料。
                 
                                  圖4  結點(diǎn)的輸入參數
    輸入數據由電廠(chǎng)設計數據直接計算而來(lái)，需要系統的管道軸側圖來(lái)計算高度，水力直徑、管道截面積、體積以及表面積，設計數據直接決定模型仿真的精度。

    2.5.2 Link的輸入參數

    Link存在有閥門(mén)與沒(méi)有閥門(mén)兩種情況，主要用于計算流導，詳細計算過(guò)程參見(jiàn)1.2.2.輸入參數主要有閥門(mén)類(lèi)型、上游壓力、下游壓力、流體速度等。帶節流孔板的flow path 選擇choke flow。沒(méi)有閥門(mén)的link閥門(mén)類(lèi)型選擇固定(fix)類(lèi)型，有閥門(mén)的link閥門(mén)類(lèi)型選擇線(xiàn)性(linear)。連接參數輸入的重點(diǎn)是將流體系統工藝流程的壓降進(jìn)行分配，這將直接決定流導的大小。

    2.5.3 泵的參數以及運行曲線(xiàn)
                  
                                    圖5  泵的輸入參數
    泵需要輸入凈吸入壓頭(NPSH)，運行曲線(xiàn)以及壓降，選取泵運行曲線(xiàn)上三個(gè)工作點(diǎn)，分別填入其壓頭以及對應的流量，所示完成泵的參數輸入。

    2.5.4 模型數據格式

    如前所述結點(diǎn)、連接、泵、閥門(mén)、熱交換器等參數輸入，實(shí)際是對模擬機模型計算方程和矩陣的參數輸入，圖形化的編程界面實(shí)際是提供了一個(gè)友好的標準參數輸入格式，便于生成標準矩陣代碼。所有的變量以及常量全部在JTOPMERET數據庫管理軟件DBM(database management system)中分配地址，在完成所有的參數輸入后就可以生成模型運行數據并鏈接可執行程序進(jìn)行調試，下面將以乏燃料池冷卻系統(sfs)為例介紹JTOPMERET生成的主要模型數據格式，子模塊在此就不作介紹。

    再完成源程序編譯后，我們得到一些主要的數據輸入塊以及控制模塊，JTOPMERET源代碼主要由C語(yǔ)言以及Fortran語(yǔ)言編譯。

    bk.sfs1 常量數據輸入塊

    bv.sfs1 變量常量數據輸入塊

    d.globla_k 常量數據庫文件

    d.global_v 變量數據庫文件

    sfs1.dbmc 常量數據庫輸入文件

    sfs1.dbmv 變量數據庫輸入文件

    sfs1.dbmb 邊界數據庫輸入文件

    sfs1.data 生成的模型輸入數據文件

    Ssfsd01t 外部參數輸入模塊

    Ssfsd02t 模型參數輸出模塊

    Sfsm16t1 控制模塊

    Sfsd11t 流導計算方程

    Sfsd12t 矩陣計算方程

    Sfsd13t 網(wǎng)絡(luò )流量計算方程

    Sfsd14t 焓、溫度、質(zhì)量、密度方程

    Sfsd15t 結點(diǎn)物理特性方程

    Sfsd18t 傳感器過(guò)程輸入方程

    3 單系統測試

    3.1 單系統調試步驟

    單系統測試是對單個(gè)系統進(jìn)行測試，沒(méi)有考慮與其他系統的接口。下面以乏燃料池冷卻系統為例，介紹模型的編譯調試步驟。

    (1)生成數據文件 SFS1.data；

    (2)生成源代碼；

    (3)在data base management 中加載常量(constant)，變量(variable)以及邊界(boundary)；

    (4)加載控制模塊以及子程序；

    (5)編譯所有控制模塊以及子程序；

    (6)生成模型可執行文件C.TESTSFS1文件；

    (7)鏈接可執行文件并啟動(dòng)S3serv；

    (8)啟動(dòng)調試窗口并加載生成的邊界文件；

    (9)運行模型并實(shí)現與JTOPMERET的通訊；

    (10)設置正確的部件開(kāi)啟次序(閥門(mén)、泵等)，對初始狀態(tài)進(jìn)行快照(snap IC)；

    (11)重啟初始條件(reset IC)；

    (12)驗證流體模型不同工況下的溫度、流量、壓力、液位等參數。

    3.2 單系統測試原則

    單系統測試的目的主要是驗證流體模型不同工況的溫度、流量、壓力、液位等重要參數以及初始條件和邊界的設置，做好記錄與電廠(chǎng)設計數據對比，如果發(fā)現模型數據與電廠(chǎng)設計數據不一致則需要重新檢查模型輸入數據以及初始條件設置，直到模型數據與電廠(chǎng)設計保持一致為止。以乏燃料池冷卻系統為例，就需要測試正常運行工況、事故工況下的系統功能，驗證正常運行工況下乏燃料池的冷卻凈化，換料腔的冷卻凈化以及反應堆冷卻劑轉移等功能，驗證溫度、壓力、流量、液位值、熱量等關(guān)鍵參數實(shí)驗值是否與設計參數保持一致。驗證事故工況下如乏燃料池沸騰，乏燃料池的補充裝量等功能。

    單系統測試作為集成測試(聯(lián)合調試)前的重要步驟具有重要意義，聯(lián)合調試的時(shí)間以及成敗很大程度上取決于單系統調試的質(zhì)量，在單系統調試中發(fā)現越多的錯誤越有助于節省聯(lián)合調試的時(shí)間，保證聯(lián)合調試的成功。

    3.3 乏燃料池冷卻系統的測試結果
                     
                                   圖6  乏燃料池冷卻系統單系統測試
    圖6是乏燃料池冷卻系統的單系統調試圖，目前正常運行工況以及事故工況下的各回路管線(xiàn)的流量、液位、水池溫度、結點(diǎn)壓力等參數與設計參數之間誤差在容許范圍內，單系統測試成功。由于目前AP1000 核島系統模型開(kāi)發(fā)還未進(jìn)入聯(lián)調階段，因此聯(lián)調內容就不在本文中做介紹。

    4 小結

    模擬機需要最大限度的反應電廠(chǎng)真實(shí)設計，模型開(kāi)發(fā)作為模擬機的核心直接決定了模擬機的仿真精度。工欲善其事，必先利其器，一套好的仿真建模工具是模型開(kāi)發(fā)工程師最好的助手，掌握好模擬機模型開(kāi)發(fā)的步驟與調試對于后續模擬機的運行維護以及技術(shù)改造升級具有重要意義。因為篇幅有限，所有不能將所有建模調試步驟一一詳述，僅希望以乏燃料池冷卻系統作為AP1000電站流體系統的例子讓更多的人了解模擬機模型的開(kāi)發(fā)過(guò)程。

    其他作者：陳森(1969-)，男，湖北人，高級工程師，哈爾濱工程大學(xué)工學(xué)學(xué)士、研究方向為全范圍模擬機開(kāi)發(fā)及應用維護。

    參考文獻：

    [1] GSE.JTOPMERET USER GUIDE，2005:30~270.

    [2] WEC.APP-SFS-M3-001_01,1版乏燃料池冷卻系統設計說(shuō)明書(shū)，2008.

    [3] WEC.APP-STS-J4-156，0版乏燃料池冷卻系統模型需求說(shuō)明書(shū)，2008.

    [4] WEC.APP-SFS-M6-001_03，3版乏燃料池冷卻系統P&ID圖，2008.

    [5] WEC. APP-SFS-M6-002_03，3版乏燃料池冷卻系統P&ID圖，2008.

    [6] WEC. APP-SFS-M6-003_03，3版乏燃料池冷卻系統P&ID圖，2008.

    [7] WEC.APP-SFS-M6K-001_01，1版乏燃料池冷卻系統過(guò)程流圖，2008.

    摘自《自動(dòng)化博覽》2010年第二期