缝洞型油藏溶洞内水驱流态研究

缝洞型碳酸盐岩油藏主要由细小的裂缝与溶洞连接而成，缝洞作为主要储集体，裂缝为流动通道的强非均质油藏[1-3]。裂缝低孔高渗的特性和溶洞高孔高渗的特性符合单孔介质的特征，使缝洞型油藏不同于常规的双孔介质油藏，开采机理与普通双孔介质油藏有所不同[4-6]。邢翠巧等[7]针对缝洞型碳酸盐岩油藏经过构造及岩溶作用之后的情况，建立了溶洞及裂缝发育的双孔模型，研究了缝洞型碳酸盐岩油藏实际储层特征和流体的流动机理。KANG等[8]和WU等[9]应用连续介质渗流理论建立了缝洞型油藏多相流问题的数值模拟方法。郑小敏等[10]利用全直径缝洞型岩芯模型，通过物理实验模拟水驱油实验研究了缝洞型碳酸盐岩油藏的渗流特点。康志宏[11]针对塔河油田4区奥陶系缝洞型碳酸盐岩油藏开展渗流机理模拟试验研究，得到碳酸盐岩油藏中裂缝(或大的溶洞)是主要渗流通道的结论。

缝洞型油藏由于储集体类型多样，微孔隙、微裂缝、裂缝和中型溶洞以及大型溶洞、大型裂缝的分布，在尺度和空间结构上都具有很大的差异，流体流动类型复杂，渗流、裂缝流和洞穴流都有存在[12-15]。赵欢等[16]利用油藏数值模拟方法，研究了不同聚合物参数对“二三结合”开发模式原射孔层和补孔层的压力、含水饱和度等渗流特征的影响。刘金玉[17]利用COMSOL Multiphysics多物理场软件对单缝概念模型进行了水驱油数值模拟，分析了油水流态、波及面积、水驱前缘及剩余油分布的基本特征。目前对缝洞型碳酸盐岩油藏水驱的研究多是从渗流角度出发，但是由于缝洞型油藏内部结构复杂，裂缝和溶洞的存在情况无法准确掌握，特别是针对一些大裂缝和大溶洞，内部属于未填充状态，多孔介质没有在内部分布，所以内部流体流动不属于渗流范畴，它的流动规律符合Navier-Stokes方程，因此需采用流体力学对其进行研究。

本文从流体力学角度出发，针对未填充缝洞型碳酸盐岩油藏进行了油水两相流动模拟，研究单裂缝-溶洞中流体流动时油藏内原油黏度、注水速度对油水界面、波及面积的影响。

1 数学模型

1.1 连续方程

连续性方程是流体运动学的基本方程，是质量守恒原理的流体力学表达式。

连续方程微分形式见式(1)。

 

(1)

式中：ρ表示流体密度，kg/s；t表示时间，s；u表示流速，m/s。

模型设定的油水两相流动为稳态流动，即：

 

(2)

油水两相都是不可压缩流体，即：

=0 。

(3)

将以上条件代入连续方程，得到油水两相流的连续方程表达式：

 

(4)

式(4)即用来描述缝洞型介质两相流动机理的连续方程，可以通过连续方程得到流体的流速，代入纳维-斯托克斯方程。

1.2 纳维-斯托克斯方程

纳维-斯托克斯方程简称N-S方程，用来描述黏性不可压缩流体动量守恒的运动方程[18]。

根据开关状态的不同,三相功率变换器共有8个输出电压矢量,V0(000),V1(100),V2(110),V3(010),V4(011),V5(001),V6(101),V7(111),如图2所示,三相功率变换器器输出的电压矢量可以表示为

REN Wenbo, CHEN Xiaofan. Study on asymmetric instability water filling in fractured-hole carbonate reservoirs [J]. Science and Technology, 2013,13 (27): 8120-8125.

 

(5)

式中：表示单位体积上的惯性力；表示单位体积上的质量力；表示单位体积上的黏滞力。

由于油水两项在大溶洞和大裂缝中的运动属于自由流动，遵循N-S定律，同时为了结合VOF方法追踪溶洞中油水界面问题，所以采用张量形式的不可压缩黏性N-S方程，用于研究未填充溶洞中的油水两相流动问题。

1.3 VOF方法及原理

