張寶巖， 盧祥國(guó)， 謝　坤， 劉義剛， 張?jiān)茖?/p>

(1.東北石油大學(xué) 提高油氣采收率教育部重點(diǎn)實(shí)驗(yàn)室，黑龍江 大慶 163318；2.中海石油有限公司天津分公司 渤海石油研究院，天津 300452)

表2　采收率實(shí)驗(yàn)數(shù)據(jù)

Table 2　Data recovery experiments

?

Cr3+聚合物凝膠與水交替注入調(diào)驅(qū)效果和機(jī)理分析
——以渤海油田儲(chǔ)層條件為例

張寶巖1， 盧祥國(guó)1， 謝坤1， 劉義剛2， 張?jiān)茖?

(1.東北石油大學(xué) 提高油氣采收率教育部重點(diǎn)實(shí)驗(yàn)室，黑龍江 大慶 163318；2.中海石油有限公司天津分公司 渤海石油研究院，天津 300452)

針對(duì)礦場(chǎng)減緩吸液剖面反轉(zhuǎn)技術(shù)需求，以渤海油藏儲(chǔ)層為模擬對(duì)象，以注入壓力、含水率和采收率為評(píng)價(jià)指標(biāo)，開(kāi)展了“Cr3+聚合物凝膠+水”交替注入調(diào)驅(qū)方式增油效果實(shí)驗(yàn)及礦場(chǎng)研究。結(jié)果表明，調(diào)驅(qū)劑進(jìn)入儲(chǔ)層中低滲透層，一方面擴(kuò)大了波及體積，另一方面增加了滲流阻力和吸液?jiǎn)?dòng)壓力，引起吸液剖面反轉(zhuǎn)。采用“Cr3+聚合物凝膠+水”交替注入調(diào)驅(qū)方式，不僅有利于增強(qiáng)Cr3+聚合物凝膠前置段塞對(duì)高滲透層的封堵作用，還可進(jìn)一步發(fā)揮后續(xù)水段塞轉(zhuǎn)向進(jìn)入中低滲透層后的驅(qū)油作用，從而減緩甚至消除調(diào)驅(qū)劑進(jìn)入儲(chǔ)層中低滲透層后引起的吸液剖面反轉(zhuǎn)現(xiàn)象。LD5-2油田A22井實(shí)施“Cr3+聚合物凝膠+水”交替注入調(diào)驅(qū)工藝后，壓力明顯提高，有利于增強(qiáng)驅(qū)替劑擴(kuò)大波及體積效果。

渤海油藏；聚合物凝膠；交替注入；增油效果；物理模擬；機(jī)理分析

近年來(lái)，在陸地油田開(kāi)發(fā)陸續(xù)進(jìn)入中高含水期以及新增探明儲(chǔ)量逐年遞減的情況下，海上油田的勘探開(kāi)發(fā)工作受到了油田科技工作者的高度重視，其中渤海油田原油儲(chǔ)量豐富，但受儲(chǔ)層非均質(zhì)性強(qiáng)、巖石膠結(jié)疏松、原油黏度高和注入水礦化度高等因素影響，水驅(qū)采收率較低，亟待進(jìn)一步提高原油采收率[1-4]。為此，渤海地區(qū)LD5-2、SZ36-1和LD10-1等普通稠油油田先后進(jìn)行了聚合物(HPAM)驅(qū)和聚合物凝膠調(diào)驅(qū)礦場(chǎng)試驗(yàn)，并取得了良好的增油效果[5-8]。但分析礦場(chǎng)實(shí)際生產(chǎn)數(shù)據(jù)發(fā)現(xiàn)，化學(xué)驅(qū)實(shí)際增油效果與預(yù)期效果仍存在一定差距，分析原因認(rèn)為，海上油田目前注入工藝為連續(xù)注入，一旦調(diào)驅(qū)劑進(jìn)入并滯留于中低滲透層中造成污染后，調(diào)驅(qū)劑滲流阻力在中低滲透層中的增加速度遠(yuǎn)大于高滲透層，造成中低滲透層的啟動(dòng)壓力升高，儲(chǔ)層吸液剖面出現(xiàn)反轉(zhuǎn)，進(jìn)而導(dǎo)致后續(xù)注入過(guò)程中調(diào)驅(qū)劑主要在剩余油飽和度較低的高滲透層中低效乃至無(wú)效循環(huán)，無(wú)法進(jìn)一步擴(kuò)大波及體積。現(xiàn)有研究表明，大慶油田利用交替注入調(diào)驅(qū)方式在一定程度上減緩了吸液剖面反轉(zhuǎn)現(xiàn)象[9-10]。為此，本文在結(jié)合Cr3+聚合物凝膠良好封堵性和單純水溶液對(duì)儲(chǔ)層無(wú)污染的基礎(chǔ)上，開(kāi)展了“Cr3+聚合物凝膠+水”交替注入調(diào)驅(qū)方式增油效果實(shí)驗(yàn)研究和機(jī)理分析，這將為渤海油田礦場(chǎng)提高原油采收率實(shí)踐技術(shù)決策提供新思路。

