田恒宇 ， 簡 季* ， 李為樂

(成都理工大學 a.地學空間信息技術(shù)國土資源部重點實驗室；b. 地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護國家重點實驗室，成都 610059)

?

基于ArcEngine的地震災(zāi)害快速評價系統(tǒng)的設(shè)計與實現(xiàn)

田恒宇a， 簡 季a*， 李為樂b

(成都理工大學 a.地學空間信息技術(shù)國土資源部重點實驗室；b. 地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護國家重點實驗室，成都 610059)

地震災(zāi)害的信息化管理是實現(xiàn)地震災(zāi)害預(yù)測、預(yù)報和預(yù)防的基礎(chǔ)。為了滿足評價過程中對不同模型數(shù)據(jù)與因子數(shù)據(jù)的管理需求，結(jié)合ArcEngine二次開發(fā)技術(shù)，設(shè)計、實現(xiàn)了集災(zāi)害數(shù)據(jù)分析、災(zāi)害評價模型管理、模型應(yīng)用與評價結(jié)果可視化等功能于一體的地震災(zāi)害快速評價系統(tǒng)，并將其應(yīng)用于蘆山地震滑坡災(zāi)害危險性評價過程。測試結(jié)果表明，該系統(tǒng)可以快速穩(wěn)定完成不同模型下數(shù)據(jù)的分析與輸出，有效地提高災(zāi)害信息提取的工作效率，為抗災(zāi)減災(zāi)贏取時間。

地震災(zāi)害評價； 滑坡； 地理信息系統(tǒng)； ArcEngine； 二次開發(fā)

0 引言

地震導(dǎo)致的山體滑坡是我國地質(zhì)災(zāi)害中最為嚴重的一類[1]，國內(nèi)、外有不少針對其成災(zāi)機理及對預(yù)測的研究工作，其評價模型經(jīng)歷了定性模型、半定量模型和定量模型的發(fā)展，定量模型基于控制災(zāi)害要素與地質(zhì)災(zāi)害之間的數(shù)學表達，采用二元或多元回歸等方法獲得地質(zhì)災(zāi)害的評價結(jié)果[2]。法國專家提出的ZERMOS理論，該理論利用兩種控制因素建立滑坡危險性分區(qū)的數(shù)學模型，并以法國局部山區(qū)為例，進行滑坡危險性分區(qū)研究[3-4]；Finlay P.J.和Fell Robin[5-6]對香港滑坡進行了災(zāi)害風險研究，建立了預(yù)測滑坡災(zāi)害水平運動距離的多元回歸模型；John R Dymond等[7]用GIS研究了暴雨條件下滑坡災(zāi)害發(fā)生的概率，建立了評價模型；李家存等[8]建立的用于區(qū)域滑坡危險性評價的概率指數(shù)模型；李為樂等[9]建立了汶川地震次生地質(zhì)敏感性評價模型。

數(shù)據(jù)處理與表達是各種評價模型應(yīng)用的重要環(huán)節(jié)，因此，GIS和計算機技術(shù)的發(fā)展為地震地質(zhì)災(zāi)害提供了平臺，可以將GIS的空間分析與預(yù)測模型相結(jié)合，完成地質(zhì)災(zāi)害各影響因子的空間疊加，快速進行災(zāi)害的危險性預(yù)測。GIS技術(shù)作為工具，主要是利用其空間數(shù)據(jù)庫、空間分析以及可視化等功能，針對輸入的各評價因子的柵格圖層作某種函數(shù)疊加運算，從而很方便的得出地質(zhì)災(zāi)害易發(fā)性、危險性和損失等方面的評價及相應(yīng)的災(zāi)害區(qū)劃[10]。

在國內(nèi)已經(jīng)有不少城市或地區(qū)針對地質(zhì)災(zāi)害建立起了相應(yīng)的管理信息系統(tǒng)[11]，但系統(tǒng)的設(shè)計實現(xiàn)大多面向于某一特定區(qū)域或者某一個模型的應(yīng)用，缺少應(yīng)用區(qū)域的擴展性和災(zāi)害模型的選擇性。作者從地震地質(zhì)災(zāi)害滑坡數(shù)據(jù)與模型管理角度對系統(tǒng)進行設(shè)計，利用GIS技術(shù)的圖形圖像功能，實現(xiàn)滑坡預(yù)測結(jié)果的直觀可視化表達，用戶可以方便地對不同模型進行編輯保存，完成模型中各種因子數(shù)據(jù)的管理操作。

