對HIS中藥品采購模塊功能的改進及一鍵采購模式的擬建立

陳淑茜，何 訓(xùn)，呂世文（金華市中心醫(yī)院藥劑科，浙江金華 321000）

對HIS中藥品采購模塊功能的改進及一鍵采購模式的擬建立

陳淑茜＊，何 訓(xùn)，呂世文（金華市中心醫(yī)院藥劑科，浙江金華 321000）

目的：提高醫(yī)院藥房制訂藥品計劃采購量的合理性和科學(xué)性。方法：分析現(xiàn)有醫(yī)院信息系統(tǒng)（HIS）中采購藥品模塊功能的不足，采用增量系數(shù)法、ABC分析法、VEN分析法改進計劃庫存量的計算公式，再通過在系統(tǒng)中設(shè)置高儲值、低儲值以及各級包裝量來限制和調(diào)整計劃庫存量；在系統(tǒng)中擬建立一鍵采購功能，并以5種典型藥品進行模擬采購的操作實踐。結(jié)果：在原有計劃庫存量的計算公式上新增代表變化趨勢的參數(shù)、增量系數(shù)k1和計劃庫存系數(shù)k2；建立了只需輸入某類藥的計劃周期天數(shù)N1和購藥安全期N2，即可得出該類藥品計劃庫存量的一鍵采購模式。經(jīng)模擬操作實踐，與改進前比較，改進后計劃采購的藥品占用資金及空間減少，采購量可根據(jù)藥品包裝自動取整。結(jié)論：對HIS中藥品采購模塊功能的改進及一鍵采購模式的擬建立使制訂的藥品計劃采購量更加合理、科學(xué)、準確，既能保障藥品的及時供應(yīng)，又避免藥品積壓，并能有效降低人力、物力和資金成本。

醫(yī)院信息系統(tǒng)；藥品；計劃庫存量；增量系數(shù)法；ABC分析法；VEN分析法；一鍵采購

隨著醫(yī)療制度的不斷改革，醫(yī)院規(guī)模的不斷擴大，現(xiàn)今的大部分綜合性醫(yī)院使用的化學(xué)藥品種往往都在千種以上，如此多的藥品不可能每次都能一次性采購?fù)耆?。對藥庫管理人員及采購員來說，多長時間采購一次，每次采購多少個品種，每個品種采購多少數(shù)量，均是一項復(fù)雜且煩瑣的任務(wù)。若無完善規(guī)范的采購流程，一方面可能會造成某些不需要采購的藥品采購過多，不但降低庫存周轉(zhuǎn)率，而且致藥品積壓或失效；另一方面可能造成藥品采購不及時，影響臨床供應(yīng)，甚至直接影響到患者的治療效果[1]。筆者試圖通過改進醫(yī)院信息系統(tǒng)（HIS）原有藥品采購模塊中計劃采購量的計算公式，完善采購流程，提高采購計劃的工作效率和準確性。

1 現(xiàn)有HIS中藥品采購模塊設(shè)計方面存在的問題

1.1 現(xiàn)有計劃庫存量計算公式不合理

我院目前采用的計劃庫存量計算公式是V＝T×N/D（V：計劃庫存量，即計劃保有藥品的庫存數(shù)量；T：過去D天消耗的數(shù)量，均以最小單包裝計；N：計劃周期的天數(shù)；D：采樣周期的天數(shù)，現(xiàn)為14 d）。采用此公式計算時，筆者發(fā)現(xiàn)取樣周期短，不能有效體現(xiàn)近段時間藥品的使用情況，易過濾掉某些該領(lǐng)用的藥品，尤其是使用頻率低、波動大的藥品，從而影響臨床科室的治療用藥。而且對于各種藥品，由于設(shè)定其計劃周期天數(shù)是相同的，沒有根據(jù)藥品的價格高低、用量大小、重要程度來分類管理，容易造成庫存積壓、資金占用大甚至藥品過期的現(xiàn)象[2]。

1.2 HIS中未設(shè)置藥品庫存的高儲值和低儲值

此處的高儲值和低儲值是指某藥品的最高庫存數(shù)量和最低庫存數(shù)量。若在系統(tǒng)中未設(shè)置此值，會影響藥品的供應(yīng)和庫存計劃的判斷與制訂。比如對搶救藥品若沒有設(shè)置低儲值，只靠人工定期觀察補貨，容易出現(xiàn)遺漏；而對包裝體積大的藥品若沒有設(shè)置高儲值又容易導(dǎo)致占用太多空間。

