Simulation of Transformation from Water-Injection Huff and Puff to Unstable Water-Flooding in Developing Fractured Tight Reservoirs
-
摘要: 为解决裂缝性致密油藏经多轮次吞吐后,产量降低过快的问题,以天然裂缝较为发育的某致密油藏M区块为例,综合考虑基质、天然裂缝和压裂裂缝的物性和压力差异,开展了数值模拟研究,分析了裂缝尖端应力场和裂缝扩展特征。在此基础上,对比分析了注水吞吐、不稳定水驱的开发效果。结果表明,地层压力随注水时间延长逐渐升高,当地层压力高于裂缝开启压力时,天然裂缝发生扩展并沟通压裂裂缝,形成复杂动态缝网。不稳定水驱开发能够充分发挥渗吸和驱替作用,改变注水量能够有效防止水窜,形成相对均匀的驱替前缘。模拟结果表明,生产时间相同时,采用周期注水时的累计产油量较采用注水吞吐方式约高18%,开发效果显著改善。将开发方式转变为不稳定水驱,能够有效提高裂缝性致密油藏水平井产量,为裂缝性致密油藏水平井开发方式优化提供理论依据。Abstract: Multiple rounds of huff and puff in tight reservoirs usually lead to a rapid decrease in production. Taking the M block of a tight reservoir with developed natural fractures as an example, a numerical simulation was conducted based on the physical properties of the matrix, natural fractures, and fractured fractures as well as the pressure difference to analyze the stress field of fracture tips and the features of fracture propagation. On this basis, a comparative analysis was then carried out to evaluate the development effect of water-injection huff and puff and unstable cyclic water injection. Results show that the formation pressure would grow with an increase in water injection time, and when the formation pressure was higher than the opening pressure of fractures, a complex and dynamic fracture network was formed with the expansion of natural fractures and the communication of fractured fractures. Unstable water-flooding can give full play to imbibition and displacement, and the change of water injection volume can effectively avoid water channeling and form relatively uniform flooding front. In addition, simulation results show that a significant increase of 18% in cumulative oil production of reservoirs with cyclic water injection compared with water-injection huff and puff. Therefore, transforming the development method into unstable water-flooding can effectively improve the oil production of horizontal wells in fractured tight reservoirs, providing a theoretical reference for optimizing the development of horizontal wells in tight reservoirs.
-
PDC钻头是石油钻探中最常用的钻头类型,中美两国近年85%以上的钻井进尺都是由PDC钻头完成的[1]。由于PDC钻头主要以切削形式破碎井底岩石,钻进时易出现粘滑振动,严重影响钻头工作寿命和破岩效率[2–4]。特别是钻遇高强度地层、砾石地层时,PDC钻头的粘滑振动会更加强烈[5]。粘滑振动还可能与其他方向的振动耦合,形成更加复杂、更具有破坏性的振动形式,对MWD等井下仪器的安全和稳定工作形成巨大威胁,易造成钻具过早疲劳、仪器失效、定向失败等复杂情况,同时还会引起钻井参数难以强化,严重影响钻井效率[6–8]。为提高PDC钻头在深部硬岩地层的切削深度并同时降低钻头的粘滑振动,从而获得更高的破岩效率与钻进稳定性,近年来国内研究者提出了联合轴向冲击和扭转冲击的复合冲击钻进新方法[9–11],并相继研制出轴向与扭转复合冲击钻井工具,在现场钻井减振方面取得了较好的应用效果[12–15]。因此,笔者针对上述PDC钻头破岩过程中的粘滑振动问题,利用室内试验,分析了钻头粘滑振动发生机理以及各参数对其振动特征的影响规律,探索了复合冲击钻井技术的减振提速机理,以期为复合冲击钻井工具的优化设计和应用选型提供理论和技术参考。
1. PDC钻头冲击破岩试验系统
钻柱系统运动极为复杂,通过现场试验测量研究钻柱运动和钻头破岩过程的成本高、影响因素多,分析难度也非常大[16–17]。通过室内试验装置模拟钻柱加载及钻头破岩过程,可以针对性开展钻头振动规律及特征研究,更加直观地了解钻头在不同加载模式下的运动状态[18],有利于指导钻头破岩理论建模和减振工具的研制。笔者利用PDC钻头冲击破岩试验装置开展了试验研究。
1.1 加载系统
PDC钻头冲击破岩试验装置主要包含加载系统、测量系统和冲击载荷发生装置,其中加载系统主要包括轴向加载系统和旋转加载系统[15],如图1所示。轴向加载系统通过试验架顶部的伺服电机驱动升降机上下运动,来模拟实际钻进中的下钻和起钻过程,并可通过设定伺服电机的转速来控制加载或卸载速度,其对钻头加载的载荷主要来自配重钻铤。钻头接触岩石前,系统的轴向负载(配重钻铤)全部由升降机承载;当钻头与岩石发生接触后,升降机负载则会逐渐降低。钻压、配重钻铤重量和升降机负载的关系为配重钻铤重量等于钻压和升降机负载的和。为模拟实际钻井过程中钻柱的轴向刚度,升降丝杠与提升机架之间通过刚度调整机构相连,利用碟簧组刚度模拟长钻柱的轴向刚度。旋转加载系统由旋转驱动电机提供动力,并通过减速机、扭转弹性杆、换向器和链轮等机构将旋转动力传递至钻头。该加载系统通过扭转弹性杆来模拟实际钻柱系统的扭转刚度,并通过链轮结构上的扭转惯性配重来模拟钻柱的扭转惯性。
1.2 测量系统
测量系统各传感器的安装方式及位置如图2所示。动态钻压传感器安装在岩心槽底部,测量钻头钻进过程中对岩石形成的压持力。动态扭矩传感器则安装在岩心槽延伸杆与固定挡块之间,通过测量岩心槽延伸杆处拉压力来计算钻头的破岩扭矩。
1.3 冲击载荷发生装置
PDC钻头冲击破岩试验装置包含轴向冲击和扭转冲击2个冲击模块,安装方式如图3所示(图3中:ω为角速度,rad/s;fai为轴向冲击载荷,N;qai,qti分别为轴向冲击载荷和扭转冲击扭矩的比例系数;
¯W(0)b 为无冲击载荷钻进时的平均钻压,N;¯M(0)t 为无冲击载荷钻进时的平均破岩扭矩,N·m;Mti为扭转冲击扭矩,N·m;f为冲击载荷的频率,Hz;t为时间,s;fti为扭转激励电机输出冲击载荷,N;Lti为扭转激励电机作用点与钻头中心的距离,m)。钻进过程中,2个冲击模块的激励源可以同时产生冲击载荷,实现轴向冲击、扭转冲击和轴扭复合冲击3种冲击条件下的破岩钻进试验。在设定2个方向冲击幅值时,首先以无冲击载荷钻进时的平均钻压
¯W(0)b 和平均破岩扭矩¯M(0)t 为参考,通过设定轴向冲击载荷幅值和平均钻压比例系数qai确定轴向激励电机的输出冲击载荷幅值,同样通过设定扭转冲击扭矩幅值与平均破岩扭矩的比例系数qti确定扭转冲击激励电机的输出冲击扭矩幅值[19],并设定载荷输出方式为正弦波形。计算公式如下:{fai=qai¯W(0)bsin(2πft)Mti=qti¯M(0)tsin(2πft) (1) 冲击扭矩是通过调节激励电机输出载荷实现的,可通过下式确定:
{fti=qti¯F(0)tsin(2πft)Mti=ftiLti (2) 其中¯F(0)t=¯M(0)tLti (3) 2. 试验方案
为了解PDC钻头破岩钻进过程中振动的形式及发展规律,利用PDC钻头冲击破岩试验系统开展了室内实钻试验,分析了钻压、转速、冲击载荷形式及幅频特性等因素对PDC钻头破岩振动的影响。试验用岩石采用尺寸130.0 mm × 130.0 mm × 100.0 mm立方体花岗岩岩样,其弹性模量为38.83 GPa,单轴抗压强度为157.98 MPa。岩样的实测单轴压缩应力−应变曲线如图4所示。
试验采用ϕ56.0 mm矿用PDC钻头,复合片直径13.0 mm,外缘布置3个切削齿,内槽布置2个切削齿,共布置5个切削齿(见图5),钻孔深度85.0 mm,采用清水代替钻井液冷却钻头和携带岩屑。
首先,利用PDC钻头冲击破岩试验系统进行了PDC钻头粘滑振动试验,通过控制送钻速度模拟钻进过程中钻压逐渐升高导致钻头粘滑振动的发展过程,并分析钻头转速的幅频,明确常规无冲击载荷作用下钻头粘滑的幅频特征。在此基础上,进行复合冲击条件下的钻头破岩试验,分析不同频率的轴扭复合冲击载荷对钻头粘滑振动的影响。最后,对比引入复合冲击载荷前后PDC钻头粘滑振动的幅频特征,分析不同频率复合冲击载荷对钻头粘滑振动影响的内在规律。
3. PDC钻头粘滑振动试验
旋转驱动电机转速为20 r/min时,进行了不同送钻速度下的破岩试验,实际钻进参数见表1,其中,ωr为旋转驱动电机转速,ωa为轴向驱动电机送钻参数。由于ωa与轴向驱动电机转速具有正比关系,因此用其表征送钻速度。
表 1 不同送钻参数下的实际钻进参数Table 1. Actual drilling parameters under different drilling parametersωa ωr/(r·min−1) 平均钻压/kN 平均扭矩/(N·m) 平均钻速/(mm·s−1) 10 20 9.29 102.29 0.14 20 20 12.73 164.19 0.28 30 20 15.16 221.37 0.42 40 20 20.12 260.78 0.63 由表1可知,随着平均钻压增大,平均钻速呈线性增加趋势,送钻速度加快了3倍,平均钻压提高了1.17倍,平均钻速提高了3.5倍,平均扭矩提高了1.55倍,平均扭矩的增大速度随着送钻速度增快有下降趋势。
斯伦贝谢公司利用粘滑比来表征粘滑振动等级,粘滑比的表达式为:
rss=ωmax−ωmin2ωr (4) 式中:rss为粘滑比;ωmax为发生粘滑时钻头的最高转速,r/min;ωmin为发生粘滑时钻头的最低转速,r/min。
不同送钻参数下钻头转速和破岩扭矩曲线如图6所示。由表1和图6可知:当送钻参数为10和20时,对应钻压分别为9.29和12.73 kN,粘滑比分别为0.22和0.48,均属于低级别粘滑振动;当送钻参数为30时,钻压为15.16 kN,粘滑比为0.97,为中级别粘滑振动;当送进参数达到40时,钻压达到20.12 kN,此时粘滑比最大为1.31,属于高级别粘滑振动,此时开始出现完全停滞现象。
对图6中钻头转速进行幅频分析,结果如图7所示。从图7可以看出:钻头转速的主振频率在0.6~0.7 Hz;送钻参数为10时,对应的主振振幅最小,为1.27 r/min;主振振幅随着送钻参数增大而升高,且升高幅度逐渐变大,送钻参数为40时对应的主振振幅为11.51 r/min。
由图1可知,破岩试验装置采用长细弹性杆模拟钻柱系统几何特征并传递破岩扭矩,同时在近钻头段安装了多个直径较大的铁盘,以保证模拟钻柱的扭转惯性,试验装置的扭矩传递部分可作为扭杆系统进行考虑,因此该装置扭转振动固有频率的计算公式为[20]:
f0=12π√kθ Im (5) 其中kθ =GJl (6) 式中:f0为扭杆系统固有频率,Hz;kθ为扭转刚度,N·m/rad;G为剪切弹性模量,Pa;J为极惯性矩,m4;l为扭转杆长度,m;Im为转动惯量,kg·m2。
破岩试验装置扭转弹性杆的长度为8 m,扭转刚度为354.43 N·m/rad,扭转配重的转动惯量为13.8 kg·m2,将这些参数代入式(4)求得破岩试验装置的扭转固有频率约为0.81 Hz。由于估算时仅考虑了扭转杆的扭转刚度和扭转惯性配重的转动惯量,忽略了减速器、换向器及其他零部件的扭转惯性,因此破岩试验装置扭转固有频率的估算值偏高。可以认为,粘滑振动的频率0.6~0.7 Hz为破岩试验装置扭转振动的一阶固有频率。
4. 轴扭复合冲击对PDC钻头粘滑的影响试验
轴扭复合冲击对PDC钻头破岩过程具有显著影响,是目前钻井提速的重要技术之一。为了解其对PDC钻头破岩过程中钻头粘滑的影响,在PDC钻头粘滑振动试验基础上,以驱动转速20 r/min、轴向送钻参数30为例,引入复合冲击载荷,进一步开展试验。试验采用频率50,100,150和200 Hz的轴向冲击载荷和扭转冲击载荷进行PDC钻头破岩试验。以无冲击钻进时的平均钻压和平均扭矩为参考,考虑冲击电机激振力量程,将轴向激励电机输出载荷幅值比例系数qai和扭转激励电机输出扭矩幅值比例系数qti均设定为10%。试验组参数设置及试验结果见表2。
表 2 冲击参数设置与实钻参数Table 2. Impact parameter settings and actual drilling parameters试验
编号fai幅值/
kNMti幅值/
(N·m)冲击
频率/Hz平均
钻压/kN平均扭
矩/(N·m)平均钻速/
(mm·s−1)无冲击 0 0 0 15.16 221.37 0.42 CC−1 1.5 22 50 13.87 165.42 0.42 CC−2 1.5 22 100 14.03 182.26 0.43 CC−3 1.5 22 150 14.16 182.68 0.44 CC−4 1.5 22 200 14.11 186.08 0.42 不同频率轴扭复合冲击钻进时测得的转速如图8所示,与图6相比,增加轴扭冲击载荷后,冲击频率在50~100 Hz时,钻头转速的波动明显降低,但随着冲击频率升高,钻头转速波动幅度呈现增大的趋势。
对图8中的钻头转速进行幅频分析,结果见图9。由图9(a)可知,由于轴扭复合冲击载荷的引入,钻头转速在各频率处均存在幅频特征,其中冲击频率50和100 Hz对应的幅频特征最明显,而冲击频率150和200 Hz对应的幅值较小。由图9(b)可知,冲击频率为50~150 Hz时,轴扭复合冲击对钻头粘滑振动主振幅值的降低作用最明显。结合图9(a)的幅频局部放大图可知,冲击载荷减小,钻头粘滑振动的幅频表现为一阶扭转振动幅值降低。
5. 结 论
1)利用PDC钻头冲击破岩试验装置研究了固定转速条件下PDC钻头粘滑振动的发生过程及特征。采用20 r/min的驱动速度,逐渐提高送钻速度,PDC钻头发生粘滑振动并逐渐加剧,转速峰值和振动范围逐渐增大;粘滑振动对PDC钻头平均扭矩和平均钻速的影响不明显,对转速和扭矩瞬时波动具有明显的影响,其主振频率为扭转系统的一阶固有频率。
2)复合冲击载荷的加入改变了系统激励条件,引起转速在冲击载荷激励频率处出现明显的幅频响应特征,而系统扭转振动一阶固有频率处幅值有不同程度减小,且其幅值减小幅度越大,系统扭转振动缓解越明显。试验参数范围内,轴扭复合冲击载荷频率在50和100 Hz时对钻头粘滑的抑制作用最为明显,随着冲击频率继续提高,扭转振动主频对应幅值呈现增大趋势。
3)钻柱振动是深井超深井、超长水平段水平井目前面临的关键难点之一,PDC钻头破岩振动是引起底部钻具组合振动的主要激励源,建议开展PDC钻头实钻过程中的振动测量和分析,在大量实钻数据基础上建立适用于超深井的钻柱振动模型,以指导减振方法及高效减振工具的研究。
-
表 1 9种周期注水方案
Table 1 9 schemes of cyclic water injection
方案 高注入量注入
时间/d低注入量注入
时间/d停注时间/
d采油时间/
d1 20 20 10 40 2 30 30 10 60 3 40 40 10 80 4 20 20 20 40 5 30 30 20 60 6 40 40 20 80 7 20 20 30 40 8 30 30 30 60 9 40 40 30 80 -
[1] 孙龙德,邹才能,贾爱林,等. 中国致密油气发展特征与方向[J]. 石油勘探与开发,2019,46(6):1015–1026. SUN Longde, ZOU Caineng, JIA Ailin, et al. Development characteristics and orientation of tight oil and gas in China[J]. Petroleum Exploration and Development, 2019, 46(6): 1015–1026.
[2] 丁士东,赵向阳. 中国石化重点探区钻井完井技术新进展与发展建议[J]. 石油钻探技术,2020,48(4):11–20. doi: 10.11911/syztjs.2020069 DING Shidong, ZHAO Xiangyang. New progress and development suggestions for drilling and completion technologies in Sinopec key exploration areas[J]. Petroleum Drilling Techniques, 2020, 48(4): 11–20. doi: 10.11911/syztjs.2020069
[3] 张映红,路保平,陈作,等. 中国陆相致密油开采技术发展策略思考[J]. 石油钻探技术,2015,43(1):1–6. ZHANG Yinghong, LU Baoping, CHEN Zuo, et al. Technical strategy thinking for developing continental tight oil in China[J]. Petroleum Drilling Techniques, 2015, 43(1): 1–6.
[4] 李阳. 中国石化致密油藏开发面临的机遇与挑战[J]. 石油钻探技术,2015,43(5):1–6. LI Yang. Opportunities and challenges for Sinopec to develop tight oil reservoirs[J]. Petroleum Drilling Techniques, 2015, 43(5): 1–6.
[5] 李国欣,覃建华,鲜成钢,等. 致密砾岩油田高效开发理论认识、关键技术与实践:以准噶尔盆地玛湖油田为例[J]. 石油勘探与开发,2020,47(6):1185–1197. LI Guoxin, QIN Jianhua, XIAN Chenggang, et al. Theoretical understandings, key technologies and practices of tight conglomerate oilfield efficient development: a case study of the Mahu Oilfield, Junggar Basin, NW China[J]. Petroleum Exploration and Development, 2020, 47(6): 1185–1197.
[6] 李洪,李治平,王香增,等. 表面活性剂对低渗透油藏渗吸敏感因素的影响[J]. 石油钻探技术,2016,44(5):100–103. LI Hong, LI Zhiping, WANG Xiangzeng, et al. The effect of surfactants on imbibition-sensitive factors of low-permeability reservoirs[J]. Petroleum Drilling Techniques, 2016, 44(5): 100–103.
[7] 赵振峰,李楷,赵鹏云,等. 鄂尔多斯盆地页岩油体积压裂技术实践与发展建议[J]. 石油钻探技术,2021,49(4):85–91. doi: 10.11911/syztjs.2021075 ZHAO Zhenfeng, LI Kai, ZHAO Pengyun, et al. Practice and development suggestions for volumetric fracturing technology for shale oil in the Ordos Basin[J]. Petroleum Drilling Techniques, 2021, 49(4): 85–91. doi: 10.11911/syztjs.2021075
[8] 马剑,黄志龙,钟大康,等. 三塘湖盆地马朗凹陷二叠系条湖组凝灰岩致密储集层形成与分布[J]. 石油勘探与开发,2016,43(5):714–722. MA Jian, HUANG Zhilong, ZHONG Dakang, et al. Formation and distribution of tuffaceous tight reservoirs in the Permian Tiaohu Formation in the Malang Sag, Santanghu Basin, NW China[J]. Petroleum Exploration and Development, 2016, 43(5): 714–722.
[9] 王友净,宋新民,田昌炳,等. 动态裂缝是特低渗透油藏注水开发中出现的新的开发地质属性[J]. 石油勘探与开发,2015,42(2):222–228. doi: 10.11698/PED.2015.02.12 WANG Youjing, SONG Xinmin, TIAN Changbing, et al. Dynamic fractures are an emerging new development geological attribute in water-flooding development of ultra-low permeability reservoirs[J]. Petroleum Exploration and Development, 2015, 42(2): 222–228. doi: 10.11698/PED.2015.02.12
[10] 严谨,程时清,郑荣臣,等. 确定压裂裂缝部分闭合的现代产量递减分析方法及应用[J]. 石油钻采工艺,2018,40(6):787–793. YAN Jin, CHENG Shiqing, ZHENG Rongchen, et al. Development and application of the modern production decline analysis method in consideration of partial closure of hydraulic fracture[J]. Oil Drilling & Production Technology, 2018, 40(6): 787–793.
[11] FAN Tianyi, SONG Xinmin, WU Shuhong, et al. A mathematical model and numerical simulation of waterflood induced dynamic fractures of low permeability reservoirs[J]. Petroleum Exploration and Development, 2015, 42(4): 541–547. doi: 10.1016/S1876-3804(15)30047-1
[12] WANG Yang, CHENG Shiqing, ZHANG Kaidi, et al. Investigation on the transient pressure response of water injector coupling the dynamic flow behaviors in the wellbore, waterflood-induced fracture and reservoir: semi-analytical modeling and a field case[J]. International Journal of Heat and Mass Transfer, 2019, 130: 668–679. doi: 10.1016/j.ijheatmasstransfer.2018.09.083
[13] WANG Yang, CHENG Shiqing, ZHANG Kaidi, et al. Pressure-transient analysis of water injectors considering the multiple closures of waterflood-induced fractures in tight reservoirs: case studies in Changqing Oilfield, China[J]. Journal of Petroleum Science and Engineering, 2019, 172: 643–653. doi: 10.1016/j.petrol.2018.07.052
[14] 许锋,姚约东,吴承美,等. 温度对吉木萨尔致密油藏渗吸效率的影响研究[J]. 石油钻探技术,2020,48(5):100–104. doi: 10.11911/syztjs.2020114 XU Feng, YAO Yuedong, WU Chengmei, et al. Effect of temperature on the imbibition efficiency of the Jimusar tight oil reservoir[J]. Petroleum Drilling Techniques, 2020, 48(5): 100–104. doi: 10.11911/syztjs.2020114
[15] 何吉祥,徐有杰,高阳,等. 裂缝性致密油藏多级压裂水平井试井模型[J]. 断块油气田,2021,28(2):241–246. HE Jixiang, XU Youjie, GAO Yang, et al. Well test model of multi-stage fractured horizontal well in fractured tight reservoirs[J]. Fault-Block Oil & Gas Field, 2021, 28(2): 241–246.
[16] 赵思远,贾自力,吴长辉,等. 低渗透油藏注水诱发裂缝实验研究:以鄂尔多斯盆地吴起吴仓堡长9油藏为例[J]. 非常规油气,2021,8(3):73–79,89. ZHAO Siyuan, JIA Zili, WU Changhui, et al. Experimental study on waterflood induced fractures simulation in low permeability reservoir: a case study from Chang 9 reservoir in Wuqi Wucangpu, Ordos Basin[J]. Unconventional Oil & Gas, 2021, 8(3): 73–79,89.
[17] 孟勇,贾庆升,张潦源,等. 东营凹陷页岩油储层层间干扰及裂缝扩展规律研究[J]. 石油钻探技术,2021,49(4):130–138. doi: 10.11911/syztjs.2021094 MENG Yong, JIA Qingsheng, ZHANG Liaoyuan, et al. Research on interlayer interference and the fracture propagation law of shale oil reservoirs in the Dongying Sag[J]. Petroleum Drilling Techniques, 2021, 49(4): 130–138. doi: 10.11911/syztjs.2021094
[18] 吴忠宝,李莉,张家良,等. 低渗透油藏转变注水开发方式研究:以大港油田孔南GD6X1区块为例[J]. 油气地质与采收率,2020,27(5):105–111. WU Zhongbao, LI Li, ZHANG Jialiang, et al. Research on transformation of waterflooding development mode in low permeability oil reservoirs: taking GD6X1 Block of Kongnan in Dagang Oilfield as an example[J]. Petroleum Geology and Recovery Efficiency, 2020, 27(5): 105–111.
[19] 康毅力,田键,罗平亚,等. 致密油藏提高采收率技术瓶颈与发展策略[J]. 石油学报,2020,41(4):467–477. doi: 10.7623/syxb202004009 KANG Yili, TIAN Jian, LUO Pingya, et al. Technical bottlenecks and development strategies of enhancing recovery for tight oil reservoirs[J]. Acta Petrolei Sinica, 2020, 41(4): 467–477. doi: 10.7623/syxb202004009
[20] 王增林,鲁明晶,张潦源,等. 东营凹陷陆相页岩油强化缝网改造生产制度优化研究[J]. 石油钻探技术,2021,49(4):71–77. doi: 10.11911/syztjs.2021074 WANG Zenglin, LU Mingjing, ZHANG Liaoyuan, et al. Production system optimization for enhanced fracture network stimulation in continental shale oil reservoirs in the Dongying Sag[J]. Petroleum Drilling Techniques, 2021, 49(4): 71–77. doi: 10.11911/syztjs.2021074
-
期刊类型引用(1)
1. 郭海伟. 基于增量学习的储气库井底压力快速计算方法. 断块油气田. 2025(02): 292-299 . 百度学术
其他类型引用(0)