，，

(1.陜西國(guó)防工業(yè)職業(yè)技術(shù)學(xué)院，陜西 戶縣 710300；2.西安石油大學(xué)，西安 710065；3.西安摩爾石油工程實(shí)驗(yàn)室，西安 710065)

·試驗(yàn)研究·

基于ISO13679熱循環(huán)試驗(yàn)要求的套管特殊螺紋強(qiáng)度分析

白松1，張建兵2，杜志杰3

(1.陜西國(guó)防工業(yè)職業(yè)技術(shù)學(xué)院，陜西 戶縣 710300；2.西安石油大學(xué)，西安 710065；3.西安摩爾石油工程實(shí)驗(yàn)室，西安 710065)

依照ISO 13679《套管和油管接頭試驗(yàn)程序推薦作法》CAL Ⅲ C系熱循環(huán)試驗(yàn)程序，采用Abaqus有限元分析軟件，對(duì)JIP ?177.8 mm×10.36 mm規(guī)格P110鋼級(jí)新型套管特殊螺紋的強(qiáng)度進(jìn)行分析評(píng)價(jià)?？紤]接觸問題計(jì)算的非線性，建立了套管特殊螺紋熱力耦合有限元模型，對(duì)接頭的C系加載路徑進(jìn)行了簡(jiǎn)化，計(jì)算了接頭在10個(gè)階段的應(yīng)力場(chǎng)與溫度場(chǎng)耦合作用下的熱應(yīng)力。結(jié)果顯示：隨著溫度由低溫到高溫的變化，特殊螺紋部分螺牙、接箍中面等塑性變形加大；而局部由于高、低溫反復(fù)作用可能會(huì)產(chǎn)生組織的變化，使先前本應(yīng)屈服的區(qū)域強(qiáng)度反而提高，但特殊螺紋整體上在C系加載完成后仍能保持結(jié)構(gòu)完整；高溫狀態(tài)，材料的屈服強(qiáng)度降低，其變形更多的處于彈性變形階段，室溫狀態(tài)時(shí)特殊螺紋的關(guān)鍵部位發(fā)生屈服的可能性更大。通過計(jì)算機(jī)仿真評(píng)價(jià)方法，按照ISO 13679的試驗(yàn)內(nèi)容和程序?qū)μ厥饴菁y進(jìn)行評(píng)價(jià)，可以大幅節(jié)約時(shí)間與經(jīng)濟(jì)成本，可在一定程度上代替接頭的實(shí)物評(píng)價(jià)。

套管；特殊螺紋；強(qiáng)度；熱-力耦合；有限元

隨著油氣開采不斷向深井、超深井、大位移井、叢式井等新型復(fù)雜井發(fā)展，油氣田對(duì)油套管的性能要求也越來越高，而油套管特殊螺紋作為連接管柱的重要環(huán)節(jié)，其性能的可靠性至關(guān)重要。

一般分析特殊螺紋性能的方法有解析法計(jì)算、實(shí)物試驗(yàn)驗(yàn)證、有限元分析等。由于特殊螺紋幾何外形的復(fù)雜性，解析法在多種假設(shè)的條件下仍很難建立精確的計(jì)算模型。試驗(yàn)法主要以ISO 13679 《套管和油管接頭試驗(yàn)程序推薦作法》[1]，采用全尺寸試驗(yàn)，對(duì)特殊螺紋進(jìn)行性能進(jìn)行評(píng)價(jià)[2]。

ISO 13679對(duì)試驗(yàn)設(shè)備要求極高，對(duì)試驗(yàn)所用試樣、試驗(yàn)條件的要求苛刻。同時(shí)，國(guó)內(nèi)僅有國(guó)家石油管材質(zhì)量監(jiān)督檢驗(yàn)中心和國(guó)家石油天然氣管材工程技術(shù)研究中心等極少數(shù)單位具備開展該評(píng)價(jià)試驗(yàn)的硬件條件和檢驗(yàn)資質(zhì)。