VOF方法基本原理是通过对网格单元中流体之间的体积比函数F来确定自有界面，追踪界面的推进情况。由Hirt和Nichols首次提出后经NOH等改进，在保证质量守恒的前提下，PUCKETT等采用方向分裂算法对输送方程进行优化，利用流体体积函数方程张量表示方式如下：

 

(6)

由于油的密度小于水的密度，两相之间存在明显界面，所以采用VOF模型对油水界面进行跟踪。根据单位网格内油水体积比F得到油水界面形状和位置。当F=0时，说明单元内全部是油相；当F=1时，说明单元内全部是水相，已经将油驱出；当0<F<1时，说明单元内存在油水两相，有油水界面存在。

党内政治生活，是党内各级组织和党员按党章及党内规章制度和优良传统习惯要求，处理政治关系、形成良性政治生态的主要平台。党的十八大以来，习近平总书记站在推进党的建设新的伟大工程战略高度，从严肃党内政治生活的重要作用、重点内容、主要途径、基本方式、衡量标准、责任主体等方面进行论述，形成了关于严肃党内政治生活的重要思想。研究探析习近平严肃党内政治生活思想，对于在新时代中国特色社会主义建设征程中推进全面从严治党向纵深发展具有重要的理论意义和实践价值。

根据流体体积函数方程张量表示方式，通过流体流动时间、速度以及流动方向上的位置可以计算出单位网格内油水体积比，利用瞬态计算方法，得到不同时刻的油水界面分布情况。

LI Aifen, ZHANG dong, GAO Chenghai.The law of water injection and oil extraction in closed volume cells[J]. Petroleum Geology and Recovery, 2012,19 (3): 94-97.

2 数值计算

2.1 初始条件

单裂缝是缝洞型油藏的基本流动单元，对流体在油藏中的流动起决定性作用。假设油和水都是不可压缩流体，密度为常数。由于溶洞三度空间长度相近，同时为了减小计算量，将溶洞简化为正方形二维模型，出入口位于左右两侧中心，左端为出口，右端为入口。入口定义为速度入口，出口定义为压力出口。

2.2 模型参数及网格划分

参数确定：直径大于100 mm属于巨洞；直径在10～100 mm属于大洞；直径在5～10 mm之间属于中洞；直径在2～5 mm之间属于小洞。由于溶洞空隙三度空间相近，且都大于2×103 μm，所以选用三度尺寸均为50 mm的溶洞，缝宽为0.6 mm。实际缝洞型碳酸盐岩油藏地层原油黏度一般为2～30 mPa·s，所以选用原油黏度为5，15和25 mPa·s。原油密度选取860 kg/m3，水的密度为996 kg/m3。

溶洞模型为边长50 mm的正方形二维模型，出入口缝宽为0.6 mm，利用gambit对模型进行六面体网格划分，见图1。根据雷诺数计算，采用Fluent中的层流模型，且为了方便追踪油水界面，采用VOF方法进行瞬态计算，收敛精度为0.000 01。

  

图1 网格划分图Fig.1 Mesh map

首先，计算入口速度为0.08 m/s、黏度为5 mPa·s的溶洞内油藏的流态，得到水驱油水分界面分布，波及面积及剩余油分布，依次改变黏度为15和25 mPa·s，得到不同黏度下的油水界面发展情况、波及面积及剩余油分布情况。改变原油黏度，得到相同入口速度下不同黏度的油水分布情况，分析原油黏度和流态之间关系；再改变入口速度，依次设定入口水速为0.08，0.16，0.24 m/s，得到原油黏度为15 mPa·s时不同入口速度下的油水分布情况，分析入口水速和流态之间的关系。

针对目前国内信息化建设普遍存在的问题，我院进行了大量的调研和学习，基于全院上下对现代医院信息化建设的深度认同，进行了科学详细的方案规划和明确的目标设定。

3 计算结果和分析

3.1 一定入口速度下黏度变化对波及面积的影响

通过图2可以看到在入口速度为0.08 m/s时完整的水驱过程波及面积变化，黏度对单缝洞型油藏中水驱波及面积的影响很大：在黏度为5 mPa·s时，水相趋向中心线，向洞口流动现象明显，油藏内形成漩涡，波及面积最小，水相占据比例为27%；当黏度增大到15 mPa·s时，开始向活塞驱替转变，没有形成漩涡，波及面积增大，水相占据比例为52%；当黏度增大到25 mPa·s时，呈现水滴状，波及面积达到最大，水相占据比例为61%。

  