1.3 HIS中未設(shè)置藥品的各級包裝量

現(xiàn)有程序不能自動計算出最接近的、方便購入和發(fā)放的整包裝量。比如亞甲藍注射液的小包裝是每盒5支，中包裝是每盒50支，若計劃采購數(shù)量是34支，目前在系統(tǒng)中自動計劃采購時的數(shù)量顯示34支，但一般為了便于整包裝發(fā)放，需人工在系統(tǒng)中修改數(shù)量為50支。

1.4 對于計劃采購量較小的藥品未設(shè)置篩除功能

使用現(xiàn)有程序時，每次計劃采購工作都需要通過人工篩除系統(tǒng)自動計劃采購但實際不需要計劃采購的藥品，增加了工作量和計劃時間。比如地特胰島素注射液現(xiàn)有庫存量為39支，而計劃庫存量設(shè)定為45支，系統(tǒng)會自動顯示計劃采購6支，但基于經(jīng)驗，現(xiàn)有庫存量已足夠，故暫時不需要采購，為了減少單次采購藥品種數(shù)，該藥的此次計劃采購單需由人工篩除。

2 對現(xiàn)有計劃庫存量計算公式的改進

2.1 通過增量系數(shù)法消除波動誤差

針對原有計算公式取樣周期短，又不能有效體現(xiàn)近期藥品消耗的變化趨勢這一缺點，現(xiàn)以取樣前2個采樣周期的消耗量T1、T2作為參考，將計劃庫存量計算公式改為V＝[（T1－T2）k1+T1]N/D（T1：前第1個采樣周期中的消耗量；T2：前第2個采樣周期中的消耗量；k1：可調(diào)整的增量系數(shù)，值越大越敏感，需通過多次試驗取得最合適的值）。采樣周期示意圖見圖1。

圖1 采樣周期示意圖Fig 1 Schematic diagram of sampling period

此算法引入了前2個采樣周期藥品消耗量的變化值（T1－T2），該值可以反映近期藥品消耗的變化趨勢：若T1－T2＞0，則表示消耗量在增加；若T1－T2=0，則表示消耗量沒有變化；若T1－T2＜0，則表示消耗量在減少。以（T1－T2）k1值來調(diào)整T1，取代原公式中單一的T值，以此制訂的V值會根據(jù)前段時間消耗量的變化趨勢而相應(yīng)調(diào)整，即對臨床需求量進行預(yù)估[3]。故與原有的計算公式相比，新公式體現(xiàn)了藥品需求的波動，也從一定程度上消除了波動誤差。

2.2 借助ABC分析法和VEN分析法確定各類藥品的計劃庫存系數(shù)

原有采購程序中，設(shè)定的各藥品的N值是相同的，未根據(jù)藥品的價格高低、用量大小、重要程度來分類管理。故針對此缺點引入ABC分析法和VEN分析法，以實現(xiàn)對不同藥品的差異性采購計劃制訂。

藥品ABC分類法的依據(jù)是庫存中各藥品每年消耗金額（某品種的年消耗量×單價）占所有藥品年消耗總金額的比值，即年消耗金額百分比。此法可用于分析1年或更短時間內(nèi)的藥品消耗量，幫助分析醫(yī)院資金是否得到了最有效的利用[4]。

ABC分析法無統(tǒng)一的標準，一般遵循以下規(guī)律：將品種百分數(shù)（某類藥品品種數(shù)占總的藥品品種數(shù)的百分數(shù)）為5%～15%、而年消耗金額百分數(shù)為60%～80%左右的前幾個藥品，確定為A類；將品種百分數(shù)為15%～25%、而年消耗金額百分數(shù)也為15%～25%的藥品，確定為B類；其余為C類，該類情況正與A類相反，其品種百分數(shù)為60%～80%，而年消耗金額百分數(shù)僅為5%～15%。根據(jù)ABC分析法的原則，少數(shù)的A類藥品占據(jù)了費用支出的絕大部分，對這類藥品應(yīng)加快資金周轉(zhuǎn)、減小每次采購量、縮短計劃周期。根據(jù)此原則，加強藥品庫存管理時，一般設(shè)定A類藥品的庫存周轉(zhuǎn)天數(shù)約為10 d，B類約為15 d，C類約為20 d。根據(jù)A、B類藥品所占年消耗藥品金額的百分比，理論上基本能滿足“85%以上藥品庫存周轉(zhuǎn)率少于10～15 d”的標準[5]。