近年來，國(guó)外的住友金屬、泰納瑞斯、GE石油天然氣集團(tuán)等[3-5]機(jī)構(gòu)的相關(guān)學(xué)者已將有限元分析技術(shù)應(yīng)用于特殊螺紋的氣密封性評(píng)價(jià)、疲勞設(shè)計(jì)、強(qiáng)度分析等方面。國(guó)內(nèi)高校學(xué)者如申昭熙、許志倩、祝效華、竇益華、高連新等[6-13]也已將計(jì)算機(jī)有限元分析用于特殊螺紋的研究工作之中。他們大多著眼于特殊螺紋的結(jié)構(gòu)對(duì)其某一方面的技術(shù)參數(shù)的影響規(guī)律的分析，這些研究對(duì)新型螺紋油套管的開發(fā)提供了很好的指導(dǎo)價(jià)值，但按照ISO 13679的試驗(yàn)內(nèi)容和程序?qū)τ吞坠芴厥饴菁y性能所進(jìn)行的系統(tǒng)性的仿真評(píng)價(jià)還基本沒有涉及。

本文以JIP ?177.8 mm×10.36 mm P110鋼級(jí)的套管特殊螺紋為分析對(duì)象，按照ISO 13679的C系熱循環(huán)試驗(yàn)程序和方法，對(duì)該套管螺紋的強(qiáng)度進(jìn)行分析。

1 C系熱循環(huán)試驗(yàn)簡(jiǎn)介

熱循環(huán)試驗(yàn)的目的是評(píng)估在服役條件下，以及拉伸和內(nèi)壓加速作用下的特殊螺紋在熱循環(huán)時(shí)泄漏的可能性，而在實(shí)際試驗(yàn)過程中特殊螺紋在熱力耦合作用下其強(qiáng)度早已發(fā)生了變化。

熱循環(huán)加載程序如圖1所示。對(duì)于試驗(yàn)高溫段，在要求溫度或更高溫度下至少保持5 min；對(duì)于試驗(yàn)低溫段，在要求溫度或更低溫度下至少應(yīng)保持5 min。對(duì)于CAL Ⅱ、Ⅲ級(jí)試驗(yàn)，試驗(yàn)高溫應(yīng)不小于135 ℃(275 ℉)；對(duì)于CAL Ⅳ級(jí)試驗(yàn)，試驗(yàn)高溫應(yīng)大于或等于180 ℃(356 ℉)。對(duì)于所有適用等級(jí)的試驗(yàn)，試驗(yàn)低溫應(yīng)小于或等于52 ℃(125 ℉)。每個(gè)循環(huán)最短時(shí)間為30 min，循環(huán)可連續(xù)進(jìn)行，也可間斷進(jìn)行。

室溫下最大機(jī)械載荷規(guī)定：①拉伸載荷取80%管體(或接箍)屈服載荷和80%特殊螺紋試驗(yàn)包絡(luò)線載荷之中的較小者，兩者都是基于材料室溫屈服強(qiáng)度得到的；②內(nèi)壓載荷取95% VME管體(或接箍)屈服載荷和95%特殊螺紋試驗(yàn)包絡(luò)線載荷之中的較小者，兩者都是基于80%拉伸載荷和室溫下材料屈服強(qiáng)度得到的。

試驗(yàn)高溫最大機(jī)械載荷規(guī)定：①內(nèi)壓與室溫下的內(nèi)壓相同；②拉伸載荷取90% VME管體(或接箍)屈服和90%特殊螺紋應(yīng)力試驗(yàn)包絡(luò)線載荷之中的較小者，兩者都是基于材料室溫屈服強(qiáng)度得到的。

2 C系載荷設(shè)計(jì)

考慮到實(shí)物試驗(yàn)操作的復(fù)雜性，試驗(yàn)周期長(zhǎng)，所以在利用有限元方法模擬C系試驗(yàn)時(shí)，簡(jiǎn)化試驗(yàn)步驟，縮短試驗(yàn)時(shí)間。將C系試驗(yàn)程序設(shè)計(jì)為室溫下的內(nèi)壓和拉伸交變作用下5次循環(huán)、在內(nèi)壓和拉伸作用一定時(shí)施加循環(huán)熱載荷50次、高溫下的內(nèi)壓和拉伸交變作用下5次循環(huán)，然后再重復(fù)熱循環(huán)50次和室溫下載荷循環(huán)5次，為方便載荷過渡，模仿實(shí)物試驗(yàn)情況設(shè)計(jì)了升溫、保載、卸載、冷卻等過程。針對(duì)JIP-7特殊螺紋的試驗(yàn)載荷如表1。