图2 不同黏度下波及面积图Fig.2 Affected area under different viscosity

出现这种现象的原因：在水驱入口速度为0.08 m/s、原油黏度从5 mPa·s增大到25 mPa·s时，流动性变差，对注入水的扩散阻力增大，相间磨损增加，水相推进速度变小，油水界面到达出口的时间变长，扩散作用变大；同时，由于原油黏度增大，油水界面向前推进的速度减慢，水相沿中心线纵向扩散趋势逐渐增强。因此，类射流现象随着原油黏度增加逐渐转换成活塞驱替，漩涡逐渐消失，波及面积逐渐增大。

3.2 一定入口速度下黏度变化对油水界面的影响

为研究水驱油过程中油水界面的发展形态，从上述3种黏度下分别取6幅不同时刻的流态图，见图3。

  

图3 不同黏度下的油水界面图Fig.3 Oil-water interface under different viscosity

根据《饮料商务》报道，中国线上酒类交易很集中，其中领头的天猫和京东占据了70%以上，这两家的总额加起来基本等于美国线上交易总额的3倍；同时，中国线上酒类交易的总额是排在其后的法国（或英国）的3倍。在这其中，成熟的专业葡萄酒零售商起到了极大的促进作用。IWSR的调查显示，大多数国家的线上酒类市场购买频率较低，只有中国和英国有超过50%的消费者至少每月购买一次或更多。

3.3 一定黏度下入口速度变化对波及面积的影响

由图4可知，入口速度对波及面积影响显著，油水界面到达右侧缝洞出口时，波及面积随着入口速度增加而减小。通过图4中完整的水驱过程波及面积图可知，入口速度对单缝洞型油藏中水驱波及面积的影响很大：在速度为0.08 m/s时，波及面积呈现水滴状，油藏内无漩涡产生，波及面积最大，水相占据比例为59%；当入口速度为0.16 m/s时，水相流动有向中心线靠近的趋势，漩涡开始形成，波及面积变小，水相占据比例为43%；当入口速度为0.24 m/s时，射流现象明显，产生漩涡，波及面积最小，油水界面推进速度明显下降。

  

图4 不同速度下的波及面积图Fig.4 Affected area under different velocity

[5] 修乃岭，熊伟，高树生，等．缝洞型碳酸盐岩油藏流动机理初探[J].钻采工艺,2008,31(1):63-65.

大学生大部分远离家乡在另外一个城市求学，因此，父母对大学生的学习生活的关注不及中小学时。虽然大学生远离家庭，但家庭作为其出生和成长的重要环境，其在外的言行仍与家庭有着千丝万缕的联系，家庭经济困难、父母感情破裂等家庭因素也是大学生恋爱危机的重要诱因。因此大学生的恋爱危机干预模式应该从学校单一的危机干预模式转为学校、家庭合作模式。当大学生发生恋爱危机时，学校应立即启动危机处理预案，相关人员要根据危机情况尽快与学生家长取得联系，做好沟通工作。对家长而言，当孩子处于恋爱危机时，最好能陪伴孩子，给予情感支持与精神抚慰，同时，家长应与学校保持沟通，反馈相关信息，与学校共谋对策。

3.4 一定黏度下入口速度变化对油水界面的影响

为了研究水驱油过程中油水界面的发展形态，从上述3种速度下分别取6幅不同时刻的流态图，见图5。

当线路短路故障引起断路器跳闸后，负荷中无法自启动的设备应自动跳闸。而负荷中允许重新带电的设备，则可以在断路器重合成功后恢复供电。重合闸时间的整定主要取决于短路点电弧的去游离时间。用户负荷中不能参加自启动的设备须自动跳闸。允许重新带电的负荷，在断路器重合后恢复供电，重合闸时间主要决定于故障点去游离时间。

  

图5 不同速度下的油水界面图Fig.5 Oil-water interface under different velocity

[8] KANG Z J，WU Y S．Modeling multiphase flow in naturally fractured vuggy petroleum reservoirs[C]/ /The 2006 SPE Annual Technical Conference and Exhibition. Texas: SPE，2006:1-10．

3.5 准确性论证

