油气二次运移数值模拟分析.pdf

第 2 9卷第 4期 -5 天然气肚 O I L＆ G A S G E O L O G Y 2 0 0 8年 8月 文章编号 0 2 5 3 9 9 8 5 2 0 0 8 0 4 0 5 2 7 0 6 油气二次运移数值模拟分析 周 波 ， 金之钧 ， 王 毅 1 ．中国石油化工股份有 限公司 石油勘探开发研 究院 ， 北京 1 0 0 0 8 3 ； 2 ．中国地质大学 ， 北京 1 0 0 0 8 3 摘要 通过分析物理实验和 数学模拟认为 ， 油气运移 的模式是 多样化 的。实 际地质 条件 的油气运移模式有 3种 稳定 式、 指进 式和优势式。不同运移模式的运移规律 并不相 同， 用任何 一种单 一的油气运移数值模拟方 法来描述 所有的油气运移 模 式是不科学的。不同的运移模式对应着不 同的数值模 拟方法 ， 其 中逾 渗方 法适 用于描述 指进 式运移和优 势 式运 移过 程， 达西渗流方程适用于描述 稳定 式运移模式。流体 势的方法 同时可 以用 于分 析 3种模 式下的运移 ， 但其模拟 结果只是 一 种定性结果。此外， 数学模拟必须以详实的地质资料为基础， 结合油气显示和地化资料， 综合分析油气运移过程， 才能 确定真正 的油气运移过程。 关键词 运移路径 ； 运移模 式； 数学模拟 ； 油气二次运 移 中图分类号 T E 1 2 2 ． 1 文献标识码 A Nu m e r i c a l s i mu l a t i o n a n a l y s i s o f s e c o n da r y h y dr o c a r b o n mi g r a t i o n Z h o u B o ， J i n Z h i j u n ， Wa n g Y i 1 ． S I N O P E C E x p l o r a t i o n a n d p r o d u c t i o n R e s e a r c h I n s t i t u t e ， B e ij i n g 1 0 0 0 8 3 ， C h i n a； 2 ． C h i n a U n i v e r s i t y o f G e o s c i e n c e s ， B e ij i n g 1 0 0 0 8 3 ， C h i na Abs t r a c t T he a n a l y s i s o f t h e p h y s i c a l e x p e rime n t a n d n u me ric a l s i mu l a t i o n r e s u l t s s h o ws v a rio u s o i l mi g r a t i o n mo d e s ． Th e r e a r e t h r e e mo d e s o f o i l mi gra t i o n u nd e r a c t u a l g e o l o g i c a l c o n d i t i o n ss t a b l e mo d e， fin g e rin g mo d e， a n d s p ring mo d e ． Di f f e r e n t mi gra t i o n mo d e s h a v e d i f f e r e n t p a t t e r ns o f mi gra t i o n， t h u s u s i n g s i n g l e n u me ric a l s i m u l a t i o n me t ho d t o s i mu l a t e all t h e mo d e s i s n o t v e r y s c i e n t i fic ． Th e r e f o r e， d i f f e r e n t mi g r a t i o n mo d e s s h all b e d e s c rib e d wi t h d i ffe r e n t n u me ric al s i mu l a t i o n me t h o ds ． I n v a s i o n p e r c o l a t i o n me t h o d i s s u i t a b l e f o r t h e de s c rip t i o n o f t h e fin g e rin g mo d e a n d t h e s p rin g mo d e， wh i l e t h e Da r c y e q u a t i o n i s a pp l i c a b l e for t h e s t a b l e mo d e ． Th e flu i d p o t e n t i al me t h o d c a n b e a p p l i e d t o all t h e t h r e e mo d e s ， b u t i t s r e s u l t s a r e q u a l i t a t i v e ． I n o r d e r t o d e t e r mi ne t h e r e al o i l m i gra t i o n p r o c e s s ， t h e m a t h e m a t i c a l s i m u l a t i o n s h al l b e b a s e d o n d e t a i l e d g e o l o gi c a l d a t a ， o i l a n d g a s s h o w s a n d g e o c h e mi c al d a t a ． Ke y wo r d s mi gra t i o n p a thwa y；mi gra t i o n mo d e；ma t h e ma t i c al s i mu l a t i o n；s e c o n d a ry h y d r o c a r b o n mi gra t i o n 随着计算机技术 的发展 ， 数学模拟技术成 了 解决物理模拟局 限性 的有利工具 ] 。早先 ， 孔 隙介质内多相流动的理论被应用于油气运移的分 析中 一 J 。在这种情况下 ， 油气的运移过程被视为 在宏观上均一 的多相流 动过程 J ， 因而可 以用 多相渗流的扩散方程和多相条件下的达西定律对 油、 气 、 水各相流体的流动过程加 以描述。为了描 述储层非均质性 的作用 ， 前人将渗透率大小代入 达西渗流方程 ， 通过数学方法解析该方程 ， 进而使 得该方法考虑了岩性对运移 的控制作用 J 。此 外 ， 流体势等概念也被用于模拟油气运移口 m ] 。 E n g l a n d等 在 Hu b b e 流体势概念的基础上考 虑了毛细管力 ， 新 的流体势方 程中考虑 了输导层 的非均质特征 ； 然而直接在流体 势方程 中加入毛 管阻力值大小的方法似乎很难体现储层非均质性 对油气运移的影响。通常毛管力和流体势数量值 收稿 日期 2 0 0 8 0 5 0 4。 第一作者简 介 周波 1 9 7 6 一 ， 男 ， 博士 ， 油气成藏 动力学和油气 数值模拟 。 基金项 目 国家重点基础研究 发展规划 9 7 3 项 目 G1 9 9 9 043 3 1 0 。 维普资讯 学兔兔 w w w .x u e t u t u .c o m 528石油与天然气地质第29卷 相差太远，因而没有多大的实用价值，故未得到推 广。目前，利用 流体势方法 预 测 油气运移的基本 原理仍然是采用H uber t的流体势公式 ⋯。 。 前人通过物理实验 发现，当两互不相溶的流 体发生 驱替时，如果驱替速 度 非常慢，这时 流体 间的排替表现出极度不均 一 的性 状。传统 的描 述宏观均 匀混相 流 动 的水 动力 学 方程并不适 用 。 ” 一 引 ，这种流动特征需要新的理论模型来解 释。研究认为，逾渗模型可以描述这种流动 ”⋯ 。 在油气运移 的过 程中，油气运移的模 式是多 样化的 扎 坩j ，不同的动 力学条件或者基本 物理 条件对应着不同的运移模式；而 目前的数学模拟 过程都没有依据油气运移 的基本条件选择合适的 运移模式，而只是 随机选择某 种运移模型进行模 拟。本文针对目前油气运移研究 中运移模型选择 存在 的问题，分析油气运移数值模拟中基本数学 模型的选择问题，同时针 对数学模型到地质模型 的转换过程中存在的问题进行探讨。 l 目前主要的数值模拟方法 目前，用来描述 油气二次运移过程的方法 基 本上分为3类。 第 一 类，基 于流体势原理 睛。，“一 ，利用构造等 高线法线的汇聚程度，对油气运移过 程进行定性 分析，公式如下 p r一 中 gz Jdp／ pg 。／ 2 1 茜 式 中中为流体势，m 。／ s 。 ；z为高程，m；g为重力加 速度，s／m ；p为压力，Pa；p为流体密度，kg ／ m 。 ；q 为流体流速，，n ／s。 第二类，基于达西渗流方程的数值 模拟 方 法 ‘“ ”。 ，主要是利用 达西渗流方程和流体流动的 连续性方程，结合物质守恒定律，通过数学方法解 析渗流方程，获取 流体流动时油相 的饱和度、压力 等数据 旧’] 。 第三类，基 于流体非均 一 流动的逾渗数值模 拟方法 。一 。。 川 。侵入逾 渗理论 中运移过程的 基本原则是，在运移的每 一 步，侵入流体选择前缘 阻力最小的吼道突破。当流体进入 一 个新 的孑 L隙 后，流体不会沿着原来 的毛细管回流。阻力的大 小主要是依据毛管力、浮力、粘滞力以及外来驱替 力来综合判断。 2 油气运移的基本模式 根据前人 的研究，在油气运移的过 程中油 相 基本上是以游离相为主。因而，本文 的所有探讨 都是以油的游离相运移为前提的。 2．1 路径和优势路径 当油气发生运移时，即使在某个优势通道 内， 油气也并不是占据了整个运移通道。有些情况 下，油气沿着优势通道 的 一 部分发生运移图 1。 这里，首先对 油气运移的路径作 一 个定义。顾名 思义，运移路径实际上就是 油气实际走过的通 道 部分。 图l和图2为笔者进行 的 一 个二维玻璃珠实 验 的结果。图 中红色部分为注入的染 色煤油，浅 绿色部分为玻璃珠，玻璃珠粒径为0 ．8 ～ 0．9m m ， 模型大小为300 m m18 0m m 。实验 中，玻璃珠被 均匀地填装入二维模型中。很显然，该模型可以 看作是 一 个均匀的运移通道。如果给该通道注入 油，这个通道 就可以被 看作是优 势通道。当油气 在该介质中发生运移的时候，油只沿着介质中间 的 一 部分范围发 生运移，其他地方仍然 为饱和水 的玻璃珠介 质图 1。 也就是说，油气运移 的路 图1油气运移路径示意图 Fi g ．1T heschematicdiagI_ amofhydmcarbon mi g ration path wa y 维普资讯 学兔兔 w w w .x u e t u t u .c o m 第4期周波，等 ．油气二次运移数值模拟分析 529 图2油气再注入时运移路径变化示意图 Fi g ．2Theschematicdiagramshowingchangesofmi grati on path w aysdur ingrecha唱ing0fhydm carbons ＆注入速率0．3mI／min． b．注入速率O．6mL ／mi n 径和优势通道不是 一 个 概 念。即使 在 优 势通道 内，油气也 不 一 定占据整个通道，而是沿着优势通 道 的全部或者 一 部分范 围发生运移。 图2a为当注入油速度为0 ．3mL ／min时的油 气运移实验 结果图。当油气运移到顶部时，不改 变其他因素，直接将注入速度增加为O ．6mL ／min ， 当油再次运移到顶部时，获得 图2b。从 图2的运 移结果来看，首 先 当油 注入速度为0 ．3mL ／min 时，油沿着运移通道的 一 部分发生运移，运移路径 呈枝权状。当提高注入速度时，油运移的路径基 本不发生变化，油仍然沿着上次运移 的路径发生 运移。也就是说，当油注入速度发生变化时在 一 定范围内变化，油运移的路径并没有发生变化， 仍然沿着原来的路径发生运移。 图3为多个运移点时 的油气运移数学模拟 图。模拟 的基本数学模型采用逾渗模型 r 2引 。图 中，A和B为油气运移的两个源，c为运移的终点。 油气首先从A点出发，沿着红色部分运移，最终到 达c点。然后，油气又从B点出发，沿着蓝色 的路 径发生运移。当运移到达A点 时，运移不再产生 新的运移路径，而是 沿着原来的路径AC发生运 移，最后到达C点。 图3油气运移优势路径示意图 Fi g ．3T heschema ticdiaF啪showingthe pred or r Iinant onmi F a tion pat h way 红色部分为从A点出发的油的运移路径，蓝色部分为从口点 出发的油 的运移路径。其中为了描述 的方便，假定A点出发 的油 的运移路径形成后，B点的油才开始运移。 维普资讯 学兔兔 w w w .x u e t u t u .c o m 5 3 0 石 油 与 天 然 气 地 质 第 2 9卷 根据图 2和 图 3的实验结果 和数学模 拟结 果， 笔者尝试对优势路径 的概念进行定义。首先 ， 对于单源的油气运移 ， 当油气运移发生后 ， 早先产 生的运移路径基本保持不变 ， 即使后来 的物理条 件发生变化 从理论上讲 ， 物理条件的变化会引起 运移路径的变化 ， 如果该变化在一定范围 内， 运 移路径仍然基本保持不变 图 2 。而对于多源的 油气运移过程 图 3 ， 如果在运移 的优势通道 内 已经形成了一个运移路径 ， 新形成 的运移路径如 果遇到早先形成 的运移路径 ， 就会沿着早先形成 的运移路径运移 ， 并且会合后的路径保持不变。 也就是说 ， 仅就单个运移过程 ， 不存在优势路 径的概念。优势路径是多次运移或者多源运移 的 结果。通常 ， 将多次运移时保持不变 的路径和多 源运移 、 多条路径合并后的路径称作优势路径。 事实上， 通道是宏观概念， 代表了一种运输油 气的能力 ； 而路径是微观定量概念 ， 是指油气实际 的运移轨迹。 造成图 2和图 3两种运移结果的主要原因可 能是早先形成的运移路径上残存着一些不连续 的 油柱 ， 这些油柱 的存在使后来发生的油 气运移过 程更容易。此外 ， 早先 的运移路径上 的孔 隙被油 长期 占据 ， 介质的润湿性向亲油型转变， 这也会促 进油气的运移。 2 ． 2 油气运移模式 在前言中， 笔者 已经提到， 油气运移 的模式是 多样化的， 任何单一 的油气运移模式或者运移原理 都不可能描述所有的运移过程。因此 ， 必须对每个 油气运移的模式进行详细探讨。此外 ， 必须指 出的 是 ， 油气运移模式的划分是以运移路径为依据的， 这也是前面提出几个概念并进行区分的原因。 L e n o r m a n d等 通过大量 的物理实验并结合 数学模拟 ， 在不考虑浮力的情况下 ， 将两相流动过 程中侵入相流体 油相 的运移模式划分为 3种 毛管指进 式 、 粘滞 力 指 进 式 和 稳定 运 移 模 式。 Me a k i n等人【 3 ． 2 。 。 针对油气运移这种两相流动的过 程 ， 认为油气运移过程应该分为毛管指进模式 和 优势式运移模式。T u k u n a g a 根据物理实验 ， 利 用两个无量纲数 ， B o n d数和 C a数 ， 将 油气运移模 式划分为稳定运移模式和非稳定运移模式。罗晓 容等Ⅲ 根据物理实验获得 了 3种油气运移模式 稳定模式 、 指进模式和优势式运移模式。综合前 人的研究成果 ， 两非混相流体 中侵入相流体 的运 移模式为 稳定模式 包含浮力稳定模式和粘滞力 稳定模式 、 粘滞力指进模式、 毛管力指进模式和 优势模式。 根据实际地质条件下油气运移时两相流体的 物理状态 ， 分析前人成果 ] ， 认 为正 如罗晓 容 等物理实验所观察到 的结果一样 ， 油气运移 应该被分为 3种运移模式 稳定运移模式 粘滞力 稳定运移模式 、 指进模式 毛管力指进模式 和优 势式。 3 油气运移模式和数学模拟方法的 对应关系 根据前文分析 ， 油气运移的模 式被分为稳定 式 、 指进式和优势式 ； 而 目前的油气运移数值模拟 方法采用的基本理论也有 3种 达西渗流理论 、 流 体势理论和逾渗理论。达西渗流理论本身就是对 孔隙介质 内流体 的宏观均一性 流动进行描述 的， 主要的 目的是分析孔隙介质中宏观流动时不同部 位流量、 饱 和度等参数 的变化。而稳定式运移 的 特点是油气整体流动 ， 正好符合达西渗流理论 的 描述范围。因此， 基于达西渗流理论的数学模拟方 法能够用于描述油气运移过程中的稳定运移模式。 图 4是指进式运移模式和优势式运移模式分 别在物理实验 中的观察结果和数学模拟结果， 数学 模拟采用逾渗数学模型 J 。将物理实验结果和数 学模拟结果所描述 的油运移路径进行对 比可以发 现 ， 二者对应运移模式的运移形态基本一致。而前 人的研究也认为， 指进式运移和优势式运移应该属 于逾渗理论的描述范围 。两者结论一致。 依上论述 ， 目前采用的 3种数值模拟方 法和 油气运移模式的对应关系如图 5所示。达西渗流 理论适用于描述稳定模式下的油气运移过程 ， 而 逾渗理论适用于指进式和优势式运移模式下的油 气运移过程描述。而流体势理论本身的依据是流 体的势能， 其原理是物理学的基本原理， 因此无论 是什么模式的油气运移 ， 都可以用流体势理论来 描述。但是流体势的原理本身并不能描述油气运 移的具体过程， 只能对流体流动 的趋势作 出定性 的判断。因此 和前面两种方法相对 比， 前 两种方 法虽然描述对象有限， 但都属于定量描述方法 ； 而 流体势方法只是一种定性讨论方法。 维普资讯 学兔兔 w w w .x u e t u t u .c o m 第4期 周波，等 ．油气二次运移数值模拟分析 531 b 物理实验条件下的指进模式a和优势模式bb数学模拟条件下的指进模式c 和 优势模式d 图4指进模式和优势模式油气运移路径物理实验结果和数学模拟结果对比 F i g ．4Compa r isonoflabexper im entandm athemat jca lsimulationresul tso ff inger i ngandspr ingmi 铲at j onm odes 图5油气运移数值模拟方法和运移模式的对应关系 Fi g ．5Thecor r espondingr elat ionshi pb etweenmi gr at i on modesandsimulations 4 实际地质条件下 油气运移数值模 拟需要注意的几个问题 上文主要针对油气运移数值模拟时 的基础理 论 问题进行探讨。然而 ，实际地质条件下的油气 运移并不是简单的数学理论就可以计算或者描述 的。下文就目前地质条件下油气运移中的几个主 要问题进行讨论。 4．1 油气运移模式的确定 油气运移的模式是多样 化 的，不同的运移模 式必须采用不同的理论方法来开 展数学模拟。那 么在实际地质条件下，如何判 断某 个疏 导 层内的 油气运移属于哪种运移模式呢 根据前人的研究，可以依据两个无 量纲数 B ond数和Ca数绘制油气运移的模式图 ，然 后 根据模式图确定实际地质条件下的运移参数在模 式图中所处的位置，从而判断运移的模式。但是， 实际地质中有些参数是很难确定的，比如地质历 史上油气运移 的速度是很难精确估计的。因此实 际模拟时，应该结合油气勘探成果及地质、地化资 料，总结油气运移的规律，判断油气运移的特点， 进而选择模拟 方案，建立地质模拟，开 展数学 模拟。 在无法确定油气运移具体模式的情况 下，也 可以同时开展几种模式的数学模拟，然后和已知 的资料相对比，并结合流体势的信息，相互印证， 最终确定出合理的运移路线。 4．2尺度变换对不同数学模拟方法的影响 鉴于 目前的计 算机模拟水平，数学模 拟的计 算力是有限的。目前数学模拟所能提供的 网格数 目远 远小于实际地质条件下的单元数目，在这种 条件下，这种从宏观 地质体到数学模拟网格的转 变会对数学模拟的结果产生 一 定的影响。 根据研究，其他条件不变时，对于指进式运移 和优势式运移，当尺度发生变化 时，油气运移路径 占据整个运移空间的比例会 发生变化 嵋⋯ 。此外， 流体势模拟方法作为 一 种宏观定性模拟方法，并 不受尺度变化的影响。 4．3 数学模拟方法的局限性 数学模拟方法提供了定量模拟油气运移的工 具。但是，任何数学模拟都是建立在实际的地质 资料之上的。高品质和详细的地质资料是油气运 移模拟可靠性 的基本前提。 由于在勘探的不同阶段，地质 资料的掌握程 维普资讯 学兔兔 w w w .x u e t u t u .c o m 5 3 2 石 油 与 天 然 气 地 质 第2 9卷 度并不相同， 进行模拟时的孔渗和构造等资料本 身也是根据已知数据推测或者 网格化的结果， 因 此 ， 数学模拟的结果本身并不一定能够完全代表 实际的运移路径 ， 需要 和油气 显示结 果、 地 化资 料 相结合进行分析。 5 结论 1 油气运移数值模拟方法有 3种 流体势方 法 、 逾渗方法 和达西渗流方法。每个方法的原理 各不相同， 所描述的油气运移规律也不相同。 2 油气运移模式包 括 稳定式 、 指进式 和优 势式。其中， 逾渗方法适用 于描述指进式运移和 优势式运移过程 ， 达西渗流方程适用于描述稳定 式运移模式 ， 二者都属于定量描述方法 。而流体 势方法虽然同时适用于 3种模式下的运移 ， 但其 模拟结果只是一种定性结果。 3 数学模拟方法是 建立在基本 的地 质资料 基础上的， 因此必须以详实的地质资料为基础， 结 合油气显示和地化资料， 综合分析油气运移过程 ， 才能确定真正的油气运移过程 。 参考文献 1 Hi r s c h L M ， T h o mp s o n A H． Mi n i mu m s a t u r a t i o n s a n d b u o y an c y i n s e c o n d a r ym i g r a t i o n [ J ] ． AA P G B u l l e t i n ， 1 9 9 5 ， 7 9 5 6 9 6 7 1 0 2 Un g e ror P． Ba s i n e v a l u a t i o n b y i n t e gra t o r t wo - d i me n s i o n al mo d e - l i n g o f h e a t t r a n s f e r ， fl u i d fl o w， h y d r o c a r b o n g e n e r a t i o n an d mig r a - t i o n [ J ] ． A A P G B u l l e t i n ， 1 9 9 0， 7 4 3 1 2 2 1 1 2 3 0 3 Me a k i n P， W a g n e r G， V e d v i k A． I n v a s i o n p e r c o l a t i o n an d s e c o n d a r y m i g r a t i o n e x p e ri me n t s and s i m u l a t i o n s [ J ] ． M a r i n e and P e t r o - l e u m G e o l o g y ， 2 0 0 0， 1 7 7 7 7 7 7 9 5 4 Un g e r e r P， Be s s i s F， C h e n e t Y， e t a1． Ge o l o g i c al an d g e o c h e mi c al m ode l s i n o i l e x p l o r a ti o n p ri n c i p l e s and p r a c t i c al e x a mp l e s [ A] ． I n De ma i s o n G， e d ． P e t r o l e u m g e o c h e mi s t r y an d b asi n e v alu a t i o n [ C ] ． A A P G M e m o ir ， 1 9 8 4， 3 5 5 3 5 7 5 S c h o walt e r T T ． Me c h an i c s o f s e c o n d a r y h y d r o c a r b o n mi gra t i o n an d e n t r a p m e n t [ J ] ． A A P G B u l l e t i n ， 1 9 7 9 ， 6 3 5 7 2 3 7 6 0 6 Hu b b e r t M K． E n t r a p me n t o f p e t r o l e u m u n d e r h y d r o d y nam i c c o n - d i t i o n s [ J ] ． A A P G B u l l e t i n ， 1 9 5 3， 3 7 8 1 9 5 42 0 2 6 7 En g l an d D A． Th e mo v e me n t e n t r a p me n t o f p e t r o l e u m flu i d i n t h e s u b s u rf a c e [ J ] ． J o u r n al o f G e o l o g i c al S oci e t y ， 1 9 8 7， l 4 4 2 3 2 73 4 7 8 Rh e a L， P e r s o n M ， Ma r s i l y G D， e t a1． Ge o s t a t i s t i c al mod e l s o f sec - o n d a r y o i l mi gra t io n wi thi n h e t e r o g e n e o u s c a r ti e r b e d s a t h e o r e t i - c al e x a m p l e [ J ] ． A A P G B u l l e ti n ， 1 9 9 4， 7 8 1 1 1 6 7 91 6 9 1 9 B e t h ke C M， Re e d J D， Ol t z D F ． L o n g r a n g e p e t r o l e u m mi gra t i o n i n t h e I l l i n o i s B asi n [ J ] ． A AP G B u l l e t i n ， 1 9 9 1 ， 7 5 5 0 9 2 5 9 4 5 1 0 Da h l b e g E C．Ap p l i e d h y d r o d y n am i c s i n p e t r o l e u m e x p l o r a t i o n [ M] ． N e w Y o r k S p ri n g e r - V e r la g ， 1 9 8 2 ． 1 1 7 1 1 1 Hi n d l e A D． P e t r o l e u m mi g r a t i o n p a t h wa y s a n d c h arg e c o n c e n t ra- t i o n a thr o e - d i m e n s i o n al mode l 【 J] ． A A P G B u ll e ti n， 1 9 9 7 ， 8 1 9 1 4 1 1 4 8 1 2 L e n o r ma n d R， Z arc o n e R， T o u b o u l E． Nu mb e r i c al mod e l s an d e x p e r i me n t s o n i m m i s c i b l e d i s p l a c e m e n t s i n p o r o u s me d i a [ J ] ． J o u r - n al o f F l u i d M e c h a n i c s ， 1 9 8 7， 1 8 9 2 1 6 51 8 7 1 3 W i l k i n s o n D． P e rco l a t i o n mod e l o f i mmi sc i b le d i s p l a c e me n t i n the p r e sen c e o f b u o y anc y f o rce s [ J ] ． P h y s ic al R e v i e w A， 1 9 8 4， 3 0 1 5 2 0 5 3 1 1 4 Wi l k i n son D． P e r c o l a t i o n e ff e c ts i n i m mi sci b l e d i s p l a c e me n t [ J ] ． P h y s i c al R e v i e w A， 1 9 8 6 ， 3 4 2 1 3 8 01 3 9 1 1 5 T u k u n a g a T， Mo g i K ， Ma t s u b ara O， e t a1． Bu o y a n c y a n d i n t e r f a c i al f o rce e ff e c t s o n t w o - p h ase d i s p l a c e me n t p a t t e r n s a n e x p e rime n t al s t u d y [ J ] ． A A P G B u ll e t i n ， 2 0 0 0， 8 4 1 6 5一 l 6 张发强， 罗晓容 ， 苗盛， 等． 石油二次运移 的模式及其 影响因素 [ J ] ． 石油实验地质 ， 2 0 0 3 ， 2 5 1 6 7～ 7 4 1 7 L u o X ， Zh a n g F， Mi a o S， e t a1． E x p e rime n t al v e rifi c a t i o n o f o i l s a t u r a t i o n and l o s s e s d u ri n g s e c o n d a r y m i gr a t i o n [ J ] ． J o u rnal o f P e - t r o l e u m Geo l o g y ， 2 0 0 4 ， 2 7 3 2 93 9 l 8 Ca t a l an L， Xi a o we n F， Ch a t z i s I ， e t a 1 ．An e x p e rime n t al s t u d y o f s e c o n d a r y o i l mi gra t i o n [ J ] ． A A P G B u l l e t i n ， 1 9 9 2， 7 6 5 6 3 8 6 4 0 l 9 侯平． 石油二次运移路径 的模式分析 [ J ] ． 中国科学 D辑 ， 2 0 0 4， 3 4 增刊 I 1 6 2～1 6 8 2 0 Me a k i n P ． I n v asi o n p e rco l a t i o n i n a d e s t a b i l i z i n g gra d i e n t [ J ] ． P h y s i c al R e v i e w A， 1 9 9 1 ， 4 6 6 3 3 5 7 3 3 6 8 21 Me h e u s t Y， L C v o l l G， Mal y K J ， e t a1． I n t e rfa c e scali n g i n a t wo - d i me n s i o n al p o r o u s me d i a u n d e r c o mb i n e d v i sco u s ，gra v i t y an d c a p i l l a r y e ff e c t s [ J ] ． P h y s i c al R e v i e w E， 2 0 0 2 ， 6 6 5 0 5 1 6 0 3 2 2 Wa g n e r G， B i rov l j e v A， M e a k i n P ， e t a1． F r a g m e n t a t i o n and m i g r a t i o n o f i n v as i o n p e rco l a t i o n c l u s t e r e x p e rime n t s an d s i mu l a t i o n s [ J ] ． P h y s i c al R e v i e w E， 1 9 9 7 ， 5 5 6 7 0 1 5 7 0 2 9 2 3 金之钧 ， 张发强． 油气运移研究现状及主要进展[ J ] ． 石油与天 然气地质， 2 0 0 5 ， 2 6 3 2 6 3～ 2 7 0 2 4 罗晓容， 张立宽， 廖前进， 等． 埕北断阶带沙河街组油气运移动 力学过程模 拟分析 [ J ] ． 石 油 与天 然气 地 质 ， 2 0 0 7 ， 2 8 2 l 9 l～l 9 7 2 5 周波 ， 金之钧 ， 罗晓容， 等． 尺度放大时逾渗模型中的油气运移 路径变化规律 探讨 [ J ] ． 石 油 与天 然气 地质 ， 2 0 0 7 ， 2 8 2 1 7 5～l 8 0 2 6 Z h o u Bo ， Di d i e r Lo g g i a， L u o Xi a o r o n g， e t a1．N u me ri c al s tu d i e s o f gra v i t y