1 需求分析

地震災(zāi)害評價預(yù)測數(shù)據(jù)處理涉及到空間分析操作，根據(jù)用戶需求，運行速度快、數(shù)據(jù)處理效率高、使用簡潔是系統(tǒng)開發(fā)的總體要求。結(jié)合系統(tǒng)用例模型圖(圖1)，總結(jié)歸納出該系統(tǒng)需要實現(xiàn)的主要功能包括：

1)因子數(shù)據(jù)讀取。其中涉及多種格式單波段柵格數(shù)據(jù)，以及用戶對柵格值的區(qū)間劃分，新值的輸入等，數(shù)據(jù)的顯示以及操作以表格的方式呈現(xiàn)。

2)因子數(shù)據(jù)重分類。系統(tǒng)需要能夠依照用戶輸入的表格及輸出路徑等參數(shù)完成重分類工作。

3)不同影響因子的空間疊加。用戶可以自由選擇因子圖層，將其添加到評價模型式中，并根據(jù)模型公式完成致災(zāi)因子疊加運算得到評價結(jié)果。

4)災(zāi)害模型的選擇、編輯、存儲。用戶在進行區(qū)域地震災(zāi)害評價過程中可能會用到多種災(zāi)害模型，不同模型會產(chǎn)生不同的評價結(jié)果，系統(tǒng)應(yīng)該滿足用戶對于多種災(zāi)害模型的選擇要求，并能使其完成模型公式的自定義。

5)地圖顯示、瀏覽、渲染、輸出。系統(tǒng)應(yīng)提供地圖顯示以及放大、縮小、平移等瀏覽功能，還有根據(jù)柵格顏色顯示自定義等數(shù)據(jù)渲染及輸出功能。

圖1 地震災(zāi)害快速評價系統(tǒng)用例模型圖Fig.1 The use case diagram of rapid earthquake disaster evaluation system

2 系統(tǒng)設(shè)計

基于ArcEngine的地震災(zāi)害快速評價系統(tǒng)屬于空間型信息系統(tǒng)，其總體設(shè)計目標是以地理信息系統(tǒng)技術(shù)為核心，使用C#語言在ArcEngine及Microsoft Visual Studio 2010平臺上進行二次開發(fā)，緊密結(jié)合地震災(zāi)害評價模型，實現(xiàn)滑坡等地質(zhì)災(zāi)害的快速評價預(yù)測，并結(jié)合地圖可視化方案呈現(xiàn)地理信息評價結(jié)果，為決策人員提供決策支持。

2.1 系統(tǒng)結(jié)構(gòu)設(shè)計

地震災(zāi)害快速評價系統(tǒng)設(shè)計用于滑坡等災(zāi)害模型的管理及應(yīng)用，利用ArcEngine中封裝的各種功能組件作為軟件的開發(fā)工具，災(zāi)害評價模型可以被寫入系統(tǒng)代碼中實現(xiàn)計算分析，幫助地震災(zāi)害評估預(yù)測。系統(tǒng)整體框架(圖2)分為三層：①基于柵格、矢量數(shù)據(jù)及模型數(shù)據(jù)的數(shù)據(jù)庫層；②系統(tǒng)功能層；③用戶層。

系統(tǒng)不同部分的交互可以描述為：功能組件從空間數(shù)據(jù)庫中讀取模型文檔，用戶對讀取的數(shù)學模型進行編輯操作，根據(jù)模型因子再從數(shù)據(jù)庫中選擇因子數(shù)據(jù)，執(zhí)行后面的一系列基于模型的操作，包括因子分類賦值、評價模型建立。結(jié)果數(shù)據(jù)包括參數(shù)修改后的模型文檔、柵格數(shù)據(jù)、圖片，這些數(shù)據(jù)都保存在空間數(shù)據(jù)庫中。最終用戶通過交互界面實現(xiàn)數(shù)據(jù)庫中地圖數(shù)據(jù)的瀏覽。

圖2 地震災(zāi)害快速評價系統(tǒng)結(jié)構(gòu)圖Fig.2 The architechure of rapid earthquake disaster evaluation system