表1 套管特殊特殊螺紋C系熱循環(huán)試驗(yàn)載荷

表1(續(xù))

3 套管特殊螺紋熱力耦合模型

采用有限元方法模擬計(jì)算熱力耦合問題有兩類方法：①間接耦合法。是指按照順序進(jìn)行多次不同場(chǎng)分析，熱力耦合中先計(jì)算在溫度場(chǎng)中隔單元的節(jié)點(diǎn)溫度，然后將節(jié)點(diǎn)溫度當(dāng)作第2次計(jì)算的初始條件，在應(yīng)力場(chǎng)分析中實(shí)現(xiàn)耦合；②直接耦合法。指利用同一計(jì)算環(huán)境中包含多個(gè)場(chǎng)變量，利用能夠?qū)崿F(xiàn)不同場(chǎng)自由度的單元實(shí)現(xiàn)耦合，通過一次計(jì)算得到分析的結(jié)果。本文計(jì)算過程中采用準(zhǔn)靜態(tài)計(jì)算，分析步驟是在恒定溫度下設(shè)定的，而不是溫度場(chǎng)和應(yīng)力場(chǎng)的解分離，溫度場(chǎng)和位移場(chǎng)之間的兩個(gè)變量相互影響。

根據(jù)以前的試驗(yàn)規(guī)律，隨溫度升高材料的屈服強(qiáng)度降低，塑性變形加強(qiáng)，故假定在高溫段管體和接箍材料的屈服強(qiáng)度下降10%，塑性擴(kuò)大2倍。材料力學(xué)性能如表2。套管特殊螺紋管體和接箍材料的密度ρ=7.85×10-6kg/mm3，引入線膨脹系數(shù)、比熱容系數(shù)、泊松比、彈性模量等參數(shù)，如圖2。

表2 套管特殊特殊螺紋與溫度相關(guān)的材料力學(xué)性能

a 泊松比

b 彈性模量

c 線膨脹系數(shù)

d 比熱系數(shù)

為提高計(jì)算效率和精度，采用四節(jié)點(diǎn)雙線性軸對(duì)稱四邊形縮減積分單元CAX4R，而CAX4R單元是能夠同時(shí)承擔(dān)位移場(chǎng)與溫度場(chǎng)自由度的單元[14]。由于接觸問題本身邊界條件復(fù)雜，具有不確定性，因此采用縮減積分單元可能會(huì)出現(xiàn)的單積分點(diǎn)線性單元受力變形沒有產(chǎn)生應(yīng)變能的情況(零能量模式)，此時(shí)單元沒有剛度，所以不能抵抗變形，故采用沙漏控制hourglass control，施加虛擬的剛度以限制沙漏模式的擴(kuò)展[15]。在網(wǎng)格劃分時(shí)，考慮到隱式求解的計(jì)算成本與自由度數(shù)的平方成正比，故只在套管特殊螺紋的螺紋牙處、密封面、臺(tái)肩處局部細(xì)化。最終得到管體部分單元數(shù)5 043個(gè)，節(jié)點(diǎn)數(shù)5 366個(gè)，接箍部分單元數(shù)5 220個(gè)，節(jié)點(diǎn)數(shù)5 569個(gè)，有限元網(wǎng)格模型如圖3所示。

圖3 套管特殊螺紋有限元網(wǎng)格模型

4 結(jié)果及分析

套管特殊螺紋結(jié)構(gòu)上發(fā)生損傷的部位發(fā)生在不完整螺紋段、特殊螺紋端面頂點(diǎn)附近、接箍中面等，故在對(duì)特殊螺紋C系熱循環(huán)分析時(shí)，重點(diǎn)討論第1、17牙螺紋、密封面和臺(tái)肩面。

套管特殊螺紋在C系熱循環(huán)下關(guān)鍵位置等效應(yīng)力值如表3，等效應(yīng)力云圖如圖4。