1　實(shí)驗(yàn)部分

1.1實(shí)驗(yàn)材料

聚合物為速溶型部分水解聚丙烯酰胺，有效含量為88%，勝利油田生產(chǎn)。

交聯(lián)劑為醋酸鉻溶液，鉻離子有效含量為2.7%。

實(shí)驗(yàn)用水為SZ36-1油田模擬注入水，總礦化度為9 047.6 mg/L，具體離子組成見(jiàn)表1。

表1　溶劑水離子組成

實(shí)驗(yàn)用油為SZ36-1油田模擬油，由渤海某油田脫氣原油與煤油混合配制而成，65 ℃時(shí)黏度為75 mPa·s。

實(shí)驗(yàn)巖心為石英砂環(huán)氧樹(shù)脂膠結(jié)而成的層內(nèi)非均質(zhì)人造巖心[11]，高、中、低滲透層的氣測(cè)滲透率分別為3 600×10-3、720×10-3、180×10-3μm2，平均滲透率1 500×10-3μm2。巖心外觀尺寸為：寬×高×長(zhǎng)=4.5 cm×4.5 cm×30 cm。

1.2實(shí)驗(yàn)儀器

采用DV-Ⅱ型布氏黏度計(jì)(美國(guó)Brookfield公司)測(cè)試調(diào)驅(qū)劑黏度，轉(zhuǎn)速為6 r/min；采用BI-200SM型廣角動(dòng)/靜態(tài)光散射儀系統(tǒng)(Brookhaven Instruments Cop，USA)測(cè)試聚合物分子線團(tuán)尺寸Dh；采用驅(qū)替實(shí)驗(yàn)裝置測(cè)試調(diào)驅(qū)劑增油降水效果，裝置由壓力傳感器、中間容器手搖泵、平流泵和巖心夾持器等組成，測(cè)試溫度為65 ℃[7]。

1.3方案設(shè)計(jì)

實(shí)驗(yàn)方案步驟均為“水驅(qū)95%+Cr3+聚合物凝膠整體或其與水交替注入+后續(xù)水驅(qū)95%”?！胺桨?-0”為整體Cr3+聚合物凝膠段塞調(diào)驅(qū)實(shí)驗(yàn)(0.1 PV，Cp=2 000 mg/L，聚合物和Cr3+的質(zhì)量比(簡(jiǎn)稱聚鉻比) =180∶1)，其為后續(xù)交替注入調(diào)驅(qū)實(shí)驗(yàn)方案增油效果的對(duì)比基礎(chǔ)?！胺桨?-1”～“方案1-4”為Cr3+聚合物凝膠與水交替注入調(diào)驅(qū)實(shí)驗(yàn)，其中Cr3+聚合物凝膠總段塞尺寸為0.1 PV，中間穿插水段塞總尺寸也為0.1 PV，依據(jù)交替注入輪次將Cr3+聚合物凝膠和水段塞尺寸進(jìn)行等分后，開(kāi)展“Cr3+聚合物凝膠+水”交替注入調(diào)驅(qū)實(shí)驗(yàn)。

2　結(jié)果分析

2.1采收率

不同交替注入輪次條件下，Cr3+聚合物凝膠與水交替注入調(diào)驅(qū)實(shí)驗(yàn)采收率數(shù)據(jù)見(jiàn)表2。

從表2中可以看出，在Cr3+聚合物凝膠和水段塞尺寸一定條件下，隨交替注入輪次數(shù)增加，Cr3+聚合物凝膠調(diào)驅(qū)采收率增幅提高。