模型采用两套网格对比，对网格无关性进行分析，第1套网格数量为16 283、节点数量为16 024，采用四面体结构网格。第2套网格数量为52 717、节点数量为54 238，采用三角形网格。将两套网格条件下模型的垂直中心轴线的铅垂线上的油相体积分布进行对比和模拟验证，见图6。由图6可知，平均相符程度为99.33%，可以认为两套网格的模拟结果相同，能做到网格无关性。基于运算量考虑，第1套网格更小，运算更快，因此后续使用第1套结构网格进行数值模拟分析。

  

图6 含油率分布曲线图Fig.6 Oil content distribution curve

4 结 论

在单缝洞型油藏水驱开采过程中，原油黏度和注水速度的改变会在很大程度上影响波及面积和油水界面，从而影响水驱效率。

1)相同注水速度下，黏度越大，油水界面推进速率越低，原油黏度较低时会在出口产生回流，浪费能量达不到水驱预期效果，波及面积随着黏度增大逐渐增大，推进速度降低。

2)原油黏度一定时，波及面积随着入口速度的增加而减小，油水界面呈对称趋势均匀推进，水速较大时会发生类射流现象并有漩涡产生，随着入口水速增大，类射流现象和漩涡现象逐渐明显。

3)由于水驱速度和原油黏度的共同影响，产生漩涡和类射流现象，在水驱过程中造成能量浪费，对于不同黏度原油需要选择合理的注入水速度，达到预期水驱效果。

由于缝洞型油藏储层地质特征复杂，溶洞内流体流动形态受到多种因素影响。本文从流体力学角度出发研究了原油黏度和注水速度对波及面积和油水界面的影响，对内部流场压力和速度等具体分布情况缺少研究，缝洞比、储层温度、溶洞所处地层深度等因素对溶洞内流体流态分布都有待进一步研究。

参考文献/References：

[1] 任文博，陈小凡．缝洞型碳酸盐岩油藏非对称不稳定注水研究[J].科学技术与工程,2013,13(27):8120-8125.

溶洞型油藏中油水两相流动属于黏性流体的不可压缩流动，由于需要研究油水界面在不同黏度原油和不同入口速度下的推进情况，所以采用张量形式的不可压缩黏性N-S方程，即：

[2] 胡蓉蓉，姚军，孙致学，等．塔河油田缝洞型碳酸盐岩油藏注气驱油提高采收率机理研究[J].西安石油大学学报(自然科学版),2015,30(2):49-53.

HU Rongrong, YAO Jun, SUN Zhixue, et al. Research on mechanism of enhanced oil recovery by gas-flooding in fractured-cavity carbonate reservoirs in Tahe oilfield [J]. Journal of Xi’an Shiyou University (Natural Science Edition), 2015, 30 (2): 49-53.

[3] 李爱芬，张东，高成海．封闭定容型缝洞单元注水替油开采规律[J].油气地质与采收率,2012,19(3):94-97.

控制了技术转移与扩散对产业升级的效果后，两组技术进步与金融结构的交叉系数方向一致，其中一组结果在5%置信程度上显著为正，表明适宜的金融结构(特别是金融市场的占比)总是能通过技术创新对产业结构升级产生正向的促进作用。这种促进一方面通过缓解企业的融资约束，促进实体企业的技术创新活动和研发投入，提高技术创新的产出；另一方面则通过有效的金融资源配置，激励“企业家精神”创建更高效率的新兴企业，从而促进技术创新加快进入生产应用，推动产业结构的高度化。

[4] 修乃岭. 缝洞型碳酸盐岩油藏流动机理研究[D].廊坊：中国科学院研究生院(渗流流体力学研究所),2008.

我挣脱开，骂他。没见过这么不要脸的人，都快老死了，还说要跟我结婚。我顺手用楼道里的煤球砸他，弄得他一身都是黑煤印，灰溜溜地走了。我不怕找不到他，这个畜生肯定还会再来的。

XIU Nailing. Study on Flow Mechanism of Fractured-hole Carbonate Reservoir [D]. Langfang：Graduate School of the Chinese Academy of Sciences (Institute of Percolation and Hydrodynamics), 2008.