采用ABC分析法管理藥品庫存時，由于只考慮藥品的出庫、消耗金額，存在一定的局限性，而在實際工作中還要考慮藥品的療效和用途。VEN分析法則是根據(jù)藥品對疾病治療的重要程度，將藥品分為V、E、N 3類。關(guān)鍵藥品（V類）：對于拯救生命或提供基本醫(yī)療保健至關(guān)重要；基本藥物（E類）：對不太嚴重但很重要的疾病有效；非基本藥物（N類）：用于輕微的或自限性疾病[6]。急救藥品、麻醉藥品、精神藥品以及一些必要的治療嚴重感染的抗菌藥物都屬于V類藥品。

將VEN分析法與ABC分析法結(jié)合后，對于A-N類藥品和C-V類藥品應(yīng)進行特別管控。A-N類藥品屬于治療輔助用藥，卻占用了大量資金，應(yīng)該降低其采購量，縮短計劃周期，故將A-N類藥品的庫存周轉(zhuǎn)天數(shù)設(shè)為5 d；C-V類藥品雖然價格低廉、周轉(zhuǎn)速度慢，但對于搶救生命至關(guān)重要，故將C-V類藥品的庫存周轉(zhuǎn)天數(shù)設(shè)為30 d。

制訂采購計劃時將N值與各類藥品的庫存周轉(zhuǎn)天數(shù)一一對應(yīng)，即可得到較合理的計劃庫存量。為了方便采購工作的實施，希望通過輸入1個N值，便可完成對各類不同藥品的各自的N值的確認。因此確定各類藥品的N值分別為A-N類5 d、A類10 d、B類15 d、C類20 d、C-V類30 d。現(xiàn)以A-N類藥品的N值為基準，轉(zhuǎn)換各類藥品相應(yīng)的計劃庫存系數(shù)k2，即確定A-N類藥品的k2值為1，則A類、B類、C類、C-V類藥品的k2值分別為2、3、4、6。因此，根據(jù)不同類藥品維持各自不同周轉(zhuǎn)天數(shù)的V的公式可改進為：V＝[（T1－T2）k1+T1]N1k2/D（N1：A-N類藥品的N值；k2：計劃庫存系數(shù)）。實際操作中，根據(jù)對AN類藥品N值的調(diào)整，可實現(xiàn)對各類藥品的N值的控制，見表1。

表1 各類藥品k2與N的關(guān)系Tab 1 Relationship between k2and N for each category

2.3 通過在系統(tǒng)中設(shè)置高儲值與低儲值來調(diào)整計劃采購量

對于難采購、波動大以及搶救藥品均在系統(tǒng)中設(shè)置低儲值L；對于包裝體積大、需要冷藏的藥品均在系統(tǒng)中設(shè)置高儲值H[7]；對于其他無特殊要求的藥品將低儲值默認設(shè)為0，高儲值默認設(shè)為∞。由于各類藥品的L和H值需根據(jù)人工經(jīng)驗判斷，那么初步計劃采購量（Qa）的取值就會出現(xiàn)以下幾種情況：（1）當V≤Va（現(xiàn)有庫存量）時，Qa＝0，即無須采購；（2）當V＞Va時，又可分為以下幾種情況：①若V≤L，則Qa＝L－Va；②若L＜V＜H，則Qa＝V－Va；③若V≥H，則Qa＝H－Va。

2.4 通過在系統(tǒng)中設(shè)置藥品的各級包裝量來調(diào)整計劃采購量

在HIS中，藥品的包裝規(guī)格信息是各不相同的。如果完全按照系統(tǒng)自動生成的采購數(shù)量來制訂計劃，會給上級部門以及供貨單位的發(fā)藥帶來不便，同時也不便于庫存部門收貨時的清點核實。如果HIS能自動識別與采購量相匹配的整包裝數(shù)量，并自動對計劃采購量進行適當調(diào)整，就能節(jié)約藥品采購與發(fā)放過程中更多的人力成本[8]。

