-
石油是一种非可再生能源,其储量随着长期不断开采逐年减少,采收率逐年降低,而开采难度和开采成本却逐年增高。如何提高石油采收率、降低开采成本成为石油开采领域亟待解决的关键问题之一[1-2]。2019年,我国石油产量为1.91亿吨,而进口石油总量达5.06亿吨,石油对外依存度高达72.6%,远超国际公认的依存度警戒线(50%)。因此,如何高效可持续开采我国剩余油藏资源,减轻开采劳动强度,保持国产原油稳产、增产已成为我国石油工业迫切需要解决的问题,也是涉及国家战略,能源安全的重要问题。为了深度挖潜国内主力油田生产能力,提高低渗透、薄、差储层动用程度,大量新型提高采收率(EOR)技术如聚合物驱油[3]、泡沫驱油[4-5]、纳米驱油[6-8]、超声驱油[9-10]、微生物[11-15]驱油等先进驱油方法被推广应用[16]。然而,大部分新型驱油方法的微观驱替机理尚不明确,传统储层模拟方式通常存在无法直观观察、可重复性差、特征结构难以提取等缺陷,因此迫切需要一种可视化、可定制、可拓展的新型三维模拟储层岩心模型。
早期室内油藏驱替实验研究中所使用的岩心模型取自天然岩心,天然岩心从地层中钻取获得,具有与地层岩心相同的物理化学性质和构造情况,可以用于模拟某一区块地层的流动环境。然而天然岩心提取困难、开采成本高,难以满足实验研究中对岩心模型快速获取的需求,因此有研究者提出采用人造岩心代替天然岩心用于实验。通过胶结剂将矿物和岩石颗粒复合可制备得到人造岩心。该种岩心可方便快速地制造获取,但是胶结颗粒随机分布,难以对内部孔隙结构进行参数化设计,且不同岩心微观结构差异较大,实验重复性差。此外,由于岩石颗粒透光率低,因此不利于观察驱油过程中孔隙结构内的驱替特点和驱油现象。
为实现岩心模型参数化、可视化设计,学者们提出数字化定制岩心模型[17-18],通过玻璃刻蚀[19]、电子束光刻[20]、微立体光刻打印[21]等制造方法进行制造并应用于实验研究中,这些制造方法均具有快速可定制制造的优点,同时由于这些制造方法都可采用透光率较高的材料进行岩心模型制造,因此能够实现岩心模型内部流动情况的实时观测。在上述制造方法中,玻璃刻蚀和电子束光刻的工艺相对复杂,并且难以制造具有复杂三维孔喉结构的岩心模型,相比之下,面投影微立体光刻技术不仅具备高精度重复制造模型能力,同时也具有快速制造复杂三维结构的能力,因此在微尺度人造岩心制造方面具有重要意义。
近年来的驱油研究中,采用微立体光刻制造的岩心模型都仅具有二维结构特征[22-23],这些岩心模型相比于实际的驱油环境忽略了包括重力在内的许多复杂因素,不利于得出可靠的驱油研究结论,因此本文提出采用微立体光刻制造具有复杂孔喉结构的三维岩心模型。该岩心制造方式,实现了百微米级三维孔隙结构的可定制制造,区别于以往研究使用的二维孔隙模型,该方法制造的岩心具有可定制、高重复性、可视化的性能,且具有三维力场结构,在该种岩心中的驱油现象更符合实际驱油特点,因此本文基于该制造方法的优势和特点对岩心模型制造展开研究。
在本研究中,首先基于面投影微立体光刻技术构建人造岩心制造系统,并通过设计光路系统和调试设置打印参数使微立体光刻系统具有制造复杂三维微结构的能力。然后将岩石微孔隙结构抽象为微球堆叠模型,通过调整模型几何参数可以实现孔隙度定量控制。最后通过分析固化成型效果对岩心模型进行评估,并提出一种适用于逐层固化微立体光刻制造的复杂三维孔喉结构岩心模型。
-
岩心模型具有复杂的微纳结构,普通的微立体光刻打印设备无法满足打印精度要求,因此本研究中自主研制了如图1所示的上成型面投影微立体光刻系统用于制造岩心模型,该系统功能组件包括隔振平台、面投影紫外光机、光阑、凸透镜、反射镜、位移台、微立体光刻成型基底、树脂槽以及控制计算机。其中微立体光刻成型基底与位移台刚性连接,控制计算机可以通过控制位移台来调整微立体光刻成型基底的位置,同时控制计算机也可通过控制光机的光照时序进行逐层微立体光刻制造。该设备可打印具有良好透光性、流动性及较高打印精度的光敏树脂耗材,可固化成型百微米级别的特征结构尺寸,同时由于光敏树脂具有良好的透光性,因此微立体光刻成型的岩心模型内部结构也可被观察。
-
岩心模型具有复杂的空间微纳结构,普通面投影立体光刻无法满足岩心模型制造所需精度要求,因此该微立体光刻制造系统采用组合光路将面投影光机成像面尺寸聚焦,以此提高光照精度并获得高精度的固化成型效果。成型件的精度主要由DMD芯片的分辨率和光路的缩放倍数共同决定,本实验采用PRO 6500 型光学引擎,其DMD芯片TI DLP6500分辨率为1920 × 1080,其微镜阵列对角线为16.51 mm,即DMD芯片尺寸为16.51 mm×9.29 mm。光路中,光阑通过控制光束的通过数量进而调节光束的强弱,使投影到成型基底的紫外光图案不受杂光的影响以提高成型质量;凸透镜作为聚焦元件控制紫外光图案的缩放倍数以缩小单像素尺寸;反射镜可以改变光传播的途径,使成像平面变为水平面以便于液态树脂光固化成型。最终,成型幅面尺寸被调整为12.4 mm×7 mm,得到最终成型面单像素尺寸为6.48 μm。
-
模型的成型精度与立体光刻系统中的光照强度、曝光时间、每层成型厚度等参数密切相关。为提高逐层固化过程中三维结构的制造精度,需要对这些参数进行精密调试。其中光照强度太小或曝光时间太短会导致固化层之间无法可靠粘结,甚至导致光敏树脂无法完全固化而变成在树脂液中悬浮的絮状固化物,而光照强度太大或曝光时间太长会导致固化层厚度过大、成型精度低。考虑到岩心模型的特征尺寸较小,故应选用较小的固化层厚来获得较高的打印精度。通过如下公式并结合固化成型效果可对光照强度、曝光时间等参数进行调试设置。
$${C_{\rm{d}}} = {D_{\rm{p}}}\ln \left(\frac{{{E_{\max }}}}{{{E_{\rm{c}}}}}\right)$$ (1) $${E_{\max }} = ISt$$ (2) 其中,Dp为光敏树脂透射深度,mm;Emax为成型面所接收到的能量,J;Ec为光敏树脂临界固化所需要光能,J;Cd为固化深度,mm;
$ I $ 为光功率密度,W·mm−2;$ S $ 为成型面积,mm2;$ t $ 为曝光时间,s。将光照强度、曝光时间、每层成型厚度等参数调试至最佳后该微立体光刻设备才可制造具有较高精度的三维岩心模型。通过多次固化成型并测量统计,得出单像素成型尺寸为6.55 μm,成型尺寸误差小于0.15 μm,相对误差小于2%,满足岩心模型制造要求。
-
岩心模型的微立体光刻制造采用逐层固化的方式进行。首先,在计算机辅助设计软件中建立模拟岩心三维结构模型,对模型进行层切片得到所需的每个打印幅面图形。其次,将切片图像导入控制计算机中,由程序控制光机的紫外光光照时序、投影图像以及微立体光刻成型基底的移动。然后,在控制程序的执行下液态微立体光刻树脂将会逐层固化并与前一固化层相互粘结堆积,最终形成三维岩心模型实体。最后,将成型的三维模型用丙酮溶液浸泡溶去表面的残留树脂,去离子水清洗、干燥后则可得到如图2所示的复杂三维模型。采用微立体光刻制造的多孔介质岩心模型具有良好的透光性,在模拟流动实验中可易于观察岩心内的流动情况,同时该方法工艺简单、模型成型速度快,因此可满足石油研究中岩心模型可定制、快速制造、可视化的需求。
-
地层岩心中的油藏存储空间由岩石颗粒的孔喉和缝隙构成[24-25],本文通过对如图3所示的岩心SEM图像进行分析,提取出岩石颗粒的堆积结构作为岩心模型的主要特征结构。由于实际岩心中的岩石颗粒都具有复杂且不规则的几何形状,难以有效地提取结构的特征参数,更难以研究这些结构参数对原油驱替效率的影响,因此本文将不规则的岩石颗粒简化为半径一致的球形颗粒,得到如图4所示的颗粒堆积结构,该结构作为主要的三维结构被用于设计、改进和制造岩心模型。
-
考虑到微球堆叠岩心模型具有复杂三维孔隙结构,能够较好地模拟地层岩石颗粒间的流动环境,因此本研究采用微球堆叠结构设计岩心模型。微球堆叠岩心模型具有如图5所示的包括简单立方堆积、体心立方堆积、面心立方堆积、六方紧密堆积等多种不同的堆积方式,不同的堆积方式也将导致设计的多孔介质岩心模型具有不同的孔隙结构。由于本研究采用逐层固化的方式制造岩心模型,通过后文分析可知简单立方堆积具有更好的成型效果,因此本研究提出采用如图6所示的简单立方堆积岩心模型进行微立体光刻制造。
-
为了实现岩心模型孔隙度可定制化,需要对岩心模型孔隙度几何参数进行分析。如图7所示立方体为简单立方堆积的微球堆叠模型基本单元,通过调整球心距s和球半径r能够定量改变岩心模型的孔隙度δ,根据以上两个结构参数计算岩心模型基本单元内的孔隙体积以及基本单元体积,计算两者之比即为该结构的孔隙度,求得孔隙度计算公式如下:
$$\delta = \frac{{8{r^3}}}{{3{s^3}}}{\text{π}} - \frac{{6{r^2}}}{{{s^2}}}{\text{π}} + \frac{\text{π}}{4} + 1$$ (3) 然而该模型的孔隙度可设计范围受到模型的结构可制造性限制,通过对球颗粒的接触关系和孔隙的联通性两方面分析考虑,可计算得出该模型的孔隙度设计范围。
考虑球颗粒的接触关系,当相邻球颗粒之间无接触时,每个球颗粒均为悬空实体而无法统一为一个刚性整体,因此该模型无法被制造。分析可知当球半径与球心距之比等于1/2时,相邻球颗粒恰好相切使得所有球颗粒相互刚性连接为一整体,此时模型具有最大的孔隙度,计算可知该模型的最大可设计孔隙度δmax为47.6%。
考虑孔隙的连通性,当处于同一平面的相邻四个球颗粒互相接触时,球颗粒间的孔隙将会被阻断,岩心模型也将由于不存在连续的孔喉结构而无法被用于流动实验。分析可知当球半径与球心距之比小于
$ \surd 2 $ /2时,模型的孔隙才可相互连通,计算可知该模型的最小可设计孔隙度δmin为3.5%。通过对简单立方堆积结构的分析,我们可以了解到简单立方堆积的微球堆叠模型可以被用于设计孔隙度处于3.5%~47.6%之间的具有三维复杂孔喉结构的岩心模型。在目前的岩石孔隙度研究中,岩浆岩孔隙度一般为0.3%~15%,沉积岩孔隙度一般为2.5%~31.7%,而常用于驱油实验研究的砂岩孔隙度一般为10%~30%,因此该模型可满足实验研究中岩心孔隙度的设计要求。
通过孔隙度分析可知球颗粒直径应为球心距的1~
$ \surd 2 $ 倍,因此本研究还尝试采用随机数生成三维直径序列用于建立颗粒堆积的非均质岩心模型。我们设计球心距为0.45 mm,考虑到微立体光刻的制造精度,颗粒直径在将在0.45~0.55 mm范围内随机生成,并以此建立5×5×5微球堆积岩心模型,如图8所示为该模型的建模过程以及制造效果,模型具有非均质特性,在三维空间内具有不规则的孔隙尺寸分布。一方面控制随机数的生成并建模可以对岩心孔隙度范围进行设计,另一方面对相同非均质岩心模型微立体光刻制造具有重复制造性。既可满足对不同油层流动环境参数的模拟需求,也可满足相同流动环境重复性试验研究的需求。 -
本文基于地层岩心中岩石颗粒的堆积结构,简化设计并通过微立体光刻技术将其制造,为了研究微立体光刻制造三维结构的可靠性,通过显微镜观察如图9所示的岩心模型实际成型情况,可发现微立体光刻岩心模型实体与建立的理想岩心模型存在差异,主要体现在岩心实体模型的表面存在与理想岩心模型光滑表面不同的阶梯状结构,通过分析微立体光刻的制造过程,可知产生这样的差异是由于采用了逐层固化成型的方式。如图10所示,该方式制造的球体模型实际由多个薄圆柱组合而成,分析实体的逐层堆叠过程可得出固化层厚d是影响成型表面阶梯结构尺寸的主要参数,随着层厚的减小,阶梯的高度也会逐渐减小,微立体光刻成形的岩心模型表面会更相似于设计的理想岩心模型表面。然而,无论层厚如何减小,阶梯结构也无法被完全消除,因此通过逐层固化方式无法制造光滑曲面,所以也无法完全一致地将设计的理想岩心模型制造出来。然而即使逐层固化产生会阶梯表面,也破坏了理想岩心模型与实体模型的一致性,但正因为制造出现了阶梯结构,使得模型表面凹凸不平,所以微立体光刻岩心模型表面实际更接近于真实岩石颗粒的不规则表面,因此在设计固化层厚时,可根据被模拟岩心的表面形貌需求进行设计,而不仅仅是取尽可能小的固化层厚。
另一方面,本文通过分析打印过程中每一层的有效固化堆积,可知如图11所示的堆积结构存在无法有效固化成型的薄圆柱,这是因为固化该层时,薄圆柱无法与已固化的实体接触粘结,所以固化后的薄圆柱会漂浮在树脂液中,并可能导致岩心模型孔喉结构堵塞。考虑到体心立方堆积、六方紧密堆积、面心立方堆积都存在有类似的结构,也存在无法有效固化粘结的缺陷层,所以最终选择采用简单立方堆积方式设计微球堆叠岩心模型,该堆积结构可解决逐层固化中可能产生缺陷层的问题,可获得较好的岩心模型成型效果。
鉴于逐层固化微立体光刻会导致模型与建模尺寸和表面形貌产生差异,因此我们通过显微镜对制造的岩心进行测量,以此标定该方法的制造精度,如图12所示为测量尺寸。我们随机测量两种成型平面上的微球直径各25个,并计算出平均值来衡量不同平面上微球的成型精度。结果显示,在x−y投影平面内尺寸误差平均值为2.36%,在x−z和y−z平面内尺寸误差平均值为0.13%。
-
本研究中使用的树脂为esun的通用刚性树脂,该树脂属于聚氨酯丙烯酸酯,在其中加入0.2%苏丹Ⅰ作为光吸收剂提高微粒体光刻精度后,可微立体光刻直径为0.5 mm、球心距为0.45 mm的微球堆叠模型球颗粒,且制造成型后采用显微镜观察可确认岩心模型孔隙相互连通。本研究中我们还使用其他树脂材料作为微立体光刻耗材,将尺寸相同的微球堆叠岩心模型重复制造,岩心模型成型效果如图13所示,其中图13(a)所示岩心模型采用zDental Model沙黄树脂成型,该树脂也属于丙烯酸树脂聚氨酯体系,可微立体成型孔隙相互连通的微球堆积岩心模型,微球成型尺寸与模型尺寸误差小,但透光性不足,无法观察岩心模型空隙内流动;图13(b)所示岩心模型采用由光聚合单体HDDA、PEGDA、TMPTA、PPTTA以及光引发剂819和光吸收剂171配置形成的光固化树脂体系固化成型,微立体成型的三维孔隙相互连通,并可成型较小尺寸孔喉,成型速度快但树脂收缩率大易产生裂纹,导致结构破坏。
-
本文基于石油开采研究需求,提出采用微立体光刻的方式进行岩心模型制造,并提出了能较好模拟地层岩心结构的微球堆叠岩心模型,同时对该模型的孔隙度特征以及微立体光刻制造效果进行研究,得出以下结论。
(1)本文的微立体光刻设备具有6.55 μm的像素点成型尺寸精度,可制造百微米级别的三维复杂孔喉结构,且基于微立体光刻制造工艺和材料的特性,该方法不仅能快速地设计制造岩心模型,而且制造的岩心模型具有较高的透光率可满足观察岩心内部流动情况的需求。
(2)微球堆叠岩心模型在结构上具有复杂的三维孔喉结构,同时立体光刻的逐层固化制造也增加了微球表面结构的复杂性,具有与地层岩性相似的宏观孔喉结构与微观表面特征。此外多种微球堆积方式中仅有简单立方堆积岩心模型可适应微立体光刻的逐层固化制造方式,避免岩心缺陷的出现。
(3)本实验微立体光刻制造的五百微米微球实体,在x−y平面上成型的直径平均误差为2.36%,在x−z平面或y−z平面上成型的直径平均误差为0.13%。
(4)简单立方堆积的微球堆叠岩心模型孔隙度由该模型的球心距与球半径之比定量控制,由模型的制造和应用分析得出该岩心模型可设计的孔隙度范围为3.5%~47.6%,该设计范围足以满足大多数岩心的孔隙度设计需求。此外,根据孔隙度设计的结论可设计非等径微球的非均质微球堆积模型,并同样可采用微立体光刻成型。
Fabrication of a three-dimensional simulated reservoir core model based on area projection micro-stereolithography
-
摘要: 首先搭建具有高精度面投影微立体光刻设备,通过理论分析和实验相结合的方法获得最优打印工艺参数,然后提出一种可用于模拟地层岩心的微球堆叠岩心模型,并通过分析岩心模型成型机理,选取具有更高成型精度的堆积方式对岩心模型进行设计。该模拟岩心制造方法具有对特殊岩心结构制造的高适应性,为实验室显微镜下研究多种EOR技术微观驱替机理提供了新思路。Abstract: Petroleum exploitation plays a very important role in national energy security. With continuous exploitation of oil fields in my country, the efficiency of conventional water injection oil production is decreasing year by year. Enhanced oil recovery (EOR) technologies, such as polymer flooding, microbial flooding, micro–nano flooding, and other flooding technologies have been proposed and developed for application. However, the microscopic displacement mechanism and displacement effect of these technologies are still unclear. Current oil displacement research needs to be verified by core displacement experiments. However, the current displacement experiments all use artificial cores, glass etching channels, photoetched microchannels, etc., as the oil displacement environment. These displacement environments are insufficient in terms of oil displacement dimensions and observation phenomena. Due to this, there is an urgent need for a core manufacturing method that is more suitable for laboratory oil displacement research. In this study, we proposed a method for manufacturing a simulated three-dimensional core structure based on micro-stereolithography technology. This method not only has the advantages of fast manufacturing speed and high forming accuracy, but is also able to realize the visualization, parameterization, and customized design of a micron structure. The core model self-searched by stereo lithography has a three-dimensional pore structure in the order of hundreds of microns and can be used to simulate the experimental study of reservoir displacement flow mechanism. In this research, a high-precision surface projection micro-stereolithography equipment was built, and the optimal printing process parameters were obtained through a combination of theoretical analysis and experiments. Then, a microsphere stacked core model was proposed that can be used to simulate formation cores. By analyzing the forming mechanism of the core model, a stacking method was selected with a higher forming accuracy to design the core model. Finally, the core of a 100-micron-sized microsphere accumulation was realized by micro-stereolithography to achieve three-dimensional molding. The simulated core manufacturing method in this study has high adaptability to special core structure manufacturing and provides a new idea for studying the microscopic displacement mechanism of various EOR technologies under a laboratory microscope.
-
[1] Dai C L, Fang J C, Jiao B L, et al. Development of the research on EOR for carbonate fractured-vuggy reservoirs in China. 中国石油大学学报, 2018, 42(6): 67 戴彩丽, 方吉超, 焦保雷, 等. 中国碳酸盐岩缝洞型油藏提高采收率研究进展. 中国石油大学学报(自然科学版), 2018, 42(6):67) [2] 朱维耀, 岳明, 刘昀枫, 等. 中国致密油藏开发理论研究进展. 工程科学学报, 2019, 41(9):1103 Zhu W Y, Yue M, Liu Y F, et al. Research progress on tight oil exploration in China. Chin J Eng, 2019, 41(9): 1103 [3] Samala R, Chaudhuri A, Vishnudas R, et al. Numerical analysis of viscous fingering and oil recovery by surfactant and polymer flooding in five-spot setup for water and oil-wet reservoirs. Geomech Geophys Geo-Energy Geo-Resour, 2020, 6(1): 3 doi: 10.1007/s40948-019-00124-1 [4] Ma K, Liontas R, Conn C A, et al. Visualization of improved sweep with foam in heterogeneous porous media using microfluidics. Soft Matter, 2012, 8(41): 10669 doi: 10.1039/c2sm25833a [5] Lv Q C, Li Z M, Li B F, et al. Study of nanoparticle-surfactant-stabilized foam as a fracturing fluid. Ind Eng Chem Res, 2015, 54(38): 9468 doi: 10.1021/acs.iecr.5b02197 [6] Olayiwola S O, Dejam M. A comprehensive review on interaction of nanoparticles with low salinity water and surfactant for enhanced oil recovery in sandstone and carbonate reservoirs. Fuel, 2019, 241: 1045 doi: 10.1016/j.fuel.2018.12.122 [7] Hendraningrat L, Li S D, Torsæter O. A coreflood investigation of nanofluid enhanced oil recovery. J Petrol Sci Eng, 2013, 111: 128 doi: 10.1016/j.petrol.2013.07.003 [8] 光新军, 豆宁辉, 贾云鹏, 等. 纳米技术在石油工程中的应用前景. 钻采工艺, 2019, 42(3):34 doi: 10.3969/J.ISSN.1006-768X.2019.03.10 Guang X J, Dou N H, Jia Y P, et al. Application prospects of nanotechnology in petroleum engineering. Drill Prod Technol, 2019, 42(3): 34 doi: 10.3969/J.ISSN.1006-768X.2019.03.10 [9] Hamida T, Babadagli T. Displacement of oil by different interfacial tension fluids under ultrasonic waves. Colloids Surf A, 2008, 316(1-3): 176 doi: 10.1016/j.colsurfa.2007.09.012 [10] 葛丹. 三元复合驱油泥稳定性及超声—破乳洗油的研究[学位论文]. 大庆: 东北石油大学, 2018. Ge D. Study on Stability of ASP Flooding Sludge and Ultrasonic-Demulsification Oil Removal [Dissertation]. Daqing: Northeast Petroleum University, 2018. [11] 刘卫. 浅谈石油开采中微生物采油技术的应用. 化工设计通讯, 2017, 43(7):63 doi: 10.3969/j.issn.1003-6490.2017.07.060 Liu W. Discussion on application of microbial oil recovery technology in oil exploitation. Chem Eng Des Commun, 2017, 43(7): 63 doi: 10.3969/j.issn.1003-6490.2017.07.060 [12] 刘杰, 王超, 孙朝阳. 微生物驱油在油田的应用. 化工设计通讯, 2019, 45(1):33 doi: 10.3969/j.issn.1003-6490.2019.01.030 Liu J, Wang C, Sun C Y. Application of microbial flooding in oil fields. Chem Eng Des Commun, 2019, 45(1): 33 doi: 10.3969/j.issn.1003-6490.2019.01.030 [13] de Araujo L L G C, Sodré L G P, Brasil L R, et al. Microbial enhanced oil recovery using a biosurfactant produced by Bacillus safensis isolated from mangrove microbiota - Part I biosurfactant characterization and oil displacement test. J Petrol Sci Eng, 2019, 180: 950 doi: 10.1016/j.petrol.2019.06.031 [14] Nazina T, Sokolova D, Grouzdev D, et al. The potential application of microorganisms for sustainable petroleum recovery from heavy oil reservoirs. Sustainability, 2020, 12(1): 15 [15] Haq B, Liu J S, Liu K Y, et al. The role of biodegradable surfactant in microbial enhanced oil recovery. J Petrol Sci Eng, 2020, 189: 106688 doi: 10.1016/j.petrol.2019.106688 [16] 董喜庆, 马晓燕. 三次采油工艺技术措施. 化工设计通讯, 2017, 43(7):42 doi: 10.3969/j.issn.1003-6490.2017.07.039 Dong X Q, Ma X Y. Technical measures for three-stage oil recovery. Chem Eng Des Commun, 2017, 43(7): 42 doi: 10.3969/j.issn.1003-6490.2017.07.039 [17] 许丽. 砂岩油藏孔隙网络模型构造及应用研究[学位论文]. 成都: 西南石油大学, 2015. Xu L. Pore Network Model Construction of Sandstone Reservoir and Application [Dissertation]. Chengdu: Southwest Petroleum University, 2015. [18] 孙泽. 基于三维CT扫描图像的致密砂岩孔隙数据提取方法. 华北科技学院学报, 2020, 17(1):6 doi: 10.3969/j.issn.1672-7169.2020.01.002 Sun Z. A method of extracting pore date of tight sandstone based on 3D CT scanning image. J North China Inst Sci Technol, 2020, 17(1): 6 doi: 10.3969/j.issn.1672-7169.2020.01.002 [19] Abgrall P, Gue A M. Lab-on-chip technologies: making a microfluidic network and coupling it into a complete microsystem—a review. J Micromech Microeng, 2007, 17(5): R15 doi: 10.1088/0960-1317/17/5/R01 [20] Xu K, Zhu P X, Huh C, et al. Microfluidic investigation of nanoparticles' role in mobilizing trapped oil droplets in porous media. Langmuir, 2015, 31(51): 13673 doi: 10.1021/acs.langmuir.5b03733 [21] Beauchamp M J, Nordin G P, Woolley A T. Moving from millifluidic to truly microfluidic sub-100-μm cross-section 3D printed devices. Anal Bioanal Chem, 2017, 409: 4311 doi: 10.1007/s00216-017-0398-3 [22] Vavra E D, Zeng Y C, Xiao S Y, et al. Microfluidic devices for characterizing pore-scale event processes in porous media for oil recovery applications. J Vis Exp, 2018, 131: e56592 [23] Li H X, Zhang T J. Imaging and characterizing fluid invasion in micro-3D printed porous devices with variable surface wettability. Soft Matter, 2019, 15(35): 6978 doi: 10.1039/C9SM01182J [24] 黄德利. 大庆油田特低渗透裂缝性油藏渗流特征研究及应用[学位论文]. 大庆: 大庆石油学院, 2010. Huang D L. Study of the Fluid Flow in Ultra-low Permeability Reservoir with Fractures and its Application at Daqing Olifield [Dissertation]. Daqing: Daqing Petroleum Institute, 2010. [25] Ding X Y. Pore Structure Characteristics of F Oil Reservoir in Daqing and the Impact on the Seepage [Dissertation]. Daqing: Northeast Petroleum University, 2013. -