2.2動(dòng)態(tài)特征

實(shí)驗(yàn)過(guò)程中注入壓力、含水率和采收率與注入PV數(shù)關(guān)系對(duì)比結(jié)果見(jiàn)圖1。

表2采收率實(shí)驗(yàn)數(shù)據(jù)

Table 2Data recovery experiments

方案編號(hào)交替注入輪次凝膠黏度/(mPa·s)含油飽和度/%采收率/%水驅(qū)最終增幅較整體段塞采收率增幅/%1-0整體16.573.526.234.78.5—1-1216.573.426.537.010.52.01-2316.573.226.338.912.64.11-3416.573.426.139.813.75.21-4516.573.526.040.914.96.4

圖1　注入過(guò)程動(dòng)態(tài)特征曲線

從圖1可以看出，隨Cr3+聚合物凝膠與水交替注入次數(shù)增加，注入壓力增加，擴(kuò)大波及體積效果增強(qiáng)，含水率降幅增大，采收率增幅提高。

2.3機(jī)理分析

Cr3+聚合物凝膠與水交替注入改善調(diào)驅(qū)效果機(jī)理分析所用地質(zhì)模型見(jiàn)圖2，模型包括高、中、低3個(gè)滲透層。

圖2　典型油藏地質(zhì)模型

從圖2可以看出，受模型自身非均質(zhì)性影響，水驅(qū)過(guò)程中高滲透層因啟動(dòng)壓力(p1)較低，吸液壓差(p-p1)較大，吸液量較多，中低滲透層因吸液壓差較小造成吸液量較少。隨著水驅(qū)過(guò)程持續(xù)進(jìn)行，高滲透層采出程度增加，水相過(guò)流通道斷面增加，滲流阻力減小，啟動(dòng)壓力(p1)降低，這又進(jìn)一步促使吸液壓差(p-p1)增加，吸液量增多。與此同時(shí)，中低滲透層吸液量逐漸降低，波及效果變差，最終影響水驅(qū)開(kāi)發(fā)效果[12]。在實(shí)施Cr3+聚合物凝膠調(diào)驅(qū)初期，調(diào)驅(qū)劑會(huì)首先進(jìn)入滲流阻力較小的高滲透層，并于其中發(fā)生滯留，造成流動(dòng)通道過(guò)流斷面減小，注入壓力(p)升高。隨注入壓力升高，中低滲透層吸液壓差(p-p2和p-p3)增加，吸液量增大，中低滲透層動(dòng)用程度增加。然而，進(jìn)入并滯留于中低滲透層中的Cr3+聚合物凝膠也會(huì)造成流動(dòng)通道過(guò)流斷面減小和滲流阻力增加，并且滲流阻力增加速度要遠(yuǎn)大于高滲透層，這便會(huì)造成中低滲透層尤其是低滲透層啟動(dòng)壓力(p2和p3)升高，吸液壓差(p-p2和p-p3)減小，吸液量降低，出現(xiàn)“吸液剖面反轉(zhuǎn)”現(xiàn)象。當(dāng)采用“Cr3+聚合物凝膠+水”交替注入方式時(shí)，可以產(chǎn)生“封堵和驅(qū)替”效應(yīng)，即凝膠前置段塞對(duì)高滲透層實(shí)施封堵，后續(xù)水段塞進(jìn)入中低滲透層驅(qū)油，減緩甚至消除了因凝膠進(jìn)入中低滲透層而引起的“吸液剖面反轉(zhuǎn)”現(xiàn)象。因此，只要交替注入段塞尺寸和輪次組合與儲(chǔ)層非均質(zhì)性相匹配，就可以最大限度地發(fā)揮“封堵和驅(qū)替”協(xié)同效應(yīng)，提高調(diào)驅(qū)增油降水效果。

進(jìn)一步分析發(fā)現(xiàn)，隨交替注入輪次增加，注入壓力增大(見(jiàn)圖1)，這與Cr3+聚合物凝膠分子線團(tuán)(聚集體)遇水膨脹特性有關(guān)。不同稀釋質(zhì)量濃度條件下，Cr3+聚合物凝膠及聚合物溶液與注入水混合稀釋前后分子線團(tuán)尺寸Dh測(cè)試結(jié)果見(jiàn)表3。