出现这种现象的原因：在原油黏度为15 mPa·s时，水驱入口速度从0.08 m/s增加到0.24 m/s，由于注入速度与压力梯度成正比，随着入口速度增加，压力梯度增大，油水界面向前推进的速度随之增加，水相推进速度变大，油水界面到达出口时间变短，从18 s减少到8 s再减少到4 s，相间扩散作用变小；同时，由于入口速度增大，水相更容易沿孔洞中间区域到达右侧出口端。因此，波及面积随着注入速度增加活塞驱替向射流发展，漩涡逐渐形成，波及面积逐渐下降。

XIU Nailing, XIONG Wei, GAO Shusheng, et al.Study on the flow mechanism of fractured-cavity carbonatite reservoir [J]. Drilling & Production Technology, 2008,31(1): 63-65.

[6] 孙玉平，修乃岭，熊伟，等．缝洞型碳酸盐岩油藏流动数学模型初探[J].石油钻探技术,2008,36(1):65-68.

SUN Yuping, XIU Nailing, XIONG Wei, et al.Study on the flow mathematical model of fractured-cavity carbonate reservoir [J]. Petroleum Drilling Technology, 2008,36 (1): 65-68.

[7] 邢翠巧,尹洪军,李兴科,等.缝洞型碳酸盐岩油藏试井分析模型研究[J].河北工业科技,2018,35(1):31-36.

XING Cuiqiao,YIN Hongjun,LI Xingke,et al.Study on well test analysis model of fracture-cavity carbonate reservoir[J].Hebei Journal of Industrial Science and Technology,2018,35(1):31-36.

由图5可知，未到达出口前，3种速度下的油水界面都是按照均匀对称的方式向前推进，到达出口附近时，由于压力变化，油水界面会产生“尖嘴状”。当速度为0.08 m/s时，水相优先于中线位置通过现象并不明显，油水界面到达出口时间为17 s；当速度为0.16 m/s时，水相沿中线向右侧出口推进的趋势增强，漩涡开始产生，油水界面均匀向前推进，到达出口时间为8.6 s；当速度为0.24 m/s时，油水界面以圆形均匀前进，到达出口时间为4 s，明显产生漩涡，部分原油被卷入漩涡中。

[9] WU Y S，DI Y，KANG Z J，et al．A multiple-continuum model for simu-lating single-phase and multiphase flow in naturally fractured vuggy reservoirs[J]．Journal of Petroleum Science and Engineering，2011，78(1):13-22．

[10] 郑小敏,孙雷,王雷,等.缝洞型碳酸盐岩油藏水驱油机理物理模拟研究[J].西南石油大学学报(自然科学版),2010,32(2):89-92.

ZHENG Xiaomin, SUN Lei, WANG Lei, et al.Physical simulation of water displacing oil mechanism for vuggy fractured carbonate rock reservoir[J]. Journal of Southwest Petroleum University(Science & Technology Edition), 2010,32 (2): 89-92.

[11] 康志宏．缝洞型碳酸盐岩油藏水驱油机理模拟试验研究[J].中国西部油气地质,2006,2(1):87-90.

在医院。我醒来时，儿子说，我已睡了很久。他说他一直在旁边看着输液管里的液体，怕它停止嘀嗒。我的眼泪一下子流了出来。

KANG Zhihong.Simulation test study on mechanism of water flooding in fractured carbonate reservoir [J]. Petroleum Geology, West China, 2006,2(1): 87-90.

[12] 杨坚，程倩，李江龙．塔里木盆地塔河4区缝洞型油藏井间连通程度[J]．石油与天然气地质，2012，33(3):484-489．

YANG Jian，CHENG Qian，LI Jianglong．Interwell-connectivity analysis for fracture-vug reservoirs in the block-4 of Tahe oilfield，Tarim Basin[J]．Oil & Gas Geology，2012，33(3):484-489．

[13] 康志江，李江龙，张冬丽．塔河缝洞型碳酸盐岩油藏渗流特征[J]．石油与天然气地质，2005，26(5):634-640．

KANG Zhijiang，LI Jianglong，ZHANG Dongli．Percolation characteristics of fractured-vuggy carbonate reservoir in Tahe oilfield[J].Oil & Gas Geology，2005，26(5):634-640．