2.2 數(shù)據(jù)庫設(shè)計

滑坡等地震災(zāi)害的評價預(yù)測工作涉及多源、多時相數(shù)據(jù)，主要包括：

1)行政區(qū)劃、交通、地形地貌等基礎(chǔ)地理數(shù)據(jù)。

2)地層巖性、斷層、河流水系等地質(zhì)水文數(shù)據(jù)。

3)地震峰值加速度等災(zāi)害數(shù)據(jù)。

4)災(zāi)害評價模型數(shù)據(jù)。

為了便于系統(tǒng)對數(shù)據(jù)的安全、高效管理，系統(tǒng)數(shù)據(jù)采用ESRI的地理數(shù)據(jù)模型Geodatabase中的個人Geodatabase進行管理[12]。Geodatabase是一種采用標準關(guān)系數(shù)據(jù)庫技術(shù)來表現(xiàn)地理信息的數(shù)據(jù)模型，它可以將空間數(shù)據(jù)(矢量、柵格、影像、三維地形等)及其相關(guān)的屬性數(shù)據(jù)統(tǒng)一存放在工業(yè)標準的數(shù)據(jù)庫管理系統(tǒng)DBM中，對于小型的GIS項目和工作組來說，個人Geodatabase是非常理想的工具。

2.3 系統(tǒng)功能設(shè)計

根據(jù)系統(tǒng)結(jié)構(gòu)設(shè)計所述，本系統(tǒng)功能模塊及一系列工具組成，所有工具都被整合進一個圖形用戶界面(GUI)。在系統(tǒng)UML序列圖(圖3)描述了用戶端與系統(tǒng)各個單元之間的信息交互，功能模塊的信息輸入包括柵格數(shù)據(jù)、模型公式。不同單元之間有不同的信息反饋，包括表格、列表、文本、圖片。

圖3 系統(tǒng)UML序列圖Fig.3 The UML sequence diagram of the system

2.3.1 因子分類賦值模塊

因子分類賦值模塊主要功能是，參照已有的經(jīng)驗?zāi)Ｐ蛯υu價因子進行分類賦值。交互界面包括載入因子數(shù)據(jù)路徑、分類屬性表格及其功能按鈕、因子分類輸出路徑、因子分類賦值輸出功能按鈕，根據(jù)用戶需求，因子分類賦值分為兩種模式：①區(qū)間賦值；②唯一值賦值。

其組織結(jié)構(gòu)圖(圖4)表示了該模塊七個單元的相互關(guān)系，其中RasterWorkspaceFactory與RasterWorkspace是ArcEngine中用于柵格數(shù)據(jù)操作的類，RasterWorkspaceFactory實現(xiàn)了IrasterWorkspaceFactory接口，RasterWorkspace類實現(xiàn)了IRasterWorkspace接口，這兩個接口定義了用于打開和設(shè)置一個數(shù)據(jù)集對象，CreateRasterDataset方法用于產(chǎn)生一個新的基于文件的柵格數(shù)據(jù)集對象，這個數(shù)據(jù)格式可以為GRID、TIFF和ERDAS Image之一，讀取的數(shù)據(jù)由GetGeodataset單元轉(zhuǎn)換成ArcEngine中的geodataset類型。RasterStatistics作為ArcEngine DataSourceRaster類庫中的普通類，封裝了用于柵格像元值統(tǒng)計的方法，這些方法通過IrasterStatistics接口實現(xiàn)，其返回結(jié)果臨時保存到表文件中，由DataGridView單元進行讀取并顯示。該表文件與geodataset類型數(shù)據(jù)作為方法參數(shù)傳遞給RasterReclassOp單元。RasterRecalssOp作為ArcEngine GeoAnalyst類庫中的普通類，封裝了用于柵格數(shù)據(jù)分析的方法，包括Lookup、Reclass、ReclassByASCIIFile、ReclassByRemap、Slice，因子分類賦值模塊的分類賦值功能由Reclass方法實現(xiàn)[13]。ToolStrip是用于表格數(shù)據(jù)編輯的輔助工具集，包括插入行、刪除行、清空表格、增加區(qū)間數(shù)、減少區(qū)間數(shù)、唯一值賦值、數(shù)據(jù)顯示等工具。