表3　Dh測(cè)試結(jié)果

由表3可以看出，隨著稀釋質(zhì)量濃度降低，Cr3+聚合物凝膠中分子線團(tuán)Dh逐漸增加，聚合物溶液中分子線團(tuán)Dh逐漸減小。在“Cr3+聚合物凝膠+水”交替注入換注水段塞時(shí)，巖心孔隙內(nèi)陽(yáng)離子質(zhì)量濃度受稀釋作用而降低，聚合物凝膠分子線團(tuán)表面原有電荷動(dòng)態(tài)平衡被打破[13-14]，聚合物分子鏈上離子基團(tuán)所帶負(fù)電荷數(shù)量增加，彼此之間的排斥力增大，進(jìn)而造成卷曲的分子鏈趨于舒展，分子線團(tuán)尺寸增大。隨聚合物凝膠分子線團(tuán)尺寸增大，一方面，在巖心孔隙中已發(fā)生捕集的聚合物分子聚集體的封堵作用得以加強(qiáng)。另一方面，部分原本處于自由運(yùn)移狀態(tài)的聚合物分子聚集體會(huì)因尺寸增大而發(fā)生捕集，從而造成孔道過(guò)流斷面進(jìn)一步減小，滲流阻力增大，注入壓力提高。

3　礦場(chǎng)試驗(yàn)效果分析

3.1試驗(yàn)區(qū)概況

LD5-2油田位于渤海遼東灣海域，發(fā)現(xiàn)于1988年6月，油田主要包括東二上、下段兩個(gè)區(qū)塊。東二上段儲(chǔ)層物性較好，孔隙度主要分布在32%～40%，滲透率一般大于1 000 mD。東二下段儲(chǔ)層孔隙度主要分布在30%～36%，滲透率分布在10～1 320 mD，主要集中在100～1 000 mD。A22井作為2號(hào)塊東二上段設(shè)計(jì)注水井，先期進(jìn)行排液生產(chǎn)，為了及時(shí)補(bǔ)充地層能量，2007年5月26日將A22井轉(zhuǎn)注，初期日注水量280 m3/d左右。A22井的注水受效井為A15、A16、A17、A23、A27、A26、A21s井。

3.2注入工藝方案設(shè)計(jì)

Cr3+聚合物凝膠與水交替注入調(diào)驅(qū)礦場(chǎng)試驗(yàn)注入工藝設(shè)計(jì)方案見(jiàn)表4。

表4　A22井注入工藝方案

3.3試驗(yàn)效果

3.3.1注入井注入壓力LD5-2油田A22井由2014年12月28日開(kāi)始按照注入施工方案進(jìn)行調(diào)驅(qū)作業(yè)，由其注入過(guò)程動(dòng)態(tài)特征曲線(見(jiàn)圖3)可見(jiàn)，實(shí)施“Cr3+聚合物凝膠+水”交替注入調(diào)驅(qū)工藝后，生產(chǎn)壓力明顯提高，有利于增強(qiáng)驅(qū)替劑擴(kuò)大波及體積效果。

3.3.2油井含水率變化和增油效果圖4為試驗(yàn)區(qū)受效井A17井綜合曲線。從圖4可以看出，調(diào)驅(qū)作業(yè)后A17井含水率上升趨勢(shì)得到了有效抑制。目前，LD5-2油田A22井調(diào)驅(qū)施工仍在進(jìn)行之中，由于作業(yè)時(shí)間較短，油井還未見(jiàn)到明顯增油降水效果。

4　結(jié)論

(1) 調(diào)驅(qū)劑進(jìn)入非均質(zhì)儲(chǔ)層中低滲透層，一方面擴(kuò)大了波及體積，另一方面增加了滲流阻力和吸液?jiǎn)?dòng)壓力，這是造成吸液剖面反轉(zhuǎn)的根本原因。

(2) 采用“Cr3+聚合物凝膠+水”交替注入方式時(shí)，可以發(fā)揮凝膠前置段塞對(duì)高滲透層封堵和后續(xù)水段塞對(duì)中低滲透層驅(qū)油作用，減緩甚至消除了因凝膠進(jìn)入中低滲透層而引起的“吸液剖面反轉(zhuǎn)”現(xiàn)象。