第1次室溫下載荷循環(huán)5次后等效應(yīng)力分布如圖4a，可見特殊螺紋外螺紋最后幾牙螺紋局部發(fā)生了屈服，但整體應(yīng)力分布比較均勻、合理；加熱升溫后等效應(yīng)力分布如圖4b，明顯發(fā)現(xiàn)接箍中部、管體端部應(yīng)力相對(duì)較高，局部單元已屈服，接箍中面有脹破的趨勢(shì)。對(duì)比圖4a與4b可見，由于該特殊螺紋設(shè)計(jì)的接箍材料強(qiáng)度低于管體，溫度升高后特殊螺紋接箍可能先發(fā)生破壞，加熱升溫過程是伴隨著卸載的，此時(shí)特殊螺紋管端附近應(yīng)力反而增大至784.5 MPa，螺紋牙應(yīng)力下降至300.0 MPa左右。分析其原因，是由于溫度升高，特殊螺紋材料屈服強(qiáng)度降低，變形增大，而在施加拉伸載荷的一側(cè)卸載，螺紋段產(chǎn)生的變形反而傳遞到管端，使密封面與臺(tái)肩面產(chǎn)生壓縮效應(yīng)，應(yīng)力增大。

第1次保載后，等效應(yīng)力分布如圖4c，保載施加拉伸載荷后特殊螺紋整體應(yīng)力水平明顯提高，在第1、16、17、18牙螺紋附近、接箍端部附近都有局部屈服發(fā)生，但特殊螺紋整體變形并不明顯。

第1次拉伸+內(nèi)壓下熱循環(huán)50次后等效應(yīng)力分布如圖4d。與圖4c對(duì)比可見，50次溫度循環(huán)下整體應(yīng)力水平反而降低，管體部分應(yīng)力水平稍有提高，且分布不均。分析可知高、低溫載荷反復(fù)作用50次后，特殊螺紋先前屈服部分可能已經(jīng)發(fā)生組織變化，部分區(qū)域強(qiáng)度已將提高，而管體由于整體尺寸較長(zhǎng)，強(qiáng)度變化區(qū)域也不一致，所以出現(xiàn)應(yīng)力分布不均現(xiàn)象。

卸載后溫度載荷并沒有變化，其應(yīng)力分布如圖4e，特殊螺紋整體應(yīng)力變化降低，只在管體端點(diǎn)附近和接箍中部表面有發(fā)生屈服的趨勢(shì)，特殊螺紋管端附近應(yīng)力增大至755.0 MPa，螺紋牙應(yīng)力下降至250.0 MPa左右，原因類似于加熱升溫過程。高溫下拉伸與內(nèi)壓載荷循環(huán)5次應(yīng)力分布如圖4f，特殊螺紋應(yīng)力分布合理，沒有屈服現(xiàn)象發(fā)生，特殊螺紋整體應(yīng)力水平都有所增大，管端附近稍微減小為748 MPa，這是由于高溫下再次施加載荷后特殊螺紋螺紋段再次承受較大載荷，螺紋牙相互之間的壓縮變形增大，傳遞至管端后密封面和臺(tái)肩面承受的壓縮變形減小，應(yīng)力水平稍有下降。

第2次保載，施加了拉伸載荷后等效應(yīng)力分布如圖4g，與圖4f對(duì)比可見，特殊螺紋應(yīng)力變化比較明顯，只在螺紋段最后幾牙附近產(chǎn)生了局部屈服，分析可知施加軸向載荷對(duì)螺紋連接的影響較大，尤其對(duì)螺紋牙的軸向剪切作用。

第2次拉伸+內(nèi)壓下熱循環(huán)50次后等效應(yīng)力分布如圖4h，與圖4d對(duì)比可見，特殊螺紋應(yīng)力整體變化水平不高，反而第2次在接箍中面的屈服趨勢(shì)更明顯，對(duì)比第1次拉伸+內(nèi)壓的熱循環(huán)50次，發(fā)現(xiàn)管端密封面與臺(tái)肩面應(yīng)力稍有下降，變形仍處于彈性階段，這是由于在多次循環(huán)載荷作用下，特殊螺紋材料處于疲勞損傷狀態(tài)。

冷卻后特殊螺紋的等效應(yīng)力分布如圖4i，冷卻的過程性相當(dāng)于卸掉包括溫度在內(nèi)的所有載荷，特殊螺紋應(yīng)力水平明顯降低。

第2次室溫下載荷循環(huán)5次后等效應(yīng)力分布如圖4k，特殊螺紋應(yīng)力水平明顯下降，各部分布都較合理，僅在接箍中部外表面出現(xiàn)應(yīng)力集中，但不會(huì)發(fā)生斷裂。

表3 套管特殊特殊螺紋C系熱循環(huán)下關(guān)鍵部位等效應(yīng)力值 MPa