故筆者在HIS中設(shè)置各藥品的最小包裝量為A1，其上一級包裝量為A2。為了自動篩除采購量極少的藥品，對計劃采購量（Qb）在0.1A1以下的藥品可不予采購。再根據(jù)四舍五入原則，最終Qb的取值可分為以下幾種情況：（1）若Qa≤0.1A1，則Qb＝0；（2）若0.1A1＜Qa＜0.5A1，則Qb＝A1；（3）若0.5A1≤Qa＜0.5A2，則對Qb＝（Qa/A1）的值小數(shù)點后四舍五入后取整，再乘以A1；（4）若Qa≥0.5A2，則對Qb＝（Qa/A2）的值小數(shù)點后四舍五入后取整，再乘以A2。

3 一鍵采購模式的擬建立

3.1 通過軟件自動分析計算藥品ABC分類

利用Excel軟件對醫(yī)院上一年度的藥品消耗進行統(tǒng)計分析，將藥品按全年消耗金額從高到低排序，計算每一品種藥品消耗金額占藥品總消耗金額的百分比、累計百分比、累計品種百分比，模擬結(jié)果見圖2。

根據(jù)ABC分析法的原則，結(jié)合以上分析得到的品種百分比和消耗金額百分比，由計算機自動確定藥品的ABC分類結(jié)果。

3.2 通過人工設(shè)置藥品VEN分類類別、各級包裝量和高、低儲值

圖2 藥品ABC分類模擬結(jié)果Fig 2 The simulation results ofABC analysis

藥品的ABC分類、藥名、規(guī)格、包裝單位可由ABC分析表直接導(dǎo)入到系統(tǒng)中；人工設(shè)置VEN分類結(jié)果，主要標注出A類中的N類藥品和C類中的V類藥品即可，其余可在系統(tǒng)顯示缺??；k2值可通過系統(tǒng)篩選A、B、C類和V、E、N類結(jié)果后進行批量設(shè)置，比如屬于C-V類的所有藥品，其k2值就可以批量設(shè)置為6；除特別的藥品由人工設(shè)置其高、低儲值外，其余各藥品的高、低儲值也可在系統(tǒng)中顯示缺省，系統(tǒng)會將這些藥品的高儲值默認為∞，低儲值默認為0；藥品各級包裝量信息可由系統(tǒng)直接導(dǎo)入，對只有1種包裝量（即平時發(fā)貨的常用包裝量）的藥品，人工將其設(shè)置為最小包裝量A1，對于其余包裝量信息缺省的藥品的A2系統(tǒng)默認為∞；k1和D這2個值需要正式啟用程序后逐步調(diào)整，直到得出最合適的值，因為所有藥品的k1、D值都取相同值，所以不需要逐個輸入，并且該項目不能缺省。程序設(shè)計時需要的數(shù)據(jù)和需實現(xiàn)的各項功能見圖3。

圖3 藥品設(shè)置數(shù)據(jù)表及各功能要求Fig 3 Drug setting data table and functional requirements

3.3 N值和購藥安全期的確定

正式采購時輸入A-N類藥品的N1，將N1代入公式V＝[（T1－T2）k1+T1]N1k2/D，可計算得出各藥品的V。但由于購藥品種數(shù)較多，為了盡量減少購藥品種數(shù)，又不至于造成缺貨的風(fēng)險，需要憑借人工經(jīng)驗確定購藥安全期（N2），以確保到下次購藥前的藥品用量足夠。此N2值可以取此購藥時間點到下次購藥時間點的間隔天數(shù)，也可取高于此的天數(shù)，使安全系數(shù)更高，但后種情況下在降低缺貨率同時可能相應(yīng)會增加庫存量。當Va≥[（T1－T2）k1+ T1]N2k2/D，即Va≥VN2/N1，也可理解成當V≤VaN1/N2時則無需采購；當V＞VaN1/N2時，則進入下一步程序。

將經(jīng)過購藥安全期篩選后的藥品V值再結(jié)合高、低儲值，可得出Qa；最后結(jié)合Qa與各級包裝量，得出最適合的整包裝倍數(shù)的Qb。

3.4 低儲警報的設(shè)置