(3) 在“Cr3+聚合物凝膠+水”交替注入調(diào)驅(qū)過(guò)程中，水段塞可稀釋巖心孔隙中陽(yáng)離子濃度，引起聚合物分子線團(tuán)尺寸增加，滲流阻力增大，注入壓力提高。

(4) LD5-2油田A22井實(shí)施“Cr3+聚合物凝膠+水”交替注入調(diào)驅(qū)工藝后，壓力明顯提高，這十分有利于擴(kuò)大波及體積效果。

圖3　A22井調(diào)驅(qū)動(dòng)態(tài)特征曲線

圖4　A17井綜合曲線

[1]王君, 范毅. 稠油油藏的開(kāi)采技術(shù)和方法[J]. 西部探礦工程, 2006(7): 84-85.

Wang Jun, Fan Yi. The distribution of hydrogen sulfide is controlled by sedim entary facies[J]. West-China Exploration Engineering, 2006(7): 84-85.

[2]張鳳久, 羅憲波, 劉英憲, 等. 海上油田叢式井網(wǎng)整體加密調(diào)整技術(shù)研究[J]. 中國(guó)工程科學(xué), 2011, 13(5): 34-40.

Zhang Fengjiu, Luo Xianbo, Liu Yingxian, et al. Research on overall encryption adjustmenttechnology of offshore oil field[J]. Engineering Sciences,2011, 13(5): 34-40.

[3]周守為, 韓明. 用于海上油田化學(xué)驅(qū)的聚合物研究[J]. 中國(guó)海上油氣, 2007, 19(1): 1-2.

Zhou Shouwei,Han Ming. Study on polymer for chemical floodingin offshore oilfield of China[J]. China Offshore Oil and Gas,2007, 19(1): 1-2.

[4]周守為, 韓明, 向問(wèn)陶, 等. 渤海油田聚合物驅(qū)提高采收率技術(shù)研究及應(yīng)用[J]. 中國(guó)海上油氣, 2006, 18(6): 386-389.

Zhou Shouwei, Han Ming, Xiang Wentao, et al. Application of EOR technology by means ofpolymer flooding in Bohai oilfields[J]. China Offshore Oil and Gas, 2006, 18(6): 386-389.

[5]姜維東. Cr3+聚合物凝膠性能及其應(yīng)用效果研究[D]. 大慶：大慶石油學(xué)院, 2009.

[6]陳欣, 盧祥國(guó), 鞠野, 等. 速溶聚合物及其Cr3+凝膠滲流特性[J]. 油田化學(xué), 2014, 31(3): 390-394.

Chen Xin, Lu Xiangguo, Ju Ye, et al. SolubleInstant polymer and its Cr3+gel seepage characteristics[J].Oilfield Chemistry, 2014, 31(3): 390-394.

[7]劉進(jìn)祥, 盧祥國(guó), 李國(guó)橋, 等. 特高礦化度Cr3+交聯(lián)聚合物溶液滲流特性及其機(jī)制[J]. 中國(guó)石油大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013, 37(6): 145-152.

Liu Jinxiang, Lu Xiangguo, Li Guoqiao, et al. Seepage characteristics and mechanism of Cr3+cross-linked polymer solution in ultra-high salinity media[J]. Journal of China University of Petroleum(Edition of Natural Science), 2013, 37(6): 145-152.

[8]鄒劍, 王榮健, 薛寶慶, 等. BZ28-2S油田Cr3+聚合物凝膠滲流特性及其調(diào)驅(qū)效果研究[J]. 海洋石油, 2014, 34(5): 46-53.

Zou Jian,Wang Rongjian,Xue Baoqing, et al. Study on the percolation property of Cr3+polymer gel and profile controlling effect in BZ28-2S oilfield [J].Offshore Oil, 2014, 34(5): 46-53.

[9]陳廣宇, 田燕春, 趙新, 等. 大慶油田二類油層復(fù)合驅(qū)注入方式優(yōu)化[J]. 石油學(xué)報(bào), 2012, 33(3): 459-464.

Chen Guangyu, Tian Yanchun, Zhao Xin, et al.Optimization of the ASP flooding injection pattern for sub-layers in Daqing oilfield[J].Acta Petrolei Sinica, 2012，33(3)：459-464.

[10]侯維虹. 聚合物驅(qū)油層吸水剖面變化規(guī)律[J]. 石油勘探與開(kāi)發(fā), 2007, 34(4): 478-482.