a 第1次室溫下載荷循環(huán)5次后等效應(yīng)力分布

b 加熱升溫后

c 第1次保載后

d 第1次拉伸+內(nèi)壓下熱循環(huán)50次后

e 卸載后

f 高溫下載荷循環(huán)5次后

g 第2次保載后

h 第2次拉伸+內(nèi)壓下熱循環(huán)50次后

i 冷卻后

k 第2次室溫下載荷循環(huán)5次后

5 結(jié)論

1) 按照ISO 13679 CALⅢ C系熱循環(huán)試驗(yàn)的實(shí)物評(píng)價(jià)程序和要求，采用有限元方法對(duì)套管特殊螺紋進(jìn)行了計(jì)算機(jī)仿真評(píng)價(jià)。接頭連接的邊界條件復(fù)雜，屬于彈塑性接觸非線性問題，本次分析又考慮了溫度場(chǎng)和應(yīng)力場(chǎng)的耦合作用，分析計(jì)算具有一定的代表性和難度。限于篇幅，本文僅進(jìn)行了熱循環(huán)試驗(yàn)評(píng)價(jià)，未進(jìn)行載荷包絡(luò)線評(píng)價(jià)，但作為一種方法上的探索，足以說明采用有限元仿真評(píng)價(jià)方法按照ISO 13679程序要求進(jìn)行復(fù)雜載荷下特殊螺紋的強(qiáng)度分析是可行的，采用仿真評(píng)價(jià)手段可以大幅地減少特殊螺紋力學(xué)評(píng)價(jià)的經(jīng)濟(jì)成本和時(shí)間。

2) 在第1、17牙螺紋、密封面、臺(tái)肩面等部位可能發(fā)生屈服，整個(gè)過程中應(yīng)力水平變化較大。隨著溫度由低溫變至高溫的變化，管體與接箍的屈服強(qiáng)度下降，部分牙螺紋、接箍中面等塑性變形加大，而局部由于高、低溫反復(fù)作用可能會(huì)產(chǎn)生局部組織的變化，使先前本應(yīng)屈服的區(qū)域強(qiáng)度反而提高，但整體上在C系加載完成后仍能保持結(jié)構(gòu)完整。說明特殊螺紋要想安全通過C系熱循環(huán)試驗(yàn)，有必要改善不完整螺紋段表面的強(qiáng)度，增強(qiáng)接箍所用材料的強(qiáng)度。

3) 套管特殊螺紋在C系熱循環(huán)作用下的高溫狀態(tài)，材料的屈服強(qiáng)度降低，其變形更多的處于彈性變形階段，室溫狀態(tài)反而特殊螺紋的關(guān)鍵部位發(fā)生屈服的可能性更大。但考慮到C系載荷循環(huán)的交替進(jìn)行，實(shí)際上反映在特殊螺紋上的塑性損傷積累作用更加明顯，特殊螺紋實(shí)際的變形也更大。

[1] API RP 5C5，Recommended Practice on Procedures for Testing Casing and Tubing Connections[S]./ISO 13679—2002，Petroleum and natural gas industries——Procedures for testing casing and tubing connections[S].

[2] 王雙來，彭娜，劉卜.高溫高壓井特殊特殊螺紋的選用與評(píng)價(jià)試驗(yàn)[J].鋼管，2016(1)：64-71.

[3] Sugino M，Nakamura K，Yamaguchi S，et al.Development of an Innovative High-performance Premium Threaded Connection for OCTG[C]//Paper OTC 20734 prepared at the 2010 Offshore Technology Conference held in Houston，Texas，USA，3-6 May 2010.

[4] Gabriel E C，Andrea Bufalini，Luis Conde，et al.Low SCF Integral Premium Connections for Use in Highly Demanding Casing and Tubing Drilling Operations[C]//Paper OTC 16566 prepared for presentation at the Offshore Technology Conference held in Houston，Texas，USA.，3-6 May 2004.

[5] Zhu B，Steen T L，Patterson W R，et al.Low Cycle Fatigue Analysis of Threaded Connectors for Subsea Drilling and Production[C]//Paper SPE 166554 prepared for presentation at the SPE Offshore Europe Oil and Gas Conference and Exhibition held in Aberdeen，UK，3-6 September 2013.