結(jié)合系統(tǒng)軟件分析的數(shù)據(jù)和人工設(shè)置的數(shù)據(jù)，并定期對其進行維護，在日常購藥時，只需要輸入本次購藥的A-N類藥品的N1和N2，就能一鍵得出所有藥品的Qb。但應(yīng)注意在得到N2過程中很可能會篩掉小部分平時用量極少的搶救藥品，所以應(yīng)對這類藥品的低儲值設(shè)置安全警報，當現(xiàn)有庫存量小于低儲值時就發(fā)出安全警報；若此藥品未出現(xiàn)在系統(tǒng)自動一鍵采購的名單中，則需單獨采購。

4 新建采購程序的模擬操作實踐

4.1 5種藥品采購相關(guān)指標數(shù)據(jù)的計算

采用Excel軟件對全院上一季度的藥品消耗進行統(tǒng)計分析，自動確定藥品的ABC分類，取其中5種典型藥品進行模擬采購。暫定k1為1，D為14，若計劃下次購藥是在3 d后，現(xiàn)取A-N類藥品的N1為5，N2為4，其他相關(guān)指標數(shù)據(jù)分析計算結(jié)果見表2。

表2 5種藥品采購相關(guān)指標數(shù)據(jù)Tab 2 Relevant data of drug purchase index for 5 categories

在表2中同時列入采用原有方法的計劃采購量（Qp）以供比較，Qp＝T×N/D－Va，其中D為14 d，T為過去14 d內(nèi)藥品的消耗量。為盡量避免各類藥品不出現(xiàn)缺貨，所有藥品的N均取值為10 d。

4.2 對5種藥品采購相關(guān)指標數(shù)據(jù)的分析

4.2.1 消疲靈顆粒 因V≤VaN1/N2，無需制訂采購計劃。

4.2.2 0.9%氯化鈉注射液 因V＞VaN1/N2，需制訂采購計劃；V≥H，則Qa＝H－Va；0.5A1≤Qa＜0.5A2（A2為∞），則Qb＝（Qa/A1），取整后小數(shù)點后四舍五入后再乘以A1。

4.2.3 地特胰島素注射液 因V≤VaN1/N2，無需制訂采購計劃；但因此藥V＞Va，在經(jīng)N2篩選后，從采購計劃中篩除此藥。

4.2.4 西吡氯銨含漱液 因V＞VaN1/N2，需制訂采購計劃；L＜V＜H，則Qa＝V－Va；0.5A1≤Qa＜0.5A2（A2為∞），則Qb＝（Qa/A1），取整后小數(shù)點四舍五入后再乘以A1。

4.2.5 亞甲藍注射液 因V＞VaN1/N2，需制訂采購計劃；L＜V＜H，則Qa＝V－Va；Qa≥0.5A2，則Qb＝（Qa/A2），取整后小數(shù)點后四舍五入后再乘以A2。

4.3 數(shù)據(jù)分析結(jié)果

由表3可看出，新程序通過k2對各藥品設(shè)置了不同權(quán)重后，降低了占用大量資金卻只是輔助藥的A-N類藥品（如消疲靈顆粒）的采購量，同時增加了周轉(zhuǎn)速度慢但至關(guān)重要的C-V類藥品（如亞甲藍注射液）的采購量。就以上5種藥品而言，采用舊公式制訂采購計劃的占用資金為14 209元，采用新公式的占用資金為6 216元，資金占用量明顯下降；H的設(shè)置避免了包裝體積大的藥品（如0.9%氯化鈉注射液）占用過多貯存空間；N2的設(shè)置對不急于采購的藥品進行了篩選，如地特胰島素注射液，在實際采購工作中可能會有大批此類藥品被篩除，從而減少了單次的購藥品種數(shù)；而藥品包裝量的設(shè)置使最終計劃采購量自動取整，便于運輸與發(fā)放。

5 結(jié)語

針對舊程序存在的問題，筆者對現(xiàn)有計劃庫存量的計算方法進行了改進，通過增量系數(shù)法、ABC分析法、VEN分析法，從消耗金額及治療價值等方面較為系統(tǒng)地在HIS中預(yù)設(shè)了各類藥品的計劃庫存量；通過設(shè)定某些特殊藥品的高、低儲值避免了人工制訂采購計劃時易出現(xiàn)的斷貨或積壓問題；通過對藥品各級包裝量的設(shè)置節(jié)約了藥品采購與發(fā)放過程中的人力成本。