圖4 因子分類賦值模塊組織結(jié)構(gòu)圖Fig.4 The organization chart of factor classification and assignment module

2.3.2 評價模型建立模塊

評價模型建立模塊包括：①圖層；②經(jīng)驗?zāi)Ｐ?；③模型編輯；④評價結(jié)果輸出與顯示四個組成部分。組織結(jié)構(gòu)圖(圖5)顯示了在程序執(zhí)行過程中各個單元的連接關(guān)系。LayerManager單元用于因子圖層的管理；Calculator包括數(shù)字鍵和基本運算符；ModelExpressionEditor是Visual Studio 2010中的Textbox控件，此單元可以讀取LayerManager中的圖層名稱和ModelManager中的模型公式，ToolStrip工具為用戶提供了添加、刪除、顯示因子圖層按鈕；RasterMapAlgebraOp單元是ArcEngine SpatialAnalyst類庫中用于柵格運算的組件類，該類封裝了ImapAlgebraOp接口，接口提供三個方法BindRaster、Execute、UnbindRaster，當用戶執(zhí)行計算輸出時，BindRaster方法將傳遞的Geodataset綁定到ModelExpression的因子名稱，然后Execute方法讀取ModelExpression開始運算，并將運算結(jié)果保存到指定路徑中。

圖5 評價模型建立模塊組織結(jié)構(gòu)圖Fig.5 The organization chart of evaluation model building module

用戶可以通過該模塊自由靈活的完成模型的新建、修改、刪除、保存，實現(xiàn)多模型管理的需求。如圖6所示，添加模型至經(jīng)驗?zāi)Ｐ椭?，可以在模型編輯中做任意調(diào)整修改，此處將該模型公式中的“因子1”刪除，并重新保存修改后的模型公式，命名為“自定義模型”。

圖6 評價模型管理系統(tǒng)界面圖Fig.6 The interface of evaluation model management

2.3.3 數(shù)據(jù)顯示與渲染模塊

數(shù)據(jù)顯示與渲染窗口為用戶提供了圖層管理、地圖瀏覽、坐標定位、距離測量、地圖渲染、圖片導(dǎo)出、地圖保存功能。為了讓用戶快速方便的對因子分類賦值、評價模型建立兩個模塊中的數(shù)據(jù)執(zhí)行顯示操作，共設(shè)置了五個功能按鈕關(guān)聯(lián)數(shù)據(jù)顯示與渲染窗口，分別用于當前載入數(shù)據(jù)顯示、分析結(jié)果數(shù)據(jù)顯示。在地圖渲染中，用戶可以對不同分類區(qū)間設(shè)置不同顏色以增強差異顯示。

3 系統(tǒng)實現(xiàn)與應(yīng)用

在本文所建立的評價系統(tǒng)中，圖層管理在圖形數(shù)據(jù)庫的管理中是一主要的特征，此種情況下，圖層為空間處理和分析的基本單元。在由ArcGIS創(chuàng)建的個人數(shù)據(jù)庫中，把空間數(shù)據(jù)按照評價因子圖層進行劃分管理，其中包括過程數(shù)據(jù)以及最終因子數(shù)據(jù)，數(shù)據(jù)類型可以是矢量，也可以是柵格。以汶川地震次生地質(zhì)敏感性評價模型為例，由于蘆山地震震后災(zāi)區(qū)影像資料匱乏，為了能大致掌握地震地質(zhì)災(zāi)害可能的分布范圍，利用該模型對“4·20”蘆山地震地質(zhì)災(zāi)害的空間分布進行快速評價，對應(yīng)急救災(zāi)有一定的指示意義。汶川地震地質(zhì)災(zāi)害敏感性多元回歸模型公式如下[9]：

(1)

式中：Ls為崩塌滑坡發(fā)生可能性，取值范圍為0～1；d為到發(fā)震斷層的距離；α為地形坡度；L為地層巖性；r為到水系距離；h為海拔高度；A為地震峰值加速度(PGA)。

在將模型評價因子載入系統(tǒng)前，數(shù)據(jù)的預(yù)處理及建立數(shù)據(jù)庫工作在ArcGIS中進行，包括坐標系定義、緩沖區(qū)分析、矢量柵格數(shù)據(jù)轉(zhuǎn)換等。