Hou Jihong. Water injection profile variation of oil layers under polymer flooding[J].Petroleum Exploration and Development, 2007，34(4)：478-482.

[11]盧祥國(guó), 高振環(huán), 閆文華, 等. 人造巖心制作及其影響因素實(shí)驗(yàn)研究[J]. 大慶石油地質(zhì)與開(kāi)發(fā), 1994, 3(4): 53-55.

Lu Xiangguo,Gao Zhenhuan,Yan Wenhua, et al. [J].Petroleum Geology & Oilfield Development in Daqing, 1994, 3(4): 53-55.

[12]盧祥國(guó), 王樹(shù)霞, 王榮健, 等. 深部液流轉(zhuǎn)向劑與油藏適應(yīng)性研究—以大慶喇嘛甸油田為例[J]. 石油勘探與開(kāi)發(fā)，2011, 38(5): 576-582.

Lu Xiangguo, Wang Shuxia, Wang Rongjian, et al. Adaptability of a deep profile control agent to reservoirs: Taking the Lamadian Oilfield in Daqing as an example[J]. Petroleum Exploration and Development，2011, 38(5): 576-582.

[13]符若文，李谷，馮開(kāi)才. 高分子物理[M]. 1版.北京：化學(xué)工業(yè)出版社，2004.

[14]陳震，趙孔雙，劉昊. 非導(dǎo)電球型粒子懸浮液寬頻介電譜—薄雙電層理論的計(jì)算機(jī)模擬[J]. 化學(xué)學(xué)報(bào)，2006，64(17)：1780-1784.

Chen Zhen，Zhao Kongshuang，Liu Hao. Broadband dielectric spectroscopy of non-conducting spherical particle suspensions:Computer simulation of thin double layer theory[J]. Acta Chimica Sinic, 2006，64(17)：1780-1784.

(編輯宋官龍)

Profile Control and Flooding Effect of Alternate Injection of Cr3+Polymer Gel and Water and Its Mechanism Analysis:Taking the Bohai Oilfield as an Example

Zhang Baoyan1, Lu Xiangguo1, Xie Kun1, Liu Yigang2, Zhang Yunbao2

(1.KeyLaboratoryofEnhancedOilandGasRecoveryofEducationMinistry，NortheastPetroleumUniversity,DaqingHeilongjiang163318,China;2.BohaiOilfieldReserchInstitute，ChinaNationalOffshoreOilfieldCorporationLimite-Tianjin,Tianjin300452,China)

Aiming at technical demand on relieving the reverse of imbibition profile, displacement effect of alternate injection of Cr3+polymer gel and water was studied, taking reservoir of Bohai oilfield as simulation object, regarding injection pressure, water content and recovery efficiency as evaluation index. Results showed that as the displacement agent flowed into the medium and low permeability layer, on the one hand, the swept volume was magnified. On the other hand, seepage resistance and start-up pressure of imbibition were increased, thus resulting in the reverse of the imbibition profile. Once the alternate injection of Cr3+polymer gel and water was adopted, polymer solution could get priority to flow into the high permeable layer and block off it, then subsequent water flowed into medium and low permeability layers to displace the oil. Therefore the phenomenon of the reverse of imbibition profile was relieved or even eliminated. Once the alternate injection of Cr3+polymer gel and water was conducted in A22 injection well in LD5-2 Oilfield, injection pressure get increased, which was good for enlarging sweep volume effect.

Bohai oilfield; Polymer gel; Alternate injection; Displacement effect; Physical simulation; Mechanism analysis

1006-396X(2016)01-0035-06

投稿網(wǎng)址：http://journal.lnpu.edu.cn

2015-08-27

2015-10-22

國(guó)家“十二五”油氣重大專項(xiàng)子課題“海上稠油保壓熱采技術(shù)示范”(2011ZX05057-005-003)。

張寶巖(1990-)，男，碩士研究生，從事提高采收率理論和技術(shù)方面的研究；E-mail: 369499907@qq.com。

盧祥國(guó)(1960-)，男，博士，教授，博士生導(dǎo)師，從事提高采收率理論和技術(shù)方面的研究； E-mail:luxiangg2003@aliyun.com。

TE357.46

Adoi:10.3969/j.issn.1006-396X.2016.01.007