本文建立的一鍵采購模式的擬建立，利用計算機網(wǎng)絡(luò)的先進性，可快速準確地獲得藥品信息；通過維護少量參數(shù)，在無需人工干預(yù)的情況下，能自動生成較為合理的計劃庫存量，使其更科學(xué)、更系統(tǒng)，不但減少了人為誤差，還可降低人力與時間成本[9]。

總之，藥品的采購是醫(yī)院藥學(xué)部門工作中非常重要的環(huán)節(jié)，承擔著為臨床和患者提供安全藥品的使命。雖然影響藥品消耗的因素多且一般無規(guī)律可循，不可能制訂出完美的計劃采購量，但是可以通過不斷地優(yōu)化算法使其更為科學(xué)合理，有效提高藥學(xué)部門管理水平，為患者提供更好的藥學(xué)服務(wù)，提高醫(yī)院的經(jīng)濟效益與社會效益。

[1] 宋冰.關(guān)于醫(yī)院藥品庫存控制的探討[J].醫(yī)學(xué)信息，2010，23（12）：4728.

[2] 陳小彪.我院合理控制藥品庫存量的持續(xù)改進[J].臨床合理用藥雜志，2014，7（2）：180-181.

[3] 楊紅斌，馮丹.藥品申領(lǐng)計劃量算法的優(yōu)化[J].醫(yī)藥導(dǎo)報，2014，33（12）：1676-1678.

[4] 聶玲，朱月皓.帕雷托分類法在我院藥品庫存管理中的運用[J].淮海醫(yī)藥，2014，32（5）：489-490.

[5] 中國醫(yī)院協(xié)會.三級綜合醫(yī)院評審標準條款評價要素與方法說明[M].北京：人民衛(wèi)生出版社，2011：184.

[6]安靜，董占軍，邱志宏.基于ABC-VEN分析法的醫(yī)院藥品管理[J].中國藥房，2015，26（4）：503-505.

[7] 劉鵬.如何提高醫(yī)院藥房藥品周轉(zhuǎn)率[J].首都醫(yī)藥，2014（20）：8-9.

[8] 吳妙蓮，沈小慶，王玨.量化決策分析法在醫(yī)院藥房自動備藥系統(tǒng)中的應(yīng)用[J].中國藥學(xué)雜志，2005，40（23）：1836-1837.

[9] 黃旭.計算機技術(shù)在藥庫規(guī)范化管理中的實踐與探討[J].中國醫(yī)藥科學(xué)，2014，4（7）：159-161.

（編輯：劉 萍）

Function Improvement of Drug Purchase Module and Quasi Establishment of One-key-getting Purchase Mode Based on HIS

CHEN Shuqian，HE Xun，LYU Shiwen（Dept.of Pharmacy，Jinhua Central Hospital，Zhejiang Jinhua 321000，China）

OBJECTIVE：To improve the rational and scientific drug purchase quantity planned by hospital pharmacy.METHODS：The deficiency in present drug purchase module of hospital information system（HIS）was analyzed；the formula of planned purchase quantity was improved by incremental coefficient method，ABC analysis method and VEN analysis method.The high and low stored values and packaging volume were set to adjust planned drug purchase quantity.Based on the computer program，a onekey-getting procedure was quasi established，and five typical drugs were selected for simulation.RESULTS：Based on previous formula，the parameters representing change trend，incremental coefficient k1and inventory coefficient k2were added in new formula. The one-key-getting procedure could calculate the quantity of all drugs by simply entering the planning cycle N1and the safe period N2of the some category drugs into the system.After simulated operation practice，compared with pre-improvement，the occupation of funds was significantly decreased and the quantity of drug was rounded off according to package in post-improvement.CONCLUSIONS：Through function improvement of drug purchase module in HIS and the quasi establishment of one-key-getting purchase mode，planned drug purchase quantity is more reasonable，scientific and accurate.It ensures the clinical drug supply in time without causing a backlog，and effectively reduces the cost of manpower，material and capital.

Hospital information system；Drug；Planned purchase quantity；Incremental coefficient method；ABC analysis method；VEN analysis method；One-key-getting purchase

R95

A

1001-0408（2017）01-0087-05

2016-04-18

2016-07-14）

＊藥師。研究方向：藥庫管理。E-mail：chen_shu_qian@163.com

DOI10.6039/j.issn.1001-0408.2017.01.23