以地形坡度因子為例[9]，按照分類賦值標準將其值分為六個區(qū)間(表1)，執(zhí)行因子分類輸出并將分類結(jié)果添加到評價模型建立模塊；按照此步驟加載其余五個評價因子數(shù)據(jù)；選擇經(jīng)驗?zāi)Ｐ凸芾韱卧械你氪Ｐ?，模型編輯單元中出現(xiàn)此模型的表達式；在圖層單元中雙擊已添加的地形坡度數(shù)據(jù)，將其加載到模型表達式中；依次加載其余五個評價因子數(shù)據(jù)；選擇評價結(jié)果輸出路徑并計算輸出執(zhí)行程序；查看分析結(jié)果并在數(shù)據(jù)顯示界面中唯一值渲染(圖7)。

表1 汶川地震次生地質(zhì)敏感性評價模型中坡度分類賦值表[9]

圖7 蘆山地震次生地質(zhì)災(zāi)害預(yù)測評價流程圖Fig.7 The flow chart of evaluation process

圖8為系統(tǒng)操作界面，為了對比系統(tǒng)輸出結(jié)果，將汶川地震地質(zhì)災(zāi)害敏感性多元回歸模型中的地形坡度、海拔高度、到發(fā)震斷層距離因子做三次參數(shù)修改并保存至經(jīng)驗?zāi)Ｐ椭?，分別命名為“自定義模型1”、“自定義模型2”、“自定義模型3”，并計算出三次結(jié)果，生成結(jié)果渲染圖(圖9)，其差異反映出三種模型對同一地區(qū)不同評價結(jié)果。從模型修改、保存到因子分類與疊加再到數(shù)據(jù)顯示與渲染，整個過程系統(tǒng)運行流暢，數(shù)據(jù)分析準確。

圖8 地震災(zāi)害快速評價系統(tǒng)界面圖Fig.8 The interface of rapid earthquake disaster evaluation system

4 結(jié)束語

作者基于ArcEngine二次開發(fā)，在Visual Studio2010平臺上利用C#語言構(gòu)建了地震災(zāi)害快速評價系統(tǒng)，該系統(tǒng)從地理信息系統(tǒng)角度，將空間數(shù)據(jù)管理、空間分析模型融入到地震災(zāi)害評價中，實現(xiàn)了對不同致災(zāi)因子數(shù)據(jù)和評價模型的管理操作，極大地簡化了在評價過程中數(shù)據(jù)處理流程，為滑坡模型實驗研究提供有效依據(jù)，能夠滿足地震災(zāi)害評價的應(yīng)用需求。但其系統(tǒng)框架是初步的，在面向不同地區(qū)的災(zāi)害風險評估指標體系的建立、評價模型的表達等方面需要不斷的完善。

圖8 三種模型輸出結(jié)果對比圖Fig.8 Comparison of output results in three models

[1] 蘇歡, 易武, 孟召平,等. 基于MAPGIS的滑坡時間預(yù)測預(yù)報系統(tǒng)框架 設(shè)計 [J]. 災(zāi)害與防治工程, 2006(02):54-58.

SU H, YI W,MENG Z P,et al. Frame Design of Time Forecasting and Predicting System for Landslide Based on MAPGIS[J]. Disaster and Control Engineering, 2006(02):54-58.(In Chinese)

[2] 陳玉, 藺啟忠, 王欽軍. 區(qū)域地質(zhì)災(zāi)害評價模型綜述 [J]. 防災(zāi)科技學院學報, 2010(04): 42-45.

CHEN Y,L Q Z,WANG Q J. Review on evaluation model of the regional geohazards [J]. Journal of Institute of Disaster Prevention 2010(04): 42-45.(In Chinese)

[3] 徐為, 胡瑞林, 吳菲,等. 淺談我國的地質(zhì)災(zāi)害風險評估 [J]. 中國地質(zhì)災(zāi)害與防治學報, 2010(04): 126-129.

XU W,HU R L,WU F,et al. Discussion on some understanding about geo-hazards risk assessment in China[J]. The Chinese Journal of Geological Hazard and Control,2010(04): 126-129.(In Chinese)

[4] LISTO F D L R, CARVALHO VIEIRA B. Mapping of risk and susceptibility of shallow-landslide in the city of S?o Paulo, Brazil [J]. Geomorphology,2012,169-170：30-44.