由图3可知，未到达出口前，3种黏度下的油水界面都是按照均匀对称的方式向前推进，到达出口后，由于压力变化，油水界面会产生“尖嘴状”。当黏度为5 mPa·s时，水流受原油阻力较小，水相优先于中线位置通过，由于水相扩散作用产生漩涡，部分原油受到漩涡波及，在漩涡作用下向右移动，水相到达出口时间最短，耗时14 s；黏度增大到15 mPa·s时，水相受到的阻力增大，水相沿垂直中线方向扩散趋势增强，产生漩涡趋势减弱，油水界面均匀向前推进，水相到达出口的时间变长，耗时17 s；黏度增大到25 mPa·s时，在到达出口前，油水界面前缘以圆形均匀前进，靠近出口时，前缘呈现“尖嘴状”，水相到达出口的时间最长，耗时19 s。

[14] 康志江，赵艳艳，张冬丽．缝洞型碳酸盐岩油藏数值模拟基础理论与应用[J]．大庆石油地质与开发，2011，30 (5):68-71．

我见过一位父亲，一心要把孩子培养成牛人，于是早早给他买《曾国藩》《胡雪岩》；另一位母亲，则忧心忡忡，认定《傲慢与偏见》《红楼梦》都只会把她的花季女儿教坏，索性根本不让女儿看“闲书”：“有那时间背背单词不好吗？”

KANG Zhijiang，ZHAO Yanyan，ZHANG Dongli．Basic theory and application of numerical simulation in fractured-vuggy carbonate reservoir[J]．Petroleum Geology and Exploitation in Daqing，2011，30(5):68-71．

[15] 李宗宇．塔河奥陶系缝洞型碳酸盐岩油藏开发对策探讨[J]．石油与天然气地质，2007，28(6):856-862．

1.1 资料来源 选取2015年5月-2017年4月于绍兴市妇幼保健院定期行孕期检查的肥胖孕妇160例。纳入标准：①单纯肥胖孕妇；②孕前体质量指数BMI≥28 kg/m2；③单胎妊娠孕妇；④骨盆外测量值正常。排出标准：①患有高血压、糖尿病、心脏病及免疫性疾病者；②习惯性流产；③年龄<20岁者。按照随机数字表法，将160例肥胖孕妇随机分为对照组和观察组，每组80例。其中，观察组年龄21～36岁，对照组年龄22～34岁。两组肥胖孕妇一般资料差异无统计学意义，具有可比性。见表1。本研究经孕妇及家属知情同意，且经过医院伦理委员会同意批准。

LI Zongyu．A discussion on development strategies for the Ordovician fractured-vuggy carbonate rock reservoirs in Tahe oilfield[J]．Oil & Gas Geology，2007，28(6):856-862．

[16] 赵欢,尹洪军,刘红,等.注聚参数对L油田B区块“二三结合”渗流特征的影响[J].河北科技大学学报,2015,36(3):324-329.

课后训练需要将作业、练习与生活实际相结合，教师在布置作业的时候尽量选取与生活相联系的题目，以此提高学生理论联系生活的能力。还可以开展有关物理的生活实践活动，在探索大自然的过程中加以学习，增强学生对物理应用价值的认识[2]。

ZHAO Huan, YIN Hongjun, LIU Hong, et al. Effects of poly condensation parameters on the seepage characteristics of “two-three combination” in block B of L oilfield[J].Journal of Hebei University of Science and Technology,2015,36(3):324-329.

[17] 刘金玉. 缝洞型介质两相流动机理的数值模拟研究[D].北京：中国石油大学,2009.

免费医学生也称农村订单式定向医学生，根据国家发展和改革委员会、原卫生部、教育部等部门在2010 年印发《关于印发开展农村订单定向医学生免费培养工作实施意见的通知》，免费医学生在入学前必须与培养学校和当地县级卫生行政部门签署定向就业协议，毕业后需要到相应的基层医疗机构服务满6年，入学后免收学费、住宿费并有一定的生活补助。而普通医学生毕业后需要进行双向选择就业，在校期间也没有相应的助学金政策[3]。

LIU Jinyu.Numerical Simulation of Two-Phase Flow Mechanism in Fractured-Cavity Media [D]. Beijing：China University of Petroleum, 2009.

[18] 张少峰，曹会敏，刘燕，等.弯管中液固两相流及壁面碰撞磨损的数值模拟[J].河北工业大学学报, 2008，37(3):48-54.

ZHANG Shaofeng, CAO Huimin, LIU Yan, et al.Numerical simulation of collision and wear between liquid-solid two-phase flow and wall in bend [J].Journal of Hebei University of Technology, 2008，37(3): 48-54.