2016, Vol. 31
油藏中聚集的原油，在经过不同的开采方式和开发阶段后，仍然滞留在油藏不同地质环境中的那一部分称作剩余油.按照剩余油“体积规模”的不同，可将剩余油的研究分为微规模、小规模、大规模和宏规模，不同规模的剩余油研究互相联系(赵军龙等, 2013).大尺度规模剩余油各种特性是微观尺度剩余油的整体表现，研究微观剩余油形成机理、赋存量以及空间分布对驱替剂驱油机理、两相及多相驱替的渗流规律等机理研究有重要的指导作用，在工程应用中可以改善开采工艺，增加可采储量，提高采收率.从储层微观孔隙结构角度研究剩余油的形成过程、数量、微观分布及其重新运移富集方式，对于正确揭示剩余油微观特性和科学开采剩余油具有重要意义.
本文将常用的微观剩余油研究分为物理实验技术和数值模拟方法.物理实验技术讨论了传统光学技术、扫描电镜技术、激光共聚焦技术、CT技术、核磁共振成像技术；数值模拟方法分析了基于网格化方法、无网格化方法、孔隙网络方法等.对各方法的优缺点加以评述，以期在勘探和开发研究微观剩余油时采用技术方法的选择作基础性的工作.
1 物理实验技术研究现状随着新的理论发现和实验仪器的发展，从最早的二维空间传统光学技术到如今的三维空间无损检测技术，已经形成了多种微观剩余油分布的物理实验技术.根据实验处理过程的不同，这些物理实验技术主要包括传统光学技术、扫描电镜技术、激光共聚焦技术、CT技术、核磁共振成像技术.
受实验条件的约束，国内物理实验技术研究微观剩余油分布长期以传统光学技术为主，由于其价格相对较低、技术手段成熟，目前已经成为各石油高校和油田研究微观剩余油的基本实验手段.国内不同油田的地质条件差异较大，储层岩性复杂，因此在薄片制备方面相对国外研究有着更加成熟的经验.尽管CT技术是近二十年才在国内应用发展，但随着高精度岩心CT成像设备的大量引进，利用CT技术研究微观剩余油逐渐成为热点.目前，国内微观剩余油研究主要依托中石化胜利油田地质科学研究院、中国石油勘探开发研究院的CT设备与配套驱替设备，两所研究院团队及中国石油大学(华东)和中国石油大学(北京)两所高校的研究学者(侯健团队等)多年来从不同岩石孔隙结构、不同驱替方式等多个角度研究微观剩余油空间分布，揭示多相流体在孔隙中流动机理，使我国在微观剩余油研究领域具有先进水平.国内利用核磁共振技术研究微观剩余油的团队较多，二十世纪九十年代以北京石油勘探开发研究院万庄分院核磁设备(BMT 47/40)为主，利用该设备对不同岩心进行驱替实验的可视化研究，得到许多提高采收率机理和效果评价的成果.随着核磁共振仪器硬件水平的提高及国内核磁共振成像仪器生产公司的增多(如苏州纽迈分析仪器股份有限公司、北京斯派克科技发展公司等)，利用核磁共振成像研究岩心中微观剩余油的团队呈现快速增长的趋势，包括中石化石油勘探开发研究院、中原油田分公司、中海油天津分公司、大连理工大学(宋永臣团队)、上海大学(狄勤丰团队)、辽宁工程技术大学(潘一山团队)等.由于各团队研究的主要方向有差异，丰富了微观剩余油在岩心尺度上形成机理、空间分布等方面的研究成果.扫描电镜技术、激光共聚焦技术存在研究尺度相对较小等缺点，国内利用相关技术研究微观剩余油的团队相对较少.
1.1 传统光学技术传统光学技术是最早研究微观剩余油形成机理、赋存量以及空间分布的实验技术，主要包括含油薄片技术、真实岩心模型和微观仿真模型，这些技术不需要高精度的仪器设备，是微观剩余油物理实验研究的常用手段.
含油薄片技术是将不受污染的含油岩心(油基泥浆、密闭、大直径取心)进行切割，利用低黏度、高强度、折光率接近玻璃的黏合剂在常温下自然渗胶、固结制作而成，它基本保持了剩余油的原始分布状况.孙国忠(2002)通过对不同类型的饱含油岩心样品薄片磨制，总结对于不同的岩心样品的磨制过程，认为岩样固结技术是制作含油薄片的关键.含油薄片技术是各个油田认识储集层微观剩余油的主要手段，不同储层、不同含水阶段微观剩余油特性既有相似性，也有鲜明的独特性(李健和李红南, 2003, 孙盈盈等，2014).王朴等(2002)用含油薄片研究胜坨油田二区沙二段砂岩高含水油层微观剩余油分布特征，根据剩余油形成的机理和在不同孔隙尺寸中分布形式，将微观剩余油划分为珠状、凝块状、环状、连片状、条带状六种形态，并讨论了剩余油分布的影响因素，主要有岩石的层理、孔隙结构、岩石的结构成分及其润湿性.含油薄片技术的最大优势在于直观真实反映剩余油在储集层孔隙中原始微观分布特点，可以评价取心井所在区域的含油性和水淹程度，也是油田研究不同储层微观剩余油空间分布最常用的技术.但它仅仅反映二维空间局部剩余油空间分布特征，只能近似代替实际赋存量，而且属于静态观测，难以进行微观剩余油形成机理研究.
真实岩心模型的常规方法是以油藏具有代表性的岩心为实验研究对象进行驱替实验，应用薄片技术分析，直接描述剩余油的形成机理、赋存量以及空间分布.宋考平(2005)对储层真实岩心进行水驱油实验，再将驱替后岩心截取三个横断面制成薄片，利用磨片荧光分析技术进行分析，研究聚合物驱的微观剩余油形成机理，统计不同类型微观剩余油的比例，给出高、中、低不同强度水淹部位的微观剩余油赋存量和空间分布规律.向丹和黄大志(2005)认为受模型材料限制，目前驱替实验不能模拟地层环境条件，通过反复实验优选配制了地层条件下模拟油，在实验室模拟水驱过程，利用驱替后制作的岩心薄片研究微观剩余油特性.相对于含油薄片技术，该技术可以进行动态驱替实验，研究微观剩余油形成机理、聚合物驱油效率，缺点是不能反映流体渗流的动态特征，而只能观察结果.另一种方法是直接对岩心磨制薄片，再进行驱替，这也被称作真实砂岩微观模型(黄书先和张超谟, 2004).它用全直径岩心磨成很薄的薄片，再用两块玻璃夹住，并用橡胶把模型薄片周围黏结起来，进行水驱油模拟实验.公言杰等(2015)以四川盆地致密油储层岩心样品为例，结合真实砂岩微观模型与核磁共振技术，研究了致密油充注过程油水流体分布及润湿性的变化模式. 赵阳(2002)对裂缝性砂岩进行水驱油的单模型和组合模型模拟实验，研究水驱油过程中注入水在孔道和裂缝中的微观水驱油机理、剩余油形成机理和裂缝对驱油效率的影响.由于采用实际岩心薄片作为水驱油物理模型，因此该方法能较为准确的模拟不同含水阶段的油水分布特征.但由于是岩心驱替后再进行薄片磨制，对于同一处薄片该过程不可重复.
微观仿真模型的主要研究对象是孔隙网络，根据制作工艺不同，主要可以分为指定孔隙网络模型和岩石刻蚀模型.这种模型并不是岩心真实的孔隙，而是根据研究需要制作指定孔隙网络模型，或者根据显微镜下储层孔隙结构特征(孔隙及喉道大小、形态、连通情况等)，在玻璃上刻画出相同孔隙特征的模型.Chatenever和Calhoun (1952)利用玻璃球进行多孔介质驱替过程的可视化研究，为动态观测孔隙中流体渗流的特征实验提供了思路.指定孔隙网络模型并不具有真实岩石的孔隙网络特征，但为微观剩余油形成机理研究提供指导. Jamaloei等(2010)在特殊规则孔隙网络模型中，比较高、低界面张力作用下微观剩余油形成机理与空间分布的不同.Al-Shalabi和Ghosh (2016)制作了低渗透率到高渗透率的六块不同微观模型，通过驱替实验确定油水分布状态，并发现低渗区域具有更高的最终采收率，认为微观孔隙不均匀和指进是这种现象的主要因素.岩石刻蚀模型制作相对复杂，但具有真实岩心孔隙网络特征，在此网络中进行水驱油实验可以近似代表岩石孔隙结构对微观剩余油的影响.1961年Mattax首次制作出刻蚀玻璃网络，随后Davis和Jones (1968)用光刻技术促进了刻蚀模型在研究微观剩余油的应用.McKellar和Wardlaw (1982)将岩石孔隙结构照片放在涂有感光材料的玻璃上，经过曝光、显影将孔隙结构图形复制在玻璃上，用氢氟酸腐蚀，使其显现孔隙结构印模，再加上盖板制成二维微观孔隙网络模型.这种技术被称作化学刻蚀法，该技术成熟，孔隙大小可控，模型透明、可视性好，比树脂浇注技术更容易控制孔隙的粗糙度和润湿性.国内广泛应用岩石刻蚀模型水驱油实验研究微观剩余油形成机理与分布特征，考虑了不同孔喉结构、原油黏度、驱替压差、驱替剂、润湿性等因素对微观剩余油的影响.王芳芳等(2013)根据微观模型水驱油实验的观察和总结，从剩余油的形成机理和赋存状态两个方面进行了研究综述，并对剩余油的赋存状态重新进行划分.目前针对微观剩余油研究成果多以直接观察和定性描述为主，定量研究较少.侯健等(2006)利用图形学概念中的形状指数定量描述微观剩余油几何形状特征，并研究了聚合物驱不同阶段各类剩余油的分布比例.马炳杰(2016)应用刻蚀模型的微观可视化物理模拟技术及图像处理技术对剩余油变化定量表征，研究注水停止前后剩余油分布的变化规律、影响因素、运移方式及受力情况，发现注水停止前后微观剩余油分布不同，而且可以分为三种运移模式.尽管微观仿真模型并不能全面反映真实岩心特征，但由于其良好的可视性及驱替实验的可重复性，该技术已经成为目前实验室动态研究微观剩余油的主要手段.
1.2 扫描电镜技术扫描电镜具有分辨率高、放大倍数大、景深大、立体感强等优点(单秀琴和严启团，1999)，广泛应用于岩石微观结构、矿物成分、孔隙类型、胶结程度等方面的研究.它基本方法是电子枪发射的电子束加速后在样品表面聚焦成极细的电子束，在双偏转线圈作用下在样品表面进行扫描，激发的各种信号被处理成像，从而在荧光屏上观察到样品表面的各种特征.但用扫描电镜进行测试时需要将样品进行干燥处理，限制了对岩心的含油性分析.在扫描电镜基础上发展起来的冷冻扫描电镜技术(Cryo-Scanning Electron Microscopy)和环境扫描电镜技术(Environmental Scanning Electron Microscopy)可以对微观剩余油分布进行研究，同时保证了扫描电镜具备的优势.
冷冻扫描电镜技术需要将样品快速冷冻，以保持样品内部的流体分布，减少冰晶的产生，再进行油水分布的扫描观测.Sutanto等(1990)研究了强水湿砂岩和混合水湿砂岩在驱替和渗吸的不同阶段微观剩余油分布的特点，利用冷冻扫描电镜技术清晰的展示了粘土分布的不同对油水在孔隙中分布的影响.Vizika等(1998)利用该技术研究了不同润湿性条件下，油气水三相流存在时微观剩余油分布.由于该技术独特的优势，可以直观显示颗粒表面润湿性薄膜，也为研究注气过程中化学参数对流体分布的影响开辟了新的途径.
环境扫描电镜技术的样品室中存在的气体阻止了表面电荷的产生，观测样品表面不需要添加导电物质，从而减少了人为因素对观测结果的影响，可以对含油样品在环境状态和低真空下进行动态实验，能最大程度地反映样品的原始状态.Combes等(1998)利用环境扫描电镜对烧结玻璃模型、砂岩、碳酸盐岩渗吸过程观察，认为该技术能够较为准确的反映孔隙中润湿性的非均质性，解释矿物组分、尺寸、分布等因素对流体物理化学行为的影响.为了更深入的认识剩余油分布及提高采收率，矿物微观表面粗糙度和润湿性的测量也尤为重要，Kareem (2015)结合原子力显微镜和环境扫描电镜技术，通过驱替实验研究了不同石英晶体粗糙面条件下润湿性的改变，发现当矿物表面粗糙时更容易形成油湿.
冷冻扫描电镜技术需要在样品表面添加涂层(比如碳、铬、金等)以增加导电性，而环境扫描电镜技术在研究动态实验过程中需要一些主观的解释，针对这些缺点，Boassen等(2006)提出一种融合冷冻扫描电镜技术和环境扫描电镜技术的方法，使得样品不需要添加涂层而且解释相对容易，通过实验观测认为该方法可以得到满意的效果和充足的信息.但扫描电镜技术只能产生样品局部表面的高分辨率图像，并不能得到较大尺度的三维空间剩余油分布，限制了其在微观剩余油研究的广泛应用.
1.3 激光共聚焦技术激光共聚焦技术(Confocal Laser Scanning Microscopy，CLSM)属于显微镜分析技术范畴，但弥补了光学显微镜与扫描电镜的不足，不仅可以获得高清晰度、高分辨率的图像，还可以进行三维立体图像重建.它基本方法是激光束射到样品的一个细微点上，避免了非照射区域产生光散射，使发射光信号通过检测针孔到达监测器.入射光源针对针孔和检测针孔的位置相对于物镜焦平面是共轭的，来自焦平面的光可以通过检测针孔被检测到，而来自焦平面上、下的光被阻挡在针孔两边，对样品进行二维图像采集，再经过处理可以得到样品的二维和三维图像.
激光共聚焦技术不仅在储层孔隙、微裂缝、有机质分布等研究应用广泛，也可以对岩石含油性、微观剩余油分布规律进行研究.白振强等(2013)通过分析激光共聚焦扫描激发荧光剩余油分布微观三维图像，得到聚合物驱替过程中微观孔隙中原油轻质和重质组分变化规律及微观分布特征，并研究了剩余油分布与粘土矿物之间的关系.Krummel等(2013)使用共聚焦技术直接对玻璃球模型进行驱替实验的三维可视化分析，在水驱油过程中，润湿相在狭窄裂隙中挤压油，形成了不连通的俘获油滴.构建油滴的三维结构模型，发现随着毛细管数(Ca)的增大，剩余油总量减小，这个现象说明了润湿相流体黏滞压力与毛管压力的相互作用可以促使油滴通过孔隙.由于该技术的空间分辨率是几百个纳米，可以直观的反映微小孔隙内和岩石粗糙表面上剩余油结构形态，但是仅能对一定厚度岩心样品的孔隙成像，规模很小，应用并不广泛.
1.4 CT技术CT技术也被称作计算机断层扫描技术(Computed Tomography，CT)，它是通过计算机图像重建，使指定层面上不同密度的材料信息以高分辨率的数字图像显示出来(于艳梅等，2012).根据射线光源的选择不同，CT技术可以分为X射线CT技术(X-ray tomographic microscopy)，中子CT技术(Neutron Computed Tomography)，超声CT技术(Ultrasonic Computed Tomography)，同步辐射CT技术(Synchrotron Computed Tomography)等.目前，在岩心微观剩余油研究中，最常用的是X射线CT技术.
X射线穿透物体时受到核外电子作用而被衰减，其质量衰减系数与材料原子序数有确定的函数关系，能在被测物体无损的条件下测量、计算被测物体断面的密度分布，以计算机断面图像的形式并通过灰度变化反映出断面材质的密度变化，直观准确地呈现出被测物体内部的结构关系、尺寸、材质的组成变化情况，因而被广泛应用到石油地质分析的多个领域(张敏和孙明霞，2002).它基本方法是将X射线信号源绕着样品旋转，对一个固定的横剖面在不同角度测量穿过样品的X射线的强度，使CT能在单个横剖面上对线性衰减系数进行空间析像.根据强度资料重新构建样品的二维横剖面，把一系列横剖面叠加在一起就形成样品的三维图像.重建的图像能够给出每个像素X射线的衰减系数，经过转换通常用CT值表示.由于不同流体线性衰减系数μ数值不同，因此可以根据这种差异区分不同流体信号，研究微观剩余油分布.
早在1982年，Wang等人尝试将X射线CT技术应用于岩心驱替过程，获得了横截面上剩余油分布图像.随后，Hove等人(1987)利用该技术对北海砂岩驱替过程可视化成像，由于受到技术的限制，成像效果仅为沿柱塞岩心轴向的二维图像，而且空间分辨率很低，为500 μm.尽管成像效果较差，但该方法为剩余油形成机理、分布和形态的研究提供了新的方向.随着空间分辨率的不断提高，成像质量的不断完善，X射线CT技术在微观剩余油分布研究中已经有了长足的发展.1998年，Coles利用该技术先对干燥砂岩进行扫描成像，再将其饱和油后进行水驱，扫描得到水驱后微观剩余油分布图像.实验过程采用的空间分辨率为30 μm，每张图像成像时间大约30 min，已经可以实现剩余油分布的三维可视化.Turner等人(2004)对气驱油后的等径堆积玻璃球模型和Berea砂岩进行成像(空间分辨率9 μm)，并对微观剩余油赋存形态进行了描述.Prodanovic′等(2006)通过注入固定体积的润湿性流体和非润湿性流体，确保驱替和渗吸过程达到不同的阶段，并进行成像分析(空间分辨率4.1 μm).为了保证驱替过程不同阶段成像质量，最重要的是保持流体在空间分布的稳定性，Youssef等(2010)结合X射线CT技术与针对驱替实验的流动微单元，设计新的实验方法固定流体分布，并描述了不同驱替阶段下孔隙几何形态与界面张力对微观剩余油的影响.Iglauer等(2012)研究了水湿砂岩和油湿砂岩经过水驱后簇状剩余油大小分布、形状和赋存量的不同，发现在其他条件相近时水湿砂岩残余油饱和度高于油湿砂岩残余油饱和度.在三维空间中，油湿砂岩残余油呈粗糙、扁平状附着在孔隙表面，或者集中在小孔隙中；水湿砂岩残余油基本存在于大孔隙的中央位置.Iglauer等(2014)又进行了三相流驱替后的CT扫描实验，它包括两类过程：一类是气直接注入原始油藏；另一类是先水驱原始油藏，再将气注入.通过对两类过程剩余油分布的比较，发现第一类过程具有较高的采收率.Bultreys等人(2015)提出了一种基于X射线CT扫描动态测量天然砂岩驱替过程中内部流体分布的新方法，每隔12 s进行一次三维可视化成像，如此短的时间间隔可以较好的显示驱替过程中流体的三维空间变化，甚至能够研究海恩斯跳跃等特殊现象.
国内应用X射线CT技术研究岩心内微观剩余油分布起步较晚，1997年胜利油田地质科学研究院引进带有微焦点X射线管的ACTIS工业CT系统(盛强等，2005)，利用驱替装置按照事先设计好的驱替流程进行驱替，同时利用CT扫描设备对岩心中的岩石骨架、流体分布进行扫描，进而获取不同驱替阶段过程岩心内部的油水分布信息，再现真实润湿状态下的油水运动规律和驱替特征.中国石油勘探开发研究院于2007年引入并研制了CT扫描岩心驱替系统(高建等，2009)，可以对岩心驱替过程进行在线CT扫描，同时采集驱替过程中的进出口流量、压力.由于实验中的油水CT值相差较小，需要添加CT增强剂区分油水信号，常用的方法是在水中加入1%~10%的碘化纳或溴化钠、碘化钾，以增强水的CT值，比如加入8%碘化钠CT增强剂可使地层水CT值为1000，而驱替油CT值为214，较好的区分油水信号；另一种方法是在油中加入5%的碘代己烷以增强油的CT值，进而区分油水在孔隙中分布.应用X射线CT技术对驱替过程可视化分析，不仅可以直观表征微观剩余油分布，还能够深入了解采油机理、监测流体分散与窜流特性、认识聚合物驱对提高波及效率影响和揭示地层伤害机理等.通过对砾岩、低渗透和特低渗透岩心水驱过程的在线CT扫描监测，可以得到驱替过程中岩心含油饱和度沿程分布信息和CT重构切面图像，进而分析砂砾岩、低渗透和特低渗透岩心的驱油机理(王瑞飞等，2008；吕伟峰等，2013；邓世冠等，2014).张顺康等(2006)利用指示克里金方法对微观CT层析图像中的剩余油图像和孔隙图像进行了分割，得到了二维CT层析图像中的孔隙以及油水分布，基于Marching Cube算法实现了不同驱替阶段中剩余油分布的三维可视化.侯健等(2014)借助CT技术，研究水驱不同驱替时刻的三维空间剩余油分布信息，建立了微观剩余油的赋存状态表征方法，对不同驱替时刻剩余油在岩心模型中微观分布状况进行定量表征. 曹永娜(2015)详细介绍了水驱与聚合物驱过程中CT扫描实验流程，研究不同驱替方式下剩余油空间分布状态，并定量分析含油饱和度沿程分布变化. Wang等人(2015)以水湿砂岩为样品，分别在水驱、聚合物驱和凝胶颗粒驱替实验条件下进行不同阶段的CT扫描，通过图像处理技术得到三维微观剩余油空间分布，进而直观认识不同驱替方式对微观剩余油的影响. X射线CT技术不仅可以研究常规微观剩余油气在孔隙中的分布特征，而且能够直接观测甲烷水合物的形成与空间分布. 胡高伟等(2014)研制了适用于沉积物中水合物微观赋存状态的CT原位探测装置，通过模拟天然气水合物生成/分解过程，直接观测水合物在沉积孔隙中的微观分布特征，为水合物形成机理和空间分布特征提供重要依据. X射线CT技术具有快速、无损和较高分辨率的优势，而且该技术目前已经高度发展和商业化，因此在微观剩余油研究中应用广泛.
除了X射线CT技术，也有学者利用中子CT技术研究微观剩余油特性.中子CT技术也是基于射线穿过物体时会发生衰减的基本原理，但与X射线CT技术相比不同之处是放射源为中子，由于中子不带电，能轻易的穿透电子层，与原子核发生核反应(裴宇阳等，2004).中子射线可以穿透大部分重元素物质，而且非常容易被氢元素吸收，对含氢物质敏感.由于油、水均为含氢物质，在实验室研究岩心微观剩余油时，需要将水替换为重水(D2O)，相对于X射线CT扫描技术可以更好的区分骨架、油、水等不同物质的空间分布，而且较好的生成流体流动过程中实时的高质量图像.De Beer和Middleton (2006)模拟在高压情况下的驱替过程，先将砂岩饱和油，再用重水进行驱油，利用中子CT技术对驱替过程进行成像分析，进而研究不同驱替阶段岩心内油水的分布情况.由于中子射线要求射线源有很高的稳定性，限制了中子CT技术的发展，目前对于微观剩余油的研究也只是起始阶段.
1.5 核磁共振成像技术核磁共振成像技术(Nuclear Magnetic Resonance Imaging，NMRI)是利用三轴梯度场技术对氢核进行空间扫描获得氢核空间分布的可视化技术，它的基本原理是基于岩石内部油、水的自旋原子核(如1H，13C)在磁场中产生共振并具有弛豫作用的特性进行岩心分析.在静磁场中叠加一个梯度磁场，建立核磁共振信号的共振频率与原子核所在位置的关系，再利用快速傅立叶变换、图像重建等技术得到各点T1、T2信号，从而获取核磁共振图像.三轴梯度线圈的设计使得任意方向的数据信息获取成为可能，比如一维剖面，二维切片甚至三维体数据.它既可以获取任选岩心切面内mm级甚至μm级分辨率下的流体分布特性，也可以直观描述岩石孔隙结构特性.由于油水的分子结构不同，而且水分子的1H弛豫率比油分子的1H弛豫率快，当水溶性Mn2+添加到水中时，水分子的1H弛豫率更强而油分子1H弛豫率不变.因此利用核磁共振成像技术可以区分油水信号，从而定量研究微观剩余油赋存量和空间分布规律.
早在1986年，Baldwin等人采用改进的医学核磁共振成像设备成功的对水湿和部分油湿Berea砂岩水驱油实验进行核磁共振成像实验，发明了在水中加入含有顺磁性离子Mn2+减少水相弛豫时间的方法，从而区分油水两相并得到岩心内剩余油分布图像.随后，Mandava等人(1990)也利用核磁共振成像方法观测岩心内两相流流动，定量得到孔隙度与饱和度剖面.由于核磁共振成像可以在任选的切面内，以毫米或者微米级的分辨率显示岩石结构或流体捕集及流体特性，这种优势吸引了更多学者进行复杂岩心内流体分布规律(Chardaire-Rivière and Roussel, 1990)、速度成像(Chang and Watson, 1999)、提高采收率机理(Romero-Zerón et al., 2009)等研究.Ersland等(2010)对含有裂缝的碳酸盐岩水驱过程进行核磁共振成像，快速获得高质量的油水分布信息，发现裂缝的存在显著提高了孤立基质岩块末端的采收率.对于核磁共振成像质量而言，高场强测量意味着获得高分辨率、高信噪比的图像.但对于天然岩心而言，内部骨架复杂成分中大多含有铁磁性物质，会对较高场强的磁场产生较大的破坏，导致成像图像重影、伪影的出现，因此在进行微观剩余油分布研究的过程中首先需要考虑岩心样品的适用性.例如，对花岗岩驱替实验过程进行核磁共振成像时(Chen and Kinzelbach, 2002)，由于花岗岩具有顺磁特性，为了获得稳定的微观剩余油空间分布图像，不能进行高场强的核磁共振成像，只能使用低场核磁共振成像.Mitchell等(2013)详细介绍了核磁共振成像技术在岩石内部流体识别领域的应用，并认为随着核磁共振仪器硬件水平的发展，高分辨率监测提高采收率过程中微观剩余油三维空间分布变化将有巨大的潜力.
国内早期已经引入核磁共振成像设备并对核磁共振成像理论进行研究，但主要局限于医学领域.80年代后期，中国科学院武汉物理研究所引进Bruker MSL-400超导NMR谱仪，肖立志(1995)利用该设备微成像附件对饱和煤油和饱和盐水的三类岩心进行成像，认为核磁共振成像技术可以对油气储集层岩石的孔隙性、渗透性和含油性及其空间分布进行定量描述，对水驱油、聚合物驱油等强化采油措施进行微观评价.1992年，北京石油勘探开发研究院万庄分院从瑞士引进磁场强度为4.7T、孔径为40 cm的核磁共振成像设备(BMT 47/40)，将核磁共振成像技术应用在岩心分析(卢祥国等，1996).利用该套设备，许多学者开展了聚合物驱油、三元复合驱油、热采驱油、周期注水驱油等实验(黄延章等, 1995, 王为民等，1997)，研究了聚合物调剖渗流机理，对多种驱油方式、剩余油分布产生新的认识.通过实验分析，王为民等(1997)认为核磁共振成像技术是模拟地层条件的可视化驱替实验研究的很好手段，能够对指进、波及效率、扫油面积、原油捕集等过程进行详尽的分析研究.
随着核磁共振成像仪器的快速发展和国内核磁共振成像设备的增加，对于利用核磁共振成像技术进行微观剩余油的研究也逐渐增长.利用核磁共振成像技术，水驱油(气)实验已经应用于砂砾岩岩心、河流相非均质储层岩心、煤层岩心、人造胶结岩心、人造填砂岩心等不同类型岩心(鲁国甫等，2001；潘一山等，2008；葛丽珍等，2014；Cheng et al., 2015；Teng et al., 2016)，得到不同的微观剩余油(气)分布规律.邓瑞健(2002)结合低场核磁共振T2弛豫测量和高场核磁共振成像(SPEC-3000型)，在水驱油的实验研究中定量地对岩心的初始饱和油分布、水驱油渗流及剩余油分布特征等进行了分析，从微观上揭示剩余油在毛细孔道中的分布特征，认为该技术实验结果能够指导实际油藏剩余油的挖潜方向.郎东江等(2006)利用核磁共振成像技术研究均质、非均质模型中复合剂驱油动态实验过程，得到了剩余油分布及复合剂驱油段塞前缘油墙的形成，认为该技术不仅可以直观地观察复合剂驱油作用时不同流体的空间分布情况，而且可以对复合剂不同配方的适配性作出评价.赵越超等(2012)总结了高温高压条件下驱替实验与核磁共振成像联测所用夹持器的发展现状，并在高压条件下多孔介质内CO2驱油进行可视化研究.
由于核磁共振成像技术仅探测岩石中流体信息，反映液相在岩石中的分布和聚集及液相与周围岩石环境的界面效应等情况，不涉及岩石基质等固体物质，该技术已经在流体分析的多个领域广泛发展.尽管在微观剩余油研究中应用核磁共振成像技术的数量相比利用X射线CT技术少，但核磁共振成像技术具有能够在相对较大尺度上研究剩余油分布和驱替过程中快速检测流体分布变化的巨大优势，使得该技术在研究微观剩余油领域依然是非常重要的手段.
2 数值模拟方法研究现状微观多孔介质物理属性模拟的基本条件需要基于多孔介质模型，而且数值模拟的正确性及适用性在很大程度上取决于所建立的岩石微观模型，只有模型的孔隙结构能够反映真实岩样的孔隙结构特征时，模拟结果才具有理论及应用价值.多孔介质模型的发展过程可分为毛管模型、随机孔隙网络模型、数值重建三维数字岩心、真实重建三维数字岩心等，目前学者构建多孔介质三维模型的主要方法有X射线CT法、FIB-SEM、顺序指示模拟(SIS)、多点地质统计学(MPS)、过程法、孔隙网络模型等(姚军等，2005；刘学锋等，2013).核磁共振成像方法也可以重建三维数字岩心模型(Minto，2014)，尽管空间分辨率低，但对于砾岩等粒径大、岩心尺寸大的岩心三维重建效果好.只有多孔介质模型重建完成后，才能进行微观剩余油数值模拟.
微观剩余油数值模拟属于流体力学中多孔介质、多相流流动模拟的范畴.随着渗流机理的逐步认识加深与计算机处理能力的大幅提升，已经出现了多种微观剩余油的数值模拟方法，根据模拟的算法类型不同，这些数值模拟方法主要包括基于网格化方法、无网格化方法、孔隙网络方法等.
尽管国内研究微观多孔介质模型的学者很多，但进一步研究微观剩余油的团队相对较少.目前国内主流方法是无网格化方法中的格子玻尔兹曼方法及孔隙网络方法，研究团队包括中国石油大学(北京)(陶果团队)、中国石油大学(华东)(姚军团队、孙建孟团队、侯健团队、徐守余团队等)、大连理工大学(宋永臣团队)等，通过改变岩石的润湿性、孔隙结构及流体性质等影响因素研究微观剩余油分布的形成机理、赋存量及空间分布等特征，丰富了微观剩余油数值模拟领域的成果.
2.1 基于网格化方法基于网格化方法(Grid-Based Methods)通常被认为是计算流体动力学(CFD)中解Navier-Stokes方程的常规方法，例如有限差分法、有限单元法和有限体积法，它们具有良好的数值模拟计算效率和出色的流体流动模拟能力.有限差分法(Finite Difference Method, FDM)将求解域划分为差分网格，用有限网格节点代替连续的求解域，以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组.有限单元法(Finite Element Method，FEM)把计算域划分为有限互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解.有限体积法(Finite Volume Method，FVM)将计算区域划分为一系列不重复的控制体积，并使每个网格点周围有一个控制体积，将待解的微分方程对每一个控制体积积分，得出一组离散方程.对于这三种网格化方法，FDM尽管易于并行编程，但由于只考虑网格点上的数值而不考虑值在网格点之间如何变化，对不规则复杂区域连续性要求高；FEM适合处理复杂区域，具有更高的自由度，计算精度高，但计算所需的内存和计算量大；FVM对复杂区域的适应性好，但计算精度相对较低.目前，对复杂多孔介质微观剩余油研究中，FEM、FVM应用较多，而FDM应用较少.
对于基于网格化方法在多孔介质两相流模拟中最大的困难是界面运动的处理.多孔介质两相流最大的特征是存在饱和度变化的相变锋面，而且对于非均质多孔介质饱和度可能会突然变化，导致流体界面变化剧烈.为了描述流体的流动状态，计算网格的尺寸必须足够精细以满足计算精度的要求，如果对计算区域采用均匀网格计算，则计算量非常大.为了确定界面位置，根据网格类型目前主要发展了三类方法(Raeini et al., 2012)：移动网格法、表面追踪法和体积追踪法.移动网格法，两相流体分别网格化，流体界面为两相网格的相交边界.而表面追踪法和体积追踪法则是对两相流体采用相同的网格化方式，界面位置用高度函数或者标记质点确定.移动网格法和表面追踪法可以很好的定义两相流体锋利的界面，但对于大的界面运动，需要重新定位，增加了计算的复杂程度.体积追踪法采用无质量粒子或指示函数对两相流界面进行标记，处理效果更好，在微观剩余油研究领域继承发展了流体体积法、水平集方法和相场模型法.
流体体积法(Volume of Fluid，VOF)广泛应用于求解两相或多相流存在而且界面处于运动状态的复杂情况，该方法首先由Nicols等人(1981)提出，该方法通过计算含标记函数的输运方程达到追踪界面的目的，在多孔介质多相流模拟算法中追踪油水界面的形状.由于VOF法追踪的是网格中的流体体积，而不是追踪流体质量的运动，所以容易实现，追踪的界面锐利性好.商业软件ANSYS Fluent®结合FVM方法和VOF可以很好的模拟两相流流动及剩余油分布情况(Piller et al., 2014)，该软件模拟剩余油分布还包含以下方法：利用分段线性插值算法追踪界面；采用连续表面张力(CSF)方法作为表面张力模型方法，在动量方程中添加源项；基于流体与固体之间的接触角模型计算黏附力；采用压力隐式算子分裂方法实现压力速度耦合；使用单压力校正方法提高计算效率.
水平集方法(Level Set Method，LSM)在两相流模拟中能够捕捉两相界面的运动，该方法首先由Osher和Sethian (1988)提出，通过定义符号距离函数，相界面可以由距离函数的零等值线确定.水平集方法不需要显式地追踪运动界面，避免了对拓扑结构变化的处理，计算稳定.Prodanovic′等(2006)发展了基于水平集方法的逐步准静态算法，能够对驱替过程逐步数值模拟，认为该方法在多孔介质两相流中有较好的模拟效果.Shabro等人(2010)利用该方法在枫丹白露砂岩和白云岩三维数字岩心数据上模拟驱替过程，得到微观剩余油分布与相对渗透率曲线，认为该方法得到的非润湿相捕获与实验观测结果一致.Peyman和Apostolos (2016)针对多种矿物组分、不同润湿性的微观多孔介质模型，提出一种结合水平集方法与流体体积法求解多相流流动的稳定算法，并研究驱替过程中润湿性的各向异性对非润湿相流体三维空间分布变化的影响.
相场模型法(Phase Field Method，PFM)也被称作扩散界面法(Diffuse Interface Method，DIM)，是van der Waals在使用随着界面变化的密度函数模拟液-气系统时提出(Rowlinson，1979)，通过相场变量φ得到界面层的信息，而不是直接追踪两种流体界面的变化，利用4阶偏微分Cahn-Hilliard方程控制相场变量的演化.商业软件COMSOL MultiphysicsTM结合FEM方法和PFM能够模拟两相流流动及剩余油分布情况(Bogdanov et al., 2011)，由于明确描述了所有重要的物理机理，模拟的结果符合实际.Amiri和Hamouda (2013)比较了LSM和PFM在多孔介质两相流模拟的效果，认为两种方法均能实现基本的流体流动过程，但PFM在复杂多孔介质中对物理现象的细节描述更加成功.
2.2 无网格化方法无网格化方法(Mesh Free Methods)是一种完全基于点近似的新兴数值方法，由于该方法不需要节点间的拓扑连接关系，完全避免了传统方法在求解时不断重分网格所带来的精度降低及计算复杂性增加等问题，在两相流流动模拟有独特的优势.无网格法具有鲜明的技术特点(张雄等，2009)：它采用基于点的近似原理，只需要节点信息而不需要将节点连成单元，与网格算法相比，可以彻底或部分地消除网格，不需要网格重构.其基本思想是在计算域上用一些离散的点来拟合场函数，从而摆脱了网格单元的限制.常规的无网格化方法包括：光滑质点流体动力学方法、分子动力学方法、耗散粒子动力学方法、移动粒子半隐式方法、格子玻尔兹曼方法等(Shakibaeinia and Jin, 2012).在多孔介质微观剩余油数值模拟中，应用最广泛的是格子玻尔兹曼方法，另外也有部分学者应用光滑质点流体动力学方法进行研究.
光滑质点流体动力学方法(Smoothed Particle Hydrodynamics，SPH)是一种单纯的Lagrangian无网格算法，由Lucy人(1977)在研究三维无边界天体物理问题时提出，它用一系列任意分布的粒子质点来代表整个连续介质流体并估计相应的偏微分方程，可以容易地追踪自由表面、移动界面以及物质的交界面.Tartakovsky等(2007)综述了SPH在多相流模拟领域的发展和优势，认为该方法在高流体密度比和高黏度比条件下模拟依然稳定，并对各向异性多孔介质的驱替过程进行模拟.Zhou等(2010)应用SPH研究了缝洞组合结构中的油水两相驱替机理，考察注水速度、壁面润湿性、裂缝宽度等一系列因素对两相流动行为特征和油的最终采收率的影响，模拟结果与相关物理实验基本相符.Kunz等人(2016)制造了二维微观仿真模型，通过驱替实验和SPH模拟结果对比，发现两者在流体流通通道和毛管压力曲线结果吻合非常好，认为SPH对于多相流体模型表面张力的表征是可行的.尽管该方法在两相流模拟方面已经有较为成熟的发展，但在多孔介质模型中的应用依然较少.
格子玻尔兹曼方法(Lattice Boltzmann Method，LBM)是20世纪80年代中期发展起来的一种全新的流场计算方法.它源于Frisch等人(1986)对格子气方法(LGA)的改进.不同于传统的流体计算方法，LBM并非采用离散化技术得到宏观连续微分方程或积分方程的数值解，而是基于微观尺度上的统计力学的Boltzmann方程，通过描述粒子的相互作用，在宏观上达到模拟复杂物理现象的目的(李学民等，2004).该方法对边界条件的处理方法比较简单，可以处理复杂的边界，已经在流体力学多个领域应用.对于多相流模拟最常用的格子玻尔兹曼模型主要有三种：颜色模型，SC模型和自由能模型.最早的多相流格子玻尔兹曼是颜色模型(Color Model)，由Gunstensen等人(1991)基于多相流格子气模型提出，该模型用红蓝两种粒子分布函数来代表不同的流体，通过“再着色”步骤实现相的分离.Ferréol和Rothman (1995)在X射线CT技术重构的枫丹白露砂岩三维孔隙结构空间的基础上应用颜色模型实现了两相流渗吸驱替过程的模拟，得到不同驱替时刻润湿相与非润湿相的三维空间分布图.Ramstad等(2012)对Bentheimer砂岩三维数字岩心进行多相流模拟，同时计算稳态和非稳态条件下相渗曲线，并结合实验结果分析相渗曲线影响因素.Jiang等人(2014)构建分辨率为3.2 μm的Berea砂岩三维数字岩心，利用GPU并行的颜色模型研究表面张力对孔隙中多相流体分布的影响，高效计算驱替过程中流体空间位置和相渗曲线，发现随着表面张力增加，相对渗透率减小.颜色模型的不足在于“再着色”步骤引起的扰动会带来伪各向异性表面张力、相界面处附近产生非物理现象，而且该模型没有考虑热动力学对复杂流体的影响.SC模型(Shan-Chen Model)是目前格子玻尔兹曼方法在多相流模拟中应用最广泛的模型，也被称作伪势模型，最初由Shan和Chen (1993)提出，该模型引入粒子间非局域相互作用势模拟两相流系统，通过流体粒子与流体粒子之间的内聚力及流体粒子与固体粒子之间的黏附力模拟表面张力、润湿性等基本界面现象，能够模拟所有的接触角.由于LBM在多孔介质两相流模拟计算量较大，Pan等人(2004)对SC模型进行优化，同时实现了两相流流渗吸驱替过程的并行化模拟，提高了计算的速度而且模拟得到的毛管压力滞后曲线与实验结果非常温吻合.Li等人(2005)利用SC模型研究了多孔介质多相流，计算得到的相对渗透率曲线与实验结果一致.Sukop等人(2008)对人造石英砂岩进行油水驱替过程的X射线CT扫描，同时采用SC模型进行两相流数值模拟，通过比较两者微观剩余油分布图像和沿Z轴定量饱和度值，揭示了LBM在描述多孔介质中多相流体分布的潜力.陈萍等(2013)利用SC模型研究孔喉道粗糙性对油水两相流动的影响，比较不同驱替模式下光滑和粗糙模型的含水饱和度等变化，认为孔喉表面粗糙性对油水流动起阻碍作用，进而影响微观剩余油的赋存量和空间分布.与颜色模型相比，SC模型稳定性相对较好，但同样没有基于热动力学的合理假设，只能应用于等温的多相流流动模拟，而且对于高密度比、高黏度比的流体应用效果差.自由能模型(Free Energy Approach)是Swift等人(1995)直接从多相流体的自由能理论出发构造的与热动力学理论一致的多相流体的格子玻尔兹曼模型，该模型根据自由能函数构造格子玻尔兹曼模型的平衡态分布函数，通过引入一个非理想流体的热力学压力张量，使得系统的总能(包括动能、内能和表面能)守恒得以满足.Huang (2011)在多孔介质中改进了自由能模型模拟两相流流动，使模型满足了伽利略不变性，而且可以准确模拟相对较大的黏度比的两相流，并将该方法与颜色模型、SC模型进行对比，认为该方法在两相流模拟的多个方面更有优势，但目前在三维复杂各向异性多孔介质中应用依然较少.另外，基于动力学理论的格子玻尔兹曼模型(He et al., 1999)也能够合理的模拟两相流流动，模型中碰撞项对局部分子间的相互作用力用Enskog理论表示，采用平均场近似的方法计算大范围的分子间相互作用力，同时考虑了分子间的吸引力和排斥力，从而更加的符合客观物理规律，但目前在三维空间多孔介质两相流模拟中缺乏相关研究.
2.3 孔隙网络方法孔隙网络方法(Pore Network Modelling，PNM)是多孔介质微观多相流模拟中发展较为成熟的一种方法，已经在石油领域和环境领域广泛应用.孔隙网络方法不仅可以求取各种条件下相渗曲线和毛管压力曲线，也可以显示微观驱替过程和微观剩余油分布.从模拟流体流动的物理机理划分，目前常用的孔隙网络方法可以分为准静态孔隙网络方法和动态孔隙网络方法(Zhang et al., 2009).准静态孔隙网络方法假设毛管压力控制着整个流动过程，不考虑黏性压降，每次运算只有一个孔隙单元或喉道单元的流体发生变化，只适用于层流模拟.动态孔隙网络方法同时考虑毛管压力和黏性压降对流体驱替的贡献，每次运算有若干个孔隙单元或喉道单元的流体发生变化，适用于界面张力非常低的两相驱替过程、含有聚合物和泡沫等高压力降的流体流动模拟.由于多相流的流动模拟的对象是孔隙网络模型，根据实验条件和研究目的不同，孔隙网络模型可以分为规则孔隙网络模型和基于图像的孔隙网络模型，尽管模型的形式不同，但模拟所用的方法基本是通用的.
规则孔隙网络模型最早由Fatt于1956年提出，他构建一个简单规则的二维毛细管规则网格模拟驱替过程，根据Young-Laplace方程确定流体进入不同半径毛管的顺序，应用等效的电阻网络测量毛管压力曲线和相渗曲线，由于模拟结果与实验结果基本吻合，因此他认为网络模型是一种可以表示真实多孔介质孔隙空间的有效模型.基于他的成果，许多学者在孔隙网络的多个方面开展研究.Chatzis和Dullien人(1977)将二维网络模型扩展到三维，并用等效的三维电阻网络进行了多相流模拟计算.由于实际孔隙空间几何形态变化多样，孔隙网络模型用规则几何形状的孔隙和喉道单元连接的系统进行表征，孔隙和喉道尺寸的大小按研究需要进行赋值.Fenwick和Blunt (1998)采用截断威布尔分布构建随机网络模型，模拟水湿条件下多孔介质油气水三相微观空间分布.胡雪涛和李允(2000)构建具有对数正态分布孔喉尺寸的随机网络模型，利用定向渗流理论模拟驱替过程，得到剩余油分布特征.另外，孔隙吼道尺寸还可以用正态分布、Haring-Greenkorn分布、压汞曲线等进行表征(王克文等，2005).逾渗理论是统计物理中描述复杂多孔介质中渗流及传输规律的方法，模拟初始时刻孔隙网络中完全饱和水，逐渐增大注入流体(油)的压力，当油水两相压力差克服孔隙体或喉道渗流阻力时，孔隙中的水被油驱替(王克文等，2006)，该理论也是国内孔隙网络模型模拟由毛管力主导准静态驱替的常用方法. 王克文等(2006)将微观剩余油分为孤立状、条带状和网络状三种不同的分布形态，利用网络模型研究了不同储层特性对聚合物驱后剩余油含量及微观分布形态的影响.Pereira等(1997)模拟油气水三相流体驱替过程，得到不同驱替条件下微观剩余油分布状态.杨永飞等(2010)基于英国赫瑞瓦特大学的三维规则网络模型模拟研究了六种不同润湿性油藏注气过程，分析油藏岩石润湿性对气驱剩余油微观分布的影响机制.
为了能够反映岩石真实的孔隙结构，逐步发展了基于图像的孔隙网络模型.这种模型与基于网格化、无网格化等多孔介质两相流模拟方法一样，数据基础都是重建的多孔介质三维模型，但不同之处在于孔隙网络模型需要提取具有三维多孔介质拓扑结构和几何特征的孔隙网络，既能保持较高的运算效率也能准确模拟多孔介质渗流性质.提取孔隙网络的方法主要包括多向扫描法、中轴线法、Voronoi多面体法和最大球法等(Silin and Patzek, 2006).目前最常用最大球法，Dong和Blunt人(2009)改进最大球法对孔隙网络进行提取，分析了孔喉尺寸、形状因子、配位数等参数，可以准确描述孔喉的尺寸和孔喉连接的拓扑结构.刘向君等人(2014)利用数学形态学算法对CT扫描数字岩心进行孔隙结构量化及表征，统计获取孔隙度、孔隙体积分布及孔径分布特征，建立了等价孔隙网络模型. 张顺康(2007)对天然岩心和人造石英砂岩进行水驱和聚合物驱实验，实验过程中进行X射线CT扫描，得到真实条件下微观剩余油三维空间分布，再对重构的数字岩心提取孔隙网络，结合网络模拟研究剩余油分布规律及成因，通过比较发现孔隙网络方法计算得到的微观剩余油分布与真实驱替过程微观剩余油分布基本一致.Hou等(2011)利用基于图像的孔隙网络模型研究润湿性、孔隙半径、形状因子、孔喉比、配位数等在水驱过程中对剩余油分布的影响，提出模糊综合评判法预测微观剩余油三维空间分布，与常规孔隙网络方法得到的结果一致.Ryazanov等(2009, 2014)利用CT扫描的不同Berea砂岩样品，提取三维孔隙网络结构进行驱替模拟，定量分析剩余油饱和度与平均接触角的关系，之后，研究了不同润湿性条件下驱替过程中微观剩余油三维空间结构的变化，得到不同孔隙结构中剩余油分布的配置关系.Al-Dhahli等(2014)利用提取的孔隙网络模型模拟在饱和油时水驱、气驱和水-气交替驱等两相流、三相流流动，比较不同驱替方式下剩余油饱和度和相渗曲线.任奕明等(2015)利用CT扫描技术建立岩心的三维映射网络模型，并对其进行水驱油模拟，基于微观驱油效率和含残余油孔隙比例随各因素的变化关系，研究了孔隙中残余油形成的影响因素及机理.
2.4 其他模拟方法除了以上三类最常用的多孔介质两相流模拟方法，还有两类方法可以研究多孔介质微观剩余油.一类是形态学方法，另一类是最大内切球方法.尽管对于这两类方法研究相对较少，但两类方法可以较快的近似获得多孔介质孔隙不同流体的空间分布情况.
形态学方法(Pore Morphology Method，PMM)是在多孔介质图像上直接计算的方法，可以在复杂孔隙结构条件下有效的计算准静态两相流驱替过程，得到微观剩余油分布及毛管压力曲线信息.Hilpert与Miller (2001)首先将数学形态法的思想融入多孔介质驱替模拟，腐蚀和膨胀运算是该方法的数学基础.初始状态孔隙空间充满润湿相，随着压力的增加，不同大小的孔隙从非润湿相注入的一端开始进行腐蚀膨胀运算，直至驱替结束.Liu等(2009)利用形态学方法模拟X射线CT扫描的数字岩心孔隙空间内流体分布，并结合有限元方法计算岩心电性变化规律.该方法具有较高的计算效率，但只能模拟完全水湿或完全油湿条件下流体空间分布.
最大内切球方法(Maximal Inscribed Spheres，MIS)最初在多孔介质中应用是提取孔隙网络，Silin和Patzek (2006)将该方法应用于两相流模拟，认为润湿相流体分布在大孔隙边缘和小孔隙中，非润湿相分布在侵入孔隙中央，并假设两种流体的界面近似球形，通过调整参数计算得到的毛管压力曲线与实验结果一致.随后，Silin等(2011)将该方法应用到CO2注入Frio砂岩和Berea砂岩的两相流分布，并与X射线CT扫描气驱实验结果进行对比，认为最大内切球方法适用于多孔介质两相流分布的模拟.该方法同样具有较高的计算效率，但没有考虑润湿性对流体分布的影响.
3 展望 3.1微观剩余油形成机理、赋存量以及空间分布的研究不仅可以改善油田开采工艺，提高采收率，而且对岩石物理研究和测井精细评价有重要的指导作用.不同微观剩余油赋存量和分布特性对岩石物理属性有较大的影响，如声波波速和Q值随含油饱和度的变化而改变(李爱兵等，1995)，岩石弹性参数受剩余油气赋存量的影响(陈信平等，2013)，核磁共振信号包含剩余油气饱和度信息(赵文杰和谭茂金，2008)，岩石电学特性通常与油水饱和度满足阿尔奇现象(孔强夫等，2015)等.这些由于微观剩余油的变化导致的岩石物理性质的改变对声、电、核磁等测井资料产生影响，因此微观剩余油也是未来较大尺度剩余油研究、精细剩余油层划分以及提高剩余油测井评价效果的研究重点.
3.2传统光学技术不能全面反映三维空间真实岩心特征，但由于实验制作方便，成本相对较低，该技术依然是油田不同层位微观剩余油特性研究主要手段.扫描电镜技术与激光共聚焦技术实验成本昂贵，视野较窄，不是研究微观剩余油的通用手段，但由于其分辨率高，成像效果好，未来在纳米尺度下研究微观剩余油形成机理、空间分布及微观剩余油与矿物颗粒相互作用领域会有很好的发展空间.
CT技术在构建数字岩心及微观剩余油研究得到了大量的成果，由于CT技术发展已经成熟和商业化，以及三维图像拥有相对较高的空间分辨率，使得CT技术在微观剩余油研究中依然具有广阔的前景.但由于分辨率和观察视野的矛盾性，要得到高分辨率的三维图像，必须减小观察视野，但是对于非均质性强的储层很难获取较小尺度下具有代表性的样品.因此，CT技术对于微观剩余油形成机理、赋存量以及空间分布规律研究、孔隙结构单一的岩心研究效果显著，而对于非均值储层的微观剩余油特性需要进一步研究，驱替过程中CT扫描快速稳定成像也是未来发展需要考虑的问题.
核磁共振成像技术对于运动的敏感性非常准确，即使在大型MRI设备上(空间分辨率小于1mm)，样品中流体移动1 μm距离，流速小于10 μm/s的条件下，流体的移动依然可以检测到，这种流体测量的准确性以及对非均质岩心(特别是全直径大尺度岩心)测量的适用性使得MRI在驱替过程的研究中有无可比拟的优势.未来MRI在微观剩余油特性研究需要基于该技术在提高分辨率方面的发展.由于信噪比的限制，核磁共振成像像素大小一般不会小于样品尺寸的百分之一.一些核磁共振成像仪器对于几毫米尺寸的岩心样品，可以获得10 μm分辨率的图像，然而这种空间分辨率能力依然低于目前X射线CT技术.较高分辨率的成像质量结合T2弛豫时间谱将会从多角度定量研究剩余油分布，一方面核磁共振成像技术可以直观获得相对较大尺度下驱替不同阶段剩余油分布图像，另一方面可以借助T2弛豫时间谱间接定量反映微观剩余油赋存量及赋存孔隙位置情况.研究多种条件下剩余油分布与孔隙空间的配置关系，更有助于应用于核磁共振测井，利用测井资料对储层含油特征进行评价.
3.3微观剩余油数值模拟前提是构建高质量的三维孔隙空间结构，比如X射线CT扫描法、核磁共振法构建三维孔隙结构质量的提高会使微观剩余油模拟效果更加真实.就数值模拟方法而言，基于网格化方法计算效率相对较慢，而且复杂孔隙结构多孔介质网格划分复杂，在微观剩余油领域研究较少.由于计算效率高，形态学方法和最大内切球方法在粗略要求计算剩余油空间分布的前提下是可行的，但单纯的数学形态等模拟很难捕捉真实流体流动复杂的变化，不能准确得到真实微观剩余油分布.无网格化中格子玻尔兹曼方法和孔隙网络方法对微观剩余油模拟研究相对成熟，由于格子玻尔兹曼是直接基于图像进行处理，得到的微观剩余油三维空间分布形态更加直观，孔隙网络方法则更倾向于微观剩余油形成机理研究，两者相互结合可以更综合的研究微观剩余油特性.对于格子玻尔兹曼，并行化快速计算、高密度比、高黏度比的两相流模拟是未来在微观剩余油分布领域发展的方向；孔隙网络方法更准确的描述孔隙结构、在非均质性强的多孔介质中的应用是未来研究的重点.
致谢 感谢审稿专家提出的修改意见和编辑部的大力支持！| [] | Al-Dhahli A, Geiger S, Van Dijke M I J .2014. Impact of pore-scale three-phase flow for arbitrary wettability on reservoir-scale oil recovery[J]. Journal of Petroleum Science and Engineering, 121 : 110–121. DOI:10.1016/j.petrol.2014.06.018 |
| [] | Al-Shalabi E W, Ghosh B .2016. Effect of pore-scale heterogeneity and capillary-viscous fingering on commingled waterflood oil recovery in stratified porous media[J]. Journal of Petroleum Engineering, 2016 : 1708929. |
| [] | Amiri H A A, Hamouda A A .2013. Evaluation of level set and phase field methods in modeling two phase flow with viscosity contrast through dual-permeability porous medium[J]. International Journal of Multiphase Flow, 52 : 22–34. DOI:10.1016/j.ijmultiphaseflow.2012.12.006 |
| [] | Bai Z Q, Wu S H, Fu Z G .2013. The distribution of microcosmic remaining oils after polymer flooding in Daqing oilfield[J]. Acta Petrolei Sinica (in Chinese), 34 (5) : 924–931. DOI:10.7623/syxb201305013 |
| [] | Baldwin B A, Yamanashi W S .1989. Detecting fluid movement and isolation in reservoir core with medical NMR imaging techniques[J]. SPE Reservoir Engineering, 4 (2) : 207–212. DOI:10.2118/14884-PA |
| [] | Boassen T, Kowalewski E, Hemmingsen P V. 2006. CRYO ESEM studies of emulsions and fluid distribution at pore scale[C].//International Symposium of the Society of Core Analysts. Trondheim, Norway. |
| [] | Bogdanov I, Kpahou J, Kamp A. 2011. Direct pore-scale modeling of two-phase flow through natural media[C].//Proceedings of the 2011 COMSOL Conference. Stuttgart, Germany. |
| [] | Bultreys T, Boone M A, Boone M N, et al .2015. Real-time visualization of Haines jumps in sandstone with laboratory-based microcomputed tomography[J]. Water Resources Research, 51 (10) : 8668–8676. DOI:10.1002/2015WR017502 |
| [] | Cao Y N .2015. Application CT scanning technology analysis micro-flooding experiments and the residual oil[J]. CT Theory and Applications (in Chinese), 24 (1) : 47–56. DOI:10.15953/j.1004-4140.2015.24.01.06 |
| [] | Chang C T P, Watson A T .1999. NMR imaging of flow velocity in porous media[J]. AIChE Journal, 45 (3) : 437–444. DOI:10.1002/(ISSN)1547-5905 |
| [] | Chardaire-Rivière C, Roussel J C .1991. Use of a high magnetic field to visualize and study fluids in porous media[J]. Magnetic Resonance Imaging, 9 (5) : 827–832. DOI:10.1016/0730-725X(91)90383-W |
| [] | Chatenever A, Calhoun Jr J C .1952. Visual examinations of fluid behavior in porous media-Part I[J]. Journal of Petroleum Technology, 4 (6) : 149–156. DOI:10.2118/135-G |
| [] | Chatzis I, Dullien F A L .1977. Modelling pore structure by 2-D and 3-D networks with application to sandstones[J]. Journal of Canadian Petroleum Technology, 16 (1) : 97–108. |
| [] | Chen P, Tao G, Dong M Z, et al .2013. The effects of the pore throat roughness on the water-oil flow in rock reservoirs[J]. Progress in Geophysics (in Chinese), 28 (2) : 824–829. DOI:10.6038/pg20130232 |
| [] | Chen Q, Kinzelbach W .2002. An NMR study of single-and two-phase flow in fault gouge filled fractures[J]. Journal of Hydrology, 259 (1-4) : 236–245. DOI:10.1016/S0022-1694(01)00599-6 |
| [] | Chen X P, Huo Q M, Lin J D, et al .2013. The relation between CBM content and the elastic parameters of CBM reservoirs:Reasoning and initial probing[J]. Geophys. (in Chinese), 56 (8) : 2837–2848. DOI:10.6038/cjg20130832 |
| [] | Cheng Y C, Di Q F, Gu C Y, et al .2015. Visualization study on fluid distribution and end effects in core flow experiments with low-field MRI method[J]. Journal of Hydrodynamics, 27 (2) : 187–194. DOI:10.1016/S1001-6058(15)60471-1 |
| [] | Coles M E, Hazlett R D, Spanne P, et al .1998. Pore level imaging of fluid transport using synchrotron X-ray microtomography[J]. Journal of Petroleum Science and Engineering, 19 (1-2) : 55–63. DOI:10.1016/S0920-4105(97)00035-1 |
| [] | Combes R, Robin M, Blavier G, et al .1998. Visualization of imbibition in porous media by environmental scanning electron microscopy:Application to reservoir rocks[J]. Journal of Petroleum Science and Engineering, 20 (3-4) : 133–139. DOI:10.1016/S0920-4105(98)00012-6 |
| [] | Davis Jr J A, Jones S C .1968. Displacement mechanisms of micellar solutions[J]. Journal of Petroleum Technology, 20 (12) : 1415–1428. DOI:10.2118/1847-2-PA |
| [] | De Beer F C, Middleton M F .2006. Neutron radiography imaging, porosity and permeability in porous rocks[J]. South African Journal of Geology, 109 (4) : 541–550. DOI:10.2113/gssajg.109.4.541 |
| [] | Deng R J .2002. Application of nuclear magnetic resonance imaging technology in water driving oil experiment[J]. Fault-Block Oil & Gas Field (in Chinese), 9 (4) : 33–37. DOI:10.3969/j.issn.1005-8907.2002.04.010 |
| [] | Deng S G, Lv W F, Liu Q J, et al .2014. Research on displacement mechanism in conglomerate using CT scanning method[J]. Petroleum Exploration and Development (in Chinese), 41 (3) : 330–335. DOI:10.11698/PED.2014.03.08 |
| [] | Dong H, Blunt M J .2009. Pore-network extraction from micro-computerized-tomography images[J]. Physical Review E, 80 (3 Pt 2) : 036307. |
| [] | Ersland G, Fernø M A, Graue A, et al .2010. Complementary imaging of oil recovery mechanisms in fractured reservoirs[J]. Chemical Engineering Journal, 158 (1) : 32–38. DOI:10.1016/j.cej.2008.11.049 |
| [] | Fatt I .1956. The network model of porous media[J]. Petroleum Transactions, AIME (207) : 160–181. |
| [] | Fenwick D H, Blunt M J .1998. Three-dimensional modeling of three phase imbibition and drainage[J]. Advances in Water Resources, 21 (2) : 121–143. DOI:10.1016/S0309-1708(96)00037-1 |
| [] | Ferréol B, Rothman D H .1995. Lattice-Boltzmann simulations of flow through Fontainebleau sandstone[J]. Transport in Porous Media, 20 (1-2) : 3–20. DOI:10.1007/BF00616923 |
| [] | Frisch U, Hasslacher B, Pomeau Y .1986. Lattice-gas automata for the Navier-Stokes equation[J]. Physical Review Letters, 56 (14) : 1505–1508. DOI:10.1103/PhysRevLett.56.1505 |
| [] | Gao J, Han D, Wang J L, et al .2009. Application of CT scanning image technique to study of oil saturation distribution in core displacement test[J]. Xinjiang Petroleum Geology (in Chinese), 30 (2) : 269–271. |
| [] | Ge L Z, Chen D X, Yang Q H .2014. Researching on remaining oil distribution in heterogeneous fluvial reservoirs by MRI[J]. China Offshore Oil and Gas (in Chinese), 26 (2) : 51–54, 60. |
| [] | Gong Y J, Liu S B, Liu K Y, et al .2015. Influence of reservoir wettability changes on oil-bearing features during tight oil accumulation:A case study of Jurassic tight oils in Sichuan Basin[J]. Petroleum Geology and Experiment (in Chinese), 37 (4) : 423–429. DOI:10.11781/sysydz201504423 |
| [] | Gunstensen A K, Rothman D H, Zaleski S, et al .1991. Lattice Boltzmann model of immiscible fluids[J]. Physical Review A, 43 (8) : 4320–4327. DOI:10.1103/PhysRevA.43.4320 |
| [] | He X Y, Chen S Y, Zhang R Y .1999. A lattice Boltzmann scheme for incompressible multiphase flow and its application in simulation of Rayleigh-Taylor instability[J]. Journal of Computational Physics, 152 (2) : 642–663. DOI:10.1006/jcph.1999.6257 |
| [] | Hilpert M, Miller C T .2001. Pore-morphology-based simulation of drainage in totally wetting porous media[J]. Advances in Water Resources, 24 (3-4) : 243–255. DOI:10.1016/S0309-1708(00)00056-7 |
| [] | Hirt C W, Nichols B D .1981. Volume of fluid (VOF) method for the dynamics of free boundaries[J]. Journal of Computational Physics, 39 (1) : 201–225. DOI:10.1016/0021-9991(81)90145-5 |
| [] | Hou J, Luo F Q, Li Z Q, et al .2014. The critical description scale study on core microscopic and reservoir macroscopic remaining oil[J]. Petroleum Geology and Recovery Efficiency (in Chinese), 21 (6) : 95–98. DOI:10.13673/j.cnki.cn37-1359/te.2014.06.024 |
| [] | Hou J, Qiu M X, Lu N, et al .2014. Characterization of residual oil microdistribution at pore scale using computerized tomography[J]. Acta Petrolei Sinica (in Chinese), 35 (2) : 319–325. DOI:10.7623/syxb201402012 |
| [] | Hou J, Zhang S K, Yuan S B .2006. Quantitative analysis of remaining oil in polymer flooding microscopic seepage images[J]. Journal of Hydrodynamics (in Chinese), 21 (1) : 41–45. |
| [] | Hou J, Zhang S K, Zhang Y H, et al .2011. Prediction of microscopic remaining oil distribution using fuzzy comprehensive evaluation[J]. Transport in Porous Media, 89 (3) : 533–545. DOI:10.1007/s11242-011-9784-z |
| [] | Hove A O, Ringen J K, Read P A .1987. Visualization of laboratory corefloods with the aid of computerized tomography of X-Rays[J]. SPE Reservoir Engineering, 2 (2) : 148–154. DOI:10.2118/13654-PA |
| [] | Hu G W, Li C F, Ye Y G, et al .2014. Observation of gas hydrate distribution in sediment pore space[J]. Chinese J. Geophys. (in Chinese), 57 (5) : 1675–1682. DOI:10.6038/cjg20140530 |
| [] | Hu X T, Li Y .2000. Study of microcosmic distribution of residual oil with stochastic simulation in networks[J]. Acta Petrolei Sinica (in Chinese), 21 (4) : 46–51. DOI:10.3321/j.issn:0253-2697.2000.04.009 |
| [] | Huang H B, Wang L, Lu X Y .2011. Evaluation of three lattice Boltzmann models for multiphase flows in porous media[J]. Computers & Mathematics with Applications, 61 (12) : 3606–3617. |
| [] | Huang S X, Zhang C M .2004. Effect of pore structural heterogeneity on residual oil distribution[J]. Journal of Jianghan Petroleum Institute (in Chinese), 26 (4) : 124–125. DOI:10.3969/j.issn.1000-9752.2004.04.046 |
| [] | Huang Y Z, Shang G H, Chen Y M .1995. A study on the mechanisms of cyclic water floodiing by nmri techniques[J]. Acta Petrolei Sinica (in Chinese), 16 (4) : 62–67. DOI:10.1007/BF02006258 |
| [] | Iglauer S, Fernø M A, Shearing P, et al .2012. Comparison of residual oil cluster size distribution, morphology and saturation in oil-wet and water-wet sandstone[J]. Journal of Colloid and Interface Science, 375 (1) : 187–192. DOI:10.1016/j.jcis.2012.02.025 |
| [] | Iglauer S, Paluszny A, Blunt M .2014. Corrigendum to "Simultaneous oil recovery and residual gas storage:A pore-level analysis using in situ X-ray micro-tomography"[Fuel 103 (2013) 905-914][J]. Fuel, 139 : 780. |
| [] | Jamaloei B Y, Asghari K, Kharrat R, et al .2010. Pore-scale two-phase filtration in imbibition process through porous media at high-and low-interfacial tension flow conditions[J]. Journal of Petroleum Science and Engineering, 72 (3-4) : 251–269. DOI:10.1016/j.petrol.2010.03.026 |
| [] | Jiang F, Tsuji T, Hu C H .2014. Elucidating the role of interfacial tension for hydrological properties of two-phase flow in natural sandstone by an improved lattice Boltzmann method[J]. Transport in Porous Media, 104 (1) : 205–229. DOI:10.1007/s11242-014-0329-0 |
| [] | Kareem R. 2015. Towards a Nanoscopic understanding of oil-sandstone wettability-implications for enhanced oil recovery[C]//18th European Symposium on Improved Oil Recovery. Dresden, Germany. |
| [] | Kong Q F, Zhou C C, Zhang Y, et al .2015. Numerical simulation methods of rock electrical properties based on digital cores:A review[J]. Progress in Geophysics (in Chinese), 30 (2) : 718–724. DOI:10.6038/pg20150232 |
| [] | Krummel A T, Datta S S, Münster S, et al .2013. Visualizing multiphase flow and trapped fluid configurations in a model three-dimensional porous medium[J]. AIChE Journal, 59 (3) : 1022–1029. DOI:10.1002/aic.v59.3 |
| [] | Kunz P, Zarikos I M, Karadimitriou N K, et al .2016. Study of Multi-phase flow in porous media:Comparison of SPH simulations with micro-model experiments[J]. Transport in Porous Media, 114 (2) : 581–600. DOI:10.1007/s11242-015-0599-1 |
| [] | Lang D J, Wang W M, Shang G H .2006. Visualization of combination corefloods and oil driving with aid of NMR imaging[J]. Computerized Tomography Theory and Applications (in Chinese), 15 (3) : 11–16. DOI:10.3969/j.issn.1004-4140.2006.03.003 |
| [] | Li A B, Ma L Y, Zhao Y X .1995. Wave velocity and quality factor Q of sandstone with porosity fluid[J]. Acta Geophysica Sinica (in Chinese), 38 (S1) : 288–294. |
| [] | Li H, Pan C, Miller C T .2005. Pore-scale investigation of viscous coupling effects for two-phase flow in porous media[J]. Physical Review E, 72 (2) : 026705. DOI:10.1103/PhysRevE.72.026705 |
| [] | Li J, Li H N .2003. The residual oil in Es2 reservoir, south area of Pucheng Oilfield[J]. Petroleum Exploration and Development (in Chinese), 30 (6) : 98–101. DOI:10.3321/j.issn:1000-0747.2003.06.030 |
| [] | Li X M, Cao J X, He X Y, et al .2004. Simulation of electric field response in the inhomogeneous medium using lattice Boltzmann method[J]. Chinese J. Geophys. (in Chinese), 47 (1) : 349–353. |
| [] | Liu X F, Sun J M, Wang H T .2009. Numerical simulation of rock electrical properties based on digital cores[J]. Applied Geophysics, 6 (1) : 1–7. DOI:10.1007/s11770-009-0001-6 |
| [] | Liu X F, Zhang W W, Sun J M .2013. Methods of constructing 3-D digital cores:A review[J]. Progress in Geophysics (in Chinese), 28 (6) : 3066–3072. DOI:10.6038/pg20130630 |
| [] | Liu X J, Zhu H L, Liang L X .2014. Digital rock physics of sandstone based on micro-CT technology[J]. Chinese J. Geophys. (in Chinese), 57 (4) : 1133–1140. DOI:10.6038/cjg20140411 |
| [] | Lu G F, Yang Y L, Feng Y, et al .2001. Study of remaining oil distribution in high aquifer stage by using a nuclear magnetic resonance imaging technique[J]. Journal of Jianghan Petroleum Institute (in Chinese), 23 (S1) : 98–99. DOI:10.3969/j.issn.1000-9752.2001.z1.038 |
| [] | Lu X G, Gao Z H, Zhao X J, et al .1996. Distribution regularities of remaining oil in heterogeneous reservoirs after polymer flooding[J]. Acta Petrolei Sinica (in Chinese), 17 (4) : 55–61. DOI:10.7623/syxb199604008 |
| [] | Lucy L B .1977. A numerical approach to the testing of the fission hypothesis[J]. Astronomical Journal, 82 : 1013–1024. DOI:10.1086/112164 |
| [] | LV W F, Leng Z P, Zhang Z B, et al .2013. Study on waterfllding mechanism in low permeability cores using CT scan ttechnology[J]. Petroleum Geology and Recovery Efficiency (in Chinese), 20 (2) : 87–90. DOI:10.3969/j.issn.1009-9603.2013.02.022 |
| [] | Ma B J .2016. The impact of water stoped on residual oil distribution in experimental study[J]. Science Technology and Engineering (in Chinese), 16 (12) : 254–259. |
| [] | Mandava S S, Watson A T, Edwards C M .1990. NMR imaging of saturation during immiscible displacements[J]. AIChE Journal, 36 (11) : 1680–1686. DOI:10.1002/(ISSN)1547-5905 |
| [] | Mattax C C, Kyte J R .1961. Ever see a water flood?[J]. Oil and Gas Journal, 59 : 115–128. |
| [] | McKellar M, Wardlaw N C .1982. A method of making two-dimensional glass micromodels of pore systems[J]. Journal of Canadian Petroleum Technology, 21 (4) . DOI:10.2118/82-04-03 |
| [] | Minto J M. 2014. Towards a rational design of gravel media water treatment filters:MRI investigation of the spatial heterogeneity in pollutant particle accumulation[Ph. D. thesis]. Glasgow:University of Glasgow. |
| [] | Mitchell J, Chandrasekera T C, Holland D J, et al .2013. Magnetic resonance imaging in laboratory petrophysical core analysis[J]. Physics Reports, 526 (3) : 165–225. DOI:10.1016/j.physrep.2013.01.003 |
| [] | Osher S, Sethian J A .1988. Fronts propagating with curvature-dependent speed:Algorithms based on Hamilton-Jacobi formulations[J]. Journal of Computational Physics, 79 (1) : 12–49. DOI:10.1016/0021-9991(88)90002-2 |
| [] | Pan C, Hilpert M, Miller C T .2004. Lattice-Boltzmann simulation of two-phase flow in porous media[J]. Water Resources Research, 40 (1) . DOI:10.1029/2003WR002120 |
| [] | Pan Y S, Tang J P, Li C Q .2008. NMRI test on two-phase transport of gas-water in coal seam[J]. Chinese J. Geophys. (in Chinese), 51 (5) : 1620–1626. DOI:10.3321/j.issn:0001-5733.2008.05.038 |
| [] | Pei Y Y, Tang G Y, Guo Z Y .2004. Progress in neutron radiography and their application[J]. Modern Instruments (in Chinese), 10 (5) : 17–22, 16. DOI:10.3969/j.issn.1672-7916.2004.05.004 |
| [] | Pereira G G, Pinczewski W V, Chan D Y C, et al .1997. Pore-scale network model for drainage-dominated three-phase flow in porous media[J]. Transport in Porous Media, 24 (2) : 167–201. |
| [] | Peyman M, Apostolos K. 2016. Modelling two-phase flow in multi-mineral porous media using coupled level set-volume of fluid (CLSVOF) method[C].//GeoConvention 2016:Optimizing Resources held in Calgary, Canada. Calgary, Canada:Calgary TELUS Convention Centre. |
| [] | Piller M, Casagrande D, Schena G, et al .2014. Pore-scale simulation of laminar flow through porous media[J]. Journal of Physics:Conference Series, 501 (1) : 012010. |
| [] | Prodanovi c ' M, Bryant S L .2006. A level set method for determining critical curvatures for drainage and imbibition[J]. Journal of Colloid and Interface Science, 304 (2) : 442–458. DOI:10.1016/j.jcis.2006.08.048 |
| [] | Prodanović M, Lindquist W B, Seright R S .2006. Porous structure and fluid partitioning in polyethylene cores from 3D X-ray microtomographic imaging[J]. Journal of Colloid and Interface Science, 298 (1) : 282–297. DOI:10.1016/j.jcis.2005.11.053 |
| [] | Raeini A Q, Blunt M J, Bijeljic B .2012. Modelling two-phase flow in porous media at the pore scale using the volume-of-fluid method[J]. Journal of Computational Physics, 231 (17) : 5653–5668. DOI:10.1016/j.jcp.2012.04.011 |
| [] | Ramstad T, Idowu N, Nardi C, et al .2012. Relative permeability calculations from two-phase flow simulations directly on digital images of porous rocks[J]. Transport in Porous Media, 94 (2) : 487–504. DOI:10.1007/s11242-011-9877-8 |
| [] | Ren Y M, Li M, Li Z C .2015. Microscopic residual oil study of water flooding based on pore-scale modeling[J]. Journal of Chongqing University of Science and Technology (Natural Sciences Edition) (in Chinese), 17 (1) : 1–5. DOI:10.3969/j.issn.1673-1980.2015.01.001 |
| [] | Romero-Zerón L B, Li L, Ongsurakul S, et al .2009. Visualization of waterflooding through unconsolidated porous media using magnetic resonance imaging[J]. Petroleum Science and Technology, 27 (17) : 1993–2009. DOI:10.1080/10916460802611267 |
| [] | Rowlinson J S .1979. Translation of J. D. van der Waals' "The thermodynamik theory of capillarity under the hypothesis of a continuous variation of density"[J]. Journal of Statistical Physics, 20 (2) : 197–200. DOI:10.1007/BF01011513 |
| [] | Ryazanov A V, Sorbi K S, Van Dijke M I J .2014. Structure of residual oil as a function of wettability using pore-network modelling[J]. Advances in Water Resources, 63 : 11–21. DOI:10.1016/j.advwatres.2013.09.012 |
| [] | Ryazanov A V, Van Dijke M I J, Sorbie K S .2009. Two-phase pore-network modelling:Existence of oil layers during water invasion[J]. Transport in Porous Media, 80 (1) : 79–99. DOI:10.1007/s11242-009-9345-x |
| [] | Shabro V, Prodanović M, Arns C H, et al. 2010. Pore scale modeling of two-phase flow[C].//18th International Conference on Computational Methods in Water Resources. Barcelona, Spain. |
| [] | Shakibaeinia A, Jin Y C. 2012. Mesh-free lagrangian modeling of mutiphase flow[C].//Proceedings of 20th Annual Conference of the CFD Society of Canada. Alberta, Canada. |
| [] | Shan X Q, Yan Q T .1999. The application summary of environment scanning electron microscope (ESEM) in oil[J]. Low Permeability Oil & Gas Field (in Chinese) (2) : 86–88, 68. |
| [] | Shan X W, Chen H D .1993. Lattice Boltzmann model for simulating flows with multiple phases and components[J]. Physical Review E, 47 (3) : 1815–1819. DOI:10.1103/PhysRevE.47.1815 |
| [] | Sheng Q, Shi X L, Liu W F, et al .2005. 3D imaging CT and multi-phase driving-substitute analytical system for rock-core and oil[J]. Computerized Tomography Theory and Applications (in Chinese), 14 (3) : 8–12. DOI:10.3969/j.issn.1004-4140.2005.03.002 |
| [] | Silin D, Patzek T .2006. Pore space morphology analysis using maximal inscribed spheres[J]. Physica A:Statistical Mechanics and Its Applications, 371 (2) : 336–360. DOI:10.1016/j.physa.2006.04.048 |
| [] | Silin D, Tomutsa L, Benson S M, et al .2011. Microtomography and pore-scale modeling of two-phase fluid distribution[J]. Transport in Porous Media, 86 (2) : 495–515. DOI:10.1007/s11242-010-9636-2 |
| [] | Song K P, Li S J, Fang W, et al .2005. Fluorescence analysis on changeable rules of microscopic remaining oil after polymer flooding[J]. Acta Petrolei Sinica (in Chinese), 26 (2) : 92–95. DOI:10.3321/j.issn:0253-2697.2005.02.021 |
| [] | Sukop M C, Huang H, Lin C L, et al .2008. Distribution of multiphase fluids in porous media:Comparison between lattice Boltzmann modeling and micro-x-ray tomography[J]. Physical Review E, 77 (2 Pt 2) : 026710. |
| [] | Sun G Z .2002. The preparation of thin section for oil-saturated rock samples[J]. Petroleum Geology & Experiment (in Chinese), 24 (6) : 566–567. DOI:10.3969/j.issn.1001-6112.2002.06.017 |
| [] | Sun Y Y, Song X M, Ma D S, et al .2014. Characteristics of remaining oil micro? distribution in Laojunmiao oilfield after waterflooding[J]. Xinjiang Petroleum Geology (in Chinese), 35 (3) : 311–314. |
| [] | Sutanto E, Davis H T, Scriven L E. 1990. Liquid distributions in porous rock examined by cryo scanning electron microscopy[C].//SPE 20518, SPE Annual Technical Conference and Exhibition. New Orleans, Louisiana:SPE, 459-474. |
| [] | Swift M R, Osborne W R, Yeomans J M .1995. Lattice Boltzmann simulation of nonideal fluids[J]. Physical Review Letters, 75 (5) : 830–833. DOI:10.1103/PhysRevLett.75.830 |
| [] | Tartakovsky A M, Ward A L, Meakin P .2007. Pore-scale simulations of drainage of heterogeneous and anisotropic porous media[J]. Physics of Fluids, 19 (10) : 103301. DOI:10.1063/1.2772529 |
| [] | Teng Y, Liu Y, Jiang L L, et al .2016. A visualization study on two-phase gravity drainage in porous media by using magnetic resonance imaging[J]. Magnetic Resonance Imaging, 34 (7) : 855–863. DOI:10.1016/j.mri.2016.03.004 |
| [] | Turner M L, Knüfing L, Arns C H, et al .2004. Three-dimensional imaging of multiphase flow in porous media[J]. Physica A:Statistical Mechanics and Its Applications, 339 (1-2) : 166–172. DOI:10.1016/j.physa.2004.03.059 |
| [] | Vizika O, Rosenberg E, Kalaydjian F .1998. Study of wettability and spreading impact in three-phase gas injection by cryo-scanning electron microscopy[J]. Journal of Petroleum Science and Engineering, 20 (3-4) : 189–202. DOI:10.1016/S0920-4105(98)00020-5 |
| [] | Wang D G, Hu Y L, Sun J J .2015. X-ray MCT based numerical analysis of residual oil pore-scale characteristics under various displacing systems[J]. Journal of Petroleum Science and Engineering, 135 : 168–176. DOI:10.1016/j.petrol.2015.09.015 |
| [] | Wang F F, Hu H G, He Z Q .2013. Remaining oil forming mechanism of occurrence research[J]. Guangdong Chemical Industry (in Chinese), 40 (3) : 66–67. DOI:10.3969/j.issn.1007-1865.2013.03.033 |
| [] | Wang K W, Guan J T, Fan Y H, et al .2005. The application of pore scale network model in seepage mechanics[J]. Advances in Mechanics (in Chinese), 35 (3) : 353–360. DOI:10.3321/j.issn:1000-0992.2005.03.006 |
| [] | Wang K W, Sun J M, Geng S C, et al .2006a. Percolation network study of shale effects on rock electrical properties under different salinity[J]. Chinese J. Geophys. (in Chinese), 49 (6) : 1867–1872. DOI:10.3321/j.issn:0001-5733.2006.06.036 |
| [] | Wang K W, Sun J M, Guan J T, et al .2006b. Network model modeling of microcosmic remaining oil distribution after polymer flooding[J]. Journal of China University of Petroleum (in Chinese), 30 (1) : 72–76. DOI:10.3321/j.issn:1000-5870.2006.01.015 |
| [] | Wang P, Cai J G, Xie Z H, et al .2002. Using oil-bearing thin section to study microscopic distribution feature of remaining oil[J]. Petroleum Geology and Recovery Efficiency (in Chinese), 9 (1) : 60–61. DOI:10.3969/j.issn.1009-9603.2002.01.019 |
| [] | Wang R F, Chen J B, Sun W .2002. Application of CT image technique to the research of oil displacement with water drive in Ultra-low permeability sandstone reservoir[J]. Progress in Geophysics (in Chinese), 23 (3) : 864–870. |
| [] | Wang S Y, Ayral S, Gryte C C .1984. Computer-assisted tomography for the observation of oil displacement in porous media[J]. Society of Petroleum Engineers Journal, 24 (1) : 53–55. DOI:10.2118/11758-PA |
| [] | Wang W M, Guo H K, Sun D Q, et al .1997. A study on polymer flooding by MRI techniques[J]. Acta Petrolei Sinica (in Chinese), 18 (4) : 54–60. DOI:10.1007/BF02951625 |
| [] | Xiang D, Huang D Z .2005. Research on microscopic distribution of water-drive remaining oil with C_5H_8O_2[J]. Chemical Engineering of Oil and Gas (in Chinese), 34 (4) : 293–295. DOI:10.3969/j.issn.1007-3426.2005.04.020 |
| [] | Xiao L Z .1995. A study on the application of NMR imaging to for[J]. Acta Petrolei Sinica (in Chinese), 16 (3) : 106–110. DOI:10.7623/syxb199503016 |
| [] | Yang Y F, Yao J, Van Dijke M I J, et al .2010. Effect of reservoir rock wettability on microcosmic distribution of residual oil after gas displacement[J]. Acta Petrolei Sinica (in Chinese), 31 (3) : 467–470. DOI:10.7623/syxb201003021 |
| [] | Yao J, Zhao X C, Yi Y J, et al .2005. The current situation and prospect on digital core technology[J]. Petroleum Geology and Recovery Efficiency (in Chinese), 12 (6) : 52–54. DOI:10.3969/j.issn.1009-9603.2005.06.017 |
| [] | Youssef S, Bauer D, Bekri S, et al. 2010. 3D in-situ fluid distribution imaging at the pore scale as a new tool for multiphase flow studies[C].//SPE Annual Technical Conference and Exhibition. Florence, Italy:SPE. |
| [] | Yu Y M, Hu Y Q, Liang W G, et al .2012. Study on pore characteristics of lean coal at different temperature by CT technology[J]. Chinese J. Geophys. (in Chinese), 55 (2) : 637–644. DOI:10.6038/j.issn.0001-5733.2012.02.027 |
| [] | Zhang M, Sun M X .2002. Applacation of X-ray CT technique to petroleum exploration[J]. Petroleum Geology & Oilfield Development in Daqing (in Chinese), 21 (2) : 15–16. DOI:10.3969/j.issn.1000-3754.2002.02.006 |
| [] | Zhang S K. 2007. The microscopic experiments and modeling for residual oil distribution after water flooding (in Chinese)[Ph. D. thesis]. Qingdao:China University of Petroleum. |
| [] | Zhang S K, Chen Y M, Hou J, et al .2006. Three dimensional visualization for CT tomographic image of microscopic flow law of rock pore[J]. Journal of Oil and Gas Technology (in Chinese), 28 (4) : 102–106. DOI:10.3969/j.issn.1000-9752.2006.04.029 |
| [] | Zhang X, Liu Y, Ma S .2009. Meshfree methods and their applications[J]. Advances in Mechanics (in Chinese), 39 (1) : 1–36. DOI:10.6052/1000-0992-2009-1-J2008-010 |
| [] | Zhang X H, Liu Q J, Lu X B. 2009. Comparisons of static, quasi-static and dynamic 3D porous media scale network models for two-phase immiscible flow in porous media[A].//Zhuang F G, Li J C eds. New Trends in Fluid Mechanics Research[M]. Berlin Heidelberg:Springer, 530. |
| [] | Zhao J L, Cao J T, Li C H, et al .2013. A review and perspective of logging evaluation technology of oil remaining in middle and high water cut stage[J]. Progress in Geophysics (in Chinese), 28 (2) : 838–845. DOI:10.6038/pg20130234 |
| [] | Zhao W J, Tan M J .2008. Review and prospect on NMR logging application in Shengli oilefield[J]. Progress in Geophysics (in Chinese), 23 (3) : 814–821. |
| [] | Zhao Y, Qu Z H, Liu Z .2002. Experimental study on water/oil displacement mechanisms in fractured reservoir by real sandstone micro-models[J]. Petroleum Exploration and Development (in Chinese), 29 (1) : 116–119. DOI:10.3321/j.issn:1000-0747.2002.01.031 |
| [] | Zhao Y C, Song Y C, Hao M, et al .2012. Application of MRI to CO2 displacement oil experiments[J]. Journal of Dalian University of Technology (in Chinese), 52 (1) : 23–28. |
| [] | Zhou G Z, Chen Z H, Ge W, et al .2010. SPH simulation of oil displacement in cavity-fracture structures[J]. Chemical Engineering Science, 65 (11) : 3363–3371. DOI:10.1016/j.ces.2010.02.029 |
| [] | 白振强, 吴胜和, 付志国.2013. 大庆油田聚合物驱后微观剩余油分布规律[J]. 石油学报, 34 (5) : 924–931. DOI:10.7623/syxb201305013 |
| [] | 曹永娜.2015. CT扫描技术在微观驱替实验及剩余油分析中的应用[J]. CT理论与应用研究, 24 (1) : 47–56. DOI:10.15953/j.1004-4140.2015.24.01.06 |
| [] | 陈萍, 陶果, 董明哲, 等.2013. 岩石孔喉道中表面粗糙度对油水两相流动的影响[J]. 地球物理学进展, 28 (2) : 824–829. DOI:10.6038/pg20130232 |
| [] | 陈信平, 霍全明, 林建东, 等.2013. 煤层气储层含气量与其弹性参数之间的关系--思考与初探[J]. 地球物理学报, 56 (8) : 2837–2848. DOI:10.6038/cjg20130832 |
| [] | 邓瑞健.2002. 核磁共振技术在水驱油实验中的应用[J]. 断块油气田, 9 (4) : 33–37. DOI:10.3969/j.issn.1005-8907.2002.04.010 |
| [] | 邓世冠, 吕伟峰, 刘庆杰, 等.2014. 利用CT技术研究砾岩驱油机理[J]. 石油勘探与开发, 41 (3) : 330–335. DOI:10.11698/PED.2014.03.08 |
| [] | 高建, 韩冬, 王家禄, 等.2009. 应用CT成像技术研究岩心水驱含油饱和度分布特征[J]. 新疆石油地质, 30 (2) : 269–271. |
| [] | 葛丽珍, 陈丹磬, 杨庆红.2014. 利用核磁共振成像技术研究河流相非均质储层剩余油分布[J]. 中国海上油气, 26 (2) : 51–54, 60. |
| [] | 公言杰, 柳少波, 刘可禹, 等.2015. 致密油充注过程中储层润湿性变化对含油性影响--以川中侏罗系致密油为例[J]. 石油实验地质, 37 (4) : 423–429. DOI:10.11781/sysydz201504423 |
| [] | 侯健, 罗福全, 李振泉, 等.2014. 岩心微观与油藏宏观剩余油临界描述尺度研究[J]. 油气地质与采收率, 21 (6) : 95–98. DOI:10.13673/j.cnki.cn37-1359/te.2014.06.024 |
| [] | 侯健, 邱茂鑫, 陆努, 等.2014. 采用CT技术研究岩心剩余油微观赋存状态[J]. 石油学报, 35 (2) : 319–325. DOI:10.7623/syxb201402012 |
| [] | 侯健, 张顺康, 袁士宝.2006. 聚合物驱微观渗流实验剩余油分布的定量分析研究[J]. 水动力学研究与进展, 21 (1) : 41–45. DOI:10.3969/j.issn.1000-4874.2006.01.007 |
| [] | 胡高伟, 李承峰, 业渝光, 等.2014. 沉积物孔隙空间天然气水合物微观分布观测[J]. 地球物理学报, 57 (5) : 1675–1682. DOI:10.6038/cjg20140530 |
| [] | 胡雪涛, 李允.2000. 随机网络模拟研究微观剩余油分布[J]. 石油学报, 21 (4) : 46–51. DOI:10.3321/j.issn:0253-2697.2000.04.009 |
| [] | 黄书先, 张超谟.2004. 孔隙结构非均质性对剩余油分布的影响[J]. 江汉石油学院学报, 26 (4) : 124–125. DOI:10.3969/j.issn.1000-9752.2004.04.046 |
| [] | 黄延章, 尚根华, 陈永敏.1995. 用核磁共振成像技术研究周期注水驱油机理[J]. 石油学报, 16 (4) : 62–67. DOI:10.1007/BF02006258 |
| [] | 孔强夫, 周灿灿, 张艳, 等.2015. 基于数字岩心岩石电性数值模拟方法综述[J]. 地球物理学进展, 30 (2) : 718–724. DOI:10.6038/pg20150232 |
| [] | 郎东江, 王为民, 尚根华.2006. 应用核磁共振成像(NMRI)技术研究复合剂驱油特征[J]. CT理论与应用研究, 15 (3) : 11–16. DOI:10.3969/j.issn.1004-4140.2006.03.003 |
| [] | 李爱兵, 马丽云, 赵永信.1995. 含孔隙液砂岩的波速和Q值[J]. 地球物理学报, 38 (S1) : 288–294. |
| [] | 李健, 李红南.2003. 濮城油田南区沙二段油藏剩余油研究[J]. 石油勘探与开发, 30 (6) : 98–101. DOI:10.3321/j.issn:1000-0747.2003.06.030 |
| [] | 李学民, 曹俊兴, 何晓燕, 等.2004. 用格子玻尔兹曼方法模拟非均匀介质中的电场响应[J]. 地球物理学报, 47 (2) : 349–353. DOI:10.3321/j.issn:0001-5733.2004.02.025 |
| [] | 刘学锋, 张伟伟, 孙建孟.2013. 三维数字岩心建模方法综述[J]. 地球物理学进展, 28 (6) : 3066–3072. DOI:10.6038/pg20130630 |
| [] | 刘向君, 朱洪林, 梁利喜.2014. 基于微CT技术的砂岩数字岩石物理实验[J]. 地球物理学报, 57 (4) : 1133–1140. DOI:10.6038/cjg20140411 |
| [] | 鲁国甫, 杨永利, 冯毅, 等.2001. 应用核磁共振成像技术研究特高含水期剩余油分布[J]. 江汉石油学院学报, 23 (S1) : 98–99. DOI:10.3969/j.issn.1000-9752.2001.z1.038 |
| [] | 卢祥国, 高振环, 赵小京, 等.1996. 聚合物驱油之后剩余油分布规律研究[J]. 石油学报, 17 (4) : 55–61. DOI:10.7623/syxb199604008 |
| [] | 吕伟峰, 冷振鹏, 张祖波, 等.2013. 应用CT扫描技术研究低渗透岩心水驱油机理[J]. 油气地质与采收率, 20 (2) : 87–90. DOI:10.3969/j.issn.1009-9603.2013.02.022 |
| [] | 马炳杰.2016. 注水停止对剩余油分布状态影响的实验研究[J]. 科学技术与工程, 16 (12) : 254–259. |
| [] | 潘一山, 唐巨鹏, 李成全.2008. 煤层中气水两相运移的NMRI试验研究[J]. 地球物理学报, 51 (5) : 1620–1626. DOI:10.3321/j.issn:0001-5733.2008.05.038 |
| [] | 裴宇阳, 唐国有, 郭之虞.2004. 中子照相技术及其应用[J]. 现代仪器, 10 (5) : 17–22, 16. DOI:10.3969/j.issn.1672-7916.2004.05.004 |
| [] | 任奕明, 李闽, 李中超.2015. 孔隙尺度水驱油模拟及残余油研究[J]. 重庆科技学院学报(自然科学版), 17 (1) : 1–5. DOI:10.3969/j.issn.1673-1980.2015.01.001 |
| [] | 单秀琴, 严启团.1999. 环境扫描电镜在石油技术领域中的应用综述[J]. 低渗透油气田 (2) : 86–88. |
| [] | 盛强, 施晓乐, 刘维甫, 等.2005. 岩心CT三维成像与多相驱替分析系统[J]. CT理论与应用研究, 14 (3) : 8–12. DOI:10.3969/j.issn.1004-4140.2005.03.002 |
| [] | 宋考平, 李世军, 方伟, 等.2005. 用荧光分析方法研究聚合物驱后微观剩余油变化[J]. 石油学报, 26 (2) : 92–95. DOI:10.3321/j.issn:0253-2697.2005.02.021 |
| [] | 孙国忠.2002. 饱含油岩心样品薄片的制作技术[J]. 石油实验地质, 24 (6) : 566–567. DOI:10.3969/j.issn.1001-6112.2002.06.017 |
| [] | 孙盈盈, 宋新民, 马德胜, 等.2014. 老君庙油田水驱后剩余油微观赋存状态[J]. 新疆石油地质, 35 (3) : 311–314. |
| [] | 王芳芳, 胡海光, 何志雄.2013. 剩余油形成机理与赋存状态研究综述[J]. 广东化工, 40 (3) : 66–67. DOI:10.3969/j.issn.1007-1865.2013.03.033 |
| [] | 王克文, 关继腾, 范业活, 等.2005. 孔隙网络模型在渗流力学研究中的应用[J]. 力学进展, 35 (3) : 353–360. DOI:10.3321/j.issn:1000-0992.2005.03.006 |
| [] | 王克文, 孙建孟, 耿生臣, 等.2006a. 不同矿化度下泥质对岩石电性影响的逾渗网络研究[J]. 地球物理学报, 49 (6) : 1867–1872. DOI:10.3321/j.issn:0001-5733.2006.06.036 |
| [] | 王克文, 孙建孟, 关继腾, 等.2006b. 聚合物驱后微观剩余油分布的网络模型模拟[J]. 中国石油大学学报(自然科学版), 30 (1) : 72–76. DOI:10.3321/j.issn:1000-5870.2006.01.015 |
| [] | 王朴, 蔡进功, 谢忠怀, 等.2002. 用含油薄片研究剩余油微观分布特征[J]. 油气地质与采收率, 9 (1) : 60–61. DOI:10.3969/j.issn.1009-9603.2002.01.019 |
| [] | 王瑞飞, 陈军斌, 孙卫.2008. 特低渗透砂岩储层水驱油CT成像技术研究[J]. 地球物理学进展, 23 (3) : 864–870. |
| [] | 王为民, 郭和坤, 孙佃庆, 等.1997. 用核磁共振成像技术研究聚合物驱油过程[J]. 石油学报, 18 (4) : 54–60. DOI:10.1007/BF02951625 |
| [] | 向丹, 黄大志.2005. 应用有机玻璃研究水驱剩余油微观分布[J]. 石油与天然气化工, 34 (4) : 293–295. DOI:10.3969/j.issn.1007-3426.2005.04.020 |
| [] | 肖立志.1995. 核磁共振成像在EOR研究中的应用[J]. 石油学报, 16 (3) : 106–110. DOI:10.7623/syxb199503016 |
| [] | 杨永飞, 姚军, VanDijke M I J, 等.2010. 油藏岩石润湿性对气驱剩余油微观分布的影响机制[J]. 石油学报, 31 (3) : 467–470. DOI:10.7623/syxb201003021 |
| [] | 姚军, 赵秀才, 衣艳静, 等.2005. 数字岩心技术现状及展望[J]. 油气地质与采收率, 12 (6) : 52–54. DOI:10.3969/j.issn.1009-9603.2005.06.017 |
| [] | 于艳梅, 胡耀青, 梁卫国, 等.2012. 应用CT技术研究瘦煤在不同温度下孔隙变化特征[J]. 地球物理学报, 55 (2) : 637–644. DOI:10.6038/j.issn.0001-5733.2012.02.027 |
| [] | 张敏, 孙明霞.2002. CT技术在油气勘探领域中的应用[J]. 大庆石油地质与开发, 21 (2) : 15–16. DOI:10.3969/j.issn.1000-3754.2002.02.006 |
| [] | 张顺康. 2007.水驱后剩余油分布微观实验与模拟[博士论文].青岛:中国石油大学. |
| [] | 张顺康, 陈月明, 侯健, 等.2006. 岩石孔隙中微观流动规律的CT层析图像三维可视化研究[J]. 石油天然气学报, 28 (4) : 102–106. DOI:10.3969/j.issn.1000-9752.2006.04.029 |
| [] | 张雄, 刘岩, 马上.2009. 无网格法的理论及应用[J]. 力学进展, 39 (1) : 1–36. DOI:10.6052/1000-0992-2009-1-J2008-010 |
| [] | 赵军龙, 曹军涛, 李传浩, 等.2013. 中高含水期剩余油测井评价技术综述与展望[J]. 地球物理学进展, 28 (2) : 838–845. DOI:10.6038/pg20130234 |
| [] | 赵文杰, 谭茂金.2008. 胜利油田核磁共振测井技术应用回顾与展望[J]. 地球物理学进展, 23 (3) : 814–821. |
| [] | 赵阳, 曲志浩, 刘震.2002. 裂缝水驱油机理的真实砂岩微观模型实验研究[J]. 石油勘探与开发, 29 (1) : 116–119. DOI:10.3321/j.issn:1000-0747.2002.01.031 |
| [] | 赵越超, 宋永臣, 郝敏, 等.2012. 磁共振成像在CO2驱油实验中应用[J]. 大连理工大学学报, 52 (1) : 23–28. |