[5] 楊義輝. 滑坡災(zāi)害風險評估方法及其應(yīng)用研究 [D].重慶：重慶交通大學,2014.

YANG Y H.Study on Landslide Hazard Risk Assessment Method and Its Application[D].Chongqing:Chongqing Jiaotong University,2014.(In Chinese).

[6] FELL R, COROMINAS J, BONNARD C, et al. Guidelines for landslide susceptibility, hazard and risk zoning for land use planning [J]. Engineering Geology,2008,102(3-4): 85-98.

[7] DYMOND J R, JESSEN M R, LOVELL L R. Computer simulation of shallow landsliding in New Zealand hill country [J]. International Journal of Applied Earth Observation and Geoinformation, 1999, 1(2): 122-31.

[8] 李家存, 李京, JASMI. 基于概率指數(shù)模型的區(qū)域滑坡危險性評價 [J]. 中國圖象圖形學報, 2007(08): 1471-1475.

LI J C,LI J,JASMI. Regional landslide susceptibility assessment based on probability lndex model [J]. Journal of Image and Graphics, 2007(08): 1471-1475.(In Chinese)

[9] 李為樂, 黃潤秋, 許強, 等. “4·20”蘆山地震次生地質(zhì)災(zāi)害預(yù)測評價 [J]. 成都理工大學學報：自然科學版, 2013(03):264-274.

LI W L,HUANG R Q,XU Q,et al. Rapid prediction of co-seismic landslides triggerd by Lushan earthquake, Sichuan, China [J]. Journal of Chengdu University of Technology:Sci&thnol,2013(03):264-274.(In Chinese)

[10]牛全福. 基于GIS的地質(zhì)災(zāi)害風險評估方法研究 [D].蘭州：蘭州大學, 2011.

NIU Q F.Study on the method of geological disaster risk assessment [D].Lanzhou: Lanzhou University,2011.(In Chinese)

[11]張宇明. 基于GIS技術(shù)的邊坡地質(zhì)災(zāi)害預(yù)測評價系統(tǒng)的研究及其應(yīng)用[D].大連：大連理工大學, 2005.

ZHANG Y M. application and study of prediction and evaluation system for slope geological hazard based on GIS technology [D].Dolian:Dalian University of Technology, 2005.(In Chinese)

[12]宋麗華, 沈明霞, 何瑞銀, 等. 基于ArcGIS的林業(yè)GIS空間數(shù)據(jù)庫建設(shè)的研究[J].計算機工程與設(shè)計,2008(19): 5117-5118+5122.

SONG L H,SHEN M X,HE R Y,et al. Study on building forestry GIS spatial database based on ArcGIS[J]. Computer Engineering and Design, 2008(19): 5117-5118+5122.(In Chinese)

[13]ESRI. ArcEngine developer's guide[Z]. 2004.

Design and realization of rapid earthquake hazard evaluation system based on ArcEngine

TIAN Heng-yua, JIAN Jia*, LI Wei-leb

(Chengdu University of Technology a.Key Laboratory of Geoscience Spatial Information Technology, Ministry of Land and Resources of the China, b.State Key Laboratory of Geohazard Prevention and Geoenvironment of Protection,Chengdu 610059,China)

The information management of earthquake hazards is the foundation of earthquake prediction, forecasting and prevention. In order to meet the administrative requirements of different models and factor data in the process of evaluation, we designed and implemented a rapid earthquake hazard evaluation system based on ArcEngine, which combined the disaster data analysis with the management and application of disaster models and the visualization of evaluation results. It was then used in the evaluation process of Lushan earthquake geological hazard. The test results prove that the system can complete the analysis and output task using different models fastly and steadily, which greatly improves the efficiency to extracting information from disasters areas.

geological hazard assessment; landslides; GIS; ArcEngine; secondary development

2015-07-13 改回日期：2016-03-21

國家自然科學青年基金資助項目(41202210)

田恒宇(1989-)，男，碩士，研究方向為地圖制圖與地理信息工程， E-mail：hengyu.hengyu@163.com。

*通信作者：簡季(1972-)男，教授，研究方向為高光譜遙感和地理信息工程，E-mail：jianji@21cn.com。

1001-1749(2016)05-0708-07

P 631.4

A

10.3969/j.issn.1001-1749.2016.05.22