[6] 申昭熙，李磊，王鵬.特殊螺紋接頭主要參數(shù)對(duì)密封性能的影響分析[J].石油礦場(chǎng)機(jī)械，2011，40(1)：38-41.

[7] 竇益華，馬亮，李明飛，等.過扭與欠扭工況下特殊螺紋油管接頭三維有限元分析[J].石油礦場(chǎng)機(jī)械，2016，45( 4)：52-56.

[8] 許志倩，閆相禎，楊秀娟.非API標(biāo)準(zhǔn)Big Omega特殊螺紋接頭連接性能數(shù)值分析[J].石油礦場(chǎng)機(jī)械，2009，38(6)：34-39.

[9] 許志倩，閆相禎，楊秀娟，等.微觀泄漏機(jī)理在非API套管特殊螺紋密封性能評(píng)估中的應(yīng)用[J].石油學(xué)報(bào)，2014(5)：963-971.

[10] 段敬黎，祝效華，童華，等.?139.7mm套管偏梯形螺紋接頭的有限元分析[J].石油礦場(chǎng)機(jī)械，2008，37(11)：59-62.

[11] 祝效華，張智，常學(xué)，等.復(fù)雜結(jié)構(gòu)井磨損套管連接螺紋的三維力學(xué)行為[J].石油學(xué)報(bào)，2015(6)：748-753.

[12] 竇益華，王軻，于洋，等.特殊螺紋油管特殊螺紋上螺紋性能三維有限元分析[J].石油機(jī)械，2015(4)：99-104.

[13] 高連新，汪華林，張毅，等.窄間隙特殊螺紋套管的使用性能分析[J].天然氣工業(yè)，2007(11)：58-60.

[14] 江丙云，孔祥宏，羅元元.ABAQUS工程實(shí)例詳解[M].北京：人民郵電出版社，2014.

[15] 莊茁，由小川，廖劍暉，等.基于ABAQUS的有限元分析和應(yīng)用[M].北京：清華大學(xué)出版社，2009：450-466.

StrengthAnalysisofCasingPremiumJointBasedonISO13679ThermalCycleTest

BAI Song1，ZHANG Jianbing2，DU Zhijie3

(1.ShaanxiInstituteofDefenseIndustryProfessionalTechnology，Huxian710300，China；2.Xi’anShiyouUniversity，Xi’an710065，China；3.Xi’anMaurerPetroleumEngineeringLaboratory，Xi’an710065，China)

In this paper，according to the“ISO 13679-casing and tubing joint test procedure recommended practice”the C series thermal cycle test procedure，abaqus finite element analysis software was used to evaluate the strength of a new type of casing special joint JIP ?177.8 mm× 10.36 mm P110.Considering the nonlinearity of the calculation of the contact problem，a finite element model of the special joint of the casing is established.The C-loading path of the joint is simplified，and the thermal stress of the joint under the coupling of the stress field and the temperature field is calculated and analyzed.The results found in this research include：As the temperature changes from low temperature to high temperature，the plastic deformation of the joint part of the thread and the middle of the coupling increases，and the local and high temperature may cause the change of the organization structure due to the high and low temperature.But on the whole，after the completion of the C series thermal cycle test is completed to maintain the integrity of the structure.In the high temperature state，the yield strength of the material decreases，and its deformation is more in the elastic deformation stage.At room temperature，the key part of the special buckle is more likely to yield.It is recognized that proper computer simulation can be applied to effectively evaluate the gas-tight seal connections in accordance with test content and procedures of the ISO 13679，which greatly saves time and economic cost and can replace the evaluation on full-scale test of connections to a certain extent.

casing；premium joint；strength；thermo-mechanical coupling；finite element

1001-3482(2017)06-0040-07

2017-06-01

國(guó)家自然科學(xué)基金項(xiàng)目“膨脹套管螺紋連接幾何、材料與接觸非線性工作特性研究”(51074126)；陜西省教育廳項(xiàng)目“石油膨脹套管螺紋連接結(jié)構(gòu)與密封完整性非線性數(shù)值模擬研究”(11JK0792)

白 松(1991-)，男，碩士，2017年畢業(yè)于西安石油大學(xué)，主要研究方向?yàn)橛?jì)算機(jī)輔助工程，E-mail：howard257@163.com。

TE931.2

A

10.3969/j.issn.1001-3482.2017.06.009