朱坤帥 ,蔡燕霞 ,魯國(guó)瑞
(1.中國(guó)科學(xué)院國(guó)家空間科學(xué)中心北京100190;2.中國(guó)科學(xué)院大學(xué)計(jì)算機(jī)與控制學(xué)院,北京100049)
空間環(huán)境模式及其技術(shù)系統(tǒng)是開(kāi)展空間環(huán)境預(yù)報(bào)和服務(wù)的重要基礎(chǔ)條件[1]??臻g環(huán)境模式建立在大量空間探測(cè)數(shù)據(jù)的基礎(chǔ)上,是對(duì)獲取到的原始數(shù)據(jù)進(jìn)行科學(xué)建模、計(jì)算,獲得可供直接使用的產(chǎn)品數(shù)據(jù)??臻g環(huán)境領(lǐng)域研究的層層深入使空間環(huán)境模式的數(shù)量持續(xù)增長(zhǎng),致使模式管理難度不斷加大,以往的單一模式業(yè)務(wù)化已不能滿足現(xiàn)實(shí)需求。
本文從空間環(huán)境模式的現(xiàn)狀出發(fā),探討開(kāi)展空間環(huán)境多模式業(yè)務(wù)化的應(yīng)用需求,并借助于“無(wú)碼編程”技術(shù)與任務(wù)調(diào)度,提出構(gòu)建基于改進(jìn)Petri網(wǎng)的業(yè)務(wù)化空間環(huán)境模式集成平臺(tái)。
歐美發(fā)達(dá)國(guó)家在空間環(huán)境模式集成方面起步早,已經(jīng)研發(fā)出一系列的空間環(huán)境模式集成系統(tǒng),如SWMF(Space Weather Modeling Framework)、CCMC(Community Coordinated Modeling Center)、SPENVIC等[2]。在NASA(美國(guó)國(guó)家航空航天局)地球科學(xué)技術(shù)部計(jì)算技術(shù)項(xiàng)目的推動(dòng)下,美國(guó)密歇根大學(xué)的空間環(huán)境建模中心CSEM(Center for Space Environment Modeling)啟動(dòng)空間環(huán)境模式集成項(xiàng)目,建成了SWMF空間環(huán)境模式計(jì)算集成框架,實(shí)現(xiàn)從太陽(yáng)表面到地球上層大氣圈各區(qū)域空間物理過(guò)程的模式集成與分析應(yīng)用[3]。NASA與美國(guó)空軍、美國(guó)國(guó)家科學(xué)基金會(huì)(National Science Foundation,United States,NSF)、美國(guó)國(guó)家海洋與氣象管理局(National Oceanic and Atmospheric Administration,NOAA)等 8個(gè)機(jī)構(gòu)合作開(kāi)發(fā)的美國(guó)協(xié)同建模中心CCMC空間環(huán)境模式集成系統(tǒng)實(shí)現(xiàn)了模式與觀測(cè)數(shù)據(jù)的融合應(yīng)用,在空間環(huán)境預(yù)報(bào)領(lǐng)域得到了廣泛的使用[4]。比利時(shí)聯(lián)邦科學(xué)研究所在ESA(歐洲航天局)的資助下,成功開(kāi)發(fā)了空間環(huán)境信息系統(tǒng)(SPENVIS),實(shí)現(xiàn)了模式的耦合計(jì)算與結(jié)果可視化表達(dá),支持空間環(huán)境分析與事件預(yù)測(cè)[5]。
近年來(lái),隨著我國(guó)自主空間環(huán)境領(lǐng)域研究的不斷深入,國(guó)內(nèi)的研究機(jī)構(gòu)開(kāi)始注重空間環(huán)境模式研究成果的應(yīng)用與推廣,空間環(huán)境模式集成技術(shù)研究也隨之進(jìn)入起步發(fā)展階段。其中,借助于插件的動(dòng)態(tài)加載和注銷(xiāo)來(lái)實(shí)現(xiàn)模式集成方法[6],集成過(guò)程復(fù)雜,缺乏靈活性;使用模式封裝模塊來(lái)實(shí)現(xiàn)模式集成方法[7],集成涉及面廣,魯棒性差;設(shè)置模型服務(wù)組件,通過(guò)模型計(jì)算包和適配器模塊分別對(duì)數(shù)據(jù)及可執(zhí)行程序進(jìn)行模式集成方法[8],但未能實(shí)現(xiàn)業(yè)務(wù)化。
構(gòu)建空間環(huán)境模式集成平臺(tái)(Spatial Environment Mode Integration Platform,SEMIP),需要空間環(huán)境模式實(shí)現(xiàn)從研究級(jí)到業(yè)務(wù)級(jí)的轉(zhuǎn)化[9]。
空間環(huán)境模式是從事空間物理學(xué)基礎(chǔ)研究的科研人員研究分析空間數(shù)理規(guī)律,從中歸納總結(jié)出來(lái)的[10-11],具有專(zhuān)業(yè)性、針對(duì)性、獨(dú)立性和多樣性的特點(diǎn)?,F(xiàn)有的模式業(yè)務(wù)化以個(gè)體為單位,獨(dú)立集成;當(dāng)增加新的模式時(shí),需執(zhí)行大量重復(fù)工作,完成整個(gè)業(yè)務(wù)過(guò)程,大大增加了工作負(fù)擔(dān)。
“無(wú)碼集成”技術(shù)通過(guò)便捷的配置只需為新模式設(shè)置相關(guān)數(shù)據(jù)信息來(lái)實(shí)現(xiàn)新模式的業(yè)務(wù)化,精簡(jiǎn)了模式集成過(guò)程,實(shí)現(xiàn)從單模式到多模式的迅速擴(kuò)展[12]。
空間環(huán)境研究正趨向于向多數(shù)據(jù)、多參數(shù)綜合應(yīng)用領(lǐng)域發(fā)展,這意味著不同模式研究之間需要更加頻繁的數(shù)據(jù)交互共享。將原有的單一的模式集成置于統(tǒng)一的平臺(tái)下,將大大提高模式運(yùn)行效率,實(shí)現(xiàn)對(duì)多模式的有效管理。
從單一的空間環(huán)境模式運(yùn)行到多模式集成業(yè)務(wù)平臺(tái),模式運(yùn)行數(shù)目不斷增加,合理調(diào)度各模式,實(shí)現(xiàn)服務(wù)器最優(yōu)分配是實(shí)現(xiàn)平臺(tái)效率及有效性的關(guān)鍵,也是實(shí)現(xiàn)空間環(huán)境模式業(yè)務(wù)化的重要環(huán)節(jié)。
圖1展示了空間環(huán)境模式集成的物理總體過(guò)程。SEMIP通過(guò)從原始數(shù)據(jù)庫(kù)獲取原始數(shù)據(jù),經(jīng)過(guò)模式集成平臺(tái)“加工”處理后,把獲得的產(chǎn)品數(shù)據(jù)送入產(chǎn)品數(shù)據(jù)庫(kù)。該過(guò)程實(shí)現(xiàn)了業(yè)務(wù)化、自動(dòng)化的空間環(huán)境模式運(yùn)作。
為了更好地描述SEMIP平臺(tái)業(yè)務(wù)化運(yùn)行過(guò)程,文中引入圖形化建模工具——Petri網(wǎng)對(duì)SEMIP進(jìn)行抽象建模[13-14]。經(jīng)典Petri網(wǎng)(如圖2)由庫(kù)所(Place)、變遷(Transition)、令牌(Token)和有向弧(Connection)等元素組成,具有自然、直觀、簡(jiǎn)單易懂和易用的特點(diǎn)[15]。
圖2 經(jīng)典Petri網(wǎng)
在經(jīng)典Petri網(wǎng)的基礎(chǔ)上提出了一種基于改進(jìn)Petri網(wǎng)的模式運(yùn)算模型。通過(guò)給令牌“著色”,并引入模式狀態(tài)評(píng)級(jí)R、模式運(yùn)行周期C、各變遷的功能函數(shù)及狀態(tài)函數(shù)等來(lái)完成構(gòu)建基于改進(jìn)Petri網(wǎng)的模式運(yùn)算模型。
模式運(yùn)算模型包括模型的靜態(tài)定義和模式運(yùn)作過(guò)程兩個(gè)部分[16],主要要素關(guān)系如下:
1)靜態(tài)定義
定義1:令牌集M={Mi∈M,i=1,2,....,n},改進(jìn)的Petri網(wǎng)給令牌賦予屬性,且每一個(gè)令牌Mi={mj∈Mi,i,j=1,2,...n}由更小的令牌子屬性組成,隨著令牌在運(yùn)作過(guò)程中不斷被轉(zhuǎn)移,令牌屬性也隨之不斷變化,變化的令牌屬性決定著變遷是否進(jìn)行下去;在一次過(guò)程中,不存在完全相同的兩個(gè)令牌,如:?x≠y且x,y∈{1,2,....,n},Mx≠M(fèi)y。
定義2:庫(kù)所集P={Pi∈P,i=1,2,...n}是變遷前后令牌的庫(kù)所,若前后一致,則表示該變遷未變更庫(kù)所。
定義 3:活動(dòng)集T={Ti∈T,i=1,2,...n}涵蓋了整個(gè)模型運(yùn)作中所有變遷活動(dòng),且每一個(gè)變遷Ti={tj∈Ti,i,j=1,2,...n}由更小的子變遷組成。同一過(guò)程中的各活動(dòng)不存在并發(fā)行為,即?x≠y且x,y∈{1,2,....,n},Tx?Ty=?。
定義 4:模式狀態(tài)評(píng)級(jí)R={1,2,...,n}表示模式狀態(tài),狀態(tài)不同,模式運(yùn)行所在的庫(kù)所不同。例如,SEMIP平臺(tái)中,評(píng)級(jí)R為服務(wù)器當(dāng)前狀態(tài),根據(jù)服務(wù)器狀態(tài)的不同,新加入的模式任務(wù)將會(huì)分配到狀態(tài)最佳服務(wù)器中去,且行為不可逆。
定義5:模式運(yùn)行周期C=c,c值為固定值,模式確定后,則不能更改。
2)運(yùn)作過(guò)程
定義6.模型運(yùn)作過(guò)程W=(M,P,T,R,C)分3步完成,分別為準(zhǔn)備前期P-(如圖3所示)、準(zhǔn)備后期P+(如圖4所示)和周期運(yùn)行Wˉ(如圖5所示)。P-=(M,P,T,R)為模式準(zhǔn)備必須的前期條件,通過(guò)變遷T的作用,令牌M屬性隨著庫(kù)所P的變更發(fā)生變化,模式狀態(tài)R作為M屬性中的一部分也發(fā)生作用。同樣,P+=(M,P,T)則為模式準(zhǔn)備必備的后期條件。模式狀態(tài)R在前期影響庫(kù)所P后就失去作用。Wˉ =(P-,P+,C)實(shí)現(xiàn)模式的周期運(yùn)作過(guò)程。
圖3 P-活動(dòng)過(guò)程
至此,完成了模式運(yùn)算模型的闡述,接下來(lái),將該模型應(yīng)用于SEMIP平臺(tái)的設(shè)計(jì)與實(shí)現(xiàn)中。
圖4 P+活動(dòng)過(guò)程
圖5 Wˉ活動(dòng)過(guò)程
圖6展示了空間環(huán)境模式集成平臺(tái)的模型架構(gòu)設(shè)計(jì)。SEMIP平臺(tái)分為“無(wú)碼集成”和“任務(wù)調(diào)度”兩大主要模塊。
圖6 模型架構(gòu)設(shè)計(jì)
根據(jù)集成的時(shí)間節(jié)點(diǎn),無(wú)碼集成又可以分為模式運(yùn)行前的輸入準(zhǔn)備和模式運(yùn)行后的輸出準(zhǔn)備,具體實(shí)現(xiàn)過(guò)程如圖7所示。
圖7 無(wú)碼集成
輸入準(zhǔn)備時(shí),首先獲取所需的模式信息,配置XML數(shù)據(jù)模板及相關(guān)原始數(shù)據(jù)的元數(shù)據(jù)表,等待指令后能夠規(guī)范自動(dòng)生成所需要的輸入數(shù)據(jù)文件,完成輸入準(zhǔn)備工作;輸出準(zhǔn)備時(shí),配置XML數(shù)據(jù)模板及相關(guān)產(chǎn)品數(shù)據(jù)元數(shù)據(jù)表,等待獲取到的輸出數(shù)據(jù)文件,完成解析,最終實(shí)現(xiàn)將產(chǎn)品數(shù)據(jù)存放入指定的產(chǎn)品數(shù)據(jù)庫(kù)中[17]。
任務(wù)調(diào)度負(fù)責(zé)服務(wù)器的初始化分配及模式的周期運(yùn)行。當(dāng)新的模式集成時(shí),SEMIP平臺(tái)會(huì)根據(jù)實(shí)時(shí)的服務(wù)器集群的運(yùn)行狀態(tài),評(píng)級(jí)最優(yōu)的服務(wù)器分配給該模式使用;將配置好的新模式開(kāi)啟,SEMIP平臺(tái)會(huì)自動(dòng)根據(jù)模式需求的運(yùn)行周期,定期啟動(dòng)模式,實(shí)現(xiàn)持續(xù)7*24業(yè)務(wù)化運(yùn)行[18]。
模型中W=(M,P,T,R,C)表示模式運(yùn)算過(guò)程。以“F10.7未來(lái)27日預(yù)報(bào)模式”模型為例,進(jìn)行模型應(yīng)用研究[19]。模式運(yùn)算模型主要定義要素如表1、表2所示。
表1 令牌及其屬性含義
表2 任務(wù)及其具體描述
首先,完成平臺(tái)的輸入準(zhǔn)備P-初始化。根據(jù)模式的物理功能,設(shè)置模式編號(hào)、起算日期;根據(jù)當(dāng)前服務(wù)器集群評(píng)級(jí)R,為新模式分配目標(biāo)服務(wù)器。配置原始數(shù)據(jù)元數(shù)據(jù)及對(duì)應(yīng)XML輸入模板。
其次,完成平臺(tái)的輸入準(zhǔn)備P+初始化。配置產(chǎn)品數(shù)據(jù)元數(shù)據(jù)及對(duì)應(yīng)XML輸出模板。
接著,開(kāi)啟任務(wù)調(diào)度。根據(jù)設(shè)定的起算日期及對(duì)應(yīng)的模式編號(hào)、輸入模板、原始數(shù)據(jù)元數(shù)據(jù),生成輸入數(shù)據(jù)文件,啟動(dòng)新模式,獲取輸出數(shù)據(jù)文件,解析文件后,根據(jù)產(chǎn)品數(shù)據(jù)元數(shù)據(jù)及輸出模板,實(shí)現(xiàn)產(chǎn)品數(shù)據(jù)入庫(kù)工作,完成模式的一次運(yùn)作。
最后,任務(wù)調(diào)度進(jìn)入Wait過(guò)程,經(jīng)24小時(shí)(F10.7未來(lái)27日預(yù)報(bào)模式運(yùn)算周期為24 h)后,任務(wù)調(diào)度重新開(kāi)啟。
SEMIP平臺(tái)基于BS架構(gòu),采用Java、JSP語(yǔ)言開(kāi)發(fā),運(yùn)行在Windows系統(tǒng)平臺(tái)上;模式程序采用多種編程語(yǔ)言實(shí)現(xiàn),如 Fortran、C、Matlab、Java,運(yùn)行在Linux/Windows系統(tǒng)平臺(tái)上;原始數(shù)據(jù)庫(kù)/產(chǎn)品數(shù)據(jù)庫(kù)采用Oracle11g數(shù)據(jù)庫(kù)。系統(tǒng)框架如圖8所示。
圖8 空間環(huán)境模式集成平臺(tái)系統(tǒng)框架
文中以F10.7未來(lái)27日預(yù)報(bào)模式為例進(jìn)行實(shí)驗(yàn)驗(yàn)證。
F10.7未來(lái)27日預(yù)報(bào)模式采用時(shí)間序列模型中的自回歸方法開(kāi)展F10.7中期預(yù)報(bào)研究。該模式利用自回歸算法處理前期若干天的F10.7觀測(cè)日值,得到未來(lái)27天的F10.7日值。目前,F(xiàn)10.7觀測(cè)日值數(shù)據(jù)存儲(chǔ)在原始數(shù)據(jù)庫(kù)中,并每天持續(xù)更新。
1)分析物理模式。對(duì)F10.7未來(lái)27日預(yù)報(bào)模式進(jìn)行分析,獲得該模式的原始數(shù)據(jù)元數(shù)據(jù)、輸入文件、輸出文件、產(chǎn)品數(shù)據(jù)元數(shù)據(jù)、運(yùn)行周期,并匹配相對(duì)應(yīng)的輸入模板(如圖9所示)、輸出模板。
圖9 F10.7觀測(cè)日值輸入模板、輸入數(shù)據(jù)
2)模式集成。通過(guò)友好界面,實(shí)現(xiàn)“無(wú)碼集成”如圖10所示。通過(guò)配置F10.7未來(lái)27日預(yù)報(bào)模式的數(shù)據(jù)源、模板、程序等信息,添加該模式到平臺(tái)任務(wù)界面如圖11所示。
圖10 無(wú)碼集成配置
圖11 平臺(tái)管理
3)任務(wù)調(diào)度。開(kāi)啟任務(wù)調(diào)度。根據(jù)模式運(yùn)行周期,定時(shí)發(fā)送任務(wù)信息至相關(guān)服務(wù)器,完成模式運(yùn)行工作,并將獲取到的產(chǎn)品數(shù)據(jù)存入產(chǎn)品數(shù)據(jù)庫(kù)。
SEMIP平臺(tái)實(shí)現(xiàn)了便捷的無(wú)碼集成和7*24的業(yè)務(wù)化運(yùn)行,完成了實(shí)驗(yàn)設(shè)想。
“無(wú)碼集成”的模式集成方式體現(xiàn)了平臺(tái)的易配置性、可重用性和擴(kuò)展靈活性。平臺(tái)在添加新模式時(shí),只需獲取到該模式所需的數(shù)據(jù)來(lái)源(元數(shù)據(jù)),不需要編碼,僅僅通過(guò)友好的可視化界面就能實(shí)現(xiàn)新模式的集成工作,并且不會(huì)影響到其他現(xiàn)有模式。通過(guò)對(duì)前期空間環(huán)境模式的調(diào)研,模式所需的輸入模板、輸出模板大致相同,即使新模式下需要新的模板,也可以通過(guò)模板配置快速生成,大大增加了可重用性和擴(kuò)展靈活性。
任務(wù)調(diào)度的模式運(yùn)行管理方式體現(xiàn)了平臺(tái)的業(yè)務(wù)化。友好的可視化界面配置實(shí)現(xiàn)了平臺(tái)的“零基礎(chǔ)使用”,7*24的業(yè)務(wù)化運(yùn)行,實(shí)現(xiàn)了模式的自動(dòng)化機(jī)制;統(tǒng)一的平臺(tái)管理界面能夠更便捷的管理多模式,更高效地完成模式運(yùn)作。
本文提出的基于改進(jìn)Petri網(wǎng)的業(yè)務(wù)化空間環(huán)境模式集成平臺(tái)有效的解決了空間環(huán)境模式業(yè)務(wù)化程度低的問(wèn)題,實(shí)現(xiàn)了空間環(huán)境模式從研究級(jí)到業(yè)務(wù)級(jí)的轉(zhuǎn)化,從單模式到多模式的擴(kuò)展,從單調(diào)控到多調(diào)度的延伸,具備易配置性、可重用性和擴(kuò)展靈活性??臻g環(huán)境模式集成平臺(tái)的構(gòu)建,為航天領(lǐng)域從事空間環(huán)境的科研人員提供了應(yīng)用性強(qiáng),使用便捷的可用產(chǎn)品數(shù)據(jù),必將大大助力科研任務(wù)的持續(xù)開(kāi)展。
參考文獻(xiàn):
[1]劉四清,師立勤,鄒自明,等.空間天氣模式轉(zhuǎn)化與系統(tǒng)集成[C]//中國(guó)氣象學(xué)會(huì)2007年年會(huì)空間天氣科技成果向業(yè)務(wù)應(yīng)用的轉(zhuǎn)化分會(huì)場(chǎng)論文集,2007.
[2]鄒自明,佟繼周,熊森林,等.大數(shù)據(jù)時(shí)代空間科學(xué)領(lǐng)域的科研信息化實(shí)踐與成果[J].大數(shù)據(jù),2016(6):83-96.
[3]Gombosi T,Toth G,Sokolov I,et al.The Space Weather Modeling Framework(SWMF):Models and Validation[C]//COSPAR Scientific Assembly.2010.
[4]Chulaki A,Bakshi S S,Berrios D,et al.Communi?ty Coordinated Modeling Center(CCMC):Providing access to space weather models and research sup?port tools[C]//AGU Fall Meeting Abstracts,2011.
[5]Heynderickx D,Quaghebeur B,Speelman E,et al.ESA'sSpaceEnvironmentInformation System(SPENVIS)-A WWW interface to models of the space environment and its effects[C]//Aerospace Sciences Meeting and Exhibit,2013.
[6]黎芳芳.空間天氣預(yù)報(bào)模式集成可視化軟件研究與實(shí)現(xiàn)[D].北京:中國(guó)科學(xué)院研究生院,2011.
[7]鄭程.空間天氣模式集成框架設(shè)計(jì)與原型實(shí)現(xiàn)[D].北京:中國(guó)科學(xué)院研究生院,2012.
[8]鄒自明,高文健,張效信,等.空間天氣定量化預(yù)報(bào)系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)[J].科技導(dǎo)報(bào),2013,31(10):18-23.
[9]吳統(tǒng)文,宋連春,劉向文,等.國(guó)家氣候中心短期氣候預(yù)測(cè)模式系統(tǒng)業(yè)務(wù)化進(jìn)展[J].應(yīng)用氣象學(xué)報(bào),2013(5):533-543.
[10]鄒自明.太陽(yáng)風(fēng)擾動(dòng)的地磁響應(yīng)與空間環(huán)境應(yīng)用模式集成[D].北京:中國(guó)科學(xué)技術(shù)大學(xué),2014.
[11]王赤.空間天氣數(shù)值預(yù)報(bào)模式研究[C]//中國(guó)空間科學(xué)學(xué)會(huì)空間物理學(xué)專(zhuān)業(yè)委員會(huì)全國(guó)日地空間物理學(xué)研討會(huì),2013.
[12]何曉亮,李國(guó)敏.無(wú)碼技術(shù)成為企業(yè)信息化新趨勢(shì)[N].科技日?qǐng)?bào),2011-01-21(003).
[13]劉聰.基于Petri網(wǎng)的港口供應(yīng)鏈業(yè)務(wù)流程仿真研究[D].大連:大連海事大學(xué),2014.
[14]李帛.Petri網(wǎng)在工作流建模中的分析與應(yīng)用[J].信息與電腦:理論版,2012(12):178-179.
[15]袁崇義.Petri網(wǎng)原理[M].北京:電子工業(yè)出版社,1998.
[16]周朝揮,蔡燕霞,魯國(guó)瑞.信牌驅(qū)動(dòng)式Web數(shù)據(jù)采集模型的應(yīng)用[J].計(jì)算機(jī)應(yīng)用,2016(S1):252-256.
[17]王宇光.基于數(shù)據(jù)模板的空間科學(xué)數(shù)據(jù)加工模型的設(shè)計(jì)與實(shí)現(xiàn)[D].北京:中國(guó)地質(zhì)大學(xué)(北京),2011.
[18]夏家莉,陳輝,楊兵.一種動(dòng)態(tài)優(yōu)先級(jí)實(shí)時(shí)任務(wù)調(diào)度算法[J].計(jì)算機(jī)學(xué)報(bào),2012(12):2685-2695.
[19]汪宏波,熊建寧,趙長(zhǎng)印,等.太陽(yáng)輻射指數(shù)F10.7的中期預(yù)報(bào)方法[J].天文學(xué)報(bào),2014,55(4):302-312.